劉芳瑛，陸　贏，王　卓，劉顯軍，胡學(xué)芳，劉旭東，徐海鵬，劉忠松

(湖南農(nóng)業(yè)大學(xué)油料作物研究所，長(zhǎng)沙 410128)

?

芥菜型油菜A9染色體BAC物理圖譜構(gòu)建及比較基因組分析

劉芳瑛，陸贏，王卓，劉顯軍，胡學(xué)芳，劉旭東，徐海鵬，劉忠松*

(湖南農(nóng)業(yè)大學(xué)油料作物研究所，長(zhǎng)沙 410128)

為獲得高精確度的芥菜型油菜A9染色體物理圖，本研究利用定位在A9染色體上的標(biāo)記對(duì)芥菜型油菜BAC文庫(kù)進(jìn)行三維混合PCR篩選，共得到陽(yáng)性BAC 3200個(gè)；構(gòu)建了由16個(gè)重疊群、5個(gè)單拷貝組成的芥菜型油菜A9染色體物理圖譜。根據(jù)最小路徑原則挑選代表整個(gè)物理圖的種子BAC 435個(gè)，并進(jìn)行雙末端測(cè)序，得到質(zhì)量良好的序列862條；按照BAC平均插入片段為126 kb，可知所有重疊群累積達(dá)52.08 M，覆蓋A9染色體的1.37倍。利用BLAST軟件(按參數(shù)值E<1e-10)將種子BAC末端序列比對(duì)到白菜和甘藍(lán)型油菜A9染色體假分子上，結(jié)果顯示該物理圖譜分別覆蓋了36.4、31.0 Mb的物理距離，相當(dāng)于全基因組組裝長(zhǎng)度的93.9%、91.8%；在862條BESs中，分別有702條(81.40%)定位在白菜A9染色體上，646條(74.9%)定位在甘藍(lán)型油菜A9染色體上，表明芥菜型油菜A9染色體同白菜A9染色體同源性高于甘藍(lán)型油菜A9染色體；芥菜型油菜A9染色體同甘藍(lán)型油菜A9染色體，相對(duì)白菜A9染色體6413031～7827177、14856710～18458808的位置發(fā)生了兩處相同的分別長(zhǎng)為1.4、3.6 Mb的倒位，白菜和芥菜型油菜A9染色體相對(duì)于甘藍(lán)型油菜15676847～16064302的區(qū)域存在長(zhǎng)為0.98 Mb的插入片段，該區(qū)域位于前述的倒位區(qū)域，而且與預(yù)測(cè)但未組裝出的甘藍(lán)型油菜著絲粒位置一致。

芥菜型油菜；A9染色體；物理圖譜；比較基因組

