小麦种质资源成株期氮效率评价及筛选

赵瑞，张旭辉，张程炀，郭泾磊，汪妤，李红霞

西北农林科技大学农学院，陕西杨凌 712100

0 引言

【研究意义】小麦是中国重要的商品粮和战略储备粮之一，对国家粮食安全和供需平衡具有重要作用。氮素作为生命元素，是作物生长和产量形成的限制因子，对中国粮食产量的贡献占 35%[1]。目前，中国氮肥使用量约占世界总用量的 31.0%，氮肥平均利用率仅为 30.0%—41.0%，且化肥的流失量占施用量的40.0%[2-4]。随着氮肥的过量使用，小麦生产成本升高，氮肥利用率降低，导致氮淋失、N2O的排放等一系列环境污染问题[5]。因此，筛选氮高效种质资源，挖掘基因潜力，是提高小麦氮素利用效率的有效途径。选育氮高效小麦品种，减少氮肥投入量，维持产量稳定，对保障粮食安全具有重大意义。【前人研究进展】已有研究表明，作物对氮素的吸收、利用效率不同，耐氮胁迫的能力存在差异[6]。董鲁明[7]研究表明，玉米在氮素吸收和利用能力上存在基因型差异。杨睿等[8]对50个甘蓝型油菜进行氮效率基因型筛选，发现品种间氮效率相差 2.5倍。氮效率鉴定是作物氮高效育种工作的一个重要环节，目前，关于作物氮效率评价筛选方面的研究，在棉花（Gossypium spp.）[9-11]、油菜（Brassica napus L.）[12]、大麦（Hordeum vulgare L.）[13-16]、谷子（Setaria italica）[17-23]、水稻（Oryza sativa）[24]、玉米（Zea mays）[25]和小麦（Triticum aestivum）[26-31]等作物上已有报道。【本研究切入点】以上研究多采用苗期水培的筛选方法，以形态、生理性状为筛选指标，对作物种质进行氮效率鉴定，关于作物成株期的氮效率指标筛选及评价却鲜见报道。【拟解决的关键问题】本研究通过田间试验的方法，设置4个氮肥梯度，选取小麦成株期株高、穗长、旗叶长、旗叶宽、茎粗、可育小穗数、穗粒数、千粒重、粒长、粒宽和单穗产量为氮效率评价指标，利用模糊隶属函数法、主成分分析法和聚类分析法对108份小麦种质资源进行基因型差异分类和氮效率评价。旨在挖掘小麦农艺和产量性状的遗传潜力，筛选氮高效品种，为小麦氮高效育种提供种质资源。

1 材料与方法

1.1 试验材料

108份小麦种质资源，包括河南品种30份、河北品种21份、陕西品种9份、北京品种20份、山西品种9份、安徽品种2份、山东品种5份、甘肃品种3份以及引进品种9份（电子附表1）。供试材料均由中国农业大学农业生物技术国家重点实验室彭惠茹教授提供。

1.2 试验设计

试验分别于2018—2019和2019—2020年小麦生长季，在陕西省杨凌示范区西北农林科技大学标本区（108.24°E，34.20°N）进行种植。试验地土壤类型为壤土，播种前试验田0—25 cm土层土壤养分状况为全氮 0.94 g·kg-1、速效钾 175 mg·kg-1、速效磷 7.42 mg·kg-1、有机质 1.29 g·kg-1。试验采取裂区设计，以施氮水平为主因素，小麦品种为副因素。设置4个施氮处理，无氮处理 N1（0 kg·hm-2）、低氮处理N2（180 kg·hm-2）、正常氮处理 N3（270 kg·hm-2）和高氮处理N4（360 kg·hm-2），尿素（N 46.4%）、磷酸二铵（P2O546%，N 18%）和氯化钾（KCl 52%）以基肥施用。采取完全随机区组设计，每个小区种植108份品种，每个品种2行，行距25 cm，行长1.5 m，每行点播30粒。4个处理，3次重复，共计12个小区。其他按照小麦栽培要求进行田间管理。

1.3 性状调查与指标测定

