43份蘆筍種質(zhì)資源遺傳多樣性評價

張冰冰， 葉艷英， 尹玉玲， 金水鳳， 羅穩(wěn)根， 周勁松， 羅紹春， 湯泳萍

(1.江西省農(nóng)業(yè)科學(xué)院蔬菜花卉研究所，江西南昌 330200； 2.江西省吉安市農(nóng)作物良種場，江西吉安 343016)

蘆筍(L.)又名石刁柏，為多年生宿根性蔬菜，富含皂苷、蕓香苷、黃酮、多糖等多種活性成分，尤其是天冬門酰胺、精氨酸硒、抗壞血酸、膽堿等含量較高，被譽(yù)為蔬菜之王。蘆筍原產(chǎn)于歐亞大陸及北非河谷及草原地帶，清朝末年傳入我國?，F(xiàn)今，我國已成為世界第一大蘆筍生產(chǎn)及出口國，年產(chǎn)量約為780萬t，約占世界總貿(mào)易量的50%，農(nóng)業(yè)產(chǎn)值已接近200億元。蘆筍產(chǎn)業(yè)已發(fā)展成為前景廣闊的朝陽產(chǎn)業(yè)。

我國蘆筍種質(zhì)資源相對貧乏，遺傳背景較狹窄，雌雄異株，育種年限長，嚴(yán)重制約了優(yōu)異品種選育。為解決我國蘆筍產(chǎn)業(yè)品種制約問題，亟須系統(tǒng)深入地進(jìn)行蘆筍新品種選育和繁殖制種技術(shù)研究，篩選適應(yīng)我國氣候條件的蘆筍品種。蘆筍種質(zhì)資源的收集、鑒定與評價、利用是蘆筍品種創(chuàng)制與基因挖掘的重要基礎(chǔ)。目前種質(zhì)資源遺傳多樣性的研究主要有形態(tài)學(xué)、細(xì)胞學(xué)、生化標(biāo)記、分子標(biāo)記等方法。形態(tài)學(xué)方法簡單直觀、易于觀察，是種質(zhì)資源鑒定與評價最直接、最快速、最高效的一種方法。為進(jìn)一步挖掘蘆筍的利用價值，選育高品質(zhì)蘆筍品種，本研究收集43份國內(nèi)外蘆筍新品種(系)，通過設(shè)施避雨栽培技術(shù)，在江西省吉安市農(nóng)作物良種場進(jìn)行比較試驗(yàn)。通過對15個主要農(nóng)藝及品質(zhì)性狀的變異分析、相關(guān)分析、通徑分析、灰色關(guān)聯(lián)度分析、主成分分析和聚類分析，研究其遺傳多樣性，以期為蘆筍資源利用和新品種選育提供參考。

1 材料與方法

1.1 試驗(yàn)材料

共收集國內(nèi)外蘆筍品種(系)43個，均為二倍體品種(表1)，其中京紫龍2號為紫蘆筍，其余均為綠蘆筍。

表1 供試試驗(yàn)材料信息

V16 Aitlas美國V38JX1510江西省V17Supreme美國V39JX1511江西省V18California Early美國V40JX1504江西省V19Purple Passion美國V41豐島2號浙江省V20Green Pride美國V42京綠1號北京市V21WB-210F1美國V43京綠龍 3號北京市V22 Atals美國

1.2 試驗(yàn)設(shè)計

吉安市農(nóng)作物良種場連體大棚試驗(yàn)地為江西省典型紅壤黏土。2017年8月中旬穴盤育苗，10月31日移栽定植，每畦栽單行，株距30 cm，行距為130～160 cm，密度1 500株/667 m。每畦定植2個蘆筍品種(系)，每材料61株，隨機(jī)區(qū)組排列，每品種(系)2個重復(fù)。參照《蘆筍設(shè)施避雨栽培技術(shù)規(guī)程》進(jìn)行田間管理。

1.3 項(xiàng)目測定及方法

1.3.1 產(chǎn)量和商品性 2019年3月1日始光頭筍時期第1次采筍，逐日記載產(chǎn)量至4月1日，并轉(zhuǎn)換成小區(qū)產(chǎn)量(，kg)/667 m進(jìn)行分析。參照Zhu等的方法進(jìn)行分級，嫩莖基部直徑大于1.2cm為優(yōu)質(zhì)筍(，kg/667 m)、大于0.6 cm 為商品性筍(，kg/667 m)。

