转录组学(meta-analysis)与功能验证鉴定长链非编码RNA(lncRNA)作为寨卡病毒(Zika virus, ZIKV)介导胶质母细胞瘤细胞系溶瘤作用的调节因子

《Cells》：Transcriptomic Meta-Analysis and Functional Validation Identify Long Non-Coding RNAs as Modulators of Zika Virus-Mediated Oncolysis in Glioblastoma Multiforme Cell Lines

　　

摘要：胶质母细胞瘤(glioblastoma multiforme, GBM)是最具侵袭性的原发性脑部恶性肿瘤，治疗选择有限且临床预后差。由于寨卡病毒(Zika virus, ZIKV)可选择性感染并杀伤快速增殖的神经细胞，将其用于GBM治疗日益受到关注。多项研究重复证实ZIKV可有效杀伤GBM细胞。研究人员旨在通过对已发表ZIKV杀伤GBM细胞的转录组研究进行meta?analysis，揭示该细胞毒效应的分子机制。研究人员整合了4套GBM数据集，并以神经母细胞瘤(neuroblastoma, NBM)数据集作为外群对照。分析鉴定出ZIKV感染的GBM细胞共有分子特征，有趣的是ZIKV杀伤的GBM细胞中TNF、NF?κB及p53信号通路等常见于增殖与转移调控的通路发生失调。采用无假设驱动设计，研究人员发现数种在ZIKV感染的GBM中一致失调的lncRNA，其中多数此前未被认识参与癌细胞死亡。在此类lncRNA中，研究人员验证了4种lncRNA在ZIKV介导溶瘤中的作用。研究人员利用siRNA敲减技术在成人GBM细胞系中对MELTF?AS1、TIPARP?AS1、NR2F1?AS1及SLC9A3?AS1进行功能检测：沉默MELTF?AS1增强ZIKV诱导的细胞死亡，而敲减TIPARP?AS1、NR2F1?AS1及SLC9A3?AS1减弱溶瘤效应，表明这些lncRNA的调控与ZIKV细胞毒性的改变相关。上述发现阐明了GBM细胞系中ZIKV溶瘤的候选机制，突出了新型lncRNA靶点，并为进一步探索以lncRNA调控增强GBM及相关恶性肿瘤溶瘤病毒疗法提供了依据。

论文解读：转录组学Meta?Analysis与功能验证鉴定长链非编码RNA作为寨卡病毒介导胶质母细胞瘤溶瘤作用的调节因子（发表于《Cells》）

一、研究背景与立题依据

胶质母细胞瘤(glioblastoma multiforme, GBM；WHO 4级，IDH?野生型成人弥漫型胶质瘤)是最常见的恶性脑肿瘤，标准治疗（最大范围手术切除＋放疗＋替莫唑胺化疗）后仍极易复发，患者中位总生存期仅约15个月，亟需创新治疗策略。溶瘤病毒(oncolytic virus)可选择性与快速增殖的肿瘤细胞结合并裂解之，同时相对 spared静止正常组织。寨卡病毒(Zika virus, ZIKV)为嗜神经病毒，其优先感染并在增殖神经细胞中复制的特性（与胎儿小头畸形相关），被提出可用于靶向杀伤高表达SOX2?整合素αvβ5及MSI1的GBM肿瘤干细胞。已有独立研究报道ZIKV在GBM及神经母细胞瘤(neuroblastoma, NBM)模型中引起细胞死亡并富集I型干扰素(interferon, IFN)信号，但均聚焦蛋白编码基因，ZIKV感染中lncRNA（长度＞200 nt且无蛋白编码潜能的非编码RNA）的表达变化及其对溶瘤过程的调节作用尚未系统研究。鉴于lncRNA可作为癌基因或抑癌基因并通过ceRNA等机制参与GBM恶性演进，研究人员假设ZIKV感染GBM后特定lncRNA发生一致性失调并参与调控溶瘤敏感性，因此开展本研究。

二、主要关键技术方法概述

研究人员检索GEO数据库筛选符合条件（人源、RNA?seq、ZIKV感染vs_mock对照、≥2生物学重复）的数据集，最终纳入2套GBM数据集（Zhu et al. 2017, Dakar株；Bulstrode et al. 2022, PE243株）与2套NBM数据集（Bonenfant et al. 2020, MR766及PRVABC59株）作为外群对照以滤除泛神经细胞共有响应。使用RefSeq（过滤＜250 bp转录本）加聚类代表性ZIKV基因组为参考，经Trimmomatic质控、bowtie2比对、salmon定量，DESeq2行差异表达分析(FDR＜0.05)，区分编码与非编码基因。对GBM共有差异表达lncRNA查阅GeneRIF注释后遴选候选靶标。体外采用U87（PTEN突变/MGMT启动子甲基化）与A172（EGFR突变/PTEN缺失）成人GBM细胞系，进行siRNA介导lncRNA敲减，单独或与MOI＝3 ZIKV（法属波利尼西亚株）共处理，CellTiter?Glo检测细胞活力，qRT?PCR验证敲减效率。

