Genomic insights into redox-driven microbial processes for carbon decomposition in thawing Arctic soils and permafrost.

Climate change is rapidly transforming Arctic landscapes where increasing soil temperatures speed up permafrost thaw. This exposes large carbon stocks to microbial decomposition, possibly worsening climate change by releasing more greenhouse gases. Understanding how microbes break down soil carbon, especially under the anaerobic conditions of thawing permafrost, is important to determine future changes. Here, we studied the microbial community dynamics and soil carbon decomposition potential in permafrost and active layer soils under anaerobic laboratory conditions that simulated an Arctic summer thaw. The microbial and viral compositions in the samples were analyzed based on metagenomes, metagenome-assembled genomes, and metagenomic viral contigs (mVCs). Following the thawing of permafrost, there was a notable shift in microbial community structure, with fermentative Firmicutes and Bacteroidota taking over from Actinobacteria and Proteobacteria over the 60-day incubation period. The increase in iron and sulfate-reducing microbes had a significant role in limiting methane production from thawed permafrost, underscoring the competition within microbial communities. We explored the growth strategies of microbial communities and found that slow growth was the major strategy in both the active layer and permafrost. Our findings challenge the assumption that fast-growing microbes mainly respond to environmental changes like permafrost thaw. Instead, they indicate a common strategy of slow growth among microbial communities, likely due to the thermodynamic constraints of soil substrates and electron acceptors, and the need for microbes to adjust to post-thaw conditions. The mVCs harbored a wide range of auxiliary metabolic genes that may support cell protection from ice formation in virus-infected cells. IMPORTANCE: As the Arctic warms, thawing permafrost unlocks carbon, potentially accelerating climate change by releasing greenhouse gases. Our research delves into the underlying biogeochemical processes likely mediated by the soil microbial community in response to the wet and anaerobic conditions, akin to an Arctic summer thaw. We observed a significant shift in the microbial community post-thaw, with fermentative bacteria like Firmicutes and Bacteroidota taking over and switching to different fermentation pathways. The dominance of iron and sulfate-reducing bacteria likely constrained methane production in the thawing permafrost. Slow-growing microbes outweighed fast-growing ones, even after thaw, upending the expectation that rapid microbial responses to dominate after permafrost thaws. This research highlights the nuanced and complex interactions within Arctic soil microbial communities and underscores the challenges in predicting microbial response to environmental change.

MetFish: a Metabolomics Pipeline for Studying Microbial Communities in Chemically Extreme Environments.

Metabolites have essential roles in microbial communities, including as mediators of nutrient and energy exchange, cell-to-cell communication, and antibiosis. However, detecting and quantifying metabolites and other chemicals in samples having extremes in salt or mineral content using liquid chromatography-mass spectrometry (LC-MS)-based methods remains a significant challenge. Here, we report a facile method based on in situ chemical derivatization followed by extraction for analysis of metabolites and other chemicals in hypersaline samples, enabling for the first time direct LC-MS-based exometabolomics analysis in sample matrices containing up to 2 M total dissolved salts. The method, MetFish, is applicable to molecules containing amine, carboxylic acid, carbonyl, or hydroxyl functional groups, and it can be integrated into either targeted or untargeted analysis pipelines. In targeted analyses, MetFish provided limits of quantification as low as 1 nM, broad linear dynamic ranges (up to 5 to 6 orders of magnitude) with excellent linearity, and low median interday reproducibility (e.g., 2.6%). MetFish was successfully applied in targeted and untargeted exometabolomics analyses of microbial consortia, quantifying amino acid dynamics in the exometabolome during community succession; in situ in a native prairie soil, whose exometabolome was isolated using a hypersaline extraction; and in input and produced fluids from a hydraulically fractured well, identifying dramatic changes in the exometabolome over time in the well. IMPORTANCE The identification and accurate quantification of metabolites using electrospray ionization-mass spectrometry (ESI-MS) in hypersaline samples is a challenge due to matrix effects. Clean-up and desalting strategies that typically work well for samples with lower salt concentrations are often ineffective in hypersaline samples. To address this gap, we developed and demonstrated a simple yet sensitive and accurate method-MetFish-using chemical derivatization to enable mass spectrometry-based metabolomics in a variety of hypersaline samples from varied ecosystems and containing up to 2 M dissolved salts.

Temporal, Spatial, and Temperature Controls on Organic Carbon Mineralization and Methanogenesis in Arctic High-Centered Polygon Soils.

