水稻雄性核不育突變體ms7的遺傳分析及基因定位

楊晉宇　白琛　丁小惠　申紅芳　王磊　應(yīng)杰政, *　鄂志國, *

水稻雄性核不育突變體的遺傳分析及基因定位

楊晉宇1, #白琛2, 3, #丁小惠2, 3申紅芳1王磊1應(yīng)杰政1, *鄂志國1, *

（1中國水稻研究所 水稻生物學(xué)國家重點(diǎn)實(shí)驗(yàn)室, 杭州 310006;2中國科學(xué)院 植物研究所 植物分子生理學(xué)重點(diǎn)實(shí)驗(yàn)室, 北京 100093;3中國科學(xué)院大學(xué)，北京 100049;#共同第一作者;*通信聯(lián)系人, E-mail: yingjiezheng@caas.cn; ezhiguo@caas.cn）

【】通過對(duì)水稻雄性不育突變體的研究，可以鑒定更多與育性或花粉發(fā)育相關(guān)的基因，有助于解析水稻雄性生殖發(fā)育的整個(gè)調(diào)控網(wǎng)絡(luò)。常規(guī)種植條件下，突變體()與對(duì)照種植于浙江富陽和海南陵水，比較它們的育性及主要農(nóng)藝性狀差異，利用混池關(guān)聯(lián)分析和圖位克隆方法進(jìn)行目標(biāo)基因定位。整個(gè)生育期，突變體生長(zhǎng)速率與野生型一致，成熟期的株高、分蘗數(shù)、葉數(shù)、葉大小、穗長(zhǎng)和每穗穎花數(shù)等性狀與野生型相比也沒有顯著差異，但結(jié)實(shí)率為0，表現(xiàn)為完全雄性不育，花藥瘦小且顏色發(fā)白，半薄切片顯示絨氈層降解推遲，花粉鏡檢呈染敗。遺傳分析表明花粉敗育受單個(gè)隱性基因控制，定位于第7染色體上BSA11與YD7045之間1.17 Mb的范圍內(nèi)。本研究為水稻雄性不育基因的克隆和功能研究打下了基礎(chǔ)。

水稻；雄性不育；混池關(guān)聯(lián)分析；基因定位

水稻是我國乃至東南亞最重要的糧食作物，在我國有著悠久的種植和馴化歷史[1]。水稻細(xì)胞核雄性不育系的發(fā)現(xiàn)及應(yīng)用，推動(dòng)了雜交水稻從三系到兩系的發(fā)展，為我國乃至世界糧食安全發(fā)揮了重要作用[2]。目前生產(chǎn)中利用的水稻核不育系多是光溫敏雄性核不育系，相應(yīng)的控制基因如[3]、/[4-5]、等也已鑒定。近年發(fā)現(xiàn)一種新型濕敏雄性核不育現(xiàn)象，突變體在低濕環(huán)境(相對(duì)濕度<60%)下因花粉不能有效防止脫水而無法完成粘附和水合作用，造成不育，而高濕環(huán)境(相對(duì)濕度>80%)下恢復(fù)育性[8-10]。除這種環(huán)境和遺傳因子互作外，還有諸多影響小孢子和花粉發(fā)育的因素，而絨氈層發(fā)生和降解異常是重要原因之一。絨氈層位于花藥壁最內(nèi)層，通過細(xì)胞程序化死亡(programmed cell death, PCD)，為小孢子的發(fā)生發(fā)育提供營養(yǎng)，分泌脂質(zhì)參與花粉外壁建成和花藥角質(zhì)層合成[11]。F-box蛋白OsADF、半胱氨酸蛋白酶OsCP1[13]、天冬氨酸蛋白酶OsAP25和OsAP37[14]等，可通過蛋白水解通路直接參與絨氈層的PCD脂肪酸還原酶DPW[15]、α整合素樣蛋白DPW3[16]、?；o酶A合成酶[17]、?；D(zhuǎn)移酶DPW2[18]和?；D(zhuǎn)移酶OsGPAT3[19]，通過參與脂質(zhì)代謝，影響絨氈層PCD和花粉外壁前體物質(zhì)合成。此外，鈣離子結(jié)合蛋白Ca。OsAGO2通過DNA甲基化直接下調(diào)己糖激酶基因表達(dá)，影響絨氈層活性氧積累，敲低或過表達(dá)都能導(dǎo)致PCD提早和花粉敗育[21]。

