Table of Contents Author Guidelines Submit a Manuscript
BioMed Research International
Volume 2018 (2018), Article ID 2936257, 7 pages
https://doi.org/10.1155/2018/2936257
Research Article

Metagenomics Biomarkers Selected for Prediction of Three Different Diseases in Chinese Population

1Wuhan National Laboratory for Optoelectronics, Key Laboratory of Information Storage System, Huazhong University of Science and Technology, Wuhan, Hubei 430000, China
2Binhai Genomics Institute, BGI-Tianjin, BGI-Shenzhen, Tianjin 300308, China
3Tianjin Translational Genomics Center, BGI-Tianjin, BGI-Shenzhen, Tianjin 300308, China
4School of Mathematics and Computer Science, Guangdong Ocean University, Zhanjiang, Guangdong 524088, China
5BGI-Shenzhen, Shenzhen 518083, China
6School of Computer Science and Technology, Huazhong University of Science and Technology, Wuhan, Hubei 430000, China

Correspondence should be addressed to Ke Zhou; nc.ude.tsuh@uohz.k

Received 17 July 2017; Revised 14 October 2017; Accepted 24 October 2017; Published 11 January 2018

Academic Editor: Clara G. de los Reyes-Gavilan

Copyright © 2018 Honglong Wu et al. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

Abstract

The dysbiosis of human microbiome has been proven to be associated with the development of many human diseases. Metagenome sequencing emerges as a powerful tool to investigate the effects of microbiome on diseases. Identification of human gut microbiome markers associated with abnormal phenotypes may facilitate feature selection for multiclass classification. Compared with binary classifiers, multiclass classification models deploy more complex discriminative patterns. Here, we developed a pipeline to address the challenging characterization of multilabel samples. In this study, a total of 300 biomarkers were selected from the microbiome of 806 Chinese individuals (383 controls, 170 with type 2 diabetes, 130 with rheumatoid arthritis, and 123 with liver cirrhosis), and then logistic regression prediction algorithm was applied to those markers as the model intrinsic features. The estimated model produced an score of 0.9142, which was better than other popular classification methods, and an average receiver operating characteristic (ROC) of 0.9475 showed a significant correlation between these selected biomarkers from microbiome and corresponding phenotypes. The results from this study indicate that machine learning is a vital tool in data mining from microbiome in order to identify disease-related biomarkers, which may contribute to the application of microbiome-based precision medicine in the future.

1. Introduction

The human microbiome plays an important role in energy harvesting, metabolism of dietary components, immunity, inflammatory bowel disease, cancer therapy, and the progression of cancers [14]. Recently, lots of studies have shown that multiple factors, such as living environment, eating habits, age, gender, and the state of health, are associated with the homeostasis of human microbiome [5, 6]. For example, in the early age, most people share common functionality of the gut microbiome, but age-associated changes in vitamin biosynthesis and metabolism were identified from 326 individuals [7, 8]. High-throughput sequencing technologies have paved the way to investigate the human microbiome in a comprehensive manner. Individual microbiome may enhance the utility of precision medicine, personalized diagnostics, and treatment modalities [911].

Machine learning algorithms have been widely introduced to tackle problems in both genomics and bioinformatics, for instance, the identification and annotation of genomic regulatory elements like promoters, transcription start sites, enhancers, splice sites, and the classification of various phenotypes [12]. Some researchers have established large-scale microbiome datasets for type 2 diabetes, liver cirrhosis, rheumatoid arthritis, and colorectal carcinoma [1317]. Metagenome-wide association studies (MWAS) based on binary classification methods have been proven to be a powerful tool to analyze the correlation between diseases and microbiome [18, 19]. To assess the risk of developing certain diseases throughout one’s life, microbiome-based multiclass biomarkers may provide a promising noninvasive diagnostic tool in large-scale population. It was also reported that many factors, such as gender, age, and body mass index (BMI), will influence the diversity of human microbiome [5, 7]. Thus, these factors were also involved in this study.