芥菜型油菜屬于十字花科蕓薹屬中芥菜的油用類型，是由白菜和黑芥通過(guò)雜交和染色體自然加倍形成的異源四倍體(AABB，n=2x=18，1068 Mb)[1]。在其進(jìn)化過(guò)程中經(jīng)歷了多次全基因組倍增事件(α，β，γ，bU復(fù)制事件)[2，3]，產(chǎn)生了大量重復(fù)序列。多倍性和重復(fù)序列使基因組增大且更復(fù)雜，從而增加了測(cè)序成本，降低了測(cè)序的準(zhǔn)確性[4]。目前，獲得基因組序列主要有兩種方法[5]：全基因組鳥(niǎo)槍法測(cè)序策略(WGS)和基于物理圖的克隆連克隆測(cè)序策略(CBCS)。新一代測(cè)序技術(shù)以高通量、時(shí)間少、低成本[6]受到科學(xué)家們的關(guān)注和青睞，但其產(chǎn)生的大量短序列不僅加大組裝過(guò)程中計(jì)算機(jī)的運(yùn)行量，而且對(duì)測(cè)序錯(cuò)誤率敏感，在基因組重復(fù)區(qū)甚至出現(xiàn)組裝錯(cuò)誤或序列丟失的現(xiàn)象[7]。CBCS選擇能覆蓋整個(gè)基因組的最低數(shù)目BAC進(jìn)行測(cè)序，一個(gè)克隆長(zhǎng)約126 kb，可以跨越重復(fù)區(qū)，從而減少空隙(GAP)數(shù)，提高了組裝質(zhì)量。完整、高質(zhì)量的基因組參考序列是進(jìn)行基因預(yù)測(cè)、基因組功能注釋、進(jìn)化及比較基因組學(xué)研究的基礎(chǔ)，而克隆連克隆測(cè)序策略是獲得精細(xì)基因組序列的前提[8，9]。 BAC重疊群的構(gòu)建主要有酶切指紋分析、分子標(biāo)記篩選、BAC末端測(cè)序[10]等方法。酶切指紋分析利用HICF[11](高信息量酶切指紋技術(shù))根據(jù)相同的指紋圖譜確定BAC之間的重疊關(guān)系，再通過(guò)指紋重疊群(FPC)或線性拓?fù)渲丿B群(LTC)兩種軟件修飾[12]組裝成重疊群。該方法可用于未知基因組測(cè)序。分子標(biāo)記法包括核酸雜交[13]和PCR篩選[14]。兩者都需要參考已知序列設(shè)計(jì)探針或引物對(duì)BAC文庫(kù)進(jìn)行篩選確定陽(yáng)性克隆之間的重疊關(guān)系。該方法可以實(shí)現(xiàn)遺傳圖譜和物理圖譜的整合。BAC末端測(cè)序則通過(guò)把末端序列與已知近緣生物基因組比對(duì)，從而把BAC錨定到染色體上[15]，重疊群側(cè)翼的BAC末端序列還可以用來(lái)設(shè)計(jì)探針或引物進(jìn)行染色體步移篩選，填補(bǔ)重疊群之間的空隙。完整的物理圖譜往往需要幾種方法綜合使用才能完成。 在蕓薹屬植物A、B、C三個(gè)基因組中，A、C基因組高度同源[16]，其中A9染色體與C8、C9染色體部分同源[17]。根據(jù)已測(cè)序的白菜[18]和甘藍(lán)型油菜[19]基因組數(shù)據(jù)顯示，在A基因組的10條染色體中，A9染色體最長(zhǎng)，白菜A9染色體組裝長(zhǎng)度為38.8 Mb，而甘藍(lán)型油菜A9染色體組裝長(zhǎng)度為33.8 Mb，另有4.1 Mb序列組裝方向未知、不均勻地分布在甘藍(lán)型油菜A9染色體上。A9染色體是易于發(fā)生重排的染色體，著絲粒區(qū)域是易于發(fā)生重排的地方[20]。Cheng等[21]根據(jù)著絲粒特異序列認(rèn)為白菜A9染色體的著絲粒位于10.4～11.3 Mb的位置，而在19.0～19.1 Mb也存在失活的著絲粒痕跡。Mason等[22]采用半四分體(half-tetrad)分析方法，認(rèn)為甘藍(lán)型油菜A9染色體的著絲粒位于15.6～15.9 Mb的區(qū)域，對(duì)應(yīng)白菜A9染色體14.8～15.4 Mb的區(qū)域。A9染色體攜帶許多控制重要性狀的基因如種子大小、種皮顏色、油脂合成、硫苷合成[23～26]等相關(guān)的基因，因而構(gòu)建A9染色體物理圖不僅具有理論意義，更具實(shí)用價(jià)值。

本研究利用定位在A9染色體的標(biāo)記對(duì)構(gòu)建的ZBjuH BAC文庫(kù)，進(jìn)行三維混合PCR 步移篩選，構(gòu)建A9染色體BAC物理圖譜，依據(jù)最小瓦片路徑(Minimal Tiling Path，MTP)原則，挑取盡可能覆蓋全基因組的BAC作為種子BAC，進(jìn)行雙末端測(cè)序，將測(cè)序成功的BES與已測(cè)序的白菜和甘藍(lán)型油菜A9染色體序列進(jìn)行BLAST，分析蕓薹屬植物A9染色體的結(jié)構(gòu)變異。

1　材料與方法

1.1芥菜型油菜BAC文庫(kù)的構(gòu)建

