董驥馳 楊 靖 郭 濤 陳立凱 陳志強 王 慧

?

基于高密度Bin圖譜的水稻抽穗期QTL定位

董驥馳**楊 靖**郭 濤 陳立凱 陳志強*王 慧*

華南農(nóng)業(yè)大學 / 國家植物航天育種工程技術(shù)研究中心, 廣東廣州 510642

以粳稻品種02428和秈稻品種玉針香進行雜交, 按單粒傳法連續(xù)自交10代, 得到包含192個株系的重組自交系(RIL)作圖群體。通過對兩親本重測序及RIL群體簡化基因組測序, 構(gòu)建了包含2711個Bin標記的高密度遺傳圖譜。該圖譜各染色體標記數(shù)在162~311個之間, 標記間平均物理距離為137.68 kb。將親本及192個株系分別于4個環(huán)境下采用隨機區(qū)組種植, 并記錄抽穗期。使用WinQTL Cartographer 2.5軟件的CIM分析方法, 進行抽穗期相關(guān)QTL檢測及定位。在4個環(huán)境下定位到影響抽穗期的QTL共14個, 分布于第1、第2、第3、第7、第8、第9和第10染色體。其中,和能在3個環(huán)境中被重復檢測到, 表型貢獻率分別為5.14%~11.15%和5.35%~16.97%, 分別能縮短抽穗期1.66 d和1.56 d, 具有聚合育種的應用價值。通過物理位置比對, 14個QTL中有11個與前人定位在相同或鄰近區(qū)域,和尚未見報道。經(jīng)對詳細分析, 在其染色體區(qū)間內(nèi)找到3個與抽穗期相關(guān)的注釋基因、和, 其中已被克隆。測序分析發(fā)現(xiàn), 這3個基因在兩親本間都存在差異, 可作為候選基因。

水稻; 抽穗期; Bin圖譜; QTL定位

水稻抽穗期決定了品種的種植地區(qū)與季節(jié)適應性[1], 而且與產(chǎn)量、品質(zhì)和抗逆性關(guān)系密切[2], 相關(guān)的研究對指導育種實踐、品種改良及品種推廣均具有重要意義[3]。自Yano等克隆了第一個水稻抽穗期QTL ()以來, 目前至少已有734個與水稻抽穗期相關(guān)的QTL被報道(http://www.grammene.org/), 用圖位克隆的方法至少克隆了14個抽穗期相關(guān)的QTL[4-17]。已有研究表明, 水稻有2個成花素基因()和(), 以及幾個抽穗抑制基因(如和), 它們形成的調(diào)控網(wǎng)絡精密控制著抽穗。至少有2個抽穗期調(diào)節(jié)通路控制著成花素基因表達, 分別是()通路和()通路。大量研究證實, 大多數(shù)抽穗期相關(guān)基因通過通路或通路來調(diào)節(jié)成花素基因的表達,和則不同, 它們獨立于和通路, 直接調(diào)節(jié)成花素基因的表達進而促進抽穗[18]。然而, 已定位和克隆的QTL還并不能完全解釋水稻抽穗期自然變異的形成, 從根本上理解水稻抽穗期的遺傳機理還需更多的深入研究, 從而更好地幫助適宜生育期品種的選育[19-20]。目前的研究大多基于傳統(tǒng)分子標記(如AFLP、RFLP、SSR等), 操作起來耗時耗力, 不能實現(xiàn)高通量操作[21]。用傳統(tǒng)分子標記構(gòu)建的遺傳圖譜密度低、定位區(qū)間過大, 導致QTL貢獻率估算偏高、基因克隆困難以及難以開發(fā)標記應用于聚合育種。此外, 標記密度低的圖譜上存在大量Gap, 在QTL定位中不能提供足夠的信息, 甚至一些研究最終陷入多態(tài)性分子標記缺乏而難以進一步定位的境地[22]。本研究中的Bin標記分布于整個基因組, 構(gòu)建的Bin圖譜標記密度高、位置精確, 達到精細定位要求, 可直接進行分子育種標記的開發(fā)和候選基因的篩選。本研究以秈稻玉針香(YZX)和粳稻02428衍生的RIL為作圖群體, 采用GBS技術(shù)構(gòu)建包含2711個Bin標記的高密度遺傳圖譜, 分別于4個環(huán)境對水稻抽穗期進行QTL定位, 以期鑒定一些新的、可穩(wěn)定遺傳的QTL位點作為育種資源, 并為水稻抽穗期調(diào)控基因克隆及分子標記輔助育種提供更多依據(jù)。

1 材料與方法

1.1 試驗材料

用秈稻玉針香(YZX)和粳稻02428雜交, 通過單粒傳法, 從F2代開始構(gòu)建玉針香/02428的重組自交系群體。群體包含192個株系, 基因型鑒定為F6世代, 表型調(diào)查為F8~F11世代。

1.2 田間種植與性狀考查

2016年早造(E1環(huán)境)、2016年晚造(E2環(huán)境)、2017年早造(E3環(huán)境)和2017年晚造(E4環(huán)境)于華南農(nóng)業(yè)大學試驗教學基地(23.17°N, 113.37°E)種植192個株系的RIL群體及兩親本。采用隨機區(qū)組設計, 田間管理(水、肥、病蟲害防治等)按當?shù)卮筇锍Ｒ?guī)栽培要求實施。E1環(huán)境, 于2016年3月5日播種育秧, 4月5日移栽。E2環(huán)境, 于2016年7月25日播種育秧, 8月8日移栽。E3環(huán)境, 于2017年2月28日播種育秧, 3月31日移栽。E4環(huán)境, 于2017年7月22日播種育秧, 8月7日移栽。每個株系6行, 每行6株, 行株距為20 cm × 20 cm, 均單本種植。

