999精品在线视频,手机成人午夜在线视频,久久不卡国产精品无码,中日无码在线观看,成人av手机在线观看,日韩精品亚洲一区中文字幕,亚洲av无码人妻,四虎国产在线观看 ?

斜帶石斑魚IL-10基因的克隆及其原核表達

2015-10-25 07:02:34肖周婷黃郁蔥簡紀常魯義善吳灶和
生物技術通報 2015年9期

肖周婷黃郁蔥簡紀常魯義善吳灶和

(1.廣東海洋大學水產學院,湛江 524088;2.廣東省水產經濟動物病原生物學及流行病學重點實驗室,湛江524088;3.廣東省教育廳水產經濟動物病害控制重點實驗室,湛江 524088;4.仲愷農業(yè)工程學院,廣州 510225)

斜帶石斑魚IL-10基因的克隆及其原核表達

肖周婷1,2,3黃郁蔥1,2,3簡紀常1,2,3魯義善1,2,3吳灶和2,3,4

(1.廣東海洋大學水產學院,湛江524088;2.廣東省水產經濟動物病原生物學及流行病學重點實驗室,湛江524088;3.廣東省教育廳水產經濟動物病害控制重點實驗室,湛江524088;4.仲愷農業(yè)工程學院,廣州510225)

白介素10(Interleukin-10,IL-10)是一個重要的多效細胞因子,主要作用是介導炎癥反應和調控一些免疫細胞的分化和增殖。利用同源克隆和RACE-PCR技術獲得全長為1 104 bp的斜帶石斑魚(Epinephelus coioides)IL-10基因。該基因ORF為564 bp,編碼187個氨基酸,預測蛋白分子量為21.7 kD,等電點為5.74。氨基酸同源性比對及進化分析結果顯示,斜帶石斑魚IL-10氨基酸序列和吉富羅非魚(O. niloticus)同源性最高,達72.25%。利用雙酶切及質粒重組技術構建IL-10原核表達載體,并轉化大腸桿菌BL21(DE3)中誘導原核表達。SDS-PAGE分析顯示,融合蛋白分子量為37.5 kD,與預期相符;在IPTG為0.02 mmol/L,37℃誘導3 h,蛋白包涵體的表達量最大。

斜帶石斑魚(Epinephelus coioides);Interleukin-10;基因克隆;重組表達

1989年,IL-10首次在鼠類T細胞中作為細胞因子合成抑制因子被報道[1]。除T細胞外,B細胞、嗜酸性粒細胞、單核巨噬細胞、NK細胞等,幾乎所有淋巴細胞均能合成IL-10[2]。IL-10主要調控炎癥反應,它不僅能通過抑制單核巨噬細胞釋放免疫介質、抗原呈遞和細胞吞噬作用,抑制CD4+T細胞的增殖和細胞因子合成,抑制呼吸爆發(fā)等,從而抑制免疫作用;也能通過刺激NK細胞、CD8+T細胞的細胞毒活性增強免疫作用[3-6]。在哺乳動物研究中表明,IL-10是以二聚體形式行使其生物學功能,通過與IL-10R1/R2受體配體復合物結合,使JAKSTAT信號通路上JAK1和Tyk2的磷酸化,從而招募和激活轉錄因子STAT3等,選擇性激活或抑制基因轉錄,發(fā)揮其抗炎生物學功能[7-10]。除JAK-STAT信號通路,IL-10還通過其他信號通路介導免疫反應。

目前為止,IL-10已在河豚、鯉、斑馬魚、黑鱸、鱈魚、金魚、虹鱒等硬骨魚中被發(fā)現(xiàn)[11-17],其中虹鱒中發(fā)現(xiàn)兩個同源性92%的IL-10剪切體。但對于斜帶石斑魚IL-10相關的研究還未見相關報道。本實驗對斜帶石斑魚IL-10基因進行克隆及其生物信息學分析,以期為進一步研究該基因的結構特點和免疫調節(jié)作用奠定基礎;同時構建IL-10原核表達載體,表達IL-10融合蛋白,旨在為進一步研究斜帶石斑魚IL-10蛋白生物學功能提供參考。