提取芥菜型油菜作圖群體親本紫葉芥S8自交系葉片總DNA，構(gòu)建芥菜型油菜紫葉芥BAC文庫(kù)(187個(gè)384孔板)。構(gòu)建步驟參照文獻(xiàn)[27]操作，由華中農(nóng)業(yè)大學(xué)羅美中教授實(shí)驗(yàn)室協(xié)助完成。

1.2BAC質(zhì)粒池的混池方法

BAC質(zhì)粒DNA的提取參照《分子克隆實(shí)驗(yàn)指南》(薩姆布魯克等，2005)進(jìn)行。為簡(jiǎn)化PCR篩選過(guò)程，本實(shí)驗(yàn)室構(gòu)建了3個(gè)水平的BAC質(zhì)粒DNA池(一級(jí)池、二級(jí)池和三級(jí)池)。一級(jí)池：由每相鄰10塊384孔板的BAC克隆質(zhì)粒DNA混和成一管，共19管。二級(jí)池：由每塊384孔板的BAC克隆質(zhì)粒DNA混和成一管，共187管。三級(jí)池分為橫向三級(jí)池和縱向三級(jí)池。橫向三級(jí)池由一個(gè)384孔板中的每相鄰2行共48個(gè)BAC克隆質(zhì)粒DNA(如002號(hào)平板的002-A/B、002-C/D、……、002-O/P)混和而成，共187×8管；縱向三級(jí)池由384孔板中的連續(xù)2列共32個(gè)BAC克隆質(zhì)粒DNA，(如004號(hào)平板的002-1/2、002-3/4、……、002-23/24)混合而成，共187×12管。在所構(gòu)建的一級(jí)池、二級(jí)池、橫向三級(jí)池、縱向三級(jí)池中都包含了整個(gè)BAC文庫(kù)。

1.3引物的設(shè)計(jì)

共設(shè)計(jì)了引物1728對(duì)，其中參考白菜基因組[18]序列設(shè)計(jì)STS引物311對(duì)，參考芥菜型油菜四川黃籽種皮的非冗余轉(zhuǎn)錄本(unigene)序列設(shè)計(jì)STS引物107對(duì)，參考芥菜型油菜Survey序列設(shè)計(jì)STS引物126對(duì)，利用BAC末端序列設(shè)計(jì)引物1184對(duì)。引物設(shè)計(jì)采用Primer premier 3或Primer premier 5軟件。設(shè)計(jì)引物時(shí)設(shè)定的標(biāo)準(zhǔn)為：產(chǎn)物大小200～500 bp，引物長(zhǎng)度控制在20～23個(gè)堿基，正反引物Tm值60℃左右，GC含量盡量大于50%。

1.4BAC單克隆的篩選

(1)用程序降落PCR篩選BAC文庫(kù)。95℃預(yù)變性5 min；94℃預(yù)變性30 s，62℃退火30 s，72℃延伸30 s，3個(gè)循環(huán)；94℃預(yù)變性30 s，59℃退火30 s，72℃延伸30 s，3個(gè)循環(huán)；94℃預(yù)變性30 s，56℃退火30 s，72℃延伸30 s，28個(gè)循環(huán)；72℃延伸10 min，22℃保存。退火溫度根據(jù)引物進(jìn)行適當(dāng)調(diào)整。

(2)BAC文庫(kù)篩選步驟。對(duì)19個(gè)一級(jí)池進(jìn)行篩選，確定陽(yáng)性克隆位于哪組相鄰的10塊384孔板中；篩選陽(yáng)性一級(jí)池包含的所有二級(jí)池，確定陽(yáng)性克隆位于哪塊384孔板中；篩選陽(yáng)性二級(jí)池對(duì)應(yīng)的橫向和縱向三級(jí)池，把陽(yáng)性克隆定位于384孔板中橫向和縱向交叉所包含的4個(gè)單克?。粚?duì)4個(gè)單克隆進(jìn)行篩選，得到單個(gè)陽(yáng)性克隆。PCR產(chǎn)物采用1.5%瓊脂糖凝膠電泳檢測(cè)。

1.5BAC末端測(cè)序和序列分析

