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.

Shifts in bacterial traits under chronic nitrogen deposition align with soil processes in arbuscular, but not ectomycorrhizal-associated trees

Nitrogen (N) deposition increases soil carbon (C) storage by reducing microbial activity. These effects vary in soil beneath trees that associate with arbuscular (AM) and ectomycorrhizal (ECM) fungi. Variation in carbon C and N uptake traits among microbes may explain differences in soil nutrient cycling between mycorrhizal associations in response to high N loads, a mechanism not previously examined due to methodological limitations. Here, we used quantitative Stable Isotope Probing (qSIP) to measure bacterial C and N assimilation rates from an added organic compound, which we conceptualize as functional traits. As such, we applied a trait-based approach to explore whether variation in assimilation rates of bacterial taxa can inform shifts in soil function under chronic N deposition. We show taxon-specific and community-wide declines of bacterial C and N uptake under chronic N deposition in both AM and ECM soils. N deposition-induced reductions in microbial activity were mirrored by declines in soil organic matter mineralization rates in AM but not ECM soils. Our findings suggest C and N uptake traits of bacterial communities can predict C cycling feedbacks to N deposition in AM soils, but additional data, for instance on the traits of fungi, may be needed to connect microbial traits with soil C and N cycling in ECM systems. Our study also highlights the potential of employing qSIP in conjunction with trait-based analytical approaches to inform how ecological processes of microbial communities influence soil functioning.

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.

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.

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.

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.

HT-SIP: a semi-automated stable isotope probing pipeline identifies cross-kingdom interactions in the hyphosphere of arbuscular mycorrhizal fungi

BACKGROUND

Linking the identity of wild microbes with their ecophysiological traits and environmental functions is a key ambition for microbial ecologists. Of many techniques that strive for this goal, Stable-isotope probing-SIP-remains among the most comprehensive for studying whole microbial communities in situ. In DNA-SIP, actively growing microorganisms that take up an isotopically heavy substrate build heavier DNA, which can be partitioned by density into multiple fractions and sequenced. However, SIP is relatively low throughput and requires significant hands-on labor. We designed and tested a semi-automated, high-throughput SIP (HT-SIP) pipeline to support well-replicated, temporally resolved amplicon and metagenomics experiments. We applied this pipeline to a soil microhabitat with significant ecological importance-the hyphosphere zone surrounding arbuscular mycorrhizal fungal (AMF) hyphae. AMF form symbiotic relationships with most plant species and play key roles in terrestrial nutrient and carbon cycling.

RESULTS

Our HT-SIP pipeline for fractionation, cleanup, and nucleic acid quantification of density gradients requires one-sixth of the hands-on labor compared to manual SIP and allows 16 samples to be processed simultaneously. Automated density fractionation increased the reproducibility of SIP gradients compared to manual fractionation, and we show adding a non-ionic detergent to the gradient buffer improved SIP DNA recovery. We applied HT-SIP to ¹³C-AMF hyphosphere DNA from a ¹³CO₂ plant labeling study and created metagenome-assembled genomes (MAGs) using high-resolution SIP metagenomics (14 metagenomes per gradient). SIP confirmed the AMF Rhizophagus intraradices and associated MAGs were highly enriched (10-33 atom% ¹³C), even though the soils' overall enrichment was low (1.8 atom% ¹³C). We assembled 212 ¹³C-hyphosphere MAGs; the hyphosphere taxa that assimilated the most AMF-derived ¹³C were from the phyla Myxococcota, Fibrobacterota, Verrucomicrobiota, and the ammonia-oxidizing archaeon genus Nitrososphaera.

CONCLUSIONS

Our semi-automated HT-SIP approach decreases operator time and improves reproducibility by targeting the most labor-intensive steps of SIP-fraction collection and cleanup. We illustrate this approach in a unique and understudied soil microhabitat-generating MAGs of actively growing microbes living in the AMF hyphosphere (without plant roots). The MAGs' phylogenetic composition and gene content suggest predation, decomposition, and ammonia oxidation may be key processes in hyphosphere nutrient cycling. Video Abstract.

Coupling plant litter quantity to a novel metric for litter quality explains C storage changes in a thawing permafrost peatland

