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

    柚皮素合成關(guān)鍵基因表達(dá)水平對(duì)目標(biāo)產(chǎn)物積累水平的量化影響

    2019-08-22 08:55:12焦亭亭周景文徐沙
    生物工程學(xué)報(bào) 2019年7期
    關(guān)鍵詞:柚皮素傳代拷貝數(shù)

    焦亭亭,周景文,徐沙

    柚皮素合成關(guān)鍵基因表達(dá)水平對(duì)目標(biāo)產(chǎn)物積累水平的量化影響

    焦亭亭1,2,周景文1,2,徐沙1,2

    1 江南大學(xué) 糧食發(fā)酵工藝與技術(shù)國(guó)家工程實(shí)驗(yàn)室,江蘇 無(wú)錫 214122 2 江南大學(xué) 生物工程學(xué)院,江蘇 無(wú)錫 214122

    柚皮素是一種天然黃酮類化合物,具有抗炎、抗氧化、抗病毒、預(yù)防動(dòng)脈粥樣硬化等多種藥理活性,也是其他黃酮類化合物合成的重要前體,具有重要的應(yīng)用價(jià)值。目前,微生物法生產(chǎn)柚皮素等黃酮類化合物由于代謝通路不平衡等原因?qū)е庐a(chǎn)量較低,在很大程度上限制了其工業(yè)應(yīng)用。文中以一株產(chǎn)柚皮素的釀酒酵母菌株Y-01為研究對(duì)象,利用啟動(dòng)子和拷貝數(shù)控制柚皮素合成代謝途徑關(guān)鍵酶4-香豆酸:CoA連接酶(4CL)、查爾酮合成酶(CHS) 和查爾酮異構(gòu)酶(CHI) 編碼基因的表達(dá)水平,考察這些基因的表達(dá)水平對(duì)目標(biāo)產(chǎn)物積累水平的量化影響。結(jié)果表明,柚皮素產(chǎn)量與4CL或CHI編碼基因的表達(dá)量之間關(guān)聯(lián)性較低,而與基因的表達(dá)量存在顯著的正相關(guān)性。通過(guò)調(diào)控基因的表達(dá)水平,獲得一株高產(chǎn)柚皮素的釀酒酵母工程菌株Y-04,產(chǎn)量較出發(fā)菌株Y-01提高了4.1倍。研究結(jié)果表明,CHS是柚皮素合成過(guò)程的關(guān)鍵調(diào)控靶點(diǎn),合理調(diào)控CHS表達(dá)可以顯著促進(jìn)釀酒酵母積累柚皮素。相關(guān)結(jié)果為采用代謝工程強(qiáng)化微生物合成柚皮素等重要黃酮類化合物提供了重要的理論參考。

    柚皮素,表達(dá)量,關(guān)鍵基因,啟動(dòng)子,拷貝數(shù)

    柚皮素 (4?,5,7-三羥基二氫黃酮) 是一種重要的二氫黃酮類化合物。柚皮素及其衍生物普遍存在于高等植物中,例如西紅柿和柑橘類水果。藥理研究表明,柚皮素及其衍生物具有抗炎、抗氧化、抗病毒、預(yù)防動(dòng)脈粥樣硬化等各種活性[1]。此外,柚皮素也是大多數(shù)黃酮類化合物的合成前體,可以用于合成具有更多生理活性和應(yīng)用價(jià)值的重要黃酮類化合物。目前,國(guó)內(nèi)外報(bào)道的柚皮素制備方法主要有植物提取法、化學(xué)合成法和微生物法。由于植物提取法存在植物生長(zhǎng)周期長(zhǎng)、含量低、代謝產(chǎn)物復(fù)雜等問(wèn)題,導(dǎo)致成本較高,嚴(yán)重阻礙了柚皮素的研究和應(yīng)用[2-3]?;瘜W(xué)合成法涉及有毒化學(xué)品和極端反應(yīng)條件,除成本較高外,還無(wú)法保證關(guān)鍵手性基團(tuán)的天然等同性[4]。由于微生物法生產(chǎn)黃酮不受季節(jié)限制,具有完全的天然等同性等特點(diǎn),是目前被認(rèn)為最有前景的發(fā)展方向。

    柚皮素的生物合成途徑從芳香族氨基酸開(kāi)始,主要由4種酶參與,即酪氨酸解氨酶 (TAL) /苯丙氨酸解氨酶 (PAL)、4-香豆酸∶CoA連接酶 (4CL)、查爾酮合成酶 (CHS) 和查爾酮異構(gòu)酶 (CHI)。其中,TAL將L-酪氨酸轉(zhuǎn)化為對(duì)香豆酸;4CL催化1分子CoA和1分子對(duì)香豆酸形成對(duì)香豆酰輔酶A;CHS負(fù)責(zé)將3分子丙二酰輔酶A和1分子對(duì)香豆酰輔酶A轉(zhuǎn)化為柚皮素查爾酮;CHI催化柚皮素查爾酮異構(gòu)化為柚皮素。為了提高目標(biāo)產(chǎn)物的產(chǎn)量,研究人員首先關(guān)注于篩選高活性酶和改善酶的性質(zhì)并獲得了一系列具有較高催化活性的酶[5-6]。大量的代謝工程實(shí)例表明,高效的黃酮生物合成途徑還需要調(diào)節(jié)高活性酶的表達(dá)水平。目前,基于不同強(qiáng)度啟動(dòng)子和不同拷貝數(shù)的模塊化優(yōu)化策略已經(jīng)用于平衡代謝流,并顯著強(qiáng)化了生松素[7]和柚皮素[8]等黃酮類物質(zhì)的積累。但是,該策略需要篩選具有梯度強(qiáng)度的啟動(dòng)子文庫(kù),并且在多基因途徑中這些啟動(dòng)子的天文數(shù)目組合在實(shí)際操作中是一個(gè)巨大的挑戰(zhàn)。因此,在構(gòu)建平衡途徑的啟動(dòng)子-關(guān)鍵基因組合之前,探究途徑中目標(biāo)物質(zhì)產(chǎn)量與關(guān)鍵酶編碼基因表達(dá)量的關(guān)系是必要的。為此,文中以一株可以積累柚皮素的釀酒酵母Y-01為研究對(duì)象,利用不同強(qiáng)度啟動(dòng)子和整合位點(diǎn),梯度過(guò)量表達(dá)酶活較高香芹菜4CL、矮牽牛CHS和紫苜宿CHI[8],結(jié)合實(shí)時(shí)熒光定量PCR (RT-PCR) 方法檢測(cè)它們的表達(dá)水平,獲得柚皮素產(chǎn)量與這3個(gè)關(guān)鍵酶表達(dá)量之間的關(guān)系。結(jié)果表明,CHS是黃酮合成過(guò)程的關(guān)鍵調(diào)控靶點(diǎn),合理調(diào)控CHS表達(dá)有利于酵母中柚皮素的積累。研究結(jié)果對(duì)后續(xù)代謝工程強(qiáng)化微生物合成柚皮素等重要黃酮類化合物具有重要的參考價(jià)值。

    1 材料與方法

    1.1 材料

    1.1.1 菌株

    大腸桿菌JM109用于質(zhì)粒構(gòu)建和保存。釀酒酵母Y-01用于整合片段及表達(dá),釀酒酵母S288C用于基因擴(kuò)增。文中使用的主要菌株如表1所示。

    表1 本研究所用的主要菌株

    from,from,from

    1.1.2 培養(yǎng)基

    LB (Luria Broth) 液體培養(yǎng)基 (/L):10 g NaCl,10 g蛋白胨,5 g酵母提取物,pH 7.0。

    YNB (Yeast Nitrogen Base) 培養(yǎng)基 (/L):200 g/L葡萄糖母液100 mL,67.4 g/L YNB (Yeast nitrogen base without amino acids) 母液100 mL,并添加所需的缺陷氨基酸 (Trp、His、Leu) 各10 mL,終濃度均為50 mg/L,加蒸餾水至1 L,過(guò)濾除菌,4 ℃冰箱保存。

    YPSe培養(yǎng)基 (/L):20 g蛋白胨,10 g酵母提取物,10 g蔗糖,10 mL無(wú)水乙醇。

    1.1.3 主要試劑

    各種抗生素和氨基酸,包括氨芐青霉素、遺傳霉素、色氨酸、組氨酸、亮氨酸,均購(gòu)自生工生物工程 (上海) 股份有限公司。柚皮素標(biāo)樣購(gòu)自Sigma-Aldrich公司。蛋白胨、酵母提取物購(gòu)自O(shè)xoid公司。乙腈 (質(zhì)譜純) 購(gòu)自Merck公司。其他本實(shí)驗(yàn)中使用的各種分析化學(xué)試劑如無(wú)特別說(shuō)明均購(gòu)自上海國(guó)藥集團(tuán)。各種一步克隆酶、連接酶、感受態(tài)制備試劑盒、pMD19-T Simple Vector購(gòu)自大連寶生物 (TaKaRa) 有限公司。細(xì)菌質(zhì)粒小量提取試劑盒購(gòu)自于生工生物工程 (上海) 股份有限公司。膠回收試劑盒購(gòu)自Fermentas公司。引物合成以及DNA測(cè)序服務(wù)均由生工生物工程 (上海) 股份有限公司提供。

    1.1.4 搖瓶培養(yǎng)條件

    從平板上挑取改造后產(chǎn)柚皮素釀酒酵母的單菌落,在含有20 mL YNB培養(yǎng)基 (添加Trp和Leu) 的250 mL搖瓶中于30 ℃、220 r/min培養(yǎng)18 h。然后以10% (/) 的接種量轉(zhuǎn)接到20 mL YPSe發(fā)酵培養(yǎng)基中,同時(shí)加入10 g/L CaCO3及終濃度500 mg/L對(duì)香豆酸,于30 ℃、220 r/min培養(yǎng)84 h。

    1.2 方法

    1.2.1 釀酒酵母基因組整合框構(gòu)建

    為了使代謝途徑基因協(xié)同表達(dá),首先確定代謝途徑關(guān)鍵基因過(guò)量表達(dá)對(duì)產(chǎn)物柚皮素含量的影響。選擇6個(gè)不同強(qiáng)度的啟動(dòng)子 (P、P、P、P、P、P) 分別表達(dá)、、基因,并選取了基因組上LEU2和核糖體DNA (rDNA) 作為整合位點(diǎn)。啟動(dòng)子P、P、P、P通過(guò)實(shí)驗(yàn)室原有質(zhì)粒為模板擴(kuò)增獲得,P、P和基因整合的上下游同源臂 (上下游同源臂各500–1 000 bp) 以S288C基因組為模板擴(kuò)增得到,以HIS5作為營(yíng)養(yǎng)缺陷標(biāo)簽。先將上下游同源臂、HIS5標(biāo)簽盒通過(guò)同源重組連接在T載上,在實(shí)驗(yàn)室已有質(zhì)粒pY26-P--P--P-、pY26-P--P-- P-、pY26-P--P--P-上分別擴(kuò)增P--T、P--T和P--T片段,通過(guò)同源重組的方法得到質(zhì)粒T- LEU2up-P--T-HIS5-LEU2down、T-LEU2up- P--T-HIS5-LEU2down、T-LEU2up-P--T-HIS5-LEU2down、T-25Sup-P--T- HIS5-5Sdown、T-25Sup-P--T-HIS5- 5Sdown、T-25Sup-P--T-HIS5-5Sdown,然后分別將另外5個(gè)啟動(dòng)子與這6個(gè)質(zhì)粒中啟動(dòng)子P通過(guò)同源重組的方法替換,得到另外30個(gè)質(zhì)粒。再通過(guò)引物Pf_LEU2和Pr_LEU2分別擴(kuò)增3個(gè)基因在LEU2位點(diǎn)上的整合片段,引物Pf_25S和Pr_5S分別擴(kuò)增3個(gè)基因在rDNA區(qū)域的整合片段。文中所用主要引物見(jiàn)表2。

    采用醋酸鋰轉(zhuǎn)化法將這36個(gè)片段分別整合到菌株Y-01的基因組上,用YNB培養(yǎng)基 (添加Trp和Leu) 篩選,30 ℃恒溫培養(yǎng)3–4 d至長(zhǎng)出單菌落。挑取單菌落至YNB液體篩選培養(yǎng)基中培養(yǎng)24 h,13 500 r/min離心5 min去上清,收集菌體提取釀酒酵母基因組,進(jìn)行PCR驗(yàn)證,得到與實(shí)際大小相符的特異性條帶,則證明菌株構(gòu)建成功。

    1.2.2 基因轉(zhuǎn)錄水平的測(cè)定

    取發(fā)酵液中培養(yǎng)至對(duì)數(shù)中期的釀酒酵母細(xì)胞,采用熱酚法提取總RNA[9]。利用PrimeScript? RT reagent Kit with gDNA Eraser試劑盒 (TaKaRa) 去除總RNA中的基因組DNA后,將總RNA反轉(zhuǎn)錄成cDNA。使用TB Green? Premix Ex? Ⅱ試劑盒處理樣品,在LightCycler 480Ⅱ熱循環(huán)儀系統(tǒng) (Roche,Mannheim,Germany) 上進(jìn)行定量PCR (qPCR) 測(cè)定。以基因作為內(nèi)參基因,用引物Pf_Actin和Pr_Actin進(jìn)行qPCR值的歸一化。、和分別用引物Pf_4CL/Pr_4CL、Pf_CHS/Pr_CHS和Pf_CHI/Pr_CHI進(jìn)行擴(kuò)增。最后采用2?ΔΔCt方法[10]獲得、和的mRNA表達(dá)水平來(lái)表征它們的相對(duì)表達(dá)量。

    表2 本研究中所用的主要引物

    1.2.3 基因整合拷貝數(shù)的測(cè)定

    酵母基因組DNA的提取按照Thermo公司GeneJET Genomic DNA Purification Kit中酵母基因組DNA提取說(shuō)明書(shū)進(jìn)行操作。

    以前的研究表明在不同條件下外源基因在基因組上單拷貝克隆具有遺傳穩(wěn)定性[11-12],而很少有研究能夠確定多拷貝克隆的遺傳穩(wěn)定性。由于一株菌株中多拷貝克隆位點(diǎn)起源于同一基因座處的多個(gè)同源重組事件,外源基因通常在基因組上串聯(lián)重復(fù)可能導(dǎo)致整合基因的不穩(wěn)定[13]。為了確定基因、和整合到酵母rDNA區(qū)域菌株的遺傳穩(wěn)定性,將以rDNA為整合位點(diǎn)的所有釀酒酵母重組菌株在YNB缺陷型平板上劃線,培養(yǎng)3–4 d至長(zhǎng)出單菌落,挑取單菌落至20 mL新鮮的液體YNB缺陷型培養(yǎng)基中,30 ℃、220 r/min培養(yǎng)24 h,然后按1%接種量轉(zhuǎn)移到新鮮的20 mL YPD中并在相同條件下溫育24 h。該過(guò)程再重復(fù)5次,總生長(zhǎng)時(shí)間為144 h,即傳代72代。

    從原始的非傳代菌株和傳代72次菌株中提取基因組DNA,通過(guò)qPCR檢測(cè)基因組中關(guān)鍵酶的整合拷貝數(shù)。對(duì)帶有目的基因、和的T載體分別進(jìn)行梯度稀釋,制備目的基因質(zhì)粒的標(biāo)準(zhǔn)曲線,取基因組DNA各1 μL為模板,用熒光引物對(duì)目的基因進(jìn)行qPCR,將得到的t值代入標(biāo)準(zhǔn)曲線中,求出DNA樣品中起始模板的拷貝數(shù)。利用相同的方法計(jì)算看家基因(編碼葡聚糖1,4-α-葡糖苷酶,在基因組以單拷貝的形式存在) 的起始模板的拷貝數(shù)。因?yàn)樵诒磉_(dá)菌株中都導(dǎo)入了質(zhì)粒pY26-P--P--P-,因此要除去質(zhì)粒中目的基因的拷貝數(shù)。利用上述同樣的方法制作該質(zhì)粒的標(biāo)準(zhǔn)曲線,以質(zhì)粒中作為目標(biāo)基因。因此,基因組中關(guān)鍵酶的整合拷貝數(shù)計(jì)算公式為:N=(N1?N2)/N0+1(N1代表基因組中目的基因的起始模板的拷貝數(shù),N2代表基因組中基因起始模板的拷貝數(shù),N0代表基因組中看家基因起始模板的拷貝數(shù),1代表基因組染色體上原有一個(gè)拷貝數(shù))。

    1.2.4 柚皮素的提取與檢測(cè)

    發(fā)酵結(jié)束后,收集發(fā)酵液并加入同體積的無(wú)水乙醇萃取柚皮素,劇烈振蕩混勻后,13 500 r/min離心5 min取上清,有機(jī)相濾膜過(guò)濾后使用Agilent 1260高效液相色譜儀進(jìn)行產(chǎn)物檢測(cè),采用C18色譜柱 (4.6 mm×250 mm) 進(jìn)行色譜分離,色譜柱溫度控制在25 ℃,進(jìn)樣量為10 μL,流動(dòng)相為超純水 (0.3%乙酸) 和乙腈 (0.3%乙酸),流速為1 mL/min。柚皮素通過(guò)乙腈/水的梯度洗脫進(jìn)行分離:0–10 min,色譜條件為10%–40%的乙腈 (體積/體積);10–15 min,色譜條件控制為40%–60%的乙腈;15–17 min,色譜條件控制在60%–10%的乙腈[14]。柚皮素的檢測(cè)波長(zhǎng)為290 nm,以標(biāo)準(zhǔn)品柚皮素濃度為橫坐標(biāo),紫外區(qū)吸收峰峰面積為縱坐標(biāo),繪制標(biāo)準(zhǔn)曲線。

    2 結(jié)果與分析

    2.1 4CL過(guò)量表達(dá)對(duì)柚皮素產(chǎn)量的影響

    以Y-01為出發(fā)菌株,將在不同強(qiáng)度啟動(dòng)子 (P、P、P、P、P、P) 控制下的基因分別整合在基因組LEU2位點(diǎn)和rDNA區(qū)域。每株菌在發(fā)酵培養(yǎng)基中發(fā)酵得到柚皮素的產(chǎn)量,qPCR得到轉(zhuǎn)錄水平上的相對(duì)表達(dá)量,從而獲得柚皮素產(chǎn)量與轉(zhuǎn)錄水平的關(guān)系。

    每株菌做3個(gè)平行,并重復(fù)3次,以確保實(shí)驗(yàn)數(shù)據(jù)的可靠性,從而獲得36組數(shù)據(jù)。結(jié)果如圖1所示,柚皮素的產(chǎn)量與的相對(duì)表達(dá)量幾乎沒(méi)有相關(guān)性 (=?0.2,=1.62×10–3),并且在基因組上單拷貝整合時(shí)有利于柚皮素的積累。以LEU2為整合位點(diǎn),在啟動(dòng)子P控制下,柚皮素的最高產(chǎn)量為165.9 mg/L (Y-02),比Y-01 (46.1 mg/L) 提高了2.6倍;在啟動(dòng)子P控制下,獲得最低產(chǎn)量81.1 mg/L (Y-03),比Y-01增加了75.9%。以rDNA區(qū)域?yàn)檎衔稽c(diǎn)過(guò)量表達(dá),柚皮素產(chǎn)量在21.5 mg/L到97.0 mg/L之間變化,最低產(chǎn)量 (Y-12) 比Y-01減少了一半。結(jié)果表明過(guò)量表達(dá)并不一定能提高柚皮素的產(chǎn)量。

    圖1 柚皮素產(chǎn)量與4cl轉(zhuǎn)錄水平的關(guān)系

    2.2 CHS過(guò)量表達(dá)對(duì)柚皮素產(chǎn)量的影響

    以上述相同的方法分析CHS過(guò)量表達(dá)對(duì)柚皮素產(chǎn)量的影響。結(jié)果如圖2所示,柚皮素的產(chǎn)量與轉(zhuǎn)錄水平有一定的線性正相關(guān)性 (=0.8,=6.24×10–10),并且單拷貝數(shù)克隆時(shí)更有利于柚皮素的積累。以LEU2為整合位點(diǎn)時(shí),在強(qiáng)啟動(dòng)子 (P、P、P) 控制下過(guò)量表達(dá),轉(zhuǎn)錄水平和柚皮素產(chǎn)量均較高,轉(zhuǎn)錄水平在0.8到1.0之間,柚皮素產(chǎn)量在185.3 mg/L到235.9 mg/L之間。當(dāng)受P啟動(dòng)子控制時(shí),柚皮素最高產(chǎn)量為235.9 mg/L (Y-04),較Y-01提高了4.1倍;當(dāng)受P啟動(dòng)子控制時(shí),最低產(chǎn)量為114.0 mg/L (Y-05),較Y-01的產(chǎn)量提高了1.5倍。以rDNA為整合位點(diǎn)過(guò)量表達(dá),柚皮素產(chǎn)量在36.5–174.9 mg/L之間變化,最低產(chǎn)量 (Y-14) 比Y-01降低了20.8%。結(jié)果表明,在一定程度上過(guò)量表達(dá)更有利于柚皮素的大量積累。

    圖2 柚皮素產(chǎn)量與chs轉(zhuǎn)錄水平的關(guān)系

    2.3 CHI過(guò)量表達(dá)對(duì)柚皮素產(chǎn)量的影響

    以上述相同的方法分析CHI過(guò)量表達(dá)對(duì)柚皮素產(chǎn)量的影響。結(jié)果如圖3所示,柚皮素的產(chǎn)量與轉(zhuǎn)錄水平的線性關(guān)系較小 (=0.4,=3.14×10–5),并且CHI在基因組單拷貝克隆時(shí)更利于柚皮素的積累。

    圖3 柚皮素產(chǎn)量與chi轉(zhuǎn)錄水平的關(guān)系

    當(dāng)LEU2作整合位點(diǎn)時(shí),在啟動(dòng)子P控制下,柚皮素最高產(chǎn)量達(dá)到168.5 mg/L (Y-06),比Y-01提高了2.7倍;在啟動(dòng)子P控制下,柚皮素最低產(chǎn)量為115.0 mg/L (Y-07),比Y-01提高了1.5倍。以rDNA為整合位點(diǎn)過(guò)量表達(dá),柚皮素的產(chǎn)量在81.9 mg/L到128.2 mg/L之間變化,最低產(chǎn)量 (Y-20) 比Y-01增加了77.7%。結(jié)果表明,過(guò)量表達(dá)提高了柚皮素的產(chǎn)量。

    2.4 多拷貝克隆菌株傳代穩(wěn)定性的驗(yàn)證

    結(jié)果如圖4所示,基因組整合重組菌Y-08、Y-09、Y-10、Y-11、Y-12及Y-13未傳代與傳代72代后拷貝數(shù)差異分別為9.0%、5.9%、23.6%、9.5%、10.1%、15.7% (圖4A);基因組整合重組菌Y-14、Y-15、Y-16、Y-17、Y-18及Y-19未傳代與傳代72代后拷貝數(shù)差異分別為12.5%、1.3%、0.5%、5.8%、1.7%、2.2% (圖4B);基因組整合重組菌Y-20、Y-21、Y-22、Y-23、Y-24及Y-25未傳代與傳代72代后拷貝數(shù)差異分別為6.4%、5.1%、7.9%、11.6%、5.7%、10.6% (圖4C)。上述結(jié)果表明,多拷貝克隆菌株傳代72代后拷貝數(shù)差異基本在15%以下,表明了在本實(shí)驗(yàn)條件下多拷貝克隆菌株具有遺傳穩(wěn)定性。

    3 討論

    文中以產(chǎn)柚皮素釀酒酵母Y-01為出發(fā)菌株,首次考察了關(guān)鍵酶4CL、CHS和CHI分別過(guò)量表達(dá)對(duì)柚皮素產(chǎn)量的影響。當(dāng)在強(qiáng)啟動(dòng)子P控制下并插入到LEU2位點(diǎn)時(shí),達(dá)到最高產(chǎn)量235.9 mg/L,比目前報(bào)道的以釀酒酵母為宿主異源合成柚皮素最高產(chǎn)量113.0 mg/L[15]提高了1.1倍。在以前的研究中,因柚皮素合成需要3分子的丙二酰輔酶A,而胞內(nèi)游離丙二酰輔酶A濃度由于脂肪酸生物合成的內(nèi)在緊密調(diào)節(jié)維持在低水平,使丙二酰輔酶A一直被認(rèn)為是黃酮合成的限制性前體因子。一般采用過(guò)量表達(dá)乙酰輔酶A羧化酶 (Acc1)[16]、敲除丙二酰輔酶A的競(jìng)爭(zhēng)途徑[17]和基于反義RNA策略[18]及CRISPRi系統(tǒng)調(diào)控[19]等方法提高胞內(nèi)丙二酰輔酶A含量,使細(xì)胞內(nèi)丙二酰輔酶A最大程度地流向目標(biāo)途徑以提高目標(biāo)化合物產(chǎn)量。但是,最新報(bào)道顯示在釀酒酵母中異源合成柚皮素的前期階段,丙二酰輔酶A并不是柚皮素合成的限速前體,關(guān)鍵酶之間的協(xié)同表達(dá)更有利于柚皮素的積累[20]。對(duì)于調(diào)控途徑基因的協(xié)同表達(dá),常用的方法是過(guò)量表達(dá)關(guān)鍵基因的同時(shí)調(diào)節(jié)其表達(dá)量,包括調(diào)節(jié)基因拷貝數(shù)、啟動(dòng)子強(qiáng)度及降低蛋白質(zhì)的降解等方式[15,21]。但這些方式所調(diào)控的范圍往往較大,目前還無(wú)法對(duì)蛋白表達(dá)進(jìn)行精確調(diào)控。近年來(lái),越來(lái)越多的計(jì)算機(jī)輔助方法應(yīng)用于代謝通量平衡分析,如OptForce[22]、OptORF[23]和RobustKnock[24]等。許多基于計(jì)算機(jī)的途徑設(shè)計(jì)策略已成功應(yīng)用于所需化合物的生產(chǎn),例如三乙酸內(nèi)酯[25]、琥珀酸[26]、脂肪酸[27]等。然而,由于代謝途徑涉及眾多的前體合成和相關(guān)修飾基因?qū)е氯狈?zhǔn)確的酶反應(yīng)參數(shù),使這些調(diào)節(jié)策略的使用受到極大的限制。

    目前,精確調(diào)控途徑基因之間的協(xié)同表達(dá)以平衡代謝通路還是一個(gè)巨大的難題。清楚地了解關(guān)鍵酶表達(dá)對(duì)產(chǎn)物柚皮素的影響有利于后期平衡代謝流的研究。據(jù)報(bào)道,CHS是燕麥葉中黃酮類化合物的限速酶[28],植物可以通過(guò)調(diào)控表達(dá)來(lái)調(diào)節(jié)黃酮類化合物的產(chǎn)生[29-30]。而以微生物為生產(chǎn)宿主異源合成柚皮素途徑尚無(wú)明確的報(bào)道。文中提出了啟動(dòng)子和拷貝數(shù)控制代謝工程,并采用qPCR方法檢測(cè)細(xì)胞中靶基因的轉(zhuǎn)錄水平,探究了關(guān)鍵酶4CL、CHS和CHI分別過(guò)量表達(dá)對(duì)柚皮素產(chǎn)量的影響。結(jié)果顯示在一定范圍內(nèi)分別過(guò)量表達(dá)4CL、CHS和CHI均能提高柚皮素的產(chǎn)量,但是柚皮素產(chǎn)量與4CL或CHI的表達(dá)量沒(méi)有密切的聯(lián)系,而與CHS的表達(dá)量存在線性正相關(guān) (=0.8,=6.24×10–10),并且柚皮素的產(chǎn)量在過(guò)量表達(dá)CHS時(shí)達(dá)到最大值,這表明CHS表達(dá)對(duì)柚皮素產(chǎn)量的影響最大,并推測(cè)CHS是該研究中柚皮素合成途徑的限速酶。

    此外,實(shí)驗(yàn)結(jié)果顯示,基因單拷貝克隆的表達(dá)水平大多高于多拷貝克隆的表達(dá)水平,造成的原因可能與釀酒酵母rDNA區(qū)域的基因沉默有關(guān)[31],外源基因插入該區(qū)域,其表達(dá)水平往往低于插入到常染色質(zhì)區(qū)域的表達(dá)水平。rDNA作為外源基因整合位點(diǎn),插入其內(nèi)部的外源基因的遺傳穩(wěn)定性還沒(méi)有明確的結(jié)論,這可能與外源基因的拷貝數(shù)和整合形式有關(guān)[32]。文中通過(guò)qPCR檢測(cè)未傳代和傳代后菌株中外源基因拷貝數(shù)并比較兩者的差異,發(fā)現(xiàn)在實(shí)驗(yàn)條件下外源基因插入rDNA區(qū)域具有遺傳穩(wěn)定性。文中確立了柚皮素產(chǎn)量與關(guān)鍵酶編碼基因相對(duì)表達(dá)量的關(guān)系,并揭示了CHS是黃酮化合物合成過(guò)程的關(guān)鍵調(diào)控位點(diǎn),為后續(xù)深入研究代謝工程強(qiáng)化微生物合成柚皮素等重要黃酮類化合物提供了重要的參考依據(jù)。

    [1] Iwashina T. Flavonoid function and activity to plants and other organisms. Biol Sci Space, 2003, 17(1): 24–44.

    [2] Fowler ZL, Koffas MA. Biosynthesis and biotechnological production of flavanones: Current state and perspectives. Appl Microbiol Biotechnol, 2009, 83(5): 799–808.

    [3] Wang YC, Chen S, Yu O. Metabolic engineering of flavonoids in plants and microorganisms. Appl Microbiol Biotechnol, 2011, 91(4): 949–956.

    [4] Xia YQ, Guo TY, Zhao HL, et al. A novel solid phase for selective separation of flavonoid compounds. J Sep Sci, 2007, 30(9): 1300–1306.

    [5] Guo L, Chen X, Li LN, et al. Transcriptome-enabled discovery and functional characterization of enzymes related to ()-pinocembrin biosynthesis fromand their application for metabolic engineering. Microb Cell Fact, 2016, 15: 27.

    [6] Zhou SH, Liu PR, Chen J, et al. Characterization of mutants of a tyrosine ammonia-lyase from. Appl Microbiol Biotechnol, 2016, 100(24): 10443–10452.

    [7] Wu JJ, Du GC, Zhou JW, et al. Metabolic engineering offor (2)-pinocembrin production from glucose by a modular metabolic strategy. Metab Eng, 2013, 16: 48–55.

    [8] Wu JJ, Zhou TT, Du GC, et al. Modular optimization of heterologous pathways forsynthesis of (2S)-naringenin in. PLoS ONE, 2014, 9(7): e101492.

    [9] Schmitt ME, Brown TA, Trumpower BL. A rapid and simple method for preparation of RNA from. Nucleic Acids Res, 1990, 18(10): 3091–3092.

    [10] Livak KJ, Schmittgen TD. Analysis of relative gene expression data using real-time quantitative PCR and the 2?ΔΔCtmethod. Methods, 2001, 25(4): 402–408.

    [11] Fukunaga T, Cha-Aim K, Hirakawa Y, et al. Designed construction of recombinant DNA at thelocus in the yeast. Yeast, 2013, 30(6): 243–253.

    [12] Ohi H, Okazaki N, Uno S, et al. Chromosomal DNA patterns and gene stability of. Yeast, 1998, 14(10): 895–903.

    [13] Zhu J, Schiestl RH. Human topoisomerase I mediates illegitimate recombination leading to DNA insertion into the ribosomal DNA locus in. Mol Genet Genomics, 2004, 271(3): 347–358.

    [14] Wu JJ. Engineeringbased on synthetic biology for production of flavonoid scaffolds [D]. Wuxi: Jiangnan University, 2015 (in Chinese).吳俊俊. 合成生物技術(shù)改造大腸桿菌生產(chǎn)黃酮骨架物質(zhì)[D]. 無(wú)錫: 江南大學(xué), 2015.

    [15] Koopman F, Beekwilder J, Crimi B, et al.production of the flavonoid naringenin in engineered. Microb Cell Fact, 2012, 11: 155.

    [16] Shin SY, Jung SM, Kim MD, et al. Production of resveratrol from tyrosine in metabolically engineered. Enzyme Microb Technol, 2012, 51(4): 211–216.

    [17] Xu P, Ranganathan S, Fowler ZL, et al. Genome-scale metabolic network modeling results in minimal interventions that cooperatively force carbon flux towards malonyl-CoA. Metab Eng, 2011, 13(5): 578–587.

    [18] Wu JJ, Yu O, Du GC, et al. Fine-tuning of the fatty acid pathway by synthetic antisense RNA for enhanced (2)-naringenin production from L-tyrosine in. Appl Environ Microbiol, 2014, 80(23): 7283–7292.

    [19] Wu JJ, Du GC, Chen J, et al. Enhancing flavonoid production by systematically tuning the central metabolic pathways based on a CRISPR interference system in. Sci Rep, 2015, 5: 13477.

    [20] Lyu XM, Ng KR, Lee JL, et al. Enhancement of naringenin biosynthesis from tyrosine by metabolic engineering of. J Agric Food Chem, 2017, 65(31): 6638–6646.

    [21] Li MJ, Kildegaard KR, Chen Y, et al.production of resveratrol from glucose or ethanol by engineered. Metab Eng, 2015, 32: 1–11.

    [22] Ranganathan S, Suthers PF, Maranas CD. OptForce: An optimization procedure for identifying all genetic manipulations leading to targeted overproductions. PLoS Comput Biol, 2010, 6(4): e1000744.

    [23] Kim J, Reed JL. Optorf: Optimal metabolic and regulatory perturbations for metabolic engineering of microbial strains. BMC Syst Biol, 2010, 4: 53.

    [24] Tepper N, Shlomi T. Predicting metabolic engineering knockout strategies for chemical production: accounting for competing pathways. Bioinformatics, 2010, 26(4): 536–543.

    [25] Chowdhury A, Zomorrodi AR, Maranas CD. k-OptForce: Integrating kinetics with flux balance analysis for strain design. PLoS Comput Biol, 2014, 10(2): e1003487.

    [26] Khodayari A, Chowdhury A, Maranas CD. Succinate overproduction: A case study of computational strain design using a comprehensivekinetic model. Front Bioeng Biotechnol, 2014, 2: 76.

    [27] Ranganathan S, Tee TW, Chowdhury A, et al. An integrated computational and experimental study for overproducing fatty acids in. Metab Eng, 2012, 14(6): 687–704.

    [28] Liu WS, Zhang B, Jiang RR. Improving acetyl-CoA biosynthesis invia the overexpression of pantothenate kinase and PDH bypass. Biotechnol Biofuels, 2017, 10: 41.

    [29] Eckermann C, Schr?der G, Eckermann S, et al. Stilbenecarboxylate biosynthesis: a new function in the family of chalcone synthase-related proteins. Phytochemistry, 2003, 62(3): 271–286.

    [30] Sun W, Meng XY, Liang LJ, et al. Molecular and biochemical analysis of chalcone synthase fromin flavonoid biosynthetic pathway. PLoS ONE, 2015, 10(3): e0119054.

    [31] Chang JS, Winston F. Spt10 and spt21 are required for transcriptional silencing in. Eukaryot Cell, 2011, 10(1): 118–129.

    [32] Fang C, Wang QH, Selvaraj JN, et al. High copy and stable expression of the xylanase XynHB inby rDNA-mediated integration. Sci Rep, 2017, 7: 8747.

    Quantitative effect of the expression level of key genes in naringenin synthesis on the accumulation level of target products

    Tingting Jiao1,2, Jingwen Zhou1,2, and Sha Xu1,2

    1 National Engineering Laboratory for Cereal Fermentation Technology, Jiangnan University, Wuxi 214122, Jiangsu, China 2 School of Biotechnology, Jiangnan University, Wuxi 214122, Jiangsu, China

    Naringenin is a natural flavonoid compound with anti-inflammatory, anti-oxidation, anti-viral, anti-atherosclerosis and other pharmacological activities. It is also an important precursor of other flavonoid synthesis and with great value of application. At present, the production of flavonoids such as naringenin by microbial methods has a low yield due to imbalance of metabolic pathways, which greatly limits its industrial application. In this study, a naringenin-producing strain ofY-01 was used in the research object. The expression levels of 4-coumaric acid: CoA ligase (4CL), chalcone synthase (CHS) and chalcone isomerase (CHI) were controlled by promoter and copy numbers to investigate the quantitative effect of key enzyme expression level on the accumulation level of target products. The results showed that the correlation between naringenin production and 4CL or CHI expression was not significant while there was a positive correlation with the expression level of CHS. Strain Y-04 with high yield of naringenin was obtained by regulating the expression level ofgene, and the yield was increased by 4.1-folds compared with the original strain Y-01. This study indicated that CHS is a key regulatory target of naringenin synthesis. Rational regulation of CHSexpression can significantly promote the accumulation of naringenin. The related results provide an important theoretical reference for the use of metabolic engineering to strengthen microbial synthesis of important flavonoids such as naringenin.

    naringenin, expression level, key gene, promoter, copy number

    December 19, 2018;

    February 28, 2019

    National Natural Science Foundation of China (No. 31670095).

    Sha Xu. Tel/Fax: +86-510-85914371; E-mail: xus1984@jiangnan.edu.cn

    國(guó)家自然科學(xué)基金 (No. 31670095) 資助。

    2019-03-19

    http://kns.cnki.net/kcms/detail/11.1998.q.20190318.1121.002.html

    焦亭亭, 周景文, 徐沙. 柚皮素合成關(guān)鍵基因表達(dá)水平對(duì)目標(biāo)產(chǎn)物積累水平的量化影響. 生物工程學(xué)報(bào), 2019, 35(7): 1256–1265.

    Jiao TT, Zhou JW, Xu S. Quantitative effect of the expression level of key genes in naringenin synthesis on the accumulation level of target products. Chin J Biotech, 2019, 35(7): 1256–1265.

    (本文責(zé)編 陳宏宇)

    猜你喜歡
    柚皮素傳代拷貝數(shù)
    柚皮素對(duì)胰脂肪酶抑制作用的研究
    嗜酸乳桿菌和雙歧桿菌遺傳性對(duì)其益生特性及貨架期菌數(shù)的影響
    線粒體DNA拷貝數(shù)變異機(jī)制及疾病預(yù)測(cè)價(jià)值分析
    塞尼卡病毒A在BHK-21細(xì)胞中培養(yǎng)條件探索
    柚皮素脂質(zhì)體對(duì)非酒精性脂肪性肝病大鼠脂質(zhì)代謝的影響及機(jī)制研究*
    胎兒染色體組拷貝數(shù)變異與產(chǎn)前超聲異常的相關(guān)性分析
    柚皮素磷脂復(fù)合物的制備和表征
    中成藥(2018年7期)2018-08-04 06:04:18
    柚皮素-銅配合物的合成、表征及其抗氧化活性研究
    微生物實(shí)驗(yàn)室菌種復(fù)蘇傳代保存流程分析
    DNA序列拷貝數(shù)變化決定黃瓜性別
    亚洲欧美精品专区久久| 九九久久精品国产亚洲av麻豆| 国产探花极品一区二区| 亚洲久久久国产精品| tube8黄色片| 女的被弄到高潮叫床怎么办| 国产爽快片一区二区三区| 人妻少妇偷人精品九色| 最后的刺客免费高清国语| 噜噜噜噜噜久久久久久91| 亚洲自偷自拍三级| 国产成人免费观看mmmm| 成人国产麻豆网| 看非洲黑人一级黄片| 国产精品秋霞免费鲁丝片| 国产欧美日韩综合在线一区二区 | 波野结衣二区三区在线| 成年人免费黄色播放视频 | 欧美亚洲 丝袜 人妻 在线| 成人毛片60女人毛片免费| 亚洲经典国产精华液单| 日韩熟女老妇一区二区性免费视频| 九草在线视频观看| 大香蕉久久网| 欧美老熟妇乱子伦牲交| 国产乱人偷精品视频| 在线观看人妻少妇| 亚洲美女搞黄在线观看| 成年av动漫网址| 久久精品国产鲁丝片午夜精品| 国产精品久久久久久久久免| 免费黄色在线免费观看| 国产精品久久久久久av不卡| 久久 成人 亚洲| 日本爱情动作片www.在线观看| 国产视频首页在线观看| 一级毛片 在线播放| 人人妻人人澡人人看| 91aial.com中文字幕在线观看| 最新的欧美精品一区二区| 欧美bdsm另类| 久久久久视频综合| 精品一区二区免费观看| 久久精品久久久久久噜噜老黄| 成年美女黄网站色视频大全免费 | 国产精品一区二区在线观看99| 午夜福利在线观看免费完整高清在| 日韩av不卡免费在线播放| 91久久精品电影网| 国产永久视频网站| 一级,二级,三级黄色视频| 国产精品成人在线| 看非洲黑人一级黄片| 国产精品人妻久久久影院| 日本猛色少妇xxxxx猛交久久| 精品久久国产蜜桃| 只有这里有精品99| 69精品国产乱码久久久| 亚洲欧洲国产日韩| 精品久久久久久电影网| 日韩三级伦理在线观看| 久久久国产一区二区| 日韩中文字幕视频在线看片| 久久人人爽人人爽人人片va| 精品一区在线观看国产| 青青草视频在线视频观看| 亚洲欧洲精品一区二区精品久久久 | 亚洲美女视频黄频| 黑人高潮一二区| 夫妻午夜视频| 亚洲国产精品一区三区| 中文字幕av电影在线播放| 激情五月婷婷亚洲| 亚洲精品国产成人久久av| 99九九在线精品视频 | 亚洲成人av在线免费| 久久久久久久久久人人人人人人| 国产av精品麻豆| 久久99一区二区三区| 永久免费av网站大全| 边亲边吃奶的免费视频| 久久99精品国语久久久| 日本vs欧美在线观看视频 | 嘟嘟电影网在线观看| 免费人成在线观看视频色| 国产伦精品一区二区三区四那| 免费看不卡的av| 色5月婷婷丁香| 免费看光身美女| 欧美日韩一区二区视频在线观看视频在线| 国产成人精品久久久久久| 男女边摸边吃奶| 最新的欧美精品一区二区| 久久精品国产a三级三级三级| 亚洲综合色惰| 国产精品一区二区性色av| 秋霞伦理黄片| 亚洲人成网站在线播| 国产成人免费无遮挡视频| 99热网站在线观看| 欧美精品人与动牲交sv欧美| 国产伦理片在线播放av一区| 中文字幕免费在线视频6| 王馨瑶露胸无遮挡在线观看| 亚洲av.av天堂| 国产免费一级a男人的天堂| 国产精品国产av在线观看| 国产69精品久久久久777片| 久久韩国三级中文字幕| 桃花免费在线播放| 久久热精品热| 日日摸夜夜添夜夜添av毛片| 三上悠亚av全集在线观看 | 欧美激情极品国产一区二区三区 | 亚洲精品第二区| 欧美日韩亚洲高清精品| 中文字幕制服av| 又爽又黄a免费视频| 黄色欧美视频在线观看| 男女无遮挡免费网站观看| 成人综合一区亚洲| 在线观看免费视频网站a站| 有码 亚洲区| 欧美高清成人免费视频www| 国产日韩欧美视频二区| 亚洲,一卡二卡三卡| 亚洲精品国产色婷婷电影| 午夜免费鲁丝| 人妻人人澡人人爽人人| 国产视频内射| 又爽又黄a免费视频| 中国美白少妇内射xxxbb| 久久久久久久久久人人人人人人| 日产精品乱码卡一卡2卡三| 2018国产大陆天天弄谢| av在线播放精品| 曰老女人黄片| 亚洲美女搞黄在线观看| 赤兔流量卡办理| 黑丝袜美女国产一区| av播播在线观看一区| 亚洲国产欧美在线一区| 黄色怎么调成土黄色| 在现免费观看毛片| 嘟嘟电影网在线观看| 青春草视频在线免费观看| 女性生殖器流出的白浆| 亚洲av中文av极速乱| 精品一区二区免费观看| 国产综合精华液| 水蜜桃什么品种好| 精品亚洲成a人片在线观看| 久久久久久久久久成人| 中文天堂在线官网| 69精品国产乱码久久久| 九色成人免费人妻av| 在线观看美女被高潮喷水网站| 亚洲av免费高清在线观看| 天堂俺去俺来也www色官网| 精品人妻一区二区三区麻豆| 久久精品国产亚洲av涩爱| 久久99一区二区三区| 又大又黄又爽视频免费| 男女边吃奶边做爰视频| av卡一久久| 在线观看www视频免费| 国产日韩欧美在线精品| 黄色日韩在线| 欧美日韩视频精品一区| a级毛片免费高清观看在线播放| 亚洲成人手机| 嘟嘟电影网在线观看| 日韩在线高清观看一区二区三区| 在线天堂最新版资源| 日本午夜av视频| 国产高清国产精品国产三级| 国精品久久久久久国模美| 女性被躁到高潮视频| 国产爽快片一区二区三区| 熟女av电影| 男女边摸边吃奶| 国产精品人妻久久久影院| 精品卡一卡二卡四卡免费| 国产精品99久久99久久久不卡 | 午夜福利视频精品| 99热6这里只有精品| 欧美一级a爱片免费观看看| 黄色配什么色好看| 久久99精品国语久久久| 高清av免费在线| 国产免费一区二区三区四区乱码| 伊人久久国产一区二区| 国产精品人妻久久久久久| 天堂中文最新版在线下载| 精品久久久久久电影网| 色视频在线一区二区三区| 国产成人aa在线观看| 伊人久久国产一区二区| 少妇被粗大的猛进出69影院 | 人人澡人人妻人| 久久久久精品性色| 亚洲av综合色区一区| 国产日韩欧美在线精品| 女性生殖器流出的白浆| 在线观看免费日韩欧美大片 | 欧美精品国产亚洲| 亚洲高清免费不卡视频| 妹子高潮喷水视频| av视频免费观看在线观看| 国产熟女欧美一区二区| 成人黄色视频免费在线看| 国产精品三级大全| 精品久久久精品久久久| 亚洲美女搞黄在线观看| 日韩成人伦理影院| 国产成人精品无人区| 爱豆传媒免费全集在线观看| 六月丁香七月| 麻豆精品久久久久久蜜桃| 久久毛片免费看一区二区三区| 高清午夜精品一区二区三区| 国产精品成人在线| 欧美另类一区| 男男h啪啪无遮挡| 久久久久久人妻| 女人久久www免费人成看片| 欧美日韩在线观看h| 国产视频首页在线观看| 色婷婷久久久亚洲欧美| 久久国产精品男人的天堂亚洲 | 色婷婷久久久亚洲欧美| 一个人免费看片子| av.在线天堂| 亚洲av成人精品一二三区| 国精品久久久久久国模美| 亚洲精品中文字幕在线视频 | 国产亚洲5aaaaa淫片| 国产亚洲最大av| 最近的中文字幕免费完整| 成年女人在线观看亚洲视频| 免费观看av网站的网址| 亚洲伊人久久精品综合| 久久99热这里只频精品6学生| 国产乱来视频区| 人妻 亚洲 视频| 国产成人精品久久久久久| 18禁裸乳无遮挡动漫免费视频| 国模一区二区三区四区视频| 亚洲人成网站在线播| 精品酒店卫生间| 在线观看免费日韩欧美大片 | 青春草亚洲视频在线观看| 女性被躁到高潮视频| 十八禁高潮呻吟视频 | 嫩草影院新地址| 热re99久久精品国产66热6| 欧美高清成人免费视频www| 亚洲av欧美aⅴ国产| 黑人高潮一二区| 我要看黄色一级片免费的| 欧美精品亚洲一区二区| 国产乱人偷精品视频| 欧美日韩综合久久久久久| 又大又黄又爽视频免费| 日本av手机在线免费观看| 亚洲成人一二三区av| 亚洲精华国产精华液的使用体验| 天天躁夜夜躁狠狠久久av| 九色成人免费人妻av| 曰老女人黄片| 国产伦精品一区二区三区四那| 永久网站在线| 99久久中文字幕三级久久日本| 大片免费播放器 马上看| 亚洲国产毛片av蜜桃av| 成人毛片a级毛片在线播放| 国产精品免费大片| 青春草国产在线视频| 日本黄色日本黄色录像| 免费播放大片免费观看视频在线观看| 午夜福利视频精品| 十八禁高潮呻吟视频 | 亚洲欧美清纯卡通| 欧美精品高潮呻吟av久久| 亚洲av日韩在线播放| 精品久久国产蜜桃| 国产69精品久久久久777片| 黄色配什么色好看| 国产伦精品一区二区三区视频9| 国产成人a∨麻豆精品| 中文在线观看免费www的网站| 三级国产精品欧美在线观看| 亚洲av福利一区| 熟女人妻精品中文字幕| 日韩视频在线欧美| 交换朋友夫妻互换小说| av一本久久久久| 韩国高清视频一区二区三区| www.av在线官网国产| www.色视频.com| 国产亚洲精品久久久com| 一级毛片aaaaaa免费看小| 亚洲激情五月婷婷啪啪| 久久综合国产亚洲精品| 国产精品蜜桃在线观看| 女人久久www免费人成看片| 国产一区二区在线观看av| 春色校园在线视频观看| 欧美97在线视频| 亚洲一区二区三区欧美精品| 国产日韩一区二区三区精品不卡 | 久久久午夜欧美精品| 亚洲情色 制服丝袜| 亚洲国产欧美日韩在线播放 | av国产精品久久久久影院| 国产亚洲av片在线观看秒播厂| 国产免费福利视频在线观看| 久久99蜜桃精品久久| 十分钟在线观看高清视频www | 少妇裸体淫交视频免费看高清| av卡一久久| 欧美xxxx性猛交bbbb| 九九爱精品视频在线观看| 久久这里有精品视频免费| 看免费成人av毛片| 中文天堂在线官网| 亚洲av成人精品一二三区| 精品一区在线观看国产| 99九九在线精品视频 | 亚洲成人av在线免费| 中文天堂在线官网| 3wmmmm亚洲av在线观看| 国产亚洲91精品色在线| 亚洲国产精品一区三区| 夜夜骑夜夜射夜夜干| 成人二区视频| 人妻 亚洲 视频| 欧美最新免费一区二区三区| 街头女战士在线观看网站| 国产成人免费观看mmmm| 国产成人a∨麻豆精品| a 毛片基地| 中文字幕人妻熟人妻熟丝袜美| 亚洲电影在线观看av| 国产精品欧美亚洲77777| 在线亚洲精品国产二区图片欧美 | 亚洲欧美一区二区三区黑人 | 在现免费观看毛片| 18+在线观看网站| 亚洲熟女精品中文字幕| 春色校园在线视频观看| 国产av一区二区精品久久| 观看av在线不卡| 性色avwww在线观看| 国产精品国产av在线观看| 午夜av观看不卡| 亚洲欧美清纯卡通| 国产精品国产av在线观看| 国产一区二区在线观看日韩| 美女cb高潮喷水在线观看| 观看免费一级毛片| a级毛片在线看网站| 欧美日韩视频高清一区二区三区二| 高清在线视频一区二区三区| 亚洲国产精品一区二区三区在线| 18禁在线播放成人免费| 乱人伦中国视频| 久久热精品热| 免费观看无遮挡的男女| 久久久久久久亚洲中文字幕| 国产精品国产三级专区第一集| 男男h啪啪无遮挡| 国产淫语在线视频| 欧美精品一区二区免费开放| 一区二区三区免费毛片| 久久狼人影院| 如日韩欧美国产精品一区二区三区 | 国产伦在线观看视频一区| 日韩制服骚丝袜av| 熟女人妻精品中文字幕| 国产一区二区三区综合在线观看 | 成人综合一区亚洲| 日本色播在线视频| 在现免费观看毛片| 日日摸夜夜添夜夜添av毛片| 搡女人真爽免费视频火全软件| 热99国产精品久久久久久7| 亚洲国产精品999| 十分钟在线观看高清视频www | 国产伦理片在线播放av一区| 亚洲综合色惰| 一级黄片播放器| 夫妻性生交免费视频一级片| 欧美少妇被猛烈插入视频| 午夜精品国产一区二区电影| 永久网站在线| 亚洲,一卡二卡三卡| 亚洲丝袜综合中文字幕| 看非洲黑人一级黄片| 亚洲国产精品成人久久小说| 伊人久久精品亚洲午夜| 国产国拍精品亚洲av在线观看| 中文字幕免费在线视频6| 91久久精品电影网| 在线观看免费高清a一片| 人人妻人人看人人澡| 国产精品国产三级国产av玫瑰| 国产精品久久久久久精品古装| 日本免费在线观看一区| 国产综合精华液| 欧美一级a爱片免费观看看| 日韩中字成人| 欧美精品亚洲一区二区| 国产黄片美女视频| 91久久精品国产一区二区三区| 亚洲av在线观看美女高潮| 免费高清在线观看视频在线观看| 最近手机中文字幕大全| 高清视频免费观看一区二区| 在线观看三级黄色| 欧美3d第一页| 国产免费一级a男人的天堂| 国产成人免费无遮挡视频| 人妻 亚洲 视频| 男人舔奶头视频| a级毛片免费高清观看在线播放| 免费在线观看成人毛片| a级一级毛片免费在线观看| 男女啪啪激烈高潮av片| 精品久久国产蜜桃| 韩国av在线不卡| 欧美老熟妇乱子伦牲交| 国产精品成人在线| av国产精品久久久久影院| 亚洲国产成人一精品久久久| 少妇高潮的动态图| 亚洲欧美成人综合另类久久久| 日日爽夜夜爽网站| 99久久人妻综合| 女性生殖器流出的白浆| 97超视频在线观看视频| 久久精品夜色国产| 亚洲精品一二三| 久久久精品94久久精品| 精品酒店卫生间| 亚洲天堂av无毛| 精品人妻熟女av久视频| 黄色日韩在线| 亚洲,一卡二卡三卡| 中文资源天堂在线| 亚洲国产欧美在线一区| 内地一区二区视频在线| 日本av免费视频播放| 秋霞在线观看毛片| 美女脱内裤让男人舔精品视频| 热re99久久国产66热| 亚洲在久久综合| 激情五月婷婷亚洲| 日韩视频在线欧美| 最近最新中文字幕免费大全7| 亚洲av在线观看美女高潮| 国产欧美日韩综合在线一区二区 | 日韩成人伦理影院| 亚洲精品日韩在线中文字幕| 一本一本综合久久| av.在线天堂| 大片免费播放器 马上看| 一个人看视频在线观看www免费| 午夜激情久久久久久久| 中文字幕制服av| 亚洲美女黄色视频免费看| 久久精品夜色国产| 制服丝袜香蕉在线| 嫩草影院入口| 亚洲在久久综合| 中国国产av一级| 赤兔流量卡办理| 久久人人爽人人爽人人片va| 丝袜喷水一区| 久久精品久久久久久噜噜老黄| 免费看av在线观看网站| 91午夜精品亚洲一区二区三区| 欧美 亚洲 国产 日韩一| 国产av精品麻豆| 日韩欧美精品免费久久| 看十八女毛片水多多多| 9色porny在线观看| 18+在线观看网站| 亚洲欧洲精品一区二区精品久久久 | 日本黄色日本黄色录像| av一本久久久久| 99久久精品一区二区三区| 亚洲av男天堂| 特大巨黑吊av在线直播| 国产精品国产三级专区第一集| a 毛片基地| 永久免费av网站大全| 久久国内精品自在自线图片| 国产男女超爽视频在线观看| 国产av国产精品国产| 欧美+日韩+精品| 偷拍熟女少妇极品色| 18+在线观看网站| 一级毛片电影观看| 国产淫语在线视频| 99热网站在线观看| 高清在线视频一区二区三区| av福利片在线| 韩国高清视频一区二区三区| 最新的欧美精品一区二区| 老司机影院成人| 国产成人精品久久久久久| 春色校园在线视频观看| 欧美三级亚洲精品| 一区二区三区免费毛片| 色94色欧美一区二区| 91久久精品国产一区二区成人| 日韩,欧美,国产一区二区三区| 久久综合国产亚洲精品| 麻豆成人av视频| 综合色丁香网| 亚洲久久久国产精品| 日韩成人av中文字幕在线观看| 伦精品一区二区三区| 青青草视频在线视频观看| 精品少妇黑人巨大在线播放| 综合色丁香网| 不卡视频在线观看欧美| 日日爽夜夜爽网站| 国产成人精品福利久久| 卡戴珊不雅视频在线播放| 亚洲在久久综合| 一二三四中文在线观看免费高清| 精品少妇黑人巨大在线播放| 天美传媒精品一区二区| 久久精品国产亚洲av涩爱| 国产综合精华液| 久久精品国产亚洲av涩爱| 少妇裸体淫交视频免费看高清| 伦精品一区二区三区| 久久人人爽av亚洲精品天堂| 亚洲国产精品999| 综合色丁香网| 欧美+日韩+精品| 黄色怎么调成土黄色| 久久久久国产网址| 校园人妻丝袜中文字幕| 亚洲国产精品999| 一级二级三级毛片免费看| 免费看光身美女| 观看美女的网站| 岛国毛片在线播放| 99久久综合免费| 高清午夜精品一区二区三区| 亚洲综合精品二区| 久久久久久久久久人人人人人人| 欧美日韩国产mv在线观看视频| 你懂的网址亚洲精品在线观看| 国产成人免费无遮挡视频| 欧美最新免费一区二区三区| videossex国产| 人体艺术视频欧美日本| 男人和女人高潮做爰伦理| 久久99一区二区三区| a 毛片基地| 日韩一本色道免费dvd| 秋霞在线观看毛片| 街头女战士在线观看网站| 久久精品久久久久久久性| 五月伊人婷婷丁香| 三上悠亚av全集在线观看 | 18禁动态无遮挡网站| 国产欧美日韩综合在线一区二区 | 国产白丝娇喘喷水9色精品| 能在线免费看毛片的网站| 成人毛片60女人毛片免费| 乱系列少妇在线播放| 只有这里有精品99| 亚洲精品中文字幕在线视频 | 午夜福利视频精品| 最近中文字幕2019免费版| 一本一本综合久久| 日韩av不卡免费在线播放| 国产亚洲一区二区精品| 精品久久久久久久久亚洲| 内地一区二区视频在线| 十分钟在线观看高清视频www | 国产免费福利视频在线观看| 国产精品99久久99久久久不卡 | 人体艺术视频欧美日本| videossex国产| 日本黄大片高清| 最近2019中文字幕mv第一页| 亚洲精品,欧美精品| 成年av动漫网址| 一级毛片 在线播放| 精品一区二区三卡| 精品一区在线观看国产| 日本黄色日本黄色录像| 在线观看美女被高潮喷水网站| 亚洲中文av在线| 亚洲va在线va天堂va国产| 男的添女的下面高潮视频| 国产精品国产三级国产专区5o| 亚洲精品日韩av片在线观看| 涩涩av久久男人的天堂| 亚洲精品色激情综合| 午夜视频国产福利| 国产精品女同一区二区软件| 色94色欧美一区二区| 热99国产精品久久久久久7| 天天操日日干夜夜撸| 丰满饥渴人妻一区二区三| 精品久久久久久久久亚洲| 日日啪夜夜爽|