收获前，每个重复小区中随机取10株参照《小麦种质资源描述规范和数据标准》[32]对株高（plant height，PH）、穗长（spike length，SL）、旗叶长（flag leaf length，FLL）、旗叶宽（flag leaf width，FLW）、茎粗（stem diameter，SD）（穗下节）、可育小穗数（fertile spikelet number，FSN）、穗粒数（grain number per spike，GNPS）和单穗产量（grain weight per spike，GWPS）等农艺性状进行考察；利用 SC-G型自动考种仪（中国 万深）测定千粒重（thousand kernel weight，TKW）、粒长（grain length，GL）和粒宽（grain width，GW）。

1.4 指标计算与数据分析

灰色关联系数 L0i(k)=（Δmin+ρΔmax）/（Δ0i(k)+ρΔmax）；

其中，Xi、CKi分别为低氮/高氮水平、正常氮水平下的性状测量值；PImin、PImax为各性状指标耐氮系数的最小值、最大值；ri为各性状指标与综合耐氮系数的相关系数；为指数权重，表示第i个性状指标在所有指标中的重要程度。Δ0i(k)为k时刻2个序列的绝对差；Δmin和Δmax分别为所有比较序列各个时刻绝对差中的最小值与最大值，Δmin=0，分辨系数ρ=0.5。

运用Excel 2010和SPSS分析软件进行生物统计、方差分析、主成分分析、聚类分析及相关性分析，R语言的GGE-Bioplot软件包作GGE双标图。

2 结果

2.1 不同氮素水平下小麦成株期氮效率相关性状的差异分析

与N3处理相比，N1处理对小麦的形态和产量构成相关性状均有不同程度的影响（表 1）。主要表现为株高降低、穗长变短、旗叶叶长和叶宽缩小、茎秆变细、可育小穗数和穗粒数减少以及单穗产量降低。下降程度依次为旗叶长（12.39%）＞穗长（8.68%）＞旗叶宽（7.28%）＞可育小穗数（5.92%）＞穗粒数（4.94%）＞株高（3.96%）＞单穗产量（3.63%）＞茎粗重（1.77%）。数据表明大多数性状在小麦生长发育过程中均受到了氮胁迫的抑制。与N2处理相比，N4处理同样对供试材料各性状产生显著影响，处理间和种间的差异均达显著水平。随着氮肥施用量的增加，株高、穗长、旗叶长、旗叶宽、可育小穗数和穗粒数等性状值随之增加，上升程度依次为旗叶宽（7.55%）＞旗叶长（5.63%）＞穗长（4.28%）＞可育小穗数（2.44%）。其他性状值（茎粗、千粒重、粒长、粒宽和单穗产量）随之降低，下降程度依次为茎粗（2.13%）＞千粒重（2.01%）＞单穗产量（1.32%）＞粒宽（1.24%）＞粒长（0.28%）。说明高氮胁迫造成了小麦旺长，营养器官变大，影响植株碳氮代谢协调，植株籽粒产量降低。

表1 不同氮水平下小麦成株期氮效率相关性状的差异Table 1 Differences of nitrogen efficiency-related of wheat traits at mature stage under different nitrogen levels

108份参试小麦品种在4个供氮水平下各性状的变异幅度都比较大，变异系数大多超过10%。变异系数越大，说明品种对氮素敏感程度越大，2年4个氮梯度环境下的各性状指标平均变异系数（氮素敏感度）依次为单穗产量（20.39%）＞旗叶长（18.10%）＞穗粒数（17.30%）＞旗叶宽（17.00%）＞株高（14.52%）＞穗长（12.25%）＞茎粗（10.86）＞千粒重（10.47）＞可育小穗数（9.75%）＞粒长（5.03%）＞粒宽（4.39%）。结果表明，氮处理效果良好，各性状表型明显。单穗产量、旗叶长、穗粒数、旗叶宽、株高、穗长、茎粗和千粒重的变异系数均超过10%，可以作为小麦耐低氮和高氮能力的主要筛选与评价指标。供试材料的氮效率相关性状具有丰富、广泛的遗传性差异，可以用于小麦氮效率品种的筛选。

2.2 耐氮系数的区间分布和相关性分析