1.3.2 農(nóng)藝性狀和抗病性 2019年初記錄各種質(zhì)見筍日期，并根據(jù)2018年底割母莖時間12月10日，轉(zhuǎn)換成見筍時間(，d)。每重復(fù)隨機(jī)連續(xù)30株插牌標(biāo)記抽筍數(shù)(，支)。插牌單株地上部生長量最大時，測單株株高(，cm)、莖粗(，cm)、生育指數(shù)(，cm)、第二主莖株高(，cm)和第一分枝高度(，cm)。其中，莖粗為距地面 5 cm 的直徑，公式：生育指數(shù)=株高×莖粗。

江西省位于長江中下游南岸，屬亞熱帶季風(fēng)氣候，雨量豐沛，多年平均年降水量1 638 mm，4—6月為降雨集中期，降水量接近全年一半，是中國多雨省區(qū)之一。蘆筍露天種植，極易受莖枯病(蘆筍癌癥)危害，植株提前枯萎，甚至全田毀滅，故采用大棚避雨栽培生產(chǎn)。整體調(diào)查每個種質(zhì)小區(qū)蘆筍褐斑病、菌核病、根腐病、莖枯病的發(fā)生，評價其抗病性()，分別標(biāo)注好、一般、差。

1.3.3 品質(zhì)性狀測定 計算優(yōu)質(zhì)筍率(，%)與商品筍率(，%)，參照《植物新品種特異性、一致性和穩(wěn)定性測試指南 蘆筍》嫩筍筍頭形狀、鱗芽包裹程度、筍身粗細(xì)對筍均勻度()和筍外觀品質(zhì)()分級，分別標(biāo)注好、一般、差。

1.4 統(tǒng)計分析

分級結(jié)果好、一般、差分別賦值為2、1、0進(jìn)行分析。采用Microsoft Excel 2010處理各性狀數(shù)據(jù)，并計算其最小值、最大值、平均值、標(biāo)準(zhǔn)差、變異幅度、變異系數(shù)；進(jìn)行10級分類，從第1級[<(-2)]到第10級[≥(+2)]，中間每級間差05，為平均值，為標(biāo)準(zhǔn)差。每一級的相對頻率用于計算Simpson多樣性指數(shù)進(jìn)行描述和評價，公式：=1-∑()，表示某性狀第個代碼值(即第級)出現(xiàn)個數(shù)，為總個體數(shù)。 采用SPSS 25.0進(jìn)行相關(guān)性分析、通徑分析、主成分分析、系統(tǒng)聚類分析。采用DPS 7.05進(jìn)行灰色關(guān)聯(lián)度分析。其中相關(guān)性分析采用皮爾遜(Pearson)相關(guān)系數(shù)表示。參考Salkind等對相關(guān)性的分級，0～0.2為極弱相關(guān)或不相關(guān)，>0.2～0.4為弱相關(guān)，>0.4～0.6為中等相關(guān)，>0.6～0.8為強(qiáng)相關(guān)，>0.8～1.0為極強(qiáng)相關(guān)。

2 結(jié)果與分析

2.1 蘆筍種質(zhì)資源不同性狀的多樣性

由表2可知，15個農(nóng)藝及品質(zhì)性狀的變異系數(shù)范圍為8.03%～53.90%，筍外觀品質(zhì)最大，為53.90%，其次是優(yōu)質(zhì)筍產(chǎn)量，為40.90%，見筍時間和商品筍率最小，均為8.03%。變異系數(shù)由大到小排序：筍外觀品質(zhì)>優(yōu)質(zhì)筍產(chǎn)量>抗病性>商品筍筍產(chǎn)量>筍粗細(xì)均勻度>總產(chǎn)量>生育指數(shù)>抽筍數(shù)>莖粗>第一分枝高度>第二主莖株高>株高>優(yōu)質(zhì)筍率>見筍時間=商品筍率。Simpson多樣性指數(shù)范圍為0.279 1～0.850 2，優(yōu)質(zhì)筍產(chǎn)量最大，為0.850 2，第二主莖株高次之，為0.845 9，筍粗細(xì)均勻度最小，為0.279 1，按由大到小排序：優(yōu)質(zhì)筍產(chǎn)量>第二主莖株高>優(yōu)質(zhì)筍率>總產(chǎn)量>莖粗>商品筍筍產(chǎn)量>第一分枝高度>生育指數(shù)>見筍時間>株高>抽筍數(shù)>抗病性>筍外觀品質(zhì)>商品筍率>筍粗細(xì)均勻度。綜上，該批蘆筍種質(zhì)資源類型豐富，遺傳多樣性較高，種質(zhì)之間存在的差異大，有利于優(yōu)異種質(zhì)的比較和篩選。

