简介

[摘要]生产转基因植物既可以研究基因功能，又可以赋予作物所需的性状。众所周知，转基因拷贝数会影响表达水平，从而影响表型。同样，转基因接合性的知识对于对表型进行定量评估和跟踪后代中转基因的遗传是可取的。自从第一批转基因植物问世以来，已经应用了多种确定拷贝数的方法，包括Southern印迹、定量实时PCR ，以及最近的测序方法；ħ H但是，每个方法具有特定的缺点，影响吞吐量，准确性，或费用。数字 PCR ( dPCR ) 将反应分成多个部分，将 PCR 的指数、模拟性质转换为线性数字信号，从而可以准确估计特定序列的出现频率。置信度随着分区数量的增加而增加；因此，能够将反应分成数万个纳米液滴的乳液技术的可用性允许在众所周知的数字液滴 PCR ( ddPCR ) 中准确确定拷贝数。ddPCR具有与其他 PCR 技术类似的低成本和可扩展性优势，但具有卓越的准确性和可靠性。


图文摘要：

数字 PCR ( dPCR ) 将反应分成多个部分，将 PCR 的指数、模拟性质转换为线性数字信号，从而可以准确评估转基因拷贝数的频率。

[背景]植物转化被用作研究基因功能以及生产基因工程作物的研究工具。有几种方法用于递送DNA导入植物细胞，其中包括直接递送到分离的原生质体，生物射弹递送采用高密度微粒，农杆菌tumefac iens介导的递送，以及最近的纳米颗粒（Altpeter等人，2016;王等人。, 2019) 。DNA递送技术可以优化解Ž编，得到的转基因品系具有单个或低更高的频率-拷贝插入; 然而，转移到植物基因组的转基因的确切数量无法控制或限制，必须单独评估转基因品系以确定转基因拷贝数（Sivamani等，2015）。因为表达水平，并因此表型，长期被公认由拷贝数的影响，准确的测定是必需的（霍布斯等人，1993）。拷贝数测定还可用于跟踪遗传和接合性。大多数植物转化方法将转基因传递到体细胞，从中再生转基因植物（称为 T 0 ），每一个都具有一个多半合的转基因基因座。在从自fertili得到的子代Ž通货膨胀，转基因座位通常遵循经典孟德尔遗传; 人虽然，非孟德尔分离也有报道（的Ahuja和Fladung ，2014）。拷贝数的准确测定可用于通过后代跟踪转基因，特别适用于具有多个转基因拷贝的品系和专性异交物种。

自从植物转化成为惯例以来的几十年里，许多不同的技术已被用于确定拷贝数。最初，拷贝数由Southern杂交分析确定(Southern，1975)。在这项技术中，基因组 DNA 用限制性内切核酸酶消化，通过电泳分离，然后转移到膜上，使用标记的 DNA 探针检测含有转基因的片段。除了费力之外，检测可能位于同一片段上的靠近插入（串联重复）的转基因也具有挑战性。此外，除非有足够的杂合度使得同源染色体在限制性消化后具有不同的片段化模式，否则合子性的检测依赖于条带强度的差异，这被认为是不可靠的（Cantsilieris等，2013）。

在 21 世纪的最初几十年，定量实时 PCR (qPCR) 成为估计植物转基因拷贝数的首选方法（Bubner和 Baldwin，2004；Li等，2004；Weng等，2004） ；Yang等人，2005a 和 2005b ；Yi 和 Hong，2019）。将每个样品中转基因的浓度与标准曲线或经过验证的参考基因（已知拷贝数的内源基因）进行比较（Pfaffl ，2001）。然而，因为每个PCR循环是加倍反应，所需的两个倍的增加来区分一个从两个拷贝，或半合子复制从纯合体，是在检测极限。即使优化解Ž通货膨胀，它可以是具有挑战性的，以实现高精度的水平所需区分单拷贝和两拷贝线（Cantsilieris等人，2013; Mieog 。等人，2013年）。

下一步-克eneration小号equencing（NGS）技术已经降低了测序成本; 然而，即使是低通浅层测序，有时也称为“脱脂测序” ，对于大型植物基因组来说也是昂贵的 –对于常规拷贝数检测而言，这是令人望而却步的（Golicz等人，2015 年；Kim等人，2016 年）。目标-捕获测序，其中将含有感兴趣的序列的基因组区域是分离的，已允许的产生与在降低成本的插入位点的序列同一性信息准确拷贝数数据一起（Guttikonda等人，2016） 。最初，短读被应用新一代测序技术（Polkoa等，2012;勒帕等人; 2013年郭等，2016） ，但最近，长期阅读的纳米孔-基于测序平台也已utili ž ED （李等人，2019;布蒂尼。等人，2020） 。然而，与基于PCR 的方法相比，即使是捕获测序技术也具有更高的成本，并且需要更多的时间进行样品制备和数据分析。

