趙衛(wèi)國　塔娜　張立堅　李保軍　趙小萍　尚麗平　王灝　穆建新

摘 要 為了解析油菜開花期性狀的遺傳機制，利用KN DH群體在冬性、半冬性和春性環(huán)境的開花期表型和KN 高密度遺傳連鎖圖譜，通過Wincart 2.5軟件的符合區(qū)間作圖法對油菜開花期性狀進行QTL定位及候選基因鑒定。結果顯示，共鑒定到119個開花期QTL，單個QTL解釋表型變異最大是qFT-13DL16-4 ?（25.96%），最小的是qFT-13ZY2-1（2.48%）。利用元分析的方法將初步鑒定的QTL整合為consensus QTL，共獲得26個環(huán)境穩(wěn)定表達QTL，包括7個開花期主效QTL。如cqFT-A2-3、cqFT-A2-4在春性環(huán)境穩(wěn)定表達，cqFT-C6-4、cqFT-C6-7、cqFT-C6-12、cqFT-C6-13在冬性和半冬性環(huán)境穩(wěn)定表達， cqFT-C6-14在冬性環(huán)境穩(wěn)定表達QTL。主效QTL置信區(qū)間共鑒定到15個與成花誘導相關的候選基因，如BnaA02g12260D（RGA1）、BnaA02g15390D（AGL12）、BnaA02g16710D（LKP2）和BnaC06g19930D（NUA）等，這些候選基因主要涉及赤霉素、光周期、生物鐘、春化作用響應和花發(fā)育等功能。可見，油菜開花期主效QTL及其候選基因的鑒定為開花期基因的精細定位和圖位克隆奠定基礎，也為培育早熟、高產油菜品種提供指導。

關鍵詞 甘藍型油菜；雙單倍體株系；開花期性狀；QTL定位；候選基因鑒定

油菜是世界上重要的油料作物之一，也是中國食用植物油的主要來源，具有很高的營養(yǎng)和經濟價值［1］。目前，中國食用植物油自給率嚴重不足，60%以上需要進口［2］。高產穩(wěn)產的油菜品種是油菜育種的主要目標。油菜開花期由多基因控制且受外界環(huán)境影響較大。在不同的氣候、生態(tài)環(huán)境，油菜開花期表現出較大的表型差異［3］。根據不同的生態(tài)環(huán)境，可將油菜分為冬性、半冬性和春性3種生態(tài)類型［4］。冬性油菜需要一定的低溫才能通過春化，春化持續(xù)時間較長；半冬性油菜對春化的要求相對較低，春化持續(xù)時間也相對較短；春性油菜開花幾乎不經過低溫春化就能夠正常開花［5］。此外，油菜作為重要的倒茬作物，開花期過長將會導致油菜成熟期的推遲，從而影響了下茬作物的正常栽培［6-7］。因此，必須嚴格控制油菜的開花期從而調節(jié)油菜的成熟期。在油菜品種選育過程中，開花期可作為油菜早熟品種選育的重要指標。

數量性狀定位（QTL）是闡明數量性狀遺傳基礎和識別QTL位點的有效途徑［8］。目前，科研工作者已經成功對玉米花期［9］、油菜株高［10］、水稻粒形、粒質量和粒長［11-12］，小麥萌發(fā)期抗旱和耐鹽［13］等重要性狀進行QTL定位研究，并且也獲得相關性狀的QTL。油菜開花期屬于數量性狀，不同的定位群體鑒定的開花期QTL數量和在染色體上的分布具有較大差異。Ferreira等［14］對以甘藍型冬、春油菜構建的DH群體開花期進行QTL定位，鑒定到3個開花期QTL分別位于A09、C02和C06 染色體上；Butruille等［15］以冬、春油菜的雜交和回交群體作為定位群體，鑒定到的開花期QTL主要位于A02、A03、A07、A08、A09、C02 和 C05 染色體上；Zhao等［16］對油菜DH群體的株高、開花期和成熟期進行QTL鑒定，檢測到7個開花期顯著性QTL；蔡長春等［17］對油菜開花期進行QTL定位，鑒定到分別位于A04和A06染色體上的2個主效QTL。這些研究結果充分表明油菜開花期性狀由不同染色上的多個位點控制，表現出開花期性狀調控的復雜性。同時，許多作物開花期QTL的定位也為本研究甘藍型油菜開花期的QTL定位提供理論基礎和技術支持。

QTL定位是獲得區(qū)間候選基因的先決條件。隨著甘藍型油菜基因組的公布，利用甘藍型油菜基因組數據庫（http：//www.genoscope.fr/ brassicanapus/），通過QTL置信區(qū)間的物理距離與油菜基因組的線性關系，可獲得QTL區(qū)間潛在基因。由于甘藍型油菜與擬南芥同源及基因的保守性，利用擬南芥基因組數據庫與甘藍型油菜候選基因進行Blast比對，確定甘藍型油菜基因在擬南芥中的同源基因及功能注釋［18］。目前，已經鑒定到油菜開花期基因BnCOA1與擬南芥CO基因同源［19］，BnFLC1與擬南芥FLC基因同源等［20］。因此，通過QTL區(qū)間及其區(qū)間的潛在基因進行鑒定和功能解析，以期獲得油菜開花期的潛在候選基因。

