Genus-Wide Transcriptional Landscapes Reveal Correlated Gene Networks Underlying Microevolutionary Divergence in Diatoms.

Marine microbes like diatoms make up the base of marine food webs and drive global nutrient cycles. Despite their key roles in ecology, biogeochemistry, and biotechnology, we have limited empirical data on how forces other than adaptation may drive diatom diversification, especially in the absence of environmental change. One key feature of diatom populations is frequent extreme reductions in population size, which can occur both in situ and ex situ as part of bloom-and-bust growth dynamics. This can drive divergence between closely related lineages, even in the absence of environmental differences. Here, we combine experimental evolution and transcriptome landscapes (t-scapes) to reveal repeated evolutionary divergence within several species of diatoms in a constant environment. We show that most of the transcriptional divergence can be captured on a reduced set of axes, and that repeatable evolution can occur along a single major axis of variation defined by core ortholog expression comprising common metabolic pathways. Previous work has associated specific transcriptional changes in gene networks with environmental factors. Here, we find that these same gene networks diverge in the absence of environmental change, suggesting these pathways may be central in generating phenotypic diversity as a result of both selective and random evolutionary forces. If this is the case, these genes and the functions they encode may represent universal axes of variation. Such axes that capture suites of interacting transcriptional changes during diversification improve our understanding of both global patterns in local adaptation and microdiversity, as well as evolutionary forces shaping algal cultivation.

Multi-trait diversification in marine diatoms in constant and warmed environments.

Phytoplankton are photosynthetic marine microbes that affect food webs, nutrient cycles and climate regulation. Their roles are determined by correlated phytoplankton functional traits including cell size, chlorophyll content and cellular composition. Here, we explore patterns of evolution in interrelated trait values and correlations. Because both chance events and natural selection contribute to phytoplankton trait evolution, we used population bottlenecks to diversify six genotypes of Thalassiosirid diatoms. We then evolved them as large populations in two environments. Interspecific variation and within-species evolution were visualized for nine traits and their correlations using reduced axes (a trait-scape). Our main findings are that shifts in trait values resulted in movement of evolving populations within the trait-scape in both environments, but were more frequent when large populations evolved in a novel environment. Which trait relationships evolved was population-specific, but greater departures from ancestral trait correlations were associated with lower population growth rates. There was no single master trait that could be used to understand multi-trait evolution. Instead, repeatable multi-trait evolution occurred along a major axis of variation defined by several diatom traits and trait relationships. Because trait-scapes capture changes in trait relationships and values together, they offer an insightful way to study multi-trait variation.

A unified theory for organic matter accumulation.

Organic matter constitutes a key reservoir in global elemental cycles. However, our understanding of the dynamics of organic matter and its accumulation remains incomplete. Seemingly disparate hypotheses have been proposed to explain organic matter accumulation: the slow degradation of intrinsically recalcitrant substrates, the depletion to concentrations that inhibit microbial consumption, and a dependency on the consumption capabilities of nearby microbial populations. Here, using a mechanistic model, we develop a theoretical framework that explains how organic matter predictably accumulates in natural environments due to biochemical, ecological, and environmental factors. Our framework subsumes the previous hypotheses. Changes in the microbial community or the environment can move a class of organic matter from a state of functional recalcitrance to a state of depletion by microbial consumers. The model explains the vertical profile of dissolved organic carbon in the ocean and connects microbial activity at subannual timescales to organic matter turnover at millennial timescales. The threshold behavior of the model implies that organic matter accumulation may respond nonlinearly to changes in temperature and other factors, providing hypotheses for the observed correlations between organic carbon reservoirs and temperature in past earth climates.

Microbial evolutionary strategies in a dynamic ocean.

Marine microbes form the base of ocean food webs and drive ocean biogeochemical cycling. Yet little is known about the ability of microbial populations to adapt as they are advected through changing conditions. Here, we investigated the interplay between physical and biological timescales using a model of adaptation and an eddy-resolving ocean circulation climate model. Two criteria were identified that relate the timing and nature of adaptation to the ratio of physical to biological timescales. Genetic adaptation was impeded in highly variable regimes by nongenetic modifications but was promoted in more stable environments. An evolutionary trade-off emerged where greater short-term nongenetic transgenerational effects (low-γ strategy) enabled rapid responses to environmental fluctuations but delayed genetic adaptation, while fewer short-term transgenerational effects (high-γ strategy) allowed faster genetic adaptation but inhibited short-term responses. Our results demonstrate that the selective pressures for organisms within a single water mass vary based on differences in generation timescales resulting in different evolutionary strategies being favored. Organisms that experience more variable environments should favor a low-γ strategy. Furthermore, faster cell division rates should be a key factor in genetic adaptation in a changing ocean. Understanding and quantifying the relationship between evolutionary and physical timescales is critical for robust predictions of future microbial dynamics.

