棉花陸海漸滲雜交群體纖維品質性狀的QTL定位

馬 君，楊延龍，師維軍，汪鵬龍，鄭巨云，郭仁松，胡文冉，楊 洋

(1.新疆農業科學院經濟作物研究所, 烏魯木齊 830091；2.新疆農業科學院庫車陸地棉試驗站，新疆庫車 842000；3.新疆農業科學院核技術生物技術研究所，烏魯木齊 830091)

0 引 言

【研究意義】陸地棉適應性廣且產量高，而海島棉高抗黃萎病且纖維品質優良，但是產量偏低、生育期長[1]。利用海島棉優良的纖維品質性狀，將其優良基因導入陸地棉，培育出適應性廣、產量高且纖維品質優良的陸海漸滲系，是對陸地棉進行優良性狀改良的有效方法[2]。在陸地棉高產的基礎上對其進行纖維品質的改良也成為當前棉花育種的重要研究方向[3]。棉花的纖維品質是由多基因控制的數量性狀，受環境因素影響較大，普遍存在基因型與環境的互作，且與產量構成因素存在負相關等，遺傳機制較復雜[4]。常規遺傳育種方法在提高陸地棉產量、改良纖維品質的過程中起到了重要作用，培育了大面積推廣的優良品種。種內雜交選育中骨干親本的集中使用和長期的人工選擇會導致陸地棉遺傳基礎狹窄[5-6]。有效發掘利用海島棉優異性狀基因，對拓寬陸地棉栽培種遺傳基礎有重要意義。【前人研究進展】分子標記技術的快速發展為標記輔助選擇和QTL聚合育種提供了條件，可以從目標性狀的分子水平上進行選擇，有利于提高表型性狀選擇的準確性，也是提高作物遺傳改良及育種進程的有效途徑[7]。通過構建連鎖遺傳圖譜，與棉花纖維品質性狀相關的QTL(Quantitative trait locus)可以在染色體組上準確定位，通過分子標記輔助選擇育種減少常規育種選擇的盲目性[2]。目前對陸地棉與海島棉優勢結合進行棉花纖維品質性狀QTL定位做了許多研究，構建了連鎖遺傳圖譜，并檢測到大量與纖維品質相關的QTL位點[2,8-14]。【本研究切入點】優良陸海漸滲系群體是基因聚合、QTL定位等研究的理想材料，需研究利用與纖維品質緊密連鎖的分子標記進行輔助選擇育種。【擬解決的關鍵問題】選用早中熟陸地棉品種新陸中60號為母本與優質抗病海島棉品種新海41號為父本進行雜交，并以新陸中60號為輪回親本構建回交群體BC1F2。利用SSR分子標記，對陸海漸滲雜交BC1F2分離群體進行基因型分型，構建陸海雜交的遺傳連鎖圖譜，種植BC1F2家系群體，測定纖維品質性狀，利用區間作圖法對纖維品質性狀進行QTL定位，檢測與棉花優質纖維品質性狀相關的穩定QTL位點，為海島棉優良基因定位和利用奠定研究基礎。

1 材料與方法

1.1 材 料

2014年在庫車陸地棉試驗站，利用高純系陸地棉品種新陸中60號為母本與優質海島棉品種新海41號為父本雜交，獲得F1種子，2015年再以新陸中60號為輪回親本與F1回交獲得BC1F1群體種子。

2018年在庫車陸地棉試驗站種植新陸中60號、新海41號及BC1F2分離群體，在棉花現蕾期取幼嫩葉片采樣后立即用液氮速凍，親本材料隨機各選5株采樣，BC1F2群體選取190株采樣。獲得151份前期已采樣的BC1F2群體。

1.2 方 法

1.2.1 纖維相關性狀的調查

于2018年10月田間取樣，在BC1F2群體中選取已標記過的151個單株收花，軋花后取15～20 g纖維樣品送農業農村部棉花品質監督測試中心進行品質檢測，檢測指標包括纖維長度、纖維強度、馬克隆值、整齊度和伸長率。

1.2.2 SSR分子標記檢測與PCR 擴增

棉花基因組DNA的提取參照宋國立等[15]改良的CTAB法，其中純化 DNA的操作步驟略作改進，用紫外分光光度計測定DNA濃度。SSR的分析流程參照李立群等[16]方法，優化體系參照許玉蘭等[17]的方法，確定最佳反應體系為15 μL體系，模板DNA(5 ng/μL)1 μL，正反引物各1 μL(濃度為10 μmol/mL)，2×PCR mix 7.5 μL(加拿大應用生物材料有限公司(abm Inc.)的2×PCR mix)，dd H2O 4.5 μL。反應程序：94℃ 5 min；95℃ 15 min；Tm 15 s；72℃15 min；30個循環；72℃10 min；退火溫度(Tm)根據引物進行調節。

1.3 數據處理

統計SSR多態性標記，在BC1F2群體中與母本新陸中60號相同的條帶記為A，與父本新海41號相同的條帶記為B，具有雙親性狀的條帶記為H，BC1F2代中有而雙親中都沒有的條帶記為D，缺失數據記為-。使用Join Map 4.0軟件構建SSR標記遺傳連鎖圖譜，連鎖的最低LOD值為3.0，最大遺傳距離為50 cM，作圖函數采用Kosambi函數。采用QTL分析軟件Windows QTL Cartographer 2.5，選用復合區間作圖(composite interval mapping，CIM)方法，進行纖維相關性狀的QTL定位和效應檢測，QTL命名方法按照 QTL(q)+性狀(英文縮寫)+染色體(chr)+QTL個數(數字)表示。采用單標記分析(Ici Mapping)方法分析纖維品質性狀相關的定位。

