基于主成分分析和聚類分析綜合評價不同品種藜麥在黃土高原地區種植的適宜性

楊釗 劉文瑜 黃杰 魏玉明 謝志軍 李琦 楊發榮

摘要：為了深入了解和利用不同藜麥品種在黃土高原地帶種植的適宜性，以6個藜麥品種（系）為材料，在觀察其物候期的基礎上對13項農藝性狀指標、產量性狀指標和品質指標進行比較分析，并采用主成分分析和聚類分析對不同藜麥品種進行適宜性種植綜合評價。結果表明：不同藜麥品種（系）在黃土高原地帶種植的生育期均有所差別，基本從開花期就初顯端倪。不同品種藜麥在黃土高原地帶種植后其性狀表現中除了水分含量指標之外，其他指標的變異系數均大于15%，并且各指標間均存在相關性，其中產量除與倒伏率間呈負相關外，與其他生長指標間均呈正相關，且與分枝數和有效分枝數間均呈顯著性正相關，與有效分枝數間的關聯程度更高。主成分分析綜合得分排名結果顯示，在黃土高原地帶種植藜麥較佳的品種（系）名稱依次分別為隴藜5號、LXM、臺灣紅藜、隴藜3號、隴藜7號、C07。聚類分析后可將6個品種（系）分為3類，分別為以矮稈為特征的隴藜5號、LXM和C07等品種（系），以株型為特征的隴藜7號和隴藜3號等品種，及以產量為代表的臺灣紅藜。各藜麥品種（系）在黃土高原地帶種植后表現出的性狀具有豐富的多樣性，在以收獲產量為目標的情況下可根據藜麥的主穗長度、主穗直徑等株型性狀進行選擇。臺灣紅藜的優異表現可作為黃土高原藜麥規模化種植及地方品種種質創新的優良品種。

關鍵詞：藜麥；黃土高原；適宜性評價；農藝性狀；產量性狀；品質指標；主成分分析；聚類分析

中圖分類號：S512.904? 文獻標志碼：A

文章編號：1002-1302（2023）24-0021-11

藜麥（Chenopodium quinoa Willd）又被稱為南美藜，是莧科藜屬一年生雙子葉草本植物，在原產地南美洲安第斯山區的秘魯、玻利維亞、厄瓜多爾等國家作為特色作物被當地人稱為“糧食之母”，自英國開始引種并小規模種植之后，到被美國、加拿大等國家陸續引進和大面積馴化種植，距今已有7 000年的栽培歷史［1-2］。20世紀80—90年代，西藏率先開展有關藜麥在我國高海拔地區種植過程中生物學特性評價、栽培育種技術及病蟲害研究等工作［3］。當前，國內藜麥的種植范圍已輻射至甘肅［4］、山西［5］、青海［6］、新疆［7］、內蒙古［8］等地，且種植面積在逐年擴大，相關產業獲得多維度、多方位發展，極大地促進了藜麥產業鏈在我國迅速發展［9］。2014年我國藜麥種植總面積約 3 333 hm2［10］，2018年就發展到1.2萬hm2［10-11］。藜麥作為一種兼具營養與生態價值的作物，在各地“調結構，轉方式，保增收”的農業政策落實和助力全面小康建設中發揮了重要作用［12-13］。

由于藜麥對干旱、寒冷、耐鹽堿等非生物脅迫的抗性較強，在世界各地的栽培面積逐年擴大［14］。但是，不同地區藜麥在株高、穗型、分枝、粒質量等生物學特征上差異較大，且易受到海拔、溫度、機械化水平等條件影響，藜麥的栽培措施具有較大的區域性差異。而黃土高原地處黃河中游，屬于半濕潤區向半干旱區過渡地帶，具有重要的生態服務功能，為區內提供糧食生產、水源涵養、土壤保持等生態系統服務的同時，也向其他區域外溢生態系統服務，尤其向黃河下游提供水源涵養、土壤保持等服務［15］。在植被類型發生變化后，黃土高原地區的土壤生產潛力得到大幅提高，而且其生物多樣性和覆蓋率一定程度上受到氣候改善和地形地貌的變動影響已處于向好向上發展［16-17］。獨特的黃土高原地區氣候和環境為種植藜麥提供了得天獨厚的地理條件。因此，本研究通過對比現有材料多個品種間的統計學差異，構建并分析了各品種在黃土高原地帶適應性種植的拓展方向，擬合出在黃土高原地帶高產量栽培的品種、高營養物質含量栽培的品種、重要性狀賦能新品種選育的品種。

1 材料與方法

1.1 試驗區概況

試驗地位于甘肅省慶陽市寧縣南義鄉馬泉村，屬半干旱氣候，典型的雨養農業區。海拔1 180 m，年降水量480～660 mm，主要集中在7—9月，年均溫10.4 ℃，≥10 ℃年積溫3 305.5 ℃·d，無霜期190 d，全生育期≥10 ℃年有效積溫2 986.75 ℃·d。前茬為玉米，土壤為黃綿土，其pH值為8.60，含有機質、全鹽、速效氮、速效磷、速效鉀分別為12.1%、0.083%、116.0 mg/g、36.0 mg/kg、189.0 mg/kg。

1.2 試驗設計

