一种新型缺氧驱动基因特征可预测胶质瘤的预后和耐药性

A novel hypoxia-driven gene signature that can predict the prognosis and drug resistance of gliomas.

作者信息Peng Ren, Jing-Ya Wang, Zhi-Rui Zeng, Nan-Xi Li, Hong-Lei Chen, Xin-Ge Peng, Ujjal K Bhawal, Wen-Zhi Guo
PMID36118887
发布时间2022-09-13
DOI10.3389/fgene.2022.976356

实验完整度

包含TCGA和CGGA多队列转录组分析、60例胶质瘤组织免疫组化验证、6种胶质瘤细胞系功能实验(CCK-8、集落形成、Western blot)及siRNA敲低机制验证。

主要模型

TCGA胶质瘤队列 CGGA胶质瘤队列 60例原发和复发胶质瘤组织 胶质瘤细胞系(U87、A172等)

重点核对

TCGA和CGGA样本量及分组 缺氧特征评分方法 免疫组化抗体及评分标准 细胞系及缺氧培养条件(1% O2) TMZ浓度及处理时间(20 μM,48 h,12天)

摘要

缺氧在胶质瘤组织内部自发形成,并调节多种基因的表达。然而,胶质瘤组织中缺氧驱动基因的状态尚未完全清楚。在本研究中,将癌症基因组图谱(TCGA)中695例胶质瘤组织的RNA-seq数据作为发现队列,用于识别缺氧驱动基因并构建新型基因特征。该特征的预后价值在TCGA和中国胶质瘤基因组图谱(CGGA)的数据中得到验证。使用免疫组织化学和受试者工作特征曲线分析缺氧驱动基因的表达和诊断价值。最后,通过蛋白质印迹、CCK-8和集落形成实验分析了缺氧驱动基因对替莫唑胺(TMZ)耐药性的影响。共鉴定出169个缺氧驱动基因,这些基因与胶质瘤患者的不良预后相关。其中,22个基因的度分数≥10,选择6个基因(WT1、HOXA2、HOXC6、MMP9、SHOX2和MYOD1)构建特征,将胶质瘤患者分为低风险组或高风险组。该特征在TCGA和CGGA中对胶质瘤患者具有显著的预后价值。HOXC6、MMP9、SHOX2和MYOD1的表达与胶质瘤组织中的缺氧程度相关,且在复发案例中具有显著的诊断价值,并与胶质瘤患者的无病生存期显著相关。此外,SHOX2在O-6-甲基鸟嘌呤-DNA甲基转移酶(MGMT)未甲基化的胶质瘤组织和TMZ耐药的胶质瘤细胞系中高表达,并与MGMT表达相关。敲低SHOX2的表达显著降低了缺氧诱导的胶质瘤细胞TMZ耐药性。最终,我们确定了六个新的缺氧驱动基因用于胶质瘤的可靠预后预测,并发现SHOX2可能是克服缺氧诱导的TMZ耐药的潜在靶点。

实验结论

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

研究问题
识别胶质瘤中的缺氧驱动基因,并探索其与预后和TMZ耐药性的关系。
核心机制
SHOX2在缺氧条件下上调,可能通过抑制p53通路、上调MGMT表达,导致胶质瘤对TMZ耐药。
主要证据
SHOX2在TMZ耐药细胞系中高表达,缺氧诱导SHOX2升高,siRNA敲低SHOX2后缺氧诱导的TMZ耐药性降低;GSEA提示SHOX2可能抑制p53通路。
研究意义
该研究提供了6个缺氧驱动基因用于胶质瘤预后预测的可靠模型,并提示SHOX2可能是克服缺氧诱导的TMZ耐药的潜在靶点。

研究路径

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

1

数据获取与预处理

获取胶质瘤转录组数据和临床信息,进行标准化处理以用于后续分析。

从TCGA下载695例胶质瘤RNA-seq数据,从CGGA下载325例胶质瘤数据,进行归一化、对数转换和基因符号转换。

2

缺氧特征评分与分组

评估每个胶质瘤组织的缺氧程度,并分为高低缺氧组。

使用GSEA数据库中的50个缺氧基因集,通过t分布随机邻域嵌入方法计算缺氧特征评分,根据评分将样本分为高或低缺氧组。

