方法文章

临床样本中RNA病毒的无偏深度测序

17.2K 次观看

⸱

DOI：

10.3791/54117

⸱

2016年7月2日

 ,  ,  ,  ,  ,  , 
1Broad Institute of MIT and Harvard, 2Harvard University

本文内容

摘要

本方案描述了一种快速且广泛适用的方法，用于对来自人类临床分离物的病毒样本进行无偏倚的RNA测序。

摘要

本文概述了一种能够实现下一代RNA测序的实验方案 从头合成 从临床和生物来源收集的病毒基因组的组装及宿主内变异检测。该方法具有无偏性和普适性，采用随机引物进行cDNA合成，无需预先了解病毒序列信息。在文库构建前，利用基于RNase H的选择性消化方法去除病毒RNA样本中不需要的RNA，包括聚(rA)载体和核糖体RNA。选择性去除可提高病毒RNA测序文库的数据质量和唯一读段数量。此外，本方案采用基于转座酶的“标签片段化”（tagmentation）步骤，可缩短文库构建的总体时间。该方案已成功实现对600多个拉沙病毒和埃博拉病毒样本（包括血液和组织分离物）的快速深度测序，且广泛适用于其他微生物基因组学研究。

引言

来自临床样本的病毒下一代测序可用于研究感染的传播途径和流行病学特征，并有助于推动新型诊断方法、疫苗和治疗手段的研发。利用随机引物进行cDNA合成，已实现对差异较大、共感染或甚至 novel 病毒的基因组检测与组装1,2。与其他无偏倚方法一样，非目标污染物会占据大量测序读长，对测序结果产生负面影响。宿主核酸和聚(rA)载体RNA是现有许多病毒样本库中存在的主要污染物。

本方案描述了一种基于无偏倚的全RNA测序（total RNA-seq）技术高效且经济地对RNA病毒基因组进行深度测序的方法。该方法利用RNase H选择性去除步骤3，以清除不需要的宿主核糖体RNA和载体RNA。选择性去除可富集病毒序列内容（图1），并提升临床样本测序数据的整体质量（图2）。此外，本方案引入了标签化片段化（tagmentation）步骤，显著缩短了文库构建时间。这些方法已被用于快速生成埃博拉病毒和拉沙病毒基因组的大规模数据集2,4,5，并可广泛应用于多种RNA病毒的研究。最后，该方法不仅限于人类样本；已有研究在感染拉沙病毒的啮齿类动物及非人灵长类疾病模型的组织样本中验证了选择性去除策略的有效性5,6。

RNA输入分析；18S rRNA拷贝数、LASV富集、临床分离株评估示意图。
图1. 总RNA含量反映了通过选择性耗竭富集拉沙病毒（Lassa Virus, LASV）的效果。 对九个不同临床分离株进行rRNA耗竭前后，初始总RNA含量（RNA输入）及独特拉沙病毒（LASV）序列读数（文库含量）的变化情况。本图改编自6。 请点击此处查看该图的放大版本。

显示碱基质量分数的测序循环图；poly(rA) 影响比较。
图 2. 去除载体 RNA 后测序质量更高。来自质控报告13的受 poly(rA) 污染的拉沙病毒文库（红色）与对照文库（文库中未检测到载体，黑色）在每个测序循环中的中位碱基质量。双端测序的 Read 1 和 Read 2 已合并至文库 BAM 文件中，每个碱基位置的质量分数均显示于此。该图已根据文献6修改。请点击此处查看此图的放大版本。

病毒RNA测序方案详细描述了直接从临床和生物样本中提取的RNA构建文库的步骤。为确保个人安全，所有病毒血清、血浆和组织样本在进行RNA提取前，应先在适当的缓冲液中灭活。某些灭活和提取试剂盒中含有载体多聚(rA) RNA；该成分将在初始的RNase H选择性耗竭步骤中被去除。根据完全回收的预期，载体RNA的浓度为100 ng/µl。在本方案中，使用110 ng/µl的寡核苷酸dT RNA（即1.1倍载体浓度）进行耗竭。如果样本中不含多聚(rA)载体，则在耗竭前不应添加寡核苷酸(dT)。

