Enhancing Metabolome Coverage in Data-Dependent LC–MS/MS Analysis through an Integrated Feature Extraction Strategy

Untargeted Metabolomics and Targeted Phytohormone Profiling of Sweet Aloes (Euphorbia neriifolia) from Guyana: An Assessment of Asthma Therapy Potential in Leaf Extracts and Latex

Background/Objectives:Euphorbia neriifolia is a succulent plant from the therapeutically rich family of Euphorbia comprising 2000 species globally. E. neriifolia is used in Indigenous Guyanese asthma therapy. Methods: To investigate E. neriifolia's therapeutic potential, traditionally heated leaf, simple leaf, and latex extracts were evaluated for phytohormones and therapeutic compounds. Full scan, data-dependent acquisition, and parallel reaction monitoring modes via liquid chromatography Orbitrap mass spectrometry were used for screening. Results: Pathway analysis of putative features from all extracts revealed a bias towards the phenylpropanoid, terpenoid, and flavonoid biosynthetic pathways. A total of 850 compounds were annotated using various bioinformatics tools, ranging from confidence levels 1 to 3. Lipids and lipid-like molecules (34.35%), benzenoids (10.24%), organic acids and derivatives (12%), organoheterocyclic compounds (12%), and phenylpropanoids and polyketides (10.35%) dominated the contribution of compounds among the 13 superclasses. Semi-targeted screening revealed 14 out of 16 literature-relevant therapeutic metabolites detected, with greater upregulation in traditional heated extracts. Targeted screening of 39 phytohormones resulted in 25 being detected and quantified. Simple leaf extract displayed 4.4 and 45 times greater phytohormone levels than traditional heated leaf and latex extracts, respectively. Simple leaf extracts had the greatest nucleotide and riboside cytokinin and acidic phytohormone levels. In contrast, traditional heated extracts exhibited the highest free base and glucoside cytokinin levels and uniquely contained methylthiolated and aromatic cytokinins while lacking acidic phytohormones. Latex samples had trace gibberellic acid levels, the lowest free base, riboside, and nucleotide levels, with absences of aromatic, glucoside, or methylthiolated cytokinin forms. Conclusions: In addition to metabolites with possible therapeutic value for asthma treatment, we present the first look at cytokinin phytohormones in the species and Euphorbia genus alongside metabolite screening to present a comprehensive assessment of heated leaf extract used in Indigenous Guyanese asthma therapy.

Toward Automated Preprocessing of Untargeted LC-MS-Based Metabolomics Feature Lists from Human Biofluids

Maximizing the extraction of true, high-quality, nonredundant features from biofluids analyzed via LC-MS systems is challenging. Here, the R packages IPO and AutoTuner were used to optimize XCMS parameter settings for the retrieval of metabolite or lipid features in both ionization modes from either faecal or urine samples from two cohorts (n = 621). The feature lists obtained were compared with those where the parameter values were selected manually. Three categories were used to compare feature lists: 1) feature quality through removing false positives, 2) tentative metabolite identification using the Human Metabolome Database (HMDB) and 3) feature utility such as analyzing the proportion of features within intensity threshold bins. Furthermore, a PCA-based approach to feature filtering using QC samples and variable loadings was also explored under this category. Overall, more features were observed after automated selection of parameter values for all data sets (1.3- to 3.7-fold), which propagated through comparative exercises. For example, a greater number of features (on average 51 vs 45%) had a coefficient of variation (CV) < 30%. Additionally, there was a significant increase (7.6-10.4%) in the number of faecal metabolites that could be tentatively annotated, and more features were present in higher intensity threshold bins. Considering the overlap across all three categories, a greater number of features were also retained. Automated approaches that guide selection of optimal parameter values for preprocessing are important to decrease the time invested for this step, while taking advantage of the wealth of data that LC-MS systems provide.

Paper Spray Ionization Mass Spectrometry Coupled with Paper-Based Three-Dimensional Tumor Model for Rapid Metabolic Gradient Profiling

The tumor microenvironment (TME), especially with its complicated metabolic characteristics, will dynamically affect the proliferation, migration, and drug response of tumor cells. Rapid metabolic analysis brings out a deeper understanding of the TME, while the susceptibility and environmental dependence of metabolites extremely hinder real-time metabolic profiling since the TME is easily disrupted. Here, we directly integrated paper spray ionization mass spectrometry with a paper-based three-dimensional (3D) tumor model, realizing the rapid capture of metabolic gradients. The entire procedure, from sample preparation to mass spectrometry detection, took less than 4 min, which was able to provide metabolic results close to real time and contributed to understanding the real metabolic processes. At present, our method successfully detected 160 metabolites; notably, over 40 significantly gradient metabolites were revealed across the six layers of the paper-based 3D tumor model. At least 22 gradient metabolites were reported to be associated with cell viability. This strategy was powerful enough to rapidly profile metabolic gradients of a paper-based 3D tumor model for revealing cell viability changes from a metabolomics perspective.

