小麦单位面积穗数和粒长主效QTL紧密连锁KASP标记的开发及其效应评价

2021-07-30范涛李治蒋庆陈姝霖欧霞陈永艳任天恒

中国农业科学 2021年14期

范涛，李治，蒋庆，陈姝霖，欧霞，陈永艳，任天恒

范涛，李治，蒋庆，陈姝霖，欧霞，陈永艳，任天恒

四川农业大学农学院/四川省植物遗传和育种省级重点实验室，成都 611130

【】小麦单位面积穗数和籽粒粒长是小麦产量相关的重要农艺性状，对其进行遗传改良有利于提高小麦的产量。通过对前期QTL定位鉴定到的提高单位面积穗数的主效QTL位点和提高籽粒粒长的主效QTL位点开发相应的KASP分子标记，并在川农18和T1208构建的RILs群体中进行验证及评价，为更好地利用这两个QTL以及分子标记辅助育种奠定基础。利用前期在川农18和T1208构建的高代自交群体中鉴定到的控制小麦单位面积穗数主效QTL位点和控制籽粒粒长主效QTL位点，结合在这两个QTL区间内的55K SNP分子标记序列，开发设计KASP分子标记，并在亲本间筛选具有多态性的KASP分子标记。将筛选到的KASP分子标记在川农18×T1208的RILs群体中分别进行基因分型和鉴定相应表型性状的高低，并分析这两个主效QTL对于其他农艺性状的影响。和在亲本中具有多态性，和在群体中的验证表明这两个分子标记分别与和连锁。和能将群体材料的基因型分为2类，按照表型划分，在3年试验中，对多穗材料的平均选择率均达到72.58%，对少穗材料的平均选择率达到71.68%；对长粒材料的平均选择率达到69.86%，对短粒基因型的平均选择率可达61.52%，表明这两个标记的可靠性。基于KASP分子标记的基因分型结果表明，这两个QTL对于株高、千粒重、粒长、粒宽、粒径比、单位面积穗数、穗粒重均具有显著性影响。在川农17×川农11的RILs群体中进行验证也表明这两个分子标记对相应性状的选择具有一定的作用。针对单位面积穗数主效QTL位点和籽粒粒长主效QTL位点分别开发了1对与之连锁的KASP分子标记，可用于相应性状的选择，与KASP标记连锁的QTL分别能显著提高单位面积穗数和籽粒粒长。对株高、千粒重、粒长、粒宽、粒径比、穗粒重是负向影响，对株高、千粒重、粒宽、粒径比和穗粒重是正向影响，但对单位面积穗数是负向影响，这两个QTL及开发的KASP标记可应用于小麦高产育种中。

小麦；单位面积穗数；粒长；QTL；KASP；农艺性状

0 引言

【研究意义】普通小麦（L.）是世界上最重要的粮食作物之一，“养育”了世界范围内大约三分之一的人口。在所有的人类食物中，小麦提供的蛋白质和卡路里约占20%[1]。但是，根据预测，到2050年粮食的产量至少每年保持2.4%的增长率才能满足需求[2]。单位面积穗数、千粒重和每穗粒数是小麦主要的产量决定因素，因此，在育种中识别和引入相关的有利等位基因或QTL，对提高小麦产量具有重要意义。小麦的产量三要素单位面积穗数、穗粒数和粒重都是复杂的农艺性状，利用表型选择进行小麦的高产育种不仅费时费力，育种周期长从而导致育种成本提高。然而，与小麦产量单位面积穗数、穗粒数和千粒重相关的有利等位基因（或QTL）紧密连锁的分子标记的开发可以快速鉴定、识别相应的基因（或QTL），节约育种成本，提高育种效率。【前人研究进展】小麦籽粒粒重主要由籽粒形态决定，包括2个主要因素，即粒长和粒宽。目前分布在小麦所有染色体上的小麦籽粒形态以及千粒重相关的QTL已有诸多报道[3]。并且，与粒重或籽粒形态相关的、、、等基因也已经被克隆[4-7]。除此之外，还开发出许多与小麦粒重相关的功能基因的分子标记，并应用于育种上[8]。迄今为止，国内外研究学者利用不同的研究群体，将小麦单位面积穗数的相关的QTL定位在1A、1B、1D、2A、2B、2D、3A、3B、4A、4B、5A、5B、5D、6A、7A、7B、7D等多条染色体上，并且开发了较多的分子标记[3,9-12]。但是相应的分子标记大多具有品种特异性，因此，并未应用到实际育种中。随着测序技术的发展，利用分子标记辅助育种作为一种传统育种方法的补充工具，越来越受到育种家们的广泛关注。刘子会等[13]通过简单关联分析鉴定到一个与小麦耐热相关的分子标记。陈泠等[14]对2个抗穗发芽的分子标记和进行了有效性验证。LI等[15]针对小麦的穗颈长主效QTL开发了一个KASP分子标记，能够有效追踪该性状。MA等[16]开发了一个与粒长、粒宽和千粒重主效QTL紧密连锁的KASP分子标记，并且在重组自交群体中验证了它的有效性。胡洋山等[10]验证了2个分别位于3A和3B染色体上的SSR标记与小麦分蘖成穗呈现极显著正相关性，并证明它们可用于小麦分蘖成穗数的筛选和育种。【本研究切入点】利用SNP开发出的KASP标记筛选多穗和长粒小麦资源的研究仍鲜见报道。本文以分蘖成穗能力强的小麦品种川农18和具有籽粒性状优良的新品系材料T1208构建的371个RILs群体为材料，在前期研究中，通过QTL定位分析发现了一个与提高单位面积穗数的主效且稳定位于2D染色体上的QTL位点，在中国春参考基因组（IWGSC RefSeq v1.0）上的物理位置是82.19—95.90 Mb，此QTL在3个环境中被稳定检测到，平均LOD值为17.20，平均解释表型变异率达到19.04%，遗传效应来源于川农18；此外，3D染色体上定位到一个提高籽粒粒长的QTL位点，在中国春参考基因组（IWGSC RefSeq v1.0）上的物理位置是463.59—518.34 Mb，此QTL平均LOD值达到15.78，平均解释表型变异率为9.40%，遗传效应来源于亲本T1208，与以往研究对比，发现这两个QTL都是控制相应性状的新QTL位点[17]。【拟解决的关键问题】本研究依据QTL位点和区间内55K SNP基因芯片中的SNP分子标记设计KASP引物，并运用川农18×T1208杂交构建的RIL群体进行验证，评估这两个QTL对产量相关农艺性状的效应，以期为小麦分子标记辅助育种和多性状聚合育种提供理论依据。

