高寒區150份豌豆種質資源農藝性狀的遺傳多樣性分析

摘要：高寒區匱乏的豆科牧草限制了其飼草的產量和品質，發展產量高、品質優的豌豆是高寒人工飼草基地建設的關鍵。為篩選高寒區優良的豌豆種質資源，本研究以國內外150份豌豆種質資源為研究對象，采用主成分分析、聚類分析方法對供試豌豆種質資源開展農藝性狀評價及遺傳多樣性分析。結果表明：供試豌豆種質資源的數量性狀遺傳變異系數大，6個質量性狀的平均變異系數與遺傳多樣性指數分別為32.22%和0.96，17個數量性狀的平均變異系數與遺傳多樣性指數分別是41.39%，1.79。17個農藝性狀簡化為6個主成分，累計貢獻率為72.25%；系統聚類分析方法將150份豌豆材料分為產草量、籽粒數量和籽粒大小3大類群。因此，本研究可為適宜高寒區種植優異豌豆種質資源篩選、選育及豌豆種質資源保護利用奠定理論基礎和提供基本依據。

關鍵詞：豌豆；種質資源；農藝性狀；遺傳多樣性

中圖分類號：S326 文獻標識碼：A 文章編號：1007-0435（2023）07-2116-12

Genetic Diversity Analysis of Agronomic Traits of 150 Pisum sativum Germplasm Resources in Alpine Region

ZHANG Peng， BAO Gen-sheng*， JIA Zhi-feng， LIU Wen-hui

（Qinghai university Qinghai Academy of Animal Science and Veterinary Medicine， Key Laboratory of Superior Forage Germplasm in the Qinghai-Tibetan Plateau， Xining， Qinghai Province 810016， China）

Abstract：Shortage of legume with high productivity and quality was considered as one of the vital limited factors in husbandry of alpine regions. Breeding legume species with high productivity and good palatability was critical factor to determine the achievement in establishment of artificial grassland. However，little studies have paid an attention to the breeding of strong adaptability and genetic diversity of Pisum sativum germplasm in alpine regions. Therefore，150 P. sativum germplasms originated from domestic and exotic were evaluated on their agronomic traits and genetic diversities through the principal component （PCA） and cluster analysis. Our results revealed that the coefficient of genetic variation of quantitative traits among the 150 examined P. sativum germplasms were relatively large，the average variation coefficient and genetic diversity index of 6 qualitative traits are 32.22% and 0.96，respectively. Similarly，the average variation coefficients and genetic diversity indexes of 17 quantitative traits were 41.39% and 1.79. Agronomic traits of 150 tested germplasms could be classify into 6 principal components，and the cumulative contribution of those components reached to 72.25% according to the PCA. Furthermore，150 tested P. sativum germplasms could be divided into three group as forage yield group，grain quantity group and seed size group based on the methodical cluster analysis. Our results would provide a theoretical foundation and basic support for selecting and breeding of strong adaptive P. sativum germplasm for establishing high quality artificial grassland in alpine regions.

Key words：Pisum sativum;Germplasm resource;Agronomic traits;Genetic diversity

豌豆（Pisum sativum L.）原產于地中海及中亞地區，是豆科豌豆屬的一年生攀援草本植物，為集糧、菜、肥、飼為一體的多用途作物^[1]，在我國廣泛種植。豌豆嫩莖葉富含蛋白質和多種維生素，已成為人類日糧的重要組成部分^[2]；同時，豌豆根部形成的根瘤能固定大氣中的氮素，將游離態的氮轉化為植物生長所需的氮素，進而有效提高貧瘠土壤的肥力；另外，豌豆植株富含粗蛋白和家畜生長所需的微量元素，成為家畜喜食的優質牧草^[3]。豌豆具有喜冷涼、耐寒旱、耐瘠薄等特點，適于在我國干旱、半干旱地區種植^[4]，其作為優良牧草和輪作倒茬作物在緩解草畜矛盾、提升土壤肥力和生態環境治理方面發揮重要作用^[5-6]。