BAC末端序列由上海生物工程技術(shù)公司采用Sanger方法進(jìn)行測(cè)序，測(cè)序引物為M13R和S2。利用chromas軟件分析BES測(cè)序峰圖以判斷測(cè)序質(zhì)量，利用BLASTN程序與BRAD、NCBI上白菜、甘藍(lán)型油菜基因組數(shù)據(jù)等進(jìn)行序列比對(duì)分析，利用比對(duì)聯(lián)配長(zhǎng)度和E值等數(shù)據(jù)初步排除假陽(yáng)性BAC，推斷BAC在染色體序列假分子上的大概位置，并設(shè)計(jì)末端引物用于驗(yàn)證和延伸。 挑選被最多標(biāo)記篩選到、代表整個(gè)物理圖譜的BAC作為種子BAC，利用BLAST軟件將種子BAC的兩個(gè)BES分別比對(duì)到甘藍(lán)型油菜和白菜A9染色體假分子上(E<1e-10)，對(duì)芥菜型油菜A9染色體BAC物理圖譜與白菜、甘藍(lán)型油菜A9染色體進(jìn)行共線性分析。

2　結(jié)果與分析

2.1芥菜型油菜A9染色體物理圖

利用PCR方法進(jìn)行染色體步移篩選，得到陽(yáng)性BAC共3200個(gè)，將其中2070個(gè)BAC進(jìn)行末端測(cè)序，獲得BAC末端序列(BES)3410條，構(gòu)建了由16個(gè)重疊群、5個(gè)單拷貝(BAC數(shù)低于5個(gè))組成的芥菜型油菜A9染色體物理圖譜。按照最小瓦片路徑原則挑選了代表整個(gè)物理圖的種子BAC 435個(gè)，并送往生工進(jìn)行末端測(cè)序，共獲得862條BES用于序列比對(duì)分析(表1)。按照每個(gè)BAC的平均插入片段為126 kb估算，所有重疊群累積共52.08 M，芥菜型油菜A9約為38 Mb，可知該物理圖譜覆蓋A9染色體的1.37倍，因此估計(jì)每個(gè)BAC延伸長(zhǎng)度約為87.6 kb。其中最長(zhǎng)的重疊群包含62個(gè)BACs，長(zhǎng)約5.4 Mb，覆蓋A9染色體的14.3%。N50值(達(dá)到總長(zhǎng)度一半時(shí)，最小的contig長(zhǎng)度)為2.8 Mb，包含32個(gè)BACs。Contig包含的BACs越多表明組裝越準(zhǔn)確。排除BAC數(shù)小于5的5個(gè)單拷貝，重疊群的平均覆蓋長(zhǎng)度約為1.68 Mb。

表1　芥菜型油菜A9染色體物理圖及比對(duì)到白菜、甘藍(lán)型A9染色體的區(qū)域Table 1　BAC physical map of Aj9 and positions mapped to Ar9、An9

注：Ar9表示白菜A9染色體，Aj9表示芥菜型油菜A9染色體，An9表示甘藍(lán)型油菜A9染色體。BAC數(shù)少于5個(gè)的被定義為單拷貝。N表示未比對(duì)到A9上，—表示錨定到A9上的區(qū)域長(zhǎng)度未知。

2.2蕓薹屬植物A9染色體結(jié)構(gòu)比較

甘藍(lán)型油菜和芥菜型油菜的A基因組都來(lái)源于白菜[1]，但它們的基因組組成不同，前者為AACC，后者為AABB。A、C基因組親緣關(guān)系近，但與B基因組親緣關(guān)系遠(yuǎn)[16]。有研究表明，種間雜種的基因型會(huì)影響減數(shù)分裂時(shí)的染色體配對(duì)聯(lián)會(huì)及以后的交換重組[28]。不同的基因組遺傳背景下，芥菜型油菜和甘藍(lán)型油菜是否發(fā)生同樣的結(jié)構(gòu)變化有待研究。本研究利用芥菜型油菜A9染色體的種子BAC末端序列，與白菜和甘藍(lán)型油菜的假分子序列進(jìn)行比較，研究A9染色體在進(jìn)化過(guò)程中的結(jié)構(gòu)變異。

2.2.1A9染色體宏觀共線性比較

