張偉民 高樹(shù)廣 胡敏杰 張春花 王瑞霞 徐東陽(yáng) 李偉峰 徐博涵 楊光宇

（1.河南省西華縣農(nóng)業(yè)科學(xué)研究所 河南 西華 466600；2.周口市農(nóng)業(yè)科學(xué)院 河南 周口 466001）

芝麻是世界上最古老的油料作物之一，也是我國(guó)的特色優(yōu)質(zhì)農(nóng)產(chǎn)品[1]，被稱為“萬(wàn)油之王”[2]。芝麻不僅富含對(duì)人體有益的不飽和脂肪酸、氨基酸和礦物質(zhì)，還含有維生素E、芝麻醇、芝麻素、芝麻蛋白、植物甾醇等功能物質(zhì)[5-8]。我國(guó)芝麻的產(chǎn)量和品質(zhì)均位居世界前列，其中河南省是中國(guó)的芝麻主產(chǎn)區(qū)，芝麻種植面積和芝麻產(chǎn)量均占我國(guó)芝麻總種植面積和芝麻總產(chǎn)量的30%[4]。近年來(lái)，隨著經(jīng)濟(jì)的發(fā)展和人民生活水平的不斷提高，人們對(duì)營(yíng)養(yǎng)保健食品的需求量越來(lái)越高，芝麻作為一種重要的藥用和食用作物，在國(guó)內(nèi)需求中有著穩(wěn)定的增長(zhǎng)[9-10]。據(jù)相關(guān)資料顯示，2018年，中國(guó)芝麻需求量達(dá)到121.8萬(wàn)t，其中進(jìn)口82.82萬(wàn)t[11]，超過(guò)我國(guó)芝麻生產(chǎn)總量的160%，表明我國(guó)芝麻產(chǎn)量遠(yuǎn)遠(yuǎn)不能滿足市場(chǎng)的需求。因此，為緩解國(guó)內(nèi)芝麻生產(chǎn)與需求的矛盾，保障國(guó)家糧油安全，選育高產(chǎn)、穩(wěn)產(chǎn)、優(yōu)質(zhì)、多抗的芝麻新品種，提高芝麻單產(chǎn)和總產(chǎn)量，實(shí)現(xiàn)農(nóng)業(yè)增產(chǎn)增效已成為當(dāng)前芝麻育種研究的首要任務(wù)。

周芝11號(hào)是周口市農(nóng)業(yè)科學(xué)院以鄭芝13號(hào)為母本、漯芝16號(hào)為父本，通過(guò)有性雜交選育的芝麻新品種，于2021年通過(guò)河南省種子管理站鑒定（鑒定編號(hào)：豫品鑒芝2021007）。該品種屬單稈型，三花四棱，莖稈綠色，茸毛中等，花色乳白，蒴果致密，成熟時(shí)蒴果輕裂，籽粒純白。為明確周芝11號(hào)的推廣應(yīng)用價(jià)值，本研究利用高穩(wěn)系數(shù)和回歸系數(shù)法分析了2018年、2019年黃淮流域芝麻新品種區(qū)域試驗(yàn)和2020年展示試驗(yàn)數(shù)據(jù)資料，進(jìn)一步明確了周芝11號(hào)的高產(chǎn)、穩(wěn)產(chǎn)、適應(yīng)性和品質(zhì)特性，為周芝11號(hào)的大規(guī)模推廣應(yīng)用提供了堅(jiān)實(shí)的理論依據(jù)。

1 材料與方法

1.1 試驗(yàn)材料和試驗(yàn)地點(diǎn)

1.1.1 試驗(yàn)材料2018年黃淮流域芝麻新品種區(qū)域試驗(yàn)參試品系19個(gè)，2019年黃淮流域芝麻新品種區(qū)域試驗(yàn)參試品系22個(gè)，2020年全國(guó)芝麻新品種展示試驗(yàn)（黃淮產(chǎn)區(qū)）參試品系8個(gè)，對(duì)照品種均為豫芝4號(hào)。1.1.2試驗(yàn)地點(diǎn)2018年、2019年區(qū)試地點(diǎn)6個(gè)，

