魏建敏，杨华连，陈莉*，卢红梅，石庆叠，张祥瑞，涂青

1(贵州省发酵工程与生物制药重点实验室，贵州 贵阳，550025)2(贵州大学 酿酒与食品工程学院，贵州 贵阳，550025)

糖尿病是一种由胰岛素分泌不足、糖脂代谢失衡而导致的以高血糖为特征，且具有遗传倾向的内分泌代谢性疾病[1-2]，其中以2型糖尿病(type 2 diabetes mellitus，T2DM)居多，占糖尿病病例的90%以上[3-4]，该病还会引发多种并发症，如心血管疾病、糖尿病、肾病、冠心病等[5]。目前，对T2DM的治疗通常采用口服或注射胰岛素等药物治疗的方式，但其存在严重的副作用[6-7]。因此，从天然植物中寻找具有显著疗效和低副作用的活性成分已成为当前糖尿病治疗领域的研究热点。

果桑是桑科桑属以结果为主、果叶兼用的落叶木本植物桑树的统称[8]，是一种“药食同源”植物。其叶、枝、果、根中均含有多种营养成分及黄酮、多糖、生物碱类等已被证实具有降糖功效的生物活性物质[9-11]。研究表明，果桑中存在天然的α-葡萄糖苷酶、α，β淀粉酶抑制剂1-脱氧野尻霉素(1-deoxynojirimycin，DNJ)和荞麦碱，桑叶DNJ可通过激活db/db小鼠骨骼肌中的胰岛素信号PI3K/AKT途径来提高胰岛素敏感性，从而达到降血糖的效果[12]。

肠道菌群是人体最大的微生态系统，人体肠道内居住着超过1014个的微生物[13]，它们在调节机体能量摄取、肠道通透性、维持肠道免疫屏障等方面发挥着重要作用[14]。近年来，很多研究表明T2DM与肠道菌群失调有关[15-17]。LARSEN等[18]研究T2DM和非T2DM人体肠道菌群组成的差异发现，与非糖尿病人相比，糖尿病人的厚壁菌门(Firmicutes)和梭状体类的比例明显降低；SEDIGHI等[19]研究发现，在T2DM患者中，乳酸菌水平明显较高，而双歧杆菌在健康人中明显较多。何雪冬等[20]通过Meta分析发现，与健康人群相比，糖尿病患者乳杆菌、拟杆菌、梭菌数量有不同程度的增加，双歧杆菌数量下降。基于之前的研究成果，本研究以桑叶、桑枝皮、桑根皮为原料，按62∶19∶19的质量比制成袋泡茶，以链脲佐菌素(streptozocin, STZ)诱导联合高脂饮食饲养建立2型糖尿病模型小鼠为对象，通过高通量测序方法，从小鼠肠道菌群方面探究果桑茶对T2DM模型小鼠肠道菌群的调整作用，为果桑降糖保健产品的开发提供理论基础。

1 材料与方法

1.1 材料与试剂

体重(20±2) g清洁级5周龄雄性C57BL/6J小鼠购于贵州医科大学实验动物中心，使用许可证号SCXX(京)2014-0004；桑叶、桑枝皮、桑根皮采自贵州开阳某桑葚基地，按62∶19∶19的质量比制成袋泡茶；STZ(纯度≥98%)，Sigma公司；羧甲基纤维素钠(carboxymethylcellulose sodium，CMC-Na)(纯度≥98%)，阿拉丁化学试剂有限公司；0.1 mol/L柠檬酸钠缓冲液(pH 4.5)(分析纯)，索莱宝生物试剂有限公司。

1.2 仪器与设备

FA2004N分析电子天平，上海菁海仪器有限公司；FW-80万能粉碎机，上海仕元科学器材有限公司；VaCo 5-Ⅱ-D真空冷冻干燥机，Zirbus technology GmbH Hilfe Gottes；SW-CJ-1FD无菌操作台，上海苏净实业有限公司；鱼跃580型血糖仪，江苏鱼跃医疗设备有限公司。

1.3 实验方法

1.3.1 试剂制备

果桑复配袋泡茶浸提液冻干粉：将袋泡茶在水温75 ℃，茶水比1∶100(g∶mL)，冲泡时间50 min的条件下冲泡1次，取其浸提液置于真空冷冻干燥机中冻干备用。

