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.

Mapping of 1 H NMR chemical shifts relationship with chemical similarities for the acceleration of metabolic profiling: Application on blood products.

One-dimensional (1D) proton-nuclear magnetic resonance (1 H-NMR) spectroscopy is an established technique for the deconvolution of complex biological sample types via the identification/quantification of small molecules. It is highly reproducible and could be easily automated for small to large-scale bioanalytical, epidemiological, and in general metabolomics studies. However, chemical shift variability is a serious issue that must still be solved in order to fully automate metabolite identification. Herein, we demonstrate a strategy to increase the confidence in assignments and effectively predict the chemical shifts of various NMR signals based upon the simplest form of statistical models (i.e., linear regression). To build these models, we were guided by chemical homology in serum/plasma metabolites classes (i.e., amino acids and carboxylic acids) and similarity between chemical groups such as methyl protons. Our models, built on 940 serum samples and validated in an independent cohort of 1,052 plasma-EDTA spectra, were able to successfully predict the 1 H NMR chemical shifts of 15 metabolites within ~1.5 linewidths (Δv1/2 ) error range on average. This pilot study demonstrates the potential of developing an algorithm for the accurate assignment of 1 H NMR chemical shifts based solely on chemically defined constraints.

1H NMR Signals from Urine Excreted Protein Are a Source of Bias in Probabilistic Quotient Normalization.

Normalization to account for variation in urinary dilution is crucial for interpretation of urine metabolic profiles. Probabilistic quotient normalization (PQN) is used routinely in metabolomics but is sensitive to systematic variation shared across a large proportion of the spectral profile (>50%). Where 1H nuclear magnetic resonance (NMR) spectroscopy is employed, the presence of urinary protein can elevate the spectral baseline and substantially impact the resulting profile. Using 1H NMR profile measurements of spot urine samples collected from hospitalized COVID-19 patients in the ISARIC 4C study, we determined that PQN coefficients are significantly correlated with observed protein levels (r2 = 0.423, p < 2.2 × 10-16). This correlation was significantly reduced (r2 = 0.163, p < 2.2 × 10-16) when using a computational method for suppression of macromolecular signals known as small molecule enhancement spectroscopy (SMolESY) for proteinic baseline removal prior to PQN. These results highlight proteinuria as a common yet overlooked source of bias in 1H NMR metabolic profiling studies which can be effectively mitigated using SMolESY or other macromolecular signal suppression methods before estimation of normalization coefficients.

peakPantheR, an R package for large-scale targeted extraction and integration of annotated metabolic features in LC-MS profiling datasets.

