《Discover Oncology》:Integrative analysis of mitochondrial-related gene profiling identifies prognostic clusters and drug resistance mechanisms in low-grade glioma
研究人员在本次研究中主要运用了以下几种关键技术方法:首先是多组学数据处理技术,从癌症基因组图谱(The Cancer Genome Atlas,TCGA)数据库获取 LGG 的 RNA 测序表达数据和临床信息,并将原始计数数据转换为每百万转录本(Transcripts Per Million,TPM),进行 log2(x + 1) 转换;同时还选取了多个公开数据集,如中国胶质瘤基因组图谱(Chinese Glioma Genome Atlas,CGGA)等进行验证。其次是共识聚类分析技术,利用 R 包 ConsensusClusterPlus 对患者进行分子亚型分类,通过多次迭代和样本子采样确定最佳聚类数。另外,还运用了突变景观分析、免疫浸润分析、基因富集分析等技术,全面探究线粒体相关基因与肿瘤各方面特征的关系。
3.1 基于线粒体相关基因的共识聚类
研究人员对 TCGA 数据集中线粒体相关基因的表达水平进行共识聚类分析,最终将患者分为三个不同的簇(Cluster A、B 和 C)。通过单样本基因集富集分析(single-sample Gene Set Enrichment Analysis,ssGSEA)计算线粒体相关通路的得分,发现簇 B 患者的生存情况最差。在多个独立验证数据集上进行同样的分析,结果都显示簇 B 患者预后较差,并且还发现簇 B 中存在一些高表达的基因,进一步证实了其独特的遗传表达特征和不良预后。
3.2 突变景观分析
对三个簇的突变数据进行深入分析,绘制瀑布图展示各簇中常见的突变基因和突变类型。通过对四个最常突变基因的差异分析,发现不同簇之间这些基因的突变频率存在显著差异。研究还分析了突变基因的共发生和相互排斥事件,构建了共突变热图,同时展示了各簇的拷贝数变异(Copy Number Variation,CNV)景观,为理解肿瘤的发生发展机制提供了重要线索。
3.3 免疫浸润分析
研究人员分析了八个免疫检查点基因在不同簇中的表达水平,发现簇 B 中这些基因的表达水平显著高于其他簇,表明其免疫抑制机制更为活跃。通过计算免疫评分、基质评分和 ESTIMATE 评分,发现簇 B 的免疫评分较高,但基质评分较低。进一步对各种免疫细胞类型进行分析,发现簇 B 中大多数免疫细胞类型的浸润程度更高,且该簇在免疫耗竭和免疫抑制相关通路中的得分较高,说明其免疫微环境具有较强的免疫逃逸能力。
通过对 TCGA 和 CGGA 数据集进行基于表达的时间序列分析,研究人员发现簇 4(在两个数据集中均为高表达组)的基因表达模式与较差的预后相关。对两个数据集的簇 4 基因集进行交集分析,确定了一些在两个数据集中均上调的基因,这些基因可能是潜在的线粒体相关预后标志物。通过跨数据集验证,进一步证实了这些基因的表达和预后意义。基于这些基因计算 MITscore,将样本分为高 MIT 和低 MIT 组,发现高 MIT 组患者预后较差,并且差异表达基因与免疫反应、细胞增殖、代谢和凋亡等生物过程相关。
3.5 与肿瘤干性特征的比较
研究人员从 MSigDB 数据库中提取肿瘤干性相关的特征基因集,计算每个样本的干性得分。结果发现簇 B 和高 MIT 组的干性得分显著更高,表明线粒体功能与预后较差患者的干性特征呈正相关。通过对线粒体标志物和干性标志物的交集分析,发现簇 B 中这些交集基因的表达水平更高,说明线粒体相关标志物与干性标志物存在重叠,线粒体功能可能通过调节肿瘤细胞的干性影响肿瘤的侵袭性、治疗抗性和预后。
3.6 化疗药物敏感性和免疫治疗分析
研究人员计算了 MIT 标记物得分与化疗药物敏感性之间的相关性,发现一些药物的敏感性与 MIT 标记物得分显著相关。选取六种高 MIT 得分组更敏感的药物进行分析,发现高 MIT 组患者对这些药物的半最大抑制浓度(IC50)值显著更低,表明他们对这些化疗药物更敏感。在免疫治疗方面,通过对 IMvigor210 和 GSE78220 免疫治疗队列的分析,发现高 MIT 得分患者在免疫治疗后的生存期显著更长,且高 MIT 得分组的肿瘤新抗原负荷(Tumor Neoantigen Burden,TNB)值更高,说明他们可能具有更强的免疫原性,对免疫治疗的反应更好。
3.7 利用 SRGA 鉴定新的 MIT 标记物
利用单区域基因关联(Single-Region Gene Association,SRGA)分析包,研究人员对 MIT 标记物与其他基因标记物进行相关性分析,确定了一些与 MIT 标记物高度相关的基因,如 KDR、LRRK2、SQSTM1、MMP9 和 ARMCX2 等。构建网络图谱展示这些基因与 MIT 标记物之间的相互关系,发现它们可能通过共同的生物学途径调节肿瘤的生物学行为。综合评估这些基因的生物学意义和临床相关性,确定它们为潜在的新 MIT 标记物,有望用于临床肿瘤的诊断和治疗监测。