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

代謝工程改造大腸桿菌合成β-丙氨酸

2017-06-21 15:10:24梁姍姍周麗張斌周哲敏
食品與發酵工業 2017年5期
關鍵詞:優化

梁姍姍,周麗,張斌,周哲敏

(江南大學 生物工程學院,工業生物技術教育部重點實驗室,江蘇 無錫,214122)

代謝工程改造大腸桿菌合成β-丙氨酸

梁姍姍,周麗,張斌,周哲敏*

(江南大學 生物工程學院,工業生物技術教育部重點實驗室,江蘇 無錫,214122)

β-丙氨酸是天然存在的唯一一種β型氨基酸,在醫藥、食品、化工等領域有重要應用。該研究考察了以重組Escherichiacoli發酵合成β-丙氨酸的可能性。在敲除副產物乙酸、甲酸、乙醇、琥珀酸、乳酸合成途徑編碼基因的EscherichiacoliCICIM B0016-050 (ackA-ptapflBadhEfrdAldhA)菌株中,考察疊加敲除β-丙氨酸的競爭途徑天冬氨酸激酶、泛酸合成酶和葡萄糖轉運蛋白(EⅡCBGlc)的編碼基因(lysC、panC、ptsG)以及過表達來源于Corynebacteriumglutamicum的panD基因對合成β-丙氨酸的影響。結果表明,疊加敲除上述基因后,β-丙氨酸的合成量依次提高了12.5%、39.3%和13.3%;過量表達panD基因,β-丙氨酸合成量提高了86.2倍;經發酵條件優化,菌株 B0016-080BB/pPL-panD搖瓶發酵可合成(5.0±0.2) g/L β-丙氨酸,體積生產強度達到(0.12±0.01) g/(L·h)。該結果為發酵法合成β-丙氨酸奠定了基礎。

β-丙氨酸;大腸桿菌;L-天冬氨酸-α-脫羧酶;代謝工程;基因敲除

β-丙氨酸即3-氨基丙酸,是自然界中唯一存在的β型氨基酸,近年來,在醫藥、化工、食品、環境等領域的作用日漸突出[1]。在醫藥領域,主要用于合成泛酸以及泛酸鈣,此外,β-丙氨酸還是合成肌肽(camosine)、帕米酸鈉(pamidronate,第二代抗骨溶解藥)、巴柳氮(balsalazide)等藥物的原料;在化工領域,用于藥劑制備中的沉淀劑,電鍍緩蝕劑,鉛中毒解毒劑等;在食品領域,可用作調味品;在環境領域用于水的凈化絮凝劑等[2]。有報道稱β-丙氨酸是未來全球12種最具開發潛力的三碳化工產品之一[3]。

國內外生產β-丙氨酸主要采用化學合成法,包括丙烯酸法、丙烯腈法、β-氨基丙腈法、β-氨基丙醇法等,這些方法合成工藝流程長,且需高溫、高壓,能耗大、生產過程對環境有一定的毒害作用。β-丙氨酸生物生產主要是采用微生物如產堿桿菌屬Alcaligenessp. OMT-MY14、氨基桿菌屬AminobacteraminobranceATCC 23314和紅球菌Rhodococcussp. G20[4],將β-氨基丙腈進行全細胞生物轉化;KONST等人還利用L-天門冬氨酸-α-脫羧酶的酶轉化法制備β-丙氨酸[5-6],成為近年來研究的熱點[6],但其生產需要昂貴的前體和大量酶,還存在著高濃度底物抑制菌生長、產率較低等問題,不適合大規模工業化生產。Chan Woo Song等[7]構建了菌株E.coliCWF4NA2 (ΔiclRΔfumCΔfumAΔfumBΔptsGΔlacIPaspA::Ptrc Pacs::Ptrc),首次采用發酵法在E.coli中合成β-丙氨酸,搖瓶發酵產量達3.94 g/L,但并未敲除有機酸、醇等副產物的合成途徑,這些副產物的積累會導致β-丙氨酸得率降低,同時影響產物的化學純度;此外,也未考察β-丙氨酸競爭代謝途徑(lysC)和分解代謝途徑(panC)的敲除對β-丙氨酸合成的影響;另外,以葡萄糖為碳源時,敲除ptsG[8]可顯著增強丙酮酸激酶(Pck)和磷酸烯醇式丙酮酸羧化酶(Ppc)的活性,從而提高草酰乙酸的積累量,而該基因在其它碳源存在時對β-丙氨酸合成的影響也未有報道。ZHOU等用PL啟動子成功地對乳酸脫氫酶[9]和丙氨酸脫氫酶[10]進行開關控制,僅通過溫控調節,在重組E.coli中實現了D-乳酸和L-丙氨酸的高產,將該方法應用于β-丙氨酸合成重組菌株,也有望大幅度提高β-丙氨酸的合成水平。