Offline Two-Dimensional Liquid Chromatography-Mass Spectrometry for Deep Annotation of the Fecal Metabolome Following Fecal Microbiota Transplantation

Biological interpretation of untargeted LC-MS-based metabolomics data depends on accurate compound identification, but current techniques fall short of identifying most features that can be detected. The human fecal metabolome is complex, variable, incompletely annotated, and serves as an ideal matrix to evaluate novel compound identification methods. We devised an experimental strategy for compound annotation using multidimensional chromatography and semiautomated feature alignment and applied these methods to study the fecal metabolome in the context of fecal microbiota transplantation (FMT) for recurrent C. difficile infection. Pooled fecal samples were fractionated using semipreparative liquid chromatography and analyzed by an orthogonal LC-MS/MS method. The resulting spectra were searched against commercial, public, and local spectral libraries, and annotations were vetted using retention time alignment and prediction. Multidimensional chromatography yielded more than a 2-fold improvement in identified compounds compared to conventional LC-MS/MS and successfully identified several rare and previously unreported compounds, including novel fatty-acid conjugated bile acid species. Using an automated software-based feature alignment strategy, most metabolites identified by the new approach could be matched to features that were detected but not identified in single-dimensional LC-MS/MS data. Overall, our approach represents a powerful strategy to enhance compound identification and biological insight from untargeted metabolomics data.

AnnoSM: An Automated Annotation Tool for Determining the Substituent Modes on the Parent Skeleton Based on a Characteristic MS/MS Fragment Ion Library

Mass spectrometry (MS) is a powerful technology for the structural elucidation of known or unknown small molecules. However, the accuracy of MS-based structure annotation is still limited due to the presence of numerous isomers in complex matrices. There are still challenges in automatically interpreting the fine structure of molecules, such as the types and positions of substituents (substituent modes, SMs) in the structure. In this study, we employed flavones, flavonols, and isoflavones as examples to develop an automated annotation method for identifying the SMs on the parent molecular skeleton based on a characteristic MS/MS fragment ion library. Importantly, user-friendly software AnnoSM was built for the convenience of researchers with limited computational backgrounds. It achieved 76.87% top-1 accuracy on the 148 authentic standards. Among them, 22 sets of flavonoid isomers were successfully differentiated. Moreover, the developed method was successfully applied to complex matrices. One such example is the extract of Ginkgo biloba L. (EGB), in which 331 possible flavonoids with SM candidates were annotated. Among them, 23 flavonoids were verified by authentic standards. The correct SMs of 13 flavonoids were ranked first on the candidate list. In the future, this software can also be extrapolated to other classes of compounds.

HDPairFinder: A data processing platform for hydrogen/deuterium isotopic labeling-based nontargeted analysis of trace-level amino-containing chemicals in environmental water

The combination of hydrogen/deuterium (H/D) formaldehyde-based isotopic methyl labeling with solid-phase extraction and high-performance liquid chromatography-high resolution mass spectrometry (HPLC-HRMS) is a powerful analytical solution for nontargeted analysis of trace-level amino-containing chemicals in water samples. Given the huge amount of chemical information generated in HPLC-HRMS analysis, identifying all possible H/D-labeled amino chemicals presents a significant challenge in data processing. To address this, we designed a streamlined data processing pipeline that can automatically extract H/D-labeled amino chemicals from the raw HPLC-HRMS data with high accuracy and efficiency. First, we developed a cross-correlation algorithm to correct the retention time shift resulting from deuterium isotopic effects, which enables reliable pairing of H- and D-labeled peaks. Second, we implemented several bioinformatic solutions to remove false chemical features generated by in-source fragmentation, salt adduction, and natural 13C isotopes. Third, we used a data mining strategy to construct the AMINES library that consists of over 38,000 structure-disjointed primary and secondary amines to facilitate putative compound annotation. Finally, we integrated these modules into a freely available R program, HDPairFinder.R. The rationale of each module was justified and its performance tested using experimental H/D-labeled chemical standards and authentic water samples. We further demonstrated the application of HDPairFinder to effectively extract N-containing contaminants, thus enabling the monitoring of changes of primary and secondary N-compounds in authentic water samples. HDPairFinder is a reliable bioinformatic tool for rapid processing of H/D isotopic methyl labeling-based nontargeted analysis of water samples, and will facilitate a better understanding of N-containing chemical compounds in water.

EISA-EXPOSOME: One Highly Sensitive and Autonomous Exposomic Platform with Enhanced in-Source Fragmentation/Annotation