聯合國糧農組織統計數據表明，全球有99個國家種植豌豆，種植面積排名前五的國家分別為加拿大、俄羅斯、中國、印度和美國^[7]。中國豌豆的種植面積、產量僅次于蠶豆，是我國第二大食用豆類作物^[8]。世界各國均特別重視對豌豆種質資源的收集、保存及利用，全球保存豌豆種質資源數量已高達34 000余份，而我國豌豆種質資源數量少，共征集到蠶、豌豆資源近萬份，尚有部分省及重點地區有待進一步征集^[9]，長期存在品種單一、優良親本匱乏、品種退化等問題^[10]。近年來，國家已將種質資源遺傳多樣性保護研究提升至國家戰略規劃，國內學者先后對豌豆資源遺傳多樣性開展了一系列研究。萬述偉等^[11]對國內外271份豌豆種質資源遺傳多樣性研究發現，初始莢位是決定遺傳多樣性的最關鍵指標，而小區產量是影響性狀變異系數的核心因素。賀晨邦等^[12]對624份國內外栽培豌豆核心種質進行了遺傳多樣性分析，發現其遺傳多樣性豐富，其中20個形態性狀的平均多樣性指數均大于1.369。李群等^[13]基于SSR標記方法對全球228份豌豆種植資源遺傳多樣性研究發現亞洲地區豌豆資源遺傳變異系數最高，而非洲地區遺傳變異系數最低。宗緒曉等^[14]利用SSR標記方法對我國1221份地方豌豆種質資源的遺傳多樣性分析發現內蒙古地區豌豆品種遺傳多樣性最高，遼寧地區最低，闡明原生境的地理、氣候因素對中國地方豌豆品種的遺傳距離存在顯著影響。其中，青海省豌豆種群存在明顯的地域獨特性^[15]，成為全國春播豌豆主要種植區。

青海地區處于青藏高原腹地，寒冷、干旱等極端逆境較多，引進紫花苜蓿（Medicago sativa）、白花草木樨（Melilotus albus）等多年生優良豆科牧草在高海拔地區不能安全越冬，致使豆科牧草匱缺，成為限制青海牧草產量和品質的關鍵限制因子^[16]。一年生禾豆混播草地是解決高寒地區牧草產量低、品質差的有效措施^[17-19]。然而，禾豆混播體系中豆科品種以箭筈豌豆（Vicia sativa）或飼用蠶豆（Vicia faba）為主，面臨可選品種少、高產穩產資源極度缺乏的問題^[20]。基于此，本研究以青海省畜牧獸醫科學院種質資源庫保存的150份豌豆材料為研究對象，采用主成分分析和聚類分析對不同豌豆材料的形態學指標進行分析，評估不同豌豆材料間的遺傳多樣性，以期為適合高寒區優異豌豆材料篩選、提高人工草地產量品質及促進青海草地畜牧業健康發展提供重要的依據。

1 材料與方法

1.1 供試材料

供試150份豌豆材料是20世紀80年代至今廣泛收集并保存于青海省畜牧獸醫科學院草原所種質資源庫，供試豌豆材料基本信息詳見表1。其中，國內品種50份，青海16份，其中地方品種7份，其余43份均為育成品種，西藏和寧夏各7份，陜西5份，北京4份，內蒙古3份，甘肅和江蘇各2份，廣東、香港、黑龍江、四川各1份；國外品種61份，皆為育成品種，其中俄羅斯22份，加拿大17份，南斯拉夫10份，匈牙利5份，美國和英國各2份，荷蘭、保加利亞、澳大利亞各1份，其余39份種質資源來源不明。上述豌豆材料每隔5年在青海省海東市民和縣中川鄉朱家嶺村（35°57′05.0″N，102°47′10.3″E，海拔1 820 m）進行擴繁以保持種子活力。

1.2 試驗地概況

試驗地位于青海省海東市樂都區壽樂鎮熊家灣村（36°33′24″N，102°24′01″E，海拔2 150 m），屬高原大陸性氣候，年均氣溫7.3℃，年均降水量330 mm，降雨集中在每年6—8月，年均無霜期120 d。種植前茬作物為油菜。土壤有機質和全氮含量分別為14.8，1.26 g·kg^-1，速效磷和速效鉀含量分別為13.2和176 mg·kg^-1。

1.3 試驗方法

2020年5月中旬進行播種，播種前對試驗田進行翻地、平整、鎮壓等處理。采用完全隨機區組設計，每個小區種植面積為2 m×5 m，且每份豌豆材料均種植3個小區，小區間設置寬1 m保護行。按照當地種植習慣豌豆種植最佳時期為4月中旬，為確保豌豆材料能正常完成生育期，采用起壟覆膜措施，每壟長4 m、寬30 cm、高25 cm，壟間距20 cm，每個小區種植4壟。采用穴播方法在壟上進行豌豆資源種植，穴間距10 cm，每穴內播種兩粒種子。播種后，壟上鋪覆寬度80 cm的白色農用薄膜塑料，四周用土壓實。播種5天后，每天觀察種子出苗情況，若豌豆子葉破土，立即撕破子葉上面薄膜塑料以防燙傷子葉。由于試驗地灌溉條件便利，分別在開花期和節莢期漫灌1次，出苗期至收獲期進行人工田間除雜3次。

