菜用大豆主要农艺性状的相关性、聚类及主成分分析

贺礼英， 尹成杰， 黄守程， 何庆元，舒英杰

(安徽科技学院 农学院， 安徽 凤阳 233100)

菜用大豆是豆科(Leguminosae)大豆属的栽培大豆种 (Glycinemax[L.]Merrill)，简称毛豆，是指R6至R7期豆荚鼓粒饱满，荚色、籽粒呈翠绿色时进行采收的大豆类型，属于特用大豆[1-2]。菜用大豆中富含蛋白、脂肪、游离氨基酸及维生素，对调节人们膳食结构和改善营养状况具有重要的作用[3-4]，因其营养价值高、口感好，被公认是无公害或少污染的安全保健食品[5]。随着人们生活水平的提高和膳食结构的改变，尤其是近年农业供给侧结构性改革后，我国大豆(包括菜用大豆)的播种面积明显增加，到2020年全国大豆种植面积将达到9.3×106hm2[6]，这对菜用大豆产业的发展是难得的机遇，但是目前市场上菜用大豆品种较多、良莠不齐，而且地方性品种不明显，因此选择适合当地种植的菜用大豆品种并进行高效栽培研究对该地区菜用大豆产业的可持续发展具有重要的理论价值和实践指导意义。

聚类和主成分分析是研究植物种质资源农艺性状常用的方法，可以对参试种质进行综合客观评价。本试验通过对沿淮地区广泛栽培的41个菜用大豆品种的12个农艺性状指标进行测定，采用聚类分析和主成分分析法对其农艺性状进行分析，以探明菜用大豆农艺性状评价的主要因素，旨在为菜用大豆品种的科学评价及筛选提供理论依据。

1 材料与方法

1.1 试验材料

本试验材料收集于江淮地区广泛种植的菜用大豆品种，其中产自辽宁省的有25份，其余省份16份，具体品种名称、编号及来源详见表1。试验材料于2017年4月3日种植于安徽科技学院种植科技园，大田内每品种直播种植一小区(小区面积6.3 m2)，行株距30 cm×25 cm。田间管理按本地的生产水平进行，灌溉、施肥与病虫害防治等保持与所选材料一致。

1.2 测定项目

测定项目有生育期(X1)、株高(X2)、茎粗(X3)、倒伏性(X4)、主茎节数(X5)、结荚高度(X6)、单株荚数(X7)、单株有效荚数(X8)、荚长(X9)、荚宽(X10)、单荚粒数(X11)、鲜豆百粒质量(X12)等12个农艺性状，具体考查方法和标准参照《大豆种质资源描述规范和数据标准》[7]。各品种选取有代表性的10株进行测定，试验数据取10株平均值作为该性状的平均值。

1.3 数据统计方法

对所测得的数据采用Excel 2010软件进行处理，应用DPS 7.05数据处理软件对菜用大豆指标进行聚类分析和主成分分析。采用欧式距离和类平均法聚类，依据方差累计贡献率达到85%以上提取主成分[8]，以各主成分对应的方差相对贡献率作为权重，构建菜用大豆的主成分农艺性状评价函数，由此确定各品种间的重要农艺性状指标[9-10]。

2 结果与分析

2.1 农艺性状的变异分析

各品种的农艺性状特征变异情况见表2，由表2可知，通过比较分析，各农艺性状均存在较大的变异，结荚高度变异系数(31.41%)最大，茎粗、株高、倒伏性、单株荚数、单株有效荚数的变异系数集中分布在11.74%～23.82%，而荚长和荚宽的变异系数最小(分别为6.69%和6.45%)。这说明，各菜用大豆栽培品种中，结荚高度、茎粗、株高、倒伏性、单株荚数、单株有效荚数具有一定的选择潜力，可通过良种选配和改善栽培等措施来提高。

2.2 农艺性状的相关性分析

菜用大豆农艺性状间相关系数见表3，结果表明，各农艺性状之间存在一定的相关性，除生育期与株高、单株荚数、荚长，荚长与荚宽间的相关系数大于0.5外，其他指标间的相关系数的绝对值均小于0.5，表明各性状之间的相关性较弱。

