- Open Access
Functional implications of microbial and viral gut metagenome changes in early stage L-DOPA-naïve Parkinson’s disease patients
- J. R. Bedarf†1, 2,
- F. Hildebrand†4,
- L. P. Coelho4,
- S. Sunagawa4, 5,
- M. Bahram9, 10,
- F. Goeser3, 11,
- P. Bork4, 6, 7, 8, 13Email author and
- U. Wüllner1, 2, 12Email author
© The Author(s). 2017
Received: 15 November 2016
Accepted: 8 April 2017
Published: 28 April 2017
The Erratum to this article has been published in Genome Medicine 2017 9:61
Parkinson’s disease (PD) presently is conceptualized as a protein aggregation disease in which pathology involves both the enteric and the central nervous system, possibly spreading from one to another via the vagus nerves. As gastrointestinal dysfunction often precedes or parallels motor symptoms, the enteric system with its vast diversity of microorganisms may be involved in PD pathogenesis. Alterations in the enteric microbial taxonomic level of L-DOPA-naïve PD patients might also serve as a biomarker.
We performed metagenomic shotgun analyses and compared the fecal microbiomes of 31 early stage, L-DOPA-naïve PD patients to 28 age-matched controls.
We found increased Verrucomicrobiaceae (Akkermansia muciniphila) and unclassified Firmicutes, whereas Prevotellaceae (Prevotella copri) and Erysipelotrichaceae (Eubacterium biforme) were markedly lowered in PD samples. The observed differences could reliably separate PD from control with a ROC-AUC of 0.84. Functional analyses of the metagenomes revealed differences in microbiota metabolism in PD involving the ẞ-glucuronate and tryptophan metabolism. While the abundances of prophages and plasmids did not differ between PD and controls, total virus abundance was decreased in PD participants. Based on our analyses, the intake of either a MAO inhibitor, amantadine, or a dopamine agonist (which in summary relates to 90% of PD patients) had no overall influence on taxa abundance or microbial functions.
Our data revealed differences of colonic microbiota and of microbiota metabolism between PD patients and controls at an unprecedented detail not achievable through 16S sequencing. The findings point to a yet unappreciated aspect of PD, possibly involving the intestinal barrier function and immune function in PD patients. The influence of the parkinsonian medication should be further investigated in the future in larger cohorts.
Idiopathic Parkinson’s disease (PD) disease is conceptualized as a progressive protein aggregation disease with the formation of neuronal cytoplasmic aggregates of misfolded α-synuclein (α-syn) and other proteins as the neuropathological hallmark (Lewy bodies [LB]) . LB are present not only in the central nervous system (CNS) but also in the enteric nervous system (ENS) of the entire gastrointestinal tract, corresponding to the clinical notion that the gastrointestinal tract is involved in PD .
Lewy, in his original thesis work in 1913, already identified the dorsal motor nucleus of the vagus as a hotspot of brain pathology and Braak et al. more recently confirmed the early involvement of the vagus and hypothesized that PD might originate in the gut and that α-syn aggregation might spread via vagal structures into the CNS and higher cortical regions . The importance of the vagus for partly bi-directional interactions of the ENS and the CNS (“gut-brain axis”) [4–6] has become even more evident since the identification of the cholinergic anti-inflammatory pathway . In line with this direct connection of gut and brain neurons, increasing evidence from cell culture and animal experiments seems to support the hypothesis of spreading or seeding of α-syn oligomers [8, 9]. The deposition of α-syn within the colon can help to distinguish PD patients from controls [10–13]. Nevertheless, when assessed with conventional immunohistochemistry, its diagnostic value as a biomarker has not been finally confirmed as α-syn staining in colonic mucosa was likewise found in PD patients and controls , which not necessarily refutes the concept of an intestinal origin of PD pathogenesis. Thus, additional yet unidentified factors beyond α-syn must be involved in the presumed PD disease process.
The gastrointestinal microbiome encompasses a vast diversity of bacterial species and may be regarded as an extracorporeal organ system, which interacts with its host in unprecedented ways just being unraveled now . Severely disturbed gut homeostasis is detrimental for the host but the effects of smaller changes or differences in species variation for nutrition [16, 17], behavior , and drug metabolism  are just beginning to emerge. Recent studies already linked an altered microbiome to PD, but most participants were well advanced and treated with L-DOPA [20–22], which affects colonic motility and may promote intestinal bacterial overgrowth .
Using metagenomic shotgun analysis, we found that early, L-DOPA-naïve PD patients carry an altered gut microbiota composition, i.e. specific taxonomic groups, among others related to intestinal barrier and immune functions, were overrepresented or underrepresented. Functional analyses also suggested differences in microbiota ẞ-glucuronate and tryptophan degrading pathways. Moreover, total virus abundance was decreased in PD participants.
Study participants and clinical characteristics
This study was approved by the local ethics committee of the University of Bonn and all participants gave informed consent (internal ethics vote 126/02). Study participants were recruited from the Department of Neurology at the University of Bonn. To reduce any potential gender effects, we included only male participants in the study. Thirty-one male PD patients (diagnosed according to the UK Brain Bank criteria ) were compared to 28 male age-matched non-parkinsonian controls.
Disease severity was measured using the Unified PD Disease Rating Scale (UPDRS part III). Gastrointestinal symptoms and presence of constipation were assessed with a modified version of an interview-based Gastrointestinal Symptom Rating Scale (GSRS) (selected items: borborygmus, abdominal distension, increased flatus, decreased passage of stools, increased passage of stools, loose stools, hard stools, urgent need for defecation and feeling of incomplete evacuation, each item was rated 0–3 according to intensity, frequency, duration or social impact depending on the respective item, Additional file 1) . To avoid alterations of gut microbiota related to either gastrointestinal dysfunction of late stage PD or L-DOPA-induced intestinal effects, we included only early stage PD participants (onset of motor symptoms and diagnosis of PD within the past year) who were naïve to L-DOPA therapy. Further exclusion criteria were: (1) chronic and inflammatory gastrointestinal diseases including chronic constipation; (2) the use of laxatives or immunosuppressive agents in the past three months; (3) atypical or secondary parkinsonism; while (4) the use of antibiotics in the past three months in principle was an exclusion criterion; however, we included three PD patients and three controls despite the intake of antibiotics for one to three days in a period of 28–34 days prior to feces sampling as the omission of those cases from the analyses did not change any result (see also Additional file 2).
Clinical characteristics and general demographic parameters of study participants
Age (years, mean ± SD)
64.8 ± 9.5
65.6 ± 10.4
UPDRS III (mean ± SD)
12.6 ± 6.9
0 ± 0
GIT symptoms incl. constipation (GSRS, mean ± SD)
3.4 ± 2.9
2.2 ± 2.0
Total serum bilirubin
0.23 ± 0.03
0.23 ± 0.03
Analysis of microbiota
Non-invasively obtained stool samples were shotgun sequenced (paired end) using an Illumina Hiseq4000 and further analyzed with the MOCAT2 pipeline . Briefly, taxonomic mapping quality-controlled generated profiles reads (minimum length cutoff 45 bp, minimum quality score cutoff 20, reads matching Illumina adapters, or human genome were removed) to a database of ten universal single-copy marker genes that were extracted from 3496 NCBI reference genomes and 263 human gut metagenomes . For abundance estimates at the species level, we used mOTU (molecular operational taxonomic units) abundances  and for genus-level and family-level abundances, individual NCBI taxonomy-annotated marker gene abundances were summed up.
Quality-controlled reads were also mapped against an annotated reference gene catalog for the human gut microbiome. On average, 88% of all reads could be mapped to the reference gene catalogue (Additional file 3). Abundances of individual genes summed up according to their KO annotation in the KEGG (Kyoto Encyclopedia of Genes and Genomes) database. KEGG and GMM (Gut specific Metabolic Modules) pathway  abundances were estimated with the same algorithm as in , re-implemented in C++, available from github.com/hildebra/Rarefaction. Briefly, we estimated for each metabolic pathway which of several alternatives had the highest coverage given our KO abundance matrix. If the coverage per pathway was higher than 30%, the median abundance of all KOs in this pathway was used to estimate pathway abundance.
To estimate the number of phages and viruses, we choose a reference gene catalog independent approach, as assembly of mobile genetic elements is notoriously error prone and imprecise. MOCAT2 quality filtered reads were mapped against the ACLAME  database using Diamond  in sensitive mode, after format conversion from fastq to fasta using sdm . Reads mapping with an e-value <1e-7 were considered valid hits against the database and the number of reads mapping to the classes of “Plasmid,” “Prophage,” and “Virus” were counted and normalized by read number.
Statistical analysis was conducted in R 3.0.0. For all univariate tests, taxa with less than five reads over all samples were excluded from this analysis to avoid artifacts, similar to the approach in . The matrix was normalized dividing each feature by the respective total sample sum and transformed with log10(x + 1), where x is the normalized feature coverage as calculated in the mOTUs algorithm. For species and mOTU level analysis, we made the filtering options more explicit to exclude spurious correlations: from the abundance matrix features were removed that were absent in more than ten samples, had less than ten accumulated or two mean read coverage. Sample count matrices were rarefied using the R implementation of the rtk toolkit . Significance between groups of samples was tested with a Kruskal–Wallis test, as implemented in R. The ordinations (non-metric multidimensional [NMDS] ordination) and subsequent statistical analysis were calculated using the R-package vegan with Bray-Curtis distance on the rarefied and log-transformed taxa abundance and visualized with custom R scripts. Intergroup differences for the microbiota were calculated using a PERMANOVA test as implemented in vegan . This test compares the intragroup distances to the intergroup distances in a permutation scheme and thus calculates a P value. For all PERMANOVA tests we used 4999 randomizations. PERMANOVA post hoc P values were corrected for multiple testing using the Benjamini–Hochberg false discovery rate (q-value) . Sample composition plots were visualized with custom R scripts (available from github.com/hildebra/PD_helpers_R).
Univariate testing for differential abundances of each taxonomic unit between two or more groups was tested using a Kruskal–Wallis test (P value), corrected for multiple testing using the Benjamini–Hochberg false discovery rate (q-value).
The Blocked “independence test” function calls (with the following options: “ytrafo = rank, teststat = scalar” for blocked WRST) were used from the COIN R package  to control for potential confounders, such as the intake of a statin.
Post hoc statistical testing for significant differences between all combinations of two groups was conducted only for taxa with a significance of P < 0.2. Wilcoxon rank-sum tests were calculated for all possible group combinations and corrected for multiple testing using Benjamini–Hochberg false discovery rate (q-value). GSRS, UPDRS III, and Bilirubin correlations to taxa were tested using a spearman correlation test; P values were corrected using Benjamini–Hochberg false discovery rate.
Further, we confirmed these results using ANCOM , a statistical test developed for microbial count data, using an R implementation version 1.1-3 and the additional parameters multcorr = 2 and sig = 0.1, that is with multiple testing correction at significance 0.1.
For generating a classifier, genera were filtered, removing any whose mean relative abundance across samples was below 0.1%. Subsequently, relative abundances were normalized by a log-transformation with a pseudo-count equal to one-tenth of the estimated detection limit, estimated as the minimal abundance of any positively detected taxa . Fitting a classifier proceeds in two steps. First, a lasso-penalized logistic regression classifier is used to select the top features (the features with the highest absolute weight are selected, Additional file 4). Second, these features are used in an unpenalized logistic regression classifier. Classifiers were based on the scikit-learn implementation . The classifier was generated using cross-validation in a leave-one-out schedule (feature selection was performed de novo at each iteration to provide unbiased estimates). Relevant taxa were extracted from a model trained on the whole data. P values were computed by the Mann–Whitney test on the prediction scores of the two classes (PD and control). Classifier performance was reported as ROC-AUC, representing the probability that the classifier will correctly label a new sample.
Structural equation model (SEM) from the same sample distribution
SEMs were used to test causality through examination of both direct and indirect effects of bacteria and medication on PD and vice versa . We also evaluated an alternative model in which medication affects bacteria and bacteria in turn affect the disease. Based on this we used an exploratory approach to approximate a significant model fit. Briefly, non-significant variables and paths were subsequently removed using backward elimination stepwise regression and new paths were added based on modification indices until significant model fit was achieved (Additional file 5). This analysis was performed using the computer program AMOS ver. 7.0 (SPSS, Chicago, IL, USA).
Microbiota in PD
Analyzing the two sample sets for differences in their bacterial composition we detected significant differences between PD and control samples (Permanova test, P < 0.001, Additional file 6).
Key taxonomic differences in the PD gut microbiota
In general, we observed Clostridiales, Bacteroidaceae, and Ruminococcaceae as the most abundant families in both sample sets accounting for 50–52% of the relative family abundance. Univariate tests revealed significant changes in family as well as in genus abundance in PD samples (Additional file 8). Verrucomicrobiaceae (genus Akkermansia), unclassified Bacteria (of the classified prokaryotes) and Firmicutes were increased, whereas Prevotellaceae (genus Prevotella) and Erysipelotrichaceae (genus Eubacterium) were markedly lowered in PD samples (Fig. 1b, Additional file 9a, and b; P < 0.05, q < 0.1; with the exception of Prevotella (genus) with a P < 0.05 and a q = 0.13, thus being n. s.). Using ANCOM, we similarly found genera Akkermansia, Prevotella, Eubacterium, unknown Bacteria, and unknown Firmicutes to be significantly different between PD and control patients.
Metagenomics provides species resolution. Among the classified bacteria, we thus identified a pattern of significantly increased key species in PD including Akkermansia muciniphila and Alistipes shahii. On the other hand, Prevotella copri, Eubacterium biforme, and Clostridium saccharolyticum were decreased (Fig. 1b and c).
Correlation of microbiota with PD clinical scores
No significant taxonomic associations were detected, neither at genus nor at species level, when microbiota abundance was correlated with clinical data (UPDRS III, GSRS, or total serum bilirubin, despite of the presumed decreased ẞ-glucuronidation in PD participants, Additional file 10). This was expected, as we had aimed to recruit a cohort of PD patients with a short duration of (motor) symptoms and subjective impairment; nevertheless, some interesting trends (q > 0.1) were observed regarding the symptom severity of PD (UPDRS III) for three different Eubacteria (E. eligens, E. rectale, and E. hallii, Fig. 1d).
Both parkinsonian and gastrointestinal symptoms severity including constipation were rather low in our sample set of early L-DOPA-naïve patients and, except for core parkinsonian features, did not differ from controls (Table 1). We could not detect any deviations pointing to a possible confounder.
Functional analysis of the PD microbiome
In order to explore differences in the metabolic potential of gut microbiota between PD and control patients, we further estimated the abundance of metabolic pathways, using our metagenomic reads mapped to functional orthologues from the KEGG and GMM databases (Additional file 11).
Interestingly, we also found two different pathways for tryptophan metabolism, which appeared to be more active in the microbiota of PD participants (tryptophan → kynurenine → 2-aminomuconate, KEGG module number M00038, Fig. 2a and b; tryptophan degradation, MF0009 (GMM module number), Fig. 2a and b, increased formate conversion, GMM module number MF0118, Fig. 2a and b). Although a lack of tryptophan and serotonin is a hallmark of PD, this finding did not reach statistical significance after multiple testing correction (all p < 0.05, q > 0.1) and needs further exploration in a larger cohort.
To determine which bacteria are involved in these pathways, we traced the contributing genes and determined their likely taxonomic origin (Fig. 2a and b). Although multiple genera contribute to all modules, only in modules MF0065 and MF0118 could we find evidence of any genera contributing significantly more reads in PD patients than control patients, after multiple testing correction (Fig. 2c). In both cases this was Eubacterium.
Influence of medication on the microbiota
Given the important influence of pharmaceuticals on the gut microbiota , we tested the effects of the concomitant PD-specific medication in our cohort (Additional file 12). No significant differences in taxa abundances (families, genera, species) were apparent in the relatively small samples grouped according to the various combinations of anti-PD medication (dopamine agonist + MAO inhibitor + amantadine: n = 10, MAO inhibitor + amantadine: n = 14, MAO inhibitor: n = 3, no therapy: n = 3; there was only one patient with dopamine agonist + MAO inhibitor, who was not included). However, patients treated with MAO inhibitors in combination with amantadine (n = 14, i.e. 45.2% of the PD participants) displayed a significantly increased richness (Additional file 12d), which did not affect the overall richness of the entire PD cohort.
Among the concomitant drugs the intake of a statin showed an influence on the gut microbiota. We identified five families significantly different between statin-treated and untreated patients (Burkholderiaceae, Propionibacteriaceae, Enterococcaceae, Actinomycetaceae, and Enterobacteriaceae), none of which contributed to the differences observed between PD participants and controls (Additional files 12a, c and 13).
Virus abundance is lowered in PD participants
We tested the fraction of the reads that could be mapped to the ACLAME database to estimate the amount of mobile elements in the metagenomes, though this method is limited to known diversity of mobile elements. While the abundances for prophages and plasmids were not different between PD and controls, total virus abundance (reflecting bacterial and archaeal phages) was decreased in PD participants (P = 6.7e-5, Additional file 14a). We also found that viruses were increased in participants treated with a statin (P = 0.0002), and therefore tested for differences between PD participants and controls, corrected for statin treatment. This again showed a decrease of the virus load in PD participants (P = 0.0009, Additional file 14b).
Gut microbes discriminate PD participants from controls
Structural equation model
Our results show a significantly altered microbial composition in early disease stage L-DOPA-naïve PD participants. Compared to earlier studies investigating microbiota in PD patients using 16S-based techniques the chosen methods allowed to detect changes at the species level and also a decreased virus load in PD.
We confirmed a decrease of Prevotella copri in PD and in addition found decreased Eubacterium biforme and Clostridium saccharolyticum and increased Akkermansia muciniphila as well as Alistipes shahii. Based on the taxonomic differences alone, logistic regression after feature selection allowed separation of PD patients in early stages from controls with good accuracy (AUC 0.84). Furthermore, our analyses point to differences in microbiota metabolism, namely in ẞ-glucuronate and tryptophan degrading pathways.
Despite several differences in study design compared to previous work [20, 21], our results strengthen the hypothesis of a PD-specific “microbial footprint.” The observed differences though are not directly comparable, since the previous studies were using 16S amplicon sequencing with its various biases , while we used metagenomics with increased precision due to usage of single copy marker genes with high confidence taxonomic assignments . Testing for causality inference of biotic and abiotic factors seem to indicate that these taxonomic changes are a consequence of the disease. However, whether the observed changes are primary changes or rather secondary, resulting from unidentified effects, and whether it is either beneficial or harmful cannot be decided at present.
The increase in Akkermansia muciniphila, which appears to be in certain consistency in regards to the findings of Keshavarzian et al., Scheperjans et al., and Unger et al. [20–22], may illustrate this dilemma: it is a common mucin degrader which has been shown to reverse diet induced pathological intestinal changes in high-fat fed mice by restoring the intestinal mucus layer and the underlying epithelium, thus being able to improve gut barrier function [44, 45]. Protective effects of extracellular vesicles derived from Akkermansia muciniphila on experimentally induced colitis further support a beneficial influence on intestinal immunity . On the other hand, inflammatory and regulatory properties have been reported for Akkermansia, probably mediated due to an increased exposure of immune cells to microbial antigens upon breaking down the mucosal mucin layer . Preliminary evidence also linked Akkermansia to multiple sclerosis . Previous studies on colonic biopsies and feces samples from treated and drug-naïve PD participants suggested an altered mucosal barrier function and PD patients exhibited significantly greater intestinal permeability than controls, paralleled by an increased mucosal staining for E. coli and α-syn . Pro-inflammatory dysbiosis may even trigger α-syn misfolding or neuronal injury from gut-derived endotoxins [21, 50]. Although we did not identify E. coli species associated with PD in our samples, the increase in Akkermansia might be associated with a yet unexplored disease related impact on mucosal barrier function.
Extending the results of Scheperjans et al. and Unger et al. [20, 22], which pointed to a relatively lower abundance of Prevotellaceae in advanced PD, Prevotella copri was markedly lowered in our samples of early stage PD. On the other side, Keshavarzian et al.  did not show this difference for fecal samples, albeit the trend was the same for mucosal-derived PD samples. However, Prevotella abundance was also reduced in Japanese multiple sclerosis patients and in autistic children, somewhat questioning the specificity of this finding [51, 52]. The Prevotella enterotype is the least prevalent in human individuals  and is related to dietary/fiber intake [54, 55]. In particular, Prevotella enrichment has been linked to non-Western and/or fiber-rich diets [56, 57]. Fibers are the primary substrate for short chain fatty acids (SCFAs) including butyrate and reductions in the latter can disrupt barrier function and promote inflammation . The fact that, in various autoimmune diseases including type 1/2 diabetes, irritable bowel disease, rheumatoid arthritis, and Behcet’s disease reduced levels of Prevotella have been found [42, 59–61], could indicate a decreased SCFA production (i.e. propionate) and in turn favor inflammatory conditions in PD .
In contrast to Scheperjans et al. , we did not find increased Ruminococcaceae (phylum Firmicutes) to compensate lower levels of Prevotella but instead an increase in unclassified Firmicutes. Interestingly, although the early, L-DOPA-naïve PD patients hold a different species pattern not yet affected by drug effects and the chronic constipation typically observed in late-stage PD, we found a certain consistency with the advanced PD patients’ pattern.
Taken together, the observed bacterial pattern in our PD samples might hint towards yet unexplored mechanisms of a disturbed intestinal and immune function in PD pathogenesis. Colonic biopsies from PD patients indeed showed enhanced pro-inflammatory cytokines and glial markers correlating with disease progression [62, 63]. Furthermore, there is evidence of α-syn contributing to neuro-inflammation by potentiating microglial or astroglial activation . In line with this, recent work highlighted the crucial role of microbiota on maturation and activity of microglia , which have been considered as one of the earliest contributors of neurodegeneration  and further supports the importance of microbial-derived mediators (gut peptides, chemokines, SCFAs) on immune regulation and CNS function [65, 67].
A possible role for SCFAs in PD
SCFAs are essential energy sources for colonocytes and reduced levels of SCFAs might not only contribute to a decreased colonic motility (i.e. constipation) but also led to an increase in intestinal barrier leakiness [68–70]. Keshavarzian et al. and Unger et al. both suggested a beneficial role for SCFAs as PD-derived feces were shown to contain less SCFA butyrate-producing bacteria, including Blautia, Roseburia, and Coprococcus  as well as Faecalibacterium prausnitzii [21, 22], which were previously attributed to exert putative anti-inflammatory effects.
While SCFA administration contra-intuitively promoted motor dysfunction and α-syn-mediated neuro-inflammation in a germ-free transgenic mouse model over-expressing α-syn, oral administration of heat killed bacteria had no effect on motor performance, indicating the putative importance for metabolically active microbiota in disease pathogenesis . Namely, when PD-derived microbiota (of treatment-naïve new onset PD donors) were orally transferred to germ-free mice, several taxa, including Roseburia, Rikenellaceae, and Enterococcus, were markedly altered in the microbial profile of the recipient mice independently of its genotype as when they received microbiota derived from healthy donors.
While inconclusive at the moment, prospective research on SFCA gene expression and metabolomic profiles of microbiota in health and disease will shed further light on this aspect.
The gut bacteria harbor a diverse phageome and virome that may contribute to function and structure of the microbiome, but evidence from comparative analyses of the human gut phageome is limited. Recently a comprehensive metagenomic analysis in 64 individuals suggested a core phageome that was shared among more than one-half of all individuals and might also exert beneficial properties as it was reduced in individuals with inflammatory bowel disease . However, based on our analysis, we could not find any differences in the abundance of prophages and plasmids between PD and control samples. In contrast, total virus amount was significantly lowered in PD participants.
Importantly, the assessment of virus and phage load is entirely dependent on the corresponding protein families being present in the ACLAME database; therefore, we only detect those phages or viruses, of which a closely related reference genome is present in the database. Testing for correlations between bacterial family and viral load abundance showed no significant correlations after multiple testing.
As viruses interact with host cells and influence immune response (i.e. prevent inflammatory conditions ), there might be various possibilities in which viruses interact in the pathogenesis of PD. Although inconclusive at the moment, exploration of the specific role of viruses in PD is a promising avenue to follow-up with more specific research.
Accumulating evidence suggests a direct impact of metabolic alterations in microbiota on human health [74, 75]. We observed a putative reduction in microbiota ẞ-glucuronidase activity in early stage PD participants. Decreased ẞ-glucuronidation in the microbiota could imply a deterioration of resistance to various pathogenic organisms . Also, microbial derived ẞ-glucuronidases affect effective dose availability of administered drugs by reactivation in the gut, which has been shown for irinotecan therapy in colorectal cancer patients . Altered metabolisms of xenobiotics or parkinsonian pharmaceuticals metabolized in the ẞ-glucuronate degrading pathway must be determined experimentally .
Our data further revealed a trend towards an increased tryptophan degradation gene copy number in PD. If one assumes that this increased genetic potential translates into an increased tryptophan metabolization, this finding is in line with previous research and is of particular interest as L-tryptophan, the precursor for serotonin, is decreased in PD patients’ brains. L-tryptophan is also metabolized to kynurenines, whereof metabolites have regulatory immune function and were described as either harmful or beneficial in PD [79–81]. Urinary metabolomics profiling demonstrated significant changes of urinary markers including an increased tryptophan metabolism, which was associated with the progression of PD . Interestingly, catabolism of serotonin also includes glucuronidation in the human intestine .
The association with these metabolic pathways point to a deeper involvement of Eubacteria with PD. Indeed Eubacteria spp. were decreased in PD (Eubacterium biforme) and other Eubacteria species (E. hallii, E. rectale, E. eligens) showed a trend towards correlation with disease severity (n.s.). Specifically colonizing the mucus layer, particularly Eubacterium rectale, might be interconnected with processes directly affecting the mucus layer due to its ability to gain access via flagella . Further, diversity of Eubacterium rectale was also reduced in an in vitro dynamic gut model (M-SHIME) of long-term colonization of the mucin layer when microbiota were derived from ulcerative colitis patients . Additionally, Eubacterium halii is viewed as a key species impacting the microbial balance due to its ability to produce several SCFAs . In turn, alterations in the abundance of different Eubaceria might contribute to the PD pathogenesis via metabolic but also direct mucosal pathways.
Lowered Eubacteria (family Erysipelotrichaceae) in mucosal as well as in fecal PD samples were similarly observed in the study of Keshavarzian et al. ; however, a correlation with disease severity was not proven.
Instead, Keshavarzian et al. found PD duration correlating with the greatest number of taxa, whereby the family Lachnospiraceae, which includes several (supposedly anti-inflammatory) butyrate producing bacteria, displayed a significant negative correlation. Scheperjans et al. further showed a significant association of Enterobacteriaceae with the postural instability and gait disorder (PIGD) PD phenotype , which was not confirmed in the work of Unger et al. .
Namely, based on our analyses, the intake of different anti-parkinsonian drugs had no overall influence on taxa abundance or microbial functions. However, in future subgroup analyses PD patients under the therapy with MAO-inhibitors and amantadine might be favorably influenced by an increased richness if assessed in a lager cohort. In this context, it is worth noting another study, which was published during the revision process of this manuscript and which demonstrates instead independent effects of different PD medications on the microbiome .
However, although the intake of a statin showed an influence on the gut microbiota with, in total, five families being different in statin-treated individuals, none of them contributed to the differences observed between PD participants and controls when controlling for statin intake with differential statistical methods. One caveat of testing for confounders in our cohort is that this result might be limited by the sample size being too small to find even small effects, which is an unavoidable inherent aspect of human cohort studies. However, future studies should address this aspect.
Our data revealed differences of colonic microbiota between PD patients and controls at an unprecedented detail not achievable through 16S sequencing: altered representation of several taxa including Eubacterium biforme, which has not been reported previously and might be limited to detection via metagenomics. The functional differences in the gut microbiota included ẞ-glucuronate and tryptophan degrading pathways. The findings point to a yet-unappreciated aspect of PD, possibly involving the intestinal barrier function and immune function in PD patients. We further show the benefits of integrating functional microbiota predictions into microbial-based profiles to discriminate health and disease that is promising as it holds the potential to identify PD patients. Furthermore, it is now evident that exploration of the PD virus populations is a promising avenue to follow up with more specific research.
This study was supported by the PD Fonds Deutschland gGmbH.
We thank members of the Bork group at the European Molecular Biology Laboratory (EMBL) for helpful discussions. We acknowledge the EMBL Genomics Core Facility for sequencing support and Y. P. Yuan and the EMBL Information Technology Core Facility for support with high-performance computing.
JRB received honoraria for scientific presentations and travel grants from IPSEN Pharma Germany and Merz pharmaceuticals Germany and is supported by the PD Fonds Deutschland gGmbH. FH received funding through the European Union’s Horizon 2020 research and innovation program under the Marie Skłodowska-Curie grant agreement no. 660375. FG is partly funded by grants of the German Center for Infection Research (DZIF). This study and LPC were partly supported by the European Molecular Biology Laboratory (EMBL). SS is supported by the Helmut Horten Foundation. MB was funded by Helge Ax:son Johnsons Stiftelse and PUT1317. PB has received funding from the European Research Council (ERC) Microbios (ERC-AdG-669830) for this project.
UW served as consultant and lecturer and on advisory boards for Boehringer-Ingelheim, Glaxo-SmithKline, Novartis and Pfizer Pharma GmbH, Lifescience KG, Lundbeck AG, Teva Pharma, UCB Pharma and Zambon and received grant/research funding from the NationalAtaxiaFoundation, Stiftung Doppelfeld, Stiftung Verum, MerckKgaA, Pharmacia&Upjohn, the Deutsche Parkinson Vereinigung (dPV), the Hans Tauber Stiftung, the International Parkinson Fonds Deutschland gGmbH, the German Research Council (DFG), the German Ministry for Research (BMBF) and the German Centre for Neurodegenerative Diseases (DZNE). He has received support from the EU/EFPIA Innovative Medicines Initiative Joint Undertaking Aetionomy [grant number 115568] for this work.
Availability of data and materials
The datasets supporting the conclusions of this article are included within the article and its additional files. All metagenomic sequences have been deposited in the European Bioinformatics Institute-Sequence Read Archive database, under accession number ERP019674.
JRB participated in the study design and patients’ recruitment, analyzed and interpreted clinical data, interpreted bioinformatics data, and drafted the manuscript. FH analyzed and interpreted all bioinformatics data, designed statistical tests to test hypothesis, and drafted the manuscript. LPC trained and evaluated the LASSO classifier. SS analyzed raw sequencing data involving MOCAT. MB performed the mathematical modeling regarding the SEM model. FG participated in the study design. PB participated in the study design and interpreted bioinformatics data. UW participated in the study design and patients’ recruitment, and interpreted clinical and bioinformatics data. All authors critically reviewed the article and approved the final version.
The authors declare that they have no competing interests.
Consent for publication
Ethics approval and consent to participate
This study has been performed in accordance with the Declaration of Helsinki and was approved by the local ethics committee of the University of Bonn and all participants gave written informed consent (126/02).
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
Open AccessThis article is distributed under the terms of the Creative Commons Attribution 4.0 International License (http://creativecommons.org/licenses/by/4.0/), which permits unrestricted use, distribution, and reproduction in any medium, provided 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 Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated.
- Luna E, Luk KC. Bent out of shape: alpha-Synuclein misfolding and the convergence of pathogenic pathways in Parkinson’s disease. FEBS Lett. 2015;589(24 Pt A):3749–59.View ArticlePubMedPubMed CentralGoogle Scholar
- Wakabayashi K, Takahashi H, Ohama E, Ikuta F. Parkinson’s disease: an immunohistochemical study of Lewy body-containing neurons in the enteric nervous system. Acta Neuropathol. 1990;79(6):581–3.View ArticlePubMedGoogle Scholar
- Braak H, de Vos RA, Bohl J, Del Tredici K. Gastric alpha-synuclein immunoreactive inclusions in Meissner’s and Auerbach’s plexuses in cases staged for Parkinson’s disease-related brain pathology. Neurosci Lett. 2006;396(1):67–72.View ArticlePubMedGoogle Scholar
- Collins SM, Surette M, Bercik P. The interplay between the intestinal microbiota and the brain. Nat Rev Microbiol. 2012;10(11):735–42.View ArticlePubMedGoogle Scholar
- Grenham S, Clarke G, Cryan JF, Dinan TG. Brain-gut-microbe communication in health and disease. Front Physiol. 2011;2:94.View ArticlePubMedPubMed CentralGoogle Scholar
- Furness JB. The enteric nervous system and neurogastroenterology. Nat Rev Gastroenterol Hepatol. 2012;9(5):286–94.View ArticlePubMedGoogle Scholar
- Tracey KJ. The inflammatory reflex. Nature. 2002;420(6917):853–9.View ArticlePubMedGoogle Scholar
- Luk KC, Kehm V, Carroll J, Zhang B, O’Brien P, Trojanowski JQ, et al. Pathological alpha-synuclein transmission initiates Parkinson-like neurodegeneration in nontransgenic mice. Science. 2012;338(6109):949–53.View ArticlePubMedPubMed CentralGoogle Scholar
- Uchihara T, Giasson BI. Propagation of alpha-synuclein pathology: hypotheses, discoveries, and yet unresolved questions from experimental and human brain studies. Acta Neuropathol. 2016;131(1):49–73.View ArticlePubMedGoogle Scholar
- Shannon KM, Keshavarzian A, Dodiya HB, Jakate S, Kordower JH. Is alpha-synuclein in the colon a biomarker for premotor Parkinson’s disease? Evidence from 3 cases. Mov Disord. 2012;27(6):716–9. doi:10.1002/mds.25020.View ArticlePubMedGoogle Scholar
- Shannon KM, Keshavarzian A, Mutlu E, Dodiya HB, Daian D, Jaglin JA, et al. Alpha-synuclein in colonic submucosa in early untreated Parkinson’s disease. Mov Disord. 2012;27(6):709–15. doi:10.1002/mds.23838.View ArticlePubMedGoogle Scholar
- Lebouvier T, Chaumette T, Damier P, Coron E, Touchefeu Y, Vrignaud S, et al. Pathological lesions in colonic biopsies during Parkinson’s disease. Gut. 2008;57(12):1741–3. doi:10.1136/gut.2008.162503.View ArticlePubMedGoogle Scholar
- Lebouvier T, Neunlist M, Bruley des Varannes S, Coron E, Drouard A, N’Guyen JM, et al. Colonic biopsies to assess the neuropathology of Parkinson's disease and its relationship with symptoms. PLoS One. 2010;5(9), e12728. doi:10.1371/journal.pone.0012728.View ArticlePubMedPubMed CentralGoogle Scholar
- Antunes L, Frasquilho S, Ostaszewski M, Weber J, Longhino L, Antony P, et al. Similar alpha-Synuclein staining in the colon mucosa in patients with Parkinson’s disease and controls. Mov Disord. 2016;31(10):1567–70.View ArticlePubMedGoogle Scholar
- Gilbert JA, Quinn RA, Debelius J, Xu ZZ, Morton J, Garg N, et al. Microbiome-wide association studies link dynamic microbial consortia to disease. Nature. 2016;535(7610):94–103.View ArticlePubMedGoogle Scholar
- Chung WS, Walker AW, Louis P, Parkhill J, Vermeiren J, Bosscher D, et al. Modulation of the human gut microbiota by dietary fibres occurs at the species level. BMC Biol. 2016;14:3.View ArticlePubMedPubMed CentralGoogle Scholar
- Sonnenburg JL, Backhed F. Diet-microbiota interactions as moderators of human metabolism. Nature. 2016;535(7610):56–64.View ArticlePubMedGoogle Scholar
- Dinan TG, Cryan JF. Microbes, immunity, and behavior: psychoneuroimmunology meets the microbiome. Neuropsychopharmacology. 2017;42(1):178–92.View ArticlePubMedGoogle Scholar
- Klaassen CD, Cui JY. Review: mechanisms of how the intestinal microbiota alters the effects of drugs and bile acids. Drug Metab Dispos. 2015;43(10):1505–21.View ArticlePubMedPubMed CentralGoogle Scholar
- Scheperjans F, Aho V, Pereira PA, Koskinen K, Paulin L, Pekkonen E, et al. Gut microbiota are related to Parkinson’s disease and clinical phenotype. Mov Disord. 2015;30(3):350–8.View ArticlePubMedGoogle Scholar
- Keshavarzian A, Green SJ, Engen PA, Voigt RM, Naqib A, Forsyth CB, et al. Colonic bacterial composition in Parkinson’s disease. Mov Disord. 2015;30(10):1351–60.View ArticlePubMedGoogle Scholar
- Unger MM, Spiegel J, Dillmann KU, Grundmann D, Philippeit H, Burmann J, et al. Short chain fatty acids and gut microbiota differ between patients with Parkinson’s disease and age-matched controls. Parkinsonism Relat Disord. 2016;32:66–72.View ArticlePubMedGoogle Scholar
- Fasano A, Visanji NP, Liu LW, Lang AE, Pfeiffer RF. Gastrointestinal dysfunction in Parkinson’s disease. Lancet Neurol. 2015;14(6):625–39.View ArticlePubMedGoogle Scholar
- Hughes AJ, Daniel SE, Kilford L, Lees AJ. Accuracy of clinical diagnosis of idiopathic Parkinson’s disease: a clinico-pathological study of 100 cases. J Neurol Neurosurg Psychiatry. 1992;55(3):181–4.View ArticlePubMedPubMed CentralGoogle Scholar
- Svedlund J, Sjodin I, Dotevall G. GSRS--a clinical rating scale for gastrointestinal symptoms in patients with irritable bowel syndrome and peptic ulcer disease. Dig Dis Sci. 1988;33(2):129–34.View ArticlePubMedGoogle Scholar
- Kultima JR, Coelho LP, Forslund K, Huerta-Cepas J, Li SS, Driessen M, et al. MOCAT2: a metagenomic assembly, annotation and profiling framework. Bioinformatics. 2016;32(16):2520–3.View ArticlePubMedPubMed CentralGoogle Scholar
- Sunagawa S, Mende DR, Zeller G, Izquierdo-Carrasco F, Berger SA, Kultima JR, et al. Metagenomic species profiling using universal phylogenetic marker genes. Nat Methods. 2013;10(12):1196–9.View ArticlePubMedGoogle Scholar
- Darzi Y, Falony G, Vieira-Silva S, Raes J. Towards biome-specific analysis of meta-omics data. ISME J. 2016;10(5):1025–8.View ArticlePubMedGoogle Scholar
- Le Chatelier E, Nielsen T, Qin J, Prifti E, Hildebrand F, Falony G, et al. Richness of human gut microbiome correlates with metabolic markers. Nature. 2013;500(7464):541–6.View ArticlePubMedGoogle Scholar
- Leplae R, Lima-Mendez G, Toussaint A. ACLAME: a CLAssification of Mobile genetic Elements, update 2010. Nucleic Acids Res. 2010;38(Database issue):D57–61.View ArticlePubMedGoogle Scholar
- Buchfink B, Xie C, Huson DH. Fast and sensitive protein alignment using DIAMOND. Nat Methods. 2015;12(1):59–60.View ArticlePubMedGoogle Scholar
- Hildebrand F, Tadeo R, Voigt AY, Bork P, Raes J. LotuS: an efficient and user-friendly OTU processing pipeline. Microbiome. 2014;2(1):30.View ArticlePubMedPubMed CentralGoogle Scholar
- Hildebrand F, Ebersbach T, Nielsen HB, Li X, Sonne SB, Bertalan M, et al. A comparative analysis of the intestinal metagenomes present in guinea pigs (Cavia porcellus) and humans (Homo sapiens). BMC Genomics. 2012;13:514. http://doi.org/10.1186/1471-2164-13-514.View ArticlePubMedPubMed CentralGoogle Scholar
- Saary P, Forslund K, Bork P, Hildebrand F. RTK: efficient rarefaction analysis of large datasets. Bioinformatics. 2017. doi:10.1093/bioinformatics/btx206.PubMedGoogle Scholar
- Anderson M. A new method for non-parametric multivariate analysis of variance. Austral Ecol. 2001;26(1):32–46.Google Scholar
- Benjamini Y, Hochberg Y. Controlling the false discovery rate: a practical and powerful approach to multiple testing. J R Stat Soc B. 1995;57(1):289–300.Google Scholar
- Hothorn T, Hornik K, van de Wiel MA, Zeileis A. A Lego system for conditional inference. Am Stat. 2006;60:257–63.View ArticleGoogle Scholar
- Mandal S, Van Treuren W, White RA, Eggesbø M, Knight R, Peddada SD. Analysis of composition of microbiomes: a novel method for studying microbial composition. Microb Ecol Health Dis. 2015;26:27663. doi:10.3402/mehd.v26.27663.PubMedGoogle Scholar
- Sunagawa S, Coelho LP, Chaffron S, Kultima JR, Labadie K, Salazar G, et al. Ocean plankton. Structure and function of the global ocean microbiome. Science. 2015;348(6237):1261359.View ArticlePubMedGoogle Scholar
- Pedregosa F, Varoquaux G, Gramfort A, Michel V, Thirion B, Grisel O. Scikit-learn: Machine learning in Python. J Mach Learn Res. 2011;12:2825–30.Google Scholar
- Grace J, Anderson TM, Olff H, Scheiner SM. On the specification of structural equation models for ecological systems. Ecol Monogr. 2010;80:67–8.View ArticleGoogle Scholar
- Forslund K, Hildebrand F, Nielsen T, Falony G, Le Chatelier E, Sunagawa S, et al. Disentangling type 2 diabetes and metformin treatment signatures in the human gut microbiota. Nature. 2015;528(7581):262–6.View ArticlePubMedPubMed CentralGoogle Scholar
- Tedersoo LAS, Bahram M, Põlme S, Riit T, Liiv I, Kõljalg U, et al. Shotgun metagenomes and multiple primer pair-barcode combinations of amplicons reveal biases in metabarcoding analyses of fungi. MycoKeys. 2015;10:1–43.View ArticleGoogle Scholar
- van Passel MW, Kant R, Zoetendal EG, Plugge CM, Derrien M, Malfatti SA, et al. The genome of Akkermansia muciniphila, a dedicated intestinal mucin degrader, and its use in exploring intestinal metagenomes. PLoS One. 2011;6(3), e16876.View ArticlePubMedPubMed CentralGoogle Scholar
- Everard A, Belzer C, Geurts L, Ouwerkerk JP, Druart C, Bindels LB, et al. Cross-talk between Akkermansia muciniphila and intestinal epithelium controls diet-induced obesity. Proc Natl Acad Sci U S A. 2013;110(22):9066–71.View ArticlePubMedPubMed CentralGoogle Scholar
- Kang CS, Ban M, Choi EJ, Moon HG, Jeon JS, Kim DK, et al. Extracellular vesicles derived from gut microbiota, especially Akkermansia muciniphila, protect the progression of dextran sulfate sodium-induced colitis. PLoS One. 2013;8(10), e76520.View ArticlePubMedPubMed CentralGoogle Scholar
- Ganesh BP, Klopfleisch R, Loh G, Blaut M. Commensal Akkermansia muciniphila exacerbates gut inflammation in Salmonella Typhimurium-infected gnotobiotic mice. PLoS One. 2013;8(9), e74963.View ArticlePubMedPubMed CentralGoogle Scholar
- Cantarel BL, Waubant E, Chehoud C, Kuczynski J, DeSantis TZ, Warrington J, et al. Gut microbiota in multiple sclerosis: possible influence of immunomodulators. J Investig Med. 2015;63(5):729–34.View ArticlePubMedPubMed CentralGoogle Scholar
- Forsyth CB, Shannon KM, Kordower JH, Voigt RM, Shaikh M, Jaglin JA, et al. Increased intestinal permeability correlates with sigmoid mucosa alpha-synuclein staining and endotoxin exposure markers in early Parkinson’s disease. PLoS One. 2011;6(12), e28032.View ArticlePubMedPubMed CentralGoogle Scholar
- Kelly LP, Carvey PM, Keshavarzian A, Shannon KM, Shaikh M, Bakay RA, et al. Progression of intestinal permeability changes and alpha-synuclein expression in a mouse model of Parkinson’s disease. Mov Disord. 2014;29(8):999–1009.View ArticlePubMedGoogle Scholar
- Kang DW, Park JG, Ilhan ZE, Wallstrom G, Labaer J, Adams JB, et al. Reduced incidence of Prevotella and other fermenters in intestinal microflora of autistic children. PLoS One. 2013;8(7), e68322.View ArticlePubMedPubMed CentralGoogle Scholar
- Miyake S, Kim S, Suda W, Oshima K, Nakamura M, Matsuoka T, et al. Dysbiosis in the gut microbiota of patients with multiple sclerosis, with a striking depletion of species belonging to Clostridia XIVa and IV clusters. PLoS One. 2015;10(9), e0137429.View ArticlePubMedPubMed CentralGoogle Scholar
- Arumugam M, Raes J, Pelletier E, Le Paslier D, Yamada T, Mende DR, et al. Enterotypes of the human gut microbiome. Nature. 2011;473(7346):174–80.View ArticlePubMedPubMed CentralGoogle Scholar
- Wu GD, Chen J, Hoffmann C, Bittinger K, Chen YY, Keilbaugh SA, et al. Linking long-term dietary patterns with gut microbial enterotypes. Science. 2011;334(6052):105–8.View ArticlePubMedPubMed CentralGoogle Scholar
- Yatsunenko T, Rey FE, Manary MJ, Trehan I, Dominguez-Bello MG, Contreras M, et al. Human gut microbiome viewed across age and geography. Nature. 2012;486(7402):222–7.PubMedPubMed CentralGoogle Scholar
- De Filippo C, Cavalieri D, Di Paola M, Ramazzotti M, Poullet JB, Massart S, et al. Impact of diet in shaping gut microbiota revealed by a comparative study in children from Europe and rural Africa. Proc Natl Acad Sci U S A. 2010;107(33):14691–6.View ArticlePubMedPubMed CentralGoogle Scholar
- David LA, Maurice CF, Carmody RN, Gootenberg DB, Button JE, Wolfe BE, et al. Diet rapidly and reproducibly alters the human gut microbiome. Nature. 2014;505(7484):559–63.View ArticlePubMedGoogle Scholar
- Kelly CJ, Zheng L, Campbell EL, Saeedi B, Scholz CC, Bayless AJ, et al. Crosstalk between microbiota-derived short-chain fatty acids and intestinal epithelial HIF augments tissue barrier function. Cell Host Microbe. 2015;17(5):662–71.View ArticlePubMedPubMed CentralGoogle Scholar
- Wang W, Chen L, Zhou R, Wang X, Song L, Huang S, et al. Increased proportions of Bifidobacterium and the Lactobacillus group and loss of butyrate-producing bacteria in inflammatory bowel disease. J Clin Microbiol. 2014;52(2):398–406.View ArticlePubMedPubMed CentralGoogle Scholar
- de Goffau MC, Fuentes S, van den Bogert B, Honkanen H, de Vos WM, Welling GW, et al. Aberrant gut microbiota composition at the onset of type 1 diabetes in young children. Diabetologia. 2014;57(8):1569–77.View ArticlePubMedGoogle Scholar
- Consolandi C, Turroni S, Emmi G, Severgnini M, Fiori J, Peano C, et al. Behcet’s syndrome patients exhibit specific microbiome signature. Autoimmun Rev. 2015;14(4):269–76.View ArticlePubMedGoogle Scholar
- Gao HM, Zhang F, Zhou H, Kam W, Wilson B, Hong JS. Neuroinflammation and alpha-synuclein dysfunction potentiate each other, driving chronic progression of neurodegeneration in a mouse model of Parkinson’s disease. Environ Health Perspect. 2011;119(6):807–14.View ArticlePubMedPubMed CentralGoogle Scholar
- Sanchez-Guajardo V, Tentillier N, Romero-Ramos M. The relation between alpha-synuclein and microglia in Parkinson’s disease: Recent developments. Neuroscience. 2015;302:47–58.View ArticlePubMedGoogle Scholar
- Devos D, Lebouvier T, Lardeux B, Biraud M, Rouaud T, Pouclet H, et al. Colonic inflammation in Parkinson's disease. Neurobiol Dis. 2013;50:42–8.View ArticlePubMedGoogle Scholar
- Erny D, Hrabe de Angelis AL, Jaitin D, Wieghofer P, Staszewski O, David E, et al. Host microbiota constantly control maturation and function of microglia in the CNS. Nat Neurosci. 2015;18(7):965–77.View ArticlePubMedGoogle Scholar
- Prinz M, Priller J. Microglia and brain macrophages in the molecular age: from origin to neuropsychiatric disease. Nat Rev Neurosci. 2014;15(5):300–12.View ArticlePubMedGoogle Scholar
- Borre YE, O’Keeffe GW, Clarke G, Stanton C, Dinan TG, Cryan JF. Microbiota and neurodevelopmental windows: implications for brain disorders. Trends Mol Med. 2014;20(9):509–18.View ArticlePubMedGoogle Scholar
- Ganapathy V, Thangaraju M, Prasad PD, Martin PM, Singh N. Transporters and receptors for short-chain fatty acids as the molecular link between colonic bacteria and the host. Curr Opin Pharmacol. 2013;13(6):869–74. doi:10.1016/j.coph.2013.08.006.View ArticlePubMedGoogle Scholar
- Singh N, Gurav A, Sivaprakasam S, Brady E, Padia R, Shi H, et al. Activation of Gpr109a, receptor for niacin and the commensal metabolite butyrate, suppresses colonic inflammation and carcinogenesis. Immunity. 2014;40(1):128–39. doi:10.1016/j.immuni.2013.12.007.View ArticlePubMedPubMed CentralGoogle Scholar
- Soret R, Chevalier J, De Coppet P, Poupeau G, Derkinderen P, Segain JP, et al. Short-chain fatty acids regulate the enteric neurons and control gastrointestinal motility in rats. Gastroenterology. 2010;138(5):1772–82. doi:10.1053/j.gastro.2010.01.053.View ArticlePubMedGoogle Scholar
- Sampson TR, Debelius JW, Thron T, Janssen S, Shastri GG, Ilhan ZE, et al. Gut microbiota regulate motor deficits and neuroinflammation in a model of Parkinson’s disease. Cell. 2016;167(6):1469–80. doi:10.1016/j.cell.2016.11.018.e12.View ArticlePubMedGoogle Scholar
- Manrique P, Bolduc B, Walk ST, van der Oost J, de Vos WM, Young MJ. Healthy human gut phageome. Proc Natl Acad Sci U S A. 2016;113(37):10400–5.View ArticlePubMedPubMed CentralGoogle Scholar
- Karst SM. Viral Safeguard: The enteric virome protects against gut inflammation. Immunity. 2016;44(4):715–8.View ArticlePubMedGoogle Scholar
- Puertollano E, Kolida S, Yaqoob P. Biological significance of short-chain fatty acid metabolism by the intestinal microbiome. Curr Opin Clin Nutr Metab Care. 2014;17(2):139–44.View ArticlePubMedGoogle Scholar
- Huda-Faujan N, Abdulamir AS, Fatimah AB, Anas OM, Shuhaimi M, Yazid AM, et al. The impact of the level of the intestinal short chain Fatty acids in inflammatory bowel disease patients versus healthy subjects. Open Biochem J. 2010;4:53–8.View ArticlePubMedPubMed CentralGoogle Scholar
- Allegretti JR, Kearney S, Li N, Bogart E, Bullock K, Gerber GK, et al. Recurrent Clostridium difficile infection associates with distinct bile acid and microbiome profiles. Aliment Pharmacol Ther. 2016;43(11):1142–53.View ArticlePubMedPubMed CentralGoogle Scholar
- Wallace BD, Wang H, Lane KT, Scott JE, Orans J, Koo JS, et al. Alleviating cancer drug toxicity by inhibiting a bacterial enzyme. Science. 2010;330(6005):831–5.View ArticlePubMedPubMed CentralGoogle Scholar
- Nayak RR, Turnbaugh PJ. Mirror, mirror on the wall: which microbiomes will help heal them all? BMC Med. 2016;14:72.View ArticlePubMedPubMed CentralGoogle Scholar
- Jenkins TA, Nguyen JC, Polglaze KE, Bertrand PP. Influence of tryptophan and serotonin on mood and cognition with a possible role of the gut-brain axis. Nutrients. 2016;8(1), E56.View ArticlePubMedGoogle Scholar
- Sanz Y, Walder K, Maes M. The role of the microbial metabolites including tryptophan catabolites and short chain fatty acids in the pathophysiology of immune-inflammatory and neuroimmune disease. Mol Neurobiol. 2016. doi:10.1007/s12035-016-0004-2.Google Scholar
- Lim CK, Essa MM, de Paula MR, Lovejoy DB, Bilgin AA, Waly MI, et al. Altered kynurenine pathway metabolism in autism: Implication for immune-induced glutamatergic activity. Autism Res. 2016;9(6):621–31.View ArticlePubMedGoogle Scholar
- Luan H, Liu LF, Tang Z, Zhang M, Chua KK, Song JX, et al. Comprehensive urinary metabolomic profiling and identification of potential noninvasive marker for idiopathic Parkinson’s disease. Sci Rep. 2015;5:13888.View ArticlePubMedPubMed CentralGoogle Scholar
- Sakakibara Y, Katoh M, Kawayanagi T, Nadai M. Species and tissue differences in serotonin glucuronidation. Xenobiotica. 2016;46(7):605–11.View ArticleGoogle Scholar
- Louis P, Flint HJ. Diversity, metabolism and microbial ecology of butyrate-producing bacteria from the human large intestine. FEMS Microbiol Lett. 2009;294(1):1–8.View ArticlePubMedGoogle Scholar
- Vermeiren J, Van den Abbeele P, Laukens D, Vigsnaes LK, De Vos M, Boon N, et al. Decreased colonization of fecal Clostridium coccoides/Eubacterium rectale species from ulcerative colitis patients in an in vitro dynamic gut model with mucin environment. FEMS Microbiol Ecol. 2012;79(3):685–96. doi:10.1111/j.1574-6941.2011.01252.x.View ArticlePubMedGoogle Scholar
- Engels C, Ruscheweyh HJ, Beerenwinkel N, Lacroix C, Schwab C. The common gut microbe Eubacterium hallii also contributes to intestinal propionate formation. Front Microbiol. 2016;7:713. doi:10.3389/fmicb.2016.00713.View ArticlePubMedPubMed CentralGoogle Scholar
- Hill-Burns EM, Debelius JW, Morton JT, Wissemann WT, Lewis MR, Wallen ZD, et al. Parkinson’s disease and Parkinson’s disease medications have distinct signatures of the gut microbiome. Mov Disord. 2017. doi:10.1002/mds.26942.PubMedGoogle Scholar