在絨氈層發(fā)育調(diào)控網(wǎng)絡(luò)中，MYB轉(zhuǎn)錄因子GAMYB[22-23]和BM1[24]，bHLH轉(zhuǎn)錄因子UDT1[25]、TIP2[26-27]、TDR[28]和EAT1/DTD[14, 29]，E3泛素連接酶OsHUB1和OsHUB2[30]，PHD鋅指蛋白PTC1/OsMS1[31-33]和凋亡抑制因子API5[34]等都起重要作用。TDR作用于GAMYB和UDT1下游，通過直接調(diào)節(jié)、和等基因的表達(dá)，激活蛋白水解通路從而引發(fā)絨氈層降解[12, 23, 28]。EAT1/DTD作用于TDR下游并能與TDR互作，通過上調(diào)和表達(dá)，正調(diào)控絨氈層PCD[14, 29]。TIP2/bHLH142能直接調(diào)控和的表達(dá)，且TIP2能與TDR互作，從而調(diào)控絨氈層PCD[26-27]。OsMS1和TIP3都具有轉(zhuǎn)錄激活活性，它們分別通過與TIP2和TDR互作，調(diào)控烏氏體發(fā)生和絨氈層降解[32-33]。OsHUB1和OsHUB2作為E3連接酶，通過組蛋白H2B單泛素化，正調(diào)控、和表達(dá)，影響絨氈層降解[30]。ATP檸檬酸裂合酶ACL催化細(xì)胞質(zhì)中的檸檬酸生成草酰乙酸和乙酰輔酶A，通過影響、、和等基因的表達(dá)，調(diào)控絨氈層PCD。ACL突變體絨氈層細(xì)胞過早降解，花藥中ATP水平和脂肪酸含量降低，不能形成有活性的花粉[35]。這些表明，在花藥壁的分化、發(fā)生和降解過程中，TIP2、TDR和EAT1這3個(gè)bHLH轉(zhuǎn)錄因子可能起核心調(diào)控作用。編碼一種AAA-ATP酶，與RAD51和DMC1互作，參與減數(shù)分裂的染色體雙鏈斷裂修復(fù)，但突變體絨氈層降解異常，是導(dǎo)致花藥異常和小孢子敗育的原因之一[36]。此外，生殖細(xì)胞也能調(diào)節(jié)絨氈層PCD。編碼一種糖蛋白，被小孢子分泌到細(xì)胞外與絨氈層細(xì)胞表面受體互作，調(diào)控絨氈層PCD和小孢子發(fā)育[37]。

我們?cè)趯?duì)秈稻中恢161進(jìn)行組培實(shí)驗(yàn)時(shí)，從其后代中鑒定出一個(gè)雄性不育突變體，也表現(xiàn)為絨氈層降解滯后和花粉敗育，其他表型在整個(gè)生育期與野生型沒有顯著差異。進(jìn)一步用潮霉素基因的分子標(biāo)記對(duì)后代不育個(gè)體進(jìn)行PCR鑒定，在無篩選標(biāo)記的個(gè)體中依然表現(xiàn)為不育，表明不育性狀并非由遺傳轉(zhuǎn)化載體結(jié)構(gòu)T-DNA插入而產(chǎn)生，推測(cè)可能是組培過程中產(chǎn)生了新的遺傳變異導(dǎo)致中恢161細(xì)胞核雄性不育。為了明確雄性不育產(chǎn)生的原因，我們用作母本，IR96671為父本構(gòu)建遺傳群體進(jìn)行遺傳分析，結(jié)果表明的不育性狀受單個(gè)隱性核基因控制。經(jīng)混池關(guān)聯(lián)分析（bulked segregant analysis, BSA）和基因定位，將定位于第7染色體上BSA11－YD7045之間，這些結(jié)果為進(jìn)一步的基因克隆和功能分析打下了基礎(chǔ)。

1?材料與方法

1.1?供試材料

本研究采用的突變體是秈型恢復(fù)系中恢161經(jīng)組織培養(yǎng)產(chǎn)生的，該突變體在杭州富陽和海南陵水均表現(xiàn)為雄性完全不育，不受溫度和日照影響。常規(guī)秈稻IR96671作為雜交親本構(gòu)建遺傳分析和定位的群體。

1.2?群體構(gòu)建及遺傳分析

突變體為母本，秈稻品種IR96671為父本，配制雜交組合，觀察F1的表型，分單株收獲種子。在F2群體中統(tǒng)計(jì)野生型和不育型的分離比例，用于遺傳分析、BSA分析和基因定位。所有材料2018?2019年種植于中國水稻研究所杭州富陽實(shí)驗(yàn)基地和海南陵水基地，正常田間管理。

1.3?花粉育性調(diào)查

抽穗期當(dāng)日開花前，田間取回野生型和不育型的穎花，在顯微鏡下觀察花藥表型并拍照。用1% I2-KI溶液進(jìn)行染色，觀察野生型和不育型花粉染色的情況并拍照。

灌漿后期觀察單株灌漿結(jié)實(shí)情況，記錄植株育性。

1.4?花藥半薄切片觀察

取野生型和不育型不同發(fā)育時(shí)期的小花，在FAA固定液（甲醛∶冰乙酸∶50%乙醇=5∶5∶90）中固定，抽真空至樣品全部沉底，換用新的FAA常溫保存?zhèn)溆?。挑選不同時(shí)期的小花剝?nèi)シf殼露出花藥后用酒精進(jìn)行梯度脫水，接著在飽和番紅水溶液中處理2～3 h，再使用Technovit 7100試劑盒進(jìn)行滲透包埋，60℃恒溫烘箱中烘烤4 d。修塊后用半薄切片機(jī)（Leica RM2265）橫切花藥，切片厚度為2 μm，經(jīng)甲苯胺藍(lán)染色后用光學(xué)顯微鏡觀察拍照。

1.5?全基因組測(cè)序

在×IR96671的F2群體中，根據(jù)育性的分離分為可育型和不育型兩組，從中分別挑選出30株，采用簡(jiǎn)易DNA提取法分單株提取葉片DNA，檢測(cè)每份DNA的品質(zhì)與濃度后等量混合制備混池。利用超聲波將DNA序列片段化形成隨機(jī)片段，對(duì)片段化的DNA進(jìn)行末端修復(fù)、3′端加A、連接測(cè)序接頭后，富集400 bp左右的隨機(jī)片段，經(jīng)PCR擴(kuò)增形成測(cè)序文庫。文庫經(jīng)質(zhì)檢合格后，用Illumina HiSeqTM平臺(tái)進(jìn)行測(cè)序，測(cè)序策略為Illumina PE150，總測(cè)序讀長(zhǎng)為300 bp。

