小麥冬前最大分蘗期根數的QTL定位

彭 濤，尹國紅，成東梅，于金林，趙偉峰，高 燕，田紀春，鄧志英

(1.山東農業大學/作物生物學國家重點實驗室，山東 泰安 271018； 2.濟源市農業科學院，河南 濟源 459002)

根系是植物吸收水分和養分的重要器官，其生長發育狀況直接影響植株地上部的發育和產量[1-4]。小麥是須根系作物，其根系由初生根和次生根組成[5-6]。初生根由種子胚根發育而成，一般于小麥種子萌發后7～10 d集中發生，功能期主要集中在出苗到拔節期。次生根的發生，一般從小麥植株分蘗期開始，到抽穗開花時停止發育[6]，其發生有2個高峰，一個是冬前最大分蘗期，另一個是拔節到抽穗期。薛麗華等[7]研究發現，冬前最大分蘗期的次生根數與小麥的有效穗數密切相關。李鴻斐等[8]研究表明，強筋小麥冬前最大分蘗期的次生根數顯著高于中筋和弱筋小麥。綜上，小麥冬前最大分蘗期的初生根數、次生根數及總根數等根系性狀對小麥的產量和品質具有重要的影響。

CARADUS[9]分析作物不同根系形態性狀的遺傳力發現，根數、根長、根體積、根表面積、根直徑、根質量、根冠比、根毛長度等衡量根系形態的指標具有較高的遺傳力，表現出典型的數量性狀遺傳。由于根系深埋在土壤中，存在工作量大、難以獲得完整準確的根系數據等問題，所以大多數根系研究主要利用水培或凝膠培養的方式在室內模擬進行根數、最大根長、根直徑、根表面積、根體積、根鮮質量、根干質量等的遺傳分析[10-15]，而在大田土壤環境條件下開展小麥根系性狀的QTL定位分析的研究較少[16]。迄今，僅見王升星等[16]在大田土壤環境條件下以196份小麥自然群體為材料，于小麥越冬前、拔節期進行次生根數測定，利用185對SSR標記進行基因型關聯分析，共關聯到32個顯著標記位點，分布于14條染色體上，其中，Barc81(1BL)、Wmc617(4DS)和Gwm190(5DS) 在3個環境下與次生根數穩定關聯。但其并未對總根數、初生根數進行QTL定位分析。為此，本研究利用在根系性狀方面差異顯著的花培3號和豫麥57及其構建的含有168個株系的加倍單倍體(DH)群體為材料，在3種土壤栽培環境條件下，對小麥冬前最大分蘗期根數(初生根數、次生根數和總根數)進行測定及QTL定位分析，旨在揭示小麥冬前最大分蘗期根數的分子遺傳基礎，以期找到有助于培育健壯根系的新的QTL，為小麥根系重要功能基因的克隆及分子標記輔助選擇育種奠定基礎。

1 材料和方法

1.1 試驗材料

試驗材料為168個DH群體，其由高產多抗的花培3號(母本)，優質、高產、綜合性狀優良的豫麥57(父本)雜交產生的F1通過花藥培養，并經染色體加倍獲得。

DH群體的遺傳圖譜由本實驗室于2008年構建[17]，含有定位于小麥21條染色體的323個位點(包括284個SSR、37個EST-SSR、1個ISSR和1個HMW-GS標記)，使用MapChart 2.1[18]軟件作圖。圖譜全長2 485.7 cM，平均2個標記間的遺傳距離是7.67 cM，形成24個連鎖群。

1.2 試驗設計

