甘蔗Saccharum hybrids spp. 是世界上最重要的糖料作物及较有发展前景的可再生绿色能源作物。在我国，蔗糖产量占国民食糖需求总量的90%以上^[1]。但随着经济的发展，甘蔗与果蔬争地矛盾日趋突出，甘蔗生产逐渐向干旱缺水的坡地和偏远山区转移，而甘蔗的种植对自然资源有一定的要求，特别是对水分和温度的要求较严格^[2]。干旱己成为制约我国蔗糖产量及品质的关键因素之一，选育和推广抗旱性强的高产、高糖品种对提高甘蔗产量、降低生产成本具有重要意义。

在干旱胁迫条件下，作物能通过产生保护性物质降低对干旱的敏感性，或产生一系列生理生化变化，减轻干旱的伤害，表现出抗逆性^[3]。研究表明甘蔗根系活力，叶片可溶性蛋白含量、可溶性糖含量、脯氨酸含量变化可作为甘蔗抗旱性鉴定和评价的重要指标^[4-6]。在干旱胁迫条件下，植物叶片中的超氧化物歧化酶(Superoxide dismutase，SOD)、过氧化物酶(Peroxidase，POD)和过氧化氢酶(Catalase，CAT)这3种酶活性快速降低，使植物的膜系统受外界不良因素的影响而损坏，这种反应与生物氧自由基有一定关联性^[7-8]。边芯等^[3]通过研究干旱胁迫对持绿性甘蔗细茎野生种叶片生理及甘蔗生物量的影响，表明在干旱胁迫条件下持绿性甘蔗细茎野生种能通过产生保护性物质降低对干旱的敏感性，通过SOD、CAT、POD活性变化，减轻干旱对植株的伤害，从而表现出更好的抗旱性，SOD、CAT、POD活性可作为衡量甘蔗抗旱能力的指标。

甘蔗受到干旱胁迫后，产生的各种生理变化错综复杂，受多种因素影响，仅根据某单一指标评定甘蔗品种的抗旱性会造成评定结果的片面性。综合运用主成分分析法和模糊隶属函数法评价作物的抗旱性可综合评定作物主要抗旱性生理生化指标，对参试作物品种的抗旱性强弱进行排序及聚类，得到较单一指标更为客观可靠的抗旱性评定结果^[9-10]。这种方法在大豆^[9]、玉米^[10]等主要农作物的抗旱性评价中已得到较好应用，在甘蔗品种抗旱性评价中尚未广泛应用，甘蔗品种抗旱性评价以单一指标或少数指标为主。本研究以华南农业大学和湛江市农业科学研究院联合培育的甘蔗10个优良新品系和甘蔗主栽品种ROC22为材料，测定甘蔗伸长初期人工干旱胁迫和正常供水2种处理甘蔗叶片及根系的8项生理参数，运用主成分分析结合聚类分析综合评价参试甘蔗材料的抗旱性，为甘蔗优良抗旱新品种选育提供科学依据。

1 材料与方法 1.1 试验材料

参试材料为华南农业大学甘蔗研究室和湛江市农业科学研究院联合培育的10个甘蔗优良新品系：A6−13115、A3−1320、A6−13111、A13−1396、A1−1305、A11−1390、A7−13120、A6−13122、A7−13104、A4−1316，以及目前我国甘蔗主栽耐旱品种ROC22(对照)。

1.2 材料种植与处理

试验于2017年2—5月在华南农业大学甘蔗育种基地温室进行。采用盆栽种植的方法，用高、上径、下径分别为32.0、33.5、28.0 cm的塑料桶，每个桶的桶底钻3~5个直径为0.4 cm的圆孔，排去土壤中多余水分。所用栽培土壤为大田土，每桶装2/3高的土壤。蔗种砍成单芽段，每桶种4~5个单芽，每份材料种6桶，下种后覆土3~5 cm，待甘蔗成苗后每桶保留3株长势相近的蔗苗。

