材料方法
結(jié)果
1.?基因組組裝
對(duì)油用亞麻品種Longya-10、纖維亞麻品種Heiya-14和pale flax分別測(cè)序得到68.2Gb,73.5Gb和 49.1Gb數(shù)據(jù),組裝出306.0Mb、303.7Mb和293.5 Mb基因組,其中contig N50/scaffold N50分別為131Kb/ 1,235Kb, 156Kb/700Kb和59Kb/384 Kb,后期利用HiC技術(shù)與遺傳圖譜輔助Longya-10基因組組裝,將434 scaffolds組裝至染色體水平,在每個(gè)基因組中預(yù)測(cè)到大約43,500個(gè)編碼基因和2600-2800個(gè)非編碼RNA;系統(tǒng)發(fā)育分析表明,栽培亞麻和白亞麻在大約232萬(wàn)年前發(fā)生了分化,自古被子植物的六倍化事件以來(lái),亞麻發(fā)生了兩次全基因組復(fù)制事件。


作物改良常常導(dǎo)致受選擇的區(qū)域核酸多態(tài)性降低,將油用亞麻群體和纖維亞麻群體分別與地方種比較分析,油用亞麻中檢測(cè)到108個(gè)選擇清除區(qū)域(長(zhǎng)度為15.5Mb,1,958個(gè)基因),纖維亞麻中有60個(gè)選擇清除區(qū)域(長(zhǎng)度為8.2Mb,1018個(gè)基因),27個(gè)選擇清除區(qū)域重疊。其中α生物素羧基載體蛋白(LuBCCP)、脂氧合酶(LuLOX)、脂肪?;鵄CP硫酯酶A(LuFatA)、脂質(zhì)轉(zhuǎn)移蛋白(LuLTP)和丙酮酸脫氫酶復(fù)合物E2組分(LuPDH-E2)等油脂相關(guān)基因以及種子大小相關(guān)基因LuBIN2和LuGW5在油用亞麻中受到選擇。纖維亞麻中,次生細(xì)胞壁生物合成相關(guān)基因(LuMYB46-1、LuXTH和LuROPGAP3)和植物株高相關(guān)基因(LuGA3ox、LuGA20ox和LuGID)受到選擇。
之前有報(bào)道纖維亞麻是從油用亞麻馴化而來(lái),將油用亞麻與纖維亞麻進(jìn)行比較分析,存在47個(gè)選擇清除區(qū)域,其中50.9%與纖維亞麻相對(duì)于地方種的選擇區(qū)域重合。這表明在纖維用亞麻改良過(guò)程中,這些區(qū)域持續(xù)受到很強(qiáng)的選擇壓力。大約一半的基因(426/867個(gè)基因)僅在油用亞麻和纖維亞麻比較中發(fā)現(xiàn),包含內(nèi)切-β-1,4-葡聚糖酶(LuKorrigan)、果膠甲酯酶(LuPME)和焦磷酸銅合成酶(LuCPS)等。作者將選擇清除區(qū)域與之前報(bào)道的QTL/GWAS位點(diǎn)比較分析,發(fā)現(xiàn)2個(gè)硬脂酸相關(guān)的QTL與油用亞麻選擇區(qū)域重疊,一個(gè)株高相關(guān)的GWAS位點(diǎn)與纖維亞麻選擇區(qū)域重疊,另外3個(gè)油合成相關(guān)的QTL與纖維亞麻選擇區(qū)域重疊,這說(shuō)明亞麻在馴化過(guò)程中受到雙重選擇。
纖維是植物中一種特殊的細(xì)胞,具有加厚的次生細(xì)胞壁。AtMYB83/MYB46是擬南芥次生細(xì)胞壁生物合成的主要調(diào)控因子。研究發(fā)現(xiàn),在亞麻、楊樹、蘋果、苜蓿和木薯中,MYB46/MYB83基因發(fā)生了物種特異性復(fù)制。8個(gè)已鑒定的LuMYB46/LuMYB83同源基因,4個(gè)在栽培種中高表達(dá)。此外,LuMYB46-1、LuMYB43-2和LuMYB83-1在栽培亞麻中發(fā)生了多種變異。與白亞麻相比,LuMYB83-1在兩個(gè)栽培種中檢測(cè)到21bp的插入,LuMYB46-1在亞麻改良期間受到了強(qiáng)烈的選擇。LuMYB46-2在兩個(gè)栽培種之間也有不同的插入/缺失變異。在亞麻中,四對(duì)MYB46/MYB83姊妹基因位于共線基因組區(qū)域,最近的一次分裂發(fā)生在最近一次WGD發(fā)生的時(shí)間,這意味著這次WGD事件導(dǎo)致了MYB46/MYB83基因在亞麻基因組中新的一次擴(kuò)增。這些結(jié)果表明,在亞麻馴化和改良過(guò)程中,LuMYB46/LuMYB83基因的擴(kuò)張和定向選擇可能在重塑亞麻形態(tài)結(jié)構(gòu)中起著至關(guān)重要的作用,qRT-PCR與RNA-Seq數(shù)據(jù)再次證明此結(jié)論。