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

    非洲栽培稻堊白粒率耐熱性QTL的定位

    2020-03-31 14:07:00曹志斌李瑤曾博虹毛凌華蔡耀輝吳曉峰袁林峰
    中國水稻科學(xué) 2020年2期
    關(guān)鍵詞:粒率表型染色體

    曹志斌 李瑤 曾博虹 毛凌華 蔡耀輝 吳曉峰,* 袁林峰,*

    非洲栽培稻堊白粒率耐熱性QTL的定位

    曹志斌1, #李瑤2, #曾博虹1毛凌華1蔡耀輝1吳曉峰1,*袁林峰1,*

    [1江西省超級水稻研究發(fā)展中心/國家水稻工程實(shí)驗(yàn)室(南昌), 南昌 330200;2江西省農(nóng)業(yè)科學(xué)院 土壤肥料與資源環(huán)境研究所/農(nóng)業(yè)部長江中下游作物生理生態(tài)與耕作重點(diǎn)實(shí)驗(yàn)室, 南昌 330200;#共2同第一作者;*通信聯(lián)系人, E-mail: yuanlf2003@sina.com; jxnkzy208@163.com]

    【】本研究旨在定位一個(gè)稻米堊白粒率高溫耐性QTL,為外觀品質(zhì)育種及解析堊白粒率高溫耐性的遺傳機(jī)制提供依據(jù)。以非洲栽培稻耐熱品種IRGC102309(Steud.)和秈稻品種R9311(L. subsp.Kato)為親本構(gòu)建的栽培稻種間染色體片段導(dǎo)入系CSIL05-23為材料構(gòu)建次級分離群體,結(jié)合人工氣候室模擬灌漿期高溫脅迫處理,采用堊白粒率高溫鈍感值為評價(jià)指標(biāo),對非洲栽培稻堊白粒率高溫耐性 QTL 進(jìn)行檢測。在BC6F2分離群體,利用單標(biāo)記分析,發(fā)現(xiàn)第5染色體上的SSR標(biāo)記RM1200與堊白粒率耐熱性狀極顯著正相關(guān)(=0.0005)。進(jìn)一步利用BC6F3和BC6F4分離群體,采用QTL Cartographer 2.5軟件和復(fù)合區(qū)間作圖法在水稻第5染色體上的SSR標(biāo)記RM1200?RM5796區(qū)間重復(fù)檢測到一個(gè)灌漿期堊白粒率耐熱性QTL, 命名為,分別解釋11.4%和17.5%表型變異。根據(jù)BC6F4分離群體的純合重組體表型分組,利用置換作圖方法將目標(biāo)QTL同樣定位在SSR標(biāo)記RM1200?RM5796之間,遺傳圖距為1.3 cM,物理圖距約為333.4 kb??刂茍装琢B誓蜔嵝缘氖且粋€(gè)能夠用于稻米外觀品質(zhì)育種的新QTL。

    水稻;堊白粒率;耐熱性;數(shù)量性狀位點(diǎn);基因定位

    水稻是全球最重要的糧食作物之一,世界上將近一半人口的能量供應(yīng)與之直接相關(guān)[1]。未來全球氣候變暖的趨勢對全球水稻生產(chǎn)不利影響將日益加劇,熱帶及亞熱帶地區(qū)作物生長季節(jié),氣溫可能會(huì)提高3~5℃,將會(huì)超過有記錄的作物生長季節(jié)性高溫極值[2]。而且高溫?zé)岷μ鞖饨?jīng)常與水稻生殖生長關(guān)鍵期相遇,在孕穗期、揚(yáng)花期以及灌漿期均有發(fā)生。因此,培育能夠適應(yīng)未來氣候變化的水稻品種,是應(yīng)對全球氣候變暖的迫切需要[3]。灌漿期對水稻籽粒產(chǎn)量與品質(zhì)形成極為關(guān)鍵。當(dāng)灌漿期遇到炎熱的氣候條件,對優(yōu)質(zhì)稻米品質(zhì)的形成很不利,當(dāng)氣溫超過35℃,籽粒的灌漿速度加快,灌漿時(shí)間縮短,出現(xiàn)高溫逼熟現(xiàn)象。產(chǎn)量方面的不利影響,主要表現(xiàn)為籽粒充實(shí)度差,千粒重顯著降低。稻米品質(zhì)方面的不利影響,主要表現(xiàn)為堊白粒率高,堊白大,整精米率下降,食味品質(zhì)變差[4-10]。

    水稻灌漿期耐熱性屬于受多基因控制的復(fù)雜數(shù)量性狀,受外界環(huán)境影響明顯[11-16]。目前已有學(xué)者利用不同作圖群體及鑒定指標(biāo)結(jié)合分子標(biāo)記技術(shù),對灌漿期耐熱性相關(guān)QTL進(jìn)行了初步定位。朱昌蘭等[17]利用Nipponbare/Kasalath//Nipponbare回交重組自交系群體,以直鏈淀粉含量為鑒定指標(biāo),在第6和9染色體上檢測到3個(gè)QTL。以膠稠度為鑒定指標(biāo),在第4、6、7、8、10和11染色體上分別定位到1個(gè)水稻灌漿期耐熱相關(guān)QTL。朱昌蘭等[18]利用Nipponbare/Kasalath//Nipponbare回交重組自交系群體,在第1、4和7染色體上分別定位了一個(gè)灌漿期粒重為耐熱指標(biāo)的QTL。Kobayashi 等[19]利用Hana-echizen/Niigatawase構(gòu)建的F3群體,以稻米外觀品質(zhì)性狀背白和腹白為灌漿期耐熱性鑒定指標(biāo),共檢測到3個(gè)QTL,其中在第4及第6染色體上分別鑒定到1個(gè)白背耐熱QTL,在第6染色上定位到一個(gè)白背和基白QTL。Shirasawa等[20]利用Kokoromachi/Tohoku 168構(gòu)建的重組群體在第1和6染色體上個(gè)檢測到1個(gè)白背耐熱性QTL。Tabata等[21]利用Koshijiwase/Chiyonishiki構(gòu)建的重組自交系群體在第1、2和4染色體上共檢測到4個(gè)稻米白背耐熱性QTL。張桂蓮等[22]以耐熱水稻品系 996 和熱敏感水稻品系4628為親本構(gòu)建的重組自交系為材料,在第1、3、6、7染色體上檢測到8個(gè)堊白粒率高溫鈍感性QTL,在第1、2、4、6、7和10染色體上檢測到12個(gè)堊白大小高溫鈍感性QTL,在第1、3、7染色體上檢測到4個(gè)堊白度高溫鈍感性QTL。

    目前這些定位結(jié)果主要是利用亞洲栽培稻秈稻種內(nèi)或秈粳亞種間雜交構(gòu)建的初級作圖群體分析得到的。本研究利用非洲栽培稻耐熱品種為供體親本,秈稻品種R9311為受體親本,獲得堊白粒率為17.8%的染色體片段導(dǎo)入系CSIL05-23。并通過進(jìn)一步構(gòu)建次級分離群體結(jié)合人工氣候室進(jìn)行表型鑒定,結(jié)合復(fù)合區(qū)間作圖及置換作圖法在第5染色體上定位到一個(gè)主效堊白粒率高溫鈍感 QTL。這一工作可為研究堊白粒率高溫鈍感性的遺傳機(jī)理,篩選與堊白粒率高溫鈍感性緊密連鎖的分子標(biāo)記以及下一步的育種應(yīng)用及圖位克隆打下基礎(chǔ)。

    1 材料與方法

    1.1 非洲栽培稻堊白粒率耐熱性QTL分析材料

    利用亞洲栽培秈稻品種R9311作為受體,非洲栽培稻品種IRGC102309作為供體,2008年夏季開始連續(xù)回交5次,自交1次獲得BC5F2種子。2010年種植BC5F2群體400株,挑選20株綜合農(nóng)藝性狀良好的單株自交套袋收種。2011年冬季將中選材料海南三亞育種試驗(yàn)站加代繁殖獲得BC5F3群體。

    2012年于南昌試驗(yàn)站中稻季節(jié)灌漿期(8月16~26日)經(jīng)歷自然高溫脅迫,日平均溫度不低于32℃,最高溫度不低于36℃。種植BC5F3群體134行,每行8株。成熟后單株收獲,脫殼調(diào)查堊白性狀,以播種時(shí)間推遲15 d,9月初開始灌漿的株系的堊白性狀值作為正常溫度下的表型對照。其中編號為CSIL05-23的染色體片段導(dǎo)入系,高溫脅迫處理下堊白粒率平均值為17.8%,而在正常溫度下堊白粒率平均值為17.1%,鈍感值為1.041,比值接近1,兩者差異不顯著,對高溫表現(xiàn)為強(qiáng)鈍感性,同時(shí)自交留種,冬季到海南繁殖得到BC5F4種子。將含有目的片段的染色體片段導(dǎo)入系CSIL05-23,分別于2013年、2014年、2015年種植在南昌高安基地,種植4行,行長1.2 m,株距16.5 cm,行距19.8 cm,隨機(jī)區(qū)組設(shè)計(jì),設(shè)2次重復(fù)。2014年夏季在南昌試驗(yàn)站種植200個(gè)BC6F2單株。2015年種植RM1200和RM5796片段雜合的BC6F2單株自交獲得的BC6F3分離群體380株,2016年種植目標(biāo)區(qū)段雜合的BC6F3單株自交獲得的BC6F4群體420株。上述分離群體單株,待幼穗分化至6期左右,帶泥移栽至盆缽中繼續(xù)生長發(fā)育。將掛牌標(biāo)識(shí)揚(yáng)花的稻穗于花后7 d開始用人工氣候室模擬高溫脅迫14 d。然后移入正常自然條件下成熟5~7 d,調(diào)查掛牌標(biāo)識(shí)稻穗堊白粒率。

    1.2 人工氣候室高溫脅迫處理

    人工氣候室(北京易盛泰和科技)模擬高溫脅迫的條件為8:00–17:00,控制在38.0℃±0.5℃,其余時(shí)間均為25.0℃±0.5℃。正常對照條件為8:00–17:00控制在32.0℃±0.5℃,其余時(shí)間均為25.0℃±0.5℃。所有處理的相對濕度均為75%±5%;光照時(shí)間為14 h,光照強(qiáng)度為20 000~25 000 lx。

    1.3 堊白粒率高溫耐性表型測定

    利用JSE-II稻米堊白觀測儀考察表型。取整精米100粒,置于觀測板上觀察,挑選有堊白(包括心白、腹白、背白)的米粒,計(jì)平均值作為堊白粒率。以堊白粒率高溫鈍感值(insensitivity value)表示堊白粒率耐熱性。堊白高溫鈍感值的計(jì)算參照鐘秉強(qiáng)等[23]方法,鈍感值(IV)=1/2,其中1,2分別為高溫脅迫下和正常溫度條件下的堊白粒率。從理論上講,鈍感值越接近1,該性狀越耐高溫。

    1.4 DNA提取,PCR擴(kuò)增及多態(tài)性標(biāo)記篩選

    應(yīng)用CTAB法[24]提取DNA,本實(shí)驗(yàn)室根據(jù)Gramene數(shù)據(jù)庫(http://www.gramene.org)選擇SSR引物進(jìn)行PCR擴(kuò)增,擴(kuò)增體系為10 μL,DNA模板l μL。

    PCR擴(kuò)增條件為:94℃下預(yù)變性5 min;94℃下變性30 s,57℃下復(fù)性30 s,72℃下延伸1 min;35個(gè)循環(huán)后,72℃下再延伸10 min。擴(kuò)增產(chǎn)物經(jīng)非變性聚丙烯酞胺凝膠電泳,凝膠濃度為6%,電泳緩沖液為0.5倍TBE,220 V恒壓電泳1 h,并銀染顯色。

    1.5 QTL分析

    根據(jù)Mapmaker/Exp 3.0b軟件[25]分析多態(tài)性標(biāo)記之間的連鎖遺傳關(guān)系, 構(gòu)建局部連鎖遺傳圖譜。設(shè)定最小LOD值為3.0, 根據(jù)Kosambi函數(shù)[26]將分子標(biāo)記重組率換算成遺傳圖距(cM), 遺傳圖距最大值為50 cM。利用Windows QTL Cartographer 2.5[27]檢測目的性狀QTL。利用復(fù)合區(qū)間作圖法(composite interval mapping, CIM)[28],取LOD值2.5為 QTL存在的閾值,利用模型6進(jìn)行分析。用SAS統(tǒng)計(jì)軟件對標(biāo)記與性狀進(jìn)行單標(biāo)記相關(guān)性分析。根據(jù)多態(tài)性分子標(biāo)記對群體單株的基因型檢測結(jié)果將數(shù)據(jù)分組。利用SAS統(tǒng)計(jì)軟件中方差相同的測驗(yàn)檢驗(yàn)組間平均數(shù)的差異,確定單標(biāo)記與性狀的相關(guān)性[29]。通過差異顯著性分析,判斷控制數(shù)量性狀的QTL是否與標(biāo)記連鎖。2是標(biāo)記能夠解釋性狀表型變異的比例,由表型性狀與標(biāo)記的回歸方程中的決定系數(shù)決定[30]。QTL命名方法參照McCouch等制定的規(guī)則[31]。

    2 結(jié)果與分析

    2.1 染色體片段導(dǎo)入系CSIL05-23高溫脅迫及常溫條件下堊白粒率表現(xiàn)

    CSIL05-23在田間高溫脅迫和人工氣候室模擬高溫脅迫條件下堊白粒率分別為17.8%和18.4%,顯著低于受體親本R9311,R9311分別為55.5%和56.8%(圖1-A)。CSIL05-23在田間高溫脅迫和人工氣候室模擬高溫脅迫條件下堊白粒率高溫鈍感值分別為1.04和1.07,R9311分別為2.25和2.30(圖1-B)。

    2.2 堊白粒率耐熱性QTL的初步鑒定

    由表1可知,BC6F2群體中堊白粒率高溫鈍感值表現(xiàn)出極顯著差異,且變異幅度較大,存在顯著的超親分離。因此推測非洲栽培稻片段含有降低堊白粒率高溫鈍感值的優(yōu)異等位基因,且這些優(yōu)異位點(diǎn)已在亞洲栽培秈稻中發(fā)生了重組并引起了群體表型的變異。偏度和峰度的絕對值均小于1,認(rèn)為堊白粒率高溫鈍感值表型均符合正態(tài)分布,呈多基因控制的數(shù)量性狀遺傳特點(diǎn)。

    從Gramene數(shù)據(jù)庫(http://www. gramene.org) 中選擇均勻覆蓋水稻基因組的430對多態(tài)性SSR引物對CSIL05-23和R9311進(jìn)行PCR擴(kuò)增,獲得3對具多態(tài)性SSR標(biāo)記,分別為RM488(Chr. 1),RM5460(Chr. 2),RM1200(Chr. 5)。利用這3個(gè)SSR標(biāo)記對200個(gè)BC6F2分離單株進(jìn)行檢測,發(fā)現(xiàn)3個(gè)標(biāo)記基因型均發(fā)生分離。根據(jù)分子標(biāo)記結(jié)果將數(shù)據(jù)分組,利用方差相同的測驗(yàn)檢驗(yàn)組間平均數(shù)的差異,確定標(biāo)記與堊白粒率高溫鈍感性狀的連鎖關(guān)系,發(fā)現(xiàn)第5染色體上RM1200與堊白粒率高溫鈍感值呈極顯著正相關(guān)(=0.0005)。該標(biāo)記位點(diǎn)的非洲栽培稻增效基因能夠顯著提高堊白粒率高溫鈍感性,解釋9.4%的表型變異。

    **差異極顯著(p<0.01), t檢驗(yàn)。

    Fig. 1. Comparison of heat tolerance of chalky grain rate of CSIL05-23 and R9311 in field and artificial climatic chambers.

    表1 親本及BC6F2、BC6F3、BC6F4群體堊白粒率高溫鈍感值的表型變異

    表2 BC6F3與BC6F4世代灌漿期高溫脅迫后堊白粒率耐熱性QTL分析

    IV indicates insensitive value of heat tolerance of chalky grain rate. X1and X2indicate chalky grain rate under heat stress condition and normal condition, respectively.

    2.3 堊白粒率耐高溫QTL復(fù)合區(qū)間作圖

    為了進(jìn)一步消除遺傳背景的影響,有效確定目標(biāo)QTL的位置和遺傳效應(yīng)。在BC6F2分離群體中利用RM488、RM5460和RM1200進(jìn)行滲入片段標(biāo)記基因型鑒定,將RM488和RM5460基因型為輪回親本純合基因型,而目標(biāo)滲入片段RM1200為雜合基因型的BC6F2單株自交,獲得BC6F3和BC6F4次級分離群體進(jìn)行分析。如表1所示,群體間差異達(dá)到極顯著水平并存在明顯的超親分離現(xiàn)象。2個(gè)世代堊白粒率耐熱性狀的鈍感值均符合正態(tài)分布,呈多基因控制的數(shù)量性狀遺傳特點(diǎn)。

    根據(jù)Gramene數(shù)據(jù)庫, 在第5染色體目標(biāo)區(qū)域附近挑選并合成25對SSR引物進(jìn)行多態(tài)性篩選, 共獲得8對多態(tài)性SSR標(biāo)記。將這些具多態(tài)性的引物進(jìn)一步分析整個(gè)BC6F3群體,應(yīng)用Mapmaker/Exp 3.0軟件構(gòu)建分子遺傳圖譜,該漸滲片段長10.3 cM,標(biāo)記間平均距離為1.5 cM(圖2)。

    應(yīng)用 QTL Cartographer 2.5軟件,采用復(fù)合區(qū)間作圖法對BC6F3單株及BC6F4群體堊白粒率高溫鈍感值,高溫脅迫處理下堊白粒率、正常溫度處理下堊白粒率分別進(jìn)行了QTL分析(表2),三個(gè)指標(biāo)性狀在2個(gè)世代相同標(biāo)記RM1200?RM5796區(qū)間都定位到目標(biāo)QTL,增效基因來源于非洲栽培稻。針對目標(biāo)QTL(a QTL for heat tolerance of chalky grain rate)以堊白粒率高溫鈍感值為鑒定指標(biāo)進(jìn)行QTL分析,發(fā)現(xiàn)在2個(gè)世代分別可以解釋11.3%,17.5%的表型變異。以高溫脅迫處理下堊白粒率X1為鑒定指標(biāo)進(jìn)行QTL分析,發(fā)現(xiàn)在2個(gè)世代分別可以解釋9.5%,16.7%的表型變異。以高溫脅迫處理下堊白粒率為鑒定指標(biāo)進(jìn)行QTL分析,發(fā)現(xiàn)在2個(gè)世代分別可以解釋9.5%,16.7%的表型變異。正常溫度處理下堊白粒率為鑒定指標(biāo)進(jìn)行QTL分析,發(fā)現(xiàn)在2個(gè)世代分別可以解釋12.6%、16.5%的表型變異。

    左右箱線圖分別表示來自BC6F3和BC6F4群體的1 LOD 和2 LOD 似然區(qū)間位置。

    Fig. 2. Genetic linkage map and likelihood intervals for QTL associated with heat tolerance of chalky grain rate.

    2.4 利用置換作圖方法定位

    利用上述8個(gè)SSR標(biāo)記對430個(gè)BC6F4單株分析,在RM1200和RM5796之間總共得到90個(gè)交換單株,通過連鎖分析構(gòu)建了目標(biāo)區(qū)域染色體片段置換圖(圖3)。根據(jù)交換單株標(biāo)記基因型信息,將90個(gè)交換單株中部分純合體按交換位置信息分為14組。利用置換作圖的方法分別利用高溫脅迫處理下堊白粒率(X1)、正常溫度處理下堊白粒率(X2)、堊白粒率高溫鈍感值(IV)三個(gè)指標(biāo)性狀對進(jìn)行定位。重組體G1~G7分別攜帶RM592到RM4777逐漸遞減的IRGC102309純合漸滲片段。G1~G6純合重組體與R9311之間比較在X1、X2和IV性狀上均存在顯著差異(<0.05),與CSIL05-23比較在X1、X2和IV性狀上差異不顯著。G7純合重組體與R9311之間比較在X1、X2和IV性狀上差異不顯著,與CSIL05-23比較在X1、X2和IV性狀上差異顯著(<0.05)。因此,根據(jù)G1~G6和G7純合重組體與對照之間的差異顯著性分析結(jié)果把堊白粒率高溫鈍感QTL界定在RM1200的下游。G8~G14分別攜帶從RM1200到RM405逐漸遞減的IRGC102309純合漸滲片段。其中,G8純合重組體與R9311比較,在X1、X2和IV性狀上差異顯著(<0.05),與CSIL05-23比較在X1、X2和IV性狀上差異不顯著。G9~G14純合重組體與R9311比較,在1、2和IV性狀上差異不顯著,與CSIL05-23比較在1、2和IV性狀上差異顯著(<0.05)。因此,G8和G9~G14純合重組體與對照親本比較結(jié)果把堊白粒率高溫鈍感QTL界定在RM5796的上游。這些研究結(jié)果表明位于SSR標(biāo)記RM1200–RM5796之間,遺傳距離為1.3 cM,物理距離約為333.4 kb。

    由380個(gè)BC6F3單株的分析數(shù)據(jù)構(gòu)建QTL區(qū)域的連鎖圖。BC6F4純合重組體的后代在灌漿期模擬高溫逆境處理后, 根據(jù)堊白粒率及鈍感值(IV=X1/X2)調(diào)查結(jié)果將qHTCGR5定位在RM1200和RM5796之間, 根據(jù)基因型將90個(gè)重組體分為14組。每組的重組體數(shù)目及與R9311和染色體片段導(dǎo)入系CSIL05-23之間的堊白粒率高溫鈍感表型差異顯著性在右邊標(biāo)出。“a”表示重組體與R9311的表型值在0.05水平無顯著差異:“b”表示重組體與CSIL05-23的表型值在0.05水平無顯著差異。

    Fig. 3. Mapping ofby a substitution mapping strategy.

    3 討論

    非洲栽培稻與亞洲栽培稻同屬AA基因組的不同種,兩者在進(jìn)化中屬于獨(dú)立起源種,并且存在明顯的生殖隔離現(xiàn)象。需要構(gòu)建染色體片段導(dǎo)入系或近等基因系來克服種間不育問題。目前通過傳統(tǒng)回交結(jié)合現(xiàn)代生物技術(shù)手段已經(jīng)能夠有效挖掘利用非洲栽培稻的種質(zhì)特性。非洲栽培稻在當(dāng)?shù)刈匀贿x擇及進(jìn)化過程中,已經(jīng)具備適應(yīng)或抵抗高溫?zé)岷?、澇害、干旱、耐鹽堿及各種病蟲害特性,對亞洲栽培稻品種改良具有重要意義。本研究利用構(gòu)建的QTL-CSIL群體為鑒定堊白粒率耐熱性QTL提供了豐富的遺傳信息。除了證實(shí)控制目的性狀的非洲栽培稻等位基因,在不同世代的遺傳效應(yīng)穩(wěn)定及方向一致外,而且通過消除背景噪音將的遺傳效應(yīng)值大小進(jìn)行準(zhǔn)確的評估。很顯然利用復(fù)合區(qū)間作圖法分析F2群體對這些QTL遺傳效應(yīng)估計(jì)的準(zhǔn)確性是不足的,實(shí)際上目前水稻耐熱QTL克隆的成功例子也是得益于高質(zhì)量近等基因系的培育[32]。

    在籽粒灌漿期遇到高溫脅迫會(huì)導(dǎo)致灌漿速率加快,導(dǎo)致籽粒堊白粒率、堊白大小、堊白度顯著增加。堊白相關(guān)性狀受多個(gè)基因/QTL 位點(diǎn)協(xié)同調(diào)控表達(dá),并且與外界環(huán)境因素存在明顯的互作效應(yīng)。目前第5染色體上堊白性狀QTL的定位研究已有報(bào)道,且都是在田間自然條件下鑒定得到的。Liu等[33]利用2個(gè)粳稻品種Koshihikari和C602雜交、自交構(gòu)建了261個(gè)重組自交系(F8)進(jìn)行堊白性狀QTL分析,在第5染色體上鑒定到一個(gè)堊白粒率QTL,定位在RM289和RM3437之間,解釋18.2%的表型變異。高方遠(yuǎn)等[34]利用秈稻保持系骨干親本崗46B與美國水稻品種Lemont多次回交培育的低堊白、窄粒株系K1075與崗46B構(gòu)建F2群體, 對堊白粒率進(jìn)行QTL分析,將第5染色體上的堊白粒率QTL定位在RM18004?RM18068之間,該QTL解釋52.55%的表型變異。上述2個(gè)堊白粒率QTL與區(qū)間標(biāo)記位置均不重疊。因此是一個(gè)新的耐熱性QTL位點(diǎn)。該位點(diǎn)在BC6F3和BC6F4兩個(gè)世代被定位在同一個(gè)標(biāo)記區(qū)間,遺傳距離為1.3 cM,并且遺傳效應(yīng)值較大。這對水稻的生殖期高溫脅迫條件下外觀品質(zhì)的改良非常有意義。

    置換作圖法最早由Paterson提出[35],并成功應(yīng)用于QTL精細(xì)定位。這種方法已經(jīng)證明置換作圖是一種區(qū)分基因多效性及排除連鎖累贅的有力工具,可以為QTL的精細(xì)定位及圖位克隆打下堅(jiān)實(shí)的基礎(chǔ)[36-37]。CSIL05-23所含漸滲區(qū)域已被證實(shí)能夠改良受體的堊白粒率高溫耐性。在本研究BC6F4群體的置換作圖中,我們對純合交換片段都進(jìn)行了分析,對不同交換片段的作用效果進(jìn)行了評估。同時(shí)SSR標(biāo)記在確定重組位置及基因滲入?yún)^(qū)域邊界是非常有效的。與受體親本相比,通過置換作圖發(fā)現(xiàn)1.3 cM的IRGC102309滲入片段上存在能夠在正常溫度下顯著降低受體親本堊白粒率X2值的QTL,在高溫脅迫下該片段上QTL也能顯著降低受體親本堊白粒率X1值,并且堊白粒率高溫鈍感值IV接近1,表現(xiàn)為堊白粒率高溫絕對鈍感(圖2)。由于本研究人工氣候室處理群體數(shù)量有限,不能用于精細(xì)定位的大規(guī)模分離群體的表型鑒定,因此本研究置換作圖分辨率仍然不能清楚地分辨出X1,X2和IV三個(gè)性狀是由于基因連鎖導(dǎo)致,還是基因的多效性導(dǎo)致這種現(xiàn)象發(fā)生。因此,下一步將設(shè)計(jì)精細(xì)定位實(shí)驗(yàn)方案,構(gòu)建高密度局部分子標(biāo)記遺傳圖譜,增加QTL精細(xì)定位分辨率,分析上述三個(gè)性狀是否由這個(gè)區(qū)域上的同一個(gè)QTL控制,還是三個(gè)QTL緊密連鎖影響,為目標(biāo)QTL的育種應(yīng)用及相關(guān)基因克隆打下堅(jiān)實(shí)的基礎(chǔ)。

    [1] Carriger S, Vallée D. More crop per drop[J]., 2007, 6(2): 10-13.

    [2] Hockley N, Gibbons J M, Edwards-Jones G. Risks of extreme heat and unpredictability[J]., 2009, 324(5924): 177-179.

    [3] Battisti D S, Naylor R L. Historical warnings of future food insecurity with unprecedented seasonal heat[J]., 2009, 323(5911): 240-244.

    [4] Fitzgerald M A, McCouch S R, Hall R D. Not just a grain of rice: the quest for quality[J]., 2009, 14(3): 133-139.

    [5] 森谷國男. 徐正進(jìn), 譯. 水稻高溫脅迫抗性遺傳育種研究概況[M]. 雜交水稻, 1992(1): 47-48.

    Sengu G N. Translated by Xu Z J. General Research on Genetic Breeding of Resistance to High Temperature Stress in Rice[M]., 1992(1): 47-48. (in Chinese with English abstract)

    [6] 李木英, 熊偉, 石慶華, 胡志紅, 潘曉華, 譚雪明. 高溫脅迫對早稻不同品種灌漿結(jié)實(shí)和稻米品質(zhì)的影響[J]. 江西農(nóng)業(yè)學(xué)學(xué)報(bào), 2006, 28(4): 483-487.

    Li M Y, Xiong W, Shi Q H, Hu Z H, Pan X H, Pan X M. Effect of high temperature stress on endosperm filling and grain quality of early rice varieties[J]., 2006, 28(4): 483-487. (in Chinese with English abstract)

    [7] 趙海燕, 姚鳳梅, 張勇, 徐賓, 袁靜, 胡亞南, 許吟隆. 長江中下游水稻開花灌漿期氣象要素與結(jié)實(shí)率和粒重的相關(guān)性分析[J]. 中國農(nóng)業(yè)科學(xué), 2006, 39(9): 1765-1771.

    Zhao H Y, Yao F M, Zhang Y, Xu B, Yuan J, Hu Y N, Xu Y L. Correlation analysis of rice seed setting rate and weight of 1000-grain and agro-meteorology over the middle and lower reaches of the Yangtze River[J]., 2006, 39(9): 1765-1771. (in Chinese with English abstract)

    [8] 李林, 沙國棟, 陸景淮. 水稻灌漿期溫光因子對稻米品質(zhì)的影響[J]. 中國農(nóng)業(yè)氣象, 1989, 10(3): 33-38.

    Li L, Sha G D, Lu J H. Effect of temperature and light on rice quality [J]., 1989, 10(3): 33-38. (in Chinese with English abstract)

    [9] 程方明, 張蒿午, 吳永常. 灌漿結(jié)實(shí)期溫度對稻米堊白形成的影響[J]. 西北農(nóng)業(yè)學(xué)報(bào), 1996, 5(2): 31-34. Cheng F M, Zhang H W, Wu Y C. Effect of high temperature stress on chalkiness at filling stage[J]., 1996, 5(2): 31-34.

    [10] Zhong L, Cheng F, Wen X, Sun Z X, Zhang G P. The deterioration of eating and cooking quality caused by high temperature during grain filling in early-season indica rice cultivars[J]., 2005, 191(3): 218-225.

    [11] Mei D Y, Zhu Y J, Yu Y H, Fan Y Y, Huang D R, Zhuang J Y. Quantitative trait loci for grain chalkiness and endosperm transparency detected in three recombinant inbred line populations of indica rice[J]., 2013, 12(1): 1-11.

    [12] 周立軍, 劉喜, 江玲, 鄭蕾娜, 陳亮明, 劉世家, 翟虎渠, 萬建民. 利用CSSL和BIL群體分析稻米堊白粒率QTL 及互作效應(yīng). 中國農(nóng)業(yè)科學(xué), 2009, 42(4): 1129-1135.

    Zhou L J, Liu X, Jiang L, Zheng L N, Chen L M, Liu S J, Zhai H Q, Wan J M. Analysis of QTL and GE effects on PGWC in rice (L.) using CSSL and BIL populations., 2009, 42(4): 1129-1135. (in Chinese with English abstract)

    [13] Liu X, Wang Y, Wang S W. QTL analysis of percentage of grains with chalkiness in japonica rice ()[J]., 2012, 11(1): 717-724.

    [14] 晁園, 馮付春, 高冠軍, 朱雪萍, 何予卿. 利用重組自交系群體定位水稻品質(zhì)相關(guān)性狀的QTL[J]. 華中農(nóng)業(yè)大學(xué)學(xué)報(bào), 2012, 31(4): 397-403.

    Chao Y, Feng F C, Gao G J, Zhu X P, He Y Q. Mapping quantitative trait loci for qualities of rice grains using a Recombinant inbred(RIL) population [J]., 2012, 31(4): 397-403. (in Chinese with English abstract)

    [15] 楊亞春, 倪大虎, 宋豐順, 李澤福, 易成新, 楊劍波. 不同生態(tài)地點(diǎn)下稻米外觀品質(zhì)性狀的QTL定位分析[J]. 中國水稻科學(xué), 2011, 25(1): 43-51.

    Yang Y C, Ni D H, Song F S, Li Z F, Yi C X, Chen J B. Identification of QTLs of rice appearance quality traits across different ecological sites[J]., 2011, 25(1): 43-51. (in Chinese with English abstract)

    [16] 王林森, 陳亮明, 王沛然, 王卓然, 鄭海, 馬宏陽, 江玲, 趙志剛, 萬建民. 利用高世代回交群體檢測水稻堊白相關(guān)性狀 QTL[J]. 南京農(nóng)業(yè)大學(xué)學(xué)報(bào), 2016, 39(2): 183-190.

    Wang L S, Chen L M, Wang P R, Wang Z R, Zheng L, Ma H Y, Jiang L, Zhao Z G, Wan J M. Detecting the QTL of rice chalkiness traits using advanced backcrossing population[J]., 2016, 39(2): 183-190. (in Chinese with English abstract)

    [17] 朱昌蘭. 稻低直鏈淀粉含量的遺傳及品質(zhì)形成對高溫耐性的QTL分析. 南京: 南京農(nóng)業(yè)大學(xué), 2004. Zhu C L. Identifying QTLs for thermo-tolerance of quality formation and inheritance of low amylose content in rice. Nanjing: Nanjing Agriculture University, 2004. (in Chinese with English abstract)

    [18] 朱昌蘭, 肖應(yīng)輝, 王春明, 江玲, 翟虎渠, 萬建民. 水稻灌漿期耐熱害的數(shù)量性狀基因位點(diǎn)分析[J]. 中國水稻科學(xué), 2005, 19(2): 117-121.

    Zhu C L, Xiao Y H, Wang C M, Jiang L, Zhai H Q, Wan J M. Mapping QTLs for heat tolerance during grain filling in rice[J]., 2005, 19(2): 117-121. (in Chinese with English abstract)

    [19] Kobayashi A, Bao G, Ye S, Tomita K. Detection of quantitative trait loci for white-back and basal white kernels under high temperature stress in japonica rice varieties[J]., 200, 57(2): 107-116.

    [20] Shirasawa, K. QTL analysis of high-temperature-stress tolerance in filling period based on rice grain quality[J]., 2006, 8(1): 155.

    [21] Tabata M, Hirabayashi H, Takeuchi Y, Ando I. Mapping of quantitative trait loci for the occurrence of white-back kernels associated with high temperatures during the ripening period of rice(L)[J]., 2007, 57(1): 47-52.

    [22] 張桂蓮, 廖斌, 唐文幫, 陳立云, 肖應(yīng)輝. 稻米堊白性狀對高溫耐性的QTL分析. 中國水稻科學(xué), 2017, 31(3): 257-264.

    Zhang GL, Liao B, Tang W B, Chen L Y, Xiao Y H. Identifying QTLs for thermo-tolerance of grain chalkiness trait in rice[J]., 2017, 31(3): 257-264. (in Chinese with English abstract)

    [23] 鐘秉強(qiáng), 楊正林, 冉啟良, 何光華. 美國水稻品種農(nóng)藝性狀和品質(zhì)性狀的溫度鈍感特性研究[J]. 中國農(nóng)學(xué)通報(bào), 2005, 21(2): 118-121.

    Zhong B Q, Yang Z L, Ran Q L, He g H. Study on temperature insensitivity of characters of agronomy and quality in the American rice variety[J]., 2005, 21(2): 118-121. (in Chinese with English abstract)

    [24] Rogers S O, Bendich A J. Extraction of DNA from plant tissues[M]. Plant molecular Biology Manual. Springer, Dordrecht, 1989: 73-83.

    [25] Lander E S, Green P, Abrahamson J, Barlow A, Daly M J, Lincoln S E, Newberg L A. MAPMAKER: An interactive computer package for constructing primary genetic linkage maps of experimental and natural populations[J]., 1987, 1(2): 174-181.

    [26] Kosambi D D. The estimation of map distances from recombination values[J]., 1943, 12(YRS 1943/5): 172-175.

    [27] Wang S C, Basten C J, Zeng Z B. Windows QTL Cartographer 2.5. Statistical Genetics, Raleigh, NC: North Carolina State, 2012.

    [28] Zeng Z B. Precision mapping of quantitative trait loci., 1994, 136: 1457-1468.

    [29] 王慧, 喻德躍, 吳巧娟, 蓋鈞鎰. 大豆對斜紋夜蛾抗生性基因的微衛(wèi)星標(biāo)記(SSR)的研究[J]. 大豆科學(xué), 2004, 23(2): 91-95.

    Wang H, Yu D Y, Wu Q J, Gai J Y. Characterization of resistance genes to cotton worm with SSR markers in soybean[J]., 2004, 23(2): 91-95. (in Chinese with English abstract)

    [30] 徐吉臣, 鄒亮星. 利用相關(guān)性分析鑒定與水稻根部性狀表達(dá)相關(guān)的分子標(biāo)記[J]. 遺傳學(xué)報(bào), 2002, 29(3): 245-249.

    Xu J C, Zou L X. Identification of molecular markers associated with rice root traits by correlation coefficient analysis[J]., 2002, 29(3): 245-249. (in Chinese with English abstract)

    [31] McCouch S R, Cho Y G, Yano M, Paul E, Blinstrub M, Morishima H, Kinoshita T. Report on QTL nomenclature[J]., 1997, 14: 11-13.

    [32] Li X M, Chao D Y, Wu Y, Huang X H, Chen K, Cui L G, Su L, Ye W W, Chen H, Chen H C, Dong N Q, Guo T, Shi M, Feng Q, Zhang P, Han B, Shan J X, Gao J P, Lin H X. Natural alleles of a proteasome α2 subunit gene contribute to thermos-tolerance and adaptation of African rice[J]., 2015, 47(7): 827.

    [33] Liu X, Wang Y, Wang S W. QTL analysis of percentage of grains with chalkiness in Japonica rice () [J]., 2012, 11(1): 717-724.

    [34] 高方遠(yuǎn), 邱玲, 陸賢軍, 任鄄勝, 吳賢婷, 任光俊, 曾禮華. 雜交秈稻骨干保持系崗46B稻谷粒形及堊白的QTL分析[J]. 中國水稻科學(xué), 2014, 28(3): 235-242.

    Gao F Y, Qiu L, Lu X J, Ren J S, Wu X T, Ren G J, Zeng L H. QTL analysis on grain shape and chalkiness of an elite maintainer line Gang 46B in hybrid rice(L.) [J]., 2014, 28(3): 235-242.

    [35] Paterson A H, Deverna J W, Lanini B, Tanksley S D. Fine mapping of quantitative trait loci using selected overlapping recombinant chromosomes, in an interspecies cross of tomato[J]., 1990, 124(3): 735-742.

    [36] Tanksley S D, Ganal M W, Martin G B. Chromosome landing: a paradigm for map-based gene cloning in plants with large genomes[J]., 1995, 11(2): 63-68.

    [37] Alpert K B, Tanksley S D. High-resolution mapping and isolation of a yeast artificial chromosome contig containing fw2.2: A major fruit weight quantitative trait locus in tomato[J]., 1996, 93(26): 15503-15507.

    QTL Mapping for Heat Tolerance of Chalky Grain Rate ofSteud.

    CAO Zhibin1, #, LI Yao2, #, ZENG Bohong1, MAO Linghua1, CAI Yaohui1, WU Xiaofeng1, *, YUAN Linfeng1, *

    (1Laboratory of Rice National Engineering, Jiangxi Research and Development Center of Super Rice, Nanchang 330200, China;2Institute of Soil and Fertilizer & Resources and Environment, Jiangxi Academy of Agricultural Sciences/Key Laboratory of Crop Ecophysiology and Farming System for the Middle and Lower Reaches of the Yangtze River, Ministry of Agriculture, Nanchang 330200, China;#These authors contributed equally to this work;*Corresponding author, E-mail: yuanlf2003@sina.com; jxnkzy208@163.com)

    【】The objective of the research is to identify a QTL for heat tolerance of chalky grain rate, so as to provide support forthe appearance quality breeding and analysis of genetic mechanism of rice high temperature tolerance of chalkygrain rate. 【】We developed a chalky grain rate heat-tolerant CSIL(chromosomal segment introgression lines), CSIL05-2, by backcrossing and marker assisted selection with African cultivated rice acc.IRGC102309 (Steud.) as donor parent and Asian cultivated riceR9311 (L. subsp.Kato) as recipient parent. And a QTL for insensitityof heat tolerance of chalky grain rate on chromosome 5 was analyzed using secondary populations from CSIL05-2.【】In a BC6F2segregation populations, the SSR marker RM1200 on chromosome 5 showed significant correlation with heat tolerance of chalky grain rateby single marker analysis(=0.0005). Using Cartographer 2.5 and the composite interval mapping with BC6F3and BC6F4populations, we further anchored the quantitative trait loci (QTL) associated with heat tolerance of chalky grain rate at filling stage to the same position within a 1.3cM interval, which was designed as qHTCGR and explained 11.4% and 17.5% of the phenotypic variances, respectively.Using homozygous recombinants screened from BC6F4populations,was also detected in a ~333.4 kb region betweenRM11633 and RM11642.【】The QTL controlling the heat tolerance of chalky grain rateis a novelQTL.

    rice; chalky grain rate; heat tolerance; QTL; gene mapping

    Q343.1+5; S511.32

    A

    1001-7216(2020)02-0135-08

    10.16819/j.1001-7216.2020.9086

    2019-07-26;

    2019-10-17。

    國家自然科學(xué)基金資助項(xiàng)目(31560382);江西省杰出青年科學(xué)基金資助項(xiàng)目(2018ACB21026);江西省自然科學(xué)基金資助項(xiàng)目(20151BAB214013; 2017BAB204020);江西省現(xiàn)代農(nóng)業(yè)科研協(xié)同創(chuàng)新專項(xiàng)(JXXTCXFY201902);國家重點(diǎn)研發(fā)計(jì)劃資助項(xiàng)目(2017YFD0301601; 2018YFD0301103);江西省重點(diǎn)研發(fā)計(jì)劃資助項(xiàng)目(20171ACF60019)。

    猜你喜歡
    粒率表型染色體
    水稻堊白主效QTL的定位與分析
    多一條X染色體,壽命會(huì)更長
    勘誤
    宿州市埇橋區(qū)2018年小麥不完善粒分析
    為什么男性要有一條X染色體?
    建蘭、寒蘭花表型分析
    能忍的人壽命長
    大豆新品種丹豆17
    GABABR2基因遺傳變異與肥胖及代謝相關(guān)表型的關(guān)系
    慢性乙型肝炎患者HBV基因表型與血清學(xué)測定的臨床意義
    日韩中字成人| 免费黄频网站在线观看国产| 大片免费播放器 马上看| 在线观看美女被高潮喷水网站| 校园人妻丝袜中文字幕| 免费av中文字幕在线| 中文精品一卡2卡3卡4更新| 在线观看三级黄色| 久久久国产一区二区| 一级a做视频免费观看| 国产一区亚洲一区在线观看| 国产成人91sexporn| 成人国产av品久久久| 国产伦理片在线播放av一区| 欧美激情国产日韩精品一区| 欧美最新免费一区二区三区| 国产av码专区亚洲av| 黄色日韩在线| 夫妻性生交免费视频一级片| av网站免费在线观看视频| freevideosex欧美| av免费观看日本| 久久久久久伊人网av| 国产亚洲欧美精品永久| 亚洲精品国产色婷婷电影| 久久99热这里只频精品6学生| 亚洲婷婷狠狠爱综合网| 国产淫片久久久久久久久| 亚洲av综合色区一区| 乱码一卡2卡4卡精品| 国产成人午夜福利电影在线观看| 人妻制服诱惑在线中文字幕| 国产成人精品久久久久久| 免费人妻精品一区二区三区视频| 菩萨蛮人人尽说江南好唐韦庄| 欧美另类一区| 毛片一级片免费看久久久久| 国产成人aa在线观看| 纯流量卡能插随身wifi吗| 精品久久久噜噜| 久久精品国产亚洲av涩爱| 亚洲第一区二区三区不卡| 黄片无遮挡物在线观看| 爱豆传媒免费全集在线观看| 国产精品女同一区二区软件| 日韩免费高清中文字幕av| 国产中年淑女户外野战色| 99久久人妻综合| 亚洲人成网站高清观看| 91狼人影院| 日本免费在线观看一区| 免费人妻精品一区二区三区视频| 一本色道久久久久久精品综合| 在线观看国产h片| 18禁裸乳无遮挡免费网站照片| 免费在线观看成人毛片| 中文在线观看免费www的网站| 97超碰精品成人国产| 色网站视频免费| 免费高清在线观看视频在线观看| 亚洲av男天堂| 国产精品精品国产色婷婷| 日韩欧美精品免费久久| 青春草国产在线视频| 极品少妇高潮喷水抽搐| 日韩成人伦理影院| 欧美xxⅹ黑人| 亚洲图色成人| 久久国产精品男人的天堂亚洲 | 国产乱人视频| 国产精品av视频在线免费观看| 一级毛片电影观看| 一区二区三区免费毛片| 欧美丝袜亚洲另类| 国产免费又黄又爽又色| 色综合色国产| 2022亚洲国产成人精品| 最近最新中文字幕免费大全7| 欧美一区二区亚洲| 日韩亚洲欧美综合| 免费黄网站久久成人精品| 赤兔流量卡办理| 黄色配什么色好看| 色网站视频免费| 2022亚洲国产成人精品| 日日撸夜夜添| 中国国产av一级| 精品久久久精品久久久| 亚洲成人一二三区av| videos熟女内射| 亚洲国产毛片av蜜桃av| 欧美精品国产亚洲| av又黄又爽大尺度在线免费看| 麻豆国产97在线/欧美| 99九九线精品视频在线观看视频| 日本爱情动作片www.在线观看| 国内精品宾馆在线| 国产淫片久久久久久久久| 蜜桃亚洲精品一区二区三区| 亚洲人成网站在线观看播放| 欧美精品人与动牲交sv欧美| 日日撸夜夜添| 啦啦啦中文免费视频观看日本| 亚洲国产精品成人久久小说| 99久久人妻综合| 久久久久久久久久成人| 大话2 男鬼变身卡| 国产真实伦视频高清在线观看| 免费不卡的大黄色大毛片视频在线观看| 激情五月婷婷亚洲| 色婷婷久久久亚洲欧美| 高清毛片免费看| 国产伦在线观看视频一区| h日本视频在线播放| 日产精品乱码卡一卡2卡三| 国产欧美另类精品又又久久亚洲欧美| 日韩成人av中文字幕在线观看| 熟女人妻精品中文字幕| 国产大屁股一区二区在线视频| 亚洲最大成人中文| 蜜臀久久99精品久久宅男| 国产成人精品婷婷| 一区二区三区乱码不卡18| 欧美日韩亚洲高清精品| 亚洲无线观看免费| 午夜福利高清视频| 日韩三级伦理在线观看| 嫩草影院入口| 91精品国产九色| 蜜桃在线观看..| 两个人的视频大全免费| 高清日韩中文字幕在线| 国产精品久久久久久av不卡| 热99国产精品久久久久久7| 又大又黄又爽视频免费| 观看av在线不卡| 国产高清不卡午夜福利| 欧美日韩一区二区视频在线观看视频在线| 亚洲av国产av综合av卡| 日本与韩国留学比较| 99热网站在线观看| 国产精品女同一区二区软件| 免费观看的影片在线观看| 蜜桃亚洲精品一区二区三区| 18+在线观看网站| www.av在线官网国产| 在线观看av片永久免费下载| 久久久午夜欧美精品| 我的女老师完整版在线观看| 免费av中文字幕在线| 国产av国产精品国产| 男人爽女人下面视频在线观看| 久久久久久久国产电影| 在线精品无人区一区二区三 | 欧美日韩视频高清一区二区三区二| av不卡在线播放| av女优亚洲男人天堂| 美女主播在线视频| 一级毛片电影观看| 久久青草综合色| 视频区图区小说| 国产精品精品国产色婷婷| 久久久色成人| 国产 一区 欧美 日韩| www.av在线官网国产| 国产黄频视频在线观看| 男女国产视频网站| 简卡轻食公司| 99热这里只有是精品50| av天堂中文字幕网| 欧美日韩精品成人综合77777| 国产精品精品国产色婷婷| 水蜜桃什么品种好| 国产日韩欧美在线精品| 熟女电影av网| 亚洲综合色惰| 男男h啪啪无遮挡| 国产v大片淫在线免费观看| 欧美人与善性xxx| 亚洲精品,欧美精品| 亚洲精品久久久久久婷婷小说| 各种免费的搞黄视频| 日日摸夜夜添夜夜爱| 性色av一级| 欧美高清成人免费视频www| 精品国产露脸久久av麻豆| 成人黄色视频免费在线看| 2018国产大陆天天弄谢| 国内精品宾馆在线| 亚洲国产毛片av蜜桃av| 99视频精品全部免费 在线| 久久精品久久久久久噜噜老黄| 午夜免费鲁丝| 人妻系列 视频| 80岁老熟妇乱子伦牲交| 1000部很黄的大片| 97超视频在线观看视频| 性色avwww在线观看| 又爽又黄a免费视频| 亚洲av二区三区四区| av播播在线观看一区| 一本色道久久久久久精品综合| 高清av免费在线| 亚洲欧美清纯卡通| 老司机影院成人| 国产日韩欧美亚洲二区| av福利片在线观看| 日韩精品有码人妻一区| 亚洲欧美日韩卡通动漫| 亚洲人与动物交配视频| 免费黄频网站在线观看国产| 插逼视频在线观看| 国产成人免费观看mmmm| 成人无遮挡网站| 亚洲美女搞黄在线观看| 欧美高清成人免费视频www| 国产成人精品久久久久久| 亚洲精品色激情综合| 99久久精品国产国产毛片| 天堂俺去俺来也www色官网| 久久久欧美国产精品| 男女啪啪激烈高潮av片| 国国产精品蜜臀av免费| 欧美xxxx性猛交bbbb| 日韩欧美一区视频在线观看 | 美女脱内裤让男人舔精品视频| 高清日韩中文字幕在线| 一级毛片我不卡| 女的被弄到高潮叫床怎么办| 欧美+日韩+精品| 深夜a级毛片| 少妇猛男粗大的猛烈进出视频| 国产爽快片一区二区三区| h日本视频在线播放| 精品一区二区免费观看| 一本久久精品| 亚洲国产毛片av蜜桃av| 日本黄大片高清| 有码 亚洲区| 秋霞在线观看毛片| 麻豆成人av视频| 深夜a级毛片| 久久精品久久精品一区二区三区| 久久久久久久大尺度免费视频| 国产在线免费精品| 男女下面进入的视频免费午夜| 久久 成人 亚洲| 亚洲不卡免费看| 最近中文字幕2019免费版| 亚洲精品色激情综合| 日韩亚洲欧美综合| 久久久精品免费免费高清| 蜜臀久久99精品久久宅男| 欧美性感艳星| 国产精品久久久久久精品古装| 一个人看视频在线观看www免费| 精品人妻视频免费看| 午夜福利视频精品| 国产大屁股一区二区在线视频| 久久人妻熟女aⅴ| 伦理电影免费视频| 久久99热6这里只有精品| 国产成人freesex在线| 日韩一本色道免费dvd| 国产男女内射视频| 午夜福利视频精品| 美女高潮的动态| 亚洲丝袜综合中文字幕| 大话2 男鬼变身卡| 99热网站在线观看| 少妇猛男粗大的猛烈进出视频| 久热久热在线精品观看| 亚洲,欧美,日韩| 午夜福利高清视频| 国产精品三级大全| 国产精品99久久99久久久不卡 | 黄片无遮挡物在线观看| 搡老乐熟女国产| 日韩免费高清中文字幕av| 免费少妇av软件| 久久久久久久久久久免费av| 国产 一区精品| 色视频在线一区二区三区| av女优亚洲男人天堂| 精品一区二区三卡| 全区人妻精品视频| 91aial.com中文字幕在线观看| 我要看日韩黄色一级片| 国产中年淑女户外野战色| 日韩电影二区| 交换朋友夫妻互换小说| a级毛片免费高清观看在线播放| 在线 av 中文字幕| 熟女av电影| 国内揄拍国产精品人妻在线| 在线 av 中文字幕| 日本黄大片高清| 能在线免费看毛片的网站| 国产一区二区三区综合在线观看 | av福利片在线观看| 亚洲高清免费不卡视频| 国产精品一区二区在线观看99| 超碰av人人做人人爽久久| 日韩欧美一区视频在线观看 | 国产片特级美女逼逼视频| 人人妻人人澡人人爽人人夜夜| 亚洲内射少妇av| videos熟女内射| 视频区图区小说| 国产精品伦人一区二区| 亚洲熟女精品中文字幕| 蜜桃亚洲精品一区二区三区| 久久国产乱子免费精品| 久热久热在线精品观看| 亚洲真实伦在线观看| 日韩 亚洲 欧美在线| 免费观看在线日韩| 能在线免费看毛片的网站| 女性生殖器流出的白浆| 国产成人a∨麻豆精品| 国产一区二区三区综合在线观看 | 99精国产麻豆久久婷婷| 精品久久久精品久久久| 国产精品久久久久久久电影| 18禁裸乳无遮挡免费网站照片| 成人亚洲精品一区在线观看 | 国产又色又爽无遮挡免| 欧美另类一区| 蜜桃在线观看..| 欧美另类一区| 欧美高清成人免费视频www| 亚洲国产欧美人成| 精品一品国产午夜福利视频| 少妇的逼好多水| 精品亚洲成a人片在线观看 | 亚洲欧美中文字幕日韩二区| 一级毛片aaaaaa免费看小| 成人毛片a级毛片在线播放| 久久ye,这里只有精品| 久久久久性生活片| 欧美+日韩+精品| 婷婷色综合大香蕉| 亚洲成人手机| 美女中出高潮动态图| 免费人妻精品一区二区三区视频| 国产av国产精品国产| 国产大屁股一区二区在线视频| 91久久精品电影网| 超碰av人人做人人爽久久| 久久精品国产鲁丝片午夜精品| 国产毛片在线视频| 天天躁夜夜躁狠狠久久av| 超碰av人人做人人爽久久| 老熟女久久久| 寂寞人妻少妇视频99o| 成人影院久久| 欧美丝袜亚洲另类| 网址你懂的国产日韩在线| 日日撸夜夜添| 午夜激情福利司机影院| 大香蕉久久网| 国产精品人妻久久久影院| 国产美女午夜福利| 国产精品不卡视频一区二区| 男女下面进入的视频免费午夜| 少妇猛男粗大的猛烈进出视频| 一级黄片播放器| 一个人免费看片子| 欧美xxⅹ黑人| 麻豆国产97在线/欧美| 男人爽女人下面视频在线观看| 观看免费一级毛片| 三级国产精品欧美在线观看| 亚洲国产欧美在线一区| 精品一区二区三区视频在线| 久久久a久久爽久久v久久| av在线蜜桃| 黑人猛操日本美女一级片| 午夜激情福利司机影院| 欧美老熟妇乱子伦牲交| 男女下面进入的视频免费午夜| 最后的刺客免费高清国语| 99国产精品免费福利视频| 色吧在线观看| 色视频www国产| 国产高清不卡午夜福利| 精品午夜福利在线看| 在线免费十八禁| 亚洲最大成人中文| 爱豆传媒免费全集在线观看| a级毛色黄片| 国产精品99久久久久久久久| 久久久久久伊人网av| 中文乱码字字幕精品一区二区三区| 日韩成人伦理影院| 免费观看a级毛片全部| 国产免费视频播放在线视频| 亚洲人成网站高清观看| 亚洲精品,欧美精品| 91精品国产国语对白视频| 国产免费又黄又爽又色| 中文欧美无线码| 简卡轻食公司| 亚洲激情五月婷婷啪啪| 久久久久国产精品人妻一区二区| 日韩人妻高清精品专区| 亚洲精品成人av观看孕妇| 日韩一本色道免费dvd| 国产淫片久久久久久久久| 亚洲第一区二区三区不卡| 久久国内精品自在自线图片| 国产免费一级a男人的天堂| 99久久精品国产国产毛片| 欧美日韩一区二区视频在线观看视频在线| 国产精品熟女久久久久浪| 免费观看在线日韩| 亚洲精品日本国产第一区| 高清黄色对白视频在线免费看 | 国内精品宾馆在线| 边亲边吃奶的免费视频| 欧美xxxx性猛交bbbb| 男女边摸边吃奶| 国产 精品1| 国产乱人视频| 国内精品宾馆在线| 亚洲色图av天堂| 狂野欧美激情性xxxx在线观看| 午夜精品国产一区二区电影| 亚洲激情五月婷婷啪啪| 免费观看的影片在线观看| 国产黄片美女视频| 欧美性感艳星| 亚洲av中文字字幕乱码综合| 精品午夜福利在线看| 亚洲国产精品国产精品| 免费观看a级毛片全部| 人妻少妇偷人精品九色| 一级毛片黄色毛片免费观看视频| 久久99热这里只频精品6学生| 欧美变态另类bdsm刘玥| 久久久成人免费电影| 国产精品蜜桃在线观看| av在线观看视频网站免费| 国产精品久久久久久精品电影小说 | 午夜日本视频在线| 夜夜骑夜夜射夜夜干| 黄色一级大片看看| 国产精品一区二区在线不卡| 少妇精品久久久久久久| 亚洲欧洲国产日韩| 热re99久久精品国产66热6| 欧美精品人与动牲交sv欧美| 中文字幕精品免费在线观看视频 | 黄片无遮挡物在线观看| 狂野欧美白嫩少妇大欣赏| 永久网站在线| 欧美成人一区二区免费高清观看| 免费大片黄手机在线观看| 在线观看美女被高潮喷水网站| 午夜免费观看性视频| av免费观看日本| 久久99热6这里只有精品| 久久久亚洲精品成人影院| 男人爽女人下面视频在线观看| 国产精品久久久久久久电影| 王馨瑶露胸无遮挡在线观看| www.av在线官网国产| 黄色怎么调成土黄色| 国产精品久久久久久av不卡| 直男gayav资源| 美女主播在线视频| 久热久热在线精品观看| 亚洲精品,欧美精品| 麻豆乱淫一区二区| 国产女主播在线喷水免费视频网站| 亚洲婷婷狠狠爱综合网| 亚洲中文av在线| 人人妻人人添人人爽欧美一区卜 | 国产黄频视频在线观看| 亚洲精品第二区| 成人毛片a级毛片在线播放| 亚洲av免费高清在线观看| 99精国产麻豆久久婷婷| 国产成人a∨麻豆精品| 黄色配什么色好看| 最近中文字幕高清免费大全6| 日韩成人伦理影院| 91久久精品国产一区二区三区| 午夜免费观看性视频| 免费久久久久久久精品成人欧美视频 | 亚洲国产欧美人成| 国产亚洲av片在线观看秒播厂| av国产免费在线观看| 国产精品国产三级国产专区5o| 男女边摸边吃奶| 亚洲国产色片| 国产免费一级a男人的天堂| 日本欧美国产在线视频| 插逼视频在线观看| 最后的刺客免费高清国语| 熟女电影av网| 18禁裸乳无遮挡免费网站照片| 黑丝袜美女国产一区| 久久午夜福利片| 国产黄片美女视频| 欧美成人a在线观看| 亚洲中文av在线| 日日摸夜夜添夜夜爱| 免费看光身美女| 91在线精品国自产拍蜜月| 亚洲欧美日韩东京热| 啦啦啦中文免费视频观看日本| 日韩国内少妇激情av| 久久青草综合色| 亚洲国产日韩一区二区| 高清不卡的av网站| 国产成人精品福利久久| 国产成人免费无遮挡视频| www.av在线官网国产| 一二三四中文在线观看免费高清| 欧美激情极品国产一区二区三区 | 在线亚洲精品国产二区图片欧美 | 亚洲性久久影院| 人人妻人人添人人爽欧美一区卜 | 国产成人91sexporn| 日日啪夜夜爽| 我的女老师完整版在线观看| 国产 一区 欧美 日韩| 日韩中文字幕视频在线看片 | 精品久久久久久久久av| 亚洲av中文字字幕乱码综合| 91久久精品国产一区二区三区| 少妇的逼水好多| 成人18禁高潮啪啪吃奶动态图 | 精品久久国产蜜桃| 男人添女人高潮全过程视频| 又粗又硬又长又爽又黄的视频| 一级毛片aaaaaa免费看小| 国国产精品蜜臀av免费| 亚洲精华国产精华液的使用体验| 国国产精品蜜臀av免费| 18+在线观看网站| 青春草亚洲视频在线观看| 麻豆成人av视频| 在线观看美女被高潮喷水网站| 成人高潮视频无遮挡免费网站| 久久99热这里只频精品6学生| 欧美激情极品国产一区二区三区 | 嘟嘟电影网在线观看| 精品久久久久久久久亚洲| 熟妇人妻不卡中文字幕| 国产午夜精品久久久久久一区二区三区| 色视频www国产| 亚洲天堂av无毛| 99热这里只有精品一区| 草草在线视频免费看| 亚洲精品乱码久久久v下载方式| 99久久精品热视频| 久久久a久久爽久久v久久| 国产男女内射视频| 中国美白少妇内射xxxbb| 在线观看美女被高潮喷水网站| 久久精品国产亚洲av涩爱| 91久久精品国产一区二区成人| 亚洲av二区三区四区| 18禁裸乳无遮挡免费网站照片| 男女免费视频国产| 在线免费十八禁| 2022亚洲国产成人精品| 国产一级毛片在线| 日本wwww免费看| tube8黄色片| 亚洲欧洲日产国产| 尾随美女入室| 高清欧美精品videossex| 免费观看a级毛片全部| h视频一区二区三区| 王馨瑶露胸无遮挡在线观看| 国产精品偷伦视频观看了| 精品久久久噜噜| 少妇猛男粗大的猛烈进出视频| 亚洲精品日韩在线中文字幕| 国产色爽女视频免费观看| 观看美女的网站| 1000部很黄的大片| 久久精品国产a三级三级三级| 欧美性感艳星| 大话2 男鬼变身卡| 高清毛片免费看| 偷拍熟女少妇极品色| 日本-黄色视频高清免费观看| 一级av片app| 国国产精品蜜臀av免费| 国产有黄有色有爽视频| 婷婷色av中文字幕| 男女边吃奶边做爰视频| 一级毛片黄色毛片免费观看视频| 国产精品蜜桃在线观看| a级毛片免费高清观看在线播放| 高清午夜精品一区二区三区| 国产精品伦人一区二区| 97在线视频观看| 肉色欧美久久久久久久蜜桃| 极品教师在线视频| 国产精品欧美亚洲77777| 在线天堂最新版资源| 欧美极品一区二区三区四区| 舔av片在线| 亚洲精品aⅴ在线观看| 伊人久久国产一区二区| 天美传媒精品一区二区| 日韩人妻高清精品专区| tube8黄色片| 美女福利国产在线 | 久久久久网色|