• <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ù)變化決定黃瓜性別
    欧美成人a在线观看| АⅤ资源中文在线天堂| 又粗又硬又长又爽又黄的视频| av视频在线观看入口| 亚洲高清免费不卡视频| 国产老妇伦熟女老妇高清| 欧美成人免费av一区二区三区| 久久久久久久国产电影| 国产精品电影一区二区三区| 天堂av国产一区二区熟女人妻| 欧美+日韩+精品| 亚洲第一区二区三区不卡| 精品人妻偷拍中文字幕| 亚洲av成人精品一二三区| www.色视频.com| 日韩强制内射视频| 午夜福利在线观看吧| 国产亚洲av嫩草精品影院| 成人二区视频| 国产精品蜜桃在线观看| 男女视频在线观看网站免费| 国产黄片视频在线免费观看| 99在线视频只有这里精品首页| 欧美丝袜亚洲另类| 听说在线观看完整版免费高清| 欧美激情久久久久久爽电影| 精品少妇黑人巨大在线播放 | 欧美成人午夜免费资源| 精品久久久久久成人av| 麻豆av噜噜一区二区三区| 国产麻豆成人av免费视频| 两个人视频免费观看高清| 亚洲国产欧美人成| 国产爱豆传媒在线观看| 免费看日本二区| 亚洲av成人精品一二三区| 国产成人freesex在线| 久久午夜福利片| 超碰av人人做人人爽久久| 亚洲欧美日韩卡通动漫| 久久精品夜色国产| 欧美极品一区二区三区四区| 最近的中文字幕免费完整| 国产片特级美女逼逼视频| av女优亚洲男人天堂| 99在线视频只有这里精品首页| 在线免费观看不下载黄p国产| 午夜爱爱视频在线播放| 97人妻精品一区二区三区麻豆| 看黄色毛片网站| 国产精品电影一区二区三区| 九九热线精品视视频播放| 国产精品三级大全| av在线蜜桃| 久久亚洲国产成人精品v| 大又大粗又爽又黄少妇毛片口| 午夜精品国产一区二区电影 | 国产精品精品国产色婷婷| 久久99精品国语久久久| 亚洲成人精品中文字幕电影| 又粗又硬又长又爽又黄的视频| 国产精品国产三级国产av玫瑰| 亚洲av成人精品一二三区| 午夜a级毛片| 噜噜噜噜噜久久久久久91| 久久精品久久精品一区二区三区| 91在线精品国自产拍蜜月| 国产精品国产三级专区第一集| 久久久久久久久大av| 国产成人精品婷婷| 亚洲av男天堂| 18禁在线播放成人免费| 亚洲高清免费不卡视频| 日本欧美国产在线视频| 成年av动漫网址| 色5月婷婷丁香| 国产欧美日韩精品一区二区| 精品一区二区免费观看| 亚洲国产精品专区欧美| 久久久久久大精品| 禁无遮挡网站| 亚洲欧美精品综合久久99| 中文字幕人妻熟人妻熟丝袜美| 日韩 亚洲 欧美在线| 女人久久www免费人成看片 | 午夜福利在线观看免费完整高清在| 国产熟女欧美一区二区| 九九在线视频观看精品| 亚洲精品日韩av片在线观看| 亚洲av免费在线观看| 嫩草影院精品99| 毛片女人毛片| 国产精品伦人一区二区| 大又大粗又爽又黄少妇毛片口| 九九在线视频观看精品| 色噜噜av男人的天堂激情| 在线观看美女被高潮喷水网站| 97热精品久久久久久| 男女边吃奶边做爰视频| 国产精品三级大全| 91久久精品国产一区二区成人| 久久人人爽人人片av| 啦啦啦韩国在线观看视频| 国产熟女欧美一区二区| 丰满人妻一区二区三区视频av| 色播亚洲综合网| 建设人人有责人人尽责人人享有的 | 老司机影院毛片| 91av网一区二区| 亚洲国产欧美人成| 亚洲aⅴ乱码一区二区在线播放| 国模一区二区三区四区视频| 麻豆乱淫一区二区| 国产av一区在线观看免费| 国产伦理片在线播放av一区| 一个人看视频在线观看www免费| 永久免费av网站大全| 久久国产乱子免费精品| 床上黄色一级片| 欧美成人a在线观看| 国产黄色视频一区二区在线观看 | 少妇丰满av| 男女国产视频网站| 国产精品久久视频播放| 干丝袜人妻中文字幕| 亚洲丝袜综合中文字幕| 国产午夜福利久久久久久| 久久久久精品久久久久真实原创| 久久精品综合一区二区三区| 亚洲精品456在线播放app| 国产精品乱码一区二三区的特点| 精品久久久久久电影网 | 亚洲激情五月婷婷啪啪| 日韩av在线大香蕉| 黄色欧美视频在线观看| 哪个播放器可以免费观看大片| 女人十人毛片免费观看3o分钟| 久久久久久久久久久丰满| 亚洲电影在线观看av| 有码 亚洲区| 色视频www国产| 伦理电影大哥的女人| 美女黄网站色视频| 久久久国产成人精品二区| 成年版毛片免费区| 欧美精品国产亚洲| 女的被弄到高潮叫床怎么办| 国产av码专区亚洲av| 亚洲内射少妇av| 超碰97精品在线观看| 久久久久久久久久久免费av| 爱豆传媒免费全集在线观看| 最近的中文字幕免费完整| 一个人看视频在线观看www免费| 99久久精品热视频| 一级黄片播放器| 伊人久久精品亚洲午夜| 久久精品久久久久久噜噜老黄 | 日日啪夜夜撸| 丰满少妇做爰视频| 91狼人影院| 国产伦精品一区二区三区视频9| 毛片女人毛片| 久久久久免费精品人妻一区二区| 热99在线观看视频| 国产美女午夜福利| 一级黄色大片毛片| 深爱激情五月婷婷| 亚洲人成网站高清观看| 婷婷色综合大香蕉| 水蜜桃什么品种好| 精品不卡国产一区二区三区| 久久久精品94久久精品| 国产在线一区二区三区精 | 欧美成人一区二区免费高清观看| 国产欧美日韩精品一区二区| 3wmmmm亚洲av在线观看| 在线免费十八禁| 欧美人与善性xxx| www.色视频.com| 精品人妻熟女av久视频| 亚洲av二区三区四区| 99在线人妻在线中文字幕| 欧美丝袜亚洲另类| 小说图片视频综合网站| 国产黄色小视频在线观看| 欧美zozozo另类| 国产精品久久电影中文字幕| 一区二区三区高清视频在线| 成人午夜高清在线视频| 天堂av国产一区二区熟女人妻| 日韩中字成人| 在线免费观看不下载黄p国产| 成年版毛片免费区| 亚洲国产欧洲综合997久久,| 免费搜索国产男女视频| 亚洲av成人精品一二三区| 亚洲第一区二区三区不卡| 国产不卡一卡二| 舔av片在线| 成人av在线播放网站| 日韩视频在线欧美| 寂寞人妻少妇视频99o| 色综合亚洲欧美另类图片| 欧美激情在线99| 国产91av在线免费观看| 99久久九九国产精品国产免费| 在线观看一区二区三区| 免费观看人在逋| 十八禁国产超污无遮挡网站| 精品久久久久久成人av| 大又大粗又爽又黄少妇毛片口| 亚洲中文字幕日韩| 日本猛色少妇xxxxx猛交久久| 亚洲欧美中文字幕日韩二区| 别揉我奶头 嗯啊视频| 人人妻人人澡人人爽人人夜夜 | 亚洲最大成人中文| 男女啪啪激烈高潮av片| 亚洲av.av天堂| 一级毛片我不卡| 99热这里只有精品一区| 人妻制服诱惑在线中文字幕| 欧美又色又爽又黄视频| 99久久中文字幕三级久久日本| 综合色av麻豆| 天堂网av新在线| 国产精品久久久久久久电影| 七月丁香在线播放| 美女cb高潮喷水在线观看| 国产亚洲av片在线观看秒播厂 | 国产大屁股一区二区在线视频| 久久精品国产亚洲av涩爱| 欧美又色又爽又黄视频| 听说在线观看完整版免费高清| 乱系列少妇在线播放| 成人午夜精彩视频在线观看| 国产女主播在线喷水免费视频网站 | 欧美一区二区精品小视频在线| 有码 亚洲区| 18禁在线播放成人免费| 我的女老师完整版在线观看| 一二三四中文在线观看免费高清| 91精品国产九色| 国产日韩欧美在线精品| 日本-黄色视频高清免费观看| 九九久久精品国产亚洲av麻豆| 18禁在线无遮挡免费观看视频| 欧美高清性xxxxhd video| 国产av码专区亚洲av| 九色成人免费人妻av| 国产淫片久久久久久久久| av专区在线播放| 波多野结衣高清无吗| 美女xxoo啪啪120秒动态图| 男人狂女人下面高潮的视频| 免费看光身美女| 亚洲怡红院男人天堂| 超碰av人人做人人爽久久| 国产免费一级a男人的天堂| 亚洲精品,欧美精品| 99久久精品热视频| 久久99热这里只频精品6学生 | 久久韩国三级中文字幕| 男人和女人高潮做爰伦理| 一级爰片在线观看| 亚洲av日韩在线播放| 少妇的逼好多水| 午夜激情福利司机影院| 看黄色毛片网站| 亚洲欧美日韩无卡精品| av卡一久久| 日本一本二区三区精品| 精品人妻偷拍中文字幕| 深爱激情五月婷婷| 国产精品美女特级片免费视频播放器| 91精品伊人久久大香线蕉| 久久久久久久久久久丰满| 亚洲av免费高清在线观看| 国产午夜精品论理片| 纵有疾风起免费观看全集完整版 | 欧美日韩国产亚洲二区| 少妇被粗大猛烈的视频| 国产片特级美女逼逼视频| 亚洲国产日韩欧美精品在线观看| 永久网站在线| 成人午夜精彩视频在线观看| 国产精品一区二区性色av| 久久久精品欧美日韩精品| 1024手机看黄色片| 国产亚洲最大av| 国产精品人妻久久久影院| 成年av动漫网址| 韩国av在线不卡| 国产高清视频在线观看网站| 六月丁香七月| 成年av动漫网址| av福利片在线观看| 国产国拍精品亚洲av在线观看| 亚洲欧美中文字幕日韩二区| 国产精品野战在线观看| av天堂中文字幕网| 国产一区二区在线av高清观看| 免费看av在线观看网站| 成年女人看的毛片在线观看| 青春草亚洲视频在线观看| 久久久久久久亚洲中文字幕| 国内少妇人妻偷人精品xxx网站| 亚洲国产精品成人综合色| 欧美高清性xxxxhd video| 国产熟女欧美一区二区| 如何舔出高潮| 99久久人妻综合| 欧美日韩一区二区视频在线观看视频在线 | 国产又黄又爽又无遮挡在线| 在线a可以看的网站| 我要看日韩黄色一级片| av线在线观看网站| 日韩欧美精品v在线| 国内精品宾馆在线| 不卡视频在线观看欧美| 久久精品综合一区二区三区| 欧美成人精品欧美一级黄| 赤兔流量卡办理| 午夜福利高清视频| 亚洲国产欧洲综合997久久,| 免费av不卡在线播放| 日韩大片免费观看网站 | 午夜亚洲福利在线播放| 午夜激情福利司机影院| 一个人看视频在线观看www免费| 2021天堂中文幕一二区在线观| 久久精品影院6| 国产精品一区www在线观看| 国产精品不卡视频一区二区| 亚洲精品乱码久久久久久按摩| 欧美zozozo另类| 成人综合一区亚洲| 日本午夜av视频| 99久久精品国产国产毛片| 看片在线看免费视频| 男女下面进入的视频免费午夜| av女优亚洲男人天堂| 不卡视频在线观看欧美| 国产精品一区二区三区四区免费观看| 国产精品1区2区在线观看.| 亚洲熟妇中文字幕五十中出| 久久精品国产亚洲av天美| 国产精品国产三级专区第一集| 99九九线精品视频在线观看视频| 国产综合懂色| 哪个播放器可以免费观看大片| 国产亚洲av嫩草精品影院| 欧美不卡视频在线免费观看| 色哟哟·www| 亚洲欧美一区二区三区国产| 中文亚洲av片在线观看爽| a级毛片免费高清观看在线播放| 青青草视频在线视频观看| 国产探花在线观看一区二区| 亚洲av.av天堂| 九九爱精品视频在线观看| 精品人妻视频免费看| 老司机影院成人| 国产伦理片在线播放av一区| 国产成人福利小说| 老女人水多毛片| 一区二区三区高清视频在线| 免费大片18禁| 亚洲在线自拍视频| 村上凉子中文字幕在线| 99在线视频只有这里精品首页| 高清视频免费观看一区二区 | 日产精品乱码卡一卡2卡三| 国产成人免费观看mmmm| 桃色一区二区三区在线观看| 1000部很黄的大片| 国产在视频线精品| 亚洲性久久影院| 国产男人的电影天堂91| 国产黄片视频在线免费观看| 成人欧美大片| 日韩三级伦理在线观看| 婷婷色av中文字幕| 国产精品一区二区性色av| 亚洲aⅴ乱码一区二区在线播放| 五月伊人婷婷丁香| 麻豆国产97在线/欧美| 青春草国产在线视频| 日韩一区二区视频免费看| 天堂中文最新版在线下载 | a级毛片免费高清观看在线播放| 精品人妻偷拍中文字幕| 久久6这里有精品| 亚洲一区高清亚洲精品| 黄色一级大片看看| 一边亲一边摸免费视频| 老司机福利观看| 日韩人妻高清精品专区| 2021少妇久久久久久久久久久| 精品酒店卫生间| 一级二级三级毛片免费看| 一级爰片在线观看| 免费看a级黄色片| 久久精品国产自在天天线| 亚洲自拍偷在线| 久久久色成人| 欧美不卡视频在线免费观看| 久久精品综合一区二区三区| 又黄又爽又刺激的免费视频.| 日本色播在线视频| 亚洲欧美中文字幕日韩二区| 黄色配什么色好看| 亚洲av一区综合| 91精品一卡2卡3卡4卡| 男女边吃奶边做爰视频| 精品久久久噜噜| 久久久久久久久久黄片| 一卡2卡三卡四卡精品乱码亚洲| 亚洲精品乱码久久久久久按摩| 日韩高清综合在线| 国产探花极品一区二区| 菩萨蛮人人尽说江南好唐韦庄 | 免费看av在线观看网站| 又爽又黄a免费视频| 青春草国产在线视频| 在线天堂最新版资源| 成年女人永久免费观看视频| 免费人成在线观看视频色| 老司机影院成人| 美女脱内裤让男人舔精品视频| 噜噜噜噜噜久久久久久91| 日韩一区二区视频免费看| 亚州av有码| 欧美极品一区二区三区四区| 美女大奶头视频| 国产精品一区www在线观看| av天堂中文字幕网| 国产亚洲午夜精品一区二区久久 | www日本黄色视频网| 亚洲人与动物交配视频| 国产久久久一区二区三区| 午夜福利在线观看免费完整高清在| 一级毛片我不卡| 一个人看的www免费观看视频| 欧美激情在线99| 国产精品乱码一区二三区的特点| 在线观看美女被高潮喷水网站| 久久久久久久久久成人| ponron亚洲| 精品欧美国产一区二区三| kizo精华| 久久久久久久午夜电影| 欧美高清成人免费视频www| 嫩草影院精品99| 欧美高清成人免费视频www| 青春草亚洲视频在线观看| 欧美日韩国产亚洲二区| 一夜夜www| 国内精品一区二区在线观看| 日日干狠狠操夜夜爽| 日韩,欧美,国产一区二区三区 | 亚洲经典国产精华液单| 欧美xxxx黑人xx丫x性爽| 欧美日韩综合久久久久久| 性色avwww在线观看| 中文亚洲av片在线观看爽| 性色avwww在线观看| 欧美三级亚洲精品| 精品欧美国产一区二区三| 身体一侧抽搐| 毛片女人毛片| 1000部很黄的大片| 久久久久九九精品影院| 亚洲电影在线观看av| 99久久精品国产国产毛片| 午夜福利在线观看免费完整高清在| 在现免费观看毛片| av在线蜜桃| 成人av在线播放网站| 成人毛片a级毛片在线播放| 成人av在线播放网站| 精品人妻偷拍中文字幕| 在线免费观看不下载黄p国产| 久久精品人妻少妇| 99热全是精品| 又粗又爽又猛毛片免费看| 韩国av在线不卡| 2021天堂中文幕一二区在线观| 少妇熟女欧美另类| 亚洲精品久久久久久婷婷小说 | 国产午夜精品论理片| 久久精品综合一区二区三区| 亚洲丝袜综合中文字幕| 少妇的逼好多水| 大又大粗又爽又黄少妇毛片口| 欧美一区二区国产精品久久精品| 中文字幕av在线有码专区| 有码 亚洲区| 18+在线观看网站| 欧美性猛交╳xxx乱大交人| 久久久欧美国产精品| 久久久精品欧美日韩精品| 色5月婷婷丁香| 校园人妻丝袜中文字幕| 日韩国内少妇激情av| 丝袜美腿在线中文| 少妇高潮的动态图| 一区二区三区乱码不卡18| 在线观看66精品国产| 亚洲最大成人手机在线| 国产免费一级a男人的天堂| 日韩 亚洲 欧美在线| 2022亚洲国产成人精品| 国产视频内射| 又粗又爽又猛毛片免费看| 欧美xxxx性猛交bbbb| 国产精品综合久久久久久久免费| 狂野欧美白嫩少妇大欣赏| 日韩欧美精品v在线| 特大巨黑吊av在线直播| 日本wwww免费看| 真实男女啪啪啪动态图| 国产成年人精品一区二区| 精品国内亚洲2022精品成人| 高清av免费在线| 97超碰精品成人国产| 岛国在线免费视频观看| 亚洲成人av在线免费| 色播亚洲综合网| 色网站视频免费| 好男人在线观看高清免费视频| 久久精品人妻少妇| 成人av在线播放网站| 国产亚洲5aaaaa淫片| 欧美3d第一页| 国产精品av视频在线免费观看| 又爽又黄a免费视频| 亚洲av成人av| 亚洲三级黄色毛片| 成人漫画全彩无遮挡| 亚洲国产精品合色在线| 嫩草影院入口| 日日摸夜夜添夜夜添av毛片| 丰满人妻一区二区三区视频av| 成人美女网站在线观看视频| 少妇的逼好多水| 国产精品人妻久久久久久| 欧美性感艳星| 中文字幕人妻熟人妻熟丝袜美| 日日撸夜夜添| 国产成人a区在线观看| 国产男人的电影天堂91| 又黄又爽又刺激的免费视频.| 91av网一区二区| 亚洲精品一区蜜桃| 亚洲国产欧洲综合997久久,| 欧美变态另类bdsm刘玥| 国产一区二区亚洲精品在线观看| 亚洲丝袜综合中文字幕| 成年女人永久免费观看视频| 国产探花在线观看一区二区| 成年女人永久免费观看视频| 内射极品少妇av片p| 国产成人a∨麻豆精品| 国产极品精品免费视频能看的| 天美传媒精品一区二区| 国产精品一区二区在线观看99 | 爱豆传媒免费全集在线观看| 久久久久久久久中文| 免费观看性生交大片5| 欧美不卡视频在线免费观看| 最近手机中文字幕大全| 久久热精品热| 午夜激情福利司机影院| 村上凉子中文字幕在线| 九九爱精品视频在线观看| 免费av毛片视频| 久久久久久久亚洲中文字幕| 日日啪夜夜撸| 高清av免费在线| 国产精品电影一区二区三区| 国产亚洲最大av| 天天一区二区日本电影三级| 亚洲国产最新在线播放| 在线天堂最新版资源| 久久午夜福利片| 永久网站在线| 久久精品夜色国产| 免费看a级黄色片| 超碰97精品在线观看| 成人国产麻豆网| 国产精品麻豆人妻色哟哟久久 | 精品久久久久久电影网 | 久久精品影院6| 日日摸夜夜添夜夜添av毛片| 色播亚洲综合网| 久久精品91蜜桃| 搡女人真爽免费视频火全软件| 好男人在线观看高清免费视频| 日本免费一区二区三区高清不卡| 99热这里只有是精品在线观看| 国产精品乱码一区二三区的特点| 欧美成人一区二区免费高清观看| 国产午夜福利久久久久久| 欧美xxxx性猛交bbbb| 成人三级黄色视频| 欧美另类亚洲清纯唯美| 亚洲av男天堂| 天堂av国产一区二区熟女人妻| 国产黄色小视频在线观看| 亚洲成人av在线免费| 国产精品一区www在线观看| 人人妻人人澡欧美一区二区| 18禁动态无遮挡网站|