Name: 使用硫酸二甲酯探测RNA结构 体 外 和细胞测序突变分析 (视频) | JoVE Journal
Uploaded: 2022-12-09T14:50:00.000+00:00
Duration: 10 min 34 s
Description: 学习使用硫酸二甲酯进行RNA结构探测的DMS-MaP方法，通过体外和细胞测序突变分析揭示RNA动态变化。该技术为生物学研究提供高通量结构预测工具。

ai_agent_IP

ai_agent

maintenance

instantTrialDurationMinutes

isBiopharma

enableNewMyJoveSectionDesign

educationFeatures

seoMetaApi

ai_copilot_institution_blacklist

ai_chat_institution_blacklist

freeAccessWindow

ssoRoleCapture

scheduledMaintenanceBanner

joveLabsFeatures

disable_autoplay_for_bots

extended_trial_duration

5 minutes

注意：有关本方案中使用的所有材料、软件、试剂、仪器和细胞的详细信息，请参见材料表。
1. 基因特异性 体外 DMS-MaP
<ol>
	<li>RNA 体外 转录
	<ol>
		<li>获取目标RNA的双链（ds）DNA序列（例如，以DNA片段、质粒或通过已有/基因组DNA进行PCR获得）。若该DNA序列含有聚合酶启动子，则直接进入步骤3。</li>
		<li>通过重叠 PCR 将 RNA 聚合酶启动子连接至目标 DNA 片段的上游（T7 聚合酶的上游引物：5' TAATACGACTCACTATAGG + 目标序列起始碱基 3'）。</li>
		<li>体外 将DNA片段转录为RNA。始终将RNA置于冰上。</li>
		<li>使用DNase消化DNA。</li>
		<li>使用基于柱的方法（步骤2.4）或乙醇沉淀法（步骤2.5）分离RNA。用适当体积洗脱，预期产量约为50 µg。</li>
		<li>通过琼脂糖凝胶电泳检测RNA完整性；上样前将RNA在70 °C下变性2-3分钟。 
		注意：缓冲液和琼脂糖可能含有降解RNA的RNase，可能会污染RNA样品。本实验室此前曾使用预制琼脂糖凝胶，但结果（尤其是RNA）有时不明确。最佳结果来自琼脂糖凝胶或PAGE凝胶。</li>
		<li>将RNA直接于-80 °C保存数月，除非解冻后可见降解现象。</li>
	</ol>
	</li>
	<li>体外 DMS修饰（在105 mM DMS条件下）
	<ol>
		<li>配制足量的复性缓冲液（0.4 M 焦碳酸钠，pH 7.2，含 6 mM MgCl₂）2). 
		注意：每个反应（终体积为 100 µL）中加入 89 µL 复性缓冲液。</li>
		<li>每个反应中，将89 µL复性缓冲液转移至指定的1.5 mL离心管中，并在化学通风橱下方的恒温振荡器中于37 °C预热。 
		注意：DMS 具有高毒性，必须始终在化学通风橱中操作，直至用还原剂淬灭。</li>
		<li>用 10 µL 无核酸酶水（NF H₂O）洗脱 1–10 pmol RNA2O）；转移至PCR管。</li>
		<li>在热循环仪中于 95 °C 孵育 1 分钟以变性 RNA。</li>
		<li>置于冰块上 立即 以避免错误折叠.</li>
		<li>将RNA样品加入含有复性缓冲液的指定离心管中，37 °C下充分混匀，孵育10–20分钟以使RNA复性。 
		注意：大多数RNA的折叠时间为毫秒到秒级，尽管存在例外情况16.</li>
		<li>向RNA样品中加入1 µL 100%（10.5 M）DMS，并在800–1,400转/分钟（rpm）振荡条件下孵育5分钟。 
		注意：此步骤中振荡（或其他混合方式）至关重要，因为DMS具有疏水性，可能无法在复性缓冲液中完全溶解。反应时间的偏差可能影响DMS反应活性的可重复性。为尽量减少移液误差，可在将DMS加入样品前先将其溶于100%乙醇中，前提是保持DMS的终浓度为1%（105 mM）。对于未处理对照组，可用二甲基亚砜（DMSO）或水替代DMS。</li>
		<li>反应5分钟后，加入60 µL 100% β-巯基乙醇（BME）终止反应，充分混匀，并立即将RNA置于冰上。 
		注意：在使用BME淬灭反应后，可将RNA安全地移出通风橱进行后续处理。然而，由于BME具有强烈气味和刺激性，仍应避免其直接暴露于环境中。</li>
		<li>使用乙酸钠-乙醇沉淀法（见步骤2.5）或柱纯化法（见步骤2.6）纯化RNA，并用10 µL水洗脱。</li>
		<li>使用分光光度计定量RNA。</li>
		<li>将修饰后的RNA直接存放在−80 °C保存。 
		注意：应避免长期保存，因为RNA在DMS处理后稳定性降低。</li>
	</ol>
	</li>
	<li>基因特异性修饰RNA的RT-PCR 
	注意：参见 图1 用于DMS处理片段的RT-PCR体系构建。
	<ol>
		<li>用 10 µL 无核酸酶（NF）水洗脱 100 ng 修饰后的 RNA2O. 转移至PCR管中。</li>
		<li>向离心管中加入 4 µL 5x 第一链缓冲液（FSB）、1 µL dNTP 混合液（各 10 mM）、1 µL 0.1 M 二硫苏糖醇（DTT）（避免冻融循环）、1 µL RNase 抑制剂、1 µL 10 µM 反向引物（单个引物或引物混合物）和 1 µL TGIRT III。 
		注意：对于引物混合物，不要将每种10 µM的引物各1 µL直接加入反转录反应；而应先将引物混合，再从混合物中取1 µL（总引物浓度为10 µM）加入。</li>
		<li>在57 °C孵育30分钟至1.5小时（通常30分钟足以生成500 nt的产物），使用PCR仪进行。</li>
		<li>加入1 µL 4 M NaOH，用移液器吹打混匀，于95 °C孵育3分钟以降解RNA。 
		注意：此步骤至关重要，因为通过降解RNA可将TGIRT从cDNA上释放。若省略该步骤，可能会影响下游的PCR反应。</li>
		<li>使用基于柱的方法进行纯化（参见步骤2.6），充分去除引物，并用10 µL NF H洗脱2O.</li>
		<li>使用专为平衡产量与保真度设计的PCR试剂盒，以每25 µL反应体系加入1 µL逆转录产物的量对cDNA进行PCR扩增。 
		注意：引物的熔解温度应约为 60 °C。</li>
		<li>取 2 µL PCR 产物在琼脂糖凝胶或预制琼脂糖凝胶上电泳，以验证 PCR 扩增结果。</li>
		<li>理想情况下，PCR 后应仅出现一条条带。若如此，使用基于柱子的方法纯化反应产物。若存在其他条带，则利用剩余的 PCR 反应产物从凝胶中切下正确条带，并用足够小的洗脱体积（例如 10 µL）进行洗脱。</li>
		<li>使用分光光度计对提取的片段进行定量。</li>
		<li>使用适合目标测序平台的方法对dsDNA片段进行测序索引。</li>
	</ol>
	</li>
