Perspective: use and reuse of NMR-based metabolomics data: what works and what remains challenging.

BACKGROUND: The National Cancer Institute issued a Request for Information (RFI; NOT-CA-23-007) in October 2022, soliciting input on using and reusing metabolomics data. This RFI aimed to gather input on best practices for metabolomics data storage, management, and use/reuse. AIM OF REVIEW: The nuclear magnetic resonance (NMR) Interest Group within the Metabolomics Association of North America (MANA) prepared a set of recommendations regarding the deposition, archiving, use, and reuse of NMR-based and, to a lesser extent, mass spectrometry (MS)-based metabolomics datasets. These recommendations were built on the collective experiences of metabolomics researchers within MANA who are generating, handling, and analyzing diverse metabolomics datasets spanning experimental (sample handling and preparation, NMR/MS metabolomics data acquisition, processing, and spectral analyses) to computational (automation of spectral processing, univariate and multivariate statistical analysis, metabolite prediction and identification, multi-omics data integration, etc.) studies. KEY SCIENTIFIC CONCEPTS OF REVIEW: We provide a synopsis of our collective view regarding the use and reuse of metabolomics data and articulate several recommendations regarding best practices, which are aimed at encouraging researchers to strengthen efforts toward maximizing the utility of metabolomics data, multi-omics data integration, and enhancing the overall scientific impact of metabolomics studies.

Mapping microhabitats of lignocellulose decomposition by a microbial consortium.

The leaf-cutter ant fungal garden ecosystem is a naturally evolved model system for efficient plant biomass degradation. Degradation processes mediated by the symbiotic fungus Leucoagaricus gongylophorus are difficult to characterize due to dynamic metabolisms and spatial complexity of the system. Herein, we performed microscale imaging across 12-µm-thick adjacent sections of Atta cephalotes fungal gardens and applied a metabolome-informed proteome imaging approach to map lignin degradation. This approach combines two spatial multiomics mass spectrometry modalities that enabled us to visualize colocalized metabolites and proteins across and through the fungal garden. Spatially profiled metabolites revealed an accumulation of lignin-related products, outlining morphologically unique lignin microhabitats. Metaproteomic analyses of these microhabitats revealed carbohydrate-degrading enzymes, indicating a prominent fungal role in lignocellulose decomposition. Integration of metabolome-informed proteome imaging data provides a comprehensive view of underlying biological pathways to inform our understanding of metabolic fungal pathways in plant matter degradation within the micrometer-scale environment.

Characterizing Families of Spectral Similarity Scores and Their Use Cases for Gas Chromatography-Mass Spectrometry Small Molecule Identification.

Metabolomics provides a unique snapshot into the world of small molecules and the complex biological processes that govern the human, animal, plant, and environmental ecosystems encapsulated by the One Health modeling framework. However, this "molecular snapshot" is only as informative as the number of metabolites confidently identified within it. The spectral similarity (SS) score is traditionally used to identify compound(s) in mass spectrometry approaches to metabolomics, where spectra are matched to reference libraries of candidate spectra. Unfortunately, there is little consensus on which of the dozens of available SS metrics should be used. This lack of standard SS score creates analytic uncertainty and potentially leads to issues in reproducibility, especially as these data are integrated across other domains. In this work, we use metabolomic spectral similarity as a case study to showcase the challenges in consistency within just one piece of the One Health framework that must be addressed to enable data science approaches for One Health problems. Here, using a large cohort of datasets comprising both standard and complex datasets with expert-verified truth annotations, we evaluated the effectiveness of 66 similarity metrics to delineate between correct matches (true positives) and incorrect matches (true negatives). We additionally characterize the families of these metrics to make informed recommendations for their use. Our results indicate that specific families of metrics (the Inner Product, Correlative, and Intersection families of scores) tend to perform better than others, with no single similarity metric performing optimally for all queried spectra. This work and its findings provide an empirically-based resource for researchers to use in their selection of similarity metrics for GC-MS identification, increasing scientific reproducibility through taking steps towards standardizing identification workflows.

Evaluating Retention Index Score Assumptions to Refine GC-MS Metabolite Identification.

