甘藍(lán)型油菜EMS誘變矮稈突變體分析

李 昕,閻晉東,楊 飄,向芙江,張 維,彭武生,卓宇紅,李新梅,趙小英*

(1.湖南大學(xué)生物學(xué)院植物功能基因組學(xué)與發(fā)育調(diào)控湖南省重點(diǎn)實(shí)驗(yàn)室,中國(guó)湖南長(zhǎng)沙410082；2.湖南大學(xué)深圳研究院,中國(guó)廣東深圳518057；3.湖南亞華種業(yè)科學(xué)研究院,中國(guó)湖南長(zhǎng)沙410001；4.長(zhǎng)沙市豐景大地農(nóng)業(yè)科技有限公司,中國(guó)湖南長(zhǎng)沙410001)

油菜是我國(guó)最主要的油料作物之一,也是我國(guó)最主要的飼用蛋白質(zhì)作物之一[1],屬于十字花科(Cruciferae)蕓薹屬(Brassica)。甘藍(lán)型油菜(Brassica napus)是油菜中最重要的栽培種。常規(guī)的油菜品種植株較高,并且隨著甘藍(lán)型油菜雜交品種的廣泛推廣,油菜株高平均增加了20 cm以上[2]。高稈油菜倒伏后對(duì)油菜的生長(zhǎng)發(fā)育影響很大,可導(dǎo)致減產(chǎn)15%～30%,嚴(yán)重時(shí)可達(dá)到60%以上[3]。油菜倒伏問題已經(jīng)成為制約油菜機(jī)械化生產(chǎn)和產(chǎn)量提高的重要因素[4]。

植物矮稈基因資源的開發(fā)和創(chuàng)造是矮化品種選育的重要基礎(chǔ)。研究表明,許多植物的矮化突變與植物激素赤霉素(gibberellin,GA)和油菜素甾醇(brassinosteroid,BR)有關(guān),而少數(shù)植物的矮化突變與生長(zhǎng)素(auxin)有關(guān)[5]。因此,植物矮稈基因也是研究植物激素生物合成和信號(hào)傳導(dǎo)及植物生長(zhǎng)發(fā)育的重要資源[6]。目前,可通過自然變異、人工誘變、雜交和生物技術(shù)等方法來產(chǎn)生和培育植物矮稈突變體,從而獲得植物矮稈基因資源[4,6～9]。

課題組前期利用甲基磺酸乙酯(ethylmethylsulfone,EMS)(0.8%)溶液處理甘藍(lán)型油菜玻里馬胞質(zhì)雄性不育系(pol cms)的保持系2B種子,篩選獲得多個(gè)矮稈突變體,將其命名為Bnd1(Brassica napus dwarf 1)、Bnd2、Bnd3 等[10]。本研究對(duì)其中 Bnd2突變體的農(nóng)藝性狀進(jìn)行了分析；利用GA和BR對(duì)突變體進(jìn)行處理,分析了突變體對(duì)植物激素的響應(yīng)；同時(shí),通過雜交對(duì)矮稈性狀進(jìn)行了遺傳分析,為油菜矮稈突變體基因的定位克隆研究提供了線索。

1 材料與方法

1.1 材料

本研究用到的野生型油菜為甘藍(lán)型油菜玻里馬胞質(zhì)雄性不育系(pol cms)的保持系2B,矮稈突變體Bnd2通過EMS誘變2B種子篩選獲得[10]。

1.2 田間種植與管理

2017年9月于湖南寧鄉(xiāng)關(guān)山油菜基地同時(shí)種植誘變親本2B、矮稈突變體Bnd2、2B與Bnd2的回交F1代和F2群體。對(duì)所有實(shí)驗(yàn)材料進(jìn)行點(diǎn)播育苗,待幼苗長(zhǎng)出兩片真葉后,根據(jù)幼苗的大小進(jìn)行分級(jí)移栽。試驗(yàn)基地每小區(qū)寬2 m、長(zhǎng)20 m、行距33 cm,每行種植10株。

1.3 農(nóng)藝性狀鑒定

