Metatranscriptomes of California grassland soil microbial communities in response to rewetting.

When very dry soil is rewet, rapid stimulation of microbial activity has important implications for ecosystem biogeochemistry, yet associated changes in microbial transcription are poorly known. Here, we present metatranscriptomes of California annual grassland soil microbial communities, collected over 1 week from soils rewet after a summer drought-providing a time series of short-term transcriptional response during rewetting.

Arbuscular mycorrhiza convey significant plant carbon to a diverse hyphosphere microbial food web and mineral-associated organic matter.

Arbuscular mycorrhizal fungi (AMF) transport substantial plant carbon (C) that serves as a substrate for soil organisms, a precursor of soil organic matter (SOM), and a driver of soil microbial dynamics. Using two-chamber microcosms where an air gap isolated AMF from roots, we 13CO2-labeled Avena barbata for 6 wk and measured the C Rhizophagus intraradices transferred to SOM and hyphosphere microorganisms. NanoSIMS imaging revealed hyphae and roots had similar 13C enrichment. SOM density fractionation, 13C NMR, and IRMS showed AMF transferred 0.77 mg C g-1 of soil (increasing total C by 2% relative to non-mycorrhizal controls); 33% was found in occluded or mineral-associated pools. In the AMF hyphosphere, there was no overall change in community diversity but 36 bacterial ASVs significantly changed in relative abundance. With stable isotope probing (SIP)-enabled shotgun sequencing, we found taxa from the Solibacterales, Sphingobacteriales, Myxococcales, and Nitrososphaerales (ammonium oxidizing archaea) were highly enriched in AMF-imported 13C (> 20 atom%). Mapping sequences from 13C-SIP metagenomes to total ASVs showed at least 92 bacteria and archaea were significantly 13C-enriched. Our results illustrate the quantitative and ecological impact of hyphal C transport on the formation of potentially protective SOM pools and microbial roles in the AMF hyphosphere soil food web.

Soil microbial community response to corrinoids is shaped by a natural reservoir of vitamin B 12.

Soil microbial communities perform critical ecosystem services through the collective metabolic activities of numerous individual organisms. Most microbes use corrinoids, a structurally diverse family of cofactors related to vitamin B 12 . Corrinoid structure influences the growth of individual microbes, yet how these growth responses scale to the community level remains unknown. Analysis of metagenome-assembled genomes suggests corrinoids are supplied to the community by members of the archaeal and bacterial phyla Thermoproteota , Actinobacteria , and Proteobacteria . Corrinoids were found largely adhered to the soil matrix in a grassland soil, at levels exceeding those required by cultured bacteria. Enrichment cultures and soil microcosms seeded with different corrinoids showed distinct shifts in bacterial community composition, supporting the hypothesis that corrinoid structure can shape communities. Environmental context influenced both community and taxon-specific responses to specific corrinoids. These results implicate corrinoids as key determinants of soil microbiome structure and suggest that environmental micronutrient reservoirs promote community stability.

Predictions of rhizosphere microbiome dynamics with a genome-informed and trait-based energy budget model.

Soil microbiomes are highly diverse, and to improve their representation in biogeochemical models, microbial genome data can be leveraged to infer key functional traits. By integrating genome-inferred traits into a theory-based hierarchical framework, emergent behaviour arising from interactions of individual traits can be predicted. Here we combine theory-driven predictions of substrate uptake kinetics with a genome-informed trait-based dynamic energy budget model to predict emergent life-history traits and trade-offs in soil bacteria. When applied to a plant microbiome system, the model accurately predicted distinct substrate-acquisition strategies that aligned with observations, uncovering resource-dependent trade-offs between microbial growth rate and efficiency. For instance, inherently slower-growing microorganisms, favoured by organic acid exudation at later plant growth stages, exhibited enhanced carbon use efficiency (yield) without sacrificing growth rate (power). This insight has implications for retaining plant root-derived carbon in soils and highlights the power of data-driven, trait-based approaches for improving microbial representation in biogeochemical models.

Astrovirology: how viruses enhance our understanding of life in the Universe.