選用適宜低海拔旱地種植的6個藜麥品種（系）為試驗材料（表1）。試驗田于2021年4月15日在一次性撒施磷酸二銨500 kg/hm2后旋耕（深度為 50 cm），并耙耱掉前茬作物的殘根，將石礫清理出田塊，最后平壓整地，做到土質緊實，墑足平整，隨即覆寬1.2 m的白膜后播種。試驗各品種小區隨機排列組合，各處理設3次重復，共18個小區，小區面積為37 m2，區間留寬1 m走道。株距30 cm，行距 30 cm，覆膜穴播，播種量4.5 kg/hm2，穴播播種深度 2～3 cm。6～8葉期時進行間苗、定苗，留苗1株/穴，保苗約82 500 株/hm2，在藜麥整個生長期間進行中耕人工清除雜草，在花蕾期灌水1次，全生育期灌水2次。其他田間管理同當地高產農田管理。

1.3 測定方法

1.3.1 物候期觀測 每小區連測5株具有代表性的植株，記錄其播種期、出苗期、顯序期、開花期、灌漿期、成熟期等各生育期時間。以田間50%以上植株達到對應生育時期的標準日期為準，其中成熟期以田間藜麥葉片80%枯黃并有部分開始脫落為準。

1.3.2 農藝性狀 每份材料選取5株，從出苗期開始，每隔7 d測量1次，直至成熟期為止。株高（cm）：從地面至植株生長最高部位的距離。分枝數：從貼近地面處數取所有的分枝。有效分枝數：從貼近地面處數取每株分枝數中結實的分枝，記為有效分枝數。主穗長度（cm）：成熟期植株主莖上生長的穗的長度，從主穗基部至穗頂的距離。側穗長度（cm）：側穗基部至穗頂的距離。主穗直徑（cm）和側穗直徑（cm）：用游標卡尺分別測量成熟期植株主莖和側莖上生長的穗的直徑。

1.3.3 產量相關性狀 倒伏率：于成熟期統計每份材料每個小區倒伏數，包括根倒伏和莖倒伏，根倒伏是植株自地表處同根系一起傾斜歪倒，根倒率=根倒伏總株數/總株數；莖倒伏是植株從基部以上某個節位折斷，莖倒率=莖稈折斷株數/總株數，重復3次，取平均值。產量：適時收獲，每株單獨收獲脫粒測產。去除邊行收獲并按小區計產，根據種植的株行距計算單位面積藜麥籽粒產量，隨機選取10株考種，脫粒、晾曬后稱量，重復3次，取平均值。千粒質量：籽粒成熟后風干，數 1 000 粒稱質量，重復3次，取平均值。

1.3.4 品質性狀 先將收獲的藜麥籽粒稱取鮮質量，烘干并測定水分含量；將藜麥籽粒粉碎，過60目篩，測定蛋白質含量和灰分含量。采用烘干法測定水分含量，參照國家標準GB/T 20264—2006《糧食、油料水分兩次烘干測定法》進行。蛋白質含量測定參考農業標準NY/T 3—1982《谷類、豆類作物種子粗蛋白測定法（半微量凱氏法）》，測定過程中使用FOSS2300型全自動定氮儀（丹麥福斯特卡托公司），其中氮轉換成蛋白質的系數為6.25。灰分含量測定參考國家標準GB/T 22510—2008《谷物、豆類及副產品灰分含量的測定》，測定過程中使用智能纖維箱式電阻爐（天津市中環實驗電爐有限公司）。所有測定由甘肅省農業科學院農業測試中心實驗室完成。

1.4 數據分析

采用Microsoft Excel 2010進行數據整理，以“平均值±標準差”表示測定結果，利用SPSS 19.0對所測數據統計分析，用Origin 2019b對統計數據作圖。不同品種（系）間同一指標進行單因素方差分析，用Duncans檢驗法進行差異顯著性比較，并對各性狀指標進行相關性分析、主成分分析及聚類分析，其中主成分分析采用降維-因子分析法，聚類方法采用系統聚類-平方歐氏距離-離差平方和法。

2 結果與分析

2.1 黃土高原地區藜麥生長物候期

表2顯示，隴藜5號、LXM、隴藜7號、C07、隴藜3號和臺灣紅藜的生育期分別為120、120、134、135、131、196 d，從開花期到成熟期分別需要59、57、72、72、68、116 d，從灌漿期到成熟期分別需要43、43、56、50、43、80 d，說明各藜麥品種從物候期上即可表現出差異。品種生育期的差別源自品種的特異性，供試藜麥品種中臺灣紅藜的生育期最長，其次為隴藜7號、C07和隴藜3號，再次為隴藜5號和LXM。除了臺灣紅藜是從顯序期開始就表現出差別之外，其他各品種生育期的差異大致是從藜麥的開花期開始就有差別，而C07和隴藜3號與其他品種相比則主要是在灌漿期出現了較大的差異。總體來看，藜麥的出苗期、顯序期和開花期出現的時間較為集中，而成熟期則較為分散。

2.2 藜麥農藝性狀、產量性狀和品質指標的描述性分析

通過比較不同品種藜麥農藝性狀評價指標的最值（表3）發現，臺灣紅藜在株高、分枝數、有效分枝數、主穗直徑、倒伏率、產量和蛋白質含量等性狀上表現最好，LXM在有效分枝數、倒伏率和千粒質量等性狀上表現最好，隴藜5號在主穗長度、側穗長度、側穗直徑等性狀上表現最好，隴藜7號、C07分別在水分含量和灰分含量等性狀上表現最好。同時，比較分析各項指標的變異系數發現，除了水分含量指標之外，不同品種藜麥在黃土高原地帶種植后其性狀表現的變異系數均大于15%，表明各品種藜麥的水分含量受黃土高原種植環境影響不明顯，但其他性狀所受的地域生長環境影響明顯。變異系數從高到低依次為倒伏率、側穗直徑、側穗長度、株高、千粒質量、產量、分枝數、主穗直徑、主穗長度、蛋白質含量、灰分含量、有效分枝數和水分含量，表明在黃土高原地帶下生長對不同品種藜麥的倒伏率、側穗長度、株高、千粒質量、產量等有較高的影響，而對水分含量、有效分枝數、灰分含量等的影響不明顯。