表2 蘆筍種質(zhì)資源農(nóng)藝與品質(zhì)性狀的多樣性分析

2.2 蘆筍種質(zhì)農(nóng)藝及品質(zhì)性狀的相關(guān)性

表3結(jié)果顯示，蘆筍總產(chǎn)量與優(yōu)質(zhì)筍產(chǎn)量、商品筍產(chǎn)量極顯著極強(qiáng)正相關(guān)；與株高極顯著強(qiáng)正相關(guān)；與莖粗、生育指數(shù)、第二主莖株高、第一分枝高度和優(yōu)質(zhì)筍率極顯著中等正相關(guān)，上述性狀與蘆筍總產(chǎn)量密切相關(guān)。而與其他性狀相關(guān)系數(shù)均小于0.3，呈弱相關(guān)或極弱相關(guān)且不顯著，說明蘆筍總產(chǎn)量與這些性狀無關(guān)。見筍時間與株高極顯著強(qiáng)負(fù)相關(guān)，與莖粗、生育指數(shù)、第二主莖株高、第一分枝高度極顯著中等負(fù)相關(guān)，說明出筍越早株高越低，進(jìn)而影響到其他株型性狀。品質(zhì)性狀中，優(yōu)質(zhì)筍率與總產(chǎn)量、商品筍產(chǎn)量極顯著中等正相關(guān)，與優(yōu)質(zhì)筍產(chǎn)量極顯著強(qiáng)相關(guān)，說明優(yōu)質(zhì)筍率是蘆筍產(chǎn)量品質(zhì)的典型代表性狀。筍粗細(xì)均勻度與抽筍數(shù)、抗病性極顯著中等正相關(guān)，筍外觀品質(zhì)與抽筍數(shù)極顯著弱正相關(guān)，與抗病性極顯著強(qiáng)正相關(guān)，表明蘆筍抽筍數(shù)、抗病性與蘆筍品質(zhì)密切相關(guān)，選育抽筍數(shù)多、抗病性強(qiáng)的材料有利于提升蘆筍品質(zhì)。

表3 蘆筍種質(zhì)農(nóng)藝與品質(zhì)性狀的相關(guān)性分析

2.3 蘆筍種質(zhì)總產(chǎn)量與各性狀的通徑分析

以蘆筍總產(chǎn)量為因變量、其他性狀為自變量進(jìn)行通徑分析。線性回歸方差分析表明，=676.832且達(dá)到極顯著水平，說明蘆筍總產(chǎn)量與各性狀間存在極顯著線性關(guān)系，可以進(jìn)行下一步通徑分析(表4)。各性狀對蘆筍總產(chǎn)量的直接作用大小排序?yàn)椋簝?yōu)質(zhì)筍產(chǎn)量(0.566)>商品筍產(chǎn)量(0.533)>生育指數(shù)(|-0.158|)>莖粗(0.125)>優(yōu)質(zhì)筍率(|-0.124|)>商品筍率(|-0.111|)>第二主莖株高(0.036)>抗病性(0.032)>株高(0.024)>第一分枝高度(|-0.023|)>筍外觀品質(zhì)(|-0.019|)>筍粗細(xì)均勻度(|-0.013|)>抽筍數(shù)(0.011)>見筍時間(|-0.008|)。直接通徑系數(shù)越大，性狀對產(chǎn)量的影響也越大，優(yōu)質(zhì)筍產(chǎn)量和商品筍產(chǎn)量對蘆筍總產(chǎn)量直接貢獻(xiàn)最大。而莖粗、第二主莖株高、抗病性、株高、抽筍數(shù)增加，蘆筍總產(chǎn)量也增加。

