Next Article in Journal
Characteristics of Clinical Symptoms, Cerebral Images and Stroke Etiology in Vertebro-Basilar Artery Fenestration-Related Infarction
Next Article in Special Issue
Mini-Review on the Possible Interconnections between the Gut-Brain Axis and the Infertility-Related Neuropsychiatric Comorbidities
Previous Article in Journal
Role of Potassium Ions Quantum Tunneling in the Pathophysiology of Phantom Limb Pain
Previous Article in Special Issue
Probiotics Alleviate the Progressive Deterioration of Motor Functions in a Mouse Model of Parkinson’s Disease
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Can Gut Microbiota Be a Good Predictor for Parkinson’s Disease? A Machine Learning Approach

1
Department of Biology, University of Rome Tor Vergata, 00133 Rome, Italy
2
UOSD Parkinson’s Center, Department of Systems Medicine, University of Rome Tor Vergata, 00133 Rome, Italy
3
Department of Systems Medicine, University of Rome Tor Vergata, 00133 Rome, Italy
*
Author to whom correspondence should be addressed.
These authors have contributed equally to this work.
Brain Sci. 2020, 10(4), 242; https://doi.org/10.3390/brainsci10040242
Submission received: 16 March 2020 / Revised: 9 April 2020 / Accepted: 16 April 2020 / Published: 19 April 2020
(This article belongs to the Special Issue Brain–Microbiome Interactions)

Abstract

:
The involvement of the gut microbiota in Parkinson’s disease (PD), investigated in several studies, identified some common alterations of the microbial community, such as a decrease in Lachnospiraceae and an increase in Verrucomicrobiaceae families in PD patients. However, the results of other bacterial families are often contradictory. Machine learning is a promising tool for building predictive models for the classification of biological data, such as those produced in metagenomic studies. We tested three different machine learning algorithms (random forest, neural networks and support vector machines), analyzing 846 metagenomic samples (472 from PD patients and 374 from healthy controls), including our published data and those downloaded from public databases. Prediction performance was evaluated by the area under curve, accuracy, precision, recall and F-score metrics. The random forest algorithm provided the best results. Bacterial families were sorted according to their importance in the classification, and a subset of 22 families has been identified for the prediction of patient status. Although the results are promising, it is necessary to train the algorithm with a larger number of samples in order to increase the accuracy of the procedure.

1. Introduction

Neurodegenerative diseases represent a heterogeneous class of neurological disorders, with a high social, sanitary and economic impact. Parkinson’s disease (PD) is one of the most common neurodegenerative diseases, with a worldwide prevalence of 0.4% that is likely to double in the next 20 years [1]. The PD phenotype is characterized by movement disorders as a result of the loss of dopaminergic neurons in the substantia nigra caused by α-synuclein (α-syn) aggregates [2]. Only 10% of PD cases are due to genetic causes [3], indicating that environmental factors like dietary habits, head injury and nicotine consumption could trigger or influence the progression of PD [4]. Among environmental factors, the role of gut microbiota and its interactions with the gut–brain axis aroused the interest of researchers worldwide [5]. Gut microbiota can interact with the vagus nerve using neuroimmune and neuroendocrine mechanisms and, at the same time, the nervous system can modulate the gut physiology and environment, affecting the gut microbiota composition [6]. PD is influenced by this bidirectional communication, and PD patients show significant comorbidity with small-intestine bacterial overgrowth, constipation and Irritable Bowel Disease (IBD) like symptoms [7]. Furthermore, α-syn aggregates are found in the enteric nervous system and can spread to the central nervous system through the vagus nerve, and gut microbiota can affect the aggregation of α-syn [8].
The role of gut microbiota was evaluated in PD patients in several studies, mainly using targeted metagenomics and sequencing the 16S rRNA gene from fecal samples [2,9,10,11,12,13,14,15,16,17]. Some common features have been found, such as the reduction of bacteria producing short-chain fatty acids from the Lachnospiraceae family [5]. Regarding other bacterial families, the results are contradictory; for example, the Bifidobacteriaceae family was reported to have a lower abundance in PD patients in some studies [2,18] and higher abundances in others [11,19]. So, although all studies indicate a correlation between microbiota and Parkinson’s disease, there is no convergence as of yet on the bacterial families identifiable as specific biomarkers. To overcome this problem, we analyzed all published data with a computational procedure that can identify taxa involved in the pathology. Machine learning (ML) algorithms are appropriate tools to create predictive models that can distinguish the pathological status of a patient using the frequencies of bacteria in the feces.
ML algorithms are currently used for building predictive models for the classification of biological data, and identify biomarkers through a training procedure [20,21]. This technology was applied to identify marker genes in breast cancer [22], and to analyze clinical data for predicting cardiovascular and diabetes risk [23,24]. Recently, ML algorithms have been used to identify biomarkers analyzing shotgun and 16S rRNA data [25].
In this study, we use three different supervised ML algorithms to analyze 16S rRNA gene sequencing data derived from six different studies [2,9,10,11,12,13] downloaded from the Sequence Read Archive database. From this analysis, we obtain a classifier that can predict the pathological status of PD patients against healthy controls (HCs), and we identify a subset of 22 bacterial families that are discriminative for the prediction.