分別是信陽(yáng)市農(nóng)業(yè)科學(xué)院、南陽(yáng)市農(nóng)業(yè)科學(xué)院、阜陽(yáng)市農(nóng)業(yè)科學(xué)院、駐馬店市農(nóng)業(yè)科學(xué)院、宿州市農(nóng)業(yè)科學(xué)院、河南省農(nóng)業(yè)科學(xué)院原陽(yáng)基地；2020年展示試驗(yàn)地點(diǎn)3個(gè)，分別是河南省農(nóng)業(yè)科學(xué)院原陽(yáng)基地、駐馬店市農(nóng)業(yè)科學(xué)院、安徽省農(nóng)業(yè)科學(xué)院。

1.2 試驗(yàn)方法

1.2.1 試驗(yàn)設(shè)計(jì)2018年、2019年區(qū)域試驗(yàn)均分為2組，采取隨機(jī)區(qū)組設(shè)計(jì)，10行區(qū)，行長(zhǎng)5 m，行距0.4 m，小區(qū)面積20 m2，條播，3次重復(fù)，重復(fù)間設(shè)走道，試驗(yàn)地四周設(shè)保護(hù)行。2020年展示試驗(yàn)，采取隨機(jī)區(qū)組設(shè)計(jì)，20行區(qū)，行長(zhǎng)5 m，行距0.4 m，小區(qū)面積40 m2，2次重復(fù)，重復(fù)間設(shè)走道，四周設(shè)保護(hù)行。

1.2.2 試驗(yàn)管理 各試驗(yàn)點(diǎn)均按照統(tǒng)一試驗(yàn)方案進(jìn)行田間管理、觀察記載和考種，每一項(xiàng)田間管理和測(cè)定項(xiàng)目均在同一天內(nèi)完成。全生育期只對(duì)蟲(chóng)害、雜草采取防控措施，不對(duì)病害進(jìn)行人工防控。

1.2.3 基因組信息測(cè)定 為評(píng)價(jià)參試品種的特異性和穩(wěn)定性，采用Illumina高通量測(cè)序技術(shù)對(duì)參試樣品進(jìn)行重測(cè)序。

1.3 數(shù)據(jù)分析

1.3.1 豐產(chǎn)性分析 采用Excel和SPSS 17.0進(jìn)行數(shù)據(jù)處理和統(tǒng)計(jì)分析[15]。采用一般通用公式計(jì)算產(chǎn)量平均值（X）和標(biāo)準(zhǔn)差，采用參試品種與對(duì)照品種相比較的方法進(jìn)行豐產(chǎn)性分析[16]，采用Duncan's新復(fù)極差法比較參試品種之間產(chǎn)量差異顯著性。

1.3.2 穩(wěn)產(chǎn)性分析 以標(biāo)準(zhǔn)差與產(chǎn)量的比值計(jì)算變異系數(shù)，用變異系數(shù)（CV）衡量品種穩(wěn)產(chǎn)性[17]；采用溫振民1994年提出的簡(jiǎn)化公式（公式1）計(jì)算高穩(wěn)系數(shù)（HSC），以高穩(wěn)系數(shù)來(lái)綜合分析高產(chǎn)性和穩(wěn)產(chǎn)性[18]。

式中，HSCi表示第i個(gè)參試品種的高穩(wěn)系數(shù)，Xi表示第i個(gè)品種的平均產(chǎn)量，Si表示第i個(gè)品種的標(biāo)準(zhǔn)差，XCK表示對(duì)照品種的平均產(chǎn)量。

1.3.3 適應(yīng)性分析 采用Eberhart-Russell[17]提出的公式計(jì)算回歸系數(shù)（公式2），用回歸系數(shù)法[19]來(lái)分析品種的適應(yīng)性。

