胶质瘤细胞系中m6A调节因子扰动后的长读长转录组分析

《Biochemistry and Biophysics Reports》：Long-read transcriptomic profiling of a glioma cell line following m6A regulator perturbations

　　

摘要：N6-甲基腺苷(m6A)是最普遍的内部mRNA修饰，且在中枢神经系统(CNS)中富集，然而其在胶质瘤中的作用仍未完全明确。研究人员利用长读长直接RNA测序(direct RNA sequencing)，在靶向敲低m6A阅读器(reader) IGF2BP2、写入器(writer) METTL3和擦除器(eraser) ALKBH5后，绘制了单个胶质瘤细胞系的转录组范围m6A修饰图谱。在多种扰动条件下，m6A的整体结构（包括转录本类别、位置富集和位点多重性）基本得以保留，而差异甲基化与基因表达的关联较弱。相比之下，m6A调节因子的扰动与广泛的异构体转换(isoform switching)以及非翻译区(UTR)、编码潜力和预测的转录本命运的变化相一致，且这些变化在很大程度上独立于批量基因表达变化。公共胶质瘤数据集进一步支持了胶质瘤中异构体特异性变化的相关性。总之，这些发现强调了与m6A调节因子扰动相关的异构体水平转录变异，并支持使用长读长、可解析异构体的方法来研究胶质瘤中的RNA调控状态。

论文解读：胶质瘤中m6A调节因子扰动引发的异构体水平重塑

研究背景与意义

超过163种RNA修饰已被鉴定，其中N6-甲基腺苷(m6A)是mRNA中最普遍的内部修饰，尤其在中枢神经系统(CNS)中含量丰富。m6A的动态平衡由写入器(甲基转移酶，如METTL3)、阅读器(结合蛋白，如IGF2BP2)和擦除器(去甲基化酶，如ALKBH5)共同维持，参与调控RNA的可变剪接、稳定性及翻译等代谢过程。胶质瘤是最常见的原发性CNS肿瘤，胶质母细胞瘤(GBM)作为其最恶性亚型，患者预后极差。尽管已有研究提示m6A与胶质瘤的进展、干性及微环境有关，但既往研究多依赖短读长测序，无法解析转录本异构体结构，且m6A在胶质瘤中是通过全局转录控制还是选择性异构体调控发挥作用仍不清楚。为此，研究人员采用长读长直接RNA测序技术，旨在从异构体层面揭示m6A调节因子扰动对胶质瘤转录组的影响。该研究成果发表在《Biochemistry and Biophysics Reports》。

关键技术方法

本研究主要采用牛津纳米孔(Oxford Nanopore)直接RNA测序技术对人胶质瘤细胞系Gli36进行转录组及m6A位点单核苷酸分辨率图谱绘制。通过siRNA介导的基因沉默技术分别敲低IGF2BP2、METTL3和ALKBH5，并设置未处理(Na?ve)和乱序siRNA(Scrambled)对照，每组包含三个技术重复。数据分析方面，利用m6Anet机器学习工具识别m6A修饰位点，使用IsoformSwitchAnalyzeR进行差异异构体使用分析，并通过定量PCR(qPCR)对特定异构体进行独立验证。此外，研究还利用了GlioVis平台和CanIsoNet数据库对临床胶质瘤样本中的m6A调节因子表达及异构体模式进行了外部验证。

研究结果

1. 胶质瘤细胞中m6A调节因子扰动的质量评估与转录组图谱分析

研究人员成功构建了靶向敲低IGF2BP2、METTL3和ALKBH5的体外胶质瘤模型。通过对15个纳米孔测序运行数据的质控分析显示，所有条件下蛋白质编码转录本占主导地位。主成分分析(PCA)表明各条件间具有可重复的条件水平结构。虽然测序产量存在差异，但全球RNA生物类型组成在各扰动条件下基本保留，为下游分析奠定了坚实基础。

2. 原始胶质瘤细胞及m6A调节因子敲低中的转录组范围m6A图谱绘制

研究人员量化了各条件下的m6A修饰情况。结果显示，原始(Na?ve)胶质瘤细胞显示出最高的m6A修饰负担（30.2%的基因被修饰）。在所有敲低(KD)和乱序对照中，检测到的m6A相对于原始细胞均有所减少，其中ALKBH5 KD的m6A修饰最低。尽管m6A丰度下降，但其分布的RNA生物类型（主要是蛋白编码转录本）在各条件下是保守的。位置分布分析显示，原始细胞和对照组在CDS区域和3'UTR近端具有较高的m6A密度，而敲低条件则表现出向3'UTR偏移的趋势，其中ALKBH5 KD细胞的位移最为显著。

