Previous Article | Next Article ![]()
Applied and Environmental Microbiology, December 2004, p. 7355-7364, Vol. 70, No. 12
0099-2240/04/$08.00+0 DOI: 10.1128/AEM.70.12.7355-7364.2004
Copyright © 2004, American Society for Microbiology. All Rights Reserved.
Oliver Bridle, and
Dave J. Scanlan*
Department of Biological Sciences, University of Warwick, Coventry, United Kingdom
Received 17 May 2004/ Accepted 3 August 2004
|
|
|---|
|
|
|---|
240 Tg of N year1 and a marine contribution of
80 Tg of N year1 (4). Trichodesmium, a filamentous nonheterocystous cyanobacterium broadly distributed throughout tropical and subtropical oceans, has long been considered to be responsible for most of this marine N2 fixation (see references 3 and 17 for recent reviews). However, recently, Zehr et al. (34) showed that unicellular cyanobacteria that are 3 to 10 µm in diameter may make a significant contribution to oceanic N2 fixation. Indeed, based on N2 fixation data for the 0.2- to 10-µm bacterioplankton size fraction and the concentrations of phycoerythrin-containing unicellular cyanobacteria, they estimated that this contribution might equal or exceed that of Trichodesmium. Several unicellular marine cyanobacterial diazotrophs belonging to the genera Synechococcus (15, 16) Cyanothece (21), and Crocosphaera (28) have previously been described. These organisms all have temporal separation of N2 fixation and photosynthesis, which is likely controlled by an endogenous circadian rhythm (22). In addition to these free-living isolates, 16S rRNA gene (rDNA) sequences derived from cyanobacterial symbionts of the marine diatom Climacodium fraeuenfeldianum show close phylogenetic relatedness, suggesting that unicellular cyanobacterial symbionts also have the potential to contribute to N2 fixation (7, 9, 10). Although several studies have now reported the presence of unicellular phycoerythrin-containing cyanobacteria that are 3 to 10 µm in diameter in geographically separated marine environments (2, 18, 27), suggesting that there is widespread distribution, the exact contribution of these unicellular cyanobacterial strains to global N2 fixation rates remains to be determined. Indeed, a means to specifically detect or analyze the diversity of these organisms is not yet available. To address this problem, we developed a 16S rDNA oligonucleotide primer that specifically recognizes the discrete marine unicellular cyanobacterial diazotroph lineage within which nearly all previously documented isolates lie. We also describe the utility of this oligonucleotide as a PCR primer for assessing the distribution of these organisms along a transect in the Arabian Sea situated in the northwest Indian Ocean.
|
|
|---|
|
View this table: [in a new window] |
TABLE 1. PCR primers used in this study
|
Phylogenetic analysis and development of the NITRO821R primer.
16S rDNA sequence alignment and phylogenetic analysis were performed by using the ARB software (13). The NITRO821R oligonucleotide primer (Table 1) was designed to specifically recognize the phylogenetic lineage in which the unicellular cyanobacterial diazotrophs Cyanothece sp. strains WH 8902, WH 8904, and PCC8801 (also known as Synechococcus sp. strain RF-1 [25]) and Crocosphaera sp. strain WH 8501 lie (Fig. 1) by using the probe design and probe match tools from the ARB program.
![]() View larger version (42K): [in a new window] |
FIG. 1. Phylogenetic tree showing the relationships of marine and freshwater cyanobacteria, including diazotrophs and nondiazotrophs, inferred from 16S rDNA sequences. The tree was constructed by the neighbor-joining method with Jukes-Cantor correction. The percentages of bootstrap replicates supporting the branching order are indicated at the nodes. Partial sequences (<1,190 nucleotides) were added to the tree by using a maximum-parsimony option within ARB. The cyanobacterium Gloeobacter violaceus PCC7421 was used as a root. The scale bar represents the equivalent of 0.1 substitution per nucleotide.
|
Sampling.
Water samples for DNA extraction were collected during the NERC-funded AMBITION cruise in the Arabian Sea from 1 to 27 September 2001 aboard the RRS Charles Darwin. Samples were obtained from discrete depths at 11 stations along a 5,500-km transect (Table 2 and Fig. 2) between Victoria, Seychelles, and Muscat, Oman, by using 20-liter Niskin bottles on a hydrographic cable. Conductivity, temperature, and barometric pressure were measured simultaneously with a CTD (model Sea-Bird 9/11). Seawater (5 to 10 liters) from each depth was filtered onto 47-mm-diameter, 0.45-µm-pore-size polysulfone filters (Supor-450; Gelman Sciences Inc., Ann Arbor, Mich.) after prefiltration through a 47-mm-diameter, 3-µm-pore-size filter (MCE MF-Millipore filters [Fisher]) with a gentle vacuum (10 mm of Hg). This allowed collection of >3- and 3- to 0.45-µm fractions. The filters were placed in a 5-ml cryovial with 3 ml of DNA lysis buffer (0.75 M sucrose, 400 mM NaCl, 20 mM EDTA, 50 mM Tris HCl [pH 9.0]) and stored at 80°C until extraction.
|
View this table: [in a new window] |
TABLE 2. Positions of the principal stations along the AMBITION cruise transect
|
![]() View larger version (30K): [in a new window] |
FIG. 2. Distribution of the UCYN2-fix lineage along the AMBITION cruise transect in the Arabian Sea. The sampling stations are indicated by dots and are numbered in the order of sampling. Representative agarose gels for stations giving detectable PCR products with the CYA359F-NITRO821R primer pair are shown on the right. The numbers above the lanes indicate the depth of sampling (in meters). Sta., station.
|
Development of the dual analytical PCR.
For specific amplification of the 16S rDNA of marine unicellular cyanobacterial diazotrophs, the specific NITRO821R primer was used in conjunction with the cyanobacterium- and plastid-specific forward primer CYA359F (19) (Table 1). For amplification of environmental nifH sequences we used the cyanobacterium-biased primer pair described by Olson et al. (20), which we designated CNF-CNR (Table 1). PCRs were carried out in 25-µl mixtures containing 0.5 µl of template (cells or environmental DNA sample), each deoxynucleoside triphosphate at a concentration of 200 µM, 2 mM MgCl2, each primer at a concentration of 0.2 µM (for 16S rDNA amplification) or 1 µM (for nifH amplification), and 0.625 U of Taq polymerase in 1x enzyme buffer (GIBCO BRL, Life Technologies Ltd., Paisley, Scotland). For amplification of environmental DNA we also included 1 mg of bovine serum albumin (Sigma) ml1. The amplification conditions comprised a denaturation step of 95°C for 5 min and then 80°C for 1 min, at which time Taq polymerase was added, followed by 30 cycles of 95°C for 1 min, 57°C for 30 s, and 72°C for 40 s and a final extension of 6 min at 72°C. Reaction mixtures were stored at 4°C prior to analysis. Products (10 of 25 µl) were resolved by gel electrophoresis on a 1.5% (wt/vol) agarose gel at 100 V. DNA was stained with ethidium bromide (0.5 µg ml1), visualized under short-wavelength UV, and photographed with a gel documentation system (UVP Inc., Upland, Calif.).
Relative quantification of PCR products.
Ethidium bromide-stained agarose gels were quantified by using the array analysis tool in the TotalLab software (Nonlinear Dynamics).
Clone libraries.
Eight clone libraries were constructed, four each for nifH and 16S rDNA, by using the CNF-CNR and CYA359F-NITRO821R primers, respectively. The libraries were constructed by using environmental DNA from the AMBITION cruise transect obtained at station 2 (depth, 10 and 25 m), station 4 (depth, 25 m), and station 11 (depth, 10 m). PCR products from each depth were ligated into the pCR2.1-TOPO cloning vector before transformation into Escherichia coli strain TOP10F (Invitrogen Corporation, San Diego, Calif.). Plasmid DNA was isolated from transformants with a QIAprep miniprep kit (QIAGEN Ltd., Crawley, West Sussex, United Kingdom). Ten clones were subsequently sequenced from each library prior to phylogenetic analysis by using the ARB software as described above.
Nucleotide sequence accession numbers.
16S rDNA and nifH sequences reported in this paper have been deposited in the GenBank database under accession numbers AY620237 to AY620241 (cultured strains) and AY621666 to AY621747 (AMBITION cruise environmental sequences).
|
|
|---|
|
View this table: [in a new window] |
TABLE 3. Specificity of the NITRO821R oligonucleotide illustrated, by alignment of the primer and target sequences with the 16S rDNA sequences of different marine unicellular diazotrophs and other cyanobacteria
|
![]() View larger version (20K): [in a new window] |
FIG. 3. Specificity of amplification with 16S rDNA primers CYA359F and NITRO821R (a) or nifH primers CNF and CNR (20) (b). Lane Av, Anabaena variabilis; lane T, Trichodesmium sp. strain WH 9601; lane A, Anabaena sp. strain PCC7118; lane Pr, Planktothrix rubescens; lane S, Synechocystis sp. strain PCC6803; lane U, Synechococcus sp. strain CC9605; lane Y, Cyanothece sp. strain WH 8902; lane C, Cyanothece sp. strain WH 8904; lane PHLI, Prochlorococcus sp. strain EQPAC1; lane PHLII, Prochlorococcus sp. strain TAK9803-2; lane PLL, Prochlorococcus sp. strain MIT9313; lane N, Nostoc sp. strain PCC7120; lane M, Microcystis sp. strain PCC7806; lane G, Gloeotrichia sp.; lane Co, control (no DNA).
|
![]() View larger version (21K): [in a new window] |
FIG. 4. Limit of detection of the CYA359F-NITRO821R PCR. The numbers of Cyanothece sp. strain WH 8902 cells in the 25-µl reaction mixtures are indicated at the top.
|
![]() View larger version (17K): [in a new window] |
FIG. 5. Relative abundance of the UCYN2-fix lineage along the AMBITION cruise transect in the Arabian Sea. The error bars indicate standard errors. The inset shows the data for stations 3, 4, and 11 in more detail.
|
Genetic diversity of cyanobacterial diazotrophs along the AMBITION cruise transect.
Clone libraries of 16S rDNA and nifH were constructed by using the CYA359F-NITRO821R and CNF-CNR primer pairs, respectively, for station 2 at depths of 10 and 25 m, for station 4 at a depth of 25 m depth, and for station 11 at a depth of 10 m. All sequenced 16S rDNA clones from these environmental libraries fell in the UCYN2-fix lineage (Fig. 6). The majority of the sequences (27 of 40 sequences), all derived from station 2 or 11, were phylogenetically most closely related to Crocosphaera sp. strain WH 8501. Within this group of sequences, those obtained from station 11 formed a distinct clade. 16S rDNA sequences from station 4 formed two separate clusters within the lineage and were generally not closely related to any known cultured unicellular diazotroph. Of the 40 environmental sequences obtained along the transect, clone 8 obtained from station 4 at a depth of 25 m and clone 5 obtained from station 11 at a depth of 10 m showed the greatest diversity (93.6% identity for the 462-bp product), suggesting that there is wide genetic diversity of these organisms in the Arabian Sea region.
![]() View larger version (25K): [in a new window] |
FIG. 6. Phylogenetic tree showing the relationships of the AMBITION cruise CYA359F-NITRO821R environmental 16S rDNA sequences. Each AMBITION cruise environmental sequence is referred to by station number, depth of sampling, and clone (cl.) number. The tree was constructed by the neighbor-joining method with Jukes-Cantor correction based on nearly full-length 16S rDNA sequences from marine and freshwater cyanobacteria. The percentages of bootstrap replicates supporting the branching order are indicated at the nodes. Partial sequences (e.g., sequences of the 462-bp product amplified by the CYA359F-NITRO821R primer set and other sequences <1,190 nucleotides long) were added to the tree by using a maximum-parsimony option within ARB. The cyanobacterium Gloeobacter violaceus PCC7421 was used as a root. The scale bar represents the equivalent of 0.1 substitution per nucleotide.
|
![]() View larger version (26K): [in a new window] |
FIG. 7. Phylogenetic tree showing the relationships of the AMBITION cruise environmental nifH sequences. Each AMBITION cruise environmental sequence is referred to by station number, depth of sampling, and clone (cl.) number. The tree was constructed by the neighbor-joining method with Jukes-Cantor correction based on partial nifH sequences (e.g., a 359-bp fragment amplified with the CNF-CNR PCR primer set used in this study). The percentages of bootstrap replicates supporting the branching order are indicated at the nodes. The cyanobacterium Symploca sp. strain PCC8002 was used as the root. The scale bar represents the equivalent of 0.1 substitution per nucleotide.
|
![]() View larger version (15K): [in a new window] |
FIG. 8. Plot of the relative abundance of the UCYN2-fix lineage along the AMBITION cruise transect as a function of seawater temperature.
|
|
|
|---|
Several 16S rDNA sequences of freshwater unicellular cyanobacterial diazotrophs which represent the genera Gloeocapsa, Gloeothece, and Aphanothece (accession numbers AB067575 to AB067581 and AB119259) and which are in the UCYN2-fix lineage described here (data not shown) have recently been deposited in the GenBank database. These sequences are also targeted with complete identity by the NITRO821R oligonucleotide, which potentially broadens its utility to freshwater systems. However, care is needed when this primer is used to enumerate this lineage in freshwater environments, given that we have shown that sequences of members of the genus Planktothrix, with a 1-bp mismatch with the NITRO821R primer, are amplified in PCRs with the CYA359F-NITRO821R primer pair. Specific size fractionation of freshwater environmental samples to remove the filamentous Planktothrix spp. would be one route to circumvent such a problem.
The utility of this molecular approach for revealing the distribution of unicellular cyanobacterial diazotrophs in natural marine systems appears to be well demonstrated, however, by the data obtained along the AMBITION cruise transect. Although the relative abundance compared to the nondiazotrophic picocyanobacterial genera Prochlorococcus and Synechococcus, which were the dominant organisms at the southern and northern stations, respectively (G. Tarran unpublished data), was low, significant numbers of the UCYN2-fix lineage were detectable, as determined by the semiquantitative analytical PCR developed in this study. Interestingly, the range of cell abundance observed along this single transect (7 x 103 ± 1.6 x 102 to 3.8 x 105 ± 8 x 104 cells liter1) is very similar to the range of average cell counts of unicellular cyanobacteria in the 2.5- to 7-µm size class reported by Falcon et al. (10) for the tropical North Atlantic and subtropical North Pacific. In that study the cell numbers in the two ocean systems were considerably different, and it was proposed that the higher Fe flux to the North Atlantic was one of the factors that were important in determining the abundance of these organisms. Although we have no data on Fe flux at the different stations along the AMBITION cruise transect, this ocean basin is known to have high and variable deposition of dust (24), which results in pulsed iron addition in this region (14). Certainly, the fact that markedly different relative abundance values for the UCYN2-fix lineage were observed along this transect shows that the abundance of these organisms can vary considerably over relatively small spatial scales. Indeed, a comparison of a suite of environmental parameters obtained from the AMBITION cruise suggests that seawater temperature (Fig. 8) accompanied by low surface SRP and nitrogen levels (generally <0.25 µM) is significant in shaping the conditions in which these organisms can proliferate. Temperature has recently been shown to be important in dictating the distribution of heterocystous and nonheterocystous cyanobacteria in oceanic systems (23). In that study, differences in the temperature dependence of oxygen flux, respiration, and N2 fixation activity were shown to be critical in explaining how Trichodesmium performs better than heterocystous species at higher temperatures, and it may be that some of these factors are important for unicellular nitrogen-fixing species too.
The construction of environmental clone libraries based on the CYA359F-NITRO821R primer pair allowed us to obtain direct insight into the diversity of the UCYN2-fix lineage along the AMBITION cruise transect. The wide genetic diversity among the derived 16S rDNA sequences for the three stations analyzed was noticeable (the lowest level of identity was 93.6% for the 462-bp product), and this diversity was much greater than the diversity that has been found for the picocyanobacterial genera Prochlorococcus and Synechococcus (11, 26), for which the lowest levels of identity between strains and environmental sequences were 96 to 97% (for sequences usually >1,000 bp long) in their respective lineages. How this translates to physiological diversity requires further culture isolation since several of the sequences, particularly those from station 4, had no closely related cultured counterparts. Within a station, particularly stations 2 and 11, the genetic diversity observed was generally lower, and interestingly, sequences from the same station formed discrete clusters or clades (Fig. 6), suggesting the potential for ecotype adaptation of members of the UCYN2-fix lineage in specific water columns. The utility of the CYA359F-NITRO821R primer pair for focusing directly on the UCYN2-fix lineage is mirrored by data for the nifH clone libraries, where the primers used are not targeted to a particular cyanobacterial lineage. Here, sequences closely related to the Trichodesmium-Katagnymene group dominated the libraries. Two sequences from station 2, however, were closely related to the nifH sequence of Crocosphaera sp. strain WH 8501, which was in good agreement with the 16S rDNA environmental sequence data from that station, where sequences closely related to this strain dominated. These sequences thus correspond to the group B nifH sequences designated by Zehr et al. (34), while the group A nifH sequences reported in that study may be related to the deeply branching 16S rDNA clusters that include the station 4 clones.
Given the wide utility of the nifH gene for investigating the diversity of marine nitrogen fixers (32, 33, 35), it is noteworthy that we also obtained two environmental sequences (from stations 2 and 11) whose closest known relatives are nifH sequences from Anabaena sp. We found no previous reports of nifH sequences from members of this genus for oligotrophic open-ocean marine environments. However, these sequences may well be equivalent to the sequence of the novel Anabaena sp. described by Carpenter and Janson (8), who reported the presence of a heterocystous cyanobacterium in the southwest Pacific Ocean and Arabian Sea at a low level (1 to 4 trichomes liter1).
The NITRO821R oligonucleotide and the analytical PCR developed here thus provide an excellent tool for rapid and sensitive screening of environmental samples to establish the presence and relative abundance of unicellular cyanobacterial diazotrophs, organisms that play important roles in oceanic new production and biogeochemistry. Use of the oligonucleotide in fluorescent in situ hybridization or quantitative PCR experiments, in which absolute quantitation of cell numbers can be obtained, or use in reverse transcription-PCR to determine the members of the community that are active in situ would further enhance its value.
This work was funded by NERC grant GST/02/2819 under the Marine and Freshwater Microbial Biodiversity initiative, as well as by research development grant SC 02009 from the University of Warwick. O.B. was the recipient of a Society for General Microbiology summer studentship.
Present address: Darling Marine Center, University of Maine, Walpole, ME 04573. ![]()
|
|
|---|
This article has been cited by other articles:
| |||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||
Copyright © 2009 by the American Society for Microbiology. For an alternate route to Journals.ASM.org, visit: http://intl-journals.asm.org | More Info»