UQCR10和TNNT1:通过整合转录组分析和机器学习识别的肌少症新型生物标志物

UQCR10 and TNNT1: novel biomarkers for sarcopenia identified through integrated transcriptomic analysis and machine learning.

作者信息Wenhan Sun, Zhitao Shangguan, Jiandong Li, Xiaoqing Ye, Qiong Lin, Gang Chen
PMID41320776
发布时间2025-12-01
DOI10.1186/s40001-025-03471-w

实验完整度

中

包含转录组数据分析、机器学习和独立队列验证,并通过qPCR和Western blot验证蛋白和mRNA表达水平,但缺乏功能验证和动物实验。

主要模型

骨骼肌组织样本(GSE226151和GSE111016) 20例肌少症患者与20例健康对照的骨骼肌样本(qPCR验证) 4例肌少症患者与4例健康对照的骨骼肌样本(Western blot验证)

重点核对

GSE226151数据集(20肌少症 vs 20健康)用于差异表达分析和WGCNA GSE111016数据集(20肌少症 vs 20健康)用于独立验证 使用LASSO、XGBoost和随机森林进行特征选择 qPCR验证:20肌少症 vs 20健康对照 Western blot验证:4肌少症 vs 4健康对照

摘要

背景 肌少症是随年龄增长而发生的肌肉质量和功能丧失,显著影响老年人的健康结局。尽管其普遍存在,但由于对其病理生理学了解不完全,分子诊断和靶向治疗仍然有限。 方法 我们分析了GSE226151数据集,比较肌少症(n=20)和健康成人(n=20)的骨骼肌转录组,使用差异表达分析和加权基因共表达网络分析(WGCNA)。通过差异表达基因与枢纽基因的交集鉴定关键基因。采用机器学习特征选择(LASSO、XGBoost和随机森林)识别最佳生物标志物,并在独立队列(GSE111016,n=40,20例肌少症对比20例健康对照)中进行验证。使用基因本体论、KEGG和Reactome数据库进行全面的功能富集分析。我们采用五折交叉验证以应对适中的样本量,性能指标通过自举法(1000次迭代)获得的95%置信区间报告。 结果 差异表达分析(97个DEG)和WGCNA(来自darkolivegreen模块的269个枢纽基因)的整合确定了肌少症中的五个关键基因:COX6C、ETFB、TNFAIP3、TNNT1和UQCR10。多方法机器学习特征选择一致将UQCR10和TNNT1列为最重要的生物标志物。在验证队列中,双生物标志物组合(UQCR10+TNNT1)在随机森林模型中取得了优越的性能(AUC=0.738 [95% CI 0.559–0.875],PRAUC=0.643 [95% CI 0.405–0.832],敏感性=0.850,特异性=0.600),显著优于单一生物标志物方法。三组分析(健康/前肌少症/肌少症)揭示了连续的剂量反应模式:所有五个基因在疾病进展中均显示出显著的线性趋势(p<0.005),其中UQCR10(线粒体)和TNNT1(收缩)表现出最强的相关性(R²=0.281和0.192),提示连续的病理生理进展而非离散的疾病阶段。功能富集分析揭示了三个主要失调的生物学过程:能量代谢和线粒体功能(UQCR10、COX6C、ETFB)、肌肉结构和收缩功能(TNNT1)以及炎症反应调节(TNFAIP3)。 结论 我们的整合转录组方法确定了用于肌少症检测的有前景的双生物标志物组合(UQCR10+TNNT1),并阐明了该病症的潜在关键分子机制。虽然这些发现是假设生成的,需要在更大、更多样化的队列中验证,但它们提供了机制见解和值得进一步开发的候选生物标志物。所识别的基因和通路为旨在保持衰老期间肌肉健康的新型治疗干预提供了潜在靶点。

实验结论

提炼研究问题、关键发现与证据,快速把握文章的核心贡献。

研究问题
如何通过整合转录组分析和机器学习识别肌少症的潜在生物标志物并阐明其分子机制?
核心机制
肌少症涉及线粒体功能障碍(UQCR10、COX6C、ETFB下调)、肌肉结构改变(TNNT1下调)和炎症反应(TNFAIP3上调)的相互作用。
主要证据
差异表达分析和WGCNA整合出五个关键基因,机器学习筛选出UQCR10和TNNT1,在独立验证队列中随机森林模型AUC=0.738,qPCR和Western blot验证两者在肌少症中表达降低。
研究意义
该双生物标志物组合可能用于肌少症的早期检测、风险分层和治疗监测,并为靶向线粒体生物能量、收缩重塑和炎症调节的干预提供潜在靶点。

研究路径

