Metabolic tradeoffs and heterogeneity in microbial responses to temperature determine the fate of litter carbon in simulations of a warmer world

Climate change has the potential to destabilize the Earth’s massive terrestrial carbon (C) stocks, but the degree to which models project this destabilization to occur depends on the kinds and complexities of microbial processes they simulate. Of particular note is carbon use efficiency (CUE), which determines the fraction of C processed by microbes that is anabolized into microbial biomass rather than lost to the atmosphere and soil as carbon dioxide and extracellular products. The temperature sensitivity of CUE is often modeled as an intrinsically fixed (homogeneous) property of the community, which contrasts with empirical data and has unknown impacts on projected changes to the soil C cycle under global warming. We used the Decomposition Model of Enzymatic Traits (DEMENT) – which simulates taxonlevel litter decomposition dynamics – to explore the effects of introducing organism-level heterogeneity into the CUE response to temperature for decomposition of leaf litter under 5 C of warming. We found that allowing the CUE temperature response to differ between taxa facilitated increased loss of litter C, unless fungal taxa were specifically restricted to decreasing CUE with temperature. Litter C loss was exacerbated by variable and elevated CUE at higher temperature, which effectively lowered costs for extracellular enzyme production. Together these results implicate a role for diversity of taxon-level CUE responses in driving the fate of litter C in a warmer world within DEMENT, which should be explored within the framework of additional model structures and validated with empirical studies.


15
Soil heterotrophs are central to the cycling and recycling of the 60 gigatons of organic carbon (C) that plants deposit onto and into the ground each year. How well these inputs are converted into relatively stable soil organic matter depends on temperature, moisture, chemical composition, and soil mineralogy, which interact to influence microbial physiology (Manzoni et al., 2012;Kallenbach et al., 2016;Oldfield et al., 2018). Predictions regarding how soil C stocks will respond to climate change are, in turn, highly sensitive to how carbon use efficiency (CUE)-or the fraction of C taken up by a cell and incorporated into 20 biomass rather than being respired--changes with temperature Wieder et al., 2013;Allison, 2014;Li et al., 2014;Sistla et al., 2014;Tang and Riley, 2015). As such, quantifying microbial decomposer CUE and its responsiveness to environmental change has been subject to intensive study (Devêvre and Horwáth, 2000;Frey et al., 2013;Blagodatskaya et al., 2014;Lee and Schmidt, 2014;Spohn et al., 2016a, b;Öquist et al., 2017;Malik et al., 2018;Geyer et al., 2019;Malik et al., 2019;Zheng et al., 2019). 25 Soil microbial communities show considerable differences in how their metabolisms respond to elevated temperatures, with their CUE increasing (Öquist et al., 2017;Zheng et al., 2019), decreasing (Devêvre and Horwáth, 2000;Frey et al., 2013;Öquist et al., 2017;Li et al., 2018;Zheng et al., 2019) or remaining unaffected by warming (Dijkstra et al., 2011b;Öquist et al., 2017;Walker et al., 2018;Zheng et al., 2019). However, models of the soil C cycle generally assume either no change Li et al., 2014;Wieder et al., 2014) or a fixed and homogeneous decrease in CUE with temperature 30 Wieder et al., 2013;Allison, 2014;Li et al., 2014). When CUE is allowed to directly increase with temperature, this temperature response is homogeneous across taxa (Frey et al., 2013;Ye et al., 2019). In other instances, CUE may be modeled as fixed within taxa, such that dynamic changes in community-level CUE with warming are the result of shifts in the dominant group or groups of organisms present as a function of their dietary preferences and/or C: nitrogen (N) ratio (Wieder et al., 2013;Sistla et al., 2014), rather than any inherent differences in the temperature sensitivity of the constituent 35 community members. This may occur, for instance, if large C-rich fungi show less-positive responses to warming than small, N-rich bacteria do (Pietikäinen et al., 2005;DeAngelis et al., 2015). Therefore, models have thus far insufficiently accounted for how the temperature sensitivity of central metabolism may differ between microbes, such that intrinsic differences in efficiency between taxa ("heterogeneity") above and beyond temperature-driven differences in substrate supply may also drive microbial community trajectories. 40 Heterogeneity in the temperature sensitivity of growth efficiency across taxa could be driven by differences in the ratelimiting step of central metabolic pathways (Dijkstra et al., 2011a), or in how well the proteins responsible for the extracellular processing and uptake of environmental nutrients are able to maintain activity as temperature increases Alster et al., 2018). For instance, there is some evidence that bacteria benefit more than fungi from an increase in temperature, as their growth rate has been observed to decrease less rapidly with temperature above its optimum than a fungal community 45 (Pietikäinen et al., 2005). The respiration rate for the two groups could not be isolated in that study, and so direct differential effects of temperature on CUE could not be parsed out. However, the CUE has been observed to differ in how sensitive it is to nutrient limitation for bacteria and fungi, such that differences in the temperature sensitivity of CUE between the two groups may play out as a consequence of changing nutrient demands (Keiblinger et al., 2010;Sinsabaugh et al., 2016). Therefore, the temperature range over which an an organism can maintain efficient growth is one important dimension of its 50 niche (Cavicchioli, 2016) and may differ between taxa.
The temperature response of CUE may also differ between taxa because of varying investments in extracellular enzyme production. Extracellular enzyme production can impose substantial metabolic costs on the cell (Allison, 2014;Malik et al., 2019), as C which could otherwise be allocated to growing the cell must instead be spent producing amino acids and ATP to synthesize the enzymes (Kaleta et al., 2013;Kafri et al., 2016). Thus if taxa need to produce more extracellular enzymes in 55 order to support rapid growth at higher temperature, then their CUE could decrease with temperature. However, some taxa may have mechanisms to adjust central metabolism so growth remains efficient at elevated temperatures, as hinted by a weak positive correlation between the Q10 of CUE and cellulolytic potential in soil communities (Zheng et al., 2019). Nonetheless, the most probable relationship between CUE temperature response and extracellular investment is currently unknown.
We explored whether interactions between the temperature sensitivity of intracellular (i.e. CUE) and extracellular (i.e. litter 60 decomposing enzyme) metabolic processes of cells can explain why CUE is observed to increase with temperature in some soils, and decrease in others. We used the litter decomposition model DEMENT (Allison, 2012) to evaluate four hypotheses: 1) allowing temperature response of CUE to differ (be heterogeneous) between taxa increases uncertainty in projected litter decomposition dynamics because more diverse phenotypic combinations exist for competitive selection (i.e. species sorting) to act upon; 2) this heterogeneity favors a community with higher CUE, in turn leading to higher microbial biomass and greater 65 litter C loss with warming; 3) forcing the temperature response of CUE to increase with the number of enzymes an organism produces causes greater litter C loss than when the two factors vary independently, because increasing CUE with temperature offsets the increased costs against CUE associated with copious enzyme production; and 4) the magnitude of litter C loss with warming is greater when the CUE of the C-rich fungal functional group increases with temperature than if only the N-rich bacterial functional group does due to higher relative C demand of the former.

