compMS2Miner: An Automatable Metabolite Identification, Visualization, and Data-Sharing R Package for High-Resolution LC-MS Data Sets

Integrated annotation prioritizes metabolites with bioactivity in inflammatory bowel disease

Microbial biochemistry is central to the pathophysiology of inflammatory bowel diseases (IBD). Improved knowledge of microbial metabolites and their immunomodulatory roles is thus necessary for diagnosis and management. Here, we systematically analyzed the chemical, ecological, and epidemiological properties of ~82k metabolic features in 546 Integrative Human Microbiome Project (iHMP/HMP2) metabolomes, using a newly developed methodology for bioactive compound prioritization from microbial communities. This suggested >1000 metabolic features as potentially bioactive in IBD and associated ~43% of prevalent, unannotated features with at least one well-characterized metabolite, thereby providing initial information for further characterization of a significant portion of the fecal metabolome. Prioritized features included known IBD-linked chemical families such as bile acids and short-chain fatty acids, and less-explored bilirubin, polyamine, and vitamin derivatives, and other microbial products. One of these, nicotinamide riboside, reduced colitis scores in DSS-treated mice. The method, MACARRoN, is generalizable with the potential to improve microbial community characterization and provide therapeutic candidates.

Progress and opportunities in microbial community metabolomics

The metabolome lies at the interface of host-microbiome crosstalk. Previous work has established links between chemically diverse microbial metabolites and a myriad of host physiological processes and diseases. Coupled with scalable and cost-effective technologies, metabolomics is thus gaining popularity as a tool for characterization of microbial communities, particularly when combined with metagenomics as a window into microbiome function. A systematic interrogation of microbial community metabolomes can uncover key microbial compounds, metabolic capabilities of the microbiome, and also provide critical mechanistic insights into microbiome-linked host phenotypes. In this review, we discuss methods and accompanying resources that have been developed for these purposes. The accomplishments of these methods demonstrate that metabolomes can be used to functionally characterize microbial communities, and that microbial properties can be used to identify and investigate chemical compounds.

HRPIF data mining based on data-dependent/independent acquisition for Rhei Radix et Rhizoma metabolite screening in rats

In traditional Chinese medicine (TCM), components with identical nuclei often share structural similarity, indicating the possibility of similar second-level mass spectrometry (MS/MS) fragments. High-resolution product-ion filter (HRPIF) technique can be utilized to identify metabolites, with similar fragments, in vivo. In principle, this technique applies to TCM; however, its application has been restricted due to the limitations of traditional MS/MS data acquisition. Therefore, a novel analysis strategy, based on data-dependent acquisition (DDA) and data-independent acquisition (DIA) datasets, has been developed for the determination of template product ions and efficient non-targeted identification of TCM-related components in vivo by HRPIF and background subtraction (BS). This DDA-DIA combination strategy, taking Rhei Radix et Rhizoma as a test case, identified 71 anthraquinone prototype components in vitro (36 of which were discovered for the first time), and 45 related components in vivo, confirming glucuronidation and sulfation as the main reactions. The developed strategy could rapidly identify TCM-related components in vivo with high sensitivity, indicating the immense importance of this novel HRPIF data mining technology in TCM analysis.

A post processing strategy to score and rank the annotation confidence of saponins in natural products by integrating MS2 spectral similarity and fragment interpretation

Tandem mass spectrometry-spectra-based annotation in natural products challenges a lot because of ambiguous structural characterization. It still lacks an efficiency method to score and rank the annotation confidence. Herein, we develop a novel approach to rank the annotation confidences of saponins. Annotations were accomplished according to fragmentation patterns. The corresponding diagnostic fragments and their abundances were recorded. Average abundances were taken as a reference spectrum, and the cosine similarity score (CSS) was calculated to measure how well the spectral matched. According to CSS values, statistic description for confidence levels can be effectively provided. Next, the fragment interpretation score (FIS) was proposed to investigate the deviators' characteristic fragmentation. FIS offset the effect from the deviators' unique fragments. Suspicious annotations involving low CSS and high FIS, may derived from the MS² spectral background interferences or co-elution. Annotations with low CSS and FIS rank as low confidences, as these annotations need more attention. Using this method, novel saccharide sequences, specific fragmentation preferences, undistinguished precursors, even new structures can also be well traced. By proposed new scoring system, confidence evaluations can be ranked, resulting in significantly enhanced annotation reliability.

Untargeted In Silico Compound Classification-A Novel Metabolomics Method to Assess the Chemodiversity in Bryophytes