Permafrost thaw is a major potential feedback source to climate change as it can drive the increased release of greenhouse gases carbon dioxide (CO₂ ) and methane (CH₄ ). This carbon release from the decomposition of thawing soil organic material can be mitigated by increased net primary productivity (NPP) caused by warming, increasing atmospheric CO₂ , and plant community transition. However, the net effect on C storage also depends on how these plant community changes alter plant litter quantity, quality, and decomposition rates. Predicting decomposition rates based on litter quality remains challenging, but a promising new way forward is to incorporate measures of the energetic favorability to soil microbes of plant biomass decomposition. We asked how the variation in one such measure, the nominal oxidation state of carbon (NOSC), interacts with changing quantities of plant material inputs to influence the net C balance of a thawing permafrost peatland. We found: (1) Plant productivity (NPP) increased post-thaw, but instead of contributing to increased standing biomass, it increased plant biomass turnover via increased litter inputs to soil; (2) Plant litter thermodynamic favorability (NOSC) and decomposition rate both increased post-thaw, despite limited changes in bulk C:N ratios; (3) these increases caused the higher NPP to cycle more rapidly through both plants and soil, contributing to higher CO₂ and CH₄  fluxes from decomposition. Thus, the increased C-storage expected from higher productivity was limited and the high global warming potential of CH₄ contributed a net positive warming effect. Although post-thaw peatlands are currently C sinks due to high NPP offsetting high CO₂ release, this status is very sensitive to the plant community's litter input rate and quality. Integration of novel bioavailability metrics based on litter chemistry, including NOSC, into studies of ecosystem dynamics, is needed to improve the understanding of controls on arctic C stocks under continued ecosystem transition.

Decreased growth of wild soil microbes after 15 years of transplant-induced warming in a montane meadow

The carbon stored in soil exceeds that of plant biomass and atmospheric carbon and its stability can impact global climate. Growth of decomposer microorganisms mediates both the accrual and loss of soil carbon. Growth is sensitive to temperature and given the vast biological diversity of soil microorganisms, the response of decomposer growth rates to warming may be strongly idiosyncratic, varying among taxa, making ecosystem predictions difficult. Here, we show that 15 years of warming by transplanting plant-soil mesocosms down in elevation, strongly reduced the growth rates of soil microorganisms, measured in the field using undisturbed soil. The magnitude of the response to warming varied among microbial taxa. However, the direction of the response-reduced growth-was universal and warming explained twofold more variation than did the sum of taxonomic identity and its interaction with warming. For this ecosystem, most of the growth responses to warming could be explained without taxon-specific information, suggesting that in some cases microbial responses measured in aggregate may be adequate for climate modeling. Long-term experimental warming also reduced soil carbon content, likely a consequence of a warming-induced increase in decomposition, as warming-induced changes in plant productivity were negligible. The loss of soil carbon and decreased microbial biomass with warming may explain the reduced growth of the microbial community, more than the direct effects of temperature on growth. These findings show that direct and indirect effects of long-term warming can reduce growth rates of soil microbes, which may have important feedbacks to global warming.

Minnesota peat viromes reveal terrestrial and aquatic niche partitioning for local and global viral populations

BACKGROUND

Peatlands are expected to experience sustained yet fluctuating higher temperatures due to climate change, leading to increased microbial activity and greenhouse gas emissions. Despite mounting evidence for viral contributions to these processes in peatlands underlain with permafrost, little is known about viruses in other peatlands. More generally, soil viral biogeography and its potential drivers are poorly understood at both local and global scales. Here, 87 metagenomes and five viral size-fraction metagenomes (viromes) from a boreal peatland in northern Minnesota (the SPRUCE whole-ecosystem warming experiment and surrounding bog) were analyzed for dsDNA viral community ecological patterns, and the recovered viral populations (vOTUs) were compared with our curated PIGEON database of 266,125 vOTUs from diverse ecosystems.

RESULTS

Within the SPRUCE experiment, viral community composition was significantly correlated with peat depth, water content, and carbon chemistry, including CH₄ and CO₂ concentrations, but not with temperature during the first 2 years of warming treatments. Peat vOTUs with aquatic-like signatures (shared predicted protein content with marine and/or freshwater vOTUs) were significantly enriched in more waterlogged surface peat depths. Predicted host ranges for SPRUCE vOTUs were relatively narrow, generally within a single bacterial genus. Of the 4326 SPRUCE vOTUs, 164 were previously detected in other soils, mostly peatlands. None of the previously identified 202,371 marine and freshwater vOTUs in our PIGEON database were detected in SPRUCE peat, but 0.4% of 80,714 viral clusters (VCs, grouped by predicted protein content) were shared between soil and aquatic environments. On a per-sample basis, vOTU recovery was 32 times higher from viromes compared with total metagenomes.