Bacterial chemotaxis to saccharides is governed by a trade-off between sensing and uptake.

To swim up gradients of nutrients, E. coli senses nutrient concentrations within its periplasm. For small nutrient molecules, periplasmic concentrations typically match extracellular concentrations. However, this is not necessarily the case for saccharides, such as maltose, which are transported into the periplasm via a specific porin. Previous observations have shown that, under various conditions, E. coli limits maltoporin abundance so that, for extracellular micromolar concentrations of maltose, there are predicted to be only nanomolar concentrations of free maltose in the periplasm. Thus, in the micromolar regime, the total uptake of maltose from the external environment into the cytoplasm is limited not by the abundance of cytoplasmic transport proteins but by the abundance of maltoporins. Here, we present results from experiments and modeling suggesting that this porin-limited transport enables E. coli to sense micromolar gradients of maltose despite having a high-affinity ABC transport system that is saturated at these micromolar levels. We used microfluidic assays to study chemotaxis of E. coli in various gradients of maltose and methyl-aspartate and leveraged our experimental observations to develop a mechanistic transport-and-sensing chemotaxis model. Incorporating this model into agent-based simulations, we discover a trade-off between uptake and sensing: although high-affinity transport enables higher uptake rates at low nutrient concentrations, it severely limits the range of dynamic sensing. We thus propose that E. coli may limit periplasmic uptake to increase its chemotactic sensitivity, enabling it to use maltose as an environmental cue.

Mechanistic model of nutrient uptake explains dichotomy between marine oligotrophic and copiotrophic bacteria.

Marine bacterial diversity is immense and believed to be driven in part by trade-offs in metabolic strategies. Here we consider heterotrophs that rely on organic carbon as an energy source and present a molecular-level model of cell metabolism that explains the dichotomy between copiotrophs-which dominate in carbon-rich environments-and oligotrophs-which dominate in carbon-poor environments-as the consequence of trade-offs between nutrient transport systems. While prototypical copiotrophs, like Vibrios, possess numerous phosphotransferase systems (PTS), prototypical oligotrophs, such as SAR11, lack PTS and rely on ATP-binding cassette (ABC) transporters, which use binding proteins. We develop models of both transport systems and use them in proteome allocation problems to predict the optimal nutrient uptake and metabolic strategy as a function of carbon availability. We derive a Michaelis-Menten approximation of ABC transport, analytically demonstrating how the half-saturation concentration is a function of binding protein abundance. We predict that oligotrophs can attain nanomolar half-saturation concentrations using binding proteins with only micromolar dissociation constants and while closely matching transport and metabolic capacities. However, our model predicts that this requires large periplasms and that the slow diffusion of the binding proteins limits uptake. Thus, binding proteins are critical for oligotrophic survival yet severely constrain growth rates. We propose that this trade-off fundamentally shaped the divergent evolution of oligotrophs and copiotrophs.

DMSP synthesis genes distinguish two types of DMSP producer phenotypes.

Dimethylsulfoniopropionate (DMSP) is an important organic carbon and sulfur source in the surface ocean that fuels microbial activity and significantly impacts Earth's climate. After three decades of research, the cellular role(s) of DMSP and environmental drivers of production remain enigmatic. Recent work suggests that cellular DMSP concentrations, and changes in these concentrations in response to environmental stressors, define two major groups of DMSP producers: high DMSP producers that contain ≥ 50 mM intracellular DMSP and low DMSP producers that contain < 50 mM. Here we show that two recently described DMSP synthesis genes (DSYB and TpMT2) may differentiate these two DMSP phenotypes. A survey of prokaryotic and eukaryotic isolates found a significant correlation between the presence of DSYB and TpMT2 genes and previous measurements of high and low DMSP concentrations, respectively. Phylogenetic analysis demonstrated that DSYB and TpMT2 form two distinct clades. DSYB and TpMT2 were also found to be globally abundant in in situ surface communities, and their taxonomic annotations were similar to those observed for isolates. The strong correlation of the DSYB and TpMT2 synthesis genes with high and low producer phenotypes establishes a foundation for direct quantification of DMSP producers, enabling significantly improved predictions of DMSP in situ.