数字 PCR ( dPCR ) 将PCR的指数、模拟性质转换为线性数字信号，允许准确估计特定序列的出现频率（Vogelstein和Kinzler ，1999）。为此，每个样品被分成大量分区，并在每个分区上单独进行PCR 反应。这是最初困难和费力的根据需要在多个区室（管）进行划分反应，但乳液技术的发展使反应是compartmentali ž在单管内纳米液滴ED（Diehl的等人，2006）。与 qPCR 协​​议一样，荧光被纳入目标扩增子中；然而，不是在反应的指数阶段监测特定循环的荧光强度，而是确定每个分区在反应终点处是否存在荧光（以及目标）。目标基因的拷贝数是通过确定相对于已知拷贝数的参考基因扩增目标基因的部分的分数来计算的。由于每个分区被打进荧光的存在，DPCR实现精确孔定量吨的转基因拷贝数的通货膨胀。通过增加种群大小（计算更多分区）提高了信心，导致通过dPCR估计拷贝数的准确性被报告为优于 qPCR 并与测序技术相比较（Abyzov等，2012）。近年来，一些商业平台，分反应到使用乳化技术的隔间数万（数滴PCR; ddPCR ）或微小号芯片已成功地应用于准确测定不同种类植物的转基因拷贝数（Glowacka等，2016；Collier等，2017；Giraldo等，2019；Cai等，2020）。下面，我们描述一个ddPCR用于确定转基因拷贝数协议在拟南芥植物（蔡等人，2020） 。

关键字：数字PCR, 微滴式数字PCR, 转基因拷贝数, 植物, 转移DNA
材料和试剂

ddPCR设置和液滴生成
PCR板热封，箔，可刺穿（Bio- R ad，目录号：1814040）
用于 QX200 TM /QX100 TM液滴发生器的DG8 TM墨盒（Bio- R ad，目录号：1864008）
液滴发生器 DG8 TM垫片（Bio- R ad，目录号：1863009）
ddPCR板，96 孔，半裙（Bio- R ad，目录号：12001925）
QubitTM的dsDNA HS测定试剂盒（赛默飞小号系统求解，目录号：Q32851）
限制性核酸内切酶在转基因内具有识别位点，但不在引物位点之间（详见程序）和兼容的消化缓冲液
QX200 TM ddPCR TM EvaGreen Supermix （Bio- R ad，目录号：1864033）
用于EvaGreen ® 的QX200 液滴生成油（Bio- R ad，目录号：1864005）

设备

ddPCR设置和液滴生成
Qubit TM 4荧光计（ThermoFisher Scientific，目录号：Q33238）
QX200 TM液滴发生器（Bio- R ad，目录号：1864002）
热循环仪（例如，Bio- R ad C1000 Touch 热循环仪，目录号：1851148）
多通道移液管小号（10-100微升）
热封机（例如，Eppendorf S100，目录号：5391000036）
DG8 筒架（Bio- R ad，目录号：1863051）
QX200 液滴阅读器（Bio- R ad，目录号：1864003）
液滴读取器板架（Bio-rad ，目录号：12006834）

软件