2. Materials and Methods

2.1. Datasets Downloadand Bioinformatic and Statistical Data Analysis

Datasets were downloaded from the Sequencing Read Archive (SRA) database [26] using the SRA Toolkit (http://ncbi.github.io/sra-tools/). Raw reads from the following BioProject IDs were selected: PRJNA510730 [9], PRJNA268515 [2], PRJEB14674 [13], PRJEB14928 [11], PRJNA381395 [12] and PRJEB27564 [10]. Non-fecal samples from PRJNA268515 and PRJNA381395 studies were removed before the analysis. The quality of raw sequencing reads was assessed with FastqC (https://www.bioinformatics.babraham.ac.uk/projects/fastqc/), and the bioinformatic data analysis was performed using the QIIME 2 pipeline [27]. Reads were quality filtered, chimera-checked and clustered in amplicon sequencing variants (ASVs) using DADA2 in the QIIME 2 pipeline [28]. The taxonomy of representative sequences was assessed using the q2-feature-classifier and the GreenGenes database [29,30]. Data normalization and filtering were performed using R 3.5.3. The dataset was loaded in R using the phyloseq package (version 1.26.1) [31]. Clusters (ASVs) with a number of reads lower than 0.05% of the total read number were removed from the analysis [32]. The read number count was normalized using DESeq2 [33]. The normalized table was summed up at the Family level using the tax_glom function in phyloseq.

2.2. Machine Learning Data Analysis

ML data analysis was performed using custom scripts in Python 3.6.7 language, using the sci-kit learn package (https://scikit-learn.org/stable/). Bacterial families were normalized across all samples using the StandardScaler method, which scales the distribution by subtracting the mean from each value and dividing the difference by the standard deviation. We used this type of standardization since, in the microbiota, some bacterial families are more abundant than others and their abundances can widely vary across samples.
The performance was evaluated using a stratified cross-validation (CV) with a K value equal to 5. The dataset was divided into training and test sets 5 times, allowing the training and the testing of the algorithms in 5 different iterations (folds). Each time, a different portion of the dataset was used in the training or the test set, in order to predict all the samples. In this way, we reduced the bias of the random division of samples in the training or test set.
We compared the performances of three different machine learning algorithms: (1) a random forest classifier (RF) [34], with 2000 estimators (trees) and a depth of 66 nodes (leaves); (2) a neural network (NN) classifier [35], with 3 hidden layers and respectively 180, 90 and 45 neurons for each layer; and (3) a support vector machine (SVM) classifier [36], with a polynomial kernel and cost parameter equal to 1.
The performance was evaluated using true positive cases (TP), false positive cases (FP), true negative cases (TN) and false negative cases (FN). A TP is a PD patient correctly classified as PD patient, while an FP is a PD patient incorrectly classified as HC. Conversely, a TN is an HC correctly classified as HC and an FP is an HC incorrectly classified as PD patient.
For each algorithm, at each fold, the number of TP, FP, TN and FN cases was computed and used to summarize the following metrics: accuracy, precision, recall, F-score and area under the curve (AUC). The accuracy is defined as the ratio of correctly predicted observations; the recall (or true positive rate, TPR) is the proportion of actual positives that are correctly identified as positive; the precision is the ratio of positive classifications identified correctly; and the F-score is the harmonic mean between precision and recall. The AUC was computed using the receiver operating characteristic (ROC) curve. The ROC curve summarizes the true positive rate and the false positive rate, and the AUC indicates the ability of the classifier to distinguish between two classes (i.e., PD or HC).
The importance of each bacterial family in the RF algorithm was evaluated using the “embedded feature selection strategy’’ analyzing the “Gini impurity decrease” [21]. The bacterial families were sorted from the most to the least relevant. The RF algorithm was then re-trained systematically using the first ‘’n’’ bacterial families, starting from the first n = 5 families and increasing this number until the AUC, the precision, the recall and the F-score values were comparable to those obtained with the whole family set (n = 52).

3. Results

3.1. Datasets Description and Supervised Machine Learning Approach

A total of 873 16S rRNA gene sequencing data of fecal samples from PD patients and HCs were downloaded from the Sequence Read Archive (SRA) database to study the association between the microbiota dysbiosis and diagnosis of PD. The distribution of PD and HC samples and the methodological approaches are reported in Table 1.
After the bioinformatic analysis, 846 samples were retained. PD patients (472 samples) represent 56% of the whole dataset. We uniformly processed the metagenomic data of all samples, using supervised ML algorithms. In detail, random forest (RF), neural network (NN) and support vector machine (SVM) were used for this evaluation since they are state-of-the-art approaches and are appropriate for this type of data [37]. The dataset was analyzed using cross-validation, randomly selecting 80% of the samples to create the training set, and the remaining 20% to create the test set and evaluate the prediction. The model was evaluated by resampling the test and training set 5 times, using a stratified cross-validation (K-fold = 5).
For each approach, the algorithm parameters were tuned through a grid search and were selected to optimize the training phase. Prediction performance was evaluated by the AUC metrics, which summarize true-positive and false-positive rates. The comparison of the three methods is reported in Figure 1A,B. The AUC is significantly higher for RF (0.80 ± 0.01) than for NN (0.67 ± 0.03) and SVM (0.54 ± 0.08) (Figure 1), indicating that RF is the most effective algorithm in distinguishing between status (PD or HC), according to the bacterial families’ frequencies in the feces. In line with this, Figure 1B indicates that the accuracy, precision, recall and F-score are higher for RF than for NN and SVM.

3.2. Feature Selection

Overall, the RF algorithm showed the best performance and has been selected to rank the importance of the bacterial families for the prediction. This process, defined as “feature selection”, identifies the most informative and relevant features in the classification of the PD status. The ranking of the bacterial families in discriminating between HC and PD patients is reported in Table 2.
In order to identify the minimal number of bacterial families that can reliably predict the pathological status, an embedded feature selection strategy was performed by re-training the RF algorithm with a subset of families and comparing the corresponding metrics. Initially, the subset included only the first 5 families in the ranking, shown in Table 2. The number of families was systematically increased to re-train the algorithm until the AUC, the precision, the recall and the F-score values were close enough to those obtained using the 52 total number of families for the training. This process allowed for the removal of bacterial families not strictly involved in the gut dysbiosis of Parkinson’s disease. The results indicate that by training the model with the first 22 bacterial families, the AUC, precision, recall the F-score values were almost identical to those obtained using the total number of families (Table S1), suggesting that this is the minimal subset that can be considered to correctly predict the pathological status.
The importance of the first 22 families in discriminating between patients and controls is plotted in Figure 2, together with their relative abundance. Interestingly, the plot shows that the importance of each family is not directly correlated with the relative family abundance in the samples. In some cases, we detected relevant species with low prevalence but high discriminative potential between “healthy” and “diseased” subjects. For example, Verrucomicrobiaceae/Akkermansiaceae and Bifidobacteriaceae were highly discriminative, although with a low average relative abundance.
It is also worth noting that not all of the 22 families identified in this analysis are cited in the previous studies investigating the role of gut microbiota in PD dysbiosis. In fact, among the first ten families identified as the most important ones in the PD diagnosis, eight were already identified in the literature, but two families—namely, Veillonellaceae and Alcaligenaceae—have never been reported before (Table 2).

4. Discussion

In this study we present, for the first time, an ML data analysis on microbiota dysbiosis in PD patients. We considered six available datasets from the SRA database, obtained from experiments carried out in different laboratories (Table 1). We downloaded and re-analyzed the datasets, uniformly processing the data using the most up-to-date bioinformatic procedures.
We initially evaluated the efficiency of three ML algorithms (RF, SVM and NN) in identifying samples belonging to HC or PD patients comparing different metrics (AUC, accuracy, precision, recall and F-score). The RF algorithm exhibited the best results providing an AUC of 80% and accuracy of 71%. This result is satisfying since we are analyzing data from studies that differ for participants’ nationality and for several methodological aspects, such as DNA extraction kit, sample transport and conservation (Table 1).
It is worth noting that, by analyzing the microbiota of diseases directly located in the gut, higher AUC values have been reported [23]. On the other hand, AUC values around 80% are in line with studies on the prediction of pathologies not directly related to the gut, such as obesity using 16S rRNA metagenomic data [25], or type-2 diabetes using shotgun metagenomic data [23]. This might indicate that the AUC value has reached the limit for Parkinson’s disease, although we cannot exclude that this value could be improved by increasing the sample size.
The relatively low accuracy (71%) could be due to various methodological approaches used in different laboratories in collecting, storing and processing data. As critically pointed out in a recent review [5], methodological inconsistencies between gut microbiome case–control studies in PD might contribute to the heterogeneity of the results. A lack of unique experimental and bioinformatic protocols prevents a direct and straightforward comparison of the data. We emphasize the importance of defining unique standards to permit a reliable comparison.
We found a subset of 22 bacterial families that provide prediction metrics almost identical to those obtained when the RF algorithm was trained with the whole microbiota (52 families). Subsets made by a lower number of relevant bacterial families (i.e., 5, 10) did not provide similar results, indicating that the combination of fewer species is insufficient to characterize the microbiota associated with this disease (Supplementary Table S1). This finding indicates the presence of a complex interplay of numerous bacterial families involved in gut dysbiosis in Parkinson’s disease. Interestingly, the rank of importance of each bacterial family is not directly correlated to its relative abundance.
Not all families identified by the RF algorithm were reported in the literature (Table 2). Indeed, eight of the first ten families in the rank (Lachnospiraceae, Ruminococcaceae, Bacteroidaceae, Verrucomicrobiaceae/Akkermansiaceae, Rikenellaceae, Bifidobacteriaceae, Porphyromonadaceae, Enterobacteriaceae) are cited, whilst two families, Veillonellaceae and Alcaligenaceae—both higher in PD patients—have never been reported before. The identification of new bacterial families that may play an important role in predicting the PD status highlights the power of a prediction analysis based on ML algorithms.
Finally, we want to point out that the values of the predictive metrics are probably too low for an immediate application of the procedure for the purpose of Parkinson’s disease’s diagnostics; however, the ranking of importance of bacterial families involved in the disease may help in its diagnosis.

5. Conclusions

In this work, we processed 846 16S rRNA microbiota data coming from six different studies, applying an ML approach. The RF algorithm provided an AUC of 80% and accuracy of 71% and identified a subset of 22 families that can be used to discriminate between PD and HC.
Unfortunately, the data deposited in the public databases are only a small fraction of the data that has been published up to now. This has prevented us from analyzing a larger number of data and increasing the population of the training and test sets. It is necessary to train the RF algorithm with a higher number of samples in order to increase the accuracy of the model and provide more robust results on the association between the gut microbiota and PD. We propose that the scientific community should build a network to share all the data produced by different laboratories, permitting the development of a fully reliable tool for the diagnosis and prognosis of this disease.

Supplementary Materials

The following are available online at https://www.mdpi.com/2076-3425/10/4/242/s1. Table S1: Random forest performance with a reduced number of bacterial families.

Author Contributions

Conception and design of the study, D.P. and A.D.; Acquisition of data, V.U. and R.C.; Bioinformatic and machine learning data analysis, D.P. and A.T.; Analysis and interpretation of data, D.P., A.T., and A.D.; Writing the article, D.P. and A.T.; Revising the article, D.P., A.T., V.U., R.C., S.B., A.S. and A.D.; All authors have read and agreed to the published version of the manuscript.

Funding

This work was partially supported by “Progetto Regione Lazio NEUROMICRO” to Alessandro Desideri.

Conflicts of Interest

The authors declare no conflicts of interest.

References

  1. Grochowska, M.; Laskus, T.; Radkowski, M. Gut Microbiota in Neurological Disorders. Arch. Immunol. Ther. Exp. (Warsz.) 2019, 67, 375–383. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  2. Keshavarzian, A.; Green, S.J.; Engen, P.A.; Voigt, R.M.; Naqib, A.; Forsyth, C.B.; Mutlu, E.; Shannon, K.M. Colonic bacterial composition in Parkinson’s disease. Mov. Disord. 2015, 30, 1351–1360. [Google Scholar] [CrossRef] [PubMed]
  3. Haikal, C.; Chen, Q.; Li, J. Microbiome changes: An indicator of Parkinson’s disease? Transl. Neurodegener. 2019, 8, 38. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  4. Singh, Y.; El-Hadidi, M.; Admard, J.; Wassouf, Z.; Schulze-Hentrich, J.M.; Kohlhofer, U.; Quintanilla-Martinez, L.; Huson, D.; Riess, O.; Casadei, N. Enriched environmental conditions modify the gut microbiome composition and fecal markers of inflammation in parkinson’s disease. Front. Neurosci. 2019, 13, 1032. [Google Scholar] [CrossRef] [Green Version]
  5. Boertien, J.M.; Pereira, P.A.B.; Aho, V.T.E.; Scheperjans, F. Increasing Comparability and Utility of Gut Microbiome Studies in Parkinson’s Disease: A Systematic Review. J. Parkinsons Dis. 2019, 9, S297–S312. [Google Scholar] [CrossRef] [Green Version]
  6. Martin, C.R.; Osadchiy, V.; Kalani, A.; Mayer, E.A. The Brain-Gut-Microbiome Axis. Cell. Mol. Gastroenterol. Hepatol. 2018, 6, 133–148. [Google Scholar] [CrossRef] [Green Version]
  7. Fasano, A.; Bove, F.; Gabrielli, M.; Petracca, M.; Zocco, M.A.; Ragazzoni, E.; Barbaro, F.; Piano, C.; Fortuna, S.; Tortora, A.; et al. The role of small intestinal bacterial overgrowth in Parkinson’s disease. Mov. Disord. 2013, 28, 1241–1249. [Google Scholar] [CrossRef]
  8. Fitzgerald, E.; Murphy, S.; Martinson, H.A. Alpha-synuclein pathology and the role of the microbiota in Parkinson’s disease. Front. Neurosci. 2019, 13, 369. [Google Scholar] [CrossRef]
  9. Pietrucci, D.; Cerroni, R.; Unida, V.; Farcomeni, A.; Pierantozzi, M.; Mercuri, N.B.; Biocca, S.; Stefani, A.; Desideri, A. Dysbiosis of gut microbiota in a selected population of Parkinson’s patients. Park. Relat. Disord. 2019, 65, 124–130. [Google Scholar] [CrossRef]
  10. Aho, V.T.E.; Pereira, P.A.B.; Voutilainen, S.; Paulin, L.; Pekkonen, E.; Auvinen, P.; Scheperjans, F. Gut microbiota in Parkinson’s disease: Temporal stability and relations to disease progression. EBioMedicine 2019, 44, 691–707. [Google Scholar] [CrossRef] [Green Version]
  11. Hopfner, F.; Künstner, A.; Müller, S.H.; Künzel, S.; Zeuner, K.E.; Margraf, N.G.; Deuschl, G.; Baines, J.F.; Kuhlenbäumer, G. Gut microbiota in Parkinson disease in a northern German cohort. Brain Res. 2017, 1667, 41–45. [Google Scholar] [CrossRef] [PubMed]
  12. Heintz-Buschart, A.; Pandey, U.; Wicke, T.; Sixel-Döring, F.; Janzen, A.; Sittig-Wiegand, E.; Trenkwalder, C.; Oertel, W.H.; Mollenhauer, B.; Wilmes, P. The nasal and gut microbiome in Parkinson’s disease and idiopathic rapid eye movement sleep behavior disorder. Mov. Disord. 2018, 33, 88–98. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  13. Hill-Burns, E.M.; Debelius, J.W.; Morton, J.T.; Wissemann, W.T.; Lewis, M.R.; Wallen, Z.D.; Peddada, S.D.; Factor, S.A.; Molho, E.; Zabetian, C.P.; et al. Parkinson’s disease and Parkinson’s disease medications have distinct signatures of the gut microbiome. Mov. Disord. 2017, 32, 739–749. [Google Scholar] [CrossRef] [PubMed]
  14. Li, F.; Wang, P.; Chen, Z.; Sui, X.; Xie, X.; Zhang, J. Alteration of the fecal microbiota in North-Eastern Han Chinese population with sporadic Parkinson’s disease. Neurosci. Lett. 2019, 707, 134297. [Google Scholar] [CrossRef]
  15. Li, W.; Wu, X.; Hu, X.; Wang, T.; Liang, S.; Duan, Y.; Jin, F.; Qin, B. Structural changes of gut microbiota in Parkinson’s disease and its correlation with clinical features. Sci. China Life Sci. 2017, 60, 1223–1233. [Google Scholar] [CrossRef]
  16. Petrov, V.A.; Saltykova, I.V.; Zhukova, I.A.; Alifirova, V.M.; Zhukova, N.G.; Dorofeeva, Y.B.; Tyakht, A.V.; Kovarsky, B.A.; Alekseev, D.G.; Kostryukova, E.S.; et al. Analysis of gut microbiota in patients with parkinson’s disease. Bull. Exp. Biol. Med. 2017, 162, 734–737. [Google Scholar] [CrossRef]
  17. Barichella, M.; Severgnini, M.; Cilia, R.; Cassani, E.; Bolliri, C.; Caronni, S.; Ferri, V.; Cancello, R.; Ceccarani, C.; Faierman, S.; et al. Unraveling gut microbiota in Parkinson’s disease and atypical parkinsonism. Mov. Disord. 2019, 34, 396–405. [Google Scholar] [CrossRef]
  18. Bedarf, J.R.; Hildebrand, F.; Coelho, L.P.; Sunagawa, S.; Bahram, M.; Goeser, F.; Bork, P.; Wüllner, U. Functional implications of microbial and viral gut metagenome changes in early stage L-DOPA-naïve Parkinson’s disease patients. Genome Med. 2017, 9, 39–52. [Google Scholar] [CrossRef]
  19. Lin, A.; Zheng, W.; He, Y.; Tang, W.; Wei, X.; He, R.; Huang, W.; Su, Y.; Huang, Y.; Zhou, H.; et al. Gut microbiota in patients with Parkinson’s disease in southern China. Park. Relat. Disord. 2018, 53, 82–88. [Google Scholar] [CrossRef]
  20. Libbrecht, M.W.; Noble, W.S. Machine learning applications in genetics and genomics. Nat. Rev. Genet. 2015, 16, 321–332. [Google Scholar] [CrossRef] [Green Version]
  21. Pasolli, E.; Truong, D.T.; Malik, F.; Waldron, L.; Segata, N. Machine Learning Meta-analysis of Large Metagenomic Datasets: Tools and Biological Insights. PLoS Comput. Biol. 2016, 12, e1004977. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  22. Tabl, A.A.; Alkhateeb, A.; ElMaraghy, W.; Rueda, L.; Ngom, A. A machine learning approach for identifying gene biomarkers guiding the treatment of breast cancer. Front. Genet. 2019, 10, 256. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  23. Kavakiotis, I.; Tsave, O.; Salifoglou, A.; Maglaveras, N.; Vlahavas, I.; Chouvarda, I. Machine Learning and Data Mining Methods in Diabetes Research. Comput. Struct. Biotechnol. J. 2017, 15, 104–116. [Google Scholar] [CrossRef] [PubMed]
  24. Weng, S.F.; Reps, J.; Kai, J.; Garibaldi, J.M.; Qureshi, N. Can Machine-learning improve cardiovascular risk prediction using routine clinical data? PLoS ONE 2017, 12, e0174944. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  25. Walters, W.A.; Xu, Z.; Knight, R. Meta-analyses of human gut microbes associated with obesity and IBD. FEBS Lett. 2014, 588, 4223–4233. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  26. Leinonen, R.; Sugawara, H.; Shumway, M. The sequence read archive. Nucleic Acids Res. 2011, 39, D19–D21. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  27. Bolyen, E.; Rideout, J.R.; Dillon, M.R.; Bokulich, N.A.; Abnet, C.C.; Al-Ghalith, G.A.; Alexander, H.; Alm, E.J.; Arumugam, M.; Asnicar, F.; et al. Reproducible, interactive, scalable and extensible microbiome data science using QIIME 2. Nat. Biotechnol. 2019, 37, 852–857. [Google Scholar] [CrossRef]
  28. Callahan, B.J.; McMurdie, P.J.; Rosen, M.J.; Han, A.W.; Johnson, A.J.A.; Holmes, S.P. DADA2: High-resolution sample inference from Illumina amplicon data. Nat. Methods 2016, 13, 581–583. [Google Scholar] [CrossRef] [Green Version]
  29. McDonald, D.; Price, M.N.; Goodrich, J.; Nawrocki, E.P.; Desantis, T.Z.; Probst, A.; Andersen, G.L.; Knight, R.; Hugenholtz, P. An improved Greengenes taxonomy with explicit ranks for ecological and evolutionary analyses of bacteria and archaea. ISME J. 2012, 6, 610–618. [Google Scholar] [CrossRef]
  30. Bokulich, N.A.; Kaehler, B.D.; Rideout, J.R.; Dillon, M.; Bolyen, E.; Knight, R.; Huttley, G.A.; Gregory Caporaso, J. Optimizing taxonomic classification of marker-gene amplicon sequences with QIIME 2’s q2-feature-classifier plugin. Microbiome 2018, 6, 90. [Google Scholar] [CrossRef]
  31. McMurdie, P.J.; Holmes, S. Phyloseq: An R Package for Reproducible Interactive Analysis and Graphics of Microbiome Census Data. PLoS ONE 2013, 8, e61217. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  32. Bokulich, N.A.; Subramanian, S.; Faith, J.J.; Gevers, D.; Gordon, J.I.; Knight, R.; Mills, D.A.; Caporaso, J.G. Quality-filtering vastly improves diversity estimates from Illumina amplicon sequencing. Nat. Methods 2013, 10, 57–59. [Google Scholar] [CrossRef] [PubMed]
  33. McMurdie, P.J.; Holmes, S. Waste not, want not: Why rarefying microbiome data is inadmissible. PLoS Comput. Biol. 2014, 10, e1003531. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  34. Chen, X.; Ishwaran, H. Random forests for genomic data analysis. Genomics 2012, 99, 323–329. [Google Scholar] [CrossRef] [Green Version]
  35. Samborska, I.A.; Alexandrov, V.; Sieczko, L.; Kornatowska, B.; Goltsev, V.; Cetner, M.D.; Kalaji, H.M. Artificial neural networks and their application in biological and agricultural research. J. NanoPhotoBioSci. 2014, 2, 14–30. [Google Scholar]
  36. Yang, Z.R. Biological applications of support vector machines. Brief. Bioinform. 2004, 5, 328–338. [Google Scholar] [CrossRef] [Green Version]
  37. Statnikov, A.; Henaff, M.; Narendra, V.; Konganti, K.; Li, Z.; Yang, L.; Pei, Z.; Blaser, M.J.; Aliferis, C.F.; Alekseyenko, A.V. A comprehensive evaluation of multicategory classification methods for microbiomic data. Microbiome 2013, 1, 11. [Google Scholar] [CrossRef]
Figure 1. (A) Average ROC curves (over 5 folds) with confidence intervals and (B) prediction performance metrics with the relative margin of error. The results for random forests are reported in green, for neural networks in purple and for support vector machines in orange, respectively.
Figure 1. (A) Average ROC curves (over 5 folds) with confidence intervals and (B) prediction performance metrics with the relative margin of error. The results for random forests are reported in green, for neural networks in purple and for support vector machines in orange, respectively.
Brainsci 10 00242 g001
Figure 2. List of the 22 bacterial families required for discriminating between HC and PD patients. For each family, the average percentage of abundance is represented by a bar, orange for HC and blue for PD patients (left scale). The importance of the family in discriminating the status is represented by a red dot (right scale).
Figure 2. List of the 22 bacterial families required for discriminating between HC and PD patients. For each family, the average percentage of abundance is represented by a bar, orange for HC and blue for PD patients (left scale). The importance of the family in discriminating the status is represented by a red dot (right scale).
Brainsci 10 00242 g002
Table 1. List of references, number of samples, methodological approaches and nationality of studies considered in this analysis.
Table 1. List of references, number of samples, methodological approaches and nationality of studies considered in this analysis.
ReferencePD SamplesHC SamplesSample TransportDNA Extraction16S RegionNationality
Method
[2]3431BD GaspakFastDNA Spin Kit for SoilV4United States
[10]6568NRPSP Spin Stool KitV3-V4Finland
[9]11682Stabilizer PSPPSP Spin Stool KitV3-V4Italy
[13]206133Ambient tempEarth microbiome project protocolV4United States
[12]2234Stabilizer PSPPSP Spin Stool KitV3-V4Germany
[11]2926Immediate freezingCustom Protocol HopfnerV4Russia
Table 2. Ranking of the importance of the bacterial families in discriminating between healthy controls and Parkinson’s disease (PD) patients.
Table 2. Ranking of the importance of the bacterial families in discriminating between healthy controls and Parkinson’s disease (PD) patients.
Bacterial FamilyRanking of ImportanceHigher (−) or Lower (+) Abundance in PD Patients from RF AlgorithmReferences in the Literature Reporting Overabundance in PD PatientsReferences in the Literature Reporting Lower Abundance in PD Patients
Lachnospiraceae1[14][2,9,10,11,13,15,16,17]
Ruminococcaceae2[2,14][13,15]
Bacteroidaceae3[2][16]
Verrucomicrobiaceae4+[2,12,13,17]
Rikenellaceae5+[19]
Bifidobacteriaceae6+[13,16,17,19][2,18]
Porphyromonadaceae7+[17]
Veillonellaceae8+
Enterobacteriaceae9+[9,15,17]
Alcaligenaceae10
Streptococcaceae11+[19]
Christensenellaceae12+[13,16,17]
Erysipelotrichaceae13+[19]
[Odoribacteraceae]14+
Prevotellaceae15[10,16,17]
Desulfovibrionaceae16+[14,19]
Coriobacteriaceae17[17]
Clostridiaceae18[2,10,16][19]
[Barnesiellaceae]19+
Lactobacillaceae20+[9,10,11,16]
[Tissierellaceae]21+
Peptostreptococcaceae22[16,17][14,15]
Methanobacteriaceae23+[19]
[Mogibacteriaceae]24
[Paraprevotellaceae]25+
Turicibacteraceae26
Pseudomonadaceae27+
Victivallaceae28
Campylobacteraceae29+
Synergistaceae30+
Pasteurellaceae31[14] [10,13,19]
Corynebacteriaceae32+
S24-733
Enterococcaceae34+[9,11,15,17]
Actinomycetaceae35+
Moraxellaceae36
Burkholderiaceae37
Comamonadaceae38+
Alcanivoracaceae39
Oxalobacteraceae40
Propionibacteriaceae41
Xanthomonadaceae42
Rhodobacteraceae43
Fusobacteriaceae44+
Staphylococcaceae45
Caulobacteraceae46+
Caldicoprobacteraceae47
Succinivibrionaceae48+
Peptococcaceae49
Flavobacteriaceae50+
[Weeksellaceae]51+
Aeromonadaceae52+

Share and Cite

MDPI and ACS Style

Pietrucci, D.; Teofani, A.; Unida, V.; Cerroni, R.; Biocca, S.; Stefani, A.; Desideri, A. Can Gut Microbiota Be a Good Predictor for Parkinson’s Disease? A Machine Learning Approach. Brain Sci. 2020, 10, 242. https://doi.org/10.3390/brainsci10040242

AMA Style

Pietrucci D, Teofani A, Unida V, Cerroni R, Biocca S, Stefani A, Desideri A. Can Gut Microbiota Be a Good Predictor for Parkinson’s Disease? A Machine Learning Approach. Brain Sciences. 2020; 10(4):242. https://doi.org/10.3390/brainsci10040242

Chicago/Turabian Style

Pietrucci, Daniele, Adelaide Teofani, Valeria Unida, Rocco Cerroni, Silvia Biocca, Alessandro Stefani, and Alessandro Desideri. 2020. "Can Gut Microbiota Be a Good Predictor for Parkinson’s Disease? A Machine Learning Approach" Brain Sciences 10, no. 4: 242. https://doi.org/10.3390/brainsci10040242

APA Style

Pietrucci, D., Teofani, A., Unida, V., Cerroni, R., Biocca, S., Stefani, A., & Desideri, A. (2020). Can Gut Microbiota Be a Good Predictor for Parkinson’s Disease? A Machine Learning Approach. Brain Sciences, 10(4), 242. https://doi.org/10.3390/brainsci10040242

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop