Amanida: an R package for meta-analysis of metabolomics non-integral data

MetaHD: A multivariate meta-analysis model for metabolomics data

MOTIVATION

Meta-analysis methods widely-used for combining metabolomics data do not account for correlation between metabolites or missing values. Within- and between-study variability are also often overlooked. These can give results with inferior statistical properties, leading to misidentification of biomarkers.

RESULTS

We propose a multivariate meta-analysis model for high-dimensional metabolomics data (MetaHD), which accommodates the correlation between metabolites, within- and between-study variances, and missing values. MetaHD can be used for integrating and collectively analysing individual-level metabolomics data generated from multiple studies as well as for combining summary estimates. We show that MetaHD leads to lower root mean square error compared to the existing approaches. Furthermore, we demonstrate that MetaHD, which exploits the borrowing strength between metabolites, could be particularly useful in the presence of missing data compared to univariate meta-analysis methods, which can return biased estimates in the presence of data missing at random.

AVAILABILITY AND IMPLEMENTATION

The MetaHD R package can be downloaded through Comprehensive R Archive Network (CRAN) repository. A detailed vignette with example datasets and code to prepare data and analyses are available on https://bookdown.org/a2delivera/MetaHD/.

SUPPLEMENTARY INFORMATION

Supplementary data are available at Bioinformatics online.

Metabolomic profiling of upper GI malignancies in blood and tissue: a systematic review and meta-analysis

OBJECTIVE

To conduct a systematic review and meta-analysis of case-control and cohort human studies evaluating metabolite markers identified using high-throughput metabolomics techniques on esophageal cancer (EC), cancer of the gastroesophageal junction (GEJ), and gastric cancer (GC) in blood and tissue.

BACKGROUND

Upper gastrointestinal cancers (UGC), predominantly EC, GEJ, and GC, are malignant tumour types with high morbidity and mortality rates. Numerous studies have focused on metabolomic profiling of UGC in recent years. In this systematic review and meta-analysis, we have provided a collective summary of previous findings on metabolites and metabolomic profiling associated with EC, GEJ and GC.

METHODS

Following the PRISMA procedure, a systematic search of four databases (Embase, PubMed, MEDLINE, and Web of Science) for molecular epidemiologic studies on the metabolomic profiles of EC, GEJ and GC was conducted and registered at PROSPERO (CRD42023486631). The Newcastle-Ottawa Scale (NOS) was used to benchmark the risk of bias for case-controlled and cohort studies. QUADOMICS, an adaptation of the QUADAS-2 (Quality Assessment of Diagnostic Accuracy) tool, was used to rate diagnostic accuracy studies. Original articles comparing metabolite patterns between patients with and without UGC were included. Two investigators independently completed title and abstract screening, data extraction, and quality evaluation. Meta-analysis was conducted whenever possible. We used a random effects model to investigate the association between metabolite levels and UGC.

RESULTS

A total of 66 original studies involving 7267 patients that met the required criteria were included for review. 169 metabolites were differentially distributed in patients with UGC compared to healthy patients among 44 GC, 9 GEJ, and 25 EC studies including metabolites involved in glycolysis, anaerobic respiration, tricarboxylic acid cycle, and lipid metabolism. Phosphatidylcholines, eicosanoids, and adenosine triphosphate were among the most frequently reported lipids and metabolites of cellular respiration, while BCAA, lysine, and asparagine were among the most commonly reported amino acids. Previously identified lipid metabolites included saturated and unsaturated free fatty acids and ketones. However, the key findings across studies have been inconsistent, possibly due to limited sample sizes and the majority being hospital-based case-control analyses lacking an independent replication group.

CONCLUSION

Thus far, metabolomic studies have provided new opportunities for screening, etiological factors, and biomarkers for UGC, supporting the potential of applying metabolomic profiling in early cancer diagnosis. According to the results of our meta-analysis especially BCAA and TMAO as well as certain phosphatidylcholines should be implicated into the diagnostic procedure of patients with UGC. We envision that metabolomics will significantly enhance our understanding of the carcinogenesis and progression process of UGC and may eventually facilitate precise oncological and patient-tailored management of UGC.