1 材料与方法

1.1 小麦田间材料及种植

供试材料一共445份，其中包括由川农18×T1208杂交构建的371个家系的重组自交系（RILs）用于KASP标记的开发，由川农17×川农11杂交构建的74个RILs群体为验证群体。川农18×T1208杂交构建的RILs群体于2015—2016、2016—2017、2017—2018、2018—2019、2019—2020年（以下简称2015、2016、2017、2018和2019）种植于四川农业大学小麦育种基地（四川省邛崃市，30°25′N，103°28′E，海拔493.3 m），川农17×川农11构建的RILs群体于2019—2020年（以下简称2019）种植于相同地点。采用随机区组设计，小区行长2 m，每行10株，行距0.25 m，每个株系种植4行，3次重复，栽培管理同大田标准化管理方法，使用除草剂和杀菌剂，以保证生长过程中田间无虫害和病害，避免对相应农艺性状的影响。在以SNP基因分型和相应性状的表型匹配中，表型选择以当年群体材料中表型平均值为标准划分，川农18×T1208的RILs群体在2015、2016和2017单位面积穗数平均值分别是291.06、280.80和258.38穗/m2，当年群体中高于此表型值的定义为多穗株系，粒长在2016、2018和2019年的平均值分别是7.27、7.22和6.82 mm，当年群体中高于此表型值的定义为长粒株系。川农17×川农11杂交构建的RILs群体中单位面积穗数和粒长平均值分别是238.50穗/m2和6.90 mm，高于此表型值的材料分别定义为多穗株系和长粒株系。

1.2 小麦各农艺性状的调查

在植株成熟后，调查株高，株高的调查是从土壤表面到植株的顶端（不包含麦芒）；单位面积穗数调查则是每行调查5株，每个株系3个重复，最后取平均值，再乘以4得到每平方米单位面积穗数；收获时每个株系收获10穗，脱粒称取重量后，计算每穗粒重，公式为每穗粒重=十穗粒重/10。采用万深SC-G自动考种分析及千粒重分析仪（由杭州万深检测科技有限公司生产）测定千粒重、粒长和粒宽。

1.3 KASP引物设计及合成

根据构建的遗传连锁图谱上的信息[18]，分别查找位于内的SNP分子标记、、、和，以及位于内的SNP分子标记和，这些SNP标记的序列如表1所示，根据它们的SNP序列，使用在线网站平台PolyMarker（http://www.polymarker.info/）针对上述SNP分子标记分别进行KASP引物设计（表2）。每对KASP引物由3条DNA序列组成，分别是正向引物1、正向引物2和反向引物。另外每条正向引物1的5′端加上FAM荧光基团接头序列，序列为5′-GAAGGTGACC AAGTTCATGCT-3′，正向引物2的5′端加上HEX荧光基团接头序列，序列为5′-GAAGGTCGGAGTCA ACGGATT-3′。引物的合成由北京擎科生物科技有限公司完成。

1.4 KASP引物多态性筛选及群体中验证

对亲本川农18和T1208随机选取10粒种子，于培养皿中发芽培养至三叶期，根据李荣华等[19]方法并稍做改进提取基因组DNA，使用1%琼脂糖凝胶电泳检测DNA质量。利用上述设计引物进行PCR扩增。PCR反应体系为4.5 μL 1×KASP Master mixture（LGC Genomics, Hoddeston, UK）、2 μL（50 ng·μL-1）基因组DNA、2 μL KASP引物mix（正向引物1﹕正向引物2﹕反向引物﹕水=6﹕6﹕15﹕23）和1.5 μL ddH2O。PCR反应程序为95℃10 min；95℃20 s，61℃40 s，循环10次，每次循环退火延伸温度下降0.6℃；95℃20 s，55℃40 s，循环30次；25℃保温，采集荧光信号。PCR反应在荧光定量PCR仪进行（BioRad，CFX-96）。利用筛选出来具有多态性的KASP标记在川农18×T1208构建的371个株系的RILs群体进行基因分型，以及在川农17×川农11杂交构建的74个RILs群体进行验证，PCR的反应条件同上。

