Large-scale Metabolomics Data Research Articles

MetMiner: A user-friendly pipeline for large-scale plant metabolomics data analysis.

The utilization of metabolomics approaches to explore the metabolic mechanisms underlying plant fitness and adaptation to dynamic environments is growing, highlighting the need for an efficient and user-friendly toolkit tailored for analyzing the extensive datasets generated by metabolomics studies. Current protocols for metabolome data analysis often struggle with handling large-scale datasets or require programming skills. To address this, we present MetMiner (https://github.com/ShawnWx2019/MetMiner), a user-friendly, full-functionality pipeline specifically designed for plant metabolomics data analysis. Built on R shiny, MetMiner can be deployed on servers to utilize additional computational resources for processing large-scale datasets. MetMiner ensures transparency, traceability, and reproducibility throughout the analytical process. Its intuitive interface provides robust data interaction and graphical capabilities, enabling users without prior programming skills to engage deeply in data analysis. Additionally, we constructed and integrated a plant-specific mass spectrometry database into the MetMiner pipeline to optimize metabolite annotation. We have also developed MDAtoolkits, which include a complete set of tools for statistical analysis, metabolite classification, and enrichment analysis, to facilitate the mining of biological meaning from the datasets. Moreover, we propose an iterative weighted gene co-expression network analysis strategy for efficient biomarker metabolite screening in large-scale metabolomics data mining. In two case studies, we validated MetMiner's efficiency in data mining and robustness in metabolite annotation. Together, the MetMiner pipeline represents a promising solution for plant metabolomics analysis, providing a valuable tool for the scientific community to use with ease.

MetaboLink: A web application for Streamlined Processing and Analysis of Large-Scale Untargeted Metabolomics Data.

The post-processing and analysis of large-scale untargeted metabolomics data face significant challenges due to the intricate nature of correction, filtration, imputation, and normalization steps. Manual execution across various applications often leads to inefficiencies, human-induced errors, and inconsistencies within the workflow. Addressing these issues, we introduce MetaboLink, a novel web application designed to process LC-MS metabolomics datasets combining established methodologies and offering flexibility and ease of implementation. It offers visualization options for data interpretation, an interface for statistical testing, and integration with PolySTest for further tests and with VSClust for clustering analysis. Fully functional tool is publicly available at https://computproteomics.bmb.sdu.dk/Metabolomics/. The source code is available at https://github.com/anitamnd/MetaboLink and a detailed description of the app can be found at https://github.com/anitamnd/MetaboLink/wiki. A tutorial video can be found at https://youtu.be/ZM6j10S6Z8Q. Supplementary data are available at Bioinformatics online.

The underappreciated diversity of bile acid modifications

The repertoire of modifications to bile acids and related steroidal lipids by host and microbial metabolism remains incompletely characterized. To address this knowledge gap, we created a reusable resource of tandem mass spectrometry (MS/MS) spectra by filtering 1.2 billion publicly available MS/MS spectra for bile-acid-selective ion patterns. Thousands of modifications are distributed throughout animal and human bodies as well as microbial cultures. We employed this MS/MS library to identify polyamine bile amidates, prevalent in carnivores. They are present in humans, and their levels alter with a diet change from a Mediterranean to a typical American diet. This work highlights the existence of many more bile acid modifications than previously recognized and the value of leveraging public large-scale untargeted metabolomics data to discover metabolites. The availability of a modification-centric bile acid MS/MS library will inform future studies investigating bile acid roles in health and disease.

COVRECON: automated integration of genome- and metabolome-scale network reconstruction and data-driven inverse modeling of metabolic interaction networks.

One central goal of systems biology is to infer biochemical regulations from large-scale OMICS data. Many aspects of cellular physiology and organismal phenotypes can be understood as results of metabolic interaction network dynamics. Previously, we have proposed a convenient mathematical method, which addresses this problem using metabolomics data for the inverse calculation of biochemical Jacobian matrices revealing regulatory checkpoints of biochemical regulations. The proposed algorithms for this inference are limited by two issues: they rely on structural network information that needs to be assembled manually, and they are numerically unstable due to ill-conditioned regression problems for large-scale metabolic networks. To address these problems, we developed a novel regression loss-based inverse Jacobian algorithm, combining metabolomics COVariance and genome-scale metabolic RECONstruction, which allows for a fully automated, algorithmic implementation of the COVRECON workflow. It consists of two parts: (i) Sim-Network and (ii) inverse differential Jacobian evaluation. Sim-Network automatically generates an organism-specific enzyme and reaction dataset from Bigg and KEGG databases, which is then used to reconstruct the Jacobian's structure for a specific metabolomics dataset. Instead of directly solving a regression problem as in the previous workflow, the new inverse differential Jacobian is based on a substantially more robust approach and rates the biochemical interactions according to their relevance from large-scale metabolomics data. The approach is illustrated by in silico stochastic analysis with differently sized metabolic networks from the BioModels database and applied to a real-world example. The characteristics of the COVRECON implementation are that (i) it automatically reconstructs a data-driven superpathway model; (ii) more general network structures can be investigated, and (iii) the new inverse algorithm improves stability, decreases computation time, and extends to large-scale models. The code is available in the website https://bitbucket.org/mosys-univie/covrecon.

Distinct metabolic features of genetic liability to type 2 diabetes and coronary artery disease: a reverse Mendelian randomization study.

Type 2 diabetes (T2D) and coronary artery disease (CAD) both have known genetic determinants, but the mechanisms through which their associated genetic variants lead to disease onset remain poorly understood. We used large-scale metabolomics data in a two-sample reverse Mendelian randomization (MR) framework to estimate effects of genetic liability to T2D and CAD on 249 circulating metabolites in the UK Biobank (N=118,466). We examined the potential for medication use to distort effect estimates by conducting age-stratified metabolite analyses. Using inverse variance weighted (IVW) models, higher genetic liability to T2D was estimated to decrease high-density lipoprotein cholesterol (HDL-C) and low-density lipoprotein cholesterol (LDL-C) (e.g. -0.05 SD; 95% CI-0.07 to-0.03, per doubling of liability), whilst increasing all triglyceride groups and branched chain amino acids (BCAAs). IVW estimates for CAD liability suggested an effect on reducing HDL-C as well as raising very-low density lipoprotein cholesterol (VLDL-C) and LDL-C. In pleiotropy-robust models, T2D liability was still estimated to increase BCAAs, but several estimates for higher CAD liability reversed and supported decreased LDL-C and apolipoprotein-B. Estimated effects of CAD liability differed substantially by age for non-HDL-C traits, with higher CAD liability lowering LDL-C only at older ages when statin use was common. Overall, our results support largely distinct metabolic features of genetic liability to T2D and CAD, illustrating both challenges and opportunities for preventing these commonly co-occurring diseases. Wellcome Trust [218495/Z/19/Z], UK MRC [MC_UU_00011/1; MC_UU_00011/4], the University of Bristol, Diabetes UK [17/0005587], World Cancer Research Fund [IIG_2019_2009].

Norm ISWSVR: A Data Integration and Normalization Approach for Large-Scale Metabolomics.

Large-scale and long-period metabolomics study is more susceptible to various sources of systematic errors, resulting in nonreproducibility and poor data quality. A reliable and robust batch correction method removes unwanted systematic variations and improves the statistical power of metabolomics data, which undeniably becomes an important issue for the quality control of metabolomics. This study proposed a novel data normalization and integration method, Norm ISWSVR. It is a two-step approach via combining the best-performance internal standard correction with support vector regression normalization, comprehensively removing the systematic and random errors and matrix effects. This method was investigated in three untargeted lipidomics or metabolomics datasets, and the performance was further evaluated systematically in comparison with that of 11 other normalization methods. As a result, Norm ISWSVR decreased the data's median cross-validated relative standard deviation (cvRSD), increased the correlation between QCs, improved the classification accuracy of biomarkers, and was well-compatible with quantitative data. More importantly, Norm ISWSVR also allows a low frequency of QCs, which could significantly decrease the burden of a large-scale experiment. Correspondingly, Norm ISWSVR favorably improves the data quality of large-scale metabolomics data.

RMisbeta: A robust missing value imputation approach in transcriptomics and metabolomics data

Transcriptomics and metabolomics data often contain missing values or outliers due to limitations of the data acquisition techniques. Most of the statistical methods require complete datasets for downstream analysis. A number of methods have been developed for missing value imputation using the classical mean and variance based on maximum likelihood estimators, which are not robust against outliers. Consequently, the performance of these methods deteriorates in the presence of outliers. Hence precise imputation of missing values and outliers handling are both concurrently important. Therefore, in this paper, we developed a robust iterative approach using robust estimators based on the minimum beta divergence method, which simultaneously impute missing values and outliers. We investigate the performance of the proposed method in a comparison with six frequently used missing value imputation methods such as Zero, KNN, robust SVD, EM, random forest (RF) and weighted least square approach (WLSA) through feature selection using both simulated and real datasets. Ten performance indices were used to explore the optimal method such as Frobenius norm (FOBN), accuracy (ACC), sensitivity (SN), specificity (SP), positive predictive value (PPV), negative predictive value (NPV), detection rate (DR), misclassification error rate (MER), the area under the ROC curve (AUC) and computational runtime. Evaluation based on both simulated and real data suggests the superiority of the proposed method over the other traditional methods in terms of various rates of outliers and missing values. The suggested approach also keeps almost equal performance in absence of outliers with the other methods. The proposed method is accurate, simple, and consumes lower computational time compared to the other methods. Therefore, our recommendation is to apply the proposed procedure for large-scale transcriptomics and metabolomics data analysis. The computational tool has been implemented in an R package, which is publicly available from https://CRAN.R-project.org/package=rMisbeta.

Metabolomics-Based Screening of Inborn Errors of Metabolism: Enhancing Clinical Application with a Robust Computational Pipeline.

Inborn errors of metabolism (IEM) are inherited conditions caused by genetic defects in enzymes or cofactors. These defects result in a specific metabolic fingerprint in patient body fluids, showing accumulation of substrate or lack of an end-product of the defective enzymatic step. Untargeted metabolomics has evolved as a high throughput methodology offering a comprehensive readout of this metabolic fingerprint. This makes it a promising tool for diagnostic screening of IEM patients. However, the size and complexity of metabolomics data have posed a challenge in translating this avalanche of information into knowledge, particularly for clinical application. We have previously established next-generation metabolic screening (NGMS) as a metabolomics-based diagnostic tool for analyzing plasma of individual IEM-suspected patients. To fully exploit the clinical potential of NGMS, we present a computational pipeline to streamline the analysis of untargeted metabolomics data. This pipeline allows for time-efficient and reproducible data analysis, compatible with ISO:15189 accredited clinical diagnostics. The pipeline implements a combination of tools embedded in a workflow environment for large-scale clinical metabolomics data analysis. The accompanying graphical user interface aids end-users from a diagnostic laboratory for efficient data interpretation and reporting. We also demonstrate the application of this pipeline with a case study and discuss future prospects.

A hierarchical approach to removal of unwanted variation for large-scale metabolomics data

Liquid chromatography-mass spectrometry-based metabolomics studies are increasingly applied to large population cohorts, which run for several weeks or even years in data acquisition. This inevitably introduces unwanted intra- and inter-batch variations over time that can overshadow true biological signals and thus hinder potential biological discoveries. To date, normalisation approaches have struggled to mitigate the variability introduced by technical factors whilst preserving biological variance, especially for protracted acquisitions. Here, we propose a study design framework with an arrangement for embedding biological sample replicates to quantify variance within and between batches and a workflow that uses these replicates to remove unwanted variation in a hierarchical manner (hRUV). We use this design to produce a dataset of more than 1000 human plasma samples run over an extended period of time. We demonstrate significant improvement of hRUV over existing methods in preserving biological signals whilst removing unwanted variation for large scale metabolomics studies. Our tools not only provide a strategy for large scale data normalisation, but also provides guidance on the design strategy for large omics studies.

Reference Standardization for Quantification and Harmonization of Large-Scale Metabolomics.

Reference standardization was developed to address quantification and harmonization challenges for high-resolution metabolomics (HRM) data collected across different studies or analytical methods. Reference standardization relies on the concurrent analysis of calibrated pooled reference samples at predefined intervals and enables a single-step batch correction and quantification for high-throughput metabolomics. Here, we provide quantitative measures of approximately 200 metabolites for each of three pooled reference materials (220 metabolites for Qstd3, 211 metabolites for CHEAR, 204 metabolites for NIST1950) and show application of this approach for quantification supports harmonization of metabolomics data collected from 3677 human samples in 17 separate studies analyzed by two complementary HRM methods over a 17-month period. The results establish reference standardization as a method suitable for harmonizing large-scale metabolomics data and extending capabilities to quantify large numbers of known and unidentified metabolites detected by high-resolution mass spectrometry methods.

A novel analysis method for biomarker identification based on horizontal relationship: identifying potential biomarkers from large-scale hepatocellular carcinoma metabolomics data.

Omics techniques develop quickly and have made a great contribution to disease study. Omics data are usually complex. How to analyze the data and mine important information has been a key part in omics research. To study the nature of disease mechanisms systematically, we propose a new data analysis method to define the network biomarkers based on horizontal comparison (DNB-HC). DNB-HC performs molecule horizontal relationships to characterize the physiological status and differential network analysis to screen the biomarkers. We applied DNB-HC to analyze a large-scale metabolomics data, which contained 550 samples from a nested case-control hepatocellular carcinoma (HCC) study. A network biomarker was defined, and its areas under curves (AUC) in the receiver-operating characteristic (ROC) analysis for HCC discrimination were larger than those defined by six efficient feature selection methods in most cases. The effectiveness was further corroborated by another nested HCC dataset. Besides, the performance of the defined biomarkers was better than that of α-fetoprotein (AFP), a commonly used clinical biomarker for distinguishing HCC from high-risk population of liver cirrhosis in other two independent metabolomics validation sets. All and 90.3% of the AFP false-negative patients with HCC were correctly diagnosed in these two sets, respectively. The experimental results illustrate that DNB-HC can mine more important information reflecting the nature of the research problems by studying the feature horizontal relationship systematically and identifying effective disease biomarkers in clinical practice. Graphical abstract.

Statistical Workflow for Feature Selection in Human Metabolomics Data.

High-throughput metabolomics investigations, when conducted in large human cohorts, represent a potentially powerful tool for elucidating the biochemical diversity underlying human health and disease. Large-scale metabolomics data sources, generated using either targeted or nontargeted platforms, are becoming more common. Appropriate statistical analysis of these complex high-dimensional data will be critical for extracting meaningful results from such large-scale human metabolomics studies. Therefore, we consider the statistical analytical approaches that have been employed in prior human metabolomics studies. Based on the lessons learned and collective experience to date in the field, we offer a step-by-step framework for pursuing statistical analyses of cohort-based human metabolomics data, with a focus on feature selection. We discuss the range of options and approaches that may be employed at each stage of data management, analysis, and interpretation and offer guidance on the analytical decisions that need to be considered over the course of implementing a data analysis workflow. Certain pervasive analytical challenges facing the field warrant ongoing focused research. Addressing these challenges, particularly those related to analyzing human metabolomics data, will allow for more standardization of as well as advances in how research in the field is practiced. In turn, such major analytical advances will lead to substantial improvements in the overall contributions of human metabolomics investigations.

WaveICA: A novel algorithm to remove batch effects for large-scale untargeted metabolomics data based on wavelet analysis

Metabolomics provides new insights into disease pathogenesis and biomarker discovery. Samples from large-scale untargeted metabolomics studies are typically analyzed using a liquid chromatography-mass spectrometry platform in several batches. Batch effects that are caused by non-biological systematic biases are unavoidable in large-scale metabolomics studies, even with properly designed experiments. The statistical analysis of large-scale metabolomics data without managing batch effects will yield misleading results. In this study, we propose a novel algorithm, called WaveICA, which is based on the wavelet transform method with independent component analysis, as the threshold processing method to capture and remove batch effects for large-scale metabolomics data. The WaveICA method uses the time trend of samples over the injection order, decomposes the original data into multi-scale data with different features, extracts and removes the batch effect information in multi-scale data, and obtains clean data. The WaveICA method was tested on real metabolomics data. After applying the WaveICA method, scattered quality control samples (QCS) and subject samples in a PCA score plot of the original data were closely clustered, respectively. The average Pearson correlation coefficients for all peaks of the QCS increased from 0.872 to 0.972. Additionally, WaveICA significantly improved the classification accuracy for metabolomics data. The method was compared with three representative methods, and outperformed all of them. To conclude, WaveICA can efficiently remove batch effects while revealing more biological information. This method can be used in large-scale untargeted metabolomics studies to preprocess raw metabolomics data.

Metabolomics in epidemiologic research: challenges and opportunities for early-career epidemiologists.

The application of metabolomics to epidemiologic studies is increasing. Here, we describe the challenges and opportunities facing early-career epidemiologists aiming to apply metabolomics to their research. Many challenges inherent to metabolomics may provide early-career epidemiologists with the opportunity to play a pivotal role in answering critical methodological questions and moving the field forward. Although generating large-scale high-quality metabolomics data can be challenging, data can be accessed through public databases, collaboration with senior researchers or participation within interest groups. Such efforts may also assist with obtaining funding, provide knowledge on training resources, and help early-career epidemiologists to publish in the field of metabolomics.

Metabolic network segmentation: A probabilistic graphical modeling approach to identify the sites and sequential order of metabolic regulation from non-targeted metabolomics data.

In recent years, the number of large-scale metabolomics studies on various cellular processes in different organisms has increased drastically. However, it remains a major challenge to perform a systematic identification of mechanistic regulatory events that mediate the observed changes in metabolite levels, due to complex interdependencies within metabolic networks. We present the metabolic network segmentation (MNS) algorithm, a probabilistic graphical modeling approach that enables genome-scale, automated prediction of regulated metabolic reactions from differential or serial metabolomics data. The algorithm sections the metabolic network into modules of metabolites with consistent changes. Metabolic reactions that connect different modules are the most likely sites of metabolic regulation. In contrast to most state-of-the-art methods, the MNS algorithm is independent of arbitrary pathway definitions, and its probabilistic nature facilitates assessments of noisy and incomplete measurements. With serial (i.e., time-resolved) data, the MNS algorithm also indicates the sequential order of metabolic regulation. We demonstrated the power and flexibility of the MNS algorithm with three, realistic case studies with bacterial and human cells. Thus, this approach enables the identification of mechanistic regulatory events from large-scale metabolomics data, and contributes to the understanding of metabolic processes and their interplay with cellular signaling and regulation processes.

Potential Impact and Study Considerations of Metabolomics in Cardiovascular Health and Disease: A Scientific Statement From the American Heart Association.

Through the measure of thousands of small-molecule metabolites in diverse biological systems, metabolomics now offers the potential for new insights into the factors that contribute to complex human diseases such as cardiovascular disease. Targeted metabolomics methods have already identified new molecular markers and metabolomic signatures of cardiovascular disease risk (including branched-chain amino acids, select unsaturated lipid species, and trimethylamine-N-oxide), thus in effect linking diverse exposures such as those from dietary intake and the microbiota with cardiometabolic traits. As technologies for metabolomics continue to evolve, the depth and breadth of small-molecule metabolite profiling in complex systems continue to advance rapidly, along with prospects for ongoing discovery. Current challenges facing the field of metabolomics include scaling throughput and technical capacity for metabolomics approaches, bioinformatic and chemoinformatic tools for handling large-scale metabolomics data, methods for elucidating the biochemical structure and function of novel metabolites, and strategies for determining the true clinical relevance of metabolites observed in association with cardiovascular disease outcomes. Progress made in addressing these challenges will allow metabolomics the potential to substantially affect diagnostics and therapeutics in cardiovascular medicine.

Sparse network modeling and metscape-based visualization methods for the analysis of large-scale metabolomics data.

Recent technological advances in mass spectrometry, development of richer mass spectral libraries and data processing tools have enabled large scale metabolic profiling. Biological interpretation of metabolomics studies heavily relies on knowledge-based tools that contain information about metabolic pathways. Incomplete coverage of different areas of metabolism and lack of information about non-canonical connections between metabolites limits the scope of applications of such tools. Furthermore, the presence of a large number of unknown features, which cannot be readily identified, but nonetheless can represent bona fide compounds, also considerably complicates biological interpretation of the data. Leveraging recent developments in the statistical analysis of high-dimensional data, we developed a new Debiased Sparse Partial Correlation algorithm (DSPC) for estimating partial correlation networks and implemented it as a Java-based CorrelationCalculator program. We also introduce a new version of our previously developed tool Metscape that enables building and visualization of correlation networks. We demonstrate the utility of these tools by constructing biologically relevant networks and in aiding identification of unknown compounds. http://metscape.med.umich.edu. Supplementary data are available at Bioinformatics online.

Normalization and integration of large-scale metabolomics data using support vector regression

Untargeted metabolomics studies for biomarker discovery often have hundreds to thousands of human samples. Data acquisition of large-scale samples has to be divided into several batches and may span from months to as long as several years. The signal drift of metabolites during data acquisition (intra- and inter-batch) is unavoidable and is a major confounding factor for large-scale metabolomics studies. We aim to develop a data normalization method to reduce unwanted variations and integrate multiple batches in large-scale metabolomics studies prior to statistical analyses. We developed a machine learning algorithm-based method, support vector regression (SVR), for large-scale metabolomics data normalization and integration. An R package named MetNormalizer was developed and provided for data processing using SVR normalization. After SVR normalization, the portion of metabolite ion peaks with relative standard deviations (RSDs) less than 30 % increased to more than 90 % of the total peaks, which is much better than other common normalization methods. The reduction of unwanted analytical variations helps to improve the performance of multivariate statistical analyses, both unsupervised and supervised, in terms of classification and prediction accuracy so that subtle metabolic changes in epidemiological studies can be detected. SVR normalization can effectively remove the unwanted intra- and inter-batch variations, and is much better than other common normalization methods.

Embracing new-generation 'omics' tools to improve drought tolerance in cereal and food-legume crops

Drought stress presents a considerable threat to the global crop production. As a dominant source of vegetarian diet, cereals and grain-legumes remain crucial to meeting the growing dietary demands worldwide. Therefore, breeding cultivars of these staple crops with enhanced drought tolerance stands to be one of the most sustainable solutions to enhance food production in changing climate. Given the context, a more focused survey of environment-defined germplasm sets is imperative to comprehend such adaptive traits. In parallel, uncovering the genetic architecture and the molecular networks that collectively contribute towards drought tolerance is urgently required through rationally combining large-scale genomics, proteomics, and metabolomics data. Also, attention needs to be directed to reasonably quantify the epistatic as well as environmental influences, thereby warranting deployment of analyses like metaquantitative trait loci (QTL) that encompass multiple environments and diverse genetic backgrounds. Further, innovative techniques like genomic selection (GS) and genome wide association study (GWAS) would help to capture the quantitative variation underlying drought tolerance. Equally importantly, integration of physiological traits-based techniques with ever-evolving ‘omics’ technologies and the new-generation phenotyping platforms will be of immense importance in advancing our existing knowledge about the genetically-complex and poorly-understood phenomena, such as plant drought response, and a deeper understanding would likely to provide a great impetus to the progress of crop breeding for drought tolerance.

Bayesian Independent Component Analysis Recovers Pathway Signatures from Blood Metabolomics Data

Interpreting the complex interplay of metabolites in heterogeneous biosamples still poses a challenging task. In this study, we propose independent component analysis (ICA) as a multivariate analysis tool for the interpretation of large-scale metabolomics data. In particular, we employ a Bayesian ICA method based on a mean-field approach, which allows us to statistically infer the number of independent components to be reconstructed. The advantage of ICA over correlation-based methods like principal component analysis (PCA) is the utilization of higher order statistical dependencies, which not only yield additional information but also allow a more meaningful representation of the data with fewer components. We performed the described ICA approach on a large-scale metabolomics data set of human serum samples, comprising a total of 1764 study probands with 218 measured metabolites. Inspecting the source matrix of statistically independent metabolite profiles using a weighted enrichment algorithm, we observe strong enrichment of specific metabolic pathways in all components. This includes signatures from amino acid metabolism, energy-related processes, carbohydrate metabolism, and lipid metabolism. Our results imply that the human blood metabolome is composed of a distinct set of overlaying, statistically independent signals. ICA furthermore produces a mixing matrix, describing the strength of each independent component for each of the study probands. Correlating these values with plasma high-density lipoprotein (HDL) levels, we establish a novel association between HDL plasma levels and the branched-chain amino acid pathway. We conclude that the Bayesian ICA methodology has the power and flexibility to replace many of the nowadays common PCA and clustering-based analyses common in the research field.