采用优化的数字PCR方法分析转基因小麦外源基因拷贝数

奥鹏网院作业

采用优化的数字PCR方法分析转基因小麦外源基因拷贝数
琚鹏举1，宁蕾1，葛林豪2，许成杰1，史华伟1，梁凯歌3，马亮4，刘陶然2，陈明2，孙黛珍1

（1山西农业大学农学院，山西太谷 030800；2中国农业科学院作物科学研究所，北京 100081；3中央民族大学，北京 100081；4石家庄市农林科学研究院，石家庄 050047）

摘要：【目的】探索基于数字PCR的小麦基因拷贝数分析技术，提高目标基因拷贝数的分析通量，促进小麦基因组研究及基因工程研究。【方法】采用中国农业科学院作物科学研究所小麦抗逆分子育种课题组创制的高抗小麦黄花叶病转nib8小麦为试验材料，使用小麦D基因组中的单拷贝纯合内源基因PINb-D1b（籽粒硬度基因）为内参基因，根据转化的抗病基因nib8序列设计4对特异性引物及相应探针，通过试验确定探针和引物的最优浓度，优化体系的退火温度，寻找最合适的模板浓度。然后用数字PCR检测转nib8小麦中外源基因的拷贝数；同时用Real-time PCR和Southern blot 2种不同方法分析的拷贝数结果，验证上述采用数字PCR方法分析拷贝数的准确性；以PINb-D1b内参基因检测Wx012（蜡质基因）和SSII（淀粉合成基因）2个内参基因的拷贝数，以Wx012和SSII为内参基因检测转nib8小麦中外源基因的拷贝数，验证以PINb-D1b内参基因的检测结果准确性；在nib8的不同区段设计特异性引物来检测转nib8株系拷贝数分析结果是否有差异。最终，建立基于数字PCR技术的高通量检测小麦目标基因拷贝数的分析方法。【结果】通过试验最终确定了基于数字PCR方法的小麦目标基因拷贝数的检测体系。试验确定引物和探针的最佳终浓度分别为500和250 nmol·L-1，确定体系反应最佳退火温度为59℃，最佳DNA模板量为40 ng。以PINb-D1b作为内参基因检测转nib8小麦中nib8的拷贝数，拷贝数检测结果显示第12、16、17、23、29、30株系中nib8拷贝数分别为7个、1个、1个、1个、1个、7个；同时通过比较发现此结果与Real-time PCR以及Southern blot结果一致；在转基因株系中，以PINb-D1b为内参基因，对Wx012和SSII 2个基因进行拷贝数进行分析，同时发现采用这3种内参基因分析转nib8小麦中nib8基因拷贝数的结果一致，说明这3个内参基因都适合用于数字PCR方法；在nib8上游、中游和下游不同区段的设计引物分析拷贝数检测结果一致。【结论】优化了基于数字PCR方法的小麦目标基因拷贝数分析方法和反应体系，确立了基于数字PCR方法的小麦目标基因拷贝数的检测体系，数字PCR分析结果稳定、可靠，检测通量明显提高，具有一定的应用前景。

关键词：小麦；基因拷贝数；数字PCR；内参基因