DEMENT background and model design
DEMENT (Allison, 2012) is a litter decomposition model designed to simulate the loss of leaf C through time. The principal advancement of DEMENT over its predecessors is that it is both microbially-and spatially-explicit. The model is able to simulate inter-and intraspecific microbial interactions, with a primary focus on the tradeoff between the ability to take up 75 and digest substrates, and the metabolic costs of creating and maintaining the machinery required to do so. Because these tradeoffs are both explicit and variable across taxa, DEMENT is an ideal model for evaluating how the physiology and ecology of microbes affects C stocks in a changing world. Furthermore, DEMENT allows for consideration of how heterogeneous responses across taxa (rather than using some homogeneous cross-taxon mean) can facilitate soil C responses to climate change. Full details about the setup and execution of DEMENT are available elsewhere (Allison, 2012(Allison, , 2014Allison and 80 Goulden, 2017); here we describe the controls on CUE in the model which are relevant to our study.
Intrinsic CUE-the maximum CUE an organism could attain under ideal temperature and stoichiometry-is calculated for each taxon as a function of the baseline CUE at 15 o C (C r ), and the number of enzymes (N e ) and uptake transporters (N u ) the taxon can produce. In turn, how much CUE is decreased due to enzyme and transporter production depends on the cost per enzyme (C e ) and cost per transporter (C u ) ( Table 1). The C used in enzyme synthesis is considered a loss from the cell, 85 and is therefore not reported as microbial biomass C. The intrinsic CUE of each taxon is adjusted for temperature, decreasing by 0.016 o C -1 by default (i.e. C t = -0.016 o C -1 ), consistent with a global meta-analysis (Qiao et al., 2019). Therefore, CUE is calculated as CU E = C i + (temperature − 288) * C t + C e * N e + C u * N u , where temperature is in Kelvin.

Modifications to DEMENT
DEMENT v0.7.2 was downloaded from GitHib (https://github.com/stevenallison/DEMENT), and modified as follows. C r was 90 adjusted downwards from its original published value of 0.58 to 0.38 at 15 o C; this not only improved model stability (Table   S1), but is also consistent with a comparative modeling study completed by Li et al. (2014), several 18 O-H 2 O based CUE measurements (Spohn et al., 2016a, b;Geyer et al., 2019), and for the structural components of litter modeled by the soil C model MIMICS (Wieder et al., 2015a). By default, the temperature sensitivity of CUE (C t ) is fixed to take on the same value (ie is homogeneous) for all taxa; therefore we modified the model so C t could vary around the mean in different ways ( Figure 1).

95
In the first set of scenarios, C t varied independent of the taxonomic identity or number of enzymes a taxon produced ( Figure 1B and hypotheses 1 and 2). In the second scenario, C t was limited to either increasing ( Figure 1C) or decreasing ( Figure 1D) as a function of the number of enzymes a given taxon had (hypothesis 3). In the third, bacteria were constrained to have a positive and fungi a negative C t ( Figure 1E), fungi a positive and bacteria a negative C t ( Figure 1F), both bacteria and fungi to have a positive C t ( Figure 1G), or both have a negative C t ( Figure 1H; hypothesis 4). In all instances, C t was selected at random from 100 a uniform distribution bounded by +/-0.022 o C -1 at the upper and/or lower limits (scenarios B-H, Fig. 1), or assigned a fixed value equivalent to the cross-taxon mean (scenario A, Fig. 1; scenario A i , A ii Fig. S1). These values are within the range of temperature sensitivities observed for both bacterial cultures in the lab and for field communities (Figure 2), as well as values inferred based on modeling CUE against mean annual temperature on a global basis (Sinsabaugh et al., 2017;Ye et al., 2019).
It was necessary to force the temperature sensitivity of CUE to take on a zero-centered uniform distribution so that simulation 105 outputs in which extracellular enzyme counts were linked to the C t could be compared to those scenarios where they were not linked, without changing the distribution of extracellular enzyme counts present in the community.

Running DEMENT
DEMENT was run on the Massachusetts Green High Performance Computing Cluster for 6,000 model days using 59 different independent starting seeds and a 100x100 grid size. "Control" runs were completed at 15 o C (equivalent to April to November 110 mean soil temperature for a northern mid-latitude temperate deciduous forest (Boose, 2001)), while "heated" runs were completed at 20 o C (Allison, 2014). The first 1000 days of each resultant output file was excluded from the analysis because of rapid shifts in the microbial community during this time. In addition, outputs were filtered to exclude any seeds where the substrate pool was two or more times greater at the end of the model run than the median during the preceding 5000 days, indicating unrealistic, unconstrained litter accumulation. R version 3.4.0 was used for all runs and analyses (R Core Team, 2016). A full set of parameters and the model used to run all these simulations can be found in the supplementary files "params.txt" and "DEMENTmodel.R", respectively.

Analysis of outputs
The model outputs of interest were litter organic matter (LOM), microbial biomass carbon (MBC), respiration rate, richness and diversity of the surviving community, median number of enzymes per taxon for taxa alive during the 5000-day 120 simulation, fungal:bacterial biomass ratio for surviving taxa, and reference and simulation temperature CUE. Two diversity metrics -richness and Shannon's H -were calculated using median daily values for the microbial community over the 5000 days of simulation and the vegan package (Oksanen et al., 2017). Biomass-weighted CUE was reported for communities at 15 o C ("reference temperature CUE") or at the simulation temperature ("simulation temperature CUE"), calculated using the formula 6000 d=1000 1000 and 6000 are the days of the simulation the outputs were examined over; 100 is the number of taxa the model was initiated with; Ne i is the number of enzyme genes taxon i has; and Nu i is the number of uptake transporter genes the taxon has; and temperature is the simulation temperature in o K.
In order to determine whether warming and model parameterization affected model outputs, we used mixed effect models with starting seed as a random effect and warming or simulation scenario as fixed effects using lmer in lme4 v 1.1-17 (Bates 130 et al., 2015). Data were visually assessed for normality and homoskedasticity using qqplots and residual plots following logtransformation. Significantly different pairwise differences were subsequently identified using emmeans v.1.3.0 (Lenth et al., 2019), with a stringent Bonferoni-corrected p-value cutoff of P < 0.0001. Warming effect sizes are plotted as the natural log ratio of model outputs in heated:control scenarios. Figures were generated using ggplot2 (Wickham, 2009), and asterisks and letters denoting significant effects of experimental factors were added in Inkscape (Inkscape'sContributors, 2003-2019).

135
3 Results and discussion LOM and MBC content were both generally higher than observed in environmental samples, leading to MBC:LOM ratios at the high end of ranges observed in the field (2-11% vs. 1-5% (Santos et al., 2012;Xu et al., 2013)). LOM and MBC values were within the range previously observed for simulations using DEMENT with daily litter inputs , but greater than those with just a single litter pulse (Allison, 2012;Allison and Goulden, 2017;Evans et al., 2017), indicating that 140 these high biomass and litter C values can be attributed to these substrate inputs.

Intertaxon variability
To evaluate the effect of intertaxon CUE variablity on LOM stocks, we ran the model at 20 o C ("heated") under two scenarios, and then compared the results to runs at 15 o C ("control"). In the first "homogeneous" scenario all taxa had an identical C t , equal to 0 o C -1 ( Figure 1A). In the second "heterogeneous" scenario, C t was assigned from a random uniform distribution bounded 145 by -0.022 and 0.022 o C -1 ( Figure 1B). Therefore, the mean C t of the heterogeneous scenario was 0 o C -1 , matching that of the homogeneous scenario.
Introducing intertaxon differences in CUE temperature response caused the characteristics of the initial microbial community (starting seed) to have a greater impact on litter decomposition than when all taxa had an identical temperature response (Table 2). This contrasts with the dampening effect proposed to explain instability in small-scale microbially-explicit models 150 compared to their macroscale counterparts (Wieder et al., 2015b). Specifically, Wieder et. al. proposed integrating diverse physiology into C-cycling models should allow different microbial subpopulations to follow distinct trajectories which average out to a more consistent community-level mean and greater certainty in model projections. However, we found that the additive effect of increased physiological diversity was to increase, rather than decrease, uncertainty in the present simulations. The median-standardized interquartile ranges of both MBC (0.25 vs. 0.16) and LOM (0.28 vs. 0.14) increased with the introduction 155 of a variable C t . Through species sorting, this heterogeneously-responding microbial community became less diverse than both the homogeneous and control communities ( Table 2). The communities characterized by heterogeneous C t maintained a higher median microbial biomass-driving two and a half times more LOM loss-than the communities characterized by a homogeneous C t ( Figure 3). Intriguingly, neither litter (r=0.16, P = 0.23) nor microbial biomass pool sizes (r=-0.44, P < 0.001) positively correlated with extracellular enzyme investment; thus, a (non-significant) 28% increase in the median enzyme count 160 is unlikely to have driven the increased decomposition under the heterogeneous scenario. Instead, increased decomposition and increased biomass are likely the consequence of elevated CUE under warming conditions. Nonetheless, the non-linearity of the model limits the degree to which causal relationships can be drawn between changes in litter C and the microbial parameters.
The homogeneous community scenario tested here is akin to the fixed "no adaptation of CUE" scenario reported in a number of other studies Li et al., 2014;Sistla et al., 2014), because the cross-taxon mean used is zero temperature 165 response. However, the effect of inter-taxon differences in C t ("heterogeneous") is less studied. Our results of reduced LOM loss in the absence of acclimation are consistent with two previous studies, but contrast with others. In an ecosystem-level model parameterized for an arctic tundra system, Sistla et al. (2014) found that greater soil organic matter (SOM) loss occurred with warming when the microbial community was able to dynamically acclimate its CN ratio (and in turn efficiency), than when the CUE was effectively fixed. Likewise, Allison (2014) found greater potential for increased LOM accumulation under warming 170 when there was greater absolute variation in CUE across taxa (i.e. C r from 0.18 to 0.58 rather than 0.38 to 0.58) (Allison, 2014), although C t was always consistent across taxa in these simulations. On the other hand, a comparison of models where taxon-level differences in CUE or C t are not possible (i.e. C t is intrinsically homogeneous at the community level) showed that soil organic matter loss increases when organisms do not adapt (Li et al., 2014). Similarly, Wieder et al. found that greater SOM loss occurred if the CUE was directly insensitive to temperature than when C t was negative (Wieder et al., 2014). These 175 microbially-explicit decomposition models vary in if and how they link CUE to microbial traits, and so our findings support the concept that nuances in how different components of CUE respond to warming is an important control on the fate of litter C (Hagerty et al., 2018).

Confirming the role of C t as an additional niche dimension
We allowed for CUE to increase with temperature for a subset of taxa in a way that most previous modeling efforts have 180 not, and so it is possible that our results deviate from those of prior studies not because of variation in C t , but rather because our simulations explore novel (positive C t ) parameter space. To facilitate comparison with previous decomposition modeling studies, we ran DEMENT simulations to test the effect of C t being homogeneous vs. heterogeneous when CUE was either always positive (homogeneous C t = 0.011 o C -1 (Figure S1A i ); heterogeneous C t = 0 to 0.022 o C -1 (Figure S1G)) or always negative (homogeneous C t = -0.011 o C -1 (Figure S1A ii ); heterogeneous C t = -0.022 to 0 o C -1 ( Figure S1H)). In contrast to 185 when C t was allowed to vary over the whole spectrum of values, introducing heterogeneity in CUE did not increase interrun uncertainty in LOM or MBC pools (Table S2). We also found that less LOM accumulated when CUE showed a variable decrease with warming than a fixed one ( Figure S2), which could be attributed to a reduction in MBC. By contrast, the homogeneous zero-centered and homogeneous positive C t scenarios, and the heterogeneous zero-centered and heterogeneous positive C t , behaved more similarly to one-another in that warming decreased LOM while increasing MBC and CUE to a 190 greater degree in the heterogeneous than homogeneous scenarios (Table S2). This finding reinforces the idea that if warming favors decomposer taxa capable of maintaining efficient growth, then soil C loss will be accelerated. Nonetheless, the strongly selected-for positive CUE response is rarely observed in complex soil communities. This indicates that additional tradeoffs with CUE temperature response are likely at play when CUE is either unaffected or decreases with temperature, but that these tradeoffs are missing in the formulation of DEMENT used in this scenario. One such tradeoff possible to explore within the 195 framework of DEMENT is the allocation of resources to extracellular enzyme activity.

Linkages between CUE temperature response and extracellular enzyme allocation
Microbes depend upon extracellular enzymes to break down substrates in the environment into digestible pieces, and enzyme activities can be, like CUE, responsive to temperature German et al., 2012;Allison et al., 2018). Soil extracellular enzymes often are active in-situ at temperatures much below their apparent activity optima (German et al., 2012;200 Pold et al., 2017;Alster et al., 2018). Therefore, warming is assumed to enable them to process substrates at a higher rate in DEMENT (Allison, 2012), increasing the supply of growth substrates to microbes. However, the affinity of enzymes for their substrates may also decrease as temperature increases (German et al., 2012;Allison et al., 2018), as is assumed in DEMENT. If this is the case, then unless enzyme V max increases faster with temperature than K m , additional resources must be diverted from growth to enzyme production to maintain microbial growth substrate supply rate. Therefore, taxa may differentially-allocate 205 resources to enzymes and so demonstrate a relationship between the temperature sensitivity of CUE and the number of enzymes they produce.
We evaluated whether litter decomposition changed its trajectory when the organisms with the greatest genomic potential to break the litter down (i.e. enzyme counts) also showed the most-or least-positive growth efficiency response to warming. In the "increase" scenario, we simulated a positive relationship between temperature sensitivity of CUE and extracellular enzymes, 210 where C t increased linearly from -0.022 o C -1 for organisms with no extracellular enzyme production potential to 0.022 o C -1 for those organisms capable of producing the model maximum of 40 enzymes ( Figure 1C). In the "decrease" scenario, we simulated a negative relationship between temperature sensitivity of CUE and extracellular enzymes, where the opposite relationship was imposed with C t decreasing with enzyme counts ( Figure 1D). These scenarios were then compared to the "heterogeneous" scenario (aka "no relation", as described above), where C t varied across the same range, but independently of the number of 215 enzymes an organism could produce. Therefore, the starting distribution of C t and enzymes per taxon was identical across scenarios, and only their relationship with one-another changed.
More taxa survived to the end of the simulation when warming was applied under the "increase" scenario than either the "decrease" or "no relation" scenario (median of 13 versus 7 and 8, respectively, median absolute deviation = 2.97 in all cases).
Under the "increase" scenario, taxa had 70% more enzymes each than the "no relation" scenario, and more than three times 220 as much as the "decrease" scenario (Table 2). This relationship caused the CUE of surviving taxa to be 20-37% lower at 15 o C for the "increase" scenario compared to the others, but this deficit was diminished at 20 o C. As a result, the "increase" scenario led to higher respiration and a greater LOM loss under warming than under the "decrease" scenario, despite an overall smaller microbial biomass pool (Figure 4).
How the relationship between C t and extracellular enzymes drives favorable trait combinations in DEMENT can also be 225 observed in Figure 5. Surviving taxa retained a median enzyme count of at least 30 and a realized CUE temperature response of no less than 0.0158 o C -1 under the "increasing" scenario (ρ=0.64, P < 0.001), but there was no relationship between realized CUE temperature response and enzyme production under either the "decrease" or "no relation" scenarios. The selection for community capable of maintaining high CUE at high temperatures was much weaker when it was associated with reduced enzyme production. When there was no relationship between C t and extracellular enzyme production costs, however, 230 communities were able to attain a high realized CUE temperature response over a much wider range of median enzyme costs.
These findings indicate that response traits-which determine how an organism reacts to changes in temperature (e.g. CUE temperature response)-and effect traits-which determine how an organism alters its environment (e.g. litter decomposition potential)-interact to determine the fate of organic C within DEMENT. However, contrary to our hypothesis, adjusting DEMENT to allow for this tradeoff did not substantially alter how community-level CUE responds to temperature. The obser-235 vation that LOM is reduced further when enzyme production is effectively cheaper contrasts with earlier work with DEMENT (Allison, 2014) showing smaller litter C pools under both ambient and elevated temperature when enzymes and transporters were cheaper to produce. However, our results are consistent in that microbial biomass was lower when enzyme costs are high, and that the microbial community was able to maintain a higher CUE under warming no matter the enzyme costs. The mechanisms underlying these phenomenological similarities differ, however, due to differences in how C t was parameterized 240 in the two sets of model simulations. Specifically, although microbes were able to attain high CUE at elevated temperatures in our simulations by balancing the benefits of elevated CUE at higher temperatures with the costs of enzyme production against CUE, CUE always decreased with temperature in earlier work with DEMENT (Allison, 2014). Furthermore, enzyme production costs varied both with and independently of enzyme counts in previous DEMENT simulations (Allison, 2014).
Within the framework of DEMENT, increased CUE is likely needed to offset the costs of extracellular enzyme production 245 that allow taxa to remain competitive at elevated temperatures, However, there is a paucity of empirical evidence regarding the hypothesized correlations between temperature sensitivity of CUE and enzyme investment in soil systems. By examining correlations between the number of enzymes an organism can produce and its CUE temperature response at the end of the DEMENT model run, we see that there is likely to be either no or a positive correlation between the two variables, rather than a negative one ( Figure 6). Limited data from bacterial isolates grown in the lab also support this, whereby CUE temperature 250 response is either positively correlated or uncorrelated with the number of enzymes produced (Pold et al., in prep). Furthermore, although the mechanisms underlying the isolate response remain unclear, it is consistent with the scenario of C t and enzyme counts being positively correlated. Specifically, we found isolates with lower CUE at 15 o C (more extracellular enzymes in DEMENT) were more likely to have a positive CUE temperature response than those with a higher CUE (fewer enzymes in DEMENT). Together, these insights support a synergism between CUE temperature response and enzyme production, rather 255 than a tradeoff within the model. Because our DEMENT simulations indicate that selection for organisms characterized by high, positive CUE temperature responses with warming can alter both the directionality and extent of projected C loss, we propose it is important for other models to explore how possible increases-rather than just decreases-in C t affect terrestrial C projections.

260
Across scenarios, fungi generally dominated the microbial biomass C pool (Table 2), as is typical for litter decomposition (Chapman et al. (2013), and references therein). This pattern occurred despite generally lower biomass-weighted CUE for surviving fungal taxa, and preferential loss of fungal taxa across most scenarios ( Table 2). The lower CUE for surviving fungi was not driven by higher enzyme costs than for bacteria, as median biomass-weighed enzyme costs were not statistically different (pairwise t-test P > 0.4) and differed by less than one enzyme for the two groups. To test whether modeled differences 265 in fungal vs. bacterial cell sizes and stoichiometry were driving this pattern, we tested how forcing differences in C t in the two groups would impact the decomposition rate.
DEMENT was run with CUE simulated to respond to temperature: 1) negatively for all fungi and positively for all bacteria (F-B+; Figure 1E); 2) negatively for all bacteria and positively for all fungi (F+B-; Figure 1F); 3) positively for all bacterial and fungal taxa (F+B+; Figure 1G) or 4) negatively for all bacterial and fungal taxa (F-B-; Figure 1F). Minimum and maximum 270 C t were set to -0.022 o C and 0.022 o C, respectively. All four scenarios were tested in order to isolate the effect of changing taxonomic domain-C t relationships from simply changing C t or taxonomic domain independently.
Litter C accumulated at higher rates with warming when fungal C t was negative, regardless of the bacterial C t (Figure 7, blue, orange). This is consistent with the observation that the CUE of surviving fungi was lower at the simulation temperature in seven out of the eleven scenarios. Because fungi have higher CNP ratios than bacteria (and thus higher C demands per 275 unit biomass), we predicted that if fungi have a negative CUE temperature response, they would be weaker competitors at higher temperature than bacterial taxa, reducing their C demand and mitigating the warming effect on SOC stocks. While this contrasts with the premise that fungi should have a higher CUE  due to their higher CN ratio (Zak et al., 1996), it is consistent with a growing body of literature indicating that substrate quality-rather than the F:B ratio correlated with it-is the underlying driver of differences in CUE between soils (Frey et al., 2013;Thiet et al., 2006;Malik et al., 2018;280 Soares and Rousk, 2019). Despite the low nutrient content of the daily inputs to the model (92:0.26:0.02 C:N:P), microbes did not show evidence for nutrient limitation as biomass CN and CP ratios were lower than are typical for soil communities (Xu et al., 2013) (4.1 and 36.7 vs. 7.6 and 42.4, respectively). Thus C was limiting, which could have further disfavored the highly C-demanding fungi when their C t was negative. The lower CUE (and increased sensitivity to warming) for fungi compared to bacteria under a given scenario was also not driven by increased metabolic costs for enzyme production in fungi, as median 285 biomass-weighted enzyme counts were statistically indistinguishable from those in bacteria.
The litter C pool decreased when fungal CUE increased with temperature (Figure 7), correlating with a smaller microbial biomass pool when bacterial CUE also decreased with temperature. By contrast, the LOM and MBC responses to warming were similar when only the C t of bacteria was changed, indicating that it is the fungal warming response which really drives changes in litter decomposition in DEMENT. This result is interesting because no a priori differences in decomposition or 290 uptake potential were imposed on the two groups, and fungal and bacterial richness was initially equivalent. Nonetheless, differences in C and nutrient translocation abilities, and in cell size, stoichiometry and turnover rates still defined the two groups. Warming decreased the enzyme costs when fungal C t was positive but bacterial C t was negative, and decreased them under the opposing scenario, as evidenced by an increase in C r in the former and decrease in the latter. Nonetheless, as long as both bacterial and fungal CUE did not both decrease their CUE with temperature, community level CUE remained higher at 295 20 o C than it was at 15 o C. Empirical evidence for high-level differences in the temperature sensitivity of CUE in bacteria and fungi is currently mixed, but indicate CUE temperature response for fungi is unlikely to be more positive than that for bacteria. Zheng et al. (2019) did not find a correlation between the lipid-based fungal:bacterial ratio and Q 10 of CUE over a range of soils. However, we (Pold et al., in prep) and our colleagues (Eric Morrison, personal comment) have found that fungi tend to show a stronger 300 negative CUE response with warming than do bacteria when examining them in isolation in the lab. This is consistent with the observation that fungal CUE decreases more strongly with warming than bacterial CUE does when C t is restricted to negative values ( Figure 7). It is also consistent with the premise that bacterial growth benefits more from elevated temperature than does fungal growth in some soils (Pietikäinen et al., 2005). Greater empirical insight into the taxonomic drivers of the temperature sensitivity of CUE will assist with constraining the parameterization and projections of microbially-explicit decomposition 305 models such as DEMENT.