Easy-Amanida: An R Shiny application for the meta-analysis of aggregate results in clinical metabolomics using Amanida and Webchem

Meta-analysis is a useful tool in clinical research, as it combines the results of multiple clinical studies to improve precision when answering a particular scientific question. While there has been a substantial increase in publications using meta-analysis in various clinical research topics, the number of published meta-analyses in metabolomics is significantly lower compared to other omics disciplines. Metabolomics is the study of small chemical compounds in living organisms, which provides important insights into an organism's phenotype. However, the wide variety of compounds and the different experimental methods used in metabolomics make it challenging to perform a thorough meta-analysis. Additionally, there is a lack of consensus on reporting statistical estimates, and the high number of compound naming synonyms further complicates the process. Easy-Amanida is a new tool that combines two R packages, "amanida" and "webchem", to enable meta-analysis of aggregate statistical data, like p-value and fold-change, while ensuring the compounds naming harmonization. The Easy-Amanida app is implemented in Shiny, an R package add-on for interactive web apps, and provides a workflow to optimize the naming combination. This article describes all the steps to perform the meta-analysis using Easy-Amanida, including an illustrative example for interpreting the results. The use of aggregate statistics metrics extends the use of Easy-Amanida beyond the metabolomics field.

Metabolomic data presents challenges for epidemiological meta-analysis: a case study of childhood body mass index from the ECHO consortium

INTRODUCTION

Meta-analyses across diverse independent studies provide improved confidence in results. However, within the context of metabolomic epidemiology, meta-analysis investigations are complicated by differences in study design, data acquisition, and other factors that may impact reproducibility.

OBJECTIVE

The objective of this study was to identify maternal blood metabolites during pregnancy (> 24 gestational weeks) related to offspring body mass index (BMI) at age two years through a meta-analysis framework.

METHODS

We used adjusted linear regression summary statistics from three cohorts (total N = 1012 mother-child pairs) participating in the NIH Environmental influences on Child Health Outcomes (ECHO) Program. We applied a random-effects meta-analysis framework to regression results and adjusted by false discovery rate (FDR) using the Benjamini-Hochberg procedure.

RESULTS

Only 20 metabolites were detected in all three cohorts, with an additional 127 metabolites detected in two of three cohorts. Of these 147, 6 maternal metabolites were nominally associated (P < 0.05) with offspring BMI z-scores at age 2 years in a meta-analytic framework including at least two studies: arabinose (Coefmeta = 0.40 [95% CI 0.10,0.70], Pmeta = 9.7 × 10-3), guanidinoacetate (Coefmeta = - 0.28 [- 0.54, - 0.02], Pmeta = 0.033), 3-ureidopropionate (Coefmeta = 0.22 [0.017,0.41], Pmeta = 0.033), 1-methylhistidine (Coefmeta = - 0.18 [- 0.33, - 0.04], Pmeta = 0.011), serine (Coefmeta = - 0.18 [- 0.36, - 0.01], Pmeta = 0.034), and lysine (Coefmeta = - 0.16 [- 0.32, - 0.01], Pmeta = 0.044). No associations were robust to multiple testing correction.

CONCLUSIONS

Despite including three cohorts with large sample sizes (N > 100), we failed to identify significant metabolite associations after FDR correction. Our investigation demonstrates difficulties in applying epidemiological meta-analysis to clinical metabolomics, emphasizes challenges to reproducibility, and highlights the need for standardized best practices in metabolomic epidemiology.

Deciphering Urogenital Cancers through Proteomic Biomarkers: A Systematic Review and Meta-Analysis