Here, we performed multiclass classification analysis on the microbiome of 806 Chinese individuals, together with phenotype information, including gender, age, and BMI. Feature-associated gene markers were identified after addressing the imbalanced datasets between classes, which could be used to establish models to predict the type of disease. The human microbiome-based knowledge that associated with the abnormal phenotypes provides cues to the future application of precision medicine [11].

2. Materials and Methods

2.1. Sample Information

We process 806 shotgun metagenomic samples from three studies, 170 type 2 diabetes (T2D), 130 liver cirrhosis, 123 rheumatoid arthritis (RA), and 383 normal controls, to identify microbiome-based biomarkers. All these normal controls are with normal values on recent screen for physical examination and other clinical testing [1517]. Those samples were filtered according to the following criteria for further analysis: (1) the sample provides complete clinical information (including age, gender, and BMI); (2) the sequencing data have been previously published and available; and (3) the clinical information matches the sequencing data.

2.2. Data Analysis and Microbiome Profiling

Sequence Read Archive (SRA) datasets or compressed files of raw reads were downloaded from NCBI or EBI. As described in Figure 1 (Figure 1), the workflow of microbiome profiling includes the following three procedures: Prepare raw sequencing reads. Paired reads were extracted from SRA files using SRA Toolkit (v2.3.2-5) with the parameters “fastq-dump –split-3 –o.” Filter the raw sequences from human genome. The host human genome sequences would be excluded by aligning them to the human reference (hg19) using Burrows-Wheeler Aligner (BWA) “MEM” module with default parameters [20]. Construct the relative abundance matrix. To calculate the relative abundance of each sequence, the remaining reads were aligned to the integrated gene catalog (IGC) reference, which consists of 9.8 million genes with BWA [9]. Finally, the relative abundance was calculated according to the following formula [15], using an in-house script:Here is the relative abundance of gene , is the number of reads that were aligned to gene of IGC reference, is the length of gene , and is the total number of genes.

Figure 1: The pipeline of data mining procedures. The whole pipeline of this study consists of preprocessing data (SRA to FASTQ, clinical information available, and discarding samples without complete clinical information), aligning to IGC and constructing the abundance matrix, feature selection and training algorithm, and biological interpretation.
2.3. Feature Selection

One key problem that we are trying to solve in this study is to identify disease-associated biomarkers in human microbiome and to build a classifier with the selected biomarkers that can classify a single sample with accurate diagnostic disease status. In this study, there are three challenges in biomarkers selection. The first one is that the selected feature list is susceptible to the training data and is sensitive to different methods due to a small sample size in high-dimensional features. The second challenge is that the process of feature selection will be performed on multiclass datasets, which is much more difficult than on binary datasets. The last challenge is the complexity of computation. Some methods have been reported, such as Max-Relevance and Min-Redundancy (mRMR) and iterative sure independence screening (ISIS), which are robust algorithms in feature selection and have been applied in T2D biomarker discovery [21]. Considering the limited computational resource and the capability of multiclass classification on higher-dimensional data, we thus introduced mRMR to be the feature selection method here. mRMR utilizes two criteria to select features. The maximum relevance means to screen out those features that characterize the labels of samples optimally, and sequentially the minimum redundancy is introduced to choose features that are maximally dissimilar to the already known ones. Top 500 features in the mRMR feature list were analyzed in the subsequent classification in this study.

2.4. Imbalanced Datasets between Classes

Referring to bioinformatics and genomics, a common problem in the application of classification algorithms is the imbalanced sample size of classes, because machine learning algorithms naturally work well on quite equal number of samples in each class. At present, datasets with skewed labels are becoming more and more frequent, and unbalanced samples have been reported to generate wrong prediction models. One approach to address this plight is to resample the dataset to offset this imbalance to generate a more robust and fair decision boundary. Random resampling of the data is a common way to solve those problems as seen in other studies [12]. While the random sampling methods have several drawbacks such as biased results and being laborious and time-consuming, some optimized methods have been developed with better performance [13, 22]. Undersampling the majority classes and oversampling the minority classes are two categories of resampling techniques.