2.3 藜麥農藝性狀、產量性狀和品質指標的顯著性分析

表4的方差齊性檢驗結果顯示，試驗各品種的株高、分枝數、有效分枝數及主穗直徑等多個指標間存在方差齊性，即可用等方差假設進一步在表2中標注各指標值的顯著性；主穗長度、側穗長度、主穗直徑和倒伏率間不存在方差齊性，則需用異方差進行進一步分析，即使用Brown-Forsythe的修正值標注各指標值的顯著性。同時表3的單因素方差分析顯示，各品種中隴藜5號、LXM、隴藜７號、C07、隴藜3號和臺灣紅藜在株高、分枝數、有效分枝數、主穗長度、側穗長度、側穗直徑、倒伏率、產量、千粒質量、蛋白質含量、水分含量及灰分含量等指標上存在極顯著差異，在主穗直徑指標上存在顯著差異。由表3可以看出，株高以臺灣紅藜最大，與其他品種間存在極顯著差異；分枝數也以臺灣紅藜最大，與其他品種間存在極顯著差異；有效分枝數以臺灣紅藜和LXM最大，與其他品種間存在極顯著差異；主穗長度以隴藜5號最大，與LXM、隴藜3號和臺灣紅藜間無顯著差異，與其他品種間存在顯著差異；側穗長度以隴藜5號最大，與其他品種間存在顯著性差異；主穗直徑以臺灣紅藜最大，除與隴藜3號有顯著性差異外，與其他各品種間無顯著性差異；側穗直徑以隴藜5號最大，與其他各品種間均存在極顯著性差異；倒伏率在各品種間無差異；產量以臺灣紅藜最大，與其他各品種間均存在極顯著性差異；千粒質量以LXM最大，與隴藜5號無顯著性差異，與其他各品種間存在極顯著性差異；蛋白質含量以臺灣紅藜最大，與其他各品種間存在顯著性差異；水分含量以隴藜7號最大，與其他各品種間存在顯著性差異；灰分含量以C07最大，與其他各品種間存在顯著性差異。

2.4 藜麥農藝性狀、產量性狀和品質指標的相關性分析

藜麥各項農藝性狀之間的相關性分析結果（表5）顯示，藜麥的株高除與側穗長度呈負相關外，與其他生長指標均呈正相關，其中與分枝數呈極顯著正相關；分枝數與側穗長度間呈負相關，與其他生長指標間均呈正相關；有效分枝數與其他指標間均呈正相關；主穗長度與其他生長指標間均呈正相關，且與主穗直徑間呈顯著性正相關；側穗長度與株高和分枝數間均呈負相關，但與有效分枝數、主穗長度、主穗直徑及側穗直徑間均呈正相關；主穗直徑、側穗直徑與其他指標間均呈正相關。

藜麥產量性狀與各農藝性狀間的相關性分析結果顯示，倒伏率與側穗長度和側穗直徑間呈正相關，與其他指標均間呈負相關；產量除與倒伏率間呈負相關外，與其他生長指標間均呈正相關，且與分枝數和有效分枝數間均呈顯著性正相關，特別是與有效分枝數間的關聯程度更高；千粒質量與株高、分枝數、主穗直徑間呈負相關，與有效分枝數、主穗長度、側穗長度、側穗直徑和倒伏率間呈正相關。

藜麥品質性狀與農藝性狀間的相關性分析結果顯示，藜麥的蛋白質含量除與側穗長度和倒伏率間呈負相關外，與株高、分枝數、有效分枝數、主穗長度、主穗直徑、側穗直徑間呈正相關，且與株高、分枝數間呈顯著性正相關，特別是與分枝數的關聯程度更高；水分含量除與株高、倒伏率間呈正相關外，與其他生長指標間均呈負相關；灰分含量與株高、分枝數、有效分枝數間呈正相關，與主穗長度、側穗長度、主穗直徑、側穗直徑和倒伏率間均呈負相關。

藜麥各項品質指標間的相關性分析結果顯示，藜麥的蛋白質含量與水分含量呈負相關，與灰分含量呈正相關；水分含量與灰分含量間呈正相關；灰分含量與蛋白質含量間的關聯程度高于其與水分含量間的關聯程度。

2.5 藜麥農藝性狀、產量性狀和品質指標的主成分分析

2.5.1 主成分提取及成分荷載矩陣 從表6可以看出，KMO統計量大于0.5，且Barlett球形度檢驗顯著性P<0.5，即檢驗陣是單位陣，說明各檢驗變量間的相關性足夠小，具備做因子分析的條件。圖1主成分碎石圖中曲線坡度的陡峭程度表明，3個主成分的特征值呈急劇下降的趨勢顯示，隨著主成分特征值的增加，從第3個主成分開始，主成分的陡峭程度變化呈緩慢下降的趨勢，表明前3個主成分因子對解釋原變量（藜麥在黃土高原生長適宜性評價的所有信息）的貢獻較大，選取3個主成分作為對原變量的分析方向。同時，由表7可知，前3個主成分的累計方差貢獻率為84.937%，即這3個主成分包含了原始變量84.937%的信息，具有較高的代表性，且其特征值分別為5.083、4.066、1.894，均大于1，表明這3個主成分基本涵蓋了不同品種藜麥在黃土高原生長適宜性評價的所有信息。從公因子方差結果可知，每一個指標的共性方差大部分在0.5以上，且大多數超過0.8，說明這3個公因子能夠較好地反映原始各項指標的大部分信息。結合上述分析結果，本研究選取前3個主成分作為藜麥在黃土高原生長適宜性評價分析的指標。

