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

    基于轉(zhuǎn)錄組測序花煙草響應(yīng)鹽脅迫活性氧清除相關(guān)基因的挖掘

    2020-12-21 09:19:40魯琳楊尚諭劉維東魯黎明
    生物技術(shù)通報 2020年12期
    關(guān)鍵詞:谷胱甘肽活性氧煙草

    魯琳 楊尚諭 劉維東 魯黎明

    (1. 四川農(nóng)業(yè)大學(xué)風(fēng)景園林學(xué)院,成都 611130;2. 四川農(nóng)業(yè)大學(xué)農(nóng)學(xué)院,成都 611130)

    土壤鹽漬化是威脅植物正常生長的最重要因素之一,并對農(nóng)業(yè)生產(chǎn)造成巨大威脅[1]。在高鹽脅迫下,Na+所產(chǎn)生的滲透脅迫和離子毒性嚴重影響了植物細胞代謝和生理功能,阻礙植物生長和發(fā)育,嚴重時會導(dǎo)致植物的死亡[2]。在長期的進化過程中,植物形成了一套復(fù)雜的抵御鹽脅迫的機制。在鹽脅迫下,植物通過活性氧(Reactive oxygen species,ROS)清除、信號轉(zhuǎn)導(dǎo)和下游基因表達,以及啟動過度鹽敏性 (Salt overly sensitive,SOS)途徑等來調(diào)節(jié)其形態(tài)、生理和生化過程,以抵御外界高鹽的環(huán)境[3]。在這些機制中,ROS清除系統(tǒng)在植物對鹽脅迫的適應(yīng)中起著重要的作用。

    植物在高鹽、干旱、寒冷、高溫和遭受病原菌的侵染時,均會產(chǎn)生大量的ROS,而過量的ROS,會被植物細胞的抗氧化防御系統(tǒng)所清除[4]。研究表明,超氧化物歧化酶(Superoxide dismutase,SOD)、過氧化物酶(Peroxidase,POD)、過氧化氫酶(Catalases,CAT)、抗壞血酸過氧化物 酶(Ascorbate peroxidase,APX)、愈 創(chuàng) 木 酚 過氧化物酶(Guaiacol peroxidase,GPOX)、谷胱甘肽還原酶(Glutathione reductase,GR)、單脫氫抗壞血酸還原酶(Monodehydroascorbate reductase,MDHAR)、脫氫抗壞血酸還原酶(Dehydroascorbate reductase,DHAR)、谷胱甘肽S-轉(zhuǎn)移酶(Glutathione S-transferases,GST)和谷胱甘肽過氧化物酶(Glutathione peroxidase,GPX)等,是植物活性氧清除的關(guān)鍵酶與主要執(zhí)行者[4]。研究表明,在鹽脅迫下,煙草幼苗的SOD、POD、CAT及APX的活性均有所增強[5-6]。同時,SbpAPX在煙草中的過表達,增強了轉(zhuǎn)基因植株對鹽及干旱的抵抗能力[7]。Cu/Zn-SOD[8-9]、Mn SOD[10-11]、POD[12]、CAT[13]、APX[14-15]及GPX[16-17]的過表達植株,均表現(xiàn)出較高的對鹽脅迫的耐受能力。

    轉(zhuǎn)錄水平的調(diào)控是植物生命活動調(diào)節(jié)的重要調(diào)控形式,轉(zhuǎn)錄組測序則是了解植物基因總體表達情況的強大工具。目前,高通量轉(zhuǎn)錄組測序技術(shù)已經(jīng)在植物響應(yīng)生物及非生物脅迫的機制研究中得到了廣泛的應(yīng)用[18-20],一大批涉及防御、物質(zhì)運輸、信號轉(zhuǎn)導(dǎo)及次生代謝的基因得到鑒定,對于揭示植物抗逆的分子機制起到了很大的推動作用。

    花煙草(Nicotiana alata)是屬于茄科煙草屬的一年生野生草本植物,由于其花色艷麗多樣、花期較長,在園林綠化中得到了較多的應(yīng)用。然而,土壤鹽漬化也嚴重影響了花煙草的正常生長。目前,對花煙草鹽脅迫響應(yīng)的相關(guān)研究尚不多見。

    本研究采用轉(zhuǎn)錄組測序(RNA-Seq)技術(shù),分析花煙草響應(yīng)鹽脅迫的基因的表達特性,重點挖掘與分析了ROS清除的相關(guān)基因,以期為花煙草耐鹽的分子機制的研究,以及耐鹽花煙草品種的選育提供參考。

    1 材料與方法

    1.1 材料

    以花煙草(Nicotiana alata)為試驗材料,由四川農(nóng)業(yè)大學(xué)煙草實驗室提供。

    1.2 方法

    1.2.1 材料處理及轉(zhuǎn)錄組測序 幼苗材料培養(yǎng):挑選籽粒飽滿且大小一致的種子,用次氯酸鈉溶液表面消毒后,用無菌水清洗數(shù)次,隨后播種于MS培養(yǎng)基上,并置于恒溫培養(yǎng)室中生長(光照:16 h光/8 h暗;溫度:25℃)。

    鹽脅迫處理及總RNA提?。捍裏熋缟L4周后,挑選36株長勢一致的煙苗,分為3組,每組煙苗為1個重復(fù),每重復(fù)12株煙苗。在無菌操作臺中,將其分別移苗至含有200 mmol/L NaCl的MS培養(yǎng)基上,然后,繼續(xù)在恒溫培養(yǎng)室中進行培養(yǎng)(光照:16 h光/8 h暗;溫度:25℃)。在處理0 h和12 h后,按重復(fù)進行整株取樣,并置于液氮中保存,用于總RNA的提取。參照魯黎明等[21]方法提取總RNA,并經(jīng)紫外分光光度計和瓊脂糖凝膠電泳檢測合格后,送北京諾禾致源科技公司進行轉(zhuǎn)錄組測序。

    1.2.2 測序數(shù)據(jù)的分析及功能注釋 對測序所獲得的原始數(shù)據(jù),按照諾禾致源公司的數(shù)據(jù)處理方法,進行原始數(shù)據(jù)的過濾,獲得高質(zhì)量的clean reads,并進行后續(xù)分析。采用BLAST軟件,在NCBI的non-redundant(NR)、UniProt/ Swiss-Prot、Gene Ontology(GO)、KEGG(Kyoto Encyclopedia of Genes and Gnomes)等數(shù)據(jù)庫進行比對,并進行差異基因的功能注釋、GO及KEGG分析。

    差異基因(Differentially expressed genes,DEGs)的篩選條件為P-value≤0.01,| log2Fold Changes|≥1。其中,F(xiàn)old Changes代表不同處理之間基因表達量的比值。該值大于1時,則該基因認定為上調(diào)表達;反之,則為下調(diào)表達。

    1.2.3 差異表達基因的驗證 為了對轉(zhuǎn)錄組測序的結(jié)果進行驗證,隨機挑選8個基因,并利用Primer 5.0軟件設(shè)計擴增引物(表1),采用qPCR的方法,對其表達水平進行驗證。將各處理的樣品總RNA反轉(zhuǎn)錄合成cDNA,以此為模板,進行qPCR反應(yīng),內(nèi)參為β-actin。qPCR的反應(yīng)體系、反應(yīng)程序以及基因相對表達量的計算均參照魯黎明等[21]方法進行。

    表1 qPCR反應(yīng)所用引物

    2 結(jié)果

    2.1 測序質(zhì)量分析

    采用高通量測序平臺,對低溫處理組及對照組的樣品進行了RNA測序。結(jié)果顯示,共獲得47 599 094-62 591 696 clean reads(表2)。就每個樣品的GC含量來看,其含量為42.97%-43.19%,偏差很小,而且基本呈隨機分布,且Q30的數(shù)值均大于90%,說明本次測序的質(zhì)量非常可靠。

    表2 各樣品測序質(zhì)量

    通過對樣品重復(fù)間數(shù)據(jù)的相關(guān)性分析,結(jié)果(圖1)表明,重復(fù)之間的決定系數(shù)R2最低為0.962,最高則達到了0.978,說明重復(fù)之間數(shù)據(jù)的重現(xiàn)性較好,可用于下一步的分析。

    2.2 差異表達基因統(tǒng)計分析

    將鹽脅迫處理12 h與0 h樣品的轉(zhuǎn)錄組數(shù)據(jù)進行對比分析,以|log2fold changes|≥1,以及P-value≤0.01為篩選條件,篩選差異表達基因。結(jié)果表明,高鹽(200 mmol/L NaCl)處理12 h后,差異表達基因的數(shù)量達到了7 239個,其中,上調(diào)表達基因4 077個,下調(diào)表達基因3 162個。

    2.3 差異表達基因的GO功能注釋分析

    對差異基因所進行的GO分析結(jié)果顯示,上調(diào)差異表達基因的功能主要集中于生物過程(Biological process)、細胞組分(Cell component)以及分子功能(Molecular function)方面。在生物過程中,差異表達基因所富集的GO的數(shù)目達108個;細胞組分以及分子功能則分別富集到7和32個GO條目(P-value≤0.01)。下調(diào)差異表達基因富集到12個GO條目,涉及到生物過程及分子功能2個部分。

    2.3.1 上調(diào)差異表達基因富集的GO 在上調(diào)差異表達基因所富集的108個GO中,差異表達基因數(shù)量超過500個的有19個GO條目。其中,Biological process最大,包括差異表達基因4 585個(76.97%),其次為Metabolic process,包括3 660個差異表達基因(61.31%)。富集上調(diào)表達基因數(shù)目較多的屬于Biological process的前30個GO條目的具體信息見表3。

    圖1 樣品各重復(fù)之間的相關(guān)性分析

    富集到Cell component部分的差異表達基因,可以歸為7個GO條目,分別為Membrane(1 747,29.26%)、Photosystem(94,1.57%)、Oxidoreductase complex(48,0.80%)、Photosystem I(30,0.50%)、Chloroplast(29,0.49%)、Plastid(29,0.49%) 和Photosystem II oxygen evolving complex(25,0.42%)。其中,有48個上調(diào)差異表達基因與氧化還原酶蛋白復(fù)合體的形成有關(guān)。

    在分子功能方面,上調(diào)差異表達基因富集到32個GO條目,其中,催化活性(Catalytic activity)基因數(shù)目最多(3 314,55.51%),其次為氧化還原酶活性(Oxidoreductase activity)基因(852,14.27%)(表4)。值得注意的是,參與氧化還原酶活性的差異基因達1 520個,占比25.46%。

    2.3.2 下調(diào)差異表達基因富集的GO 鹽脅迫下,煙草顯著下調(diào)表達的基因主要富集在生物過程(Biological process)與分子功能(Molecular function)(表5)。前者包括5個GO條目,290個基因;后者則含有7個GO條目,270個基因。生物過程的5個GO條目,全部與離子的轉(zhuǎn)運有關(guān);分子功能的GO條目包括離子轉(zhuǎn)運活性及氧化還原活性兩類。此結(jié)果顯示,花煙草在遭受高鹽逆境時,其顯著下調(diào)表達基因的數(shù)目及富集到的GO數(shù)目,均遠遠低于上調(diào)基因。說明花煙草在響應(yīng)高鹽脅迫時,是以上調(diào)表達基因所參與的生理過程為主,同時也暗示,其離子轉(zhuǎn)運及氧化還原的生理過程可能發(fā)揮較為關(guān)鍵的作用。

    2.4 差異表達基因的KEGG通路分析

    通過對差異表達基因進行KEGG通路分析(圖2)。上調(diào)差異表達基因顯著富集到7條主要的代謝通路(q value ≤0.05)。富集因子最高的為Glycerolipid metabolism途徑,其次為Glycolysis/ Gluconeogenesis和Carbon fixation in photosynthetic organisms途徑。而富集差異基因最多的途徑為Biosynthesis of secondary metabolites,富集到543個差異表達基因,其次為Carbon metabolism途徑,有162個差異表達基因被富集到該途徑。

    表3 上調(diào)差異表達基因的GO功能注釋分析(生物過程)

    下調(diào)差異表達基因顯著富集的通路數(shù)量只有5條(q value≤0.05),其中,富集因子最高的通路為Circadian rhythm-plant。在富集基因的數(shù)量方面,這5條通路所富集差異表達基因的數(shù)量,較上調(diào)基因少了很多。富集下調(diào)差異基因最多的通路為Pyrimidine metabolism,僅僅富集了58個基因。其次,為Ribosome biogenesis in eukaryotes途徑,富集的差異表達基因數(shù)為57個。

    從富集差異表達基因的代謝途徑數(shù)量的角度也可以看出,高鹽處理后,花煙草代謝途徑中,上調(diào)基因富集的代謝途徑數(shù)量遠遠高于下調(diào)基因。從而在某種程度上說明,花煙草在遭受高鹽脅迫后,其生理響應(yīng)過程主要以上調(diào)相關(guān)基因的表達為主。

    2.5 花煙草響應(yīng)高鹽脅迫的抗氧化基因分析

    植物在遭受高鹽、干旱等非生物及生物脅迫時,均會產(chǎn)生ROS,過量的ROS會對細胞的正常生命活動造成阻遏。因此,細胞會動員其抗氧化系統(tǒng)來清除這些ROS。為了探索花煙草抗氧化系統(tǒng)關(guān)鍵基因的表達對高鹽脅迫的響應(yīng),分析了差異表達基因中SOD、POD、CAT等7類ROS清除基因的表達。結(jié)果表明,在花煙草受到高鹽脅迫后,共有45個ROS清除相關(guān)基因的表達發(fā)生了顯著改變,其中,上調(diào)表達基因28個,下調(diào)表達基因17個。上調(diào)表達的倍數(shù)最小為1.10倍(GST),最高為3.10倍(CAT);下調(diào)表達的倍數(shù)最小為1.14倍(GST),最高為3.32倍(POD43)。上調(diào)表達基因中,數(shù)量最多的為GST(16),其 次 為POD(5)、SOD(4)、CAT(2)和GPX(1)(圖3)。下調(diào)表達基因中,數(shù)量最多的為GST(6),其次為POD(5)、SOD(2)、GPX(2)、APX(1)和GR(Glutathione reductase,2)。綜上所述,花煙草在高鹽脅迫下,其抗氧化的機制主要是上調(diào)抗氧化基因的表達,以啟動相應(yīng)的生理反應(yīng),從而抵御不利的環(huán)境。

    表4 上調(diào)差異表達基因的GO功能注釋分析(分子功能)

    表5 下調(diào)差異表達基因的GO功能注釋分析

    2.6 qRT-PCR驗證

    為了驗證RNA測序結(jié)果的可靠性,隨機選取了POD42、GST1、CAT1、CIPK11、ERF4、GPX4、MYB5和CIPK6等8個基因進行qPCR擴增,并與其測序結(jié)果進行對比(圖4)。從圖4可以看出,2種方法中,所選擇的8個基因的表達趨勢一致,從而說明本研究RNA測序結(jié)果較為可靠,同時也表明,對RNA測序結(jié)果所做的相關(guān)分析較為客觀。

    3 討論

    鹽漬化是植物在生長過程中所遭受到的主要的環(huán)境脅迫之一。在高鹽脅迫下,植物細胞會對自己的生理活動進行調(diào)節(jié),以抵御其對植物生長的影響。其中,活性氧清除系統(tǒng)對鹽脅迫所誘導(dǎo)產(chǎn)生的ROS的清除,是主要的生理響應(yīng)過程[4]。高等植物的酶促活性氧清除系統(tǒng)主要包括SOD、POD、CAT、APX、GPX、GST及GR等,其主要的功能,就是將鹽脅迫下在葉綠體、線粒體等細胞器中產(chǎn)生的過量的ROS進行還原,以降低其產(chǎn)生的危害。在本研究中,花煙草在鹽脅迫下,有7類共45個活性氧清除基因的表達發(fā)生了顯著變化(P <0.05),暗示鹽脅迫將導(dǎo)致花煙草眾多的生理進程發(fā)生改變。

    3.1 SOD、POD和CAT

    SOD、POD及CAT是植物最為常見的ROS清除酶類。在逆境脅迫下,植物編碼這3類酶基因表達的提高,酶活性的增強,對ROS的清除能力也大幅度提升,從而提高了植物對不良環(huán)境的抵御能力[4]。研究表明,抗鹽能力較強的基因型,其抗氧化的能力也強,在遭受高鹽脅迫時,其長勢及生物量均較不耐鹽的基因型有較大的優(yōu)勢[1]。同時,當SOD、POD和CAT在植物中過表達時,轉(zhuǎn)基因植株的耐鹽能力均顯著增強[8-13]。在本研究中,當花煙草面臨高鹽脅迫時,編碼SOD、POD及CAT 3個酶類的11個基因的發(fā)生了顯著的上調(diào)表達變化。其中,上調(diào)表達最高為3.1倍(CAT)。揭示花煙草響應(yīng)鹽脅迫過程中,SOD、POD和CAT發(fā)揮著十分重要的作用,其中,CAT的作用值得關(guān)注。

    然而,在本研究中,也有5個編碼POD的基因及2個編碼SOD的基因發(fā)生了下調(diào)表達。這種現(xiàn)象產(chǎn)生的意義需要進一步的研究,同時,此結(jié)果也說明,花煙草對鹽脅迫的響應(yīng),可能存在正向與負向兩種調(diào)控方式,其中,正向調(diào)節(jié)作用可能占主導(dǎo)地位。

    3.2 APX與GPX

    圖2 差異表達基因的KEGG分析

    抗壞血酸過氧化物酶(Ascorbate peroxidase,APX)及谷胱甘肽過氧化物酶(Glutathione peroxidase,GPX)是2個構(gòu)成植物抗氧化系統(tǒng)的關(guān)鍵酶類。APX是一類在高等植物中廣泛存在的亞鐵血紅素蛋白,包括了4個不同的亞型,即葉綠體APX(sAPX)、類囊體APX(tAPX)、微體APX(mbAPX)和胞質(zhì)APX(cAPX)[22]。APX的表達受生物及非生物脅迫因子的廣泛誘導(dǎo)。研究表明,在鹽脅迫下,堿蓬的SsAPX[23]及白樺的APX[24]的表達均顯著升高;擬南芥APX[25]及棗樹ZjAPX[26]的過表達,則顯著增強了轉(zhuǎn)基因植株的耐鹽能力。

    圖3 花煙草響應(yīng)高鹽脅迫的抗氧化基因

    圖4 RNA測序與qPCR基因表達水平的比較

    GPX也是細胞ROS清除的重要酶類,同樣是由多個同工酶組成的大家族。GPX能夠有效地清除細胞過氧化物,防止ROS對細胞的傷害,并因而在植物生長發(fā)育及逆境脅迫響應(yīng)中發(fā)揮重要作用[27]。例如,衣藻(Chlamydomonas)GPX[28]、小麥GPX[29]的過表達,均增強了轉(zhuǎn)基因植株對鹽、低溫的耐受能力。然而,植物GPX的表達較為復(fù)雜,在大多數(shù)情況下,GPX以上調(diào)表達的方式參與環(huán)境脅迫的響應(yīng)[27],但高粱(Sorghum bicolor)的SbGPX6與SbGPX7[30],在干旱脅迫下卻是下調(diào)表達。

    以上結(jié)果均表明,APX與GPX均參與了植物對鹽脅迫的響應(yīng)。在本研究中,花煙草在鹽脅迫下,有1個GPX上調(diào)表達了1.27倍,另外2個GPX及1個APX均為下調(diào)表達。揭示花煙草APX及GPX參與鹽脅迫響應(yīng)有自己的特點,而且可能主要是負向調(diào)控方式。

    3.3 GST

    谷胱甘肽-S-轉(zhuǎn)移酶(GST)也是由多基因編碼的基因家族,其功能較為廣泛,參與了細胞的解毒、抗氧化、信號轉(zhuǎn)導(dǎo)等過程,因而在植物的生長發(fā)育及逆境脅迫響應(yīng)中發(fā)揮重要的作用[31]。研究表明,植物的GST參與了高鹽、低溫、干旱、UV及ABA等多種非生物脅迫的響應(yīng)[32]。例如,擬南芥AtGSTU19參與了植物的抗旱與抗鹽脅迫[33];過表達煙草[34]及堿蓬[35]的GST,能夠顯著增強轉(zhuǎn)基因植物的耐鹽能力。然而,GST參與植物的逆境脅迫響應(yīng)的方式較為復(fù)雜,如AtGSTU17可能作為負調(diào)節(jié)因子,參與了擬南芥對干旱和鹽脅迫的應(yīng)答反應(yīng),其較低的表達量或基因敲除有助于擬南芥抗性的提高[36]。

    本研究也得出了類似的結(jié)果,在高鹽脅迫下,花煙草22個GST的表達顯著增強,其中,有16個顯著上調(diào)(最高2.81倍),6個則顯著下調(diào)(最高2.22倍)(圖3)。表明在花煙草的高鹽脅迫應(yīng)答中,GST占有舉足輕重的地位,同時,其不同成員間發(fā)揮作用的機制可能并不相同。

    3.4 GR

    谷胱甘肽還原酶(GR)與GST和GPX共同構(gòu)成了谷胱甘肽循環(huán)的重要組件,并且在ROS清除機制中發(fā)揮重要作用[31]。GR的基本功能是催化還原型谷胱甘肽(GSH)的合成,維持細胞內(nèi)的較高水平的還原力;同時,GR還能夠與SOD等酶系統(tǒng)共同作用,清除細胞內(nèi)多余的ROS[37]。所以,GR活性的增強,有利于植物抗氧化能力的提升,并最終提高了植物對生物及非生物逆境脅迫的抵御能力。植物GR在鹽脅迫應(yīng)答中的功能研究報道尚不多見,已有的研究表明,在鹽脅迫誘導(dǎo)下,水稻的GR3[38]、茶樹CsGR1[39]、油菜的GR2[40]、麻風(fēng)樹(Jatropha curcas)JcGR[41]等基因的表達明顯升高;而水稻GR3的敲除,則導(dǎo)致了突變體對鹽脅迫的敏感[42]。但是GR對脅迫的應(yīng)答的方式并不相同,如在鹽脅迫下,茶樹的CsGR2則顯著下降[39]。與此相類似,花煙草在高鹽環(huán)境下只有1個GR發(fā)生了下調(diào)表達(下調(diào)1.33倍)(圖3),說明在花煙草的鹽脅迫響應(yīng)中,GR也可能發(fā)揮著一定的作用,其作用的機制有待于進一步的探索。

    4 結(jié)論

    在鹽脅迫下,花煙草大量的基因發(fā)生了顯著的差異表達。這些差異表達基因的功能,主要集中在生物過程及金屬離子轉(zhuǎn)運等方面。其中,有7大類45個活性氧清除相關(guān)基因的表達模式變化明顯。本研究的結(jié)果說明,花煙草通過改變自身眾多的生理進程,以適應(yīng)外界的鹽脅迫環(huán)境。其中,增強對鹽脅迫所產(chǎn)生的活性氧的清除能力,是花煙草抵御鹽脅迫的主要機制之一。

    猜你喜歡
    谷胱甘肽活性氧煙草
    煙草具有輻射性?
    煙草依賴的診斷標準
    蚯蚓谷胱甘肽轉(zhuǎn)硫酶分離純化的初步研究
    煙草中茄酮的富集和應(yīng)用
    分光光度法結(jié)合抗干擾補償檢測谷胱甘肽
    TLR3活化對正常人表皮黑素細胞內(nèi)活性氧簇表達的影響
    TBBPA對美洲鰻鱺谷胱甘肽代謝相關(guān)指標的影響
    硅酸鈉處理對杏果實活性氧和苯丙烷代謝的影響
    煙草鏡頭與歷史真實
    聲屏世界(2014年6期)2014-02-28 15:18:09
    O2聯(lián)合CO2氣調(diào)對西蘭花活性氧代謝及保鮮效果的影響
    成年人午夜在线观看视频| 永久免费av网站大全| 亚洲美女视频黄频| 亚洲精品影视一区二区三区av| 热re99久久精品国产66热6| 免费少妇av软件| 国产黄片美女视频| 成人国产麻豆网| 亚洲精品国产成人久久av| 毛片女人毛片| 丝袜喷水一区| 欧美三级亚洲精品| 成人国产av品久久久| 赤兔流量卡办理| 国内精品美女久久久久久| 亚洲国产精品专区欧美| 亚洲性久久影院| 国产高潮美女av| 成人国产av品久久久| 国产精品.久久久| 精品人妻视频免费看| 精品一区在线观看国产| 青春草国产在线视频| a级一级毛片免费在线观看| 国产午夜福利久久久久久| 精品久久久久久电影网| 国产精品久久久久久av不卡| 精品亚洲乱码少妇综合久久| 夫妻性生交免费视频一级片| 日韩中字成人| 精品国产三级普通话版| 中国美白少妇内射xxxbb| av国产精品久久久久影院| 亚洲国产欧美在线一区| 一级毛片久久久久久久久女| 亚洲国产色片| 国产精品一二三区在线看| 日本一本二区三区精品| 国产日韩欧美亚洲二区| 少妇裸体淫交视频免费看高清| 精品少妇久久久久久888优播| 麻豆国产97在线/欧美| 国产精品秋霞免费鲁丝片| 婷婷色麻豆天堂久久| 波野结衣二区三区在线| 国产 精品1| 日本欧美国产在线视频| 久久97久久精品| 嫩草影院精品99| 精品一区在线观看国产| 免费av不卡在线播放| 国产爽快片一区二区三区| 超碰av人人做人人爽久久| 亚洲国产av新网站| 免费不卡的大黄色大毛片视频在线观看| 精品熟女少妇av免费看| 最近最新中文字幕大全电影3| 一级a做视频免费观看| 国产精品成人在线| 欧美老熟妇乱子伦牲交| 亚洲va在线va天堂va国产| 亚洲综合精品二区| 少妇的逼好多水| 搡老乐熟女国产| 国产一区二区亚洲精品在线观看| 天堂中文最新版在线下载 | 亚洲欧美一区二区三区黑人 | 人妻一区二区av| 欧美日韩亚洲高清精品| 少妇猛男粗大的猛烈进出视频 | av在线亚洲专区| 成人免费观看视频高清| 大话2 男鬼变身卡| 少妇高潮的动态图| 午夜精品一区二区三区免费看| 少妇熟女欧美另类| 国产精品熟女久久久久浪| 亚洲精品久久午夜乱码| 久久久国产一区二区| 成人亚洲欧美一区二区av| 两个人的视频大全免费| 99精国产麻豆久久婷婷| 极品少妇高潮喷水抽搐| 又爽又黄a免费视频| 别揉我奶头 嗯啊视频| tube8黄色片| 性色av一级| 国产精品一二三区在线看| 国产精品久久久久久av不卡| 国语对白做爰xxxⅹ性视频网站| 国产综合精华液| 在线观看美女被高潮喷水网站| 女人十人毛片免费观看3o分钟| 亚洲av免费高清在线观看| 最近2019中文字幕mv第一页| 麻豆乱淫一区二区| 国产成人精品一,二区| 国产精品国产av在线观看| 日本三级黄在线观看| 国产成人精品一,二区| 91精品一卡2卡3卡4卡| 99久久精品一区二区三区| 日韩电影二区| 成年女人在线观看亚洲视频 | 亚洲精品亚洲一区二区| 啦啦啦啦在线视频资源| av女优亚洲男人天堂| 国产精品蜜桃在线观看| 最近2019中文字幕mv第一页| 久久精品久久精品一区二区三区| 卡戴珊不雅视频在线播放| 啦啦啦在线观看免费高清www| 色哟哟·www| 另类亚洲欧美激情| 国产综合懂色| 久久精品熟女亚洲av麻豆精品| 精品国产露脸久久av麻豆| 日产精品乱码卡一卡2卡三| 国产av不卡久久| 国模一区二区三区四区视频| 亚洲国产成人一精品久久久| 久久精品国产a三级三级三级| 国产一区有黄有色的免费视频| 日本熟妇午夜| 免费观看在线日韩| av国产免费在线观看| 自拍欧美九色日韩亚洲蝌蚪91 | 尾随美女入室| 国内少妇人妻偷人精品xxx网站| 中国三级夫妇交换| 欧美成人a在线观看| 色视频www国产| 国产 一区 欧美 日韩| 久久鲁丝午夜福利片| 欧美亚洲 丝袜 人妻 在线| 18禁裸乳无遮挡动漫免费视频 | 日日摸夜夜添夜夜添av毛片| 天堂中文最新版在线下载 | 久久久久久久精品精品| 国产成人a∨麻豆精品| 各种免费的搞黄视频| 能在线免费看毛片的网站| 免费黄网站久久成人精品| 欧美人与善性xxx| 看黄色毛片网站| 永久网站在线| 国产一区二区三区av在线| 国产成人福利小说| av国产精品久久久久影院| 亚洲精品日本国产第一区| 少妇人妻精品综合一区二区| 日日啪夜夜爽| 色5月婷婷丁香| 人人妻人人澡人人爽人人夜夜| 在线a可以看的网站| 中文字幕av成人在线电影| 26uuu在线亚洲综合色| 狂野欧美激情性xxxx在线观看| 国产免费一区二区三区四区乱码| 又爽又黄无遮挡网站| 精品久久久久久久人妻蜜臀av| 不卡视频在线观看欧美| 亚洲av中文av极速乱| 男女啪啪激烈高潮av片| 国产欧美日韩一区二区三区在线 | 亚洲国产精品国产精品| 国产精品一区二区在线观看99| av国产免费在线观看| 少妇猛男粗大的猛烈进出视频 | 亚洲欧美精品专区久久| 国产亚洲91精品色在线| 国产高清国产精品国产三级 | 嫩草影院新地址| 一区二区三区免费毛片| 国产精品熟女久久久久浪| 熟妇人妻不卡中文字幕| 亚洲国产精品999| 国产探花极品一区二区| av网站免费在线观看视频| 国产精品av视频在线免费观看| 色哟哟·www| 日本一二三区视频观看| 午夜精品国产一区二区电影 | 在线天堂最新版资源| 黄片无遮挡物在线观看| 国产亚洲91精品色在线| 亚洲婷婷狠狠爱综合网| 秋霞在线观看毛片| 久久久久久久久久久丰满| 国产成人午夜福利电影在线观看| 伊人久久国产一区二区| 国产成人精品一,二区| 国产在线男女| 欧美xxⅹ黑人| 国产男人的电影天堂91| 国产成人免费无遮挡视频| 啦啦啦啦在线视频资源| 国产精品不卡视频一区二区| av天堂中文字幕网| 亚洲精品久久久久久婷婷小说| 成人特级av手机在线观看| 女人十人毛片免费观看3o分钟| 欧美精品国产亚洲| 69人妻影院| 51国产日韩欧美| 久久精品国产亚洲网站| 日韩电影二区| 99精国产麻豆久久婷婷| 亚洲综合色惰| 久久6这里有精品| 交换朋友夫妻互换小说| 18禁裸乳无遮挡动漫免费视频 | 免费在线观看成人毛片| 在线亚洲精品国产二区图片欧美 | 91精品一卡2卡3卡4卡| 网址你懂的国产日韩在线| av免费观看日本| 久久精品综合一区二区三区| 丝袜美腿在线中文| 国内精品宾馆在线| 国产乱人偷精品视频| 尾随美女入室| 欧美精品人与动牲交sv欧美| www.色视频.com| 免费大片18禁| 国产淫片久久久久久久久| av黄色大香蕉| 成人亚洲欧美一区二区av| 最新中文字幕久久久久| 亚洲婷婷狠狠爱综合网| 特大巨黑吊av在线直播| 国产伦精品一区二区三区四那| 欧美日韩综合久久久久久| 欧美日韩精品成人综合77777| 亚洲av一区综合| 又爽又黄无遮挡网站| 在线精品无人区一区二区三 | 午夜福利在线观看免费完整高清在| 亚洲综合色惰| 免费观看无遮挡的男女| 国产中年淑女户外野战色| 内地一区二区视频在线| 成人高潮视频无遮挡免费网站| av免费观看日本| 少妇猛男粗大的猛烈进出视频 | 亚洲精品成人久久久久久| 天天躁日日操中文字幕| 一本一本综合久久| 国产成人精品一,二区| 国产成人精品久久久久久| 国产精品久久久久久精品电影| 王馨瑶露胸无遮挡在线观看| 九九爱精品视频在线观看| 国产老妇女一区| 国产成人免费观看mmmm| 美女cb高潮喷水在线观看| 观看免费一级毛片| 白带黄色成豆腐渣| 国产精品久久久久久精品电影| 国产美女午夜福利| 欧美97在线视频| 日韩一本色道免费dvd| 校园人妻丝袜中文字幕| 超碰97精品在线观看| 大陆偷拍与自拍| 各种免费的搞黄视频| 国产 精品1| av在线app专区| 80岁老熟妇乱子伦牲交| 内射极品少妇av片p| 在线a可以看的网站| 亚洲精品亚洲一区二区| 亚洲精品国产成人久久av| 亚洲欧美中文字幕日韩二区| 亚洲av中文字字幕乱码综合| 日产精品乱码卡一卡2卡三| 天天躁日日操中文字幕| 国产爽快片一区二区三区| 成年女人看的毛片在线观看| 久久久国产一区二区| 黄色视频在线播放观看不卡| 搞女人的毛片| 亚洲综合色惰| 国产亚洲av嫩草精品影院| 国产综合懂色| 欧美三级亚洲精品| 久久久国产一区二区| 黄色视频在线播放观看不卡| 中国美白少妇内射xxxbb| 亚洲欧美成人综合另类久久久| 久久韩国三级中文字幕| 一级av片app| 看黄色毛片网站| 国产亚洲5aaaaa淫片| 国产高清国产精品国产三级 | 久久精品夜色国产| 精品久久久久久久人妻蜜臀av| 亚洲欧美日韩无卡精品| 色视频在线一区二区三区| 香蕉精品网在线| 人人妻人人爽人人添夜夜欢视频 | 男插女下体视频免费在线播放| 免费少妇av软件| 在线观看av片永久免费下载| 一边亲一边摸免费视频| 91精品国产九色| 亚洲,欧美,日韩| 汤姆久久久久久久影院中文字幕| 一级二级三级毛片免费看| 精品99又大又爽又粗少妇毛片| 肉色欧美久久久久久久蜜桃 | 国产成人精品一,二区| 欧美精品人与动牲交sv欧美| 欧美国产精品一级二级三级 | 一区二区av电影网| 18禁在线无遮挡免费观看视频| 久久久久久久亚洲中文字幕| 欧美xxⅹ黑人| 亚洲av一区综合| 99热这里只有精品一区| 亚洲成色77777| 欧美国产精品一级二级三级 | 亚洲av免费在线观看| 日韩亚洲欧美综合| 午夜视频国产福利| 视频区图区小说| 国产乱来视频区| av.在线天堂| 国精品久久久久久国模美| 国产在线一区二区三区精| 亚洲经典国产精华液单| 日韩三级伦理在线观看| 成人鲁丝片一二三区免费| 成人亚洲精品一区在线观看 | 白带黄色成豆腐渣| 国产精品国产三级国产av玫瑰| 精品一区二区三区视频在线| 在线播放无遮挡| 免费av毛片视频| 国内精品宾馆在线| 亚洲在久久综合| 国产成人免费观看mmmm| 国产人妻一区二区三区在| 成人免费观看视频高清| 夫妻午夜视频| 内地一区二区视频在线| 久久99精品国语久久久| 毛片一级片免费看久久久久| 国产探花在线观看一区二区| 免费观看a级毛片全部| 婷婷色综合www| 亚洲,欧美,日韩| 日日撸夜夜添| 丝袜脚勾引网站| 九草在线视频观看| 精品久久国产蜜桃| 国产精品伦人一区二区| 国产精品.久久久| 大码成人一级视频| 国产亚洲一区二区精品| 欧美人与善性xxx| 七月丁香在线播放| 大片免费播放器 马上看| 可以在线观看毛片的网站| 看非洲黑人一级黄片| 亚洲国产精品成人久久小说| 三级男女做爰猛烈吃奶摸视频| 色视频在线一区二区三区| 国产成人午夜福利电影在线观看| 欧美精品国产亚洲| 中文字幕人妻熟人妻熟丝袜美| 欧美bdsm另类| 国产精品伦人一区二区| 2022亚洲国产成人精品| 久久国内精品自在自线图片| 国产精品不卡视频一区二区| 精品国产乱码久久久久久小说| 亚洲欧美中文字幕日韩二区| 五月开心婷婷网| 国产成年人精品一区二区| 麻豆精品久久久久久蜜桃| 欧美日韩一区二区视频在线观看视频在线 | 精品国产一区二区三区久久久樱花 | 九九在线视频观看精品| 欧美zozozo另类| 欧美人与善性xxx| 大码成人一级视频| 日本黄色片子视频| 交换朋友夫妻互换小说| 久久精品久久精品一区二区三区| 秋霞伦理黄片| 亚洲av成人精品一二三区| 熟女人妻精品中文字幕| 亚洲成色77777| 久久久久久国产a免费观看| 国产日韩欧美亚洲二区| 只有这里有精品99| 久久久久网色| 人妻制服诱惑在线中文字幕| 在线观看三级黄色| 久久精品国产亚洲av天美| 中文天堂在线官网| a级毛片免费高清观看在线播放| 蜜桃亚洲精品一区二区三区| 欧美区成人在线视频| 久久久成人免费电影| 日韩 亚洲 欧美在线| 精品国产一区二区三区久久久樱花 | 交换朋友夫妻互换小说| 免费观看的影片在线观看| 久久久久国产网址| 在线观看三级黄色| 成人亚洲精品av一区二区| 亚洲精品视频女| 国产淫片久久久久久久久| 亚洲国产高清在线一区二区三| 嫩草影院精品99| 中文字幕人妻熟人妻熟丝袜美| 国产精品偷伦视频观看了| 亚洲四区av| 看十八女毛片水多多多| 校园人妻丝袜中文字幕| 亚洲av中文字字幕乱码综合| 菩萨蛮人人尽说江南好唐韦庄| 久久精品久久久久久久性| 一级毛片黄色毛片免费观看视频| 亚洲色图综合在线观看| 久热这里只有精品99| 草草在线视频免费看| 女的被弄到高潮叫床怎么办| 国产精品国产三级专区第一集| 毛片女人毛片| 日韩,欧美,国产一区二区三区| 日本猛色少妇xxxxx猛交久久| 精品亚洲乱码少妇综合久久| 黄片无遮挡物在线观看| 亚洲人与动物交配视频| 欧美精品一区二区大全| 丰满少妇做爰视频| 久久鲁丝午夜福利片| 国产高清国产精品国产三级 | 一个人看视频在线观看www免费| 国产成人精品婷婷| 国产精品蜜桃在线观看| 亚洲欧美精品专区久久| 亚洲成人久久爱视频| av天堂中文字幕网| 综合色丁香网| 看黄色毛片网站| 日韩三级伦理在线观看| 久久久久久久久久久丰满| 哪个播放器可以免费观看大片| 建设人人有责人人尽责人人享有的 | 亚洲欧美日韩另类电影网站 | 夫妻性生交免费视频一级片| 嫩草影院新地址| 国产精品一区二区在线观看99| 国产亚洲精品久久久com| 亚洲激情五月婷婷啪啪| 在线亚洲精品国产二区图片欧美 | 看黄色毛片网站| 精品少妇久久久久久888优播| av福利片在线观看| 熟女人妻精品中文字幕| 五月伊人婷婷丁香| 日本一本二区三区精品| 一区二区三区四区激情视频| 人体艺术视频欧美日本| 国产人妻一区二区三区在| 街头女战士在线观看网站| 欧美成人午夜免费资源| 国产成人精品婷婷| 欧美+日韩+精品| av国产久精品久网站免费入址| 纵有疾风起免费观看全集完整版| 又爽又黄无遮挡网站| www.色视频.com| 可以在线观看毛片的网站| 在线观看美女被高潮喷水网站| 黄色配什么色好看| 男人和女人高潮做爰伦理| 亚洲精品视频女| 中文字幕制服av| 亚洲av免费在线观看| 亚洲精品一二三| 毛片一级片免费看久久久久| 日韩欧美 国产精品| 亚洲成人一二三区av| 国产成人91sexporn| 欧美潮喷喷水| 精品少妇久久久久久888优播| 欧美激情久久久久久爽电影| 久久久色成人| 成人黄色视频免费在线看| 老师上课跳d突然被开到最大视频| 国产男女内射视频| 日韩,欧美,国产一区二区三区| 午夜福利视频1000在线观看| freevideosex欧美| 久久久久久久久久久丰满| 色哟哟·www| 人妻制服诱惑在线中文字幕| 日韩亚洲欧美综合| 亚洲天堂av无毛| 久久久a久久爽久久v久久| 综合色av麻豆| 网址你懂的国产日韩在线| 亚洲欧美中文字幕日韩二区| 黄片wwwwww| 夜夜爽夜夜爽视频| 久久99热这里只频精品6学生| 一级毛片黄色毛片免费观看视频| 青春草国产在线视频| 少妇人妻久久综合中文| 久久久欧美国产精品| 亚洲精品影视一区二区三区av| av网站免费在线观看视频| 如何舔出高潮| 五月玫瑰六月丁香| 国产精品99久久久久久久久| 国产亚洲精品久久久com| 街头女战士在线观看网站| 国产亚洲av片在线观看秒播厂| 舔av片在线| 国产老妇女一区| 麻豆乱淫一区二区| 国产精品一及| 熟女av电影| 亚洲国产最新在线播放| 2022亚洲国产成人精品| 中文字幕人妻熟人妻熟丝袜美| 搡老乐熟女国产| 精品国产一区二区三区久久久樱花 | 亚洲精品国产av成人精品| 久久99精品国语久久久| 在线观看人妻少妇| 青青草视频在线视频观看| 免费观看在线日韩| 蜜桃亚洲精品一区二区三区| 午夜激情久久久久久久| 国产成人精品福利久久| 国产精品一区www在线观看| 亚洲熟女精品中文字幕| 欧美一级a爱片免费观看看| 久久99热这里只频精品6学生| 国产精品一区二区性色av| 欧美日韩亚洲高清精品| 国产亚洲午夜精品一区二区久久 | 亚洲国产精品成人久久小说| 最新中文字幕久久久久| 小蜜桃在线观看免费完整版高清| videossex国产| 久久99精品国语久久久| 久久久久久久亚洲中文字幕| 嫩草影院新地址| 亚洲精品aⅴ在线观看| 国产精品伦人一区二区| 少妇的逼好多水| 久久久久久久久久久免费av| 亚洲精品成人av观看孕妇| 亚洲av.av天堂| 欧美+日韩+精品| 中文乱码字字幕精品一区二区三区| 少妇被粗大猛烈的视频| 午夜免费观看性视频| 大片免费播放器 马上看| 汤姆久久久久久久影院中文字幕| 久久久a久久爽久久v久久| 午夜亚洲福利在线播放| 两个人的视频大全免费| 国产久久久一区二区三区| 丰满乱子伦码专区| 国产av码专区亚洲av| 国产精品久久久久久精品古装| 国产国拍精品亚洲av在线观看| 白带黄色成豆腐渣| 人妻系列 视频| 亚洲一级一片aⅴ在线观看| 国产伦理片在线播放av一区| 热99国产精品久久久久久7| 欧美bdsm另类| 国产精品人妻久久久久久| 高清欧美精品videossex| 97人妻精品一区二区三区麻豆| 美女被艹到高潮喷水动态| 国产免费一区二区三区四区乱码| 国产精品福利在线免费观看| 欧美一级a爱片免费观看看| 国产免费视频播放在线视频| 亚洲精品久久久久久婷婷小说| 国产精品麻豆人妻色哟哟久久| 久久久久久久精品精品| 久久热精品热| 久久久a久久爽久久v久久| 亚洲,一卡二卡三卡| 男人舔奶头视频| 亚洲人成网站高清观看| 国产精品爽爽va在线观看网站| 久久久久国产网址| 免费观看的影片在线观看| 欧美bdsm另类| 91精品一卡2卡3卡4卡| 国产精品久久久久久av不卡| av在线蜜桃| 天堂俺去俺来也www色官网| 亚洲av福利一区| 99久久九九国产精品国产免费| 日韩一区二区三区影片| 中文精品一卡2卡3卡4更新| 波野结衣二区三区在线| 亚洲成人一二三区av| 欧美国产精品一级二级三级 | 一个人看的www免费观看视频| 韩国高清视频一区二区三区| 一边亲一边摸免费视频|