基于PCA方法的夏玉米穗部性狀綜合評(píng)價(jià).doc

基于PCA方法的夏玉米穗部性狀綜合評(píng)價(jià) 趙玉坤1，高根來(lái)1，田瑋瑋2，王向東1，甄勝虎1 （1山西省農(nóng)業(yè)科學(xué)院小麥研究所，山西臨汾041000；2安澤縣農(nóng)業(yè)局，山西安澤042500） 摘要：為了探究不同年際間夏玉米穗部多性狀統(tǒng)計(jì)數(shù)據(jù)的內(nèi)在關(guān)聯(lián)規(guī)律，以晉南地區(qū)品比圃鑒定2個(gè)玉米對(duì)照品種和7個(gè)雜交組合為研究對(duì)象，調(diào)查xxxx年3年間11個(gè)玉米穗部性狀因子的數(shù)據(jù)，通過(guò)主成分分析，提取34個(gè)主成分，評(píng)價(jià)各雜交組合在穗部性狀上的數(shù)量差異。結(jié)果表明，不同年際間提取的各主成分側(cè)重反應(yīng)的穗部性狀分量各有不同，但在分別構(gòu)建主成分評(píng)價(jià)函數(shù)后，主成分評(píng)價(jià)得分表現(xiàn)相似的變化趨勢(shì)；與對(duì)照品種（鄭單958、先玉335）相比，雜交組合PB6和PB7綜合主成分得分較高，表現(xiàn)出優(yōu)異的綜合穗部性狀信息，存在一定的增產(chǎn)潛力。因此，可將主成分分析作為玉米常規(guī)育種的輔助數(shù)據(jù)分析手段，提高育種效率。 關(guān)鍵詞：主成分分析；穗部性狀；評(píng)價(jià)函數(shù) ：S513.01文獻(xiàn)標(biāo)志碼：A論文編號(hào)：xx-0818 基金項(xiàng)目：山西省農(nóng)科院育種工程“優(yōu)質(zhì)多抗中早熟玉米新品種選育”(11YZGC055)。 第一作者簡(jiǎn)介：趙玉坤，男，1984年出生，山西運(yùn)城人，助理研究員，碩士，主要從事玉米遺傳育種及栽培研究。通信地址：041000山西省臨汾市幽并街33號(hào)省農(nóng)科院小麥研究所，TelE-mail：tykzhao163.。 通訊作者：高根來(lái)，男，1956年出生，山西襄汾人，研究員，主要從事玉米遺傳育種及栽培研究。E-mail：gaogenlai_1956126.。 ：xx-08-17，修回日期：xx-10-09。 0引言 玉米穗部性狀與產(chǎn)量密切相關(guān)，通過(guò)玉米穗部性狀間接篩選優(yōu)良品種的效率要遠(yuǎn)高于對(duì)產(chǎn)量的直接篩選1，針對(duì)具體穗部性狀目的明確地遺傳改良，能夠有效提高玉米田間產(chǎn)量2-4，如何借助于有效的數(shù)據(jù)統(tǒng)計(jì)分析手段，評(píng)價(jià)玉米穗部多性狀的綜合優(yōu)劣程度，是玉米育種工作的重要環(huán)節(jié)。本研究宗旨在通過(guò)對(duì)晉南地區(qū)3年間夏玉米穗部性狀數(shù)據(jù)的主成分提取及函數(shù)評(píng)價(jià)分析，揭示穗部多性狀數(shù)據(jù)在不同主成分載荷上的內(nèi)在關(guān)聯(lián)規(guī)律，以主成分函數(shù)模型評(píng)價(jià)不同玉米品種穗部性狀的綜合表現(xiàn)，為玉米高效育種提供理論參考。基于多元統(tǒng)計(jì)分析的主成分分析法(PrincipalComponentAnalysis，PCA)，通過(guò)少數(shù)幾個(gè)彼此不相關(guān)的主成分來(lái)解釋多變量的方差-協(xié)方差結(jié)構(gòu)，將多變量從高維空間映射到低維空間，簡(jiǎn)化多變量數(shù)據(jù)分析，已在模型建立5-6、指標(biāo)評(píng)價(jià)7-8中有較多應(yīng)用。 玉米穗部性狀大多是由微效多基因控制的數(shù)量性狀，受外界復(fù)雜環(huán)境因子的支配影響較大，因此，在穗部性狀研究中多采用一些生物數(shù)理統(tǒng)計(jì)方法，來(lái)完成對(duì)玉米穗部性狀的綜合評(píng)價(jià)9-11，常用的方法有相關(guān)分析、通徑分析、灰色關(guān)聯(lián)度分析等。崔震海等12研究了在玉米苗期不同供水條件下穗部性狀與產(chǎn)量的相關(guān)關(guān)系，結(jié)果表明在苗期正常供水條件下，百粒重對(duì)產(chǎn)量的影響較大，在苗期輕度干旱脅迫條件下，粒長(zhǎng)和出籽率與產(chǎn)量關(guān)系最密切，在苗期中度干旱脅迫條件下，行粒數(shù)對(duì)產(chǎn)量的影響較大。王久光等13利用通徑分析方法研究了陰雨寡照條件下玉米穗部性狀與產(chǎn)量的相關(guān)規(guī)律，表明穗行數(shù)、千粒重、行粒數(shù)對(duì)產(chǎn)量的直接作用較大，禿尖長(zhǎng)雖然直接作用比較小，但它要通過(guò)行粒數(shù)、千粒重來(lái)對(duì)產(chǎn)量造成影響，在實(shí)際育種中要注意協(xié)調(diào)選擇。岳堯海等14應(yīng)用灰色理論關(guān)聯(lián)度分析法，對(duì)夏玉米7個(gè)數(shù)量性狀與子粒產(chǎn)量進(jìn)行分析，表明各性狀與產(chǎn)量的灰關(guān)聯(lián)度依次順序?yàn)?，葉數(shù)行粒數(shù)株高穗長(zhǎng)百粒重穗行數(shù)，在玉米常規(guī)育種中要兼顧各性狀對(duì)產(chǎn)量的促進(jìn)作用，協(xié)調(diào)一致，充分發(fā)揮雜交種自身增產(chǎn)的最大潛力。主成分分析方法作為一種有效的數(shù)據(jù)統(tǒng)計(jì)分析手段，在玉米穗部性狀研究中應(yīng)用較少，且未形成成熟的系統(tǒng)理論，因此，筆者選取品比圃試驗(yàn)的2個(gè)玉米對(duì)照品種和7個(gè)玉米雜交組合為試驗(yàn)材料，通過(guò)調(diào)查記錄xxxx年3年間11個(gè)玉米穗部性狀因子的數(shù)據(jù)，利用主成分分析法進(jìn)行主成分提取和構(gòu)建函數(shù)評(píng)價(jià)，比較、分析雜交組合與對(duì)照品種間的差異，從而為玉米新品種的選育提供統(tǒng)計(jì)依據(jù)和理論借鑒。 1材料與方法 1.1試驗(yàn)材料 試驗(yàn)于xxxx年在山西省農(nóng)科院小麥研究所玉米試驗(yàn)田品比圃進(jìn)行，試驗(yàn)地前茬為休閑用地，土壤為石灰性褐土，播種前養(yǎng)分含量為：有機(jī)質(zhì)14.30g/kg，堿解氮37.80mg/kg，速效磷11.52mg/kg，速效鉀125.00mg/kg，pH7.48。 材料選取7個(gè)玉米雜交組合和2個(gè)對(duì)照品種，具體種類(lèi)如表1所示，分別于xx年5月22日、xx年5月25日和xx年5月24日播種，播種密度為49500株/hm2，寬窄行種植，寬行行距0.70m，窄行行距0.50m，株距0.33m。采用隨機(jī)區(qū)組設(shè)計(jì)，3次重復(fù)，共27個(gè)小區(qū)，每小區(qū)4行，行長(zhǎng)10m，面積為24m2左右。玉米生育期內(nèi)追肥、除草等田間管理統(tǒng)一進(jìn)行。 1.2調(diào)查項(xiàng)目及方法 在玉米成熟期時(shí)，每小區(qū)隨機(jī)選取10株，調(diào)查株高(m)和穗位高(m)，計(jì)算穗位系數(shù)（穗位高/株高）。待玉米收獲后，從小區(qū)取樣12穗中調(diào)查穗長(zhǎng)(cm)、穗粗(cm)、禿尖長(zhǎng)(cm)、穗行數(shù)、行粒數(shù)、百粒重(g)、單穗重(kg)、單穗粒重(kg)、出籽率(%)各指標(biāo)，調(diào)查參照國(guó)家玉米區(qū)域試驗(yàn)記載標(biāo)準(zhǔn)進(jìn)行。選擇上述11個(gè)穗部性狀（株高除外）作統(tǒng)計(jì)分析。 1.3數(shù)據(jù)分析 數(shù)據(jù)采用Excelxx和SPSS16.0進(jìn)行統(tǒng)計(jì)分析。玉米穗部性狀的不同參數(shù)有不同的量綱，因此有較大的數(shù)量級(jí)差別，為了排除由此對(duì)結(jié)果分析造成的干擾，需要對(duì)原始數(shù)據(jù)進(jìn)行標(biāo)準(zhǔn)化處理。采用標(biāo)準(zhǔn)差標(biāo)準(zhǔn)化，用標(biāo)準(zhǔn)化后的數(shù)據(jù)進(jìn)行主成分分析。 其中主成分評(píng)價(jià)函數(shù)Fp和綜合主成分F的構(gòu)建15采用以下公式： 2結(jié)果及分析 2.1不同穗部性狀的主成分方差分解和特征值 在SPSS軟件中進(jìn)行Factor分析，可得到提取主成分的原始特征值（表2）、主成分載荷矩陣（表3）。從表2中可以看到，在xx年，第1主成分(4.097)的方差貢獻(xiàn)率為37.249%，第2主成分(3.153)的方差貢獻(xiàn)率為28.662%，第3主成分(1.931)的方差貢獻(xiàn)率為17.558%，第4主成分(1.009)的方差貢獻(xiàn)率為9.17%，提取的前4個(gè)主成分累積貢獻(xiàn)率為92.639%，能較好地反映原先11個(gè)穗部性狀數(shù)據(jù)指標(biāo)的大部分信息及變化趨勢(shì)，符合累積貢獻(xiàn)率85%的條件，因此共提取4個(gè)主成分，完成后續(xù)數(shù)據(jù)分析。同理，對(duì)xx、xx年的11個(gè)穗部指標(biāo)分別提取4個(gè)、3個(gè)主成分完成接下來(lái)的分析。 2.2不同穗部性狀主成分載荷值 表3中，xx年在第1主成分上有較大載荷值的性狀分量分別為穗長(zhǎng)(0.933)、單穗粒重(0.885)、單穗重(0.818)、穗位高(0.738)；在第2主成分上有較大載荷值的性狀分量分別為出籽率(-0.784)、禿尖長(zhǎng)(0.707)、穗粗(0.706)；在第3主成分上有較大載荷值的性狀分量分別為穗行數(shù)(0.777)、穗位系數(shù)(0.647)；在第4主成分上有較大載荷值的性狀分量分別為行粒數(shù)(0.477)、百粒重(-0.452)；4個(gè)主成分中每個(gè)主成分均能綜合反映11個(gè)穗部性狀因子的數(shù)據(jù)信息，具體到每個(gè)主成分其優(yōu)先反映的穗部性狀指標(biāo)各不相同，而這4個(gè)主成分能代表原始11個(gè)玉米穗部性狀因子數(shù)據(jù)的綜合信息，從而簡(jiǎn)化了數(shù)據(jù)分析的工作量和避免了變量重疊的誤差。 在xx年第1主成分載荷值高的性狀分量為穗位系數(shù)(0.866)、百粒重(0.854)、單穗粒重(0.79)、穗長(zhǎng)(-0.747)、穗行數(shù)(0.726)；在第2主成分載荷值高的性狀分量為穗位高(-0.615)、出籽率(-0.585)、單穗重(0.584)、禿尖長(zhǎng)(0.578)；在第3主成分載荷值高的性狀分量為穗粗(0.817)；在第4主成分載荷值高的性狀分量為行粒數(shù)(0.652)。在xx年，穗行數(shù)(0.835)、穗位高(0.854)、單穗重(0.78)、單穗粒重(0.716)4個(gè)指標(biāo)在第1主成分有較高載荷；在第2主成分有較高載荷量的指標(biāo)分別是：禿尖長(zhǎng)(0.865)、百粒重(0.833)；第3主成分更多反映地是穗位系數(shù)(0.752)、穗粗(0.732)2個(gè)穗部性狀的綜合信息。在xxxx不同年際間，提取的主成分側(cè)重反應(yīng)的穗部性狀分量各有不同，說(shuō)明主成分分析易受外界環(huán)境影響較大，且隨著調(diào)查數(shù)據(jù)的不同而存在一定差異。 2.3主成分評(píng)價(jià)函數(shù)構(gòu)建及各主成分綜合得分排序 按前述公式，利用穗部性狀在各主成分載荷值、主成分特征值、標(biāo)準(zhǔn)化數(shù)據(jù)，即可構(gòu)建主成分（F1、F2、F3、F4）及綜合主成分F的函數(shù)表達(dá)式（表4）。將各穗部性狀數(shù)據(jù)代入函數(shù)，即可求得各主成分和綜合主成分的得分值，用所得數(shù)值大小對(duì)玉米品種及雜交組合作綜合評(píng)價(jià)（表5），幫助篩選與對(duì)照品種相比評(píng)價(jià)較高的玉米雜交組合。 由表5可知，xx年各玉米品種及雜交組合在第1主成分F1上得分較大的有PB1、PB2、PB7、鄭單958、先玉335；在第2主成分F2上得分較大的有PB1、PB4、PB6、鄭單958、先玉335；在第3主成分F3上得分較大的有PB1、PB3、PB6、PB7；在第4主成分F4上得分較大的有PB1、PB2、PB6、先玉335；在綜合主成分F上得分較大的有PB4、PB6、PB7、鄭單958、先玉335；結(jié)合各主成分及綜合主成分分析的結(jié)果，可以認(rèn)為，與對(duì)照（鄭單958和先玉335）相比，雜交組合PB4、PB6、PB7在綜合穗部性狀上存在一定優(yōu)勢(shì)。同理類(lèi)推，在xx年和xx年中，在主成分得分上與對(duì)照相比存在優(yōu)勢(shì)的雜交組合分別為PB2、PB6、PB7和PB6、PB7。通過(guò)分析不同品種及雜交組合3年穗部性狀主成分得分評(píng)價(jià)結(jié)果，雜交組合PB6和PB7在xxxx年間均表現(xiàn)出優(yōu)于對(duì)照的綜合穗部信息，而穗部性狀與產(chǎn)量密切相關(guān)，因此，雜交組合PB6（昌7-2改良系478改良系）和雜交組合PB7（先玉335變異株選系18599改良系）相比對(duì)照（鄭單958和先玉335）有一定的增產(chǎn)潛力，可作為新品種的候選組合，申請(qǐng)參加省玉米品種預(yù)備試驗(yàn)。 3結(jié)論與討論 玉米雜交組合的篩選以產(chǎn)量為主要依據(jù)，而產(chǎn)量與玉米穗部性狀密切相關(guān)。在常規(guī)育種實(shí)踐中，若穗部性狀的表觀特征只存在微量的差異，則不容易決定組合的取舍。主成分分析法(PCA)可將多變量作降維處理，且原始變量的信息能最大程度地保存，從而為玉米雜交組合穗部性狀的篩選提供了一種較為簡(jiǎn)便、準(zhǔn)確度高的數(shù)據(jù)分析方法。宋江峰16采用主成分分析法對(duì)不同甜糯玉米品種軟罐頭的主要揮發(fā)性物質(zhì)進(jìn)行表征和區(qū)分，確定了不同品種甜糯玉米軟罐頭各自特有的揮發(fā)性風(fēng)味物質(zhì)組合。傅同良17通過(guò)對(duì)33個(gè)糯玉米自交系進(jìn)行主成分分析和遺傳距離分析，初步提出了評(píng)選糯玉米自交系優(yōu)劣的遺傳主成分標(biāo)準(zhǔn)。Kroonenberg8研究了三維聚類(lèi)法和主成分分析法在玉米數(shù)量性狀分析中的應(yīng)用，結(jié)果表明，這2種方法能把變量范圍降維為幾個(gè)易于區(qū)分和計(jì)算的變量組，從而有助于對(duì)玉米新的變量特征的研究。 在本研究中，通過(guò)對(duì)xxxx年間11個(gè)玉米穗部性狀數(shù)據(jù)的主成分分析，主成分提取個(gè)數(shù)、各主成分在穗部性狀上的載荷值在3年間存在一定差異，說(shuō)明主成分?jǐn)?shù)據(jù)分析隨著數(shù)據(jù)采集地點(diǎn)、年份的不同而不同。3年間在不同玉米品種及雜交組合穗部性狀的各主成分及綜合主成分函數(shù)評(píng)價(jià)上，主成分得分表現(xiàn)出相似的變化趨勢(shì)，雜交組合PB6和PB7與2個(gè)對(duì)照品種相比，在主成分得分上表現(xiàn)優(yōu)良，可作為選育新品種的候選組合。數(shù)據(jù)經(jīng)過(guò)PCA方法處理后，11個(gè)穗部性狀指標(biāo)最終可以用34個(gè)主成分來(lái)綜合表征，簡(jiǎn)化了原始多元數(shù)據(jù)的統(tǒng)計(jì)分析，具有一定的實(shí)踐應(yīng)用價(jià)值。因此，作為一種玉米常規(guī)育種的輔助數(shù)據(jù)分析手段，主成分分析(PCA)方法能有效指導(dǎo)育種實(shí)際，縮短育種進(jìn)程，為玉米數(shù)量化育種研究提供了新的方向。 參考文獻(xiàn) 1譚巍巍，王陽(yáng)，李永祥，等.不同環(huán)境下多個(gè)玉米穗部性狀的QTL分析J.中國(guó)農(nóng)業(yè)科學(xué)，xx，44（2）：233-244. 2李召鋒，梁曉玲，阿布來(lái)提，等.不同年代主要玉米品種穗部性狀演變研究J.玉米科學(xué)，xx，18（4）：27-31. 3Lopez-ReynosoDJ，HallauerAR.Twenty-sevencyclesofdivergentmassselectionforearlengthinmaizeJ.CropScience，1998（38）：1099-1107. 4岳競(jìng)之，張莉，張春榮，等.夏玉米主要穗部性狀結(jié)構(gòu)分析和高產(chǎn)育種探索J.玉米科學(xué)，xx，13（3）：47-49，52. 5樂(lè)美旺，張冬仙，熊永鵬，等.黑芝麻數(shù)量性狀的主成份分析及其應(yīng)用J.中國(guó)油料，1995，17（3）：72-74. 6LiuXQ，Chen，XG，WuWF，etal.ProcesscontrolbasedonprincipalponentanalysisformaizedryingJ.FoodControl，xx，17（11）：894-899. 7鮑艷，胡振琪，柏玉，等.主成分聚類(lèi)分析在土地利用生態(tài)安全評(píng)價(jià)中的應(yīng)用J.農(nóng)業(yè)工程學(xué)報(bào)，xx，22（8）：87-90. 8KroonenbergRM，BasfordKE，EbskampAGM.Three-wayclust