巴青垂穗披碱草种子产量与农艺性状相关性分析

龙建廷，许赵佳，苏志刚，王明涛，苗彦军

（西藏农牧学院, 西藏 林芝 860000）

青藏高原被称为“世界屋脊”“地球第三极”和“亚洲水塔”，其在生态系统服务、生物多样性、水土保持等方面具有重要作用[1-2]。在全球气候变暖以及人类活动增加的背景下，青藏高原部分地区逐渐出现草地退化现象[3-4]。西藏自治区地处青藏高原，在其高寒的气候条件下却拥有全国面积最大的天然草地[5]。当前西藏栽培草地建植技术并不发达，栽培草地建植面积不足，所以即便坐拥我国最大的天然草地，西藏多数牧区还是出现牧草供不应求的情况[6]。垂穗披碱草(Elymus nutans)作为植被恢复的建群种，可以为土壤提供大量的凋落物，有助于提高土壤养分，一直被应用在青藏高原各地植被恢复中[7]。并且垂穗披碱草具有较高的种子繁殖能力，对土坡适应能力广，产量高，再生力强，可成为草地植物群落的建群种，一直以来都是草原地区种植栽培草地的首选草种之一[8]。

巴青垂穗披碱草(E.nutans‘Baqing’)是经过野生牧草选育出的西藏第1 个野生牧草驯化品种。该品种驯化选育始于 2004 年，经过近16 年驯化选育，于 2020 年12 月申报西藏自治区地方品种(品种登记号XZCS001)，经过近20 年的持续研究，育成的该品种适应西藏高寒环境，兼具产量高和品质好的特点，适于海拔3 000～4 800 m 区域种植。但是目前对林周县巴青垂穗披碱草种子扩繁基地种子产量的研究还处于初级阶段，针对巴青垂穗披碱草种子产量形成机制及农艺性状对其影响尚不明确，因此对巴青垂穗披碱草种子扩繁基地内种子产量与农艺性状进行相关性分析，以期为提高巴青垂穗披碱草种子产量提供理论参考。

1 材料与方法

1.1 供试材料

供试材料为西藏自治区审定品种巴青垂穗披碱草。试验材料在2021 年8 月于林周县巴青垂穗披碱草种子扩繁基地平均海拔4 200 m 播种，播种量为50 kg·hm-2。每年施用复合肥200 kg·hm-2，根据每日观测情况适量灌溉，播种当年未抽穗。试验指标测量于2022 年8 月。

1.2 试验方法与数据测量

于播种第2 年，选择一块肥力均匀、水分充足的试验地，根据垂穗披碱草根系及分蘖情况确定单株，随机选取100 株成熟巴青垂穗披碱草，为减少种子落粒导致误差在采收前提前进行套袋处理。用卷尺测定植株株高(地面至植株拉伸最高处)、旗叶长(穗下第1 片叶子长)、根长；用游标卡尺测定种子长度(不包括芒长)、芒长；用万分之一天平测定风干后单株种子重量。同时测定分蘖数(离地面10 cm的分支数)、单株穗数(有效分蘖数，即抽穗并结实的分枝)、平均穗粒数、平均拔节数(仅计算有效分蘖的平均拔节数)。计算影响垂穗披碱草种子产量指标包括单株穗数(A1)、平均穗粒数(A2)、单株产量(A3)，农艺性状指标包括种子长度(B1)、芒长(B2)、平均拔节数(B3)、株高(B4)、分蘖数(B5)、平均旗叶长(B6)、根长(B7)。其中平均拔节数、平均旗叶长为有效分蘖的平均数，种子长度与芒长为单株产出种子中随机选择20 粒种子测定的平均数。

1.3 数据处理

采用SPSS 对数据进行相关性分析。对典型变量U 和V 分别由3 个产量指标(A1～A3)和7 个农艺性状指标(B1～B7)通过线性组合构成。用逐步回归法分析建立巴青垂穗披碱草种子产量与上述农艺性状之间的关系方程。

2 结果与分析

2.1 各指标间的相关系数

在产量指标中，单株穗数(A1)与平均穗粒数(A2)呈极显著负相关关系(P＜ 0.01) (表1)，单株穗数(A1)与单株产量(A3)间呈极显著正相关关系(P＜ 0.01)。单株穗数(A1)与农艺性状相关性从大到小表现为分蘖数(B5) ＞ 平均旗叶长(B6) ＞ 株高(B4) ＞芒长(B2) ＞ 根长 (B7) ＞ 种子长度 (B1) ＞ 平均拔节数(B3)，其中单株穗数(A1)与平均拔节数(B3)呈极显著负相关关系(P＜ 0.01)，与分蘖数(B5)、平均旗叶长(B6)呈极显著正相关关系(P＜ 0.01)，与其他农艺性状无显著相关关系(P＞ 0.05)。平均穗粒数(A2)与农艺性状相关性从大到小表现为平均拔节数(B3) ＞株高(B4) ＞ 平均旗叶长(B6) ＞ 芒长(B2) ＞ 根长 (B7) ＞种子长度(B1) ＞ 分蘖数(B5)，其中平均穗粒数(A2)与平均拔节数(B3)和株高(B4)呈显著正相关关系(P＜ 0.05)，与种子长度(B1)呈显著负相关关系(P＜0.05)，与分蘖数(B5)呈极显著负相关关系(P＜ 0.01)。单株产量(A3)与农艺性状相关性从大到小表现为分蘖数(B5) ＞ 平均旗叶长(B6) ＞ 株高(B4) ＞ 芒长(B2) ＞ 种子长度(B1) ＞ 根长 (B7) ＞ 平均拔节数(B3)。其中单株产量(A3)与株高(B4)、分蘖数(B5)、平均旗叶长(B6)呈极显著正相关关系(P＜ 0.01)。