由表8可以看出， 旋轉成分載荷是主成分與變量的相關系數，其數值大小表示原有變量在降維后所構成的綜合變量中的貢獻率高低，對于一個變量而言，載荷絕對值越大的主成分與其關系越密切，也越能代表這個變量。結合表6和表7可知，主成分1的特征值和累計貢獻率分別為5.083和39.098%，與產量具有相對較高的載荷絕對值，故主成分1代表的是藜麥主要農藝性狀和品質指標中產量的構成因子。主成分2的特征值和累計貢獻率分別是4.066和70.372%，在該主成分中千粒質量的載荷絕對值最大，表明主成分2體現的是藜麥主要農藝性狀和品質指標中千粒質量的相關信息。主成分3 的特征值和累計貢獻率分別為1.894和84.937%，與主穗直徑具有相對較大的載荷絕對值，表明主成分3代表的是藜麥主要農藝性狀和品質指標中主穗直徑的相關信息。

2.5.2 主成分得分及綜合得分評價 由表9可知，根據因子得分系數矩陣及其所對應的主成分，可以計算出主成分因子得分，公式如下所示：

F1=0.222X1+0.206X2+0.026X3+0.037X4-0.119X5+0.087X6-0.007X7-0.019X8+0.124X9-0.210X10+0.205X11+0.051X12+0.075X13；

F2=0.046X1-0.022X2-0.139X3+0.284X4+0.122X5+0.283X6+0.212X7+0.060X8-0.012X9+0.09X10+0.055X11+0.118X12-0.296X13；

F3=-0.073X1+0.06X2+0.326X3-0.061X4+0.171X5-0.149X6-0.002X7+0.154X8+0.202X9+0.080X10+0.006X11-0.363X12+0.119X13。

式中：X1～X13分別表示株高、分枝數、有效分枝數、主穗長度、側穗長度、主穗直徑、側穗直徑、倒伏率、發芽率、產量、千粒質量、蛋白質含量、水分含量和灰分含量；F1～F3分別表示各主成分得分。

結合主成分因子得分公式并以各主成分方差貢獻率為權重，構建藜麥在黃土高原生長適宜性綜合評價得分函數，公式如下所示：

F=0.390 98F1+0.312 74F2+0.145 65F3。

式中：F為不同品種藜麥在黃土高原生長適宜性評價的綜合得分。

根據綜合評價得分函數可計算出不同品種藜麥在黃土高原生長適宜性的綜合得分及其排序，其中F值越大，表明該品種藜麥在黃土高原地區生長的適宜性越好，結果如表9所示。針對6個不同品種藜麥在黃土高原地區種植適宜性的綜合排名的順序而言，臺灣紅藜的評價得分最高，其次為隴藜7號、隴藜3號、隴藜5號、LXM和C07。其中，僅臺灣紅藜的綜合得分＞50，遠高于其他品種，說明這個品種相對于其他品種更適宜在黃土高原地區種植。隴藜7號、隴藜3號、隴藜5號、LXM的綜合得分均在32～36分之間，相差不大，均可視為相對適宜。而C07的綜合評分是28.14，與其他品種相比處于最低水平，說明其在黃土高原地帶種植的適應性相對較差。

2.5.3 基于主成分的二維排序分析 在圖2和圖3中，主成分1軸是原始數據矩陣方差最大方向，主成分2、3軸是方差次大方向，降維時，最大限度保留有用信息，尋找能夠使不同產地樣品得到最大分離的方向。基于不同品種在黃土高原生長適宜性評價的主成分分析結果，以產量相關信息為主要代表的主成分1值為橫坐標，以千粒質量相關信息為主要代表的主成分2值為縱坐標繪制二維排序圖（圖2），產量和千粒質量均是越高越好，符合這2個指標所在象限的只有隴藜5號和LXM，臺灣紅藜是所有品種中距離產量載荷線最近的，說明關聯程度最高，而其他相對距離產量和千粒質量成分值由近及遠的依次是隴藜麥3號、隴藜7號和C07。以產量相關信息為主要代表的主成分1值為橫坐標，以主穗長度相關信息為主要代表的主成分3值為縱坐標繪制二維排序圖（圖3），產量和主穗長度均是越高越好，雖然沒有符合這2個指標所在象限的藜麥品種，但是臺灣紅藜是所有品種中距離產量和主穗長度的載荷線最近的，說明關聯程度最高，而其他相對距離產量和主穗長度成分值由近及遠的依次是隴藜麥3號、隴藜5號、LXM、隴藜7號和C07。

同時，在圖2和圖3中以主成分1為橫軸的表達可以看出，主穗直徑、有效分枝數、主穗長度、蛋白質含量、分枝數、株高有較高的載荷值（>1）；以主成分2為縱軸的表達可以看出，千粒質量、側穗長度有較高的載荷值（>1）；以主成分3為縱軸的表達可以看出，各項指標的載荷值均小于1，差異趨勢不大。

綜上，在黃土高原地帶種植藜麥較佳的品種分別為隴藜5號、LXM、臺灣紅藜、隴藜3號；而隴藜7號和C07品種在黃土高原地區的種植適宜性不佳。

2.6 聚類分析

