Previous Article | Next Article ![]()
Applied and Environmental Microbiology, July 2002, p. 3606-3613, Vol. 68, No. 7
0099-2240/02/$04.00+0 DOI: 10.1128/AEM.68.7.3606-3613.2002
Copyright © 2002, American Society for Microbiology. All Rights Reserved.
GBFGerman Research Center for Biotechnology, Department of Environmental Microbiology, D-38124 Braunschweig, Germany
Received 10 December 2001/ Accepted 4 April 2002
|
|
|---|
|
|
|---|
The depth dependence of the structure of phyto- and zooplankton communities in lakes and oceans has been studied since the beginning of aquatic ecology and has contributed enormously to the understanding of the functioning of freshwater and marine ecosystems. Despite the recognition of the important role of bacteria in the microbial food web and the contribution of bacteria to the functioning of aquatic ecosystems, knowledge about the depth distribution of these organisms is still scarce (1, 5). Pioneering work has been done on the depth-dependent community structure of bacterioplankton in the oceans by DeLong (6), Giovannoni et al. (15), and Fuhrman et al. (13). These studies were based on analysis of 16S ribosomal DNA (rDNA) clone libraries from deep ocean samples (500 to 2,000 m) and provided insight into the depth distribution of single taxa. Field et al. (12) obtained higher resolution of the variation in bacterial communities with depth by using group-specific 16S rRNA-targeted oligonucleotide probes and environmental bulk DNA-RNA for assessment of the abundance of specific taxa belonging to the
subgroup of the Proteobacteria in the Sargasso Sea. However, community structure can be assessed most rapidly by community fingerprint methods that are based on electrophoretic separation of 16S rDNA amplicons or 5S rRNA obtained from natural microbial communities (27, 30, 34, 35). The first community fingerprinting analysis of bacterioplankton with depth was done by using amplicons of the V3 region of 16S rRNA and archaeal and bacterial primers (29). Øvreås et al. (29) found that the diversity of the bacterial community decreased from the surface to the bottom in an estuary with a strong vertical gradient of salinity and oxygen. However, a major problem of PCR-based methods is that they can be biased for some phylogenetic groups, and quantitative assessment of community diversity is difficult (39). Therefore, we have used 5S rRNA fingerprinting to obtain a direct molecular fingerprint from natural communities and have calculated, based on the relative densities of individual 5S rRNA bands, the diversity of the bacterioplankton during the year (20).
The goals of the present study were to compare changes in physical, chemical, and biological parameters with depth at a high spatial resolution (i.e., every meter) with changes in the structure of the bacterioplankton community and to extend the previous seasonal study in terms of depth (20). The ecosystem chosen was Lake Plußsee, a naturally eutrophic lake in northern Germany which has been examined extensively (28). We analyzed the 5S rRNA fraction of environmental rRNA obtained directly (i.e., without PCR) from the bacterioplankton. The structure of the bacterioplankton community as a function of depth was analyzed based on high-resolution 5S rRNA community fingerprints and partial or full sequences of individual 5S rRNA.
|
|
|---|
Biological background parameters, such as total bacterial counts, total heterotrophic nanoflagellate (HNF) counts, and total particulate chlorophyll a concentrations, were measured by standard techniques; details are given elsewhere (4).
5S rRNA fingerprints of bacterioplankton.
5S rRNA analysis of bacterioplankton started with extraction of the total RNA directly from the frozen filter sandwiches. This extraction was based on mechanical extraction with glass beads in a bead beater in combination with chemical extraction with phenol and sodium dodecyl sulfate (SDS) (8, 16). In short, total environmental RNA was extracted from the sandwiches by using 4 g of 2- and 3-mm-diameter glass beads, a high-speed cell disrupter (Microdismembrator; Braun-Diessel Corp.) with a Teflon extraction cell, and a mixture of 5 ml of RNA extraction buffer and 5 ml of buffer-saturated phenol (17). Extraction buffer consisted of 50 mM sodium acetate, 10 mM EDTA, and 1% SDS (pH 4.2). Phenol was saturated with extraction buffer lacking 1% SDS until a final pH of 4.2 to 4.4 was reached. A thin slurry (total volume, about 12 ml) of the two filters and the extraction mixture was generated by vibration for 2 min with the cell disrupter. This slurry was then centrifuged for 10 min at 8,500 x g. After repeated phenol (pH 4.2) extraction of the pellet, the supernatants were pooled and treated twice with 5 ml of chloroform. Total RNA was precipitated from the pooled supernatant by addition of ethanol (2.5 times the volume of the supernatant) and storage at -20°C overnight.
The total RNA obtained by the filter extraction procedure was 3' end labeled with cytidine 3',5'-[5'-32P]bisphosphate (specific activity, 3,000 Ci/mmol) by using RNA ligase with an efficiency of about 62% (7, 10). The radioactively labeled RNA was subjected to denaturing high-resolution electrophoresis on a 10% polyacrylamide gel (Sequicel; 0.4 mm by 38 cm by 80 cm; Bio-Rad) by using running conditions consisting of a stepwise increase over a 5-h period from 100 to 300 W. After electrophoresis, the gel was exposed to X-ray film or a storage phosphor screen (Molecular Probes Corp.) for 30 min to several hours. More details concerning RNA extraction and electrophoretic analysis are given elsewhere (8, 16, 17, 19, 40).
Comparative analysis of 5S rRNA fingerprints.
The 5S rRNA fraction was evaluated quantitatively by scanning the autoradiograms with an optical gel scanner (Hirschmann Corp.) and by using a phosphorimager (Molecular Probes Corp.). In this way the abundance of single 5S rRNA bands was quantified in terms of baseline corrected peak areas and in terms of number of pixels per band, respectively. These arbitrary units were used to calculate the relative abundance of every individual 5S rRNA band and its contribution to the total 5S rRNA.
The diversity of the community was determined by determining the number of 5S rRNA bands per sample and by calculating the Shannon diversity index (H) with the following equation:
![]() |
Sequence analysis of single environmental 5S rRNA.
After high-resolution electrophoresis, 32P-labeled 5S rRNA bands were excised from the gel, eluted with 150 µl of elution buffer (0.5 M sodium acetate, 0.1% SDS, 1 mM magnesium chloride, 100 µM EDTA; pH 7.8), and ethanol precipitated at -20°C overnight. Before application to the sequencing gel, 5S rRNA was enzymatically digested with base-specific RNases (9, 16). The cleaved 5S rRNA of a single excised band was run on a high-resolution sequencing gel (80 by 40 cm; thickness, 0.4 mm;10% acrylamide LongRanger [FMC Corp.] prepared according to the protocol of the manufacturer) to determine the positions of the bases in comparison to the bases of totally hydrolyzed 5S rRNA of Escherichia coli, which was used as the molecular size marker. This sequencing gel provided information about the positions of almost all bases of the 5S rRNA from the 5' end. Only around 2% of the molecule (i.e., two or three nucleotides at the 3' end) could not be read. These sequences could lead to classification or identification of the taxon represented by a single band provided that the reference sequence of a known bacterial species is available in the 5S rRNA sequence database (33). Also, this sequence information was used to determine the homogeneity of a band (i.e., if one or more taxa formed the band). Furthermore, it allowed comparison of 5S rRNA bands of the same length for samples from different depths and times of the year. Such comparisons were also possible with partial sequences (e.g., by using only RNase T1 to examine the positions of all guanine residues within a 5S rRNA band of a specific size).
|
|
|---|
![]() View larger version (19K): [in a new window] |
FIG. 1. Physical, chemical, and biological background data for depth profiles of Lake Plußsee during the spring (25 April 1994). (a) Temperature ( ), oxygen content ( ), and phosphorus content ( ); (b) chlorophyll a content ( ), total bacterial number ( ), and number of HNF ( ).
|
![]() View larger version (20K): [in a new window] |
FIG. 2. Physical, chemical, and biological background data for depth profiles of Lake Plußsee during the fall (19 September 1994). (a) Temperature ( ), oxygen content ( ), and phosphorus content ( ); (b) chlorophyll a content ( ), total bacterial number ( ), and number of HNF ( ).
|
![]() View larger version (66K): [in a new window] |
FIG. 3. Autoradiograms of 5S rRNA from bacterioplankton obtained at different depths. Samples of the water column were taken at different depths. For detection the 5S rRNA was 3' end labeled with 32P. (a) Depth profile for the spring (25 April 1994). Lanes M contained a molecular weight standard consisting of hydrolyzed 5S rRNA from E. coli. The intensity of the bands in the boxed area was electronically reduced for better visibility of the single bands. Each lane is labeled with the depth (in meters) at which the sample examined was obtained. (b) Depth profile for the autumn (19 September 1994): Lanes M contained a molecular weight standard consisting of hydrolyzed 5S rRNA from E. coli. Each lane is labeled with the depth (in meters) at which the sample examined was obtained.
|
The results of a detailed quantitative analysis of the amount of 5S rRNA per abundant band in the autoradiograms, obtained by using a phosphorimager for quantification, are shown in Fig. 4. In both depth profiles the number of abundant bands was rather limited, ranging from three to seven. Also, these abundant bands (i.e., bands that accounted for more than 4% of the total 5S rRNA) comprised more than 90% of the total 5S rRNA in most samples (Fig. 4B). The spring profile was slightly more variable in this respect than the fall profile due to the larger number of less abundant bands occurring in the metalimnion. Another general trend for both depth profiles was that in the epilimnion no single band exceeded a relative abundance of 34% and a small number of bands (three to five bands), each comprising 10 to 30% of the total, represented the majority of the bacterioplankton community. With increasing depth one band became dominant in both depth profiles; i.e., it had a relative abundance of more than 50%. Furthermore, in both profiles for depths below 15 m the dominant band was the same size (121 nt) despite the anoxic conditions in the fall. The two depth profiles did differ, however; when anoxic conditions were reached at depths below 6 m, the 123-nt band disappeared and a 113-nt band appeared. In general, changes in the community composition were more gradual in the spring profile than in the fall profile.
![]() View larger version (23K): [in a new window] |
FIG. 4. (A) Relative amounts of the single abundant 5S rRNA bands for the spring shown in Fig. 3a. Different colors indicate 5S rRNA bands of different sizes that account for more than 4% of the total 5S rRNA calculated as detailed in Materials and Methods. (B) Relative amounts of single abundant 5S rRNA bands for the fall shown in Fig. 3b.
|
![]() View larger version (24K): [in a new window] |
FIG. 5. (a) Total number of bands ( ) (bands containing more than 0.1% of the total 5S rRNA) and Shannon diversity index ( ) calculated from all 5S rRNA bands at each depth during the spring. The chlorophyll a content ( ) is shown for comparison. (b) Total number of bands ( ) and Shannon diversity index ( ) calculated from all 5S rRNA bands at each depth during the fall. The chlorophyll a content ( ) is shown for comparison.
|
In the fall, 5S rRNA bands at 116, 120, 121, 122, and 123 nt were compared for the upper 5 m. In addition, the 121-nt band was sequenced from all samples of the depth profile. There were no sequence differences detected in any of the fall 5S rRNA bands with depth; i.e., all bands sequenced in the upper 5 m and the 121-nt band found throughout the water column were the same. A comparison of the sequences from the spring with the sequences from the fall showed significant differences for all of the single size classes analyzed except the 116- and 122-nt bands. These bands had identical partial sequences; i.e., all As and Gs were at the same position in the 5S rRNA for both seasons. The 121-nt band, which was the most abundant 5S rRNA in the hypolimnion in both depth profiles, had different sequences in the two seasons, despite the identical sizes. These sequence data indicate that different taxa dominated the hypolimnion in the spring and in the fall.
The only 5S rRNA for which a definite phylogenetic relative could be found was the 116-nt band, which had a sequence identical to a Comamonas acidovorans sequence (accession number AJ131594). The presence of C. acidovorans has been documented in samples from Lake Plußsee before (11, 20). For all other bands except the 123-nt band either the sequence information obtained was not complete enough for alignment with sequences from the international database or no close relative could be found; the closest phylogenetic neighbor of the 123-nt band was in the genera Thiovolum and Campylobacter of the
subgroup of the Proteobacteria.
|
|
|---|
Vertical gradients of physical and chemical parameters strongly affect the biology and ecology of all pelagic ecosystems (22). This universal finding is also reflected in the structure of the bacterioplankton community. The most striking feature of both depth profiles was the much lower diversity in the hypolimnion than in the upper layers. The decrease can be seen in the reduction of the Shannon index by about one-half, as well as in the number of bands, which decreased from about 14 to 7 in both depth profiles (Fig. 5). One exception to this overall trend was the 25-m sample collected in the spring, for which the number of bands increased. Conceivably, the diatoms which dominate the phytoplankton of Lake Plußsee in early spring were still being degraded in the deep hypolimnion after they sank in late March (28, 32).
When the diversity data for the depth profiles were compared with the data from the seasonal study of bacterioplankton in Lake Plußsee, it was obvious that the low diversity, as indicated by a Shannon index of around 0.4, in the hypolimnion was encountered only during the three minima occurring in the epilimnion at special times of the year (20). During most of the year the Shannon index remained in the range found for the epilimnion for both depth profiles. The first minimum for the diversity index in the seasonal study could be explained by strong grazing by cladocerans, whereas grazing by HNF did not have a significant effect on diversity. Grazing by HNF did not directly affect the high diversity in the upper layers of depth profiles, where HNF grazing is the major mortality factor in the epilimnion (41). On the other hand, viral lysis was demonstrated to be the key mortality factor in the anoxic hypolimnion of Lake Plußsee, where all grazing was strongly reduced (41). This result corroborates the finding that only very few members of the bacterioplankton community were active in the hypolimnion because low-diversity bacterial communities can be controlled more easily by viral lysis than high-diversity communities (37). The major reason for the low diversity in the hypolimnion for both depth profiles might be the decrease in the dissolved organic carbon content, which not only decreased by about 50% with depth but also accounted for a far higher proportion of refractory organic compounds in the hypolimnion than in the epilimnion of Lake Plußsee (26).
Depth-dependent and seasonal changes in detailed structure of the bacterioplankton community.
The presence and depth distribution of C. acidovorans in Lake Plußsee have been demonstrated previously by immunofluorescence microscopy (11, 42). Additionally, a seasonal study in which 5S rRNA profiling was used showed that C. acidovorans is abundant at a depth of 1 m during all seasons (20). Both depth profiles obtained by using 5S rRNA analysis confirmed the presence of C. acidovorans in the whole profile, with an average abundance of about 6.0%. (The levels in deeper samples were around 3% and therefore are not shown in Fig. 4). This finding is in contrast to the results obtained by immunofluorescence microscopy that indicated that C. acidovorans cells disappeared at depths below 10 m in the spring of 1993 and that there was a significant lack of C. acidovorans cells in the metalimnion during the fall of 1996 (11, 42). The most likely explanation for this finding is that the monoclonal antibodies used for detection of C. acidovorans cells are specific for a certain subspecies type, whereas the 5S rRNA bands represent the whole taxon. Therefore, C. acidovorans was present in abundance in Lake Plußsee not only during all seasons but also at all depths, including the anoxic hypolimnion. Furthermore, the presence of significant amounts of rRNA indicated that this group of bacteria was metabolically active at all depths year round. This supports the hypothesis, put forward earlier (42), that C. acidovorans is a K-strategist well adapted to the pelagic environment of a eutrophic lake because of its relative resistance to protozoan grazing and its ability to utilize refractory dissolved organic matter.
Partial sequencing revealed differences in all bands except the 116- and 122-nt bands in the different seasons. The strong seasonal differences make sense considering the fact that the phytoplankton communities during the two seasons were completely different and that in the fall the hypolimnion was anoxic. Detection of members of the
subgroup of the Proteobacteria in the anaerobic water column is unusual but has been observed with clone libraries from coastal marine water samples obtained during an algal bloom and from the chemocline of a meromictic lake in the Swiss Alps (2, 21).
Comparison of bacterioplankton community structure with other depth distributions.
A variety of studies have been performed with different molecular methods to understand the depth distribution of the bacterioplankton community (14, 25). Most of these studies were based on clone libraries of 16S rDNA genes (see reference 14 for a review), followed by hybridization studies and community fingerprinting. While clone library studies provide insight into community composition and diversity, fingerprinting studies reveal the overall structure of the whole community and enable direct comparisons of a variety of samples. Fingerprinting studies of the bacterioplankton of a brackish meromictic lake and the central Baltic Sea indicated that there was higher diversity in the upper layer than in the anoxic bottom layer, as found for Lake Plußsee (18, 29). Recent studies of bacterioplankton in a deep ultraoligotrophic lake and the eastern Mediterranean Sea showed that there were distinct changes with depth but did not reveal substantial decreases in diversity (24, 38). Several studies performed with group-specific rRNA-targeted probes revealed the depth dependence of these groups in the water column of the oceans (12, 23). The amounts of some of these groups correlated well with microbial processes, such as the intensity of primary production (36). In general, the following picture appears to emerge for the structure of the bacterioplankton community with depth. Pelagic environments have specific bacterioplankton community structures at specific depths, and the overall diversity and composition might vary from layer to layer depending on the trophic status and the availability of electron acceptors.
Conclusions.
We demonstrated that there were changes in bacterioplankton community structure with depth in a eutrophic lake. These changes corresponded well with the three layers, epilimnion, metalimnion, and hypolimnion, of the lake. The more the lake was stratified, the better the three different communities were separated. Whereas the community composition generally varied with the season, some members of the community, such as C. acidovorans, were abundant at all depths during both seasons. Another common feature was the reduction in the overall diversity of the bacterioplankton with depth. This decrease in diversity means that the diversity of the bacterioplankton community in the hypolimnion was quite low compared to the diversity of the community in the epilimnion, where the diversity is low only two or three times during the year (20).
This project was supported by funds from Bundesministerium für Bildung, Wissenschaft, Forschung und Technologie grant BEO-0319433B and Deutsche Forschungsgemeinschaft grant Ho 930/2-1.
|
|
|---|
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»