本文在實驗室已有菌株E.coliB0016-050[11](敲除副產物乙酸、甲酸、乙醇、琥珀酸、乳酸合成途徑編碼基因ackA-pta、pflB、adhE、frdA、ldhA)的基礎上,敲除lysC、panC和ptsG基因,并將來源于Corynebacteriumglutamicum的panD基因克隆于溫控強啟動子PL下游,在上述重組菌株中進行高效表達,考察這些代謝途徑的改造對發酵甘油合成β-丙氨酸的影響。

圖1 E.coli菌株β-丙氨酸合成相關代謝途徑及改造策略Fig.1 Metabolic pathway and strategy for β- alanine production in E.coli注:基因及其編碼的酶: pflB, 丙酮酸甲酸裂解酶; ldhA, NAD依賴發酵型D-乳酸脫氫酶; pta, 磷酸轉乙酰酶; ackA, 乙酸激酶; adhE, 乙醇脫氫酶; frdA, 富馬酸還原酶; aspA, 天冬氨酸氨裂解酶; lysC, 天冬氨酸激酶; panC, 泛酸合成酶; ptsG, 主要的葡萄糖轉運蛋白IICB(Glc); ppc, 磷酸烯醇丙酮酸羧化酶; aspC, 天冬氨酸轉氨酶; panD, 天冬氨酸-α-脫羧酶·黑點代表本研究所刪除的基因

1 材料與方法

1.1 菌株與引物

1.1.1 菌株和質粒

本文所用到的菌株、質粒如表1所示。

E.coliB0016-050菌株為在野生型菌株E.coliCICIM B0016中敲除了基因ackA-pta、pflB、adhE、frdA、ldhA的重組菌株[11],由本實驗室構建、保藏。其他菌株由本文構建。

1.1.2 引物

本文所用到的引物如表2所示。

1.2 培養基

(1)LB培養基(g/L):胰蛋白胨10;酵母浸出粉5;NaCl 10(固體培養基,瓊脂15)。

(2)M9-1液體培養基(g/L):NaHPO4·12H2O 15.1;KH2PO43.0;NH4Cl 1.0;NaCl 0.5;(NH4)2SO413.21,補加50% (v/v)甘油 10 mL,200 g/L泛酸1 mL,0.1% (v/v) MgSO41 mol/L和0.1% (v/v)的微量元素母液。微量元素成分為(g/L):FeCl3·6H2O 2.4;CoCl2·6H2O 0.3;CuCl2·2H2O 0.15;ZnCl20.3;Na2MO4·2H2O 0.3;H3BO30.075;MnCl2·4H2O 0.495。

1.3 基因敲除方法