0 引言
【研究意义】目的基因的拷贝数一直是植物基因组研究的重要内容。外源基因整合进入受体基因组的位置和拷贝数会影响目的基因的表达以及遗传稳定性，而插入基因的拷贝数相对其插入位点而言影响更为重要[1]。当外源基因以低拷贝数（1—2个）整合到受体基因组时，通常能够稳定高效转录表达，多拷贝数的整合则会造成基因不稳定表达，甚至沉默[2]。因此，鉴定转基因产品的外源基因插入拷贝数非常重要[3]。由于小麦基因组庞大复杂，优化和建立高通量拷贝数分析方法对于对小麦转基因育种研究具有重要意义。【前人研究进展】近年来，微滴式数字 PCR（droplet-based digital PCR，ddPCR）是最新兴起的一项核酸检测技术[4-6]。该技术不依赖任何校准物，采用直接计数单个分子的原理，实现了对核酸检测的绝对定量[7-8]。ddPCR已被广泛应用于生命科学研究的多个领域，如核酸绝对定量、稀有突变检测和拷贝数检测等[9]。ddPCR具有高精确度的特点，在检测目的基因拷贝数时，利用目标基因与参考基因的双重反应，通过计算它们的比值，即可快速获得目标基因拷贝数。这种方法由于试验结果重复性好，已经成为了检测基因拷贝数的首选方法[10]。与传统方法相比较，ddPCR具有明显优势，主要体现在操作流程更加简单，无需构建标准曲线，且是对待检测的核酸分子直接进行绝对定量，试验结果更加准确可靠[11]。此外，ddPCR在线性范围、检测极限和定量极限等方面均优于qRT-PCR方法[12]。诸多研究表明，ddPCR已用于玉米、水稻等作物转基因插入拷贝数的检测[13-14]。SUN等[15]的研究表明运用 ddPCR技术可以在甘蔗复杂基因组中明确确定内源参考基因拷贝数，COLLIER等[16]研究中使用双重ddPCR获得了水稻、小麦、玉米、番茄、马铃薯和柑桔类物种的高质量转基因拷贝数测量值。其中提到 ddPCR分析产生的拷贝数测量值非常接近整数值，通常与使用Southern blot杂交得出的结论相近。并且，基于ddPCR的方法在基因组非常大的物种（如小麦和玉米）中也能很好地发挥作用，但传统的DNA印迹杂交方法在这一方面却极富挑战性。因此快速准确地确定转基因拷贝数并鉴定大量样品的半合子和纯合子的能力使 ddPCR成为植物生物学研究人员的强大工具。【本研究切入点】六倍体小麦因其基因组庞大复杂（约17 Gb），利用传统Southern blot等方法检测拷贝数费时费力，无法达到高通量检测的要求，当需要对大量的转化试验进行拷贝数分析时，传统方法无法满足要求。【拟解决的关键问题】本研究以中国农业科学院作物科学研究所小麦抗逆分子育种课题组利用农杆菌介导法获得的转nib8小麦为试验材料，以小麦PINb-D1b（籽粒硬度）为参考基因[15]，期望通过数字 PCR技术建立一种快速稳定的小麦基因拷贝数检测方法，为促进小麦基因组研究以及基因工程研究提供有力工具。

1 材料与方法
1.1 试验材料
2018年10月，中国农业科学院顺义基地种植转nib8小麦株系（由中国农科院作物科学研究所小麦抗逆分子育种课题组提供）和扬麦16（受体）。

1.2 内参基因和目标基因引物设计
选取位于小麦 D基因组中的单拷贝纯合基因PINb-D1b（籽粒硬度）作为试验内参基因[17]，通过引物的筛选和比较，确定专一性和特异性强的引物（表1）。设计Wx012（蜡质基因）和SSII（淀粉合成酶）2个内参基因的引物和探针，确保探针的特异性（表1）。

表1 数字PCR 探针引物序列
Table 1 Sequences of probes and primers used in digital PCR

pagenumber_ebook=9,pagenumber_book=1933
1.3 样本DNA的提取
结合CTAB法和全自动核酸提取仪提取法提取样本DNA[18]。

1.4 片段化基因组
用EcoRⅠ（NEB公司）对叶片基因组DNA进行核酸内切酶预处理，每微克DNA加入1个单位EcoRⅠ。37℃处理过夜，65℃ 30 min使酶失活，并将模板DNA稀释为 20 μg·μL-1。

1.5 制备反应体系及PCR反应
采用以DNA为模板的探针法。反应体系见表2，反应条件为95℃ 3 min；95℃ 30 s，57℃ 1 min，40个循环；98℃ 10 min。为保证高效的模板扩增，PCR仪的升降温速度需低于2℃·s-1。