系統聚類分析法可將受試樣本劃分為不同類群進行評價分析，其中對樣品的聚類稱為Q型聚類，對變量的聚類稱為R型聚類，其結果具有綜合性、客觀性和科學性。本試驗是在主成分分析的基礎上，采用系統聚類分析法對不同品種藜麥在黃土高原地帶種植后的各項農藝性狀和品質指標做R型聚類，同時對6個藜麥受試品種進行Q型聚類。其中，聚類方法采用組間聯接法，聚類區間為平方歐氏距離，聚類結果如圖4和圖5所示。

從圖4中可以看出，當聚類距離為2.5時，可將藜麥在黃土高原種植適宜性的農藝性狀、產量性狀和品質指標的評價指標分為3類，第1類為倒伏率、千粒質量、側穗直徑、灰分含量、主穗直徑、蛋白質含量、有效分枝數、水分含量、側穗長度、主穗長度；第2類為產量；第3類為株高；同時結合主成分載荷分析結果可知，最終選用藜麥的產量、千粒質量、主穗直徑和株高作為綜合評價藜麥在黃土高原地區種植生長適宜性的關鍵性指標。

從圖5中可以看出，當聚類距離為2.5時，可將藜麥在黃土高原適宜性種植的品種分為3類，第1類為隴藜5號、LXM和C07，第2類為隴藜3號和隴藜7號，第3類為臺灣紅藜。根據表10可知，第1類的3個藜麥品種與株高有很大的代表性（特征值為152.866 67），第2類的2個藜麥品種與主穗長度有很大的代表性（特征值為54.03），第3類的1個藜麥品種與產量有很大的代表性（特征值為148.5）。與之前的分析相比，得到大致相同的分析結果，因此采用本試驗所選的主要品質進行聚類分析是可行的。

3 討論與結論

我國引種藜麥的時間較短，在宏觀和微觀上的研究尚處于起步階段，品種的田間農藝性狀及遺傳性狀不穩定［18］。同時，藜麥作為一種耐鹽植物，有些藜麥品種可以耐受400 mmol/L NaCl溶液的鹽分脅迫［19］；藜麥耐堿性較強，適宜在肥力薄弱的土壤中種植，具有非常廣泛的適種性［20］。因地理來源、創質途徑、品種性狀等不同，在同一地區栽培表現出的差異性即可說明不同種質在表型特征方面有豐富的遺傳多樣性［21］。結合本試驗內容從現有的資料來看，在東鄉半干旱地區種植的隴藜5號折合產量為3 278.10～4 046.25 kg/hm2，比同地區參試品種臺灣紅藜和隴藜7號都高（P＞0.05），在產量結果排名當中這3個品種的名次分別為1、3和最末，與各自的倒伏率排名趨同（分別為1、2、0），其中臺灣紅藜處于產量最末且倒伏率為“0”的緣由是生育期較長導致不能成熟、無法結實，經過綜合評價，發現隴藜5號在該地區表現出良好的抗倒伏性、較好的豐產性以及優異的早熟特性，具有更高的推廣價值［22］。與隴藜5號相比，隴藜7號生育期要長 17 d （分別為119、102 d）且倒伏率高1百分點（分別為2%和1%），隴藜7號、隴藜5號的生育期時間前后相差14 d（分別為136、122 d）且產量上前者遠低于后者，這與本研究結果大致吻合。王志恒等對由主成分分析十余種藜麥品種結果顯示，隴藜3號的抗旱性最好，但抗旱性與耐鹽堿性呈負相關，故其耐鹽堿性的排名處于最末［23］。由此可見，隴藜3號較適用于在干旱半干旱地區進行引種，但還得考慮海拔和其他氣候環境條件。如在隴東地區的10個品種（系）比較結果表明，隴藜3號的主成分得分綜合排序處于第10位，聚類分析結果劃為早熟品種類，雖然產量性狀和品質指標不突出，可選擇在高海拔無霜期較短或在夏閑期短地區復播種，亦可作為冬小麥與藜麥輪作中的選擇之一［24］。而在低海拔地區——北京周邊栽植多個品種（系）的藜麥中，有結果表明隴藜3號在≥0 ℃積溫2 500 ℃以上的季節中，適播期廣，尤其是5月中下旬至6上旬均可適宜于≥500 m的山區播種，且以6月上旬播種豐產性、抗倒伏性、景觀效果最好［25］。同時，在高海拔地區——甘肅省張掖市山丹縣（海拔2 300 m左右，年平均氣溫0～2 ℃）種植的隴藜1號至4號及青藜等品種中，由于海拔、有效積溫及早、晚霜發生情況等氣候條件的制約，只有隴藜3號、青藜能夠正常成熟，生育期分別為158、163 d，折合產量分別為 2 973、2 734.5 kg/hm2［26］。但是，在甘南地區（海拔2 737 m，年均氣溫2.8 ℃）種植的隴藜1號至4號卻都能夠成熟，其中隴藜3號雖然單株穗粒數最多（為20 811粒），但是千粒質量最輕（2.65 g）、產量最低（2 167.5 kg/hm2），產量要比同品種隴藜1號低54.5%［27］。本研究結果顯示，隴藜3號的千粒質量和折合產量分別為0.32 g和 1 200 kg/hm2，更是遠低于在甘南地區種植的結果，這表明隴藜3號適合在緯度高、海拔高的地區種植。山丹縣年平均氣溫5.8 ℃，無霜期平均為138 d，日照時數平均為 2 993 h，年降水量197 mm，蒸發量為 2 246 mm［28］，可作為隴藜3號的參考種植氣候環境。LXM藜麥在天祝寒旱山區種植的生育期為168 d，折合產量為3 003.27 kg/hm2，抗倒伏性強，綜合排名為第五，因此并不建議在寒旱山區種植［29］。本研究中，LXM的折合產量為1 575 kg/hm2，低于在天祝寒旱山區種植的產量（低幅約47.6%），表明該品種在黃土高原地帶相對于寒旱山區更不太適宜種植。原產于中國臺灣的特征種臺灣紅藜，從SSR標記種質群體的親緣關系進行分析，再到主成分分析驗證，均證實臺灣紅藜的地理區域特征［30］。臺灣紅藜其秸稈纖維含量不僅較高，籽粒的粗蛋白含量還達到22.03%，因此可作為青貯飼料用［31-32］。臺灣紅藜在北京門頭溝區引種栽培后的平均產量為 2 670 kg/hm2，遠高于在黃土高原地帶種植的產量（2 227.5 kg/hm2）［33］，其差距來源可能包括經度、緯度、海拔及環境溫度等。分析四川涼山州地區的藜麥營養成分，紅藜的氨基酸、粗脂肪、多糖、多酚、黃酮含量均高于隴藜3號，表明紅藜具有更高的食用價值［34］。在本研究中，臺灣紅藜的蛋白質含量確實高于隴藜3號（高約35.8%），但水分含量和灰分含量均低于隴藜3號（分別低約8.8%和10.5%），并且在各參試品種當中臺灣紅藜在株高、分枝數、有效分枝數、主穗長度、主穗直徑、產量及蛋白質含量上均表現出顯著的優異性。