用 DATSENKO 等[12]報道的Red 重組系統基因敲除方法依次敲除B0016-050菌株染色體上的lysC、panC、和ptsG基因。其具體過程為: 以pKD13質粒為模板,分別用引物對LysC-PKD13R/LysC-PKD13F、PanC-PKD13R/PanC-PKD13F和PtsG -PKD13R/PtsG-PKD13F進行PCR 擴增,獲得兩端各具有50 bp 同源臂序列、中部為FRT-Kan-FRT 的突變盒片段。借助質粒pKD46 上的Red重組系統將突變盒整合于待突變菌株的染色體上,并用相應的驗證引物對(YlysCF/YlysCR,YpanCF/YpanCR,YptsGF/YptsGR)進行菌落PCR 驗證,最后用pCP20質粒去除重組菌中的卡那抗性基因,并進一步用相應驗證引物進行驗證。

表1 本研究中所使用的菌株和質粒

注:a:江南大學中國高校工業微生物資源和信息中心;b:耶魯大學基因保藏中心。

表2 本研究中所使用的引物

注a:小寫字母表示與pKD13質粒同源堿基序列。

1.4 搖瓶發酵培養條件

種子培養:取-80 ℃甘油保藏的菌種,在氨芐青霉素抗性的LB平板劃線,于33 ℃ 恒溫培養24 h后,轉接到 50 mL含有氨芐青霉素的LB液體培養基中,于33 ℃搖床200 r/min振蕩培養至OD600=2.5。

發酵培養:將LB種子培養基轉接至50 mL M9-1培養基中,使OD600終濃度為0.05,于33 ℃、200 r/min振蕩培養至OD600=2.5(優化時OD600取:0、0.5、1、1.5、2、2.5、3、3.5、4),補料1.5 mL 50% (v/v)甘油并用100 g/L的NaHCO3調節pH=7,同時調節溫度至42 ℃誘導,繼續200 r/min振蕩培養,每6 h補加1 mL 50% (v/v)甘油,并用100 g/L的NaHCO3調節pH=7,48 h結束搖瓶發酵,每次試驗采用3個平行并計算其平均值。

1.5 β-丙氨酸含量測定

發酵液預處理方法:發酵液樣品10 000 r/min離心2 min,取上清液加入等體積100 g/L三氯乙酸,避光放置大于3 h,適當稀釋后用0.45 μm濾膜過膜。利用鄰苯二甲醛(OPA)試劑在線衍生β-丙氨酸,用高效液相色譜(HPLC)分析β-丙氨酸的含量。色譜柱為Diamonsil C18柱(250 cm×4.6 mm, 5 μm),流動相A為:無水醋酸鈉(分析純)4.52 g,先溶于水,再加三乙胺 200 μL,四氫呋喃5 mL,水定容至1 L,醋酸調節pH 7.2;流動相B為:無水醋酸鈉4.52 g 水定容為200 mL,醋酸調節pH 7.2,加400 mL甲醇和400 mL乙腈。采用梯度洗脫:0~27 min,B液由8%上升到60%,流速為1 mL/min;20~31.5 min,B液由60%上升到100%,流速為1 mL/min ;32~34 min,B液梯度不變,流速為1.2 mL/min。34~35 min,B液由100%下降到8%,流速為1 mL/min。柱溫為40 ℃,檢測波長為338 nm;進樣量為10 μL。

2 結果和分析

2.1lysC、panC和ptsG基因的敲除

用 DATSENDO 等[12]報道的基因敲除方法,在菌株E.coliB0016-050中依次疊加敲除天冬氨酸激酶、泛酸合成酶和葡萄糖轉運蛋白(EⅡCBGlc)的編碼基因lysC、panC和ptsG(圖1),依次獲得了B0016-060A、B0016-070B和B0016-080BB菌株(表1)。用驗證引物對上述突變進行驗證,結果如圖2所示,lysC敲除前基因大小為1 714 bp,敲除后為469 bp;panC敲除前基因大小為1 334 bp,敲除后為587 bp;ptsG敲除前基因大小為1 770 bp,敲除后為340 bp。表明上述基因已成功刪除。