In plant ecology, biochemical analyses of bryophytes and vascular plants are often conducted on dried herbarium specimen as species typically grow in distant and inaccessible locations. Here, we present an automated in silico compound classification framework to annotate metabolites using an untargeted data independent acquisition (DIA)–LC/MS–QToF-sequential windowed acquisition of all theoretical fragment ion mass spectra (SWATH) ecometabolomics analytical method. We perform a comparative investigation of the chemical diversity at the global level and the composition of metabolite families in ten different species of bryophytes using fresh samples collected on-site and dried specimen stored in a herbarium for half a year. Shannon and Pielou’s diversity indices, hierarchical clustering analysis (HCA), sparse partial least squares discriminant analysis (sPLS-DA), distance-based redundancy analysis (dbRDA), ANOVA with post-hoc Tukey honestly significant difference (HSD) test, and the Fisher’s exact test were used to determine differences in the richness and composition of metabolite families, with regard to herbarium conditions, ecological characteristics, and species. We functionally annotated metabolite families to biochemical processes related to the structural integrity of membranes and cell walls (proto-lignin, glycerophospholipids, carbohydrates), chemical defense (polyphenols, steroids), reactive oxygen species (ROS) protection (alkaloids, amino acids, flavonoids), nutrition (nitrogen- and phosphate-containing glycerophospholipids), and photosynthesis. Changes in the composition of metabolite families also explained variance related to ecological functioning like physiological adaptations of bryophytes to dry environments (proteins, peptides, flavonoids, terpenes), light availability (flavonoids, terpenes, carbohydrates), temperature (flavonoids), and biotic interactions (steroids, terpenes). The results from this study allow to construct chemical traits that can be attributed to biogeochemistry, habitat conditions, environmental changes and biotic interactions. Our classification framework accelerates the complex annotation process in metabolomics and can be used to simplify biochemical patterns. We show that compound classification is a powerful tool that allows to explore relationships in both molecular biology by “zooming in” and in ecology by “zooming out”. The insights revealed by our framework allow to construct new research hypotheses and to enable detailed follow-up studies.

Salivary metabolomics for cancer detection

Introduction: Saliva is an ideal biofluid that can be collected in a noninvasive manner, enabling safe and frequent screening of various diseases. Recent studies have revealed that salivary metabolomics analysis has the potential to detect both oral and systemic cancers. Area covered: We reviewed the technical aspects, as well as applications, of salivary metabolomics for cancer detection. The topics include the effects of preconditioning and the method of sample collection, sample storage, processing, measurement, data analysis, and validation of the results. We also examined the rational relationship between salivary biomarkers and tumors distant from the oral cavity. A strategy to establish standard operating protocols for obtaining reproducible quantification data is also discussed Expert opinion: Salivary metabolomics reflects oral and systematic health status, which potently enables cancer detection. The sensitivity and specificity of each marker and their combinations have been well evaluated, but a validation study is required. Further, the standard operating protocol for each procedure should be established to obtain reproducible data before clinical usage.

MS-CleanR: A Feature-Filtering Workflow for Untargeted LC-MS Based Metabolomics