10%果桑复配茶高脂饲料：将果桑复配茶添加10%(质量分数)到高脂饮食饲料中，并压制成复合饲料备用。

STZ溶液：现用现配，临用时使用0.1 mol/L、pH 4.5柠檬酸缓冲液配制，按照给药剂量35 mg/kg配制成质量浓度为3.5 mg/mL的STZ溶液。

0.5%(质量分数)CMC-Na：称取CMC-Na 0.5 g，加入双蒸馏水至100 mL，并加热溶解至澄清溶液，置于4 ℃冰箱中保存备用。

1.3.2 动物分组及处理

小鼠进入实验室稳定喂养7 d后禁食不禁水过夜12 h，次日通过剪尾取血法测定小鼠空腹血糖和体重，剔除极值后随机选取10只小鼠作为正常对照组给予普通饲料喂养，其余小鼠采用高脂饮食饲养16周致小鼠呈肥胖状态，按小鼠体重注射35 mg/kg STZ溶液，给药量为0.1 mL/10g，7 d后血糖水平高于11.1 mmoL/L即为建模成功。

将血糖水平高于11.1 mmoL/L的小鼠稳定性饲养4周后，按照空腹血糖水平随机分组为4组，每组10只，分别为正常对照组(Control)、模型对照组(Model)、阳性对照组(Metformin)及果桑茶试药组GD-2组、GD-3组，共5组。GD-2受试组采用小鼠自由摄食的方式，正常对照组、模型对照组及GD-3受试组采取灌胃方式进行给药，给药剂量依据《药理实验方法学》中“标准体重动物剂量换算表”来确定，水提物灌胃剂量为200 mg/kg，具体分组及给药情况见表1。灌胃开始后，根据小鼠体重变化调整小鼠灌胃剂量，持续28 d。

表1 试验小鼠分组及给药情况Table 1 Grouping and administration of experimental mice

1.3.3 样品采集

连续给药28 d后禁食12 h，每组随机选3只小鼠收集粪便并于液氮中贮存。

1.3.4 DNA提取及测序

样本DNA提取并进行高可变区PCR扩增后，使用Nova Seq PE250 测序平台对各组小鼠肠道内容物16S rDNA V3+V4区进行测序。细菌16S rDNA V3+V4区的通用引物名称为341F(5′-CCTACGGGNGGCWGCAG-3′,)、805R(5′-GACTACHVGGGTATCTAATCC-3′)。测序数据使用最新的性能更优的QIIME 2分析流程，调用DADA2[21]对数据进行去噪，相当于以100%的相似度聚类，去冗余后得到特征(Feature)，然后基于Feature数据进行菌群结构组成及菌群多样性分析。

1.4 数据分析

将高通量测序得到的原始数据进行优化，得到有效数据，基于有效数据按97%相似性进行操作分类单元(operational taxonomic unit, OTU)聚类分析，对OTU序列进行物种注释。基于OTU结果，通过α多样性分析小鼠肠道菌群多样性、β多样性分析小鼠肠道菌群间的组成差异，根据绘制的样品各分类水平下的群落结构图、热图分析菌群结构组成，根据线性判别分析(linear discriminant analysis，LDA)值分布柱状图和进化分支图进行LEfSe(LDA effect size)多级物种差异判别分析，采用SPSS 19.0进行显著性差异分析，P<0.05认为差异显著。

2 结果与分析

2.1 果桑茶对T2DM模型小鼠空腹血糖值影响

如表2所示，当连续给药后，与模型组相比，其他各组小鼠的空腹血糖值明显降低，差异极显著(P<0.01)，阳性组的降糖幅度为43.60%，GDm组的降糖幅度为40.05%，其降糖作用与阳性组相比无明显差异；GDh组其降幅为25.83%。以上结果说明，果桑茶具有和二甲双胍药物一样的降血糖功效，其中GDm组降糖效果最明显。

表2 果桑茶对T2DM模型小鼠空腹血糖值影响Table 2 Effect of fruit mulberry tea on fasting blood glucose in T2DM model mice

2.2 菌群多样性分析

2.2.1 α多样性分析