式中，bi表示第i個(gè)品種的回歸系數(shù)；Yij表示第i個(gè)品種在第j個(gè)區(qū)試點(diǎn)中的平均產(chǎn)量；Ij表示第j個(gè)區(qū)試點(diǎn)的環(huán)境指數(shù)，即第j個(gè)區(qū)試點(diǎn)的平均產(chǎn)量與整個(gè)區(qū)試總平均產(chǎn)量的差值。

2 結(jié)果與分析

2.1 基因組信息測(cè)定與分析

周芝11號(hào)與豫芝11號(hào)基因組比對(duì)，mapping率為96.19%、特異SNP 446 534個(gè)、雜合SNP 331 663個(gè)、SNP雜合率為0.10%，比對(duì)效率高，數(shù)據(jù)質(zhì)量高，符合同物種基因組信息比對(duì)要求。周芝11號(hào)與100個(gè)參試樣本的遺傳距離變幅在0.126 3～0.515 6，遺傳相似性范圍為48.84%～83.34%，相似性低于90%，符合參試品種特異性要求。周芝11號(hào)基因組重測(cè)序reads與 篩 選 元 件CaMV35S、Bar、Bt和NOS進(jìn) 行 序 列 比對(duì)，表明周芝11號(hào)不含有轉(zhuǎn)基因成分。

2.2 產(chǎn)量構(gòu)成要素分析

芝麻高產(chǎn)育種要以株蒴數(shù)和蒴粒數(shù)多、千粒重高、果軸較長(zhǎng)、始蒴位低且株高適中為主要選擇目標(biāo)[20]。根據(jù)2018-2019年黃淮區(qū)芝麻新品種區(qū)域試驗(yàn)統(tǒng)計(jì)結(jié)果（表1），周芝11號(hào)平均株高178.3 cm，始蒴高度69.8 cm，黃梢尖4.0 cm，果軸104.5 cm；單株蒴數(shù)88.7個(gè)，蒴粒數(shù)61粒，千粒重3.018 g。周芝11號(hào)主要產(chǎn)量要素協(xié)調(diào)、穩(wěn)定，年際間變化小，符合高產(chǎn)芝麻品種對(duì)產(chǎn)量構(gòu)成因素的要求。

表1 2018-2019年區(qū)域試驗(yàn)中周芝11號(hào)主要產(chǎn)量構(gòu)成要素分析

2.3 豐產(chǎn)性分析

由表2可知，2018年黃淮流域芝麻新品種區(qū)域試驗(yàn)，周芝11號(hào)平均單產(chǎn)1 325.30 kg/hm2，比豫芝4號(hào)（對(duì)照）增產(chǎn)12.08%，增產(chǎn)達(dá)極顯著水平，4點(diǎn)次增產(chǎn)，增產(chǎn)點(diǎn)次占80%，居第二組參試品種第1位。2019年黃淮區(qū)芝麻新品種區(qū)域試驗(yàn)，周芝11號(hào)在6個(gè)試驗(yàn)點(diǎn)平均單產(chǎn)1 506.33 kg/hm2，比豫芝4號(hào)（對(duì)照）增產(chǎn)10.86%，6點(diǎn)次均增產(chǎn)，增產(chǎn)點(diǎn)次占100%，居第一組參試品種第2位。2年11點(diǎn)次平均單產(chǎn)1 415.82 kg/hm2，平均比豫芝4號(hào)增產(chǎn)11.43%，2年11個(gè)點(diǎn)有10個(gè)點(diǎn)比對(duì)照增產(chǎn)，增產(chǎn)點(diǎn)次占90.91%。2020年國(guó)家芝麻新品種展示試驗(yàn)（黃淮片區(qū)），周芝11號(hào)在3個(gè)試驗(yàn)點(diǎn)平均單產(chǎn)1 271.83 kg/hm2，比豫芝4號(hào)增產(chǎn)6.65%，增產(chǎn)點(diǎn)次100%，居參試品種第3位。由此表明，周芝11號(hào)是一個(gè)具有較大增產(chǎn)潛力、高產(chǎn)的優(yōu)良品種。

