长江流域甘蓝型油菜育种亲本农艺性状的遗传变异、相关及主成分分析

郑本川 崔成 李浩杰 张锦芳 柴靓 蒋俊 蒋梁材

摘要： 【目的】研究長江流域地区甘蓝型油菜育种亲本农艺性状的遗传变异，为利用不同区域甘蓝型油菜育种亲本材料的优良农艺性状、配制和选育优良品种（组合），提高油菜产量提供参考依据。【方法】以国家油菜产业技术体系提供的从长江上、中、下游区域征集的213份甘蓝型油菜育种亲本为材料，采用相关分析、通径分析和主成分分析等方法分析四川省成都平原气候条件下供试甘蓝型油菜亲本材料主要农艺性状的遗传变异及其与单株产量的关系，明确影响甘蓝型油菜育种亲本单株产量形成的主要农艺性状。【结果】长江上、中、下游区域甘蓝型油菜育种亲本农艺性状变异系数最大的均为单株产量，分别为18.29%、18.91%和18.47%;变异系数最小的均为株高，分别为6.24%、6.44%和5.95%。不同区域甘蓝型油菜育种亲本的农艺性状除上游的区域千粒重与单株产量呈负相关外，其余性状均与单株产量呈正相关，其中株高、主序有效长度、一次有效分枝数和每果粒数与单株产量呈显著（P<0.05，下同）或极显著（P<0.01，下同）正相关;上游和中游区域的主序有效角果数与单株产量呈极显著正相关，角果长度与单株产量显著正相关;中游区域的千粒重与单株产量极显著正相关。上游区域对单株产量直接作用最大的是株高，间接作用最大的是有效分枝高度;中游区域对单株产量直接作用最大的是每果粒数，间接作用最大的是株高;下游区域对单株产量直接作用最大的是主序有效长度，间接作用最大的是主序有效角果数。主成分分析可将上游和中游区域相关性状分为三大类，上游区域分别是株高性状、产量构成性状和角果性状，贡献率分别为41.353%、15.186%和12.083%;中游区域分别是株高性状、角果性状和分枝性状，贡献率分别为31.485%、22.004%和14.097%;下游区域分为四大类，分别是株高性状、角果性状、主序性状和产量构成性状，贡献率分别为32.393%、19.971%、16.035%和11.960%。【结论】在四川省成都平原气候条件下，影响长江上游区域甘蓝型油菜育种亲本材料的主要性状是株高，中游区域主要是每果粒数，下游区域主要是主序长度。

关键词： 甘蓝型油菜;农艺性状;遗传变异;相关分析;主成分分析;长江流域

中图分类号： S634.3                              文献标志码： A 文章编号：2095-1191（2019）10-2196-09

Genetic variation， correlation and principal component analysis of agronomic traits of breeding parents of Brassica napus L.

in Yangtze River

ZHENG Ben-chuan1，2， CUI Cheng1， LI Hao-jie1， ZHANG Jin-fang1，

CHAI Liang1， JIANG Jun1， JIANG Liang-cai1*

（1Crops Research Institute， Sichuan Academy of Agricultural Science， Chengdu  610066，China; 2College of Agronomy， Sichuan Agricultural University， Wenjiang， Sichuan  611130， China）

