水稻抗稻瘟病基因Pi2 的KASP 标记开发与应用

赵传超 ，丛 森，梁思怡，谢华斌，陆文宇，肖武名，黄 明，郭 涛，王加峰

(华南农业大学 国家植物航天育种工程技术研究中心/广东省植物分子育种重点实验室, 广东 广州 510642)

水稻Oryzasativa是世界上重要的粮食作物，是全球主要的食物来源之一，其产量的高低直接影响粮食安全。每年由细菌、真菌和病毒侵染造成的水稻减产高达30%～50%。其中稻瘟病对水稻生产造成的危害最为严重[1]，培育具有稻瘟病抗性的水稻品种是防治该病害最经济有效的方法[2]。然而，由于稻瘟病生理小种的迅速更替，许多抗病品种的抗病性在逐渐丧失[3]。目前已经被鉴定并克隆的稻瘟病抗性(R)基因有35 个[4]，其中Pi2基因来自对稻瘟病具有高度抗性的哥伦比亚籼稻品种5 173[5]，是一个广谱抗稻瘟病显性基因。该基因对13 个国家的36 个生理小种表现出高水平抗性[6]。针对我国华南籼稻区的稻瘟病优势生理小种，Pi2基因表现出良好的抗性[7-8]，利用Pi2基因选育具有广谱抗性的水稻品种具有重要的应用价值。

随着生物技术的发展，利用分子标记进行分子标记辅助选择(Maker-associated selection, MAS)可以大大提高育种的效率和进程[9]。目前针对Pi2基因开发并应用于育种的分子标记大多数属于第2 代分子标记，包括简单重复序列标记(Simple sequence repeats，SSR)、插入/缺失标记(Insert/Deletion，InDel)等，这类标记需要电泳技术进行检测，其操作过程较为繁琐，主要应用于基因定位和遗传图谱的绘制[10]，无法实现大批量的准确检测，严重妨碍了该基因在水稻育种中的应用。

第3 代分子标记-单核苷酸多态性(Single nucleotide polymorphism，SNP)标记[11-13]的出现为快速检测相关基因开辟了一条新的途径。SNP 是指分子水平上由单个核苷酸的改变即单个碱基的替换、插入或缺失，从而引起的DNA 序列的多态性，被认为是水稻基因组中最丰富的标记类型[14]。SNP在生物基因组中数量多，分布广泛，遗传稳定性高，与性状表达密切相关，富有代表性。另外，SNP 标记检测技术易于实现自动化高通量基因分型，可以有效缩短检测时间[15]。

近年来，竞争性等位基因聚合酶链式反应(Kompetitive allele specific PCR，KASP)技术已成为高效安全地检测SNP 的新方法。KASP 技术是一种基于SNP 标记的高通量基因分型技术，通过识别相应的荧光共振能量转移，利用基因探针对广泛的基因组DNA 样品，以及对目标SNPs 和InDels 进行精准的双等位基因分型[16]。该技术的检测过程无需电泳，可以有效地减少试验过程中有害物质的产生，现已广泛应用于各类作物的种质筛选工作[17-18]。针对Pi2基因，利用KASP 技术开发的分子标记相较于普遍使用SSR 标记具有稳定性好、检测效率高等优点，可以更加高效准确地鉴定抗病基因，为聚合多个优质基因、培育高产抗病水稻品种提供技术支持。

本研究通过对大量水稻种质资源的基因组测序数据分析，挖掘Pi2基因核心区域SNP 并开发为KASP 分子标记，进一步通过室内接种、病圃鉴定验证标记的可靠性，最终构建一套适合于Pi2基因的分型体系。研究结果将为利用稻瘟病抗性基因Pi2进行品种选育提供技术支撑。

1 材料与方法

1.1 供试材料