表1 试验中涉及的SNP标记及序列

表2 试验中涉及的KASP分子标记引物序列

下划线部分代表荧光接头序列 The underlined part represents the fluorescent junction sequence

1.5 数据统计分析

数据统计分析采用软件Microsoft Excel 2016和SPSS 22（IBM SPSS, Armonk，NY，USA）对农艺性状等数据进行方差分析、检验、相关性分析和描述性统计分析。参照李聪等[20]方法计算各农艺性状最佳线性无偏预测（best linear ubiased prediction，BLUP）值。利用CFX96荧光定量PCR仪的配套软件BioRad CFX Manager的基因分型模块，读取各样品不同荧光信号值，并分析基因型。基因分型结果与川农18相同的基因型定义为多穗基因型株系，与T1208相同的基因型定义为少穗基因型株系；基因分型结果与T1208相同的基因型定义为长粒基因型株系，与川农18相同的基因型定义为短粒基因型株系。

2 结果

2.1 川农18×T1208 RILs群体中单位面积穗数和粒长表型变异分析

在该群体中，单位面积穗数在2015、2016和2017年的平均值分别是291.06、280.80和258.38穗/m2，表型变异范围分别是185.3—442.7、179.3—387.3和120.0—442.7穗/m2；粒长在2016、2018和2019年的平均值分别是7.27、7.22和6.82 mm，表型变异范围是6.20—8.80、6.15—8.60和5.74—8.13 mm。单位面积穗数与粒长呈显著性负相关（=-0.44，＜0.001）。

2.2 KASP引物筛选结果

筛选区间内引物时发现，在亲本中具有多态性，该引物对亲本川农18分型结果是等位基因1类型，发出FAM荧光类型，T1208则是等位基因2类型，发出HEX荧光类型。其余设计的KASP引物不能很好地同时将2个亲本基因型识别出来（图1-a）。

筛选区间内引物时发现，在亲本中具有多态性，此引物对亲本川农18分型结果是等位基因1类型，发出FAM荧光类型，T1208则是等位基因2类型，发出HEX荧光类型。其余设计的KASP引物不能很好地同时将2个亲本基因型识别出来（图1-b）。

a：KASP-AX-111151907对亲本基因分型；b：KASP-AX-109962767对亲本基因分型

2.3 分子标记KASP-AX-111151907和KASP-AX- 109962767在RILs群体中分型结果

利用筛选到的多态性引物对上述371个川农18×T1208的RILs群体进行基因分型。该引物对群体中不同的材料能够很好的进行分型，各信号点荧光值较高，且不同基因型之间的夹角较大（图2-a），因此，分型效果较好。与亲本川农18相同的多穗基因型信号聚合在X轴，亲本T1208相同的少穗基因型信号聚合在Y轴。对371个株系的分型结果显示，有2个株系是杂合基因型，因此，在后续的分析中也剔除掉它们进行分析。引物鉴定到的369个株系中有158个株系是多穗基因型，这些多穗基因型中对应的多穗表型株系在2015、2016和2017年分别有99、122和123株，对应的多穗材料选择率分别是62.66%、77.22%和77.85%；剩下的211个株系是少穗基因型，其对应的少穗表型株系在2015、2016和2017年分别有148、151和155株，对应的少穗材料选择率分别是70.01%、71.56%和73.46%。在2015、2016和2017年分别有164、183和180株多穗表现型的材料，其中多穗基因型分别有100、183和180株，因此，这3年多穗材料中多穗基因型分别占60.98%、86.89%和88.33%。159个多穗基因型株系在2015、2016和2017年的表型平均值分别是311.60、303.41和290.50穗/m2，211个少穗基因型株系在2016、2018和2019年的平均值分别是275.56、263.75和234.36穗/m2。2种基因型材料在单位面积穗数上存在显著性差异（＜0.001），表明与紧密连锁（表3和图3-a）。

利用筛选到的多态性引物对川农18×T1208的RILs群体进行基因分型，结果与类似，对不同基因型材料的分型结果较好（图2-b）。该引物将与亲本T1208相同的长粒基因型信号聚合在Y轴，与亲本川农18相同的短粒基因型聚合在X轴。对371个株系基因分型结果显示，有8个株系是杂合基因型，因此，在后续的分析中剔除掉它们进行分析。鉴定到的363个株系中，有188个株系是长粒基因型，这些长粒基因型对应的长粒表型株系在2016、2018和2019年分别有127、132和135株，对应的长粒材料选择率分别是67.55%、70.21%和71.81%。剩下的175个株系是短粒基因型，其对应的短粒表型株系在2016、2018和2019年分别有104、114和105株，对应的短粒材料的选择率分别是59.43%、65.14%和60.00%。在2016、2018和2019年分别有199、194和206株长粒表型株系，而在这3年里长粒基因型分别有128、133和136株，因此，在这3年中长粒材料中长粒基因型分别占64.32%、68.56%和66.02%。188个长粒基因型株系在2016、2018和2019年的平均值分别是7.43、7.45和7.64 mm，175个短粒基因型株系在2016、2018和2019年的平均值分别是7.11、6.98和7.01 mm。2种基因型材料在粒长上存在显著性差异（＜0.001），表明与紧密连锁（表3和图3-b）。