以下方案适用于 24 个反应的 PCR 板格式（最多 250 µl 体积）。该方案的早期版本已在 Matranga 等人中报道。 等6.

方案

伦理声明：拉沙热患者是根据杜兰大学、哈佛大学、布罗德研究所、尼日利亚伊瑞瓦专科教学医院（ISTH）、凯内马政府医院（KGH）、尼日利亚伊巴丹奥约州卫生部以及塞拉利昂卫生部人类受试者委员会批准的方案招募进入本研究的。所有患者均接受同等标准的治疗护理，并被提供抗病毒药物利巴韦林，无论其是否决定参与本研究。对于拉沙热（LF）患者，利巴韦林的治疗遵循目前推荐的治疗指南，通常在高度怀疑患有拉沙热后立即给予。

由于埃博拉病毒病（EVD）疫情严重暴发，无法通过我们的标准方案获得患者的知情同意。因此，使用在塞拉利昂和哈佛大学的机构审查委员会评估并批准的EVD患者临床剩余样本。塞拉利昂伦理与科学审查委员会办公室、塞拉利昂卫生与环境卫生部以及哈佛大学人类受试者使用委员会已批准豁免知情同意，允许对在塞拉利昂埃博拉疫情中采集的患者及密切接触者样本中获得的病毒序列进行测序并公开发布。这些机构还批准使用在疫情应对期间接受治疗的所有疑似EVD患者所采集的去标识化样本的临床和流行病学数据。塞拉利昂卫生与环境卫生部还批准将非传染性、非生物性样本从塞拉利昂运往布罗德研究所和哈佛大学，用于疫情样本的基因组研究。