Warming temperatures in continuous permafrost zones of the Arctic will alter both hydrological and geochemical soil conditions, which are strongly linked with heterotrophic microbial carbon (C) cycling. Heterogeneous permafrost landscapes are often dominated by polygonal features formed by expanding ice wedges: water accumulates in low centered polygons (LCPs), and water drains outward to surrounding troughs in high centered polygons (HCPs). These geospatial differences in hydrology cause gradients in biogeochemistry, soil C storage potential, and thermal properties. Presently, data quantifying carbon dioxide (CO2) and methane (CH4) release from HCP soils are needed to support modeling and evaluation of warming-induced CO2 and CH4 fluxes from tundra soils. This study quantifies the distribution of microbial CO2 and CH4 release in HCPs over a range of temperatures and draws comparisons to previous LCP studies. Arctic tundra soils were initially characterized for geochemical and hydraulic properties. Laboratory incubations at -2, +4, and +8°C were used to quantify temporal trends in CO2 and CH4 production from homogenized active layer organic and mineral soils in HCP centers and troughs, and methanogen abundance was estimated from mcrA gene measurements. Results showed that soil water availability, organic C, and redox conditions influence temporal dynamics and magnitude of gas production from HCP active layer soils during warming. At early incubation times (2-9 days), higher CO2 emissions were observed from HCP trough soils than from HCP center soils, but increased CO2 production occurred in center soils at later times (>20 days). HCP center soils did not support methanogenesis, but CH4-producing trough soils did indicate methanogen presence. Consistent with previous LCP studies, HCP organic soils showed increased CO2 and CH4 production with elevated water content, but HCP trough mineral soils produced more CH4 than LCP mineral soils. HCP mineral soils also released substantial CO2 but did not show a strong trend in CO2 and CH4 release with water content. Knowledge of temporal and spatial variability in microbial C mineralization rates of Arctic soils in response to warming are key to constraining uncertainties in predictive climate models.

Metaphenomic Responses of a Native Prairie Soil Microbiome to Moisture Perturbations.

Climate change is causing shifts in precipitation patterns in the central grasslands of the United States, with largely unknown consequences on the collective physiological responses of the soil microbial community, i.e., the metaphenome. Here, we used an untargeted omics approach to determine the soil microbial community's metaphenomic response to soil moisture and to define specific metabolic signatures of the response. Specifically, we aimed to develop the technical approaches and metabolic mapping framework necessary for future systematic ecological studies. We collected soil from three locations at the Konza Long-Term Ecological Research (LTER) field station in Kansas, and the soils were incubated for 15 days under dry or wet conditions and compared to field-moist controls. The microbiome response to wetting or drying was determined by 16S rRNA amplicon sequencing, metatranscriptomics, and metabolomics, and the resulting shifts in taxa, gene expression, and metabolites were assessed. Soil drying resulted in significant shifts in both the composition and function of the soil microbiome. In contrast, there were few changes following wetting. The combined metabolic and metatranscriptomic data were used to generate reaction networks to determine the metaphenomic response to soil moisture transitions. Site location was a strong determinant of the response of the soil microbiome to moisture perturbations. However, some specific metabolic pathways changed consistently across sites, including an increase in pathways and metabolites for production of sugars and other osmolytes as a response to drying. Using this approach, we demonstrate that despite the high complexity of the soil habitat, it is possible to generate insight into the effect of environmental change on the soil microbiome and its physiology and functions, thus laying the groundwork for future, targeted studies.IMPORTANCE Climate change is predicted to result in increased drought extent and intensity in the highly productive, former tallgrass prairie region of the continental United States. These soils store large reserves of carbon. The decrease in soil moisture due to drought has largely unknown consequences on soil carbon cycling and other key biogeochemical cycles carried out by soil microbiomes. In this study, we found that soil drying had a significant impact on the structure and function of soil microbial communities, including shifts in expression of specific metabolic pathways, such as those leading toward production of osmoprotectant compounds. This study demonstrates the application of an untargeted multi-omics approach to decipher details of the soil microbial community's metaphenotypic response to environmental perturbations and should be applicable to studies of other complex microbial systems as well.

The MPLEx Protocol for Multi-omic Analyses of Soil Samples.

Mass spectrometry (MS)-based integrated metaproteomic, metabolomic, and lipidomic (multi-omic) studies are transforming our ability to understand and characterize microbial communities in environmental and biological systems. These measurements are even enabling enhanced analyses of complex soil microbial communities, which are the most complex microbial systems known to date. Multi-omic analyses, however, do have sample preparation challenges, since separate extractions are typically needed for each omic study, thereby greatly amplifying the preparation time and amount of sample required. To address this limitation, a 3-in-1 method for the simultaneous extraction of metabolites, proteins, and lipids (MPLEx) from the same soil sample was created by adapting a solvent-based approach. This MPLEx protocol has proven to be both simple and robust for many sample types, even when utilized for limited quantities of complex soil samples. The MPLEx method also greatly enabled the rapid multi-omic measurements needed to gain a better understanding of the members of each microbial community, while evaluating the changes taking place upon biological and environmental perturbations.