2013—2014年進行盆栽試驗，于山東農業大學網室內進行，簡稱泰安市盆栽試驗(E3)。選取DH群體168個家系及親本大小均勻一致、飽滿完整的種子各20粒，浸泡于10% H2O2中30 min，然后用蒸餾水沖洗干凈，播種于帶紗網孔的發芽板上，置入光照培養箱中。設置白天溫度為(22±2)℃，光照16 h，夜間溫度為(14±2)℃，黑暗8 h。種子萌發后，從每個家系中挑取12株生長健壯、長勢一致的幼苗，分別移栽到3個塑料盆中，重復3次。塑料盆底部有孔，盆深33 cm，上口直徑為32 cm，下口直徑為26 cm，盆栽土壤由75%農田土壤和25%細砂土混合而成，其中，農田土壤過篩去除石子、根莖等雜質。每盆裝土15 kg，裝土前拌入復合肥(含總養分45%，N∶P2O5∶K2O=15∶15∶15) 4 g。將土壤和肥料充分混勻后裝盆并澆水沉實，次日移栽。移栽日期為10月12日。在小麥幼苗生長期，及時防治雜草和病蟲害，進行正常的水肥管理。

2014—2015年進行田間栽培試驗，試驗材料種植在河南省濟源市農業科學院試驗田和山東農業大學試驗農場，簡稱濟源市田間試驗(E1)和泰安市田間試驗(E2)，完全隨機區組設計，2次重復，3行區種植，行長2 m，行距0.20 m，播種日期分別為10月5日和10月10日，基本苗150萬株/hm2。按當地小麥產量比較試驗方案進行田間管理。

1.3 測定項目及方法

對于盆栽試驗，于12月22日取樣，取樣前先用水管澆透各塑料盆土壤，再將塑料盆浸泡在水中，然后仔細清洗各單株根系。每個家系隨機選6株，統計小麥冬前最大分蘗期初生根數、次生根數和總根數。

對于田間試驗，山東省泰安市試驗田和河南省濟源市試驗田的小麥材料分別于12月20日和12月25日取樣，用鐵鏟小心挖出帶土幼苗，再把帶土幼苗浸泡在水中，然后仔細清洗各單株根系，統計小麥冬前最大分蘗期初生根數、次生根數和總根數。

1.4 數據處理和QTL分析

利用Excel 2003對數據進行匯總和整理，使用SPSS 19.0對數據進行正態分布檢驗、方差分析和相關性分析。

利用基于混合線性模型[19]的復合區間作圖軟件QTLNetwork 2.1[20]對小麥冬前最大分蘗期初生根數、次生根數和總根數等進行QTL分析。QTLNetwork 2.1軟件的掃描窗口設定為10 cM，步長設定為1 cM，以單個環境中3個重復的平均值輸入數據，以P=0.005為統計檢測閾值，即當某個根數性狀分子標記的P值小于統計檢測閾值時，認為該標記處存在1個與根數性狀有關的QTL。最后將檢測到的所有QTL以及它們之間的上位性互作整合到一個全QTL模型中，用基于Gibbs抽樣的貝葉斯方法估計遺傳效應。按照 MCINTOSH等[21]的方法命名QTL。

2 結果與分析

2.1 小麥冬前最大分蘗期根數分析

親本豫麥57的初生根數、次生根數和總根數均高于花培3號，尤其是次生根數和總根數。在DH群體中，冬前最大分蘗期初生根數、次生根數和總根數變異范圍較大。與2個田間試驗相比，盆栽條件下根數增多，變異范圍更大(表1)。不同環境條件下初生根數、次生根數和總根數3個性狀的偏度和峰度均較小，呈正態分布(表1，圖1—3)，說明該群體適合對根數進行QTL定位。

表1 小麥冬前最大分蘗期DH群體根數的分析Tab.1 Analysis of root number of wheat DH population at maximum tillering stage before winter

注：E1：濟源市田間試驗，E2：泰安市田間試驗，E3：泰安市盆栽試驗；PRN：初生根數，SRN：次生根數，RTN：總根數。下同。

Note:E1 represents field test in Jiyuan City,E2 represents field test in Tai’an City,and E3 represents pot culture test in Tai’an City;PRN represents primary root number,SRN represents secondary root number,and RTN represents total root number.The same as below.

