Construction of Viable Soil Defined Media Using Quantitative Metabolomics Analysis of Soil Metabolites

Soil metabolomics: Deciphering underground metabolic webs in terrestrial ecosystems

Soil metabolomics is an emerging approach for profiling diverse small molecule metabolites, i.e., metabolomes, in the soil. Soil metabolites, including fatty acids, amino acids, lipids, organic acids, sugars, and volatile organic compounds, often contain essential nutrients such as nitrogen, phosphorus, and sulfur and are directly linked to soil biogeochemical cycles driven by soil microorganisms. This paper presents an overview of methods for analyzing soil metabolites and the state-of-the-art of soil metabolomics in relation to soil nutrient cycling. We describe important applications of metabolomics in studying soil carbon cycling and sequestration, and the response of soil organic pools to changing environmental conditions. This includes using metabolomics to provide new insights into the close relationships between soil microbiome and metabolome, as well as responses of soil metabolome to plant and environmental stresses such as soil contamination. We also highlight the advantage of using soil metabolomics to study the biogeochemical cycles of elements and suggest that future research needs to better understand factors driving soil function and health.

Developing stable, simplified, functional consortia from Brachypodium rhizosphere for microbial application in sustainable agriculture

The rhizosphere microbiome plays a crucial role in supporting plant productivity and ecosystem functioning by regulating nutrient cycling, soil integrity, and carbon storage. However, deciphering the intricate interplay between microbial relationships within the rhizosphere is challenging due to the overwhelming taxonomic and functional diversity. Here we present our systematic design framework built on microbial colocalization and microbial interaction, toward successful assembly of multiple rhizosphere-derived Reduced Complexity Consortia (RCC). We enriched co-localized microbes from Brachypodium roots grown in field soil with carbon substrates mimicking Brachypodium root exudates, generating 768 enrichments. By transferring the enrichments every 3 or 7 days for 10 generations, we developed both fast and slow-growing reduced complexity microbial communities. Most carbon substrates led to highly stable RCC just after a few transfers. 16S rRNA gene amplicon analysis revealed distinct community compositions based on inoculum and carbon source, with complex carbon enriching slow growing yet functionally important soil taxa like Acidobacteria and Verrucomicrobia. Network analysis showed that microbial consortia, whether differentiated by growth rate (fast vs. slow) or by succession (across generations), had significantly different network centralities. Besides, the keystone taxa identified within these networks belong to genera with plant growth-promoting traits, underscoring their critical function in shaping rhizospheric microbiome networks. Furthermore, tested consortia demonstrated high stability and reproducibility, assuring successful revival from glycerol stocks for long-term viability and use. Our study represents a significant step toward developing a framework for assembling rhizosphere consortia based on microbial colocalization and interaction, with future implications for sustainable agriculture and environmental management.

The core metabolome and root exudation dynamics of three phylogenetically distinct plant species

Root exudates are plant-derived, exported metabolites likely shaping root-associated microbiomes by acting as nutrients and signals. However, root exudation dynamics are unclear and thus also, if changes in exudation are reflected in changes in microbiome structure. Here, we assess commonalities and differences between exudates of different plant species, diurnal exudation dynamics, as well as the accompanying methodological aspects of exudate sampling. We find that exudates should be collected for hours rather than days as many metabolite abundances saturate over time. Plant growth in sterile, nonsterile, or sugar-supplemented environments significantly alters exudate profiles. A comparison of Arabidopsis thaliana, Brachypodium distachyon, and Medicago truncatula shoot, root, and root exudate metabolite profiles reveals clear differences between these species, but also a core metabolome for tissues and exudates. Exudate profiles also exhibit a diurnal signature. These findings add to the methodological and conceptual groundwork for future exudate studies to improve understanding of plant-microbe interactions.

Environmental pseudotargeted metabolomics: A high throughput and wide coverage method for metabolic profiling of 1000-year paddy soil chronosequences