The evolution of trait correlations constrains phenotypic adaptation to high CO2 in a eukaryotic alga.

Microbes form the base of food webs and drive biogeochemical cycling. Predicting the effects of microbial evolution on global elemental cycles remains a significant challenge due to the sheer number of interacting environmental and trait combinations. Here, we present an approach for integrating multivariate trait data into a predictive model of trait evolution. We investigated the outcome of thousands of possible adaptive walks parameterized using empirical evolution data from the alga Chlamydomonas exposed to high CO2. We found that the direction of historical bias (existing trait correlations) influenced both the rate of adaptation and the evolved phenotypes (trait combinations). Critically, we use fitness landscapes derived directly from empirical trait values to capture known evolutionary phenomena. This work demonstrates that ecological models need to represent both changes in traits and changes in the correlation between traits in order to accurately capture phytoplankton evolution and predict future shifts in elemental cycling.

Understanding water and energy fluxes in the Amazonia: Lessons from an observation-model intercomparison.

Tropical forests are an important part of global water and energy cycles, but the mechanisms that drive seasonality of their land-atmosphere exchanges have proven challenging to capture in models. Here, we (1) report the seasonality of fluxes of latent heat (LE), sensible heat (H), and outgoing short and longwave radiation at four diverse tropical forest sites across Amazonia-along the equator from the Caxiuanã and Tapajós National Forests in the eastern Amazon to a forest near Manaus, and from the equatorial zone to the southern forest in Reserva Jaru; (2) investigate how vegetation and climate influence these fluxes; and (3) evaluate land surface model performance by comparing simulations to observations. We found that previously identified failure of models to capture observed dry-season increases in evapotranspiration (ET) was associated with model overestimations of (1) magnitude and seasonality of Bowen ratios (relative to aseasonal observations in which sensible was only 20%-30% of the latent heat flux) indicating model exaggerated water limitation, (2) canopy emissivity and reflectance (albedo was only 10%-15% of incoming solar radiation, compared to 0.15%-0.22% simulated), and (3) vegetation temperatures (due to underestimation of dry-season ET and associated cooling). These partially compensating model-observation discrepancies (e.g., higher temperatures expected from excess Bowen ratios were partially ameliorated by brighter leaves and more interception/evaporation) significantly biased seasonal model estimates of net radiation (Rn ), the key driver of water and energy fluxes (LE ~ 0.6 Rn and H ~ 0.15 Rn ), though these biases varied among sites and models. A better representation of energy-related parameters associated with dynamic phenology (e.g., leaf optical properties, canopy interception, and skin temperature) could improve simulations and benchmarking of current vegetation-atmosphere exchange and reduce uncertainty of regional and global biogeochemical models.

Evidence for contrasting roles of dimethylsulfoniopropionate production in Emiliania huxleyi and Thalassiosira oceanica.

Dimethylsulfoniopropionate (DMSP) is a globally abundant marine metabolite and a significant source of organic carbon and sulfur for marine microbial ecosystems with the potential to influence climate regulation. However, the physiological function of DMSP has remained enigmatic for >30 yr. Recent insight suggests that there are different physiological roles for DMSP based on the cellular DMSP concentrations in producers. Differential production of DMSP was tested with multiple physiological experiments that altered nitrate availability, salinity and temperature to create stressed growth and target different metabolic conditions in Emiliania huxleyi, a high DMSP producer and Thalassiosira oceanica, a low DMSP producer. Emiliania huxleyi intracellular DMSP did not respond to metabolically imbalanced conditions, while Thalassiosira oceanica intracellular DMSP was significantly correlated to stressed growth rate across all conditions tested and exhibited a plastic response on a timescale of hours in nonsteady-state. The previous assumption that proposed DMSP mechanism(s) can be universally applied to all producers is shown to be unlikely. Rather, two distinct ecological roles for DMSP likely exist that differ by producer type, where: (1) the primary role of DMSP in high producers is a constitutive compatible solute; and (2) DMSP production in low producers is a finely tuned stress response.

