主成分分析法對(duì)甜仁用杏品種 “龍王帽”種子苗優(yōu)良類(lèi)型的篩選

羅新凱,羅振興,王 利,王志遠(yuǎn),田振鐸,楊淑梅,冷志巍,馮淑華

(1.白城市林業(yè)科學(xué)研究院,吉林 白城 137000;2.通榆縣林業(yè)局,吉林 通榆 137200;3.白城市種子管理站,吉林 白城 137000)

甜仁用杏是我國(guó)的特有樹(shù)種，甜杏仁不僅口感香甜，而且營(yíng)養(yǎng)價(jià)值極高，據(jù)馮小雨研究杏仁中粗脂肪含量50.33%，總糖21.24%，蛋白質(zhì)19.25%，還含有K、Ca、Na、Mg、Fe、Zn、Mn和Cu等8種礦質(zhì)元素[1]。目前甜仁用杏主栽品種為“龍王帽”、“一窩蜂”和“優(yōu)一”等，抗寒性好，耐旱能力極強(qiáng)，年降雨量達(dá)到300 mm即可滿足其生長(zhǎng)結(jié)實(shí)的需要[2]，因此適合在我國(guó)北方干旱地區(qū)栽培。但目前我國(guó)甜仁用杏栽培面積極小，年產(chǎn)量?jī)H1 100 t(2008年統(tǒng)計(jì))，人均占有量?jī)H0.55 g[3]，每年還需對(duì)臺(tái)灣和香港大量出口，因此擴(kuò)大甜仁用杏的生產(chǎn)面積對(duì)提高人民生活水平，發(fā)展中國(guó)特色農(nóng)業(yè)經(jīng)濟(jì)意義重大。

吉林省白城市位于科爾沁草原東部，具有大面積的風(fēng)沙地和半山區(qū)丘陵，氣候干旱少雨，非常適合甜仁用杏的栽培，當(dāng)?shù)噩F(xiàn)有“龍王帽”等品種生產(chǎn)園70 hm2，即使靠自然降雨，每年也能有較高的產(chǎn)量，因此當(dāng)?shù)乜纱竺娣e發(fā)展甜仁用杏栽培。

目前甜仁用杏育種主要通過(guò)對(duì)生產(chǎn)園等資源調(diào)查來(lái)發(fā)現(xiàn)一些變異，如1998年遼寧省果樹(shù)研究所選育出“超仁”、“國(guó)仁”、“油仁”和“豐仁”[4]，2007年河北省圍場(chǎng)選育出“圍選1號(hào)”[5]等，通過(guò)雜交育種進(jìn)行新品種選育工作做得極少，白城市林業(yè)科學(xué)研究院通過(guò)播種龍王帽種子2016年從中選育出2個(gè)優(yōu)良無(wú)性系[6]。

主成分分析法可對(duì)多個(gè)變量進(jìn)行綜合分析，趙曉梅采用主成分分析法對(duì)杏品種進(jìn)行了經(jīng)濟(jì)性狀分析并進(jìn)行了優(yōu)良品種選擇[7]，郭愛(ài)華對(duì)杏各品種抗寒性進(jìn)行主成分分析[8]，在育種方面國(guó)內(nèi)對(duì)甜玉米育種和蟹爪蘭資源采用了主成分分析[9,10]，此次我們將對(duì)龍王帽種子苗通過(guò)主成分分析進(jìn)行優(yōu)良品系篩選。

1 材料和方法

2008年播種龍王帽種子，獲得60株實(shí)生苗，2013年開(kāi)始結(jié)果，其中苦仁類(lèi)型52株，甜仁類(lèi)型8株。試驗(yàn)地在白城市林業(yè)科學(xué)研究院甜仁用杏雜種圃內(nèi)，土壤為黑鈣土，pH值8.06，含鹽量0.66%，土壤有機(jī)質(zhì)2.81%；試驗(yàn)樹(shù)木為龍王帽種子苗實(shí)生苗甜仁類(lèi)型，株距2 m，行距1 m，試驗(yàn)時(shí)間為2016年，對(duì)甜仁類(lèi)型調(diào)查果實(shí)產(chǎn)量、核重、仁重、樹(shù)高和地徑，使用spss軟件對(duì)試驗(yàn)結(jié)果進(jìn)行主成分分析，確定優(yōu)良品系。

2 結(jié)果與分析

2.1 各品系主要形狀調(diào)查結(jié)果

各品系主要性狀調(diào)查結(jié)果表1。

將原始數(shù)據(jù)進(jìn)行標(biāo)準(zhǔn)化處理得到表2。

2.2 相關(guān)系數(shù)分析

對(duì)各性狀兩兩之間進(jìn)行相關(guān)性分析，分析結(jié)果見(jiàn)表3。

相關(guān)系數(shù)絕對(duì)值0～0.09為沒(méi)有相關(guān)性， 0.1～0.3為弱相關(guān)，0.3～0.5為中等相關(guān)，0.5～1.0為強(qiáng)相關(guān)。

通過(guò)性關(guān)系數(shù)分析看出，地徑與樹(shù)高之間為強(qiáng)相關(guān)性，與單仁重、單核重和果實(shí)產(chǎn)量為弱的負(fù)相關(guān)性；樹(shù)高與單核重和單仁重為弱相關(guān)，與果實(shí)產(chǎn)量為負(fù)的中等相關(guān)，即樹(shù)體越高的品系產(chǎn)量越低；果實(shí)產(chǎn)量與單核重為中等相關(guān)，與單仁重為強(qiáng)相關(guān)；單核重與單仁重之間為強(qiáng)相關(guān)。

表3 各指標(biāo)相關(guān)系數(shù)