表2供试菜用大豆品种资源简介

Table1The introduction of the tested cultivar germplasm of vegetable soybean

品种代号Code品种名称Name种子来源Source品种代号Code品种名称Name种子来源SourceYXS优系三Youxisan辽宁省铁岭市Tieling,LiaoningCJMDW超级毛豆王Chaojimaodouwang山东省济南市Jinan,ShandongZSDJ早熟大荚Zaoshudajia辽宁省铁岭市Tieling,LiaoningXS1鲜食1号XianshiNo.1河北省沧州市Canzhou,Hebei95-195-1黑龙江省尚志市Shangzhi,HeilongjiangSX15沈鲜15号ShenxianNo.15辽宁省沈阳市Shengyang,LiaoningQJN俏江南Qiaojiangnan安徽省合肥市Hefei,AnhuiLX辽鲜Liaoxian黑龙江省尚志市Shangzhi,HeilongjiangTZW特早王Tezaowang江苏省南京市Nanjing,JiangsuLX9807理想9807Lixiang9807辽宁省开原市Kaiyuan,LiaoningLX1605领鲜1605Lingxian1605辽宁省铁岭市Tieling,LiaoningSQ3苏青3号SuqingNo.3辽宁省开原市Kaiyuan,LiaoningCPTZ春棚特早Chunpengtezao辽宁省开原市Kaiyuan,LiaoningSQ1苏青1号SuqingNo.1辽宁省开原市Kaiyuan,LiaoningSZ2苏早2号SuzaoNo.2江苏省南京市Nanjing,JiangsuLXM-7理想M-7LixiangM-7辽宁省开原市Kaiyuan,LiaoningLZTZW绿洲特早王Lvzhoutezaowang辽宁省开原市Kaiyuan,LiaoningXX1新鲜1号XinxianNo.1辽宁省开原市Kaiyuan,LiaoningLLZ粒粒早Lilizao辽宁省开原市Kaiyuan,LiaoningKXXLW开心夏绿王Kaixinxialvwang辽宁省开原市Kaiyuan,LiaoningSY3苏育3号SuyuNo.3江苏省南京市Nanjing,JiangsuJPLQXF极品龙泉先锋Jipinlongquanxianfeng辽宁省开原市Kaiyuan,LiaoningMJZS美极早生Meijizaosheng江苏省南京市Nanjing,JiangsuXFTZ新福特早Xinfutezao辽宁省沈阳市Shengyang,Liaoning75-375-3黑龙江省尚志市Shangzhi,HeilongjiangDLZ大粒早Dalizao黑龙江省尚志市Shangzhi,HeilongjiangCFZ1春丰早Chunfengzao江西省新余市Xinyu,JiangxiWD15皖豆15Wandou15安徽省合肥市Hefei,AnhuiJX6京鲜6号JingxianNo.6辽宁省新民市Xinmin,LiaoningD292292黑龙江省尚志市Shangzhi,HeilongjiangSX3沈鲜3号ShenxianNo.3辽宁省沈阳市Shengyang,LiaoningRBQ日本青Ribenqing辽宁省开原市Kaiyuan,LiaoningSX11沈鲜11号ShenxianNo.11辽宁省沈阳市Shengyang,LiaoningSK1苏奎1号SukuiNo.1辽宁省铁岭市Tieling,Liaoning95ZZ95至尊Zhizun95辽宁省铁岭市Tieling,LiaoningFQ6福气6Fuqi6辽宁省铁岭市Tieling,LiaoningSN8首农8号ShounongNo.3辽宁省铁岭市Tieling,LiaoningLQSJ龙泉神剑Longquanshengjian辽宁省开原市Kaiyuan,LiaoningCFZ2春丰早Chunfengzao安徽省合肥市Hefei,AnhuiJX3H京鲜3号JingxianNo.3辽宁省沈阳市Shengyang,LiaoningLBS绿宝石Lvbaoshi江西省新余市Xinyu,Jiangxi

