Genomic fingerprints of the world's soil ecosystems.

Despite the explosion of soil metagenomic data, we lack a synthesized understanding of patterns in the distribution and functions of soil microorganisms. These patterns are critical to predictions of soil microbiome responses to climate change and resulting feedbacks that regulate greenhouse gas release from soils. To address this gap, we assay 1,512 manually curated soil metagenomes using complementary annotation databases, read-based taxonomy, and machine learning to extract multidimensional genomic fingerprints of global soil microbiomes. Our objective is to uncover novel biogeographical patterns of soil microbiomes across environmental factors and ecological biomes with high molecular resolution. We reveal shifts in the potential for (i) microbial nutrient acquisition across pH gradients; (ii) stress-, transport-, and redox-based processes across changes in soil bulk density; and (iii) greenhouse gas emissions across biomes. We also use an unsupervised approach to reveal a collection of soils with distinct genomic signatures, characterized by coordinated changes in soil organic carbon, nitrogen, and cation exchange capacity and in bulk density and clay content that may ultimately reflect soil environments with high microbial activity. Genomic fingerprints for these soils highlight the importance of resource scavenging, plant-microbe interactions, fungi, and heterotrophic metabolisms. Across all analyses, we observed phylogenetic coherence in soil microbiomes-more closely related microorganisms tended to move congruently in response to soil factors. Collectively, the genomic fingerprints uncovered here present a basis for global patterns in the microbial mechanisms underlying soil biogeochemistry and help beget tractable microbial reaction networks for incorporation into process-based models of soil carbon and nutrient cycling.IMPORTANCEWe address a critical gap in our understanding of soil microorganisms and their functions, which have a profound impact on our environment. We analyzed 1,512 global soils with advanced analytics to create detailed genetic profiles (fingerprints) of soil microbiomes. Our work reveals novel patterns in how microorganisms are distributed across different soil environments. For instance, we discovered shifts in microbial potential to acquire nutrients in relation to soil acidity, as well as changes in stress responses and potential greenhouse gas emissions linked to soil structure. We also identified soils with putative high activity that had unique genomic characteristics surrounding resource acquisition, plant-microbe interactions, and fungal activity. Finally, we observed that closely related microorganisms tend to respond in similar ways to changes in their surroundings. Our work is a significant step toward comprehending the intricate world of soil microorganisms and its role in the global climate.

Rethinking Aerobic Respiration in the Hyporheic Zone under Variation in Carbon and Nitrogen Stoichiometry.

Hyporheic zones (HZs)─zones of groundwater-surface water mixing─are hotspots for dissolved organic matter (DOM) and nutrient cycling that can disproportionately impact aquatic ecosystem functions. However, the mechanisms affecting DOM metabolism through space and time in HZs remain poorly understood. To resolve this gap, we investigate a recently proposed theory describing trade-offs between carbon (C) and nitrogen (N) limitations as a key regulator of HZ metabolism. We propose that throughout the extent of the HZ, a single process like aerobic respiration (AR) can be limited by both DOM thermodynamics and N content due to highly variable C/N ratios over short distances (centimeter scale). To investigate this theory, we used a large flume, continuous optode measurements of dissolved oxygen (DO), and spatially and temporally resolved molecular analysis of DOM. Carbon and N limitations were inferred from changes in the elemental stoichiometric ratio. We show sequential, depth-stratified relationships of DO with DOM thermodynamics and organic N that change across centimeter scales. In the shallow HZ with low C/N, DO was associated with the thermodynamics of DOM, while deeper in the HZ with higher C/N, DO was associated with inferred biochemical reactions involving organic N. Collectively, our results suggest that there are multiple competing processes that limit AR in the HZ. Resolving this spatiotemporal variation could improve predictions from mechanistic models, either via more highly resolved grid cells or by representing AR colimitation by DOM thermodynamics and organic N.

Disinfection byproducts formed during drinking water treatment reveal an export control point for dissolved organic matter in a subalpine headwater stream.