圖1 3個環境下小麥DH群體初生根數分布Fig.1 Distribution of primary root number of wheat DH population under three environments

圖2 3個環境下DH群體次生根數分布

圖3 3個環境下DH群體總根數分布

2.2 小麥冬前最大分蘗期根數性狀間的相關性分析

由表2可知，小麥冬前最大分蘗期總根數與初生根數、次生根數均呈極顯著正相關。與初生根數相比，次生根數與總根數相關性更高，在3個環境下相關系數分別為0.900、0.964、0.732，說明總根數主要由次生根數決定。次生根數與初生根數均呈正相關，但均未達到顯著水平。

表2 小麥冬前最大分蘗期DH群體根數性狀間的相關性分析(E1/E2/E3)Tab.2 Correlation analysis among root number traits of wheat DH population at maximum tillering stage before winter

注：**表示在0.01水平相關性極顯著。

Note:** represents significant correlation at 0.01 level.

2.3 小麥冬前最大分蘗期根數的QTL定位及效應分析

3個栽培環境下，共定位了7個加性效應QTL (表3和圖4) 和7對上位性互作QTL (表4)，分布在1B、2B、2D、3B、4A、4D、5D、6B、6D、7D染色體上。

對于初生根數，僅在E2環境下檢測到1個控制初生根數的加性效應QTL，位于2B染色體上，表型貢獻率為7.03%，增效等位基因來源于母本花培3號。E2環境下未檢測到上位性互作QTL。

對于次生根數，在3個環境下均檢測到控制次生根數的加性效應QTL，共4個，分別位于1B、4A、4D染色體上，可解釋表型變異的4.67%～16.56%，其中2個為主效QTL，位于1B染色體上的qSrn1B具有最大的表型貢獻率，可解釋16.56%的表型變異，其增效基因來源于父本豫麥57；位于4D染色體上的qSrn4D可解釋10.70%的表型變異，其增效基因來源于母本花培3號。在E1和E3環境下，檢測到3對控制次生根數的上位性互作QTL，分別位于染色體1B-2B、2B-4A、3B-5D上，可分別解釋6.90%、8.24%、5.24%的表型變異。

對于總根數，在E2和E3環境下，共檢測到影響總根數的2個加性效應QTL，分別位于1B和6D染色體上，可分別解釋表型變異的12.80%和9.17%，其中，位于1B染色體上的qRtn1B是主效QTL，其增效基因來源于父本豫麥57。在E1和E3環境下，檢測到4對控制總根數的上位性互作QTL，分別位于染色體6B-7D、2B-4A、2D-6D、3B-5D上，可分別解釋14.45%、9.78%、7.57%、5.65%的表型變異，其中1對位于染色體6B-7D上的上位性互作 QTL(qRtn6B和qRtn7D)的表型貢獻率達14.45%，表明上位性效應對該群體冬前最大分蘗期總根數的影響較大。

表3 小麥冬前最大分蘗期DH群體根數的加性效應QTL Tab.3 Additive QTLs for root number of wheat DH population at maximum tillering stage before winter

■代表E2環境下的初生根數，○、◎ and ●分別代表E1、E2、E3環境下的次生根數，△、▲分別代表E2、E3環境下的總根數

總體來看，位于1B染色體上影響次生根數的QTL在2個泰安市環境下均被檢測到，該位點對次生根數能解釋4.67%～16.56%的表型變異，且該位點來源于父本豫麥57，為1個典型的主效QTL位點。但該位點在E1環境下未檢測到，需增加檢測年份或地點，以確定該位點在不同環境下的表達特性。在E3環境下，分別在2B-4A、3B-5D染色體上檢測到對次生根數、總根數有共同效應的上位性互作QTL，可解釋表型變異的5.24%～9.78%。

表4 小麥冬前最大分蘗期DH群體根數的上位性互作QTL位點Tab.4 Epistatic QTLs for root number of wheat DH population at maximum tillering stage before winter

3 結論與討論