其中生育期与株高、单株荚数、荚长和荚宽呈极显著正相关，与茎粗、单株有效荚数和鲜豆百粒质量呈显著正相关；株高与荚长呈极显著正相关，与单株荚数、鲜豆百粒重呈显著正相关；茎粗与倒伏性、荚长、荚宽呈显著正相关；主茎节数与荚宽呈显著正相关；单株荚数与单株有效荚数呈极显著正相关，相关系数达0.90；荚长与荚宽呈极显著正相关。其他农艺性状间相关性不明显。

2.3 农艺性状的聚类分析

聚类分析结果表明，在欧式距离5.597 3处可将41个菜用大豆品种划分为四大类群(图1)，各类群品种主要的农艺性状详见表4。第Ⅰ类群仅包括SN8一种，占总品种数的2.4%。该类型的主要特征是生育期较短(86 d)，植株较矮小(平均株高46.20 cm)，茎秆粗壮、不易倒伏，主茎节数较多，结荚高度较高(5.66 cm)，荚较短，鲜豆百粒质量较大(87.6 g)，属直立短粗型、抗倒伏性较好的早熟大粒品种，可以作为选育大粒的早熟品种资源。第Ⅱ类群包括LX1605和LX9807等6个品种，占总品种数的14.6%。该类群的主要特征是生育期长(平均为97.67 d)，植株高大(平均株高为60.96 cm)，易倒伏，结荚高度较低，单株荚数与有效荚数均明显高于其他品种，荚较长，鲜豆百粒质量较大(87.00 g)，属多荚大粒晚熟型品种，可作为选育多荚大粒晚熟型菜用大豆品种的中间材料。第Ⅲ类群包括ZSDJ、KXXLW和WD15等3个品种，占总品种数的7.3%。该类群的主要特征是生育期适中(平均为91.83 d)，株高中等(平均株高为53.59 cm)，不易倒伏，结荚高度适中，单株荚数与有效荚数适中，荚较长，鲜豆百粒质量较大(87.30 g)，属多荚大粒中熟型品种，可作为选育多荚大粒中熟型菜用大豆品种的中间材料。第Ⅳ类群包括YXS和SZ2等31个品种，占总品种数的75.6%。该类群的主要特征是生育期短(平均为83.48 d)，株高矮(平均株高44.76 cm)，单株荚数少，荚短小，鲜豆百粒质量较小(77.30 g)，属小荚小粒早熟型品种，可作为选育小粒早熟型菜用大豆品种的中间材料。

表2菜用大豆主要农艺性状特征及变异

Table2The variations and agronomic characters of vegetable soybean

性状Characters平均值Average变幅Range标准差SD变异系数CV/%生育期Growthperiod/d86.1072.00～102.006.837.93株高Plantheight/cm47.8634.08～43.248.3917.53茎粗Stemdiameter/cm10.398.43～13.051.2211.74倒伏性Lodging1.261.00～2.300.3023.82主茎节数Nodesonmainstem3.693.20～4.200.287.48结荚高度Podheight/cm4.110.98～7.701.2931.41单株荚数Podsperplant33.5917.00～55.707.7323.01单株有效荚数Effectivepodsperplant29.1716.30～46.305.7419.68荚长Podlength/cm5.824.84～6.940.396.69荚宽Podwidth/cm1.321.14～1.460.096.45单荚粒数Grainsperpod2.382.10～3.000.177.20鲜豆百粒质量Freshweightof100grains/g78.7163.50～96.307.038.93

表3菜用大豆农艺性状间的相关系数

Table3The correlation coefficient among agronomic characters of vegetable soybean