2020年8月，每份供試豌豆材料在每個小區隨機挑選長勢一致5個單株開展數量性狀測定，測定項目包括：莢長、莢寬、主莖節數、節間長度、單株分枝數、始莢節位、單株莢數、每果節莢數、果柄長度、株高、莖粗、鮮莢籽粒數、單株籽粒重、株高、粒長、粒寬和百粒重。9月豌豆成熟期在每個小區隨機挑選5個單株分別對單株質量性狀進行測定，測定項目包括：鮮莢莢形、粒形、種子表面、粒色、臍色、鮮莢色進行測定。其中，農藝性狀觀測依據《豌豆種質資源描述規范和數據標準》進行^[21]。

1.4 數據處理

采用SPSS 22. 0軟件計算供試豌豆數量和質量性狀的平均值（X）、標準差（s）及變異系數，并對供試豌豆材料的數量和質量性狀進行遺傳多樣性分析。其中，遺傳傳多樣性指數（H′，Shannon-Wiener diversity index）的計算方法主要基于供試豌豆材料各觀測性狀的X和s值，劃分為10個等級，第1級[X_i＜（X-2s）]、第2級X₁lt;X₂=（X-1.5 s）lt;X₃=（X-1.0 s）、……第10級[X₁₀≥（X+2 s）]，即X+0 s，X±0.5 s，X±1.0 s，X±1.5 s，X±2.0 s，每0.5 s為1級。計算每一分級占資源總數的百分比為每一級的相對頻率P_i，從而得出多樣性指數^[22]。

式中，P_i為某性狀第i級別內供試豌豆份數占總供試豌豆材料份數的百分比。接著采用主成分分析（PCA，Principal component analysis）和聚類分析（CA，Cluster analysis）分別在SPSS 22. 0以及Origin 2021軟件中完成，其中聚類方法采用系統聚類法，種質間遺傳距離為歐氏距離^[10]。

2 結果與分析

2.1 豌豆種質資源質量性狀特征及遺傳多樣性

供試豌豆材料株型、莢形和種子表面形狀和顏色差異較大。其中，株型可分為：高大分枝少、高大分枝多，矮生型、莖葉細分支多和卷須型等（圖1a-e）；同時，鮮莢莢形可分為：直形、聯珠形，馬刀形，劍形和鐮刀形等（圖1 f-j）；另外，種子依據種子表面和顏色可分為：白色表面光滑、黃色表面褶皺、灰色斑紋表面光滑、灰色斑紋表面褶皺等類型種子（圖2）。可供試豌豆材料6個表型質量性狀分布類型和頻率分析發現，直形莢形占供試豌豆材料鮮莢莢形70%，扁球形、球形和柱形種子分別占供試豌豆種子粒形的36.7%，24.0%和39.3%。表面飽滿光滑種子占供試豌豆材料種子表面54%，斑紋色籽粒占供試豌豆籽粒顏色50%，黃色、黑色和灰白色分別占供試豌豆籽粒臍色總數的46.0%，46.7%和7.3%。綠色和黃色鮮莢分別占供試豌豆鮮莢的66.7%和33.3%（圖3）。豌豆6個表型質量性狀遺傳多樣性分析發現，平均遺傳多樣性指數為0.96，其中，粒色遺傳多樣性指數最高（1.18）而鮮莢色最低（0.64）（表2）；遺傳多樣性指數綜合排序為：粒色＞粒形＞種子表面＞鮮莢莢形＞臍色＞鮮莢色（表2）。另外，變異系數和遺傳多樣性指數間未表現出一致性，例如，鮮莢莢形的遺傳多樣性指數（0.96）略低于平均水平，而變異系數（42.51%）卻高于平均水平（表2）。

2.2 豌豆種質資源數量性狀變異及遺傳多樣性

