基于隸屬函數(shù)法和GGE雙標圖的飼草型小黑麥種質(zhì)適應性評價

劉晶, 宋謙, 田新會, 杜文華*, 劉漢成

(1.甘肅農(nóng)業(yè)大學草業(yè)學院，草業(yè)生態(tài)系統(tǒng)教育部重點實驗室，甘肅省草業(yè)工程實驗室，甘肅 蘭州 730070；2.青海民族大學生態(tài)環(huán)境與資源學院，青海 西寧 810000；3.甘肅省平?jīng)鍪修r(nóng)業(yè)科學院，甘肅 平?jīng)?744000)

小黑麥(Triticale)是小麥屬(Triticum)和黑麥屬(Secale)植物經(jīng)屬間有性雜交，通過染色體加倍及染色體工程育種得到的第1個人工育成的新物種，為一年生禾本科植物[1]。人工培育的小黑麥新品種通常需要進行區(qū)域試驗，以研究其在不同區(qū)域的草產(chǎn)量、營養(yǎng)品質(zhì)、抗旱性和抗寒性等，以確定最適宜種植區(qū)域，并發(fā)揮其高產(chǎn)潛力[2-4]。

目前國內(nèi)外有關小黑麥草產(chǎn)量和營養(yǎng)價值的研究主要集中在與豆科植物混播[5]、水肥對小黑麥草產(chǎn)量與品質(zhì)的影響[6-10]、播種密度和氮素營養(yǎng)對小黑麥草產(chǎn)量和品質(zhì)的影響[11]、小黑麥不同基因型的飼草品質(zhì)[12]、刈割時期對小黑麥草產(chǎn)量和品質(zhì)的影響[13-14]等方面，營養(yǎng)價值研究只分析粗蛋白、粗脂肪和中(酸)性洗滌纖維含量等單個指標[9-10]。模糊隸屬函數(shù)值是一個表示相對優(yōu)劣的綜合指數(shù)，可以反映多個構(gòu)成成分的整體水平及綜合特征，是用來評估植物營養(yǎng)品質(zhì)的有效方法[15-19]。目前尚少有應用隸屬函數(shù)法對飼草型小黑麥營養(yǎng)品質(zhì)進行評價的研究報道。GGE(基因型和基因與環(huán)境互作效應，genotype main effects plus genotype×environment interaction)雙標圖法是農(nóng)作物品種區(qū)域試驗中進行品種評價、環(huán)境評價和品種生態(tài)區(qū)劃最高效和最直觀的統(tǒng)計和圖形展示方法，在多年多點區(qū)域試驗中，此模型是一個十分理想的分析工具[20-22]。GGE雙標圖可以揭示不同因素間錯綜復雜的相互作用關系[20]，已經(jīng)廣泛用于分析作物產(chǎn)量和品質(zhì)數(shù)據(jù)[23]。柴守璽等[24]利用GGE雙標圖法研究了飼料型小黑麥基因型與環(huán)境互作效應以及種子產(chǎn)量的穩(wěn)定性。目前尚未有應用GGE雙標圖對飼草型小黑麥草產(chǎn)量和營養(yǎng)品質(zhì)進行綜合評價的報道。

本研究首先對2014-2015年4份小黑麥種質(zhì)在甘肅省不同區(qū)域的營養(yǎng)價值(粗蛋白含量，NDF含量，ADF含量)和干物質(zhì)消化率進行方差分析和多重比較，然后利用隸屬函數(shù)法計算4份小黑麥種質(zhì)及不同試驗點營養(yǎng)品質(zhì)(粗蛋白含量，NDF含量，ADF含量，干物質(zhì)消化率)的隸屬度值，并根據(jù)指標的重要性對營養(yǎng)品質(zhì)指標賦予權(quán)重，以此作為小黑麥營養(yǎng)品質(zhì)的綜合評價標準，最后利用GGE雙標圖研究小黑麥基因型與環(huán)境互作效應的穩(wěn)定性，以篩選甘肅省不同地區(qū)適宜種植的小黑麥的品種(系)和最適合種植小黑麥的區(qū)域，從而為飼草型小黑麥種質(zhì)在甘肅省的示范推廣提供科學依據(jù)。

1 材料與方法

1.1 試驗地概況

試驗于2014和2015年在甘肅省合作市、臨洮縣、瑪曲縣、夏河縣、肅南縣馬蹄鄉(xiāng)和康樂鄉(xiāng)進行。各試驗點的海拔、經(jīng)度、緯度、降水量和土壤類型見表1。

1.2 試驗材料

供試材料為甘肅農(nóng)業(yè)大學草業(yè)學院培育的飼草型小黑麥新品系P2(G1)和P4(G2)，石河子大學農(nóng)學院培育的飼草型小黑麥品種石大1號(G3)，中國農(nóng)業(yè)科學院作物科學研究所培育的飼草型小黑麥品種中飼1048(G4)。

1.3 試驗設計及田間管理

本試驗包含3個因素，分別為小黑麥種質(zhì)、試驗點和年份。小黑麥種質(zhì)有4個水平：P2(G1)，P4(G2)，石大1號(G3)，中飼1048(G4)；試驗點有6個水平：臨洮(E1)，合作(E2)，瑪曲(E3)，肅南康樂鄉(xiāng)(E4)，夏河(E5)，肅南馬蹄鄉(xiāng)(E6)；年份有2個水平：2014年(Y1)，2015年(Y2)。

