0 引言
【研究意义】据统计当今世界13.19亿hm2耕地中约有43%的土壤缺磷,中国1.07亿hm2农田中约有50%的耕地土壤速效磷含量低于5 mg·kg-1[1],而紫花苜蓿生产所需有效磷一般为10—15 mg·kg-1 [2]。磷属于非可再生资源,磷矿资源预计还可以利用50–200年[3],依靠施用大量磷肥难以满足未来农业生产的需求,而且中国磷肥的当季利用率只有15%左右,总利用率也不超过25%[4],有效磷极易与金属离子及土壤颗粒形成AlPO4、FePO4、Fe3(PO4)2、Ca3(PO4)2等难溶磷[5-6]。所以提高紫花苜蓿（Medicago sativa）磷吸收利用效率,改善磷营养遗传特性是解决土壤有效磷缺乏的有效途径。通过挖掘紫花苜蓿自身遗传潜力,筛选耐低磷品种,可以为苜蓿磷高效吸收利用机制的研究提供理论依据。【前人研究进展】磷是植物生长所必需的大量元素之一,是构成核酸、磷脂、高能磷酸化合物等的重要组分,并参与能量代谢、糖分代谢、酶促反应、光合作用等重要代谢过程[7]。磷高效植物在磷素营养供应不足时仍能维持正常生长或受低磷影响较小,如磷高效大豆表现为磷吸收能力强,再分配效率高,低磷条件下保持较多的磷,因此表现为耐低磷[8]。耐低磷植物主要通过调整根形态结构、分泌有机酸和磷酸酶、结合丛植菌根真菌、诱导高亲和磷转运蛋白[9]等措施提高土壤磷的吸收效率[10],同时通过释放液泡中磷,改造膜减少对磷脂的依赖,转移衰老组织中的磷等方式提高磷利用效率以维持体内磷浓度的平衡[3]。根系作为植物营养吸收器官,在低磷胁迫下通过增加根冠比、促进侧根和根毛的生长扩大吸收磷面积[11]。其中最典型的是白羽扇豆（Lupinus albus）在低磷下形成排根,大大提高了磷吸收效率[12]。低磷下植物根系通过分泌磷酸酶催化磷酸单酯水释放磷[13],以及释放柠檬酸、苹果酸、草酸、琥珀酸等多种有机酸活化土壤中的Al、Fe、Ca等形成的难溶性磷酸盐[14-15]。目前已有多篇文章报道了大豆（Glycine max）耐低磷品种的筛选与评价[16-18],但关于紫花苜蓿耐低磷评价及机制研究较少。高艳等[19]指出黄花苜蓿（Medicago falcata）主要通过根系分泌柠檬酸活化根际难溶磷,而蒺藜苜蓿通过维持较大的根系来适应低磷环境。任立飞[20]发现低磷胁迫下锡盟黄花苜蓿的干重、鲜重、株高、地上部磷含量、净光合速率、气孔导度与蒸腾速率等下降幅度较小,而主根长与柠檬酸释放速率的增加幅度较大,综合评价得出锡盟黄花苜蓿的耐低磷性强于呼盟黄花苜蓿。【本研究切入点】解决土壤有效磷缺乏问题只能从增施磷肥或改变植物磷吸收效率两方面入手,但施入土壤的磷肥很容易被固定并造成磷资源浪费,所以挖掘植物遗传潜力,筛选耐低磷紫花苜蓿品种是解决土壤有效磷缺乏的有效途径。【拟解决的关键问题】本研究通过模拟低磷环境,分析20个紫花苜蓿品种幼苗期耐低磷综合表现,并基于各形态生理指标与耐低磷性的关系建立苜蓿耐低磷评价模型,以期筛选出耐低磷较强的紫花苜蓿品种以及与耐低磷密切相关的指标,为紫花苜蓿耐低磷评价及磷高效吸收利用机制研究提供科学依据。1 材料与方法
1.1 供试紫花苜蓿材料
参试材料为20个国内外紫花苜蓿品种,品种名见表1。Table 1
表1
表1试验材料
Table 1Alfalfa varieties
| 编号 NO. | 品种名称 Variety name | 编号 NO. | 品种名称 Variety name | |
|---|---|---|---|---|
| 1 2 3 4 5 6 7 8 9 10 | 中苜1号 Zhongmu NO. 1 中苜2号 Zhongmu NO. 2 公农2号 Gongnong NO. 2 甘农3号 Gannong NO. 3 龙牧801 Longmu 801 新牧1号 Xinmu NO. 1 陇东 Longdong 敖汉苜蓿 Aohan 巨能Ⅱ Magnum Ⅱ 牧歌37CR Amerigraze 37CR | 11 12 13 14 15 16 17 18 19 20 | 耐盐之星 Magnum Salt MF4020 驯鹿 Caribou 皇冠 Phabulous 三得利 Sanditi WL232HQ WL168HQ 牧歌 Amerigraze 康赛 Concept 骑士T Knight T |
新窗口打开
1.2 培养条件
挑选颗粒饱满颜色纯正的紫花苜蓿种子,用75%酒精灭菌后晾干,置于含有滤纸的培养皿发芽。保持湿润条件,待幼苗生长9 d时移至36 cm×26 cm×14 cm的含有营养液的方形培养盆,置于温室培养。每盆移植60棵,温室条件为昼/夜温度为25℃/20℃,湿度为60%,光照/黑暗为14 h/10 h,每4 d换一次营养液,营养液pH 6.0–6.3。营养液成分为1 000 μmol·L-1MgSO4,250 μmol·L-1 K2SO4,250 μmol·L-1 CaCl2,1 000 μmol·L-1 NH4NO3,2 500 μmol·L-1 KNO3,500 μmol·L-1KH2PO4,100 μmol·L-1 Fe-Na2-EDTA,30 μmol·L-1 H3BO3, 5 μmol·L-1 MnSO4,1 μmol·L-1 ZnSO4,1 μmol·L-1 CuSO4,0.7 μmol·L-1 Na2MoO4。紫花苜蓿幼苗在正常磷（500 μmol·L-1 KH2PO4）条件下生长15 d后,进行低磷（5 μmol·L-1 KH2PO4）胁迫,以常磷条件为对照,连续胁迫处理30 d后分别检测低磷与常磷条件下紫花苜蓿低磷相关性状指标,每个指标检测3个生物学重复,每个重复4个单株。1.3 性状测定
取低磷条件下培养30 d的紫花苜蓿幼苗,以常磷培养为对照,3组重复,每组取4株生长均一的植株,测量株高;取紫花苜蓿根及茎叶,100℃杀青30 min后60℃烘烤48 h,分别称量茎叶及根干重,并计算根冠比;用根系扫描仪扫描完整的根系并利用万深LA-S根系分析系统计算总根长及根表面积。酸性磷酸酶活性检测：称取0.1—0.3 g新鲜根,尽量取根尖部分,液氮条件下研磨至粉末状,溶于1.2 mL的0.1 mol·L-1提取缓冲液（1 mol·L-1冰醋酸溶液 8.82 mL和0.3 mol·L-1醋酸钠溶液273.3 mL,pH 4.0）,酸性磷酸酶活性以单位时间内单位鲜重质量的紫花苜蓿根内酸性磷酸酶水解对硝基苯磷酸二钠（ρ-NPP）生成对硝基苯酚（ρ-NP）的量（μmol·g-1·h-1）表示[17]。
植株地上部分磷含量的测定：植株茎叶在100℃杀青30 min后于60℃烘烤48 h至恒重,采用钼锑抗比色法测定植株磷含量[21],并计算磷利用率[17],植株磷利用率=植株干重/植株磷吸收量。
1.4 数据统计分析
利用SAS 9.2进行相关性分析、主成分分析、聚类分析、逐步回归分析,采用隶属函数对20个品种的紫花苜蓿进行耐低磷综合评价。指标计算参考相关文献[22-23]。分别计算茎叶干重、株高、根干重、根冠比、总根长、根表面积、酸性磷酸酶活性、植株全磷含量及植株磷利用率的耐低磷系数[17]。
耐低磷系数= \(\frac{低磷处理下性状测得值}{常磷处理下性状测得值}\) （1）
不同品种紫花苜蓿综合指标的隶属函数值计算：
μ (Xj) = \(\frac{X_j-X_{min}}{X_{max}-X_{min}}\) j=1, 2, 3,…,n （2）
式中,µ (Xj)表示第j个综合指标的隶属函数值,Xj表示第j个综合指标值;Xmax表示第j个综合指标的最大值,Xmin表示第j个综合指标的最小值。
各综合指标权重的计算：
Wj = \(\frac{P_j}{\sum
m_{j=1}P_j}\) j=1, 2, 3, …, n （3）
式中,Wj表示第j个综合指标在所有综合指标中的重要程度及权重;Pj为各品种第j个综合指标的贡献率。
耐低磷综合评价值的计算：
D= \(\sum
n_{j=1}[\mu (X_j)\times X_j]\) j=1, 2, 3, …, n （4）
式中,D表示在低磷胁迫条件下各品种耐低磷能力的综合评价值。
2 结果
2.1 低磷下各性状差异分析
紫花苜蓿幼苗在常磷与低磷胁迫下各性状都存在差异。为了消除品种间固有的差异,采用了耐低磷系数（即相对值）计算方法,以便准确地反映不同品种的耐低磷性。如表2所示,紫花苜蓿茎叶干重的耐低磷系数变化范围在0.5930—1.0756,有17个品种小于1,株高的耐低磷系数变化范围为0.6454— 1.1828,其中18个品种小于1,综合表明,低磷导致茎叶干物质量及株高下降;根干重和根冠比耐低磷系数分别为0.6429—1.6120和1.0805—1.5648,表明低磷导致根系比例增大;总根长和根表面积耐低磷系数变化范围分别为0.6455—1.3110和0.6480— 1.3129,其中敖汉、耐盐之星、皇冠、驯鹿、甘农3号总根长及根表面积的耐低磷系数均大于1,有报道显示紫花苜蓿通过增大根系来适应低磷环境[19];全磷含量的耐低磷系均小于1,在0.2442—0.4471,说明低磷下全磷含量均低于常磷下的紫花苜蓿,但磷利用率的耐低磷系数均大于1,变化范围为2.2462—4.1289;酸性磷酸酶活性的耐低磷系数变化范围为0.5282—1.5567。Table 2
表2
表2低磷胁迫下不同品种各性状的耐低磷系数
Table 2Low-P tolerant coefficient of each trait of different alfalfa varieties under low-Pi stress
| 编号 NO. | 品种 Variety | 耐低磷系数 Coefficients of resistance to deficient-Pi | ||||||||
|---|---|---|---|---|---|---|---|---|---|---|
| 茎叶干重 SLW | 株高 PH | 根干重 RW | 根冠比 RS | 总根长 TRL | 根表面积 RSA | 全磷含量 TP | 磷利用率 PUE | 酸性磷酸酶活性ACPA | ||
| 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 | 中苜1号 Zhongmu NO. 1 中苜2号 Zhongmu NO. 2 公农2号 Gongnong NO. 2 甘农3号 Gannong NO. 3 龙牧801 Longmu 801 新牧1号 Xinmu NO. 1 陇东 Longdong 敖汉苜蓿 Aohan 巨能Ⅱ Magnum Ⅱ 牧歌37CR Amerigraze 37CR 耐盐之星 Magnum Salt MF4020 驯鹿 Caribou 皇冠 Phabulous 三得利 Sanditi WL232HQ WL168HQ 牧歌 Amerigraze 康赛 Concept 骑士T Knight T | 0.9340 0.6992 0.7672 0.6992 0.9258 1.0402 0.7303 1.0756 0.5930 1.0111 0.9759 0.9205 0.8733 0.9952 0.9665 0.9149 0.6619 0.7101 0.6296 0.7637 | 0.9632 0.7253 0.8526 0.7737 0.7862 1.0752 0.9295 0.9779 0.6454 0.9227 0.9834 0.8788 0.9847 0.9438 0.8917 0.7204 0.9182 0.7439 0.7379 1.1828 | 1.0607 0.9313 0.8830 0.9161 1.1436 1.2700 1.0429 1.6120 0.6429 0.9413 1.3832 1.1612 1.2046 1.4517 1.0932 0.9650 0.8944 0.9160 0.8928 1.1874 | 1.1322 1.3044 1.1846 1.3464 1.2134 1.2317 1.4006 1.4996 1.2205 1.0805 1.4287 1.2080 1.3850 1.5106 1.1247 1.1159 1.4340 1.3683 1.5074 1.5648 | 0.8172 0.9591 0.9474 1.0306 0.9539 0.8886 1.0645 1.3096 0.7736 1.0846 1.1365 0.8399 1.2044 1.3110 0.8620 0.6455 1.1895 0.9791 0.7299 1.0815 | 0.7622 0.8583 1.0452 1.1305 1.0699 0.9416 0.9140 1.2848 0.6480 0.9162 1.1656 1.0154 1.0359 1.3129 1.0318 0.8685 0.9947 1.1132 0.6578 1.0114 | 0.2634 0.2442 0.2795 0.3846 0.3895 0.4181 0.4471 0.3705 0.3582 0.3248 0.3624 0.2506 0.4083 0.3890 0.2794 0.3824 0.4007 0.3762 0.3796 0.4379 | 3.8724 4.1289 3.6083 2.6772 2.7017 2.3950 2.2462 2.7062 2.8092 3.1022 2.7769 4.1172 2.4593 2.8168 3.6220 2.6944 2.5013 2.6671 2.6771 2.3013 | 1.1416 0.7375 0.5282 1.1160 1.3806 1.4325 0.9004 0.9836 1.0737 1.3131 1.2036 1.4194 1.0053 0.8938 1.0980 1.5567 0.9389 0.8933 0.9307 0.9401 |
新窗口打开
2.2 性状相关性分析
通过分析不同性状之间的相关性可以揭示性状间是否存在依存关系及相关关系的方向与强度[24]。结果显示部分指标间存在一定相关性,即不同性状间可能存在信息重叠,单靠单项性状难以对表现复杂的耐低磷性进行综合评价。如表3,茎叶干重与株高、根干重、根表面积、酸性磷酸酶活性呈正相关（P<0.05）;根冠比与总根长、全磷含量呈极显著正相关（P<0.01）,与酸性磷酸酶活性呈显著负相关（P<0.05）;总根长与根冠比、根干重、株高、根表面积都呈极显著正相关（P<0.01）。全磷含量与磷利用率成极显著负相关（P<0.01）,相关系数为-0.9536,全磷含量与根冠比的相关系数为0.5529,呈显著正相关（P<0.05）;酸性磷酸酶活性与茎叶干重呈显著正相关（P<0.05）,与根冠比呈显著负相关（P<0.05）。Table 3
表3
表3低磷胁迫下各性状的耐低磷系数相关性
Table 3Correlation of low-Pi tolerant of each character of alfalfa in low-Pi stress
| 性状指标 Index | 耐低磷系数（低磷下性状测得值/常磷下性状测得值） Coefficients of resistance to deficient-Pi (index value under deficient-Pi/index value under sufficient-Pi) | |||||||
|---|---|---|---|---|---|---|---|---|
| 茎叶干重 SLW | 株高 PH | 根干重 RW | 根冠比 RS | 总根长 TRL | 根表面积 RSA | 全磷含量 TP | 磷利用率 PUE | |
| 株高 PH 根干重 RW 根冠比 RS 总根长 TRL 根表面积 RSA 全磷含量 TP 磷利用率 PUE 酸性磷酸酶活性 ACPA | 0.5088* 0.7676** -0.2448 0.2554 0.4852* -0.1070 0.1082 0.5042* | 0.6564** 0.3150 0.5208* 0.3783 0.2862 -0.2422 0.0292 | 0.3951 0.5903** 0.7062** 0.1898 -0.1481 0.1499 | 0.5897** 0.3581 0.5656** -0.5277 -0.4577* | 0.7439** 0.3366 -0.3243 -0.3541 | 0.1717 -0.1552 -0.0990 | -0.9848** 0.0797 | -0.0805 |
新窗口打开
2.3 主成分分析
主成分分析可在损失较少信息的基础上将大量指标转化为有效少量指标,从而达到简化指标及浓缩数据的目的[22]。由表4—5可知,前4个主成分的贡献率分别为42.17%、25.78%、17.27%、6.82%,其中第1主成分的特征值是3.7950,贡献率为42.17%,对应较大的特征向量是根干重和总根长,反映了根系因子;第2主成分特征值是2.3202,贡献率为25.78%,载荷较高的是茎叶干重;第3主成分特征值为1.5542,贡献率为17.27%,载荷较高的是全磷含量和酸性磷酸酶活性;第4主成分特征值为0.6134,贡献率为6.82%,载荷较高的是株高;前4个主成分累计贡献率达92.04%,其余贡献率可忽略不计,这样可以将9个单项指标转化为4个相互独立的综合指标,并代表了原始指标的大部分信息。Table 4
表4
表4各综合指标特征根及贡献率
Table 4The eigenvalue and contribution of each comprehensive index
| 主成分 Principal component | 特征值 Eigenvalue | 贡献率 Contributive ratio (%) | 累计贡献率 Cumulative contributive ratio (%) |
|---|---|---|---|
| 1 2 3 4 | 3.7950 2.3202 1.5547 0.6134 | 42.17 25.78 17.27 6.82 | 42.17 67.95 85.22 92.04 |
新窗口打开
Table 5
表5
表5各因子载荷矩阵
Table 5Matrix of factor loading
| 耐低磷指标 Low-P tolerant index | 主成分 Principal component | |||
|---|---|---|---|---|
| 1 | 2 | 3 | 4 | |
| 茎叶干重 SLW 株高 PH 根干重 RW 根冠比 RS 总根长 TRL 根表面积 RSA 全磷含量TP 磷利用率PUE 酸性磷酸酶活性 ACPA | 0.2301 0.3685 0.4219 0.3442 0.4305 0.3911 0.3034 -0.2878 -0.0398 | 0.5548 0.1648 0.3063 -0.374 -0.0408 0.1713 -0.379 0.3827 0.3319 | 0.1564 0.0388 -0.0241 -0.1779 -0.2931 -0.2283 0.439 -0.448 0.6405 | -0.0216 0.7733 0.0299 0.0986 -0.1466 -0.5708 -0.0569 0.1036 -0.1723 |
新窗口打开
2.4 品种综合评价
2.4.1 隶属函数分析 根据主成分分析得到的4个相互独立的综合指标,利用公式（2）计算不同品种各综合指标的隶属函数值μ(Xj)。从表6可以看出,在同一综合指标CI(1)中,敖汉的μ(X1)最大,为1.0000,巨能Ⅱ的μ(X1)最小,为0.0000,表明敖汉在CI(1)综合指标上耐低磷性最强,巨能II耐低磷性最差。根据各综合指标贡献率,利用公式（3）计算其权重Wj,计算得到4个综合指标的权重分别为0.4582、0.2801、0.1876、0.0741（表6）。进一步利用公式（4）计算耐低磷综合评价值D值,并根据D值大小对其耐低磷性强弱排序,其中耐低磷较强的有敖汉、新牧1号、耐盐之星等;耐低磷较弱的品种有巨能Ⅱ、中苜2号、康赛等（表6）。2.4.2 耐低磷品种聚类分析 为了客观地将不同耐低磷紫花苜蓿分类,对4个新的综合指标的耐低磷评价值D采用最大距离法进行了聚类分析[22,25],根据最大距离值大于1.2050可将20个紫花苜蓿品种分为3类,第一类为耐低磷较强的品种,耐低磷综合评价值为0.6534—0.7352,包括敖汉、新牧1号、耐盐之星、皇冠;第二类为中度耐低磷品种,耐低磷综合评价值为0.4404—0.6534,包括骑士T、驯鹿、牧歌37CR、龙牧801、MF4020、三得利、中苜1号、WL232HQ、陇东;第三类为耐低磷较弱的品种,耐低磷综合评价值为0.1733—0.4404,包括WL168HQ、甘农3号、牧歌、公农2号、康赛、中苜2号、巨能Ⅱ（图1）。
Table 6
表6
表6各品种的综合性状指标值、权重、μ(X)及综合评价值（D）
Table 6The value of each comprehensive index (CI), index weight, μ(X), and comprehensive evaluation value (D)
| 品种 Variety | CI (1) | CI (2) | CI (3) | CI (4) | μ (X1) | μ (X2) | μ (X3) | μ (X4) | 综合评价值 Comprehensive assessment value (D) | 排名 Order |
|---|---|---|---|---|---|---|---|---|---|---|
| 中苜1号 Zhongmu NO. 1 中苜2号 Zhongmu NO. 2 公农2号 Gongnong NO. 2 甘农3号 Gannong NO. 3 龙牧801 Longmu 801 新牧1号 Xinmu NO. 1 陇东 Longdong 敖汉苜蓿 Aohan 巨能Ⅱ Magnum Ⅱ 牧歌37CR Amerigraze 37CR 耐盐之星 Magnum Salt MF4020 驯鹿 Caribou 皇冠 Phabulous 三得利 Sanditi WL232HQ WL168HQ 牧歌 Amerigraze 康赛 Concept 骑士T Knight T 权重 Weight | -0.9800 -1.2228 -0.7423 -0.0425 0.0378 0.5970 0.4616 1.8544 -1.6613 -0.3245 1.0871 -0.7275 0.9327 1.6637 -0.5684 -0.9725 0.3480 -0.1319 -0.8729 1.2644 0.4582 | 1.2688 -0.0054 -0.0068 -0.7395 0.4256 0.7448 -1.2752 0.7125 -1.2795 1.0353 0.5898 1.6835 -0.3240 0.2455 1.2878 0.4402 -1.2750 -0.9377 -1.6145 -0.9765 0.2801 | -0.2398 -2.0016 -1.8828 0.0390 1.0440 1.8473 0.4010 -0.6802 0.5721 0.5372 0.0217 -0.3017 0.0934 -0.9268 -0.4058 2.0459 -0.2597 -0.3868 0.2906 0.1929 0.1876 | 1.7953 0.0988 0.1817 -1.5199 -1.4439 1.1452 0.5614 -0.6206 -0.2537 0.1243 -0.1594 0.0212 0.3147 -0.9114 0.0710 -1.0433 0.1135 -1.3997 0.7735 2.1513 0.0741 | 0.1938 0.1247 0.2614 0.4605 0.4833 0.6424 0.6038 1.0000 0.0000 0.3802 0.7818 0.2656 0.7378 0.9458 0.3109 0.1959 0.5715 0.4350 0.2243 0.8322 | 0.8743 0.4879 0.4875 0.2653 0.6186 0.7154 0.1029 0.7056 0.1016 0.8035 0.6684 1.0000 0.3913 0.5640 0.8800 0.6230 0.1030 0.2052 0.0000 0.1935 | 0.4353 0.0000 0.0294 0.5042 0.7525 0.9509 0.5936 0.3265 0.6359 0.6272 0.4999 0.4200 0.5176 0.2655 0.3943 1.0000 0.4304 0.3990 0.5663 0.5422 | 0.9030 0.4409 0.4635 0.0000 0.0207 0.7260 0.5669 0.2450 0.3449 0.4479 0.3706 0.4198 0.4997 0.1658 0.4333 0.1298 0.4449 0.0328 0.6247 1.0000 | 0.4823 0.2265 0.2962 0.3799 0.5374 0.7269 0.4589 0.7352 0.1733 0.5501 0.6666 0.5117 0.5818 0.6534 0.4950 0.4615 0.4044 0.3341 0.2553 0.6113 | 11 19 17 15 8 2 13 1 20 7 3 9 6 4 10 12 14 16 18 5 |
新窗口打开
显示原图|下载原图ZIP|生成PPT图120个紫花苜蓿品种的耐低磷性聚类图
-->Fig. 1The cluster of low-Pi tolerance in 20 alfalfa varieties
-->
2.5 回归模型的建立
为了检验主成分分析结果,同时建立紫花苜蓿耐低磷评价数学模型[22],本试验以耐低磷综合评价值（D）为因变量,以各单项指标的耐低磷系数作自变量进行逐步回归分析[25],得到回归方程：D=-0.7997+ 0.3856 SLW+0.2025 PH+0.3789 RW+0.1051 TRL+ 0.4188 TP+0.1347 ACPA（P<0.01）,方程的决定系数R2=0.9982,F=1211.37。由方程可知,在9个指标中有6个指标对紫花苜蓿耐低磷性有显著影响,分别是茎叶干重、根干重、株高、总根长、全磷含量、酸性磷酸酶活性,与主成分分析结果相吻合。对20个紫花苜蓿品种的耐低磷性的综合评价值与预测值进行相关性分析发现,相关系数为r=0.92674,达到极显著水平（P<0.01）,说明此方程对紫花苜蓿耐低磷性预测效果较好。3 讨论
土壤对磷具有强烈吸附作用,通常有效磷含量为0.1—10.0 μmol·L-1 [26],目前已报道了低磷（有效磷浓度≤10 μmol·L-1）对紫花苜蓿形态生理方面的影响[15,19-20]。在室内试验营养液培养环境中,一般将500 μmol·L-1 KH2PO4作为苜蓿[19-20]和大豆[27-28]正常磷水平,所以本试验以常磷500 μmol·L-1 KH2PO4为对照,在低磷（5 μmol·L-1 KH2PO4）环境下对国内外20个紫花苜蓿品种进行了耐低磷评价与分析。植物耐低磷表现是一个复杂的过程,同时受遗传背景和环境因素影响,运用多个指标综合评价才能有效反映出不同品种的耐低磷性。本试验为了消除品种间固有的差异,本试验引用了耐低磷系数进行分析。低磷胁迫下植物净光合速率、气孔导度与蒸腾速率均下降,从而引起地上生物量减少[29-30],磷含量降低[20],根冠比[31]、根系表面积、总根长、磷利用率增大,根尖酸性磷酸酶活性增强[17,20,32]等现象。低磷胁迫下,叶片光合作用下降[33],叶肉细胞变小[34],叶面积、叶片生长速率和叶片数量显著减少,从而降低了同化面积并严重限制光合作用的总碳固定量[35]。光合作用的下降主要因为有机磷循环以及RuBP（核酮糖二磷酸羧化酶）更新速率的降低,进而导致植物叶片RuBP浓度、RuBP羧化酶活性、RuBP羧化酶结合位点都有不同程度的降低[36-37]。为了适应低磷环境,植物优先将光合产物分配到地下部分,并抑制地上部分生长[11],表现为根冠比增加[32,38-39]。本试验中,低磷胁迫下大部分紫花苜蓿品种的茎叶干重、株高有所下降,所有品种表现为全磷含量下降,根冠比和磷利用率增加。
低磷下植物根系分泌酸性磷酸酶也是植物适应低磷胁迫的重要机理[32],酸性磷酸酶能酶促土壤有机磷化合物的酯磷键发生裂解,释放相应的醇和无机磷,从而提高土壤磷素的有效性[40]。高彬等[41]指出大豆缺磷时根系磷含量降低,但根系组织酸性磷酸酶活性是全磷培养植株的1.3倍。刘渊[17]发现大豆在2 μmol·L-1和1 000 μmol·L-1处理下22个品种根尖酸性磷酸酶的耐低磷系数为1.9229—35.2011。在本试验中,低磷胁迫下部分品种根组织的酸性磷酸酶活性有所增强。
关于植物耐低磷的形态及生理指标有很多,本试验初步选取了9个指标用于紫花苜蓿耐低磷评价。通过主成分分析可将9个原始的单项指标转换成4个新的相对独立的综合指标,进一步对4个新的综合指标利用隶属函数分析法得到耐低磷综合评价值,这样比较客观的展示了不同品种的耐低磷性,避免了利用单个性状指标描述耐低磷的局限性。聚类分析可以根据研究对象之间亲疏关系的统计量进行分类,从而排除人为确定标准带来的主观影响,该方法已普遍应用于植物抗逆材料的筛选研究[22,24-25],通过聚类可将20个紫花苜蓿品种客观地分为耐低磷品种、中度耐低磷品种、耐低磷较弱的品种。逐步回归得到6个显著影响耐低磷的单项指标,与主成分分析结果一致,同时得到紫花苜蓿耐低磷性评价数学模型,可以预测不同紫花苜蓿品种耐低磷性的强弱,使抗逆性的鉴定与利用研究更有预见性[23]。
4 结论
本研究对20个紫花苜蓿品种进行了耐低磷评价及筛选,根据耐低磷能力强弱可分为3类,其中耐低磷品种包括敖汉、新牧1号、耐盐之星、皇冠,中度耐低磷品种包括骑士T、驯鹿、牧歌37CR、龙牧801、MF4020、三得利、中苜1号、WL232HQ、陇东,耐低磷较弱的品种包括WL168HQ、甘农3号、牧歌、公农2号、康赛、中苜2号、巨能Ⅱ。低磷胁迫下,茎叶干重、株高、根干重、总根长、全磷含量、酸性磷酸酶活性可作为筛选耐低磷紫花苜蓿的参考指标。本试验为进一步研究紫花苜蓿磷高效吸收利用机制提供了理论依据。The authors have declared that no competing interests exist.