表1 巴青垂穗披碱草10 个变量间的相关系数Table 1 Correlation coefficient among 10 variables of Elymus nutans ‘Baqing’

2.2 典型相关分析

对巴青垂穗披碱草3 个种子产量指标为A 集合，7 个农艺性状指标为B 集合进行典型相关分析。得到3 对典型相关变量中，前两对典型变量相关系数达到0.01 极显著水平，具有统计学意义(表2)。前两对典型变量的标准化线性组合如下：U1=1.143A1- 0.025A2- 0.204A3,V1= 0.034B1+ 0.059B2-0.167B3- 0.800B4+ 0.950B5+ 0.067B6+ 0.051B7;U2= -0.149A1- 0.883A2- 0.395A3,V2= 0.485B1- 0.211B2+0.003B3- 0.705B4+ 0.340B5- 0.537B6+ 0.307B7。式中，A1～A3分别指单株穗数、平均穗粒数、单株产量，B1～B7分别指种子长度、芒长、平均拔节数、株高、分蘖数、平均旗叶长、根长。

表2 农艺性状与巴青垂穗披碱草种子产量间相关分析及典型冗余分析Table 2 Correlation analysis and typical redundancy analysis between agronomic traits and seed yield of Elymus nutans ‘Baqing’

根据典型变量线性组合公式可知，在第1 对典型变量线性组合构成中，U1和V1中分别以单株穗数和分蘖数的系数较大；在第2 对典型变量线性组合构成中，U2和V2中分别以平均穗粒数和株高的系数较大。表明第1 对典型变量主要以单株穗数和分蘖数决定，第2 对典型变量主要以平均穗粒数和株高决定。种子产量指标可以被第1、第2 典型变量(U1和U2)解释的比例分别为57.0%、35.9%；农艺性状指标可以被第1、第2 典型变量(V1和V2)解释的比例分别为17.5%和13.4% (表2)。农艺性状的典型相关结构分析表明，综合农艺性状指标变量V1和V2，分蘖数和株高结构系数相对较大，可见分蘖数与株高是决定垂穗披碱草种子产量的主要因子。

2.3 巴青垂穗披碱草与农艺性状的关系

2.3.1 单株产量与产量指标、农艺性状指标的关系

基于逐步回归分析法，建立巴青垂穗披碱草单株产量(A3)与单株穗数(A1)、平均穗粒数(A2)、种子长度(B1)、芒长(B2)、平均拔节数(B3)、株高(B4)、分蘖数(B5)、平均旗叶长(B6)、根长(B7)的关系方程如 下：YA3= - 0.980 4 + 0.265 1XA1+ 0.013 6XA2-0.318 1XB1+ 0.073 0XB2+ 0.024 6XB3+ 0.006 4XB4-0.031 3XB5+ 0.006 2XB6- 0.016 1XB7，(N=100,R2=0.909 9**)。

从上式可以看出，影响垂穗披碱草单株产量的主要指标是单株穗数(A1)、平均穗粒数(A2)，并且上述两个指标的相关系数达到了极显著水平，而其他农艺性状与单株产量(A3)的偏向相关关系不显著。这表明，单株产量会随着单株穗数和平均穗粒数的增加而增加。同时也表明，农艺性状指标对单株产量的影响不显著。

2.3.2 平均穗粒数与单株穗数、农艺性状指标的关系

基于逐步回归分析法，建立巴青垂穗披碱草平均穗粒数(A2)与单株穗数(A1)、种子长度(B1)、芒长(B2)、平均拔节数(B3)、株高(B4)、分蘖数(B5)、平均旗叶长(B6)、根长(B7)的关系方程，如下：YA2=32.166 9 - 10.277 9XA1+ 39.229 1XB1- 17.972 8XB2-13.811 1XB3+ 0.850 5XB4+ 0.580 9XB5+ 0.630 0XB6+0.285 2XB7，(N= 100,R2= 0.694 2**)。

从上式可以看出，影响垂穗披碱草平均穗粒数(A2)的主要指标是单株穗数(A1)、种子长度(B1)、平均拔节数(B3)、株高(B4)、平均旗叶长(B6)，垂穗披碱草平均穗粒数(A2)与芒长(B2)、分蘖数(B5)、根长(B7)的偏相关关系未达到显著性水平。这表明，影响平均穗粒数的主要农艺性状指标是种子长度、平均拔节数、株高、平均旗叶长，芒长、分蘖数、根长的影响则较小。