2018年4月,于油菜成熟期對(duì)矮稈突變體Bnd2的農(nóng)藝性狀進(jìn)行測(cè)量統(tǒng)計(jì),包括株高、有效分枝高度、一次有效分枝數(shù)、主花序有效長(zhǎng)度、主花序角果數(shù)、單株有效角果數(shù)、角果長(zhǎng)度、每角粒數(shù)、節(jié)間距、伸長(zhǎng)節(jié)數(shù)。隨后,收取矮稈Bnd2和親本2B植株各10株,自然風(fēng)干,用于測(cè)量千粒重和單株產(chǎn)量。

1.4 矮稈性狀遺傳分析

2016年9月在湖南寧鄉(xiāng)關(guān)山油菜基地播種親本2B和矮稈突變體Bnd2。2017年3月以2B為母本、突變體Bnd2為父本進(jìn)行回交,獲得回交F1代種子；5月份將F1代種子在青海加代繁殖,套袋自交獲得F2代種子。同年9月,在湖南寧鄉(xiāng)關(guān)山油菜基地播種親本2B、突變體Bnd2,以及F1和F2代種子,于次年成熟期測(cè)量株高。

1.5 植物激素處理

對(duì)外源激素赤霉素(GA3,購(gòu)自國(guó)藥集團(tuán)化學(xué)試劑有限公司)和油菜素甾醇(BR,購(gòu)自美侖生物技術(shù)有限公司)設(shè)定不同的濃度梯度。GA3的濃度為0 mol/L、10-7mol/L、10-6mol/L、10-5mol/L、10-4mol/L和10-3mol/L。BR的濃度為0 mol/L、10-8mol/L、10-7mol/L、10-6mol/L和10-5mol/L。取2B和Bnd2種子各30粒,置于同一鋪滿3層濾紙的培養(yǎng)瓶中,瓶中添加1/2 MS液體培養(yǎng)基10 mL(含不同濃度的GA3或BR),每個(gè)處理重復(fù)3次,播種7 d后測(cè)量各處理幼苗的下胚軸長(zhǎng)度,并拍照。

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

2 結(jié)果

2.1 矮稈突變體Bnd2的農(nóng)藝性狀分析

對(duì)田間種植的親本和矮稈突變體油菜的農(nóng)藝性狀進(jìn)行統(tǒng)計(jì)分析,結(jié)果顯示,矮稈突變體Bnd2較誘變親本2B顯著變矮(圖1,表1)。其中,誘變親本2B的株高為168.2±7.61 cm；矮稈突變體Bnd2的株高為100.65±8.09 cm。同時(shí),Bnd2的有效分枝高度(24.95±3.17 cm)、伸長(zhǎng)節(jié)數(shù)(13.6±2.11)、平均節(jié)間距(6.55±1.46 cm)、主花序有效長(zhǎng)度(52.91±7.72 cm)相對(duì)于誘變親本2B都明顯減少(表1),表明Bnd2株高降低主要是由于其有效分枝高度降低、伸長(zhǎng)節(jié)數(shù)減少、節(jié)間距和主花序縮短。此外,Bnd2的主花序角果數(shù)、一次有效分枝數(shù)、每角粒數(shù)、角果長(zhǎng)度、千粒重和單株產(chǎn)量也顯著低于2B,但單株有效角果數(shù)在兩者之間無顯著差異(表1)。

圖1 誘變親本2B和矮稈突變體Bnd2的植株表型(A)120 d苗齡的植株；(B)成熟期植株。圖中標(biāo)尺均代表10 cm。Fig.1 Phenotypes of mutagenic parent 2B and dwarf mutant Bnd2(A)120-day-old plants；(B)Plants at maturity.Size bars in(A)and(B)represent 10 cm.

2.2 矮稈突變體Bnd2的矮稈性狀遺傳分析