M-DL 5 000 DNA maker;1-lysC敲除前;2-lysC敲除后;3-panC敲除前;4-panC敲除后;5-ptsG敲除前;6-ptsG敲除后圖2 基因敲除PCR驗證電泳圖Fig.2 PCR verification of gene knockouts

2.2 pPL-panD重組質粒的構建

本實驗室前期構建的質粒pET28a-panD[14]用BglⅡ和EcoR I進行雙酶切,膠回收0.5 kb片段,克隆于pPL451質粒載體的BamH I和EcoR I酶切位點,獲得重組質粒pPL-panD,其物理圖譜如圖3A所示。重組質粒pPL-panD用EcoR I進行酶切,電泳獲得了4.8 kb 條帶(圖3B),表明質粒構建正確。

M-DL 10 000 DNA maker;1- pPL-panD digested by EcoR I圖3 pPL-panD重組質粒的物理圖譜(A)和酶切驗證電泳圖譜(B)Fig.3 The map of the recombinant plasmid pPL- panD (A) and confirmation of pPL-panD by restriction digestion (B)

2.3 泛酸添加量的優化

敲除panC基因的B0016-070B菌株為營養缺陷型,不能在無泛酸的基本培養基上生長。向M9-1培養基中添加不同濃度的泛酸(0、10-1、10-2、10-3、10-4、10-5g/L),菌體生長狀況如圖4所示。在不添加泛酸時,B0016-070B菌株基本不生長;添加量為10-4、10-5g/L時,菌體生長受到一定程度的抑制;泛酸添加量為10-2g/L時,B0016-070B菌體生長情況與未敲除panC基因的B0016-060A菌株的菌體生長情況相似,且再提高泛酸濃度,菌體生長情況不變。因此,后續實驗在培養敲除了panC基因的菌株時,基本培養基中添加10-2g/L泛酸。

圖4 B0016-070B 菌株泛酸添加量的優化Fig.4 Optimization of pantothenic acid concentration in B0016-070B

2.4 代謝途徑敲除對β-丙氨酸合成的影響

通過搖瓶發酵實驗,比較B0016-060A、B0016-070B和B0016-080BB菌株的β-丙氨酸合成水平。結果如圖5所示,B0016-050菌株的β-丙氨酸產量為5.2 mg/L,B0016-060A菌株的β-丙氨酸合成量較B0016-050提高了35.1%,說明敲除lysC基因有助于L-天冬氨酸的積累,從而提高β-丙氨酸的產量;B0016-070B菌株的β-丙氨酸合成量較B0016-060A提高了595.8%,表明敲除panC基因阻斷了β-丙氨酸的分解,從而使其合成量大幅度提高;B0016-080BB菌株的β-丙氨酸產量為49.3 mg/L,較B0016-070B提高了1.7%,說明敲除ptsG基因對β-丙氨酸的合成也有一定影響。

圖5 基因敲除對β-丙氨酸合成的影響Fig.5 Effects of gene deletions on β-alanine synthesis

2.5panD基因的異源過表達對β-丙氨酸合成的影響

將pPL-panD質粒分別轉化入菌株B0016-050、B0016-060A、B0016-070B、B0016-080BB,獲得B0016-050/pPL-panD、B0016-060A/pPL-panD、B0016-070B/pPL-panD、B0016-080BB/pPL-panD重組株,進一步進行發酵實驗檢測其β-丙氨酸合成性能。

如圖6所示,過量表達panD基因顯著提高了β-丙氨酸合成量,B0016-050/pPL-panD菌株的β-丙氨酸合成量達2.4 g/L,比B0016-050重組菌株提高了464.1倍;B0016-080BB/pPL-panD重組菌β-丙氨酸合成量最高,達4.3 g/L,比B0016-080BB重組菌β-丙氨酸合成量提高了86.2倍。敲除lysC、panC和ptsG基因的B0016-060A/pPL-panD、B0016-070B/pPL-panD和B0016-080BB/pPL-panD菌株的β-丙氨酸合成量依次升高12.5%、39.3%和13.3%。可見,敲除天冬氨酸的分解途徑關鍵基因lysC和β-丙氨酸的分解途徑關鍵基因panC顯著提高了β-丙氨酸的合成量;此外,也證明了在以甘油為碳源時,敲除ptsG基因同樣可將增加β-丙氨酸代謝流。

