本文提供了一种通过实时PCR定量沿海沉积物中细菌16S rRNA基因及其转录本的实验方案。该方法包括DNA和RNA的共提取、无DNA污染的RNA制备,以及通过反转录定量PCR(RT-q-PCR)对16S rRNA基因和转录本进行定量,同时包含标准曲线的构建。
本文提供了一种通过实时PCR定量沿海沉积物中细菌16S rRNA基因及其转录本的实验方案。该方法包括DNA和RNA的共提取、无DNA污染的RNA制备,以及通过反转录定量PCR(RT-q-PCR)对16S rRNA基因和转录本进行定量,同时包含标准曲线的构建。
实时聚合酶链式反应(Real Time Polymerase Chain Reaction),也称为定量PCR(q-PCR),是微生物生态学中广泛使用的一种工具,用于定量环境样品中分类群和功能群的基因丰度。当与逆转录反应联用时(RT-q-PCR),还可用于定量基因转录本。q-PCR利用高度灵敏的荧光检测化学方法,可在PCR反应的指数扩增阶段对PCR扩增子进行定量,从而避免了传统“终点法”PCR在平台期检测所带来的一系列偏差。本文提供了一种通过实时PCR定量沿海沉积物中细菌16S rRNA基因及其转录本的实验方案。首先,概述了从沿海沉积物中联合提取DNA和RNA的方法,包括制备无DNA污染RNA所需的额外步骤;其次,详细描述了利用q-PCR和RT-q-PCR对提取的核酸中的16S rRNA基因及转录本进行定量的分步操作流程,其中包括DNA和RNA标准曲线构建的具体细节。此外,还列出了在微生物生态学研究中应用RT-q-PCR检测时需要考虑的关键因素。
微生物是驱动生态系统功能的生物圈基石。大多数微生物尚未被培养成功1。因此,基于分子的方法对于深入理解环境中微生物的多样性和功能至关重要。这些方法的核心是从环境样品中提取核酸,并利用聚合酶链式反应(PCR)扩增目标基因。
DNA/RNA 提取的第一步旨在裂解存在的微生物群落的细胞壁,去除不需要的非核酸分子(例如,有机和无机物质),并将 DNA/RNA 保留在溶液中,以便进行后续的下游分析。在文献中报道的多种方法2,3,4,5中,包括一系列商用提取试剂盒,Griffiths 方法6被广泛采用7,8,9。该方法成本较低,特别适用于沉积物样品,因其采用珠磨法裂解细胞,并包含减少共提取 PCR 抑制剂(如腐殖酸)的步骤,同时能够同步回收 DNA 和 RNA。
第二步采用聚合酶链式反应(PCR)从提取的核酸中扩增目标基因,例如16S rRNA分类标记基因。该方法过去并持续推动着对尚未明确的微生物“黑箱”的探索10,11。然而,基于终点PCR的方法存在多种局限性,可能对微生物群落的表征造成偏差12。为了准确定量基因/转录本的丰度,必须使用实时PCR,即定量PCR(qPCR)。qPCR利用荧光报告染料系统,在每轮PCR循环后追踪扩增子的累积情况。这一点至关重要,因为这意味着定量可在PCR反应的早期指数阶段进行,而非终点阶段,此时扩增子产量仍与目标基因的初始丰度成正比。
两种报告系统常被使用:一种是嵌入式核酸染料13 以及 5' 3' DNA聚合酶的外切酶活性14由于前者报告系统会非特异性地结合所有双链DNA,因此如果出现非特异性扩增产物或引物二聚体等副产物,可能导致对目标序列的高估。为避免此问题,可能需要对扩增条件进行大量优化。在 后者 系统,通过结合5'核酸酶活性来追踪模板扩增过程 Taq 一种能够从内部探针上切割荧光团的聚合酶。由于采用了仅与引物对之间的互补靶标特异性序列结合的荧光生成探针,该特性提高了检测的特异性。在两种化学体系中,均通过确定交叉点(Cp)实现定量,即PCR扩增产物的累积量(通过荧光信号增加来测量)显著高于背景荧光水平时的循环数。
qPCR 已广泛应用于微生物生态学,以测定不同环境中基因的丰度15。此外,通过将 RNA 反转录为 cDNA 并结合 qPCR 技术,即 RT-qPCR,可实现对基因表达的定量。因此,qPCR 和 RT-qPCR 是对环境样品中基因和/或转录本数量进行快速、有效定量的重要方法。
近海沉积物中的微生物驱动着多种生态系统过程,包括有机物的矿化、污染物的降解以及氮等大量营养元素的生物地球化学循环16,17,18。要全面理解这些转化过程,必须系统掌握参与其中的微生物种群信息,包括基因和转录本丰度的定量数据。本文介绍了一系列经过充分验证、流程优化且标准化的操作方案,用于定量近海沉积物中细菌16S rRNA基因及其转录本的丰度。该方案涵盖了样品采集、DNA与RNA的同步提取、无DNA的RNA制备、提取核酸的质量检测、16S rRNA-DNA和-RNA标准品的制备,以及环境样品的定量分析。通过本方法获得的定量数据,对于揭示驱动近海生态系统功能的微生物群落具有重要意义。
1. DNA & 从海洋近海沉积物中提取 RNA
2. RNA 的制备及无 DNA RNA 的质量检测
3. 从RNA生成第一条cDNA链
4. 定量 PCR

