Abstract
Objectives Recent advances in high-resolution metabolomics platforms allow the simultaneous measurement of thousands of small molecules produced from metabolism (metabolites), providing a map of disease-related perturbations across interconnected pathways. We used high performance, untargeted metabolomics to identify metabolic disturbances and molecular events associated with Parkinson’s disease (PD) in two population-based studies.
Methods We performed a metabolome-wide association study (MWAS) on PD, using serum-based untargeted metabolomics data derived from high resolution liquid chromatography, mass spectrometry (LCMS). We used two independent, case-control populations for discovery and replication (n=642 PD patients, n=277 controls).
Results From the LCMS, 5,145 metabolites were detected across the two study populations in ≥50% of the samples (HILIC: 2913 metabolites; C18: 2063 metabolites). Using logistic regression and an FDR to correct for multiple testing, we determined 236 metabolites were associated with PD in a meta-analysis at an FDR<0.05. Of these, 110 metabolites were independently associated with PD in both discovery and replication studies at p<0.05 (187 at p<0.10), while 24 were associated with levodopa-equivalent dose among the PD patients. Intriguingly, the microbial related metabolite, p-cresol (meta-OR=1.29, 95% CI=1.13, 1.47, FDR=0.01), and it’s two metabolites, p-cresol sulfate and p-cresol glucuronide, were found at higher intensity among the PD patients relative to controls. P-cresol glucuronide was also associated with motor symptoms among patients. Pyroglutamic acid (meta OR=3.79, 95% CI=2.60, 5.54; FDR=5.30E-09), the anti-inflammatory metabolite itaconate (meta OR=0.47, 95% CI=0.36, 0.61; FDR=8.44E-06), and cysteine-S-sulfate (meta OR=1.56, 95% CI=1.32, 1.83; FDR=1.66E-05) were also among the most strongly associated metabolites. Seventeen pathways were also enriched, including several related to amino acid and lipid metabolism.
Conclusions Our results revealed PD-associated metabolites in two independent study populations, implicating individual metabolites including p-cresol and itaconate, as well as suggesting metabolic disturbances in amino acid and lipid metabolism and inflammatory processes.
Introduction
Parkinson’s disease (PD) is a complex, multi-factorial neurodegenerative disease with multi-system involvement. Pathologically, PD is defined by the loss of dopaminergic neurons in the substantia nigra and widespread intracytoplasmic aggregations of misfolded α-synuclein(1). Rare genetic mutations have been identified in early-onset, familial PD, but idiopathic PD’s complex etiopathogenesis remains unclear(2).
High-throughput technological developments over the past decade have paved the way for agnostic analysis of multiple omic measures, providing insights into disease etiology. Genome-wide association studies (GWAS), for instance, have highlighted the role of endolysosomal (vesicle trafficking, lysosomes, and autophagy) and immune pathways in PD(3). Still, biologic processes are dynamic and operate through complex interactions between gene expression, protein function, and metabolism. Investigating other principal omics, including the metabolome, may provide novel insights into biologic processes involved in PD.
Recent advances in high-resolution metabolomics platforms allow for the simultaneous measurement of thousands of metabolites. Metabolic networks shed light on the underlying biochemical activity of cells, tissues, and organs, enabling a multi-system level approach to the study of PD. Furthermore, metabolites reflect the convergence of genomic, epigenomic, transcriptomic, and proteomic action in tandem with the system’s response to environmental exposures, thus offering a readout of both physiologic and pathologic states of an individual(4). Metabolites circulating in the blood provide a wealth of information about biologic processes across different systems, including the central nervous system as metabolites can cross the blood-brain barrier(5). A growing body of work supports the use of metabolomics to provide novel information about the initiation and progression of PD(5,6). For instance, metabolites related to lipid metabolism, including glycerophospholipids and sphingolipids, mitochondrial function, and amino acids have been implicated in PD(7,8).
Here, we have performed a series of untargeted, agnostic analyses to explore metabolite signatures associated with PD. Using serum-based untargeted, high-performance liquid-chromatography mass spectrometry (LCMS), we performed a hypothesis-generating metabolome-wide association study (MWAS) with independent discovery and replication study populations. We identified individual metabolite features associated with PD, evaluated pathway enrichment, and assessed associations between PD-MWAS metabolites and symptom profiles among patients. Ultimately, identifying disrupted metabolic pathways in PD may improve our understanding of the molecular mechanisms underlying pathogenesis, paving the way for new preventative or therapeutic strategies.
Methods
Study Population
We used metabolomics data from 642 PD patients and 277 controls recruited as part of a community- based study of Parkinson’s disease (Parkinson’s Environment and Genes study, PEG). PEG is a population-based PD case-control study conducted in three Central California counties(9). Participants were recruited in two waves: PEG1, 2000-2007 and PEG2, 2011-2018. All those with serum for metabolomics were included (PEG1: n=282 PD patients, n=185 controls; PEG2: n=360 PD patients, n=90 controls). Patients were early in their disease course at enrollment (3.0 years [SD=2.6] on average from diagnosis) and all were seen by UCLA movement disorder specialists for in-person neurologic exams and confirmed as having idiopathic PD based on clinical characteristics(10). Characteristics of the PEG study subjects are shown in Supplemental Table 1. The patients were on average slightly older than the controls and a higher proportion were men, Hispanic, and never smokers.
High-Resolution Metabolomics (HRM)
HRM profiling was conducted according to established methods. Detailed methods are provided in previous publication(11). Briefly, serum samples were randomly sorted into batches of 40. Each sample was thoroughly mixed with ice-cold acetonitrile (2:1 acetonitrile to serum), placed on ice for 30 minutes, precipitated protein was removed by centrifugation, and the resulting supernatant was transferred to an autosampler vial containing a low volume insert. We analyzed all sample extracts in triplicate with a dual-column, dual-polarity approach, including hydrophilic interaction (HILIC) chromatography with positive electrospray ionization (ESI) and C18 chromatography with negative ESI, and used two types of two quality control samples (see supplement). Details on raw data extraction, feature alignment, and retention time adjustment are also provided in the supplement. Uniquely detected ions consisted of m/z, retention time (rt), and ion abundance. For analyses, we included metabolomic features with median coefficients of variation among technical replicates <75% and Pearson correlation >0.7 and features detected in >50% of all study samples, leaving 2046 C18 features and 2716 HILIC features for analysis. We log 2 transformed, quantile normalized, and batch corrected the data with ComBat after replacing zeroes with the lowest detected value. Data pre-processing, batch and technical artifact correction, and visualization are further described in the supplemental materials and Supplemental Figures 1-7.
Metabolome-Wide Association Analysis (MWAS)
To associate features with PD, we conducted unconditional logistic regression for each metabolite with PD as the outcome and age, gender, race/ethnicity, and a year of sample draw indicator included as covariates. We determined metabolite associations independently for the PEG1 and PEG2 study waves and then combined odds ratio (OR) estimates in a fixed effects meta-analysis, using a generic inverse-variance method for pooling(12). We used a false discovery rate (FDR) to correct for multiple testing. Our a priori criteria for determining the highest level of PD association were metabolites with a meta-FDR<0.05 and independent association in both discovery (PEG1) and replication (PEG2) populations at p<0.05. For metabolites which showed association with PD, using linear regression, we further tested for association with the following PD and PD symptom related phenotypes among PD patients only: levodopa equivalent daily dose (LEDD), Hoehn Yahr (HY) stage, and Unified Parkinson’s disease Rating Scale Part III (UPDRS-III) score.
Annotation and Pathway Analysis
We annotated features based on three levels. First, significant features were matched to a database of authenticated chemical standards previously characterized in the Emory laboratory, i.e., metabolites confirmed using MS/MS and authentic standards, providing the strongest level of annotation(13,14). The error tolerance was set to 5 ppm and 30s for m/z and retention time. Additional m/z feature mapping was done based on mummichog annotations and xMSannotator. mummichog is a computational algorithm which uses metabolic pathways and networks to predict functional activity from untargeted metabolite feature tables, including providing annotations of features based upon predicted ions and pathway associations(15). With xMSannotator, accurate mass m/z for adducts formed under positive/negative ESI mode was matched to HMDB, KEGG, and LipidMaps with a mass error threshold of 10 ppm(16). xMSannotator uses correlations of intensities and retention time and assigns confidence scores based on a multilevel scoring algorithm (0–3, a higher score representing higher-confidence result), ensuring annotation accuracy. Only results with scores ≥2 were considered for annotations.
For pathway enrichment analysis we used metapone, which uses a permutation-based weighted hypergeometric test with joint pathway analysis using positive and negative ion mode data to avoid double counting and account for multiple-matching uncertainty with a weighting factor(17). Metabolic pathways were compiled from KEGG, mummichog, and the small molecule pathway database (SMPDB).
Results
Our untargeted metabolome-wide association study included 4762 features (2046 C18 and 2716 HILIC). We detected 213 associated with PD at a meta-analysis FDR<0.05. Of these features, 104 met the criteria for independent replication (discovery and replication p<0.05). Figure 1A (C18) and 1B (HILIC) show the MWAS results. In total, from the C18 column, 116 features were associated at a meta-analysis FDR<0.05 (372 at p<0.05), 66 of which were independently associated in discovery and replication at p<0.05 (116 at p<0.10). From the HILIC column, 97 features were associated at a meta-analysis FDR<0.05 (441 at p<0.05), 38 of which were associated in discovery and replication at p<0.05 (54 at p<0.10). The full MWAS summary statistics are provided in Supplemental Table 2 (C18 metabolites) and Supplemental Table 3 (HILIC metabolites). Annotation based on three-layers (in-house database of metabolites, mummichog annotations, and xMSannotator high confidence matches) for all features with discovery and replication at p<0.1 is provided in Supplemental Table 4 (C18) and Supplemental Table 5 (HILIC). The full xMSannotator stage 5 annotation results are provided in Supplemental Table 6 (C18) and Supplemental Table 7 (HILIC).
Manhattan Plot of MWAS results from metabolomics LCMS data derived from both the (A) C18 negative ion channel and (B) HILIC positive ion channel. Horizontal lines are shown at FDR≤0.05. (C) Violin plots show the top three annotated metabolites from the MWAS, excluding PD medication-associated metabolites, and separated by independent study populations.
Table 1 shows the top MWAS metabolites that were associated with PD in discovery and replication. As expected, the leading PD-associated features from both columns were related to PD-medications, including medication metabolites, dopamine 3-O-sulfate and 3-O-methyldopa. In total, 11 C18 features and 16 HILIC features were associated with LEDD at an FDR<0.05 among the PD patients (Supplemental Tables 2 and 3). Additional metabolites strongly associated with PD included pyroglutamic acid (meta-OR= 3.79, 95% CI=2.6, 5.54, FDR=3.5e-9), itaconate (meta-OR= 0.47, 95% CI=0.36, 0.61, FDR=7.3e-6), and cysteine-S-sulfate (meta-OR= 1.56, 95% CI=1.32, 1.83, FDR=8.2e-6) from the C18 column (Figure 1C). In the HILIC column, a series of PD distinguishing features annotated to multiple phospholipids, including lysophosphatidylcholines (LysoPC) and lysophosphatidylethanolamines (LysoPE). A LysoPC(18:1), for example, was found at higher intensity among patients in both discovery and replication (meta-OR= 3.85, 95% CI=2.37, 6.27, FDR=1.8e-5). Features from both columns which annotated to inosine were inversely associated with PD (C18: meta-OR= 0.82, 95% CI=0.74, 0.92, FDR=1.4e-2; HLIIC: meta-OR= 0.91, 95% CI=0.87, 0.96, FDR=9.4e-3).
PD MWAS hits: Annotated features from the MWAS associated with Parkinson’s disease in both discovery and replication study populations.
P-cresol and two of its metabolites, p-cresol sulfate and p-cresol glucuronide, were also found at higher intensity among the PD patients relative to controls (p-cresol meta-OR=1.29, 95% CI=1.13, 1.47, FDR=0.01). The distributions of these metabolites by PD and across discovery and replication populations are shown in Figure 2A. The p-cresol metabolites were also correlated with age at PD diagnosis (Figure 2B) and p-cresol glucuronide was associated with a higher Hoehn Yahr (HY) stage among PD patients (beta=0.02, SE=0.007, FDR=9.5e-2; Figure 2C). In fact, of all PD-associated features (MWAS meta-p<0.05), seven were also related to HY stage among PD patients at an FDR<0.05 (115 metabolites associated at p<0.05), including, as expected, the PD medication metabolite, 3-O-methyldopa, providing confidence in the results (Supplemental Table 8). Six PD- associated features were also associated with UPDRS-III at an FDR<0.05 (100 metabolites at p<0.05; Supplemental Table 9). However, other than the PD medication metabolites and p-cresol glucuronide, the features associated with either HY stage or UPDRS-III up to FDR<0.10 could not be annotated at high confidence.
P-cresol and two p-cresol metabolites are associated with (A) Parkinson’s disease, (B) age at diagnosis among PD patients, and (C) Hoehn- Yahr Stage among PD patients.
Given the interdependent nature of metabolites, we further assessed correlation patterns between the PD-associated features. Figure 3 shows a Pearson-correlation based network of all FDR<0.05 MWAS features from both columns. Several highly correlated clusters of features are visible, including a PD medication related cluster, a phospholipids cluster, and a cluster of several features correlated with p-cresol.
|Correlations|≥0.2 shown. Features which were annotated are named to the right, while other features are shown as blank nodes.
Multiple pathways were also overrepresented in the MWAS features. Based on metabolic pathway analysis, 17 pathways were overrepresented at an FDR<0.05 (58 at p<0.05; Figure 4 and Supplemental Table 10). Glutamine and glutamate metabolism, gabaergic synapse, methionine and cysteine metabolism, glycine, serine, alanine and threonine metabolism, and leukotriene metabolism were among the most significantly overrepresented pathways. Several phospholipid pathways, including glycerophospholipid and glycosphingolipid metabolism and phospholipase d and sphingolipid signaling pathways, were also overrepresented among the PD-associated metabolites.
Based on pathway analysis of untargeted PD MWAS features using a permutation-based weighted hypergeometric test (R, Metapone: a Bioconductor package for joint pathway testing for untargeted metabolomics data). Pathways with p<0.05 are shown. Lfdr= the local FDR value for each enrichment.
Discussion
Using high resolution, untargeted serum metabolic profiling based on LC-MS, we identified 213 metabolite features associated with PD and 17 overrepresented metabolic pathways. Our agnostic metabolomics approach implicated amino acid metabolism and phospholipid pathways as important in PD, along with multiple individual metabolites with compelling links to neurodegeneration. Serum was collected from PD patients early in disease and compared to community controls. These metabolites and pathways therefore may reflect disturbances due to disease pathogenesis and progression as well as compensatory or reactive mechanisms caused by disease or treatment.
One of the strengths of our study is that we assessed metabolite profiles from a large number of PD patients recruited early in disease who were undergoing a range of treatment courses. Thus, we were able to assess the relationship between all metabolite features and levodopa equivalent daily dose to determine which features associated with PD were also associated with medication use. Predictably, PD patients differed from controls most strongly in terms of levodopa metabolites or metabolites involved in dopamine metabolism. Network analysis further showed the PD medication metabolites clustered together, as expected, but were not significantly correlated with other PD-related metabolites. Several other pathways and specific metabolites were implicated independent of levodopa medication use. Another noteworthy metabolite that attests to the validity of our analyses is inosine, which is a precursor to urate, with anti-inflammatory properties. We observed inosine at lower intensity among the patients relative to controls. This is in line with previous studies and the notion that lower uric acid level may be involved in faster PD progression, which have even resulted in inosine supplementation trials(18–21).
One of the most intriguing findings from this analysis was p-cresol and its two metabolites, p-cresol sulfate and p-cresol glucuronide. We observed higher intensities of the metabolites among PD patients relative to controls in both discovery and replication populations. Moreover, among patients, intensity of all three p- cresol metabolites was related to age at diagnosis, with older patients showing higher intensity and p-cresol glucuronide was related to higher motor symptom scores. P-cresol is an exogenous uremic toxin primarily produced by gut bacteria which express p-cresol synthesizing enzymes that are not produced by human cells. Interestingly, while this is the first study to implicate these metabolites in blood, a small study previously found higher levels of p-cresol in the CSF of PD patients(22). Multiple studies have also linked p-cresol with autism(23– 25) and altered brain dopamine metabolism in neurodevelopment(26). Furthermore, gut dysbiosis has been linked to both PD and autism, including among our own patients(27), with some research even indicating that misfolded a-synuclein retrogradely propagates from the enteric to the central nervous system(28). P-cresol and its’ metabolites therefore represent very compelling targets for future mechanistic research.
Several tricarboxylic acid (TCA) cycle metabolites were also implicated as relevant to PD in our MWAS. PD patients had higher relative abundance of oxoglutaric acid (e.g. alpha-ketoglutarate) and lower levels of itaconate. Aside from clear implications for energy metabolism, which is highly pertinent given the known mitochondrial dysfunction in PD, itaconate holds key roles in immunometabolism (e.g., changes of metabolic pathways within immune cells)(29,30). Itaconate is a mitochondrial metabolite, produced in high amounts by macrophages and monocytes by diverting aconitate away from the TCA cycle during inflammatory activation(30). The primary function appears to be anti-inflammatory, supported by human studies showing that low levels of plasma itaconate coincide with excessive inflammation(30). Inflammation and neuroinflammation are principal features of PD. Thus, it is quite interesting that in both of our patient populations we found lower relative levels of this important anti-inflammatory immunometabolite.
We further identified several amino acids as differentially abundant in PD. Glutamate and several connected metabolites, including pyroglutamic acid (PGA), had a higher relative abundance in patients’ serum relative to controls, with the glutamine and glutamate metabolism pathways significantly overrepresented. PGA is an endogenous metabolite derived from glutamate and linked to glutathione turnover(31,32). Elevated serum PGA therefore may be related to perturbed glutathione metabolism. Low levels of the antioxidant glutathione are an early neuronal biochemical finding in PD(33). But increased systemic levels of PGA may reflect an upregulation of glutathione metabolism to counter inflammatory states and oxidative stress in PD. Furthermore, the neurotransmitter glutamate itself has been linked to PD pathogenesis, with several, though not all, studies reporting increased blood-measured levels of glutamate(34,35). Other amino acid metabolic pathways, including methionine and cysteine metabolism, glycine, serine, alanine and threonine metabolism, and valine, leucine, and isoleucine degradation were also overrepresented among PD-associated features, with individual metabolites like serine, isoleucine, and tryptophan observed in higher relative abundance among the PD patients. The patients also had higher levels of cysteine-S-sulfate, a purportedly brain damaging metabolite involved in sulfite oxidase deficiency(36). Branched chain amino acids (BCAAs), including leucine, isoleucine, and valine, have been linked to PD before, as BCAAs are involved in energy metabolism, preventing oxidative damage, and regulation of protein synthesis(5).
Lipid pathways and metabolites were also strongly implicated with PD by our MWAS. Glycerophospholipid along with glycosphingolipid and sphingolipid metabolism were enriched in pathway analyses. Metabolites including glycerophosphocholine, several lysophosphotidylcholines (LysoPC), and a ceramide, were all observed at higher intensities among the patients relative to controls. The lipid profile in PD has received a great deal of interest in recent years, in part due the identification of GBA variants in PD GWAS. The glucosylceramidase-beta (GBA) gene, which encodes the lysosomal enzyme glucocerebrosidase (GCase), has directly connected lipid and sphingolipid metabolism to PD pathogenesis(37). PD pathogenic mechanisms linked to lipid metabolism include oxidative stress, inflammation and immune system signaling, pro-apoptotic processes, and interaction with a-synuclein biology, among others. Furthermore, alterations in serum, plasma, and brain measured phospholipids and sphingolipids have been widely reported in PD(38). For instance, LysoPC(18:1), implicated in our MWAS, has also been found at higher levels in the substantia nigra in animal models of PD(39). Interestingly, Cer(d18) metabolites, one of which was observed at higher intensity among our PD patients, have also been associated with physical frailty among older adults(40).
Our study is the largest untargeted high-resolution metabolomics study of PD to date, with metabolic profiles from independent discovery and replication case-control study populations allowing for validation of associated features. However, a notable limitation of the untargeted LC-MS technology is feature annotation. High-resolution LC-MS metabolomics provide metabolite features, many of which are not identified and can only be annotated based on m/z and retention time parameters from large databases (e.g. HMDB) and with consideration of feature correlation structures. While this does allow high confidence annotation, future research will be needed to identify features with certainty. For instance, para-, ortho-, and meta-cresol are all isomers. We have labeled the cresol metabolite as p-cresol due to co-occurrence with the p-cresol metabolites p-cresol sulfate and p-cresol glucuronide and because of the three exogenous metabolites, it is produced in humans via gut microbes. However, future studies will be needed to resolve the isomers. Furthermore, one-to- many matching and no matching add further uncertainty to feature annotation. Many of the features associated with PD in our MWAS, including some of the most significantly associated features, could not be annotated. Additionally, the metabolome measurements were based on a single blood-draw. Future longitudinal studies will be very informative in disentangling which if any metabolites implicated here are causally related to PD versus disease progression or reactive mechanisms.
In conclusion, based on untargeted high resolution, serum metabolic profiling from LC-MS, we have implicated over 200 individual metabolite features in PD along with multiple metabolic pathways. Several hits implicated pathways known to be disrupted in PD, including amino acid and lipid metabolism. We also present many novel findings, including for itaconate, connecting impaired anti-inflammatory signaling through immunometabolism, and three p-cresol metabolites, linking gut microbial activity to PD.
Data Availability
Data produced in the present study are available upon reasonable request to the authors
Footnotes
Conflicts of Interest: None
Funding: National Institute on Aging (K01AG07204401), National Institute of Environmental Health Sciences (grant number R21ES032593, 2R01ES010544, R21ES024356, R00ES028743, 5P30ES019776)