Temporal Dynamics of In-Field Bioreactor Populations Reflect the Groundwater System and Respond Predictably to Perturbation.

Temporal variability complicates testing the influences of environmental variability on microbial community structure and thus function. An in-field bioreactor system was developed to assess oxic versus anoxic manipulations on in situ groundwater communities. Each sample was sequenced (16S SSU rRNA genes, average 10,000 reads), and biogeochemical parameters are monitored by quantifying 53 metals, 12 organic acids, 14 anions, and 3 sugars. Changes in dissolved oxygen (DO), pH, and other variables were similar across bioreactors. Sequencing revealed a complex community that fluctuated in-step with the groundwater community and responded to DO. This also directly influenced the pH, and so the biotic impacts of DO and pH shifts are correlated. A null model demonstrated that bioreactor communities were driven in part not only by experimental conditions but also by stochastic variability and did not accurately capture alterations in diversity during perturbations. We identified two groups of abundant OTUs important to this system; one was abundant in high DO and pH and contained heterotrophs and oxidizers of iron, nitrite, and ammonium, whereas the other was abundant in low DO with the capability to reduce nitrate. In-field bioreactors are a powerful tool for capturing natural microbial community responses to alterations in geochemical factors beyond the bulk phase.

Moleculo Long-Read Sequencing Facilitates Assembly and Genomic Binning from Complex Soil Metagenomes.

Soil metagenomics has been touted as the "grand challenge" for metagenomics, as the high microbial diversity and spatial heterogeneity of soils make them unamenable to current assembly platforms. Here, we aimed to improve soil metagenomic sequence assembly by applying the Moleculo synthetic long-read sequencing technology. In total, we obtained 267 Gbp of raw sequence data from a native prairie soil; these data included 109.7 Gbp of short-read data (~100 bp) from the Joint Genome Institute (JGI), an additional 87.7 Gbp of rapid-mode read data (~250 bp), plus 69.6 Gbp (>1.5 kbp) from Moleculo sequencing. The Moleculo data alone yielded over 5,600 reads of >10 kbp in length, and over 95% of the unassembled reads mapped to contigs of >1.5 kbp. Hybrid assembly of all data resulted in more than 10,000 contigs over 10 kbp in length. We mapped three replicate metatranscriptomes derived from the same parent soil to the Moleculo subassembly and found that 95% of the predicted genes, based on their assignments to Enzyme Commission (EC) numbers, were expressed. The Moleculo subassembly also enabled binning of >100 microbial genome bins. We obtained via direct binning the first complete genome, that of "Candidatus Pseudomonas sp. strain JKJ-1" from a native soil metagenome. By mapping metatranscriptome sequence reads back to the bins, we found that several bins corresponding to low-relative-abundance Acidobacteria were highly transcriptionally active, whereas bins corresponding to high-relative-abundance Verrucomicrobia were not. These results demonstrate that Moleculo sequencing provides a significant advance for resolving complex soil microbial communities. IMPORTANCE Soil microorganisms carry out key processes for life on our planet, including cycling of carbon and other nutrients and supporting growth of plants. However, there is poor molecular-level understanding of their functional roles in ecosystem stability and responses to environmental perturbations. This knowledge gap is largely due to the difficulty in culturing the majority of soil microbes. Thus, use of culture-independent approaches, such as metagenomics, promises the direct assessment of the functional potential of soil microbiomes. Soil is, however, a challenge for metagenomic assembly due to its high microbial diversity and variable evenness, resulting in low coverage and uneven sampling of microbial genomes. Despite increasingly large soil metagenome data volumes (>200 Gbp), the majority of the data do not assemble. Here, we used the cutting-edge approach of synthetic long-read sequencing technology (Moleculo) to assemble soil metagenome sequence data into long contigs and used the assemblies for binning of genomes. Author Video: An author video summary of this article is available.

Indexing Permafrost Soil Organic Matter Degradation Using High-Resolution Mass Spectrometry.