Lacking a highly sensitive exposome screening technique is one of the biggest challenges in moving exposomic research forward. Enhanced in-source fragmentation/annotation (EISA) has been developed to facilitate molecular identification in untargeted metabolomics and proteomics. In this work, with a mixture of 50 pesticides at three concentration levels (20, 4, and 0.8 ppb), we investigated the analytical performance of the EISA technique over the well-accepted targeted MS/MS mode (TMM) in the detection and identification of chemicals at low levels using a quadrupole time-of-flight (qTOF) instrument. Compared with the TMM method, the EISA technique can recognize additional 1, 20, and 23 chemicals, respectively, at the three concentration levels (20, 4, and 0.8 ppb, respectively) investigated. At the 0.8 ppb level, intensities of precursor ions and fragments observed using the EISA technique are 30-1,154 and 3-80 times higher, respectively, than those observed at the TMM mode. A higher matched fragment ratio (MFR) between the EISA technique and the TMM method was recognized for most chemicals. We further developed a chemical annotation informatics algorithm, EISA-EXPOSOME, which can automatically search each precursor ion (m/z) in the MS/MS library against the EISA MS1 spectra. This algorithm then calculated a weighted score to rank the candidate features by comparing the experimental fragment spectra to those in the library. The peak intensity, zigzag index, and retention time prediction model as well as the peak correlation coefficient were further adopted in the algorithm to filter false positives. The performance of EISA-EXPOSOME was demonstrated using a pooled dust extract with a pesticide mixture (n = 200) spiked at 5 ppb. One urine sample spiked with a contaminant mixture (n = 50) at the 5 ppb level was also used for the validation of the pipeline. Proof-of-principal application of EISA-EXPOSOME in the real sample was further evaluated on the pooled dust sample with a modified T3DB database (n = 1650). Our results show that the EISA-EXPOSOME algorithm can remarkably improve the detection and annotation coverage at trace levels beyond the traditional approach as well as facilitate the high throughput screening of suspected chemicals.

Unlocking Potentially Therapeutic Phytochemicals in Capadulla (Doliocarpus dentatus) from Guyana Using Untargeted Mass Spectrometry-Based Metabolomics

Doliocarpus dentatus is thought to have a wide variety of therapeutic phytochemicals that allegedly improve libido and cure impotence. Although a few biomarkers have been identified with potential antinociceptive and cytotoxic properties, an untargeted mass spectrometry-based metabolomics approach has never been undertaken to identify therapeutic biofingerprints for conditions, such as erectile dysfunction, in men. This study executes a preliminary phytochemical screening of the woody vine of two ecotypes of D. dentatus with renowned differences in therapeutic potential for erectile dysfunction. Liquid chromatography-mass spectrometry-based metabolomics was used to screen for flavonoids, terpenoids, and other chemical classes found to contrast between red and white ecotypes. Among the metabolite chemodiversity found in the ecotype screens, using a combination of GNPS, MS-DIAL, and SIRIUS, approximately 847 compounds were annotated at levels 2 to 4, with the majority of compounds falling under lipid and lipid-like molecules, benzenoids and phenylpropanoids, and polyketides, indicative of the contributions of the flavonoid, shikimic acid, and terpenoid biosynthesis pathways. Despite the extensive annotation, we report on 138 tentative compound identifications of potentially therapeutic compounds, with 55 selected compounds at a level-2 annotation, and 22 statistically significant therapeutic biomarkers, the majority of which were polyphenols. Epicatechin methyl gallate, catechin gallate, and proanthocyanidin A2 had the greatest significant differences and were also relatively abundant among the red and white ecotypes. These putatively identified compounds reportedly act as antioxidants, neutralizing damaging free radicals, and lowering cell oxidative stress, thus aiding in potentially preventing cellular damage and promoting overall well-being, especially for treating erectile dysfunction (ED).

Offline Two-dimensional Liquid Chromatography-Mass Spectrometry for Deep Annotation of the Fecal Metabolome

Compound identification is an essential task in the workflow of untargeted metabolomics since the interpretation of the data in a biological context depends on the correct assignment of chemical identities to the features it contains. Current techniques fall short of identifying all or even most observable features in untargeted metabolomics data, even after rigorous data cleaning approaches to remove degenerate features are applied. Hence, new strategies are required to annotate the metabolome more deeply and accurately. The human fecal metabolome, which is the focus of substantial biomedical interest, is a more complex, more variable, yet lesser-investigated sample matrix compared to widely studied sample types like human plasma. This manuscript describes a novel experimental strategy using multidimensional chromatography to facilitate compound identification in untargeted metabolomics. Pooled fecal metabolite extract samples were fractionated using offline semi-preparative liquid chromatography. The resulting fractions were analyzed by an orthogonal LC-MS/MS method, and the data were searched against commercial, public, and local spectral libraries. Multidimensional chromatography yielded more than a 3-fold improvement in identified compounds compared to the typical single-dimensional LC-MS/MS approach and successfully identified several rare and novel compounds, including atypical conjugated bile acid species. Most features identified by the new approach could be matched to features that were detectable but not identifiable in the original single-dimension LC-MS data. Overall, our approach represents a powerful strategy for deeper annotation of the metabolome that can be implemented with commercially-available instrumentation, and should apply to any dataset requiring deeper annotation of the metabolome.

