细菌转录组分析样品制备方法
Sample Preparation for Bacterial Transcriptome Analysis   

材料与试剂

1. 50x TAE缓冲液 (Aladdi, catalog number: T197242)
2. DNA Ladder (Promega, catalog number: Z5300)
3. Agarose (中科瑞泰生物科技有限公司, catalog number: 9012-36-6)
4. RNAiso plus (TaKaRa, catalog number: T9180)
5. MICROB Express Kit (Ambion, catalog number: AM1905)
6. Glycerol (Aladdin, catalog number: G116205)
7. PrimeScriptTM 1st strand cDNA Synthesis (TaKaRa, catalog number: 6110A)
8. Synthesis (Caliper, catalog number: G2938)
9. DNA 1000 LabChip (TaKaRa, catalog number: 9151A)
10. 10x PCR Buffer (Mg2+ plus) (TaKaRa, catalog number: 4030)
11. dNTP Mixture (TaKaRa, catalog number: 2050A)
12. T4 RNA Ligase (TaKaRa, catalog number: 2050A)
13. T4 RNA Ligase Buffer (TaKaRa, catalog number: 2050A)
14. 0.1% BSA (TaKaRa, catalog number: 2050A)
15. Reverse Transcriptase XL (AMV) (TaKaRa, catalog number: 2621)
16. Glycogen (碧云天生物有限公司, catalog number: D0812)
17. AMV RT Buffer (BBI生命科学有限公司, catalog number: B610020)
18. DreamTaq DNA Polymerase (EP0704) (赛默飞世尔科技有限公司, catalog number: B300107)
19. SYTO 9 Green (Invitrogen, Eugene, Oregon, catalog number: s34854)
20. mirVana miRNA Isolation Kit (赛默飞世尔科技有限公司, catalog number: AM1561)
21. EDC粉末 (BBI生命科学有限公司, catalog number: C600433)
22. 0.1mol/L MES buffer (pH4.8) (生工生物工程股份有限公司, catalog number: C506052)
23. DEPC-Treated Water (赛默飞世尔科技有限公司, catalog number: AM9916)
24. NaOH (BBI生命科学有限公司, catalog number: A100173)
    

仪器设备

1. 生化恒温摇床器 (上海精宏实验设备有限公司, model: HZQ-QA)
2. 生化培养箱 (上海一恒科学仪器有限公司, model: SPX-250)
3. 水平电泳槽 (北京六一生物科技有限公司, model: DYCP-31DN)
4. 电泳仪 (北京六一生物科技有限公司, model: DYY-6D)
5. 微波炉 (美的集团有限公司, model: M1-L213C)
6. 冰箱 (青岛海尔股份有限公司, model: HYC-390)
7. 微量分光光度计 (Nanodrop, model: 2000c)
8. 台式冷冻离心机 (Eppendorf, model: TG16B)
9. 酶标分析仪 (Bio-Rad, model: iMark)
10. Bioanalyzer (安捷伦科技 (中国) 有限公司, model: 2100)
11. PCR核酸扩增 (Bio-Rad, model: C1000 Touch)
12. 凝胶成像分析系统 (Bio-Rad, model: GelDoc XR+)
13. 高压型电泳电源 (Bio-Rad, model: PowerPac HV)
14. 电热恒温水浴锅 (苏州江东精密仪器有限公司, model: HH-8)
15. 制冰机 (北京长流仪器有限公司, model: XD50KG)
16. 液相色谱 (安捷伦科技中国有限公司, model: 1200)
17. 紫外分光光度计 (上海元析仪器有限公司, model: JH-2100)
18. 超净工作台 (无锡一净净化设备有限公司, model: MKSD-CJ)
19. 立式压力蒸汽灭菌器 (上海博迅医疗生物仪器有限公司, model: LS-35LD )
20. 震荡仪 (海门市其林贝尔仪器制造有限公司, model: SI-0246)
21. 微型离心机 (海门市其林贝尔仪器制造有限公司, model: KST110)
22. 扫描仪 (CapitalBio, model: LuxScan 10K/A )
23. 金属浴 (卡尤迪生物科技有限公司, model: TCS10)
24. 超声仪 (上海远淮化工科技有限公司, model: KQ-3200)
25. 精密电子天平 (上海菁海仪器有限公司, model: UW820H)
26. 超低温冰箱 (青岛海尔股份有限公司, ,model: DW-86L959BP)
27. 热循环仪 (MJ Research, Watertown, MA, model: PTC-220, PTC-225)
28. 涡旋振荡器 (奥然科学技术有限公司, model: Vortex-Genie 2)
29. 超纯水系统 (Millipore, model: Milli-Q)
30. 各种量程的微量移液器 (吉尔森公司)
    

