棉花共有4個(ge) 栽培種，其中經濟價(jia) 值較高的是2個(ge) 四倍體(ti) 栽培種——陸地棉和海島棉。讓陸地棉和海島棉“聯手”，將海島棉的優(you) 質、抗性基因“引入”到高產(chan) 陸地棉中，從(cong) 而培育高產(chan) 優(you) 質多抗的棉花新品種，一直是育種家孜孜以求的目標。

　　河北農(nong) 業(ye) 大學華北作物改良與(yu) 調控國家重點實驗室教授馬峙英團隊率先組裝了陸地棉和海島棉現代品種基因組，破譯了海陸種間、陸陸種內(nei) 基因組結構變異及其規律，揭示了陸地棉大規模種質材料結構變異的遺傳(chuan) 效應，為(wei) 作物重要性狀改良提供了新的理論依據和資源。8月9日，相關(guan) 論文在線發表於(yu) 《自然—遺傳(chuan) 學》。

　　首次組裝陸海現代品種高質量基因組

　　陸地棉產(chan) 量高、適應性強，占全球棉花總產(chan) 量90%以上；海島棉產(chan) 量低、纖維品質優(you) 良、抗病性好，約占全球棉花產(chan) 量10%。長期以來，育種家在獲得多逆境抗性的同時，一直努力通過協同提高產(chan) 量和品質來培育陸地棉新品種，將海島棉的優(you) 異性狀轉育到陸地棉中，以進一步改良品質和抗性。

　　然而，論文第一作者、共同通訊作者馬峙英告訴《中國科學報》，人們(men) 對棉花現代育成品種的基因組信息、現代育種過程中基因組的結構變異仍缺乏深入了解，有關(guan) 棉花結構變異的遺傳(chuan) 效應鮮為(wei) 人知。

　　陸地棉和海島棉都是由二倍體(ti) A基因組和D基因組共同組成的異源四倍體(ti) AADD基因組，是由A基因組和D基因組通過自然雜交和染色體(ti) 加倍而形成。

　　馬峙英介紹，他們(men) 率先組裝了陸、海現代品種高質量基因組，並發現上千個(ge) 棉屬新基因。農(nong) 大棉8號（NDM8）是我國自育陸地棉現代品種，優(you) 質高產(chan) 多抗；Pima90是我國從(cong) Pima棉中係選獲得的海島棉新品係，抗病優(you) 質，長期用於(yu) 棉花抗病遺傳(chuan) 、QTL定位和基因克隆及分子育種。

　　該團隊組裝的農(nong) 大棉8號和Pima90兩(liang) 個(ge) 品種的基因組大小分別為(wei) 2.29Gb和2.21Gb，基因組錨定率分別為(wei) 99.57%和99.75%。他們(men) 在農(nong) 大棉8號中鑒定出80124個(ge) 基因，其中1499個(ge) 為(wei) 預測的新基因；在Pima90中鑒定到79613個(ge) 基因，其中1267個(ge) 為(wei) 預測的新基因。

　　破譯海陸棉基因組結構變異圖譜

　　此前，科學家已經發現了很多基因組中的單核苷酸多態性對棉花性狀產(chan) 生的影響。馬峙英告訴《中國科學報》，單核苷酸多態性是一個(ge) 堿基發生變化，實際上，影響基因表達和功能的還有染色體(ti) 結構變異，後者的影響更大。

　　染色體(ti) 結構變異的類型有堿基的缺失、重複、倒位、易位等，往往是多個(ge) 堿基同時發生變異。“這些變異對基因的轉錄翻譯產(chan) 生了很大影響，對育種家關(guan) 心的經濟性狀變化影響比較大。”馬峙英說，他們(men) 發現了一批影響基因表達的新變異。

　　論文共同通訊作者、河北農(nong) 業(ye) 大學教授王省芬告訴《中國科學報》，為(wei) 了在陸地棉現代育種中潛在有效地利用海島棉基因組變異，該團隊比對分析了Pima90與(yu) 農(nong) 大棉8號基因組。