為了確定矮稈突變體Bnd2的基因型,我們將誘變親本2B、矮稈突變體Bnd2以及2B和Bnd2的回交F1和F2代,同時(shí)播種在湖南寧鄉(xiāng)關(guān)山基地,于成熟期測(cè)量株高。株高的統(tǒng)計(jì)結(jié)果如圖2所示。親本2B的株高主要分布在150～169 cm(圖2A)；矮稈突變體Bnd2的株高主要分布在90～109 cm(圖2B)；回交F1代株高主要分布在135～144 cm(圖2C),偏向于高稈,說明矮稈突變體Bnd2為隱性突變。

F2代群體236個(gè)單株株高的統(tǒng)計(jì)結(jié)果如圖2D所示。從圖片數(shù)據(jù)可以看出,在2B和Bnd2回交F2代群體中株高發(fā)生了分離,呈明顯的雙峰分布。在拐點(diǎn)附近取臨界值110 cm進(jìn)行卡方檢驗(yàn),結(jié)果發(fā)現(xiàn)群體中高稈與矮稈的分離比符合3︰1的理論值,表明Bnd2突變體的矮稈性狀受單基因控制,且為隱性突變。

2.3 矮稈突變體Bnd2幼苗下胚軸對(duì)外源GA、BR的響應(yīng)

為了確認(rèn)Bnd2突變體的矮化是否與GA和BR信號(hào)通路有關(guān),我們檢測(cè)了Bnd2突變體幼苗下胚軸伸長(zhǎng)對(duì)外源GA和BR的響應(yīng)。將誘變親本2B和突變體Bnd2播種在含有不同濃度GA3或者不同濃度BR的1/2 MS培養(yǎng)基中,置于連續(xù)白光下培養(yǎng),7 d后測(cè)量幼苗下胚軸的長(zhǎng)度。從圖3可看出,Bnd2幼苗下胚軸比誘變親本2B的下胚軸要短,對(duì)GA3的響應(yīng)減弱,并且GA3處理不能恢復(fù)Bnd2的表型(圖3A～C)；Bnd2突變體對(duì)BR的響應(yīng)則增強(qiáng),且BR可恢復(fù)Bnd2突變體幼苗的表型(圖3D～F)。根據(jù)這些研究結(jié)果,我們推測(cè)在Bnd2突變體中GA信號(hào)通路和/或者BR合成受到影響,這為后期基因定位克隆研究提供了線索。

表1 矮稈突變體Bnd2的農(nóng)藝性狀統(tǒng)計(jì)Table 1 Agronomic characters of dwarf mutant Bnd2(±s)

表1 矮稈突變體Bnd2的農(nóng)藝性狀統(tǒng)計(jì)Table 1 Agronomic characters of dwarf mutant Bnd2(±s)

注：PH,株高；BH,有效分枝高度；NPB,一次有效分枝數(shù)；LR,主花序有效長(zhǎng)度；NSR,主花序角果數(shù)；NSP,單株有效角果數(shù)；LS,角果長(zhǎng)度；NSS,每角粒數(shù)；MIL,平均節(jié)間距；NEI,伸長(zhǎng)節(jié)數(shù)；TSW,千粒重；SYP,單株產(chǎn)量。：平均值；s：標(biāo)準(zhǔn)偏差；n=10；*：0.01＜P≤0.05,結(jié)果達(dá)到顯著水平；**：P≤0.01,結(jié)果達(dá)到極顯著水平。Notes：PH,plant height；BH,effective branch height；NPB,number of effective primary branches；LR,effective length of raceme；NSR,number of siliques on raceme；NSP,effective number of siliques per plant；LS,length of siliques；NSS,number of seeds per silique；MIL,mean internode length；NEI,number of elongated internode；TSW,thousand seed weight；SYP,seed yield per plant.：Mean；s：Standard deviation；n=10；*：0.01＜P≤0.05,the result reaches the significant level；**：P≤0.01,the result reaches the extremely significant level.