3. m6A调节因子敲低中m6A甲基化与基因表达的整合分析

研究人员聚焦于跨原始及敲低条件共有的578个转录本进行深入分析。结果显示，大多数共有转录本在敲低条件下相对于原始细胞呈低甲基化状态。值得注意的是，m6A甲基化的变化与基因表达水平的耦合度较弱。例如，ALKBH5 KD与原始细胞相比，尽管转录本甲基化水平显著降低，但基因表达谱的变化并不总是一致。功能通路分析显示，RNA结合通路（如PUF60、DDX39A）在调节因子扰动下表现出特定的甲基化和表达模式变化，但许多差异表达基因并未伴随可检测到的甲基化改变。

4. 原始胶质瘤细胞及调节因子敲低的独特m6A甲基化图谱

分析各条件独有的m6A位点发现，原始细胞拥有最多数量的独特修饰位点、转录本和基因。基序(k-mer)分析揭示了条件特异性的变化，例如GGACT在原始细胞中富集，而在敲低条件下GAACT更为普遍。区域映射显示，所有敲低条件均表现出CDS甲基化减少和3'UTR甲基化增加的趋势。相关性分析表明，独特修饰转录本的甲基化比率与基因表达呈弱正相关，这种关联在METTL3 KD中最为显著。

5. m6A调节因子敲低引发的异构体转换及结构后果

这是本研究最显著的发现之一。异构体转换分析显示，调节因子扰动导致了广泛的异构体使用变化。其中ALKBH5 KD的异构体开关数量最多（1364个），其次是IGF2BP2 KD（821个）和METTL3 KD（459个）。这些异构体层面的变化往往比基因水平的变化更明显。结构后果分析预测，这些转换常伴随着开放阅读框(ORF)丢失、结构域丢失和外显子丢失等特征。特别是ALKBH5 KD还富集了无义介导的衰变(NMD)敏感异构体。qPCR实验验证了如MET和BIRC3等基因的特定异构体在敲低条件下确实发生了显著的使用频率改变，且常与预测的NMD敏感性或外显子结构变化相关。

6. m6A调节因子、剪接因子和衰变机制的动态表达

研究人员观察到m6A调节因子及其异构体在敲低条件下的差异化表达。例如，ALKBH5-201在ALKBH5 KD中显著降低。同时，RNA剪接和衰变相关因子（如CNOT1、DCP2、EXOSC3）的表达也发生了系统性变化。这表明m6A调节因子的扰动并非孤立事件，而是引发了更广泛的RNA加工和衰变网络的协同重塑。

7. 胶质瘤中m6A调节因子表达和异构体使用的跨数据集验证

利用GlioVis平台分析公共数据发现，ALKBH5和IGF2BP2在胶质母细胞瘤(GBM)中表达最高，且与较短的总生存期相关，这与细胞系实验结果一致。通过与CanIsoNet数据库的重叠分析，研究人员在细胞系模型中发现了多个与患者胶质瘤相关的主导异构体切换事件（如TBC1D14、ZNF585A等），证实了体外模型发现的生物学相关性。

结论与讨论

本研究表明，在胶质瘤细胞模型中扰动核心m6A调节因子后，m6A的整体架构（包括转录本类型、位置分布和位点多重性）基本得以保留。然而，最突出的变化发生在转录本异构体层面，表现为广泛的异构体转换(isoform switching)，且这些转录本层面的重塑往往独立于批量基因表达变化。这解释了为何传统基于基因水平的分析可能遗漏了关键的调控信息。

研究结论强调，m6A调节因子扰动与预测的转录本结构变化（如UTR长度、编码潜力和结构域组成）密切相关，尤其是在ALKBH5敲低条件下观察到的NMD敏感异构体增加，具有重要的生物学意义。此外，m6A甲基化变化与基因表达之间仅存在微弱关联，提示m6A在胶质瘤中可能主要通过影响RNA的命运（如稳定性、定位和翻译效率）而非仅仅改变转录本总量来发挥作用。

该研究通过整合长读长测序技术，揭示了胶质瘤中m6A调控的异构体特异性机制，强调了在癌症研究中采用异构体解析方法进行转录组分析的必要性，为理解胶质瘤的表观转录组复杂性提供了新的视角。