待甘蔗生长到伸长初期时，每个材料分成2组(每组3桶)，充分淋水至桶里土壤水分饱和。一组正常供水(对照组)，另一组不供水进行干旱胁迫处理(处理组)。干旱胁迫开始后，每天定时观察2组蔗株的生长情况，当处理组大多数蔗株正3叶以下全部枯黄，且正3叶全部萎蔫时(约干旱胁迫15 d)，甘蔗达到干旱胁迫处理要求，此时分别取对照组和处理组心叶下面第1片叶片和根系进行超低温保存，次日进行生理指标测定。

1.3 测定项目及方法

甘蔗根系活力测定采用TTC法^[11]，叶片可溶性蛋白含量测定采用考马斯亮蓝G−250比色法^[12]，叶片可溶性糖含量测定采用蒽酮比色法^[13]，叶片SOD活性测定采用氮蓝四唑光化还原法^[14]，叶片CAT活性测定采用过氧化氢还原法^[14]，叶片POD活性测定用愈创木酚法^[15]，叶片脯氨酸含量测定采用酸性茚三酮显色法^[16]，叶片丙二醛(Malondialdehyde，MDA)含量测定采用硫代巴比妥酸法^[17]。除根系活力外，其余指标的试验结果都用叶片的鲜质量表示。

1.4 数据处理

抗旱系数(DC)^[18]：

$ {\rm{DC}} \!=\! {\text{干旱胁迫处理测定值}}/{\text{对照测定值}}\! \times \! 100 \;{\text{%}}{\text{ ,}} $

(1)

综合指标的隶属函数值(U_Xj)^[19]：

$ {U_{{X_j}}} = \left( {{X_j} - {X_{\min }}} \right)/\left( {{X_{\max }} - {X_{\min }}} \right)\;\;\;\;\;\;j = 1,2, \ldots ,{{n}}, $

(2)

式中，X_j表示第j个综合指标的值，X_min表示第j个综合指标的最小值，X_max表示第j个综合指标的最大值。

综合指标的权重(W_j)^[20]：

$ {W_j} = {P_j}/\sum\limits_{j = 1}^{{n}} {{P_{{j}}}} \;\;\;\;\;j = 1,2, \ldots ,{{n}}, $

(3)

式中，W_j表示第j个综合指标在所在综合指标中的重要程度，即权重；P_j为各品系第j个综合指标的贡献率。

甘蔗品种在干旱胁迫条件下抗旱性的综合评价值用D表示^[21]：

$ D = \sum\limits_{j = 1}^{{n}} {\left( {{U_{{X_j}}} \times {{{\rm{W}}_{{j}}}} } \right)} \;\;\;\;\;j = 1,2, \ldots ,{{n}}, $

(4)

采用Microsoft Excel 2013进行数据处理和作图，使用IBM SPSS 22.0软件进行方差分析和聚类分析，聚类用欧式距离作为相似性尺度。

2 结果与分析 2.1 干旱胁迫对甘蔗根系活力的影响

由图1可知，与正常供水相比，伸长初期受到干旱胁迫后，甘蔗新品系及对照种ROC22的根系活力都呈现不同程度的降低，且降幅均达显著差异水平。降幅最小的品系是A6−13122，为34.13%；降幅最大的是A4−1316，为51.47%；对照种ROC22的降幅为35.47%。

2.2 干旱胁迫对甘蔗叶片渗透调节物质含量的影响

与正常供水相比，伸长初期受到干旱胁迫后，甘蔗新品系及对照种ROC22叶片脯氨酸含量、可溶性蛋白含量和可溶性糖含量都呈现不同程度的增加趋势，且增幅均达差异显著水平(图2)。脯氨酸含量增幅最小的品系是A7−13104，为52.50%；增幅最大的品系是A6−13111，为69.48%；对照种ROC22增幅为57.34%(图2A)。可溶性蛋白含量增幅最小的品系是A4−1316，为22.83%；增幅最大的是对照种ROC22，为36.12%(图2B)。可溶性糖含量增幅最小的品系是A4−1316，为51.13%；增幅最大的品系是A3−13104，为77.93%；对照种ROC22的增幅为75.67%(图2C)。

