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

    栽培小麥Brock中抗白粉病相關基因TaRPP13-1B的克隆及功能分析

    2023-02-10 06:05:26劉曉穎王雪晴楊晨曉王光鈺卞云迪王振英
    作物學報 2023年2期
    關鍵詞:白粉白粉病抗病

    劉曉穎 張 馳 王雪晴 楊晨曉 王光鈺 卞云迪 方 芳 王 穎 王振英,*

    栽培小麥Brock中抗白粉病相關基因的克隆及功能分析

    劉曉穎1張 馳1王雪晴1楊晨曉1王光鈺1卞云迪1方 芳2王 穎2王振英1,*

    1天津師范大學生命科學學院 / 天津市動植物抗性重點實驗室, 天津 300387;2天津市寶坻區(qū)林業(yè)發(fā)展服務中心, 天津 301899

    白粉病是小麥生產(chǎn)中的主要病害之一, 發(fā)掘抗病基因并實現(xiàn)抗病基因轉(zhuǎn)育是提高作物抗病性的最經(jīng)濟有效的方法。本研究克隆了一個位于小麥染色體1B上、具有典型CC、NBS和LRR結(jié)構(gòu)域的基因。接種白粉菌后,基因在抗病小麥Brock和BJ-1中表達量雖然出現(xiàn)上下調(diào)波動, 但平均表達水平一直高于感病小麥品種京411。采用病毒誘導的基因沉默和轉(zhuǎn)基因過表達技術進行功能分析, 發(fā)現(xiàn)抑制目標基因表達, 抗病小麥品種Brock對白粉菌的抗性顯著降低; 過表達的轉(zhuǎn)基因小麥津強5號對白粉菌抗性明顯提高。說明基因參與小麥抗白粉病的防御反應過程。該研究為小麥抗病品種的選育提供了有價值的備選基因。

    小麥; NBS-LRR類基因;; 抗白粉病

    目前小麥中已克隆的抗白粉病基因中多數(shù)是NBS-LRR (nucleotide binding site leucine rich repeat, NBS-LRR)類基因, 其中又以CNL (coiled-coil NBS-LRR, CNL)類基因居多[1], 如、、、、、、及其等位基因, 攜帶這些抗性基因的小麥品種田間抗性良好, CNL類基因在小麥白粉病防治上發(fā)揮了重要作用[2-6]。(recognition of, RPP)首先在擬南芥抗霜霉病反應中被鑒定到, 具有典型的CNL結(jié)構(gòu)域, 它通過識別效應蛋白ATR13, 從而誘發(fā)抗病反應[7]。研究發(fā)現(xiàn),基因具有多態(tài)性, 其LRR結(jié)構(gòu)域的多態(tài)性影響了不同生態(tài)型擬南芥對不同病原小種的特異性識別[8-9]。Rentel等[10]、Ramachandran等[11]發(fā)現(xiàn), 霜霉病效應蛋白ATR13能夠抑制ATR13/RPP13復合體誘發(fā)超敏反應。目前, 小麥中已經(jīng)鑒定到多個同源基因, 這些基因在基因序列上都存在較大差異, 且基因功能涉及小麥生長的多個方面。例如,同源基因在小麥抗白粉病[12-14]、小麥抗銹病[15]和小麥面粉顏色方面均發(fā)揮作用[16]。在其他植物的研究中發(fā)現(xiàn),同源基因在柑橘抗白色念珠菌感染[17]、姜黃抗莖腐病[18]、番茄抗黃葉卷曲病毒[19]、葡萄抗白粉病[20]、花生抗黃曲霉病[21]和百合的雜種株高優(yōu)勢[22]也發(fā)揮作用?;蛟诘挚拐婢⒓毦⒉《厩秩竞椭参锷L發(fā)育等多個方面發(fā)揮重要作用。

    前期工作中, 課題組配制了京411的抗病近等基因系BJ-1, 對白粉菌侵染早期(12 hpi)感病小麥京411及其抗病近等基因系BJ-1應答基因轉(zhuǎn)錄組比較分析中發(fā)現(xiàn), 位于1B染色體的在小麥抗病近等基因系BJ-1中表達量上調(diào), NCBI/Blast比對分析發(fā)現(xiàn)該基因具有典型的NBS- LRR類基因結(jié)構(gòu), 為基因家族的同源基因[13,24], 命名為。本研究克隆了這個潛在的抗性基因, 并對其在不同抗、感小麥中的序列組成進行了比較; 分析了在白粉菌脅迫下的表達模式, 通過降低和過表達基因表達水平, 探究該基因在小麥抵抗白粉菌侵染過程中的貢獻, 為小麥抗病品種的選育提供了新的備選基因。

    1 材料與方法

    1.1 試驗材料

    Brock是一個抗白粉病小麥品種, 從英國引進, 在苗期和成株期均高抗小麥白粉病。京411是我國北方種植的農(nóng)藝性狀優(yōu)良的小麥品種, 但易感白粉病。本實驗室于2012年, 以Brock為抗病基因供體, 京411為輪回親本, 通過連續(xù)回交6次, 創(chuàng)制了京411抗病近等基因系BJ-1。BJ-1的產(chǎn)量、生長周期、株高、分蘗等特性與輪回親本京411一致, 它的抗病等級為1級。而抗病基因供體親本Brock為0級。BJ-1農(nóng)藝性狀上偏京411, 抗白粉病特性偏Brock[23]。抗白粉病小麥品種Brock, 感白粉病小麥品種京411和抗病近等基因系BJ-1用于小麥全長cDNA的克隆及白粉菌誘導表達分析。感病春麥品種津強5號用于過表達載體的轉(zhuǎn)化。白粉菌菌株E09 (f. sp.E09,E09)用于苗期抗病性鑒定, 由中國農(nóng)業(yè)科學院植物保護研究所提供。

    1.2 TaRPP13-1B全長cDNA的克隆

    根據(jù)的序列信息設計引物, 用Primer premier 5.0設計引物TaRPP13-1B-F和TaRPP13-1B-R。將感病小麥品種京411、抗病近等基因系BJ-1和抗病小麥品種Brock培養(yǎng)至第1片葉子完全展開后, 采用抖拂法將白粉菌(E09)接種于小麥葉片12 h后, 取葉片用RNAiso Plus (TaKaRa,中國大連)分別提取總RNA, M-MLV反轉(zhuǎn)錄試劑盒(寶生物, 中國大連)反轉(zhuǎn)錄成cDNA。分別以3個小麥品種(系)的cDNA為模板, 擴增目的基因。PCR反應體系25 μL, cDNA 1 μL、上下游引物各0.5 μL、5× GXL緩沖液5 μL、dNTP混合物2 μL、PrimeSTAR GXL DNA聚合酶0.5 μL、雙蒸水10.5 μL。反應程序為: 94℃ 5 min; 98℃ 10 s, 55℃ 15 s, 68℃ 3 min, 35個循環(huán); 68℃ 7 min。PCR產(chǎn)物經(jīng)1%瓊脂糖凝膠電泳檢測回收。按照DNA A-Tailing試劑盒方法加A尾(寶生物, 中國大連), 連接pGEM-T Easy載體(Promega, 美國)并轉(zhuǎn)化大腸桿菌DH5α, 挑取LB Amp+(氨芐青霉素, ampicilline)篩選培養(yǎng)基上的陽性克隆, 經(jīng)PCR和酶切鑒定后送蘇州金唯智生物科技有限公司測序。

    3'和5'RACE片段的擴增按SMART RACE cDNA擴增試劑盒(Clontech, 美國)的說明書進行。根據(jù)已經(jīng)掌握的部分序列, 用Primer 5.0設計特異引物RPP13-1B-GSP1F和RPP13-1B- GSP2F分別進行3'RACE反應的第1和第2輪PCR; 設計特異引物RPP13-1B-5RACE-1R和RPP13-1B- 5RACE-2R分別進行5'RACE反應的第1和第2輪PCR。將3'和5'RACE第2輪PCR擴增的片段分別克隆至pGEM-T Easy載體, 經(jīng)測序驗證后, 將3'RACE、5'RACE及部分序列進行拼接, 得到拼接后的全長cDNA序列。根據(jù)拼接的全長cDNA序列設計引物-ORF-F/R擴增的全長cDNA, 將擴增產(chǎn)物連接至pGEM-T載體并經(jīng)測序驗證。本研究所用引物序列見表1。

    表1 本研究的引物及序列

    將基因序列提交NCBI GenBank進行Blast比對;利用DNAMAN軟件對序列進行同源性分析, 構(gòu)建系統(tǒng)發(fā)育樹。利用cDNA序列的ORF在小麥聯(lián)盟網(wǎng)站(http://202.194.139.32)進行比對, 分析其染色體定位, 預測DNA序列組成。

    1.3 TaRPP13-1B基因的表達分析

    將小麥京411、Brock和BJ-1幼苗培養(yǎng)至一葉一心期時, 抖拂法接種E09 2、4、8、12、24和48 h后分別提取葉片總RNA, 反轉(zhuǎn)錄成cDNA。以SYBR Green I作為熒光染料, 利用熒光定量PCR法(quantitative real time PCR, qRT-PCR) (ABI 7500 FAST, 美國)分析白粉菌侵染后基因在苗期京411、Brock和BJ-1中的表達水平。用2–DDCT法[25]計算白粉菌侵染過程中目標基因的相對表達量變化。具體步驟按照FastStart Universal SYBR Green Master (Roche, 美國)試劑說明書進行, PCR體系20 μL, 2× Master Mix (SYBR Green I) 10 μL、上下游引物各1 μL、cDNA 1 μL。將上述反應體系混合均勻, 進行熒光定量PCR擴增。擴增程序: 50℃ 2 min; 95℃ 10 min; 95℃ 15 s, 58℃ 30 s, 40個循環(huán)。以作為內(nèi)參, 每個反應重復5次, 實驗重復3次。

    1.4 利用病毒誘導的基因沉默(virus-induced gene silencing, VIGS)技術分析TaRPP13-1B的功能

    根據(jù)的測序結(jié)果, 在基因的非保守區(qū)設計引物RPP13-1B-V-F/R, 擴增帶有I接頭的基因片段。BSMVγ:載體經(jīng)I酶切、去磷酸化后, 與沉默基因片段連接, PCR驗證后獲得基因沉默載體BSMVγ:。BSMVα和BSMVγ經(jīng)I酶切, BSMVβ經(jīng)I酶切后, 體外轉(zhuǎn)錄成BSMV病毒, 將病毒摩擦接種于抗病小麥Brock葉片表面沉默目標基因。以GKP buffer組、BSMV:組和BSMV:組為對照組。待BSMV:對照組小麥第三葉白化, 選取少量小麥葉片提取總RNA, 利用qRT-PCR法分析目標基因的沉默效率, 具體步驟同1.3。在已接種病毒的Brock葉片表面接種白粉菌, 接種白粉菌48 h、72 h和7 d后, 利用考馬斯亮藍染色, 在顯微鏡下觀察白粉菌孢子萌發(fā)狀況[26], 照相記錄接菌7 d后葉片表面的白粉菌生長狀態(tài)。明確基因功能。每組統(tǒng)計5片葉子, 每片葉子統(tǒng)計100個白粉菌孢子, 實驗重復3次。

    1.5 利用農(nóng)桿菌介導的遺傳轉(zhuǎn)化法進行TaRPP13- 1B功能鑒定

    以連接有序列的T載體為模板, 擴增獲得帶有E II和II酶切位點和接頭序列的的全長, 利用無縫克隆試劑(全式金, 中國北京)將連接至pCAMBIA1301的相應酶切位點。將連接產(chǎn)物轉(zhuǎn)至大腸桿菌DH5α, PCR驗證后, 最終獲得包含全長, 啟動子為35S的過表達載體35S:。利用本實驗室改良的農(nóng)桿菌介導的小麥成熟胚遺傳轉(zhuǎn)化法[27-28], 將含有35S:的農(nóng)桿菌C58C1轉(zhuǎn)入感病春麥津強5號成熟胚, 收獲T0代種子。利用潮霉素篩選、PCR和RT-PCR鑒定轉(zhuǎn)基因小麥株系。以轉(zhuǎn)化pCAMBIA1301質(zhì)粒的津強5號為對照組。利用潮霉素篩選35S:株系至T1代。將T1代轉(zhuǎn)基因小麥培養(yǎng)至一葉一心期, 在小麥葉片接種白粉菌E09, 逐日觀察葉片表型變化直至第7天, 拍照記錄。逐株調(diào)查材料的侵染類型, 按照病葉分級標準分級記錄, 標準如下。1級: 病斑占葉面積的10%以下; 2級: 病斑占葉面積的10%~25%; 3級: 病斑占葉面積的26%~50%; 4級: 病斑占葉面積的51%~ 80%; 5級: 病斑占葉面積的80%以上[29-30]。

    2 結(jié)果與分析

    2.1 TaRPP13-1B基因的克隆

    根據(jù)序列設計一對特異引物TaRPP13-1B-F/R, 在京411、Brock和BJ-1中擴增得到基因的cDNA序列。為確認該序列是否為全長, 進一步利用RACE技術擴增了該基因的5'和3'非編碼區(qū)。最后克隆得到基因的cDNA全長3232 bp, 包含2934 bp的ORF以及104 bp和194 bp的5'和3'非編碼區(qū)。感病小麥品種京411和抗病近等基因系BJ-1中的基因序列與Brock相同。在小麥聯(lián)盟網(wǎng)站分別進行Blast分析, 推測該基因位于1B染色體; DNA序列長3453 bp, 由2個外顯子(1~933 bp和1453~3453 bp)和1個內(nèi)含子組成(934~1452 bp) (圖1-A)。進一步的氨基酸比對分析發(fā)現(xiàn), 該基因編碼蛋白具有CC、NB-ARC和LRR結(jié)構(gòu)域, 也是一個典型的CNL類型蛋白(圖1-B)。蛋白質(zhì)的同源比對及進化樹分析發(fā)現(xiàn), 其氨基酸序列與二穗短柄草(XP_010237386.1), 烏拉爾圖小麥(EMS68464.1)和玉米(XP_00864940.1)中RPP13的同源性較高, 分別為93%、79%和66% (圖1-C)。根據(jù)染色體定位信息, 將其命名為。

    圖1 TaRPP13-1B基因的結(jié)構(gòu)域分析

    A:的基因結(jié)構(gòu); B: TaRPP13-1B蛋白的結(jié)構(gòu)域; C: 不同植物中RPP13蛋白的進化樹。

    A: predictive schematic forsequence; B: predictive domain schematic for coding TaRPP13-1B protein; C: the phylogenetic tree of diverse RPP13-like proteins family.

    將與的基因序列在IWGSC和小麥聯(lián)盟網(wǎng)站進行Blast比對發(fā)現(xiàn),和基因分別定位于染色體1B和7D, 這與最初的參考基因序列預期一致。這兩個基因均與基因家族具有較高同源性, 具有典型CNL結(jié)構(gòu)域, 但二者的氨基酸序列一致性僅有23.2% (圖2), 由此我們認為這兩個基因在小麥抵抗白粉菌侵染過程中的作用有待進一步深入研究。

    2.2 TaRPP13-1B基因的表達分析

    在抗病小麥品種Brock中, 除白粉菌侵染4 hpi,基因表達短暫下調(diào),基因表達水平一直高于抗病近等基因系BJ-1和感病小麥品種京411; BJ-1中的本底表達量與京411持平, 2 hpi時表達量逐漸上調(diào), 4 hpi達到高峰, 之后逐漸下降; 京411中基因表達模式與BJ-1一致, 但整體表達量低于抗病近等基因系BJ-1。雖然表達量遠低于BJ-1, 但表達高峰同樣出現(xiàn)在4 hpi; 之后, 表達量開始下降??傊? 抗病材料中雖然部分接菌時間段內(nèi)(4 hpi)基因的表達出現(xiàn)下調(diào), 但基因的表達水平一直高于感病品種(圖3)。故推測基因參與小麥對白粉菌的抗性反應。

    圖2 TaRPP13-1B和TaRPP13-3的同源性分析

    2.3 瞬時沉默TaRPP13-1B基因能夠顯著降低小麥對白粉菌的抗性

    2.3.1基因沉默小麥的表型分析 由于基因在抗病小麥品種Brock中表達水平較高, 在VIGS實驗中以Brock為沉默對象。待Brock的第二片葉展開, 接種重組病毒。繼續(xù)培養(yǎng)至第三片葉展開, 進行表型觀察, 結(jié)果如圖4-A所示。BSMV:對照組中小麥Brock葉片白化, 說明沉默體系有效; 利用熒光定量qRT-PCR檢測各組材料中的基因的表達量, 從圖4-C可以看出, BSMV:實驗組中目標基因表達量低于另外兩個對照組, 說明瞬時沉默基因表達成功。接種白粉菌7 d后, 觀察BSMV:實驗組、BSMV:對照組以及GKP-buffer對照組的葉片表面白粉菌生長狀態(tài)(圖4-B)。接種白粉菌7 d后, BSMV:實驗組植株的葉片上已經(jīng)出現(xiàn)大量白色菌斑, 2個對照組葉片上也出現(xiàn)菌斑, 但數(shù)量少于實驗組。這些結(jié)果說明降低抗病小麥Brock中基因表達, 能夠增加白粉菌成功侵染效率, 小麥葉片的抗病力下降。

    2.3.2基因沉默小麥葉片孢子生長情況分析 進一步觀察接種白粉菌孢子不同時間點的孢子萌發(fā)狀況。接種白粉菌孢子48 h后, 對照組葉片上的白粉菌孢子多數(shù)無法侵入, 為畸形附著胞(圖5-A, B), 而BSMV:實驗組葉片上則開始出現(xiàn)附著胞(圖5-C); 接種白粉菌72 h后, 對照組中可以發(fā)現(xiàn)少量附著胞, 但仍以纖細型或分瓣型的畸形附著胞為主(圖5-D, E), BSMV:實驗組葉片上開始出現(xiàn)少量次級菌絲(圖5-F), 實驗組葉片上白粉菌孢子生長速度明顯快于對照組, 在接菌早期就侵染成功; 接種白粉菌7 d后, 對照組和實驗組葉片上均發(fā)育出分生孢子, 但實驗組小麥葉片表面的孢子數(shù)量更多(圖5-G~I)。這些結(jié)果說明降低抗病小麥Brock中基因表達后, 小麥葉片上白粉菌孢子萌發(fā)速度加快, 分生孢子數(shù)量增多, 小麥的抗病力下降。

    圖3 TaRPP13-1B基因的白粉菌誘導表達分析

    Jing 411: 感病小麥品種京411; BJ-1: 京411抗病近等基因系; Brock: 抗病小麥品種; Student’s測驗,*:< 0.05;**:< 0.01。

    Jing 411: wheat susceptible variety; BJ-1: near isogenic line; Brock: wheat resistant variety; Student’s-test,*:< 0.05;**:< 0.01.

    2.4 農(nóng)桿菌轉(zhuǎn)化法鑒定TaRPP13-1B基因的功能

    2.4.1 感病小麥津強5號的遺傳轉(zhuǎn)化及轉(zhuǎn)基因植株的分子鑒定 利用農(nóng)桿菌介導的遺傳轉(zhuǎn)化法將重組載體轉(zhuǎn)入津強5號小麥成熟胚, 并對該基因的T0代植株進行轉(zhuǎn)基因鑒定, 獲得11株T0代植株(圖6)。培養(yǎng)轉(zhuǎn)基因苗至T1代, 利用qRT-PCR檢測目標基因表達, 發(fā)現(xiàn)部分轉(zhuǎn)基因小麥中的目標基因表達量顯著高于對照組, 平均為對照組的3.74倍, 最高的達到對照組的5倍(T1-1~T1-4號轉(zhuǎn)基因植株), 最低約為對照組的2.6倍(T1-2和T1-3號轉(zhuǎn)基因植株) (圖7)。說明成功轉(zhuǎn)入津強5號小麥并穩(wěn)定遺傳,實現(xiàn)目標基因在受體中表達。

    2.4.2T1代轉(zhuǎn)基因植株的抗病性鑒定 為了進一步驗證該基因的抗白粉病特性, 分別選取了4個過表達的T1代轉(zhuǎn)基因株系(T1-4、T1-6、T1-8和T1-11)培養(yǎng)至一葉一心期, 利用E09侵染轉(zhuǎn)基因小麥及其對照組, 觀察葉片表型。接種白粉菌48 h和72 h時, 野生型和轉(zhuǎn)基因小麥葉片表面的白粉菌孢子的生長情況相近, 無顯著差別。白粉菌侵染7 d后, 野生型和轉(zhuǎn)基因小麥葉片表面都出現(xiàn)白色菌落, 但轉(zhuǎn)基因小麥葉片表面菌落少于野生型(圖8-A)。為了進一步確定轉(zhuǎn)基因小麥對白粉病的抗病程度, 利用ImageJ對轉(zhuǎn)基因小麥及其對照組葉片病斑面積進行統(tǒng)計分析, 結(jié)果如圖8-B所示, 對照組葉片病斑面積為64.6%, 發(fā)病等級為4級(感病)。4個轉(zhuǎn)基因小麥葉片病斑面積較小, 在10%~25%之間, 發(fā)病等級均為2級。以上結(jié)果說明, 過表達提高了感病小麥津強5號的抗病性,基因是抗白粉病相關基因。

    圖4 沉默TaRPP13-1B后小麥抗病性變化

    A: Brock接種BSMV病毒后表型; B: 沉默葉片接種白粉菌E09的表型觀察; C: 各實驗組中基因的相對表達量(*:< 0.05)。

    A: phenotype of Brock leaves inoculated with BSMV; B:-knockdown plants inoculated withE09; C: the relative expression of(*:< 0.05).

    (圖5)

    LA: 分瓣型畸形附著胞; SA: 纖細型附著胞; AGT: 喙型附著胞; SH: 次生菌絲; BC: 念珠狀分生孢子。

    LA: lobed appressoria; SA: slender appressoria; AGT: appressorium germ tube; SH: the secondary hypha; BC: beaded conidia

    圖6 TaRPP13-1B T0代轉(zhuǎn)基因小麥的PCR檢測

    M: DL2000 marker; 1: 1301對照組; 2: 質(zhì)粒陽性對照; 3~14: T0代轉(zhuǎn)基因植株。

    M: DL2000 marker; 1: control; 2: positive control; 3–14: T0transgenic wheat.

    圖7 TaRPP13-1B在T1代轉(zhuǎn)基因小麥中的表達量分析

    C: 對照組; 1~11: T1-1~T1-11 轉(zhuǎn)基因單株(**:< 0.01)。

    C: control; 1–11: T1-1–T1-11: transgenic wheat (**:< 0.01).

    3 討論

    植物中克隆了大量的CNL類抗病基因。是Yahiaoui等[31]利用圖位克隆法克隆得到第1個抗白粉病基因, 包含典型的CNL結(jié)構(gòu)域。Brunner等[32-33]將及其等位基因轉(zhuǎn)化感病六倍體小麥, 所有轉(zhuǎn)基因小麥均表現(xiàn)出田間抗白粉菌特性。也是一個典型的CNL抗病基因。目前, 通過國家和省級品種審定攜帶基因小麥新品種有30多個,基因是中國小麥白粉病抗性育種利用的重要抗源[3,34-35]。[36-39]、[40]和(及其等位基因)[41-43]也都是CNL抗白粉病基因, 攜帶這些抗性基因的小麥品種均表現(xiàn)出良好的田間抗性。這些已有研究結(jié)果說明CNL類基因在小麥白粉病防治上發(fā)揮了重要作用。

    前期研究中, 我們對感病小麥京411及其抗病近等基因系BJ-1的染菌前后的轉(zhuǎn)錄組分析中發(fā)現(xiàn)了位于染色體1B的同源基因,該基因在白粉菌侵染后表達量上調(diào), 且在抗病材料上調(diào)明顯[24]。本研究以為基礎, 結(jié)合RACE技術, 克隆了Brock、BJ-1和京411中基因cDNA序列, 3個小麥品種中的cDNA序列相同, 它與二穗短柄草(XP_010237386.1)有較高同源性(93%), 與烏拉爾圖小麥(EMS68464.1)有較高同源性(79%)。雖然不同抗性小麥品種中的基因cDNA序列相同, 但目標基因在抗/感小麥品種中的表達量差異顯著, 抗病小麥中表達水平高于感病小麥。所以我們推測可能是不同小麥品種中的啟動子區(qū)或上游調(diào)控基因存在差異, 但仍需進一步研究證明。與基因一樣[13], 它們都屬于基因家族。在IWGSC進行序列比對分析后發(fā)現(xiàn), 二者位于不同的染色體, 序列一致性僅有23.2%。根據(jù)二者基因序列很難預測基因功能之間的聯(lián)系。

    圖8 過表達TaRPP13-1B基因提高小麥對白粉菌E09的抗性

    A: 轉(zhuǎn)基因小麥在白粉菌E09侵染7 d后的葉片表型; B: 白粉菌E09侵染下的小麥葉片病斑面積(**:< 0.01)。

    A: phenotype of transgenic wheat leaves inoculated withE09 for 7 days; B: the lesion area in transgenic wheat leaves (**:< 0.01).

    利用VIGS技術沉默基因, 白粉菌侵染實驗表明, 接菌48 h和72 h后, 實驗組葉片上白粉菌生長速度明顯快于對照組, 說明目標基因表達量降低有利于病原菌侵染, 目標基因在病原菌侵染早期發(fā)揮作用。接菌7 d后的實驗組和對照組葉片上都出現(xiàn)肉眼可見的白粉菌孢子, 但對照組葉片上孢子堆較少, 這表明基因表達被抑制后在一定程度上降低了抗病小麥植株的抗性。進一步構(gòu)建過表達基因的轉(zhuǎn)基因小麥, 接種白粉菌7 d, 發(fā)現(xiàn)轉(zhuǎn)基因小麥葉片表面白粉菌生長數(shù)量和狀態(tài)明顯慢于野生型津強5號。這些研究結(jié)果進一步確認,基因參與小麥-白粉菌互作過程, 發(fā)揮正向調(diào)節(jié)作用。

    4 結(jié)論

    本研究利用RACE技術克隆了小麥。沉默/過表達目標基因后的功能研究表明, 該基因參與小麥對白粉菌侵染的早期抗病防御反應, 在小麥-白粉菌互作過程中發(fā)揮正向作用, 研究結(jié)果為小麥白粉病抗性品種的選育提供備選基因。

    [1] Kourelis J, Van Der Hoorn R A L. Defended to the nines: 25 years of resistance gene cloning identifies nine mechanisms for R protein function., 2018, 30: 285–299.

    [2] Zou S H, Wang H, Li Y W, Kong Z S, Tang D Z. The NB-LRR geneconfers powdery mildew resistance in wheat., 2018, 218: 298–309.

    [3] Xing L P, Hu P, Liu J Q, Witek K, Zhou S, Xu J F, Li W H, Gao L, Huang Z P, Zhang R Q, Wang X E, Chen P D, Wang H Y, Jones D G, Karafiátová M, Vrána J, Barto? J, Dole?e l J, Tian Y C, Wu Y F, Cao A Z.fromencodes a CC-NBS-LRR protein conferring powdery mildew resistance in wheat., 2018, 11: 874–878.

    [4] Li M M, Dong L L, Li B B, Wang Z Z, Xie J Z, Qiu D, Li Y H, Shi W Q, Yang L J, Wu Q H, Chen Y X, Lu P, Guo G H, Zhang H Z, Zhang P P, Zhu K Y, Li Y W, Zhang Y, Wang R G, Yuan C G, Liu W, Yu D Z, Luo M C, Fahima T, Nevo E, Li H J, Liu Z Y. A CNL protein in wild emmer wheat confers powdery mildew resistance., 2020, 228: 1027–1037.

    [5] Xie J Z, Guo G H, Wang Y, Hu T Z, Wang L L, Li J T, Qiu D, Li Y H, Wu Q H, Lu P, Chen Y X, Dong L L, Li M M, Zhang H Z, Zhang P P, Zhu K Y, Li B B, Deal K R, Huo N X, Zhang Y, Luo M H, Liu S Z, Gu Y Q, Li H J, Liu Z Y. A rare single nucleotide variant inconfers powdery mildew resistance in common wheat., 2020, 228: 1011–1026.

    [6] Koller T, Brunner S, Herren G, Hurni S, Keller B. Pyramiding of transgenicalleles in wheat results in improved powdery mildew resistance in the field., 2018, 131: 861–871.

    [7] Hall S A, Allen R L, Baumber R E, Baxter L A, Fisher K, Bittner-Eddy P D, Rose L E, Holub E B, Beynon J L. Maintenance of genetic variation in plants and pathogens involves complex networks of gene-for-gene interactions., 2009, 10: 449–457.

    [8] Allen R L, Meitz J C, Baumber R E, Hall S A, Lee S C, Rose L E, Beynon J L. Natural variation reveals key amino acids in a downy mildew effector that alters recognition specificity by anresistance gene., 2008, 9: 511–523.

    [9] Rose L E, Bittner-Eddy P D, Langley C H, Holub E B, Beynon J L. The maintenance of extreme amino acid diversity at the disease resistance gene,, in., 2004, 166: 1517–1527.

    [10] Rentel M C, Leonelli L, Dahlbeck D, Zhao B, Staskawicz B J. Recognition of theeffector ATR13 triggers resistance against oomycete, bacterial, and viral pathogens., 2008, 105: 1091–1096.

    [11] Ramachandran S R, Yin C, Kud J, Tanaka K, Mahoney A K, Xiao F, Hulbert S H. Effectors from wheat rust fungi suppress multiple plant defense responses., 2017, 107: 75–83.

    [12] Hu W G, Wang Q H, Wang S W, Wang M M, Wang C Y, Tian Z R, Liu X L, Ji W Q, Zhang H. Gene co-expression network analysis provides a novel insight into the dynamic response of wheat to powdery mildew stress., 2020, 99: 44–45.

    [13] Liu X Y, Zhang C Q, Zhang L L, Huang J G, Dang C, Xie C J, Wang Z Y., a CC-NBS-LRR like gene located on chr 7D, promotes disease resistance to wheat powdery mildew in Brock., 2020, 168: 688–699.

    [14] 王艷紅, 肖瑩, 鄭舒揚, 劉曉穎, 王振英. 白粉菌誘導的小麥品種Brock的差異表達基因解析. 天津師范大學學報(自然科學版), 2015, 35(2): 71–76.

    Wang Y H, Xiao Y, Zheng S Y, Liu X Y, Wang Z Y. Analysis on differentially expressed genes of wheat Brock inoculated byf. sp..(Nat Sci Edn), 2015, 35(2): 71–76 (in Chinese with English abstract).

    [15] Ramachandran S R, Yin C, Kud J, Tanaka K, Mahoney A K, Xiao F, Hulbert S H. Effectors from wheat rust fungi suppress multiple plant defense responses., 2016, 107: 75–83.

    [16] Chen J H, Zhang F Y, Zhao C J, Lyu G G, Sun C W, Pan Y B, Guo X Y, Chen F. Genome-wide association study of six quality traits reveals the association of thegene with flour colour in Chinese bread wheat., 2019, 17: 2106–2122.

    [17] 劉登全, 王園秀, 崔朝宇, 秦雙林, 歐陽慧, 蔣軍喜. 晚熟溫州蜜柑 NBS-LRR類抗病基因同源序列的克隆及分析. 江西農(nóng)業(yè)大學學報, 2016, 38: 83–89.

    Liu D Q, Wang Y X, Cui C Y, Qin S L, Ou-Yang H, Jiang J X. Cloning and analysis of NBS-LRR disease-resistant gene analogs inMarc., 2016, 38: 83–89 (in Chinese with English abstract).

    [18] Prabhukarthikeyan S R, Manikandan R, Durgadevi D, Keerthana U, Harish S, Karthikeyan G, Raguchander T. Bio-suppression of turmeric rhizome rot disease and understanding the molecular basis of tripartite interaction among,and., 2017, 111: 23–31.

    [19] Chen T, Lv Y, Zhao T M, Li N, Yang Y w, Yu W, He X, Liu T, Zhang B L. Comparative transcriptome profiling of a resistant vs. susceptible tomato () cultivar in response to infection by tomato yellow leaf curl virus., 2013, 8: e80816.

    [20] Han L J, Weng K, Ma H, Xiang G Q, Li Z Q, Wang Y J, Liu G T, Xu Y. Identification and characterization of erysiphe necator-responsive MicroRNAs in Chinese wild vitis pseudoreticulata by high-throughput sequencing., 2016, 7: 621.

    [21] Khan S A, Zhang C, Ali N, Gandeka M. Highdensity SNP map facilitates fine mapping of QTLs and candidate genes discovery forresistance in peanut ()., 2020, 133: 2239–2257.

    [22] Howlader J, Robin A, Natarajan S, Biswas M K, Nou I S. Transcriptome analysis by RNA-Seq reveals genes related to plant height in two sets of parent-hybrid combinations in easter lily ()., 2020, 10: 9082–9097.

    [23] 黃建國. 小麥抗白粉病基因 SSR 分子標記篩選及AP2/ERF轉(zhuǎn)錄因子家族分析. 天津師范大學碩士學位論文, 天津, 2021.

    Huang J G. Screening for SSR Markers Linked to Powdery Mildew Resistance Gene and Bioinformatics Analysis of AP2/ERF Transcription Factor Family in Wheat. MS Thesis of Tianjin Normal University, Tianjin, China, 2021 (in Chinese with English abstract).

    [24] 劉曉穎. 小麥品種Brock響應白粉菌侵染早期抗病相關基因表達和功能分析. 中國農(nóng)業(yè)大學博士學位論文, 北京, 2021.

    Liu X Y. Expression and Functional Analysis of Genes Associated with Early Powdery Mildew Resistance in Wheat Variety Brock. PhD Dissertation of China Agricultural University, Beijing, China, 2021 (in Chinese with English abstract).

    [25] Livak K J, Schmittgen T D. Analysis of relative gene expression data using real-time quantitative PCR and the 2?ΔΔCTmethod., 2001, 25: 402–408.

    [26] Liu X Y, Wang J Y, Fan B L, Shang Y T, Sun Y F, Dang C, Xie C J, Wang Z Y, Peng Y K. Agene in wheat contributes to the early defence response against wheat powdery mildew., 2018, 166: 116–122.

    [27] 栗聰, 雒景吾, 張磊, 田增榮, 劉新倫, 吉萬全. 小麥成熟胚再生體系優(yōu)化及優(yōu)良受體基因型篩選. 麥類作物學報, 2014, 34: 583–590.

    Li C, Luo J W, Zhang L, Tian Z R, Liu X L, Ji W Q. Optimizing the regeneration system from mature embryo and screening of elite wheat genotypes., 2014, 34: 583–590 (in Chinese with English abstract).

    [28] 王艷麗, 葉興國, 劉艷鵬, 杜麗璞, 徐惠君. 農(nóng)桿菌敏感小麥基因型的篩選研究. 麥類作物學報, 2005, 25(6): 6–10.

    Wang Y L, Ye X G, Liu Y P, Du L P, Xu H J. Screening of wheat genotype sensitive toinfection., 2005, 25(6): 6–10 (in Chinese with English abstract).

    [29] 王怡, 彭永康. 抗、感白粉病春小麥感染白粉菌后抗病特性的細胞學觀察. 天津農(nóng)業(yè)科學, 2014, 20(4): 1–4.

    Wang Y, Peng Y K. Cytological analysis of spring wheat with different resistance to powdery mildew induced by Bgt., 2014, 20(4): 1–4 (in Chinese with English abstract).

    [30] 河南省植保站. 河南省主要農(nóng)作物病蟲害測報辦法. 鄭州: 河南科學技術出版社, 1995. pp 11–18.

    Plant Protection and Quarantine Station of Henan Province. Forecasting Methods of the Main Crop Diseases and Insect Pests in Henan Province. Zhengzhou: Henan Science and Technology Press, 1995. pp 11–18 (in Chinese).

    [31] Yahiaoui N, Srichumpa P, Dudler R, Keller B. Genome analysis at different ploidy levels allows cloning of the powdery mildew resistance genefrom hexaploid wheat., 2004, 37: 528–538.

    [32] Brunner S, Hurni S, Herren G, Kalinina O, Burg S V, Zeller S L, Schmid B, Winzeler M, Keller B. Transgenicwheat lines show resistance to powdery mildew in the field., 2011, 9: 897–910.

    [33] Brunner S, Stirnweis D, Quijano C D, Buesing G, Herren G, Parlange F, Barret P, Tassy C, Sautter C, Winzeler M. Transgenicmultilines of wheat show increased powdery mildew resistance in the field., 2012, 10: 398–409.

    [34] He H G, Zhu S Y, Zhao R H, Jiang Z N, Ji Y Y, Ji J, Qiu D, Li H J, Bie T D., encoding a typical CC-NBS-LRR protein, confers broad-spectrum resistance to wheat powdery mildew disease., 2018, 11: 879–882.

    [35] 王海燕, 肖進, 袁春霞, 徐濤, 于春艷, 孫昊杰, 陳佩度, 王秀娥. 攜帶抗白粉病基因的小麥–簇毛麥小片段易位染色體在不同小麥背景中的傳遞率及遺傳穩(wěn)定性. 作物學報, 2016, 42: 361–367.

    Wang H Y, Xiao J, Yuan C X, Xu T, Yu C Y, Sun H J, Chen P D, Wang X E. Transmission and genetic stability of no-homoeologous small fragment wheat–translocation chromosomes within various cultivar backgrounds of common wheat., 2016, 42: 361–367 (in Chinese with English abstract).

    [36] Sánchez-Martín J, Steuernagel B, Ghosh S, Herren G, Hurni S, Adamski N, Vrána J, Kubaláková M, Krattinger S G, Wicker T, Dole?el J, Keller B, Wulff B B H. Rapid gene isolation in barley and wheat by mutant chromosome sequencing., 2016, 17: 221.

    [37] 陳芳, 喬麟軼, 李銳, 劉成, 李欣, 郭慧娟, 張樹偉, 常利芳, 李東方, 閻曉濤, 任永康, 張曉軍, 暢志堅. 小麥新種質(zhì)CH1357抗白粉病遺傳分析及染色體定位. 作物學報, 2019, 45: 1503–1510.

    Chen F, Qiao L Y, Li R, Li C, Li X, Guo H J, Zhang S W, Chang L F, Li D F, Yan X T, Ren Y K, Zhang X J, Chang Z J. Genetic analysis and chromosomal localization of powdery mildew resistance gene in wheat germplasm CH1357., 2019, 45: 1503–1510 (in Chinese with English abstract).

    [38] 陳芳. 小偃麥衍生品系抗白粉病基因的遺傳定位及圖位克隆. 山西大學博士論文, 山西太原, 2020.

    Chen F. Genetic Mapping and Positional Cloning of Powdery Mildew Resistance Genes in the Derived Lines from a Wheat-Hybrid. PhD Dissertation of Shanxi University, Taiyuan, Shanxi, China, 2020 (in Chinese with English abstract).

    [39] 李根橋, 房體麟, 朱婕, 高亮亮, 李閃, 解超杰, 楊作民, 孫其信, 劉志勇. 普通小麥品種Brock抗白粉病基因分子標記定位.作物學報, 2009, 35: 1613–1619.

    Li G Q, Fang T L, Zhu J, Gao L L, Li S, Xie C J, Yang Z M, Sun Q X, Liu ZY. Molecular identification of a powdery mildew resistance gene from common wheat cultivar Brock., 2009, 35: 1613–1619 (in Chinese with English abstract).

    [40] Zou S H, Wang H, Li Y W, Kong Z S, Tang D Z. The NB-LRR geneconfers powdery mildew resistance in wheat., 2017, 218: 298–309.

    [41] Zhou F S, Kurth J, Wei F S, Elliott C, Valè G, Yahiaoui N, Keller B, Somerville S, Wise R, Schulze-Lefert P. Cell-autonomous expression of Barleyconfers race-specific resistance to the powdery mildew fungus via aindependent signaling pathway., 2001, 13: 337–350.

    [42] Halterman D A, Wei F, Wise R P. Powdery mildew-inducedmRNAs are alternatively spliced and contain multiple upstream open reading frames., 2003, 131: 558–567.

    [43] Chang C, Yu D S, Jiao J, Jing S J, Schulze-Lefert P, Shen Q H. Barley MLA immune receptors directly interfere with antagonistically acting transcription factors to initiate disease resistance signaling., 2013, 25: 1158–1173.

    Cloning and functional analysis ofgene related to powdery mildew resistance in wheat cultivar Brock

    LIU Xiao-Ying1, ZHANG Chi1, WANG Xue-Qing1, YANG Chen-Xiao1, WANG Guang-Yu1, BIAN Yun-Di1, FANG Fang2, WANG Ying2, and WANG Zhen-Ying1,*

    1College of Life Sciences, Tianjin Normal University, Tianjin Key Laboratory of Animal and Plant Resistance, Tianjin 300387, China;2Forestry Development Service Center of Baodi District, Tianjin 301899, China

    Powdery mildew caused byf. sp.() is a severe wheat disease in China. Cloning and pyramiding of different resistance genes to improve crop disease resistance is one of the most cost-effective methods. In this study,on chromosome 1B, which encodes the CC, NB-ARC, and LRR domains, was isolated from common wheat. The relative expression level ofin Brock and BJ-1 fluctuated afterinoculation, but the average expression levels were always higher than the susceptible wheat Jing 411. The function ofwas elucidated by virus-induced gene silencing (VIGS) and overexpression transgenic technique. Silencing ofresulted in decreased disease resistance in Brock. Overexpression ofimproved disease resistance in the transgenic wheat seedlings of Jinqiang 5 cultivar. The above results demonstrated thatwas involved in the defense response of wheat to powdery mildew, which provided valuable genetic resources for breeding of resistant varieties.

    wheat;gene;; powdery mildew resistance

    10.3724/SP.J.1006.2023.21003

    本研究由國家自然科學基金項目(31071671)和天津市科技支撐項目(18YFZCNC01100, 17JCZDJC34100)資助。

    This study was supported by the National Natural Science Foundation of China (31071671) and the Tianjin Science and Technology Program (18YFZCNC01100, 17JCZDJC34100).

    王振英, E-mail: skywangzy@tjnu.edu.cn

    E-mail: skylxy@tjnu.edu.cn

    2022-01-16;

    2022-05-05;

    2022-05-25.

    URL: https://kns.cnki.net/kcms/detail/11.1809.S.20220525.1330.002.html

    This is an open access article under the CC BY-NC-ND license (http://creativecommons.org/licenses/by-nc-nd/4.0/).

    猜你喜歡
    白粉白粉病抗病
    我國小麥基因組編輯抗病育種取得突破
    蔬菜白粉虱 秋防正當時
    一到春季就流行 蔬菜白粉病該咋防
    保護地蔬菜白粉虱的發(fā)生與防治
    bZIP轉(zhuǎn)錄因子在植物激素介導的抗病抗逆途徑中的作用
    拉薩設施月季白粉病的發(fā)生與防治
    西藏科技(2016年8期)2016-09-26 09:00:21
    葡萄新品種 優(yōu)質(zhì)又抗病
    白僵菌防治溫室大棚白粉虱藥效試驗
    黃瓜白粉病的發(fā)生與防治
    八十三團紅提葡萄白粉病的發(fā)生與防治
    日本与韩国留学比较| 国产精品1区2区在线观看.| av在线蜜桃| 高清日韩中文字幕在线| 欧美成人a在线观看| 在线免费十八禁| 中文字幕人妻熟人妻熟丝袜美| 国内精品一区二区在线观看| 亚洲欧美日韩无卡精品| 男女之事视频高清在线观看| 久久99热6这里只有精品| 免费一级毛片在线播放高清视频| 日本-黄色视频高清免费观看| 国产精品av视频在线免费观看| 一个人免费在线观看电影| 久久精品国产亚洲网站| 国产成人一区二区在线| av中文乱码字幕在线| 97在线视频观看| 亚洲最大成人手机在线| 日韩大尺度精品在线看网址| 国产精品综合久久久久久久免费| АⅤ资源中文在线天堂| 99在线视频只有这里精品首页| av在线老鸭窝| av在线蜜桃| 国产精品一区二区三区四区久久| 亚洲国产色片| 一级毛片aaaaaa免费看小| 国产一级毛片七仙女欲春2| 亚洲美女视频黄频| 午夜福利在线观看吧| 欧洲精品卡2卡3卡4卡5卡区| 国产亚洲精品久久久久久毛片| 热99在线观看视频| 99riav亚洲国产免费| 国产久久久一区二区三区| 大香蕉久久网| 我的女老师完整版在线观看| 亚洲国产色片| 中国美白少妇内射xxxbb| 日韩精品青青久久久久久| 久久精品综合一区二区三区| 免费在线观看成人毛片| 久久精品国产亚洲av涩爱 | 亚洲精品一区av在线观看| 欧美性猛交黑人性爽| 简卡轻食公司| 在线播放无遮挡| 亚洲真实伦在线观看| 久久久精品欧美日韩精品| 欧美色欧美亚洲另类二区| 中文字幕av在线有码专区| 在线观看av片永久免费下载| 国产一区二区在线av高清观看| 午夜爱爱视频在线播放| 免费一级毛片在线播放高清视频| 日韩av不卡免费在线播放| 床上黄色一级片| 亚洲欧美中文字幕日韩二区| 91在线观看av| 在线观看一区二区三区| 午夜福利视频1000在线观看| 免费人成在线观看视频色| 国产高潮美女av| 精品乱码久久久久久99久播| 午夜爱爱视频在线播放| 国产免费男女视频| 亚洲综合色惰| 99久久精品国产国产毛片| 男女那种视频在线观看| 国产乱人视频| 一级毛片电影观看 | 久久久精品94久久精品| 成人鲁丝片一二三区免费| 国产成人一区二区在线| 丰满的人妻完整版| 亚洲av免费高清在线观看| 久久久精品大字幕| 蜜臀久久99精品久久宅男| 国产精品精品国产色婷婷| 在线看三级毛片| 色5月婷婷丁香| 日韩欧美精品v在线| 天美传媒精品一区二区| 欧洲精品卡2卡3卡4卡5卡区| 国产女主播在线喷水免费视频网站 | 综合色av麻豆| 欧美中文日本在线观看视频| 变态另类成人亚洲欧美熟女| 亚洲精品粉嫩美女一区| 成人鲁丝片一二三区免费| 成人高潮视频无遮挡免费网站| 欧美bdsm另类| 免费看美女性在线毛片视频| 国产黄色视频一区二区在线观看 | 亚洲欧美精品自产自拍| 久久欧美精品欧美久久欧美| 男人舔女人下体高潮全视频| 露出奶头的视频| 日日摸夜夜添夜夜爱| 成人av在线播放网站| 精品久久久久久成人av| 免费av观看视频| 国内精品美女久久久久久| 成人三级黄色视频| 97在线视频观看| 日本一二三区视频观看| 国产精品久久久久久久久免| 国产探花在线观看一区二区| 看十八女毛片水多多多| 免费看a级黄色片| 国产av不卡久久| 天天躁日日操中文字幕| 国产单亲对白刺激| 在线国产一区二区在线| 全区人妻精品视频| or卡值多少钱| 亚洲av不卡在线观看| 网址你懂的国产日韩在线| www.色视频.com| 国内精品宾馆在线| 国产成人一区二区在线| 在线免费观看的www视频| 1024手机看黄色片| 精品少妇黑人巨大在线播放 | 国产高清视频在线播放一区| 高清日韩中文字幕在线| 少妇裸体淫交视频免费看高清| 毛片女人毛片| 免费看日本二区| 亚洲图色成人| 亚洲人成网站高清观看| 色综合色国产| 欧美日韩在线观看h| АⅤ资源中文在线天堂| 亚洲精品久久国产高清桃花| aaaaa片日本免费| 啦啦啦韩国在线观看视频| 国产老妇女一区| 美女免费视频网站| 综合色av麻豆| 日本三级黄在线观看| 日本-黄色视频高清免费观看| 欧美一区二区国产精品久久精品| 国产精品久久电影中文字幕| 日韩精品青青久久久久久| 精品免费久久久久久久清纯| 黄色欧美视频在线观看| 色噜噜av男人的天堂激情| 日韩高清综合在线| 久久久久久久亚洲中文字幕| 看片在线看免费视频| av福利片在线观看| 久久久久久久午夜电影| 别揉我奶头 嗯啊视频| 日本-黄色视频高清免费观看| 露出奶头的视频| 校园春色视频在线观看| 日本免费一区二区三区高清不卡| av卡一久久| 亚洲丝袜综合中文字幕| 国产精品av视频在线免费观看| 亚洲欧美日韩高清专用| 日韩精品中文字幕看吧| 长腿黑丝高跟| 成年av动漫网址| 不卡一级毛片| 女人十人毛片免费观看3o分钟| 观看免费一级毛片| 午夜福利18| 美女内射精品一级片tv| 国产精品野战在线观看| 国产一区二区在线观看日韩| 国产精品一区www在线观看| 亚洲欧美日韩东京热| av免费在线看不卡| 成年女人毛片免费观看观看9| 九九在线视频观看精品| 亚洲性夜色夜夜综合| 欧美高清性xxxxhd video| 麻豆久久精品国产亚洲av| 欧美一区二区亚洲| 久久精品国产亚洲av天美| 国产高潮美女av| 最近中文字幕高清免费大全6| 精品人妻一区二区三区麻豆 | 插阴视频在线观看视频| 内射极品少妇av片p| 色吧在线观看| 亚洲欧美精品综合久久99| 欧美三级亚洲精品| 亚洲aⅴ乱码一区二区在线播放| 淫妇啪啪啪对白视频| 国产三级在线视频| 女人十人毛片免费观看3o分钟| 午夜老司机福利剧场| 精华霜和精华液先用哪个| 精品午夜福利在线看| 久久午夜亚洲精品久久| 天堂网av新在线| 伦理电影大哥的女人| ponron亚洲| 免费不卡的大黄色大毛片视频在线观看 | 日韩欧美一区二区三区在线观看| 亚洲国产精品成人综合色| 欧美成人精品欧美一级黄| 乱码一卡2卡4卡精品| 中文资源天堂在线| av.在线天堂| 自拍偷自拍亚洲精品老妇| 一级毛片久久久久久久久女| 桃色一区二区三区在线观看| 欧美+日韩+精品| 亚洲精品成人久久久久久| 麻豆国产av国片精品| 亚洲av成人精品一区久久| 亚洲,欧美,日韩| 国产久久久一区二区三区| 欧美日本视频| 国产精华一区二区三区| 亚洲图色成人| 中文资源天堂在线| 国产av在哪里看| 美女xxoo啪啪120秒动态图| 97超级碰碰碰精品色视频在线观看| 亚洲精华国产精华液的使用体验 | av视频在线观看入口| 老司机午夜福利在线观看视频| 久久久久免费精品人妻一区二区| av国产免费在线观看| 国产欧美日韩一区二区精品| 三级男女做爰猛烈吃奶摸视频| 中文字幕熟女人妻在线| 大又大粗又爽又黄少妇毛片口| 一级a爱片免费观看的视频| 在线免费十八禁| 三级经典国产精品| 偷拍熟女少妇极品色| 国产色爽女视频免费观看| 搡老岳熟女国产| 简卡轻食公司| 在线a可以看的网站| 久久久精品欧美日韩精品| 又黄又爽又刺激的免费视频.| 无遮挡黄片免费观看| 99久国产av精品国产电影| 国产成人freesex在线 | 亚洲经典国产精华液单| 午夜亚洲福利在线播放| 免费av观看视频| 深夜精品福利| 老师上课跳d突然被开到最大视频| 国产精品福利在线免费观看| 狂野欧美白嫩少妇大欣赏| 九九在线视频观看精品| 欧美在线一区亚洲| 国产免费男女视频| 成人漫画全彩无遮挡| 亚洲人成网站高清观看| 国产私拍福利视频在线观看| 欧美丝袜亚洲另类| 国内少妇人妻偷人精品xxx网站| 夜夜爽天天搞| 内射极品少妇av片p| 男女边吃奶边做爰视频| 日韩大尺度精品在线看网址| 97超视频在线观看视频| 色尼玛亚洲综合影院| 亚洲精品一区av在线观看| 欧美色视频一区免费| 精品福利观看| а√天堂www在线а√下载| 欧美潮喷喷水| 哪里可以看免费的av片| 国产在线精品亚洲第一网站| 亚洲综合色惰| 亚洲第一电影网av| 色噜噜av男人的天堂激情| 麻豆乱淫一区二区| 中文字幕av成人在线电影| 插阴视频在线观看视频| 欧美日韩综合久久久久久| 午夜亚洲福利在线播放| 欧美人与善性xxx| 97在线视频观看| 99热这里只有是精品50| av免费在线看不卡| 日韩欧美国产在线观看| 久久婷婷人人爽人人干人人爱| 国产69精品久久久久777片| 精品人妻偷拍中文字幕| 99在线视频只有这里精品首页| 亚洲成人av在线免费| 成人午夜高清在线视频| 国内精品一区二区在线观看| av视频在线观看入口| 久久久久免费精品人妻一区二区| 中国国产av一级| 联通29元200g的流量卡| 亚洲熟妇熟女久久| 日韩大尺度精品在线看网址| 亚洲av成人av| 亚洲在线自拍视频| 久久久久久久久久黄片| av在线老鸭窝| 日日撸夜夜添| 一本精品99久久精品77| 久久久久久伊人网av| 欧美成人免费av一区二区三区| av中文乱码字幕在线| 精品人妻熟女av久视频| 99久久久亚洲精品蜜臀av| 亚洲欧美清纯卡通| 日本-黄色视频高清免费观看| 美女内射精品一级片tv| 在线a可以看的网站| 卡戴珊不雅视频在线播放| 国产伦一二天堂av在线观看| 18禁在线播放成人免费| 搡老熟女国产l中国老女人| 色综合站精品国产| 99久久精品一区二区三区| 国产精品一区二区三区四区久久| 成人综合一区亚洲| 国产男人的电影天堂91| 精品久久久噜噜| av免费在线看不卡| 99热这里只有是精品50| 国产男人的电影天堂91| 99精品在免费线老司机午夜| 麻豆国产97在线/欧美| 久久久精品大字幕| 亚洲av成人av| 成人欧美大片| 一级毛片电影观看 | 国产黄色视频一区二区在线观看 | 亚洲中文字幕日韩| 亚洲四区av| 熟女电影av网| 久久午夜亚洲精品久久| 两个人的视频大全免费| 能在线免费观看的黄片| 此物有八面人人有两片| 国产亚洲欧美98| 日本五十路高清| 午夜激情欧美在线| 91久久精品电影网| 一个人观看的视频www高清免费观看| 麻豆一二三区av精品| 在线免费观看的www视频| 久久99热6这里只有精品| 日本-黄色视频高清免费观看| 亚洲欧美清纯卡通| 国产成人精品久久久久久| 俺也久久电影网| 国产精品亚洲一级av第二区| 亚洲精品国产av成人精品 | 99热6这里只有精品| 一进一出抽搐gif免费好疼| 在线免费十八禁| 国产高清视频在线播放一区| 午夜激情福利司机影院| 色在线成人网| 中出人妻视频一区二区| 久久精品夜夜夜夜夜久久蜜豆| 一卡2卡三卡四卡精品乱码亚洲| 美女大奶头视频| 欧美成人精品欧美一级黄| 3wmmmm亚洲av在线观看| 直男gayav资源| 亚洲综合色惰| 亚洲欧美日韩高清在线视频| 国产精品三级大全| 亚洲精华国产精华液的使用体验 | 国产69精品久久久久777片| 午夜福利18| 国产精品伦人一区二区| 亚洲成人中文字幕在线播放| 精品久久久久久成人av| 色吧在线观看| 成人一区二区视频在线观看| 亚洲无线在线观看| 久久久久免费精品人妻一区二区| 国产高清有码在线观看视频| 欧美不卡视频在线免费观看| 日韩欧美精品免费久久| 欧美区成人在线视频| 97热精品久久久久久| 全区人妻精品视频| 久久久国产成人免费| 日本免费一区二区三区高清不卡| a级毛色黄片| 亚洲自拍偷在线| 欧美成人精品欧美一级黄| 少妇丰满av| 欧美高清性xxxxhd video| 色综合亚洲欧美另类图片| 最近最新中文字幕大全电影3| 国产成年人精品一区二区| 精品久久久噜噜| 有码 亚洲区| 小蜜桃在线观看免费完整版高清| 美女xxoo啪啪120秒动态图| 不卡视频在线观看欧美| 国产高清有码在线观看视频| 欧美日韩综合久久久久久| 国产亚洲av嫩草精品影院| 亚洲欧美精品自产自拍| 在线免费观看不下载黄p国产| 少妇被粗大猛烈的视频| 欧美高清性xxxxhd video| 国产成人91sexporn| 国产在线精品亚洲第一网站| 欧美一区二区国产精品久久精品| 国产伦一二天堂av在线观看| 男女视频在线观看网站免费| 亚洲中文字幕日韩| 国产黄色视频一区二区在线观看 | 欧美三级亚洲精品| 国产麻豆成人av免费视频| 国产男靠女视频免费网站| 亚洲精品456在线播放app| 黄片wwwwww| 亚洲国产精品国产精品| 99在线人妻在线中文字幕| 午夜日韩欧美国产| 成人午夜高清在线视频| 在线天堂最新版资源| 麻豆成人午夜福利视频| 国产精品福利在线免费观看| 亚洲精品乱码久久久v下载方式| 91麻豆精品激情在线观看国产| 黄色欧美视频在线观看| 听说在线观看完整版免费高清| 亚洲最大成人中文| 亚洲欧美清纯卡通| 男女啪啪激烈高潮av片| 久久草成人影院| 久久精品国产亚洲av香蕉五月| 十八禁国产超污无遮挡网站| 亚洲精品国产成人久久av| 国产探花极品一区二区| 精品久久久久久久末码| 亚洲成人久久性| 久久精品国产亚洲网站| 成人国产麻豆网| av在线老鸭窝| 中文字幕熟女人妻在线| 网址你懂的国产日韩在线| 俺也久久电影网| 国产午夜福利久久久久久| 国产麻豆成人av免费视频| 国产 一区精品| 欧美高清性xxxxhd video| 免费不卡的大黄色大毛片视频在线观看 | 国内揄拍国产精品人妻在线| 亚洲av成人精品一区久久| 欧美性猛交黑人性爽| 久久鲁丝午夜福利片| 18禁在线无遮挡免费观看视频 | 麻豆av噜噜一区二区三区| 在线观看一区二区三区| 国内精品久久久久精免费| 日韩强制内射视频| 伦理电影大哥的女人| 久久国内精品自在自线图片| 此物有八面人人有两片| 日韩亚洲欧美综合| 国产高清视频在线播放一区| 桃色一区二区三区在线观看| 日本 av在线| 久久久精品大字幕| 久久综合国产亚洲精品| 在线观看免费视频日本深夜| 国产欧美日韩一区二区精品| 国产精品福利在线免费观看| 亚洲国产精品久久男人天堂| 久久精品影院6| 欧美人与善性xxx| videossex国产| 国产高潮美女av| 国产单亲对白刺激| 久久久久久久久久黄片| 日日摸夜夜添夜夜添av毛片| 69av精品久久久久久| 亚洲中文字幕一区二区三区有码在线看| 久久人人爽人人片av| 成人国产麻豆网| 联通29元200g的流量卡| 91在线精品国自产拍蜜月| 欧美日韩一区二区视频在线观看视频在线 | 97热精品久久久久久| 国内少妇人妻偷人精品xxx网站| 蜜桃亚洲精品一区二区三区| 国产精品一及| 亚洲精华国产精华液的使用体验 | 天堂影院成人在线观看| 日韩强制内射视频| 国产高清不卡午夜福利| 三级毛片av免费| 在线观看一区二区三区| 美女被艹到高潮喷水动态| 狠狠狠狠99中文字幕| av.在线天堂| 99热这里只有是精品在线观看| 日本 av在线| 成人欧美大片| 免费搜索国产男女视频| 国产高潮美女av| 女的被弄到高潮叫床怎么办| 在线免费观看的www视频| 尾随美女入室| 午夜老司机福利剧场| 久久欧美精品欧美久久欧美| 久久精品国产亚洲av天美| 男女边吃奶边做爰视频| 国产黄片美女视频| 少妇人妻精品综合一区二区 | 精品久久国产蜜桃| 白带黄色成豆腐渣| 亚洲人与动物交配视频| 草草在线视频免费看| 免费人成在线观看视频色| av视频在线观看入口| 国产女主播在线喷水免费视频网站 | 1000部很黄的大片| 人妻久久中文字幕网| 变态另类丝袜制服| 亚洲色图av天堂| 国产成人aa在线观看| 色噜噜av男人的天堂激情| 啦啦啦韩国在线观看视频| av视频在线观看入口| 日韩欧美 国产精品| 国产精品久久久久久久久免| 免费高清视频大片| 国产精品国产高清国产av| 久99久视频精品免费| 成人国产麻豆网| 午夜视频国产福利| 亚洲精品影视一区二区三区av| 又粗又爽又猛毛片免费看| 日本三级黄在线观看| 国产精品爽爽va在线观看网站| 成人亚洲欧美一区二区av| 久久久久性生活片| 日日撸夜夜添| 国产伦一二天堂av在线观看| 国内精品宾馆在线| 久久精品国产自在天天线| 国产激情偷乱视频一区二区| 免费搜索国产男女视频| 此物有八面人人有两片| 国产综合懂色| 美女xxoo啪啪120秒动态图| 亚洲第一区二区三区不卡| 国产精品伦人一区二区| 亚洲内射少妇av| 婷婷精品国产亚洲av| 久久久久性生活片| 日本精品一区二区三区蜜桃| 人妻制服诱惑在线中文字幕| 国产淫片久久久久久久久| 久久精品夜夜夜夜夜久久蜜豆| 中文字幕av在线有码专区| 在线观看午夜福利视频| 亚洲精品一区av在线观看| 99热这里只有是精品50| 女人十人毛片免费观看3o分钟| 久久热精品热| 国产在视频线在精品| 免费看美女性在线毛片视频| 色哟哟哟哟哟哟| 亚洲欧美精品综合久久99| 男女边吃奶边做爰视频| 国产老妇女一区| 欧美中文日本在线观看视频| 变态另类丝袜制服| av福利片在线观看| 欧美3d第一页| 床上黄色一级片| 亚洲专区国产一区二区| 夜夜夜夜夜久久久久| 国产高清视频在线观看网站| 黄色配什么色好看| 一进一出抽搐动态| av免费在线看不卡| 看免费成人av毛片| 非洲黑人性xxxx精品又粗又长| 99热网站在线观看| 伊人久久精品亚洲午夜| 亚洲人与动物交配视频| 美女黄网站色视频| 国产精品一区二区三区四区免费观看 | 精品不卡国产一区二区三区| 国产单亲对白刺激| 久久精品综合一区二区三区| 国产熟女欧美一区二区| 亚洲av不卡在线观看| or卡值多少钱| 身体一侧抽搐| 欧洲精品卡2卡3卡4卡5卡区| 国产欧美日韩精品亚洲av| 日韩成人伦理影院| 国产成人a∨麻豆精品| 亚洲av成人精品一区久久| 欧美日韩在线观看h| 欧美国产日韩亚洲一区| 十八禁网站免费在线| 97超碰精品成人国产| 日本爱情动作片www.在线观看 | 国产精品人妻久久久影院| 99久久精品一区二区三区| 国产黄色小视频在线观看|