對(duì)測(cè)序數(shù)據(jù)進(jìn)行質(zhì)量控制，過濾掉低質(zhì)量數(shù)據(jù)。利用BWA軟件[38]將質(zhì)控后數(shù)據(jù)比對(duì)到參考基因組(IRGSP1.0, http://rice.Plantbiology.msu.edu/)上，獲得序列的位置。利用GATK軟件對(duì)位置歸屬校正，并進(jìn)行SNP檢測(cè)[39]，同時(shí)利用Samtools軟件的Vcfutils工具進(jìn)行過濾[40]，獲得高質(zhì)量的變異位點(diǎn)。采用Annovar軟件和參考基因組的基因預(yù)測(cè)信息進(jìn)行注釋?；谟H本和突變混池的數(shù)據(jù)特點(diǎn)，計(jì)算SNP指數(shù)進(jìn)行BSA關(guān)聯(lián)分析，定位目標(biāo)突變位點(diǎn)。

1.6?差異位點(diǎn)的關(guān)聯(lián)分析

對(duì)所得到的SNP和InDel位點(diǎn)進(jìn)行親本和突變混池的差異位點(diǎn)開發(fā)，計(jì)算每個(gè)差異SNP位點(diǎn)的指數(shù)，即突變基因型在突變池和野生池中所有基因型中所占的比例。由于突變位點(diǎn)與周圍標(biāo)記的連鎖效應(yīng)，在突變位點(diǎn)附近，SNP指數(shù)更接近于1，其余位點(diǎn)由于高通量測(cè)序的隨機(jī)分布，應(yīng)符合孟德爾分離比例，本研究的突變性狀屬于隱性突變，在F2群體中，突變位點(diǎn)的SNP指數(shù)為1，正常位點(diǎn)為0.5。

為了消除假陽性位點(diǎn)，利用標(biāo)記在基因組上的位置，用滑窗的方法將SNP指數(shù)進(jìn)行擬合，消除由于隨機(jī)擴(kuò)增導(dǎo)致的差異突變。計(jì)算每個(gè)變異位點(diǎn)的指數(shù)，以染色體物理位置為橫坐標(biāo)繪制曲線圖，計(jì)算所有SNP指數(shù)在95%、99%及99.9%置信水平下的閾值，挑選高于閾值的區(qū)域作為候選區(qū)間。

1.7?DNA提取和定位分析

基于BSA分析得到的初步定位區(qū)間，根據(jù)親本材料重測(cè)序結(jié)果，利用Primer 5.0軟件設(shè)計(jì)位于該區(qū)域的InDel標(biāo)記，標(biāo)記分布區(qū)域完全覆蓋BSA分析得到的關(guān)聯(lián)區(qū)間，并向左右各擴(kuò)展數(shù)兆區(qū)間。挑選多態(tài)性標(biāo)記用于驗(yàn)證BSA測(cè)序結(jié)果及篩選F2群體中的重組單株。除大田種植外，F(xiàn)2群體同時(shí)也在光溫培養(yǎng)箱中進(jìn)行水培。大田材料在移栽一周后取回葉片使用簡(jiǎn)易法提取DNA，水培材料在播種10 d后使用簡(jiǎn)易板提法提取葉片DNA。

PCR擴(kuò)增使用北京康為世紀(jì)生物科技有限公司的10 μL反應(yīng)體系：模板DNA 1 μL、10 μmol/L前后引物各0.5 μL、2×Taq MasterMix預(yù)混液5 μL、滅菌水3 μL。擴(kuò)增程序如下：94℃下預(yù)變性2 min；94℃變性30 s，55℃退火30 s（根據(jù)不同引物TM值修改），72℃延伸30 s，共計(jì)30個(gè)循環(huán)；最后72℃延伸2 min。PCR擴(kuò)增產(chǎn)物依據(jù)片段大小采用2%或2.5%瓊脂糖凝膠電泳分離后經(jīng)Gel Red染色檢測(cè)，然后拍照保存統(tǒng)計(jì)帶型。挑出目標(biāo)區(qū)域內(nèi)的重組單株后重新取樣鑒定，對(duì)最終確認(rèn)的重組單株進(jìn)行目標(biāo)區(qū)域的其他標(biāo)記檢測(cè)，移栽到大田，待灌漿完成后記錄育性。

2?結(jié)果與分析

2.1?突變體ms7表型分析

自然生長(zhǎng)條件下，突變體全生育期生長(zhǎng)速率與野生型中恢161一致，灌漿前表型沒有明顯異常，包括株高、分蘗數(shù)、葉數(shù)、葉大小、穗長(zhǎng)和每穗穎花數(shù)等性狀與野生型相比沒有顯著差異（圖1-A）。但不灌漿，結(jié)實(shí)率為0。用野生型花粉授粉后，能正常結(jié)實(shí)，表明雌蕊發(fā)育正常，該突變體為雄性不育（圖1-B）。與野生型相比，花藥瘦小且顏色變淡發(fā)白（圖1-C）。花粉鏡檢發(fā)現(xiàn)的花粉粒無法被1% I2-KI溶液染色或染色不均，為染敗，表明花粉無活性，而野生型花粉粒能全部染成深色（圖1-D）。

2.2?花藥半薄切片分析

根據(jù)張大兵等[41]將花藥發(fā)育所分成的14個(gè)時(shí)期，可通過花藥半薄切片來比較觀察野生型和突變型在花藥發(fā)育過程中存在的差異。在1～8a時(shí)期，突變體的花藥發(fā)育形態(tài)與野生型相比無明顯差異，花粉母細(xì)胞進(jìn)行正常的減數(shù)分裂并形成了二分體(圖2-A、圖2-G)。在8b時(shí)期，小孢子母細(xì)胞繼續(xù)進(jìn)行減數(shù)分裂，形成四分體。此時(shí)，野生型的絨氈層開始液泡化并伴隨著細(xì)胞程序性死亡(圖2-B)，突變體的絨氈層出現(xiàn)大量空泡(圖2-H)。第9時(shí)期，隨著胼胝質(zhì)壁的降解，自由的單倍體小孢子從四分體中釋放了出來，游離在藥室中。野生型的絨氈層中空泡逐漸被吸收，細(xì)胞質(zhì)繼續(xù)濃縮，顏色加深(圖2-C)，突變體的絨氈層仍可見空泡(圖2-I)。第10時(shí)期，小孢子液泡化，野生型的絨氈層進(jìn)一步退化(圖2-D)，突變型的絨氈層卻明顯增厚(圖2-J)。第11時(shí)期，液泡化的小孢子進(jìn)行第一次不對(duì)稱的有絲分裂，產(chǎn)生一個(gè)較小的生殖細(xì)胞和一個(gè)較大的營養(yǎng)細(xì)胞。野生型的絨氈層繼續(xù)降解(圖2-E)，而突變型的絨氈層仍無降解表現(xiàn)(圖2-K)。第12時(shí)期，小孢子內(nèi)的生殖細(xì)胞進(jìn)行了第二次有絲分裂，野生型中成熟的花粉粒由于淀粉粒和油脂的積累，細(xì)胞內(nèi)充滿儲(chǔ)存物，表皮和藥室內(nèi)壁進(jìn)一步降解，絨氈層完全消失(圖2-F)；突變體僅有極少數(shù)小孢子填充有淀粉，但缺少油脂，絕大多數(shù)小孢子表現(xiàn)為形狀不規(guī)則的空泡，無填充物(圖2-L)。

A―灌漿期的植株形態(tài)；B―ms7和野生型的結(jié)實(shí)情況；C―ms7和野生型的穎花形態(tài)；D―ms7和野生型的花粉碘化鉀染色。

Fig. 1. Phenotype ofand its wild type Zhonghui 161.

2.3?突變體ms7雄性不育性狀的遺傳分析

以突變體為母本，IR96671為父本構(gòu)建了×IR96671群體?？疾镕1植株的花粉育性，發(fā)現(xiàn)F1結(jié)實(shí)率與野生型相當(dāng)，表明突變性狀受隱性基因控制。統(tǒng)計(jì)了4個(gè)F2群體的可育型和不育型植株數(shù)，經(jīng)卡方檢驗(yàn)，≈0.43>0.05，正常育性與雄性不育個(gè)體的分離比符合3∶1（表1），說明的雄性不育性狀受單個(gè)隱性核基因控制。

A～F為野生型花藥在8a至12時(shí)期的切片；G～L為突變體花藥在8a至12時(shí)期的切片。

A?F, Wild-type anthers from stages 8a to 12; G?L,anthers from stages 8a to 12.

Ep―表皮；En―內(nèi)壁；T―絨氈層；Dy―二分體；Tds―四分體；Msp―小孢子；BP―二胞花粉；MP―成熟花粉粒；iMP―未成熟花粉粒。標(biāo)尺為20 μm。

Ep, Epidermis; En, Endothecium; T, Tapetum; Dy, Dyad cell; Tds, Tetrads; Msp, Microspore; BP, Biceullar pollen; MP, Mature pollen; iMP, Immature pollen. Bar=20 μm.

圖2?野生型中恢161和突變體不同發(fā)育時(shí)期花藥半薄切片觀察

Fig. 2. Contrast of slice of anther between themutant and its wild type.

2.4?BSA分析初定位ms7雄性不育基因

4個(gè)樣本池測(cè)序并過濾后共獲得87.59 G的數(shù)據(jù)，Q30達(dá)91.79%，樣品平均覆蓋深度43.72X，與參考基因組平均比對(duì)效率為79.49%，基因組覆蓋度96.02%。這些數(shù)據(jù)表明測(cè)序數(shù)據(jù)充足，質(zhì)量高，可用于下一步分析。

SNP分析表明，4個(gè)樣本共獲得3 905 483個(gè)SNP位點(diǎn)，其中非同義突變的SNP有1 482 935個(gè)。SNP信息統(tǒng)計(jì)顯示，不同樣品的轉(zhuǎn)換與顛換兩種類型SNP比率（i/v）大致相當(dāng)（表2）。與親本池相比，兩個(gè)子代混池樣本中純合基因型比率明顯降低，分別為38.00%和38.82%，表明雜交增加了子代的雜合SNP率。

表1 ms7中雄性不育性狀的遺傳分析

表2 SNP數(shù)據(jù)統(tǒng)計(jì)表

進(jìn)一步對(duì)兩個(gè)極端性狀混池的SNP-index在染色體上的分布情況進(jìn)行統(tǒng)計(jì)，計(jì)算兩個(gè)子代混池的|ΔSNP-index|，結(jié)合性狀進(jìn)行關(guān)聯(lián)分析后發(fā)現(xiàn)，當(dāng)置信度為0.999時(shí)，第7染色體上15.5?16.7 Mb位置被篩選為候選區(qū)域（圖3）。

曲線是變異位點(diǎn)在染色體上所對(duì)應(yīng)的Index值；紅、綠和藍(lán)三根線分別是95%、99%和99.9%置信區(qū)間閾值線。

Fig. 3. Distribution of SNP-index on the chromosomes.

2.5 雄性不育基因ms7的定位

為了驗(yàn)證測(cè)序結(jié)果并進(jìn)行下一步的精細(xì)定位，我們根據(jù)兩個(gè)親本和IR96671的重測(cè)序結(jié)果在第7染色體長(zhǎng)臂上設(shè)計(jì)了62個(gè)InDel標(biāo)記，涵蓋了BSA分析99.9%閾值對(duì)應(yīng)的關(guān)聯(lián)區(qū)間，并向兩端延伸至區(qū)間7 968 748?21 483 477。對(duì)用于混池分析的材料按單株進(jìn)行基因分型驗(yàn)證，確定基因的候選區(qū)間為15.0 Mb?18.4 Mb (15 032 966? 18 400 602)。從中篩選8個(gè)親本間多態(tài)性明顯的InDel標(biāo)記對(duì)F2群體中剩余個(gè)體進(jìn)行基因型鑒定，并結(jié)合各單株的育性表型，進(jìn)一步縮小控制育性的定位區(qū)間（表3）。

表3?8個(gè)InDel標(biāo)記的引物信息

根據(jù)約14 000株F2群體的實(shí)驗(yàn)結(jié)果可挑選出7種類型重組單株，其基因型和表型數(shù)據(jù)如圖4所示。其中，R1-R7是不同的單株，S表示該株不育，F(xiàn)表示該株可育。通過這7個(gè)單株我們將界定于第7染色體上BSA11(16.97 Mb)與YD7045 (17.99 Mb)之間約1.17 Mb的范圍內(nèi)。為了驗(yàn)證定位的結(jié)果，在BSA11?YD7045之間挑選雜合型材料發(fā)展群體在實(shí)驗(yàn)室培養(yǎng)箱種植，經(jīng)標(biāo)記檢測(cè)，群體出現(xiàn)正常的基因型分離，從中挑選不同基因型的材料種植于溫室，其結(jié)果顯示，具有純合型的單株均表現(xiàn)為不育，而雜合型和父本型均可育，此結(jié)果與定位結(jié)果一致。

2.5?GO分析

水稻第7染色體BSA11和YD7045間約1.17 Mb的物理區(qū)間內(nèi)，根據(jù)日本晴參考基因組預(yù)測(cè)，共有137個(gè)基因位點(diǎn)。其中，已克隆的僅()[41]。剩下的136個(gè)位點(diǎn)中，根據(jù)基因本體(gene ontology, GO)注釋分析，轉(zhuǎn)座子或逆轉(zhuǎn)錄轉(zhuǎn)座子合計(jì)51個(gè)，其他預(yù)測(cè)具有生物學(xué)功能的位點(diǎn)42個(gè)。從分子功能看，能與蛋白結(jié)合(GO: 0005515)的9個(gè)，具有轉(zhuǎn)移酶活性(GO: 0016740)的8個(gè)，具有水解酶活性(GO: 0016787)的7個(gè)，等等；從生物學(xué)進(jìn)程看，和預(yù)測(cè)參與花發(fā)育進(jìn)程(GO:0009908)，和預(yù)測(cè)參與細(xì)胞死亡過程(GO: 0008219)，另有參與核酸類物質(zhì)代謝(GO: 0006139)的11個(gè)，參與脂質(zhì)代謝(GO: 0006629）的4個(gè)，參與碳水化合物代謝(GO: 0005975)的2個(gè)，參與底物運(yùn)輸(GO: 0006810)的4個(gè)。

3 討論

本研究利用組培獲得了水稻雄性不育突變體，除花粉敗育外其他表型與野生型沒有差異，遺傳分析表明該性狀受單個(gè)隱性核基因控制。水稻中已鑒定的影響育性的核基因中，如這樣造成花粉完全敗育，但又不影響胚囊或其他營養(yǎng)器官發(fā)育的基因，并不多見。

圖4 重組單株的基因型組成和育性表現(xiàn)

Fig 4. Genotypic composition and fertility performance of recombinant plants.

突變體中控制雄性不育的基因定位在第7染色體上的InDel標(biāo)記BSA11和YD7045之間，該區(qū)間內(nèi)預(yù)測(cè)有基因位點(diǎn)137個(gè)。根據(jù)RiceXPro表達(dá)譜數(shù)據(jù)庫（https://ricexpro.dna.affrc. go.jp/），共有42個(gè)基因在花藥中有表達(dá)，相對(duì)表達(dá)量最高的是，預(yù)測(cè)編碼絲氨酸羧肽酶，較高的有、、、和等，未發(fā)現(xiàn)花藥特異表達(dá)基因。這些位點(diǎn)中，僅/得到鑒定。據(jù)報(bào)道，OsMSH4能與OsMSH5互作，形成異源二聚體，調(diào)控減數(shù)分裂過程中交叉的形成。突變體從終變期開始異常，導(dǎo)致雌雄配子均沒有功能，花粉和胚囊均不育[42]。由于僅雄配子異常，與不同，推測(cè)它們不是同一位點(diǎn)突變。

本文對(duì)定位的物理區(qū)間1.17 Mb，仍然偏大，當(dāng)中遺傳因子眾多，那些預(yù)測(cè)有生物學(xué)功能的位點(diǎn)，以及其他尚未注釋的位點(diǎn)，都可能是的候選基因。為了確認(rèn)基因在中是否有變異，我們將中恢161和中的DNA序列分成13段，通過PCR擴(kuò)增后送公司測(cè)序，初步分析發(fā)現(xiàn)兩者無差異。此外，為了能快速定位到目標(biāo)位點(diǎn)，我們正對(duì)中恢161和進(jìn)行第三代測(cè)序，嘗試通過長(zhǎng)序列的拼裝，結(jié)合BSA二代測(cè)序結(jié)果，希望能直接找出與野生型在目標(biāo)DNA區(qū)間上的差異。

[1] 鄂志國, 程本義, 孫紅偉, 汪玉軍, 朱練峰, 林海, 王磊, 童漢華, 陳紅旗. 近40年我國水稻育成品種分析[J]. 中國水稻科學(xué), 2019, 33(6): 523-531.

E Z G, Cheng B, Sun H, Wang Y, Zhu L, Lin Hai, Wang L, Tong H, Chen H. Analysis on Chinese improved rice varieties in recent four decades[J]., 2019, 33(6): 523-531. (in Chinese with English abstract)

[2] 范優(yōu)榮, 曹曉風(fēng), 張啟發(fā). 光溫敏雄性不育水稻的研究進(jìn)展[J]. 科學(xué)通報(bào), 2016, 61(35): 3822-3832.

Fan Y R, Cao X F, Zhang Q F. Progress on photoperiod thermo-sensitive genic male sterile rice[J]., 2016, 61(35): 3822-3832. (in Chinese with English abstract)

[3] Fan Y, Yang J, Mathioni S M, Yu J, Shen J, Yang X, Wang L, Zhang Q, Cai Z, Xu C, Li X, Xiao J, Meyers B C, Zhang Q., producing phased small-interfering RNAs, regulates photoperiod-sensitive male sterility in rice [J]., 2016, 113(52): 15144.

[4] Zhou H, Liu Q, Li J, Jiang D, Zhou L, Wu P, Lu S, Li F, Zhu L, Liu Z, Chen L, Liu Y G, Zhuang C. Photoperiod- and thermo-sensitive genic male sterility in rice are caused by a point mutation in a novel noncoding RNA that produces a small RNA[J]., 2012, 22(4): 649-660.

[5] Ding J, Lu Q, Ouyang Y, Mao H, Zhang P, Yao J, Xu C, Li X, Xiao J, Zhang Q. A long noncoding RNA regulates photoperiod-sensitive male sterility, an essential component of hybrid rice[J]., 2012, 109(7): 2654-2659

[6] Zhou H, Zhou M, Yang Y, Li J, Zhu L, Jiang D, Dong J, Liu Q, Gu L, Zhou L, Feng M, Qin P, Hu X, Song C, Shi J, Song X, Ni E, Wu X, Deng Q, Liu Z, Chen M, Liu Y G, Cao X, Zhuang C. RNase ZS1processesUbmRNAs and controls thermosensitive genic male sterility in rice [J]., 2014, 5: 4884.

[7] Pitnjam K, Chakhonkaen S, Toojinda T, Muangprom A. Identification of a deletion inand development of gene-based markers for selection[J]., 2008, 228(5): 813-822.

[8] Yu B, Liu L, Wang T. Deficiency of very long chain alkanes biosynthesis causes humidity-sensitive male sterility via affecting pollen adhesion and hydration in rice [J]., 2019, 42(12): 3340.

[9] Xue Z, Xu X, Zhou Y, Wang X, Zhang Y, Liu D, Zhao B, Duan L, Qi X. Deficiency of a triterpene pathway results in humidity-sensitive genic male sterility in rice[J]., 2018, 9: 604.

[10] Chen H, Zhang Z, Ni E, Lin J, Peng G, Huang J, Zhu L, Deng L, Yang F, Luo Q, Sun W, Liu Z, Zhuang C, Liu YG, Zhou H. HMS1 interacts with HMS1I to regulate very-long-chain fatty acid biosynthesis and the humidity-sensitive genic male sterility in rice () [J]., 2020, 225(5): 2077-2093.

[11] 王多祥, 祝萬萬, 袁政, 張大兵. 水稻雄性發(fā)育功能基因的發(fā)掘及應(yīng)用[J]. 生命科學(xué), 2016, 28(10): 1180-1188.

Wang DX, Zhu WW, Yuan Z, Zhang DB. Functional research of rice male reproduction and its utilization in breeding[J]., 2016, 28(10): 1180-1188. (in Chinese with English abstract)

[12] Li L, Li Y, Song S, Deng H, Li N, Fu X, Chen G, Yuan L. An anther development F-box (ADF) protein regulated by tapetum degeneration retardation (TDR) controls rice anther development [J]., 2015, 241(1): 157-166.

[13] Lee S, Jung K H, An G, Chung Y Y. Isolation and characterization of a rice cysteine protease gene,, using T-DNA gene-trap system[J]., 2004, 54(5): 755-765.

[14] Niu N, Liang W, Yang X, Jin W, Wilson Z A, Hu J, Zhang D. EAT1 promotes tapetal cell death by regulating aspartic proteases during male reproductive development in rice [J]., 2013, 4 : 1445.

[15] Shi J, Tan H, Yu X H, Liu Y, Liang W, Ranathunge K, Franke R B, Schreiber L, Wang Y, Kai G, Shanklin J, Ma H, Zhang D.is required for anther and microspore development in rice and encodes a fatty acyl carrier protein reductase [J]., 2011, 23(6): 2225-2246.

[16] Mondol P C, Xu D, Duan L, Shi J, Wang C, Chen X, Chen M, Hu J, Liang W, Zhang D. Defective Pollen Wall 3 (DPW3), a novel alpha integrin-like protein, is required for pollen wall formation in rice [J]., 2020, 225(2): 807-822.

[17] Yang X, Liang W, Chen M, Zhang D, Zhao X, Shi J. Rice fatty acyl-CoA synthetase OsACOS12 is required for tapetum programmed cell death and male fertility [J]., 2017, 246(1): 105-122.

[18] Xu D, Shi J, Rautengarten C, Yang L, Qian X, Uzair M, Zhu L, Luo Q, An G, Wa?mann F, Schreiber L, Heazlewood J L, Scheller H V, Hu J, Zhang D, Liang W.() encodes an acyl transferase required for rice pollen development[J]., 2017, 173(1): 240-255.

[19] Men X, Shi J, Liang W, Zhang Q, Lian G, Quan S, Zhu L, Luo Z, Chen M, Zhang D. Glycerol-3-phosphate acyltransferase 3 (OsGPAT3) is required for anther development and male fertility in rice[J]., 2017, 68(3): 513-526.

[20] Yu J, Meng Z, Liang W, Behera S, Kudla J, Tucker MR, Luo Z, Chen M, Xu D, Zhao G, Wang J, Zhang S, Kim YJ, Zhang D. A rice Ca2+binding protein is required for tapetum function and pollen formation[J]., 2016, 172(3): 1772-1786.

[21] Zheng S, Li J, Ma L, Wang H, Zhou H, Ni E, Jiang D, Liu Z, Zhuang C.controls ROS production and the initiation of tapetal PCD by epigenetically regulatingexpression in rice anthers[J]., 2019, 116(15): 7549-7558.

[22] Aya K, Ueguchi-Tanaka M, Kondo M, Hamada K, Yano K, Nishimura M, Matsuoka M. Gibberellin modulates anther development in rice via the transcriptional regulation of GAMYB[J]., 2009, 21(5): 1453.

[23] Liu Z, Bao W, Liang W, Yin J, Zhang D. Identification ofand analysis of the regulatory role ofin rice anther development[J]., 2010, 52(7): 670-678.

[24] Xiang X J, Sun L P, Yu P, Yang Z F, Zhang P P, Zhang Y X, Wu W X, Chen D B, Zhan X D, Khan R M, Abbas A, Cheng S H, Cao L Y. The MYB transcription factor Baymax1 plays a critical role in rice male fertility[J]., 2021, 134(2): 453-471.

[25] Jung K H, Han M J, Lee Y S, Kim Y W, Hwang I, Kim M J, Kim Y K, Nahm B H, An G. Riceis a major regulator of early tapetum development[J]., 2005, 17(10): 2705-2722.

[26] Ko S S, Li M J, Ku M S B, Ho Y C, Lin Y J, Chuang M H, Hsing H X, Lien Y C, Yang H T, Chang H C, Chan M T. The bHLH142 transcription factor coordinates with TDR1 to modulate the expression ofand regulate pollen development in rice[J]., 2014, 26(6): 2486-2504.

[27] Fu Z, Yu J, Cheng X, Zong X, Xu J, Chen M, Li Z, Zhang D, Liang W. The rice basic Helix-Loop-Helix transcription factor TDR INTERACTING PROTEIN2 is a central switch in early anther development., 2014, 26(4): 1512-1524.

[28] Li N, Zhang D S, Liu H S, Yin C S, Li X X, Liang W Q, Yuan Z, Xu B, Chu H W, Wang J, Wen T Q, Huang H, Luo D, Ma H, Zhang D B. The ricegene is required for tapetum degradation and anther development[J]., 2006, 18(11): 2999-3014.

[29] Ji C, Li H, Chen L, Xie M, Wang F, Chen Y, Liu YG. A novel rice bHLH transcription factor, DTD, acts coordinately with TDR in controlling tapetum function and pollen development[J]., 2013, 6(5): 1715-1718.

[30] Cao H, Li X, Wang Z, Ding M, Sun Y, Dong F, Chen F, Liu L, Doughty J, Li Y, Liu Y X. Histone H2B monoubiquitination mediated by HISTONE MONOUBI-QUITINATION1 and HISTONE MONOUBIQUITI-NATION2 is involved in anther development by regulating tapetum degradation-related genes in rice[J]., 2015, 168(4): 1389.

[31] Li H, Yuan Z, Vizcay-Barrena G, Yang C, Liang W, Zong J, Wilson Z A, Zhang D.encodes a PHD-finger protein that is required for tapetal cell death and pollen development in rice[J]., 2011, 156(2): 615-630.

[32] Yang Z, Liu L, Sun L, Yu P, Zhang P, Abbas A, Xiang X, Wu W, Zhang Y, Cao L, Cheng S.functions as a transcriptional activator to regulate programmed tapetum development and pollen exine formation in rice[J]., 2019, 99(1-2): 175-191.

[33] Yang Z, Sun L, Zhang P, Zhang Y, Yu P, Liu L, Abbas A, Xiang X, Wu W, Zhan X, Cao L, Cheng S., encoding a PHD-finger transcription factor, regulates Ubisch bodies and pollen wall formation in rice[J]., 2019, 99(5): 844-861.

[34] Li X, Gao X, Wei Y, Deng L, Ouyang Y, Chen G, Li X, Zhang Q, Wu C. Rice APOPTOSIS INHIBITOR5 coupled with two DEAD-box adenosine 5’-triphosphate-dependent RNA helicases regulates tapetum degeneration [J]., 2011, 23(4): 1416-1434.

[35] Bai W, Wang P, Hong J, Kong W, Xiao Y, Yu X, Zheng H, You S, Lu J, Lei D, Wang C, Wang Q, Liu S, Liu X, Tian Y, Chen L, Jiang L, Zhao Z, Wu C, Wan J.() encodes an ATP-citrate lyase required for tapetum programmed cell death[J]., 2019, 181(3): 1223-1238.

[36] Zhang P, Zhang Y, Sun L, Sinumporn S, Yang Z, Sun B, Xuan D, Li Z, Yu P, Wu W, Wang K, Cao L, Cheng S. The rice AAA-ATPase OsFIGNL1 is essential for male meiosis[J]., 2017, 8: 1639.

[37] Tan H, Liang W, Hu J, Zhang D.encodes a secretory fasciclin glycoprotein required for male reproductive development in rice[J]., 2012, 22(6): 1127-1137.

[38] Li H, Durbin R. Fast and accurate long-read alignment with Burrows-Wheeler transform[J]., 2009, 25(14), 1754-1760.

[39] McKenna A, Hanna M, Banks E, Sivachenko A, Cibulskis K, Kernytsky A, Garimella K, Altshuler D, Gabriel S, Daly M, DePristo M A. The genome analysis toolkit: a MapReduce framework for analyzing next-generation DNA sequencing data[J]., 2010, 20(9): 1297-1303.

[40] Li H, Handsaker B, Wysoker A, Fennell T, Ruan J, Homer N, Marth G, Abecasis G, Durbin R. The sequence alignment/map format and SAMtools[J]., 2009, 25(16): 2078-2079.

[41] Zhang D, Luo X, Zhu L. Cytological analysis and genetic control of rice anther development[J]., 2011, 38(9): 379-390.

[42] Wang C, Wang Y, Cheng Z, Zhao Z, Chen J, Sheng P, Yu Y, Ma W, Duan E, Wu F, Liu L, Qin R, Zhang X, Guo X, Wang J, Jiang L, Wan J. The role of OsMSH4 in male and female gamete development in rice meiosis[J]., 2016, 67(5): 1447-1459.

Genetic Analysis and Gene Mapping of a Male Sterile Mutantin Rice

YANG Jinyu1,#, BAI Chen2,3,#, DING Xiaohui2,3, SHEN Hongfang1, WANG Lei1, YING Jiezheng1,*, E Zhiguo1,*

(1State Key Laboratory of Rice Biology, China National Rice Research Institute, Hangzhou 310006, China;2Key Laboratory of Plant Molecular Physiology, Institute of Botany, Chinese Academy of Sciences, Beijing 100093, China;3University of Chinese Academy of Sciences, Beijing 100049, China;#These authors contributed equally to the work;*Corresponding author, E-mail: yingjiezheng@caas.cn; ezhiguo@caas.cn)

【】Genes related to pollen development and fertility can be identified through the study on male sterile mutants, which is helpful to analyze the whole regulatory network of male reproductive development in rice. 【】The fertility and agronomic traits of the() mutant and its wild type were compared under conventional planting conditions in Fuyang, Zhejiang and Lingshui, Hainan. The sterile gene was mapped by bulked segregant analysis and map-based cloning.【】During the whole growth period, the growth rate ofmutant was the same as that of wild type, and there was no significant difference in plant height, tiller number, leaf number, leaf size, panicle length and spikelets per panicle betweenand its wild type. However,exhibited slender and white anthers with inactive pollen grains which couldn’t be stained with I2-KI, and no seeds were produced. Observation results of anther cross-sections exhibited that tapetum programmed cell death (PCD) was delayed. Genetic analysis showed that pollen abortion was controlled by a single recessive nuclear gene, which was mapped to a 1.17 Mb region between BSA11 and YD7045 on chromosome 7. 【】It will contribute to the further molecular cloning and functional analysis of male sterile genein rice.

rice; male sterile; bulked segregant analysis; gene mapping

10.16819/j.1001-7216.2022.210302

2021-03-03；

2021-04-19。

浙江省自然科學(xué)基金資助項(xiàng)目(LY17G030031); 國家自然科學(xué)基金面上項(xiàng)目(32072050, 71773140); 中央級(jí)公益性科研院所基本科研業(yè)務(wù)費(fèi)專項(xiàng)(CPSIBRF-CNRRI-202129)。