一种全转录组水平RNA结构检测方法及其应用

一种全转录组水平rna结构检测方法及其应用
技术领域
1.本发明属于生物技术领域，具体涉及一种全转录组水平rna结构检测方法及其应用。


背景技术：

2.rna具有不同的功能，如：作为信使传递遗传信息，作为核酶催化反应等等。rna分子在其整个生命周期中并且在不同的亚细胞位置均受到精确调节。复杂且灵活的结构是rna分子的功能多样性和精细调节的核心。rna结构的错误折叠能够干扰可变剪接、翻译、rna修饰和编辑以及rna-蛋白质相互作用等过程，从而引起疾病。
3.rna结构检测方法利用了特异性修饰单链核苷酸的化学试剂。修饰位点能够干扰逆转录(rt)的进行，导致rt停止或突变，因此能够通过测序和生物信息学分析方法检测到修饰位点信息，从而获得rna结构信息。大多数试剂只能检测一个或两个碱基的结构信息；例如，硫酸二甲酯(dms)修饰单链胞嘧啶和腺嘌呤，乙二醛修饰单链鸟嘌呤、胞嘧啶和腺嘌呤，并且乙氧二羟丁酮修饰单链鸟嘌呤。引物延伸的选择性2-羟基酰化分析法(shape)试剂能够修饰单链区段内核糖的2'oh基团，并能获得所有四种核苷酸的结构信息。
4.全局rna结构检测研究已经揭示了功能性rna位点处往往存在结构差异，例如蛋白质和mirna结合位点，并且已有研究表明rna结构能够参与调节rna的剪接、翻译和降解过程。值得注意的是，几项研究已经表明了rna序列可以在体内与体外、在不同的亚细胞区间以及在胚胎发生的不同阶段形成不同的结构。实际上，细胞中的许多因素可以影响rna结构，包括ph、阳离子浓度、内源rna修饰(例如，甲基化、乙酰化)以及与蛋白质和/或其他rna的相互作用。因此，在其最相关的自然环境中研究rna结构对于揭示rna功能和调节机制至关重要。
5.然而，目前最先进的rna结构检测方法通常需要大量rna作为起始量，这会限制其实际应用。例如，icshape和structure-seq2的rna文库的构建需要大约107个细胞，这对于罕见的原代细胞和许多组织样品的生物学研究是难以实现的。因此，除了实验上易收集的斑马鱼早期胚胎和果蝇卵巢的一些研究之外，迄今为止的rna结构检测研究仅限于培养的细胞系。然而，细胞系中的细胞环境和由此生成的rna结构可能显著偏离原代样品，从而使得其结果不能真实反应细胞的功能状态。


技术实现要素：