本研究采用Chao1、Shannon及Simpson指数来评估小鼠肠道菌群的丰富度和多样性，Good′s Coverage反映样本的测序深度。由表3可知，5个组的Good′s Coverage值均为1，说明测序深度已饱和，实验数据合理可用。Chao1指数反映菌群的总数，Chao1值越大，菌群数越多，所测定样品中，果桑茶给药组(GDh组、GDm组)菌群的多样性均高于模型对照组和阳性对照组，但低于正常对照组。Shannon和Simpson指数反映菌群的多样性，Shannon值越大，菌群多样性越高，5组样品中正常对照组的Shannon值最大，模型对照组最小，其余3个组差别不大，均小于正常对照组，大于模型对照组；Simpson值越大，菌群多样性越低，GDh组的Simpson值最小，阳性对照组次之，其余3个组的值为0.96～0.98。以上结果表明，长期持续的高脂高糖饮食会很大程度地降低小鼠肠道菌群的多样性，而果桑茶产品对调节小鼠肠道菌群多样性有显著效果，这可能是果桑茶中的降糖成分发挥了积极作用。

表3 α多样性比较结果Table 3 α diversity comparison results

2.2.2 β多样性分析

β多样性分析可反映不同样本菌群间是否存在差异性。本试验采用主成分分析(principal component analysis，PCA)分析样本间的多样性。

如图1所示，横坐标为PCA1，贡献率为40.44%，纵坐标为PCA2，贡献率为14.04%，可以代表大部分变量信息，是样品菌群结构组成差异的主要因子。分析可知，果桑茶给药组、模型对照组和阳性对照组肠道菌群组成结构较正常对照组发生了偏移，模型对照组和阳性对照组均与GDh组在物种组成上有一定程度的相似，而GDm组与正常对照组最接近，说明GDm组在调节小鼠肠道菌群方面效果较为显著。

图1 样本细菌主成分分析Fig.1 Principal component analysis of bacteria in samples

2.3 菌群构成分析

2.3.1 菌群组成分析

各组小鼠肠道中相对丰度较高的菌群比较如表4所示。由图2、图3可知，在门水平上5个样品肠道菌群以厚壁菌门、拟杆菌门(Bacteroidetes)、变形菌门(Proteobacteria)等为主。其中，拟杆菌门是正常对照组的绝对优势菌门，其相对丰度为47.89%，与其他组差异显著(P<0.05)，模型对照组中拟杆菌门的相对丰度为15.86%，比正常对照组少32.03%，GDm组拟杆菌门的相对丰度为21.06%；厚壁菌门为正常对照组的次要优势菌门，相对丰度为39.10%，在模型对照组中其相对丰度达到了54.49%，成为绝对优势菌门，经果桑茶给药喂养后，给药组厚壁菌门数量有不同幅度降低，逐渐接近正常对照组水平，其中GDm组效果最明显，与正常对照组相比差异不显著(P<0.05)；变形菌门在正常对照组中的相对丰度为4.72%，而在模型对照组中增加至20.35%，经药物治疗后，该菌群涨幅不大，与模型对照组无显著差异(P<0.05)，说明果桑茶对此菌群的调节作用不明显。

属水平上的物种相对丰度如图4所示。Muribaculaceae_unclassified是一种有益菌属[22]，由图4、图5可知，其是正常对照组的优势菌属，相对丰度为28.89%，在模型对照组、阳性对照组、GDm组、GDh组的相对丰度分别为2.79%、0.20%、6.63%和4.25%，说明果桑茶能在一定程度上抑制有益菌属Muribaculaceae_unclassified的下降。嗜胆菌属(Bilophila)为模型对照组的优势菌属，其相对丰度为11.23%，而在正常对照组中仅占0.06%，阳性对照组、果桑茶给药组均有效地抑制了这种现象，其中以阳性对照组的效果最为显著，嗜胆菌属的相对丰度为1.56%，其次是GDh组，嗜胆菌属的相对丰度为7.86%。阳性对照组、GDm组和GDh组的优势菌属分别为Lachnospiraceae_unclassified、艾克曼菌属(Akkermansia)和毛螺菌科NK4A136群(Lachnospiraceae_NK4A136_group)。