表2 PCR扩增反应体系
Table 2 PCR amplification reaction system

pagenumber_ebook=9,pagenumber_book=1933
1.6 微滴分析
将 96孔板放入微滴分析仪中，读取微滴荧光信号，从而判断微滴的阳性或阴性。使用 QuantaSoft V1.3.2软件计算每个样品目的序列的微滴数（拷贝/μL）[19]，泊松分布的置信区间为 95%；然后根据内参基因的拷贝数计算待测目的基因的拷贝数。

2 结果
2.1 引物的特异性筛选
选用转nib8小麦株系为模板，使用表1的4个基因引物进行PCR扩增（图1）。结果发现扩增出的条带均为单一目的条带，表明引物具有特异性。

pagenumber_ebook=10,pagenumber_book=1934
图1 4对引物的特异性扩增
Fig. 1 Specific amplification results of four pairs of primers

M：DL2000Maker；1：nib8；2：PINb-D1b；3：Wx012；4：SSII；5：H2O

2.2 退火温度的筛选
通过对退火温度的梯度分析，在6个温度条件下，nib8均表现出良好的扩增图形，且随着退火温度的升高，nib8的阴性微滴和阳性微滴区分的愈发明显，扩增越来越特异（图2-A）。图2-B中内参基因随着退火温度的升高，扩增效率变化不明显。因此，为使nib8和内参基因均能表现出良好的扩增效果，最终以59℃作为最优退火温度。

2.3 灵敏度（检测低限）和准确性检测
2.3.1 DNA模板浓度的筛选 由图3可知，nib8在4个不同模板浓度下均扩增出良好的微滴分离现象，随着浓度的增加，微滴数增加。然而，nib8在160 ng模板浓度时，阳性微滴和阴性微滴已无法明显分开，在320 ng模板浓度时，该现象愈发明显。模板浓度越高，体系扩增效果越差，微滴不能明显分开，而模板为20 ng时的微滴数较少，无法明显看出扩增效果，因此，最终选定40 ng为最佳模板浓度。

2.3.2 特异性检测 使用QuantaSoft V1.3.2软件对6个转nib8小麦株系的拷贝数进行分析。由图4可以看出阳性微滴数和阴性微滴数区分明显，说明体系稳定。由表3可以看出数字PCR检测结果中有3个株系为1个拷贝，可以用于进一步基因功能分析。

2.4 体系验证及应用
2.4.1 不同检测方法的比较 由表 4中可以看出数字 PCR的检测结果与通过染料法和探针法进行的real-time PCR分析结果一致。由图5可以看出转nib8小麦第16号株系的Southern blot结果也与数字PCR结果一致，均是1个拷贝。

pagenumber_ebook=11,pagenumber_book=1935
图2 不同退火温度条件下nib8（A）和内参基因（B）的数字PCR检测
Fig. 2 Digital PCR detection of nib8 (A) and internal reference genes (B) at different annealing temperatures

蓝点为nib8阳性，绿点为内参基因阳性，灰点为阴性，红线为荧光阈值，从左往右依次为56℃、56.9℃、58.1℃、59℃、59.6℃、60℃、对照
The blue point was positive for nib8, the green point was positive for internal reference gene, and the gray point was negative. The red line was the fluorescence threshold, which was 56℃, 56.9℃, 58.1℃, 59℃, 59.6℃, 60℃ and the control from left to right

pagenumber_ebook=11,pagenumber_book=1935
图3 不同模板浓度条件下nib8（A）和内参基因（B）的数字PCR检测
Fig. 3 Digital PCR results of nib8 gene(A) and control gene (B) at different template concentrations

蓝点为内参基因阳性，绿点为nib8阳性，灰点为阴性，红线为荧光阈值限，从左往右依次为20（A01）、40（B01）、160（C01）和320 ng（D01）模板及水
The blue dots were positive for the internal reference gene, the green dots were positive for nib8, the gray dots were negative, and the red line was the fluorescence threshold. From left to right, the template and water were 20 (A01), 40 (B01), 160 (C01) and 320 ng (D01)