Ecosystem heterogeneity determines the ecological resilience of the Amazon to climate change.

Amazon forests, which store ∼ 50% of tropical forest carbon and play a vital role in global water, energy, and carbon cycling, are predicted to experience both longer and more intense dry seasons by the end of the 21st century. However, the climate sensitivity of this ecosystem remains uncertain: several studies have predicted large-scale die-back of the Amazon, whereas several more recent studies predict that the biome will remain largely intact. Combining remote-sensing and ground-based observations with a size- and age-structured terrestrial ecosystem model, we explore the sensitivity and ecological resilience of these forests to changes in climate. We demonstrate that water stress operating at the scale of individual plants, combined with spatial variation in soil texture, explains observed patterns of variation in ecosystem biomass, composition, and dynamics across the region, and strongly influences the ecosystem's resilience to changes in dry season length. Specifically, our analysis suggests that in contrast to existing predictions of either stability or catastrophic biomass loss, the Amazon forest's response to a drying regional climate is likely to be an immediate, graded, heterogeneous transition from high-biomass moist forests to transitional dry forests and woody savannah-like states. Fire, logging, and other anthropogenic disturbances may, however, exacerbate these climate change-induced ecosystem transitions.

Ecosystem heterogeneity and diversity mitigate Amazon forest resilience to frequent extreme droughts.

The impact of increases in drought frequency on the Amazon forest's composition, structure and functioning remain uncertain. We used a process- and individual-based ecosystem model (ED2) to quantify the forest's vulnerability to increased drought recurrence. We generated meteorologically realistic, drier-than-observed rainfall scenarios for two Amazon forest sites, Paracou (wetter) and Tapajós (drier), to evaluate the impacts of more frequent droughts on forest biomass, structure and composition. The wet site was insensitive to the tested scenarios, whereas at the dry site biomass declined when average rainfall reduction exceeded 15%, due to high mortality of large-sized evergreen trees. Biomass losses persisted when year-long drought recurrence was shorter than 2-7 yr, depending upon soil texture and leaf phenology. From the site-level scenario results, we developed regionally applicable metrics to quantify the Amazon forest's climatological proximity to rainfall regimes likely to cause biomass loss > 20% in 50 yr according to ED2 predictions. Nearly 25% (1.8 million km2 ) of the Amazon forests could experience frequent droughts and biomass loss if mean annual rainfall or interannual variability changed by 2σ. At least 10% of the high-emission climate projections (CMIP5/RCP8.5 models) predict critically dry regimes over 25% of the Amazon forest area by 2100.

Do dynamic global vegetation models capture the seasonality of carbon fluxes in the Amazon basin? A data-model intercomparison.

To predict forest response to long-term climate change with high confidence requires that dynamic global vegetation models (DGVMs) be successfully tested against ecosystem response to short-term variations in environmental drivers, including regular seasonal patterns. Here, we used an integrated dataset from four forests in the Brasil flux network, spanning a range of dry-season intensities and lengths, to determine how well four state-of-the-art models (IBIS, ED2, JULES, and CLM3.5) simulated the seasonality of carbon exchanges in Amazonian tropical forests. We found that most DGVMs poorly represented the annual cycle of gross primary productivity (GPP), of photosynthetic capacity (Pc), and of other fluxes and pools. Models simulated consistent dry-season declines in GPP in the equatorial Amazon (Manaus K34, Santarem K67, and Caxiuanã CAX); a contrast to observed GPP increases. Model simulated dry-season GPP reductions were driven by an external environmental factor, 'soil water stress' and consequently by a constant or decreasing photosynthetic infrastructure (Pc), while observed dry-season GPP resulted from a combination of internal biological (leaf-flush and abscission and increased Pc) and environmental (incoming radiation) causes. Moreover, we found models generally overestimated observed seasonal net ecosystem exchange (NEE) and respiration (Re ) at equatorial locations. In contrast, a southern Amazon forest (Jarú RJA) exhibited dry-season declines in GPP and Re consistent with most DGVMs simulations. While water limitation was represented in models and the primary driver of seasonal photosynthesis in southern Amazonia, changes in internal biophysical processes, light-harvesting adaptations (e.g., variations in leaf area index (LAI) and increasing leaf-level assimilation rate related to leaf demography), and allocation lags between leaf and wood, dominated equatorial Amazon carbon flux dynamics and were deficient or absent from current model formulations. Correctly simulating flux seasonality at tropical forests requires a greater understanding and the incorporation of internal biophysical mechanisms in future model developments.