实验步骤

一、细菌培养与收集。
        在相应抗性 (目的细菌菌株所含抗性) 的培养基中培养单一的目的菌种到对数中后期，须严格控制培养条件。在4 °C预冷的离心机中10,000~12,000 × g离心3~5 min沉降菌体于50 ml离心管中。留少量菌液，将菌体重悬 (全程需无菌操作)。对于厌氧细菌，除按如上操作外，在培养与菌体收集过程中需严格控制厌氧条件。

二、细菌总RNA的提取。
        目前，常用的细菌RNA提取方法主要有变性法、CTAB法、SDS-Phenol苯酚法和热硼酸法 (Chauhan et al., 2018)。在商业化的今天，国内外生物试剂公司开发了许多细菌RNA提取试剂盒。例如，Takara的RNAiso plus、Invitrogen的Trizol总RNA提取试剂盒、PurelinkTM RNA Mini试剂盒、上海生工的UNIQ-10色谱柱、Omega的细菌RNA试剂盒、HP total RNA Kit、Qiagen公司的RNeasy Protect Bacteria Mini and Midi Kits等，能够快速有效地从细胞或培养细胞中提取高质量和高纯度的RNA。经文献参考对比发现，Takara公司的RNAiso plus总RNA提取效果相对较好 (Rodríguez et al., 2019)。该试剂盒的具体步骤如下：

1. 取新鲜或者-80 °C冻存的菌体细胞，约100 mg，加入1 ml RNAiso plus裂解，振荡混匀后室温静置5 min。
2. 以每毫升RNAiso对应0.2 ml氯仿的比例加入氯仿，振荡15 s，室温静置3 min。
   
3. 12,000 × g，4 °C，离心15 min，取上清液至新的Eppendorf管中 (上层：无色水相为RNA，中层：白色为DNA，下层：红色为蛋白)。
4. 在得到的水相中加入等体积异丙醇，混匀后-20 °C放置20~30 min (提取dsRNA时，加入等体积的异丙醇和3 M NaCl混合液，异丙醇：3 M NaCl = 1:1)。
5. 12,000 × g，4 °C，离心10 min，移除上清。
6. 加入1 ml 75%乙醇，该试剂由DEPC (焦碳酸二乙酯) 和无水乙醇配置，-20 °C预冷，洗涤沉淀。
7. 8,000 × g，4 °C，离心5 min去除上清，干燥RNA沉淀。
8. 用20~30 μl DEPC水溶解RNA。
   https://www.takarabiomed.com.cn/ProductShow.aspx?m=20141220151857153056&productID=20141226160653343219以上为RNAiso plus（Takara）试剂盒说明书下载地址。
   注意：所有的试管、移液器吸头等都要求无核酸酶级别，或经DEPC处理。
   
   

三、总RNA质量检测