以單株為單位調(diào)查表型, 當單株的第1個稻穗尖露出劍葉葉鞘至少1 cm時, 記為該單株的抽穗日期, 每隔1 d調(diào)查一次。剔除異常數(shù)據(jù)后, 以株系內(nèi)全部單株從播種到抽穗所經(jīng)歷天數(shù)的均值作為該株系的抽穗期表型值。

1.3 DNA提取及高通量測序

采用本實驗室自主開發(fā)設計的高通量磁珠法制備DNA樣品。用Precellys 24研磨儀(4000′, 12 s)研磨液氮冷凍過的葉片至粉末, 迅速加入85℃預熱的LB提取液(200 mmol L–1Tris-HCl (pH 7.8), 250 mmol L–1NaCl, 25 mmol L–1ethylenediaminetetraacetic acid (EDTA), 0.5% sodium dodecyl sulfate (SDS), 2% polyvinylpyrrolidone (PVP)-40), 確保葉片粉末都浸泡在提取液中, 迅速蓋上橡膠蓋并混勻, 置于65℃水浴鍋中水浴20~25 min后, 置–20℃冰箱放置20 min。將樣品轉(zhuǎn)移至96孔深孔板, 每孔加入540 μL異丙醇和2.5 μL磁珠液(100 mg mL–1, 洛陽惠爾納米科技有限公司), 用自動核酸純化儀KingFisher Flex進行DNA分離純化。DNA最終溶于150 μL洗脫液(10 mmol L–1Tris-HCl, 1 mmol L–1EDTA (pH 8.3), 100 μg mL–1RNase A), 于–80℃保存?zhèn)溆?。以瓊脂糖凝膠電泳檢測基因組DNA的質(zhì)量, Nanodrop ND-1000微量核酸蛋白檢測儀檢測DNA的濃度與純度[23]。