本研究以华南农业大学国家植物航天育种工程技术研究中心收集的593 份水稻自然种质作为供试材料，该群体分为122 个粳稻材料和471 个籼稻材料。通过该群体的基因组测序信息筛选SNP 位点，并利用H-74 和H-78 开发KASP 标记。以该群体田间2021 年早造和晚造材料以及阳江市2020 年晚造病圃材料作为验证材料。以R1179 作为抗病对照，‘日本晴’作为感病对照。材料采用水培法[19]种植，营养液参照国际水稻研究所(International Rice Research Institute，IRRI)的培养配方配制。

1.2 试验方法

1.2.1 DNA 提取方法 本研究采用改良的十六烷基三甲基溴化铵(Cetyltrimethylammonium bromide，CTAB)法提取DNA。具体流程如下：取少量幼嫩的叶片，放入2 m L 离心管中，加入2×CTAB 提取液800 mL，用高通量组织研磨仪(Geno/Grinder)以1 500 r/min 充分研磨3 min，65 ℃水浴30 min；冷却后加入三氯甲烷∶异戊醇体积比为24∶1 的混合液至满管(约900 μL)，混合均匀，8 000 r/min 离心10 min，放入–20 ℃的冰箱中20 min，使DNA 充分沉淀；10 000 r/min 离心3 min，去除上清液，倒置1 min，加入φ为70%的乙醇溶液800 μL，10 000 r/min 离心3 min，去除上清液，干燥后加入ddH2O100 μL，使DNA 充分溶解，置于–20 ℃保存。

