花椰菜溫敏雄性不育系GS-19花藥敗育的細胞學及轉(zhuǎn)錄組分析

陶興林，侯棟，朱惠霞，劉明霞，張金文，胡立敏

?

花椰菜溫敏雄性不育系GS-19花藥敗育的細胞學及轉(zhuǎn)錄組分析

陶興林1,2，侯棟1，朱惠霞1，劉明霞1，張金文2，胡立敏1

（1甘肅省農(nóng)業(yè)科學院蔬菜研究所/農(nóng)業(yè)部園藝作物生物學與種質(zhì)創(chuàng)制西北地區(qū)科學觀測實驗站，蘭州 730070；2甘肅農(nóng)業(yè)大學/甘肅省作物遺傳改良與種質(zhì)創(chuàng)新重點實驗室，蘭州 730070）

【目的】對花椰菜溫敏雄性不育系的花器形態(tài)、花藥敗育時期和敗育分子機理進行研究，明確其不育特性發(fā)生的細胞學和分子機理，為進一步研究其雄性不育奠定理論基礎?！痉椒ā恳曰ㄒ藴孛粜坌圆挥礕S-19為試驗材料，不同溫度（15℃和20℃）處理，采用形態(tài)學、細胞學方法及高通量測序技術（Illumina HiSeq 2500），研究其形態(tài)特征、花藥敗育的細胞學及分子機理。【結(jié)果】花椰菜溫敏雄性不育系GS-19的育性轉(zhuǎn)換受溫度控制，高溫（20℃）不育，低溫（15℃）育性恢復。GS-19不育株與可育株花器差異顯著，不育株花器顯著小于可育株。花藥細胞學觀察發(fā)現(xiàn)GS-19的花藥發(fā)育過程有花粉母細胞的分化，形成正?；ǚ勰遥划a(chǎn)生花粉?；蛘弋a(chǎn)生微量的無生活力的花粉粒，花藥發(fā)育受阻于花粉母細胞到四分體時期，形成了花粉粒外壁發(fā)育異常的擬“四分體孢子”，逐漸降解，剩下花粉空殼，屬于花粉母細胞敗育類型。GS-19不育株和可育株花蕾轉(zhuǎn)錄組分析共獲得67 930個Unigene，其中Nt數(shù)據(jù)庫比對到相似序列數(shù)最多（52 191個），Nr數(shù)據(jù)庫比對到相似序列次之（46 447個），KOG數(shù)據(jù)庫比對相似序列最少（13 198個）；GO注釋到25 336個Unigene；KOG注釋到13 198個Unigene；KEGG注釋到14 778個Unigene?；虿町惐磉_分析發(fā)現(xiàn)GS-19不育株花蕾和可育株花蕾中有2 170個基因差異表達，1 078個基因上調(diào)表達和1 092個基因下調(diào)表達。【結(jié)論】花椰菜溫敏雄性不育系GS-19為高溫不育類型，不育株花器顯著小于可育株，花絲顯著縮短，花藥萎縮、干癟，花藥發(fā)育受阻于花粉母細胞到四分體時期，屬于花粉母細胞敗育類型。轉(zhuǎn)錄組測序獲得了67 930個Unigene，差異表達基因2 170個。

花椰菜；溫敏雄性不育系；花藥敗育；細胞學；轉(zhuǎn)錄組分析

0 引言

【研究意義】溫敏雄性不育系GS-19作為花椰菜育種中一種新型的雄性不育系，在兩系法花椰菜雜交育種中有十分重要的作用。明確花椰菜溫敏雄性不育系GS-19育性轉(zhuǎn)換過程中的花器形態(tài)特征、花藥發(fā)育的敗育時期、特點和基因表達差異，不僅為闡明該花椰菜溫敏雄性不育系的遺傳及分子機理奠定基礎，而且為花椰菜溫敏不育系的利用提供理論依據(jù)?！厩叭搜芯窟M展】花椰菜（）是一種非常重要的十字花科蔬菜作物，被公認為是最有營養(yǎng)的作物之一，特別是鈣、抗氧化劑、維生素A、維生素K、胡蘿卜素、核黃素及鐵的含量很豐富[1-3]，是雜種優(yōu)勢利用潛力很大的一類作物，而雄性不育途徑是培育花椰菜新品種的主要途徑之一。雄性不育主要特征是雄蕊發(fā)育受阻不能產(chǎn)生正常功能的花粉，是雜種優(yōu)勢利用中的重要材料和工具。首個水稻光敏核不育系農(nóng)墾58S的發(fā)現(xiàn)促使兩系法雜交水稻育種取得突破性進展[4]。自此以后，小麥[5]、玉米[6]、大豆[7]、水稻[8]、油菜[9]等作物中都發(fā)現(xiàn)了光（溫）敏雄性不育材料，許多學者對其敗育時期和方式進行了研究，發(fā)現(xiàn)不同作物或者同一作物不同的雄性不育系之間也有差異。對水稻[10-11]、小麥[12]、油菜[13-15]、玉米[16]、茄子[17]、辣椒[18]、甘藍[19]等作物的雄性不育系細胞學研究較多，主要有無花粉囊型、花粉母細胞敗育型和單核敗育型3種。高通量轉(zhuǎn)錄組測序技術以其速度快、測序長、通量高、準確率高等優(yōu)點已被廣泛應用于生物學[20]、醫(yī)學[21]、農(nóng)學[22-24]等研究領域。在雄性不育相關研究中，高通量測序技術應用主要集中在水稻、小麥、玉米、棉花、菊花等20多個物種中，包括細胞質(zhì)雄性不育[25-26]、細胞核雄性不育[27-28]、化學誘導雄性不育[29]、光溫敏雄性不育[30-31]等多種敗育類型，這對于雄性不育敗育機理的深入研究具有重要意義，也為相關研究提供了豐富的數(shù)據(jù)信息。【本研究切入點】2005年，筆者課題組從花椰菜種質(zhì)資源圃里發(fā)現(xiàn)了幾株溫度敏感的不育株，后經(jīng)多代自交選育成育性穩(wěn)定的溫敏不育系GS-19[32]，屬高溫不育類型，該雄性不育系已經(jīng)在生產(chǎn)中進行了應用[33]，但導致其花藥敗育的細胞學和轉(zhuǎn)錄組水平上的研究尚未見報道。【擬解決的關鍵問題】以花椰菜溫敏雄性不育系GS-19為試驗材料，不同溫度（15℃和20℃）處理，采用形態(tài)學、細胞學方法及高通量測序技術，研究花椰菜溫敏雄性不育系花器形態(tài)特征、花藥敗育的細胞學特點及轉(zhuǎn)錄組水平上的基因差異，明確不育系GS-19的花器敗育特征、花藥敗育時期、細胞學特征及差異基因表達特點，為花椰菜溫敏雄性不育進一步研究和應用奠定基礎。

