摘要: 分析影响高油大果花生品种郑花6号主要农艺性状对产量性状的直接、间接作用及其决定系数，为优化品种布局，丰富花生优质高产栽培理论提供参考。选用连续两年郑花6号参试的22个花生品种(系)为试验材料，运用基因环境互作双标图综合评价参试品种的丰产、稳产及适应性并评价18个试点的代表性和鉴别力，并利用SPSS对高油大果花生新品种郑花6号的产量及农艺性状进行相关分析和通径分析。双标图结果显示，两年产量性状的总变异值(PC1 + PC2)变异幅度为45.73%~67.28%，初步明确了密云、青岛、菏泽和濮阳、潍坊、安阳试验点对产量性状具有较强辨别力和较高的代表性，郑花6号丰产稳产性均较好，适合在河南、辽宁、山东、河北、江苏、安徽等省份大面积推广。变异分析表明单株结果数、主茎高、侧枝长的变异系数较大，出米率的变异系数最小，9个性状的变异系数在3.89~21.72%；相关分析结果显示，百果重、总分枝数、百仁重、结果枝数、出仁率、单株果数与荚果产量呈正相关，但无显著相关性(p > 0.05)，相关程度表现为百果重(r = 0.407) > 百仁重(r = 0.283) > 结果枝数(r = 0.244) > 出米率(r = 0.227) > 总分枝数(r = 0.189) > 单株果数(r = 0.175)；主茎高(r = −0.104)、侧枝长(r = −0.115)与荚果产量呈负相关，但无显著相关性(p > 0.05)。通径分析结果表明，主茎高(直接通径系数 = 1.163)，对产量的直接影响最大，其次为百果重(直接通径系数 = 0.945)、单株果数(直接通径系数 = 0.339)和出米率(直接通径系数 = 0.305)。株高、百果重和单株果数是影响郑花6号产量的主要参考性状，可作为郑花6号高产特征的指标参考，并协同出米率、结果枝数进行选择。在郑花6号栽培中应配合适宜播期，择适宜播种密度，建立合理群体结构，加强后期管理，控制株高，以实现优质高效栽培目标。
Abstract: Analyze the direct and indirect effects and determination coefficients of the main agronomic traits on yield traits of the high oil and large pod peanut variety Zhenghua 6, in order to optimize variety layout and enrich the theory of high-quality and high-yield cultivation of peanuts. 22 peanut varie-ties (lines) of Zhenghua 6 tested for two consecutive years were selected as experimental materials, and the high yield, stable yield and adaptability of the tested varieties were comprehensively eval-uated by using the gene-environment interaction biplot, and the representativeness and discrimi-nation of 18 pilot varieties were evaluated. The yield and agronomic traits of the new high-oil large-pod peanut variety Zhenghua 6 were analyzed by using SPSS. The biplot results showed that the variation range of PC1 + PC2 was 45.73%~67.28%. It was preliminarily confirmed that Miyun, Qingdao, Heze, Puyang, Weifang and Anyang test sites had strong differentiation and high repre-sentativeness for yield traits, and Zhenghua 6 had good yield stability, suitable for large-scale pro-motion in Henan, Liaoning, Shandong, Hebei, Jiangsu, Anhui and other provinces. The results of var-iation analysis showed that the variation coefficients of pod number per plant, main stem height and lateral branch length were large, and the variation coefficients of rice yield were the smallest, ranging from 3.89% to 21.72%. Correlation analysis showed that pod weight, total branch number, kernel weight, pod branch number, pod number per plant were positively correlated with pod yield, but there was no significant correlation (p > 0.05). The correlation was as follows: 100 pod weight (r = 0.407) > 100 kernel weight (r = 0.283) > number of pod branches (r = 0.244) > rice yield (r = 0.227) > total branch number (r = 0.189) > pod number per plant (r = 0.175). Main stem height(r = −0.104) and lateral branch length (r = −0.115) were negatively correlated with pod yield, but there was no significant correlation (p > 0.05). The results of diameter analysis showed that main stem height (direct diameter coefficient = 1.163) had the greatest direct effect on yield, fol-lowed by hundred pod weight (direct diameter coefficient = 0.945), pod number per plant (direct diameter coefficient = 0.339) and pod yield (direct diameter coefficient = 0.305). Plant height, pod weight and pod number per plant are the main reference traits affecting the yield of Zhenghua 6, which can be used as the index reference for the high yield characteristics of Zhenghua 6, and can be selected together with the pod yield rate and the number of pod branches. In the cultivation of Zhenghua 6, we should match suitable sowing date, select suitable sowing density, establish rea-sonable population structure, strengthen late management and control plant height, so as to achieve high quality and efficient cultivation.
1. 引言
花生(Arachis hypogaea L.)是我国重要的经济作物、油料作物 [1] ，花生种植效益较高，据2020年全国花生平均亩产值达1661元，在大宗粮棉油作物中仅次于棉花，居第二位。大力发展花生生产有助于促进农民增收、保障油料自主供给。花生品种的生态适应性及产量稳定性是花生品种合理布局和推广区划的关键，GGE (genotype + genotype-by-environment interaction)双标图 [2] [3] [4] 已广泛应用于玉米 [5] [6] 、水稻 [7] 、小麦 [8] 、马铃薯 [9] [10] 、大豆 [11] 、燕麦 [12] 、棉花 [13] 的品种区域试验评价。研究表明，农艺性状对花生产量的影响很大，百果重、出米率、百仁重等性状与产量密切相关 [14] [15] [16] [17] [18] ，但由于品种、栽培模式及地理环境的不同，花生不同农艺性状与产量之间的遗传关系具有复杂多样性。国内前人研究基于特定区域的特定品种，指导一定区域内的花生育种及生产具有意义。而不同类型花生品种影响产量的主要性状也不同，尤其是关于高油大果花生品种基因与环境互作与产量构成因素分析相结合的研究较少。本文利用GGE双标图对22份北方大果花生品种(系)进行评价，并运用通径方法对郑花6号的8个主要农艺性状与荚果产量进一步分析 [19] - [27] ，进一步明确农艺性状在产量构成因素中的影响。分析并评价郑花6号的丰产稳产性及其产量构成因素，为高油大果花生选育产量相关优良性状的选择和品种布局的优化，和丰富花生优质高产栽培理论，以及充分发挥品种的增产潜力提供参考。
2. 材料和方法
2.1. 试验材料
以2012~2013年全国北方片花生大粒组区域试验数据为基础，参试品种见表1，2012年试验点为E1~E16，2013年试验点为E1~E18 (表2)。
Table 1. Information of peanut varieties tested
表1. 参试花生品种信息
2.2. 数据分析
利用Excel对参试品种的农艺性状数据进行统计，SPSS17.0软件对郑花6号主要农艺性状作相关分析和通径分析，Genstat-biplot软件 [2] [3] [4] 对不同试点、不同品种产量作双标图分析丰产稳产性及试点代表性。
3. 结果与分析
3.1. 郑花6号的产量评价
两年多点试验共34点次，郑花6号2012年平均产量5420.11 kg/hm2，比对照增产6.20%，16点次中13点增产。2013年平均产量4901.92 kg/hm2，增产9.40%，18点次中17点增产。两年平均产量5161.02 kg/hm2，34个点次中有30个点次分别有不同程度增产，最高增幅为23.98% (E13') (图1)。
Figure 1. Yield performance of Zhenghua 6 in two years multi-site experiment. Note: Horizontal coordinates E1~E16 represent 2012 pilot, E1'~E18' represent 2013 pilot, and vertical coordinates represent output
图1. 郑花6号在两年多点试验中的产量表现。注：横坐标E1~E16代表2012年试点，E1'~E18'代表2013年试点；纵坐标代表产量
3.2. 试验点鉴评
图2环境分布可见，E4、E5、E9 (2012)和E7、E13、E18 (2013)试点与原点之间的方差向量较短，代表2012年密云、青岛和菏泽三个环境较适宜花生产量增长，2013年濮阳、潍坊、安阳三个环境品种产量表现较高。
Figure 2. Identification and representativeness of test sites. Note: The horizontal axis PC1 represents the first principal component, and the vertical axis PC2 represents the second principal component. The same applies to 2012 (left) and 2013 (right)
图2. 试验点的鉴别力和代表性。注：横坐标PC1代表第一主成分，纵坐标PC2代表第二主成分，2012年(左)及2013年(右)下同
3.3. 丰产稳产的最优基因型分析
图3中箭头前的小圆圈代表着平均环境，距离小圆圈最近的品种为最优基因型。2012年产量较稳定的是郑花6号、圣濮花1号、P10-2、P10-8、11测L3、秋乐606，最优基因型为开19-2；2013年产量较稳定的有郑花6号、泛花3号、豫花41号，最优基因型为开19-2。2012年和2013年G14和g12 (郑花6号)平均产量均位居第4，在各试点的稳产性较好。
Figure 3. Analysis of high and stable yield
图3. 丰产稳产性分析
3.4. 综合排名
以箭头为中心作等距离画圆，靠近圆心的品种综合排名更靠前，丰产稳产性越高，2012年丰产稳产较好的有P10-2、开19-2、漯花6号，郑花6号排名第6；2013年丰产稳产性最好的是P10-8，郑花6号排名第2 (图4)。
Figure 4. GGE double label chart analysis comprehensive ranking
图4. GGE双标图分析综合排名
3.5. 品种适应性分析
GGE双标图可以预测环境和基因型的相互适应性，参试品种在同一个所在区域内的试点更适宜种植。图5所示，2012年大连、潍坊和洛阳最适宜种植的材料为漯花6号和郑农花13号，石家庄、易县、菏泽、
Figure 5. Adaptability analysis of peanut varieties and classification of test sites based on GGE double plot
图5. 基于GGE双标图的花生参试品种适应性分析及试验点的类型区划分
郑州、濮阳试点最适种植品种为P10-2，包含近半数环境点(青岛、烟台、济宁、开封、漯河、徐州、固镇)最适宜品种为开19-2，郑花6号适合种植区域也在此区域，密云最适宜种植花育43号。2013年郑花6号和开19-2在包含大多数试验点(大连、密云、青岛、烟台、潍坊、济宁、菏泽、洛阳、漯河、徐州、安阳)的表现最好。
3.6. 郑花6号的产量构成
3.6.1. 主要农艺性状的变异性分析
表3可见，郑花6号主要农艺性状中单株结果数、主茎高、侧枝长的变异系数较大(20.36~21.72%)，较易受环境影响。出米率的变异系数最小(3.89%)，受环境影响较小。其他性状的变异系数在7.50~17.41%之间。
Table 3. Analysis of variation of main agronomic traits in Zhenghua 6
表3. 郑花6号主要农艺性状的变异性分析
3.6.2. 荚果产量与农艺性状的相关分析
表4可见荚果产量与其他8个农艺性状的相关程度 [14] [15] [16] [17] 是百果重 > 百仁重 > 结果枝数 > 出米率 > 总分枝数 > 单株果数 > 主茎高 > 侧枝长，其中总分枝数、百仁重、结果枝数、出仁率、单株果数、百果重与荚果产量均为正相关，侧枝长、主茎高与产量负相关。主茎高与侧枝长呈极显著正相关，两者与单株果数、百果重、百仁重呈负相关；百果重、百仁重之间极显著正相关(p < 0.01)，出米率和结果枝数与总分枝数显著正相关(p < 0.05)。说明控制主茎高与侧枝长可促进单株果数、百果重和百仁重的提高，结果枝数和总分枝数对出米率有重要影响，改良任一性状均可提高出米率。
Table 4. Correlation analysis between pod yield and main agronomic traits of Zhengzhou Flower 6
表4. 郑花6号的荚果产量与主要农艺性状的相关分析
注：*：p < 0.05；**：p < 0.01。X1，主茎高；X2，侧枝长；X3，总分枝数；X4，结果枝数；X5，单株果数；X6，百果重；X7，百仁重；X8，出米率；Y，荚果产量。表5同。
Note: *: p < 0.05；**: p < 0.01. X1, main stem height; X2, lateral branch length; X3, total branching number; X4, number of branches with pods;X5, pod number of per plant;X6, weight of 100 pods; X7, weight of 100 seeds; X8, shelling percentage; Y, pod yield. The same of Table 5.
3.6.3. 主要农艺性状的通径分析
对相关系数进行通径分析，表5表明这8个农艺性状对荚果产量的直接通径系数排序为：主茎高 > 百果重 > 单株果数 > 出米率 > 结果枝数 > 总分枝数 > 百仁重 > 侧枝长。其中主茎高的直接通径系数最大，为1.163，是影响郑花6号产量的最大因素。根据间接通径系数结果，主茎高通过对百果重、单株果数、出米率、结果枝数的正调控，对总分枝数、百仁重、侧枝长的负调控来间接影响荚果产量；且正调控大于负调控的绝对值，说明主茎高对荚果产量的影响最大。百果重对产量影响的直接通径系数也较大，为0.945，它通过对总分枝数、百果重、百仁重、出米率的正调控，对主茎高、侧枝长、结果枝数、单株果数的负调控来间接影响荚果产量。说明农艺性状之间通过互相制约和促进的关系，共同影响荚果产量。
Table 5. Path analysis of main agronomic traits on pod yield of Zhenghua 6
表5. 郑花6号主要农艺性状对荚果产量的通径分析
4. 品质性状
农业部油料及制品质量监督检验测试中心检测(表6)结果表明，郑花6号粗蛋白含量22.32%，油酸含量41.55%，亚油酸含量36.25%，油亚比1.15，平均粗脂肪含量55.74%，超过中国花生脂肪含量55%的优质标准，为高脂肪花生新品种。
Table 6. The quality performance of Zhenghua 6
表6. 郑花6号的品质表现
5. 讨论
5.1. 基因型与环境互作效应对花生品种丰产稳产性的影响
品种的多点试验是多地点、多年份的联合试验，不同的影响因素互相作用，适宜的分析方法可以为合理的品种布局提供依据。本研究花生产量双标图显示，基因型、环境互作效应占总变异的45.73%~67.28%。根据产量适应性GGE分析，郑花6号在大连、密云、青岛、烟台、潍坊、济宁、菏泽、洛阳、漯河、徐州、安阳、固镇等试验点适应性较好，即适宜在辽宁、河北、山东、河南、江苏、安徽等省份推广种植。2012年在16个试点中接近理想的参试品系为P10-2、开19-2；2013年在18个试点中接近理想的参试品系为P10-8、郑花6号。较为理想的试点为密云、青岛、菏泽、濮阳和潍坊、安阳，其代表性及鉴别力均较好。
5.2. 花生品种主要农艺性状的变异性、相关性分析
不同品种来源的花生的农艺性状差异较大。冯亚平等 [15] 对59份大花生品种的11个农艺性状进行变异分析，变异幅度为4.83~38.61%，最小的为出米率，较大的为总分枝数、单株饱果数、单株秕果数。苗建利等 [16] 对开农82的9个农艺性状及荚果产量进行变异分析，变异系数最低的为出仁率(4.17%)，较高的为主茎高、结果枝数、分枝数、结果数、饱果数。文中郑花6号的8个农艺性状和荚果产量的变异幅度为3.89~21.72%，各农艺性状的变异系数由大到小为单株结果数、主茎高、侧枝长、百果重、结果枝数、百仁重、荚果产量、总分枝数、出米率。单株结果数、主茎高、侧枝长的变异系数分别为21.72%、20.60%、20.36%，表明受环境的影响大，出米率的变异系数最小为3.89%，受环境影响较小。这与冯亚平等 [15] 关于出米率(CV，4.83%)、结果枝数(CV，15.48%)、单株结果数(CV，20.86%)、百果重(CV，15.10%)、百仁重(CV，13.96%)，苗建利等 [25] 关于出米率(CV，4.17%)、百仁重(CV，14.85%)、主茎高(CV，17.53%)、结果枝数(CV，19.45%)的研究结果相近，其余农艺性状的变异系数差异较大，可能与选取的不同研究对象及栽培环境相关。因此，品种选育鉴定关键农艺性状时，要兼顾多个性状协同选择。
大力发展花生等油料作物在近几年的中央一号文件中多次提及，2023年中央一号文件提出今年我国粮食生产目标任务重点是“两稳两扩两提”，其中包括扩油料、稳产量。出油率高、产量高的花生品种的油用经济效益表现突出，本研究以高脂肪花生品种郑花6号为例，对其产量构成要素进行了分析。通过对高产稳产的郑花6号的8个农艺性状与产量之间的相关性进行分析，相关程度表现为百果重 > 百仁重 > 结果枝数 > 出米率 > 总分枝数 > 单株果数 > 主茎高 > 侧枝长；主茎高、侧枝长与其荚果产量负相关，主茎高与侧枝长呈显著正相关，两者与单株果数、百果重、百仁重呈负相关；百果重和百仁重显著正相关，总分枝数与出米率、结果枝数显著正相关，这与前人研究结果一致 [15] [16] [17] [18] 。表明新品种选育过程中要重点考虑育种材料的百果重、结果枝数及总分枝数，并控制好株高及侧枝长的选择。
本研究发现主茎高和百果重是影响郑花6号产量的主要因素，且相关性分析和通径分析的结果并不完全一致。这与邓丽 [26] 等的研究结果相似。说明农艺性状具有复杂多样性，应综合考量各性状对产量的影响。因此在生产实践中应根据具体的品种特点，配套相应的栽培管理措施以提高品种的增产潜力。
6. 结论
6.1. 郑花6号稳产性分析及试验点鉴别
郑花6号在大连、密云、青岛、烟台、潍坊、济宁、菏泽、洛阳、漯河、徐州、安阳、固镇等试验点丰产稳产性较好，适宜在辽宁、河北、山东、河南、江苏、安徽等省份种植推广。密云、青岛、菏泽、濮阳、潍坊、安阳试点的代表性及鉴别力均较好。
6.2. 郑花6号产量构成要素
郑花6号产量构成要素(百果重、百仁重、结果枝数、出米率、总分枝数、单株果数、主茎高、侧枝长)对荚果产量的影响程度存在明显差异，其中，主茎高、百果重、单株果数是影响郑花6号产量的主要因素，可以作为郑花6号高产特征指标的参考。
6.3. 高油花生品种郑花6号适宜推广
在郑花6号栽培中应配合适宜播期，择适宜播种密度，建立合理群体结构，加强后期管理，控制株高，兼顾百果重和单株果数的增加。郑花6号脂肪含量达到54.54~56.94%，属高油花生新品种，郑花6号等高油花生品种的选育和推广可为我国花生油及其他功能产品供给提供保障，推动花生产业的发展。
基金项目
河南省重大专项(豫科计[2015]141100110600)；河南优势特色农业产业科技支撑行动计划项目：优质花生产能提升与新型推广体系建设(20230202002)；郑州市科技惠民计划项目：粮油作物新品种及病虫草害绿色防控技术示范应用(2023KJHM0042)。
参考文献
NOTES
*第一作者。
#通讯作者。