Urogenital cancers, which include prostate, bladder, and kidney malignancies, exert a substantial impact on global cancer-related morbidity and mortality. Proteomic biomarkers, emerging as valuable tools, aim to enhance early detection, prognostic accuracy, and the development of personalized therapeutic strategies. This study undertook a comprehensive systematic review and meta-analysis of the existing literature investigating the role and potential of proteomic biomarkers in plasma, tissue, and urine samples in urogenital cancers. Our extensive search across several databases identified 1879 differentially expressed proteins from 37 studies, signifying their potential as unique biomarkers for these cancers. A meta-analysis of the significantly differentially expressed proteins was executed, accentuating the findings through visually intuitive volcano plots. A functional enrichment analysis unveiled their significant involvement in diverse biological processes, including signal transduction, immune response, cell communication, and cell growth. A pathway analysis highlighted the participation of key pathways such as the nectin adhesion pathway, TRAIL signaling pathway, and integrin signaling pathways. These findings not only pave the way for future investigations into early detection and targeted therapeutic approaches but also underscore the fundamental role of proteomics in advancing our understanding of the molecular mechanisms underpinning urogenital cancer pathogenesis. Ultimately, these findings hold remarkable potential to significantly enhance patient care and improve clinical outcomes.

Untargeted serum metabolomics reveals novel metabolite associations and disruptions in amino acid and lipid metabolism in Parkinson's disease

BACKGROUND

Untargeted high-resolution metabolomic profiling provides simultaneous measurement of thousands of metabolites. Metabolic networks based on these data can help uncover disease-related perturbations across interconnected pathways.

OBJECTIVE

Identify metabolic disturbances associated with Parkinson's disease (PD) in two population-based studies using untargeted metabolomics.

METHODS

We performed a metabolome-wide association study (MWAS) of PD using serum-based untargeted metabolomics data derived from liquid chromatography with high-resolution mass spectrometry (LC-HRMS) using two distinct population-based case-control populations. We also combined our results with a previous publication of 34 metabolites linked to PD in a large-scale, untargeted MWAS to assess external validation.

RESULTS

LC-HRMS detected 4,762 metabolites for analysis (HILIC: 2716 metabolites; C18: 2046 metabolites). We identified 296 features associated with PD at FDR<0.05, 134 having a log2 fold change (FC) beyond ±0.5 (228 beyond ±0.25). Of these, 104 were independently associated with PD in both discovery and replication studies at p<0.05 (170 at p<0.10), while 27 were associated with levodopa-equivalent dose among the PD patients. Intriguingly, among the externally validated features were the microbial-related metabolites, p-cresol glucuronide (FC=2.52, 95% CI=1.67, 3.81, FDR=7.8e-04) and p-cresol sulfate. P-cresol glucuronide was also associated with motor symptoms among patients. Additional externally validated metabolites associated with PD include phenylacetyl-L-glutamine, trigonelline, kynurenine, biliverdin, and pantothenic acid. Novel associations include the anti-inflammatory metabolite itaconate (FC=0.79, 95% CI=0.73, 0.86; FDR=2.17E-06) and cysteine-S-sulfate (FC=1.56, 95% CI=1.39, 1.75; FDR=3.43E-11). Seventeen pathways were enriched, including several related to amino acid and lipid metabolism.

CONCLUSIONS

Our results revealed PD-associated metabolites, confirming several previous observations, including for p-cresol glucuronide, and newly implicating interesting metabolites, such as itaconate. Our data also suggests metabolic disturbances in amino acid and lipid metabolism and inflammatory processes in PD.

Metabolome panels as potential noninvasive biomarkers for primary glomerulonephritis sub-types: meta-analysis of profiling metabolomics studies