pagenumber_ebook=11,pagenumber_book=1935
图4 转nib8株系特异性检测
Fig. 4 Specific detection of nib8 strains was performed

从左到右为不同的转基因株系、阴性对照和空白对照
Different transgenic lines, negative control and blank control were shown from left to right

表3 转nib8株系的数字PCR检测低限测定结果
Table 3 The results of low limit determination by digital PCR were obtained by transplanting nib8 strain

pagenumber_ebook=11,pagenumber_book=1935
2.4.2 内参基因的选择 在内参基因的选择上，除小麦D基因组中的单拷贝纯合基因（PINb-D1b）外，同时还开发了其他小麦内源基因作为拷贝数检测的内参基因。

以Wx012和SSII为内参基因检测部分转nib8小麦株系目标基因的拷贝数（图6和图7），其结果与用PINb-D1b为内参基因的分析结果一致，都是1个拷贝，证明这3个基因都可作为内参基因用于小麦目标基因拷贝数数字PCR分析。

pagenumber_ebook=12,pagenumber_book=1936
图5 转nib8小麦第16号株系的Southern blot检测
Fig. 5 Southern blot results of nib8 transgenic wheat line 16

表4 探针法、染料法和数字PCR的结果对比
Table 4 Comparison between fluorescence quantitative method such as probe method and dye method, and digital PCR method

pagenumber_ebook=12,pagenumber_book=1936
pagenumber_ebook=12,pagenumber_book=1936
图6 内参基因的拷贝数分析
Fig. 6 Copy number analysis of internal reference genes

A1和B1是以PINb-D1b为内参基因时Wx012的拷贝数；C1和D1是以Wx012为内参基因时SSII的拷贝数；E1和F1是阴性对照和空白对照
A1 and B1 are the copy numbers of Wx012 when pinb-d1b is the internal reference gene.C1 and D1 are the copy numbers of SSII when Wx012 is the internal reference gene; E1 and F1 are negative controls and blank controls

pagenumber_ebook=12,pagenumber_book=1936
图7 转nib8株系的拷贝数分析
Fig.7 Copy number analysis of transferred nib8 lines

A1和B1分别以Wx012和SSII为内参的转nib8第12号株系的拷贝数分析；C1和D1分别以Wx012和SSII为内参的转nib8第17号株系的拷贝数分析；E1和F1是阴性对照和空白对照
Copy number analysis of transferred nib8 line no. 1pagenumber_ebook=12,pagenumber_book=19362 with Wx012 and SSII as internal parameters for A1 and B1pagenumber_ebook=12,pagenumber_book=1936The copy number analysis of C1 and D1 withWx012 and SSII as internal parameters, respectivelyE1 and F1 are negative controls and blank controls

2.4.3 靶基因不同区段的检测效果验证 在nib8的不同区段进行引物探针设计，分别以PINb-D1b、Wx012和SSII为内参基因检测nib8拷贝数，结果均一致（图8-A，图8-B），第12号转nib8小麦株系的拷贝数都是7（图8-C）。

pagenumber_ebook=13,pagenumber_book=1937
图8 在nib8不同位置设计引物的结果图（第12号转nib8株系）
Fig.8 Results of designing primers at different positions of nib8 gene (take 12 strains as an example)

A和B分别是用nib8不同区间引物对第12号转nib8株系内参基因和目标基因的数字PCR检测；C：第12号转nib8株系的拷贝数检测
A and B were used to detect the internal reference gene and target gene of the no.12 transferred nib8 line by digital PCR with primers of different regions of nib8; C: Copy number detection of No. 12 to nib8 line