Terrestrial and marine perspectives on modeling organic matter degradation pathways.

Organic matter (OM) plays a major role in both terrestrial and oceanic biogeochemical cycles. The amount of carbon stored in these systems is far greater than that of carbon dioxide (CO2 ) in the atmosphere, and annual fluxes of CO2 from these pools to the atmosphere exceed those from fossil fuel combustion. Understanding the processes that determine the fate of detrital material is important for predicting the effects that climate change will have on feedbacks to the global carbon cycle. However, Earth System Models (ESMs) typically utilize very simple formulations of processes affecting the mineralization and storage of detrital OM. Recent changes in our view of the nature of this material and the factors controlling its transformation have yet to find their way into models. In this review, we highlight the current understanding of the role and cycling of detrital OM in terrestrial and marine systems and examine how this pool of material is represented in ESMs. We include a discussion of the different mineralization pathways available as organic matter moves from soils, through inland waters to coastal systems and ultimately into open ocean environments. We argue that there is strong commonality between aspects of OM transformation in both terrestrial and marine systems and that our respective scientific communities would benefit from closer collaboration.

Enhancement of phytoplankton chlorophyll by submesoscale frontal dynamics in the North Pacific Subtropical Gyre.

Subtropical gyres contribute significantly to global ocean productivity. As the climate warms, the strength of these gyres as a biological carbon pump is predicted to diminish due to increased stratification and depleted surface nutrients. We present results suggesting that the impact of submesoscale physics on phytoplankton in the oligotrophic ocean is substantial and may either compensate or exacerbate future changes in carbon cycling. A new statistical tool was developed to quantify surface patchiness from sea surface temperatures. Chlorophyll concentrations in the North Pacific Subtropical Gyre were shown to be enhanced by submesoscale frontal dynamics with an average increase of 38% (maximum of 83%) during late winter. The magnitude of this enhancement is comparable to the observed decline in chlorophyll due to a warming of ~1.1°C. These results highlight the need for an improved understanding of fine-scale physical variability in order to predict the response of marine ecosystems to projected climate changes.

The fate of Amazonian ecosystems over the coming century arising from changes in climate, atmospheric CO2, and land use.

There is considerable interest in understanding the fate of the Amazon over the coming century in the face of climate change, rising atmospheric CO2 levels, ongoing land transformation, and changing fire regimes within the region. In this analysis, we explore the fate of Amazonian ecosystems under the combined impact of these four environmental forcings using three terrestrial biosphere models (ED2, IBIS, and JULES) forced by three bias-corrected IPCC AR4 climate projections (PCM1, CCSM3, and HadCM3) under two land-use change scenarios. We assess the relative roles of climate change, CO2 fertilization, land-use change, and fire in driving the projected changes in Amazonian biomass and forest extent. Our results indicate that the impacts of climate change are primarily determined by the direction and severity of projected changes in regional precipitation: under the driest climate projection, climate change alone is predicted to reduce Amazonian forest cover by an average of 14%. However, the models predict that CO2 fertilization will enhance vegetation productivity and alleviate climate-induced increases in plant water stress, and, as a result, sustain high biomass forests, even under the driest climate scenario. Land-use change and climate-driven changes in fire frequency are predicted to cause additional aboveground biomass loss and reductions in forest extent. The relative impact of land use and fire dynamics compared to climate and CO2 impacts varies considerably, depending on both the climate and land-use scenario, and on the terrestrial biosphere model used, highlighting the importance of improved quantitative understanding of all four factors - climate change, CO2 fertilization effects, fire, and land use - to the fate of the Amazon over the coming century.

Direct observations of microbial community succession on sinking marine particles.

