解码mRNA：高置信度鉴定m5C修饰及序列结构特征探究

解码mRNA：高置信度鉴定m5C修饰及序列结构特征探究一、引言1.1研究背景与意义在生命科学领域，核酸修饰的研究一直是热点话题，其中5-甲基胞嘧啶（m5C）修饰因其广泛存在于多种RNA种类中而备受关注。m5C修饰最初在转运RNA（tRNA）和核糖体RNA（rRNA）中被大量发现，近年来随着检测技术的不断进步，其在信使RNA（mRNA）中的存在也逐渐被证实。早期对于mRNA上m5C修饰的研究相对滞后，但随着亚硫酸盐测序（BS-seq）、5-氮杂胞苷交联测序（Aza-IP）、单核苷酸分辨率交联免疫沉淀测序（miCLIP-seq）等高通量技术的迅速发展，尤其是超快速亚硫酸氢盐测序（UBS-seq）和m5C-TAC-seq等新技术的出现，人们对mRNA中m5C修饰的认识有了显著进展。然而，由于mRNAm5C修饰相对较少，且存在低丰度和高背景噪声的问题，使得对转录组范围内的m5C图谱进行绘制面临诸多挑战，不同研究中估计的mRNAm5C位点数量差异较大，观察结果也不一致。例如，有的研究利用BS-seq在人类HeLa细胞和多种小鼠组织中鉴定出数千个m5C位点，而有的研究结果却与之相差甚远，这使得难以定义一组通用的mRNA底物或常见的甲基化目标序列。高置信度地鉴定mRNA上的m5C修饰位点至关重要。精确的位点鉴定是深入研究m5C修饰功能的基础。只有准确确定m5C修饰的位置，才能进一步探究其对mRNA的稳定性、翻译效率、剪接过程等方面的影响。以mRNA的稳定性为例，m5C修饰可能通过改变mRNA的二级结构，影响其与相关蛋白的相互作用，从而调控mRNA的半衰期。如果无法准确鉴定修饰位点，就无法深入理解这种调控机制。准确鉴定m5C修饰位点对于揭示其在生物过程中的作用不可或缺。在胚胎发育过程中，m5C修饰在母源mRNA上的动态变化被发现与母源-合子过渡密切相关。通过高置信度鉴定m5C修饰位点，研究人员发现母源mRNA在卵母细胞和胚胎发育早期存在大规模的m5C甲基化，随着母源-合子过渡的发生，m5C位点数目和甲基化水平急剧下降。这一发现暗示了m5C修饰在胚胎发育早期对母源mRNA的调控具有重要意义。mRNA上m5C修饰的异常与多种疾病的发生发展紧密相连，对疾病机制研究和治疗策略开发意义重大。在肿瘤领域，已有研究表明m5C修饰水平的改变与肿瘤的增殖、转移和耐药性相关。通过高置信度鉴定肿瘤细胞中mRNA的m5C修饰位点，有助于发现潜在的肿瘤生物标志物和治疗靶点。在一些癌症中，特定mRNA上的m5C修饰异常可能导致相关基因的表达失调，进而促进肿瘤的发生发展。深入研究这些异常修饰位点，有望为癌症的早期诊断和精准治疗提供新的思路和方法。在神经系统疾病等其他领域，m5C修饰也可能在疾病的病理过程中发挥作用，高置信度鉴定修饰位点将为揭示疾病机制和开发治疗手段提供关键线索。1.2国内外研究现状在mRNA的m5C修饰鉴定方法上，国内外科研人员都做出了卓越的努力。早期，液相色谱-质谱联用技术（LC-MS/MS）能够对RNA甲基化整体水平进行鉴定，实现对碱基的定性和定量分析，但无法精准定位修饰位点。随后，亚硫酸盐测序（BS-seq）技术兴起，其原理是亚硫酸盐处理使未甲基化的胞嘧啶（C）脱氨基转变为尿嘧啶（U），而m5C保持不变，结合高通量测序技术，实现了单碱基分辨率的全基因组DNA甲基化图谱绘制，为mRNAm5C修饰位点的鉴定带来了曙光。国外如美国、德国等科研团队率先运用BS-seq技术在人类细胞和小鼠组织中开展研究，国内中山大学张锐教授团队也利用该技术建立了一套基于富集的mRNA的RNA亚硫酸氢盐测序及生信分析流程，用以高置信度鉴定mRNAm5C甲基化，为后续研究奠定了重要基础。随着技术的发展，超快速亚硫酸氢盐测序（UBS-seq）利用高浓度亚硫酸盐和高反应温度，缩短反应时间减少RNA降解，同时提高C到U的转化效率，进一步优化了检测效果。2024年，北京大学伊成器、浙江大学李笑雨共同通讯发表研究，报道了m5C-TAC-seq技术，该技术结合TET辅助的m5C-f5C氧化和选择性化学标记，实现了碱基分辨率的m5C检测，且反应条件温和，可在低丰度和低序列复杂性的RNA中检测m5C，为mRNAm5C修饰鉴定提供了新的有力工具。在mRNAm5C修饰的序列结构特征研究方面，国外研究先行一步，通过大量的实验和数据分析，初步揭示了一些规律。例如，部分研究表明m5C修饰在mRNA的非翻译区（UTR）和近Argonaute蛋白结合区域存在富集现象。国内张锐教授团队则通过深入研究，发现mRNA存在两类m5C位点（TypeI和TypeII位点），具有独特的序列和结构特征。TypeI位点在多数组织和细胞类型中富集，位于茎-环结构的5’端，3’端有富含鸟嘌呤（G-rich）的基序，由NSUN2介导甲基化，其序列和结构特征与tRNA底物高度一致；TypeII位点具有较强的组织特异性，占比在不同组织中有所差异，位于茎-环结构的环区，3’端带有TCCA基序，其甲基化由NSUN6介导。这些发现进一步深化了对mRNAm5C修饰序列结构特征的认识，为后续功能研究提供了关键线索。关于mRNAm5C修饰的功能研究，国内外都取得了一定的进展。国外研究发现，m5C修饰参与了mRNA的稳定、剪接和核质穿梭等过程，对基因表达调控起着重要作用。在胚胎发育研究中，国外团队发现m5C修饰在早期胚胎发育中存在动态变化，对胚胎发育过程产生影响。国内研究也有诸多重要发现，如张锐教授团队通过构建多个物种不同发育时期的mRNAm5C图谱，意外揭示了母源mRNA的大规模m5C甲基化现象。研究发现，在卵母细胞和胚胎发育早期，m5C高度富集，随着母源-合子过渡的发生，m5C位点数目和甲基化水平急剧下降，且这种低水平甲基化维持到个体发育成熟。通过敲除果蝇中的NSUN2蛋白，导致96%的母源mRNAm5C被清除，进而造成果蝇胚胎发育迟缓，暗示了m5C在早期胚胎发育中的重要功能，这与之前国内其他团队干扰mRNAm5C读取蛋白YBX1在斑马鱼中观察到的表型相互验证，充分表明了m5C修饰在胚胎发育过程中的关键调控作用。1.3研究内容与方法本研究旨在高置信度鉴定mRNA上的m5C修饰，并深入探究其序列结构特征，具体研究内容与方法如下：mRNA的提取与纯化：选取多种细胞系和组织样本，包括人类HeLa细胞、小鼠胚胎干细胞以及小鼠的心脏、肝脏、肾脏等组织。运用Trizol试剂法提取总RNA，利用寡聚（dT）磁珠富集mRNA，以确保获取高纯度的mRNA样本，为后续实验奠定基础。m5C修饰位点的鉴定：采用超快速亚硫酸氢盐测序（UBS-seq）技术，该技术使用高浓度亚硫酸盐（约10M）和98°C的高反应温度，通过缩短反应时间减少RNA降解，并提高反应温度实现完全的C到U转化，从而更高效准确地鉴定m5C修饰位点。同时，运用m5C-TAC-seq技术进行验证，该技术结合TET辅助的m5C-f5C氧化和选择性化学标记，实现碱基分辨率的m5C检测，且反应条件温和，可在低丰度和低序列复杂性的RNA中检测m5C。对测序数据进行生物信息学分析，通过严格的筛选标准，如设定甲基化水平阈值、排除低覆盖率位点等，高置信度地识别m5C修饰位点。mRNAm5C修饰序列结构特征分析：利用生物信息学工具，对鉴定出的m5C修饰位点周围的序列进行分析，寻找保守的基序（motif）。例如，使用MEME软件进行基序搜索，探究m5C修饰位点上下游碱基组成的偏好性。借助RNA二级结构预测软件，如RNAfold，预测mRNA的二级结构，分析m5C修饰位点在二级结构中的位置分布特征，研究其是否倾向于位于茎-环结构、发夹结构等特定结构区域。m5C修饰相关蛋白的研究：通过蛋白质免疫共沉淀（co-IP）技术，使用针对m5C修饰相关蛋白（如NSUN2、NSUN6等甲基转移酶）的抗体，从细胞裂解液中富集与mRNA结合的蛋白复合物，然后进行质谱分析，鉴定与m5C修饰mRNA相互作用的蛋白。运用RNA免疫沉淀测序（RIP-seq）技术，进一步确定这些蛋白在mRNA上的结合位点，分析其与m5C修饰位点的相关性，探究它们在m5C修饰过程中的作用机制。二、mRNA上m5C修饰概述2.1m5C修饰的定义与形成机制5-甲基胞嘧啶（m5C）修饰是一种重要的RNA修饰类型，指的是在RNA中胞嘧啶碱基的第5位碳原子上添加一个甲基基团（-CH3），从而形成5-甲基胞嘧啶（m5C）。这一修饰过程是由特定的甲基转移酶催化完成的，在生物体内，S-腺苷甲硫氨酸（SAM）通常作为甲基供体，为m5C修饰提供活性甲基。甲基转移酶利用SAM提供的甲基，将其转移到RNA中特定的胞嘧啶残基上，使胞嘧啶转变为m5C。在mRNA中，m5C修饰的形成过程同样依赖于这些甲基转移酶。目前已知的RNAm5C甲基转移酶家族包括DNMT2、TRDMT家族成员、NSUN家族等多个成员。其中，NSUN家族在mRNAm5C修饰中发挥着关键作用。NSUN家族包含多个成员，如NSUN1到NSUN7以及NSUN5a/b/c等，不同成员在mRNAm5C修饰中具有不同的底物特异性和功能。NSUN2是研究较为深入的一种m5C甲基转移酶，它在多种细胞和组织中参与mRNA的m5C修饰过程。在小鼠胚胎干细胞中，NSUN2被发现参与调控一系列基因的mRNAm5C修饰，对细胞的增殖、分化等过程产生影响。m5C修饰并非一成不变，在生物体内存在着动态的调控机制。去甲基化酶可以介导m5C修饰的可逆性。TET酶家族是一类重要的去甲基化酶，在mRNAm5C修饰的调控中发挥作用。TET1-3具有作为RNA去甲基化酶的潜力，其过表达可以显著提高HEK293T细胞中RNA5hmC（5-羟甲基胞嘧啶，m5C去甲基化的中间产物）的水平，这表明TET酶家族能够参与mRNAm5C修饰的动态调控，使m5C修饰在不同的生理和病理条件下发生变化，从而影响mRNA的功能和相关生物过程。2.2m5C修饰在不同RNA中的分布特点m5C修饰在tRNA、rRNA和mRNA等不同RNA中呈现出明显的分布差异，这些差异与各类RNA的结构和功能密切相关。在tRNA中，m5C修饰较为普遍且具有特定的分布模式。研究表明，tRNA甲基化最常发生在可变环和T-stem连接区域的胞嘧啶及47-50位点的1-3个C上。这些位置的m5C修饰对tRNA的结构和功能至关重要。在大肠杆菌的tRNA中，m5C修饰位于T-ψ-C环附近，该修饰有助于稳定tRNA的L形三维结构，使得tRNA能够更准确地识别密码子，从而保证蛋白质合成的准确性。m5C修饰还可以优化密码子-反密码子配对，调节应激反应，对翻译效率和准确性产生重要影响。当细胞处于应激状态时，tRNA上特定位置的m5C修饰水平会发生变化，进而调控蛋白质合成的速率，以适应细胞的生存需求。rRNA中的m5C修饰也具有独特的分布特点。rRNA是核糖体的重要组成部分，其m5C修饰与核糖体的结构和功能密切相关。在真核生物的rRNA中，m5C修饰主要分布在一些保守区域，这些区域参与核糖体与mRNA、tRNA的相互作用。在人类的18SrRNA中，特定位置的m5C修饰被发现与核糖体的组装和功能密切相关，影响蛋白质的合成效率。m5C修饰还可能影响rRNA-tRNA-mRNA三级复合体的结构稳定性，从而在蛋白质合成过程中发挥作用。在神经胶质瘤细胞中，研究发现rRNAm5C修饰与对应激相关酶NQO1的生物活性底物的敏感性有关，表明m5C修饰在rRNA中的分布与细胞的应激反应和代谢调控存在关联。mRNA中的m5C修饰分布与tRNA和rRNA有所不同，且研究发现其具有组织特异性和区域偏好性。早期研究利用亚硫酸盐测序在人类HeLa细胞和多种小鼠组织中绘制转录组范围内的m5C图谱，发现m5C修饰主要位于编码序列（CDS）中，主要在CG环境中，以及mRNA翻译起始位点的下游区域。而另一项研究则观察到m5C位点在翻译起始密码子附近（如5'UTR末端和CDS起始处）有显著积累，在编码序列中耗竭，并在3'UTR中有混合的富集模式。中山大学张锐教授团队的研究进一步揭示了mRNA存在两类m5C位点（TypeI和TypeII位点），具有独特的序列和结构特征。TypeI位点在多数组织和细胞类型中富集，位于茎-环结构的5’端，3’端有富含鸟嘌呤（G-rich）的基序；TypeII位点具有较强的组织特异性，占比在不同组织中有所差异，位于茎-环结构的环区，3’端带有TCCA基序。这些研究表明，mRNA中m5C修饰的分布具有复杂性和多样性，可能在mRNA的不同加工过程和功能调控中发挥作用。2.3m5C修饰对mRNA功能的潜在影响m5C修饰在mRNA的生命周期中发挥着多方面的潜在影响，对mRNA的稳定性、剪接、核质穿梭和翻译等关键功能都有着重要的调控作用。在mRNA稳定性方面，m5C修饰扮演着关键的角色。研究表明，m5C修饰可以通过多种机制影响mRNA的稳定性。在斑马鱼早期胚胎发育过程中，m5C修饰的母源mRNA比未修饰的mRNA更稳定。进一步研究发现，m5C结合蛋白Ybx1通过其CSD结构域关键残基Trp45与m5C的π-π相互作用，识别m5C修饰的mRNA，进而招募Pabpc1a维持母源mRNA的稳定性，保证母源-合子转换有序进行。在肿瘤细胞中，NSUN2作为主要的m5C甲基转移酶，通过m5C修饰调控靶基因HGH1的表达，增强HGH1mRNA的稳定性，从而促进乳腺癌的发展。这些研究充分说明，m5C修饰可以通过与特定的结合蛋白相互作用，保护mRNA免受核酸酶的降解，延长其半衰期，进而影响mRNA的稳定性。mRNA的剪接过程也受到m5C修饰的影响。虽然目前对于m5C修饰在mRNA剪接中具体作用机制的研究还相对较少，但已有研究表明，m5C修饰可能通过改变mRNA的二级结构，影响剪接因子与mRNA的结合，从而调控mRNA的剪接过程。mRNA的二级结构对剪接过程至关重要，m5C修饰可能会改变mRNA的局部构象，使得剪接位点更容易或更难被剪接因子识别，进而影响剪接体的组装和剪接过程的进行。此外，m5C修饰还可能与其他RNA修饰协同作用，共同调控mRNA的剪接。在某些情况下，m5C修饰与m6A修饰可能在同一mRNA分子上同时存在，它们之间可能通过相互影响mRNA的结构和与蛋白质的相互作用，对mRNA的剪接产生综合影响。m5C修饰在mRNA的核质穿梭过程中发挥着重要作用。mRNA需要从细胞核转运到细胞质中才能进行翻译，而m5C修饰可以通过与特定的转运蛋白相互作用，促进mRNA的核输出。ALYREF是细胞核中识别的第一个mRNA阅读蛋白，具有K171的关键m5C识别位点。m5C修饰的mRNA在细胞核中被ALYREF识别后，有助于其输出到细胞质。这一过程对于保证mRNA能够及时到达细胞质进行翻译至关重要。如果m5C修饰异常，可能会导致mRNA核质穿梭受阻，进而影响蛋白质的合成。在一些疾病状态下，如肿瘤细胞中，m5C修饰的改变可能会影响ALYREF与mRNA的结合，导致mRNA在细胞核内积累，无法正常转运到细胞质中进行翻译，从而影响细胞的正常生理功能。m5C修饰对mRNA的翻译过程也具有重要影响。研究发现，m5C修饰可以影响mRNA与核糖体的结合效率，从而调控翻译起始的速率。在蛋白质合成过程中，核糖体需要准确地识别mRNA的起始密码子并与之结合，才能启动翻译过程。m5C修饰可能会改变mRNA起始密码子附近的结构，影响核糖体与mRNA的结合亲和力，进而影响翻译起始的效率。此外，m5C修饰还可能通过影响翻译延伸过程，对蛋白质的合成产生影响。在翻译延伸过程中，tRNA携带氨基酸按照mRNA的密码子顺序依次添加到多肽链上，m5C修饰可能会影响tRNA与mRNA的配对效率，从而影响翻译延伸的速率和准确性。在一些研究中，发现m5C修饰在某些基因的mRNA上的存在与蛋白质合成的速率和质量密切相关，进一步说明了m5C修饰在mRNA翻译过程中的重要作用。三、高置信度鉴定mRNA上m5C修饰的方法3.1传统鉴定方法及局限性在mRNA上m5C修饰鉴定的发展历程中，传统方法如亚硫酸盐测序（BS-seq）、5-氮杂胞苷交联测序（Aza-IP）、单核苷酸分辨率交联免疫沉淀测序（miCLIP-seq）等曾发挥了重要作用，但它们在检测mRNAm5C修饰时也暴露出诸多局限性。亚硫酸盐测序（BS-seq）是较早用于mRNAm5C修饰鉴定的经典方法之一，其原理基于亚硫酸盐处理能够使未甲基化的胞嘧啶（C）脱氨基转变为尿嘧啶（U），而m5C保持不变这一特性。通过结合高通量测序技术，理论上可以实现单碱基分辨率的全基因组DNA甲基化图谱绘制，也为mRNAm5C修饰位点的鉴定提供了可能。然而，在实际应用于mRNAm5C修饰检测时，BS-seq存在显著的局限性。一方面，该方法依赖于未修饰C的高效转换，若转换不完全，极易导致假阳性结果的出现。在实际反应过程中，由于mRNA的结构复杂性，亚硫酸盐可能无法完全接触到所有未甲基化的胞嘧啶，从而使得部分未甲基化的C未被转化为U，在后续测序分析中被误判为m5C修饰位点。另一方面，BS-seq的反应条件较为苛刻，需要较高的温度和较长的反应时间，这会不可避免地导致RNA降解。mRNA本身相对不稳定，在这种苛刻的反应条件下，更容易发生降解，这不仅限制了低起始量样本和低丰度RNA的检测，还可能导致检测到的m5C修饰位点信息不准确，因为降解的mRNA片段可能无法完整地反映真实的修饰情况。由于C转变为U后序列复杂度降低，会影响比对准确度，尤其在检测低序列复杂度RNA上的m5C时，这一问题更为突出，使得难以准确鉴定这些区域的m5C修饰位点。5-氮杂胞苷交联测序（Aza-IP）也是一种传统的m5C修饰鉴定方法。该方法利用5-氮杂胞苷（5-aza-C）这一胞嘧啶类似物，在转录过程中，RNA聚合酶会随机将其掺入到初生（nascent）RNA转录本的胞嘧啶位点。在m5C-RCMT催化结构域中，半胱氨酸残基的硫原子与目标RNA碱基的C6位置形成共价键，随后通过甲基供体S-腺苷甲硫氨酸（SAM）的作用，目标胞嘧啶的C5位置发生甲基化。在C5甲基化后，通过后续的β-消除反应断裂共价键，从而恢复自由酶。然而，由于5-氮杂胞苷中C5位置的碳原子被氮原子（N）取代，共价键得以稳定，最终导致细胞内内源性酶的耗竭，从而使RNA和DNA的甲基化水平降低。这种方法存在的问题是，5-氮杂胞苷的掺入具有随机性，可能会影响RNA的正常结构和功能，进而干扰m5C修饰的准确检测。该方法对于低丰度的mRNAm5C修饰检测灵敏度较低，容易遗漏一些真实的修饰位点，且实验操作相对复杂，对实验条件的控制要求较高。单核苷酸分辨率交联免疫沉淀测序（miCLIP-seq）则是基于甲基转移酶与内源性RNA靶标发生交联的原理来鉴定m5C修饰位点。在甲基化过程中，保守的半胱氨酸残基Cys1暂时与甲基化胞嘧啶形成共价键，而第二个半胱氨酸Cys2对于这种催化中间体的分解至关重要。在miCLIP中，Cys2突变为丙氨酸（Ala），便于捕获催化中间体，并允许甲基转移酶与内源性RNA靶标发生交联，而无需光交联，最终形成交联的表位标记的酶-底物复合体。尽管miCLIP-seq能够在一定程度上实现单核苷酸分辨率的m5C修饰位点鉴定，但它也面临着挑战。该方法需要特异性的抗体来识别交联的复合体，抗体的质量和特异性会直接影响实验结果的准确性。如果抗体的特异性不佳，可能会与非目标蛋白或RNA发生非特异性结合，导致假阳性信号的产生，从而干扰对真实m5C修饰位点的判断。miCLIP-seq的实验流程较为繁琐，需要进行多次免疫沉淀和纯化步骤，这不仅增加了实验操作的难度和时间成本，还可能在操作过程中引入误差，影响实验的重复性和可靠性。3.2新型高置信度鉴定技术原理与优势随着对mRNAm5C修饰研究的深入，传统鉴定方法的局限性日益凸显，促使科研人员不断探索创新，超快速亚硫酸氢盐测序（UBS-seq）和m5C-TAC-seq等新型高置信度鉴定技术应运而生，为mRNAm5C修饰研究带来了新的契机。超快速亚硫酸氢盐测序（UBS-seq）是一种针对传统亚硫酸氢盐测序（BS-seq）局限性进行优化的新型技术。其核心原理基于对BS反应机理的深入理解和巧妙改进。在传统BS-seq中，存在两种竞争途径，一种是期望的未甲基化胞嘧啶（C）到尿嘧啶（U）的转化途径，另一种则是导致DNA/RNA降解的途径。UBS-seq通过使用高浓度亚硫酸盐（约10M）和98°C的高反应温度来克服这些问题。高浓度的亚硫酸盐能够加快BS的转化率，从而减少DNA的降解。较高的温度不仅可以加速BS反应，还有助于使双链DNA（dsDNA）或RNA中的二级结构变性，进而在更短的时间内完成亚硫酸氢盐的完全转化。在对5-merDNA寡核苷酸AGCGA的实验中，使用UBS-1配方在98°C下孵育，MALDI-TOFMS显示在3分钟内胞嘧啶就完全转化为尿嘧啶-BS加合物，UBS条件使得BS转化速度比传统条件加速了约13倍。对5mCDNA寡核苷酸探针处理10分钟后未观察到5mC的明显反应，这表明UBS条件在10分钟内不会产生假阴性。相较于传统的BS-seq，UBS-seq具有显著的优势。在减少假阳性方面表现出色。4mC是细菌gDNA中另一种已知的胞嘧啶甲基化，传统的BS条件可介导约50%的4mC脱氨，因此使用传统的BS-seq检测5mC时容易产生假阳性。而在UBS-1条件下，4mC被定量转化为dU，4mC和C位点均被读取为T，有效消除了4mC引起的潜在假阳性。在检测背景方面，UBS-seq也远低于传统BS-seq。使用无5mC修饰的λ-DNA和有高度修饰CpG基序的pUC19质粒评估UBS-seq的假阳性率时发现，UBS-1条件下处理10min的C平均未转化率低至约0.06%，而传统BS处理在λ-DNA上的背景值高出13倍以上；对于更结构化的PUC19，在常规BS-seq中，非CpG位点的平均未转换率增加到2.3%，比UBS-1条件高出22倍以上。这些结果充分说明，UBS-seq中使用的更高温度（98°C）可以有效地使所有区域的dsDNA变性，提供更低且均匀分布的背景。当应用于RNAm5C检测时，UBS-seq大大减少了背景及假阳性。在HeLamRNA的检测中，研究人员利用UBS-seq发现了数千个m5C位点，其中约90%对NSUN2缺失敏感，证实UBS-seq产生非常低的假阳性率。m5C-TAC-seq是另一种新型的mRNAm5C修饰鉴定技术，其原理是将酶促反应与化学标记相结合。该技术利用优化后的TET酶促反应将RNAm5C氧化至f5C（5-醛基胞嘧啶），进一步使用叠氮茚二酮（AI）对于氧化产生的f5C进行特异性标记。标记产物不仅可以产生C-to-T的转变，还可以通过点击化学反应进行富集，从而实现对于m5C修饰单碱基分辨率的直接检测。这种直接检测的特性使其能够灵敏检测修饰比例低至2.5%的位点。在rRNA的检测中，m5C-TAC-seq不仅能够精准检测已知的高m5C修饰位点，还鉴定了一个新的低修饰位点m5C3683及其修饰酶，充分验证了其灵敏度。对于具有复杂二级结构且序列复杂度较低的tRNA，m5C-TAC-seq不仅能够鉴定所有已报道位点及其修饰酶，还由于保留了四元碱基组成，使得其可以实现对不同tRNAisodecoder的精准检测。m5C-TAC-seq的优势也十分突出。其反应条件温和，不会像传统BS-seq那样因苛刻的反应条件导致RNA降解，这使得它适用于低起始量样本和低丰度RNA的检测。由于是直接检测m5C修饰，不存在传统BS-seq中因C转变为U后序列复杂度降低而影响比对准确度的问题，尤其适合低序列复杂度RNA上m5C的检测。研究人员利用m5C-TAC-seq绘制了HeLa、HEK293T和mESC的单碱基分辨率m5C修饰图谱，分别鉴定了2499、765和664个m5C位点，并利用无修饰转录组文库进行对照分析验证了这些位点的可靠性。与BS-seq相比，m5C-TAC-seq在低修饰位点的检测上也表现出高精度和高灵敏度，进一步凸显了其直接检测策略的优势。m5C-TAC-seq还结合multiplexing文库构建策略，不仅可以减少样品之间的技术误差，还能够实现m5C修饰的半定量检测。得益于温和的反应条件及对碱基组成的保留，m5C-TAC-seq实现了对含有大量低复杂度序列的chromatin-associatedRNA（caRNA）上m5C位点的鉴定，为研究m5C修饰在这类特殊RNA中的分布和功能提供了有力工具。3.3案例分析：某具体研究中高置信度鉴定方法的应用张锐教授团队在mRNAm5C修饰的研究中，巧妙运用了高置信度鉴定方法，为该领域的发展做出了卓越贡献。在实验样本的选择上，团队精心挑选了多种具有代表性的样本，包括人类的HeLa细胞、多种小鼠组织，以及果蝇、斑马鱼、非洲爪蟾、热带爪蟾、小鼠、人类这六个物种不同发育时期的样本。这些样本涵盖了不同的物种、细胞类型和发育阶段，为全面研究mRNAm5C修饰提供了丰富的数据基础。在技术手段的运用上，团队主要采用了基于富集的mRNA的RNA亚硫酸氢盐测序（BS-seq）技术，并建立了一套严谨的生信分析流程。通过这一技术，他们成功地在哺乳动物细胞系与组织中构建了mRNAm5C精确图谱。在数据处理和分析过程中，团队充分考虑到mRNAm5C修饰研究中存在的假阳性问题。他们通过深入分析发现，当时报道的“m5C”位点往往成簇存在于高GC含量区域，对针对NSUN2的RNA干扰不敏感，也无法被m5C抗体所富集，极可能是实验引入的假阳性。基于这一发现，团队建立了一套新颖的生物信息学计算流程过滤噪音，得以精确地定位mRNA上的m5C。他们通过严格设定甲基化水平阈值、排除低覆盖率位点等筛选标准，高置信度地识别m5C修饰位点。在对mRNAm5C修饰位点的研究中，团队取得了一系列重大发现。他们确定了mRNA存在两类m5C位点（TypeI和TypeII位点），并揭示了其独特的序列和结构特征。TypeI位点在多数组织和细胞类型中富集，位于茎-环结构的5’端，3’端有富含鸟嘌呤（G-rich）的基序，由NSUN2介导甲基化，其序列和结构特征与tRNA底物高度一致。TypeII位点具有较强的组织特异性，占比在不同组织中有所差异，位于茎-环结构的环区，3’端带有TCCA基序，其甲基化由NSUN6介导。这些发现不仅深化了对mRNAm5C修饰序列结构特征的认识，还为后续研究修饰位点与mRNA功能之间的关系提供了重要线索。张锐教授团队还构建了目前为止最大规模的6个物种不同发育时期的mRNAm5C图谱，意外揭示了母源mRNA的大规模m5C甲基化现象。研究发现，在整个发育过程中m5C只在卵母细胞和胚胎发育早期高度富集；随着母源-合子过渡的发生，m5C位点数目和甲基化水平均出现急剧下降，并且这种低水平甲基化一直维持到个体发育成熟。与成体组织mRNA相比，母源mRNA有多达10倍到数十倍的m5C位点，单个位点的母源mRNA的甲基化水平也远远高于成体细胞中对应位点甲基化水平。以人类卵母细胞为例，研究一共发现8592个mRNAs的32941个甲基化水平大于10%的位点。为了揭示母源mRNA大规模m5C甲基化的机制，团队研究了NSUN2的表达水平和细胞定位，发现卵母细胞的NSUN2定位于细胞质，不同于成体细胞的细胞核定位。通过nocodazole处理HeLa细胞模拟卵母细胞状态，验证了位于细胞质的NSUN2蛋白可能与母源mRNA持续相互作用并导致m5C逐渐增多的假说。在功能研究方面，团队在果蝇中敲除NSUN2蛋白，导致96%的母源mRNAm5C被清除，造成果蝇胚胎发育迟缓，暗示了m5C在早期胚胎发育中的重要功能。张锐教授团队的研究成果充分展示了高置信度鉴定方法在mRNAm5C修饰研究中的关键作用。通过准确鉴定修饰位点，深入分析其序列结构特征和动态变化规律，为进一步探究mRNAm5C修饰的生物学功能和进化模式奠定了坚实基础，也为该领域的后续研究提供了重要的参考和借鉴。四、mRNA上m5C修饰的序列结构特征分析4.1一级序列特征分析4.1.1基序分析mRNA上m5C修饰位点的一级序列特征中，基序分析是关键环节，有助于揭示m5C修饰在核苷酸序列层面的偏好性和潜在规律。通过对大量已鉴定的m5C修饰位点周围核苷酸序列的深入研究，发现了一些显著的基序模式。在众多研究中，“C/G-C/G-m5C-G/A-G-G-G”基序是较为常见的一种。中国科学院上海营养与健康研究所-马普学会计算生物学伙伴研究所杨力研究组与澳大利亚国立大学约翰克汀医学研究中心ThomasPreiss研究组合作，通过构建人HeLa细胞不同多聚核糖体组分的重亚硫酸盐RNA测序（BisulfiteRNAsequencing，bsRNA-seq），建立全转录组m5C分析的生物信息学体系（bsRNA-seq-m5C），在人类HeLa细胞系鉴定了约1000个高度可信的m5C修饰位点。计算生物学分析表明，mRNA上的m5C位点在一级序列上倾向于发生在“C/G-C/G-m5C-G/A-G-G-G”基序中。这一基序中，m5C位点两侧的核苷酸种类并非随机分布，而是具有特定的偏好性。紧邻m5C位点上游的两个位置倾向于出现胞嘧啶（C）或鸟嘌呤（G），下游第一个位置倾向于出现鸟嘌呤（G）或腺嘌呤（A），随后连续三个鸟嘌呤（G）。这种碱基组成的偏好性暗示了m5C修饰与周围核苷酸之间可能存在特定的相互作用，或者该基序为m5C修饰的发生提供了某种结构基础，使得甲基转移酶更容易识别并催化该位点的甲基化反应。不同细胞类型和物种中，m5C修饰位点的基序可能存在差异。在小鼠胚胎干细胞和某些组织细胞中，对m5C修饰位点基序的分析发现，虽然“C/G-C/G-m5C-G/A-G-G-G”基序也有出现，但频率和具体分布与HeLa细胞有所不同。在小鼠胚胎干细胞中，部分m5C修饰位点周围的基序在某些位置上的碱基偏好性更加明显，可能与胚胎干细胞独特的生物学功能和基因表达调控需求有关。在不同物种间，如人类与小鼠，由于基因组序列和基因调控机制的差异，m5C修饰位点的基序也存在一定程度的变化。这些差异可能反映了m5C修饰在不同生物进化过程中为适应各自的生物学特性而发生的演变，也提示了m5C修饰在不同物种中可能具有不同的功能和调控机制。4.1.2保守序列分析保守序列分析对于理解mRNA上m5C修饰位点的进化特征和功能具有重要意义。通过对比不同物种mRNA序列中m5C修饰位点的保守序列，能够揭示m5C修饰在生物进化历程中的演变规律，以及其与mRNA功能的紧密联系。在对多种物种的研究中发现，mRNAm5C修饰位点的保守序列呈现出复杂的特征。中山大学张锐教授团队通过对人和小鼠等多个物种的研究，构建了哺乳动物mRNAm5C精确图谱。研究发现，在不同组织中mRNAm5C位点的数量与甲基化水平有很大差异，且这些m5C位点在人和老鼠之间并不保守，说明它们很可能是受到顺式调控（cis-regulation）的。这表明在进化过程中，mRNAm5C修饰位点的保守性并非普遍存在，而是在很大程度上受到物种特异性和组织特异性的影响。进一步分析发现，尽管部分m5C修饰位点在不同物种间不保守，但在某些关键基因或功能相关的mRNA区域，仍然存在保守的m5C修饰位点。在参与细胞基本代谢途径的基因mRNA中，一些m5C修饰位点在多个物种中具有较高的保守性。这暗示这些保守的m5C修饰位点可能在维持基因的基本功能、保证细胞正常代谢过程中发挥着重要作用。在进化过程中，这些关键基因的mRNAm5C修饰位点被保留下来，以确保相关生物学过程的稳定性和准确性。从进化的角度来看，mRNAm5C修饰位点保守序列的变化可能与物种的适应性进化密切相关。随着物种的进化，基因组序列发生改变，基因调控机制也不断演变。m5C修饰作为一种重要的RNA修饰方式，其修饰位点的保守序列可能会根据物种的需求发生相应的变化。在高等动物中，伴随着富集二级结构的5'UTR区域的出现，高等动物在mRNA的5’端获得了顺式元件进化介导的NSUN2依赖的m5C位点。这表明在进化过程中，mRNAm5C修饰位点的保守序列受到顺式元件进化的影响，从而产生了新的修饰位点，以适应高等动物复杂的生物学功能和基因调控需求。mRNAm5C修饰位点的保守序列还可能与mRNA的结构和功能稳定性有关。保守的m5C修饰位点可能通过影响mRNA的二级结构或与相关蛋白的相互作用，来维持mRNA的稳定性和功能。如果这些保守位点发生改变，可能会导致mRNA结构和功能的异常，进而影响细胞的正常生理过程。在某些疾病状态下，mRNAm5C修饰位点保守序列的改变与疾病的发生发展密切相关。在肿瘤细胞中，部分关键基因mRNA上m5C修饰位点保守序列的异常变化，可能导致基因表达失调，从而促进肿瘤的发生和发展。4.2二级结构特征分析4.2.1预测方法在探究mRNA上m5C修饰的二级结构特征时，准确预测mRNA的二级结构是关键的第一步。目前，RNAfold是应用广泛且经典的预测工具，其原理基于最小自由能（MFE）算法。RNA分子的二级结构形成是一个热力学过程，RNAfold通过计算不同碱基配对方式下的自由能变化，寻找自由能最小的结构，以此作为最可能的二级结构预测结果。在实际操作中，使用RNAfold预测mRNA二级结构的流程较为简便。以在线工具为例，用户只需将待分析的mRNA序列复制粘贴到RNAfoldwebserver的文本框中，点击提交按钮，系统便会进入计算程序。在计算过程中，RNAfold会考虑mRNA序列中碱基之间的相互作用，包括氢键形成、碱基堆积力等因素，从而模拟RNA分子在自然状态下的折叠情况。计算完成后，结果会以多种形式呈现，其中最直观的是二级结构图形，以点括号表示法（dot-bracketnotation）和可视化图形展示碱基配对关系和结构形状。点括号表示法中，配对的碱基用括号表示，未配对的碱基用点表示，通过这种方式可以清晰地展示mRNA二级结构中的茎-环结构、发夹结构等特征。除了RNAfold，一些新兴的预测工具也在不断发展，它们结合了深度学习等先进技术，为mRNA二级结构预测提供了新的思路和方法。DeepRNA是一种基于深度学习的RNA二级结构预测工具，它利用卷积神经网络（CNN）和循环神经网络（RNN）对大量已知二级结构的RNA序列进行学习，从而构建预测模型。与传统的基于最小自由能算法的工具相比，DeepRNA能够更好地捕捉RNA序列中的复杂模式和长程相互作用，在一些复杂RNA结构的预测上表现出更高的准确性。然而，这些新兴工具在应用中也存在一定的局限性，例如需要大量的训练数据和较高的计算资源，模型的可解释性相对较差等。在实际研究中，通常会结合多种预测工具的结果进行综合分析，以提高预测的可靠性。4.2.2结构特点通过对大量mRNA二级结构的分析，发现m5C修饰位点呈现出独特的结构特点，对mRNA的结构和功能产生重要影响。m5C修饰位点常常倾向于形成小颈环结构，这一结构在mRNA的功能调控中扮演着关键角色。中山大学张锐教授团队的研究发现，NSUN2依赖的m5C位点通常位于一个小颈环结构的底部。这种小颈环结构为m5C修饰提供了特定的微环境，可能影响m5C修饰的发生和功能。小颈环结构中的碱基配对和空间构象使得m5C修饰位点更容易被甲基转移酶识别和作用，从而促进m5C修饰的形成。小颈环结构也可能影响m5C修饰对mRNA功能的调控。研究表明，小颈环结构中的m5C修饰可以通过改变mRNA与蛋白质的相互作用，影响mRNA的稳定性和翻译效率。在某些情况下，m5C修饰的小颈环结构可以与特定的RNA结合蛋白相互作用，增强mRNA的稳定性，延长其半衰期；而在另一些情况下，这种相互作用可能会抑制mRNA的翻译过程，调节蛋白质的合成速率。mRNAm5C修饰位点的小颈环结构与tRNA底物的结构高度相似，这一发现揭示了mRNAm5C修饰序列结构特征的由来，也暗示了它们在功能上的某种联系。tRNA是蛋白质合成过程中的重要参与者，其特定的结构对于识别密码子、转运氨基酸至关重要。mRNAm5C修饰位点形成的小颈环结构与tRNA底物结构的相似性，可能意味着它们在某些功能上存在共性。在翻译过程中，mRNA上m5C修饰位点所在的小颈环结构可能像tRNA一样，参与与核糖体、翻译因子等的相互作用，从而影响mRNA的翻译效率和准确性。这种结构上的相似性还可能反映了m5C修饰在进化过程中的保守性和重要性。从进化角度来看，mRNA和tRNA在结构和功能上的这种相似性可能是在长期的进化过程中逐渐形成的，以适应生物体内复杂的基因表达调控需求。4.3案例分析：某生物mRNA上m5C修饰序列结构特征以人类HeLa细胞mRNA为例，对其m5C修饰位点的序列和结构特征进行深入剖析。在序列特征方面，通过高置信度鉴定技术，确定了大量m5C修饰位点，并对这些位点周围的序列进行基序分析。研究发现，许多m5C修饰位点倾向于存在于“C/G-C/G-m5C-G/A-G-G-G”基序中。如在某一关键基因的mRNA序列中，其m5C修饰位点周围呈现出典型的这种基序模式，紧邻m5C位点上游的两个碱基分别为鸟嘌呤（G）和胞嘧啶（C），下游第一个碱基为腺嘌呤（A），随后紧跟三个鸟嘌呤（G）。这种基序模式并非偶然出现，而是在多个不同功能的mRNA中都有一定程度的体现，暗示了该基序可能为m5C修饰的发生提供了特定的序列环境，或者与甲基转移酶的识别和催化过程密切相关。从二级结构特征来看，利用RNAfold等预测工具对含有m5C修饰位点的mRNA区域进行二级结构预测，发现m5C修饰位点常常位于茎-环结构中，且多处于小颈环结构的底部。在HeLa细胞中编码细胞周期调控蛋白的mRNA中，m5C修饰位点所在区域形成了明显的茎-环结构，m5C修饰位点正好处于小颈环结构的底部位置。这种结构特点与m5C修饰的功能紧密相关。小颈环结构为m5C修饰提供了特殊的微环境，可能影响m5C修饰对mRNA稳定性和翻译效率的调控。由于小颈环结构的存在，使得m5C修饰位点更容易与相关的RNA结合蛋白相互作用。一些RNA结合蛋白能够特异性地识别m5C修饰位点所在的小颈环结构，从而增强mRNA的稳定性，防止其被核酸酶降解。在翻译过程中，小颈环结构中的m5C修饰可能会影响核糖体与mRNA的结合效率，进而调控翻译起始的速率。如果小颈环结构中的m5C修饰发生改变，可能会导致mRNA与核糖体的结合能力下降，使得翻译起始过程受阻，从而影响蛋白质的合成。在小鼠肝脏组织mRNA中，m5C修饰位点的序列和结构特征也展现出独特之处。在序列方面，虽然也存在部分位点符合“C/G-C/G-m5C-G/A-G-G-G”基序，但相较于人类HeLa细胞，其出现的频率和具体分布存在差异。在某些与肝脏代谢功能密切相关的mRNA中，发现了一些新的基序模式。在编码参与脂肪酸代谢关键酶的mRNA中，m5C修饰位点周围呈现出“G-m5C-A-A-C”的基序，这种基序在小鼠肝脏组织mRNA中具有一定的特异性，可能与肝脏独特的代谢功能和基因调控需求相关。从二级结构上看，小鼠肝脏组织mRNA的m5C修饰位点同样倾向于位于茎-环结构中，但在结构的具体形状和稳定性上与人类HeLa细胞有所不同。在小鼠肝脏组织中参与葡萄糖代谢的mRNA中，m5C修饰位点所在的茎-环结构相对更加稳定，碱基配对更加紧密。这种结构上的差异可能导致m5C修饰在小鼠肝脏组织mRNA中的功能表现与人类HeLa细胞有所不同。由于茎-环结构的稳定性较高，m5C修饰对mRNA稳定性的影响可能更为显著，使得该mRNA在肝脏细胞中的半衰期更长，从而保证相关基因能够持续稳定地表达，以维持肝脏正常的葡萄糖代谢功能。在翻译过程中，稳定的茎-环结构可能会对核糖体的移动产生一定的阻碍作用，进而影响翻译延伸的速率，使得蛋白质的合成过程更加精细地受到调控。五、m5C修饰与生物过程及疾病的关联5.1m5C修饰在胚胎发育中的作用在胚胎发育这一复杂而有序的过程中，m5C修饰扮演着不可或缺的角色，其对母源mRNA稳定性的调控以及在胚胎发育各个阶段的动态变化，深刻影响着胚胎的正常发育进程。以斑马鱼的胚胎发育研究为例，中国科学院北京基因组研究所杨运桂研究组、中国科学院动物研究所刘峰研究组和复旦大学麻锦彪研究组展开了深入合作。研究团队凭借卓越的科研技术，首次精心绘制了斑马鱼早期胚胎发育过程中RNAm5C甲基化图谱。通过对测序数据的细致分析，一个关键发现浮出水面：在母源-合子转换过程中，m5C修饰的母源mRNA比未修饰的mRNA展现出更高的稳定性。为了进一步探究其中的分子机制，研究团队运用Oligo-pulldown结合质谱技术，成功发现了斑马鱼中新的m5C结合蛋白Ybx1。借助ITC、蛋白解构技术等先进手段，深入研究发现Ybx1倾向于通过其CSD结构域关键残基Trp45与m5C的π-π相互作用，精准识别m5C修饰的mRNA。在功能实验中，当ybx1缺失时，早期胚胎发育出现了严重阻滞，停滞在囊胚期和原肠期。研究团队并未止步于此，他们通过免疫共沉淀、GST-pulldown技术和高通量测序技术，鉴定出poly(A)结合蛋白Pabpc1a为Ybx1的互作蛋白。有趣的是，pabpc1a缺失的胚胎表现出与ybx1缺失胚胎相似的表型，且Ybx1靶向结合的mRNA的稳定性显著降低。最终，利用单基因报告系统有力地证明了Ybx1和Pabpc1a以m5C依赖的方式维持母源mRNA的稳定性，从而最终促进斑马鱼早期胚胎发育进程。这一系列研究成果清晰地揭示了m5C修饰通过其结合蛋白Ybx1招募Pabpc1a维持母源mRNA稳定性，进而保证母源-合子转换有序进行的分子机制，为全面解析斑马鱼胚胎发育过程中mRNA命运调控及细胞分化提供了重要依据。果蝇的胚胎发育研究同样为m5C修饰的作用提供了有力证据。中山大学张锐教授团队在研究中，对果蝇进行了深入探究。通过对果蝇不同发育时期的样本进行mRNABS-Seq，构建了mRNAm5C图谱。研究发现，在果蝇的卵母细胞和胚胎发育早期，m5C高度富集。随着母源-合子过渡的发生，m5C位点数目和甲基化水平急剧下降。为了深入研究m5C修饰在果蝇胚胎发育中的功能，团队进行了敲除实验。当敲除果蝇中的NSUN2蛋白时，令人惊讶的是，96%的母源mRNAm5C被清除，这一结果直接导致果蝇胚胎发育迟缓。这一实验结果充分暗示了m5C在果蝇早期胚胎发育中的重要功能，表明m5C修饰对于维持果蝇胚胎发育的正常进程至关重要。从进化的角度来看，m5C修饰在不同物种胚胎发育中的作用具有一定的保守性。斑马鱼和果蝇作为两种不同的模式生物，在胚胎发育过程中都展现出m5C修饰对母源mRNA稳定性的重要调控作用，以及m5C修饰在胚胎发育早期的动态变化特征。这表明m5C修饰在胚胎发育中的功能可能在进化过程中被保留下来，以确保不同物种胚胎发育的正常进行。这种保守性也反映了m5C修饰在胚胎发育中的重要性和普遍性，为进一步研究m5C修饰在其他物种胚胎发育中的作用提供了重要的参考和启示。5.2m5C修饰在肿瘤发生发展中的机制在肿瘤发生发展的复杂进程中，m5C修饰扮演着至关重要的角色，其背后涉及多种关键蛋白的异常以及一系列复杂的分子机制。NSUN2作为m5C修饰的关键甲基转移酶，其异常表达与肿瘤的发生发展密切相关。在食管鳞状细胞癌（ESCC）的研究中，中山大学肿瘤防治中心华南肿瘤学国家重点实验室研究团队通过RNA-BS、RNA-seq等技术揭示了其中的奥秘。研究发现，ESCC肿瘤中的m5C甲基化水平显著增加，而这一变化是由于m5C甲基转移酶NSUN2的过表达所致。进一步研究表明，NSUN2的异常表达受E2F转录因子1（E2F1）的正向调控，高表达的NSUN2预示着ESCC患者生存率低。从分子机制层面深入探究，NSUN2通过增强生长因子受体结合蛋白2（GRB2）的m5C修饰，在新型m5C介质蛋白lin-28同源物B（LIN28B）的介导下，使得GRB2mRNA更加稳定。这种稳定性的增强，使得GRB2上调，进而增加了PI3K/AKT和ERK/MAPK信号的激活，最终促进了ESCC的发生和进展。在NSUN2敲低的ESCC细胞中，GRB2的m5C修饰水平降低，mRNA稳定性下降，PI3K/AKT和ERK/MAPK信号通路受到抑制，细胞的增殖、迁移和侵袭能力也明显减弱。这一系列研究结果清晰地表明，NSUN2介导的m5C修饰通过对GRB2mRNA稳定性的调控以及相关信号通路的激活，在ESCC的发生发展中发挥着关键作用。在乳腺癌的研究中，m5C修饰同样展现出重要影响。NSUN2通过m5C修饰调控靶基因HGH1的表达，增强HGH1mRNA的稳定性。稳定的HGH1mRNA能够持续高效地翻译出相关蛋白，这些蛋白参与了乳腺癌细胞的多种生物学过程，如细胞增殖、迁移和侵袭等，从而促进了乳腺癌的发展。当通过实验手段抑制NSUN2的活性或降低其表达水平时，HGH1mRNA的m5C修饰减少，稳定性降低，乳腺癌细胞的增殖和迁移能力也随之下降。这进一步证明了m5C修饰在乳腺癌发生发展中的重要调控作用。m5C修饰还与肿瘤的耐药性密切相关。在EGFR突变型非小细胞肺癌中，异常m5C高甲基化通过NSUN2/YBX1/QSOX1轴介导对吉非替尼的固有耐药性。研究人员发现，RNAm5C高甲基化和NSUN2与EGFR-TKI的固有耐药性显著相关。NSUN2过表达导致吉非替尼耐药和肿瘤复发，而NSUN2基因干扰则导致肿瘤消退，并在体外和体内克服了对吉非替尼的固有耐药性。深入研究发现，NSUN2甲基化喹呤巯基氧化酶1（QSOX1）编码序列区域，通过m5CreaderY-box结合蛋白1（YBX1）增强QSOX1翻译。QSOX1表达的增加，使得肿瘤细胞对吉非替尼的耐药性增强，从而影响了肿瘤的治疗效果。这一发现揭示了m5C修饰在肿瘤耐药性中的关键作用，为克服肿瘤耐药提供了新的研究方向。5.3其他疾病中m5C修饰的研究进展在动脉硬化领域，m5C修饰的研究尚处于起步阶段，但已展现出其潜在的重要性。动脉粥样硬化作为动脉硬化的常见类型，是心血管疾病的主要病理基础。虽然目前关于m5C修饰与动脉粥样硬化直接关联的研究相对较少，但从相关的分子机制和细胞生物学角度来看，m5C修饰可能参与其中。mRNAm5C修饰与mRNA的稳定性、剪接和核质穿梭等过程密切相关。在动脉粥样硬化的发生发展过程中，涉及到多种细胞的功能异常，如血管内皮细胞、平滑肌细胞和巨噬细胞等。这些细胞中mRNA的正常功能对于维持血管的稳态至关重要。若mRNA的m5C修饰发生异常，可能会导致相关基因表达失调，进而影响细胞的正常功能。血管内皮细胞中某些关键基因mRNA的m5C修饰异常，可能会影响细胞的屏障功能和对炎症因子的反应，促进炎症细胞的黏附和浸润，为动脉粥样硬化的发生埋下隐患。平滑肌细胞中与增殖和迁移相关基因mRNA的m5C修饰改变，可能会导致平滑肌细胞的异常增殖和迁移，参与动脉粥样硬化斑块的形成。虽然目前还缺乏直接的证据表明m5C修饰在动脉粥样硬化中的具体作用机制，但随着研究的深入，有望揭示m5C修饰在这一疾病中的关键作用，为动脉粥样硬化的防治提供新的靶点。在自身免疫性疾病方面，m5C修饰的研究已取得了一些有价值的成果。辅助性T细胞17（Th17）是CD4+辅助T细胞的一个亚群，在自身免疫性疾病中发挥着关键作用，其功能受到多种因素的调控，其中m5C修饰的作用逐渐被揭示。中国科学院北京基因组研究所（国家生物信息中心）杨运桂、杨莹，中国科学院生物物理研究所杨鹏远和中国科学院动物研究所李伟作为共同通讯作者在NatureCommunications期刊发表的研究论文表明，在小鼠CD4+T细胞中，RNA5-甲基胞嘧啶（m5C）甲基转移酶Nsun2缺失会特异性抑制Th17细胞分化，并减轻Th17细胞诱导的结肠炎发病。从机理上讲，RoRγt能将Nsun2招募到其靶标（包括Il17a和Il17f）的染色质区域，导致转录偶联的m5C形成，从而增强mRNA的稳定性。通过构建葡聚糖硫酸钠盐（DSS）诱导和CD4+CD25-CD45RBhiT细胞过继性两种小鼠结肠炎模型，发现Nsun2缺失可减轻Th17细胞诱导的结肠炎进展，并且该表型能够通过回补Nsun2的靶标基因产物IL-17A/IL-17F蛋白逆转。这一研究成果不仅揭示了m5C修饰在Th17细胞分化和功能调控中的重要作用，也为自身免疫性疾病的治疗提供了潜在靶点。在系统性红斑狼疮、硬皮病、银屑病和类风湿性关节炎等其他自身免疫性疾病中，虽然m5C修饰的具体研究相对较少，但鉴于Th17细胞在这些疾病中的重要作用以及m5C修饰对Th17细胞的调控作用，可以推测m5C修饰可能在这些疾病的发生发展过程中也扮演着重要角色。未来的研究可以进一步深入探讨m5C修饰在不同自身免疫性疾病中的作用机制，为这些疾病的治疗提供新的思路和方法。六、结论与展望6.1研究成果总结本研究围绕mRNA上m5C修饰展开，通过运用先进的技术手段和严谨的分析方法，在多个关键方面取得了具有重要意义的研究成果。在高置信度鉴定mRNA上m5C修饰方面，成功应用超快速亚硫酸氢盐测序（UBS-seq）和m5C-TAC-seq等新型技术。UBS-seq通过使用高浓度亚硫酸盐和高反应温度，显著缩短了反应时间，减少了RNA降解，同时提高了C到U的转化效率，有效降低了检测背景和假阳性率，在HeLamRNA的检测中发现了数千个m