Primary glomerulonephritis diseases (PGDs) are known as the top causes of chronic kidney disease worldwide. Renal biopsy, an invasive method, is the main approach to diagnose PGDs. Studying the metabolome profiles of kidney diseases is an inclusive approach to identify the disease's underlying pathways and discover novel non-invasive biomarkers. So far, different experiments have explored the metabolome profiles in different PGDs, but the inconsistencies might hinder their clinical translations. The main goal of this meta-analysis study was to achieve consensus panels of dysregulated metabolites in PGD sub-types. The PGDs-related metabolome profiles from urine samples in humans were selected in a comprehensive search. Amanida package in R software was utilized for performing the meta-analysis. Through sub-type analyses, the consensus list of metabolites in each category was obtained. To identify the most affected pathways, functional enrichment analysis was performed. Also, a gene-metabolite network was constructed to identify the key metabolites and their connected proteins. After a vigorous search, among the 11 selected studies (15 metabolite profiles), 270 dysregulated metabolites were recognized in urine of 1154 PGDs and control samples. Through sub-type analyses by Amanida package, the consensus list of metabolites in each category was obtained. Top dysregulated metabolites (vote score of ≥ 4 or ≤ - 4) in PGDs urines were selected as main panel of meta-metabolites including glucose, leucine, choline, betaine, dimethylamine, fumaric acid, citric acid, 3-hydroxyisovaleric acid, pyruvic acid, isobutyric acid, and hippuric acid. The enrichment analyses results revealed the involvement of different biological pathways such as the TCA cycle and amino acid metabolisms in the pathogenesis of PGDs. The constructed metabolite-gene interaction network revealed the high centralities of several metabolites, including pyruvic acid, leucine, and choline. The identified metabolite panels could shed a light on the underlying pathological pathways and be considered as non-invasive biomarkers for the diagnosis of PGD sub-types.

Molecular and Metabolic Phenotyping of Hepatocellular Carcinoma for Biomarker Discovery: A Meta-Analysis

Identifying and translating hepatocellular carcinoma (HCC) biomarkers from bench to bedside using mass spectrometry-based metabolomics and lipidomics is hampered by inconsistent findings. Here, we investigated HCC at systemic and metabolism-centric multiomics levels by conducting a meta-analysis of quantitative evidence from 68 cohorts. Blood transcript biomarkers linked to the HCC metabolic phenotype were externally validated and prioritized. In the studies under investigation, about 600 metabolites were reported as putative HCC-associated biomarkers; 39, 20, and 10 metabolites and 52, 12, and 12 lipids were reported in three or more studies in HCC vs. Control, HCC vs. liver cirrhosis (LC), and LC vs. Control groups, respectively. Amino acids, fatty acids (increased 18:1), bile acids, and lysophosphatidylcholine were the most frequently reported biomarkers in HCC. BAX and RAC1 showed a good correlation and were associated with poor prognosis. Our study proposes robust HCC biomarkers across diverse cohorts using a data-driven knowledge-based approach that is versatile and affordable for studying other diseases.

The BinDiscover database: a biology-focused meta-analysis tool for 156,000 GC-TOF MS metabolome samples

Metabolomics by gas chromatography/mass spectrometry (GC/MS) provides a standardized and reliable platform for understanding small molecule biology. Since 2005, the West Coast Metabolomics Center at the University of California at Davis has collated GC/MS metabolomics data from over 156,000 samples and 2000 studies into the standardized BinBase database. We believe that the observations from these samples will provide meaningful insight to biologists and that our data treatment and webtool will provide insight to others who seek to standardize disparate metabolomics studies. We here developed an easy-to-use query interface, BinDiscover, to enable intuitive, rapid hypothesis generation for biologists based on these metabolomic samples. BinDiscover creates observation summaries and graphics across a broad range of species, organs, diseases, and compounds. Throughout the components of BinDiscover, we emphasize the use of ontologies to aggregate large groups of samples based on the proximity of their metadata within these ontologies. This adjacency allows for the simultaneous exploration of entire categories such as "rodents", "digestive tract", or "amino acids". The ontologies are particularly relevant for BinDiscover's ontologically grouped differential analysis, which, like other components of BinDiscover, creates clear graphs and summary statistics across compounds and biological metadata. We exemplify BinDiscover's extensive applicability in three showcases across biological domains.

Alignment and Analysis of a Disparately Acquired Multibatch Metabolomics Study of Maternal Pregnancy Samples