a：KASP-AX-111151907在RIL群体中部分基因分型结果；b：KASP-AX-109962767在RIL群体中部分基因分型结果

表3 基于KASP分子标记验证QSn.sau-2D.2和QKl.sau-3D.2对相应性状的影响

括号中的“n”表示株系的数量；基因型A表示与川农18相同，基因型B表示与T1208相同

“n” in brackets indicates the number of lines. Genotype A is the same as Chuannong18 and genotype B is the same as T1208

+：含有相应QTL的株系，-：不含有相应QTL的株系；***：在p＜0.001水平差异极显著

另外，在这371个群体中，只有39个株系含有多穗基因型和长粒基因型，有117个株系只含有多穗基因型而不含有长粒基因型，有149株只含有长粒基因型而不含有多穗基因型，既不含有多穗基因型又不含有长粒基因型的株系58株。

2.4 基于KASP分子标记基因分型评价这两个主效QTL对其他产量相关性状的影响

基于各农艺性状的BLUP值，针对上述2个分别连锁的分子标记和与连锁的分子标记分型结果，利用检验分析这两个QTL对于其他农艺性状的影响（表4）。结果表明，对于株高、千粒重、粒长、粒宽、粒径比和每穗粒重都具有显著性降低的影响（＜0.001），对于株高、千粒重、粒宽、粒径比和每穗粒重有显著性提高的作用，但是对于单位面积穗数是降低的作用（＜0.001）。

2.5 分子标记KASP-AX-111151907和KASP-AX-109962767在其他群体中的检验

为了验证本研究中在川农18×T1208的RILs群体中开发的KASP标记是否可用于其他材料中小麦单位面积穗数和粒长的筛选，利用川农17×川农11的74个RILs群体进行基因分型（图4），和能够很好地将不同基因型的小麦材料进行区分。基因分型结果显示，有43个材料是多穗基因型，其中20个材料是多穗表现型，多穗材料选择率是46.51%，有31个材料为少穗基因型，其中少穗表现型材料有14个，对于少穗材料的选择率达到45.16%。基因分型结果显示，长粒基因型有38个，其中17个为长粒表现型，长粒的选择率达到44.74%，而短粒基因型有36个，其中短粒表现型有18个，对短粒基因型的选择率达到50.00%。表明这两个分子标记对川农17×川农11杂交群体中长粒和多穗的选择具有一定的作用。

3 讨论

3.1 协调型小麦

提高产量一直以来是众多育种学家们的基本目标。单位面积穗数、穗粒数和粒重作为小麦的产量三要素，与小麦的产量直接相关。四川盆地属于中国的西南麦区，是典型的雨养农业区，阴天较多，光照少导致小麦的分蘖成穗能力弱，从而导致单位面积穗数减少，进而造成产量下降[21-22]。在过去普遍的看法是四川盆地由于光照少和暖冬而存在“生态穗容量”，本课题组根据西南麦区的气候特征、土壤类型、光照条件及栽培模式等选育的多穗协调型小麦，其代表品种有川农12、川农17和川农18[23]。该类型的小麦品种能在保持较高穗重的前提下具有较高穗数，穗数及穗重相互补偿能力较强，为打破中国西南小麦生态区“生态穗容量”的限制提供了新的思路[24]。在前期的研究中也发现川农18控制单位面积穗数的主效QTL位于2D染色体上[17]，为阐明协调型小麦品种川农18成穗能力机制奠定了基础。本研究中开发的与川农18控制单位面积穗数的主效稳定QTL位点紧密连锁的KASP标记将有助于育种人员开展多穗小麦品种的选育工作。

图4 分子标记KASP-AX-111151907（a）和KASP-AX-109962767（b）在川农17×川农11 RIL群体中验证的部分结果

表4 基于KASP分子标记验证QSn.sau-2D.2和QKl.sau-3D.2对产量相关性状的影响

***：在＜0.001水平差异极显著，基因型A表示与亲本川农18相同的基因型，基因型B表示与亲本T1208相同的基因型

***: significant at the level of<0.001, genotype A is the same genotype as parent Chuannong 18, and genotype B is the same genotype as parent T1208

3.2 多穗和长粒KASP标记的利用

