16份玉米自交系全生育期抗旱性鉴定与评价

李懿璞,苏治军,高聚林,于晓芳,王富贵,刘 剑,孙继颖,王志刚,胡树平,包海柱

(1.内蒙古农业大学农学院,内蒙古 呼和浩特 010018;2.内蒙古自治区作物栽培与遗传改良重点实验室,内蒙古 呼和浩特 010018;3.内蒙古农业大学职业技术学院,内蒙古 包头 014109)

玉米(ZeaMaysL.)是我国目前种植面积最大的粮食作物,在保障国家粮食安全方面发挥着重要作用。我国2/3的玉米主要种植在东北、华北、西北和西南等地区的丘陵或平原旱地上,由于自然降雨不足,每年因干旱造成玉米减产20%～30%,在旱情严重年份甚至可减产30%～35%[1-2]。干旱被认为是影响我国玉米总产量的首要非生物胁迫因子[3-4]。据中国气象局国家气候中心2021年预测,受大气环流异常的影响,未来10年我国发生区域性干旱的频率可能增加,并且干旱面积有扩大趋势。在全球水资源日益短缺的情况下,挖掘抗旱种质资源,选育高产抗旱品种已成为玉米育种的重要目标[5-10]。

国内外许多学者在鉴定玉米自交系抗旱性方面进行了大量研究,发现干旱胁迫对玉米自交系产量影响较大,而产量是玉米耐旱育种最为关注的性状,但仅通过田间产量评价玉米自交系的抗旱性并不科学;可运用简单相关、主成分分析、逐步回归、聚类分析和灰色关联分析方法,筛选出玉米苗期存活率、雌雄穗开花隔离时间、抽雄吐丝期叶片保水力、百粒重等作为玉米自交系抗旱性评价的主要指标[11-16],并利用这些指标筛选出一些抗旱种质材料。随着分子生物学的快速发展,抗旱QTL的定位与候选基因克隆也为玉米抗旱自交系的筛选带来了便利。近年来,一些与玉米抗旱性相关的基因相继被报道,如转录因子ZmDREB2.7[17]、蛋白磷酸酶ZmPP2C-A10[18]、转录因子ZmNAC111[19],与玉米苗期抗旱性显著相关、编码玉米液泡型氢离子焦磷酸酶的ZmVPP1基因也被成功克隆[20]。过去对玉米耐旱性研究的方法主要集中于大田鉴定和苗期胁迫,大田鉴定受到自然降水影响,试验条件难以控制;在实际农业生产中,花期干旱胁迫对玉米产量影响极大,仅通过苗期胁迫进行耐旱性鉴定所得的结果不够全面和准确。目前关于花期干旱胁迫对玉米植株形态及生理特性、灌浆、籽粒建成等相关指标影响的研究报道较少,系统的玉米全生育期抗旱性鉴定体系尚未形成。

本研究利用16份玉米自交系进行旱棚池栽试验,通过对自交系全生育期生长指标的监测和分析,明确抗旱自交系的筛选指标和科学的评价方法,为选育玉米抗旱品种提供理论依据和技术支持。

1 材料与方法

1.1 试验材料

供试材料为内蒙古玉米产区主推的7个品种的亲本和1个自育品种的亲本,共16份玉米自交系,分别为PH6WC、PH4CV、KWS49、KWS73、KW7M14、KW4M029、Mo17、E28、郑58(ZH58)、昌7-2(Chang 7-2)、DH382、DH351、Si-144、Si-287、G303、X178。

1.2 试验设计与处理

试验于2019—2020年春播期在内蒙古农业大学科技园区(包头市土默特右旗沟门镇北只图村,40°33′N,110°31′E)的旱棚中进行。棚内为正方形水泥池,边长2.2 m,深1.5 m。采用裂区试验设计,主区为水分控制处理,设置正常灌溉和干旱胁迫2个处理,正常灌溉处理为出苗后15～20 d浇水,以后每隔7 d灌溉1次,保持土壤含水量为田间最大持水量的75%左右;干旱胁迫处理为保持出苗后、开花前一周、灌浆后一周土壤达到中度干旱胁迫水平,即土壤含水量为田间最大持水量的45%左右,其他措施同正常灌溉处理。副区为16份自交系,单行区,行距60 cm,株距22 cm,种植密度75 000株·hm-2。试验设3次重复,总计96行,每池播种4行,共24个栽培池。

1.3 测定指标与方法

