Distinct fecal and oral microbiota composition in human type 1 diabetes, an observational study

PLOS ONE, Dec 2017

Objective Environmental factors driving the development of type 1 diabetes (T1D) are still largely unknown. Both animal and human studies have shown an association between altered fecal microbiota composition, impaired production of short-chain fatty acids (SCFA) and T1D onset. However, observational evidence on SCFA and fecal and oral microbiota in adults with longstanding T1D vs healthy controls (HC) is lacking. Research design and methods We included 53 T1D patients without complications or medication and 50 HC matched for age, sex and BMI. Oral and fecal microbiota, fecal and plasma SCFA levels, markers of intestinal inflammation (fecal IgA and calprotectin) and markers of low-grade systemic inflammation were measured. Results Oral microbiota were markedly different in T1D (eg abundance of Streptococci) compared to HC. Fecal analysis showed decreased butyrate producing species in T1D and less butyryl-CoA transferase genes. Also, plasma levels of acetate and propionate were lower in T1D, with similar fecal SCFA. Finally, fecal strains Christensenella and Subdoligranulum correlated with glycemic control, inflammatory parameters and SCFA. Conclusions We conclude that T1D patients harbor a different amount of intestinal SCFA (butyrate) producers and different plasma acetate and propionate levels. Future research should disentangle cause and effect and whether supplementation of SCFA-producing bacteria or SCFA alone can have disease-modifying effects in T1D.

A PDF file should load here. If you do not see its contents the file may be temporarily unavailable at the journal website or you do not have a PDF plug-in installed and enabled in your browser.

Alternatively, you can download the file locally and open with any standalone PDF reader:

http://journals.plos.org/plosone/article/file?id=10.1371/journal.pone.0188475&type=printable

Distinct fecal and oral microbiota composition in human type 1 diabetes, an observational study