表4 蘆筍種質(zhì)資源總產(chǎn)量與各性狀的通徑分析

從間接通徑系數(shù)可知，優(yōu)質(zhì)筍產(chǎn)量通過商品筍產(chǎn)量的間接通徑系數(shù)(0.546)最大，而商品筍產(chǎn)量通過優(yōu)質(zhì)筍產(chǎn)量的間接通徑系數(shù)(0.514)次之，說明優(yōu)質(zhì)筍產(chǎn)量通過商品筍產(chǎn)量對蘆筍總產(chǎn)量產(chǎn)生較大的影響；株高、莖粗、生育指數(shù)、第二主莖株高、第一分枝高度、優(yōu)質(zhì)筍率通過商品筍產(chǎn)量的間接通徑系數(shù)分別為0.321、0.267、0.314、0.313、0.299、0.290，說明這6個性狀通過商品筍產(chǎn)量對蘆筍總產(chǎn)量產(chǎn)生較大影響。結(jié)果進(jìn)一步證實(shí)了優(yōu)質(zhì)筍產(chǎn)量與商品筍產(chǎn)量是直接影響蘆筍總產(chǎn)量的重要性狀，株高、莖粗、生育指數(shù)、第二主莖株高、第一分枝高度、優(yōu)質(zhì)筍率也對蘆筍產(chǎn)量產(chǎn)生較大的影響。

2.4 蘆筍種質(zhì)總產(chǎn)量與各性狀的灰色關(guān)聯(lián)度分析

以蘆筍總產(chǎn)量為母因素，其他性狀為子因素進(jìn)行灰色關(guān)聯(lián)度分析(表5)，依據(jù)各性狀關(guān)聯(lián)度值排序?yàn)椋簝?yōu)質(zhì)筍產(chǎn)量>商品筍產(chǎn)量>優(yōu)質(zhì)筍率>株高>生育指數(shù)>第二主莖株高>第一分枝高度莖粗>莖粗>筍外觀品質(zhì)>抽筍數(shù)>抗病性>商品筍率>筍粗細(xì)均勻度>見筍時間。其結(jié)果與“2.2”“2.3”節(jié)結(jié)果基本一致，可相互驗(yàn)證。

表5 蘆筍種質(zhì)資源總產(chǎn)量與各性狀的灰色關(guān)聯(lián)度分析

2.5 蘆筍種質(zhì)性狀的主成分分析

對43份蘆筍種質(zhì)的15個農(nóng)藝及品質(zhì)性狀進(jìn)行KMO和Bartlett’s檢驗(yàn)，KMO系數(shù)為0.703，Bartlett’s球形檢驗(yàn)=0.000<0.05，說明本研究考察性狀適合進(jìn)行下一步主成分分析。根據(jù)特征值大于1的原則提取了4個因子(表6)，前4個特征根在15個特征根中累計貢獻(xiàn)率達(dá)79.121%，基本能反映全部特征。

表6 最大方差法旋轉(zhuǎn)因子載荷矩陣

K150.2290.1800.745-0.366特征值4.5933.5042.4711.300貢獻(xiàn)率(%)30.62323.35916.4728.667累計貢獻(xiàn)率(%)30.62353.98270.45479.121

第1主成分特征值為4.593，貢獻(xiàn)率為30.623%；正向載荷較高的有株高、莖粗、生育指數(shù)、第二主莖株高、第一分枝高度、見筍時間，這些性狀皆為植株農(nóng)藝性狀指標(biāo)，故第1主成分命名為蘆筍植株農(nóng)藝因子。第2主成分特征值為3.504，貢獻(xiàn)率為23.359%；正向載荷較高的性狀有蘆筍總產(chǎn)量、優(yōu)質(zhì)筍產(chǎn)量、商品筍產(chǎn)量、優(yōu)質(zhì)筍率，均為直接產(chǎn)量因子，故第2主成分命名為蘆筍產(chǎn)量因子。第3主成分特征值為2.471，貢獻(xiàn)率為16.472%，正向載荷較高的性狀有抗病性、筍外觀品質(zhì)、抽筍數(shù)、筍粗細(xì)均勻度，多數(shù)與嫩筍性狀有關(guān)，故第3主成分命名為嫩筍性狀因子。第4主成分特征值為1.300，貢獻(xiàn)率為8.667%，正向載荷較高的僅商品筍率1個指標(biāo)，因此第4主成分命名為商品筍率因子。