1 材料與方法

1.1 材料

1.1.1 實驗用魚 實驗用斜帶石斑魚采集自湛江某水產品市場。

1.1.2 載體與菌株 克隆載體pMD18-T 質粒和表達載體 pET-32a(+)質粒均購自TaKaRa公司,大腸桿菌 DH5α和BL21菌株來自本實驗室保存的菌株。

1.1.3 主要試劑 UNIQ-10柱式Trizol總RNA提取試劑盒購自上海生工生物技術有限公司;質粒提取試劑盒及DNA膠回收試劑盒購自Thermo公司;PrimeSTAR HS DNA 聚合酶、Ex Taq DNA聚合酶、限制性核酸內切酶Xho I和Bam H I、T4 DNA 連接酶、末端轉移酶TdT、dNTPs、DNA Marker、蛋白Marker等購自TaKaRa公司;M-MLV反轉錄酶購自NEB公司;用PCR的引物及序列測定均由上海生工生物技術有限公司完成。

1.2 方法

1.2.1 斜帶石斑魚總RNA的提取及cDNA第一鏈合成 斜帶石斑魚脾臟總RNA提取按照UNIQ-10柱式Trizol總RNA提取試劑盒說明書進行;按照M-MLV反轉錄酶說明書以3adapter為引物合成第一鏈cDNA,進行同源克隆及3'-RACE;以5-RECE-D1為引物合成cDNA,按照DNA膠回收試劑盒說明書對cDNA進行純化,純化后的cDNA 按照末端轉移酶TdT說明書對其加上ployA 尾,產物用于進行5'-RACE。

1.2.2 斜帶石斑魚IL-10中間片段獲得 根據(jù)GenBank上已登錄的舌齒鱸Dicentrarchus labrax(DQ821114.1)和鯽Carassius auratus(HQ259106.1)設計簡并引物IL-10-F1和IL-10-D1,同源克隆獲得IL-10的基因片段。

1.2.3 斜帶石斑魚IL-10基因3'-RACE 根據(jù)獲得的中間片段設計3'-RACE特異引物(3-RACE-F1;3-RACE-F2)。3'-RACE反應體系:第一輪PCR:cDNA 2.5 μL,3-RACE-F1(10 μmol/L)1 μL,long(0.4 μmol/L)/short(2 μmol/L)2 μL,PrimeSTAR HS DNA 聚合酶0.25 μL,5×PrimeSTAR Buffer 5 μL,dNTP mix(2.5 mmol/L)2.5 μL,加水補足25 μL。PCR反應條件:94℃ 30 s,66℃ 45 s,72℃ 2 min,5個循環(huán);94℃ 30 s,60℃ 45 s,72℃ 2 min,25個循環(huán),共計30個循環(huán)。第二輪PCR:首輪PCR產物稀釋20倍1 μL,3anchor 1 μL,3-RACE-F2(10 μmol/L)1 μL,Ex Taq DNA 聚合酶0.25 μL,10×Ex Taq Buffer 2.5 μL,dNTP mix(2.5 mmol/L)2.5 μL,加水補足25 μL。反應條件同第一輪PCR。

1.2.4 斜帶石斑魚IL-10基因5'-RACE 根據(jù)獲得的中間片段設計5'-RACE特異引物(5-RECE-D1;5-RACE-D2;5-RACE-D3)。5'-RACE反 應 體 系:第一輪PCR:加尾產物 2.5 μL,5-RACE-D2(10 μmol/L)1 μL,5oligodT-anchor 1 μL,Ex Taq DNA 聚合酶0.25 μL,10×Ex Taq Buffer 2.5 μL,dNTP mix(2.5 mmol/L)2.5 μL,加水補足25 μL。PCR反應條件:94℃ 30 s,61℃ 45 s,72℃ 1 min,5個循環(huán);94℃30 s,58℃ 45 s,72℃ 1 min,5個循環(huán);94℃ 30 s,56℃ 45 s,72℃ 1 min,18個循環(huán),共計28個循環(huán)。第二輪PCR:首輪PCR產物稀釋25倍 1 μL,5anchor 1 μL,5-RACE-D3(10 μmol/L)1 μL,Ex Taq DNA聚合酶0.25 μL,10×Ex Taq Buffer 2.5 μL,dNTP mix(2.5 mmol/L)2.5 μL,加水補足25 μL。反應條件同第一輪PCR。本實驗所用引物序列,見表1。

