Skip to main content

Application of a hierarchical enzyme classification method reveals the role of gut microbiome in human metabolism

Abstract

Background

Enzymes are known as the molecular machines that drive the metabolism of an organism; hence identification of the full enzyme complement of an organism is essential to build the metabolic blueprint of that species as well as to understand the interplay of multiple species in an ecosystem. Experimental characterization of the enzymatic reactions of all enzymes in a genome is a tedious and expensive task. The problem is more pronounced in the metagenomic samples where even the species are not adequately cultured or characterized. Enzymes encoded by the gut microbiota play an essential role in the host metabolism; thus, warranting the need to accurately identify and annotate the full enzyme complements of species in the genomic and metagenomic projects. To fulfill this need, we develop and apply a method called ECemble, an ensemble approach to identify enzymes and enzyme classes and study the human gut metabolic pathways.

Results

ECemble method uses an ensemble of machine-learning methods to accurately model and predict enzymes from protein sequences and also identifies the enzyme classes and subclasses at the finest resolution. A tenfold cross-validation result shows accuracy between 97 and 99% at different levels in the hierarchy of enzyme classification, which is superior to comparable methods. We applied ECemble to predict the entire complements of enzymes from ten sequenced proteomes including the human proteome. We also applied this method to predict enzymes encoded by the human gut microbiome from gut metagenomic samples, and to study the role played by the microbe-derived enzymes in the human metabolism. After mapping the known and predicted enzymes to canonical human pathways, we identified 48 pathways that have at least one bacteria-encoded enzyme, which demonstrates the complementary role of gut microbiome in human gut metabolism. These pathways are primarily involved in metabolizing dietary nutrients such as carbohydrates, amino acids, lipids, cofactors and vitamins.

Conclusions

The ECemble method is able to hierarchically assign high quality enzyme annotations to genomic and metagenomic data. This study demonstrated the real application of ECemble to understand the indispensable role played by microbe-encoded enzymes in the healthy functioning of human metabolic systems.

Background

Enzymes represent a significant fraction of an individual proteome [1] and catalyze a variety of specific reactions in the cellular systems [2, 3]. Hence, identification of the functions of an entire complement of enzymes in an organism helps generate the metabolic blueprint of that species. This will not only improve our understanding of defined cellular processes of individual species but also help study the metabolic interdependence of multiple species in an ecosystem such as the human gut microbiome.

The Enzyme Commission (EC) [4] has classified all enzymes based on the enzymatic reactions they catalyze. Each enzyme has an EC number, which is a hierarchical number that distinguishes enzymes by the type of reactions they catalyze. The EC groups all enzymes into six broad classes that include (1) oxidoreductases - catalyze oxidoreduction reactions; (2) transferases - catalyze the transfer of a chemical group from a donor to an acceptor; (3) hydrolases - catalyze the hydrolysis of various bonds; (4) lyases - enzymes that cleave bonds by means other than by hydrolysis; (5) isomerases - catalyze geometrical or structural changes within one molecule; and (6) ligases - catalyze the joining of two molecules coupled with hydrolysis of a pyrophosphate bond in ATP or a similar triphosphate. The EC classification system assigns a unique four-field number to each enzymatic activity (such as EC 1.2.1.3 for aldehyde dehydrogenase (NAD+)) where, the first three numbers (a.k.a. levels) of an EC number represent progressively finer description of the enzymatic reaction, while the last level mostly represents substrate specificity of a reaction [5, 6].

Experimental characterization of the enzymatic reactions of all enzymes in a genome or a metagenome is a tedious and expensive task. With the exception of a few well-characterized genomes such as Escherichia coli (E. coli) and yeast, the fraction of experimentally annotated enzymes in many sequenced genomes is very small. The problem is more pronounced in the metagenomic samples where even the species are not adequately cultured or characterized. To address this problem, computational approaches, which can build accurate models from known data to predict the unknown data, have been employed. Such models have been widely used to predict protein functions and annotate newly sequenced genomes [7–10].

Several computational methods exist for predicting enzyme classes (refer to our review [11]). Many of them use electronically inferred annotations by transferring the annotations (EC number, enzyme name and the reactants) of a known enzyme to an unknown enzyme, if the features of unknown match with the known. Most of the methods differ by the type of features they use to match proteins, which broadly include amino acid composition [12], sequence or structure homology [13–15], domain composition [16, 17] and sequence motifs [18]. A variety of machine learning (ML) and data mining algorithms, including nearest-neighbor methods [19], support vector machines (SVMs) [20, 21], Bayesian [22, 23] and ensemble approaches [24, 25] have been employed to build models for enzyme classification. The performance of these methods varies based on the classification algorithm, input datasets and features used for model building. Nevertheless, most of the existing methods fail to predict EC levels 3 and 4 due to the increasing difficulty in predicting the finer levels in the hierarchy, thus offer only limited value to enzyme annotations. Many enzyme prediction methods exhibit a lack of balance between specificity and sensitivity; EzyPred [17] and EFICAz2.5 [24] offer limited sensitivity for certain enzyme classes and/or EC levels and run extremely slow with high false positive rates. Hence, the existing methods are inadequate for annotation of high-throughput sequences generated from genomic and metagenomic projects, warranting the need for developing a new computational method.

The current study requires the identification of all enzymes encoded by the gut microbiome in human because the enzymes encoded by the gut microbiota play an extensive role in the human metabolism. Human gut microbiome is the largest and most complex of all microbial communities that harbor human body, with a gene set that is about 150 times larger than that of the human gene set [26]. Human gut microbiome alone is estimated to contain about 1000-1500 different species [26, 27], but a majority of them are yet to be characterized. These bacterial communities extensively contribute to human gut metabolism by complementing enzymes that are not encoded by the human genome, but are essential for digestion of complex polysaccharides, absorption, metabolism of amino acids and vitamins, shaping of the immunological environment and a wide range of other metabolic functions [28–30]. Changes in the composition of human microbiota have been linked to health conditions such as inflammatory bowel disease (IBD), antibiotic-resistant infections, obesity, colon cancer, symptomatic atherosclerosis and diabetes [31–33]. Hence, the identification and functional characterization of gut microbial enzymes is a very important step towards understanding the microbe-dependent component of the human metabolism.

In this study, we developed a new hierarchical enzyme classification method based on machine learning that accurately predicts if a protein sequence is an enzyme or a non-enzyme, and if an enzyme, what is the specific enzymatic reaction (class and subclasses) it carries out. We apply this method to identify the full enzyme complements of 10 sequenced genomes of model organisms, and also those of the microbial species in the metagenomic samples obtained from human gut microbiomes. The methodology developed in this project and its application to identifying the full complements of microbial enzymes has made it feasible to study the role of microbe-complemented enzymes in the human gut metabolism. To our knowledge, this study represents a novel and robust approach to studying the pathway level host-pathogen interactions in the human gut metabolism.

Results and discussion

We describe the results and discussion in two separate sections that include the method development followed by its application to study the pathways in human gut metabolism. Figure 1 shows a schematic of the method and its application. An ensemble of five different machine-learning (ML) classifiers was used to build prediction models based on protein domains that include sequence or structure-derived features. Hence, this method is named as ECemble (Enzyme Classification using ensemble approach). Bayesian network [22] model represents a set of random variables and their conditional dependencies, whereas Naïve Bayes [34] works under the strong assumption that there is independence among the features. K-nearest neighbor (KNN) [35] is 'instance based' learning algorithm. Despite its simplicity, it can offer very good performance. We tried polynomial and linear kernels in SVM [36] and polynomial kernel was found to be working better than the linear kernel. Random Forest classifiers (RFC) [37] work by generating large number of decision trees in a specific random way such that each one is de-correlated with the others. Each classifier model generates a probability distribution for all classes for each instance in the train and test set and the class with the highest probability is assigned as the predicted class.

Figure 1
figure 1

