Tracing Carbon Metabolism with Stable Isotope Metabolomics Reveals the Legacy of Diverse Carbon Sources in Soil.

Tracking the metabolic activity of whole soil communities can improve our understanding of the transformation and fate of carbon in soils. We used stable isotope metabolomics to trace 13C from nine labeled carbon sources into the water-soluble metabolite pool of an agricultural soil over time. Soil was amended with a mixture of all nine sources, with one source isotopically labeled in each treatment. We compared changes in the 13C enrichment of metabolites with respect to carbon source and time over a 48-day incubation and contrasted differences between soluble sources (glucose, xylose, amino acids, etc.) and insoluble sources (cellulose and palmitic acid). Whole soil metabolite profiles varied singularly by time, while the composition of 13C-labeled metabolites differed primarily by carbon source (R2 = 0.68) rather than time (R2 = 0.07), with source-specific differences persisting throughout incubations. The 13C labeling of metabolites from insoluble carbon sources occurred slower than that from soluble sources but yielded a higher average atom percent (atom%) 13C in metabolite markers of biomass (amino acids and nucleic acids). The 13C enrichment of metabolite markers of biomass stabilized between 5 and 15 atom% 13C by the end of incubations. Temporal patterns in the 13C enrichment of tricarboxylic acid cycle intermediates, nucleobases (uracil and thymine), and by-products of DNA salvage (allantoin) closely tracked microbial activity. Our results demonstrate that metabolite production in soils is driven by the carbon source supplied to the community and that the fate of carbon in metabolites do not generally converge over time as a result of ongoing microbial processing and recycling. IMPORTANCE Carbon metabolism in soil remains poorly described due to the inherent difficulty of obtaining information on the microbial metabolites produced by complex soil communities. Our study demonstrates the use of stable isotope probing (SIP) to study carbon metabolism in soil by tracking 13C from supplied carbon sources into metabolite pools and biomass. We show that differences in the metabolism of sources influence the fate of carbon in soils. Heterogeneity in 13C-labeled metabolite profiles corresponded with compositional differences in the metabolically active populations, providing a basis for how microbial community composition correlates with the quality of soil carbon. Our study demonstrates the application of SIP-metabolomics in studying soils and identifies several metabolite markers of growth, activity, and other aspects of microbial function.

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.

Long-read metagenomics of soil communities reveals phylum-specific secondary metabolite dynamics.

Microbial biosynthetic gene clusters (BGCs) encoding secondary metabolites are thought to impact a plethora of biologically mediated environmental processes, yet their discovery and functional characterization in natural microbiomes remains challenging. Here we describe deep long-read sequencing and assembly of metagenomes from biological soil crusts, a group of soil communities that are rich in BGCs. Taking advantage of the unusually long assemblies produced by this approach, we recovered nearly 3,000 BGCs for analysis, including 712 full-length BGCs. Functional exploration through metatranscriptome analysis of a 3-day wetting experiment uncovered phylum-specific BGC expression upon activation from dormancy, elucidating distinct roles and complex phylogenetic and temporal dynamics in wetting processes. For example, a pronounced increase in BGC transcription occurs at night primarily in cyanobacteria, implicating BGCs in nutrient scavenging roles and niche competition. Taken together, our results demonstrate that long-read metagenomic sequencing combined with metatranscriptomic analysis provides a direct view into the functional dynamics of BGCs in environmental processes and suggests a central role of secondary metabolites in maintaining phylogenetically conserved niches within biocrusts.

Characteristics of Wetting-Induced Bacteriophage Blooms in Biological Soil Crust.