三、研究结果

3.1. ZIKV?Infected GBM and NBM Show Distinct Gene Expression Profiles（ZIKV感染的GBM与NBM呈现不同基因表达谱）

对批次校正后表达矩阵作PCA显示ZIKV感染组与对照组在GBM及NBM中均清晰分离（编码与非编码基因均如此）。Venn图比较发现GBM两数据集间差异表达基因重叠高于跨肿瘤类型重叠，且无一非编码转录本在全部4个数据集共有差异表达，说明ZIKV转录响应具细胞类型特异性，故聚焦GBM内部共有变化。

3.2. Consensus Protein?Coding Expression Profile in GBM（GBM中共有蛋白编码基因表达谱）

两个GBM数据集meta?合并鉴定出294个上调与195个下调蛋白编码基因。KEGG富集显示上调基因显著富集于TNF(p＝2.9×10??)、NF?κB(p＝1.3×10??)及p53(p＝3.8×10??)信号通路；DrugBank提示与具有抗GBM活性的andrographolide靶标重叠；转录因子富集分析显示NF?κB及CREB结合位点显著富集，提示ZIKV触发促炎与抗增殖转录重编程。

3.3. Differential Expression of Long Non?Coding RNAs（长链非编码RNA的差异表达）

限定GBM两数据集共有差异表达非编码转录本，筛选出31个lncRNA（15个上调、16个下调），均无NBM共有失调，排除未表征者后获12个具文献记载功能的lncRNA，最终选取MELTF?AS1（ZIKV感染后下调，文献报道具促生存癌基因作用）、TIPARP?AS1与NR2F1?AS1（均上调，分别报道可抑制PARP?7及具抑癌/促凋亡作用）和SLC9A3?AS1（下调，文献报道在其他癌中促癌）进行功能验证。

3.4. Functional Validation of lncRNA Candidates（lncRNA候选靶点的功能验证）

在U87与A172细胞中进行siRNA敲减±ZIKV感染：

• •
  沉默MELTF?AS1单独降低细胞活力，联用ZIKV进一步降低活力，表明其具促生存作用，其下调有助于ZIKV溶瘤。
• •
  沉默TIPARP?AS1在A172细胞中（尤其120 h）显著削弱ZIKV引起的活力下降（存活率升至314%对照），U87影响较小，提示TIPARP?AS1可能通过抑制TIPARP/PARP?7解除I型IFN信号抑制而促进溶瘤，缺失削弱杀瘤效果且存细胞株依赖性。
• •
  沉默NR2F1?AS1部分削弱A172中ZIKV溶瘤（96 h存活率121%对照），U87反应相反且弱，提示其在部分GBM背景中辅助溶瘤程序。
• •
  沉默SLC9A3?AS1单独轻微降活力，但联用ZIKV反而升高A172与U87存活率（A172达303%对照），表明尽管在其他癌中为癌基因，在本GBM?ZIKV背景下它促进溶瘤，功能具组织/情境依赖性。

四、讨论与结论总结

研究人员通过转录组meta?analysis结合功能实验证明：ZIKV感染GBM引发TNF/NF?κB/p53通路激活及特定lncRNA一致性失调；其中MELTF?AS1为ZIKV溶瘤的负调因子（其敲减增强杀瘤），TIPARP?AS1、NR2F1?AS1与SLC9A3?AS1为ZIKV溶瘤的正调/辅助因子（其敲减减弱杀瘤），后者在GBM?ZIKV微环境中可表现与它在其他肿瘤中相反的生物学角色。这提示lncRNA不仅被动响应感染，还主动塑造利于病毒溶瘤的分子环境（如TIPARP?AS1?PARP?7?I型IFN轴），与蛋白编码基因共同形成"分子钳"——促死通路激活配合促生存lncRNA抑制及促溶瘤lncRNA诱导。研究局限含公共数据集少、仅用两株细胞系、未在原代细胞/体内证实因果机制，NBM仅作外群过滤而非直接类比GBM。临床转化需采用减毒/ miRNA?调控重组ZIKV并评估对非恶性星形胶质细胞影响。结论如下：

Conclusions（研究结论）：

本研究证明lncRNA表达是ZIKV介导GBM溶瘤的功能相关组分，并提出lncRNA靶向策略可优化ZIKV溶瘤治疗。转录组meta?analysis鉴定出ZIKV感染GBM中TNF、NF?κB及p53信号通路的一致失调，以及若干可调节ZIKV溶瘤效率的候选lncRNA。功能验证表明MELTF?AS1为促生存癌因子，其沉默增强ZIKV诱导细胞死亡；TIPARP?AS1、NR2F1?AS1与SLC9A3?AS1为候选促溶瘤因子，其敲减关联细胞系模型中ZIKV杀瘤效应减弱。上述发现支持进一步探究靶向lncRNA调控以增强ZIKV溶瘤效率及克服GBM治疗抵抗。