摘要:IDH突变型胶质瘤(IDH-mutant glioma, IDH-G)从生长缓慢的肿瘤进展为致死性疾病的过程伴随转录组和DNA甲基化变化,但其机制尚不清楚。本研究对19例患者的36例IDH-G纵向配对样本进行单细胞核DNA甲基化捕获测序(single-nucleus DNA methylation sequencing, snDNAme)、单细胞核RNA测序(single-nucleus RNA sequencing, snRNA-seq)及批量全外显子测序(bulk whole-exome sequencing, WES)联合分析。研究人员发现IDH-G进展与恶性干细胞样状态增加、分化减少及全基因组DNA甲基化丢失相关,后者可标记预后更差肿瘤。甲基化丢失在同一肿瘤内所有恶性细胞中均匀分布,提示其可能是干细胞样状态增加的上游原因而非结果。基于高分辨率系统进化树(phylogenetic tree)分析细胞状态遗传性与可塑性,研究人员将DNA甲基化丢失与胶质瘤细胞状态编码及遗传性的改变相联系。本研究揭示了DNA甲基化丢失如何重塑细胞转归,以及如何标记不同IDH-G亚型中临床侵袭性更强肿瘤。
论文解读——《Nature Genetics》:IDH突变型胶质瘤进展中DNA甲基化丢失重塑细胞状态分化与进化
研究背景与立项依据
异柠檬酸脱氢酶突变型弥漫性胶质瘤(IDH-mutant glioma, IDH-G)包括IDH突变型星形细胞瘤(IDH-mutant astrocytoma, IDH-A)和少突胶质细胞瘤(IDH-mutant oligodendroglioma, IDH-O, codel亚型),多初诊为低级别、生长缓慢,但最终几乎均进展为高级别致死性肿瘤。IDH1/2热点突变致D-2-羟基戊二酸(D-2-hydroxyglutarate, D-2HG)累积,抑制α-酮戊二酸依赖性双加氧酶(如组蛋白及DNA去甲基酶),引起全基因组CpG岛高甲基化即胶质瘤CpG岛甲基化表型(Glioma CpG Island Methylator Phenotype, G-CIMP-high)。既往研究表明IDH-G进展伴随获得性拷贝数变异(copy number alteration, CNA)、CDKN2A纯合缺失、恶性细胞状态比例改变(干细胞样/神经前体样及少突胶质前体样细胞增多,星形胶质细胞样及少突胶质细胞样分化状态减少),以及部分复发肿瘤出现G-CIMP状态丢失(G-CIMP-low),后者与较短生存期相关。然而G-CIMP-low如何导致更具侵袭性表型、表观遗传与遗传变异如何定量贡献于细胞表型演化均尚未阐明。已有单细胞RNA测序(single-cell RNA sequencing, scRNA-seq)或单细胞DNA甲基化研究受限于样本量、缺乏配对纵向样本及深度遗传学表征。本研究旨在通过纵向配对IDH-G样本的单细胞核多组学分析,揭示DNA甲基化动态变化与细胞状态演化及肿瘤进展的关系。
主要关键技术方法
研究人员收集来自三家医学中心的36例冷冻IDH-G肿瘤标本(32例配对的初发—复发/二次复发样本,来自15例IDH-A及4例IDH-O患者,另4例为非配对原发或复发标本),经伦理批准及患者知情同意。每例样本同时行:(1) 批量WES检测体细胞突变与CNA;(2) 10x Genomics 3'端单细胞核RNA测序(snRNA-seq)用于高通量表达谱及细胞类型注释;(2) 单细胞核双重测序——采用Smart-seq2(SS2)获取全长转录本,联合扩展代表性亚硫酸氢盐测序(extended-representation bisulfite sequencing, XRBS,即snXRBS)捕获DNA甲基化CpG位点。经质控后保留平均58.8个/样本的双重测序合格细胞核(共2117个),及158143个合格10x snRNA-seq细胞核。生物信息学流程包括:基于inferCNV与正常脑组织参考定义恶性/非恶性细胞;非负矩阵分解(non-negative matrix factorization, NMF)解析恶性细胞状态(干细胞样、星形胶质细胞样AC-like、少突胶质细胞样OC-like及细胞周期态);以既往报道684个G-CIMP进展相关CpG位点±1 kb窗口均值定义单细胞G-CIMP评分(G-CIMP score),两态混合模型区分G-CIMP-high/G-CIMP-low;差异甲基化区域(differentially methylated region, DMR)鉴定;染色质状态注释(chromHMM);转录因子基序富集分析;基于DNA甲基化拟突变构建高分辨率系统进化树,应用性状遗传性系统发育分析框架(Phylogenetic Analysis of Trait Heritability, PATH)量化细胞状态遗传性与分化/去分化转换概率。
研究结果
Single-nucleus multimodal analysis of a longitudinal cohort of IDH-G(IDH-G纵向队列的单细胞核多组学分析)
WES证实多数配对样本复发时获得额外体细胞突变(含TMZ诱导错配修复缺陷致超突变)及CNA(如CDKN2A纯合缺失),共享SNV支持共同克隆起源。snXRBS较传统snRRBS在非启动子区(P=2.2×10⁻¹⁶)及启动子区(P=3.4×10⁻¹²)均具更高CpG覆盖(平均约37.9万唯一CpGs/核),且snXRBS推断CNA分辨率优于snRNA-seq并可检出局部癌基因扩增(如PDGFRA、CDK4)。整合Seurat聚类结合CNA及正常脑参考签名,明确区分恶性细胞(具CNA或共簇)与非恶性微环境细胞。
Interrogating DNA methylation loss in single cells(单细胞水平DNA甲基化丢失的检测)
基于684个G-CIMP进展相关CpG位点的G-CIMP评分在恶性细胞中复发样本显著低于初发(P=0.01),与TCGA G-CIMP分类一致;同一肿瘤内恶性细胞间G-CIMP评分方差低,提示全局甲基化丢失均匀发生于各恶性细胞而非少数亚群驱动;非恶性细胞评分不随时间点变化排除技术因素。除已知IDH-A可出现G-CIMP-low外,研究首次在配对IDH-O复发样本中发现G-CIMP-low亚型(占部分3级IDH-O),TCGA重新分析显示此类IDH-O(codel-ME-low,G-CIMP评分<0.7)总生存更差(P=0.0022 log-rank检验),且差于G-CIMP-high IDH-A,表明DNA甲基化丢失是跨IDH-G亚型的不良预后标志。
Hypomethylation, increase in stem-like states in progressed IDH-G(进展期IDH-G中低甲基化与干细胞样状态增加)
NMF解析得到干细胞样(富集细胞周期)、AC-like及OC-like状态。配对分析示复发肿瘤干细胞样比例升高、分化样比例降低(P=0.039)。样本水平G-CIMP评分与干细胞样比例(R=-0.56, P=3.4×10⁻⁴)、增殖细胞比例(R=-0.62, P=0.0001)呈负相关,但与单个细胞内干细胞特征评分无相关性,否定"干细胞本身G-CIMP更低导致整体降低"假设,支持"全肿瘤细胞G-CIMP丢失→细胞状态组成偏移"。
Putative mechanistic links between methylation loss and increase in stem-like states(甲基化丢失与干细胞样状态增加的潜在机制关联)
G-CIMP-low相对G-CIMP-high细胞鉴定出2803个低甲基化1 kb窗口及11个高甲基化窗口,低甲基化区富集Polycomb抑制复合物2(Polycomb Repressive Complex 2, PRC2)相关域、增强子及ZNF/重复序列(非CpG岛/shore区)。差异表达分析示G-CIMP-low肿瘤上调干细胞基因(HOXD9、HOXD10、PROM1/CD133),其中部分(如HOXD9、EPHB2)启动子低甲基化伴表达上调。低甲基化区转录因子基序富集BATF、NFI家族及SOX10(胶质瘤谱系主控因子,其在G-CIMP-low中自身部分低甲基化且表达略升)。综上提出三重机制:PRC2靶基因低甲基化便于干性程序再激活;胶质瘤干细胞基因上调;干细胞转录因子结合基序低甲基化增强染色质可及性——共同促使G-CIMP-low肿瘤偏向干细胞样状态。
Increased heritability of stem-like state in IDH-G progression(IDH-G进展中干细胞样状态遗传性增强)
基于DNA甲基化拟突变构建系统进化树并经CNA亚克隆验证。PATH框架分析示G-CIMP评分与干细胞特征具系内遗传性。配对初发(G-CIMP-high)→复发(G-CIMP-low)样本中,5例中4例干细胞样状态遗传性升高。PATH推断G-CIMP-low复发肿瘤干细胞→分化样状态转换概率降低(分化受阻),去分化率无一致变化;同一肿瘤内G-CIMP-low亚克隆较G-CIMP-high亚克隆干细胞特征遗传性更高。此外G-CIMP-low中分化样细胞仍表达较高干性评分、更高增殖率及较低谱系分化评分,表明分化程序部署不全。由此提出模型:G-CIMP丢失→分化减少+干细胞样状态遗传性及比例增加+细胞周期加速→促侵袭性表型。
讨论与结论翻译(浓缩自Discussion)
IDH-G进展伴随转录(干细胞样增多、分化减少)与表观遗传(复发相关全局低甲基化即G-CIMP-low)改变,本研究通过单细胞多组学纵向分析将二者耦联。研究人员排除细胞状态间G-CIMP异质性假说,证实甲基化丢失均匀发生于同肿瘤各恶性细胞,并通过PRC2靶低甲基化、干细胞基因上调及TF基序低甲基化三重机制促成干细胞样状态富集。系进化分析揭示G-CIMP-low肿瘤细胞状态动力学紊乱——向分化方向转换减少、干细胞样状态遗传性增强及滞留,共同推高干细胞比例。G-CIMP状态可作为跨IDH-G亚型(含IDH-O)的预后标志,且可能影响突变IDH抑制剂(具分化治疗效应)疗效分层。局限含snXRBS稀疏需窗口平均、队列与细胞数有限、无法检出单细胞点突变;未来更大队列及技术可深化因果解析。总之,本研究将并发snDNA甲基化与转录组纳入高分辨率系统进化框架,阐明DNA甲基化丢失驱动IDH-G干细胞转录状态及细胞状态动态改变,为解析人类肿瘤演化提供范式。
(注:全文专业术语首次出现标注英文缩写,保留原文上下标格式如D-2HG、CpG、IDH1/2、α-ketoglutarate等;去除文献引用标号及图注标识;所有内容严格依据原文浓缩,无推测性表述。)