* 通讯作者(Corresponding author): 黄瑞冬, E-mail:r_huang@126.com
第一作者联系方式: E-mail:wangyitao-happy@163.com
采用人工气候箱内培养皿培养, PEG-6000溶液模拟干旱胁迫环境, 在萌发期以80、120、150和175 g L-1 PEG-6000水溶液处理31个高粱品种, 旨在根据高粱品种萌发期对不同干旱胁迫程度的响应, 筛选出具有抗旱能力的高粱品种并探讨高粱萌发期抗旱性鉴定的方法。通过主成分分析法(PCA)和神经网络自组织映射(SOM)聚类分析法对各高粱品种进行抗旱性综合分析与评定。PCA结果表明, 相对芽长、相对根长和相对萌发抗旱指数载荷量最大, 将其作为萌发期高粱抗旱性筛选的主要评价指标, 并对31个高粱品种抗旱性排序。通过SOM聚类分析将31个高粱品种按抗旱性强弱分为5类, 吉杂305等4个品种为高度抗旱品种, HL5等4个品种为抗旱品种, 辽杂10号等8个品种为中等抗旱品种, 锦杂103等7个品种为干旱敏感品种, 锦杂93等8个品种为高度干旱敏感品种。研究认为, 相对芽长、相对根长和相对萌发抗旱指数等可以作为高粱品种抗旱性鉴定的重要指标; SOM聚类分析可作为品种抗旱性分类的重要方法。
PEG-6000 was used to simulate drought stress environment in artificial climate chamber for studying the effects of different concentrations of PEG-6000 (80, 120, 150, and 175 g L-1) on the response of 31 sorghum cultivars at germination stage, with the objectives of screening drought-resistant cultivars and exploring methods suitable for assessing drought-resistance of sorghum cultivars. Through principal component analysis (PCA) and Self-Organizing-Map (SOM) cluster analysis, the drought-resistance of 31 sorghum cultivars was evaluated. PCA showed that relative bud length, relative root length and drought-resistance index during germination had the maximum weight in values, and were used as the main indices for evaluating and ranking the drought-resistant capabilities of 31 sorghum cultivars. Five groups of drought-resistance capabilities of 31 sorghum cultivars were detected based on SOM cluster analysis: four cultivars including Jiza 305 were highly drought-resistant, four cultivars including HL8 were drought-resistant, eight cultivars including Liaoza 10 were medium drought-resistant, seven cultivars including Jinza 103 were drought-sensitive and eight cultivars including Jinza 93 were highly drought-sensitive. In conclusion, relative bud length, relative root length and germination drought-resistance index can be used as indices for identifying drought-resistance in sorghum, and SOM cluster analysis is an important method for classifying drought-resistance.
随着全球气候的变化, 干旱成为农业生产中经常遇到的问题, 是影响农业可持续发展的重要因素[ 1, 2]。充分利用抗旱性强的作物和品种以减少干旱所带来的不利影响, 对促进我国干旱、半干旱地区农业生产具有重要意义[ 3, 4, 5]。高粱是我国北方干旱、半干旱地区的重要作物之一。尽管其综合抗旱性较强, 但不同品种间抗旱性存在较大差异[ 6]。种子萌发阶段是作物能否在干旱条件下完成生育周期的关键时期之一, 萌发期抗旱性的强弱影响出苗好坏, 进而影响作物的生长发育。近年来, 许多学者运用萌发胁迫指数法、灰色关联度法、隶属函数法和主成分分析法(PCA)等[ 7, 8, 9]对抗旱性进行了大量的研究, 但多数只侧重单一干旱处理, 具有一定的局限性, 不能反映出多水平干旱条件下作物的综合萌发状况。而近年来兴起的自组织映射(SOM)聚类分析可以较好
地从多水平、多因素着手, 建立多水平多观察值之间复杂的映射关系, 取得很好的分析效果, 已被广泛应用于各种聚类分析领域[ 10, 11]。本试验以4个不同PEG-6000溶液浓度模拟干旱胁迫, 运用PCA分析和SOM聚类分析法对31个高粱品种的萌发性状综合分析及评定, 为高粱抗旱的进一步研究和生产应用提供科学依据。
从每品种选取25粒大小均匀一致、饱满的种子, 用7%的次氯酸钠溶液消毒10 min, 再用蒸馏水冲洗5次, 用滤纸吸干后摆放在直径9 cm并铺有2层滤纸的培养皿中, 种子间距离为粒长1~2倍。干旱胁迫液质量浓度分别为80、120、150和175 g L-1的PEG-6000水溶液, 与之对应的溶液水势分别为-0.1、-0.2、-0.3和-0.4 MPa[ 12], 对照以蒸馏水替代PEG-6000溶液, 每个培养皿中加入胁迫液8 mL, 每处理3次重复。将培养皿置RXZ-1000B型人工气候箱, 其昼/夜温度为28℃/25℃, 湿度为60%, 光照/黑暗为12 h/12 h, 光照强度为134 μmol m-2 s-1, 连续培养10 d。每2 d补充适量蒸馏水, 以保持胁迫液的水势不变。在第2、第4、第6、第8和第10天调查发芽数, 并于第10天测定芽长、根长、芽干重、根干重和种子剩余干重。
发芽势=第4天发芽种子粒数/供试种子粒数×100; 发芽率=第10天发芽种子粒数/供试种子粒数×100; 萌发抗旱指数(DGRI)=干旱胁迫下种子萌发指数(PIS)/对照种子萌发指数(PIC), 式中PIS/PIC = (1.00) nd2+(0.75) nd4+(0.50) nd6+(0.25) nd8, 其中, nd2、 nd4、 nd6和 nd8分别为第2、第4、第6和第8天的萌发率[ 13]; 第10天, 从各处理随机取出高粱幼苗5株, 测定芽长和根长, 于85℃烘干后称量芽干重、根干重和籽粒剩余干重, 计算干物质转运率[ 14]。
由表2可知, 相对芽长、相对根长和相对萌发抗旱指数等8个测定指标在干旱胁迫下总体呈现不同程度的下降, 随着胁迫程度的加大, 测定指标较对照降幅增大。吉杂127 (S13)的芽长在-0.4 MPa水势下降幅最大, 相对芽长仅为0.33, 而处于同水势的HL5 (S19)为0.76, 降幅最小。龙杂9号(S21)在-0.1 MPa水势下具有较高的根长, 相对根长为1.21, 辽杂35 (S27)在-0.4 MPa水势下根长降幅最小, 相对根长为0.94。在-0.4 MPa水势下, 除吉杂127和龙杂11 (S28)相对发芽势分别为0.20和0.33外, 大多数品种相对发芽势在0.65~0.85之间, 平均值为0.74。4个水势下, 多数品种的相对发芽率在0.80~1.00之间, 龙杂9号在-0.1 MPa水势下相对发芽率最高, 为1.15, 吉杂127在-0.4 MPa水势下的相对发芽率最低, 为0.47。在-0.4 MPa水势下, 吉杂127相对芽干重和根干重均最低, 分别比对照下降65.0%和68.0%; 吉杂124 (S12)具有最高的相对芽干重, 龙杂10号(S22)具有最高的相对根干重, 二者分别仅比对照降低18.0%和8.0%。辽杂15 (S15)在-0.2 MPa水势下具有最高的相对干物质转运率。吉杂127在-0.4 MPa水势下的相对干物质转运率最低, 仅为0.26。HL8 (S20)在-0.1 MPa水势下的相对抗旱萌发指数最高, 为1.11, 而吉杂127最低, 为0.33, 其次为龙杂11 (S28), 其相对萌发抗旱指数比对照降低53.0%。上述结果表明, 低强度干旱胁迫(-0.1 MPa)对高粱种子萌发影响较小或略有提升活力作用; 高强度干旱胁迫(-0.34 ~ -0.40 MPa)足以使高粱不同抗旱性得到表现。
由表3可知, 萌发期高粱品种间、同品种不同干旱胁迫条件下测定值指标间差异显著, 并且由此引起的互作差异显著, 说明干旱胁迫对高粱萌发性状影响明显, 不同基因型间存在差异。
对31个高粱品种在4个胁迫浓度下的测定值均赋相同权重, 结果如表4所示, 各萌发性状测定指标在干旱胁迫条件下几乎都小于对照, 不同品种降幅不同。相对芽长变化的范围为0.55~0.91, 相对根长的变化范围为0.58~0.99, 相对发芽势的变化范围为0.62~1.08, 相对发芽率的变化范围为0.83~1.12, 相对芽干重的变化范围为0.55~0.90, 相对根干重的变化范围为0.61~1.08, 相对干物质转运率的变化范围为0.51~0.86, 相对萌发抗旱指数的变化范围为0.63~1.03。吉杂305 (S14)和HL5 (S19)的芽长、根长、发芽势、发芽率、芽干重、根干重、干物质转运率和萌发抗旱指数相对值较大, 表现出较强的抗旱能力。吉杂123 (S11)和吉杂124 (S12)各项指标的相对值比较小, 对干旱胁迫比较敏感, 抗旱能力相对较差。干物质转运率可以反映出各品种对种子中贮藏物质的利用速率与效率以及体内能量供应水平[ 17]。龙杂11 (S28)的相对干物质转运率最高, 为0.86, 对干旱胁迫不敏感, 吉杂319 (S18)的相对干物质转运率最低, 为0.51, 对干旱胁迫较敏感。在干旱胁迫下, 萌发抗旱指数可避免因品种本身所具有的差异对试验结果的影响, 较客观地鉴定和评价品种对干旱胁迫的各种反应[ 18]。吉杂305 (S14)相对萌发抗旱指数最高, 为1.03, 而龙杂11 (S28)的相对萌发抗旱指数最低, 为0.63。
利用双变量Pearson简单相关系数法对干旱胁迫下萌发期8个指标性状的相对值进行相关性分析(表5)[ 19], 各性状之间均呈现正相关, 多数性状间的相关性达到了显著或极显著水平。相对萌发抗旱指数与相对发芽势和相对发芽率相关系数最大, 分别达到0.844**和0.822**。相对芽长与相对根长、相对芽干重和相对干物质转运率极显著正相关, 相关系数分别为0.494**、0.810**和0.716**。相对根长与相对芽干重和相对根干重也呈极显著正相关, 相关系数分别达到0.577**和0.731**。
主成分分析是利用降维的思想, 描述隐藏在一组被测变量中无法直接测量到的隐性变量, 是综合处理多因素问题的有效方法[ 20]。本试验中, 不同胁迫强度下测得数据经加权平均后进行主成分分析, 所得结果如表6所示, 第i主成分的贡献率为47.674%。从主成分表达式可以看出, 在公式(1)中第i主成分上相对芽长和相对芽干重都有较高的载荷, 分别为0.429和0.385, 说明第i主成分由变量相对芽长决定的, 可以初步认定第i主成分反映高粱在干旱胁迫下的“芽部”萌发情况, 因此, 可以把第i主成分称之为“芽部性状因子”; 第ii主成分的贡献率为27.625%, 从公式(2)中可以看出, 第ii主成分在相对发芽势、相对发芽率和相对萌发抗旱指数上有较大的载荷, 反映的是高粱种子的萌发状况, 因此, 可以把第ii主成分称之为“萌发因子”; 第iii主成分的贡献率为12.457%, 在相对根长和相对根干重上有较大的载荷, 主要反映的是根部的生长发育状况, 因此, 可以把第iii主成分称之为“根部性状因子”。
前3个主成分的累计贡献率达87.756%, 三者足以说明该数据的变化趋势, 完全符合主成分分析的要求, 故取前3个主成分可以作为数据分析的有效成分。
经过主成分分析, 分离出的相互独立的8个主成分累计贡献率已达到100%, X1~ X8即可完全代表变量的所有变异程度。前3个主成分的累计贡献率达到87.756%, 可代表变量的绝大多数信息, 为了方便运算及形象展示, 可用前3个主成分描述变异。根据各品种的前3个主成分的因子载荷量绘图, 得到图1。各项指标由于不同的主成分载荷值而分布在不同的空间位置, 其中相对发芽势、相对发芽率和相对萌发抗旱指数之间的欧式距离较近, 说明他们之间在不同浓度处理下具有相似的响应。同样, 相对根长和相对根干重在不同浓度处理下具有相似的响应, 相对芽长、相对芽干重和相对干物质转运率在不同浓度处理下具有相似的响应。并且8个测定指标均不落在原点和坐标轴上, 说明这8个指标与这3个坐标轴所对应的因子都有关系。
表7和表8分别为各成分的因子载荷矩阵和成分向量矩阵, 据此可获得各因子得分公式。
Y1=0.429 X1+0.381 X2+0.327 X3+0.343 X4+0.385 X5+0.307 X6+0.354 X7+0.281 X8 (1)
Y2=-0.149 X1-0.158 X2+0.460 X3+0.426 X4-0.352 X5-0.207 X6-0.340 X7+0.527 X8 (2)
Y3=-0.345 X1+0.506 X2+0.018 X3-0.065 X4-0.268 X5+0.662 X6-0.329 X7-0.041×X8 (3)
根据各品种 Y1、 Y2和 Y3的值, 各品种的综合得分公式为: Y=0.47674 Y1+0.27625 Y2+0.12457 Y3 (4)
计算后 Y值见表9, Y值越大, 抗旱性越强。因此, 高粱品种萌发期抗旱性为S14> S16> S20> S21> S24> S19> S25> S23> S22> S1> S2> S3> S30> S6> S9> S27> S31> S26> S10> S15> S5> S29> S8> S17> S7> S4> S28> S18> S11> S12> S13。
将31个高粱品种在4个不同PEG-6000溶液浓度下测得的8个性状的相对值数据标准化后, 利用数据挖掘软件MATLAB进行自组织映射(SOM)聚类分析[ 21], 将标准化的数据作为网络的输入层, 进行100 000次的训练后稳定, 得到聚类分析结果如表10所示。
聚类分析结果表明, 31个高粱品种可以分为5种类型, 同时根据主成分分析结果判断, 从第I类到第V类基本遵循抗旱性由弱到强排列。第I类高粱品种抗旱性最弱, 属于高度干旱敏感品种, 含锦杂93 (S4)、四杂25 (S7)和吉杂123 (S11)等8个品种; 第II类高粱品种抗旱性较弱, 属于干旱敏感品种, 含锦杂103 (S5)、四杂29 (S8)和吉杂97 (S10)等7个品种; 第III类高粱品种抗旱性一般, 属于中等抗旱性品种, 含辽杂10号(S1)、辽杂12 (S2)和辽杂21(S3)等8个品种; 第IV类高粱品种抗旱性较强, 属于抗旱性品种, 含HL5(S19)、龙杂10号(S22)和辽杂11 (S24)等4个品种; 第V类高粱品种抗旱性最强, 属于高度抗旱性品种, 含吉杂305 (S14)、吉杂99 (S16)和HL8 (S20)等4个品种。
植物在长期进化中形成各种抗旱机能, 并表现出不同生育阶段的差异性。萌发期是植物可能遇到的首个干旱胁迫时期, 是实现苗全、苗齐、苗壮的关键时期。通过简便、快速而准确的萌发期抗旱性鉴定分析方法, 筛选出具有萌发期抗旱性的作物品种, 具有重要理论意义和实践价值。多种作物上的抗旱性研究都表明, 作物抗旱性受多种因素影响, 以单一指标评价高粱萌发期的抗旱能力具有一定的片面性, 而多指标评价又过于繁琐和低效, 通过其中的几个主要因素进行分析则可达到快速、准确的目的[ 22, 23]。因此, 本研究利用了主成分和SOM聚类分析的方法, 对31个高粱品种萌发期抗旱性进行了鉴定与分类, 并取得了比较符合实际的研究结果。主成分分析法中, 将4个胁迫浓度下测得的数据均设权重为1, 在计算中涉及到品种或测定指标的相关系数的情况下, 可能会忽略一些高粱品种在不同胁迫程度下的特征差异, 但从抗旱性的排序和分类结果上看, 与试验中的观察结果较为一致, 能够有效地对31个高粱品种的抗旱性进行评定。
本试验中, 干旱胁迫下, 大多数高粱种子芽长、根长、发芽势、发芽率、芽干重、根干重、干物质转运率和萌发抗旱指数都不同程度地下降, 说明干旱胁迫抑制了高粱萌发的生理过程, 但轻度胁迫条件下个别品种某些测定指标较对照有所升高, 可能是一定程度上干旱胁迫促进了种子的某些萌发进程, 这与前人研究的适当逆境胁迫可能会提高种子发芽率和活力具有相似之处[ 24, 25]。作物的抗旱性主要表现在对逆境的适应性和敏感性上, 强大的根系是其适应干旱胁迫的重要特征之一。本试验结果表明, 在萌发期干旱胁迫条件下, 根的生长状况与芽长显著正相关。通过主成分分析可知, 芽部性状因子对高粱萌发期抗旱性影响最大, 因此, 在田间鉴定不便于测定根系时, 也可主要通过对芽部性状因子的研究评价高粱的抗旱性。
自组织映射(SOM)神经网络是无导师自组织和自学习的网络, 具有自稳定性。SOM聚类过程无需外界定义出评价函数, 能够识别出向量空间中最有意义的特征, 将高维数据压缩到二维平面, 可以广泛应用于数据的聚类分析, 并且这一过程无需人为干涉赋予各因子权重值, 从而简化分类过程, 得到较稳定的分类结果。本试验中, SOM分类结果与主成分分析结果较为一致, 说明SOM应用于高粱抗旱性分类是可行的。SOM的缺点是压缩过程中难免有数据信息损失[ 26, 27, 28, 29]。与单一胁迫水平的聚类分析不同, 本研究中聚类分析全面考虑4个胁迫水平上31个高粱品种萌发特征的综合表现, 涉及到的数据量很大, 少量的信息损失不影响结果判定的准确性。实际应用中, 如何避免神经网络训练时出现“死神经元”造成数据损失, 以及如何在神经网络训练时确定较少的训练次数以节约时间并得出较为稳定的输出结果, 需要试验人员根据实际情况判别和遴选, 但是SOM神经网络聚类分析不失为种质资源抗旱性鉴定的一个较为简便、理想的方法。
通过主成分分析和SOM聚类分析, 将31个高粱品种按照抗旱性强弱分为5大类, 吉杂305和吉杂99等4个品种为高度抗旱品种, HL5和龙杂10号等4个品种为抗旱品种, 辽杂10号和辽杂12等8个品种为中等抗旱品种, 锦杂103和四杂29等7个品种为干旱敏感品种, 锦杂93和四杂25等8个品种为高度干旱敏感品种。相对芽长、相对根长和相对萌发抗旱指数可以作为高粱品种抗旱性鉴定的重要指标。