因根系深埋在土壤中，目前國內外對小麥根系性狀的遺傳定位研究主要基于水培條件，而本研究在土壤條件下測定小麥DH群體的根數，其結果更具有實際參考價值。

本研究在田間環境和盆栽環境條件下，以花培3號和豫麥57雜交獲得的168個株系的DH群體為試驗材料，在土壤栽培條件下對冬前最大分蘗期的根系性狀進行了QTL定位分析，共定位了小麥冬前最大分蘗期初生根數、次生根數、總根數3個根系性狀的7個加性效應QTL和7對上位性互作QTL，分布在1B、2B、2D、3B、4A、4D、5D、6B、6D、7D染色體上，單個QTL可以解釋4.67%～16.56%的遺傳變異。根數的加性效應QTL分別定位在1B、2B、4A、4D、6D染色體上，除1B染色體上的位點來自于父本豫麥57外，其余位點均來自于母本花培3號。其中，單個QTL對根數的表型貢獻率大于10%的主效QTL有3個，包括2個次生根數QTLqSrn1B(16.56%)、qSrn4D(10.70%)和1個總根數QTLqRtn1B(12.80%)。在1B染色體的標記區間Xwmc406—Xbarc156同時存在控制次生根數、總根數的QTLqSrn1B、qRtn1B，分別解釋次生根數、總根數表型變異的16.56%、12.80%，可為小麥根系的分子標記輔助選擇奠定基礎。經對比分析發現，在1B染色體上Xwmc406—Xbarc156區間控制次生根數和總根數的QTL，與利用同一DH群體定位的小麥穗部可育小穗數[22]、小麥籽粒清蛋白和球蛋白含量[23]QTL位于相同區間位置，表明控制次生根數和總根數的基因可能同時控制小麥穗部可育小穗數、籽粒清蛋白和球蛋白含量，其對小麥產量和品質均能產生影響；在4A染色體上控制次生根數的QTL與利用同一DH群體定位的籽粒產量、總小穗數、可育小穗數[22]、幼苗根長[24]QTL位于相同區間位置，表明控制次生根數的基因可能同時控制小麥籽粒產量、總小穗數、可育小穗數和幼苗根長。

本研究在1B、4A、4D染色體上檢測到在大田、盆栽環境下影響越冬前最大分蘗期次生根數的4個加性效應QTL，而王升星等[16]在1A、1B、2D、3A、3B、4A、4D、5A、5B、5D、6B、6D、7A、7B等14條染色體上關聯到大田環境下控制越冬前、拔節期次生根數的32個顯著標記位點，其中越冬前次生根數顯著關聯標記有11個，分布在1B、2D、3A、4D、5A、5D染色體上，表明1B、4A、4D等染色體上存在控制次生根數的重要位點。本研究在2B染色體上檢測到1個控制初生根數的QTL，其所在標記區間與先前定位到的控制根數的QTL標記區間一致[13-15,25-28]。

通過比較發現，一些染色體區間是控制根系性狀的熱點區間，這與在不同作物中發現相關性狀QTL常常定位于相同或相近染色體區段的研究結論一致[22,29-31]。本研究中，在染色體1B的 Xwmc406—Xbarc156標記區間同時存在控制次生根數和總根數的3個QTL，其中，2個主效位點qSrn1B和qRtn1B分別解釋次生根數和總根數表型變異的16.56%和12.80%。王升星等[16]在該區間內關聯到了2個控制次生根數的QTL；屈春艷[32]在該區間關聯到了4個控制小麥根數的QTL，其中，QRn.sdau-1B.1的表型貢獻率達13.54%；李振興等[26]也在該區間定位了2個與根數相關的磷脅迫響應QTL，其中，1個QTL的表型貢獻率達16.77%。染色體1B上具有控制次生根數和總根數的重要位點，該位點可能存在一因多效基因或緊密連鎖基因，本研究為進一步在該區段內發掘控制根系相關性狀的重要基因以及分子標記育種奠定了基礎。