在2年6個試驗點的試驗中，每個試驗點的試驗設計相同。隨機區(qū)組試驗設計，條播，行距30 cm，播種深度3～4 cm，播種量按照750萬基本苗·hm-2計算而得。4次重復，小區(qū)面積5 m×3 m=15 m2，試驗地周圍1 m種植保護行。各個試驗點的其他管理同于當?shù)卮筇锷a(chǎn)。6個試點中臨洮點有灌溉條件，分別于小黑麥出苗期和拔節(jié)期進行施肥(施肥量均為150 kg N·hm-2)和灌水(滿灌)，其他試點均無灌溉條件，為雨養(yǎng)區(qū)，試驗期間未灌水，施肥量和施肥時間同于臨洮點。

表1 試驗地概況Table 1 Introduction of the experimental sites

1.4 測定指標及方法

草產(chǎn)量：開花期進行[14]。測產(chǎn)時，齊地面刈割每個小區(qū)內(nèi)所有植株的地上部(除去邊行和地頭兩邊的50 cm部分)，稱重，得到鮮草產(chǎn)量。同時取鮮草樣500 g，在65～70 ℃烘箱中烘48 h，恒重，稱重得到干草產(chǎn)量，計算鮮干比。根據(jù)鮮干比計算每個小區(qū)的干草產(chǎn)量。

營養(yǎng)品質(zhì)：用粉碎機粉碎烘干后的草樣，過1 mm篩子，從混合均勻的草樣中隨機取4份樣品，平行測定各項指標。粗蛋白(crude protein，CP)含量測定采用凱氏定氮法[25]，中性洗滌纖維(neutral detergent fiber，NDF)和酸性洗滌纖維(acid detergent fiber，ADF)含量測定采用范氏的洗滌纖維分析法[25]，干物質(zhì)消化率(dry matter digestibility，DMD)采用人工瘤胃法[25]。

1.5 小黑麥種質(zhì)和試驗點營養(yǎng)品質(zhì)綜合評價

用隸屬函數(shù)法對4個小黑麥種質(zhì)和6個試點的營養(yǎng)品質(zhì)進行綜合評價。隸屬函數(shù)值[X(μ1)，X(μ2)]計算公式為：

X(μ1)=(X-Xmin)/(Xmax-Xmin)

(1)

X(μ2)=1-X(μ1)

(2)

式中：X為某一小黑麥種質(zhì)(或試點)某一指標的測定值；Xmax為所有小黑麥種質(zhì)(或所有試點)該指標的最大值；Xmin為所有小黑麥種質(zhì)(或試點)該指標的最小值。若所測指標與小黑麥品質(zhì)呈正相關，則采用(1)式計算隸屬值，負相關則用(2)式。

小黑麥營養(yǎng)品質(zhì)綜合評價時，根據(jù)營養(yǎng)價值和干物質(zhì)消化率在小黑麥生產(chǎn)中的重要程度[26]，各指標的權(quán)重分配如表2。根據(jù)表2的權(quán)重，利用公式(3)計算每個小黑麥種質(zhì)在每個試點的營養(yǎng)評價值。

(3)

式中：i代表試點，1代表臨洮，2代表合作，3代表瑪曲，4代表肅南康樂，5代表夏河，6代表肅南馬蹄；j代表小黑麥種質(zhì)，1代表P2，2代表P4，3代表石大1號，4代表中飼1048；k代表小黑麥營養(yǎng)品質(zhì)，1代表CP，2代表ADF，3代表NDF，4代表DMD；rij代表第i個試點第j個小黑麥種質(zhì)營養(yǎng)評價值；ξijk代表第i個試點第j個小黑麥材料第k個營養(yǎng)品質(zhì)指標對應的隸屬函數(shù)值；Wk代表營養(yǎng)品質(zhì)指標權(quán)重。

1.6 數(shù)據(jù)處理

采用Excel對小黑麥產(chǎn)量，營養(yǎng)數(shù)據(jù)，干物質(zhì)消化率進行整理，利用SPSS 19統(tǒng)計分析軟件[27]進行方差分析和多重比較，基因與環(huán)境互作采用Yan等[23]開發(fā)的GGEbiplot 6.3分析軟件。

GGE雙標圖將多點試驗得到的小黑麥種質(zhì)的產(chǎn)量性狀和品質(zhì)得分結(jié)果整理成一個小黑麥種質(zhì)-試點兩向表。每個數(shù)值是相應小黑麥種質(zhì)在相應試驗點的性狀平均值，即為表型值(P)。P的計算如下:

表2 小黑麥綜合營養(yǎng)品質(zhì)的構(gòu)建及各指標的權(quán)重Table 2 Construction on the nutrition quality of the triticale and the weight of each index

P=M+E+G+GE

(4)

式中：M為多年多點試驗性狀總平均值；E為環(huán)境主效應；G為小黑麥種質(zhì)主效應；GE為種質(zhì)-環(huán)境互作效應。通過試驗數(shù)據(jù)環(huán)境中心化，即兩向表的每個數(shù)據(jù)減去相應環(huán)境的平均值，即去掉M和E，形成環(huán)境中心化的兩向表?；诃h(huán)境中心化的小黑麥種質(zhì)-環(huán)境兩向表的雙標圖只含有與種質(zhì)評價有關的G和GE，稱GGE雙標圖[20-21]。根據(jù)作物品種或試驗點在AT軸(Average-Tester axis)上的投影長度判斷產(chǎn)量和營養(yǎng)品質(zhì)的穩(wěn)定性和試點的鑒別力。根據(jù)品種或試驗點在AT軸上的投影位置判斷參試產(chǎn)量和營養(yǎng)品質(zhì)的平均表現(xiàn)以及試驗點的代表性[21]。雙標圖解釋中所隱含的假設，就是2-D雙標圖可以充分近似其代表的兩向表數(shù)據(jù)。但是，如何知道這個假設是否滿足，首先要看雙標圖的擬合度，即前兩個主成分(PC1和PC2)所解釋的兩向表中總變異的百分數(shù)。如果擬合度較高，則雙標圖能較好地近似實際數(shù)據(jù)[20-21]。

