-
磷是植物体内核酸、蛋白质和磷脂等重要物质的必需组成成分,在核酸合成、光合与呼吸、细胞分裂增殖、生物膜的合成与稳定、酶的活化与失活、信号转导等过程发挥重要作用[1].由于磷素在土壤中易固定,移动性差,使得植物可吸收的有效磷很低[2],土壤速效磷的含量是维持喀斯特生态系统稳定的首要限制因子,对植被的机构和功能具有重要意义[3].植物在低磷胁迫下,其形态、生理、生化等方面都会发生显著变化,进而影响植物的生长与繁殖,最终影响植物的产量[4].大量研究表明,同一物种不同基因型活化土壤磷素能力不同,对磷素的吸收效率及有效利用能力不同[5-8].因此,挖掘利用磷高效植物的基因资源,是解决土壤磷素生物有效性低下问题的有效途径.
马尾松是原产于我国南方地区的针叶树种,该地区土壤胶体含有大量氢离子、铝离子,土质黏重,通透性差,土壤中有效磷很低,严重制约了马尾松生长及经济产量[9-10],因此挖掘和培育耐低磷新种质势在必行.前人通过研究筛选出许多与耐低磷相关的形态及生理指标,杨青等发现根构型的变化是马尾松高磷素吸收效率的重要机制[11].秦晓佳等发现MDA,保护酶活性及游离脯氨酸质量分数等指标是马尾松适应低磷胁迫的主要生理影响机制[10].
植物耐低磷特性是一个非常复杂的性状,单一形态或生理指标并不能全面地反映特定基因型的遗传本质,而多指标数据之间又存在一定的相关性,导致多变量之间出现信息重叠现象.前人有采用主成分分析法及灰色关联分析法进行常见作物新品种(系)的选育和评价[12-16],但迄今尚未见将这两种方法用于马尾松耐低磷种质筛选的报道.本文对7个种源地马尾松幼苗进行低磷胁迫,测定了15个形态及生理生化指标,采用主成分分析法和灰色关联分析法筛选耐低磷特性指标,并对不同种源地马尾松种质的耐低磷能力进行综合评价,以期为马尾松耐低磷遗传改良提供理论参考.
全文HTML
-
供试材料共7份,采自在全国各地生长表现较好的5个马尾松种源[17-18]:广西桐棉(GT)、福建龙岩(FL)、福建武平(FW)、江西崇义(JC)及湖南汝城(HR),以及2个中部红壤低磷地区种源:四川眉山(SM)及贵州孟关(GM).
-
选择种粒饱满无破损的种子经灭菌、催芽处理后于育苗盘中育苗,出苗35 d后移栽入塑料花盆(12 cm×11 cm)中,栽培基质为珍珠岩.每盆6株,重复3次.移栽5 d后开始浇营养液,营养液配方参考秦晓佳等[10]并稍加修改,在此基础上通过改变KH2PO4质量浓度将P浓度调节成2个梯度:1) 中磷(对照)处理(10 mg/L);2) 低磷处理(2.0 mg/L),相应改变KNO3和Ca(NO3)2·4H2O质量浓度以满足正常的N,K供应.处理在人工气候箱(上海新苗QHX-250BS-Ⅲ)中进行,光照10 h/d,光照度20 000 lx,温度18~22 ℃,湿度75%.
-
胁迫处理60 d后,小心取出整株幼苗,清洗干净,分别测定根系质量、地上部质量、根长(根颈部至主根生长点)、株高(根颈部至植株顶端生长点)等植株形态数据,并计算根冠比(根系鲜质量与地上部鲜质量的比值).
-
取各处理相同部位针叶簇及茎,测定生理生化指标.参照张志良等[19]的方法:叶绿素质量分数采用乙醇提取法测定;可溶性糖(SS)质量分数采用蒽酮比色法测定;丙二醛(MDA)质量分数采用硫代巴比妥酸(TBA)法测定;可溶性蛋白(SP)质量分数采用考马斯亮蓝G-250染料结合法测定;游离脯氨酸(Pro)质量分数采用酸性水合茚三酮显色法测定;超氧化物歧化酶(SOD)活性采用氯化硝基氮兰四唑(NBT)光还原法测定;过氧化物酶(POD)活性采用愈创木酚氧化比色法测定;过氧化氢酶(CAT)活性采用紫外吸收法测定;根系活力采用α-萘胺氧化法测定.
-
将各指标不同低磷处理的LPTC看作整体在SPSS 20.0中进行相关分析,并通过DPS 9.50处理软件选取欧氏距离采用最长距离法对其进行聚类分析.
-
式中,Cli为第i个综合指标,Clmin,Clmax为所有参试材料第i个综合指标的最小值和最大值.
-
式中,Wi表示第i个综合指标在所有综合指标中的重要程度即权重;Pi为各种源第i个综合指标的贡献率.
-
式中D值为各种源在低磷胁迫条件下综合耐低磷能力的大小(耐低磷性).
-
式中,“k”指代第k个指标,采用灰色关联分析法将低磷浓度下测定的10个生理生化指标(比较数列)与相应的D值(参考数列)进行关联分析[20-21],取分辨系数ρ=0.5计算出关联系数及关联度,并按关联度大小得出关联序,此过程在Microsoft Excel 2010中运算.
1.1. 试验材料
1.2. 试验设计
1.3. 测定项目与方法
1.3.1. 形态指标的测定
1.3.2. 生理生化指标的测定
1.4. 数据处理与统计分析
1.4.1. 耐低磷系数(low phosphorus tolerance coefficients,LPTC)
1.4.2. 隶属函数值
1.4.3. 指标权重
1.4.4. 综合值D值
1.4.5. 关联系数εi(k)及关联度γi
-
在低磷浓度胁迫下,不同种源马尾松的耐低磷系数差别较大(表 1):MDA(X9),SOD(X12),POD(X13)和CAT(X14)质量分数在低磷胁迫下均有不同程度的增加(LPTC>1);总叶绿素(X6)有所下降,耐低磷系数均小于1;其他指标的耐低磷系数因种源不同而表现较大差别.可见,各指标在马尾松耐低磷特性中所发挥的作用不尽相同,不同测定指标因种源差异所表现出的耐受性并不一致,因此,仅使用一个单一指标评价各种源的耐低磷特性是不够的.
对各形态指标和生理生化指标进行相关性分析,结果表明(表 2):低磷胁迫下,株高与游离Pro显著正相关,与SOD和POD显著负相关,与根系活力极显著正相关.主根长与根系活力极显著正相关;根系鲜质量与类胡萝卜素,SS和根系活力极显著正相关;地上部鲜质量与SOD,POD和根系活力等生化指标显著相关,根冠比与类胡萝卜素,SS和MDA极显著正相关,与CAT显著相关.不同形态指标间也存在相关性,根系鲜质量与主根长显著正相关,与株高、地上部鲜质量和根冠比极显著正相关;地上部鲜质量与主根长呈显著正相关,与株高极显著正相关(r=0.82).同样,各生理生化指标间也存在相关性,总叶绿素与SS,MDA,SOD和POD显著负相关,其中总叶绿素与POD相关性达到极显著水平;MDA与SS极显著相关;SP与游离Pro同其他指标间相关性最弱,只有游离Pro与株高显著相关.可见,各指标间存在不同程度的相关性,造成其所提供的信息存在叠加现象,影响了不同源马尾松耐低磷特性的客观性及准确性评价.
-
为避免因指标间的相关性而造成的信息重叠,客观反映不同种源在低磷胁迫下各指标的变化,本研究将5个形态性状指标(株高,主根长,根系鲜质量,地上部鲜质量,根冠比)和10个生理生化指标(总叶绿素,类胡萝卜素,SS,MDA,SP,游离Pro,SOD,POD,CAT,根系活力)在不同低磷环境下的耐低磷系数进行主成分分析.低磷胁迫下15个指标分成的15个主成分,从表 3可以看出,前4个综合指标累积贡献率达96.23%,其中第一个主成分就贡献了51.26%的信息,包括根质量,MDA,根系活力,地上部鲜质量,SS,株高,根冠比,主根长,游离Pro和类胡萝卜素等10个指标;第2个主成分由POD,SOD和CAT组成,可称为保护酶因子;第3个主成分只反应了总叶绿素质量分数的部分;第4个主成分反应了SP的信息.因此,根质量、保护酶因子、总叶绿素和可溶性蛋白(SP)可作为马尾松耐低磷种质的鉴定指标.
-
在多指标测定和主成分分析的基础上,本研究利用模糊数学中的模糊隶属函数法,对7个马尾松种源的耐低磷特性进行了综合评价,利用公式求出低磷条件下7个马尾松种源各主成分的隶属值U(i)和权重值Wi,并计算出代表各种源耐低磷能力的D值.由大到小为四川眉山(SM),福建武平(FW),福建龙岩(FL),贵州孟关(GM),江西崇义(JC),湖南汝城(HR),广西桐棉(GT). D值越大表示其忍耐低磷胁迫的能力越强,由此可知,四川眉山种源具有较强耐低磷特性,其次是福建的两个种源,武平略优于龙岩,而广西桐棉、江西崇义和湖南汝城的种源对低磷环境最为敏感,抗性最弱.
-
采用最长距离法对不同浓度下的所有D值进行系统聚类分析(图 1),在欧氏距离等于0.44的位置,7个马尾松种源按耐低磷能力的强弱被划分为3类,四川眉山种源具有最强的耐低磷特性自成一类;福建龙岩(FL)和福建武平(FW)耐低磷能力较接近聚在第二类,属于中等耐低磷类型;贵州孟关(GM)与江西崇义(JC),广西桐棉(GT)与湖南汝城(HR)这4个种源聚在第三类,属于较弱耐低磷类型.
-
本研究通过灰色关联分析法,对总叶绿素,类胡萝卜素,MDA,SS,SP,Pro,SOD,POD,CAT,根系活力等10个生理生化指标进行综合评价,旨在揭示比较序列与参考序列(D值)关联关系,从而评价各指标在耐低磷特性评定中的作用强弱.通过计算所得的关联系数如表 5,关联度由各指标关联系数的平均值来表示.根据灰色关联系统理论分析原则,关联度越大的指标,其与耐低磷特性(D值)的评价结果越接近.由表 6可以看出与马尾松耐低磷性关系最密切的是类胡萝卜素质量分数,关联度为0.566,其次是MDA,根系活力和游离脯氨酸质量分数,关联度分别为0.554,0.387和0.381.
2.1. 不同种源马尾松耐低磷系数和总相关性分析
2.2. 不同种源马尾松耐低磷特性的主成分分析
2.3. 不同种源马尾松耐低磷特性的模糊隶属函分析
2.4. 不同种源马尾松耐低磷特性的聚类分析
2.5. 灰色关联分析
-
马尾松是我国南方主要的乡土造林和工业原料树种[22],其分布范围广泛,分布区气候、土壤等生态条件存在较大差异,由于长期的适应性进化,产生了广泛的遗传变异,不同种源的马尾松对低磷环境的适应性存在差异,有必要通过一定的种源耐低磷试验,挖掘出对低磷胁迫下生长表现相对稳定的马尾松种源,同时探讨其适应低磷胁迫的应对机制.
马尾松应对低磷胁迫是一种受多基因控制的复杂生理反应,不同测定指标中所表现出的耐受性并不一致.因此,在植物耐低磷性研究中,往往采用多指标力求进行科学、客观的综合评判.但是评价指标过多时,指标间会存在一定的相关性,使得指标所提供的作物对逆境反映的信息发生交叉和重叠,且各指标在评价时因其本身之间权重差异直接导致评价结果发生偏差.主成分分析可通过降维的手段把多个复杂指标简化为少数几个综合指标,完全可以替代众多变量,并最大限度地反映原有多个变量的信息量[23],进一步求算出各个品种(系)抗逆性的综合评价值(D值).由于D值是一个无量纲的数,从而使得各个品种(系)的特性差异具有可比性[24].本研究采用主成分分析,避免重叠信息的干扰,准确分析马尾松耐低磷特性,从而挖掘出在不同程度的低磷胁迫环境下均表现优异的种源.结果表明,四川眉山的综合得分D值较大,在供试的7个马尾松种源中排名均位于前列,甚至优于前人普遍认可的磷高效型福建武平种源[10, 17].综合各成分的特征向量,地上部鲜质量、根系质量以及根系活力、保护酶系统和光合色素可以较好地反应马尾松的耐低磷特性.
灰色关联度分析是用来研究事物相互关联、相互作用等复杂因素的影响作用,使各种影响因素之间的“灰色”关系清晰化的一种统计分析方法[25],近年来成功地应用在作物新品种(系)的选育和评价中[26-27].灰色关联度分析通常把作物的多个农艺性状视为完整的一个灰色系统,通过关联度对作物各个农艺性状的重要性进行排序,也就是评价各指标在植物性状/特性评定中的作用强弱[28].本文根据灰色关联系统理论分析原则,关联度越大的指标,其与耐低磷特性的评价结果越接近.结果表明,类胡萝卜素质量分数、根系活力、游离Pro和保护酶活性等指标与马尾松耐低磷特性联系密切.
主成分分析法和加权关联分析都具有可靠的统计学基础,主成分分析可以将众多彼此相关的变量通过去除相关关系缩减为少数几个彼此无关的主成分;灰色关联分析能够从多个变量指标中找出主要决定指标,并全面、科学地评价不同品种的综合性能[28].本研究通过灰色关联分析获得的与马尾松耐低磷性密切相关的生理指标包括类胡萝卜素质量分数、根系活力、游离Pro和保护酶活性等,这些指标在主成分分析中同样被分在第一主成分中,证明两种方法在马尾松耐低磷种源筛选中同样可靠.
本文考查了5个农艺性状及10个生理生化性状,基本上能够全面描述不同种源马尾松低磷胁迫下的生长情况,但仍缺少其他相关的数据,如光合特性等.下一步研究将完善该综合评价体系,为其耐低磷良种资源的推广和应用提供参考依据.