總體而言，參試的各藜麥品種均能夠在黃土高原地區正常成熟，表明黃土高原地區氣候條件適宜飼用藜麥種植，也說明在這個生態環境下能夠獲得較全面且穩定的表型數據。從對在黃土高原地帶所種植藜麥的物候期觀察及描述性分析結果中可知，不同的藜麥品種在黃土高原地區種植后各品種間存在著不同程度的差異性（P＜0.05）。生育期時間從長到短依次為臺灣紅藜＞隴藜7號＞C07＞隴藜3號＞隴藜5號＞LXM。在藜麥各項性狀指標測定當中的最值和平均值上，臺灣紅藜占據多項與產量和品質相關的最佳值指標。同時，通過對比分析供試藜麥各品種的各項指標的變異系數發現，在黃土高原地帶下生長對不同品種藜麥在倒伏率、側穗生長、株高等方面均表現出豐富的變異，說明至少在藜麥的產量方面其遺傳多樣性指數較高，遺傳差異大且資源類型豐富，有利于特異種質材料的比較和篩選，對利用現有種質提高藜麥產量的潛力較大，但也相對反映出迭代或更換品種資源對創新當地藜麥的品質性狀或不太顯著。這與Curti等在阿根廷北部所做的藜麥品種對比分析試驗結果中有關藜麥種質的表型受地域和環境影響較大，表現為高度多樣性［35］。

從相關性分析的結果來看，在黃土高原地帶種植的藜麥株高、蛋白質含量與主穗直徑、側穗直徑間呈正相關，產量與株高呈正相關。這與王思宇等的研究結果相契合［36-37］。說明藜麥株高和直徑的增加會帶動藜麥對干物質的儲運和轉移，因為作物的莖主要起到對物質進行轉運流通的功能［38］，而株高的增加則一般會相應地增加葉片數和葉面積大小，雖然有關葉片的農藝性狀并沒有在本研究中體現，但是從其他作物的資料來看，確實有相關的證實材料［39］。此外在品質表現上，本研究的結果顯示，藜麥的水分含量與蛋白質含量呈負相關，但與灰分含量呈正相關，藜麥的蛋白質含量與灰分含量呈正相關。這與胡一波等的研究結果［40］相同，其研究中進一步指出，蛋白質含量與脂肪、淀粉及總黃酮含量間均呈負相關，而與總多酚及總皂苷的含量間呈正相關［40］。而總多酚、總皂苷和總黃酮含量的高低一般是用來評價藥用植物的藥用有效成分含量的高低［41］，進一步說明藜麥是一種全營養功能作物，其所含16種氨基酸、礦質元素及酚類和黃酮類等多種化學活性成分不僅可作為食材，也可作為高效的優質藥材［42］。

從主成分及綜合得分分析結果來看，通過對黃土高原地帶種植的6個藜麥品種13個性狀指標的適宜性表現和評價的分析中發現，前3個主成分的累計貢獻率為84.937%，其所表達的綜合信息可以用來表達全部性狀的84.937%，從而簡化各品種間比較性狀的數量，提高篩選評價效率。從綜合評價結果來看，3個主成分依次反映了藜麥的產量特征、千粒質量特征和主穗直徑特征。其中，臺灣紅藜的得分最高，綜合評價最優，說明臺灣紅藜非常適宜在黃土高原地帶種植，這可能與其生育期長的種質特性有關。在本研究中的3個隴藜系列品種間及LXM品種的綜合得分相差并不大，說明隴藜系列品種和LXM品種在黃土高原地帶均具有推廣應用價值。與其他藜麥品種相比，C07藜麥在黃土高原地帶種植后生長適宜性評價綜合表現最差，說明其在黃土高原地帶種植的利用價值不大。有研究指出，適期播種時，相應的氣候因子對藜麥的干物質積累動態、光合性狀差異和產量性狀差異等表現影響較大（P＜0.05）［43］，因此根據品種的適應性表現，調節播期，這或許是一種提高藜麥產量和品質的方法。