A comparison of feature extraction capabilities of advanced UHPLC-HRMS data analysis tools in plant metabolomics

Data analysis of ultrahigh performance liquid chromatography-high resolution mass spectrometry (UHPLC-HRMS) is an essential and time-consuming step in plant metabolomics and feature extraction is the fundamental step for current tools. Various methods lead to different feature extraction results in practical applications, which may puzzle users for selecting adequate data analysis tools to deal with collected data. In this work, we provide a comprehensive method evaluation for some advanced UHPLC-HRMS data analysis tools in plant metabolomics, including MS-DIAL, XCMS, MZmine, AntDAS, Progenesis QI, and Compound Discoverer. Both mixtures of standards and various complex plant matrices were specifically designed for evaluating the performances of the involved method in analyzing both targeted and untargeted metabolomics. Results indicated that AntDAS provide the most acceptable feature extraction, compound identification, and quantification results in targeted compound analysis. Concerning the complex plant dataset, both MS-DIAL and AntDAS can provide more reliable results than the others. The method comparison is maybe useful for the selection of suitable data analysis tools for users.

Mechanistic Understanding of the Discrepancies between Common Peak Picking Algorithms in Liquid Chromatography–Mass Spectrometry-Based Metabolomics

Inconsistent peak picking outcomes are a critical concern in processing liquid chromatography-mass spectrometry (LC-MS)-based untargeted metabolomics data. This work systematically studied the mechanisms behind the discrepancies among five commonly used peak picking algorithms, including CentWave in XCMS, linear-weighted moving average in MS-DIAL, automated data analysis pipeline (ADAP) in MZmine 2, Savitzky-Golay in El-MAVEN, and FeatureFinderMetabo in OpenMS. We first collected 10 public metabolomics datasets representing various LC-MS analytical conditions. We then incorporated several novel strategies to (i) acquire the optimal peak picking parameters of each algorithm for a fair comparison, (ii) automatically recognize false metabolic features with poor chromatographic peak shapes, and (iii) evaluate the real metabolic features that are missed by the algorithms. By applying these strategies, we compared the true, false, and undetected metabolic features in each data processing outcome. Our results show that linear-weighted moving average consistently outperforms the other peak picking algorithms. To facilitate a mechanistic understanding of the differences, we proposed six peak attributes: ideal slope, sharpness, peak height, mass deviation, peak width, and scan number. We also developed an R program to automatically measure these attributes for detected and undetected true metabolic features. From the results of the 10 datasets, we concluded that four peak attributes, including ideal slope, scan number, peak width, and mass deviation, are critical for the detectability of a peak. For instance, the focus on ideal slope critically hinders the extraction of true metabolic features with low ideal slope scores in linear-weighted moving average, Savitzky-Golay, and ADAP. The relationships between peak picking algorithms and peak attributes were also visualized in a principal component analysis biplot. Overall, the clear comparison and explanation of the differences between peak picking algorithms can lead to the design of better peak picking strategies in the future.

AntDAS-DDA: A New Platform for Data-Dependent Acquisition Mode-Based Untargeted Metabolomic Profiling Analysis with Advantage of Recognizing Insource Fragment Ions to Improve Compound Identification

Data-dependent acquisition (DDA) mode in ultra-high-performance liquid chromatography-high-resolution mass spectrometry (UHPLC-HRMS) can provide massive amounts of MS¹ and MS/MS information of compounds in untargeted metabolomics and can thus facilitate compound identification greatly. In this work, we developed a new platform called AntDAS-DDA for the automatic processing of UHPLC-HRMS data sets acquired under the DDA mode. Several algorithms, including extracted ion chromatogram extraction, feature extraction, MS/MS spectrum construction, fragment ion identification, and MS¹ spectrum construction, were developed within the platform. The performance of AntDAS-DDA was investigated comprehensively with a mixture of standard and complex plant data sets. Results suggested that features in complex sample matrices can be extracted effectively, and the constructed MS¹ and MS/MS spectra can benefit in compound identification greatly. The efficiency of compound identification can be improved by about 20%. AntDAS-DDA can take full advantage of MS/MS information in multiple sample analyses and provide more MS/MS spectra than single sample analysis. A comparison with advanced data analysis tools indicated that AntDAS-DDA may be used as an alternative for routine UHPLC-HRMS-based untargeted metabolomics. AntDAS-DDA is freely available at http://www.pmdb.org.cn/antdasdda.

Integrated Exposomics/Metabolomics for Rapid Exposure and Effect Analyses