CONCLUSIONS

Results suggest strong viral "species" boundaries between terrestrial and aquatic ecosystems and to some extent between peat and other soils, with differences less pronounced at higher taxonomic levels. The significant enrichment of aquatic-like vOTUs in more waterlogged peat suggests that viruses may also exhibit niche partitioning on more local scales. These patterns are presumably driven in part by host ecology, consistent with the predicted narrow host ranges. Although more samples and increased sequencing depth improved vOTU recovery from total metagenomes, the substantially higher per-sample vOTU recovery after viral particle enrichment highlights the utility of soil viromics. Video abstract The importance of Minnesota peat viromes in revealing terrestrial and aquatic niche partitioning for viral populations.

Active virus-host interactions at sub-freezing temperatures in Arctic peat soil

BACKGROUND

Winter carbon loss in northern ecosystems is estimated to be greater than the average growing season carbon uptake and is primarily driven by microbial decomposers. Viruses modulate microbial carbon cycling via induced mortality and metabolic controls, but it is unknown whether viruses are active under winter conditions (anoxic and sub-freezing temperatures).

RESULTS

We used stable isotope probing (SIP) targeted metagenomics to reveal the genomic potential of active soil microbial populations under simulated winter conditions, with an emphasis on viruses and virus-host dynamics. Arctic peat soils from the Bonanza Creek Long-Term Ecological Research site in Alaska were incubated under sub-freezing anoxic conditions with H218O or natural abundance water for 184 and 370 days. We sequenced 23 SIP-metagenomes and measured carbon dioxide (CO2) efflux throughout the experiment. We identified 46 bacterial populations (spanning 9 phyla) and 243 viral populations that actively took up 18O in soil and respired CO2 throughout the incubation. Active bacterial populations represented only a small portion of the detected microbial community and were capable of fermentation and organic matter degradation. In contrast, active viral populations represented a large portion of the detected viral community and one third were linked to active bacterial populations. We identified 86 auxiliary metabolic genes and other environmentally relevant genes. The majority of these genes were carried by active viral populations and had diverse functions such as carbon utilization and scavenging that could provide their host with a fitness advantage for utilizing much-needed carbon sources or acquiring essential nutrients.

CONCLUSIONS

Overall, there was a stark difference in the identity and function of the active bacterial and viral community compared to the unlabeled community that would have been overlooked with a non-targeted standard metagenomic analysis. Our results illustrate that substantial active virus-host interactions occur in sub-freezing anoxic conditions and highlight viruses as a major community-structuring agent that likely modulates carbon loss in peat soils during winter, which may be pivotal for understanding the future fate of arctic soils' vast carbon stocks. Video abstract.

The Functional Significance of Bacterial Predators

Predation structures food webs, influences energy flow, and alters rates and pathways of nutrient cycling through ecosystems, effects that are well documented for macroscopic predators. In the microbial world, predatory bacteria are common, yet little is known about their rates of growth and roles in energy flows through microbial food webs, in part because these are difficult to quantify. Here, we show that growth and carbon uptake were higher in predatory bacteria compared to nonpredatory bacteria, a finding across 15 sites, synthesizing 82 experiments and over 100,000 taxon-specific measurements of element flow into newly synthesized bacterial DNA. Obligate predatory bacteria grew 36% faster and assimilated carbon at rates 211% higher than nonpredatory bacteria. These differences were less pronounced for facultative predators (6% higher growth rates, 17% higher carbon assimilation rates), though high growth and carbon assimilation rates were observed for some facultative predators, such as members of the genera Lysobacter and Cytophaga, both capable of gliding motility and wolf-pack hunting behavior. Added carbon substrates disproportionately stimulated growth of obligate predators, with responses 63% higher than those of nonpredators for the Bdellovibrionales and 81% higher for the Vampirovibrionales, whereas responses of facultative predators to substrate addition were no different from those of nonpredators. This finding supports the ecological theory that higher productivity increases predator control of lower trophic levels. These findings also indicate that the functional significance of bacterial predators increases with energy flow and that predatory bacteria influence element flow through microbial food webs.IMPORTANCE The word "predator" may conjure images of leopards killing and eating impala on the African savannah or of great white sharks attacking elephant seals off the coast of California. But microorganisms are also predators, including bacteria that kill and eat other bacteria. While predatory bacteria have been found in many environments, it has been challenging to document their importance in nature. This study quantified the growth of predatory and nonpredatory bacteria in soils (and one stream) by tracking isotopically labeled substrates into newly synthesized DNA. Predatory bacteria were more active than nonpredators, and obligate predators, such as Bdellovibrionales and Vampirovibrionales, increased in growth rate in response to added substrates at the base of the food chain, strong evidence of trophic control. This work provides quantitative measures of predator activity and suggests that predatory bacteria-along with protists, nematodes, and phages-are active and important in microbial food webs.