圖6 panD基因的過量表達對不同菌株發酵合成β-丙氨酸的影響Fig.6 Effect of over expression of panD on the synthesis of β-alanine in different strains

2.6 B0016-080BB/pPL-panD菌株搖瓶發酵優化

2.6.1 誘導時機的優化

ZHOU等優化了PL啟動子調控的發酵過程,其最佳的菌體生長溫度為33 ℃,而最佳的產物合成溫度為42 ℃[9-10]。本文進一步對菌株B0016-080BB/pPL-panDβ-丙氨酸發酵過程中,這兩個溫度的轉折時機進行優化,結果如圖7所示。當菌體生長至OD600為2.5時,將發酵溫度由33 ℃切換為42 ℃,β-丙氨酸合成量最高,達到4.9 g/L。因此,最佳溫度誘導時機為OD600值達到2.5。

圖7 B0016-080BB/pPL-panD菌株溫度誘導時機的優化Fig.7 Optimization of thermo-induction of strain B0016-080BB/pPL-panD

2.6.2 發酵周期的優化

進一步對菌株B0016-080BB/pPL-panD的β-丙氨酸發酵周期進行優化。測定發酵液中的β-丙氨酸產量和菌體濃度,結果如圖8所示。當搖瓶發酵42 h時,β-丙氨酸合成量最高,達到5.0±0.2 g/L,體積生產強度達到(0.12±0.005) g/(L·h)。因此,菌株B0016-080BB/pPL-panD的最佳搖瓶發酵周期為42 h。

圖8 B0016-080BB/pPL-panD菌株搖瓶發酵時間的優化Fig.8 Optimization of fermentation time of strain B0016-080BB/pPL-panD

3 結論

本文在去除了有機酸、醇副產物的重組E.coliB0016-050菌株中,刪除β-丙氨酸競爭代謝途徑(lysC)和分解途徑基因(panC),可顯著提高β-丙氨酸的積累量;并證明了在以甘油為碳源時,敲除ptsG基因同樣可提高β-丙氨酸合成量。此外,用溫控型PL強啟動子調控谷氨酸棒桿菌(Corynebacteriumglutamicum)來源的L-天冬氨酸α-脫羧酶的表達,可顯著提高重組E.coli的β-丙氨酸合成水平。優化后B0016-080BB/pPL-panD菌株β-丙氨酸合成量達到(5.0±0.2) g/L,體積生產強度達到(0.12±0.005) g/(L·h)。

[1] JANG Y S,KIM B,SHIN J H,et al.Bio-based production of C2-C6 platform chemicals[J].Biotechnology and Bioengineering, 2012,109(10):2 437-2 459.

[2] LEE J W,NA D,PARK J M,et al.Systems metabolic engineering of microorganisms for natural and non-natural chemicals[J].Nature Chemical Biology,2012,8(6):536-546.

[3] STEUNENBERGP,K NST P M,SCOTT E L,et al.Polymerisation of β-alanine through catalytic ester-amide exchange[J].European Polymer Journal,2013,49(49):1 773-1 781.

[4] LIANG Lu-yi, ZHENG Yu-guo, SHEN Yin-chun.Optimization of β-alanine production from β-aminopropionitrile by resting cells ofRhodococcussp. G20 in a bubble column reactor using response surface methodology[J].Process Biochemistry,2008,43(7):758-764.

