小麥祖先種A 組/D 組與水稻基因組比較分析

許陶瑜，王長彪，韓 斌，趙興華，劉 江，任永康，牛瑜琦，唐朝暉

（1.山西大學生物工程學院，山西太原030006；2.山西省農業科學院經濟作物研究所，山西汾陽032200；3.山西省農業科學院生物技術研究中心，山西太原030031；4.山西省農業科學院，山西太原030031）

小麥祖先種烏拉爾圖小麥、粗山羊草是異源六倍體小麥基因組A 和D 的供體，是構成小麥基因的基礎。水稻（Oryza sativa）是較早進行基因測序的物種，它具有較小的基因組、高效率的遺傳轉化、高密度的遺傳圖譜和物理圖譜等優勢，使它成為了禾本科植物研究的模式作物。利用水稻作為參考基因組與小麥祖先種進行全基因組比較分析，可以更深入地了解小麥祖先種的遺傳基礎、分子結構、進化關系，為小麥育種奠定基礎。

隨著基因組測序技術的發展[1-2]，自SANGER 發明雙脫氧測序進行核苷酸序列測定以來，越來越多的物種基因完成了全基因組測序。隨著對水稻[3]、小麥[4-5]及小麥祖先種烏拉爾圖[6-7]及粗山羊草[8]基因組測序工作的完成，使得水稻和小麥祖先種基因組比對分析成為可能?；蚪M比對分析是在基因組測序基礎上，通過對核苷酸序列、氨基酸序列進行全基因組比對，從全基因組水平認識基因組之間的異同[9-10]。比較基因組研究通常有2 種方式，一種是通過比較解釋差異，另一種是通過比較揭示共性。

培育高產[11]、抗病、抗蟲、抗逆[12]、高水肥利用效率的的優質小麥、水稻品種一直是育種家的目標。過去田間育種，主要通過田間表型觀察、性狀分析來進行選擇育種，往往需要多代，選育周期長，過程復雜，工作量大。隨著基因組測序的完成，育種家將基因組和表型觀察聯系起來，可以通過基因鑒定來發掘不易鑒定的表型，將基因分析和遺傳分析聯系起來，大大地縮短了育種年限，節約了成本，提高了育種效率。通過小麥祖先種與水稻基因組比對分析，可以挖掘祖先種優異抗性基因輔助小麥分子育種。

為深入研究小麥祖先種與水稻的遺傳基礎、分子結構、進化關系，了解3 個物種同源核苷酸序列和基因序列的排列與分布[13]，本試驗構建了3 個物種的核苷酸序列和基因序列庫，利用Blast 方法篩選同源核苷酸序列和同源基因序列，然后用circos，synmap 4.2 分別對同源核苷酸和同源基因序列作圖，旨在揭示3 個物種間分子水平的共線性關系[14-15]及進化關系[16-18]。

1 材料和方法

1.1 材料

烏拉爾圖小麥序列從https://www.ncbi.nlm.nih.gov/nuccore/MKGO00000000.1 下載，數據更新公布的時間為2018 年5 月。

粗山羊草序列從https://www.ncbi.nlm.nih.gov/n uccore/MCGU00000000.1 下載。數據更新公布的時間為2017 年11 月。

水稻基因組全基因組序列從https://www.ncbi.nlm.nih.gov/genome/10 下載，數據更新時間為2018 年12 月。

1.2 方法

1.2.1 核苷酸序列全局比對 將烏拉爾圖小麥、粗山羊草作為目標基因組，水稻作為參考基因組，對核苷酸序列進行Blast N（e-value 為1e-20）比對，對結果篩選分析；使用circos 對同源核苷酸序列進行關聯作圖。

1.2.2 核苷酸序列局部比對 選取水稻1 號染色體和烏拉爾圖小麥3 號染色體的同源核苷酸序列，使用軟件mauve 對同源核苷酸序列做線性比對。

1.2.3 基因比對 構建3 個物種編碼蛋白質的基因數據庫，對它們的氨基酸序列進行Blast X（e-value為1e-10）比對，使用symap 4.2 對數據結果作圖。

2 結果與分析

2.1 核苷酸序列全局比對

結果表明，烏拉爾圖小麥中有6 786 546 條與水稻同源核苷酸序列，粗山羊草中有748 050 條與水稻同源核苷酸序列，小麥祖先種中存在大量水稻同源序列。篩選長度為4 kb 以上的同源核苷酸序列使用circos 軟件作圖（圖1），烏拉爾圖小麥篩選到2 007 條大片段（4 kb 以上）與水稻同源序列，而粗山羊草中篩選到6 252 條大片段序列。這說明粗山羊草在DNA 進化中更保守，與水稻的親緣關系近。在序列比對中，本研究發現，烏拉爾圖小麥、粗山羊草中的某條核苷酸序列可能在水稻的多條染色體中找到同源序列，水稻中的某條同源核苷酸序列也會匹配到烏拉爾圖小麥或粗山羊草的多條染色體上，這表明3 個物種均存在核苷酸序列的復制現象。

