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

    柑橘脅迫響應基因WRKY47的克隆與表達分析

    2021-03-25 13:50:00沈丹楊莉胡威匡柳青郭文芳盧婷劉德春劉勇
    江蘇農(nóng)業(yè)學報 2021年1期
    關鍵詞:非生物脅迫表達分析基因克隆

    沈丹 楊莉 胡威 匡柳青 郭文芳 盧婷 劉德春 劉勇

    摘要: 以檸檬[Citrus limon (L.) Burm.f.]、甜橙[Citrus sinensis (L.) Osbeck]、蘆柑(Citrus reticulata Blanco)、金柑[Fortunella japonica (Thunb.) Swingle]為試驗材料,基于甜橙基因組數(shù)據(jù)庫中的甜橙基因堿基序列,利用RT-PCR方法克隆獲得4個WRKY家族基因ClWRKY47、CsWRKY47、CrWRKY47和FjWRKY47的cDNA全長堿基序列。序列分析結果表明,這4個基因的cDNA全長都是1 567 bp,開放閱讀框為1 506 bp,編碼501個氨基酸。氨基酸序列和結構分析結果顯示,這4個基因編碼的蛋白質屬于Group IIa+IIb類WRKY蛋白。進化樹分析結果顯示,所克隆的4個柑橘WRKY蛋白與克萊門柚WRKY47蛋白的親緣關系最近。對CsWRKY47啟動子順式作用元件預測分析,發(fā)現(xiàn)CsWRKY47啟動子包含脫落酸響應元件(ABRE)、茉莉酸甲酯響應元件(CGTCA-motiF)、抗氧化響應元件(ARE)、參與干旱誘導的MYB結合位點(MBS)等多個與脅迫相關的順式作用元件。實時熒光定量表達分析結果表明,檸檬ClWRKY47和甜橙CsWRKY47能被高鹽、干旱、低溫脅迫誘導表達;蘆柑CrWRKY47在干旱和低溫脅迫下誘導表達,高鹽脅迫下表達量下調;金柑FjWRKY47 在高鹽脅迫下誘導表達,干旱和低溫脅迫下下調表達。為進一步開展柑橘脅迫相關基因WRKY47的功能鑒定奠定了基礎。

    關鍵詞: 柑橘;WRKY;非生物脅迫;基因克隆;表達分析

    中圖分類號: Q785 文獻標識碼: A 文章編號: 1000-4440(2021)01-0129-10

    Cloning and expression analysis of stress response gene WRKY47 in citrus

    SHEN Dan, YANG Li, HU Wei, KUANG Liu-qing, GUO Wen-fang, LU Ting, LIU De-chun, LIU Yong

    (College of Agriculture, Jiangxi Agricultural University, Nanchang 330045, China)

    Abstract: Based on the Citrus sinensis gene base sequences contained in the sweet orange genome database, the cDNAs sequences of four WRKY family genes ClWRKY47, CsWRKY47, CrWRKY47 and FjWRKY47 were cloned from lemon [Citrus limon (L.) Burm.f.], sweet orange [Citrus sinensis (L.) Osbeck], ponkan (Citrus reticulata Blanco) and kumquat [Fortunella japonica (Thunb.) Swingle] by RT-PCR, respectively. The results of sequence analysis showed that the cDNA of these four genes were 1 567 bp in length, the open reading frame (ORF) was 1 506 bp, and 501 amino acids were encoded. The results of amino acid sequence and structure analysis demonstrated that the proteins encoded by these four genes belonged to the Group IIa+IIb WRKY protein. Phylogenetic tree analysis results indicated that the four citrus WRKY proteins were closely related to the Citrus clementine WRKY47 protein. Predictive cis-acting element analysis of CsWRKY47 promoter revealed that CsWRKY47 promoter contained abscisic acid response element (ABRE), methyljasmonic acid response element (CGTCA-motiF), anaerobic response element (ARE), MYB binding site involved in drought-response (MBS) and many other stress-related cis-acting elements. Real-time fluorescence quantitative PCR analysis suggested that the expression levels of ClWRKY47 and CsWRKY47 were induced by drought, low temperature and high salt. Besides, the expression of CrWRKY47 was induced under drought and low temperature stress, but suppressed under high salt stress. However, the expression level of FjWRKY47 was upregulated under high salt stress and suppressed under drought and low temperature stress. These results lay a foundation for further functional identification of the stress-related gene WRKY47 in citrus.

    Key words: citrus;WRKY;abiotic stress;gene cloning;expression analysis

    植物的生長發(fā)育過程很大程度上受環(huán)境的影響,其中非生物脅迫主要為干旱、高鹽以及低溫。在非生物脅迫信號的感知以及生理和生化反應上,植物自身已經(jīng)進化出應對脅迫的各種復雜的響應機制。簡單的說,植物可以通過一系列脅迫信號轉導和調控相關基因來響應非生物脅迫[1-2]。這些非生物脅迫響應基因按功能可以分為兩大類,即調控型基因和功能型基因。調控基因編碼的轉錄因子(TFs)可以轉導脅迫信號并調控功能基因的表達,而這些功能基因在減輕植物遭受脅迫帶來的損害中起著直接作用[3]。

    研究發(fā)現(xiàn), bZIP、WRKY、NAC、AP2/ERF和MYB等轉錄因子家族都與植物響應逆境的過程相關 [4]。WRKY是植物中最重要的轉錄因子家族之一。該家族蛋白質包含1到2個由約60個氨基酸組成的高度保守的WRKY DNA結合結構域,還含有C2H2或C2HC鋅指結構,可與DNA序列W-box(TTGACC/T)結合,調控下游基因表達[2]。自1994年首次在甜薯中克隆出1個WRKY 家族轉錄因子基因SPF1之后[5],又陸續(xù)在高等植物以及低等植物(如苔蘚)中克隆出多個WRKY 家族基因[6]。在克隆出來的植物WRKY家族基因中,已有許多被驗證能夠響應非生物脅迫,并在這個過程中發(fā)揮調控作用 [7]。如Qiu等[8]利用Northern blotting技術,對克隆獲得的13個水稻W(wǎng)RKY基因進行分析,發(fā)現(xiàn)在這13個水稻W(wǎng)RKY基因中有10個基因能夠被非生物脅迫(高鹽、高溫、干旱和低溫)誘導表達。Jiang等[9]通過芯片譜分析法對高鹽處理后的擬南芥根系進行分析,發(fā)現(xiàn)有18個WRKY家族基因能夠被高鹽處理誘導表達。同時,單個WRKY家族基因可以被多種非生物脅迫誘導表達。例如櫻桃砧木PcWRKY1基因在滲透脅迫和高鹽脅迫下均能被誘導表達[10];芍藥PLWRKY13可以被低溫、高溫、干旱和鹽脅迫同時誘導表達,并且正向調節(jié)牡丹中番茄紅素的產(chǎn)生[11]。擬南芥AtWRKY47基因可以被硒脅迫、鎘脅迫誘導表達,并且AtWRKY47功能缺失型突變體對硒脅迫表現(xiàn)更敏感,但是對鎘脅迫更耐受,其過表達植株也對鎘脅迫表現(xiàn)出高度耐受[12-13]。前人研究結果表明,水稻OsWRKY47功能缺失型突變體對干旱脅迫敏感并且產(chǎn)量下降,而OsWRKR47過表達植株對干旱脅迫更耐受,說明OsWRKY47基因是干旱脅迫響應的正調控因子[14]。

    柑橘是世界上最重要的水果種類之一,中國柑橘種植面積和產(chǎn)量一直在不斷增長,但其產(chǎn)業(yè)發(fā)展一直受到非生物脅迫的影響和制約。目前,柑橘中的WRKY基因研究還不多,對柑橘WRKY47基因的研究未見報道。因此,本研究從檸檬、甜橙、金柑、蘆柑4個不同柑橘種類中克隆得到WRKY47基因,對克隆出的基因的堿基序列及其編碼蛋白質進行生物信息學分析,并使用實時熒光定量PCR分析這4個基因在非生物脅迫(高鹽、干旱和低溫)下的時空表達模式差異,以期為柑橘抗逆分子育種提供候選基因資源和理論依據(jù)。

    1 材料與方法

    1.1 試驗材料及材料的處理

    將檸檬[Citrus limon (L.) Burm.f.]、甜橙[Citrus sinensis (L.) Osbeck]、蘆柑(Citrus reticulata Blanco)、金柑[Fortunella japonica (Thunb.) Swingle]種子消毒催芽后分別播種于MS固體培養(yǎng)基中,置于28 ℃恒溫培養(yǎng)箱中培養(yǎng)28 d左右,再轉移至Hoagland培養(yǎng)液中預培養(yǎng)2~3 d,然后進行脅迫處理。高鹽、干旱處理時將4種幼苗分別轉移至含有250 mmol/L NaCI、20% PEG-6000的Hoagland液體培養(yǎng)基中常溫培養(yǎng)。低溫處理時幼苗轉移到干凈的Hoagland液體培養(yǎng)基中并置于4 ℃培養(yǎng)箱中培養(yǎng)。每種材料取18~24個植株進行脅迫處理,然后分別于處理0 h(CK),1 h、3 h、6 h、12 h、24 h時取其葉片,-80 ℃儲存,備用。

    1.2 葉片總RNA提取及反轉錄

    檸檬、甜橙、蘆柑、金柑葉片總RNA提取方法采用Trizol法。RNA含量用超微量分光光度計檢測, RNA提取質量用1%瓊脂糖凝膠電泳檢測。用日本TOYOBO公司生產(chǎn)的ReverTra Ace-α-TM試劑盒將提取的總RNA反轉錄成cDNA,存放于-20 ℃,備用。

    1.3 ClWRKY47、 CsWRKY47 、CrWRKY47和 FjWRKY47基因克隆

    由華中農(nóng)業(yè)大學構建的甜橙基因組數(shù)據(jù)庫中含有大量甜橙基因序列,通過序列搜索和拼接獲得包含完整開放閱讀框(ORF)的CsWRKY47轉錄因子基因序列?;谠撔蛄杏肞rimer Premier5.0軟件設計了1對基因特異性引物CsWRKY47-F和CsWRKY47-R(表1)。以獲得的4個物種的葉片cDNA為模板,通過PCR擴增獲得ClWRKY47、CsWRKY47、CrWRKY47、FjWRKY47基因堿基序列全長。PCR擴增程序參考郭文芳[15]的方法,退火溫度56 ℃ 30 s。擴增產(chǎn)物在1 %瓊脂糖凝膠中進行電泳分離,切膠回收后的PCR目的片段用DNA回收試劑盒[天根生化科技(北京)有限公司產(chǎn)品]回收并純化,然后與pMD18-T載體連接[寶日醫(yī)生物技術(北京)有限公司產(chǎn)品]并轉化到大腸桿菌中,菌液PCR鑒定為陽性克隆后送上海生物工程有限公司測序。

    1.4 生物信息學分析

    利用NCBI 中的ORF Finder工具對測序獲得的ClWRKY47、CsWRKY47、CrWRKY47、FjWRKY47基因的cDNA堿基序列進行分析,確定這4個基因的開放閱讀框和所編碼的氨基酸序列,并對其堿基序列進行比對;利用ExPASy ProtParam在線分析工具對蛋白質氨基酸序列及理化性質進行分析;使用DNAMAN軟件進行蛋白質一致性分析和不同植物蛋白質間的多序列比對分析;利用InterProScan在線軟件分析蛋白質結構功能域;采用MEGA5軟件構建4種蛋白質與其他植物蛋白質間的系統(tǒng)進化樹;通過PlantCARE在線網(wǎng)站進行基因的啟動子順式作用元件預測分析;采用TMpred預測基因編碼蛋白質的跨膜結構;通過ProtScale在線軟件預測蛋白質親疏水性;采用在線軟件PSORT II對蛋白質進行亞細胞定位分析;利用SOPMA軟件預測蛋白質的二級結構并用SWISS-MODEL軟件進行蛋白質三級結構預測建模。

    1.5 實時熒光定量PCR

    采用Primer Premier 5.0設計CsWRKY47基因的實時定量引物CsWRKY47-Fq和CsWRKY47-Rq,內參基因為柑橘CsActin基因(表1)。非生物脅迫處理后4個物種WRAKY47基因的表達量在BIO-RAD CFX9 6熒光定量PCR儀上進行檢測。qRT-PCR反應體系的總體積為20.0 μl:10.0 μl熒光染料[寶日醫(yī)生物技術(北京)有限公司產(chǎn)品]、8.0 μl ddH2O、1.0 μl F/R引物、1.0 μl總cDNA。qRT-PCR反應程序參照文獻[15]。每處理樣品設3次重復。相對表達量用2-△△Ct方法計算。

    2 結果與分析

    2.1 CsWRKY47的克隆與分析

    以甜橙為材料,通過RT-PCR擴增克隆獲得CsWRKY47基因堿基序列。測序和生物信息學分析結果顯示,該基因cDNA堿基序列全長1 567 bp,開放閱讀框(ORF)為1 506 bp,編碼501個氨基酸。ExPASy ProtParam在線分析結果顯示,CsWRKY47蛋白分子式為C2 327H3 727N699O757S25,其氨基酸組成有20種,理論等電點(PI)為8.58,PI值大于7說明該蛋白質可能為堿性蛋白質;預測相對分子質量約為54 410,脂肪系數(shù)為64.93,不穩(wěn)定指數(shù)為49.48,總平均親水性為-0.637,一般總平均親水性的值小于0為親水蛋白質,因此推測該蛋白質屬于弱穩(wěn)定性親水蛋白質。通過ProtScale在線軟件預測CsWRKY47蛋白親疏水性,進一步證明其為親水性蛋白質(圖1A)。

    采用TMpred預測蛋白質跨膜結構,結果顯示,CsWRKY47蛋白N端343~363、376~396、461~480這3個位點處具有跨膜結構。PSORT II在線軟件預測結果表明,CsWRKY47蛋白亞細胞定位于細胞核上的可能性最高。利用SOPMA軟件預測CsWRKY47蛋白的二級結構,結果顯示,該蛋白質具有27.94% α-螺旋、8.18%延伸鏈、2.59% β-折疊以及61.28%無規(guī)則卷曲(圖1B)。使用SWISS-MODEL搜索并獲得CsWRKY47蛋白同源模型,推測出該蛋白質可能的4個三級結構模型(圖1C)。通過NetPhos 2.0 Serve預測CsWRKY47蛋白磷酸化位點,結果表明,CsWRKY47蛋白含有26個絲氨酸(Ser)以及2個蘇氨酸(Thr)磷酸化位點。

    2.2 CsWRKY47基因啟動子序列預測分析

    將從甜橙基因組數(shù)據(jù)庫獲取的CsWRKY47基因起始密碼子ATG上游1 500 bp的DNA堿基序列作為該基因的啟動子序列,并通過PlantCARE在線分析網(wǎng)站對該基因啟動子序列進行順式作用元件預測分析,分析結果(表2)顯示,CsWRKY47基因啟動子除了含有16個CAAT-box和30個TATA-box 2種基本的順式作用元件外,還含有多個光響應元件,如G-Box、Sp1、TCCC-motif等;另外在DNA雙鏈上分別有1個WRKY結構域專一結合的W-box元件,同時還包含多個與脅迫相關的順式作用元件,如8個脫落酸響應元件(ABRE)、4個茉莉酸甲酯響應元件(CGTCA-motiF)、1個赤霉素響應元件(TATC-box)、1個抗氧化響應元件(ARE)、4個參與干旱誘導的MYB結合位點(MBS)等。

    2.3 不同種類柑橘WRKY47蛋白的同源性及進化分析

    為了解不同種類柑橘中WRKY47蛋白的差異和功能,通過RT-PCR擴增,分別從檸檬、蘆柑、金柑3個柑橘種類中克隆出WRKY47基因,并分別命名為ClWRKY47、 CrWRKY47和FjWRKY47。測序結果表明,這3個基因的堿基序列與CsWRKY47的cDNA堿基序列長度一致,均為1 567 bp,生物信息學分析發(fā)現(xiàn)開放閱讀框為1 506 bp并且編碼501個氨基酸;對4個基因的堿基序列進行比對發(fā)現(xiàn),這4個基因cDNA堿基序列有16個堿基存在差異,相似性達99.68%(圖2)。

    通過NCBI數(shù)據(jù)庫將ClWRKY47 、CsWRKY47、CrWRKY47、FjWRKY47的堿基序列翻譯成氨基酸序列,并用DNAman對這4條氨基酸序列進行比對,結果顯示,這4條氨基酸序列有7個位點存在差異,序列一致性達99.55%。進一步通過與其他植物WRKY47蛋白進行多序列比對分析(圖3)發(fā)現(xiàn),這4個蛋白質在N端267~333位點都含有1個高度保守的WRKY結構域,并且在該結構域的N末端有1個WRKYGQK保守氨基酸序列和1個C端鋅指結構基序C2H2(CX4-5CX22-23HX1H),說明這4個蛋白質都屬于WRKY轉錄因子家族的Group IIa+IIb類。

    為了驗證ClWRKY47、CsWRKY47、CrWRKY47和FjWRKY47的蛋白質進化關系,使用MEGA 6.0軟件中的最大似然法構建系統(tǒng)進化樹,結果(圖4)顯示,這4個新蛋白質與克萊門柚WRKY47(XP_006444621.1)蛋白親緣關系最近,并且都聚集在同一分支上,與龍眼WRKY47(AEO31517.2)蛋白親緣關系較近,與葡萄WRKY47(XP_002281194.1)蛋白親緣關系最遠。

    2.4 非生物脅迫處理對不同種類柑橘WRKY47基因表達的影響

    分別對甜橙、檸檬、蘆柑和金柑分別進行250 mmol/L NaCl、20% PEG-6000和4 ℃低溫處理,取其葉片進行實時熒光定量表達分析。由圖5a可知,250 mmol/L NaCl脅迫處理植株時,甜橙CsWRKY47在葉片中表達量逐漸升高,24 h時達到最大值,為對照的3倍左右,整體變化較大;檸檬ClWRKY47表達量在1 h時最低,此后逐漸升高,在6 h時達到最高,整體變化趨勢為先升后降;蘆柑CrWRKY47在高鹽處理后,與對照相比,整體表達是下調的,表達量表現(xiàn)出受高鹽處理的抑制;金柑FjWRKY47表達量整體呈現(xiàn)緩慢升高的趨勢,在6 h時達到最大值。

    由圖5b可知,20% PEG6000模擬干旱脅迫處理時,甜橙CsWRKY47表達量在24 h內緩慢升高,在6 h和24 h時出現(xiàn)2次表達高峰;檸檬ClWRKY47的表達量呈現(xiàn)先升高后下降的趨勢,在6 h時表達量達到最高;蘆柑CrWRKY47在24 h內表達量變化趨勢較平緩,在1 h、6 h以及12 h時表達量高于對照;金柑FjWRKY47對干旱脅迫也有響應,但其表達量低于對照,表明干旱脅迫對金柑FjWRKY47的表達有一定的抑制作用。

    由圖5c可知,4 ℃低溫處理時,甜橙CsWRKY47表達量整體呈現(xiàn)先升高后下降的趨勢,在12 h時表達量最高,為對照的2.7倍左右;檸檬ClWRKY47表達量與對照相比總體上變化不大,但在1 h和6 h時其表達量高于對照,說明低溫能誘導其表達;蘆柑CrWRKY47表達量在前12 h處于緩慢升高的趨勢,但在24 h時表達量降低到最低;金柑FjWRKY47的表達量在24 h內都低于對照,說明低溫脅迫抑制FjWRKY47的表達。

    以上結果表明,在高鹽、干旱和低溫處理下,ClWRKY47、CsWRKY47、CrWRKY47和FjWRKY47的表達量都發(fā)生了改變。檸檬ClWRKY47和甜橙CsWRKY47均能被3種非生物脅迫誘導表達;蘆柑CrWRKY47在干旱和低溫脅迫下也能被誘導表達,但在高鹽脅迫下表達量下調;金柑FjWRKY47 可以響應高鹽脅迫,但在干旱和低溫脅迫下,表達都低于對照,尤其是低溫脅迫下,F(xiàn)jWRKY47明顯下調表達。

    3 討論

    WRKY家族轉錄因子具有調控植物生長發(fā)育和響應非生物脅迫等多種功能。但是,目前有關柑橘中WRKY轉錄因子的研究較少。為了初步驗證柑橘中WRKY47基因對逆境的響應機制以及在不同柑橘種類中的差異,我們從檸檬、甜橙、蘆柑、金柑中分別克隆出ClWRKY47、CsWRKY47、CrWRKY47和FjWRKY47基因。根據(jù)所含有的保守結構域的差異,WRKY轉錄因子可分為Group I、Group II和Group III 3類[16]。Group II最初又被分為IIa、IIb、IIc、IId和IIe 5個亞組,但最近的系統(tǒng)發(fā)育分析結果表明,IIa和IIb,IId和IIe分別被合并為IIa +IIb和IId+IIe [17]。在本研究中,多序列比對分析結果表明,這4種蛋白質與其他植物WRKY47蛋白都具有高度保守的WRKY結構域和C2H2鋅指結構,這與前人研究結果[2]相吻合。同時,有研究結果表明,AtWRKY47蛋白屬于Group IIb 類[11],推測ClWRKY47、CsWRKY47、CrWRKY47和FjWRKY47蛋白屬于WRKY轉錄因子家族的Group IIa+IIb類。此外,通過系統(tǒng)進化樹可以明顯看出這4個蛋白質具有種屬特異性,與克萊門柚WRKY47蛋白這個同屬植物蛋白質親緣關系最近且都聚集在進化樹的同一個分支上,這與植物系統(tǒng)分類學的進化趨勢是一致的。另外,通過對甜橙CsWRKY47氨基酸序列的分析發(fā)現(xiàn),該蛋白質定位在細胞核上的可能性最高,符合轉錄因子的一般特征,即亞細胞定位于細胞核上[18-20],該結果還需要通過試驗進一步驗證。

    順式作用元件存在于基因的啟動子中,一般是由5~20個堿基序列組成的特異性DNA序列,具有轉錄調節(jié)功能 [21]。對玉米基因ZmCIPK10和ZmZIP71啟動子序列進行分析,發(fā)現(xiàn)兩者啟動子序列中含有大量的逆境相關元件,如脫落酸、低溫等響應元件,進一步通過非生物脅迫處理植株后,發(fā)現(xiàn)在干旱、低溫和高鹽脅迫下,ZmCIPK10和ZmZIP71基因表達量上升,表明兩者在植物應對逆境的過程中起調控作用[22-23]。在本試驗中,對CsWRKY47基因啟動子順式作用元件進行預測分析,結果表明,該基因啟動子包含多個光響應元件和植物逆境脅迫應答相關作用元件。由此推測,CsWRKY47基因可能受到光、激素及其他非生物脅迫的誘導,參與了植物非生物脅迫響應過程。歐芹在受到病原菌侵染時,PcWRKY1基因在短時間內會迅速表達且表達量顯著增加,對其啟動子序列分析發(fā)現(xiàn),其自身啟動子具有3個W-box元件, PcWRKYl蛋白不僅可以通過自身的啟動子W-box元件與其他WRKY蛋白結合以行使功能,還可以通過該元件與自身結合并形成自我調節(jié)[24] 。本研究中,CsWRKY47基因啟動子也含有W-box元件,猜測CsWRKY47蛋白可能可以與自身啟動子結合以自我調節(jié),或者與其他WRKY蛋白結合以調控其對逆境的響應表達。

    前人研究結果表明,WRKY轉錄因子在植物響應非生物脅迫中扮演著重要的角色,例如在擬南芥中,AtWRKY25和AtWRKY33的表達受到NaCl和干旱處理的誘導并顯著提高[25];蔡榮號等[26]在模式植物擬南芥中通過異源表達玉米ZmWRKY114基因發(fā)現(xiàn)轉基因植株在高鹽脅迫下綠苗率低于野生型植株,并且其根系生長也受到了抑制,ZmWRKY114基因在調控植物應對高鹽脅迫的過程中可能起負調控作用。Zou等[27]對擬南芥的研究結果證實低溫誘導了野生型擬南芥中AtWRKY34基因的表達。本研究結果表明,檸檬ClWRKY47和甜橙CsWRKY47都能被3種非生物脅迫誘導表達,但是蘆柑CrWRKY47與金柑FjWRKY47對不同的脅迫表現(xiàn)出了差異表達。蘆柑CrWRKY47在鹽脅迫后表達下調,在干旱和低溫脅迫下表達上調,金柑FjWRKY47在鹽脅迫后上調表達,在干旱和低溫脅迫后表達下調。植物在應對干旱脅迫和高鹽脅迫的過程中,有很多基因參與響應脅迫,并在這2個脅迫過程中同一個基因起著相反的作用[28-29];也有基因在調控植物應對鹽脅迫、低溫脅迫和機械損傷過程中作用相反,例如楊樹中,PtrWRKY75基因和PtrWRKY80基因的表達能被鹽脅迫誘導,從而上調表達,但這2個基因對低溫脅迫和機械損傷沒有響應;PtrWRKY61基因和PtrWRKY88基因則特異性響應鹽脅迫和機械傷害,對低溫脅迫不響應[30]。因此,蘆柑CrWRKY47與金柑FjWRKY47在不同的非生物脅迫處理后出現(xiàn)差異表達的原因可能是蘆柑和金柑這2個柑橘種類自身遺傳特性導致的,也可能是由于植物在不同的脅迫條件下,與不同柑橘種類WRKY47基因啟動子結合的上游調控因子不同導致的。

    綜上所述,本研究克隆獲得的ClWRKY47、CsWRKY47、CrWRKY47和FjWRKY47基因編碼的轉錄因子屬于WRKY轉錄因子家族的Group IIa+IIb類。順式作用元件預測發(fā)現(xiàn)CsWRKY47基因啟動子含有多種逆境相關元件。qRT-RCR表達分析結果揭示,CrWRKY47與FjWRKY47對不同非生物脅迫的響應不同。這些結果表明柑橘WRKY47可以響應多種非生物脅迫,可能在柑橘抵御非生物脅迫過程中發(fā)揮重要的作用。

    參考文獻:

    [1] HIRAYAMA T, SHINOZAKI K. Research on plant abiotic stress responses in the post-genome era: past, present and future[J]. The Plant Journal, 2010,61(6): 1041-1052.

    [2] LIU J H, PENG T, DAI W. Critical cis-acting elements and interacting transcription factors: key players associated with abiotic stress responses in plants[J]. Plant Molecular Biology Reporter, 2014, 32(2):303-317.

    [3] REN X, CHEN Z, LIU Y, et al. ABO3, a WRKY transcription factor, mediates plant responses to abscisic acid and drought tolerance in arabidopsis[J]. Plant Journal for Cell & Molecular Biology, 2010, 63(3):417-429.

    [4] XU Z S, CHEN M, LI L C,? et al. Functions and application of the AP2/ERF transcription factor family in crop improvcmcn[J].J Integr Plaint Piol,2011,53(7):570-585.

    [5] ISHIGURO S, NAKAMURA K. Characterization of a cDNA encoding a novel DNA-binding protein, SPF1, that recognizes SP8 sequences in the 5′ upstream regions of genes coding for sporamin and beta-amylase from sweet potato [J]. Molecular & General Genetics Mgg, 1994, 244(6):563-571.

    [6] EULGEM T , RUSHTON P J , ROBATZEK S , et al. The WRKY superfamily of plant transcription factors.[J]. Trends in Plant ence, 2000, 5(5):199-206.

    [7] CHEN L G, SONG Y, LI S J, et al. The role of WRKY transcription factors in plant abiotic stresses [J]. Biochim Biophys Acta, 2012, 1819 (2): 120-128.

    [8] QIU Y , JING S , FU J , et al.? Cloning and analysis of expression profile of 13 WRKY genes in rice[J]. Chinese Science Bulletin, 2004, 49(20):2159-2168.

    [9] JIANG Y, DEYHOLOS M K. Comprehensive transcriptional profiling of NaCl-stressed arabidopsis roots reveals novel classes of responsive genes[J]. BMC Plant Biol, 2006,6:25.

    [10]徐 麗,陳 新,宗曉娟,等. 櫻桃砧木PcWRKY1基因的克隆與表達分析[J].江蘇農(nóng)業(yè)學報,2018,34(3):636-641.

    [11]WANG X, LI J, GUO X, et al. PlWRKY13: a transcription factor involved in abiotic and biotic stress responses in paeonia lactiflora[J]. Int J Mol Sci, 2019,20(23):5953.

    [12]錢靚雯. WRKY47在植物硒富積及耐受中的作用研究[D].合肥:合肥工業(yè)大學,2018.

    [13]孟 云. 擬南芥WRKY47基因在鎘脅迫應答中的作用機理研究[D].合肥:合肥工業(yè)大學,2019.

    [14]RAINERI J, WANG S, PELEG Z, et al. The rice transcription factor OsWRKY47 is a positive regulator of the response to water deficit stress[J]. Plant Molecular Biology, 2015, 88(4/5):401-413.

    [15]郭文芳. 柑橘抗逆相關基因NAC83、MYB15和COR15a的克隆與表達分析[D].南昌:江西農(nóng)業(yè)大學,2015.

    [16]LLORCA C M, POTSCHIN M, ZENTGRAF U. bZIPs and WRKYs: two large transcription factor families executing two different functional strategies[J]. Frontiers in Plant Science, 2014,5:169.

    [17]RUSHTON D L, TRIPATHI P, RABARA R C, et al. WRKY transcription factors: key components inabscisic acid signaling[J]. Plant Biotechnol, 2012,10:2-11.

    [18]LIU Y, YANG T, LIN Z, et al. A WRKY transcription factor PbrWRKY53 from Pyrus betulaefolia is involved in drought tolerance and AsA accumulation[J]. Plant Biotechnol J, 2019,17:1770-1787.

    [19]GUO H, WANG Y, WANG L, et al. Expression of the MYB transcription factor gene BplMYB46 affects abiotic stress tolerance and secondary cell wall deposition in Betula platyphylla[J]. Plant Biotechnol J, 2017,15:107-121.

    [20]魏 鑫,王寒濤,魏恒玲,等. 陸地棉GhWRKY33的克隆及抗旱功能分析[J].中國農(nóng)業(yè)科學,2020,53(22):4537-4549.

    [21]郭晉艷,鄭曉瑜,鄒翠霞,等. 植物非生物脅迫誘導啟動子順式元件及轉錄因子研究進展[[J].生物技術通報,2011 (4):16-20.

    [22]趙晉鋒,余愛麗,王寒玉,等. 非生物逆境脅迫下ZmCIPK10基因表達分析[J].生物技術進展,2011,1(2):130-134.

    [23]劉彥丹,英 生,張登峰,等. 玉米逆境脅迫響應基因ZmbZIP71的克隆與表達分析[J].植物遺傳資源學報,2011,12(5):775-781.

    [24]TURCK F, ZHOU A, SOMSSICH I E. Stimulus-dependent, promoter-specific binding of transcription factor WRKY1 to Its native promoter and the defense-related gene PcPRl-1 in Parsley[J]. Plant Cell, 2004,16:2573-2585.

    [25]LI S, FU Q, CHEN L, et al. Arabidopsis thaliana WRKY25, WRKY26, and WRKY33 coordinate induction of plant thermotolerance[J]. Planta,2011, 233(6):1237-1252.

    [26]蔡榮號,李 尉,陳浩偉,等. 異源表達玉米ZmWRKY114基因增強擬南芥對鹽脅迫的敏感性[J].安徽農(nóng)業(yè)大學學報,2019,46(6):1040-1047.

    [27]ZOU C, JIANG W, YU D. Male gametophyte-specific WRKY34 transcription factor mediates cold sensitivity of mature pollen in Arabidopsis[J]. J Exp Bot, 2010,61(14):3901-3914.

    [28]XIANG Y, HUANG Y, XIONG L. Characterization of stress-responsive CIPK genes in rice for stress tolerance improvement[J]. Plant physiology, 2007, 144(3):1416.

    [29]LEE S K, KIM B G, KWON T R, et al. Overexpression of the mitogen-activated protein kinase gene OsMAPK33 enhances sensitivity to salt stress in rice (Oryza sativa L.)[J]. J Biosci, 2011,36(1):139-151.

    [30]JIANG Y Z, DUAN Y J, YIN J, et al. Genome-wide identification and characterization of the populus WRKY transcription factor family and analysis of their expression in response to biotic and abiotic stresses[J]. Journal of Experimental Botany,2014,9(22):22.

    (責任編輯:陳海霞)

    收稿日期:2020-10-13

    基金項目:國家重點研發(fā)計劃項目(2019YFD1000100);國家自然科學基金項目(31701896);江西省自然科學基金面上項目(20202BAB205001);江西省柑橘產(chǎn)業(yè)技術體系項目(JXARS-07-栽培崗位)

    作者簡介:沈 丹(1996-),女,江蘇溧陽人,碩士研究生,研究方向為果樹抗逆生理。(E-mail)sd17751769221@163.com

    通訊作者:劉德春,(E-mail)ldc873380800@163.com;劉 勇,(E-mail)liuyongjxau@163.com

    猜你喜歡
    非生物脅迫表達分析基因克隆
    白木香轉錄因子AsMYB1和AsMYB2克隆及表達分析
    三個小麥防御素基因的克隆及序列分析
    紅花生育酚環(huán)化酶基因的克隆及表達分析
    非生物脅迫對擬南芥IQM4基因表達的影響
    科技視界(2016年16期)2016-06-29 11:55:38
    水稻富亮氨酸類受體蛋白激酶參與外界脅迫應答的研究進展
    玉米紋枯病病菌y—谷氨酰轉肽酶基因克隆與表達分析
    信號分子與葉銹菌誘導下小麥病程相關蛋白1基因的表達分析
    水稻OsCATA基因相互作用蛋白的篩選和表達分析
    山桃醇腈酶基因PdHNL1的克隆、序列分析與原核表達
    紅花低分子量油體蛋白基因的克隆與系統(tǒng)進化分析
    亚洲欧美精品自产自拍| 久久青草综合色| 插逼视频在线观看| 人妻少妇偷人精品九色| 亚洲内射少妇av| 国产精品无大码| 欧美人与善性xxx| 如何舔出高潮| 在线观看三级黄色| 欧美日韩精品成人综合77777| 久久久久久久久久人人人人人人| 亚洲av不卡在线观看| 2021少妇久久久久久久久久久| 亚洲精品色激情综合| 两个人的视频大全免费| 少妇的逼水好多| 男人狂女人下面高潮的视频| 美女cb高潮喷水在线观看| 日韩三级伦理在线观看| 久久国产乱子免费精品| 高清黄色对白视频在线免费看 | 亚洲真实伦在线观看| 少妇高潮的动态图| h视频一区二区三区| 国产精品一二三区在线看| 色视频在线一区二区三区| av一本久久久久| 十分钟在线观看高清视频www | 少妇裸体淫交视频免费看高清| 国内少妇人妻偷人精品xxx网站| 蜜桃在线观看..| 91精品一卡2卡3卡4卡| av在线老鸭窝| 少妇熟女欧美另类| 男男h啪啪无遮挡| 777米奇影视久久| 久久久久久久久久人人人人人人| 欧美老熟妇乱子伦牲交| 一级爰片在线观看| 欧美日韩亚洲高清精品| 亚洲精品日本国产第一区| 美女内射精品一级片tv| 大片电影免费在线观看免费| 亚洲欧洲国产日韩| 精品国产一区二区三区久久久樱花| 两个人免费观看高清视频 | 欧美性感艳星| 久久国产乱子免费精品| 欧美精品一区二区大全| 免费大片黄手机在线观看| 久久国产精品男人的天堂亚洲 | 国产永久视频网站| 亚洲,欧美,日韩| 三上悠亚av全集在线观看 | 中文资源天堂在线| 国产亚洲一区二区精品| 热re99久久国产66热| 国产一区二区在线观看av| 2022亚洲国产成人精品| 国产有黄有色有爽视频| 人妻制服诱惑在线中文字幕| 久久这里有精品视频免费| 特大巨黑吊av在线直播| 最近2019中文字幕mv第一页| 久久国内精品自在自线图片| 色视频在线一区二区三区| 一级毛片久久久久久久久女| 亚洲精品国产色婷婷电影| 自拍欧美九色日韩亚洲蝌蚪91 | 晚上一个人看的免费电影| 三级国产精品片| 97精品久久久久久久久久精品| 精品一区二区三区视频在线| 亚洲色图综合在线观看| 一级二级三级毛片免费看| 国产成人精品婷婷| 国产国拍精品亚洲av在线观看| 一区二区av电影网| 日韩一本色道免费dvd| 天天躁夜夜躁狠狠久久av| 久久久久久久久久人人人人人人| 免费看光身美女| 人妻系列 视频| 欧美成人精品欧美一级黄| 国产爽快片一区二区三区| 99热这里只有是精品50| 另类亚洲欧美激情| a级片在线免费高清观看视频| 精品酒店卫生间| 久久人妻熟女aⅴ| 亚洲精品日韩av片在线观看| 看非洲黑人一级黄片| 极品人妻少妇av视频| 久久久久久伊人网av| 人人妻人人澡人人看| 18禁动态无遮挡网站| 日本vs欧美在线观看视频 | av免费在线看不卡| 成年美女黄网站色视频大全免费 | 日韩亚洲欧美综合| av又黄又爽大尺度在线免费看| 天堂中文最新版在线下载| 国产精品99久久久久久久久| 超碰97精品在线观看| 免费播放大片免费观看视频在线观看| 伊人亚洲综合成人网| 99re6热这里在线精品视频| 欧美精品高潮呻吟av久久| 精品人妻熟女av久视频| 免费观看性生交大片5| 国产国拍精品亚洲av在线观看| 午夜福利在线观看免费完整高清在| h日本视频在线播放| 少妇被粗大的猛进出69影院 | 国产视频首页在线观看| 久久精品久久久久久噜噜老黄| 三级国产精品欧美在线观看| 黄色一级大片看看| 97在线人人人人妻| 亚洲av不卡在线观看| 乱系列少妇在线播放| 好男人视频免费观看在线| 国产高清有码在线观看视频| 人人澡人人妻人| 9色porny在线观看| 各种免费的搞黄视频| 日韩中文字幕视频在线看片| 一本—道久久a久久精品蜜桃钙片| 国产熟女午夜一区二区三区 | 亚洲一区二区三区欧美精品| 秋霞伦理黄片| 亚洲精品一区蜜桃| 如何舔出高潮| 97精品久久久久久久久久精品| 国产精品嫩草影院av在线观看| 啦啦啦视频在线资源免费观看| 中国国产av一级| 99久久精品热视频| 人体艺术视频欧美日本| 亚洲久久久国产精品| 性色av一级| 久久人人爽人人片av| 老司机亚洲免费影院| 最后的刺客免费高清国语| 毛片一级片免费看久久久久| 久久ye,这里只有精品| 中国国产av一级| 人人澡人人妻人| 高清不卡的av网站| 大片免费播放器 马上看| 日日爽夜夜爽网站| 最近中文字幕高清免费大全6| 亚洲,一卡二卡三卡| 国产精品女同一区二区软件| av一本久久久久| 国产高清国产精品国产三级| 大码成人一级视频| 日韩大片免费观看网站| 99久久中文字幕三级久久日本| 久久久久视频综合| 国产成人91sexporn| 国产精品秋霞免费鲁丝片| 日韩大片免费观看网站| www.色视频.com| 国产精品秋霞免费鲁丝片| 少妇丰满av| 五月伊人婷婷丁香| 一区二区三区精品91| 国产极品天堂在线| a级毛片在线看网站| 国产精品一区二区性色av| av有码第一页| 免费大片黄手机在线观看| 国产精品国产av在线观看| 日本欧美视频一区| 免费黄频网站在线观看国产| 成年av动漫网址| 国产精品蜜桃在线观看| 黄色视频在线播放观看不卡| 免费观看性生交大片5| 成人影院久久| 成人影院久久| 观看美女的网站| 天堂8中文在线网| 成人特级av手机在线观看| 亚洲精品国产色婷婷电影| 成人漫画全彩无遮挡| 少妇的逼好多水| 一级毛片 在线播放| 国产精品国产三级国产av玫瑰| 午夜久久久在线观看| 国产亚洲最大av| 爱豆传媒免费全集在线观看| 国产黄频视频在线观看| 婷婷色麻豆天堂久久| 亚洲第一区二区三区不卡| 成人免费观看视频高清| 日本免费在线观看一区| 两个人的视频大全免费| 亚洲在久久综合| 啦啦啦啦在线视频资源| 欧美成人精品欧美一级黄| 水蜜桃什么品种好| 欧美另类一区| 日韩中文字幕视频在线看片| 亚洲精品国产成人久久av| 日本猛色少妇xxxxx猛交久久| 亚洲va在线va天堂va国产| 成年人午夜在线观看视频| 特大巨黑吊av在线直播| 亚洲丝袜综合中文字幕| 全区人妻精品视频| 日韩伦理黄色片| 一边亲一边摸免费视频| 天美传媒精品一区二区| av不卡在线播放| 麻豆成人午夜福利视频| 男人狂女人下面高潮的视频| 亚洲国产欧美日韩在线播放 | 80岁老熟妇乱子伦牲交| 免费看不卡的av| 波野结衣二区三区在线| 人妻夜夜爽99麻豆av| 人妻夜夜爽99麻豆av| 免费观看a级毛片全部| 七月丁香在线播放| 日本av手机在线免费观看| av天堂中文字幕网| 最新的欧美精品一区二区| 又粗又硬又长又爽又黄的视频| 国产精品女同一区二区软件| 激情五月婷婷亚洲| 好男人视频免费观看在线| 国产 精品1| 中文字幕免费在线视频6| 日韩视频在线欧美| 熟女电影av网| 91aial.com中文字幕在线观看| 日本wwww免费看| 午夜激情久久久久久久| av国产精品久久久久影院| 性色avwww在线观看| 水蜜桃什么品种好| 另类精品久久| 精品人妻一区二区三区麻豆| 国产伦理片在线播放av一区| 亚洲第一区二区三区不卡| 国产午夜精品久久久久久一区二区三区| 91久久精品国产一区二区三区| 一区二区三区免费毛片| 韩国av在线不卡| 中文在线观看免费www的网站| 日日摸夜夜添夜夜添av毛片| 美女脱内裤让男人舔精品视频| 美女主播在线视频| 日韩av免费高清视频| 99久久精品热视频| 国产一区亚洲一区在线观看| 91久久精品国产一区二区成人| 久久人人爽人人片av| 久久久久久久亚洲中文字幕| av视频免费观看在线观看| 国产伦在线观看视频一区| 久久久国产精品麻豆| 亚洲av国产av综合av卡| 国内揄拍国产精品人妻在线| 国国产精品蜜臀av免费| 亚洲欧美日韩卡通动漫| 久久久久国产精品人妻一区二区| av黄色大香蕉| 人妻 亚洲 视频| 亚洲欧美一区二区三区黑人 | 一级爰片在线观看| 五月开心婷婷网| 日韩一本色道免费dvd| 亚洲av成人精品一二三区| 国产精品人妻久久久影院| 精品少妇内射三级| 啦啦啦啦在线视频资源| 九九爱精品视频在线观看| 久久这里有精品视频免费| 成年人免费黄色播放视频 | 亚洲人成网站在线观看播放| 日韩一区二区视频免费看| 亚洲国产日韩一区二区| 日本免费在线观看一区| 欧美日韩视频精品一区| 久久综合国产亚洲精品| 亚洲综合精品二区| 国产精品久久久久久久久免| 午夜视频国产福利| 久久99一区二区三区| 国产极品粉嫩免费观看在线 | 大又大粗又爽又黄少妇毛片口| 人妻制服诱惑在线中文字幕| 亚洲美女黄色视频免费看| 日韩中文字幕视频在线看片| 毛片一级片免费看久久久久| 如日韩欧美国产精品一区二区三区 | 国产日韩欧美亚洲二区| 国产乱人偷精品视频| 国产成人精品福利久久| 人妻夜夜爽99麻豆av| 一区在线观看完整版| 亚洲精品,欧美精品| 国产成人精品婷婷| 午夜福利视频精品| 嘟嘟电影网在线观看| 极品少妇高潮喷水抽搐| av在线老鸭窝| 人人妻人人看人人澡| 免费少妇av软件| 看十八女毛片水多多多| 在线观看免费视频网站a站| 涩涩av久久男人的天堂| 三级国产精品片| 精品99又大又爽又粗少妇毛片| 一级av片app| 中文在线观看免费www的网站| 色视频www国产| 亚洲三级黄色毛片| 99视频精品全部免费 在线| 国产一区二区在线观看日韩| 精品一区在线观看国产| 最近2019中文字幕mv第一页| 两个人免费观看高清视频 | 日韩中字成人| 国产精品三级大全| 99九九在线精品视频 | 日韩熟女老妇一区二区性免费视频| 亚洲,一卡二卡三卡| 男女边摸边吃奶| 精品国产乱码久久久久久小说| 六月丁香七月| 久久6这里有精品| av专区在线播放| 这个男人来自地球电影免费观看 | 久久久国产欧美日韩av| 寂寞人妻少妇视频99o| 偷拍熟女少妇极品色| 建设人人有责人人尽责人人享有的| 在线免费观看不下载黄p国产| 在线免费观看不下载黄p国产| 18禁动态无遮挡网站| 在线观看av片永久免费下载| 大香蕉久久网| 欧美xxⅹ黑人| 成人黄色视频免费在线看| 一级毛片电影观看| 日日爽夜夜爽网站| av国产精品久久久久影院| 在线 av 中文字幕| a 毛片基地| 80岁老熟妇乱子伦牲交| 亚洲美女搞黄在线观看| 亚洲,欧美,日韩| av在线观看视频网站免费| 青青草视频在线视频观看| 亚洲av免费高清在线观看| 女的被弄到高潮叫床怎么办| 日韩成人伦理影院| 亚洲av欧美aⅴ国产| 国产真实伦视频高清在线观看| 2018国产大陆天天弄谢| 午夜福利视频精品| 久久精品久久久久久久性| 国产色婷婷99| 免费看光身美女| 精品人妻偷拍中文字幕| 亚洲美女黄色视频免费看| 18禁裸乳无遮挡动漫免费视频| 青春草亚洲视频在线观看| 91精品国产国语对白视频| 精品少妇久久久久久888优播| 一级二级三级毛片免费看| 男男h啪啪无遮挡| 午夜福利影视在线免费观看| 建设人人有责人人尽责人人享有的| av线在线观看网站| 国产伦理片在线播放av一区| 国产男人的电影天堂91| 熟妇人妻不卡中文字幕| 街头女战士在线观看网站| 热99国产精品久久久久久7| 在线播放无遮挡| 国产视频首页在线观看| 中文字幕久久专区| 草草在线视频免费看| 99热这里只有精品一区| 69精品国产乱码久久久| 久久精品夜色国产| 天天操日日干夜夜撸| 精品少妇内射三级| 99久久精品一区二区三区| 自拍偷自拍亚洲精品老妇| 黄色视频在线播放观看不卡| 人体艺术视频欧美日本| av在线播放精品| 精品少妇黑人巨大在线播放| 久久久久久久久大av| 特大巨黑吊av在线直播| 纵有疾风起免费观看全集完整版| 综合色丁香网| 日韩在线高清观看一区二区三区| 国内揄拍国产精品人妻在线| 欧美日韩亚洲高清精品| 婷婷色综合大香蕉| 九九爱精品视频在线观看| 午夜av观看不卡| 欧美 亚洲 国产 日韩一| 久久99一区二区三区| 一边亲一边摸免费视频| 成年女人在线观看亚洲视频| 爱豆传媒免费全集在线观看| 老熟女久久久| 精品一区在线观看国产| 国产亚洲一区二区精品| 永久免费av网站大全| 一区二区三区免费毛片| 日韩在线高清观看一区二区三区| 男女边吃奶边做爰视频| 美女视频免费永久观看网站| 久久久久精品久久久久真实原创| 国产无遮挡羞羞视频在线观看| 日韩熟女老妇一区二区性免费视频| 国产精品偷伦视频观看了| 国产熟女午夜一区二区三区 | 久久ye,这里只有精品| 中文精品一卡2卡3卡4更新| xxx大片免费视频| 综合色丁香网| 我要看黄色一级片免费的| 精品久久久精品久久久| 欧美精品亚洲一区二区| 久久精品国产鲁丝片午夜精品| 中文字幕久久专区| 精品少妇黑人巨大在线播放| 国产在线视频一区二区| 交换朋友夫妻互换小说| 国产极品粉嫩免费观看在线 | 欧美另类一区| 亚洲av综合色区一区| 国产精品一区二区三区四区免费观看| 26uuu在线亚洲综合色| 久久久久人妻精品一区果冻| 国产中年淑女户外野战色| 在线观看一区二区三区激情| av有码第一页| 男人舔奶头视频| 九色成人免费人妻av| 久久久精品94久久精品| 嫩草影院新地址| 成人国产av品久久久| 欧美少妇被猛烈插入视频| 老司机亚洲免费影院| 国产毛片在线视频| 国产成人精品久久久久久| 激情五月婷婷亚洲| 精品99又大又爽又粗少妇毛片| 日韩熟女老妇一区二区性免费视频| 美女视频免费永久观看网站| 亚洲激情五月婷婷啪啪| 亚洲精品成人av观看孕妇| 国产精品国产三级专区第一集| 狂野欧美激情性bbbbbb| a级毛片在线看网站| 水蜜桃什么品种好| 男的添女的下面高潮视频| 下体分泌物呈黄色| 永久网站在线| 精品久久久精品久久久| 亚洲经典国产精华液单| 欧美日韩国产mv在线观看视频| 久久女婷五月综合色啪小说| 国模一区二区三区四区视频| 久久毛片免费看一区二区三区| 日韩成人av中文字幕在线观看| 内射极品少妇av片p| 天堂8中文在线网| 亚洲人成网站在线观看播放| 国产精品国产三级国产av玫瑰| 国产淫语在线视频| 激情五月婷婷亚洲| 国产成人一区二区在线| 国产精品欧美亚洲77777| 国产69精品久久久久777片| 免费大片黄手机在线观看| 国产一区二区在线观看日韩| 一本大道久久a久久精品| 免费播放大片免费观看视频在线观看| 亚洲伊人久久精品综合| 亚洲av福利一区| 午夜激情福利司机影院| 久久婷婷青草| 久久久亚洲精品成人影院| 日韩强制内射视频| 韩国高清视频一区二区三区| 赤兔流量卡办理| 成人毛片a级毛片在线播放| 在线观看美女被高潮喷水网站| 国产视频内射| 高清视频免费观看一区二区| 91午夜精品亚洲一区二区三区| av又黄又爽大尺度在线免费看| 亚洲av免费高清在线观看| 久久婷婷青草| 国产精品成人在线| 久久久久久久精品精品| 狠狠精品人妻久久久久久综合| 久久婷婷青草| av在线观看视频网站免费| 2021少妇久久久久久久久久久| 久久亚洲国产成人精品v| 日韩制服骚丝袜av| 黄色视频在线播放观看不卡| av有码第一页| 中国国产av一级| 欧美精品高潮呻吟av久久| 欧美bdsm另类| 人妻系列 视频| 国产探花极品一区二区| 一边亲一边摸免费视频| 中文字幕制服av| 成人毛片60女人毛片免费| 亚洲色图综合在线观看| 日韩视频在线欧美| 少妇人妻久久综合中文| 国产极品粉嫩免费观看在线 | 久久人人爽人人片av| 国产伦在线观看视频一区| 国产男人的电影天堂91| 一本大道久久a久久精品| 热99国产精品久久久久久7| 美女中出高潮动态图| 国产极品粉嫩免费观看在线 | 亚洲精品自拍成人| 一级毛片久久久久久久久女| 国产深夜福利视频在线观看| 久久女婷五月综合色啪小说| 久久久精品94久久精品| 少妇精品久久久久久久| 欧美日韩精品成人综合77777| 51国产日韩欧美| 久久精品国产亚洲av涩爱| 国国产精品蜜臀av免费| 国产亚洲av片在线观看秒播厂| 国产精品无大码| 国模一区二区三区四区视频| 黄色配什么色好看| 青春草亚洲视频在线观看| 亚洲av免费高清在线观看| 欧美日本中文国产一区发布| 91在线精品国自产拍蜜月| 亚洲成人av在线免费| 国产精品久久久久成人av| 天天躁夜夜躁狠狠久久av| 国产精品国产三级国产专区5o| 午夜日本视频在线| 久久精品国产鲁丝片午夜精品| 国产亚洲91精品色在线| 色5月婷婷丁香| 久久久久久久精品精品| 插逼视频在线观看| 老司机影院成人| 一区在线观看完整版| 日本91视频免费播放| 777米奇影视久久| 精品亚洲成国产av| 亚洲自偷自拍三级| 亚洲性久久影院| 久久午夜福利片| 中文字幕人妻丝袜制服| 麻豆乱淫一区二区| 亚洲成人av在线免费| 六月丁香七月| 一本久久精品| 国产精品麻豆人妻色哟哟久久| av免费观看日本| 久久精品国产亚洲网站| 精品人妻偷拍中文字幕| av国产久精品久网站免费入址| a级片在线免费高清观看视频| 国产欧美另类精品又又久久亚洲欧美| 最近2019中文字幕mv第一页| 国产在线视频一区二区| 国产伦精品一区二区三区视频9| 免费观看的影片在线观看| 日韩成人伦理影院| 中文精品一卡2卡3卡4更新| 中文字幕制服av| 成人黄色视频免费在线看| 国产中年淑女户外野战色| 我的老师免费观看完整版| 国产免费视频播放在线视频| 看免费成人av毛片| 天堂8中文在线网| 亚洲人成网站在线观看播放| 欧美bdsm另类| 久久久久久久国产电影| 日本午夜av视频| av黄色大香蕉| 一边亲一边摸免费视频| kizo精华| 亚洲成人手机| .国产精品久久| 99热全是精品| 亚洲精品久久久久久婷婷小说| 亚洲美女搞黄在线观看| 在线观看三级黄色| 国产成人一区二区在线| 一区二区三区免费毛片| a级毛片在线看网站| 亚洲人成网站在线播| 伦理电影免费视频| 午夜福利在线观看免费完整高清在|