2 結(jié)果與分析

2.1 干草產(chǎn)量

用Kolmogorov-Smirnov正態(tài)性檢驗[27]小黑麥種質(zhì)的干草產(chǎn)量可知，P=0.34>0.05，符合正態(tài)分布，適合于進行方差分析。方差分析(表3)表明，單因素(試驗點，種質(zhì))間、試點與種質(zhì)交互作用間干草產(chǎn)量均存在極顯著差異，需要進行多重比較。根據(jù)F值大小，試點對小黑麥干草產(chǎn)量的影響最大，試點×種質(zhì)互作對干草產(chǎn)量影響次之。

表3 參試小黑麥種質(zhì)干草產(chǎn)量，CP，ADF，NDF和DMD的方差分析Table 3 Variance analysis on the hay yield, content of CP, ADF and NDF, and DMD for triticale genotypes

2.1.1單因素間小黑麥干草產(chǎn)量的差異 試點間多重比較(表4)表明，6個試點中，E6點小黑麥干草產(chǎn)量最高(14.23 t·hm-2)，除與E1和E3無顯著差異外，與其他3個試點存在顯著差異；E5點干草產(chǎn)量最低(3.45 t·hm-2)，且與其他5個試點均存在顯著差異。在干草產(chǎn)量較高的3個試點中，根據(jù)標準誤大小，E1點的小黑麥干草產(chǎn)量高且穩(wěn)定性最好。

種質(zhì)間多重比較結(jié)果(表4)可知，參試小黑麥種質(zhì)中，G1的干草產(chǎn)量(12.83 t·hm-2)顯著高于其他種質(zhì)；G2的干草產(chǎn)量(10.76 t·hm-2)次之；G4的干草產(chǎn)量(8.90 t·hm-2)顯著低于其他處理。

2.1.2試點×種質(zhì)作用間小黑麥干草產(chǎn)量的差異 試點×種質(zhì)互作效應的多重比較(表5)表明，E3G1的干草產(chǎn)量最高(18.78 t·hm-2)，除與E6G1無顯著差異外，與其他處理間均有顯著差異(P<0.05)。E5G4的干草產(chǎn)量最低(2.63 t·hm-2)，其產(chǎn)量僅為E3G1處理的15%，該處理除與E5G1、E5G2、E5G3、E2G3、E2G4處理無顯著差異外，與其他處理均存在顯著差異，所以E5點最不適合種植小黑麥。

2.2 營養(yǎng)品質(zhì)

CP、NDF、ADF和DMD測定結(jié)果經(jīng)Kolmogorov-Smirnov正態(tài)性檢驗可得，小黑麥干草營養(yǎng)品質(zhì)符合正態(tài)分布，適宜進行方差分析。方差分析(表3)可得，試點間小黑麥的營養(yǎng)品質(zhì)均存在極顯著差異，需要進行多重比較。根據(jù)F值大小，小黑麥種質(zhì)間CP含量的差異最大，試點×種質(zhì)互作效應間CP含量的差異次之；試點對小黑麥干草NDF含量的影響最大，種質(zhì)對NDF含量的影響次之；試點對小黑麥干草ADF含量的影響最大，種質(zhì)對ADF含量的影響次之；試點對小黑麥干草DMD影響最大，試點×種質(zhì)互作效應對DMD的影響次之。

表4 單因素間小黑麥干草產(chǎn)量，CP，ADF，NDF含量和DMD的多重比較和營養(yǎng)評價值Table 4 Multiple comparisons on the hay yield, content of CP, ADF and NDF, and DMD in triticale for single factor and the value of nutritional evaluation

2.2.1單因素間小黑麥營養(yǎng)品質(zhì)的差異 試點間多重比較(表4)可知，6個試點中，E3和E4點小黑麥干草的CP含量顯著高于其他試點；E5點小黑麥干草的NDF(52.36%)和ADF(37.15%)含量顯著低于其他試點；E3點小黑麥干草的DMD(59.83%)顯著高于其他試點。

種質(zhì)間多重比較(表4)可知，參試小黑麥種質(zhì)中，G1的CP含量(11.90%)顯著高于其他材料，NDF(57.57%)和ADF(38.87%)含量顯著低于G2和G3，DMD含量(53.09%)和G4相近，顯著高于G2和G3，G3的各項營養(yǎng)品質(zhì)的排名均靠后。

2.2.2試點×種質(zhì)交互作用間小黑麥營養(yǎng)品質(zhì)指標的差異 試點×種質(zhì)互作效應的多重比較(表5)表明，E3G1的CP含量為13.92%，除與E4G1干草CP含量無顯著差異之外，顯著高于其他處理的CP含量(P<0.05)。E5G1的NDF含量最低(50.26%)，與E6G1、E1G3、E2G4、E4G1、E2G3、E4G4、E4G2、E4G3有顯著性差異外，和其他處理均無顯著性差異。E5G1干草ADF含量最低為36.10%，與E4G1、E4G4、E4G2、E4G3有顯著性差異外，和其他處理均無顯著性差異。E5G1干草消化率最高，為65.36%，顯著高于其他處理，E2G4干草消化率為36.66%，顯著低于其他處理，所以G1干草CP、DMD含量顯著高于其他材料；G1干草NDF含量和ADF含量顯著低于G2、G3，這與2.2.1結(jié)論基本一致。

