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

    水稻矮桿突變體D814基因圖位克隆與功能分析

    2018-01-03 03:42:41陳新兵黃榮峰
    生物技術(shù)進(jìn)展 2017年6期
    關(guān)鍵詞:信號(hào)轉(zhuǎn)導(dǎo)突變體乙烯

    陳新兵, 黃榮峰, 王 娟

    中國農(nóng)業(yè)科學(xué)院生物技術(shù)研究所, 北京 100081

    水稻矮桿突變體D814基因圖位克隆與功能分析

    陳新兵, 黃榮峰*, 王 娟*

    中國農(nóng)業(yè)科學(xué)院生物技術(shù)研究所, 北京 100081

    水稻株高對(duì)作物產(chǎn)量有著重要的影響,在水稻整個(gè)生長發(fā)育過程中,株高受到多因素的調(diào)控,而植物激素乙烯就是重要的影響因素之一。用10 mg/m3乙烯處理水稻幼苗,對(duì)水稻突變體庫進(jìn)行篩選,獲得了3個(gè)根伸長生長對(duì)乙烯敏感性降低的突變體,其中1個(gè)突變體D814表現(xiàn)出植株矮化、分蘗數(shù)減少、千粒重下降等特征。圖位克隆將其定位在1號(hào)染色體上1 cM的區(qū)間內(nèi),該區(qū)間有6個(gè)已報(bào)道的矮桿突變基因,通過對(duì)這6個(gè)基因測(cè)序,發(fā)現(xiàn)其中1個(gè)基因OsBRI1(LOC_Os01g52050)發(fā)生了點(diǎn)突變(編碼區(qū)第1 837位G突變?yōu)門)。并在D814中分別對(duì)OsBRI1的2個(gè)同源基因(OsBRL1和OsBRL3)進(jìn)行測(cè)序,發(fā)現(xiàn)這2個(gè)基因均無突變。利用已報(bào)道的OsBRI1等位突變體gsor300084進(jìn)行乙烯處理,發(fā)現(xiàn)gsor300084與D814一樣,表現(xiàn)出根對(duì)乙烯敏感性降低。OsBRI1是植物激素油菜素內(nèi)酯(brassinosteroid,BR)的信號(hào)受體,經(jīng)檢測(cè),BR信號(hào)途徑響應(yīng)基因在D814突變體中的表達(dá)也有變化,說明D814是OsBRI1的1個(gè)等位突變體。功能分析發(fā)現(xiàn),D814參與乙烯信號(hào)轉(zhuǎn)導(dǎo)調(diào)控途徑和植物鹽脅迫應(yīng)答途徑。研究結(jié)果為探究乙烯調(diào)控水稻生長發(fā)育及耐逆性的分子機(jī)理提供了研究材料,也為進(jìn)一步探討油菜素內(nèi)酯與乙烯協(xié)同調(diào)控水稻生長發(fā)育機(jī)制奠定了理論基礎(chǔ)。

    水稻;乙烯;矮桿;圖位克??;鹽脅迫應(yīng)答

    水稻矮化是植物生長發(fā)育過程中的一個(gè)重要考察指標(biāo),對(duì)于作物產(chǎn)量、光合速率等具有重要的影響。矮化變異可作為一種穩(wěn)定性狀標(biāo)記,用于新品種培育、作物遺傳改良等生產(chǎn)實(shí)踐[1]。根據(jù)矮化程度不同,可將水稻分為:半矮桿、矮桿、極矮桿3種表型,其中矮桿是指水稻株高表型比正常株高矮一半的矮桿突變體[2]。按照慣例,矮桿和極矮桿基因統(tǒng)一用符號(hào)d來命名,而半矮桿基因用符號(hào)sd來命名。目前在國家水稻數(shù)據(jù)庫網(wǎng)站上(www.ricedate.cn)已鑒定出的矮桿基因有幾十個(gè),其中分布在1號(hào)染色體上的矮桿基因有6個(gè),分別為sd1、d2、d10、d18、d61、gid1。在這6個(gè)基因當(dāng)中,d61是較早被發(fā)現(xiàn)的一個(gè)矮桿基因,是編碼植物激素油菜素內(nèi)酯(brassinosteroid,BR)的一個(gè)受體激酶BRI1。在水稻中,OsBRI1有2個(gè)同源基因,OsBRL1與OsBRL3,OsBRI1主要通過控制細(xì)胞分裂以及細(xì)胞伸長來調(diào)控器官發(fā)育,但對(duì)器官的起始不是必需的,而OsBRL1和OsBRL3主要是在根部響應(yīng)BR信號(hào)的應(yīng)答[3]。除此之外,OsBRI1在水稻的生長發(fā)育過程中起多種作用,包括節(jié)間伸長、葉傾角、暗形態(tài)建成等[4]。近年來通過正向遺傳學(xué)等手段也發(fā)現(xiàn)了許多OsBRI1的等位突變基因,例如半矮桿突變體Fn189(d54)的點(diǎn)突變導(dǎo)致OsBRIl激酶活性降低[5]。另外一個(gè)點(diǎn)突變體gsor300084,其點(diǎn)突變導(dǎo)致第444位的色氨酸突變?yōu)榫彼?,從而?dǎo)致對(duì)BR信號(hào)的鈍感表型[6]。OsBRI1除了參與植物生長發(fā)育以外,在逆境脅迫應(yīng)答中也發(fā)揮著重要的作用。例如:鹽脅迫下,OsBRI1的表達(dá)量會(huì)顯著上升,從而導(dǎo)致植物體內(nèi)BR信號(hào)增強(qiáng);隨著氧化脅迫強(qiáng)度的上升,BR能夠通過促進(jìn)植物體內(nèi)H2O2的積累來使植物抵御這一傷害[7]。

    對(duì)于BR信號(hào)轉(zhuǎn)導(dǎo)機(jī)制,在擬南芥中已經(jīng)研究的十分清楚[8,9]。在水稻中,目前為止只鑒定出了一些與擬南芥相似的BR信號(hào)成分。OsBAK1與擬南芥中BAK1同源[10],OsBZR1與擬南芥中BZR1的功能類似[11]。而擬南芥BIN2的同源基因OsGSK2,在水稻中同樣是BR信號(hào)的負(fù)調(diào)控因子[12]。除此之外,近年來BR信號(hào)轉(zhuǎn)導(dǎo)中一些新的轉(zhuǎn)錄因子也陸續(xù)被鑒定出來:OsOFP8作為OsGSK2的作用底物參與BR信號(hào)途徑[13];OsBUL1作為一個(gè)正調(diào)控轉(zhuǎn)錄因子參與BR信號(hào)轉(zhuǎn)導(dǎo)途徑[14]。隨著生物學(xué)的發(fā)展,水稻中越來越多的BR信號(hào)途徑新組分被鑒定出來,豐富和發(fā)展著BR信號(hào)調(diào)控網(wǎng)絡(luò)。在水稻生長發(fā)育的過程中,乙烯對(duì)植物的株高以及地下部生長發(fā)育具有十分重要的影響。挖掘乙烯信號(hào)轉(zhuǎn)導(dǎo)途徑中的新組分是研究乙烯在植物體內(nèi)作用的分子基礎(chǔ),也是目前常規(guī)的流行方法。作為重要的植物激素之一,乙烯在雙子葉植物體內(nèi)的作用主要表現(xiàn)為典型的三重反應(yīng):抑制根的伸長生長、抑制下胚軸的伸長生長以及促進(jìn)頂端彎曲;而在單子葉植物中乙烯只表現(xiàn)出二重反應(yīng):抑制根的伸長生長和促進(jìn)胚芽鞘的伸長生長[15]。研究表明,乙烯在植物體內(nèi)的這種作用主要是通過乙烯信號(hào)轉(zhuǎn)導(dǎo)途徑實(shí)現(xiàn)的[16~21]。除了調(diào)控植物生長發(fā)育作用外,乙烯還參與植物鹽脅迫、病原菌侵染等逆境脅迫條件下的應(yīng)答過程[22]。其中在鹽脅迫應(yīng)答中,乙烯可以通過調(diào)節(jié)水稻根系Na+/K+的轉(zhuǎn)運(yùn)或者通過下游轉(zhuǎn)錄因子調(diào)控耐鹽相關(guān)基因的表達(dá)[23]。

    本研究用10 mg/m3乙烯處理水稻幼苗,對(duì)水稻突變體庫進(jìn)行篩選,獲得了3個(gè)根對(duì)乙烯敏感性降低突變體,其中1個(gè)突變體D814表現(xiàn)出植株矮化、分蘗數(shù)減少、千粒重下降等顯著特征。圖位克隆將其定位在1號(hào)染色體上1 cM的區(qū)間內(nèi),該區(qū)間有6個(gè)已報(bào)道的矮桿突變基因,通過對(duì)這6個(gè)基因測(cè)序發(fā)現(xiàn)其中1個(gè)基因OsBRI1(LOC_Os01g52050)發(fā)生點(diǎn)突變(第1 837位G突變?yōu)門)。同時(shí)檢測(cè)到D814中OsBRI1的2個(gè)同源基因(OsBRL1和OsBRL3)均無突變,且已報(bào)道的d61等位突變體gsor300084也表現(xiàn)出根對(duì)乙烯敏感性降低,而D814中BR響應(yīng)基因的表達(dá)也有所變化,上述結(jié)果說明該突變體很可能是OsBRI1基因點(diǎn)突變所致。進(jìn)一步的功能分析發(fā)現(xiàn),D814參與了乙烯信號(hào)轉(zhuǎn)導(dǎo)調(diào)控途徑和植物鹽脅迫應(yīng)答途徑。本研究為探究乙烯調(diào)控水稻生長發(fā)育及耐逆性的分子機(jī)理提供了研究材料,也為進(jìn)一步探討B(tài)R與乙烯協(xié)同調(diào)控水稻生長發(fā)育奠定了理論基礎(chǔ)。

    1 材料與方法

    1.1 試驗(yàn)材料與生長條件

    1.1.1試驗(yàn)材料 野生型02428和培矮64(PA64)水稻種子由本實(shí)驗(yàn)室保存,3 000份以02428為背景的甲基磺酸乙酯(ethyl methane sulfonate,EMS)誘變水稻突變體材料庫由本實(shí)驗(yàn)室創(chuàng)建,乙烯敏感性降低矮桿突變體D814由本次實(shí)驗(yàn)篩選并保存。用于D814突變體基因圖位克隆的F2群體由本實(shí)驗(yàn)室創(chuàng)建。

    1.1.2生長條件 2015年和2016年,連續(xù)2年在試驗(yàn)基地(39.6°N,116.2°E)進(jìn)行種植。土培幼苗:種植在土培盒中,在水稻培養(yǎng)溫室中(30℃,16 h白天/8 h黑夜)生長2周左右;水培幼苗:種植在鐵絲網(wǎng)架上,放入溫室培養(yǎng)4~5 d左右。

    1.1.3試驗(yàn)試劑 Trizol試劑盒購自天根生化有限公司;反轉(zhuǎn)錄試劑盒購自諾唯贊有限公司;用于熒光定量PCR的SYBR Green Mix購自翌圣生物科技有限公司。

    1.2 乙烯篩選

    選取萌發(fā)一致的水稻幼苗,種植在鐵絲網(wǎng)架上,每個(gè)鐵網(wǎng)格種植30粒水稻種子,每個(gè)株系種植3格,將網(wǎng)架放入10 L的塑料密閉培養(yǎng)箱中,加入5 L的水,使水面距離種子約1 cm處,封蓋,注入終濃度為10 mg/m3的乙烯氣體,長日照條件下培養(yǎng)4 d,進(jìn)行表型觀察并統(tǒng)計(jì)初生根長,實(shí)驗(yàn)設(shè)3次重復(fù)。

    1.3 表型分析及農(nóng)藝性狀統(tǒng)計(jì)

    在水稻幼苗期,選取土培2周左右的幼苗進(jìn)行拍照;在水稻成熟期,分別對(duì)野生型02428和突變體D814植株進(jìn)行觀察統(tǒng)計(jì):隨機(jī)選取10株成苗,測(cè)量統(tǒng)計(jì)其株高、分蘗數(shù)、各節(jié)間長度等主要農(nóng)藝性狀;等種子成熟后,分別稱取野生型02428和突變體D814千粒重,并進(jìn)行3次重復(fù)。

    1.4 基因初步定位

    將D814與PA64進(jìn)行雜交獲得F1代后加繁一代得到F2代群體。選取飽滿度大致相同的F2群體以及親本野生型02428和PA64水稻種子,萌發(fā)后土培種植,在水稻溫室中培養(yǎng)2周左右,選取矮化表型幼苗單株,根據(jù)CTAB法提取基因組DNA后,首先隨機(jī)選取10個(gè)單株F2基因組DNA,按等濃度混勻構(gòu)成DNA混池,同時(shí)分別以野生型親本02428以及PA64的基因組DNA為模板進(jìn)行PCR鑒定,PCR反應(yīng)程序?yàn)?95℃ 10 min;95℃ 15 s,55℃ 15 s,72℃ 15 s,共45個(gè)循環(huán);72℃ 5 min。PCR反應(yīng)體系為:模板1 μL,正反向引物各0.5 μL,2×PCR Mix 10 μL, 用H2O補(bǔ)齊到20 μL。PCR引物根據(jù)水稻基因組數(shù)據(jù)庫中水稻品種02428和PA64的Indel分子標(biāo)記進(jìn)行設(shè)計(jì),PCR產(chǎn)物用4%(w/V)瓊脂糖凝膠電泳。通過交換率分析確定D814突變基因所在染色體區(qū)域。再根據(jù)國家水稻數(shù)據(jù)庫網(wǎng)站(http:rapdb.dna.affrc.go.jp/)對(duì)該區(qū)域基因的注釋與分析,選取候選基因。最后設(shè)計(jì)引物,將PCR產(chǎn)物測(cè)序后確定突變基因。

    1.5 耐鹽性表型試驗(yàn)

    選取萌發(fā)一致的水稻種子,種植在裝有草炭土的種植盒中,每盒種植20粒水稻種子,每個(gè)株系種植8盒,分為對(duì)照組和處理組,封土,蓋好保鮮膜保濕,放在溫室中培養(yǎng),待種子出苗后揭開保鮮膜,繼續(xù)培養(yǎng)2周左右,選取長勢(shì)一致的水稻幼苗放入150 mmol/L NaCl溶液中生長,每天觀察,待葉片萎蔫后進(jìn)行復(fù)水處理,每天觀察并統(tǒng)計(jì)復(fù)水存活率,實(shí)驗(yàn)設(shè)3次重復(fù)。

    1.6 實(shí)時(shí)熒光定量PCR(Q-PCR)

    分別剪取生長2周左右的野生型02428和突變體D814植株幼苗葉片,按Trizol試劑盒說明書步驟操作提取水稻總RNA,再用反轉(zhuǎn)錄試劑盒反轉(zhuǎn)錄成cDNA。通過熒光定量PCR(IQ5 multicolor Real Time PCR Detection System)檢測(cè)乙烯信號(hào)途徑關(guān)鍵基因的表達(dá)。PCR反應(yīng)程序?yàn)?95℃ 10 min;95℃ 10 s,58℃ 30 s,72℃ 20 s,共45個(gè)循環(huán);以每30 s上升0.5℃的速率從55℃升至95℃來做為溶解曲線程序。PCR反應(yīng)體系為:模板1 μL,正反向引物各0.5 μL,SYBR Green Mix 10 μL, 用H2O補(bǔ)齊到20 μL。所用引物序列見表1。

    2 結(jié)果與分析

    2.1 乙烯不敏感突變體篩選

    通過對(duì)野生型02428進(jìn)行EMS誘變,獲得了3 000份水稻突變體材料,用10 mg/m3乙烯對(duì)其進(jìn)行篩選,得到3個(gè)根伸長生長對(duì)乙烯敏感性降低的突變體,分別是D848、D814以及D841(圖1,彩圖見圖版一)。用已報(bào)道的水稻乙烯不敏感突變體ein2(日本晴背景)作為對(duì)照,在黑暗條件下水培4 d。統(tǒng)計(jì)數(shù)據(jù)發(fā)現(xiàn),在10 mg/m3乙烯處理下,野生型日本晴(Nip)和ein2初生根長分別縮短為原來的35.45%和76.72%。野生型02428的初生根長在10 mg/m3乙烯處理下縮短為正常條件下的31.69%,而D848、D814和D841突變體經(jīng)同樣處理后,初生根長度分別縮短為原來的43.35%、65.50%、41.88%,該結(jié)果表明這3個(gè)突變體的根伸長生長對(duì)乙烯敏感性均有所降低,其中D814根的伸長生長對(duì)乙烯的不敏感性最明顯。

    表1 實(shí)驗(yàn)所用引物序列Table 1 Sequences of the primers used in the experiment.

    圖1 根對(duì)乙烯敏感性降低的水稻突變體篩選Fig.1 The screening of rice mutants with root decreased sensitivity to ethylenein.注:圖中為黑暗下生長4 d的幼苗在乙烯處理下的表型。其中,CK為對(duì)照組,未進(jìn)行乙烯處理;C2H4為處理組,進(jìn)行乙烯處理。(彩圖見圖版一)

    2.2 D814參與調(diào)控乙烯信號(hào)途徑

    為了進(jìn)一步研究D814突變體根對(duì)乙烯敏感性降低的分子機(jī)制,在水稻幼苗期進(jìn)行乙烯處理后,通過Q-PCR分別檢測(cè)了乙烯合成途徑和信號(hào)轉(zhuǎn)導(dǎo)途徑的下游調(diào)控基因的表達(dá)情況。對(duì)于乙烯合成家族基因,當(dāng)外源施加乙烯處理水稻幼苗時(shí),無論是野生型還是突變體,基因表達(dá)量都顯著升高(圖2A),而對(duì)于乙烯信號(hào)轉(zhuǎn)導(dǎo)途徑中的下游調(diào)控基因,當(dāng)外源施加乙烯處理水稻幼苗時(shí),OsSHR5和OsERF063的基因表達(dá)量只在野生型中有所上升,而在突變體中并未有明顯變化(圖2B)。該結(jié)果表明D814突變體阻斷了乙烯信號(hào)轉(zhuǎn)導(dǎo)途徑,而對(duì)乙烯合成途徑無影響,從而說明D814是通過乙烯信號(hào)轉(zhuǎn)導(dǎo)途徑來調(diào)控植物根對(duì)乙烯的敏感性降低。

    2.3 D814突變體表型分析

    圖2 D814調(diào)控乙烯信號(hào)途徑基因的表達(dá)Fig.2 D814 regulates the expression of ethylene signaling-related genes.注:A:乙烯合成途徑中的關(guān)鍵酶編碼基因表達(dá);B:乙烯信號(hào)下游基因OsSHR5和乙烯信號(hào)下游轉(zhuǎn)錄因子OsERF063的基因表達(dá)。CK為對(duì)照組,未進(jìn)行乙烯處理; C2H4為處理組,進(jìn)行乙烯處理。相對(duì)定量以ACTIN2為內(nèi)參基因。**表示與對(duì)照組相比差異極顯著(P<0.01)。

    d61突變體典型的表型就是株高變矮,其成熟期突變體高度相對(duì)野生型高度只有60%~70%左右,除此之外,d61突變體還會(huì)表現(xiàn)出直立葉、葉色深、突變體種子變黑等系列特征,其他的表型特征還包括葉夾角、葉鞘長度、分蘗數(shù)等[3~5]。分別對(duì)D814突變體在不同時(shí)期的表型進(jìn)行了觀察。在幼苗期,D814突變體與02428相比表現(xiàn)出株高顯著變矮(圖3A、3B,彩圖見圖版一);在成熟期,D814突變體與02428相比,在株高、花期、葉型、葉夾角、粒長、粒寬等主要農(nóng)藝性狀方面都存在顯著差異(圖3C~圖3I),尤其是突變體平均株高為野生型的50.44%,OsBRI1的等位突變體d61-2以及gsor300084等都表現(xiàn)出相似的表型。除此之外,D814也像其他OsBRI1等位突變體一樣,表現(xiàn)出典型的直立葉、葉色深黑以及褐粒等。對(duì)D814突變體成熟期的各個(gè)節(jié)間進(jìn)行統(tǒng)計(jì),發(fā)現(xiàn)D814突變體株高變矮的原因主要是第一、第三和第四節(jié)間變短造成的。同時(shí)由于D814突變體粒長顯著小于野生型、粒寬略小于野生型,其千粒重也明顯低于野生型(表2)。這些結(jié)果表明D814能夠調(diào)控水稻的生長發(fā)育,而且D814在株高及種子的發(fā)育過程中具有十分重要的作用。

    圖3 D814突變體與野生型表型比較Fig.3 Comparison of the phenotypes between D814 mutant and wild type.A:幼苗期植株;B:幼苗期根部;C:成熟期植株;D:成熟期節(jié)間長度;E:成熟期旗葉; F:成熟期穗部;G:成熟期穗頸及第一節(jié)間;H:成熟期籽長;I:成熟期粒寬(彩圖見圖版一)

    2.4 D814突變基因定位

    表2 D814突變體與野生型農(nóng)藝性狀比較Table 2 Comparison of agronomic traits between wild type and D814 mutant.

    注:*和**分別表示同組數(shù)據(jù)中突變體與野生型相比差異顯著(P<0.05)和極顯著(P<0.01)。

    利用D814的矮桿表型對(duì)其突變基因進(jìn)行定位,在D814與培矮64(PA64)雜交的F2群體中篩選出了204株具有矮桿性狀的突變體單株。提取這些單株基因組DNA后,利用圖位克隆基因定位技術(shù),首先將突變區(qū)間鎖定在1號(hào)染色體上,初步定位在2個(gè)標(biāo)記(rm1-33和rm1-34)之間,該區(qū)間的遺傳距離為1 cM(圖4A)。通過水稻數(shù)據(jù)庫(http://rice.plant)查閱發(fā)現(xiàn)該區(qū)間已報(bào)道的矮桿基因有6個(gè),分別為:調(diào)控赤霉素信號(hào)途徑中的GA20氧化酶基因sd1,參與獨(dú)角金內(nèi)酯合成途徑的細(xì)胞色素P450蛋白CYP90D2基因d2,類胡蘿卜素裂解雙加氧酶編碼基因d10,編碼赤霉素3β-羥化酶的基因d18,編碼油菜素內(nèi)酯信號(hào)轉(zhuǎn)導(dǎo)途徑中的BR受體激酶的基因d61,以及編碼赤霉素信號(hào)途徑的GA受體的基因gid1。為了確定D814突變體矮化表型是否是由于這些基因引起的,分別提取野生型02428、突變體D814以及F2群體中隨機(jī)2個(gè)矮桿單株基因組DNA進(jìn)行測(cè)序,結(jié)果表明只有d61(OsBRI1,LOC_Os01g52050)在其編碼區(qū)第1 837位堿基發(fā)生了點(diǎn)突變(由G突變?yōu)門)(圖4B),導(dǎo)致其蛋白第613位氨基酸由甘氨酸變成了纈氨酸(G613V)(圖4C)。以上結(jié)果表明,D814突變體矮化表型可能是由于1號(hào)染色體上的d61等位突變引起。

    圖4 突變基因D814在1號(hào)染色體上的定位圖Fig.4 Mapping of mutated gene D814 on chromosome 1.A、B:D814突變基因初步定位圖與測(cè)序結(jié)果;C:d61等位突變

    為了進(jìn)一步驗(yàn)證該突變體是否由OsBRI1基因突變引起的,首先利用野生型以及突變體基因組DNA對(duì)其同源基因OsBRL1與OsBRL3啟動(dòng)子及基因組全長進(jìn)行測(cè)序,發(fā)現(xiàn)只有OsBRI1基因在其CDS區(qū)第1 837位發(fā)生了堿基突變(由G突變?yōu)锳),而其同源基因并沒有發(fā)生突變。另一方面,由于突變體D814是在水稻突變體庫中用乙烯處理篩選出來的,于是利用已經(jīng)報(bào)道的以野生型松前為背景的d61水稻突變體gsor300084,與突變體D814一起進(jìn)行乙烯處理,并觀察統(tǒng)計(jì)水稻根長,結(jié)果發(fā)現(xiàn),突變體gsor300084與D814一樣,初生根都表現(xiàn)出對(duì)乙烯的敏感性降低(圖5,彩圖見圖版一)。野生型02428和松前的初生根長在乙烯處理下分別縮短為正常條件下的26.47%和22.03%,而突變體D814和gsor300084分別縮短為57.14%和64.40%。從而進(jìn)一步表明D814突變體可能是d61的等位突變體。

    圖5 d61等位突變體gsor300084乙烯表型Fig.5 Comparison of the ethylene insensitive phenotypes between D814 and d61 mutant.注:圖中為光下生長4 d的幼苗在乙烯處理下的表型。CK為對(duì)照組,未進(jìn)行乙烯處理; C2H4為處理組,進(jìn)行乙烯處理。(彩圖見圖版一)

    2.5 D814參與BR信號(hào)轉(zhuǎn)導(dǎo)途徑

    為了驗(yàn)證D814參與BR信號(hào)途徑,提取水培3 d的水稻幼苗總RNA,反轉(zhuǎn)錄成cDNA后檢測(cè)BR信號(hào)響應(yīng)基因的表達(dá)情況。結(jié)果發(fā)現(xiàn),OsBRI1及其同源基因OsBRL1與OsBRL3基因表達(dá)量相對(duì)于野生型并沒有明顯變化,BR信號(hào)傳導(dǎo)途徑基因OsBAK1,其表達(dá)量相對(duì)于野生型也沒有明顯變化。而對(duì)于BR信號(hào)下游響應(yīng)基因OsOFP8、OsBLE1和OsBUL1,在突變體D814中其基因表達(dá)水平相對(duì)于野生型02428顯著降低,而BR信號(hào)途徑負(fù)調(diào)控因子OsGSK2的基因表達(dá)略有升高(圖6)。該結(jié)果表明D814參與BR信號(hào)轉(zhuǎn)導(dǎo)途徑,也進(jìn)一步證明了D814突變體很可能就是d61等位突變體。

    2.6 D814參與植物鹽脅迫應(yīng)答

    圖6 D814調(diào)控BR信號(hào)途徑基因的表達(dá)Fig.6 D814 regulates the expression of BR signaling-related genes.注:相對(duì)定量以ACTIN2為內(nèi)參基因。*表示與野生型02428相比差異顯著(P<0.05);**表示與野生型02428相比差異極顯著(P<0.01)。

    鹽脅迫是植物非生物脅迫的一個(gè)重要方面,植物在抵御鹽脅迫過程中,一方面通過Na+/K+轉(zhuǎn)移等離子通道,另一方面通過調(diào)控體內(nèi)一些植物激素水平來減輕鹽脅迫對(duì)植物造成的傷害。之前的研究[24]表明,乙烯作為一種重要的植物激素,能夠負(fù)調(diào)控水稻耐鹽性,因此比較了突變體D814和野生型02428的耐鹽性。觀察發(fā)現(xiàn),使用150 mmol/L NaCl處理7 d后,水稻幼苗開始出現(xiàn)葉片卷曲,其中突變體發(fā)生卷曲的葉片數(shù)少于野生型(圖7A,彩圖見圖版二)。拍照統(tǒng)計(jì)后對(duì)處理組幼苗進(jìn)行復(fù)水,復(fù)水5 d后拍照并統(tǒng)計(jì)成活率。結(jié)果顯示,野生型水稻幼苗復(fù)水后成活率為9.14%,而突變體復(fù)水后成活率為56.08%(圖7B)。說明D814突變體對(duì)高鹽脅迫的耐受力明顯高于野生型,并進(jìn)一步表明D814能參與調(diào)控水稻鹽脅迫應(yīng)答。

    3 討論

    隨著對(duì)植物激素乙烯的深入研究,越來越多的乙烯信號(hào)轉(zhuǎn)導(dǎo)途徑中的新組分被挖掘出來。例如在水稻中篩選到的乙烯信號(hào)轉(zhuǎn)導(dǎo)途徑突變體mhz6等[25],這些新組分的篩選不但闡明和完善了乙烯在植物體內(nèi)的信號(hào)調(diào)控網(wǎng)絡(luò),而且為育種學(xué)家培育新品種提供了理論基礎(chǔ)。

    圖7 D814突變體與野生型鹽脅迫表型比較Fig.7 Comparison of mutants D814 and wild type in response to salt stress.注:A:野生型和突變體幼苗鹽處理表型。B:野生型和突變體幼苗鹽處理后復(fù)水成活率。**表示同組數(shù)據(jù)與野生型相比極顯著(P<0.01)。(彩圖見圖版二)

    本研究首先用10 mg/m3濃度乙烯處理2周左右的水稻幼苗,對(duì)實(shí)驗(yàn)室突變體庫進(jìn)行篩選,獲得了3個(gè)根對(duì)乙烯敏感性降低的突變體:D841、D814和D848,從而為后續(xù)探究乙烯在植物體內(nèi)的作用及分子機(jī)理提供了寶貴的材料。由于突變體D814初生根的伸長生長對(duì)乙烯的敏感性最低,因此選擇D814突變體作為進(jìn)一步的研究材料。表型觀察發(fā)現(xiàn)D814在整個(gè)植株生長發(fā)育的過程中,都表現(xiàn)出植株矮化、分蘗數(shù)減少、千粒重下降等顯著特征,與大多數(shù)OsBRI1基因的等位突變體相似,具有直立葉、葉色深等表型,說明D814可能是BR信號(hào)途徑有關(guān)的突變體,能夠參與水稻生長發(fā)育的調(diào)控,也說明D814在調(diào)控水稻生長發(fā)育過程中具有十分重要的作用,特別是在調(diào)控水稻株高方面。通過圖位克隆,將D814突變基因初步定位在1號(hào)染色體上遺傳距離為1 cM的區(qū)間內(nèi),并通過基因測(cè)序表明其中的d61基因在編碼區(qū)1 837位堿基發(fā)生了點(diǎn)突變(由G變成T),導(dǎo)致氨基酸的改變。從而說明D814突變基因很可能是d61的等位突變。之前的研究報(bào)道了d61的2個(gè)等位突變體d61-1和d61-2也都是點(diǎn)突變體(圖4)。d61編碼一個(gè)BR受體蛋白激酶,突變后會(huì)導(dǎo)致體內(nèi)BR信號(hào)轉(zhuǎn)導(dǎo)受阻,從而表現(xiàn)出節(jié)間伸長被抑制、葉片直立等典型特征[26]。在水稻中,BR不僅能夠單獨(dú)調(diào)控植物的生長發(fā)育,BR也能參與植物激素乙烯的合成與調(diào)控。例如BR與乙烯相互拮抗來共同調(diào)節(jié)下胚軸的生長[27],也可以通過增強(qiáng)乙烯合成限速酶ACS5蛋白的穩(wěn)定性來調(diào)控乙烯的生物合成[28]。BR處理后,可以提高乙烯不敏感突變體ein2在鹽脅迫下種子的發(fā)芽率[29]。而D814作為BR信號(hào)轉(zhuǎn)導(dǎo)受阻的突變體,幼苗期根部表現(xiàn)出對(duì)乙烯敏感性降低的特征,表明D814可作為一個(gè)良好的材料,對(duì)其深入研究將會(huì)為闡明油菜素內(nèi)酯與乙烯協(xié)同調(diào)控水稻生長發(fā)育機(jī)制奠定理論基礎(chǔ)。

    由于初定位的遺傳距離1 cM的區(qū)間范圍太大,在該區(qū)間內(nèi)是否還有其他基因發(fā)生突變尚不清楚,同時(shí)為了進(jìn)一步的驗(yàn)證D814突變基因確實(shí)是OsBRI1,首先通過對(duì)同源基因OsBRL1和OsBRL3啟動(dòng)子和基因組序列全長進(jìn)行測(cè)序發(fā)現(xiàn),其同源基因并沒有發(fā)生任何突變,只有OsBRI1基因在其CDs區(qū)發(fā)生了點(diǎn)突變(圖4),進(jìn)一步推測(cè),如果突變基因確實(shí)是OsBRI1,那么D814突變體就應(yīng)該和已報(bào)道的d61等其他突變體一樣具有相似的表型,于是對(duì)已報(bào)道的d61等位點(diǎn)突變體gsor300084進(jìn)行乙烯處理,觀察D814與gsor300084突變體根部表型,通過統(tǒng)計(jì)發(fā)現(xiàn)D814和gsor300084一樣,根部都表現(xiàn)出對(duì)乙烯敏感性降低(圖5),從而進(jìn)一步說明D814很可能是OsBRI1的一個(gè)新的等位突變體。最后,進(jìn)一步檢測(cè)了D814突變體中BR信號(hào)響應(yīng)基因的表達(dá)情況,發(fā)現(xiàn)突變體中這些基因的表達(dá)量與02428野生型有顯著差異(圖6),進(jìn)一步說明D814參與BR信號(hào)轉(zhuǎn)導(dǎo)途徑,但還需要進(jìn)一步的精細(xì)定位以及遺傳數(shù)據(jù)等來證明。

    乙烯功能的發(fā)揮最終都是將乙烯信號(hào)轉(zhuǎn)導(dǎo)給下游一些轉(zhuǎn)錄因子,從而啟動(dòng)乙烯相關(guān)基因的表達(dá)。在水稻中,乙烯處理能顯著誘導(dǎo)乙烯信號(hào)下游基因OsSHR5(LOC_Os08g10310)與OsERF063(LOC_Os09g11480)超高量表達(dá)[30]。而在D814突變體中,基因OsSHR5與OsERF063的表達(dá)量在乙烯處理前后并沒有顯著的差異,表明D814參與乙烯信號(hào)轉(zhuǎn)導(dǎo),從而使突變體D814表現(xiàn)出根對(duì)乙烯的敏感性下降。之前研究表明,在水稻中乙烯能夠負(fù)調(diào)控水稻的耐鹽性[31],并且發(fā)現(xiàn)D814表現(xiàn)出耐鹽性增加,進(jìn)一步表明D814通過阻斷乙烯信號(hào)途徑提高了植物的耐鹽性。以上研究對(duì)闡明乙烯在調(diào)節(jié)植物生長發(fā)育和鹽脅迫應(yīng)答中發(fā)揮功能的分子機(jī)制具有十分重要的意義。

    綜上所述,本研究鑒定了3個(gè)乙烯敏感性降低的水稻突變體,這些突變體的根都呈現(xiàn)了對(duì)乙烯的鈍感表型。通過圖位克隆方法初步鑒定D814可能為BR途徑的受體激酶OsBRI1點(diǎn)突變。進(jìn)一步研究發(fā)現(xiàn),D814突變體呈現(xiàn)耐鹽性增加,揭示了乙烯可能通過BR途徑調(diào)控根生長和鹽脅迫反應(yīng)的機(jī)制。這些研究對(duì)于揭示水稻乙烯信號(hào)轉(zhuǎn)導(dǎo)機(jī)制具有重要意義,為探究乙烯調(diào)控水稻生長發(fā)育及耐逆性的機(jī)理提供了研究材料,也為進(jìn)一步探討油菜素內(nèi)酯與乙烯協(xié)同調(diào)控水稻生長發(fā)育機(jī)制奠定了理論基礎(chǔ)。

    [1] 閔紹楷, 申宗坦, 熊振民. 水稻育種學(xué)[M]. 北京: 中國農(nóng)業(yè)科技出版杜, 1997, 12(6): 1-6.

    [2] 馬漢云, 王青林, 吳淑平, 等. 水稻株型遺傳的研究進(jìn)展[J]. 安徽農(nóng)學(xué)通報(bào), 2007, 13(23): 105-107.

    [3] Nakamura A, Fujioka S, Sunohara H,etal.. The role ofOsBRI1 and its homologous genes,OsBRL1 andOsBRL3, in rice[J]. Plant Physiol., 2006, 140(2): 580-590.

    [4] Yamamuro C, Ihara Y, Wu X,etal.. Loss of function of a ricebrassinosteroidinsensitive1 homolog prevents internode elongation and bending of the lamina joint[J]. Plant Cell, 2000, 12(9): 1591-1605.

    [5] Zhao J, Wu C, Yuan S,etal.. Kinase activity of OsBRI1 is essential for brassinosteroids to regulate rice growth and development[J]. Plant Sci., 2013, 199: 113-120.

    [6] Gao Y, Wang G, Yuan S,etal.. Phenotypic analysis and molecular characterization of an allelic mutant of thed61 gene in rice [J]. Crop J., 2014, 2(4): 175-182.

    [7] Xia X J, Wang Y J, Zhou Y H,etal.. Reactive oxygen species are involved in brassinosteroid-induced stress tolerance in cucumber[J]. Plant Physiol., 2009, 150(2): 801-814.

    [8] Wang X, Chory J. Brassinosteroids regulate dissociation of BKI1, a negative regulator of BRI1 signaling, from the plasma membrane[J]. Science, 2006, 313(5790): 1118-1122.

    [9] He J X, Gendron J M, Yang Y,etal.. The GSK3-like kinase BIN2 phosphorylates and destabilizes BZR1, a positive regulator of the brassinosteroid signaling pathway inArabidopsis[J]. Proc. Natl. Acad. Sci. USA, 2002, 99(15): 10185-10190.

    [10] Li D, Wang L, Wang M,etal.. EngineeringOsBAK1 gene as a molecular tool to improve rice architecture for high yield[J]. Plant Biotechnol. J., 2009, 7(8): 791-806.

    [11] Bai M Y, Zhang L Y, Gampala S S,etal.. Functions of OsBZR1 and 14-3-3 proteins in brassinosteroid signaling in rice[J]. Proc. Natl. Acad. Sci. USA, 2007, 104(34): 13839-13844.

    [12] Tong H, Liu L, Jin Y,etal.. Dwaftandlow-tilleringacts as a direct downstream target of a GSK3/SHAGGY-like kinase to mediate brassinosteroid responses in rice[J]. Plant Cell, 2012, 24(6): 2562-2577.

    [13] Yang C, Shen W, He Y,etal.. OVATE family protein 8 positively mediates brassinosteroidsignaling through interacting with the GSK3-like kinase in rice[J]. PLoS Genet., 2016, 12(6): e1006118.

    [14] Tanaka A, Nakagawa H, Tomita C,etal..Brassinosteroidupregulated1, encoding a helix-loop-helix protein, is a novel gene involved in brassinosteroid signaling and controls bending of the lamina joint in rice[J]. Plant Physiol., 2009, 151(2): 669-680.

    [15] Yang C, Lu X, Ma B,etal.. Ethylene signaling in rice andArabidopsis: Conserved and diverged aspects [J]. Mol. Plant, 2015, 8(4): 495-505.

    [16] Ecker J R. The ethylene signal transduction pathway in plants[J]. Science, 1995, 268(5211): 667-675.

    [17] Kieber J J, Rothenberg M, Roman G,etal..CTR1, a negative regulator of the ethylene response pathway inArabidopsis, encodes a member of the raf family of protein kinases [J]. Cell, 1993, 72(3):427-441.

    [18] Ju C, Yoon G M, Shemansky J M,etal.. CTR1 phosphorylates the central regulator EIN2 to control ethylene hormone signaling from the ER membrane to the nucleus inArabidopsis[J]. Proc. Natl. Acad. Sci. USA, 2012, 109(47): 19486-19491.

    [19] Alonso J M, Hirayama T, Roman G,etal.. EIN2, a bifunctional transducer of ethylene and stress responses inArabidopsis[J]. Science, 1999, 284(5423): 2148-2152.

    [20] Wen X, Zhang C, Ji Y,etal.. Activation of ethylene signaling is mediated by nuclear translocation of the cleaved EIN2 carboxyl terminus[J]. Cell Res., 2012, 22(11): 1613-1616.

    [21] Qiao H, Chang K, Yazaki J,etal.. Interplay between ethylene, ETP1/ETP2 F-box proteins, and degradation of EIN2 triggers ethylene responses inArabidopsis[J]. Gene. Dev., 2009, 23(4): 512-521.

    [22] Solano R, Stepanova A, Chao Q,etal.. Nuclear events in ethylene signaling: A transcriptional cascade mediated by ethylene-insensitive3 and ethylene-response-factor1[J]. Gene. Dev., 1998, 12(23): 3703-3714.

    [23] Cao W, Liu J, He X J,etal.. Modulation of ethylene responses affects plant salt-stress responses [J]. Plant Physiol., 2007, 143(2): 707-719.

    [24] Yang C, Ma B, He S,etal.. Maohuzi6/ethylene insensitive3-like1 and ethylene insensitive3-like2 regulate ethylene response of roots and coleoptiles and negatively affect salt tolerance in rice[J]. Plant Physiol., 2015, 169(1): 148-165.

    [25] Wu D, Ji J, Wang G,etal..LchERF, a novel ethylene-responsive transcription factor fromLyciumchinense, confers salt tolerance in transgenic tobacco[J]. Plant Cell Rep., 2014, 33(12): 2033-2045.

    [26] Bishop G, Yokota T. Plants steroid hormones, brassinosteroids: Current highlights of molecular aspects on their synthesis/metabolism, transport, perception and response[J]. Plant Cell Physiol., 2001, 42(2): 114-120.

    [27] Mori M, Nomura T, Ooka H,etal.. Isolation and characterization of a rice dwarf mutant with a defect in brassinosteroid biosynthesis[J]. Plant Physiol., 2002, 130(3): 1152-1161.

    [28] Deslauriers S, Larsen P. FERONIA is a key modulator of brassinosteroid and ethylene responsiveness inArabidopsishypocotyls[J]. Mol. Plant, 2010, 3(3): 626-640.

    [29] Hansen M, Chae H, Kieber J. Regulation of ACS protein stability by cytokinin and brassinosteroid[J]. Plant J., 2009, 57(4): 606-614.

    [30] Divi U, Krishna P. Overexpression of the brassinosteroid biosynthetic geneAtDWF4 inArabidopsisseeds overcomes abscisic acid-induced inhibition of germination and increases cold tolerance in transgenic seedlings[J]. J. Plant Growth Regul., 2010, 29(4): 385-393.

    [31] Ma B, He S, Duan K,etal.. Identification of rice ethylene-response mutants and characterization of MHZ7/OsEIN2 in distinct ethylene response and yield trait regulation[J]. Mol. Plant, 2013, 6(6): 1830-1848.

    Map-basedCloningandFunctionAnalysisofRiceDwarfMutantD814

    CHEN Xinbing, HUANG Rongfeng*, WANG Juan*

    BiotechnologyResearchInstitute,ChineseAcademyofAgriculturalSciences,Beijing100081,China

    Plant height confers important effects on crop yield in rice, which is regulated by many aspects in the growth and development process. Ethylene is one of important factors for regulating plant growth and stress response. In this research, we identified a mutantD814 from the mutants insensitive to ethylene in roots growth.D814 mutants showed the dwarf phenotype, decreased effective tilling numbers and reduced 1000-grain weight. Map-based cloning indicated thatD814 located in 1 cM range on chromosome 1, in which there contained six genes contributed to dwarf phenotype in the previous researches. Sequencing analyses demonstrated thatOsBRI1(LOC_Os01g52050) has a point mutation at 1 837 bp downstream of ATG, in which the G is changed to T. There were two homologous genes ofOsBRI1 in rice,OsBRL1 andOsBRL3, of which the sequences have no mutation inD814 mutants. A reported mutant ofOsBRI1,gsor300084, also showed decreased sensitivity to ethylene in root. Due to the function of OsBRI1 was a receptor in brassinosteroid (BR) signal pathway, we found the expressions of BR-response genes were affected inD814, which indicated thatD814 was an allele mutant ofOsBRI1. Moreover, our analyses showed thatD814 affected the expression levels of ethylene signaling-related genes and displayed increased salt tolerance. Thus our research provided a novel research material for exploring the molecular mechanism of ethylene regulating rice growth and stress tolerance, and laid a foundation for investigating the interplay of ethylene and brassinosteroid in growth and development.

    rice; ethylene; dwarf; map-based cloning; salt stress response

    2017-05-10;接受日期2017-07-10

    國家自然科學(xué)基金項(xiàng)目(31470366)資助。

    陳新兵,碩士研究生,研究方向?yàn)樗拘禄虬l(fā)掘與功能研究。E-mail:chenxinbing_caas@163.com。*通信作者:黃榮峰,研究員,主要從事植物激素與逆境脅迫應(yīng)答研究。E-mail:rfhuang@caas.cn;王 娟,助理研究員,主要從事作物蛋白質(zhì)功能研究。E-mail:wangjuan@caas.cn

    10.19586/j.2095-2341.2017.0040

    猜你喜歡
    信號(hào)轉(zhuǎn)導(dǎo)突變體乙烯
    Wnt/β-catenin信號(hào)轉(zhuǎn)導(dǎo)通路在瘢痕疙瘩形成中的作用機(jī)制研究
    乙烷裂解制乙烯產(chǎn)業(yè)大熱
    CLIC1及其點(diǎn)突變體與Sedlin蛋白的共定位研究
    擬南芥干旱敏感突變體篩選及其干旱脅迫響應(yīng)機(jī)制探究
    HGF/c—Met信號(hào)轉(zhuǎn)導(dǎo)通路在結(jié)直腸癌肝轉(zhuǎn)移中的作用
    Survivin D53A突變體對(duì)宮頸癌細(xì)胞增殖和凋亡的影響
    乙烯裂解爐先進(jìn)控制系統(tǒng)開發(fā)及工業(yè)應(yīng)用
    兩個(gè)基于二噻吩乙烯結(jié)構(gòu)單元雙核釕乙烯配合物的合成,表征和性質(zhì)
    2015年乙烯進(jìn)口關(guān)稅下調(diào)至1%補(bǔ)充國內(nèi)需求缺口
    中國氯堿(2014年12期)2014-02-28 01:05:10
    鈣敏感受體及其與MAPK信號(hào)轉(zhuǎn)導(dǎo)通路的關(guān)系
    三级毛片av免费| 中文字幕av电影在线播放| av网站免费在线观看视频| 999久久久国产精品视频| 亚洲国产中文字幕在线视频| 黄片大片在线免费观看| 大香蕉久久网| 天天躁狠狠躁夜夜躁狠狠躁| 成年人午夜在线观看视频| 制服人妻中文乱码| 亚洲美女黄色视频免费看| 日本黄色日本黄色录像| 99久久人妻综合| 亚洲欧美精品综合一区二区三区| 欧美日韩成人在线一区二区| 国产深夜福利视频在线观看| 日韩电影二区| 视频区欧美日本亚洲| a级毛片在线看网站| 天堂中文最新版在线下载| 久久国产精品大桥未久av| 又黄又粗又硬又大视频| 中亚洲国语对白在线视频| 欧美97在线视频| 999精品在线视频| 美女国产高潮福利片在线看| 亚洲av男天堂| 日韩欧美一区二区三区在线观看 | 丝瓜视频免费看黄片| 高潮久久久久久久久久久不卡| 人成视频在线观看免费观看| 美女视频免费永久观看网站| 午夜激情久久久久久久| 久久性视频一级片| a级毛片在线看网站| 欧美激情高清一区二区三区| 亚洲久久久国产精品| 亚洲国产毛片av蜜桃av| 女人被躁到高潮嗷嗷叫费观| 亚洲专区中文字幕在线| 亚洲欧美日韩另类电影网站| 黄色视频在线播放观看不卡| 精品欧美一区二区三区在线| 午夜老司机福利片| www.精华液| 婷婷丁香在线五月| 国产精品秋霞免费鲁丝片| 这个男人来自地球电影免费观看| 久久久久国内视频| 一区二区三区精品91| 国产精品九九99| 狠狠精品人妻久久久久久综合| 50天的宝宝边吃奶边哭怎么回事| 亚洲精品久久午夜乱码| 免费一级毛片在线播放高清视频 | 国产区一区二久久| 国产亚洲精品久久久久5区| 十八禁高潮呻吟视频| 九色亚洲精品在线播放| 日韩中文字幕欧美一区二区| 9色porny在线观看| 久久热在线av| 男女高潮啪啪啪动态图| 男女下面插进去视频免费观看| 久9热在线精品视频| 91国产中文字幕| 大型av网站在线播放| 精品一品国产午夜福利视频| 不卡一级毛片| 中文字幕制服av| 999久久久精品免费观看国产| 久久人人97超碰香蕉20202| 国产精品麻豆人妻色哟哟久久| 欧美日韩亚洲高清精品| 一本综合久久免费| 午夜福利视频精品| 久久国产亚洲av麻豆专区| 桃花免费在线播放| www.精华液| √禁漫天堂资源中文www| 狂野欧美激情性xxxx| av片东京热男人的天堂| 日本猛色少妇xxxxx猛交久久| cao死你这个sao货| 波多野结衣av一区二区av| 精品久久蜜臀av无| 18禁裸乳无遮挡动漫免费视频| 久久av网站| 国产日韩欧美视频二区| 国产亚洲av片在线观看秒播厂| 欧美黑人欧美精品刺激| 久久免费观看电影| 人妻久久中文字幕网| av电影中文网址| 国产精品一二三区在线看| 国产日韩欧美亚洲二区| 19禁男女啪啪无遮挡网站| 亚洲黑人精品在线| 精品久久蜜臀av无| 王馨瑶露胸无遮挡在线观看| 欧美 日韩 精品 国产| 精品亚洲成a人片在线观看| 91麻豆精品激情在线观看国产 | 日韩,欧美,国产一区二区三区| 老鸭窝网址在线观看| 日本av免费视频播放| 欧美日韩亚洲国产一区二区在线观看 | 亚洲午夜精品一区,二区,三区| 亚洲精品国产av蜜桃| 午夜福利视频在线观看免费| 免费在线观看完整版高清| 日韩三级视频一区二区三区| 男女下面插进去视频免费观看| 一区二区av电影网| 18禁裸乳无遮挡动漫免费视频| 午夜成年电影在线免费观看| 伦理电影免费视频| 男男h啪啪无遮挡| 99久久精品国产亚洲精品| 欧美在线一区亚洲| 美女高潮喷水抽搐中文字幕| 日韩欧美国产一区二区入口| 夜夜骑夜夜射夜夜干| 国产亚洲av高清不卡| 免费一级毛片在线播放高清视频 | 2018国产大陆天天弄谢| 搡老熟女国产l中国老女人| 少妇粗大呻吟视频| 亚洲男人天堂网一区| 一级,二级,三级黄色视频| 精品欧美一区二区三区在线| 99久久综合免费| 成人av一区二区三区在线看 | 黄色a级毛片大全视频| 九色亚洲精品在线播放| 国产成人欧美| 久久精品熟女亚洲av麻豆精品| 欧美激情极品国产一区二区三区| 老鸭窝网址在线观看| 日本猛色少妇xxxxx猛交久久| 咕卡用的链子| 久久久久久人人人人人| 国产精品国产av在线观看| 亚洲精品中文字幕一二三四区 | 三级毛片av免费| 亚洲精品av麻豆狂野| 动漫黄色视频在线观看| 日韩欧美一区视频在线观看| 伊人亚洲综合成人网| 国产成人av教育| 国产一区二区 视频在线| www日本在线高清视频| 飞空精品影院首页| 日韩欧美一区二区三区在线观看 | 亚洲五月婷婷丁香| 中文字幕av电影在线播放| 少妇精品久久久久久久| 一区在线观看完整版| 中文字幕高清在线视频| 飞空精品影院首页| 黄片小视频在线播放| 啪啪无遮挡十八禁网站| 日本vs欧美在线观看视频| 精品国产国语对白av| 少妇粗大呻吟视频| 亚洲国产日韩一区二区| 久久九九热精品免费| 国产精品欧美亚洲77777| 色播在线永久视频| 婷婷成人精品国产| 1024香蕉在线观看| 69av精品久久久久久 | 亚洲欧洲日产国产| 丁香六月欧美| 国产一级毛片在线| 在线天堂中文资源库| 黄片大片在线免费观看| 欧美另类一区| 又紧又爽又黄一区二区| 午夜免费成人在线视频| videosex国产| 99热国产这里只有精品6| 青青草视频在线视频观看| 亚洲精品成人av观看孕妇| 久久国产精品男人的天堂亚洲| 久久久国产精品麻豆| 欧美日韩一级在线毛片| 国产99久久九九免费精品| 一区二区三区激情视频| 美女高潮喷水抽搐中文字幕| 最近最新免费中文字幕在线| 熟女少妇亚洲综合色aaa.| 国产在视频线精品| 亚洲国产成人一精品久久久| 日韩,欧美,国产一区二区三区| 日本av免费视频播放| 精品第一国产精品| 天堂俺去俺来也www色官网| 亚洲国产精品999| 久久人妻熟女aⅴ| 欧美黄色淫秽网站| 亚洲av片天天在线观看| 三上悠亚av全集在线观看| 久久久国产精品麻豆| 男人操女人黄网站| 亚洲情色 制服丝袜| 欧美成人午夜精品| 精品一品国产午夜福利视频| 久久久久久久久久久久大奶| 欧美中文综合在线视频| 亚洲少妇的诱惑av| 久久香蕉激情| 9热在线视频观看99| 国产三级黄色录像| 国产精品一区二区免费欧美 | 大片电影免费在线观看免费| 97在线人人人人妻| 男人操女人黄网站| 国产亚洲av高清不卡| 国产精品自产拍在线观看55亚洲 | 国产男人的电影天堂91| 欧美+亚洲+日韩+国产| 亚洲一码二码三码区别大吗| 欧美乱码精品一区二区三区| 啦啦啦 在线观看视频| 色综合欧美亚洲国产小说| 国产无遮挡羞羞视频在线观看| 91av网站免费观看| 亚洲中文日韩欧美视频| av视频免费观看在线观看| 极品人妻少妇av视频| a 毛片基地| 在线精品无人区一区二区三| 国产片内射在线| 十八禁高潮呻吟视频| 黄片大片在线免费观看| 水蜜桃什么品种好| 黑丝袜美女国产一区| 国产片内射在线| 韩国高清视频一区二区三区| 精品一品国产午夜福利视频| 久久人妻福利社区极品人妻图片| 老司机午夜十八禁免费视频| 女人爽到高潮嗷嗷叫在线视频| 亚洲成人免费av在线播放| 欧美激情高清一区二区三区| 欧美另类亚洲清纯唯美| 午夜福利视频在线观看免费| 十八禁人妻一区二区| 高潮久久久久久久久久久不卡| 超碰97精品在线观看| 男女床上黄色一级片免费看| 亚洲五月婷婷丁香| 夫妻午夜视频| 国产精品一二三区在线看| 深夜精品福利| 久久性视频一级片| 久久久久久久精品精品| 黑人巨大精品欧美一区二区蜜桃| 久久99一区二区三区| 搡老乐熟女国产| 国产精品麻豆人妻色哟哟久久| 别揉我奶头~嗯~啊~动态视频 | 汤姆久久久久久久影院中文字幕| 手机成人av网站| 国产一区二区激情短视频 | 欧美日韩中文字幕国产精品一区二区三区 | 91大片在线观看| 亚洲va日本ⅴa欧美va伊人久久 | a级毛片在线看网站| 淫妇啪啪啪对白视频 | 中文字幕另类日韩欧美亚洲嫩草| 欧美日韩亚洲国产一区二区在线观看 | 一区福利在线观看| 精品久久蜜臀av无| 日本91视频免费播放| 欧美人与性动交α欧美精品济南到| 九色亚洲精品在线播放| 99国产精品一区二区蜜桃av | 高清在线国产一区| av欧美777| 一进一出抽搐动态| 美女主播在线视频| av不卡在线播放| 日韩三级视频一区二区三区| 亚洲欧美日韩高清在线视频 | 亚洲一码二码三码区别大吗| 两性午夜刺激爽爽歪歪视频在线观看 | 又黄又粗又硬又大视频| 91字幕亚洲| 中文字幕人妻丝袜制服| 高清在线国产一区| 免费在线观看影片大全网站| 亚洲av电影在线进入| 久久久国产一区二区| 国产精品国产三级国产专区5o| 人妻 亚洲 视频| 久久精品亚洲av国产电影网| 精品卡一卡二卡四卡免费| 日本欧美视频一区| 精品国产乱码久久久久久小说| 五月天丁香电影| 精品少妇黑人巨大在线播放| 久久精品久久久久久噜噜老黄| 女人爽到高潮嗷嗷叫在线视频| 久久精品aⅴ一区二区三区四区| 如日韩欧美国产精品一区二区三区| 亚洲国产欧美网| 久9热在线精品视频| 大陆偷拍与自拍| 少妇粗大呻吟视频| 黄色视频,在线免费观看| 日本vs欧美在线观看视频| 午夜福利免费观看在线| 亚洲国产欧美在线一区| 黄色 视频免费看| 中文字幕人妻丝袜制服| 天堂8中文在线网| 免费少妇av软件| 日本一区二区免费在线视频| 黄色 视频免费看| 国产精品久久久久久人妻精品电影 | 日本黄色日本黄色录像| 制服人妻中文乱码| 精品一区二区三卡| 色婷婷av一区二区三区视频| 老熟妇仑乱视频hdxx| 国产伦人伦偷精品视频| 国产av国产精品国产| 91麻豆精品激情在线观看国产 | 国产一卡二卡三卡精品| 国产男人的电影天堂91| 精品国产一区二区三区四区第35| 爱豆传媒免费全集在线观看| 91成人精品电影| 午夜老司机福利片| 国产高清视频在线播放一区 | www.自偷自拍.com| 99香蕉大伊视频| 亚洲熟女毛片儿| 欧美精品人与动牲交sv欧美| 成年人午夜在线观看视频| 国产精品.久久久| 色视频在线一区二区三区| 国产精品久久久人人做人人爽| 99热网站在线观看| 黑人巨大精品欧美一区二区蜜桃| 国产欧美日韩综合在线一区二区| 午夜视频精品福利| 亚洲一卡2卡3卡4卡5卡精品中文| 午夜激情av网站| 亚洲精品美女久久av网站| 久久影院123| 一区二区日韩欧美中文字幕| 高清av免费在线| 国产亚洲精品第一综合不卡| 日韩一卡2卡3卡4卡2021年| 亚洲国产欧美日韩在线播放| 麻豆国产av国片精品| 亚洲五月色婷婷综合| 美女大奶头黄色视频| 热99国产精品久久久久久7| 在线永久观看黄色视频| 男女国产视频网站| 久久天堂一区二区三区四区| 法律面前人人平等表现在哪些方面 | 免费观看人在逋| 91精品国产国语对白视频| 大香蕉久久网| 91精品伊人久久大香线蕉| 国产一区二区三区综合在线观看| 脱女人内裤的视频| 啦啦啦 在线观看视频| 50天的宝宝边吃奶边哭怎么回事| 男女床上黄色一级片免费看| 窝窝影院91人妻| 精品少妇黑人巨大在线播放| 午夜免费鲁丝| 在线天堂中文资源库| 午夜免费鲁丝| 五月开心婷婷网| 黄色毛片三级朝国网站| 美女高潮到喷水免费观看| 色视频在线一区二区三区| www.999成人在线观看| 精品一区二区三区av网在线观看 | 啦啦啦 在线观看视频| 一区二区三区四区激情视频| 日韩三级视频一区二区三区| 两个人看的免费小视频| 亚洲精品一二三| 正在播放国产对白刺激| 亚洲第一青青草原| 青春草亚洲视频在线观看| 久久亚洲国产成人精品v| 搡老乐熟女国产| 80岁老熟妇乱子伦牲交| 黑人猛操日本美女一级片| 久久精品熟女亚洲av麻豆精品| 可以免费在线观看a视频的电影网站| 91成人精品电影| 爱豆传媒免费全集在线观看| 美女中出高潮动态图| av片东京热男人的天堂| 男女午夜视频在线观看| 两人在一起打扑克的视频| 亚洲少妇的诱惑av| 99久久国产精品久久久| 伊人久久大香线蕉亚洲五| 日本一区二区免费在线视频| 丁香六月天网| 婷婷丁香在线五月| avwww免费| 中文精品一卡2卡3卡4更新| 日韩视频在线欧美| 亚洲欧洲日产国产| 久久中文看片网| av天堂久久9| 亚洲国产精品999| 黄色a级毛片大全视频| 亚洲精品粉嫩美女一区| 亚洲国产欧美网| 黑丝袜美女国产一区| 91九色精品人成在线观看| 99精国产麻豆久久婷婷| 超碰成人久久| 男女边摸边吃奶| 国产亚洲av高清不卡| 在线观看免费日韩欧美大片| 每晚都被弄得嗷嗷叫到高潮| 亚洲国产av新网站| 亚洲国产毛片av蜜桃av| 色综合欧美亚洲国产小说| 午夜福利一区二区在线看| 中文字幕最新亚洲高清| 啦啦啦 在线观看视频| 色播在线永久视频| 久久中文字幕一级| 水蜜桃什么品种好| 在线观看免费日韩欧美大片| 亚洲 国产 在线| 国产一区二区在线观看av| 波多野结衣av一区二区av| av线在线观看网站| 麻豆国产av国片精品| 久久精品国产a三级三级三级| 99热网站在线观看| 久久亚洲精品不卡| 青春草亚洲视频在线观看| 成年人免费黄色播放视频| 在线观看免费高清a一片| 大陆偷拍与自拍| 欧美激情久久久久久爽电影 | 夫妻午夜视频| 亚洲国产精品999| 啦啦啦视频在线资源免费观看| 国产成人精品在线电影| 不卡av一区二区三区| 欧美精品人与动牲交sv欧美| 亚洲精品中文字幕在线视频| 两人在一起打扑克的视频| 美女中出高潮动态图| 青春草亚洲视频在线观看| 国产一区二区激情短视频 | 国产精品 国内视频| 亚洲欧美激情在线| 欧美日韩福利视频一区二区| 亚洲成人免费av在线播放| 国产亚洲欧美精品永久| 欧美老熟妇乱子伦牲交| 黄色毛片三级朝国网站| 满18在线观看网站| 亚洲成人免费电影在线观看| 亚洲综合色网址| 日韩欧美一区视频在线观看| 人人妻人人添人人爽欧美一区卜| 人人妻人人澡人人爽人人夜夜| 夜夜骑夜夜射夜夜干| 纵有疾风起免费观看全集完整版| 建设人人有责人人尽责人人享有的| 女警被强在线播放| 男人添女人高潮全过程视频| 国产一区二区 视频在线| 亚洲精品国产精品久久久不卡| 天天躁日日躁夜夜躁夜夜| 大型av网站在线播放| 麻豆乱淫一区二区| bbb黄色大片| tube8黄色片| 成人影院久久| 国产亚洲欧美在线一区二区| 99国产精品一区二区三区| 一级,二级,三级黄色视频| 亚洲精品国产区一区二| av片东京热男人的天堂| 他把我摸到了高潮在线观看 | 亚洲国产成人一精品久久久| 日韩大片免费观看网站| 欧美激情 高清一区二区三区| 日日爽夜夜爽网站| 狠狠婷婷综合久久久久久88av| 国产精品一区二区精品视频观看| 波多野结衣一区麻豆| 午夜免费成人在线视频| 亚洲欧洲日产国产| 亚洲中文av在线| 女性被躁到高潮视频| 免费高清在线观看视频在线观看| www.自偷自拍.com| 国产精品欧美亚洲77777| 亚洲国产日韩一区二区| 午夜日韩欧美国产| 亚洲黑人精品在线| 老司机靠b影院| 桃花免费在线播放| 欧美激情 高清一区二区三区| 国产av国产精品国产| 男女免费视频国产| 精品国产超薄肉色丝袜足j| 精品一区二区三区四区五区乱码| 国产淫语在线视频| 黄色a级毛片大全视频| 丝袜美腿诱惑在线| 精品少妇一区二区三区视频日本电影| 亚洲国产av影院在线观看| 国产精品香港三级国产av潘金莲| 黄色怎么调成土黄色| 肉色欧美久久久久久久蜜桃| 日韩一区二区三区影片| 99久久国产精品久久久| 考比视频在线观看| 91精品伊人久久大香线蕉| 欧美精品人与动牲交sv欧美| 久久人人97超碰香蕉20202| 自拍欧美九色日韩亚洲蝌蚪91| 亚洲av美国av| 亚洲自偷自拍图片 自拍| 中文欧美无线码| tube8黄色片| 色视频在线一区二区三区| 51午夜福利影视在线观看| 韩国高清视频一区二区三区| 亚洲精品乱久久久久久| 亚洲,欧美精品.| 99九九在线精品视频| 一级,二级,三级黄色视频| 99久久人妻综合| 国产人伦9x9x在线观看| 国产男女超爽视频在线观看| 久久久久国产精品人妻一区二区| 欧美日韩福利视频一区二区| 日韩欧美一区视频在线观看| 精品卡一卡二卡四卡免费| 这个男人来自地球电影免费观看| 男女午夜视频在线观看| 9色porny在线观看| 香蕉丝袜av| 亚洲熟女毛片儿| 一本—道久久a久久精品蜜桃钙片| 国产91精品成人一区二区三区 | 亚洲av成人不卡在线观看播放网 | 19禁男女啪啪无遮挡网站| 免费日韩欧美在线观看| 麻豆国产av国片精品| avwww免费| 99久久综合免费| 男女边摸边吃奶| 国产日韩欧美亚洲二区| 亚洲国产精品一区三区| 久久国产精品大桥未久av| 精品少妇一区二区三区视频日本电影| 国产成人免费观看mmmm| 欧美激情 高清一区二区三区| 我要看黄色一级片免费的| 国产一区二区三区av在线| 亚洲欧美色中文字幕在线| 免费在线观看日本一区| 99国产精品一区二区蜜桃av | 建设人人有责人人尽责人人享有的| 亚洲成国产人片在线观看| 欧美黑人欧美精品刺激| 99国产精品一区二区三区| 精品熟女少妇八av免费久了| 在线永久观看黄色视频| 国产成+人综合+亚洲专区| 欧美精品一区二区免费开放| 国产在线观看jvid| 脱女人内裤的视频| 久久免费观看电影| 99久久人妻综合| 9色porny在线观看| 99国产精品一区二区蜜桃av | 欧美亚洲日本最大视频资源| 美女扒开内裤让男人捅视频| 久久久国产精品麻豆| 精品少妇一区二区三区视频日本电影| 成人国语在线视频| 十八禁网站免费在线| av免费在线观看网站| 少妇人妻久久综合中文| 国产伦理片在线播放av一区| 男人爽女人下面视频在线观看| 国产精品秋霞免费鲁丝片| kizo精华| 欧美人与性动交α欧美软件| 麻豆国产av国片精品| 国产精品久久久人人做人人爽| 高清欧美精品videossex| 国产欧美日韩综合在线一区二区| 亚洲精品一二三| 女性被躁到高潮视频| 日韩制服骚丝袜av| 国产伦理片在线播放av一区| 亚洲国产精品成人久久小说| 久久久精品国产亚洲av高清涩受|