重要胡麻栽培品种的抗旱性综合评价及指标筛选
罗俊杰1, 欧巧明1,*, 叶春雷1, 王方1, 王镛臻2, 陈玉梁1,*
1甘肃省农业科学院生物技术研究所, 甘肃兰州 730070
2甘肃农业大学资源与环境学院, 甘肃兰州 730070
* 通讯作者(Corresponding authors): 欧巧明, E-mail:ouqiaoming@163.com, Tel: 0931-7612685; 陈玉梁, E-mail:chenyl925@163.com, Tel: 0931-7616636

第一作者联系方式: E-mail:hnsljjie@163.com, Tel: 0931-7612658

摘要

以15份国内胡麻栽培品种为材料, 设自然降雨和正常灌水2个处理, 考察与抗旱性相关的7个农艺性状、8个生理生化指标及产量指标, 采用综合抗旱系数、因子分析、隶属函数、聚类分析和灰色关联度分析相结合的方法, 对其抗旱性进行综合评价、抗旱型划分和评价指标筛选。结果显示, 相关性状指标对干旱胁迫的反应及关联程度各异, 可优先选择与抗旱性关系密切的产量及其相关性状、光合作用因子、叶片抗氧化因子等相关生理生化性状; 因子分析表明, 6个公因子可代表胡麻抗旱性90.89%的原始数据信息量。基于抗旱性度量值(drought resistance comprehensive evaluation values, D值)和加权抗旱系数(weight drought resistance coefficient, WDC值)的各品种抗旱性排序相近, 位居前6位的抗旱品种相同。各品种D值与综合抗旱系数(comprehensive drought resistance coefficient, CDC值)、WDC值、产量抗旱系数(yield drought resistance coefficient, Y值)之间均呈极显著正相关, 而各品种Y值与CDC、WDC值间极显著正相关; 据D值将供试品种划分为5个抗旱级别, 可较好地反映品种的选育条件及适应地区。试验结果说明以D值为主要参数, 以WDC为辅助评价参数, 评价以产量为主要考量目标的胡麻抗旱性是适宜且准确的; 以抗旱性综合评价方法进行胡麻抗旱性综合评价、抗旱型划分、评价指标筛选是准确的。

关键词: 胡麻; 农艺性状; 生理生化指标; 抗旱性综合评价; 隶属函数
Comprehensive Valuation of Drought Resistance and Screening of Indices of Important Flax Cultivars
LUO Jun-Jie1, OU Qiao-Ming1,*, YE Chun-Lei1, WANG Fang1, WANG Yong-Zhen2, CHEN Yu-Liang1,*
1Biotechnology Institute, Gansu Academy of Agricultural Sciences, Lanzhou 730070, China
2College of Resources and Environmental Sciences, Gansu Agricultural University, Lanzhou 730070, China
Abstract

Seven main agronomic traits, eight physiological indices and yield index of 15 main flax cultivars in China were measured during maturity under the conditions of irrigation and natural rainfall. Comprehensive drought resistance coefficient (CDC value), factor analysis, subordinate function coefficients, clustering analysis, grey relational analysis were used to evaluate the drought resistance, classify drought resistance type and select evaluation indices in tested flax cultivars. The response to drought stress and correlations of tested traits and indices were different. Yield and photosynthesis factors as well as leaf antioxidant factors were closely associated with drought resistance, so could be used as priority indicator of drought resistance evaluation. Factor analysis showed that six common factors could represent 90.89% of the original information of flax drought resistance data. The ranks of drought resistance of tested flax cultivars based on drought resistance comprehensive evaluation values (D value) and weight drought resistance coefficient (WDC value) were similar, the drought resistance was the same for six flax cultivars, arranging from the first to third in drought resistance. D values of tested flax cultivars had significant and positive correlation with CDC value, WDC value and Y value. Y values of tested flax cultivars also had very significant and positive correlation with CDC value and WDC value. According to D value clustering analysis, tested cultivars were divided into five grades in drought resistance, reflecting the diffirence of cultivars in breeding condition and adaptive region growing well. Drought resistance evaluated mainly with yield by D value as index and WDC value as auxiliary index in flax were appropriate and accurate. Drought resistance comprehensive evaluation methods can be used in studying exactly drought resistance evaluation, classification of drought resistant type and screening evaluation indices in flax.

Keyword: Flax; Agronomic traits; Physiological and biochemical indices; Drought resistance comprehensive evaluation; Subordinate function coefficients

胡麻( Linum usitatissimumL.)即油用亚麻, 是中国西北、华北等年降水量300~500 mm的旱作农业区重要的油料作物[ 1, 2, 3], 具有耐瘠薄、喜阴凉、耐旱、耐寒等特点, 在抗旱资源方面具有独特研究价值。研究国内胡麻品种抗旱性及其综合评价方法和鉴定指标, 对胡麻抗旱育种、抗旱资源利用及品种的生产应用与合理布局具有重要意义。

