12月12 【成功案例】云平臺助力蕪菁花芽轉錄組比較分析解析多倍體異常減數(shù)分裂過程

ranscriptome?Analysis?of?Floral?Buds?Deciphered?an?Irregular?Course?of?Meiosis?in?Polyploid?Brassica?rapa?

云平臺助力蕪菁花芽轉錄組比較分析解析多倍體異常減數(shù)分裂過程

雜志：?Frontiers?in?Plant?Science?

影響因子：4.298

PMID:28553302

云平臺助力鄭州大學的老師們發(fā)表了一篇關于四倍體與二倍體蕪菁轉錄組比較分析的文章，這是第一個同時在細胞和轉錄組水平證明染色體多倍性對減數(shù)分裂過程有不利影響的研究，有助于全面了解二倍體和多倍體蕪菁減數(shù)分裂時花芽轉錄組的一致性和差異性。

 

研究背景：

多倍體化對于動植物進化及物種形成具有重要意義。染色體倍性增加在大多數(shù)植物中十分常見，約30-80％的被子植物在進化過程中歷經這一過程，種內基因組復制可產生同源多倍體，種間雜交可產生異源多倍體。理解多倍體對植物繁殖的影響對于多倍體育種項目具有重要意義。

本文在細胞和分子水平對同源四倍體、二倍體蕪菁的生殖組織（未成熟的花芽）進行了轉錄組比較分析。先進行細胞學分析，發(fā)現(xiàn)結果顯著異常后，進行RNA測序分析，進一步闡明這種細胞學差異的分子機制。

 

材料和方法:

材料：溫室中培育蕪菁同源四倍體和二倍體，采集1-1.5mm的花芽。

細胞學分析：花芽制片觀察染色體行為。

轉錄組分析：二倍體花芽樣本T1，T2，T3,同源四倍體樣本T4，T5，T6進行測序，篩選DEGs，進行功能注釋及相關驗證。

 

測序平臺：Illumina?HiSeq?2000

分析平臺：BMKCloud，具體分析如下：

測序原始數(shù)據(jù)去接頭、低質量序列后與參考基因組比對（TopHat2），對比對上的序列進行轉錄本重構，評估表達水平（Cufflinks）。鑒定DEGs（logFC≥2，F(xiàn)DR<0.01）并與Nr、UniProtKB\Swissprot、KOG、COG、GO、KEGG等數(shù)據(jù)庫進行比對注釋。

 

研究結果：

1.同源四倍體蕪菁異常的減數(shù)分裂過程

減數(shù)分裂過程中染色體行為分析結果顯示：二倍體中，同源染色體中期配對為二倍體，后期分離（圖1A-D）。四倍體中，在減數(shù)分裂中期I形成多倍體和單倍體，在后期I、II發(fā)生不均等分離（圖1E-L）。

圖1?四倍體蕪菁異常染色體行為

統(tǒng)計結果表明與二倍體相比，蕪菁多倍體化后整個減數(shù)分裂過程發(fā)生了不規(guī)則變化。

圖2?四倍體蕪菁減數(shù)分裂中異常染色體行為分析

2.測序結果

每個樣本文庫得到29.07?GB?clean?read，Q30≥88.55%?。同源四倍體、二倍體比對到參考基因組效率分別為74.34%，75.88%。樣本組的表達模式分析共得到40927個基因，其中有1001個新基因。

3.DEGs分析

40927個基因中，4601個基因是差異表達的（圖3A），2343個上調，2259個下調（圖3B），兩組間DEGs占比較少（(11.24%）。隨機挑選DEGs進行層次聚類，分析表達模式（圖3C），發(fā)現(xiàn)四倍體中K1、K3、K6、K9類別下調，K2、K4、K5、K7，K8類別顯著上調。

圖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%），轉錄（261個,9.98%），一般功能（515個,19.69%)，與KEGG結果一致。

GO分類中大多數(shù)DEGs聚類到細胞過程（61.97％），繁殖過程（15.49％），結合過程（42.96％）（圖5）。