从沉积物中提取高质量的DNA和RNA是定量基因和转录本丰度的第一步。成功的提取结果应显示出清晰的DNA和RNA条带,如图1中样品A-C所示,除高分子量的基因组DNA条带外,还可观察到清晰的23S和16S rRNA条带。

图1. DNA/RNA 提取。 从三份平行样品(A-B-C)各0.5 g海岸沉积物中提取DNA/RNA的典型结果。取5 µl DNA/RNA样品在1.4%琼脂糖凝胶中以85 V电泳40分钟。使用了分子量范围为10,037-200 bp的分子标记物。 请点击此处查看该图的放大版本。
制备RNA时,必须对共提取的DNA进行消化处理。随后需对RNA进行16S rRNA终点PCR检测,以确认DNA已完全去除。若DNA已被彻底清除,则仅在阳性对照中观察到条带。为确保抑制物不会影响PCR产物的生成,必须同时使用未稀释和1:10稀释的RNA进行检测。RNA现在可进行逆转录反应,将其转化为cDNA。该反应可使用基因特异性引物或随机六聚体引物进行。通常在PCR仪中进行反应,以确保最佳的温度程序。所得cDNA将作为后续qPCR反应的模板。通过定量DNA和RNA来确定每次反应所用的模板浓度。
针对16S rRNA DNA的qPCR反应方案已列出。对于16S rRNA转录本,需使用cDNA替代标准曲线和模板。由于未知样品的定量依赖于标准曲线,因此必须确保标准曲线具有良好的质量。 图2 显示了(A)标准品及环境DNA稀释液的制备,(B)标准曲线与环境样品的扩增,(C)标准曲线的线性回归转换及基因拷贝数的计算。当10倍梯度稀释的标准品被正确制备并扩增时,各标准稀释点之间会出现3.3个循环的差异(在100%扩增效率下,模板量每增加10倍所需循环数为3.3个)图2Bi)。标准曲线的参数,包括 R2 0.99 的值,且 PCR 扩增效率在 90% 至 110% 范围内(图2C)时尤为需要。如果无模板对照(NTC)出现扩增,必须报告其Cp值。此时,标准品和未知样品的Cp值截断阈值应设定为比NTC低3.3个循环即,设定为高于NTC的Cp值一个对数倍数(log fold)20 (图 2Cii)。通过沉积物DNA提取的未知模板,从原液及10倍稀释样本中进行定量-1, 10-2 和 10-3 稀释样品(B)。未稀释样品未扩增,10倍稀释样品的Cp值为10-1 至 10-3 稀释倍数分别为 24.12、26.02 和 28.40。NTC 的 Cp 值为 30.5,设定 NCT 阈值为 27.2。将基因丰度转换为 g-1 湿重沉积物产生 2.5 x 107 和 7.1 x 107 10分钟-1 和 10-2 分别为稀释范围。10-3 稀释液Cp低于NTC阈值,因此未被采用。在此情况下,10-2 被选为最佳模板稀释度。