采用WGS (Whole Genome Sequencing)和GBS (Genotyping-By-Sequencing)技術(shù)分別對兩親本及其衍生的RIL群體測序。水稻基因組經(jīng)電子酶切評估, 選擇酶切片段大小適宜, 酶切位點在基因組上分布均勻, 且能有效降低酶切片段中基因組重復序列所占比例的限制性內(nèi)切酶?；蚪M加上帶有barcode的接頭后, 對每個樣品進行擴增, 隨后混合選擇需要的片段構(gòu)建GBS文庫。利用Illumina HiSeq2500測序平臺, 進行雙末端(Paired-End)測序。測序后的數(shù)據(jù)分析流程如下: (1)數(shù)據(jù)質(zhì)量評估, 去除接頭、污染序列及低質(zhì)量reads; (2)統(tǒng)計酶捕獲的Reads數(shù)量; (3)參考基因組(http://plants.ensembl.org/Oryza_ sativa/)比對: 利用BWA軟件(http://bio-bwa.sourceforge. net/)將測得基因組與參考基因組比對, 統(tǒng)計比對率、覆蓋深度、基因組覆蓋率等; (4)群體SNP檢測; (5)遺傳標記開發(fā)?；谟H本SNP檢測結(jié)果, 開發(fā)子代標記, 并統(tǒng)計子代測序深度、覆蓋度等, 評估測序質(zhì)量。與親本玉針香相同的基因型記為“0”, 與親本02428相同的基因型記為“2”, 部分雜合基因型記為“1”。(6)遺傳標記過濾。檢查異常堿基, 檢查基因型缺失覆蓋率并過濾, 檢查基因型頻率, 過濾偏分離標記。

1.4 遺傳重組鑒定和連鎖圖譜構(gòu)建

采用ScmapV5軟件構(gòu)建遺傳圖譜。構(gòu)建原理是: (1)通過“滑動窗口”法(sliding window approach)尋找交換點; (2)將所有樣本中同一段序列上2個交換點之間緊密連鎖不發(fā)生重組的若干個SNP位點看做一個整體區(qū)塊, 即“Bin標記”, 以Bin的起點來確定其所在物理位置; (3)根據(jù)基因型, 用最大似然函數(shù)計算重組率; (4)用Kosambi作圖函數(shù)將重組率轉(zhuǎn)換為遺傳距離; (5)利用兩點檢驗法檢驗兩個標記間是否連鎖。(6)通過最短距離法進行聚類分析, 把兩兩連鎖的標記聚為一個類, 得到連鎖群; (7)基于三點測驗法, 對連鎖群內(nèi)的標記進行排序。

1.5 數(shù)據(jù)的統(tǒng)計分析

在Microsoft Excel 2010中完成表型數(shù)據(jù)處理和作圖。采用BWA軟件的默認設置進行基因組比對及相關(guān)測序數(shù)據(jù)統(tǒng)計分析。采用WinQTL Cartographer 2.5軟件定位水稻抽穗期QTL, 以LOD=2.5為閾值檢測QTL; 應用CIM算法, 掃描區(qū)間為1.0 cM, 遺傳背景控制選擇標準模型(模型6)。QTL定位結(jié)果以LOD的峰值作為該QTL的LOD值, 以LOD峰值位置的Bin標記來估計QTL效應, 以LOD值下降1.0的區(qū)域作為QTL的置信區(qū)間, 遵循McCouch等[24]的原則命名QTL。應用SPSS 19.0軟件進行方差分析及計算。應用Origin 8.0進行頻次分布圖的繪制。

2 結(jié)果與分析

2.1 遺傳連鎖圖譜構(gòu)建

親本玉針香和02428重測序Q20的比例分別為93.72%和94.22%, Q30比例分別為86.44%和87.47%。兩親本測序Reads與Nipponbare (ssp)參考基因組的比對率分別為96.29%和97.97%, 平均深度分別達到20.26×和26.21×, 1×以上覆蓋度為90.60%和96.36%, 4×以上覆蓋度為85.41%和94.53%。兩親本中共發(fā)現(xiàn)1 534 036個多態(tài)性SNP, 篩選后剩余aa×bb型1 334 454個。

利用GBS技術(shù)對192個RIL株系進行測序, 測序深度平均值為11.76×, 4×以上覆蓋度為7%。分別提取192個子代在上述的1 334 454個親本多態(tài)性標記位點的基因型, 經(jīng)過分析和篩選最終共獲得2711個均勻分布于12條染色體的Bin標記。將Bin標記錨定于連鎖群后, 各染色體上標記數(shù)量為162~311個, 各染色體的平均遺傳距離為195.3 cM, 兩標記間遺傳距離平均值為0.86 cM, 兩標記間物理距離平均值為137.36 kb, 標記整體的分布達到精細作圖的要求(圖1)。

2.2 抽穗期表型特征

方差分析表明兩親本抽穗期天數(shù)的均值存在極顯著差異, 4個環(huán)境下玉針香的抽穗期均顯著長于02428 (表1)。

表2表明, 在E1環(huán)境中, RIL群體抽穗期的平均值為95.65 d, 抽穗期最短的株系為81.33 d, 抽穗期最長的株系為117.83 d; 在E2環(huán)境中, RIL群體抽穗期的平均值為69.15 d, 抽穗期最短的株系為53.13 d, 抽穗期最長的株系為85.4 d; 在E3環(huán)境中, RIL群體抽穗期的平均值為95.29 d, 抽穗期最短的株系為81 d, 抽穗期最長的株系為110.3 d; 在E4環(huán)境中, RIL群體抽穗期的平均值為70.54 d, 抽穗期最短的株系為61.3 d, 抽穗期最長的株系為87.5 d; 四個環(huán)境中RIL群體的抽穗期峰度和偏度絕對值都小于1, 總體呈正態(tài)分布(圖2), 表現(xiàn)為數(shù)量性狀遺傳模式。

2.3 抽穗期QTL定位

在4個環(huán)境中共檢測到14個抽穗期QTL(表3和圖1)。其中, E1環(huán)境中檢測到9個, 分布于第2、第3、第7、第8和第10染色體, LOD值介于2.9~5.6之間, 貢獻率介于4.7%~9.5%之間; E2環(huán)境中檢測到5個, 分布于第2、第7和第10染色體, LOD值介于3.0~5.0之間, 貢獻率介于5.2%~9.5%之間; E3環(huán)境中檢測到3個, 分布于第1、第2和第8染色體, LOD值介于2.5~5.9之間, 貢獻率介于4.5%~11.2%之間; E4環(huán)境中檢測到5個, 分布于第3、第9和第10染色體, LOD值介于2.8~9.6之間, 貢獻率介于4.6%~17.0%之間。其中, 有2個QTL在3個環(huán)境中均檢測到。

圖1 遺傳標記在染色體上的分布及QTL位置

每條染色體上的黑色線條代表Bin標記所在位置; 紅色字體為已克隆QTL; 藍色字體為新發(fā)現(xiàn)的QTL。

Black lines represent the positions of Bin markers on each linkage group; fonts in red are cloned QTLs; fonts in blue are novel QTLs.

表1 02428與玉針香抽穗期天數(shù)比較

**表示差異在0.01水平達到極顯著。

**indicates significant difference at the 0.01 probability level.

表2 RIL群體抽穗期性狀表現(xiàn)

圖2 玉針香/02428重組自交系及親本的抽穗期表型分布特征

在第1染色體上, 僅在E3環(huán)境中被檢測到, 位于160.17 cM, 與mk186連鎖。其LOD值為2.5, 置信區(qū)間為25.95~26.05 Mb, 表型貢獻率為4.5%, 來自于02428的等位基因可縮短生育期1.3 d。

位于第2染色體, 僅在E2環(huán)境中被檢測到, 在173.61 cM LOD達到峰值3.2, 與mk484緊密連鎖, 置信區(qū)間為17.60~28.35 Mb, 能解釋5.8%的表型變異, 來自于02428的等位基因可縮短抽穗期1.3 d。同一染色體上的在3個環(huán)境中均被檢測到, 表現(xiàn)出較強的穩(wěn)定性。在E1環(huán)境檢測到的位置是178.21 cM, LOD值為3.1, 置信區(qū)間為28.25~28.65 Mb, 可解釋5.1%表型變異, 來自于02428的等位基因可縮短抽穗期1.5 d; 在E2環(huán)境檢測到的位置是179.71 cM, LOD值為5.0, 置信區(qū)間為28.45~28.65 Mb, 可解釋9.5%表型變異, 來自于02428的等位基因可縮短抽穗期1.6 d。在E3環(huán)境檢測到的位置是175.19 cM, LOD值為5.9, 置信區(qū)間為28.15~28.65 Mb, 可解釋11.2%的表型變異, 來自于02428的等位基因可縮短抽穗期1.9 d。

第3染色體上共有3個抽穗期QTL, 分別為和。在E1環(huán)境中被檢測到, 位于194.21 cM處, LOD值為2.9, 置信區(qū)間為30.05~30.65 Mb, 可解釋4.7%表型變異, 來自于02428的等位基因可縮短抽穗期1.4 d; 在E4環(huán)境中被檢測到, 位于197.3 cM處, LOD值為3.4, 置信區(qū)間為30.45~31.05 Mb, 可解釋5.6%表型變異, 來自于02428的等位基因可縮短抽穗期1.0 d。在E1環(huán)境被檢測到, 位于205.11 cM處, LOD值為3.2, 置信區(qū)間為31.85~32.45 Mb, 可解釋5.25%表型變異, 來自于02428的等位基因可縮短抽穗期1.5 d; 在E4環(huán)境被檢測到, 位于205.11 cM處, LOD值為3.0, 置信區(qū)間為31.85~32.45 Mb, 可解釋4.85%表型變異, 來自于02428的等位基因可縮短抽穗期1.0 d。僅在E1環(huán)境中被檢測到, 位于211.21 cM處, LOD值為3.2, 置信區(qū)間為33.10~33.55 Mb, 可解釋5.3%表型變異, 來自于02428的等位基因可縮短抽穗期1.5 d。

第7染色體上共有2個抽穗期QTL。只在E1環(huán)境中被檢測到, 位于48.81 cM處, 與mk1681緊密連鎖, LOD值為4.2, 置信區(qū)間為8.65~8.85 Mb, 可解釋7.0%的表型變異, 來自于02428的等位基因可推遲抽穗期1.7 d。只在E2環(huán)境中被檢測到, 位于118.41 cM處, 與mk1725緊密連鎖, LOD值為3.0, 置信區(qū)間為22.05~22.55 Mb, 可解釋5.2%的表型變異, 來自于02428的等位基因可縮短抽穗期1.2 d。第8染色體上共有2個抽穗期QTL, 即和。在E1環(huán)境檢測到的位置是42.41 cM, 與mk1823緊密連鎖, LOD值為5.6, 置信區(qū)間為4.20~4.55 Mb, 可解釋9.5%表型變異, 來自于02428的等位基因能推遲抽穗期2 d; 在E3環(huán)境檢測到的位置是42.44 cM, 與mk1823緊密連鎖, LOD值為3.9, 置信區(qū)間為4.20~4.55 Mb, 可解釋7.5%表型變異, 來自于02428的等位基因能推遲抽穗期2.4 d。僅在E1環(huán)境被檢測到, 位于56.71 cM處, 與mk1830緊密連鎖, LOD值為4.3, 置信區(qū)間為4.65~5.15 Mb, 可解釋7.4%表型變異, 來自于02428的等位基因能推遲抽穗期1.7 d。

在第9染色體上125.4 cM處, 僅在E4環(huán)境中被檢測到, 與mk2090連鎖。其LOD值為2.8, 置信區(qū)間為16.55~16.75 Mb, 表型貢獻率為4.6%, 來自于02428的等位基因可縮短生育期1.0 d。

第10染色體56.71 cM處為, 僅在E1環(huán)境被檢測到, 它與mk2245緊密連鎖, LOD值為3.5, 置信區(qū)間為16.35~16.65 Mb, 可解釋6.7%表型變異, 來自于02428的等位基因可縮短抽穗期1.6 d。同一染色體上的在E1、E2和E4環(huán)境被重復檢測到, 置信區(qū)間為16.75~17.25 Mb, E1環(huán)境的LOD值為4.5, 貢獻率為8.2%, 來自于02428的等位基因可縮短抽穗期1.4 d; E2環(huán)境的LOD值為3.2, 貢獻率為5.4%, 來自于02428的等位基因可縮短抽穗期1.5 d; E4環(huán)境的LOD值為9.6, 貢獻率為17.0%, 來自于02428的等位基因可縮短抽穗期1.7 d。第10染色體67.41 cM處是, 僅在E2環(huán)境中被檢測到, LOD值為4.2, 可解釋7.5%表型變異, 置信區(qū)間為17.75~17.95 Mb, 來自于02428的等位基因可縮短抽穗期1.4 d。

2.4 注釋基因篩選

在3個環(huán)境中可重復檢測到, 其置信區(qū)間最大時為28.15~28.65 Mb, 在該染色體區(qū)間有33個注釋基因(The Rice Annotation Project Database), 逐一篩除發(fā)現(xiàn)有3個位點很可能影響抽穗期, 分別是、和其中,與擬南芥中()同源, 可影響水稻抽穗期、胚的發(fā)育及小花器官數(shù)目。和的基因產(chǎn)物是可表達的蛋白, 其涉及的生物學進程包括花器官的發(fā)育、生殖生長等。

DNA雙向測序發(fā)現(xiàn), 玉針香和02428的3個候選基因間均有序列差異(圖3)。玉針香與02428的基因存在23處堿基差異, 共18處為單堿基的轉(zhuǎn)換, 5處為3堿基以內(nèi)的插入缺失; 其中7處在CDS區(qū)變異, 有5處發(fā)生錯義突變, 分別位于CDS+830 bp (TTT/TCT)、CDS+1008 bp (GAG/GAT)、CDS+150 bp (TTG/TTT)、CDS+2413 bp (CGT/TGT)和CDS+2 418 bp (GAC/GAG)。玉針香與02428的基因存在23處堿基差異, 其中21處為單堿基的轉(zhuǎn)換, 2處為單堿基插入缺失; 其中2處在CDS區(qū)的變異, 有1處為錯義突變, 位置是CDS+394 bp (AAT/TAT)。玉針香與02428的基因存在13處堿基差異, 其中9處為單堿基的轉(zhuǎn)換, 4處為插入缺失; 其中2處在CDS區(qū)變異, 有1處為錯義突變, 位置是CDS+883 bp (CTT/TTT)。同時, 測序比對發(fā)現(xiàn), 02428中3個候選基因的CDS序列與日本晴參考序列完全相同。

除此之外, 還在其他QTL的置信區(qū)間內(nèi)找到一些與抽穗期相關(guān)可能性較大的注釋基因, 它們都涉及到花的發(fā)育過程(表4)。

表3 4個環(huán)境下檢測到的抽穗期QTL

1)指來自于02428的等位基因造成的加性效應, 正值為延長抽穗期, 負值為縮短抽穗期。