Abstract：【Objective】This study aimed to reveal the heritable variation of agronomic traits in rapeseed （Brassica napus L.） from Yangtze River area， and provide a reference base for utilizing superior agronomic traits of parental materials from different regions， breeding fine varieties（combinations） and improving the yield of rapeseed. 【Method】A total of 213 cultivars/lines collected from the upper， middle and lower reaches of Yangtze River Basin by national rapeseed industrial technique system were used in this study. The correlation analysis， path analysis and principal component analysis were performed to analyze the relationship between the genetic variation of main agronomic characters and yield per plant under Chengdu plain climatic conditions in Sichuan， the main agronomic characters effecting the yield per plant were then confirmed. 【Result】The greatest variation coefficients of parental agronomic traits in the upper， middle and lower reaches of the Yangtze river was the yield per plant， representing 18.29%， 18.91% and 18.47%， respectively; while the lowest variation coefficient was the plant height， accounting for 6.24%， 6.44% and 5.95%， respectively. The agronomic traits in different regions were positively correlated with the yield per plant， except for the 1000-seed weight in the upper reaches of the Yangtze River， which was negatively correlated with the yield per plant. The plant height， the effective length of main inflorescence， the number of first effective branches and seeds per silique were significantly（P<0.05， the same below） or extremely（P<0.01， the same below） positively correlated with the yield per plant. The number of effective silique on main inflorescence in the upper and middle reaches of the Yangtze River was extremely positively correlated with the yield per plant， and the length of silique was significantly positively correlated with yield per plant. The 1000-seed weight in the middle reaches of the Yangtze River was positively correlated with the yield per plant. In the upper reaches of the Yangtze River， plant height had the largest direct effects on yield per plant and effective branch height had the largest indirect effects on yield per plant. In the middle reaches of the Yangtze River， seeds per silique had the largest direct effects on yield per plant and plant height had the largest indirect effects on yield per plant. In the upper reaches of Yangtze River， the length of main inflorescence had the largest direct effects on yield per plant and number of effective silique on main inflorescence had the largest indirect effects on yield per plant. The related traits in the upper and middle reaches of the Yangtze River could be clustered into three categories. In the upper reaches of the Yangtze River， the three categories were plant height traits， yield components traits and silique traits with contribution rates of 41.353%， 15.186% and 12.083%， respectively. In the middle reaches of the Yangtze River， the three were plant height traits， silique traits and branch traits， and their contribution rates were 31.485%， 22.004% and 14.097%， respectively. The main components in the lower reaches of the Yangtze River were classified in four categories， they were  plant height traits， silique traits， main inflorescence traits and yield components traits， the contribution rates were 32.393%，19.971%，16.035% and 11.960%， respectively. 【Conclusion】The most important traits affecting rapeseed breeding parents in the upper， middle and lower reaches of Yangtze River under Chengdu plain climate in Sichuan are plant height， seed number per silique and main inflorescence length， respectively.

Key words： Brassica napus L.; agronomic traits; genetic variation; correlation analysis; principal component ana-lysis; Yangtze River

0 引言

【研究意义】油菜是世界主要的油料作物，我国作为油菜的主要种植和消费国之一，常年种植面积和产量均较高，但仍无法满足国内消费需求，因此，提高油菜产量是油菜育种的首要目标（王汉中，2004;涂金星等，2007）。油菜产量除受产量构成因素如单株有效角果数、每角果粒数和千粒重等直接影响外，还受其他因素的间接影响，不同地域油菜种质含有可供育种使用的优良性状，而充分利用这些性状对于品种改良和配制优良油菜品种（组合），以及提高油菜产量均具有重要意义。【前人研究进展】前人已开展了较多有关油菜种质或品种农艺性状与产量关系的研究，且证实油菜产量与单株有效角果数、每角粒数、千粒重、主花序长度、主花序有效角果数、结角密度、分枝数及分枝部位高度等显著相关（张芳等，2012;白桂萍等，2015;李宏军等，2015;宋丰萍等，2015;倪正斌等，2017）。在油菜农艺性状变异研究方面，张文英和王凯华（2012）研究表明，甘蓝型油菜植株性状中分枝部位的变异最大、千粒重的变异最小;王健胜等（2015）研究表明，供试18个油菜亲本的11个主要农艺性状中单株产量的變异系数最大、株高最小。在研究分析方法方面，黄桃翠等（2016）通过回归分析和灰色关联分析探讨油菜农艺性状与产量的密切程度，利用Pearson相关分析探讨油菜农艺性状与产量性状间的线性关系;许进鸿等（2016）运用灰色系统理论对影响白菜型油菜产量的主要农艺性状进行分析。【本研究切入点】关于甘蓝型油菜产量与农艺性状相关性的研究已有大量报道，但主要针对本生态区主栽品种或育种亲本进行分析，遗传变异范围较窄，导致选育品种适应范围较小、适应性较差。目前针对不同生态区油菜育种亲本材料优良性状利用的研究报道较少。【拟解决的关键问题】以国家油菜产业技术体系提供的从长江上、中、下游区域征集的213份甘蓝型油菜育种亲本为材料，采用相关分析、通径分析和主成分分析等方法，探讨在四川省成都平原气候条件下供试甘蓝型油菜亲本材料主要农艺性状的遗传变异及与单株产量的关系，明确影响甘蓝型油菜育种亲本单株产量形成的主要农艺性状，为利用不同区域的甘蓝型油菜育种亲本材料的优良农艺性状、配制和选育优良品种（组合）及提高油菜产量提供参考依据。