作物抗旱性鉴定需要将形态、生理生化、产量等指标相结合, 且对各个时期的抗旱性综合评价[ 4]。因此, 简单、有效的鉴定指标及其评价方法的合理选择是抗旱性鉴定的关键[ 2]。近年来, 国内外学者从不同角度研究和提出了作物抗旱性鉴定的试验设计、指标选择、评价方法及抗旱性分级等问题[ 1, 2, 5, 6, 7, 8, 9, 10], 针对不同作物和不同育种需求的更加全面的抗旱评价方法已在小麦[ 1, 11]、大豆[ 12, 13]、油菜[ 14, 15]、谷子[ 16]、水稻[ 9]、玉米[ 17]、棉花[ 18]、胡麻[ 2, 19]等作物上应用, 如主成分分析法、隶属函数法、聚类分析法等。目前普遍认为, 多指标多方法相结合的抗旱性综合评价更加真实、可靠[ 2]

前人研究多侧重于以胡麻种质资源为材料[ 2]或针对某种农艺性状、生理生化指标单独分析, 而针对胡麻抗旱相关生理指标且采用抗旱性综合评价方法的研究鲜见报道。为了探讨国内胡麻栽培品种抗旱性及其与主要性状间的关系。本研究在分析干旱处理和灌水条件下与胡麻抗旱性相关的多个农艺性状及生理生化性状基础上, 对国内15个胡麻栽培品种的抗旱性综合评价, 试图明确其抗旱性评价中可借鉴的性状指标及简单、可靠的方法, 以期为抗旱育种及品种生产应用提供参考。

1 材料与方法
1.1 试验材料与处理

分别于2011年和2012年在甘肃省兰州市榆中县良种场种植15份目前国内主要胡麻种植区大面积应用的育成品种(表1)。试验点海拔1480 m、年平均气温6.7℃、无霜期148 d、年平均降水量 381.6 mm、年蒸发量1343.1 mm。所选试验地含有机质17.13 g kg-1、全氮0.94 g kg-1、全磷0.76 g kg-1、全钾22.98 g kg-1、碱解氮38.65 mg kg-1、速效磷24.37 mg kg-1、速效钾132.67 mg kg-1, pH 7.88。

表1 干旱胁迫处理和正常灌溉处理条件下的供试胡麻品种主要农艺性状性状测定值及其均值差异分析 Table 1 Measured value of main agronomic traits of tested flax cultivars under drought stress and irrigation and its mean variance analysis

采用随机区组排列, 分别设干旱胁迫(自然降雨, 全生育期不灌溉)和对照(全生育期灌水2次) 2个处理, 重复3次。每处理15个品种, 小区面积24 m2(6 m×4 m), 行距20 cm, 于4月7日人工开沟撒播, 种植密度900万粒 hm-2, 干旱胁迫与对照处理间隔5 m, 种植保护区, 处理间小区边界均用培土分离, 播种前施尿素300 kg hm-2和磷酸二氢铵150 kg hm-2, 按当地大田生产管理。

1.2 试验区胡麻生育期降水及土壤水分状况

图1反映了2011—2012年试验区胡麻全生育期(1月至8月)降雨量的动态趋势, 反映出本试验区降雨时空分布不均衡的特征, 说明自然降雨是胡麻生长季节干旱胁迫时水分的唯一补给来源。图2反映了2011—2012年试验区田间土壤含水量(soil water content, SWC)测试结果及其动态趋势, 相同处理不同土层深度SWC间差异明显(图2); 同一土层深度在盛花前期(6月9日) (均值差 dj = -1.344, t = 0.474, P = 0.478)和盛花后期(7月16日) (均值差 dj = -0.278, t = 0.4203, P = 0.6853) 在干旱胁迫和正常灌水处理间SWC差异不显著, 但在盛花期(6月27日) (均值差 dj = -6.578, t = 17.364, P = 0.0001)差异达极显著水平(图2), 相同处理3个测定时间的SWC差异亦达显著或极显著水平。说明干旱胁迫处理在胡麻盛花期SWC下降最明显, 水分透支严重, 胡麻受干旱胁迫亦最显著, 说明在该区域正常年份利用自然降雨控制田间水分, 实现对胡麻的干旱胁迫并加以研究是可行的。所以选择以盛花期筛选抗旱指标及评价抗旱性。

图1 2011-2012年试验区胡麻全生育期(1月至8月)旬降雨量动态Fig. 1 Dynamic of average rainfall of each ten-day at the whole growth period (from Jan. to Aug.) of flax in 2011-2012

图2 胡麻生育期不同土层深的土壤水分变化Fig. 2 Soil moisture changes of different depth in different seasons of flax growth period

1.3 考察性状与方法

在胡麻盛花期对每小区每品种材料取样, 参考邹琦[ 20]的《植物生理学实验指导》的方法, 分别测定相关生理生化性状(表3)。采用80%丙酮-无水乙醇浸提比色法(mg g-1)测定叶片叶绿素a (chlorophyll a, Chl a)、叶绿素b (chlorophyll b, Chl b)和总叶绿素(total chlorophyll, TChl)含量, 采用硫代巴比妥酸比色法(μmol mg-1)测定丙二醛含量(malondialdehyde, MDA); 采用紫外吸收法(U mg-1 min-1)、抑制NBT光还原比色法(U mg-1)、愈创木酚法(U mg-1 min-1)分别测定过氧化氢酶(catalase, CAT)活性、超氧化物歧化酶(superoxide dismutase, SOD)活性、过氧化物酶(peroxidase, POD)活性, 采用电导仪法测定细胞膜透性。每小区每品种材料取样3份, 每份样品重复测定3次, 并计算平均值, 最后取两年度平均值作为供试指标测定值。