表5 試點×種質(zhì)交互作用的小黑麥干草產(chǎn)量，CP含量，ADF含量，NDF含量和DMD的多重比較和營養(yǎng)評價值Table 5 Multiple comparisons on the hay yield, content of CP, ADF and NDF, and DMD in triticale for the interactions of experiment site and genotype and the value of nutritional evaluation

2.3 綜合評價

由于本試驗參試小黑麥種質(zhì)的3個營養(yǎng)價值指標和干物質(zhì)消化率表現(xiàn)不一致，無法對營養(yǎng)品質(zhì)進行精準評價。隸屬函數(shù)法可以消除個別指標帶來的片面性，使各小黑麥種質(zhì)營養(yǎng)品質(zhì)的差異具有真實可比性，所以本研究利用隸屬函數(shù)法對參試小黑麥種質(zhì)的營養(yǎng)品質(zhì)進行了綜合評價，結(jié)果見表4和表5。

通過對小黑麥種質(zhì)營養(yǎng)品質(zhì)綜合評價可知，6個試點中E1的營養(yǎng)評價值最高(0.51)，E3次之，為0.50(表4)，4個種質(zhì)中G1的營養(yǎng)評價值最高，為0.67，這與2.2.1結(jié)論基本一致。

試點×種質(zhì)互作效應營養(yǎng)評價值較高的處理有：E5G1>E2G1>E1G1>E3G2>E3G1=E6G4，這與2.2.1和2.2.2的結(jié)論基本一致。

綜合干草產(chǎn)量和營養(yǎng)品質(zhì)指標的方差分析和多重比較以及營養(yǎng)評價值可得： 1)6個試點中，E6、E3和E1點小黑麥干草產(chǎn)量最高，E5點的干草產(chǎn)量最低；由營養(yǎng)價值和干物質(zhì)消化率的多重比較及營養(yǎng)評價值得，E1、E3和E5點營養(yǎng)品質(zhì)綜合較好；2)4個參試小黑麥種質(zhì)中，G1的干草產(chǎn)量和營養(yǎng)品質(zhì)最好?？捎肎GE雙標圖對上述結(jié)果進行驗證。

2.4 小黑麥種質(zhì)干草產(chǎn)量和營養(yǎng)價值的平均表現(xiàn)及穩(wěn)定性

2.4.1干草產(chǎn)量的平均表現(xiàn)和穩(wěn)定性 在特定生態(tài)區(qū)域，理想品種應當既高產(chǎn)又穩(wěn)產(chǎn)且品質(zhì)較優(yōu)。方差分析和多重比較通常通過比較平均草產(chǎn)量篩選高產(chǎn)品種，根據(jù)標準差篩選穩(wěn)產(chǎn)品種。兼具“高產(chǎn)性和穩(wěn)產(chǎn)性”功能的GGE雙標圖可以同時了解品種的高產(chǎn)性和穩(wěn)產(chǎn)性[21]。本試驗GGE雙標圖(圖1)中的小圓圈代表“平均環(huán)境”；單箭頭直線表示平均環(huán)境軸，其所指方向是參試小黑麥種質(zhì)在所有環(huán)境下的近似平均草產(chǎn)量，采用聚焦品種的特征值分配(singular value partitioning，SVP=1)[21]。由圖1可知，G1的平均草產(chǎn)量最高，G4的平均草產(chǎn)量最低，這與2.1.1的結(jié)論一致。與平均環(huán)境軸垂直并通過原點、帶有雙箭頭的直線代表各小黑麥種質(zhì)與各試點相互作用的傾向性：箭頭向外指向的線段越短，說明草產(chǎn)量越穩(wěn)定。圖1中G3和G1的草產(chǎn)量較穩(wěn)定，草產(chǎn)量最不穩(wěn)定的是G2。綜合參試小黑麥種質(zhì)草產(chǎn)量的平均表現(xiàn)和穩(wěn)定性可知，小黑麥品系G1既高產(chǎn)又穩(wěn)產(chǎn)。另外，圖1中主成分(principal components, PC)PC1和PC2集中了小黑麥種質(zhì)主效應+種質(zhì)-環(huán)境互作效應(G+GE)的大部分變異信息，擬合度為90.8%，說明利用雙標圖得到的草產(chǎn)量結(jié)果比較接近真實狀況，具有較大可靠性。

2.4.2營養(yǎng)品質(zhì)的平均表現(xiàn)和穩(wěn)定性 由圖2可知，參試小黑麥種質(zhì)營養(yǎng)品質(zhì)由高到低依次為：G1>G2>G4>G3，穩(wěn)定性分別為G3>G1>G2>G4。綜合營養(yǎng)品質(zhì)和穩(wěn)定性可知，G1的營養(yǎng)品質(zhì)好而且穩(wěn)定；G4和G2的營養(yǎng)品質(zhì)相近，但G2的穩(wěn)定性較高；G3的穩(wěn)定性雖然好，但營養(yǎng)品質(zhì)最差。從圖2也可知，主成分PC1和PC2集中了G+GE的大部分變異信息(94.5%)，據(jù)此分析得到的結(jié)果可靠性較高。

圖1 基于GGE-biplot分析的參試小黑麥種質(zhì)的干草產(chǎn)量及其穩(wěn)定性Fig.1 Hay yields and their stability for different triticale genotypes based on the GGE-biplot analysis

圖2 基于GGE-biplot分析的參試小黑麥種質(zhì)營養(yǎng)品質(zhì)的表現(xiàn)及其穩(wěn)定性Fig.2 Nutritional quality and their stability for different triticale genotypes based on the GGE-biplot analysis