Untargeted metabolomics using liquid chromatography-mass spectrometry (LC-MS) is currently the gold-standard technique to determine the full chemical diversity in biological samples. However, this approach still has many limitations; notably, the difficulty of accurately estimating the number of unique metabolites profiled among the thousands of MS ion signals arising from chromatograms. Here, we describe a new workflow, MS-CleanR, based on the MS-DIAL/MS-FINDER suite, which tackles feature degeneracy and improves annotation rates. We show that implementation of MS-CleanR reduces the number of signals by nearly 80% while retaining 95% of unique metabolite features. Moreover, the annotation results from MS-FINDER can be ranked according to the database chosen by the user, which enhance identification accuracy. Application of MS-CleanR to the analysis of Arabidopsis thaliana grown in three different conditions fostered class separation resulting from multivariate data analysis and led to annotation of 75% of the final features. The full workflow was applied to metabolomic profiles from three strains of the leguminous plant Medicago truncatula that have different susceptibilities to the oomycete pathogen Aphanomyces euteiches. A group of glycosylated triterpenoids overrepresented in resistant lines were identified as candidate compounds conferring pathogen resistance. MS-CleanR is implemented through a Shiny interface for intuitive use by end-users (available at https://github.com/eMetaboHUB/MS-CleanR).

CliqueMS: a computational tool for annotating in-source metabolite ions from LC-MS untargeted metabolomics data based on a coelution similarity network

Motivation

The analysis of biological samples in untargeted metabolomic studies using LC-MS yields tens of thousands of ion signals. Annotating these features is of the utmost importance for answering questions as fundamental as, e.g. how many metabolites are there in a given sample.

Results

Here, we introduce CliqueMS, a new algorithm for annotating in-source LC-MS1 data. CliqueMS is based on the similarity between coelution profiles and therefore, as opposed to most methods, allows for the annotation of a single spectrum. Furthermore, CliqueMS improves upon the state of the art in several dimensions: (i) it uses a more discriminatory feature similarity metric; (ii) it treats the similarities between features in a transparent way by means of a simple generative model; (iii) it uses a well-grounded maximum likelihood inference approach to group features; (iv) it uses empirical adduct frequencies to identify the parental mass and (v) it deals more flexibly with the identification of the parental mass by proposing and ranking alternative annotations. We validate our approach with simple mixtures of standards and with real complex biological samples. CliqueMS reduces the thousands of features typically obtained in complex samples to hundreds of metabolites, and it is able to correctly annotate more metabolites and adducts from a single spectrum than available tools.

Availability and implementation

https://CRAN.R-project.org/package=cliqueMS and https://github.com/osenan/cliqueMS.

Supplementary information

Supplementary data are available at Bioinformatics online.

The metaRbolomics Toolbox in Bioconductor and beyond

Metabolomics aims to measure and characterise the complex composition of metabolites in a biological system. Metabolomics studies involve sophisticated analytical techniques such as mass spectrometry and nuclear magnetic resonance spectroscopy, and generate large amounts of high-dimensional and complex experimental data. Open source processing and analysis tools are of major interest in light of innovative, open and reproducible science. The scientific community has developed a wide range of open source software, providing freely available advanced processing and analysis approaches. The programming and statistics environment R has emerged as one of the most popular environments to process and analyse Metabolomics datasets. A major benefit of such an environment is the possibility of connecting different tools into more complex workflows. Combining reusable data processing R scripts with the experimental data thus allows for open, reproducible research. This review provides an extensive overview of existing packages in R for different steps in a typical computational metabolomics workflow, including data processing, biostatistics, metabolite annotation and identification, and biochemical network and pathway analysis. Multifunctional workflows, possible user interfaces and integration into workflow management systems are also reviewed. In total, this review summarises more than two hundred metabolomics specific packages primarily available on CRAN, Bioconductor and GitHub.

MolNetEnhancer: Enhanced Molecular Networks by Integrating Metabolome Mining and Annotation Tools

Metabolomics has started to embrace computational approaches for chemical interpretation of large data sets. Yet, metabolite annotation remains a key challenge. Recently, molecular networking and MS2LDA emerged as molecular mining tools that find molecular families and substructures in mass spectrometry fragmentation data. Moreover, in silico annotation tools obtain and rank candidate molecules for fragmentation spectra. Ideally, all structural information obtained and inferred from these computational tools could be combined to increase the resulting chemical insight one can obtain from a data set. However, integration is currently hampered as each tool has its own output format and efficient matching of data across these tools is lacking. Here, we introduce MolNetEnhancer, a workflow that combines the outputs from molecular networking, MS2LDA, in silico annotation tools (such as Network Annotation Propagation or DEREPLICATOR), and the automated chemical classification through ClassyFire to provide a more comprehensive chemical overview of metabolomics data whilst at the same time illuminating structural details for each fragmentation spectrum. We present examples from four plant and bacterial case studies and show how MolNetEnhancer enables the chemical annotation, visualization, and discovery of the subtle substructural diversity within molecular families. We conclude that MolNetEnhancer is a useful tool that greatly assists the metabolomics researcher in deciphering the metabolome through combination of multiple independent in silico pipelines.

Metabolomics of neonatal blood spots reveal distinct phenotypes of pediatric acute lymphoblastic leukemia and potential effects of early-life nutrition

Early-life exposures are believed to influence the incidence of pediatric acute lymphoblastic leukemia (ALL). Archived neonatal blood spots (NBS), collected within the first days of life, offer a means to investigate small molecules that reflect early-life exposures. Using untargeted metabolomics, we compared abundances of small-molecule features in extracts of NBS punches from 332 children that later developed ALL and 324 healthy controls. Subjects were stratified by early (1-5 y) and late (6-14 y) diagnosis. Mutually-exclusive sets of metabolic features - representing putative lipids and fatty acids - were associated with ALL, including 9 and 19 metabolites in the early- and late-diagnosis groups, respectively. In the late-diagnosis group, a prominent cluster of features with apparent 18:2 fatty-acid chains suggested that newborn exposure to the essential nutrient, linoleic acid, increased ALL risk. Interestingly, abundances of these putative 18:2 lipids were greater in infants who were fed formula rather than breast milk (colostrum) and increased with the mother's pre-pregnancy body mass index. These results suggest possible etiologic roles of newborn nutrition in late-diagnosis ALL.

Untargeted lipidomic features associated with colorectal cancer in a prospective cohort

BACKGROUND

Epidemiologists are beginning to employ metabolomics and lipidomics with archived blood from incident cases and controls to discover causes of cancer. Although several such studies have focused on colorectal cancer (CRC), they all followed targeted or semi-targeted designs that limited their ability to find discriminating molecules and pathways related to the causes of CRC.

METHODS

Using an untargeted design, we measured lipophilic metabolites in prediagnostic serum from 66 CRC patients and 66 matched controls from the European Prospective Investigation into Cancer and Nutrition (Turin, Italy). Samples were analyzed by liquid chromatography-high-resolution mass spectrometry (LC-MS), resulting in 8690 features for statistical analysis.

RESULTS

Rather than the usual multiple-hypothesis-testing approach, we based variable selection on an ensemble of regression methods, which found nine features to be associated with case-control status. We then regressed each selected feature on time-to-diagnosis to determine whether the feature was likely to be either a potentially causal biomarker or a reactive product of disease progression (reverse causality).

CONCLUSIONS

Of the nine selected LC-MS features, four appear to be involved in CRC etiology and merit further investigation in prospective studies of CRC. Four other features appear to be related to progression of the disease (reverse causality), and may represent biomarkers of value for early detection of CRC.

Integrating tools for non-targeted analysis research and chemical safety evaluations at the US EPA

Tens-of-thousands of chemicals are registered in the U.S. for use in countless processes and products. Recent evidence suggests that many of these chemicals are measureable in environmental and/or biological systems, indicating the potential for widespread exposures. Traditional public health research tools, including in vivo studies and targeted analytical chemistry methods, have been unable to meet the needs of screening programs designed to evaluate chemical safety. As such, new tools have been developed to enable rapid assessment of potentially harmful chemical exposures and their attendant biological responses. One group of tools, known as "non-targeted analysis" (NTA) methods, allows the rapid characterization of thousands of never-before-studied compounds in a wide variety of environmental, residential, and biological media. This article discusses current applications of NTA methods, challenges to their effective use in chemical screening studies, and ways in which shared resources (e.g., chemical standards, databases, model predictions, and media measurements) can advance their use in risk-based chemical prioritization. A brief review is provided of resources and projects within EPA's Office of Research and Development (ORD) that provide benefit to, and receive benefits from, NTA research endeavors. A summary of EPA's Non-Targeted Analysis Collaborative Trial (ENTACT) is also given, which makes direct use of ORD resources to benefit the global NTA research community. Finally, a research framework is described that shows how NTA methods will bridge chemical prioritization efforts within ORD. This framework exists as a guide for institutions seeking to understand the complexity of chemical exposures, and the impact of these exposures on living systems.

Redefining environmental exposure for disease etiology

Etiological studies of human exposures to environmental factors typically rely on low-throughput methods that target only a few hundred chemicals or mixtures. In this Perspectives article, I outline how environmental exposure can be defined by the blood exposome—the totality of chemicals circulating in blood. The blood exposome consists of chemicals derived from both endogenous and exogenous sources. Endogenous chemicals are represented by the human proteome and metabolome, which establish homeostatic networks of functional molecules. Exogenous chemicals arise from diet, vitamins, drugs, pathogens, microbiota, pollution, and lifestyle factors, and can be measured in blood as subsets of the proteome, metabolome, metals, macromolecular adducts, and foreign DNA and RNA. To conduct ‘exposome-wide association studies’, blood samples should be obtained prospectively from subjects—preferably at critical stages of life—and then analyzed in incident disease cases and matched controls to find discriminating exposures. Results from recent metabolomic investigations of archived blood illustrate our ability to discover potentially causal exposures with current technologies.

A Suspect Screening Method for Characterizing Multiple Chemical Exposures among a Demographically Diverse Population of Pregnant Women in San Francisco

BACKGROUND

In utero exposure to environmental chemicals can adversely impact pregnancy outcomes and childhood health, but minimal biomonitoring data exist on the majority of chemicals used in commerce.

OBJECTIVES

We aimed to profile exposure to multiple environmental organic acids (EOAs) and identify novel chemicals that have not been previously biomonitored in a diverse population of pregnant women.

METHODS

We used liquid chromatography-quadrupole time-of-flight mass spectrometry (LC-QTOF/MS) to perform a suspect screen for 696 EOAs, (e.g., phenols and phthalate metabolites) on the maternal serum collected at delivery from 75 pregnant women delivering at two large San Francisco Hospitals. We examined demographic differences in peak areas and detection frequency (DF) of suspect EOAs using a Kruskal-Wallis Rank Sum test or Fisher's exact test. We confirmed selected suspects by comparison with their respective reference standards.

RESULTS

We detected, on average, 56 [standard deviation (SD)]: 8) suspect EOAs in each sample (range: 32-73). Twelve suspect EOAs with DF≥60 were matched to 21 candidate compounds in our EOA database, two-thirds of which are novel chemicals. We found demographic differences in DF for 13 suspect EOAs and confirmed the presence of 6 priority novel chemicals: 2,4-Di-tert-butylphenol, Pyrocatechol, 2,4-Dinitrophenol, 3,5-Di-tert-butylsalicylic acid, 4-Hydroxycoumarin, and 2'-Hydroxyacetophenone (or 3'-Hydroxyacetophenone). The first two are high-production-volume chemicals in the United States.