农艺与产量性状：参照《NY/T 2232-2012 植物新品种特异性、一致性和稳定性测试指南—玉米》标准中的方法,田间测定出苗率、抽雄期、吐丝期、散粉期、散粉吐丝间隔期(ASI)、成熟期、株高、穗位高;收获后测定穗长、穗粗、百粒重、穗行数、行粒数、单穗粒重、秃尖长。采取定株测量法,即每个自交系测10株(不包括边行植株)并取平均值。收获期收获各小区中间2行全部果穗进行称重,用谷物水分测量仪测定玉米籽粒含水量,按照标准含水量14%计算产量。

光合叶绿素指标：玉米散粉后10 d于晴天9∶00—11∶00采用Li-6400XT光合测定系统测定穗位叶的净光合速率;散粉后10 d采用Handy PEA植物效率分析仪测定穗位叶的荧光参数值Fv/Fm和Pi;散粉后10 d于晴天9∶00—11∶00和15∶00—17∶00用SPAD仪测定穗位叶的叶绿素含量(SPAD值),以上指标每个重复选择中间10株(不包括边株)进行测定。

土壤含水量：采用TDR法测定,测定土层为0～20、20～30、30～40、40～50 cm。

1.4 统计方法

1.4.1 抗旱系数(DC) 抗旱系数=干旱胁迫某性状测定值/正常灌溉该性状测定值[21]。

1.4.2 抗旱隶属度法(SV) 综合评价品种抗旱性选用抗旱隶属度法进行[22],公式如下：

(1)

式中,Xij为i品种j性状的抗旱隶属度值,X为i品种j性状值,Xjmin为各品种j性状最小值,Xjmax为各品种j性状最大值。当指标性状值与抗旱性呈负相关时,Xij= 1-Xij。

将各自交系每个生理指标抗旱隶属度值进行累加并求平均值,获得抗旱隶属度综合值。综合值越大,抗旱性越强。根据抗旱指标隶属度综合值,品种抗旱性强弱分为4级,即1级0.80～1.00,为强抗旱型(HR);2级0.60～0.79,为抗旱型(R);3级0.30～0.59,为中抗旱型(MR);4级0.0～0.29,为弱抗旱型(S)。

1.5 数据处理

采用Microsoft Excel 2010进行数据处理、抗旱系数和隶属度的计算;采用SPSS 22.0和SAS 9.0软件进行方差分析、相关性分析和逐步回归分析。所有结果均为两个生长季数据平均值。

2 结果与分析

2.1 不同玉米自交系各性状描述性统计及方差分析

对16份玉米自交系抗旱性状进行方差分析,结果表明(表1),自交系对ASI、穗长、穗粗、秃尖长、穗行数、行粒数、百粒重、单穗粒重、Pi值共9个性状的影响达显著或极显著水平。水分条件对ASI、株高、穗位高、穗长、穗粗、秃尖长、行粒数、单穗粒重、Pi、净光合速率和SPAD等共11个性状有显著或极显著影响,表明这些指标对水分胁迫较为敏感,干旱胁迫对其影响较大,可从中进一步筛选出自交系抗旱评价的主要指标。水分条件与自交系互作效应对各性状均无显著影响。干旱胁迫下株高、穗位高和单穗粒重标准差较大,说明这些性状在不同自交系间对干旱胁迫的敏感程度存在较大差异。

2.2 玉米自交系主要性状抗旱系数的相关性分析

采用性状相对值(抗旱系数)进行抗旱性评价,选取方差分析(表1)中水分条件的影响达显著水平的11个性状进行相关分析。由表2可知,穗长、单穗粒重、行粒数、株高与产量抗旱系数显著或极显著正相关;秃尖长、ASI与产量抗旱系数显著负相关。各主要性状与产量的相关性表现为ASI>单穗粒重>穗长>行粒数>株高>秃尖长。穗粗与单穗粒重、穗粗与行粒数、穗位高与株高等性状间也存在极显著相关关系,表明不同的信息存在交叉,其对玉米自交系的抗旱性也有一定的贡献,为更加明确各性状对玉米自交系抗旱性的贡献和性状间的内在联系,需要做进一步统计分析。

2.3 抗旱性相关指标的逐步回归分析

采用逐步回归法筛选重要的抗旱性指标,建立抗旱性评价数学模型,筛选出对因变量影响最大的自变量,可对玉米自交系抗旱性进行预测。产量是评价玉米抗旱性最重要的性状,故以产量抗旱系数为因变量、与产量抗旱系数相关性达到显著水平的穗长、秃尖长、单穗粒重、行粒数、ASI和株高等6个性状作为自变量进行逐步回归分析,由表3(见92页)可确定最优回归方程为：