传统的杂交育种可以将不同的优良等位基因聚合在一起，获得理想性状，但是具有很大的局限性，对相应性状分子遗传机制不了解，育种周期长，并且需要花费大量时间和精力来鉴定表型，从而导致育种成本的提高，育种效率低下[25]。分子标记辅助育种如SSR标记、CAPS标记等已经广泛应用于小麦的育种工作中，张兆萍等[26]利用4个与穗发芽抗性相关的标记在350份不同小麦种质材料进行验证，最终筛选出2个分子标记可用于穗发芽抗性筛选，为筛选抗穗发芽小麦品种提供理论指导。Wang等[27]根据粒重相关基因的多态性序列开发了相应的CAPS标记，能够有效区分不同等位基因从而达到筛选具有更高粒重的小麦品种。SNP分子标记作为植物基因组中最丰富的一类分子标记类型，目前已经广泛应用于遗传图谱构建，遗传变异分析等[28]。但是其在分子标记辅助育种应用程度仍然不够普遍化[29]。本研究中开发的2个KASP分子标记，其中单位面积穗数相关的KASP分子标记在川农18×T1208杂交构建的RILs群体中能有效地对多穗基因型材料进行选择，在多年试验条件下，对多穗基因型材料的平均选择率可达到72.58%，对低穗基因型的材料平均选择率可达71.68%，而在多穗材料中，3年试验条件下，多穗基因型分别占了60.98%、86.89%和88.33%，而在川农17×川农11杂交构建的RIL群体中该标记对于多穗基因型的选择率可达46.51%，以上结果表明本研究开发的用于筛选多穗材料的KASP分子标记可靠性较高。值得注意的是，在本研究中单位面积穗数KASP分子标记对于多穗基因型材料的选择率分析是建立在稀播的基础上，以往的研究表明小麦品种川农18在大田生产密度下相比于对照成穗数显著性增加[23]，因此，在大田种植密度下该分子标记对于多穗基因型材料的选择率应该会有所增加。另外，本研究中开发的与粒长相关的KASP分子标记，在川农18×T1208杂交构建的RILs群体中对长粒基因型的平均选择率可达69.86%，对短粒基因型的平均选择率可达61.52%，在长粒材料中，3年试验条件下，长粒基因型分别占了64.32%、68.56%和66.02%，而在川农17×川农11杂交构建的RIL群体中该标记对于长粒基因型的选择率达44.74%，以上结果表明本研究开发的用于筛选长粒材料的KASP分子标记可靠性较高。本研究中多穗表型材料中多穗基因型占比和长粒表型材料中长粒基因型占比并未达到百分之百，可能的原因是因为这两个性状都属于数量性状由多基因控制，除了鉴定到的稳定主效QTL控制以外，还有许多微效基因共同控制相应的性状，而在之前的研究中也表明了单位面积穗数和粒长还有诸多微效QTL控制[17]。

基因聚合育种与传统聚合育种相结合已经成为了已经成为当今作物育种的新思路，并且已经取得了一定进展，如SALAMEH等[30]通过分子标记辅助育种选择将春小麦品系CM-82036的和这两个抗性主效基因导入到9个欧洲冬小麦品系中，结果发现当和聚合在一起时对赤霉病的抗性相比于单独更强。众所周知，由于营养竞争的关系，小麦的穗数增加，往往伴随着粒重的下降[31-33]。因此，如果将增加粒重的基因与增加穗数的基因聚合在一起，理论上小麦的产量将得到提高。在本研究中，川农18×T1208的RILs群体中，有39个株系同时含有增加穗数和粒长的基因型，但是它们具体的产量表现如何，还有待进一步的验证。

3.3 多穗和长粒KASP分子标记的其他用途

在小麦的产量相关性状中，各农艺性状之间往往具有显著的相关性，如小麦的千粒重与粒长、粒宽之间呈现显著性正相关，而与单位面积穗数、小穗数呈现出显著性负相关[34]。这往往归功于基因的“一因多效”，基因的“一因多效”在植物中广泛存在，如小麦的矮秆基因也多次被报道对于多个农艺性状具有影响，包括株高、千粒重、粒宽、单位面积穗数、旗叶大小等[17]。Lin等[35]定位到的2个控制小麦籽粒灌浆速度的QTL也同时对千粒重、粒长、粒宽、籽粒体积、籽粒面积、籽粒周长具有显著性影响。在本研究中，具有多穗基因型的株系在株高、千粒重、粒长、粒宽、粒径比和穗粒重方面都显著性低于少穗基因型的株系；具有长粒基因型的株系在株高、千粒重、粒宽、粒径比和穗粒重方面显著高于短粒基因型株系，在单位面积穗数上显著性低于短粒株系，这也表明本研究中开发的2个KASP分子标记对于农艺性状株高、千粒重、粒长、粒宽、粒径比、单位面积穗数和穗粒重的选择上都具有一定的作用。

4 结论

开发的与单位面积穗数相关的KASP分子标记和粒长相关的KASP分子标记能够分别区分川农18×T1208构建的RILs群体中的单位面积穗数和籽粒粒长的基因型，并有效鉴定单位面积穗数和粒长表型值的相对高低，可用于小麦高产育种。

[1] CHAVES M S, MARTINELLI J A, WESP-GUTERRES C, GRAICHEN F, BRAMMER S P, SCAGLIUSI S M, SILVA P R, WIETHOLTER P, TORRES G, LAU E Y, CONSOLI L, CAVES A. The importance for food security of maintaining rust resistance in wheat. Food Security, 2013, 5(2): 157-176.

[2] RAY D K, MUELLER N D, WEST P C, FOLEY J A. Yield trends are insufficient to double global crop production by 2050. PLoS One, 2013, 8(6): e66428.