Changes in climate, season, and vegetation can alter organic export from watersheds. While an accepted tradeoff to protect public health, disinfection processes during drinking water treatment can adversely react with organic compounds to form disinfection byproducts (DBPs). By extension, DBP monitoring can yield insights into hydrobiogeochemical dynamics within watersheds and their implications for water resource management. In this study, we analyzed temporal trends from a water treatment facility that sources water from Coal Creek in Crested Butte, Colorado. These trends revealed a long-term increase in haloacetic acid and trihalomethane formation over the period of 2005-2020. Disproportionate export of dissolved organic carbon and formation of DBPs that exceeded maximum contaminant levels were consistently recorded in association with late spring freshet. Synoptic sampling of the creek in 2020 and 2021 identified a biogeochemical hotspot for organic carbon export in the upper domain of the watershed that contained a prominent fulvic acid-like fluorescent signature. DBP formation potential analyses from this domain yielded similar ratios of regulated DBP classes to those formed at the drinking water facility. Spectrometric qualitative analyses of pre and post-reacted waters with hypochlorite indicated lignin-like and condensed hydrocarbon-like molecules were the major reactive chemical classes during chlorine-based disinfection. This study demonstrates how drinking water quality archives combined with synoptic sampling and targeted analyses can be used to identify and understand export control points for dissolved organic matter. This approach could be applied to identify and characterize analogous watersheds where seasonal or climate-associated organic matter export challenge water treatment disinfection and by extension inform watershed management and drinking water treatment.

Implications of sample treatment on characterization of riverine dissolved organic matter.

High-resolution mass spectrometry techniques are widely used in the environmental sciences to characterize natural organic matter and, when utilizing these instruments, researchers must make multiple decisions regarding sample pre-treatment and the instrument ionization mode. To identify how these choices alter organic matter characterization and resulting conclusions, we analyzed a collection of 17 riverine samples from East River, CO (USA) under four PPL-based Solid Phase Extraction (SPE) treatment and electrospray ionization polarity (e.g., positive and negative) combinations: SPE (+), SPE (-), non-SPE (-), and non-SPE (+). The greatest number of formula assignments were achieved with SPE-treated samples due to the removal of compounds that could interfere with ionization. Furthermore, the SPE (-) treatment captured the most formulas across the widest chemical compound diversity. In addition to a reduced number of assigned formulas, the non-SPE datasets resulted in altered thermodynamic interpretations that could cascade into incomplete assumptions about the availability of organic matter pools for heterotrophic microbial respiration. Thus, we infer that the SPE (-) treatment is the best single method for characterizing environmental organic matter pools unless the focus is on lipid-like compounds, in which case we recommend a combination of SPE (-) and SPE (+) to adequately characterize these molecules.

Inferring the Contribution of Microbial Taxa and Organic Matter Molecular Formulas to Ecological Assembly.

Understanding the mechanisms underlying the assembly of communities has long been the goal of many ecological studies. While several studies have evaluated community wide ecological assembly, fewer have focused on investigating the impacts of individual members within a community or assemblage on ecological assembly. Here, we adapted a previous null model ß-nearest taxon index (ßNTI) to measure the contribution of individual features within an ecological community to overall assembly. This new metric, called feature-level ßNTI (ßNTIfeat), enables researchers to determine whether ecological features (e.g., individual microbial taxa) contribute to divergence, convergence, or have insignificant impacts across spatiotemporally resolved metacommunities or meta-assemblages. Using ßNTIfeat, we revealed that unclassified microbial lineages often contributed to community divergence while diverse groups (e.g., Crenarchaeota, Alphaproteobacteria, and Gammaproteobacteria) contributed to convergence. We also demonstrate that ßNTIfeat can be extended to other ecological assemblages such as organic molecules comprising organic matter (OM) pools. OM had more inconsistent trends compared to the microbial community though CHO-containing molecular formulas often contributed to convergence, while nitrogen and phosphorus-containing formulas contributed to both convergence and divergence. A network analysis was used to relate ßNTIfeat values from the putatively active microbial community and the OM assemblage and examine potentially common contributions to ecological assembly across different communities/assemblages. This analysis revealed that P-containing formulas often contributed to convergence/divergence separately from other ecological features and N-containing formulas often contributed to assembly in coordination with microorganisms. Additionally, members of Family Geobacteraceae were often observed to contribute to convergence/divergence in conjunction with both N- and P-containing formulas, suggesting a coordinated ecological role for family members and the nitrogen/phosphorus cycle. Overall, we show that ßNTIfeat offers opportunities to investigate the community or assemblage members, which shape the phylogenetic or functional landscape, and demonstrate the potential to evaluate potential points of coordination across various community types.