根據(jù)第2主成分產(chǎn)量因子各性狀載荷量的大小研究其對總產(chǎn)量的影響，排序?yàn)椋簝?yōu)質(zhì)筍產(chǎn)量(0.929)>商品筍產(chǎn)量(0.876)>優(yōu)質(zhì)筍率(0.776)>生育指數(shù)(0.327)>株高(0.305)>莖粗(0.294)>第二主莖株高(0.292)>第一分枝高度(0.268)>筍外觀品質(zhì)(0.180)>見筍時間(|-0.073|)>商品筍率(0.066)>筍粗細(xì)均勻度(|-0.028|)>抗病性(0.017)>抽筍數(shù)(0.004)。結(jié)果顯示優(yōu)質(zhì)筍產(chǎn)量和商品筍產(chǎn)量依舊是影響蘆筍總產(chǎn)量的直接因素，優(yōu)質(zhì)筍率對產(chǎn)量的影響也較大，農(nóng)藝性狀中生育指數(shù)對產(chǎn)量的影響最大，株高、莖粗、第二主莖株高、第一分枝高度對產(chǎn)量的影響次之。

主成分綜合評價見表7，依據(jù)各種質(zhì)主成分綜合得分進(jìn)行排名，可分為4級。包括3份性狀優(yōu)良(≥1)、15份性狀較好(0.1≤<1)、3份性狀一般(0≤<0.1)、22份性狀較差(<0)。對主成分綜合評價篩選的3份性狀優(yōu)良種質(zhì)(V4、V18、V3)與其田間表現(xiàn)比較分析，具有產(chǎn)量高、株型高大、出筍早的優(yōu)質(zhì)特征，表明主成分綜合評價具有一定的代表性。

表7 蘆筍種質(zhì)資源主成分綜合評價

2.6 43份蘆筍種質(zhì)系統(tǒng)聚類分析

組間連接法對43份蘆筍種質(zhì)的15個農(nóng)藝及品質(zhì)性狀聚類分析，當(dāng)閾值為10.0時將供試種質(zhì)劃分為4個群集(圖1)，各類群的性狀平均值與變異系數(shù)見表8。第Ⅰ群集包括20份種質(zhì)，占供試種質(zhì)的46.51%，主要特征是綜合性狀一般材料，性狀特征介于第Ⅱ群集和第Ⅲ群集之間。

表8 蘆筍種質(zhì)資源各群集性狀特征表現(xiàn)

第Ⅱ群集包括8份種質(zhì)，占供試種質(zhì)的18.60%，主要特征是蘆筍總產(chǎn)量較低(58.17 kg/667 m)，株高(114.38 cm)、莖粗(0.37 cm)、生育指數(shù)(41.88 cm)、第二主莖株高(96.50 cm)和第一分枝高度(27.50 cm)，表現(xiàn)株型矮小，見筍時間晚(115.25 d)。

第Ⅲ群集包括14份種質(zhì)，占供試種質(zhì)的32.56%，主要特征是蘆筍總產(chǎn)量較高(111.42 kg/667 m)，株高(157.71 cm)、莖粗(0.53 cm)、生育指數(shù)(84.25 cm)、第二主莖株高(140.71 cm)和第一分枝高度(45.79 cm)，表現(xiàn)株型高大、見筍時間早(106.93 d)。

第Ⅳ群集包括1份種質(zhì)，占供試種質(zhì)的2.33%，主要特征是蘆筍總產(chǎn)量最高(161.41 kg/667 m)、株高(185.00 cm)、莖粗(0.57 cm)、生育指數(shù)(105.45 cm)、第二主莖高度(165.00 cm)和第一分枝高度(55.00 cm)，表現(xiàn)株型最高大壯碩，見筍時間最早(97.00 d)，即最優(yōu)種質(zhì)。

2.7 優(yōu)異種質(zhì)篩選

