第43卷第3期2021年7月
Vol.43,No.3
Jul.,2021中国糖料
Sugar Crops of China
doi：10.13570/jki.scc.2021.03.001
利用多元统计分析方法评价甘蔗新品系的抗旱性
蒋洪涛1，2，蒋夏兰1，商贺阳1，李安贵1，张木清1传闻中的七公主主题曲
（1.广西大学广西甘蔗生物学重点实验室，南宁530004；
2.福建农林大学国家甘蔗工程技术研究中心，福州350002）
摘要：评价不同甘蔗品种（系）的抗旱性，为广西不同种植地区提供适宜的耐旱品种和简单高效的鉴定方法。采用桶栽模拟干旱胁迫试验，利用主成分分析法、隶属函数法、聚类分析法和判别分析法，对17个‘GUC’系列甘蔗品系和两个对照品种‘新台糖22’（CK1）、‘福农39’（CK2）的6个抗旱性指标的抗旱性
系数进行综合评价。结果表明：水分胁迫下，不同甘蔗品种（系）的6个抗旱性指标变化程度不同，其中株高（Ph）、相对含水量（RWC）、叶绿素（Chl）含量和脯氨酸（Pro）含量4个指标相互间呈显著正相关，质膜透性（PMP）与丙二醛（MDA）含量呈显著正相关，Ph、RWC、Chl和Pro含量与PMP呈显著负相关；4个主因子代表6个抗旱指标94.12%。据隶属函数综合评价值，品系‘GUC31’、‘GUC17’、‘GUC23’和‘GUC9’的抗旱性高于CK1；其余品系低于CK1但高于CK2。聚类和判别分析将19个甘蔗品种（系）分为3类，其中‘GUC23’、‘GUC21’、‘GUC31’、‘GUC17’和‘GUC8’抗旱性较强；
‘GUC1’、‘GUC16’和‘福农39’抗旱性较差。试验中所用的生理指标和分析方法在甘蔗抗旱性评价上简单而有效，可用于未来甘蔗的抗旱性评价。
关键词：甘蔗；新品系；抗旱性；主成分分析；隶属函数；聚类分析；判别分析
中图分类号：S566.1文献标识码：A文章编号：1007-2624（2021）03-0001-07蒋洪涛，蒋夏兰，商贺阳，等.利用多元统计分析方法评价甘蔗新品系的抗旱性[J].中国糖料，2021，43(3)：1-7.
JIANG Hongtao，JIANG Xialan，SHANG Heyang，et al.Evaluation on drought-tolerance in new sugarcane clones by multiple statistics analysis[J].Sugar Crops of China,2021,43(3):1-7.
0引言
近几年随着其他作物面积的扩张，迫使甘蔗种植区域不断向旱地和坡地转移，中国旱地甘蔗所占比例达到85%以上。水资源缺乏已成为制约中国蔗糖生产的最关键因素之一，选育抗旱性强的甘蔗新品种是一种重要的经济、有效抗旱措施。抗旱性强的甘蔗新品种不仅可以保证高产稳产，而且对提高水资源利用率具有十分重要的意义。但是传统的育种周期较长，很难得到综合性能好的品种，通过引种能在较短的时间内获得适应当地推广的优良品种，而准确鉴定甘蔗品种的抗旱性首先要建立有效的甘蔗抗旱性评价体系。
作物的抗旱性与作物形态、生理生化性状密切相关。国内外学者从不同角度研究了作物抗旱性的鉴定方法和指标[1-3]。卢会文等[4]研究叶片衰老情况、株高伤害度、叶面积伤害度等形态指标和4个生理指标与甘蔗品种抗旱性的关系，指出农艺性状也是适合甘蔗抗旱性评价的指标；高三基等[5-7]研究水分胁迫下甘蔗叶绿素及叶绿素荧光与甘蔗抗旱性的关系，指出光合指标与品种抗旱性有密切关系；张木清等[8]对水分胁迫下蔗叶活性氧代谢回归分析表明，水分胁迫下丙二醛含量与质膜透性变化呈正相关，可作为鉴定甘蔗品种抗旱性较为可靠的生理生化指标。研究表明，多种作物在干旱胁迫下脯氨酸快速积累以抵御干旱胁迫的伤害[9-12]，所以该指标也可作为评价甘蔗抗旱的重要指标之一。不同植物的抗旱能力和机制不同，利用多个单一指标对作物进行抗旱性评判容易造成评价结果的片面性和不稳定性。
自制天然染发剂主成分分析可以快速准确地确定各指标权重，找出数目较少且能代表所有变量的主成分。主成分分析
对原有信息全面覆盖，同步解决了综合变量覆盖和降维问题，降低了数据对象集合的相异度和维度，实现了收稿日期：2020-04-02
基金项目：财政部和农业农村部国家现代农业（糖料）产业技术体系“抗病品种选育”（CARS170190）；福建农林大学国家甘蔗工程技术研究中心开放课题（266-PTJH1500102）资助。
第一作者：蒋洪涛（1985-），男，广西全州人，硕士研究生，研究方向为作物栽培学与耕作学，E-mail:******************。
通信作者：张木清（1966-），男，福建福州人，教授，博士生导师，研究方向甘蔗抗旱抗病性育种，E-mail:***************。
zgti.cbptki
2
中国糖料zgti.cbptki2021
数据对象集合的聚类归约[13]。研究表明，甘蔗叶片抗氧化物酶活性和渗透调节物质含量的差异是品
种间耐干旱胁迫差异的生理基础[14]。目前，利用主成分分析并结合隶属函数、聚类分析的方法已经广泛地用于多种植物的抗逆性评价。田治国等[15]对9个万寿菊品种16个生理指标的抗旱性进行综合评价，结果表明，不同品种万寿菊的生理指标变化程度不同，通过主成分分析法、隶属函数法和聚类分析法能够很好地评价万寿菊的抗旱性；王兰芬等[16]采用抗旱系数、隶属函数和抗旱指数法评价绿豆的抗旱性，结果表明，3种评价方法抗旱分级结果一致性较高，说明在抗旱性评价上抗旱系数和隶属函数法、抗旱指数法都具有很好的评价效果。边芯等[17]通过利用抗旱系数和隶属函数法对持绿性‘割手密’的抗旱性进行综合评价，探明了持绿性‘割手密’的抗旱性遗传特性，为持绿性甘蔗细茎野生种的利用提供依据；田春艳等[18]利用综合分析方法进一步明确了质膜透性（PMP）、丙二醛（MDA）含量、脯氨酸（Pro）含量、自由水/束缚水等是甘蔗抗旱性评价的关键指标，可有效评价甘蔗种质抗旱性；王天菊等[19]利用相关分析、聚类分析和模糊隶属函数分析法对26份‘割手密’的抗旱性进行评价，为‘割手密’抗旱材料筛选及抗旱基因发掘奠定基础。王焱等[20]利用相关分析、隶属函数系数、综合抗旱性系数、灰色关联分析、逐步回归和聚类分析方法对苜蓿发芽期的抗旱性进行评价，筛选出在萌发期表现出较强的抗旱性的3个苜蓿基因型。王竞红等[21]对3种观赏草、李怡等[22]对54份羊草种质资源、孙佩等[23]对杨树苗期、路之娟等[24]对不同基因型苦荞采用以上相关方法进行了抗旱性综合评价。然而目前对于引进的甘蔗品种（系）还没有抗旱性差异综合性评价。因此，本文对19个潜力较大的甘蔗品种（系）的6个抗旱性指标进行主成分分析、隶属函数分析和聚类及判别分析，从而对供试材料的抗旱性进行综合评价，旨在筛选出抗旱性强的品种（系），为今后的甘蔗育种和引种提供理论依据。建造师考试科目
1材料与方法
1.1试验材料
参试甘蔗材料为2012年广西大学在美国农业部运河点育种站按照常规杂交育种技术选择5个甘蔗组合，分别为‘HoCP01-157’בCP14-0969’、‘CP49-50’בCP96-1252’、‘CP88-1762’בCP96-1252’、‘CP00-1100’בQ209’、‘CP89-2143’בCP72-1210’；杂交所得9个甘蔗花穗，带回国后经过杂种圃和选种圃后选择出17个糖分、产量和抗逆性较好的甘蔗材料，其中‘GUC1’、‘GUC8’、‘GUC10’、‘GUC16’、‘GUC17’和‘GUC23’是‘HoCP01-157’בCP14-0969’杂交后代；‘GUC19’、‘GUC21’、‘GUC29’、‘GUC31’和‘GUC33’是‘CP49-50’בCP96-1252’杂交后代；‘GUC2’、‘GUC3’和‘GUC24’是‘CP88-1762’בCP96-1252’杂交后代；‘GUC9’和‘GUC13’是‘CP00-1100’בQ209’杂交后代；‘GUC25’是‘CP89-2143’בCP72-1210’杂交后代。在对17个甘蔗材料进行品比的同时进行抗旱性评价试验，试验对照甘蔗品种为广西本地主栽品种‘新台糖22号’（"ROC"22）和‘福农39’（‘FN39’）。
1.2试验设计
2015年7月20日，将所有试验材料砍成单芽，用多菌灵浸种处理后种植于塑料桶中，每桶装等量的种植土（种植土以沙质土壤与有机肥3:1的比例混合），每个品种种植10桶，每桶下种5～6芽，正常肥
水管理，保持所有材料的正常生长，待成苗后每桶定植3株。分蘖期移至网室大棚，每个品种选取长势均匀9桶作为研究对象，其中6桶为试验组，3桶作为对照组，试验组于10月28日开始停止供水处理，土壤相对含水量为81.9%～86.0%，停止供水20d后，自然失水到部分品种叶片干枯萎蔫，清晨无“吐水”现象时，此时土壤相对含水量为23.8%～31.2%；对照组则一直控制土壤相对含水量在81.9%～86.0%。并于早上8:00—9:00取样，每个处理组和对照组均采集6株甘蔗的+1叶混合为1份样品，放入冰壶中，带回实验室后将试验组和对照组样品分别进行随机分组，每2张叶片一组，形成3个生物学重复，进行株高（Ph）、相对含水量（RWC）、叶绿素（Chl）、丙二醛（MDA）、膜透性（PMP）、脯氨酸（Pro）含量6个生理指标的测定。
1.3测定项目与方法
株高（Ph）测定，选择成熟挺直、无病害、长势一致、具代表性的甘蔗主茎，用皮尺从地面测量至甘蔗+1叶。甘蔗叶片相对含水量（RWC）采用重量法测定，RWC=（叶片鲜重－干重）/（水饱和重－干重）×100%[25]；质膜透性（PMP）测定参照侯冬花[26]方法，用相对电导率表示，PMP=处理电导率/煮沸电导率×100%；丙二醛含量（MDA，μmol/g）测定参照HEATH等[27]方法；叶绿素含量（Chl，mg/g）测定采用混合液法[28]；脯氨酸含量（Pro，mg/g）测定参照张瑛等[29]方法。
蒋洪涛，等：利用多元统计分析方法评价甘蔗新品系的抗旱性
第43卷，第3期参照BOUSLAMA 等[30]的方法计算Ph 、RWC 、PMP 、MDA 、Chl 和Pro 含量各指标性状的抗旱系数：抗旱系数=干旱处理测定值/对照测定值×100%1.4
数据处理
采用SPSS 22.0软件进行主成分分析、隶属函数分析和判别分析，用Origin 2019软件进行聚类分析，采用欧式距离法进行聚类分析，运用的主要公式如下。1.4.1
隶属函数值
U (x i )=
x i -x min
x max -x min
kettle
（1）
式中：i =1，2，3，…，19（下同）；x i 表示第i 个因子的得分值，x min x max 分别表示第i 个因子的最小值和最大值。1.4.2
权重
W i =P i /∑i =1
n P i
（2）
式中：W i 表示第i 个公因子在所有公因子中的重要程度，P i 为各品种第i 个公因子的贡献率。1.4.3
综合评价值
学化妆盘发D i =∑i =1
n [U (x i )×W i ]
（3）
式中：D i 值表示甘蔗材料抗旱性综合评价值，U (x i )表示各品种隶属度值，W i 表示各品种权重值。2
结果与分析
2.1
各品种（系）的抗旱系数和相关性分析
干旱胁迫下，不同甘蔗品种（系）的6项抗旱相关指标呈现不同的变化（表1）。部分甘蔗在持续水分胁迫
下能保持一定的生长能力，也有少数的甘蔗MDA 含量要低于正常供水的含量；但总体来说，RWC 和Chl 含量表现为明显的下降，PMP 、MDA 含量和Pro 含量有较大幅度的增加，反映出干旱胁迫在一定程度上造成了甘蔗膜脂过氧化；而不同品种甘蔗对水分胁迫的响应表现不一样。相关分析表明，19个甘蔗品种（系）的部分生理指标之间存在显著的相关性。其中Ph 、RWC 、Chl 和Pro 与PMP 呈显著负相关；除Ph 以外，与MDA 也呈显著负相关；Ph 、RWC 、Chl 和Pro 四个指标间相互呈显著正相关；PMP 与MDA 呈显著正相关（表2）。说明，在干旱胁迫下，不同甘蔗品种的各抗旱指标之间存在一定的协同性和拮抗性，在减少水分胁迫对膜系统
表119个甘蔗品种（系）各指标的抗旱系数
Table 1
wimax是什么Drought-tolerant coefficients of different parameters of nineteen sugarcane cultivars (clones)
品种Cultivar 株高/cm Ph 相对含水量/%RWC 叶绿素含量/(mg/g)Chl 丙二醛含量/(mol/g)MDA 膜透性/%PMP 脯氨酸含量/(mg/g)Pro GUC193.2234.9454.17190.01417.901395.23GUC2103.3778.8962.99122.57344.941725.71GUC3101.5873.0566.58116.13314.611849.57GUC8112.2086.0186.1493.69199.731741.74GUC9105.9079.1481.98165.43287.401384.64GUC1098.3658.3252.40152.67364.261585.84GUC1387.6267.4654.56159.47309.931324.07GUC16100.2556.9361.10119.40272.091536.86GUC17107.6585.1483.78123.27130.001845.73GUC19109.9345.8872.31132.75375.681558.94GUC21111.9181.9681.4599.86176.682027.24GUC23111.4291.6278.99109.64176.832230.81GUC24101.0988.6561.45131.78320.471922.77GUC2596.3959.0578.99138.69250.851655.84GUC29104.8269.4283.03130.79344.771656.29GUC31112.3374.3396.86132.39186.151842.02GUC3399.6040.4121.19231.15482.411515.72"ROC"22122.1374.2075.75145.85283.251609.84FN3994.1038.0845.24109.59428.091131.53
resignNote:Ph,plant height;RWC,Relative water content;Chl,Chlorophylls;MDA,malondialdehyde;PMP,Membrane permeability;Pro,Proline,and the
same hereinafter.
3
中国糖料zgti.cbptki 2021
伤害的同时保持叶片较高的相对含水量、叶绿素含量和脯氨酸的含量，以提高甘蔗的抗旱能力，从而促进甘蔗的生长。
坏了2.2各品种（系）的主成分分析
6个抗旱指标的抗旱系数转化为4个主因子，累加贡献率达94.12%（表3），说明前4个主因子所包含的要
素信息量可以反映出6个抗旱性指标原始特征参数的大部分信息，既保留了绝大部分信息，又达到了降维的目的。主因子1中，Chl 和PMP 指标对甘蔗品种抗旱性起支配作用，反映了甘蔗抗旱性与叶片光合作用和膜透性密切相关；因子2中，RWC 和Pro 指标对甘蔗品种抗旱性起支配作用，表明渗透调节对于甘蔗抗旱至关重要；因子3中，只有Ph 具有较大载荷，反映了甘蔗抗旱性强弱直接影响甘蔗的生长；因子4中，只有MDA 含量有较大载荷，表明MDA 含量是相对较为独立的影响因子。甘蔗各生理指标的抗旱系数与其主因子系数相结合，得到每个品种（系）的4个主因子得分值CI(x i )。CI(x i )值越大，说明该品种（系）在这一主因子上的抗旱性越强，反之抗旱性越弱。由表4可以看出，‘GUC9’、‘GUC17’和‘GUC31’在Chl 和PMP 上抗旱性较强；‘GUC23’和‘GUC24’在RWC 和Pro 上的抗旱性较强；‘GUC19’和‘"ROC"22’在Ph 上的抗旱性较强；‘GUC1’、‘GUC9’和‘GUC33’在MDA 上的抗旱性较强。各CI(x i )值间的相关系数可以反映出一个因子作用性状与其他因子作用性状的间接作用。
2.3各品种（系）的隶属函数及综合评价值
依据主因子贡献率的大小，可以计算出各主因子的权重（IW ）。4个主因子的权重值不等，主因子1的权
重最大（0.307）；主因子4的权重最小（0.195）。由于4个主因子在不同甘蔗品种（系）抗旱性评价中的作用不
表3各主因子系数及贡献率
Table 3
Coefficients of principal factors and contribution rate
项目
Item
主因子1Factor 1
主因子2Factor 2主因子3Factor 3主因子4Factor 4株高Ph 0.2710.2790.907-0.117相对含水量RWC 0.5330.7200.153-0.203叶绿素Chl 0.8170.1690.389-0.303丙二醛MDA -0.300-0.225-0.1220.918膜透性PMP -0.769-0.472-0.1640.280脯氨酸Pro 0.1740.8860.291-0.190贡献率/%CR 28.9328.0518.7418.40累计贡献率/%CCR 28.93
alarm什么意思56.98
75.72
94.12
Note:CR,Contribution rate;CCR,Cumulative contribution rate
表2抗旱性指标的相关系数矩阵
Table 2
Correlative coefficient matrix of drought-tolerant indexes
项目Item Ph RWC Chl MDA PMP Pro 株高Ph
1.000
0.531**0.630**-0.365-0.532**0.554**相对含水量RWC 1.0000.651**-0.538**-0.773**0.731**叶绿素Chl 1.000**
-0.603**-0.814**0.506*丙二醛MDA 1.000
0.614**-0.453*膜透性PMP 1.000**
-0.662**脯氨酸Pro
1.000
表419个甘蔗品种（系）的4个主因子得分值
Table 4Principal factor scores of nineteen sugarcane cultivars (clones)
初三英语课文翻译品种Cultivar CI(1)CI(2)CI(3)CI(4)品种Cultivar CI(1)CI(2)CI(3)CI(4)GUC1-0.40-1.08-0.60 1.22GUC210.21 1.040.56-0.86GUC2-0.700.61-0.15-0.55GUC23-0.05 2.100.35-0.39GUC3-0.650.78-0.36-0.78GUC24-0.67 1.66-0.68-0.05GUC80.770.060.59-1.06GUC250.97-0.35-0.980.09GUC9 1.60-0.840.12 1.33GUC290.24-0.390.34-0.21GUC10-0.770.03-0.490.21GUC31 1.43-0.050.890.37GUC130.77-0.37-2.320.64GUC33-2.210.330.24 2.52GUC16-0.11-0.52-0.49-0.81"ROC"220.09-0.47 2.400.58GUC17 1.440.67-0.220.15FN39-1.26-1.99-0.67-1.87
GUC19
-0.71-1.22 1.47-0.52
4
蒋洪涛，等：利用多元统计分析方法评价甘蔗新品系的抗旱性
第43卷，第3期同，因此，在CI(x i )值的基础上利用隶属函数法进行数值转换，得到各品种的隶属函数值U(x i )和综合评价值（D i ），根据D i 值评价各品种（系）的抗旱性。由表5可知，19个甘蔗品种（系）的抗旱性强弱依次为‘GUC31’>‘GUC17’>‘GUC23’>‘GUC9’>‘"ROC"22’>‘GUC21’>‘GUC8’>‘GUC24’>‘GUC25’>‘GUC29’>‘GUC33’>‘GUC13’>‘GUC2’>‘GUC3’>‘GUC10’>‘GUC1’>‘GUC16’>‘GUC19’>‘福农39’。
2.4甘蔗不同种质材料抗旱性的聚类及判别分析
聚类分析（图1）在欧式距离5.08时，将19个甘蔗品种（系）分为3类：高抗品种（系）5个，包括‘GUC23’、
‘GUC21’、‘GUC31’、‘GUC17’和‘GUC8’；中抗品种（系）10个，包括‘"ROC"22’、‘GUC9’、‘GUC19’、‘GUC25’、‘GUC16’、‘GUC10’、‘GUC29’、‘GUC24’、‘GUC3’和‘GUC2’；不抗旱品种（系）4个，包括‘GUC33’、‘福农39’、‘GUC13’和‘GUC1’。将对甘蔗抗旱性（y ）归类有显著影响的6个抗旱性指标Ph （x 1）、RWC （x 2）、Chl （x 3）、MDA （x 4）、PMP （x 5）和Pro （x 6）数据标准化后作为判别式的变量，进行逐步判别分析，经显著性测验，所有显著指标都入选，建立如下判别函数：
y 高=3.625x 1＋1.037x 2＋3.527x 3＋0.669x 4－0.350x 5－0.968x 6－3.544y 中=-2.663x 1＋0.254
Rx 2－1.829x 3－0.654x 4＋0.465x 5＋0.715x 6－2.236y 低=-3.582x 1－3.701x 2－6.312x 3－0.481x 4－0.036x 5＋1.237x 6－
10.621
G U C
23
G U C 21G U C 31G U C 17G U C 8R O C 22G U C 9G U C 19G U C 25G U C 16G U C 10G U C 29G U C 24G U C 3G U C 2G U C 33福农39G U C 13G U
C 10
2
4
图119个甘蔗品种的抗旱性聚类分析图
Fig.1Cluster map of drought-tolerance in nineteen sugarcane cultivars (clones)
表519个甘蔗品种（系）的隶属函数值U (x i )和权重(IW )
Table 5U (x i )and IW in nineteen sugarcane cultivars (clones)
品种Cultivar U(1)U(2)U(3)U(4)综合评价值D
Comprehensive evaluation value D
GUC10.480.220.360.700.422GUC20.400.630.460.300.461GUC30.410.680.410.250.458GUC80.780.500.620.180.548GUC9 1.000.280.520.730.636GUC100.380.490.390.470.433GUC130.780.400.000.570.469GUC160.550.360.390.240.401GUC170.960.650.440.460.666GUC190.390.190.800.310.397GUC210.630.740.610.230.582GUC230.57 1.000.570.340.651GUC240.400.890.350.410.540GUC250.840.400.280.450.519GUC290.640.390.560.380.500GUC310.960.470.680.510.669GUC330.000.570.54
1.000.472"ROC"220.600.37 1.000.560.604FN390.250.000.350.000.146权重IW
0.307
0.298
0.199
0.195
–
5