目前，對油菜定位群體同時在冬性、半冬性和春性3個生態(tài)類型的開花期表型進行QTL定位，相關研究還相對較少。本研究主要通過利用油菜KN Double haploid（DH）群體在3個生態(tài)環(huán)境（冬性、半冬性和春性生態(tài)環(huán)境）連續(xù)多年的開花期表型，并結合高密度遺傳連鎖圖譜對甘藍型油菜開花期性狀的QTL進行定位、分析，在分子水平解析油菜開花期的調控機理，并獲得控制油菜開花期重要QTL及潛在候選基因。研究結果可為油菜開花期QTL的精細定位和基因的圖位克隆奠定堅實基礎，同時也為培育早熟油菜新品種提供新的基因資源和指導依據。

1 材料與方法

1.1 試驗材料

2007年以春性甘藍型油菜純系KenC-8為父本，冬性甘藍型油菜純系N53-2為母本，通過人工雜交和小孢子培養(yǎng)構建了由348份Double haploid（DH）系組成的KN DH群體。該研究以甘藍型油菜KN DH群體為試驗材料，對其不同生態(tài)環(huán)境開花期表型進行考察。

1.2 田間試驗

以 KN DH 群體及其親本作為試驗材料，2010-2013年和2015-2018年8個年份種植在冬性環(huán)境陜西省雜交油菜研究中心（大荔）試驗田（編碼為 DL，經緯度為 109.94°E，34.79°N）和陜西省雜交油菜研究中心（楊陵）試驗田（編碼為YL，經緯度為 108.07°E，34.28°N）；2012-2015年連續(xù)4 a種植在半冬性環(huán)境湖北武漢華中農業(yè)大學油菜試驗田（編碼為 WH，經緯度為 113.68°E，30.58°N）， 2010-2012 年連續(xù) 3 a種植在春性環(huán)境甘肅省張掖油菜加代試驗田（編碼為 ZY，經緯度為100.45°E，38.93°N），并記錄每年試驗的播種時間。以種植年份和地點編碼組合對各個種植環(huán)境進行命名，如2013年種植在湖北武漢命名為13WH。KN DH群體在冬性環(huán)境試驗每年3個重復，半冬性和春性環(huán)境試驗每年皆2個重復。田間試驗采用完全隨機區(qū)組設計排列，每個DH系播種2行，行距40? cm，株距20? cm，按照油菜正常的試驗管理方式進行管理。

1.3 數據統(tǒng)計與分析

油菜初花期是以該材料中50%以上株系主花序第一朵花開為標準，記錄為該材料的初花期。將初花期和播種期之間的時間進行計算，記錄為該品系的開花期。在甘肅省張掖春性環(huán)境，由于KN DH群體是由春性和冬性親本構建而來，KN DH群體母本N53-2和部分DH系對春化溫度和持續(xù)時長要求差異導致不能正常開花或者初花期推遲過長。將在油菜生長期不能正常開花或者其他DH系收獲時還沒有達到初花期的標準，以同年其他品系成熟期作為初花期來計算這部分DH系的開花期［20］。利用Excel 2010軟件對開花期表型數據進行處理。

1.4 開花期的QTL定位與QTL整合

通過KN DH群體的高密度SNP遺傳連鎖圖譜，結合各生態(tài)環(huán)境油菜的開花期表型，利用Wincart 2.5軟件進行QTL檢測分析［21］。軟件掃描以種植環(huán)境作為環(huán)境因子處理，采用復合區(qū)間作圖法（CIM），通過逐步回歸方法選擇協(xié)因子，背景標記數目設為5，窗口大小設為10 cM，掃描步長設為1.0 cM進行QTL檢測。以P=0.05分別對每個環(huán)境的表型進行1 000次排布測驗方法（permutation test）確定LOD閾值，最終將顯著性水平的LOD閾值估算為2.5。 當LOD值≥LOD閾值時的QTL稱為顯著性identified QTL［22］。Identified QTL命名以“q性狀縮寫+種植年份地點編碼+染色體編號” ［23］，若該性狀在同一種植環(huán)境同一染色體上檢測到多個顯著性QTL，將依次以“-1，-2”等來區(qū)分，如qFT.12ZY2-1代表2012年在甘肅張掖試驗點檢測到位于A02染色體上的第1個QTL。

利用BioMercator 4.2軟件對不同環(huán)境中 ?identified QTL進行整合獲得consensus QTL［24］。將consensus QTL以“cq”+“性狀縮 ?寫+染色體編號”+“在該染色體上第幾個QTL”命名，如“cqFT.A2-1”表示開花期consensus QTL位于A02染色體第1個QTL。將在兩個種植環(huán)境表型變異≥10%或者一個種植環(huán)境表型 ?變異≥20%的consensus QTL 認定為主效QTL［25］。

1.5 開花期QTL置信區(qū)間候選基因鑒定與分析