1 材料與方法

試驗于2012—2016年在甘肅省農(nóng)業(yè)科學院蔬菜研究所農(nóng)業(yè)部園藝作物生物學與種質(zhì)創(chuàng)制西北地區(qū)科學觀測實驗站進行。

1.1 試驗材料

試驗選取的花椰菜溫敏雄性不育系GS-19由甘肅省農(nóng)業(yè)科學院蔬菜研究所選育和提供。育性表現(xiàn)與溫度有關，溫度高于18℃表現(xiàn)為不育，溫度低于16℃時表現(xiàn)為育性完全恢復。

1.2 育性轉(zhuǎn)換觀察

2014年3月，在智能玻璃溫室中花盆直播不育系GS-19，共播種25盆直徑為38 cm的花盆，顯球后分別移入20℃（16 h光照/8 h黑暗）和15℃（16 h光照/8 h黑暗）的人工智能氣候箱中，開花期觀察統(tǒng)計育性情況。

1.3 花器特征觀察

在盛花期，采集在不同溫度處理下不育系GS-19開花前3—4 d的花蕾和當天完全開放的花各20朵，用游標卡尺測量7個指標，即花蕾長、花蕾寬、花冠開度、花柱長、花藥長、花絲長和花藥總長各部分大小，進行比較方差分析。

1.4 掃描電鏡觀察

取不育系GS-19的不育和可育株不同大小花蕾（0.25、0.5、0.75和1.0 cm），用手術刀橫切（或剝出花藥），確保切面整齊，將導電膠帶粘在樣品托盤上，再將試驗樣品切面向上按已標記的順序粘在膠帶上，并用鑷子進行固定，粘貼面必須保持平整，確保樣品粘貼牢固，直接放入樣品臺，用HitachiS-3400N型掃描電鏡觀察并照相。

1.5 透射電鏡觀察

選取不育系GS-19不育和可育株不同發(fā)育時期（孢原細胞、花粉母細胞、單核花粉粒和花粉粒成熟期）的花蕾，快速剝開花蕾，取出花藥。采用戊二醛和四氧化鋨雙重固定法，先用2.5%戊二醛固定6 h，再用1%四氧化鋨固定2 h，丙酮進行梯度（70%、80%、90%、100%）脫水，環(huán)氧樹脂包埋，超薄切片機切成厚度為5 μm的切片，載玻片上粘片，用乙酸雙氧鈾和檸檬酸鉛染色，TEM-1230透射電鏡觀察并照相。

1.6 轉(zhuǎn)錄組分析

RNA提?。涸诓挥礕S-19的盛花期，分別取10株待測不育株和可育株開花前發(fā)育成熟的花蕾，采用Trizol法分別提取不育和可育株花蕾總RNA，每個樣本取等量混合組成1個RNA池，用帶有Oligo(dT)的磁珠富集mRNA，將得到的mRNA逆轉(zhuǎn)錄成dsDNA，對其末端修復，加poly(A)并連接測序接頭，制備測序文庫，PCR富集測序樣本。

功能注釋：通過Illumina HiSeq 2500 的測序平臺，采用PE125的測序方法，構(gòu)建GS-19不育和可育株花蕾轉(zhuǎn)錄組文庫后進行測序并獲取相關數(shù)據(jù)。通過Blastx比對工具，將不育系GS-19與蛋白數(shù)據(jù)庫進行比對，根據(jù)基因的相似性進行功能注釋，得到與給定Unigene具有最高序列相似性的蛋白及功能注釋信息。數(shù)據(jù)庫包括Nr（Non-redundant protein database，非冗余蛋白數(shù)據(jù)庫）、Nt（NCBI non-redundant nucleotide database，非冗余核苷酸數(shù)據(jù)庫）；SwissProt（SwissProt protein database，蛋白質(zhì)序列數(shù)據(jù)庫）、COG（Cluster of Orthologous Groups，蛋白質(zhì)直系同源數(shù)據(jù)庫）、Pfam（Protein families database，蛋白質(zhì)家族域數(shù)據(jù)庫）、GO（Gene Ontology，基因本體論數(shù)據(jù)庫）、KEGG（Kyoto Encyclopedia of Genes and Genomes，東京基因與基金組百科全書）。

差異基因分析：采用TMM對read count數(shù)據(jù)進行標準化處理，再用DEG seq進行差異分析，篩選閾值為qvalue＜0.005 & |log2FoldChange|＞1；對于差異基因，如果基因的log2Foldchange＞0，則認為該差異基因是上調(diào)；反之，若log2Foldchange＜0，認為該差異基因是下調(diào)。

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

試驗結(jié)果為3次重復的平均值，用Excel 2010進行數(shù)據(jù)整理、標準偏差計算及圖表的繪制，差異顯著性用SPSS 18.0 統(tǒng)計軟件進行分析。

2 結(jié)果

2.1 溫敏雄性不育系GS-19的育性轉(zhuǎn)換分析

育性統(tǒng)計結(jié)果發(fā)現(xiàn)，2個溫度處理后的溫敏不育系GS-19的育性變化存在顯著差異。溫敏雄性不育系GS-19在2個溫度條件下可育率發(fā)生了顯著的變化，在15℃條件下，完全可育，可育率為100%，20℃條件下完全不育，可育率為0，說明溫敏不育系GS-19育性轉(zhuǎn)換受溫度控制，屬于高溫不育類型。

2.2 雄性不育系GS-19的花器特征分析

花椰菜溫敏雄性不育系GS-19的不育株和可育株的花器有顯著差異，不育株花器明顯小于可育株。方差分析表明（圖1），除雄蕊數(shù)量沒有變化外，不育和可育株的花蕾長、花蕾寬、花冠開度、花絲長、花藥總長（花絲+花藥）及花柱長都存在顯著差異，不育株花明顯小于可育株，且雄蕊很短，在開花期看不到雄蕊，雄蕊萎縮在基部，只有柱頭明顯外露（圖2），而且雄蕊也不產(chǎn)生花粉，但可育株的雄蕊基本與雌蕊高度一致，產(chǎn)生大量花粉，能完成正常的授粉。

2.3 雄性不育系GS-19的電鏡掃描觀察