[5] K NST P M,FRANSSEN M C R,SCOTT E L,et al.A study on the applicability ofL-aspartate α-decarboxylase in the biobased production of nitrogen containing chemicals[J]. Green Chemistry,2009,11(10):1 646-1 652.

[6] BORODINA I,KILDEGAARD K R,JENSEN N B,et al.Establishing a synthetic pathway for high-level production of 3-hydroxypropionic acid inSaccharomycescerevisiaevia β-alanine[J].Metabolic Engineering,2014,27:57-64.

[7] SONG C W,LEE J,KO Y S,et al.Metabolic engineering ofEscherichiacolifor the production of 3-aminopropionic acid.[J].Metabolic Engineering,2015,30(3):121-129.

[8] 呂忠成.Ptsg敲除對重組大腸桿菌碳氮代謝的影響.[D].西安:西北大學,2015.

[9] ZHOU Li,NIU Dan-dan,TIAN Kang-ming,et al.Genetically switchedD-lactate production inEscherichiacoli[J]. Metabolic Engineering,2012,14(5):560-568.

[10] ZHOU Li,DENG Can,CUI Wen-jing,et al.EfficientL-alanine production by a thermo-regulated switch inEscherichiacoli[J].Applied Biochemistry and Biotechnology,2016,178(2):324-337.

[11] 周麗,鄧璨,崔文璟,等.溫度調節基因開關調控大腸桿菌發酵合成L-丙氨[J].微生物學通報,2015,42(11):2 272-2 281.

[12] DATSENKO K A,WANNER B L.One-step inactivation of chromosomal genes inEscherichiacoliK-12 using PCR products[J].Proceedings of the National Academy of Sciences of the United States of America, 2000,97(12): 6 640-6 645.

[13] SONDEREGGER M,JEPPSSON M,HAHN-H GERDAL B,et al.Molecular basis for anaerobic growth ofSaccharomycescerevisiaeon xylose, investigated by global gene expression and metabolic flux analysis[J].Applied & Environmental Microbiology,2004,70(4):2 307-2 317.

[14] 石增秀,崔文璟,周麗,等.谷氨酸棒桿菌L-天冬氨酸α-脫羧酶基因的克隆及重組酶性質研究[J].生物技術通報,2013(4):110-115.

Metabolic engineering ofEscherichiacolifor the production of β-alanine

LIANG Shan-shan, ZHOU Li,ZHANG Bin, ZHOU Zhe-min*

(Schoool of Biotechnology and the Key Laboratory of Industrial Biotechnology, Jiangnan University, Wuxi 214122, China)

β-Alanine, the only kind of β amino acid which exists naturally, has important applications in medicine, food and chemical industry. In this study, the possibility of metabolic engineering ofEscherichiacolifor the production of β-alanine was investigated. TheEscherichiacoliCICIM B0016-050 (ackA-ptapflBadhEfrdAldhA) strain with deletions on synthetic pathways of acetic acid, formic acid, ethanol, succinic acid and lactic acid, was used as the starting strain. Effects of gene knockouts oflysC(encoding aspartokinase),panC(encoding pantothenate synthetase) andptsG(encoding the major glucose transporter IICB (Glc)) and overexpression of theCorynebacteriumglutamicumpanDgene (encodingL-aspartate-α-decarboxylase) on synthesis of β-alanine in B0016-050 were investigated. The synthesis of β-alanine was improved by 12.5%, 39.3% and 13.3%, respectively, after combined gene knockouts, and was increased by more than 100.6 times after overexpressing thepanDgene. After optimization of fermentation conditions, the yield of β-alanine in the shake flask fermentation of strain B0016-080BB/pPL-panDreached 5.0±0.2 g/L, and the volume productivity was 0.12±0.01 g/(L·h). The results laid the foundation for fermentation of β-alanine.

β-alanine;Escherichiacoli;L-aspartate-α-decarboxylase;metabolic engineering;gene knockout

10.13995/j.cnki.11-1802/ts.201705003