</ol>
<img alt="DNA测序方法示意图，显示引物、正向与反向链，每段300bp。" class="xfigimg" src="/files/ftp_upload/64820/64820fig01.jpg" data-altupdatedat="1762425657000" data-alt="图1"> 
图1：针对大尺寸DMS处理片段的RT-PCR实验设计。 在对修饰RNA进行反转录时，引物结合区域上的修饰不会被记录。因此，当片段长度超过400–500 bp时，需要设计在引物区域重叠的片段，如此处所示。片段长度取决于测序需求。使用双端150个循环测序时，片段长度不应超过300 bp。缩写：RT-PCR = 逆转录聚合酶链式反应；DMS = 二甲基硫酸酯。<a href="https://www.jove.com/files/ftp_upload/64820/64820fig01large.jpg" target="_blank">请点击此处查看此图的放大版本。</a>
2. 利用病毒感染细胞进行全基因组 DMS-MaP 分析
注意：在细胞中，DMS 处理也可与上述基因特异性扩增方法联合使用。全基因组文库需要极深的测序深度才能实现对单个基因的完全覆盖。然而，如果提取后病毒 RNA 在去核糖体 RNA 中占显著比例，则全基因组测序是合适的选择。此外，其他富集方法也可与全基因组文库构建方法结合使用。 
注意：视频方案中出于演示目的使用了未感染的细胞。
<ol>
	<li>DMS 处理
	<ol>
		<li>培养病毒感染的细胞，直至达到所需的感染阶段。</li>
		<li>将细胞容器转移至专用的通风橱中，该通风橱需符合相应生物安全等级对病毒操作的要求，并能有效处理 DMS 等试剂产生的化学烟雾。</li>
		<li>向培养基中加入 2.5% 体积的 DMS，并用封口膜密封容器（通常为 10 cm 培养皿）。 
		注意：使用 DMS 时容易出现修饰不足或过度修饰的情况。直接向细胞中添加 DMS 时，充分混匀至关重要。也可选择在 50 mL 锥形管中将新鲜培养基预热至 37 °C，然后直接加入 DMS 并剧烈震荡。弃去细胞上的旧培养基，缓慢加入含 DMS 的培养基。</li>
		<li>转移至 37 °C 培养箱中孵育 5 分钟。 
		注意：根据在培养箱外操作 DMS 所需时间的不同，5 分钟可能导致过度修饰。应尽量控制从添加 DMS 到开始孵育的时间 ≤1 分钟。若首次进行该实验，建议进行 DMS 梯度浓度测试，并调整孵育时间（3 至 10 分钟之间），以确定最佳修饰效率，并确保结果在一定浓度范围内具有可重复性。</li>
		<li>小心吸出含 DMS 的培养基（倒入适当的化学废弃物容器），并轻柔加入 10 mL 终止缓冲液（含 30% BME 的 PBS，例如 3 mL BME 与 7 mL PBS 混合）。 
		注意：DMS 和 BME 的加入可能导致贴壁性较弱的细胞从培养皿上脱落。若发生细胞脱落，可将其作为悬浮细胞处理——无需先移除含 DMS 的培养基，而是直接加入终止缓冲液，用细胞刮刀将含有 DMS 和 BME 的细胞刮下并转移至 50 mL 锥形管中。在 3,000 × g 条件下离心 3 分钟收集细胞；务必彻底清除残留的 DMS，因其可能在细胞下方形成大液滴而沉淀。若初始阶段无法完全去除 DMS 培养基，建议额外使用含 30% BME 的缓冲液清洗一次。</li>
		<li>用细胞刮刀刮取细胞，并转移至 15 mL 锥形管中。</li>
		<li>在 3,000 × g 条件下离心 3 分钟以收集细胞沉淀。</li>
		<li>弃去上清液，并用 10 mL PBS 洗涤两次。</li>
		<li>尽可能彻底地去除残余的 PBS。</li>
		<li>用适量 RNA 提取试剂溶解沉淀（例如，T75 培养瓶使用 3 mL，10 cm 培养皿使用 1 mL）。 
		注意：试剂用量不足可能影响 RNA 产量。</li>
	</ol>
	</li>
	<li>RNA 提取与核糖体 RNA (rRNA) 去除
	<ol>
		<li>向含有 1 mL RNA 提取试剂的匀浆细胞样品中加入 200 µL 氯仿，涡旋震荡 15–20 秒至溶液呈亮粉色，随后孵育最多 3 分钟，直至出现明显的相分离。 
		注意：粉色脂质相应位于下层。若未观察到此现象，则可能是涡旋时间不足。</li>
		<li>在 4 °C 条件下以最大转速（约 20,000 × g）离心 15 分钟。</li>
		<li>将上层水相转移至新离心管中。</li>
		<li>采用乙酸钠-乙醇沉淀法（见步骤 2.5）或柱纯化法（见步骤 2.6）对 RNA 进行纯化，并用足量无核酸酶水（NF H2O）洗脱。</li>
		<li>通过琼脂糖凝胶电泳检测 RNA 完整性，应可见对应于两个核糖体亚基的两条主带。</li>
		<li>采用首选方法去除 rRNA，并用适量（通常为 20–50 µL）无核酸酶水（NF H2O）洗脱。 
		​注意：后续实验建议使用约 500 ng 总 RNA，溶于 8 µL 体积中。非核糖体 RNA 通常仅占总 RNA 的 5%–10%。</li>
		<li>使用分光光度计进行定量。</li>
	</ol>
	</li>
	<li>文库构建
	<ol>
		<li>采用基因特异性 RT-PCR 或其他方法构建文库15。若使用随机六聚体引物，需增加一个低温退火步骤（Tm 为 37–42 °C），以促进六聚体与模板结合。 
		注意：也可使用标准文库构建试剂盒，只需将逆转录酶替换为 TGIRT，并将逆转录温度调整为 57 °C 即可。</li>
	</ol>
	</li>
	<li>使用 RNA Clean &amp; Concentrator 柱进行柱式 RNA 纯化 
	​注意：所有步骤均应在室温下进行。
	<ol>
		<li>向样品管中加入无核酸酶水（NF H2O），使总体积达到 50 µL。</li>
		<li>向样品中加入 100 µL 结合缓冲液和 150 µL 100% 乙醇。</li>
		<li>混匀后转移至离心柱中。</li>
		<li>在 10,000–16,000 × g 条件下离心 30 秒，弃去穿流液。</li>
		<li>加入 400 µL RNA 预洗缓冲液。</li>
		<li>在 10,000–16,000 × g 条件下离心 30 秒，弃去穿流液。</li>
		<li>加入 700 µL RNA 洗涤缓冲液。</li>
		<li>在 10,000–16,000 × g 条件下离心 30 秒，弃去穿流液。</li>
		<li>加入 400 µL RNA 洗涤缓冲液。</li>
		<li>在 10,000–16,000 × g 条件下离心 30 秒，弃去穿流液。</li>
		<li>（可选）将离心柱转移至新的收集管中，在 10,000–16,000 × g 条件下离心 2 分钟。</li>
		<li>将离心柱转移至洁净无 RNase 的离心管中，加入适量无核酸酶水（NF H2O）。</li>
		<li>在 10,000–16,000 × g 条件下离心 1 分钟。</li>
	</ol>
	</li>
	<li>酸性酚-氯仿法 RNA 提取
	<ol>
		<li>加入等体积的酸性酚:氯仿:异戊醇混合液。</li>
		<li>充分涡旋震荡，然后在 14,000 × g 条件下离心 5 分钟。</li>
		<li>若未出现相分离，加入 20 µL 2 M NaCl，再次离心。</li>
		<li>将水相转移至新的离心管中。</li>
		<li>加入 500 µL 异丙醇和 2 µL 共沉淀剂。</li>
		<li>混匀后在室温孵育 3 分钟，随后在 −80 °C 下过夜孵育。</li>
		<li>在 4 °C 条件下以最大转速（约 20,000 × g）离心 30 分钟以沉淀 RNA。</li>
		<li>用 200 µL 冰冷的 70% 乙醇洗涤沉淀。</li>
		<li>在最大转速（约 20,000 × g）下离心 5 分钟，弃去上清液。</li>
		<li>用适量无核酸酶水（NF H2O）重悬 RNA 沉淀。</li>
	</ol>
	</li>
	<li>使用 Oligo Clean and Concentrator 柱进行柱式 cDNA 纯化 
	注意：所有步骤均应在室温下进行。
	<ol>
		<li>向样品管中加入无核酸酶水（NF H2O），使总体积达到 50 µL。</li>
		<li>加入 100 µL 结合缓冲液和 400 µL 100% 乙醇。</li>
		<li>混匀后转移至离心柱中。</li>
		<li>在 10,000–16,000 × g 条件下离心 30 秒，弃去穿流液。</li>
		<li>加入 750 µL DNA 洗涤缓冲液。</li>
		<li>在 10,000–16,000 × g 条件下离心 30 秒，弃去穿流液。</li>
		<li>（可选）将离心柱转移至新的收集管中，在 10,000–16,000 × g 条件下离心 2 分钟。</li>
		<li>将离心柱转移至洁净无 RNase 的离心管中，加入适量无核酸酶水（NF H2O）。</li>
		<li>在 10,000–16,000 × g 条件下离心 1 分钟。</li>
	</ol>
	</li>
