5 OPEN ORIGINAL ARTICLE Structural modulation of gut microbiota during alleviation of type 2 diabetes with a Chinese herbal formula Jia Xu'4,Fengmei Lian2.4,Linhua Zhao2,Yufeng Zhao2,Xinyan Chen2,Xu Zhang' g',Qiang Zhou2,Zhengsheng Xue',Xiaoyan Pang', Liping Zhao Sciences ond Bntetnloe5hangaoa Medical Sciences.Beijing.People's Republic of China andMinistry of Education Key Laboratory of System Shanghai Center for Systems Biomedicine.Shanghai Jiao Tong University.Shanghai.People's The gut micr robiota is d to hav critical role in metabolic dise igh or the pl n 1)for 12 s in a do of fasting atedherTecl lobin HbA1c)co crobiota inre ment.This occe47b8 ched s BG whic HbAic and 2-hp dial blo a indi OD. es of gut m robiota are inc d by ch herbal formu 88coes6nein9sihthegarme8o8品amaeicemncn The ISME Joumal (2015)9,552-562:doi:10.1038/ismej.2014.177:published online 3 October 2014 Introduction resistance ent world which has an The e 39%a d proportion of insulin res nce(Hotamisade in abes etes among a among 2013: ation Key aboratory o vas and such as Escherichja urified for 200 bcutaneously infused into mic abundance the fecal samples of 171 diabetic patients and 174
OPEN ORIGINAL ARTICLE Structural modulation of gut microbiota during alleviation of type 2 diabetes with a Chinese herbal formula Jia Xu1,4, Fengmei Lian2,4, Linhua Zhao2 , Yufeng Zhao3 , Xinyan Chen2 , Xu Zhang1 , Yun Guo2 , Chenhong Zhang1 , Qiang Zhou2 , Zhengsheng Xue1 , Xiaoyan Pang1 , Liping Zhao1,3 and Xiaolin Tong2 1 State Key Laboratory of Microbial Metabolism, School of Life Sciences and Biotechnology, Shanghai Jiao Tong University, Shanghai, People’s Republic of China; 2 Guang’anmen Hospital, China Academy of Chinese Medical Sciences, Beijing, People’s Republic of China and 3 Ministry of Education Key Laboratory of Systems Biomedicine, Shanghai Center for Systems Biomedicine, Shanghai Jiao Tong University, Shanghai, People’s Republic of China The gut microbiota is hypothesized to have a critical role in metabolic diseases, including type 2 diabetes (T2D). A traditional Chinese herbal formula, Gegen Qinlian Decoction (GQD), can alleviate T2D. To find out whether GQD modulates the composition of the gut microbiota during T2D treatment, 187 T2D patients were randomly allocated to receive high (HD, n ¼ 44), moderate (MD, n ¼ 52), low dose GQD (LD, n ¼ 50) or the placebo (n ¼ 41) for 12 weeks in a double-blinded trial. Patients who received the HD or MD demonstrated significant reductions in adjusted mean changes from baseline of fasting blood glucose (FBG) and glycated hemoglobin (HbA1c) compared with the placebo and LD groups. Pyrosequencing of the V3 regions of 16S rRNA genes revealed a dosedependent deviation of gut microbiota in response to GQD treatment. This deviation occurred before significant improvement of T2D symptoms was observed. Redundancy analysis identified 47 GQDenriched species level phylotypes, 17 of which were negatively correlated with FBG and 9 with HbA1c. Real-time quantitative PCR confirmed that GQD significantly enriched Faecalibacterium prausnitzii, which was negatively correlated with FBG, HbA1c and 2-h postprandial blood glucose levels and positively correlated with homeostasis model assessment of b-cell function. Therefore, these data indicate that structural changes of gut microbiota are induced by Chinese herbal formula GQD. Specifically, GQD treatment may enrich the amounts of beneficial bacteria, such as Faecalibacterium spp. In conclusion, changes in the gut microbiota are associated with the antidiabetic effects of GQD. The ISME Journal (2015) 9, 552–562; doi:10.1038/ismej.2014.177; published online 3 October 2014 Introduction Type 2 diabetes (T2D), which is characterized by low-grade inflammation, insulin resistance (Shoelson, 2006) and b-cell failure (Butler et al., 2003), has become increasingly prevalent worldwide (Xu et al., 2013). The estimated proportion of diabetes among adults is 8.3% in 2010, among which T2D accounts for at least 90% (Alberti and Zimmet, 1998; Whiting et al., 2011). This proportion is projected to increase to 9.9% by 2030 (Whiting et al., 2011). Development of T2D results mostly from obesity, which has low-grade inflammation and insulin resistance (Hotamisligil, 2006). The gut microbiota may have a vital role in obesity development (Backhed et al., 2004; Collins et al., 2013; Le Chatelier et al., 2013; Zhao, 2013). For example, endotoxin produced by an opportunistic pathogen in the gut, such as Escherichia coli, induced obesity and insulin resistance when a purified form was subcutaneously infused into mice (Cani et al., 2007a). A greater abundance of opportunistic pathogens, such as Betaproteobacteria, was found in the gut of diabetic patients compared with healthy controls (Larsen et al., 2010). A more recent comparative metagenomic analysis of the fecal samples of 171 diabetic patients and 174 Correspondence: L Zhao, Ministry of Education Key Laboratory of Systems Biomedicine, Shanghai Center for Systems Biomedicine, Shanghai Jiao Tong University, Room 3-517, Biology Building, 800 Dongchuan Road, Shanghai 200240, People’s Republic of China or X Tong, Guang’anmen Hospital of China Academy of Chinese Medical Sciences, Room 432, Administration Building, 5 Beixiange Street, Xuanwu District, Beijing 100053, China. E-mail: lpzhao@sjtu.edu.cn or xiaolintong66@sina.com 4 These authors contributed equally to this work. Received 11 September 2013; revised 18 June 2014; accepted 15 August 2014; published online 3 October 2014 The ISME Journal (2015) 9, 552–562 & 2015 International Society for Microbial Ecology All rights reserved 1751-7362/15 www.nature.com/ismej
e 553 nt nuclear m nagneti ce-dp that the early onse 以o tabolism by degradati of cholin from the gut of morbidly obese and diabetic on dir ethe showing that GOD can modulate ota microbiota,particularly in hum doubl-blinded placebo-co oled clinical tria metabolic dis made to target the xamined the structura rine (Khin 28 1987 efficacy Materials and methods olled clinical tria tudy desig he was a 12-week.randor ized.double period.It was ap rove ea control for te Our h-fat ang'anmen hospita and insulin stanc anmen hosp pital,Dong Zhimen Hospita Beijing fron l partic ed in that treat bacterial diarrhea such a wa conducted in erberin might be useful treatment with the principles of the Declaration biota. Chin The inclusion and exclusion criteria of the has at can methods.Usi品n Its reco ed by SC reening,including FBG test 219).Su tand 75-02 nha atment for T2D were recruite into the stud observat eries of examinations.403 patients were excluded reduce glu s Th ia and?2 diabetic SD high (HD).mo one (MD)or low (LD)GQD.or the placebo for 12 weeks was low-dose treat Tonstd ies are clinical investigations with small sample The ISME Journal
healthy controls showed that diseased samples had lesser abundance of butyrate-producing bacteria, such as Faecalibacterium prausnitzii, but greater abundance of opportunistic pathogens, including Clostridium bolteae and Desulfovibrio sp. (Qin et al., 2012). Another study found that the early onset of high-fat-diet-induced T2D was characterized by an increased bacterial translocation from the intestine towards tissues (Amar et al., 2011). An opportunistic pathogen, Enterobacter cloacae B29, isolated from the gut of a morbidly obese and diabetic patient, induced obesity and insulin resistance in germ-free mice (Fei and Zhao, 2013). Taken together, these studies indicate that a dysbiotic gut microbiota may causatively contribute to obesity and diabetes development, and thus may serve as a potential new target for disease control. To treat obesity, T2D and other metabolic diseases, several attempts have been made to target the gut microbiota (Cani et al., 2007b, 2009; Park et al., 2013). Berberine, the major pharmacologic component of a Chinese herb Coptis chinensis (HuangLian) originally used to treat bacterial diarrhea (Khin Maung et al., 1985; Rabbani et al., 1987; Tang et al., 2009), showed clinical efficacy in treating diabetes in a multicentered, randomized, double-blinded and placebo-controlled clinical trial (Zhang et al., 2008). The herb C. chinensis has been used in traditional Chinese medicine (TCM) for diarrhea control for nearly 2000 years. Our recent study showed that berberine prevented high-fatdiet-induced obesity and insulin resistance, enriched short-chain fatty acid-producing bacteria, reduced numbers of opportunistic pathogens and alleviated inflammation in Wistar rats (Zhang et al., 2012b). Drugs that treat bacterial diarrhea, such as berberine, might be useful for T2D treatment because both diseases share a dysbiotic gut microbiota. A standardized berberine-containing Chinese herbal formula, Gegen Qinlian Decoction (GQD), has been a treatment for diarrhea in Shang Han Lun since the East Han Dynasty. Its use was recorded by the prestigious physician Zhongjing Zhang (AD 150– 219). Subsequently, GQD has been reported to have potentially beneficial effects in the treatment of diabetes in animal trials, as well as in some clinical observations. For example, GQD significantly reduced fasting blood glucose (FBG) and glycated hemoglobin (HbA1c) in streptozotocin (STZ) and high-fat-diet-induced diabetic SD rats, and the serum of SD rats that received GQD enhanced glucose consumption in 3T3-L1 adipocytes (Zhang et al., 2013). T2D patients treated with a high dose of modified GQD two times daily for 3 months showed a reduction in HbA1c of 1.79% from the initial level of 9.2%. This decrease was significantly different from that of patients receiving a low-dose treatment (Tong et al., 2011). However, these studies are either animal trials or open, non-placebo-controlled clinical investigations with small sample sizes. Moreover, the mechanism underlying GQD’s impact on glycemic efficacy has barely been elucidated. A recent nuclear magnetic resonance-based plasma metabonomic study revealed that 5 weeks of GQD treatment conspicuously modulated gut microbial metabolism by degradation of choline into methylamines, together with a decrease in FBG and an expansion of islets in STZ and high-fat-diet-induced diabetic rats (Tian et al., 2013). This finding indicates that the gut microbiota might have a pivotal role in the effect GQD has on diabetic subjects. However, there is still a lack of direct evidence showing that GQD can modulate gut microbiota, particularly in humans. In this study, we conducted a randomized, double-blinded, placebo-controlled clinical trial to evaluate the efficacy and safety of GQD in the treatment of T2D. Furthermore, we examined the structural alterations of gut microbiota in response to GQD treatment intended to alleviate T2D. Materials and methods Study design The study was a 12-week, randomized, doubleblinded and placebo-controlled clinical trial that included a 2-week washout period. It was approved by the Ethics Committee of Guang’anmen hospital of China Academy of TCM. Participants were recruited by Guang’anmen hospital, Dong Zhimen Hospital affiliated to Beijing TCM University, China–Japan Friendship Hospital or Ji Shui Tan Hospital of Beijing from August 2010 to May 2011. All participants signed informed consent forms before beginning the study. The study was conducted in accordance with the principles of the Declaration of Helsinki. The inclusion and exclusion criteria of the patients’ enrollments can be found in the Supplementary Materials and methods. Using an initial screening, including FBG test and 75-g oral glucose tolerance test, 629 recently diagnosed T2D patients who had not received prior pharmacologic treatment for T2D were recruited into the study. After a 2-week washout period and the review of a series of examinations, 403 patients were excluded for not meeting the inclusion criteria and 2 patients were excluded for other reasons. The remaining 224 patients were randomly assigned to four groups of 56 patients. Each group received one of the following treatments: high (HD), moderate (MD) or low dose (LD) GQD, or the placebo for 12 weeks. Randomization was performed centrally and was concealed and stratified in blocks of eight by the PROC PLAN process using the SAS software (SAS Institute Inc., Cary, NC, USA). After the study was completed, a total of 187 patients were included for the final analysis by the verification of data examination committee (Supplementary Figure 1). Microbiota shift in alleviation of type 2 diabetes J Xu et al 553 The ISME Journal
554 our glian (Rhiz ith bead beati ists the amount oa herb in one unit of GOD Chinese lementary Materials and methods).Each unit ugson two times daily for 12 weeks.All one described previously(Zha g et al,2012a,b).The taxonom hr.the ced process ction analysis and s Shan vels an and of OTUs undance d fo changes n serum insulin,lipids pone and R nigh were perfor d at o 4.8 and 12 weeks.Measur ments sof FBG, Natic ormed ing Microcompute CANOCO for ndows cholesterol and o-density h-density 002. instructions (Braak measurements were per Statisti ce was der the as erent genera by taxonomic signmen sequences Clinical and biochemical measurements oAhemidmentemereomearnmpit Real-time quantitative PCR high-dens ermsnitzith cholesterol an Pset of primers was use was measured by high-perormance iquid chroma LADAMS AIC HA-8160: nd the was de ibed before ody RIA (ADVIA C previ Leverkusen,Germany).ELIS Kits were Cha shen ed Min ing the GA USA)and diluted from 1x 102 t01×10 CA.USA)levels. The ISME Joural
Drug administration The TCM formula in our study was GQD, composed of four herbs, namely: Gegen (Radix Puerariae), Huangqin (Radix Scutellariae), Huanglian (Rhizoma Coptidis) and Gancao (Honey-fried Licorice Root) (Supplementary Figure 2a). Supplementary Table 1 lists the amount of each herb in one unit of GQD formula in each group. Herbs were all provided and quality controlled by Beijing Shuangqiaoyanjing Chinese herb manufacturer (Supplementary Materials and methods). The TCM intervention and placebo were given as decoction; these were prepared by Beijing Jiulong Pharmaceutical Factory according to a standard production process (Supplementary Materials and methods). Each unit of GQD formula or placebo yielded 300 ml of decoction. Each patient orally took 150 ml of the decoction two times daily for 12 weeks. All of the drugs and decoctions were quality controlled throughout the trial, and the placebo decoction was prepared by the same standardized process (Supplementary Materials and methods). Study evaluation and outcomes The following primary efficacy outcomes were used: changes in HbA1c, FBG and 2-h postprandial blood glucose (2h-PBG) levels. Secondary efficacy outcomes included changes in serum insulin, lipids levels and body mass index. Study assessments were performed at 0, 4, 8 and 12 weeks. Measurements of FBG, 2h-PBG, body mass index, waist circumference and hip circumference were taken at 0, 4, 8 and 12 weeks. Serum HbA1c, insulin, total cholesterol, triglycerides, high-density lipoprotein cholesterol and low-density lipoprotein cholesterol measurements were performed at weeks 0 and 12. Fecal samples were collected every 4 weeks until the end of trial for gut microbiota analysis. Clinical and biochemical measurements Biochemical measurements of glucose, serum lipids, HbA1c and insulin were performed in a central laboratory (Guang’anmen Hospital of China Academy of TCM, Beijing, China). Glucose, serum total cholesterol, triglycerides, high-density lipoprotein cholesterol and low-density lipoprotein cholesterol were measured by enzymatic methods (Olympus AU2700; Olympus Co. Ltd., Tokyo, Japan). HbA1c was measured by high-performance liquid chromatography using a variant hemoglobin HbA1c assay (ADAMS A1c HA-8160; Arkray Inc., Kyoto, Japan). Serum insulin was measured using a double-antibody RIA (ADVIA Centaur; Bayer Diagnostics, Leverkusen, Germany). ELISA Kits were used to measure plasma orosomucoid (Assaypro, St Charles, MO, USA), adiponectin (R&D System, Minneapolis, MN, USA), tumor necrosis factor-a (R&D System) and serum amyloid A protein (Invitrogen, Carlsbad, CA, USA) levels. Fecal DNA extraction and pyrosequencing Genomic DNA of each fecal sample was extracted by a InviMag Stool DNA Kit (Invitek, Berlin, Germany) combined with bead beating as previously published (Zhang et al., 2012b). The extracted genomic DNA was used as the template to amplify the V3 region of 16S rRNA genes. PCR reactions, pyrosequencing of the PCR amplicons and quality control of raw data were performed as described previously with minor modification (Zhang et al., 2009; Wang et al., 2011). Bioinformatics and multivariate statistics High-quality sequence alignments were performed using NAST. Sequence clustering by CD-hit and OUT delineation by DOTUR were performed as described previously (Zhang et al., 2012a, b). The representative sequences of operational taxonomy units (OTUs) with their relative abundance were used to calculate rarefaction analysis and Shannon diversity index by QIIME (Caporaso et al., 2010). In addition, the representative sequences were inserted into a pre-established phylogenetic tree of the full-length 16S rRNA gene sequences in ARB (Ludwig et al., 2004). Then, the phylogenetic tree and the relative abundance table of representative sequences of OTUs were used for UniFrac principal coordinate analysis (PCoA) (Lozupone and Knight, 2005). The statistical significance between different groups was assessed by multivariate analysis of variance in MATLAB 2010b (The MathWorks Inc., Natick, MA, USA). Redundancy analysis was performed using CANOCO for Windows 4.5 (Microcomputer Power, Ithaca, NY, USA) according to the manufacturer’s instructions (Braak and Smilauer, 2002). Statistical significance was assessed by MCPP with 499 random permutations under the full model. Ribosomal Database Project Classifier was used to assess the amounts of different genera by taxonomic assignment of all sequences. Real-time quantitative PCR of F. prausnitzii Real-time quantitative PCR (q-PCR) was used to determine the amounts of total bacteria and F. prausnitzii through detection of 16S rRNA genes. A set of universal primers was used to amplify a conserved 16S rDNA sequence in all bacteria as shown before (Wang et al., 2011). A set of specific primers was used to amplify a conserved 16S rDNA sequence in F. prausnitzii and the q-PCR reaction system and the program was described before (Balamurugan et al., 2008). A plasmid containing a F. prausnitzii full-length 16S rDNA from a previous study (Shen et al., 2006) was prepared using the EZNA Plasmid Mini Kit I (OMEGA, Doraville, GA, USA) and diluted from 1 103 to 1 109 (copies ml 1 ) to construct a standard curve for the detection of F. prausnitzii. We selected reactions with efficiencies ranging from 0.90 to 1.05 for Microbiota shift in alleviation of type 2 diabetes J Xu et al 554 The ISME Journal
四 of 555 endorf.Hamburg. cler ep Germ Spearman 到and P.value 50351 .36±0 cousnitz and clinical parameters 如 05 HD groups, Results 品e品 rom was observe in plas Spplmmlarha majo nts are shown in the ous erse events occurred uero noimg o in COD decoction ( bar-coded the structural ange9uor 1) as on sity curve altho ed with the eeks of treatment equencing. most of the diversityhac ients.The HD and MD group 1T2 eicrobiote 1.5 可i 1 GOD signif ontrol and HOMA-8 in T2D pa al Ch in FBG.(b)cha The ISME Journal
further analysis. Standard and quantified samples were performed in triplicate. PCR reactions were performed using iQ SYBR Green SuperMix (BioRad, Richmond, CA, USA) on a MasterCycler ep Realplex 4s (Eppendorf, Hamburg, Germany). Spearman’s correlation coefficient (R) and P-value were used to compare the amounts of F. prausnitzii measured by q-PCR and pyrosequencing. This coefficient was also used to evaluate the relationship between F. prausnitzii and clinical parameters using MATLAB 2010b. Results The major components of GQD decoction There were four major categories of compounds in the GQD decoction. Flavones (baicalin, puerarin, wogonoside, daidzin, liquiritin, baicalein and wogonin), alkaloids (berberine, coptisine, palmatine and jatrorrhizine) and triterpenoid sapnins (glycyrrhizin) were detected in the decoction, among which baicalin, puerarin and berberine were the major components (Supplementary Table 2). The chemical structures of these 12 components are shown in the Supplementary Figure 2c. Carbohydrates (starch, sucrose, reducing sugar and soluble dietary fiber) were also detected. Insoluble dietary fiber was undetectable in GQD decoction (Supplementary Table 3). GQD significantly improved glycemic control in T2D patients In our 12-week, randomized, double-blinded, placebo-controlled clinical trial (Supplementary Figure 1), the data of 187 participants were analyzed as shown in Supplementary Table 4. The baseline variables were not significantly different among the four groups. After 12 weeks of treatment, GQD significantly improved glycemic control in T2D patients. The HD and MD groups, when compared with the placebo and LD groups, showed significant reductions in adjusted mean changes from baseline of FBG ( 1.46±0.23 and 1.09±0.21 vs 0.16±0.22 and 0.24±0.24 mmol l 1 ; Po0.001 for HD vs LD and placebo; Po0.01 for MD vs LD and placebo). Similarly, the HD and MD groups showed significantly reduced HbA1c ( 0.88±0.14 and 0.75±0.13 vs 0.35±0.13 and 0.36±0.15%; Po0.01 for HD vs LD; Po0.05 for HD vs placebo; Po0.05, MD vs LD and placebo) (Figures 1a and b). A decrease in the mean change of 2h-PBG from baseline was also observed in the treated groups, although not reaching significant level. (Supplementary Figure 3). In addition, ANCOVA analysis showed that HOMA-b was significantly improved by HD GQD treatment compared with the placebo and LD groups (Figure 1c). Plasma orosomucoid was significantly reduced by HD GQD treatment (P ¼ 0.023) (Supplementary Figure 4a) and the HD group showed a significant reduction (P ¼ 0.034) in mean change from baseline of plasma orosomucoid compared with the LD group (Supplementary Figure 4b). No significant difference was observed in plasma adiponectin, tumor necrosis factor-a or serum amyloid A among the four groups (Supplementary Figures 5a–c). Finally, no drug-related serious adverse events occurred in this study. Overall structural modulation of gut microbiota after GQD treatment First, we used a bar-coded pyrosequencing run to analyze the structural changes of gut microbiota in the four groups before and after GQD treatment. In total, 483 304 usable raw sequences (34 753 unique sequences) and 3222 OTUs were obtained from 235 samples with an average of 2057±396 per sample. Rarefaction and Shannon diversity curves revealed that, although no rarefaction curves plateaued with the current sequencing, most of the diversity had already been captured (Supplementary Figure 6). Weighted and unweighted UniFrac PCoA analysis revealed that gut microbiota structure of the treated groups showed a dose-dependent deviation Figure 1 GQD significantly improved glycemic control and HOMA-b in T2D patients. (a) Change in FBG, (b) change in HbA1c and (c) change in HOMA-b. Placebo (n ¼ 41), LD (n ¼ 50), MD (n ¼ 52) and HD (n ¼ 44). Data are presented as mean±S.E.M. *Po0.05, **Po0.01 and ***Po0.001 vs placebo using ANCOVA; þ Po0.05, þ þ Po0.01 and þþþ Po0.001 vs LD using ANCOVA. Microbiota shift in alleviation of type 2 diabetes J Xu et al 555 The ISME Journal
556 from the rsity curves most o peTstyghadapeadypbee and PCA show seks 0.4.8 and 12 in Hp m that f its bas ne and c analysis2h-PBGaignificantly reduced.but a 6 ★★★ 00 001 0.01 0.03 -0.05 0.12 0.04 0.12 0■000▣▣▣ Placebo-wko OLD-wko MD-wk0 HD-wko ☐Placebo-wk12 LD-wk12 口MD-wk12 HD-wk12 gue Iterations of the eated with different d GOD at weeks o and 12 on th dinate (PC) 0.05 0.01 0.03 -0.07 -0.1 0.05 0.0 0.10.15 ◆■△●。△◆口 Placebo-wk0 ▲Placebo-wwk4 ◆Placebo-wk8 口Placebo-wk12 HD-wko △HD-wk4 ◆HD-wk8 HDwk12 ted wit alysis of variance (MANOVA).Each ts the score of all patients in a group at o ne point. d the e ror bar repres HD: The ISME Joural
from the baseline structure, with the HD group reaching significant level in multivariate analysis of variance test (Figures 2a and b and Supplementary Figures 7a and b). To monitor the dynamic changes of gut microbiota during GQD treatment, we analyzed the fecal samples collected at weeks 0, 4, 8 and 12 in HD and placebo groups with a second pyrosequencing run. In total, we generated 680 774 usable raw sequences (37 498 unique sequences) and 4251 OTUs from 288 samples with an average of 2364±443 per sample (one sample was excluded in later analysis because only 81 reads were obtained). Rarefaction and Shannon diversity curves revealed that most of the diversity had already been captured (Supplementary Figure 8). UniFrac PCoA and PCA showed that after 4 weeks of treatment, the gut microbiota structure of the HD group had already significantly diverged from that of its baseline and of the placebo group (Figures 3a and b; Supplementary Figure 9). At that same 4-week analysis, 2h-PBG was significantly reduced, but FBG did not reach a significant level in treated groups (Supplementary Figure 10). As the treatment Figure 2 Dose-dependent alterations of the gut microbiota in T2D patients treated with different doses of GQD at weeks 0 and 12. (a) Weighted Unifrac PCoA of gut microbiota based on the OUT data from the first pyrosequencing run. (b) Clustering of gut microbiota based on mahalanobis distances calculated with multivariate analysis of variance (MANOVA). Each point represents the mean principal coordinate (PC) score of all patients in a group at one time point, and the error bar represents the s.e.m. The sample number (n) at week 0: placebo ¼ 30, LD ¼ 28, MD ¼ 32 and HD ¼ 28. The sample number (n) at week 12: placebo ¼ 30, LD ¼ 28, MD ¼ 32 and HD ¼ 28. ***Po0.0001. Figure 3 Trajectory of the gut microbiota in T2D patients treated with HD GQD and placebo at weeks 0, 4, 8 and 12. (a) Unweighted Unifrac PCoA of gut microbiota based on the OUT data from the second pyrosequencing run. (b) Clustering of gut microbiota based on mahalanobis distances calculated with multivariate analysis of variance (MANOVA). Each point represents the mean principal coordinate (PC) score of all patients in a group at one time point, and the error bar represents the s.e.m. Placebo: n ¼ 36; HD: n ¼ 36. ***Po0.0001. Microbiota shift in alleviation of type 2 diabetes J Xu et al 556 The ISME Journal
557 diiohaaaon-baeataatten8en Figure 9).but the FBG PBG continued to faecolthocterim,nCenmiges,ifdhohocteinand Figure 10) ne study (Supplementary eroldes and usnitz and its association owed that F usnitzii is enriched by GQD,17 OTUs shov wed a d negative 0o Q sicanegativon corelae tae calibacterium was substant y enriched after HbAlc;two of wo from Gm oneom the prominent of the from antly compared with baseline( :0.011 for the LD group.p- Odoribacter pes of this bacte ar m Bacteroides (n=2, Pseu t buty P=0.0435.HD Vs MD:P=0.0120.HD Vs LD).While Bacteroide hat our group found Supplementary uppleme ntary Figure reover MCPP.P=0.002 HD Placebo 1 Axs13.6% re 4 in T2D nRe宽o上化bg The ISME Journal
progressed, the gut microbiota made no additional changes (Figures 3a and b and Supplementary Figure 9), but the FBG and 2h-PBG continued to improve until the end of the study (Supplementary Figure 10). Key phylotypes responding to the GQD treatment in T2D patients By using redundancy analysis, we identified 146 key variables that responded to the GQD treatment (Figure 4). Forty-seven OTUs were increased by GQD and 99 were decreased. Among the 47 OTUs enriched by GQD, 17 OTUs showed a significant negative correlation with FBG; the OTUs belong to Faecalibacterium (n ¼ 4), Gemmiger (n ¼ 4), Bifidobacterium (n ¼ 3), Lachnospiracea_incertae_ sedis (n ¼ 2) and Escherichia (n ¼ 1). Nine OTUs showed a significant negative correlation with HbA1c; two of which were from Faecalibacterium, two from Gemmiger, one from Bifidobacterium, one from Parasutterella and one from Escherichia (Figure 5 and Supplementary Tables 5 and 6). Among the 99 OTUs decreased by GQD treatment, 22 OTUs showed a significant positive correlation with FBG: the OTUs belong to Alistipes (n ¼ 6), Odoribacter (n ¼ 2), Parabacteroides (n ¼ 2), Bacteroides (n ¼ 2), Pseudobutyrivibrio (n ¼ 2) and one OTU to each of the following genera: Butyricimonas, Barnesiella, Oscillibacter and ClostridiumXlVa. Fourteen OTUs showed a significant positive correlation with HbA1c, four of which were from Alistipes, two from Odoribacter, one from each of the following genera: Parabacteroides, Bacteroides, Pseudobutyrivibrio, Butyricimonas, Barnesiella, Oscillibacter and ClostridiumXlVa (Figure 5 and Supplementary Tables 5 and 6). Additionally, taxon-based analysis at the genus level showed that the relative abundance of Faecalibacterium, Gemmiger, Bifidobacterium and Lachnospiracea_incertae_sedis was significantly higher after HD GQD treatment, whereas Alistipes, Parabacteroides and Pseudobutyrivibrio were significantly decreased (relative abundance 41% and Po0.01) (Supplementary Figure 11). Quantification of F. prausnitzii and its association with glycemic parameters One previous study showed that F. prausnitzii is more abundant in the gut of healthy people compared with T2D patients (Qin et al., 2012). In our study, pyrosequencing results indicated that Faecalibacterium was substantially enriched after GQD treatment. To confirm the pyrosequencing results, we quantified the relative abundance of F. prausnitzii, the prominent species of the Faecalibacterium genus, by q-PCR. All three doses of GQD treatment significantly enriched F. prausnitzii compared with baseline (P ¼ 0.004 for the HD group, P ¼ 0.024 for the MD group, P ¼ 0.011 for the LD group, P ¼ 0.217 for the placebo group) (Figure 6a). Particularly, the relative abundance of this bacterium was significantly higher in HD group than in the other three groups after 12 weeks of treatment (P ¼ 0.0002, HD vs placebo; P ¼ 0.0435, HD vs MD; P ¼ 0.0120, HD vs LD), while no significant difference was observed at baseline (Figure 6b). The relative abundance of Faecalibacterium by q-PCR and pyrosequencing showed a high significant correlation with each other, indicating that the differences of Faecalibacterium among four groups found by pyrosequencing is reliable (Supplementary Figure 12). Moreover, the relative Figure 4 Biplot of redundancy analysis (RDA) of the microbiota composition responding to GQD therapy in T2D patients. The type of treatment and time were used as environmental variables. Responding OTUs that explained more than 4% of the variability of the samples are indicated by blue arrows. Top-left, P-value was obtained by Monte Carlo permutation procedure (MCPP). Microbiota shift in alleviation of type 2 diabetes J Xu et al 557 The ISME Journal
558 Genus The ISME Joural
Microbiota shift in alleviation of type 2 diabetes J Xu et al 558 The ISME Journal
559 Discussion of a GQD tre ent is also onent of GOD.However.no gu red this p in T2D with the dose-dependent manner of of was also rep the adults nin ealthy ucces incre 出 losedepenieatmdlhationafthetmierahiote he improve T2D in humans. Moreover, the gut microbiota glucose home occurre ore sign T2D the gut GQL et nigh 6 6 4 MD HD Week o Week 12 ii am 28.p<005,p<00 are ropresented as means..m nan's and log- The color of the nan's co ween the OTU GThe family and gonus of th The ISME Journal
abundance of this bacterium negatively correlated with HbA1c, FBG and 2h-PBG, and positively correlated with HOMA-b by Spearman’s correlation coefficient (Supplementary Table 7). Discussion To our knowledge, this study is the first registered clinical trial to evaluate the dose-dependent efficacy and safety of a Chinese herbal decoction following a randomized, double-blinded and placebo-controlled design. In our study, GQD treatment provided clinically meaningful and dose-dependent reductions in FBG and HbA1c compared with placebo. HOMA-b was also significantly improved by HD GQD treatment. Compared with those in T2D patients treated with berberine in a previously published clinical trial, the HD group in our study showed similar reductions of FBG and HbA1c; however, we did not observe an improvement of dyslipidemia, which had been significantly ameliorated by berberine in the other study (Zhang et al., 2008). In addition, we found a slight alleviation of inflammation, indicated by the reduction of plasma orosomucoid in the HD group. This finding was consistent with a previous nuclear magnetic resonance-based plasma metabonomic study in T2D rats treated with GQD (Tian et al., 2013). These data show that the Chinese herbal formula is effective for glycemic control, making it a promising candidate for diabetes management. Concomitant with the improved glucose homeostasis, we observed an altered microbial composition induced by GQD. In other studies, diet (Cani et al., 2008; Turnbaugh et al., 2008; Xiao et al., 2013), drugs (Zhang et al., 2012b; Shin et al., 2013), prebiotics (Cani et al., 2007b, 2009) and probiotics (Park et al., 2013) have been used to modulate the gut microbiota in metabolic diseases, for example, the increase of Akkermansia spp. by metformin in the improvement of glucose homeostasis in high-fatdiet-induced obese mice (Shin et al., 2013). Berberine, the major pharmacologic component of a Chinese herb C. chinensis (Huang-Lian), has been shown to enrich short-chain fatty acid producers in parallel with the prevention of obesity and insulin resistance in rats (Zhang et al., 2012b). C. chinensis is also a major component of GQD. However, no gut microbiota modulation by berberine or Chinese herbs has been reported in humans. To our knowledge, this study is the first direct evidence in humans that TCMs can modulate the structure of the gut microbiota. Consistent with the dose-dependent manner of T2D amelioration, GQD also exerted a dose-dependent modulation on the gut microbiota, suggesting a strong association between the modulation of gut microbiota and T2D alleviation. Dose-dependent modulation of gut microbiota was also reported in prebiotic studies. For example, an stepwise increase of Bifidobacteria was observed by q-PCR when nine healthy adults received successively increased dosage levels of galactooligosaccharides (Davis et al., 2010). Our study is the first example that dose-dependent modulation of the gut microbiota was concomitant with a dose-dependent alleviation of T2D in humans. Moreover, the gut microbiota alterations occurred before significant improvement of T2D symptoms, supporting the hypothesis that changes in the gut microbiota induced by GQD might contribute to improved glucose homeostasis Figure 6 Relative abundance of F. prausnitzii as quantified by q-PCR. (a) The impact of different treatments on the relative abundance of F. prausnitzii before and after treatment. The sample number (n) for each group: placebo ¼ 21, LD ¼ 19, MD ¼ 19 and HD ¼ 20. (b) The differences of F. prausnitzii among four groups at week 12. The sample number (n) for each group at week 0: placebo ¼ 31, LD ¼ 23, MD ¼ 26 and HD ¼ 29. The sample number for each group at week 12: placebo ¼ 26, LD ¼ 25, MD ¼ 24 and HD ¼ 28. D Po0.05, DDPo0.05 vs before treatment using Wilcoxon’s signed rank test. *Po0.05 and **Po0.01. Placebo, LD, or MD vs HD using Mann–Whitney test. Data are represented as mean±s.e.m. Figure 5 Heat map of redundancy analysis (RDA)-identified key OTUs responding to GQD treatment and Spearman’s correlation between identified OTUs and HbA1c or FBG. The color of the spots in the left panel represents the mean relative abundance (normalized and log-transformed) of the OTU in each group. The OTUs are organized according to their phylogenetic positions. The color of the spots in the right panel represents R-value of Spearman’s correlation between the OTU and HbA1c or FBG. The family and genus names of the OTUs are shown on the right. þ Po0.05. Microbiota shift in alleviation of type 2 diabetes J Xu et al 559 The ISME Journal
560 in T2D onsequence after TCMs whic dulat thr i Baica whi onstrated n et al. 002 in diabetic rats (Waisundara eta enriched by 20O2l.howmldiosedependontynegiceectswi staphylococcus aureus and other B-lactam-tes et al. tion with nla lower abundance of F et a bloods t on gut microbiota d yo a functionally mportant phylo in ntered 2008).Additionally beer oiot-du correlated with HOMA-B 2012h mplying that tsin GOD tha d th The bohydrate compo microbiota Wa: very ifidobacteria dosage differences across the four deco ions.Othe starch easily ddit n.se These results indicate that mod of gut of "pea0d ,migh be i homeos asis bv r 。。 (Saulnier et al..2011).One ex mple from this genu microbiota-targeted diabetes 883ctgeteaeP8aicieRanwtioeal 2000 aula GQD grid BALB/c mic bacte ngss that dtsen et al., 2012.hes res um spp.in the gut.A gh t is stil of 2500 in China Conflict of Interest ties and unknown mechan The authors declare no conflict of interest. nents from being identified.Our study suggests that Acknowledgements avenue for identifying chemical components in The ISME Journal
in T2D patients rather than a mere consequence after the symptoms have been alleviated. Several putative beneficial genera that responded to GQD treatment were identified using redundancy analysis, including Faecalibacterium, Bifidobacterium and Gemmiger. Faecalibacterium, a butyrateproducing bacterial group (Duncan et al., 2002), which was profoundly enriched by GQD in our study, shows anti-inflammatory effects partly through reducing colonic cytokine synthesis and increasing anti-inflammatory cytokine secretion (Sokol et al., 2008). Diabetic obese patients have lower abundance of F. prausnitzii compared with non-diabetic obese patients, which demonstrates a negative correlation with inflammatory cytokines C-reactive protein and interleukin-6 (Furet et al., 2010). Another study also showed that the gut of T2D patients is characterized by a reduction of F. prausnitzii compared with that of healthy people (Qin et al., 2012). Furthermore, F. prausnitzii has been shown to be a functionally important phylotype because it was associated with eight urinary metabolites (Li et al., 2008). In our study, the relative abundance of F. prausnitzii was negatively correlated with HbA1c, FBG and 2h-PBG, and positively correlated with HOMA-b, suggesting that F. prausnitzii might be a pivotal phylotype associated with the improvement of T2D. Two other genera, Bifidobacterium and Gemmiger, which are also reported to confer beneficial effects (Gossling and Moore, 1975; Sokol et al., 2008; Fukuda et al., 2011), were significantly enriched by GQD in T2D patients. A selective increase of Bifidobacteria induced by a prebiotic (namely oligofructose) improved gut permeability and inflammation in ob/ob mice (Cani et al., 2009). In addition, several genera were significantly inhibited by GQD, such as Alistipes and Odoribacter. Higher abundance of several taxa of the genus Alistipes were associated with greater frequency of abdominal pain in irritable bowel syndrome pediatric patients (Saulnier et al., 2011). One example from this genus, Alistipes putredinis, was isolated from inflamed and non-inflamed intestinal tissues of children with suspected acute appendicitis (Rautio et al., 2000, 2003). Alistipes and Odoribacter were significantly increased in grid floor stress-induced BALB/c mice (Bangsgaard Bendtsen et al., 2012). These results suggest that the enrichment of beneficial bacteria, particularly Faecalibacterium spp., and reduction of pathogen-like bacteria might be involved in the amelioration of T2D by GQD. TCM formula, a form of polypharmacy, has been developed and advocated for use in the treatment of many diseases for over 2500 years in China. However, the complexities and unknown mechanisms of TCMs prevent the active chemical components from being identified. Our study suggests that gut microbiota might be involved in the effect of a widely used TCM formula, GQD. This opens an avenue for identifying chemical components in TCMs, which can modulate gut microbiota structure as a potential mechanism for disease alleviation. Baicalin, puerarin and berberine were identified as the three most abundant chemical components in GQD. Baicalin, which demonstrated antioxidant properties together with the glucose-lowering effects in STZ-induced diabetic rats (Waisundara et al., 2009), showed dose-dependent synergic effects with b-lactam antibiotics against methicillin-resistant Staphylococcus aureus and other b-lactam-resistant strains of S. aureus in vitro (Liu et al., 2000). Puerarin can dose-dependently increase glucose utilization in STZ-induced diabetic rats (Hsu et al., 2003). However, it showed poor absorption into the bloodstream after oral administration in vivo (Luo et al., 2011) and was mainly excreted via feces in its intact form (Zhu et al., 1979). It is possible that puerarin could have an impact on gut microbiota, but no study about this has been reported yet. Berberine showed significant glucose-lowering effects in a multicentered, randomized, doubleblinded and placebo-controlled clinical trial (Zhang et al., 2008). Additionally, berberine has been reported to modulate gut microbiota during prevention of high-fat-diet-induced obesity and insulin resistance in rats (Zhang et al., 2012b), implying that berberine could be one of the major active ingredients in GQD that modulated the gut microbiota in our study. The carbohydrate components of GQD likely exerted few effects on gut microbiota because insoluble dietary fiber was undetectable and soluble dietary fiber was very low. Additionally, these components showed no dosage differences across the four decoctions. Other carbohydrates, such as starch, were easily digested and absorbed before they reached large intestine. These results indicate that modulation of gut microbiota by chemical components, such as berberine, might be involved in improving glucose homeostasis by GQD, suggesting that TCMs may serve as a new source for drug leads in gut microbiota-targeted diabetes management. In conclusion, our study suggests that structural alterations of gut microbiota, induced by Chinese herbal formula GQD, are associated with the antidiabetic effects of GQD. In particular, this treatment enriched the number of beneficial bacteria, such as Faecalibacterium spp. in the gut. Although it is still unclear whether changes of gut microbiota by GQD directly contribute to the improvement of glucose homeostasis, our clinical study provides circumstantial evidence that gut microbiota might be involved. Conflict of Interest The authors declare no conflict of interest. Acknowledgements This work was supported by Project 2010CB530601 of the National Basic Research Program of China (973 Program), Microbiota shift in alleviation of type 2 diabetes J Xu et al 560 The ISME Journal
561 28529 References Alberti KG.Zin 2002. d its com n of Ama dherence and trans mensal bacteri K.Oshima K.Nakanishi Y ing wane y.Koh Cy Poitou C.Basdevant A 2008 dies in a southern Indian village.Am ICin: 643 IM D WC.Su HC.Cheng 2003 p.f200 CANOCO Refe w for Khir Mau te watery 5)Mic 291:160 omput V Yo m MA M Knauf C neiee9g85ie sity and insulin e.Diabetes 56 Cani PD,N nck AM,Fava F.Knauf C,Burcelin RG ng M.Rantalainen M.Wang S.Zhou H ove high-fa ciated Liu.Richards RM.(2000)ici and other beta 57:1470148 ane Richter I hrough a nism 110 Ya nt of gu Luo CE.Yuan M.Chen MS,Liu SM,Zhu L Huang BY et a community sequencing The ISME Journa
Project 81273720 of the National Natural Science Foundation of China (NSFC) and Key program 31330005 of the NSFC. References Alberti KG, Zimmet PZ. (1998). Definition, diagnosis and classification of diabetes mellitus and its complications. Part 1: diagnosis and classification of diabetes mellitus provisional report of a WHO consultation. Diabet Med 15: 539–553. Amar J, Chabo C, Waget A, Klopp P, Vachoux C, Bermu´dez-Humara´n LG et al. (2011). Intestinal mucosal adherence and translocation of commensal bacteria at the early onset of type 2 diabetes: molecular mechanisms and probiotic treatment. EMBO Mol Med 3: 559–572. Backhed F, Ding H, Wang T, Hooper LV, Koh GY, Nagy A et al. (2004). The gut microbiota as an environmental factor that regulates fat storage. Proc Natl Acad Sci USA 101: 15718–15723. Balamurugan R, Janardhan HP, George S, Chittaranjan SP, Ramakrishna BS. (2008). Bacterial succession in the colon during childhood and adolescence: molecular studies in a southern Indian village. Am J Clin Nutr 88: 1643–1647. Bangsgaard Bendtsen KM, Krych L, Sorensen DB, Pang W, Nielsen DS, Josefsen K et al. (2012). Gut microbiota composition is correlated to grid floor induced stress and behavior in the BALB/c mouse. PLoS One 7: e46231. Braak CT, Smilauer P. (2002). CANOCO Reference Manual and CanoDraw for Windows User’s Guide: Software for Canonical Community Ordination (version 4.5). Microcomputer Power: Ithaca, New York. Butler AE, Janson J, Bonner-Weir S, Ritzel R, Rizza RA, Butler PC. (2003). Beta-cell deficit and increased beta-cell apoptosis in humans with type 2 diabetes. Diabetes 52: 102–110. Cani PD, Amar J, Iglesias MA, Poggi M, Knauf C, Bastelica D et al. (2007a). Metabolic endotoxemia initiates obesity and insulin resistance. Diabetes 56: 1761–1772. Cani PD, Neyrinck AM, Fava F, Knauf C, Burcelin RG, Tuohy KM et al. (2007b). Selective increases of bifidobacteria in gut microflora improve high-fatdiet-induced diabetes in mice through a mechanism associated with endotoxaemia. Diabetologia 50: 2374–2383. Cani PD, Bibiloni R, Knauf C, Waget A, Neyrinck AM, Delzenne NM et al. (2008). Changes in gut microbiota control metabolic endotoxemia-induced inflammation in high-fat diet-induced obesity and diabetes in mice. Diabetes 57: 1470–1481. Cani PD, Possemiers S, Van de Wiele T, Guiot Y, Everard A, Rottier O et al. (2009). Changes in gut microbiota control inflammation in obese mice through a mechanism involving GLP-2-driven improvement of gut permeability. Gut 58: 1091–1103. Caporaso JG, Kuczynski J, Stombaugh J, Bittinger K, Bushman FD, Costello EK et al. (2010). QIIME allows analysis of high-throughput community sequencing data. Nat Methods 7: 335–336. Collins SM, Kassam Z, Bercik P. (2013). The adoptive transfer of behavioral phenotype via the intestinal microbiota: experimental evidence and clinical implications. Curr Opin Microbiol 16: 240–245. Davis LM, Martinez I, Walter J, Hutkins R. (2010). A dose dependent impact of prebiotic galactooligosaccharides on the intestinal microbiota of healthy adults. Int J Food Microbiol 144: 285–292. Duncan SH, Hold GL, Harmsen HJ, Stewart CS, Flint HJ. (2002). Growth requirements and fermentation products of Fusobacterium prausnitzii, and a proposal to reclassify it as Faecalibacterium prausnitzii gen. nov., comb. nov. Int J Syst Evol Microbiol 52: 2141–2146. Fei N, Zhao L. (2013). An opportunistic pathogen isolated from the gut of an obese human causes obesity in germfree mice. ISME J 7: 880–884. Fukuda S, Toh H, Hase K, Oshima K, Nakanishi Y, Yoshimura K et al. (2011). Bifidobacteria can protect from enteropathogenic infection through production of acetate. Nature 469: 543–547. Furet JP, Kong LC, Tap J, Poitou C, Basdevant A, Bouillot JL et al. (2010). Differential adaptation of human gut microbiota to bariatric surgery-induced weight loss: links with metabolic and low-grade inflammation markers. Diabetes 59: 3049–3057. Gossling J, Moore WEC. (1975). Gemmiger formicilis, n. gen., n. sp., an anaerobic budding bacterium from intestines. Int J Syst Bacteriol 25: 202–207. Hotamisligil GS. (2006). Inflammation and metabolic disorders. Nature 444: 860–867. Hsu FL, Liu IM, Kuo DH, Chen WC, Su HC, Cheng JT. (2003). Antihyperglycemic effect of puerarin in streptozotocin-induced diabetic rats. J Nat Prod 66: 788–792. Khin Maung U, Myo K, Nyunt Nyunt W, Aye K, Tin U. (1985). Clinical trial of berberine in acute watery diarrhoea. BMJ (Clin Res Ed) 291: 1601–1605. Larsen N, Vogensen FK, van den Berg FWJ, Nielsen DS, Andreasen AS, Pedersen BK et al. (2010). Gut microbiota in human adults with type 2 diabetes differs from non-diabetic adults. PLoS One 5: e9085. Le Chatelier E, Nielsen T, Qin J, Prifti E, Hildebrand F, Falony G et al. (2013). Richness of human gut microbiome correlates with metabolic markers. Nature 500: 541–546. Li M, Wang B, Zhang M, Rantalainen M, Wang S, Zhou H et al. (2008). Symbiotic gut microbes modulate human metabolic phenotypes. Proc Natl Acad Sci USA 105: 2117–2122. Liu IX, Durham DG, Richards RM. (2000). Baicalin synergy with beta-lactam antibiotics against methicillinresistant Staphylococcus aureus and other betalactam-resistant strains of S. aureus. J Pharm Pharmacol 52: 361–366. Lozupone C, Knight R. (2005). UniFrac: a new phylogenetic method for comparing microbial communities. Appl Environ Microbiol 71: 8228–8235. Ludwig W, Strunk O, Westram R, Richter L, Meier H, Yadhukumar et al. (2004). ARB: a software environment for sequence data. Nucleic Acids Res 32: 1363–1371. Luo CF, Yuan M, Chen MS, Liu SM, Zhu L, Huang BYet al. (2011). Pharmacokinetics, tissue distribution and relative bioavailability of puerarin solid lipid nanoparticles following oral administration. Int J Pharm 410: 138–144. Microbiota shift in alleviation of type 2 diabetes J Xu et al 561 The ISME Journal