2.3.1 不育和可育株花藥外部形態(tài)觀察 通過對不育系GS-19不育和可育的花藥進行電鏡掃描觀察，不育和可育株花藥外部形態(tài)存在顯著差異（圖3）?；ɡ匍L度為0.25 cm時，不育和可育株花藥外部形態(tài)無差異，棒狀形，光滑飽滿；花蕾長度為0.5 cm時，不育和可育株花藥變成了錐形，不育株的花藥壁出現(xiàn)輕微萎縮現(xiàn)象，表面變得不光滑，可育株花藥發(fā)育正常，仍然飽滿光滑；花蕾長度為0.75 cm時，不育株花藥變成干癟狀，明顯萎縮，表面變得粗糙不平整，中空明顯，表現(xiàn)為完全不育，可育株的花藥增大明顯，發(fā)育正常，飽滿光滑；開花期（花蕾長度為1.0 cm）不育株花藥發(fā)育只剩下開裂的空殼花藥壁，無花粉粒，可育株的花藥開裂，裂口處能看到大量飽滿的花粉粒。

2.3.2 不育和可育株花藥橫切面觀察 通過對不育系GS-19不育和可育株花蕾橫切面電鏡掃描觀察，不育和可育株花蕾發(fā)育橫切面存在顯著差異（圖4）?；ɡ匍L度為0.25 cm時，不育和可育株花藥橫切面差異不明顯，都被早期發(fā)育的小孢子填充；花蕾長度為0.5 cm時，不育株花藥發(fā)生了敗育，花藥橫切面有變空現(xiàn)象，能清楚地看到花藥壁，但可育株花藥內(nèi)部仍然被小孢子填充，沒有出現(xiàn)變空現(xiàn)象；花蕾長度為0.75 cm時，不育株花藥橫切面中空明顯，花藥壁變薄，完全不育，而可育株花藥橫切面變化不明顯，明顯看到花藥里被發(fā)育的小孢子填充，沒有中空現(xiàn)象；開花期（花蕾長度為1.0 cm左右）不育株花藥橫切電鏡掃描發(fā)現(xiàn)花藥就剩下薄薄的一層花藥壁，沒有花粉粒，但可育株的花藥被成熟花粉粒填充。

2.3.3 不育和可育株花粉粒的觀察 開花期取不育和可育株進行花藥掃描電鏡觀察（圖5），發(fā)現(xiàn)不育株（20℃）花藥干癟變形，從表面看不到花粉粒，但是剝開花藥，發(fā)現(xiàn)極少數(shù)敗育的花粉粒，敗育花粉粒為不規(guī)則形，凹陷明顯。可育株（5℃）花藥飽滿，明顯看到大量的花粉粒，花粉粒呈橢圓形，表面光滑，飽滿。

2.4 GS-19的透射電鏡觀察

通過透射電鏡觀察，發(fā)現(xiàn)不育系GS-19的花藥發(fā)育受阻于花粉母細胞到四分體時期，屬于花粉母細胞敗育類型?；ㄋ幇l(fā)育早期，不育系GS-19不育和可育株在細胞超微結(jié)構(gòu)上沒有差異，花藥發(fā)育完全正常，分化形成了表皮、藥室內(nèi)壁、中層等結(jié)構(gòu)，絨氈層和造孢細胞之間沒有明顯的分界線，形成不規(guī)則多邊形的造孢細胞，能看到細胞核，細胞質(zhì)濃厚，看不到明顯的液泡（圖6-A、E）。隨著造孢細胞分化，形成了花粉母細胞，此時可育株的花藥能明顯看到一圈排列整齊的長方形絨氈層，絨氈層細胞核較大，核仁染色深?；ǚ勰讣毎唤q氈層包裹，其細胞大，形狀不規(guī)則，核大且染色較深（圖6-B）。不育株花藥也能形成絨氈層和花粉母細胞，部分花粉母細胞沒有細胞核（圖6-F）。小孢子分化形成了被著色深的絨氈層包裹的四分體孢子，可育株的四分體孢子外壁能明顯看到深色密點狀包裹帶，其上出現(xiàn)大量的發(fā)狀物質(zhì)。四分體孢子中含有大量的細胞器，細胞內(nèi)含物逐漸增多至布滿整個細胞，染色深且均勻（圖6-C）。不育株只能形成擬“四分體孢子”，和可育株一樣，外側(cè)同樣出現(xiàn)了深色密點狀包裹帶，包裹帶上出現(xiàn)了大量的深色發(fā)狀物質(zhì)，但是四分體孢子的細胞器顯著減少，內(nèi)含物降解，出現(xiàn)液泡吞噬現(xiàn)象，導致花藥敗育（圖6-G）。四分體孢子繼續(xù)發(fā)育，可育株能形成成熟花粉粒，花粉粒近圓形，內(nèi)壁物質(zhì)顯著增厚，形成了萌發(fā)孔，內(nèi)含物豐富且染色深，絨氈層細胞完全降解消失（圖6-D）。但是不育株的擬“四分體孢子”不能正常發(fā)育，液泡解體，只剩下壁外側(cè)有發(fā)狀電子致密帶，出現(xiàn)沒有活力的“空腔”花粉?，F(xiàn)象（圖6-H）。

2.5 轉(zhuǎn)錄組分析

2.5.1 Unigene 的功能注釋相關結(jié)果分析 通過BLAST程序，對所組裝獲得不育系GS-19的Unigene進行7大數(shù)據(jù)庫的比對（表1）。結(jié)果表明，功能數(shù)據(jù)庫共注釋到67 930個Unigene，7大數(shù)據(jù)庫比對到相似序列數(shù)目存在明顯差異，比對到的相似序列數(shù)從13 198—52 191個，其中Nt數(shù)據(jù)庫比對到相似序列數(shù)最多，為52 191個，Nr數(shù)據(jù)庫比對到相似序列次之，為46 447個，KOG數(shù)據(jù)庫比對相似序列最少，為13 198個。Nr功能注釋到的匹配物種中，油菜Unigene最多（33 108），然后依次為蕪菁（8 841）、擬南芥（643）、芥藍（634）、山崳菜（626）等。從匹配物種中發(fā)現(xiàn)，不育系GS-19屬于十字花科作物，匹配率最高的也是十字花科的油菜和蕪菁，從親緣關系遠近也可以證明這一點。

表1 花椰菜溫敏雄性不育系GS-19的Unigenes功能注釋

GO注釋到25 336個Unigene（表2），分為細胞組分、分子功能、生物學過程3個本體共55個功能組；KOG注釋到13 198個Unigene（表3），分為26類，注釋到最多的是基因功能預測，共2 264個，其次為翻譯后修飾，蛋白質(zhì)折疊，分子伴侶（1 835個），最少的是未知蛋白（1個）；KEGG注釋到14 778個Unigene（表4），分細胞過程、環(huán)境信息處理、遺傳信息處理、代謝及有機系統(tǒng)5大類共19個代謝途徑，其中涉及較多的有碳水化合物（1 350個），翻譯（1 275個），折疊、分類和降解（1 016個）等，較少的為多糖生物合成與代謝（161個），膜輸送、膜轉(zhuǎn)運（60個）。