6.为了解决这一障碍，我们开发了smartshape(small amount random rt icshape，小量随机rt icshape)，一种基于icshape方法改进的新型低起始量rna二级结构检测方法。
7.因此，
8.本发明第一方面，提供一种rna结构检测方法，其中，所述方法包括：
9.1、获得包含rna的样本；2、smartshape库准备；3、rna结构检测和分析，其中，所述步骤2smartshape库准备包括：(1)、rna修饰和制备；(2)rna逆转录，去除非修饰位点引起的
逆转录终止信号(premature rt stops)，和cdna富集。
10.优选的，所述rna结构检测方法的步骤2还包括(3)、接头连接，第二链合成，和扩增。更优选的，所述接头连接包括3’接头连接和5’接头连接。
11.优选的，所述背景逆转录终止信号由非rna修饰位点导致。更优选的，所述背景逆转录终止信号可能源于内源修饰(例如m1a修饰)、局部结构(例如g-四链体)，或者源于逆转录酶的随机脱落。
12.更优选的，采用核糖核酸酶(ribonuclease，rnase)消化去除背景逆转录终止信号，更优选的，采用rnase i消化去除背景逆转录终止信号。
13.优选的，所述逆转录(rt)引物序列为5
’-
nnnnnn-3’、5
’-
nnwnnwnn-3’、5
’-
ttttttttvn-3’。优选的，利用标记试剂对rna进行修饰，更优选的，所述标记试剂为细胞膜穿透性试剂，更优选的，所述标记试剂选用硫酸二甲酯(dms)、1-甲基-7-硝基靛红酸酐(1m7)、2-甲基烟酸咪唑化物-叠氮化物(nai-n3)或乙氧二羟丁酮；更为优选的，标记试剂选用2-甲基烟酸咪唑化物-叠氮化物(nai-n3)。
14.优选的，cdna富集采用磁珠进行富集，更优选的，所述磁珠为链霉亲和素磁珠，例如myone c1磁珠。
15.优选的，所述rna结构为rna二级结构。
16.优选的，所述rna为全长rna；进一步地，rna为转录组rna。可以是长链rna，例如mrna、lncrna、rrna等，也可以包含很多小rna，例如小于200nt的小rna，蛋白质结合rna，作为dicer作用底物的rna等等。
17.优选的，所述rna可来源于任意细胞、病毒等，优选的，所述细胞包括但不仅限于实验室培养的细胞系，活体细胞，原代细胞、哺乳动物的早期胚胎、细菌、真菌以及各种感染后的细胞，例如病毒、细菌、真菌等感染后的细胞，更优选的，所述活体细胞可以是任意的体细胞、生殖细胞，例如上皮细胞、真皮细胞、腺体细胞、血液来源的细胞、骨细胞、免疫细胞(t细胞、b细胞、nk细胞、巨噬细胞等等)、受精卵等等。
18.所述rna结构检测方法还包括利用计算管道对smartshape分数进行计算处理步骤。所述计算处理步骤包括：1)去除3'接头；2)去除重复的读段；3)去除分子标签；4)将干净的读段比对到rrna标准序列；5)将未比对到rrna序列的读段比对到基因组；6)使用icshape-pipe sam2tab将sam文件转换成.tab文件；7)使用icshape-pipe calcshapenocont计算smartshape分数。
19.优选的，所述步骤7)中通过在所有外显子上以滑动窗口方式对rt停止计数进行归一化和缩尾处理来计算smartshape分数，并且将覆盖率低于100的碱基的分数定义为空(null)。
20.更优选的，所述步骤7)中的参数为：-n nai_rep1.tab,nai_rep2.tab；-size chrnamelength.txt；-out reactivity.gtab；-ijf sjdblist.fromgtf.out.tab。
21.优选的，所述检测方法在文库扩增前不包括胶回收步骤。
22.优选的，所述计算管道的文库构建中不需要对照组来去掉背景信号。
23.优选的，所述rna结构检测方法中可以少至1ng的起始量rna(104～105个细胞)对rna结构进行检测。
24.本发明还提供一种上述rna结构检测方法的应用，所述应用包括根据上述检测方
法的结果评估细胞的功能状态，研究rna对早期发育的影响，癌症的发生和发展等等。
25.优选的，所述功能状态包括各种生理及异常状态，例如，细胞炎症、损伤、缺血、免疫应激状态、早期发育过程、感染、癌症增殖等等。更优选的，所述感染由病毒、细菌、真菌等引起。
26.优选的，所述细胞来源任意组织器官，例如皮肤系统、血液淋巴系统、免疫系统、心血管系统、消化系统、呼吸系统、泌尿系统、骨骼系统、生殖系统、神经系统等等。
27.优选的，所述细胞包括免疫细胞，例如b细胞、t细胞、nk细胞、巨噬细胞等。
28.优选的，所述应用不是疾病的诊断和治疗方法。
29.本发明还提供一种细胞功能状态的评估方法，所述评估方法包括利用上述任意的检测方法对细胞的rna结构进行检测，根据检测结果评估细胞的功能状态。
30.优选的，所述细胞功能状态是细胞炎症、损伤、缺血、免疫应激状态、早期发育过程、感染、癌症增殖等等，更优选的，所述感染由病毒、细菌、真菌等引起。
31.更优选的，所述细胞功能状态是细胞的免疫应激状态。例如免疫细胞的免疫应激状态。更进一步优选的，所述免疫细胞，包括例如b细胞、t细胞、nk细胞、巨噬细胞等。
32.本发明的有益技术效果在于：
33.1、本发明去除背景逆转录终止信号，降低了背景逆转录终止信号在结构分数计算中引起的假阳性信号，从而提高检测方法的准确性。
34.2、本发明采取了不同的文库构建策略，其中我们把随机rt和珠上单链dna文库构建相结合，大大减少了由多个纯化步骤引起的损失。
35.3、smartshape需要少至1ng的起始量rna(104～105个细胞)，使得能够以非常低的样品量进行体内细胞的rna结构分析，可将其应用于任意的细胞，例如罕见的原代细胞、哺乳动物早期胚胎以及患者活检样本。
36.4、我们应用smartshape来描述来自细菌感染模型小鼠的肠道巨噬细胞的全转录组rna二级结构，其中每个样品仅有100ng总rna作为起始量。我们揭示了免疫应激后两种巨噬细胞群之间rna结构的差异，其富含免疫应答相关基因，并提供了通过rna结构调节免疫应答的证据。
37.5、本发明smartshape是一种用于研究全转录组体内rna二级结构的有效、准确和稳健的方法，只需要非常少量的rna作为起始量。我们的方法整合了随机逆转录、rnase i消化和珠上文库构建，以提高文库构建的效率并产生准确的rna结构数据。本发明结果表明，smartshape通过先rnase i消化后磁珠富集成功地去除背景逆转录终止信号，并且即使没有dmso组作为对照，也实现了优于icshape的准确度。
38.6、鉴于本发明的方法对rna起始材料的最低要求，非常有希望将smartshape应用于研究rna结构在潜在的许多其他生物环境中所起的广泛作用。例如，母体rna降解对于早期发育至关重要，并且一些研究已经报道了rna结构在斑马鱼早期胚胎发生期间在母体rna降解中起调节作用。现有技术中由于样品量有限，哺乳动物早期胚胎中的rna结构组尚未被研究，而本发明可以通过smartshape来实现。另外，已知rbp结合的失调参与了许多癌症的发生和发展，smartshape可提供一种可行手段，通过使用来自临床的罕见活检样品从rna结构角度来研究这些失调。另外，当与富集(例如，通过反义寡核苷酸或蛋白质抗体)组合使用时，预期smartshape会有助于发现并功能验证基于rna结构的调控作用，这些rna包括低水
平表达的rna(如许多lncrna)、应激颗粒中的rna种类和由rbp结合的rna片段等等。
39.以上只是概括了本发明的一些方面，不是也不应该认为是在任何方面限制本发明。除非特别说明，本发明的实践将采取细胞生物学、细胞培养、分子生物学和免疫学等的传统技术。这些技术在以下文献中进行了详细的解释。例如：
40.1、xu,h.et al.notch-rbp-j signaling regulates the transcription factor irf8 to promote inflammatory macrophage polarization.nat immunol 13,642-650,doi:10.1038/ni.2304(2012)；
41.2、li,p.,shi,r.&zhang,q.c.icshape-pipe:a comprehensive toolkit for icshape data analysis and evaluation.methods 178,96-103,doi:10.1016/j.ymeth.2019.09.020(2020)；
42.3、bolger,a.m.,lohse,m.&usadel,b.trimmomatic:a flexible trimmer for illumina sequence data.bioinformatics 30,2114-2120,doi:10.1093/bioinformatics/btu170(2014)；
43.4、langmead,b.&salzberg,s.l.fast gapped-read alignment with bowtie2.nat methods 9,357-359,doi:10.1038/nmeth.1923(2012)；
44.5、dobin,a.et al.star:ultrafast universal rna-seq aligner.bioinformatics 29,15-21,doi:10.1093/bioinformatics/bts635(2013)；
45.6、pedregosa,f.et al.scikit-learn:machine learning in python.j mach learn res 12,2825-2830(2011)；
46.7、reuter,j.s.&mathews,d.h.rna structure:software for rna secondary structure prediction and analysis.bmc bioinformatics 11,129,doi:10.1186/1471-2105-11-129(2010)；
47.8、spitale,r.c.et al.structural imprints in vivo decode rna regulatory mechanisms.nature 519,486-490,doi:10.1038/nature14263(2015)。
48.本说明书提到的所有专利和出版物都是通过参考文献作为整体而引入本发明的。本领域的技术人员应认识到，对本发明可作某些改变并不偏离本发明的构思或范围。下面的实施例进一步详细说明本发明，不能认为是限制本发明或本发明所说明的具体方法的范围。
附图说明
49.图1：smartshape文库制备示意图；
50.图2：rna打断和3’dna接头连接步骤的优化，其中图2a为nai-n3修饰或未经修饰的hek293t总rna在不同打断条件下的产率和片段分布；图2b为三种不同结构的接头示意图，包括短接头、包含10碱基分子barcode的长接头、在长接头的5’端增加一个随机核苷酸的接头；图2c为circligase和t4 dna ligase在一个合成的dna分子的3’端连接接头的连接产物。
51.图3：smartshape中经rnase i消化对背景噪音的去除，其中图3a为rnase i消化和磁珠富集去除背景噪音的示意图；图3b为28s核糖体rna中一个已知的m1a修饰所在位点，图3c为在该m1a位点上游设计引物，背景逆转录信号检测；图3d为内源性m1a或m3u修饰位点已
知的m1a修饰位点处，dmso组和nai-n3组的逆转录终止信号差异；图3e为18s核糖体rna中的一段序列，从左至右分别代表了只用nai-n3组计算的smartshape值、用nai-n3组和dmso组计算的icshape值；图3f计算了18s核糖体rna两种shape值对应的roc曲线。
52.图4：rnase i消化能够有效去除背景信号，其中图4a是合成的rna序列和结构，图4b为分别在体外折叠两条合成的rna后，用nai-n3进行修饰后进行逆转录，对逆转录产物同时进行rnase i消化和磁珠富集时，去除m1a修饰引起的背景逆转录信号；图4c为dmso组建库流程；图4d为所有核糖体rna位点的dmso组和nai-n3组的逆转录终止信号差值分布，不同直线代表核糖体rna中所有已知内源修饰位点的终止信号差异均值；图4e为背景信号异常高的位点处，不同nai-n3文库中逆转录终止信号的分布。
53.图5：使用不同的起始量rna时smartshape的覆盖率和准确率，其中图5a为四种不同起始量smartshape文库和icshape文库在rps16转录本每个位点上的逆转录终止信号；图5b为四种不同起始量smartshape文库以及icshape文库在不同测序深度下检测到的高覆盖度转录本数目；图5c为四种不同起始量smartshape文库和icshape文库在每一步处理过程中对应的读段数目；图5d为四种不同起始量smartshape文库和icshape文库在18s和28s核糖体rna中的roc曲线图5e为四种不同起始量smartshape文库和icshape文库在xbp1结构元件的auc，对应该位点的shape值。
54.图6：不同的起始量smartshape文库具有很高的可重复性和文库复杂度，其中图6a为四种不同起始量(1ng、5ng、25ng和125ng)的smartshape文库和icshape文库shape值的相关性；图6b为四种不同起始量(1ng、5ng、25ng和125ng)的smartshape文库及icshape文库中每个转录本中具有shape值的位点在不同文库技术重复间的pearson相关性的分布；图6c为不同测序深度下四种不同起始量smartshape文库中，每个转录本平均逆转录终止信号的累积分布曲线。
55.图7：smartshape文库检测到与icshape类似的结构特点，其中图7a为smartshape和icshape文库在起始密码子上游30个碱基至下游100个碱基区间及终止密码子上游100个碱基至下游30个碱基区间，每个位点的平均shape值；图7b为四种不同起始量smartshape文库及icshape文库中a、u、g、c四种不同碱基的shape值的分布；图7c为smartshape和icshape文库在m6a修饰附近，每个位点的平均shape值；图7d为smartshape文库及icshape文库中不同rna种类或区域的gini指数的分布。
56.图8：利用smartshape检测小鼠体内肠道巨噬细胞rna结构，其中图8a为小鼠巨噬细胞的分离和rna二级结构检测流程图；图8b为两种巨噬细胞smartshape文库中高覆盖率转录本的数目，即在超过80％的位点覆盖度超过100的转录本数目；图8c为两种巨噬细胞smartshape文库和icshape文库在xbp1已知结构元件的auc。
57.图9：通过流式细胞仪基于免疫相关的基因mhcii、cd45、siglecf、cd11b、cd11c、cd64和ly6c分选出ly6c
lo
组织常驻巨噬细胞和ly6c
hi
促炎巨噬细胞。
58.图10：巨噬细胞smartshape数据的准确性，其中图10a为两种巨噬细胞smartshape文库和icshape文库对于srp rna的auc；图10b为对于rfam数据库中60个已知的rna结构，分别计算了两种巨噬细胞smartshape数据和小鼠胚胎干细胞icshape数据在每个结构中的roc曲线和对应的曲线下面积，并在图中展示了每个文库对应的auc的分布。
rna clean&concentrator-5柱用于纯化。
70.2、逆转录、rnase消化、富集和3'接头连接。
71.3.5μl的rt引物混合物(50μm 5
’-
nnnnnn-3’、50μm 5
’-
nnwnnwnn-3’和6μm5
’-
ttttttttvn-3’)和3μ的5
×
第一链缓冲液(life technologies)加入到8.5μl的生物素化的rna样品中。将样品加热至85℃持续5分钟，然后缓慢冷却至4℃(每秒0.1℃)以进行引物退火和弱片段化。向具有引物的rna提供0.75μl的ribolock、1μl的100mm dtt、1μl的5
×
第一链缓冲液和1.25μl的superscript iii(life technologies)用于随机rt。cdna延伸在4℃下进行2分钟，在15℃下进行3分钟，在25℃下进行10分钟，在42℃下进行45分钟，并且在50℃下进行25分钟。向rt产物中加入5μl的rnase i(thermo fisher scientific)、3μl的10
×
tnf缓冲液和2μl的h2o，并在37℃下孵育30分钟。在cdna延伸后，样品应保持在37℃下以避免变性条件。
72.myone c1磁珠(invitrogen)(20μl/样品)通过用1ml的磁珠结合缓冲液(100mm tris-hcl ph7.0、1m nacl、10mm edta)洗涤三次并重悬于供应有1μl ribolock的10μl的磁珠结合缓冲液中来制备。rnase i消化产物与预洗涤的珠混合，并在室温下旋转孵育45分钟。在用500μl的洗涤缓冲液(100mm tris ph7.0、4m nacl、10mm edta和0.2％tween-20)洗涤五次并用500μl的1
×
pbs洗涤两次后，将与cdna样品结合的磁珠用40μl的h2o重悬。通过添加5μl的1m naoh洗脱cdna，并在混匀仪中以1000r.p.m在70℃下孵育15分钟以完全消化rna。将样品立即置于磁体上，并将45μl的cdna洗脱液移至新管中并加入5μl的1m hcl。然后用zymo dna clean&concentrator-5柱纯化洗脱液。在rnase i消化后，将dmso组直接孵育并用naoh纯化。将纯化的样品与1μl(1u)的fastap(thermo fisher scientific)、3μl的10
×
circligase ii(epicentre)和1.5μl的mncl2混合，在37℃下孵育10分钟并在95℃下孵育2分钟以进行末端修复。加入由12μl的50％peg-4000(sigma)、1.5μl的circligase ii(epicentre)和1μl的10μm 3'接头(参见表1)组成的连接混合物并通过强涡旋混合。将反应物在60℃下孵育2小时并冷却至4℃。
73.表1：3'接头体系
[0074][0075]
其中，所述seq id no.3的3’末端的c优选以dd修饰；所述seq id no.4中3’末端的tcac可选的进行硫代修饰；seq id no.6中gagat和gtgac之间可选插入索引序列。
[0076]
3、3'接头连接和第二链合成
[0077]
通过用500μl的结合缓冲液(10mm tris-hcl ph8.0、1m nacl、1mm edta、0.05％
tween-20、0.5％sds)洗涤两次并重悬于250μl的结合缓冲液中来制备myone c1磁珠(invitrogen)(20μl/样品)。将连接产物在95℃下加热2分钟，立即转移到冰上至少1分钟，并与预洗涤的磁珠在室温下旋转孵育20分钟。然后将珠用200μl的洗涤缓冲液a(10mm tris-hcl ph8.0、100mm nacl、1mm edta、0.05％tween-20、0.5％sds)洗涤一次，并用200μl的洗涤缓冲液b(10mm tris-hcl ph8.0、100mm nacl、1mm edta、0.05％tween)洗涤一次。
[0078]
用47μl的主混合物重悬磁珠，所述主混合物由40.5μl的h2o、5μl的10x等温扩增缓冲液(neb)、0.5μl的25mm dntp(thermo fisher scientific)和1μl的100μm延伸引物组成。将混合物在65℃下在混匀仪中以1000r.p.m孵育2分钟，在冰上冷却1分钟并转移到预冷却的15℃混匀仪中，然后加入3μl的bst 2.0dna聚合酶(neb)。将延伸反应物从15℃孵育至37℃(1℃/分钟)，并在混匀仪中以1500r.p.m在37℃下保持5分钟(每分钟混合15秒)。将磁珠用200μl的洗涤缓冲液a洗涤一次，并用50μl的严格洗涤缓冲液(0.1x ssc缓冲液，0.1％sds)在55℃下在混匀仪中以1500r.p.m洗涤一次(每分钟混合15秒)，并用200μl的洗涤缓冲液b洗涤一次。将磁珠用99μl的主混合物重悬，所述主混合物由86.1μl的h2o、10μl的10x tango缓冲液(thermo fisher scientific)、2.5μl的1％tween-20和0.4μl的25mm dntp和1μl的t4 dna聚合酶(thermo fisher scientific)组成。将反应物在混匀仪中以1500r.p.m在25℃下孵育15分钟(每分钟混合15秒)。如上所述，将珠洗涤三次。
[0079]
4、5'接头连接和扩增
[0080]
用98μl的主混合物重悬磁珠，所述主混合物由73.5μl的h2o、10μl的10x t4 dna连接酶缓冲液(thermo fisher scientific)、10μl的50％peg-4000(thermo fisher scientific)、2.5μl的1％tween-20和2μl的100μm双链接头(dsa)组成(参见表1)。通过在95℃下加热两个互补寡核苷酸10秒并缓慢冷却至14℃(0.1℃/秒)使dsa退火。在添加2μl(10u)的t4 dna连接酶(thermo fisher scientific)后，将连接反应物在25℃下在混匀仪中以1500r.p.m孵育1小时(每分钟混合15秒)。如上所述将珠洗涤三次，然后重悬于25μl的洗脱缓冲液(10mm tris-hcl ph8.0，0.05％tween-20)中，并在95℃下孵育10分钟。收集上清液用于扩增。
[0081]
在40μl的qpcr反应物(12μl的cdna，20μl的2x phusion hf主混合物，0.75μl的10μm p7索引引物(参见表1)，0.75μl的10μm p5引物(参见表1)，0.4μl的25x sybrgold)中扩增样品。qpcr仪编程如下：98℃ 1分钟，98℃ 15秒，65℃ 30秒，72℃ 45秒。qpcr扩增后，通过6％非变性page凝胶对样品进行大小选择(》150bp)。在用qubit(invitrogen)定量后，在hiseq x ten(illumina)上运行深度测序。
[0082]
三、用于smartshape分数计算的计算管道。
[0083]
由于插入序列大部分短于100nt，我们仅使用读段配对物1进行后续处理。使用icshape-pipe处理smartshape测序数据。处理步骤如下：1)用cutadapt去除3'接头；2)去除重复的读段；3)使用trimmomatic去除前10nt；4)使用bowtie2将干净的读段映射到人rrna；5)使用star将未映射的读段比对到人(hg38)或小鼠(mm10)基因组；6)使用icshape-pipe sam2tab将sam文件转换成.tab文件；7)使用icshape-pipe calcshapenocont计算smartshape分数，其中参数为：-n nai_rep1.tab,nai_rep2.tab；-size chrnamelength.txt；-out reactivity.gtab；-ijf sjdblist.fromgtf.out.tab。sjdblist.fromgtf.out.tab文件和chrnamelength.txt文件由star在基因组索引生成期间
产生。
[0084]
icshape-pipe基本上基于滑动窗口方案计算全基因组的smartshape值，默认窗口大小为200nt，步长为5nt，在定义窗口时跳过非编码区域，直接串联外显子。每个核苷酸计算40次，并且在计算期间仅考虑附近的核苷酸，以避免由每个转录本的不同区段的不均匀覆盖引起的偏倚。当读段的5'被比对到3'相邻位点(+1位置)时，每个位点的逆转录终止信号加一。在每个窗口内对逆转录终止信号进行归一化，并进行90％缩尾处理以获得范围为0至1的最终分数。每个碱基的最终smartshape值是包含碱基的所有窗口的平均分数。如果覆盖率低于100，则smartshape值被定义为空(null)，这意味着未能在这些位点检测到结构。
[0085]
四、rna结构分析
[0086]
用python软件包sklearn生成接受者操作特征(roc)曲线。总之，考虑到二级结构和shape分数列表(0-1)，单链碱基被认为是阳性样品，双链碱基被认为是阴性样品。如果采用shape分数的截止值将所有碱基分成阳性样品和阴性样品，则可以计算假阳性率(ftr)和真阳性率(tpr)。因此，可以通过将该截止值从0逐渐调整到1来计算roc曲线。auc是roc曲线下面积。
[0087]
rna结构建模：rna二级结构用rnastructure软件包中的fold程序建模。smartshape分数可以用作约束，斜率和截距参数设置为默认。
[0088]
实施例2：利用rnase i消化去除m1a修饰引起的背景信号
[0089]
将用nai-n3修饰的hek293t的生物素化总rna与3.5μl的特异性rt引物和3μl的5
×
第一链缓冲液混合，加热至65℃持续5分钟，并在冰上孵育2分钟。将退火的样品与0.75μl的ribolock、1μl的100mm dtt、1μl的5
×
第一链缓冲液和1.25μl的superscript iii(life technologies)混合，并在55℃下孵育30分钟。将rt产物分成5部分，其中一组省略rnase i消化和磁珠富集，并且一组直接进行磁珠富集。将其余组分别与10μl、5μl或2.5μl的rnase i在30μl的反应体系中孵育。通过myone c1磁珠富集样品，并与naoh一起孵育用于如上所述的洗脱。最后，用zymo dna clean&concentrator-5柱纯化所有样品，并通过7m尿素page分离。
[0090]
icshape和smartshape中的nai-n3修饰单链核苷酸并引起逆转录(rt)停止。然而，逆转录酶也在一些内源修饰如m1a、局部结构如g-四链体位点停止，或只是偶然地停止在未修饰的位点处。这些背景逆转录终止信号将在结构分数的计算中引起假阳性信号。因此，在以前的rna结构检测方法中，增加了dmso对照组来去掉背景信号。而在smartshape中，我们在逆转录后引入了rnase i消化步骤，来去掉非修饰位点的终止信号。如图3a所示，在逆转录过程中，一条rna上可能结合多个逆转录引物，转录出多个cdna分子。只要rna上存在一个修饰位点，就能将其上所有的cdna分子富集出来，其中可能包含非修饰位点引起的假信号。rnase i可以特异性的切割单链rna，但无法切割rna-cdna杂合链。因此，rnase i消化可以将不同的cdna分子切割成单独的片段，从而避免了背景信号的富集。理论上，在smartshape文库中捕获的所有rt信号对应于检测剂的真实修饰，因此可以省略dmso组以进一步节省起始材料、劳动力和测序成本。
[0091]
为了验证rnase i消化步骤如预期的那样起作用以去除背景逆转录终止信号，我们在人核糖体rna 28s内的已知m1a修饰位点上游设计了rt引物(图3b)。我们用nai-n3处理hek293t细胞，分离rna，并进行click-it生物素化，然后进行逆转录(具体参见实施例1)。对
于没有经过rnase i处理的样品，我们在链霉亲和素磁珠富集后，除了全长cdna之外，还观察到对应于m1a位点的强背景逆转录终止信号，rnase i消化后未能检测到该条带，说明，以nai-n3修饰的hek293t总rna为模版进行逆转录，对逆转录产物同时进行rnase i消化和磁珠富集时，能够有效去除m1a修饰引起的背景逆转录信号(参见图3c)。重要的是，rnase i处理后进行链霉亲和素珠富集消除了这种m1a位点相关的rt产物。我们用含有m1a修饰的合成rna寡核苷酸重复该分析，并观察到由m1a位点产生的rt产物也被rnase i消化和磁珠富集操作消除了(参见图4a-b)。
[0092]
为了进一步评估smartshape测序数据中背景信号的去除，我们从用nai-n3和dmso处理过的hek293t细胞中构建文库(参见图4c)。为了鉴定背景信号，我们在dmso文库的构建期间省略了rna-cdna杂合链霉亲和素珠富集步骤。我们的结果揭示了dmso组中能观察到对应于已知内源性m1a修饰位点的背景信号(参见图3d)。重要的是，这些强背景逆转录终止信号在nai-n3文库中大大减少。请注意，对于不诱导rt终止的所有其他内源修饰位点(例如am和um)，我们观察到nai-n3和dmso文库之间逆转录终止信号的平均数的差异很少，这表明rnase i消化步骤特异性地去除了背景信号(图4d)。
[0093]
实施例3：不同起始量rna的smartshape的性能
[0094]
为了评估具有不同起始量rna的smartshape的性能，我们利用1ng、5ng、25ng和125ng的rna(在rrna去除后)作为起始量构建smartshape文库，来检测hek293t细胞中的全转录组rna二级结构。所有smartshape文库在具有不同起始量的文库之间(参见图5a中的示例和图6a中的总体统计)以及在相同起始量的文库之间(参见图6b)都显示出良好的可重复性。如果超过80％的核苷酸获得有效的smartshape分数，则转录本被定义为具有“高覆盖率”。由5ng、25ng和125ng的rna作为起始量生成的文库成功地在250m测序深度下检测到超过12,000个高覆盖率转录本的二级结构，其中超过75％的转录物是mrna和lncrna，5ng、25ng和125ng smartshape文库检测到的转录本数目远高于icshape，1ng smartshape文库检测到的转录本数目与icshape相当(参见图5b，以最深测序深度为标准，从右至左依次为1ng、icshape、5ng、25ng和125ng)。因此，在这些文库中，smartshape在相同的测序深度下显示出比icshape更高的覆盖度(参见图5b)。
[0095]
为了评估不同测序深度下每个文库的复杂性，我们从每个文库的总原始测序数据中随机取样相同数量的读段(表2)，并相应地计算smartshape分数。如图5b所示，5ng、25ng和125ng文库在测序深度超过250兆时能够测得的高覆盖度转录本数目仍然快速上升，说明文库都具有很高的复杂度且尚未饱和，可以通过提高测序深度来获得更多转录本的信息。此外，这三个文库在不同测序深度下，平均逆转录终止信号的分布非常接近，表明5ng的rna起始量足以构建高度复杂的smartshape文库(参见图5b，图6c，其中，图6c中从左下往上的曲线依次代表50m至250m)。最后，尽管我们确实察觉到1ng起始量rna文库的复杂性有一定程度的降低，但我们仍然在250m测序深度下获得了超过9,000个高覆盖率转录本，这一水平与相同测序深度下的icshape相当(其需要约500ng的起始量rna)。
[0096]
表2不同文库测序深度及不同处理步骤对应读段数目
[0097][0098][0099]
我们进一步比较了每个文库中可用测序读段的比例。icshape和smartshape都使用与3'接头相邻的随机序列分子标签来标记pcr重复。pcr重复的读段和因太短而不能比对至基因组或比对到rrnas的读段对于计算rna结构分数都是无用的，需要丢弃。剩余的读段(比对到基因组的读段)被定义为可用读段。我们观察到超过60％的总测序读段在5ng、25ng和125ng文库中是可用的，而相比之下，从500ng的起始量rna生成的icshape文库中仅有约40％的可用读段，可见，5ng、25ng和125ng smartshape文库能够比对的基因组的读段数目远高于icshape文库(参见图5c)。然而，在1ng文库中仅约20％的读段是可用的，考虑到测序成本，我们建议使用超过1ng的rna作为起始量进行smartshape建库(参见图5c)。
[0100]
为了评估smartshape的准确性，我们使用计算出的smartshape值绘制18s和28s rrna中可被修饰的碱基的roc曲线。不同起始量smartshape文库18s的auc超过0.8，28s的auc超过0.7，这表明smartshape数据与已知的结构模型具有良好的一致性，且smartshape文库的准确度明显高于icshape(参见图5d)。我们还用人xbp1转录本中的已知结构元件评估了smartshape值。实际上，我们观察到smartshape值与已知的结构模型具有良好的一致性，smartshape文库的曲线下面积明显高于icshape文库(参见图5e)。
[0101]
我们还检查了smartshape库的其他质量控制参数。与先前的发现类似，smartshape数据揭示了翻译起始和终止位点处的结构特点，以及cds区段中的3-核苷酸周期性(参见图7a)。由于与cg碱基对相比，au的氢键通常较弱，因此a和u核苷酸处的smartshape值高于c和g核苷酸处的smartshape值(参见图7b)。与smartshape数据中含有相同“ggacu”基序的背景区段相比，m6a甲基化的区段显示出更高的smartshape值，这与m6a区段倾向于单链的结论一致(参见图7c)。gini指数用于定量转录本中rna结构的致密程度，gini指数越高表明双链rna结构越多。mrna和lncrna的gini指数值低于假基因、mirna和snorna的gini指数值，这与先前的发现一致(参见图7d)。
[0102]
总之，smartshape可以准确且可靠地检测不同起始量样品中的rna结构，同时仅需要其他最先进的体内rna结构检测方法所需的起始量rna的一小部分，当使用少量，例如1ng，rna作为起始量，smartshape仍然可准确检测rna结构。因此，对于大量样品材料的获取极具挑战的许多生物医学应用来说，smartshape应该非常合适。
[0103]
实施例4：用于smartshape分数计算的计算管道。
[0104]
我们开发了一种新的分析管道，用于仅基于nai-n3文库计算rna结构分数(参见实
施例1)。简言之，通过在所有外显子上以滑动窗口方式对rt终止信号进行归一化和缩尾处理来计算smartshape值，并且将覆盖率低于100的碱基的smartshape值定义为空(null)(默认窗口大小＝20nt，步长＝5nt)。我们用人核糖体rna 18s已知结构模型评估了新管道的性能(参见实施例1)。通过绘制接受者操作特征(roc)曲线，我们观察到与发表的icshape数据相比，用新管道计算出的smartshape分数的表现更好，smartshape值的曲线下面积(auc)明显高于icshape值(参见图3e-f)。这些结果进一步表明，rnase i消化和链霉亲和素珠富集步骤有效去除了背景信号，使得不再需要dmso文库作为对照。
[0105]
实施例5：smartshape测量小鼠巨噬细胞中全转录组水平的rna结构
[0106]
使鼠柠檬酸杆菌在lb肉汤中37℃下振荡生长过夜。c57bl/6j小鼠(6-8周)通过灌胃感染总体积为200μl的2
×
109cfus鼠柠檬酸杆菌，并在感染后第5天处死。取出肠道组织并置于冰冷的不含钙和镁的汉克平衡盐溶液(hbss)中。将肠道纵向切开并切成1.5厘米的片，并在含有10mm hepes、10mm edta(promega)和1mm二硫苏糖醇(dtt，fermentas)的hbss中在37℃孵育两次，持续20分钟，以除去上皮细胞和粘液。然后，在用含有10mm hepes的hbss洗涤后，将组织在含有5％热灭活的胎牛血清(fbs)、1mg/ml胶原酶ⅳ(sigma)、1mg/ml分散酶(roche)和100μg/ml dnase i(sigma)的rpmi1640(含钙和镁)中在37℃下缓慢旋转消化75分钟。通过剧烈振荡将消化的组织均质化，通过70μm细胞过滤器并重悬于40％percoll(ge health care)溶液中，然后在室温下以2,500rpm梯度密度离心20分钟。并且使用ack溶解缓冲液溶解红细胞。染色后，在facsaria4激光(bd)上分选ly6c
+
和ly6c-结肠巨噬细胞。
[0107]
先天免疫受到精确调节以有效消除病原体，同时避免由过度免疫应答引起的组织损伤。这些免疫应答的介质通常显示出瞬时表达以诱导并随后消除炎症。转录后调节对于快速抑制关键炎性介质的蛋白质表达至关重要，其中rna结构在rna降解和翻译的调节中起重要作用。例如，gait元件(哺乳动物细胞中唯一的核糖开关)通过在转换为发夹构象时募集gait复合物来阻断巨噬细胞中vegfa基因的翻译。
[0108]
为了鉴定免疫细胞中新的转录后调控rna结构元件，我们使用smartshape来检测从感染了鼠柠檬酸杆菌的小鼠分离出的肠道巨噬细胞中的rna二级结构全转录组(参见图8a和图9a)，通过向小鼠中感染小鼠柠檬酸杆菌，构建小鼠肠道炎症模型，并在五天后从肠道中分选出ly6c
lo
组织常驻巨噬细胞和ly6c
hi
促炎巨噬细胞，最后利用smartshape分别测量了两种肠道巨噬细胞中的rna二级结构。每只小鼠只有5
×
104个肠道巨噬细胞，已有的rna结构检测方法无法进行检测。值得注意的是，据我们所知，这是哺乳动物免疫细胞的第一个全局rna结构数据。
[0109]
肠道巨噬细胞对于维持肠道中免疫应答和抗原耐受性之间的平衡是必需的。具体地，从血液募集的单核细胞分化成ly6c
lo
组织常驻巨噬细胞，其通过产生抗炎细胞因子如白细胞介素(il)-10来维持肠稳态。然而，在肠道炎症期间，循环单核细胞分化成ly6c
hi
促炎巨噬细胞，其通过产生促炎细胞因子如il6、il1b和il12触发炎症。为了探索rna结构在组织常驻和促炎性巨噬细胞中的潜在差异，我们使用约100ng总rna对ly6c
lo
和ly6c
hi
巨噬细胞进行smartshape文库构建。在ly6c
lo
和ly6c
hi
巨噬细胞的smartshape数据中，我们分别获得了超过3,000个和超过2,000个高覆盖率转录本的结构信息(参见图8b)。xbp1转录本已知结构元件和srp rna的smartshape值显示出与已知结构模型良好的一致性，并且与icshape分数
相比，明显具有高得多的auc(参见图8c和图10a)。在一组结构已知的60个rna中，两种巨噬细胞smartshape值的auc均值远高于发表的小鼠胚胎干细胞icshape值的auc，这表明smartshape数据质量很高(参见图10b)。
[0110]
可见本发明的rna结构检测方法的结果可用于评估细胞的功能状态，例如，免疫应激反应。类似的，rna结构检测方法的结果可以用于评估细胞的其它功能状态，例如研究rna对早期发育的影响，癌症的发生和发展等。
[0111]
以上详细描述了本发明的优选实施方式，但是，本发明并不限于上述实施方式中的具体细节，在本发明的技术构思范围内，可以对本发明的技术方案进行多种简单变型，这些简单变型均属于本发明的保护范围。
[0112]
另外需要说明的是，在上述具体实施方式中所描述的各个具体技术特征，在不矛盾的情况下，可以通过任何合适的方式进行组合，为了避免不必要的重复，本发明对各种可能的组合方式不再另行说明。
[0113]
此外，本发明的各种不同的实施方式之间也可以进行任意组合，只要其不违背本发明的思想，其同样应当视为本发明所公开的内容。