The temperature sensitivity of soil: microbial biodiversity, growth, and carbon mineralization

Microorganisms drive soil carbon mineralization and changes in their activity with increased temperature could feedback to climate change. Variation in microbial biodiversity and the temperature sensitivities (Q₁₀) of individual taxa may explain differences in the Q₁₀ of soil respiration, a possibility not previously examined due to methodological limitations. Here, we show phylogenetic and taxonomic variation in the Q₁₀ of growth (5-35 °C) among soil bacteria from four sites, one from each of Arctic, boreal, temperate, and tropical biomes. Differences in the temperature sensitivities of taxa and the taxonomic composition of communities determined community-assembled bacterial growth Q₁₀, which was strongly predictive of soil respiration Q₁₀ within and across biomes. Our results suggest community-assembled traits of microbial taxa may enable enhanced prediction of carbon cycling feedbacks to climate change in ecosystems across the globe.

Measurement Error and Resolution in Quantitative Stable Isotope Probing: Implications for Experimental Design

Quantitative stable isotope probing (qSIP) estimates isotope tracer incorporation into DNA of individual microbes and can link microbial biodiversity and biogeochemistry in complex communities. As with any quantitative estimation technique, qSIP involves measurement error, and a fuller understanding of error, precision, and statistical power benefits qSIP experimental design and data interpretation. We used several qSIP data sets-from soil and seawater microbiomes-to evaluate how variance in isotope incorporation estimates depends on organism abundance and resolution of the density fractionation scheme. We assessed statistical power for replicated qSIP studies, plus sensitivity and specificity for unreplicated designs. As a taxon's abundance increases, the variance of its weighted mean density declines. Nine fractions appear to be a reasonable trade-off between cost and precision for most qSIP applications. Increasing the number of density fractions beyond that reduces variance, although the magnitude of this benefit declines with additional fractions. Our analysis suggests that, if a taxon has an isotope enrichment of 10 atom% excess, there is a 60% chance that this will be detected as significantly different from zero (with alpha 0.1). With five replicates, isotope enrichment of 5 atom% could be detected with power (0.6) and alpha (0.1). Finally, we illustrate the importance of internal standards, which can help to calibrate per sample conversions of %GC to mean weighted density. These results should benefit researchers designing future SIP experiments and provide a useful reference for metagenomic SIP applications where both financial and computational limitations constrain experimental scope.IMPORTANCE One of the biggest challenges in microbial ecology is correlating the identity of microorganisms with the roles they fulfill in natural environmental systems. Studies of microbes in pure culture reveal much about their genomic content and potential functions but may not reflect an organism's activity within its natural community. Culture-independent studies supply a community-wide view of composition and function in the context of community interactions but often fail to link the two. Quantitative stable isotope probing (qSIP) is a method that can link the identity and functional activity of specific microbes within a naturally occurring community. Here, we explore how the resolution of density gradient fractionation affects the error and precision of qSIP results, how they may be improved via additional experimental replication, and discuss cost-benefit balanced scenarios for SIP experimental design.

Taxon-specific microbial growth and mortality patterns reveal distinct temporal population responses to rewetting in a California grassland soil

