The nPYc-Toolbox, a Python module for the pre-processing, quality-control and analysis of metabolic profiling datasets

Adipose tissue-derived metabolite risk scores and risk for type 2 diabetes in South Asians

BACKGROUND

South Asians are at higher risk for type 2 diabetes (T2D) than many other race/ethnic groups. Ectopic adiposity, specifically hepatic steatosis and visceral fat may partially explain this. Our objective was to derive metabolite risk scores for ectopic adiposity and assess associations with incident T2D in South Asians.

METHODS

We examined 550 participants in the Mediators of Atherosclerosis in South Asians Living in America (MASALA) cohort study aged 40-84 years without known cardiovascular disease or T2D and with metabolomic data. Computed tomography scans at baseline assessed hepatic attenuation and visceral fat area, and fasting serum specimens at baseline and after 5 years assessed T2D. LC-MS-based untargeted metabolomic analysis was performed followed by targeted integration and reporting of known signals. Elastic net regularized linear regression analyses was used to derive risk scores for hepatic steatosis and visceral fat using weighted coefficients. Logistic regression models associated metabolite risk score and incident T2D, adjusting for age, gender, study site, BMI, physical activity, diet quality, energy intake and use of cholesterol-lowering medication.

RESULTS

Average age of participants was 55 years, 36% women with an average body mass index (BMI) of 25 kg/m2 and 6% prevalence of hepatic steatosis, with 47 cases of incident T2D at 5 years. There were 445 metabolites of known identity. Of these, 313 metabolites were included in the MET-Visc score and 267 in the MET-Liver score. In most fully adjusted models, MET-Liver (OR 2.04 [95% CI 1.38, 3.03]) and MET-Visc (OR 2.80 [1.75, 4.46]) were associated with higher odds of T2D. These associations remained significant after adjustment for measured adiposity.

CONCLUSIONS

Metabolite risk scores for intrahepatic fat and visceral fat were strongly related to incident T2D independent of measured adiposity. Use of these biomarkers to target risk stratification may help capture pre-clinical metabolic abnormalities.

Demonstrating the reliability of in vivo metabolomics based chemical grouping: towards best practice

While grouping/read-across is widely used to fill data gaps, chemical registration dossiers are often rejected due to weak category justifications based on structural similarity only. Metabolomics provides a route to robust chemical categories via evidence of shared molecular effects across source and target substances. To gain international acceptance, this approach must demonstrate high reliability, and best-practice guidance is required. The MetAbolomics ring Trial for CHemical groupING (MATCHING), comprising six industrial, government and academic ring-trial partners, evaluated inter-laboratory reproducibility and worked towards best-practice. An independent team selected eight substances (WY-14643, 4-chloro-3-nitroaniline, 17α-methyl-testosterone, trenbolone, aniline, dichlorprop-p, 2-chloroaniline, fenofibrate); ring-trial partners were blinded to their identities and modes-of-action. Plasma samples were derived from 28-day rat tests (two doses per substance), aliquoted, and distributed to partners. Each partner applied their preferred liquid chromatography-mass spectrometry (LC-MS) metabolomics workflows to acquire, process, quality assess, statistically analyze and report their grouping results to the European Chemicals Agency, to ensure the blinding conditions of the ring trial. Five of six partners, whose metabolomics datasets passed quality control, correctly identified the grouping of eight test substances into three categories, for both male and female rats. Strikingly, this was achieved even though a range of metabolomics approaches were used. Through assessing intrastudy quality-control samples, the sixth partner observed high technical variation and was unable to group the substances. By comparing workflows, we conclude that some heterogeneity in metabolomics methods is not detrimental to consistent grouping, and that assessing data quality prior to grouping is essential. We recommend development of international guidance for quality-control acceptance criteria. This study demonstrates the reliability of metabolomics for chemical grouping and works towards best-practice.

Mass spectrometry imaging highlights dynamic patterns of lipid co-expression with Aβ plaques in mouse and human brains

Lipids play crucial roles in the susceptibility and brain cellular responses to Alzheimer's disease (AD) and are increasingly considered potential soluble biomarkers in cerebrospinal fluid (CSF) and plasma. To delineate the pathological correlations of distinct lipid species, we conducted a comprehensive characterization of both spatially localized and global differences in brain lipid composition in AppNL-G-F mice with spatial and bulk mass spectrometry lipidomic profiling, using human amyloid-expressing (h-Aβ) and WT mouse brains controls. We observed age-dependent increases in lysophospholipids, bis(monoacylglycerol) phosphates, and phosphatidylglycerols around Aβ plaques in AppNL-G-F mice. Immunohistology-based co-localization identified associations between focal pro-inflammatory lipids, glial activation, and autophagic flux disruption. Likewise, in human donors with varying Braak stages, similar studies of cortical sections revealed co-expression of lysophospholipids and ceramides around Aβ plaques in AD (Braak stage V/VI) but not in earlier Braak stage controls. Our findings in mice provide evidence of temporally and spatially heterogeneous differences in lipid composition as local and global Aβ-related pathologies evolve. Observing similar lipidomic changes associated with pathological Aβ plaques in human AD tissue provides a foundation for understanding differences in CSF lipids with reported clinical stage or disease severity.