The totality of environmental exposures and lifestyle factors, commonly referred to as the exposome, is poorly understood. Measuring the myriad of chemicals that humans are exposed to is immensely challenging, and identifying disrupted metabolic pathways is even more complex. Here, we present a novel technological approach for the comprehensive, rapid, and integrated analysis of the endogenous human metabolome and the chemical exposome. By combining reverse-phase and hydrophilic interaction liquid chromatography (HILIC) and fast polarity-switching, molecules with highly diverse chemical structures can be analyzed in 15 min with a single analytical run as both column's effluents are combined before analysis. Standard reference materials and authentic standards were evaluated to critically benchmark performance. Highly sensitive median limits of detection (LODs) with 0.04 μM for >140 quantitatively assessed endogenous metabolites and 0.08 ng/mL for the >100 model xenobiotics and human estrogens in solvent were obtained. In matrix, the median LOD values were higher with 0.7 ng/mL (urine) and 0.5 ng/mL (plasma) for exogenous chemicals. To prove the dual-column approach's applicability, real-life urine samples from sub-Saharan Africa (high-exposure scenario) and Europe (low-exposure scenario) were assessed in a targeted and nontargeted manner. Our liquid chromatography high-resolution mass spectrometry (LC-HRMS) approach demonstrates the feasibility of quantitatively and simultaneously assessing the endogenous metabolome and the chemical exposome for the high-throughput measurement of environmental drivers of diseases.

Turning Metabolomics Data Processing from a “Black Box” to a “White Box”

Extracting thousands of metabolic features from liquid chromatography–mass spectrometry (LC–MS)–based metabolomics data is not easy. Although many feature extraction algorithms have been developed over the past few decades, automated feature extraction is still not a “white box” process. For instance, it is challenging to quickly determine the optimal parameters for the best feature extraction outcome. It is also impossible to extract every true metabolic feature. Moreover, there is contamination from false metabolic features of different sources, such as signal noise and in-source fragmentation. Our laboratory has recently developed a suite of bioinformatics tools to address these metabolic peak-picking challenges. The goal is to improve the peak-picking outcome quality, so we can effectively obtain biological information from the metabolomics data.

Integrated metabolomics and transcriptomic analysis of the flavonoid regulatory networks in Sorghum bicolor seeds

Background

The objective of this study was to reveal the flavonoid biosynthesis pathway in white (Z6), red (Z27) and black (HC4) seeds of the sweet sorghum (Sorghum bicolor) using metabolomics and transcriptomics, to identify different flavonoid metabolites, and to analyze the differentially expressed genes involved in flavonoid biosynthesis.

Results

We analyzed the metabolomics and transcriptomics data of sweet sorghum seeds. Six hundred and fifty-one metabolites including 171 flavonoids were identified in three samples. Integrated analysis of transcriptomics and metabolomics showed that 8 chalcone synthase genes (gene19114, gene19115, gene19116, gene19117, gene19118, gene19120, gene19122 and gene19123) involved in flavonoid biosynthesis, were identified and play central role in change of color. Six flavanone including homoeriodictyol, naringin, prunin, naringenin, hesperetin and pinocembrin were main reason for the color difference.

Conclusions

Our results provide valuable information on the flavonoid metabolites and the candidate genes involved in the flavonoid biosynthesis pathway in sweet sorghum seeds.

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.

MetEx: A Targeted Extraction Strategy for Improving the Coverage and Accuracy of Metabolite Annotation in Liquid Chromatography-High-Resolution Mass Spectrometry Data

Liquid chromatography-high-resolution mass spectrometry (LC-HRMS) is the most popular platform for untargeted metabolomics studies, but compound annotation is a challenge. In this work, we developed a new LC-HRMS data-targeted extraction method called MetEx for metabolite annotation. MetEx contains the retention time (t_R), MS1, and MS2 information of 30 620 metabolites from freely available spectral databases, including MoNA and KEGG. The t_R values of 95.4% of the compounds in our database were calculated by the GNN-RT model. The MS2 spectra of 39.4% compounds were also predicted using CFM-ID. MetEx was initially examined on a mixture of 634 standards, considering chemical coverage and accurate metabolite assignment, and later applied to human plasma (NIST SRM 1950), human urine, HepG2 cells, mouse liver tissue, and mouse feces. MetEx correctly assigned 252 out of 253 standards detected in our instruments. The platform also provided 8.0-44.2% more compounds in the biological samples compared to XCMS, MS-DIAL, and MZmine 2. MetEx is implemented and visualized in R and freely available at http://www.metaboex.cn/MetEx.

IDSL.IPA Characterizes the Organic Chemical Space in Untargeted LC/HRMS Data Sets

