Phenotypic differentiation of gastrointestinal microbes is reflected in their encoded metabolic repertoires
© Bauer et al. 2016
Received: 18 May 2015
Accepted: 30 September 2015
Published: 30 November 2015
The human gastrointestinal tract harbors a diverse microbial community, in which metabolic phenotypes play important roles for the human host. Recent developments in meta-omics attempt to unravel metabolic roles of microbes by linking genotypic and phenotypic characteristics. This connection, however, still remains poorly understood with respect to its evolutionary and ecological context.
We generated automatically refined draft genome-scale metabolic models of 301 representative intestinal microbes in silico. We applied a combination of unsupervised machine-learning and systems biology techniques to study individual and global differences in genomic content and inferred metabolic capabilities. Based on the global metabolic differences, we found that energy metabolism and membrane synthesis play important roles in delineating different taxonomic groups. Furthermore, we found an exponential relationship between phylogeny and the reaction composition, meaning that closely related microbes of the same genus can exhibit pronounced differences with respect to their metabolic capabilities while at the family level only marginal metabolic differences can be observed. This finding was further substantiated by the metabolic divergence within different genera. In particular, we could distinguish three sub-type clusters based on membrane and energy metabolism within the Lactobacilli as well as two clusters within the Bifidobacteria and Bacteroides.
We demonstrate that phenotypic differentiation within closely related species could be explained by their metabolic repertoire rather than their phylogenetic relationships. These results have important implications in our understanding of the ecological and evolutionary complexity of the human gastrointestinal microbiome.
KeywordsGenome-scale metabolic reconstructions Metabolic modeling Metagenomics Intestinal microbiota Evolution Ecology
Recent advances in sequencing technologies have greatly improved our knowledge about the metabolic complexity of the human microbiome and provide novel approaches to identify beneficial microbes . In particular, sequencing the (ideally) entire genomic content (i.e., metagenomic sequencing) of the intestinal microbiota has allowed the establishment of a catalog of main groups of microorganisms present in the gastrointestinal tract and potential metabolic pathways  by avoiding culturing and isolation of individual microbial organisms. In this respect, endeavors of the human microbiome project  and the MetaHIT consortium  aim at establishing comprehensive data-sets of metagenomic content, metabolic functions, and taxonomic compositions within human individuals as well as the isolation and sequencing of numerous microbial taxa.
Despite these efforts, however, we are still lacking a comprehensive mechanistic understanding of the intestinal microbiota. One major hurdle in achieving this goal is the lack of organismal system boundaries, enabling us to associate the presence of metabolic pathways in the microbiome with a specific bacterium. Inferring metabolic roles by taxonomic classification alone is difficult because phylogenetically closely related organisms might be very different in their metabolism . It may be therefore challenging to associate functional roles to entire taxonomic groups  to conjecture the biological relevance of intestinal bacteria. For instance, members of the same genus, or even of the same species, can be both probiotic and pathogenic , indicating a differential strain-specific adaptation. In this context, nutrient utilization can be a strong determinant for the adaptation to varying environments, since it can give a competitive advantage to other organisms that are metabolically less versatile. Thus, having additional metabolic functions can aid microbes in occupying further niches within the human gut. Accordingly, the functional consequences for the host change.
Current developments in systems biology allow the modeling of microbial metabolism to gain a mechanistic insight into the relationship between genotype and phenotype . Genome-scale metabolic reconstructions form the basis of such modeling efforts. A reconstruction is assembled based on the genomic sequences as well as biochemical and phenotypic data of a target organism, and accounts for metabolic genes, enzymes, and their associated reactions [10, 11]. Genome-scale metabolic reconstructions serve as a blueprint for condition-specific metabolic models [10, 11], which are obtained by the application of constraints, such as known nutrient uptake rates. The reconstruction process often includes a gap-filling procedure [12, 13], in which additional reactions are included to better model biologically relevant phenotypes, such as the formation of all known biomass precursors . Metabolic models can be studied using a variety of mathematical methods . One frequently used approach is flux balance analysis, which is applied to investigate a functional steady-state flux distribution of the modeled system, while maximizing (or minimizing) a particular cellular objective (e.g., production of biomass precursors) . This modeling approach has been used to investigate nutrient requirements , gene essentialities , and metabolic interactions  for organisms of interest, thereby providing new insights into phenotypic and metabolic properties. The reconstruction process relies on the availability of detailed phenotypic data for the target organism , which is usually not available for many of the commonly found microbes in the human gut [1, 3]. To obtain representative metabolic reconstructions for these less well-studied organisms, automatic tools have been developed in recent years, such as the Model SEED platform , to provide a valuable starting point for metabolic modeling. In fact, draft reconstructions have been used to generate hypotheses about the target organisms with subsequent experimental validation, leading to the refinement of the metabolic reconstruction [14, 20–22].
In this study, we generated automatically refined draft genome-scale metabolic models of 301 representative intestinal microbes in silico based on whole genome sequences of the human microbiome project using an established approach . We applied a combination of unsupervised machine-learning and computational modeling techniques to study individual and global differences of the metabolic models and the original genomes. Our key results include: i) divergent reactions involved in energy metabolism and membrane synthesis which are most relevant to discriminate different phylogenetic groups, ii) a linear relationship between differences in metabolic reaction potential and essential nutrients determined by flux balance analysis which indicates that the phenotype is directly correlated to the metabolic repertoire, iii) differences in metabolic reaction potential and phylogeny which exhibit an exponential relationship, suggesting an explanation as to why closely related microbes can be very different in their metabolic traits while at less-resolved phylogenetic distances only marginal differences in metabolic diversity can be observed, iv) local differences in pathway presence which can be used to further distinguish representatives of Lactobacillus, Bifidobacteria, and Bacteroides. In summary, we demonstrate the importance of the metabolic repertoire of microbes to predict their phenotypic behavior in an ecological and evolutionary context.
Results and discussion
Selected microbes as a model for the human gut microbiota
Genome and metabolic model statistics of the selected microbes
Average per taxonomic group
Our analysis also included microbes with draft genomes (Additional file 3: Table S1), requiring the assessment of the overall genome completeness and the potential impact on gene annotations and consequently on the generated metabolic reconstructions. The completeness and possible genomic contamination by other microorganisms of the individual 301 of the individual genomes was assessed using a collection of 107 universal, single-copy genes [28, 29]. In our set of 301 genomes, the average estimated genome completeness was 95 % (Table 1). We further investigated the genome size and annotated genes among the 301 organisms (Table 1). Gammaproteobacteria had generally large genomes and a high number of annotated genes, while members of the order Bacteroidia had in general larger genomes but a lower number of annotated genes. This difference could be attributed to differences in annotation efficiencies, as Proteobacteria (in particular, gut specific Escherichia species) are very well-studied and thus have more homologous genes. Consequently, the number of reactions in the constructed metabolic models was higher and the number of reactions added via gap-filling lower. In contrast, we found a higher number of gap-filled reactions and a lower number of reactions in Actinobacteria (Table 1). This bias, which is well established for metagenomic analyses , is most likely the result of having less experimental data and validated gene annotation available for Actinobacteria. The presence of this apparent annotation bias underlines the limitation in current annotation techniques affecting particularly phylogenetically distant microbes [29–31] and highlights the need for more detailed experimental biochemical studies to elucidate gene functions in phyla distant to those containing model organisms .
Global reaction differences recapitulate conserved taxonomic patterns and phenotypes
The separation between Gammaproteobacteria and Actinobacteria highlights that our reconstructions captured taxa-specific metabolic features, despite the mentioned annotation bias. Furthermore, Clostridia species showed a high metabolic diversity and overlapped with clusters of other microbial taxa (Fig. 2), which is consistent with the reported metabolic variety of these bacteria and their corresponding beneficial traits in the human gut . Erysipelotrichia representatives are closely but nonetheless distinctly placed relative to the Clostridia in the 2D principle coordinate plot (Fig. 2). Intriguingly, members of Erysipelotrichia were formerly considered as Clostridia based on the phylogeny of marker genes  but then re-assigned to a novel class based on their phylogeny and membrane composition . Similar to the Clostridia, Bacilli species were also widely spread in the 2D principle coordinate plot (Fig. 2), reflecting their metabolic versatility . In contrast, other taxa had more dense clusters, particularly Actinobacteria, reflecting more specialized roles of these bacteria, such as the conversion of polysaccharides .
Overall, we propose that metabolic reconstructions could be used, in addition to canonical approaches, to assist in the taxonomic definition of novel microbes and the re-assignment of already described microbes into better defined taxonomic groups. In particular, our approach has the advantage of considering functional characteristics, in contrast to methods solely relying on the presence and phylogeny of marker genes. As also pointed out by previous studies , functional repertoires can have a positive influence on the annotation quality of taxonomic groups. Ultimately, this could shed light onto the metabolic versatility of microbes in general or in specific habitats, such as the human gut.
Energy and membrane metabolism as markers for metabolic divergence
Following the broader characterization, we aimed to obtain a better understanding of the reactions driving the observed separation in the first two coordinates. The separation of taxonomic groups is due to reactions involved in membrane synthesis and central metabolism (Fig. 2). In particular, different types of lysophosolipase reactions exhibit the highest explanatory power (Additional file 4: Table S3). These reactions convert various phospholipid precursors (differing in their number of C-atoms) and have the same direction in the first principle coordinate, because all reactions can be carried out by single enzymes and are thus linearly dependent. Similarly, the amylomaltases catalyze multiple reactions differing in their substrates (Additional file 4: Table S3). For the enoyl-ACP reductase, we found a variety of reactions with different directions toward the first principle coordinate (Fig. 2). This variation in angle represents a potential variation in distinct yet convergent fatty acid synthesis processes involved in energy metabolism and known to be present in the human gut microbiome , thus contributing to the discrimination of the different types of bacteria. This observation is consistent with the fact that fatty acid profiles have been used to characterize microbial communities before the advent of nucleic acid-based methods . The synthesis of endotoxins was positively associated with the distribution of Gram-negative pathogenic species within the Proteobacteria and Fusobacteria, which is in accordance with previously reported correlations between various diseases and the abundance of Proteobacteria-producing endotoxins . The transport and utilization of diverse carbohydrates involved in energy metabolism, such as mannitol, mannose, and fructose, were positively associated with the location of the Bacilli cluster in the 2D principle coordinate plot (Fig. 2). This association highlights the variety of substrate consumption as represented by these reconstructions of microbial metabolism. In accordance with the literature, Bacilli are known to utilize a broad range of carbohydrates .
The differentiation of taxonomic groups based on reactions involved in energy and membrane metabolism may have important implications in understanding the evolution and heterogeneity of intestinal microbes. For instance, Gram-negative bacteria have been reported to change their membrane composition  in order to cope with environmental influences, such as antibiotics and human immune agents, many of which target bacterial membrane compounds . Additionally, ecological changes within the microbial community  can provoke a differentiation in metabolic capabilities involved in energy metabolism leading to altered interactions of the community with the human host, supporting the observed high explanatory power of metabolic reactions toward cluster separation.
The relationship between genotype, phenotype, and metabolic repertoire is non-linear
The relationship between phylogeny, metabolic repertoire, and phenotype is taxon-dependent
Summary statistics of the relationship between reaction content, phylogeny, and essential nutrients
Reaction differences reflect metabolic versatility among closely related microbes
To further investigate the metabolic divergence within closely related microbes of the same taxonomic group, we used t-distributed stochastic neighbor embedding (t-SNE)  for the two-dimensional visualization of the reaction similarities (Fig. 5, see Additional file 8: Figure S5 for point labels). t-SNE is a non-linear, non-parametric dimensionality reduction and has been used previously to reveal data-inherent cluster structures [59–61]. This method enabled us to identify fine-scale reaction differences, in addition to the principal coordinate analysis (Fig. 2). Several distinct clusters were apparent and corresponded to the different bacterial classes (Fig. 5). We further focused our analysis on the three most abundant genera in our model selection (Table 1). For Lactobacillus, we noted a widespread metabolic repertoire and thus a relatively large variability of members within this group (Fig. 5). We identified three distinct subclusters (La1, La2, and La3) within this genus. While La2 showed major overlaps with the other two clusters, La1 and La3 were distinct from each other. We investigated the differences in the reaction sets between the representatives of the different subclusters (Additional file 9: Table S4). Based on the present reaction sets, La1 corresponds to obligate homofermentative La2 to facultative homofermentative and La3 to obligate heterofermentative pathways involved in the energy metabolism of lactic acid bacteria (Fig. 5b). The pathway presence in the genomes explains why La2 overlaps with the other clusters, since the facultative homofermentative group (La2) shares reactions with the obligate homofermentative (La1) and heterofermentative group (La3) . In agreement with the literature, these subclusters correspond to known divergent pathways involved in energy metabolism in Lactobacilli . This distinction of biologically relevant phenotypic groups using predicted difference in metabolic reactions encouraged us to propose novel bacterial sub-types. Therefore, we confirmed for our choice of the number of subclusters by performing hierarchical clustering (Fig. 3) to ensure that the subclusters were substantially different. For the Bifidobacteria, we propose two distinct subclusters (Bi1 and Bi2), which differed in the reactions involved in energy metabolism and membrane biosynthesis (Fig. 5c). For the energy metabolism, numerous reactions involved in the uptake and utilization of diverse carbohydrates were observed for members of the subcluster Bi1 (Additional file 11: Table S5), corresponding to known strain-specific differences within closely related Bifidobacteria . Furthermore, we found reactions involved in the uptake and conversion of glucosamine to peptidoglycan, which could be associated with membrane composition in these two groups. To our knowledge, such pathway differentiation has not yet been proposed for Bifidobacteria. For the Bacteroidia, we could distinguish two subclusters (Ba1 and Ba2). The differences between these clusters can be attributed to the membrane biosynthesis (Fig. 5d; Additional file 12: Table S6). Members of Ba2 possess various pathway types leading to the production of varying phosphatidylglycerol compounds, whereas members of Ba1 can further process phosphatidylglycerol to myristic acid. This finding is of particular biological importance, when considering the virulence and signaling purposes of membrane lipids in Bacteroides species found in previous studies [63, 64], which links the phenotype to the synthesis of membrane compounds. Furthermore, since energy metabolism and substrate availability via the diet are major ecological driving forces within the human gut microbiota , the metabolic diversification of other closely related microbes, such as Lactobacillus spp. and Bifidobacterium spp., can be a necessary requirement to maintain a stable coexistence with each other and the host. Considering that optimal conditions for metabolic cooperation are dependent on the similarity between the metabolic repertoires of several species , this pathway analysis approach could be used to estimate cooperative as well as competitive strategies. In particular, microorganisms tend to have a higher cooperativity if they are not too similar nor too different , indicating that members of the same taxon, but different subclusters (Fig. 5b) might be able to co-exist, whereas functionally similar microbes may be more likely to compete with each other .
The requirement for a certain functional diversity to ensure a well-functioning cooperative intestinal microbiota is crucial to break down various complex dietary compounds and divide metabolic tasks among different community members . Our results complement these ideas by investigating the metabolic divergence within a model microbiota, which can be primarily distinguished by reactions involved in energy and membrane metabolism. These capabilities play important roles in shaping the interface between host and symbionts, and thereby may lead to a deeper understanding in addition to metagenomic analyses in which all microbial functions are assessed . Furthermore, the metabolic repertoire of microbes is proportional to their phenotypic properties, highlighting the importance of diversity in explaining the metabolic processes taking place within the human gut. In contrast to these properties, the metabolic repertoire exhibited an exponential relationship with phylogeny, underlining the challenges in inferring metabolic functions from phylogeny alone, in particular when using single gene-centric approaches such as via 16S rRNA gene amplicon sequencing. Moreover, this circumstance can be regarded as an important evolutionary and ecological feature of the microbiome; functional components constituting whole pathways can be very different within closely related species, whereas the metabolism in the overall metabolic repertoire is limited. In other words, by dividing the metabolic tasks between certain taxonomic groups, the microbiota can make efficient use out of a small set of functions thereby facilitating niche partitioning. This result has important implications when considering the overall species richness of the human gut microbiome in the context of different patients and diseases . Further analyses could prove these concepts by modeling interactions within bacteria and the use of the here reconstructed and refined genome-scale metabolic models.
Metabolic model selection, construction, and refinement
We selected a set of 301 microbes (Additional file 3: Table S1) representing species present in the normal gut microbiota of healthy individuals, according to previous studies [1, 3]. We retrieved the genome sequences as well as additional information about the sequencing status, oxygen requirement, taxonomic placement, and phenotype from the integrated microbial genome database . The completeness and possible genomic contamination by other microorganisms of the individual 301 genomes was assessed using a collection of 107 universal, single-copy genes . The genomic sequences were uploaded for gene annotation to the RAST server  using default parameters. Draft metabolic reconstructions were then built with these genome annotations using the Model SEED pipeline . To ensure, that the metabolic models are able to grow under anaerobic conditions, which are prevalent in their natural ecosystem, we modified, if necessary, one to five reactions to enable anaerobic growth. The reactions modified for each model are listed in Additional file 1: Table S2. For descriptive purposes, reactions in the metabolic models were translated into our in-house metabolite and reaction database. The original SEED reaction nomenclature was maintained for the growth simulation. All refined draft metabolic models are publically available in their Matlab format at http://thielelab.uni.lu/in_silico_models.
With v b as the flux through the biomass objective function, v as the vector of all reaction fluxes, v i,min as the minimal flux capacity of reaction i, and v i,max as the maximal flux capacity of reaction i. The solution (metabolic fluxes of all reactions) of this optimization problem can be obtained using linear programming. The flux through the biomass reactions can be interpreted as the growth rate of the microbe model. By setting the constraints v i,min and v i,max of exchange (transport) reactions, varying growth conditions can be simulated. Throughout this study, the maximal uptake was constrained to 10 mmol/gDW/h to estimate natural occurring conditions. The maximal achievable growth rate was calculated under these conditions by assuming that all exchange reactions are potentially active (equivalent to rich medium condition). Additionally, the absence of a particular metabolite in the medium was simulated by setting its minimal and maximal exchange reaction constraints (v i,min and v i,max ) to 0 mmol/gDW/h. By the iterative removal of each metabolite individually from the rich medium for each microbe model, different growth conditions were simulated. Essential nutrients were defined by growth rates smaller than 0.05 h−1 after removal from the medium. This cutoff was based on the estimated growth rate of microbes within the mammalian gut ; however, all calculated smaller growth rates were below 0.0001 h−1 and thus negligible.
Data mining of metabolic and genomic information
to calculate the metabolic proximity according to . Based on the obtained distance matrix of the reaction content, we used principle coordinate analysis  and t-SNE  for reducing the dimensionality from 301 to 2. The two-dimensional embeddings were visualized by scatter plots. Using principle coordinate analysis, we analyzed reaction differences between the metabolic models on a global scale by correlating each reaction to the principle coordinates and subsequently selecting the 200 reactions with the highest correlation (Additional file 4: Table S3). The t-SNE-based visualization was used to identify local differences, with a detailed analysis of cluster structures within the genera Lactobacillus, Bifidobacteria, and Bacteroides. The reaction set differences between the determined sub-types of these genera were then used to identify type specific pathways.
In addition to the determined metabolic difference, we used the phylogenetic relationships between the microbes as a measure of divergence. The phylogeny was computed with PhyloPhlAn, which uses a set of around 400 protein-coding genes for the phylogenetic placement . In addition to the 301 bacterial genomes, the genomes of the archaea Methanobrevibacter smithii ATCC35061 and Methanosphaera stadtmanae DSM 3091 were used as an out-group to root the phylogenetic tree (Additional file 5: Figure S2). The resulting phylogenetic tree was visualized using EvolView . The phylogenetic difference between the different bacteria was computed using the cophenetic distance based on the rooted tree .
Correlation between phylogeny, metabolic repertoire, and essential nutrients
to be the best fit. We complemented the exponential model with the Spearman correlation and the linear model with the Pearson correlation as a measure of association between the variables. The goodness of fit measures for the different models and subsets of the data can be found in Table 2. The fitted parameters α and β for all plots in Fig. 4 can be found in Additional file 10: Table S7.
The authors are thankful to Mrs Almut Heinken for helping with the refinement of the draft metabolic models to account for anaerobic metabolisms and Dr. Dmitry Ravcheev for providing information about aerobic and anaerobic metabolisms. This study was supported by the ATTRACT program grants (FNR/A12/01 and FNR/A09/03), the Aides à la Formation-Recherche (FNR/6783162, FNR /4964712) grant from the Luxembourg National Research Fund (FNR), and a European Union Joint Programming in Neurodegenerative Diseases grant (INTER/JPND/12/01).
Open AccessThis article is distributed under the terms of the Creative Commons Attribution 4.0 International License (http://creativecommons.org/licenses/by/4.0/), which permits unrestricted use, distribution, and reproduction in any medium, provided you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated.
- Qin J, Li Y, Cai Z, Li S, Zhu J, Zhang F, et al. A metagenome-wide association study of gut microbiota in type 2 diabetes. Nature. 2012;490(7418):55–60.PubMedView ArticleGoogle Scholar
- Flint HJ, Duncan SH, Scott KP, Louis P (2015). Links between diet, gut microbiota composition and gut metabolism. Proceedings of the Nutrition Society, 74, pp 13-22.doi:10.1017/S0029665114001463.
- Qin J, Li R, Raes J, Arumugam M, Burgdorf KS, Manichanh C, et al. A human gut microbial gene catalogue established by metagenomic sequencing. Nature. 2010;464(7285):59–65.PubMedPubMed CentralView ArticleGoogle Scholar
- Turnbaugh PJ, Ley RE, Hamady M, Fraser-Liggett C, Knight R, Gordon JI. The human microbiome project: exploring the microbial part of ourselves in a changing world. Nature. 2007;449(7164):804.PubMedPubMed CentralView ArticleGoogle Scholar
- Ehrlich SD. MetaHIT: The European Union Project on metagenomics of the human intestinal tract. Metagenomics of the Human Body. New York: Springer; 2011. p. 307–16.
- Monk JM, Charusanti P, Aziz RK, Lerman JA, Premyodhin N, Orth JD, et al. Genome-scale metabolic reconstructions of multiple Escherichia coli strains highlight strain-specific adaptations to nutritional environments. Proc Natl Acad Sci. 2013;110(50):20338–43.PubMedPubMed CentralView ArticleGoogle Scholar
- Arumugam M, Raes J, Pelletier E, Le Paslier D, Yamada T, Mende DR, et al. Enterotypes of the human gut microbiome. Nature. 2011;473(7346):174–80.PubMedPubMed CentralView ArticleGoogle Scholar
- Vebø HC, Solheim M, Snipen L, Nes IF, Brede DA. Comparative genomic analysis of pathogenic and probiotic Enterococcus faecalis isolates, and their transcriptional responses to growth in human urine. PLoS One. 2010;5(8), e12489.PubMedPubMed CentralView ArticleGoogle Scholar
- Schellenberger J, Que R, Fleming RM, Thiele I, Orth JD, Feist AM, et al. Quantitative prediction of cellular metabolism with constraint-based models: the COBRA Toolbox v2. 0. Nat Protoc. 2011;6(9):1290–307.PubMedPubMed CentralView ArticleGoogle Scholar
- Orth JD, Thiele I, Palsson BØ. What is flux balance analysis? Nat Biotechnol. 2010;28(3):245–8.PubMedPubMed CentralView ArticleGoogle Scholar
- Thiele I, Palsson BØ. A protocol for generating a high-quality genome-scale metabolic reconstruction. Nat Protoc. 2010;5(1):93–121.PubMedPubMed CentralView ArticleGoogle Scholar
- Kumar VS, Dasika MS, Maranas CD. Optimization based automated curation of metabolic reconstructions. BMC bioinformatics. 2007;8(1):212.View ArticleGoogle Scholar
- Thiele I, Vlassis N, Fleming RM. fastGapFill: efficient gap filling in metabolic networks. Bioinformatics. 2014;30(17):2529–31.PubMedPubMed CentralView ArticleGoogle Scholar
- Reed JL, Patel TR, Chen KH, Joyce AR, Applebee MK, Herring CD, et al. Systems approach to refining genome annotation. Proc Natl Acad Sci. 2006;103(46):17480–4.PubMedPubMed CentralView ArticleGoogle Scholar
- Lewis NE, Nagarajan H, Palsson BO. Constraining the metabolic genotype–phenotype relationship using a phylogeny of in silico methods. Nat Rev Microbiol. 2012;10(4):291–305.PubMedPubMed CentralGoogle Scholar
- Suthers PF, Dasika MS, Kumar VS, Denisov G, Glass JI, Maranas CD. A genome-scale metabolic reconstruction of Mycoplasma genitalium, iPS189. PLoS Comput Biol. 2009;5(2), e1000285.PubMedPubMed CentralView ArticleGoogle Scholar
- Edwards J, Palsson B. The Escherichia coli MG1655 in silico metabolic genotype: its definition, characteristics, and capabilities. Proc Natl Acad Sci. 2000;97(10):5528–33.PubMedPubMed CentralView ArticleGoogle Scholar
- Heinken A, Sahoo S, Fleming RM, Thiele I. Systems-level characterization of a host-microbe metabolic symbiosis in the mammalian gut. Gut Microbes. 2013;4(1):28–40.PubMedPubMed CentralView ArticleGoogle Scholar
- Henry CS, DeJongh M, Best AA, Frybarger PM, Linsay B, Stevens RL. High-throughput generation, optimization and analysis of genome-scale metabolic models. Nat Biotechnol. 2010;28(9):977–82.PubMedView ArticleGoogle Scholar
- Ottar R, Giuseppe P, Manuela M, Bernhard OP, Ines T. Inferring the metabolism of human orphan metabolites from their metabolic network context affirms human gluconokinase activity. Biochem J. 2013;449(2):427–35.View ArticleGoogle Scholar
- Orth JD, Palsson B. Gap-filling analysis of the iJO1366 Escherichia coli metabolic network reconstruction for discovery of metabolic functions. BMC Syst Biol. 2012;6(1):30.PubMedPubMed CentralView ArticleGoogle Scholar
- Manichaikul A, Ghamsari L, Hom EF, Lin C, Murray RR, Chang RL, et al. Metabolic network analysis integrated with transcript verification for sequenced genomes. Nat Methods. 2009;6(8):589.PubMedPubMed CentralView ArticleGoogle Scholar
- Evaldson G, Heimdahl A, Kager L, Nord C. The normal human anaerobic microflora. Scand J Infect Dis Suppl. 1981;35:9–15.Google Scholar
- Heinken A, Thiele I. Systematic prediction of health-relevant human-microbial co-metabolism through a computational framework. Gut Microbes. 2015;6(2):120–30.PubMedPubMed CentralView ArticleGoogle Scholar
- Lee TJ, Paulsen I, Karp P. Annotation-based inference of transporter function. Bioinformatics. 2008;24(13):i259–i67.PubMedPubMed CentralView ArticleGoogle Scholar
- Stearns JC, Lynch MD, Senadheera DB, Tenenbaum HC, Goldberg MB, Cvitkovitch DG, et al. Bacterial biogeography of the human digestive tract. Sci Rep. 2011;1.
- Wong JM, de Souza R, Kendall CW, Emam A, Jenkins DJ. Colonic health: fermentation and short chain fatty acids. J Clin Gastroenterol. 2006;40(3):235–43.PubMedView ArticleGoogle Scholar
- Dupont CL, Rusch DB, Yooseph S, Lombardo M-J, Richter RA, Valas R, et al. Genomic insights to SAR86, an abundant and uncultivated marine bacterial lineage. ISME J. 2011;6(6):1186–99.PubMedPubMed CentralView ArticleGoogle Scholar
- Prentice MB. Bacterial comparative genomics. Genome Biol. 2004;5(8):338.PubMedPubMed CentralView ArticleGoogle Scholar
- Carr R, Borenstein E. Comparative analysis of functional metagenomic annotation and the mappability of short reads. PLoS One. 2014;9(8), e105776.PubMedPubMed CentralView ArticleGoogle Scholar
- El Yacoubi B, de Crécy-Lagard V. Integrative Data-Mining Tools to Link Gene and Function. Gene Function Analysis. Springer; 2014. p. 43–66. http://dx.doi.org/10.1007/978-1-62703-721-1_4
- Mazumdar V. Salomon Amar, and Daniel Segre. Metabolic proximity in the order of colonization of a microbial community. PLoS One. 2013;8(10), e77617.PubMedPubMed CentralView ArticleGoogle Scholar
- Gower JC, Legendre P. Metric and Euclidean properties of dissimilarity coefficients. J Classif. 1986;3(1):5–48.View ArticleGoogle Scholar
- Garrity GM, Bell JA, Lilburn TG. Taxonomic outline of the prokaryotes. Bergey’s manual of systematic bacteriology: Springer, New York, Berlin, Heidelberg; 2004.Google Scholar
- Marchandin H, Teyssier C, Campos J, Jean-Pierre H, Roger F, Gay B, et al. Negativicoccus succinicivorans gen. nov., sp. nov., isolated from human clinical samples, emended description of the family Veillonellaceae and description of Negativicutes classis nov., Selenomonadales ord. nov. and Acidaminococcaceae fam. nov. in the bacterial phylum Firmicutes. Int J Syst Evol Microbiol. 2010;60(6):1271–9.PubMedView ArticleGoogle Scholar
- Louis P, Flint HJ. Diversity, metabolism and microbial ecology of butyrate‐producing bacteria from the human large intestine. FEMS Microbiol Lett. 2009;294(1):1–8.PubMedView ArticleGoogle Scholar
- Yutin N, Galperin MY. A genomic update on clostridial phylogeny: Gram‐negative spore formers and other misplaced clostridia. Environ Microbiol. 2013;15(10):2631–41.PubMedPubMed CentralGoogle Scholar
- Ludwig W, Schleifer K, Whitman III W, Class III. Erysipelotrichia class nov. Bergey's Manual of Systematic Bacteriology. 2009;3:1298.Google Scholar
- Adler P, Bolten CJ, Dohnt K, Hansen CE, Wittmann C. Core fluxome and metafluxome of lactic acid bacteria under simulated cocoa pulp fermentation conditions. Appl Environ Microbiol. 2013;79(18):5670–81.PubMedPubMed CentralView ArticleGoogle Scholar
- Turroni F, Ribbera A, Foroni E, van Sinderen D, Ventura M. Human gut microbiota and bifidobacteria: from composition to functionality. Antonie Van Leeuwenhoek. 2008;94(1):35–50.PubMedView ArticleGoogle Scholar
- Zhu C, Delmont TO, Vogel TM, Bromberg Y. Functional basis of microorganism classification. PLoS Comput Biol. 2015;11(8), e1004472.PubMedPubMed CentralView ArticleGoogle Scholar
- Guckert JB, Ringelberg DB, White DC, Hanson RS, Bratina BJ. Membrane fatty acids as phenotypic markers in the polyphasic taxonomy of methylotrophs within the Proteobacteria. J Gen Microbiol. 1991;137(11):2631–41.PubMedView ArticleGoogle Scholar
- Tilg H, Kaser A. Gut microbiome, obesity, and metabolic dysfunction. J Clin Invest. 2011;121(6):2126–32.PubMedPubMed CentralView ArticleGoogle Scholar
- Kandler O. Carbohydrate metabolism in lactic acid bacteria. Antonie Van Leeuwenhoek. 1983;49(3):209–24.PubMedView ArticleGoogle Scholar
- Gupta RS. Origin of diderm (Gram-negative) bacteria: antibiotic selection pressure rather than endosymbiosis likely led to the evolution of bacterial cells with two membranes. Antonie Van Leeuwenhoek. 2011;100(2):171–82.PubMedPubMed CentralView ArticleGoogle Scholar
- Bush K. Antimicrobial agents targeting bacterial cell walls and cell membranes. Rev Sci Tech. 2012;31(1):43–56.PubMedGoogle Scholar
- D'Souza G, Waschina S, Pande S, Bohl K, Kaleta C, Kost C. Less is more: selective advantages Can explain the prevalent loss of biosynthetic genes in bacteria. Evolution. 2014.
- Segata N, Börnigen D, Morgan XC, Huttenhower C. PhyloPhlAn is a new method for improved phylogenetic and taxonomic placement of microbes. Nat Commun. 2013;4.
- Sokal RR, Rohlf FJ. The comparison of dendrograms by objective methods. Taxon. 1962;11(2):33–40.View ArticleGoogle Scholar
- Zaneveld JR, Lozupone C, Gordon JI, Knight R. Ribosomal RNA diversity predicts genome diversity in gut bacteria and their relatives. Nucleic Acids Res. 2010;38(12):3869–79.PubMedPubMed CentralView ArticleGoogle Scholar
- Natale DA, Shankavaram UT, Galperin MY, Wolf YI, Aravind L, Koonin EV. Towards understanding the first genome sequence of a crenarchaeon by genome annotation using clusters of orthologous groups of proteins (COGs). Genome Biol. 2000;1(5):RESEARCH0009.PubMedPubMed CentralView ArticleGoogle Scholar
- Finn RD, Bateman A, Clements J, Coggill P, Eberhardt RY, Eddy SR, et al. Pfam: the protein families database. Nucleic Acids Res. 2013;42:gkt1223.Google Scholar
- Caspi R, Altman T, Dale JM, Dreher K, Fulcher CA, Gilham F, et al. The MetaCyc database of metabolic pathways and enzymes and the BioCyc collection of pathway/genome databases. Nucleic Acids Res. 2010;38 suppl 1:D473–D9.PubMedView ArticleGoogle Scholar
- Wilmes P, Bowen BP, Thomas BC, Mueller RS, Denef VJ, VerBerkmoes NC, et al. Metabolome-proteome differentiation coupled to microbial divergence. MBio. 2010;1(5):e00246–10.PubMedPubMed CentralView ArticleGoogle Scholar
- Plata G, Henry CS. Vitkup D. Long-term phenotypic evolution of bacteria: Nature; 2014.Google Scholar
- Smillie CS, Smith MB, Friedman J, Cordero OX, David LA, Alm EJ. Ecology drives a global network of gene exchange connecting the human microbiome. Nature. 2011;480(7376):241–4.PubMedView ArticleGoogle Scholar
- Collins M, Lawson P, Willems A, Cordoba J, Fernandez-Garayzabal J, Garcia P, et al. The phylogeny of the genus clostridium: proposal of five new genera and eleven new species combinations. Int J Syst Bacteriol. 1994;44(4):812–26.PubMedView ArticleGoogle Scholar
- Van der Maaten L, Hinton G. Visualizing data using t-SNE. J Mach Learn Res. 2008;9(2579–2605):85.Google Scholar
- Amir E-aD, Davis KL, Tadmor MD, Simonds EF, Levine JH, Bendall SC, et al. viSNE enables visualization of high dimensional single-cell data and reveals phenotypic heterogeneity of leukemia. Nat Biotechnol. 2013;31(6):545–52.PubMed CentralView ArticleGoogle Scholar
- Platzer A. Visualization of SNPs with t-SNE. PLoS One. 2013;8(2), e56883.PubMedPubMed CentralView ArticleGoogle Scholar
- Laczny CC, Pinel N, Vlassis N, Wilmes P. Alignment-free visualization of metagenomic data by nonlinear dimension reduction. Sci Rep. 2014;4.
- Lee J-H, O'Sullivan DJ. Genomic insights into bifidobacteria. Microbiol Mol Biol Rev. 2010;74(3):378–416.PubMedPubMed CentralView ArticleGoogle Scholar
- An D, Na C, Bielawski J, Hannun YA, Kasper DL. Membrane sphingolipids as essential molecular signals for Bacteroides survival in the intestine. Proc Natl Acad Sci. 2011;108(Supplement 1):4666–71.PubMedView ArticleGoogle Scholar
- Nair B, Mayberry W, Dziak R, Chen P, Levine M, Hausmann E. Biological effects of a purified lipopolysaccharide from Bacteroides gingivalis. J Periodontal Res. 1983;18(1):40–9.PubMedView ArticleGoogle Scholar
- Chiu H-C, Levy R, Borenstein E. Emergent biosynthetic capacity in simple microbial communities. PLoS Comput Biol. 2014;10(7), e1003695.PubMedPubMed CentralView ArticleGoogle Scholar
- Blaut M, Clavel T. Metabolic diversity of the intestinal microbiota: implications for health and disease. J Nutr. 2007;137(3):751S–5S.PubMedGoogle Scholar
- Le Chatelier E, Nielsen T, Qin J, Prifti E, Hildebrand F, Falony G, et al. Richness of human gut microbiome correlates with metabolic markers. Nature. 2013;500(7464):541–6.PubMedView ArticleGoogle Scholar
- Markowitz VM, Chen I-MA, Palaniappan K, Chu K, Szeto E, Grechkin Y, et al. IMG: the integrated microbial genomes database and comparative analysis system. Nucleic Acids Res. 2012;40(D1):D115–D22.PubMedView ArticleGoogle Scholar
- Aziz RK, Bartels D, Best AA, DeJongh M, Disz T, Edwards RA, et al. The RAST Server: rapid annotations using subsystems technology. BMC Genomics. 2008;9(1):75.PubMedPubMed CentralView ArticleGoogle Scholar
- Gibbons R, Kapsimalis B. Estimates of the overall rate of growth of the intestinal microflora of hamsters, guinea pigs, and mice. J Bacteriol. 1967;93(1):510.PubMedPubMed CentralGoogle Scholar
- Overbeek R, Begley T, Butler RM, Choudhuri JV, Chuang H-Y, Cohoon M, et al. The subsystems approach to genome annotation and its use in the project to annotate 1000 genomes. Nucleic Acids Res. 2005;33(17):5691–702.PubMedPubMed CentralView ArticleGoogle Scholar
- Zhang H, Gao S, Lercher MJ, Hu S, Chen W-H. EvolView, an online tool for visualizing, annotating and managing phylogenetic trees. Nucleic Acids Res. 2012;40(W1):W569–W72.PubMedPubMed CentralView ArticleGoogle Scholar
- Albertsen M, Hugenholtz P, Skarshewski A, Nielsen KL, Tyson GW, Nielsen PH. Genome sequences of rare, uncultured bacteria obtained by differential coverage binning of multiple metagenomes. Nat Biotechnol. 2013;31(6):533–8.PubMedView ArticleGoogle Scholar