</ol>
3. 测序数据分析
注意：要利用 DMS-MaP 测序数据构建 RNA 二级结构模型，必须对生成的 .fastq 文件进行多个不同步骤的处理。这些步骤可使用
<ol>
	<li>使用 TrimGalore 或 Cutadapt 去除接头序列。</li>
	<li>使用 Bowtie2 将测序读段比对至参考序列（.fasta 格式）。</li>
	<li>利用专门的 RNA 结构分析软件（例如 DREEM14、RNA-Framework17 或类似工具）进行读段计数，并生成反应性图谱。</li>
	<li>（可选）使用 DREEM14、DRACO17、DANCE-MaP18 或类似工具对读段进行聚类，以发现 RNA 的替代构象。</li>
	<li>基于反应性图谱，使用 RNAStructure12、ViennaRNA 或类似软件预测最小自由能结构。</li>
	<li>使用 VARNA（https://varna.lri.fr/）或类似工具对 RNA11 结构进行可视化。 
	注意：为提高实用性，DREEM（www.rnadreem.org）和 RNA-Framework19 等软件已将上述步骤 1–5 大幅整合至其分析流程中，从而简化了整体分析过程。然而，任何结构预测均需谨慎处理（例如，应通过验证结构与实验数据的一致性20 来确保可靠性）。</li>
</ol>

<ol><li>Kim, S. H., et al. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/pubmed?term=((Three-dimensional+tertiary+structure+of+yeast+phenylalanine+transfer+RNA%5BTitle%5D)+AND+%22Science%22%5BJournal%5D)">Three-dimensional tertiary structure of yeast phenylalanine transfer RNA.</a> Science. 185 (4149), 435-440 (1974).</li>
<li>Robertus, J. D., et al. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/pubmed?term=((Structure+of+yeast+phenylalanine+tRNA+at+3+%C3%85+resolution%5BTitle%5D)+AND+%22Nature%22%5BJournal%5D)">Structure of yeast phenylalanine tRNA at 3 Å resolution.</a> Nature. 250 (467), 546-551 (1974).</li>
<li>Zaug, A. J., Cech, T. R. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/pubmed?term=((In+vitro+splicing+of+the+ribosomal+RNA+precursor+in+nuclei+of+Tetrahymena%5BTitle%5D)+AND+%22Cell%22%5BJournal%5D)">In vitro splicing of the ribosomal RNA precursor in nuclei of Tetrahymena.</a> Cell. 19 (2), 331-338 (1980).</li>
<li>Zhao, Y., et al. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/pubmed?term=((NONCODE+2016%3A+An+informative+and+valuable+data+source+of+long+non-coding+RNAs%5BTitle%5D)+AND+%22Nucleic+Acids+Research%22%5BJournal%5D)">NONCODE 2016: An informative and valuable data source of long non-coding RNAs.</a> Nucleic Acids Research. 44, D203-D208 (2016).</li>
<li>Vandivier, L. E., Anderson, S. J., Foley, S. W., Gregory, B. D. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/pubmed?term=((The+conservation+and+function+of+RNA+secondary+structure+in+plants%5BTitle%5D)+AND+%22Annual+Review+of+Plant+Biology%22%5BJournal%5D)">The conservation and function of RNA secondary structure in plants.</a> Annual Review of Plant Biology. 67, 463(2016).</li>
<li>Jumper, J., et al. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/pubmed?term=((Highly+accurate+protein+structure+prediction+with+AlphaFold%5BTitle%5D)+AND+%22Nature%22%5BJournal%5D)">Highly accurate protein structure prediction with AlphaFold.</a> Nature. 596 (7873), 583-589 (2021).</li>
<li>Das, R. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/pubmed?term=((RNA+structure%3A+A+renaissance+begins%5BTitle%5D)+AND+%22Nature+Methods%22%5BJournal%5D)">RNA structure: A renaissance begins.</a> Nature Methods. 18 (5), 439-439 (2021).</li>
<li>Smola, M. J., Rice, G. M., Busan, S., Siegfried, N. A., Weeks, K. M. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/pubmed?term=((Selective+2%E2%80%B2-hydroxyl+acylation+analyzed+by+primer+extension+and+mutational+profiling+%28SHAPE-MaP%29+for+direct%2C+versatile+and+accurate+RNA+structure+analysis%5BTitle%5D)+AND+%22Nature+Protocols%22%5BJournal%5D)">Selective 2′-hydroxyl acylation analyzed by primer extension and mutational profiling (SHAPE-MaP) for direct, versatile and accurate RNA structure analysis.</a> Nature Protocols. 10 (11), 1643-1669 (2015).</li>
<li>Mathews, D. H., et al. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/pubmed?term=((Incorporating+chemical+modification+constraints+into+a+dynamic+programming+algorithm+for+prediction+of+RNA+secondary+structure%5BTitle%5D)+AND+%22Proceedings+of+the+National+Academy+of+Sciences+of+the+United+States+of+America%22%5BJournal%5D)">Incorporating chemical modification constraints into a dynamic programming algorithm for prediction of RNA secondary structure.</a> Proceedings of the National Academy of Sciences of the United States of America. 101 (19), 7287-7292 (2004).</li>
<li>Zuker, M., Stiegler, P. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/pubmed?term=((Optimal+computer+folding+of+large+RNA+sequences+using+thermodynamics+and+auxiliary+information%5BTitle%5D)+AND+%22Nucleic+Acids+Research%22%5BJournal%5D)">Optimal computer folding of large RNA sequences using thermodynamics and auxiliary information.</a> Nucleic Acids Research. 9 (1), 133-148 (1981).</li>
<li>Lorenz, R., et al. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/pubmed?term=((ViennaRNA+Package+2.0%5BTitle%5D)+AND+%22Algorithms+for+Molecular+Biology%22%5BJournal%5D)">ViennaRNA Package 2.0.</a> Algorithms for Molecular Biology. 6, (2011).</li>
<li>Reuter, J. S., Mathews, D. H. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/pubmed?term=((RNAstructure%3A+Software+for+RNA+secondary+structure+prediction+and+analysis%5BTitle%5D)+AND+%22BMC+Bioinformatics%22%5BJournal%5D)">RNAstructure: Software for RNA secondary structure prediction and analysis.</a> BMC Bioinformatics. 11, (2010).</li>
<li>Wells, S. E., Hughes, J. M. X., Igel, A. H., Ares, M. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/pubmed?term=((Use+of+dimethyl+sulfate+to+probe+RNA+structure+in+vivo%5BTitle%5D)+AND+%22Methods+in+Enzymology%22%5BJournal%5D)">Use of dimethyl sulfate to probe RNA structure in vivo.</a> Methods in Enzymology. , 479-493 (2000).</li>
<li>Tomezsko, P. J., et al. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/pubmed?term=((Determination+of+RNA+structural+diversity+and+its+role+in+HIV-1+RNA+splicing%5BTitle%5D)+AND+%22Nature%22%5BJournal%5D)">Determination of RNA structural diversity and its role in HIV-1 RNA splicing.</a> Nature. 582 (7812), (2020).</li>
<li>Zubradt, M., et al. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/pubmed?term=((DMS-MaPseq+for+genome-wide+or+targeted+RNA+structure+probing+in+vivo%5BTitle%5D)+AND+%22Nature+Methods%22%5BJournal%5D)">DMS-MaPseq for genome-wide or targeted RNA structure probing in vivo.</a> Nature Methods. 14 (1), (2017).</li>
<li>Woodson, S. A. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/pubmed?term=((Compact+intermediates+in+RNA+folding%5BTitle%5D)+AND+%22Annual+Reviews+in+Biophysics%22%5BJournal%5D)">Compact intermediates in RNA folding.</a> Annual Reviews in Biophysics. 39, (2010).</li>
<li>Morandi, E., et al. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/pubmed?term=((Genome-scale+deconvolution+of+RNA+structure+ensembles%5BTitle%5D)+AND+%22Nature+Methods%22%5BJournal%5D)">Genome-scale deconvolution of RNA structure ensembles.</a> Nature Methods. 18 (3), 249-252 (2021).</li>
<li>Olson, S. W., et al. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/pubmed?term=((Discovery+of+a+large-scale%2C+cell-state-responsive+allosteric+switch+in+the+7SK+RNA+using+DANCE-MaP%5BTitle%5D)+AND+%22Molecular+Cell%22%5BJournal%5D)">Discovery of a large-scale, cell-state-responsive allosteric switch in the 7SK RNA using DANCE-MaP.</a> Molecular Cell. 82 (9), 1708-1723 (2022).</li>
<li>Incarnato, D., Morandi, E., Simon, L. M., Oliviero, S. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/pubmed?term=((RNA+Framework%3A+An+all-in-one+toolkit+for+the+analysis+of+RNA+structures+and+post-transcriptional+modifications%5BTitle%5D)+AND+%22Nucleic+Acids+Research%22%5BJournal%5D)">RNA Framework: An all-in-one toolkit for the analysis of RNA structures and post-transcriptional modifications.</a> Nucleic Acids Research. 46 (16), (2018).</li>
<li>Lan, T. C. T., et al. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/pubmed?term=((Secondary+structural+ensembles+of+the+SARS-CoV-2+RNA+genome+in+infected+cells%5BTitle%5D)+AND+%22Nature+Communications%22%5BJournal%5D)">Secondary structural ensembles of the SARS-CoV-2 RNA genome in infected cells.</a> Nature Communications. 13 (1), 1128(2022).</li>
<li>Homan, P. J., et al. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/pubmed?term=((Single-molecule+correlated+chemical+probing+of+RNA%5BTitle%5D)+AND+%22Proceedings+of+the+National+Academy+of+Sciences+of+the+United+States+of+America%22%5BJournal%5D)">Single-molecule correlated chemical probing of RNA.</a> Proceedings of the National Academy of Sciences of the United States of America. 111 (38), 13858-13863 (2014).</li>
<li>Yang, S. L., et al. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/pubmed?term=((Comprehensive+mapping+of+SARS-CoV-2+interactions+in+vivo+reveals+functional+virus-host+interactions%5BTitle%5D)+AND+%22Nature+Communications%22%5BJournal%5D)">Comprehensive mapping of SARS-CoV-2 interactions in vivo reveals functional virus-host interactions.</a> Nature Communications. 12 (1), 5113(2021).</li>
<li>Manfredonia, I., et al. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/pubmed?term=((Genome-wide+mapping+of+SARS-CoV-2+RNA+structures+identifies+therapeutically-relevant+elements%5BTitle%5D)+AND+%22Nucleic+Acids+Research%22%5BJournal%5D)">Genome-wide mapping of SARS-CoV-2 RNA structures identifies therapeutically-relevant elements.</a> Nucleic Acids Research. 48 (22), 12436-12452 (2020).</li>
<li>Fischer, N. M., Polěto, M. D., Steuer, J., vander Spoel, D. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/pubmed?term=((Influence+of+Na%2B+and+Mg2%2B+ions+on+RNA+structures+studied+with+molecular+dynamics+simulations%5BTitle%5D)+AND+%22Nucleic+Acids+Research%22%5BJournal%5D)">Influence of Na+ and Mg2+ ions on RNA structures studied with molecular dynamics simulations.</a> Nucleic Acids Research. 46 (10), 4872-4882 (2018).</li>
</ol>