將種子BES與白菜scaffold序列[18]進(jìn)行比對(duì)，芥菜型油菜A9物理圖對(duì)應(yīng)白菜基因組上50個(gè)scaffold。Li等[29]利用蕓薹屬60 K Infinium SNP芯片對(duì)472份甘藍(lán)型油菜種質(zhì)進(jìn)行關(guān)聯(lián)分析，對(duì)甘藍(lán)型油菜各染色體序列支架的排列和方向也進(jìn)行了整理。將3種油菜A9染色體結(jié)構(gòu)進(jìn)行比較(圖1，封三)可知：相較于白菜，同是異源四倍體的芥菜型油菜和甘藍(lán)型油菜A9染色體上錨定了之前白菜沒(méi)有進(jìn)行定位的多個(gè)序列支架，如BraScf178、BraScf321、BraScf17_front、BraScf344、BraScf185；芥菜型油菜和甘藍(lán)型油菜又有其特異的支架，芥菜型油菜包含特異的BraScf960、BraScf70、BraScf282、BraScf460、BraScf496及插入BraScf006內(nèi)的BraScf174，甘藍(lán)型油菜包含特異的BraScf311、BraScf103、BraSc120、BraSc106，白菜包含特異的BraScf393；BraScf481在甘藍(lán)型油菜中未找到，在芥菜型油菜中發(fā)生轉(zhuǎn)座；芥菜型油菜中BraScf159和BraScf113的位置與白菜一致，與甘藍(lán)型油菜發(fā)生顛倒；BraScf66～BraScf134與甘藍(lán)型油菜一致，與白菜發(fā)生顛倒。

2.2.2A9染色體微觀共線性比較

用BLAST將862條種子BAC末端序列(BES)定位到白菜和甘藍(lán)型油菜A9染色體假分子上，分別覆蓋了36.4、31.0 Mb的區(qū)域，分別占假分子組裝長(zhǎng)度的93.9%、91.8%。在這862條BESs中，分別有702條(81.40%)定位在白菜A9染色體上、646條(74.9%)定位在甘藍(lán)型油菜A9染色體上，表明芥菜型油菜A9染色體與白菜A9染色體的同源性高于甘藍(lán)型油菜A9染色體。甘藍(lán)型油菜基因組測(cè)序時(shí)，75個(gè)總長(zhǎng)33.8 Mb的scaffolds通過(guò)遺傳標(biāo)記定向錨定在A9染色體上，另有163個(gè)總長(zhǎng)4.1 Mb的scaffolds不能通過(guò)遺傳標(biāo)記定向錨定到A9染色體因而被歸為“chrA9_random”；“chrA9_random”上的scaffolds不均勻地分布在甘藍(lán)型油菜A9染色體上[19]。上述646條定位到甘藍(lán)型油菜A9染色體的BES中，591條比對(duì)在A9染色體上，55條定位在 “chrA9_random”上(圖2紅色區(qū)域，封三)；在未定位到甘藍(lán)型油菜A9上的224條BESs中，有47條BESs定位在“chrAnn_random”上(圖2中綠色區(qū)域，封三)。與白菜A9染色體假分子[18]相比，BAC重疊群中，075M22R ～159A17L對(duì)應(yīng)白菜A9染色體6413031～7827177區(qū)域有約1.4 Mb大小的倒位，對(duì)應(yīng)白菜scaffold57的末端。185J13L ～176B07R之間對(duì)應(yīng)白菜A9染色體14888415 ～18347321區(qū)域有約3.6 Mb倒位，對(duì)應(yīng)白菜scaffold134～scaffold66。而與甘藍(lán)型油菜A9染色體假分子[19]相比，BAC重疊群中，094O07R～108C15L之間的BES均比對(duì)到chrAnn_random上，對(duì)應(yīng)甘藍(lán)型油菜15676847～15769286的區(qū)域，該區(qū)域與預(yù)測(cè)的甘藍(lán)型油菜著絲粒位置一致。因而，推測(cè)該缺失區(qū)域很可能是由于著絲粒區(qū)域序列高度重復(fù)或缺乏可利用的遺傳標(biāo)記而無(wú)法組裝到A9染色體上造成的。

3　討論

3.1BAC物理圖譜構(gòu)建中的問(wèn)題

