A Single Visualization Technique for Displaying Multiple Metabolite-Phenotype Associations

Circulating N-formylmethionine and metabolic shift in critical illness: a multicohort metabolomics study

BACKGROUND

Cell stress promotes degradation of mitochondria which release danger-associated molecular patterns that are catabolized to N-formylmethionine. We hypothesized that in critically ill adults, the response to N-formylmethionine is associated with increases in metabolomic shift-related metabolites and increases in 28-day mortality.

METHODS

We performed metabolomics analyses on plasma from the 428-subject Correction of Vitamin D Deficiency in Critically Ill Patients trial (VITdAL-ICU) cohort and the 90-subject Brigham and Women's Hospital Registry of Critical Illness (RoCI) cohort. In the VITdAL-ICU cohort, we analyzed 983 metabolites at Intensive Care Unit (ICU) admission, day 3, and 7. In the RoCI cohort, we analyzed 411 metabolites at ICU admission. The association between N-formylmethionine and mortality was determined by adjusted logistic regression. The relationship between individual metabolites and N-formylmethionine abundance was assessed with false discovery rate correction via linear regression, linear mixed-effects, and Gaussian graphical models.

RESULTS

Patients with the top quartile of N-formylmethionine abundance at ICU admission had a significantly higher adjusted odds of 28-day mortality in the VITdAL-ICU (OR, 2.4; 95%CI 1.5-4.0; P = 0.001) and RoCI cohorts (OR, 5.1; 95%CI 1.4-18.7; P = 0.015). Adjusted linear regression shows that with increases in N-formylmethionine abundance at ICU admission, 55 metabolites have significant differences common to both the VITdAL-ICU and RoCI cohorts. With increased N-formylmethionine abundance, both cohorts had elevations in individual short-chain acylcarnitine, branched chain amino acid, kynurenine pathway, and pentose phosphate pathway metabolites.

CONCLUSIONS

The results indicate that circulating N-formylmethionine promotes a metabolic shift with heightened mortality that involves incomplete mitochondrial fatty acid oxidation, increased branched chain amino acid metabolism, and activation of the pentose phosphate pathway.

Uncovering PFAS and Other Xenobiotics in the Dark Metabolome Using Ion Mobility Spectrometry, Mass Defect Analysis, and Machine Learning

The identification of xenobiotics in nontargeted metabolomic analyses is a vital step in understanding human exposure. Xenobiotic metabolism, transformation, excretion, and coexistence with other endogenous molecules, however, greatly complicate the interpretation of features detected in nontargeted studies. While mass spectrometry (MS)-based platforms are commonly used in metabolomic measurements, deconvoluting endogenous metabolites from xenobiotics is also often challenged by the lack of xenobiotic parent and metabolite standards as well as the numerous isomers possible for each small molecule m/z feature. Here, we evaluate a xenobiotic structural annotation workflow using ion mobility spectrometry coupled with MS (IMS-MS), mass defect filtering, and machine learning to uncover potential xenobiotic classes and species in large metabolomic feature lists. Xenobiotic classes examined included those of known high toxicities, including per- and polyfluoroalkyl substances (PFAS), polycyclic aromatic hydrocarbons (PAHs), polychlorinated biphenyls (PCBs), polybrominated diphenyl ethers (PBDEs), and pesticides. Specifically, when the workflow was applied to identify PFAS in the NIST SRM 1957 and 909c human serum samples, it greatly reduced the hundreds of detected liquid chromatography (LC)-IMS-MS features by utilizing both mass defect filtering and m/z versus IMS collision cross sections relationships. These potential PFAS features were then compared to the EPA CompTox entries, and while some matched within specific m/z tolerances, there were still many unknowns illustrating the importance of nontargeted studies for detecting new molecules with known chemical characteristics. Additionally, this workflow can also be utilized to evaluate other xenobiotics and enable more confident annotations from nontargeted studies.

FlexDotPlot: a universal and modular dot plot visualization tool for complex multifaceted data

Motivation

Dot plots are heatmap-like charts that provide a compact way to simultaneously display two quantitative information by means of dots of different sizes and colors. Despite the popularity of this visualization method, particularly in single-cell RNA-sequencing (scRNA-seq) studies, existing tools used to make dot plots are limited in terms of functionality and usability.