以上结果表明T2DM模型组小鼠与正常组小鼠肠道微生物群落结构存在明显差异，果桑茶可增加糖尿病小鼠中有益菌群的相对丰度，具有调节糖尿病小鼠肠道菌群失调的效果[23]。

表4 各组小鼠丰度较高的微生物比较Table 4 Microbe comparison of mice with higher abundance in each group

图2 门水平上的物种相对丰度Fig.2 Relative abundance of species at the phylum level

图3 门水平上有显著性差异的主要物种Fig.3 The main species with significant difference at the phylum level注：图中不同字母表示不同样品间物种差异显著(P<0.05)(下同)

图4 属水平上的物种相对丰度Fig.4 Relative abundance of species at the genus level

图5 属水平上有显著性差异的主要物种Fig.5 The main species with significant difference at the genus level

2.3.2 菌群结构差异性分析

将属水平丰度排名在前30的菌群，以不同组间丰度相似性聚类的方式绘制热图，根据颜色的深浅定义数值的大小，通过颜色的不同来反映各组小鼠肠道菌群结构组成的相似性及差异性。由图6分析可知，GDh组和模型对照组归为一类，GDh组中毛螺菌科NK4A136群(Lachnospiraceae_NK4A136_group)相对丰度极高，相关研究表明人体肠道共生菌群中的毛螺菌科成员可表达2种“超级抗原”，这类超级抗原可激活体内的免疫球蛋白A(immunoglobulin A，IgA)应答，在肠道稳态中发挥重要作用[24]，说明果桑茶能促进高脂饮食小鼠体内有益菌群的增长。GDm组和正常对照组归为一类，表明GDm组小鼠肠道菌群聚类分析结果与正常小鼠最为接近，GDm组调节小鼠肠道菌群的效果最好。

由以上分析可知，不同组小鼠的肠道菌群具有相似性的同时也存在一定差异性，高脂饮食造模致糖尿病的过程中小鼠肠道菌群逐渐偏离正常，而果桑茶对小鼠肠道菌群的调节具有积极作用。

图6 属水平上的相对丰度聚类热图Fig.6 Cluster heatmap of relative abundance at the genus level

2.4 LEfSe多级物种差异判别分析

LEfSe是一种用于发现和解释高维度数据生物标识的分析工具，能够在组与组之间寻找具有统计学差异的生物标志物 (biomarker)。图7是从门到种水平下各组样本的LDA值分布柱状图，物种的LDA得分用横坐标来表示，LDA值>3的物种为具有统计学差异的biomarker。图8是各组样本的进化分支图。

结合图7、图8分析可知，正常对照组中对群落结构影响较大的物种有拟杆菌纲(Bacteroidia)、拟杆菌门、拟杆菌目(Bacteroidales)、Muribaculaceae、Muribaculaceae-unclassfied、Muribaculum等。阳性对照组中对群落结构影响较大的物种有厚壁菌门、梭杆菌目(Clostridiales)、梭杆菌纲(Clostridia)、Bacteria、Ruminiciostridium-9-unclassfied、Ruminiciostridium-9等。模型对照组中对群落结构影响较大的有嗜胆菌属、unclassfied-Bilophila、理研菌科(Rikenellaceae)、另枝菌

属、Paramuribaculum、Paramuribaculum-intestinale、Alistipes-unclassfied，Alistipes是小鼠中细菌最丰富的属之一，它与调节胰岛素释放和逆转胰岛素抵抗的肠道激素的释放有关[25]。GDm组中对群落结构影响较大的物种有拟杆菌属(Bacteroides)、拟杆菌科(Bacteroidaceae)、unclassfied-Bacteroides-sp、Desulfovibrionaceae-unclassfied、副拟杆菌属(Parabcteroides)、Tannerellaceae、Parabcteroides-unclassfied。GDh组中对群落结构影响较大的物种有毛螺菌科(Lachnospiraceae)、瘤胃菌属UCG014(Ruminococcaceae-UGG-014)、Ruminococcaceae-UGG-014-unclassfied、消化链球菌科(Peptostreptococcaceae)、Romboutsia。