Mutagenic parent 2B Dwarf mutant Bnd2 Mutagenic parent 2B Dwarf mutant Bnd2 PH/cm 168.20±7.61 100.65±8.09**BH/cm 60.81±9.68 24.95±3.17**NPB 7.60±1.43 5.80±1.08*LR/cm 62.06±4.91 52.91±7.72*NSR 61.10±3.88 49.50±10.29*NSP 255.70±69.93 169.80±57.55 LS/cm 8.28±0.38 6.88±0.37**NSS 27.81±3.12 25.57±3.64**MIL/cm 8.60±0.71 6.55±1.46**NEI 19.70±1.85 13.60±2.11**TSW/g 3.36±0.08 2.67±0.18**SYP/g 18.81±6.15 9.13±1.98**

圖2 誘變親本2B(A)、突變體Bnd2(B)、回交F1代(C)及F2代(D)群體植株的株高頻率分布Fig.2 Plant height frequency distribution of mutagenic parent 2B(A),dwarf mutant Bnd2(B),backcross F1generation(C)and F2population(D)

3 討論

我國(guó)是世界上最大的油菜生產(chǎn)國(guó),油菜種植面積占世界油菜總生產(chǎn)面積的1/3[11～12]。隨著雜交油菜的利用,油菜的品質(zhì)得到顯著提高,但是雜交油菜的植株高大易倒伏,易造成油菜收獲指數(shù)的降低[2,13～14]。為了提高雜交油菜的抗倒伏性,科學(xué)家們通過誘變手段創(chuàng)制了一系列矮稈突變體。石淑穩(wěn)等[15～16]利用EMS技術(shù)獲得了兩株矮稈植株,命名為DS-1和DS-2,株高分別為106 cm和95 cm。Wang等[17]通過EMS誘變獲得一個(gè)矮稈突變體Bndwf1,該突變體苗期葉片深綠,株高僅有80～110 cm,顯著低于野生型。本研究利用EMS處理甘藍(lán)型油菜保持系2B種子,篩選獲得矮稈突變體Bnd2,其株高為100.65±8.09 cm,顯著低于親本株高(表1)。遺傳分析結(jié)果表明Bnd2的矮稈性狀為隱性性狀,且受單基因控制(圖2),這為油菜抗倒伏育種提供了新的矮稈資源。

GA、BR和生長(zhǎng)素是促進(jìn)植物生長(zhǎng)、控制株高的重要激素。目前,已報(bào)道的油菜矮稈突變體的矮化大多數(shù)與GA或者生長(zhǎng)素信號(hào)減弱有關(guān)[5,18～23]。Liu等[18]研究發(fā)現(xiàn)ds-1半矮稈突變體的DS-1基因編碼GA信號(hào)抑制子BnaA06.RGA,在突變體中該蛋白質(zhì)VHYNP模體中保守的脯氨酸(pro)突變?yōu)榱涟彼?leu)。Zhao等[19]鑒定到一個(gè)與ds-1表型相似的半矮稈突變體ds-3,DS-3基因也編碼一個(gè)在GA信號(hào)途徑中起抑制作用的DELLA蛋白(BnaC07.RGA),并且與BnaA06.rga-ds突變位點(diǎn)一樣,BnaC07.rga-ds的保守區(qū)域VHYNP中發(fā)生了一個(gè)相同的氨基酸替換(脯氨酸突變?yōu)榱涟彼?,導(dǎo)致GA信號(hào)途徑受阻,植株矮化。另外,Li和Zeng等[20～21]分別鑒定到的矮稈突變體NDF-1與bnaC.dwf,都是GA不敏感型突變體,其中NDF-1矮化是由于GA受體BnGID1的編碼基因啟動(dòng)子發(fā)生突變,導(dǎo)致基因表達(dá)降低,從而使GA信號(hào)減弱；而bnaC.dwf突變體受一對(duì)隱性基因控制,株高對(duì)外源赤霉素處理不敏感,但酶聯(lián)免疫吸附測(cè)定(enzyme-linked immunosorbent assay,ELISA)結(jié)果表明,突變體葉片與莖中內(nèi)源赤霉素含量顯著高于野生型T6(+/+),說明該突變體在GA信號(hào)轉(zhuǎn)導(dǎo)途徑中存在缺陷,但控制該性狀的候選基因目前尚未確定。最近,Cheng等[22]通過組織培養(yǎng)獲得一個(gè)矮稈突變體G7,該突變體中Bna.IAA7.C05基因所編碼蛋白質(zhì)的保守區(qū)域GWPPV突變?yōu)镋WPPV,使得生長(zhǎng)素響應(yīng)相關(guān)基因的表達(dá)下調(diào),從而導(dǎo)致植株變矮。Li等[23]鑒定了一個(gè)EMS誘變的功能獲得型半矮稈突變體sca,進(jìn)一步的研究發(fā)現(xiàn),SCA基因編碼生長(zhǎng)素信號(hào)負(fù)調(diào)控因子Aux/IAA(BnaA3.IAA7),在突變體中該蛋白質(zhì)第84位甘氨酸(G)突變?yōu)楣劝彼?E),突變蛋白質(zhì)的穩(wěn)定性增加,使得生長(zhǎng)素信號(hào)受到抑制,從而使植株矮化。本研究發(fā)現(xiàn),Bnd2突變體幼苗下胚軸對(duì)GA的響應(yīng)減弱,對(duì)BR的響應(yīng)增強(qiáng),而且BR可以恢復(fù)突變體幼苗下胚軸的表型(圖3),因此推測(cè)在Bnd2突變體中GA信號(hào)和/或者BR合成受到影響。后期對(duì)該突變位點(diǎn)的定位克隆,將對(duì)闡明GA/BR調(diào)控油菜株型和生長(zhǎng)發(fā)育的機(jī)制,以及培育矮稈油菜品種具有重要意義。