Pseudotargeted metabolomics is achieved by introducing an algorithm designed to choose ions for selected ion monitoring from identified metabolites. This method integrates the advantages of both untargeted and targeted metabolomics. In this study, environmental pseudotargeted metabolomics was established to analyze the soil metabolites, based on microwave assisted derivatization followed by gas chromatography-mass spectrometry analysis. The method development included the optimization of extraction factors and derivatization conditions, evaluation of silylation reagent types and matrix-dependent behaviors. Under the optimal conditions, the microwave oximation and silylation were completed in 5 min and 9 min. A total of 184 metabolites from 26 chemical classifications were identified in soil matrices. The method validation demonstrated excellent performance in terms of linearity (correlation coefficient > 0.99), repeatability (relative standard deviation (RSD) < 20 %), reproducibility (RSD < 25 %), stability (relative difference < 10 % within 18 h), and sensitivity (16-110 times higher signal-to-noise ratio). This developed method was applied to characterize the metabolite compositions and metabolic profiling in a 1000-year paddy soil chronosequence. The relative abundance of trehalose was highest in 6-(40.3 %), 60-(55.8 %), 300-(67.7 %)and 1000-(61.7 %)years paddy soil, respectively, but long-chain fatty acids were most abundant in marine sediment (57.4 %). Forty-two characteristic metabolites were considered as primarily responsible for discriminating and characterizing the paddy soil chronosequences development and seven major metabolic pathways were altered. In addition, GC-MS metabolite profile presented better discriminating power in paddy soil ecosystem changes than phospholipid fatty acids (PLFAs). Overall, environmental pseudotargeted metabolomics can provide a high throughout and wide coverage approach for performing metabolic profiling in the soil research.

Balancing Trade-Offs Imposed by Growth Media and Mass Spectrometry for Bacterial Exometabolomics

The bacterial exometabolome consists of a vast array of specialized metabolites, many of which are only produced in response to specific environmental stimuli. For this reason, it is desirable to control the extracellular environment with a defined growth medium composed of pure ingredients. However, complex (undefined) media are expected to support the robust growth of a greater variety of microorganisms than defined media. Here, we investigate the trade-offs inherent to a range of complex and defined solid media for the growth of soil microorganisms, production of specialized metabolites, and detection of these compounds using direct infusion mass spectrometry. We find that complex media support growth of more soil microorganisms, as well as allowing for the detection of more previously discovered natural products as a fraction of total m/z features detected in each sample. However, the use of complex media often caused mass spectrometer injection failures and poor-quality mass spectra, which in some cases resulted in over a quarter of samples being removed from analysis. Defined media, while more limiting in growth, generated higher quality spectra and yielded more m/z features after background subtraction. These results inform future exometabolomic experiments requiring a medium that supports the robust growth of many soil microorganisms. IMPORTANCE Bacteria are capable of producing and secreting a rich diversity of specialized metabolites. Yet, much of their exometabolome remains hidden due to challenges associated with eliciting specialized metabolite production, labor-intensive sample preparation, and time-consuming analysis techniques. Using our versatile three-dimensional (3D)-printed culturing platform, SubTap, we demonstrate that rapid exometabolomic data collection from a diverse set of environmental bacteria is feasible. We optimized our platform by surveying Streptomyces isolated from soil on a variety of media types to assess viability, degree of specialized metabolite production, and compatibility with downstream LESA-DIMS analysis. Ultimately, this will enable data-rich experimentation, allowing for a better understanding of bacterial exometabolomes.

Root exudates and rhizosphere soil bacterial relationships of Nitraria tangutorum are linked to k-strategists bacterial community under salt stress

