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

    干旱脅迫和復(fù)水處理后梭梭轉(zhuǎn)錄因子的轉(zhuǎn)錄組分析

    2020-04-16 03:43:54何江峰王力偉房永雨王蘊華劉紅葵
    華北農(nóng)學(xué)報 2020年1期
    關(guān)鍵詞:梭梭家族測序

    何江峰,王力偉,2,房永雨,王蘊華,王 朝,劉紅葵

    (1.內(nèi)蒙古自治區(qū)農(nóng)牧業(yè)科學(xué)院 生物技術(shù)研究中心,內(nèi)蒙古 呼和浩特 010031;2.內(nèi)蒙古大學(xué) 生命科學(xué)學(xué)院,內(nèi)蒙古 呼和浩特 010020;3.內(nèi)蒙古農(nóng)業(yè)大學(xué) 生命科學(xué)學(xué)院,內(nèi)蒙古 呼和浩特 010010)

    梭梭(Haloxylonammodendron)屬黎科(Chenopodiaceae),為超旱生小喬木,呈高大灌木狀,植株高度一般在1~4 m,個別植株可高達8~10 m,亦稱鹽木、瑣瑣樹、梭梭柴。梭梭是肉蓯蓉(CistanchedeserticolaY.C.Ma)的寄主,肉蓯蓉是名貴的中藥材,在醫(yī)學(xué)上具有抗腫瘤、保肝、補腎、通便、抗老年癡呆、抗輻射等獨特功能,其寄生在梭梭的根端[1-2]。梭梭長期生長在荒漠環(huán)境中,形成了抗旱、抗風沙、耐高溫、耐寒、耐鹽堿和耐貧瘠等特性。因此,梭梭是荒漠化地區(qū)主要的防風固沙樹種,也可作為駱駝、羊適口性飼料,是荒漠草場中主要的木本飼料樹種之一[3]。梭梭具有重要的經(jīng)濟價值、生態(tài)價值和飼用價值,是我國的二級保護植物,國內(nèi)外許多學(xué)者已從梭梭抗逆的形態(tài)學(xué)、解剖學(xué)、生理學(xué)、生態(tài)學(xué)及分子生物學(xué)等方面進行了相關(guān)研究。

    梭梭的生存環(huán)境極其嚴酷,其主要分布地區(qū)常年干旱少雨,年降水量不足50 mm,相對濕度低于40%,蒸發(fā)量高于降水量的30~100倍,生存環(huán)境的極端溫度最高可達40 ℃以上,最低為-40 ℃以下,地表溫度可達70 ℃以上[4]。為了適應(yīng)極端環(huán)境,梭梭逐步形成了自身的極端抗旱機制,但梭梭的抗旱機理、抗旱遺傳機制和抗旱相關(guān)組學(xué)的研究相對較少,前人研究發(fā)現(xiàn),部分轉(zhuǎn)錄因子及與干旱脅迫相關(guān),例如,AP2-EREBP轉(zhuǎn)錄因子(TINY、CBF1、Ptis、AtEBP、DREB1、DREB2)主要調(diào)節(jié)植物對低溫、干旱及高鹽[5]等的分子應(yīng)答反應(yīng);bZIP轉(zhuǎn)錄因子在擬南芥基因組、豆科植物基因組、水稻基因組中參與對各種逆境脅迫的響應(yīng)[6-7];擬南芥中3個NAC基因(ANAC019、ANAC055和ANAC072)的過表達能夠顯著提高植株對干旱的抗性[8];許多NAC基因過表達能夠顯著提高水稻植株的脅迫抗性,卻不會減緩植株的正常生長[9-10];bHLH基因也參與水稻應(yīng)對干旱脅迫的分子調(diào)控過程[11]。在改良植物的抗逆分子育種中,導(dǎo)入或改良轉(zhuǎn)錄因子較導(dǎo)入或改良個別功能基因在增強植物抗逆性方面更為有效。因此,本研究主要對干旱脅迫和復(fù)水處理下梭梭的轉(zhuǎn)錄因子進行分析,挖掘與干旱相關(guān)的轉(zhuǎn)錄因子,為闡明梭梭抗旱機理奠定理論基礎(chǔ),保護及開發(fā)我國抗逆植物基因資源。

    1 材料和方法

    1.1 材料及處理條件

    梭梭種子收集于內(nèi)蒙古自治區(qū)阿拉善左旗。將梭梭種植在塑料花盆中(直徑:15 cm,高:18 cm),盆內(nèi)土壤均勻,濕度及質(zhì)量均一致。梭梭幼苗在植物培養(yǎng)室種植,室內(nèi)溫度為24/20 ℃(16 h光照 / 8 h 暗處理),相對濕度為30%,光子通量密度為350 μmol/(m2·s)。本研究隨機單株種植25盆梭梭幼苗,包括1個正常澆水組(對照組,CK)、干旱脅迫組(土壤相對含水量為30%~45%,LMS;土壤相對含水量為20%~25%,MWS;土壤相對含水量為10%~15%,SWS)和旱后復(fù)水處理組(土壤相對含水量為在重度水分脅迫組上增加10%~15%,RSWS)。干旱脅迫組、復(fù)水處理組和對照組的種植均設(shè)置5個生物學(xué)重復(fù),每天在固定時間測質(zhì)量、控水、補水并記錄待幼苗生長90 d之后,分別從CK、LWS、MWS、SWS 以及重度脅迫90 d之后復(fù)水處理7 d的RSWS的植株上采集等量的根、莖和葉混合,液氮速凍,置于-80 ℃冰箱中儲存?zhèn)溆茫崛「魈幚淼目俁NA。

    1.2 總 RNA 提取

    材料總RNA的提取采用 RNAprep Pure Plant Plus KIT(DP441,TIANGEN北京生化科技有限公司)試劑盒。用 DNase酶(2212,大連寶生物工程有限公司)對樣品中的 gDNA 進行消化處理。RNA完整性和純度質(zhì)量經(jīng)1.2%非變性瓊脂糖凝膠電泳和 Agilent 2100 分析系統(tǒng)檢測。

    1.3 RNA測序和分析方法

    總RNA提取按照文獻[12]的方法,利用 Solexa/Illumina進行測序及深圳華大基因公司的Illumina系統(tǒng)軟件進行數(shù)據(jù)分析。RPKM值用于基因表達水平的分析,利用Cluster軟件[13],以歐氏距離為距離矩陣計算公式,對基因和試驗條件同時進行等級聚類分析,聚類結(jié)果用Java Treeview顯示,利用WEGO軟件對差異基因做GO功能分類統(tǒng)計[14],pathway顯著性富集分析以KEGG pathway為單位,應(yīng)用超幾何檢驗,找出與整個基因組相比后差異表達基因中顯著性富集的pathway基因和基因組的代謝途徑[15]。

    1.4 de novo 拼接及序列分析

    用Trinity[16]做轉(zhuǎn)錄組 de novo 組裝,并通過 reads overlap 關(guān)系得到 Contig。然后,將 clean reads 比對回Contig,通過paired-end reads 能確定來自同一轉(zhuǎn)錄本的不同 Contig 以及這些 Contig 之間的距離,Trinity 將 Contig 連在一起,最后得到兩端不能再延長的序列,即為 Unigene。組裝得到的 Unigene,首先使用 Tgicl 將其去冗余并進一步拼接,然后再對這些序列進行同源轉(zhuǎn)錄本聚類,最終得到非冗余 Unigene。用 Tophat v2.0.9(Broad Institute)將拼接得到的非冗余 Unigene 與已公布的同種屬基因組信息和轉(zhuǎn)錄組數(shù)據(jù)進行比對、注釋,用已獲取各測序樣本中包含的轉(zhuǎn)錄因子編碼基因信息。

    1.5 差異表達基因的篩選

    樣本中 Unigene 的表達量以公式化校正的 FPKM(Fragments per kb per million fragments)[17]值衡量。2個樣本間差異表達基因(Differential expression genes,DEGs)的篩選及分析使用DEGSeq R package軟件,以FDR(False discovery rate)≤0.001,且Log2(fold change,F(xiàn)C)≥1作為篩選閾值進行差異表達基因的篩選。

    1.6 差異表達基因 Real-time PCR 驗證

    從測序結(jié)果中隨機挑選 9 個轉(zhuǎn)錄因子編碼基因,用熒光實時定量 PCR 法驗證轉(zhuǎn)錄組測序數(shù)據(jù)的有效性和準確性。引物設(shè)計用 Primer 5 軟件,序列見表1。以18S為內(nèi)參基因,用2-ΔΔCt計算,并校正目標基因在不同處理下的表達量,每個基因3次生物學(xué)重復(fù)。10 μL反應(yīng)體系中包括5 μL SYBR Premix Ex Taq(TaKaRa,TaKaRa Clontech,Japan),5 μmol/L上游引物,5 μmol/L下游引物和20 ng 模板cDNA。所有反應(yīng)體系加樣都在96孔板進行,并在Lightcycle-480(羅氏)熒光定量PCR儀上完成,反應(yīng)程序:95 ℃ 5 min;后按95 ℃ 15 s,55 ℃ 2 min,72 ℃ 1 min,共40個循環(huán);72 ℃延伸2 min。18S基因為擴增體系的內(nèi)參基因,擴增產(chǎn)物的表達水平按 2-ΔΔCt方法計算[18]。每個樣品都有3次技術(shù)重復(fù),結(jié)果都以平均值±標準誤(s)(n= 3)計算。

    表1 實時定量PCR引物

    2 結(jié)果與分析

    2.1 轉(zhuǎn)錄組數(shù)據(jù)統(tǒng)計分析

    以梭梭幼苗為材料,對CK、LWS、MWS、SWS和RSWS處理組進行轉(zhuǎn)錄組測序和表達譜分析。對測序數(shù)據(jù)進行拼接、比對和整理,經(jīng)de novo 拼接、去冗余、同源轉(zhuǎn)錄本聚類后,最終在CK、LWS、MWS、SWS和RSWS 5個組中得到74 641條非冗余Unigene,平均長度780 bp,N50長度1 537 bp。 將CK、LWS、MWS、SWS和RSWS混合樣本拼接得到的非冗余Unigene進行比對、注釋,共檢測到56個轉(zhuǎn)錄因子家族,1 307個轉(zhuǎn)錄因子編碼基因;在CK、LWS、MWS、SWS和RSWS組分別檢索到 934,978,994,1 003,991 條轉(zhuǎn)錄因子編碼基因。LWS、MWS、SWS分別與CK組相比,RSWS與SWS相比,差異表達轉(zhuǎn)錄因子均以下調(diào)表達方式為主,CK-vs-SWS組中下調(diào)轉(zhuǎn)錄因子數(shù)目最多,達112條,約占該組差異表達轉(zhuǎn)錄因子編碼基因總數(shù)的59.3%(表2)。

    表2 差異表達轉(zhuǎn)錄因子

    2.2 差異表達轉(zhuǎn)錄因子家族統(tǒng)計分析

    獲得4個對比組轉(zhuǎn)錄因子基本數(shù)量信息后,進一步統(tǒng)計分析轉(zhuǎn)錄因子家族及其包含成員的數(shù)量關(guān)系。LWS、MWS、SWS分別與CK組對比,RSWS與SWS相比,表3分析表明,其中,AP2-EREBP、WRKY、bHLH、MYB、NAC和ABI3VP1 家族富集的DEGs 超過該組差異表達基因總數(shù)的50%以上,其中AP2-EREBP、WRKY、MYB和bHLH家族轉(zhuǎn)錄因子較多。由圖1-A-C結(jié)果顯示,AP2-EREBP家族包含的下調(diào)差異表達基因,分別為6個(CK-vs-LWS)、17個(CK-vs-MWS)和 16個(CK-vs-SWS)DEGs,約占各組DEGs總數(shù)的75%,81%,59%;WRKY家族包含的下調(diào)差異表達基因,分別為6個(CK-vs-LWS)、12個(CKvs-MWS)和 2個(CK-vs-SWS)DEGs,約占各組DEGs總數(shù)的100%,100%,33%;bHLH家族包含的下調(diào)差異表達基因,分別為5個(CK-vs-LWS)、4個(CK-vs-MWS)和 16 個(CK-vs-SWS)DEGs,約占各組DEGs總數(shù)的100%,57%,84%;MYB家族包含的下調(diào)差異表達基因,分別為5個(CK-vs-LWS)、9個(CK-vs-MWS)和 10 個(CK-vs-SWS),約占各組DEGs總數(shù)的100%,64%,45%。圖1-D結(jié)果顯示,AP2-EREBP家族、NAC家族和WRKY家族包含的上調(diào)差異表達基因分別占各自DEGs總數(shù)的67%,86%,75%。以上表達規(guī)律顯示,MYB和WRKY家族基因在LWS組中沒有上調(diào)表達,但在MWS和SWS中其上調(diào)表達數(shù)逐漸增加。綜合說明,RSWS組中偏向通過增加 AP2-EREBP、NAC家族和WRKY家族編碼基因的上調(diào)表達數(shù)來實現(xiàn)梭梭幼苗對水分刺激的應(yīng)答調(diào)節(jié);而干旱組(LWS、 MWS和 SWS)更偏向于下調(diào)AP2-EREBP和bHLH 的DEGs數(shù)量來實現(xiàn)其對干旱脅迫的響應(yīng)。

    表3 干旱和復(fù)水組各轉(zhuǎn)錄因子家族中的基因表達數(shù)

    2.3 差異表達轉(zhuǎn)錄因子表達量分析

    從4個對比組中挑選了表達量(FC)變化幅度最大的10個上調(diào)和下調(diào)基因,并對這些基因的表達規(guī)律進行統(tǒng)計分析。結(jié)果表明:差異表達基因Tify(CL5087)和NAC(Unigene28108)在干旱脅迫組(LWS、MWS和SWS)和復(fù)水組(RSWS)中均上調(diào)表達,表達量(FC)變幅為1.23~12.24,2.89~11.37;AP2-EREBP(Unigene2016)、AP2-EREBP(Unigene442)和AP2-EREBP(Unigene160)在干旱脅迫組(LWS、 MWS和SWS)中下調(diào),在復(fù)水組(RSWS)中又上調(diào)表達。表達量變幅為-4.65~12.22,-8.31~7.38,-6.99~6.85。bHLH(CL2529)和E2F-DP(CL4870)在干旱脅迫組(LWS、 MWS和SWS)中下調(diào)表達,變化幅度分別為-14.17~-3.73,-11.89~-4.56;在復(fù)水組中卻不表達。RWP-RK(CL8167)、DBP(CL6537),在干旱脅迫組(LWS、MWS和SWS)和復(fù)水組(RSWS)中均上調(diào)表達。BZIP(Unigene23554)在干旱脅迫組(LWS、 MWS和SWS)和復(fù)水組(RSWS)中均下調(diào)表達。綜合統(tǒng)計數(shù)據(jù)結(jié)果顯示,AP2-EREBP、NAC、E2F-DP、bHLH 和BZIP等5類轉(zhuǎn)錄因子對干旱和水分刺激應(yīng)答顯著(表 4)。

    A.CK-vs-LWS的差異表達基因; B.CK-vs-MWS的差異表達基因;C.CK-vs-SWS的差異表達基因; D.RSWS-vs-SWS的差異表達基因。

    表4 對比組差異表達轉(zhuǎn)錄因子表達量分析

    注:Log2(FC)≥1時,基因表達顯著。

    Note:Gene expression is significant when Log2(FC)≥1.

    2.4 Real-time PCR 驗證

    從測序數(shù)據(jù)中隨機選取了9個轉(zhuǎn)錄因子編碼基因,用Real-time PCR 方法驗證測序數(shù)據(jù)的有效性。圖2顯示,干旱和水分刺激應(yīng)答相關(guān)的AP2-EREBP、WRKY、bHLH、MYB、C2C2-Dof、RWP-RK、NAC、C2H2和ABI3VP1編碼基因在處理前后表達量變化均極顯著。干旱脅迫組和復(fù)水組測序數(shù)據(jù)與用熒光定量PCR 法得到的基因表達數(shù)據(jù)之間的相關(guān)性分別為0.947 8(P<0.01),0.967 2(P<0.01),說明干旱脅迫和復(fù)水處理后產(chǎn)生的差異表達基因數(shù)據(jù)有效。

    均值±s,樣本均為3個重復(fù);不同大寫字母表示統(tǒng)計結(jié)果在1%水平上有顯著差異。

    3 討論與結(jié)論

    本研究通過對干旱脅迫和復(fù)水處理下梭梭轉(zhuǎn)錄因子應(yīng)答的分子機理進行研究,發(fā)掘其中的關(guān)鍵轉(zhuǎn)錄因子,對培育抗旱能力強、水分利用效率高的新作物品種具有重要的意義。研究表明,轉(zhuǎn)錄因子(Transcription factors,TFs)受脅迫刺激誘導(dǎo)產(chǎn)生,并將信號傳遞、放大,能夠通過調(diào)控靶基因轉(zhuǎn)錄效率從多個層面降低脅迫對植物的傷害,對植物在逆境下的生長發(fā)育起到了至關(guān)重要的調(diào)節(jié)作用,尤其是AP2-EREBP、 E2F/DP、MYB、NAC、WRKY、bHLH和BZIP等家族,是信號轉(zhuǎn)導(dǎo)通路上的關(guān)鍵調(diào)控因子,對干旱脅迫響應(yīng)和水分刺激應(yīng)答發(fā)揮了重要的調(diào)節(jié)作用[19-23]。本試驗測序結(jié)果中,各對比組中 AP2-EREBP、MYB、bHLH、NAC、WRKY和ABI3VP1等轉(zhuǎn)錄因子編碼基因的數(shù)量占主要比例,說明這幾類轉(zhuǎn)錄因子對梭梭干旱脅迫和水分刺激應(yīng)答調(diào)控作用極顯著,其中,AP2-EREBP、WRKY、MYB和bHLH是所有56個家族中差異表達轉(zhuǎn)錄因子成員最多的蛋白家族。AP2-EREBP轉(zhuǎn)錄因子要調(diào)節(jié)植物對低溫、干旱及高鹽[23]等的分子應(yīng)答反應(yīng)。本研究中 AP2-EREBP是干旱和復(fù)水處理下數(shù)量豐富、表達量變化最顯著的一類轉(zhuǎn)錄因子,且AP2-EREBP(Unigene2016)、AP2-EREBP(Unigene442)和AP2-EREBP(Unigene160),在干旱脅迫下顯著下調(diào),而在復(fù)水組中顯著上調(diào),與 bHLH 轉(zhuǎn)錄因子類似,干旱脅迫處理下誘導(dǎo)更多的下調(diào)ERF轉(zhuǎn)錄因子表達,說明這類轉(zhuǎn)錄因子發(fā)揮調(diào)控作用的模式可能與bHLH類似。本研究結(jié)果還表明,干旱處理和復(fù)水處理后下調(diào)的bHLH編碼基因數(shù)量均大于上調(diào)數(shù),推測bHLH編碼基因主要通過下調(diào)其表達調(diào)控機體應(yīng)答脅迫過程。張子佳等[11]對水稻 22 個bHLH基因在PEG 和 ABA 脅迫下的表達譜進行了分析,結(jié)果表明,其中一些bHLH基因在響應(yīng)干旱脅迫過程中受ABA負調(diào)控。此外,不同的bHLH基因參與應(yīng)答環(huán)境脅迫的調(diào)節(jié)方式不同,表明參與 PEG 和 ABA 脅迫應(yīng)答的bHLH基因具有不同的分子路徑或模式。因此,梭梭中bHLH轉(zhuǎn)錄因子的進一步研究已在該課題組將繼續(xù)進行。

    Liao等[24]從研究的大豆基因中獲得了 156 個 MYB 家族基因,利用酵母單雜交的方法篩選出了40個左右與逆境相關(guān)的基因,將這些基因轉(zhuǎn)入擬南芥進行功能驗證,發(fā)現(xiàn)轉(zhuǎn) GmMYB177 提高了擬南芥耐旱性。丁震乾等[25]在陸地棉中克隆了 MYB 轉(zhuǎn)錄因子基因GhRAX3,發(fā)現(xiàn)其表達在干旱脅迫 0.5 h 后即表現(xiàn)為顯著上調(diào),且 48 h 內(nèi)持續(xù)高表達水平;而抑制GhRAX3的表達,則加快了棉花植株的失水率、細胞質(zhì)膜過氧化和細胞受損程度,降低了棉花對干旱脅迫的耐受性。目前,發(fā)現(xiàn)的大多數(shù) MYB 類轉(zhuǎn)錄因子基因均表現(xiàn)為提高植物抗旱能力,但亦有例外,如 Zhou 等[26]的研究任務(wù),將JcMYB001轉(zhuǎn)入擬南芥后,其過表達反而增加了擬南芥對干旱的敏感性。本研究中的MYB轉(zhuǎn)錄因子在干旱脅迫組與對照組相比,其下調(diào)差異表達基因較其他轉(zhuǎn)錄因子多,MYB家族中上調(diào)基因數(shù)隨著干旱脅迫程度的增加而增加。說明本研究結(jié)果也證實了MYB家族的轉(zhuǎn)錄因子對干旱脅迫具有正調(diào)控和負調(diào)控性。

    綜合以上結(jié)果推測,發(fā)揮主要調(diào)控作用的轉(zhuǎn)錄因子家族或成員間可能存在一定的共表達性或相互關(guān)聯(lián)性,各轉(zhuǎn)錄因子間協(xié)同作用對干旱和水分刺激做出精細應(yīng)答調(diào)控。植物對干旱脅迫和水分刺激應(yīng)答過程受眾多分子和細胞通路的調(diào)控[27]。轉(zhuǎn)錄因子能特異與順式作用元件結(jié)合激活或抑制下游靶基因表達,參與脅迫信號傳遞網(wǎng)絡(luò),其“閥門作用”的重要性不言而喻。本研究借助高通量測序技術(shù),從整體層面概要性的分析了干旱和干旱后復(fù)水處理下梭梭幼苗轉(zhuǎn)錄因子種類及初步的表達規(guī)律,結(jié)果將為干旱脅迫應(yīng)答調(diào)控網(wǎng)絡(luò)、有目的地篩選關(guān)鍵調(diào)控基因和分子設(shè)計育種等研究提供基礎(chǔ)理論數(shù)據(jù)。

    猜你喜歡
    梭梭家族測序
    杰 Sir 帶你認識宏基因二代測序(mNGS)
    新民周刊(2022年27期)2022-08-01 07:04:49
    梭梭的建筑課
    哈哈畫報(2022年4期)2022-04-19 11:11:54
    二代測序協(xié)助診斷AIDS合并馬爾尼菲籃狀菌腦膜炎1例
    傳染病信息(2021年6期)2021-02-12 01:52:58
    與生命賽跑的“沙漠植被之王”——梭梭
    HK家族崛起
    沙漠梭梭的守望者
    《小偷家族》
    電影(2019年3期)2019-04-04 11:57:18
    皿字家族
    家族中的十大至尊寶
    環(huán)保志愿者在阿拉善種梭梭固沙
    綠色中國(2015年6期)2015-07-03 01:13:10
    不卡视频在线观看欧美| 免费人妻精品一区二区三区视频| 草草在线视频免费看| 91久久精品国产一区二区三区| 国产精品久久久久久av不卡| 天堂俺去俺来也www色官网| 伊人久久国产一区二区| 欧美激情极品国产一区二区三区 | 肉色欧美久久久久久久蜜桃| 国产伦精品一区二区三区视频9| 狂野欧美激情性xxxx在线观看| 插阴视频在线观看视频| 国产深夜福利视频在线观看| 在线精品无人区一区二区三| 久久久久久久久久久久大奶| 欧美日韩视频精品一区| 成人影院久久| 国产精品99久久99久久久不卡 | 有码 亚洲区| av卡一久久| 国产精品一二三区在线看| 亚洲av成人精品一区久久| 国产成人freesex在线| 人成视频在线观看免费观看| 一级片'在线观看视频| 啦啦啦视频在线资源免费观看| 日韩 亚洲 欧美在线| 精品人妻一区二区三区麻豆| 久久女婷五月综合色啪小说| 少妇被粗大猛烈的视频| 日本av手机在线免费观看| 国产成人aa在线观看| 日本91视频免费播放| 欧美激情极品国产一区二区三区 | 一级二级三级毛片免费看| 国产精品嫩草影院av在线观看| 亚洲不卡免费看| 亚洲久久久国产精品| 国产精品 国内视频| 日日撸夜夜添| 春色校园在线视频观看| 亚洲精品一区蜜桃| videosex国产| 亚洲国产av影院在线观看| 国产综合精华液| 国产成人aa在线观看| 免费少妇av软件| 波野结衣二区三区在线| 精品国产一区二区久久| 久久午夜福利片| 婷婷成人精品国产| 伊人久久精品亚洲午夜| av福利片在线| 久久 成人 亚洲| 天堂俺去俺来也www色官网| 久久这里有精品视频免费| 亚洲成人av在线免费| 欧美bdsm另类| 毛片一级片免费看久久久久| 国产深夜福利视频在线观看| 九九在线视频观看精品| 亚洲美女视频黄频| 精品国产国语对白av| 晚上一个人看的免费电影| 亚洲精品一二三| 国产永久视频网站| 亚洲精品国产色婷婷电影| 夜夜骑夜夜射夜夜干| 女的被弄到高潮叫床怎么办| 中国美白少妇内射xxxbb| 母亲3免费完整高清在线观看 | 国产成人精品无人区| 99热网站在线观看| 婷婷色综合www| 亚洲少妇的诱惑av| 日本免费在线观看一区| 亚洲精品色激情综合| 亚洲精品av麻豆狂野| 51国产日韩欧美| 成人综合一区亚洲| 交换朋友夫妻互换小说| 久久久久人妻精品一区果冻| 久久久精品区二区三区| 不卡视频在线观看欧美| 最黄视频免费看| 人人妻人人添人人爽欧美一区卜| 妹子高潮喷水视频| 欧美激情极品国产一区二区三区 | 中文字幕av电影在线播放| 欧美丝袜亚洲另类| 免费大片黄手机在线观看| av播播在线观看一区| 久久影院123| 97精品久久久久久久久久精品| 日韩精品免费视频一区二区三区 | 大香蕉久久成人网| 观看美女的网站| a级毛色黄片| 考比视频在线观看| 天天操日日干夜夜撸| 777米奇影视久久| 丰满少妇做爰视频| 欧美xxⅹ黑人| 国产成人91sexporn| 少妇被粗大猛烈的视频| 色婷婷av一区二区三区视频| 多毛熟女@视频| 国产黄片视频在线免费观看| 久久99热这里只频精品6学生| 亚洲一区二区三区欧美精品| www.av在线官网国产| 女人久久www免费人成看片| 精品熟女少妇av免费看| 久久99精品国语久久久| 日韩精品有码人妻一区| 国产欧美另类精品又又久久亚洲欧美| 久久久精品94久久精品| 黄色毛片三级朝国网站| 国产精品嫩草影院av在线观看| 国产精品久久久久久精品古装| 新久久久久国产一级毛片| 成人午夜精彩视频在线观看| 人人妻人人澡人人爽人人夜夜| 最新的欧美精品一区二区| 久久精品国产自在天天线| 日韩中字成人| 超色免费av| 国产成人91sexporn| 亚洲欧洲精品一区二区精品久久久 | 五月开心婷婷网| 国产成人精品无人区| 一区在线观看完整版| 免费久久久久久久精品成人欧美视频 | 久久热精品热| 黑人高潮一二区| 精品一区在线观看国产| 免费黄色在线免费观看| 国产成人freesex在线| 如日韩欧美国产精品一区二区三区 | 色婷婷久久久亚洲欧美| 超碰97精品在线观看| 麻豆成人av视频| 亚洲成色77777| 欧美人与性动交α欧美精品济南到 | 久久久久久久久久久免费av| 中文精品一卡2卡3卡4更新| 简卡轻食公司| 欧美+日韩+精品| 69精品国产乱码久久久| 美女cb高潮喷水在线观看| 久久久a久久爽久久v久久| 免费观看av网站的网址| 久久ye,这里只有精品| 一边摸一边做爽爽视频免费| 岛国毛片在线播放| av又黄又爽大尺度在线免费看| a 毛片基地| 赤兔流量卡办理| 亚洲精品日韩在线中文字幕| 亚洲色图综合在线观看| 免费黄频网站在线观看国产| 91午夜精品亚洲一区二区三区| 精品一区在线观看国产| 内地一区二区视频在线| 91久久精品电影网| 天堂8中文在线网| 亚洲人成网站在线播| 久久久久久久久久久免费av| 亚洲综合色网址| av女优亚洲男人天堂| 热99久久久久精品小说推荐| 搡女人真爽免费视频火全软件| 能在线免费看毛片的网站| tube8黄色片| 成人毛片60女人毛片免费| 男女边吃奶边做爰视频| 午夜精品国产一区二区电影| 永久免费av网站大全| 多毛熟女@视频| 久久精品国产亚洲网站| 日韩熟女老妇一区二区性免费视频| 99热网站在线观看| 日韩一区二区三区影片| 人人妻人人澡人人爽人人夜夜| 18禁裸乳无遮挡动漫免费视频| 亚洲国产精品成人久久小说| 日本黄色日本黄色录像| 97超视频在线观看视频| 亚洲五月色婷婷综合| 亚洲在久久综合| 国产国语露脸激情在线看| 国产精品久久久久久av不卡| 考比视频在线观看| 美女主播在线视频| 中文乱码字字幕精品一区二区三区| 国产不卡av网站在线观看| 观看av在线不卡| 丝袜美足系列| 热99久久久久精品小说推荐| 色婷婷av一区二区三区视频| 国产精品国产三级国产专区5o| 黄色配什么色好看| av免费在线看不卡| 在线免费观看不下载黄p国产| 伊人久久精品亚洲午夜| 97超视频在线观看视频| 国产一级毛片在线| 日本黄大片高清| 久久ye,这里只有精品| 亚洲av不卡在线观看| 欧美最新免费一区二区三区| 精品午夜福利在线看| 最近2019中文字幕mv第一页| .国产精品久久| 久久av网站| 国产精品 国内视频| 免费大片18禁| av黄色大香蕉| 久久久午夜欧美精品| 亚洲欧美日韩卡通动漫| 男人添女人高潮全过程视频| 久久精品国产自在天天线| 日本欧美国产在线视频| 九草在线视频观看| 国产精品免费大片| 超碰97精品在线观看| 黑丝袜美女国产一区| 国产亚洲欧美精品永久| 18禁动态无遮挡网站| 午夜福利在线观看免费完整高清在| 欧美国产精品一级二级三级| 久久99热6这里只有精品| 女人久久www免费人成看片| 欧美精品国产亚洲| 内地一区二区视频在线| 最新的欧美精品一区二区| 91精品伊人久久大香线蕉| 天天影视国产精品| 婷婷色综合大香蕉| 最近2019中文字幕mv第一页| 少妇高潮的动态图| 久久久国产一区二区| 婷婷成人精品国产| 国产永久视频网站| 一区二区三区乱码不卡18| 99久久精品一区二区三区| 国产精品成人在线| 老司机影院成人| 午夜福利,免费看| 日本wwww免费看| 水蜜桃什么品种好| 久久国产精品大桥未久av| 男女国产视频网站| 又黄又爽又刺激的免费视频.| 国产熟女欧美一区二区| 国产极品粉嫩免费观看在线 | 欧美日本中文国产一区发布| 18禁裸乳无遮挡动漫免费视频| 99热6这里只有精品| 亚洲经典国产精华液单| 久久久久人妻精品一区果冻| 久久国产精品男人的天堂亚洲 | 精品少妇久久久久久888优播| 国产一区有黄有色的免费视频| 超色免费av| av电影中文网址| 午夜免费男女啪啪视频观看| 我的女老师完整版在线观看| 999精品在线视频| 肉色欧美久久久久久久蜜桃| 综合色丁香网| 亚洲av成人精品一区久久| 亚洲精品,欧美精品| 一本大道久久a久久精品| 久久av网站| 亚州av有码| 韩国av在线不卡| 久久女婷五月综合色啪小说| 色视频在线一区二区三区| 日日啪夜夜爽| 制服诱惑二区| 国产精品99久久99久久久不卡 | 成年人午夜在线观看视频| 久久鲁丝午夜福利片| 在线观看三级黄色| 免费黄网站久久成人精品| 日韩av在线免费看完整版不卡| 午夜免费鲁丝| 99久久精品一区二区三区| 精品国产一区二区久久| 国产一级毛片在线| 最新的欧美精品一区二区| 日韩免费高清中文字幕av| 欧美+日韩+精品| 国产男人的电影天堂91| 亚洲精品日本国产第一区| 人妻制服诱惑在线中文字幕| 国产精品免费大片| 免费黄频网站在线观看国产| 少妇人妻 视频| 日韩 亚洲 欧美在线| 99精国产麻豆久久婷婷| 亚洲情色 制服丝袜| 又黄又爽又刺激的免费视频.| 亚洲精品自拍成人| 国产视频首页在线观看| 香蕉精品网在线| 日韩欧美精品免费久久| 国产高清国产精品国产三级| 观看美女的网站| 亚洲成人一二三区av| 国产色婷婷99| 久久久亚洲精品成人影院| 搡老乐熟女国产| 欧美亚洲日本最大视频资源| 国产又色又爽无遮挡免| 午夜视频国产福利| 国产av码专区亚洲av| 国产国拍精品亚洲av在线观看| 国产69精品久久久久777片| 精品人妻在线不人妻| 亚洲av中文av极速乱| 美女国产视频在线观看| 国产69精品久久久久777片| 两个人的视频大全免费| 午夜91福利影院| 99re6热这里在线精品视频| 精品久久久久久久久av| 又黄又爽又刺激的免费视频.| 久久精品国产自在天天线| 99热6这里只有精品| 制服丝袜香蕉在线| 日韩欧美精品免费久久| 国产白丝娇喘喷水9色精品| 18禁在线播放成人免费| 少妇被粗大的猛进出69影院 | 日本色播在线视频| 51国产日韩欧美| 欧美三级亚洲精品| 一级毛片我不卡| 国产欧美亚洲国产| 日本爱情动作片www.在线观看| 久久精品熟女亚洲av麻豆精品| 日本vs欧美在线观看视频| 999精品在线视频| 亚洲av国产av综合av卡| 亚洲欧美中文字幕日韩二区| 国产免费福利视频在线观看| 一本大道久久a久久精品| 女人久久www免费人成看片| 人人澡人人妻人| 亚洲人成网站在线观看播放| 国产精品欧美亚洲77777| 少妇精品久久久久久久| 美女国产视频在线观看| 国产在线一区二区三区精| 美女主播在线视频| 91精品国产九色| 久久99精品国语久久久| 麻豆乱淫一区二区| 91午夜精品亚洲一区二区三区| 日本欧美国产在线视频| 日本-黄色视频高清免费观看| 亚洲高清免费不卡视频| videosex国产| 高清不卡的av网站| 精品人妻熟女毛片av久久网站| 伦理电影免费视频| 伦精品一区二区三区| 精品久久蜜臀av无| 亚洲av不卡在线观看| 91国产中文字幕| 水蜜桃什么品种好| 日韩人妻高清精品专区| 国产精品99久久99久久久不卡 | 婷婷成人精品国产| 欧美日韩在线观看h| 亚洲人与动物交配视频| 91午夜精品亚洲一区二区三区| 成年人午夜在线观看视频| 夜夜看夜夜爽夜夜摸| 国产成人a∨麻豆精品| h视频一区二区三区| 国产精品人妻久久久影院| 国产亚洲精品第一综合不卡 | videossex国产| 日韩熟女老妇一区二区性免费视频| 亚洲精品久久午夜乱码| 日韩,欧美,国产一区二区三区| 国产片内射在线| 成人午夜精彩视频在线观看| 欧美少妇被猛烈插入视频| 精品久久久噜噜| 各种免费的搞黄视频| 极品少妇高潮喷水抽搐| 日本wwww免费看| 自线自在国产av| 观看av在线不卡| 嘟嘟电影网在线观看| 亚洲精品久久久久久婷婷小说| 精品人妻熟女av久视频| 久久久亚洲精品成人影院| 女性被躁到高潮视频| 亚洲欧美一区二区三区黑人 | 亚洲婷婷狠狠爱综合网| 看十八女毛片水多多多| 99久久综合免费| 免费大片黄手机在线观看| 最近中文字幕高清免费大全6| 99re6热这里在线精品视频| 9色porny在线观看| 国产精品不卡视频一区二区| 激情五月婷婷亚洲| 午夜视频国产福利| 中国国产av一级| 秋霞伦理黄片| 免费av中文字幕在线| 一级毛片aaaaaa免费看小| 日韩中文字幕视频在线看片| 中文字幕久久专区| 2018国产大陆天天弄谢| 人妻 亚洲 视频| 国产精品免费大片| av专区在线播放| 男女边吃奶边做爰视频| kizo精华| 成人影院久久| 91精品三级在线观看| 在线观看国产h片| 欧美老熟妇乱子伦牲交| 在线观看美女被高潮喷水网站| 一本一本久久a久久精品综合妖精 国产伦在线观看视频一区 | 国产成人免费无遮挡视频| 不卡视频在线观看欧美| 一区二区三区四区激情视频| 国产探花极品一区二区| 免费观看无遮挡的男女| 久久国产精品大桥未久av| 午夜91福利影院| 国产有黄有色有爽视频| 2021少妇久久久久久久久久久| 国产av精品麻豆| 一区二区三区免费毛片| 国产av码专区亚洲av| 中文字幕精品免费在线观看视频 | 搡女人真爽免费视频火全软件| 天堂俺去俺来也www色官网| 全区人妻精品视频| 午夜福利影视在线免费观看| 老司机影院毛片| 一区二区三区免费毛片| 肉色欧美久久久久久久蜜桃| 女性被躁到高潮视频| 母亲3免费完整高清在线观看 | 久久久久人妻精品一区果冻| 99国产精品免费福利视频| 欧美日韩视频高清一区二区三区二| 十八禁网站网址无遮挡| 精品少妇内射三级| 99视频精品全部免费 在线| 最后的刺客免费高清国语| 国产免费现黄频在线看| 欧美激情国产日韩精品一区| 国产一区亚洲一区在线观看| 久久这里有精品视频免费| 国产亚洲最大av| 免费看不卡的av| 中文字幕最新亚洲高清| 人人妻人人爽人人添夜夜欢视频| 99久久中文字幕三级久久日本| 欧美丝袜亚洲另类| 日韩精品有码人妻一区| 午夜福利视频在线观看免费| 桃花免费在线播放| 亚洲精品一区蜜桃| 日韩制服骚丝袜av| 妹子高潮喷水视频| 日韩伦理黄色片| 啦啦啦啦在线视频资源| av网站免费在线观看视频| 亚洲国产色片| 午夜日本视频在线| 香蕉精品网在线| 日本wwww免费看| 国产色爽女视频免费观看| 亚洲精品色激情综合| 中文字幕久久专区| 男女边吃奶边做爰视频| 在线播放无遮挡| 在线观看美女被高潮喷水网站| 亚洲,欧美,日韩| 国产综合精华液| 成人免费观看视频高清| av电影中文网址| 天天操日日干夜夜撸| 丝袜脚勾引网站| av网站免费在线观看视频| 97精品久久久久久久久久精品| 免费高清在线观看视频在线观看| 日本黄色片子视频| 日本wwww免费看| 精品久久久噜噜| 亚洲精品国产色婷婷电影| 在线播放无遮挡| 国精品久久久久久国模美| 一个人看视频在线观看www免费| 国精品久久久久久国模美| 最近中文字幕2019免费版| av不卡在线播放| 国产永久视频网站| 美女视频免费永久观看网站| 日本vs欧美在线观看视频| 大又大粗又爽又黄少妇毛片口| 成人影院久久| 女性生殖器流出的白浆| 国产精品熟女久久久久浪| a级毛片在线看网站| 建设人人有责人人尽责人人享有的| 国产熟女午夜一区二区三区 | 国产伦精品一区二区三区视频9| 天天影视国产精品| 91久久精品国产一区二区三区| 不卡视频在线观看欧美| 国产高清有码在线观看视频| 中文字幕人妻丝袜制服| 99热6这里只有精品| 国产一区二区在线观看日韩| 少妇精品久久久久久久| 搡老乐熟女国产| 欧美日本中文国产一区发布| av电影中文网址| 男女免费视频国产| 熟女人妻精品中文字幕| 男人操女人黄网站| 久久精品国产鲁丝片午夜精品| 2018国产大陆天天弄谢| 99国产精品免费福利视频| av国产久精品久网站免费入址| 午夜福利在线观看免费完整高清在| 99热6这里只有精品| 最新中文字幕久久久久| tube8黄色片| 一级毛片电影观看| 国产精品一区二区三区四区免费观看| 大香蕉97超碰在线| 国产精品偷伦视频观看了| 99久久人妻综合| 国产成人精品婷婷| 亚洲在久久综合| 亚洲综合精品二区| 日韩av不卡免费在线播放| 免费av不卡在线播放| 简卡轻食公司| 欧美人与性动交α欧美精品济南到 | 久久精品国产自在天天线| 久久精品国产鲁丝片午夜精品| 午夜久久久在线观看| 欧美人与善性xxx| 一本久久精品| 日产精品乱码卡一卡2卡三| 午夜精品国产一区二区电影| 中文精品一卡2卡3卡4更新| 人人妻人人爽人人添夜夜欢视频| 能在线免费看毛片的网站| 日韩av在线免费看完整版不卡| 亚洲国产成人一精品久久久| 一本色道久久久久久精品综合| 久久久精品免费免费高清| 久久97久久精品| 国产av精品麻豆| 亚洲五月色婷婷综合| 亚洲美女视频黄频| 两个人免费观看高清视频| 91久久精品国产一区二区成人| 纵有疾风起免费观看全集完整版| 亚洲国产欧美日韩在线播放| 爱豆传媒免费全集在线观看| 亚洲精华国产精华液的使用体验| 久久韩国三级中文字幕| 22中文网久久字幕| 纯流量卡能插随身wifi吗| 国产精品人妻久久久影院| 欧美+日韩+精品| 国产精品久久久久成人av| 国产午夜精品久久久久久一区二区三区| 十八禁高潮呻吟视频| 欧美另类一区| 国产精品.久久久| av女优亚洲男人天堂| 3wmmmm亚洲av在线观看| 日韩电影二区| 国产高清有码在线观看视频| 两个人的视频大全免费| 久久久久久久大尺度免费视频| 又粗又硬又长又爽又黄的视频| 伊人久久精品亚洲午夜| 女性被躁到高潮视频| 欧美日韩精品成人综合77777| 亚洲av成人精品一区久久| 色视频在线一区二区三区| 制服丝袜香蕉在线| 丝袜在线中文字幕| 五月玫瑰六月丁香| 欧美+日韩+精品| 免费日韩欧美在线观看| 欧美人与性动交α欧美精品济南到 | 国产精品女同一区二区软件| 人体艺术视频欧美日本| 亚洲五月色婷婷综合| 国产成人免费观看mmmm| 国产毛片在线视频| 大香蕉97超碰在线| 美女主播在线视频| 日韩欧美一区视频在线观看| 人人妻人人添人人爽欧美一区卜|