　　他們(men) 發現，海島棉存在高度的基因組多樣性：插入和缺失在D亞(ya) 組上的密度明顯高於(yu) A亞(ya) 組；31296個(ge) 變異—基因對在海島棉組織中顯著特異表達；5815個(ge) 插入缺失位於(yu) 5256個(ge) 基因的外顯子區，其中3178個(ge) 變異影響了基因轉錄。

　　而且，Pima90中編碼蔗糖合酶的基因存在小片段的缺失。該基因在纖維伸長期和次生壁加厚期高表達，對海島棉纖維長度和強度的形成有重要作用。追蹤利用Pima90育成的新材料農(nong) 大373-9的基因組序列，發現該材料導入了Pima90中編碼蔗糖合酶的基因。此外，該材料還獲得了來自於(yu) Pima90的其他171個(ge) 外顯子區結構變異，表明了海島棉對陸地棉改良的可行性。

　　為(wei) 了探明陸地棉現代品種的基因組變化，他們(men) 又將農(nong) 大棉8號基因組與(yu) 早期品種TM-1基因組進行比較，發現農(nong) 大棉8號存在76568個(ge) 結構變異，其中28626個(ge) 結構變異能夠在10~1081個(ge) 重測序種質材料中檢測到。與(yu) 海島棉類似，農(nong) 大棉8號基因組的插入和缺失在D亞(ya) 組上的密度明顯高於(yu) A亞(ya) 組。然而，陸地棉基因組的插入和缺失在端粒附近具有偏好性，是染色體(ti) 其他區域的3.71倍。

　　此外，編碼肉桂酰輔酶A還原酶的基因在TM-1中存在1個(ge) 缺失，影響了抗黃萎病功能。群體(ti) 樣本分析還發現，現代品種較早期品種獲得了1128個(ge) NDM8型結構變異，反映了現代育種改良進程的選擇效應。

　　探明陸地棉結構變異遺傳(chuan) 效應

　　獲得了海島棉和陸地棉的結構變異圖譜後，他們(men) 又將結構變異對棉花重要經濟性狀的遺傳(chuan) 效應進行了分析。

　　論文共同通訊作者、河北農(nong) 業(ye) 大學教授張豔告訴《中國科學報》，為(wei) 探明基因組結構變異對重要性狀的遺傳(chuan) 效應，他們(men) 對1081份世界各地的陸地棉種質資源深度重測序，以農(nong) 大棉8號為(wei) 參考基因組獲得了304630個(ge) 結構變異，結合大規模多環境評價(jia) 獲得的纖維長度、強度、鈴重、衣分和黃萎病抗性數據分析，發現了446個(ge) 與(yu) 重要性狀顯著關(guan) 聯的結構變異。

　　研究發現，品質性狀的結構變異主要位於(yu) D亞(ya) 組，而產(chan) 量性狀的變異主要位於(yu) A亞(ya) 組。在907個(ge) 與(yu) 纖維品質、產(chan) 量相關(guan) 的基因中，84.23%的基因在纖維不同發育時期表達。

　　纖維長度是棉花最重要的經濟性狀之一。據介紹，纖維長度增加1毫米（mm），每噸皮棉售價(jia) 一般可增加300元左右。在陸地棉Dt11染色體(ti) 上，有69個(ge) 結構變異能夠使纖維長度顯著增加0.71~0.99mm，56個(ge) 結構變異能夠顯著增加1.00~1.19mm。如此可使纖維長度從(cong) 27mm或者28mm級別增至29mm級別。

　　此外，黃萎病是影響棉花產(chan) 量和品質的重要病害。他們(men) 發現，在陸地棉Dt11染色體(ti) 的3個(ge) 結構變異能夠使棉花黃萎病病情指數大大降低，讓棉花從(cong) 感病變為(wei) 耐病。

　　在發現的60個(ge) 抗黃萎病相關(guan) 基因中，編碼S—去甲烏(wu) 藥堿合成酶的基因GhNCS在抗黃萎病品種中的表達量顯著低於(yu) 感病品種，沉默該基因導致抗病性顯著增強：可以使NDM8由耐病變為(wei) 抗病，使棉花品種CCRI8從(cong) 感病變為(wei) 耐病。“這證明GhNCS是控製黃萎病抗性的一個(ge) 重要新基因。”張豔說。

　　　　記者 李晨