亞基因組的平行趨同選擇是導致蕓薹和甘藍的形態(tài)型多樣化以及趨同馴化的主要動(dòng)力
2016年8月15日,《自然-遺傳學(xué)》(Nature Genetics)雜志在線(xiàn)發(fā)表了中國農業(yè)科學(xué)院蔬菜花卉研究所王曉武研究組和北京百邁客生物科技有限公司合作的研究成果。借助重測序技術(shù)以及蕓薹屬祖先的染色體核型(translocation Proto-Calepineae Karyotype,tPCK),他們對蕓薹和甘藍的形態(tài)型的平行和趨同進(jìn)化現象進(jìn)行了深入研究,該研究最主要的成果是國際上首次發(fā)現全基因組三倍化事件是導致蕓薹類(lèi)物種的形態(tài)型多樣化以及趨同馴化的原因;此外挖掘出兩個(gè)重要的商業(yè)性狀——葉片結球性狀和塊莖形態(tài)性狀的重要調控基因,為今后的分子育種工作提供遺傳基礎。
蕓薹屬包括三個(gè)二倍體物種:蕓薹(Brassica rapa)、黑芥(B. nigra)和甘藍(B. oleracea),這三個(gè)基本種經(jīng)過(guò)兩兩雜交又產(chǎn)生了四倍體的歐洲油菜(B. napus)、非洲油菜(B. carinata)和芥菜(B. juncea),這就是“禹氏三角”。隨著(zhù)人工馴化和持續的育種工作,這些物種產(chǎn)生了高度多樣化的形態(tài)型,例如結成球狀的葉片、膨大的塊莖(或根部、頂芽、腋芽等)。有意思的是,某些物種的形態(tài)型,盡管是在不同的地區獨立培育而成,但卻表現出相似的表型特征,這便是趨同馴化的結果。
圖1.蕓薹屬禹氏三角(From Wikipedia)
導致蕓薹屬種內形態(tài)型多樣化以及種間趨同馴化的機制是什么?這是科學(xué)家們一直在探討的問(wèn)題。然而蕓薹屬的起源進(jìn)化歷史也是相當的“波折”——蕓薹屬的祖先是一個(gè)具有7條染色體的二倍體,在約1200萬(wàn)年前發(fā)生過(guò)一次全基因組三倍化事件,該事件導致一個(gè)具有42條染色體的古六倍體物種出現。隨后這個(gè)六倍體物種的基因組發(fā)生了廣泛的二倍體化——有時(shí)候一條染色體的一段區域轉移到另一條染色體上,有時(shí)候某個(gè)區域會(huì )丟失,更有時(shí)候連整條染色體都會(huì )消失。最終六倍體的又變回了二倍體的“模樣”,形成了現在的蕓薹屬二倍體物種分類(lèi)格局?!绱藦碗s的進(jìn)化歷史自然帶來(lái)不少的阻礙,所幸的是,借助于不斷發(fā)展的高通量測序技術(shù)和生物信息分析方法,科學(xué)家們能夠從基因組層面解析物種的進(jìn)化機制成為可能。
在這項研究中,王曉武研究組選擇不同形態(tài)型的蕓薹199株和甘藍119株作為實(shí)驗材料,這些材料包括不同地理區域分布的13個(gè)蕓薹及9個(gè)甘藍亞群,其中涉及葉片結球型(56份大白菜,45份卷心菜),以及塊莖膨大型(54份大頭菜及19份苤藍)。應用的是Illumina?Hiseq2000測序平臺,350bp插入片段文庫,每個(gè)個(gè)體的平均測序深度均大于8X。
圖2. 系統發(fā)育樹(shù)構建
為了調查不同形態(tài)的品種的馴化歷史,應用兩個(gè)群體共有的6707個(gè)SNP位點(diǎn)構建系統發(fā)育樹(shù),如圖2所示,蕓薹可以分成6個(gè)組,其中蕪菁(turnip)和野油菜(sarsons)等位于系統發(fā)育樹(shù)基部;大白菜(Chinese cabbage)位于最遠端,表明馴化歷史最短。甘藍群體可以劃分成7個(gè)組,其中卷心菜(cabbage)位于最遠端。從進(jìn)化樹(shù)不難看出,大白菜和卷心菜雖然具有相似的葉片結球特征,卻是不同的祖先經(jīng)過(guò)平行或趨同馴化產(chǎn)生的;相類(lèi)似地,蕪菁和苤藍具有相似的塊莖膨大特征,也同樣是平行或趨同馴化的結果。
分別對蕓薹和甘藍種內結球和不結球的群體進(jìn)行分析——主要通過(guò)ROD和PiHS兩個(gè)指數,共同篩選具有正選擇信號的基因區域。結果顯示蕓薹中有20個(gè)區域是ROD和PiHS共同定位到的,甘藍中是18個(gè)。進(jìn)一步對這些區域進(jìn)行基因注釋和富集,發(fā)現4個(gè)與植物激素應答相關(guān)的GO類(lèi)別,這些植物激素對葉片形狀和極性非常重要;此外還發(fā)現影響葉片弧度的一些基因同樣受到正選擇。
進(jìn)一步檢測古六倍體的亞基因組對蕓薹和甘藍結球性狀的平行或趨同進(jìn)化的影響,作者將兩個(gè)物種的基因組打碎成基因單元,再根據已知的祖先染色體核型,將上述基因單元重排和串接,分別構建出蕓薹的三個(gè)亞基因組一致序列和甘藍的三個(gè)亞基因組一致序列。重新進(jìn)行PiHS分析后,發(fā)現蕓薹的不同亞基因組之間存在4個(gè)平行選擇信號的區域,甘藍存在4個(gè);同時(shí)蕓薹和甘藍的相同亞基因組之間存在著(zhù)15個(gè)受到趨同選擇的區域。
隨后針對塊莖膨大性狀進(jìn)行的研究亦獲得了類(lèi)似的結果。由此證明,蕓薹族祖先的全基因組三倍化事件可能對蕓薹屬性狀平行或趨同進(jìn)化的造成深遠的影響,這一發(fā)現在國際上尚屬首次。
圖3. ?蕓薹(左)和甘藍(右)亞基因組的選擇信號檢測
參考文獻
ChengF.?et?al.?(2016).?“Subgenome?parallel?selection?is?associated?with?morphotype?diversification?and?convergent?crop?domestication?in?Brassica?rapa?and?Brassica?oleracea.”?Nat?Genet?advance?online?publication.