作者声明无任何利益冲突。

本方案描述了如何检测RNA 体外 并利用DMS突变谱分析实验在细胞中进行检测。此外，本文还提供了制备Illumina测序文库以生成基因特异性数据及分析所得.fastq文件的操作指南。同时，也可采用全基因组文库策略。然而，基因特异性的RT-PCR能够产生质量最高且最可靠的数据。因此，在不同样本之间进行比较时，必须确保各样本采用相同的测序策略制备，因为文库构建过程可能引入一定偏差。实验可重复性应始终通过设置重复样本来评估。
若干注意事项 
RNA 是一种不稳定的分子，容易因温度升高以及核糖核酸酶（RNases）的作用而降解。因此，建议采取特殊措施——使用个人防护装备（PPE）、无RNase的耗材以及RNase抑制剂。最重要的是，应尽可能将RNA保持在冰上。这一点尤其适用于甲基化RNA，因其对高温更为敏感。
必须确认目标RNA结构对DMS浓度和缓冲条件不敏感。在pH 7–7.5条件下，100 mM Tris、100 mM MOPS和100 mM HEPES等缓冲液可产生较强的信号，但可能不足以在反应过程中维持pH稳定21。由于DMS在水中会发生水解并降低pH，因此在修饰反应中使用强缓冲体系对于维持中性pH至关重要。已有研究表明，添加Bicine有助于将pH维持在微碱性状态21，但会导致鸟苷（G）和尿苷（U）的DMS修饰水平较低；这些数据可能具有信息价值，但由于其信号强度远低于腺苷（A）和胞苷（C），需单独分析，本实验方案中不再进一步讨论。
在基因特异性的 RT-PCR 中，修饰后的 RNA 被逆转录为 DNA，并通过 PCR 扩增成片段。尽管 RNA 的长度在理论上可以不受限制，但这些 PCR 片段的长度不应超过 400–500 个碱基对（bp），以避免逆转录反应过程中产生偏差。理想情况下，片段长度应处于测序读长范围之内（例如，若采用 150 × 150 循环的双端测序程序，则单一片段长度不应超过 300 bp）。当使用循环数较少的测序程序时，可使用 dsDNase 对 PCR 产物进行片段化处理。此外，由于引物序列内部不包含任何结构信息，因此当被检测的 RNA 包含 &#62;1 个片段时，这些片段必须相互重叠。逆转录反应中可包含针对不同片段的多个逆转录引物（最多可使用 10 种不同的逆转录引物）。根据序列的不同，将逆转录引物混合使用可能会降低逆转录效率，但通常仍能获得良好效果。每个 PCR 反应应单独进行。
使用DMS探测RNA时，实验条件具有额外的重要性，因为许多RNA在热力学上不稳定，并会根据温度等环境因素改变其构象。为避免出现异常，应尽可能保持实验条件恒定，反应时间也应保持一致。只要维持基本条件——缓冲能力以及单价离子（Na）和二价离子（Mg）的存在——以确保RNA正确折叠，缓冲体系在一定程度上可以互换17,20,22,23&#160;。24
关于修饰RNA的文库构建，必须考虑多个方面。首先，如前所述，修饰RNA比未修饰的RNA更不稳定，这意味着可能需要优化片段化时间以获得最佳的片段大小分布。此外，某些RNA文库构建试剂盒以及许多其他RNA测序方法在逆转录过程中使用随机引物，这可能导致参考序列的覆盖度降低，尤其是在基因的3'端，最终导致覆盖深度不足。如果某个区域的覆盖度过低，则可能需要在结构预测中剔除这些碱基。除了RT-PCR和全基因组RNA测序试剂盒外，还可采用其他文库构建方法。当使用小片段RNA或需要避免引物区域探针信息丢失时，采用将3'和/或5'接头连接至RNA的实验方案具有优势。
最后，对化学探针实验结果的分析必须始终谨慎解读。目前尚无任何软件能够仅根据序列就以高准确度预测任意RNA的二级结构。尽管化学探针数据的引入可显著提高预测准确性，但对较长的RNA（&#62;500 nt）构建可靠的结构模型仍然具有挑战性。这些模型应通过其他方法和/或定点突变实验进一步验证。RNA结构预测软件通常以最大化碱基配对数量为目标，因此对开放构象施加了显著惩罚，这可能无法准确反映RNA的真实折叠状态5。因此，所获得的结构模型应通过量化其与原始化学探针数据的一致性（例如使用AUROC指标）以及技术重复之间的一致性（例如使用mFMI指标）来进行评估，正如Lan等人的研究所示20。
理想情况下，应采用不同系统中的多个实验来验证所获得的结构模型，以加强研究假设。这些实验可包括体外（in vitro）和细胞内（in-cell）方法、补偿性突变，以及使用不同的细胞系和物种。此外，原始反应活性数据往往与结构预测同样重要，甚至更具信息量，因为它们记录了RNA折叠集合体的&#34;真实状态&#34;快照。因此，原始反应活性数据非常适合用于比较不同条件之间的结构变化，且具有很高的信息价值。重要的是，通过化学探针数据约束并结合计算预测所得到的最低自由能结构，仅应作为构建完整结构模型的初步假设。