Viruses are the most numerically abundant biological entities on Earth. As ubiquitous replicators of molecular information and agents of community change, viruses have potent effects on the life on Earth, and may play a critical role in human spaceflight, for life-detection missions to other planetary bodies and planetary protection. However, major knowledge gaps constrain our understanding of the Earth's virosphere: (1) the role viruses play in biogeochemical cycles, (2) the origin(s) of viruses and (3) the involvement of viruses in the evolution, distribution and persistence of life. As viruses are the only replicators that span all known types of nucleic acids, an expanded experimental and theoretical toolbox built for Earth's viruses will be pivotal for detecting and understanding life on Earth and beyond. Only by filling in these knowledge and technical gaps we will obtain an inclusive assessment of how to distinguish and detect life on other planetary surfaces. Meanwhile, space exploration requires life-support systems for the needs of humans, plants and their microbial inhabitants. Viral effects on microbes and plants are essential for Earth's biosphere and human health, but virus-host interactions in spaceflight are poorly understood. Viral relationships with their hosts respond to environmental changes in complex ways which are difficult to predict by extrapolating from Earth-based proxies. These relationships should be studied in space to fully understand how spaceflight will modulate viral impacts on human health and life-support systems, including microbiomes. In this review, we address key questions that must be examined to incorporate viruses into Earth system models, life-support systems and life detection. Tackling these questions will benefit our efforts to develop planetary protection protocols and further our understanding of viruses in astrobiology.

Unraveling the functional dark matter through global metagenomics.

Metagenomes encode an enormous diversity of proteins, reflecting a multiplicity of functions and activities1,2. Exploration of this vast sequence space has been limited to a comparative analysis against reference microbial genomes and protein families derived from those genomes. Here, to examine the scale of yet untapped functional diversity beyond what is currently possible through the lens of reference genomes, we develop a computational approach to generate reference-free protein families from the sequence space in metagenomes. We analyse 26,931 metagenomes and identify 1.17 billion protein sequences longer than 35 amino acids with no similarity to any sequences from 102,491 reference genomes or the Pfam database3. Using massively parallel graph-based clustering, we group these proteins into 106,198 novel sequence clusters with more than 100 members, doubling the number of protein families obtained from the reference genomes clustered using the same approach. We annotate these families on the basis of their taxonomic, habitat, geographical and gene neighbourhood distributions and, where sufficient sequence diversity is available, predict protein three-dimensional models, revealing novel structures. Overall, our results uncover an enormously diverse functional space, highlighting the importance of further exploring the microbial functional dark matter.

Rhizosphere and detritusphere habitats modulate expression of soil N-cycling genes during plant development.

IMPORTANCE: Plant roots modulate microbial nitrogen (N) cycling by regulating the supply of root-derived carbon and nitrogen uptake. These differences in resource availability cause distinct micro-habitats to develop: soil near living roots, decaying roots, near both, or outside the direct influence of roots. While many environmental factors and genes control the microbial processes involved in the nitrogen cycle, most research has focused on single genes and pathways, neglecting the interactive effects these pathways have on each other. The processes controlled by these pathways determine consumption and production of N by soil microorganisms. We followed the expression of N-cycling genes in four soil microhabitats over a period of active root growth for an annual grass. We found that the presence of root litter and living roots significantly altered gene expression involved in multiple nitrogen pathways, as well as tradeoffs between pathways, which ultimately regulate N availability to plants.

Viral but not bacterial community successional patterns reflect extreme turnover shortly after rewetting dry soils.

As central members of soil trophic networks, viruses have the potential to drive substantial microbial mortality and nutrient turnover. Pinpointing viral contributions to terrestrial ecosystem processes remains a challenge, as temporal dynamics are difficult to unravel in the spatially and physicochemically heterogeneous soil environment. In Mediterranean grasslands, the first rainfall after seasonal drought provides an ecosystem reset, triggering microbial activity during a tractable window for capturing short-term dynamics. Here, we simulated precipitation in microcosms from four distinct dry grassland soils and generated 144 viromes, 84 metagenomes and 84 16S ribosomal RNA gene amplicon datasets to characterize viral, prokaryotic and relic DNA dynamics over 10 days. Vastly different viral communities in each soil followed remarkably similar successional trajectories. Wet-up triggered a significant increase in viral richness, followed by extensive compositional turnover. Temporal succession in prokaryotic communities was much less pronounced, perhaps suggesting differences in the scales of activity captured by viromes (representing recently produced, ephemeral viral particles) and total DNA. Still, differences in the relative abundances of Actinobacteria (enriched in dry soils) and Proteobacteria (enriched in wetted soils) matched those of their predicted phages, indicating viral predation of dominant bacterial taxa. Rewetting also rapidly depleted relic DNA, which subsequently reaccumulated, indicating substantial new microbial mortality in the days after wet-up, particularly of the taxa putatively under phage predation. Production of abundant, diverse viral particles via microbial host cell lysis appears to be a conserved feature of the early response to soil rewetting, and results suggest the potential for 'Cull-the-Winner' dynamics, whereby viruses infect and cull but do not decimate dominant host populations.