CONCLUSION

Suspect screening in human biomonitoring provides a viable method to characterize a broad spectrum of environmental chemicals to prioritize for targeted method development and quantification. https://doi.org/10.1289/EHP2920.

New tools and resources in metabolomics: 2016-2017

Rapid advances in mass spectrometry (MS) and nuclear magnetic resonance (NMR)-based platforms for metabolomics have led to an upsurge of data every single year. Newer high-throughput platforms, hyphenated technologies, miniaturization, and tool kits in data acquisition efforts in metabolomics have led to additional challenges in metabolomics data pre-processing, analysis, interpretation, and integration. Thanks to the informatics, statistics, and computational community, new resources continue to develop for metabolomics researchers. The purpose of this review is to provide a summary of the metabolomics tools, software, and databases that were developed or improved during 2016-2017, thus, enabling readers, developers, and researchers access to a succinct but thorough list of resources for further improvisation, implementation, and application in due course of time.

Data acquisition workflows in liquid chromatography coupled to high resolution mass spectrometry-based metabolomics: Where do we stand?

Typical mass spectrometry (MS) based untargeted metabolomics protocols are tedious as well as time- and sample-consuming. In particular, they often rely on "full-scan-only" analyses using liquid chromatography (LC) coupled to high resolution mass spectrometry (HRMS) from which metabolites of interest are first highlighted, and then tentatively identified by using targeted MS/MS experiments. However, this situation is evolving with the emergence of integrated HRMS based-data acquisition protocols able to perform multi-event acquisitions. Most of these protocols, referring to as data dependent and data independent acquisition (DDA and DIA, respectively), have been initially developed for proteomic applications and have recently demonstrated their applicability to biomedical studies. In this context, the aim of this article is to take stock of the progress made in the field of DDA- and DIA-based protocols, and evaluate their ability to change conventional metabolomic and lipidomic data acquisition workflows, through a review of HRMS instrumentation, DDA and DIA workflows, and also associated informatics tools.