Comparison to empirical warming studies
Litter decomposition is typically observed to accelerate under warming (Lu et al., 2013). However, both the chemical composition of the litter and the identity of the living plant community at the site of decomposition are important for the magnitude of this response (Cornelissen et al., 2007;Ward et al., 2015). Consistent with these empirical studies-but inconsistent with a pre-310 vious publication using DEMENT (Allison, 2014)-we found that litter decomposition was accelerated by warming in seven of ten scenarios. The range of losses and gains of litter C we observed with warming (-62% (scenario C) to +42% (scenario A i )) approximates the -65% to +36% change observed in field experiments (Lu et al., 2013), with the upper limit only being exceeded when C t is constrained to negative values. Likewise, values for simulated litter respiration response (-5 to +6%) were within a narrow range compared to those observed for soil respiration in the field (-48 to +178%), and responses for microbial 315 biomass C were also within the observed range (-25 to +58% vs -47 to +86%) (Lu et al., 2013). Our modeled responses to warming thus suggest that one possible explanation for differences in terrestrial C pool responses to warming may be diverse temperature sensitivities of underlying decomposer communities. Nonetheless, a number of additional factors must be taken into consideration when interpreting our results in the context of global climate change, including soil mineral-mediated modulation of substrate supply (Schimel and Schaeffer, 2012;Coward et al., 2018), plant-microbe feedbacks (Melillo et al., 2011;320 Sistla et al., 2014;Suseela and Tharayil, 2018), and temporal variation in temperature. Furthermore, we modeled CUE as a somewhat inflexible property of individual taxa, rather than as the emergent property of allocation to different, biochemicallyimportant physiological processes that it is Hagerty et al. (2018). These are important features to pursue in future iterations of the model should address.