1)Additive effect from 02428’s allele. Positive (negative) values mean increasing (decreasing) heading date.

表4 QTL置信區(qū)間內(nèi)注釋基因篩選

1)Source: CHINA RICE DATA CENTER and The Rice Annotation Project Database.

圖3 玉針香與02428之間候選基因的結(jié)構(gòu)和變異

A:; B:; C:。黑色框: 外顯子; 灰色部分: 編碼序列; 紅色箭頭: SNP; 帶黑點紅色箭頭: 錯義突變; 藍色箭頭: 插入缺失。

A:; B:; C:Frames with black lines: exon; Grey boxes: protein coding sequence; Red arrow: SNP; Red arrow with a black point: missense mutation; Blue arrow: InDel.

3 討論

3.1 采用Bin圖譜定位QTL的優(yōu)勢

目前大多數(shù)研究者采用傳統(tǒng)分子標記(如SSR、RFLP等)構(gòu)建遺傳圖譜[25], 精度約為1~10 Mb, 存在標記數(shù)量少且分布不均勻的問題, 影響了QTL定位的精確度, 另外標記密度較小還可能導致某些雙交換位點的漏測[26]。Bin圖譜與傳統(tǒng)的遺傳圖譜相比, 標記密度高, 精度達100 kb, 且一個Bin內(nèi)部包含多個不發(fā)生重組的SNP, 因而雙交換能夠被精確檢測。同時, Bin圖譜是基于測序技術(shù)構(gòu)建的, 能提供準確的物理位置, 可使遺傳分析和QTL定位更準確。此外, 測序分型得到的大量SNP標記可快速應用于分子育種。本研究采用GBS技術(shù)構(gòu)建了以SNP為基礎(chǔ)的高分辨率Bin圖譜, 包含2711個Bin標記, 標記間平均物理距離137.36 kb, 標記整體的分布達到精細作圖的密度, 可直接從定位區(qū)間篩選候選基因。另外, 定位結(jié)果顯示, 在LOD值超過閾值的一段區(qū)域一般存在多個峰值, 可定位到多個QTL, 如位置相鄰的和以及和等, 說明Bin圖譜對QTL的檢測更精細, 可將位置相鄰的QTL有效地分離、解析。