1. 样品RNA的DNase处理（最多55 µl提取的总RNA，约4小时）

  1. 在生物安全柜中的冰上，按照所述方法在96孔PCR板中设置DNase反应 表1，步骤1.1（总体积， 70 µl/well）。注意：可配制预混液。
  2. 轻轻且充分涡旋振荡，然后在室温下以280 x g离心1分钟。
  3. 孵育于 37 °C 30分钟
  4. 使用 RNA 固相可逆固定（SPRI）磁珠进行纯化。
    1. 将RNA beads在RT下孵育30分钟。
    2. 轻轻摇动RNA磁珠瓶，以重悬可能已沉降的磁性颗粒。加入1.8倍体积（126 µl) RNA珠与经DNase处理的RNA的比值70 µl），用移液器吹打混合10次，在室温下孵育5分钟（孔中总体积， 196 µl).
    3. 将混合物置于磁力架上，等待溶液澄清（5 - 10 分钟）。
    4. 移除站上的上清液，用移液器吸出并弃去。在站上操作，用70%乙醇覆盖沉淀以洗涤磁珠，孵育1分钟。用移液器吸出乙醇并弃去。共重复洗涤两次。
      注意：使用精确配制的70%新鲜乙醇至关重要，因为浓度过高将导致对较小分子的洗涤效率降低，而 <70% 乙醇可能导致样品损失7.
    5. 将培养板留在工作台上并敞开放置以晾干。注意：务必让磁珠充分干燥，直至磁珠表面开始出现裂纹。
    6. 添加 55 µl 加入无核酸酶水至孔板中以洗脱RNA。从仪器上取下孔板，用移液器充分吹打混合磁珠与水。注意：也可使用较少的水量（≤ 10 µl) 以浓缩总 RNA。
    7. 将培养板放回工作站。待溶液澄清后，用移液器转移至新的螺旋盖管中，用于长期保存（-80 °C）。放置 5 µl 用于去除的新的96孔PCR板中的RNA（步骤2.4）。
    8. 可选：保存并稀释 1 µl 在 19 µl 水（1:20），用于rRNA的qRT-PCR例如，18S，28S rRNA) (表 2）和病毒标志物5.

2. 从病毒RNA样本中选择性去除核糖体RNA和载体RNA（约4小时）

  1. 按照 表1 中所述方法配制含线性丙烯酰胺载体的5倍杂交缓冲液、10倍RNase H反应缓冲液以及无核酸酶水。
  2. 在冰上，于96孔PCR板中将RNA与rRNA去除寡核苷酸（表3）及寡核苷酸(dT)混合，建立杂交反应，具体如 表1 所述。
    注意：可预先配制总反应混合液。可加入50飞克（fg）一种独特的合成RNA（ERCCs8），用于追踪病毒测序过程及潜在的索引读段交叉污染。
    1. 轻轻涡旋充分混匀，然后在室温下以280 × g离心1分钟。
      在95 °C孵育2分钟，以每秒-0.1 °C的缓慢斜率降温至45 °C。将热循环仪暂停在45 °C。
  3. 按照 表1 中所述方法在冰上配制RNase H反应混合液，然后在45 °C预热2分钟。注意：可预先配制总反应混合液。
    1. 在保持PCR板处于45 °C热循环仪中的条件下，将预热的RNase H反应液加入杂交反应体系中。
    2. 通过轻柔吹打6–8次充分混匀，于45 °C继续孵育30分钟。随后置于冰上。
  4. 按照 表1 中所述方法在冰上配制DNase反应混合液。注意：可预先配制总反应混合液。
    1. 将DNase反应液加入PCR板中的RNase H反应体系中，轻轻涡旋充分混匀，然后在室温下以280 × g离心1分钟。于37 °C孵育30分钟。
    2. 加入5 µl 0.5 M EDTA终止DNase反应。轻轻涡旋充分混匀，然后在室温下以280 × g离心1分钟。
  5. 使用RNA磁珠（见步骤1.3）进行纯化，加入1.8倍体积（144 µl）磁珠。用11 µl无核酸酶水洗脱。注意：为安全冷藏保存，建议将去除rRNA后的RNA样本在-80 °C过夜保存。

3. cDNA 合成（约 6 小时）

  1. 按照表1所述，在冰上将rRNA/载体耗尽的RNA与随机引物在96孔PCR板中混合，轻轻涡旋并充分混匀，然后在室温下以280 × g离心1分钟。
    1. 在热循环仪中于70 °C加热混合物10分钟。热变性后立即置于冰上冷却1至5分钟。在进行第一链合成反应前，RNA在冰上静置时间不得超过5分钟。
  2. 按照表1所述，在冰上配制第一链合成反应混合液。
    注意：可预先配制总混合液（master mix）。
    1. 将其加入板中的RNA/随机引物混合液中，轻轻涡旋并充分混匀，然后在室温下以280 × g离心1分钟。在22 - 25 °C孵育10分钟。
    2. 在空气孵育箱中于55 °C孵育60分钟。将板置于冰上以终止反应。注意：推荐使用空气孵育箱，以便在第一链合成反应过程中实现缓慢升温，使引物充分退火并启动第一链的延伸。
  3. 按照表1所述，在冰上配制第二链合成反应混合液。
    注意：可预先配制总混合液（master mix）。
    1. 将其加入板中的第一链合成反应体系中，轻轻涡旋并充分混匀，然后在室温下以280 × g离心1分钟。在16 °C孵育2小时（保持盖温为25 °C）。温度不得超过16 °C。
    2. 将板置于冰上，然后加入5 µl 0.5 M EDTA以终止反应，轻轻且充分混匀，再在室温下以280 × g离心1分钟。
  4. 使用DNA磁珠进行纯化（纯化步骤参见1.3），加入1.8倍体积（153 µl）的磁珠。用9 µl洗脱缓冲液（EB）洗脱。保留1 µl用于定量。后续步骤使用1 ng cDNA。若cDNA浓度过低无法检测，则使用4 µl cDNA进行片段化（参见步骤4.1）。
  5. 为安全冷藏，可将双链cDNA短期保存于4 °C过夜，或长期保存于-20 °C。

4. 文库制备——DNA文库构建（约4小时）

  1. 转移 4 µl 将cDNA加入96孔板中，并将剩余的cDNA保存，以备需要时进行第二次尝试。
  2. 在冰上按照以下所述设置转座反应 表1.
    注意：可预先配制混合液。为降低背景信号并减少总体成本，将片段化反应的总体积从 20 减少至 10 µl. 由于cDNA是限制因素，ATM的量（即， 反应中所用的转座子（transposome）用量也相应减少，以降低整合位点的数量。
    1. 将标签化混合物加入到平板中的 cDNA 中，轻轻涡旋混匀，280 × g（室温）离心 1 分钟。孵育于 55 °C 5 分钟，保持 10 °C.
    2. 一旦到达 10 °C立即加入 2.5 µl 加入中和片段化缓冲液（NT）以终止反应。用移液器吹打混匀，然后在室温下以280 × g离心1分钟。
    3. 室温孵育5分钟。
  3. 在冰上按照所述方法设置PCR扩增反应 表1.
    1. 轻轻且充分涡旋混匀，然后在室温下以 280 × g 离心 1 分钟。
    2. 使用所述条件在热循环仪中进行PCR 表1.
      注意：对于1 ng的经转座酶处理的cDNA，建议进行12个PCR循环；然而，病毒临床样本中的cDNA含量通常极低甚至无法检测。当cDNA含量较低时（<1 ng），可使用最多18个PCR循环以生成足够的文库用于测序。
  4. 文库制备——测序前的纯化与文库合并
    1. 将样品补至 50 µl 与EB一起。
    2. 使用0.6倍体积的DNA磁珠进行纯化（纯化步骤见1.3）30 µl)珠子。洗脱于 15 µl EB
    3. 确定文库浓度（图3通过使用生物分析仪软件进行区域分析（150 至 1,000 bp）9，排除引物二聚体（约120 bp）对区域分析的干扰。注意：也可采用qPCR对文库进行定量 10.
    4. 以不低于1 nM的最低摩尔浓度合并文库。若文库浓度低于1 nM，可加入少量该文库样本（约为其他文库体积的1倍）进行合并，以获取这些文库的序列信息。
    5. 使用0.7倍体积的DNA磁珠对混合物进行纯化，步骤同上（见步骤2）。洗脱于 15 µl EB。注意：磁珠体积取决于混合液的最终体积。
    6. 分析池9. 通过区域分析（150 至 1,000 bp）确定摩尔浓度9注意：也可使用qPCR来定量文库池 10.
    7. 以10 pM文库浓度上机测序，生成101 bp双端双标签测序读长11.