在胡麻成熟后每小区每品种取样20株考种, 参照《胡麻种质资源描述规范和数据标准》[ 21], 测定株高、工艺长度、有效分枝数、单株果数、不实果数、单株产量、千粒重等产量相关农艺性状, 全区收获统计小区产量, 求其平均值, 最后取两年度平均值作为供试性状测定值。

1.4 数据统计分析

参考兰巨生[ 22]、祁旭升等[ 2]、罗俊杰等[ 18]、谢小玉等[ 15]、孟庆立等[ 16]和尹利等[ 23]的方法, 以2011和2012年供试胡麻品种各小区农艺性状和生理生化性状测定值作为基础数据, 对各性状测定值进行平均数差异显著性分析, 按公式(1)和(2)分别计算单项抗旱系数(drought resistance coefficient, DC)及综合抗旱系数(comprehensive drought resistance coefficient, CDC), 式中 x i CK i表示干旱、对照处理的性状测定值。

DC=

, i=1, 2, 3, …, n (1)

CDC=

, i=1, 2, 3, …, n (2)

按公式(3)计算因子权重系数( ω i), 式中 P i为第 i个综合指标贡献率, 表示第 i个指标在所有指标中的重要程度。

ω i=

, i=1, 2, 3, …, n(3) (3)

按公式(4)计算各基因型各综合指标的隶属函数值[ μ( x i)], 式中 x i x imin x imax分别表示第 i个综合指标及第 i个综合指标的最小值、最大值。

μ( x i)=

, i=1, 2, 3, …, n(4)

根据因子权重( ω i)及隶属函数值[ μ( x i)], 按公式(5)计算抗旱性度量值(drought resistance compre-hensive evaluation values, D)。

D=

, i=1, 2, 3, …, n (5)

按公式(6)、(7)计算关联系数(ξ i)及关联度(γ i), 式中 为二级最小差, 为二级最大差, ρ为分辨系数(取值0.5)。

ξ i( k)=

, k=1, 2, 3, …, n; i=1, 2, 3, …, n (6)

γ i=

, k=1, 2, 3, …, n (7)

根据各性状关联度 γ i, 按公式(8)计算各性状指标权重系数[ ω i( γ)], 并据公式(9)计算加权抗旱系数(weight drought resistance coefficient, WDC), 式中 γ i为各基因型关联度。

ω i( γ)=

, i=1, 2, 3, …, n (8)

WDC=

, i=1, 2, 3, …, n (9)

最后针对供试胡麻品种D值, 采用欧式距离和加权配对算术平均法(weighted pair group method average, WPGMA)进行聚类分析, 划分抗旱级别; 并分别以D值、CDC值、WDC值为参考序列, 对各性状DC值进行逐步回归分析, 求取回归方程。

采用Microsoft Excel 2003和SPSS 18.0分析软件处理数据。

2 结果与分析
2.1 试验材料的代表性及其性状测定值分析

供试胡麻品种主要农艺性状两年度测定平均值的差异显著性分析(配对处理 t检验)表明(表1) , 干旱对供试品种主要农艺性状测定值均有显著影响, 除有效分枝数外( P>0.05), 处理间差异均达极显著水平( P<0.01), 品种间差异亦达显著水平( P<0.05), 品种间变异系数介于0.061~0.543。说明本试验所选胡麻品种类型较丰富, 所选主要农艺性状对干旱胁迫反应较为敏感, 干旱处理试验效果较好, 具有较好代表性。

供试品种主要生理指标2年度测定平均值的差异显著性分析表明(表2), 干旱对供试品种生理指标测定值均有明显影响, 除TChl、CAT外( P>0.05), 处理间差异均达显著( P<0.05)或极显著水平外( P<0.01), 品种间差异亦达显著水平( P<0.05), 品种间变异系数介于0.037~1.382。说明本试验干旱处理效果较好, 但所选主要生理指标对干旱胁迫反应敏感性不同, 除TChl、CAT外, 其余指标均较为敏感, 具有较好代表性。

表2 干旱胁迫处理和正常灌溉处理条件下的供试胡麻品种主要生理指标测定值及其均值差异分析 Table 2 Measured values of main physiological indicatros of tested falx cultivars under drought stress and irrigation and its mean variance analysis

另外, 表1表2还显示供试品种各性状在干旱和正常处理下的测定值相关性差异较大, 这进一步说明各性状敏感性存在较大差异, 采用各性状测定值难以直接考量其抗旱性。

2.2 单项指标与胡麻品种的抗旱性分析

根据1.4节中公式(1)获得各性状抗旱系数(DC值)(表3), 分析表明, 与正常灌水处理相比, 供试品种在生育期持续干旱处理后, 所选性状均发生了不同程度变化。同一性状各品种的DC值差异明显, 变异系数介于0.020~0.385之间, 但不同品种间DC值所反映的抗旱性不同, 且同一品种各性状的DC值存在较大差异(表3), 说明各性状指标对干旱胁迫的敏感程度各异, 且所测指标间的关系复杂, 有必要借助综合抗旱系数(CDC值)及因子分析确定其权重系数( ω)后进行抗旱性综合评价。

表3 供试胡麻品种、生育期及各农艺性状和生理生化指标的抗旱系数 Table 2 Drought resistance coefficients of agronomic traits and physiological indicatros in flax cultivars

