Chemical Discovery in the Era of Metabolomics

A pilot study for fragment identification using 2D NMR and deep learning

This paper presents a proof of concept of a method to identify substructures in 2D NMR spectra of mixtures using a bespoke image-based convolutional neural network application. This is done using HSQC and HMBC spectra separately and in combination. The application can reliably detect substructures in pure compounds, using a simple network. Results indicate that it can work for mixtures when trained on pure compounds only. HMBC data and the combination of HMBC and HSQC show better results than HSQC alone in this pilot study.

Identification of Exhaled Metabolites in Children with Cystic Fibrosis

The early detection of inflammation and infection is important to prevent irreversible lung damage in cystic fibrosis. Novel and non-invasive monitoring tools would be of high benefit for the quality of life of patients. Our group previously detected over 100 exhaled mass-to-charge (m/z) features, using on-line secondary electrospray ionization high-resolution mass spectrometry (SESI-HRMS), which distinguish children with cystic fibrosis from healthy controls. The aim of this study was to annotate as many m/z features as possible with putative chemical structures. Compound identification was performed by applying a rigorous workflow, which included the analysis of on-line MS2 spectra and a literature comparison. A total of 49 discriminatory exhaled compounds were putatively identified. A group of compounds including glycolic acid, glyceric acid and xanthine were elevated in the cystic fibrosis group. A large group of acylcarnitines and aldehydes were found to be decreased in cystic fibrosis. The proposed compound identification workflow was used to identify signatures of volatile organic compounds that discriminate children with cystic fibrosis from healthy controls, which is the first step for future non-invasive and personalized applications.

Supercritical Fluid Nanospray Mass Spectrometry

Supercritical fluids are typically electrosprayed using an organic solvent makeup flow to facilitate continuous electrical connection and enhancement of electrospray stability. This results in sample dilution, loss in sensitivity, and potential phase separation. Premixing the supercritical fluid with organic solvent has shown substantial benefits to electrospray efficiency and increased analyte charge state. Presented here is a nanospray mass spectrometry system for supercritical fluids (nSF-MS). This split flow system used small i.d. capillaries, heated interface, inline frit, and submicron emitter tips to electrospray quaternary alkyl amines solvated in supercritical CO₂ with a 10% methanol modifier. Analyte signal response was evaluated as a function of total system flow rate (0.5-1.5 mL/min) that is split to nanospray a supercritical fluid with linear flow rates between 0.07 and 0.42 cm/sec and pressure ranges (15-25 MPa). The nSF system showed mass-sensitive detection based on increased signal intensity for increasing capillary i.d. and analyte injection volume. These effects indicate efficient solvent evaporation for the analysis of quaternary amines. Carrier additives generally decreased signal intensity. Comparison of the nSF-MS system to the conventional SF makeup flow ESI showed 10-fold signal intensity enhancement across all the capillary i.d.s. The nSF-MS system likely achieves rapid solvent evaporation of the SF at the emitter point. The developed system combined the benefits of the nanoemitters, sCO₂, and the low modifier percentage which gave rise to enhancement in MS detection sensitivity.

Addressing big data challenges in mass spectrometry-based metabolomics