自从发现RNA具有结构1,2和催化3功能以来，RNA在众多生物过程中所起的重要作用及其调控功能逐渐被揭示。事实上，RNA结构对基因调控的影响正受到越来越多的关注4。与蛋白质类似，RNA具有一级、二级和三级结构，分别指核苷酸序列、碱基配对相互作用的二维图谱以及这些碱基配对结构的三维折叠。尽管确定三级结构对于理解依赖RNA过程的确切机制至关重要，但二级结构同样能提供大量关于RNA功能的信息，并且是进一步形成三维结构的基础5。
然而，利用传统方法确定RNA结构一直存在内在挑战。对于蛋白质而言，X射线晶体学、核磁共振（NMR）和冷冻电子显微镜（cryo-EM）已能够解析其多样的结构模体，使得仅从序列即可进行结构预测6，但这些方法在RNA研究中应用有限。事实上，RNA是具有高度灵活性的分子，其结构单元（核苷酸）相比氨基酸具有更强的构象和旋转自由度。此外，RNA的碱基配对相互作用也比氨基酸残基之间的相互作用更加动态且多样化。因此，经典方法仅成功应用于结构明确、高度紧凑的较小RNA分子7。
确定RNA结构的另一种方法是通过化学探针结合下一代测序技术（NGS）。该策略可提供有关RNA序列中每个碱基结合状态（即其二级结构）的信息。简而言之，RNA分子中未参与碱基配对的碱基会被小分子化学试剂选择性修饰。利用特异的逆转录酶（RTs）对这些RNA进行逆转录，可将修饰信息以突变形式整合到互补DNA（cDNA）中。随后通过聚合酶链式反应（PCR）扩增这些cDNA分子并进行测序。为了获取目标RNA各碱基处于结合或未结合状态的“&#34;状态&#34;”信息，需计算每个碱基位点的突变频率，并将这些数据作为约束条件输入结构预测软件8。基于相邻碱基对规则9和最小自由能计算10，该软件生成最符合实验数据的RNA结构模型11,12。
DMS-MaPseq 使用 DMS，后者可甲基化腺苷的 N1 氮原子 &#160;N3 氮气中&#160;以高度特异性的方式在沃森-克里克面作用于胞嘧啶13使用耐热II组内含子逆转录酶（TGIRT-III）进行逆转录可产生前所未有的高信噪比突变谱，甚至能够解析由两种或多种不同构象重叠产生的突变谱14,15此外，DMS能够穿透细胞膜和完整组织，使得在生理环境条件下进行检测成为可能。然而，获得高质量数据具有挑战性，因为操作步骤中的差异可能影响实验结果。因此，我们提供了一个详细的方案，涵盖两方面 体外 以及细胞内 DMS-MaPseq 技术，以减少偏差，并帮助初学者克服该方法可能遇到的困难。特别是在近期 SARS-CoV-2 大流行的背景下，高质量的 RNA 病毒数据是研究基因表达和寻找潜在治疗方法的重要工具。

<table><tbody><tr><td>1 Kb Plus DNA Ladder</td><td>10787018</td><td>Thermo</td><td></td></tr><tr><td>2-巯基乙醇</td><td>M6250-250ML</td><td>Sigma</td><td></td></tr><tr><td>酸性酚:氯仿，pH 4.5</td><td>AM9720</td><td>Thermo</td><td></td></tr><tr><td>Advantage PCR</td><td>639206</td><td>Takara</td><td></td></tr><tr><td>CloneAmp HiFi PCR Premix</td><td>639298</td><td>Takara</td><td></td></tr><tr><td>DMS</td><td>D186309&lt;br/&gt; &amp;nbsp;</td><td>Sigma</td><td></td></tr><tr><td>dNTPs，每种10 mM</td><td>U151B</td><td>Promega</td><td></td></tr><tr><td>E-Gel EX 琼脂糖凝胶，2%</td><td>G402022</td><td>Thermo</td><td>预制琼脂糖凝胶</td></tr><tr><td>乙醇（200度）</td><td>E7023-4X4L</td><td>Sigma</td><td></td></tr><tr><td>Falcon 管，15 mL，50 mL</td><td></td><td></td><td></td></tr><tr><td>GlycoBlue</td><td></td><td></td><td>共沉淀剂</td></tr><tr><td>HCT-8 细胞&amp;nbsp;</td><td>ATCC #CCL-244</td><td></td><td></td></tr><tr><td>Invitrogen&amp;nbsp;MgCl&lt;sub&gt;2&lt;/sub&gt;&amp;nbsp;(1 M)</td><td>AM9530G</td><td>fisherscientific</td><td></td></tr><tr><td>异丙醇</td><td>278475</td><td>Sigma</td><td></td></tr><tr><td>Megascript T7 转录试剂盒</td><td>AM1334</td><td>Thermo</td><td></td></tr><tr><td>NanoDrop 分光光度计</td><td></td><td></td><td></td></tr><tr><td>Novex TBE 凝胶，8%，10孔</td><td>EC6215BOX</td><td>Thermo</td><td></td></tr><tr><td>OC43&amp;nbsp;</td><td>ATCC #VR-1558</td><td></td><td></td></tr><tr><td>RiboRuler 低分子量 RNA Ladder</td><td>SM1831</td><td>Thermo</td><td></td></tr><tr><td>RNase H</td><td>M0297L</td><td>NEB</td><td></td></tr><tr><td>二甲胂酸钠，0.4 M，pH 7.2</td><td>102090-964</td><td>VWR</td><td></td></tr><tr><td>氢氧化钠溶液</td><td>S8263-150ML</td><td>Sigma</td><td></td></tr><tr><td>SuperScript II 逆转录酶（用于FSB和DTT）</td><td>18064014</td><td>Thermo</td><td></td></tr><tr><td>TGIRT-III 酶</td><td>TGIRT50</td><td>Ingex</td><td></td></tr><tr><td>Oligo Clean &amp; Concentrator 试剂盒</td><td>D4060</td><td>Genesee</td><td></td></tr><tr><td>RNA Clean &amp; Concentrator 试剂盒为RNA纯化试剂盒</td><td>R1016</td><td>Genesee</td><td></td></tr><tr><td>TRIzol 试剂</td><td>15596018</td><td>Thermo</td><td>RNA提取试剂</td></tr><tr><td>水（用于RNA实验）（DEPC处理，无DNase、RNase，分子生物学级）</td><td>BP561-1</td><td>fisherscientific</td><td></td></tr><tr><td>xGen 广谱RNA文库构建试剂盒，16次反应</td><td>10009865</td><td>IDT</td><td></td></tr><tr><td>Zymo RNA纯化与浓缩柱</td><td></td><td></td><td></td></tr></tbody></table>