Microbial activity increases after rewetting dry soil, resulting in a pulse of carbon mineralization and nutrient availability. The biogeochemical responses to wet-up are reasonably well understood and known to be microbially mediated. Yet, the population level dynamics, and the resulting changes in microbial community patterns, are not well understood as ecological phenomena. Here, we used sequencing of 16S rRNA genes coupled with heavy water (H218O) DNA quantitative stable isotope probing to estimate population-specific rates of growth and mortality in response to a simulated wet-up event in a California annual grassland soil. Bacterial growth and mortality responded rapidly to wet-up, within 3 h, and continued throughout the 168 h incubation, with patterns of sequential growth observed at the phylum level. Of the 37 phyla detected in the prewet community, growth was found in 18 phyla while mortality was measured in 26 phyla. Rapid growth and mortality rates were measurable within 3 h of wet-up but had contrasting characteristics; growth at 3 h was dominated by select taxa within the Proteobacteria and Firmicutes, whereas mortality was taxonomically widespread. Furthermore, across the community, mortality exhibited density-independence, consistent with the indiscriminate shock resulting from dry-down and wet-up, whereas growth was density-dependent, consistent with control by competition or predation. Total aggregated growth across the community was highly correlated with total soil CO2 production. Together, these results illustrate how previously "invisible" population responses can translate quantitatively to emergent observations of ecosystem-scale biogeochemistry.

Predictive genomic traits for bacterial growth in culture versus actual growth in soil

Relationships between microbial genes and performance are often evaluated in the laboratory in pure cultures, with little validation in nature. Here, we show that genomic traits related to laboratory measurements of maximum growth potential failed to predict the growth rates of bacteria in unamended soil, but successfully predicted growth responses to resource pulses: growth increased with 16S rRNA gene copy number and declined with genome size after substrate addition to soils, responses that were repeated in four different ecosystems. Genome size best predicted growth rate in response to addition of glucose alone; adding ammonium with glucose weakened the relationship, and the relationship was absent in nutrient-replete pure cultures, consistent with the idea that reduced genome size is a mechanism of nutrient conservation. Our findings demonstrate that genomic traits of soil bacteria can map to their ecological performance in nature, but the mapping is poor under native soil conditions, where genomic traits related to stress tolerance may prove more predictive. These results remind that phenotype depends on environmental context, underscoring the importance of verifying proposed schemes of trait-based strategies through direct measurement of performance in nature, an important and currently missing foundation for translating microbial processes from genes to ecosystems.

Stable isotope informed genome-resolved metagenomics reveals that Saccharibacteria utilize microbially-processed plant-derived carbon

BACKGROUND

The transformation of plant photosynthate into soil organic carbon and its recycling to CO₂ by soil microorganisms is one of the central components of the terrestrial carbon cycle. There are currently large knowledge gaps related to which soil-associated microorganisms take up plant carbon in the rhizosphere and the fate of that carbon.

RESULTS

We conducted an experiment in which common wild oats (Avena fatua) were grown in a ¹³CO₂ atmosphere and the rhizosphere and non-rhizosphere soil was sampled for genomic analyses. Density gradient centrifugation of DNA extracted from soil samples enabled distinction of microbes that did and did not incorporate the ¹³C into their DNA. A 1.45-Mbp genome of a Saccharibacteria (TM7) was identified and, despite the microbial complexity of rhizosphere soil, curated to completion. The genome lacks many biosynthetic pathways, including genes required to synthesize DNA de novo. Rather, it requires externally derived nucleotides for DNA and RNA synthesis. Given this, we conclude that rhizosphere-associated Saccharibacteria recycle DNA from bacteria that live off plant exudates and/or phage that acquired ¹³C because they preyed upon these bacteria and/or directly from the labeled plant DNA. Isotopic labeling indicates that the population was replicating during the 6-week period of plant growth. Interestingly, the genome is ~ 30% larger than other complete Saccharibacteria genomes from non-soil environments, largely due to more genes for complex carbon utilization and amino acid metabolism. Given the ability to degrade cellulose, hemicellulose, pectin, starch, and 1,3-β-glucan, we predict that this Saccharibacteria generates energy by fermentation of soil necromass and plant root exudates to acetate and lactate. The genome also encodes a linear electron transport chain featuring a terminal oxidase, suggesting that this Saccharibacteria may respire aerobically. The genome encodes a hydrolase that could breakdown salicylic acid, a plant defense signaling molecule, and genes to interconvert a variety of isoprenoids, including the plant hormone zeatin.

CONCLUSIONS