按研究推进顺序梳理实验设计、验证步骤与关键观察。

1

差异表达分析

识别肌少症患者与健康对照之间差异表达的基因。

使用DESeq2对GSE226151数据集进行差异表达分析,筛选标准为调整p值≤0.05且|log2FC|≥log2(1.6),共获得97个DEGs。

2

加权基因共表达网络分析

识别与肌少症相关的共表达模块及枢纽基因。

使用WGCNA构建无符号网络,软阈值选择基于尺度自由拓扑拟合指数R²>0.85,最小模块大小30,合并阈值0.25,得到与肌少症最相关的darkolivegreen模块(含1066个基因),基于MM>0.8和|GS|>0.2筛选出269个枢纽基因。

3

整合分析确定关键基因

整合差异表达基因和枢纽基因,筛选与肌少症相关的关键基因。

通过Venn图分析,取DEGs和darkolivegreen模块枢纽基因的交集,得到5个关键基因:COX6C、ETFB、TNFAIP3、TNNT1和UQCR10。

4

疾病进展分析

评估关键基因在健康、前肌少症和肌少症三组中的表达变化趋势。

使用全部60个样本进行三组分析,采用DESeq2设计公式~group,进行线性趋势分析(多项式对比),发现五个基因均呈现显著线性趋势(p<0.005)。

5

功能富集分析

揭示关键基因参与的生物学过程和通路。

使用clusterProfiler进行GO、KEGG和Reactome富集分析,采用宽松阈值(p≤0.1,q≤0.2),鉴定出能量代谢/线粒体功能、肌肉结构/功能、炎症/免疫反应等关键通路。

6

机器学习特征选择

从候选基因中筛选最佳生物标志物组合。

对五个关键基因在发现队列(GSE226151)上使用LASSO、XGBoost和随机森林进行特征选择,通过共识排序确定UQCR10和TNNT1为前两名。

7

模型构建与独立验证

评估单生物标志物和双生物标志物模型的诊断性能。

基于UQCR10和UQCR10+TNNT1配置,使用逻辑回归、XGBoost和随机森林构建六个模型,在验证队列(GSE111016)中评估性能,双生物标志物随机森林模型性能最佳(AUC=0.738)。

8

实验验证

验证UQCR10和TNNT1在肌少症中的表达水平。

采用Western blot检测蛋白表达(4肌少症 vs 4健康),qPCR检测mRNA表达(20肌少症 vs 20健康),结果显示肌少症中UQCR10和TNNT1的表达均显著降低。

研究方法

按研究目的归类文中使用的方法,便于定位所需技术。

产品清单

实验环节名称品牌货号
RIPA裂解液Elabscience Biotechnology Inc.--
BCA蛋白定量试剂盒Beijing Solabio Technology Co., Ltd.--
PVDF膜----
辣根过氧化物酶标记的山羊抗兔IgG二抗Abcamab6721
ECL化学发光底物试剂盒Shanghai Yesen Biotechnology Co., Ltd.--
ChemiDoc成像系统Bio-Rad Laboratories Inc.--
TRIzol试剂Invitrogen--
NanoDrop分光光度计Thermo Fisher Scientific--
cDNA合成试剂盒Takara Bio--
QuantStudio 5实时荧光定量PCR系统Thermo Fisher Scientific--
SYBR Green预混液Applied Biosystems--

关键环节

汇总复现实验时建议重点确认的条件及原文阅读提示。

环节核对要点
差异表达分析
数据集GSE226151(20肌少症 vs 20健康);筛选阈值:adjusted p-value≤0.05,|log2FC|≥log2(1.6)
阅读提示:Methods: Differential expression analysis
加权基因共表达网络分析
软阈值选择(R²>0.85);最小模块大小30;模块合并阈值0.25;枢纽基因阈值MM>0.8,|GS|>0.2
阅读提示:Methods: Weighted gene co-expression network analysis
机器学习特征选择
发现队列GSE226151;五折交叉验证;LASSO lambda=0.204;XGBoost和随机森林参数优化
阅读提示:Methods: Machine learning analysis
模型验证
验证队列GSE111016(20肌少症 vs 20健康);性能指标AUC、PRAUC、敏感性、特异性;置信区间通过1000次自举法计算
阅读提示:Methods: Model development and evaluation
qPCR验证
样本量:20肌少症 vs 20健康;引物序列未在文中明确说明;内参GAPDH;使用2−ΔΔCt方法
阅读提示:Methods: qPCR
Western blot验证
样本量:4肌少症 vs 4健康;上样量30μg;一抗信息未明确提供;二抗稀释比例1:2000;内参GAPDH
阅读提示:Methods: Western blotting