利用二甲基硫酸盐突变谱结合测序技术探究rna结构 体外 以及在细胞中

Probing RNA Structure with Dimethyl Sulfate Mutational Profiling with Sequencing In Vitro and in Cells

RNA结构在几乎所有生物学过程中的作用日益凸显，尤其是在过去十年中。然而，传统的RNA结构解析方法，如RNA晶体学或冷冻电镜（cryo-EM），已难以跟上该领域快速发展的步伐，也无法满足高通量分析的需求。基于测序的二甲基硫酸酯突变谱分析（DMS MaPseq）是一种通过检测碱基对DMS的反应活性来推断RNA结构的测序方法。当腺苷的N1位氮原子和胞嘧啶的N3位氮原子处于未配对状态时，DMS会对其沃森-克里克（Watson-Crick）面进行甲基化修饰。使用耐热的II型内含子逆转录酶（TGIRT-III）对经修饰的RNA进行逆转录，可使甲基化碱基在cDNA中以突变形式被引入。对所得cDNA进行测序并将其比对至参考转录本后，各碱基的相对突变率即可反映该碱基的配对状态 &#34;状态&#34; 作为成对或未成对。尽管DMS反应性具有较高的信噪比 体外 并且在细胞中，该方法对操作步骤中的偏差较为敏感。为减少此类偏差，本文提供了一种在细胞内使用DMS进行RNA处理的方案以及 体外 转录的RNA

基因特异性 体外 DMS-MaP 
为研究SARS2的5'非翻译区（5'UTR），合成了病毒前300个碱基对的gBlock序列，并设计了三对引物。其中包括用于扩增该片段的两对引物（“FW” &amp; “RV” 通过 PCR，以及一个用于连接T7启动子的序列（“FW-T7”）。这些序列可见于 表1.
<table border="1" fo:keep-together.within-page="1" fo:keep-with-next.within-page="always">
	<tbody>
		<tr>
			<td>名称</td>
			<td>序列 (5’-&gt;3’)</td>
		</tr>
		<tr>
			<td>FW</td>
			<td>ATTAAAGGTTTATACCTTCCCAGGTAAC</td>
		</tr>
		<tr>
			<td>RV</td>
			<td>GCAAACTGAGTTGGACGTGT</td>
		</tr>
		<tr>
			<td>FW-T7</td>
			<td>TAATACGACTCACTATAGG ATTAAAGGTTTATACCTTCCCAGGTAAC</td>
		</tr>
	</tbody>
