Sub-clinical detection of gut microbial biomarkers of obesity and type 2 diabetes
Genome Medicine volume 8, Article number: 17 (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
We collected fecal samples from 20 MZ Korean twins (10 twin pairs) at up to two time points each (12-44 months apart), resulting in 36 samples sequenced using metagenomic shotgun sequencing. In addition, multiple clinical parameters were measured at each sampling point, including body mass index (BMI), fasting blood sugar (FBS), cholesterol levels (LDL, HDL), fasting blood insulin (FBI), and renal and liver function (Additional file 1: Table S1), spanning the typical healthy range of these variables (Additional file 2: Figure S1). Species-level microbial abundance profiles were inferred using MetaPhlAn , and functional gene and pathway abundance profiles were generated using HUMAnN . These microbial and functional profiles were tested for statistically significant association with clinical parameters using the MaAsLin  sparse multivariate linear model (Fig. 1a).
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
Several organisms were prevalent (present in >50 % of individuals) in this cohort, although often at relatively low abundance levels. Some of these are shared with other globally surveyed populations, while others were unusually prevalent in this population (Fig. 2a). Shared organisms included Eubacterium rectale, Roseburia intestinalis, and Faecalibacterium prausnitzii, which are similarly prevalent in our cohort and in Western population (94-96 %, 83-89 %, and 96-97 %, respectively) and with similar relative abundances (5 %, approximately 1 %, and 2-4 %, respectively), confirming the similarity between this cohort and Western population.
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
We investigated the relationship between host clinical phenotype and the gut microbiome by identifying significant multivariate linear associations using MaAsLin  (Fig. 3, Additional file 6: Table S3). This model associates microbial clade or pathway abundances with metadata of interest (for example, BMI, FBS, triglyceride) while accounting for other covariates (in this case sex, smoking, age, and the twin pairing; see Methods and Additional file 7: Figure S4). The abundance of Akkermansia muciniphila was negatively correlated with BMI, FBS, and insulin levels, for example, all in gradients ranging continuously over the ranges of these clinical variables and its relative abundance. This mucin-degrading microbe has been observed to be reduced in the guts of obese mice , pregnant women [10, 11], and overweight children , but this is the first time this trend has been observed in non-pregnant adults, especially within the normal range of BMI and FBS. This suggests the organism may represent one aspect of the obese gut microbiome that may be of specifically sub-clinical significance.
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
To better understand the relation between the taxonomic compositional profiles and function, we correlated the profiles of each clade with the profiles of each module (Fig. 4, selected individual scatter plots appear in Additional file 9: Figure S6). This allows us to hypothesize which microbes contribute to, depend on, or associate with specific metabolic and biomolecular processes carried out by the gut microbiota. In particular, a positive correlation between a module and microbe can have two explanations. In some cases, the function may be encoded in the microbe’s genome (referred to as ‘directly encoded’ correlations, Additional file 10: Figure S7A). Alternatively, a microbe might correlate with a function not because it carries it itself, but because it associates with other microbes that encode it (‘indirectly associated’ correlations). The former indicates microbes that perform a particular molecular process, the latter those that depend on its presence elsewhere in the community.
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.
Remarkably, twins were not significantly more similar than unrelated in their strain composition (P = 0.15 by t-test), although (as expected) samples from the same individual over time were significantly similar (P = 1.34e-7; Fig. 5a). This suggests that either there is a genetic tendency for twins to retain broadly similar microbial compositions - but that this does not extend exert selective pressure at the strain level - or that identical strains acquired earlier in life during colonization have, by adulthood, evolved sufficiently to differ at multiple genomic elements. Only in rare cases did strains differ within individuals over time (Additional file 13: Figure S10), concordant with occasional sweeps of a replacement strain due to, for example, gene acquisition/loss or transfer from an external source.
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.
This pattern of microbial functional enrichments during inflammation has now been consistent across multiple studies regardless of their potential causality, including in other diseases like IBD (Fig. 6), suggesting a common signature of gut response to low-grade inflammation. Interestingly, several studies have examined whether inflammation can lead to obesity [45–47] and T2D [48–50] in mice, finding that inflammation drives the development of insulin resistance (potentially through the phosphorylation of insulin receptor 1 by a TNF-α mediated response ) and suggesting that particular intestinal microbial configurations can promote immune responses driving metabolic dysfunction . The extent to which the gut microbiota causes obesity is an area of active research. Many mouse models, including Lep ob, consistently demonstrate an elevated Firmicutes:Bacteroidetes (F:B) ratio in obese animals [5, 31]. In contrast, in human cohorts, the relationship has been much less consistent. The relationship between obesity and the F:B ratio has been reported as increased [5, 8], decreased , and others have reported no relationship [21, 52], indicating there is still great variability in current studies. This may arise either from technical issues, like different sample or data handling protocols, or from biological reasons, like true variation between the various cohorts.
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.
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.
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.
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.
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.
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.
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.
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.
Ley RE, Turnbaugh PJ, Klein S, Gordon JI. Microbial ecology: human gut microbes associated with obesity. Nature. 2006;444(7122):1022–3.
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.
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.
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.
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.
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.
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.
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.
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.
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.
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.
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.
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.
Human Microbiome Project Consortium. Structure, function and diversity of the healthy human microbiome. Nature. 2012;486(7402):207–14.
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.
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.
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.
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.
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.
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.
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.
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.
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.
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.
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.
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.
Ley RE. Obesity and the human microbiome. Curr Opin Gastroenterol. 2010;26(1):5–11.
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.
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.
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.
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.
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.
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.
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.
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.
Hotamisligil GS, Erbay E. Nutrient sensing and inflammation in metabolic diseases. Nat Rev Immunol. 2008;8(12):923–34.
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.
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.
Wellen KE, Hotamisligil GS. Inflammation, stress, and diabetes. J Clin Invest. 2005;115(5):1111–9.
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.
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.
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.
Pickup JC. Inflammation and activated innate immunity in the pathogenesis of type 2 diabetes. Diabetes Care. 2004;27(3):813–23.
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.
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.
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.
The authors declare that they have no competing interests.
DG, RJX, BWB, GPK, and CH initiated and designed the study. MYL, HSY, JS, YMS, KL, and GPK acquired the data. MY, TTT, EAF, XCM, DG, ESL, RJX, GPK, and CH analyzed and interpret the data. MY, XCM, ESL, RJX, GPK, and CH wrote or revised the paper. All authors read and approved the final manuscript.
The clinical markers collected at each time a sample was taken. (XLSX 36 kb)
Distributions of selected clinical variables in our data (colorful) and as collected by the human microbiome project  (gray). (PDF 3586 kb)
Prevalence and relative abundance of microbes and KEGG modules found in both this study and HMP . (XLSX 31 kb)
The average microbial prevalence (A) and abundance (B) as compared between the Human Microbiome Project  and our data. 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. (PDF 3577 kb)
Hierarchical clustering of the (A) phylogenetic and (B) functional profiles of the samples. Profiles (rows) are colored by the cluster assignment (k = 15), and samples (columns) are colored by the twin variable. The medoid profile per cluster (written in black) was chosen as a representative and shown on Fig. 2. (PDF 3597 kb)
All significant associations between the clinical markers and microbes and KEGG modules found in our data. (XLSX 35 kb)
Correlation between selected clinical variables in our data. (PDF 3598 kb)
Scatter plots of selected associations of microbial (x-axis) and functional (y-axis) profiles. Each plot represents a single entry in the correlation matrix at Fig. 4, and the R^2 value of this correlation dictates the color of the matrix entry. (PDF 3643 kb)
The genomic contribution of the microbial-to-functional associations. (A) Each entry denotes the fraction of genomic reference sequences of a specific microbe (column) that has genes from a specific functional module (row). This matrix accounts for most of the significant correlations we find on Fig. 4. (B) Taxa co-occurrence matrix in our data, using Spearman correlation. (ZIP 6726 kb)
About this article
Cite this article
Yassour, M., Lim, M.Y., Yun, H.S. et al. Sub-clinical detection of gut microbial biomarkers of obesity and type 2 diabetes. Genome Med 8, 17 (2016). https://doi.org/10.1186/s13073-016-0271-6
- Fast Blood Sugar
- Human Microbiome Project
- Bifidobacterium Longum
- Abundance Profile
- Riboflavin Biosynthesis