相关系数Correlationcoefficientx1x2x3x4x5x6x7x8x9x10x11x20.5746**x30.3487*0.0489x4-0.0064-0.06340.3390*x50.15200.05050.3046-0.1096x6-0.01620.25310.0643-0.08950.0910x70.5380**0.3746*0.20580.10600.2779-0.1357x80.3939*0.18680.24100.19060.2494-0.29750.9041**x90.5174**0.4177**0.3081*0.02900.16360.16520.22030.1608x100.4298**0.01950.3742*0.05420.3161*0.03660.23250.24910.5650**x11-0.12870.0780-0.0311-0.0637-0.17880.1269-0.1017-0.09820.0503-0.2596x120.3859*0.3692*0.03050.0682-0.01520.09680.19530.15250.07140.1066-0.0468

*和**分别表示显著水平为0.05和0.01。

* and ** meant significant difference at 0. 05 and 0.01 level， respectively.

2.3 农艺性状的主成分分析

2.3.1 特征值与贡献率

主成分的特征值与贡献率是选择主成分的依据[11]。对不同品种菜用大豆的农艺性状进行主成分分析，将原农艺性状指标转化为12个主成分信息，结果详见表5。由表5得知，第1个主成分(生育期)的特征值为3.440 5，贡献率为28.670 5%，代表全部农艺性状信息的28.670 5%，是最重要的主成分；第2个主成分(株高)的特征值是1.699 8，贡献率为14.164 6%；第3个主成分(茎粗)的特征值为1.499 4，贡献率是12.494 9%；第4个主成分(倒伏性)的特征值是1.187 0，贡献率是9.891 7%；第5个主成分(主茎节数)的特征值是1.037 0，贡献率为8.641 9；第6个主成分(结荚高度)的特征值是0.915 2，贡献率为7.626 4%；第7个主成分(单株荚数)的特征值是0.664 6，贡献率为5.538 1%。其余主成分的贡献率呈逐渐下降趋势，前7个主成分的累计贡献率为87.028 1%，已基本包含了菜用大豆农艺性状的全部信息，可将前7个主成分作为主成分分析的主要内容。

表4菜用大豆4大类群性状差异

Table4Character differences of four type vegetable soybean germplasm

类群Groups生育期Growthperiod/d株高Plantheight/cm茎粗Stemdiameter/cm倒伏性Lodging主茎节数Nodesonmainstem结荚高度Podheight/cm单株荚数Podsperplant单株有效荚数Effectivepodsperplant荚长Podlength/cm荚宽Podwidth/cm单荚粒数Grainsperpod鲜豆百粒质量Freshweightof100grains/gⅠ86.0046.2011.112.304.005.6634.7032.005.841.462.1087.60Ⅱ97.6760.9610.801.203.803.9945.5035.206.151.352.3087.00Ⅲ91.8353.5910.961.803.904.8240.1033.606.001.412.2087.30Ⅳ83.4844.7610.201.303.603.8131.3028.205.701.302.4077.30

图1 四十一个菜用大豆品种农艺性状的聚类分析图Fig.1 Agronomic characters clustering analysis of 41 germplasm resources of vegetable soybean

表5菜用大豆主要农艺性状的特征值与贡献率

Table5Eigenvalue and contribution rate of agronomic characters of vegetable soybean

主成分Principlecomponent特征值Eigenvalue贡献率Contributionrate/%累计百分率Cumulativepercentage/%13.440528.670528.670521.699814.164642.835131.499412.494955.330041.18709.891765.221751.03708.641973.863760.91527.626481.490070.66465.538187.028180.57684.807191.835290.46353.862195.6973100.25402.116397.8136110.20581.715199.5287120.05660.4713100.0000

2.3.2 主成分函数

根据累积贡献率85%的标准，选取前7个主成分及特征向量列于表6，特征向量的大小表示供试品种中各农艺性状对主成分的贡献大小。前7个主成分的函数表达式分别如下：

y1=0.440 9x1+0.167 3x2-0.103 7x3-0.015 2x4-0.146 6x5-0.196 6x6+0.061 1x7-0.325 9x8-0.220 2x9+0.566 9x10-0.440 7x11+0.180 1x12；

(1)