2.2 核苷酸序列局部比對

試驗選取水稻1 號染色體和烏拉爾圖小麥3 號染色體的Blast 結果，使用mauve 軟件對同源核苷酸序列作線性比較。結果表明，水稻1 號染色體和烏拉爾圖小麥3 號染色體存在大量同源核苷酸序列。圖中使用相同顏色表示了相同的同源核苷酸序列區塊，水稻和烏拉爾圖標識同源序列的顏色區塊一致性排列，這表明它們之間存在同源核苷酸序列宏觀的同線性和共線性（圖2）。

結果還表明，同源核苷酸序列主要分布在水稻正義鏈上，而烏拉爾圖小麥中同源序列正義鏈、反義鏈均有分布。在1 號染色體上的水稻某一同源核苷酸序列和烏拉爾圖小麥3 號染色體的多個不同位置的核苷酸序列存在同源關系。

2.3 基因比對分析

氨基酸序列比對中，本研究用烏拉爾圖小麥、粗山羊草的每一條染色體和水稻全基因組進行比較，發現每一條染色體都有80%以上的同源基因匹配到水稻全基因組中。烏拉爾圖小麥中5 號染色體與水稻的基因匹配率最高，為84%。粗山羊草中2 號染色體與水稻的匹配率最高，為88%。在共線性的區塊中發現，約50%的基因共線性區塊存在倒位現象，這表明在烏拉爾圖小麥、粗山羊草和水稻的比較中，染色體的倒位現象是普遍存在的。

從基因比較的共線性排列中可以得出，烏拉爾圖小麥1 號染色體和水稻1，5，10 號染色體存在大片段共線性排列，由水稻5，10 號染色體融合形成。烏拉爾圖小麥2 號染色體和水稻4，7 號染色體有大片段共線性排列。烏拉爾圖小麥3 號染色體和水稻1 號染色體有大片段共線性排列，由1 號染色體進化而來。烏拉爾圖小麥4 號染色體來源于水稻3，11 號染色體。烏拉爾圖小麥5 號染色體來源于水稻3，9，12 號染色體。烏拉爾圖小麥6 號染色體來源于水稻2 號染色體。烏拉爾圖小麥7 號染色體由水稻6，8 號染色體融合形成?；虮容^中還發現，烏拉爾圖小麥4，5 號染色體均和水稻3 號染色體的某一片段存在共線性，而且在烏拉爾圖小麥4 號染色體中發生了倒位（圖3）。在粗山羊草和水稻的基因組比較中也發現了類似的基因共線性排列（圖3，4）。烏拉爾圖小麥、粗山羊草的3 號與6 號染色體基因分別與水稻的1，2 號同源整體，這表明在進化中，這2 條染色體相對保守，沒有發生染色體大片段的移位和融合。

從圖3，4 可以看出，烏拉爾圖小麥和粗山羊草的差別為：和水稻共線性的基因在它們中不是相同排列的。例如，烏拉爾圖小麥7 號染色體和水稻8 號染色體只在短臂上有基因共線性，而在粗山羊草7 號染色體的短臂、長臂均與水稻8 號染色體有基因的共線性。這表明烏拉爾圖小麥和粗山羊草在分化過程中，DNA 分子也發生了選擇性進化。一些DNA 分子由于插入、缺失、轉座子插入等原因而發生了基因的丟失。而一些同源基因的DNA 分子在進化過程中由于變異的積累，形成了具有新功能基因的DNA 序列。

3 結論與討論

水稻和小麥都是禾本科單子葉植物。對小麥祖先種及水稻的基因組比較分析，有助于理解高等植物的進化與發育關系。本試驗表明，烏拉爾圖小麥、粗山羊草和水稻無論是核苷酸序列水平還是基因水平都存在高度的、整體水平的共線性，同時它們的染色體中存在著大量的倒位和復制現象。

本試驗通過對小麥祖先種和水稻全基因組比較[19]，全面地在核苷酸序列水平和基因水平對小麥祖先種和水稻進行了DNA 水平和基因水平的關聯分析，闡述了小麥祖先種染色體的進化過程，為進一步研究小麥遺傳多樣性、基因組分子進化、借鑒水稻基因組研究對小麥基因組進行注釋[20]奠定了基礎。