Results

We developed FlexDotPlot, an R package for generating dot plots from multifaceted data, including scRNA-seq data. It provides a universal and easy-to-use solution with a high versatility. An interactive R Shiny application is also available allowing non-R users to easily generate dot plots with several tunable parameters.

Availability and implementation

Source code and detailed manual are available on CRAN (stable version) and at https://github.com/Simon-Leonard/FlexDotPlot (development version). Code to reproduce figures is available at https://github.com/Simon-Leonard/FlexDotPlot_paper. A Shiny app is available as a stand-alone application within the package.

Supplementary information

Supplementary data are available at Bioinformatics Advances online.

Sex-Specific Catabolic Metabolism Alterations in the Critically Ill following High Dose Vitamin D

Pharmacological interventions are essential for the treatment and management of critical illness. Although women comprise a large proportion of the critically ill, sex-specific pharmacological properties are poorly described in critical care. The sex-specific effects of vitamin D₃ treatment in the critically ill are not known. Therefore, we performed a metabolomics cohort study with 1215 plasma samples from 428 patients from the VITdAL-ICU trial to study sex-specific differences in the metabolic response to critical illness following high-dose oral vitamin D₃ intervention. In women, despite the dose of vitamin D₃ being higher, pharmacokinetics demonstrated a lower extent of vitamin D₃ absorption compared to men. Metabolic response to high-dose oral vitamin D₃ is sex-specific. Sex-stratified individual metabolite associations with elevations in 25(OH)D following intervention showed female-specific positive associations in long-chain acylcarnitines and male-specific positive associations in free fatty acids. In subjects who responded to vitamin D₃ intervention, significant negative associations were observed in short-chain acylcarnitines and branched chain amino acid metabolites in women as compared to men. Acylcarnitines and branched chain amino acids are reflective of fatty acid B oxidation, and bioenergesis may represent notable metabolic signatures of the sex-specific response to vitamin D. Demonstrating sex-specific pharmacometabolomics differences following intervention is an important movement towards the understanding of personalized medicine.

Metabolic Impairment in Coronary Artery Disease: Elevated Serum Acylcarnitines Under the Spotlights

Coronary artery disease (CAD) remains the leading cause of death worldwide. Expanding patients' metabolic phenotyping beyond clinical chemistry investigations could lead to earlier recognition of disease onset and better prevention strategies. Additionally, metabolic phenotyping, at the molecular species level, contributes to unravel the roles of metabolites in disease development. In this cross-sectional study, we investigated clinically healthy individuals (n = 116, 65% male, 70.8 ± 8.7 years) and patients with CAD (n = 54, 91% male, 67.0 ± 11.5 years) of the COmPLETE study. We applied a high-coverage quantitative liquid chromatography-mass spectrometry approach to acquire a comprehensive profile of serum acylcarnitines, free carnitine and branched-chain amino acids (BCAAs), as markers of mitochondrial health and energy homeostasis. Multivariable linear regression analyses, adjusted for confounders, were conducted to assess associations between metabolites and CAD phenotype. In total, 20 short-, medium- and long-chain acylcarnitine species, along with L-carnitine, valine and isoleucine were found to be significantly (adjusted p ≤ 0.05) and positively associated with CAD. For 17 acylcarnitine species, associations became stronger as the number of affected coronary arteries increased. This implies that circulating acylcarnitine levels reflect CAD severity and might play a role in future patients' stratification strategies. Altogether, CAD is characterized by elevated serum acylcarnitine and BCAA levels, which indicates mitochondrial imbalance between fatty acid and glucose oxidation.

Procalcitonin metabolomics in the critically ill reveal relationships between inflammation intensity and energy utilization pathways