y2=0.292 0x1+0.436 1x2-0.317 3x3+0.004 2x4+0.013 6x5+0.049 9x6-0.311 1x7-0.333 8x8+0.255 4x9-0.012 5x10+0.587 2x11+0.053 6x12；

(2)

y3=0.277 5x1-0.095 5x2+0.390 5x3+0.379 8x4+0.137 3x5+0.196 0x6+0.219 7x7-0.461 5x8-0.416 5x9-0.310 0x10+0.156 1x11-0.068 7x12；

(3)

y4=0.086 8x1-0.269 2x2+0.066 9x3+0.736 6x4-0.152 1x5+0.178 8x6-0.263 8x7+0.114 1x8+0.400 4x9+0.265 2x10-0.067 1x11-0.003 8x12；

(4)

y5=0.221 4x1-0.105 0x2+0.323 3x3-0.415 9x4+0.233 1x5+0.492 7x6+0.197 0x7-0.147 9x8+0.521 8x9+0.140 9x10-0.108 0x11+0.029 9x12；

(5)

y6=0.001 2x1+0.524 0x2+0.209 2x3+0.036 5x4+0.089 1x5+0.512 5x6-0.334 6x7+0.382 3x8-0.337 5x9+0.060 2x10-0.144 9x11+0.114 5x12；

(6)

y7=0.415 9x1-0.234 0x2-0.320 8x3-0.104 4x4+0.220 0x5+0.108 4x6-0.178 2x7+0.189 7x8-0.165 0x9+0.012 5x10-0.070 1x11-0.705 9x12。

(7)

从函数表达式可以看出，主成分1中，系数值按大小排列，荚宽(x10)、生育期(x1)，这两个性状为系数值较大且为正值，而单株有效荚数(x8)、单荚粒数(x11)均为较大的负系数值，表明第一主成分主要反映的是菜用大豆的荚粒性状特征；主成分2中，正系数值较大的有4个，按大小排列为单荚粒数(x11)、株高(x2)、生育期(x1)、荚长(x9)，负系数值较大的有3个，分别是茎粗(x3)、单株荚数(x7)、单株有效荚数(x8)，表明第二主成分主要反映的是菜用大豆的株型豆荚性状特征；主成分3中，茎粗(x3)的正值最大，倒伏性(x4)也有相对较大的正系数值，而单株有效荚数(x8)、荚长(x9)、荚宽(x10)均为较大的负系数值，因此第三主成分所反映的是菜用大豆的茎秆性状特征；主成分4中，倒伏性(x4)的正系数值最大(0.736 6)，荚长(x9)的正系数值次之，株高(x2)、单株荚数(x7)、主茎节数(x5)的负系数值较大，所反映的是菜用大豆的倒伏性状特征；第5主成分中，荚长(x9)、结荚高度(x6)、茎粗(x3)的正系数值较大，倒伏性(x4)的负系数值最大(-0.415 9)，所反映的是菜用大豆的结荚性状特征；第6主成分中，株高(x2)、结荚高度(x6)、单株有效荚数(x8)、茎粗(x3)的正系数值较大，而单株荚数(x7)、荚长(x9)的负系数值较大，所反映的是菜用大豆的株型性状特征；第7主成分中，生育期(x1)的系数值最大，且为正值(0.415 9)，鲜豆百粒质量(x12)的负系数值也最大(-0.705 9)，表明该主成分反映的是菜用大豆的生育期性状特征。

2.3.3 主成分值分析

根据供试各菜用大豆品种的主成分值可得，第1主成分值的最小值为-3.2107(TZW)， 最大值为4.793 9(LX9807)；第2主成分值的最小值为-2.941 2(SX11)，最大值为3.657 5(WD15)；第3主成分值的最小值为-2.2120(LX1605)，最大值为2.567 5(ZSDJ)；第4主成分值的最小值为-1.614 6(FQ6)，最大值为2.592 1(LX1605)；第5主成分值的最小值为-2.017 8(SN8)，最大值为2.741 9(ZSDJ)；第6主成分值的最小值为-1.792 7(JX6)，最大值为2.222 8(XFTZ)；第7主成分值的最小值为-1.746 2(WD15)，最大值1.644 0(SK1)。这些结果与实际观察到的农艺性状基本相符。通过主成分分析，将12个农艺性状转化为7个相对独立的特征指标，降低了菜用大豆品种性状的分析难度，具有很强的表观性，由于综合信息量较大，因而代表性也很强，增加各品种性状分析的可信度[12]。