G表示參試小黑麥種質(zhì)，G1: P2，G2: P4，G3: 石大1號，G4: 中飼1048；E代表試驗點，E1：臨洮， E2：合作， E3：瑪曲，E4：肅南康樂鄉(xiāng)，E5：夏河，E6：肅南馬蹄鄉(xiāng)。下同。
G stands for triticale genotypes, G1: P2, G2: P4, G3: Shida No.1, G4: Zhongsi 1048, and E stands for the experimental sites, E1: Lintao, E2: Hezuo, E3: Maqu, E4: Kangle, Sunan, E5: Xiahe, E6: Mati, Sunan. The same below.

2.5 小黑麥種質(zhì)最佳種植區(qū)域的篩選

2.5.1干草產(chǎn)量 為了篩選各試驗點干草產(chǎn)量均較高的小黑麥種質(zhì)，把圖1中各個方向上距離最遠的點用直線連接起來，這樣，G4、G2和G1就構(gòu)成了1個三角形(圖3)，采用聚焦環(huán)境的特征值分配(singular value partitioning，SVP=2)[21]。通過三角形中心分別向3條邊做3條垂線，發(fā)現(xiàn)小黑麥種質(zhì)在3個扇區(qū)都有分布，而試點分布在2個扇區(qū)：E2和E5在同一扇區(qū)，其他4個試點(E1，E3，E4，E6)在同一扇區(qū)。說明G1在E1、E3、E4和E6的干草產(chǎn)量均較高，G2在E2和E5的草產(chǎn)量較高，G4在所有試點草產(chǎn)量表現(xiàn)均不理想，而位于多邊形內(nèi)部、靠近原點的G3的干草產(chǎn)量對試點變化不敏感。

2.5.2營養(yǎng)品質(zhì) 為了篩選各試驗點營養(yǎng)品質(zhì)均較優(yōu)的小黑麥種質(zhì)，把圖2中各個方向上距離最遠的點用直線連接起來，發(fā)現(xiàn)G1、G2、G3和G4構(gòu)成了1個四邊形(圖4)。通過中心分別向4條邊做4條垂線，將雙標圖分為4個扇區(qū)，小黑麥種質(zhì)在4個扇區(qū)均有分布，試點分布在2個扇區(qū)。位于試點扇區(qū)頂點的G1在扇內(nèi)E1、E2、E3、E4和E5點的營養(yǎng)品質(zhì)均較高，G4在扇內(nèi)E6的營養(yǎng)品質(zhì)較高，G2和G3的營養(yǎng)品質(zhì)在所有試點均不理想。

圖3 基于GGE雙標圖分析小黑麥種質(zhì)與試點互作效應的干草產(chǎn)量表現(xiàn)Fig.3 Hay yield for the interaction between the triticale genotype and experimental sites based on the GGE-biplot analysis

圖4 基于GGE雙標圖分析小黑麥種質(zhì)與試點互作效應的營養(yǎng)品質(zhì)表現(xiàn)Fig.4 Nutrition quality for the interaction between the triticale genotype and experimental sites based on the GGE-biplot analysis

2.6 小黑麥種質(zhì)與環(huán)境的優(yōu)選

利用GGE雙標圖可以方便地確定出一個理想品種(環(huán)境)的位置[20-24]。理想品種指在所有試點的草產(chǎn)量最高、營養(yǎng)品質(zhì)最好，理想環(huán)境則指對所有參試品種分辨能力最強，草產(chǎn)量和營養(yǎng)品質(zhì)在該環(huán)境下最高[28]。理想品種被定義為平均環(huán)境軸正方向上的一個點(絕對穩(wěn)定)，該點為坐標正方向上的一個點(最具有代表性)。該點距坐標原點的長度則為最長環(huán)境向量的長度(辨別力最強)?！袄硐搿逼贩N(或環(huán)境)事實上并不存在，但可作為品種(環(huán)境)理想程度比較的參照。以理想品種(或環(huán)境)為圓心做多層同心圓，根據(jù)參試小黑麥種質(zhì)(試點)與理想品種(或環(huán)境)的接近程度，可直觀地對供試種質(zhì)(或試點)優(yōu)劣進行排序。越靠近同心圓中心，則表示該品種(或環(huán)境)越理想，反之亦然。

圖5 參試小黑麥種質(zhì)與理想品種的比較Fig.5 Comparisons of the tested triticale genotypes with the ideal variety

根據(jù)圖5可得，G1距離兩個圓心均最近，說明其干草產(chǎn)量最高、營養(yǎng)品質(zhì)最好；G2次之，其干草產(chǎn)量較高、營養(yǎng)品質(zhì)次之；G3和G4的草產(chǎn)量和營養(yǎng)品質(zhì)較差。采用聚焦品種的特征值分配(SVP=1)。

根據(jù)圖6可得，6個試驗點小黑麥干草產(chǎn)量離圓心的距離由近到遠依次為：E6>E1>E3>E2>E5>E4，營養(yǎng)品質(zhì)離圓心的距離由近到遠依次為：E3>E1>E5>E4>E2>E6。E6點雖然干草產(chǎn)量最高，但營養(yǎng)品質(zhì)綜合評價值最低；E1點的干草產(chǎn)量和營養(yǎng)品質(zhì)均排名第2；E3點的干草產(chǎn)量排名第3，營養(yǎng)品質(zhì)排名第1。綜合考慮6個試驗點小黑麥種質(zhì)的草產(chǎn)量和營養(yǎng)品質(zhì)，E3和E1為最適合種植小黑麥的試點。采用聚焦環(huán)境的特征值分配(SVP=2)[21]。