As metabolomics grows into a high-throughput and high demand research field, current metrics for the identification of small molecules in gas chromatography-mass spectrometry (GC-MS) still require manual verification. Though steps have been taken to improve scoring metrics by combining spectral similarity (SS) and retention index (RI), the problem persists. A large body of literature has analyzed and refined SS scores, but few studies have explicitly studied improvements to RI scores. Here, we examined whether uninvestigated assumptions of the RI score are valid and propose ways to improve them. Query RIs were matched to library RI with a generous window of ±35 to avoid unintentional removal of valid compound identifications. Each match was manually verified as a true positive (TP), true negative, or unknown. Metabolites with at least 30 TP identifications were included in downstream analyses, resulting in a total of 87 metabolites from samples of varying complexity and type (e.g., amino acid mixtures, human urine, fungal species, and so on.). Our results showed that the RI score assumptions of normality, consistent variance across metabolites, and a mean error centered at 0 are often violated. We demonstrated through a cross-validation analysis that modifying these underlying assumptions according to empirical metabolite-specific distributions improved the TP and negative rankings. Further, we statistically determined the minimum number of samples required to estimate distributional parameters for scoring metrics. Overall, this work proposes a robust statistical pipeline to reduce the time bottleneck of metabolite identification by improving RI scores and thus minimize the effort to complete manual verification.

Gaussian Mixture Modeling Extensions for Improved False Discovery Rate Estimation in GC-MS Metabolomics.

The ability to reliably identify small molecules (e.g., metabolites) is key toward driving scientific advancement in metabolomics. Gas chromatography-mass spectrometry (GC-MS) is an analytic method that may be applied to facilitate this process. The typical GC-MS identification workflow involves quantifying the similarity of an observed sample spectrum and other features (e.g., retention index) to that of several references, noting the compound of the best-matching reference spectrum as the identified metabolite. While a deluge of similarity metrics exist, none quantify the error rate of generated identifications, thereby presenting an unknown risk of false identification or discovery. To quantify this unknown risk, we propose a model-based framework for estimating the false discovery rate (FDR) among a set of identifications. Extending a traditional mixture modeling framework, our method incorporates both similarity score and experimental information in estimating the FDR. We apply these models to identification lists derived from across 548 samples of varying complexity and sample type (e.g., fungal species, standard mixtures, etc.), comparing their performance to that of the traditional Gaussian mixture model (GMM). Through simulation, we additionally assess the impact of reference library size on the accuracy of FDR estimates. In comparing the best performing model extensions to the GMM, our results indicate relative decreases in median absolute estimation error (MAE) ranging from 12% to 70%, based on comparisons of the median MAEs across all hit-lists. Results indicate that these relative performance improvements generally hold despite library size; however FDR estimation error typically worsens as the set of reference compounds diminishes.

The Sugar Metabolic Model of Aspergillus niger Can Only Be Reliably Transferred to Fungi of Its Phylum.

Fungi play a critical role in the global carbon cycle by degrading plant polysaccharides to small sugars and metabolizing them as carbon and energy sources. We mapped the well-established sugar metabolic network of Aspergillus niger to five taxonomically distant species (Aspergillus nidulans, Penicillium subrubescens, Trichoderma reesei, Phanerochaete chrysosporium and Dichomitus squalens) using an orthology-based approach. The diversity of sugar metabolism correlates well with the taxonomic distance of the fungi. The pathways are highly conserved between the three studied Eurotiomycetes (A. niger, A. nidulans, P. subrubescens). A higher level of diversity was observed between the T. reesei and A. niger, and even more so for the two Basidiomycetes. These results were confirmed by integrative analysis of transcriptome, proteome and metabolome, as well as growth profiles of the fungi growing on the corresponding sugars. In conclusion, the establishment of sugar pathway models in different fungi revealed the diversity of fungal sugar conversion and provided a valuable resource for the community, which would facilitate rational metabolic engineering of these fungi as microbial cell factories.

Light-Dependence of Formate (C1) and Acetate (C2) Transport and Oxidation in Poplar Trees.