Inulin-type fructans and 2'fucosyllactose alter both microbial composition and appear to alleviate stress-induced mood state in a working population compared to placebo (maltodextrin): the EFFICAD Trial, a randomized, controlled trial

BACKGROUND

There is increasing interest in the bidirectional relationship existing between the gut and brain and the effects of both oligofructose and 2'fucosyllactose to alter microbial composition and mood state. Yet, much remains unknown about the ability of oligofructose and 2'fucosyllactose to improve mood state via targeted manipulation of the gut microbiota.

OBJECTIVES

We aimed to compare the effects of oligofructose and 2'fucosyllactose alone and in combination against maltodextrin (comparator) on microbial composition and mood state in a working population.

METHODS

We conducted a 5-wk, 4-arm, parallel, double-blind, randomized, placebo-controlled trial in 92 healthy adults with mild-to-moderate levels of anxiety and depression. Subjects were randomized to oligofructose 8 g/d (plus 2 g/d maltodextrin); maltodextrin 10 g/d; oligofructose 8 g/d plus 2'fucosyllactose (2 g/d) or 2'fucosyllactose 2 g/d (plus 8 g/d maltodextrin). Changes in microbial load (fluorescence in situ hybridization-flow cytometry) and composition (16S ribosomal RNA sequencing) were the primary outcomes. Secondary outcomes included gastrointestinal sensations, bowel habits, and mood state parameters.

RESULTS

There were significant increases in several bacterial taxa including Bifidobacterium, Bacteroides, Roseburia, and Faecalibacterium prausnitzii in both the oligofructose and oligofructose/2'fucosyllactose interventions (all P ≤ 0.05). Changes in bacterial taxa were highly heterogenous upon 2'fuscoyllactose supplementation. Significant improvements in Beck Depression Inventory, State Trait Anxiety Inventory Y1 and Y2, and Positive and Negative Affect Schedule scores and cortisol awakening response were detected across oligofructose, 2'fucosyllactose, and oligofructose/2'fucosyllactose combination interventions (all P ≤ 0.05). Both sole oligofructose and oligofructose/2'fuscosyllactose combination interventions outperformed both sole 2'fucosyllactose and maltodextrin in improvements in several mood state parameters (all P ≤ 0.05).

CONCLUSION

The results of this study indicate that oligofructose and combination of oligofructose/2'fucosyllactose can beneficially alter microbial composition along with improving mood state parameters. Future work is needed to understand key microbial differences separating individual responses to 2'fucosyllactose supplementation. This trial was registered at clinicaltrials.gov as NCT05212545.

Fecal microbiota transplantation plus anti-PD-1 immunotherapy in advanced melanoma: a phase I trial

Fecal microbiota transplantation (FMT) represents a potential strategy to overcome resistance to immune checkpoint inhibitors in patients with refractory melanoma; however, the role of FMT in first-line treatment settings has not been evaluated. We conducted a multicenter phase I trial combining healthy donor FMT with the PD-1 inhibitors nivolumab or pembrolizumab in 20 previously untreated patients with advanced melanoma. The primary end point was safety. No grade 3 adverse events were reported from FMT alone. Five patients (25%) experienced grade 3 immune-related adverse events from combination therapy. Key secondary end points were objective response rate, changes in gut microbiome composition and systemic immune and metabolomics analyses. The objective response rate was 65% (13 of 20), including four (20%) complete responses. Longitudinal microbiome profiling revealed that all patients engrafted strains from their respective donors; however, the acquired similarity between donor and patient microbiomes only increased over time in responders. Responders experienced an enrichment of immunogenic and a loss of deleterious bacteria following FMT. Avatar mouse models confirmed the role of healthy donor feces in increasing anti-PD-1 efficacy. Our results show that FMT from healthy donors is safe in the first-line setting and warrants further investigation in combination with immune checkpoint inhibitors. ClinicalTrials.gov identifier NCT03772899 .

Pathobionts in the tumour microbiota predict survival following resection for colorectal cancer

BACKGROUND AND AIMS

The gut microbiota is implicated in the pathogenesis of colorectal cancer (CRC). We aimed to map the CRC mucosal microbiota and metabolome and define the influence of the tumoral microbiota on oncological outcomes.

METHODS

A multicentre, prospective observational study was conducted of CRC patients undergoing primary surgical resection in the UK (n = 74) and Czech Republic (n = 61). Analysis was performed using metataxonomics, ultra-performance liquid chromatography-mass spectrometry (UPLC-MS), targeted bacterial qPCR and tumour exome sequencing. Hierarchical clustering accounting for clinical and oncological covariates was performed to identify clusters of bacteria and metabolites linked to CRC. Cox proportional hazards regression was used to ascertain clusters associated with disease-free survival over median follow-up of 50 months.

RESULTS

