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

    高粱全基因組NRT1基因鑒定、表達(dá)與DNA變異分析

    2021-11-05 02:11:30郭志強(qiáng)梁月秀朱立勛范佳麗姜曉東呂晉慧張春來(lái)
    激光生物學(xué)報(bào) 2021年5期
    關(guān)鍵詞:信號(hào)肽進(jìn)化樹(shù)硝酸鹽

    郭志強(qiáng),梁月秀,馮 凡,朱立勛,范佳麗,姜曉東,呂晉慧,張春來(lái)

    (山西農(nóng)業(yè)大學(xué)農(nóng)學(xué)院,國(guó)家功能雜糧技術(shù)創(chuàng)新中心,山西省旱作栽培與作物生態(tài)重點(diǎn)實(shí)驗(yàn)室,山西省黃土高原特色作物高效生產(chǎn)協(xié)同創(chuàng)新中心,太谷 030801)

    高粱(Sorghum bicolor)是世界上僅次于小麥、水稻、玉米和大麥的第五大糧食作物[1],主要種植在干旱和半干旱地區(qū),屬于短日照C4植物,具有抗旱、耐貧瘠和耐鹽堿等特性[2],栽培面積較為廣泛。高粱基因組DNA測(cè)序工作初步在2009年完成,并于最近得到更新,這使得從基因組水平來(lái)揭示高粱重要基因家族的功能成為可能[3-4]。

    植物氮利用率(nitrogen use efficiency,NUE)取決于植物的需求和植物對(duì)氮的吸收、運(yùn)輸、積累和再分配。氮素作為植物生長(zhǎng)發(fā)育過(guò)程中的必要元素之一,是植物體內(nèi)生物大分子如蛋白質(zhì)、葉綠素、核酸等的組成成分[5]。此外,氮素也是影響作物產(chǎn)量的因素之一[6]。氮供應(yīng)水平影響植物葉片顏色、大小以及葉片數(shù),進(jìn)而影響植物對(duì)光的吸收能力,最終影響植物的生產(chǎn)力,特別是對(duì)于多葉蔬菜[6]。硝酸鹽(NO3-)是大多數(shù)植物中最重要的氮源。提高硝酸鹽吸收能力和植物氮利用率的重要步驟是識(shí)別負(fù)責(zé)硝酸鹽運(yùn)輸?shù)倪\(yùn)輸者的特征,然后通過(guò)傳統(tǒng)的育種技術(shù)或基因工程操作來(lái)提高硝酸鹽吸收特性。近年來(lái),對(duì)擬南芥、水稻等植物的營(yíng)養(yǎng)研究表明,植物吸收轉(zhuǎn)運(yùn)硝酸鹽的過(guò)程主要由4個(gè)蛋白家族成員介導(dǎo)完成,分別為硝酸鹽轉(zhuǎn)運(yùn)蛋白1(nitrate transporter 1,NRT1)、硝酸鹽轉(zhuǎn)運(yùn)蛋白2(nitrate transporter 2,NRT2)、氯酸鹽通道家族(chloride channel,CLC)以及慢離子相關(guān)通道同系物(slow anion channel- as- sociated homolog,SLAC/SLAH)[7]。除此之外,還有許多蛋白酶、激素等參與,它們相互作用構(gòu)成了錯(cuò)綜復(fù)雜的調(diào)控網(wǎng)絡(luò)[8-10]。

    硝酸鹽轉(zhuǎn)運(yùn)蛋白(nitrate transporter,NRT),又稱(chēng)為寡肽轉(zhuǎn)運(yùn)蛋白(peptide transport,PTR),在動(dòng)物、植物、酵母、古生菌以及細(xì)菌基因組中都有發(fā)現(xiàn),是一種依賴(lài)于質(zhì)子、具有專(zhuān)一性、依賴(lài)能量的轉(zhuǎn)運(yùn)蛋白。蛋白的質(zhì)子依賴(lài)性PTR家族與ATP結(jié)合轉(zhuǎn)運(yùn)蛋白(ATP-binding cassette,ABC)不同,是通過(guò)對(duì)最近發(fā)現(xiàn)的一些肽轉(zhuǎn)運(yùn)蛋白進(jìn)行序列分析而發(fā)現(xiàn)的[11]。這些蛋白質(zhì)似乎主要參與小肽的攝取,同時(shí)還攝取質(zhì)子[12]。

    一些NRT也執(zhí)行肽的運(yùn)輸,有的還運(yùn)輸抗生素,它們通過(guò)質(zhì)子共運(yùn)發(fā)揮作用,但底物對(duì)H+的化學(xué)計(jì)量學(xué)是可變的,如高親和力的大鼠PepT2載體分別催化中性和陰離子二肽攝取兩個(gè)和三個(gè) H+,而低親和力的PepT1載體催化每個(gè)中性肽攝取一個(gè)H+[13]。

    本研究對(duì)SbNRT1的基因結(jié)構(gòu)、系統(tǒng)進(jìn)化樹(shù)、基因表達(dá)、選擇進(jìn)化及其蛋白質(zhì)的理化性質(zhì)、亞細(xì)胞定位、二級(jí)和三級(jí)結(jié)構(gòu)進(jìn)行了分析,為NRT1基因調(diào)控提供一定的理論基礎(chǔ)與試驗(yàn)依據(jù)。

    1 材料與方法

    1.1 基因檢索

    以水稻(Oryza sativa)和擬南芥(Arabidopsis thaliana)的NRT1氨基酸序列作為參考,在Phytozome(https://phytozome.jgi.doe.gov/pz/portal.html)和 NCBI(http://www.ncbi.nlm.nih.gov/)數(shù)據(jù)庫(kù)中進(jìn)行Blast搜尋和比對(duì),獲得SbNRT1基因家族的基因序列、編碼序列(coding sequence,CDS)、氨基酸序列。小麥(Triticum aestivum)與甜菜(Betavulgaris)的NRT1蛋白質(zhì)序列來(lái)自NCBI;藜麥(Chenopodium quinoa Willd)、玉米(Zea mays)、大豆(Glycine max)、谷子(Setaria italic)的NRT1蛋白質(zhì)序列來(lái)自Ensembl Plants數(shù)據(jù)庫(kù)(http://plants. ensembl. org/index.html)。

    1.2 生物信息學(xué)分析方法

    采用Gene Structure Display Sever數(shù)據(jù)庫(kù)(http://gsds.cbi.pku.edu.cn/)[14]分析基因的結(jié)構(gòu) ;采用 Ex-PASY-ProtParam tool數(shù)據(jù)庫(kù)(http://web.expasy.org/protparam)對(duì)蛋白質(zhì)的理化性質(zhì)進(jìn)行分析;采用Sig-nalP-5.0 Server數(shù)據(jù)庫(kù)(www.cbs.dtu.dk/services/SignalP/)預(yù)測(cè)蛋白質(zhì)是否含有信號(hào)肽;采用PSORT數(shù)據(jù)庫(kù)(psort1.hgc.jp/form.html)對(duì)蛋白質(zhì)進(jìn)行亞細(xì)胞定位分析預(yù)測(cè)。

    采用SOPMA數(shù)據(jù)庫(kù)(https://npsa-prabi.ibcp.fr/cgibin/npsa_automat.pl? page=npsa_sopma.html)對(duì)蛋白質(zhì)二級(jí)結(jié)構(gòu)進(jìn)行預(yù)測(cè),采用Phyre2數(shù)據(jù)庫(kù)(www.sbg.bio.ic.ac.uk/)進(jìn)行蛋白質(zhì)三級(jí)結(jié)構(gòu)模型的預(yù)測(cè),采用Ramachandran Plot Analysis數(shù)據(jù)庫(kù)(http://mordred.bioc.cam.ac.uk/)檢驗(yàn)蛋白質(zhì)三級(jí)結(jié)構(gòu)預(yù)測(cè)的準(zhǔn)確性。

    利用 MEGA 7[15]軟件對(duì)序列進(jìn)行 Clustal W[16]比對(duì),然后再使用鄰接法(neighbour joining,NJ)構(gòu)建系統(tǒng)進(jìn)化樹(shù)。利用已發(fā)表的SbNRT1基因表達(dá)Atalas數(shù)據(jù)[17-18],在phytozome(https://phytozome.jgi.doe.gov/pz/portal.html)中獲取目的基因在高粱各個(gè)生育期、不同器官部位的表達(dá)量,再利用Heml軟件作出表達(dá)量的熱圖。利用MEGA 7軟件對(duì)CDS序列進(jìn)行比對(duì),再利用TBtools軟件計(jì)算出基因的進(jìn)化選擇壓力。

    2 結(jié)果與分析

    2.1 SbNRT1基因家族基本信息

    利用Phytozome和NCBI數(shù)據(jù)庫(kù)進(jìn)行篩選,結(jié)合高粱各組織的轉(zhuǎn)錄組測(cè)序數(shù)據(jù),得到SbNRT1基因20個(gè)(表1)。SbNRT1基因長(zhǎng)度為1 762~12 364 bp,分布在第SBI-01、SBI-02、SBI-03、SBI-04、SBI-05、SBI-06、SBI-07、SBI-09、SBI-10共9條染色體上,每個(gè)SbNRT1基因都有一個(gè)或多個(gè)轉(zhuǎn)錄本。本研究選用的為初級(jí)轉(zhuǎn)錄本蛋白,具體蛋白質(zhì)的編號(hào)如表1所示。

    表1 SbNRT1基因家族基本信息Tab. 1 Basic information for the SbNRT1 gene family

    2.2 SbNRT1基因結(jié)構(gòu)分析

    SbNRT1基因結(jié)構(gòu)分析顯示:絕大多數(shù)SbNRT1均存在上下游非編碼區(qū),但其長(zhǎng)度均比較短;CDS為編碼蛋白產(chǎn)物的序列,大部分基因CDS數(shù)量很少,在3~4個(gè)之間,小部分基因CDS數(shù)量在5~8個(gè)之間 ;其中Sobic.001G133900、Sobic.001G276600、Sobic.007G044300、Sobic.010G175900等基因的內(nèi)含子較長(zhǎng),由于外顯子與內(nèi)含子接頭區(qū)存在一段高度保守的一致序列,這幾個(gè)基因相對(duì)其他SbNRT1存在更多的高度保守序列數(shù)量(圖1)。

    圖1 SbNRT1基因結(jié)構(gòu)Fig. 1 SbNRT1 gene structure

    2.3 SbNRT1蛋白質(zhì)理化性質(zhì)

    由表2可知:SbNRT1蛋白質(zhì)氨基酸數(shù)目在250~667個(gè)之間;分子量大小在26 334.30~72 337.95 Da之間;理論等電點(diǎn)在5.72~9.53之間,大部分SbNRT1蛋白質(zhì)等電點(diǎn)大于7;不穩(wěn)定系數(shù)分析表明,14個(gè)SbNRT1蛋白為穩(wěn)定蛋白,6個(gè)Sb-NRT1蛋白為不穩(wěn)定蛋白(不穩(wěn)定系數(shù)>40);脂溶系數(shù)的高低決定了蛋白質(zhì)的流動(dòng)性,系數(shù)越高,流動(dòng)性越好,除Sobic.009G049900蛋白外,其余19個(gè)SbNRT1蛋白的脂溶系數(shù)均大于90,表明SbNRT1蛋白的流動(dòng)性較好;親水性平均值都大于0,表明Sb-NRT1蛋白質(zhì)均為疏水性蛋白[19]。

    表2 SbNRT1蛋白質(zhì)理化性質(zhì)Tab. 2 The physical properties of SbNRT1 proteins

    2.4 SbNRT1蛋白質(zhì)信號(hào)肽及亞細(xì)胞定位分析

    根據(jù)蛋白質(zhì)信號(hào)肽分布,Sobic.004G260500、Sobic.006G251200、Sobic.007G044300、Sobic.007G081300、Sobic.009G049900、Sobic.010G099700這6個(gè)蛋白質(zhì)信號(hào)肽分值(signal peptide score,S)均>0.5,表明蛋白含信號(hào)肽,為分泌性蛋白;剩余的蛋白信號(hào)肽剪切位點(diǎn)分值(cut score,C)、綜合剪切位點(diǎn)分值(synthetical cut score,Y)均<0.5,表明蛋白都不含信號(hào)肽,不是分泌性蛋白。蛋白質(zhì)亞細(xì)胞定位預(yù)測(cè)顯示,有17個(gè)蛋白定位在質(zhì)膜上,概率區(qū)間為60%~82%,有3個(gè)蛋白定位在葉綠體類(lèi)囊體膜上,概率分別為94.8%、93.9%和66.4%(表3)。

    表3 SbNRT1蛋白質(zhì)亞細(xì)胞定位Tab. 3 Subcellular location of SbNRT1 proteins

    2.5 SbNRT1蛋白二級(jí)結(jié)構(gòu)分析

    SbNRT1蛋白家族的二級(jí)結(jié)構(gòu)主要由α-螺旋和無(wú)規(guī)則卷曲構(gòu)成(二者在蛋白質(zhì)中所占氨基酸數(shù)量的比例大于70%),而延伸鏈(β-折疊的組成結(jié)構(gòu))、β-轉(zhuǎn)角比例較低(表4)。由此可以推測(cè)出,α-螺旋和無(wú)規(guī)則卷曲是SbNRT1蛋白的大量結(jié)構(gòu)元件,而延伸鏈和β-轉(zhuǎn)角則散布在整個(gè)蛋白質(zhì)中。

    表4 SbNRT1蛋白二級(jí)結(jié)構(gòu)分析Tab. 4 Analysis of the secondary structure of SbNRT1 protein unit: %

    2.6 SbNRT1蛋白三級(jí)結(jié)構(gòu)模型

    采用phyre 2軟件預(yù)測(cè)高粱NRT1蛋白的三級(jí)結(jié)構(gòu),結(jié)果顯示,NRT1蛋白三級(jí)結(jié)構(gòu)主要由α-螺旋和無(wú)規(guī)則卷曲構(gòu)成,與二級(jí)結(jié)構(gòu)預(yù)測(cè)結(jié)果一致(圖2)。采用Ramachandran plot analysis軟件評(píng)估建模模型,結(jié)果顯示,氨基酸殘基所處的最佳區(qū)域比例為63.00%~88.50%,可接受區(qū)域氨基酸殘基比率為8.90%~32.80%,SbNRT1蛋白質(zhì)三級(jí)結(jié)構(gòu)模型的可信度均大于93.00%(表5)。

    表5 SbNRT1蛋白質(zhì)三級(jí)結(jié)構(gòu)模型評(píng)估Tab. 5 Evaluation of tertiary structure model of SbNRT1 protein unit: %

    圖2 SbNRT1蛋白質(zhì)三級(jí)結(jié)構(gòu)模型Fig. 2 Tertiary structure models of SbNRT1 proteins

    2.7 SbNRT1基因系統(tǒng)進(jìn)化樹(shù)分析

    SbNRT1基因的系統(tǒng)進(jìn)化樹(shù)分析表明,Sb-NRT1基因進(jìn)化樹(shù)的變異值為0.1,高粱各基因間相互存在同源性,但大部分基因同源性不高。其中Sobic.001G302800與Sobic.004G193000,Sobic.003G185100與Sobic.007G044300,Sobic.001G133900與Sobic.009G148900等旁系同源序列的同源性較高(圖3)。對(duì)高粱、擬南芥、大豆、玉米和谷子的NRT1基因家族進(jìn)行聚類(lèi)分析(圖4),不同顏色范圍代表不同亞族,結(jié)果表明:NRT1基因家族分為8類(lèi),即8個(gè)亞族,Sobic.004G260500與Zm00001eb189000,Sobic.002G091800與Zm00001eb304390,Sobic.010G099700與Zm00001eb371840、KQL10112,Sobic.003G399200與KQL15490、Zm00001eb025860,Sobic.003G107100與Zm00001eb126980互為直系同源基因,說(shuō)明高粱與谷子、玉米的NRT1基因在進(jìn)化過(guò)程中同源關(guān)系更近;進(jìn)化樹(shù)中的雙子葉植物(大豆、擬南芥)與單子葉植物(高粱、玉米和谷子)大多不能聚集在同一分支下,說(shuō)明NRT1基因在進(jìn)化過(guò)程中出現(xiàn)單雙子葉的分化。

    圖3 SbNRT1基因家族聚類(lèi)圖Fig. 3 Clustering of SbNRT1 gene family

    圖4 高粱、玉米、擬南芥、大豆、谷子的NRT1基因家族的系統(tǒng)進(jìn)化樹(shù)Fig. 4 Phylogenetic tree of NRT1 gene family in Sorghum bicolor, Zea mays L., Arabidopsis, Glycine max (Linn.) Merr, Setaria italica

    2.8 SbNRT1基因表達(dá)

    利用高粱表達(dá)Atalas數(shù)據(jù)構(gòu)建表達(dá)量熱圖(圖5),可知SbNRT1基因的表達(dá)具有時(shí)空專(zhuān)一性。Sobic.010G175900、Sobic.010G099700、Sobic.009G148900、Sobic.009G049900、Sobic.007G160900等基因在整個(gè)高粱生育期表達(dá)很低;Sobic.003G295500在整個(gè)生育期表達(dá)水平較低 ;但是Sobic.001G133900、Sobic.004G193000、Sobic.001G302800在生育期某個(gè)階段表達(dá)水平比較高,其中,Sobic.004G193000在生育期多個(gè)部位表達(dá)水平較高,尤其在左下生長(zhǎng)花期表達(dá)水平最高,Sobic.001G302800在左葉、葉鞘生長(zhǎng)花部位表達(dá)水平較高。從SbNRT1基因整個(gè)表達(dá)時(shí)期可以看出,表達(dá)水平由高到低依次為左下生長(zhǎng)花期、葉片、莖中部節(jié)間、根部。

    圖5 SbNRT1基因家族表達(dá)熱圖Fig. 5 Heat map for expression of SbNRT1 gene family

    在高粱根部Sobic.001G133900的表達(dá)受銨鹽和尿素誘導(dǎo),受硝酸鹽誘導(dǎo)表達(dá)最強(qiáng)的為Sobic.001G133900、Sobic.003G295500和Sobic.001G302800,而Sobic.003G295500似乎僅在硝酸鹽誘導(dǎo)下表達(dá)。

    2.9 SbNRT1進(jìn)化選擇分析

    在遺傳學(xué)中用非同義突變率(non-synonymous mutation rate,ka)與同義突變率(synonymous mutation rate,ks)的比值作為這個(gè)蛋白質(zhì)編碼基因是否有選擇壓力作用的標(biāo)準(zhǔn)。如果ka/ks遠(yuǎn)大于1,則存在正選擇效應(yīng);如果ka/ks=1,則認(rèn)為存在中性選擇;如果ka/ks遠(yuǎn)小于1,則存在純化選擇作用。通過(guò)對(duì)Sb-NRT1同源蛋白基因?qū)Φ倪M(jìn)化選擇壓力分析,得出SbNRT1蛋白編碼區(qū)基因的ka/ks遠(yuǎn)遠(yuǎn)小于1,說(shuō)明SbNRT1蛋白質(zhì)編碼基因存在純化選擇作用(表6)。

    表6 SbNRT1進(jìn)化選擇參數(shù)Tab. 6 SbNRT1 evolutionary selection parameters

    3 討論

    高粱抗旱性和耐鹽性強(qiáng),在貧瘠土壤上種植也能有一定的收成,是典型的“環(huán)境友好型”作物,是研究作物抗逆性的模式植物之一。由于高粱抗逆境能力比其他作物更強(qiáng),所以越來(lái)越多的植物科學(xué)家將研究的重點(diǎn)轉(zhuǎn)移到了高粱上,而其中轉(zhuǎn)運(yùn)蛋白基因的挖掘更是成為一個(gè)新的研究熱點(diǎn)。目前Sb-NRT1基因的鑒定和系統(tǒng)分析尚未見(jiàn)報(bào)道。本研究從全基因組水平上對(duì)SbNRT1基因家族進(jìn)行鑒定,并進(jìn)行了表達(dá)分析和DNA選擇分析,為研究其功能、建立分子標(biāo)記輔助育種奠定了基礎(chǔ)。

    本研究結(jié)果表明:高粱全基因組內(nèi)鑒定出SbNRT1基因共有20個(gè),分布在SBI-01、SBI-02、SBI-03、SBI-04、SBI-05、SBI-06、SBI-07、SBI-09、SBI-10共9條染色體上;SbNRT1蛋白質(zhì)均為疏水性蛋白,且理論等電點(diǎn)大于7;6個(gè)SbNRT1蛋白為不穩(wěn)定蛋白,14個(gè)SbNRT1蛋白為穩(wěn)定蛋白;75%以上的SbNRT1蛋白質(zhì)不存在信號(hào)肽,不是分泌性蛋白,這與大多數(shù)SbNRT1蛋白亞細(xì)胞定位在質(zhì)膜上的結(jié)果表現(xiàn)一致;SbNRT1蛋白質(zhì)的二級(jí)結(jié)構(gòu)主要是α-螺旋和無(wú)規(guī)則卷曲;SbNRT1蛋白質(zhì)的三級(jí)結(jié)構(gòu)比較復(fù)雜,但模型可信度均大于93%。

    系統(tǒng)進(jìn)化樹(shù)分析結(jié)果表明,SbNRT1蛋白家族分為8個(gè)亞家族,高粱與玉米、谷子的NRT1蛋白質(zhì)在進(jìn)化上同源關(guān)系更近。此結(jié)果與葉玲等[20]對(duì)谷子NRT2基因家族生物信息分析結(jié)果一致。從表達(dá)熱圖上可以分析得出,SbNRT1基因表達(dá)具有較強(qiáng)的器官、生育期特異性,尤其以Sobic.004G193000基因表達(dá)水平較高。在高粱根部受硝酸鹽誘導(dǎo)表達(dá)最強(qiáng)的基因?yàn)镾obic.001G133900、Sobic.003G295500和Sobic.001G302800,前者也受銨鹽和尿素誘導(dǎo),而Sobic.003G295500似乎僅在硝酸鹽誘導(dǎo)下表達(dá)。以上結(jié)果為研究基因功能及氮肥響應(yīng)提供了分子基礎(chǔ)。

    對(duì)SbNRT1基因的表達(dá)進(jìn)行分析有助于進(jìn)一步鑒定轉(zhuǎn)運(yùn)蛋白相關(guān)基因,為高粱相關(guān)基因的進(jìn)一步鑒定與功能研究奠定理論基礎(chǔ),為后續(xù)分子育種提供思路,也為揭示該基因家族參與植物氮素吸收利用的調(diào)控機(jī)制提供了理論依據(jù)。

    猜你喜歡
    信號(hào)肽進(jìn)化樹(shù)硝酸鹽
    嵌合信號(hào)肽提高α-淀粉酶在枯草芽孢桿菌中的分泌
    硝酸鹽并不致癌還或有益處
    中老年保健(2022年3期)2022-11-21 09:40:36
    基于心理旋轉(zhuǎn)的小學(xué)生物進(jìn)化樹(shù)教學(xué)實(shí)驗(yàn)報(bào)告
    常見(jiàn)的進(jìn)化樹(shù)錯(cuò)誤概念及其辨析*
    運(yùn)用計(jì)算機(jī)軟件預(yù)測(cè)木質(zhì)部寄生屬信號(hào)肽
    艾草白粉病的病原菌鑒定
    家畜硝酸鹽和亞硝酸鹽中毒的診斷、鑒別和防治
    內(nèi)源信號(hào)肽DSE4介導(dǎo)頭孢菌素C?;冈诋叧嘟湍钢械姆置诒磉_(dá)
    短期水分脅迫影響巴旦杏植株對(duì)硝酸鹽的吸收
    層次聚類(lèi)在進(jìn)化樹(shù)構(gòu)建中的應(yīng)用
    国产一区二区亚洲精品在线观看| 亚洲最大成人av| 各种免费的搞黄视频| 日韩三级伦理在线观看| 一级a做视频免费观看| 久久久国产一区二区| 大片电影免费在线观看免费| 国产欧美另类精品又又久久亚洲欧美| 国国产精品蜜臀av免费| 国产黄片视频在线免费观看| 卡戴珊不雅视频在线播放| 精品国产三级普通话版| 久久久精品94久久精品| 国产精品久久久久久av不卡| 日韩 亚洲 欧美在线| 美女高潮的动态| 免费黄网站久久成人精品| 97在线人人人人妻| 国产午夜精品一二区理论片| 国产探花极品一区二区| 亚洲欧美精品自产自拍| 久久国产乱子免费精品| 免费av毛片视频| 亚洲国产欧美人成| 久久韩国三级中文字幕| 亚洲最大成人手机在线| 一级av片app| 久久精品国产亚洲av天美| 91精品国产九色| 精品人妻一区二区三区麻豆| 国产成人freesex在线| 国产av码专区亚洲av| 精品人妻偷拍中文字幕| 99久久精品热视频| 亚洲自偷自拍三级| 久久精品久久精品一区二区三区| 2021少妇久久久久久久久久久| 水蜜桃什么品种好| 日本三级黄在线观看| 国产色婷婷99| 国产久久久一区二区三区| 黄色欧美视频在线观看| 成年版毛片免费区| 大话2 男鬼变身卡| 3wmmmm亚洲av在线观看| av国产免费在线观看| 少妇人妻 视频| 高清日韩中文字幕在线| 久久久久久久精品精品| 久久久久网色| 韩国高清视频一区二区三区| 国产综合懂色| 青春草国产在线视频| 欧美成人精品欧美一级黄| 人妻系列 视频| 直男gayav资源| 成人综合一区亚洲| 天天一区二区日本电影三级| 亚洲精品自拍成人| 99久久中文字幕三级久久日本| 97在线视频观看| av在线播放精品| 一区二区三区四区激情视频| 在线观看免费高清a一片| 色婷婷久久久亚洲欧美| 香蕉精品网在线| 免费黄频网站在线观看国产| 日韩国内少妇激情av| av.在线天堂| 久久久午夜欧美精品| 视频区图区小说| 国产av不卡久久| 免费观看性生交大片5| 亚洲欧美清纯卡通| av卡一久久| 亚洲欧美一区二区三区国产| 国产 精品1| 99热这里只有是精品在线观看| 色视频www国产| 免费黄频网站在线观看国产| 人人妻人人澡人人爽人人夜夜| 免费高清在线观看视频在线观看| 中国三级夫妇交换| 亚洲高清免费不卡视频| 国产成人aa在线观看| 免费黄网站久久成人精品| 欧美日本视频| 永久免费av网站大全| 日韩中字成人| 日本猛色少妇xxxxx猛交久久| 久久精品久久精品一区二区三区| 久久精品夜色国产| 男女无遮挡免费网站观看| 岛国毛片在线播放| 欧美bdsm另类| 日本av手机在线免费观看| 日韩亚洲欧美综合| 国产一级毛片在线| 亚洲人成网站在线观看播放| 99视频精品全部免费 在线| 日韩一本色道免费dvd| 日韩不卡一区二区三区视频在线| 午夜激情久久久久久久| 亚洲精品影视一区二区三区av| 熟女人妻精品中文字幕| 丝袜喷水一区| 免费看光身美女| 国产熟女欧美一区二区| 亚洲欧美精品专区久久| 亚洲欧洲日产国产| 国产亚洲av嫩草精品影院| 性色av一级| 高清在线视频一区二区三区| 国产欧美亚洲国产| 亚洲精品一二三| 精品熟女少妇av免费看| 日产精品乱码卡一卡2卡三| 久久久成人免费电影| 性色avwww在线观看| 蜜桃久久精品国产亚洲av| 久久久久久久久久久免费av| 国产黄a三级三级三级人| 国产人妻一区二区三区在| 人妻 亚洲 视频| 亚洲精品日韩在线中文字幕| 丝袜美腿在线中文| 亚洲天堂国产精品一区在线| 男女啪啪激烈高潮av片| 精品少妇久久久久久888优播| 99热全是精品| 午夜免费男女啪啪视频观看| 国内精品美女久久久久久| 日本黄大片高清| 免费播放大片免费观看视频在线观看| 日日撸夜夜添| 六月丁香七月| 久久精品久久久久久久性| 热re99久久精品国产66热6| 大又大粗又爽又黄少妇毛片口| 99久久九九国产精品国产免费| av网站免费在线观看视频| 亚洲av成人精品一二三区| 亚洲最大成人av| 嫩草影院入口| 成人鲁丝片一二三区免费| 又黄又爽又刺激的免费视频.| 少妇高潮的动态图| 亚洲天堂av无毛| 男女边吃奶边做爰视频| 亚洲欧美清纯卡通| 久久久久久久国产电影| 亚洲国产日韩一区二区| 91久久精品电影网| 久久久久久九九精品二区国产| 久久精品久久久久久噜噜老黄| 777米奇影视久久| 国产av国产精品国产| 久久久久网色| av黄色大香蕉| 国产毛片在线视频| 亚洲欧美日韩东京热| 成人毛片60女人毛片免费| 搡女人真爽免费视频火全软件| 听说在线观看完整版免费高清| 日本欧美国产在线视频| 在线观看一区二区三区激情| 天堂中文最新版在线下载 | 久久精品国产亚洲网站| 日本爱情动作片www.在线观看| 中国国产av一级| 国内精品宾馆在线| 日本欧美国产在线视频| 春色校园在线视频观看| 国产大屁股一区二区在线视频| kizo精华| 男女啪啪激烈高潮av片| 亚洲内射少妇av| 好男人视频免费观看在线| 男人狂女人下面高潮的视频| 麻豆成人av视频| 美女脱内裤让男人舔精品视频| 成人国产麻豆网| 国产高清有码在线观看视频| 狠狠精品人妻久久久久久综合| 91精品一卡2卡3卡4卡| 国产免费福利视频在线观看| 99热网站在线观看| 国产老妇伦熟女老妇高清| 国产成人精品一,二区| 亚洲经典国产精华液单| 免费播放大片免费观看视频在线观看| 在线观看免费高清a一片| 亚洲成人一二三区av| 成人午夜精彩视频在线观看| 男女啪啪激烈高潮av片| 欧美激情久久久久久爽电影| 亚洲成人中文字幕在线播放| 久久97久久精品| 人妻系列 视频| 国产精品人妻久久久久久| 纵有疾风起免费观看全集完整版| 男女边吃奶边做爰视频| 麻豆精品久久久久久蜜桃| av黄色大香蕉| 国产免费福利视频在线观看| 综合色av麻豆| 久久精品国产亚洲网站| www.av在线官网国产| 日日啪夜夜撸| 国产亚洲午夜精品一区二区久久 | 人妻制服诱惑在线中文字幕| 黄片wwwwww| 18禁在线播放成人免费| 久热久热在线精品观看| 黄色日韩在线| 久久久成人免费电影| 一本久久精品| 免费少妇av软件| 纵有疾风起免费观看全集完整版| 日韩av免费高清视频| 亚洲内射少妇av| 中文乱码字字幕精品一区二区三区| 大又大粗又爽又黄少妇毛片口| 最近中文字幕2019免费版| 久久久久久久久大av| 一本久久精品| 国产成人精品一,二区| 午夜福利网站1000一区二区三区| 午夜激情久久久久久久| 性色av一级| 久久精品国产a三级三级三级| 亚洲欧美日韩东京热| 大陆偷拍与自拍| 国产精品一及| 下体分泌物呈黄色| 一级爰片在线观看| 欧美xxⅹ黑人| 久久久精品欧美日韩精品| 日韩 亚洲 欧美在线| 熟女av电影| 久久国产乱子免费精品| 久久精品久久精品一区二区三区| 免费看日本二区| 内射极品少妇av片p| 黄色日韩在线| 日韩不卡一区二区三区视频在线| 18禁在线无遮挡免费观看视频| 97在线人人人人妻| 高清午夜精品一区二区三区| 嫩草影院精品99| 亚洲精品一区蜜桃| 国产精品国产三级专区第一集| 嫩草影院入口| 91精品一卡2卡3卡4卡| 99久久中文字幕三级久久日本| 色视频www国产| 欧美三级亚洲精品| 欧美zozozo另类| 国产淫语在线视频| 亚洲国产精品专区欧美| 人妻系列 视频| 日本免费在线观看一区| 色视频在线一区二区三区| 可以在线观看毛片的网站| 又大又黄又爽视频免费| 又爽又黄无遮挡网站| 国产午夜福利久久久久久| 观看美女的网站| 国产成人精品一,二区| 99热国产这里只有精品6| 最近2019中文字幕mv第一页| 国内精品美女久久久久久| 男男h啪啪无遮挡| 尾随美女入室| av又黄又爽大尺度在线免费看| 嘟嘟电影网在线观看| 两个人的视频大全免费| 国产免费一区二区三区四区乱码| 国产成人一区二区在线| 激情 狠狠 欧美| 亚洲人成网站高清观看| 精品久久久久久久久亚洲| 久久国内精品自在自线图片| 国产伦在线观看视频一区| 一本—道久久a久久精品蜜桃钙片 精品乱码久久久久久99久播 | 国产成人freesex在线| 亚洲国产欧美在线一区| 国产乱人偷精品视频| 亚洲精品成人久久久久久| 亚洲成人久久爱视频| 美女脱内裤让男人舔精品视频| 国产高清有码在线观看视频| 日本熟妇午夜| 亚洲av一区综合| 另类亚洲欧美激情| 精品一区在线观看国产| 欧美变态另类bdsm刘玥| 一级毛片电影观看| 午夜精品一区二区三区免费看| 亚洲伊人久久精品综合| 日本猛色少妇xxxxx猛交久久| 久久精品国产a三级三级三级| 欧美日韩在线观看h| 男插女下体视频免费在线播放| 免费大片黄手机在线观看| 综合色av麻豆| 少妇的逼水好多| 大码成人一级视频| 精品久久久噜噜| 久久久久久伊人网av| 国产欧美日韩一区二区三区在线 | 国产91av在线免费观看| 一区二区三区精品91| 水蜜桃什么品种好| 在线观看av片永久免费下载| 亚洲国产色片| 欧美精品人与动牲交sv欧美| 久久久久网色| 国产免费一级a男人的天堂| 成人无遮挡网站| 人妻制服诱惑在线中文字幕| 免费观看性生交大片5| 色5月婷婷丁香| 久久久久久国产a免费观看| 天堂网av新在线| 国产成人91sexporn| 99热这里只有精品一区| 久久久久精品久久久久真实原创| 免费黄网站久久成人精品| 日本午夜av视频| 韩国av在线不卡| 男女啪啪激烈高潮av片| 少妇人妻精品综合一区二区| 久久久亚洲精品成人影院| 99久国产av精品国产电影| 九九在线视频观看精品| 高清欧美精品videossex| 亚洲av二区三区四区| 日日摸夜夜添夜夜爱| 日本熟妇午夜| 午夜免费鲁丝| 亚洲最大成人中文| 国产高清有码在线观看视频| 99久久九九国产精品国产免费| 日本爱情动作片www.在线观看| 伦精品一区二区三区| 国产精品爽爽va在线观看网站| 人妻 亚洲 视频| 国产男女内射视频| av专区在线播放| 亚洲内射少妇av| 国产淫语在线视频| 久久久精品94久久精品| 亚洲欧美成人精品一区二区| 青青草视频在线视频观看| 中国美白少妇内射xxxbb| 一级毛片电影观看| 久久亚洲国产成人精品v| 搡老乐熟女国产| 波多野结衣巨乳人妻| 日韩av免费高清视频| 国产探花在线观看一区二区| 国产黄a三级三级三级人| 成人一区二区视频在线观看| 精品久久国产蜜桃| 人妻 亚洲 视频| 97在线视频观看| 80岁老熟妇乱子伦牲交| 一级毛片电影观看| 国产毛片在线视频| 国产高清有码在线观看视频| 国产黄片美女视频| 久久午夜福利片| 亚洲经典国产精华液单| 日韩免费高清中文字幕av| 国产有黄有色有爽视频| 国产亚洲精品久久久com| 免费观看av网站的网址| 亚洲av免费高清在线观看| 制服丝袜香蕉在线| 国产精品99久久99久久久不卡 | 国产久久久一区二区三区| 草草在线视频免费看| 成年av动漫网址| 日韩不卡一区二区三区视频在线| 精品人妻一区二区三区麻豆| 五月玫瑰六月丁香| 久热这里只有精品99| 少妇的逼好多水| 日韩电影二区| 久久久精品欧美日韩精品| 日日摸夜夜添夜夜添av毛片| 欧美日韩视频高清一区二区三区二| 亚洲无线观看免费| 日日摸夜夜添夜夜爱| 三级经典国产精品| 最后的刺客免费高清国语| 亚洲国产高清在线一区二区三| 亚洲综合精品二区| 永久网站在线| 97精品久久久久久久久久精品| freevideosex欧美| 国产v大片淫在线免费观看| 精品国产乱码久久久久久小说| 欧美高清性xxxxhd video| 搡女人真爽免费视频火全软件| 国产乱人偷精品视频| 18禁裸乳无遮挡动漫免费视频 | 18禁在线播放成人免费| 国产精品一区二区三区四区免费观看| 成人无遮挡网站| 亚洲第一区二区三区不卡| 国产精品熟女久久久久浪| 少妇高潮的动态图| 春色校园在线视频观看| 日韩人妻高清精品专区| 一级av片app| 成人国产麻豆网| 中文在线观看免费www的网站| xxx大片免费视频| 日韩一区二区视频免费看| 国产精品国产三级国产av玫瑰| 国产高清三级在线| 精品久久久久久久久av| 另类亚洲欧美激情| av在线播放精品| av福利片在线观看| 亚洲四区av| 国产高潮美女av| 久久精品夜色国产| 精品久久久久久久久亚洲| 亚洲精品一区蜜桃| 日韩av免费高清视频| 免费av观看视频| 久久99热这里只频精品6学生| www.av在线官网国产| 欧美成人一区二区免费高清观看| 免费观看在线日韩| 亚洲欧美一区二区三区国产| 国产 一区精品| 久久精品国产亚洲av天美| 97超视频在线观看视频| 国产又色又爽无遮挡免| 黄色怎么调成土黄色| 大香蕉久久网| 亚洲欧美清纯卡通| 亚洲精品日韩在线中文字幕| 中文字幕制服av| 99久久中文字幕三级久久日本| 成人一区二区视频在线观看| 一级毛片久久久久久久久女| 夫妻性生交免费视频一级片| 国产精品嫩草影院av在线观看| 精品人妻偷拍中文字幕| 如何舔出高潮| 99热网站在线观看| 80岁老熟妇乱子伦牲交| 亚洲色图综合在线观看| 肉色欧美久久久久久久蜜桃 | 少妇丰满av| 丝袜喷水一区| 激情 狠狠 欧美| 黄色一级大片看看| 久久久色成人| 国产精品福利在线免费观看| 蜜桃久久精品国产亚洲av| 爱豆传媒免费全集在线观看| 亚洲精品一二三| 亚洲aⅴ乱码一区二区在线播放| 在线天堂最新版资源| 国产精品久久久久久精品电影小说 | 永久免费av网站大全| tube8黄色片| 亚州av有码| 国产成人午夜福利电影在线观看| 免费观看无遮挡的男女| 身体一侧抽搐| 国产精品久久久久久精品古装| 久久久成人免费电影| 日韩电影二区| 最近中文字幕高清免费大全6| 亚洲成色77777| 亚洲国产精品成人久久小说| 街头女战士在线观看网站| 狂野欧美白嫩少妇大欣赏| 97人妻精品一区二区三区麻豆| videos熟女内射| 亚洲内射少妇av| 久久亚洲国产成人精品v| 欧美激情在线99| 九九在线视频观看精品| 一级黄片播放器| eeuss影院久久| 最近最新中文字幕大全电影3| 国产亚洲午夜精品一区二区久久 | 午夜福利视频1000在线观看| 亚洲真实伦在线观看| 精品99又大又爽又粗少妇毛片| 各种免费的搞黄视频| 日韩在线高清观看一区二区三区| 91在线精品国自产拍蜜月| 肉色欧美久久久久久久蜜桃 | 亚洲伊人久久精品综合| 赤兔流量卡办理| av国产精品久久久久影院| 天堂中文最新版在线下载 | av网站免费在线观看视频| 蜜桃亚洲精品一区二区三区| 最近最新中文字幕免费大全7| 免费看不卡的av| av卡一久久| 在线观看av片永久免费下载| 亚洲久久久久久中文字幕| 2018国产大陆天天弄谢| 精品少妇黑人巨大在线播放| 亚洲最大成人手机在线| 深夜a级毛片| 涩涩av久久男人的天堂| 国产亚洲最大av| 人人妻人人澡人人爽人人夜夜| 熟女人妻精品中文字幕| 国产精品国产三级专区第一集| 成年免费大片在线观看| 黄片无遮挡物在线观看| 在线观看免费高清a一片| 国产爽快片一区二区三区| 国产午夜福利久久久久久| 五月天丁香电影| 亚洲av国产av综合av卡| 亚洲欧美日韩另类电影网站 | 天堂网av新在线| 丰满人妻一区二区三区视频av| 性色avwww在线观看| 日本-黄色视频高清免费观看| 亚洲欧美清纯卡通| 永久网站在线| 一区二区三区免费毛片| .国产精品久久| 午夜亚洲福利在线播放| 日韩国内少妇激情av| 综合色丁香网| 久久精品国产亚洲av涩爱| 啦啦啦中文免费视频观看日本| 人人妻人人看人人澡| 爱豆传媒免费全集在线观看| 国产成人a区在线观看| 亚洲自偷自拍三级| 在线免费观看不下载黄p国产| 成人亚洲精品一区在线观看 | 久久国产乱子免费精品| 如何舔出高潮| 最新中文字幕久久久久| 免费黄网站久久成人精品| 国产精品av视频在线免费观看| 亚洲精品乱码久久久v下载方式| 日本wwww免费看| 精品亚洲乱码少妇综合久久| av免费在线看不卡| 男女无遮挡免费网站观看| 综合色av麻豆| 久久久久久久久久成人| 国精品久久久久久国模美| 亚洲国产日韩一区二区| 国精品久久久久久国模美| 人人妻人人爽人人添夜夜欢视频 | 大香蕉97超碰在线| 一个人看的www免费观看视频| 天天躁日日操中文字幕| h日本视频在线播放| 18禁裸乳无遮挡免费网站照片| av天堂中文字幕网| 午夜福利视频1000在线观看| 日韩欧美 国产精品| 99热6这里只有精品| 久久国产乱子免费精品| 内地一区二区视频在线| 高清日韩中文字幕在线| 黄色欧美视频在线观看| 狂野欧美激情性xxxx在线观看| 麻豆成人av视频| 一个人观看的视频www高清免费观看| 九九在线视频观看精品| 卡戴珊不雅视频在线播放| 亚洲成人一二三区av| 日韩伦理黄色片| 亚洲精品一区蜜桃| 日本猛色少妇xxxxx猛交久久| 少妇的逼水好多| 啦啦啦中文免费视频观看日本| 三级经典国产精品| 国产精品久久久久久精品电影| 麻豆成人午夜福利视频| 亚洲精品亚洲一区二区| 国产亚洲一区二区精品| 日韩中字成人| 亚洲av在线观看美女高潮| 亚洲精品色激情综合| 亚洲人成网站高清观看| 国产国拍精品亚洲av在线观看| 联通29元200g的流量卡| 久久久久久久午夜电影| 欧美+日韩+精品| 久久97久久精品| 国产精品一区二区在线观看99| 亚洲欧美清纯卡通| 纵有疾风起免费观看全集完整版| 日韩 亚洲 欧美在线| 亚洲自拍偷在线| 别揉我奶头 嗯啊视频| 777米奇影视久久| 久久亚洲国产成人精品v| 91精品伊人久久大香线蕉| 在线看a的网站| 欧美3d第一页| 亚洲av男天堂|