Although apparent light inhibition of leaf day respiration is a widespread reported phenomenon, the mechanisms involved, including utilization of alternate respiratory pathways and substrates and light inhibition of TCA cycle enzymes are under active investigation. Recently, acetate fermentation was highlighted as a key drought survival strategy mediated through protein acetylation and jasmonate signaling. Here, we evaluate the light-dependence of acetate transport and assimilation in Populus trichocarpa trees using the dynamic xylem solution injection (DXSI) method developed here for continuous studies of C1 and C2 organic acid transport and light-dependent metabolism. Over 7 days, 1.0 L of [13C]formate and [13C2]acetate solutions were delivered to the stem base of 2-year old potted poplar trees, while continuous diurnal observations were made in the canopy of CO2, H2O, and isoprene gas exchange together with δ13CO2. Stem base injection of 10 mM [13C2]acetate induced an overall pattern of canopy branch headspace 13CO2 enrichment (δ13CO2 +27‱) with a diurnal structure in δ13CO2 reaching a mid-day minimum followed by a maximum shortly after darkening where δ13CO2 values rapidly increased up to +12‱. In contrast, 50 mM injections of [13C]formate were required to reach similar δ13CO2 enrichment levels in the canopy with δ13CO2 following diurnal patterns of transpiration. Illuminated leaves of detached poplar branches pretreated with 10 mM [13C2]acetate showed lower δ13CO2 (+20‱) compared to leaves treated with 10 mM [13C]formate (+320‱), the opposite pattern observed at the whole plant scale. Following dark/light cycles at the leaf-scale, rapid, strong, and reversible enhancements in headspace δ13CO2 by up to +60‱ were observed in [13C2]acetate-treated leaves which showed enhanced dihydrojasmonic acid and TCA cycle intermediate concentrations. The results are consistent with acetate in the transpiration stream as an effective activator of the jasmonate signaling pathway and respiratory substrate. The shorter lifetime of formate relative to acetate in the transpiration stream suggests rapid formate oxidation to CO2 during transport to the canopy. In contrast, acetate is efficiently transported to the canopy where an increased allocation towards mitochondrial dark respiration occurs at night. The results highlight the potential for an effective integration of acetate into glyoxylate and TCA cycles and the light-inhibition of citrate synthase as a potential regulatory mechanism controlling the diurnal allocation of acetate between anabolic and catabolic processes.

Root exudate composition reflects drought severity gradient in blue grama (Bouteloua gracilis).

Plant survival during environmental stress greatly affects ecosystem carbon (C) cycling, and plant-microbe interactions are central to plant stress survival. The release of C-rich root exudates is a key mechanism plants use to manage their microbiome, attracting beneficial microbes and/or suppressing harmful microbes to help plants withstand environmental stress. However, a critical knowledge gap is how plants alter root exudate concentration and composition under varying stress levels. In a greenhouse study, we imposed three drought treatments (control, mild, severe) on blue grama (Bouteloua gracilis Kunth Lag. Ex Griffiths), and measured plant physiology and root exudate concentration and composition using GC-MS, NMR, and FTICR. With increasing drought severity, root exudate total C and organic C increased concurrently with declining predawn leaf water potential and photosynthesis. Root exudate composition mirrored the physiological gradient of drought severity treatments. Specific compounds that are known to alter plant drought responses and the rhizosphere microbiome mirrored the drought severity-induced root exudate compositional gradient. Despite reducing C uptake, these plants actively invested C to root exudates with increasing drought severity. Patterns of plant physiology and root exudate concentration and composition co-varied along a gradient of drought severity.

Re-routing of Sugar Catabolism Provides a Better Insight Into Fungal Flexibility in Using Plant Biomass-Derived Monomers as Substrates.

The filamentous ascomycete Aspergillus niger has received increasing interest as a cell factory, being able to efficiently degrade plant cell wall polysaccharides as well as having an extensive metabolism to convert the released monosaccharides into value added compounds. The pentoses D-xylose and L-arabinose are the most abundant monosaccharides in plant biomass after the hexose D-glucose, being major constituents of xylan, pectin and xyloglucan. In this study, the influence of selected pentose catabolic pathway (PCP) deletion strains on growth on plant biomass and re-routing of sugar catabolism was addressed to gain a better understanding of the flexibility of this fungus in using plant biomass-derived monomers. The transcriptome, metabolome and proteome response of three PCP mutant strains, ΔlarAΔxyrAΔxyrB, ΔladAΔxdhAΔsdhA and ΔxkiA, grown on wheat bran (WB) and sugar beet pulp (SBP), was evaluated. Our results showed that despite the absolute impact of these PCP mutations on pure pentose sugars, they are not as critical for growth of A. niger on more complex biomass substrates, such as WB and SBP. However, significant phenotypic variation was observed between the two biomass substrates, but also between the different PCP mutants. This shows that the high sugar heterogeneity of these substrates in combination with the high complexity and adaptability of the fungal sugar metabolism allow for activation of alternative strategies to support growth.

Preoperative Metabolic Signatures of Prostate Cancer Recurrence Following Radical Prostatectomy.