Microbial community dynamics on sinking particles control the amount of carbon that reaches the deep ocean and the length of time that carbon is stored, with potentially profound impacts on Earth's climate. A mechanistic understanding of the controls on sinking particle distributions has been hindered by limited depth- and time-resolved sampling and methods that cannot distinguish individual particles. Here, we analyze microbial communities on nearly 400 individual sinking particles in conjunction with more conventional composite particle samples to determine how particle colonization and community assembly might control carbon sequestration in the deep ocean. We observed community succession with corresponding changes in microbial metabolic potential on the larger sinking particles transporting a significant fraction of carbon to the deep sea. Microbial community richness decreased as particles aged and sank; however, richness increased with particle size and the attenuation of carbon export. This suggests that the theory of island biogeography applies to sinking marine particles. Changes in POC flux attenuation with time and microbial community composition with depth were reproduced in a mechanistic ecosystem model that reflected a range of POC labilities and microbial growth rates. Our results highlight microbial community dynamics and processes on individual sinking particles, the isolation of which is necessary to improve mechanistic models of ocean carbon uptake.

Confronting model predictions of carbon fluxes with measurements of Amazon forests subjected to experimental drought.

Considerable uncertainty surrounds the fate of Amazon rainforests in response to climate change. Here, carbon (C) flux predictions of five terrestrial biosphere models (Community Land Model version 3.5 (CLM3.5), Ecosystem Demography model version 2.1 (ED2), Integrated BIosphere Simulator version 2.6.4 (IBIS), Joint UK Land Environment Simulator version 2.1 (JULES) and Simple Biosphere model version 3 (SiB3)) and a hydrodynamic terrestrial ecosystem model (the Soil-Plant-Atmosphere (SPA) model) were evaluated against measurements from two large-scale Amazon drought experiments. Model predictions agreed with the observed C fluxes in the control plots of both experiments, but poorly replicated the responses to the drought treatments. Most notably, with the exception of ED2, the models predicted negligible reductions in aboveground biomass in response to the drought treatments, which was in contrast to an observed c. 20% reduction at both sites. For ED2, the timing of the decline in aboveground biomass was accurate, but the magnitude was too high for one site and too low for the other. Three key findings indicate critical areas for future research and model development. First, the models predicted declines in autotrophic respiration under prolonged drought in contrast to measured increases at one of the sites. Secondly, models lacking a phenological response to drought introduced bias in the sensitivity of canopy productivity and respiration to drought. Thirdly, the phenomenological water-stress functions used by the terrestrial biosphere models to represent the effects of soil moisture on stomatal conductance yielded unrealistic diurnal and seasonal responses to drought.

Mechanistic constraints on the trade-off between photosynthesis and respiration in response to warming.

Phytoplankton are responsible for half of all oxygen production and drive the ocean carbon cycle. Metabolic theory predicts that increasing global temperatures will cause phytoplankton to become more heterotrophic and smaller. Here, we uncover the metabolic trade-offs between cellular space, energy, and stress management driving phytoplankton thermal acclimation and how these might be overcome through evolutionary adaptation. We show that the observed relationships between traits such as chlorophyll, lipid content, C:N, and size can be predicted on the basis of the metabolic demands of the cell, the thermal dependency of transporters, and changes in membrane lipids. We suggest that many of the observed relationships are not fixed physiological constraints but rather can be altered through adaptation. For example, the evolution of lipid metabolism can favor larger cells with higher lipid content to mitigate oxidative stress. These results have implications for rates of carbon sequestration and export in a warmer ocean.

Identification of microbial metabolic functional guilds from large genomic datasets.

Heterotrophic microbes play an important role in the Earth System as key drivers of major biogeochemical cycles. Specifically, the consumption rate of organic matter is set by the interaction between diverse microbial communities and the chemical and physical environment in which they reside. Modeling these dynamics requires reducing the complexity of microbial communities and linking directly with biogeochemical functions. Microbial metabolic functional guilds provide one approach for reducing microbial complexity and incorporating microbial biogeochemical functions into models. However, we lack a way to identify these guilds. In this study, we present a method for defining metabolic functional guilds from annotated genomes, which are derived from both uncultured and cultured organisms. This method utilizes an Aspect Bernoulli (AB) model and was tested on three large genomic datasets with 1,733-3,840 genomes each. Ecologically relevant microbial metabolic functional guilds were identified including guilds related to DMSP degradation, dissimilatory nitrate reduction to ammonia, and motile copiotrophy. This method presents a way to generate hypotheses about functions co-occurring within individual microbes without relying on cultured representatives. Applying the concept of metabolic functional guilds to environmental samples will provide new insight into the role that heterotrophic microbial communities play in setting rates of carbon cycling.