PD patients have higher calprotectin and lower SCFA levels in stool and lower CXCL8 levels in plasma compared to controls
Patient and control groups were similar with regard to basic demographics such as age, sex, and body mass index (BMI); however, as we knew from previous analyses of the same subjects (31, 32), the groups differed regarding medications, medical history, and various symptom scores. A higher percentage of controls reported a history of stroke and use of medications for high blood pressure and cholesterol, and PD patients scored higher on scales of non-motor symptoms, gastrointestinal problems, and constipation (Table 1). Contrasting inflammatory markers and SCFAs between PD patients and control subjects, patients had lower levels of butyric and propionic acid and higher levels of calprotectin in their stool (Additional File 3) and lower levels of CXCL8 in plasma (Additional File 3). When the data were stratified by sex, the differences were particularly prominent for butyric acid in males and plasma CXCL8 and stool calprotectin in females (Fig. 1, Additional File 3).
Table 1. Demographic and Clinical Details of Subjects
|
Control subjects
|
PD patients
|
p-value
|
Number of Subjects
|
56
|
55
|
|
Sex (% Male)
|
48.21
|
52.73
|
0.706
|
Age at stool collection (mean ± SD)
|
66.38 ± 6.73
|
67.63 ± 5.21
|
0.293
|
Body Mass Index (mean ± SD)
|
26.7 ± 3.56
|
27.63 ± 4.74
|
0.261
|
History of TIA/ischemic stroke (%)
|
37.50
|
5.56
|
< 0.001
|
Medication: ACE-inhibitor or
AT1 antagonist (% yes)
|
48.21
|
29.09
|
0.051
|
Medication: calcium channel blocker (% yes)
|
19.64
|
5.45
|
0.042
|
Medication: statin (% yes)
|
48.21
|
18.18
|
0.001
|
Levodopa equivalent daily dose
(LEDD mg; mean ± SD)
|
0 ± 0
|
608.61 ± 301.46
|
< 0.001
|
Non-Motor Symptoms Scale (NMSS) total (mean ± SD)
|
7.09 ± 6.35
|
48.36 ± 36.55
|
< 0.001
|
Rome III constipation subscore
(items 9-15; mean ± SD)
|
2.46 ± 3.16
|
7.45 ± 5.22
|
< 0.001
|
Rome III irritable bowel syndrome
criteria fulfilled (%)
|
7.14
|
36.36
|
< 0.001
|
p-values reflect Fisher’s exact test; SD–standard deviation, TIA–transient ischemic attack, ACE–Angiotensin-Converting Enzyme, AT1–angiotensin II type 1
Inflammatory markers in plasma and stool are not correlated, but stool SCFAs, zonulin, and stool and plasma cytokines are related to PD onset and symptom severity
The inflammatory markers measured from the same sample material (plasma or stool) were largely intercorrelated, but we found no correlations between measurements of the same markers in stool and plasma (Fig. 2). PCAs for combining the most correlated markers resulted in one principal component (PC) for stool (summarizing IL-10, IL-12p70, IL-13, TNF, IL-6, and IFNγ), and two plasma PCs (PC1: IL-10, IL-4, TNF; PC2: IL-1β, IL-2, IL-12p70, IL-13) (Additional File 2). In stool, the inflammatory markers neutrophil gelatinase-associated lipocalin (NGAL) and calprotectin and the gut permeability marker zonulin were significantly correlated (Fig. 2). In general, calprotectin and zonulin showed fewer correlations with other stool markers in the PD group as compared to the control group. Except for correlations with stool CXCL8 and IL-1β in the control group, none of the three most abundant SCFAs was significantly correlated with any of the markers in stool or plasma (Fig. 2).
Contrasting clinical variables against SCFAs and inflammatory markers suggested inverse correlations between stool SCFAs and several PD-related clinical variables such as the Non-Motor Symptoms Scale (NMSS) total score, the Rome III constipation/defecation subscore, and the Geriatric Depression Scale-15 (GDS15) (Fig. 3A, Additional File 4A). A similar pattern was seen for the Rome III constipation/defecation subscore, Rome III irritable bowel syndrome (IBS) criteria, and NMSS score and plasma CXCL8 (Fig. 3B, Additional File 4A). These correlations were significant for the entire cohort but not when evaluated only within PD patients. Butyric acid levels were correlated, however, with the age of onset for motor and non-motor symptoms of PD, while CXCL8 and IL-1β levels in stool were inversely correlated with the age of motor symptom onset (Fig. 3A, 3C, Additional File 4B). In PD patients, calprotectin levels were inversely linked to IBS symptoms, and stool zonulin was inversely correlated with variables related to PD severity, including the Unified Parkinson’s Disease Rating Scale (UPDRS) score, GDS15, and NMSS score (Fig. 3C, Additional File 4B). Symptom severity-related correlations remained significant after correcting for disease duration (Additional File 2).
Microbial alpha diversity is inversely associated with stool inflammatory and permeability markers and SCFAs
Microbial alpha diversity indices, which indicate richness and evenness of bacterial taxa, were inversely correlated with total SCFAs, butyric and propionic acid, as well as with NGAL, calprotectin, zonulin, and IL-1β in stool (Fig. 4). Correlations between alpha diversity and SCFAs, NGAL, and zonulin were most evident in control subjects while correlations with calprotectin were most apparent in PD patients. In patients only, a significant inverse correlation between the Shannon diversity index and stool CXCL8 was also found (Fig. 4). Linear regression modeling corrected for sex and PD/control status suggested that the relationship between alpha diversity and stool inflammatory markers was independent of sex except for calprotectin, for which an inverse relationship was observed in females but not in males (Fig. 5A, Additional File 5A-C). There was a significant interaction between PD/control status and alpha diversity for propionic acid, showing an inverse correlation between microbial alpha diversity and propionic acid concentration in controls but not in PD patients (Fig. 5B). The total SCFAs variable showed a similar trend (Additional File 5B). The interaction in the propionic acid model remained significant after confounder correction (Additional File 5D).
Microbial beta diversity is related to SCFA levels and inflammatory and permeability markers in stool
Microbial beta diversity describes the degree to which different communities of bacteria differ from one another in composition. We found significant associations between beta diversity and total SCFAs, acetic, butyric, and propionic acid as well as stool NGAL, zonulin, IL2, and the stool marker PC in full data, with or without correction for PD/control status, sex, constipation, and BMI (Additional File 6A). The effect of calprotectin, while not statistically significant, was worth noting (p = 0.058 for full data, p = 0.079 when corrected for confounders). Control-only comparisons showed significant effects for SCFAs, NGAL, and zonulin; only butyric acid and zonulin were significant in PD-only comparisons (Additional File 6A, 7). These two variables also had a significant interaction with PD/control status when tested with the variable split into two categories by median (Additional File 6B).
The Prevotella enterotype is associated with higher levels of butyric acid and lower levels of NGAL and zonulin in stool
One broad characterization method for microbiota composition is defining enterotypes based on certain indicator species in bacterial communities (36). In our subjects, the Firmicutes enterotype was more common and the Prevotella enterotype less common among PD patients compared to control subjects (Additional File 8A). Stool zonulin, NGAL, propionic acid, and butyric acid had differences in concentrations between enterotypes (Additional File 8B). Butyric acid levels were higher in PD patients with the Prevotella enterotype than the other two enterotypes, and control subjects with this enterotype had the lowest concentrations of NGAL and zonulin (Fig. 6). The findings regarding NGAL (full data and controls) and zonulin (controls only) remained significant after FDR correction (Additional File 8B).
Specific bacterial taxa are associated with levels of stool SCFAs and inflammatory and permeability markers
Guided by the associations with beta diversity, we explored the associations of specific bacterial taxa with the three most abundant SCFAs and with stool NGAL, calprotectin, and zonulin. The abundances of Butyricicoccus, Clostridium sensu stricto, and Roseburia were positively correlated with levels of SCFAs while the abundances of Akkermansia, Escherichia/Shigella, Flavonifractor, Intestinimonas, Phascolarctobacterium, and Sporobacter decreased with increasing SCFA concentrations (Fig. 7, Additional File 9, 10A). Similar patterns could be seen for the parent families of these taxa, such as a positive association between butyric acid and Lachnospiraceae and a negative one for Enterobacteriaceae and Verrucomicrobiaceae (Additional File 10A). A positive association between SCFAs and Bacteroides was observed only in controls, while a negative association with Bifidobacterium was observed only in PD patients (Fig. 7, Additional File 9, 10A).
Microbiota associations with the stool markers differed between PD patients and control subjects. Stool NGAL levels were positively associated with Bifidobacterium in PD but with Bacteroides and Roseburia in controls (Fig. 7, Additional File 10B, 11A). Relative abundance of Prevotella was negatively associated with both NGAL and zonulin in controls only (Fig. 7, Additional File 10B, 11). The only microbiota-immune association found to be significant in both PD patients and controls was the negative relationship between stool zonulin and the genus Coprococcus (Fig. 7, Additional File 10B, 11B), a taxon known to have a strong influence on gut permeability (41) and stool zonulin levels (42). Other associations with low-abundance taxa were statistically significant but were driven by a small number of data points, precluding confident interpretation (Additional File 11).