A subset of viruses thrives following microbial resuscitation during rewetting of a seasonally dry California grassland soil.

Viruses are abundant, ubiquitous members of soil communities that kill microbial cells, but how they respond to perturbation of soil ecosystems is essentially unknown. Here, we investigate lineage-specific virus-host dynamics in grassland soil following "wet-up", when resident microbes are both resuscitated and lysed after a prolonged dry period. Quantitative isotope tracing, time-resolved metagenomics and viromic analyses indicate that dry soil holds a diverse but low biomass reservoir of virions, of which only a subset thrives following wet-up. Viral richness decreases by 50% within 24 h post wet-up, while viral biomass increases four-fold within one week. Though recent hypotheses suggest lysogeny predominates in soil, our evidence indicates that viruses in lytic cycles dominate the response to wet-up. We estimate that viruses drive a measurable and continuous rate of cell lysis, with up to 46% of microbial death driven by viral lysis one week following wet-up. Thus, viruses contribute to turnover of soil microbial biomass and the widely reported CO2 efflux following wet-up of seasonally dry soils.

The predictive power of phylogeny on growth rates in soil bacterial communities.

Predicting ecosystem function is critical to assess and mitigate the impacts of climate change. Quantitative predictions of microbially mediated ecosystem processes are typically uninformed by microbial biodiversity. Yet new tools allow the measurement of taxon-specific traits within natural microbial communities. There is mounting evidence of a phylogenetic signal in these traits, which may support prediction and microbiome management frameworks. We investigated phylogeny-based trait prediction using bacterial growth rates from soil communities in Arctic, boreal, temperate, and tropical ecosystems. Here we show that phylogeny predicts growth rates of soil bacteria, explaining an average of 31%, and up to 58%, of the variation within ecosystems. Despite limited overlap in community composition across these ecosystems, shared nodes in the phylogeny enabled ancestral trait reconstruction and cross-ecosystem predictions. Phylogenetic relationships could explain up to 38% (averaging 14%) of the variation in growth rates across the highly disparate ecosystems studied. Our results suggest that shared evolutionary history contributes to similarity in the relative growth rates of related bacteria in the wild, allowing phylogeny-based predictions to explain a substantial amount of the variation in taxon-specific functional traits, within and across ecosystems.

Terabase-Scale Coassembly of a Tropical Soil Microbiome.

Petabases of environmental metagenomic data are publicly available, presenting an opportunity to characterize complex environments and discover novel lineages of life. Metagenome coassembly, in which many metagenomic samples from an environment are simultaneously analyzed to infer the underlying genomes' sequences, is an essential tool for achieving this goal. We applied MetaHipMer2, a distributed metagenome assembler that runs on supercomputing clusters, to coassemble 3.4 terabases (Tbp) of metagenome data from a tropical soil in the Luquillo Experimental Forest (LEF), Puerto Rico. The resulting coassembly yielded 39 high-quality (>90% complete, <5% contaminated, with predicted 23S, 16S, and 5S rRNA genes and ≥18 tRNAs) metagenome-assembled genomes (MAGs), including two from the candidate phylum Eremiobacterota. Another 268 medium-quality (≥50% complete, <10% contaminated) MAGs were extracted, including the candidate phyla Dependentiae, Dormibacterota, and Methylomirabilota. In total, 307 medium- or higher-quality MAGs were assigned to 23 phyla, compared to 294 MAGs assigned to nine phyla in the same samples individually assembled. The low-quality (<50% complete, <10% contaminated) MAGs from the coassembly revealed a 49% complete rare biosphere microbe from the candidate phylum FCPU426 among other low-abundance microbes, an 81% complete fungal genome from the phylum Ascomycota, and 30 partial eukaryotic MAGs with ≥10% completeness, possibly representing protist lineages. A total of 22,254 viruses, many of them low abundance, were identified. Estimation of metagenome coverage and diversity indicates that we may have characterized ≥87.5% of the sequence diversity in this humid tropical soil and indicates the value of future terabase-scale sequencing and coassembly of complex environments. IMPORTANCE Petabases of reads are being produced by environmental metagenome sequencing. An essential step in analyzing these data is metagenome assembly, the computational reconstruction of genome sequences from microbial communities. "Coassembly" of metagenomic sequence data, in which multiple samples are assembled together, enables more complete detection of microbial genomes in an environment than "multiassembly," in which samples are assembled individually. To demonstrate the potential for coassembling terabases of metagenome data to drive biological discovery, we applied MetaHipMer2, a distributed metagenome assembler that runs on supercomputing clusters, to coassemble 3.4 Tbp of reads from a humid tropical soil environment. The resulting coassembly, its functional annotation, and analysis are presented here. The coassembly yielded more, and phylogenetically more diverse, microbial, eukaryotic, and viral genomes than the multiassembly of the same data. Our resource may facilitate the discovery of novel microbial biology in tropical soils and demonstrates the value of terabase-scale metagenome sequencing.