采用综合抗旱系数公式(2)获得供试品种CDC值, CDC值介于0.770~1.024之间, 变异系数为0.079。根据CDC值大小对供试品种进行抗旱性排序(表3表6), 15个供试品种的抗旱性强弱顺序为天亚9号、晋亚7号、坝亚12、伊亚4号、轮选2号、宁亚19、陇亚杂2号、陇亚10号、陇亚8号、陇亚11、定亚23、晋亚10号、晋亚11、坝选3号。

通过各性状DC值的连续变数次数分布统计分析[ 2], 表明同一区间各性状DC值分布次数和频率相差较大, DC值在0.9表4)。说明千粒重、有效分枝数对干旱胁迫的反应迟钝, 株高、工艺长度、MDA、Chl b次之, 而其他性状则较为敏感。

表4 供试胡麻品种各性状指标的抗旱系数及其在不同区间的分布 Table 4 Drought resistance coefficient of traits indices and diffirent distributions of tested flax cultivars
2.3 因子分析及因子载荷与权重系数评价

采用因子分析对各性状指标DC值进一步分析, 获得特征向量、因子载荷和贡献率(表5), 以便进一步确定和评价抗旱鉴定的重要指标。结果表明, 各因子特征值中前6个因子的累积贡献率已达到90.89%, 其特征根 λ>0.6371。抽取前6个因子, 将具有相同本质的变量归为一类, 可将原来各单项性状指标转换为6个新的相互独立的综合指标(公因子, 用F1、F2、F3、F4、F5、F6表示)。公因子载荷矩阵的载荷系数绝对值表明(表2), F1在单株果实、单株产量、千粒重、Chl a、CAT上有较高荷载量; F2在株高、工艺长度、TChl上有较高载荷量; F3在MDA、Chl b、SOD、产量上有较高载荷量; F4在膜透性、POD上有较高载荷量; F5、F6分别在不实果数、有效分枝数上有较高载荷量。

表5 因子载荷矩阵及方差贡献率 Table 5 Rotated factor pattern and cumulative variance contribution

在获得各公因子特征向量、贡献率的基础上, 据公式(3)计算各公因子权重系数(ω)(表5), 用以计算加权隶属函数值。

2.4 隶属函数分析及其品种的综合抗旱性评价

在获得各公因子特征向量的基础上, 利用模糊隶属函数法, 据公式(4)计算各因子隶属函数值( μ), 并据各因子权重系数( ω), 对各因子隶属函数值( μ)赋予相应权重系数( ω), 按公式(5)计算加权隶属函数值, 作为综合抗旱性度量值(D值), 据此可较准确地评价供试品种的抗旱性, D值越大抗旱性越强(表6)。

表6 不同胡麻品种抗旱性评价的D值、CDC值、WDC值及其预测值 Table 6 Subordinate function value,D value,CDC value,WDC value,and their predict value on flax cultivars based of drought resistance evaluation

供试品种D值介于0.366~0.738之间, 平均值0.488, 变异系数为0.191。根据D值大小对供试胡麻品种进行抗旱性排序(表6), 其顺序为天亚9号、晋亚7号、轮选2号、伊亚4号、坝亚12、宁亚19、陇亚10号、陇亚杂2号、晋亚11、陇亚8号、晋亚10号、定亚23、陇亚11、陇亚杂1号、坝选3号。这与基于CDC值的供试品种的抗旱性强弱评价结果基本吻合。而抗旱性最强的天亚9号D值较所有品种平均值高出11.998%, 远高于其他品种, 这与各品种在生产应用中的多年抗旱性表现相吻合。

2.5 不同性状与D值的灰色关联分析

按照公式(6)和(7), 对所有性状单项抗旱系数(比较数列)与其D值(参考数列)进行灰色关联度分析, 获得各性状DC值与D值间的关联度( γD), 结果显示各性状与D值的密切程度依次为单株果实、千粒重、CAT、Chl a、产量、单株产量、MDA、有效分枝、不实果数、SOD、Chl b、膜透性、TChl、株高、工艺长度、POD (表7), 这和各性状DC值与Y值的关联程度相近。

表7 供试胡麻品种各性状抗旱系数与D值、WDC值的关联度及各性状指标权数 Table 7 Correlation degree between drought resistance coefficients of tested traits and D value, WDC value and traits
weight of flax cultivars

上述结果还表明, 与胡麻品种抗旱性关系最为密切的是与产量相关的性状(如单株果实、千粒重、单株产量、产量), 其次是与叶片抗氧化和光合作用有关的生理生化性状(CAT、MDA、Chl a), 它们对干旱胁迫反映最为直接, 且所受影响最大, 因而与抗旱性关联度最大; 而株高、工艺长度与抗旱性关联较弱。

对各性状关联度( γD)归一化处理后获得性状指标权重系数[ ω i( γ)], 进而获得供试品种加权抗旱系数(WDC值)(表4), 并基于WDC值获得各性状关联度( γWDC)(表7), 用以辅助评价各性状与品种抗旱性的关联度及各品种抗旱性。各性状关联度γWDC与各性状 γD达极显著正相关( r=0.93), 说明各性状同WDC值的密切程度与各性状同D值的密切程度相近, 用WDC值辅助进行品种抗旱性评价是适宜且准确的。