Biological soil crusts (biocrusts) are photosynthetic "hot spots" in deserts and cover ∼12% of the Earth's terrestrial surface, and yet they face an uncertain future given expected shifts in rainfall events. Laboratory wetting of biocrust communities is known to cause a bloom of Firmicutes which rapidly become dominant community members within 2 days after emerging from a sporulated state. We hypothesized that their bacteriophages (phages) would respond to such a dramatic increase in their host's abundance. In our experiment, wetting caused Firmicutes to bloom and triggered a significant depletion of cyanobacterial diversity. We used genome-resolved metagenomics to link phage to their hosts and found that the bloom of the genus Bacillus correlated with a dramatic increase in the number of Caudovirales phages targeting these diverse spore-formers (r = 0.762). After 2 days, we observed dramatic reductions in the relative abundances of Bacillus, while the number of Bacillus phages continued to increase, suggestive of a predator-prey relationship. We found predicted auxiliary metabolic genes (AMGs) associated with sporulation in several Caudovirales genomes, suggesting that phages may influence and even benefit from sporulation dynamics in biocrusts. Prophage elements and CRISPR-Cas repeats in Firmicutes metagenome-assembled genomes (MAGs) provide evidence of recent infection events by phages, which were corroborated by mapping viral contigs to their host MAGs. Combined, these findings suggest that the blooming Firmicutes become primary targets for biocrust Caudovirales phages, consistent with the classical "kill-the-winner" hypothesis.IMPORTANCE This work forms part of an overarching research theme studying the effects of a changing climate on biological soil crust (biocrust) in the Southwestern United States. To our knowledge, this study was the first to characterize bacteriophages in biocrust and offers a view into the ecology of phages in response to a laboratory wetting experiment. The phages identified here represent lineages of Caudovirales, and we found that the dynamics of their interactions with their Firmicutes hosts explain the collapse of a bacterial bloom that was induced by wetting. Moreover, we show that phages carried host-altering metabolic genes and found evidence of proviral infection and CRISPR-Cas repeats within host genomes. Our results suggest that phages exert controls on population density by lysing dominant bacterial hosts and that they further impact biocrust by acquiring host genes for sporulation. Future research should explore how dominant these phages are in other biocrust communities and quantify how much the control and lysis of blooming populations contributes to nutrient cycling in biocrusts.

Learning representations of microbe-metabolite interactions.

