People with type 1 diabetes are at risk for developing micro- and macrovascular complications. Little is known about the gut microbiome in long-standing type 1 diabetes. We explored differences in the gut microbiome of participants with type 1 diabetes compared with healthy control subjects and associated the gut microbiome with diabetes-related complications.
Microbiome data of 238 participants with type 1 diabetes with an average disease duration of 28 ± 15 years were compared with 2,937 age-, sex-, and BMI-matched individuals. Clinical characteristics and fecal samples were collected, and metagenomic shotgun sequencing was performed. Microbial taxonomy was associated with type 1 diabetes–related characteristics and vascular complications.
No significant difference in the α-diversity of the gut microbiome was found between participants with type 1 diabetes and healthy control subjects. However, 43 bacterial taxa were significantly depleted in type 1 diabetes, while 37 bacterial taxa were significantly enriched. HbA1c and disease duration explained a significant part of the variation in the gut microbiome (R2 > 0.008, false discovery rate [FDR] <0.05), and HbA1c was significantly associated with the abundance of several microbial species. Additionally, both micro- and macrovascular complications explained a significant part of the variation in the gut microbiome (R2 > 0.0075, FDR < 0.05). Nephropathy was strongly associated with several microbial species. Macrovascular complications displayed similar associations with nephropathy.
Our data show that the gut microbiome is altered in people with (long-standing) type 1 diabetes and is associated with glycemic control and diabetes-related complications. As a result of the cross-sectional design, the causality of these relationships remains to be determined.
Introduction
Chronic hyperglycemia in type 1 diabetes drives the development of various microvascular complications, including retinopathy, neuropathy, and nephropathy (1), and accelerates macrovascular complications consisting of coronary artery disease, peripheral arterial disease, and stroke (2).
Several studies have suggested that the gut microbiome plays a role in the development of immune-mediated diseases, including type 1 diabetes (3–5). These include many cross-sectional studies in humans (4,6) and some studies showing that microbial changes are present before the onset of an immune-mediated disease (7). The gut microbiome and immune system are tightly connected, and their interaction is homeostatic. The microorganisms in the gut are a major source of antigen variation to which the immune system should control its response (8). The myelopoiesis and the epigenome of immune cells are affected by signals from the gut microbiome. Dysbiosis in the gut microbiome is associated with systemic effects on the immune system and can contribute to the development of autoimmune diseases (8,9).
Some reports have also suggested an association between the gut microbiome and the development of vascular complications (10). Alterations in gut microbial composition result in changes in the metabolites produced in the gut that are suggested to be partly responsible for the induction of chronic inflammation and to thereby promote cardiovascular disease (CVD) (10,11).
Only a few studies have investigated the role of the gut microbiome in long-standing type 1 diabetes. One study reported differences in the gut microbiome of 53 participants with type 1 diabetes with an average disease duration of 9 years compared with healthy control subjects (12). However, studies focusing on differences in the gut microbial composition in long-standing type 1 diabetes in larger cohorts and on the potential associations with diabetes-associated complications are lacking. Here, we characterized the gut microbiome using whole-genome shotgun sequencing in a cohort of patients with long-standing type 1 diabetes with a large range of glucose control, diabetes duration, and complications. We compared these results to a large, matched, control group to identify whether the microbiome composition is associated with type 1 diabetes–related traits.
Research Design and Methods
Cohorts and Metadata Collection
The study included two different cohorts: 240 participants with type 1 diabetes recruited in the Radboud University Medical Center, the Netherlands, and 2,937 age-, sex-, and BMI-matched healthy participants from the Lifelines Dutch Microbiome Project cohort.
Participants with type 1 diabetes were included from the outpatient diabetes clinic of the Radboud University Medical Center, the Netherlands. All participants were diagnosed with type 1 diabetes based on accepted clinical criteria, and data include age of onset, severity of symptoms, glucose variability, and absolute insulin dependency, sometimes in combination with documented positivity for autoantibodies (mainly anti-GAD). Exclusion criteria were age <18 years, pregnancy, and the use of antibiotics in the 4 weeks before inclusion. Tests were rescheduled in case of fever in the week before inclusion. This cohort is part of the Human Functional Genomics Project (13).
Besides general characteristics, in-depth information related to type 1 diabetes was collected, including HbA1c (mean HbA1c was used, based on 35 measurements on average), disease duration, and use of medication. In addition, we determined the presence of diabetes-associated complications, including retinopathy, neuropathy, nephropathy, myocardial infarction, stroke, and peripheral arterial disease (14,15). One participant had celiac disease and was on a long-term gluten-free diet. We did not observe any significant differences in α-diversity, β-diversity, or relative abundances of individual taxa or pathways between the participant and the remainder of the cohort, and he was therefore not excluded. As metformin is known to change the composition of the gut microbiome, we removed 2 metformin users from the cohort, resulting in 238 participants with type 1 diabetes being included in our analysis.
Lifelines is a multidisciplinary prospective population-based cohort study with a three-generation design that is examining the health and health-related behaviors of 167,729 people living in the north of the Netherlands (16). During the first Lifelines data collection visit, all participants were invited to participate in a parallel project—The Dutch Microbiome Project (DMP)—on a voluntary basis. The DMP’s goal is to evaluate the impact of different exposures and lifestyles on gut microbiota composition. A total of 10,000 participants were recruited in the DMP, and 8,208 were retained for downstream analysis after stringent quality control. DMP participants who did not self-report any serious or chronic diseases (any cancer, any form of diabetes, celiac disease, inflammatory bowel disease [Crohn disease or ulcerative colitis] and irritable bowel syndrome [diagnosed based on Rome III questionnaires]) were age-, sex-, and BMI-matched to the participants with type 1 diabetes using the MatchIt version (v)3.0.2 package for R (17). The ratio of control subjects to case subjects was selected as the maximal number of control subjects for which no significant difference in age, sex, or BMI could be identified between case subjects and control subjects using a Mann-Whitney U test, which resulted in 2,937 control subjects.
Determination of Gut Microbiome
Fecal samples from the two cohorts were collected and processed using a uniform workflow, as described previously (6). Briefly, fecal samples were collected by participants, frozen within 15 min after production, transported on dry ice, and stored at −80°C. Microbial DNA extraction was performed using QIAamp Fast DNA Stool Mini Kit (Qiagen, Hilden, Germany) according to the manufacturer’s instructions, and DNA library preparation and sequencing were performed at Novogene, Beijing, China. Library preparation of samples with a total DNA yield of <200 ng (as determined by Qubit 4 Fluorometer) was performed using NEBNext Ultra DNA Library Prep Kit for Illumina, while libraries for other samples were prepared using NEBNext Ultra II DNA Library Prep Kit for Illumina. Sequencing was performed using the Illumina HiSEq 2000 platform to generate ∼8 gigabase (Gb) of 150 base pair paired-end reads per sample (mean 7.9 Gb, SD 1.2 Gb).
Metagenomic sequencing data were analyzed using standardized methods we have used previously (6). In brief, quality control was performed by removing low-quality reads (Phred quality ≤30), adapters, and human DNA (reference Genome Reference Consortium Human build [GRCh] 37/hg19) using KneadData tools v0.5.1. Relative abundances of microbial taxa were determined using MetaPhlAn2 v2.7.2 (18), and the functional potential of the microbiome was profiled using the HUMAnN2 pipeline v0.11.1 (19) integrated with the UniRef90 v0.1.1 database of proteins (20), the ChocoPhlAn pangenome database, and the DIAMOND alignment tool v0.8.22 (21). In this study, we focused on bacterial and archaea taxa and pathways present in at least 1% of samples with a relative abundance >0.001%. We removed redundant pathways by discarding one pathway among highly correlated (Spearman R2 > 0.95) pathways.
Case-Control Analysis of Microbiome Signatures in Type 1 Diabetes
Microbiome α-diversity was calculated on the level of microbial species and genera using the Shannon diversity index, while β-diversity was calculated as the Bray-Curtis dissimilarity on the level of species. Principal coordinate analysis was used to project the Bray-Curtis dissimilarity matrix into a limited number of dimensions (n = 5) for exploratory data analysis, visualization of data, and clustering. Comparisons between centroids of groups (participants with type 1 diabetes vs. healthy control subjects) were performed on principal coordinates using the Mann-Whitney U test.
To quantify the proportion of microbiome variance explained by the disease and participant phenotypes, including age, sex, BMI, and medication use, we performed permutational multivariate ANOVA using distance matrices (ADONIS) analysis on the Bray-Curtis dissimilarity matrix using 10,000 permutations.
To determine microbiome signatures of type 1 diabetes, we compared microbiome diversity and the relative abundances of bacterial taxa and pathways between participants with type 1 diabetes and matched healthy control subjects by constructing Gaussian family multivariate generalized linear models, with age, sex, BMI, smoking, sequencing read depth, and proton pump inhibitor (PPI) use included as covariates to correct for the potential confounding effects of these variables. Relative abundances of microbiome taxa and pathways were normalized by centered log-ratio transformation. As none of the participants reported use of antibiotics up to 3 months before the fecal sampling, antibiotic use was not considered in the analysis. Results with a false discovery rate (FDR) of <0.05 were deemed statistically significant.
Analysis of Microbiome Signatures of Type 1 Diabetes Phenotypes
Analysis of the microbiome signatures of phenotypes related to type 1 diabetes was performed within the type 1 diabetes cohort using methods analogous to the case–control analysis described above in Case-Control Analysis of Microbiome Signatures in Type 1 Diabetes. We considered the following phenotypes in the study: HbA1c (mmol/mol), insulin sensitivity (derived from insulin need in insulin units/kg/day), duration of diabetes (years), glycemic load (defined as HbA1c ∗ disease duration), microvascular complications (including nephropathy, neuropathy, and retinopathy), and macrovascular complications (including peripheral arterial disease and myocardial infarction). Of these, HbA1c and insulin sensitivity were found to have highly skewed distributions and were log-transformed before analysis to normalize the distribution. All analysis included corrections for age, sex, BMI, read depth, smoking, and PPI use.
Statistical Tools and Software
Calculations were performed using R 3.6.0 for Windows. Generalized multivariate models were implemented using the MaAsLin2 (v2.1.0) package for R and the glm function from the R base package. ADONIS analysis and calculations of α- and β-diversity were implemented using the vegan (v2.5.6) package for R. The Benjamini-Hochberg correction was used to control for multiple testing.
Results
Baseline Characteristics
Characteristics of both cohorts are summarized in Supplementary Table 1. People with type 1 diabetes and healthy control subjects were well matched with respect to age, sex, and BMI. Participants with type 1 diabetes had an average age of 52 ± 16 years, an average HbA1c of 8.0% (63.6 mmol/mol), and an average disease duration of 28 ± 15 years. There were 40 participants (16.7%) with a history of macrovascular complications (myocardial infraction, stroke, or peripheral arterial disease), and 161 participants (67.3%) affected by one or more microvascular complications (retinopathy, neuropathy, or nephropathy).
Patients With Type 1 Diabetes Have Increased α-Diversity of Microbial Pathways
Data analysis revealed no differences in α-diversity of microbial taxa between participants with type 1 diabetes and healthy control subjects (Fig. 1). There was a small but significant (FDR = 0.001) increase in α-diversity of microbial biochemical pathways in participants with type 1 diabetes (Fig. 1G).
Gut Microbiome Composition of Patients With Type 1 Diabetes Is Significantly Different
Next, we compared the microbial species and pathways between participants with type 1 diabetes and healthy control subjects. Significant differences in various species and pathways were observed between both groups. A total of 37 bacterial species were significantly enriched in participants with type 1 diabetes, whereas 43 bacterial taxa were significantly depleted (Supplementary Table 2). Depleted species included Alistipes putredinis (FDR = 1.1 ∗ 10−11), Prevotella copri (FDR = 7.7 ∗ 10−10), and Bifidobacterium longum (FDR = 3.5 ∗ 10–5) (Fig. 2A–C). Enriched species in participants with type 1 diabetes included multiple species from families Ruminococcaceae (FDR = 2.5 ∗ 10−9), Clostridiaceae (FDR = 8.3 ∗ 10−8) (Fig. 2D), Clostridiales (FDR = 3.2 ∗ 10−5), and genus Oscillibacter (FDR = 0.01). Interestingly, enriched species consisted mainly of opportunistic pathogens (Clostridiales, Oscillibacter), whereas commensal species were significantly depleted in participants with type 1 diabetes (Supplementary Table 2). Regarding microbial pathways, 28 pathways were significantly depleted in participants with type 1 diabetes, whereas 11 were enriched (Supplementary Table 3). Depleted pathways included the pyrimidine deoxyribonucleotides de novo biosynthesis III pathway (PWY.6545) (FDR = 2.3 ∗ 10−8) and the adenosine deoxyribonucleotides de novo biosynthesis II pathway (PWY.7220) (FDR = 1.0 ∗ 10–4) (Fig. 2E and G). Enriched pathways included the l-isoleucine biosynthesis II pathway (PWY.5101) (FDR = 1.0 ∗ 10–4) and the stachyose degradation pathway (PWY.6527) (FDR = 4.2 ∗ 10−4).
Diabetes-Related Characteristics Explain a Significant Amount of Gut Microbiome Variance
To define the impact of diabetes-associated characteristics on the gut microbiome composition, microbiome variance per phenotype was determined. This approach was used to determine what factors would explain the variance in the gut microbiome between participants with type 1 diabetes and healthy control subjects and within participants with type 1 diabetes. In addition to looking at the individual species level, we also looked at enriched microbial pathways. When participants with type 1 diabetes were compared with healthy control subjects, beyond having type 1 diabetes, HbA1c explained a significant part of the variance in species (explained variance R2 > 0.004, FDR = 0.0017) and pathways (explained variance R2 > 0.003, FDR = 0.0014) (Fig. 3A and B). The diabetes-associated characteristics that impacted the microbiome variance included duration of diabetes, age of onset, HbA1c levels/glycemic load, and medication use. Age explained the largest amount of variation in microbiome composition (explained variance R2 = 0.013, FDR = 0.0066) (Fig. 3C), followed by the use of PPIs or platelet aggregation inhibitors (explained variance R2 = 0.011, FDR = 0.0092 and R2 = 0.011, FDR = 0.0092, respectively) and HbA1c (explained variance R2 = 0.010, FDR = 0.0092) (Fig. 3C). Within the cohort with type 1 diabetes, none of the diabetes-related clinical characteristics impacted the microbial pathway variation significantly (Fig. 3D).
Diabetes-Related Characteristics Show Significant Associations With Microbial Species and Pathways
Since we observed a significant effect of diabetes-related characteristics on the microbiome variance, we set out to determine species-specific associations. Glycemic control, as reflected by HbA1c, showed a significant association with several microbial species, with the strongest positive association found with Dorea formicigenerans (FDR = 0.028) and the strongest negative associations found with genus Fecalibacterium (P = 0.0016), which was dominated by Fecalibacterium prausnitzii (P = 0.0016). Disease duration and HbA1c showed positive associations with bacteria from family Ruminococcaceae (P < 0.05), and disease duration showed a negative association with Roseburia intestinalis (P = 0.01). HbA1c also showed significant negative and positive associations with several microbial pathways. The highest negative FDRs were found for the galacturonate degradation pathway I (FDR = 0.008) and the super pathway of l-serine and glycine biosynthesis I (FDR = 0.008). Again, disease duration and glycemic load displayed similar associations, showing significant positive associations with both the flavin biosynthesis I bacteria and plants pathway and the super pathway of pyrimidine ribonucleotides de novo biosynthesis pathway.
Microbial Species and Pathways Are Significantly Associated With Diabetes-Related Complications
The presence of micro- and macrovascular complications explained a significant amount of the variance in gut microbiome composition on species level (R2 > 0.0075, FDR < 0.05) (Fig. 3). Overall, the different microvascular complications showed similar associations, whereas different macrovascular complications displayed different associations.
Diabetic nephropathy, retinopathy, and neuropathy were strongly associated with several microbial species (Fig. 4C), with the most significant associations found for diabetic nephropathy, which was positively associated with seven Clostridium species (FDR < 0.05), genus Bacteroides (FDR = 0.042), and family Bacteroidaceae (FDR = 0.03). Retinopathy showed a positive association with 15 species, mainly Clostridium species, and a negative association with 7 species, including R. intestinalis (P = 0.015). Diabetic neuropathy was positively associated with 14 species, of which the strongest association was with Anaerostipes (P = 0.0007), but there was also an association with two Ruminococcaceae species (P < 0.05). Diabetic nephropathy was also associated with a variety of microbial pathways, whereas fewer associations were found for other complications. The strongest positive associations were found for the folate transformation II pathway (PWY.3841) (FDR = 0.001), the formyl tetrahydrofolate biosynthesis pathway (FDR = 0.0022), and the super pathway of pyridoxal 5 phosphate biosynthesis and salvage (PWY.0.845; FDR = 0.0039).
Presence of peripheral arterial disease and myocardial infarction were strongly associated with the same microbial species, but peripheral arterial disease was associated with more species (Fig. 4C). Although not FDR significant, these complications were positively associated with Slackia (P < 0.05), order Clostridiales, and five Clostridium species. One negative association was found for both macrovascular peripheral arterial disease and myocardial infarction, which was with Bifidobacterium adolescentis (P < 0.05). Macrovascular complications were positively associated with seven pathways, including the pyridoxal 5 phosphate biosynthesis I pathway (Fig. 4D). Macrovascular complications were negatively associated with three pathways, including the super pathway of pyridoxal 5 phosphate biosynthesis and salvage (PWY.0.845; P = 0.025) (Fig. 4D).
Conclusions
The current study identified clear differences in gut microbiome composition between participants with long-standing type 1 diabetes and healthy control subjects, with a significant part of the variation seemingly related to glycemic control. Gut microbiome composition was also associated with the presence of micro- and macrovascular complications, particularly diabetic nephropathy, independent of age, sex, BMI, PPI use, and smoking. Overall, our study provides evidence for a difference in gut microbiome composition in participants with long-standing type 1 diabetes compared with healthy control subjects that associates with glycemic control and disease-related complications. Given the cross-sectional design, this study cannot determine cause-and-effect relationships.
Previous research on the role of the gut microbiome in participants with type 1 diabetes primarily focused on the pathogenesis of the disease and found reduced diversity in gut microbiome composition in children with type 1 diabetes and children at risk for developing type 1 diabetes (22). These data are suggestive for the involvement of decreased microbial diversity in the autoimmune process that drives type 1 diabetes development. Interestingly, we did not observe significant differences in microbial diversity between participants with long-standing type 1 diabetes and healthy control subjects, a finding consistent with two previous studies in participants with long-standing type 1 diabetes (12,23). These results suggest that a recovery in diversity occurs over time.
We did identify differences in gut microbial species between participants with type 1 diabetes and healthy control subjects. Previous studies using 16S sequencing revealed conflicting results on this topic, finding both enrichment and depletion of Bifidobacteria in long-standing type 1 diabetes (12,24). Alongside a smaller cohort size, these differences between previous studies could be explained by the average disease duration of 10 years and the younger age of the participants (12,24). In our study, where participants had an average diabetes duration of 28 years, we observed a reduction in Bifidobacteria. Bifidobacteria are producers of the short-chain fatty acid butyrate, which has anti-inflammatory effects through suppression of nuclear factor-κB signaling, induces T regulatory cells in the gut mucosa, and is important in tight-junction expression (25–27). The reduction in Bifidobacteria may be related to the more proinflammatory status in type 1 diabetes driving complications. Alongside a reduction in Bifidobacteria, we also observed a reduction in the commensal species D. longicatena. This species is highly prevalent in healthy individuals and has therefore been associated with a healthy gut microbiome (28). Interestingly, a lower abundance of D. longicatena in individuals with type 1 diabetes was previously observed to result in an increased risk of developing heart failure (29). Besides a lower abundance of several species, we found increased levels of Clostridiales and Oscillibacter in type 1 diabetes. These species are opportunistic pathogens that have previously been linked to type 2 diabetes and obesity in humans and mice (30–32). Overall, our data suggest that the differences in the gut microbiome could be caused by a decrease in species linked to a healthy gut microbiome and an increase in species associated with disease.
Using the diabetes-associated characteristics collected for our study cohort, we tried to identify the driving force behind the changes in the microbiome composition in participants with type 1 diabetes. Hyperglycemia is an important risk factor in the development of several diabetes-related complications (33), and poor glycemic control is associated with systemic low-grade inflammation (33,34). HbA1c explained a significant amount of the variation found in the gut microbiome of participants with type 1 diabetes. At the species level, we found a significant positive association between HbA1c and D. formicigenerans. This species has previously been linked to the presences of glycoside hydrolases involved in mucin degradation and can thereby affect gut permeability (35). The positive association between HbA1c and D. formicigenerans may impact gut permeability in participants with poor glycemic control, and increased gut permeability can contribute to a higher systemic inflammatory status (11).
Disease duration and poor glycemic control were also positively correlated with Ruminococcaceae, which were previously found to be correlated with clinical characteristics such as HbA1c, blood pressure, and lysophosphatidylcholine (LPC)/phosphatidylcholine (PC) metabolites in patients with type 2 diabetes (36). The association with LPC/PC metabolites was found to be strongly age dependent (36). Since we found an association with disease duration, we suggest that Ruminococcaceae, and thereby LPC/PC metabolites, increase over time and might contribute to the development of CVD.
Our results suggest a link between the level of glycemic control and gut microbiome composition but cannot provide insight into causality. We hypothesize that poor glycemic control could result in an altered gut microbiome composition that may subsequently contribute to the development of systemic inflammation and an increased risk of complications. Conversely, the altered gut microbiome composition may affect glycemic control and thereby drive the risk of developing complications. Previous research in people with prediabetes reported clear associations between the composition of the microbiome and postprandial glucose excursions in response to carbohydrate ingestion (37). Further research would be needed to determine the causal relation between gut microbiome composition and diabetes-related characteristics.
The presence of diabetes results in an increased risk of developing micro- and macrovascular complications. CVD, which is linked to chronic inflammation (38), is the leading cause of death in participants with type 1 diabetes, and recent studies suggest that alterations in the gut microbiome could contribute to chronic inflammation and CVD risk (10). This is in line with our results, also revealing an association between complications and gut microbiome composition. The strongest associations were found between microvascular complications and the gut microbiome, particularly with nephropathy, which were mainly with family Clostridioforme (C. bolteae and C. clostridioforme) and Bacteroides, known opportunistic pathogens (39,40). There were fewer associations between other microvascular complications and microbial species. Nephropathy is considered as one of the most severe complications associated with diabetes. Most participants with type 1 diabetes had several complications, so although the number of participants with nephropathy was not high, many had additional complications. This may partially explain the more pronounced association. Furthermore, previous literature has shown an important role for the gut-kidney axis in the development of kidney disease (41), which may explain the stronger associations between nephropathy and the gut microbiome in our study.
We observed that the presence of microvascular complications was associated with increased levels in Clostridia, Dialister, and Anaerotruncus and decreases in butyrate producers from genera Roseburia and Alistipes. Butyrate producers are crucial for colon health due to their role in mucin production and their anti-inflammatory properties (42). We found similar associations between Clostridioforme (C. boltee and C. clostridioforme) and Bacteroides in participants with macrovascular complications. Furthermore, we found an association between micro- and macrovascular complications and Ruminococcus. Higher levels of Ruminococcus have previously been found in diabetes-prone rats. A common metabolite of Ruminococcus is butyrate, which is important for colon health (42). In a healthy gut, besides mucin-producing species, mucin-degrading species, such as Prevotella, are also important for mucin production (43). We observed that Prevotella were reduced in participants with type 1 diabetes, which could contribute to a reduction in mucin production, thereby increasing gut permeability.
Microbial species also produce metabolites, with one interesting metabolite being trimethylamine-N-oxide (TMAO). TMAO’s precursor (trimethylamine) is mainly produced by Firmicutes, and trimethylamine is oxidized to TMAO by the hepatic enzyme flavin monooxygenase 3 (FMO3). Plasma levels of TMAO have been shown to be predictive for the development of CVD and atherosclerosis (44). P. copri, which was increased in participants with type 1 diabetes in our cohort, has been linked to increased TMAO production (45). Furthermore, several Ruminococcaceae, part of the phylum Firmicutes, show a positive association with microvascular neuropathy and macrovascular peripheral arterial disease in our data. Overall, these data support the hypothesis that alterations in the gut microbiome can increase the risk of developing micro- and macrovascular complications via microbially produced metabolites such as butyrate and TMAO.
Our study has limitations. The cross-sectional design makes it impossible to determine causality. In addition, while our cohort consisted of a reasonable number of participants, it is relatively small for this type of study. The correlation structure in our data (HbA1c is correlated with complications and complications are correlated with other complications) and the limited number of cases of individual complications limited our power to detect specific complication-microbiome associations. Yet, our cohort included carefully selected people with type 1 diabetes with a wide range of disease duration, glycemic control (HbA1c), and presence of complications, which renders our findings generalizable to the majority of the population with established type 1 diabetes.
The study also has strengths, including the detailed microbiome analyses, the large matched control population, and the large amount of clinical information obtained for the cohort with type 1 diabetes.
In conclusion, the composition of the gut microbiome in participants with type 1 diabetes differs significantly from the microbiome of healthy control subjects, while diversity is not affected. Changes in the gut microbiome are associated with glycemic control and disease-related complications. Altogether, these results suggest that the gut microbiome is not only important in the context of type 1 diabetes development but may also be involved in the development of diabetes-associated complications.
This article contains supplementary material online at https://doi.org/10.2337/figshare.20032256.
J.I.P.v.H., and R.G. share first authorship. L.A.B.J., and C.J.T. share senior authorship.
This article is featured in a podcast available at diabetesjournals.org/journals/pages/diabetes-core-update-podcasts.
Article Information
Acknowledgments. The authors thank all participants for their contribution. The authors thank Anouk Janssen, Department of Internal Medicine, Radboudumc, and Martin Jaeger, Department of Internal Medicine, Radboudumc, for the sample collection and Anneke Hijmans and Kiki Schraa, Department of Internal Medicine, Radboudumc, for processing the samples. The authors would like to thank the Center for Information Technology of the University of Groningen (Rijksuniversiteit Groningen [RUG]) for their support and for providing access to the Peregrine high-performance computing cluster, the Genomics Coordination Center (University Medical Center Groningen and RUG) for their support and for providing access to Calculon and Boxy high-performance computing clusters, and the MOLGENIS (Molecular Genetics Information System) team for data management and analysis support. The authors also thank K. McIntyre, Department of Genetics, University of Groningen, for English and content editing.
Funding. Sequencing of the cohort was funded by a grant from CardioVasculair Onderzoek Nederland (CVON 2012-03) to J.F. and A.Z.; J.I.P.v.H., R.S., L.A.B.J., C.J.T., R.G., and R.K.W. are supported by the collaborative TIMID project (LSHM18057-SGF) financed by the public-private-partnership allowance made available by Top Sector Life Sciences & Health to Samenwerkende Gezondheidsfondsen (SGF) to stimulate public–private partnerships and cofinanced by health foundations that are part of the SGF. R.K.W. is supported by the Seerave Foundation and the Dutch Digestive Foundation (16–14). A.Z. is supported by European Research Council (ERC) Starting Grant 715772, Netherlands Organization for Scientific Research (NWO) VIDI grant 016.178.056, CVON grant 2018–27, and NWO Gravitation grant ExposomeNL 024.004.017. J.F. is supported by the Dutch Heart Foundation IN-CONTROL (CVON2018–27), the ERC Consolidator grant (grant agreement no. 101001678), NWO-VICI grant VI.C.202.022, and the Netherlands Organ-on-Chip Initiative, an NWO Gravitation project (024.003.001) funded by the Ministry of Education, Culture and Science of the government of the Netherlands.
Duality of Interest. No potential conflicts of interest relevant to this article were reported.
Author Contributions. J.I.P.v.H., and R.G. wrote the first draft of the manuscript. R.G. analyzed the data. R.S. reviewed the first draft of the manuscript. R.S., R.K.W., H.J.M.H., L.A.B.J., and C.J.T. conceived and planned the study. J.F., A.Z., and H.J.M.H. contributed to the bioinformatical framework and the DArm Gezondheid (DAG) 3 cohort. All authors provided critical feedback and helped to shape the research and the manuscript. L.A.B.J. and C.J.T. are the guarantors of this work and, as such, had full access to all the data in the study and take responsibility for the integrity of the data and the accuracy of the data analysis.