用于EBOV分析的cDNA文库与引物二聚体的琼脂糖凝胶电泳。
图3. 从埃博拉病毒临床样本构建的文库。 4个代表性埃博拉病毒（EBOV）文库的凝胶图像。图中显示了文库片段和引物二聚体的区域。 请点击此处查看该图的放大版本。

结果

本方案能够从低起始量的病毒RNA样本中生成高质量的测序读段，同时富集独特的病毒序列内容。如图1所示，在所有样本中（与未去除的对照相比），该方案对拉沙病毒独特序列的富集倍数至少达到五倍，即使样本中仅含有约一百万个18S rRNA拷贝（约100 pg总RNA）。同样，测序成功率也与特定样本中的病毒含量呈正相关。以qRT-PCR作为病毒量的替代指标，含有约1,000个或更多病毒基因组拷贝的样本通常能够完成完整的基因组组装（数据未显示）。此外，去除聚(rA)载体可减少文库中A和T的同聚物序列，从而获得更洁净的文库制备产物，并确保更高质量的测序读段（图2）。来自低起始量病毒临床样本的最终文库，其片段长度通常较宽，范围在150至1,000 bp之间（图3）。

测序完成后，为降低样本误识别以及文库池中各文库之间的交叉干扰12，仅保留去复用过程中索引读段碱基质量分数为25（q25）且无任何错配的读段。病毒基因组的组装采用专为高度变异病毒设计的生物信息学分析流程2,4-6。这些工具可通过 https://github.com/broadinstitute/viral-ngs 获取，或通过商业云平台获得4。