2.3 干旱胁迫对甘蔗叶片抗氧化酶活性的影响

与正常供水相比，伸长初期受到干旱胁迫后，甘蔗新品系及对照种ROC22叶片SOD、POD和CAT活性都呈现不同程度的增加趋势，且增幅均达差异显著水平。SOD活性增幅最小的品系是A1−1305，为34.17%；增幅最大的品系是A7−13120，为59.55%；对照种ROC22增幅为49.67%(图2D)。POD活性增幅最小的品系是A7−13120，为36.36%；增幅最大的是对照种ROC22，为64.75%(图2E)。CAT活性增幅最小的品系是A7−13120，增幅为41.82%；增幅最大的品系是A6−13115，为84.21%；对照种ROC22的增幅为71.59%(图2F)。

2.4 干旱胁迫对甘蔗叶片丙二醛含量的影响

与正常供水相比，伸长初期受到干旱胁迫后，甘蔗新品系及对照种ROC22叶片MDA含量都呈现不同程度的增加趋势，且增幅均达显著差异水平。增幅最小的品系是A7−13120，为29.46%；增幅最大的品系是A6−13115，为86.53%；对照种ROC22的增幅为70.87%(图3)。

2.5 干旱胁迫下甘蔗主要生理指标的主成分分析

对根系活力、脯氨酸含量、可溶性蛋白含量、可溶性糖含量、MDA含量、SOD活性、POD活性、CAT活性8个生理指标的抗旱系数(表1)进行主成分分析后，将8个指标转化成为相互独立的综合新指标CI₁~CI₈，前2个综合指标的贡献率分别为68.62%、12.64%，累积贡献率达到81.26%(表2)，因此忽略CI₃~CI₈这6个综合指标。由表2可以看出，在第一主成分CI₁中，根系活力、脯氨酸含量、可溶性蛋白含量、可溶性糖含量、SOD活性、POD活性、CAT活性的系数较大，占了全部信息的68.62%，因此该成分代表了根系活力、脯氨酸含量、可溶性蛋白含量、可溶性糖含量、SOD活性、POD活性、CAT活性；在第二主成分CI₂中，MDA含量的系数较大，反映了全部信息的12.64%，因此该成分代表了MDA含量。

通过公式(2)可得出2个主成分综合指标相应的隶属函数值(U_Xj)，根据2个指标的贡献率通过公式(3)求出权重(W_j)，最后通过公式(4)得到不同品系抗旱性的综合评价值(D)，并对D值进行排序(表3)。D值越大，说明甘蔗的综合抗旱能力越强。A3−1320的D值最大，为0.886，评价为最抗旱；A4−1316的D值最小，为0.115，评价为最不抗旱；其他品系介于二者之间；对照种ROC22的D值为0.603，综合抗旱能力排在第5位。

2.6 甘蔗抗旱性综合评价值的聚类分析

通过对参试甘蔗材料抗旱性综合评价值(D)进行聚类分析，将10个甘蔗新品系及对照种ROC22分成抗旱性强、中度抗旱和抗旱性弱3个类型(图4)。抗旱性强的品系有2个：A6−13115、A3−1320，所占比例为18.18%；中度抗旱的品系(品种)有3个：ROC22、A13−1396、A6−13122，所占比例为27.27%；抗旱性弱的品系有6个：A6−13111、A1−1305、A11−1390、A7−13120、A7−13104、A4−1316，所占比例为54.55%。

3 讨论与结论