3.2 本研究與前人定位QTL的比較

本研究定位的在第1染色體上, 僅在E3環(huán)境中被檢測到, 其LOD值為2.5, 置信區(qū)間為25.95~26.05 Mb, 表型貢獻率為4.5%, 在前人的研究中未見報道。

位于第2染色體27.95 Mb處, 置信區(qū)間為17.60~28.35 Mb, 表型貢獻率為5.8%, LOD值為3.2, 在該區(qū)間有多篇相關(guān)的報道[27-32], 其中Thomson等[31]報道的表型貢獻率為4.4%, LOD為4.06, 在5個環(huán)境中檢測到1次, 與位置也最近, 它們可能為同一QTL。在第2染色體上, E1環(huán)境于28.55 Mb位置檢測到1個QTL, 其置信區(qū)間為28.25~28.65 Mb; E2環(huán)境于28.65 Mb位置檢測到1個QTL, 置信區(qū)間為28.45~28.65 Mb; E3環(huán)境于28.35 Mb位置檢測到1個QTL, 置信區(qū)間為28.15~28.65 Mb。這3個QTL位置十分相近, 置信區(qū)間十分相似, 因此視為同一QTL, 命名為, 它在前人的研究中未見報道, 很可能是與水稻抽穗期相關(guān)的新位點。

第3染色體上檢測到3個QTL和。在2個環(huán)境中被重復檢測到, 在E1環(huán)境定位的位置是30.55 Mb, 置信區(qū)間分別為30.55~30.65 Mb; E4環(huán)境定位的位置是30.95 Mb, 置信區(qū)間為30.45~31.05 Mb。也在2個環(huán)境中被重復檢測到, 定位的位置是32.35 Mb, 置信區(qū)間為31.85~32.45 Mb。僅在E1環(huán)境被檢測到, 位于33.35 Mb, 置信區(qū)間為33.10~ 33.55 Mb。這3個QTL所在區(qū)域已有相關(guān)報道[31,33-34]。

位于第7染色體8.85 Mb處, 置信區(qū)間為8.65~8.85 Mb, 已精細定位的包含了該區(qū)間[35]。同一染色體上的位于22.35 Mb處, 置信區(qū)間為22.05~22.55 Mb, 該區(qū)間位于Thomson等[31]定位的所在區(qū)間RM125–RM336內(nèi)部, Jiang等[36]定位的與的區(qū)間也有部分重疊。

位于第8染色體4.45 Mb處, 置信區(qū)間為4.20~ 4.55 Mb, 與已克隆的[39]位置相吻合, 都表現(xiàn)為長日照延長抽穗期。位于第8染色體5.15 Mb處, 置信區(qū)間為4.65~5.15 Mb, 在該區(qū)段也有多個相關(guān)報道[37-40]。

僅在E4環(huán)境被檢測到, 位于第9染色體16.75 Mb處, 與mk2090緊密連鎖, 置信區(qū)間為16.55~16.75 Mb,未發(fā)現(xiàn)相關(guān)報道。

