突变Aerolysin单体在检测RNA碱基序列及RNA修饰方面的应用

摘要

本发明涉及基因检测技术领域，具体涉及突变Aerolysin单体在检测RNA碱基序列及RNA修饰方面的应用。本发明人证实了使用突变Aerolysin纳米孔检测mRNA片段中的m6A修饰，可以从混合物中将甲基化与未甲基化的RNA链区分开。因此，本发明提供了一种方法表征靶标RNA以及RNA上的碱基修饰，其包括但不限于m6A，m1A,m7G,poly(A)尾，U尾,该方法包括：(a)利用CRISPR‑Cas13a和crRNA特异性识别包含碱基修饰的目标RNA序列，并利用RNA内切酶将长链RNA序列切割成短的包含修饰碱基的短RNA序列。修饰碱基在序列上的位置特定，如核苷酸长度为14，该RNA链从3’端数，第4/11/14个为修饰的碱基。

说明书

技术领域

本发明涉及基因检测技术领域，涉及纳米孔检测RNA修饰方法，具体为突变Aerolysin单体在检测RNA 碱基序列及RNA修饰方面的应用。

技术背景

目前开发应用范围跨度广、快速且廉价的核酸鉴别技术是必要的。现有技术检测慢、价格昂贵的主要原因是它们需要扩增技术以产生大量体积的核酸，并且需要大量的用于信号检测的专业荧光化学试剂。纳米孔作为直接的电生物传感器，可降低所需核苷酸和试剂的数量，在快速、廉价地实现核酸检测方面具有巨大的潜力。

纳米孔传感技术通过检测待测物分子与孔道的相互作用产生的特征电流阻断信号来实现单分子水平的分析。当在纳米孔上施加电压时，当存在待测物，例如核酸的情况下，会产生电流的变化。通过分析物的特征电流标志，例如阻断电流大小和持续时间来揭露分析物的特性。使用纳米孔进行核酸检测的两个基本组分是(1)控制通过孔的核酸移动和(2)在核酸聚合物通过孔移动时鉴别核苷酸。

CRISPR是目前流行的强大的基因编辑工具，一些Cas酶例如Cas9, 靶向DNA，而另一些Cas酶，例如Cas 13a（又称C2c2）可以与CRISPR RNA(cr-RNA)重编程，以特异性识别RNA。识别其RNA靶标后，激活Cas13a参与附近非靶标RNA的“附带”切割。该方法可以通过设计cr-RNA序列，根据纳米孔对非特异性切割的RNA片段的读出，判断靶RNA的存在，且具有较高的灵敏度。

目前，纳米孔技术可以检测的合适的核酸序列通常长度较短，因此从长RNA序列上截取适合纳米孔分析的包含修饰位点的目标片段的样品制备方法是十分重要的。Rnase H酶是一种RNA内切酶，其不能切割单链DNA或单链RNA，但可以切割DNA与RNA形成的双链。因此当已知靶RNA的序列，拟利用纳米孔技术进行RNA修饰检测时，可先设计DNA序列，与包含修饰位点的目标RNA片段的周围序列互补配对成双链，经Rnase H酶切割可得适合纳米孔检测长度的，包含修饰位点的目标RNA片段。

RNA的表观遗传修饰是RNA调节基因表达的基础，RNA经过化学修饰后会影响其活性、定位和稳定性，常见的RNA表观遗传修饰包括m

发明内容

本发明所要解决的技术问题是：提供突变Aerolysin单体在检测RNA 碱基序列及RNA修饰方面的应用，发现突变Aerolysin单体制备的纳米孔可以检测mRNA片段中的m

为了解决上述技术问题，本发明采用的技术方案如下：提供一种突变Aerolysin单体在检测RNA 碱基序列及RNA修饰方面的应用；所述突变Aerolysin单体包含T232K/K238Q两个氨基酸位点的突变。

优选的，所述突变Aerolysin单体还包含N226Q、G270K、S228K、T232K、K238Q中至少一个的氨基酸位点的突变。

其中，所述Aerolysin单体在检测RNA 碱基序列及其修饰方面的应用，其特征在于，所述RNA修饰为RNA上的碱基修饰，其包括但不限于m