在對6份藜麥種質材料13個性狀指標進行聚類后發現，可將6份藜麥品種聚為3個群類，從R型聚類的結果并結合主成分載荷分析結果最終選用藜麥的產量、千粒質量、主穗直徑和株高作為綜合評價藜麥在黃土高原地區種植生長適宜性的關鍵性指標。這4個指標從農藝學角度可簡單概括為產量性狀和株型性狀，與張亞萍的研究結論［44-45］相一致。進一步說明藜麥的產量主要積累于主枝穗的大小，故藜麥主枝穗的大小是選育高產藜麥品種的重要指標［46］。同時，R型聚類的結果中4個指標沒有與主成分分析中綜合表現相近的指標聚為一類，這與魏玉明等在隴東地區的試驗結果［24］相類似，證明藜麥這類異型雜交（異花授粉）作物非常容易受外來基因的干擾，異型雜交可高達17.4%［47］，其個體差異大于品種間的差異。

從Q型聚類的結果可知，6份材料可以分為3類種質材料。以株高為代表（特征值為 152.866 67）的第1類藜麥包括隴藜5號、LXM和C07，根據長勢可以看出其為矮稈作物，但千粒質量高（即籽粒飽滿），說明這3個藜麥品種（系）資源可以嘗試以特異種質在藜麥密植矮稈型品種（系）或旅游區景觀植物品種（系）選育目標中作為的親本之一。其次是以主穗長度為代表（特征值為54.03）的第2類藜麥包括隴藜3號和隴藜7號，這類藜麥主穗長度處于適中程度，并與主穗直徑、株高等藜麥生長的長勢狀況均呈正相關關系，也從側面說明這2個藜麥品種在黃土高原地帶具有廣泛的適種性，可以考慮將其作為秸稈飼料作物或藥材作物。最后就是以產量為代表（特征值為148.5）的第3類藜麥包括臺灣紅藜，這類藜麥除了側穗直徑、側穗長度、倒伏率、千粒質量、水分含量和灰分含量等方面表現不是最突出外，其他方面均表現出與其他同試品種（系）的差異性較大，這一點在圖3中可以得出參照。可見，臺灣紅藜不僅非常適宜于在黃土高原地帶推廣種植，而且還可以在以遺傳改良當地更為普適性種植、高附加值加工產出的藜麥品種中發揮優質藜麥親本材料的作用。

由于藜麥的開發應用涉及原糧制造、食品加工、藥用開發與保健品、秸稈飼用、景觀綠化等眾多領域［18］，所以更好地了解種質資源性狀的遺傳相似性、差異性，并根據各類材料的性狀表現和特點進行相應的選擇，可提高藜麥在產業升級、鄉村振興等方面的綜合利用價值。

參考文獻：

［1］Manjarres-Hernández E H，Morillo-Coronado A C. Genetic diversity of colombian quinoa （Chenopodium quinoa Willd.）：implications for breeding programs［J］. Genetic Resources and Crop Evolution，2022，69：2447-2458.

［2］Christensen S A，Pratt D B，Pratt C，et al. Assessment of genetic diversity in the USDA and CIP-FAO international nursery collections of quinoa （Chenopodium quinoa Willd.） using microsatellite markers［J］. Plant Genetic Resources，2007，5（2）：82-95.

［3］貢布扎西，旺 姆，張崇璽，等. 南美藜在西藏的生物學特性表現［J］. 西南農業學報，1994（3）：54-62.

［4］黃 杰，楊發榮. 藜麥在甘肅的研發現狀及前景［J］. 甘肅農業科技，2015（1）：49-52.

［5］劉鎖榮，范文虎. 促進山西藜麥種植規模化及產業鏈形成的建議［J］. 山西農業科學，2011，39（7）：767-769.

［6］劉 洋，熊國富，閆殿海，等. “糧食之母”、“超級食物”——藜麥“落戶”青海［J］. 青海農林科技，2014（4）：95-98.

［7］周 英，魏明鋒. 新疆北疆地區藜麥栽培技術［J］. 新疆農墾科技，2017，40（7）：20-21.

［8］任永峰，王志敏，趙沛義，等. 內蒙古陰山北麓區藜麥生態適應性研究［J］. 作物雜志，2016（2）：79-82，2.

［9］李娜娜，丁漢鳳，郝俊杰，等. 藜麥在中國的適應性種植及發展展望［J］. 中國農學通報，2017，33（10）：31-36.

［10］任貴興，楊修仕，么 楊. 中國藜麥產業現狀［J］. 作物雜志，2015（5）：1-5.

［11］張琴萍，邢 寶，周幫偉，等. 藜麥飼用研究進展與應用前景分析［J］. 中國草地學報，2020，42（2）：162-168.

［12］楊修仕. 第三屆中國藜麥產業（濟南）高峰論壇召開［J］. 作物雜志，2017（5）：2.

［13］首屆中國藜麥產業（長春）高峰論壇在吉林省長春市舉行［J］. 作物雜志，2015（5）：2.

［14］Jacobsen S E，Mujica A，Jense C R. The worldwide potential for quinoa （Chenopodium quinoa Willd.）［J］. Food Reviews International，2003，19（1/2）：167-177.

［15］趙雪雁，馬平易，李文青，等. 黃土高原生態系統服務供需關系的時空變化［J］. 地理學報，2021，76（11）：2780-2796.

［16］王 建，趙牡丹，樊 藝，等. 黃土高原人類活動與生物多樣性的演變及關聯性［J］. 水土保持研究，2022，29（6）：154-160.

