maplet: an extensible R toolbox for modular and reproducible metabolomics pipelines

AutoFocus: a hierarchical framework to explore multi-omic disease associations spanning multiple scales of biomolecular interaction

Recent advances in high-throughput measurement technologies have enabled the analysis of molecular perturbations associated with disease phenotypes at the multi-omic level. Such perturbations can range in scale from fluctuations of individual molecules to entire biological pathways. Data-driven clustering algorithms have long been used to group interactions into interpretable functional modules; however, these modules are typically constrained to a fixed size or statistical cutoff. Furthermore, modules are often analyzed independently of their broader biological context. Consequently, such clustering approaches limit the ability to explore functional module associations with disease phenotypes across multiple scales. Here, we introduce AutoFocus, a data-driven method that hierarchically organizes biomolecules and tests for phenotype enrichment at every level within the hierarchy. As a result, the method allows disease-associated modules to emerge at any scale. We evaluated this approach using two datasets: First, we explored associations of biomolecules from the multi-omic QMDiab dataset (n = 388) with the well-characterized type 2 diabetes phenotype. Secondly, we utilized the ROS/MAP Alzheimer's disease dataset (n = 500), consisting of high-throughput measurements of brain tissue to explore modules associated with multiple Alzheimer's Disease-related phenotypes. Our method identifies modules that are multi-omic, span multiple pathways, and vary in size. We provide an interactive tool to explore this hierarchy at different levels and probe enriched modules, empowering users to examine the full hierarchy, delve into biomolecular drivers of disease phenotype within a module, and incorporate functional annotations.

The HuMet Repository: Watching human metabolism at work

Metabolism oscillates between catabolic and anabolic states depending on food intake, exercise, or stresses that change a multitude of metabolic pathways simultaneously. We present the HuMet Repository for exploring dynamic metabolic responses to oral glucose/lipid loads, mixed meals, 36-h fasting, exercise, and cold stress in healthy subjects. Metabolomics data from blood, urine, and breath of 15 young, healthy men at up to 56 time points are integrated and embedded within an interactive web application, enabling researchers with and without computational expertise to search, visualize, analyze, and contextualize the dynamic metabolite profiles of 2,656 metabolites acquired on multiple platforms. With examples, we demonstrate the utility of the resource for research into the dynamics of human metabolism, highlighting differences and similarities in systemic metabolic responses across challenges and the complementarity of metabolomics platforms. The repository, providing a reference for healthy metabolite changes to six standardized physiological challenges, is freely accessible through a web portal.

Serum and CSF metabolomics analysis shows Mediterranean Ketogenic Diet mitigates risk factors of Alzheimer's disease

Alzheimer's disease (AD) is influenced by a variety of modifiable risk factors, including a person's dietary habits. While the ketogenic diet (KD) holds promise in reducing metabolic risks and potentially affecting AD progression, only a few studies have explored KD's metabolic impact, especially on blood and cerebrospinal fluid (CSF). Our study involved participants at risk for AD, either cognitively normal or with mild cognitive impairment. The participants consumed both a modified Mediterranean Ketogenic Diet (MMKD) and the American Heart Association diet (AHAD) for 6 weeks each, separated by a 6-week washout period. We employed nuclear magnetic resonance (NMR)-based metabolomics to profile serum and CSF and metagenomics profiling on fecal samples. While the AHAD induced no notable metabolic changes, MMKD led to significant alterations in both serum and CSF. These changes included improved modifiable risk factors, like increased HDL-C and reduced BMI, reversed serum metabolic disturbances linked to AD such as a microbiome-mediated increase in valine levels, and a reduction in systemic inflammation. Additionally, the MMKD was linked to increased amino acid levels in the CSF, a breakdown of branched-chain amino acids (BCAAs), and decreased valine levels. Importantly, we observed a strong correlation between metabolic changes in the CSF and serum, suggesting a systemic regulation of metabolism. Our findings highlight that MMKD can improve AD-related risk factors, reverse some metabolic disturbances associated with AD, and align metabolic changes across the blood-CSF barrier.

A Modified Mediterranean Ketogenic Diet mitigates modifiable risk factors of Alzheimer's Disease: a serum and CSF-based metabolic analysis