为减小品种间遗传及生物学差异造成的误差，采用耐低氮系数PI1（N1/N3）和耐高氮系数PI2（N4/N2）表示性状对低氮和高氮胁迫的敏感程度。将各性状指标的耐低氮系数PI1和耐高氮系数PI2分别以0.5为组距，分成 3个区间，统计数值的分布频数情况（表2），不同性状指标的PI1和PI2在同一区间或不同区间的分布频率差别较大。2年数据结果显示，在低氮胁迫条件下，穗长耐低氮系数0≤PI1<1区间的品种占71.30%和78.10%、旗叶长0≤PI1<1区间的品种占73.15%和66.67%、茎粗0≤PI1<1区间的品种占69.44%和70.48%、旗叶宽0≤PI1<1区间的品种占77.78%和54.29%、穗粒数0≤PI1<1区间的品种占63.89%和63.81%、单穗产量0≤PI1<1区间的品种占67.31%和59.05%、株高0≤PI1<1区间的品种占65.74%和61.90%、可育小穗数0≤PI1<1区间的品种占50.93%和71.43%、说明大多数小麦的穗长、旗叶长、茎粗、旗叶宽、穗粒数、单穗产量、株高和可育小穗数在低氮胁迫后降低，且对低氮胁迫响应明显。在高氮胁迫条件下，穗长耐高氮系数PI2≥1区间的品种占74.07%和67.31%、旗叶长PI2≥1区间的品种占68.52%和62.50%、可育小穗数PI2≥1区间的品种占61.11%和58.65%，说明穗长、旗叶长和可育小穗数对高氮胁迫响应敏感，耐肥性较强。

表2 小麦种质各性状指标耐氮系数PI的区间分布Table 2 The interval distribution of nitrogen tolerance coefficient PI of each trait in wheat germplasm

各性状指标耐低氮系数相关性分析显示，多数性状间的相关性达到显著或极显著水平（图1-A），株高与穗长、旗叶长、可育小穗数、穗粒数、千粒重、粒宽和单穗产量极显著相关，其中，与旗叶长的相关性达 0.62。穗长与旗叶长、可育小穗数、穗粒数、茎粗、粒长和单穗产量极显著相关，与千粒重显著相关。旗叶长与可育小穗数、穗粒数、茎粗、粒长和单穗产量极显著相关，与千粒重和粒宽显著相关。旗叶宽与千粒重、粒长和粒宽相关，与其他不相关。可育小穗数与穗粒数、茎粗和单穗产量极显著相关，与其他不相关。穗粒数与茎粗、千粒重和单穗产量极显著相关，与粒宽显著相关。茎粗与千粒重显著相关，与粒长负相关。千粒重与单穗产量极显著相关，与粒长负相关。粒长与粒宽极显著相关，与单穗产量显著相关。粒宽与单穗产量极显著相关。除了旗叶宽之外，各性状指标与单穗产量均显著相关，其中穗粒数和千粒重与其相关性最大，分别为0.81和0.79。

各性状指标耐高氮系数相关性分析显示（图1-B），株高与旗叶长和茎粗极显著相关，与可育小穗数、千粒重、粒长和粒宽负相关。穗长与旗叶长极显著相关，与旗叶宽和可育小穗数显著相关。旗叶长与旗叶宽和穗粒数呈显著相关，与茎粗极显著相关。旗叶宽与可育小穗数、穗粒数和单穗产量极显著相关，与千粒重、粒长和粒宽显著相关。可育小穗数与穗粒数和单穗产量极显著相关，与千粒重、粒宽和粒长显著相关。穗粒数与茎粗和单穗产量极显著相关。茎粗与单穗产量极显著相关，与粒长呈负相关。千粒重与粒长、粒宽和单穗产量极显著相关。粒长与粒宽和单穗产量极显著相关。粒宽与单穗产量极显著相关。在高氮胁迫下，除了株高和穗长之外，其他指标均与单穗产量呈显著相关或极显著相关。相关性分析说明各性状指标间均存在不同程度的相关性，且提供的耐氮信息交叉重叠。

2.3 小麦各性状指标的主成分分析