步骤 1.1：DNase 反应
试剂每反应体积（µl)
10x DNase 缓冲液7
无核酸酶水6
提取的病毒RNA55
DNase（2 U/µl)2
总体积70
步骤 2.1：5× 杂交缓冲液
试剂1 ml 体积（µl)
5 M NaCl200
1 M Tris-HCl (pH 7.4)500
无核酸酶水300
总体积1,000
第2.1步：10x RNase H反应缓冲液
试剂1 ml 体积（µl)
5 M NaCl200
1 M Tris-HCl（pH 7.5）500
1 M MgCl2200
无核酸酶水500
总体积1,000
第2.1步：含线性丙烯酰胺的水
试剂1 ml 缓冲液的体积（µl)
无核酸酶水992
线性丙烯酰胺（5 mg/ml）8
总体积1,000
步骤 2.2：用于选择性去除的杂交反应
试剂每反应体积（µl)
5x 杂交缓冲液2
rRNA 耗竭寡核苷酸混合物（100 µM)1.22
Oligo(dT)（550 ng/μL）µl)1
经DNase处理的总RNA 最多 5
外源添加的RNA（可选）0.5
水（含线性丙烯酰胺）加至总体积 10
总体积10
第2.3步：用于选择性去除的RNase H反应
试剂每反应体积（µl)
10x RNase H反应缓冲液       2
水（含线性丙烯酰胺）5
耐热RNase H（5 U/µl)3
总体积10
步骤 2.4：选择性耗竭后的 DNase 反应
试剂每反应体积（µl)
10x DNase 缓冲液7.5
水（含线性丙烯酰胺）44.5
RNase 抑制剂（20 U/µl)1
RNase-free DNase I (2.72 U/µl) 2
总反应体积（含 RNase H 反应）75
步骤 3.1：  cDNA合成，随机引物杂交
试剂每反应体积（µl)
去除rRNA/载体的RNA10
3 µg 随机引物1
总体积11
第3.2步：第一链cDNA合成反应
试剂体积µl)
5x 第一链反应缓冲液4
0.1 M DTT2
10 mM dNTP 混合液1
RNase 抑制剂（20 U/µl)1
逆转录酶（最后添加）1
总体积（含上述 RNA）20
步骤 3.3：第二链 cDNA 合成反应
试剂体积µl)
无RNase水43
10x 第二链反应缓冲液8
10 mM dNTP 混合液3
E. coli DNA连接酶（10 U/μL）µl)1
E. coli DNA聚合酶I（10 U/µl)4
E. coli RNase H（2 U/μL）µl)1
总体积（含1st 链反应80
步骤 4.2：  标签化反应
试剂体积µl)
扩增子片段化混合液（ATM）1
Tagment DNA 缓冲液（TD）5
总体积（含 cDNA）10
第4.3步：文库PCR反应
试剂体积µl)
PCR预混液（NPM）7.5
Index 1 引物（i7）2.5
Index 2 引物（i5）2.5
总反应体积（含经标签化处理的cDNA）25
第4.3.2步：文库PCR反应条件
72 °C，3 分钟
95 °C，30秒
最多18个循环——95℃下10秒 °C，55℃下30秒 °C，72℃下30秒 °C
72 °C，5 分钟
10 °C，永远

表1：反应体系设置与缓冲液。 包含所有缓冲液和反应混合物成分的逐步说明表格。

引物名称序列 (5' 至 3')
Ebola KGH FWGTCGTTCCAACAATCGAGCG
Ebola KGH RVCGTCCCGTAGCTTTRGCCAT
Ebola KULESH FWTCTGACATGGATTACCACAAGATC
Ebola KULESH RVGGATGACTCTTTGCCGAACAATC
Lassa SL FWGTA AGC CCA GCD GYA AAB CC
Lassa SL RVAAG CCA CAG AAA RCT GGS AGC A
18S rRNA FWTCCTTTAACGAGGATCCATTGG
18S rRNA RVCGAGCTTTTTAACTGCAGCAACT

