Establishing a framework for best practices for quality assurance and quality control in untargeted metabolomics

BACKGROUND

Quality assurance (QA) and quality control (QC) practices are key tenets that facilitate study and data quality across all applications of untargeted metabolomics. These important practices will strengthen this field and accelerate its success. The Best Practices Working Group (WG) within the Metabolomics Quality Assurance and Quality Control Consortium (mQACC) focuses on community use of QA/QC practices and protocols and aims to identify, catalogue, harmonize, and disseminate current best practices in untargeted metabolomics through community-driven activities.

AIM OF REVIEW

A present goal of the Best Practices WG is to develop a working strategy, or roadmap, that guides the actions of practitioners and progress in the field. The framework in which mQACC operates promotes the harmonization and dissemination of current best QA/QC practice guidance and encourages widespread adoption of these essential QA/QC activities for liquid chromatography-mass spectrometry.

KEY SCIENTIFIC CONCEPTS OF REVIEW

Community engagement and QA/QC information gathering activities have been occurring through conference workshops, virtual and in-person interactive forum discussions, and community surveys. Seven principal QC stages prioritized by internal discussions of the Best Practices WG have received participant input, feedback and discussion. We outline these stages, each involving a multitude of activities, as the framework for identifying QA/QC best practices. The ultimate planned product of these endeavors is a "living guidance" document of current QA/QC best practices for untargeted metabolomics that will grow and change with the evolution of the field.

Small Intestinal Permeability and Metabolomic Profiles in Feces and Plasma Associate With Clinical Response in Patients With Active Psoriatic Arthritis Participating in a Fecal Microbiota Transplantation Trial: Exploratory Findings From the FLORA Trial

OBJECTIVE

We investigated intestinal permeability and fecal, plasma, and urine metabolomic profiles in methotrexate-treated active psoriatic arthritis (PsA) and how this related to clinical response following one sham or fecal microbiota transplantation (FMT).

METHODS

This exploratory study is based on the FLORA trial cohort, in which 31 patients with moderate-to-high peripheral PsA disease activity, despite at least 3 months of methotrexate-treatment, were included in a 26-week, double-blind, 1:1 randomized, sham-controlled trial. Participants were randomly allocated to receive either one healthy donor FMT (n = 15) or sham (n = 16) via gastroscopy. The primary trial end point was the proportion of treatment failures through 26 weeks. We performed a lactulose-to-mannitol ratio (LMR) test at baseline (n = 31) and at week 26 (n = 26) to assess small intestinal permeability. Metabolomic profiles in fecal, plasma, and urine samples collected at baseline, weeks 4, 12, and 26 were measured using 1 H Nuclear Magnetic Resonance.

RESULTS

Trial failures (n = 7) had significantly higher LMR compared with responders (n = 19) at week 26 (0.027 [0.017-0.33]) vs. 0.012 [0-0.064], P = 0.013), indicating increased intestinal permeability. Multivariate analysis revealed a significant model for responders (n = 19) versus failures (n = 12) at all time points based on their fecal (P < 0.0001) and plasma (P = 0.005) metabolomic profiles, whereas urine metabolomic profiles did not differ between groups (P = 1). Fecal N-acetyl glycoprotein GlycA correlated with Health Assessment Questionnaire Disability Index (coefficient = 0.50; P = 0.03) and fecal propionate correlated with American College of Rheumatology 20 response at week 26 (coefficient = 27, P = 0.02).

CONCLUSION

Intestinal permeability and fecal and plasma metabolomic profiles of patients with PsA were associated with the primary clinical trial end point, failure versus responder.

The postprandial secretion of peptide YY1-36 and 3-36 in obesity is differentially increased after gastric bypass versus sleeve gastrectomy

OBJECTIVES

Peptide tyrosine tyrosine (PYY) exists as two species, PYY1-36 and PYY3-36 , with distinct effects on insulin secretion and appetite regulation. The detailed effects of bariatric surgery on PYY1-36 and PYY3-36 secretion are not known as previous studies have used nonspecific immunoassays to measure total PYY. Our objective was to characterize the effect of sleeve gastrectomy (SG) and Roux-en-Y gastric bypass (RYGB) on fasting and postprandial PYY1-36 and PYY3-36 secretion using a newly developed liquid chromatography-tandem mass spectrometry (LC-MS/MS) assay.

DESIGN AND SUBJECTS

Observational study in 10 healthy nonobese volunteers and 30 participants with obesity who underwent RYGB (n = 24) or SG (n = 6) at the Imperial Weight Centre [NCT01945840]. Participants were studied using a standardized mixed meal test (MMT) before and 1 year after surgery. The outcome measures were PYY1-36 and PYY3-36 concentrations.

RESULTS

Presurgery, the fasting and postprandial levels of PYY1-36 and PYY3-36 were low, with minimal responses to the MMT, and these did not differ from healthy nonobese volunteers. The postprandial secretion of both PYY1-36 and PYY3-36 at 1 year was amplified after RYGB, but not SG, with the response being significantly higher in RYGB compared with SG.

CONCLUSIONS

There appears to be no difference in PYY secretion between nonobese and obese volunteers at baseline. At 1 year after surgery, RYGB, but not SG, is associated with increased postprandial secretion of PYY1-36 and PYY3-36 , which may account for long-term differences in efficacy and adverse effects between the two types of surgery.

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.

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.

A proposed framework to evaluate the quality and reliability of targeted metabolomics assays from the UK Consortium on Metabolic Phenotyping (MAP/UK)

Targeted metabolite assays that measure tens or hundreds of pre-selected metabolites, typically using liquid chromatography-mass spectrometry, are increasingly being developed and applied to metabolic phenotyping studies. These are used both as standalone phenotyping methods and for the validation of putative metabolic biomarkers obtained from untargeted metabolomics studies. However, there are no widely accepted standards in the scientific community for ensuring reliability of the development and validation of targeted metabolite assays (referred to here as 'targeted metabolomics'). Most current practices attempt to adopt, with modifications, the strict guidance provided by drug regulatory authorities for analytical methods designed largely for measuring drugs and other xenobiotic analytes. Here, the regulatory guidance provided by the European Medicines Agency, US Food and Drug Administration and International Council for Harmonisation of Technical Requirements for Pharmaceuticals for Human Use are summarized. In this Perspective, we have adapted these guidelines and propose a less onerous 'tiered' approach to evaluate the reliability of a wide range of metabolomics analyses, addressing the need for community-accepted, harmonized guidelines for tiers other than full validation. This 'fit-for-purpose' tiered approach comprises four levels-discovery, screening, qualification and validation-and is discussed in the context of a range of targeted and untargeted metabolomics assays. Issues arising with targeted multiplexed metabolomics assays, and how these might be addressed, are considered. Furthermore, guidance is provided to assist the community with selecting the appropriate degree of reliability for a series of well-defined applications of metabolomics.

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.

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.