When plants are subjected to various biotic and abiotic stresses, the root system responds actively by secreting different types and amounts of bioactive compounds, while affects the structure of rhizosphere soil bacterial community. Therefore, understanding plant-soil-microbial interactions, especially the strength of microbial interactions, mediated by root exudates is essential. A short-term experiment was conducted under drought and salt stress to investigate the interaction between root exudates and Nitraria tangutorum rhizosphere bacterial communities. We found that drought and salt stress increased rhizosphere soil pH (9.32 and 20.6%) and electrical conductivity (1.38 and 11 times), respectively, while decreased organic matter (27.48 and 31.38%), total carbon (34.55 and 29.95%), and total phosphorus (20 and 28.57%) content of N. tangutorum rhizosphere soil. Organic acids, growth hormones, and sugars were the main differential metabolites of N. tangutorum under drought and salt stress. Salt stress further changed the N. tangutorum rhizosphere soil bacterial community structure, markedly decreasing the relative abundance of Bacteroidota as r-strategist while increasing that of Alphaproteobacteria as k-strategists. The co-occurrence network analysis showed that drought and salt stress reduced the connectivity and complexity of the rhizosphere bacterial network. Soil physicochemical properties and root exudates in combination with salt stress affect bacterial strategies and interactions. Our study revealed the mechanism of plant-soil-microbial interactions under the influence of root exudates and provided new insights into the responses of bacterial communities to stressful environments.

Artificial Soils Reveal Individual Factor Controls on Microbial Processes

Soil matrix properties influence microbial behaviors that underlie nutrient cycling, greenhouse gas production, and soil formation. However, the dynamic and heterogeneous nature of soils makes it challenging to untangle the effects of different matrix properties on microbial behaviors. To address this challenge, we developed a tunable artificial soil recipe and used these materials to study the abiotic mechanisms driving soil microbial growth and communication. When we used standardized matrices with varying textures to culture gas-reporting biosensors, we found that a Gram-negative bacterium (Escherichia coli) grew best in synthetic silt soils, remaining active over a wide range of soil matric potentials, while a Gram-positive bacterium (Bacillus subtilis) preferred sandy soils, sporulating at low water potentials. Soil texture, mineralogy, and alkalinity all attenuated the bioavailability of an acyl-homoserine lactone (AHL) signaling molecule that controls community-level microbial behaviors. Texture controlled the timing of AHL sensing, while AHL bioavailability was decreased ~10⁵-fold by mineralogy and ~10³-fold by alkalinity. Finally, we built artificial soils with a range of complexities that converge on the properties of one Mollisol. As artificial soil complexity increased to more closely resemble the Mollisol, microbial behaviors approached those occurring in the natural soil, with the notable exception of organic matter. IMPORTANCE Understanding environmental controls on soil microbes is difficult because many abiotic parameters vary simultaneously and uncontrollably when different natural soils are compared, preventing mechanistic determination of any individual soil parameter's effect on microbial behaviors. We describe how soil texture, mineralogy, pH, and organic matter content can be varied individually within artificial soils to study their effects on soil microbes. Using microbial biosensors that report by producing a rare indicator gas, we identify soil properties that control microbial growth and attenuate the bioavailability of a diffusible chemical used to control community-level behaviors. We find that artificial soils differentially affect signal bioavailability and the growth of Gram-negative (Escherichia coli) and Gram-positive (Bacillus subtilis) microbes. These artificial soils are useful for studying the mechanisms that underlie soil controls on microbial fitness, signaling, and gene transfer.

Substrate Utilization and Competitive Interactions Among Soil Bacteria Vary With Life-History Strategies

Microorganisms have evolved various life-history strategies to survive fluctuating resource conditions in soils. However, it remains elusive how the life-history strategies of microorganisms influence their processing of organic carbon, which may affect microbial interactions and carbon cycling in soils. Here, we characterized the genomic traits, exometabolite profiles, and interactions of soil bacteria representing copiotrophic and oligotrophic strategists. Isolates were selected based on differences in ribosomal RNA operon (rrn) copy number, as a proxy for life-history strategies, with pairs of “high” and “low” rrn copy number isolates represented within the Micrococcales, Corynebacteriales, and Bacillales. We found that high rrn isolates consumed a greater diversity and amount of substrates than low rrn isolates in a defined growth medium containing common soil metabolites. We estimated overlap in substrate utilization profiles to predict the potential for resource competition and found that high rrn isolates tended to have a greater potential for competitive interactions. The predicted interactions positively correlated with the measured interactions that were dominated by negative interactions as determined through sequential growth experiments. This suggests that resource competition was a major force governing interactions among isolates, while cross-feeding of metabolic secretion likely contributed to the relatively rare positive interactions observed. By connecting bacterial life-history strategies, genomic features, and metabolism, our study advances the understanding of the links between bacterial community composition and the transformation of carbon in soils.