Development of a Liquid Chromatography-High Resolution Mass Spectrometry Metabolomics Method with High Specificity for Metabolite Identification Using All Ion Fragmentation Acquisition

High-resolution mass spectrometry (HRMS)-based metabolomics approaches have made significant advances. However, metabolite identification is still a major challenge with significant bottleneck in translating metabolomics data into biological context. In the current study, a liquid chromatography (LC)-HRMS metabolomics method was developed using an all ion fragmentation (AIF) acquisition approach. To increase the specificity in metabolite annotation, four criteria were considered: (i) accurate mass (AM), (ii) retention time (RT), (iii) MS/MS spectrum, and (iv) product/precursor ion intensity ratios. We constructed an in-house mass spectral library of 408 metabolites containing AMRT and MS/MS spectra information at four collision energies. The percent relative standard deviations between ion ratios of a metabolite in an analytical standard vs sample matrix were used as an additional metric for establishing metabolite identity. A data processing method for targeted metabolite screening was then created, merging m/z, RT, MS/MS, and ion ratio information for each of the 413 metabolites. In the data processing method, the precursor ion and product ion were considered as the quantifier and qualifier ion, respectively. We also included a scheme to distinguish coeluting isobaric compounds by selecting a specific product ion as the quantifier ion instead of the precursor ion. An advantage of the current AIF approach is the concurrent collection of full scan data, enabling identification of metabolites not included in the database. Our data acquisition strategy enables a simultaneous mixture of database-dependent targeted and nontargeted metabolomics in combination with improved accuracy in metabolite identification, increasing the quality of the biological information acquired in a metabolomics experiment.