1. 总量：微量分光光度计测260 nm吸收值计算。
2. 纯度：微量分光光度计测260 nm/230 nm (>1.8) 吸收值的比值，用于评估有机溶剂残留情况；260 nm/280 nm吸收值的比值 (1.8~2.1)，用于评估蛋白质污染比例情况 (Schroeder et al., 2006)。
3. 完整性。在进行bioanalyzer分析之前，可先进行琼脂糖凝胶电泳，对总RNA质量进行初步分析。电泳显示有无基因组DNA污染，有无RNA降解 (23S rRNA，16S rRNA条带)；通过目测23S rRNA和16S rRNA条带的比例可以初步评价总RNA的质量。一般认为23S rRNA:16S rRNA ≥ 2可以初步判定总RNA完整性较好 (Schroeder et al., 2006)。电泳RNA应出现两条清晰的条带 (23S/16S rRNA)，有时会有第三条条带 (5S rRNA)，最上方不可出现DNA污染条带。具体评价标准见结果示例分析。若结果较为理想即可进行下一步分析。以Agilent Bioanalyzer进行毛细管电泳（Capillary Electrophoresis），并以软件的RIN (RNA Integrity Number) 分数评估，10为RNA完整性最好，0为最差，推荐使用RIN值在8.0以上的RNA进行建库和测序 (Schroeder et al., 2006)。Agilent bioanalyzer 2100具体操作步骤简图 (图1) 如下，按Agilent 2100芯片操作SOP依次将胶、染料混合液、marker、Ladder、样品等加入对应的芯片孔内，随即放入2100 Bioanalyzer进行分析。RNA质量评价标准见结果分析示例。
   
   
   图1. RNA上样简图 (Agilent 2100说明书)
   
   

四、mRNA的富集
        mRNA只占原核细胞总RNA的1~5% (图2)，因此建议在转录组测序前对mRNA进行富集。虽然mRNA的富集不是绝对必要的，但它可以显著增加转录组覆盖率，从而提高所得到的转录组图谱的分辨率 (Kumar et al., 2016)。未经富集的低覆盖率样品可以通过测序更多的cDNA来弥补。但是从经济角度出发，这可能会大大增加其实验成本。目前，已知的细菌mRNA的分离纯化方法主要有以下几种：rRNA消减杂交、5’单核苷酸依赖的外切酶处理法、选择性引物扩增法、依赖于双链特异核酸酶的cDNA均一化法、大肠杆菌 poly(A)聚合酶加尾法、与RNA结合蛋白Hfq免疫共沉淀法 (Kumar et al., 2016)。此处将对rRNA消减杂交法进行详细说明，其余方法详细步骤及原理可在附件查看。

        

        图2. 细胞内RNA分布图

        rRNA消减杂交法。目前应用此原理的试剂盒主要有MICROB Express Ki (Ambion)、RiboMinus bacteria transcriptome isolation Kit (Invitrogen)、Ribo-Zero rRNA removal Kit (Epicentre) 等。本文以MICROB Express Kit (Ambion) 试剂盒为例进行说明，该试剂盒适用于绝大多数细菌种属 (古菌、支原体等不可用)。其原理为：16S和23S rRNA与总RNA样本进行杂交，使用与rRNA互补的寡核苷酸探针进行杂交后，带有16S和23S rRNA的磁珠在磁力的作用下，移至管子的一侧 (Stewar et al., 2010)。上清液中浓缩的RNA被转移至新的试管中。洗涤磁珠，以保证磁珠上没有RNA的残留。最后用乙醇沉淀RNA。此过程产生的RNA包含mRNA、tRNA、5S rRNA和其他小RNA (图3)。实验具体操作步骤如下：

1. 将2~10 μg总RNA添加到200 μl结合缓冲液中。
2. 加入4 μl Capture Oligo Mix。
3. 加热至70 °C，10 min。
4. 37 °C，持续15 min。
5. 每个样品取出50 μl Oligo MagBeads到1.5 ml管中。
6. 使用磁力试管架捕获Oligo MagBeads，并小心地移除和丢弃上清液。
7. 用等量的无核酸酶水清洗Oligo MagBeads。
8. 通过与相同体积的缓冲液结合，以平衡Oligo MagBeads。
9. 重新悬浮Oligo MagBeads至相同体积的缓冲液中，并使缓冲液达到37 °C。
10. 同时将洗涤液预热至37 °C。
11. 在RNA预制液中，加入50 μl制备的Oligo MagBeads，混匀，在37 °C下孵育15 min。
12. 捕获Oligo MagBeads，并将富集的mRNA上清液移至收集管。
13. 用100 μl洗涤液在37 °C下洗涤，并回收洗涤液，从而回收Oligo MagBeads中残留的mRNA，与前一步骤的上清液合并。
14. 乙醇沉淀富集的mRNA。
15. 将富集的mRNA重新悬浮在适当的缓冲液中。
    注意：MICROB Express Kit (Ambion) 试剂盒详细说明书可以在以下网址下载。https://www.thermofisher.com/order/catalog/product/AM1905#/AM1905 。
    
    
    图3. rRNA的杂交捕获 (MICROB Express Kit说明书)
    
    