Melanization slows the rapid movement of fungal necromass carbon and nitrogen into both bacterial and fungal decomposer communities and soils.

Microbial necromass contributes significantly to both soil carbon (C) persistence and ecosystem nitrogen (N) availability, but quantitative estimates of C and N movement from necromass into soils and decomposer communities are lacking. Additionally, while melanin is known to slow fungal necromass decomposition, how it influences microbial C and N acquisition as well as elemental release into soils remains unclear. Here, we tracked decomposition of isotopically labeled low and high melanin fungal necromass and measured 13C and 15N accumulation in surrounding soils and microbial communities over 77 d in a temperate forest in Minnesota, USA. Mass loss was significantly higher from low melanin necromass, corresponding with greater 13C and 15N soil inputs. A taxonomically and functionally diverse array of bacteria and fungi was enriched in 13C and/or 15N at all sampling points, with enrichment being consistently higher on low melanin necromass and earlier in decomposition. Similar patterns of preferential C and N enrichment of many bacterial and fungal genera early in decomposition suggest that both microbial groups co-contribute to the rapid assimilation of resource-rich soil organic matter inputs. While overall richness of taxa enriched in C was higher than in N for both bacteria and fungi, there was a significant positive relationship between C and N in co-enriched taxa. Collectively, our results demonstrate that melanization acts as a key ecological trait mediating not only fungal necromass decomposition rate but also necromass C and N release and that both elements are rapidly co-utilized by diverse bacterial and fungal decomposers in natural settings. IMPORTANCE Recent studies indicate that microbial dead cells, particularly those of fungi, play an important role in long-term carbon persistence in soils. Despite this growing recognition, how the resources within dead fungal cells (also known as fungal necromass) move into decomposer communities and soils are poorly quantified, particularly in studies based in natural environments. In this study, we found that the contribution of fungal necromass to soil carbon and nitrogen availability was slowed by the amount of melanin present in fungal cell walls. Further, despite the overall rapid acquisition of carbon and nitrogen from necromass by a diverse range of both bacteria and fungi, melanization also slowed microbial uptake of both elements. Collectively, our results indicate that melanization acts as a key ecological trait mediating not only fungal necromass decomposition rate, but also necromass carbon and nitrogen release into soil as well as microbial resource acquisition.

A standardized quantitative analysis strategy for stable isotope probing metagenomics.

Stable isotope probing (SIP) facilitates culture-independent identification of active microbial populations within complex ecosystems through isotopic enrichment of nucleic acids. Many DNA-SIP studies rely on 16S rRNA gene sequences to identify active taxa, but connecting these sequences to specific bacterial genomes is often challenging. Here, we describe a standardized laboratory and analysis framework to quantify isotopic enrichment on a per-genome basis using shotgun metagenomics instead of 16S rRNA gene sequencing. To develop this framework, we explored various sample processing and analysis approaches using a designed microbiome where the identity of labeled genomes and their level of isotopic enrichment were experimentally controlled. With this ground truth dataset, we empirically assessed the accuracy of different analytical models for identifying active taxa and examined how sequencing depth impacts the detection of isotopically labeled genomes. We also demonstrate that using synthetic DNA internal standards to measure absolute genome abundances in SIP density fractions improves estimates of isotopic enrichment. In addition, our study illustrates the utility of internal standards to reveal anomalies in sample handling that could negatively impact SIP metagenomic analyses if left undetected. Finally, we present SIPmg, an R package to facilitate the estimation of absolute abundances and perform statistical analyses for identifying labeled genomes within SIP metagenomic data. This experimentally validated analysis framework strengthens the foundation of DNA-SIP metagenomics as a tool for accurately measuring the in situ activity of environmental microbial populations and assessing their genomic potential. IMPORTANCE Answering the questions, "who is eating what?" and "who is active?" within complex microbial communities is paramount for our ability to model, predict, and modulate microbiomes for improved human and planetary health. These questions can be pursued using stable isotope probing to track the incorporation of labeled compounds into cellular DNA during microbial growth. However, with traditional stable isotope methods, it is challenging to establish links between an active microorganism's taxonomic identity and genome composition while providing quantitative estimates of the microorganism's isotope incorporation rate. Here, we report an experimental and analytical workflow that lays the foundation for improved detection of metabolically active microorganisms and better quantitative estimates of genome-resolved isotope incorporation, which can be used to further refine ecosystem-scale models for carbon and nutrient fluxes within microbiomes.

