Previous Article | Next Article ![]()
Applied and Environmental Microbiology, January 2006, p. 212-220, Vol. 72, No. 1
0099-2240/06/$08.00+0 doi:10.1128/AEM.72.1.212-220.2006
Copyright © 2006, American Society for Microbiology. All Rights Reserved.
Limnology/Department of Ecology and Evolution, Evolutionary Biology Centre, Uppsala University, Norbyvägen 20, 75236 Uppsala, Sweden,1 Département des Sciences Biologiques, Université du Québec à Montréal (UQÀM), Montréal, Québec, Canada2
Received 16 May 2005/ Accepted 5 October 2005
|
|
|---|
|
|
|---|
The assembly of communities at the local scale is influenced by both local environmental conditions (abiotic and biotic factors) and regional factors (e.g., climate, migration, and speciation), as shown, e.g., for zooplankton (7, 8). Similar processes have been suggested to be important in the regulation of bacterial community composition (32, 54). One view concerning the regulation of local bacterial community structure goes back to the famous statement "everything is everywherethe environment selects" (1). Accordingly, based on the high local but low global diversity of microorganisms, the active community is selected from the total species pool by the prevailing environmental conditions (18). The idea that microorganisms are ubiquitously distributed is supported by studies of protozoan communities (17, 20) and studies reporting on the global distribution of closely related freshwater bacteria (e.g., see references 25, 29, 30, and 55). On the other hand, it is widely recognized that ecological differentiation and adaptation to certain environmental conditions occur even among major bacterial groups, since, e.g., beta-proteobacteria are dominant constituents of freshwater bacterioplankton but very rare in the ocean (24, 45, 55). Further on, the "everything is everywhere" concept is questioned by studies of soils and salt marshes showing increasing numbers of species with increasing sampled habitat sizes, in analogy with established species-area relationships demonstrated for higher organisms (26, 33).
Another view is that the local bacterial community structure is regulated by the size and diversity of the surrounding regional community, or "metacommunity" (11). Curtis and Sloan (11) also arguebased on the neutral theory of biodiversity (34)that local bacterial community composition is a product of random events in connection to the recruitment of functionally equivalent bacterial taxa from the "meta- or source community." Hence, this concept does not assume that there is a one-to-one match between environment (or niche availability) and community structure.
The objective of this study was to compare these contrasting concepts of how bacterial communities are shaped and how bacterial community structure and function are related by using an experimental approach. Different bacterial communities were inoculated into incubation vessels with the same sterile medium, hence making it possible to investigate the importance of the origin of the source community (i.e., the inoculum) versus the environmental conditions (i.e., the medium) for the composition and functional performance of the emerging bacterial communities. Previous studies using dilution cultures often found much stronger effects of the medium than of the inoculum on functional parameters (21, 22, 53), which can be explained either by similarly composed communities being selected by identical media or by differently composed bacterial communities maintaining similar functions (40).
According to the "everything is everywherethe environment selects" concept, initially different local communities are expected to converge towards similar composition upon growth under identical conditions. According to the "metacommunity concept," however, deviating communities consisting of functionally equivalent taxa should be enriched even under identical conditions.
|
|
|---|
|
View this table: [in a new window] |
TABLE 1. Characteristics of habitats from which water for the preparation of inocula was retrieved
|
The sampling sites are up to 400 km apart. All water samples were taken on 19 April 2004, except for those from Lake Eckarfjärden, which was sampled on 16 April 2004. Two liters of water was sampled from a depth of
1m, depending on the depth of the water body. The water was stored at +4°C until the setup of the experiment (on 22 April). Physicochemical parameters (temperature, pH, conductivity, alkalinity, total phosphorus, total nitrogen, and chlorophyll a) were measured for each of the eight waters according to standard procedures. The DOC concentration was analyzed according to Langenheder et al. (39).
Experimental setup.
Bacterial cultures were set up in autoclavable 840-ml polycarbonate tubes. Each incubation vessel was closed without headspace with a polycarbonate piston with silicone sealings. The piston had two stainless steel pipes (one inlet and one outlet) with an inner diameter of 4 mm that were connected to Teflon tubing and sterile three-way valves. For sample collection, the outlet valve was opened and water was pressed out by moving the piston with the help of a hand lever.
The cultures were filled with autoclaved artificial lake water medium prepared according to the method of Bastviken et al. (3), using Nordic reservoir natural organic matter from the International Humic Substances Society (1R108N), diluted to a final DOC concentration of approximately 10 mg liter1. A stock solution with a DOC concentration of 520 mg liter1 was prepared from the freeze-dried powder and filtered twice through a sterile 0.2-µm Mini Filter polysulfone capsule (Gelman Sciences). Prior to addition to the cultures, the stock solution was filtered again through sterile 0.2-µm Supor syringe filters (Gelman Sciences). All filter units were extensively rinsed with autoclaved Milli-Q water prior to use to avoid contamination with external DOC. The pH in the cultures was 7.0 ± 0.03 (n = 3), and the alkalinity was 0.34 ± 0.02 meq liter1 (n = 3).
For preparation of the inocula, bacteria from each of the eight sites described earlier were filtered through GF/F filters (Whatman; precombusted for 8 h at 450°C). Since the concentration of bacteria was substantially lower in the groundwater than in the other samples, groundwater bacteria were concentrated. To achieve this, the water volume of the GF/F filtrate was decreased by filtration through autoclaved 0.2-µm-pore-size 47-mm polycarbonate filters, resulting in an approximately sixfold increase in bacterial abundance compared to the value given in Table 1. Ten milliliters of inoculum was added to each incubation vessel, and the resulting batch cultures were incubated at 19°C in the dark. Bacterial growth in the cultures was monitored, and the experiment was stopped when no further net increase in bacterial abundance was observed. Hence, the experiment was stopped between 242 and 505 h after inoculation, depending on the treatment. Three controls with only artificial lake water and natural organic matter extract were prepared and remained sterile throughout the entire experiment. Moreover, by use of microscopy, flagellates were confirmed to be virtually absent from the filtered water samples at the final stage of the experiment.
Generally, all glass- and plasticware coming into contact with the samples was soaked for several hours in 1 M HCl and rinsed with excessive amounts of Milli-Q water afterwards.
Oxygen consumption.
The net change in the dissolved oxygen concentration in the cultures was measured according to a spectrophotometric modification of the Winkler titrimetric method (48) in 30-ml Winkler flasks. The absorbance at 450 nm was measured and transformed into dissolved oxygen concentration using a calibration curve given by Mille-Lindblom and Tranvik (46). Bacterial respiration (BR) was defined as O2 utilization in the cultures and was transformed into carbon units by assuming a respiratory quotient of 1 (12).
Bacterial abundance, biomass, growth, and growth efficiency.
Subsamples of 10 ml were taken at 20- to 26-h intervals (up to 48-h intervals at later stages of the experiment), preserved and stained as described previously (15), and analyzed flow cytometrically using a FACScan flow cytometer (Becton Dickinson), with minor modifications of the original protocol (13). The cell volume and carbon content per cell were determined and calculated as described earlier (15). The bacterial biomass yield (YB) was obtained by multiplying the obtained cellular carbon content and the respective value of maximum cell abundance in the cultures. The bacterial growth efficiency (BGE) was calculated according to the following formula: BGE = [YB/(YB + BR)] x 100. The maximum intrinsic growth rate (µ) was calculated as the slope of the linear regression curve within the interval of linear increase of ln-transformed bacterial abundance plotted as a function of time, including at least three data points. The lag phase was defined as the period between inoculation and the observable increase in bacterial cell numbers in the cultures.
Activities of ectoenzymes.
The activities of two extracellular enzymes, leucine aminopeptidase and ß-glucosidase, were measured at the end of the incubation period (381 h for the SKOT and BOG inocula and 216 h for all other treatments). For leucine aminopeptidase activity, L-leucine 7-amido-4-methylcoumarin hydrochloride (Sigma) was dissolved in 5 mM bicarbonate buffer (pH 8), added to 200 µl of a sample from the cultures at a final concentration of 0.2 mM, and incubated at room temperature. At this concentration, the Vmax of the enzyme was measured (50). The release of the fluorescent product 7-amino-4-methylcoumarin (AMC) was measured once per hour for a total period of 10 h, using a MicroWell plate reader connected to a FluroMax-2 instrument (ISA Horiba Group) at a 380-nm excitation wavelength and a 440-nm emission wavelength. Enzyme activities were calculated from the increase in fluorescence for the time interval with a linear increase and converted into AMC concentrations with the help of a standard curve. ß-Glucosidase activity was measured following a similar procedure. 4-Methylumbelliferyl-ß-D-glucopyranoside (Sigma) was dissolved in methylcellosolve (high-performance liquid chromatography grade; Sigma) and added to 200-µl subsamples at a final concentration of 0.2 mM. The release of the fluorescent product 4-methylumbelliferyl was followed for 10 h at 365-nm excitation and 445-nm emission, and extracellular enzyme activities were calculated as described above. ß-Glucosidase activity decreased during the first 5 h (including the sterile controls) and started to increase linearly after this, perhaps because the 4-methylumbelliferyl was not completely dissolved when samples were added. ß-Glucosidase activities were calculated for the period of linear increase after 5 h of incubation, but the results should generally be treated with care.
[14C]benzoic acid uptake.
The utilization of radioactively labeled benzoic acid was measured at the end of the incubation period (381 h for the SKOT and BOG inocula and 216 h for all other treatments). Triplicate subsamples of 5 ml plus a blank where bacteria were killed by the addition of formaldehyde were incubated in test tubes sealed with rubber septum plugs at 20°C for 6 h with 14C-ring benzoic acid (8.0 mCi mmol1; Sigma) at a final concentration of 1.25 µM. The incubation was terminated by adding 250 µl of formaldehyde through the rubber septum with a syringe. 14CO2 was driven out of the samples by adding 250 µl of 4 M HCl and purging with air for 3 min, and the 14CO2 was collected in 2.5 ml of Carbo-Sorb S (Packard BioScience) in a scintillation vial, as described by Tranvik (52). After that, 2.5 ml of Permafluor E+ (Perkin-Elmer) was added, and radioactivity was counted in a Packard Tri-Carb 2100TR liquid scintillation counter. Particulate 14C was collected on 0.2-µm cellulose nitrate filters (Sartorius, Göttingen, Germany) that were presoaked in benzoic acid solution (1 g liter1), and the filters were subsequently rinsed five times with 2 ml of the same benzoic acid solution. Radioactivity was assessed by liquid scintillation counting after the addition of 3.5 ml Filtercount (Packard BioScience). Total uptake was calculated as the sum of both parameters.
Statistical analysis of functional parameters.
Differences in functional parameters depending on the inoculum source were tested using one-way multivariate analysis of variance (MANOVA; Pillai's trace test), and a subsequent Tukeys honestly significant difference (HSD) test was performed on results from univariate ANOVAs for each single functional parameter. Two different MANOVAs were performed: in the first analysis data from all treatments were included, whereas in the second one data from three treatments (LJUS, SKOT, and BOG) were excluded because they showed retarded growth and differed clearly from the rest. Hence, a new MANOVA was run including only the five treatments that were similar in their overall growth patterns. Principal component analysis (PCA) was performed on a correlation matrix of normalized data. Log transformation was generally used to normalize data, but values for BGE were arcsine square-root transformed. All statistical analyses were done using Statistica 6.0 (StatSoft, Inc.).
DNA extraction and terminal restriction fragment polymorphism length analysis.
At the end of the experiment, community compositions were analyzed by T-RFLP analysis as described previously (40). Briefly, DNAs from samples of 4.5 ml, collected by centrifugation, were extracted using a DNeasy kit (QIAGEN, Hilden, Germany), and PCRs were set up and performed as described earlier (40), with minor modifications. Three replicate PCR mixtures for each sample were pooled, treated with mung bean nuclease (14), and purified and concentrated using a QIAquick PCR purification kit (QIAGEN, Hilden, Germany). PCR products (
75 ng) were digested using the restriction enzymes HhaI, HaeIII, and AluI (Invitrogen, Carlsbad, Calif.), and hexachlorofluorescein-labeled fragments were separated and detected with an ABI 3700 96-capillary sequencer running in GeneScan mode (Applied Biosystems). T-RFLP electropherograms were inspected with the free software GenScan View 4 (CRIBI group [http://www.bmr.cribi.unipd.it]). The T-RFLP patterns obtained from the three different enzymes were pooled for each sample and further processed using a program in Visual Basic for Microsoft Excel (51). The data set for subsequent multivariate statistical analysis was constructed taking into account peaks with a size of 40 to 520 base pairs and a relative peak height of >1% of the total signal, and peaks that were fewer than 0.5 bases from a larger peak were merged. To account for small differences in running time among samples, peaks with <0.5-bp differences were considered to be of the same length. These settings were obtained by analyzing five replicates of the same sample, for which independent DNA extractions, PCRs, and restriction digests were made.
T-RFLP patterns were compared through calculation of a similarity matrix, using the Jaccard similarity coefficient, which calculates the similarity between a pair of samples based on binary data as follows: Sjxy = nxy/(nxy +
xy), in which nxy is the number of peaks common to both samples x and y and
xy is the number of peaks found only in x or in y. The matrix was calculated using the programs LecPCR and DistAFLP, which are part of the ADE-4 package (http://pbil.univ-lyon1.fr/ADE-4/microb/), and was analyzed by using the UPGMA (unweighted-pair group method using average linkages) clustering algorithm of Statistica 6.0 (StatSoft, Inc.).
|
|
|---|
![]() View larger version (16K): [in a new window] |
FIG. 1. UPGMA analysis based on T-RFLP patterns obtained from eight ambient sampling sites (A) and the final stages of the culture experiment (B). See Table 1 for assignments of treatment abbreviations to the different habitats from which the source communities were obtained. For the culture experiment, three or four replicates were included for each treatment.
|
![]() View larger version (23K): [in a new window] |
FIG. 2. Bacterial abundances in batch cultures inoculated with different starting communities. Note the different scaling of the x axis for SKOT and BOG curves. All symbols represent mean values ± standard errors (SE) (n = 4).
|
![]() View larger version (11K): [in a new window] |
FIG. 3. Length of lag phase (adaptation time) is dependent on the pH at the sampling sites from which the inocula were retrieved. All symbols represent mean values ± standard deviations (n = 4).
|
![]() View larger version (15K): [in a new window] |
FIG. 4. PCA with functional parameters. Parameters included in the analysis were bacterial abundance (abund), biomass yield (YB), bacterial respiration (BR), bacterial growth efficiency (BGE), length of lag phase (lag), maximum intrinsic growth rate (µ), beta-glucosidase activity (beta), amino-leucine peptidase activity (amino), and benzoic acid uptake (benz). All values were log(x + 1) transformed, except for µ (untransformed) and BGE, which was arcsine square-root transformed. (a) Distribution of treatments (mean ± SE; n = 4) along the first two principal components. (b) Loadings of parameters included in the analysis.
|
|
View this table: [in a new window] |
TABLE 2. Results from one-way ANOVA to test the effect of the origin of the bacterial inoculum on functional parametersa
|
![]() View larger version (34K): [in a new window] |
FIG. 5. Functional parameters measured during or at the end of the experiments. Only results for ROCK, GROW, ECKA, LUMP, and VATT are shown. YB, bacterial biomass yield; BR, bacterial respiration; BGE, bacterial growth efficiency; µ, maximum intrinsic growth rate. All bars represent mean values ± SE calculated from four replicate cultures. A post hoc test revealed no significant differences between treatments for YB, BR, and BGE. For VATT, µ was significantly lower than that for LUMP.
|
|
View this table: [in a new window] |
TABLE 3. Activities of ß-glucosidase, leucine-aminopeptidase, and 14C-ring benzoic acid uptake in ROCK, GROW, ECKA, LUMP, and VATTa
|
|
|
|---|
Previous experiments provided convincing support of the "everything is everywherethe environment selects" concept for eukaryotic microbes (17). These studies included a range of enrichment strategies and incubations extended up to roughly 4 months in order to provoke the growth of very rare cryptic species. In contrast, our experiments lasted 242 to 505 h. Based on the assumption that rare or cryptic populations have the potential to become dominant under the right circumstances, a longer incubation time would possibly also reflect a higher degree of ubiquity for the bacteria in our experiment, if ubiquitous species in our inocula were rare and needed extensive enrichment to reach detectable quantities. It is possible to gauge the hypothetical time (t) it takes for such rare populations to reach 1% of the total bacterial concentration (assumed to be 109 cells liter1), i.e., a fraction of the total community that is likely to be easily detected by T-RFLP. This is done by resolving the exponential growth equation for t, i.e., Nt = Ni x eµt, where µ is the exponential growth rate, Ni is the initial abundance, and Nt is the number of cells at which the population becomes detectable. Thus, for a rare population with an initial abundance of 100 cells liter1 in a culture growing with a doubling time of 11 h (corresponding to the average community growth rate found in our experiment), it would take roughly 180 h for the population to reach a detectable size of 107 cells liter1. The actual incubation times in the experiment were 242 to 505 h, suggesting that rare populations would in fact be able to reach detectable cell densities unless they were growing substantially slower than the average community growth rate. In general terms, the calculations above illustrate that the successful proliferation of rare or cryptic species enabling them to reach high levels of abundance requires substantial time under suitable conditions. Hence, their window of opportunity may be restricted to environments with stable environmental conditions at scales that are relevant for bacteria. Hence, in frequently disturbed communities such as pelagic environments, time might simply not be sufficient for cryptic species to become dominant members of the community. Thus, the composition of a local bacterial community should not exclusively be regulated either by selection by the local environmental condition or by "metacommunity processes." In contrast, the two mechanisms are not mutually exclusive and occur simultaneously but apply to different time scales, and their relative importance might vary between habitats differing with regard to their disturbance regimen (5, 41).
At a rough glance, the growth of differently composed communities under identical conditions seems to contradict the previously suggested global distribution of freshwater bacteria (25, 29, 30, 55, 56). However, due to the lack of information about the identity of community members enriched in the culture experiments, our results cannot be used to deduce that freshwater planktonic bacteria are not globally dispersed. Even though communities seemed to be different from each other depending on the origin of the source community, they might still contain highly related, though different, representatives belonging to the same bacterial cluster. Accordingly, studies of sediments and soils suggest that the degree of ubiquity of microbial taxa depends on the level of phylogenetic resolution (6, 33). Additionally, it has also been shown that bacteria with identical 16S rRNA gene sequences (which by operational definition would belong to the same species) can have clearly different ecophysiological and phenotypic traits (30, 35). Hence, even if there is a global dispersal of bacterial species, the local adaptation of genotypes within these species might lead to various responses by the same species (as defined by 16S rRNA sequences) to changes in environmental conditions.
There was some resemblance between the clustering patterns between culture communities and ambient sites. The most apparent was that the GROW (groundwater) and BOG (bog pool) samples had the highest fractions of unique peaks in the original samples as well as in the cultures. Groundwater might be unique because it lacks direct contact with airborne bacteria, and its connectivity to other habitats is very limited. The bog pool provides hostile conditions (pH < 4) for many bacteria and might therefore also harbor a highly adapted community. Hence, GROW and BOG cultures developed clearly distinct communities because they were inoculated with clearly distinct source communities. More unexpected was the observation that source communities from two lakes (LJUS and SKOT) also resulted in the growth of quite distinct bacterial communities (Fig. 1B). It appeared, however, that pH was generally important in regulating the overall growth patterns of selected bacterial communities (Fig. 2). Hence, all bacterial source communities from acidic sites (LJUS, SKOT, and BOG) had extended lag phases of several days (in the case of LJUS and SKOT; Fig. 2) or even weeks (BOG) before growth could be observed. Possibly, a long apparent lag phase is due to a "shock response" of bacteria exposed to strong deviations in pH, resulting in either a very long adaptation time or a very small effective inoculum of specific bacteria that are able to grow under the culture conditions.
Although the variability among replicate cultures with the same inoculum was always less than the variability between cultures with different inocula (Fig. 1B), there was a high fraction (30% or more) of peaks that were unique to each replicate culture. Recently, Hewson and Fuhrman (31) could show for another fingerprinting method (automated ribosomal intergenic spacer analysis) that it is highly reproducible as long as only the largest peaks, i.e., those comprising >0.5% of the total signal, are taken into account. Assuming that the same accounts for T-RFLP analysis, it seems likely that the observed among-replicate variability in BCC resembles real differences in taxon composition and is not a methodological artifact. The size of the source community (i.e., the total number of individuals), its diversity, and the distribution of individuals among taxa are crucial for community assemblage patterns (11). Highly similar BCC among replicate cultures requires that replicate inocula (i.e., the replicate 10-ml subsamples of the source community) also show a high degree of similarity. Deviating patterns in BCC among replicate cultures might thus simply reflect that the inocula were not identical. Long and Azam (43) demonstrated that variations in BCC are minor from the microliter to the liter scale, but this accounted only for the most abundant community members. Considering that some of the taxa that grew in the cultures may have been rare in the source community, it seems likely that differences in the presence or absence of rare species and the absolute abundance of taxa caused the observed variability in community assembly patterns among replicate cultures. Together with size-selective removal of cells during filtration and cultivation-induced shifts in community structure and diversity (16, 44), the random exclusion of rare taxa from the inocula also might have limited the potential of the different source communities to converge in BCC upon growth under identical conditions.
Coupling between community structure and function.
The second major aim of this study was to test if there is a link between BCC and the functioning of the system. The pH-affected communities (LJUS, SKOT, and BOG) were functionally clearly different from those within the remaining treatment groups, i.e., ROCK, GROW, ECKA, LUMP, and VATT (Fig. 2 and 4). Similar effects of pH on bacterial community composition and functioning have been observed previously (40, 47) and agree with findings from field studies suggesting that pH is a strong regulating factor of BCC in lakes (42, 54) and a major threat to the survival of allochthonous bacteria in aquatic systems (2). On the other hand, broad-scale functions such as bacterial biomass production and respiration were similar in treatments ROCK, GROW, ECKA, LUMP, and VATT, independently of BCC, suggesting functional redundancy. This agrees with the results of a previous full-factorial switch experiment (40), where we investigated the importance of the source of the growth medium versus the inoculum as shaping factors of community composition and broad-scale functioning and found a loose connection between community composition and function. Other culture studies also found that functioning was rather independent of the origins of inocula (21, 22, 53).
On the other hand, there are studies showing that community structure matters for function (4, 37, 38). It is interesting to observe, however, that these studies mostly measured "narrow" functions, such as pesticide-induced or semilabile DOC degradation and enzyme activities, while studies supporting that function is independent of BCC typically address "broad" functions, such as respiration and biomass yield (e.g., see reference 40). In this study, different bacterial communities showed at least partly differing ectoenzyme activities (Table 3) and various potentials to utilize benzoic acid (Table 3), which was used as a model for aromatic ring structures. This indicates that different components within the total DOM pool were degraded in cultures with different inocula. Thus, specific functions such as enzyme activities differed in treatments with different BCC, whereas there were no observable effects on broader functions such as biomass production and community respiration. This idea is compatible with the finding of functional equivalency of different fungal species with respect to carbon mineralization (49) and that well-defined narrow niche functions are more sensitive to a reduction in diversity (27) or perturbation (23) than broad-scale functions carried out by a wide range of organisms.
To summarize, differently composed bacterial communities developed from different source communities under identical conditions. These structural differences, however, were only partly reflected in differences in community function, and the functions that appeared sensitive to community composition were relatively specific, while aggregated functions (e.g., community biomass and respiration) were insensitive to community composition. Very similar findings have been made for zooplankton communities (36), suggesting that the overall mechanisms behind community structure and function might be similar for microbes and higher organisms.
This work was supported by grants from the Swedish Research Council to L.J.T. and the Helge Ax:son Johnson and Malméns foundations to S.L.
|
|
|---|
| |||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||
Copyright © 2009 by the American Society for Microbiology. For an alternate route to Journals.ASM.org, visit: http://intl-journals.asm.org | More Info»