五、纯化后mRNA质量的测定。
        方法同上文总RNA质量检测。

六、片段化处理。
        mRNA纯化之后的文库构建通常有两种思路，一种是首先用oligo (dT)引物反转录mRNA，再进行cDNA的片段化；另外一种则是先将mRNA打断，再结合随机引物进行反转录。先针对mRNA进行打断再进行反转录获得测序序列主要是针对基因本体。而先转录再进行片段化的方法，尤其是结合oligo (dT)进行反转录获得的测序结果对转录本3'端具有比较强的偏好性。所以在mRNA-Seq中建议采用先对mRNA打断再进行反转录文库构建的方法。以上两种方法具体步骤如下：

1. mRNA片段化。
           建议使用专门试剂Fragmentation Reagent (BioVendor) 将纯化的mRNA片段化。mRNA片段化完成后，再结合随机引物进行反转录。 mRNA片段化处理方法主要包括碱处理法、金属离子 (Mg²⁺、Zn²⁺) 溶液处理法、酶 (RNaseIII) 处理法等 (Kumar et al., 2012)。mRNA片段化之后应立即进行第一条cDNA链的合成（详细步骤见“七. cDNA文库的构建”），因为mRNA在该体系下非常容易降解。选择金属离子 (Mg²⁺、Zn²⁺ ) 溶液处理法打断时需根据需要的文库大小选择合适的片段化温度和时间 (常用设置条件：150~200 bp，94 °C，15 min； 200~300 bp，94 °C，10 min；250~550 bp，94 °C，5 min)。本文以Hieff NGS® MaxUpTM II Dual-mode mRNA Library Prep Kit for Illumina®为例进行详细说明。具体操作步骤如下。试剂盒详细说明书可在此地址下载https://www.bio-equip.com/show1equip.asp?equipid=4505343。
   
   1.1

   将mRNA Capture Beads从2~8 °C取出，静置使其温度平衡至室温，约30 min。
   1.2

   准备一个Nuclease free离心管，取0.1~4 μg总RNA，用Nuclease free水将体积补至50 μl，冰上放置备用。
   1.3

   颠倒或旋涡振荡混匀磁珠，吸取50 μl磁珠悬液加入至50 μl总RNA样品中，用移液器吹打6次，使其充分混匀。
   1.4

   将磁珠与RNA的混合物置于PCR仪中，65 °C，5 min；4 °C，静置，使得RNA变性。
   1.5

   室温孵育5 min，使mRNA与磁珠完全结合。
   1.6

   将样品置于磁力架中，室温静置5 min，使mRNA与总RNA分离，小心移除上清。
   1.7

   将样品从磁力架上取出，用200 μl Beads Wash Buffer重悬磁珠，移液器反复吹打6次以彻底混匀。将样品置于磁力架中，室温静置5 min，小心移除上清。
   1.8

   重复步骤1.7，共洗涤两次。
   1.9

   将样品从磁力架上取出，加入50 μl Tris Buffer重悬磁珠，用移液器反复吹打6次以彻底混匀。
   1.10

   将样品置于PCR仪中，80 °C，2 min；25 °C，静置，将mRNA洗脱下来。
   1.11

   将样品从PCR仪中取出，加入50 μl Beads Binding Buffer，用移液器反复吹打6次以彻底混匀。
   1.12

   室温放置5 min，使mRNA结合到磁珠上。
   1.13

   将样品置于磁力架中，室温静置5 min，小心移除上清。
   1.14

   将样品从磁力架上取出，用200 μl Beads Wash Buffer重悬磁珠，移液器反复吹打6次以彻底混匀，将样品重新放回至磁力架中，室温静置5 min，吸掉全部上清。