另外, 除株高、工艺长度、不实果数、膜透性、TChl、POD外, 其余性状DC值与CDC值、WDC值之间均呈极显著正相关性( P<0.01), 这和所有性状DC值与D值、Y值均呈极显著正相关的结果相当吻合。说明上述6个性状指标对干旱胁迫的敏感程度相对较小, 但这并不足以认定它们完全不宜作为供试品种抗旱性评价的指标。

2.6 聚类分析及抗旱级别的划分

图3表明参试15份胡麻品种按抗旱性在0.60处分为5类, 其中I级抗旱型1份, II级抗旱型4份, Ⅲ级抗旱型3份, IV级抗旱型6份, V级抗旱型1份。

图3 基于D值的供试胡麻品种抗旱性系统聚类图(WPGMA法)Fig. 3 Fuzzy clustering dendrogram (WPGMA) of drought resistance of tested flax cultivars based on D values

按照供试胡麻品种的抗旱性聚类结果及抗旱级别, 对供试品种及其隶属函数( μ)、抗旱性综合评价指标CDC值、D值及WDC值分级显示(表8), 大部分性状指标的隶属函数值( μ)及上述3个评价指标均随抗旱级别的升高而降低, 但个别性状亦存在异常变化。上述3个评价指标在不同抗旱级别上的差异较为明显, 可为其他胡麻品种抗旱级别划分及不同研究者的相关研究结果之间提供可比较的依据。

表8 供试胡麻品种不同抗旱性综合评价指标的分级及其评价数据平均值 Table 8 Classification of drought resistance comprehensive evaluation indices of tested flax cultivars and their
average of evaluation data
2.7 逐步回归分析及抗旱性预测评价

分别利用供试胡麻品种的D值、CDC值、WDC值与所有性状指标抗旱系数进行逐步回归分析, 得到回归方程(表9), 进而据回归方程(1)、(2)、(3)对各指标隶属函数的预测D值与D值、预测CDC值与CDC值、预测WDC值与WDC值做相关性分析, 两两间均达极显著正相关( P < 0.01); 统计分析显示模型方程(1)、(2)、(3)的决定系数 R2≈ 1, F检验均达极显著水平( P < 0.01), 说明预测值与实际值之间拟合度好(Durbin-Watson统计量 d ≈ 2, P < 0.01), 回归方程最优, 其解释能力强, 预测精度高, 用这3个方程进行胡麻品种抗旱性评价预测效果好。

表9 不同胡麻品种抗旱性预测模型及不同综合评价指标间的相关性 Table 9 Predict model of drought resistance and correlation of diffirent comprehensive valuation indices of flax cultivars

据回归方程(1)可知, 在胡麻品种抗旱性鉴定中, 有选择地测定与D值密切相关的指标, 如单株果实、不实果数、MP、TChl、CAT等5个指标, 可有效鉴定胡麻品种的抗旱性, 从而使鉴定工作简单化。

2.8 不同抗旱性综合评价方法及度量指标的比较

基于不同综合评价指标的抗旱性评价结果, 基于WDC值和基于D值, 供试胡麻品种抗旱性排序相近, 与位居前6位的抗旱品种相同。相关性分析显示(表9), 供试品种D值、WDC值、Y值两两之间均呈极显著正相关( P<0.01), 说明以D值为主要抗旱性综合评价指标, 以WDC值作为辅助综合评价指标, 进行以产量为主要考量目标的胡麻抗旱性评价是准确的。

另外, 相关性分析还显示, CDC值与D值、Y值亦达极显著正相关( P<0.01), 但考虑到其在计算过程中未考量各性状指标的敏感程度及权重, 仅据此将其作为抗旱性综合评价指标有待进一步研究。

3 讨论
3.1 关于抗旱评价指标的选择问题

性状指标的合理选择是作物抗旱性鉴定的关键, 对此目前国内外进行了大量研究[ 1, 2, 7, 8, 9, 10, 11, 12, 13, 14, 15, 16, 17, 18, 19], 并从不同角度提出了抗旱性鉴定的相关性状。Zhang等[ 24]和孟庆立等[ 16]对谷子成株期抗旱性研究认为, 产量性状对干旱胁迫最敏感。Chernyad’ev等[ 25]和王士强等[ 1]认为小麦生理生化性状的影响大于农艺性状; Tambussi等[ 26]和李贵全等[ 12]研究表明大豆农艺性状的综合指标与品种的抗旱性显著相关。然而, 不同鉴定指标对干旱胁迫的敏感程度、不同指标间有何关系, 对抗旱性的贡献等问题, 目前尚无定论, 且多数研究未将农艺性状和生理指标结合起来考虑。

关联度反映的是构成该系统的各性状组成的比较数列和参考数列间的密切程度[ 1, 23]。本研究选取与胡麻抗旱性相关的16个关键农艺性状及生理生化性状为指标, 其结合起来, 从宏观和微观层面对抗旱性的内在机制认识和分析, 结果显示, 与抗旱性度量值(D值)的关联程度较高的6个性状依次为单株果实、千粒重、CAT、Chl a、产量抗旱系数、单株产量, 这和各性状抗旱系数与加权抗旱系数(WDC值)的密切程度相吻合, 说明基于D值和WDC值的性状关联度评价是适宜且准确的, 进而提示我们, 除以产量等核心农艺性状作为胡麻抗旱性筛选指标外, 还应重视CAT、MDA、Chl a等生化指标的分析, 以提高判断的准确性, 这与前人的研究结果是一致的[ 9, 15]。而与小麦等作物明显不同的是, 胡麻株高、工艺长度等性状与其抗旱性关联较弱, 不宜作为品种抗旱性选择的核心指标。