2.5.2 差異表達分析 本試驗中發(fā)現(xiàn)不育系GS-19不育株的上調(diào)表達基因表達水平高于可育株的基因表達水平。如圖7所示，在TMS（GS-19不育株花蕾）和TMF（GS-19可育株花蕾）之間發(fā)現(xiàn)的67 931 Unigenes中有2 170個基因差異表達，其中1 078個基因上調(diào)表達和1 092個基因下調(diào)表達。每個比較組合都會得到一個差異基因集，將所有比較組合的差異基因集的并集在每個試驗組/樣品中的的FPKM（Fragment Per Kilobase of exon model per Million mapped reads）值，就得到了TMS和TMF的聚類分析結(jié)果（圖8）。

表2 花椰菜溫敏雄性不育系GS-19 Unigene的GO分類

表3 花椰菜溫敏雄性不育系GS-19 Unigene的KOG系統(tǒng)分類

3 討論

本研究發(fā)現(xiàn)花椰菜溫敏雄性不育系GS-19的育性受溫度控制，高溫（20℃）條件下表現(xiàn)為不育，低溫（15℃）條件下完全可育，與前人在甘藍型油菜溫敏核不育系SP2S[14]、TE5A[15]、細胞質(zhì)雄性不育系392A[13]的育性轉(zhuǎn)換研究結(jié)果相同，都表現(xiàn)為高溫不育、低溫育性恢復。但在油菜溫敏雄性不育系Huiyou50S的育性轉(zhuǎn)換研究結(jié)果正好與本研究結(jié)果相反，表現(xiàn)出低溫不育、高溫可育的現(xiàn)象[34]。為進一步對溫敏雄性不育系GS-19利用，在制種過程中要找到適宜區(qū)域，注意監(jiān)控溫度條件，避免影響雜種種子的純度。

花器形態(tài)特征是雄性不育和可育區(qū)別的直接表型，大部分研究表明雄性不育的花明顯小于可育株的花，不育花器形態(tài)有花冠退化型、雄蕊萎縮型、花粉敗育型3種[35-36]。本研究發(fā)現(xiàn)，高溫條件下溫敏不育系GS-19的花器顯著小于低溫下溫敏核不育系GS-19，尤其是花的發(fā)育后期，不育株花器的雄蕊都萎縮于雌蕊基部，柱頭外露，花藥干癟，與前人研究中的雄蕊萎縮型表型一致。童旭宏等[37]觀察了陸地棉的G22A、H109A和1793A 3個不育系，均表現(xiàn)花瓣小、花絲短，柱頭明顯外露，敗育徹底，無花粉；張體德等[38]發(fā)現(xiàn)芝麻細胞核雄性不育系ms86-1的花器也是不育株小于可育株，這些研究結(jié)果都與本研究結(jié)果相同。

隨著花椰菜雄性不育系的發(fā)現(xiàn)，花藥發(fā)育的細胞學研究也顯得非常重要。目前，雄性不育敗育研究主要集中在油菜、小麥、水稻等大田作物中，主要有無花粉囊型、花粉母細胞敗育型和單核敗育型3種，已經(jīng)研究發(fā)現(xiàn)的典型無花囊型油菜不育系有533S[39]、宜3A[40]、212A[41]、Shaan-GMS[42]、9012[43]、BnaC.Tic40[44]、BnMs3[45]等。本研究發(fā)現(xiàn)在花蕾發(fā)育早期（≤1 mm）時，低溫條件下的溫敏雄性不育系GS-19處于孢原時期和造孢時期，能分化形成造孢細胞和壁細胞，該時期高溫條件下的GS-19的花藥發(fā)育也正常。隨著花蕾的繼續(xù)發(fā)育（＞1 mm），高溫條件下溫敏雄性不育系GS-19不育株的花粉母細胞和可育株一樣，也能進行正常減數(shù)分裂，但形成的四分體孢子不正常，出現(xiàn)擬四分體孢子現(xiàn)象，細胞器顯著減少，內(nèi)含物降解，出現(xiàn)液泡吞噬現(xiàn)象，導致敗育，屬于花粉母細胞敗育類型，與溫敏核不育系小麥BNS[12]和BNS366[46]、甘藍Ogura CMS[47]的敗育類型一致，此外，水稻[45]、小麥[48]、高粱[49]、大豆[50]、紅菜薹[51]等作物的雄性不育系也是因為減數(shù)分裂過程出現(xiàn)問題，導致花粉母細胞不能形成四分體，也與本試驗不育類型一致。但是一些其他的溫敏雄性不育水稻、小麥、甘藍型油菜的花藥發(fā)育的細胞學觀察發(fā)現(xiàn)，這些敗育類型屬于單核敗育型[10,34,52]。推斷GS-19是一個具有重要利用價值的花椰菜溫敏雄性不育系，其不育發(fā)生的細胞學研究結(jié)果，對于揭示不育發(fā)生的細胞學機制和雜交制種具有重要的指導意義。

表4 花椰菜溫敏雄性不育系GS-19Unigene的KEGG代謝途徑分析

通過基因功能注釋，可以明確基因的功能，也可以了解基因的保守結(jié)構(gòu)域，基因自身表達條件和潛在的代謝途徑[53]。本研究通過GO數(shù)據(jù)庫，對溫敏雄性不育系GS-19 Unigene進行基因生物學特征功能分類，發(fā)現(xiàn)37 220個GO條目歸屬于細胞組分中的20個功能組；29 992個GO條目歸屬于分子功能中的10個功能組，62 204個GO條目歸屬于生物學過程中的25個功能組，明顯看出GO分析雄性不育相關轉(zhuǎn)錄組數(shù)據(jù)量大，相關差異表達基因集聚在代謝中的多個通路，初步證明花椰菜溫敏雄性不育產(chǎn)生機理的復雜性。前人研究發(fā)現(xiàn)植物代謝過程中的花粉壁形成、類黃酮合成等過程都與植物雄性不育形成過程中的花粉發(fā)育密切相關[54]，細胞程序性死亡也與植物雄性不育關系密切[55]，這些研究結(jié)果與本研究結(jié)果相似。相反，大豆的細胞質(zhì)雄性不育研究中，發(fā)現(xiàn)參與胚發(fā)育過程的基因數(shù)量最多[56]；大蒜細胞質(zhì)雄性不育的花青素代謝相關基因在差異表達基因中數(shù)量較多[57]，說明敗育方式差異會導致雄性不育表現(xiàn)出特有差異基因富集的功能通路。