表1 本實驗所用引物

1.2.5 斜帶石斑魚IL-10生物信息學分析 利用NCBI在線軟件ORF finder(http://www.ncbi.nlm.nih. gov/gorf/orfig.cgi)確定開放閱讀框(ORF)及蛋白氨基酸序列;運用ExPASy ProtParam(http://web. expasy.org/protparam/)預測蛋白分子量(Mw)和理論等電點(p I);采用SignalP 4.1 Server(http://www. cbs.dtu.dk/services/SignalP/)預測對信號肽進行預測;通過TMHMM Server v.2.0(http://www.cbs.dtu. dk/services/TMHMM/)分析是否存在跨膜結構域;功能位點分析使用ExPASy PROSITE在線軟件(http:// prosite.expasy.org/);蛋白結構預測采用PredictProtein(https://www.predictprotein.org/)、Phyre 2 server(http://www.sbg.bio.ic.ac.uk/phyre2/ htm l/page. cgi?id=index)和SWISS Model(http://swissmodel.ex pasy.org/);用DNAMAN6.0軟件對斜帶石斑魚IL-10氨基酸序列與其他物種IL-10進行同源比對分析及序列拼接;利用MEGA6軟件進行系統(tǒng)進化樹的構建。

1.2.6 斜帶石斑魚IL-10原核表達載體pET-32a-IL10構建及其原核表達條件的優(yōu)化 根據(jù)已得到的IL-10基因cDNA 全長設計含酶切位點的引物擴增IL-10去信號肽核心片段,引物為IL10F(Bam H I)及IL10R(Xho I)。將用Ex Taq酶PCR擴增獲得的含酶切位點的基因片段連接到pMD18-T 載體,對測序正確的陽性克隆和pET-32a(+)載體進行用Bam H I和Xho I雙酶切,1.5%瓊脂糖電泳后回收相應的目的條帶,用T4 DNA連接酶進行連接后轉化至大腸桿菌BL21感受態(tài)細胞中,對陽性克隆進行測序。用終濃度為0.1 mmol/L IPTG誘導4 h進行原核表達。融合蛋白表達的優(yōu)化按黃瑜等[18]的方法進行。

2 結果

2.1 斜帶石斑魚IL-10基因克隆

結果(圖1)顯示,獲得全長1 104 bp的IL-10基因,ORF為564 bp,編碼187個氨基酸。GenBank登錄號為KJ741852。蛋白理論分子量(Mw)為21.7 kD,理論等電點(p I)為5.74,存在一個22個氨基酸的信號肽,無跨膜區(qū)。ExPASy預測斜帶石斑魚IL-10蛋白存在4個蛋白激酶C磷酸化位點(24、73-75、135-137、177-179),3個酪蛋白激酶II磷酸化位點(69-72、85-88、150-153),1個N端酰基化位點(146-149)及1個N端糖基化位點(156-161)。

2.2 斜帶石斑魚IL-10氨基酸同源性及進化分析

多序列比對結果(圖2)顯示,斜帶石斑魚IL-10的氨基酸序列與吉富羅非魚的同源性較高,為72.25%,而與人和猴的僅有26.18%和25.13%。在IL-10多肽鏈的30、80、110和116位存在4個非常保守的半胱氨酸。運用MEGA6 軟件,以N-J法構建IL-10的系統(tǒng)進化樹,結果(圖3)顯示,斜帶石斑魚IL-10與維多利亞湖慈鯛(P. nyererei)、斑馬宮麗魚(M. zebra)和吉富羅非魚(O. niloticus)聚為一支。

2.3 斜帶石斑魚IL-10的結構分析

軟件分析顯示IL-10蛋白二級結構中有117處α-螺旋,占總二級結構的62.57%;42處無規(guī)則卷曲,占總二級結構的22.46%;14處β-轉角和14處共有延伸鏈,各占總二級結構的7.49%。PredictProtein和SWISS MODEL的結果(圖4)顯示,斜帶石斑魚IL-10蛋白存在6個α螺旋,三維結構以二聚體形式存在,其空間結構與人IL-10相似。

圖1 IL-10 cDNA核苷酸及其編碼的氨基酸序列

2.4 斜帶石斑魚IL-10原核表達載體pET-32a-IL10構建及原核表達條件優(yōu)化

對測序正確的陽性菌與pET-32a(+)質粒的空載體菌分別在0.1 mmol/L IPTG的濃度下誘導4 h進行預誘導,結果(圖5-A)顯示分子量為37.5 kD的融合蛋白在大腸桿菌BL21中表達。在IPTG為0.1 mmol/L濃度時,融合蛋白在3 h后表達量不再增加(圖5-B);在IPTG濃度為0.02、0.06、0.1、0.4、0.6及1 mmol/L,誘導4 h后,蛋白在0.02 mmol/L和0.6 mmol/L時表達量最大(圖5-C);在IPTG濃度為0.1 mmol/L時,溫度為16℃、28℃和37℃,誘導4 h后,37℃時蛋白在沉淀中表達量最大,而28℃時,上清中融合蛋白較多(圖5-D)。

3 討論

本研究成功克隆了全長為1 104 bp的斜帶石斑魚IL-10基因,ORF為564 bp,編碼187個氨基酸,在3'端非編碼區(qū)存在2個mRNA不穩(wěn)定基序(ATTTA),表明該基因可能在炎癥反應和免疫應答中較為活躍[19]。在斜帶石斑魚IL-10氨基酸序列中存在LLDESVEESFKSPFGCHAMNSILEFYLNTVL和GLFKAMGDLDVLFNYIE兩個屬于IL-10特征基序L-[FILMV]-X3-[ILV]-X3-[FILMV]-X5-C-X5-[ILMV]和G-X2-KA-X2-[ED]-X-D-[ILV]-[FLY]-[FILMV]-X2-[ILMV]-[EKQZ]的序列。4個非常保守的半胱氨酸(Cys-30、Cys-80、Cys-110和Cys-116)在斜帶石斑魚IL-10多肽鏈中被發(fā)現(xiàn),在人和哺乳動物IL-10研究表明,IL-10存在對蛋白空間結構的形成非常重要的兩個二硫鍵,對蛋白行使生物學功能不可或缺[8,9]。不同于人和哺乳動物IL-10,斜帶石斑魚IL-10成熟肽還存在一個額外的半胱氨酸(Cys-26),在其他已報道的魚類IL-10,如斑馬魚[13]、鯉[12]、虹鱒[17]等,也存在類似的現(xiàn)象,IL-10通常以二聚體行使生物學功能,這一可能的二硫鍵形成位點,是否對IL-10的蛋白具有生物學意義還有待研究。哺乳動物IL-10通過磷酸化行使其一些生物學功能[20],在鼠IL-10發(fā)現(xiàn)一個潛在的糖基化位點,但不是蛋白功能所必需的[21]。斜帶石斑魚IL-10蛋白分析顯示其存在4個蛋白激酶C磷酸化位點、3個酪蛋白激酶II磷酸化位點以及1個N端酰基化位點和1個N端糖基化位點,這些修飾位點對蛋白功能是否有影響可以通過蛋白點突變來驗證。

圖2 IL-10氨基酸序列的同源性比對

圖3 IL-10氨基酸的系統(tǒng)進化樹

圖4 預測斜帶石斑魚IL-10(A)與人的IL-10(B)的空間結構結構比較

圖5 IL-10蛋白的原核表達

氨基酸同源性比對顯示斜帶石斑魚IL-10的氨基酸序列與吉富羅非魚的同源性最高,為72.25%,而與人類和猴的同源性較低。斜帶石斑魚IL-0蛋白結構預測顯示斜帶石斑魚IL-10含有6個α螺旋,以二聚體形式存在于生物體中,與對人和哺乳動物的研究一致。其三級結構與人IL-10空間結構類似,蛋白的結構與蛋白功能之間的關系非常密切,斜帶石斑魚IL-10可能和哺乳動物功能相同,以二聚體形式與其受體相結合,通過Jak-Sata信號通路來實現(xiàn)其對細胞的免疫調控[7,10,20]。Grayfer等[16,22]基于對金魚和斑馬魚IL-10及其受體IL10R特征和功能的研究,也認為IL-10的功能在進化上可能高度保守。

本實驗成功表達出分子量為37.5 kD的斜帶石斑魚IL-10融合蛋白。選擇pET-32a(+)質粒構建原核表達載體是因為pET系列擁有T7強啟動子的載體,能短時間內大量表達目的蛋白,而表達的融合蛋白帶有6個組氨酸標簽有利于蛋白后續(xù)的分離與純化,且多數(shù)情況下表達的融合蛋白不影響蛋白的免疫原性和免疫反應性。原核表達條件優(yōu)化顯示,在IPTG濃度為0.02 mmol/L,37℃誘導3 h,蛋白包涵體的表達量最大。這些結果為后續(xù)進行IL-10蛋白的富集及抗體的制備奠定基礎。而斜帶石斑魚IL-10融合蛋白在低溫誘導時,上清也有明顯表達,可以省去包涵體蛋白復性等相關復雜操作,獲取活性蛋白。本研究的后續(xù)實驗將進行IL-10的分離純化,參照Grayfera等[16]的實驗操作,研究其對細胞呼吸爆發(fā)的影響,以及IL-10對細胞內一些免疫相關基因表達的調控,探討IL-10在石斑魚免疫反應中的作用與哺乳動物是否類似。

4 結論

本研究從斜帶石斑魚成功克隆了IL-10基因。該基因全長1 104 bp,編碼187個氨基酸,理論分子量為21.7 kD,理論等電點為5.74,同源分析顯示與吉富羅非魚親緣關系較近。并對斜帶石斑魚IL-10的重組表達和表達條件進行了探索,其優(yōu)化條件為:0.02 mmol/L的IPTG,37℃誘導3 h。

[1]Fiorentino DF, Bond MW, Mosmann TR. Two types of mouse T helper cell. IV. Th2 clones secrete a factor that inhibits cytokine production by Th1 clones[J]. Exp Med, 1989, 170(6):2081-2095.

[2]Saraia M, O' Garra A. The regulation of IL-10 production by immune cells[J]. Nat Rev Immunol, 2010, 10(3):170-181.

[3]Bogdan C, Vodovotz Y, Nathan C. Macrophage deactivation by interleukin10[J]. Exp Med, 1991, 174(6):1549-1555.

[4]Oswald IP, Wynn TA, Sher A, et al. Interleukin 10 inhibits macrophage microbicidal activity by blocking the endogenous production of tumor necrosis factor alpha required as a costimulatory factor for interferon gamma-induced activation[J]. Proc Natl AcadSci U S A, 1992, 89(18):8676-8680.

[5]Ding L, Shevach EM. IL-10 inhibits mitogen-induced T cell proliferation by selectively inhibiting macrophage costimulatory function[J]. Immunol, 1992, 148(10):3133-3139.

[6]Lang R, Patel D, Morris JJ, et al. Shaping gene expression in activated and resting primary macrophages by IL-10[J]. Immunol,2002, 169(5):2253-2263.

[7]Zdanov A, Schalk-Hihi C, Gustchina A, et al. Crystal structure of interleukin-10 reveals the functional dimer with an unexpected topological similarity to interferon gamma[J]. Structure, 1995, 3(6):591-601.

[8]W illiam TW, Rosalinda S, Anthony T, et al. Disulfide bond assignments and secondary structure analysis of human and marine interleukin 10[J]. Biochemistry, 1993, 32(34):8807-8815.

[9] Syto R, Murgolo NJ, Braswell EH, et al. Structural and biological stability of the human interleukin 10 homodimer[J]. Biochemistry,1998, 37(48):16943-16951.

[10]Weber-Nordt RM, Riley JK, Greenlund AC, et al. Stat3 recruitment by two distinct ligand-induced, tyrosinephosphorylated docking sites in the interleukin-10 receptor intracellular domain[J]. Biol Chem, 1996, 271(44):27954-27961.

[11]Zou J, Clark MS, Secombes CJ. Characterisation, expression and promoter analysis of an interleukin 10 homologue in the puffer fish,F(xiàn)ugu rubripes[J]. Immunogenetics, 2003, 55(5):325-335.

[12]Savan R, Igawa D, Sakai M. Cloning, characterization and expression analysis of interleukin-10 from the common carp,Cyprinus carpio L. [J]. Biochem, 2003, 270(23):4647-4654.

[13] Zhang DC, Shao YQ, Huang YQ, et al. Cloning, characterization and expression analysis of interleukin-10 from the zebrafish(Danio rerion)[J]. Biochem Mol Bio, 2005, 38(5):571-576.

[14] Pinto RD, Nascimento DS, Reis MR, et al. Molecular characterization, 3D modelling and expression analysis of sea bass(Dicentrarchus labrax L.)interleukin-10[J]. Mol Immunol, 2007, 44(8):2056-2065.

[15] Ingerslev HC, Ronneseth A, Pettersen EF, et al. Differential expression of immune genes in Atlantic salmon(Salmo salar L.)challenged intraperitoneally or by cohabitation with IPNV. Scand.[J]. Immunol, 2009, 69(2):90-98.

[16] Grayfer L, Hodgkinson J, Hitchen SJ, et al. Characterization and functional analysis of goldfish(Carassius auratus L.)interleukin-10[J]. Molecular Immunology, 2011, 48(4):563-571.

[17]Inoue Y, Kamota S, Ito K, et al. Molecular cloning and expression analysis of rainbow trout(Oncorhynchus mykiss)interleukin-10 cDNAs[J]. Fish Shellfish Immunol, 2005, 18(4):335-344.

[18]黃瑜, 張雪利, 魯義善, 等. 紅笛鯛 tdt 基因融合蛋白原核表達條件的優(yōu)化及純化[J]. 廣東海洋大學學報, 2013, 33(1):44-49.

[19]Gay E, Babajko S. AUUUA sequences compromise human insulinlike growth factor binding protein-1 mRNA stability[J]. Biochemical and Biophysical Research Communications, 2000, 267(2):509-515.

[20] Sabat R, Grütz G, Warszawsk K, et al. Biology of interleukin-10[J]. Cytokine & Growth Factor Reviews, 2010, 21(5):331-344.

[21] Moore KW, O’Garra A, de Waal Malefyt R, et al. Interleukin-10[J]. Annu Rev Immunol, 1993, 11:165-190.

[22] Grayfer L, Belosevic M. Identification and molecular characterization of the interleukin-10 receptor 1 of the zebrafish(Danio rerio)and the goldfish(Carassius auratus L.)[J]. Developmental and Comparative Immunology, 2012, 36(2):408-417.

(責任編輯 馬鑫)

Cloning and Prokaryotic Expression of IL-10 Gene of Epinephelus coioides

Xiao Zhouting1,2,3Huang Yucong1,2,3Jian Jichang1,2,3Lu Yishan1,2,3Wu Zaohe2,3,4
(1. College of Fisheries,Guangdong Ocean University,Zhanjiang524088;2. Guangdong Provincial Key Laboratory of Pathogenic Biology and Epidemiology for Aquatic Economic Animals,Zhanjiang524088;3. Key Laboratory of Diseases Controlling for Aquatic Economic Animals of Guangdong Higher Education Institutions,Zhanjiang524088;4. Zhongkai University of Agriculture and Engineering,Guangzhou510225)

Interleukin 10(IL-10)is the critical multiple-functional cytokine playing the major role of mediating inflammatory responses and regulating the differentiation and proliferation of some immune cells. The full-length cDNA of IL-10 of Epinephelus coioides was cloned using homological cloning and rapid amplification of cDNA ends(RACE-PCR). Results showed the full-length of IL-10’s cDNA was 1 104 bp,containing an open reading frame of 564 bp encoding 187 amino acids with an estimated molecular weight of 21.7 kD and an estimated isoelectric point of 5.74, and it had a signal peptide of 22 amino acid residues and no transmembranous region. Amino acid homology analysis showed that the amino acid sequences of IL-10 with Oreochromis niloticus had the highest homology and up to 72.25%. SDS-PAGE indicated that the molecular weight of the fusion protein was 37.5 kD which was in accord with the estimated. The optimal condition of inducible expression for IL-10 protein in Escherichia coli was at 37℃ for 3 h with 0.02 mmol/L of IPTG.

Epinephelus coioides;Interleukin-10;gene cloning;prokaryotic expression

10.13560/j.cnki.biotech.bull.1985.2015.09.018

2014-12-07

廣東教育廳科技創(chuàng)新重點項目(2012CXZD0026),十二五國家科技支撐計劃(2012BAD17B03)

肖周婷,女,碩士研究生,研究方向:水產動物病害防治;E-mail:xihu1tingyu@163.com

簡紀常,教授,研究方向:水產動物免疫學及病害控制;E-mail:jianjc@gmail.com

主站蜘蛛池模板: 97成人在线视频| 国产精品吹潮在线观看中文| 精品国产电影久久九九| 成人福利在线观看| 日韩无码白| 中美日韩在线网免费毛片视频| 欧美日韩在线亚洲国产人| 国产欧美日韩一区二区视频在线| 91久久性奴调教国产免费| 国产成人欧美| 天堂网亚洲系列亚洲系列| 在线观看国产精品第一区免费 | 国产乱码精品一区二区三区中文| 中文字幕在线视频免费| 99久久精品国产麻豆婷婷| 青草91视频免费观看| 中文字幕在线视频免费| 亚洲精品无码日韩国产不卡| 精品伊人久久大香线蕉网站| 日本伊人色综合网| 日韩精品一区二区三区大桥未久 | 在线精品欧美日韩| 欧美亚洲中文精品三区| 一级爆乳无码av| 无码福利视频| 五月天丁香婷婷综合久久| 国产精品yjizz视频网一二区| 亚洲AⅤ无码日韩AV无码网站| 亚洲国产成人精品无码区性色| 国产女人水多毛片18| 在线网站18禁| 欧美成人国产| 2021精品国产自在现线看| 国产麻豆精品手机在线观看| 国产精品成人免费视频99| 久久大香香蕉国产免费网站| 亚洲第一视频区| 国产在线视频导航| 国产亚洲一区二区三区在线| 国产在线观看91精品亚瑟| 91国语视频| 国产网友愉拍精品| 青草国产在线视频| 国产午夜人做人免费视频中文 | 中文字幕日韩久久综合影院| 99国产精品国产| 国产精品人成在线播放| 欧美日韩中文国产va另类| 国产黄色爱视频| 无码精品一区二区久久久| 一级毛片基地| 风韵丰满熟妇啪啪区老熟熟女| 中文无码伦av中文字幕| 亚洲女同一区二区| 中国毛片网| 人妻丰满熟妇AV无码区| 久操中文在线| 999国产精品| 六月婷婷综合| 久久综合色天堂av| 国产一区亚洲一区| 欧美在线导航| 中国国产高清免费AV片| 六月婷婷激情综合| 成人在线观看一区| 91热爆在线| 一区二区理伦视频| 99视频精品在线观看| 国产无套粉嫩白浆| 中国特黄美女一级视频| 91精品啪在线观看国产91九色| 亚洲三级电影在线播放| 18禁高潮出水呻吟娇喘蜜芽| 国产欧美在线观看一区| 欧美19综合中文字幕| 国产呦视频免费视频在线观看| 国产精品香蕉在线观看不卡| 亚洲国产亚洲综合在线尤物| 久久无码高潮喷水| 欧美国产综合视频| 99久久精品国产综合婷婷| 免费无码网站|