• <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
    不同品種大花萱草抗寒性比較
    国产一区二区在线观看日韩| 国产精品一区二区三区四区免费观看| 午夜老司机福利剧场| 亚洲成av人片在线播放无| 亚洲av不卡在线观看| 美女xxoo啪啪120秒动态图| 少妇高潮的动态图| 一进一出抽搐动态| 国产高清激情床上av| 亚洲人成网站在线播放欧美日韩| 99久国产av精品| 51国产日韩欧美| 婷婷亚洲欧美| 性色avwww在线观看| 一区二区三区四区激情视频 | 97热精品久久久久久| 黄色视频,在线免费观看| or卡值多少钱| 一区福利在线观看| 日本成人三级电影网站| 欧美xxxx性猛交bbbb| 精品久久久久久久久亚洲| 亚洲无线观看免费| 女人被狂操c到高潮| 国产精品爽爽va在线观看网站| 26uuu在线亚洲综合色| 日韩欧美 国产精品| 少妇熟女aⅴ在线视频| 丝袜喷水一区| 99九九线精品视频在线观看视频| 久久99蜜桃精品久久| 看十八女毛片水多多多| 十八禁国产超污无遮挡网站| 亚洲va在线va天堂va国产| 中文精品一卡2卡3卡4更新| 精品人妻视频免费看| 久久久久久久久久久免费av| 欧美高清成人免费视频www| 国产91av在线免费观看| 免费看av在线观看网站| 国产麻豆成人av免费视频| 国产人妻一区二区三区在| 乱系列少妇在线播放| av在线播放精品| 免费不卡的大黄色大毛片视频在线观看 | 亚洲图色成人| 亚洲精品久久久久久婷婷小说 | 一本一本综合久久| 天天躁夜夜躁狠狠久久av| 欧美成人a在线观看| 欧美性感艳星| 一个人观看的视频www高清免费观看| 免费不卡的大黄色大毛片视频在线观看 | 国产精品一区二区在线观看99 | 久久久久国产网址| 联通29元200g的流量卡| 国产精品久久电影中文字幕| 高清在线视频一区二区三区 | 丝袜美腿在线中文| 卡戴珊不雅视频在线播放| 岛国在线免费视频观看| 69人妻影院| 国产精品久久久久久av不卡| 91久久精品国产一区二区成人| 成年女人永久免费观看视频| 久久久久久久久久黄片| eeuss影院久久| 神马国产精品三级电影在线观看| 亚洲经典国产精华液单| 国产熟女欧美一区二区| 两性午夜刺激爽爽歪歪视频在线观看| 国产91av在线免费观看| 亚洲成a人片在线一区二区| 嫩草影院精品99| 亚洲成人精品中文字幕电影| 国产极品精品免费视频能看的| 伊人久久精品亚洲午夜| 亚洲av中文字字幕乱码综合| 高清在线视频一区二区三区 | 免费看a级黄色片| 中文亚洲av片在线观看爽| 在线a可以看的网站| 搞女人的毛片| 亚洲国产精品国产精品| 1024手机看黄色片| 好男人在线观看高清免费视频| 欧美成人a在线观看| 午夜福利在线观看吧| 性色avwww在线观看| 亚洲精品久久国产高清桃花| 春色校园在线视频观看| 国产精品国产三级国产av玫瑰| 亚洲成人久久性| 成人高潮视频无遮挡免费网站| 亚洲国产日韩欧美精品在线观看| 啦啦啦观看免费观看视频高清| 精品久久久久久久人妻蜜臀av| 午夜精品国产一区二区电影 | 亚洲国产精品国产精品| 亚洲久久久久久中文字幕| 亚洲av电影不卡..在线观看| 成熟少妇高潮喷水视频| 国产乱人偷精品视频| 一级毛片我不卡| 搡老妇女老女人老熟妇| 亚洲自拍偷在线| a级一级毛片免费在线观看| 高清毛片免费看| 51国产日韩欧美| 淫秽高清视频在线观看| 深夜a级毛片| 日本-黄色视频高清免费观看| 男人舔奶头视频| 精品不卡国产一区二区三区| 精品熟女少妇av免费看| 成人毛片60女人毛片免费| 国产真实伦视频高清在线观看| 一区福利在线观看| 久久精品国产清高在天天线| 久久精品国产亚洲av香蕉五月| 少妇的逼好多水| 天堂√8在线中文| 久久精品国产亚洲av香蕉五月| 亚洲经典国产精华液单| 简卡轻食公司| 国产淫片久久久久久久久| 亚洲精品色激情综合| 日韩一区二区三区影片| 国产精品久久视频播放| 99热全是精品| 老师上课跳d突然被开到最大视频| 国产亚洲精品av在线| 国产精品电影一区二区三区| 最近最新中文字幕大全电影3| 99久久人妻综合| 美女高潮的动态| 能在线免费观看的黄片| 中文字幕精品亚洲无线码一区| 插阴视频在线观看视频| 成人漫画全彩无遮挡| 亚洲精品粉嫩美女一区| 床上黄色一级片| 黑人高潮一二区| 国产精品女同一区二区软件| 欧美成人一区二区免费高清观看| 老司机影院成人| 亚洲欧美精品综合久久99| 成人亚洲精品av一区二区| 国产中年淑女户外野战色| 成年免费大片在线观看| 简卡轻食公司| 久久久久九九精品影院| 深爱激情五月婷婷| 久久综合国产亚洲精品| 久久综合国产亚洲精品| 能在线免费看毛片的网站| 久99久视频精品免费| 美女内射精品一级片tv| 在线观看午夜福利视频| 在现免费观看毛片| 热99re8久久精品国产| 成人鲁丝片一二三区免费| 久久精品国产鲁丝片午夜精品| 免费人成在线观看视频色| 一个人看视频在线观看www免费| h日本视频在线播放| 一进一出抽搐gif免费好疼| 色噜噜av男人的天堂激情| 99在线人妻在线中文字幕| 亚洲精品日韩在线中文字幕 | 亚洲av中文av极速乱| 成人毛片60女人毛片免费| 国产毛片a区久久久久| 精品久久久久久久久久免费视频| 老女人水多毛片| 国产精品一二三区在线看| 国产亚洲精品久久久com| 插阴视频在线观看视频| 久久久久网色| 成人性生交大片免费视频hd| 国产精品久久久久久久久免| 国内精品久久久久精免费| 特级一级黄色大片| 寂寞人妻少妇视频99o| 亚洲乱码一区二区免费版| 看非洲黑人一级黄片| 亚洲精品粉嫩美女一区| 一区二区三区四区激情视频 | 成人综合一区亚洲| av在线亚洲专区| 不卡视频在线观看欧美| 日韩精品青青久久久久久| 日韩欧美国产在线观看| 亚洲欧美日韩高清在线视频| 麻豆一二三区av精品| 久久精品国产亚洲av涩爱 | 丰满人妻一区二区三区视频av| 深夜精品福利| 欧美精品国产亚洲| 亚洲国产欧美人成| 变态另类丝袜制服| 99热只有精品国产| 亚洲最大成人av| 色播亚洲综合网| 精品人妻熟女av久视频| 内地一区二区视频在线| 精品久久久久久久人妻蜜臀av| 国产熟女欧美一区二区| 黄片无遮挡物在线观看| 国产精品免费一区二区三区在线| 18禁黄网站禁片免费观看直播| 日韩欧美精品v在线| 欧美日韩在线观看h| www.色视频.com| 最近视频中文字幕2019在线8| 成人无遮挡网站| 国产精华一区二区三区| 亚洲成a人片在线一区二区| 老师上课跳d突然被开到最大视频| 国产精品不卡视频一区二区| 国产精品一区二区在线观看99 | 国产午夜精品久久久久久一区二区三区| eeuss影院久久| 欧美一区二区国产精品久久精品| 欧美性感艳星| 国产人妻一区二区三区在| av在线老鸭窝| 亚洲精品456在线播放app| 亚洲不卡免费看| 美女高潮的动态| 成人欧美大片| 真实男女啪啪啪动态图| 欧美日本视频| 小蜜桃在线观看免费完整版高清| 美女被艹到高潮喷水动态| 亚洲欧美成人精品一区二区| 亚洲在线观看片| 午夜视频国产福利| 国产精品人妻久久久久久| 99视频精品全部免费 在线| 少妇被粗大猛烈的视频| 国产伦精品一区二区三区视频9| 一级毛片电影观看 | 国产三级在线视频| 久久久久久久午夜电影| 日本欧美国产在线视频| 亚洲av电影不卡..在线观看| a级一级毛片免费在线观看| 精品不卡国产一区二区三区| 国产一区二区三区av在线 | 国内久久婷婷六月综合欲色啪| 国产毛片a区久久久久| 久久精品91蜜桃| 永久网站在线| 日本黄大片高清| 2021天堂中文幕一二区在线观| 亚洲av中文av极速乱| 久久99蜜桃精品久久| av免费在线看不卡| 色综合站精品国产| 岛国在线免费视频观看| 成人性生交大片免费视频hd| 午夜精品国产一区二区电影 | 午夜精品在线福利| 好男人视频免费观看在线| 欧美色视频一区免费| 少妇熟女欧美另类| 久久精品国产鲁丝片午夜精品| 亚洲天堂国产精品一区在线| 成人综合一区亚洲| 国产白丝娇喘喷水9色精品| 国产极品天堂在线| 国产成人福利小说| 美女cb高潮喷水在线观看| 欧美高清性xxxxhd video| 看免费成人av毛片| av国产免费在线观看| 日韩成人av中文字幕在线观看| 99久国产av精品| 欧美日韩国产亚洲二区| 欧美区成人在线视频| 人人妻人人澡欧美一区二区| 久久精品国产自在天天线| 国产一区二区亚洲精品在线观看| 久久午夜福利片| 人妻久久中文字幕网| 日韩国内少妇激情av| 欧美+亚洲+日韩+国产| 亚洲综合色惰| 国产精品女同一区二区软件| 麻豆av噜噜一区二区三区| 亚洲精华国产精华液的使用体验 | 亚洲国产欧美人成| 男人的好看免费观看在线视频| 精品99又大又爽又粗少妇毛片| 日韩,欧美,国产一区二区三区 | 久久久久久九九精品二区国产| 亚洲自偷自拍三级| 春色校园在线视频观看| 乱码一卡2卡4卡精品| 精品久久久久久久久久久久久| 亚洲人成网站高清观看| 亚洲人成网站高清观看| 搞女人的毛片| 亚洲色图av天堂| 日韩欧美精品v在线| 男人舔女人下体高潮全视频| 26uuu在线亚洲综合色| 久久中文看片网| 色综合站精品国产| 最后的刺客免费高清国语| 99热这里只有精品一区| 欧美最黄视频在线播放免费| 精品不卡国产一区二区三区| 精品一区二区免费观看| 男女视频在线观看网站免费| 亚洲天堂国产精品一区在线| 国产精品久久久久久久电影| 免费不卡的大黄色大毛片视频在线观看 | 久久久久九九精品影院| 久久久精品大字幕| .国产精品久久| 99久久久亚洲精品蜜臀av| 国产私拍福利视频在线观看| 欧美成人a在线观看| 成人av在线播放网站| 一级av片app| 久久久精品大字幕| 麻豆久久精品国产亚洲av| 久久精品国产亚洲av香蕉五月| 免费大片18禁| 非洲黑人性xxxx精品又粗又长| 国产av麻豆久久久久久久| 麻豆久久精品国产亚洲av| 国产精品免费一区二区三区在线| 欧美成人一区二区免费高清观看| 欧美高清性xxxxhd video| 少妇熟女aⅴ在线视频| 97人妻精品一区二区三区麻豆| 噜噜噜噜噜久久久久久91| 美女cb高潮喷水在线观看| 午夜久久久久精精品| 欧美最新免费一区二区三区| avwww免费| 免费搜索国产男女视频| 晚上一个人看的免费电影| 国产成人午夜福利电影在线观看| 亚洲自拍偷在线| 久久韩国三级中文字幕| 亚洲国产精品国产精品| 国产精品.久久久| 免费观看a级毛片全部| 嫩草影院入口| 精品久久久久久久久久免费视频| 日韩欧美 国产精品| 国产精品爽爽va在线观看网站| 色视频www国产| 国产伦在线观看视频一区| 亚洲成人精品中文字幕电影| 综合色av麻豆| 国产老妇伦熟女老妇高清| 精品人妻一区二区三区麻豆| 国产精品99久久久久久久久| 欧美三级亚洲精品| 婷婷六月久久综合丁香| 国产精品嫩草影院av在线观看| 国产av麻豆久久久久久久| 22中文网久久字幕| 国产精品日韩av在线免费观看| 午夜激情欧美在线| 国产伦一二天堂av在线观看| 久久亚洲精品不卡| 午夜精品在线福利| 日日摸夜夜添夜夜添av毛片| 国内揄拍国产精品人妻在线| 日本欧美国产在线视频| 国产在线男女| 男女下面进入的视频免费午夜| 午夜免费激情av| 欧美高清成人免费视频www| av国产免费在线观看| 男女视频在线观看网站免费| 精品久久久噜噜| 欧美+亚洲+日韩+国产| 国产高清有码在线观看视频| 成人二区视频| 国内精品久久久久精免费| 亚洲欧美精品自产自拍| 成人国产麻豆网| 一个人观看的视频www高清免费观看| 免费在线观看成人毛片| 欧美高清成人免费视频www| 成年女人看的毛片在线观看| 久久这里只有精品中国| 日本免费a在线| 精品日产1卡2卡| 免费av观看视频| 波多野结衣高清无吗| 美女xxoo啪啪120秒动态图| 久久99热6这里只有精品| 黄色一级大片看看| 国产蜜桃级精品一区二区三区| 欧美色视频一区免费| 特级一级黄色大片| 亚洲av男天堂| 51国产日韩欧美| 久久久久久久久久久免费av| 麻豆国产97在线/欧美| 如何舔出高潮| 99视频精品全部免费 在线| 国国产精品蜜臀av免费| 美女高潮的动态| 免费av不卡在线播放| 亚洲国产欧美在线一区| 三级毛片av免费| 欧美高清成人免费视频www| 国产三级中文精品| 亚洲精品色激情综合| 欧美一区二区亚洲| 中文字幕精品亚洲无线码一区| 欧美潮喷喷水| 此物有八面人人有两片| 久久综合国产亚洲精品| 午夜精品国产一区二区电影 | 国产黄色视频一区二区在线观看 | 99久久精品一区二区三区| 国产老妇女一区| 成人国产麻豆网| 午夜精品在线福利| 亚洲欧美成人精品一区二区| 26uuu在线亚洲综合色| 高清毛片免费观看视频网站| 久久久久久久久久久丰满| .国产精品久久| 亚州av有码| 一进一出抽搐动态| 在线观看66精品国产| 久久久久久久午夜电影| 99久久人妻综合| 日本五十路高清| 亚洲精品久久国产高清桃花| 精品久久久久久久久久免费视频| 2021天堂中文幕一二区在线观| 亚洲欧美清纯卡通| 国产黄片视频在线免费观看| 噜噜噜噜噜久久久久久91| 99久久成人亚洲精品观看| 日韩人妻高清精品专区| 日日摸夜夜添夜夜爱| 国产伦精品一区二区三区四那| 午夜a级毛片| 三级毛片av免费| 中文亚洲av片在线观看爽| 一个人免费在线观看电影| 亚洲人成网站在线播| av天堂中文字幕网| 国产一区二区三区在线臀色熟女| 国产日本99.免费观看| 久久久久久久久久成人| 搡女人真爽免费视频火全软件| 欧美xxxx黑人xx丫x性爽| 成人永久免费在线观看视频| 婷婷色综合大香蕉| 少妇猛男粗大的猛烈进出视频 | 久久人人爽人人片av| ponron亚洲| 一本久久精品| 亚洲在久久综合| 99热网站在线观看| 狂野欧美白嫩少妇大欣赏| 国产老妇女一区| 噜噜噜噜噜久久久久久91| 直男gayav资源| 国产亚洲5aaaaa淫片| 久久午夜亚洲精品久久| 少妇熟女aⅴ在线视频| 九九热线精品视视频播放| 精品99又大又爽又粗少妇毛片| 波多野结衣高清无吗| 国产精华一区二区三区| 欧美日本视频| 国产精品一区二区三区四区久久| 亚洲久久久久久中文字幕| 在线观看av片永久免费下载| 国产精品久久视频播放| 精品久久久久久久人妻蜜臀av| 国产成人91sexporn| 如何舔出高潮| 美女大奶头视频| 性插视频无遮挡在线免费观看| 国产在线精品亚洲第一网站| 欧美+亚洲+日韩+国产| 天美传媒精品一区二区| 成人综合一区亚洲| 日韩中字成人| 老司机影院成人| 白带黄色成豆腐渣| 嫩草影院精品99| 在线播放国产精品三级| 日本五十路高清| 亚洲欧美日韩高清专用| 老熟妇乱子伦视频在线观看| 久久久久久大精品| 国产色爽女视频免费观看| eeuss影院久久| 韩国av在线不卡| 不卡一级毛片| 亚洲婷婷狠狠爱综合网| 亚洲性久久影院| 亚洲av成人av| 亚洲久久久久久中文字幕| 欧美+亚洲+日韩+国产| 毛片女人毛片| 大香蕉久久网| 我要看日韩黄色一级片| 美女被艹到高潮喷水动态| 在线观看一区二区三区| 国产精品美女特级片免费视频播放器| 搞女人的毛片| 亚洲精品粉嫩美女一区| 久久草成人影院| 网址你懂的国产日韩在线| 美女xxoo啪啪120秒动态图| 久久久a久久爽久久v久久| 大型黄色视频在线免费观看| 性插视频无遮挡在线免费观看| 国产成人一区二区在线| 综合色av麻豆| 欧美最黄视频在线播放免费| 黑人高潮一二区| 日韩欧美国产在线观看| 成人毛片a级毛片在线播放| 女人十人毛片免费观看3o分钟| 亚洲精华国产精华液的使用体验 | 国产黄色小视频在线观看| 精品久久久久久久久久久久久| 国产 一区精品| 桃色一区二区三区在线观看| 久久热精品热| 1024手机看黄色片| 日本黄色视频三级网站网址| 天天躁日日操中文字幕| 国产精品不卡视频一区二区| 亚洲无线在线观看| 99久久精品一区二区三区| 少妇的逼好多水| 免费观看精品视频网站| 日韩三级伦理在线观看| 老女人水多毛片| 亚洲一级一片aⅴ在线观看| 国产久久久一区二区三区| 精品久久久久久久久亚洲| 久久草成人影院| 久久人人精品亚洲av| 国产亚洲欧美98| 日日啪夜夜撸| 一级黄片播放器| 久久99精品国语久久久| 日韩一区二区三区影片| 久久久午夜欧美精品| 久久精品夜色国产| 男人舔奶头视频| 桃色一区二区三区在线观看| 国产午夜福利久久久久久| 久久韩国三级中文字幕| 老熟妇乱子伦视频在线观看| 欧美性感艳星| 麻豆一二三区av精品| 99在线人妻在线中文字幕| 国产大屁股一区二区在线视频| 黄片无遮挡物在线观看| 美女高潮的动态| 中文字幕av成人在线电影| 美女cb高潮喷水在线观看| 一级黄色大片毛片| 丝袜美腿在线中文| 欧美日本视频| 中文字幕熟女人妻在线| 精品久久久噜噜| 色吧在线观看| 国产又黄又爽又无遮挡在线| 美女大奶头视频| 国产中年淑女户外野战色| 亚洲精品亚洲一区二区| 亚洲成a人片在线一区二区| 亚洲欧美精品专区久久| 91精品国产九色| 午夜a级毛片| 级片在线观看| 丰满人妻一区二区三区视频av| 久久久精品大字幕| 欧美精品国产亚洲| 国产激情偷乱视频一区二区| 久久人人爽人人爽人人片va| 亚洲丝袜综合中文字幕| 有码 亚洲区| 人人妻人人看人人澡| 国产精品精品国产色婷婷| 内射极品少妇av片p| 国产高清有码在线观看视频| 国产精品.久久久| 日韩欧美精品免费久久| 九九在线视频观看精品| 好男人在线观看高清免费视频| 精品人妻偷拍中文字幕| 韩国av在线不卡| 亚洲成人精品中文字幕电影| 久久久久性生活片| 欧美日本亚洲视频在线播放| 国产单亲对白刺激| 久久精品久久久久久噜噜老黄 | 亚洲精品国产av成人精品| 欧美三级亚洲精品| 亚洲最大成人手机在线| 国产爱豆传媒在线观看|