Network Marker Selection for Untargeted LC-MS Metabolomics Data

Bayesian functional analysis for untargeted metabolomics data with matching uncertainty and small sample sizes

Untargeted metabolomics based on liquid chromatography-mass spectrometry technology is quickly gaining widespread application, given its ability to depict the global metabolic pattern in biological samples. However, the data are noisy and plagued by the lack of clear identity of data features measured from samples. Multiple potential matchings exist between data features and known metabolites, while the truth can only be one-to-one matches. Some existing methods attempt to reduce the matching uncertainty, but are far from being able to remove the uncertainty for most features. The existence of the uncertainty causes major difficulty in downstream functional analysis. To address these issues, we develop a novel approach for Bayesian Analysis of Untargeted Metabolomics data (BAUM) to integrate previously separate tasks into a single framework, including matching uncertainty inference, metabolite selection and functional analysis. By incorporating the knowledge graph between variables and using relatively simple assumptions, BAUM can analyze datasets with small sample sizes. By allowing different confidence levels of feature-metabolite matching, the method is applicable to datasets in which feature identities are partially known. Simulation studies demonstrate that, compared with other existing methods, BAUM achieves better accuracy in selecting important metabolites that tend to be functionally consistent and assigning confidence scores to feature-metabolite matches. We analyze a COVID-19 metabolomics dataset and a mouse brain metabolomics dataset using BAUM. Even with a very small sample size of 16 mice per group, BAUM is robust and stable. It finds pathways that conform to existing knowledge, as well as novel pathways that are biologically plausible.

An integrated deep learning framework for the interpretation of untargeted metabolomics data

Untargeted metabolomics is gaining widespread applications. The key aspects of the data analysis include modeling complex activities of the metabolic network, selecting metabolites associated with clinical outcome and finding critical metabolic pathways to reveal biological mechanisms. One of the key roadblocks in data analysis is not well-addressed, which is the problem of matching uncertainty between data features and known metabolites. Given the limitations of the experimental technology, the identities of data features cannot be directly revealed in the data. The predominant approach for mapping features to metabolites is to match the mass-to-charge ratio (m/z) of data features to those derived from theoretical values of known metabolites. The relationship between features and metabolites is not one-to-one since some metabolites share molecular composition, and various adduct ions can be derived from the same metabolite. This matching uncertainty causes unreliable metabolite selection and functional analysis results. Here we introduce an integrated deep learning framework for metabolomics data that take matching uncertainty into consideration. The model is devised with a gradual sparsification neural network based on the known metabolic network and the annotation relationship between features and metabolites. This architecture characterizes metabolomics data and reflects the modular structure of biological system. Three goals can be achieved simultaneously without requiring much complex inference and additional assumptions: (1) evaluate metabolite importance, (2) infer feature-metabolite matching likelihood and (3) select disease sub-networks. When applied to a COVID metabolomics dataset and an aging mouse brain dataset, our method found metabolic sub-networks that were easily interpretable.

Metapone: a Bioconductor package for joint pathway testing for untargeted metabolomics data

MOTIVATION

Testing for pathway enrichment is an important aspect in the analysis of untargeted metabolomics data. Due to the unique characteristics of untargeted metabolomics data, some key issues have not been fully addressed in existing pathway testing algorithms: (1) matching uncertainty between data features and metabolites; (2) lacking of method to analyze positive mode and negative mode LC/MS data simultaneously on the same set of subjects; (3) the incompleteness of pathways in individual software packages.

RESULTS

We developed an innovative R/Bioconductor package: metabolic pathway testing with positive and negative mode data (metapone), which can perform two novel statistical tests that take matching uncertainty into consideration - (1) a weighted GSEA-type test, and (2) a permutation-based weighted hypergeometric test. The package is capable of combining positive and negative ion mode results in a single testing scheme. For comprehensiveness, the built-in pathways were manually curated from three sources: KEGG, Mummichog, and SMPDB.

AVAILABILITY

The package is available at https://bioconductor.org/packages/devel/bioc/html/metapone.html.

SUPPLEMENTARY INFORMATION

Supplementary data are available at Bioinformatics online.

Lipid Profiles of Human Serum Fractions Enhanced with CD9 Antibody-Immobilized Magnetic Beads

Blood samples are minimally invasive and can be collected repeatedly, but they are far from the site of disease and the target molecules are diluted by the large amount of blood. Therefore, we performed lipidomics using immunoprecipitation as a method to enrich specific fractions of serum. In this study, a CD9 antibody was immobilized on magnetic beads to enrich CD9-containing components in the serum for lipidomics. The percentages of phospholipids recovered from serum by methanol and isopropanol extractions were not significantly different, but triglycerides were barely recovered from serum by methanol extraction, requiring the use of isopropanol. However, once the serum was enriched with CD9 magnetic beads, triglycerides, and phospholipids were recovered at similar levels in both methanol and isopropanol extractions. Therefore, it is possible that the triglyceride fraction of the whole serum and the triglyceride fraction were enriched in CD9 magnetic beads differ in localization and properties. In addition, the variation per disease was small in general serum lipidomics; however, the difference per disease appeared larger when CD9 magnetic bead enrichment was employed.