1 材料与方法

1. 1 试验材料

供试的213份甘蓝型油菜育种亲本材料由国家油菜产业技术体系提供，其中，来源于长江上游区域的亲本材料41份、中游区域144份、下游区域28份。供试甘蓝型油菜亲本材料信息见表1。

1. 2 试验方法

1. 2. 1 试验设计 试验材料分别于2012—2013和2013—2014年度在四川省成都市新都区泰兴镇四川省农业科学院科技创新示范园区内种植，株行距0.25 m×0.25 m，完全随机试验设计，3次重复，每个材料种5行，每行8株，田间管理同大田生产。

1. 2. 2 测定项目及方法 于油菜成熟期每个材料取4株进行农艺性状考察，主要指标有株高、有效分枝高度、主序有效长度、一次有效分枝数、主序有效角果数、角果长度、每果粒数、千粒重和单株产量。指标测定参考郑本川等（2013）的方法。

1. 3 统计分析

试验数据采用Excel 2010进行统计整理，应用SPSS 22.0计算各农艺性状的平均值、标准差和变异系数，相关性分析采用Pearson双侧T检验。通径分析和主成分分析参考陶爱芬等（2008）、杜家菊和陈志伟（2010）、周丽艳等（2011）、杜鹃（2012）、宋小园等（2016）的方法。

2 结果与分析

2. 1 不同区域甘蓝型油菜育种亲本农艺性状变异分析结果

由表2可知，长江上、中、下游区域甘蓝型油菜亲本农艺性状变异系数最大的均为单株产量，分别为18.29%、18.91%和18.47%;变异系数最小的均为株高，分别为6.24%、6.44%和5.95%。上游区域甘蓝型油菜亲本农艺性状中有效分枝高度、一次有效分枝数、主序有效角果数和千粒重的变异系数分布在11.06%～11.99%;中游区域甘蓝型油菜亲本农艺性状中有效分枝高度、一次有效分枝数、角果长度和千粒重的变异系数分布在10.27%～14.68%;下游区域甘蓝型油菜亲本农艺性状中有效分枝高度、主序有效长度、主序有效角果数和千粒重的变异系数分布在9.23%～12.27%。说明单株产量、分枝高度、千粒重、一次有效分枝数和主序有效角果数在育种选育中具有一定的改进余地和利用价值。

2. 2 不同区域甘蓝型油菜育种亲本农艺性状与单株产量的相关性分析结果

由表3可知，不同区域甘蓝型油菜亲本的农艺性状除上游区域的千粒重与单株产量呈负相关外，其余性状均与单株产量呈正相关，其中株高、主序有效长度、一次有效分枝数和每果粒数与单株产量呈显著（P<0.05，下同）或极显著（P<0.01，下同）正相关;上游和中游区域的主序有效角果数与单株产量呈极显著正相关，角果长度与单株产量显著正相关;中游区域的千粒重与单株产量极显著正相关。表明株高、主序有效长度、一次有效分枝数和每果粒数等的增加对单株产量的提高具有显著促进作用。