圖6 6個試點與理想試點的比較Fig.6 Comparisons of the tested experimental sites with the ideal site

3 討論

隨著我國草食畜牧業(yè)不斷發(fā)展，飼草供求不平衡、家畜日糧不足等問題日益凸顯。我國天然草原由于長期超載過牧，生產(chǎn)力急劇下降，已不能滿足草食畜牧業(yè)可持續(xù)發(fā)展的要求[29]。為了畜牧業(yè)穩(wěn)步健康發(fā)展，尋找優(yōu)質(zhì)人工牧草是解決草食畜牧業(yè)可持續(xù)發(fā)展的重要措施。本研究將方差分析、隸屬函數(shù)法和GGE雙標圖法相結(jié)合，篩選適合于甘肅省不同氣候區(qū)域種植的草產(chǎn)量較高、營養(yǎng)品質(zhì)較好的小黑麥種質(zhì)，可以為解決牧草短缺問題提供理論依據(jù)。

3.1 高產(chǎn)優(yōu)質(zhì)小黑麥種質(zhì)的篩選

由于年份是不可重復且不可控制的隨機因素，所有與年份互作分析也都是沒有意義的，所以只需對試點、種質(zhì)、試點×種質(zhì)進行方差分析。根據(jù)參試小黑麥種質(zhì)在不同試點的干草產(chǎn)量和營養(yǎng)品質(zhì)多重比較，品系P2的干草產(chǎn)量最高、營養(yǎng)品質(zhì)最優(yōu)，具有推廣價值；P4的干草產(chǎn)量和營養(yǎng)品質(zhì)較優(yōu)；中飼1048干草產(chǎn)量最低，石大1號的營養(yǎng)品質(zhì)最差。根據(jù)GGE雙標圖結(jié)果，小黑麥品系P2的干草產(chǎn)量最高、營養(yǎng)品質(zhì)最好，穩(wěn)定性最好；品系P4的干草產(chǎn)量較高、營養(yǎng)品質(zhì)較高，穩(wěn)定性較好；中飼1048小黑麥的營養(yǎng)品質(zhì)較低、干草產(chǎn)量最低，且穩(wěn)定性也較差；石大1號小黑麥的草產(chǎn)量雖高，但營養(yǎng)評價值和穩(wěn)定性最差。通過分析可知，參試小黑麥種質(zhì)和外界條件以及它們的互作是造成草產(chǎn)量和營養(yǎng)品質(zhì)存在顯著差異的主要原因。在甘肅中南部和高寒牧區(qū)，銹病是危害小黑麥生長發(fā)育的主要病害[30]。小黑麥感染銹病后對草產(chǎn)量和飼草營養(yǎng)品質(zhì)，尤其是粗蛋白含量和干物質(zhì)消化率影響較大[31-32]。石大1號小黑麥雖然草產(chǎn)量中等，但由于其高感銹病，植株感染銹病后，葉片枯黃，對營養(yǎng)品質(zhì)影響較大。中飼1048小黑麥本身草產(chǎn)量較低，再加上其高感銹病，所以在4個種質(zhì)中表現(xiàn)最差。P2和P4是甘肅農(nóng)業(yè)大學草業(yè)學院以高產(chǎn)優(yōu)質(zhì)和抗銹病為育種目標培育的小黑麥新品系，P4雖然感病，但感病率較低，對草產(chǎn)量和營養(yǎng)品質(zhì)的影響較?。籔2由于分蘗性能強，CP含量高，而且高抗銹病，所以其草產(chǎn)量和營養(yǎng)品質(zhì)均最佳。

3.2 小黑麥種植區(qū)域的篩選

正確選擇某一作物的種植區(qū)域?qū)Πl(fā)揮其生產(chǎn)潛力具有重要意義。根據(jù)不同試點參試小黑麥種質(zhì)的干草產(chǎn)量和營養(yǎng)品質(zhì)，臨洮和瑪曲點小黑麥的干草產(chǎn)量和營養(yǎng)評價值明顯優(yōu)于其他試點。根據(jù)GGE雙標圖小黑麥種質(zhì)與環(huán)境的優(yōu)選結(jié)果，臨洮和瑪曲點是小黑麥最理想的生存環(huán)境，這與宋謙[32]和趙丹[33]的研究結(jié)果一致。這主要是因為，小黑麥性喜冷涼濕潤氣候條件，瑪曲點降水量較多(表1)，而且主要分布在6-9月[34]，和小黑麥的生長節(jié)律相吻合。臨洮點的降水量雖然居中，但其有灌水條件，可以滿足小黑麥對水分的需求。

4 結(jié)論

參試的4個小黑麥種質(zhì)中，小黑麥新品系P2在甘肅中部和高寒牧區(qū)的草產(chǎn)量最高而且穩(wěn)定、營養(yǎng)品質(zhì)最優(yōu)，有待于進一步示范推廣。

參考文獻References：

[1] Sun M, Guo Y. Biology property, nutritive value and use foreground of triticale. Journal of Shanxi Agricultural University, 2003, (3): 200-204.

孫敏, 郭媛. 小黑麥生物學特性、營養(yǎng)價值及利用前景. 山西農(nóng)業(yè)大學學報, 2003, (3): 200-204.

[2] Chen L X, Tian X H, Du W H. Studies on the seed yield and yield components of new triticale line in the winter wheat grown area of Longzhong Region, Gansu Province. Grassland and Turf, 2015, 35(5): 70-74.

陳麗霞, 田新會, 杜文華. 小黑麥新品系在甘肅省隴中冬小麥區(qū)的種子產(chǎn)量及產(chǎn)量構(gòu)成因素的研究. 草原與草坪, 2015, 35(5): 70-74.