Improved quantitation of short-chain carboxylic acids in human biofluids using 3-nitrophenylhydrazine derivatization and liquid chromatography with tandem mass spectrometry (LC-MS/MS)

Short-chain carboxylic acids (SCCAs) produced by gut microbial fermentation may reflect gastrointestinal health. Their concentrations in serum and urine are indicative of specific metabolic pathway activity; therefore, accurate quantitation of SCCAs in different biofluids is desirable. However, it is often challenging to quantitate SCCAs since matrix effects, induced by the presence of a vast variety of other compounds other than SCCAs in complex biofluids, can suppress or enhance signals. Materials used for sample preparation may introduce further analytical challenges. This study reports for the first time a LC-MS/MS-based method to quantitate ten SCCAs (lactate, acetate, 2-hydroxybutyrate, propionate, isobutyrate, butyrate, 2-methylbutyrate, isovalerate, valerate and hexanoate) and evaluates the matrix effects in five human biofluids: serum, urine, stool, and contents from the duodenum and intestinal stoma bags. The optimized method, using 3-Nitrophenylhydrazone as a derivatization agent and a Charge Surface Hybrid reverse phase column, showed clear separation for all SCCAs at a concentration range of 0.1-100 µM, in a 10.5 min run without carry-over effects. The validation of the method showed a good linearity (R2 > 0.99), repeatability (CV ≤ 15%) assessed by intra- and inter-day monitoring. The lowest limit of detection (LLOD) was 25 nM and lowest limit of quantitation (LLOQ) was 50 nM for nine SCCA except acetate at 0.5 and 1 µM, respectively. Quantitative accuracy in all biofluids for most compounds was < ±15%. In summary, this methodology has the advantages over other techniques for its simple and fast sample preparation and a high level of selectivity, repeatability and robustness for SCCA quantification. It also reduced interferences from the matrix or sample containers, making it ideal for use in high-throughput analyses of biofluid samples from large-scale studies.

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.

NMRpQuant: an automated software for large scale urinary total protein quantification by one-dimensional 1H NMR profiles

SUMMARY

1H nuclear magnetic resonance (NMR) spectroscopy is an established bioanalytical technology for metabolic profiling of biofluids in both clinical and large-scale population screening applications. Recently, urinary protein quantification has been demonstrated using the same 1D 1H NMR experimental data captured for metabolic profiling. Here, we introduce NMRpQuant, a freely available platform that builds on these findings with both novel and further optimized computational NMR approaches for rigorous, automated protein urine quantification. The results are validated by interlaboratory comparisons, demonstrating agreement with clinical/biochemical methodologies, pointing at a ready-to-use tool for routine protein urinalyses.

AVAILABILITY AND IMPLEMENTATION

NMRpQuant was developed on MATLAB programming environment. Source code and Windows/macOS compiled applications are available at https://github.com/pantakis/NMRpQuant, and working examples are available at https://doi.org/10.6084/m9.figshare.18737189.v1.

SUPPLEMENTARY INFORMATION

Supplementary data are available at Bioinformatics online.

Evidence of altered fatty acid metabolism in dogs with naturally occurring valvular heart disease and congestive heart failure

INTRODUCTION

Myxomatous mitral valve disease (MMVD) is the most common cardiac condition in adult dogs. The disease progresses over several years and affected dogs may develop congestive heart failure (HF). Research has shown that myocardial metabolism is altered in cardiac disease, leading to a reduction in β-oxidation of fatty acids and an increased dependence upon glycolysis.

OBJECTIVES

This study aimed to evaluate whether a shift in substrate use occurs in canine patients with MMVD; a naturally occurring model of human disease.

METHODS

Client-owned dogs were longitudinally evaluated at a research clinic in London, UK and paired serum samples were selected from visits when patients were in ACVIM stage B1: asymptomatic disease without cardiomegaly, and stage C: HF. Samples were processed using ultra-performance liquid chromatography mass spectrometry and lipid profiles were compared using mixed effects models with false discovery rate adjustment. The effect of disease stage was evaluated with patient breed entered as a confounder. Features that significantly differed were screened for selection for annotation efforts using reference databases.

RESULTS

Dogs in HF had altered concentrations of lipid species belonging to several classes previously associated with cardiovascular disease. Concentrations of certain acylcarnitines, phospholipids and sphingomyelins were increased after individuals had developed HF, whilst some ceramides and lysophosphatidylcholines decreased.

CONCLUSIONS