結合甘藍型油菜“Darmor-bzh”參考基因組，將與QTL連鎖的50-mer SNP探針序列在NCBI數據庫進行Blast，同源位點作為QTL區(qū)域候選基因。KN高密度SNP遺傳連鎖圖譜上所有的標記用于確定遺傳連鎖圖和物理圖的線性關系［26］。如果1個SNP比對到參考基因組多個位點區(qū)域，僅選擇與該位點相一致的連鎖群位點，如果1個SNP在同一連鎖群上有多個同源位點，通過該位點上下游物理位置手動進行確認［27］。將與QTL置信區(qū)間對應一致的基因組區(qū)域內的所有基因認定為該QTL區(qū)域內的候選基因［28］。結合擬南芥數據庫對QTL區(qū)間的候選基因進行同源基因鑒定和功能注釋分析。

2 結果與分析

2.1 油菜KN DH群體在3個生態(tài)環(huán)境中開花期的表型鑒定

對KN DH群體、父本KenC-8和母本N53-2在冬性、半冬性和春性生態(tài)環(huán)境多年的開花期表型進行統(tǒng)計分析（表1）。KN DH群體的父本KenC-8在3個生態(tài)環(huán)境的開花期早于母本N53-2，父本KenC-8的開花期皆短于母本N53-2。在冬性環(huán)境的陜西大荔和楊陵、半冬性環(huán)境的湖北武漢，父本KenC-8和母本N53-2的開花期差異較小，父本KenC-8比母本N53-2提前0.2～ ?3.0 d開花。在春性環(huán)境兩親本開花期差異較大，父本KenC-8比母本N53-2提前6.5～49.0 d開花。這是由于母本N53-2的冬性較強，需要較長時間的低溫才能完全通過春化，從而導致母本N53-2初花期推遲，開花期延長。

同時，對KN DH群體在冬性、半冬性和春性3個生態(tài)環(huán)境的開花期表型分析，該群體在冬性生態(tài)環(huán)境的開花期時間最長（189.7～197.6? d），其次是半冬性環(huán)境（149.5～176.7 d），開花期最短的是在春性環(huán)境（92.2～94.7 d）。油菜在冬環(huán)境和半冬性環(huán)境，第一年的9-10月播種，第二年的3-4月開花，導致這兩個生態(tài)環(huán)境開花期時間較長，而在春性環(huán)境當年4月份種植，7月初就開始開花，因此在春性環(huán)境開花期較短。KN DH群體在春性環(huán)境的開花期變幅最大，變幅為 ?73.0 d（63.0～136.0 d），其次在冬性和半冬性的開花期變幅相差不大，分別為44.3 d（161.0～205.3 d）和45.5 d（138.5～184.0 d）。油菜KN群體開花期表型與兩親本開花期相比較，KN群體的開花期表現出超親的遺傳現象。此外，對KN DH群體3個生態(tài)環(huán)境開花期表型分布分析，發(fā)現油菜開花期表型基本服從正態(tài)或者近似的正態(tài)分布，表現出開花期屬于數量性狀的遺傳特點 ?（圖1）。

2.2 油菜開花期性狀的QTL定位與分析

利用WinCart 2.5軟件對油菜KN群體3個生態(tài)環(huán)境的15個種植環(huán)境開花期表型數據，結合KN高密度SNP遺傳連鎖圖譜對油菜開花期性狀性狀進行QTL定位于分析。首選是將開花期表型和高密度遺傳連鎖圖譜通過WinCart 2.5軟件進行1 000次排布確定顯著性QTL的LOD閾值為2.5。通過WinCart 2.5軟件的初步掃描，共鑒定到119個identified QTL，主要分布在A03、A06、C06等11條染色體上（表2，圖2），其中69個identified QTL分布在A基因組上，50個分布在C基因組。C06染色體上的開花期QTL數量最多（32個），其次是A03染色體（28個），在A09和C05染色體上僅檢測到1個QTL。進一步分析發(fā)現開花期QTL的表型變異最高的是位于C06染色體上的qFT-13DL16-4，其表型變異達到25.96%；其次是位于A02染色體上的qFT-13ZY2-3，其表型變異為23.84%；表型變異最小的則是位于A02染色體上的qFT-13ZY2-1，僅為2.48%。

2.3 油菜開花期性狀的QTL整合與分析

在不同環(huán)境檢測到部分identified QTL置信區(qū)間具有重疊區(qū)域，通過利用BioMercator4.2軟件對identified QTL進行整合。119個開花期 ?identified QTL被整合為76個consensus QTL（表2，圖3），大大縮小QTL置信區(qū)間，提高了QTL的準確性。在consensus QTL中，26個consensus QTL能夠在多個種植環(huán)境中檢測到，屬于環(huán)境穩(wěn)定表達QTL；50個consensus僅在1個種植環(huán)境表達，屬于環(huán)境特異表達QTL。在26個環(huán)境穩(wěn)定表達QTL中，14 個consensus QTL能夠在多個生態(tài)環(huán)境穩(wěn)定表達（表2），如 cqFT-A3-1和cqFT-A6-6在冬性、半冬性和春性環(huán)境都能夠穩(wěn)定表達，cqFT-A3-2和cqFT-A10-2在半冬性和春性環(huán)境穩(wěn)定表達，cqFT-C3-2和cqFT-C6-3在冬性和半冬性環(huán)境穩(wěn)定表達；12 個consensus QTL僅在一個生態(tài)環(huán)境穩(wěn)定表達，如 cqFT-A2-3和cqFT-A2-4僅在春性環(huán)境穩(wěn)定表達，cqFT-A6-7和cqFT-A6-8只在半冬性環(huán)境穩(wěn)定表達，僅有cqFT-A6-5在冬性環(huán)境穩(wěn)定表達。