Procalcitonin is a biomarker of systemic inflammation and may have importance in the immune response. The metabolic response to elevated procalcitonin in critical illness is not known. The response to inflammation is vitally important to understanding metabolism alterations during extreme stress. Our aim was to determine if patients with elevated procalcitonin have differences in the metabolomic response to early critical illness. We performed a metabolomics study of the VITdAL-ICU trial where subjects received high dose vitamin D3 or placebo. Mixed-effects modeling was used to study changes in metabolites over time relative to procalcitonin levels adjusted for age, Simplified Acute Physiology Score II, admission diagnosis, day 0 25-hydroxyvitamin D level, and the 25-hydroxyvitamin D response to intervention. With elevated procalcitonin, multiple members of the short and medium chain acylcarnitine, dicarboxylate fatty acid, branched-chain amino acid, and pentose phosphate pathway metabolite classes had significantly positive false discovery rate corrected associations. Further, multiple long chain acylcarnitines and lysophosphatidylcholines had significantly negative false discovery rate corrected associations with elevated procalcitonin. Gaussian graphical model analysis revealed functional modules specific to elevated procalcitonin. Our findings show that metabolite differences exist with increased procalcitonin indicating activation of branched chain amino acid dehydrogenase and a metabolic shift.

Metabolic View on Human Healthspan: A Lipidome-Wide Association Study

As ageing is a major risk factor for the development of non-communicable diseases, extending healthspan has become a medical and societal necessity. Precise lipid phenotyping that captures metabolic individuality could support healthspan extension strategies. This study applied ‘omic-scale lipid profiling to characterise sex-specific age-related differences in the serum lipidome composition of healthy humans. A subset of the COmPLETE-Health study, composed of 73 young (25.2 ± 2.6 years, 43% female) and 77 aged (73.5 ± 2.3 years, 48% female) clinically healthy individuals, was investigated, using an untargeted liquid chromatography high-resolution mass spectrometry approach. Compared to their younger counterparts, aged females and males exhibited significant higher levels in 138 and 107 lipid species representing 15 and 13 distinct subclasses, respectively. Percentage of difference ranged from 5.8% to 61.7% (females) and from 5.3% to 46.0% (males), with sphingolipid and glycerophophospholipid species displaying the greatest amplitudes. Remarkably, specific sphingolipid and glycerophospholipid species, previously described as cardiometabolically favourable, were found elevated in aged individuals. Furthermore, specific ether-glycerophospholipid and lyso-glycerophosphocholine species displayed higher levels in aged females only, revealing a more favourable lipidome evolution in females. Altogether, age determined the circulating lipidome composition, while lipid species analysis revealed additional findings that were not observed at the subclass level.

Metabolomic differences between critically Ill women and men

Metabolism differs in women and men at homeostasis. Critically ill patients have profound dysregulation of homeostasis and metabolism. It is not clear if the metabolic response to critical illness differs in women compared to men. Such sex-specific differences in illness response would have consequences for personalized medicine. Our aim was to determine the sex-specific metabolomic response to early critical illness. We performed a post-hoc metabolomics study of the VITdAL-ICU trial where subjects received high dose vitamin D3 or placebo. Using mixed-effects modeling, we studied sex-specific changes in metabolites over time adjusted for age, Simplified Acute Physiology Score II, admission diagnosis, day 0 25-hydroxyvitamin D level, and 25-hydroxyvitamin D response to intervention. In women, multiple members of the sphingomyelin and lysophospholipid metabolite classes had significantly positive Bonferroni corrected associations over time compared to men. Further, multiple representatives of the acylcarnitine, androgenic steroid, bile acid, nucleotide and amino acid metabolite classes had significantly negative Bonferroni corrected associations over time compared to men. Gaussian graphical model analyses revealed sex-specific functional modules. Our findings show that robust and coordinated sex-specific metabolite differences exist early in critical illness.

One-Year Effects of Omega-3 Treatment on Fatty Acids, Oxylipins, and Related Bioactive Lipids and Their Associations with Clinical Lipid and Inflammatory Biomarkers: Findings from a Substudy of the Vitamin D and Omega-3 Trial (VITAL)