Untargeted liquid chromatography-mass spectrometry metabolomics studies are typically performed under roughly identical experimental settings. Measurements acquired with different LC-MS protocols or following extended time intervals harbor significant variation in retention times and spectral abundances due to altered chromatographic, spectrometric, and other factors, raising many data analysis challenges. We developed a computational workflow for merging and harmonizing metabolomics data acquired under disparate LC-MS conditions. Plasma metabolite profiles were collected from two sets of maternal subjects three years apart using distinct instruments and LC-MS procedures. Metabolomics features were aligned using metabCombiner to generate lists of compounds detected across all experimental batches. We applied data set-specific normalization methods to remove interbatch and interexperimental variation in spectral intensities, enabling statistical analysis on the assembled data matrix. Bioinformatics analyses revealed large-scale metabolic changes in maternal plasma between the first and third trimesters of pregnancy and between maternal plasma and umbilical cord blood. We observed increases in steroid hormones and free fatty acids from the first trimester to term of gestation, along with decreases in amino acids coupled to increased levels in cord blood. This work demonstrates the viability of integrating nonidentically acquired LC-MS metabolomics data and its utility in unconventional metabolomics study designs.

iDMET: network-based approach for integrating differential analysis of cancer metabolomics

BACKGROUND

Comprehensive metabolomic analyses have been conducted in various institutes and a large amount of metabolomic data are now publicly available. To help fully exploit such data and facilitate their interpretation, metabolomic data obtained from different facilities and different samples should be integrated and compared. However, large-scale integration of such data for biological discovery is challenging given that they are obtained from various types of sample at different facilities and by different measurement techniques, and the target metabolites and sensitivities to detect them also differ from study to study.

RESULTS

We developed iDMET, a network-based approach to integrate metabolomic data from different studies based on the differential metabolomic profiles between two groups, instead of the metabolite profiles themselves. As an application, we collected cancer metabolomic data from 27 previously published studies and integrated them using iDMET. A pair of metabolomic changes observed in the same disease from two studies were successfully connected in the network, and a new association between two drugs that may have similar effects on the metabolic reactions was discovered.

CONCLUSIONS

We believe that iDMET is an efficient tool for integrating heterogeneous metabolomic data and discovering novel relationships between biological phenomena.

Urine NMR Metabolomics for Precision Oncology in Colorectal Cancer

Metabolomics is a fundamental approach to discovering novel biomarkers and their potential use for precision medicine. When applied for population screening, NMR-based metabolomics can become a powerful clinical tool in precision oncology. Urine tests can be more widely accepted due to their intrinsic non-invasiveness. Our review provides the first exhaustive evaluation of NMR metabolomics for the determination of colorectal cancer (CRC) in urine. A specific search in PubMed, Web of Science, and Scopus was performed, and 10 studies met the required criteria. There were no restrictions on the query for study type, leading to not only colorectal cancer samples versus control comparisons, but also prospective studies of surgical effects. With this review, all compounds in the included studies were merged into a database. In doing so, we identified up to 100 compounds in urine samples, and 11 were found in at least three articles. Results were analyzed in three groups: case (CRC and adenomas)/control, pre-/post-surgery, and combining both groups. When combining the case-control and the pre-/post-surgery groups, up to twelve compounds were found to be relevant. Seven down-regulated metabolites in CRC were identified, creatinine, 4-hydroxybenzoic acid, acetone, carnitine, d-glucose, hippuric acid, l-lysine, l-threonine, and pyruvic acid, and three up-regulated compounds in CRC were identified, acetic acid, phenylacetylglutamine, and urea. The pathways and enrichment analysis returned only two pathways significantly expressed: the pyruvate metabolism and the glycolysis/gluconeogenesis pathway. In both cases, only the pyruvic acid (down-regulated in urine of CRC patients, with cancer cell proliferation effect in the tissue) and acetic acid (up-regulated in urine of CRC patients, with chemoprotective effect) were present.