KEGG功能注釋結果表明，4,601個DEGs中有1,453個可歸到50個生物途徑。高度富集的通路為：氨基酸合成（48個，3.3％），植物激素信號轉導（61個，4.2％）和蛋白質加工內質網途徑（47個，3.2％）（圖6）。

圖4?COG分類

圖5?新DEGs的GO分類

圖6?KEGG富集通路

4.減數(shù)分裂相關DEGs分析

為闡述四倍體蕪菁減數(shù)分裂過程改變相關的遺傳信息，選擇以下明顯與減數(shù)分裂相關的COG類別進行細分：（1）復制，重組和修復;（2）染色質結構和動力學;（3）細胞周期控制，細胞分裂，染色體分區(qū);（4）細胞骨架。在4,601個位點中，共鑒定出288個與減數(shù)分裂相關的基因（圖7A）。
應用擬南芥數(shù)據(jù)庫（TAIR）的BLASTN搜索，進一步鑒定蕪菁中減數(shù)分裂直系同源基因，11個兩組間差異表達的已知減數(shù)分裂基因。如圖7所示，顯著富集的L組（復制，重組和修復）包含121個上調、102個下調基因，這組聚類結果（圖8B）顯示已知減數(shù)分裂基因（包括RAD54，DMC1和RPA）在四倍體蕪菁中表達顯著下調。D組（細胞周期控制，細胞分裂和染色體分配）中，10個下調，13個上調，已知的減數(shù)分裂基因DIF1/SYN1和CyclinA1-2/TAM上調（圖7C）。在B組（染色質結構和動態(tài)變化）（圖7D），有7個基因上調，包括已知調節(jié)開花時間和植物發(fā)育的NF-YB8基因，14個基因下調，包括組蛋白相關基因（H2AV和H4）。Z組（細胞骨架）中，8個基因上調，7個基因下調（圖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蛋白質引起的DSBs的精確修復十分重要，作為調節(jié)關鍵介質的DMC1，RAD54和RPA在這個途徑中顯著下調。Cyclin-A1-2?/?TAM是DNA復制途徑中唯一與細胞周期蛋白有關的基因，表達上調。這些結果表明染色體集在同源四倍體中增加了一倍，與減數(shù)分裂有關的大多數(shù)基因可能通過上調或至少正常表達以滿足減數(shù)分裂期間基因組含量增加的需求。否則，減數(shù)分裂關鍵基因下調可能會干擾四倍體減數(shù)分裂中染色體行為。

qRT-PCR?對24個減數(shù)分裂相關基因進行驗證，20多個基因表達與測序結果一致。對擬南芥中直系同源基因分析結果顯示兩個物種減數(shù)分裂相關基因的總體表達模式相近。

 

結果：

本研究細胞學分析發(fā)現(xiàn)，合成多倍體整體減數(shù)分裂過程是顯著不規(guī)則的。為從分子水平闡明這一過程的遺傳基礎，本文進行了同源四倍體和二倍體蕪菁花芽間的遺傳調節(jié)差異的轉錄組比較分析。在40927個表達基因中，同源四倍體蕪菁花芽中鑒定出4601個DEGs，其中288個與減數(shù)分裂相關。DMC1（已知減數(shù)分裂特異性基因，與DSBs同源染色體依賴性修復相關）在同源四倍體蕪菁中表達顯著下調，可能與減數(shù)分裂I期的異常有關。某些RNA解旋酶，細胞周期、體細胞DNA修復相關的DEG在基因組復制后表達上調，減數(shù)分裂DSB修復有關基因表達顯著下調。蕪菁和擬南芥中減數(shù)分裂相關基因的總體表達模式相近。

 

創(chuàng)新點：

這是第一個在細胞和轉錄組水平證明染色體多倍性對減數(shù)分裂過程有不利影響的研究，有助于全面了解二倍體和多倍體蕪菁減數(shù)分裂時花芽轉錄組的一致性和差異性。

 

參考文獻：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.