采用各性状指标的耐低氮系数PI-L（2年PI1的平均值）进行主成分分析（表3），提取4个主成分，贡献率分别为 39.766%、16.661%、9.361%和9.275%，累积贡献率达75.064%。第1主成分贡献率最大，主要由单穗产量、可育小穗数和穗长决定，此主成分主要反映低氮胁迫下小麦穗部性状的适应性变化。第 2主成分主要由粒长、粒宽和千粒重决定，其中，粒宽的特征向量最大（0.660），茎粗的特征向量最小（-0.383），说明小麦通过茎秆的变细来增加籽粒的饱满度和容重从而适应低氮环境。第3主成分中旗叶宽的特征向量最大（0.712），株高的特征向量最小（-0.425），说明小麦通过降低植株的高度来增加旗叶宽，此主成分反映小麦叶片的适应性变化。第4主成分中株高的特征向量最大（0.465），单穗产量的特征向量最小（-0.345），说明小麦通过降低单穗产量来保证小麦的正常生长，此主成分反映小麦植株的整体适应性变化。

表3 小麦氮效率相关性状的主成分分析Table 3 Principal component analysis of nitrogen efficiency related wheat traits

对各性状指标的耐高氮系数PI-G（2年PI2的平均值）进行主成分分析，将原来各单项性状指标转换成 4个新的相互独立的综合指标，累积贡献率达74.037%。主成分 PC1的特征值为 3.448，贡献率为31.348%，主要由单穗产量和千粒重决定，此主成分主要反映高氮胁迫下小麦籽粒性状的适应性变化。主成分PC2的特征值为2.243，贡献率为20.387%，其中旗叶长的特征向量最大（0.742），千粒重的特征向量最小（-0.490），说明高氮环境造成了小麦旗叶长度的增加，同时导致千粒重的降低。主成分PC3的特征值为1.370，贡献率为12.452%，其中穗粒数的特征向量最大（0.506），穗长的特征向量最小（-0.413），说明小麦通过减小穗长来增加穗粒数，此主成分反映小麦穗部的适应性变化。

2.4 小麦种质的氮效率综合评价

2.4.1 耐低氮性评价 利用性状指标的耐低氮系数，采用模糊隶属函数法计算隶属函数值U(X)，得出各性状指标的权重。2019年株高、穗长、旗叶长、旗叶宽、可育小穗数、穗粒数、茎粗、千粒重、粒长、粒宽和单穗产量的权重分别为0.083、0.120、0.100、0.115、0.111、0.114、0.077、0.072、0.026、0.050 和 0.131。2020年为 0.046、0.073、0.107、0.058、0.051、0.117、0.071、0.145、0.058、0.116和0.162。分别计算出2019年和 2020年每个小麦品种的耐低氮性综合评价值（D），D值越大表明其耐低氮能力越强。以2年的耐低氮性综合评价值D为评价指标，采用欧氏最长距离法对108份小麦品种进行聚类分析，将其分为5个类群（图2，电子附表1）。第5类为低氮敏感型，包括中种麦10号、泛麦5号、冀麦585、郑农46、农大3753、周麦25、洛麦23、兰考矮早8、晋47等31个品种，占29.25%；第4类为低氮较敏感型，包括烟农23、周麦24、洛旱7号、农大211、兰考926、京冬8号、陕麦 139和中麦 175等 21个品种，占总数的19.81%；第3类为中间型，包括许科316、育德1号、轮选 987、中旱 110、郑麦 9962、临汾 138、洛夫林10号、石麦12、西农585、新科8号、西农979和山农0917等47个品种，占44.34 %；第2类为耐低氮型，包括温麦19、西农529、石4185、陇麦212和丰抗2号5个品种，占4.72%；第1类为强耐低氮型，包括中麦875和西农158，占1.89%。