Multi-Omics Analysis to Generate Hypotheses for Mild Health Problems in Monkeys

Certain symptoms associated with mild sickness and lethargy have not been categorized as definitive diseases. Confirming such symptoms in captive monkeys (Macaca fascicularis, known as cynomolgus monkeys) can be difficult; however, it is possible to observe and analyze their feces. In this study, we investigated the relationship between stool state and various omics data by considering objective and quantitative values of stool water content as a phenotype for analysis. By examining the food intake of the monkeys and assessing their stool, urine, and plasma, we attempted to obtain a comprehensive understanding of the health status of individual monkeys and correlate it with the stool condition. Our metabolomics data strongly suggested that many lipid-related metabolites were correlated with the stool water content. The lipidomic analysis revealed the involvement of saturated and oxidized fatty acids, metallomics revealed the contribution of selenium (a bio-essential trace element), and intestinal microbiota analysis revealed the association of several bacterial species with the stool water content. Based on our results, we hypothesize that the redox imbalance causes minor health problems. However, it is not possible to make a definite conclusion using multi-omics alone, and other hypotheses could be proposed.

DNLC: differential network local consistency analysis

BACKGROUND

The biological network is highly dynamic. Functional relations between genes can be activated or deactivated depending on the biological conditions. On the genome-scale network, subnetworks that gain or lose local expression consistency may shed light on the regulatory mechanisms related to the changing biological conditions, such as disease status or tissue developmental stages.

RESULTS

In this study, we develop a new method to select genes and modules on the existing biological network, in which local expression consistency changes significantly between clinical conditions. The method is called DNLC: Differential Network Local Consistency. In simulations, our algorithm detected artificially created local consistency changes effectively. We applied the method on two publicly available datasets, and the method detected novel genes and network modules that were biologically plausible.

CONCLUSIONS

The new method is effective in finding modules in which the gene expression consistency change between clinical conditions. It is a useful tool that complements traditional differential expression analyses to make discoveries from gene expression data. The R package is available at https://cran.r-project.org/web/packages/DNLC.

Inborn Errors of Metabolism in the Era of Untargeted Metabolomics and Lipidomics

Inborn errors of metabolism (IEMs) are a group of inherited diseases with variable incidences. IEMs are caused by disrupting enzyme activities in specific metabolic pathways by genetic mutations, either directly or indirectly by cofactor deficiencies, causing altered levels of compounds associated with these pathways. While IEMs may present with multiple overlapping symptoms and metabolites, early and accurate diagnosis of IEMs is critical for the long-term health of affected subjects. The prevalence of IEMs differs between countries, likely because different IEM classifications and IEM screening methods are used. Currently, newborn screening programs exclusively use targeted metabolic assays that focus on limited panels of compounds for selected IEM diseases. Such targeted approaches face the problem of false negative and false positive diagnoses that could be overcome if metabolic screening adopted analyses of a broader range of analytes. Hence, we here review the prospects of using untargeted metabolomics for IEM screening. Untargeted metabolomics and lipidomics do not rely on predefined target lists and can detect as many metabolites as possible in a sample, allowing to screen for many metabolic pathways simultaneously. Examples are given for nontargeted analyses of IEMs, and prospects and limitations of different metabolomics methods are discussed. We conclude that dedicated studies are needed to compare accuracy and robustness of targeted and untargeted methods with respect to widening the scope of IEM diagnostics.

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.

Missing value imputation for LC-MS metabolomics data by incorporating metabolic network and adduct ion relations

Motivation

Metabolomics data generated from liquid chromatography-mass spectrometry platforms often contain missing values. Existing imputation methods do not consider underlying feature relations and the metabolic network information. As a result, the imputation results may not be optimal.

Results

We proposed an imputation algorithm that incorporates the existing metabolic network, adduct ion relations even for unknown compounds, as well as linear and nonlinear associations between feature intensities to build a feature-level network. The algorithm uses support vector regression for missing value imputation based on features in the neighborhood on the network. We compared our proposed method with methods being widely used. As judged by the normalized root mean squared error in real data-based simulations, our proposed methods can achieve better accuracy.

Availability and implementation

The R package is available at http://web1.sph.emory.edu/users/tyu8/MINMA.

Contact

jiankang@umich.edu or tianwei.yu@emory.edu.

Supplementary information

Supplementary data are available at Bioinformatics online.

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.