图2. 来自海岸沉积物的标准品和环境DNA的16S rRNA基因qPCR扩增。(A)i)DNA标准曲线的制备和ii)用于qPCR扩增的环境DNA稀释液的制备;(B)i)DNA标准曲线和ii)环境样品的qPCR扩增,显示每个样品的Cp值;(C)i)标准曲线的线性回归及其标准曲线参数;ii)由环境样品计算基因丰度。NTC:阴性模板对照。请点击此处查看该图的放大版本
| 试剂 | 用量 |
| CTAB | 5 g |
| K2HPO4x3H2O | 2.58 g |
| KH2PO4 | 0.1 g |
| NaCl | 2.05 g |
| dH2O | 加至 100 ml |
表1. CTAB-磷酸盐缓冲液的配制。
| 试剂 | 用量 |
| PEG 6000 | 30 g |
| NaCl | 9.35 g |
| dH2O | 加至 100 ml |
| 注意:在中等加热并搅拌 50 ml dH2O 的条件下,缓慢加入 PEG 6000。 | |
表2. PEG-NaCl 沉淀溶液的配制。
| 试剂 | 用量 |
| 未稀释 RNA / 1:10 稀释 RNA | 1.0 µl |
| 10x PCR 缓冲液 | 5.0 µl |
| dNTPs 10 mM | 1.0 µl |
| 63F 正向引物 | 1.0 µl |
| CAG GCC TAA CAC ATG CAA GTC | |
| 518R 反向引物 | 1.0 µl |
| ATT ACC GCG GCT GCT GG | |
| Taq 聚合酶 | 0.4 µl |
| 无 RNA 酶的 dH2O | 40.6 µl |
表3. 无DNA的RNA质量检测PCR预混液。
| 试剂 | 用量 |
| RNA | 0.5-5 μg |
| dNTPs 10 mM | 1 µl |
| 逆转录引物 10 μM / 随机六聚体 50 µM | 1 µl |
| DEPC-H2O | 加至 13 µl |
表4. cDNA合成反应混合液A。
| 试剂 | 用量 |
| 5x 缓冲液 | 4 µl |
| DTT | 1 µl |
| RNA酶抑制剂 | 1 µl |
| 逆转录酶 | 1 µl |
表5. cDNA合成反应体系B。
| 靶标 | 引物名称 | 序列 | 退火温度 (°C) | 参考文献 |
| 16S rRNA 细菌 | Bact1369F | CGG TGA ATA CGT TCY CGG | 56 | Suzuki et al. 2000 19 |
| Prok1492R | GGW TAC CTT GTT ACG ACT T | |||
| TM1389F (探针) | CTT GTA CAC ACC GCC CGT C |
表6. 16S rRNA 定量PCR检测所用的引物和探针。
| 试剂 | 用量 |
| 2x Master Mix | 10 μl |
| 正向引物 (10 μM) | 0.8 μl |
| 反向引物 (10 μM) | 0.8 μl |
| 探针 (10 μM) | 0.4 μl |
| 水 | 7.0 μl |
| 终体积 | 19 μl |
表7. 实时PCR反应体系。
| 描述符 | 意义 |
| 相关系数 (R2) | 标准曲线线性程度的度量。理想情况下应为 R2=1,R2=0.98–0.99 的值可接受。 |
| 斜率 (α) | 反应效率的度量。理想情况下应等于 -3.32,可接受范围为 -3.58 至 -3.1。 |
| 扩增效率 (=10(-1/斜率)-1) | 若为 100%,则在指数扩增期间每个热循环后模板数量翻倍。良好的扩增效率范围为 90% 至 110%。 |
| Y 轴截距 (β) | 反应的理论检测极限。尽管不用于量化反应灵敏度,但在比较同一靶标的不同的标准曲线时具有重要意义。 |
表8. qPCR反应参数描述。
DNA/RNA 提取与qPCR相结合,为从多种环境样品(如近海沉积物)中灵敏定量基因和转录本丰度提供了一种快速、准确且相对经济的方法。
初始核酸提取是确保获得具有代表性的微生物群落组成的关键步骤22. 提取方案存在若干局限性需加以考虑:a)细胞完全裂解的实现,b)抑制性化合物的共提取例如,腐殖酸或多酚类物质),c)RNA组分中的DNA污染,d)提取的核酸在储存过程中迅速降解。必须采取预防措施以避免这些局限性。例如,需特别注意确保核酸提取方法针对特定样本类型进行了优化例如沉积物、土壤、废水 等等。)。通过开展预实验,可显著提高DNA和RNA的核酸得率和质量23. 为尽量减少抑制性化合物对基于 PCR 的应用的影响24,对提取的 DNA 和 RNA 进行一系列稀释度测试。为避免因多次冻融循环导致提取的 DNA/RNA 快速降解,并防止可能的遗传信息丢失,应将多个小体积等分试样于 -80 °C 保存。
当精心设计时,qPCR 是一种稳健、高度可重复且灵敏的方法。值得注意的是,本方案中所述的扩增及标准曲线构建方法可适用于任何目标基因,包括其他系统发育标记即,古菌16S rRNA,真菌18S rRNA)或环境中重要功能相关的基因。qPCR技术已知的局限性包括:a)难以生成可重复的高质量标准曲线以实现绝对定量;b)引物/探针的选择及qPCR检测条件的优化;c)使用低质量或剪切的核酸;d)DNA/RNA工作稀释度的选择不当可能导致抑制效应。此外,需注意qPCR技术提供的是基因或转录本的丰度数据,未必等同于实际的细胞数量:当靶标为16S和18S rRNA基因时尤其如此,因为不同微生物在其基因组中核糖体RNA基因的拷贝数存在差异25.
标准曲线质量不佳将导致目标基因的定量结果不准确。建议将高浓度标准品分装成小份储存,以便每次使用时制备新鲜的标准曲线。切勿储存已稀释的标准曲线,每次均应从最高浓度的储备液中新鲜配制稀释系列。为准确测定环境样品,标准曲线的浓度范围应覆盖未知样品预期的Cp值范围。在定量转录本时,应使用RNA而非双链DNA构建标准曲线。在可能的情况下,为避免不同批次间的变异,需直接比较的样品应于同一次实验中完成定量。20这并不总是可行的。因此,为了比较不同检测之间获得的基因丰度,建议在不同检测之间对部分随机选取的样本进行重复检测。目前已有多种针对微生物类群和功能类群的qPCR引物和探针组合可供选择。15. 在选择这些参数时需仔细考量,以确保对目标群体具有最大程度的覆盖性和特异性。如果qPCR检测的反应效率不理想,应首先通过测试不同的热循环条件(如退火时间和/或温度)和/或反应条件来进行故障排查。 例如,调整引物-探针浓度。在q-PCR检测方法及反应条件优化后,始终需对一系列DNA/cDNA稀释浓度进行初步测试,以确定合适的模板浓度。选择产生最高拷贝数的稀释范围作为后续检测的最优模板稀释度。
目前,新一代测序技术能够高效揭示多种环境中微生物群落的结构与功能26,27,28。然而,这些数据集通常基于终点PCR扩增子文库,因此仅能对特定分类单元的丰度进行半定量评估。因此,基于实时PCR的技术能够靶向特定的分类标记(从高阶分类域到菌株水平),从而有效验证新一代测序所获得的结果。此外,qPCR已成功与其他微生物生态学分子方法结合使用,例如稳定同位素探针技术(SIP)或系统发育/功能基因芯片。与前者结合时,qPCR可用于量化代谢活跃的微生物群落29,30。当与基因芯片分析结合时,qPCR可为基于系统发育标记和功能基因的环境调查提供关键的定量解读31,32。
因此,无论单独使用,还是与其他(通常是基于过程的)生态系统功能评估方法结合使用,定量PCR都是微生物生态学家探索微生物群落与生态系统功能之间难以捉摸的联系所必需的工具。
作者无任何利益冲突需要披露。
本出版物所涉及的研究工作获得了自然环境研究理事会(NERC)资助号 NERC NE/JO11959/1 以及爱尔兰科学基金会的财政支持 & 玛丽·居里行动COFUND项目(资助号:11/SIRG/B2159)资助予CJS,以及东部学术研究联盟(Eastern ARC)
| 姓名 | 公司 | 目录编号 | 评论 |
|---|---|---|---|
| 焦碳酸二乙酯 | Sigma | 40718 | 有毒,需在通风橱中开启 |
| 十六烷基三甲基溴化铵(CTAB) | Sigma | 52365 | 刺激性,需在通风橱中开启 |
| 磷酸氢二钾 | Sigma | RES20765 | |
| 磷酸二氢钾 | VWR | P9791 | |
| 苯酚:氯仿:异戊醇 pH 8 | Sigma | P2069 | 使用前调节至 pH 8 |
| 氯仿:异戊醇 | Sigma | 25666 | |
| 氯化钠 | VWR | 1.06404.0500 | |
| 聚乙二醇 6000 | VWR | 528877-100 | |
| 分子级乙醇 | Sigma | E7148 | |
| Lysing Matrix E 管 | MP Biomedical | 116914050 | |
| Turbo DNA酶 | Ambion | AM1907 | |
| Taq 聚合酶 | Sigma | D1806 | |
| dNTPs | Bioline | BIO-39028 | |
| SuperScript III | Life Technologies | 18080044 | |
| RNase 抑制剂 | Life Technologies | 10777019 | |
| 无RNA酶/DNA酶 0.2 ml PCR管 | Sarstedt | 72.985.002 | |
| SureCleanPlus | Bioline | BIO-37047 | |
| pGEM Easy T 载体 | Promega | A1360 | |
| E. coli JM109 感受态细胞 | Promega | L2005 | |
| 胰蛋白胨 | Sigma | T7293 | |
| 酵母提取物 | Sigma | 70161 | |
| 氨苄青霉素 | Sigma | 59349 | |
| X-Gal | Bioline | BIO-37035 | |
| IPTG | Bioline | BIO-37036 | |
| 琼脂 | VWR | 20768.235 | |
| 质粒中量提取试剂盒 | Qiagen | 12143 | |
| Qubit | Life Technologies | Q33216 | |
| Quant-IT DNA HS 检测试剂 | Life Technologies | Q-33120 | |
| Quant-IT RNA HS 检测试剂 | Life Technologies | Q32855 | |
| MEGAshortscript 试剂盒 | Ambion | AM1354 | |
| TaqMan SensiFast Probe Lo-ROX 试剂盒 | Bioline | BIO-84002 | |
| qPCR 仪 |