Gut microbiota in children with type 1 diabetes differs from that in healthy children: a case-control study
BMC Medicine volume 11, Article number: 46 (2013)
A recent study using a rat model found significant differences at the time of diabetes onset in the bacterial communities responsible for type 1 diabetes modulation. We hypothesized that type 1 diabetes in humans could also be linked to a specific gut microbiota. Our aim was to quantify and evaluate the difference in the composition of gut microbiota between children with type 1 diabetes and healthy children and to determine the possible relationship of the gut microbiota of children with type 1 diabetes with the glycemic level.
A case-control study was carried out with 16 children with type 1 diabetes and 16 healthy children. The fecal bacteria composition was investigated by polymerase chain reaction-denaturing gradient gel electrophoresis and real-time quantitative polymerase chain reaction.
The mean similarity index was 47.39% for the healthy children and 37.56% for the children with diabetes, whereas the intergroup similarity index was 26.69%. In the children with diabetes, the bacterial number of Actinobacteria and Firmicutes, and the Firmicutes to Bacteroidetes ratio were all significantly decreased, with the quantity of Bacteroidetes significantly increased with respect to healthy children. At the genus level, we found a significant increase in the number of Clostridium, Bacteroides and Veillonella and a significant decrease in the number of Lactobacillus, Bifidobacterium, Blautia coccoides/Eubacterium rectale group and Prevotella in the children with diabetes. We also found that the number of Bifidobacterium and Lactobacillus, and the Firmicutes to Bacteroidetes ratio correlated negatively and significantly with the plasma glucose level while the quantity of Clostridium correlated positively and significantly with the plasma glucose level in the diabetes group.
This is the first study showing that type 1 diabetes is associated with compositional changes in gut microbiota. The significant differences in the number of Bifidobacterium, Lactobacillus and Clostridium and in the Firmicutes to Bacteroidetes ratio observed between the two groups could be related to the glycemic level in the group with diabetes. Moreover, the quantity of bacteria essential to maintain gut integrity was significantly lower in the children with diabetes than the healthy children. These findings could be useful for developing strategies to control the development of type 1 diabetes by modifying the gut microbiota.
Type 1 diabetes is a worldwide problem, mainly in children, and it is associated with a significant burden, mostly related to the development of vascular complications . Type 1 diabetes is the result of a complex interaction between different degrees of genetic susceptibility and environmental factors [2–4]. The intestinal microbiota is one of these environmental factors currently under study, partly as a result of observations in both non-obese diabetic (NOD) mice and BioBreeding diabetes-prone rats, where the use of antibiotics was shown to prevent the onset of diabetes [5, 6]. Moreover, a recent study using NOD mice suggested that the development of type 1 diabetes can be prevented through modulation of the intestinal microbiota . Newly, Vaarala et al. suggested that the interaction between the intestinal environment, the barrier function and the immune system are crucial in the onset of type 1 diabetes . Using a rat model, Roesch et al. found significant differences at the time of diabetes onset in the bacterial communities responsible for type 1 diabetes modulation . Moreover, other studies have shown that beneficial bacteria, such as probiotic bacteria, have a protective effect in rodent models by delaying or preventing the onset of type 1 diabetes [9, 10]. With respect to mechanisms of action, Wen et al. found that the gut microbiome of NOD mice lacking an adaptor for multiple innate immune receptors responsible for recognizing microbial stimuli correlates with the disease onset, revealing a relationship between gut microbiota and the immune system . Recent studies have demonstrated that commensal bacteria are crucial for maturation and function of the mucosal immune system. The balance between two major effector T cell populations in the intestine, IL-17+ T helper 17 cells and Foxp3+ regulatory T cells, requires signals from commensal bacteria and is dependent on the composition of the intestinal microbiota [12–14]. In addition, increased gut permeability has been observed in patients with type 1 diabetes as well as in NOD mouse and BioBreeding rat models [15–18]. It has been suggested that this increased gut permeability (commonly called leaky gut) may affect the absorption of antigens that can attack and damage pancreatic beta cells . Because gut microbes can affect intestinal permeability, the gut ecology may play a role in the development of type 1 diabetes .
Only a few studies have evaluated the ecology of intestinal microbiota in autoimmune children who were not yet diabetic [21, 22]. These studies used a very low number of participants (four patients and four controls) and neither of them have controlled for such an important factor as the mode of delivery (natural birth or Cesarean) or the type and time of infant feeding (formula-fed or breast-fed), both of which determine the gut microbial composition during infancy [23, 24].
The aim of the present study, therefore, was to characterize the composition of fecal microbiota in children with type 1 diabetes as compared with children without diabetes (controlling for such factors as mode of delivery and breastfeeding time) using PCR-denaturing gradient gel electrophoresis (DGGE) and real-time quantitative PCR (qPCR) analysis. This was to determine whether there were significant differences in the gut microbiota composition between these groups and, if so, to quantify the differences and determine the possible relation of the gut microbiota of children with type 1 diabetes with their glycemic level.
Study participants and design
The case-control study included 16 Caucasian children with type 1 diabetes, aged 7.16 ±0.72 years, and 16 healthy Caucasian children, aged 7.48 ±0.87 years. Type 1 diabetes was diagnosed following the criteria of the American Diabetes Association [25, 26] and the appearance of at least two persistent, confirmed anti-islet autoantibodies (insulin autoantibodies, glutamic acid decarboxylase autoantibodies or tyrosine phosphatase autoantibodies). The patients with diabetes were treated and monitored according to a standard medical protocol. Patients were excluded if they had any other acute or chronic inflammatory diseases or infectious diseases at study entry. The study participants received no antibiotic treatment, probiotics, prebiotics or any other medical treatment influencing intestinal microbiota during the 3 months before the start of the study. The selected healthy children were all type 1 diabetes autoantibody negative and they were matched to the children with diabetes for age, gender, race, mode of delivery and duration of breastfeeding. The parents of the patients and controls completed a structured interview to obtain the following data: health status, lifestyle aspects (such as living environment and physical activity) and dietary habit. The dietary intake patterns in patients and controls were determined from a food frequency questionnaire that allowed us to assess the consumption of groups of foods. The written guardian or parental consents of the children were obtained. The sampling and experimental processes were performed with the approval of the local Ethics Committee of Ciudad de Jaen hospital. Stool samples were collected by parents at home and delivered to the storage area for frozen storage at -80ºC within one hour .
Body weight and height were measured according to standardized procedures .
Fasting venous blood samples were collected. The serum was separated in aliquots and immediately frozen at -80ºC. Serum biochemical parameters were measured in duplicate. Serum glucose, cholesterol and triglycerides were measured using a standard enzymatic method (Randox Laboratories Ltd., Antrim, UK). The quantitative detection of autoantibodies to islet cell antigens was done using the Elisa RSR GADAb Kit, Elisa RSR IA-2Ab Kit and RIA RSR IAA Kit (RSR Limited, Cardiff, UK).
DNA extraction from fecal samples
Fecal samples were immediately kept after collection at -80°C and stored until analyzed. DNA extraction from 200 mg of stools was done using the QIAamp DNA Stool Mini Kit (Qiagen, Hilden, Germany) following the manufacturer's instructions. The DNA concentration was determined by absorbance at 260 nm (A260), and the purity was estimated by determining the A260 to A280 ratio with a Nanodrop spectrophotometer (Nanodrop Technologies, Wilmington, DE, USA).
Analysis of fecal microbiota by PCR-DGGE
Fecal samples from each participant were examined by determining PCR-DGGE profiles as recently published by us . The V2 to V3 regions of the 16S rRNA genes (positions 339 to 539 in the Escherichia coli gene) of bacteria in the fecal samples were amplified by primers HDA1-GC (5′-CGC CCG CCG CGC GCG GCG GGC GGG GCG GGG GCA CGG GGG GCC TAC GGG AGG CAG CAG T-3′; (the GC clamp is in boldface)) and HDA2 (5′-GTA TTA CCG CGG CTG CTG GCA C-3′) generating a 200 bp product. Aliquots (2 μL) of DNA were amplified by real-time PCR (20 μL final volume) in a 7500 Fast Real-Time PCR Systems instrument using Fast SYBR Green Master Mix and 200 nM of each of the universal primers HDA1-GC or HDA2 with the following amplification program: initial denaturation at 95°C for 20 s; amplification using 45 cycles including denaturation at 95°C for 3 s; annealing at 55°C for 30 s; and extension at 72°C for 1 min. Negative controls without a DNA template were included in each analysis.
After real-time PCR, 15 μL of products were mixed with 6 μL of loading dye before loading. Electrophoresis was performed with a DCode Universal Mutation Detection System instrument (Bio-Rad Laboratories, S.A, Madrid, Spain). Six percent polyacrylamide gels were prepared and electrophoresed with 1× TAE buffer prepared from 50× TAE buffer (2 M Tris base, 1 M glacial acetic acid, 50 mM ethylenediaminetetraacetic acid (EDTA)). The denaturing gradient was formed by using two 6% acrylamide (acrylamide to bisacrylamide ratio 37.5:1) stock solutions (Bio-Rad). The gels contained a 20% to 80% gradient of urea and formamide that increases in the direction of electrophoresis. Electrophoretic runs were in a TAE buffer (40 mmol/L Tris, 20 mmol/L acetic acid, and 1 mmol/L EDTA, pH 7.4) at 130 V and 60°C for 4.5 h. Electrophoresis was stopped when a xylene cyanol dye marker reached the bottom of a gel. Gels were stained with ethidium bromide (0.5 mg/L) for 5 min, rinsed with deionized water, viewed by UV transillumination and photographed with Gelcapture image acquisition software (DNR Bio-Imaging Systems Ltd, Mahale HaHamisha, Jerusalen, Israel). All the samples were analyzed on the same DGGE run to avoid the possible influence of variations in electrophoretic conditions between different runs. No band was observed in the negative controls. Similarities between banding patterns in the DGGE profile were calculated based on the presence and absence of bands and expressed as a similarity coefficient. Gels were analyzed using BioNumerics software (Applied Maths, Sint-Martens-Latem, Belgium). Normalized banding patterns were used for cluster analysis. The Dice similarity coefficient was used to calculate pairwise comparisons of the DGGE fingerprint profiles obtained. A similarity coefficient value of 100% indicates that DGGE profiles are identical while completely different profiles result in a similarity coefficient value of 0%. The unweighted pair group method with arithmetic mean algorithm was used for construction of dendrograms.
Sequencing of selected bands from DGGE gels
Bands were excised from DGGE gels with a sterile razor, placed in 40 μL sterile water and incubated at 4°C for diffusion of DNA into the water. DNA were used in a second PCR with HDA1/2 primers without a GC-clamp (initial denaturation at 95°C for 20 s, followed by 45 cycles including denaturation at 95°C for 3 s, annealing at 55°C for 15 s and extension at 72°C for 10 s). Subsequently, the PCR products were directly cloned into pCR 4-TOPO (Invitrogen, Carlsbad, CA, USA) according to the manufacturer's instructions. Plasmid DNA was isolated from the cells using the Qiagen Mini Spin Prep Kit (Qiagen), and subjected to PCR (HDA1/2-GC) as earlier described. PCR products were diluted until 20 ng/μL, purified with ExoSAP-IT (USB Corporation, Cleveland, OH, USA) and sequenced in an ABI 3130 (Applied Biosystems, Inc., Foster City, CA, USA) using the BigDie-Kit-Standard. Nucleotide sequence data obtained were analyzed using MicroSeqID v2.1.1 software (Applied Biosystems).
Microbial quantification by real-time qPCR
Specific primers targeting different bacterial genera were used to characterize the fecal microbiota by quantitative real-time qPCR (Table 1) [30–36]. Briefly, quantitative PCR experiments were performed with a LightCycler 2.0 PCR sequence detection system using the FastStart DNA Master SYBR Green kit (Roche Diagnostics, Indianapolis, IN, USA). All PCR tests were carried out in duplicate, with a final volume of 20 μL containing 1 μL of each fecal DNA preparation and 200 nM of each primer (Table 1). The thermal cycling conditions used were as follows: an initial DNA denaturation step at 95°C for 10 min; 45 cycles of denaturation at 95°C for 10 s; primer annealing at optimal temperature for 20 s; extension at 72°C for 15 s. Finally, melt curve analysis was performed by slowly cooling the PCRs from 95°C to 60°C (0.05°C per cycle) with simultaneous measurement of the SYBR Green I signal intensity. Melting-point-determination analysis allowed the confirmation of the specificity of the amplification products. Each participant's extracted DNA was subjected to a human β-globin PCR to ensure that amplifiable DNA was successfully extracted from the sample and to monitor for PCR inhibitors with the reaction conditions described previously . The bacterial concentration from each sample was calculated by comparing the threshold cycle values obtained from the standard curves with the LightCycler 4.0 software. Standard curves were constructed for each experiment using serial 10-fold dilutions of bacterial genomic DNA (of known concentration) from pure cultures, corresponding to 101 to 1010 copies per gram of feces.
The different strains used were obtained from the Spanish Collection of Type Cultures (CECT) (Bacteroides vulgatus NCTC 11154, Fusobacterium varium NCTC 10560, Enterococcus faecalis CECT 184, Enterobacter cloacae CECT 194, Clostridium perfringens CECT 376) and the American Type Culture Collection (ATCC) (Bifidobacterium bifidum ATCC 15696, Lactobacillus casei ATCC 334D-5, Prevotella intermedia ATCC 25611D-5, Ruminococcus productus ATCC 27340D-5 and Veillonella dispar ATCC 17745). Standard curves were normalized to the copy number of the 16S rRNA gene for each species. For species for which the copy number of 16S rRNA operon was not published, the copy number was calculated by averaging the operon numbers of the closest bacterial taxa from the ribosomal RNA database rrnDB . Negative controls containing all the elements of the reaction mixture except template DNA were performed in every analysis and no product was ever detected. The data presented are the mean values of duplicate real-time qPCR analyses. The amplification efficiency of the qPCR for all primer pairs was determined using the linear regression slope of a dilution series based on the following equation E = 10(-1/slope). We found that for 13 primer pairs the efficiency ranged from 98% (E = 1.96) to 100% (E = 2) with slopes values in the range of -3.4 to -3.32.
Results are expressed as mean values and standard deviations. The statistical analysis was performed with SPSS 15.0 software (SPSS Inc., Chicago, IL, USA). The sample size was calculated to obtain a difference in the mean bacterial number between the healthy children and those with type 1 diabetes of at least 2 × 105 copies per gram of feces. With a power of 80%, an alpha error of 0.05 and an estimated standard deviation between group of 1.13 × 105 copies per gram of feces (data obtained from Wu et al. ), six children were needed in each group. However, we increased the number of participants to 16 children and 16 controls. The bacterial copy number values were converted into logarithmic values before the statistical analysis. Given the low number of participants analyzed, the Mann-Whitney U test was used to check changes in bacterial number and biochemical variables between the two groups. The Spearman correlation coefficient was calculated to estimate the linear correlations between variables. A multivariate regression analysis was performed to identify individual bacteria as independent predictors for plasma glucose level. Statistical significance was set at a P value of <0.05. All data are presented in the text as the mean ± SD.
The healthy children and those with diabetes all had similar physical activity and dietary habits. The analysis of the food frequency questionnaires showed no significant differences in the consumption patterns of rice, wheat, vegetables, fish or meat between the two study groups, although the children with diabetes had a fast carbohydrate restriction (foods made with white flour and refined sugar).
Anthropometric and biochemical measurements
The anthropometric and biochemical variables of the healthy children and those with diabetes are shown in Table 2. Apart from the levels of glucose and HbA1c, which were significantly higher in the children with diabetes, no other significant differences were seen between the groups in the anthropometric and biochemical variables. In addition, because the healthy children and the children with diabetes were matched for breastfeeding time and mode of delivery, no significant differences were noted in these variables.
PCR-DGGE and bacterial band identification
Variations were found in the presence or absence (qualitative) and intensity (quantitative) of the bands between the healthy children and the children with diabetes in the host-specific fingerprints generated. DGGE band profiles showed differences in band richness between the two groups. Analysis of the diversity of the microbiota showed that the mean of the DGGE bands was 13.85 ±3.87 for the healthy children and 11.63 ±3.64 for the children with diabetes, though the difference was not significant. Some bands were seen in fingerprints from all the children (in different lanes but at the same position), indicating that specific species of the predominant microbiota were common to all the children.
The Dice similarity coefficient was used to calculate the similarity index between DGGE band profiles related to sampling of healthy children and those with diabetes. The mean similarity index was 47.39% for the healthy children and 37.56% for the children with diabetes. The mean similarity index between the groups was 26.69%, lower than the intra-group similarity (Table 3). The DGGE gel and the results of the cluster analysis are shown in Figure 1. The cluster analysis showed that the intra-group similarity for the diabetic and the healthy groups was significantly higher than the inter-group similarity. These results demonstrate that the dominant microbiota in the healthy group was different from that of the diabetic group.
All the bands from the profiles of all the healthy children and the children with diabetes were cloned and sequenced to identify the dominant microbiota. The sequence similarity matches for bands were analyzed by MicroSeqID v2.1.1 software. Bacterial identification showed that the majority of bacteria represented in the fingerprints obtained corresponded to five phyla (Table 4). Most of the sequences belonged to Firmicutes and Bacteroidetes, with the rest distributed among Actinobacteria, Fusobacteria and Proteobacteria. Nevertheless, we observed important differences between the healthy children and the children with diabetes in the distribution ratio of the different genera within Bacteroidetes, Firmicutes and Actinobacteria phyla. In the children with diabetes, we found an increase in the Clostridium, Bacteroides, Veillonella, Eggerthella and Bacillus frequencies and a disappearance of Prevotella and Bifidobacterium as compared with the healthy children (Table 4).
Comparative analysis of gut microbiota communities in healthy children and children with diabetes
Changes in the bacterial population abundance were assessed in the fecal samples of both groups. The results obtained in the real-time qPCR experiments with the different primers are shown in Tables 5 and 6. Relevant differences were found in the bacteria number of three phyla between the diabetic and the healthy children. The number of Actinobacteria, Firmicutes and Bacteroidetes were significantly different between groups whilst the quantity of Proteobacteria and Fusobacteria were similar between the groups. In the children with diabetes, the bacterial number of Actinobacteria and Firmicutes was significantly decreased while that of Bacteroidetes was significantly increased with respect to the healthy children. Moreover, the Firmicutes to Bacteroidetes ratio was significantly lower in the children with diabetes than the healthy children.
Within Firmicutes, the quantity of Veillonella was significantly higher and the number of bacteria from the Blautia coccoides-Eubacterium rectale group was significantly lower in the children with diabetes compared with the healthy children. The Lactobacillus number was significantly lower and Clostridium levels significantly higher in the children with diabetes. However, no significant differences were found in Enterococcus levels between the two groups. Within Bacteroidetes, the quantity of Bacteroides was significantly higher whereas the number of Prevotella was significantly lower in the children with diabetes compared with the healthy children. Finally, within Actinobacteria, the number of Bifidobacterium was significantly lower in the children with diabetes.
Relationship between gutmicrobiota composition in children with type 1 diabetes and glycemic level
In the children with diabetes, we found a significant univariate correlation between the amount of specific bacterial groups and the plasma glucose levels (Bifidobacterium r = -0.797, P = 0.008; Clostridium r = 0.676, P <0.05; Lactobacillus r = -0.698, P <0.05; and Firmicutes to Bacteroidetes ratio r = -0.473, P <0.05) and HbA1c levels (Bifidobacterium r = -0.573, P <0.05; Clostridium r = 0.452, P <0.05; Firmicutes r = -0.559, P <0.05; Firmicutes to Bacteroidetes ratio r = -0.765, P = 0.012). A multivariate regression analysis that included all the bacterial groups analyzed showed that only the reduction in the number of Bifidobacterium and Lactobacillus was associated with the plasma glucose level (P <0.05, β = -0.476, R2 = 0.587; and P = 0.012, β = -0.687, R2 = 0.539, respectively) whereas the higher HbA1c level was associated with the decrease in the Firmicutes to Bacteroidetes ratio (P <0.001, β = -1.047, R2 = 0.781) and the increase in the number of Clostridium (P = 0.016, β = 0.867, R2 = 0.499).
In the present study we found significant differences in the fecal microbial composition between healthy children and children with type 1 diabetes. We are unaware of any other similar studies in children with type 1 diabetes using simultaneously DGGE molecular profiling, unweighted pair group method with arithmetic mean algorithm dendrogram construction, sequencing and real-time qPCR analysis. To determine the characteristics of the gut microbiota based on the condition of just type 1 diabetes, we excluded the influence of physiological factors such as age, gender, dietary habits and race. In addition, we also controlled for the mode of delivery at birth and the duration of breastfeeding. This was because the first year of life has a crucial impact on gut microbiota composition and epidemiological studies in humans at genetic risk for type 1 diabetes have suggested that a short duration of breastfeeding and early feeding in infancy with complex dietary proteins such as cow's milk proteins can modulate the development of beta cell autoimmunity, clinical type 1 diabetes, or both [40–42]. No significant differences were found between the two groups of children (type 1 diabetes and controls).
The DGGE analysis of the fecal microbiota revealed a significantly lower intra-group similarity index in children with diabetes than in healthy children. In other words, the DGGE profiles in healthy children were more similar to each other, whereas in children with diabetes they were less similar. A similar result was found by Giongo et al. . These data suggest that diabetic status may influence specific bacterial groups of the gut microbiota community.
Sequence analysis of the DGGE bands cloned enable the association of specific bacterial genotypes with health or diabetes situations. Consistent with previous human and animal studies [11, 21, 39, 37, 43], the gut microbiota of healthy children and children with diabetes was predominately composed of Firmicutes and Bacteroidetes and the main difference lies in the proportion of genus-division bacteria within this two phyla and the Actinobacteria phylum between both the group with diabetes and the healthy group. These results suggest that the dominant microbiota genera are different in children with type 1 diabetes compared with healthy children. Recently, three robust clusters, referred to as "enterotypes", which are not nation or continent specific have been identified. Assignment of an individual microbiome into a given enterotype is based upon the relative enrichment of that microbiome in one of three genera: Bacteroides (enterotype 1), Prevotella (enterotype 2) or Ruminococcus (enterotype 3) . In this study, within Bacteroidetes, the Bacteroides genus was prevalent in the diabetic group, whereas the Prevotella genus was associated with the healthy group. Thus, the type 1 diabetic gut microbiomes could be classified into enterotype 1 and the healthy microbiomes could be classified into enterotype 2.
As DGGE is considered a semiquantitative tool for monitoring the dynamics of the predominant bacterial species of fecal microbiota, an additional analysis with real-time qPCR was performed to obtain a quantitative estimation of the changes found in the gut microbiota between children with diabetes and healthy children. We noted significant quantitative differences between the major microbial phyla present in the feces of healthy children and those with diabetes. In contrast to the situation in healthy children, we found a significant increase in the quantity of Bacteroidetes and a significant decrease in the number of Firmicutes and Actinobacteria in children with type 1 diabetes. Our data showed a significantly lower Firmicutes to Bacteroidetes ratio in children with type 1 diabetes compared with healthy children. Moreover, we saw a negative correlation between this ratio and both the glucose and the HbA1C levels in children with diabetes, which could help to explain the significantly higher glycemic level in this group. In agreement with this, Giongo et al. observed that the Firmicutes to Bacteroidetes ratio in study participants with type 1 diabetes was changing during the first 6 months after birth before the development of the autoimmune disease. These authors showed a successive decline in Firmicutes and an increase in Bacteroidetes number in the gut microbiome over time until the children became diabetic . Moreover, this imbalance observed at the phylum level between Bacteroidetes and Firmicutes has been previously described in several human disorders. A decline in the Firmicutes to Bacteroidetes ratio compared with controls has been described in human type 2 diabetes , whereas in Crohn´s disease, both Bacteroidetes and Firmicutes seem to decline . The opposite happens in obesity, where the imbalance is due to the increase in the Firmicutes to Bacteroidetes ratio , indicating that obesity and diabetes are associated with different groups of intestinal microbiota.
However, the major difference between the two groups was found in the number of bacteria at genus-division level. The most remarkable result was the significant increase in the number of Clostridium, Bacteroides and Veillonella in the children with diabetes, whereas the number of Lactobacillus, Bifidobacterium, the Blautia coccoides/Eubacterium rectale group and Prevotella genus were all significantly decreased in children with diabetes. Our findings concerning the microbiota of children with diabetes are in line with observations in other animal studies. Roesch et al. found higher levels of Lactobacillus and Bifidobacterium in BioBreeding diabetes-resistant rats whereas Bacteroides and Clostridium were more abundant in BioBreeding diabetes-prone rats . In contrast with this, however, Brown et al. found that Lactobacillus and Bifidobacterium were more abundant in participants with type 1 diabetes than in healthy participants .
The significant decrease in the number of Lactobacillus and Bifidobacterium observed in children with type 1 diabetes in our study was associated with their higher levels of plasma glucose, as indicated by the negative correlation found. Also, the regression analysis showed that the decrease in the number of Lactobacillus and Bifidobacterium could be associated with the plasma glucose level in the children with diabetes. In previous studies, the levels of Bifidobacterium have also been related to improved glucose metabolism, insulin resistance and low-grade inflammation [48, 49]. Moreover, Valladares et al. determined that the administration of Lactobacillus johnsonii isolated from BioBreeding diabetes-resistant rats delays or inhibits the onset of type 1 diabetes in BioBreeding diabetes-prone rats .
Both Lactobacillus and Bifidobacterium have members with probiotic characteristics and these have been associated with positive effects for the host in the large intestine . In addition, both bacterial groups have the capacity to produce the beneficial organic acid lactate, which is converted into butyrate by butyrate-producing bacteria in the gut . Barcenilla et al.  showed that most of the butyrate-producing isolates from human fecal samples are related to the Blautia coccoides-Eubacterium rectale group. Previous studies have shown that butyrate induces mucin synthesis (a glycoprotein produced by the host that could maintain the integrity of the gut epithelium) , decreases bacterial transport across the epithelium , and improves gut integrity by increasing tight junction assembly . In addition, the genera Prevotella are responsible for the degradation of this mucin ; thus, the significant decline in the numbers of the Blautia coccoides-Eubacterium rectale group and Prevotella that we found in children with type 1 diabetes compared with healthy children could indicate a reduction in mucin synthesis by the host and a lack of this mucin on the epithelial layer of the gut, which would lead to a significant alteration in intestinal permeability. Other studies have described an association between type 1 diabetes and compromised barrier permeability in humans and both the NOD mouse and BioBreeding rat models [16–18, 20].
The significant increase in the number of Clostridium, Bacteroides and Veillonella in the children with diabetes with respect to the healthy children was accompanied by a significant positive correlation between both the plasma levels of glucose and HbA1c and the quantity of Clostridium. These bacteria are able to ferment glucose and lactate to propionate, acetate and succinate. However, these short fatty acids do not induce mucin synthesis . This situation would, though, reduce the tight junction assembly, generating an increase in the gut permeability in children with type 1 diabetes .
Finally, we propose a possible mechanism to explain the relationship we have found between the gut microbiota present in children with type 1 diabetes and the glycemic levels observed. The short-chain fatty acids (such as butyrate and propionate) formed by this gut microbiota have a role in the regulation of the levels of gut hormones such as glucose-dependent insulinotropic polypeptide, glucagon-like peptide 1 and ghrelin. These hormones have important effects on carbohydrate metabolism , thus allowing gut microbiota to affect glycemic levels. In addition, Huml et al. have previously demonstrated an altered secretion pattern of gut hormones in children with type 1 diabetes that may impact on the metabolic control of diabetes in these patients . Further studies will be necessary to demonstrate this proposed mechanism.
A limitation of the 16S rRNA gene-based method is that the function of the identified bacteria is unknown. Future studies using a microbial metagenomic sequencing analysis will be carried out to obtain information about the functional diversity of the bacterial community analyzed here.
This is the first study showing that type 1 diabetes is associated with compositional changes in gut microbiota. Our results show that gut microbiota found in children with type 1 diabetes differed significantly from that found in healthy children. The gut microbiota in the children with diabetes was less similar than the gut microbiota in the healthy children. The significant differences between the diabetic and the healthy children in the number of Bifidobacterium, Lactobacillus and Clostridium and the Firmicutes to Bacteroidetes ratio could be implicated in the glycemic level of the children with diabetes. In addition, the numbers of lactic acid-producing bacteria, butyrate-producing bacteria and mucin-degrading bacteria, essential to maintain gut integrity, were significantly lower in the children with diabetes than the healthy children. These bacterial differences could be responsible for the altered gut permeability previously described in patients with type 1 diabetes. These findings could be useful for developing strategies to control the development of type 1 diabetes by modifying the gut microbiota.
denaturing gradient gel electrophoresis
non-obese diabetic mice
polymerase chain reaction
quantitative polymerase chain reaction
buffer with Tris base, glacial acetic acid and ethylenediaminetetraacetic acid.
Marcovecchio ML, Tossavainen PH, Dunger DB: Prevention and treatment of microvascular disease in childhood type 1 diabetes. Br Med Bull. 2010, 94: 145-164. 10.1093/bmb/ldp053.
Ehehalt S, Dietz K, Willasch AM, Neu A, Baden-Württemberg : Diabetes Incidence Registry (DIARY) Group. Epidemiological perspectives on type 1 diabetes in childhood and adolescence in Germany: 20 years of the Baden-Württemberg Diabetes Incidence Registry (DIARY). Diabetes Care. 2010, 3: 338-340.
Patterson CC, Dahlquist G, Soltesz G, Green A, Grp EAS: Is childhood onset Type I diabetes a wealth-related disease? An ecological analysis of European incidence rates. Diabetologia. 2001, 4: 9-16.
Vaarala O, Atkinson MA, Neu J: The ''perfect storm'' for type 1 diabetes - the complex interplay between intestinal microbiota, gut permeability, and mucosal immunity. Diabetes. 2008, 57: 2555-2562. 10.2337/db08-0331.
Brugman S, Klatter FA, Visser JT, Wildeboer-Veloo AC, Harmsen HJ, Rozing J, Bos NA: Antibiotic treatment partially protects against type 1 diabetes in the Bio-Breeding diabetes-prone rat. Is the gut flora involved in the development of type 1 diabetes?. Diabetologia. 2006, 49: 2105-2108. 10.1007/s00125-006-0334-0.
Schwartz RF, Neu J, Schatz D, Atkinson MA, Wasserfall C: Comment on: Brugman S et al. (2006) Antibiotic treatment partially protects against type 1 diabetes in the Bio-Breeding diabetes-prone rat. Is the gut flora involved in the development of type 1 diabetes? Diabetologia 49:2105-2108. Diabetologia. 2007, 50: 220-221.
King C, Sarvetnick N: The incidence of type-1 diabetes in NOD mice is modulated by restricted flora not germ-free conditions. Plos One. 2011, 6: e17049-10.1371/journal.pone.0017049.
Roesch LFW, Lorca GL, Casella G, Giongo A, Naranjo A, Pionzio AM, Li N, Mai V, Wasserfall CH, Schatz D, Atkinson MA, Neu J, Triplett EW: Culture-independent identification of gut bacteria correlated with the onset of diabetes in a rat model. Isme J. 2009, 3: 536-548. 10.1038/ismej.2009.5.
Lai KK, Lorca GL, Gonzalez CF: Biochemical properties of two cinnamoyl esterases purified from a Lactobacillus johnsonii strain isolated from stool samples of diabetes-resistant rats. Appl Environ Microbiol. 2009, 75: 5018-5024. 10.1128/AEM.02837-08.
Valladares R, Sankar D, Li N, Williams E, Lai KK, Abdelgeliel AS, Gonzalez CF, Wasserfall CH, Larkin J, Schatz D, Atkinson MA, Triplett EW, Neu J, Lorca GL: Lactobacillus johnsonii N6.2 mitigates the development of type 1 diabetes in BB-DP rats. PLoS One. 2010, 5: e10507-10.1371/journal.pone.0010507.
Wen L, Ley RE, Volchkov PY, Stranges PB, Avanesyan L, Stonebraker AC, Hu C, Wong FS, Szot GL, Bluestone JA, Gordon JI, Chervonsky AV: Innate immunity and intestinal microbiota in the development of Type 1 diabetes. Nature. 2008, 45: 1109-1113.
Ivanov II, Atarashi K, Manel N, Brodie EL, Shima T, Karaoz U, Wei D, Goldfarb KC, Santee CA, Lynch SV, Tanoue T, Imaoka A, Itoh K, Takeda K, Umesaki Y, Honda K, Littman DR: Induction of intestinal Th17 cells by segmented filamentous bacteria. Cell. 2009, 139: 485-498. 10.1016/j.cell.2009.09.033.
Ivanov II, Littman DR: Segmented filamentous bacteria take the stage. Mucosal Immunol. 2010, 3: 209-212. 10.1038/mi.2010.3.
Romano-Keeler J, Weitkamp JH, Moore DJ: Regulatory properties of the intestinal microbiome effecting the development and treatment of diabetes. Curr Opin Endocrinol Diabetes Obes. 2012, 19: 73-80.
Kuitunen M, Saukkonen T, Ilonen J, Akerblom HK, Savilahti E: Intestinal permeability to mannitol and lactulose in children with type 1 diabetes with the HLA-DQB1*02 allele. Autoimmunity. 2002, 35: 365-368. 10.1080/0891693021000008526.
Watts T, Berti I, Sapone A, Gerarduzzi T, Not T, Zielke R, Fasano A: Role of the intestinal tight junction modulator zonulin in the pathogenesis of type I diabetes in BB diabetic-prone rats. Proc Natl Acad Sci. 2005, 102: 2916-2921. 10.1073/pnas.0500178102.
Bosi E, Molteni L, Radaelli MG, Folini L, Fermo I, Bazzigaluppi E, Piemonti L, Pastore MR, Paroni R: Increased intestinal permeability precedes clinical onset of type 1 diabetes. Diabetologia. 2006, 49: 2824-2827. 10.1007/s00125-006-0465-3.
Lee AS, Gibson DL, Zhang Y, Sham HP, Vallance BA, Dutz JP: Gut barrier disruption by an enteric bacterial pathogen accelerates insulitis in NOD mice. Diabetologia. 2010, 53: 741-748. 10.1007/s00125-009-1626-y.
Vehik K, Dabelea D: The changing epidemiology of type 1 diabetes: why is it going through the roof?. Diabetes Metab Res Rev. 2011, 27: 3-13. 10.1002/dmrr.1141.
Mathis D, Benoist C: The influence of the microbiota on type-1 diabetes: on the threshold of a leap forward in our understanding. Immunol Rev. 2012, 245: 239-249. 10.1111/j.1600-065X.2011.01084.x.
Giongo A, Gano KA, Crabb DB, Mukherjee N, Novelo LL, Casella G, Drew JC, Ilonen J, Knip M, Hyöty H, Veijola R, Simell T, Simell O, Neu J, Wasserfall CH, Schatz D, Atkinson MA, Triplett EW: Toward defining the autoimmune microbiome for type 1 diabetes. Isme J. 2011, 5: 82-91. 10.1038/ismej.2010.92.
Brown CT, Davis-Richardson AG, Giongo A, Gano KA, Crabb DB, Mukherjee N, Casella G, Drew JC, Ilonen J, Knip M, Hyöty H, Veijola R, Simell T, Simell O, Neu J, Wasserfall CH, Schatz D, Atkinson MA, Triplett EW: Gut microbiome metagenomics analysis suggests a functional model for the development of autoimmunity for type 1 diabetes. PLoS One. 2011, 6: e25792-10.1371/journal.pone.0025792.
Penders J, Thijs C, Vink C, Stelma FF, Snijders B, Kummeling I, van den Brandt PA, Stobberingh EE: Factors influencing the composition of the intestinal microbiota in early infancy. Pediatrics. 2006, 118: 511-521. 10.1542/peds.2005-2824.
Musso G, Gambino R, Cassader M: Obesity, diabetes, and gut microbiota: the hygiene hypothesis expanded?. Diabetes Care. 2010, 3: 2277-2284.
Report of the Expert Committee on the diagnosis and classification of diabetes mellitus. Diabetes Care. 1997, 20: 1183-1197.
American Diabetes Association: Diagnosis and classification of diabetes mellitus. Diabetes Care. 2010, 3: S62.
Roesch LF, Casella G, Simell O, Krischer J, Wasserfall CH, Schatz D, Atkinson MA, Neu J, Triplett EW: Influence of sample storage on bacterial community diversity in fecal samples. Open Microbiol J. 2009, 3: 40-46. 10.2174/1874285800903010040.
Standardization of anthropometric measurements. The Airlie (VA) Consensus Conference. Edited by: Loham T, Roche A, Martorel R. 1988, Champaign, IL: Human Kinetics, 20-37.
Queipo-Ortuño MI, Boto-Ordóñez M, Murri M, Gomez-Zumaquero JM, Clemente-Postigo M, Estruch R, Cardona Diaz F, Andrés-Lacueva C, Tinahones FJ: Influence of red wine polyphenols and ethanol on the gut microbiota ecology and biochemical biomarkers. Am J Clin Nutr. 2012, 95: 1323-1334. 10.3945/ajcn.111.027847.
Guo X, Xia X, Tang R, Zhou J, Zhao H, Wang K: Development of a real-time PCR method for Firmicutes and Bacteroidetes in faeces and its application to quantify intestinal population of obese and lean pigs. Lett Appl Microbiol. 2008, 47: 367-373. 10.1111/j.1472-765X.2008.02408.x.
Delroisse JM, Boulvin AL, Parmentier I, Dauphin RD, Vandenbol M, Portetelle D: Quantification of Bifidobacterium spp. and Lactobacillus spp. in rat fecal samples by real-time PCR. Microbiol Res. 2008, 163: 663-670. 10.1016/j.micres.2006.09.004.
Friswell MK, Gika H, Stratford IJ, Theodoridis G, Telfer B, Wilson ID, McBain AJ: Site and strain-specific variation in gut microbiota profiles and metabolism in experimental mice. PLoS One. 2010, 5: e8584-10.1371/journal.pone.0008584.
Stach JE, Maldonado LA, Ward AC, Goodfellow M, Bull AT: New primers for the class Actinobacteria: application to marine and terrestrial environments. Environ Microbiol. 2003, 5: 828-841. 10.1046/j.1462-2920.2003.00483.x.
Matsuki T, Watanabe K, Fujimoto J, Takada T, Tanaka R: Use of 16S rRNA gene-targeted group-specific primers for real-time PCR analysis of predominant bacteria in human feces. Appl Environ Microbiol. 2004, 70: 7220-7228. 10.1128/AEM.70.12.7220-7228.2004.
Bekele AZ, Koike S, Kobayashi Y: Genetic diversity and diet specificity of ruminal Prevotella revealed by 16S rRNA gene-based analysis. FEMS Microbiol Lett. 2010, 305: 49-57. 10.1111/j.1574-6968.2010.01911.x.
Rinttilä T, Kassinen A, Malinen E, Krogius L, Palva A: Development of an extensive set of 16S rDNA-targeted primers for quantification of pathogenic and indigenous bacteria in faecal samples by real-time PCR. J Appl Microbiol. 2004, 97: 1166-1177. 10.1111/j.1365-2672.2004.02409.x.
Fredricks DN, Fiedler TL, Thomas KK, Oakley BB, Marrazzo JM: Targeted PCR for detection of vaginal bacteria associated with bacterial vaginosis. J Clin Microbiol. 2007, 45: 3270-3276. 10.1128/JCM.01272-07.
Lee ZMP, Bussema C, Schmidt TM: rrnDB: documenting the number of rRNA and tRNA genes in bacteria and archaea. Nucleic Acids Res. 2009, 37: 489-493. 10.1093/nar/gkn689.
Wu X, Ma C, Han L, Nawaz M, Gao F, Zhang X, Yu P, Zhao C, Li L, Zhou A, Wang J, Moore JE, Millar BC, Xu J: Molecular characterisation of the faecal microbiota in patients with type II diabetes. Curr Microbiol. 2010, 61: 69-78. 10.1007/s00284-010-9582-9.
Stuebe A: The risks of not breastfeeding for mothers and infants. Rev Obstet Gynecol. 2009, 2: 222-231.
Vaarala O: The gut as a regulator of early inflammation in type 1 diabetes. Curr Opin Endocrinol Diabetes Obes. 2011, 18: 241-247. 10.1097/MED.0b013e3283488218.
Knip M, Virtanen SM, Becker D, Dupré J, Krischer JP, Åkerblom HK, TRIGR Study Group: Early feeding and risk of type 1 diabetes: experiences from the Trial to Reduce Insulin-dependent diabetes mellitus in the Genetically at Risk (TRIGR). Am J Clin Nutr. 2011, 94: 1814-1820. 10.3945/ajcn.110.000711.
Gill SR, Pop M, Deboy RT, Eckburg PB, Turnbaugh PJ, Samuel BS, Gordon JI, Relman DA, Fraser-Liggett CM, Nelson KE: Metagenomic analysis of the human distal gut microbiome. Science. 2006, 312: 1355-1359. 10.1126/science.1124234.
Arumugam M, Raes J, Pelletier E, Le Paslier D, Yamada T, Mende DR, Fernandes GR, Tap J, Bruls T, Batto JM, Bertalan M, Borruel N, Casellas F, Fernandez L, Gautier L, Hansen T, Hattori M, Hayashi T, Kleerebezem M, Kurokawa K, Leclerc M, Levenez F, Manichanh C, Nielsen HB, Nielsen T, Pons N, Poulain J, Qin J, Sicheritz-Ponten T, Tims S, et al: Enterotypes of the human gut microbiome. Nature. 2011, 473: 174-180. 10.1038/nature09944.
Larsen N, Vogensen FK, van den Beg FWL, Nielsen DS, Andreasen AS, Pedersen BK, Al-Soud WA, Sørensen SJ, Hansen LH, Jakobsen M: Gut microbiota in human adults with type 2 diabetes differs from non-diabetic adults. PLoS One. 2010, 5: e9085-10.1371/journal.pone.0009085.
Willing B, Halfvarson J, Dicksved J, Rosenquist M, Järnerot G, Engstrand L, Tysk C, Jansson JK: Twin studies reveal specific imbalances in the mucosa-associated microbiota of patients with ileal Crohn's disease. Inflamm Bowel Dis. 2009, 15: 653-660. 10.1002/ibd.20783.
Armougom F, Henry M, Vialettes B, Raccah D, Raoult D: Monitoring bacterial community of human gut microbiota reveals an increase in Lactobacillus in obese patients and methanogens in anorexic patients. PLoS One. 2009, 4: e7125-10.1371/journal.pone.0007125.
Philippe D, Favre L, Foata F, Adolfsson O, Perruisseau-Carrier G, Vidal K, Reuteler G, Dayer-Schneider J, Mueller C, Blum S: Bifidobacterium lactis attenuates onset of inflammation in a murine model of colitis. World J Gastroenterol. 2011, 17: 459-469. 10.3748/wjg.v17.i4.459.
Cani PD, Possemiers S, Van de Wiele T, Guiot Y, Everard A, Rottier O, Geurts L, Naslain D, Neyrinck A, Lambert DM, Muccioli GG, Delzenne NM: Changes in gut microbiota control inflammation in obese mice through a mechanism involving GLP-2-driven improvement of gut permeability. Gut. 2009, 58: 1091-1103. 10.1136/gut.2008.165886.
Tzounis X, Rodriguez-Mateos A, Vulevic J, Gibson GR, Kwik-Uribe C, Spencer JP: Prebiotic evaluation of cocoa-derived flavanols in healthy humans by using a randomized, controlled, double-blind, crossover intervention study. Am J Clin Nutr. 2011, 93: 62-72. 10.3945/ajcn.110.000075.
Barcenilla A, Pryde SE, Martin JC, Duncan SH, Stewart CS, Henderson C, Flint HJ: Phylogenetic relationships of butyrate-producing bacteria from the human gut. Appl Environ Microbiol. 2000, 6: 1654-1661.
Burger-van Paassen N, Vincent A, Puiman PJ, van der Sluis M, Bouma J, Boehm G, van Goudoever JB, van Seuningen I, Renes IB: The regulation of intestinal mucin MUC2 expression by short-chain fatty acids: implications for epithelial protection. Biochem J. 2009, 420: 211-219. 10.1042/BJ20082222.
Lewis K, Lutgendorff F, Phan V, Soderholm JD, Sherman PM, McKay DM: Enhanced translocation of bacteria across metabolically stressed epithelia is reduced by butyrate. Inflamm Bowel Dis. 2010, 16: 1138-1148. 10.1002/ibd.21177.
Peng LY, Li Z, Green RS, Holzman IR, Lin J: Butyrate enhances the intestinal barrier by facilitating tight junction assembly via activation of AMP-activated protein kinase in Caco-2 cell monolayers. J Nutr. 2009, 139: 1619-1625. 10.3945/jn.109.104638.
Wright DP, Knight CG, Parker SG, Christie DL, Roberton AM: Cloning of a mucin-desulfating sulfatase gene from Prevotella strain RS2 and its expression using a Bacteroides recombinant system. J Bacteriol. 2000, 182: 3002-3007. 10.1128/JB.182.11.3002-3007.2000.
Lin HV, Frassetto A, Kowalik EJ, Nawrocki AR, Lu MM, Kosinski JR, Hubert JA, Szeto D, Yao X, Forrest G, Marsh DJ: Butyrate and propionate protect against diet-induced obesity and regulate gut hormones via free fatty acid receptor 3-independent mechanisms. PLoS One. 2012, 7: e35240-10.1371/journal.pone.0035240.
Huml M, Kobr J, Siala K, Varvařovská J, Pomahačová R, Karlíková M, Sýkora J: Gut peptide hormones and pediatric type 1 diabetes mellitus. Physiol Res. 2011, 60: 647-658.
The pre-publication history for this paper can be accessed here:http://www.biomedcentral.com/1741-7015/11/46/prepub
The authors wish to thank all the participants for their collaboration and the Plataforma de Secuenciacion y Genotipado and the Unidad Central de Biología Molecular (Servicio de PCR a tiempo real) of FIMABIS for their help in laboratory assays. We also gratefully acknowledge the help of Ian Johnstone for his expertise in preparing this manuscript. The research group belongs to the Centros de Investigación en Red (CIBEROBN, CB06/03/0018 and CIBERDEM, CB07/08/0019) of the Instituto de Salud Carlos III. This work was partially funded by a grant from CIBER, CB06/03/0018 of the Instituto de Salud Carlos III to MM, FC, FJT and MIQO; the Instituto de Salud Carlos III, Madrid, Spain (CP07/0095) to FJT; and the Servicio Andaluz de Salud, Andalucía, Spain (PI0696/2010) to FJT. The funding agencies had no role in the design and performance of the study, the interpretation of the data or the writing of the manuscript.
The authors declare that they have no competing interests.
IL, FC, FS, FJT and MIQO conceived the study and developed the experimental design. IL, FS, FJT and MIQO were responsible for acquisition and selection of all samples utilized in this study. MM, IL, FC, JMGZ and MIQO performed all laboratory assays. MM, IL, JMGZ, FC and MIQO compiled the database and performed statistical analysis and data interpretation. MM, IL, JMGZ, FC, FS, FJT and MIQO wrote the paper. FS, FJT and MIQO provided critical revision. All the authors have read and approved the final manuscript.
Mora Murri, Isabel Leiva contributed equally to this work.
Authors’ original submitted files for images
Below are the links to the authors’ original submitted files for images.
About this article
Cite this article
Murri, M., Leiva, I., Gomez-Zumaquero, J.M. et al. Gut microbiota in children with type 1 diabetes differs from that in healthy children: a case-control study. BMC Med 11, 46 (2013). https://doi.org/10.1186/1741-7015-11-46
- butyrate-producing bacteria
- glycemic level
- gut integrity
- gut microbiota
- gut permeability
- HbA1c level
- lactic acid-producing bacteria
- mode of delivery
- type 1 diabetes