2. 3 不同区域甘蓝型油菜育种亲本农艺性状与单株产量的通径分析结果

2. 3. 1 长江上游区域甘蓝型油菜育种亲本农艺性状与单株产量的通径分析结果 由表4可知，甘蓝型油菜各主要农艺性状对单株产量的直接作用排序为株高（1.104）>有效分枝高度（[-0.530]）>主序有效长度（[-0.333]）>每果粒数（0.271）>千粒重（[-0.200]）>一次有效分枝数（[-0.135]）>主序有效角果数（0.091）>角果长度（0.043）。间接综合效应作用排序为有效分枝高度（0.761）>主序有效长度（0.760）>一次有效分枝数（0.460）>株高（[-0.450]）>主序有效角果数（0.433）>角果长度（0.290）>每果粒数（0.278）>千粒重（0.183）。其中，株高、主序有效长度、主序有效角果数和每果粒数与单株产量呈极显著正相关，说明这些性状的改变对单株产量的提高具有显著促进作用。

2. 3. 2 长江中游区域甘蓝型油菜育种亲本农艺性状与单株产量的通径分析结果 由表5可知，甘蓝型油菜各主要性状对单株产量的直接作用排序为每果粒数（0.286）>一次有效分枝数（0.216）>千粒重（0.191）>主序有效长度（0.183）>主序有效角果数（0.151）>株高（0.146）>有效分枝高度（[-0.062]）>角果长度（[-0.020]）。间接综合效应作用排序为株高（0.282）>主序有效长度（0.225）>角果长度（0.205）>主序有效角果数（0.201）>有效分枝高度（0.149）>千粒重（0.080）>一次有效分枝数（0.073）>每果粒数（0.043）。其中，株高、主序有效角果数、每果粒数、一次有效分枝数、千粒重和主序有效长度与单株产量呈极显著正相关，说明这些性状的改变对提高单株产量具有促进作用。

2. 3. 3 长江下游区域甘蓝型油菜育种亲本农艺性状与单株产量的通径分析结果 由表6可知，甘蓝型油菜各性状对单株产量的直接作用排序为主序有效长度（0.627）>一次有效分枝数（0.417）>主序有效角果数（[-0.397]）>株高（0.207）>千粒重（0.164）>角果长度（0.109）>每果粒数（0.054）>有效分枝高度（0.041）。间接综合作用排序为主序有效角果数（0.481）>每果粒数（0.381）>株高（0.332）>千粒重（0.236）>主序有效长度（[-0.187]）>角果长度（0.177）>一次有效分枝数（0.109）>有效分枝高度（0.089）。其中，主序有效长度和每果粒数与单株产量呈显著正相关，一次有效分枝数和株高与单株产量呈极显著正相关，说明这些性状的改变对提高单株产量具有促进作用。

2. 4 不同区域甘蓝型油菜育种亲本主要性状的主成分分析结果

由表7可知，长江上游和中游区域均提取出3个主成分，其累计贡献率分别为68.622%和67.586%，下游区域提取出4个主成分，其累计贡献率为80.359%。

长江上游区域第一主成分（PC1）中特征值较大且为正值的是株高（0.910），因此第一主成分综合为株高性状，其贡献率为41.353%;第二主成分（PC2）中特征值较大且为正值的是千粒重（0.649），因此第二主成分综合为产量构成性状，其贡献率为15.186%;第三主成分（PC3）中特征值较大且为正值的是角果长度（0.681），因此第三主成分综合为角果性状，其贡献率为12.083%。

长江中游区域第一主成分（PC1）中特征值较大且为正值的是株高（0.882），因此第一主成分综合为株高性状，其贡献率为31.485%;第二主成分（PC2）中特征值较大且为正值的是角果长度（0.832），因此第二主成分综合为角果性状，其贡献率为22.004%;第三主成分（PC3）中特征值较大且为正值的是一次有效分枝数（0.609），因此第三主成分可综合为分枝性状，其贡献率为14.097%。