位于第10染色體16.45 Mb處, 置信區(qū)間為16.35~16.65 Mb, 與Thomson等[31]定位的區(qū)間有部分相重疊。在第10染色體上, E1環(huán)境于17.25 Mb位置檢測到1個QTL, 其置信區(qū)間為16.75~17.25 Mb; E2環(huán)境于16.75 Mb位置檢測到1個QTL, 置信區(qū)間為16.75~ 17.25 Mb。這2個QTL位置相近、區(qū)間相同, 因此視為同一QTL, 命名為, Doi等[7]報道的的位置在其區(qū)間內(nèi)。同在第10染色體上的位于17.85 Mb處, 置信區(qū)間為17.75~17.95 Mb, Mei等[32]報道了與抽穗期相關(guān)的RG241a–CDO98染色體區(qū)段在的區(qū)間內(nèi), Zhou等[28]報道的與區(qū)間基本相同。

3.3 檢測到的QTL應用價值及研究價值

本研究在4個環(huán)境中共檢測到14個影響抽穗期的QTL, 分布于第1、第2、第3、第7、第8、第9和第10染色體上。除、、、和能在2個及以上環(huán)境被檢測到外, 其他QTL只在單個環(huán)境被檢測到。其中,和所在的位置分別是克隆基因和。

Hori等[20]同時用12個群體進行抽穗期相關(guān)的QTL定位, 在其中3個群體中能檢測到, 表型貢獻率為7.9%~74.0%; 在其中2個群體中能檢測到, 表型貢獻率為18.2%~18.3%。Cheng等[2]同時用3個群體進行抽穗期相關(guān)的QTL定位, 在其中1個群體中能檢測到, 表型貢獻率為3.5%; 在3個群體中均能檢測到, 表型貢獻率為3.6%~33.7%。而本研究中檢測到的表型貢獻率為7.5%~9.5%;的表型貢獻率為5.4%~ 17.0%。這些研究表明, 對于不同群體, 由于材料遺傳背景影響, QTL貢獻率的變化很大; 其次, 圖譜的精度對QTL貢獻率的估算影響也較大, 傳統(tǒng)分子標記構(gòu)建的圖譜精度通常約為4 Mb, 兩標記之間經(jīng)常會存在多個相關(guān)基因或QTL簇, 粗定位到的QTL貢獻率會較高。因此我們認為, 重點關(guān)注QTL在不同環(huán)境或不同群體中的重復定位情況對聚合育種的意義更大。

和在定位到的QTL中表型貢獻率相對較小, 表達也較不穩(wěn)定, 只能2次被重復檢測到, 因而利用的可能性小。在兩年的早造中被重復檢測到, 僅在早造環(huán)境特異表達, 表現(xiàn)為延長生育期, 這些都與在長日條件下延長抽穗期的特性相對應, 它對于將早熟品種改良成半晚熟的高產(chǎn)品種具有一定價值。和都能在3個環(huán)境中被重復檢測到, 且都能夠縮短生育期, 說明它們受環(huán)境影響較小、表達較穩(wěn)定, 對光照和溫度變化不敏感, 它們的表型貢獻率最低時分別為5.1%和5.4%, 表型貢獻率最高時分別為11.2%和17.0%, 因此對于QTL聚合改良品種生育期具有重要價值。

此外,的置信區(qū)間在3次定位中略有變化, 在E2環(huán)境判定的置信區(qū)間較其他2個環(huán)境中小,不在此置信區(qū)間內(nèi), 但此時對表型的貢獻率仍然達到9.5%。因此和更有可能是導致抽穗期表型變異的新QTL, 具體情況還需深入研究。

[1] Fujino K, Sekiguchi H. Mapping of quantitative trait loci controlling heading date among rice cultivars in the northernmost region of Japan., 2008, 58: 367–373

[2] Cheng L R, Wang J M, Ye G, Luo C G, Xu J L, Li Z K. Identification of stably expressed QTL for heading date using reciprocal introgression line and recombinant inbred line populations in rice., 2012, 94: 245–253

[3] 戴高興, 楊占烈, 鄧國富, 張迎信, 王會民, 翟榮榮, 曹立勇, 程式華. 超級雜交稻協(xié)優(yōu)9308重組自交系主莖葉片數(shù)的動態(tài)QTL分析. 中國水稻科學, 2012, 26: 291–296 Dai G X, Yang Z L, Deng G F, Zhang Y X, Wang H M, Zhai R R, Cao L Y, Cheng S H. QTL analysis for leaf number on main stem in RILs of super hybrid rice Xieyou 9308.i, 2012, 26: 291–296 (in Chinese with English abstract)

[4] Yano M, Katayose Y, Ashikari M, Yamanouchi U, Monna L, Fuse T, Baba T, Yamamoto K, Umehara Y, Nagamura Y. Hd1, a major photoperiod sensitivity quantitative trait locus in rice, is closely related to the Arabidopsis flowering time gene., 2000, 12: 2473–2484

[5] Takahashi Y, Shomura A, Sasaki T, Yano M. Hd6, a rice quantitative trait locus involved in photoperiod sensitivity, encodes the α subunit of protein kinase CK2., 2001, 98: 7922–7927

[6] Kojima S, Takahashi Y, Kobayashi Y, Monna L, Sasaki T, Araki T, Yano M. Hd3a, a rice ortholog of the Arabidopsis FT gene, promotes transition to flowering downstream of Hd1 under short-day conditions., 2002, 43: 1096–1105

[7] Doi K, Izawa T, Fuse T, Yamanouchi U, Kubo T, Shimatani Z, Yano M, Yoshimura A., a B-type response regulator in rice, confers short-day promotion of flowering and controlsgene expression independently of Hd1., 2004, 18: 926–936

[8] Xue W, Xing Y, Weng X, Zhao Y, Tang W, Wang L, Zhou H, Yu S, Xu C, Li X, Zhang Q. Natural variation inis an important regulator of heading date and yield potential in rice., 2008, 40: 761–767