Expression of macromolecular organic nitrogen degrading enzymes identifies potential mediators of soil organic N availability to an annual grass.

Nitrogen (N) is frequently limiting to plant growth, in part because most soil N is present as polymeric organic compounds that are not readily taken up by plants. Microbial depolymerization of these large macromolecular N-substrates gradually releases available inorganic N. While many studies have researched and modeled controls on soil organic matter formation and bulk N mineralization, the ecological-spatial, temporal and phylogenetic-patterns underlying organic N degradation remain unclear. We analyzed 48 time-resolved metatranscriptomes and quantified N-depolymerization gene expression to resolve differential expression by soil habitat and time in specific taxonomic groups and gene-based guilds. We observed much higher expression of extracellular serine-type proteases than other extracellular N-degrading enzymes, with protease expression of predatory bacteria declining with time and other taxonomic patterns driven by the presence (Gammaproteobacteria) or absence (Thermoproteota) of live roots and root detritus (Deltaproteobacteria and Fungi). The primary chitinase chit1 gene was more highly expressed by eukaryotes near root detritus, suggesting predation of fungi. In some lineages, increased gene expression over time suggests increased competitiveness with rhizosphere age (Chloroflexi). Phylotypes from some genera had protease expression patterns that could benefit plant N nutrition, for example, we identified a Janthinobacterium phylotype and two Burkholderiales that depolymerize organic N near young roots and a Rhizobacter with elevated protease levels near mature roots. These taxon-resolved gene expression results provide an ecological read-out of microbial interactions and controls on N dynamics in specific soil microhabitats and could be used to target potential plant N bioaugmentation strategies.

Thermophilic Dehalococcoidia with unusual traits shed light on an unexpected past.

Although the phylum Chloroflexota is ubiquitous, its biology and evolution are poorly understood due to limited cultivability. Here, we isolated two motile, thermophilic bacteria from hot spring sediments belonging to the genus Tepidiforma and class Dehalococcoidia within the phylum Chloroflexota. A combination of cryo-electron tomography, exometabolomics, and cultivation experiments using stable isotopes of carbon revealed three unusual traits: flagellar motility, a peptidoglycan-containing cell envelope, and heterotrophic activity on aromatics and plant-associated compounds. Outside of this genus, flagellar motility has not been observed in Chloroflexota, and peptidoglycan-containing cell envelopes have not been described in Dehalococcoidia. Although these traits are unusual among cultivated Chloroflexota and Dehalococcoidia, ancestral character state reconstructions showed flagellar motility and peptidoglycan-containing cell envelopes were ancestral within the Dehalococcoidia, and subsequently lost prior to a major adaptive radiation of Dehalococcoidia into marine environments. However, despite the predominantly vertical evolutionary histories of flagellar motility and peptidoglycan biosynthesis, the evolution of enzymes for degradation of aromatics and plant-associated compounds was predominantly horizontal and complex. Together, the presence of these unusual traits in Dehalococcoidia and their evolutionary histories raise new questions about the timing and selective forces driving their successful niche expansion into global oceans.

Nutrients strengthen density dependence of per-capita growth and mortality rates in the soil bacterial community.

Density dependence in an ecological community has been observed in many macro-organismal ecosystems and is hypothesized to maintain biodiversity but is poorly understood in microbial ecosystems. Here, we analyze data from an experiment using quantitative stable isotope probing (qSIP) to estimate per-capita growth and mortality rates of bacterial populations in soils from several ecosystems along an elevation gradient which were subject to nutrient addition of either carbon alone (glucose; C) or carbon with nitrogen (glucose + ammonium-sulfate; C + N). Across all ecosystems, we found that higher population densities, quantified by the abundance of genomes per gram of soil, had lower per-capita growth rates in C + N-amended soils. Similarly, bacterial mortality rates in C + N-amended soils increased at a significantly higher rate with increasing population size than mortality rates in control and C-amended soils. In contrast to the hypothesis that density dependence would promote or maintain diversity, we observed significantly lower bacterial diversity in soils with stronger negative density-dependent growth. Here, density dependence was significantly but weakly responsive to nutrients and was not associated with higher bacterial diversity.

