数据获取与差异表达分析
识别骨肉瘤组织中差异表达的失巢凋亡相关基因(ARGs)
从TCGA获取86例骨肉瘤样本,从GEO获取GSE21257(53例)和GSE33382(3正常,84肿瘤)。将794个ARGs与GSE33382差异基因取交集。
Develop a Novel Signature to Predict the Survival and Affect the Immune Microenvironment of Osteosarcoma Patients: Anoikis-Related Genes.
研究结合TCGA和GEO数据构建预后模型,并进行免疫微环境、药物敏感性分析及RT-qPCR验证,但缺乏体内功能实验。
TCGA骨肉瘤队列(86例) GEO骨肉瘤队列(GSE21257,53例) 人骨肉瘤细胞系143B和U2OS 人正常成骨细胞系hFOB1.19
训练队列:TCGA数据库86例骨肉瘤患者;验证队列:GEO GSE21257数据库53例;差异分析数据集:GSE33382(3正常,84肿瘤) 分组:按风险评分中位数划分,训练队列高风险43例,低风险43例;验证队列高风险27例,低风险26例 LASSO回归使用最优惩罚参数确定11个ARGs用于模型构建 RT-qPCR验证:使用143B和U2OS细胞系,对照为hFOB1.19,实验独立重复至少三次 统计显著性阈值P < 0.05,药物敏感性分析阈值P < 0.001
提炼研究问题、关键发现与证据,快速把握文章的核心贡献。
按研究推进顺序梳理实验设计、验证步骤与关键观察。
识别骨肉瘤组织中差异表达的失巢凋亡相关基因(ARGs)
从TCGA获取86例骨肉瘤样本,从GEO获取GSE21257(53例)和GSE33382(3正常,84肿瘤)。将794个ARGs与GSE33382差异基因取交集。
探索差异ARGs的功能通路并构建风险预后模型
对145个差异ARGs进行GO和KEGG富集分析。使用单变量Cox回归和LASSO回归确定11个与预后相关的ARGs并计算风险评分。
验证风险模型的预测能力和独立性
绘制风险曲线、生存曲线、ROC曲线,进行独立预后分析、PCA和t-SNE分析。
建立临床预测工具并探索不同风险组的通路活性
基于风险评分和临床特征构建列线图,并绘制校准曲线。使用GSEA分析高、低风险组中的活跃通路。
评估风险模型与免疫细胞浸润及免疫功能的关系
使用ESTIMATE计算基质评分和免疫评分,通过ssGSEA比较免疫细胞和功能差异,利用CIBERSORT分析22种免疫细胞相关性,并评估免疫检查点基因差异。
识别对骨肉瘤潜在有效的药物
使用pRRophetic等R包评估药物敏感性,筛选在高、低风险组间敏感性有显著差异的药物。
验证11个ARGs在骨肉瘤细胞系中的mRNA表达水平
提取143B、U2OS和hFOB1.19细胞总RNA,进行RT-qPCR检测ARGs表达,以GAPDH为内参,使用2−ΔΔCt方法计算相对表达。
按研究目的归类文中使用的方法,便于定位所需技术。
| 实验环节 | 名称 | 品牌 | 货号 |
|---|---|---|---|
| RNA提取 | TRIzol试剂 | GenStar | P118-05 |
| RT-qPCR | SYBR Green Master Mix | TOPSCIENCE | C0006 |
| 引物对 | Accurate Biology | -- | |
| 数据获取 | TCGA数据库 | -- | -- |
| GEO数据库 | -- | -- |
汇总复现实验时建议重点确认的条件及原文阅读提示。
| 环节 | 核对要点 |
|---|---|
| 数据获取和差异分析 | 训练队列:TCGA数据库86例;验证队列:GSE21257 53例;差异分析数据集:GSE33382(3正常,84肿瘤);筛选阈值P<0.05和|logFC|≥0.5 阅读提示:参见2.1节和3.1节 |
| 预后模型构建 | 单变量Cox回归P<0.01筛选12个基因;LASSO回归选取最优lambda确定11个基因;风险评分公式 阅读提示:参见2.4节和3.3节 |
| 模型验证 | 训练和验证队列的生存曲线、ROC曲线(1、3、5年)、独立预后分析(单因素和多因素Cox)、PCA和t-SNE 阅读提示:参见2.5、2.6节和3.4、3.5节 |
| 列线图 | 基于风险评分和临床特征(性别、转移状态、原发部位等)构建,评估1、3、5年生存概率,校准曲线 阅读提示:参见2.7节和3.6节 |
| 免疫微环境分析 | ESTIMATE算法计算基质、免疫和综合评分;ssGSEA比较免疫细胞和功能;CIBERSORT分析22种免疫细胞,排除P>0.05样本;免疫检查点差异分析 阅读提示:参见2.9-2.13节和3.8-3.12节 |
| 药物敏感性分析 | 使用pRRophetic等包,筛选阈值P<0.001,识别8种敏感药物 阅读提示:参见2.14节和3.13节 |
| RT-qPCR验证 | 细胞系(143B、U2OS、hFOB1.19);RNA提取和RT-qPCR条件;内参GAPDH;2−ΔΔCt方法;实验重复至少三次 阅读提示:参见2.15、2.16节和3.14节 |