Schematic representation of the ECemble method and its application.

This is a hierarchical prediction method that predicts enzymes and enzyme classes at 5 levels (henceforth referred to as L0, L1, L2, L3 and L4), where the first model at L0 predicts if a protein is an enzyme or a non-enzyme, and the subsequent models from L1 to L4 predict specific class and subclass of an enzyme in the EC number hierarchy. The predictions are selected by a consensus approach, i.e., only when at least two of the three top-performing ML classifiers show consistent predictions. To demonstrate its usefulness, we applied this method to predict the full complement of enzymes from ten sequenced proteomes. In addition, we also tested about 2.5 million protein sequences that were obtained from metagenomic sequencing and assembly of the human gut microbiome [26]. Using the full complements of enzymes from human and human gut microbiomes, we further investigated the role of microbe-derived enzymes in the human gut metabolism.

Development of ECemble method

Feature selection and optimization

We collected the known enzyme and non-enzyme protein sequences as positive and negative datasets, respectively. Sequences were clustered at 70% identity to generate 64,948 enzyme and 128,292 non-enzyme sequences, where each sequence has mapping to at least one of the following databases. These sequences were mapped against three domain databases that include Pfam [38], Superfamily [39] and Prosite [40] to extract the enzyme-specific and non-enzyme-specific features (as discussed in the Methods section). A protein domain is a conserved part of a protein's sequence or structure that can evolve, function, and exist independently of the rest of the protein chain. We used three domain databases that include Pfam, Superfamily and Prosite to represent the functional, structural and motif or active site regions of protein domains, respectively. These domains are used as features, where, in combination they offer a comprehensive feature set for ML methods to discriminate between enzymes and non-enzymes. We extracted 5,045 and 9,196 overlapping features for enzyme and non-enzyme datasets, respectively, from three domain databases (Table 1).

Table 1 Distribution of enzyme and non-enzyme features.

To facilitate the building of an effective ML-based classifier, it is essential to represent the domain features of each sequence in a binary format (feature vector), using '1' for the presence and '0' for the absence of a domain. The size of the feature vector matrix is N × ∑ (P, Q, R), where N is the total number of enzyme and non-enzyme sequences, and P, Q and R represent the total number of mapped domains in Pfam, Superfamily and Prosite, respectively (Table 2). In theory, the dimensionality of this vector gets very large; however, since each protein sequence contains only a few domains, we used the sparse-formatting option by storing the occurrence of features with their locations in the feature space to significantly speed up the model building process.

Table 2 Feature vectors and dimensionality for the dataset.

Hierarchical design of prediction models

We applied five different machine-learning (ML) algorithms to best exploit features from the training dataset and optimize at each level. These include Naïve Bayes, k-Nearest Neighbor (KNN) classifier, Support Vector Machine (SVM), Decision Stump (DS) [41] and Random Forest (tree-based) classifiers (RFC). ML algorithms are employed to learn discriminative features of classes from the training data, build models, and test how related the unknown instances (testing data) are to these models. We used the WEKA (Waikato Environment for Knowledge Analysis) [42] framework to build prediction models in an iterative fashion at 5 different levels (L0 to L4). As part of the enzyme identification step, the first model at L0 predicts if a protein sequence is an enzyme or not. Only the sequences predicted as enzymes at L0 are forwarded to build models for predicting enzyme classes and subclasses at L1-L4, sequentially. At L0, there are only two classes (enzyme vs. non-enzyme) and similarly, at L1, there are only 6 enzyme classes; hence, one model is sufficient to predict classes at these two levels. However, the six enzyme classes at L1 are further divided into 51 subclasses at L2. As a result, 6 prediction models are constructed to predict 51 subclasses at L2. Similarly, due to the increasing number of subclasses at each subsequent level, 51 and 169 models are constructed to predict 169 and 1,921 subclasses at L3 and L4, respectively (Table 3). This hierarchical design of prediction models ensures that subclasses of a superclass are not predicted outside of that superclass and hence minimize false positives. For instance, members of EC 1.x.x.x superclass are never predicted as members of EC 2.x.x.x, and so on and so forth.

Table 3 Overall prediction accuracy of ECemble method.

Evaluation of prediction performance

All feature vectors were randomly divided into 80 and 20 percent subsets for training and testing, respectively. Since the datasets are unbalanced across classes (and subclasses), class distributions are approximately preserved at all EC levels using stratified partitioning for training and testing sets. We used a two-step validation procedure that include determining a10-fold cross-validation accuracy on the training set, and the testing accuracy using the testing dataset that is not a part of the training data. We also report standard performance measures over each class at each level, including true positive rate (TPR), false positive rate (FPR), and receiver operating characteristic (ROC) curves and the area under the curve (AUC). Please refer to the Methods section for more details.

Table 4 shows the 10-fold accuracy and the testing accuracy for each of the five ML algorithms (Additional file 1: Figure S1). Overall, these two accuracies are consistent across the five classifiers indicating that models are not over trained. With the exception of decision stump (DS) classifier, the other four classifiers achieved at least 92.5% and 95.9% testing accuracies at L0 and L1, respectively. The three top performing classifiers are KNN, SVM and RFC, where testing accuracies reached at least 94.4% and 97.3% at L0 and L1, respectively. At the same time the false positive rates are very low; at L0 and L1, the FPR for the top three methods ranges from 0.054 to 0.061 and 0.011 to 0.026, respectively.

Table 4 Ten-fold cross validation and testing accuracy for enzyme identification and enzyme classification.

Figure 2 illustrates ROC curves that show the relationship between TPR (sensitivity) and FPR (1-specificity) for a single class. An ideal ROC curve heads straight up on the Y-axis and then to the right parallel to the X-axis of the graph; thus maximizing the area under the curve (AUC). Such curves indicate that the classifier is predicting maximum true positives with minimum false positives, with AUC values closer to one. Figure 2 shows fairly consistent ROC curves for the top performing three ML methods at L0 and L1, respectively. At L1, enzymes from EC6 class are consistently the best predicted, while those from EC2 class showed relatively the least performance. This is probably because EC2 is one of the largest classes with most divergent subclass distribution, while EC6 has the least number of subclasses. The minimum AUC values at L0 and L1 are 0.945 and 0.989, respectively, indicating the superior performance of these three classifiers; hence, we chose only these three classifiers (KNN, SVM, RFC) for further use in this study.

Figure 2
figure 2

ROC curve at Level-0 and Level-1 using top 3 performing classifiers. (A) Testing at Level-0 using KNN, (B) Testing at level-0 using RFC, (C) Testing at Level-0 using SVM, (D) Testing at Level-1 using KNN, (E) Testing at level-1 using RFC, (F) Testing at Level-1 using SVM. Due to the high accuracies, False Positive Rate (X-axis) is shown till 0.5 for Level-0 and till 0.05 for Level-1.

Consensus-based ensemble approach

A consensus approach adds confidence to the prediction accuracy and drastically reduces the false positives; hence, we implemented it by considering only those predictions, where the same enzyme class is predicted for a sequence by at least two out of the top three classifiers. Table 3 shows the overall prediction accuracy for testing data at each level (L0-L4) after using the consensus approach. We achieved an overall accuracy of 97.7% at L0 for identifying enzymes and non-enzymes, and at least 98.7% accuracy using models at L1 to L4 for predicting enzyme classes and subclasses. These results are very promising despite the fact that the size of the training data per model kept diminishing as the number of models increase at the lower levels (L2-L4). It can also be noted that at L0, 96.7% of the correctly predicted instances are consistent across the top three classifiers, while at L1-L4, over 99% of the correctly predicted instances are consistent (Additional file 2: Table S1).

Effect of sequence identity in the training dataset