2. cDNA的片段化及文库构建。
           先用随机引物反转录mRNA，再进行cDNA的片段化。首先得到长段的一链cDNA，之后反转录为双链cDNA（详细步骤见“七. cDNA文库的构建”），再利用DNA片段化的方法对其进行片段化处理 (酶切法或机械法)，随后纯化得到cDNA文库。酶切法构建DNA文库，相较于传统的超声法和转座酶法具有以下优势。首先，无偏好片段化酶：具有稳定的片段化效果，无需复杂的机械片段化过程。酶切产物片段大小同样本类型、Input、DNA量和GC含量等均无相关性，仅与酶切时间有关。其次，建库流程精简。一步完成片段化、末端修复、加A、纯化等。末端修复加A (30 min)-接头连接 (15 min)-纯化、分选 (30 min)-文库扩增 (15 min)-纯化分选 (30 min)，常规建库预计耗时约2 h， PCR-free建库预计耗时1 h。最后，宽泛的样本投入范围。相比较固定投入量，固定酶切时间的TN5建库，酶切法建库适用500~1 μg，满足个性化建库。本文以Hieff NGS® OnePot DNA Library Prep Kit for Illumina®为例，对酶切法制备cDNA文库进行详细说明。试剂盒说明书下载地址为https://www.yeasen.com/products/detail/932。该试剂盒具体步骤如下：
   
   2.1

   DNA 片段化、末端修复、dA 尾添加 (DNA Fragment/End Preparation/dA-Tailing)该步骤将基因组DNA片段化，同时进行末端修复及 dA 尾添加。
   

   1)

   将Input DNA 、Smearase® Mix 、10 ddH2O解冻后，颠倒混匀，置于冰上备用。

   2)

   于冰上配制反应体系。加入DNA X μl、Smearase® Mix 10 μl，补水ddH2O至60 μl。

   3)

   使用移液器轻轻吹打或低速振荡混匀，并短暂离心将反应液离心至管底。

   4)

   将上述 PCR 管置于 PCR 仪，反应程序为热盖，105 °C ； 4 °C ，1 min； 30 °C ，3~20 min； 72 °C ，20 min ；4 °C ，静置。进行DNA片段化，末端修复及 dA 尾添加反应。
   2.2

   接头连接 (Adapter Ligation) 该步骤将2.1步骤的产物末端，连接特定的Illumina®接头。
   

   1)

   根据Input DNA量按试剂盒说明稀释Adapter至合适浓度。

   2)

   将dA-tailed DNA、Ligation Enhancer、Fast T4 DNA Ligase、 DNA Adapter解冻后颠倒混匀，置于冰上备用。

   3)

   于PCR管中配制表反应体系。具体为：dA-tailed DNA， 60 μl； Ligation Enhancer ，30 μl；Fast T4 DNA Ligase ，5 μl； DNA Adapter 5 μl。

   4)

   使用移液器轻轻吹打或振荡混匀，并短暂离心将反应液收集至管底。

   5)

   将 PCR 管置于 PCR 仪中，反应程序为热盖，105 °C； 20 °C， 15 min； 4 °C，静置。进行接头连接反应。
   2.3

   连接产物磁珠纯化 (Post Ligation Clean Up) 该步骤使用磁珠对上一步骤的产物进行纯化或分选。纯化可除去未连接的 Adapter 或 Adapter Dimer 等无效产物。
   

   1)

   将 Hieff NGS® DNA Selection Beads 磁珠由冰箱中取出，室温平衡至少 30 min。配制 80%乙醇。

   2)

   涡旋振荡或充分颠倒磁珠以保证充分混匀。

   3)

   吸取 60 μl Hieff NGS® DNA Selection Beads (0.6x, Beads:DNA=0.6:1) 至 Adapter Ligation 产物中，室温孵育 5 min。

   4)

   将 PCR 管短暂离心并置于磁力架中分离磁珠和液体，待溶液澄清后 (约 5 min)，小心移除上清。

   5)

   保持 PCR 管始终置于磁力架中，加入 200 μl新鲜配制的 80%乙醇漂洗磁珠，室温孵育 30 s 后，小心移除上清。

   6)

   重复步骤 5)，总计漂洗两次。

   7)

   保持 PCR 管始终置于磁力架中，开盖空气干燥磁珠至刚刚出现龟裂 (不超过 5 min)。

   8)

   将 PCR 管从磁力架中取出，进行洗脱。
   2.4

   文库扩增 (Library Amplification) 该步骤将对纯化或长度分选后的接头连接产物进行 PCR 扩增富集。
   

   1)

   将2x Super Canace® II High-Fidelity Mix、Primer Mix、Adapter Ligated DNA (上一步骤产物) 解冻后颠倒混匀，置于冰上备用。

   2)

   于无菌 PCR 管中配制，其反应体系为2x Super Canace® II High-Fidelity Mix，25 μl；Primer Mix，5 μl；Adapter Ligated DNA，20 μl。

   3)

   使用移液器轻轻吹打或振荡混匀，并短暂离心将反应液收集至管底。

   4)

   将 PCR 管置于PCR仪中，反应程序为98 °C ，1 min ； 98 °C， 10 s；60 °C ，30 s； 72 °C ，30 s；72 °C， 5 min；4 °C，静置。进行 PCR 扩增。
   2.5

   扩增产物磁珠纯化。同上步骤中纯化操作步骤。使用 Hieff NGS® DNA Selection Beads (0.9x，Beads:DNA=0.9:1) 纯化文库扩增产物。
   2.6

   文库质量控制。详细步骤见“八. cDNA文库质量检测”。
   
   