引物设计软件，例如引物 3 ( https://primer3.ut.ee/ )
序列分析软件例如，Benchling （https://benchling.com）
QuantaSoft软件（Bio- R ad，目录号：1863007）
能够读取逗号分隔值 (.csv) 文件的软件（例如Microsoft Excel）
Droplet Digital TM PCR 应用指南
http://www.bio-rad.com/webroot/web/pdf/lsr/literature/Bulletin_6407.pdf

程序

在此过程中，通过对每个样品进行两次 PCR 反应来确定转基因拷贝数。第一个扩增转基因区域（靶基因），第二个扩增已知拷贝数的内源基因（参考基因）区域。另一种方法是使用水解探针检测目标基因和参考基因的多重扩增；之前已经报道了几种植物参考基因的引物和探针（Collier等，2017）。

限制性内切酶的引物设计与选择
使用引物设计软件（例如，Primer 3）选择可扩增转基因目标基因和参考基因的 100-120 碱基对区域的引物。引物应最佳地已退火〜60的温度下℃下，类似的估计结合动力学，以及最小dimeri ž通货膨胀。
检查每个引物对是否在测试 PCR 反应中有效地扩增了所需的目标。包括与 ( i ) 来自野生型非转基因植物的基因组 DNA 和 (ii) 无 DNA 的对照反应。组装反应如表 1 所示。

表 1. 用于测试ddPCR引物对的PCR 反应

使用表 2 中提供的试剂制造商推荐的循环条件孵育 PCR 反应。

表 2. 用于测试ddPCR引物对的PCR 循环条件。对所有步骤使用2°C/s 的升温速率。

Visuali ž ËPCR产物通过琼脂糖凝胶电泳，从而确保这两组引物获得了预期大小的单条带清晰。
在扩增前消化基因组 DNA 将确保，即使有两个转基因副本串联插入，每个副本也将位于单独的 DNA 分子上，因此位于单独的液滴中。ANALY ž e是转基因序列为限制性内切酶（RE）使用识别序列的存在的序列分析软件例如，Benchling 。选择在转基因中至少有一个识别位点的n RE，但在目标或参考基因的引物之间没有位点。避免使用对甲基化敏感的限制性内切核酸酶。

ddPCR设置
ddPCR需要高质量的基因组 DNA 。甲十六烷基三甲基溴化铵（CTAB）为基础的方法或市售的植物DNA提取试剂盒将提供足够纯的基因组DNA。每个ddPCR反应通常需要 5-50 ng 基因组 DNA。
注意小号：
为了准确孔定量泰特DNA浓度，我们建议使用一个量子位荧光计。
在每个反应中所用的基因组DNA的确切数量必须是优化解ž的生物体的基因组大小编辑。
使用选定的 RE 消化每个基因组 DNA 样本。一个例子限制性消化反应示于表3.不同量的酶和缓冲液将根据制造商和浓度需要的酶。目标是每 g 基因组 DNA包含 ~10-30 U 酶。
注意：EvaGree n ®产品说明表明，也可以通过向PCR 反应中加入 1 μl (2-5 U) 所需的限制性内切酶来消化基因组 DNA 。

表 3. 基因组 DNA 限制性内切酶消化的示例反应

通过加热灭活 RE。温育的精确温度和长度将取决于在特定的酶和制造商。灭活后，消化的基因组 DNA 可直接用于下一步，无需进一步清理。
按照表 4 中所示的顺序组装 PCR 反应。除了样品反应之外，还设置了与 ( i ) 来自野生型非转基因植物的基因组 DNA 和 (ii) 无 DNA 的对照反应。在进行液滴生成之前，组装好的ddPCR反应必须充分混合，因为每个液滴必须包含相同浓度的每种成分。

表 4. ddPCR反应示例

注意：如果多个样品正在ANALY Ž版，作为其他的PCR方法，建议向组装使主混合物和稀浓缩DNA之前ddPCR反应。

液滴生成
按下中间的闩锁，打开 DG8 TM墨盒固定器。
将 DG8 TM Cartridge 插入支架并通过将两半压在一起来关闭支架（图 1）。
图 1.液滴生成的工作流程。PCR 反应和油被添加到墨盒中，墨盒插入支架并由垫圈固定到位。液滴生成后，乳化反应被转移到微量滴定板进行热循环。

加载每个20 -在先前步骤中成DG8的单个样品池组装微升反应TM使用多道移液器墨盒。
将 70 μl QX200 用于EvaGreen ® 的液滴生成油加入油井中。
用 DG8 TM垫圈盖住滤芯。将垫圈的孔绕在墨盒支架上的四个挂钩上，以将垫圈固定到位。
打开 DG8 液滴发生器并插入墨盒支架。按下按钮关闭盖子。发电机将自动启动。
将整个体积的液滴乳液（通常为 40 μl ）转移到半裙边 96 孔ddPCR板的所需孔中（图 1）。盖上刺穿箔片热封的板和附连使用根据一个板密封器的制造商的说明。旋转板并重复此步骤以确保箔完全贴合。
注意：为避免破坏乳化液滴，请注意移液速度非常缓慢。可以使用大口径移液器吸头（或使用刀片去除末端的移液器吸头）。
密封板放置到热循环仪和我们Ë完全相同的循环条件为那些在测试PCR中使用（见表2）。在最初的热循环过程中，每个液滴周围的油都会热固。

液滴计数
PCR 循环完成后，打开 QX200 TM液滴阅读器并打开QuantaSoft TM软件。
在QuantaSoft TM软件中，开始一个新实验（模板 > 新建）。
选择要检测的孔并输入样品名称或记录每个样品的目标和参考基因的孔位置。
在井编辑器中选择关键参数如下：
实验：ABS（绝对孔定量牛逼通货膨胀）。
注：QuantaSoft TM软件提供了三种默认实验：绝对孔定量牛逼通货膨胀（ABS），稀有事件检测（RED）和拷贝数变异（CNV）。ABS 和 CNV 均可用于确定拷贝数。对于 CNV 模式，需要提供已知的参考基因拷贝数。ABS模式更加灵活。例如，拷贝数的不需要参考基因，并且可以计算目标和参照的相对比率。
Supermix ：QX200 TM ddPCR EvaGreen Supermix 。
目标 1 类型： Ch1 Unknown 。
目标 2 类型： Ch2 Unknown 。
名称：é NTER样品名称例如，转基因行x，参照基因。这有助于区分哪些孔用于转基因，哪些孔用于参考基因。
单击“确定”将所选参数应用于所选孔。
从热循环仪中取出 PCR 板并插入读板器支架的底座。盖上盖子上的板，然后按所述闩锁保持器。
打开 QX200 TM液滴读取器，插入读板器支架，然后合上盖子。
检查“电源”、“瓶位”和“盘子就位”的指示灯是否为绿色。
单击QuantaSoft TM软件中的“运行” 。在“运行选项”对话框中选择EvaGreen 。单击“确定”开始阅读。

数据分析

阅读完成后，单击“分析”。
单击“一维振幅”按钮以查看正负液滴的分布。示例分布如图 2 所示。

图 2. 示例图显示了八个样品 (A01-H01) 的阳性（蓝色）和阴性（灰色）液滴。X 轴（事件编号）表示整个实验中测量的液滴数量，Y 轴表示荧光幅度。

注意：单个反应应产生 12,000-20,000 个液滴。的拷贝数的计算依赖于正的比率到负液滴。如果阳性液滴饱和（大多数液滴是阳性的），这意味着每个液滴都有目标分子，很可能引物是非特异性的，或者反应中使用了过多的 DNA。为了获得准确的结果，应检查引物特异性和/或用较少的 DNA 重复反应。
的阈值，这区别于负液滴正液滴，以及浓度值，则自动地确定; ħ H但是，如果需要的话，该阈值可以手动调节。
单击“浓度”按钮。将浓度数据导出为逗号分隔值(.csv) 文件。
在打开.csv文件的适当的软件（例如，Microsoft Excel）中。QuantaSoft TM软件中提供的任何样品名称都将出现在“样品”列中。在“ ç oncentration”列包含每PCR扩增子的数量微升需要计算拷贝数。
靶基因的拷贝数是通过计算（ç目标/ Ç REF ）* Ñ REF ，其中Ç目标是转基因的浓度，Ç REF是内源基因的浓度，和Ñ REF是的拷贝数参考基因。在表5提供的例子中，参考基因是拟南芥基因组中的单拷贝纯合基因；因此，在每个提取 DNA 的二倍体体细胞中都有两个拷贝。在第 1 行中，目标 ( C target )的浓度为 115 拷贝/ μl ，参考浓度 ( C ref ) 为 227 拷贝/ μl 。因此，转基因的拷贝数是（227分之115）* 2 = 1.0，这意味着该线1具有转基因（半合子）的单个副本。在第 2 行中，C target = 784 拷贝/ μl和C ref = 173 拷贝/ μl 。因此，转基因的拷贝数是（173分之784）* 2 = 9.1 ，这意味着该线2具有转基因的9份。在第 3 行中，C target = 11.7 拷贝/ μl和C ref = 333 拷贝/ μl 。因此，对拷贝数的目标转基因是（11.7 / 333）* 2 = 0.07份，这意味着线3是不太可能包含转基因。
注意：在确定多倍体植物中的转基因拷贝数时，许多通常表现为单拷贝基因的内源基因可能存在于多个基因组中。例如，在四倍体植物的二倍体细胞中，典型的单拷贝参考基因可能有四个拷贝；因此，单拷贝半合子转基因的浓度值将是参考基因的四分之一。上述方法也可用于确定后代中的转基因接合度；ħ H但是，因为该ddPCR协议将只提供原始值，父行（多个）的拷贝数的知识是必需的。例如，如果估计后代植物的二倍体细胞含有六个拷贝，则需要了解任何亲本系的拷贝数以确定植物是否具有六个独立的半合子位点或三个纯合子位点。

表 5.来自三个转基因植物（第 1-3 行）的示例ddPCR数据，显示了扩增转基因（目标）和已知拷贝的内源基因（参考）时获得的阳性和阴性液滴的数量。显示浓度（每微升PCR 扩增子数）的列用于计算拷贝数。

致谢

我们感谢 Oleg Raitskin的技术帮助和 Simon Foster 的设备维护。我们感谢英国生物技术和生物科学以及工程和物理研究委员会（BBSRC 和EPSRC ）合成生物学增长计划（OpenPlant合成生物学研究中心 BB / L 014130/1）的财政支持，以及 BBS/E/ T/000PR9816 和 BB/R021554/1。该协议旨在生成 Cai (2020) 中提供的转基因拷贝数数据。

利益争夺

没有要申报的竞争利益。