Ecological theory applied to environmental metabolomes reveals compositional divergence despite conserved molecular properties.

Stream and river systems transport and process substantial amounts of dissolved organic matter (DOM) from terrestrial and aquatic sources to the ocean, with global biogeochemical implications. However, the underlying mechanisms affecting the spatiotemporal organization of DOM composition are under-investigated. To understand the principles governing DOM composition, we leverage the recently proposed synthesis of metacommunity ecology and metabolomics, termed 'meta-metabolome ecology.' Applying this novel approach to a freshwater ecosystem, we demonstrated that despite similar molecular properties across metabolomes, metabolite identity significantly diverged due to environmental filtering and variations in putative biochemical transformations. We refer to this phenomenon as 'thermodynamic redundancy,' which is analogous to the ecological concept of functional redundancy. We suggest that under thermodynamic redundancy, divergent metabolomes can support equivalent biogeochemical function just as divergent ecological communities can support equivalent ecosystem function. As these analyses are performed in additional ecosystems, potentially generalizable concepts, like thermodynamic redundancy, can be revealed and provide insight into DOM dynamics.

Representing Organic Matter Thermodynamics in Biogeochemical Reactions via Substrate-Explicit Modeling.

Predictive biogeochemical modeling requires data-model integration that enables explicit representation of the sophisticated roles of microbial processes that transform substrates. Data from high-resolution organic matter (OM) characterization are increasingly available and can serve as a critical resource for this purpose, but their incorporation into biogeochemical models is often prohibited due to an over-simplified description of reaction networks. To fill this gap, we proposed a new concept of biogeochemical modeling-termed substrate-explicit modeling-that enables parameterizing OM-specific oxidative degradation pathways and reaction rates based on the thermodynamic properties of OM pools. Based on previous developments in the literature, we characterized the resulting kinetic models by only two parameters regardless of the complexity of OM profiles, which can greatly facilitate the integration with reactive transport models for ecosystem simulations by alleviating the difficulty in parameter identification. The two parameters include maximal growth rate (µmax) and harvest volume (Vh) (i.e., the volume that a microbe can access for harvesting energy). For every detected organic molecule in a given sample, our approach provides a systematic way to formulate reaction kinetics from chemical formula, which enables the evaluation of the impact of OM character on biogeochemical processes across conditions. In a case study of two sites with distinct OM thermodynamics using ultra high-resolution metabolomics datasets derived from Fourier transform ion cyclotron resonance mass spectrometry analyses, our method predicted how oxidative degradation is primarily driven by thermodynamic efficiency of OM consistent with experimental rate measurements (as shown by correlation coefficients of up to 0.61), and how biogeochemical reactions can vary in response to carbon and/or oxygen limitations. Lastly, we showed that incorporation of enzymatic regulations into substrate-explicit models is critical for more reasonable predictions. This result led us to present integrative biogeochemical modeling as a unifying framework that can ideally describe the dynamic interplay among microbes, enzymes, and substrates to address advanced questions and hypotheses in future studies. Altogether, the new modeling concept we propose in this work provides a foundational platform for unprecedented predictions of biogeochemical and ecosystem dynamics through enhanced integration with diverse experimental data and extant modeling approaches.

Using Community Science to Reveal the Global Chemogeography of River Metabolomes.

River corridor metabolomes reflect organic matter (OM) processing that drives aquatic biogeochemical cycles. Recent work highlights the power of ultrahigh-resolution mass spectrometry for understanding metabolome composition and river corridor metabolism. However, there have been no studies on the global chemogeography of surface water and sediment metabolomes using ultrahigh-resolution techniques. Here, we describe a community science effort from the Worldwide Hydrobiogeochemistry Observation Network for Dynamic River Systems (WHONDRS) consortium to characterize global metabolomes in surface water and sediment that span multiple stream orders and biomes. We describe the distribution of key aspects of metabolomes including elemental groups, chemical classes, indices, and inferred biochemical transformations. We show that metabolomes significantly differ across surface water and sediment and that surface water metabolomes are more rich and variable. We also use inferred biochemical transformations to identify core metabolic processes shared among surface water and sediment. Finally, we observe significant spatial variation in sediment metabolites between rivers in the eastern and western portions of the contiguous United States. Our work not only provides a basis for understanding global patterns in river corridor biogeochemical cycles but also demonstrates that community science endeavors can enable global research projects that are unfeasible with traditional research models.