论文部分内容阅读
目的:Graves病(GD)的发病机制目前尚未完全明确,但可以确定的是在遗传易感性基础上,环境因素导致的免疫功能紊乱已被证明是重要的组成部分,而环境因素与肠道菌群的组成和结构密切相关。本研究旨在分析甲巯咪唑治疗前后GD患者肠道菌群的多样性变化,以及探讨肠道菌群变化与甲状腺相关抗体的相关性。方法:选择遵义医科大学附属医院内分泌科初诊GD患者(UGD组,18例),经过甲巯咪唑(赛治)治疗后甲功接近或恢复正常(TGD组,10例),和同期遵义医科大学附属医院体检健康人群(NC组,11例)做为研究对象。采集以上研究对象病史资料,测量血压、身高、体重、甲状腺功能及其相关抗体、血常规、肝肾功能、糖脂代谢等常规及生化指标,计算体质指数(BMI),并收集新鲜粪便样本,提取肠道菌群DNA,16S rRNA基因扩增,Illumina平台测序,对测序结果进行物种注释、多样性及物种差异分析。采用Spearman秩检验对NC组及UGD组双歧杆菌属相对丰度与TRAb、TgAb、TPOAb进行相关性分析。结果:(1)Illumina测序共测得1562445条高质量序列,其中NC组中获得477940条序列;UGD组中获得了722725条序列;TGD组中获得了361780条序列;每个样本平均含有(40063±6355)条有效序列。根据获得的OTU丰度矩阵,使用R软件计算3组OTU的数量,并通过绘制OTUs的Venn图直观地呈现3组所共有和独有OTU所占数量。NC组特有409个OTU,UGD组特有525个OTU,TGD组特有191个OTU,其中3组共有2177个OTU相互重叠。3组间Chao1指数、ACE指数、Shannon指数及Simpson指数没有统计学差异(P>0.05),但GD组均低于NC组。基于非度量多维尺度(NMDS)的β多样性分析可以将NC组样本与UGD组和TGD组样本区分开来。(2)3组样本厚壁菌门(Firmicutes)、拟杆菌门(Bacteroidetes)、放线菌门(Actinobacteria)及蛋白菌门(Proteobacteria)为优势菌门,占细菌总数99%以上。UGD组Actinobacteria、TM7和Cyanobacteria丰度显著高于NC组,而Firmicutes丰度低于NC组(P<0.05);TGD组Proteobacteria和TM7丰度高于NC组,而[Thermi]丰度低于NC组(P<0.05);UGD组Actinobacteria丰度高于NC组(P<0.05)。(3)所有样本共检出19个丰度>1%属水平物种,其中拟杆菌属(Bacteroides)丰度最高。UGD组双歧杆菌属(Bifidobacterium)和柯林斯氏菌属(Collinsella)高于NC组,罗斯氏菌属(Roseburia)和小杆菌属(Dialister)低于NC组(P<0.05);TGD组Bacteroides高于NC组,普雷沃菌属(Prevotella)和Dialister低于NC组(P<0.05);UGD组Prevotella和Collinsella高于TGD组(P<0.05);(4)线性判别分析(LDA)结果表明UGD组中科水平上Bifidobacteriaceae、Clostridiaceae及Carnobacteriaceae相对于NC组优势丰度分布,而Porphyromonadaceae、Rikenellacaeae、_Paraprevotellaceae_、Veillonellaceae、Desulfovibrionaceae在NC组优势丰度分布。TGD组富集Xanthomonadates、Xanthomonadaceae、_Chromatiaceae_和Alteromonadales,而Leuconostocaceae主要富集于UGD组(LDA值>2,P<0.05)。(5)NC组及UGD组以Bifidobacterium为因变量的Spearman相关性分析显示:Bifidobacterium丰度与TRAb(r=0.588,P=0.002),TgAb(r=0.463,P=0.023),TPOAb(r=0.578,P=0.002)呈正相关。结论:GD患者肠道菌群生物多样性降低,正常人及甲巯咪唑治疗前后肠道菌群确有差异,主要表现为Collinsella等某些致病菌,Roseburia等产短链脂肪酸有益菌丰度的改变,而Bifidobacterium丰度的改变与TRAb、TgAb、TPOAb呈正相关,提示可能与GD的免疫机制有关。因此,我们推测肠道菌群紊乱可能参与了GD的发病。