Original Research ARTICLE
Carbon Availability Modifies Temperature Responses of Heterotrophic Microbial Respiration, Carbon Uptake Affinity, and Stable Carbon Isotope Discrimination
- Department of Ecology and Evolutionary Biology, Kansas Biological Survey, University of Kansas, Lawrence, KS, USA
Microbial transformations of organic carbon (OC) generate a large flux of CO2 into the atmosphere and influence the C balance of terrestrial and aquatic ecosystems. Yet, inherent heterogeneity in natural environments precludes direct quantification of multiple microbial C fluxes that underlie CO2 production. Here we used a continuous flow bioreactor coupled with a stable C isotope analyzer to determine the effects of temperature and C availability (cellobiose concentration) on C fluxes and 13C discrimination of a microbial population growing at steady-state in a homogeneous, well-mixed environment. We estimated C uptake affinity and C use efficiency (CUE) to characterize the physiological responses of microbes to changing environmental conditions. Temperature increased biomass-C specific respiration rate and C uptake affinity at lower C availability, but did not influence those parameters at higher C availability. CUE decreased non-linearly with increasing temperature. The non-linear, negative relationship between CUE and temperature was more pronounced under lower C availability than under relatively high C availability. We observed stable isotope fractionation between C substrate and microbial biomass C (7~12‰ depletion), and between microbial biomass and respired CO2 (4~10‰ depletion). Microbial discrimination against 13C-containing cellobiose during C uptake was influenced by temperature and C availability, while discrimination during respiration was only influenced by C availability. Shifts in C uptake affinity with temperature and C availability may have modified uptake-induced 13C fractionation. By stressing the importance of C availability on temperature responses of microbial C fluxes, C uptake affinity, CUE, and isotopic fractionation, this study contributes to a fundamental understanding of C flow through microbes. This will help guide parameterization of microbial responses to varying temperature and C availability within Earth-system models.
Heterotrophic microorganisms break down organic carbon (OC) into assimilable molecules, which are subsequently incorporated into biomass or returned to the environment as exudates or respired CO2. In doing so, microbes influence the quantity and biochemical composition of OC in terrestrial and aquatic systems as well as atmospheric CO2 concentration. Because of the important role they play in regulating OC in terrestrial and aquatic systems and atmospheric CO2, ecosystem scientists have explicitly incorporated microbial dynamics into Earth system models (Allison et al., 2010; Wieder et al., 2013; Hagerty et al., 2014; Tang and Riley, 2015). However, it remains unclear what values should be used for microbial parameters [e.g., C use efficiency (CUE), biomass-C specific respiration] and how these parameters change with environmental conditions (Luo et al., 2016). This is because heterogeneity in natural environments, especially in soils, makes it nearly impossible to empirically determine many of these parameters in situ (Werth and Kuzyakov, 2010; Billings and Ballantyne, 2013; Billings et al., 2015).
Assessing microbial OC transformations in controlled, experimental systems can help us overcome some of the difficulties ecosystem scientists face when using natural samples to investigate microbial C flows (Billings et al., 2015). For example, chemostats, continuous flow bioreactors for the cultivation of microorganisms, have been widely employed to investigate fundamental microbial physiology during OC transformations by removing the influence of natural heterogeneity (Rhee and Gotham, 1981; Cajal-Medrano and Maske, 2005; Cotner et al., 2006; Chrzanowski and Grover, 2008; Ferenci, 2008). Recently, Lehmeier et al. (2016) used a coupled chemostat and C isotope analyzer to improve our understanding of C flows from substrate through microbial biomass and into respired CO2, which is of great interest at multiple scales and levels of organization, from microbes to ecosystems (Blair et al., 1985; Hagström et al., 1988; Henn and Chapela, 2000; Šantrůčková et al., 2000; Henn et al., 2002; Werth and Kuzyakov, 2010; Brüggemann et al., 2011; Dijkstra et al., 2011a,b; Penger et al., 2014). In particular, quantifying the effects of temperature on biomass-C specific respiration rate, CUE, and microbial 13C discrimination for a ubiquitous microbe (Lehmeier et al., 2016) provides a mechanistic basis and justification for the often assumed temperature dependence of CUE in Earth-system models.
Temperature responses of microbial OC transformations in Lehmeier et al. (2016) were quantified using a single C substrate, cellobiose, at a single concentration, 10 mM. However, it remains unknown how the temperature responses of microbial respiration, CUE, and 13C discrimination may be influenced by C availability which varies among ecosystems, vertically and laterally throughout soil profiles, and as microbes’ substrate landscapes shift over time (Kuzyakov, 2010; Richter and Markewitz, 2013). Furthermore, the absolute and relative abundance of readily assimilable C can also vary with temperature (Bárta et al., 2013; Lehmeier et al., 2013; Min et al., 2014), and thus can influence the temperature responses of microbial OC transformations (Rhee and Gotham, 1981; Pomeroy and Wiebe, 2001; Hall and Cotner, 2007; Manzoni et al., 2012; Billings and Ballantyne, 2013; Frey et al., 2013; Xu et al., 2014).
When C availability decreases, microorganisms can adjust their ability to take up C as a compensatory measure (Death et al., 1993; Williams et al., 1994; Ferenci, 1999; Gresham and Hong, 2015). The ability to take up C, or C uptake affinity, is often defined as the ratio of the maximum uptake rate (Vmax) to the half-saturation concentration of C (km) of the Michaelis-Menten function (Healey, 1980; Aksnes and Egge, 1991; Button, 1993; Reay et al., 1999). Because C uptake affinity defined this way is equivalent to the slope of the Michaelis–Menten function at a C concentration of zero, it reflects the sensitivity of the uptake response to changing C availability at low absolute concentrations. This suggests that microbes may allocate resources to a more energetically costly, high affinity C uptake strategy at relatively low C availability (Patzer and Hantke, 1998), resulting in enhanced respiratory costs. Furthermore, changing C uptake affinity may alter 13C discrimination during uptake and subsequently influence δ13C of remaining OC pools, as well as downstream pools such as biomass and respired CO2 (Lehmeier et al., 2016). Understanding any such C uptake affinity-induced shift in fractionation would be critical to correctly interpret δ13C of terrestrial C pools and fluxes.
We grew Pseudomonas fluorescens, a cosmopolitan Gram-negative bacterium, in a chemostat system coupled to a 13CO2/12CO2 analyzer and quantified how the sizes of C pools in the system and their C stable isotope ratios varied with temperature and C availability. Combining the measurements with a dynamic model of substrate-C and biomass-C, we estimated how C uptake and C uptake affinity varied with changing temperature (11.5–25.5°C) and C availability (1 and 20 mM cellobiose). The cellobiose concentrations of 1 and 20 mM yielded medium C:N of 1 and 20, respectively, to span the observed range of C:N of bacterial biomass (5–10) reported in Nagata (1986) and Cleveland and Liptzin (2007). Our primary goal was to quantify how changing C substrate availability alters the temperature dependence of C flow through a common microbe and associated C stable isotope fractionations. As such, we asked: (1) how will temperature and C availability influence biomass-C specific respiration, C uptake affinity, and CUE during OC transformations? (2) how will changes in these microbial parameters with temperature and C availability, if any, modify 13C discrimination among substrate, biomass and respired CO2?
Materials and Methods
The chemostat system coupled to a continuous flow 13CO2/12CO2 gas analyzer employed in this study is described in full detail in Lehmeier et al. (2016). Here we summarize the protocol of running a chemostat experiment and expand Lehmeier et al.’s approach by providing estimates for C uptake affinity and a simple model for CUE.
Microorganisms and Nutrient Solution With Two Contrasting C Concentrations
We used Pseudomonas fluorescens (ATCC ®13525, Carolina Biological Supply, USA), because it is one of the most common bacteria inhabiting both soil and water. As growth medium for P. fluorescens, we prepared nutrient solution modified from Abraham et al. (1998), containing 10 mM NH4Cl, 1.6 mM KNO3, 2.6 mM K2HPO4, 1.0 mM KH2PO4, 0.8 mM MgSO4, 0.2 mM CaCl2, 0.1 mM CuCl2, 0.04 mM FeSO4, 0.03 mM MnCl2 and 0.02 mM ZnSO4. After adjustment to pH 6.5, the nutrient solution was autoclaved and stored under UV light. We added cellobiose (Sigma–Aldrich, USA) to the sterile nutrient solution as the only organic C source, inoculated it with P. fluorescens, and grew cultures at 10°C. These pre-cultures of P. fluorescens served as inoculum for the chemostat bioreactor. The concentration of cellobiose in the nutrient solution was adjusted to either 1 or 20 mM, which corresponded to resource C to N ratios of 1 and 20, respectively. Cellobiose δ13C values ranged from –24.4 to –24.9‰, depending on the lot from which the supply was issued. For each individual chemostat run, we used cellobiose derived from only one lot to ensure a constant δ13C value of microbial substrate for each run (Supplementary Table 1).
Continuous Flow Chemostat System
The chemostat system consisted of two 1.9 L Mason jars (Ball, USA), which were connected via flexible tubing (Masterflex Tygon E-LFL tubing, Cole-Parmer, USA) and situated on stir plates in separate incubators. One jar functioned as a bioreactor, into which we poured 1 L of fresh nutrient solution and injected inoculum (5 mL of P. fluorescens pre-culture) at the start of each individual chemostat run; the other jar functioned as a reservoir for fresh nutrient solution and was regularly refilled during a chemostat run.
In the continuous culture mode, a peristaltic pump supplied fresh nutrient solution from the reservoir to the bioreactor at a constant rate, providing new resources to the growing microbial population. The pump simultaneously removed well-mixed bioreactor medium (including microbial cells) at the same rate as the nutrient solution supply rate via a waste line, such that the total volume of liquid in the bioreactor remained constant. The dilution rate of the chemostat bioreactor averaged 0.128 h–1 (Figure 1A). By definition, the dilution rate equals the growth rate of the microbial population in the chemostat reactor under steady-state conditions (Hoskisson and Hobbs, 2005; Bull, 2010). Steady state in this sense means that the microbial population in the bioreactor reaches a state of constant size and activity under the prevailing environmental conditions, and that the number of cells leaving the bioreactor via dilution is counterbalanced by cell division.
FIGURE 1. Biomass-C specific rates of growth (A), respiration (B), and uptake (C) with varying temperature for P. fluorescens grown in chemostats with cellobiose concentrations in the nutrient solution of either 1 mM (open circles) or 20 mM (closed circles). We imposed a regression line on data points when there was a significant temperature effect on rates (see Supplementary Table 2). Error bars in (B) represent ± 1SE.
Analogous to the continuous flow of liquid, gas also flowed continuously through the bioreactor. A 13CO2/12CO2 analyzer (G2101-i, Picarro, USA) removed headspace air from the bioreactor at a rate of 25 mL min–1 on average and continuously measured the concentration and the δ13C of headspace CO2. The air volume removed was instantaneously replaced with CO2-free air, bubbled into the nutrient solution at the same rate that the 13CO2/12CO2 analyzer withdrew headspace air via a mass flow controller (MC-50SCCM, Alicat Scientific, USA). In this way, bioreactor headspace air pressure remained constant and the CO2-free air provided the microbial population with oxygen.
Immediately after the inoculation and closure of the chemostat bioreactor, CO2 concentration in the headspace air represented lab air conditions (~400 ppm CO2). Soon after closure of the bioreactor, the CO2 concentration approached zero because of the bubbling of CO2-free air into the bioreactor and little microbial activity. Subsequently, the size of the microbial population and hence the rate of respiratory CO2 production increased, and after the switch to the continuous culture mode, it achieved steady state conditions. Prior to the steady state, the 13CO2/12CO2 measurement of the headspace air did not solely reflect microbial respiration, because the bioreactor medium acted as a sink for respired CO2 (in the form of H2CO3 and ; Stumm and Morgan, 1981; Lehmeier et al., 2016). When the microbial population reached steady state, both the CO2 production rate and the δ13C of the CO2 became constant, indicating that chemical and isotopic equilibria of the carbonate system were achieved. However, even in the steady-state phase, not all of the respired CO2 entering the bioreactor medium was released into the bioreactor headspace. This is because some respired CO2 entered the carbonate system and was removed via the waste line due to dilution. We estimated that not more than 6~8% of respired CO2 was lost from the bioreactor via dilution, hence, we may have underestimated microbial respiration rates by about 6–8%. Importantly, as both the supply rates of fresh medium and the steady-state CO2 concentrations in the bioreactor headspace (the two most important factors determining the sink strength of the supplied fresh medium) were very similar across all chemostat runs, the degree of underestimating microbial respiration rate was very similar across chemostat runs. Therefore, we considered the headspace CO2 representing the concentration and δ13C of respired CO2 (for more detailed information, please refer to the Supplementary Material in Lehmeier et al., 2016). We thus used measurements of bioreactor headspace air when the microbial population was at steady state to quantify microbial respiration rates and δ13C of respired CO2, allowing us to draw conclusions about microbial behavior without concern for interactions between headspace CO2 and inorganic C pools in the bioreactor medium.
For each cellobiose concentration, chemostat experiments were conducted at six different bioreactor temperatures, randomly ordered: at 1 mM of cellobiose, reactor temperature was maintained at 11.8, 13, 15.5, 17, 21, or 25.5°C; at 20 mM of cellobiose, it was 11.5, 13, 14.7, 17.5, 20.8, or 25.5°C. This experimental temperature range covers an ecologically relevant range of temperatures in many ecosystems. Below 11.5°C, microorganisms’ steady-state growth was unsustainable at the dilution rate of 0.128 h–1. This is consistent with findings in Höfle (1979), which reported that the maximum growth rate of P. fluorescens at 10°C is 0.13 h–1. We conducted one repeated chemostat run at 17°C and 1 mM of cellobiose. We acknowledge that one replication is not enough to quantify an overall uncertainty of the measurements, but similar values of microbial OC transformations for these replicated runs, in conjunction with clear temperature-related trends (see Results), suggest that our measurements were reliable (Figures 1–4).
Quantification of Microbial Biomass and Respiration at Steady-State
At steady state, we collected and filtered the bioreactor medium from the waste line on dry, pre-weighed 0.22 μm polyethersulfone filters (Pall, USA). We dried the wet filters at 75°C for 48 h in an oven and reweighed them to determine dry microbial biomass. Carbon and N in microbial biomass were quantified using an elemental analyzer coupled to an isotope-ratio mass spectrometer, either at the Keck Paleoenvironmental and Environmental Stable Isotope Laboratory at the University of Kansas or at the Stable Isotope Laboratory at the University of Arkansas. In both laboratories, the standard deviation of repeated lab standard measurements was a maximum of 0.2‰.
We calculated microbial respiration rates by multiplying the average reactor headspace CO2 concentration measured over 5 h at steady state with the molar air flow rate through the bioreactor. After steady state CO2 measurements for each chemostat run, we measured a laboratory standard gas, which had a CO2 concentration of 1015 ppm and a δ13C of –48.9‰ as calibrated against secondary CO2 standards from KU’s Keck Paleoenvironmental and Environmental Stable Isotope Laboratory. This served to check the concentrations and to correct the isotopic signatures of bioreactor headspace CO2 measurements.
Carbon Dynamics in the Chemostat Bioreactor
In a chemostat bioreactor provided with cellobiose as the sole OC substrate there are two OC pools, cellobiose-C and biomass-C, if we assume exudation is negligible (Blair et al., 1985; Hagström et al., 1988). The temporal dynamics of the two OC pools can be described with two ordinary differential equations (Smith and Waltman, 1995).
where S is cellobiose-C in the chemostat bioreactor, I is cellobiose-C input to the chemostat bioreactor, d is the dilution rate (which equals the biomass-C-specific population growth rate in a steady-state chemostat; Hoskisson and Hobbs, 2005; Bull, 2010), B is biomass-C in the chemostat reactor, f(S) is the function describing biomass-C specific uptake rate of cellobiose-C, ρ is the biomass-C specific respiration rate, 𝜀 is the biomass-C specific exudation rate, and t is time. At steady state, denoted by ∧, biomass-C and cellobiose-C concentrations are both constant . This means that the biomass-C specific uptake rate, f(S), equals the sum of d, ρ and 𝜀. If we assume 𝜀 is negligible, f(S) is estimated as the sum of d and ρ. The steady state conditions in the chemostat allow biomass-C specific rates to be computed by dividing gross C fluxes by C pool sizes. For example, ρ was estimated by dividing the total rate of respiration of the bioreactor by microbial biomass-C at steady state, which we obtained via elemental analysis of dried biomass.
Estimates of C Uptake Affinity and CUE
Although the steady state measurements of d and ρ allow us to estimate f(S), this alone does not provide the physiological basis for a change in uptake fluxes. However, if we employ a standard assumption of f(S) saturating as a function of C availability (Button, 1993; Ferenci, 1999), we can begin to probe temperature-dependent physiological shifts. We can express f(S), the function describing the biomass-C specific uptake rate of cellobiose-C (see equation 1), with the standard Michaelis-Menten formulation for cellobiose-C uptake as
This approach allows us to relate changes in C uptake physiology, namely in Vmax and km, to changing temperature, such that a change in C uptake physiology would be reflected in a change in either or both with temperature and C availability.
Assuming that 𝜀 is negligible and combining equations 2 and 3, we can express steady-state concentrations of biomass as
We can rearrange equation 4 to consolidate all our measurements and pre-specified quantities on the right hand side:
For convenience, we replace the right side of equation 5 with A and rearrange terms to arrive at
Here, 1/A represents the slope of Vmax over km and equals . If we assume Vmax is greatly larger than the sum of d and ρ (Schulze and Lipe, 1964; Button, 1985; Kovárová-Kovar and Egli, 1998), then 1/A approaches C uptake affinity, defined as the ratio of Vmax over km (Button, 1993; Nedwell, 1999; Reay et al., 1999). Thus, hereafter we refer to 1/A as estimated cellobiose-C uptake affinity for P. fluorescens at steady state. Our estimates of C uptake affinity are valuable because they provide insight to how P. fluorescens responds physiologically to changing temperature and C availability, without independently quantifying the kinetic parameters Vmax and km.
We estimated CUE by substituting the apparently linear relationship between biomass-C specific respiration rate and temperature (see Statistics below) via equation 7:
where i corresponds to either 1 or 20 mM cellobiose concentrations in the nutrient medium, αi is the y-intercept of the regression line relating biomass-C specific respiration rate to temperature at each C availability, and βi is the slope of the line at each C availability. Using errors obtained from the modeled temperature-dependence of biomass-C specific respiration rates, we calculated 95% confidence envelopes surrounding these CUE estimates. Because the dilution rate of the chemostat bioreactor was relatively constant at 0.128 h–1, changes in CUE estimates were driven by temperature-dependent changes in biomass-C specific respiration rates.
We used general linear models to relate biomass-C specific respiration rate and C uptake affinity to temperature and cellobiose concentration in the nutrient solution. Both models consider temperature as a continuous covariate, and cellobiose concentration as a categorical predictor. Because biomass-C specific respiration rates were a key feature of our estimates of biomass-C specific uptake rates and CUE (see Carbon Dynamics in the Chemostat Bioreactor), we did not fit separate models for these two responses to test the effects of temperature and cellobiose concentration. For the C stable isotope data, we treated temperature as a continuous covariate, and cellobiose concentration and C pool (biomass and CO2) as categorical predictors. To correct any effect of varying δ13C of cellobiose across experiments on δ13C of both biomass and CO2, we used the differences between δ13C of cellobiose and δ13C of biomass, and between δ13C of cellobiose and δ13C of CO2, as dependent variables. We used the Akaike Information Criterion (AIC) to inform model selection for each of the three response variables. All statistical analyses were performed with R (R version 3.1.2. R development core team).
Microbial Biomass-C Specific Rates: Respiration, Growth, and Uptake
Increasing temperature significantly enhanced biomass-C specific respiration at 1 mM of cellobiose (p < 0.001), while temperature did not significantly influence biomass-C specific respiration at 20 mM of cellobiose (Figure 1B; Supplementary Table 2). Because our estimates of biomass-C specific uptake were computed by summing biomass-C specific rates of growth and respiration (see Carbon Dynamics in the Chemostat Bioreactor), and biomass-C specific growth rate was relatively constant across all chemostat runs (Figure 1A), temperature responses of biomass-C specific uptake rates were similar to those of biomass-C specific respiration at each C availability (Figure 1C). Biomass-C specific uptake rates at both C availabilities were similar at lower temperatures, but diverged with increasing temperature. At 25.5°C, biomass-C specific uptake rate of the microbes at 1 mM cellobiose was almost twice as high as the biomass-C specific uptake rates at 20 mM cellobiose. We calculated that P. fluorescens took up, on average across temperatures, 17.4 ± 1.9 and 1.3 ± 0.1% of the cellobiose-C substrate supplied per hour at 1 and 20 mM, respectively.
Estimates of C Uptake Affinity and CUE
The interactive effect of temperature and C availability on biomass-C specific respiration rate influenced the estimates of C uptake affinity and CUE. As temperature increased from 11.5 to 25.5°C, C uptake affinity increased from 1.31 to 4.01 mL h–1 mg–1 C at 1 mM of cellobiose, but temperature did not influence C uptake affinity at 20 mM of cellobiose (Figure 2; Supplementary Table 2). Carbon uptake affinity was always higher at 1 mM of cellobiose than at 20 mM, and exhibited no convergence at the lowest temperature in contrast to the biomass-C specific uptake and respiration rates.
FIGURE 2. Temperature responses of C uptake affinity, the slope of Vmax over km (see Estimates of C Uptake Affinity and CUE), estimated from the two-pool model of P. fluorescens grown in chemostats with cellobiose concentrations in the nutrient solution of either 1 mM (open circles) or 20 mM (closed circles). Model parameters are available in Supplementary Table 2.
Microbial CUE ranged between 0.3 and 0.8, with similar maximum values for both C availabilities (Figure 3). Within the experimental temperature range, microbial CUE decreased non-linearly with increasing temperature. The non-linear nature of the relationship between temperature and CUE estimates (see equation 7) was much more pronounced at 1 mM cellobiose than at 20 mM cellobiose. This was a consequence of the much greater temperature response of biomass-C specific respiration and of the associated difference between biomass-C specific respiration and growth rates at 1 mM cellobiose.
FIGURE 3. Estimated microbial CUE is plotted as a function of temperature. P. fluorescens was grown in chemostats with cellobiose concentrations in the nutrient solution of either1 mM (dashed line) or 20 mM (dotted line). Shaded areas represent 95% confidence envelopes of estimated CUE.
δ13C of C Pools and Isotopic Fractionations
Temperature and C availability independently influenced δ13C of microbial biomass and respired CO2 (Table 1; Supplementary Table 2). When cellobiose was incorporated into microbial biomass and eventually respired as CO2, large fractionations were apparent at both 1 and 20 mM cellobiose (Figure 4). The value of δ13C for biomass was always more negative than δ13C of the cellobiose, and δ13C of respired CO2 was even more negative than δ13C of the biomass at both C availabilities. The degree to which temperature increased both δ13C of biomass and respired CO2 was statistically indistinguishable for both cellobiose concentrations, as illustrated by the common slope coefficient depicted in Figure 4 (slope = 0.23‰ increase in δ13C per °C, p = 0.007). The values of δ13C of biomass and respired CO2 were always more negative at 1 mM cellobiose than at 20 mM cellobiose. Assuming that P. fluorescens did not exude significant amounts of C, standard isotopic mass balance dictates that the δ13C of C taken up by the bacteria must have been more negative than the δ13C of the cellobiose substrate (Supplementary Table 1). This apparent discrimination against 13C-containing cellobiose molecules during uptake was influenced by temperature and C availability. More specifically, the lower the C availability and the lower the temperature, the greater the fractionation during uptake (the differences in δ13C of cellobiose and δ13C of microbial biomass in Figure 4). In contrast, fractionation during respiration (the differences in δ13C of microbial biomass and δ13C of respired CO2) did not vary significantly with temperature, as depicted by parallel regression lines in Figure 4.
FIGURE 4. δ13C values of provided cellobiose (the sole OC source in the system; triangles), and of microbial biomass (circles), and respired CO2 (squares) of P. fluorescens in chemostats with cellobiose concentrations in the nutrient solution of either 1 mM (A) or 20 mM (B) are plotted against bioreactor temperature. Regression lines represent statistically significant, common temperature effects on δ13C of microbial biomass and respired CO2 (see Supplementary Table 2). As indicated by arrows, uptake fractionation is the difference in δ13C values between cellobiose and biomass, and respiratory fractionation is the difference in δ13C values between biomass and respired CO2.
TABLE 1. Analysis of covariance results testing the effects of microbial C pools, temperature, and C availability on δ13C of cellobiose, biomass, and respired CO2 during cellobiose transformations.
We used chemostats and stable C isotope analyses to quantify rates and isotopic signatures of C flowing from a single C substrate through microbial biomass into CO2. Using a two C pool dynamic model for substrate and biomass coupled with measurements of pool sizes and fluxes, we estimated how microbial C uptake, C uptake affinity, and CUE change with temperature and C availability. The combination of modeling and measurements allowed us to infer changes with environmental conditions in underlying microbial physiology. Our results highlight how temperature and C availability influence biomass-C specific rates and microbial 13C discrimination during OC transformations. Below, we discuss the inferred, underlying mechanistic basis for the observed responses to varying temperature and C availability, and implications for interpreting the isotopic signature of ecosystem respiration.
Interactive Effects of Temperature and C Availability on C Uptake Affinity, Biomass-C Specific Rates, and CUE
We compared microbial C fluxes at 1 and 20 mM cellobiose (current study) with those at 10 mM cellobiose as reported in Lehmeier et al. (2016). Microbial biomass-C specific rates, C uptake affinity, and biomass-C specific C uptake rates at 10 mM cellobiose fall between those reported here at 1 and 20 mM cellobiose, suggesting gradual changes in microbial temperature responses to C availability. The shifts in C uptake affinity and C uptake rates from 1 mM to 10 mM to 20 mM are consistent with the notion that microorganisms can sense their environment and adjust their ability to acquire C as substrate availability varies (Williams et al., 1994; Ferenci, 1999; Gresham and Hong, 2015). Indeed, microorganisms possess multiple C transporter systems, each exhibiting varying affinity, and these systems’ relative expression can change with environmental C concentrations (Williams et al., 1994; Schlösser and Schrempf, 1996; Ferenci, 1999; Kajikawa and Masaki, 1999; Koning et al., 2001; Gresham and Hong, 2015). It is beyond the scope of our study to investigate the exact mechanisms driving shifts in C uptake affinity. However, enhanced C uptake affinity can occur through various mechanisms, including changes in transporter density, microbial surface to volume ratio, membrane thickness, and gene expression of a suite of transporters with different affinity (Aksnes and Egge, 1991; Button, 1993; Kovárová-Kovar and Egli, 1998; Liu and Ferenci, 2001; Gresham and Hong, 2015).
Our system is different from natural systems in that one population of microorganisms grew with surplus C at steady state without interactions with other microbial populations. This distinction allows us to make multiple observations about microbial C flows with only temperature and C availability varying, which cannot be addressed directly using environmental samples. First, our work reinforces the idea that CUE declines with rising temperature (Steinweg et al., 2008; Allison et al., 2010; Tucker et al., 2013; Wieder et al., 2013; Schindlbacher et al., 2015), though this has not always been inferred in studies using environmental samples (Crowther and Bradford, 2013; Hagerty et al., 2014). The maximum, convergent CUE around 0.7-0.8 in the current study and in Lehmeier et al. (2016) is consistent with inferred CUE values reported in aquatic systems (Hall and Cotner, 2007; Manzoni et al., 2012), soils (Dijkstra et al., 2011b; Manzoni et al., 2012; Frey et al., 2013; Blagodatskaya et al., 2014; Lee and Schmidt, 2014), and batch culture (Keiblinger et al., 2010). The ceiling value of CUE suggests that inherent limitations might keep microorganisms from achieving higher CUE values. At relatively low temperatures, microbial respiratory costs for growth, maintenance, and foraging are minimized (Rivkin and Legendre, 2001; Apple et al., 2006), likely yielding relatively high CUE. However, we expect that the negative relationship between CUE and temperature breaks down below a certain temperature limit: when microbes enter dormancy at extremely low temperature (i.e., growth rate is null), then CUE defined as growth rate/(growth rate + respiration rate) approaches zero. Thus, we highlight that current conceptions and estimates of CUE are well-defined only when C is flowing into microbial cells and that the negative CUE-temperature relationship may only hold true at temperatures at which microbes are growing.
Second, our results demonstrate that assuming a linear decline of CUE with temperature (Allison et al., 2010; Wieder et al., 2013) may not be appropriate. Indeed, CUE as defined here and in multiple other studies (see Manzoni et al., 2012) must decline with temperature in a non-linear fashion if respiration increases in a linear fashion as observed in this (Figure 1B) and a related study (Lehmeier et al., 2016). Our observation of a non-linear temperature dependence of CUE estimates emphasizes the need for further research to advance CUE estimates in Earth system models (Ballantyne and Billings in review).
Third, greater temperature responses of respiration and CUE at relatively low C availability, as observed in this study, support the idea that lower resource C:N may be generally linked to higher temperature sensitivity of OC transformations (Fierer et al., 2003; Rumpel and Kögel-Knabner, 2011; Billings and Ballantyne, 2013). We acknowledge that the meaning of resource C:N in our chemostat system (as determined by the ratio of cellobiose-C to mineral N) may be different from that in natural conditions, where microbes have varying access to both organic and inorganic mineral nutrients and resources can be protected physically and chemically. Furthermore, we do not know if changes in C availability caused nutrient limitation or a shift in relative nutrient limitation, common in natural conditions. However, if these results are robust and are not masked by the myriad other features at play in natural systems, we would predict that environments exhibiting relatively lower C concentrations such as deep, relatively old soil OC and fertilized soils/aquatic systems would be especially vulnerable to OC losses via microbial respiration with warming.
Temperature and C Availability Effects on Microbial C Isotope Discrimination
Our isotope data suggest that C availability can significantly alter microbial uptake as well as uptake discrimination against 13C-containing compounds. Bigger δ13C differences between cellobiose and microbial biomass at low C availability can be partly explained by greater C uptake affinity and associated biomass-C specific uptake rates. Shifts from a low C affinity uptake system to a high C uptake affinity system at lower C availability may lead to greater uptake fractionation, given that transport pathways with different affinities can exhibit varying degrees of C isotope fractionation (Henn and Chapela, 2000). However, even if different C uptake affinity systems have the same isotopic fractionation effect, greater uptake flux rates through a higher C uptake affinity system can enhance the magnitude of 13C discrimination during protein-mediated reactions (Tcherkez et al., 2011). Although these scenarios provide potential clues explaining why stable C uptake fractionation was higher at low C availability, they cannot explain why the degree of microbial 13C discrimination during uptake apparently was not influenced by changes in temperature (Figure 4). Thus, further studies need to examine the degree of coupling between changing C uptake affinity/uptake rate and uptake 13C fractionation under different environmental conditions.
Contrary to Lehmeier et al. (2016), we observed similar respiratory fractionation across temperatures for both cellobiose concentrations, as illustrated by two parallel regression lines in Figure 4. If we assume that this similar respiratory fractionation can be used to interpret δ13C values of microbial biomass and respired CO2 in natural samples, our results may help tease apart the effects of multiple, confounding factors driving apparent respiratory C isotope fractionation in environmental samples. In natural samples, respiratory fractionation is a net result of kinetic isotope effects, preferential substrate utilization, and heterogeneity in microbial community composition (Werth and Kuzyakov, 2010). By controlling the type and supply of substrate, and employing one population of microorganisms, our results approach intrinsic kinetic isotope effects as much as possible, with little influence of other factors. Thus, comparison of respiratory fractionation in this study to that in natural samples can be useful for quantifying the magnitude and direction of the confounding factors influencing C isotope fractionations during microbial C mineralization. Direct comparisons should be made with caution, because the range of respiratory 13C discrimination in this study (from ~5 to ~11‰, the difference between δ13C of biomass and respired CO2) is higher than apparent respiratory 13C discrimination reported in other studies of natural samples (Šantrůčková et al., 2000; Werth and Kuzyakov, 2010; Brüggemann et al., 2011). This discrepancy may arise from the relatively high C availability and decoupled C and N sources in our chemostat system. High absolute C availability in chemostats may induce microorganisms to enhance discrimination against 13C during OC transformations, given reduced discrimination against heavy stable isotopes when resources are in relatively scant supply (Fry, 2006). However, our work provides baseline data describing the potential degree of discrimination against 13C during microbially induced C flows and how it may or may not change with environmental conditions.
Ecosystem scientists frequently assume that heterotrophic microbial 13C discrimination during OC transformations is minimal and that the δ13C value of soil-respired CO2 closely reflects the δ13C of the C source (Ehleringer et al., 2000; Bowling et al., 2008). This assumption may be valid for some scenarios. However, the results reported here demonstrate that microorganisms can discriminate substantially against 13C during uptake and respiration, and that the degree of discrimination during OC uptake may change with C availability and temperature. This is relevant for interpretations of δ13C of OC substrates in natural environments, microorganisms themselves (their live biomass, and their necromass’ contribution to OC), and respired CO2. These data serve as a cautionary tale for those using δ13C signatures to infer underlying mechanisms driving those values, or those projecting ecosystem C pools and fluxes.
KM, CL, FB, and SB contributed to the development of the idea. KM and CL performed the experiment. KM, FB, and SB analyzed the data and all were involved in the data interpretation. KM wrote a draft and all contributed to revising it.
This research was supported by NSF DEB-0950095 and EAR-1331846.
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 appreciate laboratory assistance from Yanjun Chen, Mitch Sellers, Samantha Elledge, Carl Heroneme, and Dr. Kate Buckeridge in operating the chemostat experiments. We thank two reviewers for their time and input. KM also thanks Saengseon Gileum and Poohy Min for editing assistance.
The Supplementary Material for this article can be found online at: https://www.frontiersin.org/article/10.3389/fmicb.2016.02083/full#supplementary-material
Apple, J. K., del Giorgio, P. A., and Kemp, W. M. (2006). Temperature regulation of bacterial production, respiration, and growth efficiency in a temperate salt-marsh estuary. Aquat. Microb. Ecol. 43, 243–254.
Bárta, J., Šlajsová, P., Tahovská, K., Picek, T., and Šantrůčková, H. (2013). Different temperature sensitivity and kinetics of soil enzymes indicate seasonal shifts in C, N and P nutrient stoichiometry in acid forest soil. Biogeochemistry 117, 525–537. doi: 10.1007/s10533-013-9898-1
Billings, S. A., and Ballantyne, F. (2013). How interactions between microbial resource demands, soil organic matter stoichiometry, and substrate reactivity determine the direction and magnitude of soil respiratory responses to warming. Glob. Change Biol. 19, 90–102.
Billings, S. A., Tiemann, L. K., Ballantyne, F. IV, Lehmeier, C. A., and Min, K. (2015). Investigating microbial transformations of soil organic matter: synthesizing knowledge from disparate fields to guide new experimentation. Soil 1, 313–330. doi: 10.5194/soil-1-313-2015
Blagodatskaya, E., Blagodatsky, S., Anderson, T. H., and Kuzyakov, Y. (2014). Microbial growth and carbon use efficiency in the rhizosphere and root-free soil. PLoS ONE 9:e0093282. doi: 10.1371/journal.pone.0093282
Brüggemann, N., Gessler, A., Kayler, Z., Keel, S. G., Badeck, F., Barthel, M., et al. (2011). Carbon allocation and carbon isotope fluxes in the plant-soil-atmosphere continuum: a review. Biogeosciences 8, 3457–3489. doi: 10.5194/bg-8-3457-2011
Cajal-Medrano, R., and Maske, H. (2005). Growth efficiency and respiration at different growth rates in glucose-limited chemostats with natural marine bacteria populations. Aquat. Microb. Ecol. 38, 125–133. doi: 10.3354/ame038125
Chrzanowski, T. H., and Grover, J. P. (2008). Element content of Pseudomonas fluorescens varies with growth rate and temperature: a replicated chemostat study addressing ecological stoichiometry. Limnol. Oceanogr. 53, 1242–1251. doi: 10.4319/lo.2008.53.4.1242
Death, A., Notley, L., and Ferenci, T. (1993). Derepression of LamB protein facilitates outer membrane permeation of carbohydrates into Escherichia coli under conditions of nutrient stress. J. Bacteriol. 175, 1475–1483. doi: 10.1128/jb.175.10.3245.1993
Dijkstra, P., Blankinship, J. C., Selmants, P. C., Hart, S. C., Koch, G. W., Schwartz, E., et al. (2011a). Probing carbon flux patterns through soil microbial metabolic networks using parallel position-specific tracer labeling. Soil Biol. Biochem. 43, 126–132. doi: 10.1016/j.soilbio.2010.09.022
Dijkstra, P., Thomas, S. C., Heinrich, P. L., Koch, G. W., Schwartz, E., and Hungate, B. A. (2011b). Effect of temperature on metabolic activity of intact microbial communities: evidence for altered metabolic pathway activity but not for increased maintenance respiration and reduced carbon use efficiency. Soil Biol. Biochem. 43, 2023–2031. doi: 10.1016/j.soilbio.2011.05.018
Ehleringer, J. R., Buchmann, N., and Flanagan, L. B. (2000). Carbon isotope ratios in belowground carbon cycle processes. Ecol. Appl. 10, 412–422. doi: 10.1890/1051-0761 (2000)010[0412:CIRIBC]2.0.CO;2
Ferenci, T. (1999). “Growth of bacterial cultures” 50 years on: towards an uncertainty principle instead of constants in bacterial growth kinetics. Res. Microbiol. 150, 431–438. doi: 10.1016/S0923-2508(99)00114-X
Fierer, N., Allen, A. S., Schimel, J. P., and Holden, P. (2003). Controls on microbial CO2 production: a comparison of surface and subsurface soil horizons. Glob. Change Biol. 9, 1322–1332. doi: 10.1046/j.1365-2486.2003.00663.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. Change 4, 903–906. doi: 10.1038/nclimate2361
Hagström,Å., Azam, F., Andersson, A., Wikner, J., and Rassoulzadegan, F. (1988). Microbial loop in an oligotrophic pelagic marine ecosystem: possible roles of cyanobacteria and nanoflagellates in the organic fluxes. Mar. Ecol. Prog. Ser. 49, 171–178. doi: 10.3354/meps049171
Hall, E. K., and Cotner, J. B. (2007). Interactive effect of temperature and resources on carbon cycling by freshwater bacterioplankton communities. Aquat. Microb. Ecol. 49, 35–45. doi: 10.3354/ame01124
Henn, M. R., and Chapela, I. H. (2000). Differential C isotope discrimination by fungi during decomposition of C(3)- and C(4)-derived sucrose. Appl. Environ. Microbiol. 66, 4180–4186. doi: 10.1128/AEM.66.10.4180-4186.2000
Henn, M. R., Gleixner, G., and Chapela, I. H. (2002). Growth-dependent stable carbon isotope fractionation by basidiomycete fungi?: del13 C pattern and physiological process. Appl. Environ. Microbiol. 68, 4956–4964. doi: 10.1128/AEM.68.10.4956
Keiblinger, K. M., Hall, E. K., Wanek, W., Szukics, U., Hammerle, I., Ellersdorfer, G., et al. (2010). The effect of resource quantity and resource stoichiometry on microbial carbon-use-efficiency. FEMS Microbiol. Ecol. 73, 430–440. doi: 10.1111/j.1574-6941.2010.00912.x
Koning, S. M., Koning, S. M., Elferink, M. G. L., Elferink, M. G. L., Driessen, A. J. M., and Driessen, A. J. M. (2001). Cellobiose uptake in the hyperthermophilic archaeon. Microbiology 183, 4979–4984. doi: 10.1128/JB.183.17.4979
Lehmeier, C. A., Ballantyne, F. IV, Min, K., and Billings, S. A. (2016). Temperature-mediated changes in microbial carbon use efficiency and 13C discrimination. Biogeosciences 13, 3319–3329. doi: 10.5194/bgd-12-17367-2015
Lehmeier, C. A., Min, K., Niehues, N. D., Ballantyne, F. IV, and Billings, S. A. (2013). Temperature-mediated changes of exoenzyme-substrate reaction rates and their consequences for the carbon to nitrogen flow ratio of liberated resources. Soil Biol. Biochem. 57, 374–382. doi: 10.1016/j.soilbio.2012.10.030
Liu, X., and Ferenci, T. (2001). An analysis of multifactorial influences on the transcriptional control of ompF and ompC porin expression under nutrient limitation. Microbiology 147, 2981–2989. doi: 10.1099/00221287-147-11-2981
Luo, Y., Ahlstrõm, A., Allison, S. D., Batjes, N. H., Brovkin, V., Carvalhais, N., et al. (2016). Toward more realistic projections of soil carbon dynamics by Earth system models. Glob. Biogeochem. Cycle 30, 40–56. doi: 10.1002/2015GB005239
Manzoni, S., Taylor, P., Richter, A., Porporato, A., and Ågren, G. I. (2012). Environmental and stoichiometric controls on microbial carbon use efficiency in soils. New Phytol. 196, 79–91. doi: 10.1111/j.1469-8137.2012.04225.x
Min, K., Lehmeier, C. A., Ballantyne, F., Tatarko, A., and Billings, S. A. (2014). Differential effects of pH on temperature sensitivity of organic carbon and nitrogen decay. Soil Biol. Biochem. 76, 193–200. doi: 10.1016/j.soilbio.2014.05.021
Nedwell, D. (1999). Effect of low temperature on microbial growth: lowered affinity for substrates limits growth at low temperature. FEMS Microbiol. Ecol. 30, 101–111. doi: 10.1111/j.1574-6941.1999.tb00639.x
Penger, J., Conrad, R., and Blaser, M. (2014). Stable carbon isotope fractionation of six strongly fractionating microorganisms is not affected by growth temperature under laboratory conditions. Geochim. Cosmochim. Acta 140, 95–105. doi: 10.1016/j.gca.2014.05.015
Reay, D. S., Nedwell, D. B., Priddle, J., and Ellis-Evans, J. C. (1999). Temperature dependence of inorganic nitrogen uptake: reduced affinity for nitrate at suboptimal temperatures in both algae and bacteria. Appl. Environ. Microbiol. 65, 2577–2584.
Rhee, G., and Gotham, I. J. (1981). The effect of environmental factors on phytoplankton growth?: temperature and the interactions of temperature with nutrient limitation. Limnol. Oceanogr. 26, 635–648. doi: 10.1111/j.1462-2920.2011.02498.x
Šantrůčková, H., Bird, M. I., and Lloyd, J. (2000). Microbial processes and carbon-isotope fractionation in tropical and temperate grassland soils. Funct. Ecol. 14, 108–114. doi: 10.1046/j.1365-2435.2000.00402.x
Schindlbacher, A., Schnecker, J., Takriti, M., Borken, W., and Wanek, W. (2015). Microbial physioloby and soil CO2 eﬄux after 9 years of soil warming in a temperate forest – no indications for thermal adaptations. Glob. Change Biol. 21, 4265–4277. doi: 10.1111/gcb.12996
Schlösser, A., and Schrempf, H. (1996). A lipid-anchored binding protein is a component of an ATP-dependent cellobiose/cellotriose-transport system from the cellulose degrader Streptomyces reticuli. Eur. J. Biochem. 242, 332–338. doi: 10.1111/j.1432-1033.1996.0332r.x
Schulze, K. L., and Lipe, R. S. (1964). Relationship between substrate concentration, growth rate, and respiration rate of Escherichia coli in continuous culture. Arch. Mikrobiol. 48, 1–20. doi: 10.1007/BF00406595
Steinweg, J. M., Plante, A. F., Conant, R. T., Paul, E. A., and Tanaka, D. L. (2008). Patterns of substrate utilization during long-term incubations at different temperatures. Soil Biol. Biochem. 40, 2722–2728. doi: 10.1016/j.soilbio.2008.07.002
Tucker, C. L., Bell, J., Pendall, E., and Ogle, K. (2013). Does declining carbon-use efficiency explain thermal acclimation of soil respiration with warming? Glob. Change Biol. 19, 252–263. doi: 10.1111/gcb.12036
Werth, M., and Kuzyakov, Y. (2010). 13C fractionation at the root-microorganisms-soil interface: a review and outlook for partitioning studies. Soil Biol. Biochem. 42, 1372–1384. doi: 10.1016/j.soilbio.2010.04.009
Williams, S. G., Greenwood, J. A., and Jones, C. W. (1994). The effect of nutrient limitation on glycerol uptake and metabolism in continuous cultures of Pseudomonas aeruginosa. Microbiology 140, 2961–2969. doi: 10.1099/13500872-140-11-2961
Xu, X., Schimel, J. P., Thornton, P. E., Song, X., Yuan, F., and Goswami, S. (2014). Substrate and environmental controls on microbial assimilation of soil organic carbon: a framework for Earth system models. Ecol. Lett. 17, 547–555. doi: 10.1111/ele.12254
Keywords: carbon use efficiency, respiration, isotope, resource, microbial growth, chemostat, nutrient, stoichiometry
Citation: Min K, Lehmeier CA, Ballantyne F IV and Billings SA (2016) Carbon Availability Modifies Temperature Responses of Heterotrophic Microbial Respiration, Carbon Uptake Affinity, and Stable Carbon Isotope Discrimination. Front. Microbiol. 7:2083. doi: 10.3389/fmicb.2016.02083
Received: 14 October 2016; Accepted: 08 December 2016;
Published: 26 December 2016.
Edited by:Michael Sauer, University of Natural Resources and Life Sciences, Vienna, Austria
Reviewed by:Stefano Manzoni, Stockholm University, Sweden
Andreas Schindlbacher, Austrian Research Centre for Forests (BFW), Austria
Copyright © 2016 Min, Lehmeier, Ballantyne and Billings. 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: Sharon A. Billings, email@example.com
†Present address: Ford Ballantyne, Odum School of Ecology, University of Georgia, Athens, GA, USA