不同數量性狀間差異大，其中，數量性狀的平均變異系數為41.39%，鮮莢籽粒數變異系數最大（74.51%），其次是單株莢數（68.46%），而鮮莢長變異系數最小（18.84%）（表3）。同時，供試150份豌豆種質資源的17個數量性狀的變異系數均大于10%。另外，150份豌豆種質資源數量性狀的平均遺傳多樣性指數為1.79。其中，鮮莢長的遺傳多樣性指數最高（2.18），主莖節數和株高次之（分別為2.06和2.04），而單株分枝數的遺傳多樣性指數最低（0.94），其余數量性狀的遺傳多樣性指數介于1～2。綜合分析數量性狀的遺傳多樣性指數排序為：鮮莢長＞主莖節數＞株高＞節間長度＞草產量＞始莢節位＞百粒重＞鮮莢籽粒數＞籽粒產量＞每果結莢數＞單株莢數＞果柄長度＞鮮莢寬＞粒寬＞粒長＞莖粗＞單株分枝數。

2.3 豌豆主要性狀相關性分析

150份豌豆種質資源的17個農藝性狀相關性分析發現，大多農藝性狀間存在顯著或極顯著的正相關或負相關關系（表4）。其中，草產量與鮮莢長、鮮莢寬呈顯著負相關（Plt;0.05，表4），與主莖節數、節間長度、單株分枝數、單株莢數、每果結莢數、果柄長度以及株高呈極顯著正相關（Plt;0.01，表4）；籽粒產量與鮮莢籽粒數呈顯著正相關（Plt;0.05，表4），與單株分枝數、百粒重、草產量呈極顯著正相關（Plt;0.01，表4）。

2.4 豌豆種質資源17個數量性狀的主成分分析

對150份豌豆材料的17個數量性狀數據進行主成分分析，PCA分析使150份豌豆材料的17個數量性狀之間的影響更為顯著（圖4）。按照特征值大于1的原則，提取了6個主成分，方差貢獻率為72.25％，其中，第一和第二主成分分別解釋了21.5%和15.8%的方差（圖4）。表5列出了各主成分特征向量，第一、第二主成分特征值為分別為3.65，2.69。第一主成分中鮮莢長、鮮莢寬、始莢節位、粒長、粒寬及百粒重的載荷為負值，其余數量性狀均為正值。其中，載荷較高性狀（株高、草產量、單株分枝數、主莖節數和每果結莢數）均與籽粒產量正相關，說明第一主成分主要反映數量性狀對種子產量增產有明顯的意義。第二主成分中果柄長度、鮮莢寬、莖粗和鮮莢長載荷值均與豌豆籽粒性狀呈正相關，說明第二主成分主要反映豌豆鮮莢數量和大小有關（圖4）。

2.5 豌豆種質資源農藝性狀的聚類分析

150份豌豆種質資源的17個數量性狀、6個質量性狀進行系統聚類分析，將供試豌豆種質資源劃分為3個類群（圖5）。由表6知，第Ⅰ、Ⅱ類群的鮮莢莢形均以直形為主，第Ⅲ類群的鮮莢莢形以馬刀形為主；第Ⅰ、Ⅲ類群種子表面以光滑為主，第Ⅱ類群的種子表面以光滑、凹坑為主；第Ⅰ、Ⅱ、Ⅲ類群的粒色均已斑紋為主；第Ⅰ、Ⅲ類群臍色主要以黑、黃為主，第Ⅱ類群的臍色主要以黑色為主；第Ⅰ、Ⅲ類群鮮莢色主要以綠色為主，第Ⅱ類群的鮮莢色以綠、黃為主；第Ⅰ類群粒形以扁球形為主，第Ⅱ、Ⅲ類群粒形以柱形為主（表6）。此外，鮮莢長、粒長、粒寬、百粒重、莖粗、始莢節位、鮮莢籽粒數在第Ⅰ、Ⅱ、Ⅲ個類群兩兩間均沒有顯著性差異；第Ⅰ類群和第Ⅱ、Ⅲ類群的鮮莢寬無顯著性差異，第Ⅱ、Ⅲ類群的鮮莢寬有顯著性差異（Plt;0.05）；第Ⅰ、Ⅲ類群的籽粒產量、節間長度、株高、主莖節數、單株分枝數均無顯著性差異，第Ⅱ和Ⅰ、Ⅲ類群的籽粒產量、節間長度、株高、主莖節數、單株分枝數均有顯著性差異（Plt;0.05）；第Ⅰ類群和第Ⅱ、Ⅲ類群的果柄長度有顯著性差異，第Ⅱ、Ⅲ類群的鮮莢寬無顯著性差異，草產量、單株莢數、每果結莢數在第Ⅰ、Ⅱ、Ⅲ個類群兩兩間均有顯著性差異（Plt;0.05）。

