Previous Article | Next Article ![]()
Applied and Environmental Microbiology, May 2008, p. 3014-3021, Vol. 74, No. 10
0099-2240/08/$08.00+0 doi:10.1128/AEM.01809-07
Copyright © 2008, American Society for Microbiology. All Rights Reserved.

Meredith A. J. Hullar,2,
David A. Stahl,2 and
Louis A. Kaplan3
Department of Microbiology, Miami University, Oxford, Ohio 45056,1 Civil and Environmental Engineering, University of Washington, Seattle, Washington 98108,2 Stroud Water Research Center, Avondale, Pennsylvania 193113
Received 3 August 2007/ Accepted 20 March 2008
|
|
|---|
|
|
|---|
This study examined microbial communities from streambed sediments in low-order, forested streams within three biomes, geographic areas distinguished by climate and their predominant terrestrial vegetation. Forested headwater streams derive most of their organic energy from allochthonous sources, either directly from the adjacent terrestrial vegetation or as products of decaying vegetation modified by soil diagenesis. Our study was designed to investigate the general question of how similar are stream microbial communities within and among biomes and, more specifically, do heterotrophic bacteria within streambed communities exhibit biogeographic patterns at the biome level? Nine streams, three located in each of three biomes, were assayed for bacterial abundance, microbial biomass, and microbial and bacterial community structures using a combination of classical, biochemical, and molecular methods. Multivariate statistical analyses were performed to compare the patterns of community structure within and among biomes.
|
|
|---|
Environmental variables.
Concomitant with benthic sampling, temperature was measured with a hand-held field thermometer and water was collected for analyses of dissolved organic carbon (DOC), anions, cations, pH, and conductivity. Samples for DOC analyses were collected in borosilicate glassware rendered organic C free by combustion (500°C for 6 h), and then filtered through precombusted glass fiber filters (Whatman GF/F) with a syringe and syringe-type filter holder (32). DOC concentration was determined by Pt-catalyzed persulfate oxidation with either an OI 700 or OI 1010 total organic carbon analyzer. Anions and cations were measured with a Dionex DX 500 ion chromatography system equipped with an ED40 electrochemical detector after passing each sample through a sterile 0.22-µm Pall Gelman HF Tuffryn acrodisc filter. Conductivity was determined with a YSI model 32 conductance meter, and pH was measured with a Fisher Scientific pH probe and meter.
Sampling of benthic microbial communities.
Sediments were delimited with a 2- by 10-cm (height by inside diameter) Plexiglas ring that was inserted 2 cm deep into the sediment. Plexiglas plates were slipped under and over the ring, effectively trapping the sediments and allowing them to be removed from the stream without disturbance. Sediments in the upper 2 mm were transferred with a methanol-washed spatula to a clean plastic weigh boat and thoroughly mixed. Sediments were subsampled for bacterial abundance, lipid analysis, and molecular analysis (approximately 0.5, 1.0, and 0.5 g [wet weight] of sediment, respectively), and each subsample was placed in separate 2-ml microcentrifuge tubes. This procedure, beginning with delimiting of the sediment, was repeated twice to yield three true replicate samples from each stream. Tubes receiving subsamples for molecular analysis contained 0.5 g of prebaked (180°C for at least 2 h) zirconium beads (0.1-mm diameter; Biospec). Subsamples were preserved (2.5% formaldehyde for bacterial abundance and freezing for lipid and molecular analyses) and shipped to the appropriate laboratory for further analysis. Several subsamples designated for molecular analysis were lost in shipping or analysis.
Epifluorescence microscopic counts.
Bacterial abundance was estimated by epifluorescence direct microscopic counts (28). Briefly, preserved sediments were treated with 0.1 mM sterile tetrasodium pyrophosphate (56) for 10 min, and cells were removed by sonication (75 W, 60 s). Suspensions were vortexed, and subsamples (1 ml) were transferred into sterile vials containing 1 ml of tetrasodium pyrophosphate, thoroughly vortexed, and sonicated (75 W, 30 s). Samples were mixed with sterile glycerol (30% final concentration) in 2.5% formaldehyde, vortexed, and spun at low speed (3 min, 25 x g) to pellet bacterium-free particles. After staining with propidium iodide, bacterial cells in the supernatant fluid were filtered onto a black 0.2-µm-pore-size Nuclepore filter and 20 randomly selected fields or 300 to 500 cells were enumerated (4). One filter was counted per sediment sample, with duplicate filters counted for 10% of the samples.
Phospholipid analysis.
Microbial biomass and community structure were determined using phospholipid analysis following the methods of Findlay (19). Briefly, lipids were extracted from frozen sediment samples by dichloromethane-methanol. The extraction mixture was partitioned into aqueous and organic fractions through the addition of dichloromethane and water, and the organic fraction containing the lipids was recovered. The organic fraction was subsampled for total microbial biomass, determined by phospholipid phosphate (PLP) (21). Phospholipid fatty acids (PLFAs) were recovered from the remaining lipid by differential elution from silicic acid columns (J. T. Baker) and were analyzed as their methyl esters (19). Fatty acid methyl esters were quantified by gas chromatography, with identification based on relative retention times, coelution with known standards, and mass spectrometry analysis. Using polyenoic fatty acids as indicators of microeukaryotes, total biomass was divided between prokaryotic and microeukaryotic organisms (20) and the results were expressed as percentages. To investigate bacterial community structure, PLFAs assigned a priori to the functional group microeukaryotes (19) and those known to be common to both bacteria and microeukaryotes were removed from the PLFA profiles (58), leaving only PLFAs of bacterial origin, and weight percentages were recalculated. Prokaryotic biomass was converted into bacterial abundance using a conversion factor of 100 nmol PLP = 4 x 109 cells (3).
DNA extraction.
Phosphate and MT buffers from the Fast DNA spin kit for soils (Qbiogene, Carlsbad, CA) were added to the sample tubes containing the frozen sediment and zirconium beads. The samples were shaken in a Bio101 FastPrep (Qbiogene, Carlsbad, CA) at speed 4.5 for 15 s (twice) and placed on ice for 1 min between mechanical disruptions. The samples were centrifuged at 4°C at 15,000 x g for 5 min. Supernatants were placed into new tubes, and samples were processed according to the manufacturer's protocol. Genomic DNA was eluted in 50 µl diethyl pyrocarbonate (Ambion, Austin, TX)-treated water and stored at –80°C until further analysis. To evaluate quality, the DNA was resolved on 0.8% (wt/vol) high-melt agarose gels with 1x TAE (Tris-acetate-EDTA) buffer and visualized by ethidium bromide staining (48).
tRFLP analysis.
For terminal restriction fragment length polymorphism (tRFLP) analysis, bacterial 16S rRNA gene sequences were amplified using primer S-D-Bact-1512-a-A-21 (35) and a second primer (S-D-Bact-008-a-S-17) 5' end labeled with 6-carboxyfluorescein (35). Each 50-µl PCR mixture contained primers at 0.7 µmol liter–1, 1.0 µmol liter–1 Tris-HCl (pH 8.8), 5.0 µmol liter–1 KCl, and 0.15 µmol liter–1 MgCl2; 125 µmol liter–1 of each deoxynucleoside triphosphate (Invitrogen, Carlsbad, CA); 25 µg of bovine serum albumin; and 2.5 U Taq DNA polymerase (Invitrogen, Carlsbad, CA). Each reaction mixture was incubated on a PTC-100 thermal cycler (MJ Research, Waltham, MA) using the following "touchdown" parameters: initial denaturation at 94°C for 1 min, followed by four cycles of denaturation at 94°C for 30 s, annealing at 65°C for 30 s, and extension at 72°C for 30 s. The annealing temperature for the reaction was subsequently reduced by 1°C per cycle until a 50°C annealing temperature was reached, cycled five more times at a 50°C annealing temperature, and terminated with a final extension at 72°C for 5 min. The PCR products were visualized on a 0.8% (wt/vol) high-melt agarose gel stained with ethidium bromide and quantified based on relative intensity using a quantitative DNA ladder (Invitrogen, Carlsbad, CA) and NIH image v 1.63 (http://rsb.info.nih.gov/nih-image).
PCR products, derived from environmental samples, were digested overnight in the dark at 37°C using 10 U of the restriction enzyme HaeIII, which has a 4-bp recognition site, in a standard restriction buffer (One-Phor-All Plus; Amersham-Pharmacia). Digested samples were ethanol precipitated (48), dried, and analyzed using an ABI 377 DNA sequencer (Applied Biosystems, Inc., Fremont, CA). At least 100 fmol of each digested sample was loaded onto a 4% polyacrylamide gel for fragment analysis in Gene Scan mode. The Genescan TAMRA 500 size standard (HaeIII) was added to each sample lane. ASCII files of electropherograms were analyzed using Dax analysis software (van Mierlo, Inc., Holland). Peaks representing tRFLP fragments were measured as integrated peak areas. Peaks less than 2 bp apart were binned and normalized as the relative percentage of total peak area. Peaks less than 1% of the total area or found in less than three replicate profiles were excluded from further analysis.
Statistical analysis.
Comparisons of water chemistry, microbial biomass, and bacterial abundance were made using nested analyses of variance (streams within and among biomes) (Systat 10 and SAS). Relationships among variables were investigated using linear regression (Microsoft Excel 2004 for Mac 11.2). Fatty acid profiles were subjected to principal component analysis (PCA) after log transformation [ln(x + 1)] of weight percent fatty acid data (SPSS, version 7.5, and Minitab, release 12.23); PLFA profiles were interpreted with a functional group approach (19).
Nonmetric multidimensional scaling, which does not require the assumption of linearity, was used to compare tRFLP fragment composition among streams (PC-ORD version 4) (38) following the calculation of Sorensen's distance (34). The raw data set consisted of 15 samples made up of 16 fragment lengths or OTU. To reduce skew, the square root of the relative percent area was arcsine transformed (52). For the ordination, the autopilot option was set to "slow and thorough." Choosing the axis that minimizes stress and maximizes the interpretation of the data assessed the amount of variability explained in the data set. Stress, an indicator of goodness of fit, measured the inverse of the fit between the original data matrix and the reduced dimension ordination matrices. Stress values between 0 and 15 give a good approximation of the data in multivariable space with a low risk of drawing false inferences (38). A plot between the final stress and the number of iterations was used to assess the stability of the solutions. Monte Carlo tests identified the dimensions that gave solutions that were significantly different from those due to random chance. The proportion of the variance of the distances between the original matrix and the ordination space explained was described for each axis.
|
|
|---|
0 to 21°C), while the tropical evergreen streams had nearly constant temperatures (range, 19 to 22°C). Biome-level average DOC concentrations at baseflow were lowest in the tropical evergreen streams (0.7 mg C/liter), twofold higher in the deciduous forest streams (1.5 mg C/liter), and highest in the coniferous forest streams (8.6 mg C/liter). |
View this table: [in a new window] |
TABLE 1. Summary of chemical characterizations of stream water and sedimentsa
|
10- to 160-fold greater in southeastern coniferous biome streams than in eastern deciduous and tropical evergreen streams (P < 0.001); there were no significant differences among streams within biomes (P = 0.22) (Table 2). The range of total microbial biomass for streams within the eastern deciduous biome overlapped the range observed from streams within the tropical evergreen biome. Prokaryotes comprised between 39 and 88% of total microbial biomass, and in seven of the nine streams, bacteria accounted for greater than 60% of total biomass. |
View this table: [in a new window] |
TABLE 2. Microbial biomass and bacterial abundance in sediments of nine streams
|
Microbial community structure.
PCA of PLFA profiles based on a data set of 48 PLFAs separated the streams within the three biomes from each other at the biome scale, with principal components 1 and 2 (PC1 and -2, respectively) accounting for 66% of the variance (Fig. 1). The samples from all three tropical evergreen streams were tightly clustered, both within and among streams. All three streams in the eastern deciduous biome displayed low within-stream variation but showed high among-stream variation, as West Creek samples clearly separated from samples from the other two streams. We observed tight clustering of all samples from McDonalds Branch in the southeastern coniferous biome, which also clustered with single samples from Shinns Branch and Mt. Misery. The latter two streams, however, displayed high within-stream variation.
![]() View larger version (10K): [in a new window] |
FIG. 1. PCA of variation in sedimentary microbial community composition in the nine streams by PLFA analysis. Each abbreviation represents an individual sample: E, eastern deciduous; T, tropical evergreen; S, southern coniferous; and 1, 2, and 3, individual streams within the biome. The percent variation explained by each axis is indicated in parentheses. Identified fatty acids had component loadings of >0.5 or <–0.5 and exerted strong influence on the pattern of variation among samples along the respective component axes.
|
![]() View larger version (9K): [in a new window] |
FIG. 2. Relationship between PCA factor 1 score and the calculated percentage that microeukaryotes contribute to total microbial biomass for all stream samples.
|
![]() View larger version (11K): [in a new window] |
FIG. 3. Patterns of variation in sedimentary bacterial community composition in the nine streams by PLFA analysis after removal of fatty acids assigned a priori to the functional group microeukaryotes (19) and those known to be common to both bacteria and microeukaryotes (58) from the PLFA profiles. Figure layout and abbreviations are as described in the legend to Fig. 1.
|
![]() View larger version (7K): [in a new window] |
FIG. 4. Patterns of variation in sedimentary microbial community composition in the nine streams by molecular analysis. Shown is a scatter plot of sample scores from nonmetric multidimensional scaling analysis, with each abbreviation representing an individual sample. The percent variation explained by each axis is represented in parentheses. Abbreviations are as described in the legend to Fig. 1. Identified fragment lengths had Pearson's correlation coefficients of >0.5 or <–0.5 and exerted strong influence on the pattern of variation among samples along the respective component axes. Samples not plotted were not sampled or had been sampled but lost in transport or processing.
|
|
View this table: [in a new window] |
TABLE 3. Monte Carlo test of stress in relation to dimensionality
|
|
|
|---|
3, 20:5
3) (11) in describing the variation in microbial community structure among samples from these streams (Fig. 1). Similarly, the separation of White Clay Creek (E1) and Birch Run (E2) from West Creek in the eastern deciduous biome is consistent with the twofold-greater eukaryote contribution to microbial biomass in E1 and E2 and the importance of a lipid marker for algal chlorophytes (18:3
3) (11) in describing the variation in microbial community structure quantified by PC1 (Fig. 1). A high correlation of the ratio of prokaryotic to eukaryotic biomass to variation in microbial community structure quantified by PCA has been reported in microbial communities from shallow, subtidal marine sediments (22), stream sediments (36, 53), and freshwater reservoir sediments (51), but those studies were limited to comparisons within a single ecosystem where changes in community structures were attributed to seasonality. This is the first study, to our knowledge, in which the sampling regimen eliminated the seasonal component of variation and the prokaryotic/eukaryotic biomass ratio influenced microbial community structure at the biome level. Despite the positive correlation between algal biomass and bacterial abundance reported for other stream communities (45) and observed in our study, it is clear that in our study streams algal biomass did not control bacterial community structure. In fact, when microeukaryotic lipids were removed from the analysis, the within- and among-stream variation for each biome was greatly reduced for those biomes with higher algal contributions (southern coniferous and eastern deciduous) and remained mostly unchanged for the highly shaded tropical evergreen biome streams with low levels of algal biomass. Clearly, the major variation in sedimentary bacterial community structure occurred at the biome level, and the variation among streams within a biome was comparable to the variation observed within individual streams. This pattern was evident in both the tRFLP and bacterial PLFA descriptions of community structure and supports the contention, based on a study of bacterial communities within coastal marine waters, that biogeographical patterns are generated by processes that are comparable across major global biomes (42).
Several studies have used PLFA profiles to investigate longitudinal variation within streams (53) or the impacts of anthropogenic stress on microbial community structure (e.g., see references 5 and 36), but we are unaware of any other studies comparing variation in sedimentary microbial and bacterial community structure across streams. Two studies have used the PLFA method to investigate spatial variation in community structure in marine sediments. Federle et al. (18) found significant spatial variation at the centimeter and 100-m levels within subtropical estuarine sediments. Baird et al. (2), using a replicated sampling scheme unprecedented for deep-sea research, found a similar pattern at a deep-sea site subjected to high-energy currents. Several studies have been conducted using molecular approaches to examine spatial variation in benthic bacterial community composition. Gao et al. (25) examined a single stream in each of nine states and reported high within-stream (upstream versus downstream) and among-stream variation. In some cases, variation could be correlated with differences in environmental conditions, but no overall geographic pattern was reported. Braker et al. (7) found distinct bacterial community structures, based upon unique distributions of 16S rRNA gene tRFLP, in marine sediments from Puget Sound and the Washington margin, although the study is based on four cores, limiting the information available to regional spatial scales. A study focused on denitrifying bacterial communities in marine sediments found the most variation in community composition at the kilometer scale (the largest examined), whereas relative abundance of similar OTU became important for differentiating communities at the centimeter scale (49). In our study, within-stream samples characterize variation at the meter scale and among-stream samples characterize variation at the kilometer scale. Clearly, the occurrence of the principal component of variation at the biome level indicates that the greatest variation in bacterial community structure in streams occurs at the largest spatial scales and suggests the development of biogeographic patterns within these communities.
The contrasting patterns of microbial and bacterial community structure suggest that microeukaryotes and bacteria are responding to different environmental determinants. The three geographic locations differ in important physical, chemical, and biological parameters, including terrestrial plant communities, which, in turn, lead to differences in the terrestrially derived, dissolved organic matter quantity and quality in stream water (24, 33). Which of these determinants are the proximate causes of the observed bacterial distribution is not known at this time.
The total microbial biomass for the eastern deciduous and tropical evergreen stream sediments was similar to total microbial biomass for several other temperate streams (36, 53) and within the range reported for freshwater sediments (12). Sedimentary biomass in White Clay Creek sediments was similar to that reported 15 years prior, while biomass in West Creek was elevated in the current study (6). In contrast, total microbial biomass for the southeastern coniferous stream sediments was higher than the range previously reported by Dobbs and Findlay (12) and comparable to that in sediments from a seasonally anaerobic, depositional zone within a freshwater reservoir (51). Bacterial abundances were also greatest in the highly organic sediments of the southeastern coniferous streams, and these values were similar to those observed for benthic organic matter in streams (17, 23). The discrepancy between direct microscopy and biochemical measurements of abundance was likely due to two issues: the necessity of converting biomass, measured as moles of phospholipid phosphate, to number of cells and the difficulty in obtaining a clean, debris-free preparation for microscopic counts. The conversion from biomass to cell number was based on cultured bacteria, particularly Escherichia coli (3). Assuming the average size of the cultured bacteria used in the study by Balkwill et al. (3) was a rod 1.5 µm in length and 1 µm in diameter, a 2.5 times reduction in size to 0.6 µm by 0.4 µm (length by diameter) would yield a biomass-to-cell number conversion factor that would our bring estimates of abundance for eastern deciduous and tropical evergreen stream sediments to unity. The lack of correlation between the two measures of abundance for southeastern coniferous stream sediments is likely due to the high organic content of the sediments causing difficulty in obtaining clean, debris-free preparations for microscopic counts.
Biogeography of aquatic bacteria is in its infancy, and a globally consistent pattern, if one exists, has yet to emerge. Different communities have been observed in a comparison of freshwater lakes and pelagic marine habitats (13), as well as for lakes and streams within an Arctic tundra catchment (9). Within coastal marine sediments, spatially isolated assemblages occur over scales of a few kilometers, with environmental conditions selecting for the individual assemblages (27). A review of studies from a wide range of terrestrial and aquatic environments concluded that free-living bacteria exhibit biogeographical patterns that are the result of both environmental selection and the legacies of historical events (37). It has been suggested that different groups of microbes may show very different biogeographies (14). In all of the recent investigations, streams have been largely ignored. The emergence of clear biome-level patterns in our study suggests that further investigations into the microbial communities present in low-order streams may provide clues to the physical, chemical, and biotic factors influencing the biogeography of bacteria.
This research was supported by NSF DEB 9904047 and DEB 0096276.
Published ahead of print on 31 March 2008. ![]()
Present address: Department of Natural Resources and Water, Brisbane, Queensland 4072, Australia. ![]()
Present address: Fred Hutchinson Cancer Research Center, Seattle, WA 98112. ![]()
|
|
|---|
| |||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||
Copyright © 2009 by the American Society for Microbiology. For an alternate route to Journals.ASM.org, visit: http://intl-journals.asm.org | More Info»