Original Research ARTICLE
Microbes as Engines of Ecosystem Function: When Does Community Structure Enhance Predictions of Ecosystem Processes?
- 1Institute of Arctic and Alpine Research, University of Colorado Boulder, Boulder, CO, USA
- 2Biological Sciences Division, Pacific Northwest National Laboratory, Richland, WA, USA
- 3US Department of Energy, Joint Genome Institute, Walnut Creek, CA, USA
- 4Department of Forest Ecology, Federal Research and Training Centre for Forests, Bundesforschungs- und Ausbildungszentrum für Wald, Vienna, Austria
- 5Department of Soil Science, University of Saskatchewan, Saskatoon, SK, Canada
- 6Helmholtz Centre for Environmental Research – Centre for Environmental Biotechnology, Leipzig, Germany
- 7Department of Natural Resources and Environmental Sciences, University of Illinois at Urbana-Champaign, Urbana, IL, USA
- 8Life and Environmental Sciences and Sierra Nevada Research Institute, University of California – Merced, Merced, CA, USA
- 9School of Medicine, Flinders University, Adelaide, SA, Australia
- 10Institut National de la Recherche Agronomique – Agroecology, Dijon, France
- 11Institute of Biological and Environmental Sciences, University of Aberdeen, Aberdeen, UK
- 12Irstea, UR MALY, Centre de Lyon-Villeurbanne, Villeurbanne, France
- 13Department of Biogeography and Global Change, Museo Nacional de Ciencias Naturales, Consejo Superior de Investigaciones Científicas, Madrid, Spain
- 14Environment Centre Wales, Bangor University, Gwynedd, UK
- 15Department of Microbiology and Ecosystem Science, University of Vienna, Vienna, Austria
- 16Häme University of Applied Sciences, Hämeenlinna, Finland
- 17School of Freshwater Sciences, University of Wisconsin-Milwaukee, Milwaukee, WI, USA
- 18Department of Surface Waters, Eawag: Swiss Federal Institute of Aquatic Science and Technology, Kastanienbaum, Switzerland
- 19Centre for Carbon, Water and Food, The University of Sydney, Sydney, NSW, Australia
- 20Institute of Landscape Ecology, University of Münster, Münster, Germany
- 21Department of Environmental and Biological Sciences, University of Eastern Finland, Kuopio, Finland
- 22Natural Resources Institute, Vantaa, Finland
- 23Institute of Soil Science and Site Ecology, Technische University, Dresden, Germany
- 24Institute of Aquatic Ecology, Facultat de Ciències, University of Girona, Girona, Spain
- 25Institute for Sustainability Sciences – Agroscope, Zurich, Switzerland
- 26CSIRO Agriculture Flagship, Crace, ACT, Australia
- 27Department of Biology, University of Texas at Tyler, Tyler, TX, USA
- 28EDF R&D, National Hydraulics and Environmental Laboratory, Chatou, France
- 29Department of Microbiology and Ecosystem Science, University of Vienna, Vienna, Austria
- 30Division of Terrestrial Ecosystem Research, Department of Microbiology and Ecosystem Science, University of Vienna, Vienna, Austria
- 31Department of Ecosystem and Conservation Sciences, University of Montana, Missoula, MT, USA
- 32Centro de Investigación y Docencia Económicas – Consejo Superior de Investigaciones Científicas, Valencia, Spain
- 33Department of Biological Science, Virginia Institute of Marine Science, Gloucester Point, VA, USA
- 34AES School of Natural Resources Sciences, North Dakota State University, Fargo, ND, USA
- 35LEPABE - Laboratory for Process Engineering, Environmental, Biotechnology and Energy, Faculdade de Engenharia da Universidade do Porto, Porto, Portugal
- 36Southern California Coastal Water Research Project Authority, Costa Mesa, CA, USA
- 37UMR, Interactions Sol Plante Atmosphère, INRA Bordeaux, Villenave d’Ornon, France
- 38Department of Forest Mycology and Plant Pathology, Swedish University of Agricultural Sciences, Uppsala, Sweden
- 39Department of Biological Sciences, Virginia Polytechnic Institute, State University, Blacksburg, VA, USA
- 40Centre Tecnològic Forestal de Catalunya, Solsona, Spain
- 41Centre de Recerca Ecològica i Aplicacions Forestals, Cerdanyola del Vallès, Barcelona, Spain
- 42School of Civil and Environmental Engineering, Yonsei University, Seoul, South Korea
- 43Department of Applied Biological Chemistry, The University of Tokyo, Tokyo, Japan
- 44Department of Biomedicine, Biotechnology and Public Health, University of Cadiz, Puerto Real, Spain
- 45Research Centre for Agrobiology and Pedology, Florence, Italy
- 46Department of Ecology and Genetics/Limnology, Uppsala University, Uppsala, Sweden
- 47Vale Living with Lakes Centre and Department of Biology, Laurentian University, Sudbury, ON, Canada
- 48Biology Department, Duke University, Durham, NC, USA
Microorganisms are vital in mediating the earth’s biogeochemical cycles; yet, despite our rapidly increasing ability to explore complex environmental microbial communities, the relationship between microbial community structure and ecosystem processes remains poorly understood. Here, we address a fundamental and unanswered question in microbial ecology: ‘When do we need to understand microbial community structure to accurately predict function?’ We present a statistical analysis investigating the value of environmental data and microbial community structure independently and in combination for explaining rates of carbon and nitrogen cycling processes within 82 global datasets. Environmental variables were the strongest predictors of process rates but left 44% of variation unexplained on average, suggesting the potential for microbial data to increase model accuracy. Although only 29% of our datasets were significantly improved by adding information on microbial community structure, we observed improvement in models of processes mediated by narrow phylogenetic guilds via functional gene data, and conversely, improvement in models of facultative microbial processes via community diversity metrics. Our results also suggest that microbial diversity can strengthen predictions of respiration rates beyond microbial biomass parameters, as 53% of models were improved by incorporating both sets of predictors compared to 35% by microbial biomass alone. Our analysis represents the first comprehensive analysis of research examining links between microbial community structure and ecosystem function. Taken together, our results indicate that a greater understanding of microbial communities informed by ecological principles may enhance our ability to predict ecosystem process rates relative to assessments based on environmental variables and microbial physiology.
The links between complex environmental microbial communities and ecosystem processes remain unclear (Carney and Matson, 2005; Prosser et al., 2007; van der Heijden et al., 2008; Petersen et al., 2012; Wallenstein and Hall, 2012; Graham et al., 2014; Martiny et al., 2015), and an emerging field of research has begun to investigate the utility of microbial data for improving predictions of carbon (C) and nitrogen (N) cycling beyond estimates based solely on environmental data (Todd-Brown et al., 2012; Wieder et al., 2013; Reed et al., 2014; Powell et al., 2015). While researchers have attempted to enhance ecosystem process models by parameterizing microbial physiological properties such as drought tolerance (Manzoni et al., 2014), growth efficiency (Hagerty et al., 2014), dormancy (Wang et al., 2015), and turnover rates (Wieder et al., 2013), these models often fail to consider variation in microbial community structure that may regulate ecosystem process rates (Bouskill et al., 2012; Kaiser et al., 2014). As such, we still lack an integrated understanding of the interactions between microbial communities and ecosystem function, and a central question in ecosystem science remains: ‘Under what circumstances does information on microbial communities add to our predictive power of ecosystem processes?’ Addressing this question is essential not only for improving knowledge on how critical biogeochemical cycles may respond to current and impending environmental change, but also for allowing us to identify factors that determine microbial community structure and activity in space and time.
The specific conditions in which microbial community structure – broadly defined here as information on diversity and/or abundance of taxa within a community – may improve predictions of ecosystem process rates beyond models based on environmental and physiological attributes varies with a myriad of biotic and abiotic factors (Knelman and Nemergut, 2014; Nemergut et al., 2014). Research has demonstrated global patterns in microbial communities that correlate with environmental factors such as salinity (Lozupone and Knight, 2007; Auguet et al., 2010), pH (Lauber et al., 2009), and habitat type (Dinsdale et al., 2008; Nemergut et al., 2011; Fierer et al., 2012), and links between microbial community and ecosystem processes have been observed within numerous individual studies (van der Heijden et al., 1998; Torsvik and Øvreås, 2002; Bardgett and van der Putten, 2014; Wagg et al., 2014). However, these relationships and the underlying ecological principles that generate them may vary among environments, relevant microbial traits, and processes of interest (Wallenstein and Hall, 2012; Nemergut et al., 2014). For example, stochastic assembly processes (Hubbell, 2001; Rosindell et al., 2011; Stegen et al., 2012; Nemergut et al., 2013), phenotypic plasticity (DeWitt et al., 1998), and priority effects (Fukami, 2004; Fukami et al., 2010) can decouple microbial community structure from environmental conditions, and under such conditions, microbial community structure should be central to explaining ecosystem process rates. A unifying perspective that accounts for the dynamic relationships between microbial community structure and ecosystem function is vital for improving predictions of ecosystem process rates, yet recent work has been largely unable to decipher coherent trends in such relationships across studies (Rocca et al., 2014; Bier et al., 2015).
Thus, to assess the importance of microorganisms in explaining ecosystem processes, we used statistical modeling to evaluate the value of environmental variables and microbial community structure both alone and in conjunction for explaining rates of C and N cycling processes within global datasets. Although the literature reporting environmental variables, microbial community data, and biogeochemical processes in a single study is sparse (Rocca et al., 2014), we collected and re-analyzed 82 datasets spanning a multitude of diverse environments. We hypothesized that while the environment is a strong driver of most ecosystem processes, microbial community structure adds to our understanding of processes under certain circumstances, and we discuss patterns across ecological, biogeochemical, and phylogenetic subsets of data. Our analysis generates new insights about the relationship between microbial community structure and a variety of ecosystem processes.
Materials and Methods
We completed an exhaustive literature review for studies measuring environmental variables, microbial community structure, and biogeochemical processes across a range of ecosystems (Supplementary Table S1). We then constructed statistical models for three predictor sets and their additive combinations – environmental data, microbial community data, and microbial biomass measurements (Supplementary Figure S1). Because the measurement of microbial biomass (defined as microbial C and/or N content) was not a criterion for dataset selection, only 28 datasets included biomass measurements. To incorporate differences in microbial community structure into our models, we used ordination-based techniques to condense multivariate community data into vectors that we included as predictors in our models (Supplementary Methods). The Shannon diversity index of each sample was also included as a measure of alpha diversity. In total, these taxonomic metrics are referred to as ‘community diversity’ when applied to 16S rRNA genes or PLFAs and ‘functional diversity’ when applied to functional genes. Functional gene abundances were included as measures of total gene abundance per gram of soil when possible (n = 17), although five datasets only listed gene abundances as normalized relative to 16S rRNA gene abundance.
Within each type of data, methodology to collect environmental, microbial, and process data, as well as variables collected, inevitably varied. To minimize potential error introduced by these differences, we re-analyzed data within each dataset using a multimodel inference approach with the ‘MuMIn’ package (Barton, 2011) in R software (R Core Team, 2014) and subsequently compared results from our analyses. Multimodel inference is a broad regression-based, model-averaging statistical approach designed to reduce errors in model selection, and this method has the advantage of standardizing our approach across studies while accounting for a lack of a priori system-specific knowledge (Burnham and Anderson, 2002). This statistical approach has been used in other small scale studies with similar objectives (Powell et al., 2015). We validated the accuracy of MuMIn in our dataset by comparing models for several datasets with expert-built regression models, which yielded comparable results (more details available in Supplementary Methods).
Using the ‘dredge’ command in the ‘MuMIn’ package to fit and evaluate the explanatory power of all possible combinations of variables within a predictor set on a process rate, we selected a set of best fit models for each predictor set consisting of all models with a delta AICc value no more than four greater than the model with the lowest AICc value (Burnham and Anderson, 2004). We generated an averaged final model from this model set using Akaike’s weights, implemented with the ‘model.avg’ command in R (Supplementary Methods). Final models from different predictor sets were compared for statistically significant differences at a delta AICc value of four to provide conservative estimates of model improvements (Burnham and Anderson, 2004). Models were also evaluated for ecologically relevant improvement, defined by an increase in adjusted R2 value greater than 10% of the environmental model adjusted R2 value. This criterion was implemented to remove artifacts from datasets in which statistically different models, based on AICc values, yielded similar adjusted R2 values. Only models that showed both statistical and ecological improvement were considered to be improved.
Finally, model results were synthesized across studies within ecologically relevant subsets of data. We examined results within the full dataset and within biogeochemical process [C mineralization (referred to here as ‘respiration’), nitrification, denitrification, N mineralization], microbial data (PLFA, tRFLP, ARISA, DGGE, qPCR, next generation sequencing) and environment types (natural soil, sediment, agricultural soil) with sufficient replication (n ≥ 12). We report the mean adjusted R2 of models and the increase in adjusted R2 value as measures of effect size as well as the percent of models statistically improved by the addition of microbial data relative to models constructed with only environmental parameters. Differences among the explanatory power of models with environmental, microbial, or both environmental and microbial data were assessed by comparing the distribution of model adjusted R2 values within each predictor set using unpaired one-sided Mann–Whitney U-tests for non-parametric data. We also examined correlations between predictor sets to determine the extent to which environmental variables explained variation in microbial community structure and biomass and to which biomass explained variation in microbial community structure. We analyzed these relationships with redundancy analysis (RDA), utilizing the ‘ordistep’ function in the ‘vegan’ package (Oksanen et al., 2013) in R to automate forward model selection, and we report the average R2 values of correlations between predictor sets as a measure of effect size (Supplementary Methods).
Overall, models based on environmental variables explained significantly more variation in processes than models based on microbial community structure (n = 82, average adjusted R2 0.56 vs. 0.31, Mann–Whitney U, p < 0.0001, Figure 1). On average, models with both environmental and microbial predictors explained more variation in processes than environmental models (n = 82, average adjusted R2 0.65 vs. 0.56, Mann–Whitney U, p = 0.046); however, only 29% of datasets were significantly improved by adding information on microbial community structure, by an average of 0.08 increase in adjusted R2 within all models. Data on microbial community structure from ARISA, tRFLP, qPCR, and next generation technology did not differ in their explanatory power of process rates, but DGGE and PLFA had significantly lower explanatory power than other metrics (Supplementary Figure S2A). All microbial data types showed weak correlations to environmental variables but gene abundance data displayed a trend for higher correlation (n = 22, average adjusted R2 0.36) than community (n = 55, average adjusted R2 0.20) or functional diversity (n = 5, average adjusted R2 0.21) metrics (Supplementary Figure S2B).
FIGURE 1. Within the full dataset (n = 82), environmental variables explained more variation in ecosystem process rates than microbial community structure. Microbial predictors alone had low explanatory power on processes but improved model explanatory power on average. Error bars denote standard errors, and letters represent significantly different groups (Mann–Whitney U, p < 0.05).
We also examined each process individually for those functions for which we had the most data: nitrification (n = 14), denitrification (n = 17), N mineralization (n = 12), and respiration (n = 26). Overall, variation in N mineralization and denitrification rates was less well-explained than other processes by any predictor set (Figure 2). Only 17% of datasets examining N mineralization were improved by microbial community data – an average increase in adjusted R2 of 0.008 within all N mineralization datasets – and the average explanatory power of microbial community structure alone on N mineralization rates was 0.21, compared to 0.31 in the full dataset. Within N mineralization studies that assayed microbial biomass, we found that microbial biomass correlated with environmental variables (n = 5, average adjusted R2 0.58) but not with process rates; none of the N mineralization models were improved with the addition of data on biomass. No denitrification studies examined microbial biomass or PLFAs. However, despite low replication, it is notable that in three of four studies, models of denitrification were improved with the addition of 16S rRNA gene diversity data (0.13 average increase in adjusted R2), while only three of eleven models were improved with data on functional gene abundance (0.04 average increase in adjusted R2).
FIGURE 2. Average explanatory power varied by process and microbial data type, and we present the explanatory power of functional genes (abundance and diversity) and community diversity for each N cycling process. Environmental model bars slightly differ within each figure as they represent only the environmental models for each data type, and functional genes are not depicted for N mineralization due to low sample size (n = 2). Error bars represent standard errors. (A) N mineralization (n = 12) and (B) denitrification (n = 17) models had relatively low explanatory by all predictor sets. However, community diversity metrics (n = 4) provided more added value for denitrification rates than for other processes. (C) By contrast, rates of nitrification (n = 14) were well-explained by both environmental and functional gene predictor sets (n = 6) and were more likely to be improved by functional data.
Rates of nitrification were well-explained by all predictor sets (Figure 2). Nitrification models based on microbial data yielded the highest average adjusted R2 of any process (0.38 across all data types and 0.50 if PLFA data are removed) with gene abundance (n = 4) data showing higher adjusted R2 than community diversity data (n = 8; 0.61 vs. 0.21, Mann–Whitney U, p = 0.008). Moreover, 50% of environmental models of nitrification were improved with data on functional gene abundance or functional diversity (n = 6), while only 30% of models were improved with data on 16S rRNA gene diversity (n = 6), and PLFA data only improved 20% of models (n = 2). As with N mineralization, biomass was correlated with environmental variables (n = 3, average adjusted R2 0.62) and did not add explanatory power to nitrification rates.
Respiration rates were well-explained by the environment, with an average adjusted R2 of 0.66 (Figure 3). Similar to the full dataset, microbial-only respiration models showed an average adjusted R2 of 0.29 and only 23% of respiration models were improved with the addition of 16S rRNA gene diversity or PLFA data (no functional gene abundance datasets existed), by an average increase in 0.06 in adjusted R2 across all respiration datasets. Interestingly, microbial biomass, when measured, improved 35% of models of respiration rates (n = 17, 0.09 average increase in adjusted R2), while 53% of models were improved with data on biomass and community structure in combination (0.15 average increase in adjusted R2). Microbial biomass was correlated to environmental variables with an average adjusted R2 of 0.56 but not to microbial community structure (average adjusted R2 0.11).
FIGURE 3. Rates of respiration (n = 26) were well-explained by all predictor sets relative to the full dataset. In this figure, we present microbial biomass in addition to microbial community structure, with standard errors represented as error bars. Microbial community structure and biomass independently had low explanatory power compared to environmental variables; however, they were largely uncorrelated with each other and provided greater model improvement when added to environmental variables in conjunction than when added independently.
Although most of our datasets were from natural soil studies (n = 47), we had sufficient replication of sediment (n = 12) and agricultural soil (n = 12) datasets to examine these independently; other environments (n < 12) were excluded from this analysis. While statistically significant differences were not found, possibly due to the large variation in our datasets, we observed some notable trends across environments. Datasets from sediment samples were less well-explained by environmental or microbial predictor sets than other datasets (Figure 4). Relationships between community diversity and environmental conditions were also weak in sediment systems (n = 8, average adjusted R2 0.17), and community diversity metrics provided more added value (0.12 average increase in adjusted R2) in sediments than other systems. Agricultural systems had high explanatory power of process by both environmental (average adjusted R2 0.59) and microbial predictors (average adjusted R2 0.38) relative to other systems, and metrics of community diversity were more strongly correlated with the environment (n = 7, average adjusted R2 0.24) and had lower added value (average adjusted R2 0.06) than in other ecosystem types.
FIGURE 4. We analyzed natural soil (n = 47), sediment (n = 12), and agricultural soil (n = 12) environments separately. Bar height represents the average explanatory power of each predictor set, with error bars denoting standard errors. While predictor sets in natural soils nearly matched those within the full dataset, environmental and microbial predictor sets had lower explanatory power within sediments and higher explanatory power within agricultural soils. Microbial community diversity was less tightly coupled to the environment and provided more model improvement in sediments, with the opposite relationship in agricultural soils.
Despite the vast amount of variation in our datasets, our results indicate that data on microbial community structure may strengthen predictions of C and N cycling rates. Within our full dataset, environmental variables alone had greater explanatory power than either metrics of microbial community structure or microbial biomass but left 44% of variation in process rates unexplained on average (Figure 1). While we acknowledge that factors such as unmeasured variables, macroinvertebrates (Ferris et al., 2004), fungi (Talbot et al., 2008), plant communities (Hanson et al., 2000), and non-linearity could reduce the explanatory power of our models, a number of ecological dynamics can disassociate environmental conditions from microbial communities mediating C and N cycle processes and alter the relative value of environmental and microbial data for explaining rates of ecosystem processes (Knelman and Nemergut, 2014; Nemergut et al., 2014). When such factors strongly impact microbial communities, data on microbial community structure should enhance our predictions of ecosystem processes beyond those based solely on environmental data.
We observed stronger trends in the explanatory power of microbial predictors on ecosystem process rates when we reduced noise by examining subsets of data containing a single process or environment type. However, the lack of consistent trends within the full dataset is not surprising due to a variety of ecological and methodological factors. For instance, microorganisms can enter a state of lowered activity or dormancy in response to unfavorable environmental conditions (Jones and Lennon, 2010; Lennon and Jones, 2011), and research suggests that the percent of soil microbial communities catalyzing chemical reactions in soils at a given point in time can vary widely with resource history and disturbances (Blagodatskaya and Kuzyakov, 2013). These factors also induce spatially and temporally variable relationships between microbial communities and ecosystem processes, and the same microbial community may display different relationships to ecosystem function when sampled through time. Moreover, the extent to which functional traits are phylogenetically conserved varies with the trait of interest (Martiny et al., 2013), and thus, some processes may exhibit inherently stronger correlations with taxonomic metrics than others. Although, rapidly improving transcriptomic and proteomic sequencing technology may enhance predictions of ecosystem processes by identifying active segments of microbial communities (Abram, 2015), too few datasets presently exist to conduct a robust analysis with these data. Studies investigating the extent to which these dynamics are masked within observations of microbial community structure and function are imperative to deciphering and improving the current ability of microbial data to predict ecosystem process rates.
As well, discrepancies in methodology may also confound the relationship we observed in the full dataset between microbial communities and biogeochemical processes (Sinsabaugh et al., 2015). Microbial DNA is typically extracted from a fraction of gram of soil at a single time point and methods detect both dormant and active organisms, while corresponding process data are measured in several grams of soil over longer timeframes. Further, process rate measures are either collected as net measurements that aggregate multiple chemical reactions, potential rates that represent microbial response to substrate addition, or gross rates that trace individual chemical transformations (Bier et al., 2015). These approaches inevitably introduce data aggregation issues (Clark et al., 2011) and rely on assumptions that microbial and biogeochemical data are spatially and temporally representative. Recently, studies have attempted to reduce uncertainty associated with scalar differences by pairing spatially and temporally explicit field studies with ecological modeling approaches that interpolate data across scales (Lira-Noriega et al., 2013; Walker and Wardle, 2014). Leveraging these approaches in conjunction requires substantial resource investment but is increasingly feasible and presents a promising avenue for maximizing the value of microbial data.
Despite inherent limitations in microbial data, our data analysis suggests that microbial data types differ in their ability to explain processes that are phylogenetically broad vs. narrow and that are products of facultative vs. obligate metabolisms (Figure 2). Nitrification is the product of an obligate chemolithotrophic metabolism, and functional markers for this process are narrowly constrained within ammonia-oxidizing bacteria and archaea that catalyze the rate-limiting step in nitrification, putatively representing a likely dataset for detecting relationships between microbial data and process rates (Schimel, 1995; Schimel and Gulledge, 1998). While both functional gene abundance and community diversity metrics had similar average explanatory power within the full dataset (an effect that also did not vary among microbial technologies, Supplementary Figure S2A), models based on functional gene abundances or functional diversity explained more variation in process rates than community diversity within datasets measuring rates of nitrification (Figure 2C). Although these metrics were more strongly correlated with environmental factors than community diversity, they also provided greater explanatory power in combination with environmental variables than community diversity metrics, which include variation in many microbial guilds external to nitrification. Notably, models based on environmental variables alone also had higher explanatory power for nitrification than for denitrification or N mineralization, suggesting that the environment may also be more strongly coupled to nitrification rates than other N cycle processes. Taken together, these results suggest that despite a relatively tight connection between the environment, microbial community structure, and process data, the microbial guilds driving nitrification rates are partially decoupled at a functional level from environmental drivers of nitrification. Importantly, although nitrification was the only narrow process in our dataset with sufficient replication to examine independently, the high explanatory power of functional genes, both alone and in combination with environmental variables, for describing nitrification rates relative to other processes also conveys the potential for the inclusion of functional metrics in improving predictions of ecosystem processes that are mediated by phylogenetically narrow functional guilds.
By contrast, denitrification is a broad process, catalyzed by diverse facultative anaerobic microorganisms, and community diversity metrics provided more added value to denitrification rates than functional gene diversity or abundances (Figure 2B). While we acknowledge that genetic markers may not denote all functional genes involved in a process (Jones et al., 2014; Verbaendert et al., 2014), we observed more added value from community diversity metrics than functional genes despite several studies assaying up to three genes to explain denitrification rates, as opposed to one or two functional markers measured for nitrification. For broad processes such as denitrification, representations of niche complementarity or resource availability reflected in diversity metrics (Cadotte et al., 2011) rather than functional gene abundances that do not account for interactions between individuals (Barberán et al., 2012; Gagic et al., 2015) or functional redundancy (Allison and Martiny, 2008), may be important for explaining process rates. Thus, it is plausible that denitrifying communities are decoupled from environmental variables when measured at a broad taxonomic level and that such taxonomic metrics encompass more variation in within community interactions affecting denitrification rates than functional approaches. Although the sample size of datasets measuring denitrification was small, our analysis suggests that measurements of community diversity may be more likely than functional metrics to improve predictions of ecosystem processes that are catalyzed by facultative metabolisms and/or phylogenetically broad suites of organisms.
Similarly, the ability to respire carbon is widely distributed among microorganisms, and community diversity metrics (no datasets measured functional markers) described rates of respiration better than all other processes except denitrification (Figure 3). In fact, environmental and microbial variables, both alone and in combination, yielded higher average explanatory power of respiration rates than all other processes we examined independently with the exception of nitrification. Importantly, while recent work has suggested that microbial biomass (Schimel and Weintraub, 2003) or physiological properties such as drought tolerance (Manzoni et al., 2014) and growth efficiency (Wieder et al., 2013) can improve C cycling models, our results indicate that information on microbial community structure may further enhance our understanding of ecosystem C cycling. Microbial community structure and biomass appeared to explain differing portions of variation in process, as they exhibited extremely weak correlations with each other and, together, provided additive value for explaining rates of respiration. Thus, biomass and structure may jointly contribute to explaining rates of respiration as biomass may serve as a proxy for unmeasured environmental variables that regulate activity, such as soil structure (Gupta and Germida, 1988), while community structure may relate to the genetic capability of a community to respire carbon. Regardless, our results suggest that these data are at least partially independent and that measurements of microbial community structure could reduce uncertainty within Earth Systems models.
Lastly, the relationship between microbial community structure and ecosystem processes in our dataset varied by environment type (Figure 4). Sediment systems had weaker explanatory power of process rates by both environmental variables and microbial community structure than natural or agricultural soils. Sediment processes are influenced by both sediment and porewater chemistry (Middelburg and Levin, 2009), and these datasets may be more likely to have unmeasured variables or errors due to spatial discrepancies than other systems. In particular, dissolved oxygen concentration is a strong determinant of redox potential and ecosystem processes in sediments (Abell et al., 2011) but was only measured in three of 12 datasets. Despite the low explanatory power of sediment processes by all predictor sets, microbial community diversity added more value to explaining process rates and was more decoupled from environmental variables in sediments than in natural soil or agricultural systems, indicating that measurements of microbial community structure may improve predictions of sediment processes. Conversely, in agricultural soils, microbial community structure explained more variation in ecosystem processes but provided less improvement over environmental variables than other systems. Environmental variables were also better predictors of process in agricultural soils, and microbial community diversity was more tightly correlated with the environment than in sediments or natural soils. Agricultural systems homogenize soil structure and decrease soil organic matter (Dick, 1992; Chan and Heenan, 1996), reducing the number and complexity of microbial niches in agricultural soils. Thus, relatively coarse-scale measures of environmental conditions may more directly correlate with both community structure and ecosystem process rates than in more variable systems. Overall, these results indicate that ecosystem-specific dynamics may be crucial to understanding the value of microbial community structure for explaining ecosystem processes and emphasize the importance of future investigations into understudied biomes in enhancing our understanding of global relationships between microbial communities and ecosystem processes.
Here, we present the first extensive investigation into the relationship between environmental conditions, microbial community structure, and ecosystem function by re-analyzing 82 datasets collected from an international team of collaborators, and we demonstrate that nuanced relationships between the environment and microbial communities can influence the importance of microbial community structure for explaining ecosystem-level processes. Our analysis provides an empirical basis for future hypothesis-testing on the roles of assembly, dormancy, redundancy and phenotypic plasticity in microbial community structure and function, and despite complexity in our datasets, the trends we observed represent an encouraging step for linking microbial community data to ecosystem function. As a whole, our results indicate that a greater understanding of microbial community structure informed by ecological principles may further our ability to accurately predict rates of ecosystem processes beyond environmental variables and bulk physiological characterizations of microbial communities.
All authors collected and contributed datasets for analysis as well as participated in the conceptual drafting and revision of this manuscript. EG conducted all data analysis and was the primary author in writing and revising the manuscript. DN and JK contributed significantly in manuscript development and revision, and DN provided financial support for this work.
This work was supported by NSF grant DEB-1221215 to DN, as well as grants supporting the generation of our datasets as acknowledged in their original publications and in Supplementary Table S1.
Conflict of Interest Statement
The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.
We thank the USGS Powell Center ‘Next Generation Microbes’ working group, anonymous reviews, Brett Melbourne, and Alan Townsend for valuable feedback on this project.
The Supplementary Material for this article can be found online at: https://www.frontiersin.org/article/10.3389/fmicb.2016.00214
Abell, G. C., Banks, J., Ross, D. J., Keane, J. P., Robert, S. S., Revill, A. T., et al. (2011). Effects of estuarine sediment hypoxia on nitrogen fluxes and ammonia oxidizer gene transcription. FEMS Microbiol. Ecol. 75, 111–122. doi: 10.1111/j.1574-6941.2010.00988.x
Barberán, A., Bates, S. T., Casamayor, E. O., and Fierer, N. (2012). Using network analysis to explore co-occurrence patterns in soil microbial communities. ISME J. 6, 343–351. doi: 10.1038/ismej.2011.119
Barton, K. (2011). MuMIn: Multi-Model Inference. R Package Version 1.0. 0. Vienna, Austria: R Foundation for Statistical Computing. Available at: http://CRAN.R-project.org/package=MuMIn
Bier, R. L., Bernhardt, E. S., Boot, C. M., Graham, E. B., Hall, E. K., Lennon, J. T., et al. (2015). Linking microbial community structure and microbial processes: an empirical and conceptual overview. FEMS Microbiol. Ecol. 91:fiv113. doi: 10.1093/femsec/fiv113
Bouskill, N. J., Tang, J., Riley, W. J., and Brodie, E. L. (2012). Trait-based representation of biological nitrification: model development, testing, and predicted community composition. Front. Microbiol. 3:364. doi: 10.3389/fmicb.2012.00364
Cadotte, M. W., Carscadden, K., and Mirotchnick, N. (2011). Beyond species: functional diversity and the maintenance of ecological processes and services. J. Appl. Ecol. 48, 1079–1087. doi: 10.1111/j.1365-2664.2011.02048.x
Carney, K. M., and Matson, P. A. (2005). Plant communities, soil microorganisms, and soil carbon cycling: does altering the world belowground matter to ecosystem functioning? Ecosystems 8, 928–940. doi: 10.1007/s10021-005-0047-0
Chan, K., and Heenan, D. (1996). The influence of crop rotation on soil structure and soil physical properties under conventional tillage. Soil Tillage Res. 37, 113–125. doi: 10.1016/0167-1987(96)01008-2
Clark, J. S., Bell, D. M., Hersh, M. H., Kwit, M. C., Moran, E., Salk, C., et al. (2011). Individual-scale variation, species-scale differences: inference needed to understand diversity. Ecol. Lett. 14, 1273–1287. doi: 10.1111/j.1461-0248.2011.01685.x
Ferris, H., Venette, R., and Scow, K. (2004). Soil management to enhance bacterivore and fungivore nematode populations and their nitrogen mineralisation function. Appl. Soil Ecol. 25, 19–35. doi: 10.1016/j.apsoil.2003.07.001
Fierer, N., Leff, J. W., Adams, B. J., Nielsen, U. N., Bates, S. T., Lauber, C. L., et al. (2012). Cross-biome metagenomic analyses of soil microbial communities and their functional attributes. Proc. Natl. Acad. Sci. U.S.A. 109, 21390–21395. doi: 10.1073/pnas.1215210110
Fukami, T., Dickie, I. A., Paula Wilkie, J., Paulus, B. C., Park, D., Roberts, A., et al. (2010). Assembly history dictates ecosystem functioning: evidence from wood decomposer communities. Ecol. Lett. 13, 675–684. doi: 10.1111/j.1461-0248.2010.01465.x
Gagic, V., Bartomeus, I., Jonsson, T., Taylor, A., Winqvist, C., Fischer, C., et al. (2015). Functional identity and diversity of animals predict ecosystem functioning better than species-based indices. Proc. R. Soc. Biol. Sci. U.S.A. 282:20142620. doi: 10.1098/rspb.2014.2620
Graham, E. B., Wieder, W. R., Leff, J. W., Weintraub, S. R., Townsend, A. R., Cleveland, C. C., et al. (2014). Do we need to understand microbial communities to predict ecosystem function? A comparison of statistical models of nitrogen cycling processes. Soil Biol. Biochem. 68, 279–282. doi: 10.1016/j.soilbio.2013.08.023
Gupta, V., and Germida, J. (1988). Distribution of microbial biomass and its activity in different soil aggregate size classes as affected by cultivation. Soil Biol. Biochem. 20, 777–786. doi: 10.1016/0038-0717(88)90082-X
Hagerty, S. B., van Groenigen, K. J., Allison, S. D., Hungate, B. A., Schwartz, E., Koch, G. W., et al. (2014). Accelerated microbial turnover but constant growth efficiency with warming in soil. Nat. Clim. Chang. 4, 903–906. doi: 10.1038/nclimate2361
Hanson, P., Edwards, N., Garten, C., and Andrews, J. (2000). Separating root and soil microbial contributions to soil respiration: a review of methods and observations. Biogeochemistry 48, 115–146. doi: 10.1023/A:1006244819642
Jones, C. M., Spor, A., Brennan, F. P., Breuil, M.-C., Bru, D., Lemanceau, P., et al. (2014). Recently identified microbial guild mediates soil N2O sink capacity. Nat. Clim. Chang. 4, 801–805. doi: 10.1038/nclimate2301
Knelman, J. E., and Nemergut, D. R. (2014). Changes in community assembly may shift the relationship between biodiversity and ecosystem function. Front. Microbiol. 5:424. doi: 10.3389/fmicb.2014.00424
Lauber, C. L., Hamady, M., Knight, R., and Fierer, N. (2009). Pyrosequencing-based assessment of soil pH as a predictor of soil bacterial community structure at the continental scale. Appl. Environ. Microbiol. 75, 5111–5120. doi: 10.1128/AEM.00335-09
Manzoni, S., Schaeffer, S., Katul, G., Porporato, A., and Schimel, J. (2014). A theoretical analysis of microbial eco-physiological and diffusion limitations to carbon cycling in drying soils. Soil Biol. Biochem. 73, 69–83. doi: 10.1016/j.soilbio.2014.02.008
Nemergut, D. R., Costello, E. K., Hamady, M., Lozupone, C., Jiang, L., Schmidt, S. K., et al. (2011). Global patterns in the biogeography of bacterial taxa. Environ. Microbiol. 13, 135–144. doi: 10.1111/j.1462-2920.2010.02315.x
Nemergut, D. R., Schmidt, S. K., Fukami, T., O’Neill, S. P., Bilinski, T. M., Stanish, L. F., et al. (2013). Patterns and processes of microbial community assembly. Microbiol. Mol. Biol. Rev. 77, 342–356. doi: 10.1128/MMBR.00051-12
Oksanen, J., Blanchet, F. G., Kindt, R., Legendre, P., Minchin, P. R., O’Hara, R., et al. (2013). Package ‘Vegan’. Community Ecology Package, Version 2. Available at: http://CRAN.R-project.org/package=vegan
Petersen, D. G., Blazewicz, S. J., Firestone, M., Herman, D. J., Turetsky, M., and Waldrop, M. (2012). Abundance of microbial genes associated with nitrogen cycling as indices of biogeochemical process rates across a vegetation gradient in Alaska. Environ. Microbiol. 14, 993–1008. doi: 10.1111/j.1462-2920.2011.02679.x
Powell, J. R., Welsh, A., and Hallin, S. (2015). Microbial functional diversity enhances predictive models linking environmental parameters to ecosystem properties. Ecology 96, 1985–1993. doi: 10.1890/14-1127.1
Prosser, J. I., Bohannan, B. J., Curtis, T. P., Ellis, R. J., Firestone, M. K., Freckleton, R. P., et al. (2007). The role of ecological theory in microbial ecology. Nat. Rev. Microbiol. 5, 384–392. doi: 10.1038/nrmicro1643
Reed, D. C., Algar, C. K., Huber, J. A., and Dick, G. J. (2014). Gene-centric approach to integrating environmental genomics and biogeochemical models. Proc. Natl. Acad. Sci. U.S.A. 111, 1879–1884. doi: 10.1073/pnas.1313713111
Rocca, J. D., Hall, E. K., Lennon, J. T., Evans, S. E., Waldrop, M. P., Cotner, J. B., et al. (2014). Relationships between protein-encoding gene abundance and corresponding process are commonly assumed yet rarely observed. ISME J. 9, 1693–1699. doi: 10.1038/ismej.2014.252
Schimel, J. (1995). “Ecosystem consequences of microbial diversity and community structure,” in Arctic and Alpine Biodiversity: Patterns, Causes and Ecosystem Consequences, ed. C. Krner (Berlin: Springer), 239–254.
Schimel, J. P., and Weintraub, M. N. (2003). The implications of exoenzyme activity on microbial carbon and nitrogen limitation in soil: a theoretical model. Soil Biol. Biochem. 35, 549–563. doi: 10.1016/S0038-0717(03)00015-4
Sinsabaugh, R. L., Shah, J. J. F., Findlay, S. G., Kuehn, K. A., and Moorhead, D. L. (2015). Scaling microbial biomass, metabolism and resource supply. Biogeochemistry 122, 175–190. doi: 10.1007/s10533-014-0058-z
Stegen, J. C., Lin, X., Konopka, A. E., and Fredrickson, J. K. (2012). Stochastic and deterministic assembly processes in subsurface microbial communities. ISME J. 6, 1653–1664. doi: 10.1038/ismej.2012.22
Talbot, J., Allison, S., and Treseder, K. (2008). Decomposers in disguise: mycorrhizal fungi as regulators of soil C dynamics in ecosystems under global change. Funct. Ecol. 22, 955–963. doi: 10.1111/j.1365-2435.2008.01402.x
Todd-Brown, K. E., Hopkins, F. M., Kivlin, S. N., Talbot, J. M., and Allison, S. D. (2012). A framework for representing microbial decomposition in coupled climate models. Biogeochemistry 109, 19–33. doi: 10.1007/s10533-011-9635-6
van der Heijden, M. G. A., Bardgett, R. D., and van Straalen, N. M. (2008). The unseen majority: soil microbes as drivers of plant diversity and productivity in terrestrial ecosystems. Ecol. Lett. 11, 296–310. doi: 10.1111/j.1461-0248.2007.01139.x
van der Heijden, M. G., Klironomos, J. N., Ursic, M., Moutoglis, P., Streitwolf-Engel, R., Boller, T., et al. (1998). Mycorrhizal fungal diversity determines plant biodiversity, ecosystem variability and productivity. Nature 396, 69–72. doi: 10.1038/23932
Verbaendert, I., Hoefman, S., Boeckx, P., Boon, N., and De Vos, P. (2014). Primers for overlooked nirK, qnorB, and nosZ genes of thermophilic Gram-positive denitrifiers. FEMS Microbiol. Ecol. 89, 162–180. doi: 10.1111/1574-6941.12346
Wagg, C., Bender, S. F., Widmer, F., and van der Heijden, M. G. (2014). Soil biodiversity and soil community composition determine ecosystem multifunctionality. Proc. Natl. Acad. Sci. U.S.A. 111, 5266–5270. doi: 10.1073/pnas.1320054111
Wallenstein, M. D., and Hall, E. K. (2012). A trait-based framework for predicting when and where microbial adaptation to climate change will affect ecosystem functioning. Biogeochemistry 109, 35–47. doi: 10.1007/s10533-011-9641-8
Wang, G., Jagadamma, S., Mayes, M. A., Schadt, C. W., Steinweg, J. M., Gu, L., et al. (2015). Microbial dormancy improves development and experimental validation of ecosystem model. ISME J. 9, 226–237. doi: 10.1038/ismej.2014.120
Keywords: microbial diversity, functional gene, statistical modeling, microbial ecology, ecosystem processes, respiration, nitrification, denitrification
Citation: Graham EB, Knelman JE, Schindlbacher A, Siciliano S, Breulmann M, Yannarell A, Beman JM, Abell G, Philippot L, Prosser J, Foulquier A, Yuste JC, Glanville HC, Jones DL, Angel R, Salminen J, Newton RJ, Bürgmann H, Ingram LJ, Hamer U, Siljanen HMP, Peltoniemi K, Potthast K, Bañeras L, Hartmann M, Banerjee S, Yu R-Q, Nogaro G, Richter A, Koranda M, Castle SC, Goberna M, Song B, Chatterjee A, Nunes OC, Lopes AR, Cao Y, Kaisermann A, Hallin S, Strickland MS, Garcia-Pausas J, Barba J, Kang H, Isobe K, Papaspyrou S, Pastorelli R, Lagomarsino A, Lindström ES, Basiliko N and Nemergut DR (2016) Microbes as Engines of Ecosystem Function: When Does Community Structure Enhance Predictions of Ecosystem Processes? Front. Microbiol. 7:214. doi: 10.3389/fmicb.2016.00214
Received: 09 December 2015; Accepted: 09 February 2016;
Published: 24 February 2016.
Edited by:Gary M. King, Louisiana State University, USA
Reviewed by:Steffen Kolb, Landscape Biogeochemistry – Leibniz Centre for Agricultural Landscape Research, Germany
Hongchen Jiang, Miami University, USA
Kristen M. DeAngelis, University of Massachusetts Amherst, USA
Copyright © 2016 Graham, Knelman, Schindlbacher, Siciliano, Breulmann, Yannarell, Beman, Abell, Philippot, Prosser, Foulquier, Yuste, Glanville, Jones, Angel, Salminen, Newton, Bürgmann, Ingram, Hamer, Siljanen, Peltoniemi, Potthast, Bañeras, Hartmann, Banerjee, Yu, Nogaro, Richter, Koranda, Castle, Goberna, Song, Chatterjee, Nunes, Lopes, Cao, Kaisermann, Hallin, Strickland, Garcia-Pausas, Barba, Kang, Isobe, Papaspyrou, Pastorelli, Lagomarsino, Lindström, Basiliko and Nemergut. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) or licensor are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.
*Correspondence: Emily B. Graham, firstname.lastname@example.org