Sequence redundancy in the training datasets often results in overtraining and inaccurate estimation of prediction accuracy. To test the effect of sequence identity on the prediction accuracy, we created four datasets at 70%, 60%, 50% and 40% sequence identities using the CD-HIT clustering algorithm [43] and accordingly labeled as cdh70, cdh60, cdh50, and cdh40. At lower sequence identity thresholds like 40%, more sequences got removed resulting in a fewer number of sequences in each enzyme class compared to datasets with higher sequence identities. Similarly, the number of enzyme classes containing minimum number of sequences (ten) for model building started to go down from thresholds 70 to 40 percent identity (Additional file 3: Figure S2A). We generated the enzyme prediction models for EC levels L0-L4 for all the 4 datasets (shown in Additional file 4: Table S2). We needed at least 10 sequences in each subclass to build models using 10-fold cross validation. Accuracy is the highest for the cdh70 dataset (70% sequence identity) compared to all other datasets (cdh60, cdh50, cdh40) (Additional file 3: Figure S2B); hence, we used this dataset for model building. The 70% sequence identity is considered an optimal threshold in many other ML datasets, because the enzyme function starts to diverge quickly when the sequence identity is below 70% [5].

Comparison of ECemble with other methods

We compared the performance of ECemble with that of two existing methods, BLAST and EFICAz [24]. We chose these two methods because a number of methods are BLAST-based homology searching methods and EFICAz is a popular open-source tool. We used the same train and test datasets against these three methods to compare their performance. In the first step, a BLAST database was created with train dataset and the test data was queried against it for enzyme identification. Because BLAST generates a number of hits for each query with varying levels of confidence, only the top hit was considered (with a minimum E-value threshold of 10-5 for blastp) as the prediction for each sequence in the test set. To predict the enzyme classes and subclasses, we performed a second query against a BLAST database created using only the enzyme sequences and used the same procedure as in the first step. Similarly, EFICAz method was trained and tested for both enzyme identification and classification. We also used four different clustered datasets (using CD-HIT) that were described earlier. For enzyme identification, ECemble reports the highest accuracy (94.9%) compared to BLAST (89%) and EFICAz (88.9%) using cdh70 dataset (Figure 3A). It can also be seen that the accuracy goes down as the percent identity in the datasets goes down from 70 to 40; however, this effect is more pronounced in the BLAST method compared to ECemble and EFICAz suggesting that reduced identity has minimal effect on these two methods. Similarly for enzyme classification at L1, the overall accuracy of ECemble (98.4%) is better than BLAST (28.43%) and EFICAz (93.36%) methods using cdh70 dataset (Figure 3B). We compared four different sequence identity thresholds (70%, 60%, 50% and 40%) of this dataset and our method performed better than others irrespective of the dataset used. These results convincingly demonstrate that the ECemble method consistently performed better than BLAST and EFICAz methods at all identity thresholds (70%, 60%, 50%, 40%), and that it is also suitable for accurate annotation of protein sequences with low sequence identity. Hence, in the next section, we used ECemble to identify and annotate complete enzyme complements of the sequenced genomes and metagenomes and applied this method to study the role of gut microbiome-derived enzymes in human metabolism.

Figure 3
figure 3

Comparison of ECemble with BLAST and EFICAz methods.

Application of ECemble method

Annotation of full complements of enzymes in sequenced genomes

The ability of ECemble to predict solely based on a protein sequence enables it to annotate full complements of enzymes from sequenced genomes as well as from mixed genomes such as metagenomic samples. We applied ECemble to annotate a large number of unknown or unclassified enzyme sequences from 10 proteomes (that include reviewed and unreviewed proteins from UniProt database [44]) of both eukaryotic and prokaryotic model organisms (Additional file 5). As seen in Figure 4, E. coli and yeast proteomes have a high fraction of known enzymes (27% and 24%, respectively) compared to less than 5% in most of the other species. In contrast, the fractions of ECemble predicted enzymes are smaller in E. coli and yeast (around 4.5%), but reach up to 12% in chicken and zebrafish proteomes. In human and mouse, the fractions of unknown or unclassified enzymes predicted by ECemble account to 7.4% and 8.1% of their proteomes, respectively. These data underscore an important and generic application of our method that is to annotate a large number of unknown enzymes and enzyme classes in the sequenced genomes.

Figure 4
figure 4

Fractions of known and ECemble predicted enzymes in the proteomes of 10 model organisms from UniProt. The data are sorted from highest to lowest known enzyme fractions. Both reviewed and unreviewed sequences from UniProt were used.

We further investigated about 200 newly predicted enzymes from the reviewed human proteome sequences (from the SwissProt database [45]), that were either not known as enzymes or had partially annotated enzyme subclasses. These 200 enzymes map to 119 unique enzyme reactions (with unique EC numbers), of which, only 30 are not previously known to be in human. Table 5 lists these 30 enzymes with the current annotation in UniProt and the predicted annotation by ECemble. Note that despite the well-characterized nature of the reviewed proteins, ECemble method is able to complement or correct the existing enzyme annotations (marked with * in Table 5). For instance, human genes such as FOXRED1 (FAD-dependent oxidoreductase domain-containing protein 1) and SCCPDH (saccharopine dehydrogenase-like oxidoreductase) have been broadly known as oxidoreductases and accordingly labeled as 1.-.-.-, in the UniProt database. ECemble method predicted the specific reactions of these enzymes as malate dehydrogenase (quinone) [EC 1.1.5.4] and Saccharopine dehydrogenase (NADP(+), L-glutamate-forming) [EC 1.5.1.10], respectively. In certain cases (marked with ** in Table 5), ECemble predictions differ only at the subclass levels (L2-L4), while in a small number of cases (marked with *** in Table 5), the predictions also differ from the current annotations at the class level (L1). Hence, experimental validation of these predictions is worth pursuing in the future. Because ECemble prediction models are built only from reviewed enzyme sequences, the accuracy and coverage of prediction by our model will continue to improve as the newly characterized enzymes from experimental studies become available. These results prove that ECemble method that solely uses protein sequences for prediction, is highly promising for the identification and classification of full complements of enzymes in the sequenced genomes.

Table 5 ECemble predicted enzymes from reviewed human proteome.

Prediction of enzymes from the gut microbiome

To understand the role played by the most densely colonized human microbiome (gut microbiome) in human metabolism, we would require a catalogue of all the enzymes that can potentially exist in the human gut environment. There are several hundreds to thousands of microbial species that inhabit human gut and many of them are unknown, under characterized or not culturable in laboratory conditions, which would have made our goal impossible to accomplish. However, Next-Generation Sequencing (NGS) data from metagenomic samples could be used for de novo assembly of microbial genomes, and consequently for the prediction of the translated proteomes from the assembled scaffolds. We applied the ECemble method on 2.5 million proteins from metagenomic sequences (discussed in the Methods section) and assigned 213,313 (8.53%) sequences to 513 distinct enzymatic reactions (Additional file 6). Of these, 276 reactions are also encoded in human, leaving 237 reactions that are exclusively encoded by human gut microbial genomes. Of these, 222 enzymes are known gut bacterial enzymes in SwissProt database, while the remaining 15 enzymes are newly predicted by the ECemble method.

Role of gut microbial enzymes on human metabolism

Application of ECemble to identify the enzymes and enzyme classes in human and gut microbial species has enabled us to ask the following questions. Which gut microbial enzymes are involved in human metabolism? Which human pathways are partly or fully driven by microbe-derived enzymes? Are there any previously unknown microbial enzymes that are involved in human gut metabolism? Even though, gut microbiome can contains bacteria, fungi and other small eukaryotes, about 98% of the gut microbial enzymes used in this study originated from bacteria [26], hence we limited this study only to bacteria-derived enzymes. To answer these questions, we created 4 different sets of enzymes. (i) Known human enzyme reactions (with unique EC numbers) (1,271), (ii) Unknown human enzyme reactions that are predicted by ECemble (19), (iii) Known bacterial enzyme reactions in the gut microflora predicted by ECemble (222) and (iv) Unknown gut bacterial enzyme reactions that are predicted by ECemble (15). For human and bacteria, known enzymes were obtained from SwissProt and KEGG (Kyoto Encyclopedia of Genes and Genomes) [46] databases (Additional file 7 shows known human enzyme reactions).