表2：qRT-PCR引物序列。 用于检测宿主（18S rRNA）和病毒（埃博拉病毒和拉沙病毒）含量的引物。“KGH”指塞拉利昂的凯内马政府医院，埃博拉病毒引物在此进行了测试2。“Kulesh”为设计该引物组的研究人员14。

表3：核糖体RNA（rRNA）去除寡核苷酸。195条长度为50个核苷酸 、与人类rRNA互补的序列，用于选择性去除步骤6。请点击此处下载该文件。

讨论

该方案具有稳健性、通用性和快速测序的特点，在2014年埃博拉病毒暴发期间已被用于病毒测序2,4。通过将选择性耗竭和cDNA合成与基于标签片段化（tagmentation）的文库构建相结合，整体流程时间比此前的接头连接方法缩短了约2天。最近，该方案已被国际合作者及其他研究团队成功应用15,16，并将被推广至西非的实验室，以支持当地的基因组学研究和诊断工作17。

本方案采用随机引物制备用于病毒RNA测序文库的cDNA。与以往的病毒RNA测序方法不同，该方法无需预先了解序列信息，也无需针对特定病毒或进化枝设计复杂且耗时的引物。该方法可应用于任何病毒RNA样本。例如，该方法已成功用于埃博拉病毒和拉沙病毒样本的病毒组分分析6。此外，该方案也可用于宿主转录组、宏基因组及病原体发现测序项目1。

该实验方案的关键步骤是靶向RNase H消化，这是一种高通量、低成本的方法，用于从病毒样本中去除不需要的载体和宿主RNA。方案中的选择性耗竭步骤涉及多种组分，需要较高的操作技巧和准确性。在初始设置阶段应额外投入时间并谨慎操作。

由于大多数临床血清和血浆样本中的核酸物质含量极低，污染和样本损失较为常见。为避免这些问题，在使用本方案时应格外谨慎。首先，RNA极易降解，因此所有操作区域应保持清洁且无核酸酶。其次，为筛选适用于本方案的样本，应采用针对宿主RNA和病毒的qRT-PCR检测进行定量5,6。在将输入量与本方案的测序结果进行比较时发现，测序成功（即 获得足够数据以完成病毒全基因组组装）与样本中总RNA含量至少达到100 pg且病毒拷贝数至少为1,000相关。第三，应避免样本接触环境中外源性核酸来源。此处所述方案应在生物安全柜中进行，以确保操作安全并减少环境污染物的引入。此外，我们课题组及其他研究团队已注意到，商用酶类可能是低输入样本中污染性细菌核酸的另一来源6,18。使用洁净的工作空间（例如，PCR操作台、生物安全柜）以及设置阴性对照（例如，水或缓冲液）将分别有助于减轻污染并追踪污染来源。对于总RNA含量<100 pg的样本，仅应去除poly(rA)载体RNA，而不应去除rRNA，以确保获得高质量的测序结果，同时尽量减少样本损失。对于输入量极低的样本，cDNA扩增方法可能更为适用19，但应在cDNA合成前去除poly(rA)载体。

去除宿主rRNA可富集测序文库中的病毒成分，适用于多种样本类型，包括血清或血浆，以及啮齿类动物和非人灵长类动物的多种组织5,6。在非人生物中，去除后仍存在比对到28S rRNA的读长，表明28S rRNA在人类与其他物种之间保守性较低6,20。当将该方法用于非人来源的样本时，可能需要补充针对特定宿主差异性rRNA序列的DNA寡核苷酸3,21。

