Application of Differential Network Enrichment Analysis for Deciphering Metabolic Alterations

DNEA: an R package for fast and versatile data-driven network analysis of metabolomics data

BACKGROUND

Metabolomics is a high-throughput technology that measures small molecule metabolites in cells, tissues or biofluids. Analysis of metabolomics data is a multi-step process that involves data processing, quality control and normalization, followed by statistical and bioinformatics analysis. The latter step often involves pathway analysis to aid biological interpretation of the data. This approach is limited to endogenous metabolites that can be readily mapped to metabolic pathways. An alternative to pathway analysis that can be used for any classes of metabolites, including unknown compounds that are ubiquitous in untargeted metabolomics data, involves defining metabolite-metabolite interactions using experimental data. Our group has developed several network-based methods that use partial correlations of experimentally determined metabolite measurements. These were implemented in CorrelationCalculator and Filigree, two software tools for the analysis of metabolomics data we developed previously. The latter tool implements the Differential Network Enrichment Analysis (DNEA) algorithm. This analysis is useful for building differential networks from metabolomics data containing two experimental groups and identifying differentially enriched metabolic modules. While Filigree is a user-friendly tool, it has certain limitations when used for the analysis of large-scale metabolomics datasets.

RESULTS

We developed the DNEA R package for the data-driven network analysis of metabolomics data. We present the DNEA workflow and functionality, algorithm enhancements implemented with respect to the package's predecessor, Filigree, and discuss best practices for analyses. We tested the performance of the DNEA R package and illustrated its features using publicly available metabolomics data from the environmental determinants of diabetes in the young. To our knowledge, this package is the only publicly available tool designed for the construction of biological networks and subsequent enrichment testing for datasets containing exogenous, secondary, and unknown compounds. This greatly expands the scope of traditional enrichment analysis tools that can be used to analyze a relatively small set of well-annotated metabolites.

CONCLUSIONS

The DNEA R package is a more flexible and powerful implementation of our previously published software tool, Filigree. The modular structure of the package, along with the parallel processing framework built into the most computationally extensive steps of the algorithm, make it a powerful tool for the analysis of large and complex metabolomics datasets.

Metabolomic Alterations Associated with Phthalate Exposures among Pregnant Women in Puerto Rico

Although phthalate exposure has been linked with multiple adverse pregnancy outcomes, their underlying biological mechanisms are not fully understood. We examined associations between biomarkers of phthalate exposures and metabolic alterations using untargeted metabolomics in 99 pregnant women and 86 newborns [mean (SD) gestational age = 39.5 (1.5) weeks] in the PROTECT cohort. Maternal urinary phthalate metabolites were quantified using isotope dilution high-performance liquid chromatography-tandem mass spectrometry (LC-MS/MS), while metabolic profiles in maternal and cord blood plasma were characterized via reversed-phase LC-MS. Multivariable linear regression was used in metabolome-wide association studies (MWAS) to identify individual metabolic features associated with elevated phthalate levels, while clustering and correlation network analyses were used to discern the interconnectedness of biologically relevant features. In the MWAS adjusted for maternal age and prepregnancy BMI, we observed significant associations between specific phthalates, namely, di(2-ethylhexyl) phthalate (DEHP) and mono(3-carboxypropyl) phthalate (MCPP), and 34 maternal plasma metabolic features. These associations predominantly included upregulation of fatty acids, amino acids, purines, or their derivatives and downregulation of ceramides and sphingomyelins. In contrast, fewer significant associations were observed with metabolic features in cord blood. Correlation network analysis highlighted the overlap of features associated with phthalates and those identified as differentiating markers for preterm birth in a previous study. Overall, our findings underscore the complex impact of phthalate exposures on maternal and fetal metabolism, highlighting metabolomics as a tool for understanding associated biological processes. Future research should focus on expanding the sample size, exploring the effects of phthalate mixtures, and validating identified metabolic features in larger, more diverse populations.

CorrelationCalculator and Filigree: Tools for Data-Driven Network Analysis of Metabolomics Data

A significant challenge in the analysis of omics data is extracting actionable biological knowledge. Metabolomics is no exception. The general problem of relating changes in levels of individual metabolites to specific biological processes is compounded by the large number of unknown metabolites present in untargeted liquid chromatography-mass spectrometry (LC-MS) studies. Further, secondary metabolism and lipid metabolism are poorly represented in existing pathway databases. To overcome these limitations, our group has developed several tools for data-driven network construction and analysis. These include CorrelationCalculator and Filigree. Both tools allow users to build partial correlation-based networks from experimental metabolomics data when the number of metabolites exceeds the number of samples. CorrelationCalculator supports the construction of a single network, while Filigree allows building a differential network utilizing data from two groups of samples, followed by network clustering and enrichment analysis. We will describe the utility and application of both tools for the analysis of real-life metabolomics data.

Body mass index associates with amyotrophic lateral sclerosis survival and metabolomic profiles

INTRODUCTION/AIMS

Body mass index (BMI) is linked to amyotrophic lateral sclerosis (ALS) risk and prognosis, but additional research is needed. The aim of this study was to identify whether and when historical changes in BMI occurred in ALS participants, how these longer term trajectories associated with survival, and whether metabolomic profiles provided insight into potential mechanisms.

METHODS

ALS and control participants self-reported body height and weight 10 (reference) and 5 years earlier, and at study entry (diagnosis for ALS participants). Generalized estimating equations evaluated differences in BMI trajectories between cases and controls. ALS survival was evaluated by BMI trajectory group using accelerated failure time models. BMI trajectories and survival associations were explored using published metabolomic profiling and correlation networks.

RESULTS