[3] CAO S, XU D, HANIF M, XIA X, HE Z. Genetic architecture underpinning yield component traits in wheat. Theoretical and Applied Genetics, 2020, 133(6): 1811-1823.

[4] YANG J, ZHOU Y, WU Q, CHEN Y, ZHANG P, ZHANG Y, HU W, WANG X, ZHAO H, DONG L, HAN J, LIU Z, CAO T. Molecular characterization of a novelallele and its association with grain length in wheat (L.). Theoretical and Applied Genetics, 2019, 132(6): 1799-1814.

[5] SU Z, HAO C, WANG L, DONG Y, ZHANG X. Identification and development of a functional marker ofassociated with grain weight in bread wheat (L.). Theoretical and Applied Genetics, 2011, 122: 211-223.

[6] ZHANG P, HE Z, TIAN X, GAO F, XU D, LIU J, WEN W, FU L, LI G, SUI X, XIA X, WANG C, CAO S. Cloning ofassociated with grain weight in bread wheat and development of functional marker. Molecular Breeding, 2017, 37: 78.

[7] QIN L, HAO C, HOU J, WANG Y, LI T, WANG L, MA Z, ZHANG X. Homologous haplotypes, expression, genetic effects and geographic distribution of the wheat yield gene.BMC Plant Biology, 2014, 14(1): 1-36.

[8] 张福彦, 范家霖, 陈晓杰, 陈锋, 齐红志, 王嘉欢, 程仲杰, 杨保安, 张建伟. 小麦粒重相关基因的遗传定位和分子标记辅助育种进展. 植物遗传资源学报, 2020, 21(3): 507-516.

ZHANG F Y, FAN J L, CHEN X J, CHEN F, QI H Z, WANG J H, CHENG Z J, YANG B A, ZHANG J W. Genetic localization and marker assisted breeding of grain weight-related genes in common wheat. Journal of Plant Genetic Resources, 2020, 21(3): 507-516. (in Chinese)

[9] GAO F, WEN W, LIU J, RASHEED A, YIN G, XIA X, WU X, HE Z. Genome-wide linkage mapping of QTL for yield components, plant height and yield-related physiological traits in the Chinese wheat cross Zhou 8425B/Chinese Spring. Frontiers in Plant Science, 2015, 6: 1099.

[10] 胡洋山, 汤颖子, 李治, 晏本菊, 任正隆, 任天恒. 小麦分蘖成穗数相关分子标记在重组自交系(RIL)群体中的有效性验证及实用性评价. 麦类作物学报, 2018, 38(1): 8-15.

HU Y S, TANG Y Z, LI Z, YAN B J, REN Z L, REN T H. Evaluation and validation of molecular markers associated with maximum tiller number and spike number per unit area of wheat in a RIL population. Journal of Triticeae Crops, 2018, 38(1): 8-15. (in Chinese)

[11] FAN X, CUI F, JI J, ZHANG W, ZHAO X, LIU J, MENG D, TONG Y, WANG T, LI J. Dissection of pleiotropic QTL regions controlling wheat spike characteristics under different nitrogen treatments using traditional and conditional QTL mapping. Frontiers in Plant Science, 2019, 10: 187.

[12] MOHLER V, ALBRECHT T, CASTELL A, DIETHELM M, SCHWEIZER G, HARTL L. Considering causal genes in the genetic dissection of kernel traits in common wheat. Journal of Applied Genetics, 2016, 57: 467-476.

[13] 刘子会, 刘畅, 彭义峰, 李国良, 张华宁, 张红梅, 段硕楠, 郭秀林. 小麦耐热分子标记的鉴定及分析. 植物遗传资源学报, 2021, 22(1): 95-101.

LIU Z H, LIU C, PENG Y F, LI G L, ZHANG H N, ZHANG H M, DUAN S N, GUO X L. Identification and analysis of heat-tolerant molecular marker

[14] 陈泠, 高春保, 王翠, 朱展望, 佟汉文, 刘易科, 张宇庆, 邹娟, 何伟杰. 穗发芽抗性相关分子标记和在红白粒小麦中的有效性验证. 湖北农业科学, 2018, 57(24): 66-69.

CHEN L, GAO C B, WANG C, ZHU Z W, TONG H W, LIU Y K, ZHANG Y Q, ZOU J, HE W J. Validation ofandassociated with pre-harvest sprouting tolerance in red-grained and white-grained wheat varieties. Hubei Agricultural Sciences, 2018, 57(24): 66-69. (in Chinese)

[15] Li C, TANG H P, LUO W, ZHANG X M, MU Y, DENG M, LIU Y X, JIANG Q T, CHEN G Y, WANG J R, QI P F, PU Z E, JIANG Y F, WEI Y M, ZHENG Y L, LAN X J, MA J. A novel, validated, and plant height-independent QTL for spike extension length is associated with yield-related traits in wheat. Theoretical and Applied Genetics, 2020, 133: 3381-3393.

[16] MA J, ZHANG H, LI S, ZOU Y Y, LI T, LIU J J, DING P Y, MU Y, TANG H P, DENG M, LIU Y X, JIANG Q T, CHEN G Y, KANG H Y, LI W, PU Z E, WEI Y M, ZHENG Y L, LAN X J. Identification of quantitative trait loci for kernel traits in a wheat cultivar Chuannong16. BMC Genetics, 2019, 20(1): 77.