[9] Wei X, Xu J, Guo H, Jiang L, Chen S, Yu C, Zhou Z, Hu P, Zhai H, Wan J. DTH8 suppresses flowering in rice, influencing plant height and yield potential simultaneously., 2010, 153: 1747–1758

[10] Bian X F, Liu X, Zhao Z G, Jiang L, Gao H, Zhang Y H, Zheng M, Chen L M, Liu S J, Zhai H Q. Heading date gene,controlled late flowering inSteud. by down-regulating., 2011, 30: 2243–2254

[11] Yan W H, Wang P, Chen H X, Zhou H J, Li Q P, Wang C R, Ding Z H, Zhang Y S, Yu S B, Xing Y Z. A major QTL, Ghd8, plays pleiotropic roles in regulating grain productivity, plant height and heading date in rice., 2011, 4: 319–330

[12] Matsubara K, Ogisotanaka E, Hori K, Ebana K, Ando T, Yano M. Natural variation in Hd17, a homolog of Arabidopsis ELF3 that is involved in rice photoperiodic flowering., 2012, 53: 709–716

[13] Gao H, Zheng X, Fei G, Chen J, Jin M, Ren Y, Wu W, Zhou K, Sheng P, Zhou F, Jiang L, Wang J, Zhang X, Guo X, Wang J, Cheng Z, Wu C, Wang H, Wan J. Ehd4 encodes a novel and Oryza-Genus-Specific regulator of photoperiodic flowering in rice.2013, 9: e1003281

[14] Hori K, Ogisotanaka E, Matsubara K, Yamanouchi U, Ebana K, Yano M., a gene for casein kinase I, is involved in the control of rice flowering time by modulating the day-length response., 2013, 76: 36–46

[15] Koo B H, Yoo S C, Park J W, Kwon C T, Lee B D, An G, Zhang Z Y, Li J J, Li Z C, Paek N C. Natural variation in OsPRR37 regulates heading date and contributes to rice cultivation at a wide range of latitudes., 2013, 6: 1877–1888

[16] Ogisotanaka E, Matsubara K, Yamamoto S, Nonoue Y, Wu J, Fujisawa H, Ishikubo H, Tanaka T, Ando T, Matsumoto T. Natural variation of thecontributes to flowering time divergence in rice., 2013, 8: e75959

[17] Wu W, Zheng X M, Lu G, Zhong Z, Gao H, Chen L, Wu C, Wang H J, Wang Q, Zhou K, Wang J, Wu F, Zhang X, Guo X, Cheng Z, Lei C, Lin Q, Jiang L, Wang H, Ge S, Wan J. Association of functional nucleotide polymorphisms at DTH2 with the northward expansion of rice cultivation in Asia., 2013, 110: 2775–2780

[18] Sun B, Zhan X D, Lin Z C, Wu W X, Yu P, Zhang Y X, Sun L P, Cao L Y, Cheng S H. Fine mapping and candidate gene analysis of qHD5, a novel major QTL with pleiotropism for yield-related traits in rice (L.).2016, 130: 1–12

[19] Shen G, Xing Y. Two novel QTLs for heading date are identified using a set of chromosome segment substitution lines in rice (L.)., 2014, 41: 659–662

[20] Hori K, Nonoue Y, Ono N, Shibaya T, Ebana K, Matsubara K, Ogisotanaka E, Tanabata T, Sugimoto K, Taguchishiobara F. Genetic architecture of variation in heading date among Asian rice accessions., 2015, 15: 115

[21] Xie W, Feng Q, Yu H, Huang X, Zhao Q, Xing Y, Yu S, Han B, Zhang Q. Parent-independent genotyping for constructing an ultrahigh-density linkage map based on population sequencing., 2010, 107: 10578–10583

[22] Yu H, Xie W, Wang J, Xing Y, Xu C, Li X, Xiao J, Zhang Q. Gains in QTL detection using an ultra-high density SNP map based on population sequencing relative to traditional RFLP/SSR markers.e, 2012, 6: e17595

[23] Chen L, Gao W, Guo T, Huang C, Huang M, Wang J, Xiao W, Yang G, Liu Y, Wang H, Chen Z. A genotyping platform assembled with high-throughput DNA extraction, codominant functional markers, and automated CE system to accelerate marker-assisted improvement of rice., 2016, 36: 123, https://doi.org/10.1007/s11032-016-0547-y

[24] McCouch S R. Gene nomenclature system for rice., 2008, 1: 72–84

[25] Golicz A A, Bayer P E, Edwards D. Skim-based genotyping by sequencing., 2015, 1245: 257–270

[26] Yu H, Xie W, Li J, Zhou F, Zhang Q. A whole-genome SNP array (RICE6K) for genomic breeding in rice.2014, 12: 28–37

[27] Zou J H, Pan X B, Chen Z X, Xu J Y, Lu J F, Zhai W X, Zhu L H. Mapping quantitative trait loci controlling sheath blight resistance in two rice cultivars (L.).2000, 101: 569–573

[28] Zhou Y, Li W, Wu W, Chen Q, Mao D, Worland A J. Genetic dissection of heading time and its components in rice., 2001, 102: 1236–1242

[29] Li Z K, Yu S B, Lafitte H R, Huang N, Courtois B, Hittalmani S, Vijayakumar C H M, Liu G F, Wang G C, Shashidhar H E, Zhuang J Y, Zheng K L, Singh V P, Sidhu J S, Srivantaneeyakul S, Khush G S. QTL×environment interactions in rice: I. Heading date and plant height., 2003, 108: 141–153

[30] Xiao J, Li J, Yuan L, Tanksley S D. Identification of QTLs affecting traits of agronomic importance in a recombinant inbred population derived from a subspecific rice cross., 1996, 92: 230–244

