0 引言
疲劳是一种常见的生理与病理交叉性症状。世界卫生组织将其描述为个体在体力活动、脑力运转或情绪调节过程中,出现的持续性功能衰退状态,表现为躯体乏力、精神萎靡、注意力涣散及工作效率显著下降,且无法通过短期休息完全缓解
[1]。随着社会节奏加快、工作压力增大及生活方式改变,我国疲劳人群规模持续扩大,并且疲劳对个体与社会的多重影响已成为不容忽视的公共卫生问题
[2-3]。疲劳病因复杂,缺乏针对性治疗,尽管现代医学在抗疲劳方面取得了一定进展,但现有药物普遍存在副作用大、疗效不稳定等问题
[4-5]。因此,亟需探索安全、有效且针对性强的治疗手段。
肠道菌群是环境因素影响人类健康与疾病的核心枢纽,其组成或功能的失调是引发多种病理状态的直接推动因素
[6-7]。近年来,研究者发现肠道菌群与疲劳之间存在密切关系。研究结果显示,乳酸菌与双歧杆菌通过良性调节肠道微生态、维护屏障完整性,协同改善肠道症状和机体疲劳
[8-11]。此外,肠道中厚壁菌门数目减少,变形杆菌门增加,肠道屏障被破坏引起患者产生疲劳感
[12-13]。由此可知,肠道菌群失调可引发肠道局部代谢异常与物理屏障破坏,促使菌群及其代谢产物易位入血;进而触发系统性炎症与免疫失衡,引起疲劳
[14]。因此,通过靶向调节肠道菌群可能是缓解疲劳的有效治疗手段。
随着孟德尔随机化、网络药理学等生物信息学方法的快速发展,为疾病诊断、药物研发、靶点筛选等领域提供了高效的技术支撑,尤其在复杂疾病机制解析与创新治疗策略探索中发挥着不可替代的作用。孟德尔随机化(Mendelian randomization, MR)是以遗传学原理为基础,利用单核苷酸多态性(single nucleotide polymorphisms, SNPs)作为工具变量(instrumentalvariables, IVs)进行因果推断的流行病学方法
[15-16]。其核心优势在于能通过遗传工具变量模拟随机对照试验,从而在关联估计中有效避免混杂偏倚和反向因果。本研究以肠道菌群为暴露因素,疲劳作为结局变量,分析两者的因果关联。网络药理学可构建化合物-靶点-疾病的相互作用网络,整合多数据库信息筛选肠道菌群抗疲劳的关键靶点。
因此,本研究旨在利用孟德尔随机化方法确定具有抗疲劳作用的肠道菌群,其次利用网络药理学深入探究肠道菌群在抗疲劳中的作用机制,为靶向肠道菌群以缓解疲劳的潜在干预策略提供了初步的理论依据。
1 资料与方法
1.1 孟德尔随机化研究设计
两样本 MR分析的3个关键假设,①相关性假设: IVs与肠道菌群密切相关;②独立性假设:IVs独立于已知的混杂因素;③排他性假设:IVs不能直接影响结局,只能通过与肠道菌群的关联间接影响疲劳。
1.2 孟德尔随机化数据来源
本研究中以肠道菌群作为暴露因素,疲劳作为结局变量进行MR分析。肠道菌群遗传基因数据来自荷兰微生物组项目(Dutch Microbiome Project),队列包含7 738名参与者,测定了 207个微生物分类群。所有数据在NHGRI-EBI GWAS Catalog上的编号为GCST90027446-GCST90027857。疲劳GWAS数据来源于IEU OpenGWAS数据库,数据集为finn-b-R18_MALAI_FATIG,该研究共有180 362个样本,包括7 363例疲劳人群和172 999例对照,所有样本为欧洲血统。
1.3 工具变量选择
为了满足MR分析的三个假设,通过质控筛选与暴露因素密切相关并符合标准的IVs。筛选标准:① 使用p <1×进行IVs筛选。② 采用R包ieugwasr的ld clump函数去除连锁不平衡的工具变量,疲劳参数设置:r2=0.001;kb=10 000,肠道菌群参数设置:r2=0.001;kb=10 000。③ 采用F统计量测量IVs强度,去除F统计值小于10的弱工具变量。
1.4 分析方法
为了保证结果的可靠性,本研究使用Inverse variance weighted (IVW)、MR Egger、Weighted median、Simple mode、Weighted mode 这5种方法评估潜在的因果效应,以 IVW方法为主。通过Egger回归等方法进行敏感性分析,以确认结果的可靠性。
1.5 菌群特征代谢物及靶点筛选
利用gutMGENE数据库获取与疲劳呈负相关肠道菌群的代谢产物及作用靶点。
1.6 疲劳相关靶点筛选
以“fatigue”为关键词,在drugbank、Genecard、OMIM数据库中筛选相关疾病靶点,获得疲劳相关靶点。之后,利用R语言绘制疾病靶点与代谢产物靶点Venn图,获得交集靶点。
1.7 特征代谢产物抗疲劳核心靶点的筛选
将特征代谢产物作用靶点与疲劳共有靶点输入STRING数据库,获得PPI网络。然后将数据导入Cytoscape.3.10.3,运用centiscape插件进行拓扑分析。依据度中心值(Degree>48.94)、紧密中心值(Closeness>0.0031)、中介中心值(Betweenness>297.16)等3个指标筛选核心靶点。
1.8 富集分析
基于DAVID数据库对交集基因进行基因本体(Gene Ontology, GO)功能分析和京都基因和基因组百科全书(Kyoto Encyclopedia of Genes and Genomes, KEGG)通路富集分析,所得结果利用微生信平台(
https://www.bioinformatics.com.cn/)进行可视化。
2 结果
2.1 MR分析
为判断肠道菌群与疲劳之间的因果关系,进行两样本MR分析。IVW结果显示,12种肠道菌群与疲劳之间存在因果关系(
表1)。其中,其中包括8种保护因素,即OR<1,包括
g_Akkermansia(OR:0.877,95%CI:0.795
0.967)、
s_Bifidobacterium_bifidum (OR:0.940,95%CI:0.889
0.995)、
s_Akkermansia_muciniphila (OR:0.863,95%CI:0.780
0.955)、
s_Roseburia_hominis(OR:0.821,95%CI:0.700
0.964)等,与疲劳呈负相关;以及4种危险因素,即OR>1,包括
g_Flavonifractor(OR:1.149,95%CI:1.027
1.285)、
s_Bilophila_wadsworthia(OR:1.178,95%CI:1.017
1.363)、
s_Bacteroides_plebeius(OR:1.122,95%CI:1.038
1.213)、
s_Lachnospiraceae_bacterium_8_1_57FAA(OR:1.064,95%CI:1.009
1.122)等,与疲劳呈正相关。此外,5种不同分析方法中,上述12种肠道菌群的所有β值方向基本一致(
图1)。
2.2 敏感性分析
为确认上述结果的可靠性,进行敏感性分析,包括三个检验:异质性检验、水平多效性检验、逐一剔除检验
[17-19]。异质性检验旨在识别群体间因果效应的不一致性,结果如
表2所示,各菌群与疲劳之间的异质性检验pval>0.05,表明不存在异质性。水平多效性检验用于评价因果效应估计偏差,结果如
表3所示,肠道菌群的MR分析水平多效性pval>0.05,说明没有水平多效性效应。逐一剔除检验用于评估单个SNP是否对整体MR产生了过度影响,结果如
图2所示,通过逐步剔除每个SNP,其剩余SNP对结局变量的效应变化不大。由此可知,本研究中不存在混杂因素,且结果可靠、稳定性良好。
2.3 菌株特征代谢物及作用靶点确定
基于孟德尔随机化结果,以降低疲劳风险的肠道菌群进行深入探究,在gutMGENE数据库筛选特征代谢产物及作用靶点。其中,
c_Verrucomicrobiae,f_Verrucomicrobiaceae,o_Verrucomicrobiales,p_Verrucomicrobia四个水平微生物未查到对应特征代谢物,
g_Akkermansia,
s_Akkermansia_muciniphila,
s_Bifidobacterium_bifidum,
s_Roseburia_hominis分别查找到4,2,2,1种代谢产物,其次,分别查找特征代谢产物对应基因靶点,共40个(
图3)。由图可知,这40个靶点主要是Acetate、Proplonate、Trimethylamine oxide三种代谢产物的作用靶点。
2.4 疲劳靶点筛选
以fatigue为关键词筛选疲劳靶点。在drugbank、Genecard、OMIM数据库中筛选疾病靶点,共得到3 571个,去除重复基因后,得到疲劳疾病靶点3 553个。
2.5 核心靶点筛选
利用R语言绘制代谢产物靶点和疲劳靶点Venn图。由
图4可知,代谢产物和疲劳共有靶点30个。取交集靶点上传至STRING构建PPI网络,结果显示共有30个节点,186条边,平均度值12.4。经拓扑分析,得到9个核心靶点(
图5)。基于
图2,推测这9个靶点可能是
s_Akkermansia_muciniphila和
s_Roseburia_hominis抗疲劳的关键靶点。
2.6 富集分析
利用DAVID数据库对交集基因进行GO功能富集分析和KEGG通路富集分析。GO分析包括生物过程(BP)、细胞组分(CC)和分子功能(MF)三个方面。结果如
图6所示,在BP中,交集基因主要富集于炎症反应(inflammatory response)、基因表达的正向调控(positive regulation of gene expression)等过程;在CC 中,交集基因主要富集于胞外区(extracellular region)、细胞外间隙(extracellular space)等组分;在MF中,交集基因主要富集于细胞因子活性(cytokine activity)、蛋白连接(protein binding)等功能。以Fold Enrichment为依据对交集基因富集的KEGG通路进行排名,前20的通路结果如
图7所示,关键通路包括疟疾(Malria)、炎症性肠病(Inflammatory bowel disease)、类风湿性关节炎(Rheumatoid arthritis)、Toll样受体信号通路(Toll-like receptor signaling pathway)等。
3 讨论
随着生活节奏加快、竞争压力增大,疲劳成为常见的亚健康状态,其特征是疲惫乏力、能量不足等,严重影响人们的日常生活
[20-21]。近年来,以肠道菌群为靶点探究抗疲劳方式的研究日益增多。早期研究显示疲劳患者肠道菌群多样性降低,组成及丰富度改变,随着研究的深入发现,疲劳还与肠道菌群代谢功能紊乱有关
[22-23]。为进一步研究肠道菌群抗疲劳作用机制,本文通过MR分析和网络药理学,系统研究了肠道菌群抗疲劳的作用机制,为靶向肠道菌群抗疲劳提供了新的研究思路。结果表明,肠道中
s_Akkermansia_muciniphila等微生物的特征代谢产物可作用于核心靶点,进一步调节相关信号通路发挥抗疲劳作用。
首先,利用MR分析,阐释肠道菌群与疲劳之间的因果关系。通过GWAS数据,发现
g_Akkermansia,
s_Bifidobacterium_bifidum,
s_Akkermansia_muciniphila,
s_Roseburia_hominis等12种与疲劳呈负相关的肠道菌群,以及
g_Flavonifractor,
s_Bilophila_wadsworthia,
s_Bacteroides_plebeius,
s_Lachnospiraceae_bacterium_8_1_57FAA等4种与疲劳呈正相关的肠道菌群,这一结果从遗传学角度识别了与疲劳风险相关的肠道微生物。
s_Bifidobacterium_bifidum属于双歧杆菌属,具有调节肠道免疫、改善肠道屏障损伤等功能
[24-25],虽然尚未见直接报道与抗疲劳相关的内容,但其作为传统益生菌具有较大的潜力。
s_Roseburia_hominis是一种新型益生菌,当前研究主要集中于改善肠炎,可作为改善肠道疾病的潜在靶点
[26-27]。除了这三种微生物之外,其余微生物是疣微菌门(
p_Verrucomicrobia)到嗜黏蛋白阿克曼菌(
s_Akkermansia_muciniphila)不同分类水平的物种,表明嗜黏蛋白阿克曼菌可能是抗疲劳的关键菌种。
s_Akkermansia_muciniphila是肠道中调节宿主代谢功能的潜在益生菌,在缓解多种疾病中具有重要作用
[28]。MR分析确立了肠道菌群与疲劳间的因果关联。该方法有效控制了因素混杂偏倚,弥补传统观察性研究的固有局限。
其次,通过gutMGENE数据库获得降低疲劳风险菌群的特征代谢产物9种及40个靶点,利用多种数据库挖掘得到3 553个疲劳疾病靶点,并确定了与特征代谢产物的交集靶点30个。进一步利用PPI网络分析得到特征代谢产物发挥作用的关键靶点,包括:IL6、TNF、CXCL8、FOXP3、IL1β、TLR2、CCL2、JUN、CD163。其中,IL6、TNF、IL1β等炎症因子在疲劳患者体内显著升高,其原因可能是肠道中致病菌增加,损害肠屏障,致使全身炎症激活,产生疲劳
[29-30]。经GO功能富集分析以及KEGG通路富集分析,发现这些靶点在Toll样受体信号通路、IL 17通路等经典信号通路发挥着重要作用,阐明了肠道菌群特征代谢产物抗疲劳的多靶点及多通路的协同作用机制。
综上所述,本研究运用孟德尔随机化及网络药理学方法,深入探讨了肠道菌群抗疲劳的作用机制,为治疗疲劳提供了新的科学证据。未来,我们将继续深入探索这些靶点的具体作用机制,以期开发出更加安全、有效的抗疲劳药物。
本研究的发现虽初步阐释了肠道菌群抗疲劳的分子机制,但存在一定局限性。首先,MR分析所采用的肠道菌群与疲劳GWAS数据均来自欧洲人群(荷兰微生物组项目与FinnGen 研究),这可能限制结论向其他人群的外推——遗传背景、饮食及环境等因素均可影响肠道菌群构成及其与宿主表型的关联,因此本研究结果在不同人群中的普适性尚需进一步验证。需指出的是,目前大规模、高质量的非欧洲人群肠道菌群GWAS数据仍十分有限,本研究基于欧洲人群数据进行因果推断,是在现有最佳证据基础上的有益探索。此外,受限于菌群及其特征代谢物筛选的广度,本研究在机制阐释的全面性上仍有不足。未来应整合多组学技术,基于更完整的微生物及代谢物谱,系统构建肠道菌群调控疲劳的作用网络。
山西省基础研究计划青年项目(202403021222289)
山西中医药大学博士启动基金(2023BK29)
山西中医药大学优秀博士毕业生来晋工作奖励经费科研启动基金(2023BKS18)
山西省卫生健康委科研课题(2024164)
中药生命组学与创新药物研发研究室(zyyyjs2024019)