本文要点
- m1A 单碱基定位:TGIRT 逆转录突变谱(A-to-T/C/G 特征峰)+ AlkB/ALKBH3 去甲基化阴性对照,区分 SNP 与二级结构假象
- m1A 的 Dimroth 重排陷阱:pH>8 或加热会转成假 m6A,样本全程保持 pH 6.0–7.2 低温
- m5C 经典方案 RNA-BS-seq 损伤大(高温强酸致 RNA 降解);Aza-IP、miCLIP、DRAM-seq 实现无损单碱基解析
- m7G 用 NaBH4 还原+苯胺裂解(TRAC-seq/AlkAniline-seq)单碱基定位,前置 DXO/TAP 脱帽避免帽端信号压制
- 纳米孔直接 RNA 测序(DRS)保留长程多修饰共存信息,是单分子检测前沿
- m1A 单碱基定位:TGIRT 逆转录突变谱(A-to-T/C/G 特征峰)+ AlkB/ALKBH3 去甲基化阴性对照,区分 SNP 与二级结构假象
- m1A 的 Dimroth 重排陷阱:pH>8 或加热会转成假 m6A,样本全程保持 pH 6.0–7.2 低温
- m5C 经典方案 RNA-BS-seq 损伤大(高温强酸致 RNA 降解);Aza-IP、miCLIP、DRAM-seq 实现无损单碱基解析
- m7G 用 NaBH4 还原+苯胺裂解(TRAC-seq/AlkAniline-seq)单碱基定位,前置 DXO/TAP 脱帽避免帽端信号压制
- 纳米孔直接 RNA 测序(DRS)保留长程多修饰共存信息,是单分子检测前沿
序幕:后 m6A 时代的表观转录组版图重构
在真核生物转录组多样性的化学图景中,核糖核酸链条上的共价修饰构成了调控转录本命运与生物学功能的核心轴心。围绕 N6-甲基腺苷(m6A)建立的甲基转移酶、去甲基化酶与结合蛋白调控范式不断扩展,推动表观转录组学研究视野迅速延伸至其他具备独特物理化学特征的核苷修饰位点。分子结构层面的细微甲基化取代,赋予了不同核苷酸截然不同的生化行为与空间构象。腺嘌呤第一位氮原子发生的甲基化(m1A),在生理 pH 介质中携带局域正电荷,通过阻断经典 Watson-Crick 氢键配对边缘,对RNA局部双螺旋堆积以及核糖体翻译延伸造成显著的空间位阻与静电排斥。胞嘧啶第五位碳原子的甲基化(m5C),在维持 Watson-Crick 配对几何构象的同时,改变了碱基堆积自由能,调节局部二级结构刚性,并特异性招募细胞核输出适配体与转录本稳定调控因子。鸟嘌呤第七位氮原子的甲基化(m7G),既占据真核信使RNA(mRNA)5′端封闭帽子结构的起始位点,又以内部转录后修饰的形式分布于非编码小RNA及特定编码转录本中,其携带的离域正电荷使嘌呤环对亲核开环反应呈现出高度的化学敏感性。解析这些修饰在全转录组范围内的精确丰度与化学定位,构成了揭示转录后分子调控网络的前提。早期基于抗体交联富集的方法虽拓展了修饰图谱的宏观视界,却常受限于非特异性交叉反应与抗体表位识别宽泛度,引发了关于位点丰度与计量学的持续学术探讨。近期基于特异化学反应、逆转录酶通读阻断突变谱学以及内切酶定向裂解等技术的交融突破,正推动 m1A、m5C 与 m7G 向着单碱基分辨率与绝对化学计量学的精确高度迈进。
m1A 检测方案:空间电荷阻障、酶促校正与单碱基突变特征谱
生物学调控轴线与拓扑学争议
哺乳动物细胞体系中,tRNA 甲基转移酶复合体 TRMT6/TRMT61A 催化特定 tRNA 发卡结构中第58位腺嘌呤发生甲基化修饰,该复合体亦能靶向细胞质部分 mRNA 中呈现 GUUCRA 保守基序的茎环结构并沉积 m1A。在线粒体基质内部,线粒体专属的 TRMT10C 与 TRMT61B 则分别对线粒体编码的 ND5 等转录本引入位点特异性 m1A 修饰。在去甲基化机制维度,AlkB 家族同源蛋白 ALKBH1 与 ALKBH3 表现出针对单链 RNA 上 m1A 去甲基化的铁依赖性双加氧酶催化活性。围绕细胞质 mRNA 内部 m1A 丰度与化学计量,学术界曾出现过深刻的技术辨析。早期抗体富集联合测序报告了转录组中数千个高丰度内部位点。然而后续利用单碱基分辨生化图谱与去甲基化酶平行校正的研究澄清了相关认识:细胞质 mRNA 内部仅包含极少数由 TRMT6/TRMT61A 低效催化的超低计量修饰位点,早期海量富集峰主要源自抗体对 mRNA 5′端帽子附近 cap+1 位点甲基化修饰的非特异吸附、线粒体转录本错配比对以及高丰度 tRNA 片段污染。
逆转录错配突变谱与延伸截断特征
化学结构上,腺嘌呤 N1 位的甲基化占据了经典碱基配对的氢键供体位点,逆转录酶在合成互补脱氧核糖核酸(cDNA)过程中遭遇该位点时,遭遇明显的空间碰撞障碍。这种酶学受阻现象产生两种互为表里的测序读段特征,一方面表现为逆转录延伸提前截断(RT truncation/arrest),另一方面表现为强行跨越位点时的非同源碱基错误掺入(RT misincorporation)。常规禽成髓细胞瘤病毒(AMV)或莫洛尼氏鼠白血病病毒(M-MLV)逆转录酶在跨越 m1A 时极易发生解离,产生大量非特异截断末端,易受 RNA 自身复杂二级结构干扰。为了获得信噪比更高的单碱基定位信号,采用兼具高保真通读能力与模板跳跃活性的逆转录酶(如嗜热栖热菌内含子编码的 TGIRT 酶)展开逆转录反应,显著提高了跨越修饰位点的通读比例。在测序文库中,TGIRT 跨越 m1A 位点时主要诱导形成 A-to-T、A-to-C 或 A-to-G 的特征突变突发峰,突变频谱展现出独特的化学指纹。借助算法从错配率与突变类型分布出发,研究者可以直接推算潜在的修饰坐标。
酶促去甲基化阴性对照的设计逻辑
单纯依赖逆转录突变谱仍面临两大生信假阳性挑战:基因组内部的单核苷酸多态性(SNP)位点,以及复杂 RNA 高级结构折叠引起的酶促停滞假象。确立 m1A 真实性的核心实验设计,在于引入高活性重组去甲基化酶的平行对照反应。在体外反应体系中,加入经工程化改造的大肠杆菌 AlkB 蛋白或人类重组 ALKBH3,在 Fe2+ 与 2-氧戊二酸(2-OG)共存条件下,选择性氧化切除 N1 位的甲基基团,将 m1A 转化为未修饰腺嘌呤。由于 AlkB 催化过程未破坏天然磷酸二酯骨架与常规腺苷序列,真实存在的 m1A 位点经 AlkB 处理后,其逆转录截断信号与特定突变比例发生明显坍缩;相反,SNP 位点或序列本征结构造成的延伸停顿,在 AlkB 处理前后保持稳定的突变频率。这一双轨消除机制构成了确证全转录组低丰度 m1A 位点不可或缺的技术基准。
化学重排防护:Dimroth 重排陷阱与质控策略
由于 m1A 具备独特的化学活性,样品制备过程需严格规避碱性诱导的 Dimroth 重排反应。在 pH > 8.0 伴随加热的环境中,腺嘌呤嘧啶环的 C2 位易遭氢氧根亲核进攻,引发嘌呤环局部开环与氨基旋转,最终自发异构化重排为热力学更为稳定的 N6-甲基腺苷(m6A)。发生非预期的 Dimroth 重排后,不仅原本的 m1A 信号在逆转录中彻底隐匿,还会人为抬高下游 m6A 谱学的背景噪音。因此,全套裂解、片段化与沉淀操作均需维持在严格的中性至微弱酸性缓冲液体系(pH 6.0–7.2)及低温环境中进行。
m5C 检测方案:重亚硫酸盐脱氨机制与免转化无损策略
生物学调控轴线与细胞功能
哺乳动物 RNA 胞嘧啶 C5 位甲基化主要由 NOL1/NOP2/SUN 结构域(NSUN)家族甲基转移酶(如 NSUN2、NSUN6)以及 DNA 甲基转移酶同源物 TRDMT1(DNMT2)催化介导。在细胞代谢网络中,NSUN2 针对 tRNA 可变环及反密码子环特定胞嘧啶(如 C34、C48)沉积 m5C,维系其构象完整性并防御血管生成素(angiogenin)引发的异常位点切割。在信使核糖核酸层面,细胞质 mRNA 内部由 NSUN2 引入的 m5C 位点能够被特异性核转运适配蛋白 ALYREF 直接结合,加速成熟 mRNA 的核质转运进程。转录本 3′UTR 区域的 m5C 修饰亦可招募冷休克结构域蛋白 YBX1,借由稳定 mRNA 链条参与恶性肿瘤的浸润转移调控。
重亚硫酸盐转录组测序(RNA-BS-seq)的化学反应机制与技术短板
RNA 重亚硫酸盐转化测序(RNA-BS-seq)代表了长久以来解析 m5C 单碱基分布的经典实验路径。其化学反应依托单链 RNA 中胞嘧啶对于重亚硫酸根离子的三步亲核加成历程:Cytidine + HSO3− ⇌ Cytidine-6-sulfonate。该中间体随后脱氨基形成尿嘧啶-6-磺酸盐,并在弱碱性温育下水解脱除磺酸基团,转化为尿嘧啶(在后续逆转录产物中呈现为胸腺嘧啶 T)。相比之下,由于胞嘧啶 C5 位甲基基团的存在,增强了环系给电子共轭效应,并对 C6 位的亲核加成构成强烈的立体位阻,致使 m5C 对重亚硫酸盐诱导的脱氨作用具备显著抵抗性,在逆转录建库与测序输出端依然保留原始胞嘧啶(C)的读段特征。尽管 RNA-BS-seq 赋予了位点定量的直接途径,但其固有的理化局限性引发了持续的分析争议。首先,长时间处于高温(60–70°C)及高浓度亚硫酸氢盐(> 3.5 M)环境促使 RNA 骨架中的核糖-磷酸二酯键频繁自发水解,造成严重的模板降解,要求起始材料通常多达数微克纯化 Poly(A) RNA 才能满足建库需求。其次,RNA 分子自发折叠成的稳定发卡或茎环双链结构显著抑制了亚硫酸氢根对 C6 位的亲核进攻,容易造成未完全脱氨的假阳性位点。即便采用二甲基亚砜(DMSO)进行变性预热,结构复杂区域的残留未转化背景仍对信噪比构成干扰。最后,全转录组几乎全部非修饰 C 均转变为 U,导致 cDNA 测序读段退化为近似三碱基字母表,在跨越富含重复序列或可变剪接结合区时易引发错配与假阳性注释。
酶促捕获与交联免疫共沉淀(Aza-IP 与 miCLIP)
为了绕开重亚硫酸盐诱导的剧烈化学断裂,学界开发了基于修饰酶共价捕获特性的单碱基分析策略。在 5-氮杂胞苷介导的免疫共沉淀(Aza-IP)中,通过向细胞培养基中脉冲掺入核苷类似物 5-氮杂胞苷(5-azacytidine),当内源甲基转移酶(如 NSUN2、DNMT2)利用其活化位点的亲核半胱氨酸残基进攻胞嘧啶类似物的 C6 位进行催化时,活化中间体无法完成 β-消除过程,从而在酶蛋白与受体 RNA 之间形成不可逆的共价交联锁定。经由免疫沉淀拉下目标靶标后,交联位点处发生的逆转录通读往往伴随高频的 C-to-G 碱基颠换现象,提供单碱基分辨率的酶作用靶点图谱。类似地,甲基转移酶交联突变免疫沉淀(miCLIP)利用将甲基转移酶催化口袋的关键半胱氨酸突变为丙氨酸(如 NSUN2 C321A)的工程化突变体,在缺乏外源类似物的情况下使底物 RNA 滞留并交联沉淀,经逆转录阻断精准锁定靶标坐标。
免亚硫酸盐与脱氨酶融合追踪(DRAM 与无损测序新途径)
为了在微量生物样本中实现低损伤单碱基 m5C 解析,新型酶学编辑测序技术相继涌现,其代表为脱氨酶与阅读器辅助甲基化分析技术(DRAM)。DRAM 系统避免了苛刻的外源化学试剂,将脱氨酶活性结构域(如胞苷脱氨酶 APOBEC1 或脱氧腺苷脱氨酶突变体 TadA-8e)与 m5C 专属阅读蛋白(如 ALYREF 或 YBX1)进行融合表达。当修饰阅读模块精准附着于内源转录本的 m5C 位点附近时,人工偶联的脱氨酶催化模块便在空间邻近效应驱动下,对修饰位点两侧单链区域特定的未修饰胞嘧啶或腺嘌呤施加高效单碱基定点脱氨修饰,产生 C-to-U 或 A-to-I 的点突变特征。在对照体系中引入丧失 m5C 结合活性的突变阅读器(如引入 ALYREF K171A 或缺失 YBX1 结合核心基序的 DRAMmut),能够精准扣除胞苷脱氨酶固有的非特异性背景编辑。这种无损原位编辑策略为解析神经元及胚胎发育组织中极微量 RNA 的 m5C 分布开辟了可行路径。
m7G 检测方案:硼氢化钠特异还原、苯胺链切断裂与帽子结构分离
生物学调控轴线与空间分布特征
作为RNA修饰领域中横跨两类空间拓扑的代表,m7G 分布于真核转录本两套截然不同的生化微环境中。在 5′ 帽子结构定位上,转录起始阶段由三磷酸鸟苷转移酶和鸟嘌呤-N7-甲基转移酶(RNMT-RAM 复合体)协同装配,形成经典的 m7GpppN 结构,赋予转录本抵抗外切核酸酶降解的能力,并协助招募真核翻译起始复合物。在内部位点特异性修饰上,甲基转移酶样蛋白 1 与 WD 重复结构域蛋白 4 组成的 METTL1-WDR4 催化复合物,专一性催化特定 tRNA(如第46位)及 microRNA 前体形成内部 m7G 修饰。该修饰对维持 tRNA 核心三级折叠构象及翻译保真度展现出显著影响,其表达失调与小头畸形原始侏儒症以及恶性肿瘤转化密切相关。此外,小核糖体亚基 rRNA 内部(如酵母 1575 位或人 1639 位)亦由 Bud23-TRMT112 催化甲基化以维系核糖体生物发生的正常进程。
帽子结构与内部位点的酶学分离前置策略
由于 mRNA 5′ 端帽子的 m7G 丰度在摩尔比上远超转录本内部位点数个数量级,若缺乏彻底的分离纯化,5′ 端的强信号将掩盖内部转录后动态调控信号。前置纯化主要依赖特异性酶学拆分路径,利用双功能焦磷酸酶(如非典型脱帽酶 DXO 或烟草酸性焦磷酸酶 TAP)定向水解 5′ 端的三磷酸桥连键,生成单磷酸 5′ 末端并与下游富集步骤隔离,随后方可针对留存的链内结构展开特异化学标记与测序建库。
特异性亲核还原与苯胺诱导的链切断裂(TRAC-seq 与 AlkAniline-seq)
由于鸟嘌呤第7位氮原子带有局域正电荷,大幅增强了嘌呤环咪唑部分的亲电反应倾向。这一独特性质奠定了基于化学选择性裂解测序方案的基石。在 TRAC-seq(tRNA 还原与裂解测序)及 AlkAniline-seq 实验流程中,化学断裂严密遵循两个连续阶段。首先进入选择性硼氢化钠还原阶段。在冰浴及弱碱性缓冲介质中,使用新配制的硼氢化钠(NaBH4)处理纯化后的 RNA 链,NaBH4 专一性进攻带正电荷的 m7G 咪唑环 C8 位,将 C8=N7 双键还原为 7,8-二氢鸟苷衍生物;该还原产物自发水解脱去碱基,形成带有活性的无碱基脱嘌呤糖环(abasic site)。随后进入苯胺催化裂解阶段,向体系中加入酸性苯胺缓冲液(冰醋酸与水及苯胺的配比将体系维持在 pH ≈ 4.5),苯胺分子亲核进攻该无碱基糖环,催化连续的 β-消除与 δ-消除反应,断裂该位点两侧的核糖-磷酸二酯键。断裂反应在切口下游暴露出均一的 5′-单磷酸末端,使得后续流程中能够利用 T4 RNA 连接酶实现特异性 5′ 接头的定向连接,在随后的二代测序截断分析中产生单碱基分辨率的起始读段高峰。
预处理消除多重化学干扰
在针对小 RNA(如 tRNA)的检测中,除 m7G 外,序列内部通常富含 3-甲基胞苷(m3C)、双氢尿嘧啶(D)及 5-羟基胞苷(ho5C)等不稳定修饰,这些结构在碱性或苯胺处理时亦呈现出弱自发裂解倾向。为降低非特异性干扰,完善的 TRAC-seq 体系在进行 NaBH4 还原前,引入重组 AlkB 蛋白预先进行酶促去甲基化处理。由于天然野生型 AlkB 能够高效识别并擦除单链中的 m1A 与 m3C,却对嵌入嘌呤环内的 m7G 缺乏催化活性,这种选择性预消除步骤显著净化了背景杂质,使最终检测到的苯胺裂解峰聚焦于 METTL1 催化的 m7G 靶标。
逆转录突变分析图谱(m7G-MaP-seq 与 BoRed-seq)
除直接裂解 RNA 骨架之外,化学修饰导致的逆转录酶通读错配同样被拓展应用于内部 m7G 的精确捕获。在 m7G-MaP-seq 工作流中,经 NaBH4 温和还原后形成的无碱基糖环结构,在特定的二价阳离子(如 Mn2+ 替代 Mg2+)与工程化逆转录酶环境下,促使逆转录酶在跨越损伤位点时高频掺入非同源 dNTP,产生高达 15%–25% 的高对比度单核苷酸错配率,直接由序列突变图谱反映其精确化学修饰坐标。
核心检测方案横向多维对比
为在实验体系搭建时评估适宜的技术路线,针对 m1A、m5C 与 m7G 的主流检测方案的技术边界、起始量门槛与假阳性风险梳理如下:
| 修饰 | 主流方案 | 核心反应原理 | 单碱基分辨 | 推荐起始量 | 主要背景/假阳性来源 | 质控对冲策略 |
|---|---|---|---|---|---|---|
| m1A | m1A-MAP / TGIRT-seq | N1位正电荷空间位阻诱导逆转录突变与截断 | 具备 | 1–5 μg Poly(A) RNA | 5′端非模板掺入、tRNA片段污染、SNP混淆 | 平行引入 AlkB / ALKBH3 去甲基化酶对照 |
| m1A | MeRIP-seq(富集法) | m1A 表位单克隆抗体物理富集 | 缺乏(约百nt富集峰) | >10 μg Total RNA | 抗体对 5′Cap+1/Cap+2 交叉反应、非特异吸附 | 弱酸性反应规避 Dimroth 重排,结合去甲基化酶验证 |
| m5C | RNA-BS-seq | 亚硫酸氢盐使未修饰 C 脱氨为 U,m5C 抵抗 | 具备 | 2–10 μg Poly(A) RNA | 高温酸碱致 RNA 降解、结构区脱氨不完全 | 体外转录未修饰 Spike-in 监测未转化率 |
| m5C | Aza-IP / miCLIP | 5-氮杂胞苷共价结合催化酶诱发 C-to-G 颠换 | 具备 | 细胞量 >1×107 | 代谢类似物干扰、低表达转录本交联效率低 | 催化突变体(如 NSUN2 C321A)平行交联对照 |
| m5C | DRAM-seq | 识别蛋白融合脱氨酶诱导邻近碱基定点突变 | 具备(邻位标记) | 微量细胞 / 100 ng Total RNA | 游离脱氨酶脱靶、识别蛋白结合非甲基化序列 | 结合缺陷型突变体(如 ALYREF K171A)扣除噪音 |
| m7G | TRAC-seq | NaBH4 还原+酸性苯胺诱导骨架特异断裂 | 具备 | 10–50 μg Total RNA(需富集小RNA) | 其他修饰化学自发裂解、末端水解假象 | 还原前用野生型 AlkB 去除 m1A 与 m3C |
| m7G | m7G-MaP-seq | NaBH4 还原形成无碱基位点诱导逆转录错配 | 具备 | 2–5 μg Total RNA | 5′帽子高丰度信号压制、酶延伸过度终止 | 双功能脱帽酶(如 DXO)前置剥离端点修饰 |
技术跃迁前瞻:单分子直接测序与多重修饰动态互作
随着第三代单分子测序技术的发展,牛津纳米孔技术(ONT)直接 RNA 测序(direct RNA sequencing, DRS)为表观转录组研究提供了全新的观察维度。原生单链 RNA 分子在跨膜电压驱动下穿过蛋白质纳米孔,由于不同核苷酸碱基及甲基取代团在空间体积与表面电荷上的差异,引起流经微孔的离子电流产生皮安(pA)级别的特征性电平阻滞。这种物理传感机制展现出超越传统短读长测序的显著特征:摆脱了体外逆转录为 cDNA 及后续 PCR 扩增的生化流程,避免了逆转录酶延伸停滞与序列扩增偏倚带来的系统误差;完整保留了单条天然转录本上的长程序列信息,支持在数千碱基长度的物理分子上同时记录多个修饰位点的共存状态。纳米孔直接测序在精准分辨 m1A、m5C 与 m7G 时面临的核心难点集中在信号解码层面。由于纳米微孔感受区同时容纳约 5 个连续核苷酸(5-mer),相邻碱基的排列对基线电流构成明显的背景叠合效应,化学性质接近的修饰基团对电流幅值与驻留时间(dwell time)造成的扰动特征可能发生重叠。针对这一瓶颈,建立体外转录未修饰对照(IVT)基准模型、整合修饰酶敲除株多维数据,结合基于深度学习的残差卷积网络解码算法,正在逐步降低对复杂微观修饰位点的判读模糊度。与此同时,细胞生理调控中不同修饰之间表现出紧密的协同与动态制衡。内部 m7G 与非编码 tRNA 的精细修饰直接调节翻译机器读取三联体密码子的延伸速率,而成熟信使 RNA 上由 NSUN2 催化的 m5C 与由 METTL3 催化的 m6A,在介导转录本核输出与胞质液-液相分离凝聚体分配中展现出交替配合的调控倾向。从传统基于抗体富集的粗放识别,演进至基于特征化学反应与单碱基突变图谱的精准解析,再推进到单分子直接传感技术的开发应用,表观转录组学研究正迈向高分辨率与绝对计量的全新阶段,持续为揭示复杂生命活动中RNA分子的精细调控法则提供生化基石。
相关检测产品
本文涉及的 m5C、m1A、m6A、m7G 四种 RNA 修饰,常誉生物提供 RayBiotech 全线检测产品,按研究目的分为三条路线(完整对比见 RNA 修饰研究专题):
整体水平定量 — ELISA(EIA 系列)
- m1A ELISA 试剂盒(EIA-m1A) — 标准曲线定量
- m5C ELISA 试剂盒(EIA-m5C) — 标准曲线定量
- m6A ELISA 试剂盒(EIA-m6A) — 标准曲线定量
- m7G ELISA 试剂盒(EIA-m7G) — 标准曲线定量
组合半定量 — Dot Blot(DB 系列)
- 四合一 RNA 甲基化 Dot Blot 组合试剂盒(DB-4-RNA) — m1A/m5C/m6A/m7G 一次检测
- 六合一 RNA 修饰 Dot Blot 组合试剂盒(DB-6-RNA) — 增加 ac4C 与 A-I 编辑
- m1A Dot Blot 试剂盒(DB-m1A) · m5C Dot Blot 试剂盒(DB-m5C) · m6A Dot Blot 试剂盒(DB-m6A) · m7G Dot Blot 试剂盒(DB-m7G)
富集与定量 — MeRIP 系列(QMERIP / QMEDIP)
- m1A RNA 富集与定量试剂盒(QMERIP-M1A) — MeRIP + qPCR
- m6A RNA 富集与定量试剂盒(QMERIP-M6A) — MeRIP + qPCR
- m7G RNA 富集与定量试剂盒(QMERIP-M7G) — MeRIP + qPCR
- m5C 富集与定量试剂盒(QMEDIP-M5C) — MeRIP + qPCR
以上产品均仅供科研使用(Research Use Only),具体规格、报价与货期请联系常誉生物。
参考文献
- Dominissini, D., et al. (2016). The dynamic N1-methyladenosine methylome in eukaryotic messenger RNA. Nature, 530(7591), 441–446. DOI: 10.1038/nature16998
- Safra, M., et al. (2017). The m1A landscape on cytosolic and mitochondrial mRNA at single-base resolution. Nature, 551(7679), 251–255. DOI: 10.1038/nature24456
- Li, X., et al. (2017). Base-resolution mapping reveals distinct m1A methylome in nuclear- and mitochondrial-encoded transcripts. Molecular Cell, 68(5), 993–1005. DOI: 10.1016/j.molcel.2017.10.019
- Squires, J. E., et al. (2012). Widespread occurrence of 5-methylcytosine in human coding and non-coding RNA. Nucleic Acids Research, 40(11), 5023–5033. DOI: 10.1093/nar/gks144
- Khoddami, V., & Cairns, B. R. (2013). Identification of direct targets and modified bases of RNA cytosine methyltransferases. Nature Biotechnology, 31(5), 458–464. DOI: 10.1038/nbt.2566
- Yang, X., et al. (2017). 5-methylcytosine promotes mRNA export—NSUN2 as the methyltransferase and ALYREF as an m5C reader. Cell Research, 27(5), 606–625. DOI: 10.1038/cr.2017.55
- Chen, X., et al. (2019). 5-methylcytosine promotes pathogenesis of bladder cancer through stabilizing mRNAs. Nature Cell Biology, 21(8), 978–990. DOI: 10.1038/s41556-019-0361-y
- Lin, S., et al. (2018). Mettl1/Wdr4-mediated m7G tRNA methylome is required for normal mRNA translation and embryonic stem cell self-renewal and differentiation. Molecular Cell, 71(2), 244–255. DOI: 10.1016/j.molcel.2018.06.001
- Marchand, V., et al. (2018). AlkAniline-Seq: Profiling of m7G and m3C RNA Modifications at Single Nucleotide Resolution. Angewandte Chemie International Edition, 57(51), 16785–16790. DOI: 10.1002/anie.201810946
- Zhang, L. S., et al. (2019). Transcriptome-wide mapping of internal N7-methylguanosine methylation in mammalian mRNA. Molecular Cell, 74(6), 1304–1316. DOI: 10.1016/j.molcel.2019.03.036