[31] Thomson M J, Tai T H, Mcclung A M, Lai X H, Hinga M E, Lobos K B, Xu Y, Martinez C P, McCouch S R. Mapping quantitative trait loci for yield, yield components and morphological traits in an advanced backcross population betweenand thecultivar Jefferson., 2003, 107: 479–493

[32] Mei H W, Luo L J, Ying C S, Wang Y P, Yu X Q, Guo L B, Paterson A H, Li Z K. Gene actions of QTLs affecting several agronomic traits resolved in a recombinant inbred rice population and two testcross populations., 2003, 107: 89–101

[33] Takeuchi Y, Hayasaka H, Chiba B, Tanaka I, Shimano T, Yamagishi M, Nagano K, Sasaki T, Yano M. Mapping quantitative trait loci controlling cool-temperature tolerance at booting stage in temperate japonica rice., 2001, 51: 191–197

[34] Sarma R N, Gill B S, Sasaki T, Galiba G, Sutka J, Laurie D A, Snape J W. Comparative mapping of the wheat chromosome 5A Vrn-A1 region with rice and its relationship to QTL for flowering time., 1998, 97: 103–109

[35] Lin H, Liang Z W, Sasaki T, Yano M. Fine mapping and characterization of quantitative trait loci Hd4 and Hd5 controlling heading date in rice., 2003, 53: 51–59

[36] Jiang L, Xu J, Wei X, Wang S, Tang J, Zhai H, Wan J. The inheritance of early heading in the rice variety USSR5., 2007, 34: 46–55

[37] 國廣泰史, 錢前, 佐藤宏之, 滕勝, 曾大力, 藤本寬, 朱立煌. 水稻紋枯病抗性QTL分析. 遺傳學報, 2002, 29: 50–55 Kunihiro Y, Qian Q, Sato H, Teng S, Zeng D L, Fujimoto K, Zhu L H. QTL analysis of sheath blight resistance in rice (L.)., 2002, 29: 50–55 (in Chinese with English abstract)

[38] Wang C M, Yasui H, Yoshimura A, Wan J M, Zhai H Q. Identification of quantitative trait loci controlling F2sterility and heading date in rice., 2002, 29: 339–342

[39] Lu C, Shen L, Tan Z, Xu Y, He P, Chen Y, Zhu L. Comparative mapping of QTLs for agronomic traits of rice across environments using a doubled haploid population., 1996, 93: 1211–1217

[40] 譚震波, 沈利爽, 況浩池, 陸朝福, 陳英, 周開達, 朱立煌. 水稻上部節(jié)間長度等數(shù)量性狀基因的定位及其遺傳效應分析. 遺傳學報, 1996, 23: 439–446Tan Z B, Shen L S, Kuang H C, Lu C F, Chen Y, Zhou K D, Zhu L H. Identification of QTLs for lengths of the top internodes and other traits in rice and analysis of their genetic effects., 1996, 23: 439–446 (in Chinese with English abstract)

QTL Mapping for Heading Date in Rice Using High-density Bin Map

DONG Ji-Chi**, YANG Jing**, GUO Tao, CHEN Li-Kai, CHEN Zhi-Qiang*, and WANG Hui*

National Engineering Research Centre of Plant Space Breeding / South China Agricultural University, Guangzhou 510642, Guangdong, China

A recombination inbred lines (RIL) population including 192 lines derived from an inter-subspecific cross betweenrice ‘Yuzhenxiang’ andrice ‘02428’ was used in the experiment. The two parent varieties and RIL population were separately sequenced by Whole Genome Sequencing (WGS) and Genotyping-By-Sequencing (GBS) to construct a genetic linkage map with 2711 recombination Bin markers. The number of markers on 12 chromosomes ranged from 162 to 311, and the average physical distance between two markers was 137.68 kb. The WinQTL Cartographer 2.5 was used to analysis QTLs associated with heading date in four different environments. A total of 14 QTLs associated with heading date were detected on chromosomes 1, 2, 3, 7, 8, 9, and 10. Among them,and, which explained 5.14%-11.15% and 5.35%-16.97% of the total phenotypic variation for heading date separately, could be detected in three environments and shorting heading date about 1.66 days and 1.56 days on average. Thetwo QTLs could inherit stablely, having a good potential to be applied in QTL pyramiding. After comparing the physical positions of these QTLs with those previously reported, we found 11 QTLs were located in the same or near position, among them,,andwere newly reported. Furthermore, we found one cloned geneand two annotated genesandin the genomic region of, might be related to heading date. DNA sequence comparison between YZX and 02428 revealed that all the three genes could be candidate genes.

rice; heading date; Bin map; QTL mapping

2017-12-14;

2018-03-25;

2018-04-16.

10.3724/SP.J.1006.2018.00938

王慧, E-mail: wanghui@scau.edu.cn; 陳志強, E-mail: chenlin@scau.edu.cn

**同等貢獻(Contributed equally to this work)

董驥馳, E-mail: Johnkeatsno.1@gmail.com; 楊靖, E-mail: 1349643559@qq.com

URL: http://kns.cnki.net/kcms/detail/11.1809.S.20180416.0843.004.html

本研究由國家現(xiàn)代農(nóng)業(yè)產(chǎn)業(yè)技術(shù)體系建設專項(CARS-01-12), 國家重點研發(fā)計劃項目(2016YFD0102102)和廣東省應用型研發(fā)項目(2015B020231011)資助。

This study was supported by the China Agriculture Research System (CARS-01-12), the National Key Research and Development Program of China (2016YFD0102102), and the Research and Development Program for Application in Guangdong Province (2015B020231011).