Sub-clinical detection of gut microbial biomarkers of obesity and type 2 diabetes
© Yassour et al. 2016
Received: 9 December 2015
Accepted: 25 January 2016
Published: 17 February 2016
Obesity and type 2 diabetes (T2D) are linked both with host genetics and with environmental factors, including dysbioses of the gut microbiota. However, it is unclear whether these microbial changes precede disease onset. Twin cohorts present a unique genetically-controlled opportunity to study the relationships between lifestyle factors and the microbiome. In particular, we hypothesized that family-independent changes in microbial composition and metabolic function during the sub-clinical state of T2D could be either causal or early biomarkers of progression.
We collected fecal samples and clinical metadata from 20 monozygotic Korean twins at up to two time points, resulting in 36 stool shotgun metagenomes. While the participants were neither obese nor diabetic, they spanned the entire range of healthy to near-clinical values and thus enabled the study of microbial associations during sub-clinical disease while accounting for genetic background.
We found changes both in composition and in function of the sub-clinical gut microbiome, including a decrease in Akkermansia muciniphila suggesting a role prior to the onset of disease, and functional changes reflecting a response to oxidative stress comparable to that previously observed in chronic T2D and inflammatory bowel diseases. Finally, our unique study design allowed us to examine the strain similarity between twins, and we found that twins demonstrate strain-level differences in composition despite species-level similarities.
These changes in the microbiome might be used for the early diagnosis of an inflamed gut and T2D prior to clinical onset of the disease and will help to advance toward microbial interventions.
The human gut microbiota plays an important role in health and disease [1, 2] and can be viewed as a mirror into the host physiology. One of the primary roles of the microbiota is energy harvest; thus, it is not surprising that microbial dysbiosis has been associated with various metabolic disorders, including type 2 diabetes (T2D) [3, 4] and obesity [5–7]. T2D is often a consequence of obesity. As the diagnosis is threshold-based, risk of developing T2D in the near future correlates with high levels of two biomarkers, fasting blood sugar (FBS) and HbA1c, even when they do not meet the clinical criteria (HbA1c >6.5 % or FBS >125). However, the microbial changes that occur in the sub-clinical state, prior to the onset of disease, have never been examined, but may potentially be used for early diagnosis and intervention.
Previous profiles of the gut microbiome during clinical T2D have found compositional changes between patients and healthy controls [3, 4], including an obesity-related change in the abundance ratio of Bacteroidetes:Firmicutes [5, 8], and a decreased abundance of mucin-degrading Akkermansia muciniphila in overweight children  and pregnant women [10, 11]. However, there is no strong consensus across studies in taxa changing in obese versus lean individuals . The causes for this inconsistency may be either technical or biological. From a technical standpoint, a lack of consistent standard operating procedures for sample preparation and sequencing can lead to great variance between different labs and studies . Biologically, the specific composition of the community may be much less important than its overall functional capability.
Indeed, there is greater consensus between these studies when microbial functional dysbioses are considered rather than microbial composition . The gut communities of T2D patients showed increased capacity for oxidative stress resistance, and a decreased capacity for flagellar assembly and riboflavin metabolism [3, 4]. Interestingly, oxidative stress resistance was also enriched in the guts of patients with inflammatory bowel diseases (IBD) , potentially indicating that the microbiome is generally stressed by low-level inflammation and immune activation, which may be present at the sub-clinical state of T2D as well.
Despite recent studies associating the microbiome with T2D [3, 4] and obesity [6, 8], all previous work has examined individuals with well-established disease. These data may be further influenced by additional factors, such as decreased subject mobility, and it is difficult to conclude from study design whether the observed microbial changes preceded the onset of disease. Furthermore, these studies have rarely taken into account the various genetic backgrounds of the patients. We have addressed these issues by performing the first metagenomic profile of the gut microbiome of monozygotic (MZ) twins, spanning the entire healthy range of T2D clinical indicators, including body mass index (BMI) and fasting blood sugar (FBS). Identifying gradient-like associations between these parameters and gut microbiome features in the sub-clinical state of these diseases will open the way to discover potential markers for early diagnosis of T2D and obesity.
We found several taxa associated with sub-clinical changes in BMI, blood pressure, sugar, and triglycerides, including enrichment of the Roseburia genus and depletion of the Akkermansia muciniphila species. Additionally, riboflavin and NAD biosynthesis were metagenomically enriched in participants with high blood pressure and BMI values. Interestingly, similar functional enrichments are shared with other gut inflammatory conditions such as IBD  and clinically-established T2D , suggesting shifts in the gut microbial population prior to full disease onset that may be either causal or an early correlative indicator.
Finally, this cohort included a unique combination of MZ twins and longitudinal sampling, which allowed us to identify the degree to which specific microbial strains were shared between the guts of siblings and maintained over time. Despite the small size of this targeted cohort, the deep metagenomic sequencing (mean 3.5 Gnt per sample) combined with a focus on taxa of high relative abundance (see Methods) enabled us to determine strain similarity with high resolution in these data. MZ twins have also been previously observed to share a greater proportion of gut microbes than unrelated individuals , and some strains appear to be maintained within the guts of isolated individuals for months to years . Surprisingly, we observed that while twins in our MZ cohort indeed share a substantial subset of microbial species, strains within these species differ between related twins. Thus, the gut microbial similarities of twins may arise from sources such as genetic pressure to acquire certain species (but not specific strains), or from early colonization by the same strains, with subsequent genetic divergence over the course of a lifetime.
Written informed consent was obtained from each participant. The study protocol was approved by the institutional review board (IRB) of Samsung Medical Center, Busan Paik Hospital, and Seoul National University (IRB No. 144-2011-07-11).
Study participants and specimen collection
The participants were MZ twins who enrolled for the Healthy Twin study in Seoul and Busan, South Korea. The zygosity of twins was determined using AmpFlSTR Identifier Kit with 16 short tandem repeat markers (15 autosomal STR markers + one sex determining marker) or a questionnaire with a validated accuracy of >90 % . Details on methodology of this cohort have been previously described .
A total of 36 fecal samples from the participants were collected: samples from two twin pairs were taken once and those from eight twin pairs were obtained twice with an average interval of 2 years. Twins were in the age range of 30-48 years at the first sampling point. Fecal samples were taken in conjunction with a health examination and immediately stored at -25 °C. They were subsequently transported to the two central clinics and stored at -80 °C until DNA extraction. Blood samples were drawn by vein puncture after an overnight fast and sent to a central laboratory to measure biochemical factors.
During each visit, individuals also completed a questionnaire recording life style, medication, and dietary habits. Anthropometrical measurements (height, weight, waist circumference, and so on) and biochemical tests (glucose, hsCRP, total cholesterol, HDL-C, LDL-C, triglyceride, and so on) were also conducted (Additional file 1: Table S1). The derived homeostasis model assessment (HOMA) index uses the fasting blood sugar and insulin to predict the insulin resistance of patients  and was calculated as standard (insulin * glucose)/405, both measured after fasting and glucose levels measured in mg/dL .
Nucleic acid extraction and metagenomic shotgun sequencing
Total DNA was extracted from each fecal sample using the MoBio Power Soil DNA Isolation kit (MoBio, Solana Beach, CA, USA) according to the manufacturer’s instructions and stored at -80 °C until subsequent analysis. All samples were sequenced using the Illumina Hiseq2000 instrument, which produced paired-end reads of 101 nt, yielding average 3.5 Gnt per each of fecal samples.
Metagenomic shotgun sequences analysis
The gut microbial composition of each sample was profiled using MetaPhlAn . MetaPhlAn uses a unique set of markers for each species (and higher level clades) to estimate the abundance of species in each sample according to the number of mapped reads to its markers. The relative abundances of the gut microbial functional pathways from metagenomically sequenced communities were determined using HUMAnN . HUMAnN maps the sequenced reads to a non-redundant set of genes extracted from the KEGG database  and estimates the abundance level of each functional module by the number of matches to member genes fully compatible with it being carried out by one or more microbes.
Testing for significant associations with the clinical metadata variables
To identify significant associations between microbial and phenotypic variables, we applied a linear multivariate regression model specifically adapted to microbiome data (MaAsLin, Multivariate microbial Association by Linear models ). MaAsLin constructs boosted, additive general linear models to associate metadata and transformed microbial taxonomic or functional relative abundances. Since microbial community profiles are typically high-dimensional, boosting is used for feature selection over potential covariates to identify those most associated with each microbial feature. Selected metadata are then used in a general linear model with metadata as predictors and arcsin-square root transformed microbial relative abundances as the responses. In this study, model covariates of interest comprised of clinical variables included in Additional file 1: Table S1, and each model also included age, smoking status, sex, and twin as potential confounders (the latter as a random effect to accommodate repeated longitudinal measures).
Comparing strains between samples
We performed taxonomic profiling with MetaPhlAn . Briefly, MetaPhlAn operates by mapping raw sequence reads to a database of pre-defined clade-specific marker genes. Markers are those genes occurring in isolates from a particular clade but not outside of that clade. After mapping reads to clade-specific marker genes, the resulting raw counts are normalized for total marker gene length and outliers, yielding profiles of: (1) clade relative abundance; (2) marker gene presence/absence; and (3) marker gene abundance (in reads-per-kb (RPK) units, where 10 RPK would correspond to about ×1 coverage, given our 100 bp reads). Due to gene gain and loss events, an individual strain will not necessarily carry all of the markers associated with its corresponding species. A specific pattern of marker presence and absence can therefore be used as a molecular ‘barcode’ to identify a strain across samples. We next compared the marker gene abundance profiles of various samples (unrelated, twin or self; with median marker abundance >5 RPK) using a Bray-Curtis distance.
Generating and analyzing the taxon-function correlation matrix
Spearman correlation was calculated between the profiles of each microbe and each function to generate the taxon-function correlation matrix. The KEGG database  was used to identify the ‘encoded’ correlations, by calculating the fraction of its reference sequences that include sufficient genes for any given module. The microbial co-occurrence matrix was calculated using spearman correlation between all taxa profiles to identify the ‘associated’ correlations.
Comparing Korean and Western microbial populations
The prevalence and average abundances of all clades were calculated within our cohort and the HMP , and these were compared using Pearson correlation. Prevalence was defined as percent of samples with >0.001 relative abundance for each species, and average abundance was calculated only for samples passing that criterion.
Sequence accession numbers and availability
Sequences generated in this study are publicly available at the European nucleotide archive (ERP002391).
Monozygotic twin cohort and longitudinal metagenomic profiles
As is typical in the gut, the taxa that vary most widely in our cohort were members of the Bacteroidetes, Firmicutes, and at lower abundances the Esherichia, Methanobrevibacter, and Bifidobacterium genera (Fig. 1b). The x-axis depicts mostly the Firmicutes (left) to Bacteroidetes (ratio), as observed before , with a modest positive association between Bacteroides and BMI (see below). The y-axis is dominated by less prevalent or more variable genera, like Methanobrevibacter (25 % prevalence) and Bifidobacterium (2.7-3.3 coefficient of variation; Additional file 3: Table S2). The Esherichia and Prevotella genera are prevalent in our cohort (92 and 81 %, respectively, see Additional file 3: Table S2), but only a few individuals have high abundances of one or both genera, driving their contribution to population variability. Overall, we found that the prevalence and abundance profiles of the various taxa in our data are consistent with those measured in Western populations ; P value = 0.0001 by Pearson correlation; Additional file 4: Figure S2).
In our cohort, clinical parameters such as BMI (Fig. 1c) and microbial community composition (Fig. 1d) were both more similar between twins than between unrelated individuals, and both were self-similar over time. BMI is most stable between ‘self’ samples, especially since no individual has become obese during this study, and indeed, twins are more concordant on BMI compared to unrelated . When comparing their microbial composition we found that, as expected, twins were somewhat more similar than unrelated , but self-samples were significantly more stable , indicating that for both clinical and microbial phenotypes, longitudinal samples were more similar than twins, which were in turn more similar than unrelated.
Phylogenetic and functional diversity in the Korean gut microbiome
Less prevalent organisms included Methanobrevibacter smithii, the predominant archaeon in the human gut  with the important role of reducing the hydrogen via methanogenesis. M. smithii here mirrored its relatively low prevalence as seen in Western populations (25 %, compared to 21 % in the Human Microbiome Project (HMP)) , but with higher average relative abundance when present (9 %, compared to 0.3 %) . Despite the major differences in the Korean diet and environment, M. smithii is still the major microbe likely responsible for this function.
Several organisms, including Bifidobacterium longum, Escherichia coli, Prevotella copri, and Bacteroides plebeius, were significantly more prevalent in our population compared to that of the HMP (Additional file 3: Table S2). It has been previously shown that Prevotella is more common in international cohorts , and indeed we found P. copri present in 81 % of our samples, compared to 16 % in the HMP . Moreover, we found lower abundance of P. copri in our cohort when present (approximately 10 % relative abundance), unlike the high 37 % relative abundance found in the individuals who carry it in the HMP’s Western population , although the overall average of P. copri relative abundance in all individuals is similar between cohorts (6-8 %).
Another difference between the Korean and Western populations was Bifidobacterium longum, a microbe that ferments sugars into lactic acid noted as one of the first colonizers of the infant gut . It was carried at an unusually high prevalence of 94 % (with average relative abundance of 2.5 %), in contrast to its presence in the HMP of 59 % prevalence and 0.4 % relative abundance . Bacteroides plebeius was likewise enriched here, with 97 % prevalence in our cohort and only 9 % in the HMP . P. plebius has been previously found in Japanese populations, likely due to its capability to break down complex carbohydrates specific to seaweed , and it may play a similar role in the guts of Koreans, as this is a major staple of their diet.
As has been previously observed , despite variability in the composition of the microbiome among these individuals, the distribution of microbial metabolic processes remained relatively stable (Fig. 2b, Additional file 5: Figure S3). One of the most variable modules was the transport system of autoinducer-2 (AI-2), a quorum-sensing signaling molecule traditionally associated with the Enterobacteriaceae and Vibrionaceae  and recently characterized in some Bifidobacterium species . We also see striking variability in the biosynthesis of GABA, a major neurotransmitter in the central nervous system that has also been recently shown to be produced by some Bifidobacterium species . Both these modules had a striking 3.8 coefficient of variance and approximately 18 % prevalence (compared to 94 % and 60 % in Western population, respectively) and shared similar abundance profiles in our data, suggesting a potential link between the two processes, perhaps through carriage by specific Bifidobacterium strains. As a control, we also see several ubiquitous ‘housekeeping’ processes such as the ribosome and translation, glycolysis, and gluconeogenesis were present at high levels with low variability among individuals (>4 % average relative abundance and 100 % prevalence, Additional file 3: Table S2 and Additional file 5: Figure S3).
Host factors such as BMI associate with some microbes and processes in a graded fashion
Other continuous associations of clades with phenotypes included a positive correlation between the Bacteroides genus and BMI, again spanning a range of the latter outside of clinical obesity. The Bacteroidetes:Firmicutes ratio is one of the earliest features of the gut microbiome suggested to associate with obesity in mice , but in this and other studies [11, 32, 33] associating an increase in Bacteroidetes with obesity rather than an increase in Firmicutes. Both positive and negative associations with the Bacteroidetes have been found in human populations , suggesting that this finding is not generalizable and depends greatly on factors that may include the underlying demographics, diet, sample preparation, and analysis (see Discussion).
We found multiple associations between microbial molecular function and clinical phenotypes (Fig. 3, Additional file 6: Table S3), including an increase in riboflavin-, NAD-, and tetrahydrofolate-biosynthesis and a decrease in pyruvate ferredoxin oxidoreductase accompanying increasing BMI levels (Fig. 3). Riboflavin and NAD are both required for the biosynthesis of the reduced form of glutathione (GSH) , an important antioxidant that alleviates the damage done by reactive oxygen species, and indeed we found that glutathione biosynthesis also increased at higher BMI levels (Additional file 6: Table S3). The direction of these associations suggests that the gut microbiota is producing more glutathione, and potentially processing more from the host, to relieve the increased oxidative stress at high levels of BMI.
Besides the gradient associations discussed above, we also identified two unusual threshold-like associations: BMI with Eubacterium siraeum and blood pressure with Methanobrevibacter smithii. In both these cases, the microbe is present only below a certain threshold (22.5 BMI, and 70 diastolic blood pressure, respectively). Interestingly, these trends do not appear in data collected for the human microbiome project (Additional file 8: Figure S5) and maybe unique to our cohort. Such discrete associations can be very interesting to investigate further, to potentially reveal the mechanism underlying the microbe’s sensing of host conditions.
Contribution of specific microbes to overall gut community function
As expected, many ‘encoded’ correlations are found in our data and each of them induces others that are ‘associated’, based on the species co-occurrence network as computed using Spearman correlation (Additional file 10: Figure S7B). One interesting ‘encoded’ example is a set of archaeal functions, such as coenzyme M biosynthesis and archaeal RNA polymerase and ribosome, correlated with the archaeal species Methanobrevibacter smithii and Methanosphaera stadtmanae (Fig. 4 Box A). Correspondingly, we found ‘associated’ correlations between these archaeal functions and the abundance of Deltaproteobacteria, Akkermansia muciniphila and Eubacterium siraeum, which co-occur with the archaeal species in our data (Fig. 4 Box A; Additional file 10: Figure S7B Box A). We do not yet know why particularly these microbes would tend to share the archaeal environment, but it is interesting to examine if this relationship holds in other cohorts as well.
Additionally, three metal transport systems (TSs) were correlated with specific taxa. The manganese/zinc/iron TS M00319 is an ABC transporter, comprising four proteins, found originally in Treponema pallidum , and this specific TS is encoded only in the Veillonella species in our data. In addition to ‘encoded’ correlations between this TS and three Veillonella species, we identified two ‘associated’ correlations with Streptococcus salivarius and Haemophilus parainfluenzae, co-occurring species with the Veillonella (Fig. 4 Box B; Additional file 10: Figure S7B Box B). Similarly, the cobalt/nickel TS M00245 is another four-protein ABC transporter, estimated to be the most widespread uptake system for the two metals . This module is encoded in many microbial species, and specifically in Eubacterium hallii and Ruminococcus obeum  in our communities, resulting in those encoded correlations. These were accompanied by additional ‘associated’ correlations with co-occurring species including Dorea longicatena (Fig. 4 Box C; Additional file 10: Figure S7B Box C). Lastly, two modules involving nickel TSs (cobalt/nickel M00245 and peptide/nickel M00239), together with sugar and amino acid metabolism and TS modules, were ‘encoded’ correlations found in the Bifidobacterium species. Methionine degradation M00035 was also encoded by these organisms, which generates S-adenosyl-L-methionine (SAM), a major methyl donor in the cell , and Bifidobacterium is used outside of the gut as a source of SAM in the functional food industry .
This cohort’s microbial co-occurrence network (Additional file 10: Figure S7B) can also explain some of the negative associations found in our data. For example, the abundance of Bacteroides xylanisolvens is negatively correlated with several taxa, including Methanosphaera stadtmanae and Ruminococcus obeum (Additional file 10: Figure S7B Box D), resulting in a negative correlation between these species and modules encoded by B. xylanisolvens, like NAD-, tetrahydrofolate-, and biotin biosynthesis (Fig. 4 Boxes D, E, F). Another example is the negative correlation between Roseburia intestinalis and E. coli (Additional file 10: Figure S7B Box E), resulting in the negative associations between E. coli and the Shikimate pathway, encoded in R. intestinalis genome sequence (Fig. 4 Boxes G, H). Such negative associations can also be the result of conflicting functionalities between certain microbes and metabolic functions not present in their genomes.
Interestingly, some correlations are neither ‘encoded’ nor ‘associated’, and we can only hypothesize as to their cause. One such example is the correlation between the archaeal functions mentioned above and abundance of Bacteroides fragilis, Odoribacter splanchnicus and Parabacteroides merdae (Fig. 4 Box I). Although these functions are not encoded in any of these genomes, the gene comB (2-phosphosulfolactate phosphatase), which belongs to the coenzyme M biosynthesis module, is encoded in the genome of Parabacteroides merdae , potentially explaining this association. Such a correlation might arise due to these organisms’ metabolic dependence on a function encoded by diverse organisms in different hosts.
Another interesting example is the positive correlation between riboflavin biosynthesis and the abundance of Bacteroides xylanisolvens (Fig. 4 Box E). As discussed above, we found several compositional and functional associations with BMI, including Bacteroides, riboflavin-, NAD-, and tetrahydrofolate- biosynthesis. Indeed, the latter two are ‘encoded’ correlations with B. xylanisolvens, suggesting that the increase in this species abundance in higher levels of BMI is contributing to the increase in NAD- and tetrahydrofolate- biosynthesis. Additionally, riboflavin biosynthesis may also be an ‘encoded’ association of B. xylanisolvens and we failed to identify it as such, potentially due to the incompleteness of current functional databases.
Finally, additional ‘encoded’ correlations were found between various amino acid transport systems and E. coli (Fig. 4 Box K). Many of these were detected in at least some strains of E. coli, due to the combination of E. coli’s very large pan-genome and the extent to which its strain variation space is well-covered by the many available reference genomes. Many of these correlations, like all of this cohort’s microbe-function correlations, were also found in the HMP  (yellow marks in Fig. 4, and Additional file 11: Figure S8). This suggests both that simple ‘encoded’ correlations recur across populations, as expected, and that more subtle ‘associated’ microbial dependencies may be consistent among diverse gut ecologies.
Microbial species, but typically not strains, are shared between twins
Several studies have observed that related individuals, and particularly twins, carry more similar microbial communities than do unrelated individuals , and we reproduce this finding in our cohort (Fig. 1d). However, it has not been previously determined whether this similarity is due to ecological pressures that select for similar microbes among individuals, dispersal effects that cause the acquisition of identical microbes, or other factors. In cases where twins in this study shared similar taxonomic profiles and identical species, we thus tested whether these species were of the same strain. Defining microbial clades at genus-, species-, or strain- level is a difficult task [11, 41], and here we chose to define a strain as a combination of genomic markers, allowing us to identify dominant, near-clonal populations. Microbes were strain-typed within samples by identifying conserved or differential patterns of unique mobile element loss and gain using MetaPhlAn (see Methods), which has previously been successful in differentiating strains among individuals and over time [21, 42]. In particular, individuals were shown in these previous analyses to often carry a single dominant strain of most species  and for that strain to be significantly stable over time [14, 21] (Additional file 12: Figure S9). Applying this method to our data allowed us to determine: (1) when twins shared the same strain, in addition to the same species; and (2) when these strains were retained within an individual over time.
We quantified the degree to which each microbial species was represented by the same or different strains among twins, specifically by calculating for each microbe the mean distance across genomic elements of all sample pair-wise comparisons (Additional file 14: Figure S11, see Methods). While some species retained very similar strains among twins (for example, Methanobrevibacter smithii, Fig. 5b), others consistently had more distinct strains (for example, Prevotella copri, Fig. 5c). While this detailed analysis is only possible in microbes that are particularly abundant in multiple samples, it raises the intriguing possibility that identical strains acquired early in life persist in the gut, but evolve rapidly through the gain and loss of genetic elements.
Here, we analyzed metagenomes from 36 fecal samples drawn from healthy Korean MZ twins over time, identifying associations between T2D-related biomarkers (for example, BMI, FBS) and microbial clades and functions. We found, among other examples, that BMI was negatively correlated with the abundance of Akkermansia muciniphila and positively correlated with riboflavin and NAD biosynthesis. These associations occurred over both the pre- and post-onset range of T2D clinical markers, suggesting that the microbiota may contribute to or react to changes in the host environment prior to the onset of disease. Furthermore, functional changes in the gut microbiome at higher sub-clinical values of BMI, FBS, and triglycerides resembled the signatures found in patients with established IBD or T2D, suggesting a shared response to oxidative stress in the gut, induced even at low levels of inflammation or immune activation. Finally, we found that while twins were more similar than unrelated individuals in microbial composition, they often carried different strains of these species. The computational framework presented here can be easily applied to other MZ twin cohorts, identifying early microbial markers of various other diseases, even in their sub-clinical phase.
Obesity and metabolic syndrome have long been associated with chronic, low-grade inflammation . For example, macrophages of obese individuals accumulate in adipose tissue, where they express pro-inflammatory cytokines such as TNFa, IL6, and INOS , and the gut microbiota can initiate the inflammation and insulin resistance associated with obesity . Interestingly, although there was no indication of host inflammation in our data (as measured by hsCRP), we observed the microbiome responding to this in the form of decreased abundance of Akkermansia muciniphila and increased NAD- and riboflavin-biosynthesis. This collection of functional changes together specifically enables the recharging of glutathione to its reduced form, promoting redox homeostasis in microbes potentially exposed to an increasingly hostile, inflammatory, oxidatively stressed environment the gut.
It is likely that additional compositional and functional shifts accompany this low-grade inflammation in T2D and related conditions, which will be better detected in other, larger cohorts capturing an even broader range of phenotypes and disease states. Finding common microbial changes is an important step towards understanding the cross talk between the gut microbiota and the diseased host, but whether these shifts are causal, responsive, or both remains an open question. For example, the microbial response to redox stress is more likely to be reactive, but microbes that are robust to this environment may promote its maintenance and thus contribute to immune activation or obesity. A combination of interventions in model systems and longitudinal prospective cohort studies of high-risk individuals, identifying the microbial changes that occur before the onset and during the early progression of the diseases, will enable us to determine whether the microbial shifts trigger host symptoms, or vice versa, and potentially by what specific molecular mechanisms.
Many microbial genes, particularly housekeeping genes, are transcribed at a basal level, and thus their measured DNA and RNA levels are well-correlated. In contrast, other classes of genes such as vitamin and amino acid biosynthesis are much more dynamically regulated, so DNA and RNA levels are less correlated . While this study measured only DNA abundance of genes, concurrent examination of the subclinical biomarker meta-transcriptome would be an informative extension of this work.
Finally, a unique dataset comprising shotgun metagenomes of MZ twins over time enabled us to find that although twins are more similar in their species composition, they often harbor different strains. This unexpected discrepancy between the microbial population structures of strains versus species in the human microbiome should be further explored, as it can be explained by a variety of very distinct ecological and molecular hypotheses. Hosts with similar genetic profiles may exert a modest but continuous selective pressure for the acquisition and maintenance of similar species in the gut, which could be tested by collecting data on the degree of shared early life and persistent environment versus genetics. Alternatively, initially identical strains acquired from a shared environment, possibly in early life, may be maintained but diverge through fixation of genetic drift and laterally transferred elements over time. Occasional strain differences within the same individual over time suggest a fast divergence rate; however, a larger study, with temporally dense sampling of both adults and infants, will be needed to address this question.
To conclude, this study provides evidence of low-grade inflammation of the gut with increasing values of obesity- and T2D-related biomarkers. Compositional and functional microbial signatures indicate the presence of sub-clinical inflammation in adults increasingly at risk of these conditions, even before they are reflected by clinical markers. If these microbial shifts play a causal role in the onset of obesity or T2D, they may represent not only novel markers for early diagnosis, but also a target for preventative therapeutic intervention. Even if these shifts are not ultimately the primary causal agents behind their associated diseases, microbial dysbioses may still be manipulated to avert disease onset, and their specifics are likely to improve our mechanistic understanding of host-microbiota interaction and its role in disease prevention and treatment.
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.
- Clemente JC, Ursell LK, Parfrey LW, Knight R. The impact of the gut microbiota on human health: an integrative view. Cell. 2012;148(6):1258–70.View ArticlePubMedGoogle Scholar
- Lim MY, Rho M, Song YM, Lee K, Sung J, Ko G. Stability of gut enterotypes in Korean monozygotic twins and their association with biomarkers and diet. Sci Rep. 2014;4:7348.PubMed CentralView ArticlePubMedGoogle Scholar
- 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.View ArticlePubMedGoogle Scholar
- Karlsson FH, Tremaroli V, Nookaew I, Bergstrom G, Behre CJ, Fagerberg B, et al. Gut metagenome in European women with normal, impaired and diabetic glucose control. Nature. 2013;498(7452):99–103.View ArticlePubMedGoogle Scholar
- Turnbaugh PJ, Ley RE, Mahowald MA, Magrini V, Mardis ER, Gordon JI. An obesity-associated gut microbiome with increased capacity for energy harvest. Nature. 2006;444(7122):1027–31.View ArticlePubMedGoogle Scholar
- Turnbaugh PJ, Hamady M, Yatsunenko T, Cantarel BL, Duncan A, Ley RE, et al. A core gut microbiome in obese and lean twins. Nature. 2009;457(7228):480–4.PubMed CentralView ArticlePubMedGoogle Scholar
- Verdam FJ, Fuentes S, de Jonge C, Zoetendal EG, Erbil R, Greve JW, et al. Human intestinal microbiota composition is associated with local and systemic inflammation in obesity. Obesity (Silver Spring). 2013;21(12):E607–15.View ArticleGoogle Scholar
- Ley RE, Turnbaugh PJ, Klein S, Gordon JI. Microbial ecology: human gut microbes associated with obesity. Nature. 2006;444(7122):1022–3.View ArticlePubMedGoogle Scholar
- Karlsson CL, Onnerfalt J, Xu J, Molin G, Ahrne S, Thorngren-Jerneck K. The microbiota of the gut in preschool children with normal and excessive body weight. Obesity (Silver Spring). 2012;20(11):2257–61.View ArticleGoogle Scholar
- Santacruz A, Collado MC, Garcia-Valdes L, Segura MT, Martin-Lagos JA, Anjos T, et al. Gut microbiota composition is associated with body weight, weight gain and biochemical parameters in pregnant women. Br J Nutr. 2010;104(1):83–92.View ArticlePubMedGoogle Scholar
- Collado MC, Isolauri E, Laitinen K, Salminen S. Distinct composition of gut microbiota during pregnancy in overweight and normal-weight women. Am J Clin Nutr. 2008;88(4):894–9.PubMedGoogle Scholar
- Finucane MM, Sharpton TJ, Laurent TJ, Pollard KS. A taxonomic signature of obesity in the microbiome? Getting to the guts of the matter. PLoS One. 2014;9(1):e84689.PubMed CentralView ArticlePubMedGoogle Scholar
- Morgan XC, Tickle TL, Sokol H, Gevers D, Devaney KL, Ward DV, et al. Dysfunction of the intestinal microbiome in inflammatory bowel disease and treatment. Genome Biol. 2012;13(9):R79.PubMed CentralView ArticlePubMedGoogle Scholar
- Faith JJ, Guruge JL, Charbonneau M, Subramanian S, Seedorf H, Goodman AL, et al. The long-term stability of the human gut microbiota. Science. 2013;341(6141):1237439.PubMed CentralView ArticlePubMedGoogle Scholar
- Song YM, Lee D, Lee MK, Lee K, Lee HJ, Hong EJ, et al. Validity of the zygosity questionnaire and characteristics of zygosity-misdiagnosed twin pairs in the Healthy Twin Study of Korea. Twin Res Hum Genet. 2010;13(3):223–30.View ArticlePubMedGoogle Scholar
- Sung J, Cho SI, Lee K, Ha M, Choi EY, Choi JS, et al. Healthy Twin: a twin-family study of Korea--protocols and current status. Twin Res Hum Genet. 2006;9(6):844–8.View ArticlePubMedGoogle Scholar
- Matthews DR, Hosker JP, Rudenski AS, Naylor BA, Treacher DF, Turner RC. Homeostasis model assessment: insulin resistance and beta-cell function from fasting plasma glucose and insulin concentrations in man. Diabetologia. 1985;28(7):412–9.View ArticlePubMedGoogle Scholar
- Segata N, Waldron L, Ballarini A, Narasimhan V, Jousson O, Huttenhower C. Metagenomic microbial community profiling using unique clade-specific marker genes. Nat Methods. 2012;9(8):811–4.PubMed CentralView ArticlePubMedGoogle Scholar
- Abubucker S, Segata N, Goll J, Schubert AM, Izard J, Cantarel BL, et al. Metabolic reconstruction for metagenomic data and its application to the human microbiome. PLoS Comput Biol. 2012;8(6):e1002358.PubMed CentralView ArticlePubMedGoogle Scholar
- Kanehisa M, Goto S, Furumichi M, Tanabe M, Hirakawa M. KEGG for representation and analysis of molecular networks involving diseases and drugs. Nucleic Acids Res. 2010;38(Database issue):D355–60.PubMed CentralView ArticlePubMedGoogle Scholar
- Human Microbiome Project Consortium. Structure, function and diversity of the healthy human microbiome. Nature. 2012;486(7402):207–14.View ArticleGoogle Scholar
- Koren O, Knights D, Gonzalez A, Waldron L, Segata N, Knight R, et al. A guide to enterotypes across the human body: meta-analysis of microbial community structures in human microbiome datasets. PLoS Comput Biol. 2013;9(1), e1002863.PubMed CentralView ArticlePubMedGoogle Scholar
- Samuel BS, Hansen EE, Manchester JK, Coutinho PM, Henrissat B, Fulton R, et al. Genomic and metabolic adaptations of Methanobrevibacter smithii to the human gut. Proc Natl Acad Sci U S A. 2007;104(25):10643–8.PubMed CentralView ArticlePubMedGoogle Scholar
- De Filippo C, Cavalieri D, Di Paola M, Ramazzotti M, Poullet JB, Massart S, et al. Impact of diet in shaping gut microbiota revealed by a comparative study in children from Europe and rural Africa. Proc Natl Acad Sci U S A. 2010;107(33):14691–6.PubMed CentralView ArticlePubMedGoogle Scholar
- Schell MA, Karmirantzou M, Snel B, Vilanova D, Berger B, Pessi G, et al. The genome sequence of Bifidobacterium longum reflects its adaptation to the human gastrointestinal tract. Proc Natl Acad Sci U S A. 2002;99(22):14422–7.PubMed CentralView ArticlePubMedGoogle Scholar
- Hehemann JH, Correc G, Barbeyron T, Helbert W, Czjzek M, Michel G. Transfer of carbohydrate-active enzymes from marine bacteria to Japanese gut microbiota. Nature. 2010;464(7290):908–12.View ArticlePubMedGoogle Scholar
- Surette MG, Miller MB, Bassler BL. Quorum sensing in Escherichia coli, Salmonella typhimurium, and Vibrio harveyi: a new family of genes responsible for autoinducer production. Proc Natl Acad Sci U S A. 1999;96(4):1639–44.PubMed CentralView ArticlePubMedGoogle Scholar
- Sun Z, He X, Brancaccio VF, Yuan J, Riedel CU. Bifidobacteria exhibit LuxS-dependent autoinducer 2 activity and biofilm formation. PLoS One. 2014;9(2):e88260.PubMed CentralView ArticlePubMedGoogle Scholar
- Barrett E, Ross RP, O’Toole PW, Fitzgerald GF, Stanton C. gamma-Aminobutyric acid production by culturable bacteria from the human intestine. J Appl Microbiol. 2012;113(2):411–7.View ArticlePubMedGoogle Scholar
- Everard A, Belzer C, Geurts L, Ouwerkerk JP, Druart C, Bindels LB, et al. Cross-talk between Akkermansia muciniphila and intestinal epithelium controls diet-induced obesity. Proc Natl Acad Sci U S A. 2013;110(22):9066–71.PubMed CentralView ArticlePubMedGoogle Scholar
- Ley RE, Backhed F, Turnbaugh P, Lozupone CA, Knight RD, Gordon JI. Obesity alters gut microbial ecology. Proc Natl Acad Sci U S A. 2005;102(31):11070–5.PubMed CentralView ArticlePubMedGoogle Scholar
- Schwiertz A, Taras D, Schafer K, Beijer S, Bos NA, Donus C, et al. Microbiota and SCFA in lean and overweight healthy subjects. Obesity (Silver Spring). 2010;18(1):190–5.View ArticleGoogle Scholar
- Zhang H, DiBaise JK, Zuccolo A, Kudrna D, Braidotti M, Yu Y, et al. Human gut microbiota in obesity and after gastric bypass. Proc Natl Acad Sci U S A. 2009;106(7):2365–70.PubMed CentralView ArticlePubMedGoogle Scholar
- Ley RE. Obesity and the human microbiome. Curr Opin Gastroenterol. 2010;26(1):5–11.View ArticlePubMedGoogle Scholar
- Mulherin DM, Thurnham DI, Situnayake RD. Glutathione reductase activity, riboflavin status, and disease activity in rheumatoid arthritis. Ann Rheum Dis. 1996;55(11):837–40.PubMed CentralView ArticlePubMedGoogle Scholar
- Hazlett KR, Rusnak F, Kehres DG, Bearden SW, La Vake CJ, La Vake ME, et al. The Treponema pallidum tro operon encodes a multiple metal transporter, a zinc-dependent transcriptional repressor, and a semi-autonomously expressed phosphoglycerate mutase. J Biol Chem. 2003;278(23):20687–94.View ArticlePubMedGoogle Scholar
- Rodionov DA, Hebbeln P, Gelfand MS, Eitinger T. Comparative and functional genomic analysis of prokaryotic nickel and cobalt uptake transporters: evidence for a novel group of ATP-binding cassette transporters. J Bacteriol. 2006;188(1):317–27.PubMed CentralView ArticlePubMedGoogle Scholar
- Markowitz VM, Chen IM, Palaniappan K, Chu K, Szeto E, Pillay M, et al. IMG 4 version of the integrated microbial genomes comparative analysis system. Nucleic Acids Res. 2014;42(Database issue):D560–7.PubMed CentralView ArticlePubMedGoogle Scholar
- Chiang PK, Gordon RK, Tal J, Zeng GC, Doctor BP, Pardhasaradhi K, et al. S-Adenosylmethionine and methylation. FASEB J. 1996;10(4):471–80.PubMedGoogle Scholar
- Kim J-Y, Suh J-W, Ji G-E. Evaluation of S-Adenosyl-L-Methionine Production by Bifidobacterium bifidum BGN4. Food Sci Biotechnol. 2008;17(1):184–7.Google Scholar
- Ciccarelli FD, Doerks T, von Mering C, Creevey CJ, Snel B, Bork P. Toward automatic reconstruction of a highly resolved tree of life. Science. 2006;311(5765):1283–7.View ArticlePubMedGoogle Scholar
- Scher JU, Sczesnak A, Longman RS, Segata N, Ubeda C, Bielski C, et al. Expansion of intestinal Prevotella copri correlates with enhanced susceptibility to arthritis. Elife (Cambridge). 2013;2:e01202.Google Scholar
- Hotamisligil GS, Erbay E. Nutrient sensing and inflammation in metabolic diseases. Nat Rev Immunol. 2008;8(12):923–34.PubMed CentralView ArticlePubMedGoogle Scholar
- Weisberg SP, McCann D, Desai M, Rosenbaum M, Leibel RL, Ferrante Jr AW. Obesity is associated with macrophage accumulation in adipose tissue. J Clin Invest. 2003;112(12):1796–808.PubMed CentralView ArticlePubMedGoogle Scholar
- Engstrom G, Hedblad B, Stavenow L, Lind P, Janzon L, Lindgarde F. Inflammation-sensitive plasma proteins are associated with future weight gain. Diabetes. 2003;52(8):2097–101.View ArticlePubMedGoogle Scholar
- Wellen KE, Hotamisligil GS. Inflammation, stress, and diabetes. J Clin Invest. 2005;115(5):1111–9.PubMed CentralView ArticlePubMedGoogle Scholar
- Cani PD, Amar J, Iglesias MA, Poggi M, Knauf C, Bastelica D, et al. Metabolic endotoxemia initiates obesity and insulin resistance. Diabetes. 2007;56(7):1761–72.View ArticlePubMedGoogle Scholar
- Xu H, Barnes GT, Yang Q, Tan G, Yang D, Chou CJ, et al. Chronic inflammation in fat plays a crucial role in the development of obesity-related insulin resistance. J Clin Invest. 2003;112(12):1821–30.PubMed CentralView ArticlePubMedGoogle Scholar
- Uysal KT, Wiesbrock SM, Marino MW, Hotamisligil GS. Protection from obesity-induced insulin resistance in mice lacking TNF-alpha function. Nature. 1997;389(6651):610–4.View ArticlePubMedGoogle Scholar
- Pickup JC. Inflammation and activated innate immunity in the pathogenesis of type 2 diabetes. Diabetes Care. 2004;27(3):813–23.View ArticlePubMedGoogle Scholar
- Kau AL, Ahern PP, Griffin NW, Goodman AL, Gordon JI. Human nutrition, the gut microbiome and the immune system. Nature. 2011;474(7351):327–36.PubMed CentralView ArticlePubMedGoogle 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.PubMed CentralView ArticlePubMedGoogle Scholar
- Franzosa EA, Morgan XC, Segata N, Waldron L, Reyes J, Earl AM, et al. Relating the metatranscriptome and metagenome of the human gut. Proc Natl Acad Sci U S A. 2014;111(22):E2329–38.PubMed CentralView ArticlePubMedGoogle Scholar