3 讨论
3.1 不同拷贝数检测方法的比较
目前，小麦拷贝数检测方法很多，但用数字PCR检测转基因小麦外源基因拷贝数报道较少。经典的外源基因拷贝数检测方法主要有 qRT-PCR和 Southern blot。其中qRT-PCR分析方法的步骤较为繁琐，首先需要准备标准 DNA样品，然后利用此样品构建标准曲线，而标准曲线的建立又涉及体系的摸索与优化，需要耗费大量的时间和精力，所需试验周期较长[20]；另外，借助标准曲线进行定量本身就是一种相对定量方法，所以试验结果并不准确[21]。而传统的Southern blot方法不但对DNA需求量大，试验操作步骤复杂，而且对人员操作的技术要求比较高[22]。与定量PCR不同，数字PCR无需标准曲线即可实现对DNA样品绝对定量的灵敏、准确。一个PCR反应被分成许多独立的反应，每个反应都有一个阳性或阴性的信号。应用泊松统计量，直接从阳性反应和阴性反应的数量计算出原始样品中DNA分子的数量。本研究中数字PCR模板要求低，而且相比Real-time PCR的灵敏性（0.1 pg·μL-1）高，数字 PCR 灵敏性能达到 0.001 pg·μL-1，是Real-time PCR的100多倍。相比Southern blot，数字PCR试验周期短，操作简单，检测通量高。

3.2 体系准确性的比较
一步多重PCR（multiplex PCR）又称多重引物PCR或复合PCR，是指在同一PCR反应体系中添加2对以上引物，同时扩增出多个核酸片段的PCR反应，且每对引物只对应1个[23]靶标片段，与常用的RAPD技术（1对引物对应多个扩增产物）[24]全然不同。考虑到二重微滴式数字PCR不仅能避免单重ddPCR定量因二次取样所产生的样品内的误差，而且可以消除不同样品 DNA因取样不一致而造成的样品间的误差，所以本研究采用二重PCR反应扩增，选用了3个不同的小麦单拷贝基因Wx012、PIN-b1和SSII作内参基因，其中以PIN-b1为内参基因时，Wx012和SSII的拷贝数为1（图6）。又以这3个基因为内参基因检测转基因株系12的拷贝数时，结果也相同。为了检验表达是否一致，在nib8序列的不同区段设计引物，结果显示一致。为了试验结果的准确性，本研究选择小麦基因组里的其他2个单拷贝基因Wx012和SSII作为内参基因检测PIN-b1[25]。以此来检测该体系的准确性。

3.3 数字PCR检测小麦拷贝数方法的应用
本研究建立的检测体系短时间内可以检测几百个样品，大大提高了工作效率。本研究建立的拷贝数方法不仅可以检测外源基因。也可以检测内源基因的拷贝数。因此，本方法具有快速、灵活、高效、所需样本量少等优点，而且已经在基因表达研究、基因组拷贝数鉴定、癌症标志物稀有突变检测、转基因成分鉴定等诸多领域显示出广阔的应用前景[26-28]。BURNS等[29]应用数字 PCR来评估几种转基因分析方法的检测限。

DOBNIK等[30]应用数字 PCR技术，针对欧盟授权的 12种转基因玉米品系建立了转基因成分多重定量分析方法。目前，该技术已经快速应用于各个领域。

4 结论
建立了针对小麦双重微滴式数字 PCR检测基因拷贝数的方法。确定了引物和探针的最佳终浓度分别为500和250 nmol·L-1，退火温度为59℃。该检测方法可以选用3个内参基因（Wx012、PIN-b1和SSII）。

References

[1] WANG J W, WANG J J, ZHAO S J, LI N, LI S J, ZHANG C.Detecting copy number of foreign gene in human alpha-lactalbum in transgenic cloned cattle. China Animal Husbandry & Veterinary Medicine, 2013, 40(9): 183-186.

[2] HEYRIES K A, TROPINI C, VANINSBERGHE M, DOOLIN C,PETRIV O I, SINGHAL A, LEUNG K, HUGHESMAN C B,HANSEN C L. Megapixel digital PCR. Nature Methods, 2011(8):649-651.

[3] MARTíNEZ-GARCíA E, DE LORENZO V. The quest for the minimal bacterial genome. Current Opinion in Biotechnology, 2016,42: 216-224.