长江下游区域第一主成分（PC1）中特征值较大且为正值的是株高（0.901），因此第一主成分可综合为株高性状，其贡献率为32.393%;第二主成分（PC2）中特征值较大且为正值的是角果长度（0.610），因此第二主成分可综合为角果性状，其贡献率为19.971%;第三主成分（PC3）中特征值较大且为正值的是主序长度（0.509），因此第三主成分可综合为主序性状，其贡献率为16.035%;第四主成分（PC4）中特征值较大且为正值的是千粒重（0.677），因此第四主成分可综合为产量构成性状，其贡献率为11.960%。

综上所述，不同区域甘蓝型油菜育种亲本第一主成分均是株高性状，但上游区域贡献率大于中、下游区域;上游区域的第二主成分和下游区域第四主成分为产量构成性状，上游区域第三主成分和中、下游区域的第二主成分均为角果性状，中、下游区域第三主成分分别为分枝性状和主序性状。

3 讨论

高产是油菜育种工作一直追求的目标。油菜产量除受产量构成性状的直接影响外，还受其他性状的间接作用。采用统计学方法对油菜种质资源农艺性状和产量关系进行鉴定和评价，可充分利用优良农艺性状进行亲本改良和组合配制。为此，本研究以长江上、中、下游区域不同油菜主产区征集的213份甘蓝型油菜育种亲本为材料，选取8个与单株产量相关的性状进行遗传变异分析、相关分析、通径分析和主成分分析。

对213份甘蓝型油菜育种亲本主要农艺性状的变异分析结果表明，不同区域甘蓝型油菜亲本农艺性状变异系数最大的为单株产量、最小的为株高，与王健胜等（2015）的研究结果一致，因此，提高单株产量对油菜增产具有重要意义。而张文英和王凯华（2012）研究表明，甘蓝型油菜植株性状中分枝部位变异最大，千粒重变异最小。究其原因可能是研究材料选取不同所致，本研究选用的是育种亲本，而张文英和王凯（2011）选择的是杂交种或常规品种。此外，上游区域甘蓝型油菜的有效分枝高度、一次有效分枝数、主序有效角果数和千粒重的变异系数分布在11.06%～11.99%，中游区域甘蓝型油菜的有效分枝高度、一次有效分枝数、角果长度和千粒重的变异系数分布在10.27%～14.68%，下游区域甘蓝型油菜的有效分枝高度、主序有效长度、主序有效角果数和千粒重的变异系数分布在9.23%～12.27%，说明供试不同区域甘蓝型油菜育种亲本材料在可改良的农艺性状方面略有差异，可进一步改造利用。

关于油菜产量与农艺性状的相关性研究，杨安中和彭春华（2006）、张文英和王凯华（2012）研究表明，單株产量与单株有效角果数、一次有效分枝数及千粒重等呈显著正相关;杨安中和彭春华（2006）研究表明，单株产量与每角果粒数呈极显著负相关;关周博等（2013）研究表明，单株产量除了与单株有效角果数、每角粒数等产量构成性状呈显著正相关外，还与株高呈显著正相关。本研究结果表明，不同区域甘蓝型油菜育种亲本的农艺性状除上游区域的千粒重与单株产量呈负相关外，其余性状均与单株产量呈正相关，其中株高、主序有效长度、一次有效分枝数和每果粒数均与单株产量呈显著（或极显著）正相关，上游和中游区域的油菜主序有效角果数与单株产量极显著正相关，角果长度与单株产量呈显著正相关，中游区域的油菜千粒重与单株产量呈极显著正相关。这一结果与前人的研究存在异同，主要是因为本研究选取的研究对象是长江流域不同区域部分育种亲本，前人研究选取的是杂交种或常规品种，研究对象的特性存在一定差异。