In this study, NearMiss (version 3) and SMOTEENN were applied to solve unbalanced datasets. NearMiss uses the -nearest neighbor (KNN) classifier to achieve undersampling. It selects a given number of the closest majority samples for each minority sample to guarantee that every minority sample is surrounded by some majority samples. On the other hand, SMOTEENN, the method of oversampling, tends to remove examples from both classes by filtering the misclassified example’s three nearest neighbors from the training dataset. Thus, SMOTEENN shows an improved performance on datasets with a small number of positive samples [13].

2.5. Accuracy Assessment

More than one measurement can be used to assess the performance of machine learning algorithms, such as score and the receiver operating characteristic curve (ROC) [23]. The score can be interpreted as a weighted average of precision and recall. Correspondingly, the score reaches its best value at 1 and worst value at 0. The relative contributions of precision and recall to the score are equal. The formula for the score is

In the multiclass and multilabel case, the weighted average of the score of each class is computed. To evaluate the precision of the algorithm, the AUC was also calculated for comparison with that in previous work.

2.6. Accuracy Estimation of Selected Biomarkers

The prediction performance of the selected biomarkers should be measured by training a classifier on the data that are restricted to the selected biomarkers. Seven popular classifiers, -nearest neighbors (KNN), logistic regression (LR), random forest (RF), support vector machine (SVM), gradient boosting decision tree (GBDT), stochastic gradient descent (SGD), and adaptive boosting (ADA), were performed on the selected features. Based on the score, the best algorithm of classification on this dataset was selected in this analysis.

The classification of multiclasses can be implemented as known one-versus-all strategy [24]. This strategy consists of fitting one classifier per class, with the samples of that class as positives and other samples as negatives. Then these basic classifiers are used to predict a single sample by aggregating their decisions.

Scikit-learn method is a simple but efficient tool for data mining and subsequent data analysis. It brings machine learning to nonspecialists with a general-purpose high-level language, where the API of seven algorithms was provided [25].

2.7. Cross-Validation Method

We utilized the 5-fold cross-validation approach to evaluate the performance of prediction models. Following this method, the dataset was randomly divided into five equal-sized partitions. Each time we fitted a model on four partitions as the training data and test it on the remaining partition. The process was repeated five times and the average result was used to generate the estimation.

2.8. Biological Interpretation

The integrated gene catalog that has been annotated is a comprehensive resource for metagenomics analyses [9]. To understand these distinct feature-related biomarker genes, we annotated these genes with an in-house script and picked out the items from the profile dataset of genes, genus, KEGG, and eggNOG, which can be downloaded from the following website: http://meta.genomics.cn/meta/dataTools. Finally, we interpreted the biomarkers with details in multiple levels, such as Phylum, genus, KEGG, and eggNOG.