[4] CILLONI D, PETITI J, ROSSO V, ANDREANI G, DRAGANI M,FAVA C, SAGLIO G. Digital PCR in myeloid malignancies: Ready to replace quantitative PCR?. International Journal of Molecular Sciences, 2019, 20(9): 2249.

[5] LI H, BAI R, ZHAO Z, TAO L, MA M, JI Z. Application of droplet digital PCR to detect the pathogens of infectious diseases. Bioscience Reports, 2018, 38(6): 237-245.

[6] POMARI E, PIUBELLI C, PERANDIN F, BISOFFI Z. DIGITAL PCR: A new technology for diagnosis of parasitic infections. Clinical Microbiology and Infection, 2019, 25(12): 1510-1516.

[7] 柳方方, 张玲, 王晶. 数字PCR测定DNA含量及测量结果不确定度评定. 化学分析计量, 2013, 22(1): 18-22.LIU F F, ZHANG L, WANG J. Determination of DNA content by digital PCR and evaluation of uncertainty of measurement results.Chemometrics, 2013, 22(1): 18-22. (in Chinese)

[8] MILAVEC M, DOBNIK D, YANG L. GMO quantification: Valuable experience and insights for the future. Analytical and Bioanalytical Chemistry, 2014, 406(26): 6485-6497.

[9] COLLIER R, DASGUPTA K, XING Y P, HERNANDEZ B T, SHAO M, ROHOZINSKI D, KOVAK E, LIN J, LUIZA M, OLIVRIRA D,STOVER E, MCCUE K F, HARMON F G, BLECHL A, THOMSON J G, THOMSON, THILMONY R. High-throughput droplet digital PCR system for absolute quantitation of DNA copy number.Analytical Chemistry, 2011, 83(22): 8604-8610.

[10] SEDLAK R H, JEROME K R. Viral diagnostics in the era of digital PCR. Diagnostic Microbiology and Infectious Disease, 2012, 75(1):1-4.

[11] DHONT A, ISON D, ALIX K. Determination of basic chromosome numbers in the genus Saccharum by physical mapping of ribosomal RNA genes. Genome, 1998, 41(2): 221-225.

[12] 陈嘉茵. 逆转录微滴数字PCR技术检测蔬果中GⅠ型及GⅡ型诺如病毒[D]. 广州: 暨南大学, 2018.CHEN J Y. Reverse transcription droplet digital PCR technique to detect fruits and vegetables in G Ⅰ type and G Ⅱ type, such as virus[D]. Guangzhou: Jinan University, 2018. (in Chinese)

[13] 夏清燕. 柑橘溃疡病微滴数字PCR与定量PCR检测方法的比较研究[D]. 重庆: 重庆大学, 2017.XIA Q Y. Comparative study of quantitative PCR and digital PCR for detection of citrus canker disease [D]. Chongqing: Chongqing University, 2017. (in Chinese)

[14] 周圆. 运用微滴式数字 PCR技术检测转基因玉米品系的研究. 安徽农业科学, 2018, 591(14): 175-178.ZHOU Y. Study on detection of transgenic maize strains by micro-drop digital PCR. Anhui Agricultural Science, 2018, 591(14):175-178. (in Chinese)

[15] SUN Y, JOYCE P A. Application of droplet digital PCR to determine copy number of endogenous genes and transgenes in sugarcane. Plant Cell Reports, 2017, 36: 1775-1783.