[3] Dong Y Q, Cheng Z J. Breeding of Qianzhong 7, a new allo-octoploid triticale variety. Guizhou Agricultural Sciences, 2005, 33(5): 7-9.

董永琴, 程志軍. 異源八倍體小黑麥黔中7號新品種的選育. 貴州農(nóng)業(yè)科學, 2005, 33(5): 7-9.

[4] Kong G C. New type high yield and good quality forage triticale breeding and application of No.2. Seed, 2005, 24(11): 76-77.

孔廣超. 高產(chǎn)優(yōu)質(zhì)飼草型小黑麥新小黑麥2號的選育與應用. 種子, 2005, 24(11): 76-77.

[5] Vasiljevi? D M S, Karagi? D, Mihailovi? V,etal. Breeding in a world of scarcity. Switzerland Cham: Springer, 2016: 215-218.

[6] Verloop J, Hilhorst G J, Oenema J,etal. Soil N mineralization in a dairy production system with grass and forage crops. Nutrient Cycling in Agroecosystems, 2014, 98(3): 267-280.

[7] Li X N, Wu J Y, Sun W Y.etal. Influence of reclaimed water irrigation on forage triticale quality. Journal of Triticeae Crops, 2012, 32(3): 460-464.

李曉娜, 武菊英, 孫文元, 等. 再生水灌溉對飼用小黑麥品質(zhì)的影響. 麥類作物學報, 2012, 32(3): 460-464.

[8] Dong Z R, Tian L Z, Zhao B,etal. Effects of nitrogen application on the yield and quality ofTriticalecereal. Pratacultural Science, 2008, 25(5): 64-67.

董召榮, 田靈芝, 趙波, 等. 小黑麥牧草產(chǎn)量與品質(zhì)對施氮的響應. 草業(yè)科學, 2008, 25(5): 64-67.

[9] Wang C H, Du X K, Zhang D,etal. Effect of nitrogen fertilization technology on yield and quality ofTriticalecereale. Journal of Northeast Agricultural University, 2012, 43(11): 20-23.

王春宏, 杜曉坤, 張迪, 等. 氮肥施用技術對小黑麥產(chǎn)量及飼用品質(zhì)的影響. 東北農(nóng)業(yè)大學學報, 2012, 43(11): 20-23.

[10] Li J, Zu W, Ji B,etal. Effect of usage of nitrogen on group development and feeding quality of triticale Dongnong 96026. Chinese Agricultural Science Bulletin, 2009, 25(7): 141-144.

李晶, 祖?zhèn)? 吉彪, 等. 氮用量對小黑麥東農(nóng)96026群體生長及飼用品質(zhì)的影響. 中國農(nóng)學通報, 2009, 25(7): 141-144.

[11] Li T. Effect of density and nitrogen nutrition on yield and quality of triticale. Haerbin: Northeast Agricultural University, 2008: 46-50.

李陶. 密度和氮素營養(yǎng)對小黑麥產(chǎn)量與品質(zhì)的影響. 哈爾濱: 東北農(nóng)業(yè)大學, 2008: 46-50.

[12] Zhu X K, Sun J Y, Guo W S,etal. Characteristics of forage yield and quality of different triticale varieties. Barley and Cereal Sciences, 2010, (3): 1-7.

朱新開, 孫建勇, 郭文善, 等. 不同類型小黑麥飼草產(chǎn)量和品質(zhì)特性研究. 大麥與谷類科學, 2010, (3): 1-7.

[13] Yang X H. Effect of yield and quality with N regulation on triticale of different harvest time. Haerbin: Northeast Agricultural University, 2011: 39-42.

楊小輝. 氮素調(diào)控對不同收獲期小黑麥產(chǎn)量及品質(zhì)的影響. 哈爾濱: 東北農(nóng)業(yè)大學, 2011: 39-42.

[14] Zhao Y J, Tian X H, Du W H. Studies on the optimal cutting period of forage triticale in Dingxi area. Pratacultural Science, 2015, 32(7): 1143-1149.

趙雅姣, 田新會, 杜文華. 飼草型小黑麥在定西地區(qū)的最佳刈割期. 草業(yè)科學, 2015, 32(7): 1143-1149.

[15] Zhang B C, Yuan H L, Liu C Y. Analysis and evaluation of nutritional quality for radish variety resources in Anhui. Journal of Crop Variety Resources, 1999, 34(2): 41-42.

張部昌, 袁華玲, 劉才宇. 安徽蘿卜品種資源營養(yǎng)品質(zhì)分析與評價. 作物品種資源, 1999, 34(2): 41-42.

[16] Hou X J, Jiang B. Application of fuzzy mathematics and CI_GSMCM on quality evaluation for tomatoes. Journal of Anhui Agricultural Sciences, 2013, 41(9): 4068-4070.

侯曉靜, 姜波. 模糊數(shù)學和CI_GSMCM在加工番茄品質(zhì)綜合評價中的應用. 安徽農(nóng)業(yè)科學, 2013, 41(9): 4068-4070.

[17] Yang L. The nutrient analysis and quality assessment of characterization and utilization ofHemerocallis. Changchun: Jilin Agricultural University, 2014: 8-10.

楊利. 萱草屬植物營養(yǎng)成分分析及品質(zhì)評價. 長春: 吉林農(nóng)業(yè)大學, 2014: 8-10.

[18] Yue S J, Wang X M, Zhang B. Application of fuzzy mathematics to taste quality evaluation of super sweet corn. Journal of Maize Science, 2003, 11(2): 96-98.

