Previous Article |
Table of Contents
| Next Article 
BIOLOGICAL SCIENCES / ECOLOGY
Diversity predicts stability and resource use efficiency in natural phytoplankton communities
*Norwegian Institute for Water Research, Gaustadalléen 21, 0349 Oslo, Norway; European Commission, Joint Research Centre, Institute for Environment and Sustainability, I-21020 Ispra, Italy; Department of Biology, University of Oslo, P.O. Box 1066, Blindern, 0316 Oslo, Norway; ¶Finnish Environment Institute, P.O. Box 140, FIN-00251, Helsinki, Finland; and ||Swedish University of Agricultural Sciences, P.O. Box 7070, SE-750 07 Uppsala, Sweden
Edited by Paul G. Falkowski, Rutgers, The State University of New Jersey, New Brunswick, NJ, and approved February 5, 2008 (received for review September 3, 2007)
 |
Abstract |
|---|
 Top Abstract ResultsDiscussionMaterials and MethodsAcknowledgments.References |
|---|
50% to the global primary production. We show here that phytoplankton resource use, and thus carbon fixation, is directly linked to the diversity of phytoplankton communities. Datasets from freshwater and brackish habitats show that diversity is the best predictor for resource use efficiency of phytoplankton communities across considerable environmental gradients. Furthermore, we show that the diversity requirement for stable ecosystem functioning scales with the nutrient level (total phosphorus), as evidenced by the opposing effects of diversity (negative) and resource level (positive) on the variability of both resource use and community composition. Our analyses of large-scale observational data are consistent with experimental and model studies demonstrating causal effects of microbial diversity on functional properties at the system level. Our findings point at potential linkages between eutrophication and pollution-mediated loss of phytoplankton diversity. Factors reducing phytoplankton diversity may have direct detrimental effects on the amount and predictability of aquatic primary production.
biodiversity | carbon cycle | ecosystem functioning
In view of globally-increasing species losses, the importance of diversity for ecosystem processes, such as resource use and productivity, have recently received considerable attention. Studies on terrestrial plants were the first to show that levels of primary production generally increase with increasing diversity and that more diverse communities are more resistant to extreme events (4). Meanwhile, there is accumulating evidence that this positive effect on productivity exists across various habitats and trophic levels (5–7). Cardinale et al. (8) and Gross and Cardinale (9) have shown that this positive effect may largely be attributed to more efficient resource use in more diverse communities.
Phototrophic bacteria and protists in the surface layers of lakes and oceans (phytoplankton) account for 50% of the global primary production, making this polyphyletic group a major component of the global carbon cycle (10). A few milliliters of sea water may contain tens to hundreds of species from very different taxonomic groups. Despite their obvious importance, however, our knowledge about the functional role of phytoplankton diversity (and microbial diversity in general) is very limited (11, 12). In addition to conceptual problems regarding the diversity of unicellular organisms (such as the species concept in organisms with largely asexual reproduction), it is generally unknown whether microbial diversity relates to ecosystem functioning in ways similar to those demonstrated for higher organisms.
The enormous dispersal potential of microbes and many species' being found across vast spatial scales (13, 14) have been the basis for the assumption that the spatial distribution of unicellular organisms is not dispersal limited, which should imply that the number of species present in a phytoplankton community becomes only a matter of local environmental conditions. With respect to ecosystem functioning, it has been proposed that the "local species richness [of microbes] is always sufficient to drive ecosystem functions" (13).
An increasing number of studies question the view of a global distribution of microbes (15, 16). Microbes show biogeographical patterns that do not always correlate to environmental gradients (e.g., refs. 17 and 18). For example, the number of diatom species found in a given lake correlates with the number of surrounding lakes (19). The existence of spatial diversity patterns alone cannot prove, however, that ecosystem function is related to microbial diversity.
Experiments with artificial communities of phytoplankton and other microbes indicate that processes carried out by microbes might be related to diversity (20–22) and that the relationships are similar to those known from the macroscopic realm. Such artificial communities, however, usually consist of arbitrary assemblages taken from easily cultivatable taxa and contain far fewer species than those found in natural assemblages, raising a serious concern as to what extent such "synthetic" communities exhibit natural responses.
Here, we investigate whether the variation in natural phytoplankton diversity spans a range that affects ecosystem functions and whether such patterns are comparable to those found for higher organisms. In particular, we test (i) whether phytoplankton diversity enhances resource use efficiency and (ii) whether diversity dampens variability in resource use and community composition. We use >3,000 phytoplankton samples from Scandinavian lakes and from the Baltic Sea. [See supporting information (SI) for a map with the sampling stations.]
Phosphorus is generally considered the limiting resource for phytoplankton primary production in temperate lakes (23, 24) and frequently represents the (co-) limiting nutrient in the Northern Baltic Sea (25). Moreover, phosphorus is rarely available in excess even in nitrogen limited systems (26, 27). Therefore, we use total phosphorus (TP) as a proxy for potential system productivity and express phytoplankton resource use efficiency (RUE) by the ratio between phytoplankton biomass and TP. To increase robustness of results, we use two independent measures of algal biomass, namely chlorophyll-a and algal carbon content as derived from microscopical cell counts, yielding two measures of resource use efficiency (RUEchl and RUEcarb). As a proxy for phytoplankton diversity, we use genus richness (G). Because both diversity and system productivity may affect system stability (28), we include TP as a covariable into our analysis.
|
Results |
|---|
|
TopAbstractResultsDiscussionMaterials and MethodsAcknowledgments.References |
|---|
|
|
= 0.16, P < 0.001). RUEchl is consistently high along the upper flank of this relationship but variable where high TP is met by low G (Fig. 2A).
|
Variability of RUE and Community Composition. Temporal and spatial variability of RUE were analyzed separately. In a heteroscedasticity analysis (see Materials and Methods), we tested whether G and TP predict variability of RUE among sites (i.e., spatial variability). In the combined dataset and in most of the separate datasets, error variances showed significant negative relationships with G (Table 2). Moreover, although TP did not have a consistent effect on absolute RUE (Table 1), the heteroscedasticity analysis revealed that variability of RUE increased with TP in the Norwegian, Swedish, and combined datasets (Table 2).
|
|
|
|
Discussion |
|---|
|
TopAbstractResultsDiscussionMaterials and MethodsAcknowledgments.References |
|---|
The patterns appear to be consistent among the different datasets. Moreover, despite considerable differences in regional diversity and environmental gradients (brackish vs. freshwater and continental vs. oceanic climate), diversity was the best predictor for RUE in all datasets (see Materials and Methods and SI). This indicates a general and strong dependency between phytoplankton diversity and resource use. The positive relationship between diversity and RUE is much more pronounced in areas with a low average diversity (Norwegian lakes and the Baltic Sea) compared with those with higher diversity (Finnish and Swedish lakes), indicating that diversity effects are stronger for less diverse communities.
Although species richness is often analyzed as a function of productivity (e.g., 29, 30), the direction of causality has been seriously challenged (9, 8). Experimental studies in which species richness is manipulated give considerable support that species richness enhances resource use and thus productivity (20–22, 31). Moreover, a recent mathematical model explicitly shows that positive correlations between productivity and diversity may not be interpreted such that productivity drives diversity (9).
A causal relationship between diversity and RUE cannot be established based on observational data alone. The relationships between diversity and RUE shown here, however, are very consistent with those known from higher organisms (4, 6), and supports the view of general scaling rules for microbes and higher organisms (15). The data are also in good agreement with experimental data from artificial microbial communities (see above) but represent the first of its kind from natural species-rich microbial communities that have not been subject to experimental manipulation. Moreover, stabilizing effects of diversity and destabilizing effects of enrichment or trophic state ("paradox of enrichment;" ref. 28) have been hitherto reported independently from each other. To the best of our knowledge, however, this is the first study showing that stabilizing effects of diversity interact with resource levels in natural systems.
Our analysis indicates that RUE and community composition are both stabilized by diversity. A coupling between fluctuations in RUE and fluctuations of community composition seems plausible, because fluctuating RUE implies temporal windows where resources are not used efficiently. Such situations provide available niche-space for new species to become established in the system, pointing to higher risk of species invasion and/or resource monopolization (algal blooms) in species-poor communities. The inverse relationship between community turnover and diversity reported here (Fig. 3) is consistent with previous results for communities of higher organisms, including fast dispersing groups, such as macroalgae and zooplankton (32, 33). Therefore, we suggest that the stabilizing effect of diversity may be a universal rule.
A tight coupling between microbial diversity and ecosystem function implies that factors impairing microbial diversity are likely to affect efficiency and predictability of ecosystem processes. For example, pollution stress, especially from toxic substances, is often manifested as biodiversity loss. Our results indicate that pollution-mediated losses of microbial diversity may have direct effects on ecosystem processes, pointing to important linkages between toxic pollutants, nutrient enrichment, and biodiversity. This conclusion is supported by two recent studies on microbial communities in soils (34, 35) in which diversity of natural communities was artificially reduced by toxic compounds, and the metabolic rates were measured. A reduction of microbial diversity lowered community function and made the communities less resistant to further stress (34, 35).
We still know very little about the diversity of phytoplankton and other microbial communities and how it affects ecosystem functioning. This is due in part to methodological difficulties in assessment of microbial diversity (12). General ecological theories have only recently been applied to microbial communities, but the emerging patterns encourage a closer look into diversity of microbial assemblages as a driver of ecosystem function.
|
Materials and Methods |
|---|
|
TopAbstractResultsDiscussionMaterials and MethodsAcknowledgments.References |
|---|
For all datasets, only surface samples from the period between July and September were used to limit seasonal variation and to exclude periods with unfavorable weather conditions except for the time series analysis. Most lakes with multiple observations per year were sampled monthly; thus, we had to use samples from May to September to make this analysis meaningful. We excluded observations with TP concentrations <2 µg·liter–1 (n = 15) because they are close to the detection limit, resulting in very imprecise measurements.
Abundances and biovolumes of phytoplankton species and genera were estimated with the Utermöhl method. Phytoplankton carbon biomass (PPC) was estimated from biovolumes, applying a constant conversion factor of 14% of algal wet weight (36, 37).
Phytoplankton biomass divided by the limiting resource (total phosphorus) gave the yield or resource use efficiency (RUE) for the community. We estimated RUE both in terms of chlorophyll-a (RUEchl = chlorophyll-a:TP) and carbon biomass (RUEcarb = PPC:TP) (see also Validation of Underlying Assumptions). For the Baltic Sea, chlorophyll-a data were missing for many samples, causing a considerably lower number of observations for RUEchl (n = 99) compared with RUEcarb (n = 512).
Data Analysis. Regression analysis of RUE. Initial analyses indicated nonconstant error variances in regression models predicting RUE. We applied nonparametric bootstrapping (1,000 replications) to obtain robust estimates for linear regression model parameters (38). Statistical inferences on regression coefficients were based on their bootstrap distributions.
To avoid bias due to varying number of observations per locality, observations from localities with multiple observations were down-weighted by the inverse of the corresponding number of observations. Furthermore, for lakes with >10 observations, 10 observations were selected randomly from the total of observations.
Heteroscedasticity of the regression models was analyzed by using the Breusch–Pagan test (38), which tests whether the variability of the response variable scales with the predictors, i.e., whether the predictors affect the predictability of the relationship (38). Using the regression coefficients estimated from the bootstrapped regressions (see above), we calculated the log-transformed absolute residuals (v) and analyzed heteroscedasticity by v = a + b ln(G) + c ln(TP).
Species turnover and temporal variability of RUE.
Variability in community composition and RUE were analyzed for all lakes and years ("lake year") with a minimum of four observations per year. We calculated Bray–Curtis dissimilarity 
for all pairs of observations within a lake year, using square-root transformed genus-level data. Within a lake year, species turnover between two samples was a saturating function of the corresponding time interval (
t). Dividing by ln(t) gave a standardized s, which was not correlated with time. We calculated the averages of s, ln(G), ln(TP), and the standard deviation (SD) of ln(RUE) for each lake year. Using these data, we then estimated the average effects of G and TP on community turnover (s) and on variability (SD) of RUE in linear regressions. We did not estimate turnover rates for stations in the Baltic Sea, because sampling stations in an open system do not represent distinct systems, making changes in species composition influenced by horizontal advection rather than "local" dynamics.
Validation of Underlying Assumptions. There is a risk of confounding true effects with spurious correlations when analyzing relationships between two parameters (RUE and diversity) that both reflect cellular growth. Environmental factors acting upon phytoplankton growth might synchronize patterns seen in productivity and diversity. To minimize such effects, we excluded sites with extreme conditions from the analysis (see above). Using multiple regressions with diversity and important environmental factors as predictors of RUE, we double-checked that effects of diversity seen in the analyses do not represent spurious correlations driven by environmental factors. These regressions revealed that diversity is the best predictor of RUE in our datasets, whereas the explanatory power of other factors, such as pH, temperature, salinity, and lake morphometry, is rather poor (see SI).
Our definition of RUE depends on phosphorus being either directly limiting or at least not being in excess compared with the actual limiting nutrient. We are aware that phytoplankton is rarely limited by phosphorus alone, particularly at high TP-levels (24). Elser et al. (26) recently showed how synergistic effects of both phosphorus and nitrogen additions fuel the pelagic productivity of lakes in a consistent manner across a large number of studies. Because phosphorus is rarely available in great excess compared with nitrogen (27), TP should give a reasonable surrogate for potential productivity even in nitrogen limited systems. This assumption may, however, be questioned for eutrophic systems with predominant light limitation. To safeguard against potential biases that might emerge from including eutrophic sites in our analysis, we performed the analysis shown in Table 1 for the same dataset with observations >0.5 µmol of TP per liter being excluded (n = 1570 of 2535; temperate lakes with TP-levels <0.5 µmol of P per liter were shown to be predominantly P-limited (24)). The results of this analysis compare very well to the results shown for the entire dataset in the SI.
|
Acknowledgments. |
|---|
|
TopAbstractResultsDiscussionMaterials and MethodsAcknowledgments.References |
|---|
|
Footnotes |
|---|
To whom correspondence should be addressed. E-mail: robert.ptacnik{at}niva.noFreely available online through the PNAS open access option.
Author contributions: R.P. designed research; P.B., L.L., and E.W. performed research; R.P., A.G.S., and T.A. analyzed data; and R.P., A.G.S., T.A., T.T., and S.R. wrote the paper.
The authors declare no conflict of interest.
This article is a PNAS Direct Submission.
This article contains supporting information online at www.pnas.org/cgi/content/full/0708328105/DCSupplemental.
© 2008 by The National Academy of Sciences of the USA
|
References |
|---|
|
TopAbstractResultsDiscussionMaterials and MethodsAcknowledgments.References |
|---|
on MA, Bloem J, Romkens PFAM, de Ruiter PC (2006) Functional stability of microbial communities in contaminated soils near a Zinc smelter (Budel, the Netherlands) Ecotoxicology 15:187–197.[CrossRef][ISI][Medline]
CiteULike 
Complore 
Connotea 
Del.icio.us 
Digg What's this?
| ||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||