根據(jù)性狀鑒定結(jié)果，綜合考慮矮稈(株高≤100 cm)，高稈(株高≥170 cm)，出筍早(見筍時間<100 d)，高產(chǎn)(總產(chǎn)量≥120 kg/667 m)等性狀，篩選上述4類優(yōu)異種質(zhì)共9份(表9)，其中矮稈種質(zhì)2份，為JX1506(江西)、京紫龍2號(北京)，隸屬于聚類群集Ⅱ；高稈種質(zhì)4份，分別為Dogeras、Jersey knight、UC115、Imperial，均收集于美國，隸屬于聚類群集Ⅲ；高產(chǎn)種質(zhì)2份，為巨豐(山東)、JX1516(江西)，隸屬于聚類群集Ⅲ；兼具高稈、出筍早、高產(chǎn)最優(yōu)種質(zhì)1份，為California early(美國)，隸屬于群集聚類Ⅳ，有利于進(jìn)一步生產(chǎn)與育種利用。

表9 蘆筍優(yōu)異種質(zhì)性狀特征表現(xiàn)

3 結(jié)論與討論

種質(zhì)資源遺傳多樣性研究對資源引選、親本組配、雜種優(yōu)勢利用及新品種創(chuàng)制具有重要的指導(dǎo)意義。農(nóng)藝性狀與作物生產(chǎn)直接相關(guān)，易識別鑒定，是作物種質(zhì)研究和遺傳分類的重要依據(jù)。但農(nóng)藝性狀復(fù)雜多樣，易受環(huán)境影響，選擇難度大，因此科學(xué)方法的運(yùn)用顯得尤為重要。本研究搜集43份國內(nèi)外蘆筍新品種(系)，對15個主要農(nóng)藝及品質(zhì)性狀進(jìn)行變異分析、Simpson 多樣性指數(shù)、相關(guān)分析、通徑分析、灰色關(guān)聯(lián)度分析、主成分分析、系統(tǒng)聚類分析，研究其遺傳多樣性，結(jié)果可為育種實(shí)踐提供參考。

3.1 蘆筍種質(zhì)遺傳多樣性

遺傳變異系數(shù)是遺傳變異潛力大小的標(biāo)志，表示群體中直接選擇的范圍。對43份蘆筍種質(zhì)的主要農(nóng)藝及品質(zhì)性狀分析發(fā)現(xiàn)，變異系數(shù)最大為53.90%，最小為8.03%，平均26.33%，Simpson多樣性指數(shù)最大為0.850 2，最小為0.279 1，平均0.749 0，說明該批蘆筍種質(zhì)資源性狀的資源類型豐富，農(nóng)藝性狀差異明顯，具有豐富的遺傳多樣性。優(yōu)質(zhì)筍產(chǎn)量、商品筍產(chǎn)量、莖粗、生育指數(shù)、第一分枝高度、 抽筍數(shù)、 抗病性、筍外觀品質(zhì)等8個性狀的變異系數(shù)均在20%以上，且均大于0.6，說明上述性狀指標(biāo)具有較大的選擇潛力，可利用雜交配組進(jìn)行高產(chǎn)、高抗、品優(yōu)的選育。

3.2 蘆筍種質(zhì)產(chǎn)量與各性狀的多重分析