Thirteen mucosal microbiota clusters were identified, of which five were significantly different between tumour and paired normal mucosa. Cluster 7, containing the pathobionts Fusobacterium nucleatum and Granulicatella adiacens, was strongly associated with CRC (P_FDR = 0.0002). Additionally, tumoral dominance of cluster 7 independently predicted favourable disease-free survival (adjusted p = 0.031). Cluster 1, containing Faecalibacterium prausnitzii and Ruminococcus gnavus, was negatively associated with cancer (P_FDR = 0.0009), and abundance was independently predictive of worse disease-free survival (adjusted p = 0.0009). UPLC-MS analysis revealed two major metabolic (Met) clusters. Met 1, composed of medium chain (MCFA), long-chain (LCFA) and very long-chain (VLCFA) fatty acid species, ceramides and lysophospholipids, was negatively associated with CRC (P_FDR = 2.61 × 10^-11); Met 2, composed of phosphatidylcholine species, nucleosides and amino acids, was strongly associated with CRC (P_FDR = 1.30 × 10^-12), but metabolite clusters were not associated with disease-free survival (p = 0.358). An association was identified between Met 1 and DNA mismatch-repair deficiency (p = 0.005). FBXW7 mutations were only found in cancers predominant in microbiota cluster 7.

CONCLUSIONS

Networks of pathobionts in the tumour mucosal niche are associated with tumour mutation and metabolic subtypes and predict favourable outcome following CRC resection. Video Abstract.

Discriminatory Changes in Circulating Metabolites as a Predictor of Hepatocellular Cancer in Patients with Metabolic (Dysfunction) Associated Fatty Liver Disease

Introduction

The burden of metabolic (dysfunction) associated fatty liver disease (MAFLD) is rising mirrored by an increase in hepatocellular cancer (HCC). MAFLD and its sequelae are characterized by perturbations in lipid handling, inflammation, and mitochondrial damage. The profile of circulating lipid and small molecule metabolites with the development of HCC is poorly characterized in MAFLD and could be used in future studies as a biomarker for HCC.

Methods

We assessed the profile of 273 lipid and small molecule metabolites by ultra-performance liquid chromatography coupled to high-resolution mass spectrometry in serum from patients with MAFLD (n = 113) and MAFLD-associated HCC (n = 144) from six different centers. Regression models were used to identify a predictive model of HCC.

Results

Twenty lipid species and one metabolite, reflecting changes in mitochondrial function and sphingolipid metabolism, were associated with the presence of cancer on a background of MAFLD with high accuracy (AUC 0.789, 95% CI: 0.721-0.858), which was enhanced with the addition of cirrhosis to the model (AUC 0.855, 95% CI: 0.793-0.917). In particular, the presence of these metabolites was associated with cirrhosis in the MAFLD subgroup (p < 0.001). When considering the HCC cohort alone, the metabolic signature was an independent predictor of overall survival (HR 1.42, 95% CI: 1.09-1.83, p < 0.01).

Conclusion

These exploratory findings reveal a metabolic signature in serum which is capable of accurately detecting the presence of HCC on a background of MAFLD. This unique serum signature will be taken forward for further investigation of diagnostic performance as biomarker of early stage HCC in patients with MAFLD in the future.

The gut microbiota and metabolome are associated with diminished COVID-19 vaccine-induced antibody responses in immunosuppressed inflammatory bowel disease patients

BACKGROUND

Patients with inflammatory bowel disease (IBD) treated with anti-TNF therapy exhibit attenuated humoral immune responses to vaccination against SARS-CoV-2. The gut microbiota and its functional metabolic output, which are perturbed in IBD, play an important role in shaping host immune responses. We explored whether the gut microbiota and metabolome could explain variation in anti-SARS-CoV-2 vaccination responses in immunosuppressed IBD patients.

METHODS

Faecal and serum samples were prospectively collected from infliximab-treated patients with IBD in the CLARITY-IBD study undergoing vaccination against SARS-CoV-2. Antibody responses were measured following two doses of either ChAdOx1 nCoV-19 or BNT162b2 vaccine. Patients were classified as having responses above or below the geometric mean of the wider CLARITY-IBD cohort. 16S rRNA gene amplicon sequencing, nuclear magnetic resonance (NMR) spectroscopy and bile acid profiling with ultra-high-performance liquid chromatography mass spectrometry (UHPLC-MS) were performed on faecal samples. Univariate, multivariable and correlation analyses were performed to determine gut microbial and metabolomic predictors of response to vaccination.

FINDINGS

