Thank you for visiting nature.com. You are using a browser version with limited support for CSS. To obtain the best experience, we recommend you use a more up to date browser (or turn off compatibility mode in Internet Explorer). In the meantime, to ensure continued support, we are displaying the site without styles and JavaScript.
npj Parkinson's Disease volume 9, Article number: 165 (2023 ) Cite this article citric acid
Gut microbial proteolytic metabolism has been reportedly altered in Parkinson’s disease (PD). However, the circulating aromatic amino acids (AAA) described in PD are inconsistent. Here we aimed to investigate plasma AAA profiles in a large cohort of PD patients, and examine their correlations with clinical severity and gut microbiota changes. We enrolled 500 participants including 250 PD patients and 250 neurologically normal controls. Plasma metabolites were measured using liquid chromatography mass spectrometry. Faecal samples were newly collected from 154 PD patients for microbiota shotgun metagenomic sequencing combined with data derived from 96 PD patients reported before. Data were collected regarding diet, medications, and motor and non-motor symptoms of PD. Compared to controls, PD patients had higher plasma AAA levels, including phenylacetylglutamine (PAGln), p-cresol sulfate (Pcs), p-cresol glucuronide (Pcg), and indoxyl sulfate (IS). Multivariable linear regression analyses, with adjustment for age, sex, and medications, revealed that the plasma levels of PAGln (coefficient 4.49, 95% CI 0.40–8.58, P = 0.032) and Pcg (coefficient 1.79, 95% CI 0.07–3.52, P = 0.042) positively correlated with motor symptom severity but not cognitive function. After correcting for abovementioned potential confounders, these AAA metabolites were also associated with the occurrence of constipation in PD patients (all P < 0.05). Furthermore, plasma levels of AAA metabolites were correlated with the abundance of specific gut microbiota species, including Bacteroides sp. CF01-10NS, Bacteroides vulgatus, and Clostridium sp. AF50-3. In conclusion, elevated plasma AAA metabolite levels correlated with disease characteristics in PD, suggesting that upregulated proteolytic metabolism may contribute to the pathophysiology of PD.
Parkinson’s disease (PD) is a common neurodegenerative disease characterized by various progressive motor and non-motor symptoms1. Neuropathological evidence reveals that, at least in some patients, the pathognomonic PD pathology termed the Lewy body may initiate from the enteric nervous system, and then propagate to the brain via cell-to-cell transmission2. The onset of gastrointestinal dysfunction, especially constipation, may precede the motor symptoms of PD by decades3. Many reports have shown altered gut microbiomes in PD patients compared to unaffected controls4. These observations suggest that gut microenvironment alterations may trigger the PD process though the complicated gut–brain axis5. The gut microbiota produces thousands of small molecules and metabolites that accumulate in the gastrointestinal system or reach distant organs through systemic circulation. Intestinal microorganisms interact with the host through secreted metabolites that act as signaling molecules in host–microbiome crosstalk, potentially modulating immune responses, endocrine secretion, and neurotransmission5. Understanding of how gut metabolites affect gut microbiome–host interactions in the PD process could shed light for future developing gut microbiome-based therapeutic interventions for PD.
The gut microbiota impact host metabolism by contributing enzymes, which are not encoded by the human genome, for the breakdown of dietary components. For example, the gut microbiota can fermentation carbohydrates to short chain fatty acids, proteins to aromatic amino acid (AAA) metabolites, lipids to bile acids and can synthesize certain vitamins, notably vitamin K and B. AAA metabolites—including phenylalanine (PHE), tyrosine, and tryptophan, have been reported to regulate immune, inflammatory, metabolic, and neuronal responses in the gut and brain6. Through redox pathways, PHE can be catabolized into phenylacetic acid (PAA), phenylpropionic acid (PPA), and phenyllactic acid (PLA)6. Furthermore, absorbed PAA can be conjugated by glutamine in the liver to form phenylacetylglutamine (PAGln). The dopamine precursor tyrosine is degraded by anaerobic bacteria to form p-cresol, which is metabolized into p-cresol sulfate (Pcs) and p-cresol glucuronide (Pcg) in the liver. These molecules are mainly excreted in urine, and are considered uremic toxins that promote inflammation and generate oxidative stress6. In contrast, indole derivatives produced from tryptophan by colonic bacteria generally have beneficial anti-inflammatory and anti-oxidative effects—with the exception of indoxyl sulfate (IS), which is a uremic toxin6. An aberrant degradative pathway of AAA is associated with several systemic disorders, including chronic renal disease, cardiovascular disorder, and neurodegeneration, including PD7,8. Some studies show higher concentrations of AAA metabolites, including PAGln and p-cresol, in the plasma and cerebrospinal fluid (CSF) of PD patients9,10; however, contradictory results exist11,12. Of note, a recent study observed reduced plasma AAA profiles in PD patients compared to controls and, furthermore, the levels were negatively correlated with the motor symptom severity of PD11.
In addition, some AAA metabolites, such as phenylalanine and phenylacetic acid, also have direct dietary sources and not solely derived from gut microbiota fermentation. Phenylalanine is present in a variety of protein-rich foods, including milk, eggs, and meat. It is also available as a dietary supplement and the artificial sweetener aspartame also contains phenylalanine. Phenylacetic acid, characterized by its sweet and floral flavors, can be detected in various food items such as hyssop, cowpea, endive, and shea tree. Hence, in the current study, we comprehensively examined plasma AAA levels in a large cohort of PD patients and unaffected controls matched for age, sex and diet habit. We further investigated how these AAA levels were associated with clinical severity (both motor and non-motor aspects), and gut microbiota changes, to delineate amino acid metabolism in the PD process.
This study enrolled a total of 500 participants, including 250 PD patients (67.4 ± 7.7 years old, 65.6% male) and 250 normal controls (67.4 ± 7.5 years old, 66.0% male). The participants’ clinical characteristics are listed in Table 1. The groups did not significantly differ in demographic distribution or laboratory results, including age, sex, body mass index, dietary pattern, medical co-morbidities, serum creatinine levels and dietary intake patterns measured by using the FFQ (Supplementary Table 1). MMSE scores were higher in control participants compared to PD patients (P < 0.01). Constipation was significantly more frequent in PD patients than controls (52% vs 41%, P = 0.01). In the PD group, the average Hoehn-Yahr stage was 3.2 ± 0.9 and mean MDS-UPDRS part III motor score was 35.0 ± 9.9. All PD patients received dopaminergic medication, with an average LEDD of 649.2 ± 206.0 mg/day.
Compared to control participants, PD patients exhibited higher plasma concentrations of PHE oxidative metabolites, including PAA (113.7 ± 82.9 vs 81.9 ± 106.1 ng/mL, P < 0.001) and the downstream metabolite PAGln (1024.3 ± 1025.6 vs 849.1 ± 2778.0 ng/mL, P < 0.001); other PHE-derived products, including PLA (30.3 ± 18.7 vs 26.7 ± 15.1 ng/mL, P = 0.015) and PPA (61.5 ± 122.7 vs 28.8 ± 40.6 ng/mL, P < 0.001); tyrosine-related metabolites, including Pcs (3265.1 ± 2899.3 vs 2110.7 ± 2603.1 ng/mL, P < 0.001) and Pcg (41.5 ± 61.5 vs 25.2 ± 101.8 ng/mL, P < 0.001); and the tryptophan derivative IS (1132.2 ± 1380.6 vs 1059.8 ± 2838.7, P < 0.001) (Table 1, Fig. 1A). We found that age was not correlated with levels of the individual AAA subtypes in either the control or PD group (Supplementary Fig. 1A, B).
A Violin plot shows plasma concentrations of microbial metabolites, with the solid and dashed horizontal lines representing the median value and interquartile range, respectively. B The accuracy of PD prediction with various models including age, sex, and/or different microbial metabolites. AUC, area under the receiver operating characteristic curves. ns, non-specific, *P < 0.05, **P < 0.01. ***P < 0.001.
We next examined whether plasma concentrations of AAA metabolites could distinguish between PD patients and normal controls. Multivariable logistic regression models, adjusted for age and sex, revealed that PD occurrence was associated with higher plasma levels of PAA (OR 4.00, 95% CI 2.52–6.53, P < 0.001), PAGln (OR 3.73, 95% CI 2.43–5.84, P < 0.001), PLA (OR 3.03, 95% CI 1.25–7.50, P = 0.015), PPA (OR 1.75, 95% CI 1.36–2.27, P < 0.001), Pcs (OR 2.29, 95% CI 1.74–3.07, P < 0.001), Pcg (OR 2.07, 95% CI 1.65–2.61, P < 0.001), and IS (OR 3.02, 95% CI 1.79–5.24, P < 0.001) (Supplementary Table 2). ROC curve analysis was performed to evaluate the accuracy of predicting PD occurrence based on plasma AAA metabolites. The prediction accuracy, expressed as the AUC, improved from 0.51 in the basic model (age and sex; 95% CI 0.46 – 0.56, P = 0.745), and improved to 0.69 (95% CI 0.55–0.65, P < 0.001) with addition of the PHE-associated metabolites PHE, PAGln, PLA, and PPA (Fig. 1B). The AUC further improved to 0.71 (95% CI 0.66–0.75, P < 0.001) in the full model containing age, sex, and all AAA proteolytic metabolites, including PHE, PAGln, PLA, PPA, Pcg, and IS (Fig. 1B).
Next, we examined whether AAA-related metabolites in systemic circulation reflected disease severity. Compared to patients with early-stage PD, those with advanced-stage PD exhibited higher plasma concentrations of PAGln (1100.0 ± 1100.2 vs 688.3 ± 467.7 ng/mL, P = 0.002), Pcs (3518.1 ± 3006.3 vs 2143.0 ± 2040.8 ng/mL, P = 0.003), Pcg (46.3 ± 65.8 vs 19.8 ± 29.1 ng/mL, P = 0.006), and IS (1177.5 ± 1466.4 vs 931.4 ± 892.1 ng/mL, P = 0.018) (Fig. 2).
A Violin plot shows the plasma levels of microbial metabolites in patients with early-stage and advanced-stage PD. Solid and dashed horizontal lines represent the median value and interquartile range, respectively. B Violin plots illustrate the comparison of plasma levels of microbial metabolites between PD patients with and without constipation. Solid and dashed horizontal lines indicate the median value and interquartile range, respectively. ns, non-specific. The comparison was conducted employing the Student’s t-test. *P < 0.05, **P < 0.01. ***P < 0.001.
Since antiparkinsonian medications were reported to affect the gut microbiota structure and microbial metabolite concentrations13,14, we examined the potential effects of antiparkinsonian medication on plasma levels of AAA metabolites. Spearman’s rank correlation test showed that most AAA metabolites had some degree of positive association with LEDD (Supplementary Fig. 2A). We further explored which types of antiparkinsonian medications affected plasma AAA levels. Among PD patients, users of dopaminergic agents and trihexyphenidyl exhibited lower plasma levels of PAGln, Pcs, Pcg, and IS compared to non-users (all P < 0.05), while users of catechol-O-methyltransferase inhibitor and amantadine had increased levels of PAGln, PAA, Pcs, Pcg, and IS compared to non-users (all P < 0.001) (Supplementary Table 3). Therefore, we adjusted for antiparkinsonian medication usage in the following further analyses.
Individual multivariable linear regression models analysing plasma levels of individual AAA metabolites and clinical severity adjusted for age, sex, occurrence of constipation, and anti-parkinsonian medication usage revealed that MDS-UPDRS part III motor score positively correlated with plasma concentrations of PAGln (coefficient 4.49, 95% CI 0.40–8.58, P = 0.032) and Pcg (coefficient 1.79, 95% CI 0.07–3.52, P = 0.042) (Supplementary Table 4 and summarized in Table 2). Consistently, further scatter plot analyses revealed that MDS-UPDRS part III motor scores significantly correlated with plasma levels of PAGln (ρ = 0.183, P = 0.004), PAA (ρ = 0.135, P = 0.033), Pcs (ρ = 0.165, P = 0.009), Pcg (ρ = 0.163, P = 0.010) and IS (ρ = 0.181, P = 0.004) (Supplementary Fig. 2B). However, spearman’s rank correlation analysis only showed that MMSE score tended to be inversely associated with plasma levels of PAGln (ρ = −0.124, P = 0.05) and IS (ρ = − 0.145, P = 0.02). After adjusting for age, sex, the occurrence of constipation, and anti-parkinsonian medication usage revealed that MMSE scores were not correlated with plasma concentrations of AAA metabolites (Table 2). The effects of other covariates, including age, sex, constipation, and anti-parkinsonism medications in the multiple linear regression analyses for individual AAA metabolite level and clinical disease severity were shown in the Supplementary Table 4.
Constipation is the most common non-motor symptom in PD, and may affect gut microbiota composition and the gut metabolites in systemic circulation. In the control group, we found that the plasma levels of AAA metabolites did not significantly differ between constipated and non-constipated participants (Supplementary Fig. 3). However, PD patients with constipation had higher plasma levels of individual AAA metabolites than PD patients without constipation (all P < 0.01, Fig. 2B). After adjustment for age, sex, disease severity and LEDD, constipation remained significantly associated with higher plasma concentrations of PHE (P = 0.02), PAGln (P = 0.03), PLA (P = 0.02), Pcs (P = 0.008) and IS (P = 0.03) within the PD group (Supplementary Table 5).
We next generated heat-maps depicting the Spearman correlations between the relative abundance of different bacterial species and individual AAA subtypes in plasma. We found that PAGln level was positively associated with the abundance of Bacteroides sp. CF01-10NS, Bacteroides vulgatus, and Clostridium sp. AF50-3, and negatively correlated with the abundance of Alistipes species (Fig. 3A). Plasma PAA, which is in the same pathway as PAGln, had similar associations with abundances of similar bacterial species, and was also positively correlated with abundance of Faecalibacterium sp. OM04-11BH. Another PHE reductive product, PPA, was positively related to abundance of Clostridium perfringens and Clostridium spiroforme. Plasma Pcs and Pcg are reportedly produced by similar bacterial species as PAGln15. Consistently, we found that plasma levels of Pcs and Pcg were positively correlated with the abundance of Bacteroides vulgatus and Bacteroides finegoldii, and negatively correlated with various Alistipes species. Bacteroides species were also positively correlated with plasma IS level (Fig. 3B). Together, these findings indicated that altered abundance of specific gut microbiota species were correlated with plasma levels of AAA metabolites, which were associated with PD clinical severity.
A Individual phenylalanine-related microbial metabolites, including phenylacetic acid (PAA), phenylpropionic acid (PPA) and phenylacetylglutamine (PAGln), are correlated with specific bacterial species. B Tryptophan and indole derivative-metabolites, including p-cresol sulfate (Pcs), p-cresol glucuronide (Pcg) and Indoxyl sulfate (IS), are correlated with specific bacterial species. The r values are represented by gradient colors, with red cells indicating positive correlations, and blue cells indicating negative correlations. Gut microbiota analysis was only performed in patients with PD. Among 250 PD patients, the fecal shotgun metagenomic data derived from 96 PD patients was reported before35. The correlations between the relative abundance of gut microbiota and plasma AAA levels were performed using Spearman’s rank sum tests with Benjamini and Hochberg adjustment for multiple tests. *P < 0.05, **P < 0.01. ***P < 0.001.
In a large cohort of participants, we performed simultaneous analyses of various AAA metabolites in systemic circulation, and examined their associations with gut microbiota detected by shotgun metagenomic sequencing, to comprehensively delineate the alterations of microbial AAA metabolites in relationship to PD occurrence and severity. We found that PD patients had higher plasma concentrations of AAA metabolites compared to controls. Moreover, among PD patients, these AAA metabolite levels were higher in patients with constipation than those without constipation. After adjustment for confounders, the increased plasma levels of PAGln and Pcg were correlated with motor symptom severity, and associated with specific gut microbiota changes.
PAGln level was increased in PD patients compared to controls, and was correlated with motor symptom severity after adjustment for confounders. The PAGln precursor PAA was also elevated in PD patients. PAGln is a conjugated product of glutamate and PAA produced from PHE through gut microbial proteolytic processes associated with the commensal bacteria Clostridium, which was detected in our cohort16. Elevated plasma PAGln levels are reported in systemic and neurological diseases, including chronic renal disease, cardiovascular syndrome, and diabetic neuropathy16. In line with our findings, two recent untargeted metabolomics studies with large multi-cohort study designs also revealed elevated plasma levels of PAGln in PD patients compared to controls and other neurological diseases9,10. However, some studies have not detected different plasma PAGln levels between PD patients and controls11,12. These inconsistent results may be explained by diet habits and medical co-morbidities, since PAGln is a deleterious biomarker for cardiovascular disorder, DM, and chronic renal failure16. Therefore, a control group well-matched for age, sex, and these other potential confounders is critical for elucidating the role of PAGln in the PD process. In our current study, the PD and control groups were matched for these confounders, and anti-parkinsonian medications were considered as an adjustment covariate in the analyses. Our results suggested that PAGln and its related pathway may be a mechanistic link in the PD process. Notably, recent functional studies show that PAGln—which is upregulated in PD patients and promotes PD neuropathology through neuroinflammation in PD rodent models—exacerbates inflammatory cascades by activating toll-like receptor 4-associated pathways17. Furthermore, PAGln administration in a mouse model exacerbates colitis by upregulating coagulation-related biological processes, and plasma PAGln level is elevated in patients with Crohn’s disease18,19. Crohn’s disease is associated with PD, since both diseases shared a common genetic risk factor, leucine-rich repeat kinase 2 (LRRK2), as discovered in genome-wide association studies20. A recent epidemiology meta-analysis study demonstrated that Crohn’s disease patients have an 28% increased risk of PD compared to controls21. Furthermore, PAGln can cross the blood–brain barrier and exert neurotoxicity, thereby impairing neuronal activity22. These observations and our present results suggest that plasma PAGln level is a novel gut-oriented biosignature for PD. Further studies are warranted to investigate the molecular mechanism of PAGln in PD pathophysiology.
Our results also showed elevated plasma levels of Pcs and Pcg in PD patients compared to controls, which is in line with a recent untargeted metabolomic study23. Pcs and Pcg are both conjugated metabolites of p-cresol, a deleterious bacterial degradative product from tyrosine6. Elevated Pcs levels in CSF and plasma have been reported in several neurological disorders, including PD, Alzheimer’s disease, and multiple sclerosis24. Pcs impairs neuronal physiological function and exerts cytotoxic effect through generation of oxidative stress and upregulation of pro-inflammatory cytokines22. Pcg and Pcs have different biological functions. Some studies have reported that Pcg can prevent lipopolysaccharide-induced blood–brain barrier disruption25, while others show that Pcg suppresses mitochondrial function and aggravates the Pcs-induced oxidative burst of leukocytes, synergistically promoting inflammatory responses26. We observed that plasma Pcg level positively correlated with motor symptom severity in PD patients. Functional studies are needed to explore the molecular mechanisms underlying how Pcs and Pcg participate in the PD disease mechanism.
Constipation is among the most common non-motor and prodromal PD symptoms, and exhibits increasing prevalence with PD progression27. A major cause of constipation in PD is prolonged colonic transition time, likely resulting from α-synuclein accumulation in the enteric nervous system27. Impaired colonic peristalsis with delayed transit time is associated with gut dysbiosis, which results in increased colonic protein fermentation and a shift of the colonic microbiota metabolism away from carbohydrate fermentation and toward proteolysis15. In previous studies, more abundant Ruminoccoccaceae, Oscillospira, Lachnospiraceae, and Clostridiales, and higher plasma PAGln and Pcs levels, are observed in patients with impaired colonic motility, and are correlated with severity of constipation in PD9,15. Accordingly, we found that plasma PAGln and Pcs levels were higher in constipated than non-constipated PD patients, after adjustment for potential confounders, reinforcing that the colonic metabolism was skewed towards proteolysis in PD. In our cohort, these AAA proteolytic metabolites were associated with specific microbial species that have pro-inflammatory effects. Higher levels of PAGln, Pcs, Pcg, and IS were all associated with increased abundance of Bacteroides species, including Bacteroides vulgatus, which is a commensal bacteria that induces colitis in rodent models and is positively correlated with severity of inflammatory bowel disease28,29. More abundant Bacteroides species have also been reported in PD, although though are some conflicting results30. Low-grade intestinal inflammation could lead to intestinal epithelial barrier dysfunction and increased infiltration of luminal microbial metabolites, including AAA derivatives, into systemic circulation, which may partially explain the increased plasma AAA metabolites in our cohort31. Although we observed elevated plasma AAA levels in PD patients compared to controls, especially among those with advanced-stage PD and those with constipation symptoms, we did not find that plasma AAA concentrations were correlated with MMSE scores. One possible reasons is that the cognitive function of most of our enrolled PD patients did not reach dementia status, with a mean MMSE of >26. Another possible explanation may reflect the diverse pathophysiologic processes of motor and cognitive decline in PD32. Further studies enrolling PD patients with different severity of cognitive dysfunction are needed to explore the role of plasma AAA levels in PD dementia.
To our knowledge, our study is the largest cohort to comprehensively investigate the composition of AAA proteolytic metabolites in a cohort of PD patients and well-matched controls, and delineate their associations with gut microbiota. However, our study has several limitations. First, the cross-sectional study design limits our inference of causal relationship of AAA metabolite changes and PD. In addition, constipation symptoms in the current study may come from heterogeneous etiologies, including as a prodromal symptom of PD process, a co-morbidity due to motor dysfunction and gastrointestinal side effects of anti-parkinsonism medications (ex: anti-cholinergic agents and dopamine agonizts). A future prospective follow-up study including patients in the prodromal stage and drug naive patients are needed to correlate the changes of gut microbiota, plasma AAA levels and the symptoms of constipation. Further in vitro and in vivo assays are also needed to provide molecular insights into AAA metabolites in PD pathophysiology. Second, gut microbiota analysis was only performed in PD patients in this study. In this regard, as we did not have shotgun metagenomic data from control participants, which hampers us from performing functional analysis to compare the pathway differences which are highly correlated with AAAs. Further studies enrolling participants with concomitant information of fecal metagenomic shotgun sequencing data and plasma metabolite details are warranted to identify the specific bacterial species, which are highly correlated with AAAs, and the relevant functional pathways in the bacterial genomes. Third, although we showed that increased AAA metabolites were associated with gut microbiota, we cannot determine whether the correlations resulted from increased microbial production, enteric barrier dysfunction, or other mechanisms. Future research should include a concomitant leaky gut test. Fourth, we did not examine the inflammatory markers in feces or plasma to confirm our hypothesis that the AAA metabolites were associated with pro-inflammatory status in the gut and systemic circulation. Finally, we assessed cognitive function only with MMSE, a simple measurement of global cognitive function. Future studies should include detailed neuropsychological tests to evaluate individual cognitive domains, and participants with diverse cognitive dysfunction, to further assess the correlation between plasma AAA levels and cognitive decline in patients with PD.
In conclusion, our results revealed that AAA proteolytic metabolites in systemic circulation were upregulated in PD patients compared to controls. These metabolites were also associated with specific colonic bacterial species, and with motor symptom severity after adjusting for confounders, especially PAGln and Pcg. These findings suggest that dysregulation of microbial-associated proteolytic metabolism and AAA metabolites may contribute to PD pathophysiology.
PD patients and neurologically normal controls were enrolled from the movement disorder clinic of National Taiwan University Hospital. PD was diagnosed according to the UK PD Society Brain Bank Clinical Diagnostic Criteria33. As controls, we enrolled neurologically unaffected persons accompanying the PD patients. Participants were excluded if they had renal or liver disease; history of inflammatory bowel disease, irritable bowel syndrome, colitis, or colon cancer; had used antibiotics or probiotic supplements within 3 months of enrollment; or were vegetarians. Constipation was diagnosed based on the standard Rome IV gastrointestinal disorders diagnostic criteria34. Among the enrolled participants, 96 patients with PD had undergone fecal shotgun metagenomic sequencing, which was previously published35.
The research protocol was reviewed by the Institutional Research Board Committee at National Taiwan University Hospital. All participants provided written informed consent before entering the study.
A comprehensive dietary history was collected using the Food Frequency Questionnaire (FFQ)36. Motor symptom severity was assessed during the “on” phase of PD using Hoehn and Yahr staging37 and Movement Disorder Society-Unified Parkinson’s Disease Rating Scale (MDS-UPDRS) part III motor scores38. Hoehn and Yahr stage of <3 was considered early-stage PD, and ≥3 as advanced-stage PD. Cognitive function was examined using the Mini-Mental State Examination (MMSE). Patients’ anti-PD medication dosage during the study period was converted into a levodopa equivalent daily dosage (LEDD).
After participants fasted for at least 8 h, 10 µL venous blood was collected, centrifuged, and stored at −80 °C before further analysis. To quantify microbial metabolites, 500 ng/mL internal standards of all AAA-derived gut microbial metabolites, in 200 µL methanol, was added to 50 µL calibration standard solution and to 50 µL plasma, for protein precipitation. Samples were mixed with Geno/Grinder 2020 (OPS Diagnostics, LLC, NJ, USA) at 1000 rpm for 3 min, and then centrifuged at 15000 rcf for 5 min at 4 °C. We then filtered 200 µL of supernatant through a 0.2-µm filter (RC-4, Sartorius, Gottingen, Germany), followed by analysis with liquid chromatography-mass spectrometry (LC-MS).
LC separation was performed using an Agilent 1290 UHPLC system and Waters ACQUITY UPLC HSS T3 column (2.1 × 100 mm, 1.8 μm). The mobile phase was water with 2 mM ammonium acetate as solvent A, and water:acetonitrile (2:3, v-v) as solvent B. The flow rate was 0.3 mL/min, and the injection volume was 5 µL. The LC gradient was as follows: 0–1.5 min, 0–25% B; 1.5–4.5 min, 25% B; 4.5–4.6 min, 25–100% B; and 4.6–7 min, 100% B. Using an Agilent 6495 triple quadrupole mass spectrometer (Agilent Technologies, Santa Clara, CA) with an ESI source, we analysed the metabolites in negative ion mode, with a capillary voltage of 3500 V and nozzle voltage of 2000 V. For the dry gas and sheath gas, respectively, the temperatures were 290 °C and 325 °C, and the flow rates were 11 L/min and 10 L/min. The nebulizer was set to 40 psi. Spectra were acquired in multiple reaction monitoring mode. Analyte concentrations were calculated from calibration curves, with the peak area ratio of each analyte to its corresponding internal standard.
For data imputation, we imputed the missing values based on the average level of the individual AAA metabolite across all participants. If the mean metabolite measure across samples was >100,000, then zero was imputed, otherwise one half of the minimum measure for that sample was imputed39.
Gut microbiota analysis was only performed in patients with PD. Among 250 PD patients, the fecal shotgun metagenomic data derived from 96 PD patients was reported before with plasma samples collected at the same time35. The remaining 154 PD patients were newly recruited for assessing gut microbiota and plasma AAA metabolites in the current study. There was no significant batch effect as the beta diversity between these two batches of microbiota analyses did not reveal significant differences (Supplementary Fig. 4). Fecal sample collection, fecal DNA extraction, and gut microbiota analysis were performed35.
Fresh fecal samples were obtained from each participant and immediately preserved in stool specimen collection tubes containing DNA stabilizer (Sarstedt). These samples were promptly flash-frozen on dry ice and stored at −80 °C until further analysis. The total fecal DNA was extracted using the QIAamp DNA Stool Mini Kit (Qiagen, Hilden, Germany). To prepare the DNA for shotgun metagenomic sequencing, several steps were taken. First, the purity of the DNA was assessed using an Epoch Microplate Spectrophotometer (BioTek, USA). Then, the quantity of DNA was measured using the Quant-iT PicoGreen dsDNA Assay (ThermoFisher Scientific, USA). Subsequently, libraries were constructed using Illumina DNA Prep kits (Illumina, USA), targeting insert sizes of ~350 base pairs. The prepared libraries, with the specified insert sizes, were subjected to sequencing on an Illumina NovaSeq 6000 sequencer, utilizing S4 flow cells. For both case and control samples, the metagenomic Illumina paired-end reads underwent trimming, which was performed using Trimmomatic v0.39 with the following parameters: “PE -phred33 ILLUMINACLIP:TruSeq3-PE.fa:2:30:10 LEADING:12 TRAILING:12 SLIDINGWINDOW:4:15 MINLEN:36.”40.
The trimmed paired-end reads from the samples were subsequently subjected to a prokaryotic species profiling analysis using Kraken2 v2.0.6-beta 4941. To enhance the efficiency of this search, we employed the maxikraken2 database, which has a total size of 140 GB. This specialized database, constructed by Daniel Fischer at the Natural Resources Institute Finland and accessible at https://lomanlab.github.io/mockcommunity/mc_databases.html, was used in place of the default kraken2 database. Prior to investigating the relationship between gut microbiota and plasma AAA levels, we implemented a filtering process to exclude rare taxa, defined as those with a prevalence <10% among all samples and a mean abundance of <0.001%. Correlations between the relative abundance of gut microbiota and plasma AAA levels were assessed through Spearman’s rank sum tests, with Benjamini and Hochberg adjustment applied to account for multiple tests. Taxa with a P-value < 0.01 were selected for further analysis. To visualize the correlation between the relative abundance of gut microbiota and plasma levels of AAA, we utilized the R package “pheatmap”. All analyses were conducted using R software 3.4.142.
The codes that were used in the R software were shown as below.
pd500 <- read.csv(“pd500_clean.csv”, header=TRUE, sep = “,”) ##read AAA data
kraken_pd500 <- read.csv(“kraken_pd500.csv”, header=TRUE, row.names=1, sep = “,”) %>% as.matrix() ##read abundance table
## filter out rare taxa (prevelance < 10%, mean abundance < 0.001%)
prev_idx <- sapply(1:1974,function(x) sum(kraken_pd500[x,]>0)) > 0.1*ncol(kraken_pd500)
abun_idx <- sapply(1:1974,function(x) mean(kraken_pd500[x,])) > 0.00001
kraken_pd500_f <- kraken_pd500[(prev_idx & abun_idx),]
metabolite <- colnames(pd500)[3:11]
rho <- data.frame()
p_val <- data.frame()
for (j in 1: nrow(kraken_pd500_f)) {
test <- cor.test(log(pd500[,i + 2]),log(kraken_pd500_f[j,]),method = “spearman”)
rho[j,i] <- test$estimate
p_val[j,i] <- test$p.value
p_adj <- apply(p_val,2,function(x) p.adjust(x,“BH”))
colnames(rho) <- metabolite
colnames(p_val) <- paste0(metabolite,”_p”)
colnames(p_adj) <- paste0(metabolite,”_p_adj”)
result_simple <- rho
colnames(result_simple) <- metabolite
rownames(result_simple) <- rownames(kraken_pd500_f)
result <- cbind(rho,p_val,p_adj)
rownames(result) <- rownames(kraken_pd500_f)
result_pc <- result_simple[,6:9]
asterisk_pc <- cut(unlist(p_val[,6:9]), breaks=c(-Inf, 0.001, 0.01, 0.05, Inf),
label=c(“***“, “**“, “*“, “”)) %>% matrix(., ncol=4)
result_phe <- result_simple[,1:5]
asterisk_phe <- cut(unlist(p_val[,1:5]), breaks=c(-Inf, 0.001, 0.01, 0.05, Inf),
label=c(“***“, “**“, “*“, “”)) %>% matrix(., ncol=5)
idx_p_pc <- (p_val[,6] < 0.01) | (p_val[,7] < 0.01) | (p_val[,8] <0.01) | (p_val[,9] < 0.01)
result_pc <- result_pc[,2:4]
asterisk_pc <-asterisk_pc[,2:4]
pc_names <- lapply(
function(x) bquote(italic(.(x))))
pheatmap(result_pc[idx_p_pc,], color = colorRampPalette(c(“#33A1C9”, “white”, “brown3”))(50),
fontsize_row=11, cluster_cols=FALSE, display_numbers = asterisk_pc[idx_p_pc,],
labels_row = as.expression(pc_names))
idx_p_phe <- (p_val[,1] < 0.01) | (p_val[,2] < 0.01) | (p_val[,3] <0.01) | (p_val[,4] < 0.01) | (p_val[,5] < 0.01)
desired_order <- c(“PHE”,“PAA”, “PAGln”, “PPA”, “PLA”) ## Lin’s order 20230626
desired_res_phe <- result_phe[,desired_order]
pheatmap(result_phe[idx_p_phe,], color = colorRampPalette(c(“#33A1C9”, “white”, “brown3”))(50),
fontsize_row=7, cluster_cols=FALSE, display_numbers = asterisk_phe[idx_p_phe,])
#labels_col could not change order
function(x) bquote(italic(.(x))))
pheatmap(desired_res_phe[idx_p_phe,], color = colorRampPalette(c(“#33A1C9”, “white”, “brown3”))(50),
cluster_cols=FALSE, display_numbers = asterisk_phe[idx_p_phe,],
labels_row = as.expression(newnames)) # change spp names to italic
Continuous variables are expressed as mean ± standard deviation, and categorical factors as numbers and percentages. We examined homogeneity of variances with Levene test. Variables were compared using Student’s t-test or analysis of variance if fulfilling Gaussian distribution, and using the nonparametric Mann-Whitney U test or Fisher’s exact test if the variables violated homoscedasticity or assumption of normality. The microbial AAA metabolite concentrations were log-transformed to achieve normal distribution for further analyses. A logistic regression model was applied to examine the relationship between AAA levels and PD occurrence, with adjustment for age and sex. We used the area under the receiver operating characteristic (ROC) curve (AUC) to quantify the model’s diagnostic performance for exploring the ability of AAA to distinguish between PD patients and controls. Spearman’s rank correlation was used to examine correlations between plasma AAA levels and disease severity. To examine associations between AAA and motor or cognitive symptom severity, we used multivariable linear regression models, with MDS-UPDRS part III and MMSE scores as dependent factors. The effect size of the model was estimated with Cohen’s f2. Independent variables included the log-transformed AAA levels, age, sex, disease duration, and LEDD. The correlations between the relative abundance of gut microbiota and plasma AAA levels were performed using Spearman’s rank sum tests with Benjamini and Hochberg adjustment for multiple tests. The statistical analyses were majorly performed in MedCalc software version 19.0.3 (MedCalc Software bvba, Ostend, Belgium) while the shotgun sequencing and microbiota related bioinformatic analysis were conducted using R software (3.4.1, R Foundation for Statistical Computing, Vienna, Austria). P < 0.05 was considered significant.
The raw anonymized microbiome sequencing data were uploaded to the European Nucleotide Archive database (https://www.ebi.ac.uk/ena/browser/home) with the accession number PRJEB57770. Anonymized data not published within this article will be made available upon request.
All code for data analysis associated with the current study is provided in the methods.
Kalia, L. V. & Lang, A. E. Parkinson’s disease. Lancet Lond. Engl. 386, 896–912 (2015).
Braak, H. et al. Staging of brain pathology related to sporadic Parkinson’s disease. Neurobiol. Aging 24, 197–211 (2003).
Savica, R. et al. Medical records documentation of constipation preceding Parkinson disease: a case-control study. Neurology 73, 1752–1758 (2009).
Article CAS PubMed PubMed Central Google Scholar
Chen, S. J. & Lin, C. H. Gut microenvironmental changes as a potential trigger in Parkinson’s disease through the gut-brain axis. J. Biomed. Sci. 29, 54 (2022).
Article CAS PubMed PubMed Central Google Scholar
Tan, A. H., Lim, S. Y. & Lang, A. E. The microbiome-gut-brain axis in Parkinson disease - from basic research to the clinic. Nat. Rev. Neurol. 18, 476–495 (2022).
Liu, Y., Hou, Y., Wang, G., Zheng, X. & Hao, H. Gut microbial metabolites of aromatic amino acids as signals in host-microbe interplay. Trends Endocrinol. Metab. 31, 818–834 (2020).
Article CAS PubMed Google Scholar
Vanholder, R., Schepers, E., Pletinck, A., Nagler, E. V. & Glorieux, G. The uremic toxicity of indoxyl sulfate and p-cresyl sulfate: a systematic review. J. Am. Soc. Nephrol. 25, 1897–1907 (2014).
Article CAS PubMed PubMed Central Google Scholar
Sankowski, B. et al. Higher cerebrospinal fluid to plasma ratio of p-cresol sulfate and indoxyl sulfate in patients with Parkinson’s disease. Clin. Chim. Acta. 501, 165–173 (2020).
Article CAS PubMed Google Scholar
Cirstea, M. S. et al. Microbiota composition and metabolism are associated with gut function in Parkinson’s disease. Mov. Disord. 35, 1208–1217 (2020).
Article CAS PubMed Google Scholar
Willkommen, D. et al. Metabolomic investigations in cerebrospinal fluid of Parkinson’s disease. PLoS One 13, e0208752 (2018).
Article CAS PubMed PubMed Central Google Scholar
Zhang, Y. et al. Plasma branched-chain and aromatic amino acids correlate with the gut microbiota and severity of Parkinson’s disease. NPJ Park. Dis. 8, 48 (2022).
Yan, Z. et al. Alterations of gut microbiota and metabolome with Parkinson’s disease. Microb. Pathog. 160, 105187 (2021).
Article CAS PubMed Google Scholar
Hertel, J. et al. Integrated analyses of microbiome and longitudinal metabolome data reveal microbial-host interactions on sulfur metabolism in Parkinson’s disease. Cell Rep. 29, 1767–1777.e1768 (2019).
Article CAS PubMed PubMed Central Google Scholar
Hill-Burns, E. M. et al. Parkinson’s disease and Parkinson’s disease medications have distinct signatures of the gut microbiome. Move. Disord. 32, 739–749 (2017).
Roager, H. M. et al. Colonic transit time is related to bacterial metabolism and mucosal turnover in the gut. Nat. Microbiol. 1, 16093 (2016).
Article CAS PubMed Google Scholar
Nemet, I. et al. A cardiovascular disease-linked gut microbial metabolite acts via adrenergic receptors. Cell 180, 862–877.e822 (2020).
Article CAS PubMed PubMed Central Google Scholar
Fu, H., Kong, B., Zhu, J., Huang, H. & Shuai, W. Phenylacetylglutamine increases the susceptibility of ventricular arrhythmias in heart failure mice by exacerbated activation of the TLR4/AKT/mTOR signaling pathway. Int. Immunopharmacol. 116, 109795 (2023).
Article CAS PubMed Google Scholar
Feng, R. et al. Gut microbiome generated phenylacetylglutamine from dietary protein is associated with Crohn’s Disease and exacerbates colitis in mouse model possibly via platelet activation. J. Crohns. Colitis. 15, 774–795 (2023).
Ma, R. et al. A novel serum metabolomic panel for the diagnosis of Crohn’s disease. Inflamm. Bowel Dis. 88, 55–115 (2023).
Hui, K. Y. et al. Functional variants in the LRRK2 gene confer shared effects on risk for Crohn’s disease and Parkinson’s disease. Sci. Transl. Med. 10, 345–422 (2018).
Zhu, F. et al. The risk of Parkinson’s disease in inflammatory bowel disease: a systematic review and meta-analysis. Dig. Liver Dis. 51, 38–42 (2019).
Ntranos, A. et al. Bacterial neurotoxic metabolites in multiple sclerosis cerebrospinal fluid and plasma. Brain 145, 569–583 (2022).
Poesen, R. et al. Metabolism, protein binding, and renal clearance of microbiota-derived p-Cresol in patients with CKD. Clin. J. Am. Soc. Nephrol. 11, 1136–1144 (2016).
Article CAS PubMed PubMed Central Google Scholar
Kalecky, K., German, D. C., Montillo, A. A. & Bottiglieri, T. Targeted metabolomic analysis in Alzheimer’s disease plasma and brain tissue in non-hispanic whites. J. Alzheimers Dis. 86, 1875–1895 (2022).
Article CAS PubMed PubMed Central Google Scholar
Stachulski, A. V., Knausenberger, T. B., Shah, S. N., Hoyles, L. & McArthur, S. A host-gut microbial amino acid co-metabolite, p-cresol glucuronide, promotes blood-brain barrier integrity in vivo. Tissue Barriers 11, 2073175 (2023).
Zhu, S., Rong, Y. & Kiang, T. K. L. Effects of p-Cresol on oxidative stress, glutathione depletion, and necrosis in HepaRG Cells: comparisons to other uremic toxins and the role of p-Cresol glucuronide formation. Pharmaceutics 13, 52–88 (2021).
Knudsen, K., Krogh, K., Ostergaard, K. & Borghammer, P. Constipation in parkinson’s disease: subjective symptoms, objective markers, and new perspectives. Mov. Disord. 32, 94–105 (2017).
Bloom, S. M. et al. Commensal bacteroides species induce colitis in host-genotype-specific fashion in a mouse model of inflammatory Bowel disease. Cell Host. Microb. 9, 390–403 (2011).
Mills, R. H. et al. Multi-omics analyses of the ulcerative colitis gut microbiome link bacteroides vulgatus proteases with disease severity. Nat. Microbiol. 7, 262–276 (2022).
Article CAS PubMed PubMed Central Google Scholar
Keshavarzian, A. et al. Colonic bacterial composition in Parkinson’s disease. Mov. Disord. 30, 1351–1360 (2015).
Article CAS PubMed Google Scholar
Yang, X. et al. Parkinson’s disease is associated with impaired gut-blood barrier for short-chain fatty acids. Mov. Disord. 37, 1634–1643 (2022).
Article CAS PubMed Google Scholar
Walker, L., Stefanis, L. & Attems, J. Clinical and neuropathological differences between Parkinson’s disease, Parkinson’s disease dementia and dementia with Lewy bodies - current issues and future directions. J. Neurochem. 150, 467–474 (2019).
Article CAS PubMed Google Scholar
Hughes, A. J., Daniel, S. E., Kilford, L. & Lees, A. J. Accuracy of clinical diagnosis of idiopathic Parkinson’s disease: a clinico-pathological study of 100 cases. J. Neurol. Neurosurg. Psychiatr. 55, 181–184 (1992).
Drossman, D. A. Functional gastrointestinal disorders: history, pathophysiology, clinical features and Rome IV. Gastroenterology 92, e80091 (2016).
Chen, S. J. et al. Association of fecal and plasma levels of short-chain fatty acids with gut microbiota and clinical severity in patients with Parkinson disease. Neurology 98, e848–e858 (2022).
Article CAS PubMed PubMed Central Google Scholar
Chiu, T. H. et al. Relative validity and reproducibility of a quantitative FFQ for assessing nutrient intakes of vegetarians in Taiwan. Public Health Nutr. 17, 1459–1466 (2014).
Hoehn, M. M. & Yahr, M. D. Parkinsonism: onset, progression and mortality. Neurology 17, 427–442 (1967).
Article CAS PubMed Google Scholar
Goetz, C. G. et al. Movement disorder dociety-sponsored revision of the unified Parkinson’s disease rating scale (MDS-UPDRS): scale presentation and clinimetric testing results. Mov. Disord. 23, 2129–2170 (2008).
Sreekumar, A. et al. Metabolomic profiles delineate potential role for sarcosine in prostate cancer progression. Nature 457, 910–914 (2009).
Article CAS PubMed PubMed Central Google Scholar
Bolger, A. M. et al. Trimmomatic: a flexible trimmer for Illumina sequence data. Bioinformatics 30, 2114–2120 (2014).
Article CAS PubMed PubMed Central Google Scholar
Wood, D. E. et al. Improved metagenomic analysis with Kraken 2. Genom. Biol. 20, 257 (2019).
Team, R.C. R: A Language and Environment for Statistical Computing (Foundation for Statistical Computing Vienna, 2020).
We thank the patients for participating in this study. We also thank the staff of the Second Core Lab, Department of Medical Research, National Taiwan University Hospital, for technical support during the study. We are grateful to the funding supports by grants from National Taiwan University (NTU-110-A-CC-5400-64841) and National Health Research Institutes (NHRI-EX112-11136NI).
These authors contributed equally: Szu-Ju Chen, Yu-Jun Wu.
Department of Neurology, National Taiwan University Hospital, College of Medicine, National Taiwan University, Taipei, Taiwan
Szu-Ju Chen & Chin-Hsien Lin
Department of Neurology, National Taiwan University Hospital Bei-Hu Branch, Taipei, Taiwan
Graduate Institute of Clinical Medicine, College of Medicine, National Taiwan University, Taipei, Taiwan
Szu-Ju Chen & Chieh-Chang Chen
School of Pharmacy, College of Medicine, National Taiwan University, Taipei, Taiwan
Division of Gastroenterology and Hepatology, Department of Internal Medicine, National Taiwan University Hospital, College of Medicine, National Taiwan University, Taipei, Taiwan
Chieh-Chang Chen, Jyh-Ming Liou & Ming-Shiang Wu
Graduate Institute of Biomedical Informatics, College of Medical Science and Technology, Taipei Medical University, Taipei, Taiwan
Yu-Wei Wu & Ching-Hua Kuo
The Metabolomics Core Laboratory, NTU Centers of Genomic and Precision Medicine, National Taiwan University, Taipei, Taiwan
Department of Pharmacy, National Taiwan University Hospital, Taipei, Taiwan
Institute of Molecular Medicine, College of Medicine, National Taiwan University, Taipei, Taiwan
You can also search for this author in PubMed Google Scholar
You can also search for this author in PubMed Google Scholar
You can also search for this author in PubMed Google Scholar
You can also search for this author in PubMed Google Scholar
You can also search for this author in PubMed Google Scholar
You can also search for this author in PubMed Google Scholar
You can also search for this author in PubMed Google Scholar
You can also search for this author in PubMed Google Scholar
Study concept and design: C.-H.K., C.-H.L. Acquisition of data: S.-J.C., Y.-J.W., C.-C.C., J.-M.L., M.-S.W, C.-H.K., C.-H.L. Analysis and interpretation of data: S.-J.C., Y.-J.W., C.-C.C. Drafting of the manuscript: S.J.C. Critical revision of the manuscript for important intellectual content: C.-H.K., C.-H.L. Obtained funding: C.-H.K., C.-H.L.
Correspondence to Ching-Hua Kuo or Chin-Hsien Lin.
The authors declare no competing interests.
Publisher’s note Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made. The images or other third party material in this article are included in the article’s Creative Commons license, unless indicated otherwise in a credit line to the material. If material is not included in the article’s Creative Commons license and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this license, visit http://creativecommons.org/licenses/by/4.0/.
Chen, SJ., Wu, YJ., Chen, CC. et al. Plasma metabolites of aromatic amino acids associate with clinical severity and gut microbiota of Parkinson’s disease. npj Parkinsons Dis. 9, 165 (2023). https://doi.org/10.1038/s41531-023-00612-y
DOI: https://doi.org/10.1038/s41531-023-00612-y
Anyone you share the following link with will be able to read this content:
Sorry, a shareable link is not currently available for this article.
Provided by the Springer Nature SharedIt content-sharing initiative
npj Parkinson's Disease (npj Parkinsons Dis.) ISSN 2373-8057 (online)
corn starch food Sign up for the Nature Briefing: Translational Research newsletter — top stories in biotechnology, drug discovery and pharma.