圖3 矮稈突變體Bnd2幼苗下胚軸的伸長(zhǎng)生長(zhǎng)對(duì)外源GA3和BR的響應(yīng)(A)對(duì)照組(-GA)及10-4mol/L GA3處理組(+GA)的2B與Bnd2幼苗表型；(B)不同濃度GA3處理后2B和Bnd2的下胚軸長(zhǎng)度；(C)不同濃度GA3處理后2B和Bnd2中下胚軸長(zhǎng)度與對(duì)照組(0 mol/L)下胚軸長(zhǎng)度的比值；(D)對(duì)照組(-BR)及10-5mol/L BR處理組(+BR)的2B與Bnd2幼苗表型；(E)不同濃度BR處理后2B和Bnd2的下胚軸長(zhǎng)度；(F)不同濃度BR處理后2B和Bnd2中下胚軸長(zhǎng)度與對(duì)照組(0 mol/L)下胚軸長(zhǎng)度的比值。圖(B)～(C)與(E)～(F)的數(shù)據(jù)代表平均值±標(biāo)準(zhǔn)差,n=20；*：0.01＜P≤0.05,結(jié)果達(dá)到顯著水平；**：P≤0.01,結(jié)果達(dá)到極顯著水平。圖(A)、(D)中的標(biāo)尺均代表10 mm。Fig.3 Responses of hypocotyl elongation of dwarf mutant Bnd2 seedlings to exogenous GA3and BR(A)Phenotypes of 2B and Bnd2 seedlings in control group(-GA)and treatment group with 10-4mol/L GA3(+GA)；(B)Hypocotyl lengths of 2B and Bnd2 treated with different concentrations of GA3；(C)Ratio of hypocotyl length of 2B and Bnd2 treated with different concentrations of GA3to that of the control group(0 mol/L)；(D)Phenotypes of 2B and Bnd2 seedlings in control group(-BR)and treatment group with 10-5mol/L BR(+BR)；(E)Hypocotyl lengths of 2B and Bnd2 treated with different concentrations of BR；(F)Ratio of hypocotyl length of 2B and Bnd2 treated with different concentrations of BR to that of the control group(0 mol/L).Data shown in(B)～(C)and(E)～(F)are represented as mean ± SD,n=20；*：0.01＜P≤0.05,the result reaches the significant level；**：P≤0.01,the result reaches the extremely significant level.Size bars shown in(A)and(D)are represented as 10 mm.