更為重要的是在A02和C06染色體上鑒定到7個開花期主效QTL，包括cqFT-A2-3、cqFT-A2-4、cqFT-C6-4、cqFT-C6-7、cqFT-C6-12、cqFT-C6-13和 cqFT-C6-14，其在2個種植環(huán)境的QTL表型變異達到10%以上或在一個種植環(huán)境達到20%以上（表2）。位于A02染色體上主效QTL cqFT-A2-3和cqFT-A2-4（圖4），其加性效應都為正值，表明該QTL的貢獻率主要來自母本N53-2，并發(fā)現這2個主效QTL僅在春性環(huán)境穩(wěn)定表達。位于C06染色體上主效QTL（圖4），其加性效應都為負值，表明該QTL的貢獻率主要來自父本KenC-8。主效QTL cqFT-C6-4、cqFT-C6-7、cqFT-C6-12和cqFT-C6-13 在冬性和半冬性環(huán)境都能夠穩(wěn)定表達，cqFT-C6-14僅在冬性環(huán)境穩(wěn)定表達。主效QTL的獲得及其QTL的加性效應貢獻親本的確定對于下一步選擇正確親本進行開花期近等基因系的構建和該位點重要基因的圖位克隆，以及利用該位點進行適宜花期材料的創(chuàng)制具有重要的指導作用。

2.4 油菜開花期主效QTL置信區(qū)間候選基因鑒定

根據甘藍型油菜KN DH群體的高密度遺傳連鎖圖與甘藍型油菜“Darmor-bzh”基因組間的共線性關系，從7個油菜開花期主效QTL置信區(qū)間中鑒定出15個與成花誘導相關的候選基因（表3）。這些基因涉及油菜不同的開花路徑，包括赤霉素、光周期、生物鐘、春化作用響應和花發(fā)育等。主效QTL cqFT.A2-3和cqFT.A2-4在春性環(huán)境穩(wěn)定表達，其置信區(qū)間均鑒定到 BnaA02g12130D（FT） 和BnaA02g12260D ?（RGA1），這2個基因均參與調控開花整合子和赤霉素信號轉導途徑。在主效QTL cqFT.A2-4的置信區(qū)間鑒定到BnaA02g15390D（AGL12）、BnaA02g16710D（LKP2）、Bna02g20160D（SNZ），這3個基因主要參與了油菜的光周期和生物鐘路徑中的開花誘導過程。位于C06 染色體的主效QTL cqFT.C6-4、 cqFT.C6-7、cqFT.C6-12和cqFT.C6-13在冬性和半冬性環(huán)境能夠穩(wěn)定表達， cqFT.C6-14能夠在冬性環(huán)境穩(wěn)定表達（圖4）。除了cqFT.C6-13，其他主效QTL區(qū)間均鑒定到BnaC06g24850D（MMP）和BnaC06g27090D （TSF），BnaC06g24850D的擬南芥同源基因主要調控花的發(fā)育，BnaC06g27090D的擬南芥同源基因主要受到光周期的調控來調節(jié)擬南芥的開花誘導；在cqFT.C6-4、cqFT.C6-13和cqFT.C6-14的置信區(qū)間鑒定到BnaC06g32640D （AGL12），其擬南芥同源基因功能涉及光周期和生物鐘誘導途徑；cqFT.C6-4、cqFT.C6-7和 cqFT.C6-14的置信區(qū)間還鑒定到BnaC06g19930D（NUA），其擬南芥同源基因功能主要是對春化作用的響應。總之，油菜開花期在春性環(huán)境主要受到赤霉素的誘導途徑，在冬性環(huán)境和半冬性環(huán)境主要受到光周期、生物鐘和春化作用響應的誘導開花路徑。此外，在這些QTL置信區(qū)間還鑒定到如BnaC06g25500D（AP1）、BnaC06g19420D（SAR3）和BnaC06g27170D ?（ARF8） 等參與花發(fā)育調控路徑。

3 討? 論