第Ⅰ類群包含96份材料，占供試豌豆資源總數的64%。該類群果柄長度在3個類群中最高為11.15 mm，變異系數為25.40%；其余16個數量性狀均值均在3個類群中處于中等水平；鮮莢長、粒長、粒寬以及籽粒產量的變異系數均最高，其余數量性狀的變異系數均介于第Ⅱ類群和第Ⅲ類群之間；從來源看，此類群中來自國內的豌豆材料有22份，國外豌豆材料55份，未知豌豆材料19份（圖5）。

第Ⅱ類群包含42份材料，占資源總數的28%，其中鮮莢長、鮮莢寬、粒長、粒寬、始莢節位以及鮮莢籽粒數最高，分別為59.69，11.89，7.24，6.06 mm，3.51個，其余10個數量性狀在此類群中均最低（表6）。鮮莢長、粒長、百粒重和粒寬、單株籽粒產量的變異系數介于3個類群中間和低于其他2個類群，為18.70%，18.95%，37.49%和15.59%，36.10%，其余11個數量性狀的變異系數在3個類群中均最高（表6）。此外，從來源看，此類群中來自國內的豌豆材料有16份，來自國外的豌豆材料有6份，未知來源的資源有20份（圖5）。

第Ⅲ類群僅有12份豌豆種質資源，占資源總數的8%。單株籽粒產量、節間長度、株高、莖粗、草產量、單株莢數、每果結莢數、主莖節數、單株分枝數均值在此類群中的值最高，分別為44.44 g，94.77 mm，181.94 cm，5.35 mm，197.24 g，201.64個，11.24個，28.54個，4.19個，其余7個數量性狀均值較低（表6）；粒寬、果柄長度、主莖節數和百粒重、籽粒產量的變異系數介于其余兩個類群之間和最高，為19.97%，29.72%，17.81%，46.98%，44.12%，其余10個數量性狀的變異系數在3個類群中均最低（表6）；從來源看，該類群的種質資源均來自國內，其中青海8份，西藏4份，該類群共包含12份種質材料（圖5）。

3 討論

遺傳多樣性是物種多樣性和生態多樣性的基礎，也是生物多樣性的核心，主要反映同一種群內部和不同種群間的遺傳變異程度^[23-24]。遺傳多樣性高的物種，其所含的基因變化豐富，具有對環境壓力較強的適應能力以及較大的進化潛力^[25]。豌豆種質資源改良的關鍵是種質資源的有效利用，農藝性狀的鑒定和描述是研究種質資源最基本的方法和途徑^[26]。目前，采用形態性狀來檢測遺傳多樣性是最經典和可行方法^[27]，即在作物的生長發育期間對農藝性狀進行測定和評價。由于表型性狀是基因多樣性的外在表現，因此種質資源間種間差異性在某種程度上也能表現出來^[28]。豌豆基因變化豐富，種質資源農藝性狀多樣，其表現出極多的數量性狀和質量性狀。

3.1 豌豆種質資源表型性狀的遺傳多樣性分析

本試驗發現150份豌豆種質資源的6個質量性狀中鮮莢莢形以直形居多，鮮莢色以綠色為主，籽粒顏色以斑紋為主，籽粒臍色以灰白色資源占比較高。同時，種子表面飽滿光滑居多，但出現了皺褶的種子。原因可能是：（1）種子收獲期誤判，導致豌豆收獲期提前；（2）青海暖季較短造成豌豆生育期延遲，收獲期種子未完全發育成熟，以致于出現皺褶；（3）試驗地缺乏灌溉條件、土壤養分貧瘠，在一定程度上會影響豌豆結莢和籽粒形成過程，造成部分種子表面未完全成熟，種子表面出現皺褶^[29]。豌豆資源具有豐富的遺傳變異，這孫敏杰等^[30]的研究結果較一致。鮮莢莢形、粒形、種子表面、粒色、臍色、鮮莢色等6個農藝性狀變異幅度為25.93%～47.60%，平均變異系數為32.22%，種子表面的變異系數最大，為47.60%；莢長、莢寬、主莖節數、節間長度、單株分枝數、初莢節位、單株莢數、每果節莢數、果柄長度、株高、莖粗、鮮莢籽粒數、單株籽粒重、株高、粒長、粒寬和百粒重等17個數量性狀變異幅度為18.84%～74.51%，平均變異系數為41.39%，鮮莢籽粒數的變異系數最大，為74.51%；這與張佳欣等^[31]研究結果相似，值略高于孫敏杰等^[30]的研究。這種差異的原因可能與資源種類數量及觀測農藝性狀數量有關，部分觀測性狀對整個群體的變異系數產生誤差。單株分枝數、草產量以及株高變異系數較高，表明可以通過一定的技術手段獲得牧草高產品種。