</table>
表1：SARS-CoV2 5'UTR的DMS-MaP RT-PCR引物序列 此处需要 FW-T7 和 RV 来生成 DNA 模板 体外 转录过程中，RV用于逆转录，FW-RV引物对用于后续cDNA的PCR扩增。正向引物（FW）结合于SARS-CoV2基因组的起始位置，反向引物（RV）结合于目标区域下游紧邻的序列。缩写：DMS-MaP = 二甲基硫酸盐结合测序的突变谱分析；RT-PCR = 逆转录聚合酶链式反应；SARS-CoV2 = 严重急性呼吸综合征冠状病毒2；UTR = 非翻译区；RV = 反向引物；FW = 正向引物。
为了从gBlock片段生成RNA，根据图2A所示方案，使用PCR预混液通过重叠PCR连接T7聚合酶启动子序列。从延长后的片段出发，利用T7转录试剂盒生成RNA。随后使用DNase消化DNA模板，并通过RNA纯化与浓缩柱分离RNA。
质量控制的 体外 转录结果通过将RNA产物与单链RNA标记物一起在1%琼脂糖凝胶上电泳分析。由于仅观察到一条条带， 体外 DMS 探针法和 RT-PCR 实验已进行（参见 图2B).
为了验证PCR反应的成功，使用dsDNA Ladder在2%琼脂糖凝胶上对样品进行电泳。加接索引序列后，同一凝胶上条带应比原来高约150 bp，以计入索引引物的长度。
<img alt="PCR and IVT diagram with gel electrophoresis results for RNA synthesis and nucleotide separation." class="xfigimg" src="/files/ftp_upload/64820/64820fig2v2.jpg" data-altupdatedat="1762425657000" data-alt="Figure 2"> 
图2: 体外 DNA模板的转录 (A为了 体外 将不含固有RNA聚合酶启动子的DNA模板进行转录时，必须首先通过重叠PCR将启动子序列连接至模板。该过程使用一条上游引物（正向引物），其序列包含TAATACGACTCACTATAGG（以T7 RNA聚合酶为例），该序列位于与目标片段重叠的起始碱基上游。此处下划线标示的碱基代表聚合酶的转录起始位点。一旦启动子序列连接至双链DNA片段，即可被T7 RNA聚合酶转录。需要特别注意的是，该聚合酶以与上述启动子序列互补的链（蓝色）作为模板，从而生成与启动子序列下游紧邻区域序列相同（红色）的RNA。B1% 琼脂糖凝胶，含 ssRNA Ladder（第1泳道）以及 体外 300 nt 处的转录 RNA 产物（泳道 2）。(C2% 琼脂糖凝胶电泳图，第1道为GeneRuler 1 kb Plus Ladder，第2道为RT-PCR后在300 bp处的PCR产物，第3道为文库构建后在470 bp处的加接头片段。缩写：RT-PCR = 逆转录聚合酶链式反应；DMS = 二甲基硫酸酯；nt = 核苷酸；dsDNA = 双链DNA；ssRNA = 单链RNA。 <a href="https://www.jove.com/files/ftp_upload/64820/64820fig2v2large.jpg" target="_blank">请点击此处以查看此图的放大版本。</a>
全基因组 体内 使用病毒感染细胞进行 DMS-MaP 
在DMS处理之前，HCT-8细胞已感染OC43。感染后4天（dpi）观察到细胞病变效应（CPE）（如图所示） 图3A），这些细胞经过处理后提取RNA并进行了核糖体RNA去除。将总RNA在琼脂糖凝胶上电泳时，可见两条明亮条带，分别对应核糖体的40S和60S亚基，二者约占总RNA质量的95%（参见 图 3B）。当RNA提取失败或发生降解（例如，因多次冻融循环）时，RNA降解产物会在凝胶底部可见（参见 图3C，第二泳道）。此外，rRNA 耗尽后，两条明亮的条带消失，仅留下泳道中的弥散条带（见 图3C，第三条带）。最后，在文库构建完成后，各样本的片段大小分布不同，在最终的PAGE凝胶上呈现为弥散条带。根据计划用于分析这些文库的150 × 150双端测序方案，切取200核苷酸（nt）至500 nt之间的条带。最重要的是，位于约150 nt处的接头二聚体已被分离去除（见 图3D).
<img alt="Static equilibrium; gel electrophoresis results; ribosomal RNA analysis; diagrams showing subunit bands." class="xfigimg" src="/files/ftp_upload/64820/64820fig3v2.jpg" data-altupdatedat="1762425657000" data-alt="Figure 3"> 
图3：检查点 体内 病毒感染细胞中的 DMS-MaP。 (A) 病毒感染的HCT-8细胞在4 dpi时的光学显微镜图像。为了从总RNA中获得尽可能高的病毒RNA产量，同时尽量减少细胞死亡带来的不利影响，应在出现细胞病变效应（CPE）时或甚至更早时加入DMS，如图所示。B) 含六份总RNA样品（每份1 µg）的1%琼脂糖凝胶。每条泳道中均可见两条明亮条带，分别对应40S和60S亚基，因为核糖体RNA约占总RNA的95%。注意：细胞内DMS处理会导致部分RNA断裂和条带拖尾，但两条rRNA条带仍应可见。修饰后的轻度断裂是可以接受的，因为在活细胞中进行DMS孵育时，携带甲基化标记的信息已被生成并反映了RNA的结构状态。C1% 琼脂糖凝胶电泳图，泳道1为GeneRuler 1 kb plus DNA标记物，泳道2为在−80 °C保存6个月的总RNA，泳道3为去核糖体RNA。当RNA长期储存并经历多次冻融循环后，RNA开始降解，可能不适用于探针实验。此外，总RNA经去核糖体处理后，代表核糖体40S和60S亚基的两条明亮条带变淡，残留RNA的弥散条带开始显现。DGeneRuler 1 kb plus DNA 标记物（第1道）和全基因组RNA文库样本的PAGE凝胶电泳图。应根据测序需求切取凝胶中的目标片段。对于双端测序（从两端各进行150个循环），应切取300 bp至500 bp之间的片段。需将接头二聚体（在170 bp处迁移）分离去除。缩写：DMS-MaP = 二甲基硫酸盐结合测序的突变谱分析；dpi = 感染后天数；CPE = 细胞病变效应。 <a href="https://www.jove.com/files/ftp_upload/64820/64820fig3v2large.jpg" target="_blank">请点击此处查看此图的放大版本。</a>
测序完成后，将 .fastq 文件与一个 .fasta 参考序列文件一起提交至 DREEM 网络服务器（http://rnadreem.org/）进行分析。服务器生成的输出结果包括由 fastqc（https://www.bioinformatics.babraham.ac.uk/projects/fastqc/）和 TrimGalore（https://www.bioinformatics.babraham.ac.uk/projects/trim_galore/）产生的质控文件，以及包含群体平均突变频率的其他输出文件。除了一张以交互式 .html 格式展示突变频率的图表（见图 4A）、每个碱基的原始反应性数据的 .csv 文件和一个可被多种 RNA 结构预测软件读取的 struct_constraint.txt 文件外，还包括一个 bitvector.txt 文件，用于记录每条测序读段的突变情况。基于这些数据，将 .fasta 文件和 struct_constraint.txt 文件提交至 RNAfold 网络服务器（http://rna.tbi.univie.ac.at/cgi-bin/RNAWebSuite/RNAfold.cgi），以计算群体平均结构。该服务器利用 ViennaRNA 软件，基于最小自由能原理生成结构预测结果，预测结构可在线查看，也可下载为 ct 或 Vienna 格式。为了构建 RNA 结构模型，将上述可下载的文件提交至 VARNA 软件（https://varna.lri.fr/，见图 4B）。最后，可使用 DREEM 的稳定版本（https://codeocean.com/capsule/6175523/tree/v1）分析 bitvector.txt 文件，以搜索替代性的 RNA 构象。为了利用 DREEM 获得高质量的结构模型，建议每个碱基的测序深度达到 10,000 条读段；而对于聚类分析，每个碱基可能需要高达 100,000 条读段。整个工作流程的概览见图 4C。
<img alt="RNA 反应性分析：图表、二级结构图和测序工作流程图。" class="xfigimg" src="/files/ftp_upload/64820/64820fig4v2.jpg" data-altupdatedat="1762425657000" data-alt="Figure 4"> 
图 4：SARS-CoV2 5'非翻译区（5'UTR）化学探针实验获得的代表性数据。（A）SARS-CoV2 基因组前 300 个碱基的反应性图谱，按碱基着色（A：红色，C：蓝色，U：绿色，G：黄色）。原始反应性值通过绝对突变频率除以测序深度计算得出。具有开放构象的碱基反应性值较高；参与碱基配对的碱基反应性值较低。U 和 G 不被 DMS 修饰，其低反应性值源于聚合酶的错配错误。预测由 DREEM 网络服务器完成。（B）基于反应性值使用 VARNA 预测的 SARS-CoV2 5'UTR 结构模型。高反应性值的碱基以红色表示；低反应性值的碱基以白色表示。（C）DMS-MaP 分析的工作流程，起始于测序获得的 .fastq 文件。可使用 fastqc 进行质量控制；利用 TrimGalore 去除接头序列后，使用 Bowtie2 将序列比对至参考序列。从获得的 .bam 文件中，DREEM 统计每条读段中的突变，生成突变图谱或 .bitvector.txt 文件。这些文件以位置依赖的方式报告每条读段的突变情况，据此可构建群体平均反应性图谱。此外，也可使用 DREEM 对 bitvector 进行聚类，以寻找 RNA 的替代构象。最后，使用软件（例如 VARNA）对获得的结构模型进行可视化展示。缩写：DMS-MaP = 使用二甲基硫酸盐的突变谱型测序；SARS-CoV2 = 严重急性呼吸综合征冠状病毒 2。 <a href="https://www.jove.com/files/ftp_upload/64820/64820fig4v2large.jpg" target="_blank">请点击此处查看该图的放大版本。</a>

该方案提供了使用二甲基硫酸酯对RNA进行修饰以开展突变谱分析实验的操作指导。内容包括 体外 和 体内 采用两种不同的文库构建方法进行检测。

利用二甲基硫酸盐突变谱结合测序技术探究RNA结构 体外 以及在细胞中

RNA 结构在几乎所有生物过程中的作用在过去十年中变得愈发明显。然而，经典……