七、cDNA文库的构建。
        为了获得覆盖整个转录组的序列数据，必须从整个RNA样本中随机产生小cDNA分子。这通常是通过随机六聚体启动的反转录来实现的。由于常规cDNA文库构建方法需进行两条链的合成，所以其正向和反向转录组信息相同。这就使得该方法特异性较差。而与常规cDNA文库构建方法相比，仅使用第一条链来构建序列库，就可以保持关于转录方向的信息。或者，使用RNA片段化、cDNA片段化方法，通过定向的、特异的、逐步连接的接头序列来保留转录模板链上的信息。此处主要对常规cDNA文库构建方法 (图4)进行详细说明，其他方法详见附件。

1. cDNA第一条链合成。在随机六聚体引物 (Random Hexamer Primer)、逆转录酶的作用下，以mRNA为模版合成一链cDNA。具体过程按照PrimeScriptTM 1st strand cDNA Synthesis Kit说明书进行PCR反应。
2. cDNA第二链合成并删除mRNA，产生双链cDNA。在cDNA第一条链混合物中加入51 μl超纯水、20 μl 5x SecondStrand Buffer 和3 μl dNTP mix (10 mM)，dNTP mix中用dUTP代替dTTP，使cDNA第二链中仅包含A/U/T/G。在冰上放置5 min，然后加入1 μl RNase H (2U/μl) 和5 μl DNA Pol聚合酶 (10 U/μl) 混匀，然后置于16 °C，2.5 h。
3. 双链cDNA的末端修复。取纯化DNA溶液35 μl，加入50 μl T4 DNA连接酶缓冲液 (2x)、4 μl dNTPs mix、T4 DNA聚合酶、T4 Polynucleotide Kinase、Klenow DNA聚合酶，在20 °C下反应30 min，进行末端修复。之后用QIAquick-PCR纯化试剂盒对cDNA进行纯化。
4. 连接poly(A)端。将前一步纯化的DNA中加入Klenow缓冲液、dATP和Klenow exo，在37 °C下反应30 min，用PCR纯化试剂盒纯化DNA，最后将样品溶解在Ethidium Bromide溶液中。每个连接序列都有一个Index序列 (6bp)，可用于不同的库构建。
5. 添加接头序列。将纯化后的DNA加入T4-DNA连接酶缓冲液、接头序列寡聚ligation solution MasterMix和T4 DNA连接酶。然后用PCR纯化试剂盒进行纯化。
6. 回收并纯化产物凝胶进行电泳。将250~300 bp (或自己实验所需的其他大小的目的片段) 的片段回收并放入Eppendorf管中。按QIAquick凝胶提取试剂盒说明书对凝胶进行纯化回收。
7. PCR扩增前，使用UNG酶消化第二链cDNA ，文库中只含有第一条。用PCR方法扩增富集cDNA文库。
8. 使用凝胶回收试剂盒对PCR产物进行回收纯化，待用。具体步骤参照试剂盒说明书。
   
           
   

                图4. cDNA文库构建方法简图