To address the problem whether these biomarkers are specifically associated with a particular phenotype, Venn diagram was plotted with the online program (http://genevenn.sourceforge.net/).

3. Results and Discussions

3.1. Sample Information

In this study, 806 samples with different phenotypes were analyzed. Those samples were from three large available metagenome-wide association studies at present, which focused on type 2 diabetes (T2D), liver cirrhosis, and rheumatoid arthritis (RA) (Table S1) [7, 1417]. Of these 806 samples, 383 were with normal phenotypes while 423 with different disease status were defined as abnormal ones and were labeled with different tags (normal status was labeled as class 0, type 2 diabetes as class 1, rheumatoid arthritis as class 2, and liver cirrhosis as class 3) as shown in Table 1. To investigate the impact of age and BMI, we used the one-sample Kolmogorov-Smirnov test to test the normality of the age and BMI distribution. And then one-way analysis of variance (ANOVA) from “scipy.stats” package was adopted to perform significance test. In previous studies, aging was found to have a global impact on the physiology of human digestion system, and the process of aging can affect the composition of the human gut microbiome [8, 26]. In this study, age is associated with phenotype (ANOVA test, ) among these four groups. However, in line with previous findings, the body mass index (BMI) shows no significantly direct correlation with phenotype (ANOVA test, ) [27], and BMI is used as a screening tool to estimate the health status instead of as a diagnostic tool for disease risk. To concern these related factors on phenotype, age, gender, and BMI were included. The number of samples labeled as 0 is more than other classes, so imbalanced datasets would be one special case for classification problem in this study (Table 1).

Table 1: Statistics on sample information.
3.2. Data Alignment and Matrix Construction

A total of 383 samples with normal phenotype were sequenced to generate an average of 9.30 (±4.02) Gb data, and 92% of these were aligned to the IGC reference. In addition, 423 samples with abnormal phenotype generated a total of 11.27 (±6.95) Gb data also with an average mapping rate of 92% (Table S2). The strategy of mapping to the IGC reference was adopted. The catalog of IGC may reach saturated coverage of core gene content and the pipeline of assembly and annotation would cost more time and computing resource [9]. Finally, an 806 × 9,879,896 relative abundance matrix was constructed by calculating the number of reads that aligned to the IGC catalog and was normalized with the length of gene.

3.3. Feature Selection and Accuracy of the Algorithm

The relative matrix was preprocessed at the initial stage of data analysis, and the frequency of gene less than 90% in population scale was removed, and then 13,990 genes remained in the next stage. Following this, we applied mRMR methods to features selection. The top 500 features in the mRMR feature list were selected as candidate biomarkers. After two criteria of mRMR process, we iteratively repeated the procedure of prediction with seven algorithms to increase the number of biomarkers, which initiated with 100 features and finally increased by 20 features (Table S3). The performance of SOMTEENN was overall better than NearMiss3 on all these seven algorithms, especially the logistic regression and support vector machine. The largest value of optimized combination achieved 0.9142 on 300 features with the logistic regression (Table 2). In the cross-validation, the average value was 0.92 (±0.01), supporting the fact that this model is stable and accurate (Figure 2).

Table 2: The evaluation of algorithms based on 1 score.
Figure 2: The ROC plot of 5-fold cross-validations. For repeated cross-validation, multiple curves were plotted with green color, where each class was repeated five times, and the mean curve was plotted with blue color. The confidence interval was fulfilled with grey color.

Compared with previous studies, we plotted the ROC and calculated the area under the ROC curve (AUC). The AUC that we obtained with optimized features and algorithm was 0.85 in normal samples, 0.99 in type 2 diabetes, 0.96 in rheumatoid arthritis, and 0.99 in liver cirrhosis, respectively (Figure 3). The AUC value of the control group was lower than the other groups and may be attributed to the diversity of the human microbiome in healthy people. The AUC of Chinese T2D metagenome from previous study was 0.81, which was calculated based on the species and metagenomics clusters (MGC). The AUCs of rheumatoid arthritis and liver cirrhosis were published as 0.94 and 0.83, respectively. The better performance of our biomarkers can be explained by the relatively complete reference, more features than previous studies, and integrated clinical information [1517].

Figure 3: The ROC plot and AUC. Multicurves were plotted on the same figure, (normal class = blue, type 2 diabetes = green, rheumatoid arthritis = red, and liver cirrhosis = aqua-blue); the numbers in parenthesis were values of AUC.

The predictive power of these microbiome biomarkers is promising to be applied to disease diagnostics, especially disease screening within large-scale population. The multiclass classifier is a potential tool in the personalized medicine in a wide spectrum of phenotypes.

3.4. Biological Interpretability

ANOVA was performed to investigate the difference among these four groups, and these biomarkers showed significant difference (Table S4). These biomarkers were annotated as 67 Bacteroidetes, 207 Firmicutes, 3 Proteobacteria, and 21 biomarkers that were not classified into phylum. The highest proportion of biomarkers was Firmicutes, which was correlated with the fat storage and energy harvest [28].

It is well accepted that many factors synergistically shape the diversity of human microbiome, so phenotype-specific biomarkers are important in the evaluation of health status. The top biomarker in the ranking list according to the P value from ANOVA is DOF003_GL0053139, which is a gene from the phylum of Firmicutes and the genus of Clostridium. The relative abundance of this biomarker in RA is higher than that of others, which is consistent with the previous report [29]. The biomarker 469590.BSCG_05503 is overexpressed in liver cirrhosis, which belongs to Bacteroides with its function being well studied [30].

To comprehensively investigate the phenotype-related microbiome, we performed pairwise comparison of the relative abundance of one phenotype against other phenotypes. A total of 117 biomarkers were significantly correlated with T2D disease, 257 biomarkers with rheumatoid arthritis, and 220 biomarkers with liver cirrhosis, respectively (Table S5).

In the T2D group, 48 biomarkers were enriched and 81% of them were Bacteroides involved in the process of protein and carbohydrate breakdown. However, 69 of 117 biomarkers were significantly reduced in T2D disease and 78% of them were Firmicutes known to enhance the absorption of fat. All these results support the previous finding that the ratio of Bacteroides/Firmicutes was altered in T2D patients [31].

In the early stage of RA, Gram-positive bacteria were enriched and Gram-negative bacteria were depleted in the microbiome of RA patients [17]. However, two biomarkers of Veillonella, one biomarker of Haemophilus, and 74% of 73 depleted biomarkers of Bacteroides are Gram-negative bacteria in this study. Besides, we also identified 184 RA-enriched biomarkers including 16 Clostridium, 16 Faecalibacterium, and 10 Blautia as well as Coprococcus, Enterococcus, Eubacterium, and Roseburia. The component of biomarkers found in RA patients was also consistently reported from previous results.

In the liver cirrhosis group, 66 biomarkers were found to be enriched in liver cirrhosis, including 8 Streptococcus and 2 Veillonella, in line with previous work. It indicated that these two genera should play important roles in the development of liver cirrhosis. Interestingly, another important genus, Faecalibacterium, an anti-inflammatory commensal bacterium, was also found [16, 32].

In this study, 154 biomarkers were depleted in the microbiome of liver cirrhosis patients, which is similar to the RA enriched biomarkers, including 12 with Clostridium, 13 with Faecalibacterium, and 10 with Blautia as well as Coprococcus, Enterococcus, Eubacterium, Roseburia, and Ruminococcus. The phenomenon may be attributed to interindividual variation in the human microbiome.

To address the problem of whether these biomarkers are specific to a specific phenotype, Venn diagram shows that no biomarker was shared by three phenotypes. T2D shared 78 biomarkers with RA and 39 biomarkers with liver cirrhosis. RA shared 179 biomarkers with liver cirrhosis. No biomarker was specific to T2D and RA. It is noteworthy that two biomarkers were specific to liver cirrhosis. These two biomarkers are 1000570.HMPREF9966_1928 and 1000570.HMPREF9966_1926, which belong to Streptococcus (Figure 4). Previous study has reported that Streptococcus is associated with liver disease [33].

Figure 4: The Venn diagram of phenotype-specific biomarkers. Each circle stands for one phenotype and the number stands for biomarkers.

To characterize the functional role of microbiome in phenotype, we annotated each biomarker by the KEGG database. The pathways that include more than 10 biomarkers are related to membrane transport, genetic information processing, translation, and cellular processes and signaling (Table S4). Particularly, the membrane transporters were found to be enriched in type 2 diabetes [14, 15]. We also found that some biomarkers were enriched in the pathways of carbohydrate metabolism, lipid metabolism, and amino acid metabolism, which are crucial fundamental physiological processes for living life.

4. Discussions and Conclusions

In this study, we improved the performance of classification considering the problem of multiclass and imbalanced datasets by using the technology of machine learning. A total of 300 biomarkers were selected from 13,990 features including clinical information and the matrix of relative gene abundance from 806 human microbiomes through using logistic regression classifier. The phenotype-specific biomarkers were interpreted comprehensively.

Our study pinpoints the potential role of human microbiome which may lead to the research and development of microbiome knowledge-based personalized precision medicine by monitoring and modulating the diversity of microbiome [11]. While in this manuscript we selected biomarkers from Chinese individuals, all these selected biomarkers should be validated on other populations and then can be exploited in the future personal medicine.

This study uncovered some interesting phenomena. First, the ratio of Bacteroides/Firmicutes was altered in T2D patients, which is in line with previous study. Second, these biomarkers related to RA were found to be depleted in Gram-negative bacteria and enriched in Gram-positive bacteria. Third, Streptococcus and Veillonella were found to be enriched in liver cirrhosis and, particularly, another bacterium named Faecalibacterium was also found. Finally, some bacteria were enriched in one phenotype; however, they were depleted in anther phenotype. All these findings suggest the complexity of microbiome and give us cues to treat microbiome according to the status of one person.

Disease-specific biomarkers were analyzed in this study and two liver cirrhosis specific biomarkers that belong to the genus of Streptococcus were found. Other biomarkers were shared by two types of disease. Many intriguing biomarkers play important role in the core metabolism process, such as O2.UC32-1_GL0019091 that belongs to K03088, which is thought to play a role in enhancing transcriptional specificity in low-G+C Gram-positive bacteria [34].

Incorporated with clinical information and the profiling of human microbiome, the healthy status of one individual can be predicted. However, those previously published microbiome-based biomarkers were almost specific to the population studies; it means that all these biomarkers should be demonstrated in other independent study cohorts. Other pitfalls and challenges such as the diversity of dietary, standardization in sample collection and the treatment, and different microbiome analytic tools play a vital, important role in the utility of microbiome. Thus, comprehensive information of human and standardization in the pipeline of treatment will be the promise, which will ensure the machine learning to tap its full potential.

Finally, this study reveals that gut microbiome biomarkers are able to distinguish abnormal cases from controls with a higher level of specificity than previous results [7, 12, 14, 16, 17]. Our method can also be extended from the features to the type of abnormal phenotype, which will score the possibility of the specific disorder and make the mode of monitoring gut health become a reality. The technology of machine learning will accelerate the speed of the application of human microbiome.

Conflicts of Interest

The authors declare that they have no conflicts of interest.

Authors’ Contributions

Honglong Wu and Lihua Cai contributed equally to this work.

Acknowledgments

The authors gratefully acknowledge support from the Science and Technology Program of Tianjin, China (no. 13ZCZDSY02500).

Supplementary Materials

S1 Table. The detailed clinical information of 806 samples. S2 Table. The statistical information of data analysis. S3 Table. Evaluation of the seven algorithms. S4 Table. The annotated and ANOVA test information of 300 biomarkers. S5 Table. Pairwise comparison of biomarkers. (Supplementary Materials)

References

  1. Q. Feng, S. Liang, H. Jia et al., “Gut microbiome development along the colorectal adenoma-carcinoma sequence,” Nature Communications, vol. 6, article 6528, 2015. View at Publisher · View at Google Scholar · View at Scopus
  2. N. Kamada, S.-U. Seo, G. Y. Chen, and G. Núñez, “Role of the gut microbiota in immunity and inflammatory disease,” Nature Reviews Immunology, vol. 13, no. 5, pp. 321–335, 2013. View at Publisher · View at Google Scholar · View at Scopus
  3. E. Le Chatelier, T. Nielsen, and J. Qin, “Richness of human gut microbiome correlates with metabolic markers,” Nature, vol. 500, no. 7464, pp. 541–546, 2013. View at Publisher · View at Google Scholar
  4. L. Zitvogel, M. Ayyoub, B. Routy, and G. Kroemer, “Microbiome and anticancer immunosurveillance,” Cell, vol. 165, no. 2, pp. 276–287, 2016. View at Publisher · View at Google Scholar · View at Scopus
  5. The Human Microbiome Project Consortium, “Structure, function and diversity of the healthy human microbiome,” Nature, vol. 486, no. 7402, pp. 207–214, 2012. View at Publisher · View at Google Scholar
  6. J. Zhang, Z. Guo, Z. Xue et al., “A phylo-functional core of gut microbiota in healthy young Chinese cohorts across lifestyles, geography and ethnicities,” The ISME Journal, vol. 9, no. 9, pp. 1979–1990, 2015. View at Publisher · View at Google Scholar · View at Scopus
  7. F. Bäckhed, J. Roswall, Y. Peng, Q. Feng, H. Jia, P. Kovatcheva-Datchary et al., “Dynamics and stabilization of the human gut microbiome during the first year of life,” Cell Host & Microbe, vol. 17, no. 5, pp. 690–703, 2015. View at Publisher · View at Google Scholar
  8. T. Yatsunenko, F. E. Rey, M. J. Manary et al., “Human gut microbiome viewed across age and geography,” Nature, vol. 486, no. 7402, pp. 222–227, 2012. View at Publisher · View at Google Scholar · View at Scopus
  9. J. Li, H. Jia, X. Cai, H. Zhong, Q. Feng, S. Sunagawa et al., “An integrated catalog of reference genes in the human gut microbiome,” Nature Biotechnology, vol. 32, no. 8, pp. 834–841, 2014. View at Publisher · View at Google Scholar
  10. J. Qin, R. Li, J. Raes, M. Arumugam, K. S. Burgdorf, C. Manichanh et al., “A human gut microbial gene catalogue established by metagenomic sequencing,” Nature, vol. 464, no. 7285, pp. 59–65, 2010. View at Publisher · View at Google Scholar
  11. N. Zmora, D. Zeevi, T. Korem, E. Segal, and E. Elinav, “Taking it personally: personalized utilization of the human microbiome in health and disease,” Cell Host & Microbe, vol. 19, no. 1, pp. 12–20, 2016. View at Publisher · View at Google Scholar · View at Scopus
  12. M. W. Libbrecht and W. S. Noble, “Machine learning applications in genetics and genomics,” Nature Reviews Genetics, vol. 16, no. 6, pp. 321–332, 2015. View at Publisher · View at Google Scholar
  13. G. E. A. P. A. Batista, R. C. Prati, and M. C. Monard, “A study of the behavior of several methods for balancing machine learning training data,” ACM SIGKDD Explorations Newsletter, vol. 6, no. 1, pp. 20–29, 2004. View at Publisher · View at Google Scholar
  14. F. H. Karlsson, V. Tremaroli, I. Nookaew et al., “Gut metagenome in European women with normal, impaired and diabetic glucose control,” Nature, vol. 498, no. 7452, pp. 99–103, 2013. View at Publisher · View at Google Scholar · View at Scopus
  15. J. J. Qin, Y. R. Li, Z. M. Cai et al., “A metagenome-wide association study of gut microbiota in type 2 diabetes,” Nature, vol. 490, no. 7418, pp. 55–60, 2012. View at Publisher · View at Google Scholar · View at Scopus
  16. N. Qin, F. Yang, and A. Li, “Alterations of the human gut microbiome in liver cirrhosis,” Nature, vol. 513, no. 7516, pp. 59–64, 2014. View at Publisher · View at Google Scholar
  17. X. Zhang, D. Zhang, H. Jia et al., “The oral and gut microbiomes are perturbed in rheumatoid arthritis and partly normalized after treatment,” Nature Medicine, vol. 21, no. 8, pp. 895–905, 2015. View at Publisher · View at Google Scholar
  18. J. A. Gilbert, R. A. Quinn, J. Debelius et al., “Microbiome-wide association studies link dynamic microbial consortia to disease,” Nature, vol. 535, no. 7610, pp. 94–103, 2016. View at Publisher · View at Google Scholar · View at Scopus
  19. J. Wang and H. Jia, “Metagenome-wide association studies: fine-mining the microbiome,” Nature Reviews Microbiology, vol. 14, no. 8, pp. 508–522, 2016. View at Publisher · View at Google Scholar
  20. H. Li and R. Durbin, “Fast and accurate short read alignment with Burrows-Wheeler transform,” Bioinformatics, vol. 25, no. 14, pp. 1754–1760, 2009. View at Publisher · View at Google Scholar · View at Scopus
  21. L. Cai, H. Wu, D. Li, K. Zhou, and F. Zou, “Type 2 diabetes biomarkers of human gut microbiota selected via iterative sure independent screening method,” PLoS One, vol. 10, no. 10, Article ID e0140827, 2015. View at Publisher · View at Google Scholar
  22. J. P. Zhang and I. Mani, “KNN approach to unbalanced data distributions: a case study involving information extraction,” in Proceedings of the Proceeding of International Conference on Machine Learning, Workshop on Learning from Imbalanced Data Sets (ICML '03), Washington DC, Wash, USA, August 2003.
  23. D. M. W. Powers, “Evaluation: From Precision, Recall and F-Measure to ROC, Informedness, Markedness & Correlation,” 2011.
  24. C. M. Bishop, Pattern Recognition and Machine Learning, Springer, New York, NY, USA, 2006. View at MathSciNet
  25. F. E. A. Pedregosa, “Scikit-learn: machine learning in python,” Journal of Machine Learning Research, vol. 12, pp. 2825–2830, 2011. View at Google Scholar · View at MathSciNet
  26. E. Biagi, L. Nylund, M. Candela et al., “Through ageing, and beyond: gut microbiota and inflammatory status in seniors and centenarians,” PLoS ONE, vol. 5, Articale ID e10667, 2010. View at Publisher · View at Google Scholar · View at Scopus
  27. D. Rojo, A. Hevia, R. Bargiela et al., “Ranking the impact of human health disorders on gut metabolism: systemic lupus erythematosus and obesity as study cases,” Scientific Reports, vol. 5, article 8310, 2015. View at Publisher · View at Google Scholar · View at Scopus
  28. P. J. Turnbaugh, R. E. Ley, M. A. Mahowald, V. Magrini, E. R. Mardis, and J. I. Gordon, “An obesity-associated gut microbiome with increased capacity for energy harvest,” Nature, vol. 444, no. 7122, pp. 1027–1031, 2006. View at Publisher · View at Google Scholar · View at Scopus
  29. A. C. Essenmacher, N. Khurram, and G. T. Bismack, “A case of reactive arthritis due to Clostridium difficile colitis,” Journal of Community Hospital Internal Medicine Perspect, vol. 6, no. 1, p. 30151, 2016. View at Publisher · View at Google Scholar
  30. H. Lu, G. Qian, Z. Ren et al., “Alterations of Bacteroides sp., Neisseria sp., Actinomyces sp., and Streptococcus sp. populations in the oropharyngeal microbiome are associated with liver cirrhosis and pneumonia,” BMC Infectious Diseases, vol. 15, no. 1, article 239, 2015. View at Publisher · View at Google Scholar
  31. N. Larsen, F. K. Vogensen, F. W. J. Van Den Berg et al., “Gut microbiota in human adults with type 2 diabetes differs from non-diabetic adults,” PLoS ONE, vol. 5, no. 2, Article ID e9085, 2010. View at Publisher · View at Google Scholar · View at Scopus
  32. H. Sokol, B. Pigneur, L. Watterlot et al., “Faecalibacterium prausnitzii is an anti-inflammatory commensal bacterium identified by gut microbiota analysis of Crohn disease patients,” Proceedings of the National Acadamy of Sciences of the United States of America, vol. 105, no. 43, pp. 16731–16736, 2008. View at Publisher · View at Google Scholar · View at Scopus
  33. A. Gonzlez-Quintela, C. Martinez-Rey, J. F. Castroagudin, M. C. Rajo-Iglesias, and M. J. Dominguez-Santalla, “Prevalence of liver disease in patients with Streptococcus bovis bacteraemia,” Journal of Infection, vol. 42, no. 2, pp. 116–119, 2001. View at Publisher · View at Google Scholar
  34. X. Xue, J. Tomasch, H. Sztajer, and I. Wagner-Döbler, “The delta subunit of RNA polymerase, RpoE, is a global modulator of Streptococcus mutans environmental adaptation,” Journal of Bacteriology, vol. 192, no. 19, pp. 5081–5092, 2010. View at Publisher · View at Google Scholar · View at Scopus