Alzheimer's disease (AD) is influenced by a variety of modifiable risk factors, including a person's dietary habits. While the ketogenic diet (KD) holds promise in reducing metabolic risks and potentially affecting AD progression, only a few studies have explored KD's metabolic impact, especially on blood and cerebrospinal fluid (CSF). Our study involved participants at risk for AD, either cognitively normal or with mild cognitive impairment. The participants consumed both a modified Mediterranean-ketogenic diet (MMKD) and the American Heart Association diet (AHAD) for 6 weeks each, separated by a 6-week washout period. We employed nuclear magnetic resonance (NMR)-based metabolomics to profile serum and CSF and metagenomics profiling on fecal samples. While the AHAD induced no notable metabolic changes, MMKD led to significant alterations in both serum and CSF. These changes included improved modifiable risk factors, like increased HDL-C and reduced BMI, reversed serum metabolic disturbances linked to AD such as a microbiome-mediated increase in valine levels, and a reduction in systemic inflammation. Additionally, the MMKD was linked to increased amino acid levels in the CSF, a breakdown of branched-chain amino acids (BCAAs), and decreased valine levels. Importantly, we observed a strong correlation between metabolic changes in the CSF and serum, suggesting a systemic regulation of metabolism. Our findings highlight that MMKD can improve AD-related risk factors, reverse some metabolic disturbances associated with AD, and align metabolic changes across the blood-CSF barrier.

Plasma metabolomics profiling of 580 patients from an Early Detection Research Network prostate cancer cohort

Prostate cancer is the second most common cancer in men and affects 1 in 9 men in the United States. Early screening for prostate cancer often involves monitoring levels of prostate-specific antigen (PSA) and performing digital rectal exams. However, a prostate biopsy is always required for definitive cancer diagnosis. The Early Detection Research Network (EDRN) is a consortium within the National Cancer Institute aimed at improving screening approaches and early detection of cancers. As part of this effort, the Weill Cornell EDRN Prostate Cancer has collected and biobanked specimens from men undergoing a prostate biopsy between 2008 and 2017. In this report, we describe blood metabolomics measurements for a subset of this population. The dataset includes detailed clinical and prospective records for 580 patients who underwent prostate biopsy, 287 of which were subsequentially diagnosed with prostate cancer, combined with profiling of 1,482 metabolites from plasma samples collected at the time of biopsy. We expect this dataset to provide a valuable resource for scientists investigating prostate cancer metabolism.

Comparative brain metabolomics reveals shared and distinct metabolic alterations in Alzheimer's disease and progressive supranuclear palsy

Metabolic dysregulation is a hallmark of neurodegenerative diseases, including Alzheimer's disease (AD) and progressive supranuclear palsy (PSP). While metabolic dysregulation is a common link between these two tauopathies, a comprehensive brain metabolic comparison of the diseases has not yet been performed. We analyzed 342 postmortem brain samples from the Mayo Clinic Brain Bank and examined 658 metabolites in the cerebellar cortex and the temporal cortex between the two tauopathies. Our findings indicate that both diseases display oxidative stress associated with lipid metabolism, mitochondrial dysfunction linked to lysine metabolism, and an indication of tau-induced polyamine stress response. However, specific to AD, we detected glutathione-related neuroinflammation, deregulations of enzymes tied to purines, and cognitive deficits associated with vitamin B. Taken together, our findings underscore vast alterations in the brain's metabolome, illuminating shared neurodegenerative pathways and disease-specific traits in AD and PSP.

Urine-based multi-omic comparative analysis of COVID-19 and bacterial sepsis-induced ARDS

BACKGROUND

Acute respiratory distress syndrome (ARDS), a life-threatening condition during critical illness, is a common complication of COVID-19. It can originate from various disease etiologies, including severe infections, major injury, or inhalation of irritants. ARDS poses substantial clinical challenges due to a lack of etiology-specific therapies, multisystem involvement, and heterogeneous, poor patient outcomes. A molecular comparison of ARDS groups holds the potential to reveal common and distinct mechanisms underlying ARDS pathogenesis.

METHODS

We performed a comparative analysis of urine-based metabolomics and proteomics profiles from COVID-19 ARDS patients (n = 42) and bacterial sepsis-induced ARDS patients (n = 17). To this end, we used two different approaches, first we compared the molecular omics profiles between ARDS groups, and second, we correlated clinical manifestations within each group with the omics profiles.

RESULTS

The comparison of the two ARDS etiologies identified 150 metabolites and 70 proteins that were differentially abundant between the two groups. Based on these findings, we interrogated the interplay of cell adhesion/extracellular matrix molecules, inflammation, and mitochondrial dysfunction in ARDS pathogenesis through a multi-omic network approach. Moreover, we identified a proteomic signature associated with mortality in COVID-19 ARDS patients, which contained several proteins that had previously been implicated in clinical manifestations frequently linked with ARDS pathogenesis.

CONCLUSION

In summary, our results provide evidence for significant molecular differences in ARDS patients from different etiologies and a potential synergy of extracellular matrix molecules, inflammation, and mitochondrial dysfunction in ARDS pathogenesis. The proteomic mortality signature should be further investigated in future studies to develop prediction models for COVID-19 patient outcomes.

rox: A Statistical Model for Regression with Missing Values

High-dimensional omics datasets frequently contain missing data points, which typically occur due to concentrations below the limit of detection (LOD) of the profiling platform. The presence of such missing values significantly limits downstream statistical analysis and result interpretation. Two common techniques to deal with this issue include the removal of samples with missing values and imputation approaches that substitute the missing measurements with reasonable estimates. Both approaches, however, suffer from various shortcomings and pitfalls. In this paper, we present "rox", a novel statistical model for the analysis of omics data with missing values without the need for imputation. The model directly incorporates missing values as "low" concentrations into the calculation. We show the superiority of rox over common approaches on simulated data and on six metabolomics datasets. Fully leveraging the information contained in LOD-based missing values, rox provides a powerful tool for the statistical analysis of omics data.

Microbes Contribute to Chemopreventive Efficacy, Intestinal Tumorigenesis, and the Metabolome

Bacteria are believed to play an important role in intestinal tumorigenesis and contribute to both gut luminal and circulating metabolites. Celecoxib, a selective cyclooxygenase-2 inhibitor, alters gut bacteria and metabolites in association with suppressing the development of intestinal polyps in mice. The current study sought to evaluate whether celecoxib exerts its chemopreventive effects, in part, through intestinal bacteria and metabolomic alterations. Using ApcMin/+ mice, we demonstrated that treatment with broad-spectrum antibiotics (ABx) reduced abundance of gut bacteria and attenuated the ability of celecoxib to suppress intestinal tumorigenesis. Use of ABx also impaired celecoxib's ability to shift microbial populations and gut luminal and circulating metabolites. Treatment with ABx alone markedly reduced tumor number and size in ApcMin/+ mice, in conjunction with profoundly altering the metabolite profiles of the intestinal lumen and blood. Many of the metabolite changes in the gut and circulation overlapped and included shifts in microbially derived metabolites. To complement these findings in mice, we evaluated the effects of ABx on circulating metabolites in patients with colon cancer. This showed that ABx treatment led to a shift in blood metabolites, including several that were of bacterial origin. Importantly, changes in metabolites in patients given ABx overlapped with alterations found in mice that also received ABx. Taken together, these findings suggest a potential role for bacterial metabolites in mediating both the chemopreventive effects of celecoxib and intestinal tumor growth.

PREVENTION RELEVANCE

This study demonstrates novel mechanisms by which chemopreventive agents exert their effects and gut microbiota impact intestinal tumor development. These findings have the potential to lead to improved cancer prevention strategies by modulating microbes and their metabolites.

Hopomics: Humulus lupulus Brewing Cultivars Classification Based on LC-MS Profiling and Nested Feature Selection

Omics approaches in plant analysis find many different applications, from classification to new bioactive compounds discovery. Metabolomics seems to be one of the most informative ways of describing plants’ phenotypes, since commonly used methods such as liquid chromatography–mass spectrometry (LC-MS) and nuclear magnetic resonance spectroscopy (NMR) could provide a huge amount of information about samples. However, due to high efficiency, many disadvantages arise with the complexity of the experimental design. In the present work, we demonstrate an untargeted metabolomics pipeline with the example of a Humulus lupulus classification task. LC-MS profiling of brewing cultivars samples was carried out as a starting point. Hierarchical cluster analysis (HCA)-based classification in combination with nested feature selection was provided for sample discrimination and marker compounds discovery. Obtained metabolome-based classification showed an expected difference compared to genetic-based classification data. Nine compounds were found to have the biggest classification power during nested feature selection. Using database search and molecular network construction, five of them were identified as known hops bitter compounds.

metabolomicsR: a streamlined workflow to analyze metabolomic data in R

Summary

metabolomicsR is a streamlined, flexible and user-friendly R package to preprocess, analyze and visualize metabolomic data. metabolomicsR includes comprehensive functionalities for sample and metabolite quality control, outlier detection, missing value imputation, dimensional reduction, batch effect normalization, data integration, regression, metabolite annotation and visualization of data and results. In this application note, we demonstrate the step-by-step use of the main functions from this package.

Availability and implementation

The metabolomicsR package is available via CRAN and GitHub (https://github.com/XikunHan/metabolomicsR/). A step-by-step online tutorial is available at https://xikunhan.github.io/metabolomicsR/docs/articles/Introduction.html.

Supplementary information

Supplementary data are available at Bioinformatics Advances online.

Multi-omic comparative analysis of COVID-19 and bacterial sepsis-induced ARDS

BACKGROUND

Acute respiratory distress syndrome (ARDS), a life-threatening condition characterized by hypoxemia and poor lung compliance, is associated with high mortality. ARDS induced by COVID-19 has similar clinical presentations and pathological manifestations as non-COVID-19 ARDS. However, COVID-19 ARDS is associated with a more protracted inflammatory respiratory failure compared to traditional ARDS. Therefore, a comprehensive molecular comparison of ARDS of different etiologies groups may pave the way for more specific clinical interventions.

METHODS AND FINDINGS

In this study, we compared COVID-19 ARDS (n = 43) and bacterial sepsis-induced (non-COVID-19) ARDS (n = 24) using multi-omic plasma profiles covering 663 metabolites, 1,051 lipids, and 266 proteins. To address both between- and within- ARDS group variabilities we followed two approaches. First, we identified 706 molecules differently abundant between the two ARDS etiologies, revealing more than 40 biological processes differently regulated between the two groups. From these processes, we assembled a cascade of therapeutically relevant pathways downstream of sphingosine metabolism. The analysis suggests a possible overactivation of arginine metabolism involved in long-term sequelae of ARDS and highlights the potential of JAK inhibitors to improve outcomes in bacterial sepsis-induced ARDS. The second part of our study involved the comparison of the two ARDS groups with respect to clinical manifestations. Using a data-driven multi-omic network, we identified signatures of acute kidney injury (AKI) and thrombocytosis within each ARDS group. The AKI-associated network implicated mitochondrial dysregulation which might lead to post-ARDS renal-sequalae. The thrombocytosis-associated network hinted at a synergy between prothrombotic processes, namely IL-17, MAPK, TNF signaling pathways, and cell adhesion molecules. Thus, we speculate that combination therapy targeting two or more of these processes may ameliorate thrombocytosis-mediated hypercoagulation.

CONCLUSION

We present a first comprehensive molecular characterization of differences between two ARDS etiologies-COVID-19 and bacterial sepsis. Further investigation into the identified pathways will lead to a better understanding of the pathophysiological processes, potentially enabling novel therapeutic interventions.

Multi-omic comparative analysis of COVID-19 and bacterial sepsis-induced ARDS

Background

Acute respiratory distress syndrome (ARDS), a life-threatening condition characterized by hypoxemia and poor lung compliance, is associated with high mortality. ARDS induced by COVID-19 has similar clinical presentations and pathological manifestations as non-COVID-19 ARDS. However, COVID-19 ARDS is associated with a more protracted inflammatory respiratory failure compared to traditional ARDS. Therefore, a comprehensive molecular comparison of ARDS of different etiologies groups may pave the way for more specific clinical interventions.

Methods and Findings

In this study, we compared COVID-19 ARDS (n=43) and bacterial sepsis-induced (non-COVID-19) ARDS (n=24) using multi-omic plasma profiles covering 663 metabolites, 1,051 lipids, and 266 proteins. To address both between- and within-ARDS group variabilities we followed two approaches. First, we identified 706 molecules differently abundant between the two ARDS etiologies, revealing more than 40 biological processes differently regulated between the two groups. From these processes, we assembled a cascade of therapeutically relevant pathways downstream of sphingosine metabolism. The analysis suggests a possible overactivation of arginine metabolism involved in long-term sequelae of ARDS and highlights the potential of JAK inhibitors to improve outcomes in bacterial sepsis-induced ARDS. The second part of our study involved the comparison of the two ARDS groups with respect to clinical manifestations. Using a data-driven multi-omic network, we identified signatures of acute kidney injury (AKI) and thrombocytosis within each ARDS group. The AKI-associated network implicated mitochondrial dysregulation which might lead to post-ARDS renal-sequalae. The thrombocytosis-associated network hinted at a synergy between prothrombotic processes, namely IL-17, MAPK, TNF signaling pathways, and cell adhesion molecules. Thus, we speculate that combination therapy targeting two or more of these processes may ameliorate thrombocytosis-mediated hypercoagulation.

Conclusion

We present a first comprehensive molecular characterization of differences between two ARDS etiologies - COVID-19 and bacterial sepsis. Further investigation into the identified pathways will lead to a better understanding of the pathophysiological processes, potentially enabling novel therapeutic interventions.

Urine-based multi-omic comparative analysis of COVID-19 and bacterial sepsis-induced ARDS

Acute respiratory distress syndrome (ARDS), a life-threatening condition during critical illness, is a common complication of COVID-19. It can originate from various disease etiologies, including severe infections, major injury, or inhalation of irritants. ARDS poses substantial clinical challenges due to a lack of etiology-specific therapies, multisystem involvement, and heterogeneous, poor patient outcomes. A molecular comparison of ARDS groups holds the potential to reveal common and distinct mechanisms underlying ARDS pathogenesis. In this study, we performed a comparative analysis of urine-based metabolomics and proteomics profiles from COVID-19 ARDS patients (n = 42) and bacterial sepsis-induced ARDS patients (n = 17). The comparison of these ARDS etiologies identified 150 metabolites and 70 proteins that were differentially abundant between the two groups. Based on these findings, we interrogated the interplay of cell adhesion/extracellular matrix molecules, inflammation, and mitochondrial dysfunction in ARDS pathogenesis through a multi-omic network approach. Moreover, we identified a proteomic signature associated with mortality in COVID-19 ARDS patients, which contained several proteins that had previously been implicated in clinical manifestations frequently linked with ARDS pathogenesis. In summary, our results provide evidence for significant molecular differences in ARDS patients from different etiologies and a potential synergy of extracellular matrix molecules, inflammation, and mitochondrial dysfunction in ARDS pathogenesis. The proteomic mortality signature should be further investigated in future studies to develop prediction models for COVID-19 patient outcomes.

Angiopoietin 2 Is Associated with Vascular Necroptosis Induction in Coronavirus Disease 2019 Acute Respiratory Distress Syndrome

Vascular injury is a well-established, disease-modifying factor in acute respiratory distress syndrome (ARDS) pathogenesis. Recently, coronavirus disease 2019 (COVID-19)-induced injury to the vascular compartment has been linked to complement activation, microvascular thrombosis, and dysregulated immune responses. This study sought to assess whether aberrant vascular activation in this prothrombotic context was associated with the induction of necroptotic vascular cell death. To achieve this, proteomic analysis was performed on blood samples from COVID-19 subjects at distinct time points during ARDS pathogenesis (hospitalized at risk, N = 59; ARDS, N = 31; and recovery, N = 12). Assessment of circulating vascular markers in the at-risk cohort revealed a signature of low vascular protein abundance that tracked with low platelet levels and increased mortality. This signature was replicated in the ARDS cohort and correlated with increased plasma angiopoietin 2 levels. COVID-19 ARDS lung autopsy immunostaining confirmed a link between vascular injury (angiopoietin 2) and platelet-rich microthrombi (CD61) and induction of necrotic cell death [phosphorylated mixed lineage kinase domain-like (pMLKL)]. Among recovery subjects, the vascular signature identified patients with poor functional outcomes. Taken together, this vascular injury signature was associated with low platelet levels and increased mortality and can be used to identify ARDS patients most likely to benefit from vascular targeted therapies.

metaboprep: an R package for preanalysis data description and processing

MOTIVATION

Metabolomics is an increasingly common part of health research and there is need for preanalytical data processing. Researchers typically need to characterize the data and to exclude errors within the context of the intended analysis. Whilst some preprocessing steps are common, there is currently a lack of standardization and reporting transparency for these procedures.

RESULTS

Here, we introduce metaboprep, a standardized data processing workflow to extract and characterize high quality metabolomics datasets. The package extracts data from preformed worksheets, provides summary statistics and enables the user to select samples and metabolites for their analysis based on a set of quality metrics. A report summarizing quality metrics and the influence of available batch variables on the data are generated for the purpose of open disclosure. Where possible, we provide users flexibility in defining their own selection thresholds.

AVAILABILITY AND IMPLEMENTATION

metaboprep is an open-source R package available at https://github.com/MRCIEU/metaboprep.

SUPPLEMENTARY INFORMATION

Supplementary data are available at Bioinformatics online.