表6供试菜用大豆品种的规格化特征向量

Table6Normalized eigenvector of vegetable soybean

主成分PrincipleComponentx1x2x3x4x5x6x7x8x9x10x11x1210.44090.1673-0.1037-0.0152-0.1466-0.19660.0611-0.3259-0.22020.5669-0.44070.180120.29200.4361-0.31730.00420.01360.0499-0.3111-0.33380.2554-0.01250.58720.053630.2775-0.09550.39050.37980.13730.19600.2197-0.4615-0.4165-0.31000.1561-0.068740.0868-0.26920.06690.7366-0.15210.1788-0.26380.11410.40040.2652-0.0671-0.003850.2214-0.10500.3233-0.41590.23310.49270.1970-0.14790.52180.1409-0.10800.029960.00120.52400.20920.03650.08910.5125-0.33460.3823-0.33750.0602-0.14490.114570.4159-0.2340-0.3208-0.10440.22000.1084-0.17820.1897-0.16500.0125-0.0701-0.7059

3 小结与讨论

3.1 选择农艺性状变异系数较高的品种

通过研究分析，结荚高度的变异系数(31.41%)最大，茎粗、株高、倒伏性、单株荚数、单株有效荚数的变异系数同样较大，变异幅度在11.74%～23.82%，而荚长和荚宽的变异系数最小(分别为6.69%和6.45%)。因此，进行品种选育时，首先要考虑对结荚高度、茎粗、株高、倒伏性、单株荚数、单株有效荚数等性状的选择，因其具有一定的选择潜力，可通过良种选配和改善栽培等措施来提高，此外，荚长和荚宽也有一定的改进余地，这样才能比较容易获得综合性优良菜用大豆品种。

3.2 选配亲本组合时应在不同的类群

聚类分析的依据是遗传距离，遗传距离大的品种分为不同的类，遗传距离小的品种分在同一类[13]。遗传距离的计算依赖于各农艺性状的表现型，而表现型受基因型和外部环境条件的共同影响，只有在外部环境一致的情况下，聚类分析结果才能更准确。分析结果表明，来自同一地区亲缘关系较近的品种大多数被分在同一类群，绝大部分为辽、黑地区品种，只有少数被分在不同的类群，这与罗珊等[14]的研究结果一致。本试验41个品种被分为4大类群，各类型的农艺性状差异明显，其中第Ⅰ类群和第Ⅱ类群所表现的性状特征较为优良，可用于指导菜用大豆的亲本选育实践，具有一定的实际意义。另外，将农艺性状聚类分析结果与用分子标记方法的分类结果结合起来，综合评价菜用大豆种质资源实践效果可能会更好。

3.3 集中考察综合性状因子来提高育种效率

主成分分析法由于各主成分是一个相对独立的指标体系，各主成分之间不存在相关性，并且数据直观、易于分析[15]。从主成分分析来看，其性状信息主要集中在前7个主成分，累计贡献率达87.0281%，其所表达的综合信息基本可以用来表达全部农艺性状的全部信息，育种工作往往不可能同时考虑所有性状，做到对逐个性状的改良，所以在菜用大豆育种中，可以尝试集中考察7个综合性状因子，重点关注影响菜豆的生育期、株高和茎粗等主要性状，这对菜用大豆的生产应用尤为重要，忽略其次要性状，以提高菜用大豆育种效率，此结果与石惠等[16]的研究结果一致。

3.4 充分利用性状间的相关性提高选择效率

