论文部分内容阅读
摘要： 以7份杭白菊为试材，采用人工海水作为筛选渗透压，在测定一系列生长指标和水分生理指标的基础上，通过主成分分析、模糊数学隶属函数法、聚类分析和灰色关联等方法，对杭白菊种质的耐盐性进行综合评价，筛选出合适的评价指标。结果表明：生物量、膜透性、SOD（超氧化物歧化酶）、可溶性总糖、POD（过氧化物酶）、游离脯氨酸等指标可作为杭白菊耐盐性快速鉴定评价的指标；聚类分析法将7份杭白菊种质的耐盐性分为强、中、弱3类，耐盐性强的种质为单瓣4号，耐盐性中的种质为NJK1、管花2号、NJK2、NJK3，耐盐性弱的种质为大花4号和红心2号；模糊数学隶属函数法评价杭白菊种质耐盐性由强到弱依次为单瓣4号>NJK1>NJK2>红心2号>NJK3>管花2号>大花4号。主成分分析、隶属函数法和聚类分析对7份杭白菊种质的耐盐性评价结果相对较为一致，单瓣4号、NJK1、NJK2可作为杭白菊耐盐育种的优良种质资源。
关键词： 杭白菊；人工海水；主成分分析；综合评价；耐盐性；生长指标；生理指标
中图分类号： S682 1 10 1 文献标志码： A
文章编号：1002-1302（2015）08-0234-05
收 稿日期：2014-08-20
菊花为菊科植物菊（Chrysanthemum morifolium Ramat ）的干燥头状花序，是传统常用中药，具有散风清热、平肝明目、清热解毒的功能。近年来，随着国内外食品、医药和饮料加工业的迅猛发展，作为药食兼用的菊花用量激增。江苏省射阳县东临黄海，土地盐碱程度参差不齐，可开发沿海滩涂面积巨大。滨海盐土的最大特点是土壤和地下水盐分组成与海水较为类似，模拟海水胁迫，探索海水对植物生长影响的机制，可为沿海滩涂盐碱地种植农作物提供理论依据 [1]。
江苏省射阳县是我国杭白菊的主产区之一，筛选耐盐能力较强的杭白菊种质对其育种和生产具有十分重要的意义。目前，已有许多学者开展菊耐盐性及小麦等海水耐盐筛选研究 [2-7]，但有关药用杭白菊的耐盐性评价研究鲜见报道 [8-9]，采用数学统计方法进行耐盐性评价的就更少。为此，本试验采用人工海水胁迫处理的方法，综合运用主成分分析、模糊数学隶属函数法、聚类分析和灰色关联等对7份杭白菊种质材料的耐盐性进行综合评价，以期建立较为完善的人工海水胁迫条件下杭白菊耐盐性鉴定方法与评价指标，为杭白菊的耐盐性评价及耐盐品种的选育提供一定的理论依据。
1 材料与方法
1 1 材料
试验于2013年3—6月在南京农业大学中药材研究所实验室进行，以NJK3、NJK2、大花4号、单瓣4号、管花2号、NJK1、红心2号共7份优良杭白菊种质为试材。
1 2 方法
1 2 1 试验设计 2013年3月下旬，选取7份种质无病虫害、健壮枝条进行扦插；4月下旬，选取根系良好、生长一致的扦插幼苗定植于放有穴盘、规格为50 cm×34 cm×16 cm 的水培箱中，每个穴盘有140个小孔，幼苗采用“之”字形间隔种植，培养液为1/2 Hoagland营养液，pH值控制在5 5～7 0之间，用空气压缩泵给水培箱营养液间歇通入空气，30 min/h； 预培养5 d，待植株恢复正常生长，以40%人工海水进行胁迫处理，以1/2 Hoagland营养液处理作为空白对照（CK），每处理重复3次。为避免盐冲击效应，初始采用浓度为10%的人工海水，后每天递增15%逐渐提高人工海水浓度，于第3天达到设定浓度后即维持不变。每隔3 d更换1次营养液。
1 2 2 相关性状及生理指标测定 试验处理前后，分别采用电子数显卡尺测定植株株高，采用分析天平称量植株的鲜质量；处理后3周，随机选取各杭白菊种质相同部位的叶片测定各个生理指标，可溶性总糖、可溶性蛋白、SOD（超氧化物歧化酶）活性、CAT（过氧化氢酶）活性、POD（过氧化物酶）活性、MDA（丙二醛）含量、游离脯氨酸（Pro）含量分别采用蒽酮比色法、Folin酚法、氮蓝四唑（NBT）法、高锰酸钾滴定法、愈创木酚法、硫代巴比妥酸（TBA）法、磺基水杨酸法测定，相对含水量、光合色素和超氧阴离子（ O-2[KG- 2]· ）采用高俊凤等的方法 [10]测定，膜透性采用美国EXTECH EC400笔式电导率TDS盐度计测定。
1 3 数据处理
采用Excel 2007和统计软件SPSS 19 0进行数据标准化处理和主成分分析，利用隶属函数法和主成分赋予权重法分别求出耐盐性的综合评价值，并进行耐盐性排序；对7份杭白菊种质材料进行耐盐性聚类分析，采用LSD法检验各处理的差异显著性。
2 结果与分析
2 1 人工海水胁迫对杭白菊叶片抗氧化酶活性、超氧阴离子、MDA含量和膜透性的影响
由表1可见，40%人工海水胁迫下，7种杭白菊种质的SOD、POD和CAT活性及超氧阴离子、MDA含量和膜透性均出现不同程度的增加，且各供试材料的差异达到显著水平；单瓣4号SOD活性增加最多，为对照的2 50倍，而NJK3的SOD活性仅为对照的0 83倍；单瓣4号CAT活性增加也最多，为对照的2 99倍，而管花2号增加最少，仅为对照的132倍；管花2号POD活性增加最多，为相应对照的2 08倍，而红心2号仅为对照的1 15倍；大花4号超氧阴离子含量最高，为1 275 0 μg/g，单瓣4号相对最小，为513 8 μg/g；NJK1的MDA含量增加最为明显，为对照的3 28倍，而单瓣4号增加幅度相对最小，为对照的1 78倍；大花4号的膜透性增加最多，为对照的2 22倍，而NJK3仅为对照的1 21倍。
2 2 人工海水胁迫对杭白菊叶片中光合色素含量的影响
由表2可见，40%人工海水胁迫下，各供试材料叶绿素a、叶绿素b、类胡萝卜素和总叶绿素含量与对照相比，总体均呈降低趋势，且各材料间的差异达到显著水平；管花2号的叶绿素a含量降低最多，较对照降低了51 81%，而NJK1降幅最小，较对照仅降低了10 84%；NJK2的叶绿素b含量降低最多，较对照降低了37 25%，而红心2号降幅最小，较对照仅降低了12 20%；管花2号的类胡萝素含量降低最多，较对照降低了61 54%，而NJK1、单瓣4号、NJK3的类胡萝卜素含量却较对照有不同程度增加；管花2号的总叶绿素含量降低最多，较对照降低了52 30%，红心2号降幅最小，较对照仅降低了13 42%；人工海水胁迫下，NJK2、单瓣4号、NJK3、红心2号的总叶绿素含量相对较高，分别为1 65、145、1 35、1 29 mg/g。 2 3 人工海水胁迫对杭白菊生长和渗透调节物质的影响
由表3可见，40%人工海水胁迫下，7种供试材料的生长量、生物量、相对含水量较对照均有不同程度的下降，可溶性总糖、可溶性蛋白和脯氨酸含量较对照均有不同程度的上升，且各材料间的差异达到显著水平；大花4号的生长量降低最多，较对照降低了72 73%，单瓣4号降低最少，较对照仅降低了38 46%；大花4号的生物量降低最多，较对照降低了7557%，红心2号降低最少，较对照仅降低了39 39%；管花2号的相对含水量降低最多，较对照降低了27 98%，单瓣4号降低最少，较对照仅降低了12 63%；单瓣4号的可溶性总糖、可溶性蛋白增加最多，分别为对照的3 93、3 33倍，而管花2号仅分别为对照的1 12、1 20倍；管花2号的脯氨酸含量增加最多，为对照的34 81倍，而NJK2仅为对照的6 51倍；人工海水胁迫处理下，单瓣4号的可溶性总糖、NJK2的可溶性蛋白及管花2号的脯氨酸含量最高，分别为0 59 mg/g、0 14 mg/g、248 9 μg/g。
2 4 耐盐性综合评价
2 4 1 主成分分析 试验测定的16个指标间部分存在显著或极显著差异，所提供的信息也发生重叠，直接利用这些指标对杭白菊幼苗进行耐盐性评价会产生较大偏差。主成分分析法充分考虑了各指标之间的重叠信息，将原来个数较多且彼此相关的指标，转换成新的个数较少且彼此独立的综合指标，可较为全面地对不同品种的抗逆性进行评价。对人工海水胁迫处理的7个杭白菊种质资源16个指标的耐盐系数进行主成分分析，结果（表4）表明，有5个主成分的累计贡献率达到92 932%，可以较充分地概括原有指标的信息，确定提取这5个主成分进行分析，把原有的16个单项指标转化为5个相互独立的综合指标。
根据各综合指标的贡献率大小，可以知道各指标的相对重要性。根据16个指标在盐胁迫下平均值的标准化值及各综合指标的相关系数（表5），求出每一种杭白菊种质资源的5个综合指标值（表6），根据主成分分析综合得分公式，可以求[CM（25]得7种供试材料所对应的主成分综合得分。由表7可见，单瓣4号、NJK1、NJK2耐盐能力较强，而大花4号和NJK3对盐胁迫较为敏感。
2 4 2 模糊数学隶属函数评价 采用16个单指标评价7种供试材料的耐盐性强弱时，评价结果存在较大差异。为更好地全面反映杭白菊不同材料真实的耐盐能力，需对各项指标进行综合考虑，而隶属函数综合评价法提供了可行性，避免了各测试指标间信息的重叠性。隶属函数综合评价以主成分分析提取的5个主成分综合评价值为指标，计算每一材料各综合指标的隶属函数值，编号分别为D1、D2、D3、D4、D5，公式为：μ（xj）=（xj-xmin）/（xmax-xmin） [10]，式中，xj代表各材料所对应的每一个综合评价值。根据综合指标贡献率的大小，计算各隶属函数值对应的权重 [11]，公式为：Wj=Pj/ P，式中，Wj表示第j个综合评价值在5个主成分的贡献率之和。由表6可见，各隶属函数值所对应的权重分别为0 262 4、0 261 4、0 177 1、0 170 9、0 128 2。从而进一步计算各材料的隶属函数综合评价值，公式为：D=∑（μj×Wj）。由表7可见，单瓣4号、NJK1、NJK2耐盐能力较强，而大花4号和管花2号对盐胁迫较为敏感。
2 4 3 聚类分析 聚类分析作为一种无管理模式的识别方法，试验结果很大程度上取决于方法所采用的距离、类间距等参数 [12]。以主成分分析提取的5个主成分综合评价值为指标，对数据进行标准化处理，以欧氏距离的平方为相似尺度，采用组间连接的方法进行系统聚类分析，将供试材料耐盐性分为强、中、弱3类。由图1可见，耐盐性强的种质为单瓣4号，耐盐性中的种质为NJK1、管花2号、NJK2、NJK3，耐盐性弱的种质为大花4号和红心2号。
2 4 4 灰色关联分析 聚类分析、主成分分析和模糊数学隶属函数法分别对杭白菊供试材料的耐盐能力进行了分门别类和综合评价排序，但对杭白菊耐盐筛选的具体量化指标缺乏一个明确的定性，不能方便有效地用于杭白菊其他材料耐盐能力的快速判断，采用更为有效的方法对耐盐指标进行明确和简化很有必要。按照灰色系统理论， 将7种供试材料的耐盐综合评价值及16个指标视为一个整体，耐盐综合评价值为母序列，其他16个量化指标为子序列。利用SPSS 19 0对原始数值进行标准化处理，使其变成无量纲数值；利用关联系数和关联度求解公式，分别求得各子序列和母序列之间的关联系数和关联度并排序，结果表明，SOD、可溶性总糖、生物量、膜透性、POD、Pro、类胡萝素、超氧阴离子、CAT、叶绿素b、MDA、生长量、叶绿素总量、可溶性蛋白、相对含水量、叶绿素a关联度系数由高到低分别为0 99、0 77、0 76、0 73、0 69、0 69、0 68、0 67、0 64、0 63、0 60、0 60、0 59、0 57、0 57和0 56，SOD、可溶性总糖、生物量、膜透性、POD、Pro等指标与杭白菊耐盐性评价有较好的关联度。
3 结论与讨论
盐胁迫是全球普遍存在的一种非生物胁迫，是盐碱地上限制植物生长和发育的主要环境因素。不同物种对盐渍环境的适应能力和适应机制均会存在明显的差异，同一物种的不同栽培类型的耐盐性也会有所不同。何先元等研究发现，大白菊和长瓣菊耐盐性较弱，红心菊和小白菊耐盐性较强；红心菊和小白菊适于在盐碱地种植 [8]。王康才等研究发现，在低浓度盐胁迫下，红心菊耐盐性高于黄菊，而在高盐胁迫下，黄菊耐盐性高于红心菊 [9]。本试验结果表明，7份杭白菊种质在人工海水胁迫下，其耐盐性有较大差别，总体表现为抗氧化酶、可溶性总糖、Pro等渗透调节物质及超氧阴离子、MDA、膜透性等指标均出现不同程度的增加，而光合色素含量、生长量、生物量及相对含水量等指标较对照存在不同程度的下降。当然也有例外，NJK1、单瓣4号、NJK3在遭受胁迫时，类胡萝素含量与对照相比没有降低，反而有不同程度的增加，这说明一定程度的胁迫作用会促进部分材料的某些生理进程，并表现出对盐胁迫具有较强的适应性。基因决定表现型，杭白菊不同材料耐盐性不同，也许还与杭白菊本身基因高度杂合、染色体数从43条到55条不等有关，这也说明了杭白菊的耐盐能力不能仅凭借几个简单的指标就可以评价。 植物耐盐性是由多因素相互作用而形成的一个十分复杂的生理生态响应过程，不同植物不同时期的耐盐性都会所不同，任何单一指标的测定与研究都不可能有效、准确地评价植物的耐盐性，而评价指标过多，不仅增加研究难度，也会造成资源的浪费。因此，找出那些起关键作用的指标，可以达到以较少指标准确评价植物耐盐性的目的，而综合利用形态生理和生化等指标，对植物的耐盐性进行综合评价是比较科学合理的。管志勇等运用主成分分析，提炼出叶片伤害率、叶绿素a含量、叶绿素总量、净光合速率、丙二醛含量 5个关键筛选指标， 对大岛野路菊等5种菊花的耐盐性进行了综合评价 [3]；李丰先等通过模糊数学隶属函数法和主成分分析法，对各高粱品种进行耐碱性综合评定，筛选出根长、发芽指数、叶干重等指标用于大量高粱品种萌发期耐碱性的筛选，并对结果进行聚类分析，对高粱品种进行了耐碱性分类 [13]；杨升等利用相关分析、主成分分析和隶属函数法，对3个树种的耐盐性进行综合评价，并通过灰色关联分析，最终筛选出株高、相对生长量、膜透性、脯氨酸、根茎叶中的Na 和K 含量等作为耐盐性评价的重要指标 [14]。本研究结合前人的研究成果，采用主成分分析、模糊数学隶属函数法、聚类分析、灰色关联分析方法，利用16个指标对7份杭白菊种质材料的耐盐性进行综合评价，结果发现，主成分分析和模糊数学隶属函数法对NJK2、大花4号、单瓣4号、NJK1、红心2号的耐盐性评价较为一致，仅对NJK3和管花2号的评价稍有不同；将系统聚类分析与主成分分析、模糊数学隶属函数法得分排序的结果对比发现，除红心2号和NJK3外，其他供试材料3种评价方法结果较为一致，尤其盐能力较强杭白菊种质单瓣4号和对盐胁迫较为敏感的种质大花4号，3种评价方法的结果完全一致。聚类分析法与主成分分析、模糊数学隶属函数法之所以存在耐盐性评价差异，这是由于聚类分析法是一种等权重的综合评价方法，而实际上不同指标在评价耐盐性上的权重是不相同的。
灰色关联分析结果显示，SOD、可溶性总糖、生物量、膜透性、POD、Pro等指标与杭白菊耐盐性评价有较大的关联度，这些指标可反映植物在盐胁迫环境下的生长情况、膜系统稳定性、渗透调节能力和抗氧化酶膜透性的大小，相应反映了植物叶片的受损程度，膜透性越大，叶片受损越严重，必然导致叶绿素含量降低。Mohamed等认为，膜透性可以作为植物耐盐能力评价的一个重要指标，而脯氨酸是植物在逆境环境条件下的一种重要渗透调节物质，在植物体内不仅作为渗透调节剂调节渗透压，还起着改善植物微环境、保证其他生化反应正常进行的作用 [15]。Wided等研究发现，盐胁迫下叶片中的脯氨酸会显著积累，这可以作为一个耐盐的评价指标 [16]；可溶性糖是很多非盐生植物的主要渗透调节剂 [17-18]，它也是合成其他有机溶质的碳架和能量来源 [19]，对细胞膜和原生质胶体有稳定作用，还可在细胞内无机离子浓度高时起保护酶类的作用；SOD、POD是膜系统的抗氧化酶，可以减轻膜系统的氧化损伤；脯氨酸、可溶性糖、MDA含量、SOD、POD活性可以作为鉴定番茄幼苗耐盐性的生理指标 [7]。本试验将灰色关联分析通过与主成分分析中因子载荷矩阵进行对比发现，可溶性总糖、POD、Pro与主成分高度相关，因子载荷矩阵虽然不能直接反映相关指标对耐盐性的影响，但却可以通过影响主成分并最终对杭白菊耐盐性产生重大影响。因此，利用SOD、可溶性总糖、生物量、膜透性、POD、Pro这6个指标对杭白菊耐盐性的快速评价较为合适。同时，综合评价认为，杭白菊种质单瓣4号、NJK1、NJK2耐盐能力较强，可以作为杭白菊耐盐育种的优良种质资源。
需说明的是，本研究是在模拟海水的盐分条件下在室内所做的水培试验，与实际滨海盐碱地条件还有一定差别，自然条件下盐碱地植物会受各种因素的影响，各种逆境条件相互作用、相互影响，植物在不同时期、不同环境下对盐胁迫的反应会不尽相同；因此，大田试验和不同生育期对比试验将是下一步的研究方向。
参考文献：
[1] 唐奇志，刘兆普，陈铭达，等 海水处理对向日葵幼苗生长及叶片一些生理特性的影响[J] 植物学通报，2004，21（6）：667-672
[2]刘爱荣，张远兵，方园园，等 盐胁迫对金盏菊生长、抗氧化能力和盐胁迫蛋白的影响[J] 草业学报，2011，20（6）：52-59
[3]管志勇，陈发棣，滕年军，等 5种菊花近缘种属植物的耐盐性比较[J] 中国农业科学，2010，43（4）：787-794
[4]管志勇，陈素梅，陈发棣，等 32个菊花近缘种属植物耐盐性筛选[J] 中国农业科学，2010，43（19）：4063-4071
[5]薛建平，张爱民，高 翔，等 安徽药菊耐盐突变体的筛选[J] 中国中药杂志，2004，29（9）：834-837
[6]刘妍妍，吴纪中，许璋阳，等 人工海水胁迫下小麦芽期和苗期的耐盐性鉴定方法[J] 植物生理学报，2014，50（2）：214-222
[7]费 伟，陈火英，曹 忠，等 盐胁迫对番茄幼苗生理特性的影响[J] 上海交通大学学报：农业科学版，2005，23（1）：5-9， 30
[8]何先元，郭巧生，罗庆云，等 药用白菊花苗期耐盐性鉴定[J] 中国中药杂志，2003，28（6）：499-503
[9]王康才，黄 莺，汤兴利，等 药用杭白菊和黄菊及其杂交F1代耐盐特性研究[J] 中国中药杂志，2011，36（17）：2321-2324
[10] 高俊凤 植物生理学实验指导[M] 北京：高等教育出版社，2006
[11]张 耿，高洪文，王 赞，等 偃麦草属植物苗期耐盐性指标筛选及综合评价[J] 草业学报，2007，16（4）：55-61
[12] 宋江峰，李大婧，刘春泉，等 甜糯玉米软罐头主要挥发性 物质主成分分析和聚类分析[J] 中国农业科学，2010，43（10）：2122-2131
[13]李丰先，周宇飞，王艺陶，等 高粱品种萌发期耐碱性筛选与综合鉴定[J] 中国农业科学，2013，46（9）：1762-1771
[14]杨 升，刘正祥，张华新，等 3个树种苗期耐盐性综合评价及指标筛选[J] 林业科学，2013，49（1）：91-98
[15]Mohamed M F M，Karima H A S Cellular basis of salinity tolerance in plants[J] Environmental and Experimental Botany，2004，52（2）：113-122
[16]Wided M，Nader B A，Ahmed D，et al Salt tolerance of the annual halophyte Cakile maritima as affected by the provenance and the developmental stage[J] Acta Physiologiae Plantarum，2007，29（4）：375-384
[17]赵可夫 植物抗盐生理[M] 北京：中国科学技术出版社，1993
[18]刘祖祺，张石城 植物抗性生理学[M] 北京：中国农业出版社，1994：369
[19]张海燕，赵可夫 盐分和水分胁迫对盐地碱蓬幼苗渗透调节效应的研究[J] 植物学报：英文版，1998，40（1）：61-66
关键词： 杭白菊；人工海水；主成分分析；综合评价；耐盐性；生长指标；生理指标
中图分类号： S682 1 10 1 文献标志码： A
文章编号：1002-1302（2015）08-0234-05
收 稿日期：2014-08-20
菊花为菊科植物菊（Chrysanthemum morifolium Ramat ）的干燥头状花序，是传统常用中药，具有散风清热、平肝明目、清热解毒的功能。近年来，随着国内外食品、医药和饮料加工业的迅猛发展，作为药食兼用的菊花用量激增。江苏省射阳县东临黄海，土地盐碱程度参差不齐，可开发沿海滩涂面积巨大。滨海盐土的最大特点是土壤和地下水盐分组成与海水较为类似，模拟海水胁迫，探索海水对植物生长影响的机制，可为沿海滩涂盐碱地种植农作物提供理论依据 [1]。
江苏省射阳县是我国杭白菊的主产区之一，筛选耐盐能力较强的杭白菊种质对其育种和生产具有十分重要的意义。目前，已有许多学者开展菊耐盐性及小麦等海水耐盐筛选研究 [2-7]，但有关药用杭白菊的耐盐性评价研究鲜见报道 [8-9]，采用数学统计方法进行耐盐性评价的就更少。为此，本试验采用人工海水胁迫处理的方法，综合运用主成分分析、模糊数学隶属函数法、聚类分析和灰色关联等对7份杭白菊种质材料的耐盐性进行综合评价，以期建立较为完善的人工海水胁迫条件下杭白菊耐盐性鉴定方法与评价指标，为杭白菊的耐盐性评价及耐盐品种的选育提供一定的理论依据。
1 材料与方法
1 1 材料
试验于2013年3—6月在南京农业大学中药材研究所实验室进行，以NJK3、NJK2、大花4号、单瓣4号、管花2号、NJK1、红心2号共7份优良杭白菊种质为试材。
1 2 方法
1 2 1 试验设计 2013年3月下旬，选取7份种质无病虫害、健壮枝条进行扦插；4月下旬，选取根系良好、生长一致的扦插幼苗定植于放有穴盘、规格为50 cm×34 cm×16 cm 的水培箱中，每个穴盘有140个小孔，幼苗采用“之”字形间隔种植，培养液为1/2 Hoagland营养液，pH值控制在5 5～7 0之间，用空气压缩泵给水培箱营养液间歇通入空气，30 min/h； 预培养5 d，待植株恢复正常生长，以40%人工海水进行胁迫处理，以1/2 Hoagland营养液处理作为空白对照（CK），每处理重复3次。为避免盐冲击效应，初始采用浓度为10%的人工海水，后每天递增15%逐渐提高人工海水浓度，于第3天达到设定浓度后即维持不变。每隔3 d更换1次营养液。
1 2 2 相关性状及生理指标测定 试验处理前后，分别采用电子数显卡尺测定植株株高，采用分析天平称量植株的鲜质量；处理后3周，随机选取各杭白菊种质相同部位的叶片测定各个生理指标，可溶性总糖、可溶性蛋白、SOD（超氧化物歧化酶）活性、CAT（过氧化氢酶）活性、POD（过氧化物酶）活性、MDA（丙二醛）含量、游离脯氨酸（Pro）含量分别采用蒽酮比色法、Folin酚法、氮蓝四唑（NBT）法、高锰酸钾滴定法、愈创木酚法、硫代巴比妥酸（TBA）法、磺基水杨酸法测定，相对含水量、光合色素和超氧阴离子（ O-2[KG- 2]· ）采用高俊凤等的方法 [10]测定，膜透性采用美国EXTECH EC400笔式电导率TDS盐度计测定。
1 3 数据处理
采用Excel 2007和统计软件SPSS 19 0进行数据标准化处理和主成分分析，利用隶属函数法和主成分赋予权重法分别求出耐盐性的综合评价值，并进行耐盐性排序；对7份杭白菊种质材料进行耐盐性聚类分析，采用LSD法检验各处理的差异显著性。
2 结果与分析
2 1 人工海水胁迫对杭白菊叶片抗氧化酶活性、超氧阴离子、MDA含量和膜透性的影响
由表1可见，40%人工海水胁迫下，7种杭白菊种质的SOD、POD和CAT活性及超氧阴离子、MDA含量和膜透性均出现不同程度的增加，且各供试材料的差异达到显著水平；单瓣4号SOD活性增加最多，为对照的2 50倍，而NJK3的SOD活性仅为对照的0 83倍；单瓣4号CAT活性增加也最多，为对照的2 99倍，而管花2号增加最少，仅为对照的132倍；管花2号POD活性增加最多，为相应对照的2 08倍，而红心2号仅为对照的1 15倍；大花4号超氧阴离子含量最高，为1 275 0 μg/g，单瓣4号相对最小，为513 8 μg/g；NJK1的MDA含量增加最为明显，为对照的3 28倍，而单瓣4号增加幅度相对最小，为对照的1 78倍；大花4号的膜透性增加最多，为对照的2 22倍，而NJK3仅为对照的1 21倍。
2 2 人工海水胁迫对杭白菊叶片中光合色素含量的影响
由表2可见，40%人工海水胁迫下，各供试材料叶绿素a、叶绿素b、类胡萝卜素和总叶绿素含量与对照相比，总体均呈降低趋势，且各材料间的差异达到显著水平；管花2号的叶绿素a含量降低最多，较对照降低了51 81%，而NJK1降幅最小，较对照仅降低了10 84%；NJK2的叶绿素b含量降低最多，较对照降低了37 25%，而红心2号降幅最小，较对照仅降低了12 20%；管花2号的类胡萝素含量降低最多，较对照降低了61 54%，而NJK1、单瓣4号、NJK3的类胡萝卜素含量却较对照有不同程度增加；管花2号的总叶绿素含量降低最多，较对照降低了52 30%，红心2号降幅最小，较对照仅降低了13 42%；人工海水胁迫下，NJK2、单瓣4号、NJK3、红心2号的总叶绿素含量相对较高，分别为1 65、145、1 35、1 29 mg/g。 2 3 人工海水胁迫对杭白菊生长和渗透调节物质的影响
由表3可见，40%人工海水胁迫下，7种供试材料的生长量、生物量、相对含水量较对照均有不同程度的下降，可溶性总糖、可溶性蛋白和脯氨酸含量较对照均有不同程度的上升，且各材料间的差异达到显著水平；大花4号的生长量降低最多，较对照降低了72 73%，单瓣4号降低最少，较对照仅降低了38 46%；大花4号的生物量降低最多，较对照降低了7557%，红心2号降低最少，较对照仅降低了39 39%；管花2号的相对含水量降低最多，较对照降低了27 98%，单瓣4号降低最少，较对照仅降低了12 63%；单瓣4号的可溶性总糖、可溶性蛋白增加最多，分别为对照的3 93、3 33倍，而管花2号仅分别为对照的1 12、1 20倍；管花2号的脯氨酸含量增加最多，为对照的34 81倍，而NJK2仅为对照的6 51倍；人工海水胁迫处理下，单瓣4号的可溶性总糖、NJK2的可溶性蛋白及管花2号的脯氨酸含量最高，分别为0 59 mg/g、0 14 mg/g、248 9 μg/g。
2 4 耐盐性综合评价
2 4 1 主成分分析 试验测定的16个指标间部分存在显著或极显著差异，所提供的信息也发生重叠，直接利用这些指标对杭白菊幼苗进行耐盐性评价会产生较大偏差。主成分分析法充分考虑了各指标之间的重叠信息，将原来个数较多且彼此相关的指标，转换成新的个数较少且彼此独立的综合指标，可较为全面地对不同品种的抗逆性进行评价。对人工海水胁迫处理的7个杭白菊种质资源16个指标的耐盐系数进行主成分分析，结果（表4）表明，有5个主成分的累计贡献率达到92 932%，可以较充分地概括原有指标的信息，确定提取这5个主成分进行分析，把原有的16个单项指标转化为5个相互独立的综合指标。
根据各综合指标的贡献率大小，可以知道各指标的相对重要性。根据16个指标在盐胁迫下平均值的标准化值及各综合指标的相关系数（表5），求出每一种杭白菊种质资源的5个综合指标值（表6），根据主成分分析综合得分公式，可以求[CM（25]得7种供试材料所对应的主成分综合得分。由表7可见，单瓣4号、NJK1、NJK2耐盐能力较强，而大花4号和NJK3对盐胁迫较为敏感。
2 4 2 模糊数学隶属函数评价 采用16个单指标评价7种供试材料的耐盐性强弱时，评价结果存在较大差异。为更好地全面反映杭白菊不同材料真实的耐盐能力，需对各项指标进行综合考虑，而隶属函数综合评价法提供了可行性，避免了各测试指标间信息的重叠性。隶属函数综合评价以主成分分析提取的5个主成分综合评价值为指标，计算每一材料各综合指标的隶属函数值，编号分别为D1、D2、D3、D4、D5，公式为：μ（xj）=（xj-xmin）/（xmax-xmin） [10]，式中，xj代表各材料所对应的每一个综合评价值。根据综合指标贡献率的大小，计算各隶属函数值对应的权重 [11]，公式为：Wj=Pj/ P，式中，Wj表示第j个综合评价值在5个主成分的贡献率之和。由表6可见，各隶属函数值所对应的权重分别为0 262 4、0 261 4、0 177 1、0 170 9、0 128 2。从而进一步计算各材料的隶属函数综合评价值，公式为：D=∑（μj×Wj）。由表7可见，单瓣4号、NJK1、NJK2耐盐能力较强，而大花4号和管花2号对盐胁迫较为敏感。
2 4 3 聚类分析 聚类分析作为一种无管理模式的识别方法，试验结果很大程度上取决于方法所采用的距离、类间距等参数 [12]。以主成分分析提取的5个主成分综合评价值为指标，对数据进行标准化处理，以欧氏距离的平方为相似尺度，采用组间连接的方法进行系统聚类分析，将供试材料耐盐性分为强、中、弱3类。由图1可见，耐盐性强的种质为单瓣4号，耐盐性中的种质为NJK1、管花2号、NJK2、NJK3，耐盐性弱的种质为大花4号和红心2号。
2 4 4 灰色关联分析 聚类分析、主成分分析和模糊数学隶属函数法分别对杭白菊供试材料的耐盐能力进行了分门别类和综合评价排序，但对杭白菊耐盐筛选的具体量化指标缺乏一个明确的定性，不能方便有效地用于杭白菊其他材料耐盐能力的快速判断，采用更为有效的方法对耐盐指标进行明确和简化很有必要。按照灰色系统理论， 将7种供试材料的耐盐综合评价值及16个指标视为一个整体，耐盐综合评价值为母序列，其他16个量化指标为子序列。利用SPSS 19 0对原始数值进行标准化处理，使其变成无量纲数值；利用关联系数和关联度求解公式，分别求得各子序列和母序列之间的关联系数和关联度并排序，结果表明，SOD、可溶性总糖、生物量、膜透性、POD、Pro、类胡萝素、超氧阴离子、CAT、叶绿素b、MDA、生长量、叶绿素总量、可溶性蛋白、相对含水量、叶绿素a关联度系数由高到低分别为0 99、0 77、0 76、0 73、0 69、0 69、0 68、0 67、0 64、0 63、0 60、0 60、0 59、0 57、0 57和0 56，SOD、可溶性总糖、生物量、膜透性、POD、Pro等指标与杭白菊耐盐性评价有较好的关联度。
3 结论与讨论
盐胁迫是全球普遍存在的一种非生物胁迫，是盐碱地上限制植物生长和发育的主要环境因素。不同物种对盐渍环境的适应能力和适应机制均会存在明显的差异，同一物种的不同栽培类型的耐盐性也会有所不同。何先元等研究发现，大白菊和长瓣菊耐盐性较弱，红心菊和小白菊耐盐性较强；红心菊和小白菊适于在盐碱地种植 [8]。王康才等研究发现，在低浓度盐胁迫下，红心菊耐盐性高于黄菊，而在高盐胁迫下，黄菊耐盐性高于红心菊 [9]。本试验结果表明，7份杭白菊种质在人工海水胁迫下，其耐盐性有较大差别，总体表现为抗氧化酶、可溶性总糖、Pro等渗透调节物质及超氧阴离子、MDA、膜透性等指标均出现不同程度的增加，而光合色素含量、生长量、生物量及相对含水量等指标较对照存在不同程度的下降。当然也有例外，NJK1、单瓣4号、NJK3在遭受胁迫时，类胡萝素含量与对照相比没有降低，反而有不同程度的增加，这说明一定程度的胁迫作用会促进部分材料的某些生理进程，并表现出对盐胁迫具有较强的适应性。基因决定表现型，杭白菊不同材料耐盐性不同，也许还与杭白菊本身基因高度杂合、染色体数从43条到55条不等有关，这也说明了杭白菊的耐盐能力不能仅凭借几个简单的指标就可以评价。 植物耐盐性是由多因素相互作用而形成的一个十分复杂的生理生态响应过程，不同植物不同时期的耐盐性都会所不同，任何单一指标的测定与研究都不可能有效、准确地评价植物的耐盐性，而评价指标过多，不仅增加研究难度，也会造成资源的浪费。因此，找出那些起关键作用的指标，可以达到以较少指标准确评价植物耐盐性的目的，而综合利用形态生理和生化等指标，对植物的耐盐性进行综合评价是比较科学合理的。管志勇等运用主成分分析，提炼出叶片伤害率、叶绿素a含量、叶绿素总量、净光合速率、丙二醛含量 5个关键筛选指标， 对大岛野路菊等5种菊花的耐盐性进行了综合评价 [3]；李丰先等通过模糊数学隶属函数法和主成分分析法，对各高粱品种进行耐碱性综合评定，筛选出根长、发芽指数、叶干重等指标用于大量高粱品种萌发期耐碱性的筛选，并对结果进行聚类分析，对高粱品种进行了耐碱性分类 [13]；杨升等利用相关分析、主成分分析和隶属函数法，对3个树种的耐盐性进行综合评价，并通过灰色关联分析，最终筛选出株高、相对生长量、膜透性、脯氨酸、根茎叶中的Na 和K 含量等作为耐盐性评价的重要指标 [14]。本研究结合前人的研究成果，采用主成分分析、模糊数学隶属函数法、聚类分析、灰色关联分析方法，利用16个指标对7份杭白菊种质材料的耐盐性进行综合评价，结果发现，主成分分析和模糊数学隶属函数法对NJK2、大花4号、单瓣4号、NJK1、红心2号的耐盐性评价较为一致，仅对NJK3和管花2号的评价稍有不同；将系统聚类分析与主成分分析、模糊数学隶属函数法得分排序的结果对比发现，除红心2号和NJK3外，其他供试材料3种评价方法结果较为一致，尤其盐能力较强杭白菊种质单瓣4号和对盐胁迫较为敏感的种质大花4号，3种评价方法的结果完全一致。聚类分析法与主成分分析、模糊数学隶属函数法之所以存在耐盐性评价差异，这是由于聚类分析法是一种等权重的综合评价方法，而实际上不同指标在评价耐盐性上的权重是不相同的。
灰色关联分析结果显示，SOD、可溶性总糖、生物量、膜透性、POD、Pro等指标与杭白菊耐盐性评价有较大的关联度，这些指标可反映植物在盐胁迫环境下的生长情况、膜系统稳定性、渗透调节能力和抗氧化酶膜透性的大小，相应反映了植物叶片的受损程度，膜透性越大，叶片受损越严重，必然导致叶绿素含量降低。Mohamed等认为，膜透性可以作为植物耐盐能力评价的一个重要指标，而脯氨酸是植物在逆境环境条件下的一种重要渗透调节物质，在植物体内不仅作为渗透调节剂调节渗透压，还起着改善植物微环境、保证其他生化反应正常进行的作用 [15]。Wided等研究发现，盐胁迫下叶片中的脯氨酸会显著积累，这可以作为一个耐盐的评价指标 [16]；可溶性糖是很多非盐生植物的主要渗透调节剂 [17-18]，它也是合成其他有机溶质的碳架和能量来源 [19]，对细胞膜和原生质胶体有稳定作用，还可在细胞内无机离子浓度高时起保护酶类的作用；SOD、POD是膜系统的抗氧化酶，可以减轻膜系统的氧化损伤；脯氨酸、可溶性糖、MDA含量、SOD、POD活性可以作为鉴定番茄幼苗耐盐性的生理指标 [7]。本试验将灰色关联分析通过与主成分分析中因子载荷矩阵进行对比发现，可溶性总糖、POD、Pro与主成分高度相关，因子载荷矩阵虽然不能直接反映相关指标对耐盐性的影响，但却可以通过影响主成分并最终对杭白菊耐盐性产生重大影响。因此，利用SOD、可溶性总糖、生物量、膜透性、POD、Pro这6个指标对杭白菊耐盐性的快速评价较为合适。同时，综合评价认为，杭白菊种质单瓣4号、NJK1、NJK2耐盐能力较强，可以作为杭白菊耐盐育种的优良种质资源。
需说明的是，本研究是在模拟海水的盐分条件下在室内所做的水培试验，与实际滨海盐碱地条件还有一定差别，自然条件下盐碱地植物会受各种因素的影响，各种逆境条件相互作用、相互影响，植物在不同时期、不同环境下对盐胁迫的反应会不尽相同；因此，大田试验和不同生育期对比试验将是下一步的研究方向。
参考文献：
[1] 唐奇志，刘兆普，陈铭达，等 海水处理对向日葵幼苗生长及叶片一些生理特性的影响[J] 植物学通报，2004，21（6）：667-672
[2]刘爱荣，张远兵，方园园，等 盐胁迫对金盏菊生长、抗氧化能力和盐胁迫蛋白的影响[J] 草业学报，2011，20（6）：52-59
[3]管志勇，陈发棣，滕年军，等 5种菊花近缘种属植物的耐盐性比较[J] 中国农业科学，2010，43（4）：787-794
[4]管志勇，陈素梅，陈发棣，等 32个菊花近缘种属植物耐盐性筛选[J] 中国农业科学，2010，43（19）：4063-4071
[5]薛建平，张爱民，高 翔，等 安徽药菊耐盐突变体的筛选[J] 中国中药杂志，2004，29（9）：834-837
[6]刘妍妍，吴纪中，许璋阳，等 人工海水胁迫下小麦芽期和苗期的耐盐性鉴定方法[J] 植物生理学报，2014，50（2）：214-222
[7]费 伟，陈火英，曹 忠，等 盐胁迫对番茄幼苗生理特性的影响[J] 上海交通大学学报：农业科学版，2005，23（1）：5-9， 30
[8]何先元，郭巧生，罗庆云，等 药用白菊花苗期耐盐性鉴定[J] 中国中药杂志，2003，28（6）：499-503
[9]王康才，黄 莺，汤兴利，等 药用杭白菊和黄菊及其杂交F1代耐盐特性研究[J] 中国中药杂志，2011，36（17）：2321-2324
[10] 高俊凤 植物生理学实验指导[M] 北京：高等教育出版社，2006
[11]张 耿，高洪文，王 赞，等 偃麦草属植物苗期耐盐性指标筛选及综合评价[J] 草业学报，2007，16（4）：55-61
[12] 宋江峰，李大婧，刘春泉，等 甜糯玉米软罐头主要挥发性 物质主成分分析和聚类分析[J] 中国农业科学，2010，43（10）：2122-2131
[13]李丰先，周宇飞，王艺陶，等 高粱品种萌发期耐碱性筛选与综合鉴定[J] 中国农业科学，2013，46（9）：1762-1771
[14]杨 升，刘正祥，张华新，等 3个树种苗期耐盐性综合评价及指标筛选[J] 林业科学，2013，49（1）：91-98
[15]Mohamed M F M，Karima H A S Cellular basis of salinity tolerance in plants[J] Environmental and Experimental Botany，2004，52（2）：113-122
[16]Wided M，Nader B A，Ahmed D，et al Salt tolerance of the annual halophyte Cakile maritima as affected by the provenance and the developmental stage[J] Acta Physiologiae Plantarum，2007，29（4）：375-384
[17]赵可夫 植物抗盐生理[M] 北京：中国科学技术出版社，1993
[18]刘祖祺，张石城 植物抗性生理学[M] 北京：中国农业出版社，1994：369
[19]张海燕，赵可夫 盐分和水分胁迫对盐地碱蓬幼苗渗透调节效应的研究[J] 植物学报：英文版，1998，40（1）：61-66