Omega-3 (n-3) treatment may lower cardiovascular risk, yet its effects on the circulating lipidome and relation to cardiovascular risk biomarkers are unclear. We hypothesized that n-3 treatment is associated with favorable changes in downstream fatty acids (FAs), oxylipins, bioactive lipids, clinical lipid and inflammatory biomarkers. We examined these VITAL200, a nested substudy of 200 subjects balanced on demographics and treatment and randomly selected from the Vitamin D and Omega-3 Trial (VITAL). VITAL is a randomized double-blind trial of 840 mg/d eicosapentaenoic acid (EPA) + docosahexaenoic acid (DHA) vs. placebo among 25,871 individuals. Small polar bioactive lipid features, oxylipins and FAs from plasma and red blood cells were measured using three independent assaying techniques at baseline and one year. The Women's Health Study (WHS) was used for replication with dietary n-3 intake. Randomized n-3 treatment led to changes in 143 FAs, oxylipins and bioactive lipids (False Discovery Rate (FDR) < 0.05 in VITAL200, validated (p-values < 0.05)) in WHS with increases in 95 including EPA, DHA, n-3 docosapentaenoic acid (DPA-n3), and decreases in 48 including DPA-n6, dihomo gamma linolenic (DGLA), adrenic and arachidonic acids. N-3 related changes in the bioactive lipidome were heterogeneously associated with changes in clinical lipid and inflammatory biomarkers. N-3 treatment significantly modulates the bioactive lipidome, which may contribute to its clinical benefits.

Metabolomic basis for response to high dose vitamin D in critical illness

BACKGROUND & AIMS

It is unclear if intervention can mitigate the dramatic alterations of metabolic homeostasis present in critical illness. Our objective was to determine the associations between increased 25-hydroxyvitamin D levels following high dose vitamin D₃ and more favorable metabolomic profiles in critical illness.

METHODS

We performed a post-hoc metabolomics study of the VITdAL-ICU randomized double-blind, placebo-controlled trial. Trial patients from Medical and Surgical Intensive Care Units at a tertiary university hospital with 25-hydroxyvitamin D level ≤20 ng/mL received either high dose oral vitamin D₃ (540,000 IU) or placebo. We performed an analysis of 578 metabolites from 1215 plasma samples from 428 subjects at randomization (day 0), day 3 and 7. Using mixed-effects modeling, we studied changes in metabolite profiles in subjects receiving intervention or placebo relative to absolute increases in 25-hydroxyvitamin D levels from day 0 to day 3.

RESULTS

55.2% of subjects randomized to high dose vitamin D₃ demonstrated an absolute increase in 25-hydroxyvitamin D ≥ 15 ng/ml from day 0 to day 3. With an absolute increase in 25-hydroxyvitamin D ≥ 15 ng/ml, multiple members of the sphingomyelin, plasmalogen, lysoplasmalogen and lysophospholipid metabolite classes had significantly positive Bonferroni corrected associations over time. Further, multiple representatives of the acylcarnitine and phosphatidylethanolamine metabolite classes had significantly negative Bonferroni corrected associations over time with an absolute increase in 25-hydroxyvitamin D ≥ 15 ng/ml. Changes in these highlighted metabolite classes were associated with decreased 28-day mortality.

CONCLUSIONS

Increases in 25-hydroxyvitamin D following vitamin D₃ intervention are associated with favorable changes in metabolites involved in endothelial protection, enhanced innate immunity and improved mitochondrial function.

Blood Metabolite Signature of Metabolic Syndrome Implicates Alterations in Amino Acid Metabolism: Findings from the Baltimore Longitudinal Study of Aging (BLSA) and the Tsuruoka Metabolomics Cohort Study (TMCS)

Rapid lifestyle and dietary changes have contributed to a rise in the global prevalence of metabolic syndrome (MetS), which presents a potential healthcare crisis, owing to its association with an increased burden of multiple cardiovascular and neurological diseases. Prior work has identified the role that genetic, lifestyle, and environmental factors can play in the prevalence of MetS. Metabolomics is an important tool to study alterations in biochemical pathways intrinsic to the pathophysiology of MetS. We undertook a metabolomic study of MetS in serum samples from two ethnically distinct, well-characterized cohorts-the Baltimore Longitudinal Study of Aging (BLSA) from the U.S. and the Tsuruoka Metabolomics Cohort Study (TMCS) from Japan. We used multivariate logistic regression to identify metabolites that were associated with MetS in both cohorts. Among the top 25 most significant (lowest p-value) metabolite associations with MetS in each cohort, we identified 18 metabolites that were shared between TMCS and BLSA, the majority of which were classified as amino acids. These associations implicate multiple biochemical pathways in MetS, including branched-chain amino acid metabolism, glutathione production, aromatic amino acid metabolism, gluconeogenesis, and the tricarboxylic acid cycle. Our results suggest that fundamental alterations in amino acid metabolism may be central features of MetS.