碩士研究生(周哲敏教授為通訊作者,E-mail:zhmzhou@jiangnan.edu.cn)。

國家自然科學基金項目(31300087);基本科研-重點項目(JUSRP51611A)

2017-01-05,改回日期:2017-02-09

猜你喜歡
優化
超限高層建筑結構設計與優化思考
房地產導刊(2022年5期)2022-06-01 06:20:14
PEMFC流道的多目標優化
能源工程(2022年1期)2022-03-29 01:06:28
民用建筑防煙排煙設計優化探討
關于優化消防安全告知承諾的一些思考
一道優化題的幾何解法
由“形”啟“數”優化運算——以2021年解析幾何高考題為例
圍繞“地、業、人”優化產業扶貧
今日農業(2020年16期)2020-12-14 15:04:59
事業單位中固定資產會計處理的優化
消費導刊(2018年8期)2018-05-25 13:20:08
4K HDR性能大幅度優化 JVC DLA-X8 18 BC
幾種常見的負載均衡算法的優化
電子制作(2017年20期)2017-04-26 06:57:45
主站蜘蛛池模板: 午夜国产在线观看| 欧美成人精品一级在线观看| 欧美视频二区| 久久99久久无码毛片一区二区 | 国产精品久久久久久久久kt| 亚洲三级视频在线观看| 91视频青青草| 精品人妻无码中字系列| 亚洲欧美成人影院| 日日碰狠狠添天天爽| 72种姿势欧美久久久久大黄蕉| 午夜久久影院| 激情无码字幕综合| 久久人人爽人人爽人人片aV东京热 | 麻豆国产精品| 中文精品久久久久国产网址| 激情综合婷婷丁香五月尤物| 日本欧美成人免费| 亚洲高清在线播放| 欧美日韩另类在线| 国产aaaaa一级毛片| 亚洲精品午夜天堂网页| 天天躁夜夜躁狠狠躁躁88| 美女国内精品自产拍在线播放 | 中文字幕天无码久久精品视频免费| 天天躁日日躁狠狠躁中文字幕| 国产粉嫩粉嫩的18在线播放91| 午夜视频免费一区二区在线看| 国产成人欧美| 97无码免费人妻超级碰碰碰| 第一区免费在线观看| 亚洲性影院| 国产二级毛片| 精品国产自在现线看久久| 精品欧美一区二区三区久久久| 青青草原国产| 成人欧美日韩| 91系列在线观看| 人人91人人澡人人妻人人爽 | 精品视频91| 国模粉嫩小泬视频在线观看| 中文字幕永久在线看| 麻豆国产精品| 亚洲色图综合在线| 四虎综合网| 91免费精品国偷自产在线在线| 青青青国产视频手机| 精品国产成人高清在线| 亚卅精品无码久久毛片乌克兰| 在线观看亚洲天堂| 青青青国产视频手机| 毛片久久网站小视频| 日韩人妻精品一区| 精品国产成人三级在线观看| 久久夜色精品国产嚕嚕亚洲av| AV在线天堂进入| 91精品啪在线观看国产91九色| 丰满少妇αⅴ无码区| 亚洲男人的天堂久久香蕉| 亚洲久悠悠色悠在线播放| AV无码一区二区三区四区| 中文字幕久久精品波多野结| 国国产a国产片免费麻豆| 国产成人精品一区二区三区| 日韩少妇激情一区二区| 国内老司机精品视频在线播出| 国产簧片免费在线播放| 无码电影在线观看| 国产一级毛片在线| 国产91精品久久| 91福利一区二区三区| 欧美日韩另类国产| 人妖无码第一页| 久久综合伊人77777| 久久人妻系列无码一区| 在线观看国产黄色| 欧美在线观看不卡| 国产毛片高清一级国语| 久久久久久尹人网香蕉| 国产精品人莉莉成在线播放| 一级毛片高清| 国产激情在线视频|