各品种中生育期与株高、单株荚数、荚长和荚宽呈极显著正相关，因此可根据安徽不同地区的具体情况，选择生育期合适的品种，在生育期许可的范围内，可适当选择生育期较长的菜用大豆品种，有利于提高单株产量；品种生育期越长，营养生长期也就越长，荚大且荚数多，为高产奠定一定的基础，即生育期长、结荚多、籽粒多、株型高大的多为高产类型。此结果与梁镇林等[17]的研究结果相似。株高的增加会引起结荚高度、单株荚数的相应增加，鲜豆百粒重也会有明显的增加，此结果与陈学珍等[18]的研究结果不太一致，可能与当地种植条件有关。茎粗、主茎节数均与荚长荚宽正相关，增加茎粗的同时茎秆越抗倒伏，但增加主茎节数可能会造成鲜豆百粒重的下降，此结果与杨守臻等[19]的研究结果基本相同。各性状间相互影响、相互制约，因此想要选育高产优质的菜用大豆品种，须考虑各性状特征间的相互协调。

[1] 李鹏. 菜用大豆品种遗传多样性的聚类分析[J]. 种子世界，2017(1): 16-19.

LI P. Clustering analysis of genetic diversity of vegetable soybean varieties[J].SeedWorld，2017(1):16-19. (in Chinese)

[2] 盖钧镒，王明军，陈长之.中国毛豆生产的历史渊源与发展[J]. 大豆科学，2002，21(1): 7-13.

GAI J Y， WANG M J， CHEN C Z. Historical origin and development of Maodou production in china[J].SoybeanScience， 2002，21(1): 7-13. (in Chinese with English abstract)

[3] 张德纯.蔬菜史话·菜用大豆[J]. 中国蔬菜，2009(17):40.

ZHANG D C. Vegetables soybeans[J].ChineseVegetables， 2009(17): 40. (in Chinese)

[4] 赵璇，王玉岭，李占军，等. 菜用大豆的发展现状及展望[J]. 河北农业科学，2012，16(7):39-41.

ZHAO X， WANG Y L， LI Z J， et al. Development status and prospect of vegetable soybean[J].HebeiAgriculturalSciences， 2012，16(7): 39-41. (in Chinese with English abstract)

[5] 顾卫红，郑洪建. 菜用大豆的国际需求及科研生产动态(综述)[J]. 上海农业学报， 2002， 18(2):45-48.

GU W H， ZHENG H J. International demand and quality of vegetable soybean research and production dynamics (review)[J].JournalofShanghaiAgricultural， 2002， 18 (2): 45-48. (in Chinese with English abstract)

[6] 农业部关于印发《全国种植业结构调整规划(2016-2020年)》的通知.[2016-04-28]. http://www. moa.gov.cn/ zwllm/tzgg/tz/201604/t20160428_5110638.htm.

[7] 常汝镇，孙建英，邱丽娟. 中国大豆种质资源研究进展[J]. 作物杂志， 1998(3): 7-9.

CHANG R Z， SUN J Y， QIU L J. Research progress of soybean germplasm resource in China[J].Crops， 1998(3): 7-9.(in Chinese)

[8] 刘国伟，任艳云，高园园，等. 大蒜品种资源生育和产量性状主成分聚类分析及评价[J]. 浙江农业学报，2016， 28(9): 1508-1513.

LIU G W， REN Y Y， GAO Y Y， et al. Cluster analysis and evaluation of agronomic characters in garlic variety resources based on principal components[J].ActaAgriculturaeZhejiangensis， 2016， 28(9): 1508-1513.(in Chinese with English abstract)

[9] 徐臣善，高东升. 基于主成分分析的设施桃果实品质综合评价[J]. 食品工业科技， 2014， 35(23): 84-88，94.

XU C S， GAO D S. Comprehensive evaluation on fruit quality of peach cultivars in greenhouse based on prin-cipa component analysis[J].Science&TechnologyofFoodIndustry， 2014， 35(23): 84-88.94). (in Chinese with English abstract)

[10] 鲍江峰，夏仁学，邓秀新，等. 用主成分分析法选择纽荷尔脐橙品质的评价因素[J]. 华中农业大学报， 2004，23(6): 663-666.