基因組和轉(zhuǎn)錄組數(shù)據(jù)的缺乏大大限制花椰菜溫敏雄性不育系的分子機理方面的相關研究。目前，對溫敏雄性不育的轉(zhuǎn)錄組分析研究發(fā)現(xiàn)，遺傳背景、敗育時期和方式等多方面因素都會影響基因的表達水平，Mei等[58]在蘿卜2個核背景存在差異的同一細胞質(zhì)雄性不育材料的轉(zhuǎn)錄組分析中，發(fā)現(xiàn)與其相對應的保持系差異表達基因數(shù)分別為3 843和2 041個，這些基因都是與不育特性密切相關的基因；An等[59]發(fā)現(xiàn)甘藍型油菜細胞質(zhì)雄性不育系Pol不育株中的基因上調(diào)表達數(shù)目顯著高于可育株。本研究發(fā)現(xiàn)，在不育株花蕾和可育株花蕾中有2 170個基因差異表達，其中1 078個上調(diào)表達和1 092個下調(diào)表達，下調(diào)表達主要出現(xiàn)在不育材料中，上調(diào)表達主要出現(xiàn)在可育材料中，與Qu等[60]在油菜細胞核雄性不育系中的研究結(jié)果相一致，可育材料中的差異基因以上調(diào)表達為主。但是與Mei等[58]和An等[59]的研究結(jié)果有一定差異，可能與基因型有很大的關系，這些研究都說明了基因的表達量及表達的上、下調(diào)關系與遺傳背景和差異基因的功能密切相關，因此，可以通過對基因表達水平變化的分析，來預測和研究其基因功能的變化。

4 結(jié)論

花椰菜溫敏雄性不育系GS-19育性轉(zhuǎn)換受溫度控制。GS-19不育與可育株花器差異顯著，不育株的花器明顯小于可育株。超微結(jié)構(gòu)觀察發(fā)現(xiàn)GS-19不育株的花藥發(fā)育過程有造胞細胞和花粉母細胞的分化，形成正?；ǚ勰?，不產(chǎn)生花粉粒或者產(chǎn)生微量的無生活力的花粉粒，花藥發(fā)育受阻于花粉母細胞到四分體時期，形成了花粉粒外壁發(fā)育異常的擬“四分體孢子”，隨著小孢子發(fā)育，擬“四分體孢子”逐漸降解，只剩下花粉空殼，屬于花粉母細胞敗育類型。轉(zhuǎn)錄組分析發(fā)現(xiàn)，GS-19共獲得了67 930個Unigene，GO注釋到25 336個Unigene；KOG注釋到13 198個Unigene；KEGG注釋到14 778個Unigene。GS-19不育株花蕾和可育株花蕾中有2 170個基因差異表達，1 078個上調(diào)表達和1 092個下調(diào)表達。

References

[1] VALLEJO F, GARCIA-VIGUERA C, TOMAS-BARBERAN F A. Changes in broccoli () health promoting compounds with inflorescence development., 2003, 51: 3776-3782.

[2] ABDEL-WAHHAB M A, ALY S E. Antioxidants and radical scavenging properties of vegetable extracts in rats fed aflatoxin- contaminated diet., 2003, 51: 2409-2414.

[3] 李杰, 賈豪語, 頡建明, 郁繼華, 楊萍. 生物肥部分替代化肥對花椰菜產(chǎn)量、品質(zhì)、光合特性及肥料利用率的影響. 草業(yè)學報, 2015, 24(1): 47-55.

LI J, JIA H Y, XIE J M, YU J H, YANG P. Effects of partial substitution of mineral fertilizer by bio-fertilizer on yield, quality, photosynthesis and fertilizer utilization rate in broccoli., 2015, 24(1): 47-55. (in Chinese)

[4] 石明松. 對光照長度敏感的隱性雄性不育水稻的發(fā)現(xiàn)與初步研究.中國農(nóng)業(yè)科學, 1985, 18(2): 44-48.

SHI M S. The discovery and study of the photo-sensitive recessive male-sterile rice., 1985, 18(2): 44-48.(in Chinese)

[5] GUO R X, SUN D F, TAN Z B, RONG D F, LI C D. Two recessive genes controlling thermo-photoperiod-sensitive male sterility in wheat., 2006, 112(7): 1271-1276.

[6] TANG J H, FU Z Y, HU Y M, LI J S, SUN L L, JI H Q. Genetic analyses and mapping of a new thermo-sensitive genic male sterile gene in maize., 2006, 113: 11-15.

[7] FRASCH R M, WEIGAND C, PEREZ P T, PALMER R G, SANDHU D. Molecular mapping of 2 environmentally sensitive male sterile mutants in soybean., 2010, 102: 11-16.

[8] DING J H, LU Q, YANG Y D, MAO H L, ZHANG P B, YAO J L, XU C G, LI X H, XIAO J H, ZHANG Q F. A long noncoding RNA regulates photoperiod-sensitive male sterility, an essential component of hybrid rice., 2012, 109(7): 2654-2659.

[9] ZENG X H, LI W P, WU Y H, LIU F, LUO J L, CAO Y L, ZHU L, LI Y J, LI J, YOU Q B, WU G. Fine mapping of a dominant thermo-sensitive genic male sterility gene (BntsMs) in rapeseed () with AFLP-and-derived PCR markers., 2014, 127: 1733-1740.

[10] SHI Y, ZHAO S, YAO J. Premature tapetum degeneration:a major cause of abortive pollen development in photo-period sensitive genic male sterility in rice., 2009, 51(8): 774-781.

[11] ZHOU H, LIU Q, LI J, JIANG D G, ZHOU L Y, WU P, LU S, ZHU L Y, LIU Z L, CHEN L T, LIU G Y, ZHUANG C X. Photoperiod- and thermo-sensitive genic male sterility in rice are caused by a point mutation in a novel non-coding RNA that produces a small RNA．, 2012, 22(4): 649-660.

[12] 李東霄, 鄧小莉, 李淦, 徐龍龍, 茹振鋼. 溫敏核不育小麥可育花藥和敗育花藥發(fā)育觀察. 中國細胞生物學學報, 2012, 34(9): 880-885.

LI D X, DENG X L, LI G, XU L L, RU Z G. Observation in fertile and sterile anthers of a thermo-sensitive genic male-sterile wheat (L)., 2012, 34(9): 880-885. (in Chinese)

[13] 聶明建, 王國槐, 朱衛(wèi)平. 甘藍型油菜 3 種類型雄性不育系花藥敗育的細胞學研究. 中國農(nóng)業(yè)科學, 2007, 40(7):1543-1549.

NIE M J, WANG G H, ZHU W P. Cytology research on the anther abortion of three male sterility lines in rapeseed (L.)., 2007, 40(7): 1543-1549. (in Chinese)

