DiabetesCare Metformin Is Associated With acobo de lo Cuesta-Zuluaga, Higher Relative abundance of tes-Agudelo, Eliang p velasquez-Meiig Mucin-Degrading akkermansia 从d muciniphila and Several Short- Chain Fatty Acid-Producing Microbiota in the Gut D0t:10.2337/dc16-1324 OBJECTIVE t the beneficial effects of metformin on glu memicroballymeddWemnd tection of type2dabetes metformin,and gut microbiota in community-dwelling Colombian adults.On the basis of previous research,we hypoth obiogaofshortehn sized that metformin is associated with fatty acid(SCFA)-producing and m cin-degrading RESEARCH DESIGN AND METHODS Participants were selected from a larger cohort of 459 participants.The present analyses focus on the 28 participants diagnosed with diabetes-14 taking metformin- and the 84 participants without diabetes (3-to-1)to P pants with di betes by s age.am We measured de analyze the composition and structure of the gut microbiota. RESULTS We found an association between diabetes and gut microbiota that was modified by metformin use nts with participants known for m rin de adation and s biota known for production of SCFAs,including Butyrivib Bifid acterium bifidum,Megasph and an ope rati n out diabetes pa 20 June 2016 and accepted 27 Septem 02dos and a distinct pperational taxonomis unit of Prevptell anda lowe r206 abundance of Enterococcus casseliflavus. CONCLUSIONS 2016 by the A Our results support the hypothesis that metformin shifts gut microbiota compo 北as ng as the wo n through the enr n-degrading A.n ld as w for profit.。 shifts mediate met min's glyc emic and anti-infammatory p le ot http://www. Diabetes Care Publish Ahead of Print,published online November 14,2016
Metformin Is Associated With Higher Relative Abundance of Mucin-Degrading Akkermansia muciniphila and Several ShortChain Fatty Acid–Producing Microbiota in the Gut DOI: 10.2337/dc16-1324 OBJECTIVE Recent studies suggest the beneficial effects of metformin on glucose metabolism may be microbially mediated. We examined the association of type 2 diabetes, metformin, and gut microbiota in community-dwelling Colombian adults. On the basis of previous research, we hypothesized that metformin is associated with higher levels of short-chain fatty acid (SCFA)–producing and mucin-degrading microbiota. RESEARCH DESIGN AND METHODS Participants were selected from a larger cohort of 459 participants. The present analyses focus on the 28 participants diagnosed with diabetesd14 taking metformind and the 84 participants without diabetes who were matched (3-to-1) to participants with diabetes by sex, age, and BMI. We measured demographic information, anthropometry, and blood biochemical parameters and collected fecal samples from which we performed 16S rRNA gene sequencing to analyze the composition and structure of the gut microbiota. RESULTS We found an association between diabetes and gut microbiota that was modified by metformin use. Compared with participants without diabetes, participants with diabetes taking metformin had higher relative abundance of Akkermansia muciniphila, a microbiota known for mucin degradation, and several gut microbiota known for production of SCFAs, including Butyrivibrio, Bifidobacterium bifidum, Megasphaera, and an operational taxonomic unit of Prevotella. In contrast, compared with participants without diabetes, participants with diabetes not taking metformin had higher relative abundance of Clostridiaceae 02d06 and a distinct operational taxonomic unit of Prevotella and a lower abundance of Enterococcus casseliflavus. CONCLUSIONS Our results support the hypothesis that metformin shifts gut microbiota composition through the enrichment of mucin-degrading A. muciniphila as well as several SCFA-producing microbiota. Future studies are needed to determine if these shifts mediate metformin’s glycemic and anti-inflammatory properties. 1 VidariumdNutrition, Health and Wellness Research Center, Grupo Empresarial Nutresa, Medellin, Colombia 2 Department of Epidemiology, Johns Hopkins Bloomberg School of Public Health, Baltimore, MD 3 Dinamica I.P.S. ´ dEspecialista en Ayudas Diagnosticas, Medellin, Colombia ´ 4 EPS y Medicina Prepagada Suramericana S.A., Medellin, Colombia Corresponding author: Juan S. Escobar, jsescobar@ serviciosnutresa.com. Received 20 June 2016 and accepted 27 September 2016. This article contains Supplementary Data online at http://care.diabetesjournals.org/lookup/ suppl/doi:10.2337/dc16-1324/-/DC1. © 2016 by the American Diabetes Association. Readers may use this article as long as the work is properly cited, the use is educational and not for profit, and the work is not altered.More information is available at http://www.diabetesjournals .org/content/license. Jacobo de la Cuesta-Zuluaga,1 Noel T. Mueller,2 Vanessa Corrales-Agudelo,1 Eliana P. Velasquez-Mej ´ ´ıa,1 Jenny A. Carmona,3 Jose M. Abad, ´ 4 and Juan S. Escobar1 Diabetes Care 1 EMERGING TECHNOLOGIES AND THERAPEUTICS Diabetes Care Publish Ahead of Print, published online November 14, 2016
2 Type 2 Diabetes,Metformin,and Gut Microbiota Diabetes Car (11-13).a mucin-degrading bacteria nt.and individuals diae their associated genes(microbiome) with Alzheimer disease,Parkinson dis ute the T our env In numans,parti r any other ne rative d ved to deterministic role in had similar abundance of Subdoliaran and gastrointestinal diseases (crohn dis lum and Akkerm ansi ease.ulcerative colitis short bo wel syn of type .or diah globe,has been linked in nonhuman (1) may help ameliorate a type 2 diabetes dance with the principles of the Dec ation of Helsink the nost rerent hu takir elative abundance of Adlercreutzia(17). and 008430 of 1993).All of the partici don was bn a poten metage ants were th nt (6) and pron e assured of anonymity and conf Metformi 1.1 imethylbigua dia ontrast tiality nfor wa th ith with had Th (7).and findings from recent studi dance of fubacte nmittee of Sede de cio t ma o pre 53 17 min has nleiot ic effe vet the Subdoli anulum and a and an ved the r s de echa nist cluster of butyrate- oducing Clostr 588 ings pro alters hepatic glucose production via to the antidiabe effects of metformin metric,Clinical,and Dietary Evalu through 03 that mucin we BM (kg)/heigh esi the gut (11 n this lean (18.5 BMI ke/m 509% eralizability of ght B/m ng t inte inal an -30 21 ation of metformin is 300 times that of of HDI IDI VIDI total choles dysbiosis in a colombian adult popula iglycerides,apolipoprotein B asting the plasma (16),mal on.Given the conside ariation in etes and that the obtained (collection and m administration of metformin lombians is differen to tha of othe nt explained in Supplemen y Data a14 518 hyp n the 24-h 。e throug ns a Data) al me ormin a rogates in s ferent to those obs and 2 ction and enc effect of metfo nin enrichment of mucin-degrading and scFA Each participant collected their own fe- ela producing micro ealed,ste the RESEARCH DESIGN AND METHODS receptacle ffect than imm r ext Study Design erated in h sehold freezers and brough e July and N ember 2 ty in ea xposure than the other metformin for ld with BMI =18 5 kg/m2 living in sively in the morning (6 AM-12 PM) mulations(15 the Colom of M n,bt ted the anima 11 2,1 16171 nlled in the dav ay that metformin may partially restore sured by the health insu rance prov ide samples were stored on dry ice and ith type next-da formin treatment increased the relativ individual who consu anti evaluated by trained lab abundance of Akkermansia muciniphilo biotics or antiparasitics <3 months prior oratory technicians
Microbial communities (microbiota) and their associated genes (microbiome) constitute the interface of our environs and our cells, and their composition is believed to play a deterministic role in human health and disease. In particular, the development of type 2 diabetes, a disease rising in prevalence around the globe, has been linked in nonhuman (1) and human (2–5) studies to imbalances in microbiota of the intestinal tract (gut). However, the most recent human study on this topic found that the association was modulated in a potentially beneficial manner by metformin treatment (6). Metformin (1,1-dimethylbiguanide hydrochloride) is the most frequent medication used to treat type 2 diabetes (7), and findings from recent studies suggest it may also prevent cancer (7) and cardiovascular events (8). Metformin has pleiotropic effects, yet the majority of mechanistic studies have focused on changes in liver function (7,9,10). Although metformin certainly alters hepatic glucose production via effects on AMPK, there is growing evidence that the genesis of its action is in the gut (11–15). Metformin is ;50% bioavailable, allowing for near-equal intestinal and plasma exposure, but intestinal accumulation of metformin is 300 times that of the plasma (16), making the gut the primary reservoir for metformin in humans. Unlike oral administration, intravenous administration of metformin in humans does not improve glycemia (14). Moreover, in mice, oral administration of a broad-spectrum antibiotic cocktail with oral metformin abrogates metformin’s glucose-lowering effect (12). Providing yet further evidence that the glucoselowering effect of metformin may originate in the lower bowel, a delayed-release oral metformin, which targets the ileum, had a similar or greater glucose-lowering effect than immediate-release or extendedrelease metformin, despite the delayedrelease metformin having lower systemic exposure than the other metformin formulations (15). Recent studies in animals (11,12,13) and humans (6,17) provide evidence that metformin may partially restore gut dysbiosis associated with type 2 diabetes. In mice fed a high-fat diet, metformin treatment increased the relative abundance of Akkermansia muciniphila (11–13), a mucin-degrading bacteria that has been shown to reverse metabolic disorders (1,12). In humans, participants with diabetes taking metformin had similar abundance of Subdoligranulum and, to some extent, Akkermansia compared with control subjects without diabetes, suggesting that metformin may help ameliorate a type 2 diabetes– associated gut microbiome (6). It has also been shown that people with diabetes taking metformin had a higher relative abundance of Adlercreutzia (17), and metagenomic functional analyses demonstrated significantly enhanced butyrate and propionate production in people with diabetes using metformin (6). In contrast, people with diabetes who were not treated with metformin had a higher abundance of Eubacterium and Clostridiaceae SMB53 (17) and lower levels of shortchain fatty-acid (SCFA)–producers, such as Roseburia, Subdoligranulum, and a cluster of butyrate-producing Clostridiales (6). These findings provide evidence that gut microbes may contribute to the antidiabetes effects of metformin through pathways that include mucin degradation and SCFA production. In this study we aimed to test the generalizability of previous observations concerning the influence of metformin on the association of type 2 diabetes and gut dysbiosis in a Colombian adult population. Given the considerable variation in the microbiota associated with type 2 diabetes and that the gut microbiota of Colombians is different to that of other populations (18), we hypothesized that the microbial taxa involved in the type 2 diabetes dysbiosis of Colombians are different to those observed in Chinese and European populations (2,3) but that the effect of metformin is similar, i.e., through enrichment of mucin-degrading and SCFAproducing microbiota. RESEARCH DESIGN AND METHODS Study Design Between July and November 2014, we enrolled 459 men and women 18–62 years old, with BMI $18.5 kg/m2 , living in the Colombian cities of Medellin, Bogota, Barranquilla, Bucaramanga, and Cali. All participants enrolled in the study were insured by the health insurance provider EPS y Medicina Prepagada Suramericana S.A. (EPS SURA). We excluded pregnant women, individuals who consumed antibiotics or antiparasitics ,3 months prior to enrollment, and individuals diagnosed with Alzheimer disease, Parkinson disease, or any other neurodegenerative diseases; current or recent cancer (,1 year); and gastrointestinal diseases (Crohn disease, ulcerative colitis, short bowel syndrome, diverticulosis, or celiac disease). This study was conducted in accordance with the principles of the Declaration of Helsinki, as revised in 2008, and had minimal risk according to the Colombian Ministry of Health (Resolution 008430 of 1993). All of the participants were thoroughly informed about the study and procedures. Participants were assured of anonymity and confi- dentiality. Written informed consent was obtained from all the participants before beginning the study. The Bioethics Committee of Sede de Investigacion´ UniversitariadUniversity of Antioquia reviewed the protocol and the consent forms and approved the procedures described here (approbation act 14–24–588 dated 28 May 2014). Anthropometric, Clinical, and Dietary Evaluations We calculated BMI as weight (kg)/height squared (m2 ) to classify participants as lean (18.5 # BMI , 25.0 kg/m2 ), overweight (25.0 # BMI , 30.0 kg/m2 ), or obese (BMI $30 kg/m2 ). In addition, values of HDL, LDL, VLDL, total cholesterol, triglycerides, apolipoprotein B, fasting glucose, glycated hemoglobin (HbA1c), fasting insulin, adiponectin, and hs-CRP were obtained (collection and measurement explained in Supplementary Data). Dietary intake was evaluated through 24-h dietary recalls (see Supplementary Data). DNA Extraction and Sequence Analysis Each participant collected their own fecal sample in a hermetically sealed, sterile receptacle provided by the research team. Samples were immediately refrigerated in household freezers and brought to an EPS SURA facility in each city within 12 h; receipt of samples occurred exclusively in the morning (6 A.M.–12 P.M.). As such, stools were collected between the evening of the day before and the morning of the day of sample receipt. Fecal samples were stored on dry ice and sent to a central laboratory via next-day delivery. Before DNA extraction, stool consistency was evaluated by trained laboratory technicians. 2 Type 2 Diabetes, Metformin, and Gut Microbiota Diabetes Care
care.diabete de la Cuesta-Zuluaga and Associates sing 14 with ty (Qiagen,Hilden,Germany)following metformin (T2D-met"),14 with type 2 di- plemented in the Vegan package of R(22). the es (ND (T2-met) We ine be ingstep with the lysis buffer (0at agnostically identify microbial bic 15 Hz using a stainless stee with a stical An markers.LEfS e uses the nonparam d clinical variables ified DNA e compared across study groups us e noDrop spectrophotometer (Ny ing ANOVA and t tests after check ing fo differential abundance amon ote ns,Fra nce)and sent it to th par then p of Michigan Medical School (Ann Art ants using the unpaired wil MI).The V4 hypervariable region of the lity tests sum test,and final lly uses LD est transformed (natural log for unco and R806 (S'-GGACTACHVGGGTWTC ined variables or ar quare roo of lEfse com with sta statis -3)p rs and equenced sing for p rati tica approaches is that, y and the Statistical analyses were formed ing strategy(19).n addition t with Rv..3.2.2(21】 een each OTU and the grou g cat we sequ the QIamp DNA S 's EB el median 28,699).To limit the effects o nicrobial biomarkers in ou tio an hhe da etained f they had a A 0.05 and 0g10 D M e 12D-me na ment's run. quences were curated ugh Finally,across groups,we tested for dif foll ving th eg sta a erences in e the likeli 13620 of detecting false positives.especially maior butyrate. microbia seq tedat the NCBIand mic units a,including Butyrivibni can be th he a nd com as assessed by quantifying and ssified in each of thes ion of Type 2 iab and preting (4 for A.muciniphil 11 for that had typedi;6self-reported and the erence using ANOVA nd t te te t n We compared the nine whet (fasting bloc 2126 mg/d ne analysis。 bunda of the 28 only ir formin treatment,14 not (1 was Diversity using Results from lEfSe and from the pooled treat with insutin alo with gli d g mu of ph otype e corrected fo logical treat nent for =0 5)calcr including the 2 participants unaware of ated with the GUniF ge ofR(23 age of R(25).Tests were considered betes status)(supplementary a P value. We matched each participant with di nce per OTU.then obtained a dis with thr matrix from pairwis (to th finally eda rela cf20nerP0iemtheohg ndrNSteristiG difference between case and control joining phylogenetic tree using Mothur pants.There were no statistically signif mong en es (a ght or ohese)This left us with atotal sing the adonis function (ANOVA using cal narameters hetw en T2D-met*and
Total microbial DNA was extracted using the QIAamp DNA Stool Mini Kit (Qiagen, Hilden, Germany) following the manufacturer’s instructions, with a slight modification consisting in a beadbeating step with the lysis buffer (20 s at 15 Hz using a stainless steel bead with a 5-mm diameter). After extraction, we quantified DNA concentration using a NanoDrop spectrophotometer (Nyxor Biotech, Paris, France) and sent it to the Microbial Systems Laboratory, University of Michigan Medical School (Ann Arbor, MI). The V4 hypervariable region of the 16S rRNA gene was amplified using the F515 (59-CACGGTCGKCGGCGCCATT-39) and R806 (59-GGACTACHVGGGTWTC TAAT-39) primers and sequenced using the Illumina MiSeq sequencing platform with V2 chemistry and the dual-index sequencing strategy (19). In addition to DNA from fecal samples, we sequenced negative controls (ultrapure water and the QIAamp DNA Stool Mini Kit’s EB elution buffer), a DNA extraction blank, and a mock community (HM-782D, BEI Resources, Manassas, VA) in each instrument’s run. Sequences were curated following the MiSeq standard operating procedure implemented by Mothur v.1.36 (20) (see Supplementary Data). Raw sequences were deposited at the NCBI and can be accessed through the BioProject (accession number PRJNA325931). Definition of Type 2 Diabetes and Selection of Control Subjects We identified 28 participants in our study that had type 2 diabetes; 26 self-reported physician-diagnosed diabetes prior to the beginning of the study and 2 were diagnosed through laboratory testing (fasting blood glucose $126 mg/dL and HbA1c$6.5%). Of the 28 participants with type 2 diabetes, 14 were under metformin treatment, 14 were not (1 was treated with insulin alone, 2 with glibenclamide, and 11 were under no pharmacological treatment for type 2 diabetes, including the 2 participants unaware of their diabetes status) (Supplementary Table 1). We matched each participant with diabetes with three participants without diabetes in our study based on sex, age (to the closest possible age; maximum difference between case and control subjects 6 years; mean 1.5 years; median 1 year), and BMI category (lean, overweight, or obese). This left us with a total analytic sample of 112 study participants comprising 14 with type 2 diabetes using metformin (T2D-met+ ), 14 with type 2 diabetes not using metformin (T2D-met2), and 84 without diabetes (ND). Statistical Analyses Anthropometric and clinical variables were compared across study groups using ANOVA and t tests after checking for homoscedasticity and normal distribution of residuals (using Fligner-Killeen tests of homogeneity of variances and Shapiro-Wilk normality tests). When necessary, variables were appropriately transformed (natural log for unconstrained variables or arcsin square root for proportions). Sex ratio and stool consistency were compared using x2 tests. Statistical analyses were performed with R v.3.2.2 (21). Curated DNA sequences ranged from 69 to 102,660 sequences per sample (median 28,699). To limit the effects of uneven sampling, we rarefied the data set to 4,091 sequences per sample, resulting in the exclusion of one T2D-met2 participant with 69 reads. Although rarefaction may lead to missing low-abundance data, it is a powerful way to reduce the likelihood of detecting false positives, especially among those operational taxonomic units (OTUs) with very low abundance. The gut microbiota structure and composition was assessed by quantifying and interpreting similarities based on intraand intergroup diversity analyses (a and b diversity, respectively). For a diversity, we calculated Good’s coverage and the number of OTUs of each sample using Mothur and constructed rarefaction curves.We compared these indices among groups of participants using analysis of similarity (ANOSIM) with 1,000 permutations using the Vegan package of R (22). b Diversity was assessed using phylogeny-based generalized UniFrac distances (with the a parameter controlling weight on abundant lineages = 0.5) calculated with the GUniFrac package of R (23). For this, we first reduced the alignment and the OTU table to one representative sequence per OTU, then obtained a distance matrix from uncorrected pairwise distances between aligned sequences, and finally constructed a relaxed neighborjoining phylogenetic tree using Mothur and Clearcut. Comparisons among groups of participants were performed using the adonis function (ANOVA using distance matrices) of the permutational multivariate ANOVA (PERMANOVA) implemented in the Vegan package of R (22). We next used linear discriminant analysis (LDA) effect size (LEfSe) (24) to agnostically identify microbial biomarkers. LEfSe uses the nonparametric factorial Kruskal-Wallis sum rank test to detect individual OTUs with significant differential abundance among groups of participants, then performs a set of pairwise tests among groups of participants using the unpaired Wilcoxon rank sum test, and finally uses LDA to estimate the effect size of each differentially abundant OTU (24). The strength of LEfSe compared with standard statistical approaches is that, in addition to providing P values, it provides an estimation of the magnitude of the association between each OTU and the grouping categories (e.g., metformin, type 2 diabetes) through the LDA score. For stringency, microbial biomarkers in our study were retained if they had a P , 0.05 and a (log10) LDA score $3, i.e., one order of magnitude greater than LEfSe’s default. Finally, across groups, we tested for differences in relative abundance of the mucin-degrading A. muciniphila and major butyrate-producing microbial genera, including Butyrivibrio, Roseburia, Subdoligranulum, and Faecalibacterium. For this analysis, we pooled all OTUs classified in each of these phylotypes (4 for A. muciniphila, 11 for Butyrivibrio, 4 for Roseburia, 10 for Subdoligranulum, and 5 for Faecalibacterium) and tested for differences using ANOVA and t tests on arcsin square root transformed relative abundances. This pooling served to examine whether differences in relative abundance of these groups of bacteria occurred across all OTUs or only in specific OTUs. Results from LEfSe and from the pooled analysis of phylotypes were corrected for multiple testing using the Bayesian approach implemented in the qvalue package of R (25). Tests were considered significant if they had a P value # 0.05 and a q value #0.1. RESULTS In Table 1 we present the characteristics of T2D-met+ , T2D-met2, and ND participants. There were no statistically significant differences (all P values . 0.10) in demographic, anthropometric, or clinical parameters between T2D-met+ and care.diabetesjournals.org de la Cuesta-Zuluaga and Associates 3
4 Type 2 Diabetes,Metformin,and Gut Microbiota Diabetes Care 18g8 8昌888838昌888 3日888B言昌君目8昌日 西始器然路1111 用Eea 888器111 日g 合器器阳9 woy sanjeA d llV.'OS ueaw se pal
Table 1 —General characteristics of T2D-met+, T2D-met2, and ND participants among community-dwelling Colombian adults Group P* T2D-met+ T2D-met2 ND T2D-met+ vs. T2D-met2 T2D-met+ vs. ND T2D-met2 vs. ND n 14 14 84 d dd Age (years) 50 6 10 44 6 9 47 6 9 0.11 0.33 0.24 Sex (F/M) 0.36 0.50 0.43 0.70 0.83 0.84 Anthropometry BMI (kg/m2 ) 31.88 6 4.63 32.15 6 6.36 31.11 6 4.53 0.98 0.56 0.62 Body fat (%) 0.41 6 0.05 0.41 6 0.03 0.40 6 0.04 0.96 0.68 0.58 Waist circumference (cm) 104.6 6 9.8 102.8 6 14.2 102.0 6 11.3 0.70 0.39 0.85 Clinical parameters Total cholesterol (mg/dL) 178 6 51 208 6 44 187 6 30 0.11 0.54 0.10 HDL (mg/dL) 40 6 11 38 6 6 44 6 12 0.77 0.21 0.0136 LDL (mg/dL) 105 6 35 128 6 39 115 6 28 0.11 0.31 0.26 VLDL (mg/dL) 38 6 40 52 6 52 31 6 15 0.19 0.64 0.0471 Apolipoprotein B (mg/dL) 102 6 30 108 6 30 97 6 26 0.56 0.55 0.18 Triglycerides (mg/dL) 176 6 202 244 6 253 154 6 75 0.15 0.97 0.07 Fasting glucose (mg/dL) 127 6 47 145 6 76 90 6 10 0.57 0.0047 0.0060 HbA1c [% (mmol/mol)] 6.9 6 1.4 (52.0 6 15.3) 7.1 6 1.7 (54.0 6 18.6) 5.6 6 0.3 (38.0 6 3.3) 0.78 0.0026 0.0047 Fasting insulin (mU/mL) 22.24 6 12.58 24.24 6 12.86 15.20 6 8.79 0.49 0.11 0.0029 Insulin resistance index 2.9 6 1.5 3.7 6 3.0 1.9 6 1.1 0.39 0.0420 0.0018 Leptin (ng/mL) 7.39 6 6.09 8.34 6 5.40 7.89 6 6.14 0.38 0.91 0.25 Adiponectin (mg/mL) 4.59 6 1.97 4.96 6 3.22 6.79 6 3.90 0.91 0.0195 0.07 hs-CRP (mg/L) 2.70 6 2.53 3.35 6 2.76 4.04 6 5.02 0.61 0.21 0.62 Dietary intake Energy intake (calories) 1,843 6 268 1,865 6 584 1,945 6 572 0.75 0.61 0.57 Carbohydrate (g) 250 6 37 260 6 86 267 6 88 0.92 0.52 0.68 Protein (g) 74.7 6 9.1 69.7 6 13.1 74.1 6 13.8 0.26 0.85 0.26 Fat (g) 60.7 6 13.0 60.2 6 21.3 62.9 6 17.1 0.64 0.76 0.51 Cholesterol (mg) 336 6 31 330 6 41 347 6 38 0.67 0.28 0.18 Dietary fiber (g) 17.9 6 4.8 18.6 6 6.4 17.5 6 4.8 0.88 0.73 0.64 Stool consistency [n (%)] 0.38 0.58 0.50 Hard 4 (28) 2 (14) 13 (15) d d d Normal 8 (57) 9 (64) 54 (64) d dd Mushy 2 (14) 1 (7) 13 (15) d dd Diarrheic 0 (0) 2 (14) 4 (5) d dd Values presented as mean 6 SD. *All P values from t tests except in sex and stool consistency (x2 tests). 4 Type 2 Diabetes, Metformin, and Gut Microbiota Diabetes Care
care.diabetesjo nals.org de la Cuesta-Zuluaga and Associates 5 T2D R2=0019 had higher fasting glucose,HbAand in- p=0.34 rent T20 16S rRNA Gene Sequencing e ND 009 g (mean Good's cov e±5D=09901 PCo1(10.34% 0.001);9%of OTUs were detected at ast bych biota.We next tested for ween participants with 20 diabetes and ND participants(A OSIM 2519 ipants (ANOSIM statistic R=-0.018,P= 0.557 The numbe of obse ved O met"and ND than bet n T2D-met (Supplementary Fig ce- PCo1(10.34% ANOSIM statistic R=0.018,P=0.348; :ANOSIM statistic R- We c s (PE 2D-me A) with dia nd ND partic =0.009,P=0.416)(Fg18) ever,the com significance (R2=0.013,P=0.036)(Fig. C).the PCo1(10.34%) De0036 0.01 T2D-met ND Participants (R 0.008,P=0.9 43).These results suggested and partici s not takir D the -me (T t-and VA.(A high-quality color represe e in the onine issue.) We next use to examine microbiota at the OTU level.Note that (represented by OTUs with [og10]LDA relevant.OTUs (19 displayed in Fig.2 and
T2D-met2 participants. Compared with ND participants, T2D-met+ participants had higher fasting glucose, HbA1c, and insulin resistance than ND participants and lower levels of the insulin-sensitizing hormone adiponectin (P , 0.05). No other demographic, anthropometric, or clinical parameters were statistically different. 16S rRNA Gene Sequencing Gut microbiota communities were specific to each participant with marked intersubject differences (Fig. 1A) (overall interindividual generalized UniFrac distance = 0.720 6 0.009). We found high coverage across all groups of participants (mean Good’s coverage 6 SD = 0.990 6 0.001); 99% of OTUs were detected at least by two DNA reads, demonstrating thorough sampling of the gut microbiota. We next tested for differences in the number of observed OTUs across the groups of participants. We found no differences between participants with diabetes and ND participants (ANOSIM statistic R = 0.005, P value = 0.425) or between T2D-met+ and T2D-met2 participants (ANOSIM statistic R = 20.018, P = 0.557). The number of observed OTUs tended to be more similar between T2Dmet+ and ND than between T2D-met2 and ND participants (Supplementary Fig. 1); however, these differences were not statistically significant (T2D-met+ vs. ND: ANOSIM statistic R = 0.018, P = 0.348; T2D-met2 vs. ND: ANOSIM statistic R = 0.009, P = 0.409). We observed no significant differences in b diversity estimates among the three groups of participants (PERMANOVA: R2 = 0.019, P = 0.335) (Fig. 1A) or between participants with diabetes and ND participants (R2 = 0.009, P = 0.416) (Fig. 1B). However, the comparison between metformin and nonmetformin users reached significance (R2 = 0.013, P = 0.036) (Fig. 1C), demonstrating differences in the bacterial community structure associated with metformin use. The difference was also significant when comparing T2Dmet+ and ND participants (R2 = 0.015, P = 0.036) but not when comparing T2D-met2 and ND participants (R2 = 0.008, P = 0.943). These results suggested the microbial communities of T2D-met+ versus T2D-met2 were modestly phylogenetically dissimilar. We next used LEfSe to examine differences in the relative abundance of gut microbiota at the OTU level. Note that we were only interested in OTUs displaying strong associations in the LDA (represented by OTUs with [log10] LDA scores $3); such stringency resulted in fewer retained, but more likely biologically relevant, OTUs (19 displayed in Fig. 2 and PCo2 (6.03%) PCo1 (10.34%) ND R2 = 0.019 P = 0.34 T2D-met+ T2D-metPCo2 (6.03%) R2 = 0.009 P = 0.42 ND T2D PCo1 (10.34%) PCo2 (6.03%) metR2 = 0.013 P = 0.036 T2D-met+ PCo1 (10.34%) A B C Figure 1—Principal coordinates analysis based on generalized UniFrac. A: Comparison among the three groups of participants. B: Comparison between participants with diabetes and ND participants. C: Comparison between T2D-met+ and participants not taking metformin (T2D-met2 and ND). Ellipses encompass 75% of data distribution in each group of participants. R2 and P values from PERMANOVA. (A high-quality color representation of this figure is available in the online issue.) care.diabetesjournals.org de la Cuesta-Zuluaga and Associates 5