本研究使用A9特異性標(biāo)記及BAC末端序列標(biāo)記步移篩選法，構(gòu)建了芥菜型油菜A9染色體BAC物理圖譜。在物理圖譜的構(gòu)建過(guò)程中，BAC文庫(kù)的篩選是最關(guān)鍵也是最耗時(shí)的步驟，三級(jí)DNA池的逐步篩選法大量減少了文庫(kù)篩選過(guò)程，而且每個(gè)陽(yáng)性克隆都經(jīng)過(guò)了4次驗(yàn)證，結(jié)果可靠。但是，文庫(kù)篩選過(guò)程中由于基因組中的重復(fù)序列仍存在假陽(yáng)性現(xiàn)象，因而在引物設(shè)計(jì)過(guò)程中可以與白菜或甘藍(lán)型油菜A9假分子序列比對(duì)，使引物唯一比對(duì)在預(yù)測(cè)目的區(qū)域。另外可以增加引物長(zhǎng)度、提高退火溫度、采用PCR產(chǎn)物驗(yàn)證等方法減少假陽(yáng)性出現(xiàn)的概率。目前所構(gòu)建的芥菜型油菜A9染色體物理圖仍然存在一些空隙，且多對(duì)應(yīng)白菜scaffold的末端。scaffold的末端一般是因?yàn)槲挥谥貜?fù)序列而無(wú)法將相鄰scaffold連接起來(lái)。 因此空隙存在的地方很可能是重復(fù)區(qū)域，但也可能是由于HindIII內(nèi)切酶的偏好型，導(dǎo)致某些區(qū)域難以形成126 kb左右的片段，而導(dǎo)致BAC文庫(kù)中沒(méi)有代表該區(qū)域的陽(yáng)性克隆。因此，在BAC文庫(kù)構(gòu)建過(guò)程中通常需要采用多酶切的方法。

3.2蕓薹屬植物A9染色體比較分析

通過(guò)芥菜型油菜種子BAC末端序列與白菜、甘藍(lán)型油菜基因組序列比對(duì)，發(fā)現(xiàn)芥菜型油菜和甘藍(lán)型油菜A9染色體在異源四倍體的形成過(guò)程中發(fā)生了兩處相同的染色體倒位。預(yù)測(cè)的甘藍(lán)型油菜A9染色體著絲粒跨越其中一處倒位，同Parkin等[20]認(rèn)為的A9染色體是易于發(fā)生重排的染色體，著絲粒是易于發(fā)生重排的區(qū)域結(jié)論一致。芥菜型油菜A9染色體物理圖和白菜基因組包含了預(yù)測(cè)的甘藍(lán)型油菜著絲粒區(qū)域，表明針對(duì)大而復(fù)雜基因組測(cè)序時(shí)，基于物理圖的克隆連克隆測(cè)序策略優(yōu)于全基因組鳥(niǎo)槍法測(cè)序策略。芥菜型油菜A9染色體與白菜A9染色體的同源性高于甘藍(lán)型油菜A9染色體。甘藍(lán)型油菜全基因組序列顯示A、C亞基因組之間，特別是A1與C1，A2與C2，A9與C9之間的同源交換非常頻繁[19]，因此推測(cè)這有可能是因?yàn)楦仕{(lán)型油菜在進(jìn)化過(guò)程中，A基因組與高度同源的C基因組之間的同源重組使甘藍(lán)型油菜A9產(chǎn)生了異于白菜和芥菜型油菜的變異，也可能是甘藍(lán)型油菜組裝時(shí)，遺傳標(biāo)記數(shù)目有限導(dǎo)致scaffolds無(wú)法特異定位在A9上。

4　結(jié)論

本研究構(gòu)建了由16個(gè)重疊群、5個(gè)單拷貝組成的芥菜型油菜A9染色體物理圖，選擇其中435個(gè)BAC作為A9染色體測(cè)序的候選克隆，通過(guò)BAC末端序列與Ar9、An9假分子比對(duì)，發(fā)現(xiàn)Aj9與Ar9同源性高于An9；Aj9同An9，相對(duì)Ar9發(fā)生了兩處相同的分別長(zhǎng)為1.4、3.6 Mb的倒位；Aj9物理圖覆蓋了An9未組裝出的著絲粒區(qū)域。

[1]UN.Genome analysis inBrassicawith special reference to the experimental formation ofB.napusand peculiar mode of fertilization[J].Japanese Journal of Botany，1935，7：389-452.

[2]Bowers JE，Chapman BA，Rong J，et al.Unravelling angiosperm genome evolution by phylogenetic analysis of chromosomal duplication events[J].Nature，2003，422：433-438.

[3]Bancroft I，Bancroft I，Botanik HC.Genetics and Genomics of theBrassicaceae[M].New York：Springer，2011.

[4]Jackson SA，Iwata A，Lee SH，et al.Sequencing crop genomes：approaches and applications[J].New Phytologist，2011，191：915-925.

[5]Catherine F，Leach JE，Jane R，et al.Crop genome sequencing：lessons and rationales[J].Trends in Plant Science，2011，16(2)：77-88.

[6]Shendure J，Ji H.Next-generation DNA sequencing[J].Nature Biotechnology，2008，26(10)：1135-1145.

[7]Paszkiewicz K，Studholme DJ.De novo assembly of short sequence reads[J].Briefings in Bioinformatics，2010，11(5)：457-472.

[8]Vivien M.Next-generation sequencing：The genome jigsaw[J].Nature，2013，501：263-268.

[9]Pan Y，Deng Y，Lin H，et al.Comparative BAC-based physical mapping ofOryzasativassp.indicavar.93-11 and evaluation of the two rice reference sequence assemblies[J].Plant Journal，2014，77(5)：795-805.

[10]湯飛宇，張?zhí)煺?重疊群物理圖譜的構(gòu)建及其應(yīng)用[J].基因組學(xué)與應(yīng)用生物學(xué)，2009，28(1)：195-201.