A Defined Medium for Cultivation and Exometabolite Profiling of Soil Bacteria

Exometabolomics is an approach to assess how microorganisms alter, or react to their environments through the depletion and production of metabolites. It allows the examination of how soil microbes transform the small molecule metabolites within their environment, which can be used to study resource competition and cross-feeding. This approach is most powerful when used with defined media that enable tracking of all metabolites. However, microbial growth media have traditionally been developed for the isolation and growth of microorganisms but not metabolite utilization profiling through Liquid Chromatography Tandem Mass Spectrometry (LC-MS/MS). Here, we describe the construction of a defined medium, the Northen Lab Defined Medium (NLDM), that not only supports the growth of diverse soil bacteria but also is defined and therefore suited for exometabolomic experiments. Metabolites included in NLDM were selected based on their presence in R2A medium and soil, elemental stoichiometry requirements, as well as knowledge of metabolite usage by different bacteria. We found that NLDM supported the growth of 108 of the 110 phylogenetically diverse (spanning 36 different families) soil bacterial isolates tested and all of its metabolites were trackable through LC–MS/MS analysis. These results demonstrate the viability and utility of the constructed NLDM medium for growing and characterizing diverse microbial isolates and communities.

Nutrient Availability Shifts the Biosynthetic Potential of Soil-Derived Microbial Communities

Secondary metabolites produced by microorganisms are the main source of antimicrobials and other pharmaceutical drugs. Soil microbes have been the primary discovery source for these secondary metabolites, often producing complex organic compounds with specific biological activities. Research suggests that secondary metabolism broadly shapes microbial ecological interactions, but little is known about the factors that shape the abundance, distribution, and diversity of biosynthetic gene clusters in the context of microbial communities. In this study, we investigate the role of nutrient availability on the abundance of biosynthetic gene clusters in soil-derived microbial consortia. Soil microbial consortia enriched in high sugar medium (150 mg/L of glucose and 200 mg/L of trehalose) had more biosynthetic gene clusters and higher inhibitory activity than those enriched in low sugar medium (15 mg/L of glucose + 20 mg/L of trehalose). Our results demonstrate that experimental microbial communities are a promising tool to study the ecology of specialized metabolites.

Biofilm Interaction Mapping and Analysis (BIMA) of Interspecific Interactions in Pseudomonas Co-culture Biofilms

Pseudomonas species are ubiquitous in nature and include numerous medically, agriculturally and technologically beneficial strains of which the interspecific interactions are of great interest for biotechnologies. Specifically, co-cultures containing Pseudomonas stutzeri have been used for bioremediation, biocontrol, aquaculture management and wastewater denitrification. Furthermore, the use of P. stutzeri biofilms, in combination with consortia-based approaches, may offer advantages for these processes. Understanding the interspecific interaction within biofilm co-cultures or consortia provides a means for improvement of current technologies. However, the investigation of biofilm-based consortia has been limited. We present an adaptable and scalable method for the analysis of macroscopic interactions (colony morphology, inhibition, and invasion) between colony-forming bacterial strains using an automated printing method followed by analysis of the genes and metabolites involved in the interactions. Using Biofilm Interaction Mapping and Analysis (BIMA), these interactions were investigated between P. stutzeri strain RCH2, a denitrifier isolated from chromium (VI) contaminated soil, and 13 other species of pseudomonas isolated from non-contaminated soil. One interaction partner, Pseudomonas fluorescens N1B4 was selected for mutant fitness profiling of a DNA-barcoded mutant library; with this approach four genes of importance were identified and the effects on interactions were evaluated with deletion mutants and mass spectrometry based metabolomics.

Nucleoside Metabolism Is Induced in Common Bean During Early Seedling Development