We mapped both human- and gut bacteria-encoded enzymes to the KEGG reference pathways to identify 48 human metabolic pathways, where each pathway contains both human- and bacteria-encoded enzymes plus at least one of them is predicted by ECemble. We refer to them henceforth as gut microbe complemented (GMC) pathways. The first set (39 pathways) contains human enzymes that are complemented by known gut bacterial enzymes. This set serves to validate the known role of gut bacteria in human metabolism (Additional file 8: Figure S3). The second set (9 pathways) is same as the first one; in addition, contains predicted enzymes (previously unknown) from gut bacterial species. These nine GMC pathways reveal the role of newly discovered gut bacterial enzymes in human gut metabolism, which is made possible with the ECemble method (Additional file 9: Figure S4). All the 48 pathways are mapped with enzymes (EC numbers) using a color-coded format. Light red colored enzymes represent known human enzymes, pink represents unknown human enzymes that are predicted by ECemble, light green represents known bacterial enzyme reactions in the gut microflora predicted by ECemble and light blue represents unknown gut bacterial enzyme reactions that are predicted by ECemble.

Functional distribution of GMC pathways is shown in Figure 5A and a full list of 48 pathways is given in Additional file 10. GMC pathways contribute to the metabolism of a variety of nutrients that primarily include carbohydrates, amino acids, vitamins and cofactors. Thus, the metabolic role of GMC pathways closely matches with the known nutritional requirements of humans. Our results show that gut microbial enzymes also play a role in lipid and energy metabolism and also in the metabolism of terpenoids, polyketides and derived amino acids such as taurine, D-glutamine, etc. Our method predicted 27 Carbohydrate active enzymes (CAZymes) (Additional file 11) in the human gut microbiome that include carbohydrate esterases, glycoside hydrolases, glycosyl transferases and polysaccharide lyases, which are primarily involved in carbohydrate metabolic pathways. As shown in Figure 5B, human- and microbe-derived enzymes complement different categories of human metabolic pathways. Note that 64, 60 and 30 exclusively microbe-derived enzymes complement the amino acid, carbohydrate and vitamin/cofactor metabolic pathways, respectively, which underscores their pivotal role as well as the dependency of human pathways on microbe-derived enzymes. Microbes not only complement, but also supplement some of the enzyme functions that are common to both human and microbes. For instance, we found that 518 enzymatic reactions are shared between human and gut-bacterial species; and we hypothesize that some of these common enzymes that are secreted by bacteria could supplement to perform human metabolic functions, and vice-versa. These common enzymes are highlighted in blue color in the Additional file 12. These results empirically support the argument that symbiotic gut microbiome has coevolved with human (or the coelomate animals) and they play a huge role in the metabolic interactions between host-gut microbiota [32].

Figure 5
figure 5

Distribution of pathways and enzymes in each of the major pathway categories. (A) Distribution of mapped pathways in each of the 8 major pathway categories (B) Distribution of human and bacterial enzymes in each of the 8 major pathway categories.

Humans lack the enzymes needed to degrade oxalate [47]. It has been shown that the bacterial degradation of oxalate is carried by Oxalobacter formigenes in the human intestinal tract and the absence of Oxalobacter formigenes is considered a risk factor for urolithiasis [48]. Similarly Choline, an essential dietary nutrient, is found to be metabolized in the liver by the gut microbial enzymes. Thus conversion of dietary choline is used as a metabolic hallmark for liver and cardiovascular diseases [49]. Similarly, nitrate reductase that converts nitrate into nitrite and nitric oxide, is synthesized only by gut microbiome; elevated levels of nitric oxide have been associated with both IBD and obesity-induced insulin resistance [31]. Gut microbiome also plays a crucial role in the metabolism of xenobiotics; at least thirty commercially available drugs are shown to be metabolized as substrates by bacterial enzymes [50, 51]. On the other hand, gut microbiome is also a source of inflammatory molecules such as lipopolysaccharide and peptidoglycan that may contribute to metabolic diseases [52].

Gut microbe-complemented pathways for dietary carbohydrate metabolism

Carbohydrates are a major component of the human diet that includes starch (amylose and amylopectin) and disaccharides such as sucrose, lactose, and maltose. Human gut bacteria produce a vast panel of CAZymes to degrade components of dietary fiber into metabolisable monosaccharides and disaccharides [33]. The human genome encodes at best 17 digestive enzymes [53]; for ex. lactase, α-amylase, maltase, isomaltase and sucrose. It has been known that human enzymes can hydrolyze disaccharides (sucrose, lactose and maltose, etc.) and starch, but not other complex polysaccharides [54]. Hence, our ability to digest dietary plant carbohydrates resides entirely in our gut, where gut microbe-derived enzymes can hydrolyze complex dietary carbohydrates by producing a variety of CAZymes [55]. Thirteen gut bacterial enzymes predicted by ECemble were mapped in starch and sucrose metabolism pathway as shown in Figure 6. Enzymes responsible for the conversion of sucrose to glucose and bacterial degradation of pectin (a common component of dietary fibers) and xylan (polysaccharides in plant cell walls) are shown in the Table 6. These enzymes are predicted by our ECemble method from the gut microbial metagenomic data.

Figure 6
figure 6

Starch and sucrose metabolism pathway. The following color coding scheme is used for the pathways: Known human enzyme reactions (Light Red), Unknown human enzyme reactions that are predicted by ECemble (Pink), Known bacterial enzyme reactions in the gut microflora predicted by ECemble (Light green) and Unknown gut bacterial enzyme reactions that are predicted by ECemble (Light Blue).

Table 6 Enzymes and their function in starch and sucrose metabolism pathway.

Another GMC pathway, fructose and mannose metabolism, explains how bacterial enzymes complement human enzymes to metabolize dietary sugars (Figure 7). Fructose occurs as a free monosaccharide and an isomer of glucose. In Figure 7, predicted bacteria-encoded enzymes and known human-encoded enzymes are shown, where D-Fructose (fructose) is catalyzed by bacterial enzymes, Protein-N (pi) -phosphohistidine-sugar phosphotransferase (2.7.1.69) and Fructokinase (2.7.1.4) into D-Fructose-1 Phosphate and β-D-Fructose-6 Phosphate, respectively. β-D-Fructose-6P is metabolized to Glyceraldehyde-3P using human-encoded enzymes Phosphofructokinase (2.7.1.11) and Fructose-biphosphate aldolase (4.1.2.13). The Glyceraldehyde-3P compound is a part of the glycolysis (normal metabolism of sugars) pathway, which is the main energy generating mechanism in the body. This pathway demonstrates that the bacteria- and human-encoded enzymes complement and work in unison in the digestion and energy metabolism pathways.

Figure 7
figure 7

Fructose and mannose metabolism pathway. The following color coding scheme is used for the pathways: Known human enzyme reactions (Light Red), Unknown human enzyme reactions that are predicted by ECemble (Pink), Known bacterial enzyme reactions in the gut microflora predicted by ECemble (Light green) and Unknown gut bacterial enzyme reactions that are predicted by ECemble (Light Blue).

Analysis of enzyme profiles of obese and IBD subjects