从抗旱性状因子分析的结果可以看出, 单株果实、单株产量、千粒重、Chl a、CAT占据了首要位置, 表明在干旱条件下产量相关性状是胡麻不同基因型间抗旱性的重要体现, 这与祁旭升等[ 2]的研究结果相近; 而干旱条件下的地上部分光合作用和抗氧化能力与其抗旱性密切相关。可见, 除丰产性是胡麻抗旱性重要外在表现外, 良好的逆境生理调控能力、繁茂的生物积累量和流畅的根系物质转储能力[ 27]是影响胡麻抗旱性的重要因素, 这与孟庆立等[ 16]的研究结论相近。

3.2 作物抗旱性的分析方法

近年国内外学者针对不同作物和不同育种需求从不同角度研究和提出了更加全面的抗旱评价鉴定方法[ 2, 18, 28, 29, 30, 31]。目前普遍认为, 多指标多方法相结合的抗旱性综合评价比较可靠[ 2]。本研究认为各性状指标权重及最终核心综合评价指标的确定往往是抗旱性综合评价的关键。

本研究采用D值、综合抗旱系数(CDC值)、WDC值、Y值等综合评价其抗旱性, 确定D值及WDC值为综合评价指标, 进行较确切的聚类分析和抗旱型划分, 并结合灰色关联分析、逐步回归分析筛选关键性状指标, 建立了拟合度较好的回归方程。结果显示, 以D值为主要综合评价指标、以WDC作为辅助综合评价指标, 评价以产量为主要考量目标的胡麻抗旱性是适宜且准确的。虽然CDC值与D值亦达极显著正相关, 但在其计算过程中未考量各性状指标的敏感程度、重要性及权重, 是否能将其作为抗旱性综合评价指标有待进一步研究确认。另外, 本研究依据D值进行抗旱级别的划分, 并提供相应抗旱级别的D值、CDC值和WDC值, 可为其他胡麻品种抗旱级别划分及不同研究者的相关研究结果之间提供可比较的依据。

本研究结果也说明产量相关性状仍是决定胡麻抗旱性的根本因素, 这与孟庆立等[ 16]对谷子成株期抗旱性研究结论相同。而利用本研究所建立最优回归方程可有效预测所测品种抗旱性的强弱, 提高抗旱育种效率。故笔者认为, 基于D值的抗旱性鉴定和抗旱型划分能很好地反映各品种抗旱性、抗旱特点及选育和应用区域, 进而提示我们采用隶属函数法及综合/加权抗旱系数法, 并结合灰色关联分析法评价胡麻品种抗旱性, 结果客观、可靠。

4 结论

采用抗旱性综合评价方法, 确定抗旱性度量值(D值)及加权抗旱系数(WDC值)为综合评价指标, 进行聚类分析和抗旱型划分, 并结合灰色关联分析、逐步回归分析筛选出关键性状指标, 建立了回归方程, 能有效反映各参试品种抗旱性、抗旱特点, 为参试品种在育种及适应区域的合理利用及其他品种(系)的抗旱性鉴定提供有价值的参考。

The authors have declared that no competing interests exist.

作者已声明无竞争性利益关系。