SUMMARY: Untargeted liquid chromatography-mass spectrometry (LC-MS) profiling assays are capable of measuring thousands of chemical compounds in a single sample, but unreliable feature extraction and metabolite identification remain considerable barriers to their interpretation and usefulness. peakPantheR (Peak Picking and ANnoTation of High-resolution Experiments in R) is an R package for the targeted extraction and integration of annotated features from LC-MS profiling experiments. It takes advantage of chromatographic and spectral databases and prior information of sample matrix composition to generate annotated and interpretable metabolic phenotypic datasets and power workflows for real-time data quality assessment. AVAILABILITY AND IMPLEMENTATION: peakPantheR is available via Bioconductor (https://bioconductor.org/packages/peakPantheR/). Documentation and worked examples are available at https://phenomecentre.github.io/peakPantheR.github.io/ and https://github.com/phenomecentre/metabotyping-dementia-urine. SUPPLEMENTARY INFORMATION: Supplementary data are available at Bioinformatics online.

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.

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

SUMMARY: As large-scale metabolic phenotyping studies become increasingly common, the need for systemic methods for pre-processing and quality control (QC) of analytical data prior to statistical analysis has become increasingly important, both within a study, and to allow meaningful inter-study comparisons. The nPYc-Toolbox provides software for the import, pre-processing, QC and visualization of metabolic phenotyping datasets, either interactively, or in automated pipelines. AVAILABILITY AND IMPLEMENTATION: The nPYc-Toolbox is implemented in Python, and is freely available from the Python package index https://pypi.org/project/nPYc/, source is available at https://github.com/phenomecentre/nPYc-Toolbox. Full documentation can be found at http://npyc-toolbox.readthedocs.io/ and exemplar datasets and tutorials at https://github.com/phenomecentre/nPYc-toolbox-tutorials.

Metabolic phenotyping for enhanced mechanistic stratification of chronic hepatitis C-induced liver fibrosis.

OBJECTIVES: The invasive nature of biopsy alongside issues with categorical staging and sampling error has driven research into noninvasive biomarkers for the assessment of liver fibrosis in order to stratify and personalize treatment of patients with liver disease. Here, we sought to determine whether a metabonomic approach could be used to identify signatures reflective of the dynamic, pathological metabolic perturbations associated with fibrosis in chronic hepatitis C (CHC) patients. METHODS: Plasma nuclear magnetic resonance (NMR) spectral profiles were generated for two independent cohorts of CHC patients and healthy controls (n=50 original and n=63 validation). Spectral data were analyzed and significant discriminant biomarkers associated with fibrosis (as graded by enhanced liver fibrosis (ELF) and METAVIR scores) identified using orthogonal projection to latent structures (O-PLS). RESULTS: Increased severity of fibrosis was associated with higher tyrosine, phenylalanine, methionine, citrate and, very-low-density lipoprotein (vLDL) and lower creatine, low-density lipoprotein (LDL), phosphatidylcholine, and N-Acetyl-α1-acid-glycoprotein. Although area under the receiver operator characteristic curve analysis revealed a high predictive performance for classification based on METAVIR-derived models, <40% of identified biomarkers were validated in the second cohort. In the ELF-derived models, however, over 80% of the biomarkers were validated. CONCLUSIONS: Our findings suggest that modeling against a continuous ELF-derived score of fibrosis provides a more robust assessment of the metabolic changes associated with fibrosis than modeling against the categorical METAVIR score. Plasma metabolic phenotypes reflective of CHC-induced fibrosis primarily define alterations in amino-acid and lipid metabolism, and hence identify mechanistically relevant pathways for further investigation as therapeutic targets.

Comparison of phenomics and cfDNA in a large breast screening population: the Breast Screening and Monitoring Study (BSMS).

To assess their roles in breast cancer diagnostics, we aimed to compare plasma cell-free DNA (cfDNA) levels with the circulating metabolome in a large breast screening cohort of women recalled for mammography, including healthy women and women with mammographically detected breast diseases, ductal carcinoma in situ and invasive breast cancer: the Breast Screening and Monitoring Study (BSMS). In 999 women, plasma was analyzed by nuclear magnetic resonance (NMR) and Ultra-Performance Liquid Chromatography-Mass Spectrometry (UPLC-MS) and then processed to isolate and quantify total cfDNA. NMR and UPLC-MS results were compared with data for 186 healthy women derived from the AIRWAVE cohort. Results showed no significant differences between groups for all metabolites, whereas invasive cancers had significantly higher plasma cfDNA levels than all other groups. When stratified the supervised OPLS-DA analysis and total cfDNA concentration showed high discrimination accuracy between invasive cancers and the disease/medication-free subjects. Furthermore, comparison of OPLS-DA data for invasive breast cancers with the AIRWAVE cohort showed similar discrimination between breast cancers and healthy controls. This is the first report of agreement between metabolomics and plasma cfDNA levels for discriminating breast cancer from healthy subjects in a true screening population. It also emphasizes the importance of sample standardization. Follow on studies will involve analysis of candidate features in a larger validation series as well as comparing results with serial plasma samples taken at the next routine screening mammography appointment. The findings here help establish the role of plasma analysis in the diagnosis of breast cancer in a large real-world cohort.

Data-driven approach for metabolite relationship recovery in biological 1H NMR data sets using iterative statistical total correlation spectroscopy.

Statistical total correlation spectroscopy (STOCSY) is a well-established and valuable method in the elucidation of both inter- and intrametabolite correlations in NMR metabonomic data sets. Here, the STOCSY approach is extended in a novel Iterative-STOCSY (I-STOCSY) tool in which correlations are calculated initially from a driver peak of interest and subsequently for all peaks identified as correlating with a correlation coefficient greater than a set threshold. Consequently, in a single automated run, the majority of information contained in multiple STOCSY calculations from all peaks recursively correlated to the original user defined driver peak of interest are recovered. In addition, highly correlating peaks are clustered into putative structurally related sets, and the results are presented in a fully interactive plot where each set is represented by a node; node-to-node connections are plotted alongside corresponding spectral data colored by the strength of connection, thus allowing the intuitive exploration of both inter- and intrametabolite connections. The I-STOCSY approach has been here applied to a (1)H NMR data set of 24 h postdose aqueous liver extracts from rats treated with the model hepatotoxin galactosamine and has been shown both to recover the previously deduced major metabolic effects of treatment and to generate new hypotheses even on this well-studied model system. I-STOCSY, thus, represents a significant advance in correlation based analysis and visualization, providing insight into inter- and intrametabolite relationships following metabolic perturbations.

Impact of Pelvic Radiation Therapy for Prostate Cancer on Global Metabolic Profiles and Microbiota-Driven Gastrointestinal Late Side Effects: A Longitudinal Observational Study.

PURPOSE: Radiation therapy to the prostate and pelvic lymph nodes (PLNRT) is part of the curative treatment of high-risk prostate cancer. Yet, the broader influence of radiation therapy on patient physiology is poorly understood. We conducted comprehensive global metabolomic profiling of urine, plasma, and stools sampled from patients undergoing PLNRT for high-risk prostate cancer. METHODS AND MATERIALS: Samples were taken from 32 patients at 6 timepoints: baseline, 2 to 3 and 4 to 5 weeks of PLNRT; and 3, 6, and 12 months after PLNRT. We characterized the global metabolome of urine and plasma using 1H nuclear magnetic resonance spectroscopy and ultraperformance liquid chromatography-mass spectrometry, and of stools with nuclear magnetic resonance. Linear mixed-effects modeling was used to investigate metabolic changes between timepoints for each biofluid and assay and determine metabolites of interest. RESULTS: Metabolites in urine, plasma and stools changed significantly after PLNRT initiation. Metabolic profiles did not return to baseline up to 1 year post-PLNRT in any biofluid. Molecules associated with cardiovascular risk were increased in plasma. Pre-PLNRT fecal butyrate levels directly associated with increasing gastrointestinal side effects, as did a sharper fall in those levels during and up to 1 year postradiation therapy, mirroring our previous results with metataxonomics. CONCLUSIONS: We showed for the first time that an overall metabolic effect is observed in patients undergoing PLNRT up to 1 year posttreatment. These metabolic changes may effect on long-term morbidity after treatment, which warrants further investigation.

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.

Statistical total correlation spectroscopy editing of 1H NMR spectra of biofluids: application to drug metabolite profile identification and enhanced information recovery.

Here we present a novel method for enhanced NMR spectral information recovery, utilizing a statistical total correlation spectroscopy editing (STOCSY-E) procedure for the identification of drug metabolite peaks in biofluids and for deconvolution of drug and endogenous metabolite signals. Structurally correlated peaks from drug metabolites and those from closely related drug metabolite pathways are first identified using STOCSY. Subsequently, this correlation information is utilized to scale the biofluid (1)H NMR spectra across these identified regions, producing a modified set of spectra in which drug metabolite contributions are reduced and, thus, facilitating analysis by pattern recognition methods without drug metabolite interferences. The application of STOCSY-E is illustrated with two exemplar (1)H NMR spectroscopic data sets, posing various drug metabolic, toxicological, and analytical challenges viz. 800 MHz (1)H spectra of human urine (n = 21) collected over 10 h following dosing with the antibiotic flucloxacillin and 600 MHz (1)H NMR spectra of rat urine (n = 27) collected over 48 h following exposure to the renal papillary toxin 2-bromoethanamine (BEA). STOCSY-E efficiently identified and removed the major xenobiotic metabolite peaks in both data sets, providing enhanced visualization of endogenous changes via orthogonal to projection filtered partial least-squares discriminant analysis (OPLS-DA). OPLS-DA of the STOCSY-E spectral data from the BEA-treated rats revealed the gut bacterial-mammalian co-metabolite phenylacetylglycine as a previously unidentified surrogate biomarker of toxicity. STOCSY-E has a wide range of potential applications in clinical, epidemiology, toxicology, and nutritional studies where multiple xenobiotic metabolic interferences may confound biological interpretation. Additionally, this tool could prove useful for applications outside of metabolic analysis, for example, in process chemistry for following chemical reactions and equilibria and detecting impurities.