大负荷运动后骨骼肌炎症反应的关键基因时序分析及信号通路鉴定.pdfVIP

中国运动医学2025年1月第44卷第1期ChinJSportsMed,Jan.2025Vol.44,No.129

大负荷运动后骨骼肌炎症反应的关键基因时序分析及信号通路鉴定

张燕！梁隆钰！夏雨²钱艳1丁海丽1.2

1成都体育学院运动医学与健康学院（四川成都610041）

2成都体育学院运动医学与健康研究所（四川成都610041）

摘要目的：基于生物信息学探讨大负荷运动后骨骼肌炎症反应基因的时间窗效应，鉴定运动性骨骼

肌损伤（exercise-inducedmuscledamage，EIMD）炎症反应的关键基因与信号通路。方法：48只SD大鼠

随机分为对照组（C组，n=8）和运动组（E组，n=40），分别于运动后即刻（E0）、12小时（E12）、24小时

（E24）、48小时（E48）、72小时（E72)取腓肠肌，通过转录组测序，筛选差异表达基因，对其基因本体论

(geneontology，GO）和京都基因与基因组百科全书代谢通路数据库（kyotoencyclopediaofgenesand

genomes,KECG）的注释进行富集分析。从数据库提取炎症相关基因，与差异表达基因（differentially

expressedgenes，DEGs）相交识别获取差异表达炎症相关基因（differentiallyexpressedinflammation-

relatedgenes,DEIRGs），基于Mfuzz算法对DEIRGs进行时间序列聚类获得具有相似表达特征的子集，并

进行CO和KEGG分析以及蛋白互作网络分析，筛选获得关键DEIRGs，二次功能富集，并进行关键基因表

达时间变化分析，确定关键信号通路。结果：大负荷运动后的骨骼肌炎症反应基因经Mfuzz时序聚类获得

7个DEIRGs子集（cluster），cluster5的差异基因表达总体下调，cluster7的差异基因表达总体上调，

cluster3和4的差异基因表达在E0表达上调后在E12快速下调，cluster2和6的差异基因表达在E0表达

下调后在E12快速上调，cluster1的差异基因表达在E0上调后在E12快速下调,E24再次上调，筛选鉴定

出关键差异表达炎症基因TP53、STAT3、CD44、AKT1、KDR、GJA1、CYCS、HIF1A、IQGAP3、FASN、TFRC，富

集在低氧诱导因子1（hypoxia-induciblefactor-1,HIF-1）、细胞凋亡、铁死亡、丝裂原活化蛋白激酶

(mitogen-activatedproteinkinase,MAPK）、血管内皮生长因子（vascularendothelial-derivedgrowth

factorB

，VEGF）、磷脂酰肌醇3-激酶-蛋白激酶

AktO(forkheadboxO,FoxOAMP-activatedprotein

）、胰岛素抵抗、叉头盒转录因子）、腺苷酸激活酶（

kinase，AMPK）和Janus激酶-信号转导与转录激活子(Januskinasesignaltransductionandactivatorof

transcription,JAK-STAT)信号通路。结论：炎症相关基因呈时间动态变化参与EIMD，运动后12小时呈现

明显的时间窗效应，关键靶点STAT3、AKT1、HIF-1A经JAK-STAT、PI3K-Akt、HIF-1、VECF信号通路响应

运动诱发的骨骼肌损伤，促进肌组织修复。

关键词运动；骨骼肌损伤；差异表达基因；炎症反应；时序分析

Timeseriesanalysisofkeygenesandidentificationofsignalingpathwaysinskeletalmuscle

infla