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

    干旱脅迫下大苞萱草WRKY 基因家族成員鑒定及生物信息學分析

    2023-04-07 07:52:40陳麗飛李嘉峻劉云怡慧周蘊薇
    河南農(nóng)業(yè)科學 2023年2期
    關鍵詞:萱草基序結構域

    陳麗飛,李嘉峻,劉云怡慧,白 云,孟 緣,李 江,周蘊薇

    (1.吉林農(nóng)業(yè)大學 林學與草學學院,吉林 長春 130118;2.定西市水土保持科學研究所,甘肅 定西 743000;3.吉林農(nóng)業(yè)大學 園藝學院,吉林 長春 130118)

    萱 草(Hemerocallisspp.)屬 阿 福 花 科(Asphodelaceae)萱草屬(Hemerocallis),多年生草本植物,民間多稱其“忘憂草”[1-2]。萱草屬植物分布廣泛、適應性強、品種繁多、觀賞價值高,被稱為中國的“母親花”。大苞萱草(Hemerocallis middendorffii)是現(xiàn)代萱草的原始親本之一,主要分布在我國東北、華北等地。大苞萱草抗旱能力強、花色艷麗,而且具有二次開花特性,可廣泛應用于各類園林綠地或進行盆栽裝飾。

    干旱不僅限制著植物的正常生長發(fā)育,還會影響植物體的代謝功能,降低光合速率、蒸騰速率以及抑制植物生長[3]。干旱是植物非生物脅迫的主要原因之一,迫使植物發(fā)生形態(tài)和生理變化[4]。植物對干旱脅迫的響應主要表現(xiàn)為根系和葉片結構的改變,生理水平主要表現(xiàn)為光合作用、滲透調(diào)節(jié)代謝、抗氧化代謝及激素物質等方面的變化,而分子水平則是調(diào)控抗旱調(diào)節(jié)基因和功能基因的表達,改變功能性蛋白和有機分子合成酶的編碼,進而影響干旱信號轉導以及激素、滲透物質的合成和積累[5]。

    WRKY 轉錄因子的主要特征在于其包含60 個左右氨基酸DNA 結合區(qū)的WRKY 域,其中包含一個高度保守的七肽序列(WRKYGQK)。WRKY 轉錄因子在植物防御生物及非生物脅迫中發(fā)揮著關鍵作用[6]。超表達的WRKY 家族基因能夠幫助植物在干旱環(huán)境中維持較好的水分狀況,緩解因干旱所造成的萎蔫及生長抑制[7]。在大豆研究中,過表達的GmWRKY54可能通過激活ABA 以及Ca2+信號通路來控制氣孔的關閉,增強大豆的抗旱能力[8]。在干旱脅迫下,植物會受到諸如活性氧(ROS)損傷等次生脅迫損傷。VlWRKY48能提高葡萄過氧化氫酶、過氧化物酶和超氧化物歧化酶等抗氧化酶活性,清除活性氧來提高其抗旱性[9]。

    目前,對于大苞萱草的開發(fā)利用還不充分,關于大苞萱草在分子層面響應干旱脅迫的研究鮮見報道,不利于深入研究大苞萱草的抗旱機制。鑒于此,以大苞萱草為試驗材料,基于轉錄組數(shù)據(jù),采用生物信息學分析方法對干旱脅迫下的大苞萱草WRKY 基因(HmWRKY)進行鑒定、構建進化樹、分析和預測蛋白質理化性質以及蛋白質二級結構等,根據(jù)HmWRKY基因的表達特性,明確候選基因,為進一步研究HmWRKY基因功能奠定基礎。

    1 材料和方法

    1.1 材料及處理

    供試材料大苞萱草取自吉林農(nóng)業(yè)大學野生花卉資源圃,選取田間長勢一致且植株健壯的大苞萱草幼苗移栽至高13 cm、直徑10 cm 的塑料花盆中,單株單盆,栽培基質為草炭、珍珠巖按體積比3∶1混合配制成的營養(yǎng)土,進行正常的盆栽管理。緩苗后,將供試材料分為對照組(CK 組)和干旱脅迫處理組(D組)。干旱處理前使各組均處于飽和含水狀態(tài),處理后,CK 組每3 d 澆一次水,D 組采取自然失水脅迫的方式處理。于5 d 后的10:00 對CK 組和D組進行取樣,每組處理設置3個重復,將嫩葉和新根放入液氮中帶回并存放在-80 ℃冰箱保存,備用。RNA 的提取、質控、建庫及轉錄組的三代測序均由北京諾禾致源科技股份有限公司完成。

    1.2 HmWRKY基因家族成員鑒定

    以大苞萱草轉錄組測序結果為材料。在TAIR網(wǎng) 站(https://www.arabidopsis.org/)下 載 擬 南 芥WRKY 基因家族的轉錄因子序列和蛋白質序列,共72 個。通過本地blast+進行比對,篩選出大苞萱草的WRKY 轉錄因子。利用Smart(http://smart.emblheidelberg.de/)在線軟件進行鑒定以確保含有WRKY 保守域,通過NCBI(https://blast.ncbi.nlm.nih.gov)在線軟件進行結構域完整性的預測。

    1.3 HmWRKY基因家族成員生物信息學分析

    1.3.1 HmWRKY 基因家族聚類分析及多重序列比對 利用MEGA 7 對HmWRKY 轉錄因子蛋白構建進 化 樹(Neighbor-joining,NJ),Bootstrap 設 置 為1 000次重復,其余參數(shù)選擇默認值[10]。利用序列分析軟件DNAMAN 6.0 對HmWRKY 蛋白序列進行多重序列比對,截取包含WRKY 結構域與鋅指結構的氨基酸序列,通過TBtools 軟件的SeqLogo 繪制HmWRKY轉錄因子保守結構域的logo圖。

    1.3.2 HmWRKY 基因家族的理化性質分析 利用ExPASy ProtParam tool(https://web.expasy.org/protparam/)對HmWRKY蛋白的理化性質進行預測,其中包括氨基酸序列長度、分子質量、蛋白質等電點、總平均親水性、不穩(wěn)定系數(shù)以及脂溶性系數(shù)。使 用CELLO v.2.5(http://cello.life.nctu.edu.tw/)對HmWRKY 家族成員進行植物細胞的亞細胞定位,整理數(shù)據(jù)進行匯總。

    1.3.3 HmWRKY 蛋白二級結構預測 利用SOPMA(https://npsa-prabi.ibcp.fr/cgi-bin/npsa_automat.pl?page=npsa-sopma.html)在線軟件對大苞萱草WRKY蛋白質二級結構進行預測。

    1.3.4 HmWRKY 蛋白保守基序分析 通過在線軟件MEME(https://meme-suite.org/meme/tools/meme)對HmWRKY 蛋白保守基序進行搜索,搜索參數(shù)設置為最大搜索數(shù)量20 個,ZOOPS 參數(shù)設置為Zero or one occurrence per sequence,其余參數(shù)使用網(wǎng)頁默認設置。最后使用TBtools將HmWRKY轉錄因子蛋白保守基序可視化。

    1.3.5 HmWRKY 基因家族在葉和根中的表達分析 利用北京諾禾致源科技股份有限公司提供的售后平臺工具NovoMagic(https://magic.novogene.com/customer/main#/login)分別制作處理組與對照組葉和根中HmWRKY基因家族41個成員的表達模式熱聚圖。

    1.4 HmWRKY 基因實時熒光定量PCR(qRTPCR)分析

    將公司測序后剩余的樣品進行液氮研磨,使用RNApure Plant Kit(DNase Ⅰ)植物RNA 提取試劑盒(DNase Ⅰ)提取葉片總RNA,使用NovoScript?Plus All-in-one 1st Strand cDNA Synthesis(gDNA Purge)反轉錄試劑盒合成cDNA,使用NovoStart?SYBR qPCR SuperMix Plus 試劑盒進行qRT-PCR,使用Primer 5.0進行qRT-PCR 引物設計(表1),內(nèi)參基因為HfEF-1α(GenBank No.MT096368)[11]。qRT-PCR反應體系為2×SYBR Mix 主混合物10 μL,上下游引物各1 μL,模板(cDNA)2 μL,ddH2O 6 μL。反應程序:預變性95 ℃1 min;95 ℃20 s、60 ℃20 s、72 ℃30 s,40 個循環(huán)。試驗進行3 次生物學重復,采用2-ΔΔCt方法計算基因的相對表達量[12]。

    表1 HmWRKY基因qRT-PCR引物序列Tab.1 qRT-PCR primer sequences for HmWRKY gene

    2 結果與分析

    2.1 HmWRKY基因家族成員鑒定

    通過本地blast+與擬南芥WRKY 轉錄因子進行比對,從轉錄組測序結果中共篩選出182 個大苞萱草WRKY 轉錄因子。利用在線軟件Smart鑒定以確保含有WRKY 保守域,通過NCBI 在線軟件預測結構域的完整性。最終得到41 個大苞萱草WRKY 轉錄因子,并將這些基因命名為HmWRKY1—HmWRKY41。

    2.2 HmWRKY基因家族成員生物信息學分析

    2.2.1 HmWRKY 基因家族聚類分析及多重序列比對 將41 個HmWRKY 蛋白與72 個擬南芥WRKY蛋白進行多序列分析并通過軟件構建進化樹(圖1)。WRKY DNA 結構域(DBD)的特征是在其N 端有一個不變的七肽WRKYGQK 氨基酸基序,在其C端包含一個鋅指結構。根據(jù)WRKY 結構域的數(shù)量和鋅指狀基序的特征將WRKY 轉錄因子分為Ⅰ、Ⅱ、Ⅲ三類,其中Ⅱ類分為5個亞類。

    圖1 大苞萱草和擬南芥WRKY轉錄因子聚類分析Fig.1 Cluster analysis of WRKY transcription factors in Hemerocallis middendorffii and Arabidopsis thaliana

    WRKY 結構域高度保守,多重序列比對結果如圖2 和圖3 所示。WRKYGQK 的氨基酸序列存在替換的現(xiàn)象。Ⅰ類共有8個HmWRKY轉錄因子,含有2 個WRKY 結構域,鋅指結構為CX7CX24~25HXH 型。Ⅱ類共有30 個HmWRKY 轉錄因子,只含有1 個WRKY 結構域且鋅指結構為CX7CX25HXH 型,Ⅰ類和Ⅱ類轉錄因子均有不同程度的七肽結構變異和鋅指結構缺失的現(xiàn)象。Ⅲ類有4 個HmWRKY 轉錄因子,含有1 個WRKY 結構域,鋅指結構為CX7CX25HXC 型,Ⅲ類HmWRKY 轉錄因子相對保守,七肽結構不發(fā)生變異且鋅指結構完整。

    圖2 HmWRKY轉錄因子家族保守結構域分析Fig.2 Conservative domain analysis of HmWRKY transcription factor family

    圖3 HmWRKY轉錄因子家族保守結構域分析Fig.3 Fig.3 Conservative domain analysis of HmWRKY transcription factor family

    2.2.2 HmWRKY 基因家族的理化性質分析HmWRKY 基因家族成員的理化性質差異性較大(表2),HmWRKY 轉錄因子的氨基酸數(shù)目在92(HmWRKY11)~641(HmWRKY7),分 子 質 量 在10 551.76(HmWRKY11)~69 972.1(HmWRKY7)ku,理論等電點在4.97(HmWRKY26)~9.99(HmWRKY24),脂溶性指數(shù)在47.08(HmWRKY26)~75.93(HmWRKY6),41 個HmWRKY 轉錄因子的親水性平均值均為負值,判定為親水性蛋白。除HmWRKY2、HmWRKY4、HmWRKY5、HmWRKY27、HmWRKY41 為穩(wěn)定蛋白外,其余均為不穩(wěn)定蛋白(不穩(wěn)定系數(shù)>40)。HmWRKY 家族成員亞細胞定位結果如表3所示,除HmWRKY41 定位在線粒體中外,其余的HmWRKY轉錄因子均定位在細胞核中。

    表2 HmWRKY轉錄因子蛋白的理化性質Tab.2 Physicochemical properties of HmWRKY transcription factor proteins

    續(xù)表2 HmWRKY轉錄因子蛋白的理化性質Tab.2(Continued)Physicochemical properties of HmWRKY transcription factor proteins

    2.2.3 HmWRKY 蛋白的二級結構預測 41 個HmWRKY 轉錄因子蛋白二級結構均含有α-螺旋、β-折疊、延伸鏈以及無規(guī)卷曲,無其他二級結構(表4)。絕大多數(shù)HmWRKY 蛋白都表現(xiàn)出無規(guī)卷曲數(shù)>α-螺旋數(shù)≥延伸鏈數(shù)>β-折疊數(shù)的特點,但少數(shù)蛋白質,如HmWRKY11、HmWRKY16、HmWRKY17等表現(xiàn)出無規(guī)卷曲數(shù)>延伸鏈數(shù)>α-螺旋數(shù)>β-折疊數(shù)的特點。

    表4 HmWRKY轉錄因子蛋白二級結構Tab.4 Secondary structure of HmWRKY transcription factor proteins %

    續(xù)表4 HmWRKY轉錄因子蛋白二級結構Tab.4(Continued)Secondary structure of HmWRKY transcription factor proteins %

    2.2.4 HmWRKY 蛋白保守基序分析 HmWRKY蛋白序列的保守基序預測結果見圖4。Motif 1 是WRKY 轉錄因子的特征基序且組成了WRKY 結構域。除HmWRKY24、HmWRKY26、HmWRKY27 和HmWRKY30外,其他成員均包含Motif 1。Motif 2和Motif 5 共同構成C2H2型的鋅指結構域,存在于大多數(shù)HmWRKY 轉錄因子中。Ⅰ類成員幾乎全包含Motif 5、Motif 17 以 及 在N 端、C 端 分 別 含 有1 組Motif 1 和Motif 2。Ⅱ類成員的Motif種類繁多,但呈現(xiàn)出在進化樹中親緣關系較近的各亞類成員,包含Motif 種類也基本一致的規(guī)律,如Ⅱa 亞類與Ⅱb 亞類 除HmWRKY2 外,均 包 含Motif 19;Ⅱc 亞 類HmWRKY33 和HmWRKY39 的Motif 種類完全相同且成員多數(shù)包含Motif 4;Ⅱd 亞類HmWRKY14 和HmWRKY37包含的Motif種類以及排列順序完全相同,其他成員幾乎都包含Motif 7。Ⅲ類成員中存在的Motif種類比較單一,只包含Motif 1和Motif 2。

    圖4 HmWRKY轉錄因子家族保守基序Fig.4 Conserved motif of HmWRKY transcription factor family

    2.2.5 HmWRKY 基因家族在葉和根中的表達分析 根據(jù)轉錄組數(shù)據(jù)的分析,對比HmWRKY 基因家族的41 個成員在不同處理下葉和根中的表達見圖5,在葉中有10 個上調(diào)基因,29 個下調(diào)基因,其中HmWRKY25和HmWRKY32在 葉 中 無 表 達。HmWRKY22、HmWRKY9上 調(diào) 較 明 顯,HmWRKY2、HmWRKY11、HmWRKY12、HmWRKY20、HmWRKY26、HmWRKY30、HmWRKY34和HmWRKY36下 調(diào) 較 明顯。41 個HmWRKY基因在根中均有表達,其中有20 個上調(diào)基因,21 個下調(diào)基因。HmWRKY4、HmWRKY8、HmWRKY14、HmWRKY15、HmWRKY22、HmWRKY23、HmWRKY24、HmWRKY35在受到干旱脅迫時,表達量上調(diào)較為明顯。HmWRKY2、HmWRKY16、HmWRKY27、HmWRKY30、HmWRKY41在下調(diào)基因中表達量變化較為明顯。

    圖5 大苞萱草葉片和根中HmWRKY基因在不同處理下的表達模式熱聚圖Fig.5 Thermopolygram of HmWRKY gene expression pattern under different treatments in leaves and roots of Hemerocallis middendorffii

    HmWRKY差異基因的表達變化如圖6 所示,在葉中共發(fā)現(xiàn)10 個差異基因,1 個上調(diào)基因,9 個下調(diào)基因。大苞萱草受到干旱脅迫后,只有HmWRKY9的表達量升高。下調(diào)的差異基因中HmWRKY2、HmWRKY26和HmWRKY30的表達量變化顯著。在根中發(fā)現(xiàn)4 個差異基因,其中3 個上調(diào)基因,1 個下調(diào)基因。HmWRKY4和HmWRKY24上調(diào)趨勢比較明顯,可能在根系中起到積極的調(diào)控作用。此外,通過葉和根的差異基因比對發(fā)現(xiàn),HmWRKY2在2個組織內(nèi)均有表達且均為負表達,表達量變化十分明顯。

    圖6 大苞萱草葉片和根中HmWRKY差異基因在不同處理下的表達水平Fig.6 Expression of HmWRKY differential gene under different treatments in leaves and roots of Hemerocallis middendorffii

    2.3 候選差異基因的qRT-PCR驗證

    為了驗證候選差異基因HmWRKY9在葉中表達量的正確性,進行qRT-PCR 分析。結果如圖7 所示,HmWRKY9在受到干旱脅迫后,表達量升高,正向響應了干旱脅迫。HmWRKY9的qRT-PCR 結果與轉錄組測序的變化趨勢基本一致,從而驗證了轉錄組數(shù)據(jù)的準確性。

    3 結論與討論

    WRKY 基因家族在植物體中響應生長發(fā)育、生物脅迫以及非生物脅迫等過程,在非生物脅迫中對干旱脅迫反應的調(diào)控網(wǎng)絡發(fā)揮重要作用。CmWRKY10通過ABA 途徑來調(diào)控菊花的干旱耐受性[13],水稻在干旱脅迫條件下過表達OsWRKY11可增強其耐旱性[14],過表達AtWRKY30轉基因植株能夠顯著提高小麥的耐旱能力[15]。以上研究均揭示了WRKY轉錄因子家族成員可提高植物抗旱能力。

    本研究基于轉錄組測序,共挖掘出41 個HmWRKY 轉錄因子。將HmWRKY 基因家族與已公布的擬南芥WRKY 基因共同構建進化樹并進行聚類分析,HmWRKY 基因家族的分類與擬南芥WRKY 基因家族的分類一致,表明本試驗轉錄組數(shù)據(jù)具有可靠性,HmWRKY 基因家族在進化上相對保守。WRKY 轉錄因子的七肽和鋅指結構是WRKY 基因行使功能的兩大基團,隨著植物進化,這兩大基團的氨基酸組成也有所變化。Ⅰ類HmWRKY 轉錄因子中只有HmWRKY41N 端的七肽結構發(fā)生變異,變異為“WRKYGPK”,Ⅱ類HmWRKY 轉錄因子的七肽結構易發(fā)生變異,其中變異為“WKKYGQK、WRKYGKK、WRKYGEK”[16-17]。Ⅲ類HmWRKY 轉錄因子的結構域相比前兩類保守,沒有發(fā)生變異情況。ZHANG 等[16]針對真核生物WRKY 基因的研究發(fā)現(xiàn),WRKY 家族成員有19 個WRKY 結構域的變異,WRKYGEK 和WRKYGKK 是2 個常見的變異。EULGEM 等[17]在香蕉和楊樹中也發(fā)現(xiàn)七肽變異類型為WRKYGKK 的WRKY 轉錄因子蛋白,其中的變異類型與本研究結果一致。YU等[18]認為高度保守的WRKY 結構域可以特異識別靶基因啟動子區(qū)內(nèi)的W-BOX 元件從而調(diào)控該基因的活性,進一步發(fā)揮其生物學功能。因此,具有變異的WRKY 結構域的HmWRKY 轉錄因子可能會影響其發(fā)揮生物學功能。保守基序分析表明,絕大多數(shù)HmWRKY 蛋白含有的特征基序構成不同,如Ⅰ類成員多數(shù)含有2 組Motif 1、Motif 2 和Motif 5,Ⅱa和Ⅱb 亞類的特征基序是Motif 19,Ⅱc 亞類特征基序為Motif 4,Ⅱd和Ⅱe亞類的特征基序均為Motif 7,而Ⅲ類成員中只存在2 種Motif,保守基序種類不同說明這些成員在進化過程中發(fā)生了不同程度的變異,這可能與基因家族的分類及功能調(diào)控有關。

    通過對干旱脅迫下HmWRKY基因的表達分析發(fā)現(xiàn),41個HmWRKY基因在葉和根中的表達量變化均有不同,少數(shù)HmWRKY基因如HmWRKY6、HmWRKY29等在干旱脅迫前后葉和根中的表達量變化不明顯,而HmWRKY22、HmWRKY36等在處理前后葉和根中的表達量變化相對顯著。另外,HmWRKY6在干旱處理前不表達,處理后開始表達,HmWRKY21和HmWRKY38則相反,說明HmWRKY6、HmWRKY21和HmWRKY38可能與大苞萱草的抗旱性相關。已有研究表明,擬南芥中部分WRKY 轉錄因子可在植物抵御干旱脅迫的過程中起到重要調(diào)控作用。AtWRKY28和AtWRKY57均能提高植株的耐旱性[19-20]。過表達AtWRKY29可使擬南芥植株變小且葉片顏色變深,增強其抗旱能力[21]。AtWRKY33可通過直接調(diào)控CesA8 參與干旱脅迫[22]。AtWRKY40可以通過增強植株的滲透調(diào)節(jié)能力和抗氧化能力來維持細胞吸水,從而提高植株對干旱脅迫的抵抗力[23]。過表達AtWRKY57的擬南芥耐旱性顯著提高[19]。本研究通過大苞萱草轉錄組測序結果得到13個差異基因,部分差異基因與上述擬南芥抗旱基因具有較近的親緣關系,如HmWRKY2、HmWRKY3、HmWRKY4、HmWRKY9、HmWRKY30和HmWRKY34。這 也 間 接 表 明 ,HmWRKY2、HmWRKY3、HmWRKY4、HmWRKY9、HmWRKY30和HmWRKY34可能通過促進或抑制自身的表達來參與大苞萱草響應干旱脅迫的過程。大苞萱草受到干旱脅迫的前后,HmWRKY2均存在于葉和根中,HmWRKY2的表達量變化均為顯著下降且HmWRKY2與AtWRKY40同 源 性 較 高。因 此,HmWRKY2可能通過抑制其表達在大苞萱草抵御干旱脅迫的過程中發(fā)揮重要作用。

    本研究基于大苞萱草的轉錄組數(shù)據(jù),首次對大苞萱草WRKY 基因家族進行鑒定,初步確定了41個HmWRKY基因。通過對HmWRKY 轉錄因子的結構特征及在不同處理下葉和根中表達量的綜合分析,揭示了大苞萱草WRKY 基因家族成員可能與大苞萱草抗旱方面相關,為HmWRKY 轉錄因子家族的功能研究以及抗性基因篩選提供參考。確定HmWRKY9為候選基因,后續(xù)將開展HmWRKY9的功能驗證工作,探討其是否與大苞萱草的抗旱性有關。

    猜你喜歡
    萱草基序結構域
    EPIYA 基序與幽門螺桿菌感染相關胃病關系的研究進展
    新醫(yī)學(2023年10期)2023-12-09 15:04:51
    帶TRS基序突變的新型冠狀病毒威脅更大
    芥藍Aux/IAA家族基因生物信息學與表達分析
    四月節(jié)令百草之萱草 誰言四月芳菲盡 萱草葵花覺昨非
    紫禁城(2020年5期)2021-01-07 02:13:34
    蛋白質結構域劃分方法及在線服務綜述
    重組綠豆BBI(6-33)結構域的抗腫瘤作用分析
    北堂種萱草
    海峽姐妹(2016年7期)2016-02-27 15:21:17
    組蛋白甲基化酶Set2片段調(diào)控SET結構域催化活性的探討
    秋葉有情,萱草有愛
    人間(2015年19期)2016-01-04 12:46:44
    不同品種大花萱草抗寒性比較
    亚洲精华国产精华精| 老司机福利观看| 欧美国产日韩亚洲一区| 欧美绝顶高潮抽搐喷水| 色在线成人网| 久久中文看片网| 一个人观看的视频www高清免费观看 | 在线播放国产精品三级| 午夜福利影视在线免费观看| 国产亚洲精品第一综合不卡| 国产精品av久久久久免费| 视频在线观看一区二区三区| 久久伊人香网站| 亚洲最大成人中文| 一区福利在线观看| 俄罗斯特黄特色一大片| 亚洲欧美激情综合另类| 国产熟女xx| 久久欧美精品欧美久久欧美| 黄色视频,在线免费观看| 女性生殖器流出的白浆| 欧美成狂野欧美在线观看| 久久这里只有精品19| 久久这里只有精品19| 国产成人免费无遮挡视频| 搞女人的毛片| 搡老岳熟女国产| 国产精品精品国产色婷婷| 人人妻人人澡人人看| 亚洲 欧美一区二区三区| 动漫黄色视频在线观看| 深夜精品福利| 久99久视频精品免费| 一级毛片女人18水好多| 女人被躁到高潮嗷嗷叫费观| 精品乱码久久久久久99久播| 欧美人与性动交α欧美精品济南到| 欧美久久黑人一区二区| 12—13女人毛片做爰片一| 男女床上黄色一级片免费看| 日韩有码中文字幕| 黄色a级毛片大全视频| 日本黄色视频三级网站网址| 欧美日本中文国产一区发布| 黄色片一级片一级黄色片| 窝窝影院91人妻| 12—13女人毛片做爰片一| 亚洲av成人av| 欧美乱妇无乱码| 国产在线精品亚洲第一网站| 悠悠久久av| www.999成人在线观看| 色精品久久人妻99蜜桃| www.熟女人妻精品国产| 身体一侧抽搐| 变态另类丝袜制服| 久99久视频精品免费| 国产成人免费无遮挡视频| 999久久久国产精品视频| 天堂√8在线中文| 午夜精品在线福利| 怎么达到女性高潮| 女人被躁到高潮嗷嗷叫费观| 亚洲九九香蕉| 色综合欧美亚洲国产小说| 欧美大码av| 国产精品亚洲美女久久久| 国产野战对白在线观看| 久久久久久久精品吃奶| 国产aⅴ精品一区二区三区波| 国产免费男女视频| 亚洲视频免费观看视频| 99香蕉大伊视频| 在线播放国产精品三级| 一进一出抽搐gif免费好疼| 亚洲男人天堂网一区| 两性夫妻黄色片| 免费久久久久久久精品成人欧美视频| 久久精品91蜜桃| 中文字幕人妻丝袜一区二区| aaaaa片日本免费| 久久精品国产综合久久久| 国产男靠女视频免费网站| 国产精品自产拍在线观看55亚洲| 亚洲男人的天堂狠狠| av免费在线观看网站| 日韩欧美一区二区三区在线观看| 一区在线观看完整版| 国产精品电影一区二区三区| 久久久国产欧美日韩av| 日本免费a在线| 老司机午夜福利在线观看视频| www日本在线高清视频| 国产又爽黄色视频| 又大又爽又粗| 啦啦啦韩国在线观看视频| 好男人电影高清在线观看| 国产精品爽爽va在线观看网站 | 国产亚洲欧美精品永久| 美女国产高潮福利片在线看| 中文字幕人妻丝袜一区二区| 国产亚洲欧美98| 亚洲va日本ⅴa欧美va伊人久久| 中文字幕人成人乱码亚洲影| 久久久久国内视频| 国产成人免费无遮挡视频| 人人妻,人人澡人人爽秒播| 99国产综合亚洲精品| 性欧美人与动物交配| 欧美激情久久久久久爽电影 | 亚洲 欧美一区二区三区| 亚洲国产欧美日韩在线播放| 欧美黄色片欧美黄色片| 老鸭窝网址在线观看| 午夜福利视频1000在线观看 | 国产精品乱码一区二三区的特点 | 久久久久久久久久久久大奶| 18禁黄网站禁片午夜丰满| 色综合欧美亚洲国产小说| 看片在线看免费视频| 一级毛片女人18水好多| 天天一区二区日本电影三级 | 亚洲欧美激情综合另类| 中文字幕av电影在线播放| 午夜激情av网站| 国产精品影院久久| 成人免费观看视频高清| 视频在线观看一区二区三区| 在线观看免费午夜福利视频| 伦理电影免费视频| 精品日产1卡2卡| 三级毛片av免费| 亚洲人成77777在线视频| 少妇裸体淫交视频免费看高清 | 亚洲狠狠婷婷综合久久图片| 99精品久久久久人妻精品| 最近最新中文字幕大全电影3 | www.精华液| 又黄又粗又硬又大视频| 美女国产高潮福利片在线看| 亚洲伊人色综图| 亚洲一区二区三区不卡视频| 无限看片的www在线观看| 亚洲av成人不卡在线观看播放网| 长腿黑丝高跟| 大香蕉久久成人网| 国产麻豆成人av免费视频| 久久热在线av| 久久久久国产精品人妻aⅴ院| 国产精品亚洲一级av第二区| 久久人人97超碰香蕉20202| 久热这里只有精品99| 欧美日本中文国产一区发布| 日韩大码丰满熟妇| 女人被狂操c到高潮| 可以在线观看毛片的网站| 侵犯人妻中文字幕一二三四区| 精品人妻1区二区| 亚洲无线在线观看| 日韩av在线大香蕉| 午夜影院日韩av| 99在线人妻在线中文字幕| 成人手机av| 免费久久久久久久精品成人欧美视频| 久久精品国产综合久久久| 久久国产精品男人的天堂亚洲| 电影成人av| 国产日韩一区二区三区精品不卡| 日本黄色视频三级网站网址| 亚洲成av片中文字幕在线观看| 最近最新免费中文字幕在线| 久久久国产成人精品二区| 国产亚洲精品一区二区www| 久久香蕉精品热| 久久久久亚洲av毛片大全| 精品一品国产午夜福利视频| 91成年电影在线观看| 久久人妻熟女aⅴ| 亚洲色图 男人天堂 中文字幕| 高清在线国产一区| 欧美精品啪啪一区二区三区| 久久国产乱子伦精品免费另类| 亚洲片人在线观看| 大陆偷拍与自拍| 琪琪午夜伦伦电影理论片6080| 久久精品人人爽人人爽视色| 亚洲熟妇中文字幕五十中出| 在线天堂中文资源库| 午夜免费激情av| 亚洲国产欧美一区二区综合| 一二三四社区在线视频社区8| 国产精品爽爽va在线观看网站 | 亚洲熟女毛片儿| 亚洲全国av大片| 天堂动漫精品| 99国产精品一区二区蜜桃av| 欧美色欧美亚洲另类二区 | 99国产综合亚洲精品| 69精品国产乱码久久久| 国产精品一区二区在线不卡| 757午夜福利合集在线观看| 人人妻,人人澡人人爽秒播| 久久久久久国产a免费观看| 大香蕉久久成人网| 一卡2卡三卡四卡精品乱码亚洲| 午夜福利影视在线免费观看| 久久久久久人人人人人| 日本精品一区二区三区蜜桃| 一级片免费观看大全| 亚洲国产中文字幕在线视频| 韩国av一区二区三区四区| 日韩欧美一区视频在线观看| 日本撒尿小便嘘嘘汇集6| 女人爽到高潮嗷嗷叫在线视频| 日韩中文字幕欧美一区二区| 欧美另类亚洲清纯唯美| 人人澡人人妻人| 亚洲欧美激情在线| 他把我摸到了高潮在线观看| av超薄肉色丝袜交足视频| 老司机午夜福利在线观看视频| 欧美不卡视频在线免费观看 | www.精华液| 一区二区三区激情视频| 成人亚洲精品一区在线观看| 成人国语在线视频| 又黄又爽又免费观看的视频| 乱人伦中国视频| 日本a在线网址| 亚洲成a人片在线一区二区| 久久国产亚洲av麻豆专区| 国产乱人伦免费视频| av视频免费观看在线观看| 免费观看人在逋| 性少妇av在线| 色综合婷婷激情| 91字幕亚洲| 久久久国产精品麻豆| 日本vs欧美在线观看视频| 丝袜美足系列| 久久久久亚洲av毛片大全| 欧美激情极品国产一区二区三区| 午夜福利视频1000在线观看 | 91麻豆精品激情在线观看国产| 啪啪无遮挡十八禁网站| 亚洲avbb在线观看| 熟妇人妻久久中文字幕3abv| 巨乳人妻的诱惑在线观看| 日本免费a在线| 亚洲第一欧美日韩一区二区三区| 在线国产一区二区在线| 国产精品影院久久| 欧美激情高清一区二区三区| 色av中文字幕| 国产精品免费视频内射| 可以免费在线观看a视频的电影网站| 国产aⅴ精品一区二区三区波| 国产精品国产高清国产av| 村上凉子中文字幕在线| 99国产精品免费福利视频| 日韩 欧美 亚洲 中文字幕| 精品国产超薄肉色丝袜足j| 亚洲美女黄片视频| 亚洲无线在线观看| 后天国语完整版免费观看| 50天的宝宝边吃奶边哭怎么回事| 午夜两性在线视频| 人人妻人人爽人人添夜夜欢视频| 十八禁人妻一区二区| 精品第一国产精品| 国产精品亚洲一级av第二区| 欧美老熟妇乱子伦牲交| a在线观看视频网站| 亚洲国产精品久久男人天堂| av免费在线观看网站| 日韩欧美一区二区三区在线观看| 亚洲人成电影观看| 欧美最黄视频在线播放免费| 欧美+亚洲+日韩+国产| 欧美日韩乱码在线| 天堂影院成人在线观看| 成人特级黄色片久久久久久久| 国产午夜精品久久久久久| 久久国产精品男人的天堂亚洲| 一区福利在线观看| 男女下面插进去视频免费观看| 日本在线视频免费播放| 999久久久精品免费观看国产| 嫩草影院精品99| 欧美成人免费av一区二区三区| 欧美精品亚洲一区二区| av有码第一页| 成熟少妇高潮喷水视频| 亚洲 欧美 日韩 在线 免费| 色综合站精品国产| 大码成人一级视频| 精品国产亚洲在线| videosex国产| 别揉我奶头~嗯~啊~动态视频| 搞女人的毛片| 国产私拍福利视频在线观看| av片东京热男人的天堂| 国语自产精品视频在线第100页| 99在线人妻在线中文字幕| 操美女的视频在线观看| 亚洲五月色婷婷综合| 久久天躁狠狠躁夜夜2o2o| 日本撒尿小便嘘嘘汇集6| 亚洲午夜精品一区,二区,三区| www.999成人在线观看| 中文字幕色久视频| 国产亚洲av高清不卡| 欧美日韩中文字幕国产精品一区二区三区 | 久9热在线精品视频| 中文字幕精品免费在线观看视频| 久久中文看片网| 一进一出抽搐动态| 美女国产高潮福利片在线看| 91精品三级在线观看| 可以免费在线观看a视频的电影网站| 中文亚洲av片在线观看爽| 黄色视频,在线免费观看| 制服人妻中文乱码| 亚洲 国产 在线| 大陆偷拍与自拍| 日本三级黄在线观看| 国产欧美日韩一区二区三区在线| 亚洲av第一区精品v没综合| 熟女少妇亚洲综合色aaa.| 免费无遮挡裸体视频| 欧美黑人欧美精品刺激| 国产欧美日韩一区二区精品| 亚洲av成人一区二区三| 国产在线观看jvid| 亚洲国产欧美一区二区综合| 国产精品一区二区免费欧美| 男人舔女人的私密视频| 一级,二级,三级黄色视频| 中文字幕高清在线视频| 国产精品电影一区二区三区| 国产亚洲精品一区二区www| 亚洲片人在线观看| 两性夫妻黄色片| 亚洲中文日韩欧美视频| 精品国产一区二区三区四区第35| www.www免费av| 午夜日韩欧美国产| 精品人妻在线不人妻| 久久午夜亚洲精品久久| 久久精品亚洲精品国产色婷小说| av视频免费观看在线观看| 嫁个100分男人电影在线观看| 真人做人爱边吃奶动态| 久久国产乱子伦精品免费另类| av天堂在线播放| 亚洲,欧美精品.| 国产国语露脸激情在线看| 香蕉久久夜色| 亚洲精品久久成人aⅴ小说| 国产欧美日韩一区二区精品| 免费在线观看视频国产中文字幕亚洲| 国产蜜桃级精品一区二区三区| 国产成人一区二区三区免费视频网站| 9热在线视频观看99| 变态另类丝袜制服| 大码成人一级视频| 999久久久精品免费观看国产| 午夜亚洲福利在线播放| 亚洲情色 制服丝袜| 亚洲国产精品999在线| 一区二区三区精品91| 欧美乱色亚洲激情| 变态另类丝袜制服| 亚洲精品国产色婷婷电影| 涩涩av久久男人的天堂| 青草久久国产| 久久精品亚洲精品国产色婷小说| 久久久水蜜桃国产精品网| www.熟女人妻精品国产| 91成人精品电影| 麻豆av在线久日| 九色亚洲精品在线播放| 午夜福利视频1000在线观看 | 1024视频免费在线观看| 久久香蕉精品热| 亚洲五月婷婷丁香| 老司机靠b影院| av福利片在线| 国产野战对白在线观看| 一卡2卡三卡四卡精品乱码亚洲| 午夜a级毛片| 丰满人妻熟妇乱又伦精品不卡| 一进一出好大好爽视频| 日韩大尺度精品在线看网址 | 老司机深夜福利视频在线观看| 亚洲av第一区精品v没综合| 欧美激情久久久久久爽电影 | 国产欧美日韩一区二区三区在线| 极品人妻少妇av视频| 日本五十路高清| 国产97色在线日韩免费| 一区福利在线观看| 9色porny在线观看| 久久久久久久午夜电影| 女人爽到高潮嗷嗷叫在线视频| 黄频高清免费视频| 亚洲专区国产一区二区| 精品欧美国产一区二区三| 精品电影一区二区在线| 首页视频小说图片口味搜索| www日本在线高清视频| 国产精品久久久人人做人人爽| 18禁国产床啪视频网站| 亚洲少妇的诱惑av| av网站免费在线观看视频| 怎么达到女性高潮| 纯流量卡能插随身wifi吗| 成人手机av| 动漫黄色视频在线观看| 久久久久久久久久久久大奶| 精品第一国产精品| 久久久国产欧美日韩av| 欧美人与性动交α欧美精品济南到| 电影成人av| 国产精品美女特级片免费视频播放器 | 免费在线观看完整版高清| 国产精品香港三级国产av潘金莲| 给我免费播放毛片高清在线观看| 91在线观看av| www.999成人在线观看| 日本三级黄在线观看| 精品久久久久久,| 一进一出抽搐动态| 身体一侧抽搐| 午夜福利在线观看吧| 亚洲最大成人中文| 一个人免费在线观看的高清视频| xxx96com| av视频免费观看在线观看| 成年女人毛片免费观看观看9| 黄色视频不卡| 高潮久久久久久久久久久不卡| 日韩三级视频一区二区三区| 91九色精品人成在线观看| 国产伦人伦偷精品视频| 国产精品1区2区在线观看.| 男人的好看免费观看在线视频 | 色综合婷婷激情| 性色av乱码一区二区三区2| 国产成人一区二区三区免费视频网站| 亚洲欧美日韩无卡精品| 亚洲国产精品久久男人天堂| 欧美黑人欧美精品刺激| 9191精品国产免费久久| 久久久久久久久久久久大奶| 亚洲中文字幕一区二区三区有码在线看 | 亚洲av五月六月丁香网| 成年人黄色毛片网站| 成人国产一区最新在线观看| 精品欧美一区二区三区在线| 国产一卡二卡三卡精品| 久久久久久久久中文| 日日干狠狠操夜夜爽| 少妇 在线观看| 好男人在线观看高清免费视频 | 中文字幕另类日韩欧美亚洲嫩草| 国语自产精品视频在线第100页| 两性夫妻黄色片| 久久欧美精品欧美久久欧美| 成熟少妇高潮喷水视频| 中文字幕人妻丝袜一区二区| 韩国av一区二区三区四区| 色播亚洲综合网| 亚洲三区欧美一区| 国产麻豆69| 成人三级黄色视频| 国产日韩一区二区三区精品不卡| 亚洲av日韩精品久久久久久密| 亚洲va日本ⅴa欧美va伊人久久| 在线观看午夜福利视频| 日韩免费av在线播放| 久久草成人影院| 国产欧美日韩一区二区三| 两个人看的免费小视频| 国产精品,欧美在线| 亚洲欧洲精品一区二区精品久久久| 久久热在线av| 村上凉子中文字幕在线| 亚洲欧美激情综合另类| 欧美一区二区精品小视频在线| 波多野结衣一区麻豆| 亚洲精品在线美女| 日韩高清综合在线| 女同久久另类99精品国产91| 日本 欧美在线| 国产乱人伦免费视频| 老鸭窝网址在线观看| 国产乱人伦免费视频| 久久亚洲真实| 国产精品一区二区精品视频观看| 久久久久久久精品吃奶| 色老头精品视频在线观看| 亚洲一码二码三码区别大吗| 欧美国产精品va在线观看不卡| 少妇裸体淫交视频免费看高清 | 91麻豆精品激情在线观看国产| 正在播放国产对白刺激| 丰满的人妻完整版| 久久人妻av系列| 在线免费观看的www视频| 国产精品一区二区免费欧美| 欧美黑人欧美精品刺激| 国内精品久久久久精免费| 美女高潮喷水抽搐中文字幕| 看片在线看免费视频| 人成视频在线观看免费观看| 亚洲电影在线观看av| 久9热在线精品视频| ponron亚洲| 久久精品aⅴ一区二区三区四区| 欧美人与性动交α欧美精品济南到| 不卡一级毛片| 青草久久国产| 国产精品美女特级片免费视频播放器 | 嫩草影院精品99| 两个人看的免费小视频| 中文字幕人成人乱码亚洲影| 亚洲免费av在线视频| 波多野结衣av一区二区av| 女人精品久久久久毛片| 久久人人精品亚洲av| av片东京热男人的天堂| 男女午夜视频在线观看| e午夜精品久久久久久久| 村上凉子中文字幕在线| 精品久久蜜臀av无| 午夜福利欧美成人| 别揉我奶头~嗯~啊~动态视频| 免费在线观看亚洲国产| 亚洲少妇的诱惑av| 成人国语在线视频| 午夜a级毛片| 国产aⅴ精品一区二区三区波| 十八禁网站免费在线| 麻豆一二三区av精品| 国产av精品麻豆| 757午夜福利合集在线观看| 老鸭窝网址在线观看| 精品一区二区三区视频在线观看免费| 国产熟女午夜一区二区三区| 亚洲精品在线观看二区| 99久久综合精品五月天人人| 国产高清激情床上av| 多毛熟女@视频| 两个人免费观看高清视频| 97超级碰碰碰精品色视频在线观看| 日韩免费av在线播放| 99久久国产精品久久久| 人人妻人人澡欧美一区二区 | 亚洲欧美日韩另类电影网站| a级毛片在线看网站| 97人妻天天添夜夜摸| 又大又爽又粗| 丝袜在线中文字幕| 精品无人区乱码1区二区| 亚洲美女黄片视频| 校园春色视频在线观看| 欧美午夜高清在线| 国产精品电影一区二区三区| 午夜亚洲福利在线播放| av免费在线观看网站| 两性夫妻黄色片| 自线自在国产av| www国产在线视频色| 国产主播在线观看一区二区| 久久影院123| av有码第一页| 免费在线观看视频国产中文字幕亚洲| 久久久久久久久久久久大奶| 日本黄色视频三级网站网址| 免费在线观看黄色视频的| 午夜福利18| 日本三级黄在线观看| 中文字幕人妻丝袜一区二区| 在线国产一区二区在线| 大型av网站在线播放| 精品国产亚洲在线| 精品无人区乱码1区二区| 在线观看一区二区三区| 精品久久久久久久久久免费视频| 亚洲中文av在线| 欧美 亚洲 国产 日韩一| 国产成人精品无人区| 久久午夜综合久久蜜桃| 12—13女人毛片做爰片一| 精品熟女少妇八av免费久了| 最好的美女福利视频网| 国产亚洲欧美精品永久| 天堂动漫精品| 久久人人爽av亚洲精品天堂| 日韩三级视频一区二区三区| 色综合婷婷激情| 精品免费久久久久久久清纯| 亚洲午夜精品一区,二区,三区| 久久久久久久精品吃奶| 国产精品免费视频内射| 黄片大片在线免费观看| 男人舔女人的私密视频| 91麻豆av在线| 黄色女人牲交| 黑人欧美特级aaaaaa片| 黄色视频不卡| 香蕉久久夜色| 啪啪无遮挡十八禁网站| 夜夜看夜夜爽夜夜摸|