［17］魏建洲，牛雪娜，史戰紅，等. 黃土高原草地凈初級生產力變化的驅動因素分析［J］. 環境科學與技術，2022，45（6）：203-211.

［18］楊發榮，黃 杰，魏玉明，等. 藜麥生物學特性及應用［J］. 草業科學，2017，34（3）：607-613.

［19］姜奇彥，牛風娟，胡 正，等. 金藜麥耐鹽性分析及營養評價［J］. 植物遺傳資源學報，2015，16（4）：700-707.

［20］Villa D Y G，Russo L，Kerbab K，et al. Chemical and nutritional characterization of Chenopodium pallidicaule （caihua） and Chenopodium quinoa （quinoa） seeds［J］. Emirates Journal of Food and Agriculture，2014，26 （7）：609-615.

［21］南 銘，趙桂琴，柴繼寬. 黃土高原半干旱區飼用燕麥種質表型性狀遺傳多樣性分析及綜合評價［J］. 草地學報，2017，25（6）：1197-1205.

［22］黃 杰，劉文瑜，楊發榮，等. 不同藜麥品種在東鄉半干旱區的適應性表現［J］. 甘肅農業科技，2022，53（6）：46-50.

［23］王志恒，黃思麒，李成虎，等. 13種藜麥萌發期抗逆性綜合評價［J］. 西北農林科技大學學報（自然科學版），2021，49（1）：25-36.

［24］魏玉明，楊發榮，劉文瑜，等. 隴東旱塬區復種不同藜麥品種（系）的適應性初步評價［J］. 西北農業學報，2020，29（5）：675-686.

［25］梅 麗. 北京藜麥適應性栽培研究進展及展望［J］. 作物雜志，2022（6）：14-22.

［26］王 元，張忠福，王建梅. 山丹縣藜麥引種試驗總結［J］. 農業科技與信息，2018，（7）：11-12.

［27］徐冬麗，周蘭蘭，王國平，等. 高原藜麥新品種（系）引進比較試驗初報［J］. 西藏農業科技，2021，43（1）：29-32.

［28］陳正華，王 建. 利用遙感技術建立干旱半干旱地區草地生態健康模型［J］. 遙感技術與應用，2005（6）：558-562.

［29］李良斌，王 耀，雷成軍，等. 寒旱山區藜麥引種試驗初報［J］. 甘肅農業科技，2021，52（4）：64-67.

［30］孫夢涵，邢 寶，崔宏亮，等. 藜麥種質資源遺傳多樣性SSR標記分析［J］. 植物遺傳資源學報，2021，22（3）：625-637.

［31］石振興. 國內外藜麥品質分析及其減肥活性研究［D］. 北京：中國農業科學院，2016.

［32］趙小雪，郭鳳根，羅富成. 貯藏時間及纖維素酶對臺灣紅藜青貯品質的影響［J］. 飼料研究，2022，45（10）：76-80.

［33］喻 泓. 門頭溝區臺灣紅藜引種栽培技術［J］. 農業科技通訊，2022（1）：257-258，271.

［34］鄧 杰，趙 鋼，徐漪沙，等. 四川地區2種藜麥營養成分的比較分析［J］. 糧食與油脂，2021，34（2）：43-46.

［35］Curti R N，Andrade A J，Bramardi S，et al. Ecogeographic structure of phenotypic diversity in cultivated populations of quinoa from Northwest Argentina［J］. Annals of Applied Biology，2012，160（2）：114-125.

［36］王思宇，左文博，朱凱莉，等. 71份藜麥品種資源的農藝性狀及營養品質分析與評價［J］. 作物雜志，2022（3）：63-72.

［37］王艷青，李勇軍，李春花，等. 藜麥主要農藝性狀與單株產量的相關和通徑分析［J］. 作物雜志，2019（6）：156-161.

［38］韓 霄. 從作物的源流庫理論展望新型育種技術［J］. 生物技術通報，2015，31（4）：34-39.

［39］邱正高，楊 華，祁志云，等. 航空誘變的糯玉米矮稈突變體株高與葉片生長分析［J］. 西南農業學報，2009，22（4）：901-904.

［40］胡一波，楊修仕，陸 平，等. 中國北部藜麥品質性狀的多樣性和相關性分析［J］. 作物學報，2017，43（3）：464-470.[HJ]

［41］顏小捷，谷陟欣，盧鳳來，等. 裸花紫珠中總酚、總黃酮和總皂苷的含量測定［J］. 時珍國醫國藥，2013，24（12）：2852-2854.

［42］王麗娜，任翠梅，王明澤，等. 中國藜麥種質資源分布及研究現狀［J］. 黑龍江農業科學，2020（12）：142-145.

［43］任永峰，梅 麗，楊亞東，等. 播期對藜麥農藝性狀及產量的影響［J］. 中國生態農業學報，2018，26（5）：643-656.

［44］張亞萍，王致和，張秀華，等. 14個藜麥品種（系）在祁連山區的農藝性狀表現及其與產量的關系分析［J］. 山東農業科學，2021，53（8）：17-21.

［45］黃 杰，劉文瑜，呂 瑋，等. 38份藜麥種質資源農藝性狀與產量的關系分析［J］. 甘肅農業科技，2018（12）：72-75.

［46］黃 杰，楊發榮，李敏權，等. 13個藜麥材料在甘肅臨夏旱作區適應性的初步評價［J］. 草業學報，2016，25（3）：191-201.

［47］王艷青，李春花，盧文潔，等. 135份國外藜麥種質主要農藝性狀的遺傳多樣性分析［J］. 植物遺傳資源學報，2018，19（5）：887-894.