宋 慧 劉金榮* 王素英 閆宏山 王 濤 邢曉寧 邢 璐 付 楠
(1.安陽市農(nóng)業(yè)科學(xué)院 谷子研究所,河南 安陽 455000; 2.安陽工學(xué)院 生物與食品工程學(xué)院,河南 安陽 455000)
谷子(Setariaitalica(L.) Beauv.)為禾本科狗尾草屬作物,起源于中國,在中國農(nóng)耕文明中占舉足輕重的地位[1-2]。在深化農(nóng)業(yè)供給側(cè)結(jié)構(gòu)性改革和實現(xiàn)鄉(xiāng)村振興戰(zhàn)略政策的引導(dǎo)下,培育適應(yīng)性廣、高產(chǎn)、穩(wěn)產(chǎn)和適合機械化作業(yè)等突破性谷子新品種,是當(dāng)前谷子產(chǎn)業(yè)發(fā)展的方向之一。據(jù)國家谷子糜子產(chǎn)業(yè)技術(shù)體系產(chǎn)業(yè)經(jīng)濟崗位調(diào)查結(jié)果[3],2013年‘豫谷18’全國種植面積5 333 hm2,2014年約 1.3萬 hm2,2015年新疆華大谷業(yè)科技有限公司僅在新疆推廣約 1.3萬 hm2,2013—2017年連續(xù)被農(nóng)業(yè)部列為華北夏谷區(qū)谷子生產(chǎn)上主推品種,推廣應(yīng)用前景良好。以華北夏谷區(qū)組試驗為例對該品種進行豐產(chǎn)穩(wěn)產(chǎn)和適應(yīng)性科學(xué)評價,找到基因型與環(huán)境互作效應(yīng)可重復(fù),提高谷子新品種的豐產(chǎn)穩(wěn)產(chǎn)性和適應(yīng)性育種的選擇效率,對持續(xù)提升谷子產(chǎn)業(yè)綜合競爭力和可持續(xù)發(fā)展能力具有重要意義。
GGE雙標(biāo)圖(GGE,genotype main effects and genotype-environment interaction effects)通過分析基因型主效和基因型與環(huán)境互作效應(yīng),以模型的方式科學(xué)描述和量化表達[4-5]。目前國內(nèi)外已被廣泛應(yīng)用于棉花[6-9]、小麥[10]、玉米[11-12]、水稻[13]、燕麥[14]、小豆[15]、糜子[16]、花生[17]和蕎麥[18]等作物的品種豐產(chǎn)性與穩(wěn)產(chǎn)性分析、品種適宜種植區(qū)域劃分[19],種質(zhì)資源評價[20]和試驗環(huán)境評價[21-24]等研究。以谷子為例,用GGE雙標(biāo)圖評價品種的豐產(chǎn)穩(wěn)產(chǎn)性和適應(yīng)性,并為品種推薦適宜種植區(qū)域,在谷子新品種推廣和應(yīng)用中報道尚少。鑒于此,本研究應(yīng)用2010—2011年國家谷子品種區(qū)域試驗華北夏谷區(qū)組‘豫谷18’所在組別參試品種的谷子小區(qū)產(chǎn)量數(shù)據(jù),通過GGE雙標(biāo)圖的“豐產(chǎn)性與穩(wěn)定性”功能圖、“理想品種”功能圖和“適宜種植區(qū)域劃分”功能圖對‘豫谷18’等品種的豐產(chǎn)性與穩(wěn)產(chǎn)性、適宜種植區(qū)域劃分進行分析,并采用“成對比較”功能圖分析‘豫谷18’與對照品種‘冀谷19’在目標(biāo)區(qū)域的適應(yīng)性表現(xiàn),旨在為‘豫谷18’的合理利用提供理論依據(jù),以期為其他作物品種的高產(chǎn)穩(wěn)產(chǎn)性與適應(yīng)性評價提供參考方法。
研究數(shù)據(jù)來源于2010—2011年華北夏谷區(qū)‘豫谷18’所在組別谷子新品種區(qū)域試驗的產(chǎn)量數(shù)據(jù)(表1)。2010年參試品種10個,2011年參試品種15個,‘冀谷19’為對照品種,2010年試驗在河北省石家莊市高新區(qū)(簡稱為SJZ,下同)、河北省唐山市灤縣(LX)、河北省邯鄲市武安市(WA)、河北省保定市蓮池區(qū)(BD)、河北省滄州市運河區(qū)(CZ)、河北省衡水市桃城區(qū)(HS),河南省安陽市北關(guān)區(qū)(AY)、河南省洛陽市洛龍區(qū)(LY)、河南省鄭州市金水區(qū)(ZZ),山東省濟南市歷城區(qū)(JN)和北京市昌平區(qū)(BJ)共11個試點。2011年新增河北省邯鄲市邯山區(qū)(HD)、山東省臨沂市莒縣(LX),遼寧省錦州市太和區(qū)(JZ),山東省德州市德城區(qū)(DZ)。河北省武安市不再作為參試點,總計14個試點。不同試點經(jīng)緯度、海拔等地理信息,見表2。
試驗采用隨機區(qū)組排列,3次重復(fù),6~8行,行距0.4 m,小區(qū)面積20 m2,留苗密度60萬株/hm2,收獲時按實收小區(qū)產(chǎn)量。
采用GGE-biplot軟件(基因和基因×環(huán)境雙標(biāo)圖)參照《農(nóng)作物品種試驗數(shù)據(jù)管理分析》[25]對參試品種進行豐產(chǎn)性和穩(wěn)定性、理想品種、適宜種植區(qū)域劃分和成對比較評價。
表3可知, 2010—2011年國家谷子區(qū)試華北夏谷區(qū)組品種產(chǎn)量多環(huán)境方差分析均達到差異極顯著水平。其中2010年基因型主效應(yīng)、環(huán)境主效應(yīng)和基因×環(huán)境主效應(yīng)的變異分別占處理總變異的11.6%、43.56%和27.3%,2011年分別占總變異的10.8%、60.9%和17.7%??梢姯h(huán)境主效應(yīng)是產(chǎn)量變異來源,基因型×環(huán)境互作效應(yīng)2年試驗中都大于基因型效應(yīng),采用GGE模型有效剔除環(huán)境效應(yīng)的影響,科學(xué)評價基因型、基因型×環(huán)境互作效應(yīng),推薦各參試品種的適宜種植區(qū)域。
由圖1(a)可知,2010年國家谷子品種區(qū)試華北夏谷區(qū)組中第1主成分(PC1)效應(yīng)37.3%,第2主成分(PC2)的效應(yīng)為21.6%,基因型和環(huán)境互作效應(yīng)58.9%。各品種的豐產(chǎn)性表現(xiàn)為‘豫谷18’(YG18)>‘206058’>‘豫谷19’(YG19)>‘200131’>‘保213’(B213)>‘濟0515’(J0515) >‘滄372’(C372) >‘衡2001’(H2001) >‘冀谷19’(JG19) >‘鄭06-6’(Z06-6) >‘A2×測351’(A2C),‘豫谷18’最接近AEA軸垂足的正方向,與其他品種距離較遠,表明‘豫谷18’的豐產(chǎn)性優(yōu)于其他參試品種,與品種間差異多重比較結(jié)果一致(表1)。參試品種穩(wěn)產(chǎn)性表現(xiàn)為(圖1(a)),‘豫谷19’和‘豫谷18’在AEC軸的投影距離短,穩(wěn)產(chǎn)指數(shù)低,穩(wěn)產(chǎn)性最好,‘濟0515’、‘衡2001’、‘冀谷19’、‘鄭06-6’穩(wěn)產(chǎn)性較好,‘206058’、‘200131’、‘滄372’穩(wěn)產(chǎn)性較差,‘保213’穩(wěn)產(chǎn)性最差。由表1和圖1(b)可知,2011年參試品種的豐產(chǎn)性表現(xiàn)為:‘豫谷18’ >‘京安7505’ >‘濟0506’>‘豫谷19’>‘200475’>‘M1508’>‘K1174’>‘保213’>‘鄭07-1’>‘冀谷19’>‘石38’>‘A2’ב夏父1’>‘滄372’>‘A2×4170’>‘滄369’>‘保200302’,其中‘豫谷18’的豐產(chǎn)性優(yōu)于其他參試品種。各參試品種的穩(wěn)產(chǎn)性表現(xiàn)為,有‘豫谷18’和‘京安7505’穩(wěn)產(chǎn)性最好,‘豫谷19’、‘M1508’、‘保213’,以‘A2×4170’和‘200475’穩(wěn)產(chǎn)性最差。綜合分析表明,‘豫谷18’在2年區(qū)域試驗中豐產(chǎn)性和穩(wěn)產(chǎn)性均表現(xiàn)良好。
表1 2010—2011年國家谷子品種區(qū)域試驗華北夏谷區(qū)組參試品種信息表Table 1 Discription of foxtail millet cultivars in the national foxtail millet region trials in the summer sowing region of North China in 2010-2011
注:不同小寫和大寫字母分別表示在5%和1%水平上差異顯著。
Note:Yield data in each year followed by different lowercases and capital letters are significantly different at 5% and 1% probability levels, respectively.
表2 2010—2011 年國家谷子區(qū)域試驗華北夏谷區(qū)組環(huán)境的地理因子Table 2 Geographical factors and planting regions of test locations of national foxtail millet regional trials in the summer sowing region of North China
表3 2010—2011年國家谷子區(qū)域試驗華北夏谷區(qū)谷子品種產(chǎn)量的聯(lián)合方差分析Table 3 Combined analysis of variance for foxtail millet yield from the national foxtail regional trials in the summer sowing region of North China in 2010-2011
注:**表示差異極顯著(P<0.01).
Note:**indicates significant difference at 0.01 level.
PC1,第一主成分;PC2,第二主成分。帶“*”表示品種,“+”表示參試點。下同?!皁”為各試點環(huán)境坐標(biāo)的平均值,代表目標(biāo)環(huán)境;通過原點和平均環(huán)境的直線是平均環(huán)境軸,AEA軸的箭頭方向為正方向,豐產(chǎn)性由雙標(biāo)圖中參試品種在AEA軸上相對位置對應(yīng),越靠近正方向說明豐產(chǎn)性越好;穩(wěn)定性是由參試品種向AEA軸做垂線,距離越遠穩(wěn)定性越差,反之則穩(wěn)定性好。
C372,滄372;B213, 保213;v206058,206058;H2001,衡2001;v200131,200131;YG19,豫谷19;YG18,豫谷18;J0515,濟0515;Z06-6,鄭06-6;A2C,A2×測351變;JG19,冀谷19;M1,M1508;K1,K1174;z07,鄭07-1;v20,200475;B20,保200302;J05,濟0506;C369,滄369;S38,石3839;JA,京安7505;A24,A2×4170;A2X,A2×夏父1。下同。
PC1, first principal component; PC2, second principal component. The same below. The plus sign “+” shows the test location. “*” is cultivar. The small circle in the
Figure is the average value of the coordinates of the position of the test location, representing the target environment. The straight line through the origin and the average environment is the average environment axis, and the arrow direction of the AEA axis is in the positive direction. The relative position of the tested varieties on the AEA axis corresponds to the high yield and the closer they are to the positive direction, the better the high yield. The stability is represented by the vertical line from tested varieties to AEA axis. The farther the distance the worse the stability, and vice versa, the better the stability. C372, Cang 372; B213, Bao 213; v206058, 206058; H2001, Heng 2001; v200131, 200131; YG19, Yugu 19; YG18, Yug 18; J0515, Ji 0515; Z06-6, Zheng 06-6; A2C, A2×Ce 351 Bian; JG19, Jigu 19; M1, M1508; K1, K1174; z07, Zheng 07-1; v20, 200475; B20, Bao 200302; J05, Ji 0506; C369, Cang 369; S38, Shi 3839; JA, Jingan 7505; A24, A2×4170; A2X, A2×Xiafu 1. The same below.
圖1 2010(a)和2011年(b)華北夏谷區(qū)組谷子產(chǎn)量GGE雙標(biāo)圖的豐產(chǎn)性與穩(wěn)產(chǎn)性功能圖
Fig.1 Mean and stability view of the GGE biplot for millet yield data from the millet regional trials in the summer sowing region in north of China in 2010 (a) and 2011 (b)
圖2(a)和表1可知,2010年各參試品種中‘豫谷18’距離理想品種最近,理想指數(shù)最小,說明該品種的豐產(chǎn)性和穩(wěn)產(chǎn)性綜合表現(xiàn)最好;‘豫谷19’和‘206058’,‘200131’、‘滄372’、‘濟0515’、‘保213’較為理想;‘衡2001’、‘鄭06-6’表現(xiàn)較差;‘A2×測351’綜合表現(xiàn)最差。2011年理想指數(shù)最小的品種是‘豫谷18’和‘京安7505’,其次是‘豫谷19’、‘保283’、‘M1508’、‘K1174’、‘200475’,‘A2×4170’和‘保200302’表現(xiàn)差(表1和圖2(b))?!ス?8’在2年區(qū)域試驗中豐產(chǎn)性和穩(wěn)產(chǎn)性綜合表現(xiàn)顯著優(yōu)于其他參試品種和對照品種。
“o”是理想品種的圖標(biāo),以理想品種為圓心畫出若干個同心圓,越接近圓心的品系理想程度越好。
“o” at the center of concentric circle is the mark of the ideal cultivar. Draw several concentric circles with the ideal cultivar as the center of the circle. The closer to the center of the circle is the better the cultivar.
圖2 2010(a)和2011年(b)國家谷子區(qū)域試驗華北夏谷區(qū)組谷子產(chǎn)量GGE雙標(biāo)圖的“理想品種”功能圖
Fig.2 Ideal cultivar view of the GGE biplot for millet yield data from the millet regional trials in the summer sowing region in north of China in 2010 (a) and 2011 (b)
由圖3(a)可知,2010年‘豫谷18’在河北省石家莊市高新區(qū)、衡水市桃城區(qū)、滄州市運河區(qū),河南省安陽市北關(guān)區(qū)、洛陽市洛龍區(qū)、鄭州市金水區(qū)和山東省臨沂市莒縣等7個試點綜合表現(xiàn)良好,‘206058’在山東省濟南市歷城區(qū)、河北省保定市蓮池區(qū)、唐山市灤縣表現(xiàn)最好。由圖3(b)可知,2011年‘豫谷18’,‘京安7505’、‘M1508’在河南省安陽市北關(guān)區(qū)、洛陽市洛龍區(qū),山東省德州市德城區(qū)、濟南市歷城區(qū),臨沂市莒縣,河北省衡水桃城區(qū)、石家莊市高新區(qū)、唐山市灤縣、滄州市運河區(qū)、邯鄲市邯城區(qū),北京市昌平區(qū)等11個試點表現(xiàn)好,適宜種植,‘206058’僅在鄭州市1個試點表現(xiàn)好,‘A2×4170’在錦州試點表現(xiàn)最好。綜述之,‘豫谷18’能夠適應(yīng)華北夏谷區(qū)組大部分試點,為該區(qū)組適應(yīng)性最廣的品種。
圖4(a)可知,2010年‘豫谷18’和‘冀谷19’位于等值線的兩側(cè),試驗環(huán)境主要位于等值線右側(cè)與‘豫谷18’同側(cè),‘豫谷18’與等值線之間的距離大于‘冀谷19’與等值線的距離,說明‘豫谷18’在不同試點中的豐產(chǎn)性綜合表現(xiàn)優(yōu)于‘冀谷19’。2011年有相似的表現(xiàn)(圖4(b))?!ス?8’在華北夏谷區(qū)域試點中比對照品種更具有推廣優(yōu)勢。
在作物區(qū)域試驗中產(chǎn)量受基因型主效、環(huán)境主效×基因型與環(huán)境互作效應(yīng)的共同影響,同一個品種在不同的參試環(huán)境中的產(chǎn)量表現(xiàn)的差異會很大,對選育適應(yīng)性廣、高產(chǎn)穩(wěn)產(chǎn)品種的選育帶來困擾?;蛐汀镰h(huán)境互作效應(yīng)對品種變異的影響通常大于基因型主效[20],是影響品種穩(wěn)定性和適應(yīng)性的主要因素。在區(qū)域試驗中科學(xué)評價和利用可重復(fù)的基因型環(huán)境互作效應(yīng),對提高育種效率,選育適應(yīng)性廣高產(chǎn)穩(wěn)產(chǎn)品種起關(guān)鍵作用[26]。目前,對作物品種豐產(chǎn)和穩(wěn)產(chǎn)性進行評價應(yīng)用最多的是AMMI 模型,但是這些分析方法主要通過分析基因型與環(huán)境互作效應(yīng),對基因型進行全面評價,忽視基因型主效、基因與環(huán)境主效綜合評價。而GGE雙標(biāo)圖模型將基因型主效和基因型與環(huán)境互作效應(yīng)相結(jié)合[27],比AMMI模型能更多地解釋基因型和基因×環(huán)境互作效應(yīng),完成對品種的豐產(chǎn)性和穩(wěn)產(chǎn)性同步評價和選擇。本研究采用GGE雙標(biāo)圖的“豐產(chǎn)性和穩(wěn)定性”和“理想品種”功能圖對2010—2011年國家谷子區(qū)域試驗華北夏谷區(qū)組中的‘豫谷18’等參試品種的豐產(chǎn)性和穩(wěn)產(chǎn)性分析表明,在2年區(qū)域試驗中,‘豫谷18’都在AEA軸正方向的最遠處,居參試品種第一位,在AEC軸的投影相對較短,穩(wěn)產(chǎn)性優(yōu)良;與參試品種相比,‘豫谷18’理想指數(shù)最小,豐產(chǎn)性和穩(wěn)產(chǎn)性綜合評價第一,是高產(chǎn)穩(wěn)產(chǎn)的優(yōu)良新品種。本研究科學(xué)評價‘豫谷18’的適應(yīng)廣、豐產(chǎn)和穩(wěn)產(chǎn)性特點,為該品種科學(xué)合理推廣應(yīng)用提供了理論依據(jù),同時也進一步證實GGE雙標(biāo)圖模型在谷子品種高產(chǎn)穩(wěn)產(chǎn)性同步評價方面的可行性。
同一方向上距離原點最遠的品種依次連接形成一個多邊形,由雙標(biāo)圖原點向多邊形各邊作垂線,形成多個不同扇區(qū),分布在同一個扇區(qū)的品種屬于相同環(huán)境組合,每個扇區(qū)中位于多邊形頂部的品種為該環(huán)境組合中最優(yōu)品種。
SJZ,河北省石家莊市;LX,河北省灤縣市;WA,河北省武安市;BD,河北省保定市;CZ,河北省滄州市;HS,河北省衡水市;HD,河北省邯鄲市;LX,河北省莒縣;AY,河南省安陽市;LY,河南省洛陽市;ZZ,河南省鄭州市;JN,山東省濟南市;DZ,山東省德州市;BJ,北京市;JZ,遼寧省錦州市。
The varieties farthest from the origin in the same direction are connected in turn to form a polygon. The vertical line is drawn from the origin of GGE biplot to each side of the polygon to form many different sectors. Varieties distributed in the same sector belong to the same environmental combination. In each sector, the varieties located in polygons are the best environmental combination. SJZ, Shijiazhuang, Hebei; LX, Luanxian, Hebei; WA, Wuan, Hebei; BD, Baoding, Hebei; CZ, Cangzhou, Hebei; HS, Hengshui, Hebei; HD, Handan, Hebei; LX, Juxian, Hebei; AY, Anyang, Henan; LY, luoyang, Henan; ZZ, Zhengzhou, Henan; JN, Jinan, Shandong; DZ, Dezhou, Shandong; BJ, Beijing; JZ, Jinzhou, Liaoning.
圖3 2010(a)和2011年(b)華北夏谷區(qū)組谷子產(chǎn)量GGE雙標(biāo)圖的“適宜種植區(qū)域劃分”功能圖
Fig.3 Which-Won-Where view of the GGE biplot for millet yield data from the milletregional trials in the summer sowing region in north of China in 2010 (a) and 2011 (b)
過雙標(biāo)圖原點對直線做垂線也稱等值線,品種在等值線上的任意環(huán)境中數(shù)值相等,距離等值線越遠在目標(biāo)環(huán)境表現(xiàn)越好。
Connect the icons of two comparative varieties with a strait directly and draw a vertical line from the origin of GGE biplot to this strait. This vertical line is called an isoline and the values of varieties are equal in any environment on the isoline. The farther the cultivar is from the isoline, the better the performance of this cultivar in the target environment.
圖4 2010(a)和2011年(b)國家谷子區(qū)域試驗華北夏谷區(qū)組谷子產(chǎn)量GGE雙標(biāo)圖的“成對比較”功能圖
Fig.4 Pairwise comparison view of the GGE biplot for millet yield data from the millet regional trials in the summer sowing region in north of China in 2010 (a) and 2011 (b)
適應(yīng)性廣、豐產(chǎn)和穩(wěn)產(chǎn)是評價農(nóng)作物品種的最重要的經(jīng)濟指標(biāo)和產(chǎn)量育種目標(biāo)[30]。20世紀(jì)80年代‘豫谷1號’的選育打破了谷子適應(yīng)性窄、地域性強的傳統(tǒng)理論??稍谠ァ⒓?、魯、遼和陜等省大面積推廣種植,獲國家技術(shù)發(fā)明二等獎,到目前生產(chǎn)上主推品種‘豫谷18’,都為適應(yīng)性廣谷子新品種的選育積累了優(yōu)秀的育種材料。因此,筆者總結(jié)團隊多年的育種經(jīng)驗[28-29],分別針對谷子的穩(wěn)產(chǎn)性、豐產(chǎn)性和廣適應(yīng)性狀的育種策略進行探討。對谷子穩(wěn)產(chǎn)性狀的選擇主要依據(jù)品種抗逆性、抗病性和適應(yīng)性,低世代關(guān)注單株抗倒性、垂直抗病性的選擇,淘汰易倒伏、感病的株系,高世代關(guān)注群體抗倒性、水平抗病性,對重點株系進行接種,設(shè)置抗逆抗病鑒定圃,確保品種抗逆性的遺傳。對谷子豐產(chǎn)性狀選擇通過提高單位面積有效穗數(shù)和單穗粒重的目標(biāo),低世代通過測定單株的穗長、穗粗、穗碼數(shù)、碼粒數(shù)等產(chǎn)量性狀達到提高單穗粒重目的,一般谷穗中等偏長、粗穗、穗碼數(shù)適中、碼粒數(shù)多的類型為主。高世代以評測群體為主,設(shè)置高密度選種圃,選擇成穗率高、結(jié)實性好、抗倒性強的株系。這與小麥高產(chǎn)高效品種選育有相似之處[30]。參照程順和等[31]提出高代品系進行不同生態(tài)區(qū)多點試驗提高適應(yīng)性育種效率的觀點,在谷子的選育過程中首先以光溫不敏感的材料/品種作親本配制雜交組合,在安陽市農(nóng)業(yè)科學(xué)院的安陽市和海南省三亞市試驗基地進行跨生態(tài)區(qū)交替選擇,由于海南省三亞市崖城鎮(zhèn)特殊氣候條件下,低世代材料不輕易淘汰,以防漏選,高代材料布置多點試驗鑒定,結(jié)合綜合表現(xiàn),重點選擇高產(chǎn)穩(wěn)產(chǎn)、抗病、一致性好的株系,同時考慮成穗率高、邊際效應(yīng)弱、耐密植株系,育成品種自身調(diào)節(jié)能力強,利于穩(wěn)產(chǎn)和廣適應(yīng)。
本研究采用GGE雙標(biāo)圖,為‘豫谷18’劃定優(yōu)勢種植區(qū)域,在華北夏谷區(qū)大部分試點適合該品種的推廣,同時利用 GGE 雙標(biāo)圖的“成對比較”功能圖將‘豫谷18’與對照品種‘冀谷19’在華北夏谷區(qū)組的產(chǎn)量表現(xiàn)進行相互比較。結(jié)果表明,‘豫谷18’在2年試驗中表現(xiàn)出對華北夏谷區(qū)域最廣泛的適應(yīng)性,在與對照品種的成對比較中也體現(xiàn)出絕對的區(qū)域適應(yīng)性優(yōu)勢。明確‘豫谷18’是華北夏谷區(qū)廣適性的品種,結(jié)合其高產(chǎn)和穩(wěn)產(chǎn)的特征,說明‘豫谷18’是針對華北夏谷區(qū)廣適性育種的成功例證,同時也說明 GGE 雙標(biāo)圖對品種的適宜種植區(qū)域劃分的有效性和實用性。
中國農(nóng)業(yè)大學(xué)學(xué)報2020年1期