[14] YU C, GUO Y, GE J, HU Y M, DONG J G, DONG Z S. Characterization of a new temperature-sensitive male sterile line SP2S in rapeseed (L.)., 2015, 206(2): 473-485.

[15] 李可琪, 曾新華, 袁榮, 閆曉紅, 吳剛. 甘藍型油菜溫敏細胞核雄性不育系 TE5A花藥發(fā)育的細胞學研究. 中國農(nóng)業(yè)科學, 2016, 49(12): 2408-2417.

LI K Q, ZENG X H, YUAN R, YAN X H, WU G. Cytological researches on the anther development of a thermo- sensitive genic male sterile line TE5A in., 2016, 49(12): 2408-2417. (in Chinese)

[16] TIMOFEJEVA LS, KIBBS D S, LEE S, GOLUBOVSKAYA L, WANG R, HARPER L, WALBOT V, CANDE W Z. Cytological characterization and allelism testing of anther developmental mutants identified in a screen of maize male sterile lines, 2013, 3(2):231-249.

[17] 陳雪平, 羅雙霞, 軒淑欣, 張成合, 申書興. 茄子雄性不育系花粉敗育的細胞學觀察. 植物遺傳資源學報, 2013, 14(4): 754-758.

CHEN X P, LUO S X, XUAN S X, ZHANG C H, SHEN S X. Cytological observation of pollen abortion in eggplant male-sterile lines., 2013, 14(4): 754-758. (in Chinese)

[18] 何長征, 劉志敏, 熊興耀, 鄒學校, 蕭浪濤. 辣椒細胞質(zhì)雄性不育系9704A花藥發(fā)育的細胞觀察. 園藝學報, 2008, 35(4): 521-528.

HE C Z, LIU Z M, XIONG X Y, ZHOU X X, XIAO L T. Cytologic observations on anther development of 9704A a cytoplasmic male sterile line inL., 2008, 35(4): 521-528. (in Chinese)

[19] 王濤濤, 蔡曉峰, 張俊紅, 李漢霞, 葉志彪. 芥菜型油菜雄性不育系與甘藍遠緣雜交胚培養(yǎng)及早代育性鑒定. 園藝學報, 2010, 37(10): 1661-1666.

WANG T T, CAI X F, ZHANG J H, LI H X, YE Z B. The culture and early male sterile identification of distant hybrid embryos derived fromL. and male sterile line in., 2010, 37(10): 1661-1666. (in Chinese)

[20] SULTAN M, SCHULZ M H, RICHARD H, MAGEN A, KLINGENHOFF A, SCHERF M, SEIFERT M, BORODINA T, SOLDATOV A, PARKHOMCHUK D, SCHMIDT D, O'KEEFFE S, HAAS S, VINGRON M, LEHRACH H, YASPO M L. A global view of gene activity and alternative aplicing by deep sequencing of the human transcriptome., 2008, 21: 956-960.

[21] MAHER C A, SINHAL C K, CAO X H, KALYANA-SUNDARAM S, HAN B, JING X J, SAM L, BARRETTE T, PALANISAMY N, CHINNAIYAN A M. Transcriptome sequencing to detect gene fusions in cancer., 2009, 458: 97-101.

[22] FILICHKIN S A, PRIEST H D, SCOTT A, GIVAN S A, SHEN R K, BRYANT D W, FOX S E, WONG W K, MOCKLER T C. Genome-wide mapping of alternative splicing in, 2010, 20: 45-58.

[23] SHI C Y, YANG H, WEI C L, YU O, ZHANG Z Z, JIANG C J, SUN J, LI Y Y, CHEN Q, XIA T, WAN X C. Deep sequencing of the Camellia sinensis transcriptome revealed candidate genes for major metabolic path-ways of tea-specific compounds., 2011, 12: 131.

[24] WEI W L, QI X Q, WANG L H, ZHANG Y X, HUA W, LI D H, Lü H X, ZHANG X R. Characterization of the sesame (L.) global transcriptome using illumina paired-end sequencing and development of EST-SSR markers., 2011, 12: 451.

[25] AN H, YANG Z H, YI B, WEN J, SHEN J X, TU J X, MA C Z, FU T D. Comparative transcript profiling of the fertile and sterile flower buds ofCMS in., 2014, 15: 258.

[26] LI J J, HAN S H, DING X L, HE T T, DAI J Y, YANG S P, GAI J Y. Comparative transcriptome analysis between the cytoplasmic male sterile line NJCMS1A and its maintainer NJCMS1B in soybean ((L.) Merr.)., 2015, 10(5): e126771.

[27] JEONG H J, KANG J H, ZHAO M A, KWON J K, CHOI H S, BAE J H, LEE H A, JOUNG Y H, CHOI D, KANG B C. Tomato male sterile 1035 is essential for pollen development and meiosis in anthers., 2014, 65(22): 6693-6709.

[28] CHEN C M, CHEN G J, CAO B H, LEI J J. Transcriptional profiling analysis of genic male sterile-fertilereveal candidate genes for pollen development and maturation by RNA-Seq technology., 2015, 122(2): 465-476.

[29] ZHU Q D, SONG Y L, ZHANG G S, JU L, ZHANG J, YU Y A, NIU N, WANG J W, MA S C.assembly and transcriptome analysis of wheat with male sterility induced by the chemical hybridizing agent SQ-1., 2015, 10(4): e0123556.

[30] OMIDVAR V, MOHORIANU I, DALMAY T, FELLNER M. Identification of miRNAs with potential roles in regulation of anther development and male-sterility inmale-sterile tomato mutant., 2015, 16(1): 183.

[31] TANG Z H, ZHANG L P, XU C G, YUAN S H, ZHANG F T, ZHENG Y L, ZHAO C P. Uncovering small RNA-mediated responses to cold stress in a wheat thermo-sensitive genic male-sterile line by deep sequencing., 2012, 159(2): 721-738.

[32] 陶興林, 胡立敏, 侯棟, 朱惠霞. 花椰菜溫敏雄性不育系GS-19的選育. 農(nóng)業(yè)科技通訊, 2010(11): 59-60.

TAO X L, HU L M, HOU D, ZHU H X. Breeding of the cauliflower thermo-sensitive male sterile line GS-19., 2010(11): 59-60. (in Chinese)

[33] 胡立敏, 陶興林, 朱惠霞, 劉明霞. 花椰菜新品種‘圣雪4號’. 園藝學報, 2015, 42(5): 1007-1008.

HU L M, TAO X L, ZHU H X, LIU M X. A new cauliflower cultivar ‘Shengxue 4’., 2015, 42(5): 1007-1008. (in Chinese)