关于油菜产量的影响因素研究，张锦芳等（2007）、关周博等（2013）认为，对单株产量的直接、间接影响较大的性状为单株有效角果数、每角粒数;宋稀等（2010）認为不同种植密度品种对单株产量的影响不同，常规密度品种对产量影响最大的是单株有效角果数，其次是株高，而耐密植油菜对单株产量直接作用最大的是结角密度，其次是分枝高度和株高;李宏军等（2015）认为，在低密度条件下对产量影响最大的是株高，而在高密度条件下对产量影响最大的是每果粒数;黄益国等（2017）认为，对早熟油菜品种产量直接影响最大的是千粒重，间接影响最大的是主花序长，对中熟油菜品种产量直接影响最大的是株高，间接影响最大的是单株有效角果数。本研究中，不同区域甘蓝型油菜亲本对单株产量的影响因素不同，长江上游区域育种亲本对单株产量直接作用最大的是株高，其次是有效分枝高度和主序有效长度，间接作用最大的是有效分枝高度，其次是主序有效长度和一次有效分枝数;长江中游区域对单株产量直接作用最大的是每果粒数，其次是一次有效分枝数和千粒重，间接作用最大的是株高，其次是主序有效长度和角果长度;长江下游区域对单株产量直接作用最大的是主序有效长度，其次是一次有效分枝数和主序有效角果数，间接作用最大的是主序有效角果数，其次是每果粒数和株高。研究结果不同的原因主要是不同研究中所选取的考察性状不同，本研究中作为产量构成性状的单株有效角果数没有进行考种，因此在进行分析时未给出单株有效角果数对单株产量的贡献。

主成分分析是将原来的多个性状因子归纳为少数几个性状进行分析。本研究结果表明，长江上游和中游区域均提取出3个主成分，其累计贡献率分别为68.621%和67.586%，下游区域提取出4个主成分，其累计贡献率为80.359%。各区域第一主成分均为株高性状，且上游区域的贡献率最高;上游区域的第二主成分是产量构成性状，中、下游区域为角果性状;上游区域第三主成分为角果性状，中游区域为分枝性状，下游区域为主序性状;下游区域第四主成分为产量构成性状。

4 结论

在四川省成都平原气候条件下，影响长江上游区域甘蓝型油菜育种亲本材料的主要性状是株高，中游区域主要是每果粒数，下游区域主要是主序长度。

参考文献：

白桂萍，刘克钊，谭永强，尹羽丰，余华强，王汉中. 2015. 油菜高产群体各农艺性状对产量的影响[J]. 作物杂志，（6）：33-38. [Bai G P，Liu K Z，Tan Y Q，Yin Y F，Yu H Q，Wang H Z. 2015. Effect of agronomic traits on seed yield in high-yielding rapeseed populations[J]. Crops，（6）：33-38.]

杜家菊，陈志伟. 2010. 使用SPSS线性回归实现通径分析的方法[J]. 生物学通报，45（2）：4-6. [Du J J，Chen Z W. 2010. Method for path analysis using SPSS linear regression[J]. Biology Bulletin，45（2）：4-6.]

杜鹃. 2012. 通径分析在Excel和SPSS中的实现[J]. 陕西气象，（1）：15-18. [Du J. 2012. Path analysis in Excel and SPSS[J]. Journal of Shaanxi Meteorology，（1）：15-18.]

关周博，田建华，郑磊，姚雪雁，韦世豪，李少钦，董育红，李殿荣. 2013. 适宜机械化栽培的甘蓝型油菜农艺性状与单株产量的相关性分析及耐密油菜育种探讨[J]. 中国农学通报，29（18）：79-83. [Guan Z B，Tian J H，Zheng L，Yao X Y，Wei S H，Li S Q，Dong Y H，Li D R. 2013. Correlation analysis of agronomic traits and yield per plant in Brassica napus L. of suitable mechanization planting and study of breeding of tolerance close planting canola[J]. Chinese Agricultural Science Bulletin，29（18）：79-83.]