2.4.2 耐高氮性评价 利用性状指标的耐高氮系数，采用模糊隶属函数法计算出隶属函数值U(X)，得出各性状指标的权重。2019年株高、穗长、旗叶长、旗叶宽、可育小穗数、穗粒数、茎粗、千粒重、粒长、粒宽和单穗产量的权重分别为 0.044、0.103、0.114、0.113、0.119、0.135、0.083、0.058、0.046、0.047 和0.138。2020年分别为 0.055、0.089、0.086、0.108、0.091、0.134、0.070、0.080、0.048、0.084 和 0.156。以2年耐高氮性综合评价值D为依据，采用欧氏最长距离法对108份小麦品种进行聚类分析，将其分为4个类群（图3，电子附表1）。第1类为耐高氮型，包括兰考矮早8、CA9553、中原18、尧麦16、良星99、农大179、豫农9901、兰考926和郑农46等9个品种，占 8.65%；第 2类为中间型，包括陕麦 139、洛旱 7号、石4185、陇麦212、育德1号、西农979、许科316、郑麦366、山农0917、金丰7183等53个品种，占50.96%。第3类为高氮较敏感型，包括许农5号、偃展4110、农大211、河麦58-3、京冬8号、农大3492、科农2009、京411和CA9722等13个品种，占总数的12.50%；第4类为高氮敏感型，包括西农585、温麦19、轮选987、西农158、农大3753、西农529、中麦415、烟农23、中旱110等29个品种，占27.88%。

2.5 氮效率相关性状与综合评价D值的相关性分析

2.5.1 灰色关联度分析 灰色关联度分析显示（表4），在低氮胁迫条件下，各性状指标与综合评价D值的关联度最小值为0.748，最大值为0.815。关联程度由大到小依次为单穗产量>穗长>千粒重>旗叶长>株高>粒宽>可育小穗数>粒长>茎粗>穗粒数>旗叶宽，结果显示，与小麦品种耐低氮性关系最密切的指标性状为单穗产量、穗长、千粒重、旗叶长、株高。在高氮胁迫条件下，各性状指标与综合评价值的关联度最小值为0.663，最大值为0.765。关联程度由大到小依次为单穗产量>旗叶长>千粒重>粒宽>茎粗>穗粒数>株高>穗长>旗叶宽>粒长>可育小穗数，结果表明，与小麦品种耐高氮性关系最密切的指标性状为单穗产量、旗叶长、千粒重、粒宽和茎粗。

表4 各性状指标耐氮系数与综合评价D值的灰色关联度分析Table 4 Grey relational analysis between nitrogen tolerance coefficient and comprehensive evaluation value D

2.5.2 GGE双标图分析 利用GGE双标图对小麦品种各性状指标耐氮系数与综合评价D值进行相关性分析，通过品种和性状的相互关系对不同小麦品种的耐氮性进行分组，并揭示组内耐氮能力最强的品种。将同一方向上距离原点最远的品种连接构成多边形，多边形各边的垂线把双标图分成若干个扇形区，不同的小麦品种包括相关性状分布在不同的扇区内。每个扇区内“顶角”的品种为所包含性状表现最好的品种，位于多边形内部或者靠近原点的品种则为表现较差的品种。结果如图4-A所示，以耐低氮品种新麦19（93）为代表的扇区内，包含温9519（89）、石4185（73）、陇麦212（83）和郑麦9962（61）等耐低氮性强的品种，耐低氮性状为穗长、株高、穗粒数、茎粗、粒长和粒宽。在耐低氮品种西农 158（84）为代表的扇区内，有中麦 875（82）、丰抗 2号（78）、许科 316（94）和育德1号（97）等耐低氮性强的品种，耐低氮性状为千粒重和小穗数。以低氮敏感型品种泛麦 5号（52）为代表的扇区内没有耐低氮性状。由图 4-B可知，以耐高氮品种良星99（39）为代表的扇区内，包含 CA9553（49）和许农5号（95）等品种，耐高氮性状为穗长和株高。在耐高氮品种兰考矮早8（105）为代表的扇区内，有兰考926（104）、郑农46（100）、郑麦9962（61）和尧麦16（67）等品种，耐高氮性状为粒长、粒宽千粒重和穗粒数。以敏感型品种石麦12（74）为代表的扇区内没有耐低氮性状。结果表明，不同基因型小麦氮素吸收利用效率存在差异，且耐氮机制也不相同。

2.6 籽粒产量和氮效率分类