抗旱性是作物在干旱胁迫下正常生长、繁殖或生存，胁迫解除后能够迅速恢复生长的能力^[22]。抗旱性鉴定是通过筛选、评价和归类，提供优异抗旱种质资源的过程^[23]。植物遭受干旱胁迫时会造成细胞膜系统的破坏和代谢紊乱，从而导致大量活性氧积累，SOD、POD、CAT是植物活性氧的清除剂，这3种酶的活性与植物抗旱性密切相关^[24-25]。本研究中，甘蔗受干旱胁迫后，与对照组(正常供水)相比，甘蔗叶片SOD、POD和CAT活性均显著增加，上述3种酶的活性增幅最小的品系依次为A1−1305、A7−13120和A7−13120，抗旱性聚类均为抗旱性弱类型，CAT活性增幅最大的品系为A6−13115，抗旱性聚类属于抗旱性强类型，验证了Devi等^[24]和Huseynova等^[25]的研究结果。作物根系生理活性与根系主动吸水有关，根系越发达，作物吸水效率越高，可缓解干旱对植株的影响^[26-27]。本研究中，甘蔗受干旱胁迫后，与正常供水相比，甘蔗根系活力均显著降低。根系活力降幅最小的品系为A6−13122，降幅最大的品系为A4−1316，抗旱性聚类表明前者属于中度抗旱类型，后者为抗旱性弱类型，这也验证了上述观点^[26-27]。王川等^[28]研究发现小麦植株受干旱胁迫后，抗旱性强的小麦品种总能保持较高的可溶性糖含量，因此认为可溶性糖含量可以作为小麦抗旱性评价的指标，本研究也得到类似结论。张明生等^[29]发现甘薯处于逆境时，与抗逆性相关的基因被激活，诱导合成新的蛋白，导致植株可溶性蛋白含量增加，增强细胞的吸水能力从而提高植株抗旱性，因此认为可溶性蛋白含量可作为甘薯抗旱性评价指标。本研究中，甘蔗受干旱胁迫后，叶片可溶性蛋白质含量均显著增加，增幅最小的品系为A4−1316，最大的为对照种ROC22，抗旱性评价聚类表明前者为抗旱性弱类型，后者为中度抗旱类型，这也表明可溶性蛋白含量与甘蔗的抗旱性相关。

参试甘蔗品系抗旱性综合评价结果表明，A6−13115抗旱性强，A7−13104、A4−1316、A1−1305和A7−13120抗旱性弱。作物抗旱性是各生理生化代谢途径对干旱胁迫响应综合作用的结果，对不同指标作物的抗旱反应不一定相同，单一指标评价作物抗旱性往往具有片面性，评价结果可靠性差，采用多个指标综合分析能使评价结果更加可靠和科学。本研究采用抗旱系数对各指标进行数据处理，使其具有可比性，结合主成分分析法可以在不损失或很少损失原有信息的前提下，将原来个数较多而且彼此相关的指标转换成少数几个新的彼此独立的综合指标，并结合其隶属函数值及权重，得到甘蔗抗旱性综合评价值D值。D值是一个无量纲的数，且位于[0, 1]，从而使得甘蔗品种的抗旱性差异具有可比性^[30]。主成分分析法在大豆、水稻、小麦的抗旱性鉴定中已得到较好的应用^{[9, 21, 23]}，因此采用主成分分析结合聚类分析的综合评价方法可提高甘蔗抗旱性鉴定结果的可靠性。本研究中，A3−1320鉴定为抗旱性强品系，A13−1396鉴定为中度抗旱品系，与田间抗旱表现一致。对照品种ROC22鉴定为中度抗旱的甘蔗品种，这与ROC22多年田间种植表现一致，但与檀小辉等^[31]对54个广西区试甘蔗品种进行抗旱性评价的结果不一致，可能是抗旱评价标准不同导致。

采用主成分和聚类方法对多个指标综合评价，能提高甘蔗抗旱评价的可靠性和准确性。本研究对参试的甘蔗新品系和对照品种ROC22作抗旱性评价，将其分为抗旱性强、中度抗旱和抗旱性弱3个类型。抗旱性强的品系有A6−13115、A3−1320；中度抗旱的品系有A13−1396、A6−13122及ROC22；抗旱性弱的品系有A6−13111、A1−1305、A11−1390、A7−13120、A7−13104、A4−1316。