[34] 葛娟, 郭英芬, 于澄宇, 張國云, 董軍剛, 董振生. 甘藍型油菜光、溫敏雄性不育系Huiyou50S花粉敗育的細胞學觀察. 作物學報, 2012, 38(3): 541-548.

GE J, GUO Y F, YU C Y, ZHANG G Y, DONG J G, DONG Z S. Cytological observation of anther development of photoperiod/thermo sensitive male sterile line Huiyou50S inL., 2012, 38(3): 541-548. (in Chinese)

[35] NIVISON H T, HANSON M R. Identification of a mitochondrial protein associated with cytoplasmic male sterility in petunia, 1989, 1(11): 1121-1130.

[36] 周長久, 張友良.蘿卜雄性不育的幾種特性研究. 園藝學報, 1994, 21(1): 13-16.

ZHOU C J. ZHANG Y L. Studies on several properties of radish male-sterility.1994, 21(1): 13-16. (in Chinese)

[37] 童旭宏, 秦利, 孔廣超, 祝水金. 新型陸地棉細胞質(zhì)雄性不育系花器形態(tài)學和細胞學觀察. 棉花學報, 2009, 21(2): 83-88.

TONG X H, QIN L, KONG G C, ZHU S J. Studies on the morphological and anatomic characters of a novel CMS inL.2009, 21(2): 83-88. (in Chinese)

[38] 張體德, 張海洋, 鄭永戰(zhàn), 衛(wèi)雙玲, 梅鴻獻. 芝麻細胞核雄性不育系ms86-1的花器形態(tài)及小孢子發(fā)育的細胞學觀察. 中國農(nóng)學通報, 2010, 26(23): 64-68.

ZHANG T D, ZHANG H Y, ZHENG Y Z, WEI S L, MEI H X. The cytological observation on the morphology of the flowers and the microspore development in sesame nuclear male-sterile lines ‘ms86-1’., 2010, 26(23): 64-68. (in Chinese)

[39] 董軍剛, 董振生, 劉絢霞, 劉創(chuàng)社, 李紅兵. 甘藍型油菜生態(tài)雄性不育系533S花藥發(fā)育的細胞學研究. 西北農(nóng)林科技大學學報(自然科學版), 2004, 32(7): 61-66.

DONG J G, DONG Z S, LIU X X, LIU C S, LI H B. Cytological studies on anther development of ecological male sterile line 533S inL., 2004, 32(7): 61-66. (in Chinese)

[40] 楊光圣, 瞿波, 傅廷棟. 三個甘藍型油菜隱性細胞核雄性不育系小孢子發(fā)生的細胞學研究. 華中農(nóng)業(yè)大學學報, 1999, 18(6): 520-523.

Yang G S, QU B. Fu T D. Anatomical studies on the anther’s development of the dominant genic male sterile line YI-3A inL., 1999, 18(5): 405-408. (in Chinese)

[41] 劉燕, 董振生, 張改文, 董軍剛, 劉創(chuàng)社, 劉絢霞, 李宏兵. 甘藍型油菜CMS 212A 花藥發(fā)育的細胞學研究. 西北農(nóng)業(yè)學報, 2005, 14(1): 33-37.

LIU Y, DONG Z S, ZHANG G W, DONG J G, LIU C S, LIU X X, LI H B. Cytological study on growth of anther of CMS 212A inL., 2005, 14(1): 33-37. (in Chinese)

[42] XIAO Z D, XIN X Y, CHEN H Y, HU S W. Cytological investigation of anther development in DGMS line Shaan-GMS inL., 2013, 49: 16-23.

[43] WAN L L, XIA X Y, HONG D F, LI J, YANG G S. Abnormal vacuolization of the tapetum during the tetrad stage is associated with male sterility in the recessive genic male sterileL. line 9012A., 2010, 53: 121-133.

[44] Dun X L, Zhou Z F, Xia S Q, WEN J, YI B, SHEN J X, MA C Z, TU J X, FU T D. Bna C.Tic40, a plastid inner membrane translocon originating from, is essential for tapetal function and microspore development in.2011, 68: 532-545.

[45] ZHOU Z F, DUN X L, XIA S Q, SHI D Y, QIN M M, YI B, WEN J, SHEN J X, MA C Z, TU J X, FU T D. BnMs3 is required for tapetal differentiation and degradation, microspore separation, and pollen wall biosynthesis in.,2012, 63: 2041-2058.

[46] 賀曉敏, 周美蘭, 余傳啟, 蔣敏敏, 茹振鋼. 溫敏雄性不育小麥BNS366花粉敗育的細胞學觀察. 麥類作物學報, 2014, 34(4): 460-466.

HE X M, ZHOU M L, YU C Q, JIANG M M, RU Z G. Cytological observation on pollen male fertility abortion of thermo-sensitive male sterile wheat line BNS366.,2014, 34(4): 460-466. (in Chinese)

[47] MA Y, KANG J G, WU J, ZHU Y G, WANG X W. Identification of tapetum-specific genes by comparing global gene expression of four different male sterile lines in.,2015, 87(6): 541-554.

[48] 樊建青, 張立平, 趙昌平, 許晨光, 王靈云, 苑少華. 光溫敏核雄性不育小麥 BS366 花粉母細胞減數(shù)分裂的細胞學研究. 中國細胞生物學學報, 2011, 33(6): 622-628.

FAN J Q, ZHANG L P, ZHAO C P, XU C G, WANG L Y, YUAN S H. Studies on mieosis pollen mother cells in photoperiod-temperature sensitive genic male sterile wheat line BS366.,2011, 33(6): 622-628. (in Chinese)

[49] 梁小紅, 儀治本, 趙威軍, 段運平, 崔貴梅, 孫毅. 高粱A2型細胞質(zhì)雄性不育系小孢子發(fā)生的細胞學觀察和減數(shù)分裂染色體行為分析. 作物學報, 2006, 32(8): 1107-1110.

LIANG X H, YI Z B, ZHAO W J, DUAN Y P, CUI G M, SUN Y. Cytological observation of microsporogenesis and its chromosomal behavior in meiosis of A2 cytoplasmic-male sterile line in sorghum., 2006, 32(8): 1107-1110.(in Chinese)

[50] 王芳, 衛(wèi)保國, 李貴全, 李艷花. 大豆光敏雄性不育株88-428BY-827小孢子母細胞的細胞學觀察. 中國農(nóng)業(yè)科學, 2004, 37(8): 1110-1113.

WANG F, WEI B G, LI G Q, LI Y H. A cytological observation of the pollen mother cells of the photoperiod- sensitive male sterile soybean plant of 88-428BY-827., 2004, 37(8): 1110-1113. (in Chinese)

[51] 田福發(fā), 徐躍進, 袁黎, 陳建軍. 紅菜薹雄性不育系花藥敗育的細胞形態(tài)學觀察. 武漢植物學研究, 2004, 22(3): 269-272.

