Previous Article | Next Article ![]()
Applied and Environmental Microbiology, September 2003, p. 5157-5169, Vol. 69, No. 9
0099-2240/03/$08.00+0 DOI: 10.1128/AEM.69.9.5157-5169.2003
Copyright © 2003, American Society for Microbiology. All Rights Reserved.
Laboratoire d'Algologie, de Mycologie et de Systématique Expérimentale, Institut de Botanique B22,1 Centre d'Ingénierie des Protéines, Institut de Chimie B6, Université de Liège, B-4000 Liège, Belgium,2 Deutsche Sammlung von Mikroorganismen und Zellkulturen GmbH, D-38124 Braunschweig, Germany,3 Laboratoire d'Océanographie Biologique, UMR 5805 CNRS and Université Bordeaux 1, F-33120 Arcachon, France4
Received 14 March 2003/ Accepted 9 June 2003
|
|
|---|
|
|
|---|
Molecular ecology studies have greatly enhanced our ability to detect and identify microorganisms in nature (1). Moreover, the need for molecular tools to study the genotypic relationships of cyanobacteria, to reconstruct their evolution, and to improve their taxonomy has long been recognized (69). Garcia-Pichel et al. (17) pointed out that the cyanobacterial genotypic diversity and the diversity inferred from identification of morphotypes may or may not coincide depending on the specific case analyzed and suggested that the correlation between morphotypes and phylotypes became gradually better as more complex forms were considered.
One of the questions that now can be addressed with molecular tools is the geographic distribution of the cyanobacteria and the existence of endemic taxa in Antarctica. Currently, there is no consensus concerning the extent of endemism of cyanobacteria in polar environments (62, 63). Arguments against the existence of endemic cyanobacterial species in Antarctica are related to the efficient dispersal abilities of the organisms and the relatively young age of the ice-free area. Moreover, most of the Antarctic cyanobacteria identified to date based on morphological features have apparently cosmopolitan distributions (62). However, Komárek (27) attributed the assignment of Antarctic cyanobacterial taxa to cosmopolitan species to the use of taxonomic keys developed for temperate or tropical microflora. In support of the hypothesis that there are endemic cyanobacteria in polar environments (63) are the fact that Antarctica has been more isolated than other parts of the world for several million years, the fact that dispersal processes which favor local species are more efficient than long-range dispersal processes, and the observation that there has probably been environmental selection for adaptative strategies. Accordingly, on the basis of morphology, Komárek (27) determined that about 60% of the species in various microbiotopes of ice-free areas of King George Island are probably endemic. However, it is necessary to obtain sequence data because there is evidence that morphological features do not reflect the real genetic and physiological divergence (63). Recent sequence information (37, 44, 52) has shown that there are clusters containing only Antarctic sequences. This could support the hypothesis that there is cyanobacterial endemism. Hence, molecular techniques can provide new insights into the genotypic diversity of the Antarctic cyanobacterial communities and their affinities with Arctic and nonpolar communities. These techniques provide a different way of sampling microbial diversity that is complementary to cultivation and microscopy (65).
This study was part of a comprehensive biodiversity study of a microbial mat from Lake Fryxell performed in the framework of the joint research project MICROMAT (Biodiversity of Microbial Mats in Antarctica [http://www.nerc-bas.ac.uk/public/mlsd/micromat/]). The same mat sample was shared by the different research groups and was used for determinations of bacterial and archaeal diversity (8, 59, 61), as well as an inoculum for culturing mat communities with a benthic gradient chamber (BGC) (11). The artificial communities cultured in the BGC at 5°C were clearly cold adapted, as oxygenic photosynthesis was optimal at 10°C, and they shared several features with the field mat, such as pinnacle formation and pigment composition (11, 43).
To describe the morphological and genotypic diversity of cyanobacteria from a Lake Fryxell mat sample, we used a polyphasic approach. The morphological diversity was assessed by light microscopy, and the genotypic diversity was studied by examining clone libraries and by performing denaturing gradient gel electrophoresis (DGGE) based on rRNA gene amplicons. The internal transcribed spacer (ITS) has proven to be useful for taxonomic studies of cyanobacterial cultures (3, 7, 13), as well as for bacterial community studies (5). Therefore, we retrieved ITS sequences from the cyanobacterial communities to complement the information obtained from the 16S rRNA gene. Using the same approach, we analyzed the cyanobacteria from the BGC mat, which originated from Lake Fryxell and showed functional similarities with the natural community. Finally, by performing a detailed phylogenetic analysis we compared the cyanobacteria observed in this study with previously described polar cyanobacteria to study possible endemism in Antarctica and to make comparisons between polar regions.
|
|
|---|
Cultivation of artificial mats.
An artificial microbial mat was cultured in a BGC with opposing oxygen and sulfide gradients at 5°C for 8.5 months at the University of Bordeaux. Further details of the cultivation procedure and a structural description of the mat have been provided by Buffan-Dubau et al. (11). The artificial mat grew mainly above the sand surface and reached a final thickness of 2.5 to 5 mm. The surface of the mat was dark green with a few light green and brown patches. Pinnacle formation similar to that described for Antarctic mats was observed on the surface. At the end of the culture period the BGC was opened, and the artificial microbial mat was collected. A subsample was stored in 50% ethanol and was sent to the University of Liege for subsequent analyses.
Microscopic observations.
Immediately upon receipt, subsamples were fixed with 4% formaldehyde. The major cyanobacterial taxa were observed with a Wild MS-20 microscope equipped with a screw micrometer (68). The diacritical morphological traits used for botanical species descriptions were considered (cell shape for both intercalary and end cells; width and length of intercalary cells; presence or absence of constriction at the cross wall, of necridic cells, and of a sheath; color of the sheath; number of trichomes per filament; presence or absence of heterocysts; width and length of heterocysts). For each biometrical character, 30 measurements were obtained from cells, heterocysts, and filaments sampled at random. The taxonomic works of Geitler (18), Anagnostidis and Komárek (2), and Komárek and Anagnostidis (26), as well as descriptions in the Antarctic literature (9, 10, 33), were used.
Isolation of DNA.
The DNA extraction protocol used in this study was derived from that of Smalla et al. (51). A 0.5-g mat sample was suspended in 0.5 ml of SNT solution (500 mM Tris-HCl, 100 mM NaCl, 25% saccharose) supplemented with 26 µl of freshly added lysozyme (25%, wt/vol). The resulting suspension was shaken and incubated for 30 min at 37°C. After incubation, 0.5 ml of solution II (500 mM Tris base, 500 mM EDTA, 1% sodium dodecyl sulfate, 6% phenol) and 0.25 g of glass beads (diameter, 0.17 to 0.18 mm; Braun Biotech) were added to the sample and shaken for 60 s in a bead beater (Braun Biotech). The resulting suspension was placed on ice for 1 h and vortexed every 10 min. After incubation, the suspension was centrifuged for 10 min at 720 x g (MinifugeT; Heraeus), and 1 ml of the aqueous phase was mixed with an equal volume of phenol, after which it was centrifuged for 5 min at 13,600 x g. The supernatant was then transferred into new tubes, extracted with equal volumes of phenol-chloroform-isoamyl alcohol (25:24:1), and reextracted with equal volumes of chloroform-isoamyl alcohol (24:1). Then a standard Na acetate-ethanol precipitation was performed, and the dried pellet was resuspended in 100 µl of TE buffer (10 mM Tris-Cl, 1 mM EDTA; pH 8).
Cloning of the ribosomal DNA and screening by amplified ribosomal DNA restriction analysis (ARDRA).
The crude DNA preparations were subjected to several cesium chloride (CsCl)- and potassium acetate-based purifications steps by using the protocol of Smalla et al. (51), except for a few modifications. Two hundred milligrams of CsCl was added and mixed, and the tubes were centrifuged for 10 min at 13,600 x g after 2 h of incubation at 20°C (with mixing every 30 min). Three volumes of water and 0.6 volume of 80% isopropanol were added to the supernatant, and the tubes were centrifuged for15 min at 13,600 x g after 30 min of incubation at 20°C. After this the supernatant was discarded, and the pellet was dried. The DNA was dissolved by adding 100 µl of TE buffer, followed by addition of 20 µl of 8 M potassium acetate and 1.5 h of incubation at 20°C. After 15 min of centrifugation at 13,600 x g, 0.6 volume of 80% isopropanol was added to the supernatant, mixed, and incubated for 30 min at 20°C. After centrifugation for 15 min at 13,600 x g and 4°C, the supernatant was discarded, and the pellet was dried. The DNA was dissolved in 200 µl of TE buffer.
PCR amplification of a cyanobacterial 16S rRNA gene plus the ITS was performed in a 50-µl (total volume) reaction mixture containing 0.5 µl of mat DNA, 1x Super Taq Plus PCR buffer, each deoxynucleoside triphosphate at a concentration of 0.2 mM, 0.5 µM primer 16S27F (Table 1), 0.5 µM primer 23S30R (Table 1), 1 mg of bovine serum albumin (Sigma Chemical Co., St. Louis, Mo.) ml-1, and 1 U of Super Taq Plus polymerase with proofreading activity (HT Biotechnology, Cambridge, United Kingdom). Amplification was carried out with a Gene Cycler (Bio-Rad, Hercules, Calif.) as follows: one cycle of 5 min at 94°C; 10 cycles of 45 s at 94°C, 45 s at 57°C, and 2 min at 68°C; 25 cycles of 45 s at 92°C, 45 s at 54°C, and 2 min at 68°C; and a final elongation step of 7 min at 68°C. PCR products were purified with a Quantum Prep PCR Kleen Spin column (Bio-Rad). Poly(A) extension for 20 min at 72°C with dATP and the Goldstar polymerase (Eurogentec, Seraing, Belgium) was performed.
|
View this table: [in a new window] |
TABLE 1. Primer sequences and target sites
|
Plasmid DNAs were extracted with a Quantum Prep Plasmid Miniprep kit (Bio-Rad) by following the manufacturer's instructions. The inserted 16S rRNA gene plus ITS was reamplified with primers 16S27F and 23S30R as described above and subjected to ARDRA to screen the clone libraries. For the ARDRA, all steps were performed as previously described (49), except for the following changes: MboI and HpaII (Gibco Life Sciences) were used as restriction enzymes, and electrophoresis was performed at a constant voltage of 3 V cm-1 for 235 min.
For each ARDRA type, partial 16S rRNA gene sequences that were at least 400 bp long were determined for one to three clones by using sequencing primer 16S378F and/or primer 16S784R (Table 1). Sequencing was carried out by GenomeExpress (Paris, France) with an ABI PRISM system 377 (PE Applied Biosystems, Foster City, Calif.). Complete 16S rRNA gene sequences (Escherichia coli positions 27 to 1492) were determined for one DNA strand for one clone selected at random from among the clones identified as members of the same phylotype. We define phylotype as a group of sequences that exhibit more than 97.5% similarity with each other.
We set the threshold at 97.5% for two reasons. First, for complete 16S rRNA gene sequences, a binary similarity value less than 97.5% likely corresponds to a DNA-DNA hybridization value less than 70%, and therefore the sequences probably correspond to two different species (56). Second, there is a good correlation between complete and partial 16S rRNA gene sequence similarity values. We calculated binary similarity values for all complete 16S rRNA gene sequences of cyanobacterial strains available from GenBank (indels were not taken into account). The 53 strains with binary similarity values between 96 and 98% were selected, and their partial 16S rRNA gene sequence similarities were calculated (for E. coli positions 378 to 784). The average similarities were 97.05 and 97.01% for complete and partial sequences, respectively. Hence, on average, partial sequence similarity reflected the complete sequence similarity well, although the standard deviation was 1.24.
The sequencing was carried out with primers 16S1494R and 16S784R (Table 1). The consensus sequences were obtained with the software AlignIR (LI-COR). In addition, the complete ITS sequences were determined by using primer 23S30R and/or primer 16S1407F (Table 1) for the 16 clones with complete 16S rRNA gene sequences, as well as for nine closely related clones. Six clones with very similar ITS and different 16S rRNA genes were subjected to an additional control by sequencing ca. 750 bp with primer 16S1114F (Table 1).
DGGE analysis.
The crude DNA preparations were subjected to a purification step by using the Wizard DNA Clean-up system (Promega) and following the manufacturer's instructions.
16S rRNA gene fragments that were 422 bp long were generated by seminested PCR. The primers used for the first PCR were 16S378F and 23S30R (Table 1). The PCR conditions were similar to those used when the forward primer was 16S27F, but the reaction was carried out as follows: incubation for 5 min at 94°C, followed by 30 cycles of 45 s at 94°C, 45 s at 54°C, and 2 min at 68°C and then a final elongation step of 7 min at 68°C. The resulting PCR products (0.5 µl) served as templates for the second PCR, which was performed with forward primer 16S378F and reverse primers 16S781R(a) and 16S781R(b) (Table 1). A 38-nucleotide GC-rich sequence was attached to the 5' end of each of the reverse primers. The reaction conditions were the same as those described above except that amplification was carried out as follows: incubation for 5 min at 94°C, followed by 35 cycles of 1 min at 94°C, 1 min at 60°C, and 1 min at 68°C and then a final elongation step of 7 min at 68°C. Two distinct reactions were performed for each reverse primer. The negative control for the first PCR was used in the second PCR to check for contamination.
DGGE was carried out with a Dcode gene system (Bio-Rad) and was performed as described by Nübel et al. (40), with the following modifications. The PCR products obtained with primers 16S781R(a) and 16S781R(b) were applied separately to the polyacrylamide gel. The gel contained a linear 40 to 65% denaturant gradient, the pH of the TAE buffer was adjusted to 7.4, and electrophoresis was performed for 16 h at 45 V and 60°C.
The DGGE bands were excised with a surgical scalpel. Each small gel block was placed in 100 µl of sterile water for 2 h at room temperature. Each solution was used as a template for PCR amplification as described above. The PCR products were then electrophoresed under the DGGE conditions described above to confirm their positions relative to the bands from which they were excised and to detect potential heteroduplexes. A partial ca. 350-bp 16S rRNA gene sequence from primer 16S378F or 16S784R was determined.
Analysis of sequence data.
The sequences were initially analyzed by a similarity search by using the BLAST software, widely available on the Internet, and chimera detection was performed by using Check Chimera in the Ribosomal Database Project (30). The sequences determined in this study were included in the database of the ARB software package (29) available at http://www.arb-home.de and were aligned with the cyanobacterial sequences available from GenBank.
A distance tree was constructed with the software package TREECON for Windows 1.3b (60). The dissimilarity values were corrected for multiple substitutions by the method of Jukes and Cantor (25) and were used to calculate a distance matrix, and a tree was constructed by the neighbor-joining method (48). Aligned partial 16S rRNA gene sequences corresponding to E. coli sequence positions 405 to 780 were used, but the indels were not taken into account. A bootstrap analysis was performed that involved construction of 500 resampled trees. The tree comprised all sequences determined in this study (clones and DGGE bands) and all the sequences of polar cyanobacteria available in databases, as well as the related sequences (as indicated by BLAST) that exhibited more than 97.5% similarity with the sequences mentioned above. However, large groups of sequences that exhibited more than 97.5% similarity were represented by a limited number of selected sequences (e.g., Nostoc sp.). In addition, we included in the tree the two nearest neighbors indicated by BLAST when no sequences were related at a similarity level of more than 97.5%. If these two sequences were from uncultured clones, we added the sequences of the two closest cultured strains indicated by BLAST. In addition, we included at least one sequence for each of the clusters defined by Wilmotte and Herdmann (67).
To check the robustness of this analysis, other phylogenetic trees were constructed by neighbor-joining, maximum-likelihood, and maximum-parsimony methods by using the ARB software as described by Garcia-Pichel et al. (16).
Nucleotide sequence accession numbers.
New sequence data were deposited in the GenBank database. Sixteen almost complete 16S rRNA gene sequences were deposited under accession numbers AY151721 to AY151736, 33 partial 16S rRNA gene sequences were deposited under accession numbers AY151737, AY151738, and AY151740 to AY151770, and 24 ITS sequences were deposited under accession numbers AF547626 to AF547637, AF547639 to AF547649, and AY188584.
|
|
|---|
|
View this table: [in a new window] |
TABLE 2. Cyanobacterial morphotypes distinguished in the field sample and in the BGC sample
|
![]() View larger version (175K): [in a new window] |
FIG. 1. Diversity of cyanobacterial morphotypes identified in the artificial and natural microbial mats. (A) Nostoc sp.; (B) Hydrocoryne cf. spongiosa; (C) Nodularia cf. harveyana; (D) Leptolyngbya sp. 1; (E) Leptolyngbya sp. 2; (F) Phormidium cf. autumnale; (G) Oscillatoria cf. subproboscidea; (H) Schizothrix sp.
|
The eight morphotypes were observed in the field sample, but only three of them (morphotypes E, F, and H) were found in the BGC sample (Table 2). The major differences between samples were the absence of the three heterocystous taxa (Nostoc sp., Nodularia cf. harveyana, and Hydrocoryne cf. spongiosa) and the great abundance of Phormidium cf. autumnale in the BGC sample. Phormidium cf. autumnale, which was a minor component in the original sample, was the dominant morphotype in the BGC. This species is one of the most frequently reported taxa observed in Antarctic biotopes, and it comprises numerous morphotypes and ecotypes (27).
Clone libraries.
Altogether, 87 clones with an insert of the correct size were obtained; 45 of these clones produced 28 different ARDRA patterns for the field sample, and 42 clones produced 14 different ARDRA patterns for the BGC sample.
To assign clones to broad taxonomic clusters, 42 partial 16S rRNA gene sequences were determined for one representative of each ARDRA pattern. In addition, two cyanobacterial sequences obtained by Brambilla et al. (8) from a clone library of bacterial sequences for the field sample were added. Five sequences of chimeric origin, one plastid sequence, and three sequences resulting from cloning rearrangements were detected and excluded from the analysis. Finally, 35 partial cyanobacterial sequences were obtained.
To evaluate the levels of discrimination of the ARDRA with the restriction enzymes MboI and HpaII, partial 16S rRNA gene sequences were determined for nine additional clones selected at random from among seven ARDRA patterns represented by several clones. On average, 16S rRNA gene sequences (E. coli positions 405 to 780) of two clones with identical ARDRA patterns showed only a single substitution (i.e., 99.75% similarity). Hence, in the rest of our analysis we considered that clones with identical ARDRA patterns had virtually identical DNA partial sequences.
The 23 sequences (38 clones) obtained for the field sample and the 12 sequences (40 clones) obtained for the BGC sample were grouped in 15 phylotypes (Table 3) by using a threshold of 97.5% similarity (E. coli positions 405 to 780). Eleven phylotypes were detected only in the original sample and two phylotypes were detected only in the BGC sample, whereas two cyanobacterial phylotypes were present in both samples.
|
View this table: [in a new window] |
TABLE 3. Summary of the data obtained from the clone libraries and DGGE
|
The total number of phylotypes, the Shannon-Wiener diversity index, and the Berger-Parker dominance index (Table 4) showed that the genotypic diversity was twice as high in the field sample as in the BGC sample and that the BGC sample was dominated by a smaller subset of taxa.
|
View this table: [in a new window] |
TABLE 4. Richness. coverage, Shannon-Wiener, and Berger-Parker indices
|
For each of the 15 phylotypes, one or two complete 16S rRNA gene sequences were determined. Three complete sequences determined to be chimeras were discarded. This was also the case for clone Fr285, and phylotype 6 was therefore represented by only the partial sequence of Fr285 (E. coli positions 405 to 780). Hence, 16 complete 16S rRNA gene sequences that belonged to 14 phylotypes were obtained.
For the first time in a study of cyanobacterial diversity with environmental samples, ITS sequences were exploited to complement the information obtained from the 16S rRNA gene. ITS sequences of 26 clones were determined, including 2 sequences that were determined to be chimeras and discarded.
An alignment of the ITS sequences (http://www.ulg.ac.be/cingprot/ITS%20alignment.pdf) was constructed on the basis of conserved domains (24) and tRNAs. Groups of ITS sequences in which the alignment seemed meaningful (72) were distinguished, and we defined nine different ITS types (Table 3). ITS sequences belonging to the same ITS type exhibited between 94.3 and 100% similarity (average, 97.9%). Both tRNA genes were found in ITS types 1 to 5, as well as in the ITS sequences of clones Fr032 and Fr121, whereas the ITS sequence of clone Fr147 had only the tRNAIle gene and ITS type 6 did not have any tRNA genes. Table 3 shows the ITS types that were observed for the different phylotypes. For five phylotypes, the ITS sequences of several clones were obtained. Four of these phylotypes comprised a number of clones with different ITS types (Phy4/ITS4-5-6, Phy9/ITS1-2, Phy10/ITS1-3, Phy11/ITS2-3), while only a single ITS type (4) was found for phylotype 1. Unexpectedly, the same ITS types were found for clones that belonged to different phylotypes (ITS1/Phy9-10-11, ITS2/Phy9-11, ITS3/Phy10-11, ITS4/Phy1-4-13-15, ITS5/Phy1-3-4-8-10, ITS6/Phy4-5).
DGGE.
A total of 43 bands were observed in the DGGE patterns (data not shown). Theoretically (14, 45), if bands at the same position had identical sequences, these bands corresponded to 35 genotypes (21 for the field sample and 22 for the BGC sample, with 8 genotypes identical for the two samples).
Altogether, five DGGE band sequences (FrF1, FrF1bis, FrF3, FrF4, FrF5) belonging to four phylotypes were obtained for the field sample (Table 3). Two phylotypes (2 and 8) were found only in the clone library of the field sample, one phylotype (5) was found in both libraries, and interestingly, the last phylotype (9) was found in only the BGC clone library. Based on the number of clones in each phylotype (Table 3), we noticed that our DGGE band sequences did not belong to the phylotypes (1, 4, 10, and 11) which comprised the highest numbers of clones. No DGGE band sequences were obtained from the BGC patterns despite several attempts.
Phylogenetic analysis.
As shown in the distance tree in Fig. 2, all polar sequences (including the sequences obtained in this study) were distributed in 22 lineages; two of the lineages (clusters X and XIV) were novel, and seven (clusters VIII, IX, XI, XII, XV, XIX, and XXII) contained Antarctic representatives exclusively. Table 5 shows the origins of the polar sequences. The sequences obtained in this study were distributed in 11 lineages, which included the 2 novel lineages and 3 lineages that were exclusively composed of Antarctic representatives. The same novel and Antarctic clades were obtained with different methods of tree construction.
![]() View larger version (58K): [in a new window] |
FIG.2. Neighbor-joining tree based on partial 16S rRNA gene sequences corresponding to E. coli positions 405 to 780, with the exception of the sequences of the uncultured Antarctic cyanobacteria CLEAR-10, PENDANT-2, -4, and QSSC (E. coli positions 519 to 780). The tree includes the 49 sequences of clones and DGGE bands determined in the present study (boldface, italic, underlined type), 125 previously published sequences (Antarctic sequences are in boldface italic type, and Arctic sequences are in boldface roman type), and the E. coli sequence used as an outgroup. Bootstrap values equal to or greater than 70% are indicated at the nodes. The evolutionary distance between two sequences is obtained by adding the lengths of the horizontal branches connecting them and using the scale bar (0.1 mutation per position). Signatures are shown next to the cluster numbers. Abbreviations: Uncult., uncultured; cyanobact., cyanobacterium; Ant., Antarctic; str., strain.
|
|
View this table: [in a new window] |
TABLE 5. Origins of the polar sequences
|
Antarctic clusters.
Clusters IX, XI, and XII comprised sequences of clones isolated from the field sample from Lake Fryxell and other sequences of uncultured Antarctic cyanobacteria. Sequences belonging to these three clusters exhibited at least 8% dissimilarity with their closest relatives from nonpolar environments available in the databases. Therefore, we propose that these clusters are Antarctic lineages. The minimum level of sequence similarity within cluster IX was 91.5%, and sequences in this cluster belonged to phylotypes 7 and 15. Clusters VIII, XV, XIX, and XXII comprised only previously published Antarctic sequences. The closest nonpolar relatives exhibited at least 6% dissimilarity with the sequences belonging to these clusters. In addition, two previously described Antarctic strains, Cyanothece aeruginosa NIVA-CYA 258.2 and the LPP group cyanobacterium QSSC3cya, exhibited at least 9 and 12% dissimilarity, respectively, with all other database sequences.
Other clusters.
The rest of the sequences obtained in this study fell into clusters I, III, V, VII, XIII, and XXI. Six clone sequences from the BGC sample and one clone sequence from the field sample fell in cluster I and were grouped with polar and nonpolar sequences. Three of the Antarctic sequences were from the psychrophilic organisms Phormidium sp. strains Ant-Lunch and Ant-Orange and Oscillatoria sp. strain ANT-G16. Cluster III comprised two clone sequences from the field sample, as well as the sequence of the uncultured cyanobacterium cl lichen 4 (46). The levels of sequences similarity within this cluster ranged from 94.9 to 96.9%, and the clone sequences belonged to phylotypes 12 and 13. Furthermore, the most closely related strain exhibited at least 8% sequence dissimilarity. Interestingly, cluster V comprised five BGC clone sequences and only one sequence from the field sample obtained from a DGGE band. A single sequence from a nonpolar organism fell in this cluster; this was the sequence of Geitlerinema sp. strain PCC 9222 (35), which exhibited levels of similarity ranging from 98 to 99% with our sequences. Cluster VII corresponded to the order Nostocales. One clone sequence and two DGGE band sequences from the field sample fell in this cluster and were grouped together with other polar and nonpolar sequences of Nostoc strains. Cluster XIII contained one clone sequence from the field sample obtained with bacterial primers (8), the Arctic sequences of Phormidium sp. strain E18, and nonpolar sequences. Cluster XXI comprised the sequences of six BGC clones, the uncultured Antarctic clone LB3-1, the LPP group cyanobacterium QSSC8cya, and the nonpolar uncultured cyanobacterium TAF-B22 (41), which exhibited levels of similarity ranging from 96.5 to 97.8% with the Antarctic sequences of this cluster. In clusters II, IV, VI, XVI, XVII, XVIII, and XX, previously described Antarctic and/or Arctic sequences were grouped with nonpolar sequences.
Nadeau et al. (37) pointed out that there is a rare 11-nucleotide insertion (5'-AGTTGTGAAAG-3') in the 16S rRNA genes of the Antarctic isolates Ant-Lunch and Ant-Orange and the Arctic isolate NIVA-CYA 203. We found this insertion in the 16S rRNA genes of clones BGC-Fr025, BGC-Fr032, BGC-Fr068, BGC-Fr072, BGC-Fr078, and BGC-Fr080 in cluster I. In addition, the Microcoleus vaginatus PCC9802 sequence and many sequences recently deposited in the databases also have this insertion (7, 16, 46; D. A. Casamatta and M. L. Vis, unpublished data). We found that there was a good correlation between the highly variable region corresponding to E. coli positions 463 to 468, where this insertion was located, and the clusters of the tree (Fig. 2). Indeed, for 17 clusters, this region could serve as a signature, because all the sequences in these clusters were identical or nearly identical to each other at E. coli positions 463 to 468. However, the signatures of clones Fr121 and Fr127, which were distantly related to the other clones in cluster IX, were AGAAC and TAA, respectively. In addition, three different signatures were found among the sequences in cluster I (AAGTTGTGAAAGCAARCC, ACGCAN, and ACAAAA).
|
|
|---|
Based on morphological data, five of the eight species identified in the field sample by using morphological criteria (morphotypes) have been found in Canada Stream, which flows into Lake Fryxell close to our sampling site (9, 33). Even if a comparison based on species inventory had limits, these species were not observed in other mat samples or in the plankton of Lake Fryxell (34, 50, 53, 66).
When geographical distribution was considered, with the exception of Oscillatoria cf. subproboscidea, a species apparently endemic to Antarctica (27), the species found in this study have been found in various biotopes outside Antarctica and appear to have cosmopolitan distributions. Thus, on their own, the morphological results support the idea that endemism is rare among Antarctic cyanobacteria.
Molecular tools, however, revealed an entirely different pattern. Indeed, a higher number of phylotypes (15) than morphotypes (8) was detected. In addition, microdiversity, commonly found in molecular ecology studies (15), was observed. On this basis and considering the small size (0.5 g) of the mat sample studied, we assume that the cyanobacterial diversity in microbial mats of Lake Fryxell is quite high. Brambilla et al. (8) have described the enormous complexity of the eubacterial community of these mats on the basis of 16S rRNA gene sequences of clones and isolates.
The phylogenetic analysis showed that sequences determined in this study along with the polar sequences available from the databases were distributed in 22 lineages, 2 of which were novel and 7 of which were exclusively Antarctic. The sequences obtained in this study were distributed in 11 lineages, which included the 2 novel lineages and 3 lineages that are exclusively Antarctic. Furthermore, five Antarctic sequences determined in the present study and 13 sequences available in the databases were distantly related (less than 97.5% similarity) to the other sequences in their clusters. Therefore, if we considered numbers of phylotypes rather than numbers of lineages, our sequences belonged to 15 phylotypes, 4 of which were closely related to Antarctic sequences only and 7 of which were novel phylotypes. These results suggest that a previously undiscovered diversity was found and that Antarctic endemic species are more abundant than estimated previously on the basis of morphological features. In contrast, the three Nostoc sequences determined in this study belonged to cluster VII, which comprised other Nostoc sequences from polar and nonpolar environments. It seems, therefore, that the genus Nostoc has a more cosmopolitan distribution. However, these conclusions must be considered with care. Indeed, 16S rRNA gene sequence databases have grown very fast during the last decade, but they are certainly not exhaustive.
Nadeau et al. (37) suggested a bipolar distribution for several oscillatorians based on a rare 11-nucleotide insertion (5'-AGTTGTGAAAG-3') in the 16S rRNA genes of several Antarctic and Arctic strains. However, sequences from nonpolar environments also have this insertion, and a strict bipolar distribution was not observed in our analysis. In contrast to Antarctic sequences, Arctic sequences always clustered with nonpolar sequences. However, as only eight Arctic sequences were available in the databases, the apparent lack of Arctic endemism must be considered with care.
Interestingly, Gordon et al. (20) demonstrated, by using oligonucleotide probes specific for 16S rRNA gene sequences of uncultured Antarctic cyanobacteria isolated from the ice cover of Lake Bonney, that the lake ice microbial community was dominated by organisms that originated elsewhere in the surrounding region (e.g., terrestrial microbial mats). In agreement with this conclusion, some of our sequences clustered with sequences with which the probes were designed (clusters IX and XII).
Finally, as explained above, our findings support the conclusion that there is not a univocal relationship between morphological diversity and molecular diversity. The molecular results displayed greater diversity, and decisive arguments concerning the worldwide geographical distribution of cyanobacteria should therefore include genotypic data. However, as remarked by other authors, it cannot be presumed that the molecular approach provides a complete view of diversity. In the present study, for example, sequences related to the genera Nodularia and Hydrocoryne were not detected, although these genera were observed by light microscopy. To obtain as complete a picture of the diversity as possible, molecular methods and the more traditional methods based on microscopy should be combined.
Molecular tools: DGGE versus clone libraries, ITS sequences, and PCR primers.
The problem of DGGE band reamplification precluded a complete comparison of DGGE and clone libraries in terms of the sequences obtained. However, we observed that the clone libraries provided a more complete picture of the cyanobacterial diversity in Antarctic microbial mats than DGGE provided as they yielded a higher number of different sequences or phylotypes. In addition, Ferris et al. (14) indicated that the DGGE technique could be used to detect the most abundant organisms, but this was not the case in our study.
Since the ITS is less subject to selection pressure, it may accumulate more mutations than functional genes (39) and is therefore a variable marker. Consequently, it has been used successfully in several studies based on restriction fragment length polymorphism (49) or direct sequencing (3, 7, 13) to discriminate cyanobacterial strains at the intra- or interspecific level. In addition, it has proven to be useful in bacterial community studies (5). Here we evaluated the potential of its high variability, first in combination with the 16S rRNA gene to screen our clone libraries by the ARDRA technique and second to reveal finer-scale variations in the genotypic community composition by sequencing. However, the existence of multiple rRNA operons with ITS that may differ in sequence, length, secondary structure, and the presence or absence of tRNA genes in a single organism (7) must be taken into consideration when data are interpreted.
Our analysis showed that different ITS types, which could not be aligned meaningfully, were retrieved from clones belonging to the same phylotype (16S rRNA gene partial sequences that were up to >99.5% similar). In addition, we did not expect to find clones belonging to different phylotypes with identical ITS types (levels of sequence similarity ranging from 94.3 to 100%). As we resequenced the junctions between the 16S rRNA gene and the ITS, this observation could be explained only by the existence of lateral gene transfer or by the presence of undetected chimeras among our clones. The latter hypothesis is more credible as the formation of chimeras is a well-known artifact of PCR (45, 54), and we excluded several such artifacts during this study.
For the first time in community analysis, a PCR primer pair that allowed amplification of the 16S rRNA gene and the ITS that is specific for the cyanobacterial sequences was used. Primer 16S27F is a universal primer for bacteria. Primer 23S30R and its previous version (primer 23S26R) have been used with success with a wide variety of strains since 1993 (28, 38, 70, 71, 72). The high number of cyanobacterial sequences and the absence of other bacterial sequences in our clone libraries confirm the applicability of this primer combination, which performs better than the universal bacterial primers. Indeed, with the latter primers, only five cyanobacterial sequences (including two chimeras) and one plastid sequence were obtained among a total of 325 sequences for a subsample of the same microbial mat from Lake Fryxell (8).
Comparison of the field mat and the cultured artificial mat, linking the community structure of the artificial mat with its functional properties.
Our analysis of the cyanobacterial diversity in a field mat sample from Lake Fryxell was enlarged to include a derived cultured mat from a BGC, since the BGC mat has been extensively studied by using microsensor techniques and pigment analysis (11, 43), which permit linkage of its community structure directly to its functional properties.
In terms of pigment content and primary productivity, Buffan-Dubau et al. (11) have suggested that artificial mats have some characteristics of natural polar mats, although some differences were found.
Not unexpectedly, the present study showed that the BGC culture conditions selected a cyanobacterial mat community which represented some of the diversity observed in the original field mat used as the inoculum. With the exception of one phylotype, the morphotypes and phylotypes retrieved from the BGC sample were present in the field sample. The most notable differences between the two mats were the high number of sequences closely related to psychrophilic oscillatorians, the dominance of Phormidium cf. autumnale, and the apparent absence of Nostocales in the BGC sample. Based on the 16S rRNA gene, psychrophilic oscillatorian strains appear to be restricted to a single monophyletic clade (37), while most polar cyanobacteria are psychrotolerant (58) and the order Oscillatoriales is polyphyletic (23, 38, 67, 69). This suggests that the BGC clones belonging to cluster I represented possibly psychrophilic cyanobacteria. The numerical importance of these clones may explain the observation that oxygenic photosynthesis in the artificial cold-adapted mats was optimal at 10°C, while no oxygenic photosynthesis was detected at temperatures above 15°C. The high concentration of hydrogen sulfide, which was more than 100 µM at the end of the 6-h dark period below a depth of 1 mm (43), along with the presence of oscillatorians able to migrate in the mats depending on the light conditions (4, 36), suggests that Phormidium cf. autumnale (morphotype F) migrated to the surface when the H2S concentration was too high and was therefore selected in the BGC. The apparent absence of heterocystous taxa was also reflected by the scytonemin contents, which were 2 orders of magnitudes lower in the BGC mat than in the field mat. Scytonemin plays an important role in photoprotection against the deleterious effects of UV radiation in Antarctic mats (21) but does not provide a real advantage under culture conditions where artificial light is less harmful (11). This may suggest that members of the Nostocales were overgrown by fast-growing oscillatorians. Most likely, the absence of heterocystous cyanobacteria in the BGC was related to the availability of high levels of bound nitrogen (NH4Cl), as well as the greater sensitivity of the organisms to high sulfide concentrations (22).
This study provided genotypic confirmation of the existence of endemic cyanobacteria in Antarctica, as well as the existence of cosmopolitan taxa that have Antarctic representatives. Additional studies of cyanobacteria living in microbial mats in different regions of Antarctica would be interesting to verify the geographical conclusions of this study and to obtain more information about the ecological range of the cyanobacterial genotypes. Other techniques, like fluorescent in situ hybridization, could also help workers assign morphotypes to clusters containing only sequences of uncultured organisms.
We thank Cathy Welsch, who collected the mat material, and the Long Term Ecosystem Research Program, under whose auspices the material was collected. Johanna Laybourn-Parry (University of Nottingham, Nottingham, United Kingdom) organized the transport of samples. Dominic Hodgson (British Antarctic Survey) and Warwick Vincent (Laval University, Quebec, Canada) are acknowledged for stimulating discussions, for reading versions of the manuscript, and for contributing many helpful comments and suggestions. We thank the anonymous reviewers for their corrections.
|
|
|---|
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»