Microbial degradation of soil organic matter (SOM) is a key process for terrestrial carbon cycling, although the molecular details of these transformations remain unclear. This study reports the application of ultrahigh resolution mass spectrometry to profile the molecular composition of SOM and its degradation during a simulated warming experiment. A soil sample, collected near Barrow, Alaska, USA, was subjected to a 40-day incubation under anoxic conditions and analyzed before and after the incubation to determine changes of SOM composition. A CHO index based on molecular C, H, and O data was utilized to codify SOM components according to their observed degradation potentials. Compounds with a CHO index score between -1 and 0 in a water-soluble fraction (WSF) demonstrated high degradation potential, with a highest shift of CHO index occurred in the N-containing group of compounds, while similar stoichiometries in a base-soluble fraction (BSF) did not. Additionally, compared with the classical H:C vs O:C van Krevelen diagram, CHO index allowed for direct visualization of the distribution of heteroatoms such as N in the identified SOM compounds. We demonstrate that CHO index is useful not only in characterizing arctic SOM at the molecular level but also enabling quantitative description of SOM degradation, thereby facilitating incorporation of the high resolution MS datasets to future mechanistic models of SOM degradation and prediction of greenhouse gas emissions.

A pan-Arctic synthesis of CH4 and CO2 production from anoxic soil incubations.

Permafrost thaw can alter the soil environment through changes in soil moisture, frequently resulting in soil saturation, a shift to anaerobic decomposition, and changes in the plant community. These changes, along with thawing of previously frozen organic material, can alter the form and magnitude of greenhouse gas production from permafrost ecosystems. We synthesized existing methane (CH4 ) and carbon dioxide (CO2 ) production measurements from anaerobic incubations of boreal and tundra soils from the geographic permafrost region to evaluate large-scale controls of anaerobic CO2 and CH4 production and compare the relative importance of landscape-level factors (e.g., vegetation type and landscape position), soil properties (e.g., pH, depth, and soil type), and soil environmental conditions (e.g., temperature and relative water table position). We found fivefold higher maximum CH4 production per gram soil carbon from organic soils than mineral soils. Maximum CH4 production from soils in the active layer (ground that thaws and refreezes annually) was nearly four times that of permafrost per gram soil carbon, and CH4 production per gram soil carbon was two times greater from sites without permafrost than sites with permafrost. Maximum CH4 and median anaerobic CO2 production decreased with depth, while CO2 :CH4 production increased with depth. Maximum CH4 production was highest in soils with herbaceous vegetation and soils that were either consistently or periodically inundated. This synthesis identifies the need to consider biome, landscape position, and vascular/moss vegetation types when modeling CH4 production in permafrost ecosystems and suggests the need for longer-term anaerobic incubations to fully capture CH4 dynamics. Our results demonstrate that as climate warms in arctic and boreal regions, rates of anaerobic CO2 and CH4 production will increase, not only as a result of increased temperature, but also from shifts in vegetation and increased ground saturation that will accompany permafrost thaw.

Stoichiometry and temperature sensitivity of methanogenesis and CO2 production from saturated polygonal tundra in Barrow, Alaska.

Arctic permafrost ecosystems store ~50% of global belowground carbon (C) that is vulnerable to increased microbial degradation with warmer active layer temperatures and thawing of the near surface permafrost. We used anoxic laboratory incubations to estimate anaerobic CO2 production and methanogenesis in active layer (organic and mineral soil horizons) and permafrost samples from center, ridge and trough positions of water-saturated low-centered polygon in Barrow Environmental Observatory, Barrow AK, USA. Methane (CH4 ) and CO2 production rates and concentrations were determined at -2, +4, or +8 °C for 60 day incubation period. Temporal dynamics of CO2 production and methanogenesis at -2 °C showed evidence of fundamentally different mechanisms of substrate limitation and inhibited microbial growth at soil water freezing points compared to warmer temperatures. Nonlinear regression better modeled the initial rates and estimates of Q10 values for CO2 that showed higher sensitivity in the organic-rich soils of polygon center and trough than the relatively drier ridge soils. Methanogenesis generally exhibited a lag phase in the mineral soils that was significantly longer at -2 °C in all horizons. Such discontinuity in CH4 production between -2 °C and the elevated temperatures (+4 and +8 °C) indicated the insufficient representation of methanogenesis on the basis of Q10 values estimated from both linear and nonlinear models. Production rates for both CH4 and CO2 were substantially higher in organic horizons (20% to 40% wt. C) at all temperatures relative to mineral horizons (<20% wt. C). Permafrost horizon (~12% wt. C) produced ~5-fold less CO2 than the active layer and negligible CH4 . High concentrations of initial exchangeable Fe(II) and increasing accumulation rates signified the role of iron as terminal electron acceptors for anaerobic C degradation in the mineral horizons.