Technological advances in mass spectrometry (MS), liquid chromatography (LC) separations, nuclear magnetic resonance (NMR) spectroscopy, and big data analytics have made possible studying metabolism at an "omics" or systems level. Here, we applied a multiplatform (NMR + LC-MS) metabolomics approach to the study of preoperative metabolic alterations associated with prostate cancer recurrence. Thus far, predicting which patients will recur even after radical prostatectomy has not been possible. Correlation analysis on metabolite abundances detected on serum samples collected prior to surgery from prostate cancer patients ( n = 40 remission vs n = 40 recurrence) showed significant alterations in a number of pathways, including amino acid metabolism, purine and pyrimidine synthesis, tricarboxylic acid cycle, tryptophan catabolism, glucose, and lactate. Lipidomics experiments indicated higher lipid abundances on recurrent patients for a number of classes that included triglycerides, lysophosphatidylcholines, phosphatidylethanolamines, phosphatidylinositols, diglycerides, acyl carnitines, and ceramides. Machine learning approaches led to the selection of a 20-metabolite panel from a single preoperative blood sample that enabled prediction of recurrence with 92.6% accuracy, 94.4% sensitivity, and 91.9% specificity under cross-validation conditions.

Ambient mass spectrometry in metabolomics.

Since the introduction of desorption electrospray ionization (DESI) mass spectrometry (MS), ambient MS methods have seen increased use in a variety of fields from health to food science. Increasing its popularity in metabolomics, ambient MS offers limited sample preparation, rapid and direct analysis of liquids, solids, and gases, in situ and in vivo analysis, and imaging. The metabolome consists of a constantly changing collection of small (<1.5 kDa) molecules. These include endogenous molecules that are part of primary metabolism pathways, secondary metabolites with specific functions such as signaling, chemicals incorporated in the diet or resulting from environmental exposures, and metabolites associated with the microbiome. Characterization of the responsive changes of this molecule cohort is the principal goal of any metabolomics study. With adjustments to experimental parameters, metabolites with a range of chemical and physical properties can be selectively desorbed and ionized and subsequently analyzed with increased speed and sensitivity. This review covers the broad applications of a variety of ambient MS techniques in four primary fields in which metabolomics is commonly employed.

13C Metabolomics: NMR and IROA for Unknown Identification.

ABSTRACT BACKGROUND: Isotopic Ratio Outlier Analysis (IROA) is an untargeted metabolomics method that uses stable isotopic labeling and LC-HRMS for identification and relative quantification of metabolites in a biological sample under varying experimental conditions. OBJECTIVE: We demonstrate a method using high-sensitivity 13C NMR to identify an unknown metabolite isolated from fractionated material from an IROA LC-HRMS experiment. METHODS: IROA samples from the nematode Caenorhabditis elegans were fractionated using LC-HRMS using 5 repeated injections and collecting 30 sec fractions. These were concentrated and analyzed by 13C NMR. RESULTS: We isotopically labeled samples of C. elegans and collected 2 adjacent LC fractions. By HRMS, one contained at least 2 known metabolites, phenylalanine and inosine, and the other contained tryptophan and an unknown feature with a monoisotopic mass of m/z 380.0742 [M+H]+. With NMR, we were able to easily verify the known compounds, and we then identified the spin system networks responsible for the unknown resonances. After searching the BMRB database and comparing the molecular formula from LC-HRMS, we determined that the fragments were a modified anthranilate and a glucose modified by a phosphate. We then performed quantum chemical NMR chemical shift calculations to determine the most likely isomer, which was 3'-O-phospho-ß-D-glucopyranosyl-anthranilate. This compound had previously been found in the same organism, validating our approach. CONCLUSION: We were able to dereplicate previously known metabolites and identify a metabolite that was not in databases by matching resonances to NMR databases and using chemical shift calculations to determine the correct isomer. This approach is efficient and can be used to identify unknown compounds of interest using the same material used for IROA.

An overview of methods using (13)C for improved compound identification in metabolomics and natural products.