黄桃翠，李加纳，唐世义. 2016. 2014—2015年度长江中游组区试油菜稳定性探析[J]. 中国油料作物学报，（4）：423-430. [Huang T C，Li J N，Tang S Y. 2016. Stability of oilseed rape varieties in 2014-2015 national rapeseed variety test in regions along Midway of Yangtze River[J]. Chinese Journal of Oil Crop Sciences，（4）：423-430.]

黄益国，张学昆，张毅，雷利琴. 2017. 不同熟期油菜品种农艺性状与产量的相关分析[J]. 作物研究，31（3）：260-264.[Huang Y G，Zhang X K，Zhang Y，Lei L Q. 2017. Correlation analysis about agronomic traits and yield of rape varieties with different maturity stages[J]. Crop Research，31（3）：260-264.]

李宏军，杨鸿，朱传霞，胡金荣，薛高尚，黄琳，朱俊子，罗晓玲，张洋，文奕峰，魏廷龙. 2015. 湖南省油菜品种农艺性状分析[J]. 作物杂志，（3）：41-44. [Li H J，Yang H，Zhu C X，Hu J R，Xue G S，Huang L，Zhu J Z，Luo X L，Zhang Y，Wen Y F，Wei T L. 2015. Analysis on agronomic traits of rape varieties in Hunan Province[J]. Crops，（3）：41-44.]

倪正斌，孙红芹，万林生. 2017. 甘蓝型油菜产量与主要农艺性状的灰色关联度分析[J]. 浙江农业科学，58（7）：1146-1149. [Ni Z B，Sun H Q，Wan L S. 2017. Grey relational analysis of yield and main agronomic traits of Brassica napus[J]. Journal of Zhejiang Agricultural Sciences，58（7）：1146-1149.]

宋豐萍，蒙祖庆，窦胜玮，刘丹. 2015. 春播半冬性甘蓝型油菜光温因子与产量及农艺性状的典型相关分析[J]. 中国生态农业学报，23（8）：987-993. [Song F P，Meng Z Q，Dou S W，Liu D. 2015. Canonical correlations of light and temperature with yield and agronomic traits of semi-winter rapeseed（Brassica napus L.） sowed in spring[J]. Chinese Journal of Eco-Agriculture，23（8）：987-993.]

宋稀，刘凤兰，郑普英，张学昆，陆光远，付桂萍，程勇. 2010. 高密度种植专用油菜重要农艺性状与产量的关系分析[J]. 中国农业科学，43（9）：1800-1806. [Song X，Liu F L，Zheng P Y，Zhang X K，Lu G Y，Fu G P，Cheng Y. 2010. Correlation analysis between agronomic traits and yield of rapeseed（Brassica napus L.） for high-density plan-ting[J]. Scientia Agricultura Sinica，43（9）：1800-1806.]

宋小园，朱仲元，刘艳伟，赵宏瑾. 2016. 通径分析在SPSS 逐步线性回归中的实现[J]. 干旱区研究，33（1）：108-112.[Song X Y，Zhu Z Y，Liu Y W，Zhao H J. 2016. Application of path analysis in stepwise linear regression SPSS[J]. Arid  Zone Research，33（1）：108-112.]

陶爱芬，祁建民，林培青，方平平，吴建梅，林荔辉. 2008. 红麻优异种质产量和品质性状主成分聚类分析与综合评价[J]. 中国农业科学，41（9）：2859-2867.[Tao A F，Qi J M，Lin P Q，Fang P P，Wu J M，Lin L H. 2008. Cluster ana-lysis and evaluation of elite kanaf germplasm based on principal components[J]. Scientia Agricultura Sinica，41（9）：2859-2867.]

涂金星，张冬晓，张毅，傅廷栋. 2007. 我国油菜育种目标及品种审定问题的商榷[J]. 中国油料作物学报，29（3）：350-352.[Tu J X，Zhang D X，Zhang Y，Fu T D. 2007. Discussion on some standards of variety registration and breeding goals of Brassica napus in China[J]. Chinese Journal of Oil Crop Sciences，29（3）：350-352.]