Rhizosphere Saccharibacteria likely depend on other bacteria for basic cellular building blocks. We propose that isotopically labeled CO₂ is incorporated into plant-derived carbon and then into the DNA of rhizosphere organisms capable of nucleotide synthesis, and the nucleotides are recycled into Saccharibacterial genomes.

Depth matters: effects of precipitation regime on soil microbial activity upon rewetting of a plant-soil system

Changes in frequency and amplitude of rain events, that is, precipitation patterns, result in different water conditions with soil depth, and likely affect plant growth and shape plant and soil microbial activity. Here, we used ¹⁸O stable isotope probing (SIP) to investigate bacterial and fungal communities that actively grew or not upon rewetting, at three different depths in soil mesocosms previously subjected to frequent or infrequent watering for 12 weeks (equal total water input). Phylogenetic marker genes for bacteria and fungi were sequenced after rewetting, and plant-soil microbial coupling documented by plant ¹³C-CO₂ labeling. Soil depth, rather than precipitation pattern, was most influential in shaping microbial response to rewetting, and had differential effects on active and inactive bacterial and fungal communities. After rewetting, active bacterial communities were less rich, more even and phylogenetically related than the inactive, and reactivated throughout the soil profile. Active fungal communities after rewetting were less abundant and rich than the inactive. The coupling between plants and soil microbes decreased under infrequent watering in the top soil layer. We suggest that differences in fungal and bacterial abundance and relative activity could result in large effects on subsequent soil biogeochemical cycling.

Growth and death of bacteria and fungi underlie rainfall-induced carbon dioxide pulses from seasonally dried soil

The rapid increase in microbial activity that occurs when a dry soil is rewetted has been well documented and is of great interest due to implications of changing precipitation patterns on soil C dynamics. Several studies have shown minor net changes in microbial population diversity or abundance following wet-up, but the gross population dynamics of bacteria and fungi resulting from soil wet-up are virtually unknown. Here we applied DNA stable isotope probing with H218O coupled with quantitative PCR to characterize new growth, survival, and mortality of bacteria and fungi following the rewetting of a seasonally dried California annual grassland soil. Microbial activity, as determined by CO2 production, increased significantly within three hours of wet-up, yet new growth was not detected until after three hours, suggesting a pulse of nongrowth activity immediately following wet-up, likely due to osmo-regulation and resuscitation from dormancy in response to the rapid change in water potential. Total microbial abundance revealed little change throughout the seven-day post-wet incubation, but there was substantial turnover of both bacterial and fungal populations (49% and 52%, respectively). New growth was linear between 24 and 168 hours for both bacteria and fungi, with average growth rates of 2.3 x 10(8) bacterial 16S rRNA gene copies x [g dry mass](-1) x h(-1) and 4.3 x 10(7) fungal ITS copies x [g dry mass](-1) x h(-1). While bacteria and fungi differed in their mortality and survival characteristics during the seven-day incubation, mortality that occurred within the first three hours was similar, with 25% and 27% of bacterial and fungal gene copies disappearing from the pre-wet community, respectively. The rapid disappearance of gene copies indicates that cell death, occurring either during the extreme dry down period (preceding five months) or during the rapid change in water potential due to wet-up, generates a significant pool of available C that likely contributes to the large pulse in CO2 associated with wet-up. A dynamic assemblage of growing and dying organisms controlled the CO2 pulse, but the balance between death and growth resulted in relatively stable total population abundances, even after a profound and sudden change in environment.

Evaluating rRNA as an indicator of microbial activity in environmental communities: limitations and uses

Microbes exist in a range of metabolic states (for example, dormant, active and growing) and analysis of ribosomal RNA (rRNA) is frequently employed to identify the 'active' fraction of microbes in environmental samples. While rRNA analyses are no longer commonly used to quantify a population's growth rate in mixed communities, due to rRNA concentration not scaling linearly with growth rate uniformly across taxa, rRNA analyses are still frequently used toward the more conservative goal of identifying populations that are currently active in a mixed community. Yet, evidence indicates that the general use of rRNA as a reliable indicator of metabolic state in microbial assemblages has serious limitations. This report highlights the complex and often contradictory relationships between rRNA, growth and activity. Potential mechanisms for confounding rRNA patterns are discussed, including differences in life histories, life strategies and non-growth activities. Ways in which rRNA data can be used for useful characterization of microbial assemblages are presented, along with questions to be addressed in future studies.