2.3 主成分分析

首先進(jìn)行KMO檢驗(yàn)，5個(gè)性狀指標(biāo)KMO值為0.386，KMO值小于0.5時(shí)不適合做主成分分析，通過(guò)上面相關(guān)系數(shù)分析我們知道地徑只與樹(shù)高呈強(qiáng)相關(guān)，與各經(jīng)濟(jì)指標(biāo)為弱的負(fù)相關(guān)，因此去掉該項(xiàng)指標(biāo)對(duì)其余4個(gè)性狀指標(biāo)重新做KMO檢驗(yàn)，結(jié)果KMO值=0.541，可以進(jìn)行主成分分析。

通過(guò)各指標(biāo)相關(guān)矩陣可求出4個(gè)主成分，并得出各主成分的方差貢獻(xiàn)率，各主成分初始特征值及對(duì)應(yīng)的方差貢獻(xiàn)率見(jiàn)表4。根據(jù)表4，第一和第二主成分值大于1，其方差貢獻(xiàn)率累計(jì)達(dá)到93.117%，保存了各指標(biāo)的絕大多數(shù)信息，因此提出這2個(gè)主成分，記作Y1和Y2。

表4 各主成分方差貢獻(xiàn)率

求出各性狀指標(biāo)的成分系數(shù)得分矩陣，見(jiàn)表5。

表5 各指標(biāo)的成分系數(shù)得分矩陣

通過(guò)表5，得出2個(gè)主成分由各性狀變量的線性組合公式為：

(代數(shù)式1)

可以看出第一主成分Y1主要受單仁重和單核重影響，果實(shí)產(chǎn)量影響也較大，第二主成分Y2主要受樹(shù)高的影響和果實(shí)產(chǎn)量的負(fù)影響，單仁重和單核重得分系數(shù)小，對(duì)Y2影響不大。

將表2的數(shù)據(jù)代入代數(shù)式1中，得出各品系第1和第2主成分得分，用 FAC1和FAC2表示，結(jié)果見(jiàn)表6。

表6 各品系主成分得分

第一主成分“23號(hào)”得分最高，其次是“24號(hào)”和“5號(hào)”，這三個(gè)品系仁重、核重等經(jīng)濟(jì)指標(biāo)高于其他品系。

將表6中各品系2個(gè)主成分得分作出散點(diǎn)圖，見(jiàn)圖1，可更直觀表示各品系的性狀。

圖1 各品系主成分得分散點(diǎn)圖

“2號(hào)”、“5號(hào)”和“23號(hào)”處于第一象限，F(xiàn)AC1>0，因此經(jīng)濟(jì)性狀好，但“2號(hào)”FAC1值太小，與“5號(hào)”相差太大，經(jīng)濟(jì)性狀不夠理想；FAC2>0，因此“5號(hào)”和“23號(hào)”樹(shù)體具有高度優(yōu)勢(shì)，其生長(zhǎng)勢(shì)略強(qiáng)?！?4號(hào)”位于第四象限，F(xiàn)AC1>0并且超過(guò)“5號(hào)”，因此經(jīng)濟(jì)性狀非常好，F(xiàn)AC2<0，說(shuō)明其不具備高度優(yōu)勢(shì)，生長(zhǎng)勢(shì)弱。

其余品系處于第二、第三象限和縱坐標(biāo)軸上，F(xiàn)AC1<=0，經(jīng)濟(jì)性狀都很差。

通過(guò)上述分析，可篩選出“5號(hào)”、“23號(hào)”和“24號(hào)”3個(gè)經(jīng)濟(jì)性狀優(yōu)良的品系。

3 結(jié)論

(1)雜種后代的篩選需要考慮包括經(jīng)濟(jì)性狀在內(nèi)的多個(gè)性狀，在初選階段如果對(duì)各個(gè)性狀逐個(gè)比較分析往往顧此失彼，既抓不住重點(diǎn)，又會(huì)遺漏一些問(wèn)題，因此采用主成分分析法降維，可將所有的性狀指標(biāo)轉(zhuǎn)化成幾個(gè)綜合指標(biāo)，在損失較少信息的前提下分析問(wèn)題既全面又簡(jiǎn)化[11]；筆者試驗(yàn)第一和第二主成分方差累計(jì)貢獻(xiàn)率達(dá)到93.117%，可以表達(dá)全部性狀信息。

(2)第一主成分在單仁重方面的得分系數(shù)最大，單核重略次之，在果實(shí)產(chǎn)量方面得分系數(shù)也較大，因此第一主成分代表經(jīng)濟(jì)性狀，得分高的品系即為優(yōu)良品系，得分最高的前3個(gè)品系“23號(hào)”、“24號(hào)”和“5號(hào)”可確定為優(yōu)良品系；第二主成分主要代表樹(shù)高這一生長(zhǎng)勢(shì)性狀，“23號(hào)”和“5號(hào)”第二主成分得分>0，因此具有一定的生長(zhǎng)優(yōu)勢(shì)。

(3)筆者試驗(yàn)在進(jìn)行主成分分析前進(jìn)行了KMO測(cè)度，KNO測(cè)度是判斷變量(指標(biāo))是否適合做因子分析的統(tǒng)計(jì)檢驗(yàn)方法之一，只有KMO值大于0.5時(shí)，才可進(jìn)行因子分析，在樣品數(shù)量較少情況下，變量過(guò)多也很難滿足KMO值>0.5。筆者試驗(yàn)在開(kāi)始KMO值過(guò)低的情況下刪除了地徑這一相關(guān)性不強(qiáng)的變量，使KMO值符合了因子分析的要求，從而保證分析過(guò)程正當(dāng)合理，結(jié)果正確可信。