A recent study published the metagenomic profiles of obese, lean and Inflammatory Bowel Disease (IBD) subjects [26]. This study also reported the translated gene products of gut microbiome from 124 metagenomic subject samples. We used these protein sequences to predict all the enzymes in each subject with our ECemble method, and analyzed the enzyme profiles of 'obese versus lean' (42 obese/82 lean) and 'IBD versus non-IBD' (25 IBD/99 Non-IBD) subjects. We identified 237 unique bacterial enzymes that are not encoded in human from the metagenomic samples of obesity, lean, IBD and non-IBD subjects. These include 222 known and 15 previously unknown enzymes in gut bacterial species. Details on how these enzyme reactions maps to KEGG human pathways are shown in Additional files 13 and 14. The taxonomic distribution of bacterial species from metagenomic samples is also shown in Additional file 15. The frequencies of enzymes present in the subjects were normalized based on the number of subjects in the obese/lean and IBD/non-IBD comparison groups and a Fisher's exact-test (P-value <0.05) using R [56] was conducted to determine the significantly enriched or depleted enzymes in the obese and IBD subjects (Additional files 16 and 17). Of the obesity-enriched enzymes, the most significant enzyme (P-value, 9.84E-04) is polygalacturonase (EC 3.2.1.15), which is encoded by Bacteroides and Prevotella species, and carries out pentose and glucuronate interconversions in starch and sucrose metabolism (Additional file 16). In contrary, urease (EC: 3.5.1.5) encoding bacteria are found in fewer number of obese subjects compared to lean subjects (obese/lean ratio = 0.88, P-value = 0.0179), suggesting that the loss or absence of this enzyme may be associated with obesity. In our analysis of IBD bacterial enzymes (Additional file 17), we found that non-IBD subjects predominantly host bacterial populations that contain an enzyme, Glucose-1-phosphate thymidylyltransferase (EC: 2.7.7.24, IBD/Non-IBD ratio= 0.06; P-value = 2.9E-12) compared to the IBD group. This enzyme is involved in the biosynthesis of L-rhamnose in bacteria. While there is no direct evidence to link lower levels of L-rhamnose to IBD, increased lactulose/L-rhamnose permeability ratio or decreased L-rhamnose in human intestinal permeability is found to be associated with IBD [57, 58]. Hence, the consequences of the absence of bacterial populations that produce L-rhamnose in the IBD patients is worth investigating by experimental studies. The distribution of significant (p-value <0.05) enzymes in obese/lean and IBD/Non-IBD groups is presented in Figure 8. The comprehensive list of enzymes in each category (obese/lean and IBD/Non-IBD) is given in Additional files 16 and 17.

Figure 8
figure 8

Number of Obesity and IBD enzymes in each of the major KEGG pathway categories. (A) Number of significant (p-value < 0.05; Fisher Exact Test) enzymes from obese vs. lean subjects in each major pathway categories (B) Number of significant (p-value < 0.05; Fisher Exact Test) enzymes from IBD vs. Non-IBD subjects in each major pathway categories.

Conclusions

A consensus-based ensemble method, ECemble, was implemented to identify enzymes from non-enzymes, and to hierarchically predict the class and subclasses of an enzyme up to L4 of the EC number. Comparison against BLAST and EFICAz methods showed superior performance of ECemble, both in coverage and accuracy. The superior coverage can be attributed to the generic protein domain feature space used in this method, while the improved accuracy resulted from the stringent consensus-based ensemble approach. Application of ECemble to predict full complements of enzymes from 10 sequenced genomes of model organisms has generated new annotations for unknown enzymes as well as full annotations for undercharacterized enzymes. Similarly, ECemble method enabled us to predict bacterial enzymes present in the gut microbiome and consequently use them to study the dependence of human metabolism on gut microbe-derived enzymes. Mapping of human and microbe-derived enzymes to KEGG metabolic pathways revealed that gut microbe-derived enzymes, especially those involved in the digestion of dietary nutrients, are essential components of a number of human pathways. Further application of this method to study the profiles of gut microbe-derived enzymes in lean versus obese, and IBD versus non-IBD subjects showed that certain enzymes were significantly enriched or depleted in these comparison groups, warranting further studies to understand the role of these enzymes on certain disease conditions. Two important merits of ECemble are that it can predict solely based on the protein sequence and also fully annotates enzymes by hierarchically assigning classes and subclasses up to L4. Hence, it can be a valuable tool for accurately annotating the entire enzyme complements of individual genomes as well as the mixed genomes from metagenomic studies. As evident from this study, ECemble can be effectively used to study the metabolic interactions between the host and microbes or those among the members of a community in a microflora. Being a generic method, it can be applied to study the systems level pathway interactions in any organism or microbial community.

Methods

ECemble dataset

The dataset include enzyme sequences taken from the EXPASY enzyme database [59], which is built based on SwissProt database [45] that contains reviewed and experimentally determined annotations on enzymes. We applied the following filters to obtain high-quality data for testing and training our program: duplicate sequences, sequences annotated as 'fragments' and those shorter than 50 residues in length were removed; sequences having more than one EC number were also removed. To avoid redundancy in the dataset, no more than 70% sequence identity was allowed among sequences within a class for each of the six EC classes, using the CD-HIT program [43]. Consequently, 64,950 non-redundant enzyme sequences (positive dataset) were collected that broadly cover all the known enzyme classes and subclasses. Similarly, a negative dataset (non-enzymes) of 128,475 sequences were collected from a SwissProt database by following the same filtering steps as mentioned above. A class-wise distribution of known enzyme sequences is listed in Table 7. Sequence to domain mapping was done for enzyme and non-enzyme sequences based on the Pfam [38] and Superfamily [39] databases using HMMSCAN [60], and the Prosite database using PSSCAN [61].

Table 7 Class-wise statistics on the number of enzyme sequences, and the subclasses at each level.

Feature databases

The most important features of enzymes that differentiate them from non-enzymes are their structure, function and catalytic sites. The feature set we used for machine learning in this study was based on the functional and structural domains, and sequence motifs. Structural domains define evolutionarily conserved region of proteins that can fold independently, while functional domains define evolutionarily conserved regions that can independently perform a specific biological function. Similarly, short sequence motifs define the catalytic or binding sites of enzymes. Machine learning algorithms exploit these different sets of features from enzyme and non-enzyme data in differentiating one from other.

Three popular domain databases - Pfam, Superfamily and Prosite - were used to extract the domain features of enzyme and non-enzyme sequences for building the ECemble method. Pfam (Protein Family) database contains functional domain information for protein sequences; Prosite database is a collection of biologically important sites, patterns and sequence motifs associated with protein functions; and the Superfamily database contains sets of homologous proteins that conserve structural and active site features. Together, these three databases provide a comprehensive coverage of functional, active site and structure-based features of enzyme and non-enzyme protein sequences.

Machine learning methods

We selected five diverse and most popular machine learning classifiers; Naïve Bayes [34], k-Nearest Neighbor (KNN) [35] classifier, Support Vector Machine (SVM) [36], Decision Stump (DS) [41] and Random Forest (tree-based) classifiers (RFC) [37], to build models for enzyme identification and classification. We used the WEKA 3.7.5 [42] package, which is an open-source, Java-based framework to build classification models using different ML techniques.

Ten proteome datasets

We downloaded the proteomes of ten model organisms from the SwissProt database [45], ranging from animal, plant, fungal and bacterial species such as Saccharomyces cerevisiae (yeast), Caenorhibditis elegans (C. elegans), Drosophila melanogaster (fruitfly), Danio rerio (zebrafish), Gallus gallus (chicken), Mus musculus (mouse), Homo sapiens (human), Escherichia coli (E. coli), Oryza sativa (rice) and Arabidopsis thaliana (Arabidopsis). We filtered out the experimentally known enzyme sequences from each of the ten proteomes and tested the unclassified sequences with ECemble.

Human gut metagenomic samples

We used a study on the human gut microbiome from the Beijing Genomics Institute's (BGI) Metagenomic Sequencing Project [26] that describes a set of about 3.3 million microbial genes sequenced and assembled from fecal samples of 124 individuals (Additional file 18 represent the patient profiles and health status), to study the impact of human gut microbiome on human metabolism. As part of the filtering steps, sequences mapped to viruses, archaea, eukaryota, unclassified and unknown sequences were removed. About 2.5 million translated protein sequences from the assembled scaffolds of bacterial genomes were predicted by our ECemble method to identify all the bacterial enzymes in the human gut microbiome. The diversity of taxonomic groups of bacteria at each level is given in Additional file 19: Figure S5.