遺傳多樣性指數的高低代表各性狀的分布趨勢^[32]，指數越低說明性狀的分布越集中，指數越高說明性狀呈均衡化趨勢分布^[33]。6個質量性狀和17個數量性狀的遺傳多樣性指數平均值分別為0.96，1.79，這與萬述偉等^[11]研究結果相似。本研究供試150份豌豆種質資源遺傳變異范圍較大，具有豐富的遺傳變異，這與李玲^[34]對642份國內豌豆種質資源遺傳多樣性研究發現中國豌豆種質資源具有更高的形態多樣性的結論相一致。另外，賀晨幫等^[35]對豌豆資源20個農藝性狀進行評價，發現國內外624份不同地理來源的豌豆資源群間的遺傳變異大。鮮莢長、主莖節數以及株高3個性狀的遺傳多樣性指數均大于2.00，表明這3個性狀的分散度好，在優異資源選育過程中可利用性更高，為我國豌豆新品種的選育能提供優異的種質基礎。

3.2 豌豆種質資源農藝性狀的相關性分析

農藝性狀間相關性分析也證明，豌豆大多農藝性狀間存在顯著的相關關系；對17個數量性狀的相關性分析表明，主莖節數、節間長度、單株莢數、單株分枝數、每果節莢數、果柄長度、株高與草產量具有極顯著的正相關性，鮮莢長、鮮莢寬與草產量具有顯著的負相關性，說明鮮莢的大小制約牧草產量，因此在篩選牧草產量高的品種時應重點關注主莖節數、節間長度、單株莢數、單株分枝數、每果節莢數、果柄長度、株高等這些性狀，忽略鮮莢大小；單株分枝數、百粒重與鮮莢籽粒數和籽粒產量具有極顯著和顯著的正相關性，表明單株分枝數、百粒重、鮮莢籽粒數可以促進豌豆高產，在篩選籽粒產量的品種應重點關注單株分枝數、百粒重、鮮莢籽粒數這三個性狀；此外，單株分枝數與草產量和籽粒產量均有極顯著的正相關性，在選用糧飼兼用的豌豆品種時，可篩選單株分枝數高的豌豆品種進行培育；另外，李莉等^[36]研究結果表明，豌豆資源呈極顯著正相關的性狀有主莖節數、株高、果節數、單株莢數，豌豆的單株產量與單株分枝數、單株莢數極顯著正相關，這與本研究結果相似。

對17個數量性狀進行主成分分析，其中前6個主成分的累積貢獻率為72.25%，6個主成分包含的性狀信息具有一定相關性，能充分涵蓋供試150份豌豆資源17個性狀的關鍵信息。同時，6個主成分特征向量間具有相關性，第一主成分說明這些特征向量對產草量的貢獻較大，尤其對株高影響最顯著；第三、四、六主成分與籽粒數量性狀密切相關，進而間接影響籽粒產量；第二、五主成分主要反映對籽粒大小的貢獻情況。由此可見，供試150份豌豆資源可通過產草量、籽粒數量和籽粒大小這三類關鍵性狀對現有豌豆種質綜合性狀進行優劣評價^[37]。

3.3 豌豆種質資源農藝性狀的聚類分析