[17] REN T H, FAN T, CHEN S L, LI C S, CHEN Y Y, OU X, JIANG Q, REN Z L, TAN F Q, LUO P G, CHEN C, LI Z. Utilization of a Wheat55K SNP array-derived high-density genetic map for high-resolution mapping of quantitative trait loci for important kernel-related traits in common wheat. Theoretical and Applied Genetics, 2020, https://doi.org/10.1007/s00122-020-03732-8.

[18] REN T H, HU Y S, TANG Y Z, LI C S, YAN B J, REN Z L, TAN F Q, TANG Z X, FU S L, LI Z. Utilization of a wheat55K SNP Array for mapping of major QTL for memporal expression of the tiller number. Frontiers in Plant Science, 2018, 9: 333.

[19] 李荣华, 夏岩石, 刘顺枝, 孙莉丽, 郭培国, 缪绅裕, 陈健辉. 改进的CTAB提取植物DNA方法. 实验室研究与探索, 2009, 28(9): 14-16.

LI R H, XIA Y S, LIU S Z, SUN L L, GUO P G, MIAO S Y, CHEN J H. CTAB-improved method of DNA extraction in plant. Research and Exploration in Laboratory, 2009, 28(9): 14-16. (in Chinese)

[20] 李聪, 马建, 刘航, 丁浦洋, 杨聪聪, 张涵, 秦娜娜, 兰秀锦. 基于小麦55K SNP芯片检测小麦穗长和株高性状QTL. 麦类作物学报, 2019, 39(11): 1284-1292.

LI C, MA J, LIU H, DING F Y, YANG C C, ZHANG H, QIN N N, LAN X J. Detection of QTLs for spike length and plant height in wheat based on 55K SNP array. Journal of Triticeae Crops, 2019, 39(11): 1284-1292. (in Chinese)

[21] 任正隆. 雨养农业区的小麦育种. 北京: 科学出版社, 2011: 534.

REN Z L.. Beijing: Science Press, 2011: 534. (in Chinese)

[22] 任正隆. 中国南方小麦优质高效生产的若干问题. 四川农业大学学报, 2002, 20(3): 302.

REN Z L. Several limiting factors of wheat production in south area of china and the new approach of wheat breeding. Journal of Sichuan Agricultural University, 2002, 20(3): 302. (in Chinese)

[23] 谭飞泉, 张怀琼, 任正隆. “协调型”小麦新品种的产量潜力及其构成的研究. 四川农业大学学报, 2003, 21(3): 189-192.

TAN F Q, ZHANG H Q, REN Z L. Study on yield potential and its components of the "Coordination-type" wheat new cultivars. Journal of Sichuan Agricultural University, 2003, 21(3): 189-192. (in Chinese)

[24] 郝艳玲, 张紫晋, 粟永英, 张怀渝, 任正隆. 西南麦区高产多穗型小麦单株分蘖特征研究. 核农学报, 2016, 30(11): 2248-2257.

HAO Y L, ZHANG Z J, SU Y Y, ZHANG H Y, REN Z L. Morphological characteristics of tillers per plant in high-yield and multi-spike type wheat in southwest china. Journal of Nuclear Agricultural Sciences, 2016, 30(11): 2248-2257. (in Chinese)

[25] SUN C W, DONG Z D, ZHAO L, REN Y, ZHANG N, CHEN F. The Wheat 660K SNP array demonstrates great potential for marker-assisted selection in polyploid wheat. Plant Biotechnology Journal, 2020, 18(6): 1354-1360.

[26] 张兆萍, 周丽敏, 宋晓朋, 连俊方, 孙道杰. 小麦穗发芽抗性鉴定及相关分子标记的有效性验证. 麦类作物学报, 2015, 35(3): 300-305.

ZHANG Z P, ZHOU L M, SONG X P, LIAN J F, SUN D J. Identification of PHS tolerance in wheat varieties and validation of molecular markers associated with PHS tolerance. Journal of Triticeae Crops, 2015, 35(3): 300-305. (in Chinese)

[27] WANG S S, ZHANG X F, CHEN F, CUI D Q. A Single-Nucleotide polymorphism ofgene revealed its association with kernel weight in chinese bread wheat. Frontiers in Plant Science, 2015, 6(1166): 1166.

[28] CUI F, ZHANG N, FAN X L, ZHANG W, ZHAO C H, YANG L J, PAN R Q, CHEN M, HAN J, ZHAO X Q, JI J, TONG Y P, ZHANG H X, JIA J Z, ZHAO G Y, LI J M. Utilization of a wheat 660K SNP array-derived high-density genetic map for high-resolution mapping of a major QTL for kernel number. Scientific Reports, 2017, 7(1): 3788.

[29] 魏广辉, 李执, 陈强, 李阳, 陈诗豪, 裴英, 周勇, 程梦萍, 唐豪, 王际睿, 魏育明, 刘登才, 陈黎, 郑有良, 蒲至恩. 人工合成小麦SHW-L1高硒含量KASP分子标记开发及其应用. 中国农业科学, 2020, 53(20): 4103-4112.

