差异表达分析
识别肌少症患者与健康对照之间差异表达的基因。
使用DESeq2对GSE226151数据集进行差异表达分析,筛选标准为调整p值≤0.05且|log2FC|≥log2(1.6),共获得97个DEGs。
UQCR10 and TNNT1: novel biomarkers for sarcopenia identified through integrated transcriptomic analysis and machine learning.
包含转录组数据分析、机器学习和独立队列验证,并通过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健康对照
提炼研究问题、关键发现与证据,快速把握文章的核心贡献。
按研究推进顺序梳理实验设计、验证步骤与关键观察。
识别肌少症患者与健康对照之间差异表达的基因。
使用DESeq2对GSE226151数据集进行差异表达分析,筛选标准为调整p值≤0.05且|log2FC|≥log2(1.6),共获得97个DEGs。
识别与肌少症相关的共表达模块及枢纽基因。
使用WGCNA构建无符号网络,软阈值选择基于尺度自由拓扑拟合指数R²>0.85,最小模块大小30,合并阈值0.25,得到与肌少症最相关的darkolivegreen模块(含1066个基因),基于MM>0.8和|GS|>0.2筛选出269个枢纽基因。
整合差异表达基因和枢纽基因,筛选与肌少症相关的关键基因。
通过Venn图分析,取DEGs和darkolivegreen模块枢纽基因的交集,得到5个关键基因:COX6C、ETFB、TNFAIP3、TNNT1和UQCR10。
评估关键基因在健康、前肌少症和肌少症三组中的表达变化趋势。
使用全部60个样本进行三组分析,采用DESeq2设计公式~group,进行线性趋势分析(多项式对比),发现五个基因均呈现显著线性趋势(p<0.005)。
揭示关键基因参与的生物学过程和通路。
使用clusterProfiler进行GO、KEGG和Reactome富集分析,采用宽松阈值(p≤0.1,q≤0.2),鉴定出能量代谢/线粒体功能、肌肉结构/功能、炎症/免疫反应等关键通路。
从候选基因中筛选最佳生物标志物组合。
对五个关键基因在发现队列(GSE226151)上使用LASSO、XGBoost和随机森林进行特征选择,通过共识排序确定UQCR10和TNNT1为前两名。
评估单生物标志物和双生物标志物模型的诊断性能。
基于UQCR10和UQCR10+TNNT1配置,使用逻辑回归、XGBoost和随机森林构建六个模型,在验证队列(GSE111016)中评估性能,双生物标志物随机森林模型性能最佳(AUC=0.738)。
验证UQCR10和TNNT1在肌少症中的表达水平。
采用Western blot检测蛋白表达(4肌少症 vs 4健康),qPCR检测mRNA表达(20肌少症 vs 20健康),结果显示肌少症中UQCR10和TNNT1的表达均显著降低。
按研究目的归类文中使用的方法,便于定位所需技术。
| 实验环节 | 名称 | 品牌 | 货号 |
|---|---|---|---|
| 蛋白质提取 | RIPA裂解液 | Elabscience Biotechnology Inc. | -- |
| 蛋白质定量 | BCA蛋白定量试剂盒 | Beijing Solabio Technology Co., Ltd. | -- |
| Western blot | PVDF膜 | -- | -- |
| 辣根过氧化物酶标记的山羊抗兔IgG二抗 | Abcam | ab6721 | |
| ECL化学发光底物试剂盒 | Shanghai Yesen Biotechnology Co., Ltd. | -- | |
| ChemiDoc成像系统 | Bio-Rad Laboratories Inc. | -- | |
| RNA提取 | TRIzol试剂 | Invitrogen | -- |
| RNA质量检测 | NanoDrop分光光度计 | Thermo Fisher Scientific | -- |
| cDNA合成 | cDNA合成试剂盒 | Takara Bio | -- |
| qPCR | 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 |