• <tr id="yyy80"></tr>
  • <sup id="yyy80"></sup>
  • <tfoot id="yyy80"><noscript id="yyy80"></noscript></tfoot>
  • 99热精品在线国产_美女午夜性视频免费_国产精品国产高清国产av_av欧美777_自拍偷自拍亚洲精品老妇_亚洲熟女精品中文字幕_www日本黄色视频网_国产精品野战在线观看 ?

    基于高密度SNP遺傳圖譜的粳稻芽期耐低溫QTL鑒定

    2020-07-21 08:07:20姜樹坤王立志楊賢莉母偉杰董世晨車韋才李忠杰遲力勇李明賢張喜娟李文華
    作物學(xué)報 2020年8期
    關(guān)鍵詞:冷性自交系染色體

    姜樹坤 王立志 楊賢莉 李 波 母偉杰 董世晨 車韋才 李忠杰 遲力勇 李明賢 張喜娟 姜 輝 李 銳 趙 茜 李文華,*

    基于高密度SNP遺傳圖譜的粳稻芽期耐低溫QTL鑒定

    姜樹坤1,*王立志1楊賢莉1李 波2母偉杰3董世晨3車韋才3李忠杰1遲力勇1李明賢1張喜娟1姜 輝2李 銳1趙 茜1李文華2,*

    1黑龍江省農(nóng)業(yè)科學(xué)院耕作栽培研究所/ 黑龍江省寒地作物生理生態(tài)重點實驗室/ 黑龍江省農(nóng)作物低溫冷害工程技術(shù)研究中心, 黑龍江哈爾濱 150086;2黑龍江省農(nóng)業(yè)科學(xué)院, 黑龍江哈爾濱 150086;3哈爾濱師范大學(xué)生命科學(xué)與技術(shù)學(xué)院, 黑龍江哈爾濱 150025

    水稻直播由于省時、省工和節(jié)約成本而備受農(nóng)戶關(guān)注。然而, 芽期耐冷性不強致使現(xiàn)行推廣的許多優(yōu)良水稻品種不適于直播生產(chǎn)。因此, 挖掘鑒定芽期耐冷位點, 為后續(xù)的輔助育種提供基因資源就日益受到重視。本研究利用麗江新團黑谷和沈農(nóng)265構(gòu)建的重組自交系群體及其重測序構(gòu)建的包含2818個bin標(biāo)記的遺傳圖譜對水稻芽期的耐冷性進行QTL定位分析。共檢測到5個芽期耐冷QTL, 分布在水稻的1號、3號、9號和11號染色體上, 增效等位基因均來自耐冷親本麗江新團黑谷。這些QTL的LOD值的范圍從3.05到24.01, 表型貢獻率為8.0%~53.5%。其中, 表型貢獻率最大的主效QTL是, 位于11號染色體長臂端的21.24 Mb~22.03 Mb之間, 物理圖譜區(qū)間為790 kb。隨后利用“選擇作圖”的策略進行了QTL驗證和累加效應(yīng)分析, 明確了可以通過QTL的累加聚合實現(xiàn)芽期耐冷能力的遺傳改良, 聚合的增效QTL越多, 耐冷能力提升越明顯。上述研究結(jié)果不僅可以增強人們對芽期水稻耐冷能力遺傳基礎(chǔ)的認識和理解, 也可以為后續(xù)直播品種的遺傳改良提供理論依據(jù)和技術(shù)指導(dǎo)。

    粳稻; 芽期耐冷性; 重測序; 遺傳圖譜; 數(shù)量性狀基因座

    作為世界上最重要的糧食作物之一, 水稻起源于熱帶、亞熱帶, 相較小麥、大麥等作物, 對低溫更加敏感[1]。冷害是全球性自然災(zāi)害, 在世界許多國家均有發(fā)生, 尤以澳大利亞南部、日本、巴西、韓國、朝鮮為甚, 大規(guī)模冷害的發(fā)生會造成數(shù)十上百億公斤糧食的損失[2-6]。我國多數(shù)稻作區(qū)均有低溫冷害發(fā)生, 在華南和長江中下游的雙季稻區(qū), 春季的“倒春寒”和秋季的“寒露風(fēng)”常常引起低溫冷害[7]。東北地區(qū)平均3~4年就會遭遇一次大范圍的冷害, 造成水稻大量減產(chǎn), 危害糧食安全[8]。水稻萌發(fā)期、芽期、苗期、孕穗期、開花期和灌漿期易受低溫冷害影響[9-10], 與產(chǎn)量有直接或間接關(guān)系。20世紀90年代以來, 我國的水稻生產(chǎn)發(fā)生了較大的變化。一方面, 水稻生產(chǎn)目標(biāo)從追求產(chǎn)量向降低成本、提高經(jīng)濟效益轉(zhuǎn)變; 另一方面, 隨著我國經(jīng)濟的快速發(fā)展, 農(nóng)村主要勞動力向城鎮(zhèn)轉(zhuǎn)移, 農(nóng)村勞動力緊缺已成為水稻生產(chǎn)發(fā)展的主要限制因素。在這樣的背景下, 水稻直播由于省時、省工和節(jié)約成本, 深得農(nóng)民的歡迎, 并迅速推廣。據(jù)不完全統(tǒng)計, 水稻直播面積在安徽省已達56.7萬公頃, 湖北省為49萬公頃, 江西省為43.3萬公頃, 江蘇省為69.3萬公頃, 浙江省為37.3萬公頃, 廣東省為10萬公頃, 黑龍江省的直播稻面積為40萬公頃[11-12]。華南雙季稻區(qū)早稻直播和東北高寒稻區(qū)的單季稻直播往往會遇到春季低溫天氣, 造成缺苗斷壟, 導(dǎo)致嚴重減產(chǎn)[9-12]。這就要求用于直播的水稻種子有較強的芽期耐低溫能力。而上述稻作區(qū)一直以育秧移栽為主, 芽期耐低溫能力長期不作為育種目標(biāo), 致使許多高產(chǎn)、優(yōu)質(zhì)的水稻品種由于芽期耐低溫能力弱而不能用作直播稻, 限制了直播稻的生產(chǎn)。因此, 水稻芽期耐冷性是其生長發(fā)育過程中不可忽視的重要性狀。

    自從20世紀30年代, 科研人員開始研究水稻冷害以來, 已在冷害發(fā)生的氣象學(xué)原因、水稻耐冷的生理機制、遺傳解析和基因克隆等方面積累了大量的研究成果[1,13-19]。但與萌發(fā)期、苗期和孕穗期冷害相比, 人們在水稻芽期耐冷性方面的研究稍顯不足, 目前主要開展的是芽期耐冷基因的QTL定位研究。嚴長杰等[20]利用南京11和巴利拉的DH群體, 在7號染色體上鑒定出1個芽期冷害QTL。Zhang等[21]利用Lemont和特青的重組自交系群體, 在3號、7號和11號染色體上檢測到3個控制水稻芽期耐冷性的QTL。喬永利等[22]利用密陽23和吉冷1號的F2:3群體, 在2號、4號和7號染色體上檢測到3個芽期耐冷QTL。陳瑋等[23]利用Lemont和特青的重組自交系群體, 在1號、3號、7號和11號染色體上檢測到4個控制水稻芽期耐冷性的QTL。張露霞等[24]利用Asominori和IR24的重組自交系群體和相應(yīng)的染色體片段置換系, 在5號和12號染色體上檢測到3個芽期耐冷QTL。鞏迎軍等[25]利用越光和Kasalath的回交重組自交系群體, 在4號、6號和11號染色體上檢測到4個芽期耐冷QTL。林靜等[26]利用以秈稻9311為受體、粳稻日本晴為供體的染色體片段置換系, 在5號和7號染色體上鑒定出4個芽期耐冷QTL。Baruah等[27]利用栽培稻A58和野生稻W(wǎng)107的重組自交系群體, 在1號、11號和12號染色體上檢測到了3個芽期耐冷QTL。Ji等[28]利用TN1和春江06的DH群體, 在2號、4號和11號染色體上鑒定出3個芽期耐冷QTL。周勇等[29]利用秈稻9311為受體、粳稻日本晴為供體的染色體單片段置換系, 在除1號和11號染色體以外的10條染色體上檢測出18個芽期耐冷QTL。楊洛淼等[30]利用空育131和東農(nóng)422的重組自交系群體, 在4號染色體上檢測到1個芽期耐冷QTL。朱金燕等[31]利用廣陸矮4號為受體, 日本晴為供體的染色體片段置換系, 在1號、6號、8號、9號和10號染色體上檢測到8個芽期耐冷QTL。Yang等[32]利用華粳秈74為受體, 南洋占為供體的染色體片段置換系, 在5號和6號染色體上鑒定了2個芽期耐冷QTL。Zhang等[33]利用249份秈稻種質(zhì)資源材料結(jié)合5K水稻芯片進行芽期耐冷的GWAS分析, 在1號、3號、4號、5號、6號、10號和12號染色體上共檢測到13個芽期耐冷相關(guān)區(qū)域。

    雖然人們對芽期耐冷性進行了較多的QTL鑒定分析, 而且在水稻的12條染色體上均檢測到了QTL的存在, 但鑒定到的QTL多為微效QTL, 進行精細定位和后續(xù)研究難度較大。同時, 由于研究所用的遺傳群體和試驗環(huán)境等方面的不同, 研究結(jié)果不盡相同, 因此有必要采用與以往不同的遺傳群體和更高密度的遺傳圖譜來對已報道的芽期耐冷QTL進行驗證, 并檢測新的芽期耐冷位點。本研究利用麗江新團黑谷和沈農(nóng)265構(gòu)建的重組自交系群體, 結(jié)合重測序技術(shù)構(gòu)建的高密度bin遺傳圖譜, 挖掘控制芽期耐冷性的QTL, 以期為直播稻的耐冷設(shè)計育種提供有用的“基因資源”, 同時為闡明水稻芽期耐冷的遺傳和分子機制提供參考。

    1 材料與方法

    1.1 試驗材料

    2004年, 以芽期耐冷性強的云南地方粳稻品種麗江新團黑谷(LTH)為母本, 芽期耐冷性弱的北方超級稻品種沈農(nóng)265 (SN265)為父本配制雜交組合, 2005—2013年分別在遼寧沈陽(2005—2006, 2008—2009)、海南三亞(2008)和黑龍江哈爾濱(2011—2015)三地, 采用單粒傳法(single seed descent, SSD)套袋自交11代獲得包含144個株系的穩(wěn)定遺傳重組自交系群體(recombinant inbred lines, RILs)。

    參照Huang等[34]報道的重測序手段, 采用“滑動窗口”法構(gòu)建了麗江新團黑谷和沈農(nóng)265的重組自交系群體bin遺傳圖譜。該bin連鎖圖譜共包含58,738個重組熱點, 平均每個株系405個重組熱點。最終產(chǎn)生了2818個重組bin標(biāo)記(圖1-A), 這些標(biāo)記覆蓋了大部分的重組熱點(圖1-B), 平均每條染色體上覆蓋235個bin標(biāo)記, 標(biāo)記區(qū)間平均為128.8 kb[35]。

    1.2 芽期耐冷性評價方法

    親本麗江新團黑谷(LTH)、沈農(nóng)265 (SN265)和144份重組自交系群體材料存放于種子儲藏柜中待用, 儲藏條件為溫度5°C, 相對濕度為10%。評價芽期耐冷性參照韓龍植等[36]的方法略作修改, 每個株系3次重復(fù), 利用3次的平均值進行后續(xù)分析。首先將種子從儲藏柜中取出, 在50°C恒溫干燥箱內(nèi)高溫處理48 h, 使其充分干燥并打破種子休眠。每次重復(fù)挑選飽滿粒100粒, 經(jīng)0.1%的氯化汞溶液消毒10 min, 自來水沖洗3~4次, 去離子水洗滌3次, 于墊雙層濾紙的直徑9 cm培養(yǎng)皿中, 加入10 mL去離子水。蓋上培養(yǎng)皿蓋, 在30°C恒溫培養(yǎng)箱(上海一恒MGC-250)內(nèi)催芽2~3 d。然后從恒溫箱中取出, 用去離子水沖洗1~2次。從發(fā)芽的種子中精心挑選芽長約1~2 cm的種子50粒, 種植于邊長2 cm的營養(yǎng)缽盆中, 缽盆中事先放置配制好的營養(yǎng)育苗土。28°C光照培養(yǎng)箱(上海一恒MGC-450HP)內(nèi)放置1 d后, 轉(zhuǎn)入大型人工氣候箱(日本三江平原項目援建霜凍害實驗室)內(nèi)處理48 h, 溫度條件為2°C。處理后, 將試驗材料轉(zhuǎn)移到溫度為20~30°C, 并有陽光的地方, 使其恢復(fù)生長, 每天適量澆水。恢復(fù)10 d后, 調(diào)查秧苗生長狀態(tài)。按照圖2的標(biāo)準(zhǔn)給處理后的秧苗恢復(fù)狀態(tài)賦值, 共分為5級: 恢復(fù)后幼芽枯死的賦值為0; 恢復(fù)后幼芽停止發(fā)育, 出現(xiàn)畸形的賦值為1; 恢復(fù)后有一定生長能力, 葉色較淡, 但發(fā)育非常緩慢的賦值為2; 恢復(fù)后生長能力較強, 葉色淡, 發(fā)育較緩慢的賦值為3; 恢復(fù)后葉色較深, 能夠快速生長的賦值為4。并按照不同級別進行加權(quán)計分。芽期耐冷級別=(0×0級的株數(shù)+1×1級的株數(shù)+2×2級的株數(shù)+3×3級的株數(shù)+4×4級的株數(shù))/處理的材料總株數(shù)。

    圖1 麗江新團黑谷-沈農(nóng)265重組自交系群體的物理圖譜(A)和基因型(B)

    物理圖譜位置基于日本晴參考基因組(MSU Rice Genome Annotation Project Release 7), 黃色表示沈農(nóng)265; 藍色表示麗江新團黑谷。

    Physical position is based on MSU Rice Genome Annotation Project Release 7 sequence. Yellow: SN265 genotype; Blue: LTH genotype.

    圖2 水稻芽期冷害評價等級標(biāo)準(zhǔn)

    1.3 芽期耐冷性的QTL鑒定和QTL效應(yīng)分析

    采用R/qtl的CIM方法進行QTL定位, 采用命令進行排列組合1000次的LOD閾值(α=0.05)確定, 當(dāng)實際求得的LOD值大于LOD閾值時, 就認為該區(qū)段存在1個QTL, 其置信區(qū)間為LOD峰值向下1個LOD值單位的區(qū)間[37]。

    為了進一步驗證本試驗定位結(jié)果的準(zhǔn)確性, 同時評價篩選到的QTL累加效應(yīng), 利用“選擇作圖”策略進行QTL驗證和累加效應(yīng)分析[38]。首先將控制芽期耐冷的其他QTL位點根據(jù)分子標(biāo)記數(shù)據(jù)進行非增效基因的固定, 以此來避免其他QTL對后續(xù)分析的影響, 然后在RILs群體中篩選目標(biāo)QTL區(qū)間的重組交換單株, 進而對QTL進行驗證和累加效應(yīng)分析。

    2 結(jié)果與分析

    2.1 親本的芽期耐冷性比較及群體的芽期耐冷性分布

    課題組前期的研究表明, 麗江新團黑谷與沈農(nóng)265在多個不同時期的耐冷相關(guān)性狀上存在顯著差異, 麗江新團黑谷具有更優(yōu)秀的綜合耐冷特性。本研究中, 2個親本在2°C低溫處理2 d并恢復(fù)10 d后, 沈農(nóng)265的平均芽期耐冷等級僅為0.813, 而麗江新團黑谷高達3.188 (圖3和圖4), 差異達到極顯著水平。重組自交系群體的芽期耐冷性等級分布范圍為0.313~3.688, 平均值為1.699, 標(biāo)準(zhǔn)差為0.767。整體上呈現(xiàn)近似正態(tài)分布, 并表現(xiàn)出較強的超親分離(圖4)。

    圖3 親本麗江新團黑谷和沈農(nóng)265芽期不經(jīng)冷處理(A)和冷處理后(B)的比較

    圖4 重組自交系群體的芽期耐冷性分布

    2.2 群體芽期耐冷性的QTL分布與效應(yīng)

    利用R/qtl軟件共檢測到5個控制水稻芽期耐冷性的QTL, 分布在水稻的1號、3號、9號和11號染色體上, 這些QTL的增效等位基因均來自于麗江新團黑谷(表1和圖5)。表型貢獻率最大的主效QTL是(LOD值為24.01), 位于11號染色體長臂端的21.24 Mb~22.03 Mb之間, 區(qū)間大小0.79 Mb。另外兩個貢獻率大于15%的主效QTL是(LOD值為5.36)和(LOD值為8.58), 分別位于水稻第9染色體上19.90 Mb~22.30 Mb之間和11號染色體上8.53 Mb~9.93 Mb之間, 區(qū)間大小分別為2.40 Mb和1.40 Mb。兩個微效QTL是(LOD值為3.26)和(LOD值為3.05), 分別位于水稻1號染色體上0.43 Mb~0.93 Mb之間和3號染色體上29.40 Mb~32.87 Mb之間, 區(qū)間大小分別為0.50 Mb和0.61 Mb。沒能從沈農(nóng)265中檢測到芽期耐冷性的增效等位QTL。

    表1 利用麗江新團黑谷-沈農(nóng)265的重組自交系群體檢測的芽期耐冷QTL信息

    a在水稻物理圖譜上的位置(單位: Mb)。

    aPositions in physical linkage map (unit: Mb). LTH: Lijiangxintuanheigu.

    圖5 控制芽期耐冷性的QTL在染色體上的位置分布

    : 低溫誘導(dǎo)的鋅指蛋白[46];: 參與冷脅迫的E3泛素連接酶[47];,: 低溫誘導(dǎo)的AP2/ EREBP轉(zhuǎn)錄因子基因[48];: 芽期耐冷QTL[27]。

    : cold inducible zinc finger protein[46];: a rice E3-Ubiquitin Ligase in the modulation of cold stress response[47];: cold inducible AP2/EREBP transcription factor gene[48];: a QTL for cold tolerance at bud stage[27].

    2.3 芽期耐冷性QTL的驗證和累加效應(yīng)分析

    我們對檢測到的芽期耐冷QTL進行了驗證并評價了這些QTL的互作和累加效應(yīng)。首先, 將控制芽期耐冷性的其他QTL利用分子標(biāo)記輔助進行非增效基因的固定, 然后在重組自交系群體中篩選目標(biāo)區(qū)間(待評價QTL區(qū)間)為增效等位基因的單株, 進而對這5個QTL進行驗證。隨后, 利用類似的策略, 通過標(biāo)記輔助選擇出不同QTL組合的單株, 進行QTL的互作分析和累加效應(yīng)評價。不導(dǎo)入任何芽期耐冷QTL株系的芽期低溫抗性等級為0.80 (圖6和圖7-A, X), 含有1個QTL株系的芽期低溫抗性等級從1.03到2.33 (圖6和圖7-C~F), 含有的株系具有最高的芽期低溫抗性等級(圖6和圖7-G)。含有2個QTL株系的芽期低溫抗性等級從1.67到2.82 (圖6和圖7-H~M), 含有株系的芽期低溫抗性等級普遍較高(2.44~2.82), 2個QTL之間的作用總體上表現(xiàn)為簡單的加性效應(yīng), 互作效應(yīng)較弱。由于受群體規(guī)模的影響, 并沒有覆蓋到全部的3個QTL聚合類型, 篩選到的6個類型含有3個QTL株系的芽期低溫抗性等級從1.33到2.98 (圖6和圖7-N~S), 與兩個QTL聚合的規(guī)律類似, 含有株系的芽期低溫抗性等級普遍較高(2.52~2.98), 3個QTL之間的作用總體上也表現(xiàn)為簡單的加性效應(yīng)。4個QTL聚合的株系檢測到3個類型, 均含有和, 它們的芽期低溫抗性等級從3.00到3.06 (圖6和圖7-T~V)。共檢測到同時聚合5個QTL的株系8個, 它們的芽期低溫抗性等級從2.94到3.50 (圖6和圖7-W)。通過上述分析, 基本確定了芽期耐低溫QTL在重組自交系群體中聚合時主要表現(xiàn)為加性效應(yīng)。

    圖6 芽期耐冷QTL的驗證和聚合效應(yīng)分析使用的重組單株基因型

    “-”: 缺失。“-”: missing.

    3 討論

    3.1 水稻芽期耐冷性評價方法的比較

    水稻芽期耐冷能力是決定水稻直播成苗率和整齊度的重要農(nóng)藝性狀, 尤其是在北方稻區(qū)和南方雙季稻早稻進行直播栽培的重要抗逆性狀。近年來, 隨著直播技術(shù)的改良, 越來越多的地區(qū)在水稻直播時選擇“芽種直播”或“盲谷播種”, 使得芽期耐冷特性成為決定品種是否適宜直播栽培的關(guān)鍵因素。因此, 準(zhǔn)確評價試驗材料的芽期耐冷特性就變得越來越重要。芽期耐冷鑒定方法中比較常用的是利用低溫培養(yǎng)箱處理萌芽后的種子, 并利用秧苗成活率作為評價指標(biāo), 這種方法簡單、直觀、結(jié)果可靠[22-33]。本研究對這種方法進行了進一步的改進, 主要包括以下幾個方面: 一是低溫處理的溫度和時間, 前人報道中, 多采用5°C 5 d[22-33], 本研究采用了2°C 2 d,既滿足了低溫處理的低溫條件, 又縮短了處理時間。二是低溫處理后的秧苗培養(yǎng)方式, 前人多使用全程水培[22-25], 本研究選擇土培, 這樣更接近水稻直播生產(chǎn)實際。三是芽期冷害等級的評價方式, 前人多是簡單地將處理緩苗后的秧苗成活率作為芽期耐冷性鑒定的指標(biāo)[22-25], 本研究為了更全面地反映每一個株系的綜合抗性, 采用了不同級別加權(quán)合計的方式評價某一株系整體的耐冷水平。通過上述改進, 不僅縮短了處理時間, 而且還提高了鑒定的精準(zhǔn)度。

    3.2 重測序技術(shù)和高密度bin標(biāo)記遺傳圖譜在水稻QTL分析中的應(yīng)用

    隨著重測序成本的大幅下降和測序數(shù)據(jù)利用技術(shù)的日益成熟, 利用重測序技術(shù)結(jié)合“滑動窗口”法構(gòu)建bin遺傳圖譜已成為水稻[39-41]、玉米[42-43]、高粱[44]等作物農(nóng)藝性狀遺傳研究的重要技術(shù)手段。與傳統(tǒng)的分子標(biāo)記相比, 重測序獲得的SNP標(biāo)記具有高基因組覆蓋度和高密度的優(yōu)勢, 它不僅顯著提高了QTL鑒定和定位的準(zhǔn)確性和精確性, 而且能夠為后續(xù)的基因功能研究提供相對較小的染色體區(qū)間參考。本研究利用包含2818個bin標(biāo)記的遺傳圖譜分別將控制水稻芽期耐冷性的QTL定位到了0.50 Mb~2.40 Mb的區(qū)間內(nèi), 最小的物理區(qū)間分別是500 kb ()和790 kb (), 相對較大的物理區(qū)間也僅分別是1.40 Mb ()、2.40 Mb ()和3.47 Mb (), 這遠遠小于傳統(tǒng)分子標(biāo)記高達5~10 cM的定位區(qū)間(約4 Mb~10 Mb)。500 kb左右的QTL定位區(qū)間雖然離直接完成基因克隆尚有一定的差距[45], 但是也為后續(xù)的精細定位提供了足夠小的物理區(qū)間。此外, 傳統(tǒng)的QTL分析作圖群體主要是秈-粳亞種間雜交得到的, 少數(shù)為秈-秈交群體, 很少使用兩個粳稻作為親本材料的遺傳群體, 主要是因為粳稻亞種內(nèi)的多態(tài)性遠遠小于秈稻亞種內(nèi), 與秈-粳亞種間的多態(tài)性相比, 更是差距極大[46]。重測序技術(shù)成為解決這一問題的有效手段, 使得利用粳-粳交群體也可以構(gòu)建高密度的遺傳圖譜, 為研究粳稻食味品質(zhì)、孕穗期障礙型冷害等要求雙親均為粳稻的實驗提供了極大的便利。

    圖7 控制芽期耐冷性的QTL效應(yīng)及在LTH-SN265重組自交系群體中的聚合效果分析

    A: 親本沈農(nóng)265 (SN265); B: 親本麗江新團黑谷(LTH); C~G:、、、和的芽期耐冷效應(yīng); H:和的芽期耐冷聚合效應(yīng); I:和的芽期耐冷聚合效應(yīng); J:和的芽期耐冷聚合效應(yīng); K:和的芽期耐冷聚合效應(yīng); L:和的芽期耐冷聚合效應(yīng); M:和的芽期耐冷聚合效應(yīng); N:、和的芽期耐冷聚合效應(yīng); O:、和的芽期耐冷聚合效應(yīng); P:、和的芽期耐冷聚合效應(yīng); Q:、和的芽期耐冷聚合效應(yīng); R:、和的芽期耐冷聚合效應(yīng); S:、和的芽期耐冷聚合效應(yīng); T:、、和的芽期耐冷聚合效應(yīng); U:、、和的芽期耐冷聚合效應(yīng); V:、、和的芽期耐冷聚合效應(yīng); W:、、、和的芽期耐冷聚合效應(yīng); X: 無QTL的株系表型。

    A: Shennong 265 (SN265); B: Lijiangxintuanheigu (LTH); C–G: Cold tolerance effect of,,,, andat bud stage; H: Cold tolerance cluster effect ofandat bud stage; I: Cold tolerance cluster effect ofandat bud stage; J: Cold tolerance cluster effect ofandat bud stage; K: Cold tolerance cluster effect ofandat bud stage; L: Cold tolerance cluster effect ofandat bud stage; M: Cold tolerance cluster effect ofandat bud stage; N: Cold tolerance cluster effect of,,andat bud stage; O: Cold tolerance cluster effect of,,andat bud stage; P: Cold tolerance cluster effect of,,andat bud stage; Q: Cold tolerance cluster effect of,,andat bud stage; R: Cold tolerance cluster effect of,,andat bud stage; S: Cold tolerance cluster effect of,,andat bud stage; T: Cold tolerance cluster effect of,,,andat bud stage; U: Cold tolerance cluster effect of,,andat bud stage; V: Cold tolerance cluster effect of,,,andat bud stage; W: Cold tolerance cluster effect of,,,,andat bud stage; X: none QTL line.

    3.3 耐冷QTL全部來自LTH的可能原因及利用RIL群體進行QTL聚合評估的可行性

    本研究中發(fā)現(xiàn)SN265和LTH構(gòu)建的重組自交系群體存在芽期耐冷性的超親分離, 但檢測到的耐冷QTL全部來自LTH, 其原因主要是QTL檢測時, 為了盡可能降低假陽性, 使用了相對較高的LOD閾值(3.0)。如果將LOD閾值降低到2.4, 能檢測到一個可能的QTL-, 其LOD值僅為2.45, 表型貢獻率很小, 僅為5.8%, 加性效應(yīng)值為?0.13。該QTL的耐冷增效等位基因來自沈農(nóng)265, 但由于表型貢獻率太小, 還需要構(gòu)建該QTL的高世代回交群體進行真?zhèn)悟炞C。此外, 對檢測的QTL進行效應(yīng)評價, 尤其是對QTL進行聚合效果分析是QTL用于水稻育種實踐的重要前提。但是, 一般而言, 由于重組自交系中導(dǎo)入片段相對較多, 排除遺傳背景的噪音影響難度很大, 很可能影響QTL累加效應(yīng)評估。所以, 重組自交系較少用于QTL的累加效應(yīng)評估。本研究采用了“選擇作圖”的方法進行QTL的累加效應(yīng)分析, 該方法的關(guān)鍵就是充分發(fā)揮已構(gòu)建的高密度物理圖譜和詳細的群體基因型數(shù)據(jù)將目標(biāo)QTL以外的其他QTL進行基因型固定, 以盡可能排除遺傳背景的影響。我們之前也成功利用這一方法分析了水稻根系的QTL聚合效果[38]。當(dāng)然, 以選擇作圖分析QTL聚合效果的能力肯定不如近等基因系或染色體片段代換系, 但在沒有更高級的遺傳群體時, 仍不失為一種有效手段。

    3.4 本研究鑒定的水稻芽期耐低溫QTL與其他研究結(jié)果的比較和育種應(yīng)用潛力

    隨著農(nóng)村勞動力的大規(guī)模轉(zhuǎn)移和農(nóng)業(yè)生產(chǎn)成本的不斷攀升, 水稻直播生產(chǎn)越來越受到黑龍江省地廣人稀的農(nóng)業(yè)主產(chǎn)區(qū)農(nóng)戶的重視。但是, 由于獨特的寒地生態(tài)環(huán)境, 黑龍江水稻直播不得不面對春季的低溫環(huán)境。水稻芽期耐冷性是在低溫下水稻幼芽細胞維持生活的能力, 從外觀上表現(xiàn)為幼芽誘發(fā)和維持綠苗的能力, 是水稻直播必具備的重要農(nóng)藝性狀[22]。本文從芽期耐冷性強的云南地方粳稻品種麗江新團黑谷中, 鑒定出5個控制芽期耐低溫能力的QTL, 同QTL數(shù)據(jù)庫(Q-TARO: http://qtaro.abr.affrc. go.jp/)搜集整理的耐冷相關(guān)QTL比對發(fā)現(xiàn),附近有一個已知的受冷誘導(dǎo)的鋅指蛋白基因, 該基因過量表達可以增強水稻耐冷性, 同時增加脯氨酸含量[47]。區(qū)間內(nèi)包含一個已知的耐冷相關(guān)基因, 它通過與互作來參與水稻的冷脅迫反應(yīng)[48]。區(qū)間內(nèi)包含2個成簇存在的耐冷相關(guān)基因, 分別是和, 二者均受冷誘導(dǎo), 參與水稻的冷脅迫反應(yīng)[49]。附近沒有發(fā)現(xiàn)已經(jīng)鑒定的冷害相關(guān)QTL和基因, 是一個新發(fā)現(xiàn)的位點。主效QTL-與Baruah等[27]鑒定的芽期耐冷QTL-位于相同的區(qū)間,來源于北海道地方品種A58。綜上可以看出, 控制水稻芽期耐冷性的QTL在不同遺傳背景條件下具有一定的重演性, 而且第11染色體上鑒定的主效位點具有較大的后續(xù)研究價值。我們也利用“選擇作圖”的方法評價了這些QTL的育種潛力和聚合效果, 發(fā)現(xiàn)芽期耐低溫QTL主要表現(xiàn)為加性效應(yīng), 非常適于水稻的分子輔助育種。而且主效QTL-不僅自身具有很大的遺傳效應(yīng), 而且在與其他QTL聚合時, 常常表現(xiàn)出更加突出的改良效果, 極具育種應(yīng)用價值。

    4 結(jié)論

    構(gòu)建了麗江新團黑谷和沈農(nóng)265重組自交系群體包含2818個bin標(biāo)記的遺傳圖譜。利用此群體和bin標(biāo)記圖譜, 共檢測到5個芽期耐冷QTL, 分布在水稻的1號、3號、9號和11號染色體上, 表型貢獻率分別為8.0%、8.2%、16.6%、24.4%和53.5%, 增效等位基因均來自耐冷親本麗江新團黑谷。主效QTL-位于11號染色體長臂端的17.97 Mb~23.03 Mb之間。利用“選擇作圖”策略進行了QTL驗證和累加效應(yīng)分析, 明確了可以通過QTL的累加聚合實現(xiàn)芽期耐冷能力的遺傳改良, 聚合的增效QTL越多, 耐冷能力等級越高。

    [1] 劉次桃, 王威, 毛畢剛, 儲成才. 水稻耐低溫逆境研究: 分子生理機制及育種展望. 遺傳, 2018, 40: 171–185. Liu C T, Wang W, Mao B G, Chu C C. Cold stress tolerance in rice: physiological changes, molecular mechanism, and future prospects., 2018, 40: 171–185 (in Chinese with English abstract).

    [2] 矯江, 許顯濱, 孟英. 黑龍江省水稻低溫冷害及對策研究. 中國農(nóng)業(yè)氣象, 2004, 25(2): 27–29. Jiao J, Xu X B, Meng Y. Analysis of rice chilling injury and countermeasures in Heilongjiang province., 2004, 25(2): 27–29 (in Chinese with English abstract).

    [3] 聶元元, 蔡耀輝, 顏滿蓮, 李永輝, 毛凌華, 顏龍安, 楊曉莉. 水稻低溫冷害分析研究進展. 江西農(nóng)業(yè)學(xué)報, 2011, 23(3): 63–66. Nie Y Y, Cai Y H, Yan M L, Li Y H, Mao L H, Yan L A, Yang X L. Research advances in chilling injury to rice., 2011, 23(3): 63–66 (in Chinese with English abstract).

    [4] Cruz R P, Sperotto R A, Cargnelutti D, Adamski J M, de FreitasTerra T, Fett J P. Avoiding damage and achieving cold tolerance in rice plants., 2013, 2: 96–119.

    [5] Borjas A H, Leon T B D, Subudhi P K. Genetic analysis of germinating ability and seedling vigor under cold stress in US weedy rice., 2016, 208: 1–14.

    [6] Satoh T, Tezuka K, Kawamoto T, Matsumoto S, Satoh-Nagasawa N, Ueda K, Sakurai K, Watanabe A, Takahashi H, Akagi H. Identification of QTLs controlling low-temperature germination of the east European rice variety Maratteli., 2016, 207: 245–254.

    [7] 吳立, 霍治國, 姜燕, 張蕾, 于彩霞. 氣候變暖背景下南方早稻春季低溫災(zāi)害的發(fā)生趨勢與風(fēng)險. 生態(tài)學(xué)報, 2016, 36: 1263–1271. Wu L, Huo Z G, Jiang Y, Zhang L, Yu C X. Trends and risk of spring low-temperature damage to early rice in southern China against the background of global warming., 2016, 36: 1263–1271 (in Chinese with English abstract).

    [8] 郭麗穎, 耿艷秋, 金峰, 宋微, 邵璽文. 寒地水稻低溫冷害防御栽培技術(shù)研究進展. 作物雜志, 2017, (4): 7–14. Guo L Y, Geng Y Q, Jin F, Song W, Shao X W. Research advances about low temperature, cold damage defense cultivation techniques of rice in cold region of china., 2017, (4): 7–14 (in Chinese with English abstract).

    [9] 戴陸園, 葉昌榮, 余騰瓊, 徐福榮.水稻耐冷性研究: I. 稻冷害類型及耐冷性鑒定評價方法概述.西南農(nóng)業(yè)學(xué)報, 2002, 15(1): 41–45. Dai L Y, Ye C R, Yu T Q, Xu F R. Studies on cold tolerance of rice,L. : I. Description on types of cold injury and classifications of evaluation methods on cold tolerance in rice., 2002, 15(1): 41–45 (in Chinese with English abstract).

    [10] 韓龍植, 曹桂蘭, 芮鐘斗, 安永平, 喬永利, 黃興九, 高熙宗. 水稻芽期耐冷性與其他耐冷性狀的相關(guān)關(guān)系. 作物學(xué)報, 2004, 30: 990–995. Han L Z, Cao G L, Yea J D, An Y P, Qiao Y L, Huang X J, Koh H J. Relationship between cold tolerance at the bud bursting period and other traits related to cold tolerance in rice., 30: 990–995 (in Chinese with English abstract).

    [11] 陳品, 陸建飛. 長江中下游地區(qū)直播稻的生理生態(tài)特性及其栽培技術(shù)的研究進展. 核農(nóng)學(xué)報, 2013, 27: 487–494. Chen P, Lu J F. Research advances on the physiological and ecological characteristics and cultivation techniques of direct seeding rice in the Middle-Lower Yangtze area., 2013, 27: 487–494 (in Chinese with English abstract).

    [12] 張喜娟, 來永才, 孟英, 唐傲, 董文軍, 冷春旭, 王立志, 姜樹坤, 姜輝, 丁國華. 水直播對寒地粳稻產(chǎn)量和品質(zhì)性狀的影響. 中國稻米, 2016, 22(2): 43–46. Zhang X J, Lai Y C, Meng Y, Tang A, Dong W J, Leng C X, Wang L Z, Jiang S K, Jiang H, Ding G H. Effects of water direct seeding on yield and quality ofrice in cold area., 2016, 22(2): 43–46 (in Chinese with English abstract).

    [13] 余會康, 郭建平. 氣候變化下東北水稻冷害時空分布變化. 中國生態(tài)農(nóng)業(yè)學(xué)報, 2014, 22: 594?601. Yu H K, Guo J P. Variation in spatial and temporal distribution of chilling injury of rice under climate change in Northeast China., 2014, 22: 594?601 (in Chinese with English abstract).

    [14] Fujino K, Sekiguchi H, Matsuda Y, Sugimoto K, Ono K, Yano M. Molecular identification of a major quantitative trait locus,, controlling low-temperature germinability in rice., 2008, 105: 12623–12628.

    [15] Liu C T, Ou S, Mao B, Tang J, Wang W, Wang H, Cao S, Schl?ppi M, Zhao B, Xiao G. Early selection of bZIP73 facilitated adaptation ofrice to cold climates., 2018, 9: 3302.

    [16] Saito K, Hayano-Saito Y, Kuroki M, Sato Y. Map-based cloning of the rice cold tolerance gene., 2010, 179: 97–102.

    [17] Zhang Z Y, Li J J, Pan Y H, Li J L, Zhou L, Shi H L, Zeng Y W, Guo H F, Yang S M, Zheng W W, Yu J P, Sun X M, Li G L, Ding Y L, Ma L, Shen S Q, Dai L Y, Zhang H L, Yang S H, Guo Y, Li Z C. Natural variation in CTB4a enhances rice adaptation to cold habitats., 2017, 8: 14788.

    [18] Ma Y, Dai X Y, Xu Y Y, Luo W, Zheng X M, Zeng D L, Pan Y J, Lin X L, Liu H H, Zhang D J, Xiao J, Guo X Y, Xu S J, Niu Y D, Jin J B, Zhang H, Xu X, Li L L, Wang W, Qian Q, Ge S, Chong K. COLD1 confers chilling tolerance in rice., 2015, 160: 1209–1221.

    [19] Mao D H, Xin Y Y, Tan Y J, Hu X J, Bai J J, Liu Z Y, Yu Y L, Li L Y, Peng C, Fan T, Zhu Y X, Guo Y L, Wang S H, Li D P, Xing Y Z, Yuan L P, Chen C Y. Natural variation in thegene confers chilling tolerance in rice and allowed adaptation to a temperate climate., 2019, 116: 3494–3501.

    [20] 嚴長杰, 李欣, 程祝寬, 于恒秀, 顧銘洪, 朱立煌. 利用分子標(biāo)記定位水稻芽期耐冷性基因. 中國水稻科學(xué), 1999, 13: 134–138. Yan C J, Li X, Cheng Z K, Yu H X, Gu M H, Zhu L H. Identification of QTL for cold tolerance at early seeding stage in rice () via RFLP markers., 1999, 13: 134–138 (in Chinese with English abstract).

    [21] Zhang Z H, Li S, Li W, Chen W, Zhu Y G. A major QTL conferring cold tolerance at the early seedling stage using recombinant inbred lines of rice (L.)., 2005, 168: 527–534.

    [22] 喬永利, 韓龍植, 安永平, 張媛媛, 曹桂蘭, 高熙宗. 水稻芽期耐冷性QTL的分子定位. 中國農(nóng)業(yè)科學(xué), 2005, 38: 217–221. Qiao Y L, Han L Z, An Y P, Zhang Y Y, Cao G L, Koh H J. Molecular mapping of QTLs for cold tolerance at the bud bursting period in rice., 2005, 38: 217–221 (in Chinese with English abstract).

    [23] 陳瑋, 李煒. 水稻RIL群體芽期耐冷性基因的分子標(biāo)記定位. 武漢植物學(xué)研究, 2005, 23: 116–120. Chen W, Li W. Mapping of QTL conferring cold tolerance at early seeding stage of rice by molecular markers., 2005, 23: 116–120 (in Chinese with English abstract).

    [24] 張露霞, 王松鳳, 江鈴, 萬建民. 利用重組自交系群體檢測水稻芽期耐冷性QTL. 南京農(nóng)業(yè)大學(xué)學(xué)報, 2007, 30(4): 1–5. Zhang L X, Wang S F, Jiang L, Wan J M. QTL analysis of cold tolerance at the bud bursting period in rice (L.) by using recombinant inbred lines., 2007, 30(4): 1–5 (in Chinese with English abstract).

    [25] 鞏迎軍, 阮雯君, 荀星, 董彥君, 林冬枝, 葉勝海, 張小明. 水稻芽性狀耐冷性的QTL分析. 分子植物育種, 2009, 7: 273–278. Gong Y J, Ruan W J, Xun X, Dong Y J, Lin D Z, Ye S H, Zhang X M. QTL analysis of cold tolerance for two bud traits in rice., 2009, 7: 273–278 (in Chinese with English abstract).

    [26] 林靜, 朱文銀, 張亞東, 朱鎮(zhèn), 趙凌, 陳濤, 趙慶勇, 周麗慧, 方先文, 王艷平. 利用染色體片段置換系定位水稻芽期耐冷性QTL. 中國水稻科學(xué), 2010, 24: 233–236. Lin J, Zhu W Y, Zhang Y D, Zhu Z, Zhao L, Chen T, Zhao Q Y, Zhou L H, Fang X W, Wang Y P. Detection of quantitative trait loci for cold tolerance at the bursting stage by using chromosome segment substitution lines in rice ()., 2010, 24: 233–236 (in Chinese with English abstract).

    [27] Baruah A R, Ishigo-Oka N, Adachi M, Oguma Y, Tokizono Y, Onishi K, Sano Y. Cold tolerance at the early growth stage in wild and cultivated rice., 2009, 165: 459–470.

    [28] Ji Z J, Zeng Y X, Zeng D L, Ma L Y, Li X M, Liu B X, Yang C D. Identification of QTLs for rice cold tolerance at plumule and 3-leaf-seedling stages by using QTL Network software., 2010, 17: 282–287.

    [29] 周勇, 朱孝波, 袁華, 鄭英, 欽鵬, 魏應(yīng)海, 王玉平, 黃世君, 李仕貴. 水稻單片段代換系芽期和苗期耐冷性分析及耐冷性QTL鑒定. 中國水稻科學(xué), 2013, 27: 381–388. Zhou Y, Zhu X B, Yuan H, Zheng Y, Qin P, Wei Y H, Wang Y P, Huang S J, Li S G. Characterization of cold tolerance and identification of cold tolerance QTLs for rice single segment substitution lines at plumule and seedling stages., 2013, 27: 381–388 (in Chinese with English abstract).

    [30] 楊洛淼, 王敬國, 劉化龍, 孫健, 鄭洪亮, 鄒德堂. 寒地粳稻發(fā)芽期和芽期的耐冷性QTL定位. 作物雜志, 2014, (6): 44–51. Yang L S, Wang J G, Liu H L, Sun J, Zheng H L, Zou D T. QTL mapping of cold tolerance for the germination period and the bud bursting period ofin cold area., 2014, (6): 44–51 (in Chinese with English abstract).

    [31] 朱金燕, 楊梅, 嵇朝球, 王軍, 楊杰, 范方軍, 李文奇, 王芳權(quán), 梁國華, 周勇, 仲維功. 利用染色體單片段代換系定位水稻芽期耐冷QTL. 植物學(xué)報, 2015, 50: 338–345. Zhu J Y, Yang M, Ji C Q, Wang J, Yang J, Fan F J, Li W Q, Wang F Q, Liang G H, Zhou Y, Zhong W G. Identification of cold tolerance at the plumule stage quantitative trait loci with single segment substituted lines in rice., 2015, 50: 338–345 (in Chinese with English abstract).

    [32] Yang T F, Zhang S H, Zhao J L. Liu Q, Huang Z H, Mao X X, Dong J F, Wang X F, Zhang G Q, Liu B. Identification and pyramiding of QTLs for cold tolerance at the bud bursting and the seedling stages by use of single segment substitution lines in rice (L.)., 2016, 36: 96.

    [33] Zhang M C, Ye J, Xu Q, Feng Y, Yuan X P, Yu H Y, Wang Y P, Wei X H, Yang Y L. Genome-wide association study of cold tolerance of Chineserice varieties at the bud burst stage., 2018, 37: 529–539.

    [34] Huang X, Feng Q, Qian Q, Zhao Q, Wang L, Wang A, Guan J, Fan D, Weng Q, Huang T. High-throughput genotyping by whole-genome resequencing., 2009, 19: 1068–1076.

    [35] Jiang S, Yang C, Xu Q, Wang L, Yang X, Song X, Wang J, Zhang X, Li B, Li H, Li Z, Li W. Genetic dissection germinability under low temperature by building a resequencing linkage map inrice., 2020, 21: 1284.

    [36] 韓龍植, 魏興華. 水稻種質(zhì)資源描述規(guī)范和數(shù)據(jù)標(biāo)準(zhǔn). 北京: 中國農(nóng)業(yè)出版社, 2006. p 107. Han L Z, Wei X H. Descriptors and Data Standard for Rice (L.). Beijing: China Agriculture Press, 2006. p 107.

    [37] Arends D, Prins P, Jansen R, Broman K. R/QTL: high-throughput multiple QTL mapping., 2010, 26: 2990–2992.

    [38] 姜樹坤, 張鳳鳴, 白良明, 孫世臣, 王彤彤, 丁國華, 姜輝, 張喜娟. 水稻移栽后新生根系相關(guān)性狀的QTL分析. 中國水稻科學(xué), 2014, 28: 598–604. Jiang S K, Zhang F M, Bai L M, Sun S C, Wang T T, Ding G H, Jiang H, Zhang X J. QTL analysis on new root traits after rice transplanting., 2014, 28: 598–604 (in Chinese with English abstract).

    [39] 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., 2011, 6: e17595.

    [40] Gao Z Y, Zhao S C, He W M, Guo L B, Peng Y L, Wang J J, Guo X S, Zhang X M, Rao Y C, Zhang C, Dong G J, Zheng F Y, Lu C X, Hu J, Zhou Q, Liu H J, Wu H Y, Xu J, Ni P X, Zeng D L, Liu D H, Tian P, Gong L H, Ye C, Zhang G H, Wang J, Tian F K, Xue D W, Liao Y, Zhu L, Chen M S, Li J Y, Cheng S H, Zhang G Y, Wang J, Qian Q. Dissecting yield-associated loci in super hybrid rice by resequencing recombinant inbred lines and improving parental genome sequences., 2013, 110: 14492–14497.

    [41] Jiang N, Shi S, Shi H, Khanzada H, Wassan G M, Zhu C, Peng X, Yu Q, Chen X, He X, Fu J, Hu L, Xu J, Ou-Yang L, Sun X, Zhou D, He H, Bian J. Mapping QTL for seed germinability under low temperature using a new high-density genetic map of rice., 2017, 8: 1223.

    [42] Liu H, Niu Y, Gonzalez-Portilla P J, Zhou H, Wang L, Zuo T, Qin C, Tai S, Jansen C, Shen Y, Lin H, Lee M, Ware D, Zhang Z, Lübberstedt T, Pan G. An ultra-high-density map as a community resource for discerning the genetic basis of quantitative traits in maize., 2015, 16: 1078.

    [43] Song W, Wang B, Hauck A L, Dong X, Li J, Lai J. Genetic dissection of maize seedling root system architecture traits using an ultra-high density bin-map and a recombinant inbred line population., 2016, 58: 266–279.

    [44] Zou G, Zhai G, Feng Q, Yan S, Wang A, Zhao Q, Shao J, Zhang Z, Zou J, Han B, Tao Y. Identification of QTLs for eight agronomically important traits using an ultra-high-density map based on SNPs generated from high-throughput sequencing in sorghum under contrasting photoperiods., 2012, 63: 5451–5462.

    [45] 姜樹坤, 徐正進, 陳溫福. 水稻QTL圖位克隆的特征分析. 遺傳, 2008, 30: 1121–1126. Jiang S K, Xu Z J, Chen W F. Analysis of features of 15 successful positional cloning of QTL in rice.(Beijing), 2008, 30: 1121–1126 (in Chinese with English abstract).

    [46] Jiang S K, Huang C, Zhang X J, Wang J Y, Xu Z J, Chen W F. Development of a high informative microsatellite markers (SSRs) framework for genotyping of rice (L.)., 2010, 9: 1697–1704.

    [47] Liu K M, Wang L, Xu Y Y, Chen N, Ma Q B, Li F, Chong K. Overexpression of, a putative cold inducible zinc finger protein, increased tolerance to chilling, salt and drought, and enhanced proline level in rice., 2007, 226: 1007–1016.

    [48] Louren?o T, Sapeta H, Figueiredo D D, Rodrigues M, Cordeiro A, Abreu I A, Saibo N J M, Oliveira M M. Isolation and characterization of rice (L.) E3-ubiquitin ligasegene in the modulation of cold stress response., 2013, 83: 351–363.

    [49] Mao D H, Chen C Y. Colinearity and similar expression pattern of ricereveal their functional conservation in the cold-responsive pathway., 2012, 7: e47275.

    Detection of QTLs controlling cold tolerance at bud bursting stage by using a high-density SNP linkage map inrice

    JIANG Shu-Kun1,*, WANG Li-Zhi1, YANG Xian-Li1, LI Bo2, MU Wei-Jie3, DONG Shi-Chen3, CHE Wei-Cai3, LI Zhong-Jie1, CHI Li-Yong1, LI Ming-Xian1, ZHANG Xi-Juan1, JIANG Hui2, LI Rui1, ZHAO Qian1, and LI Wen-Hua2,*

    1Institute of Crop Cultivation and Tillage, Heilongjiang Academy of Agricultural Sciences / Heilongjiang Provincial Key Laboratory of Crop Physiology and Ecology in Cold Region / Heilongjiang Provincial Engineering Technology Research Center of Crop Cold Damage, Harbin 150086, Heilongjiang, China;2Heilongjiang Academy of Agricultural Sciences, Harbin 150086, Heilongjiang, China;3Life Science and Technology College, Harbin Normal University, Harbin 150025, Heilongjiang, China

    Direct seeded rice (DSR) has received much attention because of its time- and labour-saving and low-input demand. However, the long-term cultivation method of seedling-transplantation has led to loss of some low-temperature-tolerant genes expressed at the bud stage. It has made many currently popular rice varieties unsuitable for direct seeding production. Therefore, it is important to identify cold-tolerance genes at the bud stage and to provide genes for subsequent molecular marker assistant breeding. In this study, we used a recombinant inbred line population constructed by cross of Lijiangxintuanheigu (LTH) and Shennong 265 (SN265) and its linkage map containing 2818 markers to detect cold tolerance QTLs at the bud stage. A total of five QTLs were detected on rice chromosomes 1, 3, 9, and 11. All the cold tolerance alleles were from the cold-tolerant parent LTH. The LOD values of these QTLs ranged from 3.05 to 24.01, and the phenotypic variations ranged from 8.0% to 53.5%. Among them, the major QTLwas located in a 790 kb interval of 21.24 Mb to 22.03 Mb on the long arm of chromosome 11. Subsequently, the “selective mapping” strategy was used for QTL verification and pyramiding effect analysis. Genetic improvement of cold tolerance at the bud stage would be achieved through pyramiding more QTLs. These results not only promote people’s understanding of the genetic basis for cold tolerance at the bud stage in rice but also provide theoretical basis and technical guidance for genetic improvement of DSR varieties.

    rice; cold tolerance at bud bursting stage; re-sequencing; genetic map; QTLs

    10.3724/SP.J.1006.2020.92066

    本研究由國家重點研發(fā)計劃項目(2018YFD0300105-5-2), 國家自然科學(xué)基金項目(31661143012), 國家重點研發(fā)項目黑龍江省級資助(768001), 黑龍江省自然科學(xué)基金優(yōu)秀青年基金項目(YQ2019C020)和黑龍江省農(nóng)業(yè)科學(xué)院長期性基礎(chǔ)性項目(2018CQJC002, 2019CQJC002)資助。

    This study was supported by the National Key Research and Development Program of China(2018YFD0300105-5-2), the National Natural Science Foundation of China(31661143012), the Provincial Funding for the National Key Research and Development Program in Heilongjiang Province (768001), the Natural Science Foundation of Heilongjiang Province of China (YQ2019C020), and the Heilongjiang Province Agricultural Science and Technology Innovation Project (2018CQJC002, 2019CQJC002).

    姜樹坤, E-mail:sk_jiang@126.com, Tel: 023-68251264, 李文華, E-mail: nkylwh@163.com

    E-mail: sk_jiang@126.com

    2019-12-08;

    2020-03-24;

    2020-04-10.

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

    猜你喜歡
    冷性自交系染色體
    優(yōu)質(zhì)水稻新品種芽期耐冷性鑒定
    北方水稻(2020年2期)2020-04-21 09:00:28
    多一條X染色體,壽命會更長
    玉米自交系京92遺傳改良研究
    為什么男性要有一條X染色體?
    不同來源玉米自交系穗粒性狀的多樣性分析與改良
    轉(zhuǎn)檉柳冷適應(yīng)蛋白基因水稻的獲得及抗冷性分析
    能忍的人壽命長
    SRAP結(jié)合SSR標(biāo)記分析油菜自交系的遺傳多樣性
    干旱脅迫對4份玉米自交系生理與光合特性的影響
    黃瓜果實耐冷性與CsHSFs基因表達關(guān)系的研究
    天美传媒精品一区二区| 伦精品一区二区三区| 看十八女毛片水多多多| 2021天堂中文幕一二区在线观| 国产精品一区二区三区四区久久| 日韩亚洲欧美综合| 亚洲在线自拍视频| 亚洲第一区二区三区不卡| 婷婷六月久久综合丁香| 久久韩国三级中文字幕| 内地一区二区视频在线| 国产精品1区2区在线观看.| www日本黄色视频网| 欧美一区二区亚洲| 夜夜夜夜夜久久久久| 亚洲国产色片| 99热这里只有是精品50| 欧美区成人在线视频| 国内揄拍国产精品人妻在线| 亚洲婷婷狠狠爱综合网| 亚洲国产精品国产精品| 全区人妻精品视频| 国产精品一区二区三区四区久久| 亚洲欧美清纯卡通| 亚洲在线自拍视频| 熟女电影av网| 国产精品一区二区三区四区久久| 国内精品宾馆在线| 美女被艹到高潮喷水动态| 国内精品美女久久久久久| 91av网一区二区| 偷拍熟女少妇极品色| 51国产日韩欧美| 日韩大尺度精品在线看网址| 高清午夜精品一区二区三区 | 日本三级黄在线观看| 少妇熟女aⅴ在线视频| 国产精品福利在线免费观看| 亚洲七黄色美女视频| 女人被狂操c到高潮| 中国美白少妇内射xxxbb| 亚洲精品乱码久久久久久按摩| 久久精品影院6| 亚洲精品国产av成人精品| 在线播放国产精品三级| 亚洲第一区二区三区不卡| 99热只有精品国产| 成人av在线播放网站| 天堂网av新在线| 在线天堂最新版资源| 免费av毛片视频| 一区二区三区免费毛片| 国产在视频线在精品| 日本成人三级电影网站| 五月玫瑰六月丁香| 少妇熟女欧美另类| 国产精品不卡视频一区二区| 欧美日本亚洲视频在线播放| 亚洲aⅴ乱码一区二区在线播放| 亚洲欧美清纯卡通| 色尼玛亚洲综合影院| 精品久久久久久久人妻蜜臀av| 熟妇人妻久久中文字幕3abv| 国产 一区精品| 亚洲av一区综合| 欧美3d第一页| 亚洲最大成人中文| 日韩中字成人| 一级毛片我不卡| 天堂影院成人在线观看| 一区二区三区高清视频在线| 哪里可以看免费的av片| 夜夜夜夜夜久久久久| 伦精品一区二区三区| 国产一区二区三区av在线 | 变态另类成人亚洲欧美熟女| 国产男人的电影天堂91| 乱系列少妇在线播放| 麻豆一二三区av精品| 中文字幕av成人在线电影| 国产一区二区亚洲精品在线观看| 成人二区视频| 看十八女毛片水多多多| 久久人人爽人人爽人人片va| 国产老妇伦熟女老妇高清| 又黄又爽又刺激的免费视频.| 日本黄色视频三级网站网址| 精品久久久久久久久av| 桃色一区二区三区在线观看| 成年版毛片免费区| 五月玫瑰六月丁香| 日本与韩国留学比较| 男女边吃奶边做爰视频| 中出人妻视频一区二区| 国产亚洲精品久久久久久毛片| 国产精品,欧美在线| 亚洲最大成人手机在线| 亚洲精华国产精华液的使用体验 | 高清毛片免费看| 三级男女做爰猛烈吃奶摸视频| 高清日韩中文字幕在线| 亚洲最大成人手机在线| 国产成人一区二区在线| 久久婷婷人人爽人人干人人爱| 一边摸一边抽搐一进一小说| 国产精品麻豆人妻色哟哟久久 | 国语自产精品视频在线第100页| 深夜a级毛片| 91在线精品国自产拍蜜月| 蜜桃久久精品国产亚洲av| 人人妻人人澡人人爽人人夜夜 | 九九在线视频观看精品| 一级毛片电影观看 | 别揉我奶头 嗯啊视频| 国产老妇女一区| 黑人高潮一二区| 99在线视频只有这里精品首页| 国产av麻豆久久久久久久| 男人舔奶头视频| 我要搜黄色片| 老女人水多毛片| 人妻夜夜爽99麻豆av| 国产一区二区激情短视频| 久久久久性生活片| www.av在线官网国产| 亚洲欧美成人综合另类久久久 | 国产色爽女视频免费观看| 悠悠久久av| 青春草亚洲视频在线观看| 女的被弄到高潮叫床怎么办| 日本-黄色视频高清免费观看| 国产成人精品婷婷| 久久精品国产99精品国产亚洲性色| 人妻制服诱惑在线中文字幕| 久久人人爽人人爽人人片va| 亚洲欧美精品综合久久99| av天堂中文字幕网| 久久久久久久午夜电影| 久久久久性生活片| 在线观看免费视频日本深夜| 国产一区二区激情短视频| 干丝袜人妻中文字幕| 白带黄色成豆腐渣| 男人舔奶头视频| 69人妻影院| 国产在视频线在精品| 欧美精品一区二区大全| 卡戴珊不雅视频在线播放| 午夜久久久久精精品| 成人毛片a级毛片在线播放| 激情 狠狠 欧美| 精品日产1卡2卡| 欧美日韩综合久久久久久| 国产亚洲91精品色在线| 成熟少妇高潮喷水视频| 91久久精品电影网| 日韩av在线大香蕉| 又黄又爽又刺激的免费视频.| 亚洲激情五月婷婷啪啪| 亚洲国产精品成人综合色| 真实男女啪啪啪动态图| 91aial.com中文字幕在线观看| 亚洲aⅴ乱码一区二区在线播放| 观看免费一级毛片| 亚洲激情五月婷婷啪啪| 日韩高清综合在线| 国产精品不卡视频一区二区| 爱豆传媒免费全集在线观看| 99热这里只有是精品在线观看| 老熟妇乱子伦视频在线观看| 免费不卡的大黄色大毛片视频在线观看 | 美女大奶头视频| 在线观看午夜福利视频| 国产亚洲精品av在线| 成人国产麻豆网| 欧美色视频一区免费| 国产成人一区二区在线| 国产男人的电影天堂91| 精品久久久久久成人av| 亚洲av中文av极速乱| 卡戴珊不雅视频在线播放| 日韩大尺度精品在线看网址| 最新中文字幕久久久久| 最近的中文字幕免费完整| 成人高潮视频无遮挡免费网站| 国产精品,欧美在线| 99热全是精品| 精品熟女少妇av免费看| av视频在线观看入口| 久久久久久久亚洲中文字幕| eeuss影院久久| 好男人在线观看高清免费视频| 色5月婷婷丁香| 99热网站在线观看| 51国产日韩欧美| 老熟妇乱子伦视频在线观看| 国产高清视频在线观看网站| 国产片特级美女逼逼视频| 国产欧美日韩精品一区二区| 美女大奶头视频| 中国国产av一级| 精品久久久久久成人av| 最好的美女福利视频网| 一本久久精品| 国产精品乱码一区二三区的特点| 国产精品日韩av在线免费观看| 蜜桃亚洲精品一区二区三区| 一进一出抽搐gif免费好疼| 天美传媒精品一区二区| 赤兔流量卡办理| 亚洲欧美成人综合另类久久久 | 午夜视频国产福利| 美女国产视频在线观看| 青春草国产在线视频 | 日日干狠狠操夜夜爽| 亚洲欧美日韩高清在线视频| 精品国产三级普通话版| 色视频www国产| 插逼视频在线观看| 69av精品久久久久久| 联通29元200g的流量卡| 国产蜜桃级精品一区二区三区| 美女 人体艺术 gogo| h日本视频在线播放| 国国产精品蜜臀av免费| 观看免费一级毛片| 久久这里只有精品中国| 韩国av在线不卡| 精品不卡国产一区二区三区| 女人十人毛片免费观看3o分钟| 成人一区二区视频在线观看| 国内揄拍国产精品人妻在线| 国产单亲对白刺激| 大型黄色视频在线免费观看| 国产亚洲精品久久久com| 亚洲性久久影院| 美女黄网站色视频| 99久久人妻综合| 亚洲精品乱码久久久久久按摩| 中国美女看黄片| 久久亚洲国产成人精品v| 日韩欧美三级三区| 亚洲精品乱码久久久v下载方式| 91午夜精品亚洲一区二区三区| 成人无遮挡网站| 亚洲性久久影院| 波多野结衣高清作品| 精品日产1卡2卡| 男女啪啪激烈高潮av片| 日本欧美国产在线视频| 免费av不卡在线播放| 中国国产av一级| 国产日韩欧美在线精品| 久久精品国产亚洲av香蕉五月| 麻豆精品久久久久久蜜桃| 亚洲中文字幕日韩| 一区二区三区免费毛片| 婷婷亚洲欧美| 69人妻影院| 亚洲人成网站在线观看播放| 九九久久精品国产亚洲av麻豆| 夜夜夜夜夜久久久久| 色哟哟哟哟哟哟| 69av精品久久久久久| 亚洲成av人片在线播放无| 99热只有精品国产| 亚洲天堂国产精品一区在线| 成人国产麻豆网| 国产亚洲欧美98| 国产女主播在线喷水免费视频网站 | 91aial.com中文字幕在线观看| 精品一区二区三区视频在线| 亚洲成人av在线免费| 欧美性感艳星| 日韩欧美在线乱码| 亚洲国产欧洲综合997久久,| av在线老鸭窝| av免费在线看不卡| 成人二区视频| 国产精品福利在线免费观看| 亚洲人与动物交配视频| 久久久久久九九精品二区国产| 好男人视频免费观看在线| 国产伦理片在线播放av一区 | 99久久成人亚洲精品观看| 国产 一区精品| 中文字幕精品亚洲无线码一区| 国产极品精品免费视频能看的| 亚洲aⅴ乱码一区二区在线播放| 18+在线观看网站| 亚洲熟妇中文字幕五十中出| 成人美女网站在线观看视频| 热99在线观看视频| 亚洲一区高清亚洲精品| 日日啪夜夜撸| 非洲黑人性xxxx精品又粗又长| 国产精品一区二区性色av| 六月丁香七月| 成年版毛片免费区| 色播亚洲综合网| 久久亚洲国产成人精品v| 波多野结衣高清作品| 在线播放无遮挡| 国产精品,欧美在线| 国产精品蜜桃在线观看 | 99国产精品一区二区蜜桃av| 中文字幕av在线有码专区| 我的女老师完整版在线观看| 久久久久久久久久久丰满| 只有这里有精品99| 国产黄色视频一区二区在线观看 | 美女大奶头视频| 简卡轻食公司| 身体一侧抽搐| 一本久久中文字幕| 国语自产精品视频在线第100页| 嫩草影院精品99| 啦啦啦韩国在线观看视频| 少妇被粗大猛烈的视频| 日本黄色片子视频| 久久99精品国语久久久| 中文字幕精品亚洲无线码一区| 精品熟女少妇av免费看| 最近手机中文字幕大全| 你懂的网址亚洲精品在线观看 | 日韩 亚洲 欧美在线| 免费av不卡在线播放| 我要搜黄色片| 啦啦啦韩国在线观看视频| 国产蜜桃级精品一区二区三区| 99久国产av精品| 亚洲真实伦在线观看| 在线免费十八禁| 中文亚洲av片在线观看爽| 日日摸夜夜添夜夜爱| 精品少妇黑人巨大在线播放 | 亚洲高清免费不卡视频| www.av在线官网国产| 在线观看美女被高潮喷水网站| 国产亚洲av片在线观看秒播厂 | 婷婷精品国产亚洲av| а√天堂www在线а√下载| 尾随美女入室| or卡值多少钱| 亚洲av中文字字幕乱码综合| 精品久久久久久久久av| 99热精品在线国产| 99久久精品热视频| 午夜福利成人在线免费观看| 一区福利在线观看| 校园春色视频在线观看| 一区福利在线观看| 99热这里只有精品一区| 国产蜜桃级精品一区二区三区| 欧美高清性xxxxhd video| 五月伊人婷婷丁香| 亚洲不卡免费看| 色噜噜av男人的天堂激情| 人妻久久中文字幕网| 哪个播放器可以免费观看大片| 亚洲成a人片在线一区二区| 日韩一本色道免费dvd| 国产精品美女特级片免费视频播放器| 日本免费一区二区三区高清不卡| 能在线免费观看的黄片| ponron亚洲| 国产真实乱freesex| 1000部很黄的大片| 午夜激情欧美在线| 国产午夜精品久久久久久一区二区三区| 久久99热6这里只有精品| 亚洲精品日韩在线中文字幕 | 婷婷色综合大香蕉| 99热这里只有精品一区| 99热网站在线观看| 成人美女网站在线观看视频| 国产一区二区激情短视频| 免费一级毛片在线播放高清视频| 综合色av麻豆| 亚洲无线观看免费| 亚洲av第一区精品v没综合| 国产高清三级在线| 美女国产视频在线观看| 国产又黄又爽又无遮挡在线| 一本—道久久a久久精品蜜桃钙片 精品乱码久久久久久99久播 | 亚洲av中文av极速乱| 在线免费观看不下载黄p国产| 91午夜精品亚洲一区二区三区| 午夜激情欧美在线| 国产一区二区三区av在线 | 国产精品人妻久久久影院| 国产成人精品一,二区 | www.av在线官网国产| 亚洲精品国产成人久久av| 老女人水多毛片| 日韩欧美精品v在线| 春色校园在线视频观看| www.av在线官网国产| 超碰av人人做人人爽久久| 99久久精品一区二区三区| 色综合亚洲欧美另类图片| 91精品国产九色| 国产精品福利在线免费观看| 一进一出抽搐动态| 三级经典国产精品| 69人妻影院| 特级一级黄色大片| 成人漫画全彩无遮挡| 国产在线男女| 有码 亚洲区| 亚洲在线自拍视频| 久久精品国产亚洲av香蕉五月| 免费在线观看成人毛片| 精品熟女少妇av免费看| 一边摸一边抽搐一进一小说| www.av在线官网国产| 九九热线精品视视频播放| 亚洲国产精品久久男人天堂| 国产精品久久久久久亚洲av鲁大| 中国国产av一级| 亚洲精品乱码久久久v下载方式| 91久久精品电影网| 18禁裸乳无遮挡免费网站照片| 日韩成人伦理影院| 亚洲18禁久久av| 欧美高清性xxxxhd video| 午夜福利在线在线| 精品欧美国产一区二区三| 国产精品一区二区三区四区免费观看| 亚洲国产欧美人成| 又黄又爽又刺激的免费视频.| 免费观看a级毛片全部| 乱系列少妇在线播放| 99视频精品全部免费 在线| 中国国产av一级| 可以在线观看毛片的网站| 久久久久久久久大av| 国产精品蜜桃在线观看 | 亚洲熟妇中文字幕五十中出| 亚洲人成网站在线播| 久久99精品国语久久久| 精品一区二区免费观看| 日韩在线高清观看一区二区三区| 亚洲久久久国产精品| 国产一区二区三区av在线| 亚洲av免费高清在线观看| 五月开心婷婷网| 国产亚洲欧美精品永久| 久久精品国产a三级三级三级| 成人毛片60女人毛片免费| 国产精品国产三级专区第一集| 啦啦啦在线观看免费高清www| 久热久热在线精品观看| 婷婷色麻豆天堂久久| 亚洲精品日本国产第一区| 亚洲国产精品999| 99热国产这里只有精品6| 男人爽女人下面视频在线观看| 亚洲性久久影院| 亚洲精品乱久久久久久| 欧美 日韩 精品 国产| 交换朋友夫妻互换小说| 卡戴珊不雅视频在线播放| 黄色欧美视频在线观看| 26uuu在线亚洲综合色| 国产精品一区www在线观看| 一级毛片黄色毛片免费观看视频| 又粗又硬又长又爽又黄的视频| 另类精品久久| 成人毛片a级毛片在线播放| 一二三四中文在线观看免费高清| 日韩成人伦理影院| 欧美另类一区| 国产成人精品一,二区| 亚洲伊人久久精品综合| 天堂俺去俺来也www色官网| 纵有疾风起免费观看全集完整版| xxxhd国产人妻xxx| 国产日韩欧美在线精品| 亚洲av.av天堂| 亚洲av不卡在线观看| 22中文网久久字幕| 日韩免费高清中文字幕av| 国产精品人妻久久久久久| 精品久久久久久电影网| 飞空精品影院首页| 国产精品国产三级专区第一集| 免费日韩欧美在线观看| 高清av免费在线| 免费日韩欧美在线观看| 欧美日韩综合久久久久久| 亚洲av电影在线观看一区二区三区| 国产免费又黄又爽又色| 久久久久久伊人网av| 99久久人妻综合| 一区在线观看完整版| 一级黄片播放器| 伊人亚洲综合成人网| av不卡在线播放| 夫妻性生交免费视频一级片| 日韩成人av中文字幕在线观看| 国产成人午夜福利电影在线观看| 欧美少妇被猛烈插入视频| 夜夜爽夜夜爽视频| 又大又黄又爽视频免费| 亚洲精品456在线播放app| 久久久亚洲精品成人影院| 精品一区二区三区视频在线| 大陆偷拍与自拍| 卡戴珊不雅视频在线播放| 国产一区二区三区综合在线观看 | 亚洲国产av影院在线观看| 女的被弄到高潮叫床怎么办| 制服诱惑二区| 97精品久久久久久久久久精品| 亚洲欧美一区二区三区黑人 | 国产女主播在线喷水免费视频网站| 亚洲欧美清纯卡通| 成人亚洲欧美一区二区av| 卡戴珊不雅视频在线播放| 黄色欧美视频在线观看| 久久精品国产自在天天线| 麻豆精品久久久久久蜜桃| 九九在线视频观看精品| a级毛片免费高清观看在线播放| 亚洲少妇的诱惑av| 亚洲精品aⅴ在线观看| 国产成人免费无遮挡视频| 熟女av电影| 欧美bdsm另类| 久久久久精品性色| 久久久国产精品麻豆| 欧美精品国产亚洲| 亚洲成人一二三区av| 国产永久视频网站| 最近手机中文字幕大全| 久久久久国产网址| 欧美性感艳星| 成人午夜精彩视频在线观看| 国产精品国产三级国产av玫瑰| 亚洲精品第二区| 人妻夜夜爽99麻豆av| 一级毛片我不卡| 看十八女毛片水多多多| 9色porny在线观看| 一个人免费看片子| 极品人妻少妇av视频| 黑人巨大精品欧美一区二区蜜桃 | 精品国产一区二区久久| 亚洲国产成人一精品久久久| av福利片在线| 国产爽快片一区二区三区| 最新中文字幕久久久久| 内地一区二区视频在线| 美女主播在线视频| 亚洲国产成人一精品久久久| 街头女战士在线观看网站| 精品少妇内射三级| 日韩视频在线欧美| 亚洲av男天堂| 男女无遮挡免费网站观看| 免费久久久久久久精品成人欧美视频 | av福利片在线| 91精品国产九色| 黄片无遮挡物在线观看| 老司机影院成人| 高清黄色对白视频在线免费看| 欧美日韩av久久| 国产精品 国内视频| 狂野欧美激情性bbbbbb| 国产永久视频网站| 啦啦啦啦在线视频资源| 婷婷色综合大香蕉| 自线自在国产av| 各种免费的搞黄视频| 精品视频人人做人人爽| 有码 亚洲区| a级片在线免费高清观看视频| 在线观看美女被高潮喷水网站| 国产在线一区二区三区精| 午夜福利网站1000一区二区三区| 国产欧美另类精品又又久久亚洲欧美| 两个人免费观看高清视频| 日日爽夜夜爽网站| 人妻夜夜爽99麻豆av| 亚洲精品第二区| 国产 一区精品| 久久久久久久大尺度免费视频| 26uuu在线亚洲综合色| 国产高清三级在线| 日韩强制内射视频| 九草在线视频观看| 91精品国产国语对白视频| 国产高清不卡午夜福利| 一本大道久久a久久精品| 亚洲av欧美aⅴ国产| 久久精品国产亚洲网站| 街头女战士在线观看网站| 亚洲不卡免费看| 国产精品国产三级国产av玫瑰| 日本欧美视频一区| 日韩中文字幕视频在线看片| 亚洲欧美色中文字幕在线| 有码 亚洲区| 在线播放无遮挡| 日本爱情动作片www.在线观看| 日韩熟女老妇一区二区性免费视频| 一本大道久久a久久精品| 日韩一本色道免费dvd| 建设人人有责人人尽责人人享有的| 国产精品一国产av| 高清在线视频一区二区三区| 亚洲国产欧美在线一区| 国产精品国产av在线观看| 午夜免费鲁丝|