WEI G H, LI Z, CHEN Q, LI Y, CHEN S H, PEI Y, ZHOU Y, CHENG M P, TANG H, WANG J R, WEI Y M, LIU D C, CHEN L, ZHENG Y L, PU Z E. Development and utilization of KASP marker for Se concentration in synthetic wheat SHW-L1. Scientia Agricultura Sinica, 2020, 53(20): 4103-4112. (in Chinese)

[30] SALAMEH A, BUERSTMAYR M, STEINER B, NEUMAYER A, LEMMENS M, BUERSTMAYRET H. Effects of introgression of two QTL for fusarium head blight resistance from Asian spring wheat by marker-assisted backcrossing into european winter wheat on fusarium head blight resistance, yield and quality traits. Molecular Breeding, 2011, 28(4): 485-494.

[31] GRIFFITHS S, WINGEN L, PIETRAGALLA J, GARCIA G, HASAN A, MIRALLES D, CALDERINI D F, ANKLESHWARIA J B, WAITE M L, SIMMONDS J, SNAPE J, REYNOLDS M. Genetic dissection of grain size and grain number trade-offs in CIMMYT wheat germplasm. PLoS One, 2015, 10(3): e0118847.

[32] PINTHUS M J, MILLET E. Interactions among number of spikelets, number of grains and grain weight in the spikes of wheat (L.). Annals of Botany, 1978, 42(4): 839-848.

[33] CALDERINI D F, SAVIN R, ABELEDO L G, REYNOLDS M P, SLAFER G A. The importance of the period immediately preceding anthesis for grain weight determination in wheat. Euphytica, 2001, 119(1/2): 199-204.

[34] LI F J, WEN W E, HE Z H, LIU J D, JIN H, CAO S H, GENG H W, YAN J, ZHANG P Z, WAN Y X, XIA X C. Genome-wide linkage mapping of yield-related traits in three Chinese bread wheat populations using high-density SNP markers.Theoretical and Applied Genetics, 2018, 131: 1-22.

[35] LIN Y, JIANG X J, TAO Y, YANG X L, WANG Z Q, WU F K, LIU S H, LI C X, DENG M, MA J, CHEN G D, WEI Y M, ZHENG Y L, LIU Y X. Identification and validation of stable quantitative trait loci for grain filling rate in common wheat (L.). Theoretical and Applied Genetics, 2020, 133: 2377-2385.

Development and Effect Evaluation of KASP Markers Closely Linked to Major QTLs of Spike Number Per Unit Area and Grain Length in Wheat

FAN Tao, LI Zhi, JIANG Qing, CHEN ShuLin, OU Xia, CHEN YongYan, REN TianHeng

College of Agronomy, Sichuan Agricultural University/Provincial Key Laboratory of Plant Genetics and Breeding of Sichuan Province, Chengdu 611130

】Spike numbers per unit area (SN) and kernel length (KL) are both important yield-related traits of wheat. Genetic improvement on SN and KL will help increase the yield of wheat. The corresponding KASP molecular markers were developed for theand the QTL, which were identified in previous study for increasing SN and KL, respectively. Both KASP markers were verified and evaluated in the RIL population crossed by Chuannong 18 and T1208, which laid a foundation for better utilization of these two QTLs.【】The major QTLfor SN and the major QTLfor KL were previously identified in the RIL population crossed by Chuannong 18 and T1208. The SNP molecular marker sequence of Wheat 55K SNP array within the two QTLs were used to develop and design the KASP molecular markers. The polymorphic KASP markers were screened with wheat parents and then verified in the RIL population. The selected KASP molecular markers were genotyped and used to identify the corresponding phenotypic traits in Chuannong 18×T1208 RILs population. The effects of the two major QTLs on other agronomic traits were analyzed.【】andhave polymorphic between the wheat parents. Subsequently, these two KASP markers were verified in the RIL population linkage withand, respectively.andcould divide the genotypes into two groups. According to the phenotype, the average selection rate offor multi spike lines was 72.58%, the average selection rate offor few spike lines was 71.68%, the average selection rate offor long kernel lines was 69.86%, and the average selection rate offor short kernel lines was 61.52%, indicating the reliability of the two KASP markers. Moreover, the genotyping results based on KASP molecular markers showed that the two QTLs had significant effects on plant height, 1000 kernel weight, kernel length, kernel width, kernel diameter ratio, spike numbers per unit area and kernel weight per spike, respectively. The validation in the RIL population of Chuannong17×Chuannong11 also indicated that these two KASP markers had a certain effect on the selection of the corresponding traits.【】This study developed two KASP markers for two major QTLand, respectively. Both KASP markers can be used for the selection of corresponding traits. These two QTLs which are tightly linked to these two KASP markers could significantly improve the SN and KL, respectively. In addition,had negative effects on plant height, 1000 kernel weight, kernel length, kernel width, kernel diameter ratio and kernel weight per spike. On the other hand,had positive effect on plant height, 1000 kernel weight, kernel width, kernel diameter ratio and kernel weight per spike, but had negative effects on SN. These two QTLs and the developed KASP markers can be used in the future high yield wheat breeding program.

wheat; spike numbers per unit area; kernel length; QTL; KASP; agronomic traits