图8 各组样本的进化分支图Fig.8 Evolutionary branch diagram of each group of samples注：小圆圈:图中由内至外辐射的圆圈代表了由门至属的分类级别。 不同分类级别上的每一个小圆圈代表该水平下的一个分类，小圆圈 直径大小与相对丰度大小呈正比；颜色：无显著差异的物种统一着色 为黄色，差异显著的物种Biomarker跟随组别进行着色，红色节点表示 在红色组别中起到重要作用的微生物类群，绿色节点表示在绿色组别中 起到重要作用的微生物类群；未能在图中显示的Biomarker对应的物种 名会展示在右侧，字母编号与图中对应

3 讨论

肠道菌群的状态可以通过粪便菌群的变化反映出来，高通量测序技术现已广泛应用于肠道菌群研究等多个领域，与传统检测方法相比，具有更准确定量、实时检测的优势。本研究通过高通量测序技术对小鼠粪便进行检测，分析探讨果桑茶对2型糖尿病模型小鼠肠道菌群的调节效果。

菌群α多样性分析表明，模型对照组小鼠肠道菌群的多样性与正常对照组相差较大，正常对照组与给药组的Chao1和Shannon值均大于模型对照组，模型对照组的Simpson值大于GDm组和正常组而小于GDh组和阳性对照组，说明糖尿病的产生会降低小鼠肠道菌群的多样性，而果桑茶可以提高肠道菌群的丰度，其中GDm组效果最好。β多样性分析发现正常对照组与GDm组、GDh组与阳性对照组群落构成差异较小，GDm组小鼠的肠道菌群构成最接近正常对照组。

菌群构成分析发现，在门和属水平上，各组小鼠肠道菌群存在明显差异。在门分类水平上，拟杆菌门和厚壁菌门是正常对照组小鼠肠道菌群中的优势菌门，这与BEZIRTZOGLOU等[26]，QIN等[27]的研究结果相似，厚壁菌门和拟杆菌门能共同促进对糖的分解利用[28]，相关研究表明厚壁菌门与拟杆菌门比值与个体体重呈正相关，经高脂饮食后，模型对照组中厚壁菌门相对丰度升高而拟杆菌门相对丰度降低，给药治疗后，拟杆菌门相对丰度有所增加，GDm组效果最显著。在属分类水平上，Muribaculaceae_unclassified是正常组中的优势菌属，其在模型对照组和阳性对照组中的丰度最低，与之相比在果桑茶给药组中其丰度有明显提高。嗜胆菌属是一种2型糖尿病的标志物，其中的沃氏嗜胆菌Bilophilawadsworthia能引起和加重炎症[29]，与正常对照组相比，模型对照组中嗜胆菌属呈明显上升趋势，而给药组有效地抑制了这种现象。GDm组中的优势菌属艾克曼菌属可降解黏蛋白生成短链脂肪酸，与肥胖、糖尿病等疾病呈负相关[30]。GDh组优势菌属毛螺菌科NK4A136群对肠道黏膜损伤具有一定的修复作用，在调节和稳定肠道菌群方面发挥着重要作用[31]。说明果桑茶能在调节2型小鼠糖尿病肠道菌群方面发挥积极作用。

LEfSe多级物种差异判别分析表明，与模型对照组相比，GDm组中的嗜胆菌属等有害菌群的丰度降低，而有益菌群拟杆菌属等的丰度增高；此外，GDm组小鼠的肠道菌群多样性和物种丰富度较模型对照组显著提高，与正常对照组小鼠接近，这些结果表明果桑茶在将糖尿病小鼠肠道菌群紊乱调节到接近正常状态方面具有一定的功效，这可能与果桑茶中的多糖、黄酮等生物活性物质有关[32]。

4 结论

本研究采用高通量测序对果桑茶对2型糖尿病模型小鼠肠道菌群的影响进行研究。结果不仅验证了2型糖尿病会破坏小鼠肠道菌群平衡，降低菌群丰度，同时证明了果桑茶能提高小鼠肠道菌群的丰度，并在一定程度上增加有益菌群的数量，降低高血糖特征菌群及有害菌群的影响，促进肠道菌群良性发展。本研究的结果可为果桑资源的综合应用奠定理论基础，同时为果桑茶应用于2型糖尿病的干预治疗提供一定的科学依据。