開花期作為作物的重要農藝性狀之一，研究油菜開花期對培育優(yōu)質早熟油菜品種具有重大意義。許多研究者通過構建DH系、F2：3家系等不同類型群體對油菜開花期進行QTL定位，鑒定到許多主要QTL，為進一步解析開花期的遺傳機理奠定基礎［29］。本研究利用一個包含348個株系的KN DH群體，通過將該群體在冬性、半冬性和春性生態(tài)環(huán)境的多年開花期表型數據進行分析，發(fā)現KN DH群體及親本在不同生態(tài)條件下，油菜開花期的表型差異較大。在冬性和半冬性環(huán)境兩親本開花期差異較小，在春性環(huán)境父本墾C-8早于母本N53-2提前開花，且開花期差異較大。KN DH群體在冬性和半冬性環(huán)境下由于經過寒冷的冬季延長了植株的營養(yǎng)生長時間導致開花期較長，不同DH系之間開花期變幅差異較小，KN DH群體開花期在春性環(huán)境相對較短，變幅相差較大（73.0 d）。油菜開花期表現出超親的遺傳特點，而且對不同DH系的開花期頻數分布分析表現出正態(tài)或近似正態(tài)分布的特點。這些結果表明開花期屬于數量性狀遺傳，且受到種植環(huán)境的影響，該結論與前人研究結果一致［30-31］。

QTL定位是研究數量性狀最有效的方法之一［32］。關于甘藍型油菜重要農藝性狀的QTL定位和分析的研究已有大量研究，其中對開花期QTL定位是油菜作物研究的熱點之一［33］。Mei等［29］利用F2：3群體和遺傳連鎖圖獲得了6個開花期QTL，其中位于C03染色體上的2個主效QTL的表型變異超過30%。Wei等［34］對油菜開花期的QTL進行定位，共檢測到17個QTL，其中位于A10和C01染色體上的3個主效QTL解釋的表型變異均超過20%。魏雅琪等［31］以油菜的F2：3群體為材料進行開花期QTL定位，獲得了5個位于A02、A06、C02、C03、C06的半冬性環(huán)境穩(wěn)定表達QTL。Luo等［35］對以春油菜為材料的DH群體進行開花期QTL定位，鑒定到12個開花期QTL。這些研究的定位群體只是包含了兩個生態(tài)環(huán)境，而且由于定位群體和遺傳連鎖圖譜的差異，開花期QTL數量和分布差異較大。在本研究中，油菜KN DH群體的種植地點的生態(tài)類型包含了冬性、半冬性和春性環(huán)境，油菜開花期表型數據全面，得到的結果準確可靠。通過對KN DH群體在冬性、半冬性和春性環(huán)境近10年的開花期性狀進行QTL定位解析，共鑒定到119個QTL，解釋表型變異在2.48%～25.96%。本研究獲得的位于A02、C02、C03和C06染色體上的QTL與前人研究鑒定的QTL所在染色體結果一致［33，36］。然而，本研究與前人研究獲得的QTL是否屬于相同位點，下一步工作還需要借助油菜參考基因組進行比對分析才能獲得結果。

環(huán)境穩(wěn)定QTL，特別是主效QTL的鑒定對于下一步通過QTL位點聚合和分子標記輔助選擇培育廣適性的早熟油菜新品種（系）具有重要指導作用，同時也有利于為下一步開花期主效QTL位點的近等基因系構建和圖位克隆奠定基礎。本研究鑒定到76個consensus QTL，其中26個 consensus QTL屬于環(huán)境穩(wěn)定表達QTL，包括cqFT-A2-3和cqFT-A2-4僅在春性環(huán)境穩(wěn)定表達，cqFT-A3-1和cqFT-A6-6、cqFT-C6-4在3個生態(tài)環(huán)境都能夠穩(wěn)定表達，cqFT-A3-2和cqFT-A10-2在半冬性和春性環(huán)境穩(wěn)定表達， cqFT-A3-9、 cqFT-C6-2、cqFT-C6-3、 cqFT-C6-7和cqFT-C6-13等在冬性和半冬性環(huán)境穩(wěn)定表達。沒有檢測到在冬性、春性環(huán)境穩(wěn)定表達QTL。此外，對于不同穩(wěn)定表達QTL，其加性效應來源不同，如cqFT-A2-3、cqFT-A3-1、cqFT-A6-6等QTL的加性效應為正值，其效應主要來自母本N53-2，位于C06染色體上的環(huán)境穩(wěn)定表達QTL 如cqFT-C6-3、 cqFT-C6-4、 cqFT-C6-7、 cqFT-C6-10、 cqFT-C6-12和cqFT-C6-13等QTL的加性效應皆為負值，其效應主要來自父本KenC-8。因此，在利用某個QTL位點進行近等基因系構建時，要特別注意其效應來自哪個親本，以免搞錯供體親本和受體親本耽誤科研進程。

