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

    谷氨酸棒狀桿菌葡萄糖代謝阻斷工程菌的構(gòu)建

    2017-02-08 07:42:35韓武洋劉金雷杜紅燕王北辰李天明
    食品科學 2017年2期
    關(guān)鍵詞:碳源基因組質(zhì)粒

    韓武洋,劉金雷,杜紅燕,王北辰,儀 宏,李天明,*

    谷氨酸棒狀桿菌葡萄糖代謝阻斷工程菌的構(gòu)建

    韓武洋1,劉金雷1,杜紅燕1,王北辰2,儀 宏1,李天明1,*

    (1.河北科技大學生物科學與工程學院,河北 石家莊 050018;2.威斯康星大學農(nóng)業(yè)與生命學院,美國 麥迪遜 53706)

    采用反向代謝工程的策略,以谷氨酸棒桿菌ATCC 13032野生型為出發(fā)菌株,利用無抗性標記的同源重組方法,敲除了編碼葡萄糖pts系統(tǒng)關(guān)鍵轉(zhuǎn)運蛋白基因ptsG、ptsH-ptsI和葡萄糖轉(zhuǎn)運系統(tǒng)關(guān)鍵轉(zhuǎn)運蛋白基因abc、abc2和iolt1,得到了5 株逐次敲除了相應基因的突變株。結(jié)果表明:當以葡萄糖為唯一碳源培養(yǎng)時,CGΔptsG菌株的葡萄糖的消耗是野生型的50%,菌體OD值為1.473;CGΔptsH-ptsI菌株的葡萄糖的消耗是野生型的39.5%,菌體OD值為1.226;CGΔabc菌株的葡萄糖的消耗是野生型的36%,菌體OD值為1.092;CGΔabc2菌株的葡萄糖的消耗是野生型的26.2%,菌體OD值為0.486;CGΔiolt1葡萄糖的消耗和菌體生長OD值為0,實現(xiàn)了谷氨酸棒狀桿菌葡萄糖代謝的阻斷,說明ptsG、ptsH-ptsI、abc、abc2和iolt1所編碼的轉(zhuǎn)運蛋白具有葡萄糖轉(zhuǎn)運功能。

    谷氨酸棒狀桿菌;葡萄糖;基因敲除

    谷氨酸棒狀桿菌(Corynebacterium glutamicum)是重要的工業(yè)生產(chǎn)菌之一,在50多年的發(fā)展過程中,主要被用于工業(yè)生產(chǎn)風味增強劑和動物飼料中的添加劑,如谷氨酸、賴氨酸、谷氨酰胺等[1-3]。C. glutamicum可利用多種糖類和蜜糖作為原料進行工業(yè)發(fā)酵生產(chǎn)[4],葡萄糖是C. glutamicum工業(yè)生產(chǎn)氨基酸等產(chǎn)品中使用最多的碳源,雖然關(guān)于C. glutamicum葡萄糖代謝及其調(diào)控的生理研究已經(jīng)有了一些報道和進展,但在葡萄糖運輸系統(tǒng)及分子水平上調(diào)節(jié)的研究相對較少[5]。近年來,C. glutamicum ATCC 13032全基因組測序的完成[6-7]和其他組學數(shù)據(jù)的積累[8],使科研人員可以利用反向代謝工程的技術(shù)手段[9-10],進一步探明葡萄糖運輸系統(tǒng)和分子水平上的調(diào)節(jié),為利用代謝工程改良菌種提供理論指導,以實現(xiàn)提高葡萄糖的利用效率、促進菌體生長的目標,或者實現(xiàn)不利用葡萄糖與其他利用葡萄糖微生物共發(fā)酵生產(chǎn)的目標。

    綜合目前的國內(nèi)外研究,認為谷氨酸棒狀桿菌中存在著兩條不同的葡萄糖轉(zhuǎn)運途徑[11]。一是當谷氨酸棒狀桿菌以葡萄糖為唯一碳源生長時,大部分葡萄糖會在磷酸轉(zhuǎn)移酶系統(tǒng)(phosphotransferase system,pts)的作用下磷酸化成6-磷酸葡萄糖,然后進入糖酵解途徑[12];二是少部分的葡萄糖會先以葡萄糖分子的形態(tài)進入細胞內(nèi),再在葡萄糖磷酸激酶的作用下生成6-磷酸葡萄糖進入糖酵解途徑[13-14]。根據(jù)文獻及基因組信息,可能的涉及轉(zhuǎn)運相關(guān)的基因及運輸途徑如圖1所示,但是對于有些基因的功能還不十分清楚,需要進一步研究確證。

    圖11 C. glutamiiccuumm ATCC 13032葡萄糖代謝途徑[[2200]]Fig. 1 Schematic illustration of the glucose metabolic pathways of C. glutamicum ATCC 13032

    本實驗以C. glutamicum ATCC 13032野生型為出發(fā)菌株,利用無抗性標記的同源重組方法,逐個敲除了調(diào)控葡萄糖pts和葡萄糖轉(zhuǎn)運系統(tǒng)關(guān)鍵酶基因ptsG、ptsH-ptsI、abc、abc2和iolt1,實現(xiàn)C. glutamicum葡萄糖代謝系統(tǒng)的阻斷,為在分子水平研究葡萄糖轉(zhuǎn)運與代謝提供參考依據(jù)。同時,也為C. glutamicum與其他葡萄糖利用微生物共發(fā)酵提供了基礎(chǔ)遺傳資源。

    1 材料與方法

    1.1 材料

    1.1.1 菌株與質(zhì)粒

    本實驗所用的菌株、質(zhì)粒分別如表1所示。

    表1 本實驗所用的菌株、質(zhì)粒Table1 Strains and plasmids used in this study

    1.1.2 試劑

    表2 本實驗所用引物Table2 PCR primers used in this study

    High-Fidelity DNA聚合酶 北京全式金生物技術(shù)有限公司;限制性內(nèi)切酶XhoⅠ、Hind Ⅲ、BamHⅠ和EcoRⅠ 美國NEB公司;T4 DNA 連接酶 大連寶生物工程公司;基因組提取試劑盒、質(zhì)粒小提中量試劑盒、超薄DNA產(chǎn)物純化試劑盒、瓊脂糖凝膠DNA回收試劑盒 天根生化科技(北京)有限公司;聚合酶鏈式反應(polymerase chain reaction,PCR)引物由生工生物工程(上海)股份有限公司合成(表2)。

    1.1.3 培養(yǎng)基

    LB培養(yǎng)基:胰蛋白胨10 g/L、酵母提取物5 g/L、NaCl 10 g/L,調(diào)pH值為6.5左右,需要時加入卡那霉素至終質(zhì)量濃度30 mg/L。固體培養(yǎng)基則添加瓊脂粉至20 g/L。

    成分確定(chemical difined,CD)培養(yǎng)基:KH2PO41 g/L、(NH4)2SO43 g/L、NH4NO31 g/L、MgSO40.5 g/L、CaCl20.2 g/L、Fe2(SO4)30.01 g/L、MnSO40.01 g/L、CuSO40.002 g/L、ZnSO40.001 g/L、生物素0.002 g/L,pH 6.5左右。

    孵育培養(yǎng)基:CD培養(yǎng)基、乳酸鈉50 g/L、乙酸鈉50 g/L,pH 6.5左右。

    1.2 儀器與設(shè)備

    HWS型培養(yǎng)箱 寧波江南儀器廠;SW-CJ型超凈工作臺 蘇凈集團安泰公司制造;PCR儀、電轉(zhuǎn)儀美國Bio-Rad公司;Neofuge23R冷凍離心機 力新儀器(上海)有限公司;pH211酸度計 意大利Hanna公司;Jl003A型電子天平 常熟雙杰測試儀器廠;DZKW-D水浴鍋 河北黃驊航天儀器廠;HYG型搖床 上海欣蕊自動化設(shè)備有限公司;F1-F2凝膠成像儀上海天能科技有限公司;PoweWave HT酶標儀 美國伯騰儀器有限公司;BCD-191W/HC冰箱 上海海爾股份有限公司;G80W23YCsL-03(RO)微波爐 上海申安醫(yī)療器械廠;Finnpipette移液槍 美國Thermo公司。

    1.3 方法

    1.3.1 菌株培養(yǎng)

    C. glutamicum ATCC 13032在30 ℃、200 r/min條件下培養(yǎng),E. coli TransT1在37 ℃、200 r/min條件下培養(yǎng)。

    1.3.2 打靶質(zhì)粒的構(gòu)建

    ptsG基因敲除打靶質(zhì)粒構(gòu)建方法如下:首先以C. glutamicum ATCC 13032野生型菌株基因組為模板,利用PCR擴增出待敲除基因ptsG上游和下游同源臂片段,利用重疊PCR技術(shù)將ptsG基因的上下游同源臂融合,將所得片段純化后與pK18mobsacB空質(zhì)粒利用XhoⅠ和EcoRⅠ進行雙酶切,酶切體系(50 μL)為:XhoⅠ 1.5 μL、EcoRⅠ 1.5 μL、10×buffer 5 μL、模板2 000 ng/μL和超純水,37 ℃條件下反應3 h。目的片段與pK18mobsacB空質(zhì)粒酶切產(chǎn)物通過超薄DNA產(chǎn)物回收試劑盒回收,再通過T4連接酶16 ℃過夜連接。連接體系(25 μL)為:目的片段酶切產(chǎn)物5 μL、pK18mobsacB空質(zhì)粒酶切產(chǎn)物1 μL、10×buffer 1.25 μL、T4連接酶0.5 μL、超純水17.25 μL。連接過程應在冰上操作并充分混勻。連接產(chǎn)物轉(zhuǎn)化E. coli TransT1感受態(tài)中,轉(zhuǎn)化子提取質(zhì)粒,經(jīng)酶切驗證后獲得打靶質(zhì)粒pK18mobsacB/ΔptsG。構(gòu)建示意圖如圖2所示。

    圖2 pK18mobsaaccBB/ΔppttssGG構(gòu)建示意圖Fig. 2 Schematic illustration of construction of pK18mobsacB/ΔptsG

    本實驗所涉及的敲除ptsH-ptsI、abc、abc2和iolt1基因打靶質(zhì)粒利用類似方法構(gòu)建,分別獲得打靶質(zhì)粒pK18mobsacB/ΔptsH-ptsI、pK18mobsacB/Δabc、pK18mobsacB/Δabc2和pK18mobsacB/Δiolt1。

    1.3.3 C. glutamicum感受態(tài)的制備及轉(zhuǎn)化

    C. glutamicum感受態(tài)細胞制備方法參考Jang等[15],將過夜培養(yǎng)的菌株按體積分數(shù)10%接種于三角瓶中培養(yǎng),當OD600nm達到0.35~0.45時,于4℃條件下4 000 r/min離心5 min收集菌體,4 ℃條件下利用10%甘油吹懸漂洗2 次,取約5 μg打靶質(zhì)粒與40 μL感受態(tài)細胞輕輕混合,轉(zhuǎn)入電擊杯,靜置5 min后進行電擊(電擊條件:2.5 kV/cm電壓、250 Ω電阻、25 μF電容、2 mm電擊杯)。電擊完成后,30 ℃、180 r/min條件下孵育培養(yǎng)4 h,涂布于含有卡那霉素抗性的固體培養(yǎng)基上,30 ℃條件下培養(yǎng)約40 h,固體培養(yǎng)基平板上長出均勻單菌落,挑取單菌落進行篩選。

    1.3.4 工程菌株的篩選及鑒定

    將單菌落擴培至含卡那霉素抗性的LB固體平板上,經(jīng)菌落PCR驗證,確定打靶質(zhì)粒經(jīng)同源重組整合到染色體上。選取驗證正確的轉(zhuǎn)化子,在無抗性壓力條件下振蕩培養(yǎng)6 h,稀釋涂布于含有10%蔗糖的LB固體平板上,利用蔗糖致死基因sacB反向篩選出發(fā)生第2次同源重組的轉(zhuǎn)化子,挑取蔗糖平板上單菌落,分別轉(zhuǎn)接至卡那霉素抗性平板及10%蔗糖平板上,能在蔗糖平板上生長但不能在卡那霉素抗性平板生長的菌株即為發(fā)生第2次重組的菌株,但可能是突變菌株和恢復突變菌株兩種可能,通過設(shè)計相對應的驗證引物來進行菌落PCR驗證,從而獲得陽性工程菌。

    1.3.5 工程菌株搖瓶發(fā)酵

    從新鮮活化的LB固體平板上刮取菌株于種子培養(yǎng)基中,30 ℃條件下培養(yǎng)14 h,將種子液4 000 r/min離心2 min收集菌體,并用0.9%生理鹽水洗滌1~2 次,接種于含2%葡萄糖的CD培養(yǎng)基中,轉(zhuǎn)接前后OD600nm保持不變。培養(yǎng)溫度30 ℃,搖床轉(zhuǎn)速180 r/min,pH 7.0左右,發(fā)酵過程中間歇性調(diào)節(jié)pH值,使發(fā)酵液pH值保持在7.0左右,并于不同時間檢測菌體濃度并取樣保留,利用葡萄糖測定儀測定樣品中殘留的葡萄糖濃度。

    2 結(jié)果與分析

    2.1 ptsG基因缺失菌株構(gòu)建及葡萄糖利用特性

    2.1.1 ptsG基因敲除打靶質(zhì)粒構(gòu)建

    圖3 PCR產(chǎn)物(A)和打靶質(zhì)粒pK18mobsaaccBB/ΔppttssGG酶切驗證(BB)Fig. 3 PCR products (A) and identification of targeting plasmid pK18mobsacB/ΔptsG by enzymatic digestion (B)

    ptsG是葡萄糖pts特異酶EⅡ的編碼基因。以C. glutamicum ATCC 13032基因組DNA為模板,經(jīng)PCR擴增獲得上下同源臂,利用重疊PCR融合獲得目的片段,大小為2 900 bp,結(jié)果如圖3A所示。將目的片段與空質(zhì)粒pK18mobsacB連接,獲得重組質(zhì)粒pK18mobsacB/ΔptsG,經(jīng)酶切驗證,得到5 600 bp載體片段和2 900 bp目的片段兩條條帶,結(jié)果如圖3B所示,說明重組質(zhì)粒pK18mobsacB/ΔptsG構(gòu)建成功。

    2.1.2 工程菌株CGΔptsG的篩選及驗證

    圖4 同源交換原理示意圖Fig. 4 Schmatic illustration of the principle of homologous exchange

    將重組質(zhì)粒pK18mobsacB/ΔptsG電轉(zhuǎn)化至C. glutamicum ATCC 13032感受態(tài)細胞中,利用兩次同源重組方法和sacB基因反向篩選技術(shù)來實現(xiàn)基因的無痕敲除,基因重組過程如圖4所示。

    圖5 ptsGptsG缺失基因PCR驗證PCRFig. 5 Identification of deletion of ptsG by PCR

    按照C. glutamicum ATCC 13032基因組序列中待敲除基因ptsG的基因序列設(shè)計敲除驗證引物QPTSGYZF/ QPTSGYZR,分別以CGΔptsG菌株與C. glutamicum ATCC 13032菌株的基因組DNA為模板,進行PCR驗證,結(jié)果如圖5所示,C. glutamicum ATCC 13032菌株基因組能擴增出大小為2 502 bp的ptsG基因片段,與理論值大小一致,而CGΔptsG菌株基因組沒有擴增出相應大小的片段,并且空白陰性對照(超純水為模版)無條帶,說明CGΔptsG工程菌基因組中的ptsG基因已被敲除。

    2.1.3 工程菌株CGΔptsG的生長狀況研究

    CGΔptsG和C. glutamicum ATCC 13032在以葡萄糖為唯一碳源的條件下CD培養(yǎng)基中搖瓶發(fā)酵48 h,測定生長曲線和葡萄糖濃度,結(jié)果如圖6所示,敲除ptsG基因的突變株與野生型相比在延滯期生長速率開始降低,生長速率緩慢,而葡萄糖濃度升高,表明敲除ptsG基因使菌株對葡萄糖利用率降低,但仍能利用葡萄糖繼續(xù)生長。工程菌株CGΔptsG的生長速率、菌液OD600nm和葡萄糖利用率分別是野生菌株的74.4%、77.7%和50%。

    圖6 C. glutamicum ATCC 13032 和C GppttssGG在以葡萄糖為碳源培養(yǎng)基上的生長曲線(A)及葡萄糖利用情況(BB)Fig. 6 Growth curves and glucose consumption of C. glutamicum ATCC 13032 and CGΔptsG on glucose as the sole carbon source

    2.2 ptsH、ptsI基因簇缺失菌株構(gòu)建及葡萄糖利用特性

    2.2.1 ptsH-ptsI基因簇敲除打靶質(zhì)粒構(gòu)建

    圖7 PCR產(chǎn)物(A)和打靶質(zhì)粒pK18mobsacB/ΔppttssHH-ppttssII酶切驗證(BB)Fig. 7 PCR products (A) and identification of targeting plasmid pK18mobsacB/ΔptsH-ptsI by enzymatic digestion (B)

    葡萄糖pts還與特異性蛋白E1(由ptsH編碼)和搬運蛋白HPr(由ptsI編碼)相關(guān),在C. glutamicum ATCC 13032基因組中,p t s H和p t s I基因是一個基因簇。以C. glutamicum ATCC 13032基因組DNA為模板,分別擴增出ptsH-ptsI基因簇的上下游同源臂,利用重疊PCR融合獲得目的片段,大小為2 700 bp,結(jié)果如圖7A所示。將目的片段與空質(zhì)粒pK18mobsacB連接,獲得重組質(zhì)粒pK18mobsacB/ΔptsH-ptsI,經(jīng)酶切驗證,得到5 600 bp載體片段和2 700 bp目的片段兩條條帶,結(jié)果如圖7B所示,說明重組質(zhì)粒pK18mobsacB/ΔptsH-ptsI構(gòu)建成功。

    2.2.2 工程菌株CGΔptsH-ptsI的篩選及驗證

    將重組質(zhì)粒pK18mobsacB/ΔptsH-ptsI電轉(zhuǎn)化至CGΔptsG感受態(tài)細胞中,利用兩次同源重組方法和sacB基因反向篩選技術(shù)來實現(xiàn)基因的無痕敲除,基因重組過程與圖4相似。

    按照C. glutamicum ATCC 13032基因組序列中待敲除基因ptsH-ptsI的基因序列設(shè)計敲除驗證引物QPTSHIYZF/ QPTSHIYZR,分別以CGΔptsH-ptsI菌株與CGΔptsG菌株的基因組DNA為模板,進行PCR驗證,結(jié)果如圖8所示,CGΔptsG菌株基因組能擴增出大小為1 700 bp的ptsH-ptsI基因片段,與理論值大小一致,而CGΔptsH-ptsI菌株基因組沒有擴增出相應大小的片段,并且空白陰性對照(超純水為模版)無條帶,說明CGΔptsH-ptsI工程菌基因組中的ptsH-ptsI基因已被敲除。

    圖8 ptsHptsH-ptsIptsI缺失基因PCR驗證PCRFig. 8 Identification of deletion of ptsH-ptsI by PCR

    2.2.3 工程菌株CGΔptsH-ptsI的生長狀況研究

    圖9 C. glutamicum ATCC 13032和CGptsH--ppttssII在以葡萄糖為碳源培養(yǎng)基上的生長曲線(A)及葡萄糖利用情況(BB)Fig. 9 Growth curves and glucose consumption of C. glutamicum ATCC 13032 and CGΔptsH-ptsI on glucose as the sole carbon source

    CGΔptsH-ptsI和C. glutamicum ATCC 13032在以葡萄糖為唯一碳源的條件下CD培養(yǎng)基中搖瓶發(fā)酵48 h,測定菌株生長曲線和殘留的葡萄糖,結(jié)果如圖9所示,敲除ptsH-ptsI基因的突變株與野生型相比在延滯期生長速率開始降低,生長速率緩慢,而葡萄糖殘量升高,表明敲除ptsH-ptsI基因使菌株對葡萄糖利用率降低。CGΔptsH-ptsI的生長速率、菌液OD600nm和葡萄糖利用率分別是野生型菌株的59.7%、66.8%和39.5%。

    2.3 abc基因缺失菌株構(gòu)建及葡萄糖利用特性

    2.3.1 abc基因敲除打靶質(zhì)粒構(gòu)建

    圖10 PCR產(chǎn)物(A)和打靶質(zhì)粒pK18mobsacB/Δaabbcc酶切驗證(BB)Fig. 10 The PCR products (A) and identification of targeting plasmid pK18mobsacB/Δabc by digestion (B)

    2.3.2 工程菌株CGΔabc的篩選及驗證

    圖11 abcabc缺失基因PCR驗證PCRFig. 11 Identification of deletion of abc by PCR

    將重組質(zhì)粒pK18mobsacB/Δabc電轉(zhuǎn)化至CGΔptsH-ptsI感受態(tài)細胞中,利用兩次同源重組方法和sacB基因反向篩選技術(shù)來實現(xiàn)基因的無痕敲除,基因重組過程與圖4相似。

    按照C. glutamicum ATCC 13032基因組序列中待敲除基因abc的基因序列設(shè)計敲除驗證引物QABCYZF/ QABCYZR,分別以CGΔabc菌株與CGΔptsH-ptsI菌株的基因組DNA為模板,進行PCR驗證,結(jié)果如圖11所示,CGΔptsH-ptsI菌株基因組能擴增出大小為1 000 bp的abc基因片段,與理論值大小一致,而CGΔabc菌株基因組沒有擴增出相應大小的片段,并且空白陰性對照(超純水為模版)無條帶,說明CGΔabc工程菌基因組中的abc基因已被敲除。

    2.3.3 工程菌株CGΔabc的生長狀況研究

    CGΔabc和C. glutamicum ATCC 13032在以葡萄糖為唯一碳源的條件下CD培養(yǎng)基中搖瓶發(fā)酵48 h,測定菌株生長曲線和殘留的葡萄糖,結(jié)果如圖12所示,CGΔabc的生長速率、菌液OD600nm和葡萄糖利用率分別是野生型菌株的53.6%、60%和36%,說明敲除abc基因后,工程菌株對葡萄糖的利用率進一步降低。

    圖12 C. glutamiiccuumm ATCC 13032和CCGGΔaabbcc在以葡萄糖為碳源培養(yǎng)基上的生長曲線(A)及葡萄糖利用情況(BB)Fig. 12 Growth curves and glucose consumption of C. glutamicum ATCC 13032 and CGΔabc on glucose as the sole carbon source

    2.4 abc2基因缺失菌株構(gòu)建及葡萄糖利用特性

    2.4.1 abc2基因敲除打靶質(zhì)粒構(gòu)建

    圖13 PCR產(chǎn)物(A)和打靶質(zhì)粒pK18mobsacB/Δabc2abc2酶切驗證(B)Fig. 13 PCR products (A) and identification of targeting plasmid pK18mobsacB/Δabc2 by enzymatic digestion (B)

    abc2基因是C. glutamicum ATCC 13032另一個abc型糖轉(zhuǎn)運系統(tǒng)的相關(guān)基因。以C. glutamicum ATCC 13032基因組DNA為模板,分別擴增出abc基因的上下游同源臂基因片段,利用重疊PCR融合獲得目的片段,大小為2 600 bp,結(jié)果如圖13A所示。將目的片段與空質(zhì)粒pK18mobsacB連接,獲得重組質(zhì)粒pK18mobsacB/Δabc2,經(jīng)酶切驗證,得到5 600 bp載體片段和2 600 bp目的片段兩條條帶,結(jié)果如圖13B所示,說明重組質(zhì)粒pK18mobsacB/Δabc2構(gòu)建成功。

    2.4.2 工程菌株CGΔabc2的篩選及驗證

    圖14 ptsabc2sabc2缺失基因PCR驗證PCRFig. 14 Identification of deletion of ptsabc2 by PCR

    將重組質(zhì)粒pK18mobsacB/Δabc2電轉(zhuǎn)化至CGΔabc感受態(tài)細胞中,利用兩次同源重組方法和sacB基因反向篩選技術(shù)來實現(xiàn)基因的無痕敲除,基因重組過程與圖4相似。

    按照C. glutamicum ATCC 13032基因組序列中待敲除基因abc2的基因序列設(shè)計敲除驗證引物QABC2YZF/ QABC2YZR,分別以CGΔabc2菌株與CGΔabc菌株的基因組DNA為模板,進行PCR驗證,結(jié)果如圖14所示,CGΔabc菌株基因組能擴增出大小為1 130 bp的abc2基因片段,與理論值大小一致,而CGΔabc2菌株基因組沒有擴增出相應大小的片段,并且空白陰性對照(超純水為模版)無條帶,說明CGΔabc2工程菌基因組中的abc2基因已被敲除。

    2.4.3 工程菌株CGΔabc2的生長狀況研究

    CGΔabc2和C. glutamicum ATCC 13032在以葡萄糖為唯一碳源的條件下CD培養(yǎng)基中搖瓶發(fā)酵48 h,測定菌株生長曲線和殘留的葡萄糖,結(jié)果如圖15所示,CGΔabc2的生長速率、菌液OD600nm和葡萄糖利用率分別是野生型菌株的10.3%、30.1%和26.2%,說明abc2基因?qū)τ谡{(diào)控葡萄糖轉(zhuǎn)運有重要作用,敲除abc2基因后,工程菌株對葡萄糖的利用率大大降低。

    圖15 C. glutamicum ATCC 13032和CGaabbcc22在以葡萄糖為碳源培養(yǎng)基上的生長曲線(A)及葡萄糖利用情況(BB)Fig. 15 Growth curves and glucose consumption of C. glutamicum ATCC 13032 and CGΔabc2 on glucose as the sole carbon source

    2.5 iolt1基因缺失菌株構(gòu)建及葡萄糖利用特性

    2.5.1 iolt1基因敲除打靶質(zhì)粒構(gòu)建

    圖16 PCR產(chǎn)物(A)和打靶質(zhì)粒pK18mobsacB/Δiioolltt11酶切驗證(BB)Fig. 16 PCR products (A) and identification of targeting plasmid pK18mobsacB/Δiolt1 by enzymatic digestion (B)

    在基因組中iolt1基因是編碼肌醇轉(zhuǎn)運的基因,但在Krings等[16]報道中,其可以轉(zhuǎn)運葡萄糖。以C. glutamicum ATCC 13032基因組DNA為模板,分別擴增出abc基因的上下游同源臂基因片段,利用重疊PCR融合獲得目的片段,大小為1 800 bp,結(jié)果如圖16A所示。將目的片段與空質(zhì)粒pK18mobsacB連接,獲得重組質(zhì)粒pK18mobsacB/ Δiolt1,經(jīng)酶切驗證,得到5 600 bp載體片段和1 800 bp目的片段兩條條帶,結(jié)果如圖16B所示,說明重組質(zhì)粒pK18mobsacB/Δiolt1構(gòu)建成功。

    2.5.2 工程菌株CGΔiolt1的篩選及驗證

    圖17 ptsGptsG缺失基因PCR驗證PCRFig. 17 Identification of deletion of ptsGptsG by PCR

    將重組質(zhì)粒pK18mobsacB/Δiolt1電轉(zhuǎn)化至CGΔabc2感受態(tài)細胞中,利用兩次同源重組方法和sacB基因反向篩選技術(shù)來實現(xiàn)基因的無痕敲除,基因重組過程與圖4相似。

    按照C. glutamicum ATCC 13032基因組序列中待敲除基因iolt1的基因序列設(shè)計敲除驗證引物QIOLT1YZF/ QIOLT1YZR,分別以CGΔiolt1菌株與CGΔabc2菌株的基因組DNA為模板,進行PCR驗證,結(jié)果如圖17所示,CGΔabc2菌株基因組能擴增出大小為1 480 bp的iolt1基因片段,與理論值大小一致,而CGΔiolt1菌株基因組沒有擴增出相應大小的片段,并且空白陰性對照(超純水為模版)無條帶,說明CGΔiolt1工程菌基因組中的iolt1基因已被敲除。

    2.5.3 工程菌株CGΔiolt1在以葡萄糖為碳源培養(yǎng)基上的生長情況

    圖18 C. glutamicum ATCC 13032和CGiioolltt11在以葡萄糖為碳源培養(yǎng)基上的生長曲線(A)及葡萄糖利用情況(BB)Fig. 18 Growth curves and glucose consumption of C. glutamicum ATCC 13032 and CGΔiolt1 on glucose as the sole carbon source

    CGΔiolt1和C. glutamicum ATCC 13032在以葡萄糖為唯一碳源的條件下CD培養(yǎng)基中搖瓶發(fā)酵48 h,測定菌株生長曲線和殘留的葡萄糖,結(jié)果如圖18所示,CGΔiolt1幾乎無法利用葡萄糖生長,生長速率和殘留葡萄糖量接近于0。

    2.6 工程菌株CGΔptsG、CGΔptsH-ptsI、CGΔabc、CGΔabc2、CGΔiolt1與野生型菌株在以葡萄糖為碳源培養(yǎng)基的生長情況

    圖19 固體培養(yǎng)基上生長情況Fig. 19 Growth of C. glutamicum ATCC 13032 and engineered strains on glucose-containing solid minimal medium

    將工程菌株CGΔptsG、CGΔptsH-ptsI、CGΔabc、CGΔabc2、CGΔiolt1與野生型菌株依次接種于在僅以葡萄糖為碳源的固體CD培養(yǎng)基上,生長情況如圖19所示,原始菌株生長良好,隨著葡萄糖轉(zhuǎn)運系統(tǒng)的相關(guān)基因ptsG、ptsH-ptsI、abc、abc2的依次敲除,菌株生長越來越弱,直至基因iolt1敲除后,工程菌株不生長。

    圖20 8 h和24 h菌液ODD60000nnmmFig. 20 OD600nmof fermentation broths at 8 and 24 h

    工程菌株與野生型菌株同時接種于液體CD培養(yǎng)基中進行發(fā)酵培養(yǎng),培養(yǎng)過程中,分別在8、24 h取樣測OD600nm,結(jié)果如圖20所示,基因ptsG和ptsH、ptsI基因簇敲除后OD600nm降低,生長速率下降,但并沒有完全阻斷葡萄糖代謝,說明在C. glutamicum ATCC 13032中除葡萄糖轉(zhuǎn)運除了pts,還存在其他方式的轉(zhuǎn)運系統(tǒng)。基因abc和abc2敲除后OD600nm降低幅度大,生長速率下降幅度大;基因iolt1敲除后,菌株不生長,在敲除pts的基礎(chǔ)上,abc、abc2和iolt1 3個基因的敲除,完全阻斷了葡萄糖代謝,說明C. glutamicum ATCC 13032中,abc、abc2和iolt1 3個基因是負責葡萄糖分子轉(zhuǎn)運的關(guān)鍵基因。

    3 討 論

    據(jù)報道,C. glutamicum ATCC 13032中含有4 組pts,包括葡萄糖pts、果糖pts、蔗糖pts和一個未知pts[17-20]。葡萄糖pts由特異性蛋白E1(由ptsH編碼)、搬運蛋白HPr(由ptsI編碼)和酶Ⅱ(EⅡs)組成[21]。與C. glutamicum ATCC 13032不同,另一株野生型菌株C. glutamicum R含有兩組β-葡糖苷pts,他們負責β-葡糖苷的攝取,例如水楊苷、熊果苷和葡萄糖[22]。本研究獲得的ptsG和ptsI-ptsH基因缺失突變株,在以葡萄糖為唯一碳源培養(yǎng)基上的生長和葡萄糖的消耗都受到影響,說明在C. glutamicum ATCC 13032中,葡萄糖轉(zhuǎn)運的pts除了ptsG編碼的蛋白EⅡGlc,還包括特異性蛋白E1(由ptsH編碼)和搬運蛋白HPr(由ptsI編碼),與現(xiàn)有的文獻報道一致[21]。但是突變株仍保留了運輸葡萄糖糖的能力,說明葡萄糖的轉(zhuǎn)運除了pts外還存在其他轉(zhuǎn)運系統(tǒng)。

    據(jù)報道C. glutamicum能夠有效的使用葡萄糖酸和核糖作為唯一碳源,但他們的攝取是由非pts介導。葡萄糖酸由gntP編碼的特定葡萄糖酸透酶GntP調(diào)控吸收,再經(jīng)gntK編碼的葡萄糖酸激酶GntK磷酸化成6-磷酸葡萄糖酸進入戊糖磷酸途徑?;騡ntP或gntK的缺失都會影響在葡萄糖上的生長[22]。還有報道,核糖先通過rbsACBD基因編碼的特定轉(zhuǎn)運系統(tǒng)RbsACBD轉(zhuǎn)運到細胞內(nèi),再通過rbsK1和rbsK2基因編碼的核糖激酶RbsK1和RbsK2磷酸化形成5-磷酸核糖進入戊糖磷酸途徑,rbsK1和rbsK2雙基因缺失會影響在核糖上的生長的[24-25],但沒有文獻具體報道負責葡萄糖原型轉(zhuǎn)運的基因。本研究根據(jù)C. glutamicum ATCC 13032的KEGG數(shù)據(jù)庫,找到了類似葡萄糖酸透酶和abc型糖轉(zhuǎn)運系統(tǒng)的相關(guān)基因abc1和abc2基因(圖1)。為了確認這兩個基因?qū)ζ咸烟欠肿拥霓D(zhuǎn)運作用,實驗組利用兩次重組技術(shù),逐次刪除了abc1和abc2基因,得到工程菌株CGΔabc1和CGΔabc2。葡萄糖的消耗速率和菌體的生長速率都大幅度下降,首次發(fā)現(xiàn)abc1和abc2基因編碼的糖類蛋白對葡萄糖的轉(zhuǎn)運起到重要作用。

    Ikeda等[1]發(fā)現(xiàn)C. glutamicum中編碼肌醇轉(zhuǎn)運的兩個基因iolt1和iolt2可能介導葡萄糖的攝取,負責編碼葡萄糖通透酶。鑒于本研究得到的4 個葡萄糖轉(zhuǎn)運相關(guān)基因缺失的突變株仍具有葡萄糖代謝功能,本研究參考Ikeda等[1]的研究,在CGΔabc2突變株基礎(chǔ)上,敲除了iolt1基因,得到的工程菌株CGΔiolt1在葡萄糖為唯一碳源的培養(yǎng)基中不再生長,表明C. glutamicum ATCC 13032的葡萄糖代謝途徑已經(jīng)被完全阻斷,說明在基因組中標注為編碼肌醇轉(zhuǎn)運的蛋白質(zhì),同樣負責葡萄糖的轉(zhuǎn)運,與其報道一致,但iolt2對葡糖糖轉(zhuǎn)運的作用還有待進一步的研究。

    [1] IKEDA M, TAKENO S. Amino acid production by Corynebacterium glutamicum[M]//Corynebacterium glutamicum. Berlin: Springer, 2013: 107-147. DOI:10.1007/978-3-642-29857-8_4.

    [2] BLOMBACH B, SEIBOLD G M. Carbohydrate metabolism in Corynebacterium glutamicum and applications for the metabolic engineering of L-lysine production strains[J]. Applied Microbiology and Biotechnology, 2010, 86(5): 1313-1322. DOI:10.1007/s00253-010-2537-z.

    [3] HOU Xiaohu, CHEN Xinde, ZHANG Yue, et al. L-Valine production with minimization of by-products’ synthesis in Corynebacterium glutamicum and Brevibacterium favum[J]. Amino Acids, 2012, 43(6): 2301-2311. DOI:10.1007/s00726-012-1308-9.

    [4] XU J, ZHANG J, GUO Y, et al. Improvement of cell growth and L-lysine production by genetically modified Corynebacterium glutamicum during growth on molasses[J]. Journal of Industrial Microbiology and Biotechnology, 2013, 40(12): 1423-1432. DOI:10.1007/s10295-013-1329-8.

    [5] KIEFER P, HEINZLE E, ZELDER O, et al. Comparative metabolic flux analysis of lysine-producing Corynebacterium glutamicum cultured on glucose or fructose[J]. Applied and Environmental Microbiology, 2004, 70(1): 229-239. DOI:10.1128/AEM.70.1.229-239.2004.

    [6] KALINOWSKI J, BATHE B, BARTELS D, et al. The complete Corynebacterium glutamicum ATCC 13032 genome sequence and its impact on the production of L-aspartate-derived amino acids and vitamins[J]. Journal of Biotechnology, 2003, 104: 5-25. DOI:10.1016/ S0168-1656(03)00154-8.

    [7] MENTZ A, NESHAT A, PFEIFER-SANCAR K, et al. Comprehensive discovery and characterization of small RNAs in Corynebacterium glutamicum ATCC 13032[J]. BMC genomics, 2013, 14: 714. DOI:10.1186/1471-2164-14-714.

    [8] WENDISCH V F, POLEN T. Transcriptome/proteome analysis of Corynebacterium glutamicum[M]//Corynebacterium glutamicum. Berlin: Springer, 2013: 173-216. DOI:10.1007/978-3-642-29857-8_6.

    [9] 羅玉常, 竇文芳, 張曉梅, 等. 谷氨酸棒桿菌ilvE基因的敲除對相關(guān)氨基酸合成的影響[J]. 生物技術(shù)通報, 2012(11): 185-191.

    [10] LU Dongmei, LIU Jianzhong, MAO Zongwan. Engineering of Corynebacterium glutamicum to enhance L-ornithine production by gene knockout and comparative proteomic analysis[J]. Chinese Journal of Chemical Engineering, 2012, 20(4): 731-739. DOI:10.1016/S1004-9541(11)60242-5.

    [11] TERAM OTO H, INUI M. Regulation of sugar uptake, glycolysis, and the pentose phosphate pathway in Corynebacterium glutamicum[M]// Corynebacterium glutamicum. Berlin: Springer, 2013: 263-279. DOI:10.1007/978-3-642-29857-8_9.

    [12] PARCHE S, BURKOVSKI A, SPRENGER G A, et al. Corynebacterium glutamicum: a dissection of the PTS[J]. Journal of Molecular Microbiology and Biotechnology, 2001, 3(3): 423-428.

    [13] SHI F, LI K, HUAN X J, et al. Expression of NAD(H) kinase and glucose-6-phosphate dehydrogenase improve NADPH supply and L-isoleucine biosynthesis in Corynebacterium glutamicum ssp. lactofermentum[J]. Applied Biochemistry and Biotechnology, 2013, 171(2): 504-521. DOI:10.1007/s12010-013-0389-6.

    [14] PARK S Y, KIM H K, YOO S K, et al. Characterization of glk, a gene coding for glucose kinase of Corynebacterium glutamicum[J]. FEMS Microbiology Letters, 2000, 188(2): 209-215. DOI:10.1111/j.1574-6968.2000.tb09195.x.

    [15] JANG K, BRITZ M. Improved electrotransformation frequencies of Corynebacterium glutamicum using cell-surface mutants[J]. Biotechnology Letters, 2000, 22(7): 539-545. DOI:10.1023/ A:1005629224109.

    [16] KRINGS E, KRUMBACH K, BATHE B, et al. Characterization of myo-inositol utilization by Corynebacterium glutamicum: the stimulon, identification of transporters, and influence on L-lysine formation[J]. Journal of Bacteriology, 2006, 188(23): 8054-8061. DOI:10.1128/ JB.00935-06.

    [17] MOON M W, PARK S Y, CJOI S K, et al. The phosphotransferase system of Corynebacterium glutamicum: features of sugar transport and carbon regulation[J]. Journal of Molecular Microbiology and Biotechnology, 2007, 12(1/2): 43-50. DOI:10.1159/000096458.

    [18] ZHOU Z H, WANG C, XU H T, et al. Increasing succinic acid production using the PTS-independent glucose transport system in a Corynebacterium glutamicum PTS-defective mutant[J]. Journal of Industrial Microbiology and Biotechnology, 2015, 42(7): 1073-1082. DOI:10.1159/S10295-015-1630-9.

    [19] DOMINGUEZ H, ROLLIN C, GUYONVARCH A, et al. Carbon-flux distribution in the central metabolic pathways of Corynebacterium glutamicum during growth on fructose[J]. European Journal of Biochemistry, 1998, 254(1): 96-102. DOI:10.1046/j.1432-1327.

    [20] WITTMANN C, KIEFER P, ZELDER O. Metabolic fluxes in Corynebacterium glutamicum during lysine production with sucrose as carbon source[J]. Applied and Environmental Microbiology, 2004, 70(12): 7277-7287. DOI:10.1128/AEM.70.12.7277-7287.2004.

    [21] IKEDA M. Sugar transport systems in Corynebacterium glutamicum: features and applications to strain development[J]. Applied Microbiology and Biotechnology, 2012, 96(5): 1191-1200. DOI:10.1007/S00253-012-4488-z.

    [22] KOTRBA P, INUI M, YUKAWA H. A single V317A or V317M substitution in enzyme Ⅱ of a newly identified β-glucoside phosphotransferase and utilization system of Corynebacterium glutamicum R extends its specificity towards cellobiose[J]. Microbiology, 2003, 149: 1569-1580. DOI:10.1099/mic.0.26053-0.

    [23] LETEK M, VALBUENA N, RAMOS A, et al. Characterization and use of catabolite-repressed promoters from gluconate genes in Corynebacterium glutamicum[J]. Journal of Bacteriology, 2006, 188(2): 409-423. DOI:10.1128/JB.188.2.409-423.2006.

    [24] NENTWICH S S, BRINKROLF K, GAIGALAT L, et al. Characterization of the LacI-type transcriptional repressor RbsR controlling ribose transport in Corynebacterium glutamicum ATCC 13032[J]. Microbiology, 2009, 155(1): 150-164. DOI:10.1099/ mic.0.020388-0.

    [25] BRINKROLF K, PL?GER S, SOLLE S, et al. The LacI/GalR family transcriptional regulator UriR negatively controls uridine utilization of Corynebacterium glutamicum by binding to catabolite-responsive element (CRE)-like sequences[J]. Microbiology, 2008, 154: 1068-1081. DOI:10.1099/mic.0.2007/014001-0.

    Construction of Engineered Strain Blocking Glucose Metabolism in Corynebacterium glutamicum

    HAN Wuyang1, LIU Jinlei1, DU Hongyan1, WANG Beichen2, YI Hong1, LI Tianming1,*
    (1. College of Biological Science and Engineering, Hebei University of Science and Technology, Shijiazhuang 050018, China; 2. College of Agriculture and Life Science, University of Wisconsin, Madison 53706, USA)

    In this study, we used the wild-type Corynebacterium glutamicum strain ATCC 13032 as a starting strain to obtain f ve mutants with knockout of the ptsG and ptsH-ptsI genes as well as the abc, abc2 and iolt1 genes by homologous recombination without resistance marker using reverse metabolic engineering strategies. Our experimental results showed that compared to the wild-type strain, the glucose consumption of the mutant CGΔptsG was 50% using glucose as the sole carbon source, giving an OD value of 1.473, the glucose consumption of the mutant CGΔptsH-ptsI was 39.5%, giving an OD value of 1.226, the glucose consumption of the mutant CGΔabc was 36%, giving an OD value of 1.09, and the glucose consumption of CGΔabc2 was 26.2%, giving an OD value of 0.486, while the mutant CGΔiolt1 could not utilize glucose to grow, suggesting that glucose metabolism of C. glutamicum was blocked. It turned out that the glucose transporter function was controlled by the ptsG, ptsH-ptsI, abc, abc2 and iolt1 genes, encoding transporter proteins.

    Corynebacterium glutamicum; glucose; gene knockout

    10.7506/spkx1002-6630-201702011

    Q789

    A

    1002-6630(2017)02-0065-10

    韓武洋, 劉金雷, 杜紅燕, 等. 谷氨酸棒狀桿菌葡萄糖代謝阻斷工程菌的構(gòu)建[J]. 食品科學, 2017, 38(2): 65-74. DOI:10.7506/spkx1002-6630-201702011. http://www.spkx.net.cn

    HAN Wuyang, LIU Jinlei, DU Hongyan, et al. Construction of engineered strain blocking glucose metabolism in Corynebacterium glutamicum[J]. Food Science, 2017, 38(2): 65-74. (in Chinese with English abstract)

    10.7506/ spkx1002-6630-201702011. http://www.spkx.net.cn

    2016-03-28

    國家高技術(shù)研究發(fā)展計劃(863計劃)項目(2014AA022102)

    韓武洋(1991—),男,碩士研究生,研究方向為合成生物學與代謝工程。E-mail:924084885@qq.com

    *通信作者:李天明(1986—),男,助理研究員,碩士,研究方向為合成生物學與代謝工程。E-mail:iamltm2000@hotmail.com

    猜你喜歡
    碳源基因組質(zhì)粒
    緩釋碳源促進生物反硝化脫氮技術(shù)研究進展
    牛參考基因組中發(fā)現(xiàn)被忽視基因
    不同碳源對銅溜槽用鋁碳質(zhì)涂抹料性能的影響
    昆鋼科技(2021年6期)2021-03-09 06:10:20
    短乳桿菌天然質(zhì)粒分類
    食品科學(2018年10期)2018-05-23 01:27:28
    四甘醇作碳源合成Li3V2(PO4)3正極材料及其電化學性能
    重組質(zhì)粒rAd-EGF構(gòu)建并轉(zhuǎn)染hDPSCs對其增殖的影響
    Survivin-siRNA重組質(zhì)粒對人神經(jīng)母細胞瘤細胞SH-SY5Y的作用
    BAG4過表達重組質(zhì)粒構(gòu)建及轉(zhuǎn)染結(jié)腸癌SW480細胞的鑒定
    基因組DNA甲基化及組蛋白甲基化
    遺傳(2014年3期)2014-02-28 20:58:49
    有趣的植物基因組
    世界科學(2014年8期)2014-02-28 14:58:31
    亚洲精品粉嫩美女一区| 久久精品91无色码中文字幕| 99久久99久久久精品蜜桃| 一级,二级,三级黄色视频| 巨乳人妻的诱惑在线观看| 一级毛片女人18水好多| 一边摸一边抽搐一进一小说 | 国产精品久久久久久人妻精品电影| 欧美亚洲日本最大视频资源| 人妻久久中文字幕网| 免费女性裸体啪啪无遮挡网站| 美女视频免费永久观看网站| 国产激情久久老熟女| 国产成人精品久久二区二区91| 欧美黄色片欧美黄色片| 免费黄频网站在线观看国产| 王馨瑶露胸无遮挡在线观看| 18禁黄网站禁片午夜丰满| 午夜福利在线免费观看网站| 欧美老熟妇乱子伦牲交| 悠悠久久av| 99久久99久久久精品蜜桃| 久久人人爽av亚洲精品天堂| 久久久国产成人精品二区 | 亚洲伊人色综图| 精品少妇一区二区三区视频日本电影| 中文字幕av电影在线播放| 夜夜夜夜夜久久久久| 18禁裸乳无遮挡免费网站照片 | 欧美乱色亚洲激情| 欧美日韩一级在线毛片| 中亚洲国语对白在线视频| a在线观看视频网站| 欧美在线黄色| 一级毛片高清免费大全| 国产av又大| 日本精品一区二区三区蜜桃| 亚洲精品中文字幕在线视频| 国产aⅴ精品一区二区三区波| 亚洲欧美一区二区三区黑人| 久久久国产欧美日韩av| 久久精品国产综合久久久| 一进一出抽搐动态| 一级a爱片免费观看的视频| 日韩免费av在线播放| 亚洲成人手机| 多毛熟女@视频| 久久精品国产99精品国产亚洲性色 | 亚洲中文日韩欧美视频| 自拍欧美九色日韩亚洲蝌蚪91| 一区在线观看完整版| 天天添夜夜摸| 国产一区有黄有色的免费视频| 欧美日韩黄片免| 欧美激情久久久久久爽电影 | a级毛片黄视频| 欧美另类亚洲清纯唯美| 精品福利永久在线观看| 欧美老熟妇乱子伦牲交| 久久婷婷成人综合色麻豆| 看免费av毛片| 很黄的视频免费| 一边摸一边抽搐一进一小说 | 国产亚洲欧美精品永久| 一进一出抽搐动态| 交换朋友夫妻互换小说| 国产一区二区三区在线臀色熟女 | 亚洲第一av免费看| a级毛片黄视频| 久久ye,这里只有精品| 老司机午夜福利在线观看视频| tocl精华| 1024视频免费在线观看| 黄色视频,在线免费观看| 欧美+亚洲+日韩+国产| 黑人欧美特级aaaaaa片| 日本vs欧美在线观看视频| 国产免费现黄频在线看| 波多野结衣av一区二区av| 色婷婷久久久亚洲欧美| 国产人伦9x9x在线观看| 欧美日韩视频精品一区| 久久久国产精品麻豆| 无遮挡黄片免费观看| 色老头精品视频在线观看| 国产免费男女视频| 欧美丝袜亚洲另类 | 国产亚洲精品久久久久久毛片 | 午夜视频精品福利| 啦啦啦视频在线资源免费观看| 99久久综合精品五月天人人| 成年人黄色毛片网站| 满18在线观看网站| 欧美日韩亚洲综合一区二区三区_| 中文字幕人妻丝袜一区二区| 婷婷精品国产亚洲av在线 | 人人妻人人澡人人爽人人夜夜| 亚洲欧美色中文字幕在线| 三上悠亚av全集在线观看| 国产精品一区二区免费欧美| bbb黄色大片| 日韩 欧美 亚洲 中文字幕| 国产精品乱码一区二三区的特点 | 欧美精品啪啪一区二区三区| 国产高清videossex| 欧美日韩一级在线毛片| 欧美激情高清一区二区三区| 在线永久观看黄色视频| 精品人妻1区二区| 久久久精品区二区三区| 中出人妻视频一区二区| 嫁个100分男人电影在线观看| 久久久久精品人妻al黑| 精品视频人人做人人爽| 久久国产精品男人的天堂亚洲| 涩涩av久久男人的天堂| 国产成人系列免费观看| 丝袜在线中文字幕| 国产精品1区2区在线观看. | 亚洲欧美精品综合一区二区三区| 99久久99久久久精品蜜桃| 又黄又粗又硬又大视频| 91九色精品人成在线观看| 在线观看舔阴道视频| 欧美性长视频在线观看| 欧美性长视频在线观看| 王馨瑶露胸无遮挡在线观看| 9191精品国产免费久久| 亚洲成a人片在线一区二区| 午夜福利在线观看吧| 啦啦啦 在线观看视频| 国产亚洲精品第一综合不卡| 免费在线观看影片大全网站| 国产精品久久久久成人av| 黄色片一级片一级黄色片| 国产日韩欧美亚洲二区| 操美女的视频在线观看| 一边摸一边抽搐一进一小说 | 麻豆av在线久日| 亚洲精品美女久久av网站| 欧美黑人精品巨大| 成人av一区二区三区在线看| 亚洲欧美一区二区三区久久| 免费av中文字幕在线| 视频区图区小说| 亚洲一码二码三码区别大吗| 岛国毛片在线播放| 久久精品aⅴ一区二区三区四区| 18禁观看日本| 午夜激情av网站| 国产精品久久久久成人av| 一区二区三区国产精品乱码| 欧美乱码精品一区二区三区| 中文字幕av电影在线播放| 欧美另类亚洲清纯唯美| 成人国语在线视频| 亚洲成国产人片在线观看| 一边摸一边抽搐一进一小说 | 国产精品久久久久久人妻精品电影| 国产欧美日韩一区二区三| 侵犯人妻中文字幕一二三四区| 国产精品电影一区二区三区 | 欧美丝袜亚洲另类 | 久久精品国产亚洲av高清一级| 精品少妇一区二区三区视频日本电影| 亚洲欧美一区二区三区久久| 久久精品国产99精品国产亚洲性色 | 大型av网站在线播放| 日本黄色日本黄色录像| 精品一品国产午夜福利视频| 久久久久久人人人人人| 欧美乱色亚洲激情| 人人妻人人添人人爽欧美一区卜| 人妻久久中文字幕网| 欧美日韩一级在线毛片| 欧美日韩成人在线一区二区| 一区二区三区精品91| 叶爱在线成人免费视频播放| 黄色女人牲交| 午夜免费成人在线视频| 亚洲成人免费电影在线观看| 一本综合久久免费| 国产成人精品久久二区二区免费| 王馨瑶露胸无遮挡在线观看| 久久青草综合色| 精品高清国产在线一区| 午夜成年电影在线免费观看| 首页视频小说图片口味搜索| 亚洲avbb在线观看| 大码成人一级视频| 国产亚洲欧美98| 欧美黄色淫秽网站| 日本wwww免费看| e午夜精品久久久久久久| 啪啪无遮挡十八禁网站| 久久精品成人免费网站| 99国产精品一区二区蜜桃av | 欧美亚洲 丝袜 人妻 在线| 国产伦人伦偷精品视频| 男女下面插进去视频免费观看| 极品教师在线免费播放| 人人妻人人澡人人爽人人夜夜| 亚洲熟妇中文字幕五十中出 | 老司机在亚洲福利影院| 亚洲av成人一区二区三| 亚洲av片天天在线观看| 丝袜美足系列| 精品国产国语对白av| 国产成人精品久久二区二区免费| 国产成人啪精品午夜网站| 香蕉丝袜av| 极品教师在线免费播放| 99国产精品免费福利视频| 日韩制服丝袜自拍偷拍| 婷婷成人精品国产| 老司机靠b影院| 国产精品久久久av美女十八| 国产高清激情床上av| 国产一区有黄有色的免费视频| av中文乱码字幕在线| 女人高潮潮喷娇喘18禁视频| 天天躁狠狠躁夜夜躁狠狠躁| 大型黄色视频在线免费观看| 精品国产超薄肉色丝袜足j| 成年人午夜在线观看视频| 亚洲精品av麻豆狂野| 亚洲欧美日韩高清在线视频| 色尼玛亚洲综合影院| √禁漫天堂资源中文www| 777米奇影视久久| 国产麻豆69| 国产成人一区二区三区免费视频网站| 老司机午夜十八禁免费视频| 99久久国产精品久久久| 久久精品亚洲av国产电影网| 久久久水蜜桃国产精品网| 黄片小视频在线播放| 国产真人三级小视频在线观看| tube8黄色片| 新久久久久国产一级毛片| 成人国产一区最新在线观看| 飞空精品影院首页| 亚洲精品在线美女| 亚洲成人免费电影在线观看| 色老头精品视频在线观看| 色综合婷婷激情| 精品国内亚洲2022精品成人 | 丰满迷人的少妇在线观看| 韩国av一区二区三区四区| 夜夜躁狠狠躁天天躁| 午夜日韩欧美国产| 99re在线观看精品视频| 欧美精品高潮呻吟av久久| 一区二区日韩欧美中文字幕| 中文字幕人妻丝袜一区二区| 欧美乱妇无乱码| 天天躁狠狠躁夜夜躁狠狠躁| 亚洲精品国产精品久久久不卡| 欧美精品人与动牲交sv欧美| 国产区一区二久久| 精品少妇久久久久久888优播| 日韩免费高清中文字幕av| 日韩免费高清中文字幕av| av免费在线观看网站| 手机成人av网站| 99国产精品99久久久久| 国产亚洲一区二区精品| 国产99白浆流出| 午夜精品国产一区二区电影| 亚洲精品成人av观看孕妇| 两性午夜刺激爽爽歪歪视频在线观看 | 99香蕉大伊视频| 精品熟女少妇八av免费久了| 91麻豆av在线| av片东京热男人的天堂| 亚洲成国产人片在线观看| 成人影院久久| 欧美另类亚洲清纯唯美| 高清在线国产一区| 女警被强在线播放| 成人黄色视频免费在线看| 十八禁人妻一区二区| 免费在线观看亚洲国产| 如日韩欧美国产精品一区二区三区| 欧美色视频一区免费| 久久久久久久久久久久大奶| 精品国产乱码久久久久久男人| 亚洲一区二区三区欧美精品| 亚洲五月天丁香| 亚洲熟女精品中文字幕| 美女 人体艺术 gogo| 久久久国产一区二区| 成人国语在线视频| 欧美中文综合在线视频| 成人精品一区二区免费| 高清欧美精品videossex| 成年人午夜在线观看视频| 动漫黄色视频在线观看| 亚洲欧美激情在线| 国产精品久久久av美女十八| 日本欧美视频一区| 精品一区二区三卡| 精品一品国产午夜福利视频| 日韩精品免费视频一区二区三区| 午夜免费观看网址| 亚洲精华国产精华精| 69精品国产乱码久久久| 99热网站在线观看| 高潮久久久久久久久久久不卡| 制服人妻中文乱码| aaaaa片日本免费| 夜夜夜夜夜久久久久| 中文字幕色久视频| 狂野欧美激情性xxxx| 国产精品久久电影中文字幕 | 91成年电影在线观看| 久久午夜亚洲精品久久| 三级毛片av免费| 老司机午夜福利在线观看视频| 校园春色视频在线观看| 50天的宝宝边吃奶边哭怎么回事| 国产激情欧美一区二区| 高清av免费在线| 欧美中文综合在线视频| 在线播放国产精品三级| 黑人欧美特级aaaaaa片| 搡老乐熟女国产| 熟女少妇亚洲综合色aaa.| 亚洲午夜理论影院| 久久精品亚洲熟妇少妇任你| 一本综合久久免费| 黑人操中国人逼视频| 亚洲va日本ⅴa欧美va伊人久久| 身体一侧抽搐| 免费黄频网站在线观看国产| 国产三级黄色录像| 在线观看午夜福利视频| 极品教师在线免费播放| 欧美日韩视频精品一区| 人人妻人人澡人人看| 国产男女内射视频| 国产亚洲精品一区二区www | 侵犯人妻中文字幕一二三四区| 亚洲专区国产一区二区| 成人特级黄色片久久久久久久| 999久久久国产精品视频| 国产亚洲精品第一综合不卡| 天天躁夜夜躁狠狠躁躁| 深夜精品福利| 老司机影院毛片| 国产精品 国内视频| 国产黄色免费在线视频| 嫁个100分男人电影在线观看| 熟女少妇亚洲综合色aaa.| 丰满迷人的少妇在线观看| 国产精品电影一区二区三区 | 12—13女人毛片做爰片一| 亚洲自偷自拍图片 自拍| 下体分泌物呈黄色| 亚洲欧美色中文字幕在线| avwww免费| 一个人免费在线观看的高清视频| 久久久水蜜桃国产精品网| 亚洲av第一区精品v没综合| 亚洲黑人精品在线| 高清在线国产一区| xxxhd国产人妻xxx| 激情在线观看视频在线高清 | videos熟女内射| 80岁老熟妇乱子伦牲交| a级毛片在线看网站| 最新在线观看一区二区三区| 国产成人啪精品午夜网站| xxxhd国产人妻xxx| 两性午夜刺激爽爽歪歪视频在线观看 | 欧美激情 高清一区二区三区| 久久中文字幕一级| 久久香蕉国产精品| 18禁裸乳无遮挡免费网站照片 | 亚洲精品国产一区二区精华液| 一区二区三区国产精品乱码| 亚洲午夜精品一区,二区,三区| 亚洲国产欧美日韩在线播放| 国产一区二区激情短视频| 1024视频免费在线观看| 首页视频小说图片口味搜索| 人人妻人人添人人爽欧美一区卜| 成人av一区二区三区在线看| 80岁老熟妇乱子伦牲交| 一个人免费在线观看的高清视频| av超薄肉色丝袜交足视频| 一区在线观看完整版| 亚洲欧洲精品一区二区精品久久久| 91九色精品人成在线观看| 9191精品国产免费久久| 免费少妇av软件| 人人妻人人爽人人添夜夜欢视频| 高清黄色对白视频在线免费看| 亚洲国产中文字幕在线视频| 久久国产亚洲av麻豆专区| 亚洲午夜精品一区,二区,三区| 免费看a级黄色片| 久久久国产精品麻豆| 黄色视频,在线免费观看| 国产av又大| 校园春色视频在线观看| 欧美国产精品一级二级三级| 老司机靠b影院| 日本黄色日本黄色录像| 黄频高清免费视频| 久久午夜亚洲精品久久| 夫妻午夜视频| 国产成人啪精品午夜网站| 97人妻天天添夜夜摸| 国产一区二区三区视频了| 老司机靠b影院| tube8黄色片| 黄色视频不卡| 俄罗斯特黄特色一大片| 欧美激情 高清一区二区三区| 人人妻人人澡人人爽人人夜夜| 天天操日日干夜夜撸| 欧美乱妇无乱码| 狂野欧美激情性xxxx| 亚洲精品久久午夜乱码| 精品福利观看| 亚洲精品在线观看二区| 欧美老熟妇乱子伦牲交| 叶爱在线成人免费视频播放| 别揉我奶头~嗯~啊~动态视频| 丰满人妻熟妇乱又伦精品不卡| 亚洲国产欧美网| 美女国产高潮福利片在线看| 村上凉子中文字幕在线| 天堂动漫精品| 久久亚洲真实| 成年女人毛片免费观看观看9 | 超色免费av| 老司机福利观看| 精品久久久精品久久久| 亚洲专区字幕在线| 纯流量卡能插随身wifi吗| 精品人妻在线不人妻| 在线天堂中文资源库| 女性生殖器流出的白浆| 久久九九热精品免费| 国产精品永久免费网站| 国产成人精品久久二区二区免费| av免费在线观看网站| 超碰97精品在线观看| 深夜精品福利| 久久精品91无色码中文字幕| 国产精品免费大片| 亚洲自偷自拍图片 自拍| 午夜影院日韩av| 久久精品人人爽人人爽视色| 欧美人与性动交α欧美精品济南到| 最新美女视频免费是黄的| 亚洲av成人一区二区三| 午夜精品国产一区二区电影| 国产精品二区激情视频| a在线观看视频网站| 国产单亲对白刺激| 亚洲精品av麻豆狂野| 亚洲欧美精品综合一区二区三区| 大型黄色视频在线免费观看| 久久中文看片网| 一边摸一边抽搐一进一小说 | 91国产中文字幕| 欧美中文综合在线视频| netflix在线观看网站| 欧美不卡视频在线免费观看 | 精品无人区乱码1区二区| 黑人欧美特级aaaaaa片| 国产在线精品亚洲第一网站| xxxhd国产人妻xxx| 国产精品久久久久成人av| 手机成人av网站| 国产一区在线观看成人免费| 久久久国产精品麻豆| 99香蕉大伊视频| 法律面前人人平等表现在哪些方面| 99国产精品99久久久久| 亚洲五月色婷婷综合| 日本欧美视频一区| 欧美激情高清一区二区三区| 精品福利永久在线观看| 18禁黄网站禁片午夜丰满| 久久国产亚洲av麻豆专区| 欧美成人午夜精品| 精品人妻在线不人妻| 免费观看人在逋| 久久久国产成人免费| 欧美黄色淫秽网站| 美女福利国产在线| 久久精品国产综合久久久| 在线观看免费视频日本深夜| 午夜福利一区二区在线看| 满18在线观看网站| 在线av久久热| 欧美日韩国产mv在线观看视频| 欧美亚洲日本最大视频资源| 91精品三级在线观看| 69精品国产乱码久久久| 夜夜躁狠狠躁天天躁| 999久久久精品免费观看国产| 人妻 亚洲 视频| 天堂中文最新版在线下载| 国产91精品成人一区二区三区| 亚洲黑人精品在线| 国产成人精品久久二区二区91| 狂野欧美激情性xxxx| 欧美日韩国产mv在线观看视频| 亚洲全国av大片| 人人妻人人澡人人爽人人夜夜| 在线观看免费午夜福利视频| 免费在线观看黄色视频的| 国产精品一区二区在线不卡| 热99国产精品久久久久久7| 欧美成人午夜精品| 午夜福利视频在线观看免费| 91成年电影在线观看| 一进一出抽搐动态| 高潮久久久久久久久久久不卡| 在线观看66精品国产| 无遮挡黄片免费观看| 90打野战视频偷拍视频| 久久午夜亚洲精品久久| 欧美中文综合在线视频| 国产精品乱码一区二三区的特点 | 亚洲五月婷婷丁香| 久久久久久久精品吃奶| 中文字幕色久视频| 韩国av一区二区三区四区| 国产精品综合久久久久久久免费 | 人妻 亚洲 视频| 亚洲精华国产精华精| 亚洲第一欧美日韩一区二区三区| netflix在线观看网站| 日韩欧美国产一区二区入口| 黑丝袜美女国产一区| 一区二区三区精品91| 国产人伦9x9x在线观看| 一边摸一边做爽爽视频免费| 国产一区在线观看成人免费| 成人亚洲精品一区在线观看| 欧美日韩视频精品一区| 亚洲欧洲精品一区二区精品久久久| 成年版毛片免费区| 不卡一级毛片| 亚洲熟妇熟女久久| 老司机午夜十八禁免费视频| 午夜91福利影院| 亚洲一区二区三区不卡视频| 精品视频人人做人人爽| 精品欧美一区二区三区在线| 女人高潮潮喷娇喘18禁视频| 后天国语完整版免费观看| 国产精品98久久久久久宅男小说| 少妇裸体淫交视频免费看高清 | 国产亚洲精品一区二区www | 91老司机精品| 精品国产美女av久久久久小说| 久久草成人影院| 精品亚洲成a人片在线观看| 国产成人av教育| 成人av一区二区三区在线看| 国产av又大| 水蜜桃什么品种好| 狠狠婷婷综合久久久久久88av| 欧美精品一区二区免费开放| 黄色丝袜av网址大全| 黄色 视频免费看| 亚洲中文av在线| 啪啪无遮挡十八禁网站| 亚洲国产中文字幕在线视频| 亚洲欧美激情综合另类| 夫妻午夜视频| 国产亚洲av高清不卡| 黑丝袜美女国产一区| 两个人看的免费小视频| 999精品在线视频| 中文字幕人妻丝袜制服| 免费久久久久久久精品成人欧美视频| www日本在线高清视频| 午夜福利,免费看| 久久这里只有精品19| 亚洲国产精品合色在线| 亚洲欧美日韩另类电影网站| 亚洲精品国产区一区二| 天堂√8在线中文| 国产极品粉嫩免费观看在线| xxxhd国产人妻xxx| 自线自在国产av| 最近最新免费中文字幕在线| 国产99久久九九免费精品| 在线视频色国产色| 日韩欧美免费精品| 一级毛片女人18水好多| 一本综合久久免费| 午夜免费观看网址| 欧美丝袜亚洲另类 | 美国免费a级毛片| av线在线观看网站| 69av精品久久久久久| 久久久国产成人精品二区 | 极品教师在线免费播放| 校园春色视频在线观看| 国产亚洲一区二区精品| 国产精品久久久人人做人人爽| 久久这里只有精品19| 欧美日韩精品网址| 国产精品九九99| 亚洲成国产人片在线观看| 搡老熟女国产l中国老女人|