基于主成分分析和聚类分析的大豆高产群体性状研究.pdf

江西农业学报 2023，35（06）：2429ActaAgriculturaeJiangxiDOI:10.19386/ki.jxnyxb.2023.06.004基于主成分分析和聚类分析的大豆高产群体性状研究孙 影，叶 录*，王凤梧，廉培勇，郑成忠，梅 雪，张子臻，王千军，徐振朋，王小明，黄文娟（乌兰察布市农林科学研究所，内蒙古 乌兰察布 012000）摘 要:综合分析了20份大豆品种的产量性状和农艺性状，得出了影响大豆产量形成的关键因子和高产大豆群体的性状特征，结果表明：大豆的单株无效荚数和有效分枝数的变异系数大，分别为0.83、0.77，说明这2个指标受环境的影响大，在遗传性状改良时可作为主要改良的目标性状；大豆产量与株高、主茎节数、底荚高度、单株有效荚数、单株总粒数呈显著正相关；主成分分析得出单株有效荚数、单株总粒数、株高、百粒重、主茎节数的载荷量较高，分别为0.774、0.766、0.705、0.863和0.661，可作为产量的主要影响因子；高产群体的农艺和产量性状表现为株高高、主茎节数多、底荚高度高、粒重大、单株总粒数和有效荚数多，说明在引进或者选育高产品种时，应首要考虑对这些指标进行科学、合理的选择。关键词:大豆；产量性状；农艺性状；高产群体；主成分分析；聚类分析 中图分类号：S565.1 文献标志码：A 文章编号：1001-8581（2023）06-0024-06Study on High-yield Population Traits of Soybean Based on Principal Component Analysis and Cluster Analysis SUNYing,YELu*,WANGFeng-wu,LIANPei-yong,ZHENGCheng-zhong,MEIXue,ZHANGZi-zhen,WANGQian-jun,XUZhen-peng,WANGXiao-ming,HUANGWen-juan(UlanqabInstituteofAgriculturalandForestryScience,Ulanqab012000,China)Abstract:Throughcomprehensiveanalysisoftheyieldandagronomictraitsof20soybeanvarieties,thekeyfactorsaffectingsoybeanyieldformationandthetraitcharacteristicsofhigh-yieldsoybeanpopulationsweresummarized.Theresultsshowedthatthecoefficientofvariationoftheineffectivepodsperplantandtheeffectivebranchesofsoybeanwaslarge,whichwas0.83and0.77respectively,indicatingthatthesetwoindicatorsweregreatlyaffectedbytheenvironment,andcouldbeusedasthemainimprovedtraitsintheimprovementofgenetictraits.Theyieldofsoybeanwassignificantlypositivelycorrelatedwithplantheight,nodenumberofmainstem,heightofbottompod,effectivepodnumberperplantandtotalgrainnumberperplant.Theprincipalcomponentanalysisshowedthattheloadofeffectivepodsperplant,totalgrainsperplant,plantheightand100-grainweightwere0.774,0.766,0.705,0.863and0.661respectively,whichcouldbeusedasthemainfactorsaffectingtheyield.Theagronomicandyieldtraitsofhigh-yieldpopulationaretallplant,morenodesofmainstem,highbottompodheight,largegrainweight,moretotalgrainsperplantandmoreeffectivepods,indicatingthattheseindicatorsshouldbeconsideredfirstwhenintroducingorbreedinghigh-yieldvarieties.Key words:Soybean;Yieldtrait;Agronomictrait;High-yieldpopulation;Principalcomponentanalysis;Clusteranalysis大豆作为重要的蛋白质和油脂来源之一，在我国已有2000多年的栽培历史1。大豆产量性状和相关农艺性状多为数量性状，容易受环境因素的影响，而且农艺性状之间还会相互影响，导致分析产量与相关农艺性状之间的关系比较困难。在特定的生态条件下，研究特定性状与产量的关系，可以直观地筛选出高产的指标性状2。大豆的产量受多个基因控制，在育种上很难直接选择某个性状，因此，通过选择与产量相关的农艺性状已成为作物高产育种的一个重要途径3。主成分分析、聚类分析与其他分析方法联合使用，在品种评价、分类、选育、建立理想的高产株型方面均有独特的优势4-7。目前，关于大豆的研究多集中在高产栽培技术收稿日期：2023-04-28基金项目：乌兰察布市农林科学研究所青年项目（2022WNLQN01）。作者简介：孙影（1989），江苏盱眙人，农艺师，硕士，从事燕麦等杂粮杂豆栽培与育种研究。*通信作者：叶录。6 期孙影等：基于主成分分析和聚类分析的大豆高产群体性状研究25方面，对大豆品种主要农艺性状与产量之间的关系及其变化规律已有一些研究8-10，但采用主成分分析与聚类分析相结合的方法对不同品种相关农艺性状进行量化分析的研究鲜见报道。本研究采用多元统计分析方法，通过对20种大豆品种的产量和农艺性状进行分析，探讨了它们之间的量化关系，并对各参试品种的性状进行了聚类分析和综合评价，旨在明确大豆高产群体的特征，为内蒙古乌兰察布市大豆的生产、引种及良种选育提供有价值的参考。1 材料与方法1.1 试验地点与供试材料本试验于2022年在乌兰察布市农林科学研究所平地泉基地进行，地理位置113.1196 E、40.9232 N，海拔1342m，生育期平均温度16.92，降水量278.9mm，日照时数1084.7h。土壤质地为砂壤土，地力良好、平衡，前茬作物为马铃薯。供试大豆品种：内豆4号、蒙豆343、蒙豆160、蒙豆170、蒙豆42、蒙豆43、蒙豆46、蒙豆48、蒙豆49、蒙豆50、蒙豆15、登科5号、登科1号、蒙豆53、蒙豆54、蒙豆57、蒙豆58、蒙豆1137、蒙豆912、东农690共20个大豆品种，均由呼伦贝尔市农牧科学研究所提供。1.2 试验设计2022年5月20日播种，以品种为处理，采取随机区组排列，设置3个重复，合计60个小区，每个小区面积为16.5m2（5.5m3.0m），株行距18cm25cm，保护行2行，小区间观察道宽0.5m，小区中间过道宽1.0m。1.3 测定内容与方法记录各参试大豆品种的生育期、收获期；在每个小区随机挑选大豆植株10株，并分别测定其株高、主茎节数、结荚高度、分枝数等农艺性状，以及单株有效荚数、单株无效荚数、单株总粒数、百粒重等产量构成性状。每个小区收割20穴，脱粒，晒干，测定20株的总产量。同时，每个小区全部收割，混合脱粒，测定每个小区的实际产量。1.4 统计分析运用Excel2010和SPSS19.0软件对试验数据进行处理和统计分析。2 结果与分析2.1 生育期由表1可以看出，内豆4号的生育期最短，为92d；蒙豆912的生育期最长，为122d。除此之外，其他参试大豆品种的生育期在108118d之间。总体来说，除内豆4号外，参试大豆品种的生育期较为接近，在内蒙古乌兰察布地区种植均能够正常成熟。表1 参试大豆品种的生育期情况品种生育期/d收获期/(月/日)品种生育期/d收获期/(月/日)内豆4号928/20蒙豆501139/10蒙豆531089/5登科1号1139/10蒙豆541089/5蒙豆581139/10蒙豆3431099/6东农6901139/10蒙豆1601099/6蒙豆481149/11登科5号1099/6蒙豆151149/11蒙豆571099/6蒙豆421159/12蒙豆1701109/7蒙豆431189/15蒙豆461109/7蒙豆11371189/15蒙豆491119/8蒙豆9121229/192.2 产量和农艺性状测定参试品种的单株有效荚数、无效荚数、总粒数、有效分枝数，以及百粒重、产量、株高、主茎节数、底荚高度等指标，结果(表2)表明：不同品种的各个农艺性状和产量性状之间的差异较大，变异系数最大的是单株无效荚数，为0.83；其次是单株有效分枝数，为0.77，这说明无效荚数和有效分枝数会随外界环境的改变而发生较大的变化。在大豆的实际生产中，可以通过调控农艺和产量性状来达到提高产量的目的。主茎节数的变异系数最小，为0.12；其他产量和农艺性状的变异系数较小，集中在0.130.23之间，说明这些指标受遗传因素的影响较大，有一定的稳定性。2.3 产量和农艺性状之间的相关性分析对参试大豆品种的农艺性状和产量性状进行Person相关性分析，结果如表3所示，由此可以看出：小区产量与株高、主茎节数、底荚高度、单株有效荚数、单株总粒数呈显著正相关，与20株产量呈极显著正相关；20株产量与株高、单株有效荚数呈显著正相关，与单株总粒数呈极显著正相关；百粒重与单株有效荚数、单株总粒数呈显著负相关；单株总粒数与单株有效荚数呈极显著正相关；单株无效荚数与底荚高度呈极显著负相关；株高与单株有效荚数、主茎节数呈显著正相关。江 西 农 业 学 报35 卷26表2 参试大豆品种主要产量和农艺性状的表现品种产量性状农艺性状单株有效荚数单株无效荚数单株总粒数百粒重/g20株产量/g小区产量/kg株高/cm主茎节数底荚高度/cm单株有效分枝数内豆4号17.81.343.018.56.83.565.09.37.50.3蒙豆34324.00.553.520.98.94.256.111.710.01.0蒙豆16039.00.881.220.813.75.995.013.38.91.0蒙豆17027.04.061.318.67.83.742.410.03.51.0蒙豆4232.20.576.519.413.15.668.912.310.01.3蒙豆4334.50.883.017.914.46.485.012.710.20.0蒙豆4622.00.556.320.47.53.668.511.38.20.3蒙豆4825.51.347.321.88.64.272.111.713.30.0蒙豆4926.00.881.818.711.95.359.211.311.71.7蒙豆5028.31.073.520.212.85.576.811.711.31.3蒙豆1524.00.855.321.89.55.166.713.314.52.0登科5号25.01.565.521.112.35.467.115.710.50.7登科1号24.02.059.319.510.05.283.014.010.71.0蒙豆5320.80.354.818.99.34.756.711.78.81.0蒙豆5418.51.350.320.68.33.969.310.79.81.0蒙豆5722.00.064.020.912.15.465.510.712.50.3蒙豆5820.51.858.321.011.55.267.510.310.70.0蒙豆113732.50.573.320.414.06.278.912.312.20.3蒙豆91225.30.359.021.39.65.188.712.311.50.3东农69038.30.8101.310.38.64.076.912.710.82.0平均值26.361.0464.9319.6510.544.9170.4711.9510.330.83标准差6.080.8714.522.482.370.8812.281.472.330.63变异系数0.230.830.220.130.220.180.170.120.230.77表3 参试大豆品种产量性状、农艺性状间的相关性性状株高主茎节数底荚高度单株有效分枝数单株有效荚数单株无效荚数单株总粒数百粒重20株产量小区产量株高1主茎节数0.474*1底荚高度0.4000.4001单株有效分枝数-0.2000.2000.1001单株有效荚数0.509*0.4000.1000.3001单株无效荚数-0.400-0.200-0.584*0.000-0.1001单株总粒数0.4000.4000.1000.4000.856*-0.2001百粒重0.0000.0000.200-0.400-0.454*0.000-0.635*120株产量0.459*0.4000.300-0.1000.557*-0.3000.582*0.2001小区产量0.535*0.510*0.450*-0.1000.485*-0.3000.480*0.2000.956*1注：*、*分别表示在0.05、0.01水平(双侧)上显著相关。2.4 参试大豆品种的主成分分析对试验数据进行Bartlett球形度检验和KMO检验(表4)，根据检验要求，如果KMO0.8，则说明非常适合进行分析；如果KMO值介于0.70.8之间，则说明比较适合进行分析；如果KMO值介于0.60.7之间，则说明可以进行分析；如果KMO值0.6，则说明不适合进行分析。本试验符合KMO0.6且P0.05的标准，这说明适合进行主成分分析。表4 Bartlett球形检验和 KMO检验结果KMO值Bartlett球形度检验近似卡方dfP值0.617134.053450.0006 期孙影等：基于主成分分析和聚类分析的大豆高产群体性状研究27对参试品种的农艺和产量性状做主成分分析，从表5可以看出：通过主成分分析提取出3个主成分，第1主成分的特征值为4.057，贡献率为40.569%；第2主成分的特征值为2.230，贡献率为22.302%；第3主成分的特征值为1.197，贡献率为11.970%；前3个主成分的累计贡献率为74.841%。表5 参试大豆品种主成分分析结果主成分特征值方差解释率/%累计贡献率/%14.05740.56940.56922.23022.30262.87131.19711.97074.841表6反映出了各主成分的载荷量和共同度，由此可以看出：第1主成分包括小区产量(载荷量0.847)、20株产量(0.838)、单株有效荚数(0.774)、单株总粒数(0.766)、株高(0.705)、主茎节数(0.661)和底荚高度(0.508)；第2主成分包括百粒重(0.863)；第3主成分包括单株无效荚数(0.548)。表6 载荷量与共同度性状载荷量公因子方差(共同度)主成分1主成分2主成分3株高0.7050.2170.0640.548主茎节数0.6610.042-0.1380.458底荚高度0.5080.481-0.5770.823单株有效分枝数0.173-0.605-0.5470.695单株有效荚数0.774-0.4720.2110.867单株无效荚数-0.455-0.3560.5480.635单株总粒数0.766-0.5810.0640.928百粒重-0.1490.8630.1370.78620株产量0.8380.2190.3350.863小区产量0.8470.3220.2480.8822.5 聚类分析对20个品种的10个性状进行系统聚类分析，结果如图1所示，按照聚类结果在最大距离为10时，可将20个品种划分为4个大类：类群包括蒙豆343、蒙豆53、蒙豆46、蒙豆54、内豆4号、蒙豆57、蒙豆58、蒙豆48、蒙豆912；类群包括蒙豆42、蒙豆50、蒙豆49、蒙豆43、蒙豆1137、蒙豆160、登科5号、登科1号、蒙豆15；类群包括东农690；类群包括蒙豆170。11图1 参试大豆品种的聚类分析结果根据聚类分析结果，对4个大类的产量和农艺性状进行系统分析，总结出各大类群的性状特征。结果表明(表7)：类群品种的株高、主茎节数、底荚高度、单株有效荚数、20株产量和小区产量适中，有效分枝数和单株无效荚数少，百粒重高。类群品种的株高、主茎节数、底荚高度、百粒重、20株产量和小区产量高，单株总粒数和单株有效荚数较多，单株有效分枝数适中，单株无效荚数少，该类型是典型的高产类型，具有植株高大，单株总粒数和单株总荚数较多且颗粒大的特点。类群品种的株高、主茎节数、单株有效分枝数、单株有效荚数、单株总粒数高，底荚高度、20株产量和小区产量适中，单株无效荚数和百粒重低，该类型是典型的株高粒小型。类群品种的株高、主茎节数、底荚高度、20株产量和小区产量低，单株有效分枝数、单株有效荚数、单株总粒数、百粒重适中，单株无效荚数多，该类型植株矮小，产量低，容易受环境的影响而产生空瘪粒。表7 参试大豆品种4个类群的性状特征类群株高/cm主茎节数底荚高度/cm单株有效分枝数单株有效荚数单株无效荚数单株总粒数百粒重/g20株产量/g小区产量/kg67.711.110.20.521.80.854.020.5183.74.475.613.011.11.029.50.972.120.0248.45.676.912.710.82.038.30.8101.310.3172.24.042.410.03.51.027.04.061.318.6156.03.7江 西 农 业 学 报35 卷283 讨论3.1 产量、农艺性状的变异性和相关性分析大豆植株性状的变异系数反映了不同大豆品种随着环境的改变而发生改变的可能性，韩秉进等11认为除了百粒重、主茎节数、每荚粒数和经济系数的变异较小，改良空间较小外，大豆的其他性状都有一定的改良空间；4粒荚数、分枝数的变异系数最大，可通过品种的更换或改良、施肥的改善等栽培措施，使其具有较大的改良空间，其也可以作为育种选择的性状。本研究发现，大豆的单株无效荚数和单株有效分枝数的变异系数最大，有很大的改良空间，而百粒重、主茎节数、株高等性状的变异系数小，受遗传因素的影响较大，这与前人的研究结果相似。符小发等12认为单株粒数、单株荚数、生物产量、株高等性状能为大豆产量的积累提供足够的“库”，这些指标与产量呈显著正相关；潘相文等13亦认为在选择大豆高产品种时，应侧重选择这些性状为重要指标。王金陵等14-15认为大豆的株高与单株产量呈正相关，而孟庆喜等16认为大豆的株高与单株产量在肥力高的田块呈负相关，在肥力低的田块则呈正相关。本研究结果表明：与大豆产量呈显著正相关的性状指标有株高、主茎节数、底荚高度、单株有效荚数、单株总粒数，这与前人的研究结果相似。3.2 主成分分析与聚类分析众多学者认为单株粒重是提高产量的关键因素，即单株粒重对产量的贡献最大17-21。符小发等12研究发现不同农艺性状对不同大豆品种单株产量的影响表现为单株荚数单株粒数生物产量有效分枝数百粒重生育期株高。任海龙等22得出大豆单株产量与农艺性状的关联强度表现为单株粒数单株荚数有效分枝数株髙主茎节数生育期百粒重。王彩洁等6认为分枝数、单株荚数是影响大豆产量的关键因素，因而是育种过程中选择的主要目标性状。本研究对20份参试大豆材料进行主成分分析后，得出了3个主成分，3个主成分的累计贡献率为74.841%。第1主成分的特征值为4.057，贡献率为40.569%，说明在本试验中该类主成分对产量形成有重要的影响。除去产量因素，还有百粒重、单株有效荚数、单株总粒数、株高和主茎节数的载荷量较高，分别为0.863、0.774、0.766、0.705和0.661，说明这几个农艺性状对产量的影响较大，因此在育种中应该优先选择具有这类优良性状的大豆材料，并在高产栽培模式创建过程中对这些性状进行重点调控。本研究通过聚类分析发现产量最高的是类群，说明高产大豆品种的特点为株高、主茎节数、底荚高度、百粒重高，单株总粒数和单株有效荚数较多，单株有效分枝数适中，单株无效荚数少。由此可知，选育高产大豆品种时不能一味地追求单株总荚数和总粒数，而是在保证这2项指标的基础上，适当提高株高，增加百粒重。这与陈桂娟等23的研究结果相似。4 结论乌兰察布市气候干燥，昼夜温差大，光能资源极为丰富，日平均日照时数为7.68.4h，可以为大豆产量的积累提供足够的“源”。本试验选取的20个大豆品种均能在本地区正常生长成熟。综合聚类分析、主成分分析和相关性分析结果，大豆高产群体的农艺性状和产量性状表现为株高高、主茎节数多、底荚高度高、粒重大、单株总粒数和有效荚数多，说明在选育高产品种时，要首先考虑这些指标。变异性分析表明：单株无效荚数和单株有效分枝数的变异系数最大，说明在品种遗传改良时这2个性状可作为优先考虑的农艺性状。参考文献：1 JianYY.SituationofsoybeanproductionandresearchinChina J.TropicalAgricultureResearchSeries,1984,17:67-72.2 DonaldCM.Competitiveplants,communalplants,andyieldinwheatcrops M/EvansLTP.Wheatsciencetodayandtomorrow.Cambridge:CambridgeUniversityPress,1981:223-247.3 王海杰,任海龙,林力,等.灰色关联分析在春大豆产量分析中的应用 J.广东农业科学,2013,40(24):11-14.4 郭数进,杨凯敏,霍瑾,等.大豆不同性状与产量的相关性及主成分分析 J.山西农业科学,2015,43(5):505-508.5 赵银月,耿智德,保丽萍,等.云南省大豆地方品种资源的主成分分析及聚类分析 J.湖南农业大学学报:自然科学版,2007(s1):120-122.6 王彩洁,李连华,李伟,等.大豆品种产量与主要性状的主成分分析 J.山东农业科学,2008(1):5-6.7 张玉革,胡绪彬.基于主成分和聚类分析的大豆品种生6 期孙影等：基于主成分分析和聚类分析的大豆高产群体性状研究29物学性状的比较研究 J.大豆科学,2004,23(3):178-183.8 KumudiniS,HumeDJ,ChuG.Geneticimprovementinshortseasonsoybeans:.Drymatteraccumulation,partiti-oning,andleafareaduration J.CropScience,2001,41(2):391-398.9 VoldengHD,CoberER,HumeDJ,etal.Fifty-eightyearsofgeneticimprovementofshort-seasonsoybeancultivarsinCanada J.CropScience,1997,37(2):428-431.10 杨守臻,陈怀珠,李初英,等.木薯主要农艺性状的遗传变异、相关性和主成分分析 J.中国农学通报,2006,22(7):232-234.11 韩秉进,潘相文,金剑,等.大豆农艺及产量性状的主成分分析 J.大豆科学,2008,27(1):67-73.12 符小发,高强,任海龙,等.南繁大豆农艺性状与产量的主成分和聚类分析 J.广东农业科学,2018,45(4):6-13.13 潘相文,李文滨,李艳华,等.主成分分析在大豆抗旱性评价上的应用 J.大豆科学,2006,25(4):379-384,389.14 王金陵.大豆农艺性状的遗传传递规律与大豆的杂交育种 J.中国农业科学,1962(12):1-5.15 马育华,盖钧镒.江淮下游地区大豆地方品种的初步研究:数量性状的遗传变异 J.遗传学报,1979,6(3):331-338.16 孟庆喜,杨庆凯,高凤兰,等.不同生态条件对大豆品种产量性状的影响 J.东北农学院学报,1981(3):15-21.17 NimmakayalaP,AbburiVR,SaminathanT,etal.Genome-widedivergenceandlinkagedisequilibriumana-lysesforCapsicum baccatumrevealedbygenome-anchoredsinglenucleotidepolymorphisms J.FrontiersinPlantScience,2016,7:1646.18 AndersenJR,LubberstedtT.FunctionalmarkersinplantsJ.TrendsinPlantScience,2003,8(11):554-560.19 KimJE,OhSK,LeeJH,etal.Genome-wideSNPcallingusingnextgenerationsequencingdataintomato J.MoleculesandCells,2014,37(1):36-42.20 ZhangT,YouJ,YuGL,etal.Genemappingandcandidategeneanalysisofaberrant-floral spikelet1(afs1)inrice(Oryza sativaL.)J.JournalofIntegrativeAgriculture,2020,19(4):921-930.21 宋健.大豆种皮色相关基因的图位克隆及功能解析D.北京:中国农业科学院,2019:31-45.22 任海龙,徐麟,刘国志,等.新疆大豆在低纬地区单株产量相关性状的灰色分析 J.内蒙古农业大学学报:自然科学版,2014,35(2):23-27.23 陈桂娟,吴国华.夏大豆产量与主要农艺性状相关性分析研究 J.农业科技通讯,2014(9):123-124.（责任编辑：曾小军）