Life history strategies among soil bacteria-dichotomy for few, continuum for many.

Study of life history strategies may help predict the performance of microorganisms in nature by organizing the complexity of microbial communities into groups of organisms with similar strategies. Here, we tested the extent that one common application of life history theory, the copiotroph-oligotroph framework, could predict the relative population growth rate of bacterial taxa in soils from four different ecosystems. We measured the change of in situ relative growth rate to added glucose and ammonium using both 18O-H2O and 13C quantitative stable isotope probing to test whether bacterial taxa sorted into copiotrophic and oligotrophic groups. We saw considerable overlap in nutrient responses across most bacteria regardless of phyla, with many taxa growing slowly and few taxa that grew quickly. To define plausible life history boundaries based on in situ relative growth rates, we applied Gaussian mixture models to organisms' joint 18O-13C signatures and found that across experimental replicates, few taxa could consistently be assigned as copiotrophs, despite their potential for fast growth. When life history classifications were assigned based on average relative growth rate at varying taxonomic levels, finer resolutions (e.g., genus level) were significantly more effective in capturing changes in nutrient response than broad taxonomic resolution (e.g., phylum level). Our results demonstrate the difficulty in generalizing bacterial life history strategies to broad lineages, and even to single organisms across a range of soils and experimental conditions. We conclude that there is a continued need for the direct measurement of microbial communities in soil to advance ecologically realistic frameworks.

Initial soil organic carbon stocks govern changes in soil carbon: Reality or artifact?

Changes in soil organic carbon (SOC) storage have the potential to affect global climate; hence identifying environments with a high capacity to gain or lose SOC is of broad interest. Many cross-site studies have found that SOC-poor soils tend to gain or retain carbon more readily than SOC-rich soils. While this pattern may partly reflect reality, here we argue that it can also be created by a pair of statistical artifacts. First, soils that appear SOC-poor purely due to random variation will tend to yield more moderate SOC estimates upon resampling and hence will appear to accrue or retain more SOC than SOC-rich soils. This phenomenon is an example of regression to the mean. Second, normalized metrics of SOC change-such as relative rates and response ratios-will by definition show larger changes in SOC at lower initial SOC levels, even when the absolute change in SOC does not depend on initial SOC. These two artifacts create an exaggerated impression that initial SOC stocks are a major control on SOC dynamics. To address this problem, we recommend applying statistical corrections to eliminate the effect of regression to the mean, and avoiding normalized metrics when testing relationships between SOC change and initial SOC. Careful consideration of these issues in future cross-site studies will support clearer scientific inference that can better inform environmental management.

Diversifying and perennializing plants in agroecosystems alters retention of new C and N from crop residues.

Managing soils to retain new plant inputs is key to moving toward a sustainable and regenerative agriculture. Management practices, like diversifying and perennializing agroecosystems, may affect the decomposer organisms that regulate how new residue is converted to persistent soil organic matter. Here we tested whether 12 years of diversifying/perennializing plants in agroecosystems through extended rotations or grassland restoration would decrease losses of new plant residue inputs and, thus, increase retention of carbon (C) and nitrogen (N) in soil. We tracked dual-labeled (13 C and 15 N), isotopically enriched wheat (Triticum aestivum) residue in situ for 2 years as it decomposed in three agroecosystems: maize-soybean (CS) rotation, maize-soybean-wheat plus red clover and cereal rye cover crops (CSW2), and spring fallow management with regeneration of natural grassland species (seven to 10 species; SF). We measured losses of wheat residue (Cwheat and Nwheat ) in leached soil solution and greenhouse gas fluxes, as well as how much was recovered in microbial biomass and bulk soil at 5-cm increments down to 20 cm. CSW2 and SF both had unique, significant effects on residue decomposition and retention dynamics that were clear only when using nuanced metrics that able to tease apart subtle differences. For example, SF retained a greater portion of Cwheat in 0-5 cm surface soils (155%, p = 0.035) and narrowed the Cwheat to Nwheat ratio (p < 0.030) compared to CS. CSW2 increased an index of carbon-retention efficiency, Cwheat retained in the mesocosm divided by total measured, from 0.18 to 0.27 (49%, p = 0.001), compared to CS. Overall, we found that diversifying and extending the duration of living plants in agroecosystems can lead to greater retention of new residue inputs in subtle ways that require further investigation to fully understand.