八、 cDNA文库质量检测。
        文库容量、重组率及插入片段的大小是鉴定cDNA文库质量的重要指标 (Croucher et al., 2010)。先对第二轮PCR产物进行琼脂糖凝胶电泳初步分析，再利用Agligent 2100 bioanalyzer仪进行检测，检测流程参照其说明书。Bioanalyzer分析按照DNA 1000 LabChip试剂盒所提供的说明进行。每1试剂盒含有25块芯片和以下试剂： DNA长度测定梯度标准品 (sizing ladder)、注射器、凝胶基质、染料浓缩液、DNA相对分子质量标准 (marker) 以及旋转过滤器。凝胶染料混合物的准备：将25 μl DNA染料与凝胶基质混合后，6,000 × g离心，在芯片指定位置分别加入9 μl凝胶染料混合物以及5 μl maker和1 μl ladder (1,000)，取7个PCR产物各1 μl分别加入1~7号样品孔中，最后将芯片涡旋混匀后放入仪器中进行自动检测分析。

九、若所测cDNA文库质量较好，满足建库要求，则可以进行cDNA上机测序。

结果与分析

        RNA样本总量及完整性是评判样品质量的关键点，其中RNA完整性评估依靠RIN值、23S/16S (原核生物) 以及bioanalyzer检测峰图基线是否平整来综合评判。对于降解样品，实验难以获取完整的转录本信息，这样就会影响数据质量及其完整性。当RNA总量较低时，会导致建库成功率低，或数据Duplication rate高等问题[9]。基于此，本文将展示部分较好的结果示例就如何对RNA、cDNA质量进行评价作进一步的详细说明。

1. RNA质量评价标准——凝胶电泳图。
           RNA质量评价首先要通过凝胶电泳图来对其进行初步分析。主要需查看以下几方面：目的条带是否明亮清晰；泳道内是否有明显的降解弥散区；胶孔处有无蛋白污染；有无DNA污染；有无外源物种核酸污染等 (Han et al., 2017)。由下文结果示例可以看出，该总RNA的凝胶电泳图目的条带明亮清晰，泳道无弥散区，无蛋白和DNA的污染。由此得出结论该样品质量较好，如图5所示 (图5为bioanalyzer获得)。
   
   
   图5. RNA 质量较好的示例 (Han et al., 2017)
   
   
2. RNA质量评价标准——电泳图及RIN值。
           在凝胶电泳检测完成后，需要对其进行bioanalyzer分析。图6为bioanalyzer (Bio-Rad)原核生物典型电泳图。最左侧峰为Ladder峰，中间峰为16S rRNA，最右侧峰为23S rRNA。整体样品峰图基线较为平整，峰型较好，无拖尾峰、前沿峰、包裹峰及多余的降解峰。通过mRNA富集前后对比电泳图，可以清晰地看到mRNA富集程度 (即16S rRNA、23S rRNA去除程度)。
   
   
   图6. mRNA富集前后对比 (Han et al., 2017)
   
           若样本采用Agilent bioanalyzer 2100质检。则需判别其RIN值，RIN值即为RNA完整值（RNA integrity number），上样完成后，质检报告中会显示其具体数值。RIN数值范围为1~10，将生物总RNA质量进行分类，其中1代表降解最严重的情况，10代表最完整，建议采用RIN数值8以上的样品进行cDNA文库的构建。通过这种检测方法，有助于解释电泳图，有利于样品之间的相互比较，同时确保实验的可重复性。