Abundance of microbial genes associated with nitrogen cycling as indices of biogeochemical process rates across a vegetation gradient in Alaska

Nitrification and denitrification processes are crucial to plant nutrient availability, eutrophication and greenhouse gas production both locally and globally. Unravelling the major environmental predictors for nitrification and denitrification is thus pivotal in order to understand and model environmental nitrogen (N) cycling. Here, we sampled five plant community types characteristic of interior Alaska, including black spruce, bog birch, tussock grass and two fens. We assessed abundance of functional genes affiliated with nitrification (bacterial and archaeal amoA) and denitrification (nirK/S and nosZ) using qPCR, soil characteristics, potential nitrification and denitrification rates (PNR and PDR) and gross mineralization rates. The main chemical and biological predictors for PNR and PDR were assigned through path analysis. The potential N cycling rates varied dramatically between sites, from some of the highest (in fens) to some of the lowest (in black spruce) measured globally. Based on path analysis, functional gene abundances were the most important variables to predict potential rates. PNR was best explained by bacterial amoA gene abundance followed by ammonium content, whereas PDR was best explained directly by nosZ gene abundance and indirectly by nirK/S gene abundance and nitrate. Hence, functional gene abundance is a valuable index that integrates recent environmental history and recent process activity, and therefore is a good predictor of potential rates. The results of this study contribute to our understanding of the relative importance of different biological and chemical factors in driving the potential for nitrification and denitrification across terrestrial ecosystems.

Metagenomic analysis of a permafrost microbial community reveals a rapid response to thaw

Permafrost contains an estimated 1672 Pg carbon (C), an amount roughly equivalent to the total currently contained within land plants and the atmosphere. This reservoir of C is vulnerable to decomposition as rising global temperatures cause the permafrost to thaw. During thaw, trapped organic matter may become more accessible for microbial degradation and result in greenhouse gas emissions. Despite recent advances in the use of molecular tools to study permafrost microbial communities, their response to thaw remains unclear. Here we use deep metagenomic sequencing to determine the impact of thaw on microbial phylogenetic and functional genes, and relate these data to measurements of methane emissions. Metagenomics, the direct sequencing of DNA from the environment, allows the examination of whole biochemical pathways and associated processes, as opposed to individual pieces of the metabolic puzzle. Our metagenome analyses reveal that during transition from a frozen to a thawed state there are rapid shifts in many microbial, phylogenetic and functional gene abundances and pathways. After one week of incubation at 5 °C, permafrost metagenomes converge to be more similar to each other than while they are frozen. We find that multiple genes involved in cycling of C and nitrogen shift rapidly during thaw. We also construct the first draft genome from a complex soil metagenome, which corresponds to a novel methanogen. Methane previously accumulated in permafrost is released during thaw and subsequently consumed by methanotrophic bacteria. Together these data point towards the importance of rapid cycling of methane and nitrogen in thawing permafrost.

Dynamics of ¹⁸O incorporation from H₂ ¹⁸O into soil microbial DNA

Heavy water (H(2)(18)O) has been used to label DNA of soil microorganisms in stable isotope probing experiments, yet no measurements have been reported for the (18)O content of DNA from soil incubated with heavy water. Here we present the first measurements of atom% (18)O for DNA extracted from soil incubated with the addition of H(2)(18)O. Four experiments were conducted to test how the atom% (18)O of DNA, extracted from Ponderosa Pine forest soil incubated with heavy water, was affected by the following variables: (1) time, (2) nutrients, (3) soil moisture, and (4) atom% (18)O of added H(2)O. In the time series experiment, the atom% (18)O of DNA increased linearly (R(2) = 0.994, p < 0.01) over the first 72 h of incubation. In the nutrient addition experiment, there was a positive correlation (R(2) = 0.991, p = 0.006) between the log(10) of the amount of tryptic soy broth, a complex nutrient broth, added to soil and the log(10) of the atom% (18)O of DNA. For the experiment where soil moisture was manipulated, the atom% (18)O of DNA increased with higher soil moisture until soil moisture reached 30%, above which (18)O enrichment of DNA declined as soils became more saturated. When the atom% (18)O for H(2)O added was varied, there was a positive linear relationship between the atom% (18)O of the added water and the atom% (18)O of the DNA. Results indicate that quantification of (18)O incorporated into DNA from H(2)(18)O has potential to be used as a proxy for microbial growth in soil.