Integrating multiomics datasets is critical for microbiome research; however, inferring interactions across omics datasets has multiple statistical challenges. We solve this problem by using neural networks (https://github.com/biocore/mmvec) to estimate the conditional probability that each molecule is present given the presence of a specific microorganism. We show with known environmental (desert soil biocrust wetting) and clinical (cystic fibrosis lung) examples, our ability to recover microbe-metabolite relationships, and demonstrate how the method can discover relationships between microbially produced metabolites and inflammatory bowel disease.

Untargeted Soil Metabolomics Using Liquid Chromatography-Mass Spectrometry and Gas Chromatography-Mass Spectrometry.

The molecular composition of soil organic matter (SOM) sets the foundation for terrestrial microbial community structures and carbon cycling dynamics. However, the specific chemical constituents of SOM are underexplored. In this chapter we present a protocol for the extraction of small molecule metabolites from soil followed by compound detection and identification using liquid chromatography-mass spectrometry and gas chromatography-mass spectrometry. There are options within the protocol to assess either the extracellular pool of metabolites or the total pool (including intracellular) and either polar or nonpolar metabolites, depending on the reader's research interests. These methods can be followed individually for a more targeted analysis or all methods can be combined to obtain a more comprehensive understanding of SOM metabolite composition (such as amino acids, nucleobases, organic acids, fatty acids, carbohydrates, secondary metabolites, and antibiotics).

Linking soil biology and chemistry in biological soil crust using isolate exometabolomics.

Metagenomic sequencing provides a window into microbial community structure and metabolic potential; however, linking these data to exogenous metabolites that microorganisms process and produce (the exometabolome) remains challenging. Previously, we observed strong exometabolite niche partitioning among bacterial isolates from biological soil crust (biocrust). Here we examine native biocrust to determine if these patterns are reproduced in the environment. Overall, most soil metabolites display the expected relationship (positive or negative correlation) with four dominant bacteria following a wetting event and across biocrust developmental stages. For metabolites that were previously found to be consumed by an isolate, 70% are negatively correlated with the abundance of the isolate's closest matching environmental relative in situ, whereas for released metabolites, 67% were positively correlated. Our results demonstrate that metabolite profiling, shotgun sequencing and exometabolomics may be successfully integrated to functionally link microbial community structure with environmental chemistry in biocrust.

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

Exometabolomics enables analysis of metabolite utilization of low molecular weight organic substances by soil bacteria. Environmentally-based defined media are needed to examine ecologically relevant patterns of substrate utilization. Here, we describe an approach for the construction of defined media using untargeted characterization of water soluble soil microbial metabolites from a saprolite soil collected from the Oak Ridge Field Research Center (ORFRC). To broadly characterize metabolites, both liquid chromatography mass spectrometry (LC/MS) and gas chromatography mass spectrometry (GC/MS) were used. With this approach, 96 metabolites were identified, including amino acids, amino acid derivatives, sugars, sugar alcohols, mono- and di-carboxylic acids, nucleobases, and nucleosides. From this pool of metabolites, 25 were quantified. Molecular weight cut-off filtration determined the fraction of carbon accounted for by the quantified metabolites and revealed that these soil metabolites have an uneven quantitative distribution (e.g., trehalose accounted for 9.9% of the <1 kDa fraction). This quantitative information was used to formulate two soil defined media (SDM), one containing 23 metabolites (SDM1) and one containing 46 (SDM2). To evaluate the viability of the SDM, we examined the growth of 30 phylogenetically diverse soil bacterial isolates from the ORFRC field site. The simpler SDM1 supported the growth of 13 isolates while the more complex SDM2 supported 15 isolates. To investigate SDM1 substrate preferences, one isolate, Pseudomonas corrugata strain FW300-N2E2 was selected for a time-series exometabolomics analysis. Interestingly, it was found that this organism preferred lower-abundance substrates such as guanine, glycine, proline and arginine and glucose and did not utilize the more abundant substrates maltose, mannitol, trehalose and uridine. These results demonstrate the viability and utility of using exometabolomics to construct a tractable environmentally relevant media. We anticipate that this approach can be expanded to other environments to enhance isolation and characterization of diverse microbial communities.

Cycloxaprid insecticide: nicotinic acetylcholine receptor binding site and metabolism.

Cycloxaprid (CYC) is a novel neonicotinoid prepared from the (nitromethylene)imidazole (NMI) analogue of imidacloprid. In this study we consider whether CYC is active per se or only as a proinsecticide for NMI. The IC50 values (nM) for displacing [(3)H]NMI binding are 43-49 for CYC and 2.3-3.2 for NMI in house fly and honeybee head membranes and 302 and 7.2, respectively, in mouse brain membranes, potency relationships interpreted as partial conversion of some CYC to NMI under the assay conditions. The 6-8-fold difference in toxicity of injected CYC and NMI to house flies is consistent with their relative potencies as in vivo nicotinic acetylcholine receptor (nAChR) inhibitors in brain measured with [(3)H]NMI binding assays. CYC metabolism in mice largely involves cytochrome P450 pathways without NMI as a major intermediate. Metabolites of CYC tentatively assigned are five monohydroxy derivatives and one each of dihydroxy, nitroso, and amino modifications. CYC appears be a proinsecticide, serving as a slow-release reservoir for NMI with selective activity for insect versus mammalian nAChRs.

Aldehyde oxidase importance in vivo in xenobiotic metabolism: imidacloprid nitroreduction in mice.

Aldehyde oxidase (AOX) metabolizes many xenobiotics in vitro, but its importance in vivo is usually unknown relative to cytochrome P450s (CYPs) and other detoxification systems. Currently, the most important insecticides are neonicotinoids, which are metabolized in vitro by AOX on reduction of the nitroimino group and by CYPs via oxidation reactions. The goal of this study was to establish the relative importance of AOX and CYPs in vivo using the mouse model. The procedure was to reduce liver AOX activity by providing tungsten or hydralazine in the drinking water or to use the AOX-deficient DBA/2 mouse strain. None of these approaches reduced CYP activity measured in vitro with an isozyme nonspecific substrate. Liver AOX activity was reduced by 45% with tungsten and 61% with hydralazine and 81% in AOX-deficient mice relative to controls. When mice were treated ip with the major neonicotinoid imidacloprid (IMI), metabolism by CYP oxidation reactions was not appreciably affected, whereas the AOX-generated nitrosoguanidine metabolite was decreased by 30% with tungsten and 56% with hydralazine and 86% in the AOX-deficient mice. The other IMI nitroreduction metabolite, desnitro-IMI, was decreased by 55%, 65%, and 81% with tungsten, hydralazine, and in the AOX-deficient mice, respectively. Thus, decreasing liver AOX activity by three quite different procedures gave a corresponding decrease for in vivo reductive metabolites in the liver of IMI-treated mice. Possible AOX involvement in IMI metabolism in insects was evaluated using AOX-expressing and AOX-deficient Drosophila, but no differences were found in IMI nitroreduction or sensitivity between the two strains. This is the first study to establish the in vivo relevance of AOX in neonicotinoid metabolism in mammals and one of the first for xenobiotics in general.

Neonicotinoid formaldehyde generators: possible mechanism of mouse-specific hepatotoxicity/hepatocarcinogenicity of thiamethoxam.

Thiamethoxam (TMX), an important insecticide, is hepatotoxic and hepatocarcinogenic in mice but not rats. Studies of Syngenta Central Toxicology Laboratory on species specificity in metabolism established that TMX is a much better substrate for mouse liver microsomal CYPs than the corresponding rat or human enzymes in forming desmethyl-TMX (dm-TMX), which is also hepatotoxic, and clothianidin (CLO), which is not hepatotoxic or hepatocarcinogenic. They proposed that TMX hepatotoxicity/hepatocarcinogencity is due to dm-TMX and a further metabolite desmethyl-CLO (dm-CLO) (structurally analogous to a standard inducible nitric oxide synthase inhibitor) acting synergistically. The present study considers formation of formaldehyde (HCHO) and N-methylol intermediates as an alternative mechanism of TMX hepatotoxicity/hepatocarcinogenicity. Comparison of neonicotinoid metabolism by mouse, rat and human microsomes with NADPH showed two important points. First, TMX and dm-TMX yield more HCHO than any other commercial neonicotinoid. Second, mouse microsomes give much higher conversion than rat or human microsomes. These observations provide an alternative hypothesis of HCHO and N-methylol intermediates from CYP-mediated oxidative oxadiazinane ring cleavage as the bioactivated hepatotoxicants. However, the proposed mono-N-methylol CYP metabolites are not observed, possibly further reacting in situ.

Neonicotinoid insecticides: oxidative stress in planta and metallo-oxidase inhibition.

Neonicotinoids not only control insect pests but also sometimes independently alter plant growth and response to stress. We find that imidacloprid, thiacloprid, acetamiprid, thiamethoxam, and clothianidin but not nitenpyram and dinotefuran induce foliar lesions and peroxidative damage in soybean ( Glycine max ) seedlings assayed with the 3,3'-diaminobenzidine stain. The chloropyridinyl-carboxylic acid (COOH) but not the -carboxaldehyde (CHO) metabolites induce peroxidative damage but in a different pattern. Surprisingly, the chlorothiazolyl -CHO and -COOH metabolites induce chlorosis but no clear superimposable peroxidative damage or cell death. Four metallo-oxidases known to modulate reactive oxygen species were not sensitive in vitro to the parent neonicotinoid itself but were to several CHO and COOH metabolites and related compounds, with a sensitivity order of CHO > COOH and tyrosinase > xanthine oxidase and aldehyde oxidase > catalase. Although metallo-oxidase inhibition does not correlate overall with lesion formation, it may play an as yet unknown role in plant response to neonicotinoids.