Compound identification is a major bottleneck in metabolomics studies. In nuclear magnetic resonance (NMR) investigations, resonance overlap often hinders unambiguous database matching or de novo compound identification. In liquid chromatography-mass spectrometry (LC-MS), discriminating between biological signals and background artifacts and reliable determination of molecular formulae are not always straightforward. We have designed and implemented several NMR and LC-MS approaches that utilize (13)C, either enriched or at natural abundance, in metabolomics applications. For LC-MS applications, we describe a technique called isotopic ratio outlier analysis (IROA), which utilizes samples that are isotopically labeled with 5% (test) and 95% (control) (13)C. This labeling strategy leads to characteristic isotopic patterns that allow the differentiation of biological signals from artifacts and yield the exact number of carbons, significantly reducing possible molecular formulae. The relative abundance between the test and control samples for every IROA feature can be determined simply by integrating the peaks that arise from the 5 and 95% channels. For NMR applications, we describe two (13)C-based approaches. For samples at natural abundance, we have developed a workflow to obtain (13)C-(13)C and (13)C-(1)H statistical correlations using 1D (13)C and (1)H NMR spectra. For samples that can be isotopically labeled, we describe another NMR approach to obtain direct (13)C-(13)C spectroscopic correlations. These methods both provide extensive information about the carbon framework of compounds in the mixture for either database matching or de novo compound identification. We also discuss strategies in which (13)C NMR can be used to identify unknown compounds from IROA experiments. By combining technologies with the same samples, we can identify important biomarkers and corresponding metabolites of interest.

Metabolomics and Natural-Products Strategies to Study Chemical Ecology in Nematodes.

This review provides an overview of two complementary approaches to identify biologically active compounds for studies in chemical ecology. The first is activity-guided fractionation and the second is metabolomics, particularly focusing on a new liquid chromatography-mass spectrometry-based method called isotopic ratio outlier analysis. To illustrate examples using these approaches, we review recent experiments using Caenorhabditis elegans and related free-living nematodes.

(13)C NMR Metabolomics: INADEQUATE Network Analysis.

The many advantages of (13)C NMR are often overshadowed by its intrinsically low sensitivity. Given that carbon makes up the backbone of most biologically relevant molecules, (13)C NMR offers a straightforward measurement of these compounds. Two-dimensional (13)C-(13)C correlation experiments like INADEQUATE (incredible natural abundance double quantum transfer experiment) are ideal for the structural elucidation of natural products and have great but untapped potential for metabolomics analysis. We demonstrate a new and semiautomated approach called INETA (INADEQUATE network analysis) for the untargeted analysis of INADEQUATE data sets using an in silico INADEQUATE database. We demonstrate this approach using isotopically labeled Caenorhabditis elegans mixtures.

¹³C NMR metabolomics: applications at natural abundance.

(13)C NMR has many advantages for a metabolomics study, including a large spectral dispersion, narrow singlets at natural abundance, and a direct measure of the backbone structures of metabolites. However, it has not had widespread use because of its relatively low sensitivity compounded by low natural abundance. Here we demonstrate the utility of high-quality (13)C NMR spectra obtained using a custom (13)C-optimized probe on metabolomic mixtures. A workflow was developed to use statistical correlations between replicate 1D (13)C and (1)H spectra, leading to composite spin systems that can be used to search publicly available databases for compound identification. This was developed using synthetic mixtures and then applied to two biological samples, Drosophila melanogaster extracts and mouse serum. Using the synthetic mixtures we were able to obtain useful (13)C-(13)C statistical correlations from metabolites with as little as 60 nmol of material. The lower limit of (13)C NMR detection under our experimental conditions is approximately 40 nmol, slightly lower than the requirement for statistical analysis. The (13)C and (1)H data together led to 15 matches in the database compared to just 7 using (1)H alone, and the (13)C correlated peak lists had far fewer false positives than the (1)H generated lists. In addition, the (13)C 1D data provided improved metabolite identification and separation of biologically distinct groups using multivariate statistical analysis in the D. melanogaster extracts and mouse serum.

Isotopic ratio outlier analysis global metabolomics of Caenorhabditis elegans.

We demonstrate the global metabolic analysis of Caenorhabditis elegans stress responses using a mass-spectrometry-based technique called isotopic ratio outlier analysis (IROA). In an IROA protocol, control and experimental samples are isotopically labeled with 95 and 5% (13)C, and the two sample populations are mixed together for uniform extraction, sample preparation, and LC-MS analysis. This labeling strategy provides several advantages over conventional approaches: (1) compounds arising from biosynthesis are easily distinguished from artifacts, (2) errors from sample extraction and preparation are minimized because the control and experiment are combined into a single sample, (3) measurement of both the molecular weight and the exact number of carbon atoms in each molecule provides extremely accurate molecular formulas, and (4) relative concentrations of all metabolites are easily determined. A heat-shock perturbation was conducted on C. elegans to demonstrate this approach. We identified many compounds that significantly changed upon heat shock, including several from the purine metabolism pathway. The metabolomic response information by IROA may be interpreted in the context of a wealth of genetic and proteomic information available for C. elegans . Furthermore, the IROA protocol can be applied to any organism that can be isotopically labeled, making it a powerful new tool in a global metabolomics pipeline.