所述突变Aerolysin单体在检测RNA 碱基序列的具体方法如下：(a) 设计crRNA序列，利用CRISPR-Cas13a和crRNA特异性识别包含碱基修饰的靶RNA序列，并将长链RNA序列切割成短的包含修饰碱基的短RNA序列，（c）设计DNA序列与靶RNA片段的周围序列互补配对；(d) 用核酸内切酶切割DNA-RNA双链，得到靶RNA序列；（e）使用Aerolysin单体制备的纳米孔检测靶RNA片段的信号，从而确认靶RNA的存在及其修饰。

其中，上述步骤（e）使用Aerolysin单体制备的纳米孔检测靶RNA片段的信号具体方法为：使所述靶RNA片段相对于Aerolysin单体制备的纳米孔中移动，测量所述短RNA的穿孔电流信号，根据数学变换解析出该电流信号的瞬时频率从而表征短RNA片段。

本发明的有益效果如下：

使用本发明所述突变Aerolysin单体检测RNA 碱基序列及RNA修饰，使用Aerolysin单体制备的纳米孔检测靶RNA片段的信号，其所述孔道已具有控制核酸移动速度的功能，无需现有纳米孔RNA测序中所使用的聚合酶的参与，无需特殊化学修饰或荧光标记，通过纳米孔电信号即可在单分子水平上对样品中是否含有目的RNA序列进行检测以及对修饰的RNA链进行区分。

附图说明

图1为本发明的CRISPR-Cas13a检测靶RNA存在示意图。

图2为本发明的Rnase H切割DNA-RNA双链获取目标RNA链示意图。

图3为本发明的Aerolysin突变孔用于检测RNA修饰示意图。包括T232KK238Q,N226QK238Q, N226QG270K, N226QS228K Aerolysin突变孔。

图4A为本发明的实施例1所述RNA样品相对于孔的移动的电流信号示意图，RNA1序列：5’- CAGACUCUGUAGCG -3’和RNA2序列：5’- CAGm

具体实施方式

本文所提到的剪切酶均指CRISPR-Cas 13a，其用于纳米孔检测RNA修饰的样品前处理，特异性识别靶RNA序列并剪切，而非指纳米孔测序中使用的聚合酶，本发明中纳米孔检测方法采用无酶的方案，突变Aerolysin纳米孔实现对RNA的控速。下述是本发明技术方案中所用到的几个基本原理

一、基于CRISPR-Cas13的纳米孔检测靶RNA

本发明提供了用CRISPR-Cas 13a特异性识别样品中是否存在靶RNA序列的纳米孔检测的方法。包含靶向靶RNA的crRNA的设计与合成，CRISPR-Cas 13a系统的特异性识别和剪切。

靶向含有碱基修饰位点的目标RNA的crRNA的设计原则如下：

（1）crRNA包括间隔序列（spacer）和直接重复序列（Direct Repeat, DR），序列格式为：5’-与Cas 13a蛋白结合的直接重复序列-crRNA间隔序列-3’，直接重复序列需要根据Cas13 a蛋白来确定，使其能够与选定的Cas 13a蛋白匹配并结合；

（2）crRNA的间隔序列长度为22~28个碱基序列；

（3）crRNA的直接重复序列大于24个碱基序列；

（4）crRNA的直接重复序列应包含茎环（stem loop）结构；

（5）在crRNA的间隔序列中间为seed区，与靶序列结合时不能发生错配。

CRISPR-Cas 13系统特异性识别靶RNA。所述CRISPR-Cas 13a核酸酶结合crRNA可以被目标核酸序列特异性激活，并剪切该序列。CRISPR-Cas核酸酶还可以包括LbaCas13,LbuC13a, LwaCas13a, AspCas13a, BzoCas13a, CcaCas13a, AspCas13b, BzoCas13b,CcaCas13b, PsmCas13b, PinCas13b, PbuCas13b, PguCas13b, PigCas13b, PsaCas13b,RanCas13b, PspCas13b, EsCas13d, RspCas13d。该酶可以包括野生型、改造过的、密码子优化过的、进化过的、嗜热的、工程化的和/或一种以上Cas蛋白的混合物。CRISPR-Cas核酸酶优选地为LwaCas 13a。

二 Rnase H酶剪切获得包含修饰位点的目标RNA链

本发明提供了样品制备的方法，用于纳米孔检测RNA的碱基修饰。该方法包含设计DNA序列与目标RNA片段的周围序列互补配对，Rnase H剪切DNA-RNA双链从而得到包含修饰位点的适合纳米孔检测的目标RNA链。Rnase H起核酸内切酶的作用需要二价阳离子的参与，其包括Mg