2.4 穩(wěn)產(chǎn)性分析

由表2可知，在2年黃淮區(qū)芝麻新品種區(qū)域試驗(yàn)和國(guó)家芝麻新品種展示試驗(yàn)（黃淮流域）中，周芝11號(hào)變異系數(shù)分別為18.19%、15.87%、5.32%，平均變異系數(shù)為13.13%，相對(duì)較低，說(shuō)明周芝11號(hào)產(chǎn)量受環(huán)境影響不大，穩(wěn)產(chǎn)性較好。周芝11號(hào)HSC值分別為83.35%、84.79%、91.80%，在參試品系中分別排第3、2、2位，平均高穩(wěn)系數(shù)為86.65%，相對(duì)較高，充分說(shuō)明該品種遺傳基礎(chǔ)較好。

2.5 適應(yīng)性分析

由表2可知，2018年、2019年黃淮流域芝麻新品種區(qū)域試驗(yàn)和2020年國(guó)家芝麻新品種展示試驗(yàn)（黃淮流域）中，周芝11號(hào)回歸系數(shù)分別為1.358 6、1.143 9、0.933 3，平均回歸系數(shù)為1.145 3，趨近于1，說(shuō)明周芝11號(hào)具有廣泛的適應(yīng)性。

表2 周芝11號(hào)豐產(chǎn)性、穩(wěn)產(chǎn)性和適應(yīng)性分析

2.6 品質(zhì)分析

周芝11號(hào)籽粒純白，經(jīng)農(nóng)業(yè)農(nóng)村部農(nóng)產(chǎn)品質(zhì)量監(jiān)督檢驗(yàn)測(cè)試中心（鄭州）測(cè)定（表3），周芝11號(hào)平均粗脂肪含量55.4%、平均粗蛋白質(zhì)含量20.9%。

表3 2018-2019年區(qū)域試驗(yàn)周芝11號(hào)品質(zhì)分析

2.7 抗性表現(xiàn)

2018年、2019年在芝麻枯萎病和莖點(diǎn)枯病病圃進(jìn)行人工接種鑒定（表4），周芝11號(hào)2018年枯萎病和莖點(diǎn)枯病的病情指數(shù)分別為5.07和15.68，高抗枯萎病和莖點(diǎn)枯??；2019年枯萎病和莖點(diǎn)枯病的病情指數(shù)分別為10.09和44.58，抗枯萎病，高感莖點(diǎn)枯病；枯萎病和莖點(diǎn)枯病2年平均病情指數(shù)分別為7.58和30.13，周芝11號(hào)高抗枯萎病、感莖點(diǎn)枯病。

3 討論與結(jié)論