[11]Luo MC，Thomas C，You FM，et al.High-throughput fingerprinting of bacterial artificial chromosomes using the snapshot labeling kit and sizing of restriction fragments by capillary electrophoresis[J].Genomics，2003，82(3)：378-389.

[12]Barabaschi D，Magni F，Volante A，et al.Physical mapping of bread wheat chromosome 5A：An integrated approach[J].Plant Genome，2015，8(3)：1-24.

[13]Mozo T，Dewar K，Dunn P，et al.A complete BAC-based physical map of theArabidopsisthalianagenome[J].Nature Genetics，1999，22(3)：271-275.

[14]Vu GT，Caligari PD，Wilkinson MJ.A simple，high throughput method to locate single copy sequences from Bacterial Artificial Chromosome (BAC) libraries using High Resolution Melt analysis[J].BMC Genomics，2010，11(10)：301.

[15]Mun J，Kwon S，Seol Y，et al.Sequence and structure ofBrassicarapachromosome A3[J].Genome Biology，2010，11(9)：79-82.

[16]栗茂騰，張椿雨，劉列釗，等.蕓薹屬A，B和C基因組之間關(guān)系研究進(jìn)展[J].遺傳，2005，27(4)：671-676.

[17]Panjabi P，Jagannath A，Bisht NC，et al.Comparative mapping ofBrassicajunceaandArabidopsisthalianausing Intron Polymorphism (IP) markers：homoeologous relationships，diversification and evolution of the A，B and CBrassicagenomes[J].BMC Genomics，2008，9(4)：113.

[18]Wang X，Wang H，Wang J，et al.The genome of the mesopolyploid crop speciesBrassicarapa[J].Nature Genetics，2011，43(10)：1035-1039.

[19]Boulos C，F(xiàn)rance D，Shengyi L，et al.Plant genetics.Early allopolyploid evolution in the post-NeolithicBrassicanapusoilseed genome[J].Science，2014，345：950-953.

[20]Parkin I.Chasing ghosts：comparative mapping in theBrassicaceae[J].Plant Genetics & Genomics Crops & Models，2010，9：153-170.

[21]Cheng F，Mandáková T，Wu J，et al.Deciphering the diploid ancestral genome of the mesohexaploidBrassicarapa[J].Plant Cell，2013，25(5)：1541-1554.