Y=0.408+0.442X1+1.177X2+0.520X3(R2=0.623,F=12.097,P=0.004)

(2)

其中,X1、X2、X3分别为单穗粒重、ASI和穗长。因此,单穗粒重、ASI和穗长可作为玉米自交系抗旱性鉴定的次级性状筛选指标。

表3 回归系数的标准化检验Table 3 Standardized test of regression coefficients

2.4 不同玉米自交系抗旱性综合评价

采用两种方法进行抗旱性综合评价,即通过逐步回归筛选主要性状指标抗旱系数和计算旱情发生时期关键指标抗旱系数,以增强自交系抗旱性评价结果的可靠性。根据由逐步回归法确定的3个次级性状筛选指标(单穗粒重、ASI、穗长),基于3个性状抗旱系数计算各自交系的隶属度值,根据小区产量计算各自交系产量抗旱系数的隶属度值,通过平均隶属度评价自交系的抗旱性,按照模糊隶属函数法原理,隶属度越大,抗旱性越强。通过分析次级性状单穗粒重、ASI、穗长和产量抗旱系数的平均隶属度(表4,见92页),筛选出强抗旱型自交系5份,分别为Si-287、E28、Chang 7-2、PH6WC和Si-144,抗旱型玉米自交系3份,中抗旱型自交系7份,弱抗旱型1份(X178)。

表4 16份玉米自交系主要性状抗旱性评价结果Table 4 Drought resistance identified by subordinate function value method for 16 maize inbred lines

考虑玉米主产区旱情发生规律,为了检验逐步回归筛选的次级性状指标进行抗旱性综合评价的可靠性,分别对自交系苗期、开花期、收获期关键性状指标出苗率、ASI和产量的抗旱系数进行计算,综合评价自交系在不同生理时期的抗旱性,依据3个指标抗旱系数的平均隶属度进行分级评价(表5),共选出强抗旱型自交系4份,分别为Si-287、PH6WC、Chang 7-2和Si-144,抗旱型玉米自交系4份,中抗旱型自交系7份,弱抗旱型1份(X178)。与主要性状指标的抗旱性综合评价结果进行比较发现,完全一致的有11份,强抗旱型的4份与弱抗旱型的1份评级结果完全一致,抗旱型和中抗旱型有5份评级结果存在一定差异。

表5 16份玉米自交系不同生育时期重要性状抗旱性评价结果Table 5 Drought resistance identified by subordinate function value of vital traits of different growth stages for 16 maize inbred lines

基于两种评价方法的平均隶属度,利用平方欧氏距离法对16份参试自交系进行聚类分析,当欧氏距离等于4时,可以将参试自交系分成4类(图1)：第1类包括KWS49、KW7M14、DH382、G303、DH351和ZH58共6份抗旱型自交系;第2类包括1份弱抗旱型自交系(X178);第3类包括Si-287、PH6WC、Chang 7-2、Si-144和E28共5份强抗旱型自交系;第4类包括PH4CV、KWS73、Mo17和KW4M029共4份中抗旱型自交系。

图1 16份玉米自交系聚类分析Fig.1 Cluster analysis of 16 maize inbred lines

3 讨 论

3.1 自交系抗旱指标、鉴定综合评价

抗旱性是复杂的数量性状,了解玉米自交系的抗旱性是提高玉米抗旱性和选育抗旱品种的重要基础。要想全面确定自交系的抗旱性,应利用多指标综合评价。于永涛等[23]认为衡量不同玉米品种的抗旱性强弱时,产量是最重要的指标。浦军等[24]研究表明,ASI、保绿性的抗旱系数与产量抗旱系数呈极显著正相关,是田间鉴定玉米抗旱性的有效指标。Cairns等[25]认为,干旱胁迫条件下,开花前生物量与玉米自交系的产量密切相关[26]。师亚琴等[26]研究认为,干旱胁迫下,植株茎秆强度、花后干物质积累量、穗位叶SPAD值降低,且其抗旱系数与产量抗旱系数呈显著正相关,这3个性状可作为玉米自交系抗旱性鉴定的次级指标。张赛楠等[27]研究认为,玉米自交系的ASI、百粒重和穗位高可作为玉米抗旱性鉴定的次级性状筛选指标。本研究首先通过方差分析筛选出受水分条件影响显著的指标,并对各指标及产量抗旱系数进行相关分析,从中选择与产量相关性显著的指标进行逐步回归分析,建立较优回归方程,最终筛选出单穗粒重、ASI和穗长作为玉米自交系抗旱性鉴定的次级性状指标,与前人的研究结果基本一致。ASI在多个研究中被认为是鉴定玉米自交系抗旱性的有效指标,在筛选中需要高度关注。