三、表征靶RNA

本发明提供了表征靶RNA的方法。已知靶RNA序列和碱基修饰发生的位置，截取碱基长度为14的RNA链，使该碱基修饰位点位于从3’ 端开始数的第11个碱基。甲基化和未甲基化的RNA截取方式相同。获得该混合物RNA后使用Aerolysin突变孔表征。

本发明样品制备非常简单，因为纳米孔能够检测单个RNA分子，不需扩增靶RNA或互补核酸。本发明通常不涉及聚合酶链式反应或逆转录。这大大减少了表征靶RNA的工作流程。

四、膜

本发明的纳米孔需要插入膜中构成一个体系。膜优选地是两亲层。两亲层是由两亲性分子形成的两亲膜，所述两亲性分子如磷脂，具有亲水性和亲脂性两种特性。两亲分子可以是合成的或天然存在的。

下面结合说明书附图对本发明作进一步说明。

实施例1 RNA特征序列的检测

使用Aerolysin单体制备的纳米孔对RNA特征序列的表征：，待测RNA序列为天然牛催乳素NCBI NM_173953.2的3’UTR序列的一段片段，其RNA序列和crRNA序列见表1。Cas13a蛋白 (2 mg/mL, 4 μL)使用前用122.5 μL 储液重悬。将1 μL待测RNA样品，2 μL酶缓冲溶液（400 mM Tris pH 7.4），9.6 μL DEPC处理水，2 μL Cas13a蛋白溶液，1 μL cr-RNA(10ng/μL)，SUPERnase核糖核酸抑制剂1 μL，MgCl

表1 待测靶RNA, crRNA和DNA序列

实施例2 RNA修饰的检测

Rnase H剪切目标RNA用于纳米孔RNA碱基修饰检测：第一步，退火。将待测RNA样品30 μL（约100 μM，含无碱基修饰和有碱基修饰序列）与30 μL的DNA1（100 μM）和30 μL的DNA2（100 μM）混合，95°C变性5min，后慢慢冷却到室温，形成DNA-RNA双链。第二步，将第一步得到的样品90 μL（约30 μM）与22.5 μL的5×Rnase H酶缓冲溶液（500 mM Tris-HCl；750mM KCl；30 mM MgCl

使用T232K/K238Q Aerolysin突变单体制备的纳米孔检测上述制备获得的靶RNA片段的碱基修饰：制备T232K/K238Q Aerolysin纳米孔，用于检测RNA短片段的T232K/K238QAerolysin突变孔如图3所示，对Aerolysin纳米孔β桶内232位点和238位点氨基酸进行突变，可以延长RNA链在孔道内的停留时间，根据甲基化修饰和未甲基化修饰的RNA链与纳米孔道内氨基酸残基的相互作用不同产生的纳米孔电信号不同来进行检测区分。具体原理如下：在室温下用胰蛋白酶酶切Proaerolysin蛋白6小时，得到有活性的Aerolysin蛋白。实验所用的缩醛树脂检测池由分别称为cis端和trans端的两个检测池组成，两个检测池中分别含有1mL的缓冲溶液（1.0 M /3.0 M KCl，10 mM Tris，pH 8.0），该缓冲溶液由DEPC处理水配制以防止RNA的降解，磷脂双分子层趴附在检测池中间50 µm的小孔上，在cis端检测池中加入一定量的T232K/K238Q Aerolysin突变单体蛋白溶液，T232K/K238Q Aerolysin突变单体自组装成七聚体并插入到磷脂双分子层上形成纳米孔道，用于控制RNA的移动速度，所有的纳米孔电化学分析实验温度均控制在20 ± 2 ℃。将经Rnase H剪切得到的目标RNA样品（体积为40 μL，浓度为20 μM）加入到cis端检测池中，检测不同电压下的电流信号。使用膜片钳放大器(Axopatch 200B) 放大和采集Aerolysin纳米孔道内的电流信号，并通过模数转换器(Digidata 1440A)转换成数字信号。采样率为 100 kHz，滤波频率为5 kHz; 通过Clampex 10.4软件(Molecular Devices, Forest City, CA)获取电流信号并使用MOSAIC和Origin软件对所得数据进行统计分析。仪器所采集到的每个数据点的电流值，包括的内容有：所处时间点下纳米孔道的电流基线、电流阻断以及电流噪音的综合数据。由于某些信号阻断时间较短，处理时被认为是噪声而丢失，可以通过计算阻断数据中的所有极小极大值点来进行信号的识别。