2.3.3 有效分蘖与农艺性状的关系

基于逐步回归分析法，建立巴青垂穗披碱草单株穗数(A1)、种子长度(B1)、芒长(B2)、平均拔节数(B3)、株高(B4)、分蘖数(B5)、平均旗叶长(B6)、根长(B7)的关系方程，如下：YA1= - 0.546 2 + 3.430 6XB1+2.270 0XB2- 1.254 9XB3+ 0.022 9XB4- 0.001 5XB5+0.127 3XB6+ 0.007 3XB7，（N=100,R2= 0.827 0**）。

从上式可以看出，影响垂穗披碱草单株穗数(A1)的主要指标是种子长度(B1)、芒长(B2)、平均拔节数(B3)、平均旗叶长(B6)，垂穗披碱草单株穗数(A1)与株高(B4)、分蘖数(B5)、根长(B7)的偏相关系数未达显著水平。这表明，单株穗数与农艺性状指标中种子长度、芒长、平均拔节数、平均旗叶长之间具有密切关系，与株高、分蘖数、根长的关系相对较弱。

3 讨论

牧草在完成生育期的过程中，会遇到各种因素，气温、雨水、光照等都对牧草种子产量影响较大[9-10]。农艺性状与种子产量的相关分析一直是育种专家学者关注的重点问题之一[11]。史关燕等[12]对谷子(Setaria italica)杂交种产量与农艺性状分析表明，单株穗数与单株粒重之间存在显著正相关关系。本研究中相关性分析结果表明，影响巴青垂穗披碱草单株产量的主要因素可能是单株穗数，这与史关燕等[12]的研究结果有相同之处。垂穗披碱草在播种两年间，种子产量与播种年限呈正相关关系，同时也有研究表明其他多年生牧草播种两年后，种子产量随着年限增加而降低，因此垂穗披碱草种子产量可能会随着种植年限的进一步增加，产量有所变化[13-14]。根据产量指标与农艺性状指标相关性结果分析可知，减少拔节数或提高分蘖数、旗叶长可增加单株穗数；减少种子长度和分蘖数或提高拔节数和株高可增加单株穗粒数；提高株高、分蘖数、旗叶长可增加单株产量。其中旗叶长、拔节数与部分研究[13-14]有所不同，可能是播种牧草品种、播种密度和试验区域不同导致的。

刘婷娜等[15]研究表明，垂穗披碱草种子产量与其构成因素的相关性不完全一致，单株干重相关系数最大。雷雄等[16]研究表明，喷施多效唑(PP333)可增加垂穗披碱草单株穗数，提高种子质量。彭珍等[17]研究表明，单株穗数和单穗种子数与垂穗披碱草种子产量呈极显著正相关关系(P＜ 0.01)。上述结果表明，垂穗披碱草种子产量的确是多个农艺性状共同作用的结果。当前关于垂穗披碱草种子产量与农艺性状相关的研究多围绕每生殖枝小穗数与每小穗种子数等指标，本研究针对巴青垂穗披碱草的产量指标与农艺性状指标分析发现，单株产量受单株穗数与平均穗粒数影响最大，这与前人[15-17]的研究基本一致。当前关于农艺性状对平均穗粒数与单株穗数影响的研究较少，其主要原因可能是相关学者的研究中平均穗粒数与单株穗数仅作为单株产量相关分析的指标，平均穗粒数与单株穗数对单株产量间接影响的研究较少。本研究结果表明平均穗粒数显著受种子长度、株高、平均旗叶长的影响，单株穗数显著受种子长度、芒长、平均旗叶长的影响。该结论为巴青垂穗披碱草种子产量的提高提供参考。

目前关于禾本科植物种子产量、农艺性状的研究主要应用简单相关[18-19]或主成分分析等方法[20-21]，多局限于两个因子间或集合内不同因子间的关系分析。相关研究发现[22-23]，种子产量与株高、分蘖数、旗叶宽等农艺性状有直接关系，其中株高影响种子产量的结论与本研究结果相近。

4 结论

垂穗披碱草是多年生禾本科牧草，对其分蘖数、株高、拔节数等性状的分析可有效解决种子扩繁中遇到的部分问题。本研究结果表明，最直接影响垂穗披碱草单株产量的因素是单株穗数与平均穗粒数；影响平均穗粒数的主要农艺性状指标则分别是种子长度、平均拔节数、株高、平均旗叶长；影响单株穗数的主要农艺性状指标分别是种子长度、芒长、平均拔节数、平均旗叶长。单株穗数和单株产量都与分蘖数、平均旗叶长呈极显著正相关(P＜0.01)，但平均穗粒数与分蘖数呈极显著负相关关系(P＜ 0.01)。因此或可通过增加旗叶长来增加单株穗数与单株产量，进而提高巴青垂穗披碱草种子产量。