KEGG database

We used Kyoto Encyclopedia of Genes and Genomes (KEGG) [46], the most comprehensive database source that integrates genomic, proteomic and systemic functional information for pathway analysis. The KEGG Pathway suite is a collection of manually drawn maps demonstrating the existing knowledge on the molecular interaction and reaction networks. KEGG pathways were used as reference pathways to map human and gut metagenomic enzymes for analysis of gut microbe-complemented human metabolic pathways.

Performance measurements

We used standard evaluation metrics that include 10-fold cross validation and ROC curves. In 10-fold cross-validation, sequences at each level are divided into ten parts, models are built using nine parts, and predictions are generated and evaluated on the data contained in one part. This procedure is repeated ten times, where each part is tested against the models built from nine other parts. The average performance of the ten models is considered as an unbiased estimate of the training model performance. After cross-validation, we assessed the performance of the fully trained classifier models using the test set (20% of original data) that were hidden from the classifiers. We report standard performance measures over each enzyme class including the following: true positives (TP) as the number of sequences that are correctly identified in a class that belongs to them; false negatives (FN) as the number of sequences that are not identified in a class that belongs to them; true negatives (TN) as the number of sequences that are not found in a class that does not belong to them; false positives (FP) as the number of sequences that are identified in a class that does not belong to them; sensitivity as the proportion of true positives that are predicted as positives; specificity is the proportion of true negatives that are predicted as negatives. The sensitivity and specificity are given by, sensitivity = TP/ (TP + FN); specificity = TN/ (TN + FP).

We also report accuracy in a class as the ratio of the number of correctly predicted enzyme sequences to the total number of sequences in that class. We optimize and validate the accuracy of ECemble by selecting the optimal model(s) that has maximum true positive rate (sensitivity) and minimum false positive rate (1-specificity). Finally, we show ROC (Receiver Operating Characteristic) curves as a graphical means of measuring the performance for each class at each level of the prediction process, and the area under the curve (AUC) as a numeric measure of performance depicted by ROC curves.

Abbreviations

AUC:

Area Under Curve

DS:

Decision Stump

EC:

Enzyme Commission

IBD:

Inflammatory Bowel Disease

KNN:

K-Nearest Neighbor

KEGG:

Kyoto Encyclopedia of Genes and Genomes

ML:

Machine Learning

NBC:

Naïve Bayes Classifier

RFC:

Random Forest Classifier

ROC:

Receiver Operating Characteristic

SF:

Superfamily

SVM:

Support Vector Machine

WEKA:

Waikato Environment for Knowledge Analysis