樂素菊, 王曉明, 張璧. 模糊數(shù)學在超甜玉米食用品質(zhì)綜合評價上的應用. 玉米科學, 2003, 11(2): 96-98.

[19] Chen P, Jiang T, Yu C M,etal. Comprehensive evaluation on nutritive quality of 33 ramie varieties by subordinate function. Hubei Agricultural Sciences, 2015, 54(10): 2435-2439.

陳平, 姜濤, 喻春明, 等. 應用隸屬函數(shù)法評價33個苧麻資源的營養(yǎng)品質(zhì). 湖北農(nóng)業(yè)科學, 2015, 54(10): 2435-2439.

[20] Yan W K, Sheng Q L, Hu Y G,etal. GGE bipot-an ideal tool for studying genotype by environment interaction of regional yield trial data. Acta Agronomica Sinica, 2001, 27(1): 21-28.

嚴威凱, 盛慶來, 胡躍高, 等. GGE疊圖法-分析品種×環(huán)境互作模式的理想方法. 作物學報, 2001, 27(1): 21-28.

[21] Yan W K. Optimal use of biplots in analysis of multi-location variety test data. Acta Agronomica Sinica, 2010, 36(11): 1805-1819.

嚴威凱. 雙標圖分析在農(nóng)作物品種多點試驗中的應用. 作物學報, 2010, 36(11): 1805-1819.

[22] Yan W K, Hunt L A, Sheng Q L,etal. Cultivar evaluation and mega-environment investigation based on the GGE biplot. Crop Science, 2000, 40(3): 597-605.

[23] Yan W K, Hunt L A. Biplot analysis of diallel data. Crop Science, 2002, 42(1): 21-30.

[24] Chai S X, Chang L, Yang R J,etal. Analysis on genotype×environment interaction and yield stability of triticale. Journal of Nuclear Agricultural Sciences, 2011, 25(1): 155-161.

柴守璽, 常磊, 楊蕊菊, 等. 小黑麥基因型與環(huán)境互作效應及產(chǎn)量穩(wěn)定性分析. 核農(nóng)學報, 2011, 25(1): 155-161.

[25] Yang S. Feed and feed quality detection technology. Beijing: China Agricultural University Press, 1998: 330-338.

楊勝. 飼料分析及飼料質(zhì)量監(jiān)測技術. 北京: 中國農(nóng)業(yè)大學出版社, 1998: 330-338.

[26] Li D M, Tian X H, Du W H. Studies on seed yield and yield components of triticale lines. Acta Agrestia Sinica, 2016, 24(1): 241-244.

李冬梅, 田新會, 杜文華. 5個小黑麥新品系的種子產(chǎn)量及產(chǎn)量構(gòu)成因素分析. 草地學報, 2016, 24(1): 241-244.

[27] Xu X H, He M Z. Experimental design and application of design-expert and SPSS. Beijing: Science Press, 2010: 9-19.

徐向宏, 何明珠. 試驗設計與Design-Expert、SPSS應用. 北京: 科學出版社, 2010: 9-19.

[28] Chang L, Chai S X. Application of GGE biplot in spring wheat yield stability analysis in rainfed areas of China. Chinese Journal of Eco-Agriculture, 2010, 18(5): 988-994.

常磊, 柴守璽. GGE雙標圖在我國旱地春小麥穩(wěn)產(chǎn)性分析中的應用. 中國生態(tài)農(nóng)業(yè)學報, 2010, 18(5): 988-994.

[29] He J F, Zhao M L. Forage characteristics and application perspective of Triticale in grassland ecology. Chinese Journal of Grassland, 2012, 34(1): 1101-1107.

何江峰, 趙萌莉. 小黑麥的飼用特性及其在草地生態(tài)中的應用前景. 中國草地學報, 2012, 34(1): 1101-1107.

[30] Zhang R. Survey and identification of the alpine grassland major fungal diseases in southern Gansu province. Lanzhou: Gansu Agricultural University, 2009: 40-43.

張蓉. 甘南高寒草地植物主要真菌病害調(diào)查與鑒定. 蘭州: 甘肅農(nóng)業(yè)大學, 2009: 40-43.

[31] Song Q, Tian X H, Du W H. Studies on production performance of new triticale lines in alpine pastoral areas of Gansu. Pratacultural Science, 2016, 33(7): 1367-1374.

宋謙, 田新會, 杜文華. 甘肅省高寒牧區(qū)小黑麥新品系的生產(chǎn)性能. 草業(yè)科學, 2016, 33(7): 1367-1374.

[32] Song Q. Evaluation on cold resistance and forage production of new triticale lines in alpine pastoral area of Gansu province. Lanzhou: Gansu Agricultural University, 2016: 39-40.

宋謙. 飼草型小黑麥新品系的抗寒性評價及在甘肅省高寒牧區(qū)的生產(chǎn)性能研究. 蘭州: 甘肅農(nóng)業(yè)大學, 2016: 39-40.

[33] Zhao D. ISSR analysis and studies on the productivity of triticale new lines. Lanzhou: Gansu Agricultural University, 2014: 68-69.

趙丹. 小黑麥新品系遺傳多樣性的ISSR分析及生產(chǎn)性能研究. 蘭州: 甘肅農(nóng)業(yè)大學, 2014: 68-69.

[34] Tian J S, Tian X H, Ka Z J,etal. Climate dynamic in Maqu county from 1981 to 2010. Grassland and Turf, 2015, 35(1): 84-88.

田久勝, 田新會, 卡召加, 等. 甘肅省瑪曲縣1981-2010年氣候變化動態(tài). 草原與草坪, 2015, 35(1): 84-88.