12月12 【成功案例】云平臺助力蕪菁花芽轉(zhuǎn)錄組比較分析解析多倍體異常減數(shù)分裂過程
ranscriptome?Analysis?of?Floral?Buds?Deciphered?an?Irregular?Course?of?Meiosis?in?Polyploid?Brassica?rapa?
云平臺助力蕪菁花芽轉(zhuǎn)錄組比較分析解析多倍體異常減數(shù)分裂過程
雜志：?Frontiers?in?Plant?Science?
影響因子：4.298
PMID:28553302
云平臺助力鄭州大學的老師們發(fā)表了一篇關于四倍體與二倍體蕪菁轉(zhuǎn)錄組比較分析的文章，這是第一個同時在細胞和轉(zhuǎn)錄組水平證明染色體多倍性對減數(shù)分裂過程有不利影響的研究，有助于全面了解二倍體和多倍體蕪菁減數(shù)分裂時花芽轉(zhuǎn)錄組的一致性和差異性。
研究背景：
多倍體化對于動植物進化及物種形成具有重要意義。染色體倍性增加在大多數(shù)植物中十分常見，約30-80％的被子植物在進化過程中歷經(jīng)這一過程，種內(nèi)基因組復制可產(chǎn)生同源多倍體，種間雜交可產(chǎn)生異源多倍體。理解多倍體對植物繁殖的影響對于多倍體育種項目具有重要意義。
本文在細胞和分子水平對同源四倍體、二倍體蕪菁的生殖組織（未成熟的花芽）進行了轉(zhuǎn)錄組比較分析。先進行細胞學分析，發(fā)現(xiàn)結(jié)果顯著異常后，進行RNA測序分析，進一步闡明這種細胞學差異的分子機制。
材料和方法:
材料：溫室中培育蕪菁同源四倍體和二倍體，采集1-1.5mm的花芽。
細胞學分析：花芽制片觀察染色體行為。
轉(zhuǎn)錄組分析：二倍體花芽樣本T1，T2，T3,同源四倍體樣本T4，T5，T6進行測序，篩選DEGs，進行功能注釋及相關驗證。
測序平臺：Illumina?HiSeq?2000
分析平臺：BMKCloud，具體分析如下：
測序原始數(shù)據(jù)去接頭、低質(zhì)量序列后與參考基因組比對（TopHat2），對比對上的序列進行轉(zhuǎn)錄本重構(gòu)，評估表達水平（Cufflinks）。鑒定DEGs（logFC≥2，F(xiàn)DR<0.01）并與Nr、UniProtKB\Swissprot、KOG、COG、GO、KEGG等數(shù)據(jù)庫進行比對注釋。
研究結(jié)果：
1.同源四倍體蕪菁異常的減數(shù)分裂過程
減數(shù)分裂過程中染色體行為分析結(jié)果顯示：二倍體中，同源染色體中期配對為二倍體，后期分離（圖1A-D）。四倍體中，在減數(shù)分裂中期I形成多倍體和單倍體，在后期I、II發(fā)生不均等分離（圖1E-L）。
圖1?四倍體蕪菁異常染色體行為
統(tǒng)計結(jié)果表明與二倍體相比，蕪菁多倍體化后整個減數(shù)分裂過程發(fā)生了不規(guī)則變化。
圖2?四倍體蕪菁減數(shù)分裂中異常染色體行為分析
2.測序結(jié)果
每個樣本文庫得到29.07?GB?clean?read，Q30≥88.55%?。同源四倍體、二倍體比對到參考基因組效率分別為74.34%，75.88%。樣本組的表達模式分析共得到40927個基因，其中有1001個新基因。
3.DEGs分析
40927個基因中，4601個基因是差異表達的（圖3A），2343個上調(diào)，2259個下調(diào)（圖3B），兩組間DEGs占比較少（(11.24%）。隨機挑選DEGs進行層次聚類，分析表達模式（圖3C），發(fā)現(xiàn)四倍體中K1、K3、K6、K9類別下調(diào)，K2、K4、K5、K7，K8類別顯著上調(diào)。
圖3?DEGs表達分析
與COG、GO、KOG、Nr、KEGG、Swiss-Prot比對后，97%的DEGs至少能與一個數(shù)據(jù)庫比對上（表1）
表1?差異表達Unigene注釋
COG分類中，2615個DEGs可歸到25個COG類別（圖4），其中包含最多DEGs的類別為：復制重組與修復（222個,8.49%），轉(zhuǎn)錄（261個,9.98%），一般功能（515個,19.69%)，與KEGG結(jié)果一致。
GO分類中大多數(shù)DEGs聚類到細胞過程（61.97％），繁殖過程（15.49％），結(jié)合過程（42.96％）（圖5）。
KEGG功能注釋結(jié)果表明，4,601個DEGs中有1,453個可歸到50個生物途徑。高度富集的通路為：氨基酸合成（48個，3.3％），植物激素信號轉(zhuǎn)導（61個，4.2％）和蛋白質(zhì)加工內(nèi)質(zhì)網(wǎng)途徑（47個，3.2％）（圖6）。
圖4?COG分類
圖5?新DEGs的GO分類
圖6?KEGG富集通路
4.減數(shù)分裂相關DEGs分析
為闡述四倍體蕪菁減數(shù)分裂過程改變相關的遺傳信息，選擇以下明顯與減數(shù)分裂相關的COG類別進行細分：（1）復制，重組和修復;（2）染色質(zhì)結(jié)構(gòu)和動力學;（3）細胞周期控制，細胞分裂，染色體分區(qū);（4）細胞骨架。在4,601個位點中，共鑒定出288個與減數(shù)分裂相關的基因（圖7A）。
應用擬南芥數(shù)據(jù)庫（TAIR）的BLASTN搜索，進一步鑒定蕪菁中減數(shù)分裂直系同源基因，11個兩組間差異表達的已知減數(shù)分裂基因。如圖7所示，顯著富集的L組（復制，重組和修復）包含121個上調(diào)、102個下調(diào)基因，這組聚類結(jié)果（圖8B）顯示已知減數(shù)分裂基因（包括RAD54，DMC1和RPA）在四倍體蕪菁中表達顯著下調(diào)。D組（細胞周期控制，細胞分裂和染色體分配）中，10個下調(diào)，13個上調(diào)，已知的減數(shù)分裂基因DIF1/SYN1和CyclinA1-2/TAM上調(diào)（圖7C）。在B組（染色質(zhì)結(jié)構(gòu)和動態(tài)變化）（圖7D），有7個基因上調(diào)，包括已知調(diào)節(jié)開花時間和植物發(fā)育的NF-YB8基因，14個基因下調(diào)，包括組蛋白相關基因（H2AV和H4）。Z組（細胞骨架）中，8個基因上調(diào)，7個基因下調(diào)（圖7E），且沒有已知減數(shù)分裂相關基因。已知減數(shù)分裂基因MMD1和HOP2在一般功能組（R）中富集（圖4）。HOP2，XRI1，ZYP1a和ZYP1b未歸到任何COG組。
KEGG比對后將288個中的108個DEGs分配到28個途徑，包括減數(shù)分裂相關同源重組（2.8％），DNA修復和重組（6.5％），DNA復制（3.7％），染色體和相關蛋白（4.6％）。
圖7?COG分類中減數(shù)分裂相關基因分布及表達分析
同源重組對于真核生物減數(shù)分裂中SPO11蛋白質(zhì)引起的DSBs的精確修復十分重要，作為調(diào)節(jié)關鍵介質(zhì)的DMC1，RAD54和RPA在這個途徑中顯著下調(diào)。Cyclin-A1-2?/?TAM是DNA復制途徑中唯一與細胞周期蛋白有關的基因，表達上調(diào)。這些結(jié)果表明染色體集在同源四倍體中增加了一倍，與減數(shù)分裂有關的大多數(shù)基因可能通過上調(diào)或至少正常表達以滿足減數(shù)分裂期間基因組含量增加的需求。否則，減數(shù)分裂關鍵基因下調(diào)可能會干擾四倍體減數(shù)分裂中染色體行為。
qRT-PCR?對24個減數(shù)分裂相關基因進行驗證，20多個基因表達與測序結(jié)果一致。對擬南芥中直系同源基因分析結(jié)果顯示兩個物種減數(shù)分裂相關基因的總體表達模式相近。
結(jié)果：
本研究細胞學分析發(fā)現(xiàn)，合成多倍體整體減數(shù)分裂過程是顯著不規(guī)則的。為從分子水平闡明這一過程的遺傳基礎，本文進行了同源四倍體和二倍體蕪菁花芽間的遺傳調(diào)節(jié)差異的轉(zhuǎn)錄組比較分析。在40927個表達基因中，同源四倍體蕪菁花芽中鑒定出4601個DEGs，其中288個與減數(shù)分裂相關。DMC1（已知減數(shù)分裂特異性基因，與DSBs同源染色體依賴性修復相關）在同源四倍體蕪菁中表達顯著下調(diào)，可能與減數(shù)分裂I期的異常有關。某些RNA解旋酶，細胞周期、體細胞DNA修復相關的DEG在基因組復制后表達上調(diào)，減數(shù)分裂DSB修復有關基因表達顯著下調(diào)。蕪菁和擬南芥中減數(shù)分裂相關基因的總體表達模式相近。
創(chuàng)新點：
這是第一個在細胞和轉(zhuǎn)錄組水平證明染色體多倍性對減數(shù)分裂過程有不利影響的研究，有助于全面了解二倍體和多倍體蕪菁減數(shù)分裂時花芽轉(zhuǎn)錄組的一致性和差異性。
參考文獻：Braynen?J,?Yang?Y,?Wei?F,?et?al.?Transcriptome?Analysis?of?Floral?Buds?Deciphered?an?Irregular?Course?of?Meiosis?in?Polyploid?Brassica?rapa[J].?Frontiers?in?Plant?Science,?2017,?8:768.