Department of Marine Sciences, University of Georgia, Athens, Georgia 30602-3636
Received 20 October 2003/ Accepted 20 February 2004
| ABSTRACT |
|---|
|
|
|---|
| INTRODUCTION |
|---|
|
|
|---|
85-g liter1 salinity and pH 9.8, Mono Lake is a midpoint on the gradient of saline lakes, and as such, its chemistry and the general features of its ecology are similar to those of numerous other saline lakes and similar environments in the western United States and elsewhere in the world (21). Summer thermal stratification and the high rates of primary production characteristic of many temperate saline lakes result in seasonal bottom water anoxia. Saline lakes are also prone to longer periods of anoxia during meromixis (persistent chemical stratification). During the past two decades, two extended periods (>5 year) of meromixis were initiated at Mono Lake by unusually high runoff. Reduced species, especially sulfide, ammonia, and methane, accumulate to high concentrations beneath the chemocline (14, 17). Because of the strong redox gradients that are generated when the lake is stratified, chemolithoautotrophy may be an important trophic mode in Mono Lake, particularly during meromixis.
A common carbon fixation pathway for chemolithoautotrophic microorganisms is via the reductive pentose phosphate or Calvin-Benson-Bassham (CBB) cycle. One of the key enzymes in the CBB pathway, ribulose-1,5-bisphosphate carboxylase/oxygenase (RuBisCO) (EC 4.1.1.39), catalyzes two competing reactions: the carboxylation or the oxygenolysis of ribulose-1,5-bisphosphate. RuBisCO exists in two forms. The form I enzyme is made up of eight large catalytic subunits and eight small subunits, with the cbbL (rbcL in plants) gene encoding the large subunit and the cbbS (rbcS in plants) gene encoding the small subunit (16). This structural form is typically found in all terrestrial plants, algae, cyanobacteria, and most photo- and chemoautotrophic bacteria (26).
Form II RuBisCO is composed of two large catalytic subunits that are encoded by the cbbM gene and arranged in a homodimeric structure. Although the active site residues are conserved in the two forms, the overall level of similarity between the forms I and II large subunits is only 25 to 30% (24). The cbbM gene has been found in several photosynthetic bacteria (27), aerobic and facultatively anaerobic chemoautotrophic bacteria (4, 11, 16, 22, 25), and dinoflagellates (18).
Due to its relatively large size, degree of conservation, widespread distribution, functional significance, and an increasing number of published sequences from chemoautotrophs, the RuBisCO gene is particularly useful for studying chemoautotrophic populations, in addition to and for comparison with other genes, such as the rRNA gene (28). We investigated the use of RuBisCO as an indicator of the diversity of carbon-fixing microorganisms in the water column of Mono Lake. In this study, we analyzed the distribution and phylogenetic richness of the cbbL and cbbM genes in four samples selected to represent the range of water column redox conditions in Mono Lake.
| MATERIALS AND METHODS |
|---|
|
|
|---|
PAR sensor, a WetLabs 10-cm-path-length transmissometer, and a WetLabs WetStar fluorometer. Oxygen profiles were taken with a YSI polarographic oxygen sensor equipped with a Clark-type electrode. Water samples were collected at discrete depths using a Niskin water sampler. Samples were screened as necessary to remove brine shrimp, placed in dark plastic bottles, and held in coolers until they were processed (within 4 h of collection).
DNA collection and extraction.
Methods of DNA collection, extraction, and purification were as described previously (9, 13). Water samples (500 to 1,000 ml) were pressure filtered (
50 kPa) through Millipore Sterivex filter cartridges (0.22-µm pore size) to collect microbial biomass for subsequent DNA extraction. Excess water was expelled, and then the cartridges were filled with extraction buffer, capped, frozen on dry ice, shipped to the laboratory, and stored at 70°C until processed. Blanks were prepared with each set of samples using Sterivex cartridges through which no water had been filtered.
PCR and cloning.
Samples from 2 m (mixolimnion), 17.5 m (base of the oxycline, also corresponding to a fluorescence maximum/transparency minimum), 23 m (top of the chemocline), and 35 m (monimolimnion) were chosen for analysis. These depths were chosen because they represent the range of water column redox conditions (see Fig. 1). Although only four samples taken from one station on one day were analyzed in detail, extensive analysis of the distribution of Mono Lake bacteria ribotypes using survey techniques (PCR and denaturing gradient gel electrophoresis) indicates that the composition of these samples is representative of other depths, stations and times (12; J. T. Hollibaugh, N. Bano, S. Humayoun, and G. LeCleir, unpublished data). These are also the same samples that were used in an analysis of the diversity of 16S rRNA genes of Mono Lake bacteria (13).
|
Phylogenetic analysis.
All sequences were obtained from an automatic sequencer operated by Molecular Genetics Instrumentation Facility. A total of 57 clones (11 from 2 m, 17 from 17.5 m, 15 from 23 m, and 14 from 35 m) were selected randomly; DNA was extracted from them, purified, and sequenced using the plasmid primer SP6, T7, or both. Sequence reads were
650 bp. Most clones were sequenced in only one direction; however, a subset of clones representing major clades were sequenced from both ends to obtain complete sequences of the inserts. Nucleotide sequences were aligned and converted to inferred amino acid sequences by using programs in the Wisconsin package (version 10.0) Genetics Computer Group Inc. Sequences (nucleotide and amino acid) were compared to known sequences using BLAST (2). Phylogenetic trees of RuBisCO genes were constructed using the ARB program. The topology of the nucleotide tree was established using full-length sequences (772 bp after primer sequences were removed) and the neighbor-joining method. Bootstrap analysis (5,000 replicates) was performed, and then partial sequences (including some of the reference sequences) were inserted into the consensus tree. The inferred amino acid tree used the neighbor-joining method and evolutionary (Kimura) distances and was constructed with PHYLIP (8), using 228 positions. Bootstrap values for the amino acid tree are based on 100 replicates. The final trees were unrooted, with a Chlorella-like algal form I RuBisCO large subunit sequence used as the outgroup.
Nucleotide sequence accession numbers.
Sequences have been deposited in GenBank under accession numbers AY291478 to AY291533 and AY293400 as indicated on the figures.
| RESULTS AND DISCUSSION |
|---|
|
|
|---|
Phylogenetic affiliations of RuBisCO genes.
Amplicons of the form I RuBisCO gene (cbbL) were obtained from all samples. The diversity of cbbL genes we retrieved from our Mono Lake samples was lower than we expected. Coverage estimates (10) were 82, 83, 93, and 93% for libraries from 2, 17.5, 23, and 35 m, respectively. Fifty-two of the fifty-seven sequences we obtained are classified as form IA RuBisCO, and five form 1B sequences were from cyanobacteria. Nucleotide (Fig. 2) and inferred amino acid (Fig. 3) tree topologies were generally congruent, with one exception (discussed below). The nucleotide sequences fell into nine distinct clades, with the largest clade containing 42% of the sequences. One clade of four closely related sequences merged with the largest clade when inferred amino acid sequences were used for phylogenetic analysis, leaving only eight distinct clades, with the largest containing 52% of the sequences. Most of the sequences we obtained are not closely affiliated (nucleotide similarities, <90%) with RuBisCO sequences in the databases; however, some of them were closely related to sequences obtained from Mono Lake isolates or enrichments.
|
|
Clade A contains three sequences from the oxic mixolimnion and one from the monimolimnion that were most similar to RuBisCO sequences from T. denitrificans, A. vinosum, and M. capsulatus.
Clade B is dominated by sequences retrieved from chemocline and monimolimnion samples (Fig. 2). Sequences in Clade B were most similar to RuBisCO sequences from Nitrosomonas sp., Nitrosospira sp., N. winogradskyi, and T. intermedia (T. intermedius).
Clade C sequences were found in chemocline and monimolimnion samples and were most similar to those of Mono Lake isolate strain MLHE-1 (20) and to clade B (nucleotide sequence, Fig. 2). We recently retrieved RuBisCO sequences with >99% similarity (431 of 435 bp) to those of clade C from an arsenate-reducing, sulfide-oxidizing enrichment culture derived from a Mono Lake sample.
Clade D contains form 1B sequences that are >99% similar to the RuBisCO sequence from a Cyanobium-like unicellular cyanobacterium isolated from Mono Lake (C. Budinoff and J. T. Hollibaugh, unpublished data). The Cyanobium is abundant in subchemocline waters (
108 cells liter1; Hollibaugh, unpublished data), and although it grows photoautotrophically in the laboratory (C. Budinoff, unpublished data), it has not yet been observed in surface water samples from Mono Lake. It is unclear whether the monimolimnion population is actively growing, represents a relict population from the last time the lake completely mixed, or is maintained by a sinking flux (for example, brine shrimp fecal pellets).
Branches E and H each contain one sequence retrieved from the mixolimnion sample. These sequences were most closely related to cyanobacterial genes (nucleotide similarities of 77 and 82%, respectively).
Clades F and G grouped together on the second major branch of the Mono Lake RuBisCO tree along with the reference sequence from Hydrogenovibrio marinus (nucleotide similarities of 89 and 87%, respectively). Clades F1 and F2 account for more than half of the sequences we retrieved. Their nucleotide sequences are >95% similar; however, their amino acid sequences are >99% similar, because differences in the nucleotide sequences correspond to variability in the third base pair of codons for the same amino acids. No Clade F sequences were retrieved from the surface water sample; most were retrieved from the micro-oxic or anoxic region between the oxycline and the chemocline (17.5- and 23-m samples), where gradients of reduced substrates were strongest (Fig. 1). Sequences that were identical to those in clade F1 were retrieved from dilution cultures of Mono Lake water amended with nitrate and sulfide (Hollibaugh et al., unpublished). Clade G contains sequences from the mixolimnion and the oxycline.
No PCR products were obtained from any samples when primers specific for the form II RuBisCO gene, cbbM, were used for amplification, despite successful amplification of the gene from the positive control (R. rubrum, ATCC 9791). We also failed to obtain RuBisCO sequences when samples were amplified with RuBisCO primers used in a study of chemoautotrophs in groundwater (1).
Given the variety and strength of geochemical gradients in the Mono Lake water column, we expected a greater diversity of chemoautotrophs and thus of form 1a and form 2 RuBisCO genes in our samples. The low diversity of RuBisCO sequences may reflect a low diversity of organisms; low diversity of RuBisCO genes due to physicochemical conditions in Mono Lake; inefficient extraction of DNA from some groups of organisms, for example, methanotrophs; or restricted specificity of the primer set.
Elsaied and Naganuma (7) report identifying a total of 29 different cbbL sequences (operational taxonomic units [OTUs])from 300 cloned genes. Their samples covered a range of habitats associated with a hydrothermal vent site, including sediment, water, and symbionts. The richness of cbbL genes they detected (0.10 OTU clone1) is comparable to the richness we observed (0.12 OTU clone1), despite the differences in habitat diversity. However, they detected cbbM genes in their samples (24 OTUs in 250 cbbM clones) where we detected no cbbM genes in Mono Lake samples. Alfreider et al. (1) recovered 24 cbbL OTUs in 100 clones (0.24 OTU clone1) and 6 cbbM OTUs in 40 clones (0.15 OTU clone1) in groundwater samples. Nanba et al. (19) retrieved at least 7 different OTUs from 10 clones of a microbial mat sample and at least 23 different OTUs from 110 clones retrieved from volcanic soil samples. These comparisons are not entirely satisfactory, because each study used a different primer set. Furthermore, the heterogeneity of these environments is greater than that of the Mono Lake water column, but there are as yet no other comparable studies from other aquatic environments.
RuBisCO gene diversity might be low because Mono Lake presents a very stable environment with regard to RuBisCO substrates. Because of the elevated pH (9.8) and high dissolved inorganic carbon concentrations of Mono Lake (360 mM; W.-J. Cai and Y. Wang, unpublished data), Mono Lake partial CO2 pressure and alkalinity are elevated (512 and 593 µatm and 661 and 634 mM, respectively, measured in two samples by the procedure described in reference 3 [Cai and Wang, unpublished data]) relative to those of offshore surface seawater in equilibrium with the atmosphere (
365 µatm and 2.2 mM). Oxygen concentrations are low and constant (
0) anywhere in the lake below the oxycline. Thus, the Mono Lake environment may be less variable with regard to the conditions under which RuBisCO functions than other aquatic environments, which may be reflected in lower diversity of RuBisCO genes. Alternatively, chemoautotrophic populations may be diverse, reflecting the diversity of physiochemical conditions in the lake's highly stratified water column, but Mono Lake chemoautotrophs may rely on other pathways for carbon fixation or the RuBisCO primer set we used may not adequately sample this diversity.
Another possible reason for the apparent low richness of RuBisCO genes in Mono Lake samples may be that the primers we used recognize a subset of cbbL sequences that are not well represented in Mono Lake chemoautotrophs. The composition of the Mono Lake chemoautotroph population may be significantly different from those of populations of chemoautotrophs in the other environments that have been studied. This speculation is supported by the observation that none of the sequences we obtained was closely related (nucleotide similarities of <90%) to environmental sequences obtained by Alfreider et al. (1) or Elsaied and Naganuma (7).
| ACKNOWLEDGMENTS |
|---|
We are grateful to G. LeCleir for help with sample collection and for providing cell count and DOC data and to C. R. Budinoff for his help with data analysis. S. B. Joye and S. Carini supplied sulfide, ammonia, and methane data. Field laboratory facilities, limnological data processing, and other survey and logistical support were provided by R. Jellison, S. Roll, and the Sierra Nevada Aquatic Research Laboratory with support from the National Science Foundation (MCB-9977901). We are also grateful to two anonymous reviewers for comments that greatly helped improve the manuscript.
| FOOTNOTES |
|---|
| REFERENCES |
|---|
|
|
|---|
This article has been cited by other articles:
| |||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||
| J. Bacteriol. | Microbiol. Mol. Biol. Rev. | Eukaryot. Cell | All ASM Journals |
|---|