油菜開花期的調控機制較為復雜，內部基因調控和外部環(huán)境信號刺激共同作用才能促使油菜正常開花。在以往的油菜花期性狀QTL定位研究中，許多研究者只是對油菜開花期進行了顯著性QTL檢測［32-33，35］，而對于其QTL置信區(qū)間候選基因未進行進一步鑒定。利用甘藍型油菜基因組數據庫，通過與KN高密度SNP遺傳連鎖圖的線性比對，獲得了位于7個主效QTL區(qū)間的15個候選基因及擬南芥同源基因。油菜和擬南芥具有共同的祖先，且擬南芥作為蕓薹屬十字花科的模式植物，由于擬南芥基因功能的保守性，在甘藍型油菜中的基因可能具有相似的功能。這15個候選基因主要涉及調控開花過程中赤霉素、光周期、生物鐘、花發(fā)育、春化作用等開花調控誘導途徑。半冬性環(huán)境穩(wěn)定表達QTL? cqFT.C6-7對應的候選基因為NUA（BnaC06g19930D）、MMP（BnaC06g24850D）和AGL12（BnaC06g32640D）等，參與調控花發(fā)育、春化和光周期等途徑。AGL12是SOC，FLC和LFY的上游正向調節(jié)子，可有效的促進植物提早開花［36］。這些研究結果充分表明油菜開花期的遺傳機理較為復雜，受到多條路徑的共同調控。

4 結? 論

中國是世界上第二大油菜生產國。作為水稻、玉米等重要倒茬作物，解析油菜開花期性狀對于培育早熟油菜品種具有重要意義。本研究分析近10個年份的3個生態(tài)環(huán)境的油菜KN DH群體開花期表型，油菜開花期表現出數量性狀的遺傳特點，受到生態(tài)環(huán)境條件的影響。共檢測到油菜開花期119個identified QTLs，表型變異最大的是qFT-13DL16-4（25.96%）；QTL整合后獲得26個環(huán)境穩(wěn)定表達consensus QTLs，其中包括在春性環(huán)境穩(wěn)定表達的主效QTL cqFT-A2-3和cqFT-A2-4，在冬性和半冬性環(huán)境穩(wěn)定表達的主效QTL cqFT-C6-4、cqFT-C6-7、cqFT-C6-12和cqFT-C6-13，在冬性環(huán)境穩(wěn)定表達主效QTL cqFT-C6-14。7個主效QTL區(qū)間的15個候選基因［BnaA02g12260D（RGA1），BnaA02g15390D（AGL12）、BnaA02g16710D（LKP2）和BnaC06g19930D（NUA）等］，其同源基因功能主要涉及了赤霉素、光周期、生物鐘、春化作用響應和花發(fā)育等。 本研究下一步將對主效QTL區(qū)域的候選基因進行功能驗證，以期挖掘出控制開花期的關鍵基因，為利用該基因位點加快培育優(yōu)質早熟油菜品種，縮短育種周期奠定堅實基礎和理論支持。

參考文獻 Reference：

［1]CAI G，YANG Q，YANG Q，et al.Identification of candidate genes of QTLs for seed weight in Brassica napus through comparative mapping among Arabidopsis and Brassica species［J］.BMC Genetics，2012，13：105.

［2]王漢中，殷 艷.我國油料產業(yè)形勢分析與發(fā)展對策建議［J］.中國油料作物學報，2014，36（3）：414-421.

WANG H ZH，YIN Y.Analysis and strategy for oil crop industry in China［J］.Chinese Journal of Oil Crop Sciences，2014，36（3）：414-421.

［3]PUTTERILL J，LAURIE R，MACKNIGHT R.Its time to flower：the genetic control of flowering time ［J］.Bioessays，2004，26（4）：363-73.

［4]官春云，王國槐，李 栒，等.油菜光溫生態(tài)特性的研究和應用［J］.作物研究，1995（1）：12-14.

GUAN CH Y，WANG G H，LI L，et al.Study and application of light temperature ecological characteristics of rape ［J］.Crop Research，1995（1）：12-14.

［5]LIU L，ZHANG J，ADRIAN J，et al.Elevated levels of MYB30 in the phloem accelerate flowering in Arabidopsis through the regulation of? flowering locust［J］.PLoS One，2014，9（2）：e89799.

［6]LI L，LONG Y，ZHANG L，et al.Genome wide analysis of flowering time trait in multiple environments via high-throughput genotyping technique in Brassica napus L ［J］.PLoS One，2015，10（3）：e0119425.

［7]高永同.甘藍型早熟油菜的花期選擇［J］.湖北農業(yè)科學，1979（11）：20.

GAO? Y T.Flowering period selection of early maturing Brassica napus ［J］.Hubei Agricultural Sciences，1979（11）：20.

［8]WANG X，CHEN L，WANG A，et al.Quantitative trait loci analysis and genome-wide comparison for silique related traits in Brassica napus ［J］.BMC Plant Biology，2016，16：71.

［9]蘭進好.玉米開花期相關性狀的QTL分析［J］.西北植物學報，2010，20（3）：471-480.

LAN J H.QTL analysis of flowering related traits in maize ［J］.Acta Botanica Boreali-Occidentalia Sinica，2010， ?20（3）：471-480.

［10][ZK（#]魏麗娟，劉瑞影，張 莉，等.甘藍型油菜莖高QTL定位及株高相關位點整合［J］.作物學報，2019，45（6）：818-828.

WEI L J，LIU R Y，ZHANG L，et al.QTL mapping of stem height and integration of plant height related loci in Brassica napus ［J］.Acta Crops Sinica，2019，45 （6）：818-828.

［11]林荔輝，吳為人.水稻粒型和粒重的QTL定位分析［J］.分子植物育種，2003，1（3）：337-342.