2.6.1 不同氮水平处理对小麦籽粒产量的影响 比较不同施氮水平下小麦品种的单穗产量（图 5）。结果表明，在低氮胁迫下，2年单穗产量均降低的品种有48份，其中降低幅度小于5.00%的品种有3份，刑531、洛旱7号和石新633。降幅超过20.00%的品种有郑麦98、矮杆早、新乡大穗、冀729、丰舞981、金丰7183、京411、泛麦5号、京冬17、CA9553和YX11-57。降幅超过 30.00%的品种有中种麦 10号、益科麦5号、剡麦98和河麦58-3。2年单穗产量均增加的品种有16份，其中，升高幅度超过20.00%的品种有西农585、石4185、温麦19、农大3492和陇麦212，升幅超过30.00%的品种有育德1号、许科316、西农158和丰抗2号。

在高氮胁迫下，2年单穗产量均降低的品种有26份，其中，降低幅度小于5.00%的品种有3份，中鉴49、刑531和山农0917。降幅超过20.00%的品种有西农158、新乡大穗、临抗22、周麦25、豫农201和洛麦23。2年单穗产量均增加的品种有20份，其中，升高幅度超过20.00%的品种包括长4738、CA9553、农大211、剡麦98和长丰1号，增幅超过30.00%的品种有兰考矮早8和许农5号。

2.6.2 基于单穗籽粒产量的氮效率分类 作物在氮胁迫下的产量常被用作氮效率评价指标。以耐低氮性综合评价D值和低氮处理下的单穗籽粒产量进行散点作图，将供试小麦品种分为4类（图6-A）。高产耐低氮型（Ⅰ区），此类品种的单穗籽粒产量和D值均高于供试品种的平均值，包括新科8号（77）、西农158（84）、陇麦 212（83）、温麦 19（89）、新麦 19（93）、西农529（90）、西农979（92）、许科316（94）和育德1号（97）等；高产低氮敏感型（Ⅱ区），此类品种的单穗籽粒产量高于平均值，D值低于平均值。包括平凉44号（28）、陇辐 2号（29）、L54223（32）和剡麦 98（31）等；低产低氮敏感型（Ⅲ区），此类品种的单穗籽粒产量和D值均低于平均值。包括中种麦10号（42）、新乡大穗（44）、金丰 7183（55）、衡 4399（45）、河麦58-3（46）和益科麦5号（58）等；低产耐低氮型（Ⅳ区），此类品种的单穗籽粒产量低于平均值，而D值高于平均值，包括丰抗2号（78）、轮选987（79）、农大3492（21）和山农0917（41）等。

以耐高氮综合评价D值和高氮处理下的单穗籽粒产量进行散点作图，将供试小麦品种分为 4类（图6-B）。高产耐高氮型（Ⅰ区），此类品种的单穗籽粒产量和D值均高于供试品种的平均值，包括兰考矮早8（105）、中原18（102）、许农5号（95）、兰考926（104）和剡麦98（31）等；高产高氮敏感型（Ⅱ区），此类品种的单穗籽粒产量高于平均值，D值低于平均值。包括西农158（84）、新麦19（93）、陇辐2号（29）、LS6109（33）、温麦19（89）、郑麦98（60）和冀729（37）等；低产高氮敏感型（Ⅲ区），此类品种的单穗籽粒产量和D值均低于平均值。有石麦 12（74）、丰舞 981（87）、中育 9398（101）、新乡大穗（44）、周麦25（64）、临抗22（66）和丰抗2号（78）等；低产耐高氮型（Ⅳ区），此类品种的单穗产量低于平均值，而D值高于平均值，包括矮杆早（9）、良星99（39）、尧麦16（67）、京冬8号（10）、有芒白4号（23）、农大3492（21）、京411（16）和郑农46（100）。

综合2年结果，根据低、高氮胁迫下的单穗籽粒产量差异，将108份小麦品种分为4种氮效率类型。（1）双高效型，低氮或高氮胁迫下的籽粒产量高于平均值。包括陇麦 212（83）、西农 979（92）、兰考926（104）、西农158（84）、存麦12号（85）、温麦19（89）、剡麦98（31）和陇辐2号（29）等21个品种；（2）低氮高效型，在低氮胁迫下的籽粒产量高于平均值，高氮胁迫下籽粒产量低于平均值，包括石4185（73）、西农585（91）和郑农46（100）等7个品种；（3）高氮高效型，在低氮胁迫下的籽粒产量低于平均值，高氮胁迫下则相反，包括科农2009（38）、长丰1号（13）、邯农2312（36）、郑麦366（99）和中种麦10号（42）等10个品种；（4）双低效型，在低氮和高氮胁迫下的籽粒产量均低于平均值，包括农大 3492（21）、冀 38（71）、金丰 7183（55）、泛麦5号（52）、京冬8号（10）和轮选987（79）等24个品种。

