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

    煙草NtWRKY40在植物應(yīng)答病毒侵染過程中的作用

    2016-11-09 02:22:50劉晶晶程春玲席玉珍魏書
    生物技術(shù)通報 2016年10期
    關(guān)鍵詞:前體侵染轉(zhuǎn)基因

    劉晶晶 程春玲 席玉珍 魏書

    (安徽農(nóng)業(yè)大學(xué)茶樹生物學(xué)與資源利用國家重點實驗室,合肥 230036)

    煙草NtWRKY40在植物應(yīng)答病毒侵染過程中的作用

    劉晶晶 程春玲 席玉珍 魏書

    (安徽農(nóng)業(yè)大學(xué)茶樹生物學(xué)與資源利用國家重點實驗室,合肥 230036)

    WRKY轉(zhuǎn)錄因子家族在植物的抗病、抗逆反應(yīng)中具有重要功能。已有研究表明,煙草花葉病毒(TMV)的侵染顯著地誘導(dǎo)煙草NtWRKY的表達,有必要進一步探明該基因在植物應(yīng)答病毒侵染過程中的作用。采用PCR的方法克隆獲得NtWRKY cDNA,生物信息學(xué)分析結(jié)果顯示,該基因?qū)儆赪RKY Ⅱa 亞族成員,與絨毛狀煙草NtoWRKY40高度同源,命名為NtWRKY40。以此建立了過表達該基因的轉(zhuǎn)基因煙草,并以TMV為毒源進行了轉(zhuǎn)基因煙草和野生煙草的侵染實驗,以觀察NtWRKY40在煙草應(yīng)答病毒侵染過程中的作用。實驗結(jié)果表明,野生煙草在TMV侵染后9 d,NtWRKY40的表達量顯著升高,而NtWRKY40過表達轉(zhuǎn)基因煙草在病毒侵染后,病毒相關(guān)基因的表達高于野生型對照,與染病程度成正相關(guān),說明過表達NtWRKY40增加了植株對病毒的敏感性,該基因為負調(diào)控因子。此外,為探索應(yīng)用人工miRNA的抗病毒技術(shù),以煙草天然miR167前體為骨架、馬鈴薯Y病毒(PVY)外殼蛋白基因的一段反向互補序列為成熟序列,構(gòu)建了amiR167-PVY植物表達載體并轉(zhuǎn)化煙草,以抑制PVY。對amiRNA轉(zhuǎn)基因植株進行抗病毒實驗的結(jié)果顯示,amiR167-PVY能夠部分抑制病毒基因的表達,轉(zhuǎn)基因植株具有一定的抗病毒能力。

    煙草花葉病毒;馬鈴薯Y病毒;WRKY轉(zhuǎn)錄因子;人工miRNA;抗病毒

    DOI:10.13560/j.cnki.biotech.bull.1985.2016.10.023

    煙草花葉病毒(Tobacco mosaic virus,TMV)和馬鈴薯Y病毒(Potato virus Y,PVY)均為RNA病毒,寄主廣泛,可危害煙草、馬鈴薯、番茄和辣椒等多種茄科作物,造成巨大經(jīng)濟損失[1]。

    WRKY(WRKY transcription factors,WRKY TFs)等多種轉(zhuǎn)錄因子參與植物對生物逆境的防衛(wèi)[2]。WRKY轉(zhuǎn)錄因子的氨基端包含保守的WRKY結(jié)構(gòu)域,該結(jié)構(gòu)域約含有60個氨基酸,由1個七肽WRKYGQR序列和一個鋅指結(jié)構(gòu)組成,WRKY結(jié)構(gòu)域能夠與靶基因啟動子上的W-box元件[(C/T)TGAC(C/T)序列]特異性結(jié)合,實現(xiàn)對靶基因的調(diào)控作用[3]。根據(jù)WRKY 結(jié)構(gòu)域數(shù)目和鋅指結(jié)構(gòu)特征,將WRKY 轉(zhuǎn)錄因子家族分為類型I、IIa+ IIb、IIc、IId+IIe和III[4]。在擬南芥中的研究表明,WRKY40是一個負調(diào)控植物白粉病抗性的轉(zhuǎn)錄因子,它的表達會抑制抗病信號通路EDSl和抗病物質(zhì)亞麻薺素(camalexin)的合成[5]。在煙草中過表達苜蓿的WRKY基因,能夠顯著誘導(dǎo)NtPR2(endo-1,3-beta-glucanase)基因的表達,增強轉(zhuǎn)基因煙草對煙草花葉病毒抗性反應(yīng)[6]。在煙草中,WRKY8與植物的免疫相關(guān)[7,8],WRKY3和WRKY6參與JA(茉莉酸)信號調(diào)控的防御反應(yīng),將這兩個基因單獨或同時沉默,茉莉酸的合成顯著下降,加重?zé)煵萏於甑那趾Γ?]。也有報道指出,煙草的WRKY1和WRKY2基因能夠迅速響應(yīng)機械損傷和真菌侵害[10,11]。研究表明,病原物相關(guān)分子模式觸發(fā)的免疫(PTI)和效應(yīng)子觸發(fā)的免疫(ETI)可以誘導(dǎo)大部分WRKY轉(zhuǎn)錄因子在轉(zhuǎn)錄水平上的表達,例如煙草的WRKY7、8、9、11[12]。我們曾發(fā)現(xiàn),當煙草花葉病毒(Tobacco mosaic virus,TMV)、馬鈴薯Y病毒(Potato virus Y,PVY)和黃瓜花葉病毒(Cucu-mber mosaic virus,CMV)侵染普通煙草后,煙草幼葉的一種WRKY基因的表達激增[13]。因此,進一步研究該基因在煙草應(yīng)答病毒侵染過程中的作用,對闡明植物抗病信號轉(zhuǎn)導(dǎo)途徑及其作用機制具有重要作用。

    近年來發(fā)展的人工miRNA(artificial miRNA,amiRNA)的技術(shù),以生物體內(nèi)源miRNA前體作為天然骨架,以人為選定的靶基因的特定片段取代天然miRNA成熟序列,構(gòu)建amiRNA以達到抑制生物體內(nèi)靶基因表達的目的,實現(xiàn)特定基因的沉默,特別是用于共同抑制多個基因的表達[14,15]。這種技術(shù)以其特異性、高效性受到了普遍的關(guān)注,在動植物代謝工程、生長發(fā)育調(diào)控以及抵抗外源病毒等方面起著重要的作用[16]。但是由于難以明確所有amiRNA在體內(nèi)都能按預(yù)期被加工產(chǎn)生而具有顯著的抑制功能。因此,有必要對相關(guān)因素進行進一步的研究。

    為此,本研究在克隆獲得NtWRKY cDNA之后,建立該基因過表達轉(zhuǎn)基因煙草,并以TMV為毒源,以明確該基因在煙草應(yīng)答病毒侵染過程中的作用。此外,為探索應(yīng)用人工miRNA的抗病毒技術(shù),以煙草天然miR167e前體為骨架、PVY病毒外殼蛋白基因的一段反向互補序列為成熟序列,構(gòu)建35S啟動子驅(qū)動的植物表達載體35S∷amiRNA-PVY并轉(zhuǎn)化煙草,以抑制PVY病毒,旨在為植物基因表達的調(diào)控及其應(yīng)用技術(shù)提供科學(xué)依據(jù)。

    1 材料與方法

    1.1材料

    1.1.1植物材料 實驗所用煙草品種為云煙85,為本實驗室培養(yǎng)的無菌苗。

    1.1.2病毒材料 煙草花葉病毒(TMV)和馬鈴薯Y病毒普通株系(PVYO系)繁殖并保存于云煙草85上,由安徽農(nóng)業(yè)大學(xué)江彤副教授提供。

    1.2方法

    1.2.1DNA提取 煙草葉片DNA的提取采用CTAB法,DNA樣品經(jīng)RNase I 酶(TaKaRa)處理以去除RNA污染。DNA濃度用NanoDrop 2000超微量分光光度計(Thermo)測定。

    1.2.2RNA提取及cDNA的合成 煙草葉片總RNA(Total RNA)和小RNA(Small RNA)提取分別使用RNAiso Plus試劑盒(TaKaRa)和RNAiso Plus for Small RNA試劑盒(TaKaRa)并按使用說明提取??俁NA樣品經(jīng)DNase I 酶(TaKaRa)處理以去除DNA污染。RNA濃度用NanoDrop 2000超微量分光光度計(Thermo)測定,RNA的完整性通過1%瓊脂糖凝膠電泳檢測。cDNA第一鏈合成分別使用PrimeScript RT reagent Kit(TaKaRa)和SYBR PrimeScript miRNA RT-PCR Kit(TaKaRa)并按說明書進行。

    1.2.3生物信息學(xué)分析 在擬南芥信息資源網(wǎng)站TAIR上下載擬南芥WRKY家族的部分基因編碼蛋白的序列,包括AtWRKY1(AT2G04880)、AtWRKY4(AT1G13960)、AtWRKY6(AT1G62300)、AtWRKY9(AT1G68150)、AtWRKY10(AT1G55600)、AtWRKY14(AT1G30650)、AtWRKY15(AT2G23320)、AtWRKY17(AT2G24570)、AtWRKY36(AT1G69810)、AtWRKY40(AT1G80840)、AtWRKY56(AT1G64000)、AtWRKY57(AT1G69310)、AtWRKY60(AT2G25000)、AtWRKY61(AT1G18860)、AtWRKY63(AT1G66600)、AtWRKY64(AT1G66560)、AtWRKY65(AT1G29280)、AtWRKY67(AT1G66550)、AtWRKY71(AT1G29860),在NCBI數(shù)據(jù)庫下載其他植物WRKY家族基因編碼蛋白的序列,其中包括棉花的GhWRKY1(KF669831)、GhWRKY2(KF669759)、GhWRKY3(KF669771)、GhWRKY4(KF669822)、GhWRKY5(KF669781)、GhWRKY12(KF669853)、GhWRKY19(KF669784)、GhWRKY34(KF669842)、GhWRKY40(AGX24945.1),大豆的GmWRKY25(ABS18427.1)、GmWRKY40(ABC26914.1)、GmWRKY53(ABC26915.1)。運用MEGA5.10 軟件基于鄰位歸并法(NJ 法)構(gòu)建系統(tǒng)進化樹。

    運用DNAMAN軟件,對普通煙草(Nicotiana tobacum)NtWRKY40(AB024510.1)、絨毛狀煙草(N. tomentosiformis)NtoWRKY40(XM_009612054.1),美花煙草(N. sylvestris)NsWRKY40(XM_0097821 27.1)和番茄(Solanum lycopersicum)SlWRKY40(NM_001317914.1)進行多重序列比對。應(yīng)用WoLF PSORT(http://wolfpsort.org/)在線軟件對NtWRKY40進行亞細胞定位預(yù)測。

    1.2.4質(zhì)粒的構(gòu)建 以煙草cDNA為模板,NtWRKY40基因上下游引物進行克隆,獲得煙草NtWRKY40基因全長序列,通過引物上的Xba I和Sac I酶切位點,連接到植物雙元表達載體pBI121上,獲得pBI121-NtWRKY40植物表達載體。

    以病毒外殼蛋白為靶基因,以煙草天然miR167e為骨架[17],在網(wǎng)站W(wǎng)MD3[18](Web miRNA Designer,http://wmd3.weigelworld.org/)上設(shè)計amiRNA引物,通過重疊PCR(overlapping PCR)的方法,以煙草DNA為模板,先以amiR167-F與 amiR167-I、amiR167-II與amiR167-III和amiR167IV與amiR167-R這3對引物進行分段擴增,回收3段目的片段,以這3段回收片段為模板,用引物amiR167-F和amiR167-R進行融合PCR擴增,即得到amiR167-PVY序列。然后通過BamH I和Sac I酶切位點,連接到植物雙元表達載體pBI121上,獲得pBI121-amiR167-PVY植物表達載體。

    1.2.5基因的定量分析 檢測NtWRKY40基因、病毒相關(guān)基因和amiRNA前體(pre-amiR167-PVY)基因表達的實時熒光定量PCR是以總RNA反轉(zhuǎn)錄產(chǎn)物為模板,所用基因?qū)R恍砸?,見?。內(nèi)參基因為Ntubc2[19],熒光定量PCR用TransStart Top Green qPCR Mix試劑盒(TransGen)并按照說明書配制反應(yīng)體系,在Bio-Rad熒光定量PCR儀上擴增。反應(yīng)條件為:95℃預(yù)變性30 s;95℃變性5 s,55℃退火25 s,72℃延伸20 s,40個循環(huán);72℃延伸5min。數(shù)據(jù)讀取由實時熒光定量PCR儀自動完成,供試基因的擴增效率均在85%以上?;虮磉_水平計算采用2-ΔΔCt方法。試驗進行3次技術(shù)重復(fù),采用t檢驗進行顯著性差異分析。

    實時熒光定量PCR檢測amiRNA成熟序列。以小RNA反轉(zhuǎn)錄產(chǎn)物為模板,專一性引物:一條為試劑盒內(nèi)的通用引物Uni-miR qPCR Primer,另一條為amiR167-PVY的成熟序列(轉(zhuǎn)基因煙草)或煙草天然miR167成熟序列(野生型煙草),內(nèi)參基因為煙草miR159成熟序列[20]。熒光定量PCR用SYBR PrimeScript miRNA RT-PCR Kit試劑盒(TaKaRa)并按照說明書配制反應(yīng)體系,在Bio-Rad熒光定量PCR儀上擴增。反應(yīng)條件為:95℃預(yù)變性10 s;95℃變性5 s,60℃退火20 s,72℃延伸20 s,40個循環(huán);72℃延伸5 min。數(shù)據(jù)讀取由實時熒光定量PCR儀自動完成,供試基因的擴增效率均在85%以上?;虮磉_水平計算采用2-ΔΔCt方法。試驗進行3次技術(shù)重復(fù),采用t檢驗進行顯著性差異分析。

    1.2.6農(nóng)桿菌介導(dǎo)轉(zhuǎn)化煙草獲得轉(zhuǎn)基因植株 將構(gòu)建正確的植物雙元表達載體轉(zhuǎn)化到農(nóng)桿菌GV3101中,挑取單菌落搖菌,經(jīng)PCR驗證正確后以1∶100的比例擴大培養(yǎng)。參照Horsc等[21]的方法進行農(nóng)桿菌介導(dǎo)轉(zhuǎn)化煙草,獲得轉(zhuǎn)基因植株。

    1.2.7病毒侵染 將保存的帶有病毒毒源的煙草葉片研磨,加入PBS緩沖液稀釋(1 g∶10 mL),對煙草的倒數(shù)第2、3片葉進行摩擦接毒,取樣部位為侵染葉的非接毒區(qū)和非侵染葉(倒數(shù)第4片葉)。在病毒侵染前期,我們以葉片中病毒外殼蛋白基因(Coat protein,CP)和移動蛋白基因(Move protein,MP)的表達水平判斷病毒含量;在侵染后期,以煙草植株上的病毒癥狀判斷染病情況。

    表1 用于構(gòu)建質(zhì)粒和定量分析的引物

    2 結(jié)果

    2.1煙草NtWRKY40對病毒侵染的響應(yīng)

    2.1.1NtWRKY40的生物信息學(xué)分析 克隆獲得的受病毒誘導(dǎo)的煙草WRKY基因cDNA長777 bp,編碼258個氨基酸,在NCBI比對相似度最高的為絨毛狀煙草NtoWRKY40(XM_009612054.1)(95%),其次是美花煙草NsWRKY40(XM_009782127.1)(93.0%),番茄SlWRKY40(NM_001317914.1)(83.0%),故將該基因命名為NtWRKY40,GenBank登錄號為AB024510.1。

    參照模式植物擬南芥WRKY家族的分類、進化關(guān)系,NtWRKY40屬于WRKY 家族基因Ⅱa 亞族成員(圖1-A),多重序列比對顯示NtWRKY40具有該家族完整的結(jié)構(gòu),包括一個典型的WRKY結(jié)構(gòu)域和一個假定的鋅指結(jié)構(gòu)(C-X4-5-C-X22-23-H-X1-H,圖1-B)。在線軟件WoLF PSORT預(yù)測NtWRKY40定位在細胞核中,序列分析表明,NtWRKY40含有3個可能的核定位信號(82RKRK85,104PKRPREI110,163PVKKKVQ169,圖1-B)。

    2.1.2NtWRKY40過表達轉(zhuǎn)基因煙草植株的獲得 為了進一步研究煙草NtWRKY40在應(yīng)答病毒侵染中的作用,通過雙酶切的方法將目的基因插入植物表達載體pBI121,并置于花椰菜花葉病毒35S啟動子(35S)之后。再通過農(nóng)桿菌GV3101介導(dǎo),獲得35S∷NtWRKY40轉(zhuǎn)基因煙草植株。經(jīng)過PCR鑒定,獲得7個不同的轉(zhuǎn)基因株系(W-1-W-7)。NtWRKY40在不同獨立的轉(zhuǎn)基因株系中的表達水平與野生型煙草相比,有不同程度的提高(圖2-A)。

    2.1.3NtWRKY40過表達植株對植物病毒的抗病分析 選擇NtWRKY40表達水平較高(W-5)和較低(W-3)的兩個株系各20株接種TMV和PVY兩種病毒,野生型煙草(WT)作為對照。觀察發(fā)現(xiàn),過表達植株比野生型植株早2 d出現(xiàn)病毒癥狀,在第9天時,轉(zhuǎn)基因煙草的斑駁花葉病癥較野生型更為明顯(圖3),表明NtWRKY40過表達株系對病毒的敏感性增加。TMV病毒侵染實驗(圖2-B)表明,野生型煙草的NtWRKY40表達在侵染后6 h和3 d均有升高,但兩者沒有明顯差異。但在侵染后9 d,NtWRKY40表達量顯著升高至148倍。NtWRKY40過表達轉(zhuǎn)基因植株W-3在病毒侵染6 h和3 d時,NtWRKY40的表達量均比野生型高。但在侵染后的9 d時,其表達水平則稍低于野生型。W-5轉(zhuǎn)基因植株在病毒侵染后6 h和3 d時,NtWRKY40的表達水平遠高于野生型,而9 d時雖高于野生型(W-5為184,WT為148),但二者沒有顯著性差異。從轉(zhuǎn)基因植株來看,9 d時NtWRKY40的表達水平均比6 h和3 d略有降低。

    圖1 NtWRKY40基因系統(tǒng)進化樹(A)和氨基酸序列分析(B)

    圖2 35S∷NtWRKY40過表達煙草中NtWRKY40的表達(A)及病毒侵染后的變化(B)

    為了進一步研究病毒侵染后病毒含量的動態(tài)變化,以確定NtWRKY40基因與感病程度的關(guān)系,利用定量PCR的方法分析了NtWRKY40過表達煙草與野生煙草的非侵染葉中TMV病毒外殼蛋白基因(coat protein,CP)和移動蛋白基因(move protein,MP)的表達水平。在接種TMV病毒后6 h,即可檢測到CP和MP的表達(圖4)。分別以6 h時野生型煙草中CP和MP的表達水平為1,計算此后不同時間點的病毒基因表達的相對量。接種病毒后3 d(圖4-A),野生植株CP和MP的表達水平比6 h時高了1-2倍,而W-3和W-5中CP和MP升高了4-5倍;接種病毒后9 d(圖4-B),野生植株CP和MP的表達水平升高了5 000多倍,轉(zhuǎn)基因植株W-3升高了1×104多倍,而W-5的表達水平則升高了近1.5× 105倍。

    圖3 TMV病毒侵染后9 d野生型煙草(WT)與NtWRKY40轉(zhuǎn)基因煙草株系W-3和W-5病癥表型

    綜合圖2和圖4的結(jié)果,NtWRKY40的表達水平與CP和MP的表達水平呈正相關(guān),且與葉片的感病程度一致,說明過表達NtWRKY40使植物對病毒的侵染更加敏感,表明NtWRKY40是一個負調(diào)控因子。

    圖4 病毒外殼蛋白基因(CP)和移動蛋白基因(MP)在TMV病毒侵染后3 d(A)和9 d(B)的表達變化

    2.2amiRNA抗病毒技術(shù)

    2.2.1amiR167-PVY轉(zhuǎn)基因煙草植株的獲得 為進一步探索應(yīng)用人工miRNA抑制病毒病害的技術(shù),以煙草miR167e為骨架,以PVY病毒外殼蛋白為抑制對象,應(yīng)用WMD3程序獲得反義基因序列,進而構(gòu)建出35S∷amiR167-PVY重組質(zhì)粒,通過農(nóng)桿菌介導(dǎo)法建立35S∷amiR167-PVY轉(zhuǎn)基因植株的3個株系PVY-1、PVY-2和PVY-3。在35S啟動子的驅(qū)動下,3個株系的amiR167-PVY成熟序列均比天然miR167有上百倍的表達(圖5-A)。以野生型煙草中miR167前體序列的表達水平為對照1,以此計算轉(zhuǎn)基因植株中amiR167-PVY前體序列(preamiR167-PVY)的轉(zhuǎn)錄水平,結(jié)果(圖5-B)顯示,3個株系的前體序列的表達水平也有明顯的提高,且成熟序列的表達水平與前體基因的表達水平呈正相關(guān)。

    圖5 35S∷amiR167-PVY轉(zhuǎn)基因煙草植株中amiR167-PVY成熟序列(A)和前體(B)的定量分析

    2.2.2amiR167-PVY轉(zhuǎn)基因煙草植株抗病毒實驗 將35S∷amiR167-PVY轉(zhuǎn)基因植株的3個株系擴繁后接種PVY病毒,侵染30 d后對野生型煙草及PVY-1、PVY-2和PVY-3三個株系的抗病能力進行了統(tǒng)計。結(jié)果(表2)顯示,amiR167-PVY的3個株系均出現(xiàn)了病毒癥狀,但與野生型相比,轉(zhuǎn)基因植株發(fā)病率較野生型低,病癥發(fā)展較慢,對病毒有一定的抑制作用。

    表2 野生型煙草和轉(zhuǎn)基因煙草的抗病能力比較

    圖6 病毒侵染6 h和3 d后病毒外殼蛋白PVY-CP基因的定量分析

    圖7 野生煙草(CK)和過表達amiR167-PVY的轉(zhuǎn)基因煙草(3-9,3-1,3-8)病毒侵染1個月后的病癥表型

    為了追蹤病毒在侵染過程中的動態(tài)變化,對PVY-3號株系中的10棵煙草和野生型煙草中PVY病毒外殼蛋白基因(PVY-CP)進行定量分析。以侵染后6 h野生型煙草PVY-CP的表達水平為1,計算侵染后6 h和3 d轉(zhuǎn)基因植株非侵染葉中PVY-CP的轉(zhuǎn)錄水平,結(jié)果(圖6)顯示,6 h時除3-9中的表達水平與野生型相似外,其他植株的表達較低或者檢測不到,3 d時野生型植株的表達升高至6 h的3倍,而轉(zhuǎn)基因植株的表達顯著低于野生型,對病毒具有抑制效果。

    病毒侵染1個月后,總體上轉(zhuǎn)基因煙株與野生株相比,病癥顯著減輕,但仍有個別轉(zhuǎn)基因植株的病癥與野生株相近(圖7)。為了進一步了解病癥與人工amiR167-PVY、病毒外殼蛋白基因表達的關(guān)系,選擇amiR167-PVY-3轉(zhuǎn)基因植株F2分離后代中病癥較重的3-9、病癥較輕的3-1和抗病的3-8(圖7),進行上述基因的定量分析。

    以野生型煙草中miR167成熟序列的表達水平為1,計算轉(zhuǎn)基因植株中amiR167-PVY成熟序列的轉(zhuǎn)錄水平。結(jié)果(圖8-A)顯示,病癥較輕植株的表達水平高于病癥重者(3-8>3-1>3-9)。以侵染后6 h的野生型煙草中PVY-CP的表達水平為1,計算轉(zhuǎn)基因植株病葉和非病葉中PVY-CP的轉(zhuǎn)錄水平,結(jié)果(圖8-B)顯示,與6 h的野生型煙草相比,染病植株中PVY-CP的表達水平均有不同程度的升高,且與病癥的嚴重程度相關(guān),即病癥嚴重的植株P(guān)VYCP表達水平高,病癥較輕的則表達水平低,而在未染病植株3-8中始終未檢測到病毒相關(guān)基因的表達。在轉(zhuǎn)基因植株中,PVY-CP的表達水平與成熟序列的表達水平表現(xiàn)出負相關(guān),說明amiRAN抑制了病毒基因的表達,起到了抗病毒作用。

    圖8 侵染病毒一個月后染病植株中成熟序列(A)和病毒外殼蛋白(B)的定量分析

    3 討論

    3.1WRKY家族基因在病毒防御信號中的作用

    對影響植物防衛(wèi)反應(yīng)的一些WRKY家族基因的研究表明,當植物受到不同病原體或是防御信號物質(zhì)誘導(dǎo)時,WRKY基因的轉(zhuǎn)錄、蛋白合成以及結(jié)合活性都會產(chǎn)生明顯變化,從而對不同的抗性反應(yīng)作出調(diào)節(jié)[5,22]。在病原菌感染以及用病原菌誘導(dǎo)子或水楊酸(salicylic acid,SA)處理時,許多植物中的WRKY基因迅速被誘導(dǎo)[23,24],隨后WRKY蛋白通過特異性的識別啟動子上的W-box,激發(fā)一系列與抗病相關(guān)基因的表達,例如,PR基因等[25,26]。一些研究對特定的WRKY蛋白參與植物防御反應(yīng)提供了直接證據(jù),例如,西芹中WRKY1作為轉(zhuǎn)錄激活子,通過與W-box結(jié)合來調(diào)節(jié)真菌誘導(dǎo)子誘導(dǎo)的基因表達[27]。有研究報道,植物對機械損傷和病原菌侵染的響應(yīng)有著共同的信號傳導(dǎo)元件,都能夠誘導(dǎo)WRKY轉(zhuǎn)錄因子的表達[5]。當前對WRKY家族的研究主要集中該轉(zhuǎn)錄因子調(diào)節(jié)防衛(wèi)反應(yīng)中下游抗病基因的表達機制,以及確立它所參與的信號轉(zhuǎn)導(dǎo)途徑中一些上游信號傳遞元件[28]。在擬南芥中的研究發(fā)現(xiàn),當病原菌侵染擬南芥時,AtWRKY40主要通過調(diào)控JZS的表達進而影響一些相關(guān)信號途徑,如對JA信號途徑的影響使得擬南芥對霉粉菌的入侵更加敏感[5]。在煙草中過表達NtWRKY-R1基因,導(dǎo)致JA信號途徑關(guān)鍵防御基因PDF1.2表達量降低[29],說明NtWRKY-R1 基因影響JA 信號途徑,參與逆境信號調(diào)節(jié)作用。

    目前利用農(nóng)桿菌介導(dǎo)的轉(zhuǎn)基因技術(shù)并結(jié)合植物的逆境處理,可以有效揭示出WRKY轉(zhuǎn)錄因子所參與和調(diào)控的信號傳導(dǎo)途徑[5,22]。本研究中所用的35S啟動子是組成型啟動子,其本身不受病毒或機械損傷的影響,而植物內(nèi)源的NtWRKY40受病毒誘導(dǎo),從病毒侵染后的發(fā)病癥狀來看,NtWRKY40過表達煙草對病毒的反應(yīng)更加敏感,感病時效更短,病毒的表達與NtWRKY40的表達呈現(xiàn)出正相關(guān),因此我們認為NtWRKY40是一個負調(diào)控因子,這有可能是病毒侵染6 h和3 d時,過表達株系中NtWRKY40的表達水平顯著高于野生型的主要原因,而野生煙草在病毒侵染9 d時才受到顯著誘導(dǎo)(圖2-B)。我們在取樣過程中對植株造成一定的機械損傷,可能引起了PRs和JZS基因的表達變化,使得NtWRKY40過表達植株對TMV的敏感性增加。另外,有研究顯示[13],煙草WRKY的表達水平與葉片發(fā)育程度呈負相關(guān),即幼葉中表達水平高,老葉中表達水平低,本研究中,NtWRKY40的轉(zhuǎn)基因煙草在病毒侵染后9 d,NtWRKY40的表達水平低于6 h和3 d,可能是因為轉(zhuǎn)基因煙草對病毒的侵染更加敏感,導(dǎo)致其衰老程度加快,因而9 d時NtWRKY40的表達水平下降。

    3.2amiRNA骨架和病毒靶基因位點的選擇

    到目前為止,一般amiRNA均以天然miRNA作為骨架,一方面可以避免amiRNA轉(zhuǎn)錄后的降解碎片可能引起的內(nèi)源基因的表達受到意外抑制,另一方面,可以使amiRNA在生物體內(nèi)得以像天然miRNA一樣順利產(chǎn)生成熟的短鏈amiRNA片段[30,31]。不同的miRNA骨架對靶基因的抑制效果也不盡相同,Ai等[32]在擬南芥上的研究發(fā)現(xiàn),在選用相同的靶基因和成熟序列的前體下,miR159a,miR167b和miR171a三種不同前體對PVX(Potato virus X,馬鈴薯X病毒)和PVY病毒的抑制效果有所不同,以miR159a為骨架對病毒的抑制效率達到100%,而miR167b和miR171a則分別為75%和50%。這可能是由于植物在產(chǎn)生amiRNA時,miR159前體(pre-miR159a)的莖環(huán)結(jié)構(gòu)更容易被DCL1(Dicer-like protein 1)識別,從而更有效的產(chǎn)生miRNA:miRNA*復(fù)合體(*代表反義鏈,與miRNA成熟序列互補);也可能與miRNA前體的自由能有關(guān)(pre-miR159a △G=-73 kcal·mol-1,premiR167b △G=-46.9 kcal·mol-1,pre-miR171a △G= -54.2 kcal·mol-1),自由能越低莖環(huán)結(jié)構(gòu)越穩(wěn)定,從而前體轉(zhuǎn)錄折疊成二級結(jié)構(gòu)的時間變短,加快形成RISC復(fù)合體(RNA-induced silencing complex),實現(xiàn)更有效的基因沉默[32]。本研究以煙草miR167e為骨架構(gòu)建amiRNA,根據(jù)在線網(wǎng)站RNAfold WebServer(http://rna.tbi.univie.ac.at/)預(yù)測,煙草miR167e的自由能為-31.35 kcal·mol-1,自由能偏高可能是amiR167-PVY不能完全抑制PVY病毒的原因之一。

    此外,為數(shù)不多的有關(guān)amiRNA中人工靶基因片段序列選擇的研究表明,選擇高效率的小RNAs(small RNAs)和具有高度專一性的靶基因是大多數(shù)在線amiRNA序列分析軟件設(shè)計amiRNA的原則,如WMD3[15]和PsRNAtarget[33](http://plantgrn. noble.org/psRNATarget/)等。本研究中選用的amiRNA成熟序列便是基于WMD3網(wǎng)站預(yù)測獲得,選取供試序列中的最佳序列。有文獻報道,WMD3網(wǎng)站平臺是基于序列的互補配對和amiRNA與靶基因的結(jié)合能量進行預(yù)測,為每個靶基因提供出上百條候選amiRNA序列,而這些amiRNA在生物體內(nèi)的表達、加工機器成熟序列的產(chǎn)生尚難預(yù)測[34]。amiRNA與靶基因的互作受到生物體多種因素的影響,包括靶基因mRNA的序列和mRNA結(jié)合蛋白等[35,36]。因此,用于基因沉默的最佳amiRNA僅從目標反向互補的序列上判斷,仍有不足,尚需進一步的研究。

    本研究及此前的研究均表明,RNA病毒如TMV、PVY的侵染能夠極顯著地誘導(dǎo)煙草NtWRKY40的表達,而該基因在干旱、低溫等逆境條件下的表達卻無顯著變化[13]。這些結(jié)果表明,該基因的表達控制元件對RNA病毒有較為專一的反應(yīng)。因此,利用這一特點可以構(gòu)建以抗病毒人工amiRNA和NtWRKY40表達控制元件為主體的表達系統(tǒng),有可能實現(xiàn)對病毒侵染自動應(yīng)答式的控制。

    4 結(jié)論

    煙草NtWRKY40是WRKY家族基因Ⅱa 亞族成員,與擬南芥AtWRKY40基因同源,具有該家族的典型結(jié)構(gòu)域,它的過量表達使植物對于病毒的侵染更加敏感,是抗病負調(diào)控因子,同時它的過表達可能刺激了一些與植物生長有關(guān)的基因的表達,我們推斷NtWRKY40可能是參與機械損傷脅迫和病原菌侵染的重要組分。

    另外,以馬鈴薯Y病毒基因特定反相互補片段序列構(gòu)建的amiR167-PVY在轉(zhuǎn)基因煙草的抗病毒試驗中,表現(xiàn)出對病毒有一定的抑制作用。

    [1]Chen W, HuangT, Dai J, et al. Evaluations of tobacco cultivars resistance to Tobacco mosaic virus and Potato virus Y[J]. Plant Pathology Journal, 2014, 13(1):37-43.

    [2]Yan DH, Fenning T, Tang S, et al. Genome-wide transcriptional response of Populus euphratica to long-term drought stress[J]. Plant Sci, 2012, 195:24-35.

    [3]Zhao MY, Zhang ZB, Chen SY, et al. WRKY transcription factor superfamily:Structure, origin and functions[J]. African Journal of Biotechnology, 2012, 11(32):8051-8059.

    [4]Schmutz J, Cannon SB, Schlueter J, et al. Genome sequence of the palaeopolyploid soybean[J]. Nature, 2010, 463(7278):178-183.

    [5] Pandey SP, Roccaro M, Sch?n M, et al. Transcriptional reprogramming regulated by WRKY18 and WRKY40 facilitates powdery mildew infection of Arabidopsis[J]. Plant J, 2010, 64(6):912-923.

    [6]Naoumkina M, He X, Dixon R. Elicitor-induced transcription factors for metabolic reprogramming of secondary metabolism in Medicago truncatula[J]. BMC Plant Biol, 2008, 8:132.

    [7] Popescu SC, Popescu GV, Bachan S, et al. MAPK target networks in Arabidopsis thaliana revealed using functional protein microarrays[J]. Genes Dev, 2009, 23(1):80-92.

    [8]Ishihama N, Yamada R, Yoshioka M, et al. Phosphorylation of the Nicotiana benthamiana WRKY8 transcription factor by MAPK functions in the defense response[J]. Plant Cell, 2011, 23(3):1153-1170.

    [9] Skibbe M, Qu N, Galis I, et al. Induced plant defenses in the natural environment:Nicotiana attenuata WRKY3 and WRKY6 coordinate responses to herbivory[J]. Plant Cell, 2008, 20(7):1984-2000.

    [10]Nishiuchi T, Shinshi H, Suzuki K. Rapid and transient activation of transcription of the ERF3 gene by wounding in tobacco leaves -Possible involvement of NtWRKYs and autorepression[J]. J Biol Chem, 2004, 279(53):55355-55361.

    [11]Yamamoto S, Nakano T, Suzuki K, Shinshi H. Elicitorinduced activation of transcription via W box-related cis-acting elements from a basic chitinase gene by WRKY transcription factors in tobacco[J]. Biochim Biophys Acta, 2004, 1679(3):279-287.

    [12]Adachi H, Nakano T, Miyagawa N, et al. WRKY transcription factors phosphorylated by MAPK regulate a plant immune NADPH oxidase in Nicotiana benthamiana[J]. Plant Cell, 2015, 27(9):2645-2663.

    [13] 周漢琛, 席玉珍, 魏書. RNA病毒誘導(dǎo)的煙草WRKY轉(zhuǎn)錄因子基因的應(yīng)急表達[J]. 安徽農(nóng)業(yè)大學(xué)學(xué)報, 2014, 41(6):1020-1026.

    [14]Carbonell A, Fahlgren N, Mitchell S, et al. Highly specific gene silencing in a monocot species by artificial microRNAs derived from chimeric miRNA precursors[J]. Plant J, 2015, 82(6):1061-1075.

    [15] Niu QW, Lin SS, Reyes JL, et al. Expression of artificial micro-RNAs in transgenic Arabidopsis thaliana confers virus resistance[J]. Nat Biotechnol, 2008, 24(11):1420-1428.

    [16]Qu J, Ye J, Fang R. Artificial microRNA-mediated virus resistance in plants[J]. J Virol, 2007, 81(12):6690-6699.

    [17]Tang S, Wang Y, Li Z, et al. Identification of wounding and topping responsive small RNAs in tobacco(Nicotiana tabacum)[J]. BMC Plant Biol, 2012, 12:28.

    [18]Ossowski S, Schwab R, Weigel D. Gene silencing in plants using artificial microRNAs and other small RNAs[J]. Plant J, 2008,53(4):674-690.

    [19] Gregor WS, Sven KD. Stable internal reference genes for normalization of real-time RT-PCR in tobacco(Nicotiana tabacum)during development and abiotic stress[J]. Mol Genet Genomics,2010, 283(3):233-241.

    [20]Lang QL, Jin CZ, Lai LY, et al. Tobacco microRNAs prediction and their expression infected with Cucumber mosaic virus and Potato virus X[J]. Mol Biol Rep, 2011, 38(3):1523-1531.

    [21]Horsch RB, Fry JE, Hoffmann NL, et al. A simple and general method for transferring genes into plant[J]. Science, 1985, 227(4691):1229-1231.

    [22]Abbruscato P, Nepusz T, Mizzd L, et al. OsWRKY22, a monocot WRKY gene, plays a role in the resistance response to blast[J]. Mol Plant Pathol, 2012, 13(8):828-841.

    [23]Pan LJ, Jiang L. Identification and expression of the WRKY transcription factors of Carica papaya in response to abiotic andbiotic stresses[J]. Mol Biol Rep, 2014, 41(3):1215-25.

    [24]Rushton PJ, Torres JT, Parniske M, et al. Interaction of elicitorinduced DNA-binding proteins with elicitor response elements in the promoters of parsley PR1 genes[J]. EMBO J, 1996, 15(20):5690-700.

    [25]Chen C, Chen Z. Isolation and characterization of two pathogenand salicylic acid-induced genes encoding WRKY DNA-binding proteins from tobacco[J]. Plant Mol Biol, 2000, 42(2):387-96.

    [26]Dong J, Chen C, Chen Z. Expression profiles of the Arabidopsis WRKY gene superfamily during plant defense response[J]. Plant Mol Biol, 2003, 51(1):21-37.

    [27]Turck F, Zhou A, Somssich IE. Stimulus-dependent, promoterspecific binding of transcription factor WRKY1 to Its native promoter and the defense-related gene PcPR1-1 in Parsley[J]. Plant Cell, 2004, 16(10):2573-2585.

    [28]Chi Y, Yang Y, Zhou Y, et al. Protein-protein interactions in the regulation of WRKY transcription factors[J]. Mol Plant, 2013, 6(2):287-300.

    [29]胡曼, 宋丹丹, 楊金淼, 等. 煙草WRKY-R1基因的克隆及瞬時表達分析[J]. 山東農(nóng)業(yè)科學(xué), 2014, 46(9):12-16.

    [30]Tiwari M, Sharma D, Trivedi PK. Artificial microRNA mediated gene silencing in plants:progress and perspectives[J]. Plant Mol Biol, 2014, 86(1-2):1-18.

    [31]Sablok G, Pérez-Quintero AL, Hassan M, et al. Artificial microRNAs(amiRNAs)engineering - On how microRNA-based silencing methods have affected current plant silencing research[J]. Biochem Biophys Res Commun, 2011, 406(3):315-319.

    [32]Ai T, Zhang L, Gao Z, et al. Highly efficient virus resistance mediated by artificial microRNAs that target the suppressor of PVX and PVY in plants[J]. Plant Biol(Stuttg), 2011, 13(2):304-316.

    [33]Dai X, Zhao PX. psRNATarget:a plant small RNA target analysis server[J]. Nucleic Acids Res, 2011, 39(Web Server issue):W155-9.

    [34]Li J, Chung HS, Niu Y, et al. Comprehensive protein-based artificial microRNA screens for effective gene silencing in plants[J]. Plant Cell, 2013, 25(5):1507-1522.

    [35]Fabian MR, Sonenberg N, Filipowicz W. Regulation of mRNA translation and stability by microRNAs[J]. Annu Rev Biochem,2010, 79(1):351-379.

    [36]Pasquinelli AE. MicroRNAs and their targets:Recognition,regulation and an emerging reciprocal relationship[J]. Nat Rev Genet, 2012, 13(4):271-28.

    (責(zé)任編輯李楠)

    Roles of Nicotiana tobacum NtWRKY40 in Plant Responding to Virus Infection

    LIU Jing-jing CHENG Chun-ling XI Yu-zhen WEI Shu
    (State Key Laboratory of Tea Plant Biology and Utilization,Anhui Agricultural University,Hefei 230036)

    Some of WRKY transcription factors have involved in plant abiotic defense. A WRKY gene,a homologue of Arabidopsis WRKY40,was found being dramatically induced by plant virus and further characterized in this study. To verify its function,NtWRKY40 overexpressed tobacco plants were established. After virus infection,the transcription levels of NtWRKY40 and virus related genes were quantified using quantitative real time PCR at three stages(6 h,3 d,9 d after infection). Results showed that tobacco mosaic virus infection led to 150 fold higher of NtWRKY40 transcript in non-transgenic wild type tobacco than non-infected control. No significant change in NtWRKY40 transcription levels was noted in the NtWRKY40 over-expressed tobacco after virus infection. NtWRKY40 over-expressed tobaccos were more sensitive to TMV infection than non-transgenic wild type. The data suggested that NtWRKY40 was a negative regulator in responding to TMV. In addition,an artificial miRNA vector amiR167-PVY was constructed,for silencing the coat protein gene of Potato virus Y,using the miR167 precursor of Nicotiana tabacum as backbone and cauliflower mosaic virus 35S promoter. Over-expression of amiR167-PVY exhibited inhibitory effect on the virus. The results provide the basis for the further research on amiRNA-mediated virus resistance mechanism,and contributes to antiviral transgenic plants cultivation with the optimal target sequence.

    Tobacco mosaic virus;Potato virus Y;WRKY transcription factor;artificial microRNA;virus resistance

    2016-06-01

    國家自然科學(xué)基金項目(31070614,31370687),教育部博士學(xué)科點專項科研基金項目(20123418110002)

    劉晶晶,女,博士研究生,研究方向:生物化學(xué)與分子生物學(xué),E-mail:xiaoayatou@163.com;程春玲同為本文第一作者

    魏書,男,教授,博士生導(dǎo)師,研究方向:茶樹次生代謝;E-mail:weishu@ahau.edu.cn

    猜你喜歡
    前體侵染轉(zhuǎn)基因
    探秘轉(zhuǎn)基因
    轉(zhuǎn)基因,你吃了嗎?
    揭示水霉菌繁殖和侵染過程
    N-末端腦鈉肽前體與糖尿病及糖尿病相關(guān)并發(fā)癥呈負相關(guān)
    N-端腦鈉肽前體測定在高血壓疾病中的應(yīng)用研究
    蕓薹根腫菌侵染過程及影響因子研究
    甘藍根腫病菌休眠孢子的生物學(xué)特性及侵染寄主的顯微觀察
    天然的轉(zhuǎn)基因天然的轉(zhuǎn)基因“工程師”及其對轉(zhuǎn)基因食品的意蘊
    茶葉香氣前體物研究進展
    茶葉通訊(2014年2期)2014-02-27 07:55:40
    煙草靶斑病(Thanatephorus cucumeris)侵染特性研究
    真人一进一出gif抽搐免费| 人人妻,人人澡人人爽秒播| 伦理电影免费视频| 大型av网站在线播放| 亚洲国产精品999在线| 欧美精品亚洲一区二区| 亚洲欧美精品综合一区二区三区| 国产亚洲精品综合一区在线观看 | 999精品在线视频| 侵犯人妻中文字幕一二三四区| 亚洲中文日韩欧美视频| 91麻豆精品激情在线观看国产| 大型黄色视频在线免费观看| 亚洲国产欧美网| 国产又黄又爽又无遮挡在线| 在线观看免费日韩欧美大片| 18禁国产床啪视频网站| 12—13女人毛片做爰片一| 亚洲国产精品sss在线观看| 国产三级黄色录像| 久久久久九九精品影院| 国产色视频综合| 午夜福利一区二区在线看| 国产精品免费视频内射| 久久久久久久久中文| 欧美大码av| 99精品久久久久人妻精品| 宅男免费午夜| 一a级毛片在线观看| 欧美日韩瑟瑟在线播放| 国产激情欧美一区二区| √禁漫天堂资源中文www| 亚洲第一欧美日韩一区二区三区| 老司机午夜福利在线观看视频| 久久久久国内视频| 久久中文字幕一级| 国产91精品成人一区二区三区| 伊人久久大香线蕉亚洲五| 成年版毛片免费区| 看黄色毛片网站| 欧美绝顶高潮抽搐喷水| 欧美+亚洲+日韩+国产| 伊人久久大香线蕉亚洲五| 一区二区三区高清视频在线| 非洲黑人性xxxx精品又粗又长| 亚洲狠狠婷婷综合久久图片| 免费一级毛片在线播放高清视频| 亚洲aⅴ乱码一区二区在线播放 | 亚洲在线自拍视频| 丝袜美腿诱惑在线| 久久精品夜夜夜夜夜久久蜜豆 | 99热这里只有精品一区 | 欧美成狂野欧美在线观看| 丁香欧美五月| 好看av亚洲va欧美ⅴa在| 久久这里只有精品19| 香蕉丝袜av| 黄频高清免费视频| 久久亚洲真实| 黄频高清免费视频| 他把我摸到了高潮在线观看| 国产爱豆传媒在线观看 | 波多野结衣av一区二区av| 国产亚洲精品av在线| svipshipincom国产片| 午夜福利视频1000在线观看| 两个人免费观看高清视频| 在线观看www视频免费| 九色国产91popny在线| 欧美黑人巨大hd| 色综合亚洲欧美另类图片| 亚洲熟妇中文字幕五十中出| 国产精品久久久人人做人人爽| 国产精品野战在线观看| 亚洲欧美激情综合另类| 大型av网站在线播放| av欧美777| 亚洲精品一卡2卡三卡4卡5卡| 天天一区二区日本电影三级| 国产亚洲精品av在线| 亚洲色图 男人天堂 中文字幕| 亚洲精品在线美女| 亚洲欧洲精品一区二区精品久久久| 两个人视频免费观看高清| 丝袜美腿诱惑在线| or卡值多少钱| 日本免费一区二区三区高清不卡| www日本在线高清视频| 一区二区日韩欧美中文字幕| 男男h啪啪无遮挡| 麻豆国产av国片精品| 国产精品免费一区二区三区在线| 国产野战对白在线观看| 日本 欧美在线| 亚洲国产欧美日韩在线播放| 国产亚洲精品久久久久久毛片| 欧美日韩瑟瑟在线播放| 少妇 在线观看| 麻豆成人av在线观看| 欧美乱码精品一区二区三区| 色尼玛亚洲综合影院| 国产日本99.免费观看| www.www免费av| 99热只有精品国产| 免费av毛片视频| 国产欧美日韩一区二区三| 久久精品亚洲精品国产色婷小说| 中文字幕精品亚洲无线码一区 | 一a级毛片在线观看| 亚洲成人免费电影在线观看| 成人永久免费在线观看视频| 亚洲欧美日韩高清在线视频| 狠狠狠狠99中文字幕| 欧美黑人精品巨大| 国产一级毛片七仙女欲春2 | 婷婷亚洲欧美| 欧美精品啪啪一区二区三区| 日本三级黄在线观看| 97超级碰碰碰精品色视频在线观看| 色精品久久人妻99蜜桃| 欧美日韩亚洲综合一区二区三区_| 国产成人系列免费观看| 一级a爱视频在线免费观看| 亚洲av中文字字幕乱码综合 | 亚洲色图av天堂| 国产精品乱码一区二三区的特点| 久久性视频一级片| 一区福利在线观看| 日韩中文字幕欧美一区二区| 窝窝影院91人妻| 人妻久久中文字幕网| 久久这里只有精品19| 少妇的丰满在线观看| 国产精品久久视频播放| 亚洲第一青青草原| 999久久久国产精品视频| 美国免费a级毛片| 色综合婷婷激情| 黑人欧美特级aaaaaa片| 色综合欧美亚洲国产小说| 亚洲五月婷婷丁香| 男女床上黄色一级片免费看| 久久婷婷人人爽人人干人人爱| 成人av一区二区三区在线看| 亚洲最大成人中文| 国产精品亚洲美女久久久| netflix在线观看网站| 国内精品久久久久精免费| 国产真实乱freesex| 成人免费观看视频高清| 在线视频色国产色| 身体一侧抽搐| 亚洲国产高清在线一区二区三 | 久久久国产成人免费| 国产精品一区二区三区四区久久 | 男男h啪啪无遮挡| 欧美国产日韩亚洲一区| 色av中文字幕| 中出人妻视频一区二区| 两性夫妻黄色片| or卡值多少钱| 亚洲av成人av| 色在线成人网| 老汉色∧v一级毛片| 中文资源天堂在线| 欧美绝顶高潮抽搐喷水| 欧美日韩瑟瑟在线播放| 侵犯人妻中文字幕一二三四区| 精品久久久久久久久久免费视频| 久久人妻av系列| 欧美另类亚洲清纯唯美| 日韩有码中文字幕| 国产成人影院久久av| 夜夜爽天天搞| 亚洲aⅴ乱码一区二区在线播放 | 国产亚洲精品综合一区在线观看 | 午夜a级毛片| 亚洲成人久久性| 老熟妇仑乱视频hdxx| 熟女电影av网| 国产成人欧美| 国产亚洲av高清不卡| 亚洲欧美精品综合一区二区三区| 一二三四在线观看免费中文在| 色婷婷久久久亚洲欧美| 91在线观看av| 757午夜福利合集在线观看| 一级毛片高清免费大全| 中文资源天堂在线| 国产激情偷乱视频一区二区| 成人一区二区视频在线观看| 十分钟在线观看高清视频www| 1024视频免费在线观看| 久久香蕉国产精品| 琪琪午夜伦伦电影理论片6080| 日本a在线网址| 99riav亚洲国产免费| a在线观看视频网站| 在线观看免费午夜福利视频| 久久香蕉激情| 欧美色欧美亚洲另类二区| 欧美大码av| 欧美中文综合在线视频| 精品国产乱子伦一区二区三区| 色婷婷久久久亚洲欧美| 亚洲中文字幕日韩| 久久香蕉精品热| 激情在线观看视频在线高清| 99精品在免费线老司机午夜| 久久青草综合色| 淫妇啪啪啪对白视频| 黑人欧美特级aaaaaa片| 人人澡人人妻人| 色综合婷婷激情| 一级a爱视频在线免费观看| 脱女人内裤的视频| 久久精品国产综合久久久| 亚洲av电影不卡..在线观看| 国产成人一区二区三区免费视频网站| 在线看三级毛片| 校园春色视频在线观看| 亚洲在线自拍视频| 国产亚洲精品av在线| 国产精品久久电影中文字幕| 校园春色视频在线观看| 免费看日本二区| 在线观看日韩欧美| 真人一进一出gif抽搐免费| 成年人黄色毛片网站| 又大又爽又粗| 国产色视频综合| 两性夫妻黄色片| av片东京热男人的天堂| 精品一区二区三区四区五区乱码| 在线观看日韩欧美| 成年人黄色毛片网站| 中亚洲国语对白在线视频| 免费在线观看成人毛片| 日韩欧美在线二视频| 亚洲国产精品久久男人天堂| www.熟女人妻精品国产| 国产成人欧美在线观看| 两人在一起打扑克的视频| 久久精品国产99精品国产亚洲性色| 国产精品久久电影中文字幕| 男男h啪啪无遮挡| 国产99白浆流出| 欧美激情久久久久久爽电影| 欧美激情久久久久久爽电影| 亚洲全国av大片| 欧美乱码精品一区二区三区| 欧美性长视频在线观看| 波多野结衣巨乳人妻| 在线十欧美十亚洲十日本专区| 女同久久另类99精品国产91| 一进一出抽搐gif免费好疼| 欧美黑人巨大hd| 欧美日韩福利视频一区二区| 看免费av毛片| 日韩免费av在线播放| 琪琪午夜伦伦电影理论片6080| 啪啪无遮挡十八禁网站| 黄色丝袜av网址大全| 亚洲三区欧美一区| 亚洲激情在线av| 一本久久中文字幕| 欧美国产日韩亚洲一区| 精品国内亚洲2022精品成人| 国产精品99久久99久久久不卡| 国产99白浆流出| 亚洲av成人不卡在线观看播放网| 最好的美女福利视频网| 香蕉av资源在线| 国产欧美日韩精品亚洲av| 嫩草影院精品99| 亚洲精品中文字幕一二三四区| 国产色视频综合| 国产真实乱freesex| 一卡2卡三卡四卡精品乱码亚洲| 中文在线观看免费www的网站 | 啦啦啦免费观看视频1| 国产av在哪里看| 久久狼人影院| 亚洲第一av免费看| 视频区欧美日本亚洲| 国产精品 国内视频| 非洲黑人性xxxx精品又粗又长| 亚洲人成网站高清观看| 色老头精品视频在线观看| 男男h啪啪无遮挡| 久热爱精品视频在线9| 99久久久亚洲精品蜜臀av| 丝袜在线中文字幕| 天天躁狠狠躁夜夜躁狠狠躁| 久久精品夜夜夜夜夜久久蜜豆 | 欧美一级a爱片免费观看看 | 亚洲精华国产精华精| 亚洲成av人片免费观看| 色尼玛亚洲综合影院| 一区二区三区国产精品乱码| 天天一区二区日本电影三级| 欧美+亚洲+日韩+国产| 国产又爽黄色视频| 免费高清在线观看日韩| 国产麻豆成人av免费视频| 国产成人啪精品午夜网站| 精品卡一卡二卡四卡免费| 一级a爱视频在线免费观看| 欧美日本亚洲视频在线播放| 18禁国产床啪视频网站| 免费av毛片视频| 琪琪午夜伦伦电影理论片6080| 男女视频在线观看网站免费 | 免费高清视频大片| 成年人黄色毛片网站| 最近最新中文字幕大全电影3 | 日本免费一区二区三区高清不卡| 久久国产精品男人的天堂亚洲| 欧美日韩亚洲综合一区二区三区_| 国内精品久久久久久久电影| netflix在线观看网站| 午夜免费成人在线视频| 在线视频色国产色| 黄片小视频在线播放| 51午夜福利影视在线观看| 久久久精品欧美日韩精品| 日本免费a在线| 99在线视频只有这里精品首页| 久久这里只有精品19| 人人澡人人妻人| 手机成人av网站| 成人国产综合亚洲| 此物有八面人人有两片| 免费电影在线观看免费观看| 男女床上黄色一级片免费看| 欧美性猛交╳xxx乱大交人| 日日干狠狠操夜夜爽| xxxwww97欧美| 老司机深夜福利视频在线观看| 国产精品亚洲一级av第二区| 婷婷精品国产亚洲av在线| 亚洲av成人av| 中国美女看黄片| 色播亚洲综合网| a在线观看视频网站| 国产精品亚洲av一区麻豆| 久久 成人 亚洲| 亚洲精品在线观看二区| 少妇被粗大的猛进出69影院| 波多野结衣高清作品| 俄罗斯特黄特色一大片| 99热6这里只有精品| www国产在线视频色| 中文字幕久久专区| 亚洲一区二区三区不卡视频| 日韩免费av在线播放| 露出奶头的视频| 亚洲avbb在线观看| 国产成人系列免费观看| 午夜免费激情av| 久久 成人 亚洲| www.自偷自拍.com| 免费在线观看完整版高清| 中国美女看黄片| 两性夫妻黄色片| 看免费av毛片| 婷婷六月久久综合丁香| 免费看美女性在线毛片视频| 精品一区二区三区四区五区乱码| 可以在线观看毛片的网站| 亚洲一区二区三区色噜噜| 中亚洲国语对白在线视频| 无限看片的www在线观看| 国产欧美日韩精品亚洲av| 国产国语露脸激情在线看| 国产亚洲精品综合一区在线观看 | 男女做爰动态图高潮gif福利片| а√天堂www在线а√下载| 国产亚洲欧美98| 欧美最黄视频在线播放免费| 中文亚洲av片在线观看爽| 免费av毛片视频| 国产精品亚洲av一区麻豆| 久久婷婷成人综合色麻豆| 久久久国产欧美日韩av| 极品教师在线免费播放| 亚洲成人久久爱视频| 午夜久久久久精精品| 欧美日韩瑟瑟在线播放| 久热这里只有精品99| 最新在线观看一区二区三区| 成人国产一区最新在线观看| 1024手机看黄色片| 久热爱精品视频在线9| 哪里可以看免费的av片| 好男人在线观看高清免费视频 | 国产激情久久老熟女| 激情在线观看视频在线高清| 久久人人精品亚洲av| 精品乱码久久久久久99久播| 久久天躁狠狠躁夜夜2o2o| 国产v大片淫在线免费观看| 欧美日韩瑟瑟在线播放| 99re在线观看精品视频| 欧美激情高清一区二区三区| 中文字幕精品亚洲无线码一区 | 久久天躁狠狠躁夜夜2o2o| 国产亚洲欧美98| 国产精品野战在线观看| 日韩欧美 国产精品| 听说在线观看完整版免费高清| 国产人伦9x9x在线观看| avwww免费| 精品一区二区三区四区五区乱码| 老司机靠b影院| 欧美日韩亚洲综合一区二区三区_| 怎么达到女性高潮| av免费在线观看网站| av在线播放免费不卡| 琪琪午夜伦伦电影理论片6080| 无人区码免费观看不卡| 精品卡一卡二卡四卡免费| 亚洲欧洲精品一区二区精品久久久| 最近最新中文字幕大全电影3 | 淫秽高清视频在线观看| 亚洲av成人一区二区三| 免费看美女性在线毛片视频| 久久久久国产精品人妻aⅴ院| 久久亚洲真实| 在线观看www视频免费| 国产色视频综合| 一区二区三区高清视频在线| 国产视频内射| 国产精品香港三级国产av潘金莲| 亚洲欧洲精品一区二区精品久久久| 久久精品aⅴ一区二区三区四区| x7x7x7水蜜桃| 欧美亚洲日本最大视频资源| 久久中文字幕一级| 级片在线观看| 亚洲人成网站在线播放欧美日韩| 亚洲全国av大片| 色播亚洲综合网| 国产亚洲精品av在线| 欧美激情高清一区二区三区| 中文字幕av电影在线播放| 欧美激情 高清一区二区三区| 久久久久久亚洲精品国产蜜桃av| 亚洲一区中文字幕在线| 欧美乱码精品一区二区三区| 在线观看www视频免费| 亚洲国产精品999在线| av视频在线观看入口| 夜夜躁狠狠躁天天躁| 精品久久久久久,| 久久精品成人免费网站| 国产精品亚洲美女久久久| 亚洲精品美女久久av网站| 欧美国产日韩亚洲一区| 久久精品影院6| 脱女人内裤的视频| 日本免费一区二区三区高清不卡| 亚洲精品国产一区二区精华液| 黄片播放在线免费| 日韩 欧美 亚洲 中文字幕| 国产精品美女特级片免费视频播放器 | 欧美色视频一区免费| 欧美国产精品va在线观看不卡| 听说在线观看完整版免费高清| 又黄又爽又免费观看的视频| 免费人成视频x8x8入口观看| 亚洲专区字幕在线| АⅤ资源中文在线天堂| 国产成人啪精品午夜网站| 久久精品人妻少妇| 免费高清在线观看日韩| 久久伊人香网站| a在线观看视频网站| 欧美中文综合在线视频| 欧美日本亚洲视频在线播放| √禁漫天堂资源中文www| 天堂动漫精品| 久久 成人 亚洲| 琪琪午夜伦伦电影理论片6080| 亚洲人成网站高清观看| 国产又爽黄色视频| 丝袜人妻中文字幕| 欧美国产精品va在线观看不卡| 国产av不卡久久| 久久午夜综合久久蜜桃| 久久人人精品亚洲av| 真人一进一出gif抽搐免费| 9191精品国产免费久久| 国产主播在线观看一区二区| 夜夜夜夜夜久久久久| 亚洲中文字幕一区二区三区有码在线看 | 国产亚洲精品综合一区在线观看 | 久久精品国产亚洲av香蕉五月| 国产亚洲欧美98| 老司机靠b影院| 日韩大码丰满熟妇| 欧美日韩瑟瑟在线播放| 天天添夜夜摸| 久久国产精品人妻蜜桃| 两个人免费观看高清视频| 久久草成人影院| 亚洲精品一区av在线观看| 国产极品粉嫩免费观看在线| 精品不卡国产一区二区三区| 午夜亚洲福利在线播放| 午夜两性在线视频| 国内精品久久久久久久电影| 国产精品久久久人人做人人爽| 久久人妻av系列| 亚洲第一青青草原| 国产人伦9x9x在线观看| 国产亚洲精品第一综合不卡| 在线观看免费午夜福利视频| 亚洲午夜精品一区,二区,三区| 黄色视频,在线免费观看| 久久久久久久久中文| 两个人免费观看高清视频| 久久这里只有精品19| 亚洲五月色婷婷综合| 欧美黑人巨大hd| 亚洲最大成人中文| 国产精品精品国产色婷婷| 香蕉久久夜色| 一区二区三区精品91| 2021天堂中文幕一二区在线观 | 欧美激情 高清一区二区三区| 成人18禁高潮啪啪吃奶动态图| 久久中文看片网| 欧美日韩亚洲国产一区二区在线观看| 欧美日本亚洲视频在线播放| 无人区码免费观看不卡| 日韩免费av在线播放| 国产熟女午夜一区二区三区| 免费在线观看影片大全网站| 精品久久久久久久末码| 欧美绝顶高潮抽搐喷水| 午夜福利免费观看在线| 看免费av毛片| 波多野结衣高清作品| 久久热在线av| 99在线人妻在线中文字幕| www日本在线高清视频| 最近在线观看免费完整版| 巨乳人妻的诱惑在线观看| 免费在线观看日本一区| 99国产精品一区二区蜜桃av| 国语自产精品视频在线第100页| 国产精品久久视频播放| 黑人巨大精品欧美一区二区mp4| 日本 欧美在线| 一二三四在线观看免费中文在| 69av精品久久久久久| 国产成+人综合+亚洲专区| 久久久精品欧美日韩精品| netflix在线观看网站| 免费看美女性在线毛片视频| 美国免费a级毛片| 琪琪午夜伦伦电影理论片6080| 精品久久久久久久人妻蜜臀av| 国产精品免费一区二区三区在线| 中文字幕最新亚洲高清| 美女高潮喷水抽搐中文字幕| 两个人看的免费小视频| 99热只有精品国产| 国产成人av激情在线播放| 国产av不卡久久| 免费一级毛片在线播放高清视频| a级毛片a级免费在线| 在线国产一区二区在线| 欧美日韩亚洲综合一区二区三区_| 欧美乱色亚洲激情| 国产伦人伦偷精品视频| 日本熟妇午夜| 91av网站免费观看| 国产精品乱码一区二三区的特点| 国产黄色小视频在线观看| 成年版毛片免费区| 大型av网站在线播放| 亚洲色图av天堂| 亚洲精品久久国产高清桃花| 999久久久精品免费观看国产| 成年女人毛片免费观看观看9| 1024香蕉在线观看| 日韩大尺度精品在线看网址| 亚洲av日韩精品久久久久久密| 色综合婷婷激情| 免费女性裸体啪啪无遮挡网站| 亚洲自偷自拍图片 自拍| 国产黄色小视频在线观看| www日本在线高清视频| 人成视频在线观看免费观看| 国产精品98久久久久久宅男小说| 日韩大码丰满熟妇| 欧美三级亚洲精品| 久99久视频精品免费| 国产av在哪里看| 午夜久久久久精精品| 露出奶头的视频| 亚洲国产精品合色在线| 丝袜人妻中文字幕| 国产99白浆流出| 午夜老司机福利片| 亚洲精品一区av在线观看| 成年免费大片在线观看| 色在线成人网| 99热这里只有精品一区 | 亚洲,欧美精品.| 日日爽夜夜爽网站| 久久中文字幕一级| 超碰成人久久|