参考文献
[1] 王士强, 胡银岗, 佘奎军, 周琳璘, 孟凡磊. 小麦抗旱相关农艺性状和生理生化性状的灰色关联度分析. 中国农业科学, 2007, 40: 2452-2459
Wang S Q, Hu Y G, She K J, Zhou L L, Meng F L. Gray relational grade analysis of agronomical and physi-biochemical traits related to drought tolerance in wheat. Sci Agric Sin, 2007, 40: 2452-2459 (in Chinese with English abstract) [本文引用:6] [CJCR: 1.889]
[2] 祁旭升, 王兴荣, 许军, 张建平, 米君. 胡麻种质资源成株期抗旱性评价. 中国农业科学, 2010, 43: 3076-3087
Qi X S, Wang X R, Xu J, Zhang J P, Mi J. Drought-resistance evaluation of flax germplasm at adult plant stage. Sci Agric Sin, 2010, 43: 3076-3087 (in Chinese with English abstract) [本文引用:12] [CJCR: 1.889]
[3] 赵利, 党占海, 张建平, 关天霞, 田彩萍. 不同类型胡麻品种资源品质特性及其相关性研究. 干旱地区农业研究, 2008, 26(5): 6-9
Zhao L, Dang Z H, Zhang J P, Guan T X, Tian C P. Study on quality characters and correlation of different types of flax germplasm. Agric Res Arid Areas, 2008, 26(5): 6-9 (in Chinese with English abstract) [本文引用:1]
[4] 张木清, 陈如凯. 作物抗旱分子生理与遗传改良. 北京: 科学出版社, 2005. pp22-23
Zhang M Q, Chen R K. Molecular Physiological and Genetic Improvement of Crop Drought Tolerance. Beijing: Science Press, 2005. pp22-23(in Chinese) [本文引用:1]
[5] 张正斌. 作物抗旱节水的生理遗传育种基础. 北京: 科学出版社, 2003. pp67-209
Zhang Z B. Fundamentals of Physiological and Genetic Breeding Basis in Crop Drought Resistance and Water Saving. Beijing: Science Press, 2003. pp67-209(in Chinese) [本文引用:1]
[6] Kamoshita A, Babu R C, Boopathi N M, Fukai S. Phenotypic and genotypic analysis of drought-resistance traits for development of rice cultivars adapted to rainfed environments. Field Crops Res, 2008, 109: 1-23 [本文引用:1] [JCR: 2.474]
[7] Upadhyaya H D. Variability for drought resistance related traits in the mini core collection of peanut. Crop Sci, 2005, 45: 1432-1440 [本文引用:2] [JCR: 1.513]
[8] Hura T, Hura K, Grzesiak S. Physiological and biochemical parameters for identification of QTLs controlling the winter triticale drought tolerance at the seedling stage. Plant Physiol Biochem, 2009, 47: 210-214 [本文引用:2]
[9] 王贺正, 马均, 李旭毅, 李艳, 张荣萍, 汪仁全. 水稻开花期一些生理生化特性与品种抗旱性的关系. 中国农业科学, 2007, 40: 399-404
Wang H Z, Ma J, Li X Y, Li Y, Zhang R P, Wang R Q. Relationship between some physiological and biochemical characteristics and drought tolerance at rice flowering stage. Sci Agric Sin, 2007, 40: 399-404 (in Chinese with English abstract) [本文引用:4] [CJCR: 1.889]
[10] 武斌, 李新海, 肖木辑, 谢传晓, 郝转芳, 李明顺, 张世煌. 53份玉米自交系的苗期耐旱性分析. 中国农业科学, 2007, 40: 665-676
Wu B, Li X H, Xiao M J, Xie C X, Hao Z F, Li M S, Zhang S H. Genetic variation in fifty three maize inbred lines in relation to drought tolerance at seedling stage. Sci Agric Sin, 2007, 40: 665-676 (in Chinese with English abstract) [本文引用:2] [CJCR: 1.889]
[11] 白志英, 李存东, 孙红春, 赵金锋. 小麦代换系抗旱生理指标的主成分分析及综合评价. 中国农业科学, 2008, 41: 4264-4272
Bai Z Y, Li C D, Sun H C, Zhao J F. Principal component analysis and comprehensive evaluation on physiological indices of drought resistance in wheat substitution. Sci Agric Sin, 2008, 41: 4264-4272 (in Chinese with English abstract) [本文引用:2] [CJCR: 1.889]
[12] 李贵全, 张海燕, 季兰, 赵二开, 刘建兵, 李玲, 张家蓉. 不同大豆品种抗旱性综合评价. 应用生态学报, 2006, 17: 2408-2412
Li G Q, Zhang H Y, Ji L, Zhao E K, Liu J B, Li L, Zhang J R. Comprehensive evaluation on drought-resistance of different soybean varieties. Chin J Appl Ecol, 2006, 17: 2408-2412 (in Chinese with English abstract) [本文引用:3] [CJCR: 1.742]
[13] 孔照胜, 武云帅, 岳爱琴, 李贵全, 彭永康. 不同大豆品种抗旱性生理指标综合分析. 华北农学报, 2001, 16(3): 40-45
Kong Z S, Wu Y S, Yue A Q, Li G Q, Peng Y K. Comprehensive analysis of physiological indexes for drought resistance in different soybean varieties. Acta Agric Boreali-Sin, 2001, 16(3): 40-45 (in Chinese with English abstract) [本文引用:2] [CJCR: 0.951]
[14] 朱宗河, 郑文寅, 张学昆. 甘蓝型油菜耐旱相关性状的主成分分析及综合评价. 中国农业科学, 2011, 44: 1775-1787
Zhu Z H, Zheng W Y, Zhang X K. Principal component analysis and comprehensive evaluation on morphological and agronomic traits of drought tolerance in rapeseed (Brassica napus L. ). Sci Agric Sin, 2011, 44: 1775-1787 (in Chinese with English abstract) [本文引用:2] [CJCR: 1.889]
[15] 谢小玉, 张霞, 张兵. 油菜苗期抗旱性评价及抗旱相关指标变化分析. 中国农业科学, 2013, 46: 476-485
Xie X Y, Zhang X, Zhang B. Evaluation of drought resistance and analysis of variation of relevant parameters at seedling stage of rapeseed (Brassica napus L. ). Sci Agric Sin, 2013, 46: 476-485 (in Chinese with English abstract) [本文引用:4] [CJCR: 1.889]
[16] 孟庆立, 关周博, 冯佰利, 柴岩, 胡银岗. 谷子抗旱相关性状的主成分与模糊聚类分析. 中国农业科学, 2009, 42: 2667-2675
Meng Q L, Guan Z B, Feng B L, Chai Y, Hu Y G. Principal component analysis and fuzzy clustering on drought tolerance related traits of foxtail millet (Setaria italica). Sci Agric Sin, 2009, 42: 2667-2675 (in Chinese with English abstract) [本文引用:6] [CJCR: 1.889]
[17] 徐蕊, 王启柏, 张春庆, 吴承来. 玉米自交系抗旱性评价指标体系的建立. 中国农业科学, 2009, 42: 72-84
Xu R, Wang Q B, Zhang C Q, Wu C L. Drought resistance evaluation system of maize inbred. Sci Agric Sin, 2009, 42: 72-84 (in Chinese with English abstract) [本文引用:2] [CJCR: 1.889]
[18] 罗俊杰, 石有太, 陈玉梁, 王红梅, 刘新星. 甘肃不同色彩陆地棉抗旱指标筛选及评价研究. 核农学报, 2012, 26: 952-959
Luo J J, Shi Y T, Chen Y L, Wang H M, Liu X X. Screening and wvaluation of drought tolerant indices of colored upland cotton (Gossypium hirsutum L. ) in Gansu. Acta Agric Nucl Sin, 2012, 26: 952-959 (in Chinese with English abstract) [本文引用:4]
[19] 石仓吉. 亚麻品种抗旱性评价研究. 干旱地区农业研究, 2008, 26(5): 1-5
Shi C J. Research on assessment of flax drought resistance. Agric Res Arid Areas, 2008, 26(5): 1-5 (in Chinese with English abstract) [本文引用:2]
[20] 邹琦. 植物生理学实验指导. 北京: 中国农业出版社, 2000. pp41-184
Zou Q. Guidebook of Plant Physiology Experiments. Beijing: China Agriculture Press, 2000. pp41-184(in Chinese) [本文引用:1]
[21] 王玉富, 粟建光. 亚麻种质资源描述规范和数据标准. 北京: 中国农业出版社, 2006. pp1-99
Wang Y F, Su J G. Descriptors and Data Stand ard for Flax (Linum usitatissimum L. ). Beijing: China Agriculture Press, 2006. pp1-99(in Chinese) [本文引用:1]
[22] 兰巨生. 农作物综合抗旱性评价方法的研究. 西北农业学报, 1998, 7(3): 85-87
Lan J S. Comparison of evaluating methods for agronomic drought resistance in crops. Acta Agric Boreali-OccidentSin, 1998, 7(3): 85-87 (in Chinese with English abstract) [本文引用:1] [JCR: 1.44]
[23] 尹利, 逯晓萍, 傅晓峰, 李美娜, 郭建. 高丹草杂交种灰色关联分析与评判. 中国草地学报, 2006, 28: 21-25
Yin L, Lu X P, Fu X F, Li M N, Guo J. The grey relation analysis and evaluation of hybrid pacesetter. Chin J Grassl, 2006, 28: 21-25 (in Chinese with English abstract) [本文引用:2]
[24] 张霞, 谢小玉. PEG胁迫下甘蓝型油菜种子萌发期抗旱鉴定指标的研究. 西北农业学报, 2012, 21(2): 72-77
Zhang X, Xie X Y. Studies on identification indexes of drought resistance by PEG during seed germination of rapeseed (Brassica napus L. ). Acta Agric Boreali-Occident Sin, 2012, 21(2): 72-77 (in Chinese with English abstract) [本文引用:1]
[25] Chernyad’ev I I, Monakhova O F. Effects of cytokinin preparations on the pools of pigments and proteins of wheat cultivars differing in their tolerance to water stress. Appl Biochem Micro, 2003, 39: 524-531 [本文引用:1]
[26] Tambussi E A, Nogués S, Araus J L. Ear of durum wheat under water stress: water relations and photosynthetic metabolism. Planta, 2005, 221: 446-458 [本文引用:1] [JCR: 3.347]
[27] Manschadi A M, Hammer G L, Christopher J T, de Voil P. Genotypic variation in seedling root architectural traits and implications for drought adaptation in wheat (Triticum aestivum L. ). Plant Soil, 2008, 303: 115-129 [本文引用:1] [JCR: 2.638]
[28] Condon A G, Richards R A, Rebetzke G J, Farquhar G D. Breeding for high water use efficiency. J Exp Bot, 2004, 55: 2447-2460 [本文引用:1] [JCR: 5.242]
[29] Subrahmanyam D, Subash N, Haris A, Sikka A K. Influence of water stress on leaf photosynthetic characteristics in wheat cultivars differing in their susceptibility to drought. Photosyn thetica, 2006, 44: 125-129 [本文引用:1]
[30] 陈玉梁, 石有太, 罗俊杰, 王蒂, 厚毅清, 李忠旺, 张秉贤. 甘肃彩色棉花抗旱性农艺性状指标的筛选鉴定. 作物学报, 2012, 38: 1680-1687
Chen Y L, Shi Y T, Luo J J, Wang D, Hou Y Q, Li Z W, Zhang B X. Screening of drought tolerant agronomic trait indices of colored cotton varieties (lines) in Gansu Province. Acta Agron Sin, 2012, 38: 1680-1687 (in Chinese with English abstract) [本文引用:1] [CJCR: 1.667]
[31] 陈玉梁, 石有太, 罗俊杰, 李忠旺, 厚毅清, 王蒂. 干旱胁迫对彩色棉花农艺、品质性状及水分利用效率的影响. 作物学报, 2013, 39: 2074-2082
Chen Y L, Shi Y T, Luo J J, Li Z W, Hou Y Q, Wang D. Effect of drought stress on agronomic traits, quality, and WUE in different colored upland cotton varieties (lines). Acta Agron Sin, 2013, 39: 2074-2082 (in Chinese with English abstract) [本文引用:1] [CJCR: 1.667]