Forty-three infliximab-treated patients with IBD were recruited (30 Crohn's disease, 12 ulcerative colitis, 1 IBD-unclassified; 26 with concomitant thiopurine therapy). Eight patients had evidence of prior SARS-CoV-2 infection. Seventeen patients (39.5%) had a serological response below the geometric mean. Gut microbiota diversity was lower in below average responders (p = 0.037). Bilophila abundance was associated with better serological response, while Streptococcus was associated with poorer response. The faecal metabolome was distinct between above and below average responders (OPLS-DA R²X 0.25, R²Y 0.26, Q² 0.15; CV-ANOVA p = 0.038). Trimethylamine, isobutyrate and omega-muricholic acid were associated with better response, while succinate, phenylalanine, taurolithocholate and taurodeoxycholate were associated with poorer response.

INTERPRETATION

Our data suggest that there is an association between the gut microbiota and variable serological response to vaccination against SARS-CoV-2 in immunocompromised patients. Microbial metabolites including trimethylamine may be important in mitigating anti-TNF-induced attenuation of the immune response.

FUNDING

JLA is the recipient of an NIHR Academic Clinical Lectureship (CL-2019-21-502), funded by Imperial College London and The Joyce and Norman Freed Charitable Trust. BHM is the recipient of an NIHR Academic Clinical Lectureship (CL-2019-21-002). The Division of Digestive Diseases at Imperial College London receives financial and infrastructure support from the NIHR Imperial Biomedical Research Centre (BRC) based at Imperial College Healthcare NHS Trust and Imperial College London. Metabolomics studies were performed at the MRC-NIHR National Phenome Centre at Imperial College London; this work was supported by the Medical Research Council (MRC), the National Institute of Health Research (NIHR) (grant number MC_PC_12025) and infrastructure support was provided by the NIHR Imperial Biomedical Research Centre (BRC). The NIHR Exeter Clinical Research Facility is a partnership between the University of Exeter Medical School College of Medicine and Health, and Royal Devon and Exeter NHS Foundation Trust. This project is supported by the National Institute for Health Research (NIHR) Exeter Clinical Research Facility. The views expressed are those of the authors and not necessarily those of the NIHR or the UK Department of Health and Social Care.

Assessing the clinical value of faecal bile acid profiling to predict recurrence in primary Clostridioides difficile infection

BACKGROUND

Factors influencing recurrence risk in primary Clostridioides difficile infection (CDI) are poorly understood, and tools predicting recurrence are lacking. Perturbations in bile acids (BAs) contribute to CDI pathogenesis and may be relevant to primary disease prognosis.

AIMS

To define stool BA dynamics in patients with primary CDI and to explore signatures predicting recurrence METHODS: Weekly stool samples were collected from patients with primary CDI from the last day of anti-CDI therapy until recurrence or, otherwise, through 8 weeks post-completion. Ultra-high performance liquid chromatography-mass spectrometry was used to profile BAs. Stool bile salt hydrolase (BSH) activity was measured to determine primary BA bacterial deconjugation capacity. Multivariate and univariate models were used to define differential BA trajectories in patients with recurrence versus those without, and to assess faecal BAs as predictive markers for recurrence.

RESULTS

Twenty (36%) of 56 patients (median age: 57, 64% male) had recurrence; 80% of recurrences occurred within the first 9 days post-antibiotic treatment. Principal component analysis of stool BA profiles demonstrated clustering by recurrence status and post-treatment timepoint. Longitudinal faecal BA trajectories showed recovery of secondary BAs and their derivatives only in patients without recurrence. BSH activity increased over time only among non-relapsing patients (β = 0.056; likelihood ratio test p = 0.018). A joint longitudinal-survival model identified five stool BAs with area under the receiver operating characteristic curve >0.73 for predicting recurrence within 9 days post-CDI treatment.

CONCLUSIONS

Gut BA metabolism dynamics differ in primary CDI patients between those developing recurrence and those who do not. Individual BAs show promise as potential novel biomarkers to predict CDI recurrence.

Diet Patterns Are Associated with Circulating Metabolites and Lipid Profiles of South Asians in the United States

BACKGROUND

South Asians are at higher risk for cardiometabolic disease than many other racial/ethnic minority groups. Diet patterns in US South Asians have unique components associated with cardiometabolic disease.

OBJECTIVES

We aimed to characterize the metabolites associated with 3 representative diet patterns.

METHODS

We included 722 participants in the Mediators of Atherosclerosis in South Asians Living in America (MASALA) cohort study aged 40-84 y without known cardiovascular disease. Fasting serum specimens and diet and demographic questionnaires were collected at baseline and diet patterns previously generated through principal components analysis. LC-MS-based untargeted metabolomic and lipidomic analysis was conducted with targeted integration of known metabolite and lipid signals. Linear regression models of diet pattern factor score and log-transformed metabolites adjusted for age, sex, caloric intake, and BMI and adjusted for multiple comparisons were performed, followed by elastic net linear regression of significant metabolites.

RESULTS

There were 443 metabolites of known identity extracted from the profiling data. The "animal protein" diet pattern was associated with 61 metabolites and lipids, including glycerophospholipids phosphatidylethanolamine PE(O-16:1/20:4) and/or PE(P-16:0/20:4) (β: 0.13; 95% CI: 0.11, 0.14) and N-acyl phosphatidylethanolamines (NAPEs) NAPE(O-18:1/20:4/18:0) and/or NAPE(P-18:0/20:4/18:0) (β: 0.13; 95% CI: 0.11, 0.14), lysophosphatidylinositol (LPI) (22:6/0:0) (β: 0.14; 95% CI: 0.12, 0.17), and fatty acid (FA) (22:6) (β: 0.15; 95% CI: 0.13, 0.17). The "fried snacks, sweets, high-fat dairy" pattern was associated with 12 lipids, including PC(16:0/22:6) (β: -0.08; 95% CI: -0.09, -0.06) and FA (22:6) (β: 0.14; 95% CI: -0.17, -0.10). The "fruits, vegetables, nuts, and legumes" pattern was associated with 5 metabolites including proline betaine (β: 0.17; 95% CI: 0.09, 0.25) (P < 0.0002).

CONCLUSIONS

Three predominant dietary patterns in US South Asians are associated with circulating metabolites differentiated by lipids including glycerophospholipids and PUFAs and the amino acid proline betaine.

Metabolome-wide association study on ABCA7 indicates a role of ceramide metabolism in Alzheimer's disease

Genome-wide association studies (GWASs) have identified genetic loci associated with the risk of Alzheimer's disease (AD), but the molecular mechanisms by which they confer risk are largely unknown. We conducted a metabolome-wide association study (MWAS) of AD-associated loci from GWASs using untargeted metabolic profiling (metabolomics) by ultraperformance liquid chromatography-mass spectrometry (UPLC-MS). We identified an association of lactosylceramides (LacCer) with AD-related single-nucleotide polymorphisms (SNPs) in ABCA7 (P = 5.0 × 10^-5 to 1.3 × 10^-44). We showed that plasma LacCer concentrations are associated with cognitive performance and genetically modified levels of LacCer are associated with AD risk. We then showed that concentrations of sphingomyelins, ceramides, and hexosylceramides were altered in brain tissue from Abca7 knockout mice, compared with wild type (WT) (P = 0.049-1.4 × 10^-5), but not in a mouse model of amyloidosis. Furthermore, activation of microglia increases intracellular concentrations of hexosylceramides in part through induction in the expression of sphingosine kinase, an enzyme with a high control coefficient for sphingolipid and ceramide synthesis. Our work suggests that the risk for AD arising from functional variations in ABCA7 is mediated at least in part through ceramides. Modulation of their metabolism or downstream signaling may offer new therapeutic opportunities for AD.

Metabolomic profiling in small vessel disease identifies multiple associations with disease severity

Cerebral small vessel disease is a major cause of vascular cognitive impairment and dementia. There are few treatments, largely reflecting limited understanding of the underlying pathophysiology. Metabolomics can be used to identify novel risk factors to better understand pathogenesis and to predict disease progression and severity. We analysed data from 624 patients with symptomatic cerebral small vessel disease from two prospective cohort studies. Serum samples were collected at baseline and patients underwent MRI scans and cognitive testing at regular intervals with up to 14 years of follow-up. Using ultra-performance liquid chromatography-mass spectrometry and nuclear magnetic resonance spectroscopy, we obtained metabolic and lipidomic profiles from 369 annotated metabolites and 54 764 unannotated features and examined their association with respect to disease severity, assessed using MRI small vessel disease markers, cognition and future risk of all-cause dementia. Our analysis identified 28 metabolites that were significantly associated with small vessel disease imaging markers and cognition. Decreased levels of multiple glycerophospholipids and sphingolipids were associated with increased small vessel disease load as evidenced by higher white matter hyperintensity volume, lower mean diffusivity normalized peak height, greater brain atrophy and impaired cognition. Higher levels of creatine, FA(18:2(OH)) and SM(d18:2/24:1) were associated with increased lacune count, higher white matter hyperintensity volume and impaired cognition. Lower baseline levels of carnitines and creatinine were associated with higher annualized change in peak width of skeletonized mean diffusivity, and 25 metabolites, including lipoprotein subclasses, amino acids and xenobiotics, were associated with future dementia incidence. Our results show multiple distinct metabolic signatures that are associated with imaging markers of small vessel disease, cognition and conversion to dementia. Further research should assess causality and the use of metabolomic screening to improve the ability to predict future disease severity and dementia risk in small vessel disease. The metabolomic profiles may also provide novel insights into disease pathogenesis and help identify novel treatment approaches.

Associations of NAFLD with circulating ceramides and impaired glycemia

AIM

Determine the association of circulating ceramides with NAFLD and glycemic impairment.

METHODS

Sample: 669 participants in the Mediators of Atherosclerosis in South Asians Living in America (MASALA) cohort aged 40-84 years without cardiovascular disease, cirrhosis, or significant alcohol intake.

CLINICAL MEASURES

Computed tomography scans at baseline for hepatic attenuation. Fasting serum specimens at baseline and after 5 years. Lipidomics: LC-MS-based analysis of 19 known ceramide signals.

STATISTICAL ANALYSIS

Linear and logistic regression models of log-transformed ceramides, hepatic attenuation and glucose adjusted for age, sex, calories, study site, BMI, exercise, diet quality, alcohol, saturated fat, lipid-lowering medications and fasting glucose.

RESULTS

Average age was 55 years, 44% were women, mean BMI was 25.9 kg/m2, and 8% had NAFLD. In adjusted models, Cer(d16:1/20:0) and Cer(d18:1/18:0) were associated with lower mean hepatic attenuation (increased liver fat) (β -4.29; 95% CI [-5.98, -2.59]) and (β -3.40; 95% CI [-5.11, -1.70]), and LacCer(d18:1/16:0) with higher attenuation (β 4.44; 95% CI [2.15, 6.73]). All three ceramides partially mediated the relationship between hepatic attenuation and fasting glucose by 16%, 11% and 5%, respectively, after 5-years.

CONCLUSIONS

Three circulating ceramides were strongly associated with NAFLD and fasting glucose after 5 years, and partially mediated this association.

Antiviral metabolite 3'-deoxy-3',4'-didehydro-cytidine is detectable in serum and identifies acute viral infections including COVID-19

BACKGROUND

There is a critical need for rapid viral infection diagnostics to enable prompt case identification in pandemic settings and support targeted antimicrobial prescribing.

METHODS

Using untargeted high-resolution liquid chromatography coupled with mass spectrometry, we compared the admission serum metabolome of emergency department patients with viral infections (including COVID-19), bacterial infections, inflammatory conditions, and healthy controls. Sera from an independent cohort of emergency department patients admitted with viral or bacterial infections underwent profiling to validate findings. Associations between whole-blood gene expression and the identified metabolite of interest were examined.

FINDINGS

3'-Deoxy-3',4'-didehydro-cytidine (ddhC), a free base of the only known human antiviral small molecule ddhC-triphosphate (ddhCTP), was detected for the first time in serum. When comparing 60 viral with 101 non-viral cases in the discovery cohort, ddhC was the most significantly differentially abundant metabolite, generating an area under the receiver operating characteristic curve (AUC) of 0.954 (95% CI: 0.923-0.986). In the validation cohort, ddhC was again the most significantly differentially abundant metabolite when comparing 40 viral with 40 bacterial cases, generating an AUC of 0.81 (95% CI 0.708-0.915). Transcripts of viperin and CMPK2, enzymes responsible for ddhCTP synthesis, were among the five genes most highly correlated with ddhC abundance.

CONCLUSIONS

The antiviral precursor molecule ddhC is detectable in serum and an accurate marker for acute viral infection. Interferon-inducible genes viperin and CMPK2 are implicated in ddhC production in vivo. These findings highlight a future diagnostic role for ddhC in viral diagnosis, pandemic preparedness, and acute infection management.

FUNDING

NIHR Imperial BRC; UKRI.

Machine Learning at the Interface of Polymer Science and Biology: How Far Can We Go?

This Perspective outlines recent progress and future directions for using machine learning (ML), a data-driven method, to address critical questions in the design, synthesis, processing, and characterization of biomacromolecules. The achievement of these tasks requires the navigation of vast and complex chemical and biological spaces, difficult to accomplish with reasonable speed. Using modern algorithms and supercomputers, quantum physics methods are able to examine systems containing a few hundred interacting species and determine the probability of finding them in a particular region of phase space, thereby anticipating their properties. Likewise, modern approaches in chemistry and biomolecular simulation, supported by high performance computing, have culminated in producing data sets of escalating size and intrinsically high complexity. Hence, using ML to extract relevant information from these fields is of paramount importance to advance our understanding of chemical and biomolecular systems. At the heart of ML approaches lie statistical algorithms, which by evaluating a portion of a given data set, identify, learn, and manipulate the underlying rules that govern the whole data set. The assembly of a quality model to represent the data followed by the predictions and elimination of error sources are the key steps in ML. In addition to a growing infrastructure of ML tools to address complex problems, an increasing number of aspects related to our understanding of the fundamental properties of biomacromolecules are exposed to ML. These fields, including those residing at the interface of polymer science and biology (i.e., structure determination, de novo design, folding, and dynamics), strive to adopt and take advantage of the transformative power offered by approaches in the ML domain, which clearly has the potential of accelerating research in the field of biomacromolecules.

Direct on-swab metabolic profiling of vaginal microbiome host interactions during pregnancy and preterm birth

The pregnancy vaginal microbiome contributes to risk of preterm birth, the primary cause of death in children under 5 years of age. Here we describe direct on-swab metabolic profiling by Desorption Electrospray Ionization Mass Spectrometry (DESI-MS) for sample preparation-free characterisation of the cervicovaginal metabolome in two independent pregnancy cohorts (VMET, n = 160; 455 swabs; VMET II, n = 205; 573 swabs). By integrating metataxonomics and immune profiling data from matched samples, we show that specific metabolome signatures can be used to robustly predict simultaneously both the composition of the vaginal microbiome and host inflammatory status. In these patients, vaginal microbiota instability and innate immune activation, as predicted using DESI-MS, associated with preterm birth, including in women receiving cervical cerclage for preterm birth prevention. These findings highlight direct on-swab metabolic profiling by DESI-MS as an innovative approach for preterm birth risk stratification through rapid assessment of vaginal microbiota-host dynamics.

Integration of global metabolomics and lipidomics approaches reveals the molecular mechanisms and the potential biomarkers for postoperative recurrence in early-stage cholangiocarcinoma

Background

Cholangiocarcioma (CCA) treatment is challenging because most of the patients are diagnosed when the disease is advanced, and cancer recurrence is the main problem after treatment, leading to low survival rates. Therefore, our understanding of the mechanism underlying CCA recurrence is essential in order to prevent CCA recurrence and improve patient outcomes.

Methods

We performed ¹H-NMR and UPLC-MS-based metabolomics on the CCA serum. The differential metabolites were further analyzed using pathway analysis and potential biomarker identification.

Results

At an early stage, the metabolites involved in energy metabolisms, such as pyruvate metabolism, and the TCA cycle, are downregulated, while most lipids, including TGs, PCs, PEs, and PAs, are upregulated in recurrence patients. This metabolic feature has been described in cancer stem-like cell (CSC) metabolism. In addition, the CSC markers CD44v6 and CD44v8-10 are associated with CD36 (a protein involved in lipid uptake) as well as with recurrence-free survival. We also found that citrate, sarcosine, succinate, creatine, creatinine and pyruvate, and TGs have good predictive values for CCA recurrence.

Conclusion

Our study demonstrates the possible molecular mechanisms underlying CCA recurrence, and these may associate with the existence of CSCs. The metabolic change involved in the recurrence pathway might be used to determine biomarkers for predicting CCA recurrence.

Supplementary Information

The online version contains supplementary material available at 10.1186/s40170-021-00266-5.

Statistical analysis in metabolic phenotyping

Metabolic phenotyping is an important tool in translational biomedical research. The advanced analytical technologies commonly used for phenotyping, including mass spectrometry (MS) and nuclear magnetic resonance (NMR) spectroscopy, generate complex data requiring tailored statistical analysis methods. Detailed protocols have been published for data acquisition by liquid NMR, solid-state NMR, ultra-performance liquid chromatography (LC-)MS and gas chromatography (GC-)MS on biofluids or tissues and their preprocessing. Here we propose an efficient protocol (guidelines and software) for statistical analysis of metabolic data generated by these methods. Code for all steps is provided, and no prior coding skill is necessary. We offer efficient solutions for the different steps required within the complete phenotyping data analytics workflow: scaling, normalization, outlier detection, multivariate analysis to explore and model study-related effects, selection of candidate biomarkers, validation, multiple testing correction and performance evaluation of statistical models. We also provide a statistical power calculation algorithm and safeguards to ensure robust and meaningful experimental designs that deliver reliable results. We exemplify the protocol with a two-group classification study and data from an epidemiological cohort; however, the protocol can be easily modified to cover a wider range of experimental designs or incorporate different modeling approaches. This protocol describes a minimal set of analyses needed to rigorously investigate typical datasets encountered in metabolic phenotyping.

Prediction of response of methotrexate in patients with rheumatoid arthritis using serum lipidomics

Methotrexate (MTX) is a common first-line treatment for new-onset rheumatoid arthritis (RA). However, MTX is ineffective for 30–40% of patients and there is no way to know which patients might benefit. Here, we built statistical models based on serum lipid levels measured at two time-points (pre-treatment and following 4 weeks on-drug) to investigate if MTX response (by 6 months) could be predicted. Patients about to commence MTX treatment for the first time were selected from the Rheumatoid Arthritis Medication Study (RAMS). Patients were categorised as good or non-responders following 6 months on-drug using EULAR response criteria. Serum lipids were measured using ultra‐performance liquid chromatography–mass spectrometry and supervised machine learning methods (including regularized regression, support vector machine and random forest) were used to predict EULAR response. Models including lipid levels were compared to models including clinical covariates alone. The best performing classifier including lipid levels (assessed at 4 weeks) was constructed using regularized regression (ROC AUC 0.61 ± 0.02). However, the clinical covariate based model outperformed the classifier including lipid levels when either pre- or on-treatment time-points were investigated (ROC AUC 0.68 ± 0.02). Pre- or early-treatment serum lipid profiles are unlikely to inform classification of MTX response by 6 months with performance adequate for use in RA clinical management.

Circulating metabolites and lipids are associated with glycaemic measures in South Asians

BACKGROUND

South Asians are at higher risk for diabetes (DM) than many other racial/ethnic groups. Circulating metabolites are measurable products of metabolic processes that may explain the aetiology of elevated risk. We characterized metabolites associated with prevalent DM and glycaemic measures in South Asians.

METHODS

We included 717 participants from the Mediators of Atherosclerosis in South Asians Living in America (MASALA) study, aged 40-84 years. We used baseline fasting serum for metabolomics and demographic, behavioural, glycaemic data from baseline and at 5 years. We performed LC-MS untargeted metabolomic and lipidomic analysis with targeted integration of known signals. Individual linear and ordinal logistic regression models were adjusted for age, sex, BMI, diet, exercise, alcohol, smoking and family history of DM followed by elastic net regression to identify metabolites most associated with the outcome.

RESULTS

There were 258 metabolites with detectable signal in >98% of samples. Thirty-four metabolites were associated with prevalent DM in an elastic net model. Predominant metabolites associated with DM were sphingomyelins, proline (OR 15.86; 95% CI 4.72, 53.31) and betaine (OR 0.03; 0.004, 0.14). Baseline tri- and di-acylglycerols [DG (18:0/16:0) (18.36; 11.79, 24.92)] were positively associated with fasting glucose and long-chain acylcarnitines [CAR 26:1 (-0.40; -0.54, -0.27)] were inversely associated with prevalent DM and HbA_1c at follow-up.

DISCUSSION

A metabolomic signature in South Asians may help determine the unique aetiology of diabetes in this high-risk ethnic group. Future work will externally validate our findings and determine the effects of modifiable risk factors for DM.

Representing the Metabolome with High Fidelity: Range and Response as Quality Control Factors in LC-MS-Based Global Profiling

Liquid chromatography-mass spectrometry (LC-MS) is a powerful and widely used technique for measuring the abundance of chemical species in living systems. Its sensitivity, analytical specificity, and direct applicability to biofluids and tissue extracts impart great promise for the discovery and mechanistic characterization of biomarker panels for disease detection, health monitoring, patient stratification, and treatment personalization. Global metabolic profiling applications yield complex data sets consisting of multiple feature measurements for each chemical species observed. While this multiplicity can be useful in deriving enhanced analytical specificity and chemical identities from LC-MS data, data set inflation and quantitative imprecision among related features is problematic for statistical analyses and interpretation. This Perspective provides a critical evaluation of global profiling data fidelity with respect to measurement linearity and the quantitative response variation observed among components of the spectra. These elements of data quality are widely overlooked in untargeted metabolomics yet essential for the generation of data that accurately reflect the metabolome. Advanced feature filtering informed by linear range estimation and analyte response factor assessment is advocated as an attainable means of controlling LC-MS data quality in global profiling studies and exemplified herein at both the feature and data set level.

Urinary metabolic phenotyping for Alzheimer's disease

Finding early disease markers using non-invasive and widely available methods is essential to develop a successful therapy for Alzheimer’s Disease. Few studies to date have examined urine, the most readily available biofluid. Here we report the largest study to date using comprehensive metabolic phenotyping platforms (NMR spectroscopy and UHPLC-MS) to probe the urinary metabolome in-depth in people with Alzheimer’s Disease and Mild Cognitive Impairment. Feature reduction was performed using metabolomic Quantitative Trait Loci, resulting in the list of metabolites associated with the genetic variants. This approach helps accuracy in identification of disease states and provides a route to a plausible mechanistic link to pathological processes. Using these mQTLs we built a Random Forests model, which not only correctly discriminates between people with Alzheimer’s Disease and age-matched controls, but also between individuals with Mild Cognitive Impairment who were later diagnosed with Alzheimer’s Disease and those who were not. Further annotation of top-ranking metabolic features nominated by the trained model revealed the involvement of cholesterol-derived metabolites and small-molecules that were linked to Alzheimer’s pathology in previous studies.

A Python-Based Pipeline for Preprocessing LC-MS Data for Untargeted Metabolomics Workflows

Preprocessing data in a reproducible and robust way is one of the current challenges in untargeted metabolomics workflows. Data curation in liquid chromatography–mass spectrometry (LC–MS) involves the removal of biologically non-relevant features (retention time, m/z pairs) to retain only high-quality data for subsequent analysis and interpretation. The present work introduces TidyMS, a package for the Python programming language for preprocessing LC–MS data for quality control (QC) procedures in untargeted metabolomics workflows. It is a versatile strategy that can be customized or fit for purpose according to the specific metabolomics application. It allows performing quality control procedures to ensure accuracy and reliability in LC–MS measurements, and it allows preprocessing metabolomics data to obtain cleaned matrices for subsequent statistical analysis. The capabilities of the package are shown with pipelines for an LC–MS system suitability check, system conditioning, signal drift evaluation, and data curation. These applications were implemented to preprocess data corresponding to a new suite of candidate plasma reference materials developed by the National Institute of Standards and Technology (NIST; hypertriglyceridemic, diabetic, and African-American plasma pools) to be used in untargeted metabolomics studies in addition to NIST SRM 1950 Metabolites in Frozen Human Plasma. The package offers a rapid and reproducible workflow that can be used in an automated or semi-automated fashion, and it is an open and free tool available to all users.

SMolESY: an efficient and quantitative alternative to on-instrument macromolecular 1H-NMR signal suppression

One-dimensional (1D) proton-nuclear magnetic resonance (1H-NMR) spectroscopy is an established technique for measuring small molecules in a wide variety of complex biological sample types. It is demonstrably reproducible, easily automatable and consequently ideal for routine and large-scale application. However, samples containing proteins, lipids, polysaccharides and other macromolecules produce broad signals which overlap and convolute those from small molecules. NMR experiment types designed to suppress macromolecular signals during acquisition may be additionally performed, however these approaches add to the overall sample analysis time and cost, especially for large cohort studies, and fail to produce reliably quantitative data. Here, we propose an alternative way of computationally eliminating macromolecular signals, employing the mathematical differentiation of standard 1H-NMR spectra, producing small molecule-enhanced spectra with preserved quantitative capability and increased resolution. Our approach, presented in its simplest form, was implemented in a cheminformatic toolbox and successfully applied to more than 3000 samples of various biological matrices rich or potentially rich with macromolecules, offering an efficient alternative to on-instrument experimentation, facilitating NMR use in routine and large-scale applications.