《Non-Coding RNA》：A Transcriptome-Wide Analysis Nominates DIO3OS and C9orf139 as Survival-Associated Long Non-Coding RNAs in Glioblastoma
编辑推荐：
背景：长链非编码RNA（lncRNA）被广泛认为可决定胶质母细胞瘤（GBM）对替莫唑胺（TMZ）的反应，但许多提名依赖于分析者定义的生存终点、未根据O6-甲基鸟嘌呤-DNA甲基转移酶（MGMT）启动子甲基化校正、以及用筛选数据本身评估模型。研究人员追问：当终点
背景：长链非编码RNA（lncRNA）被广泛认为可决定胶质母细胞瘤（GBM）对替莫唑胺（TMZ）的反应，但许多提名依赖于分析者定义的生存终点、未根据O6-甲基鸟嘌呤-DNA甲基转移酶（MGMT）启动子甲基化校正、以及用筛选数据本身评估模型。研究人员追问：当终点定义、混杂校正与外部验证被显式处理时，哪些lncRNA–结局关联仍然成立。方法：研究人员分析来自癌症基因组图谱（TCGA）的94例IDH野生型、TMZ治疗原发性GBM患者，采用TCGA泛癌临床数据资源（TCGA-CDR）的规则化无进展间隔（PFI）与总生存（OS）作为独立终点。OS发现在TCGA与同样处理的CPTAC-3队列（n=282，209例死亡）间联合进行，校正MGMT启动子甲基化状态、细胞组成与表达亚型；CPTAC-3因无IDH注释而未按IDH选择。提名转录本在两个独立中国胶质瘤基因组图谱（CGGA）队列外部验证，增量区分度以乐观校正后C指数变化（ΔC-index）评估，并与转录组范围搜索本身产生的增益比较。结果：无lncRNA在PFI上达到转录组范围错误发现率（FDR）<0.05。联合OS发现提名DIO3OS为唯一达转录组范围显著性的转录本（每标准差风险比1.38，95%置信区间1.20–1.58，FDR 0.022）与C9orf139（FDR 0.065）；二者在各发现队列独立显著，且对MGMT状态、细胞组成与表达亚型校正稳健。但二者在两组CGGA中均为零效应且点估计反向，与PFI无关，仅增加有限乐观校正区分度（ΔC-index 0.07–0.11），接近转录组搜索自身产生的增益。复制仅见于共享处理流程的两个队列。结论：研究人员提名DIO3OS与C9orf139作为GBM中候选生存相关lncRNA，需在具备完整IDH注释与进展终点的队列中前瞻性评价，并展示终点定义、选择乐观性与外部验证如何决定哪些lncRNA关联经得起检验。
研究背景方面，胶质母细胞瘤（glioblastoma，GBM；世界卫生组织WHO 4级）是成人最常见且致死性最高的原发性恶性脑肿瘤，采用最大安全切除联合放疗与替莫唑胺（temozolomide，TMZ）的Stupp方案后中位总生存（overall survival，OS）仍仅约14.6个月。O6-甲基鸟嘌呤-DNA甲基转移酶（O6-methylguanine-DNA methyltransferase，MGMT）启动子甲基化是TMZ反应最佳确立的预测标志物，但MGMT为二分类而TMZ反应连续，仍存在大量异质性。长链非编码RNA（long non-coding RNA，lncRNA；定义为超过200核苷酸的非蛋白编码转录本）被反复报道调控GBM的TMZ敏感性，但该领域存在四类方法学弱点：分析者从不完整临床记录推导进展终点；先非校正筛选再校正使得无法知道全转录组中真正独立于MGMT的关联；在同一批筛选数据上报告模型表现而不做选择步骤的重采样校正；仅凭相关网络断言竞争内源RNA（competing endogenous RNA，ceRNA）机制而不设匹配零模型与化学计量检验。研究人员因此开展一项以终点校正、转录组范围发现、选择乐观校正与独立外部验证为核心的研究，论文发表于《Non-Coding RNA》。
主要关键技术方法方面，研究人员使用TCGA GBM的94例IDH野生型TMZ治疗原发肿瘤、CPTAC-3同流程GENCODE v36 STAR计数队列、以及CGGA mRNAseq_693与mRNAseq_325两个独立中国队列；采用TCGA-CDR规则化无进展间隔（progression-free interval，PFI）与OS双终点；全lncRNA转录组范围MGMT校正Cox回归与逆方差荟萃；bootstrap内重复候选选择以作乐观校正C-index；标志基因打分估算髓系/基质/增殖/缺氧等区室组成；Ensembl别名检索判断CGGA可量化性；ENCORI三元相关加置换匹配零模型与每细胞分子数估算检验ceRNA化学计量。
研究结果方面，2.1节终点审计显示，由GDC临床字段推导的进展变量因recurrence字段全缺而与OS数值完全相同（r=1.000），TCGA-CDR PFI则记录71事件、中位8.4月，证明必须采用规则化终点。2.2节PFI转录组范围分析显示3785个表达lncRNA中无一个达FDR<0.05，最小FDR 0.33，说明该样本量下无单个lncRNA可检出与进展相关。2.3节TCGA+CPTAC-3联合OS发现中，DIO3OS唯一达转录组范围显著（HR 1.38/SD，FDR 0.022），C9orf139次之（FDR 0.065），二者在两GDC/GENCODE队列同向独立显著；AC068888.1因MGMT校正后不显著未被推进；MROCKI仅在TCGA显著、CPTAC-3不复制。2.4节TCGA内稳健性显示，DIO3OS与C9orf139在MGMT状态、标志细胞组成、表达亚型校正下仍显著，留一患者jackknife不依赖个别样本，且不属DLK1-DIO3印记域整体信号。2.5节区分度显示，基线MGMT+年龄+性别C=0.566，加DIO3OS校正ΔC 0.071、C9orf139 0.113；但把整个选择步骤放入bootstrap后，搜索自身产生校正ΔC约0.071，说明候选增量并不高于搜索噪声。2.6节CGGA外部验证中，DIO3OS与C9orf139均为零效应且HR点估计反向（如DIO3OS 0.86、0.95；C9orf139 0.83、0.86），CGGA中DIO3OS与MEG3/MEG9共表达反而更强，说明并非注释错误，而是人群、治疗时代或共享流程导致的未明 discordance。2.7节常规非校正筛选提名16个lncRNA，在MGMT校正PFI模型下无一显著，LINC00707、ZFPM2-AS1、H19均不复制；用GDC衍生终点才会显得强，证明终点定义是假阳性来源。2.8节MGMT校正本身几乎不改变估计（HR变动≤3.1%），表达不按MGMT分层，交互未达校正显著，说明本队列中筛出候选脆弱性不是MGMT混杂造成。2.9节6月PFI组间差异表达找到ZFPM2-AS1、LINC00707等14个lncRNA，但二分类丢失时间/删失信息，仅作警示不作结论。2.10节区室打分显示ZFPM2-AS1与基质分、H19与增殖/缺氧/基质、LINC00707与髓系相关，说明候选可能来自非恶性区室，不能默认肿瘤细胞内表达。2.11节H19 ceRNA网络在置换零模型下仅边际富集（经验p=0.022），且H19约115分子/细胞，而竞争靶位点池约34000–44000转录本/细胞，占比约0.3%，两个miR（hsa-miR-519c-3p、hsa-miR-519d-3p）在TCGA-GBM近零读数，故只能说相关网络而非海绵机制。2.12节18个经典胶质瘤lncRNA中14个可检测，PFI下无一FDR<0.05，仅MALAT1名义p=0.023，敲低实验与群体预后预测不可等同。
讨论部分总结方面，研究人员强调正结果为DIO3OS与C9orf139在两家GENCODE流程队列中与较短OS相关且多种校正下稳健，但负结果同样关键：PFI无lncRNA可检出、常规筛选候选不复制、MGMT混杂假设被直检否定、H19相关网络机制薄弱。四大因素决定关联能否存活：终点定义、同一数据选择与评估、假定混杂而不检验、相关网络冒充机制。DIO3OS位于14q32的DLK1-DIO3印记区但不是全域信号；C9orf139几乎未表征；两候选都缺IDH选择性验证、都无PFI关联、都在CGGA反向归零，所以只是前瞻候选而非生物标志物。所有患者均接受TMZ，因此无法区分TMZ特异预测与总预后，无治疗对照则不能称TMZ反应标志物。局限性包括PFI事件少、TCGA-CDR不能分假性进展、CPTAC-3无MGMT/IDH/治疗注释、切除范围缺失、miRNA为前体水平、相关分析仅43例等。
结论部分翻译为：一项终点校正、转录组范围的胶质母细胞瘤分析提名DIO3OS，以及次级的C9orf139，作为在两个独立GENCODE注释队列中表达与较短总生存相关的lncRNA——TCGA队列经选择为IDH野生型，CPTAC-3队列未按IDH选择——并对MGMT启动子甲基化状态、细胞组成与转录亚型的校正稳健。两者均不与进展相关，在两个CGGA队列中均不复制，且增加的区分度不超过转录组范围搜索偶然产生的增益；因此它们应作为前瞻性评价候选而非生物标志物。在同一分析中，无任何lncRNA与无进展间隔相关，常规非校正筛选提名的转录本不能存活，而为这种脆弱性所设想的MGMT混杂解释在直接检验时未被支持。研究人员同时报告阳性与阴性发现及产生它们的框架，因为其所暴露的陷阱——退化为总生存的分析者衍生终点、在同一数据内提名与评估、从未被检验的假定混杂因子、仅由相关就断言机制——在该文献中极易出现，且标准流程输出中不可见。