References

  1. Jeremy MB, John LT, Lubert S: Biochemistry. 1988, W.H. Freeman, 444-445.

    Google Scholar 

  2. Schmidt S: Metabolites: a helping hand for pathway evolution?. Trends Biochem Sci. 2003, 28: 336-341. 10.1016/S0968-0004(03)00114-2.

    Article  CAS  PubMed  Google Scholar 

  3. Shah I, Hunter L: Predicting enzyme function from sequence: a systematic appraisal. Proc International Conf Intell Syst Mol Biol- ISMBInternational Conf Intell Syst Mol Biol. 1997, 5: 276-283.

    CAS  Google Scholar 

  4. Webb EC: Enzyme Nomenclature. 1992, San Diego: Academic Press

    Google Scholar 

  5. Tian W, Skolnick J: How well is enzyme function conserved as a function of pairwise sequence identity?. J Mol Biol. 2003, 333: 863-882. 10.1016/j.jmb.2003.08.057.

    Article  CAS  PubMed  Google Scholar 

  6. Espadaler J, Eswar N, Querol E, Avilés FX, Sali A, Marti-Renom M a, Oliva B: Prediction of enzyme function by combining sequence similarity and protein interactions. BMC Bioinformatics. 2008, 16: 249-

    Article  Google Scholar 

  7. Cheng BYM, Carbonell JG, Klein-Seetharaman J: Protein classification based on text document classification techniques. Proteins. 2005, 58: 955-970. 10.1002/prot.20373.

    Article  CAS  PubMed  Google Scholar 

  8. McCallum a, Nigam K: A comparison of event models for naive bayes text classification. Work Learn Text Categ. 1998, 41-48.

    Google Scholar 

  9. Mitchell TM: Machine Learning Boston. 1997, 4: 432-

    Google Scholar 

  10. Duda R, Hart P, Stork D: Pattern Classification. 2001, 680-

    Google Scholar 

  11. Mohammed A, Guda C: Computational Approaches for Automated Classification of Enzyme Sequences. J Proteomics Bioinforma. 2011, 04: 147-152.

    Article  CAS  Google Scholar 

  12. Chou KC: Using amphiphilic pseudo amino acid composition to predict enzyme subfamily classes. Bioinformatics. 2005, 21: 10-19. 10.1093/bioinformatics/bth466.

    Article  CAS  PubMed  Google Scholar 

  13. Otto TD, Guimarães ACR, Degrave WM, de Miranda AB: AnEnPi: identification and annotation of analogous enzymes. BMC Bioinformatics. 2008, 16: 544-

    Article  Google Scholar 

  14. Röttig M, Rausch C, Kohlbacher O: Combining structure and sequence information allows automated prediction of substrate specificities within enzyme families. PLoS Comput Biol. 2010, 6: e1000636-10.1371/journal.pcbi.1000636.

    Article  PubMed Central  PubMed  Google Scholar 

  15. Almonacid DE, Babbitt PC: Toward mechanistic classification of enzyme functions. Curr Opin Chem Biol. 2011, 15: 435-442. 10.1016/j.cbpa.2011.03.008.

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  16. Lu L, Qian Z, Cai YD, Li Y: ECS: an automatic enzyme classifier based on functional domain composition. Comput Biol Chem. 2007, 31: 226-32. 10.1016/j.compbiolchem.2007.03.008.

    Article  CAS  PubMed  Google Scholar 

  17. Shen HB, Chou KC: EzyPred: a top-down approach for predicting enzyme functional classes and subclasses. Biochem Biophys Res Commun. 2007, 364: 53-9. 10.1016/j.bbrc.2007.09.098.

    Article  CAS  PubMed  Google Scholar 

  18. Weingart U, Lavi Y, Horn D: Data mining of enzymes using specific peptides. BMC Bioinformatics. 2009, 10: 446-10.1186/1471-2105-10-446.

    Article  PubMed Central  PubMed  Google Scholar 

  19. Nasibov E, Kandemir-Cavas C: Efficiency analysis of KNN and minimum distance-based classifiers in enzyme family prediction. Comput Biol Chem. 2009, 33: 461-4. 10.1016/j.compbiolchem.2009.09.002.

    Article  CAS  PubMed  Google Scholar 

  20. Shi R, Hu X: Predicting enzyme subclasses by using support vector machine with composite vectors. Protein Pept Lett. 2010, 17: 599-604. 10.2174/092986610791112710.

    Article  CAS  PubMed  Google Scholar 

  21. Wang YC, Wang Y, Yang ZX, Deng NY: Support vector machine prediction of enzyme function with conjoint triad feature and hierarchical context. BMC Syst Biol. 2011, 5 (Suppl 1): S6-10.1186/1752-0509-5-S1-S6.

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  22. Borro LC, Oliveira SRM, Yamagishi MEB, Mancini AL, Jardine JG, Mazoni I, Santos EH Dos, Higa RH, Kuser PR, Neshich G: Predicting enzyme class from protein structure using Bayesian classification. Genet Mol Res. 2006, 5: 193-202.

    CAS  PubMed  Google Scholar 

  23. Hung SS, Wasmuth J, Sanford C, Parkinson J: DETECT--a density estimation tool for enzyme classification and its application to Plasmodium falciparum. Bioinformatics. 2010, 26: 1690-8. 10.1093/bioinformatics/btq266.

    Article  CAS  PubMed  Google Scholar 

  24. Arakaki AK, Huang Y, Skolnick J: EFICAz2: enzyme function inference by a combined approach enhanced by machine learning. BMC Bioinformatics. 2009, 10: 107-10.1186/1471-2105-10-107.

    Article  PubMed Central  PubMed  Google Scholar 

  25. Syed U, Yona G: Enzyme function prediction with interpretable models. Methods Mol Biol. 2009, 541: 373-420.

    Article  CAS  PubMed  Google Scholar 

  26. Qin J, Li R, Raes J, Arumugam M, Burgdorf KS, Manichanh C, Nielsen T, Pons N, Levenez F, Yamada T, Mende DR, Li J, Xu J, Li S, Li D, Cao J, Wang B, Liang H, Zheng H, Xie Y, Tap J, Lepage P, Bertalan M, Batto JM, Hansen T, Le Paslier D, Linneberg A, Nielsen HB, Pelletier E, Renault P, et al: A human gut microbial gene catalogue established by metagenomic sequencing. Nature. 2010, 464: 59-65. 10.1038/nature08821.

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  27. Bäckhed F: Host responses to the human microbiome. Nutr Rev. 2012, 70 (Suppl 1): S14-7.

    Article  PubMed  Google Scholar 

  28. Cummings JH: Microbial digestion of complex carbohydrates in man. Proc Nutr Soc. 1984, 43: 35-44. 10.1079/PNS19840025.

    Article  CAS  PubMed  Google Scholar 

  29. Hooper L V, Midtvedt T, Gordon JI: How host-microbial interactions shape the nutrient environment of the mammalian intestine. Annu Rev Nutr. 2002, 22: 283-307. 10.1146/annurev.nutr.22.011602.092259.

    Article  PubMed  Google Scholar 

  30. Kelly D, Campbell JI, King TP, Grant G, Jansson EA, Coutts AGP, Pettersson S, Conway S: Commensal anaerobic gut bacteria attenuate inflammation by regulating nuclear-cytoplasmic shuttling of PPAR-gamma and RelA. Nat Immunol. 2004, 5: 104-112. 10.1038/ni1018.

    Article  CAS  PubMed  Google Scholar 

  31. Greenblum S, Turnbaugh PJ, Borenstein E: Metagenomic systems biology of the human gut microbiome reveals topological shifts associated with obesity and inflammatory bowel disease. Proc Natl Acad Sci USA. 2012, 109: 594-9. 10.1073/pnas.1116053109.

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  32. Nicholson JK, Holmes E, Kinross J, Burcelin R, Gibson G, Jia W, Pettersson S: Host-gut microbiota metabolic interactions. Science. 2012, 336: 1262-7. 10.1126/science.1223813.

    Article  CAS  PubMed  Google Scholar 

  33. Lepage P, Leclerc MC, Joossens M, Mondot S, Blottière HM, Raes J, Ehrlich D, Doré J: A metagenomic insight into our gut's microbiome. Gut. 2013, 62: 146-58. 10.1136/gutjnl-2011-301805.

    Article  PubMed  Google Scholar 

  34. Rish I: An empirical study of the naive Bayes classifier. IJCAI 2001 Work Empir Methods Artif Intell. 2001, 3: 41-46.

    Google Scholar 

  35. Cover T, Hart P: Nearest neighbor pattern classification. IEEE Trans Inf Theory. 1967, 13: 21-27.

    Article  Google Scholar 

  36. Cortes C, Vapnik V: Support-vector networks. Mach Learn. 1995, 20: 273-297.

    Google Scholar 

  37. Breiman L: Random Forests. Mach Learn. 2001, 45: 5-32. 10.1023/A:1010933404324.

    Article  Google Scholar 

  38. Bateman A, Birney E, Durbin R, Eddy SR, Finn RD, Sonnhammer EL: Pfam 3.1: 1313 multiple alignments and profile HMMs match the majority of proteins. Nucleic Acids Res. 1999, 27: 260-262. 10.1093/nar/27.1.260.

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  39. Wilson D, Madera M, Vogel C, Chothia C, Gough J: The SUPERFAMILY database in 2007: families and functions. Nucleic Acids Res. 2007, 35: D308-D313. 10.1093/nar/gkl910.

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  40. Bairoch A, Bucher P, Hofmann K: The PROSITE database, its status in 1997. Nucleic Acids Res. 1997, 25: 217-221. 10.1093/nar/25.1.217.

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  41. Iba W, Langley P: Induction of One-Level Decision Trees. ML92 Proc Ninth Int Conf Mach Learn Aberdeen, Scotland, 1-3 July 1992. 1992, San Francisco, CA: Morgan Kaufmann, 233-240.

    Google Scholar 

  42. Hall M, Frank E, Holmes G, Pfahringer B, Reutemann P, Witten IH: The WEKA data mining software: an update. SIGKDD Explor. 2009, 11: 10-18. 10.1145/1656274.1656278.

    Article  Google Scholar 

  43. Li W, Godzik A: Cd-hit: a fast program for clustering and comparing large sets of protein or nucleotide sequences. Bioinformatics. 2006, 22: 1658-1659. 10.1093/bioinformatics/btl158.

    Article  CAS  PubMed  Google Scholar 

  44. Magrane M, Consortium U: UniProt Knowledgebase: a hub of integrated protein data. Database (Oxford). 2011, 2011: bar009-

    Article  Google Scholar 

  45. Bairoch A, Apweiler R: The SWISS-PROT protein sequence database and its supplement TrEMBL in 2000. Nucleic Acids Res. 2000, 28: 45-48. 10.1093/nar/28.1.45.

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  46. Kanehisa M, Goto S: KEGG: kyoto encyclopedia of genes and genomes. Nucleic Acids Res. 2000, 28: 27-30. 10.1093/nar/28.1.27.

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  47. Abratt VR, Reid SJ: Oxalate-degrading bacteria of the human gut as probiotics in the management of kidney stone disease. Adv Appl Microbiol. 2010, 72: 63-87.

    Article  CAS  PubMed  Google Scholar 

  48. Blaut M, Clavel T: Metabolic diversity of the intestinal microbiota: implications for health and disease. J Nutr. 2007, 137: 751S-5S.

    CAS  PubMed  Google Scholar 

  49. Dumas ME, Barton RH, Toye A, Cloarec O, Blancher C, Rothwell A, Fearnside J, Tatoud R, Blanc V, Lindon JC, Mitchell SC, Holmes E, McCarthy MI, Scott J, Gauguier D, Nicholson JK: Metabolic profiling reveals a contribution of gut microbiota to fatty liver phenotype in insulin-resistant mice. Proc Natl Acad Sci USA. 2006, 103: 12511-6. 10.1073/pnas.0601056103.

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  50. Clayton TA, Baker D, Lindon JC, Everett JR, Nicholson JK: Pharmacometabonomic identification of a significant host-microbiome metabolic interaction affecting human drug metabolism. Proc Natl Acad Sci USA. 2009, 106: 14728-14733. 10.1073/pnas.0904489106.

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  51. Sousa T, Paterson R, Moore V, Carlsson A, Abrahamsson B, Basit AW: The gastrointestinal microbiota as a site for the biotransformation of drugs. Int J Pharm. 2008, 1-25.

    Google Scholar 

  52. Cani PD, Amar J, Iglesias MA, Poggi M, Knauf C, Bastelica D, Neyrinck AM, Fava F, Tuohy KM, Chabo C, Waget A, Delmée E, Cousin B, Sulpice T, Chamontin B, Ferrières J, Tanti JF, Gibson GR, Casteilla L, Delzenne NM, Alessi MC, Burcelin R: Metabolic endotoxemia initiates obesity and insulin resistance. Diabetes. 2007, 56: 1761-1772. 10.2337/db06-1491.

    Article  CAS  PubMed  Google Scholar 

  53. Turnbaugh PJ, Henrissat B, Gordon JI: Viewing the human microbiome through three-dimensional glasses: integrating structural and functional studies to better define the properties of myriad carbohydrate-active enzymes. Acta Crystallogr Sect F Struct Biol Cryst Commun. 2010, 66 (Pt 10): 1261-1264.

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  54. Devaraj S, Hemarajata P, Versalovic J: The human gut microbiome and body metabolism: implications for obesity and diabetes. Clin Chem. 2013, 59: 617-28. 10.1373/clinchem.2012.187617.

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  55. Cantarel BL, Lombard V, Henrissat B: Complex carbohydrate utilization by the healthy human microbiome. PLoS One. 2012, 7: e28742-10.1371/journal.pone.0028742.

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  56. R Development Core Team: R: A Language and Environment for Statistical Computing. R Found Stat Comput Vienna Austria. 2013, {ISBN} 3-900051-07-0

    Google Scholar 

  57. Sanderson IR, Boulton P, Menzies I, W-SJ : Improvement of abnormal lactulose / rhamnose permeability in active Crohn 's disease of the small bowel by an elemental diet. Gut. 1987, 28 (February): 1073-1076.

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  58. Desai D, Faubion WA, Sandborn WJ: Review article: biological activity markers in inflammatory bowel disease. Aliment Pharmacol Ther. 2007, 25: 247-255.

    Article  CAS  PubMed  Google Scholar 

  59. Bairoch A: The ENZYME database in 2000. Nucleic Acids Res. 2000, 28: 304-305. 10.1093/nar/28.1.304.

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  60. Johnson LS, Eddy SR, Portugaly E: Hidden Markov model speed heuristic and iterative HMM search procedure. BMC Bioinformatics. 2010, 11: 431-10.1186/1471-2105-11-431.

    Article  PubMed Central  PubMed  Google Scholar 

  61. de Castro E, Sigrist CJA, Gattiker A, Bulliard V, Langendijk-Genevaux PS, Gasteiger E, Bairoch A, Hulo N: ScanProsite: detection of PROSITE signature matches and ProRule-associated functional and structural residues in proteins. Nucleic Acids Res. 2006, 34: W362-W365. 10.1093/nar/gkl124.

    Article  PubMed Central  CAS  PubMed  Google Scholar 