[22]Mason AS，Rousseau-Gueutin M，Morice J，et al.Centromere locations inBrassicaA and C genomes revealed through half-tetrad analysis[J].Genetics，2015，115：183-210.

[23]劉顯軍，袁謀志，官春云，等.芥菜型油菜黃籽性狀的遺傳、基因定位和起源探討[J].作物學(xué)報(bào)，2009，35(5)：839-847.

[24]Ji F，Yan L，Lei S，et al.Characterization of metabolite quantitative trait loci and metabolic networks that control glucosinolate concentration in the seeds and leaves ofBrassicanapus[J].New Phytologist，2012，193(1)：96-108.

[25]Hui W，Jian W，Silong S，et al.Glucosinolate biosynthetic genes inBrassicarapa[J].Gene，2011，487(2)：135-142.

[26]Zhao J，Huang J，Chen F，et al.Molecular mapping ofArabidopsisthalianalipid-related orthologous genes inBrassicanapus[J].Theoretical & Applied Genetics，2012，124(2)：407-421.

[27]Luo M，Wing RA.An improved method for plant BAC library construction[J].Methods in Molecular Biology，2003，236：3-20.

[28]Mason AS，Huteau V，Eber F，et al.Genome structure affects the rate of autosyndesis and allosyndesis in AABC，BBAC，CCAB[J].Chromosome Research，2010，18(6)：655-666.

[29]Li F，Chen B，Xu K，et al.Genome-wide association study dissects the genetic architecture of seed weight and seed quality in rapeseed (BrassicanapusL.)[J].DNA Research，2014，21(4)：1-13.

Construction of BAC Physical Map on Chromosome A9 in Brassica juncea and Comparative Genome Analysis

LIU Fangying，LU Ying，WANG Zhuo，LIU Xianjun，HU Xuefang，LIU Xudong，XU Haipeng，LIU Zhongsong*

(Oil Crops Research Institute，Hunan Agricultural University，Changsha，Hunan 410128，China)

To abtain a precise physical map ofBrassicajunceachromosome A9，the markers located on the A9 chromosome were used in this study to screenBrassicajunceaBAC library by three-dimensional mixed PCR screening.A total of 3200 positive BAC were screened out；constructing a physical map composed of 16 contigs and 4 singletons.According to the principle of the minimum path，435 BACs were chosen as seed BACs to represent the whole physical map and were subject to end sequencing，at last，862 qualified sequences were received.As the average insert size of BAC is 126 kb，the contigs accounted up to 52.08 M，covering 1.37 times of chromosome A9.By BLAST，the seed BAC end sequences were mapped to chromosome A9 pseudomolecule ofBrassicarapaandBrassicanapusrespectively，results showed that the physical map covered 36.4 Mb and 31.0 Mb，equivalent to 93.9%、91.8% of the whole genome assemble length.Among the 862 BAC end sequences，702 sequences (81.4%) were located onBrassicarapaA9 chromosome and 646 sequences(74.9%)toBrassicanapusA9 chromosome，which meansBrassicajunceaA9 chromosome have higher homology withBrassicarapaA9 chromosome thanBrassicanapusA9 chromosome.Compared withBrassicanapusA9 pseudomolecule and BAC end sequence ofBrassicajunceaA9 chromosome，there are two chromosome inversions about 1.4 Mb、3.6 Mb in length at the position of 6413031-7827177 and 14856710-18458808 inBrassicarapaA9.While compared with BAC end sequence ofBrassicajunceaA9，there is a inserted mutation about 0.98 Mb in length at the position of 15676847-16064302 ofBrassicanapusA9 pseudomolecule，which lie in one of the two chromosome inversions and was consist with the predicted position of centromere ofBrassicanapusA9 chromosome remained to assembled.

Brassicajuncea；chromosome A9；physical map；compare genomics

2016-03-07

劉芳瑛(1991-)，女，碩士研究生，Email：365698105@qq.com。*通信作者：劉忠松，博士，教授，Email：zsliu48@sohu.com。

國(guó)家自然科學(xué)基金(30971799)。

S565.403.2

A

1001-5280(2016)04-0347-07

10.16848/j.cnki.issn.1001-5280.2016.04.01