TIAN F F, XU Y J, YUAN L, CHEN J J. The cytom orphology study of anther abortion in male-sterile lines of Hong Cai-tai (L. ssp.LLee.).,2004, 22(3): 269-272. (in Chinese)

[52] 賀國良, 賀浩華, 劉宜柏. 光敏核不育水稻農(nóng)墾58S細胞學研究. 江西農(nóng)業(yè)大學學報, 2001, 23(1): 24-27.

HE G L, HE H H, LIU Y B. The cytological ctudy of photoperiod- sensitive genic sterile rice Nongken 58S.,2001, 23(1): 24-27. (in Chinese)

[53] HUANG L, YANG X, SUN P, TONG W, HU S. First Illumina-based de novo transcriptome sequencing and analysis of safflower flowers.2012, 7: e38653. doi: 10.1371/journal.pone.0038653 PMID: 22723874.

[54] YONEKURA-SAKAKIBARA K, NAKABAYASHI R, SUGAWARA S, TOHGE T, ITO T, KOYANAGI M, KITAJIMA M, TAKAYAMA H,SAITO K. A flavonoid 3--glucoside: 2"--glucosyltransferase responsible for terminal modification of pollen-specific flavonols in, 2014, 79(5): 769-782.

[55] CHEN L T, LIU Y G. Male sterility and fertility restoration in crops., 2014, 65(1): 579-606.

[56] LI J J, HAN S H, DING X L, HE T T, DAI J Y, YANG S P, GAI J Y. Comparative transcriptome analysis between the cytoplasmic male sterile line NJCMS1A and its maintainer NJCMS1B in soybean ((L.) Merr.)., 2015, 10(5): e126771.

[57] SHEMESH-MAYER E, BEN-MICHAEL T, ROTEM N, RABINOWITCH H D, DORON-FAIGENBOIM A, KOSMALA A, PERLIKOWSKI D, SHERMAN A, KAMENETSKY R. Garlic (L.) fertility: transcriptome and proteome analyses provide insight into flower and pollen development., 2015, 6: 271.

[58] MEI S Y, LIU T M, WANG Z W. Comparative transcriptome profile of the cytoplasmic male sterile and fertileral buds of radish (L),2016, 17(1): 42.

[59] AN H, YANG Z H, YI B, WEN J, SHEN J X, TU J X, MA C Z, FU T D. Comparative transcript profiling of the fertile and sterile flower buds ofCMS in.., 2014, 15(1): 258.

[60] QU C M, FU F Y, LIU M, ZHAO H Y, LIU C, LI J N, TANG Z L, XU X F, QIU X, WANG R, LU K. Comparative transcriptome analysis of recessive male sterility (RGMS) in sterile and fertilelines., 2015, 10(12): e0144118.

（責任編輯 趙伶俐）

Transcriptome and Cytological Researches on the Anther Abortion of a Thermo-sensitive Genic Male Sterile Line GS-19 in Cauliflower

TAO XingLin1,2, HOU Dong1, ZHU HuiXia1, LIU MingXia1, ZHANG JinWen2, HU LiMin1

(1Institute of Vegetable, Gansu Academy of Agricultural Science/Lanzhou Research Station of Horticultural Crop Biology and Germplasm Enhancement, Ministry of Agriculture, Lanzhou 730070;2Gansu Key Lab of Crop Improvement & Germplasm Enhancement/Gansu Agricultural University, Lanzhou 730070)

【Objective】In order to clarify the mechanism of cellular and molecular sterility characteristics, the flower morphology, anther abortion period and molecular mechanism of the thermos-sensitive male sterile line of cauliflower were studied, thus laying a theoretical foundation for further study of the male sterility. 【Method】The thermo-sensitive male sterile line GS-19 was used as the experimental material, using morphological, cytological and high-throughput sequencing techniques, the morphological, cytological and molecular mechanism were studied with treatment at different temperatures (15℃ and 20℃). 【Result】The fertility conversion of the thermo-sensitive male sterile line GS-19 was controlled by temperature, sterility at high temperature (20℃), and fertility at low temperature (15℃).The flower difference between sterile and fertile plants reached a significant level in GS-19, the flower was significantly smaller than that of fertile plants. The anthers abortive stages and manners were different between sterile and fertile plants of GS-19by the microstructureobservation.The anther development of GS-19had differentiation of pollenmother cell, formed normal pollen sac, did not produce pollen grains or produce trace ofpollen grains withoutviability, the developmentof antherswas blocked from pollen mother cell to the tetradstage. The pollen mother cells could not pass the meiosis, with no dyads and tetrads formed in GS-19, insteading of some “pseudo microspores”, belong to an abortive type of pollen mother cell. The pollen exines of the “pseudo microspores” were abnormal. The “pseudo microspores” were decayed gradually and only empty shells remained at last. A total of 67 930 unigenes were generated in GS-19 alabastrum by sequence, Nt database compared to the largest number of similar sequences, for the 52 191, Nr database to similar sequence of times (46 447), KOG database than similar sequence of the least (13 198). By GO database the 25 336 unigenes were divided into 3 categories containing 55 function groups; by KOG databases the 13 198 unigenes were grouped into 26 functional categories; by KEGG database the 14 778 unigenes were divided into 19 metabolism pathways.Gene differential expression analysis revealed that there were 2 170 differentially expressed genes in the bud of GS-19 male sterile line and fertile plants, 1 078 genes up-regulated and 1 092 genes down regulated. 【Conclusion】The fertility conversion of the thermo-sensitive male sterile line GS-19 is controlled by temperature, sterility at high temperature. The flower is significantly smaller than that of fertile plants. The filaments are significantly shortened, the anthers shrank, and dried. The developmentof anthersis blocked from pollen mother cell to the tetradstage, belonging to an abortive type of pollen mother cell. A total of 67 930 unigenes were generated in GS-19 alabastrum by transcriptome analysis. Gene differential expression analysis reveals that there are 2 170 differentially expressed genes in the bud of GS-19 male sterile and fertile plants.

cauliflower (L.); thermo-sensitive genic male sterile; anther abortion; cytology; transcripiome analysis

2016-12-15；接受日期：2017-04-13

國家自然科學基金（31460519）、農(nóng)業(yè)部園藝作物生物學與種質(zhì)創(chuàng)制西北地區(qū)科學觀測試驗站項目（2015-A2621-620321-G1203-066）、現(xiàn)代農(nóng)業(yè)產(chǎn)業(yè)技術體系（CARS-25-G-48）、甘肅省農(nóng)科院科技創(chuàng)新工程學科團隊（2014GAAS07）

陶興林，E-mail：taoxinglin77@126.com