325
Our results indicate that accounting for heterogeneous temperature response increases uncertainty regarding future litter C stocks, but only when C t does not differ from zero on average. However, by combining simulations, empirical studies, and literature searches, we can conclude that microbes with high enzyme costs are likely to have larger increases in intrinsic CUE with temperature; that taxa can sort on a CUE temperature response axis; and that fungi are more likely to increase CUE with warming than bacteria. The simulations meeting each or all of these criteria lead to loss of litter C under warming, and so 330 DEMENT favors the prediction that litter will to become a net atmospheric C source in a warmer world. Nonetheless, our results also indicate that soil C stocks may change heterogeneously with warming because of diversity in how decomposers' CUE responds to temperature. We encourage models functioning on larger scales to explore the effect of including heterogeneity in the temperature response of CUE in order to determine the robustness of our conclusions to other model structures. However, ultimately increased integration of the growing body of literature on the temperature sensitivity of CUE must be explored for 335 root causes of heterogeneity in temperature sensitivity of CUE in taxa under in situ conditions. Dijkstra, P., Thomas, S. C., Heinrich, P. L., Koch, G. W., Schwartz, E., and Hungate, B. A.: Effect of temperature on metabolic activity of intact microbial communities: evidence for altered metabolic pathway activity but not for increased maintenance respiration and reduced   to control conditions. CUE temperature response was forced to increase, decrease, or remain independent of the number of enzymes a taxon could produce. Values above the zero line indicate warming increased the value, and values below indicate a decrease with warming.
Boxplots denote 1st to 3rd quartiles with the median. Asterisks denote significant warming effect at P < 0.0001 after correcting for multiple comparisons using the Bonferoni method. Letters denote warmed scenarios which are significantly different from one-another by the same criteria. "Reference temperature CUE" denotes the ratio of CUE for the surviving community at 15 o C in warmed to the ratio of CUE for the surviving community at 15 o control scenarios, while "Simulation temperature CUE" denotes the ratio of CUE in heated scenarios run at 20 o C to control scenarios run at 15 o C. , and values below indicate a decrease with warming. Boxplots denote 1st to 3rd quartiles with the median. Asterisks denote significant warming effect at P < 0.0001 after correcting for multiple comparisons using the Bonferoni method. Letters denote warmed scenarios which are significantly different from one-another by the same criteria. "Reference temperature CUE" denotes the ratio of CUE for the surviving community at 15 o C in warmed to the ratio of CUE for the surviving community at 15 o control scenarios, while "Simulation temperature CUE" denotes the ratio of CUE in heated scenarios run at 20 o C to control scenarios run at 15 o C. Values for "simulation temperature CUE" are also shown for bacteria and fungi separately.  Table 2. Median (or median-standardized interquartile range (IQRm)) output values for various DEMENT model runs, marked according to warming effect (+/-) and model structure effects (letters) determined using Bonferoni-corrected post-hoc tests following linear mixed effect models. Symbols: "+" warming increased value; "-" warming decreased value). Letters: differences between warmed scenarios. All heated scenarios were compared to values in the "control" column to determine the warming effect. Only values within boxes defined by vertical lines were evaluated for significance differences of warming scenario, as different seeds needed to be excluded for failing to constrain litter accumulation in the two boxes. NA indicates that fungi died out completely in many instances (i.e. median fungal biomass of zero), so the parameter output could not be determined. The number in brackets denotes the median excluding the scenarios where all the fungi died out. "Excluded runs" denotes the number of runs removed from analysis due to unconstrained litter C accumulation