3

差异表达基因分析

鉴定高低缺氧组之间的差异表达基因,作为缺氧驱动基因。

使用EdgeR包分析DEGs,通过火山图和热图展示。

4

预后相关DEGs筛选与PPI网络构建

筛选与不良预后相关的缺氧驱动基因,并构建蛋白互作网络。

通过DEGs分析比较预后好和差的胶质瘤患者,取交集,将169个缺氧驱动基因导入STRING数据库构建PPI网络,用Cytoscape可视化并计算度分数。

5

预后风险模型构建

建立基于缺氧驱动基因的风险评分模型,以预测胶质瘤患者预后。

对22个重要基因进行单因素Cox回归和LASSO Cox回归,再通过多因素Cox回归选择6个基因,构建风险评分公式。

6

风险模型验证

在TCGA和CGGA队列中验证风险模型的预后价值。

根据中位风险评分将患者分为高低风险组,进行Kaplan-Meier生存分析和ROC曲线评估。

7

临床样本验证

在胶质瘤组织中验证缺氧驱动基因表达与缺氧程度及复发的关系。

对60例胶质瘤组织进行免疫组化染色,检测Hypoxyprobe和6个基因的表达,分析其与缺氧水平及原发/复发状态的关系。

8

TMZ耐药性分析

评估缺氧驱动基因与TMZ耐药性的关联。

通过OncoPredict计算TMZ分数,分析基因表达与TMZ分数的关系;在细胞系中检测IC50,并检测SHOX2表达;缺氧处理下检测SHOX2表达,及敲低SHOX2对TMZ耐药性的影响。

研究方法

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

产品清单

实验环节名称品牌货号
缺氧探针哌莫硝唑(与Mab结合;1:200;货号HP1-100试剂盒)Hypoxyprobe, Inc.HP1-100
HOXC6抗体Thermo Fisher ScientificPA5-41479
MMP9抗体Proteintech10375-2-AP
WT1抗体AbconalA16319
SHOX2抗体Abcamab229851
MYOD1抗体Proteintech18943-1-AP
HOXA2抗体Abcamab229960
DMEM培养基----
胎牛血清----
靶向SHOX2的小干扰RNAiGene Biotechnology Co., Ltd.--
Lipo2000转染试剂Sangon--
RIPA裂解液Servicebio--
SDS-PAGE凝胶Meilune--
PVDF膜Thermo Fisher Scientific--
抗SHOX2抗体Proteintech16366-1-AP
抗ACTB抗体AbconalAC026
ECL化学发光试剂Proteintech--
CCK-8试剂----
结晶紫----
三气培养箱Thermo Fisher Scientific--

关键环节

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

环节核对要点
数据获取与预处理
TCGA和CGGA样本量、归一化方法、R软件版本
阅读提示:Materials and methods: Data acquisition and preprocessing
缺氧特征评分
50个缺氧基因集、t-SNE方法、分组阈值
阅读提示:Materials and methods: Identification of hypoxia characteristics in glioma tissues
差异表达分析
DEGs阈值 |LogFC|>1, adj p<0.05
阅读提示:Materials and methods: Analysis of DEGs
PPI网络
STRING数据库、Cytoscape、度分数阈值≥10
阅读提示:Materials and methods: Protein-protein interaction network
预后模型构建
单因素Cox p<0.05,LASSO交叉验证,多因素Cox所选基因
阅读提示:Materials and methods: Prognostic gene screening and prognostic signature
模型验证
中位风险评分分组、Kaplan-Meier、ROC曲线
阅读提示:Results: Prognostic risk model validation
免疫组化
抗体稀释比、抗原修复条件、DAB染色条件
阅读提示:Materials and methods: Immunohistochemistry staining
细胞培养与缺氧模拟
细胞系来源、培养基、缺氧条件(1% O2)、siRNA序列及转染试剂
阅读提示:Materials and methods: Cell culture and transfection
TMZ敏感性检测
细胞密度、TMZ浓度、处理时间、CCK-8试剂孵育时间、IC50计算
阅读提示:Materials and methods: Cell count kit-8, Colony formation assay