Download references

Acknowledgements

This research has been partly supported by National Institutes of Health [1R01GM086533-01A1 to CG], startup funds to CG from UNMC, and the graduate fellowship from UNMC to AM. The publication costs for this article have been funded by the startup funds to CG from UNMC.

This article has been published as part of BMC Genomics Volume 16 Supplement 7, 2015: Selected articles from The International Conference on Intelligent Biology and Medicine (ICIBM) 2014: Genomics. The full contents of the supplement are available online at http://www.biomedcentral.com/bmcgenomics/supplements/16/S7.

Author information

Authors and Affiliations

Authors

Corresponding author

Correspondence to Chittibabu Guda.

Additional information

Competing interests

The authors declare that they have no competing interests.

Authors' contributions

AM developed the ECemble method using machine-learning approaches and applied it to analyze gut metabolic pathways, and drafted the manuscript. CG conceived the idea of the project, provided overall guidance and coordination, suggested improvement and corrected the manuscript. All authors read and approved the final version of the manuscript.

Electronic supplementary material

12864_2015_7210_MOESM1_ESM.pdf

Additional file 1: Figure S1. 10-fold cross validation and testing accuracies. A) For enzyme identification at EC Level-0 using ML classifiers Decision Stump (DS), Naïve Bayes Classifier (NBC), K-Nearest Neighbor (KNN), Support Vector Machine (SVM), and Random Forest Classifier (RFC). B) For enzyme classification at EC Level-1 using ML classifiers Decision Stump (DS), Naïve Bayes Classifier (NBC), K-Nearest Neighbor (KNN), Support vector Machine (SVM), and Random Forest Classifier (RFC). (PDF 335 KB)

Additional file 2: Table S1. Overall predictions of 2 classifiers vs. 3 classifiers. (PDF 174 KB)

12864_2015_7210_MOESM3_ESM.pdf

Additional file 3: Figure S2. Accuracy and distribution of enzyme sequence and class. A) Distribution of enzyme sequence and class coverage for cdh70, cdh60, cdh50 and cdh40 datasets. B) Accuracy at each EC level for cdh70, cdh60, cdh50 and cdh40 datasets. (PDF 135 KB)

Additional file 4: Table S2. Accuracy of ECemble method using cdh70, cdh60, cdh50 and cdh40 datasets (PDF 204 KB)

12864_2015_7210_MOESM5_ESM.xlsx

Additional file 5: Statistics for Known and ECemble predicted enzyme sequences from ten proteomes. The proteomes contains both reviewed and unreviewed sequences from UniProt (Data Sheet 'ReviewedAndUnreviewed'). The ECemble predictions for the reviewed sequences from 10 proteomes are given in Data Sheet 'PredictionsReviewedSet'. (XLSX 51 KB)

Additional file 6: ECemble predicted gut microbial enzymes. (XLSX 24 KB)

Additional file 7: Known human enzyme reactions from SwissProt and KEGG. (XLSX 56 KB)

Additional file 8: Figure S3. GMC pathways that validate the role of gut bacteria in human metabolism (PDF 1 MB)

12864_2015_7210_MOESM9_ESM.pdf

Additional file 9: Figure S4. GMC pathways that reveal the newly discovered role of gut bacteria in human metabolism (PDF 493 KB)

Additional file 10: Mapped human pathways in major pathway categories. (XLSX 11 KB)

Additional file 11: CAZymes found in metagenomic samples using ECemble. (XLSX 10 KB)

Additional file 12: Known enzymes common in both human and bacteria. (XLSX 30 KB)

Additional file 13: Known bacterial enzyme reactions in the gut microflora predicted by ECemble (XLSX 17 KB)

Additional file 14: Unknown gut bacterial enzyme reactions that are predicted by ECemble. (XLSX 13 KB)

Additional file 15: ECemble predicted gut bacterial enzymes with subject information. (XLSX 5 MB)

12864_2015_7210_MOESM16_ESM.xlsx

Additional file 16: Obesity/Lean enzymes from metagenomic samples using Fisher-exact test. Obese_Normalized is calculated as the ratio of number of obese subjects having a particular enzyme to the total number of obese subjects (42). Similarly, Lean_Normalized is calculated as the ratio of number of lean subjects having a particular enzyme to the total number of lean subjects (82). (XLSX 19 KB)

12864_2015_7210_MOESM17_ESM.xlsx

Additional file 17: IBD/Non-IBD enzymes from metagenomic samples using Fisher-exact test. IBD_Normalized is calculated as the ratio of number of IBD subjects having a particular enzyme to the total number of IBD subjects (25). Similarly, Non-IBD_Normalized is calculated as the ratio of number of Non-IBD subjects having a particular enzyme to the total number of Non-IBD subjects (99). (XLSX 22 KB)

Additional file 18: Metagenomic samples profile. (XLSX 13 KB)

Additional file 19: Figure S5. Taxonomic distribution of bacterial species from metagenomic samples. (PDF 159 KB)

Rights and permissions

This article is published under license to BioMed Central Ltd. This is an Open Access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/4.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited. 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.

Reprints and permissions

About this article

Check for updates. Verify currency and authenticity via CrossMark

Cite this article

Mohammed, A., Guda, C. Application of a hierarchical enzyme classification method reveals the role of gut microbiome in human metabolism. BMC Genomics 16 (Suppl 7), S16 (2015). https://doi.org/10.1186/1471-2164-16-S7-S16

Download citation

  • Published:

  • DOI: https://doi.org/10.1186/1471-2164-16-S7-S16

Keywords