Nucleoside hydrolases (NSH; nucleosidases) catalyze the cleavage of nucleosides into ribose and free nucleobases. These enzymes have been postulated as key elements controlling the ratio between nucleotide salvage and degradation. Moreover, they play a pivotal role in ureidic legumes by providing the substrate for the synthesis of ureides. Furthermore, nucleotide metabolism has a crucial role during germination and early seedling development, since the developing seedlings require high amount of nucleotide simultaneously to the mobilization of nutrient in cotyledons. In this study, we have cloned two nucleosidases genes from Phaseolus vulgaris, PvNSH1 and PvNSH2, expressed them as recombinant proteins, and characterized their catalytic activities. Both enzymes showed a broad range of substrate affinity; however, PvNSH1 exhibited the highest activity with uridine, followed by xanthosine, whereas PvNSH2 hydrolyses preferentially xanthosine and shows low activity with uridine. The study of the regulation of nucleosidases during germination and early postgerminative development indicated that nucleosidases are induced in cotyledons and embryonic axes just after the radicle emergence, coincident with the induction of nucleases activity and the synthesis of ureides in the embryonic axes, with no remarkable differences in the level of expression of both nucleosidase genes. In addition, nucleosides and nucleobase levels were determined as well in cotyledons and embryonic axes. Our results suggest that PvNSH1 and PvNSH2 play an important role in the mobilization of nutrients during this crucial stage of plant development.

Ecophysiological Study of Paraburkholderia sp. Strain 1N under Soil Solution Conditions: Dynamic Substrate Preferences and Characterization of Carbon Use Efficiency

We used time-resolved metabolic footprinting, an important technical approach used to monitor changes in extracellular compound concentrations during microbial growth, to study the order of substrate utilization (i.e., substrate preferences) and kinetics of a fast-growing soil isolate, Paraburkholderia sp. strain 1N. The growth of Paraburkholderia sp. 1N was monitored under aerobic conditions in a soil-extracted solubilized organic matter medium, representing a realistic diversity of available substrates and gradient of initial concentrations. We combined multiple analytical approaches to track over 150 compounds in the medium and complemented this with bulk carbon and nitrogen measurements, allowing estimates of carbon use efficiency throughout the growth curve. Targeted methods allowed the quantification of common low-molecular-weight substrates: glucose, 20 amino acids, and 9 organic acids. All targeted compounds were depleted from the medium, and depletion followed a sigmoidal curve where sufficient data were available. Substrates were utilized in at least three distinct temporal clusters as Paraburkholderia sp. 1N produced biomass at a cumulative carbon use efficiency of 0.43. The two substrates with highest initial concentrations, glucose and valine, exhibited longer usage windows, at higher biomass-normalized rates, and later in the growth curve. Contrary to hypotheses based on previous studies, we found no clear relationship between substrate nominal oxidation state of carbon (NOSC) or maximal growth rate and the order of substrate depletion. Under soil solution conditions, the growth of Paraburkholderia sp. 1N induced multiauxic substrate depletion patterns that could not be explained by the traditional paradigm of catabolite repression.IMPORTANCE Exometabolomic footprinting methods have the capability to provide time-resolved observations of the uptake and release of hundreds of compounds during microbial growth. Of particular interest is microbial phenotyping under environmentally relevant soil conditions, consisting of relatively low concentrations and modeling pulse input events. Here, we show that growth of a bacterial soil isolate, Paraburkholderia sp. 1N, on a dilute soil extract resulted in a multiauxic metabolic response, characterized by discrete temporal clusters of substrate depletion and metabolite production. Our data did not support the hypothesis that compounds with lower energy content are used preferentially, as each cluster contained compounds with a range of nominal oxidation states of carbon. These new findings with Paraburkholderia sp. 1N, which belongs to a metabolically diverse genus, provide insights on ecological strategies employed by aerobic heterotrophs competing for low-molecular-weight substrates in soil solution.

Root morphology and exudate availability are shaped by particle size and chemistry in Brachypodium distachyon