2 結果與分析

2.1 遺傳圖譜的構建與纖維品質性狀的QTL位點

研究表明，共獲得105對差異明顯的多態性引物，其中有23對引物不能檢測出BC1F2代的基因分離，多態性頻率為42.17%。構建了一個包含52個標記、14個連鎖群的遺傳圖譜，該圖譜總長824 cM，覆蓋棉花基因組的18.5%，最長的連鎖群為150.3 cM，包含6個標記，最短的為0.3 cM，包含2個標記。圖1

2.2 親本及BC1F2家系纖維品質性狀

研究表明，兩親本間各性狀差異明顯、范圍適中，海陸漸滲系性狀值(纖維長度, 斷裂比強度和伸長率)均高于親本新陸中60號，而馬克隆值和纖維整齊度相對偏低。各項纖維指標在BC1F2家系群體中均呈連續分布，這些性狀都受多基因控制的數量性狀，其中，纖維長度略超親本，存在超親分離。纖維品質各項指標除伸長率外，偏度都小于1，符合典型的數量性狀正態分布特征。表1

圖1 遺傳圖譜與纖維品質性狀的 QTL 位點Fig.1 Linkage map and QTL location of fiber quality traits

表1 親本及BC1F2家系纖維品質性狀Table 1 Analysis of fiber quality in parent and BC1F2 population

2.3 BC1F2群體纖維品質性狀的相關性

研究表明，BC1F2群體纖維品質各性狀間存在著復雜的相關性，纖維長度與纖維整齊度呈極顯著正相關，相關性系數為0.552；纖維長度與馬克隆值呈極顯著負相關，相關性系數為-0.534；斷裂比強度與纖維整齊度和伸長率呈極顯著正相關，相關性系數分別為0.872與0.983；纖維整齊度和伸長率呈極顯著正相關，相關性系數為0.908。存在正相關的各纖維品質之間可相互促進提高，存在負相關的各纖維品質之間相互制約不可能同步提高。表2

2.4 纖維品質性狀的QTL定位

研究表明，檢測到與纖維品質性狀相關的QTL位點1個，該位點與纖維長度相關，位于Chr14的32.31 cM處，標記區間位于BNL1905-BNL3502，LOD值為2.0，加性效應-0.96，顯性效應0，解釋表型變異8.59%；檢測到與馬克隆值相關的單標記位點1個，引物名稱為NAU1119，位于Chr26上，LOD值為3.0，加性效應-0.15，顯性效應3.7，解釋表型變異8.88%。表2

表2 BC1F2群體纖維品質性狀相關性Table 2 Analysis of correlation among fiber quality in BC1F2 population

3 討 論

海島棉和陸地棉是2個主要的異源四倍體栽培種，在產量、纖維品質和生態適應性等方面存在可互補性[18]。研究選用高純系主推的陸地棉品種新陸中60號為母本與優質海島棉品種新海41號為父本進行雜交構建作圖群體，將二者優良特性優勢結合，實現產量和纖維品質的同步提高。選取249對分布于棉花全基因組的SSR分子標記進行多態性篩選，共獲得105對差異明顯的多態性引物。當前開發的棉花分子標記引物多來自海陸群體，遺傳連鎖圖譜的加密主要依賴于開發的引物數量和親本遺傳基礎的差異大小。棉花海陸群體QTL定位研究趨于飽和。張正圣等[19]利用陸地棉種間雜交(渝棉1號×T586)F2群體，構建了一張包括173個標記位點(159個SSR,11個AFLP和3個形態標記)的遺傳連鎖圖，覆蓋1 360 cM，約占棉花基因組的3.06%。Ullao等[20]利用陸地棉種內雜交(MD5678ne-×Prema)F2群體構建的連鎖圖，包括17個連鎖群，81個RFLP 位點，標記間的平均距離8.7 cM，覆蓋700.7 cM，約占棉花基因組的15.8%。

因多效現象在自然界生物進化遺傳中普遍存在，在多環境多群體中獲得的 QTL位點往往有較高的穩定性和精確度。所獲得的1個位點在相關文獻中已有報道，qFL-14-1與qFM-14-1[21-22]、q FS-14-1[22-24]及q FU-14-1[2]在同一位點，效應不同，該位點對應的等位基因有可能連鎖遺傳，或者是一因多效[24]。在不同群體不同環境下定位到的QTL，能更準確地反映控制棉花纖維品質性狀的等位基因的位置，環境因素對數量性狀位點的檢測影響很大。

4 結 論

共獲得105對差異明顯的多態性引物。篩選的與優質性狀相關SSR標記可應用于棉花優質分子標記輔助選擇上。檢測到與纖維長度相關的QTL位點1個，位于Chr14的32.31 cM處，解釋表型變異8.59%。檢測到與馬克隆值相關的位點1個，位于Chr26上，解釋表型變異8.88%。