Previous Article | Next Article ![]()
Applied and Environmental Microbiology, March 2008, p. 1856-1868, Vol. 74, No. 6
0099-2240/08/$08.00+0 doi:10.1128/AEM.01719-07
Copyright © 2008, American Society for Microbiology. All Rights Reserved.
,
Ellen B. Duffy,1,
Vincenzo Milano,1,
Maggie McLellan,1,¶
Brett M. Stevens,1,||
Charles W. Boylen,1,3 and
Sandra A. Nierzwicki-Bauer1,3*
Darrin Fresh Water Institute, Rensselaer Polytechnic Institute, 5060 Lake Shore Drive, Bolton Landing, New York 12814,1 Skidaway Institute of Oceanography, 10 Ocean Science Circle, Savannah, Georgia 31411,2 Department of Biology, Rensselaer Polytechnic Institute, 110 8th Street, Troy, New York 121803
Received 25 July 2007/ Accepted 23 January 2008
|
|
|---|
|
|
|---|
Before culture-independent molecular genetics-based techniques were available, freshwater microbial communities were believed to be more similar to marine and soil communities than different from marine and soil communities (41). However, more recent examinations of bacterial communities in freshwater environments have suggested that they are distinguishable from marine communities, largely by the dominance of members of the Betaproteobacteria class and representatives of the Actinobacteria class (16, 33, 58, 59, 64). Representation of these groups appears to be intermediate in estuarine environments (6, 7, 40). Investigations of bacterial communities in Adirondack lakes in the state of New York were among the first to recognize unique freshwater bacterial lineages of freshwater Betaproteobacteria and Actinobacteria (22, 33). Additionally, recent studies are revealing diversity of freshwater picoeukaryote communities that rival the diversity seen in marine environments (42). Less is known about freshwater Archaea communities, although several studies have investigated these communities (8, 26, 39, 53).
The importance and numerical abundance of these "typical" freshwater bacterial groups have been confirmed in studies employing fluorescence in situ hybridization approaches (17). Interestingly, the abundance and distribution of Actinobacteria and Betaproteobacteria suggest that these two groups are influenced differently by lake hydrologic type, nutrient conditions, seasonality, and grazing pressures (9, 29, 61). Recently, Hahn and others have succeeded in isolating representatives of some of the freshwater Actinobacteria and Betaproteobacteria species and have begun to demonstrate the ecological significance of some phylogenetically distinct groups of these freshwater bacteria (19, 20, 38, 47, 59).
One of the longest and most comprehensive investigations of freshwater systems has focused on lakes and ponds in the Adirondack Park located in northeastern section of New York State in the United States. The Adirondack Park is a 6-million-acre forested upland and mountainous region, which contains approximately 2,770 lakes and ponds in six different watersheds (11, 25). The Adirondack region receives an average of 100 to 150 cm of precipitation annually, and due to the low acid-neutralizing capacity of most Adirondack waters, the region has been severely impacted by acid precipitation (25). In recent years, Adirondack lakes have been on a trajectory to recovery from acidification in response to decreasing levels of air pollution (11, 25). Lakes in the Adirondacks are characterized hydrologically and physiochemically into distinct types defined by their geomorphology, drainage, and chloride, calcium, and dissolved organic carbon (DOC) concentrations (25, 37). Most Adirondack lakes are small in terms of surface area and volume, and they have moderate drainage areas, resulting in short hydraulic retention times. Detailed hydrographic and chemical characteristics of many Adirondack lakes can be obtained from the Adirondack Lake Survey Corporation at the website www.adirondacklakessurvey.org. From 1994 to 2007, under the auspices of the Adirondack Effects Assessment Program (AEAP) supported by the U.S. Environmental Protection Agency, the potential recovery of Adirondack lakes from acidification has been investigated. The AEAP has focused on the long-term monitoring of water chemistry and biota, including microorganisms, in 30 representative lakes. In this study we report on bacterial communities in 17 of these lakes, as well as an 18th Adirondack lake, Lake George, that has not been significantly impacted by acid deposition.
Lake George is a 62.5-km-long glacial lake located in the southeastern corner of the Adirondack Park. It is relatively large with respect to other lakes in the region, with a surface area of 114 km2. On the basis of nutrient concentrations, Lake George is classified as an oligo- to mesotrophic lake with moderate alkalinity (35). Continuous long-term monitoring of water chemistry in Lake George has been ongoing since 1979 (35) (http://www.rpi.edu/dept/DFWI/).
The deleterious impact of atmospheric acid deposition on the biota of freshwater environments, particularly macroorganisms (plants and animals), has been shown to be dramatic, particularly with respect to large invertebrate and vertebrate predators (21, 44). Most acutely in the northeastern portion of North America, acid deposition in freshwater environments has resulted in the loss of aquatic species, including algae, fish, submerged vegetation, amphibians, and aquatic invertebrates (44). In recent years, primarily due to reduced atmospheric sulfur deposition, lakes, rivers, and streams across Europe and North America have experienced increasing pH with concomitant decreases in SO4 and other chemical indicators of acidic precipitation (49, 51, 57). In some, but not all, cases, the biotic communities (submerged vegetation, invertebrates, and fish) are also recovering, although generally there has been a substantial lag period between chemical and biological recovery (11, 24, 48). Because relatively little is known about the response, potential recovery, or changes of bacterial communities in association with acid deposition processes, we investigated the composition and diversity of bacterial communities in Adirondack lakes recovering from acidification (along with one lake that has not been impacted). We also examined the relationship between bacterioplankton community structure and lake physiochemical and hydrological properties to elucidate some of the factors that may influence microbial diversity in an effort to develop specific testable hypotheses concerning environmental controls on the composition of freshwater bacterial communities.
|
|
|---|
|
View this table: [in a new window] |
TABLE 1. Locations and physical characteristics of the 18 study lakes
|
Sample collection.
With the exception of Lake George, water samples from lakes that exhibited thermal stratification were collected from the epilimnion (well-mixed surface layer) and the hypolimnion (1 meter above the lake bottom). Not all sample lakes contained a hypolimnion layer during sampling (see Table 3). Of the 17 lakes that were sampled, 8 had both layers, 9 did not have a hypolimnetic layer developed, and 1 did not have an epilimnion. Depth-integrated epilimnetic water samples were collected by lowering a 2.54-cm (1-in.)-diameter hose to 1 m above the thermocline, plugging the top with a rubber cork to create a vacuum, and retrieving the water into a well-rinsed bucket. Water samples from the hypolimnion were collected using a Wildco Alpha VanDorn style water bottle. Two liters of water from each sample was transferred to sterile darkened containers, stored on ice, and transported to the lab. All samples were processed within 12 h of collection. In Lake George, three epilimnion samples were collected from the microsurface layer (top 60 µm) and from the epilimnion (top 10 m). The microsurface layer was sampled as previously described (36) using a custom-modified towed skimming sled. Integrated water samples from the top 10 m (epilimnion) were collected using a modified zooplankton pump by the method of Frischer et al. (14). Samples were kept on ice, transported to the lab, and processed within 3 h of collection.
|
View this table: [in a new window] |
TABLE 3. Distribution of clones in each lake sample
|
DNA purification, PCR amplification, cloning, and sequencing.
AEAP lakes were sampled for bacterioplankton from 10 to 24 September 2002. Lake George was sampled in November 2000. For AEAP lakes, epilimnetic and hypolimnetic subsamples (900 to 1,000 ml each) were collected and prefiltered (5 µm), and bacterioplankton communities were then collected on 0.22-µm mixed cellulose ester filters (Millipore) by gentle vacuum using Nalgene MityvacII hand pumps and stored at –80°C until further processed. For all lakes, total genomic DNA was purified using the soil DNA extraction kit (Q-Biogene) with a bead-beating step using filters cut into small pieces as the "soil." Purified DNA was quantified by spectrophotometry and analyzed for quality by visual assessment after agarose gel electrophoresis on 1% gels. Total genomic DNA was used for PCR amplification of the 16S rRNA genes using universal forward 8F (5' AGA GTT TGA TCM TGG CTT CAG) and reverse 1492R (5' GGT TAC CTT GTT ACG ACT T) primers (1). Amplification products were resolved by gel electrophoresis and subcloned into the pCR 2.1 vector using the TA Cloning kit (Invitrogen). After the amplification product was cloned into the plasmid vector, the 16S rRNA gene insert was reamplified with primers M13 Forward (–20) (5' GTA AAA CGA CGG CCA GTG) and M13 Reverse (–27) (5' GGA AAC AGC TAT GAC CAT G) targeting the M13 cloning site. The amplified insert was verified visually by agarose gel electrophoresis (1% gels) and purified using the Qiagen MinElute purification kit following the manufacturer's instruction (Qiagen). Purified PCR amplicon (typically 5 to 125 ng µl–1) was eluted in 10 µl EB buffer (10 mM Tris-Cl, pH 8.5). Sequencing was accomplished using the Big Dye Terminator v3.0 cycle sequencing kit (Applied Biosystems) in-house starting with 1 to 2 µl of purified PCR amplicon using the M13 primers. Automated sequencing was performed using an ABI 3700 sequencer at the Center for Functional Genomics at the State University of New York at Albany. The identity of each clone sequence was determined by identifying its nearest neighbor using the Sequence Match and Classifier tools at the Ribosomal Database Project (RDP) (http://rdp.cme.msu.edu/) on 18 May 2006 (RDP update 39).
Clone identification.
A total of 1,907 clones from 25 libraries produced from 17 AEAP study lakes and 228 clones from Lake George in six libraries were isolated and analyzed. Because all of the clones from the Lake George libraries were derived from epilimnetic samples, they were pooled and analyzed as a single library. Initially, 96 clones were picked from each library but depending on downstream processing, fewer clones were usually analyzed. For the AEAP lakes, the number of clones successfully analyzed from each library ranged from 25 (Big Moose Lake, epilimnion) to 97 (Hell Diver Pond, hypolimnion; 97 because there was low recovery from the first library, a second library was also processed). The average number of clones analyzed per lake was 76 ± 19, and the specific number of clones analyzed in each lake is presented below (see Table 3). Initially, amplified rRNA gene inserts from clones were verified visually by agarose gel electrophoresis (1% gels) to confirm the presence of full-length rRNA gene inserts. Clone inserts were preliminarily identified on the basis of restriction digestion profiling of amplified 16S rRNA gene fragments (amplified rRNA gene restriction analysis [ARDRA]) essentially as previously described by Vergin et al. (56) except that the isoschizomer of BsuRI, HaeIII was utilized. ARDRA patterns were assessed after visualization by agarose gel electrophoresis on 3% agarose gels. Restriction profiles were visually sorted into patterns differentiated by the number and size of fragments. Pattern recognition was further facilitated by digital image analysis using the Kodak 1D Image Analysis Software Package v3.6 (Eastman Kodak, NY). ARDRA samples were first pooled by lake layer and then between epilimnion and hypolimnion layers of the same lake. Each unique ARDRA profile was assigned an identification code, and on average, 62% of the representative clones of each ARDRA pattern were sequenced. In the case of the same ARDRA pattern being discovered in different lakes, representative clones from each lake were always sequenced, which provided identical taxonomic classifications at the class and subclass levels as determined using the RDP Classifier tool. However, in some cases when the same ARDRA pattern was observed in a single lake, but in different layers, representatives from both layers were not always sequenced. Of the 1,907 clones obtained from the AEAP lakes and classified into 86 unique ARDRA profiles, 431 were sequenced. There was evidence of microdiversity among clones that had been classified by ARDRA (i.e., multiple ARDRA patterns were identified as belonging to the same class or subclass and very minor sequence differences between some clones were observed). From the Lake George libraries, 228 clones were obtained, all of which were sequenced. In 10 instances, when the ARDRA patterns were very distinct, less than 20% of the clones with identical ARDRA patterns were sequenced.
Each clone was phylogenetically classified to the taxonomic level of class and subclass according to the hierarchical taxonomy proposed by Garrity et al. (15) using the RDP Classifier tool (http://rdp.cme.msu.edu/classifier/classifier.jsp) release 9.0. The taxonomic classification of subclass is an intermediate taxonomic rank between class and order used by Garrity et al. (15) and adopted by the RDP. Although the majority of sequences were reliably classified to much lower taxonomic ranks, it was not possible to consistently assign all 2,135 clone identities to less than a class designation. Therefore, the generic term subclass was used to define classifications below the rank of class but higher than order. A comprehensive phylogenetic analysis of the Adirondack lake bacterial clones obtained in this study was beyond the scope of the current study and will be presented elsewhere. Likewise, since all clones could not be assigned to taxonomic ranks directly comparable to the recently proposed ecologically significant classifications of the freshwater betaproteobacteria (19) and Actinobacteria (58), comparison on the basis of these rankings was also considered outside the scope of this study.
Database organization.
All clone information is maintained in an Excel spreadsheet and identified by a unique clone identification number consisting of an ADK (Adirondack) designation code, lake identifier, date (year) identifier, and clone number. Epilimnetic and hypolimnetic samples are designated by an e or h, respectively, following the lake identifier. The lake identifier codes are shown in Table 1. For example, the clone identifier ADK-DTe02-58 uniquely identifies clone 58 obtained from the epilimnion of Dart Lake in 2002. Complete clone information, sequences, and classification are maintained on the Darrin Fresh Water Institute website (www.rpi.edu/dept/DFWI/index.html) and are available for research purposes with permission.
Statistical analyses.
Richness at the taxonomic levels of class and subclass in all clone libraries was estimated using the software package EstimateS v7.5.1 (R. K. Colwell, EstimateS [http://purl.oclc.org/estimates]), which gives a statistical estimation of species richness and shared species from samples. The richness estimator Chao2 without bias correction is reported, since this metric is recommended as among the most robust for sample sets with many rare species (5) although the Chao1, ICE, and ACE indices were also calculated and resulted in similar richness estimates. Richness of the complete Adirondack data set was estimated on the basis of the total number (2,135 clones) of phylogenetically classified clones (Table 1). Sampling sufficiency of each library (lake sample) was determined as described by Kemp and Aller (27) using the "Large Enough" estimator available online at http://www.aslo.org/lomethods/free/2004/0114a.html. The Large Enough software provides a convenient calculator for the Chao1 and ACE richness estimates based on the number of uniquely classified clones in each clone library (see Table 3). The actual clone numbers for each lake sample can be calculated by multiplying the number of clones per sample (n) by the percent occurrence of each bacterial class. Diversity indexes including the Shannon diversity index, the Brillouin index, Simpson's index, McIntosh's index, and the Berger-Parker index were all calculated using the BIO-DAP software package (G. M. Thomas [http://nhsbig.inhs.uiuc.edu]), a biodiversity analysis package, although only the Shannon diversity index is reported. Diversity indices were calculated on the basis of the number of uniquely classified clones in each clone library. The overall bacterioplankton community composition was explored by principal component analysis (PCA) using the Unscrambler software package v9.2 (CAMO, Norway). Principal component models were built with full leverage correction and based on percent normalized (by lake sample) abundance of clones in each library (see Table 3). The statistical significance of PCAs were assessed by analysis of variance of PCA scores between lake samples. The significance of correlations between bacterial classes and the 31 physical and chemical properties was assessed by correlation analysis between each chemical and physical parameter and frequency of occurrence of each phylogenetic bacterial class in each lake sample. In cases where normality and equality of variance of least square residuals assumptions were satisfied, the parametric Pearson correlation procedure was used. In cases where normality and equality of variance assumptions were not met, the nonparametric Spearman rank procedure was utilized. All analysis of variance and correlation analyses were facilitated using the software package SigmaStat v3.0 (SPSS Inc.).
Nucleotide sequence accession numbers.
The sequences of all unique clones obtained in this study have been deposited in GenBank (accession numbers EF520350 to EF520639).
|
|
|---|
![]() View larger version (20K): [in a new window] |
FIG. 1. Map of New York State with the Adirondack Park showing the 18 lakes sampled during this study. The locations of the 18 lakes sampled during the study are indicated (). The inset indicates the location of New York State within the United States.
|
|
View this table: [in a new window] |
TABLE 2. Classified rRNA gene clones from Adirondack lakes
|
Of the 476 Actinobacteria clones, 207 were most similar (98.3% nucleotide similarity) to an uncultured Actinobacterium clone (GenBank accession number AJ575501) belonging to the acI-A clade as first described by Warnecke et al. (58) and originally detected in the
imov Reservoir, a eutrophic dimictic reservoir located in southern Bohemia. The second most abundant group (116 clones) was most similar (92.5% nucleotide similarity) to an uncultured Actinobacterium clone FukuS81 (GenBank accession number AJ290054). This clone (FukuS81) was described by Glöckner et al. (17) and originally identified from a water sample collected from Lake Fuchskuhle, a meso- to acidotrophic forest lake in the Brandeburg Mecklenburg lake district (Germany). According to the phylogenetic analysis presented by Warnecke et al. (58), this organism belongs to the acI-B cluster of the freshwater Actinobacteria. At the subclass level, both of these clones were classified as unclassified Actinomycetales using the RDP classifier tool release 9.0 (http://rdp.cme.msu.edu/classifier/classifier.jsp). The remaining 153 Actinobacteria clones were distributed among seven additional subclasses with nearest neighbors to 41 different uncultured clone sequences. However, none of these subclasses were represented by more than 16 clones, and in most cases, these subclasses were represented by three or fewer clones. GenBank accession numbers of all unique clones recovered and sequenced in this study are provided in Table 2.
Richness, diversity, and similarity of Adirondack lake bacterial communities.
The overall richness of bacterial communities in the study lakes was estimated using the Chao2 estimator of richness (4) (Fig. 2). Total richness at the phylogenetic class and subclass levels was estimated to be 27 ± 6.5 and 149 ± 12, respectively. In both analyses, saturation of the Chao2 estimator was reached. Several other richness estimators, including the ACE, ICE, Chao1, and Jack estimators were also calculated with similar results. Richness estimates at the class level ranged from 23 to 31 and from 120 to 173 at the subclass level in these analyses.
![]() View larger version (8K): [in a new window] |
FIG. 2. Richness estimates of 16S rRNA gene fragment clone collections from AEAP lakes and Lake George at the phylogenetic level of class (a) and subclass (b). A total of 2,135 individual clones were classified using the RDP Classifier tool (http://rdp.cme.msu.edu/classifier/classifier.jsp) from 26 samples derived from 18 different Adirondack lakes. The richness estimator Chao2 was calculated as described by Chao (4) and implemented in the EstimateS software package v7.5.1 (http://purl.oclc.org/estimates). The Chao2 richness estimator was utilized without bias correction.
|
![]() View larger version (32K): [in a new window] |
FIG. 3. Sampling sufficiency estimates of clone libraries from AEAP study lakes and Lake George. The "Large Enough" calculator (27) was used to determine whether individual clone libraries were sampled sufficiently at the class level. If the estimated phylotype richness reached an asymptote, we inferred that the library was large enough to yield a stable estimate of phylotype richness. With the exception of the libraries derived from the Big Moose Lake and Sagamore Lake epilimnion samples (b), all other samples appeared to have been sufficiently sampled.
|
![]() View larger version (9K): [in a new window] |
FIG. 4. pH versus estimated diversity (a and c) and richness (b) of bacterial community composition based on 16S rRNA gene fragment clone collections from AEAP lakes and Lake George at the phylogenetic level of class (a and b) and subclass (c). Richness was not estimated at the subclass level because sampling at this level was not sufficient. Significant correlation between diversity and pH (r = 0.387, P = 0.05) (a) and between richness and pH (r = 0.5, P = 0.01) (b) were observed. Trends at the subclass level were similar, but statistical correlations were nonsignificant (r = 0.3, P = 0.13). Diversity was estimated using the Shannon index as implemented in the BIO-DAP software (http://nhsbig.inhs.uiuc.edu). Richness for each lake sample was estimated using the Chao1 statistic as implemented using the "Large Enough" estimator (27). Richness estimates for libraries derived from the epilimnion of Big Moose Lake and Sagamore Lake were omitted because analysis suggested that these libraries were insufficiently sampled. A total of 2,135 individual clones were classified using the RDP Classifier tool (http://rdp.cme.msu.edu/classifier/classifier.jsp) from 26 samples derived from 18 Adirondack lakes. RDP classifications were assigned on 18 May 2006 (RDP update 39).
|
![]() View larger version (19K): [in a new window] |
FIG. 5. (a) Principal component analysis of bacterial community composition based on 16S rRNA gene fragment clone collections from AEAP lakes and Lake George at the phylogenetic level of class. The Actinobacteria, Betaproteobacteria, and Gammaproteobacteria accounted for significant differences (P = 0.008) between samples. A total of 75% of the variation was explained by the first two principal components (the first principal component [PC1] and the second principal component [PC2]). (b) Average distribution of the Actinobacteria, betaproteobacteria (β-proteobacteria), and gammaproteobacteria ( -proteobacteria), and all other bacteria in each PCA-discriminated group. A total of 2,135 individual clones were classified using the RDP Classifier tool (http://rdp.cme.msu.edu/classifier/classifier.jsp) from 26 samples derived from 18 Adirondack lakes. RDP classifications were assigned on 18 May 2006 (RDP update 39).
|
![]() View larger version (8K): [in a new window] |
FIG. 6. Principal component analysis of bacterial community composition based on 16S rRNA gene fragment clone collections from AEAP lakes and Lake George at the phylogenetic level of subclass. At the phylogenetic level of subclass, most differences could be attributed to the unclassified betaproteobacteria and an unclassified Actinobacteria subclass, but differences were not significant (P > 0.05). The first two principal components (the first principal component [PC1]) and the second principal component [PC2]) explained 57% of the overall variation. The clone library from the epilimnion of Middle Settlement Lake was omitted from the analysis as an outlier. A total of 2,135 individual clones were classified using the RDP Classifier tool (http://rdp.cme.msu.edu/classifier/classifier.jsp) from 26 samples derived from 18 Adirondack lakes. RDP classifications were assigned on 18 May 2006 (RDP update 39).
|
![]() View larger version (8K): [in a new window] |
FIG. 7. Relationship between the percentage of Alphaproteobacteria and Bacteroidetes in each lake clone library and pH. The percentage of 16S rRNA gene clones in each lake sample clone library is shown. (a) Relationship between the percentage of alphaproteobacteria ( -proteobacteria) and pH (r = –0.78, P = 0.0002); (b) relationship between the percentage of Bacteroidetes and pH (r = 0.48, P = 0.02).
|
0.05) of bacterial classes with the physiochemical properties of the Adirondack lakes in this study. Of the most abundant bacterial classes, the Alphaproteobacteria, Actinobacteria, and Betaproteobacteria classes were significantly correlated with lake physiochemical properties. In addition to these properties, the abundance (proportion) of Alphaproteobacteria was also significantly correlated with total phosphorus (r = 0.5, P = 0.01), molybdate-reactive phosphorus (MRP) (r = 0.6, P = 0.03), chlorophyll a (r = 0.6, P = 0.0007), total monomeric aluminum (r = 0.4, P = 0.05), and nonlabile monomeric aluminum (r = 0.4, P = 0.04). The abundance of Actinobacteria was positively correlated to MRP (r = 0.4, P = 0.01). The abundance of Betaproteobacteria was negatively correlated (r = –0.5, P = 0.05) with lake hydraulic retention times. However, since most of the lakes had relatively short retention times, it is difficult to interpret these results. Second to properties directly related to acidity (pH and acid-neutralizing capacity), the concentrations of reactive phosphorus appeared to most significantly explain the abundance and distribution of the Actinobacteria and Alphaproteobacteria within the Adirondack lakes. In the complete data set, there were 81 significant correlations between lake physiochemical properties and the abundance of specific groups of bacteria (see Table S2 in the supplemental material), but the full exploration of these correlations are beyond the scope of this report. |
|
|---|
Bacterial communities in the Adirondack lake set were, in many ways, similar to bacterioplankton communities described in other freshwater systems. Communities were dominated by "typical" freshwater assemblages consisting of Betaproteobacteria and Actinobacteria lineages (54, 55, 60). The dominant groups of these bacteria observed in the Adirondack lakes have been previously described and phylogenetically classified, and some have even been isolated (19, 33, 38, 58). In 15 of the 26 lake samples examined, the Betaproteobacteria were the most abundant group with the largest number of clones (291 clones) belonging to the P. necessarius cluster. Of these clones, 152 (52%) were classified within the PnecC subcluster. The most abundant group within the Actinobacteria was the Actinomycetales subclass belonging to the acI-A clade (391 of 476 clones) as originally described by Warnecke et al. (58). Interestingly, seven lake samples were distinguished by the abundance of Gammaproteobacteria and in two samples (Cascade Pond epilimnion and Ice House Lake hypolimnion), Gammaproteobacteria were the most dominant group. To the best of our knowledge, this is the first report of the dominance of Gammaproteobacteria in any freshwater lake, although Gammaproteobacteria are common in marine systems (16). Two clone libraries (Brooktrout Lake epilimnion and hypolimnion) were dominated by representatives of the Gammaproteobacteria. Dominance of the Alphaproteobacteria is also more typical in marine systems (18, 40), but unlike marine systems that are dominated by the SAR and Roseobacter clades of the Alphaproteobacteria, the Alphaproteobacteria most commonly observed in the Adirondack lakes belonged to representatives of the Caulobacter (22%), Methylocystaceae (19%), and Rhizobiales (14%) taxa.
Several recent studies have investigated the relationship between environmental factors and the composition of bacterioplankton in freshwater lakes and rivers based on direct analysis of amplified rRNA gene and rRNA gene clone libraries (64). These studies have examined a variety of lake types from the tropics to polar regions and have included lakes of various sizes, hydrological types, and trophic status, but only a few have compared communities between large suites of lakes (30, 31, 62). Although these studies included acidity (pH) as a variable, none of them were specifically designed to investigate the influence of pH on bacterioplankton communities. On the basis of these investigations, it has been reported that lake size and hydraulic retention time (30, 31) and regional/landscape factors, including location and hydrological type (32, 62), appear to be primary variables responsible for shaping the composition of freshwater prokaryotic communities.
The association of pH with distribution of freshwater bacterial groups has also been reported (31, 50, 55, 62). Lindström et al. (31) examined 15 diverse lakes in northern Europe and identified that the pH, temperature, and hydraulic retention time were most strongly related to variations in the distribution of bacterial taxa. Similarly, in a study of four shallow eutrophic lakes in Belgium, a significant relationship between pH and the composition of lake bacterioplankton communities was reported (54).
In the Adirondack lakes, the diversity and richness of bacterial communities were significantly correlated with pH, mirroring similar relationships observed between higher organisms and pH in lakes (12, 23, 43, 46). However, pH was not significantly correlated with the overall structure of bacterioplankton communities as investigated by PCA, since the most dominant groups of bacteria (Betaproteobacteria and Actinobacteria) were not correlated with pH. Although these most abundant bacterial classes were not correlated with pH, the Alphaproteobacteria and Bacteroidetes taxa were. In the study of Lindström et al. (31), a relationship between these groups and pH was not reported, but they observed that two actinobacterial groups were correlated with pH. The ACK-M1 group was associated with high pH; Sta2-30 was associated with low pH. A third actinobacterial group (Urk-014) was not correlated with pH. These groups correspond to the acI-A, acI-B, and acIV actinobacterial lineages described by Warnecke et al. (58), respectively. In the present study, in which we recovered 476 Actinobacteria clones, the acI-A and acI-B groups exhibited essentially the same distribution among the lakes and did not correlate with pH. Lindström et al. (31) also reported an association between low pH and betaproteobacteria belonging to the P. necessarius group, which was not observed in the present study.
Possible relationships between other physiochemical parameters and bacterial communities were also examined. Significant correlations between the abundance of different bacterial classes and 26 of the 31 parameters examined were observed (see Table S2 in the supplemental material). Of the more abundant taxa, the Alphaproteobacteria, Bacteroidetes, and Verrucomicrobiae were correlated with a number of lake parameters other than acidity. The Bacteroidetes (largely composed of the unclassified Rikenellaceae subclass) was significantly correlated with 12 parameters (pH, dissolved oxygen, temperature, acid-neutralizing capacity, NH4, dissolved inorganic carbon, SO4, total P, MRP, total filterable P, total N, and Fe). The Adirondack Bacteroidetes clones were most similar to bacteria recovered from other acidic environments, including Lake Fuchskule, a mesotrophic acidotrophic forest lake (17). This group of organisms may be of particular interest for investigating bacteria associated with aquatic systems with low pH. The Verrucomicrobiae taxon was also significantly correlated with four parameters (water clarity, DOC, TFP, and total aluminum). This group of organisms has been reported from a wide variety of environments (10), including lakes and rivers, acidic forest soils, Antarctic sediments, and manure leakage, but little is known about the physiology of the Verrucomicrobiae.
The roles of regional and landscape factors have been suggested to be the most important factors influencing the composition of freshwater bacterial communities (31, 32, 62). Specifically, the influence of geographic region and the landscape level factors of lake type (seepage versus drainage) and hydraulic retention time. In this study, the lakes examined were within close proximity to each other and were predominantly drainage lakes; therefore, it was not possible to investigate regional differences or lake hydrological type. However, the most dominant group of bacteria observed in the Adirondack lakes, the Betaproteobacteria, was correlated with hydraulic retention time, supporting previous reports.
The correlation of some bacterial groups, particularly the Alphaproteobacteria, may be useful as biological indicators of recovery from acidification in Adirondack lakes. Initially, it was hypothesized that because of their short generation times and small size, bacterioplankton might be excellent as rapid biological indicators of recovery from acidification. However, relatively weak and inconsistent correlations between bacterial community and lake physiochemical properties suggest that bacterioplankton communities are indirectly rather than directly influenced by acid deposition. For example, bacterial growth rates may be influenced by alterations in populations of primary producer communities and chemical effects on nutrient cycling that affect the availability of dissolved organic and inorganic nutrients (13). Alternatively, the composition of bacterial communities may be selectively controlled by the composition of bactivore communities (2). If this were the case, then the restoration of Adirondack bacterioplankton communities to conditions that resemble similar but nonimpacted freshwater lakes might be interpreted as an end point metric in lake recovery rather than an early indicator of ongoing biological recovery.
Preliminary comparisons between a small data set collected in 1994 (22, 33, 34) and the 2002 studies described here (comparison of five lakes) are potentially useful for this purpose. Over this 8-year period, the AEAP study lakes experienced an average pH increase of 0.6 ± 0.3, while the pH of Lake George remained circumneutral (S. A. Nierzwicki-Bauer, unpublished). Although a robust analysis between these data sets is not possible because of methodological differences and the small sample size of the 1994 data set, several trends were observed. Bacterioplankton communities in the AEAP study lakes changed significantly from 1994 to 2002, while those in Lake George did not. Further, in the AEAP lake set, the abundance of Alphaproteobacteria changed, although this change was not significantly correlated with acidity. Interestingly, a significant negative correlation (r = –0.4, P = 0.02) between the abundance of Gammaproteobacteria and pH was observed in the lakes studied in 1994 and 2002 (data not shown), which was not suggested by examination of the 2002 data alone (18 study lakes). The significance of these preliminary observations with respect to recovery from acidification is not yet clear, but ongoing larger-scale studies are under way to explore this process.
Although it is difficult to reach conclusions regarding the functional significance or mechanisms responsible for relationships between acidity and the composition of bacterioplankton communities based only on ribosomal markers, the results of this study indicate that at least some groups of freshwater bacterioplankton are influenced by parameters associated with acid deposition. We observed that it was not the most abundant groups of bacterioplankton that were affected by pH, but several less dominant groups. In the future, it will be important to investigate the physiological and functional roles of these organisms to make predictions concerning potential biogeochemical and ecological significance. Perhaps in freshwater environments, the Betaproteobacteria and Actinobacteria can be considered ecological generalists, while other groups are more specialized. If this is the case, loss or gain of specialized communities may be of importance to particular biogeochemical or other processes.
Although this work has been funded by the U.S. Environmental Protection Agency, it has not been subjected to the agency's review and, therefore, does not necessarily reflect the views of the agency, and no official endorsement should be inferred.
We are greatly indebted to our colleagues in the AEAP that have made significant contributions to the collection, analysis and/or interpretation of samples and results. In particular, Robert Bombard, Larry Eichler, Laurie Ahrens-Franklin, James Harrison, Jim Sutherland, Alexey Vepritskiy, and David Winkler made significant contributions. We also gratefully acknowledge the outstanding contributions of numerous undergraduate students in carrying out the molecular studies. Katie Opper, Mike Schaaff, Christy Stagnar, and Claudia Weigl particularly made outstanding contributions. Anna Boyette (Skidaway Institute of Oceanography) prepared the figures. We thank the useful and constructive comments of two anonymous reviewers.
Published ahead of print on 1 February 2008. ![]()
Supplemental material for this article may be found at http://aem.asm.org/. ![]()
Present address: Molecular and Microbiology Department, Regeneron Pharmaceuticals, 81 Columbia, Turnpike, Rensselaer, NY 12144. ![]()
Present address: Albany Medical College, Center for Immunology and Microbial Disease, 47 New Scotland Ave., MC-151, Albany, NY 12208. ![]()
¶ Present address: Department of Chemistry and Chemical Biology, Rensselaer Polytechnic Institute, 110 8th Street, Troy, NY 12180. ![]()
|| Present address: Department of Pharmacology and Physiology, University of Rochester Medical Center, 601 Elmwood Ave., Box 633, Rochester, NY 14642. ![]()
|
|
|---|
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»