在高產(chǎn)穩(wěn)產(chǎn)芝麻新品種的選育過(guò)程中，品種區(qū)域試驗(yàn)是必不可少的環(huán)節(jié)，該環(huán)節(jié)通過(guò)對(duì)產(chǎn)量的構(gòu)成要素、品質(zhì)及抗性進(jìn)行分析評(píng)價(jià)，綜合各性狀進(jìn)行全面鑒定[21，24]，從而評(píng)價(jià)參試品種的豐產(chǎn)性、穩(wěn)定性和適應(yīng)性[21]，是檢驗(yàn)品種是否優(yōu)良的基本依據(jù)[22-23]。其中，用目標(biāo)品種與對(duì)照品種的產(chǎn)量相比較的方法分析目標(biāo)品種的豐產(chǎn)性，具有直接、簡(jiǎn)單、明了、清晰的優(yōu)點(diǎn)，但該方法不能反映目標(biāo)品種的穩(wěn)產(chǎn)性。而目標(biāo)品種產(chǎn)量的標(biāo)準(zhǔn)差和變異系數(shù)可以反映出品種的穩(wěn)產(chǎn)性，但不能分析其豐產(chǎn)性。因此本研究選用計(jì)算簡(jiǎn)便、結(jié)果易于量化和排序的高穩(wěn)系數(shù)來(lái)綜合分析新品種豐產(chǎn)性和穩(wěn)產(chǎn)性，選用能體現(xiàn)環(huán)境對(duì)品種產(chǎn)量影響的回歸系數(shù)反映品種的適應(yīng)性[25]。品種的穩(wěn)產(chǎn)性主要以變異系數(shù)（CV）和高穩(wěn)系數(shù)（HSC）來(lái)衡量，其中變異系數(shù)主要衡量品種穩(wěn)產(chǎn)性，變異系數(shù)越小，反映品種的產(chǎn)量在不同環(huán)境下變化越小，表明該品種穩(wěn)產(chǎn)性越好[26]；高穩(wěn)系數(shù)衡量豐產(chǎn)性，并兼顧穩(wěn)產(chǎn)性，高穩(wěn)系數(shù)越大，表明品種豐產(chǎn)穩(wěn)產(chǎn)性越好。作物產(chǎn)量是遺傳因素和外部環(huán)境共同作用的結(jié)果，以各試驗(yàn)點(diǎn)的平均產(chǎn)量與整個(gè)區(qū)試總平均產(chǎn)量的差值（環(huán)境指數(shù)）為自變量，以每個(gè)參試品種在相應(yīng)各試驗(yàn)點(diǎn)的產(chǎn)量為依變量計(jì)算回歸系數(shù)b，以回歸系數(shù)b反映品種對(duì)不同環(huán)境的敏感程度。當(dāng)b=1時(shí)，表示該品種具有廣泛適應(yīng)性；當(dāng)b<1時(shí)表示該品種高于平均適應(yīng)性；當(dāng)b>1時(shí)，表示該品種的適應(yīng)性低于平均適應(yīng)性[27]。

本研究應(yīng)用高通量測(cè)序技術(shù)、變異系數(shù)、高穩(wěn)系數(shù)、回歸系數(shù)、人工接種鑒定等方法，對(duì)芝麻新品種周芝11號(hào)的物種特異性、產(chǎn)量、品質(zhì)及抗性進(jìn)行了綜合分析與評(píng)價(jià)?；蚪M信息比對(duì)表明周芝11號(hào)不含有轉(zhuǎn)基因成分，且符合芝麻同物種基因組信息要求、參試品種特異性要求。區(qū)域試驗(yàn)、展示試驗(yàn)（黃淮流域）中，芝麻新品種周芝11號(hào)的平均單產(chǎn)分別為1 415.82 kg/hm2、1 271.83 kg/hm2， 較豫芝4號(hào)（對(duì)照）分別增產(chǎn)11.41%、6.65%。在區(qū)域試驗(yàn)和展示試驗(yàn)（黃淮流域）中，周芝11號(hào)平均變異系數(shù)為13.13%，較各參試品系相對(duì)較低，說(shuō)明周芝11號(hào)產(chǎn)量受環(huán)境影響不大，穩(wěn)產(chǎn)性好；平均高穩(wěn)系數(shù)為86.65%，較各參試品系相對(duì)較高，說(shuō)明該品種遺傳基礎(chǔ)較好；回歸系數(shù)平均值為1.145 3，接近1，說(shuō)明周芝11號(hào)具有廣泛的適應(yīng)性。經(jīng)人工接種鑒定，周芝11號(hào)高抗枯萎病、感莖點(diǎn)枯病。

綜合分析表明，周芝11號(hào)主要產(chǎn)量要素協(xié)調(diào)、穩(wěn)定，年際間變化小，是一個(gè)綜合性狀好、豐產(chǎn)性、穩(wěn)產(chǎn)性和適應(yīng)性強(qiáng)的芝麻新品種，可在黃淮流域芝麻產(chǎn)區(qū)大面積推廣種植，但在苗期和盛花期應(yīng)及時(shí)做好莖點(diǎn)枯病的防治工作。