聚類分析能對種質資源的形態性狀的綜合表現進行聚類，能粗略的反映出種質資源間的親緣關系^[38]。本研究對150份豌豆種質資源的17個數量性狀和6個質量性狀進行聚類分析發現，可將23個農藝性狀相似的150份豌豆種質資源分為3大類群，Ⅰ類群包括96份材料，主要以飼草與籽粒產量居中的豌豆種質資源為主，此類群豌豆材料屬于糧飼兼用型豌豆資源。Ⅱ類群包含42份材料，此類群植株矮小、結莢數少但籽粒性狀優，表明此類群可作為優異種用型雜交材料。Ⅲ類群包括12份材料，此類群籽粒產量和草產量最高且遺傳性狀穩定，表明該類群12份材料適合在可作為高寒區糧飼兼用型種植豆科資源。此外，從3個類群來源地來看，Ⅰ類群國內豌豆材料22份，國外豌豆材料55份，未知來源豌豆材料19份；Ⅱ類群國內豌豆材料16份，國外豌豆材料6份，未知來源20份；Ⅲ類群種質資源均來自國內，其中，青海省8份，西藏藏族自治區4份。可見，地理來源相同的豌豆材料被歸入相同類群，這與王亞飛等^[39-40]的研究結果相似，認為基本按照地理來源分布聚類，部分材料有交叉；同時，Ⅲ類群豌豆資源的飼草和籽粒產量高較高，屬糧飼一體的優質豌豆材料，這類材料可作為高寒區人工飼草基地建植的首選高蛋白型牧草。

4 結論

150份豌豆種質資源在高寒地區開展農藝性狀測定，發現豌豆的23個農藝性狀遺傳變異系數和遺傳多樣性指數均較高。同時，農藝性狀間相關性分析也證明，豌豆大多農藝性狀間存在顯著的相關關系，且主成分分析發現將17個農藝性狀可降維為6個主成分，累計貢獻率為72.25%。此外，系統聚類分析方法將150份豌豆種質資源分為3個類群，地理來源相同的豌豆材料被歸入相同類群。因此，本研究可為適宜高寒區種植優異豌豆種質資源篩選、選育及豌豆種質資源保護利用奠定理論基礎和提供基本數據。

參考文獻

[1] 王娜，馬紹英，馬蕾，等. 肉桂酸和棕櫚酸對豌豆種子萌發和幼苗生長的化感效應[J]. 植物生理學報，2021，57（8）：1657-1667

[2] 向娟，潘紹坤，馬躍洲，等. 我國豌豆尖產業現狀、發展與優勢種植區[J]. 四川農業科技，2021（12）：72-75

[3] CAZZOLA F，BERMEJO C J，GUINDON M F，et al. Speed breeding in pea （Pisum sativum L.），an efficient and simple system to accelerate breeding programs[J]. Euphvtica，2020，216（11）：1-11

[4] DAHL W J，FOSTER L M，TYLER R T. Review of the health benefits of peas （Pisum sativum L.）[J]. British Journal of Nutrition，2012，108（S1）：S3-S10

[5] 王軻. 簡述青海省飼草輪作及加工模式[J]. 養殖與飼料，2021，20（9）：82-83

[6] 何亞慧. 豌豆根瘤菌抗逆、促生性能及固氮效果研究[D]. 蘭州：甘肅農業大學，2011：42-56

[7] LEBAS F，COUDERT P，ROUVIER R，et al. The rabbit：husbandry，health，and production[M]. Rome：Food and Agriculture Organization of the United Nations，1997：1-15

[8] 葛長軍，閆良，徐麗榮，等. 食莢豌豆品種比較試驗[J]. 湖南農業科學，2021，7：11-13

[9] 楊俊品，羅菊枝. 中國蠶、豌豆育種進展[J]. 西南農業學報，1996（S1）：142-146

[10]喬玲，陳紅霖，王麗俠，等. 國外綠豆種質資源農藝性狀的遺傳多樣性分析[J]. 植物遺傳資源學報，2015，16（5）：986-993

[11]萬述偉，宋鳳景，郝俊杰，等. 271份豌豆種質資源農藝性狀遺傳多樣性分析[J]. 植物遺傳資源學報，2017，18（1）：10-18

[12]賀晨邦. 國內外栽培豌豆核心種質形態性狀遺傳多樣性研究[D]. 北京：中國農業科學院，2008：17-18

[13]李群，王棟，張文蘭，等. 基于SSR標記的世界豌豆種質遺傳多樣性分析[J]. 植物遺傳資源學報，2021，22（3）：684-691

[14]宗緒曉，關建平，王述民，等. 中國豌豆地方品種SSR標記遺傳多樣性分析[J]. 作物學報，2008，34（8）：1330-1338

[15]馬進福. 發展青海高原特色旱作農業的比較優勢與思路[J]. 青海科技，2009，16（1）：17-18

[16]DONG S K，SHANG Z H，GAO J X，et al. Enhancing sustainability of grassland ecosystem through ecological restoration and grazing management in an era of climate change on Qinghai-Tibetan Plateau[J]. Agriculture，Ecosystem and Environment，2020，287：106684

