Previous Article | Next Article ![]()
Applied and Environmental Microbiology, May 2005, p. 2260-2266, Vol. 71, No. 5
0099-2240/05/$08.00+0 doi:10.1128/AEM.71.5.2260-2266.2005
Copyright © 2005, American Society for Microbiology. All Rights Reserved.
Department of Biological Oceanography, Royal Netherlands Institute for Sea Research (NIOZ), P.O. Box 59, 1790 AB Den Burg, Texel, The Netherlands
Received 17 August 2004/ Accepted 22 November 2004
|
|
|---|
43% of the variation in BGE. Our results indicate that despite the observed shifts in the community composition, the main function of the bacterioplankton, the remineralization of dissolved organic carbon to CO2, is rather stable. |
|
|---|
The extent to which variations in bacterial growth yield, and particularly in respiration, are related to shifts in the bacterioplankton community composition has not been studied in detail yet. Thus, the link between the phylogenetic composition of the bacterioplankton community and its major function in carbon cycling, i.e., the remineralization of DOC, is unclear, although shifts in the bacterioplankton community composition on a seasonal scale have been reported and related to the occurrence of specific phytoplankton blooms (1). Consequently, microbial communities may respond to changing environmental conditions with physiological adaptation or by shifts in the community composition (2, 11, 37, 54).
Naeem and Li (35) and Morin and McGrady-Steed (34) presented evidence that species loss changes important aspects of ecosystem functioning. In mesocosm experiments, increased biodiversity decreased the variation in the measured CO2 flux (30). However, experiments linking biodiversity with ecosystem functioning are usually performed with the aim to maximize experimental control and therefore are likely to miss the spatial and temporal heterogeneity in natural environments (reference 5 and references therein).
For this study, we tested the hypothesis that bacterial production and respiration are largely independent of the bacterioplankton community composition. For this purpose, we conducted six cruises in the southern North Sea, covering a full seasonal cycle. We measured bacterial production (BP), bacterial respiration (BR), and the bacterial community composition by terminal restriction fragment length polymorphism (T-RFLP) and compared the patterns of bacterioplankton richness with these bacterial activity measurements.
|
|
|---|
![]() View larger version (24K): [in a new window] |
FIG. 1. Map of study area. Dots indicate the individual stations occupied during the six cruises. The months of sampling are indicated with different symbols.
|
Bacterial production.
Bacterial production in the 0.8-µm-pore-size filtered seawater was measured by [14C]leucine incorporation (specific activity, 0.295 Ci mmol1; final concentration, 10 nmol liter1). Two samples and one blank were incubated in the dark. The blank was fixed immediately with concentrated formaldehyde (final concentration, 4% [vol/vol]) 10 min prior to adding the tracer. After incubation of the samples and the blank at the in situ temperature for 60 min, the samples were fixed with formaldehyde (4% final concentration), filtered with 0.45-µm nitrocellulose filters (25-mm diameter; Millipore HA), and rinsed twice with 5 ml ice-cold 5% trichloroacetic acid (Sigma Chemicals) for 5 min. The filters were dissolved in 1 ml ethyl acetate, and after 10 min, 8 ml of scintillation cocktail (Insta-Gel Plus; Canberra Packard) was added. The radioactivity incorporated into bacterial cells was counted in a liquid scintillation counter (model 1212; LKB Wallac). The amount of leucine incorporated into the bacterial biomass was converted to carbon production by using the empirical conversion factor 0.07 x 1018 cells mol1 Leu (42) and assuming a C content of Bacteria of 20 fg C cell1 (25). Applying this conversion factor resulted in bacterial production estimates similar to the theoretical factor of 1.55 kg C mol1 Leu, assuming there was no isotope dilution (48; data not shown). In the text below, the abbreviation BP is used for bacterioplankton production measured in filtered (0.8 µm) seawater.
Bacterial respiration.
The filtrate (filtered through a 0.8-µm filter) was carefully transferred to calibrated borosilicate glass BOD bottles with a nominal volume of 120 ml by a sipper system to avoid the introduction of air bubbles. For determinations of the initial O2 concentration (T0), samples were fixed immediately with Winkler reagents and incubated together with the live samples in a water bath in the dark at the in situ temperature (±1°C) for 12 to 24 h, after which the incubations were terminated (T1). Triplicate bottles were used for determinations of the initial and final O2 concentrations. All glassware was washed with 10% HCl and thoroughly rinsed with Milli-Q water prior to use. The oxygen concentrations in the T0 and T1 sample bottles were measured spectrophotometrically in one run (36, 45), essentially according to the standard protocol for the determination of oxygen by Winkler titration (7). The amount of total iodine was determined at a wavelength of 456 nm. Measurements were done at 20°C on a Hitachi U-1100 spectrophotometer using a 1-cm flowthrough cuvette. To increase the sensitivity of the absorbance readings, we connected a four-digit voltmeter (Metex M4650) to the spectrophotometer. Calibration was performed by standard additions of iodate to distilled water, resulting in an empirical coefficient of 0.54455 nM cm1 (G. Kraay, personal communication). The samples were withdrawn from the BOD bottles with a Teflon tube and a peristaltic pump (Gilson Minipuls) and directly fed to the flowthrough cuvette of the spectrophotometer. The end of the tube was placed near the bottom of the bottles to avoid a possible loss of volatile iodine. The spectrophotometer was zeroed against Milli-Q water. The coefficient of variation of the oxygen determinations was <0.5%.
T-RFLP sampling and analysis.
A total of 19 liters of seawater was filtered through 0.8-µm-pore-size filters (Isopore ATTP, 142-mm diameter; Millipore), using a stainless steel filter holder and an air pressure pump, to exclude most of the nonbacterioplankton particles from the analysis. The bacterioplankton fraction (<0.8 µm) was concentrated to a final volume of 400 to 550 ml by tangential-flow filtration with a 0.22-µm-pore-size filter cassette (Pellicon PTGVPPC05; Millipore). Subsequently, the bacterial concentrate was filtered onto 0.22-µm-pore-size filters (Isopore GVWP, 100-mm diameter; Millipore) and the filters were frozen in liquid nitrogen and subsequently stored at 80°C until further processing in the lab. The filtration units were thoroughly rinsed with sample water prior to use and soaked in 1 N HCl between samplings.
Extraction of the nucleic acids from the filters was performed as previously described. Briefly, four freeze-thaw cycles (196°C to 37°C) and a subsequent treatment with lysozyme (Sigma) and proteinase K (Fluka) in 1% sodium dodecyl sulfate were conducted. The liquid phase was extracted with a mixture of phenol, chloroform, and isoamyl alcohol and the nucleic acids were precipitated with ethanol overnight at 20°C (52). The resulting pellet was redissolved in 100 µl ultrapure water (Sigma). The nucleic acids in 50 µl of this solution were further purified by use of a QIAEX II gel extraction kit (QIAGEN) as recommended by the manufacturer for DNA fragments larger than 10 kbp. The nucleic acids were recovered in a final volume of 20 µl of elution buffer (QIAGEN) and used for subsequent PCR amplification. The integrity of the DNA was checked by agarose gel electrophoresis.
PCR conditions and chemicals were applied as described by Moeseneder et al. (33). Briefly, 1 to 2 µl of the cleaned nucleic acid extract was used as a template in a 50-µl PCR mixture. The Bacteria-specific primer 27F and the universal primer 1492R (24) were used to amplify a ca. 1,480-bp fragment of the bacterial 16S rRNA gene. The forward primer was fluorescently 5' end labeled with the phosphoramidite fluorochrome 5-carboxy-fluorescein, and the reverse primer was labeled with 6-carboxy-4',5'-dichloro-2',7'-dimethoxyfluorescein (both from Interactiva, Germany). After PCR amplification, excess fluorescently labeled primers were removed by ethanol precipitation and subsequent gel purification in a 1% agarose gel using 6-carboxytetramethylrhodamine loading dye (Perkin Elmer-Applied Biosystems). The PCR fragments were recovered from the gel with a QIAquick gel extraction kit (QIAGEN) in a final volume of 20 µl of elution buffer (QIAGEN). The concentrations of the PCR fragments were estimated with a 1% agarose gel containing a standard (Smart-Ladder; Eurogentec).
Restriction digestion of 50 ng of the purified PCR fragments was performed with 20 U of the restriction enzyme HhaI together with the recommended buffer (both from Amersham Pharmacia) in a total volume of 100 µl. Restriction was performed at 37°C for 12 h to ensure complete digestion. The DNA fragments from the restriction digests were recovered in a final volume of 2 µl ultrapure water (Sigma) by linear polyacrylamide precipitation (32). T-RFLP analysis was performed on an automated capillary sequencer (ABI Prism 310; Perkin Elmer-Applied Biosystems) as previously described (33).
Statistical analyses.
The T-RFLP patterns were analyzed by recording the number of peaks (presence versus absence). Due to the variability in the discriminative power of the forward and the reverse primers, subsequent analysis was performed on the combined data set from both primers, serving as a relative measure of bacterial community richness (32).
The similarities of the T-RFLP patterns between the different stations were assessed with the Bray-Curtis similarity index, which is similar to Sorenson's similarity index when applied to presence/absence data (8). The resulting similarity matrix was assessed by nonparametric multidimensional scaling (MDS), which is a powerful tool for assessing community profiles obtained by molecular fingerprinting techniques (26, 46, 53). For MDS analysis, the recommendations of Clarke and Warwick (8) were followed. With MDS, the complexity of a given similarity matrix is reduced by plotting the data two-dimensionally. The more similar the samples are, the closer they appear in the plot. The stress factor is an indication of the goodness of fit of the raw versus plotted data, with values of <0.1 representing a good ordination of the similarity matrix with little risk of misinterpretation. Higher stress factors were checked by cluster analysis.
Dendrograms of the obtained T-RFLP patterns were constructed by using either the Bray-Curtis similarity matrix or the bootstrapped Jaccard or Simple Match similarity matrix calculated from the same data set. Clusters were constructed with the unweighted-pair group method using average linkages. MDS calculations were done with the software package Primer 5 from Primer-E, and other statistical analyses were done with the software package Statistica from Statsoft.
DOC measurement.
Samples for DOC measurements were filtered through rinsed 0.2-µm polycarbonate filters and sealed in combusted (450°C for 4 h) glass ampoules after the addition of 50 µl of 40% phosphoric acid. Subsequently, the samples were stored frozen at 20°C. DOC concentrations were determined by the high-temperature combustion method using a Shimadzu TOC-5000 analyzer (4). Standards were prepared with potassium hydrogen phthalate (Nacalai Tesque, Inc., Kyoto, Japan). Ultrapure Milli-Q blanks were run before and after the sample analysis. The blank value was, on average, 16.3 ± 6.8 µM, and the mean of triplicate injections was calculated for each sample. The average analytical precision of the instrument was <3%.
Chl a measurement.
One-liter samples were gently filtered through 47-mm Whatman GF/F filters and stored at 60°C until analysis within 4 weeks. Chlorophyll a (Chl a) was extracted in 10 ml of 90% acetone at 40°C in the dark for 48 h. Subsequently, the filters were sonicated on ice for 1 min (Branson model 3200) and centrifuged to remove particles. The Chl a concentration in the supernatant was determined fluorometrically with a Hitachi F-2000 fluorometer (21).
|
|
|---|
40% ± 20% compared to unfiltered seawater samples.
Seasonal dynamics in physicochemical parameters.
Low salinity values were recorded in April but varied only within a rather narrow range during the rest of the year. The temperature changed significantly between months (least significant difference test; P < 0.05), with the lowest temperatures in spring and winter and an average maximum temperature in September. Monthly averages for DOC concentrations showed no clear pattern and varied around an annual average of 137 ± 69 mmol C m3. Chl a concentrations peaked in April and September and were low during the other months (Table 1).
|
View this table: [in a new window] |
TABLE 1. Monthly means and SD for selected physicochemical and biological parameters of the southern North Sea
|
1.0 x 106 ± 0.5 x 106 ml1. BP decreased significantly from the late spring towards the winter (least significant difference test; P < 0.0001). The highest BP level was measured in April, with 1.7 mmol C m3 day1, while in December, BP was only 0.03 mmol C m3 day1. BR was variable over the seasonal cycle, with an average of 1.4 ± 1.0 mmol C m3 day1. Bacterial growth efficiencies (BGEs) ranged from 42% in April to 5% in December, with an overall average of
21% ± 14% (Table 1). Monthly averages in BGEs were variable but decreased significantly overall from spring to winter (analysis of variance [ANOVA]; P < 0.05).
Spatial versus seasonal heterogeneity in the richness of the bacterial community.
As indicated by the dendrogram based on the Bray-Curtis similarity index (Fig. 2), bacterial communities were more closely related within individual months than between different months. An MDS plot highlights the seasonal heterogeneity of the bacterial community, but spatial dissimilarities are also visible (Fig. 3). Major changes in the bacterioplankton community composition from April to June were apparent, while between July and December these changes were much smaller (Fig. 3).
![]() View larger version (24K): [in a new window] |
FIG. 2. Cluster analysis of T-RFLP patterns from an annual cycle in the southern North Sea by the unweighted-pair group method using average linkages. The distance matrix was calculated by using the Bray-Curtis similarity index. The main groups are shaded and correspond to the months of sampling, and the branches represent stations. Subgroups are indicated with dotted lines for relative similarities of 75 and 83%. Numbers next to the lines denote the groups used for the analysis depicted in Fig. 7.
|
![]() View larger version (27K): [in a new window] |
FIG. 3. Two-dimensional MDS plot of T-RFLP patterns of the bacterioplankton community of the southern North Sea. The configuration found by MDS was based on a Bray-Curtis similarity matrix for the bacterial richness between the stations. The months are indicated by different symbols. The arrows inside the panel indicate the seasonal sequence from spring to winter. The arrows on dimensions 1 and 2 indicate the directions in which scores increase. Spearman rank correlations were calculated between the scores of the dimension found for the different stations and bacterial richness, bacterial production (BP), bacterial respiration (BR), chlorophyll a (Chl a), dissolved organic carbon (DOC), and temperature (T). SF, stress factor; *, P < 0.01.
|
|
View this table: [in a new window] |
TABLE 2. Spearman rank correlations of biological and selected physicochemical parameters
|
52 operational taxonomic units), cell-specific BP leveled off (Fig. 4a). The decline in cell-specific BR was less pronounced and more variable (Fig. 4b).
![]() View larger version (23K): [in a new window] |
FIG. 4. Relationship between bacterioplankton richness and cell-specific bacterial production (BP; fmol C cell1 day1) (a) or cell-specific respiration (BR; fmol C cell1 day1) (b). Exponential decreases are illustrated by a power regression model fit to the raw data.
|
80 operational taxonomic units (one-way ANOVA; r2 = 0.68; P < 0.0001), while BGE significantly decreased from April to December, from
42% in the spring to 5% in the winter (one-way ANOVA; r2 = 0.45; P < 0.0001) (Fig. 5a). As a consequence, BGE was negatively related to bacterioplankton richness, which explained
43% of the variation in BGE (Fig. 5b).
![]() View larger version (25K): [in a new window] |
FIG. 5. Dynamics of bacterial growth efficiency (BGE) and richness during the seasonal cycle in the southern North Sea. (a) Monthly averages of BGE and bacterioplankton richness measured by T-RFLP. Error bars indicate standard deviations of the means (n = 8 to 19). (b) Relationship between BGE and bacterioplankton richness, with months indicated by different symbols.
|
|
|
|---|
Evidence that biodiversity can influence rates of ecosystem processes has come from experimental studies (20, 27, 30, 35). Usually, the opposing view that the ecosystem function affects species diversity is reported. However, relationships between diversity and productivity are, in fact, two-dimensional projections of a three-dimensional relationship among factors, e.g., site productivity, metabolic rates, and species diversity (18, 47). Thus, conclusions based on effects of changes in bacterial richness on ecosystem properties seem appropriate.
There is much debate about the relationship of species diversity to productivity (27, 47). Chl a is commonly used as a simple indicator of phytoplankton productivity, and some studies have interpreted changes in the bacterial community composition as a consequence of gradients in Chl a concentrations (23, 38). However, analyzing the patterns of
- and ß-proteobacteria and Bacteroidetes richness individually, Horner-Devine et al. (22) found that the richness within the individual groups changed differently along a Chl a gradient but that the overall bacterial richness showed no trend. Also, in our study, only a weak correlation was found between Chl a and bacterioplankton richness (Table 2). Due to the effects of statistical averaging, diverse communities may appear more stable than they actually are, as pointed out by Doak (14) and Tilman (50). This effect might mask productivity versus richness patterns in whole-community analyses.
With seawater cultures, Covert and Moran (10) and Carlson et al. (6) found that shifts in the bacterial community were related to the quality of the DOC. In the southern North Sea, DOC concentrations showed no clear seasonality (Table 1) and, consequently, no substantial correlation with bacterial richness (Table 2). However, a negative correlation in the MDS plot was apparent (Fig. 3). Bulk DOC measurements are usually not a good indicator of the direct availability of a substrate to bacteria since the DOC pool consists of >80% semilabile to refractory compounds (3). Nevertheless, the weak gradient detected suggests that changes in the DOC concentration and/or composition play a role in structuring the composition of bacterial communities in the southern North Sea. This was also evident in a previous study performed at a coastal site in the North Sea for which the wax and wane of a phytoplankton bloom were monitored (1).
Generally, low BP is measured in the winter. In the early summer, BP and BR are high (41) but community richness is low in the southern North Sea (Fig. 5a). From published reports, Mittelbach et al. (31) deduced that both hump-shaped and positive relationships of species diversity to productivity are common in nature. In many studies, however, primary production is compared to the diversity of animals, i.e., complex feedback interactions between productivity and diversity are to be expected. Thus, direct productivity estimates of the studied organisms should be used (19). The bacterial carbon demand (BCD), calculated as the sum of BP and BR, reveals the carbon flux mediated by the bacterial community. We found a broad negative relationship between bacterial richness and the BCD (Fig. 6). This would argue for a highly active bacterial assemblage in spring with few community members and for reduced productivity with increasing bacterial richness.
![]() View larger version (16K): [in a new window] |
FIG. 6. Relationship between bacterioplankton richness and bacterial carbon demand (BCD; mmol C m3 day1).
|
Generally, stability-diversity hypotheses argue that with increasing diversity or species richness, community responses tend to become more stable (16, 29). We assessed the stability-diversity relationship for North Sea bacterioplankton, using bacterioplankton richness and activity, by constructing a dendrogram based on the similarity matrix and then intersecting the resulting tree at relative similarities of 75 and 83% (Fig. 2) in a way that allowed calculations of averages and standard deviations. The tendency of variance to increase with the mean is a major source of bias in variability estimates (9). McArdle et al. (28) argued that the standard deviation of log-transformed raw data [standard deviation (SD) of log(Xi)] is less sensitive to highly skewed distributions, as is often the case for aggregate community properties, than the coefficient of variation.
When the standard deviations of log-transformed cell-specific BP obtained from the individual months were plotted against the respective averages of bacterioplankton richness, a significant increase in variability with increasing richness became apparent (r2 = 0.70; P < 0.002) (Fig. 7a). This suggests that increasing richness leads to an increased instability of cell-specific BP. However, the variability of cell-specific BR barely increased along the richness gradient (Fig. 7b). This is in contrast to the results of McGrady-Steed et al. (30), who measured the CO2 flux in mesocosms supplied with increasing species richness of a manipulated aquatic microbial assemblage. They found that with increasing biodiversity, ecosystem respiration becomes more predictable, i.e., the standard deviation for replicates with similar species richness decreased with increasing richness. Our data on natural bacterioplankton assemblages, however, suggest that bacterioplankton respiration is rather independent of shifts in bacterioplankton community richness. Thus, the main function of heterotrophic bacterioplankton, i.e., the remineralization of DOC to CO2, is rather stable over the seasonal cycle.
![]() View larger version (15K): [in a new window] |
FIG. 7. Relationship between standard deviations (SD) of log-transformed raw data [SD of log(Xi)] and averages of bacterioplankton richness calculated for the subgroups with >83% relative similarity, as indicated in the dendrogram in Fig. 2. The graphs show the relationship between richness and the SD of log BP/cell (a) or the SD of log BR/cell (b).
|
This work was supported by the Dutch Science Foundation (NWO-ALW) project 811.33.002 and the European Commission through the COMET and AIRWIN projects to G.J.H. This work was done in partial fulfillment of the requirements for a Ph.D. degree from the University of Groningen by T.R.
|
|
|---|
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»