由于该实验方案无偏向性，病毒序列读长可能仅占文库总内容的一小部分。尽管核糖体RNA（rRNA）是宿主RNA中最丰富的种类，且rRNA读长仅有一小部分被保留（<1%）在选择性去除后仍存在，所有其他宿主RNA（例如，mRNA）在耗尽后仍会残留，并可能导致样本中产生大量测序读段。因此 "过采样" (即为了获得足够的测序深度以进行病毒基因组组装和变异位点检测，需要对单个文库进行高通量测序（oversequencing）。在我们的研究中，每个样本力求测定约2000万条读长，以确保对病毒基因组及其相关变异和宏基因组内容的分析具有足够的测序深度。2,5在宏基因组学和病原体发现研究中，需要注意的是，污染性的宿主DNA会通过DNase消化被去除。因此，含有DNA基因组的病毒及其他病原体可能在此过程中丢失，但其RNA中间产物仍可能被测序。

披露

作者无竞争性财务利益。

致谢

本工作部分由美国国立卫生研究院主任办公室创新奖（编号：DP2OD06514）（PCS）的联邦资金资助，以及由美国卫生与公共服务部下属的美国国立卫生研究院过敏与传染病研究所根据合同（编号：HHSN272200900018C、HHSN272200900049C 和 U19AI110818）提供资助。

材料

本文使用的材料清单
姓名公司目录编号评论
96孔PCR板VWR47743-953
八联管盖条VWR47745-512
无核酸酶水AmbionAM993750 ml 瓶装
TURBO DNaseAmbionAM2238RNA提取后步骤，2 U/µl，含缓冲液
PCR仪任意PCR仪 
Agencourt RNAClean XP SPRI磁珠 Beckman Coulter GenomicsA63987用于RNA纯化的磁珠
实时定量PCR系统任意系统
DynaMag-96侧裙磁板Invitrogen12027
70%乙醇现配
qRT-PCR引物IDT DNA见 表2
5 M NaCl AmbionAM9760G
1 M Tris-HCl pH 7.4 SigmaT2663-1L
1 M Tris-HCl pH 7.5 Invitrogen15567-027
1 M MgCl2 AmbionAM9530G
线性聚丙烯酰胺 Ambion AM9520
覆盖完整rRNA区域的DNA寡核苷酸IDT DNA见 表3，订购实验室就绪型，浓度为100 µM
Oligo (dT)IDT DNA40 nt 长，脱盐处理
Hybridase耐热RNase H EpicentreH39100
无RNase的DNase试剂盒 Qiagen79254选择性去除步骤后使用
SUPERase-In RNase抑制剂AmbionAM2694
随机引物 Invitrogen48190-011主要为六聚体
10 mM dNTP混合液New England BiolabsN0447L
SuperScript III反转录酶 Invitrogen18080-093含第一链合成缓冲液、DTT
空气孵育仪任意空气孵育仪 
NEBNext第二链合成（无dNTP）反应缓冲液New England BiolabsB6117S10x
E. coli DNA连接酶New England BiolabsM0205L10 U/μl
E. coli DNA聚合酶I New England BiolabsM0209L10 U/μl
E. coli RNase HNew England BiolabsM0297L2 U/μl
0.5 M EDTAAmbionAM9261
Agencourt AMPure XP SPRI磁珠Beckman Coulter GenomicsA63881用于DNA纯化的磁珠
洗脱缓冲液 Qiagen10 mM Tris HCl，pH 8.5
Quant-iT dsDNA HS检测试剂盒 InvitrogenQ32854
Qubit荧光计 InvitrogenQ32857
Nextera XT DNA样本制备试剂盒 IlluminaFC-131-1096
Nextera XT DNA索引试剂盒 IlluminaFC-131-1001
Tapestation 2200AgilentG2965AA
高灵敏度D1000试剂Agilent5067-5585
高灵敏度D1000 ScreenTapeAgilent5067-5584
BioAnalyzer 2100 AgilentG2939AA
高灵敏度DNA试剂Agilent5067-4626
文库定量完整试剂盒（通用型）Kapa BiosystemsKK4824可替代Tapestation或BioAnalyzer用于文库定量

参考文献

  1. Stremlau, M. H., et al. Discovery of novel rhabdoviruses in the blood of healthy individuals from West Africa. PLoS Negl Trop Dis. 9, e0003631(2015).
  2. Gire, S. K., et al. Genomic surveillance elucidates Ebola virus origin and transmission during the 2014 outbreak. Science. 345, 1369-1372 (2014).
  3. Morlan, J. D., Qu, K., Sinicropi, D. V. Selective depletion of rRNA enables whole transcriptome profiling of archival fixed tissue. PLoS One. 7, e42882(2012).
  4. Park, D. J., et al. Ebola Virus Epidemiology, Transmission, and Evolution during Seven Months in Sierra Leone. Cell. 161, 1516-1526 (2015).
  5. Andersen, K. G., et al. Clinical Sequencing Uncovers Origins and Evolution of Lassa Virus. Cell. 162, 738-750 (2015).
  6. Matranga, C. B., et al. Enhanced methods for unbiased deep sequencing of Lassa and Ebola RNA viruses from clinical and biological samples. Genome Biol. 15, 519(2014).
  7. Tang, F., et al. RNA-Seq analysis to capture the transcriptome landscape of a single cell. Nat Protoc. 5, 516-535 (2010).
  8. Jiang, L., et al. Synthetic spike-in standards for RNA-seq experiments. Genome Res. 21, 1543-1551 (2011).
  9. Agilennt Technologies. , Available from: http://www.agilent.com/cs/library/usermanuals/Public/G2946-90004_Vespucci_UG_eBook_(NoSecPack).pdf. (2015).
  10. Kapa Biosystems. , Available from: https://www.kapabiosystems.com/product-applications/products/next-generation-sequencing-2/library-quantification/ (2015).
  11. Illumina Technologies. , Available from: https://support.illumina.com/content/dam/illumina-support/documents/documentation/system_documentation/miseq/preparing-libraries-for-sequencing-on-miseq-15039740-d.pdf. (2015).
  12. Kircher, M., Sawyer, S., Meyer, M. Double indexing overcomes inaccuracies in multiplex sequencing on the Illumina platform. Nucleic Acids Res. 40, 3(2012).
  13. Andrews, S. Babraham Bioinformatics. , Available from: http://www.bioinformatics.babraham.ac.uk/projects/fastqc/ (2010).
  14. Trombley, A. R., et al. Comprehensive panel of real-time TaqMan polymerase chain reaction assays for detection and absolute quantification of filoviruses, arenaviruses, and New World hantaviruses. Am J Trop Med Hyg. 82, 954-960 (2010).
  15. Hu, Y., et al. Serial high-resolution analysis of blood virome and host cytokines expression profile of a patient with fatal H7N9 infection by massively parallel RNA sequencing. Clin Microbiol Infect. 21, e1-4 713(2015).
  16. Simon-Loriere, E., et al. Distinct lineages of Ebola virus in Guinea during the 2014 West African epidemic. Nature. 524, 102-104 (2015).
  17. Folarin, O. A., Happi, A. N., Happi, C. T. Empowering African genomics for infectious disease control. Genome Biol. 15, 515(2014).
  18. Blainey, P. C., Quake, S. R. Digital MDA for enumeration of total nucleic acid contamination. Nucleic Acids Res. 39, 19(2011).
  19. Malboeuf, C. M., et al. Complete viral RNA genome sequencing of ultra-low copy samples by sequence-independent amplification. Nucleic Acids Res. 41, 13(2013).
  20. Gonzalez, I. L., Sylvester, J. E., Smith, T. F., Stambolian, D., Schmickel, R. D. Ribosomal RNA gene sequences and hominoid phylogeny. Mol Biol Evol. 7, 203-219 (1990).
  21. Adiconis, X., et al. Comparative analysis of RNA sequencing methods for degraded or low-input samples. Nat Methods. 10, 623-629 (2013).

重印与许可

标签

RNA病毒测序病毒基因组学cDNA合成RNase H耗竭随机引物标签化反应文库构建病毒监测
临床样本中RNA病毒的无偏深度测序
联系我们
向图书馆推荐

隐私政策

使用条款

政策

版权所有 © 2026 MyJoVE Corporation。保留所有权利。