[16] COLLIER R, DASGUPTA K, XING Y P, HINDSON B J, NESS K D,MASQUELIER D A, BELGRADER P, HEREDIA N, MAKAREWICZ A J, BRIGHT I J, LUCERO M Y, HIDDESSEN A L, LEGLER T C,KITANO T K, HODEL M R, PETERSEN J F, WYATT P W,STEENBLOCK E R, SHAH P H, BOUSSE L, JTROUP C B,MELLEN J C, WITTMANN D K, ERNDT N G, CAULEY T H,KOEHLER R T, SO A P, DUBE S, ROSE K A, MONTESCLAROS L,WANG S L, STUMBO D P, HODGES S P, ROMINE S,MILANOVICH F P, WHITE H E, REGAN J F, KARLIN G A,HINSON C M, SAXONOV S, COLSTON B W. Accurate measurement of transgene copy number in crop plants using droplet digital PCR. The Plant Journal, 2017, 90(5): 1014-1025.

[17] 朱鹏宇, 张亮亮, 杜智欣. 双重数字 PCR在转基因水稻检测中的应用. 植物检疫, 2018, 32(4): 52-56.ZHU P Y, ZHANG L L, DU Z X. Application of double digital PCR in detection of transgenic rice. Plant Quarantine, 2018, 32(4): 52-56. (in Chinese)

[18] LI Z, HANSEN J L, YING L. Using real-time PCR to determine transgene copy number in wheat. China Biotechnology, 2010, 22(2):179-188.

[19] COLLIE R, DASGUPTA K, XING Y P. Accurate measurement of transgene copy number in crop plants using droplet digital PCR. The Plant Journal, 2017, 90(5): 1014-1025.

[20] WHALE A S, HUGGETT J F, COWEN S. Comparison of microfluidic digital PCR and conventional quantitative PCR for measuring copy number variation. Nucleic Acids Research, 2012,40(11): e82.

[21] 姜羽, 胡佳莹, 杨立桃. 利用微滴数字 PCR分析转基因生物外源基因拷贝数. 农业生物技术学报, 2014, 22(10): 1298-1305.JIANG Y, HU J Y, YANG L T. Analysis of exogenous gene copy number in genetically modified organisms by microdroplet digital PCR. Journal of Agricultural Biotechnology, 2014, 22(10): 1298-1305.(in Chinese)

[22] YANG L, DING J, ZHANG C. Estimating the copy number of transgenes in transformed rice by real-time quantitative PCR. Plant Cell Reports, 2005, 23(11): 759-763.

[23] CHAMBERLAIN J S, GIBBS R A, RAINER J E. Deletion screening of the Duchenne muscular dystrophy locus via multiplex DNA amplification. Nucleic Acids Research, 1988, 16(23): 11141-11156.

[24] EDWARDS M C, GIBBS R A. Multiplex PCR: Advantages,development, and applications. Genome Research, 1994, 3(4): 65-75.

[25] YANG L, QUAN S, ZHANG D. Endogenous reference genes and their quantitative real-time PCR assays for genetically modified bread wheat (Triticum aestivum L.) detection. Methods in Molecular Biology,2017, 1679: 259-268.

[26] OTTESEN E A, HONG J W, QUAKE S R. Microfluidic digital PCR enables multigene analysis of individual environmental bacteria.Science, 2006, 314(5804): 1464-1467.

[27] WHITE T B, MCCOY A M, STREVA V A. A droplet digital PCR detection method for rare L1 insertions in tumors. Mobile DNA, 2014,5(1): 30.

[28] YUNG T K F, CHAN K C A, MOK T S K. Single-molecule detection of epidermal growth factor receptor mutations in plasma by microfluidics digital PCR in non-small cell lung cancer patients.Clinical Cancer Research, 2009, 15(6): 2076-2084.

[29] BURNS M J, BURRELL A M, FOY C A. The applicability of digital PCR for the assessment of detection limits in GMO analysis.European Food Research and Technology, 2010, 231(3): 353-362.

[30] DOBNIK D, BJORN S, ALEXANDRA B. Multiplex quantification of twelve EU authorized GM maize lines with droplet digital PCR.Analytical Chemistry, 2015, 87(16): 8218-8226.

Analysis of Foreign Gene Copy Number in Transgenic Wheat by Optimized Digital PCR