BAO J F， XIA R X， DENG X X， et al. The quality evaluation factors selection of Newhall Orange by the principal component analysis[J].JournalofHuazhongAgriculturalUniversity， 2004，23(6): 663-666.(in Chinese with English abstract)

[11] 周丽艳，郭振清，马玉玲，等. 春小麦品种农艺性状的主成分分析与聚类分析[J]. 麦类作物学报， 2011， 31(6): 1057-1062.

ZHOU L Y， GUO Z Q， MA Y L，et al. Principal component and cluster analysis of different spring wheat cultivars based on agronomic traits[J].JournalofTriticeaeCrops， 2011， 31(6): 1057-1062.(in Chinese with English abstract)

[12] 张运锋，黄光和，李希国，等. 普通菜豆农艺性状相关性分析及主成分分析[J]. 安徽农学通报， 2008，14(7): 97-98,137.

ZHANG Y F，HUANG G H，LI X G，et al. Correlation analysis and principal component analysis of agronomic traits of common[J].AuhuiAgriculturalScienceBulletin， 2008，14(7): 97-98,137.(in Chinese with English abstract)

[13] 李卫华，伏河山，徐红军，等. 42个春小麦品种数量性状的聚类分析[J]. 新疆农业科学， 2006，43(5): 357-361.

LI W H，FU H S，XU H J，et al. Cluster analysis of quantitative traits in 42 spring wheat varieties[J].XinjiangAgriculturalSciences， 2006， 43(5):357-361.(in Chinese with English abstract)

[14] 罗珊，康玉凡，濮绍京，等. 黑河地区55份大豆品种资源农艺性状和营养成分的聚类分析[J]. 大豆科学， 2009， 28(3): 421-425.

LUO S，KANG Y F，PU S J，et al. Cluster analysis of agronomic characters and nutrition composition of 55 soybean accessions from Heihe area[J].SoybeanScience， 2009， 28(3): 421-425. (in Chinese with English abstract)

[15] 李桂萍，张根生，巴青松，等. 杂种小麦品质性状的性状相关和主成分分析[J]. 浙江农业学报， 2016， 28(9): 1447-1453.

LI G P，ZHANG G S，BA Q S，et al. Correlation analysis and principal component analysis on quality traits in hybrid wheat[J].ActaAgriculturaeZhejiangensis， 2016， 28(9): 1447-1453.(in Chinese with English abstract)

[16] 石惠，许海涛.大豆主要农艺性状的遗传变异及相关性和主成分分析[J]. 黑龙江农业科学， 2008(2): 29-32.

SHI H，XU H T. Genetic variation correlation and principal component analysis on major agronomic trait of soybean[J].HeilongjiangAgriculturalScience， 2008(2): 29-32. (in Chinese with English abstract)

[17] 梁镇林，梁颖. 间作大豆产量与主要经济性状的相关及选择[J]. 大豆科学，1997，16(1):54-59.

LIANG Z L，LIANG Y. Correlation and selection of intercropped soybean yield with main economic traits[J].SoybeanScience， 1997，16(1): 54-59.(in Chinese with English abstract)

[18] 陈学珍，谢皓，郑晓宇，等. 菜用大豆的农艺及品质性状评价与相关性分析[J]. 北京农学院学报，2005，20(1):23-26.

CHEN X Z，XIE H，ZHENG X Y，et al. Relation analysis and evaluation on agronomic and quality characters in vegetable voybean[J].JournalofBeijingAgriculturalCollege， 2005，20(1): 23-26.(in Chinese with English abstract)

[19] 杨守臻，孙祖东，陈怀珠，等. 毛豆品种的农艺性状鉴定及相关性分析[J]. 中国农学通报，2005，21(3):169-171.

YANG S Z，SUN Z D，CHEN H Z，et al. The agronomic trait and correlation analysis of soybean varieties[J].ChineseAgriculturalScienceBulletin， 2005，21(3): 169-171.(in Chinese with English abstract)