3 讨论

关于氮效率筛选方法的选择，针对不同的作物常采用苗期水培试验[33]、盆栽试验[34]和田间试验[35]，或两两结合的方法控制氮环境条件进行氮高效种质的筛选。相比前两者，田间试验将作物的氮利用效率与产量联系起来，使其更为精准地筛选氮高效种质。选择合适的性状指标是作物氮效率鉴定的关键。冯洋等[35]利用叶片SPAD值、氮素积累量和产量相关性状筛选出不同水稻氮效率种质。崔文芳等[36]研究认为玉米杂交种氮高效筛选的有效指标为叶绿素SPAD值、茎叶总氮量和穗位叶氨基酸浓度。本研究采用大田试验，对108份小麦品种的耐氮性进行全生育期筛选。利用小麦产量性状指标的相对值控制环境因素和遗传背景造成的差异，依据氮胁迫系数和变异系数反映性状指标受氮胁迫的大小和不同品种间受氮胁迫的程度，筛选出单穗产量、旗叶长、穗粒数、旗叶宽、株高、穗长、茎粗和千粒重 8个指标（CV>10%）可以作为小麦成株期耐氮性的评价指标。

根据耐氮性不同将品种划分成不同的类型，有利于品种间的比较。张旭[37]利用氮素吸收效率和氮肥农学效率将 14个小麦品种划分氮高效、中效和低效 3种类型。黄永兰等[38]对氮效率相关性状进行分析，依据不同施氮水平下的产量差异将 45份水稻种质分为双高效型、高氮高效型、低氮高效型和双低效型。崔文芳等[39]以耐低氮胁迫指数和高氮下产量将 27份玉米自交材料划分成4种氮效率类型，高产氮高效型、低产氮高效型、高产氮低效型、低产氮低效型。顾炽明等[40]以苗期生物量的平均值为依据，将162份甘蓝型油菜氮效率分为双高效型、双低效型、中间型。本研究利用相关分析、主成分分析、隶属函数分析和聚类分析等方法，对108份小麦种质成株期的耐氮性进行综合评价，依据耐低氮性综合评价D值将供试小麦品种分为低氮敏感型、低氮较敏感型、中间型、耐低氮型和强耐低氮型5类；以耐高氮性综合评价值将其分为4类，高氮敏感型、高氮较敏感型、中间型和耐高氮型。最终基于籽粒产量和氮综合评价D值将108份供试小麦品种分为4种氮效率类型，双高效型、低氮高效型、高氮高效型和双低效型。对不同氮水平下小麦单穗籽粒产量分析发现，耐低氮型小麦（石4185、丰抗2、西农158和丰抗2）的单穗产量连续2年均增产超过20%，低氮敏感型小麦（郑麦98、金丰7183、中种麦10号、益科麦5号、剡麦98和河麦58-3）的单穗产量连续2年均减产超过20%，结果进一步证实筛选指标的代表性和评价方法的可行性，以及耐氮性类型划分结果的可靠性。该研究结果对小麦氮高效育种具有重要的参考价值，但由于试验无其他生态区的多点实验佐证，因此，具有一定局限性，后续会在黄淮麦区不同点进行鉴定。

4 结论

氮胁迫对小麦的形态和产量构成相关性状指标均有不同程度的影响，其中单穗产量对氮胁迫的敏感程度最大。单穗产量、旗叶长、穗粒数、旗叶宽、株高、穗长、茎粗和千粒重 8个指标（CV>10%）适合作为小麦成株期氮效率评价的初选指标。并将108份小麦品种划分为双高效型、低氮高效型、高氮高效型和双低效型4种氮效率类型，其中，陇麦212、西农979、兰考926、西农158、温麦19等21个小麦品种为双高效型，可作为小麦氮高效育种的种质资源。