Ten-year BMI trends differed between ALS and controls, with BMI loss in the 5 years before diagnosis despite BMI gains 10 to 5 years beforehand in both groups. An overall 10-year drop in BMI associated with a 27.1% decrease in ALS survival (P = .010). Metabolomic networks in ALS participants showed dysregulation in sphingomyelin, bile acid, and plasmalogen subpathways.

DISCUSSION

ALS participants lost weight in the 5-year period before enrollment. BMI trajectories had three distinct groups and the group with significant weight loss in the past 10 years had the worst survival. Participants with a high BMI and increase in weight in the 10 years before symptom onset also had shorter survival. Certain metabolomics profiles were associated with the BMI trajectories. Replicating these findings in prospective cohorts is warranted.

Metabolomics identifies shared lipid pathways in independent amyotrophic lateral sclerosis cohorts

Amyotrophic lateral sclerosis (ALS) is a fatal neurodegenerative disease lacking effective treatments. This is due, in part, to a complex and incompletely understood pathophysiology. To shed light, we conducted untargeted metabolomics on plasma from two independent cross-sectional ALS cohorts versus control participants to identify recurrent dysregulated metabolic pathways. Untargeted metabolomics was performed on plasma from two ALS cohorts (cohort 1, n = 125; cohort 2, n = 225) and healthy controls (cohort 1, n = 71; cohort 2, n = 104). Individual differential metabolites in ALS cases versus controls were assessed by Wilcoxon, adjusted logistic regression and partial least squares-discriminant analysis, while group lasso explored sub-pathway level differences. Adjustment parameters included age, sex and body mass index. Metabolomics pathway enrichment analysis was performed on metabolites selected using the above methods. Additionally, we conducted a sex sensitivity analysis due to sex imbalance in the cohort 2 control arm. Finally, a data-driven approach, differential network enrichment analysis (DNEA), was performed on a combined dataset to further identify important ALS metabolic pathways. Cohort 2 ALS participants were slightly older than the controls (64.0 versus 62.0 years, P = 0.009). Cohort 2 controls were over-represented in females (68%, P < 0.001). The most concordant cohort 1 and 2 pathways centred heavily on lipid sub-pathways, including complex and signalling lipid species and metabolic intermediates. There were differences in sub-pathways that were enriched in ALS females versus males, including in lipid sub-pathways. Finally, DNEA of the merged metabolite dataset of both ALS and control cohorts identified nine significant subnetworks; three centred on lipids and two encompassed a range of sub-pathways. In our analysis, we saw consistent and important shared metabolic sub-pathways in both ALS cohorts, particularly in lipids, further supporting their importance as ALS pathomechanisms and therapeutics targets.

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.

Comparing the Fasting and Random-Fed Metabolome Response to an Oral Glucose Tolerance Test in Children and Adolescents: Implications of Sex, Obesity, and Insulin Resistance

As the incidence of obesity and type 2 diabetes (T2D) is occurring at a younger age, studying adolescent nutrient metabolism can provide insights on the development of T2D. Metabolic challenges, including an oral glucose tolerance test (OGTT) can assess the effects of perturbations in nutrient metabolism. Here, we present alterations in the global metabolome in response to an OGTT, classifying the influence of obesity and insulin resistance (IR) in adolescents that arrived at the clinic fasted and in a random-fed state. Participants were recruited as lean (n = 55, aged 8–17 years, BMI percentile 5–85%) and overweight and obese (OVOB, n = 228, aged 8–17 years, BMI percentile ≥ 85%). Untargeted metabolomics profiled 246 annotated metabolites in plasma at t0 and t60 min during the OGTT. Our results suggest that obesity and IR influence the switch from fatty acid (FA) to glucose oxidation in response to the OGTT. Obesity was associated with a blunted decline of acylcarnitines and fatty acid oxidation intermediates. In females, metabolites from the Fasted and Random-Fed OGTT were associated with HOMA-IR, including diacylglycerols, leucine/isoleucine, acylcarnitines, and phosphocholines. Our results indicate that at an early age, obesity and IR may influence the metabolome dynamics in response to a glucose challenge.

Advantages of Studying the Metabolome in Response to Mixed-Macronutrient Challenges and Suggestions for Future Research Designs

Evaluating the postprandial response to a dietary challenge containing all macronutrients-carbohydrates, lipids, and protein-may provide stronger insights of metabolic health than a fasted measurement. Metabolomic profiling deepens the understanding of the homeostatic and adaptive response to a dietary challenge by classifying multiple metabolic pathways and biomarkers. A total of 26 articles were identified that measure the human blood metabolome or lipidome response to a mixed-macronutrient challenge. Most studies were cross-sectional, exploring the baseline and postprandial response to the dietary challenge. Large variations in study designs were reported, including the macronutrient and caloric composition of the challenge and the delivery of the challenge as a liquid shake or a solid meal. Most studies utilized a targeted metabolomics platform, assessing only a particular metabolic pathway, however, several studies utilized global metabolomics and lipidomics assays demonstrating the expansive postprandial response of the metabolome. The postprandial response of individual amino acids was largely dependent on the amino acid composition of the test meal, with the exception of alanine and proline, 2 nonessential amino acids. Long-chain fatty acids and unsaturated long-chain acylcarnitines rapidly decreased in response to the dietary challenges, representing the switch from fat to carbohydrate oxidation. Studies were reviewed that assessed the metabolome response in the context of obesity and metabolic diseases, providing insight on how weight status and disease influence the ability to cope with a nutrient load and return to homeostasis. Results demonstrate that the flexibility to respond to a substrate load is influenced by obesity and metabolic disease and flexibility alterations will be evident in downstream metabolites of fat, carbohydrate, and protein metabolism. In response, we propose suggestions for standardization between studies with the potential of creating a study exploring the postprandial response to a multitude of challenges with a variety of macronutrients.