LIN L H，WU W R.QTL mapping analysis of grain type and grain weight in rice ［J］.Molecular Plant Breeding，2003，1（3）：337-342.

［12]鄭躍濱，李 智，趙海燕，等.水稻粒長QTL定位與主效基因的遺傳分析［J］.西北植物學報，2020，40（4）：58-64.

ZHENG Y B，LI ZH，ZHAO H Y，et al.Mapping quantitative trait loci associated with grain length and genetic? ?analysis of major quantitative loci in rice ［J］.Acta Botanica Boreali-Occidentalia Sinica，2020，40（4）：58-64.

［13] 王秀華，于鴻翔，孫福來，等.小麥萌發(fā)期抗旱和耐鹽性狀的QTL分析［J］.西北植物學報，2020，40（2）：62-73.

WANG X H，YU H X，SUN F L，et al.QTL mapping of drought and soult tolerant traits in wheat at germination stage ［J］.Acta Botanica Boreali-Occidentalia Sinica，2020，40（2）：62-73.

［14]FERREIRA M E，SATAGOPAN J，YANDELL B S，et al.Mapping loci controlling vernalization requirement and flowering time in Brassica napus ［J］.Theoretical and Applied Genetics，1995，90（5）：727-732.

［15]BUTRUILLE D V，GURIES R P，OSBORN T C.Linkage analysis of molecular markers and quantitative trait loci in populations of inbred backcross lines of Brassica napus L ［J］.Genetics，1999，153（2）：949-964.

［16]ZHAO J Y，BECKER H C，DING H D，et al.QTL of three agronomically important traits and their interactions with environment in a European x Chinese rapeseed population［J］.Journal of Genetics and Genomics，2005， ?32（9）：969-978.

［17]蔡長春，傅廷棟，陳寶元，等.甘藍型油菜遺傳圖譜的構建及開花期的QTL分析［J］.中國油料作物學報，2007， ?29（1）：1-8.

CAI CH CH，FU T D，CHEN B Y，et al.Construction of genetic map and QTL analysis of flowering period in Brassica napus ［J］.Acta Oilseed Crops Sinica，2007，29 （1）：1-8.

［18]CHALHOUB B，DENOEUD F，LIU? S，et al.Plant genetics.Early allopolyploid evolution in the post-Neolithic Brassica napus oilseed genome ［J］.Science，2014，345（6199）：950-953.

［19]ROBERT L S ，ROBSON F ，SHARPE A ，et al.Conserved structure and function of the Arabidopsis flowering time gene CONSTANS in Brassica napus ［J］.Plant Molecular Biology，1998，37（5）：763-772.

［20]TADEGE M ，SHELDON C C ，HELLIWELL C A ，et al.Control of flowering time by FLC orthologues in Brassica napus ［J］.Plant Journal，2001，28（5）：545-553.

［21]WANG S，BASTEN C，ZENG Z.Windows QTL Cartographer 2.5［M］.Raleigh，NC：North Carolina State University，2007.

［22]BURNS M J，BARNES S R，BOWMAN J G，et al.QTL? ?analysis of an intervarietal set of substitution lines in Brassica napus：（i） Seed oil content and fatty acid composition ［J］.Heredity （Edinb），2003，90（1）：39-48.

［23]WANG X，WANG H，LONG Y，et al.Identification of QTLs associated with oil content in a high-oil Brassica napus cultivar and construction of a high-density consensus map for QTLs comparison in B.napus ［J］.PLoS One，2013，8（12）：e80569.

［24]ARCADE A，LABOURDETTE A，FALQUE M，et al.Bio-Mercator：integrating genetic maps and QTL towards discovery of candidate genes ［J］.Bioinformatics，2004， ?20（14）：2324-2326.

［25]SHI J Q，LI R Y，QIU D，et al.Unraveling the complex trait of crop yield with quantitative trait loci mapping in Brassica napus ［J］.Genetics，2009，182：851-861.

［26]LU K，XIAO Z，JIAN H，et al.A combination of genome-wide association and transcriptome analysis reveals candidate genes controlling harvest index-related traits in Brassica napus ［J］.Scientific Reports，2016，6：36452.

［27]CHAO H，WANG H，WANG X，et al.Genetic dissection of seed oil and protein content and identification of networks associated with oil content in Brassica napus ［J］.Scientific Reports，2017，7：46295.

［28]SHI J，ZHAN J，YANG Y，et al.Linkage and regional association analysis reveal two new tightly-linked major-QTLs for pod number and seed number per pod in rapeseed （Brassica napus L.） ［J］.Scientific Report，2015，5：14481.

［29]MEI D S，WANG H Z，HU Q，et al.QTL analysis on plant height and flowering time in Brassica napus ［J］.Plant Breeding，2010，128（5）：458-465.

［30]JIAN H，ZHANG A，MA J，et al.Joint QTL mapping and transcriptome sequencing analysis reveal candidate flowering time genes in Brassica napus L ［J］.BMC Genomics，2019，20（1）：21.

［31]魏雅琪，危文亮，劉道敏，等.油菜開花期QTL定位和候選基因分析［J］.中國油料作物學報，2019，41（5）：679-687.