王汉中. 2004. 中国油菜品种改良的中长期发展战略[J]. 中国油料作物学报，26（3）：98-101. [Wang H Z. 2004. Medium and long-term development strategy of rape variety improvement in China[J]. Chinese Journal of Oil Crop Sciences，26（3）： 98-101.]

王健胜，梁亚红，侯桂玲，韦世豪，赵干卿. 2015. 油菜育种亲本主要性状相关性及主成分分析[J]. 中国农学通报，31（27）：147-152. [Wang J X，Liang Y H，Hou G L，Wei S H，Zhao G Q. 2015. Correlation and principal components analysis of the main agricultural traits of rapeseed bree-ding parents[J]. Chinese Agricultural Science Bulletin，31（27）：147-152.]

许进鸿，安克杰，次德吉，蒙祖庆，旦巴. 2016. 西藏部分地方白菜型油菜产量与主要农艺性状灰色关联分析[J]. 西藏科技，（1）：72-75. [Xu J H，An K J，Ci D J，Meng Z Q，Dan B. 2016. Grey correlation analysis between the yield of cabbage rapeseed and main agronomic characters in Tibet[J]. Tibet Science and Technology，（1）：72-75.]

杨安中，彭春华. 2006. 油菜单株产量与若干农艺性状的相关分析[J]. 安徽农学通报，12（2）：33-34. [Yang A Z，Peng C H. 2006. Correlation aralysis between agrieultural cha-racters and yield per plant in rape[J]. Anhui Agricultural Science Bulletin，12（2）：33-34.]

张芳，赵永国，谷铁城，张冬晓，刘凤兰，郭瑞星，付桂萍，张学昆. 2012. 2001—2010年国家审定冬油菜品种的产量与主要性状分析[J]. 中国油料作物学报，34（3）：239-244.[Zhang F，Zhao Y G，Gu T C，Zhang D X，Liu F L，Guo R X，Fu G P，Zhang X K. 2012. Yield and agronomic traits of winter rapeseed cultivars registered in China from 2001 to 2010[J]. Chinese Journal of Oil Crop Scien-ces，34（3）：239-244.]

张锦芳，蒲晓斌，李浩杰，张启行，蒋梁材. 2007. 不同来源甘蓝型油菜主要农艺性状与产量的相关分析[J]. 西南农业学报，20（4）：587-590. [Zhang J F，Pu X B，Li H J，Zhang Q X，Jiang L C. 2007. Correlation analysis between major agronomil characters and yield per plant in rapeseed（Brassica nanpus L.） from different source[J]. Southwest China Journal of Agricultural Science，20（4）：587-590.]

张文英，王凯华. 2012. 甘蓝型油菜农艺性状相关性及主成分分析[J]. 湖北农业科学，51（11）：2183-2185. [Zhang W Y，Wang K H. 2012. Correlation and principle component analysis on agronomic trait of Brassica napus L.[J]. Hubei Agricultural Sciences，51（11）：2183-2185.]

郑本川，张锦芳，李浩杰，蒲晓斌，崔成，柴靓，蒋俊，牛应泽，蒋梁材. 2013. 甘蓝型油菜生育期天数与产量构成性状的相关分析[J]. 中国油料作物学报，35（3）： 240-245.[Zheng B C，Zhang J F，Li H J，Pu X B，Cui C，Chai L，Jiang J，Niu Y Z，Jiang L C. 2013. Correlation between duration of growth periods and yield components of Brassica napus L.[J]. Chinese Journal of Oil Crop Scien-ces，35（3）： 240-245.]

周麗艳，郭振清，马玉玲，东方阳，林小虎. 2011. 春小麦品种农艺性状的主成分分析与聚类分析[J]. 麦类作物学报，31（6）：1057-1062. [Zhou L Y，Guo Z Q，Ma Y L，Dongfang Y，Lin X H. 2011. Principal component and cluster analysis of different spring wheat cultivars based on agronomic traits[J]. Journal of Triticeae Crops，31（6）：1057-1062.]

（责任编辑 麻小燕）