1.2.2 KASP 标记开发 在美国国家生物技术信息中心(National Center for Biotechnology Information，NCBI)下载Pi2序列(Genbank 登录号：DQ352453)以及与Pi2互为等位基因的Pi9(Genbank 登录号：DQ285630.1) 和Pigm序列(Genbank 登录号：KU904633.2)。分别利用程序VISTA tools(https://genome.lbl.gov/vista/index.shtml) 与Softberry(http://www.softberry.com/berry.phtml)在线分析，获得序列组成和基因结构信息。对H-74 和H-78 材料进行单倍型分析，获取Pi2抗性位点核心区域和重要SNP 差异位点，并根据样品材料的相似性判断目标基因是否存在。结合单倍型分析结果获得的Pi2基因抗性位点核心区域的物理区间，利用程序Primer-BLAST(https://www.ncbi.nlm.nih.gov/tools) 和软件Premier5 设计普通PCR 引物。H-74 和H-78 这2 个材料各设置2 个重复，对目标物理位置进行扩增并测序。对上述获得的测序结果进行比对，分析差异SNP，优先选择单碱基突变、SNP 突变位点位于编码区且周边没有其他的差异SNP 位点，利用这些位点基于KASP 技术设计成若干KASP 分子标记。

1.2.3 KASP 基因分型 将待检测的DNA 样品稀释至50 mg/L，加入96 孔PCR 板中。配制KASP扩增反应体系20 μL(DNA 5 μL、引物Mix 0.4 μL、2×KASP Master Mix 4.6 μL，矿物油10 μL)。KASP扩增反应程序：95 ℃预变性180 s，循环1 次；95 ℃变性10 s，63.4 ～ 56.2 ℃退火/延伸30 s(–0.8 ℃/循环)，循环9 次；95 ℃变性10 s，55 ℃退火/延伸30 s，循环29 次。待PCR 扩增程序结束后，用荧光检测仪(FLUOstar Omega SNP)进行荧光检测，获得基因分型结果。

1.2.4 稻瘟病菌接种鉴定 使用开发的KASP 标记进行基因分型后，选择检测结果为阳性和阴性的材料各3 份，以及已知具有Pi2基因的材料R1179作为阳性对照和无Pi2基因的材料‘日本晴’作为阴性对照。将种子浸种催芽处理，待种子浸种露白后，以8×8 的规格种植于30 cm×20 cm 的盆中，每个品种种植4 株，重复2 次，待水稻生长至3～4 叶期，参照韩美丽等[20]的方法离体单菌接种稻瘟病菌GDYJ7，6～7 d 后调查。

自然诱发病圃设置在广东省阳江市阳江农业科学研究所白沙试验场，该病圃常年处于高温高湿的气候环境中，具有发病的地理环境条件，是广东省具有代表性的稻瘟病区。试验材料按单本种植，每份材料种5 行，每行6 株，行距23.3 cm，株距16.7 cm，四周种植诱发品种‘Co39’。病圃全生育期防虫不防病。在成熟期，以‘Co39’充分发病为准，对供试材料进行穗颈瘟调查，按IRRI 0～9 级的分级标准进行记录(0 级：无病；1 级：抗；3 级：中抗；5 级：中感；7 级：感；9 级：高感)。

2 结果与分析

2.1 Pi2 基因核心SNP 筛选

基于593 份自然群体的全基因组重测序结果，结合Pi2基因序列，利用MEGA X64 构建系统进化树，根据进化树分布密度，分为7 个部分(图1A)。从中选择2 个抗性不同且亲缘关系较远的材料H-74 和H-78，其中H-74 为抗病材料，与R1179 等抗病材料在III 区；H-78 为感病材料，与‘日本晴’、‘丽江新团黑谷’等感病材料在II 区。Pi2基因第2 外显子上第787～788 位密码子位置的核苷酸变异与稻瘟病抗性密切相关[21-24]。因此对H-74 和H-78 这2 个材料，在Pi2基因区段以及基因上下游100 kb 区域进行单倍型分析，确定Pi2基因抗性位点核心区域的物理位置在10 386 510～10 389 466 bp的范围内。单倍型分析结果中发现该区域内含有99 个重要的SNP 突变位点。根据材料测序结果与分析结果中位点突变情况的相似性判断H-74 材料含有Pi2基因，H-78 材料不含有Pi2基因。进一步利用593 份自然群体的全基因组重测序结果，分析区域范围内的所有差异S N P 位点，共筛选出175 个SNP 位点，并获得这些位点的碱基分布情况(图1B)。

图1 系统进化树及差异SNP 位点分析Fig. 1 Phylogenetic tree and analysis of different SNP sites

根据物理位置10 386 510～10 389 466 bp，利用软件Primer5 设计引物(表1)，分别利用3 对引物对H-74 和H-78 材料的目标片段进行扩增和测序，测序结果进行序列比对分析，以获得更加精确的SNP 位点及其分布和碱基突变情况(图1C，图2)。结合593 份材料中发现的175 个SNP 位点与单倍型分析获得的99 个SNP 位点，过滤数据，在10 386 510～10 389 466 bp 的区间内共筛选出28 个关键SNP 位点，其中9 个位于内含子，19 个位于外显子。

表1 测序引物Table 1 Sequencing primers

图2 Pi2 基因差异位点序列比对Fig. 2 Sequence alignment analysis of the different sites of Pi2 gene

2.2 KASP 标记开发

以不含有Pi2基因的‘日本晴’基因组作为参考序列，将H-74 材料与‘日本晴’基因组序列进行比对，对上述的28 个关键SNP 位点进一步筛选，剔除相同的突变SNP 位点。最后获得5 个关键的SNP 位点，分别在10 388 566、10 388 841、10 389 007、10 389 256 和10 389 341 bp 发生SNP 变异，利用这5 个SNP 位点设计KASP 引物(表2)。

表2 KASP 引物Table 2 KASP primers

以含有Pi2基因的水稻材料R1179 为阳性对照，不含Pi2基因的‘日本晴’为阴性对照，ddH2O 为空白对照，在同一批593 份自然群体材料中随机选取84 份，利用上述引物对其进行KASP基因分型检测。待测样品扩增产物的荧光信号显示为蓝色圆点，且靠近X轴，则表明该材料包含Pi2基因，即T/T 基因型；如果待测的扩增产物荧光信号显示为红色圆点，且靠近Y轴，则表明该材料不包含Pi2基因，即A/A 基因型。

如图3 所示，引物Pi2-C1、Pi2-C2 和Pi2-C3 的基因分型中有着分群聚集的趋势，引物Pi2-C4 和Pi2-C5 分型结果不够明确(图3)。为进一步确认这5 个引物的可行性，利用3 个已知包含Pi2基因的样品材料和3 个已知不包含Pi2基因的样品材料作为阴、阳对照，同样以ddH2O 作为空白对照，对其进行KASP 基因分型检测，结果如图4 所示。引物Pi2-C3 可以成功将阴、阳对照及空白对照和检测样品材料分开。基因分型结果和预期结果相匹配，初步判断引物Pi2-C3 的可行性。然而利用引物Pi2-C3 的分型结果较为分散，为获得较为直观和清晰的图片，在原有的检测PCR 程序基础上加4 个循环，获得新的KASP 基因分型结果(图5)。

图3 KASP 引物基因分型效果Fig. 3 Genotyping effects of KASP primers

图4 对照KASP 基因分型结果Fig. 4 KASP genotyping results for control samples

图5 利用引物Pi2-C3 增加PCR 循环数后KASP 基因分型结果Fig. 5 KASP genotyping results after increasing PCR cycles by primer Pi2-C3

2.3 KASP 标记鉴定效果分析

为评估KASP 标记Pi2-C3 的检测效果，本研究基于全基因组重测序结果选取了6 个基因型为T/T 的材料作为阳性对照、4 个基因型为A/A 的材料作为阴性对照，ddH2O 作为空白对照，对84 个田间水稻材料，分别利用设计的KASP 分子标记Pi2-C 3 和已报道的S S R 标记M-P i 2(P i 2-F：5′-CAGCGATGGTATGAGCACAA-3′，Pi2-R：5′-CGTTCCTATACTGCCACATCG-3′)进行检测，并对两者的检测结果进行比较。SSR 标记扩增后，含有Pi2基因的片段大小为450 bp，不含Pi2基因的片段大小为282 bp[25]。结果如图6A 所示，KASP 标记的基因分型结果中，有50 个材料的荧光信号点为蓝色，38 个材料的荧光信号点为红色，4 个材料的荧光信号点为绿色，前后进行2 次检测，分型结果一致，表明KASP 标记Pi2-C3 具有较好的重复性和稳定性。将上述结果与SSR 标记的电泳结果(图6B)相比较，两者有95%的相似性，并且在两者的结果中，10 个对照材料的检测结果与前期测序结果一致。因此可初步判断，KASP 标记Pi2-C3 可以成功区分纯合型T/T、A/A 和杂合型T/A 这3 种基因型。

图6 84 个田间材料Pi2-C3 KASP 分子标记的验证Fig. 6 Verification of Pi2-C3 KASP molecular marker for 84 field samples

从593 份材料中对该位点进行筛选，挑选出基因型为T/T 的2 份材料607-557 和RT1031-69；基因型为A/T 的22 份材料，包括京农大‘W030’、‘早粳占’、CT15672-3-2-3-2-2-M、‘华粳籼74’、白水糯、雪华新占、海南糯、1 641 选、远野4 号、‘黄广华占1 号’、黑糯47、Getu、海南黑、本地糯(阳春三甲)、黑节早、广粳籼165、观音稻、IRRI129、黄广软占、‘大丰糯’、‘鹏稻2 号’、‘凤新丝苗-1’；其余材料基因型为A/A。从基因型为A/A 的材料中随机挑选出60 个材料，再加上挑选出的基因型为T/T和A/A 的24 份材料，对这84 个材料进行KASP基因分型检测，进一步验证KASP 标记的可靠性。如图7A 所示，a-1 和a-2 对应基因型为T/T的2 个材料，与检测结果保持一致。a-3 至b-12 对应基因型为A/T 的22 个材料，基因型数量比T/T∶A/A＝13∶6，初步判断这22 个材料后代发生性状分离。c-1 至g-12 对应基因型为A/A 的60 个材料，54 个检测为A/A，6 个检测为T/T，准确率为93%(图7B)。

图7 593 份材料中部分材料KASP 基因分型结果Fig. 7 KASP genotyping results of partial materials of the 593 materials

2.4 抗稻瘟病表型鉴定

为进一步验证基因分型结果的准确性，以R1179与‘日本晴’分别作为阳性对照与阴性对照，对基因分型中确定含有Pi2基因的种质(H-11、H-12 和H-94) 与基因分型中确定不含有Pi2基因的种质(H-14、H-15 和H-151)进行稻瘟病离体打孔接种试验，5～7 d 后调查接种表型，根据叶片病斑大小分析材料的抗病性。结果(图8)表明，H-11、H-12 和H-94 对稻瘟病具有抗性而H-14、H-15 和H-151 对稻瘟病表现高感。表明抗性表型与基因分型结果一致，Pi2基因KASP 标记Pi2-C3 具有较高的可靠性。

图8 叶片接种表型Fig. 8 Leaf inoculation phenotype

利用开发的Pi2-C3 抗性标记对2020 年晚造中阳江市病圃的84 份材料进行KASP 基因分型，确定供试材料中Pi2基因的存在情况，同时进行田间病圃抗性评价(表3)。从阳江市病圃材料中，检测到携带Pi2抗性基因的材料46 份，占54.76%，未携带Pi2抗性基因的材料38 份，占45.24%。结合病圃表型观察，抗性材料68 份，包括高抗材料16 份、中抗材料25 份；感病材料16 份，包括高感材料5 份、中感材料5 份。在后续对病圃检测结果的验证中表明，携带Pi2基因的材料对稻瘟病表现为抗性，其中5 份材料表现为高抗，7 份表现为中抗。

表3 广东阳江病圃材料Pi2 基因检测结果及抗性评价Table 3 Pi2 gene detection results and resistance evaluation of materials from the Yangjiang disease nursery in Guangdong

3 讨论与结论

3.1 Pi2 基因分子标记的开发和应用

培育和种植抗稻瘟病水稻品种可以有效减少稻瘟病造成的损失，利用与抗稻瘟病基因连锁的分子标记进行辅助选择可以提高品种选育的效率。针对抗稻瘟病基因Pi2，先后有不同的研究人员开发了不同类型的分子标记用于该基因的分型。Wu等[26]利用Pi2基因的SSR 标记AP22，对277 份优质中国水稻亲本材料进行了检测；高利军等[25]基于Pi2基因与日本晴基因序列的差异开发共显性分子标记M-Pi2；李潜龙[27]基于Pi2基因1 543～1 545处3 个碱基的缺失，开发Indel 功能标记Pi2-1，使用功能标记Pi2-1 辅助选择育种改良文香1 号的稻瘟病抗性，抗性表型的选择准确率在90%以上；李冰等[28]利用Pi2位点Nbs4-Pi2 及其复等位基因Pi9、Pigm位点Pigm-R2 DNA 序列比对结果，开发了Indel 标记Pi2-CM1，该标记可有效区分Pi2与其等位基因Pi9和Pigm。除了传统的二代标记，已有研究人员针对Pi2基因开发了一系列SNP 标记。华丽霞等[23]开发的dCAPS 标记，在PCR 扩增后需要通过酶切和电泳进行检测，操作复杂，而金名捺等[24]开发的HRM 标记，对PCR 扩增体系要求较高，不适合进行大批量筛选。基于KASP 的基因分型技术无需电泳，稳定性好，可以有效解决以上问题，杨义强等[29]利用高抗稻瘟病品种‘黄广油占’与高感品种‘广陆矮4 号’Pi2基因第787 和788 密码子上的突变，开发了Pi2基因特异性KASP 标记W-Pi2，通过对40 个品种和部分F2分离群体的检测发现该标记可以实现对Pi2基因的准确分型。邓钊等[30]利用基于KASP 技术开发的针对16 个稻瘟病抗性基因的标记组合对隆两优和晶两优系列杂交稻品种及其亲本进行基因型检测，Pia、Pita、Pi2、Pi5和Piz基因在品种中的检出频率相对较高，均达到50%以上，其中Pi2在父本中的检出率达到了71.2%。

本研究在593 份自然群体材料中筛选出2 个表型不同、亲缘关系较远的材料H-74 和H-78，进一步对这2 个材料进行单倍型分析，获得Pi2基因抗性位点的核心区域和99 个重要SNP 位点。结合593 份材料的全基因组重测序结果，筛选出关键SNP 位点，最后成功开发出KASP 特异性分子标记Pi2-C3。该标记被应用于已知基因型的84 个自然群体材料的检测中，显示出较高的准确度。同时，从阳江市病圃材料中，选出84 份材料进行KASP基因分型和抗性评价，发现46 份材料携带Pi2抗性基因，占总体的54.76%，而38 份材料未携带Pi2抗性基因，占45.24%。本研究开发的KASP 标记专门针对Pi2基因进行检测。然而，在被检测的病圃材料中未携带Pi2抗性基因的材料可能由于携带Pi1、Pita及Piz-t等其他抗性基因而表现出稻瘟病抗性。与以往的研究相比，本研究开发的KASP 标记Pi2-C3 不需要进行电泳，准确度较高，并且易于实现高通量和自动化。另外，相较于其他KASP 分子标记，本研究基于Pi2抗性位点的核心区域进行标记的开发，确保了标记的真实性和可靠性。

3.2 KASP 基因分型技术的应用前景

目前，SSR 标记鉴定体系已被广泛应用，显著提高了品种鉴定的效率和准确性。然而，在水稻基因组中，具有多个等位基因的位点增加了构建DNA 指纹库的难度，需要频繁更新数据以避免转换错误。一些通量化的SSR 标记检测技术缺乏稳定性。相比之下，SNP 标记仅有3 个基因型，不受扩增片段大小的影响，可广泛用于DNA 指纹图谱的构建，推动了检测技术的发展。基于PCR 原理的KASP 基因分型技术具有直观的检测结果，适用于高通量检测。与同为SNP 检测技术的高分辨率熔解曲线(High resolution melting，HRM) 相比，KASP 技术成本更低，操作更简便。本研究开发的KASP 标记Pi2-C3 和SSR 标记M-Pi2 分别对包含阴阳对照的96 个材料进行检测，2 种标记检测结果高达95%的一致性。由于KASP 和SSR 标记是基于不同的原理和标记位点设计的，SSR 不是基于基因的标记，其多态性不能完全反映基因的存在，而KASP 是基于SNP 差异开发的，具有更高的准确性。此外，Pi2区域的某些位点可能在个别个体中呈现差异，这也可能导致KASP 和SSR 标记的检测结果不一致。另一方面，利用KASP 标记进行基因分型可以减少电泳点样的步骤，从而降低试验误差。此外， KASP 标记的应用避免了使用有毒化学物质如核酸染料，减少了对人体和环境的潜在危害。KASP 的结果更加直观，可以通过颜色区分不同基因型，直接定位到相应的材料。

3.3 结论

本研究针对稻瘟病抗性基因Pi2的核心区域，开发了基于KASP 基因分型技术的特异性分子标记Pi2-C3，并建立了一套水稻Pi2基因的KASP 基因分型体系。该标记能够准确检测水稻材料中Pi2基因，将不同水稻材料的Pi2位点分为T/T、T/A 和A/A 基因型，其中T/T 为抗病基因型。相较于其他分子标记，操作流程更加简便，检测结果更加高效。KASP 基因分型检测技术为水稻抗稻瘟病种质资源的鉴定和抗性群体的高通量筛选提供了技术支持，为水稻新时代抗性育种奠定了基础。