Previous Article | Next Article ![]()
Applied and Environmental Microbiology, February 2003, p. 1030-1042, Vol. 69, No. 2
0099-2240/03/$08.00+0 DOI: 10.1128/AEM.69.2.1030-1042.2003
Copyright © 2003, American Society for Microbiology. All Rights Reserved.
Department of Marine Sciences, University of Georgia, Athens, Georgia 30602-3636,1 Agricultural Research Service, U.S. Department of Agriculture, Watkinsville, Georgia 306772
Received 27 June 2002/ Accepted 14 November 2002
|
|
|---|
- and
-Proteobacteria (6 and 10%, respectively), Cytophaga-Flexibacter-Bacteroides (19%), high-G+C-content gram-positive organisms (Actinobacteria; 25%), and low-G+C-content gram-positive organisms (Bacillus and Clostridium; 19%). Twelve percent were identified as chloroplasts. The remaining 9% represented ß- and
-Proteobacteria, Verrucomicrobiales, and candidate divisions. Mixolimnion and oxycline samples had low microbial diversity, with only 9 and 12 distinct phylotypes, respectively, whereas chemocline and monimolimnion samples were more diverse, containing 27 and 25 phylotypes, respectively. The compositions of microbial assemblages from the mixolimnion and oxycline were not significantly different from each other (P = 0.314 and 0.877), but they were significantly different from those of chemocline and monimolimnion assemblages (P < 0.001), and the compositions of chemocline and monimolimnion assemblages were not significantly different from each other (P = 0.006 and 0.124). The populations of sequences retrieved from the mixolimnion and oxycline samples were dominated by sequences related to high-G+C-content gram-positive bacteria (49 and 63%, respectively) distributed in only three distinct phylotypes, while the population of sequences retrieved from the monimolimnion sample was dominated (52%) by sequences related to low-G+C-content gram-positive bacteria distributed in 12 distinct phylotypes. Twelve and 28% of the sequences retrieved from the chemocline sample were also found in the mixolimnion and monimolimnion samples, respectively. None of the sequences retrieved from the monimolimnion sample were found in the mixolimnion or oxycline samples. Elevated diversity in anoxic bottom water samples relative to oxic surface water samples suggests a greater opportunity for niche differentiation in bottom versus surface waters of this lake. |
|
|---|
With prolonged meromixis, anoxic bottom water at a high pH can attain high levels of toxic inorganic compounds, including sulfide and ammonia (33). These conditions have a profound effect on the ecology of the lake (30, 44) and may be expected to affect the composition and activity of the microbial community of the lake. Stratification affects the flux of reduced metabolites into the surface layer from bottom water and lake bottom sediments. Reduced substances accumulate to high concentrations in bottom water, and biota in surface layers may be deprived of essential nutrients. High concentrations of biologically generated reduced substances (sulfide, for example) may directly affect geochemical processes (i.e., metal precipitation or speciation) and thus indirectly affect biota by creating secondary nutrient limitations. Sulfide toxicity has also been shown to affect microbial processes directlyfor example, denitrification (31) and ammonia oxidation (34). Sulfide toxicity may result in the exclusion of zooplankton grazers from anoxic or hypoxic regions of the water column. In Mono Lake, this situation leads to the accumulation of bacterial biomass at depth, to an increase in average cell size (large bacterioplankton are observed in surface water before the annual bloom of the brine shrimp Artemia monica, the dominant zooplankter in Mono Lake), and to the formation of a deep chlorophyll maximum at the base of the oxycline.
Compared to more thoroughly investigated marine and freshwater environments, relatively little is known about the microbial ecology of saline lakes. Interest in soda lake microbiology has centered primarily on the isolation and characterization of individual microorganisms with potential industrial applications (32). Little attention has been paid to the community as a whole. A recent study examined the phylogenetic diversity of bacteria isolated from soda lakes (15). However, little work has been done on the in situ phylogenetic diversity of bacterial populations in soda lakes or on ascribing aspects of the biogeochemical cycles of soda lakes to these bacteria. Since soda lakes are environments of extreme physiochemical conditions, organisms living in these lakes may possess novel adaptations to these environmentsfor example, transport systems, osmoregulatory compounds, or ectoenzymes. Some of these adaptations may involve enzymes or other physiological properties with potential commercial significance (24, 28).
Soda lakes are also anomalous environments in that they have high concentrations of dissolved organic carbon (DOC), dissolved inorganic phosphorus, and high standing crops of bacteria. These characteristics suggest that the microbial loop processes that have been studied more thoroughly for marine and freshwater environments (3, 4) are somehow decoupled in soda lakes, allowing DOC to accumulate. Understanding the reasons for this situation will help improve understanding of the consequences of the microbial loop for the ecology and biogeochemistry of other aquatic environments and provide insight into the processes that contribute to the storage and deposition of organic carbon (e.g., see references 14 and 45).
In a previous study (27), denaturing gradient gel electrophoresis (DGGE) was used to investigate seasonal changes in the vertical structure of Mono Lake microbial communities during a period of holomixis. This study indicated that there were significant differences in the compositions of bacterial assemblages with depth under stratified conditions but did not provide a detailed description of the phylogenetic compositions of these communities. The study presented here was conducted to compare the phylogenetic compositions and diversity of bacterial assemblages at four depths chosen to represent a range of water column redox conditions. Compared to samples from other aquatic environments, the Mono Lake microbial assemblages that we sampled displayed low diversity. We found distinctly different patterns of diversity in oxic versus anoxic waters of the lake, with surface water microbial communities being characterized by lower microbial diversity dominated by high-G+C-content gram-positive bacteria (Actinobacteria), while bottom water microbial communities had higher diversity dominated by low-G+C-content gram-positive bacteria (Bacillus and Clostridium).
|
|
|---|
sensor), turbidity (WetLabs C-Star transmissometer, 10-cm path length), and in vivo fluorescence (WetLabs WetStar fluorometer) were obtained with a SeaBird SeaCat profiler. Oxygen profiles were obtained with a polarographic oxygen sensor (YSI) equipped with a Clarke-type electrode.
Bacterial counts.
Water (10 ml) for bacterioplankton enumeration was preserved with buffered formaldehyde at a 2% final concentration. Samples were stored at 4°C in the dark for no more than 1 month until enumeration was done. Bacterial cells in 0.1 to 1 ml of lake water were mixed with filtered lake water as needed to dilute and disperse cells; the mixtures were filtered onto black 0.2-µm-pore-size polycarbonate membrane filters, rinsed once with 2 ml of 10% acetic acid, stained with 2 ml of 4',6'-diamidino-2-phenylindole (DAPI) (1 µg ml-1 final concentration), and enumerated by epifluorescence microscopy (54).
DNA extraction.
Water samples (125 to 900 ml; see Table 2) were pressure (ca. 50 kPa) filtered through Millipore Sterivex filter cartridges (0.22-µm-pore size) to collect microbial biomass for subsequent DNA extraction. Excess water was expelled. The cartridges then were filled (1.8 ml) with a buffer containing 50 mM Tris (pH 8.3), 40 mM EDTA, and 0.75 M sucrose; capped; frozen on dry ice; shipped to the laboratory; and stored at -70°C until processed. DNA extraction and purification were carried out essentially as described by Ferrari and Hollibaugh (19). Briefly, 40 µl of lysozyme (50 mg ml-1) was added to each cartridge, and the cartridges were incubated for 60 min at 37°C. Fifty microliters of proteinase K (20 mg ml-1) and 100 µl of a 20% (wt/vol) solution of sodium dodecyl sulfate were added to each cartridge, and the cartridges were incubated at 55°C for 2 h. DNA was purified from 800 µl of the lysate by sequential extraction with 800 µl of phenol-chloroform-isoamyl alcohol (25:24:1), chloroform-isoamyl alcohol (24:1), and finally n-butanol. The aqueous phase was removed, placed in a Centricon-100 concentrator (Amicon), mixed with 500 µl of a buffer containing 10 mM Tris and 1 mM EDTA (pH 8.0), and centrifuged at 1,000 x g for 10 min. Next, 500 µl of TE was added to the Centricon-100 concentrator, and the mixture was centrifuged for another 10 min.
|
View this table: [in a new window] |
TABLE 2. Comparison of the compositions of libraries from samples collected at different depths in Mono Lakea
|
Clone libraries.
The PCR product (50 ng) was ligated into pGEM-T Easy Vector (Promega) and transformed into competent E. coli JM109 cells. The transformed cells were plated on Luria-Bertani plates containing 100 µg of ampicillin ml-1, 80 µg of 5-bromo-4-chloro-3-indolyl-ß-D-galactopyranoside (X-Gal) ml-1, and 0.5 mM isopropyl-ß-D-thiogalactopyranoside (IPTG) as recommended by the manufacturer and incubated overnight at 37°C.
Sequencing.
Clones containing an insert were selected randomly from each library for sequencing. We used only one library from each sample for the analysis presented here; however, 10 other libraries were constructed from different samples and used to identify DGGE bands. Although clones from these other libraries were selected for sequencing because they gave DGGE bands that matched bands of interest from the samples, rather than being selected randomly, they contained the same sequences with approximately the same frequency distribution as that which we determined by random sampling.
All sequences were obtained with an automatic sequencer operated at the Molecular Genetics Instrumentation Facility. A total of 239 clones (55 from 2 m, 61 from 17.5 m, 57 from 23 m, and 66 from 35 m) were selected, purified, and screened by sequencing in one direction through a hypervariable region (V3) of the 16S rRNA gene by using primer 356F (forward, 5'-CCTACGGGAGGCAGCAG-3'), yielding readable sequences of
600 bp. Ten of the sequences were ambiguous and were discarded. The remaining sequences were tested for chimeras by using the Ribosomal Database Project CHIMERA-CHECK program. Chimeras were also detected by generating phylogenetic trees with different regions of the sequences. Seventeen of these clones were identified as containing chimeric sequences and were excluded from further analysis. The remaining 212 cloned sequences (51 from the oxic mixolimnion, 2 m; 59 from the base of the oxycline, 17.5 m; 46 from the top of the chemocline, 23 m; and 56 from the monimolimnion, 35 m) were used in phylogenetic analyses. Forty of the these clones, representing all of the major phylotypes, were also sequenced with plasmid primers Sp6 and T7 to obtain nearly complete sequences of the 16S rRNA gene.
Phylogenetic analyses.
Sequences were assigned to major groups, and reference sequences were chosen, based on BLAST (1) similarities. For further phylogenetic analyses, sequences were aligned by using the PILEUP tool of the Wisconsin package, version 10.0 (Genetics Computer Group Inc.). Phylogenetic trees were constructed by using positions 411 to 854 (Proteobacteria), 491 to 882 (high-G+C-content gram-positive organisms), 541 to 972 (low-G+C-content gram-positive organisms), and 540 to 993 (Cytophaga-Flexibacter-Bacteroides [CFB]). Phylogenetic trees were inferred, and bootstrap analysis (100 replicates) was performed with the PHYLIP package (17) by using evolutionary distances (i.e., Jukes-Cantor distances) and the neighbor-joining method. Another phylogenetic tree was constructed by using the nearly full-length sequences (positions 110 to 1404). Sequence similarity values cited below are based on GAP analysis (Wisconsin package, version 10.0) of the nearly full-length sequences representing each phylotype, when they were available, or the shorter, "screening" sequences when they were not.
Statistical analysis and sequence population diversity.
For the purposes of assigning sequences to distinct phylotypes, sequences with similarities of greater than 97% were considered to represent the same phylotypes. In most studies in which similar analyses were performed, the discriminator values used were greater (i.e., 98 to 99% similarity) (36, 47, 59). We chose to use a relatively conservative value to be inclusive (i.e., we tended to "lump" rather than "split") in consideration of possible errors in the sequences due to PCR (70) and cloning (51, 65). We also wanted to avoid focusing on microdiversity (11, 20). As a result, the taxa (phylotypes) that we define for the purposes of this analysis may not represent individual "species" of bacteria, although 97% 16S rRNA gene sequence similarity was proposed previously as a discriminator of species versus strain differences (25, 66).
We then assumed that the relative abundance of phylotypes defined as described above reflects the relative abundance of the bacteria possessing those sequences in the population of bacteria in our sample. It is not necessary to make this assumption for the purposes of the subsequent analyses, but the analyses are devoid of ecological significance without it. We recognize that PCR biases, differential extraction efficiencies, multiple copies of the rrn operon per genome, and similar artifacts (16, 56, 67, 70) may affect the distribution of sequences in our libraries. We also recognize the possibility that not all of the bacteria present in Mono Lake when we sampled are represented in our libraries, as highly divergent sequences may not have been amplified by the primers used. While the assumption made above is a major one, analysis of the same and similar samples by DGGE with a different set of PCR primers produced qualitatively similar results (27; J. T. Hollibaugh, unpublished data), and the distributions determined by direct counting of readily identifiable organisms (a Synechococcus-like organism and the chloroplast of Picocystis salinarum) agreed well with the distributions obtained from analyses of the clone libraries.
Coverage (C) was calculated by using the following formula: C = 1 - (n1/N), where n1 is the number of phylotypes that occurred only once in the clone library and N is the total number of clones analyzed (48). Rarefaction curves (26) were produced by using software available online at http://www.uga.edu/
strata/software.html. The phylogenetic compositions of libraries were compared by using the Sorensen similarity index, Cs = 2j/(a + b), where j is the number of phylotypes common to both samples and a and b are the numbers of phylotypes in libraries A and B, respectively (43). The diversity of the libraries and the total number of species in our samples were estimated by using the approach presented by Curtis et al. (10), which assumes a lognormal distribution of species abundance and uses the relative abundance of phylotypes in clone libraries and the total abundance of organisms in the samples (here, DAPI abundance estimates times sample volume) to estimate the total number of phylotypes present. Sequences identified as originating from chloroplast 16S rRNA genes were not included in any of these analyses.
The statistical significance of differences in the compositions of pairs of libraries was tested by using the LIBSHUFF program (60). Sequences (
500 bp) were aligned as described above, and evolutionary distances were calculated by using the Jukes-Cantor algorithm. The distance matrix was then used to determine the significance of differences in the compositions of pairs of libraries.
Nucleotide sequence accession numbers.
The sequences of the cloned inserts were deposited in GenBank (see Fig. 2 for accession numbers). Only one representative of sequences that were more than 99% similar was deposited.
![]() ![]() ![]() ![]() View larger version (161K): [in a new window] |
FIG. 2. Neighbor-joining trees (partial sequences, 400 bp) showing phylogenetic relationships of 16S rRNA gene sequences cloned from Mono Lake to closely related sequences from GenBank. Clones from this study are coded as follows, with ML623J-31 as an example: ML, Mono Lake; 6, station 6; 23, sample depth in meters; J, month in which samples were collected (July); and 31, number assigned to the clone in that library (thus, ML623J is the library identifier). Scale bars indicate Jukes-Cantor distances. Bootstrap values of >50% (for 100 iterations) are shown. The trees are unrooted, with Halobacterium salinarum as the outgroup. GenBank accession numbers are given in parentheses. pos, positive.
|
|
|
|---|
![]() View larger version (27K): [in a new window] |
FIG. 1. Depth distributions of selected physicochemical and biological variables in the Mono Lake water column at station 6 on 20 July 2000. (A) PAR, oxygen, temperature, and excess density. µE, microeinsteins. (B) Ammonia (NH3) and sulfide (S2-). (C) Bacterial abundance, DOC, fluorescence, and chlorophyll. Note that some variable values are scaled as indicated (by numbers following slashes) to fit axis ranges.
|
(i) Chloroplast-like sequences
A total of 18 (29%) of the sequences retrieved from the 17.5-m sample, 5 (11%) retrieved from the 23-m sample, and 1 each (2%) retrieved from the 2- and 35-m samples were from the chloroplast 16S rRNA gene of the chlorophyte alga P. salinarum. Two other sequences retrieved from the libraries of the 2- and 17.5-m samples were related (96% similarity) to the chloroplast 16S rRNA gene of the diatom Odontella sinensis.
(ii)
-Proteobacteria.
Eleven of the 13 sequences that fell into the
-Proteobacteria group clustered together (Fig. 2A and 3) and were found in libraries generated from the 2-, 17.5-, and 23-m samples. The nearly full-length sequence (1,391 bp) of one of these clones (ML602J-43) was 97 to 98% similar to Rhodobaca bogoriensis (46), Roseonatronobacter thiooxidans (64), and Natronohydrobacter thiooxidans. It was also 97% similar to a bacterium (ML-168a) previously isolated from Mono Lake sediment (GenBank description). R. bogoriensis is an anoxygenic, phototrophic, alkaliphilic purple nonsulfur bacterium isolated from soda lakes in the African Rift Valley (GenBank description). It is capable of both phototrophic growth and chemoautotrophic growth on a wide variety of organic compounds (46). R. thiooxidans is a strictly aerobic and obligately heterotrophic alkaliphile isolated from Russian soda lakes (64). An identical sequence was found in a DGGE band in a previous study of Mono Lake (27). The only
-Proteobacteria sequence obtained from the monimolimnion was related (96% similarity over 1,385 bp) to a Roseobacter sp.
![]() View larger version (46K): [in a new window] |
FIG. 3. Neighbor-joining tree of nearly full-length sequences ( 1,300 bp) of representative examples from each of the major clades from Fig. 2 showing phylogenetic relationships to 16S rRNA gene sequences from GenBank. The scale bar indicates Jukes-Cantor distances. Bootstrap values of >50% (for 100 iterations) are shown. The trees are unrooted, with Halobacterium salinarum as the outgroup. GenBank accession numbers are given in parentheses. pos, positive.
|
(iv)
-Proteobacteria.
Twenty-one sequences (only 1 from the monimolimnion) clustered with the
-Proteobacteria (Fig. 2A and 3). Thirteen of these sequences (9 from the mixolimnion) clustered together (cluster 1). The nearly full-length sequence (1,462 bp) of one of these clones (ML602J-47) was distantly related (93 to 94% similarity) to Alcalilimnicola halodurans, Alkalispirilum mobile, and Nitrococcus mobilis. It was also 92% similar to a bacterium (MLHE-1) (49) isolated from Mono Lake. A. halodurans is an alkaliphilic, extremely halotolerant bacterium isolated from the sediment of soda-depositing Lake Natron, Tanzania (72). A. mobile is a moderately halophilic, alkaliphilic, gram-negative, obligate aerobe (58). Previous studies of Mono Lake reported sequences that were 99% similar to these sequences in oxycline samples (27, 68). Two sequences (one from the chemocline and one from the monimolimnion) fell into cluster 2 and were related (96% similarity) to an obligately autotrophic, sulfur-oxidizing, alkaliphilic bacterium (Arh1) of the genus Thioalkalivibrio. This organism was reported to use thiocyanate as a sole source of energy under highly alkaline conditions (63). A sequence found only in the clone library from the chemocline sample was related (98% similarity) (Fig. 3) to Thioalcalivibrio jannaschii, a haloalkaliphilic, obligately chemolithoautotrophic, sulfur-oxidizing bacterium recently isolated from Mono Lake (61). Three additional sequences related most closely to chemoautotrophic sulfur oxidizers (Thiomicrospira thyasirae, 98% similarity; and Thioalcalomicrobium sibericum, 96% similarity) (Fig. 2A) were recovered from the oxycline sample. T. sibericum is a member of a newly described genus isolated from soda lakes in Siberia and Kenya (62).
(v)
-Proteobacteria.
One sequence from the chemocline clone library was distantly related (90% similarity) to a clone designated GOUTA17 (Fig. 2A) and obtained from groundwater samples contaminated with chlorobenzene (GenBank description) as well as to sulfate-reducing bacteria (<90% similarity) obtained from salt marsh sediments and deep-sea environments. A sequence from the monimolimnion library was related (95% similarity) (Fig. 2A) to Desulfonatronovibrio hydrogenovorans, an alkaliphilic, sulfate-reducing bacterium (75).
(vi) High-G+C-content gram-positive bacteria (Actinobacteria).
A total of 53 sequences clustered with the high-G+C-content gram-positive bacteria (Fig. 2B and 3). These sequences were associated with the oxygenated surface water of Mono Lake, where they were the dominant sequences in clone libraries: only three of these sequences were obtained from the chemocline library, and none were obtained from the monimolimnion library. Twenty-three sequences formed a separate lineage (cluster 1) characterized by only one previously described sequence, from an isolate identified as Chibacore 1500 (GenBank accession no. AB038407; 94% similarity). Chibacore 1500 is an alkaliphilic bacterium isolated from a depth of 1,500 m and showing only sparse growth on any culture medium; it was subsequently lost (GenBank description; H. Kaneko, personal communication). Another 27 sequences (26 from the mixolimnion and oxycline libraries) fell into a second cluster and were 94 to 96% similar (Fig. 3) to Rathayibacter tritici and Microbacterium arborescens and to an environmental isolate, FukuN101, cloned from a sample collected from an acidic forest lake in Germany (23). Three sequences that fell into a third cluster were distantly related (93% similarity) to a clone designated HTG5 and obtained from microbial populations associated with Metallogenium morphotypes and amorphous particles in a freshwater reservoir (GenBank description).
(vii) CFB group.
Forty sequences (7 from the mixolimnion, 6 from the oxycline, 10 from the chemocline, and 17 from the monimolimnion libraries) fell into the CFB clade (Fig. 2C and 3). Fifteen monimolimnion and 5 chemocline sequences formed a separate lineage (cluster 1) and were not affiliated with any previously documented sequences. Two other sequences from the monimolimnion library were distantly related (92%) to the partial sequence (
900 bp) of a clone designated BURTON-42 and obtained from the anoxic sediment of a meromictic lake in eastern Antarctica (5). Three deeply branching sequences from the chemocline sample (ML623J-25, -29, and -58) were distantly related (89% similarity over positions 514 to 1236) to the partial sequence of clone kpc103f (GenBank accession no. AF195431) (not shown in Fig. 2C and 3), found to co-occur with the brown-tide phytoplankter Aureococcus anophagefferens (GenBank description). Eight sequences from the mixolimnion and oxycline libraries that clustered together were related (95% similarity) to the partial sequence (
1,000 bp) of a clone, designated MPD-32, obtained from a selenium-contaminated, hypersaline evaporation pond (12). They were 90% similar to the partial sequence (
900 bp) of a clone, designated MT70, obtained from hydrogen sulfide-rich, black mud from a coastal marine environment (GenBank description). No other previously documented sequences fell into this cluster. Three sequences from the mixolimnion library were distantly related (86 to 88% similarity) to the sequence of a clone designated WCHB1-53 and obtained from a hydrocarbon- and chlorinated-solvent-contaminated aquifer (13); to a clone designated MERTZ_21CM_238 and obtained from continental shelf sediments off the coast of Antarctica (GenBank description); and to the 16S rRNA gene sequence of bacterium JTB132 (GenBank accession no. AB015260) (not shown in Fig. 2C and 3), isolated from a cold seep in the Japan Trench (40). A previous study (27) of Mono Lake reported the presence of a similar sequence in mixolimnion samples.
(viii) Low-G+C-content gram-positive bacteria (Bacillus and Clostridium spp.).
Low-G+C-content gram-positive bacterial sequences dominated libraries from the anoxic water of Mono Lake. Forty sequences (2 from the mixolimnion, 9 from the chemocline, and 29 from the monimolimnion) fell into the low-G+C-content gram-positive clade (Fig. 2D and 3). None of these sequences was closely related to any 16S rRNA gene sequence in the GenBank database. Twenty-seven sequences (19 from the monimolimnion and 8 from the chemocline) formed a separate lineage (cluster 1). Two sequences were distantly related (86%) to Anaerobranca bogoriae, a recently described thermoalkaliphile (55). This organism is an obligate anaerobe capable of growth on a variety of carbohydrates and does not use sulfate as an electron acceptor. Two sequences were distantly related (85%) to Natronoanaerobium salstagnum, an obligately anaerobic alkaliphile isolated from salterns of alkaline, hypersaline Lake Magadi, Kenya (32). The only two sequences obtained from oxic water (mixolimnion library) and related to the low-G+C-content gram-positive bacteria were distantly related (83%) to Enterococcus spp., facultative anaerobes that are part of the normal gut microflora, and thus might have originated from fecal material produced by the numerous birds that feed on the brine shrimp population of Mono Lake during the summer.
(ix) Spirochaetales.
Sequences of three clones obtained from the chemocline library (Fig. 2D) were related to moderately halophilic, alkaliphilic Spirochaeta spp. isolated from various soda lakes (74), including Mono Lake (GenBank description).
(x) Others.
We also retrieved eight sequences from the chemocline and monimolimnion libraries that clustered with the Verrucomicrobiales (Fig. 2D). Three sequences clustered with candidate divisions (Fig. 2D). A short sequence obtained from the monimolimnion library was related (97% similarity) to the marine cyanobacteria Synechococcus and Prochlorococcus (data not shown).
Distribution of bacterial groups.
The relative abundances of different phylogenetic groups in the libraries (excluding chloroplast sequences) were calculated for all four samples (Fig. 4). The phylogenetic compositions of our clone libraries changed completely with depth. A total of 49 and 63% of the sequences retrieved from the mixolimnion (2-m) and oxycline (17.5-m) libraries, respectively, were affiliated with high-G+C-content gram-positive bacteria, whereas 52% of the sequences retrieved from the monimolimnion (35-m) library were affiliated with the low-G+C-content gram-positive group. A total of 22% of the sequences retrieved from the chemocline (23-m) library were affiliated with low-G+C-content gram-positive bacteria, and 7% were affiliated with high-G+C-content gram-positive bacteria. Sequences affiliated with the CFB group were less abundant in the mixolimnion library than in the chemocline and monimolimnion libraries (14% versus 24 and 31%, respectively) and, as mentioned above, the mixolimnion and monimolimnion libraries were dominated by different CFB.
![]() View larger version (75K): [in a new window] |
FIG. 4. Frequencies of phylotypes affiliated with major phylogenetic groups in libraries from 2-, 17.5-, 23-, and 35-m samples. Chloroplast-like sequences were omitted from the data set used to calculate percentages. , , and , Proteobacteria; Verruco, Verrucomicrobiales; Spiro, Spirochaetales.
|
|
View this table: [in a new window] |
TABLE 1. Properties of the distributions of phylotypes in clone libraries from Mono Lake samplesa
|
![]() View larger version (17K): [in a new window] |
FIG. 5. Rarefaction curves generated for 16S rRNA genes in clone libraries from samples collected at 2, 17.5, 23, and 35 m.
|
|
|
|---|
Even more striking was the dramatic difference in the pattern of diversity in oxic versus anoxic samples (Fig. 2 and 3). Relative to what has been found for other environments that have been characterized (7, 21, 36, 42, 48), the diversity of Bacteria in Mono Lake is low, as indicated by the coverage, rarefaction, and diversity analyses that we performed. However, within the lake, diversity was higher, with more deeply diverging taxa in anoxic deep water than in oxic surface water. The mixolimnion and oxycline libraries contained 9 and12 phylotypes in four and five clades that accounted for 78 and 85%, respectively, of all of the sequences that we retrieved from these libraries. In contrast, the chemocline and monimolimnion libraries contained 25 and 27 phylotypes in five and six major clades that accounted for only 44 and 48%, respectively, of all of the sequences retrieved from these libraries. This difference is also reflected in the coverage estimates for these libraries, the estimated total number of phylotypes in the samples, and their diversity indices (Table 1) and by their rarefaction curves (Fig. 5).
Note that the oxycline sample was identified as the least diverse by the analytical approach described by Curtis et al. (10), whereas the mixolimnion sample was identified as the least diverse by coverage estimates and rarefaction curves. Curtis et al. (10) assume a lognormal distribution of species abundance. We tested this assumption with our data sets and found that it was met (P < 0.05) only by data from the monimolimnion sample. However, the data distributions were highly skewed to the right tail of the curves (toward the most abundant phylotypes), rather than being distributed over the whole range of species abundances (i.e., rare phylotypes were not enumerated), greatly reducing the power of the statistical test of lognormality. As pointed out by Curtis et al. (10), correcting this bias would involve sequencing of a much larger number of clones from each library.
Our results are in contrast to those of a previous study of Mono Lake microbial assemblages based on DGGE analysis (27). That study concluded that diversity was higher in oxic surface water samples than in anoxic deep water samples. A similar DGGE-based study of another stratified lake (50) also reported higher diversity in oxic surface water than in anoxic deep water. A likely explanation for this discrepancy, at least with regard to Mono Lake, is that most of the phylotypes found in deep water were present at a very low relative abundance and thus formed faint bands that were not readily detected by DGGE. This explanation is supported by DGGE analysis of these samples (Hollibaugh, unpublished) by using technology superior to that available to Hollibaugh et al. (27).
Although it is an anthropocentric view, the bottom water of Mono Lake might be considered to be more extreme than the surface water because of the potential toxicity of the high concentrations of sulfide and ammonia that have accumulated after 8 years of meromixis. Current ecological theory would lead to the expectation that the species (here, sequence) diversity of populations in the extreme environment of bottom water would be low relative to the diversity of populations from surface water, where less extreme conditions prevail, yet we found the opposite to be true.
We considered the possibility that heterogeneity of the 16S rRNA gene in different rrn operons of the same genome might contribute to an artifactually high apparent diversity of Mono Lake bottom water (16, 35), since our analyses indicated that these samples were dominated by low-G+C-content gram-positive bacteria related to Bacillus and Clostridium, which are known to contain multiple copies of rrn operons. We tried to minimize the effect of intragenomic heterogeneity of rrn operons on our analysis by using a low similarity value (97%) to decide whether or not two sequences were the same.
If intragenomic heterogeneity is not the cause of the high diversity in the chemocline and monimolimnion samples, then environmental factors must be contributing to the diversity that we found. At least for the monimolimnion sample, it is unlikely that this diversity can be ascribed to mixing of disparate populations because these waters have been isolated for several years. Sinking Artemia fecal pellets may also have transported surface water microorganisms into deep water. This process may account for the presence in the chemocline library (and thus partly for the greater diversity of that library) of sequences related to the cluster 1 low-G+C-content gram-positive bacteria that were dominant in the mixolimnion and oxycline libraries (Fig. 2A). However, these sequences were not found in the monimolimnion library. Similarly, resuspension of bacterial assemblages from sediments seems unlikely because the sediments in this portion of the lake are isolated from wind-driven surface currents (41), and there is no bioturbation.
Although the strong physicochemical gradients in Mono Lake no doubt contribute to the overall diversity of the bacterioplankton in the lake by supporting a range of redox environments (niches), it is unlikely that this element of habitat diversity is responsible for the elevated diversity in the monimolimnion, as these gradients are weak there relative to other locations in the water column (the oxycline, for example). Microscope observations and transmissometer data indicated elevated particle and floc concentrations in Mono Lake bottom water, which may contribute to habitat diversity and thus elevated phylogenetic diversity. Another hypothesis to explain the elevated diversity of monimolimnion bacterial populations is that it may be possible to maintain a higher diversity of energetic pathways in anoxic environments because organisms whose metabolism is based on oxygen as an electron acceptor do not out-compete organisms that use other, less energetically efficient, forms of metabolism; such a situation would lead to the retention of higher metabolic diversity and thus ecological diversity under anaerobic conditions than under aerobic conditions. Madrid et al. (42) also found higher diversity in the deep (1,310-m) anoxic water of the Cariaco Basin than in the overlying oxic oceanic water. Regardless of the reason, these results suggest that the ecological forces that act to structure aerobic microbial communities are fundamentally different from those that act to structure anaerobic microbial communities.
The differences in the phylogenetic affinities of dominant phylotypes from mixolimnion versus monimolimnion samples are consistent with the differences in the physicochemical milieus of these two habitats and with what is known from other studies of the physiology, ecology, and distribution of representatives of these groups. Although they are not closely related to the Mono Lake high-G+C-content gram-positive bacteria, previous studies have shown that marine planktonic high-G+C-content gram-positive bacteria are generally found in the photic zone (22). The frequency of Picocystis chloroplast 16S rRNA gene sequences in our libraries is consistent with the vertical distribution of Picocystis cells in the lake. The bacterial community that we found in the monimolimnion sample is not typical of the water column but rather resembles that of sediment or soil. For example, low-G+C-content gram-positive bacteria are characteristic of deep sea sediments (40). The largest clade of sequences retrieved from sulfide-rich, anoxic sediment from a meromictic eastern Antarctica lake clustered with low-G+C-content gram-positive bacteria (5). Madrid et al. (42) also reported that the diversity and composition of 16S rRNA gene sequences from their deep anoxic water sample resembled those of a sediment or soil sample; however, low-G+C-content gram-positive bacteria were not abundant in this library, which was dominated by sequences most closely affiliated with the
-Proteobacteria.
Sequences from mixolimnion and oxycline libraries that were related to the CFB group were affiliated with the order Cytophaga-Flavobacteria, which is dominated by aerobic bacteria. CFB sequences from the monimolimnion and some of the chemocline CFB sequences fell within the order Bacteroidales, which is dominated by anaerobes that are commonly members of gut microflora.
Although the compositions of the mixolimnion and oxycline libraries were similar, two mixolimnion phylotypes were not found in the oxycline samples and four oxycline phylotypes were not found in the mixolimnion samples. One of the unique oxycline phylotypes was closely related (98% similarity) to Thiomicrospira spp., chemolithoautotrophic bacteria that use reduced sulfur compounds, such as sulfide, thiosulfate, and sulfur, as an energy source. The oxycline is a logical place in the physicochemical gradient of the lake to find sulfide-oxidizing bacteria like Thiomicrospira spp. It was shown that Thiomicrospira spp. increased at an interface between the oxidized surface and the anaerobic layer of black mud, assimilating the reduced product from the lower anoxic zone (6). Another unique oxycline phylotype was distantly related (88% similarity) to Hahella chejuensis, a facultative anaerobe capable of using several different carbohydrates as a carbon source (38). The two other unique oxycline phylotypes were related to Cytophaga spp. The planktonic alga P. salinarium was abundant in the oxycline samples, coinciding with peaks in fluorescence, chlorophyll, and DOC, suggesting that algal exudates provided the major source of organic carbon in this environment. Members of the CFB have been identified as being associated with the degradation of complex organic substrates (53, 57). Yager et al. (71) found that Cytophaga spp. dominated during the peak of an algal bloom. The CFB phylotypes that appeared to be unique to the oxycline samples may have been involved in degrading organic matter produced by Picocystis, which has an unusual composition (39).
If we accept the criterion that two organisms with 16S rRNA gene sequences that are 93% similar are related at the level of genus (48), 75% of our cloned sequences, representing 50 phylotypes, belonged to taxonomic groups that are distinct at the genus level or higher and that have not been previously described. None of the sequences that we obtained and that were related to CFB, low-G+C-content gram-positive bacteria, and candidate division groups were similar (>97% sequence similarity) to 16S rRNA gene sequences of cultivated organisms or to environmental sequences recovered from other sites. The sequences found in this study are also different from sequences previously retrieved from other alkaline environments. Most previous studies of alkaline environments have been based on classical isolation techniques (15),which may not completely sample the diversity present in a habitat (32). It is estimated that >99% of microorganisms observable in nature typically are not cultivated by standard techniques (2, 18). Thus, the data presented in this report significantly advance the understanding of the microbial ecology of alkaline, hypersaline environments.
We are grateful to G. Lecleir for help with sample collection and for providing cell count and DOC data. S. B. Joye and S. Carini supplied sulfide and ammonia data, and A. Burd helped with diversity analyses. E. Biers helped with suggestions and excellent cakes.
|
|
|---|
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»