WEI Y Q，WEI W L，LIU D M，et al.QTL mapping and candidate gene analysis of flowering stage in Brassica napus ［J］.Acta Oil Seed Crops Sinica，2019，41（5）：679-687.

［32]LONG Y，SHI J，QIU D，et al.Flowering time QTL analysis of oilseed Brassica in multiple environments and genome-wide alignment with Arabidopsis ［J］.Genetics，2008，177（4）：2433-2444.

［33]ZHAO J Y ，BECKER H C ，DING H D ，et al.QTL of three agronomically important traits and their interactions with environment in a European x Chinese rapeseed population ［J］.Journal of Genetics and Genomics，2005， ?32（9）：969-978.

［34]WEI D，MEI J，FU Y，et al.Quantitative trait loci analyses for resistance to sclerotinia sclerotiorum and flowering time in Brassica napus ［J］.Molecular Breeding，2014， ?34（4）：1797-1804.

［35]LUO Y，LUO D，DU D Z，et al.Quantitative trait analysis of flowering time in spring rapeseed （B.napus L.） ［J］.Euphytica，2014，200：321-335.

［36]TAPIA-LOPEZ R，GARCIA-PONCE B，DUBROVSKY J G，et al.An AGAMOUS-related MADS-box gene，XAL1 （AGL12），regulates root meristem cell proliferation and flowering transition in Arabidopsis ［J］.Plant Physiology，2008，146（3）：1182-1192.

QTL Mapping and Candidate Gene Identification of Brassica napus at Flowering Stage

ZHAO Weiguo1，TA Na1， ZHANG Lijian2， LI Baojun2， ZHAO Xiaoping2， SHANG Liping2， WANG Hao2 and? MU Jianxin2

Abstract This study aimed to investigate the genetic mechanism of Brassica napus at flowering period using? the double haploid （DH） population and KN population KN high-density genetic linkage map.The phenotype data at flowering period under winter， semi-winter and spring environments were analyzed by using composite interval mapping （CIM） of Wincart 2.5 software to locate the QTL for flowering-related traits and identify candidate genes.The results revealed a total of 119 associated with flowering stage，with qFT-13DL16-4 explaining up to 19.63% of the phenotypic variation （PV） and qFT-13ZY2-1 exhibiting the minimum PV （2.48%）.Furthermore， these QTLs were integrated into a consensus QTL by meta-analysis， resulting in the identification of 26 environment-stable expression QTLs， including seven major QTLs.For example， cqFT-A2-3 and cqFT-A2-4 were stably expressed under the spring ecological environment， cqFT-C6-4， cqFT-C6-7， cqFT-C6-12， and cqFT-C6-13 were stably expressed in winter and semi-winter ecological regions， and cqFT-C6-14 was stably expressed in the winter ecological region.In conclusion， this study provides a foundation for fine mapping and map-based cloning of causal genes and offers guidance for the cultivation of new early maturing and high-yield rapeseed varieties in the future.Fifteen candidate genes involved in flowering regulation were identified within QTL confidence intervals，including BnaA02g12260D（RGA1）， BnaA02g15390D（AGL12）， BnaA02g16710D（LKP2） andBnaC06g19930D（NUA）， and these candidate genes are involved in the regulatory pathways， such as gibberellin， photoperiod， circadian clock， vernalization response and flower development.In conclusion，this study lays a foundation for fine mapping and map-based cloning of causal genes and offers guidance for the cultivation of new early maturing and high-yield rapeseed varieties in the future.

Key words Brassica napus; Double haploid population; Flowering trait; QTL analysis; Candidate gene

Received2021-12-30Returned 2022-08-30

Foundation item The Key Research and Research Plan of Shaanxi Province （No.2020ZDLNY04-01）; The Youth Innovation Team of Shaanxi University and Collaborative Innovation Center Project of Shaanxi Tea Industry（No.23JY001）;High Level? Talent Project of? Ankang University （No.2021AYQDZR12，No.2016AYQDZR12）.

First author ZHAO Weiguo， male， associate professor.Research area：biotechnology and molecular breeding of B.napus.E-mail：Zhaoweiguo0517@126.com

Corresponding?? author WANG Hao， male， research fellow.Research area：biotechnology and molecular breeding of B.napus.E-mail：wangzhuoyuan846@sohu.com

（責任編輯：成 敏 Responsible editor：CHENG Min）

收稿日期：2021-12-30修回日期：2022-08-30

基金項目：陜西省重點研發(fā)計劃項目（2020ZDLNY04-01）；陜西高校青年創(chuàng)新團隊和陜西省茶產業(yè)協(xié)同創(chuàng)新中心項目（23JY001）;安康學院高層次人才引進項目（2021AYQDZR12，2016AYQDZR12）。

第一作者：趙衛(wèi)國，男，副教授，研究方向為油菜生物技術與分子育種。E-mail：zhaoweiguo0517@126.com

通信作者：王 灝，男，研究員，研究方向為油菜生物技術與分子育種。E-mail：wangzhuoyuan846@sohu.com