Generating comprehensive and high-fidelity metabolomics data matrices from LC/HRMS data remains to be extremely challenging for population-scale large studies (n > 200). Here, we present a new data processing pipeline, the Intrinsic Peak Analysis (IDSL.IPA) R package (https://ipa.idsl.me), to generate such data matrices specifically for organic compounds. The IDSL.IPA pipeline incorporates (1) identifying potential ¹²C and ¹³C ion pairs in individual mass spectra; (2) detecting and characterizing chromatographic peaks using a new sensitive and versatile approach to perform mass correction, peak smoothing, baseline development for local noise measurement, and peak quality determination; (3) correcting retention time and cross-referencing peaks from multiple samples by a dynamic retention index marker approach; (4) annotating peaks using a reference database of m/z and retention time; and (5) accelerating data processing using a parallel computation of the peak detection and alignment steps for larger studies. This pipeline has been successfully evaluated for studies ranging from 200 to 1600 samples. By specifically isolating high quality and reliable signals pertaining to carbon-containing compounds in untargeted LC/HRMS data sets from larger studies, IDSL.IPA opens new opportunities for discovering new biological insights in the population-scale metabolomics and exposomics projects. The package is available in the R CRAN repository at https://cran.r-project.org/package=IDSL.IPA.

JPA: Joint Metabolic Feature Extraction Increases the Depth of Chemical Coverage for LC-MS-Based Metabolomics and Exposomics

Extracting metabolic features from liquid chromatography-mass spectrometry (LC-MS) data has been a long-standing bioinformatic challenge in untargeted metabolomics. Conventional feature extraction algorithms fail to recognize features with low signal intensities, poor chromatographic peak shapes, or those that do not fit the parameter settings. This problem also poses a challenge for MS-based exposome studies, as low-abundant metabolic or exposomic features cannot be automatically recognized from raw data. To address this data processing challenge, we developed an R package, JPA (short for Joint Metabolomic Data Processing and Annotation), to comprehensively extract metabolic features from raw LC-MS data. JPA performs feature extraction by combining a conventional peak picking algorithm and strategies for (1) recognizing features with bad peak shapes but that have tandem mass spectra (MS²) and (2) picking up features from a user-defined targeted list. The performance of JPA in global metabolomics was demonstrated using serial diluted urine samples, in which JPA was able to rescue an average of 25% of metabolic features that were missed by the conventional peak picking algorithm due to dilution. More importantly, the chromatographic peak shapes, analytical accuracy, and precision of the rescued metabolic features were all evaluated. Furthermore, owing to its sensitive feature extraction, JPA was able to achieve a limit of detection (LOD) that was up to thousands of folds lower when automatically processing metabolomics data of a serial diluted metabolite standard mixture analyzed in HILIC(−) and RP(+) modes. Finally, the performance of JPA in exposome research was validated using a mixture of 250 drugs and 255 pesticides at environmentally relevant levels. JPA detected an average of 2.3-fold more exposure compounds than conventional peak picking only.

Fungi population metabolomics and molecular network study reveal novel biomarkers for early detection of aflatoxigenic Aspergillus species

Mycotoxins threaten global food safety, public health and cause huge socioeconomic losses. Early detection is an effective preventive strategy, yet efficient biomarkers for early detection of aflatoxigenic Aspergillus species are lacking. Here, we proposed to use untargeted metabolomics and machine learning to mine biomarkers of aflatoxigenic Aspergillus species. We systematically delineated metabolic differences across 568 extensive field sampling A. flavus and performed biomarker analysis. Versicolorin B, 11-hydroxy-O-methylsterigmatocystin et.al metabolites shown a high correlation (from 0.71 to 0.95) with strains aflatoxin-producing capacity. Molecular networking analysis deciphered the connection of aflatoxins and biomarkers as well as potential emerging mycotoxins. We then developed a model using the biomarkers as variables to discern aflatoxigenic Aspergillus species with 97.8% accuracy. A validation dataset and metabolome from other 16 fungal isolates confirmed the robustness and specificity of these biomarkers. We further demonstrated the solution feasibility in agricultural products by early detection of biomarkers, which predicted aflatoxin contamination risk 35-47 days in advance. A developed operable decision rule by the XGBoost algorithm help regulators to intuitively assess the risk prioritization with 87.2% accuracy. Our research provides novel insights into global food safety risk assessment which will be crucial for early prevention and control of mycotoxins.

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.

Deciphering Microbial Metal Toxicity Responses via Random Bar Code Transposon Site Sequencing and Activity-Based Metabolomics

To uncover metal toxicity targets and defense mechanisms of the facultative anaerobe Pantoea sp. strain MT58 (MT58), we used a multiomic strategy combining two global techniques, random bar code transposon site sequencing (RB-TnSeq) and activity-based metabolomics. MT58 is a metal-tolerant Oak Ridge Reservation (ORR) environmental isolate that was enriched in the presence of metals at concentrations measured in contaminated groundwater at an ORR nuclear waste site. The effects of three chemically different metals found at elevated concentrations in the ORR contaminated environment were investigated: the cation Al³⁺, the oxyanion CrO₄^2-, and the oxycation UO₂²⁺. Both global techniques were applied using all three metals under both aerobic and anaerobic conditions to elucidate metal interactions mediated through the activity of metabolites and key genes/proteins. These revealed that Al³⁺ binds intracellular arginine, CrO₄^2- enters the cell through sulfate transporters and oxidizes intracellular reduced thiols, and membrane-bound lipopolysaccharides protect the cell from UO₂²⁺ toxicity. In addition, the Tol outer membrane system contributed to the protection of cellular integrity from the toxic effects of all three metals. Likewise, we found evidence of regulation of lipid content in membranes under metal stress. Individually, RB-TnSeq and metabolomics are powerful tools to explore the impact various stresses have on biological systems. Here, we show that together they can be used synergistically to identify the molecular actors and mechanisms of these pertubations to an organism, furthering our understanding of how living systems interact with their environment. IMPORTANCE Studying microbial interactions with their environment can lead to a deeper understanding of biological molecular mechanisms. In this study, two global techniques, RB-TnSeq and activity metabolomics, were successfully used to probe the interactions between a metal-resistant microorganism, Pantoea sp. strain MT58, and metals contaminating a site where the organism can be located. A number of novel metal-microbe interactions were uncovered, including Al³⁺ toxicity targeting arginine synthesis, which could lead to a deeper understanding of the impact Al³⁺ contamination has on microbial communities as well as its impact on higher-level organisms, including plants for whom Al³⁺ contamination is an issue. Using multiomic approaches like the one described here is a way to further our understanding of microbial interactions and their impacts on the environment overall.

Cross platform solutions to improve the zebrafish polar metabolome coverage using LC-QTOF MS: Optimization of separation mechanisms, solvent additives, and resuspension solvents

Untargeted metabolomics has been widely used for studies with zebrafish embryos. Until now, the number of analytical approaches to determine metabolites in zebrafish is limited, and there is a lack of consensus on the best platforms for comprehensive metabolomics analysis of zebrafish embryos. In addition, the capacity of these methods to detect metabolites is unsatisfactory and the confidence level for identifying compounds is relatively low. To improve the metabolome coverage, we mainly focused on the optimization of separation mechanisms, mobile phase additives, and resuspension solvents based on liquid chromatography (LC) coupling to high-resolution mass spectrometry (HRMS) techniques. Moreover, the procedures for optimizing methods were assessed when taking metabolite profiles in both positive and negative ionization modes into account. Four LC columns were studied: C18, T3, PFP, and HILIC. In positive ionization mode, it was strongly recommended to employ the HILIC approach operated at the neutral condition, which led to the presence of more than 4700 features and the annotation of 151 metabolites, mainly zwitterionic and basic compounds, in comparison to reverse phase (RP)-based methods with less than 1000 features. In negative ionization mode, the PFP column operated at 0.02% acetic acid showed the best performance in terms of metabolite coverage: 3100 metabolic features were detected and 218 metabolites were annotated in zebrafish embryos. Metabolite profiles mainly contained acidic and zwitterionic compounds. HILIC-based platforms were complementary to RP columns when analyzing highly polar metabolites. Additionally, it was preferable to reconstitute zebrafish extracts in 100% water for analysis of metabolites on RP columns, with a 20-30% increase in the number of identified metabolites compared to a 50% water in methanol solution. However, water/methanol (1:9, v/v), as resuspension solution, was advantageous over water/methanol (1:1, v/v) for HILIC analysis showing an 8-15% increase in detected metabolites. In total 336 polar metabolites were annotated by the combination of the optimized HILIC (positive) and PFP (negative) approaches. The largest metabolome coverage of polar metabolites in zebrafish embryos was obtained when three approaches were combined (negative PFP and HILIC, and HILIC positive) resulting in more than 420 annotated compounds.

Global-Scale Metabolomic Profiling of Human Hair for Simultaneous Monitoring of Endogenous Metabolome, Short- and Long-Term Exposome

Hair is a unique biological matrix that adsorbs short-term exposures (e. g., environmental contaminants and personal care products) on its surface and also embeds endogenous metabolites and long-term exposures in its matrix. In this work, we developed an untargeted metabolomics workflow to profile both temporal exposure chemicals and endogenous metabolites in the same hair sample. This analytical workflow begins with the extraction of short-term exposures from hair surfaces through washing. Further development of mechanical homogenization extracts endogenous metabolites and long-term exposures from the cleaned hair. Both solutions of hair wash and hair extract were analyzed using ultra-high-performance liquid chromatography-high-resolution mass spectrometry (UHPLC-HRMS)-based metabolomics for global-scale metabolic profiling. After analysis, raw data were processed using bioinformatic programs recently developed specifically for exposome research. Using optimized experimental conditions, we detected a total of 10,005 and 9,584 metabolic features from hair wash and extraction samples, respectively. Among them, 274 and 276 features can be definitively confirmed by MS² spectral matching against spectral library, and an additional 3,356 and 3,079 features were tentatively confirmed as biotransformation metabolites. To demonstrate the performance of our hair metabolomics, we collected hair samples from three female volunteers and tested their hair metabolic changes before and after a 2-day exposure exercise. Our results show that 645 features from wash and 89 features from extract were significantly changed from the 2-day exposure. Altogether, this work provides a novel analytical approach to study the hair metabolome and exposome at a global scale, which can be implemented in a wide range of biological applications for a deeper understanding of the impact of environmental and genetic factors on human health.

Enhancing the power of liquid chromatography-Mass spectrometry for chemical fingerprinting of phytotoxins in the environment

Phytotoxins are plant secondary metabolites. They have recently been considered as chemicals of emerging concern (CECs) and there is a growing interest in their environmental fate and potential threat to public health. Dedicated target and non-target screening (NTS) analysis of phytotoxins in environmental samples are sparse, meanwhile phytotoxins are rarely detected in NTS-based analysis due to lack of an efficient methodology. Development of new analytical measurement methods is therefore highly needed. In this study, we for the first time investigated key parameters of reversed phase liquid chromatography-high resolution mass spectrometry (RPLC-HRMS) for five major classes of phytotoxins (alkaloids, steroids, terpenoids, flavonoids and aromatic polyketides) in environmental matrices; the investigation included analytical conditions which have not yet been explored by others, e.g. ionization at alkaline pH above 9. As the outcome we established a new analytical method for target/suspect screening and NTS of phytotoxins in the environment, which significantly improved the detection sensitivity with up to 40 times compared to previous methods, and enabled the discovery of over 30 phytotoxins in a NTS-based environmental study. We also observed that the negative ionization of phenols could be facilitated by the number of hydroxyl groups on the ring rather than their position of substitution. This study is of interest for a better fundamental understanding of the behavior of phytotoxins in LC-MS. Dedicated target/suspect screening and NTS methods will facilitate a better risk characterization of phytotoxins in the environment and stimulate implementation of new public regulation on phytotoxins.

A graph density-based strategy for features fusion from different peak extract software to achieve more metabolites in metabolic profiling from high-resolution mass spectrometry

In metabolomics study, it is not easy to extract the metabolites from data of ultra high-performance liquid chromatography-high-resolution mass spectrometry, especially for those with low abundance. Different software for peak recognition and matching use different algorithms, leading to different extract results. Therefore, integration of results from different software can obtain richer metabolome information, but the redundant features should be removed. In this study, an integrated strategy of fusing features and removing redundancy based on graph density (FRRGD) was proposed. A graph is used to cover the ion features generated by two open access software (XCMS, MZmine 2) and a software (SIEVE) from an instrument vendor, and redundant features were removed by searching the maximal complete sub-graphs. A standard mixture containing 41 metabolites and a spontaneous urine were utilized to develop the method and demonstrate its usefulness. For the standard mixture, 19, 19 and 27 metabolites were extracted by XCMS, MZmine 2 and SIEVE, respectively. After fusion by FRRGD, 37 metabolites were obtained. For the diluted spontaneous urine sample, 1103, 1500 and 387 metabolites were extracted by XCMS, MZmine 2 and SIEVE, respectively, FRRGD produced 1619 metabolites which were much more than individual software, significantly increasing metabolome coverage. The proposed FRRGD shows a great prospect as a new data processing strategy for metabolomics study.

From Samples to Insights into Metabolism: Uncovering Biologically Relevant Information in LC-HRMS Metabolomics Data

Untargeted metabolomics (including lipidomics) is a holistic approach to biomarker discovery and mechanistic insights into disease onset and progression, and response to intervention. Each step of the analytical and statistical pipeline is crucial for the generation of high-quality, robust data. Metabolite identification remains the bottleneck in these studies; therefore, confidence in the data produced is paramount in order to maximize the biological output. Here, we outline the key steps of the metabolomics workflow and provide details on important parameters and considerations. Studies should be designed carefully to ensure appropriate statistical power and adequate controls. Subsequent sample handling and preparation should avoid the introduction of bias, which can significantly affect downstream data interpretation. It is not possible to cover the entire metabolome with a single platform; therefore, the analytical platform should reflect the biological sample under investigation and the question(s) under consideration. The large, complex datasets produced need to be pre-processed in order to extract meaningful information. Finally, the most time-consuming steps are metabolite identification, as well as metabolic pathway and network analysis. Here we discuss some widely used tools and the pitfalls of each step of the workflow, with the ultimate aim of guiding the reader towards the most efficient pipeline for their metabolomics studies.