Advancements in computer science and software engineering have greatly facilitated mass spectrometry (MS)-based untargeted metabolomics. Nowadays, gigabytes of metabolomics data are routinely generated from MS platforms, containing condensed structural and quantitative information from thousands of metabolites. Manual data processing is almost impossible due to the large data size. Therefore, in the "omics" era, we are faced with new challenges, the big data challenges of how to accurately and efficiently process the raw data, extract the biological information, and visualize the results from the gigantic amount of collected data. Although important, proposing solutions to address these big data challenges requires broad interdisciplinary knowledge, which can be challenging for many metabolomics practitioners. Our laboratory in the Department of Chemistry at the University of British Columbia is committed to combining analytical chemistry, computer science, and statistics to develop bioinformatics tools that address these big data challenges. In this Feature Article, we elaborate on the major big data challenges in metabolomics, including data acquisition, feature extraction, quantitative measurements, statistical analysis, and metabolite annotation. We also introduce our recently developed bioinformatics solutions for these challenges. Notably, all of the bioinformatics tools and source codes are freely available on GitHub (https://www.github.com/HuanLab), along with revised and regularly updated content.

Using MetaboAnalyst 5.0 for LC-HRMS spectra processing, multi-omics integration and covariate adjustment of global metabolomics data

Liquid chromatography coupled with high-resolution mass spectrometry (LC-HRMS) has become a workhorse in global metabolomics studies with growing applications across biomedical and environmental sciences. However, outstanding bioinformatics challenges in terms of data processing, statistical analysis and functional interpretation remain critical barriers to the wider adoption of this technology. To help the user community overcome these barriers, we have made major updates to the well-established MetaboAnalyst platform ( www.metaboanalyst.ca ). This protocol extends the previous 2011 Nature Protocol by providing stepwise instructions on how to use MetaboAnalyst 5.0 to: optimize parameters for LC-HRMS spectra processing; obtain functional insights from peak list data; integrate metabolomics data with transcriptomics data or combine multiple metabolomics datasets; conduct exploratory statistical analysis with complex metadata. Parameter optimization may take ~2 h to complete depending on the server load, and the remaining three stages may be executed in ~60 min.

AI/ML-driven advances in untargeted metabolomics and exposomics for biomedical applications

Metabolomics describes a high-throughput approach for measuring a repertoire of metabolites and small molecules in biological samples. One utility of untargeted metabolomics, unbiased global analysis of the metabolome, is to detect key metabolites as contributors to, or readouts of, human health and disease. In this perspective, we discuss how artificial intelligence (AI) and machine learning (ML) have promoted major advances in untargeted metabolomics workflows and facilitated pivotal findings in the areas of disease screening and diagnosis. We contextualize applications of AI and ML to the emerging field of high-resolution mass spectrometry (HRMS) exposomics, which unbiasedly detects endogenous metabolites and exogenous chemicals in human tissue to characterize exposure linked with disease outcomes. We discuss the state of the science and suggest potential opportunities for using AI and ML to improve data quality, rigor, detection, and chemical identification in untargeted metabolomics and exposomics studies.

Improvement of Biosynthetic Ansamitocin P-3 Production Based on Oxygen-Vector Screening and Metabonomics Analysis

A novel approach involving exogenous oxygen vectors was developed for improving the production of biosynthetic Ansamitocin P-3 (AP-3). Four types of oxygen vectors including soybean oil, n-dodecane, n-hexadecane, and Tween-80 were applied to explore the effect of exogenous oxygen vectors on AP-3 yield. It was observed that soybean oil exhibited a better ability for promoting AP-3 generation than the other three oxygen vectors. Based on the results of the single-factor experiment, response surface methodology was employed to obtain the optimal soybean oil addition method. The optimum soybean oil concentration was 0.52%, and the addition time was 50 h. Under this condition, the yield of AP-3 reached 106.04 mg/L, which was 49.48% higher than that of the control group without adding oxygen vectors. To further investigate the influence of dissolved oxygen on precious orange tufts actinomycetes variety A. pretiosum strain metabolism and AP-3 yield, metabolomics analysis was carried out by detecting strain intermediate metabolites at various stages under different dissolved oxygen levels. Moreover, differential metabolite screening and metabolic pathway enrichment analysis were combined to exploit the effect mechanism of soybean oil on AP-3 production. Results suggested that primary metabolic levels of the TCA cycle and amino acid metabolism increased with the increase in dissolved oxygen level, which was beneficial to the life activities of bacteria and the synthesis of secondary metabolic precursors, thus increasing the production of AP-3.

N-Acetyl-Cysteinylated Streptophenazines from Streptomyces

Here, we describe two N-acetyl-cysteinylated streptophenazines (1 and 2) produced by the soil-derived Streptomyces sp. ID63040 and identified through a metabolomic approach. These metabolites attracted our interest due to their low occurrence frequency in a large library of fermentation broth extracts and their consistent presence in biological replicates of the producer strain. The compounds were found to possess broad-spectrum antibacterial activity while exhibiting low cytotoxicity. The biosynthetic gene cluster from Streptomyces sp. ID63040 was found to be highly similar to the streptophenazine reference cluster in the MIBiG database, which originates from the marine Streptomyces sp. CNB-091. Compounds 1 and 2 were the main streptophenazine products from Streptomyces sp. ID63040 at all cultivation times but were not detected in Streptomyces sp. CNB-091. The lack of obvious candidates for cysteinylation in the Streptomyces sp. ID63040 biosynthetic gene cluster suggests that the N-acetyl-cysteine moiety derives from cellular functions, most likely from mycothiol. Overall, our data represent an interesting example of how to leverage metabolomics for the discovery of new natural products and point out the often-neglected contribution of house-keeping cellular functions to natural product diversification.

MAFFIN: Metabolomics Sample Normalization Using Maximal Density Fold Change with High-Quality Metabolic Features and Corrected Signal Intensities

MOTIVATION

Post-acquisition sample normalization is a critical step in comparative metabolomics to remove the variation introduced by sample amount or concentration difference. Previously reported approaches are either specific to one sample type or built on strong assumptions on data structure, which are limited to certain levels. This encouraged us to develop MAFFIN, an accurate and robust post-acquisition sample normalization workflow that works universally for metabolomics data collected on mass spectrometry (MS) platforms.

RESULTS

MAFFIN calculates normalization factors using maximal density fold change (MDFC) computed by a kernel density-based approach. Using both simulated data and 20 metabolomics data sets, we showcased that MDFC outperforms four commonly used normalization methods in terms of reducing the intragroup variation among samples. Two essential steps, overlooked in conventional methods, were also examined and incorporated into MAFFIN. (1) MAFFIN uses multiple orthogonal criteria to select high-quality features for normalization factor calculation, which minimizes the bias caused by abiotic features or metabolites with poor quantitative performance. (2) MAFFIN corrects the MS signal intensities of high-quality features using serial quality control (QC) samples, which guarantees the accuracy of fold change calculations. MAFFIN was applied to a human saliva metabolomics study and led to better data separation in principal component analysis (PCA) and more confirmed significantly altered metabolites.

AVAILABILITY AND IMPLEMENTATION

The MAFFIN algorithm was implemented in an R package named MAFFIN. Package installation, user instruction, and demo data are available at https://github.com/HuanLab/MAFFIN.

SUPPLEMENTARY INFORMATION

Supplementary data are available at Bioinformatics online.

MINE 2.0: Enhanced biochemical coverage for peak identification in untargeted metabolomics

SUMMARY

Although advances in untargeted metabolomics have made it possible to gather data on thousands of cellular metabolites in parallel, identification of novel metabolites from these datasets remains challenging. To address this need, Metabolic in silico Network Expansions (MINEs) were developed. A MINE is an expansion of known biochemistry which can be used as a list of potential structures for unannotated metabolomics peaks. Here, we present MINE 2.0, which utilizes a new set of biochemical transformation rules that covers 93% of MetaCyc reactions (compared to 25% in MINE 1.0). This results in a 17-fold increase in database size and a 40% increase in MINE database compounds matching unannotated peaks from an untargeted metabolomics dataset. MINE 2.0 is thus a significant improvement to this community resource.

AVAILABILITY AND IMPLEMENTATION

The MINE 2.0 website can be accessed at https://minedatabase.ci.northwestern.edu. The MINE 2.0 web API documentation can be accessed at https://mine-api.readthedocs.io/en/latest/. The data and code underlying this article are available in the MINE-2.0-Paper repository at https://github.com/tyo-nu/MINE-2.0-Paper. MINE 2.0 source code can be accessed at https://github.com/tyo-nu/MINE-Database (MINE construction), https://github.com/tyo-nu/MINE-Server (backend web API), and https://github.com/tyo-nu/MINE-app (web app).

SUPPLEMENTARY INFORMATION

Supplementary data are available at Bioinformatics online.

Resolving isobaric interferences in direct infusion tandem mass spectrometry

RATIONALE

The co-fragmentation of precursors in direct infusion (DI) tandem high-resolution mass spectrometry (HRMS) can complicate the fragment spectra and consequently lead to false hits during compound identification.

METHODS

The method herein described, termed IQAROS (incremental quadrupole acquisition to resolve overlapping spectra), modulates the intensities of precursors and fragments by stepwise movement of the quadrupole isolation window over the mass-to-charge (m/z) range of the precursors. The modulated signals are then deconvoluted by a linear regression model to reconstruct the fragment spectra with less interference. The hardware to demonstrate the use of IQAROS was an orbitrap with electrospray ionization (ESI) or secondary electrospray ionization (SESI), although the method can also be applied to other ionization techniques or mass analyzers.

RESULTS

Assessing the performance of IQAROS with isobaric standards revealed that the reconstructed fragment spectra match with spectra acquired from the pure standards and that more compounds were correctly identified compared with the classical approach with the quadrupole centered at the m/z value of the precursor of interest. Moreover, the strength of IQAROS is exemplified by the identification of two isobaric biomarkers directly from a breath sample with SESI-HRMS.

CONCLUSIONS

With IQAROS, cleaner fragment spectra of co-fragmenting isobars during DI-HRMS analysis can be obtained. IQAROS can easily be set up by the standard graphical user interface of the instrument. Therefore, it facilitates the characterization of features of interest in samples analyzed by DI-HRMS, for example, in high-throughput or real-time metabolomics.

Improving confidence in lipidomic annotations by incorporating empirical ion mobility regression analysis and chemical class prediction

MOTIVATION

Mass spectrometry-based untargeted lipidomics aims to globally characterize the lipids and lipid-like molecules in biological systems. Ion mobility increases coverage and confidence by offering an additional dimension of separation and a highly reproducible metric for feature annotation, the collision cross-section (CCS).

RESULTS

We present a data processing workflow to increase confidence in molecular class annotations based on CCS values. This approach uses class-specific regression models built from a standardized CCS repository (the Unified CCS Compendium) in a parallel scheme that combines a new annotation filtering approach with a machine learning class prediction strategy. In a proof-of-concept study using murine brain lipid extracts, 883 lipids were assigned higher confidence identifications using the filtering approach, which reduced the tentative candidate lists by over 50% on average. An additional 192 unannotated compounds were assigned a predicted chemical class.

AVAILABILITY AND IMPLEMENTATION

All relevant source code is available at https://github.com/McLeanResearchGroup/CCS-filter.

SUPPLEMENTARY INFORMATION

Supplementary data are available at Bioinformatics online.

Global Tracking of Transformation Products of Environmental Contaminants by 2H-Labeled Stable Isotope-Assisted Metabolomics

Stable isotope-assisted metabolomics (SIAM) enables global tracking of isotopic labels in nontargeted metabolomics in living organisms. However, its application in tracking transformation products (TPs, as metabolites of contaminants) of environmental contaminants is still a challenge due to limits in methodology, unmatured algorithms, and the high cost of ¹³C-labeled contaminants. Therefore, we developed a ²H-SIAM pipeline coupled with a highly flexible algorithm ²H-SIAM(1.0) (https://github.com/kechen1984/2H-SIAM), facilitating tracking TPs of contaminants in the environmental matrix. A detailed discussion illustrates the theory, behavior, and prospect of ²H-SIAM. We demonstrate that the proposed ²H-SIAM pipeline has unique advantages over ¹³C-SIAM, for example, cost-effective ²H-labeled contaminants, easy synthesis of ²H-labeled emerging contaminants, and providing more structural information. A pyrene soil degradation study further confirmed its high performance. It efficiently discarded 99% of noise signals and extracted 52 features from the nontargeted high resolution mass spectrometry (HRMS) data. Among them, 13 features were annotated as TPs of pyrene with identification confidence from Level 2a to Level 5, and 5 TPs were reported for the first time. In conclusion, the proposed ²H-SIAM pipeline is powerful in tracking potential TPs of environmental contaminants with unique advantages.

The Human Meconium Metabolome and Its Evolution during the First Days of Life

Meconium represents the first newborn stools, formed from the second month of gestation and excreted in the first days after birth. As an accumulative and inert matrix, it accumulates most of the molecules transferred through the placenta from the mother to the fetus during the last 6 months of pregnancy, and those resulting from the metabolic activities of the fetus. To date, only few studies dealing with meconium metabolomics have been published. In this study, we aimed to provide a comprehensive view of the meconium metabolic composition using 33 samples collected longitudinally from 11 healthy newborns and to analyze its evolution during the first 3 days of life. First, a robust and efficient methodology for metabolite extraction was implemented. Data acquisition was performed using liquid chromatography coupled to high-resolution mass spectrometry (LC-HRMS), using two complementary LC-HRMS conditions. Data preprocessing and treatment were performed using the Workflow4Metabolomics platform and the metabolite annotation was performed using our in-house database by matching accurate masses, retention times, and MS/MS spectra to those of pure standards. We successfully identified up to 229 metabolites at a high confidence level in human meconium, belonging to diverse chemical classes and from different origins. A progressive evolution of the metabolic profile was statistically evidenced, with sugars, amino acids, and some bacteria-derived metabolites being among the most impacted identified compounds. Our implemented analytical workflow allows a unique and comprehensive description of the meconium metabolome, which is related to factors, such as maternal diet and environment.

Paramounter: Direct Measurement of Universal Parameters To Process Metabolomics Data in a "White Box"

Choosing appropriate data processing parameters is critical in processing liquid chromatography-mass spectrometry (LC-MS)-based untargeted metabolomics data. The conventional design of experiments (DOE) approach is time-consuming and provides no intuitive explanation why the selected parameters generate the best results. After studying commonly used metabolomics data processing software, this work summarized a set of universal parameters, including mass tolerance, peak height, peak width, and instrumental shift. These universal parameters are shared among different feature extraction programs and are critical to metabolic feature extraction. We then developed Paramounter, an R program that automatically measures these universal parameters from raw LC-MS-based metabolomics data prior to metabolic feature extraction. This is made possible through novel concepts of rank-based intensity sorting, zone of interest, and many others. Paramounter also translates universal parameters to software-specific parameters for data processing in different programs. Applying Paramounter is demonstrated to provide a threefold increase in the extracted metabolites compared to using default parameters in MS-DIAL-based feature extraction. Furthermore, the comparison between Paramounter, AutoTuner, and IPO showed that Paramounter generates 3.7- and 1.6-fold more true positive features than AutoTuner and IPO, respectively. Further validation of Paramounter on 11 datasets covering different sample types, data acquisition modes, and MS vendors proved that Paramounter is a convenient and robust program. Overall, the proposed universal parameters and the development of Paramounter address a critical need in metabolomics data processing, transforming metabolomics feature extraction from a "black box" to a "white box." Paramounter is freely available on GitHub (https://github.com/HuanLab/Paramounter).

Reproducible untargeted metabolomics workflow for exhaustive MS2 data acquisition of MS1 features

Unknown features in untargeted metabolomics and non-targeted analysis (NTA) are identified using fragment ions from MS/MS spectra to predict the structures of the unknown compounds. The precursor ion selected for fragmentation is commonly performed using data dependent acquisition (DDA) strategies or following statistical analysis using targeted MS/MS approaches. However, the selected precursor ions from DDA only cover a biased subset of the peaks or features found in full scan data. In addition, different statistical analysis can select different precursor ions for MS/MS analysis, which make the post-hoc validation of ions selected following a secondary analysis impossible for precursor ions selected by the original statistical method. Here we propose an automated, exhaustive, statistical model-free workflow: paired mass distance-dependent analysis (PMDDA), for reproducible untargeted mass spectrometry MS2 fragment ion collection of unknown compounds found in MS1 full scan. Our workflow first removes redundant peaks from MS1 data and then exports a list of precursor ions for pseudo-targeted MS/MS analysis on independent peaks. This workflow provides comprehensive coverage of MS2 collection on unknown compounds found in full scan analysis using a “one peak for one compound” workflow without a priori redundant peak information. We compared pseudo-spectra formation and the number of MS2 spectra linked to MS1 data using the PMDDA workflow to that obtained using CAMERA and RAMclustR algorithms. More annotated compounds, molecular networks, and unique MS/MS spectra were found using PMDDA compared with CAMERA and RAMClustR. In addition, PMDDA can generate a preferred ion list for iterative DDA to enhance coverage of compounds when instruments support such functions. Finally, compounds with signals in both positive and negative modes can be identified by the PMDDA workflow, to further reduce redundancies. The whole workflow is fully reproducible as a docker image xcmsrocker with both the original data and the data processing template.

Comparative metabolomics with Metaboseek reveals functions of a conserved fat metabolism pathway in C. elegans

Untargeted metabolomics via high-resolution mass spectrometry can reveal more than 100,000 molecular features in a single sample, many of which may represent unidentified metabolites, posing significant challenges to data analysis. We here introduce Metaboseek, an open-source analysis platform designed for untargeted comparative metabolomics and demonstrate its utility by uncovering biosynthetic functions of a conserved fat metabolism pathway, α-oxidation, using C. elegans as a model. Metaboseek integrates modules for molecular feature detection, statistics, molecular formula prediction, and fragmentation analysis, which uncovers more than 200 previously uncharacterized α-oxidation-dependent metabolites in an untargeted comparison of wildtype and α-oxidation-defective hacl-1 mutants. The identified metabolites support the predicted enzymatic function of HACL-1 and reveal that α-oxidation participates in metabolism of endogenous β-methyl-branched fatty acids and food-derived cyclopropane lipids. Our results showcase compound discovery and feature annotation at scale via untargeted comparative metabolomics applied to a conserved primary metabolic pathway and suggest a model for the metabolism of cyclopropane lipids.

Untargeted mass spectrometry-based metabolomics can reveal new biochemistry, but data analysis is challenging. Here, the authors develop Metaboseek, an open-source software that facilitates metabolite discovery, and apply it to characterize fatty acid alpha-oxidation in C. elegans.

Ultrahigh-Resolution NMR with Water Signal Suppression for a Deeper Understanding of the Action of Antimetabolic Drugs on Diffuse Large B-Cell Lymphoma

Ultrahigh-resolution NMR has recently attracted considerable attention in the field of complex samples analysis. Indeed, the implementation of broadband homonuclear decoupling techniques has allowed us to greatly simplify crowded ¹H spectra, yielding singlets for almost every proton site from the analyzed molecules. Pure shift methods have notably shown to be particularly suitable for deciphering mixtures of metabolites in biological samples. Here, we have successfully implemented a new pure shift pulse sequence based on the PSYCHE method, which incorporates a block for solvent suppression that is suitable for metabolomics analysis. The resulting experiment allows us to record ultrahigh-resolution 1D NOESY ¹H spectra of biofluids with suppression of the water signal, which is a crucial step for highlighting metabolite mixtures in an aqueous phase. We have successfully recorded pure shift spectra on extracellular media of diffuse large B-cell lymphoma (DLBCL) cells. Despite a lower sensitivity, the resolution of pure shift data was found to be better than that of the standard approach, which provides a more detailed vision of the exo-metabolome. The statistical analyses carried out on the resulting metabolic profiles allow us to successfully highlight several metabolic pathways affected by these drugs. Notably, we show that Kidrolase plays a major role in the metabolic pathways of this DLBCL cell line.

Molecular networking in infectious disease models

Small molecule metabolites are the product of many enzymatic reactions. Metabolomics thus opens a window into enzyme activity and function, integrating effects at the post-translational, proteome, transcriptome and genome level. In addition, small molecules can themselves regulate enzyme activity, expression and function both via substrate availability mechanisms and through allosteric regulation. Metabolites are therefore at the nexus of infectious diseases, regulating nutrient availability to the pathogen, immune responses, tropism, and host disease tolerance and resilience. Analysis of metabolomics data is however complex, particularly in terms of metabolite annotation. An emerging valuable approach to extend metabolite annotations beyond existing compound libraries and to identify infection-induced chemical changes is molecular networking. In this chapter, we discuss the applications of molecular networking in the context of infectious diseases specifically, with a focus on considerations relevant to these biological systems.

Strategies for structure elucidation of small molecules based on LC–MS/MS data from complex biological samples

LC–MS/MS is a major analytical platform for metabolomics, which has become a recent hotspot in the research fields of life and environmental sciences. By contrast, structure elucidation of small molecules based on LC–MS/MS data remains a major challenge in the chemical and biological interpretation of untargeted metabolomics datasets. In recent years, several strategies for structure elucidation using LC–MS/MS data from complex biological samples have been proposed, these strategies can be simply categorized into two types, one based on structure annotation of mass spectra and for the other on retention time prediction. These strategies have helped many scientists conduct research in metabolite-related fields and are indispensable for the development of future tools. Here, we summarized the characteristics of the current tools and strategies for structure elucidation of small molecules based on LC–MS/MS data, and further discussed the directions and perspectives to improve the power of the tools or strategies for structure elucidation.

Metabolomics Insights into Chemical Convergence in Xanthomonas perforans and Metabolic Changes Following Treatment with the Small Molecule Carvacrol

Microbes are natural chemical factories and their metabolome comprise diverse arrays of chemicals. The genus Xanthomonas comprises some of the most important plant pathogens causing devastating yield losses globally and previous studies suggested that species in the genus are untapped chemical minefields. In this study, we applied an untargeted metabolomics approach to study the metabolome of a globally spread important xanthomonad, X. perforans. The pathogen is difficult to manage, but recent studies suggest that the small molecule carvacrol was efficient in disease control. Bacterial strains were treated with carvacrol, and samples were taken at time intervals (1 and 6 h). An untreated control was also included. There were five replicates for each sample and samples were prepared for metabolomics profiling using the standard procedure. Metabolomics profiling was carried out using a thermo Q-Exactive orbitrap mass spectrometer with Dionex ultra high-performance liquid chromatography (UHPLC) and an autosampler. Annotation of significant metabolites using the Metabolomics Standards Initiative level 2 identified an array of novel metabolites that were previously not reported in Xanthomonas perforans. These metabolites include methoxybrassinin and cyclobrassinone, which are known metabolites of brassicas; sarmentosin, a metabolite of the Passiflora-heliconiine butterfly system; and monatin, a naturally occurring sweetener found in Sclerochiton ilicifolius. To our knowledge, this is the first report of these metabolites in a microbial system. Other significant metabolites previously identified in non-Xanthomonas systems but reported in this study include maculosin; piperidine; β-carboline alkaloids, such as harman and derivatives; and several important medically relevant metabolites, such as valsartan, metharbital, pirbuterol, and ozagrel. This finding is consistent with convergent evolution found in reported biological systems. Analyses of the effect of carvacrol in time-series and associated pathways suggest that carvacrol has a global effect on the metabolome of X. perforans, showing marked changes in metabolites that are critical in energy biosynthesis and degradation pathways, amino acid pathways, nucleic acid pathways, as well as the newly identified metabolites whose pathways are unknown. This study provides the first insight into the X. perforans metabolome and additionally lays a metabolomics-guided foundation for characterization of novel metabolites and pathways in xanthomonad systems.

Modifying Chromatography Conditions for Improved Unknown Feature Identification in Untargeted Metabolomics

Untargeted metabolomics is an essential component of systems biology research, but it is plagued by a high proportion of detectable features not identified with a chemical structure. Liquid chromatography-tandem mass spectrometry (LC-MS/MS) experiments produce spectra that can be searched against databases to help identify or classify these unknowns, but many features do not generate spectra of sufficient quality to enable successful annotation. Here, we explore alterations to gradient length, mass loading, and rolling precursor ion exclusion parameters for reversed phase liquid chromatography (RPLC) and hydrophilic interaction liquid chromatography (HILIC) that improve compound identification performance for human plasma samples. A manual review of spectral matches from the HILIC data set was used to determine reasonable thresholds for search score and other metrics to enable semi-automated MS/MS data analysis. Compared to typical LC-MS/MS conditions, methods adapted for compound identification increased the total number of unique metabolites that could be matched to a spectral database from 214 to 2052. Following data alignment, 68.0% of newly identified features from the modified conditions could be detected and quantitated using a routine 20-min LC-MS run. Finally, a localized machine learning model was developed to classify the remaining unknowns and select a subset that shared spectral characteristics with successfully identified features. A total of 576 and 749 unidentified features in the HILIC and RPLC data sets were classified by the model as high-priority unknowns or higher-importance targets for follow-up analysis. Overall, our study presents a simple strategy to more deeply annotate untargeted metabolomics data for a modest additional investment of time and sample.

Resistance to a nonselective 4-hydroxyphenylpyruvate dioxygenase-inhibiting herbicide via novel reduction-dehydration-glutathione conjugation in Amaranthus tuberculatus

Metabolic resistance to 4-hydroxyphenylpyruvate dioxygenase (HPPD)-inhibiting herbicides is a threat in controlling waterhemp (Amaranthus tuberculatus) in the USA. We investigated resistance mechanisms to syncarpic acid-3 (SA3), a nonselective, noncommercial HPPD-inhibiting herbicide metabolically robust to Phase I oxidation, in multiple-herbicide-resistant (MHR) waterhemp populations (SIR and NEB) and HPPD inhibitor-sensitive populations (ACR and SEN). Dose-response experiments with SA3 provided ED₅₀ -based resistant : sensitive ratios of at least 18-fold. Metabolism experiments quantifying parent SA3 remaining in excised leaves during a time course indicated MHR populations displayed faster rates of SA3 metabolism compared to HPPD inhibitor-sensitive populations. SA3 metabolites were identified via LC-MS-based untargeted metabolomics in whole plants. A Phase I metabolite, likely generated by cytochrome P450-mediated alkyl hydroxylation, was detected but was not associated with resistance. A Phase I metabolite consistent with ketone reduction followed by water elimination was detected, creating a putative α,β-unsaturated carbonyl resembling a Michael acceptor site. A Phase II glutathione-SA3 conjugate was associated with resistance. Our results revealed a novel reduction-dehydration-GSH conjugation detoxification mechanism. SA3 metabolism in MHR waterhemp is thus atypical compared to commercial HPPD-inhibiting herbicides. This previously uncharacterized detoxification mechanism presents a unique opportunity for future biorational design by blocking known sites of herbicide metabolism in weeds.

Applications of quantitative metabolomics to revolutionize early diagnosis of inborn errors of metabolism in India

Inborn errors of metabolism (IEMs) are a group of disorders caused by disruption of metabolic pathways, which leads to accumulation, decreased circulating levels, or increased excretion of metabolites as a consequence of the underlying genetic defects. These heterogeneous groups of disorders cause significant neonatal and infant mortality across the whole world and it is of utmost concern for developing countries like India owing to lack of awareness and standard preventive strategies like newborn screening (NBS). Though the predictive cumulative incidence of IEMs is said to be ∼1:800 newborns, data pertaining to the true prevalence of individual IEMs is not available in the context of Indian population. There is a need for a large population-based study to get a clear picture of the prevalence of different IEMs. One of the best ways to screen for IEMs is by applying advanced liquid chromatography-mass spectrometry (LC-MS) technology using a quantitative metabolomics approaches such as selected or multiple reaction monitoring (SRM or MRM). Recent developments in LC-MS/MRM based quantification of marker metabolites in newborns have opened a novel opportunity to screen multiple disorders simultaneously from a minuscule volume of biological fluids. In this review article, we have highlighted how LC-MS/MRM based metabolomics approach with its high sensitivity and diagnostic capability can make an impact on the nation's public health through NBS programs.

MassGenie: A Transformer-Based Deep Learning Method for Identifying Small Molecules from Their Mass Spectra

The 'inverse problem' of mass spectrometric molecular identification ('given a mass spectrum, calculate/predict the 2D structure of the molecule whence it came') is largely unsolved, and is especially acute in metabolomics where many small molecules remain unidentified. This is largely because the number of experimentally available electrospray mass spectra of small molecules is quite limited. However, the forward problem ('calculate a small molecule's likely fragmentation and hence at least some of its mass spectrum from its structure alone') is much more tractable, because the strengths of different chemical bonds are roughly known. This kind of molecular identification problem may be cast as a language translation problem in which the source language is a list of high-resolution mass spectral peaks and the 'translation' a representation (for instance in SMILES) of the molecule. It is thus suitable for attack using the deep neural networks known as transformers. We here present MassGenie, a method that uses a transformer-based deep neural network, trained on ~6 million chemical structures with augmented SMILES encoding and their paired molecular fragments as generated in silico, explicitly including the protonated molecular ion. This architecture (containing some 400 million elements) is used to predict the structure of a molecule from the various fragments that may be expected to be observed when some of its bonds are broken. Despite being given essentially no detailed nor explicit rules about molecular fragmentation methods, isotope patterns, rearrangements, neutral losses, and the like, MassGenie learns the effective properties of the mass spectral fragment and valency space, and can generate candidate molecular structures that are very close or identical to those of the 'true' molecules. We also use VAE-Sim, a previously published variational autoencoder, to generate candidate molecules that are 'similar' to the top hit. In addition to using the 'top hits' directly, we can produce a rank order of these by 'round-tripping' candidate molecules and comparing them with the true molecules, where known. As a proof of principle, we confine ourselves to positive electrospray mass spectra from molecules with a molecular mass of 500Da or lower, including those in the last CASMI challenge (for which the results are known), getting 49/93 (53%) precisely correct. The transformer method, applied here for the first time to mass spectral interpretation, works extremely effectively both for mass spectra generated in silico and on experimentally obtained mass spectra from pure compounds. It seems to act as a Las Vegas algorithm, in that it either gives the correct answer or simply states that it cannot find one. The ability to create and to 'learn' millions of fragmentation patterns in silico, and therefrom generate candidate structures (that do not have to be in existing libraries) directly, thus opens up entirely the field of de novo small molecule structure prediction from experimental mass spectra.

Metabolite discovery through global annotation of untargeted metabolomics data

Liquid chromatography-high-resolution mass spectrometry (LC-MS)-based metabolomics aims to identify and quantify all metabolites, but most LC-MS peaks remain unidentified. Here we present a global network optimization approach, NetID, to annotate untargeted LC-MS metabolomics data. The approach aims to generate, for all experimentally observed ion peaks, annotations that match the measured masses, retention times and (when available) tandem mass spectrometry fragmentation patterns. Peaks are connected based on mass differences reflecting adduction, fragmentation, isotopes, or feasible biochemical transformations. Global optimization generates a single network linking most observed ion peaks, enhances peak assignment accuracy, and produces chemically informative peak-peak relationships, including for peaks lacking tandem mass spectrometry spectra. Applying this approach to yeast and mouse data, we identified five previously unrecognized metabolites (thiamine derivatives and N-glucosyl-taurine). Isotope tracer studies indicate active flux through these metabolites. Thus, NetID applies existing metabolomic knowledge and global optimization to substantially improve annotation coverage and accuracy in untargeted metabolomics datasets, facilitating metabolite discovery.

HERMES: a molecular-formula-oriented method to target the metabolome

Comprehensive metabolome analyses are essential for biomedical, environmental, and biotechnological research. However, current MS1- and MS2-based acquisition and data analysis strategies in untargeted metabolomics result in low identification rates of metabolites. Here we present HERMES, a molecular-formula-oriented and peak-detection-free method that uses raw LC/MS1 information to optimize MS2 acquisition. Investigating environmental water, Escherichia coli, and human plasma extracts with HERMES, we achieved an increased biological specificity of MS2 scans, leading to improved mass spectral similarity scoring and identification rates when compared with a state-of-the-art data-dependent acquisition (DDA) approach. Thus, HERMES improves sensitivity, selectivity, and annotation of metabolites. HERMES is available as an R package with a user-friendly graphical interface for data analysis and visualization.

Liquid Chromatography-Mass Spectrometry for Clinical Metabolomics: An Overview

Metabolomics is a discipline that offers a comprehensive analysis of metabolites in biological samples. In the last decades, the notable evolution in liquid chromatography and mass spectrometry technologies has driven an exponential progress in LC-MS-based metabolomics. Targeted and untargeted metabolomics strategies are important tools in health and medical science, especially in the study of disease-related biomarkers, drug discovery and development, toxicology, diet, physical exercise, and precision medicine. Clinical and biological problems can now be understood in terms of metabolic phenotyping. This overview highlights the current approaches to LC-MS-based metabolomics analysis and its applications in the clinical research.

A guide to interrogating immunometabolism

The metabolic charts memorized in early biochemistry courses, and then later forgotten, have come back to haunt many immunologists with new recognition of the importance of these pathways. Metabolites and the activity of metabolic pathways drive energy production, macromolecule synthesis, intracellular signalling, post-translational modifications and cell survival. Immunologists who identify a metabolic phenotype in their system are often left wondering where to begin and what does it mean? Here, we provide a framework for navigating and selecting the appropriate biochemical techniques to explore immunometabolism. We offer recommendations for initial approaches to develop and test metabolic hypotheses and how to avoid common mistakes. We then discuss how to take things to the next level with metabolomic approaches, such as isotope tracing and genetic approaches. By proposing strategies and evaluating the strengths and weaknesses of different methodologies, we aim to provide insight, note important considerations and discuss ways to avoid common misconceptions. Furthermore, we highlight recent studies demonstrating the power of these metabolic approaches to uncover the role of metabolism in immunology. By following the framework in this Review, neophytes and seasoned investigators alike can venture into the emerging realm of cellular metabolism and immunity with confidence and rigour.

Recognizing Contamination Fragment Ions in Liquid Chromatography-Tandem Mass Spectrometry Data

Tandem mass spectral (MS/MS) data in liquid chromatography-tandem mass spectrometry (LC-MS/MS) analysis are often contaminated as the selection of precursor ions is based on a low-resolution quadrupole mass filter. In this work, we developed a strategy to differentiate contamination fragment ions (CFIs) from true fragment ions (TFIs) in an MS/MS spectrum. The rationale is that TFIs should coelute with their parent ions, but CFIs should not. To assess coelution, we performed a parallel LC-MS/MS analysis in data-independent acquisition (DIA) with all-ion-fragmentation (AIF) mode. Using the DIA (AIF) data, peak-peak correlation (PPC) score is calculated between the extracted ion chromatogram (EIC) of the fragment ion using the MS/MS scans and the EIC of the precursor ion using the MS1 scans. A high PPC score is an indication of TFIs, and a low PPC score is an indication of CFIs. Tested using metabolomics data generated by high resolution QTOF and Orbitrap MS from various vendors in different LC-MS configurations, we found that more than 70% of the fragment ions have PPC scores < 0.8 and identified three common sources of CFIs, including (1) solvent contamination, (2) adjacent chemical contamination, and (3) undetermined signals from artifacts and noise. Combining PPC scores with other precursor and fragment ion information, we further developed a machine learning model that can robustly and conservatively predict CFIs. Incorporating the machine learning model, we created an R program, MS2Purifier, to automatically recognize CFIs and clean MS/MS spectra of metabolic features in LC-MS/MS data with high sensitivity and specificity.

A Workflow to Perform Targeted Metabolomics at the Untargeted Scale on a Triple Quadrupole Mass Spectrometer

The thousands of features commonly observed when performing untargeted metabolomics with quadrupole time-of-flight (QTOF) and Orbitrap mass spectrometers often correspond to only a few hundred unique metabolites of biological origin, which is in the range of what can be assayed in a single targeted metabolomics experiment by using a triple quadrupole (QqQ) mass spectrometer. A major benefit of performing targeted metabolomics with QqQ mass spectrometry is the affordability of the instruments relative to high-resolution QTOF and Orbitrap platforms. Optimizing targeted methods to profile hundreds of metabolites on a QqQ mass spectrometer, however, has historically been limited by the availability of authentic standards, particularly for "unknowns" that have yet to be structurally identified. Here, we report a strategy to develop multiple reaction monitoring (MRM) methods for QqQ instruments on the basis of high-resolution spectra, thereby enabling us to use data from untargeted metabolomics to design targeted experiments without the need for authentic standards. We demonstrate that using high-resolution fragmentation data alone to design MRM methods results in the same quantitative performance as when methods are optimized by measuring authentic standards on QqQ instruments, as is conventionally done. The approach was validated by showing that Orbitrap ID-X data can be used to establish MRM methods on a Thermo TSQ Altis and two Agilent QqQs for hundreds of metabolites, including unknowns, without a dependence on standards. Finally, we highlight an application where metabolite profiling was performed on an ID-X and a QqQ by using the strategy introduced here, with both data sets yielding the same result. The described approach therefore allows us to use QqQ instruments, which are often associated with targeted metabolomics, to profile knowns and unknowns at a comprehensive scale that is typical of untargeted metabolomics.

Unraveling mosquito metabolism with mass spectrometry-based metabolomics

Nearly half a million people die annually due to mosquito-borne diseases. Despite aggressive mosquito population-control efforts, current strategies are limited in their ability to control these vectors. A better understanding of mosquito metabolism through modern approaches can contribute to the discovery of novel metabolic targets and/or regulators and lead to the development of better mosquito-control strategies. Currently, cutting-edge technologies such as gas or liquid chromatography-mass spectrometry-based metabolomics are considered 'mature technologies' in many life-science disciplines but are still an emerging area of research in medical entomology. This review primarily discusses recent developments and progress in the application of mass spectrometry-based metabolomics to answer multiple biological questions related to mosquito metabolism.

Simultaneous 3-Nitrophenylhydrazine Derivatization Strategy of Carbonyl, Carboxyl and Phosphoryl Submetabolome for LC-MS/MS-Based Targeted Metabolomics with Improved Sensitivity and Coverage

Metabolomics is a powerful and essential technology for profiling metabolic phenotypes and exploring metabolic reprogramming, which enables the identification of biomarkers and provides mechanistic insights into physiology and disease. However, its applications are still limited by the technical challenges particularly in its detection sensitivity for the analysis of biological samples with limited amount, necessitating the development of highly sensitive approaches. Here, we developed a highly sensitive liquid chromatography tandem mass spectrometry method based on a 3-nitrophenylhydrazine (3-NPH) derivatization strategy that simultaneously targets carbonyl, carboxyl, and phosphoryl groups for targeted metabolomic analysis (HSD^ccp-TM) in biological samples. By testing 130 endogenous metabolites including organic acids, amino acids, carbohydrates, nucleotides, carnitines, and vitamins, we showed that the derivatization strategy resulted in significantly improved detection sensitivity and chromatographic separation capability. Metabolic profiling of merely 60 oocytes and 5000 hematopoietic stem cells primarily isolated from mice demonstrated that this method enabled routine metabolomic analysis in trace amounts of biospecimens. Moreover, the derivatization strategy bypassed the tediousness of inferring the MS fragmentation patterns and simplified the complexity of monitoring ion pairs of metabolites, which greatly facilitated the metabolic flux analysis (MFA) for glycolysis, the tricarboxylic acid (TCA) cycle, and pentose phosphate pathway (PPP) in cell cultures. In summary, the novel 3-NPH derivatization-based method with high sensitivity, good chromatographic separation, and broad coverage showed great potential in promoting metabolomics and MFA, especially in trace amounts of biospecimens.

ISFrag: De Novo Recognition of In-Source Fragments for Liquid Chromatography-Mass Spectrometry Data

In-source fragmentation (ISF) is a naturally occurring phenomenon during electrospray ionization (ESI) in liquid chromatography-mass spectrometry (LC-MS) analysis. ISF leads to false metabolite annotation in untargeted metabolomics, prompting misinterpretation of the underlying biological mechanisms. Conventional metabolomic data cleaning mainly focuses on the annotation of adducts and isotopes, and the recognition of ISF features is mainly based on common neutral losses and the LC coelution pattern. In this work, we recognized three increasingly important patterns of ISF features, including (1) coeluting with their precursor ions, (2) being in the tandem MS (MS²) spectra of their precursor ions, and (3) sharing similar MS² fragmentation patterns with their precursor ions. Based on these patterns, we developed an R package, ISFrag, to comprehensively recognize all possible ISF features from LC-MS data generated from full-scan, data-dependent acquisition, and data-independent acquisition modes without the assistance of common neutral loss information or MS² spectral library. Tested using metabolite standards, we achieved a 100% correct recognition of level 1 ISF features and over 80% correct recognition for level 2 ISF features. Further application of ISFrag on untargeted metabolomics data allows us to identify ISF features that can potentially cause false metabolite annotation at an omics-scale. With the help of ISFrag, we performed a systematic investigation of how ISF features are influenced by different MS parameters, including capillary voltage, end plate offset, ion energy, and "collision energy". Our results show that while increasing energies can increase the number of real metabolic features and ISF features, the percentage of ISF features might not necessarily increase. Finally, using ISFrag, we created an ISF pathway to visualize the relationships between multiple ISF features that belong to the same precursor ion. ISFrag is freely available on GitHub (https://github.com/HuanLab/ISFrag).

Untargeted metabolomics approach to discriminate mistletoe commercial products

Mistletoe (Viscum album L.) is used in German-speaking European countries in the field of integrative oncology linking conventional and complementary medicine therapies to improve quality of life. Various companies sell extracts, fermented or not, for injection by subcutaneous or intra-tumoral route with a regulatory status of anthroposophic medicinal products (European Medicinal Agency (EMA) assessment status). These companies as well as anthroposophical physicians argue that complex matrices composed of many molecules in mixture are necessary for activity and that the host tree of the mistletoe parasitic plant is the main determining factor for this matrix composition. The critical point is that parenteral devices of European mistletoe extracts do not have a standard chemical composition regulated by EMA quality guidelines, because they are not drugs, regulatory speaking. However, the mechanism of mistletoe's anticancer activity and its effectiveness in treating and supporting cancer patients are not fully understood. Because of this lack of transparency and knowledge regarding the matrix chemical composition, we undertook an untargeted metabolomics study of several mistletoe extracts to explore and compare their fingerprints by LC-(HR)MS(/MS) and 1H-NMR. Unexpectedly, we showed that the composition was primarily driven by the manufacturer/preparation method rather than the different host trees. This differential composition may cause differences in immunostimulating and anti-cancer activities of the different commercially available mistletoe extracts as illustrated by structure-activity relationships based on LC-MS/MS and 1H-NMR identifications completed by docking experiments. In conclusion, in order to move towards an evidence-based medicine use of mistletoe, it is a priority to bring rigor and quality, chemically speaking.

Ultra-high-performance liquid chromatography high-resolution mass spectrometry variants for metabolomics research

Ultra-high-performance liquid chromatography high-resolution mass spectrometry (UHPLC-HRMS) variants currently represent the best tools to tackle the challenges of complexity and lack of comprehensive coverage of the metabolome. UHPLC offers flexible and efficient separation coupled with high-sensitivity detection via HRMS, allowing for the detection and identification of a broad range of metabolites. Here we discuss current common strategies for UHPLC-HRMS-based metabolomics, with a focus on expanding metabolome coverage.

Marine Actinomycetes, New Sources of Biotechnological Products

The Actinomycetales order is one of great genetic and functional diversity, including diversity in the production of secondary metabolites which have uses in medical, environmental rehabilitation, and industrial applications. Secondary metabolites produced by actinomycete species are an abundant source of antibiotics, antitumor agents, anthelmintics, and antifungals. These actinomycete-derived medicines are in circulation as current treatments, but actinomycetes are also being explored as potential sources of new compounds to combat multidrug resistance in pathogenic bacteria. Actinomycetes as a potential to solve environmental concerns is another area of recent investigation, particularly their utility in the bioremediation of pesticides, toxic metals, radioactive wastes, and biofouling. Other applications include biofuels, detergents, and food preservatives/additives. Exploring other unique properties of actinomycetes will allow for a deeper understanding of this interesting taxonomic group. Combined with genetic engineering, microbial experimental evolution, and other enhancement techniques, it is reasonable to assume that the use of marine actinomycetes will continue to increase. Novel products will begin to be developed for diverse applied research purposes, including zymology and enology. This paper outlines the current knowledge of actinomycete usage in applied research, focusing on marine isolates and providing direction for future research.

Secretory products of the corpus luteum and preeclampsia

BACKGROUND

Despite significant advances in our understanding of the pathophysiology of preeclampsia (PE), there are still many unknowns and controversies in the field. Women undergoing frozen-thawed embryo transfer (FET) to a hormonally prepared endometrium have been found to have an unexpected increased risk of PE compared to women who receive embryos in a natural FET cycle. The differences in risk have been hypothesized to be related to the absence or presence of a functioning corpus luteum (CL).

OBJECTIVE AND RATIONALE

To evaluate the literature on secretory products of the CL that could be essential for a healthy pregnancy and could reduce the risk of PE in the setting of FET.

SEARCH METHODS

For this review, pertinent studies were searched in PubMed/Medline (updated June 2020) using common keywords applied in the field of assisted reproductive technologies, CL physiology and preeclampsia. We also screened the complete list of references in recent publications in English (both animal and human studies) on the topics investigated. Given the design of this work as a narrative review, no formal criteria for study selection or appraisal were utilized.

OUTCOMES

The CL is a major source of multiple factors regulating reproduction. Progesterone, estradiol, relaxin and vasoactive and angiogenic substances produced by the CL have important roles in regulating its functional lifespan and are also secreted into the circulation to act remotely during early stages of pregnancy. Beyond the known actions of progesterone and estradiol on the uterus in early pregnancy, their metabolites have angiogenic properties that may optimize implantation and placentation. Serum levels of relaxin are almost undetectable in pregnant women without a CL, which precludes some maternal cardiovascular and renal adaptations to early pregnancy. We suggest that an imbalance in steroid hormones and their metabolites and polypeptides influencing early physiologic processes such as decidualization, implantation, angiogenesis and maternal haemodynamics could contribute to the increased PE risk among women undergoing programmed FET cycles.

WIDER IMPLICATIONS

A better understanding of the critical roles of the secretory products of the CL during early pregnancy holds the promise of improving the efficacy and safety of ART based on programmed FET cycles.

Development of analytical protocol for the investigation of transformation products of pre-harvest fungicides in fruits using LC-MS/MS methods

Chemical protection of plants is critical to permit all-year-long availability of plant products. These chemical agents are transformed both biotically and abiotically after spraying. Our purpose was to develop a workflow which is suitable to investigate those transformation products from plant matrices. Two field trials were set up in two years in two different plant matrices (apples, cherries) to develop a workflow to map the transformation products (TP) of three selected pre-harvest fungicides (boscalid, fluopyram, pyraclostrobin) in the fruits. Modified QuEChERS extraction method was applied for the extraction of TPs from the fruit matrices. We used liquid chromatograph-mass spectrometers to identify and confirm the transformation products of fungicides. LC-QTOF-MS method was suitable to map the key fragmentation routes of parent fungicides. Based on fragmentation pathways, MRM (multiple-reaction monitoring) transitions of fungicide-metabolites mentioned in the literature were predicted. HPLC-QTRAP-MS in target mode was successfully applied to monitor trace level of metabolites and in some cases their isomers. For confirmation of the identified metabolites, LC-QTOF-MS was used. Five earlier-documented as well as one novel transformation product (deschloro-FLP) were found in the investigated fruit samples, the latter has not been reported in plant matrices so far. Area-normalisation method was used to follow the relative concentration of the transformation products in the fruits as a function of time.

SteroidXtract: Deep Learning-Based Pattern Recognition Enables Comprehensive and Rapid Extraction of Steroid-Like Metabolic Features for Automated Biology-Driven Metabolomics

Despite the vast amount of metabolic information that can be captured in untargeted metabolomics, many biological applications are looking for a biology-driven metabolomics platform that targets a set of metabolites that are relevant to the given biological question. Steroids are a class of important molecules that play critical roles in many physiological systems and diseases. Besides known steroids, there are a large number of unknown steroids that have not been reported in the literature. The ability to rapidly detect and quantify both known and unknown steroid molecules in a biological sample can greatly accelerate a broad range of steroid-focused life science research. This work describes the development and application of SteroidXtract, a convolutional neural network (CNN)-based bioinformatics tool that can recognize steroid molecules in mass spectrometry (MS)-based untargeted metabolomics using their unique tandem MS (MS²) spectral patterns. SteroidXtract was trained using a comprehensive set of standard MS² spectra from MassBank of North America (MoNA) and an in-house steroid library. Data augmentation strategies, including intensity thresholding and Gaussian noise addition, were created and applied to minimize data overfitting caused by the limited number of standard steroid MS² spectra. The CNN model embedded in SteroidXtract was further compared with random forest and XGBoost using nested cross-validations to demonstrate its performance. Finally, SteroidXtract was applied in several metabolomics studies to demonstrate its sensitivity, specificity, and robustness. Compared to conventional statistics-driven metabolomics data interpretation, our work offers a novel automated biology-driven approach to interpreting untargeted metabolomics data, prioritizing biologically important molecules with high throughput and sensitivity.

metabCombiner: Paired Untargeted LC-HRMS Metabolomics Feature Matching and Concatenation of Disparately Acquired Data Sets

LC-HRMS experiments detect thousands of compounds, with only a small fraction of them identified in most studies. Traditional data processing pipelines contain an alignment step to assemble the measurements of overlapping features across samples into a unified table. However, data sets acquired under nonidentical conditions are not amenable to this process, mostly due to significant alterations in chromatographic retention times. Alignment of features between disparately acquired LC-MS metabolomics data could aid collaborative compound identification efforts and enable meta-analyses of expanded data sets. Here, we describe metabCombiner, a new computational pipeline for matching known and unknown features in a pair of untargeted LC-MS data sets and concatenating their abundances into a combined table of intersecting feature measurements. metabCombiner groups features by mass-to-charge (m/z) values to generate a search space of possible feature pair alignments, fits a spline through a set of selected retention time ordered pairs, and ranks alignments by m/z, mapped retention time, and relative abundance similarity. We evaluated this workflow on a pair of plasma metabolomics data sets acquired with different gradient elution methods, achieving a mean absolute retention time prediction error of roughly 0.06 min and a weighted per-compound matching accuracy of approximately 90%. We further demonstrate the utility of this method by comprehensively mapping features in urine and muscle metabolomics data sets acquired from different laboratories. metabCombiner has the potential to bridge the gap between otherwise incompatible metabolomics data sets and is available as an R package at https://github.com/hhabra/metabCombiner and Bioconductor.

Interlaboratory Comparison of Untargeted Mass Spectrometry Data Uncovers Underlying Causes for Variability

Despite the value of mass spectrometry in modern natural products discovery workflows, it remains very difficult to compare data sets between laboratories. In this study we compared mass spectrometry data for the same sample set from two different laboratories (quadrupole time-of-flight and quadrupole-Orbitrap) and evaluated the similarity between these two data sets in terms of both mass spectrometry features and their ability to describe the chemical composition of the sample set. Somewhat surprisingly, the two data sets, collected with appropriate controls and replication, had very low feature overlap (25.7% of Laboratory A features overlapping 21.8% of Laboratory B features). Our data clearly demonstrate that differences in fragmentation, charge state, and adduct formation in the ionization source are a major underlying cause for these differences. Consistent with other recent literature, these findings challenge the conventional wisdom that electrospray ionization mass spectrometry (ESI-MS) yields a simple one-to-one correspondence between analytes in solution and features in the data set. Importantly, despite low overlap in feature lists, principal component analysis (PCA) generated qualitatively similar PCA plots. Overall, our findings demonstrate that comparing untargeted metabolomics data between laboratories is challenging, but that data sets with low feature overlap can yield the same qualitative description of a sample set using PCA.

A map of mass spectrometry-based in silico fragmentation prediction and compound identification in metabolomics

Metabolomics, the comprehensive study of the metabolome, and lipidomics-the large-scale study of pathways and networks of cellular lipids-are major driving forces in enabling personalized medicine. Complicated and error-prone data analysis still remains a bottleneck, however, especially for identifying novel metabolites. Comparing experimental mass spectra to curated databases containing reference spectra has been the gold standard for identification of compounds, but constructing such databases is a costly and time-demanding task. Many software applications try to circumvent this process by utilizing cutting-edge advances in computational methods-including quantum chemistry and machine learning-and simulate mass spectra by performing theoretical, so called in silico fragmentations of compounds. Other solutions concentrate directly on experimental spectra and try to identify structural properties by investigating reoccurring patterns and the relationships between them. The considerable progress made in the field allows recent approaches to provide valuable clues to expedite annotation of experimental mass spectra. This review sheds light on individual strengths and weaknesses of these tools, and attempts to evaluate them-especially in view of lipidomics, when considering complex mixtures found in biological samples as well as mass spectrometer inter-instrument variability.

Approaches in metabolomics for regulatory toxicology applications

Innovative methodological approaches are needed to conduct human health and environmental risk assessments on a growing number of marketed chemicals. Metabolomics is progressively proving its value as an efficient strategy to perform toxicological evaluations of new and existing substances, and it will likely become a key tool to accelerate chemical risk assessments. However, additional guidance with widely accepted and harmonized procedures is needed before metabolomics can be routinely incorporated in decision-making for regulatory purposes. The aim of this review is to provide an overview of metabolomic strategies that have been successfully employed in toxicity assessment as well as the most promising workflows in a regulatory context. First, we provide a general view of the different steps of regulatory toxicology-oriented metabolomics. Emphasis is put on three key elements: robustness of experimental design, choice of analytical platform, and use of adapted data treatment tools. Then, examples in which metabolomics supported regulatory toxicology outputs in different scenarios are reviewed, including chemical grouping, elucidation of mechanisms of toxicity, and determination of points of departure. The overall intention is to provide insights into why and how to plan and conduct metabolomic studies for regulatory toxicology purposes.

Targeting unique biological signals on the fly to improve MS/MS coverage and identification efficiency in metabolomics

When using liquid chromatography/mass spectrometry (LC/MS) to perform untargeted metabolomics, it is common to detect thousands of features from a biological extract. Although it is impractical to collect non-chimeric MS/MS data for each in a single chromatographic run, this is generally unnecessary because most features do not correspond to unique metabolites of biological relevance. Here we show that relatively simple data-processing strategies that can be applied on the fly during acquisition of data with an Orbitrap ID-X, such as blank subtraction and well-established adduct or isotope calculations, decrease the number of features to target for MS/MS analysis by up to an order of magnitude for various types of biological matrices. We demonstrate that annotating these non-biological contaminants and redundancies in real time during data acquisition enables comprehensive MS/MS data to be acquired on each remaining feature at a single collision energy. To ensure that an appropriate collision energy is applied, we introduce a method using a series of hidden ion-trap scans in an Orbitrap ID-X to find an optimal value for each feature that can then be applied in a subsequent high-resolution Orbitrap scan. Data from 100 metabolite standards indicate that this real-time optimization of collision energies leads to more informative MS/MS patterns compared to using a single fixed collision energy alone. As a benchmark to evaluate the overall workflow, we manually annotated unique biological features by independently subjecting E. coli samples to a credentialing analysis. While credentialing led to a more rigorous reduction in feature number, on-the-fly annotation with blank subtraction on an Orbitrap ID-X did not inappropriately discard unique biological metabolites. Taken together, our results reveal that optimal fragmentation data can be obtained in a single LC/MS/MS run for >90% of the unique biological metabolites in a sample when features are annotated during acquisition and collision energies are selected by using parallel mass spectrometry detection.

Metabolomic Analysis of Aqueous Humor Identifies Aberrant Amino Acid and Fatty Acid Metabolism in Vogt-Koyanagi-Harada and Behcet's Disease

To investigate aqueous metabolic profiles in Vogt-Koyanagi-Harada (VKH) and Behcet's disease (BD), we applied ultra-high-performance liquid chromatography equipped with quadrupole time-of flight mass spectrometry in aqueous humor samples collected from these patients and controls. Metabolite levels in these three groups were analyzed by univariate logistic regression. The differential metabolites were subjected to subsequent pathway analysis by MetaboAnalyst. The results showed that both partial-least squares discrimination analysis and hierarchical clustering analysis showed specific aqueous metabolite profiles when comparing VKH, BD, and controls. There were 28 differential metabolites in VKH compared to controls and 29 differential metabolites in BD compared to controls. Amino acids and fatty acids were the two most abundant categories of differential metabolites. Furthermore, pathway enrichment analysis identified several perturbed pathways, including pantothenate and CoA biosynthesis when comparing VKH with the control group, and D-arginine and D-ornithine metabolism and phenylalanine metabolism when comparing BD with the control group. Aminoacyl-tRNA biosynthesis was altered in both VKH and BD when compared to controls. Our findings suggest that amino acids metabolism as well as two fatty acids, palmitic acid and oleic acid, may be involved in the pathogenesis of BD and VKH.