DMS-MaP 是一种基于测序的方法，可让我们获取 RNA 结构的快照，并了解其在不同条件下的变化情况。与传统的结构测定方法（如晶体学和电子显微镜技术）相比，DMS-MaP 可在细胞内使用，并能解析 RNA 结构集合体。由于 RNA 结构在多种生物学现象中具有重要意义，我们的方法几乎可为所有生物学研究领域提供机制层面的深入见解。每个终体积为 100 微升的反应，先将 89 微升复性缓冲液转移至指定的 1.5 毫升离心管中。将离心管置于化学通风橱下的恒温振荡器中，在 37 摄氏度预热。向一个 PCR 管中加入 10 微升无核酸酶水，并将 1 至 10 皮摩尔的 RNA 转移至其中。将样品置于PCR仪中，在95摄氏度下孵育1分钟，以使RNA变性。随后立即将离心管放置在冰块上，防止RNA错误折叠。接下来，为了使RNA重新折叠，将样品加入指定的含有37摄氏度重折叠缓冲液的离心管中，充分混匀后，在37摄氏度下孵育10至20分钟。接下来，向含有RNA样品的离心管中加入1 µL浓度为10.5 M的100%二甲基硫酸盐（DMS），并将反应混合物在800至1,400 rpm的振荡条件下孵育5分钟。反应5分钟后，加入60 µL 100% β-巯基乙醇终止反应，充分混匀，涡旋振荡后立即将RNA置于冰上。随后纯化RNA，并用10 µL无核酸酶的水洗脱。使用分光光度计定量RNA。最后，将修饰后的RNA储存在-80摄氏度。将反应混合物加入PCR管后，将10微升无核酸酶水中至少100纳克的洗脱修饰RNA转移至PCR管中。将混合物置于热循环仪中，在 57 摄氏度下孵育 30 分钟至 1.5 小时。30 分钟足以生成含有 500 个核苷酸的产物。接着，向其中加入 1 µL 4 M 氢氧化钠，用移液器充分混匀，并在 95 摄氏度下孵育 3 分钟，以降解 RNA 并使 TGIRT 从互补 DNA 上释放。采用基于层析柱的方法去除引物，纯化互补DNA，并通过PCR进行扩增。在进行文库构建之前，取2 µL PCR产物在琼脂糖凝胶上电泳，以验证PCR扩增是否成功。然后使用分光光度计对提取的片段进行定量，再对片段进行索引以用于测序。将培养至所需感染阶段的病毒感染细胞转移至具备相应生物安全级别的专用通风橱中。向含有细胞的培养基中加入2.5%体积的DMS，并用封口膜密封容器。立即将容器在37摄氏度下孵育5分钟。孵育结束后，小心地用移液器吸出含DMS的培养基，并弃入适当的化学废弃物容器中。然后轻轻向细胞中加入10毫升终止缓冲液，刮下细胞，并将其转移至15毫升锥形管中。将离心管在 3,000 g 条件下离心 3 分钟，使细胞沉淀。弃去上清液，并用 10 毫升 PBS 将细胞沉淀洗涤两次。每次洗涤后，尽可能彻底地移除残留的 PBS。将沉淀溶解于适量的RNA提取试剂中，以进行RNA提取。接着，向含有1毫升匀浆细胞的RNA提取试剂中加入200微升氯仿，涡旋振荡15至20秒，直至细胞溶液变为亮粉色。等待分相完成，此时可见粉红色的脂质相沉于底部。在4°C条件下，以约20,000 g的转速离心15分钟，然后将上层水相转移至新离心管中。纯化RNA，并用足量无核酸酶的水洗脱。随后采用选定的方法进行核糖体RNA去除，再用适量无核酸酶的水洗脱RNA，并使用分光光度计对其进行定量。纯化的RNA可再次用于RT-PCR。利用上述网络服务器，可方便地分析DMS-MaP数据。将测序读段比对至参考序列，并统计每个碱基位点的突变数。所得的群体平均文件可用作已知RNA结构预测软件中的约束条件。所得的最小自由能结构可使用VARNA等软件进行可视化。目前，只有DREAM的稳定版本能够解析RNA结构集合。然而，目前正在进行相关工作，以使该功能对整个科研群体更加易于使用。通过在gBlock片段上连接T7聚合酶启动子进行体外转录，获得目标RNA，在1%琼脂糖凝胶上显示为一条约300个核苷酸的单一条带。对经DMS修饰的RNA进行基因特异性RT-PCR，成功的结果表现为在2%琼脂糖凝胶上出现一条约300个碱基对的单一条带。如预期所示，带索引的片段位置高出约150个碱基对。在病毒感染后第四天出现细胞病变效应的HCT-8细胞上进行了DMS处理。提取的总RNA在琼脂糖凝胶上显示为两条明亮条带，分别对应40S和60S亚基。核糖体RNA耗尽后，两条明亮的条带消失，相应泳道中仅留下弥散条带。文库构建后，各样本的片段大小分布各异，均呈现为弥散条带。切取200至500个核苷酸之间的条带，并将接头二聚体分离去除。新冠病毒基因组的结构模型显示，处于开放构象的碱基比参与碱基配对的碱基具有更高的反应活性。碱基的反应活性图谱表明，尿嘧啶和鸟嘌呤的反应活性值较低，且未被DMS修饰。通过我们的方法，能够以高通量的方式在细胞内预测RNA结构，且不受片段大小限制。对照指标非常重要，因为它们能够提供有关结构预测准确性的信息。为了进一步验证预测的准确性，应开展更多干扰或改变结构的实验。

probing-rna-structure-with-dimethyl-sulfate-mutational-profiling-with

Research

科研

JoVE Journal

Biochemistry

生物化学

使用硫酸二甲酯探测RNA结构 体 外 和细胞测序突变分析

本方案描述了用于Entamoeba histolytica...

体外 tRNA 甲基化实验：使用溶组织内阿米巴 DNA 和 tRNA 甲基转移酶 Dnmt2 (Ehmeth) 酶

结合特定RNA序列的蛋白质在基因表达中发挥关键作用。对这些结合位点进行详细表征，对于我们理解基因调控机制至关重要。本文描述了一种针对RNA中蛋白质结合位点的单步饱和突变方法。该方法适用于RNA中所有蛋白质结合位点的研究。...

一种新型饱和突变方法：利用硫代磷酸酯一步法表征RNA中调控蛋白结合位点

核糖核苷酸是在真核生物细胞核DNA复制过程中掺入基因组中最丰富的非经典核苷酸之一。如果未能被及时清除，核糖核苷酸可能导致DNA损伤和诱变。本文介绍了两种用于评估DNA中核糖核苷酸掺入丰度及其诱变效应的实验方法。...

研究核糖核苷酸的掺入：酵母基因组中核糖核苷酸的链特异性检测及核糖核苷酸诱导突变的测量

本文介绍了一种对酵母RNA修饰酶Mod5进行生化表征的实验方案，并讨论了该方案如何应用于其他RNA修饰酶的研究。...

一个 体外 检测tRNA-异戊烯基转移酶活性的实验方法

定点突变是一种用于在脱氧核糖核酸（DNA）中引入特定突变的技术。本实验方案介绍了基于两步法和三步法聚合酶链式反应（PCR）的定点突变操作方法，该方法适用于任何感兴趣的DNA片段。...

位点定向诱变在体外和体内实验中的应用——以大肠杆菌（Escherichia coli）中的RNA相互作用为例

本文介绍了在存在靶向RNA的小分子的情况下，通过引物延伸分析法（SHAPE）进行体外和细胞内选择性2'-羟基乙酰化的详细实验方案，以确定pre-mRNA目标序列的二级结构。...

利用体外和细胞内SHAPE技术研究小分子诱导的pre-mRNA结构变化

本文介绍一种依赖脱氧核酶（DNAzyme）的RNA切割实验方案。该方法可实现对RNA...

DNA酶依赖性分析rRNA的2'-O-甲基化

本方案描述了单分子荧光原位杂交（smFISH）技术的一种应用，用于测量体内mRNA合成与降解的动力学过程。...

在 Escherichia coli 中利用双色单分子荧光原位杂交技术研究 mRNA 的时空动力学

RNA二级结构主要通过结构探针技术在成熟RNA中进行观察。共转录结构追踪测序技术（CoSTseq）将已用于研究聚合酶在新生RNA上位置的核内延伸技术与结构探针技术相结合，从而实现了对处于活跃转录过程中的RNA的二级结构的观测。...

新生转录本与成熟转录本的RNA结构比较分析 酿酒酵母（Saccharomyces cerevisiae）

本方案描述了一种利用分裂型 Broccoli RNA 报告系统进行的可视化检测方法，用于检测和定量 RNA 簇中的分子间碱基配对 体外.

分子间RNA碱基配对的可视化与定量分析 体外 利用分裂型 Broccoli RNA 报告系统检测 RNA 簇

学习使用硫酸二甲酯进行RNA结构探测的DMS-MaP方法，通过体外和细胞测序突变分析揭示RNA动态变化。该技术为生物学研究提供高通量结构预测工具。

利用二甲基硫酸盐突变谱结合测序技术探究RNA结构体外以及在细胞中

探索更多视频

利用二甲基硫酸盐突变谱结合测序技术探究RNA结构 体外 以及在细胞中

探索更多视频

利用二甲基硫酸盐突变谱结合测序技术探究RNA结构体外以及在细胞中