JU PengJu 1, NING Lei 1, GE LinHao 2, XU ChengJie 1, SHI HuaWei1, LIANG KaiGe 3, MA Liang 4,LIU TaoRan 2, CHEN Ming 2, SUN DaiZhen 1
(1 College of Agriculture, Shanxi Agricultural University, Taigu 030800, Shanxi; 2 Institute of Crop Science, Chinese Academy of Agricultural Sciences, Beijing 100081; 3 Minzu University of China, Beijing 100081; 4 Shijiazhuang Academy of Agricultural and Forestry Sciences, Shijiazhuang 050047)

Abstract: 【Objective】In order to explore the analysis technology of wheat gene copy number based on digital PCR, and improve the analysis efficiency of target gene copy number, and promote the research of wheat genome and gene engineering.【Method】Using the transgenic wheat with high resistance to wheat yellow mosaic disease transformed by nib8 gene created by the wheat stress resistance molecular breeding group in the Crop Science Institute of Chinese Academy of Agricultural Sciences as the experimental material, and used the single copy homozygous endogenous gene PIN-D1b (grain hardness gene) in wheat D genome as the internal reference gene, and four pairs of specific primers and corresponding probes were designed according to the sequence of the transformed disease resistance gene nib8, and the optimal concentration of probe and primer, the optimal annealing temperature of the system, to find the most appropriate template concentration were determined through experiments. Then, the copy number of nib8 genes in the nib8 transgenic wheat was detected by digital PCR; the accuracy of the above-mentioned copy number assay results were verified through using methods including the real time PCR and Southern blot; the copy number of the Wx012(waxy gene) and SSII (starch synthesis gene) were detected by using the PINb-D1b as reference gene, and copy number of nib8 were detected by using the Wx012 and SSII as reference genes in the nib8 transgenic wheat, to verify the accuracy of the results using the PINb-D1b as reference gene; the specific primers were designed in different regions of the nib8 gene to compare the difference of the copy number analysis results of the nib8 transgenic wheat. Finally, a high-throughput assay method based on digital PCR was established to detect the copy number of target gene in wheat genome. 【Result】Finally, the detection system of target gene copy number in wheat genome based on digital PCR was determined. The optimal final concentration of primer and probe was 500 and 250 nmol·L-1, respectively. The optimal annealing temperature was 59 and the optimal amount of DNA template was 40 ng.℃PINb-D1b gene was used as the internal reference gene to detect the copy number of nib8 gene in the nib8 transgenic wheat. The results showed that the copy number of nib8 gene in the 12th, 16th, 17th, 23rd, 29th and 30th lines was 7, 1, 1, 1, 1, 1, 7, respectively.At the same time, the results were consistent with the real-time PCR and Southern blot results. In the nib8 transgenic lines, the PINb-D1b was used as the internal reference gene to detect copy number of the Wx012 and SSII, and we found that the results of copy number of nib8 gene in nib8 transgenic wheat by using these three kinds of internal reference gene analysis were consistent,which indicated that these three internal reference genes were all suitable for digital PCR method; the primers were designed in the upstream, middle and downstream regions of nib8 gene, and the results of copy number using different primer were consistent.【Conclusion】The method and reaction system to detect target gene copy number in wheat genome based on digital PCR method were optimized, and the detection system to detect target gene copy number in wheat genome based on digital PCR method was established. The results of digital PCR analysis are stable and reliable, and the detection efficiency is significantly improved, which has a certain application prospect.

Key words: wheat; copy number of genes; the digital PCR; the reference genes

开放科学（资源服务）标识码（OSID）：pagenumber_ebook=7,pagenumber_book=1931

收稿日期：2019-09-26；接受日期：2019-12-19

基金项目：国家重点研发计划（2016YFD0101802）、谷子氮素高效利用重要基因克隆及价值评估（2018ZX08009-17B）

联系方式：琚鹏举，E-mail：13610642179@163.com。

通信作者孙黛珍，E-mail：sdz64@126.com。

通信作者陈明，E-mail：chenming02@caas.cn

（责任编辑 李莉）