The canine metabolome appears to change as MMVD progresses. Findings from this study suggest that in HF myocardial metabolism may be characterised by reduced β-oxidation. This proposed explanation warrants further research.

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 ¹H nuclear magnetic resonance (NMR) spectroscopy is employed, the presence of urinary protein can elevate the spectral baseline and substantially impact the resulting profile. Using ¹H 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 (r² = 0.423, p < 2.2 × 10^-16). This correlation was significantly reduced (r² = 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 ¹H NMR metabolic profiling studies which can be effectively mitigated using SMolESY or other macromolecular signal suppression methods before estimation of normalization coefficients.

Mass recalibration for desorption electrospray ionization mass spectrometry imaging using endogenous reference ions

BACKGROUND

Mass spectrometry imaging (MSI) data often consist of tens of thousands of mass spectra collected from a sample surface. During the time necessary to perform a single acquisition, it is likely that uncontrollable factors alter the validity of the initial mass calibration of the instrument, resulting in mass errors of magnitude significantly larger than their theoretical values. This phenomenon has a two-fold detrimental effect: (a) it reduces the ability to interpret the results based on the observed signals, (b) it can affect the quality of the observed signal spatial distributions.

RESULTS

We present a post-acquisition computational method capable of reducing the observed mass drift by up to 60 ppm in biological samples, exploiting the presence of typical molecules with a known mass-to-charge ratio. The procedure, tested on time-of-flight and Orbitrap mass spectrometry analyzers interfaced to a desorption electrospray ionization (DESI) source, improves the molecular annotation quality and the spatial distributions of the detected ions.

CONCLUSION

The presented method represents a robust and accurate tool for performing post-acquisition mass recalibration of DESI-MSI datasets and can help to increase the reliability of the molecular assignment and the data quality.

Mass recalibration for desorption electrospray ionization mass spectrometry imaging using endogenous reference ions

BACKGROUND

Mass spectrometry imaging (MSI) data often consist of tens of thousands of mass spectra collected from a sample surface. During the time necessary to perform a single acquisition, it is likely that uncontrollable factors alter the validity of the initial mass calibration of the instrument, resulting in mass errors of magnitude significantly larger than their theoretical values. This phenomenon has a two-fold detrimental effect: (a) it reduces the ability to interpret the results based on the observed signals, (b) it can affect the quality of the observed signal spatial distributions.

RESULTS

We present a post-acquisition computational method capable of reducing the observed mass drift by up to 60 ppm in biological samples, exploiting the presence of typical molecules with a known mass-to-charge ratio. The procedure, tested on time-of-flight and Orbitrap mass spectrometry analyzers interfaced to a desorption electrospray ionization (DESI) source, improves the molecular annotation quality and the spatial distributions of the detected ions.

CONCLUSION

The presented method represents a robust and accurate tool for performing post-acquisition mass recalibration of DESI-MSI datasets and can help to increase the reliability of the molecular assignment and the data quality.

Reference materials for MS-based untargeted metabolomics and lipidomics: a review by the metabolomics quality assurance and quality control consortium (mQACC)

INTRODUCTION

The metabolomics quality assurance and quality control consortium (mQACC) is enabling the identification, development, prioritization, and promotion of suitable reference materials (RMs) to be used in quality assurance (QA) and quality control (QC) for untargeted metabolomics research.

OBJECTIVES

This review aims to highlight current RMs, and methodologies used within untargeted metabolomics and lipidomics communities to ensure standardization of results obtained from data analysis, interpretation and cross-study, and cross-laboratory comparisons. The essence of the aims is also applicable to other 'omics areas that generate high dimensional data.

RESULTS

The potential for game-changing biochemical discoveries through mass spectrometry-based (MS) untargeted metabolomics and lipidomics are predicated on the evolution of more confident qualitative (and eventually quantitative) results from research laboratories. RMs are thus critical QC tools to be able to assure standardization, comparability, repeatability and reproducibility for untargeted data analysis, interpretation, to compare data within and across studies and across multiple laboratories. Standard operating procedures (SOPs) that promote, describe and exemplify the use of RMs will also improve QC for the metabolomics and lipidomics communities.

CONCLUSIONS

The application of RMs described in this review may significantly improve data quality to support metabolomics and lipidomics research. The continued development and deployment of new RMs, together with interlaboratory studies and educational outreach and training, will further promote sound QA practices in the community.

Circulating Metabolome and White Matter Hyperintensities in Women and Men

BACKGROUND

White matter hyperintensities (WMH), identified on T2-weighted magnetic resonance images of the human brain as areas of enhanced brightness, are a major risk factor of stroke, dementia, and death. There are no large-scale studies testing associations between WMH and circulating metabolites.

METHODS

We studied up to 9290 individuals (50.7% female, average age 61 years) from 15 populations of 8 community-based cohorts. WMH volume was quantified from T2-weighted or fluid-attenuated inversion recovery images or as hypointensities on T1-weighted images. Circulating metabolomic measures were assessed with mass spectrometry and nuclear magnetic resonance spectroscopy. Associations between WMH and metabolomic measures were tested by fitting linear regression models in the pooled sample and in sex-stratified and statin treatment-stratified subsamples. Our basic models were adjusted for age, sex, age×sex, and technical covariates, and our fully adjusted models were also adjusted for statin treatment, hypertension, type 2 diabetes, smoking, body mass index, and estimated glomerular filtration rate. Population-specific results were meta-analyzed using the fixed-effect inverse variance-weighted method. Associations with false discovery rate (FDR)-adjusted P values (PFDR)<0.05 were considered significant.

RESULTS

In the meta-analysis of results from the basic models, we identified 30 metabolomic measures associated with WMH (PFDR<0.05), 7 of which remained significant in the fully adjusted models. The most significant association was with higher level of hydroxyphenylpyruvate in men (PFDR.full.adj=1.40×10-7) and in both the pooled sample (PFDR.full.adj=1.66×10-4) and statin-untreated (PFDR.full.adj=1.65×10-6) subsample. In men, hydroxyphenylpyruvate explained 3% to 14% of variance in WMH. In men and the pooled sample, WMH were also associated with lower levels of lysophosphatidylcholines and hydroxysphingomyelins and a larger diameter of low-density lipoprotein particles, likely arising from higher triglyceride to total lipids and lower cholesteryl ester to total lipids ratios within these particles. In women, the only significant association was with higher level of glucuronate (PFDR=0.047).

CONCLUSIONS

Circulating metabolomic measures, including multiple lipid measures (eg, lysophosphatidylcholines, hydroxysphingomyelins, low-density lipoprotein size and composition) and nonlipid metabolites (eg, hydroxyphenylpyruvate, glucuronate), associate with WMH in a general population of middle-aged and older adults. Some metabolomic measures show marked sex specificities and explain a sizable proportion of WMH variance.

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.

Finding Correspondence between Metabolomic Features in Untargeted Liquid Chromatography-Mass Spectrometry Metabolomics Datasets

Integration of multiple datasets can greatly enhance bioanalytical studies, for example, by increasing power to discover and validate biomarkers. In liquid chromatography–mass spectrometry (LC–MS) metabolomics, it is especially hard to combine untargeted datasets since the majority of metabolomic features are not annotated and thus cannot be matched by chemical identity. Typically, the information available for each feature is retention time (RT), mass-to-charge ratio (m/z), and feature intensity (FI). Pairs of features from the same metabolite in separate datasets can exhibit small but significant differences, making matching very challenging. Current methods to address this issue are too simple or rely on assumptions that cannot be met in all cases. We present a method to find feature correspondence between two similar LC–MS metabolomics experiments or batches using only the features’ RT, m/z, and FI. We demonstrate the method on both real and synthetic datasets, using six orthogonal validation strategies to gauge the matching quality. In our main example, 4953 features were uniquely matched, of which 585 (96.8%) of 604 manually annotated features were correct. In a second example, 2324 features could be uniquely matched, with 79 (90.8%) out of 87 annotated features correctly matched. Most of the missed annotated matches are between features that behave very differently from modeled inter-dataset shifts of RT, MZ, and FI. In a third example with simulated data with 4755 features per dataset, 99.6% of the matches were correct. Finally, the results of matching three other dataset pairs using our method are compared with a published alternative method, metabCombiner, showing the advantages of our approach. The method can be applied using M2S (Match 2 Sets), a free, open-source MATLAB toolbox, available at https://github.com/rjdossan/M2S.

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.

Automated Annotation of Untargeted All-Ion Fragmentation LC-MS Metabolomics Data with MetaboAnnotatoR

Untargeted metabolomics and lipidomics LC-MS experiments produce complex datasets, usually containing tens of thousands of features from thousands of metabolites whose annotation requires additional MS/MS experiments and expert knowledge. All-ion fragmentation (AIF) LC-MS/MS acquisition provides fragmentation data at no additional experimental time cost. However, analysis of such datasets requires reconstruction of parent-fragment relationships and annotation of the resulting pseudo-MS/MS spectra. Here, we propose a novel approach for automated annotation of isotopologues, adducts, and in-source fragments from AIF LC-MS datasets by combining correlation-based parent-fragment linking with molecular fragment matching. Our workflow focuses on a subset of features rather than trying to annotate the full dataset, saving time and simplifying the process. We demonstrate the workflow in three human serum datasets containing 599 features manually annotated by experts. Precision and recall values of 82-92% and 82-85%, respectively, were obtained for features found in the highest-rank scores (1-5). These results equal or outperform those obtained using MS-DIAL software, the current state of the art for AIF data annotation. Further validation for other biological matrices and different instrument types showed variable precision (60-89%) and recall (10-88%) particularly for datasets dominated by nonlipid metabolites. The workflow is freely available as an open-source R package, MetaboAnnotatoR, together with the fragment libraries from Github (https://github.com/gggraca/MetaboAnnotatoR).

The Metabolomic Effects of Tripeptide Gut Hormone Infusion Compared to Roux-en-Y Gastric Bypass and Caloric Restriction

CONTEXT

The gut-derived peptide hormones glucagon-like peptide-1 (GLP-1), oxyntomodulin (OXM), and peptide YY (PYY) are regulators of energy intake and glucose homeostasis and are thought to contribute to the glucose-lowering effects of bariatric surgery.

OBJECTIVE

To establish the metabolomic effects of a combined infusion of GLP-1, OXM, and PYY (tripeptide GOP) in comparison to a placebo infusion, Roux-en-Y gastric bypass (RYGB) surgery, and a very low-calorie diet (VLCD).

DESIGN AND SETTING

Subanalysis of a single-blind, randomized, placebo-controlled study of GOP infusion (ClinicalTrials.gov NCT01945840), including VLCD and RYGB comparator groups.

PATIENTS AND INTERVENTIONS

Twenty-five obese patients with type 2 diabetes or prediabetes were randomly allocated to receive a 4-week subcutaneous infusion of GOP (n = 14) or 0.9% saline control (n = 11). An additional 22 patients followed a VLCD, and 21 underwent RYGB surgery.

MAIN OUTCOME MEASURES

Plasma and urine samples collected at baseline and 4 weeks into each intervention were subjected to cross-platform metabolomic analysis, followed by unsupervised and supervised modeling approaches to identify similarities and differences between the effects of each intervention.

RESULTS

Aside from glucose, very few metabolites were affected by GOP, contrasting with major metabolomic changes seen with VLCD and RYGB.

CONCLUSIONS

Treatment with GOP provides a powerful glucose-lowering effect but does not replicate the broader metabolomic changes seen with VLCD and RYGB. The contribution of these metabolomic changes to the clinical benefits of RYGB remains to be elucidated.

Structure Elucidation and Mitigation of Endogenous Interferences in LC-MS-Based Metabolic Profiling of Urine

Liquid chromatography-mass spectrometry (LC-MS) is the main workhorse of metabolomics owing to its high degree of analytical sensitivity and specificity when measuring diverse chemistry in complex biological samples. LC-MS-based metabolic profiling of human urine, a biofluid of primary interest for clinical and biobank studies, is not widely considered to be compromised by the presence of endogenous interferences and is often accomplished using a simple "dilute-and-shoot" approach. Yet, it is our experience that broad obscuring signals are routinely observed in LC-MS metabolic profiles and represent interferences that lack consideration in the relevant metabolomics literature. In this work, we chromatographically isolated the interfering metabolites from human urine and unambiguously identified them via de novo structure elucidation as two separate proline-containing dipeptides: N,N,N-trimethyl-l-alanine-l-proline betaine (l,l-TMAP) and N,N-dimethyl-l-proline-l-proline betaine (l,l-DMPP), the latter reported here for the first time. Offline LC-MS/MS, magnetic resonance mass spectrometry (MRMS), and nuclear magnetic resonance (NMR) spectroscopy were essential components of this workflow for the full chemical and spectroscopic characterization of these metabolites and for establishing the coexistence of cis and trans isomers of both dipeptides in solution. Analysis of these definitive structures highlighted intramolecular ionic interactions as responsible for slow interconversion between these isomeric forms resulting in their unusually broad elution profiles. Proposed mitigation strategies, aimed at increasing the quality of LC-MS-based urine metabolomics data, include modification of column temperature and mobile-phase pH to reduce the chromatographic footprint of these dipeptides, thereby reducing their interfering effect on the underlying metabolic profiles. Alternatively, sample dilution and internal standardization methods may be employed to reduce or account for the observed effects of ionization suppression on the metabolic profile.

The impact of bariatric surgery on serum tryptophan-kynurenine pathway metabolites

This study aims to explore the immediate effects of bariatric surgery on serum tryptophan–kynurenine pathway metabolites in individuals with type 2 diabetes and BMI > 30. With the goal of providing insight into the link between tryptophan pathway metabolites, type 2 diabetes, and chronic obesity-induced inflammation. This longitudinal study included 20 participants. Half were diagnosed with type 2 diabetes. 11 and 9 underwent RYGB and SG respectively. Blood samples were obtained at pre-operative and 3 months post-operative timepoints. Tryptophan and downstream metabolites of the kynurenine pathway were quantified with an ultrahigh-performance liquid chromatography tandem mass spectrometry with electrospray ionisation method. At 3 months post-operation, RYGB led to significant reductions in tryptophan, kynurenic acid and xanthurenic acid levels when compared to baseline. Significant reductions of the same metabolites after surgery were also observed in individuals with T2D irrespective of surgical procedure. These metabolites were significantly correlated with serum HbA1c levels and BMI. Bariatric surgery, in particular RYGB reduces serum levels of tryptophan and its downstream kynurenine metabolites. These metabolites are associated with T2D and thought to be potentially mechanistic in the systemic processes of obesity induced inflammation leading to insulin resistance. Its reduction after surgery is associated with an improvement in glycaemic control (HbA1c).

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.

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

 

Untargeted 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

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.

Roux-en-Y gastric bypass-induced bacterial perturbation contributes to altered host-bacterial co-metabolic phenotype

BACKGROUND

Bariatric surgery, used to achieve effective weight loss in individuals with severe obesity, modifies the gut microbiota and systemic metabolism in both humans and animal models. The aim of the current study was to understand better the metabolic functions of the altered gut microbiome by conducting deep phenotyping of bariatric surgery patients and bacterial culturing to investigate causality of the metabolic observations.

METHODS

Three bariatric cohorts (n = 84, n = 14 and n = 9) with patients who had undergone Roux-en-Y gastric bypass (RYGB), sleeve gastrectomy (SG) or laparoscopic gastric banding (LGB), respectively, were enrolled. Metabolic and 16S rRNA bacterial profiles were compared between pre- and post-surgery. Faeces from RYGB patients and bacterial isolates were cultured to experimentally associate the observed metabolic changes in biofluids with the altered gut microbiome.

RESULTS

Compared to SG and LGB, RYGB induced the greatest weight loss and most profound metabolic and bacterial changes. RYGB patients showed increased aromatic amino acids-based host-bacterial co-metabolism, resulting in increased urinary excretion of 4-hydroxyphenylacetate, phenylacetylglutamine, 4-cresyl sulphate and indoxyl sulphate, and increased faecal excretion of tyramine and phenylacetate. Bacterial degradation of choline was increased as evidenced by altered urinary trimethylamine-N-oxide and dimethylamine excretion and faecal concentrations of dimethylamine. RYGB patients' bacteria had a greater capacity to produce tyramine from tyrosine, phenylalanine to phenylacetate and tryptophan to indole and tryptamine, compared to the microbiota from non-surgery, normal weight individuals. 3-Hydroxydicarboxylic acid metabolism and urinary excretion of primary bile acids, serum BCAAs and dimethyl sulfone were also perturbed following bariatric surgery.

CONCLUSION

Altered bacterial composition and metabolism contribute to metabolic observations in biofluids of patients following RYGB surgery. The impact of these changes on the functional clinical outcomes requires further investigation. Video abstract.

A Computationally Lightweight Algorithm for Deriving Reliable Metabolite Panel Measurements from 1D 1H NMR

Small Molecule Enhancement SpectroscopY (SMolESY) was employed to develop a unique and fully automated computational solution for the assignment and integration of 1H nuclear magnetic resonance (NMR) signals from metabolites in challenging matrices containing macromolecules (herein blood products). Sensitive and reliable quantitation is provided by instant signal deconvolution and straightforward integration bolstered by spectral resolution enhancement and macromolecular signal suppression. The approach is highly efficient, requiring only standard one-dimensional 1H NMR spectra and avoiding the need for sample preprocessing, complex deconvolution, and spectral baseline fitting. The performance of the algorithm, developed using >4000 NMR serum and plasma spectra, was evaluated using an additional >8800 spectra, yielding an assignment accuracy greater than 99.5% for all 22 metabolites targeted. Further validation of its quantitation capabilities illustrated a reliable performance among challenging phenotypes. The simplicity and complete automation of the approach support the application of NMR-based metabolite panel measurements in clinical and population screening applications.

Unaltered Liver Regeneration in Post-Cholestatic Rats Treated with the FXR Agonist Obeticholic Acid

In a previous study, obeticholic acid (OCA) increased liver growth before partial hepatectomy (PHx) in rats through the bile acid receptor farnesoid X-receptor (FXR). In that model, OCA was administered during obstructive cholestasis. However, patients normally undergo PHx several days after biliary drainage. The effects of OCA on liver regeneration were therefore studied in post-cholestatic Wistar rats. Rats underwent sham surgery or reversible bile duct ligation (rBDL), which was relieved after 7 days. PHx was performed one day after restoration of bile flow. Rats received 10 mg/kg OCA per day or were fed vehicle from restoration of bile flow until sacrifice 5 days after PHx. Liver regeneration was comparable between cholestatic and non-cholestatic livers in PHx-subjected rats, which paralleled liver regeneration a human validation cohort. OCA treatment induced ileal Fgf15 mRNA expression but did not enhance post-PHx hepatocyte proliferation through FXR/SHP signaling. OCA treatment neither increased mitosis rates nor recovery of liver weight after PHx but accelerated liver regrowth in rats that had not been subjected to rBDL. OCA did not increase biliary injury. Conclusively, OCA does not induce liver regeneration in post-cholestatic rats and does not exacerbate biliary damage that results from cholestasis. This study challenges the previously reported beneficial effects of OCA in liver regeneration in cholestatic rats.

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.

Metabolic phenotyping reveals a reduction in the bioavailability of serotonin and kynurenine pathway metabolites in both the urine and serum of individuals living with Alzheimer's disease

BACKGROUND

Both serotonergic signalling disruption and systemic inflammation have been associated with the pathogenesis of Alzheimer's disease (AD). The common denominator linking the two is the catabolism of the essential amino acid, tryptophan. Metabolism via tryptophan hydroxylase results in serotonin synthesis, whilst metabolism via indoleamine 2,3-dioxygenase (IDO) results in kynurenine and its downstream derivatives. IDO is reported to be activated in times of host systemic inflammation and therefore is thought to influence both pathways. To investigate metabolic alterations in AD, a large-scale metabolic phenotyping study was conducted on both urine and serum samples collected from a multi-centre clinical cohort, consisting of individuals clinically diagnosed with AD, mild cognitive impairment (MCI) and age-matched controls.

METHODS

Metabolic phenotyping was applied to both urine (n = 560) and serum (n = 354) from the European-wide AddNeuroMed/Dementia Case Register (DCR) biobank repositories. Metabolite data were subsequently interrogated for inter-group differences; influence of gender and age; comparisons between two subgroups of MCI - versus those who remained cognitively stable at follow-up visits (sMCI); and those who underwent further cognitive decline (cMCI); and the impact of selective serotonin reuptake inhibitor (SSRI) medication on metabolite concentrations.

RESULTS

Results revealed significantly lower metabolite concentrations of tryptophan pathway metabolites in the AD group: serotonin (urine, serum), 5-hydroxyindoleacetic acid (urine), kynurenine (serum), kynurenic acid (urine), tryptophan (urine, serum), xanthurenic acid (urine, serum), and kynurenine/tryptophan ratio (urine). For each listed metabolite, a decreasing trend in concentrations was observed in-line with clinical diagnosis: control > MCI > AD. There were no significant differences in the two MCI subgroups whilst SSRI medication status influenced observations in serum, but not urine.

CONCLUSIONS

Urine and serum serotonin concentrations were found to be significantly lower in AD compared with controls, suggesting the bioavailability of the neurotransmitter may be altered in the disease. A significant increase in the kynurenine/tryptophan ratio suggests that this may be a result of a shift to the kynurenine metabolic route due to increased IDO activity, potentially as a result of systemic inflammation. Modulation of the pathways could help improve serotonin bioavailability and signalling in AD patients.

A targeted ultra performance liquid chromatography - Tandem mass spectrometric assay for tyrosine and metabolites in urine and plasma: Application to the effects of antibiotics on mice

Tyrosine plays a key role in mammalian biochemistry and defects in its metabolism (e.g., tyrosinemia, alkaptonuria etc.) have significant adverse consequences for those affected if left untreated. In addition, gut bacterially-derived p-cresol and its metabolites are of interest as a result of various effects on host xenobiotic metabolism. A fit-for-purpose quantitative ultra-performance liquid chromatography tandem mass spectrometry (UPLC-MS/MS) assay was developed to target and quantify tyrosine and eleven metabolites in urine and plasma. Dansylation, using dansyl chloride, was used to improve chromatographic and mass spectral properties for tyrosine and nine phenolic metabolites, with detection using positive electrospray ionisation (ESI). The sulfate and glucuronide conjugates of p-cresol, where the phenol group was blocked, were quantified intact, using negative ESI via polarity switching during the same run. Sample preparation for urine and plasma involved deproteinization by solvent precipitation (of acetonitrile:isopropyl alcohol (1:1 v/v)) followed by in situ dansylation in 96 well plates. To minimize sample and solvent usage, and maximize sensitivity, analysis was performed using microbore reversed-phase gradient UPLC on a C8 phase with a 7.5 min. cycle time. The coefficients of variation obtained were <15%, with lower limits of quantification ranging from 5 to 250 nM depending upon the analyte. The method was applied to plasma and urine samples obtained from mice placed on a high tyrosine diet with one subgroup of animals subsequently receiving antibiotics to suppress the gut microbiota. Whilst plasma profiles were largely unaffected by antibiotic treatment clear reductions in the amount of p-cresol sulfate and p-cresol glucuronide excreted in the urine were observed for these mice.

Dissemination and analysis of the quality assurance (QA) and quality control (QC) practices of LC-MS based untargeted metabolomics practitioners

INTRODUCTION

The metabolomics quality assurance and quality control consortium (mQACC) evolved from the recognized need for a community-wide consensus on improving and systematizing quality assurance (QA) and quality control (QC) practices for untargeted metabolomics.

OBJECTIVES

In this work, we sought to identify and share the common and divergent QA and QC practices amongst mQACC members and collaborators who use liquid chromatography-mass spectrometry (LC-MS) in untargeted metabolomics.

METHODS

All authors voluntarily participated in this collaborative research project by providing the details of and insights into the QA and QC practices used in their laboratories. This sharing was enabled via a six-page questionnaire composed of over 120 questions and comment fields which was developed as part of this work and has proved the basis for ongoing mQACC outreach.

RESULTS

For QA, many laboratories reported documenting maintenance, calibration and tuning (82%); having established data storage and archival processes (71%); depositing data in public repositories (55%); having standard operating procedures (SOPs) in place for all laboratory processes (68%) and training staff on laboratory processes (55%). For QC, universal practices included using system suitability procedures (100%) and using a robust system of identification (Metabolomics Standards Initiative level 1 identification standards) for at least some of the detected compounds. Most laboratories used QC samples (>86%); used internal standards (91%); used a designated analytical acquisition template with randomized experimental samples (91%); and manually reviewed peak integration following data acquisition (86%). A minority of laboratories included technical replicates of experimental samples in their workflows (36%).

CONCLUSIONS

Although the 23 contributors were researchers with diverse and international backgrounds from academia, industry and government, they are not necessarily representative of the worldwide pool of practitioners due to the recruitment method for participants and its voluntary nature. However, both questionnaire and the findings presented here have already informed and led other data gathering efforts by mQACC at conferences and other outreach activities and will continue to evolve in order to guide discussions for recommendations of best practices within the community and to establish internationally agreed upon reporting standards. We very much welcome further feedback from readers of this article.

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.

Performance of metabonomic serum analysis for diagnostics in paediatric tuberculosis

We applied a metabonomic strategy to identify host biomarkers in serum to diagnose paediatric tuberculosis (TB) disease. 112 symptomatic children with presumptive TB were recruited in The Gambia and classified as bacteriologically-confirmed TB, clinically diagnosed TB, or other diseases. Sera were analysed using ¹H nuclear magnetic resonance (NMR) spectroscopy and mass spectrometry (MS). Multivariate data analysis was used to distinguish patients with TB from other diseases. Diagnostic accuracy was evaluated using Receiver Operating Characteristic (ROC) curves. Model performance was tested in a validation cohort of 36 children from the UK. Data acquired using ¹H NMR demonstrated a sensitivity, specificity and Area Under the Curve (AUC) of 69% (95% confidence interval [CI], 56-73%), 83% (95% CI, 73-93%), and 0.78 respectively, and correctly classified 20% of the validation cohort from the UK. The most discriminatory MS data showed a sensitivity of 67% (95% CI, 60-71%), specificity of 86% (95% CI, 75-93%) and an AUC of 0.78, correctly classifying 83% of the validation cohort. Amongst children with presumptive TB, metabolic profiling of sera distinguished bacteriologically-confirmed and clinical TB from other diseases. This novel approach yielded a diagnostic performance for paediatric TB comparable to that of Xpert MTB/RIF and interferon gamma release assays.

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.

Plasma lipid profiles discriminate bacterial from viral infection in febrile children

Fever is the most common reason that children present to Emergency Departments. Clinical signs and symptoms suggestive of bacterial infection are often non-specific, and there is no definitive test for the accurate diagnosis of infection. The 'omics' approaches to identifying biomarkers from the host-response to bacterial infection are promising. In this study, lipidomic analysis was carried out with plasma samples obtained from febrile children with confirmed bacterial infection (n = 20) and confirmed viral infection (n = 20). We show for the first time that bacterial and viral infection produces distinct profile in the host lipidome. Some species of glycerophosphoinositol, sphingomyelin, lysophosphatidylcholine and cholesterol sulfate were higher in the confirmed virus infected group, while some species of fatty acids, glycerophosphocholine, glycerophosphoserine, lactosylceramide and bilirubin were lower in the confirmed virus infected group when compared with confirmed bacterial infected group. A combination of three lipids achieved an area under the receiver operating characteristic (ROC) curve of 0.911 (95% CI 0.81 to 0.98). This pilot study demonstrates the potential of metabolic biomarkers to assist clinicians in distinguishing bacterial from viral infection in febrile children, to facilitate effective clinical management and to the limit inappropriate use of antibiotics.

Systematic Isolation and Structure Elucidation of Urinary Metabolites Optimized for the Analytical-Scale Molecular Profiling Laboratory

Annotation and identification of metabolite biomarkers is critical for their biological interpretation in metabolic phenotyping studies, presenting a significant bottleneck in the successful implementation of untargeted metabolomics. Here, a systematic multistep protocol was developed for the purification and de novo structural elucidation of urinary metabolites. The protocol is most suited for instances where structure elucidation and metabolite annotation are critical for the downstream biological interpretation of metabolic phenotyping studies. First, a bulk urine pool was desalted using ion-exchange resins enabling large-scale fractionation using precise iterations of analytical scale chromatography. Primary urine fractions were collected and assembled into a “fraction bank” suitable for long-term laboratory storage. Secondary and tertiary fractionations exploited differences in selectivity across a range of reversed-phase chemistries, achieving the purification of metabolites of interest yielding an amount of material suitable for chemical characterization. To exemplify the application of the systematic workflow in a diverse set of cases, four metabolites with a range of physicochemical properties were selected and purified from urine and subjected to chemical formula and structure elucidation by respective magnetic resonance mass spectrometry (MRMS) and NMR analyses. Their structures were fully assigned as tetrahydropentoxyline, indole-3-acetic-acid-O-glucuronide, p-cresol glucuronide, and pregnanediol-3-glucuronide. Unused effluent was collected, dried, and returned to the fraction bank, demonstrating the viability of the system for repeat use in metabolite annotation with a high degree of efficiency.

Breaking the cycle: Reversal of flux in the tricarboxylic acid cycle by dimethyl fumarate

Objective

To infer molecular effectors of therapeutic effects and adverse events for dimethyl fumarate (DMF) in patients with relapsing-remitting MS (RRMS) using untargeted plasma metabolomics.

Methods

Plasma from 27 patients with RRMS was collected at baseline and 6 weeks after initiating DMF. Patients were separated into discovery (n = 15) and validation cohorts (n = 12). Ten healthy controls were also recruited. Metabolomic profiling using ultra-high-performance liquid chromatography mass spectrometry (UPLC-MS) was performed on the discovery cohort and healthy controls at Metabolon Inc (Durham, NC). UPLC-MS was performed on the validation cohort at the National Phenome Centre (London, UK). Plasma neurofilament concentration (pNfL) was assayed using the Simoa platform (Quanterix, Lexington, MA). Time course and cross-sectional analyses were performed to identify pharmacodynamic changes in the metabolome secondary to DMF and relate these to adverse events.

Results

In the discovery cohort, tricarboxylic acid (TCA) cycle intermediates fumarate and succinate, and TCA cycle metabolites succinyl-carnitine and methyl succinyl-carnitine increased 6 weeks following treatment (q < 0.05). Methyl succinyl-carnitine increased in the validation cohort (q < 0.05). These changes were not observed in the control population. Increased succinyl-carnitine and methyl succinyl-carnitine were associated with adverse events from DMF (flushing and abdominal symptoms). pNfL concentration was higher in patients with RRMS than in controls and reduced over 15 months of treatment.

Conclusion

TCA cycle intermediates and metabolites are increased in patients with RRMS treated with DMF. The results suggest reversal of flux through the succinate dehydrogenase complex. The contribution of succinyl-carnitine ester agonism at hydroxycarboxylic acid receptor 2 to both therapeutic effects and adverse events requires investigation.

Ultrahigh-Performance Liquid Chromatography Tandem Mass Spectrometry with Electrospray Ionization Quantification of Tryptophan Metabolites and Markers of Gut Health in Serum and Plasma-Application to Clinical and Epidemiology Cohorts

A targeted ultrahigh-performance liquid chromatography tandem mass spectrometry with electrospray ionization (UHPLC-ESI-MS/MS) method has been developed for the quantification of tryptophan and its downstream metabolites from the kynurenine and serotonin pathways. The assay coverage also includes markers of gut health and inflammation, including citrulline and neopterin. The method was designed in 96-well plate format for application in multiday, multiplate clinical and epidemiology population studies. A chromatographic cycle time of 7 min enables the analysis of two 96-well plates in 24 h. To protect chromatographic column lifespan, samples underwent a two-step extraction, using solvent protein precipitation followed by delipidation via solid-phase extraction (SPE). Analytical validation reported accuracy of each analyte <20% for the lowest limit of quantification and <15% for all other quality control (QC) levels. The analytical precision for each analyte was 2.1–12.9%. To test the applicability of the method to multiplate and multiday preparations, a serum pool underwent periodic repeat analysis during a run consisting of 18 plates. The % CV (coefficient of variation) values obtained for each analyte were <15%. Additional biological testing applied the assay to samples collected from healthy control participants and two groups diagnosed with inflammatory bowel disease (IBD) (one group treated with the anti-inflammatory 5-aminosalicylic acid (5-ASA) and one group untreated), with results showing significant differences in the concentrations of picolinic acid, kynurenine, and xanthurenic acid. The short analysis time and 96-well plate format of the assay makes it suitable for high-throughput targeted UHPLC-ESI-MS/MS metabolomic analysis in large-scale clinical and epidemiological population studies.

The effects of kisspeptin on β-cell function, serum metabolites and appetite in humans

AIMS

To investigate the effect of kisspeptin on glucose-stimulated insulin secretion and appetite in humans.

MATERIALS AND METHODS

In 15 healthy men (age: 25.2 ± 1.1 years; BMI: 22.3 ± 0.5 kg m-2 ), we compared the effects of 1 nmol kg-1 h-1 kisspeptin versus vehicle administration on glucose-stimulated insulin secretion, metabolites, gut hormones, appetite and food intake. In addition, we assessed the effect of kisspeptin on glucose-stimulated insulin secretion in vitro in human pancreatic islets and a human β-cell line (EndoC-βH1 cells).

RESULTS

Kisspeptin administration to healthy men enhanced insulin secretion following an intravenous glucose load, and modulated serum metabolites. In keeping with this, kisspeptin increased glucose-stimulated insulin secretion from human islets and a human pancreatic cell line in vitro. In addition, kisspeptin administration did not alter gut hormones, appetite or food intake in healthy men.

CONCLUSIONS

Collectively, these data demonstrate for the first time a beneficial role for kisspeptin in insulin secretion in humans in vivo. This has important implications for our understanding of the links between reproduction and metabolism in humans, as well as for the ongoing translational development of kisspeptin-based therapies for reproductive and potentially metabolic conditions.

Quantitative Lipoprotein Subclass and Low Molecular Weight Metabolite Analysis in Human Serum and Plasma by 1H NMR Spectroscopy in a Multilaboratory Trial

We report an extensive 600 MHz NMR trial of quantitative lipoprotein and small-molecule measurements in human blood serum and plasma. Five centers with eleven 600 MHz NMR spectrometers were used to analyze 98 samples including 20 quality controls (QCs), 37 commercially sourced, paired serum and plasma samples, and two National Institute of Science and Technology (NIST) reference material 1951c replicates. Samples were analyzed using rigorous protocols for sample preparation and experimental acquisition. A commercial lipoprotein subclass analysis was used to quantify 105 lipoprotein subclasses and 24 low molecular weight metabolites from the NMR spectra. For all spectrometers, the instrument specific variance in measuring internal QCs was lower than the percentage described by the National Cholesterol Education Program (NCEP) criteria for lipid testing [triglycerides <2.7%; cholesterol <2.8%; low-density lipoprotein (LDL) cholesterol <2.8%; high-density lipoprotein (HDL) cholesterol <2.3%], showing exceptional reproducibility for direct quantitation of lipoproteins in both matrixes. The average relative standard deviations (RSDs) for the 105 lipoprotein parameters in the 11 instruments were 4.6% and 3.9% for the two NIST samples, whereas they were 38% and 40% for the 37 commercially sourced plasmas and sera, respectively, showing negligible analytical compared to biological variation. The coefficient of variance (CV) obtained for the quantification of the small molecules across the 11 spectrometers was below 15% for 20 out of the 24 metabolites analyzed. This study provides further evidence of the suitability of NMR for high-throughput lipoprotein subcomponent analysis and small-molecule quantitation with the exceptional required reproducibility for clinical and other regulatory settings.

Guidelines and considerations for the use of system suitability and quality control samples in mass spectrometry assays applied in untargeted clinical metabolomic studies

BACKGROUND

Quality assurance (QA) and quality control (QC) are two quality management processes that are integral to the success of metabolomics including their application for the acquisition of high quality data in any high-throughput analytical chemistry laboratory. QA defines all the planned and systematic activities implemented before samples are collected, to provide confidence that a subsequent analytical process will fulfil predetermined requirements for quality. QC can be defined as the operational techniques and activities used to measure and report these quality requirements after data acquisition.

AIM OF REVIEW

This tutorial review will guide the reader through the use of system suitability and QC samples, why these samples should be applied and how the quality of data can be reported.

KEY SCIENTIFIC CONCEPTS OF REVIEW

System suitability samples are applied to assess the operation and lack of contamination of the analytical platform prior to sample analysis. Isotopically-labelled internal standards are applied to assess system stability for each sample analysed. Pooled QC samples are applied to condition the analytical platform, perform intra-study reproducibility measurements (QC) and to correct mathematically for systematic errors. Standard reference materials and long-term reference QC samples are applied for inter-study and inter-laboratory assessment of data.

The pathophysiology of human obstructive cholestasis is mimicked in cholestatic Gold Syrian hamsters

Obstructive cholestasis causes liver injury via accumulation of toxic bile acids (BAs). Therapeutic options for cholestatic liver disease are limited, partially because the available murine disease models lack translational value. Profiling of time-related changes following bile duct ligation (BDL) in Gold Syrian hamsters revealed a biochemical response similar to cholestatic patients in terms of BA pool composition, alterations in hepatocyte BA transport and signaling, suppression of BA production, and adapted BA metabolism. Hamsters tolerated cholestasis well for up to 28days and progressed relatively slowly to fibrotic liver injury. Hepatocellular necrosis was absent, which coincided with preserved intrahepatic energy levels and only mild oxidative stress. The histological response to cholestasis in hamsters was similar to the changes seen in 17 patients with prolonged obstructive cholestasis caused by cholangiocarcinoma. Hamsters moreover upregulated hepatic fibroblast growth factor 15 (Fgf15) expression in response to BDL, which is a cytoprotective adaptation to cholestasis that hitherto had only been documented in cholestatic human livers. Hamster models should therefore be added to the repertoire of animal models used to study the pathophysiology of cholestatic liver disease.

Improving Visualization and Interpretation of Metabolome-Wide Association Studies: An Application in a Population-Based Cohort Using Untargeted 1H NMR Metabolic Profiling

¹H NMR spectroscopy of biofluids generates reproducible data allowing detection and quantification of small molecules in large population cohorts. Statistical models to analyze such data are now well-established, and the use of univariate metabolome wide association studies (MWAS) investigating the spectral features separately has emerged as a computationally efficient and interpretable alternative to multivariate models. The MWAS rely on the accurate estimation of a metabolome wide significance level (MWSL) to be applied to control the family wise error rate. Subsequent interpretation requires efficient visualization and formal feature annotation, which, in-turn, call for efficient prioritization of spectral variables of interest. Using human serum ¹H NMR spectroscopic profiles from 3948 participants from the Multi-Ethnic Study of Atherosclerosis (MESA), we have performed a series of MWAS for serum levels of glucose. We first propose an extension of the conventional MWSL that yields stable estimates of the MWSL across the different model parameterizations and distributional features of the outcome. We propose both efficient visualization methods and a strategy based on subsampling and internal validation to prioritize the associations. Our work proposes and illustrates practical and scalable solutions to facilitate the implementation of the MWAS approach and improve interpretation in large cohort studies.

Parkinson's disease and Parkinson's disease medications have distinct signatures of the gut microbiome

BACKGROUND

There is mounting evidence for a connection between the gut and Parkinson's disease (PD). Dysbiosis of gut microbiota could explain several features of PD.

OBJECTIVE

The objective of this study was to determine if PD involves dysbiosis of gut microbiome, disentangle effects of confounders, and identify candidate taxa and functional pathways to guide research.

METHODS

A total of 197 PD cases and 130 controls were studied. Microbial composition was determined by 16S rRNA gene sequencing of DNA extracted from stool. Metadata were collected on 39 potential confounders including medications, diet, gastrointestinal symptoms, and demographics. Statistical analyses were conducted while controlling for potential confounders and correcting for multiple testing. We tested differences in the overall microbial composition, taxa abundance, and functional pathways.

RESULTS

Independent microbial signatures were detected for PD (P = 4E-5), participants' region of residence within the United States (P = 3E-3), age (P = 0.03), sex (P = 1E-3), and dietary fruits/vegetables (P = 0.01). Among patients, independent signals were detected for catechol-O-methyltransferase-inhibitors (P = 4E-4), anticholinergics (P = 5E-3), and possibly carbidopa/levodopa (P = 0.05). We found significantly altered abundances of the Bifidobacteriaceae, Christensenellaceae, [Tissierellaceae], Lachnospiraceae, Lactobacillaceae, Pasteurellaceae, and Verrucomicrobiaceae families. Functional predictions revealed changes in numerous pathways, including the metabolism of plant-derived compounds and xenobiotics degradation.

CONCLUSION

PD is accompanied by dysbiosis of gut microbiome. Results coalesce divergent findings of prior studies, reveal altered abundance of several taxa, nominate functional pathways, and demonstrate independent effects of PD medications on the microbiome. The findings provide new leads and testable hypotheses on the pathophysiology and treatment of PD. © 2017 International Parkinson and Movement Disorder Society.