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

    基于RNA-seq數(shù)據(jù)發(fā)掘響應(yīng)干旱的谷子C2H2型鋅指蛋白基因

    2019-08-13 08:54:10齊帆范興陳義
    江蘇農(nóng)業(yè)科學 2019年2期
    關(guān)鍵詞:差異表達基因干旱脅迫谷子

    齊帆 范興 陳義

    摘要:C2H2型鋅指蛋白在植物響應(yīng)逆境脅迫以及激素信號轉(zhuǎn)導過程中起著重要作用,為了發(fā)掘響應(yīng)干旱脅迫的谷子C2H2型鋅指蛋白基因(SiC2H2),利用轉(zhuǎn)錄組測序技術(shù)構(gòu)建谷子聚乙二醇(PEG)處理前后轉(zhuǎn)錄組文庫,根據(jù)表達譜數(shù)據(jù)鑒定出11個SiC2H2差異表達基因,對其進行理化性質(zhì)、啟動子元件等部分生物信息學分析和進化樹及基因結(jié)構(gòu)分析。實時熒光定量逆轉(zhuǎn)錄PCR(RT-PCR)驗證結(jié)果表明,隨機選取的4個基因在PEG脅迫前后的表達特點與表達譜數(shù)據(jù)基本一致。研究結(jié)果為后續(xù)進一步研究C2H2型鋅指蛋白在谷子抗旱性中的作用以及谷子抗旱分子育種提供了理論依據(jù)。

    關(guān)鍵詞:谷子;C2H2型鋅指蛋白;干旱脅迫;差異表達基因

    中圖分類號: S515.01 ?文獻標志碼: A ?文章編號:1002-1302(2019)02-0041-05

    干旱是我國糧食生產(chǎn)的重要制約因素之一,且干旱的頻率和程度日趨嚴重,因此選育抗旱品種已成為育種的重要目標。谷子[Setaria italica (L.) P. Beauv.]屬禾本科狗尾草屬,具有耐旱、耐瘠薄、適應(yīng)性強、水分利用率高等特點,在作物抗旱種質(zhì)資源利用方面具有其獨特的研究價值[1]。谷子全基因組測序的完成為研究谷子功能基因組提供了平臺[2-3],發(fā)掘與鑒定谷子抗旱基因資源、探索谷子抗旱機制對于促進其基因資源的利用和抗旱育種有著極為重要的理論和實踐意義。

    鋅指蛋白是植物轉(zhuǎn)錄調(diào)節(jié)子中最大的家族之一,具有指狀結(jié)構(gòu)且能結(jié)合鋅離子[4],其中的C2H2型鋅指蛋白具有一段高度保守的氨基酸特征序列QALGGH,含有由大約30個氨基酸組成的鋅指結(jié)構(gòu),該結(jié)構(gòu)的保守序列為X2CX2~4CX12HX2~8H[5]。C2H2型鋅指蛋白廣泛參與調(diào)控植物的生長和發(fā)育[6-9],而且越來越多的研究表明,C2H2型鋅指蛋白在植物響應(yīng)逆境脅迫過程中起著非常重要的作用,許多與植物脅迫耐受相關(guān)的C2H2型鋅指蛋白已得到鑒定[10-12]。DNA芯片分析結(jié)果顯示,擬南芥C2H2型鋅指蛋白中Cl-2i類(Zat家族)的一些成員在各種生物和非生物脅迫下上調(diào)表達,其中2個被廣泛研究的成員Zat10(STZ)和Zat12,能夠被干旱、鹽、滲透、溫度、傷害、氧化和高光等脅迫所誘導;遺傳分析結(jié)果也表明,Zat10和Zat12參與調(diào)節(jié)了活性氧清除基因的表達進而提高了植物對干旱、鹽漬以及氧化脅迫的忍耐力[6]。水稻中C2H2型鋅指蛋白的部分成員,如ZFP245、ZFP252和ZFP179也被證明參與了水稻對干旱、鹽和氧化脅迫的應(yīng)激反應(yīng)[13]。水稻ZFP182和ZFP36參與了脫落酸(ABA)誘導的抗氧化防護,在提高植物對水分脅迫以及氧化脅迫的忍耐力方面起到了重要的調(diào)節(jié)作用[14-15]。在煙草中過量表達矮牽牛C2H2型鋅指蛋白基因ZPT2-3,明顯提高了轉(zhuǎn)基因煙草對干旱脅迫的耐受能力[16]。目前已經(jīng)在許多植物,如擬南芥、水稻、小麥、煙草中分離到近百個C2H2型鋅指蛋白基因。

    谷子作為我國重要的農(nóng)作物之一,近年來已逐漸成為開展C4植物和抗旱耐逆研究的模式植物[17]。目前人們已經(jīng)從谷子中鑒定出124個SiC2H2,它們大都參與非生物逆境脅迫響應(yīng)[18],然而具體哪些成員參與了干旱脅迫響應(yīng)過程還不清楚,在這一過程中的具體作用和功能還有待于進一步研究。本研究擬通過轉(zhuǎn)錄組測序表達譜數(shù)據(jù),從谷子中鑒定響應(yīng)干旱脅迫的C2H2型鋅指蛋白,篩選參與干旱脅迫應(yīng)答的關(guān)鍵SiC2H2基因,以期為闡明C2H2型鋅指蛋白在谷子抗旱性中的作用以及利用分子生物學手段提高作物的抗旱性奠定基礎(chǔ)。

    1 材料與方法

    試驗于2016年9月至2017年9月在山西農(nóng)業(yè)大學農(nóng)業(yè)生物工程研究所、雜糧基因與分子育種山西省高等學校重點實驗室內(nèi)完成。

    1.1 試驗材料

    以谷子JF16為試驗材料,其種子由山西農(nóng)業(yè)大學農(nóng)業(yè)生物工程研究所提供。將谷子在人工氣候室中培養(yǎng),采用體積比為3 ∶ 1的營養(yǎng)土、蛭石混合土種植,培養(yǎng)3周后,對谷子幼苗進行20% PEG-6000(聚乙二醇-6000)溶液模擬干旱脅迫處理,以蒸餾水處理作為對照,處理0.5 h后分別剪取幼苗植株葉片,迅速放于液氮中冷凍,用于轉(zhuǎn)錄組測序;同上,處理0、0.5、1、2、4 h后取材,用于實時熒光定量PCR驗證。

    1.2 SiC2H2差異表達基因的篩選

    基因表達量用TPM值(transcripts per million,即每百萬條reads的轉(zhuǎn)錄本)表示,運用R語言繪制基因的表達值熱圖。根據(jù)對照組和干旱處理組中基因的TPM值來比較不同處理間的基因差異表達。本研究中差異基因篩選標準為錯誤發(fā)現(xiàn)率(false discovery rate,簡稱FDR)≤0.001且差異倍數(shù)(fold change,簡稱FC)在2倍以上,即∣log2(fold change)∣≥1。

    1.3 SiC2H2基因生物信息學分析

    從谷子基因組數(shù)據(jù)庫phytozome v 12.0(https://phytozome.jgi.doe.gov/pz/portal.html)中找出所有C2H2型鋅指蛋白基因,并從轉(zhuǎn)錄組數(shù)據(jù)中查找到相對應(yīng)的基因,提取其在基因組數(shù)據(jù)庫中的基本信息;利用DNAMAN6.0預(yù)測蛋白質(zhì)分子量和等電點等理化性質(zhì);用Plant CARE(http:bioinformatics.psb.ugent.be/webtools/plantcare/html/)在線分析SiC2H2基因啟動子(起始密碼子ATG上游1.5 kb的DNA序列)順式作用元件;使用MEGA6.0構(gòu)建SiC2H2氨基酸序列系統(tǒng)進化樹;利用GSDS在線工具(http://gsds.cbi.pku.edu.cn/)制作SiC2H2基因結(jié)構(gòu)圖。

    1.4 實時熒光定量RT-PCR分析

    參照RNAiso Plus、PrimeScriptTM RT reagent Kit(TaKaRa)使用說明,從各處理樣品中提取總RNA以及合成cDNA。根據(jù)基因序列信息,采用軟件Primer 5.0進行real-time RT-PCR特異性引物設(shè)計(表1)。使用美國Bio-Rad CFX96TM Touch Real-Time PCR System對基因進行相對定量RT-PCR分析。兩步法PCR反應(yīng)程序:95 ℃預(yù)變性30 s;95 ℃ 5 s,60 ℃ 30 s,40個循環(huán)。每個樣品重復(fù)3次,以β-actin作為內(nèi)參,用2-ΔΔCT的方法計算基因的相對表達量,基因在處理 0 h 時的相對表達量設(shè)定為1。

    2 結(jié)果與分析

    2.1 SiC2H2差異表達基因篩選

    利用高通量轉(zhuǎn)錄組測序技術(shù)從表達譜數(shù)據(jù)中共得到44個SiC2H2基因,將這些基因在正常澆水對照組(control)與PEG模擬干旱處理組(PEG)中的表達量進行了標準化熱圖分析(圖1),其中18個基因上調(diào)表達,25個下調(diào)表達,1個無差異表達。

    為篩選鑒定響應(yīng)干旱的關(guān)鍵SiC2H2基因,對正常澆水和干旱處理的表達譜進行比較,根據(jù)基因表達量(TPM值)計算該基因在同一品種不同處理間的差異表達倍數(shù),以FDR≤0.001、∣log2(FC)∣≥1(即倍數(shù)差異在2倍以上)為篩選標準,共獲得11個差異表達基因(圖2),其中編號為Si017660m、Si022296m、Si023013m、Si003051m、Si037343m、Si037672m的基因在干旱脅迫下與對照相比上調(diào)表達,而編號為Si016736m、Si018078m、Si029540m、Si037164m、Si029897m的基因下調(diào)表達。

    2.2 SiC2H2基本信息與理化性質(zhì)

    針對獲得的11個差異表達基因,對其基本信息與理化性質(zhì)進行了分析,如表2所示,它們分別位于第1、2、3、5、9號染色體上,基因序列長度處于844~7 484 bp之間,平均長度為2 652 bp;氨基酸序列長度處于184~585個之間,平均長度為332個,分子量處于19.4~60.2 ku之間,平均分子量為 352 ku;理論等電點為6.05~ 9.17,平均為7.12。

    2.3 SiC2H2進化樹及基因結(jié)構(gòu)分析

    谷子C2H2型鋅指蛋白家族成員分成5個亞家族,本研究使用MEGA6.0軟件中的Neighbor-Joining算法對11個SiC2H2氨基酸序列構(gòu)建系統(tǒng)進化樹,發(fā)現(xiàn)它們分別屬于Ⅱ、Ⅲ、Ⅳ、Ⅴ亞家族,基因結(jié)構(gòu)分析顯示,同一亞族的基因具有相似的外顯子-內(nèi)含子結(jié)構(gòu)(圖3)。

    2.4 SiC2H2基因啟動子元件分析

    啟動子是基因表達與調(diào)控的“分子開關(guān)”,對啟動子順式作用元件的分析將有助于研究該基因的表達調(diào)控機制及其生物學功能。如表3所示,11個SiC2H2基因啟動子區(qū)包含多種激素響應(yīng)順式作用元件,如核心序列為CACGTG的脫落酸應(yīng)答元件ABRE、水楊酸應(yīng)答元件TCA-element、乙烯應(yīng)答元件ERE、赤霉素應(yīng)答元件GARE-motif、茉莉酸甲酯應(yīng)答元件CGTCA-motif及生長素應(yīng)答元件TGA-element;同時包含多種參與植物非生物逆境脅迫響應(yīng)的順式作用元件,如干旱誘導相關(guān)的MYB結(jié)合位點MBS、熱激應(yīng)答元件HSE、參與低溫響應(yīng)的順式作用元件LTR以及參與防御和脅迫響應(yīng)的元件 TC-rich repeats。每個基因至少含有6個與激素及脅迫應(yīng)答相關(guān)的順式作用元件,并且所有基因都含有與干旱脅迫響應(yīng)相關(guān)的ABRE作用元件或者MBS作用元件,表明這些SiC2H2基因極有可能參與調(diào)控谷子生長發(fā)育及脅迫應(yīng)激反應(yīng),尤其是干旱脅迫響應(yīng)。

    2.5 SiC2H2基因熒光定量RT-PCR驗證

    為了進一步驗證表達譜中響應(yīng)干旱的SiC2H2基因表達特性,以PEG處理0、0.5、1、2、4 h的谷子幼苗為材料,隨機篩選4個基因Si037343m、Si023013m、Si022296m、Si037164m,使用實時熒光定量PCR技術(shù)對它們的基因表達情況進行驗證。如圖4所示,4個基因在PEG干旱脅迫處理條件下的表達趨勢與表達譜數(shù)據(jù)結(jié)果基本吻合,說明表達譜數(shù)據(jù)可靠,可用于后續(xù)研究。

    3 討論與結(jié)論

    干旱是世界范圍內(nèi)嚴重的環(huán)境脅迫因子,為保證全球人口日益增加的糧食供應(yīng),在干旱條件下最大限度地提高作物產(chǎn)量潛力、保持作物產(chǎn)量穩(wěn)定是十分重要的。植物感知和應(yīng)對干旱脅迫的機制非常復(fù)雜,在耐旱脅迫反應(yīng)過程中許多干旱響應(yīng)基因被激活,其中植物C2H2型鋅指蛋白作為植物轉(zhuǎn)錄因子大家族之一,廣泛參與了植物生長發(fā)育與逆境脅迫信號轉(zhuǎn)導,許多C2H2型鋅指蛋白介導的干旱脅迫響應(yīng)過程被發(fā)現(xiàn)。最近的研究表明,擬南芥ZAT18在干旱脅迫下相應(yīng)基因表達量顯著上調(diào),過表達植株與對照相比具有更強的耐旱性,ZAT18在干旱脅迫響應(yīng)過程中起到正調(diào)控作用[19],而大豆中的1個具有QALGGH保守結(jié)構(gòu)域的C2H2型鋅指蛋白GmZFP3則通過ABA依賴途徑負調(diào)控植物對干旱脅迫的響應(yīng)[20]。甘薯IbZFP1通過參與調(diào)控ABA信號轉(zhuǎn)導、脯氨酸合成及活性氧清除進而增強了轉(zhuǎn)基因擬南芥對鹽脅迫與干旱脅迫的耐性[21]。

    隨著全基因組測序的完成,目前在谷子中發(fā)現(xiàn)124個C2H2型鋅指蛋白家族成員[18],本研究通過分析谷子PEG處理葉片轉(zhuǎn)錄組,進一步研究了其中11個響應(yīng)干旱的差異表達基因,其中6個在干旱脅迫下上調(diào)表達,5個下調(diào)表達,并對這11個應(yīng)答干旱脅迫的SiC2H2基因進行了理化性質(zhì)等生物信息學分析。這些SiC2H2基因啟動子序列區(qū)中含有多種與逆境以及激素相關(guān)的順式作用元件,其中大多數(shù)基因都具有脫落酸應(yīng)答元件ABRE、干旱誘導相關(guān)的MYB結(jié)合位點MBS以及茉莉酸甲酯應(yīng)答元件CGTCA-motif,說明它們極有可能通過依賴ABA途徑或者非ABA依賴途徑信號轉(zhuǎn)導過程來調(diào)控植物對干旱脅迫的響應(yīng)。使用熒光定量RT-PCR對4個隨機選擇的基因進行表達情況檢測,它們在PEG脅迫下的表達特性總體趨勢與表達譜測序結(jié)果基本是吻合的,但在具體表達量上調(diào)或下調(diào)倍數(shù)以及出現(xiàn)峰值的時間點上存在差異,這可能與所取材料、數(shù)據(jù)分析方法及檢測靈敏度等不同有關(guān)。

    本研究基于PEG干旱脅迫谷子葉片轉(zhuǎn)錄組測序數(shù)據(jù),鑒定出44個SiC2H2基因,并對其中11個差異表達基因進行生物信息學分析及實時熒光定量RT-PCR驗證,研究結(jié)果為從轉(zhuǎn)錄組水平揭示谷子響應(yīng)干旱脅迫的分子調(diào)控機制提供了理論依據(jù)。

    參考文獻:

    [1]Diao X M,Schnable J,Bennetzen J L,et al. Initiation of Setaria as a model plant[J]. Frontiers of Agricultural Science and Engineering,2014,1(1):16-20.

    [2]Zhang G,Liu X,Quan Z,et al. Genome sequence of foxtail millet (Setaria italica) provides insights into grass evolution and biofuel potential[J]. Nature Biotechnology,2012,30(6):549-554.

    [3]Bennetzen J L,Schmutz J,Wang H,et al. Reference genome sequence of the model plant Setaria[J]. Nature Biotechnology,2012,30(6):555-561.

    [4]Takatsuji H. Zinc-finger proteins:the classical zinc finger emerges in contemporary plant science[J]. Plant Molecular Biology,1999,39(6):1073-1078.

    [5]Pabo C O,Peisach E,Grant R A. Design and selection of novel Cys2His2 zinc finger proteins[J]. Annual Review of Biochemistry,2001,70(70):313-340.

    [6]Kiebowicz-Matuk A. Involvement of plant C2H2-type zinc finger transcription factors in stress responses[J]. Plant Science,2012,185-186(4):78-85.

    [7]Sakamoto H,Maruyama K,Sakuma Y,et al. Arabidopsis Cys2/his2-type zinc-finger proteins function as transcription repressors under drought,cold,and high-salinity stress conditions[J]. Plant Physiology,2004,136(1):2734-2746.

    [8]Xiao H,Tang J,Li Y,et al. STAMENLESS 1,encoding a single C2H2 zinc finger protein,regulates floral organ identity in rice[J]. Plant Journal,2009,59(5):789-801.

    [9]Wang Y,Dou D,Wang X,et al. The PsCZF1 gene encoding a C2H2 zinc finger protein is required for growth,development and pathogenesis in Phytophthorasojae[J]. Microbial Pathogenesis,2009,47(2):78-86.

    [10]Tian Z D,Zhang Y,Liu J,et al. Novel potato C2H2-type zinc finger protein gene,StZFP1,which responds to biotic and abiotic stress,plays a role in salt tolerance[J]. Plant Biology,2010,12(5):689-697.

    [11]Sun S J,Guo S Q,Yang X,et al. Functional analysis of a novel Cys2/His2-type zinc finger protein involved in salt tolerance in rice[J]. Journal of Experimental Botany,2010,61(10):2807-2818.

    [12]Kodaira K S,Qin F,Lam-Son P T,et al. Arabidopsis Cys2/His2 zinc-finger proteins AZF1 and AZF2 negatively regulate abscisic acid-repressive and auxin-inducible genes under abiotic stress conditions[J]. Plant Physiology,2011,157(2):742-756.

    [13]Huang J,Sun S,Xu D,et al. A TFIIIA-type zinc finger protein confers multiple abiotic stress tolerances in transgenic rice (Oryza sativa L.)[J]. Plant Molecular Biology,2012,80(3):337-350.

    [14]Zhang H,Ni L,Liu Y P,et al. The C2H2-type Zinc Finger Protein ZFP182 is involved in abscisic acid-Induced antioxidant defense in rice[J]. Journal of Integrative Plant Biology,2012,54(7):500-510.

    [15]Zhang H,Liu Y,Wen F,et al. A novel rice C2H2-type zinc finger protein,ZFP36,is a key player involved in abscisic acid-induced antioxidant defence and oxidative stress tolerance in rice[J]. Journal of Experimental Botany,2014,65(20):5795-5809.

    [16]Sugano S,Kaminaka H,Rybka Z,et al. Stress-responsive zinc finger gene ZPT2-3 plays a role in drought tolerance in petunia[J]. Plant Journal,2003,36(6):830-841.

    [17]賈冠清,刁現(xiàn)民. 谷子[Setaria italica (L.) P. Beauv.]作為功能基因組研究模式植物的發(fā)展現(xiàn)狀及趨勢[J]. 生命科學,2017(3):292-301.

    [18]Muthamilarasan M,Bonthala V S,Mishra A K,et al. C2H2 type of zinc finger transcription factors in foxtail millet define response to abiotic stresses[J]. Functional and Integrative Genomics,2014,14(3):531-543.

    [19]Yin M,Wang Y,Zhang L,et al. The arabidopsis Cys2/his2 zinc finger transcription factor ZAT18 is a positive regulator of plant tolerance to drought stress[J]. Journal of Experimental Botany,2017,68(11):2991-3005.

    [20]Zhang D,Tong J,Xu Z,et al. Soybean C2H2-type zinc finger protein GmZFP3 with conserved QALGGH motif negatively regulates drought responses in transgenic Arabidopsis[J]. Frontiers in Plant Science,2016,7(1237):325.

    [21]Wang F,Tong W,Zhu H,et al. A novel Cys2/His2 zinc finger protein gene from sweetpotato,IbZFP1,is involved in salt and drought tolerance in transgenic Arabidopsis[J]. Planta,2016,243(3):783-797.

    猜你喜歡
    差異表達基因干旱脅迫谷子
    打谷子
    當代陜西(2021年21期)2022-01-19 02:00:32
    生物信息學分析患有乳腺癌的乳腺球樣本中與自我更新相關(guān)的關(guān)鍵基因
    一氧化氮參與水楊酸對玉米幼苗根系抗旱性的調(diào)控
    干旱脅迫對扁豆生長與生理特性的影響
    不同水分條件下硫肥對玉米幼苗葉片光合特性的影響
    干旱脅迫對金花茶幼苗光合生理特性的影響
    條斑紫菜優(yōu)良品系的基因芯片表達譜分析
    谷子栽培技術(shù)
    高溫脅迫下草坪草高羊茅差異表達基因的分子研究
    HOW TO CORRECT GRAMMATICAL MISTAKES OF CHINESE ESL STUDENTS IN THEIR COMPOSITONS
    黑人巨大精品欧美一区二区蜜桃| 9热在线视频观看99| 久久久欧美国产精品| 国产精品久久久久成人av| 两个人看的免费小视频| 久久99一区二区三区| 亚洲av美国av| 亚洲精品国产av蜜桃| 99热全是精品| 国产精品秋霞免费鲁丝片| 日本五十路高清| 我的亚洲天堂| 99香蕉大伊视频| 亚洲国产欧美网| 欧美黄色淫秽网站| 久久女婷五月综合色啪小说| 日本91视频免费播放| 亚洲中文日韩欧美视频| 国产97色在线日韩免费| 日韩视频在线欧美| 久久ye,这里只有精品| 老汉色av国产亚洲站长工具| 一级黄色大片毛片| 久久久久久久久久久久大奶| 电影成人av| 亚洲欧美色中文字幕在线| 两个人看的免费小视频| 黄色a级毛片大全视频| 成人国产av品久久久| 日韩电影二区| 精品高清国产在线一区| 婷婷丁香在线五月| 99热国产这里只有精品6| 大陆偷拍与自拍| 午夜福利免费观看在线| 精品国产乱码久久久久久男人| 久久精品人人爽人人爽视色| 国产欧美亚洲国产| 亚洲成人国产一区在线观看 | 日韩人妻精品一区2区三区| 日韩av在线免费看完整版不卡| 欧美xxⅹ黑人| 中文字幕另类日韩欧美亚洲嫩草| 亚洲中文av在线| av视频免费观看在线观看| 久久亚洲精品不卡| 丝瓜视频免费看黄片| 国产欧美日韩一区二区三区在线| 欧美日韩国产mv在线观看视频| 日本猛色少妇xxxxx猛交久久| 国产主播在线观看一区二区 | 国产免费视频播放在线视频| 国产在线免费精品| 亚洲视频免费观看视频| 极品少妇高潮喷水抽搐| 亚洲国产成人一精品久久久| 一级黄片播放器| 97人妻天天添夜夜摸| 美女国产高潮福利片在线看| 久久99精品国语久久久| 青草久久国产| 国产一区有黄有色的免费视频| 国产日韩欧美在线精品| 少妇精品久久久久久久| 国产精品国产av在线观看| 国产精品二区激情视频| 国产在线免费精品| 国产日韩欧美亚洲二区| 亚洲av综合色区一区| 国产深夜福利视频在线观看| 精品亚洲乱码少妇综合久久| 午夜日韩欧美国产| 日韩一卡2卡3卡4卡2021年| 天天操日日干夜夜撸| 国产一卡二卡三卡精品| 亚洲av综合色区一区| 精品国产国语对白av| av天堂久久9| 国产麻豆69| www.av在线官网国产| 丁香六月天网| 亚洲精品一区蜜桃| 80岁老熟妇乱子伦牲交| 免费在线观看日本一区| 中文字幕最新亚洲高清| 久久影院123| 天天躁日日躁夜夜躁夜夜| 久久精品亚洲熟妇少妇任你| 老司机亚洲免费影院| 女人爽到高潮嗷嗷叫在线视频| 女性被躁到高潮视频| 97人妻天天添夜夜摸| 午夜福利视频精品| 精品高清国产在线一区| 黑人欧美特级aaaaaa片| 久久久国产一区二区| 国产成人欧美在线观看 | 亚洲五月婷婷丁香| www.999成人在线观看| 久久国产精品男人的天堂亚洲| 欧美日韩黄片免| 久久久久久久国产电影| 亚洲中文日韩欧美视频| 在线观看免费高清a一片| 欧美精品啪啪一区二区三区 | 精品少妇久久久久久888优播| 亚洲精品久久午夜乱码| 久久综合国产亚洲精品| 黄色怎么调成土黄色| 亚洲国产中文字幕在线视频| 久久久精品94久久精品| 如日韩欧美国产精品一区二区三区| 国产精品成人在线| 亚洲天堂av无毛| 亚洲国产av新网站| 制服诱惑二区| 操出白浆在线播放| 成人手机av| 99国产综合亚洲精品| 五月开心婷婷网| 大话2 男鬼变身卡| e午夜精品久久久久久久| 69精品国产乱码久久久| 成年女人毛片免费观看观看9 | 热99久久久久精品小说推荐| av欧美777| 波多野结衣av一区二区av| 欧美精品高潮呻吟av久久| 亚洲精品一二三| 午夜免费男女啪啪视频观看| 看免费成人av毛片| 亚洲成色77777| 51午夜福利影视在线观看| 赤兔流量卡办理| 美女视频免费永久观看网站| 美女视频免费永久观看网站| 97人妻天天添夜夜摸| 热99国产精品久久久久久7| 午夜激情av网站| 美女视频免费永久观看网站| 丝袜脚勾引网站| 无限看片的www在线观看| 国产精品一区二区精品视频观看| 老熟女久久久| 女人高潮潮喷娇喘18禁视频| 乱人伦中国视频| www.自偷自拍.com| 每晚都被弄得嗷嗷叫到高潮| 色播在线永久视频| 欧美乱码精品一区二区三区| 黄色毛片三级朝国网站| 久久精品亚洲熟妇少妇任你| 不卡av一区二区三区| 亚洲国产欧美日韩在线播放| 男人添女人高潮全过程视频| 一二三四社区在线视频社区8| 精品国产乱码久久久久久男人| av网站免费在线观看视频| 这个男人来自地球电影免费观看| 成人黄色视频免费在线看| 各种免费的搞黄视频| 夫妻午夜视频| 脱女人内裤的视频| 成年女人毛片免费观看观看9 | 国产黄频视频在线观看| 国产伦人伦偷精品视频| 美女主播在线视频| 国产高清视频在线播放一区 | 七月丁香在线播放| 老鸭窝网址在线观看| 亚洲国产精品一区二区三区在线| 成人亚洲欧美一区二区av| 宅男免费午夜| 五月开心婷婷网| 欧美日韩亚洲综合一区二区三区_| 视频在线观看一区二区三区| 国产国语露脸激情在线看| 亚洲av电影在线观看一区二区三区| 欧美 亚洲 国产 日韩一| 亚洲 国产 在线| 国产不卡av网站在线观看| 天堂中文最新版在线下载| 国产熟女午夜一区二区三区| 看十八女毛片水多多多| 国产一区二区三区av在线| xxx大片免费视频| 国产av国产精品国产| 精品福利观看| 午夜免费成人在线视频| 亚洲av美国av| 一级,二级,三级黄色视频| 操美女的视频在线观看| 搡老乐熟女国产| 成年人午夜在线观看视频| a级毛片黄视频| 人人妻,人人澡人人爽秒播 | 亚洲五月婷婷丁香| 国产成人欧美在线观看 | 涩涩av久久男人的天堂| 看免费成人av毛片| 日韩制服骚丝袜av| 亚洲色图综合在线观看| 少妇粗大呻吟视频| 看免费成人av毛片| 大话2 男鬼变身卡| 国产精品免费大片| 一边摸一边抽搐一进一出视频| 国产日韩欧美亚洲二区| 免费观看av网站的网址| 成人国产av品久久久| av线在线观看网站| 国产男女超爽视频在线观看| 国产国语露脸激情在线看| 最近手机中文字幕大全| 十八禁高潮呻吟视频| 久久久精品免费免费高清| 又大又爽又粗| 日本av手机在线免费观看| 国产欧美日韩一区二区三区在线| 午夜福利免费观看在线| 久久这里只有精品19| 日韩一卡2卡3卡4卡2021年| 一本综合久久免费| 亚洲精品久久成人aⅴ小说| 在线观看免费日韩欧美大片| 日日爽夜夜爽网站| 少妇精品久久久久久久| 视频在线观看一区二区三区| 久久久精品94久久精品| 国产精品久久久久久精品古装| 久久亚洲精品不卡| 一级片免费观看大全| 欧美精品一区二区大全| 一级片'在线观看视频| 亚洲自偷自拍图片 自拍| 国产日韩欧美在线精品| 1024视频免费在线观看| 亚洲综合色网址| 手机成人av网站| 精品亚洲成国产av| 国产高清不卡午夜福利| 又粗又硬又长又爽又黄的视频| 国产精品一区二区精品视频观看| 亚洲国产看品久久| 天堂中文最新版在线下载| 亚洲综合色网址| 国产日韩欧美视频二区| 一二三四社区在线视频社区8| 欧美人与性动交α欧美精品济南到| 亚洲 欧美一区二区三区| 午夜视频精品福利| 99精国产麻豆久久婷婷| 亚洲成国产人片在线观看| av不卡在线播放| 在线观看一区二区三区激情| 天堂中文最新版在线下载| 久久久精品94久久精品| 少妇粗大呻吟视频| 天天添夜夜摸| 日韩中文字幕视频在线看片| 久久综合国产亚洲精品| 一边摸一边做爽爽视频免费| 亚洲伊人色综图| 久久久精品免费免费高清| 亚洲欧美一区二区三区黑人| 两人在一起打扑克的视频| 国产激情久久老熟女| 老司机深夜福利视频在线观看 | 麻豆国产av国片精品| 日本wwww免费看| 成人国语在线视频| 在线观看免费日韩欧美大片| 在线观看国产h片| 老汉色∧v一级毛片| 午夜福利视频在线观看免费| 女性生殖器流出的白浆| 国产成人av教育| 久久精品亚洲熟妇少妇任你| 国产欧美日韩综合在线一区二区| 午夜老司机福利片| 国产av国产精品国产| 久久人人97超碰香蕉20202| av在线老鸭窝| 日韩av在线免费看完整版不卡| 黑人猛操日本美女一级片| 男女国产视频网站| 99国产精品一区二区蜜桃av | 天天躁狠狠躁夜夜躁狠狠躁| 午夜福利免费观看在线| 亚洲熟女精品中文字幕| 欧美中文综合在线视频| 午夜91福利影院| 每晚都被弄得嗷嗷叫到高潮| 亚洲精品美女久久av网站| 婷婷色麻豆天堂久久| 亚洲国产精品一区二区三区在线| 人成视频在线观看免费观看| 首页视频小说图片口味搜索 | 久久久久久久大尺度免费视频| 尾随美女入室| 别揉我奶头~嗯~啊~动态视频 | 中文字幕av电影在线播放| 亚洲精品国产区一区二| 国产精品.久久久| 欧美亚洲 丝袜 人妻 在线| 日韩av免费高清视频| 手机成人av网站| 男人添女人高潮全过程视频| av网站免费在线观看视频| av天堂久久9| 国产精品久久久久久精品古装| 国产精品国产av在线观看| 精品福利永久在线观看| 亚洲国产av影院在线观看| 一区福利在线观看| 欧美成人精品欧美一级黄| 久久精品熟女亚洲av麻豆精品| 亚洲国产精品一区二区三区在线| 伦理电影免费视频| 亚洲精品一卡2卡三卡4卡5卡 | 亚洲欧美一区二区三区国产| 一边摸一边做爽爽视频免费| 亚洲第一青青草原| 一级毛片黄色毛片免费观看视频| 99精品久久久久人妻精品| 50天的宝宝边吃奶边哭怎么回事| 国产老妇伦熟女老妇高清| 欧美在线黄色| 国产一区二区三区av在线| 久久国产精品人妻蜜桃| 人人澡人人妻人| 天堂俺去俺来也www色官网| 好男人电影高清在线观看| 一本一本久久a久久精品综合妖精| 天堂中文最新版在线下载| 高清av免费在线| 97人妻天天添夜夜摸| 人人妻,人人澡人人爽秒播 | 看十八女毛片水多多多| 国产伦人伦偷精品视频| 又大又黄又爽视频免费| www.精华液| 亚洲精品日本国产第一区| 久久精品国产亚洲av涩爱| 一区二区av电影网| 美女国产高潮福利片在线看| 黄片小视频在线播放| 国产成人av教育| 中文字幕av电影在线播放| 欧美成人精品欧美一级黄| 尾随美女入室| 国产精品偷伦视频观看了| 免费人妻精品一区二区三区视频| 50天的宝宝边吃奶边哭怎么回事| 欧美精品亚洲一区二区| av线在线观看网站| 亚洲精品日韩在线中文字幕| 人人妻人人添人人爽欧美一区卜| 美女视频免费永久观看网站| 亚洲人成77777在线视频| 女人精品久久久久毛片| 欧美性长视频在线观看| www.自偷自拍.com| 夫妻性生交免费视频一级片| 在线av久久热| 欧美性长视频在线观看| 男女床上黄色一级片免费看| 爱豆传媒免费全集在线观看| 日韩 亚洲 欧美在线| 91麻豆精品激情在线观看国产 | 永久免费av网站大全| 免费不卡黄色视频| 国产精品久久久久成人av| 久久国产精品男人的天堂亚洲| 久久人妻熟女aⅴ| 一级片'在线观看视频| 国产在视频线精品| av在线app专区| 亚洲男人天堂网一区| 成年人黄色毛片网站| 看免费av毛片| 欧美国产精品va在线观看不卡| 国产无遮挡羞羞视频在线观看| 国产真人三级小视频在线观看| 99精国产麻豆久久婷婷| 国产淫语在线视频| 纵有疾风起免费观看全集完整版| 中国美女看黄片| 50天的宝宝边吃奶边哭怎么回事| 亚洲成人免费av在线播放| 亚洲一卡2卡3卡4卡5卡精品中文| 91九色精品人成在线观看| 男女边摸边吃奶| 免费久久久久久久精品成人欧美视频| 麻豆乱淫一区二区| 在线精品无人区一区二区三| 高清黄色对白视频在线免费看| svipshipincom国产片| 免费看av在线观看网站| 男女床上黄色一级片免费看| av片东京热男人的天堂| 在线 av 中文字幕| 国产精品二区激情视频| 女警被强在线播放| 在线 av 中文字幕| 国产不卡av网站在线观看| 亚洲一卡2卡3卡4卡5卡精品中文| www日本在线高清视频| 国产伦理片在线播放av一区| 亚洲天堂av无毛| 美女福利国产在线| 欧美精品亚洲一区二区| www.999成人在线观看| 别揉我奶头~嗯~啊~动态视频 | 国产精品一国产av| 亚洲欧美一区二区三区黑人| 欧美日韩福利视频一区二区| 国产免费视频播放在线视频| 天天影视国产精品| 亚洲天堂av无毛| av不卡在线播放| 午夜免费成人在线视频| 一级,二级,三级黄色视频| 男男h啪啪无遮挡| 自线自在国产av| 国产精品久久久久久人妻精品电影 | 黄色a级毛片大全视频| 美女国产高潮福利片在线看| 黄片小视频在线播放| 最黄视频免费看| 真人做人爱边吃奶动态| 国产一区有黄有色的免费视频| 汤姆久久久久久久影院中文字幕| 国产99久久九九免费精品| 丁香六月天网| 一级片免费观看大全| 我的亚洲天堂| 亚洲国产欧美日韩在线播放| 美女高潮到喷水免费观看| 狂野欧美激情性xxxx| 在线 av 中文字幕| 久久狼人影院| 久久99精品国语久久久| 国产免费一区二区三区四区乱码| 久久天堂一区二区三区四区| 99热网站在线观看| 在现免费观看毛片| 亚洲第一青青草原| 99热国产这里只有精品6| 黄色怎么调成土黄色| 每晚都被弄得嗷嗷叫到高潮| 男女免费视频国产| 丝瓜视频免费看黄片| 日本a在线网址| 欧美黄色淫秽网站| 中文乱码字字幕精品一区二区三区| 精品亚洲成国产av| 丝袜人妻中文字幕| 99国产精品99久久久久| 国产成人系列免费观看| 亚洲色图 男人天堂 中文字幕| 久久99热这里只频精品6学生| 日韩电影二区| 少妇人妻久久综合中文| 国产视频一区二区在线看| 国产精品一区二区在线观看99| 国产精品一国产av| 美国免费a级毛片| 两个人免费观看高清视频| 日本色播在线视频| 欧美乱码精品一区二区三区| h视频一区二区三区| 亚洲天堂av无毛| 一边摸一边抽搐一进一出视频| 亚洲中文av在线| 国产高清国产精品国产三级| 好男人电影高清在线观看| 在线看a的网站| 亚洲色图综合在线观看| 欧美精品高潮呻吟av久久| 亚洲国产欧美一区二区综合| 久久青草综合色| 精品亚洲成国产av| 婷婷色麻豆天堂久久| av国产久精品久网站免费入址| 久久精品国产a三级三级三级| 色网站视频免费| 亚洲综合色网址| 亚洲自偷自拍图片 自拍| 久久毛片免费看一区二区三区| 人人妻人人澡人人看| 亚洲熟女毛片儿| 亚洲av电影在线进入| 中国国产av一级| av天堂在线播放| 久久99精品国语久久久| 女人爽到高潮嗷嗷叫在线视频| 亚洲国产毛片av蜜桃av| 国产色视频综合| 在线观看www视频免费| 女性生殖器流出的白浆| 男女国产视频网站| 亚洲国产毛片av蜜桃av| 国产黄频视频在线观看| 精品少妇内射三级| 精品人妻在线不人妻| 一区在线观看完整版| 亚洲午夜精品一区,二区,三区| 19禁男女啪啪无遮挡网站| 国产精品免费大片| 精品人妻1区二区| 日韩一本色道免费dvd| 99久久99久久久精品蜜桃| 黄片小视频在线播放| xxx大片免费视频| 少妇 在线观看| 波野结衣二区三区在线| 一级毛片电影观看| 丝瓜视频免费看黄片| 99精品久久久久人妻精品| 欧美日韩av久久| 成人亚洲欧美一区二区av| 大片电影免费在线观看免费| 亚洲精品第二区| 叶爱在线成人免费视频播放| 黑人猛操日本美女一级片| 极品少妇高潮喷水抽搐| 亚洲熟女精品中文字幕| 曰老女人黄片| 人人妻人人澡人人爽人人夜夜| 91老司机精品| 考比视频在线观看| 69精品国产乱码久久久| 欧美性长视频在线观看| 午夜福利在线免费观看网站| 久久国产精品人妻蜜桃| 两个人免费观看高清视频| 中文精品一卡2卡3卡4更新| 亚洲国产精品成人久久小说| 成人免费观看视频高清| 亚洲欧洲精品一区二区精品久久久| 欧美av亚洲av综合av国产av| 香蕉国产在线看| 在线观看免费高清a一片| 老汉色∧v一级毛片| 两性夫妻黄色片| 国产高清国产精品国产三级| 久久久国产欧美日韩av| 国产成人av教育| 欧美日韩黄片免| 青草久久国产| 悠悠久久av| 免费看十八禁软件| 人妻人人澡人人爽人人| 久久人妻福利社区极品人妻图片 | 91成人精品电影| 亚洲国产av影院在线观看| 亚洲专区中文字幕在线| 欧美精品一区二区大全| 高清欧美精品videossex| 91精品三级在线观看| 国产视频首页在线观看| 在线观看免费视频网站a站| av网站在线播放免费| 丝袜喷水一区| 久久久精品国产亚洲av高清涩受| 久久国产精品影院| 日韩av在线免费看完整版不卡| 国产精品久久久久成人av| 成在线人永久免费视频| 亚洲色图综合在线观看| 熟女av电影| 两个人看的免费小视频| 欧美黑人精品巨大| 日韩中文字幕视频在线看片| 亚洲成人免费av在线播放| 国产一区二区在线观看av| 18禁观看日本| 老鸭窝网址在线观看| 国产精品久久久久久精品电影小说| 大型av网站在线播放| 王馨瑶露胸无遮挡在线观看| 久久天堂一区二区三区四区| 成人影院久久| 午夜影院在线不卡| 一边亲一边摸免费视频| 精品人妻1区二区| 桃花免费在线播放| 少妇精品久久久久久久| 成人国语在线视频| 亚洲国产毛片av蜜桃av| 最新在线观看一区二区三区 | 日本五十路高清| 国产成人91sexporn| 999久久久国产精品视频| 69精品国产乱码久久久| 国产爽快片一区二区三区| 精品国产一区二区久久| 成年av动漫网址| 两人在一起打扑克的视频| 女人高潮潮喷娇喘18禁视频| 爱豆传媒免费全集在线观看| 欧美黑人精品巨大| 在线观看国产h片| 日本黄色日本黄色录像| www.熟女人妻精品国产| 欧美老熟妇乱子伦牲交| 免费在线观看日本一区| 女人久久www免费人成看片| 欧美日韩国产mv在线观看视频| 人体艺术视频欧美日本| 男人爽女人下面视频在线观看| a级毛片在线看网站| 亚洲精品第二区| 波多野结衣一区麻豆|