盐碱地是全球农业可持续发展的主要障碍之一,世界盐碱地面积占陆地面积的6%,我国各类盐碱地面积之和大约为1亿hm
2,主要分布在北方干旱半干旱地区
[1-2]。盐碱地分为以硫酸盐-氯化物为主的盐土、含有苏打的苏打盐土和pH>9、碱化度大于30%的碱土
[3-4]。水文地质、气候、地形和植被等原因使易溶于水的盐在土壤表层聚集,导致土壤含盐量高于正常水平,这是盐碱地形成的自然因素;人类不合理的灌溉与自然资源的不合理利用使得土壤盐碱化是次生盐碱地形成的重要原因
[5-6]。随着全球人口的增加,为了解决人多地少导致的粮食不足问题,合理利用盐碱化土地是非常必要的。但是盐碱混合胁迫引起的渗透胁迫、氧化胁迫、离子毒性和高pH值胁迫会破坏植物的细胞膜结构,使酶活性丧失,从而引起植物细胞中的离子失衡和代谢紊乱
[7-9]。因此,筛选耐盐碱作物对培育耐盐碱品种和盐碱地的改良利用具有重要的意义。
紫花苜蓿(
Medicago sativa L.)被誉为“牧草之王”,蛋白含量高,适口性好,是家畜喜食的优质豆科牧草,也可作为生态修复草,兼具经济效益与环境效益
[10]。紫花苜蓿的种植区域主要为干旱和半干旱地区,也是盐碱地分布区,虽然紫花苜蓿具有一定的耐盐碱能力,但日益加剧的土壤盐碱化严重影响了紫花苜蓿的生长发育及生理调控过程,降低了其产量和品质
[11]。在萌发期,盐碱胁迫会使紫花苜蓿种子吸水困难从而引起生理性干旱,细胞膜透性改变,造成离子堆积烧伤种子从而导致种子停止萌发
[12];在幼苗期,盐碱胁迫可显著降低紫花苜蓿的生长指标,造成渗透胁迫和氧化胁迫损伤,同时Na
+过度积累造成离子毒害
[13]。有研究用100 mmol/L硫酸钠盐碱溶液对4种紫花苜蓿种子进行胁迫,结果发现WL343HQ品种的抗盐碱能力最强
[14];高荣等
[15]用75 mmol/L Na
2CO
3-NaHCO
3(Na
2CO
3与NaHCO
3物质的量比=1∶2)对5种紫花苜蓿进行盐碱胁迫,结果发现Magna-salt品种苗期耐盐碱能力最强;崔禄
[16]采用
m(NaCl)∶
m(NaHCO
3)∶
m(Na
2SO
4)∶
m(Na
2CO
3)=1∶9∶9∶1的混合盐碱溶液对23种苜蓿进行盐碱胁迫处理,结果发现新牧一号的耐盐碱性最强,阿尔冈金的耐盐碱性最弱。尽管紫花苜蓿的耐盐碱特性已有不少研究报道,但目前关于紫花苜蓿种质资源对复合盐碱胁迫的响应差异和耐盐碱性综合评价还不够全面。因此,本研究结合宁夏地区盐碱土地类型,以175 mmol/L NaCl-Na
2SO
4模拟盐碱胁迫,对52份紫花苜蓿种质资源进行处理,测定其生长指标与生理指标,利用相关性分析、主成分分析、聚类分析和逐步回归分析等方法,综合评价了紫花苜蓿苗期的耐盐碱性,并筛选了评价的关键指标,旨在为后续种质资源的合理利用、抗性基因挖掘和耐盐碱品种的培育提供重要的理论支撑。
1 材料与方法
1.1 试验材料
52份紫花苜蓿种质资源均由宁夏大学林业与草业学院种子库保存。其中,三得利来自北京百绿种业有限公司;贝多叶、金星、青塔、冰峰、泰山,来自北京赤锐达生态科技有限公司;中苜2号和中苜3号,来自北京畜牧兽医研究所;先行者,来自北京正道种业有限公司;赤峰(地方种质资源)和扶风(地方种质资源)来自宁夏大学;渭南(地方种质资源)、罗马尼亚、拉达克、甘农6号、波兰、土库曼、甘农2号、东德、甘农8号、荷兰向阳、罗佑码、澳大利亚、礼县(地方种质资源)、雷纳伊、陕北子州、陕北(地方种质资源)和加拿大,来自甘肃农业大学;甘农5号,来自甘肃省农业科学院;阿尔冈金、雷达克之星、皇冠、金皇后、巨能6、巨能7、MF4020、耐盐之星、牧歌37CR和维多利亚,来自克劳沃生态科技有限公司;草原3号,来自内蒙古农业大学;康赛、骑士2和阿迪娜,来自宁夏草原站;中苜4号、赛沃5、牧歌耐盐和威龙,来自宁夏西贝农林牧生态科技有限公司;路特斯、秘鲁、印第安、伊犁(地方种质资源)和安斯塔,来自中国农科院草原研究所。
1.2 试验方法
基于前期预试验结果,本试验用175 mmo/L NaCl-Na2SO4(NaCl与Na2SO4物质的量比=9∶4),利用NaHCO3调节pH为8.3)模拟盐碱胁迫,对苗期紫花苜蓿进行胁迫处理。选取饱满的紫花苜蓿种子,首先用75%(体积分数)的酒精消毒15 s,蒸馏水冲洗干净后再用5%(体积分数)的次氯酸钠溶液消毒6 min,用蒸馏水冲洗直至无味后均匀地点播在培养皿中,于4 ℃冰箱中倒置春化2 d,然后放于人工气候培养箱(昼/夜温度25 ℃/20 ℃,相对湿度60%,光照/黑暗时间16 h/8 h,光照强度5 000 lx)中倒置生长2 d。将发芽一致的种子去除种皮后移栽到蛭石中,其间浇灌霍格兰营养液,待其生长至30 d时开始进行盐碱处理,为了防止应激,第1天用50 mmol/L的盐碱胁迫,第2天用100 mmol/L的盐碱胁迫,第3天进行175 mmol/L盐碱胁迫处理,共计处理30 d,每份种质资源选取生长一致的植株进行指标测定。
1.3 指标测定
形态及生长指标测定:利用卷尺测量株高(plant height,PH),用游标卡尺测量茎粗(stem thickness,ST),利用叶面积仪测量叶面积(leaf area,LA);使用万分之一分析天平称量52份紫花苜蓿地上部分的鲜质量(fresh weight,FW);测量后将地上部分分别放入烘箱中,105 ℃杀青30 min,65 ℃烘干至恒质量,然后测量干质量(dry weight,DW),计算每份种质资源的相对含水量(relative water content,RWC)。对每份种质资源的节数(number of nodes,NN)和分枝数(number of branches,NB)进行记录。
生理指标测定:称取0.2 g叶片,剪碎后放入20 mL蒸馏水中过夜,利用电导率仪测量初始电导率(C1),随后将其放入沸水浴锅中煮沸30 min测量煮沸后的电导率(C2),计算相对电导率(relative electrical conductivity,REC):REC=C1/C2×100%。将紫花苜蓿植株在黑暗处放置30 min,随后利用叶绿素荧光成像系统捕捉光合系统Ⅱ(PSⅡ)的最大量子产量(Fv/Fm)。利用手持叶绿素仪测量叶绿素含量(SPAD)。
1.4 数据分析
使用Excel 365软件对数据进行整理,使用SPSS 24进行分析。根据以下公式计算52份紫花苜蓿种质资源耐盐碱性的综合评价值(D),并结合聚类分析与逐步回归分析对紫花苜蓿种质资源的耐盐碱性进行归类,并筛选关键耐盐碱性指标。
(1)隶属函数值:
。
(2)各综合指标权重:
(3)综合评价值:
式中:u(Xj)为各指标隶属函数值,Xj 为指标相对值,Xmax和Xmin为所有参试种质资源某一指标的最大值和最小值,Wj 为第j个指标占所有指标的权重,Vj 表示各种质资源第j个指标的贡献率,D表示在盐胁迫下各种质资源耐盐能力的综合评价值。
2 结果与分析
2.1 盐碱胁迫下紫花苜蓿苗期各指标变异分析
由
表1可知,52种紫花苜蓿的株高为18.44~38.82 cm,平均为29.03 cm,变异系数为16.20%;茎粗为0.46~1.86 mm,平均为1.11 mm,变异系数为21.20%;节数为7.40~11.20节/株,平均为9.45节/株,变异系数为9.80%;分枝数为3.00~4.80枝/株,平均为3.81枝/株,变异系数为11.90%;叶面积为96.20~298.00 mm
2/片,平均为186.08 mm
2/片,变异系数为18.50%;鲜质量为1.30~3.28 g/株,平均为2.27 g/株,变异系数为18.70%;干质量为0.29~0.86 g/株,平均为0.60 g/株,变异系数为21.90%;相对含水量为62%~89%,平均为76%,变异系数为7.80%;相对电导率为41%~92%,平均为66%,变异系数为18.10%;
Fv/
Fm为0.61~0.82,平均为0.71,变异系数为7.30%;SPAD为36.38~72.22,平均为60.07,变异系数为12.70%。11个指标的变异系数在7.30%~21.90%,其中
Fv/
Fm的变异系数最小(7.30%),干质量的变异系数最大(21.90%),表明盐碱胁迫处理对紫花苜蓿苗期各个指标产生的影响差异较大。
2.2 盐碱胁迫下紫花苜蓿苗期各指标相关性分析
对盐碱胁迫下紫花苜蓿苗期的11个指标进行相关性分析,结果见
图1。
由
图1可知,株高与茎粗、节数、叶面积、鲜质量和干质量正相关,其中与茎粗、节数、叶面积、鲜质量的相关性达到极显著水平(
P<0.01),与干质量的相关性达显著水平(
P<0.05)。茎粗与干质量显著正相关(
P<0.05),与鲜质量极显著正相关(
P<0.01);节数与干质量显著正相关(
P<0.05),与叶面积和鲜质量极显著正相关(
P<0.01);分枝数与相对电导率显著负相关(
P<0.05);叶面积与鲜质量和干质量显著正相关(
P<0.05);鲜质量与干质量极显著正相关(
P<0.01);相对含水量与相对电导率显著负相关(
P<0.05);相对电导率与
Fv/
Fm极显著负相关(
P<0.01)。鲜质量与干质量的正相关系数最大(0.842),相对电导率与
Fv/
Fm的负相关系数最大(-0.508)。结果表明,盐碱胁迫下紫花苜蓿苗期各个指标之间存在一定的相关性,无法利用单独指标判断紫花苜蓿的耐盐碱性。
2.3 盐碱胁迫下紫花苜蓿苗期各指标主成分分析
以累计方差贡献率超过85%为标准,对盐碱胁迫下的紫花苜蓿苗期11个指标进行主成分(principal component,PC)分析,结果见
表2。由
表2可知,试验提取到PC1~PC7共7个主成分,其贡献率依次为28.889%,17.443%,12.463%,9.584%,8.667%,6.158%和5.343%,累计贡献率达到88.546%。由
表2可以看出,PC1中株高、茎粗、节数、叶面积、鲜质量和干质量的载荷系数较高,PC2中相对含水量和
Fv/
Fm的载荷系数较高,PC3中分枝数和SPAD的载荷系数较高,PC4中茎粗和SPAD的载荷系数较高,PC5中节数和叶面积的载荷系数较高,PC6中相对含水量的载荷系数较高,PC7中节数的载荷系数较高。
2.4 紫花苜蓿苗期耐盐碱性综合评价
根据主成分分析的结果,进一步计算52份紫花苜蓿种质资源PC1~PC7的隶属函数值,进而计算
D值,结果如
表3所示。由
表3可知,52份紫花苜蓿种质资源的
D值为0.332~0.845。
D值越大,表明紫花苜蓿种质资源的耐盐碱性越强,根据
D值大小对52份紫花苜蓿种质资源的耐盐碱性进行排名和分级,结果表明,耐盐碱性强的前5个种质资源分别为罗马尼亚(
D值为0.845)、贝多叶(
D值为0.784)、康赛(
D值为0.720)、中苜4号(
D值为0.708)和巨能6号(
D值为0.676);最不耐盐碱的种质资源为加拿大(
D值为0.332),其次是耐盐之星(
D值为0.409)、扶风(
D值为0.420)、安斯塔(
D值为0.430)和牧歌37CR(
D值为0.434)。
2.5 紫花苜蓿苗期耐盐碱性聚类分析
通过系统聚类分析方法,在欧氏距离为3.0处,可以将52份紫花苜蓿种质资源分为5个耐盐碱级别,结合耐盐碱性
D值排序结果,第1级别为强耐盐碱级别,包括2份种质资源,分别为罗马尼亚和贝多叶;第2级别为较强耐盐碱级别,包括7份种质资源,分别为康赛、中苜4号、巨能6、印第安、甘农6号、土库曼和拉达克;第3级别为中等耐盐碱级别,包括27份种质资源;第4级别为弱耐盐碱级别,包括15份种质资源;第5级别为不耐盐碱级别,共1份种质资源,为加拿大(
图2)。同时根据
D值的大小结合聚类分析结果,将52份紫花苜蓿种质资源划分成5个耐盐碱等级,
D值在0.784~0.845的种质资源为强耐盐碱品种,在0.656~0.720的为较强耐盐碱品种,在0.514~0.641的为中等耐盐碱品种,在0.409~0.502的为弱耐盐碱品种,
D≤0.332的为不耐盐碱品种。
2.6 逐步回归分析
以
D值作为因变量,将各个种质资源的指标平均值归一化处理后作为自变量,采用逐步回归分析进一步探究各指标与耐盐碱性的关系,根据
R2>0.9和
P值无限接近于0,建立最优回归方程:
D=0.246+0.222
X1+0.103
X2+0.111
X3+0.079
X4-0.092
X5+0.055
X6+0.118
X7+0.034
X8(
R2=0.990,
F=610.958,
P=0.000),式中
X1、
X2、
X3、
X4、
X5、
X6、
X7、
X8分别代表茎粗、干质量、
Fv/Fm、株高、相对电导率、分枝数、鲜质量和SPAD,表明这8个指标对紫花苜蓿种质资源的耐盐碱性影响显著。接着利用最优回归方程计算52份种质资源各自的回归
D值,结果见
表4。由
表4可以看出,回归方程计算的回归
D值在0.379~0.809,拟合误差在0.000~0.023,估计精度平均为0.987,证明逐步回归分析能准确地描述各指标与
D值间的关系,实现精准预测。结果表明,茎粗、干质量、
Fv/
Fm、株高、相对电导率、分枝数、鲜质量和SPAD等8个指标可作为评价紫花苜蓿种质资源耐盐碱性的关键指标。
3 讨 论
紫花苜蓿是多年生豆科牧草,其优势众多,不仅生态适应性强,而且营养价值高、适口性好,同时还具备固氮能力,是畜牧业优质饲料来源,目前已得到广泛栽培,但紫花苜蓿的产量和质量受到一系列因素的影响,包括各种非生物胁迫等外部环境因素
[17]。盐胁迫和碱胁迫在自然界中经常同时发生,盐碱胁迫会造成紫花苜蓿离子毒害、渗透胁迫以及抗氧化系统损伤,但目前大多数研究停留在盐或碱单一胁迫上,关于高盐与高pH伴随的盐碱混合胁迫对紫花苜蓿生长发育和品质影响的研究较少
[7]。紫花苜蓿是一种具有复杂的基因组和自交不亲和性的多年生杂交豆科植物,难以获得纯合自交系,对其进行耐盐碱性综合评价并筛选关键评价指标,可为耐盐碱苜蓿品种的选育及推广应用提供科学的理论依据和实践指导
[17-18]。
植物的外部形态和生长状况是其是否具有耐盐碱性最直观的表现,例如生物量以及地上部形态。混合盐碱胁迫下,植物株高、叶片数、茎长及地上部分干质量等均有所下降
[19],本研究也有类似的发现。对苗期各指标的变异分析发现,干质量的变异系数最大(21.90%),茎粗(21.20%)和鲜质量(18.70%)次之。同时,王宁等
[20]和张会丽等
[21]的研究也证实,盐胁迫会极显著抑制苜蓿的鲜质量、干质量及根长,这些发现均表明,植物生物量积累与形态建成过程对盐碱胁迫高度敏感。此外,盐碱胁迫还会减弱植物的光合作用、扰乱植物体内离子的动态平衡并造成细胞膜透性增加
[22]。本研究也对相关指标进行了分析,发现盐碱胁迫下紫花苜蓿的光合与水分代谢指标(
Fv/
Fm、SPAD、相对含水量)变异系数较小,总体来看这些生理功能在胁迫下比较稳定,这可能与紫花苜蓿固有的渗透调节能力有关
[23],但相对电导率的变异较大,表明盐碱胁迫会影响紫花苜蓿的细胞膜透性。相关性分析结果表明,鲜质量与干质量呈极显著正相关(
r=0.842),且二者均与株高、茎粗、叶面积显著关联,表明形态发育与生物量积累协同响应盐碱胁迫,这与李雪等
[24]的研究结果相符。值得注意的是,相对电导率与
Fv/Fm呈极显著负相关(
r=
-0.508),说明膜脂过氧化加剧直接损伤光系统Ⅱ反应中心,这为利用电导率快速鉴定耐盐碱性提供了生理依据
[25]。
通过主成分分析将11个指标转换为7个独立的综合指标变量,接着利用隶属函数法和聚类分析综合评价52份紫花苜蓿种质资源的耐盐碱性,得到2份强耐盐碱(罗马尼亚、康赛)和1份不耐盐碱(加拿大)的种质资源。这种多维指标降维方法有效地避免了单一指标的片面性,与姚兴洁等
[23]对61个苜蓿品种萌发期耐盐性评价的策略一致。吕艳贞等
[26]研究发现,在紫花苜蓿种子萌发期,根长、成苗率等指标在碱胁迫下对
D值贡献最大,而本研究发现苗期的干质量和茎粗等指标对
D值影响较大,说明不同生育期需采用差异化的评价体系
[21]。本研究将52份紫花苜蓿种质资源分为5个耐盐碱等级,其中强耐盐碱等级的种质资源仅占3.8%(2份),而中等及以下级别的种质资源占82.7%,说明多数引进种质资源的耐盐碱能力有限。张莹莹
[27]利用隶属函数法对130份紫花苜蓿种质资源进行综合评价时发现,加拿大为碱敏感性材料,本研究也发现加拿大为不耐盐碱种质材料,未来研究可将其作为解析盐碱敏感机制的对照材料
[25],可结合转录组学和代谢组学等多组学技术联合挖掘耐盐碱基因。本研究也对各个指标进行了逐步回归分析,发现茎粗、干质量、
Fv/
Fm、株高、相对电导率、分枝数、鲜质量和SPAD为评价紫花苜蓿苗期耐盐碱性的关键指标。已有研究发现,下调
SPL8基因可促进苜蓿的分枝数增多,显著提高苜蓿的生物量以及耐盐性
[28],表明分枝数和生物量与植物的耐盐性息息相关,这与本研究的结果一致;过表达
MsSPL12基因作为正调节因子可提高紫花苜蓿的耐盐性,同时使转基因紫花苜蓿的相对电导率降低
[29],同样表明相对电导率也可作为评价植物耐盐碱性的关键指标之一。除此之外,植物遭受盐碱胁迫时,光合作用效率会下降
[22],本研究也发现耐盐碱等级越低的紫花苜蓿,
Fv/Fm和SPAD的平均值也较低。综上,对紫花苜蓿种质资源进行耐盐碱胁迫评价时,不能根据单一的指标进行主观评价,需要根据各个指标的权重以及与耐盐碱性的相关程度进行综合评价,以构建合适且科学的评价体系。
4 结 论
本研究对52个紫花苜蓿种质资源苗期的耐盐碱性进行了综合评价,揭示了紫花苜蓿苗期耐盐碱性的显著差异。从52份紫花苜蓿种质资源中筛选出2份强耐盐碱种质资源(罗马尼亚和贝多叶),7份耐盐碱性较强的种质资源,27份中等耐盐碱的种质资源,15份耐盐碱性较弱的种质资源,1份不耐盐碱的种质资源(加拿大)。茎粗、干质量、Fv/Fm、株高、相对电导率、分枝数、鲜质量和SPAD等8个指标是评价紫花苜蓿苗期耐盐碱性的关键指标。
国家自然科学基金项目(U24A20426)