December Distinct fecal and oral microbiota composition in human type 1 diabetes, an observational study Pieter F. de Groot 0 1 2 Clara Belzer 0 2 OÈ mruÈ m Aydin 0 1 2 Evgeni Levin 0 1 2 Johannes H. Levels 0 1 2 Steven Aalvink 0 2 Fransje Boot 0 1 2 Frits Holleman 0 1 2 DanieÈ l H. van Raalte 0 2 Torsten P. Scheithauer 0 1 2 Suat Simsek 0 2 Frank G. Schaap 0 2 Steven W. M. Olde Damink 0 2 Bart O. Roep 0 2 Joost B. Hoekstra 0 1 2 Willem M. de Vos 0 2 Max Nieuwdorp 0 1 2 0 Competing interests: WM de Vos is on the Scientific Advisory Board of Caelus Health , and M 1 Department of Internal and Vascular Medicine, Academic Medical Center±University of Amsterdam, Amsterdam, the Netherlands, 2 Laboratory of Microbiology, Wageningen University, Wageningen, the Netherlands, 3 Department of Internal medicine, VU University Medical Center , Amsterdam , The Netherlands, 4 ICAR , VU University Medical Center , Amsterdam , The Netherlands , 5 Department of Internal Medicine, Medisch Centrum Alkmaar, Alkmaar, the Netherlands, 6 Department of Surgery, Maastricht University , Maastricht , The Netherlands , 7 NUTRIM School of Nutrition and Translational Research in Metabolism, Maastricht, the Netherlands, 8 Department of General, Visceral and Transplantation Surgery, RWTH University Hospital Aachen , Aachen, Germany , 9 Department of Immunohaematology & Blood Transfusion, Leiden University Medical Center, Leiden, the Netherlands, 10 Beckman Research Institute, DMRI, City of Hope , Duarte, CA , United States of America , 11 RPU Immunobiology , University of Helsinki , Helsinki, Finland, 12 Wallenberg Laboratory , Department of Molecular and Clinical Medicine, Sahlgrenska Academy, University of Gothenburg , Gothenburg , Sweden 2 Editor: Nalini Kumar Vudattu, Pfizer Global Research and Development , UNITED STATES - Data Availability Statement: All relevant data are within the paper and its Supporting Information files. Funding: M Nieuwdorp is supported by a ZONMWVIDI grant 2013 (016.146.327) of the Netherlands Organization for Scientific Research (NWO, www. nwo.nl) and CVON Young Talent grant 2012. WM de Vos is supported by the NWO Gravitation Grant (SIAM 024.002.002) and the Spinoza Award. Objective Environmental factors driving the development of type 1 diabetes (T1D) are still largely unknown. Both animal and human studies have shown an association between altered fecal microbiota composition, impaired production of short-chain fatty acids (SCFA) and T1D onset. However, observational evidence on SCFA and fecal and oral microbiota in adults with longstanding T1D vs healthy controls (HC) is lacking. Research design and methods We included 53 T1D patients without complications or medication and 50 HC matched for age, sex and BMI. Oral and fecal microbiota, fecal and plasma SCFA levels, markers of intestinal inflammation (fecal IgA and calprotectin) and markers of low-grade systemic inflammation were measured. Results Oral microbiota were markedly different in T1D (eg abundance of Streptococci) compared to HC. Fecal analysis showed decreased butyrate producing species in T1D and less butyryl-CoA transferase genes. Also, plasma levels of acetate and propionate were lower in T1D, with similar fecal SCFA. Finally, fecal strains Christensenella and Subdoligranulum correlated with glycemic control, inflammatory parameters and SCFA. Nieuwdorp is on the Scientific Advisory Board of Caelus Health. The other authors have nothing to disclose. Conclusions We conclude that T1D patients harbor a different amount of intestinal SCFA (butyrate) producers and different plasma acetate and propionate levels. Future research should disentangle cause and effect and whether supplementation of SCFA-producing bacteria or SCFA alone can have disease-modifying effects in T1D. Introduction Type 1 Diabetes (T1D) is the second most frequent autoimmune disease in childhood and its incidence has tripled in the last thirty years. Notwithstanding decades of intensive research in animals, the environmental factors driving T1D are still unknown and therapeutic strategies have invariably failed to halt disease progression. As the increased T1D incidence is primarily observed in subjects who are not genetically predisposed, environmental factors including altered diet[ 1 ], infections with concomitant antibiotic use[ 2 ] as well as mode of birth[ 3 ], have been suggested to play a role. As autoimmune beta cell inflammation is one of the hallmarks of T1D, this insulitis may originate from an innate immune response to intestinal pathogens[4]. Accordingly, an increased amount of pathogenic bacterial species has been observed in fecal samples of T1D patients around time of diagnosis[ 5 ] and an altered composition of the fecal microbiota composition was observed in adolescent T1D patients[ 5 ]. Interestingly, this altered fecal microbiota is already present before the clinical onset of T1D[ 6 ] and is related to islet autoantibodies[ 6 ]. Moreover, a decrease in fecal short-chain fatty acid (SCFA) producing bacteria was observed in small studies of T1D subjects[ 7 ], while dietary acetate and butyrate supplementation elicits favorable immunological effects and protection from T1D in NOD mice [ 8 ]. Most human studies have been performed in children and adolescents with T1D[5,7], whereas patients with longstanding T1D have only been studied in a small group[ 9 ] without taking several confounding factors including gender and diet into account. More recently, oral microbiota have been implicated in development of type 2 diabetes[ 10 ] and liver disease[ 11 ], but these have never been studied in T1D. We therefore set out to study both fecal and oral microbiota composition, as well as SCFA metabolism in a large matched case-control study of T1D patients and control subjects. Methods Recruitment and study visit Outpatient clinics of six medical centers in the Amsterdam region in the Netherlands were screened for eligible T1D patients. From a total amount of 854 T1D subjects, 82 eligible subjects were selected by their treating physician based on our inclusion criteria and approached for participation. The large majority of noneligible patients were excluded due to cormorbidity affecting glucose control or gut microbiota (eg recent antibiotic use or use of proton pump inhibitors) or because end organ damage was present. In a few cases (N = 5) other reasons were given (e.g. expected nonadherence, mental performance, language barrier). Finally, from 82 eligible patients, 53 subjects agreed to participate. In total, 83 control subjects were recruited by advertisement of which 50 were eligible for participation. Subjects were matched for age, sex and BMI. Written informed consent was obtained from all subjects. The study was conducted at the Academic Medical Center (Amsterdam) in the period between November 2013 2 / 14 and April 2014, in accordance with the Declaration of Helsinki (updated version 2013). The study was approved by the ethics review board of the Academic Medical Center. Inclusion criteria for T1D patients were of Western European descent, age 18±65 years, normal BMI (18.5±25 kg/m2), and a Western dietary pattern. Exclusion criteria were known determinants of altered microbiota composition including medication use including statins and proton pump inhibitors, use of antibiotics three months prior to inclusion, use of probiotic-containing food, unusual dietary habits (e.g. vegan diet) and medical conditions believed to affect glucose metabolism or gut microbiota (e.g. cholecystectomy, celiac disease and irritable bowel syndrome)[ 12 ]. As we aimed to investigate uncomplicated longstanding T1D patients, we excluded poorly regulated subjects (HbA1c > 10% or 86 mmol/mol) as well as subjects with microvascular complications of T1D (nephropathy, neuropathy or retinopathy). Both T1D and controls were allowed to continue their diet and were asked to fill out an online nutritional diary (www.voedingscentrum.nl) for the duration of one week before the study visit to monitor caloric intake including the amount of dietary carbohydrates, fat, proteins and fibers. Anthropometric and physiological measurements including blood pressure, length, weight, hip and waist circumference) were recorded. After an overnight fast, blood was drawn and fresh morning feces were collected. Participants were asked to refrain from tooth brushing in the twelve hours before collection of oral swabs, that were taken taken between the upper lip and the front row of teeth for oral microbiota analysis as previously published[ 10 ]. All samples were stored at -80ÊC until analyzed. Plasma samples Total cholesterol, low density lipoprotein cholesterol (LDLc), high density lipoprotein choles terol (HDLc) and triglycerides (TG) were measured by using commercially available enzymatic assays (Randox, USA and Daiichi, Japan). All analyses were performed using a Cobas Mira autoanalyzer (Horiba, France). C-Reactive Protein (CRP, Roche Diagnostics), LPS-bind ing protein (Human LBP ELISA kit, Hycult biotech, catalog no. HK315), fecal calprotectin (Buhlmann), and fecal IgA (human IgA ELISA Kit, eBioscience) were determined. Plasma SCFA were determined using liquid chromatography-mass spectrometry according to van Eijk et al[13], with minor modifications. Fecal samples Fecal SCFA were measured using high-performance liquid chromatography with UV detection as previously published by De Baere et al.[ 14 ]. DNA was extracted from non-thawed fecal samples using a bead-beating protocol and analysis of microbiota diversity and composition was performed by Illumina Miseq sequencing (Illumina, San Diego, CA, USA) of 16S rRNA genes from extracted DNA with primers 27F-DegS and 338R[ 15 ]. The average number of reads was 12158 with a lower limit of 1187. Sequencing data were analyzed using NG-tax[ 16 ] and for follow up analyses the QIIME software package (available at http://qiime.sourceforge. net/), the Canoco 5 software package (Biometris, Wageningen, the Netherlands) and R-studio were used. Several fecal and saliva samples were not collected, properly stored or lost during the DNA extraction or processing steps. Hence, fecal microbiota were assessed in 45 T1D vs 35 controls. Oral microbiota were assessed in 51 T1D vs 42 controls. Fecal SCFA were assessed in 43 T1D vs 47 controls. All other parameters including plasma SCFA were assessed in all participants. Butyrate production capacity was assessed by qPCR targeting the Butyryl-coenzyme A (CoA)-CoA transferase (ButCoA) gene using 25 ng of template DNA and 0.5μM primer concentration in each reaction as previously described[ 17 ]. qPCR amplifications were performed 3 / 14 in triplicate using MX3005P Real-Time PCR System (Stratagene, La Jolla, CA, USA) in a volume of 25 μl. Each reaction was amplified by using 5 μl of HOT FIREPol EvaGreen qPCR Mix Plus, no ROX (Solis BioDyne, Tartu, Estonia). The number of copies of the butyryl-CoA transferase-acetate-CoA transferase gene in feces was measured and expressed as a percentage of the total amount of 16S rRNA gene copies, which allows to make an estimation of the abundance of all butyrate producing bacteria. Statistics To test significance between group differences, unpaired Student t-test or the Mann-Whitney test were used, depending on the distribution of the data. Accordingly, data are expressed as mean ± the standard deviation or as median with interquartile range. For correlation analyses, Spearman's Rank test was used (as all parameters shown had non-normal distribution). A p value < 0.05 was considered significant. Differences in fecal and oral microbiota composition between groups (beta diversity) were assessed visually using RDA-plots and per species with Wilcoxon's signed rank test. A false-discovery-rate corrected p-value (q-value <0.05) was considered statistically significant; `qvalues' package in R was used for testing. For testing alpha diversity (the species diversity within an individual or mean individual diversity within a group) we used Shannon's diversity index. Intestinal bacterial species discriminating between controls and the T1D group were selected using the elastic net algorithm[ 18,19 ]. A randomization test was conducted to evaluate the statistical validity of the results obtained via elastic net algorithm. We followed the procedure where the outcome variable (e.g T1D vs. control) was randomly reshuffled while the corresponding microbial profiles were kept intact. This was repeated up to 100 times and Receiver-Operating-Characteristics Area-Under-Curve (ROC AUC) scores were computed each time. The performance measure used for a binary classification task is an ROC AUC. The ROC can be understood as a plot of the probability of correctly classifying T1D vs. control subjects. Cross-validation within the dataset was accomplished by randomly hiding 20% of the subjects from the model and evaluating the prediction quality on that group. The ROC AUC score measures the predictive accuracy of the classification model with 0.5 AUC corresponding to a random result. A significance value of 0.05 was defined and the true AUC of the original dataset was compared with this value. All statistical analyses were performed using Numerical Python, R, and MATLAB. Results Baseline characteristics of both study groups are depicted in Table 1. T1D and matched controls did not differ in age, sex or BMI. Total caloric intake (proteins, fat or carbohydrate) was not significantly different between groups, in spite of a trend towards higher fiber intake observed in T1D subjects (Fig 1A). Although fecal SCFA levels were not different between both groups (Table 1), in fasting plasma of T1D we observed lower acetate and propionate levels (Fig 1B±1D). While plasma acetate and butyrate were significantly correlated in T1D (r = +0.357, p = 0.009), a negative correlation between plasma acetate and propionate was observed in the controls (r = -0.432, p = 0.002) (S1 Fig). Markers of chronic inflammation (e.g. LBP and CRP) as well as fecal calprotectin and fecal IgA were not significantly different between the groups (see Table 1). However, a correlation between fecal calprotectin and plasma LPS-binding protein (LBP) levels (r = +0.23, p = <0.05), and plasma LBP and CRP (r = 0.6, p<0.001) was noted in both groups. Furthermore, in the T1D group HbA1c was directly correlated to plasma LBP (r = +0.3, p<0.05), whereas fecal IgA correlated strongly with fecal butyrate in the T1D group (r = +0.65, p<0.001). Finally, although age range-dependent sex differences have been reported in relation to the incidence of type 1 diabetes[ 20,21 ], no difference was found in 4 / 14 Group characteristics expressed as mean ± SD or median [IQR] with p-value using independent T-test for parametric and Mann-Whitney test for nonparametric data. inflammatory parameters or short-chain fatty acids between males and females in either group. Fecal and oral microbiota composition A difference in fecal microbiota beta diversity was seen as depicted by RDA plot (Fig 2A), while there was no significant in alpha diversity (Shannon's diversity index) in T1D (N = 45) versus controls (N = 35). Regarding fecal microbiota, on the level of individual taxa only one genus and one family were significantly different (p<0.05) between the groups; Christensenella (p = 0.015) and Rhodospirillales (p = 0.009) whereas the adjusted q-values did not remain significant; When applying the elastic net algorithm to identify intestinal bacterial species most discriminating as a group between T1D and controls, we could identify several fecal and oral species with high sensitivity (AUC 0.88) that were predictive to belong to either the T1D or the HC group (Table 2 and Fig 3A and 3B). In feces, Christensenella associated with T1D and cor related negatively with fecal acetate (Spearman r = -0.416, p<0.05), whereas Subdoligranulum correlated significantly with plasma markers of endotoxemia (LBP, r = 0.519, p<0.001) and inflammation (CRP, r = 0.362, p<0.05, S2A Fig) in T1D. Also, butyryl-CoA:acetate-CoAtransferase gene ratio was significantly decreased in T1D fecal samples (3.1± 1.1% vs 5.8±0.6%, p = 0.03), which is in line with our findings of lower Roseburia, a known butyrate producer in T1D (Table 2 and Fig 3A and 3B). Finally, these key fecal microbiota species from the elastic 5 / 14 Fig 1. Plasma SCFA and dietary fiber intake. The figure shows that (A) despite higher fiber intake, we (B) found decreased plasma levels of acetate and (C) propionate and (D) similar levels of plasma butyrate in controls vs T1D subjects. net analyses were correlated with clinical and metabolic parameters, which is summarized in S2A and S2B Fig. With respect to oral microbiota, beta diversity also differed between the groups (Fig 2B). T1D subjects were characterized by a significantly higher abundance of taxa belonging to the phyla Actinobacteria and Firmicutes, including Streptococcus spp., Actinomyces spp. and Rothia spp (q values <0.05) (see Table 3 and Fig 4A and 4B). In contrast, Bacteroidetes and Proteobacteria phyla were significantly increased in controls including Pasteurellaceae (q = 0.032). In contrast with fecal microbiota, we found no significant correlations between oral microbiota and glycemic control, inflammatory parameters or short-chain fatty acids. Discussion Recent studies have linked altered fecal microbiota composition to the recent rise in many immunological disorders including T1D. Indeed, the intestinal microbiota composition seems to differ between new onset adolescent T1D subjects and controls[ 5 ]. In this regard, antibiotic administration has been associated with T1D incidence both in mice[ 22 ]and humans[ 23 ]. In patients with longstanding T1D, we now report that oral as well as fecal microbiota are altered compared to healthy matched controls. While previous studies reported a decrease in alpha diversity in young adolescent T1D[ 5 ], we only found significant changes in beta diversity in both fecal and oral microbiota in our well-controlled longstanding T1D patients, which is in line with a previous smaller study in longstanding T1D[ 9 ]. Moreover, on the genus level we found that fecal samples of T1D were enriched in Christensenella and Bifidobacteria (Table 2 and Fig 3A and 3B), which suggests that low fecal Bifidobacterium levels reported around the time of diagnosis[ 7 ] have recovered in longstanding T1D patients. Our finding of increased abundance of Christensenella, which has been previously linked to lower obesity risk,[ 24 ] is a novel finding in T1D. In line with lower plasma SCFA levels, decreased fecal levels of the butyryl-CoA:acetate-CoA transferase gene as well as SCFA producers like Roseburia, were seen in the T1D group. These findings are 6 / 14 Fig 2. Fecal microbiota composition. Differences in (A) fecal and (B) oral microbiota composition as depicted by a biplot of Redundancy Analysis (RDA axis 1 vs. axis 2) constrained by T1D or controls. in line with a recent study in non-obese diabetic (NOD) mice that showed an inverse correlation with key features of T1D and plasma blood SCFA concentrations, whereas diet containing 7 / 14 Mean relative abundance, percentage of subjects harboring species and ratio of most discriminating fecal microbiota from our predictive models analysis. Taxa are ordered by weight in the predictive model from highest positive to lowest negative. *Taxa shown are discriminatory between cases and controls as a group and not on the individual level. Nevertheless, calculated p and q values (using Mann-Whitney U-test and q-value package in R) are included for completeness. **ªT1Dº type 1 diabetes, ªHCº healthy controls, ªinfº infinite. Fig 3. Most discriminating fecal species. Relative abundances of most discriminating fecal microbiota resulting from our elastic net analysis for (A) T1D patients and (B) healthy controls (B). 8 / 14 Mean relative abundance, percentage of subjects harboring species and ratio of most discriminating oral microbiota from our predictive models analysis. Taxa are ordered by weight in the predictive model from highest positive to lowest negative. *Taxa shown are discriminatory between cases and controls as a group and not on the individual level. Nevertheless, calculated p and q values (using Mann-Whitney U-test and q-value package in R) are included for completeness. **ªT1Dº type 1 diabetes, ªHCº healthy controls, ªinfº infinite. SCFA boosted the number and function of regulatory T cells, enhanced gut integrity and decreased chronic low grade inflammation[ 8 ]. In this study various in vivo and in vitro experiments using transfer of T-cell subsets and labeling show that oral administration of acetate changes the composition of B cell subsets, thereby affecting the proliferation of T-cells in the spleen and pancreatic lymph nodes whereas infiltration of pancreatic islets was decreased, thereby protecting from the development of autoimmune diabetes[ 8 ]. In contrast, no protective effects of propionate against T1D development were observed[ 8 ]. Nevertheless literature on role of SCFAs in T1D or its effects on the adaptive immune system in general are scarce and warrants further study. Interestingly, we also found changes in bacteria in the proximal part of the intestine with Streptococcus in the oral cavity being positively associated with T1D, whereas fecal Streptococcus showed an inverse relationship (Table 2 and Fig 3A and 3B). This might be clinically relevant since the immune system in the small intestine of T1D subjects seems to reacts differently than that of healthy controls to Streptococcus which are abundant throughout the entire alimentary tract[ 15 ]. In line with our data, Subdoligranulum has been reported to be associated with poorer metabolic control and chronic inflammation[ 25 ]. Besides altered function of the adaptive (T-cell) immune system[ 26 ], T1D has been linked to an altered innate immune function[ 27 ] including changes in Toll-like receptor signaling[ 28 ] and increased bacterial translocation[ 29 ]. In line, a recent study showed that gut-derived bacterial lipopolysaccharide from Bacteroides may preclude proper immune system regulation[30]. Furthermore, our observed correlation between fecal IgA and fecal butyrate levels suggests that SCFA play a role in differential IgA targeting of specific intestinal bacteria in (de novo) T1D (e.g. Bacteroides spp)[ 31 ]. While differences in altered butyrate production between T1D and controls have been reported, fecal SCFA levels did not differ between groups. However, SCFA levels poorly reflect biologically active SCFA levels, as only 5±10% is found in feces and plasma[ 32 ], and the majority of microbiota-derived butyrate is consumed by the intestinal epithelium. Moreover after SCFAs have been formed in the intestine, they are subjected to a high intestinal inter-conversion due to microbial cross-feeding and the differential pathways via which these SCFA are 9 / 14 Fig 4. Most discriminating oral species. Relative abundances of most discriminating oral microbiota resulting from our elastic net analysis for (A) T1D patients and (B) healthy controls (B). Please note that the yaxes differ between figures and that Streptococcus abundance was relatively very large (0.53 in T1D and 0.43 in healthy controls) and did not fit in the figure. further metabolized have been only partly unveiled24. Nevertheless, the decreased amounts of plasma acetate and propionate in T1D might be of interest, since propionate has been linked to endogenous hepatic glucose production[ 33 ], which was previously reported to be higher in T1D[34]. The discordant correlations between propionate and acetate in both T1D and con trols (S2 Fig) may thus point towards relevant differences in acetate metabolism. In this regard it is noteworthy that acetate can protect from the development of anti-islet cell autoantibodies [ 35 ], while serum acetate directly influences CD8+ memory T cell function[ 36 ]. Finally, it has previously been recognized that altered oral health is observed in T1D and is invariably related to glucose regulation[ 37 ]. In contrast to fecal microbiota, we observed a much stronger oral microbiota signature (including Streptococcus) in T1D subjects. Interestingly, pancreatic duct infusion with Streptococci is linked to CD43+ accumulation and beta cell inflammation in animals[ 4 ]. Since oral microbiota are claimed to highly resemble the small 10 / 14 intestinal microbiome[ 15 ], it is tempting to speculate that these oral bacterial strains may indeed contribute to T1D development. Our study was performed in T1D patients without complications. Since this comprises only one third of all T1D subjects, it precludes direct extrapolation of our findings to other subgroups of T1D patients. Another caveat involves the notion that diet can influence gut[ 38 ] and oral microbiota composition. In spite of close monitoring of caloric intake in both groups, a trend towards different fiber intake in T1D was noted that might introduce some bias. Finally, there might be a reciprocal relation between hyperglycemia and altered microbiota composition[ 39 ]. However, our population had good glycemic control without diabetes complications and no correlations were seen between significantly different oral or fecal taxa and HbA1c. Finally, as we did not perform HLA typing, we cannot provide data on the potential relation between HLA type and intestinal microbiota composition In T1D. In conclusion, our study reveals that intestinal microbiota are different in well-controlled subjects with longstanding T1D compared to matched healthy controls. Moreover, correlations between markers of gut inflammation, endotoxemia and glycemic control in relation to SCFA were observed in T1D subjects suggesting potential pathophysiological links in T1D pathogenesis. In this respect, we are currently executing a fecal microbiota transplantation study in new onset T1D which may provide more insight into potential causality of (small) intestinal microbiota and their metabolites[ 40 ] in this autoimmune disease. Especially the observed differences in oral microbiota are of interest in view of recent human studies suggesting a link between oral microbiota and microbiota in the small intestine, which is believed to be the culprit site in T1D development[ 41 ]. Supporting information S1 Fig. SCFA correlations. The figure shows correlations between plasma and fecal SCFA's in T1D subjects and healthy controls. Of note, acetate is differentially correlated to propionate and butyrate in T1D vs healthy controls. <0.05, <0.01. (PPTX) S2 Fig. Correlations of key microbes with other parameters. Correlations of key microbes are depicted with clinical and inflammatory parameters and SCFA for T1D (A) and healthy controls (B). (+) indicates increased abundance in T1D, (-) indicates reduced abundance in T1D. Taxa Erysipelotrichaceae and Subdoligranulum were not significantly different between groups, but correlated significantly with a large number of clinical parameters. (PDF) S1 Database File. Excel database file containing raw data. (CSV) Acknowledgments The authors would like to thank Cees Rustemeijer, Teri Brouwer, Sytze van Dam, Victor Gerdes, Eelco Meesters and Frank Stam for inclusion of patients, Alinda Schimmel for techni cal assistance in the lab work, Lennart Kleinjans for advice on microbiota analyses and all type 1 diabetes patients and healthy control subjects for their enthusiastic participation. Author Contributions Conceptualization: Pieter F. de Groot, Frits Holleman, Willem M. de Vos, Max Nieuwdorp. Data curation: Pieter F. de Groot, OÈmruÈm Aydin, Suat Simsek, Max Nieuwdorp. 11 / 14 Formal analysis: Pieter F. de Groot, Evgeni Levin, Johannes H. Levels, Steven Aalvink, Fransje Boot, Torsten P. Scheithauer, Frank G. Schaap, Steven W. M. Olde Damink. Funding acquisition: Willem M. de Vos, Max Nieuwdorp. Investigation: Pieter F. de Groot, OÈ mruÈm Aydin. Methodology: Pieter F. de Groot, Willem M. de Vos, Max Nieuwdorp. Project administration: Pieter F. de Groot, Max Nieuwdorp. Supervision: Pieter F. de Groot, Clara Belzer, Willem M. de Vos, Max Nieuwdorp. Validation: Pieter F. de Groot. Visualization: Pieter F. de Groot, Evgeni Levin. Writing ± original draft: Pieter F. de Groot. Writing ± review & editing: Clara Belzer, Frits Holleman, DanieÈl H. van Raalte, Frank G. Schaap, Bart O. Roep, Joost B. Hoekstra, Willem M. de Vos, Max Nieuwdorp. 12 / 14 13 / 14 1. Virtanen SM , Nevalainen J , Kronberg-KippilaÈ C , Ahonen S , Tapanainen H , Uusitalo L , et al. Food consumption and advanced β cell autoimmunity in young children with HLA-conferred susceptibility to type 1 diabetes: a nested case-control design . Am J Clin Nutr . 2012 ; 95 : 471 ±8. https://doi.org/10.3945/ajcn. 111.018879 PMID: 22237062 2. Bach J-F . The effect of infections on susceptibility to autoimmune and allergic diseases . N Engl J Med . Massachusetts Medical Society; 2002 ; 347 : 911 ± 20 . https://doi.org/10.1056/NEJMra020100 PMID: 12239261 3. Clausen TD , Bergholt T , Eriksson F , Rasmussen S , Keiding N , Løkkegaard EC . Prelabor Cesarean Section and Risk of Childhood Type 1 Diabetes: A Nationwide Register-based Cohort Study . Epidemiology . 2016 ; 27 : 547 ± 55 . https://doi.org/10.1097/EDE.0000000000000488 PMID: 27031040 4. Korsgren S , Molin Y , Salmela K , Lundgren T , Melhus Å , Korsgren O . On the Etiology of Type 1 Diabetes: A New Animal Model Signifying a Decisive Role for Bacteria Eliciting an Adverse Innate Immunity Response . Am J Pathol . 2012 ; 181 : 1735 ± 1748 . https://doi.org/10.1016/j.ajpath. 2012 . 07 .022 PMID: 22944599 5. Kostic AD , Gevers D , Siljander H , Vatanen T , HyoÈtylaÈinen T , HaÈmaÈlaÈinen A-M , et al. The Dynamics of the Human Infant Gut Microbiome in Development and in Progression toward Type 1 Diabetes . Cell Host Microbe . 2015 ; 17 : 260 ± 273 . https://doi.org/10.1016/j.chom. 2015 . 01 .001 PMID: 25662751 6. Endesfelder D , zu Castell W , Ardissone A , Davis-Richardson AG , Achenbach P , Hagen M , et al. Compromised gut microbiota networks in children with anti-islet cell autoimmunity . Diabetes . 2014 ; 63 : 2006 ±14. https://doi.org/10.2337/db13-1676 PMID: 24608442 7. de Goffau MC , LuopajaÈrvi K , Knip M , Ilonen J , Ruohtula T , HaÈrkoÈnen T , et al. Fecal microbiota composition differs between children with β-cell autoimmunity and those without . Diabetes . American Diabetes Association; 2013 ; 62 : 1238 ± 44 . https://doi.org/10.2337/db12-0526 PMID: 23274889 8. Marino E , Richards JL , McLeod KH , Stanley D , Yap YA , Knight J , et al. Gut microbial metabolites limit the frequency of autoimmune T cells and protect against type 1 diabetes . Nat Immunol . Nature Publishing Group, a division of Macmillan Publishers Limited . All Rights Reserved. ; 2017;advance on . Available: http://dx.doi.org/10.1038/ni.3713 9. Stewart CJ , Nelson A , Campbell MD , Walker M , Stevenson EJ , Shaw JA , et al. Gut microbiota of Type 1 diabetes patients with good glycaemic control and high physical fitness is similar to people without diabetes: an observational study . Diabet Med . 2016 ; https://doi.org/10.1111/dme.13140 PMID: 27100052 10. Koren O , Spor A , Felin J , Fåk F , Stombaugh J , Tremaroli V , et al. Human oral, gut, and plaque microbiota in patients with atherosclerosis . Proc Natl Acad Sci U S A . 2011 ; 4592 ±8. https://doi.org/10.1073/ pnas.1011383107 PMID: 20937873 11. Qin N , Yang F , Li A , Prifti E , Chen Y , Shao L , et al. Alterations of the human gut microbiome in liver cirrhosis . Nature . 2014 ; 513 : 59 ± 64 . https://doi.org/10.1038/nature13568 PMID: 25079328 12. Kootte RS , Vrieze A , Holleman F , Dallinga-Thie GM , Zoetendal EG , de Vos WM , et al. The therapeutic potential of manipulating gut microbiota in obesity and type 2 diabetes mellitus . Diabetes Obes Metab . 2012 ; 14 : 112 ± 20 . https://doi.org/10.1111/j.1463- 1326 . 2011 . 01483 . x PMID : 21812894 13. van Eijk HMH , Bloemen JG , Dejong CHC . Application of liquid chromatography-mass spectrometry to measure short chain fatty acids in blood . J Chromatogr B Analyt Technol Biomed Life Sci . 2009 ; 877 : 719 ± 24 . https://doi.org/10.1016/j.jchromb. 2009 . 01 .039 PMID: 19230798 14. De Baere S , Eeckhaut V , Steppe M , De Maesschalck C , De Backer P , Van Immerseel F , et al. Development of a HPLC-UV method for the quantitative determination of four short-chain fatty acids and lactic acid produced by intestinal bacteria during in vitro fermentation . J Pharm Biomed Anal . 2013 ; 80 : 107 ± 15 . https://doi.org/10.1016/j.jpba. 2013 . 02 .032 PMID: 23542733 15. van den Bogert B , Erkus O , Boekhorst J , de Goffau M , Smid EJ , Zoetendal EG , et al. Diversity of human small intestinal Streptococcus and Veillonella populations . FEMS Microbiol Ecol . 2013 ; 85 : 376 ± 88 . https://doi.org/10.1111/ 1574 - 6941 .12127 PMID: 23614882 16. Ramiro-Garcia J , Hermes GDA , Giatsis C , Sipkema D , Zoetendal EG , Schaap PJ , et al. NG-Tax, a highly accurate and validated pipeline for analysis of 16S rRNA amplicons from complex biomes . F1000Research . 2016 ; 5 : 1791 . https://doi.org/10.12688/f1000research.9227.1 17. Louis P , Flint HJ . Development of a Semiquantitative Degenerate Real-Time PCR-Based Assay for Estimation of Numbers of Butyryl-Coenzyme A (CoA) CoA Transferase Genes in Complex Bacterial Samples . Applied and Environmental Microbiology . 2007 . pp. 2009 ± 2012 . https://doi.org/10.1128/ AEM.02561-06 PMID: 17259367 18. Zou H , Hastie T. Regularization and variable selection via the elastic net . J R Stat Soc Ser B (Statistical Methodol . Blackwell Publishing Ltd; 2005 ; 67 : 301 ± 320 . https://doi.org/10.1111/j.1467- 9868 . 2005 . 00503 .x 19. Tibshirani R . Regression shrinkage and selection via the lasso: a retrospective . J R Stat Soc Ser B (Statistical Methodol . Blackwell Publishing Ltd; 2011 ; 73 : 273 ± 282 . https://doi.org/10.1111/j.1467- 9868 . 2011 . 00771 .x 20. Gale EA , Gillespie KM , Dabelea D , Mayer-Davis EJ , Saydah S , Imperatore G , et al. Time trends in the incidence of type 1 diabetes in Finnish children: a cohort study . Lancet (London, England). Germany; 2014 ; 311 : 1777 ± 82 . https://doi.org/10.1007/s001250051573 21. Dabelea D , Mayer-Davis EJ , Saydah S , Imperatore G , Linder B , Divers J , et al. Prevalence of type 1 and type 2 diabetes among children and adolescents from 2001 to 2009 . JAMA. United States; 2014 ; 311 : 1778 ± 1786 . https://doi.org/10.1001/jama. 2014 .3201 PMID: 24794371 22. Hansen CHF , Krych L , Nielsen DS , Vogensen FK , Hansen LH , Sørensen SJ , et al. Early life treatment with vancomycin propagates Akkermansia muciniphila and reduces diabetes incidence in the NOD mouse . Diabetologia . 2012 ; 55 : 2285 ± 94 . https://doi.org/10.1007/s00125-012 -2564-7 PMID: 22572803 23. Clausen TD , Bergholt T , Bouaziz O , Arpi M , Eriksson F , Rasmussen S , et al. Broad-Spectrum Antibiotic Treatment and Subsequent Childhood Type 1 Diabetes: A Nationwide Danish Cohort Study . Bouchama A, editor . PLoS One . 2016 ; 11 : e0161654. https://doi.org/10.1371/journal.pone. 0161654 PMID: 27560963 24. Goodrich JK , Waters JL , Poole AC , Sutter JL , Koren O , Blekhman R , et al. Human genetics shape the gut microbiome . Cell. NIH Public Access; 2014 ; 159 : 789 ± 99 . https://doi.org/10.1016/j.cell. 2014 . 09 .053 PMID: 25417156 25. Chumpitazi BP , Hollister EB , Oezguen N , Tsai CM , McMeans AR , Luna RA , et al. Gut microbiota influences low fermentable substrate diet efficacy in children with irritable bowel syndrome . Gut Microbes . United States; 2014 ; 5 : 165 ± 175 . https://doi.org/10.4161/gmic.27923 PMID: 24637601 26. Horwitz MS , Bradley LM , Harbertson J , Krahl T , Lee J , Sarvetnick N. Diabetes induced by Coxsackie virus: initiation by bystander damage and not molecular mimicry . Nat Med . 1998 ; 4 : 781 ± 5 . Available: http://www.ncbi.nlm.nih.gov/pubmed/9662368 PMID: 9662368 27. Wen L , Ley RE , Volchkov PY , Stranges PB , Avanesyan L , Stonebraker AC , et al. Innate immunity and intestinal microbiota in the development of Type 1 diabetes . Nature . Nature Publishing Group; 2008 ; 455 : 1109 ± 13 . https://doi.org/10.1038/nature07336 PMID: 18806780 28. Burrows MP , Volchkov P , Kobayashi KS , Chervonsky A V . Microbiota regulates type 1 diabetes through Toll-like receptors . Proc Natl Acad Sci . 2015 ; 112 : 9973 ± 9977 . https://doi.org/10.1073/pnas. 1508740112 PMID: 26216961 29. Costa FRC , FrancËozo MCS , de Oliveira GG , Ignacio A , Castoldi A , Zamboni DS , et al. Gut microbiota translocation to the pancreatic lymph nodes triggers NOD2 activation and contributes to T1D onset . J Exp Med . 2016 ; 213 : 1223 ± 1239 . https://doi.org/10.1084/jem.20150744 PMID: 27325889 30. Vatanen T , Kostic AD , d'Hennezel E , Siljander H , Franzosa EA , Yassour M , et al. Variation in Microbiome LPS Immunogenicity Contributes to Autoimmunity in Humans . Cell. Elsevier; 2016 ; 165 : 842 ± 853 . https://doi.org/10.1016/j.cell. 2016 . 04 .007 PMID: 27133167 31. Planer JD , Peng Y , Kau AL , Blanton L V. , Ndao IM , Tarr PI , et al. Development of the gut microbiota and mucosal IgA responses in twins and gnotobiotic mice . Nature . Nature Publishing Group; 2016 ; 534 : 263 ± 266 . https://doi.org/10.1038/nature17940 PMID: 27279225 32. Canfora EE , Jocken JW , Blaak EE . Short-chain fatty acids in control of body weight and insulin sensitivity . Nat Rev Endocrinol. Nature Research ; 2015 ; 11 : 577 ± 591 . https://doi.org/10.1038/nrendo. 2015 .128 PMID: 26260141 33. De Vadder F , Kovatcheva-Datchary P , Goncalves D , Vinera J , Zitoun C , Duchampt A , et al. Microbiotagenerated metabolites promote metabolic benefits via gut-brain neural circuits . Cell . 2014 ; 156 : 84 ± 96 . https://doi.org/10.1016/j.cell. 2013 . 12 .016 PMID: 24412651 34. PRATO SDEL , NOSADINI R , TIENGO A , TESSARI P , AVOGARO A , TREVISAN R , et al. Insulin-Mediated Glucose Disposal in Type I Diabetes: Evidence for Insulin Resistance . J Clin Endocrinol Metab . 1983 ; 57 : 904 ± 910 . https://doi.org/10.1210/jcem-57-5-904 PMID: 6352727 35. Endesfelder D , Engel M , Davis-Richardson AG , Ardissone AN , Achenbach P , Hummel S , et al. Towards a functional hypothesis relating anti-islet cell autoimmunity to the dietary impact on microbial communities and butyrate production . Microbiome . 2016 ; 4 : 17 . https://doi.org/10.1186/s40168-016 - 0163-4 PMID: 27114075 36. Balmer ML , Ma EH , Bantug GR , GraÈhlert J , Pfister S , Glatter T , et al. Memory CD8(+) T Cells Require Increased Concentrations of Acetate Induced by Stress for Optimal Function . Immunity . 2016 ; https:// doi.org/10.1016/j.immuni. 2016 . 03 .016 PMID: 27212436 37. Lalla E , Papapanou PN . Diabetes mellitus and periodontitis: a tale of two common interrelated diseases . Nat Rev Endocrinol . 2011 ; 7 : 738 ± 48 . https://doi.org/10.1038/nrendo. 2011 .106 PMID: 21709707 38. Zeevi D , Korem T , Zmora N , Israeli D , Rothschild D , Weinberger A , et al. Personalized Nutrition by Prediction of Glycemic Responses . Cell . 2015 ; 163 : 1079 ± 1094 . https://doi.org/10.1016/j.cell. 2015 . 11 .001 PMID: 26590418 39. Hartstra A V. , Bouter KEC , Backhed F , Nieuwdorp M. Insights Into the Role of the Microbiome in Obesity and Type 2 Diabetes . Diabetes Care . 2014 ; 38 : 159 ± 165 . https://doi.org/10.2337/dc14-0769 PMID: 25538312 40. Heintz-Buschart A , May P , Laczny CC , Lebrun LA , Bellora C , Krishna A , et al. Integrated multi-omics of the human gut microbiome in a case study of familial type 1 diabetes . Nat Microbiol . Nature Publishing Group; 2016 ; 2 : 16180 . https://doi.org/10.1038/nmicrobiol. 2016 .180 PMID: 27723761 41. Vaarala O. Is the origin of type 1 diabetes in the gut[quest] . Immunol Cell Biol . Australasian Society for Immunology Inc .; 2012 ; 90 : 271 ± 276 . Available: http://dx.doi.org/10.1038/icb. 2011 .115 PMID: 22290506


This is a preview of a remote PDF: http://journals.plos.org/plosone/article/file?id=10.1371/journal.pone.0188475&type=printable

Pieter F. de Groot, Clara Belzer, Ömrüm Aydin, Evgeni Levin, Johannes H. Levels, Steven Aalvink, Fransje Boot, Frits Holleman, Daniël H. van Raalte, Torsten P. Scheithauer, Suat Simsek, Frank G. Schaap, Steven W. M. Olde Damink, Bart O. Roep, Joost B. Hoekstra, Willem M. de Vos, Max Nieuwdorp. Distinct fecal and oral microbiota composition in human type 1 diabetes, an observational study, PLOS ONE, 2017, DOI: 10.1371/journal.pone.0188475