-
低温是限制植物生长发育和影响植物地理分布的主要环境因素之一[1]。由于纬度和地形的差异,低温限制成为优良林木品种引种和培育迫切需要解决的难题[2-3]。植物通过自身的生理和形态的可塑性来提高在低温条件下的适应性[1,4],如抗氧化酶系统(超氧物歧化酶、过氧化物酶、过氧化氢酶)将在低温条件下被激活以抵御超氧自由基的毒害[5-6];植物组织中可溶性蛋白、非结构性碳水化合物以及游离脯氨酸含量的增加也可提高细胞液渗透势,降低冰点,从而有效地避免细胞液结冰[7-8]。由于树种对低温的耐受能力不同,了解目标树种的耐寒极限在林木培育和良种引进过程中尤为重要。当植物受到低温胁迫时,电解质因细胞膜被破坏而大量外渗,通过测定组织电导率可判定植物受低温伤害程度[9-11]。然而,判断引进树种是否适合引种地栽培仍然需要根据生理指标进行综合判定[12]。
雷竹(Phyllostachys praecox C. D. Chu et C. S. Chao ‘Prevernalis’)隶属禾本科(Gramineae)刚竹属(Phyllostachys Sieb. et Zucc.),是早竹的栽培变型。雷竹易栽培、产量高、竹笋品质好,是我国长江以南地区重要的笋用竹种,具有较高的经济价值和生态服务价值[13];同时,具有特殊的秆形、叶形、秆色的雷竹变型也极具观赏价值;但选育优良观赏雷竹品种并推广栽培需要对雷竹的生长习性、抗逆能力进行综合评价。目前,关于观赏竹对环境因子的适应性研究主要包括温度[14-15]、水分[16-17]、养分[18-19]等方面,通过综合评价的方法筛选出一些抗逆性较好的观赏竹品种,如耐旱性较好的黄秆京竹(Ph. aureosulcata McClure ‘Aureocarlis’),耐寒性较好的平安竹(Pseudosasa japonica var. tsusumiana Yanagita)和金镶玉竹(Ph. aureosulcata McClure ‘Spectabilis’)等。关于雷竹适应性的研究主要包括经营措施[20]、养分[21]、水分[22]、土壤盐碱度[23]等方面,而针对观赏雷竹品种耐寒性的研究还相对匮乏。因此,本研究对5个具有良好观赏性表型的雷竹变型进行低温胁迫试验,采用Logistic方程拟合、主成分分析等方法对其耐寒性进行综合评价,以期为观赏性雷竹的引种和培育提供参考。
HTML
-
试验材料取自浙江省杭州临安太湖源观赏竹种园有限公司竹种园(119°95′ E,29°48′ N),该地地处浙江省西北部,亚热带季风气候区南缘,温暖湿润,光照充足,年均降水量1613.9 mm,年均气温13~22℃,无霜期平均为237 d,受台风、寒潮和冰雹等灾害性天气影响。本研究的5个观赏雷竹变型均于2016年引种栽培,基本信息见表1。每个品种以同样的种植密度定植于20 m2的样地内。采用同样的田间管理方式(灌溉、除草等),以保证各品种生长条件基本一致。
变种
Varieties引种来源
Source表型
Phenotype株高
Plant height/
m地径
Basal diameter/
mm胸径
Diameter at breast
height/mm枝下高
Under branch
height/m观赏性表型
Ornamental phenotype黄条雷竹
Ph. praecox f. Notata浙江临安
Lin’an, Zhejiang5.28 ± 1.17 34.7 ± 7.87 33.5 ± 0.75 2.16 ± 0.82 竹秆主体为绿色,节间分枝一侧沟槽内具黄色条纹
The main body of bamboo stalk is green with yellow stripe in one side groove of internode branch黄皮雷竹
Ph. praecox f. Chrysoderma浙江临安
Lin’an, Zhejiang2.30 ± 0.44 25.9 ± 5.13 22.0 ± 6.49 1.14 ± 0.15 竹秆和竹枝为黄色,基部节间偶有绿色纵条纹
Bamboo stalks and branches are yellow, with occasional green longitudinal stripes between basal nodes弯秆雷竹
Ph. praecox cv. Linanensis浙江临安
Lin’an, Zhejiang3.03 ± 0.75 26.3 ± 0.11 24.1 ± 12.2 1.32 ± 0.45 竹节节间呈一定幅度弯曲
The internodes of bamboo are bent at a certain degree花秆雷竹
Ph. praecox f. Viridisulcata浙江临安
Lin’an, Zhejiang5.29 ± 90.5 41.3 ± 9.75 37.1 ± 11 2.18 ± 0.7 竹秆黄色,间分枝一侧沟槽为绿色,且沟槽外有绿色细纵条纹;部分叶片淡白色纵条纹;竹箨淡黄褐色
The bamboo stalk is yellow, the groove on the side of the branch is green, and the bamboo pole except the groove has green and fine longitudinal stripes; some leaves are light white and longitudinal stripes; the bamboo stem is yellowish brown细叶乌头雷竹
Ph. praecox ‘Xiyewutouleizhu’浙江临安
Lin’an, Zhejiang2.74 ± 0.49 39.7 ± 9.47 38.88 ± 9.95 1.47 ± 0.38 叶细而窄呈瓦状上翘,边缘略反转,笋壳乌褐色
The leaves are thin, narrow and upturned like tiles. The upward reversed slightly at the edge of the leaf. The shell of bamboo shoot is dark brown.Table 1. Basic information of 5 varieties
-
于2019年9月随机采集5个观赏雷竹变型标准木的功能叶10 g,置于0~2℃冰盒中带回实验室。将叶片用清水洗净并用吸水纸擦干表面水分备用。试验设置0、−5、−10、−15、−20、−25、−30℃ 7个温度梯度,将等量的叶片分为 7个试验组分别进行不同温度处理,处理时间24 h,每个试验组3次重复。用超低温冷冻循环仪(DX-4015,北京长流科学仪器有限公司)以5℃·h−1逐渐降温到设置温度,处理完毕后,以5℃·h−1逐渐升温至0℃,然后在0℃冰箱内保存待测,待测时间小于1 h。
叶片细胞膜损害程度测定:(1)相对电导率(REC):取待测叶片用打孔器制成大小相同的圆片,随机选取形状完整的20片,放置于试管底部,加入20 mL去离子水,然后将试管放入真空干燥箱中用真空抽气3次,每次20 min,常温下放置6 h,期间每30 min 振荡一次试管,用电导率仪(STARTER3000C,浙江纳德科学仪器有限公司)测定电导率值S1,再将试管置于水浴锅中沸水浴40 min 后冷却,在常温下稳定30 min后测定电导率值S2。以去离子水电导率值S0作为对照。根据公式REC=(S1−S0)/(S2−S0)计算相对电导率。(2)丙二醛(MDA)含量用硫代巴比妥酸(TBA)显色法测定。
叶片抗氧化酶活性测定:超氧化物歧化酶(SOD)采用氮蓝四唑法(NBT法)测定;超氧化物酶(POD)的活性采用愈创木酚法测定;过氧化氢酶(CAT)采用紫外吸收法测定。
渗透调节物质测定:分别测定可溶性蛋白和脯氨酸含量分析雷竹在低温条件下的生理可塑性。可溶性蛋白含量采用考马斯亮蓝染色法测定;脯氨酸含量采用酸性茚三酮比色法测定。以上生理指标均参照李合生的植物生理试验方法[24]。
-
运用SPSS17.0软件将相对电导率拟合Logistic方程y = k/(1+ae−bx),后将获得的参数计算后得出低温半致死温度(LT50)。公式中x为处理温度,y代表相对电导率,k、a、b为参数,k为y的最大极限值,b反映了x与y之间的对应关系,a表示曲线对原点的相对位置,如果y为累计频率,则常数k为1[25]。当方程曲线出现拐点时的温度即为引起膜不可逆伤害的临界点,可用来确定植物组织的LT50,从而定量判断5个雷竹变型的抗寒能力[26]。为了方便对 a、b 求解,对原拟合方程进行线性化:令y1 = ln[(k−y)/y],则y1 = −bx+lna,从而将Logistic方程转化为 y1与处理温度x之间的直线方程,再运用SPSS 17.0软件进行线性拟合可求出 a、b值以及相关系数R,而 LT50=(lna)/b[27]。
-
根据模糊数学的隶属函数法计算5个雷竹变型耐寒性的隶属函数值,从而得到综合排名。隶属函数值可根据公式求出:当指标与耐寒性呈正相关时(如抗氧化酶活性、可溶性蛋白和脯氨酸含量),Uijk= (Xijk−Xmin
)/(Xmax−Xmin );而指标与耐寒性呈负相关时(如相对电导率和MDA含量),Uijk = 1−(Xijk−Xmin)/(Xmax−Xmin),其中,i代表竹种,j代表温度,k代表指标,Uijk为第i个竹种在第j个温度处理下第k项指标的隶属函数值,Xmin 表示所有竹种所有温度下第k项指标的最小值,Xmax表示第k项指标的最大值[2]。 -
对所有生理指标用SPSS 17.0软件进行双因素方差分析(ANOVA),采用Canoco 5.0软件进行主成分分析(PCA)和冗余分析(RDA),确定生理指标变化的主分量,计算其载荷值,并分析各生理指标与温度间的相关性。根据主成分分析结果计算生理参数的主成分载荷量,从而确定影响雷竹耐寒性的主要生理参数。运用Origin 9.0软件进行图表绘制。
1.1. 试验材料
1.2. 试验方法
1.3. 数据处理
1.3.1. Logistic方程的建立和LT50确定
1.3.2. 隶属函数值的计算
1.3.3. 差异性分析(ANOVA)及主成分分析(PCA)
-
由图1可知:5个雷竹变型叶片相对电导率(REC)均随着温度的降低呈先上升后下降的趋势;但弯秆雷竹叶片相对电导率在−20℃达到峰值,而其余4种雷竹则均在−25℃达到峰值。另外,5个雷竹变型相对电导率的峰值差异显著,其中,花秆雷竹相对电导率最高,而细叶乌头雷竹相对电导率最低。根据相对电导率的结果拟合Logistic方程,曲线的拟合效果均达到极显著水平(P < 0.01)。根据a、b值分别计算出5个雷竹变型的半致死温度(LT50)。结果显示:弯秆雷竹的LT50最低(−18.73℃),而花秆雷竹的LT50最高(−14.57℃)(表2)。
-
由图2可知:5个雷竹变型叶片超氧化物歧化酶(SOD)活性随处理温度的降低均呈先上升后下降的趋势。花秆雷竹SOD活性在−10℃时达到峰值,而其余4种雷竹均在−15℃时达到峰值。5个雷竹变型SOD活性上升幅度为:黄条雷竹(590.72)> 黄皮雷竹(548.52)> 弯秆雷竹(471.17)> 细叶雷竹(428.97)> 花秆雷竹(393.81)。双因素方差分析发现,SOD活性在不同温度间和不同变种间差异极显著(P < 0.01)(表3)。
变种
varietieslogistic方程 a b LT50/℃ R2 黄条雷竹
Ph. praecox f. Notatay = 100/(1+2.35e−0.138 7) 10.53 −0.138 7 −16.98 0.927** 黄皮雷竹
Ph. praecox f. Chrysodermay = 100/(1+6.63e−0.112 5) 6.63 −0.112 5 −16.82 0.880** 弯秆雷竹
Ph. praecox cv. Linanensisy = 100/(1+4.62e−0.081 7) 4.62 −0.081 7 −18.73 0.725** 花秆雷竹
Ph. praecox f. Viridisulcatay = 100/(1+6.63e−0.129 9) 6.63 −0.129 9 −14.57 0.901** 细叶乌头雷竹
Ph. praecox ‘Xiyewutouleizhu’y = 100/(1+6.59e−0.110 3) 6.59 −0.110 3 −17.09 0.943** 注:**表示相关性在P < 0.01水平上显著。
Note: ** represents a significant correlation at 0.01 level.Table 2. Regression equation of conductivity under low temperature treatment and LT50 of 5 varieties
参数
Parameters处理
Treatments变种
Varieties处理 × 变种
Treatments × VarietiesF-value P-value F-value P-value F-value P-value 相对电导率 REC 1 087.349 < 0.01 12.725 < 0.01 9.718 < 0.01 超氧化物歧化酶 SOD 47.793 < 0.01 17.765 < 0.01 0.717 0.818 超氧化物酶 POD 19.793 < 0.01 0.926 0.454 1.025 0.449 过氧化氢酶 CAT 50.826 < 0.01 3.874 < 0.01 0.595 0.922 可溶性蛋白 Soluble protein 30.314 < 0.01 18.696 < 0.01 0.428 0.989 丙二醛 MDA 225.034 < 0.01 33.078 < 0.01 0.926 0.568 脯氨酸 Proline 86.603 < 0.01 43.230 < 0.01 0.576 0.934 Table 3. Two-way variance analysis on varieties and temperatures
超氧化物酶(POD)和过氧化氢酶(CAT)活性随温度变化的趋势与SOD的类似。花秆雷竹的POD活性在−10℃时达到极值,而其余4种雷竹均在−15℃时达到峰值。5个雷竹变型POD活性上升幅度为:黄皮雷竹(2 577.78)> 黄条雷竹(2 266.67)> 细叶乌头雷竹(2 155.56)> 弯秆雷竹(2 133.33)> 花秆雷竹(1 555.56)(图3)。5个雷竹变型CAT活性均在−15℃时达到峰值,其上升幅度为:弯秆雷竹(540.74)> 黄条雷竹(511.11)> 细叶乌头雷竹(488.89)> 黄皮雷竹(481.48)> 花秆雷竹(437.04)(图4)。
-
由图5可知:5个雷竹变型叶片中丙二醛(MDA)含量随处理温度的降低均呈先上升后下降的趋势。花秆雷竹的MDA含量在−15℃时达到峰值,而其余4种雷竹均在−20℃时达到峰值。5个雷竹变型MDA含量上升幅度为:花秆雷竹(28.71)> 黄条雷竹(24.9)> 细叶乌头雷竹(24.8)> 弯秆雷竹(24.27)> 黄皮雷竹(23.59)。双因素方差分析发现,MDA含量在不同温度间和不同变种间差异极显著(P < 0.01),但5个雷竹变型叶片MDA含量在达到峰值后的下降程度很小,其在−30℃时的含量仍远高于0℃时的含量。
-
由图6可知:5个雷竹变型叶片中可溶性蛋白含量均随处理温度的降低均呈先上升后下降的趋势。弯秆雷竹可溶性蛋白含量在−20℃时达到峰值,而其余4种雷竹均在−15℃时达到峰值。5个雷竹变型可溶性蛋白含量的上升幅度为:黄皮雷竹(0.236)> 细叶乌头雷竹(0.229)> 弯秆黄皮雷竹(0.221)> 黄条雷竹(0.219)> 花秆雷竹(0.179)。双因素方差分析发现,可溶性蛋白含量在不同温度间和不同变种间差异极显著(P < 0.01)。
由图7可知:5个雷竹变型叶片中脯氨酸含量均随处理温度的降低均呈先上升后下降的趋势,且均在−15℃时达到峰值。5个雷竹变型脯氨酸含量的上升幅度为:弯秆雷竹(56.08)> 黄条雷竹(55.35)> 黄皮雷竹(54.35)> 细叶乌头雷竹(50.56)> 花秆雷竹(47.09)。双因素方差分析发现,可溶性蛋白含量在不同温度间和不同变种间差异极显著(P < 0.01)。
-
对低温胁迫下的5个雷竹变型生理特征的主效应分析(PCA)结果(表4)显示:第1主成分可反映生理指标变化59.38%的贡献率,其中,可溶性蛋白和脯氨酸这两个体现细胞渗透调节响应的指标载荷量较高;第2主成分可解释约27.8%的信息量,其中,体现细胞抗氧化酶调节能力的指标SOD和POD载荷量较高。从温度对生理指标的影响看,可将0℃、−5℃归为一类,−10℃、−15℃和−20℃归为一类,−25℃和−30℃归为一类,同类的处理对生理指标产生的影响相近(图8)。冗余分析(RDA)结果(图9)表明:3种抗氧化酶活性、可溶性蛋白以及脯氨酸间表现出显著正相关关系,而MDA和相对电导率间表现出显著正相关关系。温度的改变对相对电导率、脯氨酸、MDA以及抗氧化酶活性的影响较大,而对可溶性蛋白的影响相对较小。运用隶属函数法将与耐寒性相关(正相关和负相关)的生理指标综合分析比单独分析一个或一类指标更具科学性。根据隶属函数值得到5个雷竹变型耐寒性综合排名为:弯秆雷竹 > 细叶乌头雷竹 > 黄皮雷竹 > 黄条雷竹 > 花秆雷竹(表5)。
生理指标
Physiological indices第1主成分
PC1第2主成分
PC2相对电导率 REC −0.974 0.078 超氧化物歧化酶 SOD 0.620 0.753 超氧化物酶 POD 0.733 0.604 过氧化氢酶 CAT 0.793 −0.462 丙二醛 MDA −0.942 0.325 可溶性蛋白 Soluble protein 0.947 −0.247 脯氨酸 Proline 0.925 0.077 Table 4. Component matrix of physiological indices
Figure 8. Principal component analysis(PCA)of physiological parameters of 5 varieties under different treatments
Figure 9. Redundant analysis of temperature response by physiological indexes of Ph. praecox ‘Prevernalis’
变种
Varieties隶属函数值
Fuzzy mathematics总评分
Score综合排名
Comprehensive
ranking相对电导率
REC超氧化物歧化酶
SOD超氧化物酶
POD过氧化氢酶
CAT丙二醛
MDA可溶性蛋白
Soluble protein脯氨酸
Proline黄条雷竹
Ph. praecox f. Notata0.547 9 0.464 7 0.420 2 0.471 4 0.493 2 0.581 4 0.494 9 0.496 2 4 黄皮雷竹
Ph. praecox f. Chrysoderma0.532 1 0.604 5 0.454 8 0.487 5 0.473 2 0.471 0 0.492 5 0.502 2 3 弯秆雷竹
Ph. praecox cv. Linanensis0.565 5 0.636 1 0.449 6 0.476 8 0.558 2 0.642 8 0.691 0 0.574 3 1 花秆雷竹
Ph. praecox f. Viridisulcata0.475 1 0.466 2 0.378 2 0.378 6 0.295 7 0.265 8 0.369 6 0.375 6 5 细叶乌头雷竹
Ph. praecox ‘Xiyewutouleizhu’0.538 9 0.472 2 0.392 9 0.539 3 0.572 8 0.576 7 0.547 6 0.520 1 2 Table 5. The membership function and the comprehensive ranking on cold tolerance of 5 varieties