抗旱系数反映了不同自交系对水分胁迫的敏感程度,可较为客观地评价不同自交系的抗旱性。隶属度可消除单个指标评价带来的片面性,将各自交系的抗旱指标扩展到[0,1]闭区间上,使得各自交系抗旱性的差异具有可比性。通过2 a的旱棚试验,采用抗旱系数结合隶属度综合评价法评价的16份自交系的抗旱性结果可靠。两种方法评级均为强抗旱型的Si-287和Si-144为‘吉单27’的亲本,而‘吉单27’的抗旱性和丰产性表现都非常优异[28];前人研究表明,骨干自交系‘昌7-2’具有较强的抗旱性[29],本研究与其评价结果一致。这3个自交系可作为优质抗旱种质资源创制的供体材料。

3.2 抗旱鉴定时期与方法

根据玉米的生长发育特征,可将其抗旱鉴定分为萌发期抗旱鉴定、苗期抗旱鉴定、开花期抗旱鉴定、灌浆期抗旱鉴定和进行整体综合评价的全生育期抗旱鉴定。苗期抗旱鉴定方法简单易行,可进行大批量材料的鉴定,但其不涉及产量性状,所以只适合抗旱指标的初步筛选和抗旱性状的QTL定位与基因克隆。全生育期抗旱鉴定主要采用旱棚、人工气候室鉴定和田间直接鉴定等方法,旱棚、人工气候室鉴定法下田间持水量可控,结果较可靠,但样本数量有限;田间直接鉴定法工作量大,且结果易受环境影响,但其简单易行,多用于玉米杂交种抗旱性的鉴定。本研究采用旱棚全生育期鉴定法,评价结果全面可靠,但受制于干旱池的面积,只选用了16个具有代表性的自交系进行抗旱性鉴定。后期利用分子标记辅助选育的杂交组合还需要通过田间直接鉴定法进行多年多点的抗旱性鉴定。

多指标抗旱性测定结果虽准确,但工作量大而繁杂,某些指标测定成本高。本研究同时对旱情发生时期的抗旱系数进行了综合评价,与筛选出次级性状进行抗旱性评价的强抗旱型和弱抗旱型结果基本一致,只有1份自交系结果存在差异,所以根据旱情发生时期抗旱系数进行评价的方法简单实用,适宜小型育种企业推广使用。

我国西北春玉米的不同主产区降水量差异大且同一地区降水量年际之间波动幅度大,玉米全生育期均有可能受到干旱胁迫。张文英等[30]和朴明鑫等[31]认为在玉米抗旱育种实际应用过程中可依据当地的气候条件和育种目标对抗旱指标进行取舍。本研究在后期试验中将结合当地旱情的发生情况对自交系及杂交种进行综合评价,其中对杂交种产量评估时将增加产量抗旱指数作为评价参数。分子标记辅助育种结合常规育种会大大加快育种进程,目前玉米已有多个抗旱基因被克隆,后期研究将结合与抗旱基因(ZmNAC11,ZmVPP1)共分离的分子标记对自交系的抗旱性进行评价,结合基因芯片对遗传背景进行选择,进一步提高抗旱性评价的可靠性和抗旱育种新材料创制的效率。另外,将对筛选出的强抗旱型和弱抗旱型自交系进行抗旱机理研究,以期揭示玉米抗旱的分子机制。

4 结 论

运用方差分析结合逐步回归分析的方法,对玉米全生育期11个性状的抗旱系数与产量抗旱系数的相关性进行分析,表明单穗粒重、散粉吐丝间隔期、穗长3个性状与抗旱性密切相关,可作为玉米自交系抗旱性鉴定的次级性状指标。基于这3个性状抗旱系数的平均隶属度和旱情发生时期抗旱系数的平均隶属度两种方法对16份材料进行聚类分析,发现5份强抗旱型玉米自交系(Si-287、PH6WC、昌7-2、Si-144、E28)、6份抗旱型玉米自交系(KWS49、KW7M14、DH382、G303、DH351、郑58)、4份中抗旱型自交系(PH4CV、KWS73、Mo17、KW4M029)和1份弱抗旱型自交系(X178)。

共同第一作者贡献说明：李懿璞和苏治军共同完成了本试验的实施和论文写作等过程,其中苏治军主要完成试验材料收集、试验设计与田间管理,李懿璞主要完成试验数据采集与分析、论文写作。