论文部分内容阅读
摘 要：夏季高温已成为制约绣球引种与园林應用的首要限制因子，为构建绣球品种耐热性评价体系并筛选耐热品种，本研究以8个绣球品种为材料，对其高温半致死温度（HLT50）及其高温胁迫下热害指数（HII）进行观测，通过测定44 ℃/33 ℃（昼/夜）高温胁迫6 d后的14项生理生化指标，采用多元统计分析方法对其耐热性进行综合评价。结果表明：通过主成分分析法将原有14项指标转化为4个彼此独立的综合指标，利用隶属函数法得到各品种的综合评价值（D值），通过聚类分析将8个绣球品种分为3个类群，耐热型包括‘花手鞠’与‘头花’;中等耐热型包括‘爱莎’‘含羞叶’‘小町’和‘纱织小姐’;不耐热型包括‘银边’与‘灵感’。在此基础上采用逐步回归法建立了以游离脯氨酸（Pro）含量、叶绿素（Chla+Chlb）含量、类胡萝卜素（Car）含量、叶绿素a/b值（Chla/Chlb）、超氧化物歧化酶（SOD）活性、相对含水量（RWC）与可溶性糖（SS）含量7个指标为变量的耐热性评价模型。此外，HLT50与HII也可作为鉴定绣球品种耐热性的可靠方法。本研究可为绣球品种在南方城乡园林中的应用推广以及耐热型绣球品种选育提供理论依据。
关键词：绣球品种;耐热性;生理生化;隶属函数法;综合评价体系
Abstract： High temperature is one of the main limiting factors for the introduction and landscape application of Hydrangea in summer. In order to construct the heat resistance evaluation system of Hydrangea cultivars and select heat resistant cultivars， eight Hydrangea cultivars were used to observe the high temperature semi-lethal temperature （HLT50） and the heat injury indexes （HII） under high temperature stress， and 14 physiological and biochemical indexes were measured after six days of high temperature stress at 44 ℃/33 ℃， after which heat resistance was evaluated by multivariate statistical analysis method. The original 14 indexes were transformed into four independent comprehensive indexes by the principal component analysis， and the comprehensive evaluation value （D value） of each variety was obtained by the subordinate function method. According to D value， the eight Hydrangea cultivars were divided into three groups. The moderate heat resistance cultivars included ‘Essar’ ‘Elbtal’ ‘Komachi’ and ‘Shakespeare’. The heat sensitive cultivars included ‘Tricolor’ and ‘Inspiration’. The heat resistance cultivars included ‘Stockings’ and ‘Corsage’. On the basis of D value， the stepwise regression method was used to set up the variable heat resistance evaluation model with the content of proline （Pro）， chlorophyll （Chla+Chlb）， carotenoid （Car）， chlorophyll a/b value （Chla/Chlb）， superoxide dismutase （SOD） activity， relative water content （RWC） and soluble sugar （SS） contents as indicators. In addition， HLT50 and HII of high temperature could also be used as a reliable method to identify the heat resistance of Hydrangea cultivars. This study would provide a theoretical basis for the application and popularization of Hydrangea cultivars in southern urban and rural gardens as well as the breeding of heat-resistant hydrangea cultivars.
Keywords： Hydrangea cultivars; heat resistance; physiology and biochemistry; subordination function method; comprehensive evaluation system 绣球（Hydrangea macrophylla）是中国传统观赏花木之一，又名八仙花、紫阳花，在我国有着悠久的栽培历史，不仅花色丰富、枝繁叶茂、观赏期长，具有极高的观赏价值及园林应用价值，而且其繁殖容易、抗性强，为建设生态型、节约型园林十分理想的绿化材料，因而具有广阔的市场前景，在国际市场上备受园艺爱好者的青睐[1]。绣球品种繁多，现已有约500余个，并不断有新品种被培育出来[2]。目前虽然有部分品种已逐渐应用在园林绿化美化中，但由于大部分绣球品种的原产地为温带气候，一些优良品种引种至夏季高温地区后，生长发育与观赏价值均受到严重影响，高温成为制约其园林推广与应用的首要限制因子，而全球气候变暖与城市热岛效应也为引种栽培工作带来更为严峻的挑战。若要在高温地区更好地引种栽培较为合适的绣球品种，耐热性是品种选择的重要指标之一。但是，目前国内外对绣球的研究多集中在资源调查[3-4]、栽培繁殖[5-6]及育种技术[7]等方面，对其抗性生理虽偶有报道，却大多集中在干旱胁迫[8-9]、铝盐胁迫[10-11]等方面，针对绣球品种的耐热性研究相对较少[12-13]。
相关研究表明，植物的耐热性作为多基因控制的复杂性状，往往会受到多方面复杂因素的影响[14]。传统上常用单一的指标对植物的耐热性进行鉴定，然而由于不同植物对高温环境的抗性机制不同，单一的指标通常无法将准确完整的信息反映出来[15-16]。评价植物耐热能力不仅需要选择适宜的评价指标，更需要依靠科学、准确的评价方法。目前常用多元统计分析方法对植物的耐热性进行综合评价，主要采用主成分分析法结合其他数据统计方法，如Logistic回归分析法、隶属函数分析法及聚类分析法等。该评价方法目前已经广泛应用于芍药属（Paeonia）[15-17]、铁线莲属（Clematis）[18]、杜鹃属（Rhododendron）[19]及蔷薇属（Rosa）[20]等多种植物的耐热性评价，这种多元建模分析技术的评价结果与田间热害观察结果一致。然而，目前有关绣球品种在高温环境下的形态表现、生理生化响应的研究鲜见报道，耐热性综合评价体系尚未建立。彭勇政等[20]在对绣球进行耐热品种筛选时发现，相对电导率、可溶性蛋白与脯氨酸可作为绣球品种耐热性鉴定的重要指标。辛雅芬等[12]研究了相对含水量、细胞膜透性、叶绿素含量、丙二醇含量与绣球耐热性的关系，研究发现不同生理指标在高温胁迫下变化显著，并且各指标之间存在一定的联系。虽然前人对绣球耐热性已经展开了一些研究，但大多集中在高温胁迫对部分生理指标的影响，对其耐热能力的鉴定也仅用单项指标来衡量，难以全面、准确地反映绣球品种的耐热性，仍缺乏更为系统和深入的评价体系研究[12-13]。
因此，本研究对人工模拟43 ℃/33 ℃（昼/夜）高温胁迫下8个绣球品种进行热害指数观测及生理生化指标的测定，以各项指标的耐热系数作为耐热性评价依据，利用主成分分析法、隶属函数法与聚类分析法对其耐热能力进行综合评价，最后利用逐步回归分析建立最优方程，筛选出与绣球品种耐热性密切关联的相关指标，旨在建立科学的绣球品种耐热性综合评价体系并简化鉴定工作，为绣球耐热品种的选育工作及在南方高温地区城乡园林绿化中的应用与推广提供理论依据。
1 材料与方法
1.1 材料
供试材料选择目前生产上具有一定规模且观赏性较高的8个绣球园艺品种，分别为：‘头花’（‘Corsage’）、‘花手鞠’（‘Stockings’）、‘灵感’（‘Inspiration’）、‘纱织小姐’（‘Shakespeare’）、‘含羞叶’（‘Elbtal’）、‘爱莎’（‘Essar’）、‘小町’（‘Komachi’）、‘银边’（‘Tricolor’）。所选用的8个品种近些年在我国南方地区苗圃中普遍栽培繁殖，尤其以云南昆明、浙江萧山、福建泉州等地栽培、应用数量相对较多。试验材料均为生长健壮、长势一致的两年生扦插苗，每个品种5株，共40株，栽培于福建农林大学（福建省福州市）森林兰苑遮阴棚内，常规水肥管理。
1.2 方法
1.2.1 高温半致死温度测定 为避免自然高温影响，于2020年9月20日（当日温度21～26 ℃）采集样叶。每个品种选取5株，采集1片叶片，采摘后自封袋保存迅速带回实验室，用去离子水洗净叶片后擦干，按照品种将所取叶片剪成0.5 cm2的碎叶，混匀后每次称取0.2 g碎叶置于装有20 mL去离子水的试管中，分别在25、35、40、45、50、55、60、65、70 ℃条件下水浴20 min，取出后静置2 h，期间自然冷却至室温，用电导仪测定电导值（Ct），之后全部放置在100 ℃水浴锅中沸水浴15 min，取出静置冷却至室温后测定电导值（Cm）。以室温下叶片的电导值（Cck）为对照。每个温度处理重复3次。根据公式求细胞伤害率：细胞伤害率=（Ct–Cck）/（Cm–Cck）×100%[21]。
采用Logistic方程进行拟合计算，Logistic方程的表达方式为：y=k/（1+ae–bt），其中y为细胞的伤害率，t为不同的处理温度，k为细胞伤害率的饱和容量，本试验消除了本底干扰k为100%。方程经线性化处理后采用直线回归法求得方程参数a、b，代入公式T=lna/b，T值即为高温半致死温度（HLT50）。
1.2.2 高温处理与采样 胁迫试验于2020年9月开始，实验前将材料统一放置在人工气候室内进行适应性栽培7 d，將气候箱的光照强度控制在240 μmol/（m2·s），光周期为13 h/11 h（昼/夜），温度控制在25 ℃/18 ℃（昼/夜），相对湿度80%，高温处理前统一浇水，保持各盆土壤湿度基本一致。于第8天模拟高温条件开始胁迫试验，将人工气候室温度提升至43 ℃/33 ℃（昼/夜），高温处理时间为8：00—17：00，光照强度和时间与湿度条件均不变。试验过程中为减轻高温引发的干旱胁迫，每隔1 d加一定量的水以补充基质水分。于高温胁迫处理第6天上午8：00—9：00进行采样，高温胁迫的第7天恢复正常温度，用于观察恢复状况。一部分样叶剪下后迅速用液氮冷冻，放置于–80 ℃超低温冰箱内冷藏，用于抗氧化酶等生理指标的测定，另一部分样叶于采集当天进行相对含水量、细胞膜透性、光合色素含量的测定。采样时选择同一叶位的当年生成熟功能叶进行测定相关指标，每个品种处理5株。 1.2.3 热害指数观测 植物热害等级评价方法参照张佳平[22]的方法做一定调整，将热害表现分为6个等级。0级：茎叶正常生长、几乎无热害表现;1级：少于1/4的茎叶萎蔫或焦尖;2级：1/4～1/2的茎叶萎蔫或焦边;3级：1/2～3/4的茎叶枯焦或穿孔;4级：3/4以上茎叶枯焦;5级：整株枯萎或死亡。根据公式计算出热害指数：
热害指数=[∑各株级数/（最高级数×总株数）]×100%。
1.2.4 生理指标测定 相对含水量（RWC）的测定采用烘干饱和称重法;细胞膜透性（cell membrane permeability，CMP）的测定采用电导率法;丙二醇（malondialdehyde，MDA）含量的测定采用硫代巴比妥酸TBA比色法;叶绿素a（chlorophyll a，Chla）含量、叶绿素b（chlorophyll b，Chlb）含量、叶绿素a/b值（Chla/Chlb）、叶绿素（chlorophyll content，Chla+Chlb）总含量、类胡萝卜素（carotenoids，Car）含量的测定采用酒精浸提法;超氧化物歧化酶（superoxide dismutase，SOD）活性的测定采用氮蓝四唑NBT法;过氧化物酶（peroxidase，POD）活性的测定采用愈创木酚法;过氧化氢酶（catalase，CAT）活性的测定采用紫外吸收法;可溶性糖（soluble sugar，SS）含量的测定采用蒽酮比色法;可溶性蛋白（soluble protein，SP）含量的测定采用考马斯亮蓝G-250试剂染色法;游离脯氨酸（proline，Pro）含量的测定采用酸性茚三酮显色法。试验方法均参考自李合生[21]的《植物生理生化实验原理和技术》，在浓度配比、反应时间及浸提、煮沸、冷却与烘干时长上适当改进，所有指标的测定均重复3次。
1.2.5 耐热性综合评价 （1）对各项指标进行耐热系数转化，公式如下：耐热系数（α）=高温处理测定值/对照测定值× 100% （1）
（2）综合评价方法。隶属函数法（subordination function method，SF method）和主成分分析法（principal component analysis，PCA）相结合的方法进行综合评价[23-24]，首先将数据输入SPSS 20.0进行相关性分析，并获得特征值、贡献率与累积贡献率，计算主成分得分（principal component score，PCS）并利用隶属函数公式进行定量转换，公式如下：
将8个品种各自的PCS代入公式（2）。Xj为某品种第j个主成分的得分，Xmin为8个品种第j个主成分得分内的最小值，Xmax为8个品种第j个主成分得分内的最大值，X（PCSj）为第j个主成分得分的隶属函数值。将主成分得分的隶属函数值代入公式（3），Wj为第j个主成分的贡献率比率，Pj为第j个主成分的贡献率。所获得的隶属函数值之和[公式（4）]即为该品种耐热性的综合得分（comprehensive score，CS），即耐热性综合评价值（D值），综合得分越高则耐热性越强。根据D值进行聚类分析和多元逐步回归分析。
1.3 数据处理
利用Microsoft Office Excel 2016软件进行数据整理与图表绘制。在SPSS 20.0软件中进行差异显著性检验、多重比较及双变量相关性分析、
主成分分析、隶属函数计算、系统聚类分析及多元逐步回归分析。
2 结果与分析
2.1 不同绣球品种的高温半致死温度
耐热性研究方法有很多，通过电导法配合Logistic回归方程的拟合计算出待测植物的高温半致死温度（HLT50）来评价植物耐热性是其中的一种[25]。通过对不同温度处理下各绣球品种叶片细胞伤害率的测定结果进行Logistic方程拟合，拟合度R2在0.804～0.968之间（表1），均大于决定系数显著性临界值，说明拟合结果符合Logistic方程的变化规律且准确度较高。
HLT50是用于评价植物耐热性的重要指标之一，与绣球耐热性呈正相关。所测8个品种的HLT50范围在55.10～59.02 ℃之间（表1），平均值为57.41 ℃。根据HLT50数值，参试8个绣球品种耐热性排序依次为：‘头花’>‘花手鞠’>‘小町’>‘爱莎’>‘含羞叶’>‘纱织小姐’>‘银边’>‘灵感’。
2.2 不同绣球品种的热害指数观测
热害指数（HII）反映了不同绣球品种的耐热能力，热害指数越低则植株耐热性相对越强，此方法简单易行，可作为评价植物耐热性强弱的指标之一[22]。高温胁迫对绣球品种茎叶生长势影响显著，各品种的热害指数均显著升高。结果表明，43 ℃高温胁迫6 d后，8个品种形态与对照组相比均发生了较为明显的变化，各品种叶片均表现出不同程度的叶缘卷曲、褐变失绿、失水皱破、萎蔫干枯、脱落等热害症状，严重影响其观赏价值，其中品种‘纱织小姐’与‘银边’全株超过3/4的叶片发生褐变、老叶枯萎甚至脱落，受热害程度最严重;品种‘花手鞠’与‘头花’受热害程度相对较轻，表现出较强的耐热性。各品种根据热害指数由到高到低次序为：‘纱织小姐’>‘银边’>‘灵感’>‘小町’>‘爱莎’>‘含羞叶’>‘花手鞠’≈‘头花’，热害指数分别为：14.4%、13.6%、12.8%、8.8%、7.2%、6.4%、5.6%與5.6%。
2.3 绣球品种的耐热性综合评价
2.3.1 各项指标的耐热系数及相关性分析 将43 ℃/33 ℃高温胁迫6 d后的各项指标测定结果进行耐热系数α值转化（表2）。由表2可知，各绣球品种经高温胁迫后的CMP、MDA含量、POD活性、CAT活性、SS含量、SP含量与Pro含量与对照相比均有提升（α>100%）;RWC与对照相比有所下降（α<100%）;Chla含量、Chlb含量、Chla+Chlb含量、Chla/Chlb值、Car含量与SOD活性在品种‘含羞叶’‘小町’‘爱莎’‘银边’‘纱织小姐’和‘灵感’中与对照相比有所下降，在品种‘头花’‘花手鞠’中与对照相比有所提升。此外，14项耐热指标在经过耐热系数转化后变化幅度差异显著，其中变化幅度最大的指标为POD活性，变化幅度最小的为RWC，并且不同品种在同一指标内的变化差异较大。 从表3可见，RWC与SOD活性之间存在极显著相关关系，与Pro含量之间存在显著相关关系;Chla含量与Chlb含量、Chla/Chlb及Chla+Chlb含量之间存在极显著相关关系，与SOD活性、SS含量及Pro含量之间存在显著相关关系;Chlb含量与Chla+Chlb含量之间存在极显著相关关系;Chla/Chlb与Chla+Chlb含量、SS含量及Pro含量之间存在显著相关关系;Chla+Chlb含量与SOD活性、SS含量及Pro含量之间存在显著相关关系;Car含量与SP含量及Pro含量之间存在极显著相关关系，与SOD活性的之间存在显著相关关系;SOD活性与Pro含量之间存在极显著相关关系;POD活性与SS含量之间存在显著相关关系;SS含量与Pro含量之间存在显著相关关系;SP含量与Pro含量之间存在显著相关关系。其余各项指标之间关系不显著（表3）。可见，各单项指标之间普遍存在信息重叠现象，可能会在一定程度上影响绣球品种耐热性鉴定结果。
2.3.2 不同绣球品种各单项指标的主成分分析 主成分分析法能通过降维的方式将多个相互关联的指标转化为少数几个包含原变量大部分信息且相互独立的综合指标。特征值能在一定程度上反映主成分的影响力，当特征值大于1则可作为较好的评价依据。通过对高温胁迫下绣球品种各项指标的耐热系数进行主成分分析，前4个主成分的特征值大于1，把原本相互关联的13个单项指标重新降维成4个新的互不干扰的综合指标，依次用Z1～Z4表示，4个综合指标的特征值分别为7.299、1.987、1.589和1.195，贡献率分別为52.134%、14.189%、11.351%和8.533%，累计贡献率达86.207%（当累积贡献率超过80%则具有较强数据代表性），说明所提取的这4个综合指标基本上能代表不同绣球品种耐热性原始数据的绝大部分信息（表4）。
主成分中各项指标的荷载值大小能反映出该指标在主成分中的关键程度。对各综合指标进行特征向量分析，其中Z1中Chla含量、Chla+Chlb含量和Pro含量有较大的正载荷量;Z2中SP含量、Car含量和CMP有较大的正载荷量，Chlb含量和CAT活性有较大负荷载量;Z3中MDA含量、CAT活性与SP含量有较大的正载荷量;Z4中POD活性、Chla含量与Chlb含量有较大的正载荷量，RWC与SS含量有较大负荷载量。综合指标的贡献率的大小能反映不同综合指标的相对重要性，
主成分Z1的贡献率最大（表4），说明主成分Z1所涵盖的信息能更全面地反映8个绣球品种的耐热能力。
2.3.3 不同绣球品种耐热性综合分析 将8个绣球品种的各项生理指标的耐热系数与其主成分的荷载值相结合，根据线性组合方程计算结果得出4个主成分的得分值（comprehensive index，CI），得分值越高说明该品种在这一综合指标中的耐热表现越好。由表5可知，品种‘头花’与‘花手鞠’在第1主成分上的耐热性最强;品种‘含羞叶’与‘小町’在第2主成分上的耐热性最强;品种‘含羞叶’与‘纱织小姐’在第3主成分上。
2.3.4 不同绣球品种耐热性聚类分析 聚类分析法是一种能根据对象性质划分为不同集合的样本归类方法[26]。由图1可知，本研究根据各个材料的耐热指标综合评价值（D值），采用平方欧氏距离，用离差平方和法对不同绣球品种的耐热性进行了系统聚类分析，聚类结果在欧式距离5处可将绣球耐热性分为3类：中等耐热型品种、不耐热型品种、耐热型品种。其中，中等耐热型品种包括‘爱莎’‘含羞叶’‘小町’和‘纱织小姐’;不耐热型包括品种‘银边’与‘灵感’;耐热型包括品种‘花手鞠’与‘头花’。8个绣球品种的耐热性鉴定结果与热害指数鉴定结果、高温半至死温度鉴定结果大体一致。
2.3.5 耐热鉴定指标的筛选 为探究各项指标与耐热性之间的关系，在保证评价结果准确性的同时精简评价指标，建立可用于绣球品种耐热性评价的数学模型，将8个绣球品种的综合评价D值作为因变量，各单项指标的耐热系数作为自变量进行逐步回归分析，得出预测绣球品种耐热能力
的最优回归方程：D=–0.707–0.000068X1+ 0.001408X2+0.004131X3+0.003451X4+0.003746X5–0.000398X6–0.000038X7（R2=1，P=0.0001）。由方程可知，原有的14个指标被精简为X1～X7 7个指标，它们分别代表Pro含量、Chla+Chlb含量、Car含量、Chla/Chlb、SOD活性、RWC与SS含量的耐热系数α，可用于鉴定不同绣球品种的耐热性，即在相同条件下测定其他品种的这7个指标，并计算其耐热系数，进而利用该方程预测相应品种的耐热性。
用该回归方程对8个品种的耐热性进行预测，即把8个品种的上述7个指标的耐热系数（α值）代入该方程，得到各品种的预测值：‘含羞叶’为0.435，‘头花’为0.828，‘小町’为0.308、‘爱莎’为0.437，‘银边’为0.126，‘花手鞠’为0.665，‘纱织小姐’为0.293，‘灵感’为0.088。将各品种的预测值与综合评价值进行相关性分析，二者相关系数R2=1，在0.1水平上显著相关。说明此方程能够很好地量化各品种的综合耐热能力与各指标值间的因果关系，对绣球品种间耐热性的预测效果较好。
此外，将各品种的耐热性综合评价D值与高温半致死温度HLT50进行相关性分析，发现二者成极显著正相关（R2=0.868）;将各品种的耐热性综合评价D值与热害指数HII进行相关分析，发现二者呈显著负相关（R2=–0.743），即测定HLT50与HII也可作为鉴定绣球品种耐热性的可靠方法。
3 讨论
相关研究表明，植物细胞膜的热稳定性与其耐热能力密切相关[27]，通过电导率法结合Logistic方程拟合估算HLT50，能较为准确地反映植物耐性差异。通过该方法评价植物耐热性既便捷又高效[28-29]，现已应用在杜鹃花科（Ericaceae）[30]、蔷薇科（Rosaceae）[31]及景天科（Crassulaceae）[32]等多种植物材料的耐热性鉴定中，并已有学者验证该方法在多个经济作物与观赏花卉品种耐热性差异评价中的可行性[33-34]。但也有学者认为由于该实验方法在离体条件下进行，因而并不能准确反映出品种间的耐热差异[33]。本研究对8个绣球品种的HLT50进行测定，温度范围在55.10～59.02 ℃之间，测定结果与耐热性综合评价结果呈极显著正相关（R2=0.868），验证了其作为绣球品种耐热评价指标的可靠性。刘婉迪等[30]测定出9个杜鹃品种（Rhododendron simsii）的HLT50范围在46.92～53.84 ℃之间，刘志高等[18]对20个铁线莲品种的HLT50进行测定，HLT50最高的品种仅为42.17 ℃，其耐热能力均不如本研究参试的8个绣球品种。从耐热性角度考虑，8个绣球品种的半致死温度均超过了55 ℃，说明这8个品种都有一定的高温耐受能力，是华南地区园林绿化材料更新潜在的植物资源，但其耐热能力与凌霄（Campsis grandiflora）、地锦（Parthenocissus inserta）、牛轭草（Murdannia loriformis）和细竹篙草（Murdannia simplex）等平均HLT50在60 ℃以上的观赏植物相比仍存在差距[35-36]。这一方面体现了不同植物种类间的耐热性差异，另一方面反应出绣球品种对高温的耐受能力有限，在实际应用中需要尽量配置在通风良好的遮荫环境下。此外，各品种的耐热性综合评价D值与HLT50呈极显著正相关（R2=0.868），说明绣球的高温半致死温度能较好地反映其耐热性，可作为其耐热性评价的关键指标。 持续高温会造成绣球的叶片失绿褐变、萎蔫下垂甚至枯萎，这是由于高温引起植物体内的酶代谢失活最终导致叶绿素合成减少，或是叶片细胞含水量降低的表现[16，37]。HII通过将高温胁迫下植物的表型特征数量化直观地反映受热害程度，被作为可靠评价指标广泛地应用在耐热性评价研究中[17-18]。本研究结果表明，高温胁迫下各绣球品种皆表现出不同程度的热害症状，其中‘花手鞠’与‘头花’在高温环境下形态表现较好，为耐热品种;‘纱织小姐’与‘银边’受热害程度较高，为热敏品种，该鉴定结果与耐热综合评价D值聚类鉴定结果基本一致。此外，不同耐热型绣球品种的HII差异显著，各品种的耐热性综合评价D值与热害指数HII呈显著负相关（R2=–0.743），进一步验证了以HII作为绣球耐热性评价关键指标的可靠性，为通过表型观测对品种耐热性进行预判提供了依据。
植物能利用自身体内的一系列生理代谢途径来抵御或适应高温胁迫，比如通过产生兼容溶质组织蛋白质与细胞结构，通过渗透调节保持细胞充盈，通过抗氧化系统重建细胞氧化还原平衡与内稳态等[38]。不同植物种类对高温胁迫的动态响应机制存在较大差异，植物的耐热能力不仅与高温胁迫的程度与胁迫时间有关，也受到植物的遗传背景、表型性状、生长发育阶段与生理代谢水平等多方面因素的影响，是在植物与各环境因子之间相互作用的复杂过程中形成的[39]。本研究采用人工气候室模拟高温胁迫环境，有效地避免了干旱、强光、高湿等复合逆境的影响，便于更加客观地对绣球品种的耐热性进行评价。可用于观赏植物耐热性鉴定的指标，不同指标对高温胁迫的敏感程度表现不一，且各指标间存在着复杂的关联性[40-41]。本研究选用的14项生理生化指标在高温胁迫下的响应方式各异，其相关性分析结果也表明各单项指标所提供的耐热响应信息发生不同程度的交叉重叠，这与上述观点一致，说明根据其中任何一项指标对绣球品种的耐热性进行评价都将得到不同的评价结果，因此无法直接通过相互矛盾的单一指标对绣球各品种的耐热性进行评价，耐热性作为复杂的综合性状需要依靠多元的统计分析方法，建立高效、科学的绣球耐热性综合评价体系[42]。
为此，本研究采用主成分分析法结合隶属函数分析法进一步将各指标进行综合分析，以期弥补单项指标评价耐热性的缺陷，得到更为准确的耐热性评价结果。通过测定对照及高温胁迫处理下不同绣球品种的14个单项指标，以各项指标的耐热系数α为基础，利用主成分分析法将原有的14个单项指标降维成4个互不干扰的综合指标，再利用隶属函数分析法获得不同绣球品种的耐热性综合评价D值。由于D值不仅与各指标间的关联性相关，也与各指标的权重相关，所以利用其对绣球耐热性评价更具有系统性和科学性。根据综合评价D值进行聚类分析，可将各参试绣球品种划分为耐热型、中等耐热型和不耐热型3个类群，耐热型品种包括‘花手鞠’与‘头花’，其耐热表现较好，是应用在南方夏季高温地区的潜在优良品种;中等耐热型品种包括‘爱莎’‘含羞叶’‘小町’与‘纱织小姐’;不耐热型品种包括‘银边’与‘灵感’，其耐热表现较差，更适合栽植在温度不高的温带地区。
此外，本研究利用逐步回归法建立了能够快速、准确鉴定绣球品种耐热性评价模型，筛选出Pro含量、Chla+Chlb含量、Car含量、Chla/Chlb、SOD活性、RWC与SS含量等7个显著影响绣球品种耐热性单项指标。经过预测值验证，结果表明由这7项指标构建的评价模型能能夠很好地量化各品种的综合耐热能力与各指标值间的因果关系，进而高效、准确地预测不同绣球品种耐热能力的强弱，可作为绣球品种耐热性鉴定的关键指标，进一步为绣球品种耐热性鉴定工作简化了工作量。
在本研究的基础上，更多研究可以进一步展开，如通过解剖技术进一步探究绣球叶片显微结构与其耐热性的联系;选择耐热性强的品种作为亲本培育新种质展开育种工作;随着分子生物技术的成熟，对植物的耐热性研究也已达到分子水平，选择耐热与不耐热绣球品种构建基因文库，筛选与耐热性直接相关的关键基因，通过基因与表达产物进一步揭示绣球对高温的响应机制，导入外源基因进行耐热性育种是未来研究的方向之一。此外，绣球由于其极高的观赏价值与经济价值，在园林美化、盆花销售与切花生产中扮演越来越重要的角色，但其生产栽培仍受不同地区气候差异的限制，故除研究品种本身的抗逆性强弱和响应机制外，探索更多有助于提升其抗性的方法，从而降低栽培管理成本并扩大生产应用区域也是较好的切入点。
参考文献
乔 谦， 王江勇， 陶吉寒. 绣球属植物研究进展[J]. 农学学报， 2020， 10（4）： 60-64.
李艳香. 19个绣球属种质资源的SRAP分析[D]. 长沙： 湖南农业大学， 2009.
彭尽晖， 周朴华， 周红灿， 等. 湖南省绣球属植物资源调查[J]. 湖南农业大学学报（自然科学版）， 2008（5）： 563-567.
孙 强， 虞秀明， 姚红军. 绣球属品种资源收集及扦插生根能力比较[J]. 北方园艺， 2016（11）： 71-73.
张国华， 陈旦旦， 曹受金， 等. 绣球花扦插繁殖技术研究[J]. 现代园艺， 2019， 42（23）： 24-25.
王红梅， 刘春风， 孟紫莲， 等. 圆锥绣球‘北极熊’组培快繁体系的建立[J]. 分子植物育种， 2021， 19（7）： 2305-2311.
Uemachi T， Nishio T. Relationship between the inflorescence structure and setting of decorative flowers on lacecap cultivars of Hydrangea macrophylla （Thunb.） Ser. and H. serrata （Thunb.） Ser[J]. Horticultural Research， 2005， 4（4）： 435-438. 章 毅， 韦孟琪， 孙欧文， 等. 不同绣球品种对干旱胁迫的生理响应及抗旱机制研究[J]. 西北林学院学报， 2018， 33（1）： 90-97.
蔡建国， 章 毅， 孙欧文， 等. 绣球抗旱性综合评价及指标体系构建[J]. 应用生态学报， 2018， 29（10）： 3175-3182.
李叶华， 陈 爽， 赵 冰. 铝胁迫对绣球组培苗生长及生理特性的影响[J]. 浙江农林大学学报， 2020， 37（6）： 1064-1070.
汪 霖. NaCl胁迫对绣球生理指标的影响[J]. 北方园艺， 2011（24）： 80-83.
辛雅芬， 石玉波. 高温对八仙花生理指标的影响[J]. 安徽农业科学， 2009， 37（11）： 5204-5206.
蒋昌华， 李 健. 高温胁迫对八仙花品种生理特性的影响[C]//中国植物园（第十五期）， 西安： 中国林业出版社， 2011： 196-198.
Gulen H， Eris A. Effect of heat stress on peroxidase activity and total protein content in strawberry plants[J]. Plant Science， 2003， 166（3）： 739-744.
Tian Z G， Yang Y， Wang F. A comprehensive evaluation of heat tolerance in nine cultivars of marigold[J]. Horticulture Environment and Biotechnology， 2015， 56（6）： 749-755.
Stirbet A， Govindjee. On the relation between the Kautsky effect （chlorophyll a fluorescence induction） and photosystem II： Basics and applications of the OJIP fluorescence transient[J]. Journal of Photochemistry and Photobiology，2011， 104（1-2）： 236-257.
張佳平， 李丹青， 聂晶晶， 等. 高温胁迫下芍药的生理生化响应和耐热性评价[J]. 核农学报， 2016， 30（9）： 1848-1856.
刘志高， 邵伟丽， 申亚梅， 等. 铁线莲品种耐热性分析及评价指标筛选[J]. 核农学报， 2020， 34（1）： 203-213.
李 辉. 杜鹃耐热资源评价及其生理机制研究[D]. 扬州： 扬州大学， 2019.
彭勇政， 刘智媛， 朱晓非， 等. 5个月季品种高温处理后生理指标变化及其耐热性评价[J]. 上海交通大学学报（农业科学版）， 2019， 37（5）： 53-58.
李合生. 植物生理生化实验原理和技术[M]. 北京： 高等教育出版社， 2004.
张佳平. 芍药在杭州栽培的耐热评价及地下芽休眠机理研究[D]. 杭州： 浙江大学， 2015.
张旻桓. 湖南牡丹资源遗传多样性及耐热性研究[D]. 长沙： 中南林业科技大学， 2019.
田治国. 万寿菊属植物耐热性与抗旱性的评价及生长生理特性的研究[D]. 杨凌： 西北农林科技大学， 2012.
赵 伟. 三种国产唇柱苣苔属野生花卉的组培、扦插和幼苗半致死温度研究[D]. 上海： 上海师范大学， 2015.
尹立为. 函数数据降维方法的分析与研究[D]. 长春： 吉林大学， 2014.
Ashraf M， Saeed M M， Qureshi M J. Tolerance to high temperature in cotton （Gossypium hirsutum L.） at initial growth stages[J]. Environmental and Experimental Botany， 1994， 34（3）： 275-283.
史彦江， 宋锋惠， 俞 涛， 等. 高温胁迫下平欧杂种榛品种的耐热性研究[J]. 中国农学通报， 2015， 31（28）： 23-28.
Castor T P， Tyler T A. Determination of taxol in taxus media needles in the presence of interfering components[J]. J Liq Chromatogr Relat Technol， 1993， 16（3）： 723-731.
刘婉迪， 王 威， 谢 倩， 等. 9个杜鹃品种的高温半致死温度与耐热性评价[J]. 西北林学院学报， 2018， 33（5）： 105-110， 136.
吕晓倩. 梅花杂交育种与后代耐热性评价研究[D]. 北京： 北京林业大学， 2014.
张燕利， 高捍东， 吴锦华. 4种景天科植物耐热性测定[J]. 西南林业大学学报， 2010， 30（6）： 52-54.
Yeh D M， Lin H F. Thermostability of cell membranes as a measure of heat tolerance and relationship to flowering delay in Chrysanthemum[J]. Journal of the American Society for Horticultural Science American Society for Horticultural Science， 2003， 128（5）： 656-660. 薛盼盼， 张 强， 魏学智. 模拟干旱胁迫对酸枣种子萌发及活力的影响[J]. 西北农业学报， 2016， 25（12）： 1837-1843.
徐玉芬， 杨 静， 邓 磊， 等. 应用Logistic方程测定4种鸭跖草科屋顶绿化植物的耐热性[J]. 仲恺农业工程学院学报， 2018， 31（4）： 52-55， 70.
袁雪涛， 夏 冬， 许永利， 等. 应用Logistic方程测定四种藤本植物的耐热性[J]. 北方园艺， 2017（14）： 99-102.
Wang Q， Cheng T R， Yu X N， et al. Physiological and biochemical responses of six herbaceous peony cultivars to cold stress[J]. South African Journal of Botany， 2014， 94： 140-148.
Valliyodan B， Nguyen H T. Understanding regulatory networks and engineering for enhanced drought tolerance in plants[J]. Current Opinion in Plant Biology， 2006， 9（2）： 189-195.
Khanna-Chopra R， Chauhan S. Wheat cultivars differing in heat tolerance show a differential response to oxidative stress during monocarpic senescence under high temperature stress[J]. Protoplasma， 2015， 252（5）： 1241-1251.
Zhang Y T， Du H M. Differential accumulation of proteins in leaves and roots associated with heat tolerance in two kentucky bluegrass genotypes differing in heat tolerance[J]. Acta Physiologiae Plantarum， 2016， 38（9）： 1-15.
Cheng Z Y， Sun L， Wang X J， et al. Ferulic acid pretreatment alleviates heat stress in blueberry seedlings by inducing antioxidant enzymes proline and soluble sugars[J]. Biologia Plantarum， 2018， 62（3）： 534-542.
Patton A J， Cunningham S M， Volenec J J， et al. Differences in freeze tolerance of zoysiagrasses： II. carbohydrate and proline accumulation[J]. Crop Science， 2007， 47（5）： 2170-2181.
責任编辑：沈德发
关键词：绣球品种;耐热性;生理生化;隶属函数法;综合评价体系
Abstract： High temperature is one of the main limiting factors for the introduction and landscape application of Hydrangea in summer. In order to construct the heat resistance evaluation system of Hydrangea cultivars and select heat resistant cultivars， eight Hydrangea cultivars were used to observe the high temperature semi-lethal temperature （HLT50） and the heat injury indexes （HII） under high temperature stress， and 14 physiological and biochemical indexes were measured after six days of high temperature stress at 44 ℃/33 ℃， after which heat resistance was evaluated by multivariate statistical analysis method. The original 14 indexes were transformed into four independent comprehensive indexes by the principal component analysis， and the comprehensive evaluation value （D value） of each variety was obtained by the subordinate function method. According to D value， the eight Hydrangea cultivars were divided into three groups. The moderate heat resistance cultivars included ‘Essar’ ‘Elbtal’ ‘Komachi’ and ‘Shakespeare’. The heat sensitive cultivars included ‘Tricolor’ and ‘Inspiration’. The heat resistance cultivars included ‘Stockings’ and ‘Corsage’. On the basis of D value， the stepwise regression method was used to set up the variable heat resistance evaluation model with the content of proline （Pro）， chlorophyll （Chla+Chlb）， carotenoid （Car）， chlorophyll a/b value （Chla/Chlb）， superoxide dismutase （SOD） activity， relative water content （RWC） and soluble sugar （SS） contents as indicators. In addition， HLT50 and HII of high temperature could also be used as a reliable method to identify the heat resistance of Hydrangea cultivars. This study would provide a theoretical basis for the application and popularization of Hydrangea cultivars in southern urban and rural gardens as well as the breeding of heat-resistant hydrangea cultivars.
Keywords： Hydrangea cultivars; heat resistance; physiology and biochemistry; subordination function method; comprehensive evaluation system 绣球（Hydrangea macrophylla）是中国传统观赏花木之一，又名八仙花、紫阳花，在我国有着悠久的栽培历史，不仅花色丰富、枝繁叶茂、观赏期长，具有极高的观赏价值及园林应用价值，而且其繁殖容易、抗性强，为建设生态型、节约型园林十分理想的绿化材料，因而具有广阔的市场前景，在国际市场上备受园艺爱好者的青睐[1]。绣球品种繁多，现已有约500余个，并不断有新品种被培育出来[2]。目前虽然有部分品种已逐渐应用在园林绿化美化中，但由于大部分绣球品种的原产地为温带气候，一些优良品种引种至夏季高温地区后，生长发育与观赏价值均受到严重影响，高温成为制约其园林推广与应用的首要限制因子，而全球气候变暖与城市热岛效应也为引种栽培工作带来更为严峻的挑战。若要在高温地区更好地引种栽培较为合适的绣球品种，耐热性是品种选择的重要指标之一。但是，目前国内外对绣球的研究多集中在资源调查[3-4]、栽培繁殖[5-6]及育种技术[7]等方面，对其抗性生理虽偶有报道，却大多集中在干旱胁迫[8-9]、铝盐胁迫[10-11]等方面，针对绣球品种的耐热性研究相对较少[12-13]。
相关研究表明，植物的耐热性作为多基因控制的复杂性状，往往会受到多方面复杂因素的影响[14]。传统上常用单一的指标对植物的耐热性进行鉴定，然而由于不同植物对高温环境的抗性机制不同，单一的指标通常无法将准确完整的信息反映出来[15-16]。评价植物耐热能力不仅需要选择适宜的评价指标，更需要依靠科学、准确的评价方法。目前常用多元统计分析方法对植物的耐热性进行综合评价，主要采用主成分分析法结合其他数据统计方法，如Logistic回归分析法、隶属函数分析法及聚类分析法等。该评价方法目前已经广泛应用于芍药属（Paeonia）[15-17]、铁线莲属（Clematis）[18]、杜鹃属（Rhododendron）[19]及蔷薇属（Rosa）[20]等多种植物的耐热性评价，这种多元建模分析技术的评价结果与田间热害观察结果一致。然而，目前有关绣球品种在高温环境下的形态表现、生理生化响应的研究鲜见报道，耐热性综合评价体系尚未建立。彭勇政等[20]在对绣球进行耐热品种筛选时发现，相对电导率、可溶性蛋白与脯氨酸可作为绣球品种耐热性鉴定的重要指标。辛雅芬等[12]研究了相对含水量、细胞膜透性、叶绿素含量、丙二醇含量与绣球耐热性的关系，研究发现不同生理指标在高温胁迫下变化显著，并且各指标之间存在一定的联系。虽然前人对绣球耐热性已经展开了一些研究，但大多集中在高温胁迫对部分生理指标的影响，对其耐热能力的鉴定也仅用单项指标来衡量，难以全面、准确地反映绣球品种的耐热性，仍缺乏更为系统和深入的评价体系研究[12-13]。
因此，本研究对人工模拟43 ℃/33 ℃（昼/夜）高温胁迫下8个绣球品种进行热害指数观测及生理生化指标的测定，以各项指标的耐热系数作为耐热性评价依据，利用主成分分析法、隶属函数法与聚类分析法对其耐热能力进行综合评价，最后利用逐步回归分析建立最优方程，筛选出与绣球品种耐热性密切关联的相关指标，旨在建立科学的绣球品种耐热性综合评价体系并简化鉴定工作，为绣球耐热品种的选育工作及在南方高温地区城乡园林绿化中的应用与推广提供理论依据。
1 材料与方法
1.1 材料
供试材料选择目前生产上具有一定规模且观赏性较高的8个绣球园艺品种，分别为：‘头花’（‘Corsage’）、‘花手鞠’（‘Stockings’）、‘灵感’（‘Inspiration’）、‘纱织小姐’（‘Shakespeare’）、‘含羞叶’（‘Elbtal’）、‘爱莎’（‘Essar’）、‘小町’（‘Komachi’）、‘银边’（‘Tricolor’）。所选用的8个品种近些年在我国南方地区苗圃中普遍栽培繁殖，尤其以云南昆明、浙江萧山、福建泉州等地栽培、应用数量相对较多。试验材料均为生长健壮、长势一致的两年生扦插苗，每个品种5株，共40株，栽培于福建农林大学（福建省福州市）森林兰苑遮阴棚内，常规水肥管理。
1.2 方法
1.2.1 高温半致死温度测定 为避免自然高温影响，于2020年9月20日（当日温度21～26 ℃）采集样叶。每个品种选取5株，采集1片叶片，采摘后自封袋保存迅速带回实验室，用去离子水洗净叶片后擦干，按照品种将所取叶片剪成0.5 cm2的碎叶，混匀后每次称取0.2 g碎叶置于装有20 mL去离子水的试管中，分别在25、35、40、45、50、55、60、65、70 ℃条件下水浴20 min，取出后静置2 h，期间自然冷却至室温，用电导仪测定电导值（Ct），之后全部放置在100 ℃水浴锅中沸水浴15 min，取出静置冷却至室温后测定电导值（Cm）。以室温下叶片的电导值（Cck）为对照。每个温度处理重复3次。根据公式求细胞伤害率：细胞伤害率=（Ct–Cck）/（Cm–Cck）×100%[21]。
采用Logistic方程进行拟合计算，Logistic方程的表达方式为：y=k/（1+ae–bt），其中y为细胞的伤害率，t为不同的处理温度，k为细胞伤害率的饱和容量，本试验消除了本底干扰k为100%。方程经线性化处理后采用直线回归法求得方程参数a、b，代入公式T=lna/b，T值即为高温半致死温度（HLT50）。
1.2.2 高温处理与采样 胁迫试验于2020年9月开始，实验前将材料统一放置在人工气候室内进行适应性栽培7 d，將气候箱的光照强度控制在240 μmol/（m2·s），光周期为13 h/11 h（昼/夜），温度控制在25 ℃/18 ℃（昼/夜），相对湿度80%，高温处理前统一浇水，保持各盆土壤湿度基本一致。于第8天模拟高温条件开始胁迫试验，将人工气候室温度提升至43 ℃/33 ℃（昼/夜），高温处理时间为8：00—17：00，光照强度和时间与湿度条件均不变。试验过程中为减轻高温引发的干旱胁迫，每隔1 d加一定量的水以补充基质水分。于高温胁迫处理第6天上午8：00—9：00进行采样，高温胁迫的第7天恢复正常温度，用于观察恢复状况。一部分样叶剪下后迅速用液氮冷冻，放置于–80 ℃超低温冰箱内冷藏，用于抗氧化酶等生理指标的测定，另一部分样叶于采集当天进行相对含水量、细胞膜透性、光合色素含量的测定。采样时选择同一叶位的当年生成熟功能叶进行测定相关指标，每个品种处理5株。 1.2.3 热害指数观测 植物热害等级评价方法参照张佳平[22]的方法做一定调整，将热害表现分为6个等级。0级：茎叶正常生长、几乎无热害表现;1级：少于1/4的茎叶萎蔫或焦尖;2级：1/4～1/2的茎叶萎蔫或焦边;3级：1/2～3/4的茎叶枯焦或穿孔;4级：3/4以上茎叶枯焦;5级：整株枯萎或死亡。根据公式计算出热害指数：
热害指数=[∑各株级数/（最高级数×总株数）]×100%。
1.2.4 生理指标测定 相对含水量（RWC）的测定采用烘干饱和称重法;细胞膜透性（cell membrane permeability，CMP）的测定采用电导率法;丙二醇（malondialdehyde，MDA）含量的测定采用硫代巴比妥酸TBA比色法;叶绿素a（chlorophyll a，Chla）含量、叶绿素b（chlorophyll b，Chlb）含量、叶绿素a/b值（Chla/Chlb）、叶绿素（chlorophyll content，Chla+Chlb）总含量、类胡萝卜素（carotenoids，Car）含量的测定采用酒精浸提法;超氧化物歧化酶（superoxide dismutase，SOD）活性的测定采用氮蓝四唑NBT法;过氧化物酶（peroxidase，POD）活性的测定采用愈创木酚法;过氧化氢酶（catalase，CAT）活性的测定采用紫外吸收法;可溶性糖（soluble sugar，SS）含量的测定采用蒽酮比色法;可溶性蛋白（soluble protein，SP）含量的测定采用考马斯亮蓝G-250试剂染色法;游离脯氨酸（proline，Pro）含量的测定采用酸性茚三酮显色法。试验方法均参考自李合生[21]的《植物生理生化实验原理和技术》，在浓度配比、反应时间及浸提、煮沸、冷却与烘干时长上适当改进，所有指标的测定均重复3次。
1.2.5 耐热性综合评价 （1）对各项指标进行耐热系数转化，公式如下：耐热系数（α）=高温处理测定值/对照测定值× 100% （1）
（2）综合评价方法。隶属函数法（subordination function method，SF method）和主成分分析法（principal component analysis，PCA）相结合的方法进行综合评价[23-24]，首先将数据输入SPSS 20.0进行相关性分析，并获得特征值、贡献率与累积贡献率，计算主成分得分（principal component score，PCS）并利用隶属函数公式进行定量转换，公式如下：
将8个品种各自的PCS代入公式（2）。Xj为某品种第j个主成分的得分，Xmin为8个品种第j个主成分得分内的最小值，Xmax为8个品种第j个主成分得分内的最大值，X（PCSj）为第j个主成分得分的隶属函数值。将主成分得分的隶属函数值代入公式（3），Wj为第j个主成分的贡献率比率，Pj为第j个主成分的贡献率。所获得的隶属函数值之和[公式（4）]即为该品种耐热性的综合得分（comprehensive score，CS），即耐热性综合评价值（D值），综合得分越高则耐热性越强。根据D值进行聚类分析和多元逐步回归分析。
1.3 数据处理
利用Microsoft Office Excel 2016软件进行数据整理与图表绘制。在SPSS 20.0软件中进行差异显著性检验、多重比较及双变量相关性分析、
主成分分析、隶属函数计算、系统聚类分析及多元逐步回归分析。
2 结果与分析
2.1 不同绣球品种的高温半致死温度
耐热性研究方法有很多，通过电导法配合Logistic回归方程的拟合计算出待测植物的高温半致死温度（HLT50）来评价植物耐热性是其中的一种[25]。通过对不同温度处理下各绣球品种叶片细胞伤害率的测定结果进行Logistic方程拟合，拟合度R2在0.804～0.968之间（表1），均大于决定系数显著性临界值，说明拟合结果符合Logistic方程的变化规律且准确度较高。
HLT50是用于评价植物耐热性的重要指标之一，与绣球耐热性呈正相关。所测8个品种的HLT50范围在55.10～59.02 ℃之间（表1），平均值为57.41 ℃。根据HLT50数值，参试8个绣球品种耐热性排序依次为：‘头花’>‘花手鞠’>‘小町’>‘爱莎’>‘含羞叶’>‘纱织小姐’>‘银边’>‘灵感’。
2.2 不同绣球品种的热害指数观测
热害指数（HII）反映了不同绣球品种的耐热能力，热害指数越低则植株耐热性相对越强，此方法简单易行，可作为评价植物耐热性强弱的指标之一[22]。高温胁迫对绣球品种茎叶生长势影响显著，各品种的热害指数均显著升高。结果表明，43 ℃高温胁迫6 d后，8个品种形态与对照组相比均发生了较为明显的变化，各品种叶片均表现出不同程度的叶缘卷曲、褐变失绿、失水皱破、萎蔫干枯、脱落等热害症状，严重影响其观赏价值，其中品种‘纱织小姐’与‘银边’全株超过3/4的叶片发生褐变、老叶枯萎甚至脱落，受热害程度最严重;品种‘花手鞠’与‘头花’受热害程度相对较轻，表现出较强的耐热性。各品种根据热害指数由到高到低次序为：‘纱织小姐’>‘银边’>‘灵感’>‘小町’>‘爱莎’>‘含羞叶’>‘花手鞠’≈‘头花’，热害指数分别为：14.4%、13.6%、12.8%、8.8%、7.2%、6.4%、5.6%與5.6%。
2.3 绣球品种的耐热性综合评价
2.3.1 各项指标的耐热系数及相关性分析 将43 ℃/33 ℃高温胁迫6 d后的各项指标测定结果进行耐热系数α值转化（表2）。由表2可知，各绣球品种经高温胁迫后的CMP、MDA含量、POD活性、CAT活性、SS含量、SP含量与Pro含量与对照相比均有提升（α>100%）;RWC与对照相比有所下降（α<100%）;Chla含量、Chlb含量、Chla+Chlb含量、Chla/Chlb值、Car含量与SOD活性在品种‘含羞叶’‘小町’‘爱莎’‘银边’‘纱织小姐’和‘灵感’中与对照相比有所下降，在品种‘头花’‘花手鞠’中与对照相比有所提升。此外，14项耐热指标在经过耐热系数转化后变化幅度差异显著，其中变化幅度最大的指标为POD活性，变化幅度最小的为RWC，并且不同品种在同一指标内的变化差异较大。 从表3可见，RWC与SOD活性之间存在极显著相关关系，与Pro含量之间存在显著相关关系;Chla含量与Chlb含量、Chla/Chlb及Chla+Chlb含量之间存在极显著相关关系，与SOD活性、SS含量及Pro含量之间存在显著相关关系;Chlb含量与Chla+Chlb含量之间存在极显著相关关系;Chla/Chlb与Chla+Chlb含量、SS含量及Pro含量之间存在显著相关关系;Chla+Chlb含量与SOD活性、SS含量及Pro含量之间存在显著相关关系;Car含量与SP含量及Pro含量之间存在极显著相关关系，与SOD活性的之间存在显著相关关系;SOD活性与Pro含量之间存在极显著相关关系;POD活性与SS含量之间存在显著相关关系;SS含量与Pro含量之间存在显著相关关系;SP含量与Pro含量之间存在显著相关关系。其余各项指标之间关系不显著（表3）。可见，各单项指标之间普遍存在信息重叠现象，可能会在一定程度上影响绣球品种耐热性鉴定结果。
2.3.2 不同绣球品种各单项指标的主成分分析 主成分分析法能通过降维的方式将多个相互关联的指标转化为少数几个包含原变量大部分信息且相互独立的综合指标。特征值能在一定程度上反映主成分的影响力，当特征值大于1则可作为较好的评价依据。通过对高温胁迫下绣球品种各项指标的耐热系数进行主成分分析，前4个主成分的特征值大于1，把原本相互关联的13个单项指标重新降维成4个新的互不干扰的综合指标，依次用Z1～Z4表示，4个综合指标的特征值分别为7.299、1.987、1.589和1.195，贡献率分別为52.134%、14.189%、11.351%和8.533%，累计贡献率达86.207%（当累积贡献率超过80%则具有较强数据代表性），说明所提取的这4个综合指标基本上能代表不同绣球品种耐热性原始数据的绝大部分信息（表4）。
主成分中各项指标的荷载值大小能反映出该指标在主成分中的关键程度。对各综合指标进行特征向量分析，其中Z1中Chla含量、Chla+Chlb含量和Pro含量有较大的正载荷量;Z2中SP含量、Car含量和CMP有较大的正载荷量，Chlb含量和CAT活性有较大负荷载量;Z3中MDA含量、CAT活性与SP含量有较大的正载荷量;Z4中POD活性、Chla含量与Chlb含量有较大的正载荷量，RWC与SS含量有较大负荷载量。综合指标的贡献率的大小能反映不同综合指标的相对重要性，
主成分Z1的贡献率最大（表4），说明主成分Z1所涵盖的信息能更全面地反映8个绣球品种的耐热能力。
2.3.3 不同绣球品种耐热性综合分析 将8个绣球品种的各项生理指标的耐热系数与其主成分的荷载值相结合，根据线性组合方程计算结果得出4个主成分的得分值（comprehensive index，CI），得分值越高说明该品种在这一综合指标中的耐热表现越好。由表5可知，品种‘头花’与‘花手鞠’在第1主成分上的耐热性最强;品种‘含羞叶’与‘小町’在第2主成分上的耐热性最强;品种‘含羞叶’与‘纱织小姐’在第3主成分上。
2.3.4 不同绣球品种耐热性聚类分析 聚类分析法是一种能根据对象性质划分为不同集合的样本归类方法[26]。由图1可知，本研究根据各个材料的耐热指标综合评价值（D值），采用平方欧氏距离，用离差平方和法对不同绣球品种的耐热性进行了系统聚类分析，聚类结果在欧式距离5处可将绣球耐热性分为3类：中等耐热型品种、不耐热型品种、耐热型品种。其中，中等耐热型品种包括‘爱莎’‘含羞叶’‘小町’和‘纱织小姐’;不耐热型包括品种‘银边’与‘灵感’;耐热型包括品种‘花手鞠’与‘头花’。8个绣球品种的耐热性鉴定结果与热害指数鉴定结果、高温半至死温度鉴定结果大体一致。
2.3.5 耐热鉴定指标的筛选 为探究各项指标与耐热性之间的关系，在保证评价结果准确性的同时精简评价指标，建立可用于绣球品种耐热性评价的数学模型，将8个绣球品种的综合评价D值作为因变量，各单项指标的耐热系数作为自变量进行逐步回归分析，得出预测绣球品种耐热能力
的最优回归方程：D=–0.707–0.000068X1+ 0.001408X2+0.004131X3+0.003451X4+0.003746X5–0.000398X6–0.000038X7（R2=1，P=0.0001）。由方程可知，原有的14个指标被精简为X1～X7 7个指标，它们分别代表Pro含量、Chla+Chlb含量、Car含量、Chla/Chlb、SOD活性、RWC与SS含量的耐热系数α，可用于鉴定不同绣球品种的耐热性，即在相同条件下测定其他品种的这7个指标，并计算其耐热系数，进而利用该方程预测相应品种的耐热性。
用该回归方程对8个品种的耐热性进行预测，即把8个品种的上述7个指标的耐热系数（α值）代入该方程，得到各品种的预测值：‘含羞叶’为0.435，‘头花’为0.828，‘小町’为0.308、‘爱莎’为0.437，‘银边’为0.126，‘花手鞠’为0.665，‘纱织小姐’为0.293，‘灵感’为0.088。将各品种的预测值与综合评价值进行相关性分析，二者相关系数R2=1，在0.1水平上显著相关。说明此方程能够很好地量化各品种的综合耐热能力与各指标值间的因果关系，对绣球品种间耐热性的预测效果较好。
此外，将各品种的耐热性综合评价D值与高温半致死温度HLT50进行相关性分析，发现二者成极显著正相关（R2=0.868）;将各品种的耐热性综合评价D值与热害指数HII进行相关分析，发现二者呈显著负相关（R2=–0.743），即测定HLT50与HII也可作为鉴定绣球品种耐热性的可靠方法。
3 讨论
相关研究表明，植物细胞膜的热稳定性与其耐热能力密切相关[27]，通过电导率法结合Logistic方程拟合估算HLT50，能较为准确地反映植物耐性差异。通过该方法评价植物耐热性既便捷又高效[28-29]，现已应用在杜鹃花科（Ericaceae）[30]、蔷薇科（Rosaceae）[31]及景天科（Crassulaceae）[32]等多种植物材料的耐热性鉴定中，并已有学者验证该方法在多个经济作物与观赏花卉品种耐热性差异评价中的可行性[33-34]。但也有学者认为由于该实验方法在离体条件下进行，因而并不能准确反映出品种间的耐热差异[33]。本研究对8个绣球品种的HLT50进行测定，温度范围在55.10～59.02 ℃之间，测定结果与耐热性综合评价结果呈极显著正相关（R2=0.868），验证了其作为绣球品种耐热评价指标的可靠性。刘婉迪等[30]测定出9个杜鹃品种（Rhododendron simsii）的HLT50范围在46.92～53.84 ℃之间，刘志高等[18]对20个铁线莲品种的HLT50进行测定，HLT50最高的品种仅为42.17 ℃，其耐热能力均不如本研究参试的8个绣球品种。从耐热性角度考虑，8个绣球品种的半致死温度均超过了55 ℃，说明这8个品种都有一定的高温耐受能力，是华南地区园林绿化材料更新潜在的植物资源，但其耐热能力与凌霄（Campsis grandiflora）、地锦（Parthenocissus inserta）、牛轭草（Murdannia loriformis）和细竹篙草（Murdannia simplex）等平均HLT50在60 ℃以上的观赏植物相比仍存在差距[35-36]。这一方面体现了不同植物种类间的耐热性差异，另一方面反应出绣球品种对高温的耐受能力有限，在实际应用中需要尽量配置在通风良好的遮荫环境下。此外，各品种的耐热性综合评价D值与HLT50呈极显著正相关（R2=0.868），说明绣球的高温半致死温度能较好地反映其耐热性，可作为其耐热性评价的关键指标。 持续高温会造成绣球的叶片失绿褐变、萎蔫下垂甚至枯萎，这是由于高温引起植物体内的酶代谢失活最终导致叶绿素合成减少，或是叶片细胞含水量降低的表现[16，37]。HII通过将高温胁迫下植物的表型特征数量化直观地反映受热害程度，被作为可靠评价指标广泛地应用在耐热性评价研究中[17-18]。本研究结果表明，高温胁迫下各绣球品种皆表现出不同程度的热害症状，其中‘花手鞠’与‘头花’在高温环境下形态表现较好，为耐热品种;‘纱织小姐’与‘银边’受热害程度较高，为热敏品种，该鉴定结果与耐热综合评价D值聚类鉴定结果基本一致。此外，不同耐热型绣球品种的HII差异显著，各品种的耐热性综合评价D值与热害指数HII呈显著负相关（R2=–0.743），进一步验证了以HII作为绣球耐热性评价关键指标的可靠性，为通过表型观测对品种耐热性进行预判提供了依据。
植物能利用自身体内的一系列生理代谢途径来抵御或适应高温胁迫，比如通过产生兼容溶质组织蛋白质与细胞结构，通过渗透调节保持细胞充盈，通过抗氧化系统重建细胞氧化还原平衡与内稳态等[38]。不同植物种类对高温胁迫的动态响应机制存在较大差异，植物的耐热能力不仅与高温胁迫的程度与胁迫时间有关，也受到植物的遗传背景、表型性状、生长发育阶段与生理代谢水平等多方面因素的影响，是在植物与各环境因子之间相互作用的复杂过程中形成的[39]。本研究采用人工气候室模拟高温胁迫环境，有效地避免了干旱、强光、高湿等复合逆境的影响，便于更加客观地对绣球品种的耐热性进行评价。可用于观赏植物耐热性鉴定的指标，不同指标对高温胁迫的敏感程度表现不一，且各指标间存在着复杂的关联性[40-41]。本研究选用的14项生理生化指标在高温胁迫下的响应方式各异，其相关性分析结果也表明各单项指标所提供的耐热响应信息发生不同程度的交叉重叠，这与上述观点一致，说明根据其中任何一项指标对绣球品种的耐热性进行评价都将得到不同的评价结果，因此无法直接通过相互矛盾的单一指标对绣球各品种的耐热性进行评价，耐热性作为复杂的综合性状需要依靠多元的统计分析方法，建立高效、科学的绣球耐热性综合评价体系[42]。
为此，本研究采用主成分分析法结合隶属函数分析法进一步将各指标进行综合分析，以期弥补单项指标评价耐热性的缺陷，得到更为准确的耐热性评价结果。通过测定对照及高温胁迫处理下不同绣球品种的14个单项指标，以各项指标的耐热系数α为基础，利用主成分分析法将原有的14个单项指标降维成4个互不干扰的综合指标，再利用隶属函数分析法获得不同绣球品种的耐热性综合评价D值。由于D值不仅与各指标间的关联性相关，也与各指标的权重相关，所以利用其对绣球耐热性评价更具有系统性和科学性。根据综合评价D值进行聚类分析，可将各参试绣球品种划分为耐热型、中等耐热型和不耐热型3个类群，耐热型品种包括‘花手鞠’与‘头花’，其耐热表现较好，是应用在南方夏季高温地区的潜在优良品种;中等耐热型品种包括‘爱莎’‘含羞叶’‘小町’与‘纱织小姐’;不耐热型品种包括‘银边’与‘灵感’，其耐热表现较差，更适合栽植在温度不高的温带地区。
此外，本研究利用逐步回归法建立了能够快速、准确鉴定绣球品种耐热性评价模型，筛选出Pro含量、Chla+Chlb含量、Car含量、Chla/Chlb、SOD活性、RWC与SS含量等7个显著影响绣球品种耐热性单项指标。经过预测值验证，结果表明由这7项指标构建的评价模型能能夠很好地量化各品种的综合耐热能力与各指标值间的因果关系，进而高效、准确地预测不同绣球品种耐热能力的强弱，可作为绣球品种耐热性鉴定的关键指标，进一步为绣球品种耐热性鉴定工作简化了工作量。
在本研究的基础上，更多研究可以进一步展开，如通过解剖技术进一步探究绣球叶片显微结构与其耐热性的联系;选择耐热性强的品种作为亲本培育新种质展开育种工作;随着分子生物技术的成熟，对植物的耐热性研究也已达到分子水平，选择耐热与不耐热绣球品种构建基因文库，筛选与耐热性直接相关的关键基因，通过基因与表达产物进一步揭示绣球对高温的响应机制，导入外源基因进行耐热性育种是未来研究的方向之一。此外，绣球由于其极高的观赏价值与经济价值，在园林美化、盆花销售与切花生产中扮演越来越重要的角色，但其生产栽培仍受不同地区气候差异的限制，故除研究品种本身的抗逆性强弱和响应机制外，探索更多有助于提升其抗性的方法，从而降低栽培管理成本并扩大生产应用区域也是较好的切入点。
参考文献
乔 谦， 王江勇， 陶吉寒. 绣球属植物研究进展[J]. 农学学报， 2020， 10（4）： 60-64.
李艳香. 19个绣球属种质资源的SRAP分析[D]. 长沙： 湖南农业大学， 2009.
彭尽晖， 周朴华， 周红灿， 等. 湖南省绣球属植物资源调查[J]. 湖南农业大学学报（自然科学版）， 2008（5）： 563-567.
孙 强， 虞秀明， 姚红军. 绣球属品种资源收集及扦插生根能力比较[J]. 北方园艺， 2016（11）： 71-73.
张国华， 陈旦旦， 曹受金， 等. 绣球花扦插繁殖技术研究[J]. 现代园艺， 2019， 42（23）： 24-25.
王红梅， 刘春风， 孟紫莲， 等. 圆锥绣球‘北极熊’组培快繁体系的建立[J]. 分子植物育种， 2021， 19（7）： 2305-2311.
Uemachi T， Nishio T. Relationship between the inflorescence structure and setting of decorative flowers on lacecap cultivars of Hydrangea macrophylla （Thunb.） Ser. and H. serrata （Thunb.） Ser[J]. Horticultural Research， 2005， 4（4）： 435-438. 章 毅， 韦孟琪， 孙欧文， 等. 不同绣球品种对干旱胁迫的生理响应及抗旱机制研究[J]. 西北林学院学报， 2018， 33（1）： 90-97.
蔡建国， 章 毅， 孙欧文， 等. 绣球抗旱性综合评价及指标体系构建[J]. 应用生态学报， 2018， 29（10）： 3175-3182.
李叶华， 陈 爽， 赵 冰. 铝胁迫对绣球组培苗生长及生理特性的影响[J]. 浙江农林大学学报， 2020， 37（6）： 1064-1070.
汪 霖. NaCl胁迫对绣球生理指标的影响[J]. 北方园艺， 2011（24）： 80-83.
辛雅芬， 石玉波. 高温对八仙花生理指标的影响[J]. 安徽农业科学， 2009， 37（11）： 5204-5206.
蒋昌华， 李 健. 高温胁迫对八仙花品种生理特性的影响[C]//中国植物园（第十五期）， 西安： 中国林业出版社， 2011： 196-198.
Gulen H， Eris A. Effect of heat stress on peroxidase activity and total protein content in strawberry plants[J]. Plant Science， 2003， 166（3）： 739-744.
Tian Z G， Yang Y， Wang F. A comprehensive evaluation of heat tolerance in nine cultivars of marigold[J]. Horticulture Environment and Biotechnology， 2015， 56（6）： 749-755.
Stirbet A， Govindjee. On the relation between the Kautsky effect （chlorophyll a fluorescence induction） and photosystem II： Basics and applications of the OJIP fluorescence transient[J]. Journal of Photochemistry and Photobiology，2011， 104（1-2）： 236-257.
張佳平， 李丹青， 聂晶晶， 等. 高温胁迫下芍药的生理生化响应和耐热性评价[J]. 核农学报， 2016， 30（9）： 1848-1856.
刘志高， 邵伟丽， 申亚梅， 等. 铁线莲品种耐热性分析及评价指标筛选[J]. 核农学报， 2020， 34（1）： 203-213.
李 辉. 杜鹃耐热资源评价及其生理机制研究[D]. 扬州： 扬州大学， 2019.
彭勇政， 刘智媛， 朱晓非， 等. 5个月季品种高温处理后生理指标变化及其耐热性评价[J]. 上海交通大学学报（农业科学版）， 2019， 37（5）： 53-58.
李合生. 植物生理生化实验原理和技术[M]. 北京： 高等教育出版社， 2004.
张佳平. 芍药在杭州栽培的耐热评价及地下芽休眠机理研究[D]. 杭州： 浙江大学， 2015.
张旻桓. 湖南牡丹资源遗传多样性及耐热性研究[D]. 长沙： 中南林业科技大学， 2019.
田治国. 万寿菊属植物耐热性与抗旱性的评价及生长生理特性的研究[D]. 杨凌： 西北农林科技大学， 2012.
赵 伟. 三种国产唇柱苣苔属野生花卉的组培、扦插和幼苗半致死温度研究[D]. 上海： 上海师范大学， 2015.
尹立为. 函数数据降维方法的分析与研究[D]. 长春： 吉林大学， 2014.
Ashraf M， Saeed M M， Qureshi M J. Tolerance to high temperature in cotton （Gossypium hirsutum L.） at initial growth stages[J]. Environmental and Experimental Botany， 1994， 34（3）： 275-283.
史彦江， 宋锋惠， 俞 涛， 等. 高温胁迫下平欧杂种榛品种的耐热性研究[J]. 中国农学通报， 2015， 31（28）： 23-28.
Castor T P， Tyler T A. Determination of taxol in taxus media needles in the presence of interfering components[J]. J Liq Chromatogr Relat Technol， 1993， 16（3）： 723-731.
刘婉迪， 王 威， 谢 倩， 等. 9个杜鹃品种的高温半致死温度与耐热性评价[J]. 西北林学院学报， 2018， 33（5）： 105-110， 136.
吕晓倩. 梅花杂交育种与后代耐热性评价研究[D]. 北京： 北京林业大学， 2014.
张燕利， 高捍东， 吴锦华. 4种景天科植物耐热性测定[J]. 西南林业大学学报， 2010， 30（6）： 52-54.
Yeh D M， Lin H F. Thermostability of cell membranes as a measure of heat tolerance and relationship to flowering delay in Chrysanthemum[J]. Journal of the American Society for Horticultural Science American Society for Horticultural Science， 2003， 128（5）： 656-660. 薛盼盼， 张 强， 魏学智. 模拟干旱胁迫对酸枣种子萌发及活力的影响[J]. 西北农业学报， 2016， 25（12）： 1837-1843.
徐玉芬， 杨 静， 邓 磊， 等. 应用Logistic方程测定4种鸭跖草科屋顶绿化植物的耐热性[J]. 仲恺农业工程学院学报， 2018， 31（4）： 52-55， 70.
袁雪涛， 夏 冬， 许永利， 等. 应用Logistic方程测定四种藤本植物的耐热性[J]. 北方园艺， 2017（14）： 99-102.
Wang Q， Cheng T R， Yu X N， et al. Physiological and biochemical responses of six herbaceous peony cultivars to cold stress[J]. South African Journal of Botany， 2014， 94： 140-148.
Valliyodan B， Nguyen H T. Understanding regulatory networks and engineering for enhanced drought tolerance in plants[J]. Current Opinion in Plant Biology， 2006， 9（2）： 189-195.
Khanna-Chopra R， Chauhan S. Wheat cultivars differing in heat tolerance show a differential response to oxidative stress during monocarpic senescence under high temperature stress[J]. Protoplasma， 2015， 252（5）： 1241-1251.
Zhang Y T， Du H M. Differential accumulation of proteins in leaves and roots associated with heat tolerance in two kentucky bluegrass genotypes differing in heat tolerance[J]. Acta Physiologiae Plantarum， 2016， 38（9）： 1-15.
Cheng Z Y， Sun L， Wang X J， et al. Ferulic acid pretreatment alleviates heat stress in blueberry seedlings by inducing antioxidant enzymes proline and soluble sugars[J]. Biologia Plantarum， 2018， 62（3）： 534-542.
Patton A J， Cunningham S M， Volenec J J， et al. Differences in freeze tolerance of zoysiagrasses： II. carbohydrate and proline accumulation[J]. Crop Science， 2007， 47（5）： 2170-2181.
責任编辑：沈德发