相關(guān)性分析結(jié)果顯示，各性狀至少與其他1個性狀存在極顯著或顯著相關(guān)。優(yōu)質(zhì)筍產(chǎn)量、商品筍產(chǎn)量與蘆筍總產(chǎn)量極顯著強(qiáng)正相關(guān)，是影響其總產(chǎn)量的直接因素?？偖a(chǎn)量與株高極顯著強(qiáng)相關(guān)，與莖粗、生育指數(shù)、第二主莖株高、第一分枝高度極顯著中等相關(guān)，與優(yōu)質(zhì)筍率極顯著中等相關(guān)，說明這些性狀的變化直接影響蘆筍總產(chǎn)量的變化，其結(jié)論在通徑分析、灰色關(guān)聯(lián)度分析和主成分因子載荷量中都得到相互驗(yàn)證，表明上述性狀的改良有利于快速提升蘆筍品質(zhì)。值得注意的是，見筍時間與產(chǎn)量指標(biāo)呈弱負(fù)相關(guān)且不顯著，與株高、莖粗、生育指數(shù)、第二主莖株高、第一分枝高度呈中等或強(qiáng)負(fù)相關(guān)且極顯著，說明見筍時間通過影響蘆筍株形而影響產(chǎn)量。但蘆筍最重要的采收時間為春季光頭筍時期，出筍越早，越有利于提前上市，產(chǎn)品價格越高。據(jù)調(diào)查發(fā)現(xiàn)，1、2月份春節(jié)期間上市的鮮蘆筍比平時價格高20～30元/kg。因此在創(chuàng)制新品種時應(yīng)注意育種目標(biāo)與生產(chǎn)實(shí)際結(jié)合，選育貼合市場的優(yōu)良品種。陳河龍等認(rèn)為，株高、第一分枝高度、莖粗是蘆筍的主要農(nóng)藝性狀，株高與產(chǎn)量密切相關(guān)。王培等認(rèn)為，蘆筍品種間株高和莖數(shù)差異顯著，與產(chǎn)量密切相關(guān)。李霞等認(rèn)為，優(yōu)質(zhì)筍率是評價蘆筍質(zhì)量的重要指標(biāo)。本研究結(jié)果與前人結(jié)論基本一致，但是主成分累計貢獻(xiàn)率為79.12%，說明還有一些性狀未被考慮在內(nèi)，如一級分枝數(shù)、一級分枝間距、筍鱗片數(shù)量、長寬度等。

3.3 蘆筍種質(zhì)主成分綜合評價與系統(tǒng)聚類分析

主成分分析可將較多的性狀進(jìn)行降維分類，以提煉出較少的比較獨(dú)立的主成分，用于作物農(nóng)藝性狀的評價和篩選，既能把握其綜合性狀表現(xiàn)，又能簡化選擇程序，且更具有科學(xué)性。本研究通過主成分分析，將15個性狀降為4個主成分，分別反映了蘆筍的植株農(nóng)藝因子、產(chǎn)量因子、嫩筍性狀因子與商品筍率因子，累計貢獻(xiàn)率79.121%。通過主成分綜合評價，將43份種質(zhì)劃分為4類，其中3份性狀優(yōu)良，15份性狀較好，3份性狀一般，22份性狀較差。育種工作者可根據(jù)上述4個主要特征及材料群集，有針對性地發(fā)掘育種中間材料，進(jìn)而有效利用資源。

系統(tǒng)聚類分析不僅可揭示品種類群間的遺傳差異，還可了解類群內(nèi)品種的遺傳相似性，為作物新品種選育提供科學(xué)依據(jù)。本研究當(dāng)閾值為10時，將43份蘆筍種質(zhì)劃分為4類，分別為綜合性狀一般種質(zhì)，低產(chǎn)、株型矮小、晚出筍種質(zhì)，高產(chǎn)、株型高大、早出筍種質(zhì)與最優(yōu)種質(zhì)。與主成分綜合評價結(jié)果比較分析，發(fā)現(xiàn)3份性狀優(yōu)良種質(zhì)分別隸屬于第Ⅲ和第Ⅳ群集，性狀較好的15份種質(zhì)除V10和V11隸屬于第Ⅰ群集外，其余均隸屬第Ⅲ群集。表明雖算法不同，但其對種質(zhì)的最終評價基本可以相互印證。

3.4 蘆筍優(yōu)異種質(zhì)篩選

選擇優(yōu)良品種是發(fā)展蘆筍產(chǎn)業(yè)、保障經(jīng)濟(jì)效益的關(guān)鍵。通過聚類分析篩選獲得14份高產(chǎn)種質(zhì)和1份最優(yōu)種質(zhì)。Bussell等研究發(fā)現(xiàn)，通過兩季測產(chǎn)評價蘆筍多年產(chǎn)量，其方法的可靠性被Scholten等在品種選育和生產(chǎn)中驗(yàn)證。因此上述15份種質(zhì)可用于下一步品種生產(chǎn)推廣綜合評價。

除高產(chǎn)種質(zhì)用于推廣生產(chǎn)外，其他優(yōu)異種質(zhì)可用于基因的挖掘和新品種的創(chuàng)制。依據(jù)對蘆筍性狀改良實(shí)際，從43份蘆筍種質(zhì)中篩選了高稈、矮稈、出筍早、高產(chǎn)等優(yōu)異種質(zhì)9份，可作為蘆筍育種的親本材料。