3. cDNA质量评价标准——凝胶电泳图。
           通过PCR进行第二链合成后，直接取5 μl PCR产物进行琼脂糖凝胶电泳检测。示例较好的结果如图7所示，双链cDNA片段弥散分布，大小主要集中在200~500 bp (实验所需cDNA片段大小，片段化条件控制文库DNA分子的长度分布具体方法可参考所用试剂盒，比如TruSeq® RNA Sample Preparation v2 Guide试剂盒)，说明实验所需mRNA均成功反转录并得到了有效扩增，其cDNA合成质量较高，基本包含全部cDNA，可满足文库构建的需要，可用于后续筛选目的基因、大规模测序等研究。
   
   
   图7. cDNA模拟胶图 (Han et al., 2017)
   
   
4. cDNA质量评价标准——电泳图。
           cDNA质量同样需要进行bioanalyzer来进一步确定其质量。图8为bioanalyzer (Bio-Rad) 原核生物cDNA典型电泳图。最左侧峰为1,700 bp内标物，中间峰为cDNA，最右侧峰为50 bp内标物（内标物可根据实验所需进行设置，包括数量、大小等）。由图可知，该样品整体峰图基线较为平整，峰型较为平滑，没有出现拖尾峰、前沿峰、包裹峰及其他的杂乱峰。cDNA峰较为明显且峰面积较大。说明文库质量较高，可满足文库构建的需要。
   
   
   图8. cDNA电泳图 (Han et al., 2017)
   

致谢

        感谢中国农业科学院科技创新工程对本研究的资助。使用本实验方案已发表的文章有：Han, D., Link, H. and Liesack, W. (2017). Response of Methylocystis sp. strain SC2 to salt stress: physiology, global transcriptome, and amino acid profiles. Appl Environ Microbiol 83(20): e00866-17.

参考文献

1. Chauhan, A., Sharma, J. N., Modgil, M. and Siddappa, S. (2018). Comparison of various RNA extraction methods， cDNA preparation and isolation of calmodulin gene from a highly melanized isolate of apple leaf blotch fungus Marssonina coronaria. J Microbiol Methods 151: 7-15.
2. Rodríguez, A. and Vaneechoutt, M. (2019). Comparison of the efficiency of different cell lysis methods and different commercial methods for RNA extraction from Candida albicans stored in RNAlater. BMC Microbiol 19(1): 94.
3. Schroeder, A., Mueller, O., Stocker, S., Salowsky R., Leiber M., Gassmann M., Lightfoot S., Menzel W., Granzow M. and Ragg, T. (2006) The RIN: an RNA integrity number for assigning integrity values to RNA measurements. BMC Molecular Biol 7: 3.
4. Kumar, N., Lin, M., Zhao, X., Ott, S., Santana-Cruz, I., Daugherty, S. and Dunning Hotopp, J. C. (2016). Efficient enrichment of bacterial mRNA from host-bacteria total RNA samples.Sci Rep 6: 34850.
5. Kumar, N., Lin, M., Zhao, X., Ott, S., Santana-Cruz, I., Daugherty, S. Rikihisa, Y., Sadzewicz, L., Tallon, L. J., Fraser, C. M. and Dunning Hotopp, J. C. (2016). Efficient enrichment of bacterial mRNA from host-bacteria total RNA samples.Sci Rep 6: 34850.
6. Stewart, F. J., Ottesen, E. A. and DeLong, E. F. (2010). Development and quantitative analyses of a universal rRNA-subtraction protocol for microbial metatranscriptomics.ISME J 4(7): 896-907.
7. Kumar, R., Ichihashi, Y., Kimura, S., Chitwood, D. H., Headland, L. R., Peng, J. and Sinha, N. R. (2012). A high-throughput method forillumina RNA-Seq library preparation. Front Plant Sci 3: 202.
8. Croucher, N. J. and Thomson, N. R. (2010). Studying bacterial transcriptomes using RNA-Seq. Curr Opin Microbiol 13(5): 619-624.
9. Hayden, C. A., Wheeler, T. J. and Jorgensen, R. A. (2005). Evaluating and improving cDNA sequence quality with cQC. Bioinformatics 21(24): 4414-4415.
10. Han, D., Link, H. and Liesack, W. (2017). Response of Methylocystis sp. strain SC2 to salt stress: physiology global transcriptome and amino acid profiles. Appl Environ Microbiol 83(20): e00866-17.