[17]WANG Z K，ZHANG X M，MA Q H，et al. Seed mixture of oats and common vetch on fertilizer and water-use reduction in a semi-arid alpine region[J]. Soil and Tillage Research，2022，219：105329

[18]LI R，ZHANG Z Z，TANG W，et al. Common vetch cultivars improve yield of oat row intercropping on the Qinghai-Tibetan plateau by optimizing photosynthetic performance[J]. European Journal Agronomy，2020，117：126088

[19]LI R，ZHANG Z Z，TANG W，et al. Effect of row configuration on yield and radiation use of common vetch-oat strip intercropping on the Qinghai-Tibetan plateau[J]. European Journal Agronomy，2021，128：126290

[20]南志標，王彥榮，聶斌，等. 春箭筈豌豆新品種‘蘭箭3號’選育與特性評價[J]. 草業學報，2021，30（4）：111-120

[21]宗緒曉，王志剛，關建平. 豌豆種質資源描述規范和數據標準[M]. 北京：中國農業出版社，2005：11-14

[22]梁國玲，劉文輝，馬祥. 590份皮燕麥種質資源穗部性狀遺傳多樣性分析[J]. 草地學報，2021，29（3）：495-503

[23]謝文輝，趙文武，王雷挺，等. 22份百脈根種質資源表型數量性狀的遺傳多樣性分析[J]. 草地學報，2023，31（1）：173-179

[24]郭麗芬，張躍，徐寧生，等. 紅花種質資源形態性狀遺傳多樣性分析[J]. 熱帶作物學報，2015，36（1）：83-91

[25]涂敏，王云月，盧寶榮，等. 云南省不同地理位置水稻品種遺傳多樣性分析[J]. 熱帶作物學報，2011，32（6）：998-1003

[26]劉金，關建平，徐東旭，等. 小扁豆種質資源形態標記遺傳多樣性分析[J]. 植物遺傳資源學報，2008，9（2）：173-179

[27]陳珊珊，周明芹. 淺析遺傳多樣性的研究方法[J]. 長江大學學報，2010，7（3）：54-57，65

[28]賀熙勇，倪書邦，陳國云，等. 澳洲堅果種質資源形態性狀的遺傳多樣性分析[J]. 中國農學通報，2010，26（3）：206-215

[29]孫孟園，白云，李振宇，等. ⁶⁰Co-γ射線誘變大豆‘桂夏7號’突變體篩選[J]. 四川農業大學學報，2021，36（6）：737-744

[30]孫敏杰，張天靜，沈寶宇，等. 菜用豌豆種質資源形態性狀遺傳多樣性分析[J]. 種子，2017，36（8）：60-62，68

[31]張佳欣，李靜妍，朱少文，等. 引進豌豆種質資源遺傳多樣性分析[J]. 北方園藝，2021（2）：1-9

[32]曾翰國，岳佳銘，張曙恒，等. 76份老芒麥種質資源的農藝性狀綜合評價[J]. 草地學報，2022，30（11）：3046-3055

[33]于海天，楊峰，呂梅媛，等. 伊朗鷹嘴豆種質資源農藝性狀遺傳多樣性分析及綜合評價[J]. 南方農業學報，2021，52（3）：769-778

[34]李玲. 國內豌豆種質資源形態性狀多樣性分析[D]. 北京：中國農業科學院，2009：58-61

[35]賀晨幫，宗緒曉. 豌豆種質資源形態標記遺傳多樣性分析[J]. 植物遺傳資源學報，2011，12（1）：42-48

[36]李莉，萬正煌，焦春海，等. 外引豌豆資源的鑒定及主要數量性狀的主成分分析[J]. 湖北農業科學，2014，53（23）：5643-5648

[37]馬育華. 植物育種的數量遺傳學基礎[M]. 南京：江蘇科技出版社，1996：76-83

[38]夏利娟，蔡鯤鵬，馬莉娟，等. 小麥品質相關分析和聚類分析[J]. 農學學報，2021，11（8）：1-7

[39]王亞飛，李世景，徐萍，等. 黃淮和長江中下游冬麥區小麥品種（系）農藝性狀及其聚類分析[J]. 中國生態農業學報，2020，28（3）：395-404

[40]尹淑英，VLADIMI F C，NICKOLAY I D，等. 37份山羊豆種質資源的遺傳多樣性SSR分析[J]. 中國草地學報，2016，38（6）：29-35

（責任編輯 彭露茜）