肠道菌群在糖尿病合并衰弱的发生发展中起到了重要的作用,因此研究中国中老年糖尿病衰弱人群肠道菌群的结构及功能特征性,具有重要意义。本研究收集30例社区中老年糖尿病患者粪便进行宏基因测序,结合临床数据,分析糖尿病衰弱及衰弱前期与糖尿病非衰弱两组人群菌群结构及功能差异性。结果表明:糖尿病衰弱及衰弱前期人群菌群α多样性有所降低,柯林斯氏菌属(Collinsella)、丁酸弧菌属(Butyricimonas)在衰弱人群中相对丰度较高。糖尿病衰弱及衰弱前期人群富集上调的功能基因分布在组氨酸代谢(Histidine metabolism)、EB病毒感染(Epstein-Barr virus infection)、硫代谢(Sulfur metabolism)和Ⅱ型聚酮生物合成(Biosynthesis of type Ⅱ polyketide products)等,而丁酸代谢(Butanoate metabolism)以及苯丙氨酸代谢(Phenylalanine metabolism)等在该人群中下调。本研究为糖尿病衰弱人群的肠道菌群机制探索奠定了理论基础,同时也为糖尿病衰弱及衰弱前期人群的预防与干预提供了一定的科学依据。
引用本文: 彭旭超, 赵艳莉, 林泰平, 舒哓宇, 侯利莎, 高浪丽, 王慧, 葛宁, 岳冀蓉. 中老年糖尿病衰弱患者与肠道菌群结构及功能相关性研究. 生物医学工程学杂志, 2021, 38(6): 1126-1133. doi: 10.7507/1001-5515.202101083 复制
引言
全球老龄化问题正日益凸显,中国老龄化进程明显快于世界平均水平。2010年我国65岁及以上人群占比8.87%,预计到2050年,这一数字将提高至34%[1]。衰弱是老龄化突出问题之一,是指老年人生理储备下降导致机体易损性增加和抗应激能力减退的非特异性状态,其核心是老年人生理储备下降或异常,外界较小刺激即可引起临床事件的发生[2]。衰弱人群更容易出现各种不良预后,包括跌倒、失能、认识功能下降、再入院率增加甚至死亡风险等[3]。
糖尿病与衰弱两者可相互作用。一方面,高血糖状态会抑制骨骼肌生长,严重时会造成肌肉萎缩[4],而胰岛素抵抗会抑制骨骼肌的能量代谢造成肌肉再生功能障碍[5],机体衰弱发生风险增加[6]。另一方面,衰弱促使胰岛素抵抗增加,血糖呈高水平状态,增加糖尿病发生风险[7]。我国的流行病学调研显示,我国65岁以上糖尿病合并衰弱人群占糖尿病人群比例为5%~48%[8],该人群出现低血糖及死亡风险高[9],预后较单纯糖尿病老年人群差,极大地增加了家庭照护和社会经济负担。
微生物在人体内广泛分布,其中以胃肠道菌群数量最多。既往研究表明,衰弱与肠道菌群α多样性呈负相关[10],衰弱程度与Ruminococcacea菌、F. prausnitzii菌及Lachnospiraceae菌丰度呈负相关,与Eubacteriumdolichum菌和Eggerthellalenta菌丰度呈正相关。van Tongeren等[11]采用衰弱指数(frailty index,FI)评价衰弱,发现非衰弱人群,菌群特点表现为普氏菌属(Prevotella)与拟杆菌属(Bacteroides)的比值较高;而在衰弱人群中,菌群特征表现为瘤胃球菌属(Ruminococcus)及奇异菌属(Atopobium)占比增加,乳酸杆菌属(Lactobacillus)占比减少。同时,Biagi团队[12]也发现慢性疾病越多的老年人,肠道内克里斯滕森菌科(Christensenellaceae)、嗜黏蛋白阿克曼菌(Akkermansia)及双歧杆菌属(Bifidobacterium)等数量减少,但这些细菌数量级与衰弱并非呈线性关系,个体差异性较大。
肠道菌群可从多个方面影响衰弱。首先,肠道微生物组成失调会影响蛋白质代谢和生物利用度[13-14],影响骨骼肌合成代谢所需的维生素B12、叶酸、核黄酸等物质的产生[14],加速肌细胞线粒体受损累积[15],增加胰岛素抵抗风险[16],这些因素会导致骨骼肌数量与质量的下降;其次,随着年龄增长,微生物组成向着促炎表型发展,总体表现为微生物多样性减少,个体变异增加[17],该变化与握力降低显著相关[18];另一方面,肠道菌群结构会随着地理环境、饮食结构及居住条件、年龄等因素发生改变[19],而目前尚无文献报道中国中老年糖尿病衰弱人群肠道菌群结构及功能特征。本研究将采用因美纳(Illumina)高通量测序等方式,对肠道菌群进行定性定量分析,与糖尿病非衰弱人群对比,探索糖尿病衰弱人群肠道菌群结构及功能特征,揭示菌群与糖尿病衰弱之间的相互关系。
1 材料与方法
1.1 对象
本研究于2017年5月至2017年7月在四川大学华西医院老年科门诊招募患者。纳入标准:① 已确诊为2型糖尿病的患者;② 愿意参加该研究并配合试验留取粪便。排除标准:① 终末期疾病,远期预后不超过半年;② 严重肝肾功能不全。衰弱评估采用FRAIL量表,该量表包括了五个条目:疲劳(Fatigue)、抗阻力(Resistance)、步行(Ambulation)、疾病(Illness)、体重下降(Loss weight),其中 ≥ 3项满足条件为衰弱,1~2项满足条件为衰弱前期,0项为没有衰弱。本研究试验组为糖尿病衰弱及衰弱前期患者,对照组按1∶1比例匹配相同年龄、性别纳入糖尿病非衰弱患者,收集临床信息及患者粪便。
本研究严格遵守赫尔辛基宣言,并通过四川大学华西医院伦理委员会审查,在参与任何实验活动之前,每位受试者均已签署知情同意书。
1.2 样品采集
所有入组患者均正常饮食,采用自然排便方式,于清晨使用粪便取样器截取样品中段里部,采取新鲜粪便标本,样品于2 h内收集完成,放置于 − 80 ℃冰箱保存,粪便样本统一由北京量化健康科技有限公司(干冰运输)进行粪便宏基因测序。
1.3 DNA提取
采用QIAamp PowerFecal Pro DNA Kit试剂对粪便样本DNA进行提取,并采用琼脂糖凝胶电泳(agarose gel electrophoresis,AGE)的方式分析DNA纯度和完整性;采用Qubit对DNA浓度进行精确定量。
1.4 文库构建及测序
采用KAPA HyperPlus PCR-free(96rxn)试剂盒进行文库构建,检测合格的DNA样品采用Covaris超声破碎仪随机打断成长度约350 bp的片段,经末端修复、加A尾、加测序接头、纯化、PCR扩增等步骤完成文库制备,文库构建完成后,使用Qubit2.0进行初步定量,稀释文库至2 ng/μL,随后使用普瑞凯Qsep 100对文库的嵌入部分(insert size)进行检测,符合预期后,再采用Q-PCR方法对文库有效浓度进行精确定量。把不同文库按照有效浓度及目标下机数据量需求,合并进行测序。
1.5 测序分析
所有原始宏基因组测序数据将通过MOCAT2软件进行质量控制[20],采用Cutadapt软件(v1.14,参数:-M30)对原始序列进行接头处理,使用SolexaQA package过滤掉长度小于30 bp的序列,质控后得到较为干净的序列。再使用SOAPaligner软件(v2.21,参数:-M4 -I30 -V10)比对宿主基因组去除污染的宿主基因序列,得到高质量的干净数据。采用SOAPdenovo软件对序列进行从头组装(v2.04,参数:all -D1 -M3 -L500),将干净序列打断成K-mer构建德布鲁因图(de Bruijin),寻找欧拉途径(Eulerian Path)组装成重叠群(contigs),根据双端测序(pair-end reads)位置关系连接成支架序列(scaffold),并从中挑选出连续的重叠群,得到不含氮端的序列片段(scaftigs)。
基因组装后获得的序列采用MetaGeneMark进行基因结构预测[21],预测后用CD-HIT的方式进行聚类去冗余[22],构建基因集,采用DIAMOND软件,将所获得的非冗余参考基因对比到京都基因组百科全书基因库(KEGG Orthology,KO),对所有预测出的基因进行功能注释。
在物种及其他高级分类层面上,采用MetaPhlan2.0获得菌群相对丰度,并用来自7 500个物种的100多万个标记基因做预测,获得不同分类水平下的相对丰度[23]。同时采用MOCAT2基于宏基因组操作分类单位(metagenomic operational taxonomic units,MetaOTUs)对基因组序列未知物种进行分析,获得MetaOTUs的相对丰度[24]。
对于KO的相对丰度,采用R/bioconductor的GAGE包进行KEGG Pathway富集分析[25]。
1.6 物种累积曲线及多样性分析
处理好后的菌群基因,采用R软件中vegan-specaccum函数计算物种累积曲线,对物种丰度进行预测;在属水平使用vegan-diversity函数计算香农指数(Shannon Index)以及辛普森指数(Simpson Index),用于评估物种阿尔法(alpha,α)多样性;采用基于布雷柯蒂斯距离(Bray-Curtis)的主坐标分析(principal co-ordinates analysis,PCoA)法,在二维坐标上显示不同样本之间的相似性及差异性,同时采用相似百分比分析法(similar percentage,SIMPER)确定主要差异菌群。
为比较两组间肠道菌群基因、功能和分类群差异,本研究通过威尔科克森符号秩检验(Wilcoxon)比较两组间相对丰度来得到显著差异单元,通过本贾米尼-霍克伯格法(Benjamini-Hochberg,HB)方法来矫正P值;并采用线性判别分析(linear discriminant analysis,LDA)对数据进行降维,评估差异显著物种的影响力(LDA Score),得到线性判别分析效果大小结果图(Linear discriminant analysis Effect Size,LEfSe)。
2 结果
2.1 纳入人群基本特征
本研究共纳入糖尿病衰弱或衰弱前期患者15名(A组),按照年龄、性别、衰弱程度等进行1∶1配对后纳入糖尿病非衰弱患者15名(B组)。具体情况见表1。
2.2 粪便样本测序基本统计结果
采用Illumina测序平台进行测序,共获得了193.9 GB的原始测序数据,样本的序列数量在34 956 712~59 919 586之间,平均为46 275 356,所有样本纯净序列占未处理序列的95%以上,单样本序列丰度 > 0.1%。
2.3 肠道微生物多样性
通过对两组共30例患者的粪便进行宏基因测序,共得到16门、127属、298种细菌。其中,糖尿病衰弱及衰弱前期组共提取出12门、116属、264种微生物;糖尿病非衰弱组共提取出16门、100属、243种微生物。通过物种累积曲线验证微生物测序数据量在多样性分析中是否充足,是反映样本物种多样性指标之一。横坐标代表样本量,纵坐标代表被检测的物种数,结果显示末端曲线趋于平缓,表明样本量足够(见图1)。
香农指数常用于反映α多样性指数,香农指数值越大,说明菌落多样性越高。对两组人群进行香农指数计算,如图2所示,A组较B组香农指数略小,但经秩和检验分析,两组α多样性差异无统计学意义。再采用辛普森指数进行计算,A、B两组间菌群多样性差异无统计学意义(见图3)。
本研究使用PCoA来分析两组人群粪便微生物β多样性(见图4),该方法通过计算Bray-Curtis距离,对菌群进行排序、降维,形成相似度及相异度可视化图表。结果显示,两组人群肠道菌群分布具有差异性,但各自聚集不明显。使用非参数多元方差分析检验分组因素对菌群差异的影响,结果显示衰弱对糖尿病患者菌群结构总体变化影响不显著,无统计学意义。
2.4 肠道微生物组成特征
如图5所示,在门水平上,糖尿病衰弱及衰弱前期组肠道菌群相对丰度较高的为厚壁菌门(Firmicutes)、拟杆菌门(Bacteroidetes)、变形菌门(Proteobacteria)、放线菌门(Actinobacteria);在科水平上,糖尿病衰弱及衰弱前期组肠道菌群相对丰度较高的为瘤胃菌科(Ruminococcaceae)、拟杆菌科(Bacteroidaceae)、韦荣氏球菌科(Veillonellaceae)、毛螺菌科(Lachnospir-aceae)、优杆菌科(Eubacteriaceae)、普雷沃氏菌科(Prevotellaceae)等;在属水平上,糖尿病衰弱及衰弱前期组患者肠道菌群相对丰度较高的为粪杆菌属(Faecalibacterium)、拟杆菌属(Bacteroides)、优杆菌属(Eubacterium);在种水平上,糖尿病衰弱组患者肠道菌群相对丰度较高的为普拉梭菌(Faecalibacterium prausnitzii)、直肠真杆菌(Eubacterium rectale)、Prevotellacopri菌、长双歧杆菌(Bifidobacterium longum)以及大肠杆菌(Escherichia coli)。
2.5 肠道菌群组间差异性分析
本研究采用秩和检验对两组患者肠道菌群各个水平相对丰度进行统计分析,如见表2所示。在门水平,糖尿病衰弱及衰弱前期组的拟杆菌门(Bacteroidetes,红色)相对丰度高于糖尿病非衰弱组(P = 0.016)。在纲水平,糖尿病衰弱及衰弱前期组的拟杆菌纲(Bacteroidia)相对丰度低于糖尿病非衰弱组(P = 0.02)。在目水平,糖尿病衰弱及衰弱前期Bacteroidales相对丰度低于糖尿病非衰弱组(P = 0.016)。在科水平,糖尿病衰弱及衰弱前期组Bacteroidales noname高于糖尿病非衰弱组(P < 0.01),而Bacteroidaceae在糖尿病非衰弱组中相对丰度较高(P = 0.03)。在属水平上,糖尿病衰弱及衰弱前期组Bacteroides_noname(P = 0.01)、Collinsella(P = 0.022)、Butyricimonas(P = 0.03)相对丰度高于糖尿病非衰弱组,而Parasutterella(P = 0.01)、Bacteroides(P = 0.03)在糖尿病非衰弱组中相对丰度较高。在种水平丰度上,衰弱及衰弱前期组人群Bacteroidales-bacterium-ph8(P < 0.01)、Butyricimonas_synergistica(P = 0.03)、Bacteroides_ovatus(P = 0.01)相对丰度显著高于糖尿病非衰弱组患者,而Parasutterella_excrementihominis(P = 0.02)在糖尿病非衰弱组中相对丰度较高。LEfSe分析结果显示糖尿病衰弱及衰弱前期患者肠道菌群贡献最大的菌为Butyricimonas菌属,而糖尿病非衰弱组人群关键菌包括Bacteroidetes门bacteroides属等(见图6、图7)。
2.6 差异功能分析
本研究从30例样本中共提取及构建了1 264 594非冗余基因集。其中,糖尿病衰弱及衰弱前期组人群共有1 147 676个基因集,糖尿病非衰弱组人群共有1 044 831个基因集。其中,共有基因集占有927 913个,占总基因集73.4%。两组人群微生物功能基因涉及多条功能通路,通过将糖尿病衰弱及衰弱前期和糖尿病非衰弱组人群差异性表达的KO注释到KEGG模块水平,发现糖尿病衰弱及衰弱前期的人群富集上调的功能基因分布在组氨酸代谢(histidine metabolism)、EB病毒感染(Epstein-Barr virus infection)、硫代谢(sulfur metabolism)以及Ⅱ型聚酮生物合成(biosynthesis of type Ⅱ polyketide products)等;而磷酸转移酶系统(phosphotransferase system,pts)、生物膜形成(biofilm formation)、ABC转运蛋白(ABC transporters)、鞭毛组装(flagellar assembly)、丁酸代谢(butanoate metabolism)以及苯丙氨酸代谢(phenylalanine metabolism)等在糖尿病衰弱及衰弱前期人群中出现下调(见表3)。
3 讨论
本研究调查了糖尿病衰弱人群的肠道微生物,发现糖尿病衰弱与肠道菌群α多样性并无相关性,两组人群β多样性呈现差异性分布,但各自微生物群聚集不明显。对差异菌群进行功能分析,发现衰弱人群肠道微生物在氨基酸代谢、EB病毒感染等方面上调,而在生物膜形成、丁酸代谢、苯丙氨酸代谢等方面下调。
3.1 菌群结构特征分析
肠道微生物会随着年龄增长而出现多样性下降或生态失衡,个体变异性增加,促使寄主生理紊乱显著[26],而这种变化与肠道通透性增加以及慢性低度炎症具有相关性[27]。已有证据表明衰弱状态伴随着系统性炎性[28]。既往研究中,衰弱人群的肠道微生物最一致的特征是多样性的丧失[27],本研究发现糖尿病衰弱及衰弱前期组的α多样性少于糖尿病非衰弱组人群,但两组并没有统计学意义,原因可能是糖尿病对肠道菌群的影响[29],导致在没有衰弱状态时就缺乏微生物多样性,β多样性分析也提示两组人群差异并没有统计学意义。
两组人群的菌群组成具有明显差异性。在糖尿病衰弱及衰弱前期组中,柯林斯氏菌属(Collinsella)和丁酸弧菌属(Butyricimonas)富集。既往研究提示Collinsella丰度与空腹甘油三酯、总胆固醇浓度呈正相关[30],促进体内脂肥胖形成,进一步促进体内炎症因子形成[31]。在给无菌小鼠肠道喂养Collinsella菌后发现肠细胞紧密连接蛋白表达降低,肠漏及菌血症的发生风险显著升高[32]。因此我们推测,Collinsella菌属可能通过改变体内慢性炎症环境,导致衰弱发病率上升,是潜在的糖尿病衰弱致病菌群。然而,目前尚不清楚该菌的富集是衰弱的原因还是结果,其潜在机制有待进一步研究。
在糖尿病非衰弱组中,副萨特氏菌属(Parasutterella)及卵形拟杆菌(Bacteroides ovatus)富集,目前的研究认为副萨特氏菌属是肠道的核心菌种,在胆汁酸维持和胆固醇代谢中具有重要作用[33]。在一项随机试验中,副萨特氏菌属与低密度脂蛋白(low density lipoprotein,LDL)降低具有相关性,而低水平LDL与较低的衰弱发病率具有相关性[34],因此我们推测副沙特氏菌是糖尿病衰弱的保护菌群。卵形拟杆菌在肠道中参与了多种碳水化合物的消化分解,提供能量[35],目前有研究已证明卵形拟杆菌能够诱导肠道产生免疫球蛋白A[36],限制细菌与细菌源性毒素进入肠上皮细胞,减少肠道感染风险[37],同时进食含有卵形拟杆菌的益生菌,能够改善社区人群腹泻等症状[38]。因此我们推测,卵形拟杆菌数量增加可能会降低糖尿病人群衰弱发生的风险,而相关数量级的确定,还需要更多的研究来支持。
3.2 代谢功能特征性分析
通过与KEGG代谢通路数据库进行比对,发现两组人群菌群代谢通路具有显著差异,在糖尿病衰弱及衰弱前期组人群中,上调的通路包括组氨酸代谢通路、EB病毒感染通路、硫代谢通路,以及Ⅱ型聚酮生物合成通路。组氨酸代谢通路中,咪唑丙酸是其重要产物,通过激活雷帕霉素靶蛋白敏感性复合体1通路(mammalian target of rapamycin complex 1,mTORC1),降解胰岛素受体底物1(Insulin Receptor Substrate 1,IRS1)和胰岛素受体底物2(Insulin Receptor Substrate 2,IRS2),抑制胰岛素信号通路,造成胰岛素抵抗,可能是2型糖尿病的病因之一[39];另一方面,胰岛素抵抗可以通过抑制骨骼肌细胞能量代谢造成肌肉收缩障碍[5]。这提示我们,组氨酸代谢的相关细菌所引起的胰岛素抵抗,可能在老年衰弱的发生发展中起到一定的作用。同时EB病毒感染通路能激活Toll样受体信号通路(Toll-like receptor signaling pathway,TLR pathway),促进体内白细胞介素6(interleukin-6,IL-6)及肿瘤细胞坏死因子α(tumor necrosis factor-α,TNF-α)的产生[40],而这些炎性标志物与衰弱具有相关性,因此推测衰弱老年人与EB病毒感染通路激活相关,而具体机制还需要更多研究支持。
在糖尿病衰弱及衰弱前期人群菌群通路中,下调的通路包括磷酸转移酶系统、生物膜形成、ABC转运蛋白、鞭毛组装及丁酸代谢。其中丁酸代谢在长寿人群当中具有重要意义[41]。同时丁酸通过结肠细胞上的单羧酸转运蛋白(monocarboxylate transporter,MCT-1)和溶质转运蛋白(solute carrier family 5 member 8,SLC5A8)转运,作为组蛋白脱乙酰酶抑制剂或G蛋白偶联受体分子,在抗炎、维持肠黏膜屏障以及肠黏膜免疫方面发挥重要作用[42]。糖尿病衰弱及衰弱前期人群肠道微生物对丁酸利用减少,上述丁酸作用相对减弱,前期的研究中发现,肠道免疫系统紊乱是衰弱发生发展的重要因素[43],因此我们推测丁酸可能在预防衰弱发生或改善衰弱程度等方面具有一定的作用,具体机制还需要更多研究证明。
3.3 创新性
既往的研究多集中在欧美等地国家,研究地点多选择在病房以及养老院等[10-11, 44],而关于社区衰弱人群的肠道菌群差异性在国内未见报道。由于肠道菌群受宿主种族、地域及饮食等的密切影响[45],不同地点人群肠道菌群特点值得研究。本研究采用人群年龄、性别完全配比的方式控制了混杂因素,发现柯林斯氏菌以及EB病毒在衰弱发生中起到重要作用,拓宽了以往的研究结果,为衰弱的干预提供了新的理论依据。
4 结论
本研究通过比较糖尿病衰弱及衰弱前期和非衰弱人群的肠道菌群差异,揭示了糖尿病衰弱及衰弱前期人群的肠道菌群结构特征。本研究采用宏基因测序对糖尿病衰弱及衰弱前期患者肠道菌群结构及功能特征性进行研究,结果表明衰弱及衰弱前期对肠道菌群整体结构影响不显著,但仍具有特异性菌群。通过对菌群基因功能分析,表明肠道菌群可通过胰岛素抵抗和肠道免疫等途径对糖尿病人群的衰弱产生一定影响。本研究的结果为衰弱人群的肠道菌群机制探索奠定了理论基础,同时也为糖尿病衰弱及衰弱前期人群的预防与干预提供了一定的科学依据。
利益冲突声明:本文全体作者均声明不存在利益冲突。
引言
全球老龄化问题正日益凸显,中国老龄化进程明显快于世界平均水平。2010年我国65岁及以上人群占比8.87%,预计到2050年,这一数字将提高至34%[1]。衰弱是老龄化突出问题之一,是指老年人生理储备下降导致机体易损性增加和抗应激能力减退的非特异性状态,其核心是老年人生理储备下降或异常,外界较小刺激即可引起临床事件的发生[2]。衰弱人群更容易出现各种不良预后,包括跌倒、失能、认识功能下降、再入院率增加甚至死亡风险等[3]。
糖尿病与衰弱两者可相互作用。一方面,高血糖状态会抑制骨骼肌生长,严重时会造成肌肉萎缩[4],而胰岛素抵抗会抑制骨骼肌的能量代谢造成肌肉再生功能障碍[5],机体衰弱发生风险增加[6]。另一方面,衰弱促使胰岛素抵抗增加,血糖呈高水平状态,增加糖尿病发生风险[7]。我国的流行病学调研显示,我国65岁以上糖尿病合并衰弱人群占糖尿病人群比例为5%~48%[8],该人群出现低血糖及死亡风险高[9],预后较单纯糖尿病老年人群差,极大地增加了家庭照护和社会经济负担。
微生物在人体内广泛分布,其中以胃肠道菌群数量最多。既往研究表明,衰弱与肠道菌群α多样性呈负相关[10],衰弱程度与Ruminococcacea菌、F. prausnitzii菌及Lachnospiraceae菌丰度呈负相关,与Eubacteriumdolichum菌和Eggerthellalenta菌丰度呈正相关。van Tongeren等[11]采用衰弱指数(frailty index,FI)评价衰弱,发现非衰弱人群,菌群特点表现为普氏菌属(Prevotella)与拟杆菌属(Bacteroides)的比值较高;而在衰弱人群中,菌群特征表现为瘤胃球菌属(Ruminococcus)及奇异菌属(Atopobium)占比增加,乳酸杆菌属(Lactobacillus)占比减少。同时,Biagi团队[12]也发现慢性疾病越多的老年人,肠道内克里斯滕森菌科(Christensenellaceae)、嗜黏蛋白阿克曼菌(Akkermansia)及双歧杆菌属(Bifidobacterium)等数量减少,但这些细菌数量级与衰弱并非呈线性关系,个体差异性较大。
肠道菌群可从多个方面影响衰弱。首先,肠道微生物组成失调会影响蛋白质代谢和生物利用度[13-14],影响骨骼肌合成代谢所需的维生素B12、叶酸、核黄酸等物质的产生[14],加速肌细胞线粒体受损累积[15],增加胰岛素抵抗风险[16],这些因素会导致骨骼肌数量与质量的下降;其次,随着年龄增长,微生物组成向着促炎表型发展,总体表现为微生物多样性减少,个体变异增加[17],该变化与握力降低显著相关[18];另一方面,肠道菌群结构会随着地理环境、饮食结构及居住条件、年龄等因素发生改变[19],而目前尚无文献报道中国中老年糖尿病衰弱人群肠道菌群结构及功能特征。本研究将采用因美纳(Illumina)高通量测序等方式,对肠道菌群进行定性定量分析,与糖尿病非衰弱人群对比,探索糖尿病衰弱人群肠道菌群结构及功能特征,揭示菌群与糖尿病衰弱之间的相互关系。
1 材料与方法
1.1 对象
本研究于2017年5月至2017年7月在四川大学华西医院老年科门诊招募患者。纳入标准:① 已确诊为2型糖尿病的患者;② 愿意参加该研究并配合试验留取粪便。排除标准:① 终末期疾病,远期预后不超过半年;② 严重肝肾功能不全。衰弱评估采用FRAIL量表,该量表包括了五个条目:疲劳(Fatigue)、抗阻力(Resistance)、步行(Ambulation)、疾病(Illness)、体重下降(Loss weight),其中 ≥ 3项满足条件为衰弱,1~2项满足条件为衰弱前期,0项为没有衰弱。本研究试验组为糖尿病衰弱及衰弱前期患者,对照组按1∶1比例匹配相同年龄、性别纳入糖尿病非衰弱患者,收集临床信息及患者粪便。
本研究严格遵守赫尔辛基宣言,并通过四川大学华西医院伦理委员会审查,在参与任何实验活动之前,每位受试者均已签署知情同意书。
1.2 样品采集
所有入组患者均正常饮食,采用自然排便方式,于清晨使用粪便取样器截取样品中段里部,采取新鲜粪便标本,样品于2 h内收集完成,放置于 − 80 ℃冰箱保存,粪便样本统一由北京量化健康科技有限公司(干冰运输)进行粪便宏基因测序。
1.3 DNA提取
采用QIAamp PowerFecal Pro DNA Kit试剂对粪便样本DNA进行提取,并采用琼脂糖凝胶电泳(agarose gel electrophoresis,AGE)的方式分析DNA纯度和完整性;采用Qubit对DNA浓度进行精确定量。
1.4 文库构建及测序
采用KAPA HyperPlus PCR-free(96rxn)试剂盒进行文库构建,检测合格的DNA样品采用Covaris超声破碎仪随机打断成长度约350 bp的片段,经末端修复、加A尾、加测序接头、纯化、PCR扩增等步骤完成文库制备,文库构建完成后,使用Qubit2.0进行初步定量,稀释文库至2 ng/μL,随后使用普瑞凯Qsep 100对文库的嵌入部分(insert size)进行检测,符合预期后,再采用Q-PCR方法对文库有效浓度进行精确定量。把不同文库按照有效浓度及目标下机数据量需求,合并进行测序。
1.5 测序分析
所有原始宏基因组测序数据将通过MOCAT2软件进行质量控制[20],采用Cutadapt软件(v1.14,参数:-M30)对原始序列进行接头处理,使用SolexaQA package过滤掉长度小于30 bp的序列,质控后得到较为干净的序列。再使用SOAPaligner软件(v2.21,参数:-M4 -I30 -V10)比对宿主基因组去除污染的宿主基因序列,得到高质量的干净数据。采用SOAPdenovo软件对序列进行从头组装(v2.04,参数:all -D1 -M3 -L500),将干净序列打断成K-mer构建德布鲁因图(de Bruijin),寻找欧拉途径(Eulerian Path)组装成重叠群(contigs),根据双端测序(pair-end reads)位置关系连接成支架序列(scaffold),并从中挑选出连续的重叠群,得到不含氮端的序列片段(scaftigs)。
基因组装后获得的序列采用MetaGeneMark进行基因结构预测[21],预测后用CD-HIT的方式进行聚类去冗余[22],构建基因集,采用DIAMOND软件,将所获得的非冗余参考基因对比到京都基因组百科全书基因库(KEGG Orthology,KO),对所有预测出的基因进行功能注释。
在物种及其他高级分类层面上,采用MetaPhlan2.0获得菌群相对丰度,并用来自7 500个物种的100多万个标记基因做预测,获得不同分类水平下的相对丰度[23]。同时采用MOCAT2基于宏基因组操作分类单位(metagenomic operational taxonomic units,MetaOTUs)对基因组序列未知物种进行分析,获得MetaOTUs的相对丰度[24]。
对于KO的相对丰度,采用R/bioconductor的GAGE包进行KEGG Pathway富集分析[25]。
1.6 物种累积曲线及多样性分析
处理好后的菌群基因,采用R软件中vegan-specaccum函数计算物种累积曲线,对物种丰度进行预测;在属水平使用vegan-diversity函数计算香农指数(Shannon Index)以及辛普森指数(Simpson Index),用于评估物种阿尔法(alpha,α)多样性;采用基于布雷柯蒂斯距离(Bray-Curtis)的主坐标分析(principal co-ordinates analysis,PCoA)法,在二维坐标上显示不同样本之间的相似性及差异性,同时采用相似百分比分析法(similar percentage,SIMPER)确定主要差异菌群。
为比较两组间肠道菌群基因、功能和分类群差异,本研究通过威尔科克森符号秩检验(Wilcoxon)比较两组间相对丰度来得到显著差异单元,通过本贾米尼-霍克伯格法(Benjamini-Hochberg,HB)方法来矫正P值;并采用线性判别分析(linear discriminant analysis,LDA)对数据进行降维,评估差异显著物种的影响力(LDA Score),得到线性判别分析效果大小结果图(Linear discriminant analysis Effect Size,LEfSe)。
2 结果
2.1 纳入人群基本特征
本研究共纳入糖尿病衰弱或衰弱前期患者15名(A组),按照年龄、性别、衰弱程度等进行1∶1配对后纳入糖尿病非衰弱患者15名(B组)。具体情况见表1。
2.2 粪便样本测序基本统计结果
采用Illumina测序平台进行测序,共获得了193.9 GB的原始测序数据,样本的序列数量在34 956 712~59 919 586之间,平均为46 275 356,所有样本纯净序列占未处理序列的95%以上,单样本序列丰度 > 0.1%。
2.3 肠道微生物多样性
通过对两组共30例患者的粪便进行宏基因测序,共得到16门、127属、298种细菌。其中,糖尿病衰弱及衰弱前期组共提取出12门、116属、264种微生物;糖尿病非衰弱组共提取出16门、100属、243种微生物。通过物种累积曲线验证微生物测序数据量在多样性分析中是否充足,是反映样本物种多样性指标之一。横坐标代表样本量,纵坐标代表被检测的物种数,结果显示末端曲线趋于平缓,表明样本量足够(见图1)。
香农指数常用于反映α多样性指数,香农指数值越大,说明菌落多样性越高。对两组人群进行香农指数计算,如图2所示,A组较B组香农指数略小,但经秩和检验分析,两组α多样性差异无统计学意义。再采用辛普森指数进行计算,A、B两组间菌群多样性差异无统计学意义(见图3)。
本研究使用PCoA来分析两组人群粪便微生物β多样性(见图4),该方法通过计算Bray-Curtis距离,对菌群进行排序、降维,形成相似度及相异度可视化图表。结果显示,两组人群肠道菌群分布具有差异性,但各自聚集不明显。使用非参数多元方差分析检验分组因素对菌群差异的影响,结果显示衰弱对糖尿病患者菌群结构总体变化影响不显著,无统计学意义。
2.4 肠道微生物组成特征
如图5所示,在门水平上,糖尿病衰弱及衰弱前期组肠道菌群相对丰度较高的为厚壁菌门(Firmicutes)、拟杆菌门(Bacteroidetes)、变形菌门(Proteobacteria)、放线菌门(Actinobacteria);在科水平上,糖尿病衰弱及衰弱前期组肠道菌群相对丰度较高的为瘤胃菌科(Ruminococcaceae)、拟杆菌科(Bacteroidaceae)、韦荣氏球菌科(Veillonellaceae)、毛螺菌科(Lachnospir-aceae)、优杆菌科(Eubacteriaceae)、普雷沃氏菌科(Prevotellaceae)等;在属水平上,糖尿病衰弱及衰弱前期组患者肠道菌群相对丰度较高的为粪杆菌属(Faecalibacterium)、拟杆菌属(Bacteroides)、优杆菌属(Eubacterium);在种水平上,糖尿病衰弱组患者肠道菌群相对丰度较高的为普拉梭菌(Faecalibacterium prausnitzii)、直肠真杆菌(Eubacterium rectale)、Prevotellacopri菌、长双歧杆菌(Bifidobacterium longum)以及大肠杆菌(Escherichia coli)。
2.5 肠道菌群组间差异性分析
本研究采用秩和检验对两组患者肠道菌群各个水平相对丰度进行统计分析,如见表2所示。在门水平,糖尿病衰弱及衰弱前期组的拟杆菌门(Bacteroidetes,红色)相对丰度高于糖尿病非衰弱组(P = 0.016)。在纲水平,糖尿病衰弱及衰弱前期组的拟杆菌纲(Bacteroidia)相对丰度低于糖尿病非衰弱组(P = 0.02)。在目水平,糖尿病衰弱及衰弱前期Bacteroidales相对丰度低于糖尿病非衰弱组(P = 0.016)。在科水平,糖尿病衰弱及衰弱前期组Bacteroidales noname高于糖尿病非衰弱组(P < 0.01),而Bacteroidaceae在糖尿病非衰弱组中相对丰度较高(P = 0.03)。在属水平上,糖尿病衰弱及衰弱前期组Bacteroides_noname(P = 0.01)、Collinsella(P = 0.022)、Butyricimonas(P = 0.03)相对丰度高于糖尿病非衰弱组,而Parasutterella(P = 0.01)、Bacteroides(P = 0.03)在糖尿病非衰弱组中相对丰度较高。在种水平丰度上,衰弱及衰弱前期组人群Bacteroidales-bacterium-ph8(P < 0.01)、Butyricimonas_synergistica(P = 0.03)、Bacteroides_ovatus(P = 0.01)相对丰度显著高于糖尿病非衰弱组患者,而Parasutterella_excrementihominis(P = 0.02)在糖尿病非衰弱组中相对丰度较高。LEfSe分析结果显示糖尿病衰弱及衰弱前期患者肠道菌群贡献最大的菌为Butyricimonas菌属,而糖尿病非衰弱组人群关键菌包括Bacteroidetes门bacteroides属等(见图6、图7)。
2.6 差异功能分析
本研究从30例样本中共提取及构建了1 264 594非冗余基因集。其中,糖尿病衰弱及衰弱前期组人群共有1 147 676个基因集,糖尿病非衰弱组人群共有1 044 831个基因集。其中,共有基因集占有927 913个,占总基因集73.4%。两组人群微生物功能基因涉及多条功能通路,通过将糖尿病衰弱及衰弱前期和糖尿病非衰弱组人群差异性表达的KO注释到KEGG模块水平,发现糖尿病衰弱及衰弱前期的人群富集上调的功能基因分布在组氨酸代谢(histidine metabolism)、EB病毒感染(Epstein-Barr virus infection)、硫代谢(sulfur metabolism)以及Ⅱ型聚酮生物合成(biosynthesis of type Ⅱ polyketide products)等;而磷酸转移酶系统(phosphotransferase system,pts)、生物膜形成(biofilm formation)、ABC转运蛋白(ABC transporters)、鞭毛组装(flagellar assembly)、丁酸代谢(butanoate metabolism)以及苯丙氨酸代谢(phenylalanine metabolism)等在糖尿病衰弱及衰弱前期人群中出现下调(见表3)。
3 讨论
本研究调查了糖尿病衰弱人群的肠道微生物,发现糖尿病衰弱与肠道菌群α多样性并无相关性,两组人群β多样性呈现差异性分布,但各自微生物群聚集不明显。对差异菌群进行功能分析,发现衰弱人群肠道微生物在氨基酸代谢、EB病毒感染等方面上调,而在生物膜形成、丁酸代谢、苯丙氨酸代谢等方面下调。
3.1 菌群结构特征分析
肠道微生物会随着年龄增长而出现多样性下降或生态失衡,个体变异性增加,促使寄主生理紊乱显著[26],而这种变化与肠道通透性增加以及慢性低度炎症具有相关性[27]。已有证据表明衰弱状态伴随着系统性炎性[28]。既往研究中,衰弱人群的肠道微生物最一致的特征是多样性的丧失[27],本研究发现糖尿病衰弱及衰弱前期组的α多样性少于糖尿病非衰弱组人群,但两组并没有统计学意义,原因可能是糖尿病对肠道菌群的影响[29],导致在没有衰弱状态时就缺乏微生物多样性,β多样性分析也提示两组人群差异并没有统计学意义。
两组人群的菌群组成具有明显差异性。在糖尿病衰弱及衰弱前期组中,柯林斯氏菌属(Collinsella)和丁酸弧菌属(Butyricimonas)富集。既往研究提示Collinsella丰度与空腹甘油三酯、总胆固醇浓度呈正相关[30],促进体内脂肥胖形成,进一步促进体内炎症因子形成[31]。在给无菌小鼠肠道喂养Collinsella菌后发现肠细胞紧密连接蛋白表达降低,肠漏及菌血症的发生风险显著升高[32]。因此我们推测,Collinsella菌属可能通过改变体内慢性炎症环境,导致衰弱发病率上升,是潜在的糖尿病衰弱致病菌群。然而,目前尚不清楚该菌的富集是衰弱的原因还是结果,其潜在机制有待进一步研究。
在糖尿病非衰弱组中,副萨特氏菌属(Parasutterella)及卵形拟杆菌(Bacteroides ovatus)富集,目前的研究认为副萨特氏菌属是肠道的核心菌种,在胆汁酸维持和胆固醇代谢中具有重要作用[33]。在一项随机试验中,副萨特氏菌属与低密度脂蛋白(low density lipoprotein,LDL)降低具有相关性,而低水平LDL与较低的衰弱发病率具有相关性[34],因此我们推测副沙特氏菌是糖尿病衰弱的保护菌群。卵形拟杆菌在肠道中参与了多种碳水化合物的消化分解,提供能量[35],目前有研究已证明卵形拟杆菌能够诱导肠道产生免疫球蛋白A[36],限制细菌与细菌源性毒素进入肠上皮细胞,减少肠道感染风险[37],同时进食含有卵形拟杆菌的益生菌,能够改善社区人群腹泻等症状[38]。因此我们推测,卵形拟杆菌数量增加可能会降低糖尿病人群衰弱发生的风险,而相关数量级的确定,还需要更多的研究来支持。
3.2 代谢功能特征性分析
通过与KEGG代谢通路数据库进行比对,发现两组人群菌群代谢通路具有显著差异,在糖尿病衰弱及衰弱前期组人群中,上调的通路包括组氨酸代谢通路、EB病毒感染通路、硫代谢通路,以及Ⅱ型聚酮生物合成通路。组氨酸代谢通路中,咪唑丙酸是其重要产物,通过激活雷帕霉素靶蛋白敏感性复合体1通路(mammalian target of rapamycin complex 1,mTORC1),降解胰岛素受体底物1(Insulin Receptor Substrate 1,IRS1)和胰岛素受体底物2(Insulin Receptor Substrate 2,IRS2),抑制胰岛素信号通路,造成胰岛素抵抗,可能是2型糖尿病的病因之一[39];另一方面,胰岛素抵抗可以通过抑制骨骼肌细胞能量代谢造成肌肉收缩障碍[5]。这提示我们,组氨酸代谢的相关细菌所引起的胰岛素抵抗,可能在老年衰弱的发生发展中起到一定的作用。同时EB病毒感染通路能激活Toll样受体信号通路(Toll-like receptor signaling pathway,TLR pathway),促进体内白细胞介素6(interleukin-6,IL-6)及肿瘤细胞坏死因子α(tumor necrosis factor-α,TNF-α)的产生[40],而这些炎性标志物与衰弱具有相关性,因此推测衰弱老年人与EB病毒感染通路激活相关,而具体机制还需要更多研究支持。
在糖尿病衰弱及衰弱前期人群菌群通路中,下调的通路包括磷酸转移酶系统、生物膜形成、ABC转运蛋白、鞭毛组装及丁酸代谢。其中丁酸代谢在长寿人群当中具有重要意义[41]。同时丁酸通过结肠细胞上的单羧酸转运蛋白(monocarboxylate transporter,MCT-1)和溶质转运蛋白(solute carrier family 5 member 8,SLC5A8)转运,作为组蛋白脱乙酰酶抑制剂或G蛋白偶联受体分子,在抗炎、维持肠黏膜屏障以及肠黏膜免疫方面发挥重要作用[42]。糖尿病衰弱及衰弱前期人群肠道微生物对丁酸利用减少,上述丁酸作用相对减弱,前期的研究中发现,肠道免疫系统紊乱是衰弱发生发展的重要因素[43],因此我们推测丁酸可能在预防衰弱发生或改善衰弱程度等方面具有一定的作用,具体机制还需要更多研究证明。
3.3 创新性
既往的研究多集中在欧美等地国家,研究地点多选择在病房以及养老院等[10-11, 44],而关于社区衰弱人群的肠道菌群差异性在国内未见报道。由于肠道菌群受宿主种族、地域及饮食等的密切影响[45],不同地点人群肠道菌群特点值得研究。本研究采用人群年龄、性别完全配比的方式控制了混杂因素,发现柯林斯氏菌以及EB病毒在衰弱发生中起到重要作用,拓宽了以往的研究结果,为衰弱的干预提供了新的理论依据。
4 结论
本研究通过比较糖尿病衰弱及衰弱前期和非衰弱人群的肠道菌群差异,揭示了糖尿病衰弱及衰弱前期人群的肠道菌群结构特征。本研究采用宏基因测序对糖尿病衰弱及衰弱前期患者肠道菌群结构及功能特征性进行研究,结果表明衰弱及衰弱前期对肠道菌群整体结构影响不显著,但仍具有特异性菌群。通过对菌群基因功能分析,表明肠道菌群可通过胰岛素抵抗和肠道免疫等途径对糖尿病人群的衰弱产生一定影响。本研究的结果为衰弱人群的肠道菌群机制探索奠定了理论基础,同时也为糖尿病衰弱及衰弱前期人群的预防与干预提供了一定的科学依据。
利益冲突声明:本文全体作者均声明不存在利益冲突。