采用上述方法检测甲基化和未甲基化RNA混合物的电流信号，其相对于孔的移动的电流轨迹如图4A所示，若样品中含有特异性的RNA链，则可以用Cas13a酶切割，产生纳米孔阻断电流信号，若样品中不含特异性的RNA链，则不会被Cas13a酶切割，不会产生纳米孔阻断电流信号。

实施例3本实例描述用于纳米孔RNA特征序列表征和将其剪切成5个核苷酸长度（RNA3：5’-GGACU-3’和RNA4: 5’-GGm

采用RNA序列为m

表2 待测靶RNA, crRNA和DNA序列

Rnase H剪切目标RNA用于纳米孔RNA碱基修饰检测：第一步，退火。将待测RNA样品20 μL（约100 μM，含无碱基修饰和有碱基修饰序列）与20 μL的DNA3（100 μM）混合，95°C变性5min，后慢慢冷却到室温，形成DNA-RNA双链。第二步，将第一步得到的样品90 μL（约30 μM）与10 μL的5×Rnase H酶缓冲溶液（500 mM Tris-HCl；750 mM KCl；30 mM MgCl

采用T232K/K238Q Aerolysin突变纳米孔检测上述制备获得的靶RNA片段的碱基修饰：制备T232K/K238Q Aerolysin纳米孔，包括将质粒Proaerolysin进行T232K/K238Q两个氨基酸位点的突变并将其表达和纯化，用于检测RNA的T232K/K238Q Aerolysin突变孔如图3所示，用于检测RNA的T232KK238Q Aerolysin突变孔如图3左上所示，第238位点带正电荷的赖氨酸突变成不带电荷的谷氨酰胺，延长核酸的过孔时间，第232位点的苏氨酸突变成带正电荷的赖氨酸，增加灵敏位点。其他突变孔如N226QK238Q/同样起着延长时间，灵敏探测的作用。在室温下用胰蛋白酶酶切Proaerolysin蛋白6小时，得到有活性的Aerolysin蛋白。实验所用的缩醛树脂检测池由分别称为cis端和trans端的两个检测池组成，两个检测池中分别含有1mL的缓冲溶液（1.0 M /3.0 M KCl，10 mM Tris，pH 8.0），该缓冲溶液由DEPC处理水配制以防止RNA的降解，磷脂双分子层趴附在检测池中间50 µm的小孔上，在cis端检测池中加入一定量的T232K/K238Q Aerolysin突变单体蛋白溶液，T232K/K238QAerolysin突变单体自组装成七聚体并插入到磷脂双分子层上形成纳米孔道，用于控制RNA的移动速度，所有的纳米孔电化学分析实验温度均控制在20 ± 2 ℃。将经Rnase H剪切得到的目标RNA样品（体积为20 μL，浓度为50 μM）加入到cis端检测池中，检测不同电压下的电流信号。使用膜片钳放大器(Axopatch 200B) 放大和采集Aerolysin纳米孔道内的电流信号，并通过模数转换器(Digidata 1440A)转换成数字信号。采样率为 100 kHz，滤波频率为5 kHz; 通过Clampex 10.4软件(Molecular Devices, Forest City, CA)获取电流信号并使用MOSAIC和Origin软件对所得数据进行统计分析。仪器所采集到的每个数据点的电流值，包括的内容有：所处时间点下纳米孔道的电流基线、电流阻断以及电流噪音的综合数据。由于某些信号阻断时间较短，处理时被认为是噪声而丢失，可以通过计算阻断数据中的所有极小极大值点来进行信号的识别。

采用上述方法检测甲基化和未甲基化RNA混合物的电流信号，其相对于孔的移动的电流轨迹如图4B所示，与未甲基化的RNA相比，甲基化的RNA阻断程度更大，阻断时间更长，可能是与孔道内R220区域的相互作用更强。

以上显示和描述了本发明的基本原理、主要特征和本发明的优点。本行业的技术人员应该了解，本发明不受上述实施例的限制，上述实施例和说明书中描述的仅为本发明的优选例，并不用来限制本发明，在不脱离本发明精神和范围的前提下，本发明还会有各种变化和改进，这些变化和改进都落入要求保护的本发明范围内。本发明要求保护范围由所附的权利要求书及其等效物界定。