Root morphology and exudation define a plants' sphere of influence in soils. In turn, soil characteristics influence plant growth, morphology, root microbiome, and rhizosphere chemistry. Collectively, all these parameters have significant implications on the major biogeochemical cycles, crop yield, and ecosystem health. However, how plants are shaped by the physiochemistry of soil particles is still not well understood. We explored how particle size and chemistry of growth substrates affect root morphology and exudation of a model grass. We grew Brachypodium distachyon in glass beads with various sizes (0.5, 1, 2, 3 mm), as well as in sand (0.005, 0.25, 4 mm) and in clay (4 mm) particles and in particle-free hydroponic medium. Plant morphology, root weight, and shoot weight were measured. We found that particle size significantly influenced root fresh weight and root length, whereas root number and shoot weight remained constant. Next, plant exudation profiles were analyzed with mass spectrometry imaging and liquid chromatography-mass spectrometry. Mass spectrometry imaging suggested that both, root length and number shape root exudation. Exudate profiles were comparable for plants growing in glass beads or sand with various particles sizes, but distinct for plants growing in clay for in situ exudate collection. Clay particles were found to sorb 20% of compounds exuded by clay-grown plants, and 70% of compounds from a defined exudate medium. The sorbed compounds belonged to a range of chemical classes, among them nucleosides, organic acids, sugars, and amino acids. Some of the sorbed compounds could be desorbed by a rhizobacterium (Pseudomonas fluorescens WCS415), supporting its growth. This study demonstrates the effect of different characteristics of particles on root morphology, plant exudation and availability of nutrients to microorganisms. These findings further support the critical importance of the physiochemical properties of soils when investigating plant morphology, plant chemistry, and plant-microbe interactions.

Local Phenomena Shape Backyard Soil Metabolite Composition

Soil covers most of Earth's continental surface and is fundamental to life-sustaining processes such as agriculture. Given its rich biodiversity, soil is also a major source for natural product drug discovery from soil microorganisms. However, the study of the soil small molecule profile has been challenging due to the complexity and heterogeneity of this matrix. In this study, we implemented high-resolution liquid chromatography-tandem mass spectrometry and large-scale data analysis tools such as molecular networking to characterize the relative contributions of city, state and regional processes on backyard soil metabolite composition, in 188 soil samples collected from 14 USA States, representing five USA climate regions. We observed that region, state and city of collection all influence the overall soil metabolite profile. However, many metabolites were only detected in unique sites, indicating that uniquely local phenomena also influence the backyard soil environment, with both human-derived and naturally-produced (plant-derived, microbially-derived) metabolites identified. Overall, these findings are helping to define the processes that shape the backyard soil metabolite composition, while also highlighting the need for expanded metabolomic studies of this complex environment.

Evaluation of an untargeted nano-liquid chromatography-mass spectrometry approach to expand coverage of low molecular weight dissolved organic matter in Arctic soil

Characterizing low molecular weight (LMW) dissolved organic matter (DOM) in soils and evaluating the availability of this labile pool is critical to understanding the underlying mechanisms that control carbon storage or release across terrestrial systems. However, due to wide-ranging physicochemical diversity, characterizing this complex mixture of small molecules and how it varies across space remains an analytical challenge. Here, we evaluate an untargeted approach to detect qualitative and relative-quantitative variations in LMW DOM with depth using water extracts from a soil core from the Alaskan Arctic, a unique system that contains nearly half the Earth's terrestrial carbon and is rapidly warming due to climate change. We combined reversed-phase and hydrophilic interaction liquid chromatography, and nano-electrospray ionization coupled with high-resolution tandem mass spectrometry in positive- and negative-ionization mode. The optimized conditions were sensitive, robust, highly complementary, and enabled detection and putative annotations of a wide range of compounds (e.g. amino acids, plant/microbial metabolites, sugars, lipids, peptides). Furthermore, multivariate statistical analyses revealed subtle but consistent and significant variations with depth. Thus, this platform is useful not only for characterizing LMW DOM, but also for quantifying relative variations in LMW DOM availability across space, revealing hotspots of biogeochemical activity for further evaluation.