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

    霍亂弧菌新發(fā)菌株中vipA,vipB和clpV對(duì)T6SS轉(zhuǎn)錄調(diào)控

    2022-02-28 12:21:10吳文軒張幸鼎趙文婧
    關(guān)鍵詞:硫胺素霍亂弧菌基因簇

    吳文軒,郭 揚(yáng),張幸鼎,趙文婧

    (中山大學(xué)醫(yī)學(xué)院,廣東深圳 518107)

    霍亂是一種烈性傳染病,2010 年海地爆發(fā)霍亂疫情,相對(duì)于El Tor 經(jīng)典株海地變種El Tor 霍亂弧菌不僅在臨床上表現(xiàn)出更快的傳播力和更高的死亡率,也在動(dòng)物模型上表現(xiàn)出更強(qiáng)的定植力和毒力[1-2]。已有研究表明,六型分泌系統(tǒng)(Type Ⅵsecretion system,T6SS)對(duì)霍亂弧菌的生存和致病起著關(guān)鍵作用,例如其效應(yīng)因子VgrG-1 靶向真核生物,可以引起巨噬細(xì)胞或捕食性阿米巴細(xì)胞毒性肌動(dòng)蛋白交聯(lián),從而逃避捕食[3-5]。霍亂弧菌T6SS 基因簇由一個(gè)編碼結(jié)構(gòu)蛋白的主基因簇(VC0105-VC0124)和至少3 個(gè)附軸(VCA0017-VCA0021、VCA1415-VCA1419 和 VCA0284-VCA0286)構(gòu)成[6]。有研究報(bào)道通過(guò)轉(zhuǎn)錄組測(cè)序的方法,將第七次霍亂流行的El Tor 菌株N16961、C6706 和海地變異菌株H1 的T6SS 轉(zhuǎn)錄水平進(jìn)行比較,H1 株T6SS有著更高的轉(zhuǎn)錄水平;尤其是T6SS 鞘蛋白(VipA,VipB)和鞘水解蛋白(ClpV)[5]。我們猜測(cè)基因vipA,vipB和clpV中是否存在轉(zhuǎn)錄調(diào)控因子調(diào)控T6SS轉(zhuǎn)錄。雖然霍亂弧菌是研究T6SS 調(diào)控的模式菌株,但是霍亂弧菌會(huì)隨著環(huán)境的變化會(huì)快速產(chǎn)生適應(yīng)性突變。在2016 也門爆發(fā)近代史上最大的霍亂疫情,也門霍亂弧菌同樣為El Tor 霍亂弧菌的變種,同樣表現(xiàn)出更強(qiáng)的持續(xù)時(shí)間和嚴(yán)重程度[7-8]。探索新發(fā)疫情細(xì)菌H1 的T6SS 調(diào)控機(jī)制對(duì)于深刻了解T6SS 對(duì)于病原菌致病性,以及在細(xì)菌進(jìn)化中的作用,有著重要意義。為了探索霍亂弧菌新發(fā)菌株H1 的T6SS 轉(zhuǎn)錄調(diào)控機(jī)制,本研究以轉(zhuǎn)錄組測(cè)序?yàn)榛A(chǔ),分析H1 株分別由于vipA敲除(H1ΔvipA),vipB敲除(H1ΔvipB)和clpV敲除(H1ΔclpV)時(shí)T6SS基因簇的轉(zhuǎn)錄水平,并對(duì)其差異基因進(jìn)行了基因本體論(Gene ontology,GO)富集和京都基因與基因組百科全書(shū)(Kyoto Encyclopedia of Genes and Genomes,KEGG)富集。旨在分析H1 株中結(jié)構(gòu)相關(guān)基因與T6SS 轉(zhuǎn)錄調(diào)控的相關(guān)性,完善霍亂弧菌T6SS 調(diào)控圖譜;H1 株的T6SS 調(diào)控對(duì)將T6SS 作為靶標(biāo)來(lái)防治霍亂弧菌等病原菌的傳害提供重要理論基礎(chǔ)。

    1 材料與方法

    1.1 實(shí)驗(yàn)菌株

    根據(jù)實(shí)驗(yàn)需求將霍亂弧菌H1 株wt 和突變株培養(yǎng)于添加終濃度為100 μg/mL 的鏈霉素的LB 培養(yǎng)基,37 ℃,1.4×g的搖床中培養(yǎng);添加100 μg/mL的鏈霉素和1%瓊脂粉的LB 培養(yǎng)基,37 ℃靜置培養(yǎng)。實(shí)驗(yàn)操作嚴(yán)格遵守《生物安全保障承諾書(shū)》中的相關(guān)規(guī)定,按照標(biāo)準(zhǔn)操作展開(kāi)實(shí)驗(yàn)活動(dòng),確保實(shí)驗(yàn)室安全。

    1.2 主要試劑

    載體(pWM91 Vector,實(shí)驗(yàn)室自存);限制性內(nèi)切酶(NotI-HF 和SacI,NEB 公司);連接酶(T4 DNA Ligase,NEB 公司);Taq 酶(PrimeSTAR Max Premix,Takara 公司);感受態(tài)(SM10 和DH5α λ pir,實(shí)驗(yàn)室自存);TRIzol 試劑(TRIzol Reagent,Thermo公司);RNA 提取試劑盒(Pure Link RNA Mini Kit,Invitrogen 公司);DNA 清除試劑盒(Ambionogen-RNA Mini Kit,Invitrogen 公司);RNA-seq 建庫(kù)試劑盒(Next? Ultra Ⅱ?Directional RNA Library Prep Kit for Illumina,NEB 公司);反轉(zhuǎn)錄試劑盒(Evo MMLV RT Kit with gDNA Clean for qPCR,艾科瑞公司);熒光定量試劑盒(SYBR? Green Premix Pro Taq HS qPCR Kit,艾科瑞公司)。

    1.3 突變菌株構(gòu)建

    選取目標(biāo)基因上下游各800 bp 的序列分別進(jìn)行擴(kuò)增,將上下游片段的黏性末端相連組成1 600 bp片段克隆到pWM91自殺質(zhì)粒中。將重組質(zhì)粒轉(zhuǎn)化至SM10 感受態(tài)中,并與H1 wt 于37 ℃共培養(yǎng)2 h;保證重組質(zhì)粒整合至H1菌株的基因組中完成第1 次同源重組。將完成第一次同源重組的H1 菌株在含有6%蔗糖的平板上涂布,由于重組質(zhì)粒攜帶sacB基因。迫于蔗糖的選擇壓力,此時(shí)發(fā)生第2 次同源重組,目標(biāo)基因被敲除突變菌株構(gòu)建完成。按照以上方法在H1 wt基礎(chǔ)上構(gòu)建H1ΔvipA,H1ΔvipB和H1ΔclpV。

    1.4 細(xì)菌總RNA提取

    于37 ℃培養(yǎng)過(guò)夜,取109CFU 細(xì)菌于1 mL 的TRIzol 中,加入250 mg 的玻璃砂,勻漿2 次每次1 min。加入400 μL 的氯仿于樣品中混勻,在4 ℃預(yù)冷離心機(jī)中13 200 ×g離心30 min。取水相上清,按照Pure Link RNA Mini Kit 說(shuō)明書(shū)提取總RNA,隨后根據(jù)Ambionink RNA Mini Kit 的說(shuō)明書(shū)去除樣品中的DNA,得到純化的總RNA 樣品。利用Nanodrop One,1%瓊脂糖電泳(120V,30 min)和Agilent 4200檢測(cè)RNA 樣品的完整性、純度和總量,合格者進(jìn)行接下來(lái)的測(cè)序分析。每個(gè)菌株重復(fù)3 次,作為3次生物學(xué)重復(fù)。

    1.5 轉(zhuǎn)錄組建庫(kù)、測(cè)序及數(shù)據(jù)分析

    將合格的RNA 樣品按照Next? Ultra Ⅱ?Directional RNA Library Prep Kit for Illumina 的說(shuō)明書(shū)進(jìn)行cDNA 雙鏈合成和測(cè)序文庫(kù)構(gòu)建,利用Agilent 2100 對(duì)文庫(kù)質(zhì)量進(jìn)行評(píng)估;合格的文庫(kù)通過(guò)llumina Hiseq 平臺(tái)進(jìn)行RNA-seq 測(cè)序得到RNA-seq 原始數(shù)據(jù)。CLC Workbench 12.0 對(duì)數(shù)據(jù)進(jìn)行質(zhì)量檢測(cè)和過(guò)濾。將過(guò)濾后的數(shù)據(jù)與從NCBI(https://www.ncbi.nlm.nih.gov/)下載H1 wt 菌株的參考基因組文件進(jìn)行比對(duì),得到測(cè)序片段計(jì)數(shù)矩陣并運(yùn)用TMM 的方法進(jìn)行均一化。最后以H1 wt 菌株轉(zhuǎn)錄組作為內(nèi)參獲得差異基因的表達(dá)矩陣,并用BH 方法對(duì)P-value 進(jìn)行矯正;以FDR<0.01 且|log2FC|≥1為標(biāo)準(zhǔn)定義差異基因。

    1.6 逆轉(zhuǎn)錄和熒光定量PCR檢測(cè)

    利用熒光定量PCR 方法,對(duì)T6SS 基因簇中部分基因的RNA-seq 結(jié)果進(jìn)行qPCR 驗(yàn)證。RNA 樣品1 μg 按照Evo M-MLV RT Kit with gDNA Clean for qPCR 說(shuō)明書(shū)進(jìn)行逆轉(zhuǎn)錄得到cDNA。實(shí)時(shí)定量PCR 反應(yīng)體系如下:2×SYBR? Green Pro Taq HS Premix 10 μL、cDNA 2 μL、上下游特異性引物(10 μmol/L)各0.4 μL、ROX Reference Dye 0.4 μL、RNase free water 6.8 μL,引物序列均在NCBI 中驗(yàn)證。擴(kuò)增程序:95 ℃運(yùn)行30 s循環(huán)1次,95 ℃運(yùn)行5 s,60 ℃運(yùn)行30 s循環(huán)40次。每個(gè)樣本重復(fù)測(cè)3次,以gyrB作為管家基因,以H1 wt的CT值作為樣本間標(biāo)準(zhǔn)化內(nèi)參,使用log2FC計(jì)算基因的相對(duì)表達(dá)量。

    1.7 統(tǒng)計(jì)學(xué)方法

    使用Hiplot(https://hiplot.com.cn)對(duì)RNA-seq的差異基因進(jìn)行聚類分析、可視化分析和交叉關(guān)系展示。將H1 wt 基因組序列與eggNOG(http://eggnog-mapper.embl.de/)數(shù)據(jù)庫(kù)比對(duì),得到有GO和KO功能注釋的差異基因,并使用clusterProfiler R 進(jìn)行GO 和KEGG 富集。采用GraphPad Prism 8 統(tǒng)計(jì)軟件進(jìn)行數(shù)據(jù)分析,計(jì)量資料以均數(shù)±標(biāo)準(zhǔn)誤(Mean±SEM)表示。

    2 結(jié)果

    2.1 H1ΔvipA,H1ΔvipB 和H1ΔclpV 的T6SS 基因簇轉(zhuǎn)錄水平不同

    利用RNA-seq 技術(shù),將H1 wt 轉(zhuǎn)錄組作為樣本間標(biāo)準(zhǔn)化的內(nèi) 參,得 到H1ΔvipA,H1ΔvipB和H1ΔclpV的T6SS 基因簇log2FC 趨勢(shì)圖(圖1)。以|log2FC|≥1 為基準(zhǔn),認(rèn)為有差異。由圖1 可知,H1ΔvipA株T6SS 主 軸vipB(log2FC=4.18)、hsiF(log2FC=4.28)、vasA(log2FC=3.59)、vasB(log2FC=3.71)、vasC(log2FC=3.91)、vasD(log2FC=3.81)、vasE(log2FC=3.41)、vasF(log2FC=3.58)、clpV(log2FC=3.39)、vasH(log2FC=3.38)、vasI(log2FC=3.71)、vasJ(log2FC=3.40)、vasK(log2FC=3.31)、vasL(log2FC=2.60)和vgrG-3(log2FC=1.51),以及T6SS 附軸hcp-1(log2FC=3.59)、vgrG-1(log2FC=1.62)、hcp-2(log2FC=3.45)和vgrG-2(log2FC=2.21)的轉(zhuǎn)錄水平都呈增高趨勢(shì);僅vipA(log2FC=-11.78)轉(zhuǎn)錄水平由于敲除呈降低趨勢(shì)。H1ΔvipB株tsiV-3(log2FC=-1.23)、VC1420(log2FC=-1.06)和vipB(log2FC=-9.02)的轉(zhuǎn)錄水平呈降低趨勢(shì),其余基因|log2FC|均小于1。H1ΔclpV株tsiH(log2FC=1.32)轉(zhuǎn)錄水平呈增高趨勢(shì);僅clpV(log2FC=-8.45)轉(zhuǎn)錄水平由于敲除呈降低趨勢(shì)。

    圖1 H1ΔvipA,H1ΔvipB和H1ΔclpV株的T6SS基因簇轉(zhuǎn)錄水平Fig.1 The expression level of genes of T6SS clusters across H1ΔvipA,H1ΔvipB and H1ΔclpV

    2.2 H1ΔvipA,H1ΔvipB和H1ΔclpV轉(zhuǎn)錄組分析

    以H1 wt 轉(zhuǎn)錄組為樣本間標(biāo)準(zhǔn)化的內(nèi)參,尋找由于vipA,vipB和clpV基因敲除引起的轉(zhuǎn)錄差異。以FDR<0.01 且|log2FC|≥1 為標(biāo)準(zhǔn)定義差異表達(dá)基因,并且對(duì)差異基因進(jìn)行聚類分析(圖2A)。H1ΔvipA的差異基因共有22 個(gè),其中19 個(gè)基因轉(zhuǎn)錄水平上調(diào),3 個(gè)基因轉(zhuǎn)錄水平下調(diào)(圖2B);H1ΔvipB的差異基因有44 個(gè),其中17 個(gè)基因轉(zhuǎn)錄水平上調(diào),27 個(gè)基因轉(zhuǎn)錄水平下調(diào)(圖2C);H1ΔclpV的差異基因有86 個(gè),其中49 個(gè)基因轉(zhuǎn)錄水平上調(diào),35 個(gè)基因轉(zhuǎn)錄水平下調(diào)(圖2D)。并且H1ΔvipB和H1ΔclpV擁有相同的23 個(gè)差異基因;H1ΔvipB和H1ΔvipA有4 個(gè)相同的差異基因(hapA,hapR,vipB,KW3_RS0115980);H1ΔvipA和H1ΔclpV有4 個(gè)相同的差異基因(hapA,hapR,clpV,KW3_RS0115980);H1ΔvipA,H1ΔvipB和H1ΔclpV共同擁有的差異基因有3 個(gè)(hapA,hapR,KW3_RS0 115980;圖2E)。

    圖2 H1ΔvipA,H1ΔvipB和H1ΔclpV差異基因分析Fig.2 Differentially expressed genes across H1ΔvipA,H1ΔvipB and H1ΔclpV

    2.3 H1ΔvipA,H1ΔvipB 和H1ΔclpV 的差異表達(dá)基因的GO和KEGG分析

    以H1 wt 的GO 注釋作為基因背景,對(duì)H1ΔvipA,H1ΔvipB和H1ΔclpV的差異基因進(jìn)行GO 富集。H1ΔvipA的差異基因GO 富集通路主要涉及細(xì)胞外區(qū)域(extracellular region)、多機(jī)體過(guò)程(multi-organism process)、發(fā)病機(jī)制(pathogenesis)和物種間相互作用(interspecies interaction)等(圖3A)。H1ΔvipB的差異基因GO 富集通路主要涉及二甲基酸運(yùn)輸(dicarboxylic acid transport)、C4-二甲基運(yùn)輸(C4-dicarboxylate transport)、有機(jī)酸運(yùn)輸(organic acid transport)和有機(jī)離子運(yùn)輸(organic anion transport)等(圖3B)。H1ΔclpV的差異基因GO 富集通路主要涉及結(jié)構(gòu)分子活性(structural molecule activity)、核糖體的結(jié)構(gòu)成分(structural constituent of ribosome)、肽代謝過(guò)程(peptide metabolic process)和肽生物合成過(guò)程(peptide biosynthetic process)等(圖3C)。

    圖3 H1ΔvipA,H1ΔvipB和H1ΔclpV差異表達(dá)基因的GO富集Fig.3 GO enrichment of DEGs across H1ΔvipA,H1ΔvipB and H1ΔclpV

    以H1 wt 的KEGG 注釋作為基因背景,對(duì)H1ΔvipA,H1ΔvipB和H1ΔclpV的差異基因進(jìn)行KEGG 富集。H1ΔvipA的差異基因KEGG 富集通路主要涉及霍亂弧菌生物膜形成(biofilm formation-Vibrio cholerae,圖4A)。H1ΔvipB的差異基因KEGG富集通路主要涉及布塔諾酸鹽代謝(butanoate metabolism)、雙組件系統(tǒng)(two-component system)、霍亂弧菌生物膜形成(biofilm formation-Vibrio cholerae)和酮體的合成和降解等(synthesis and degradation of ketone bodies,圖4B)。H1ΔclpV的差異基因KEGG 富集通路主要涉及硫銨代謝(thiamine metabolism)等(圖4C)。

    圖4 H1ΔvipA,H1ΔvipB和H1ΔclpV差異表達(dá)基因的KEGG富集Fig.4 KEGG enrichment of DEGs across H1ΔvipA,H1ΔvipB and H1ΔclpV

    2.4 熒光定量PCR驗(yàn)證

    挑選T6SS 基因簇中部分基因以及hapR對(duì)H1ΔvipA,H1ΔvipB和H1ΔclpV的轉(zhuǎn)錄組結(jié)果進(jìn)行qPCR 驗(yàn) 證。H1ΔvipA的決定性系數(shù)為0.801 3;H1ΔvipB的決定性系數(shù)為0.816 5;H1ΔclpV的決定性系數(shù)為0.800 0,各個(gè)菌株T6SS 部分基因和hapR轉(zhuǎn)錄水平變化趨勢(shì)與RNA-seq結(jié)果基本一致(圖5)。

    3 討論

    為了探究霍亂弧菌H1 新發(fā)株T6SS 結(jié)構(gòu)蛋白編碼基因是否調(diào)控T6SS 轉(zhuǎn)錄,我們構(gòu)建了多個(gè)T6SS 敲除株。通過(guò)H1ΔvipA,H1ΔvipB和H1ΔclpV的T6SS 轉(zhuǎn)錄水平比 較,發(fā)現(xiàn) 僅vipA缺 失,H1 株T6SS 下游基因整體轉(zhuǎn)錄水平顯著性增高,包括下游的vipB和clpV;與其共編碼VipA/VipB 鞘蛋白復(fù)合體的vipB缺失和編碼鞘蛋白水解酶的clpV缺失,T6SS 轉(zhuǎn)錄水平卻無(wú)整體變化。同時(shí)RNA-seq 的結(jié)果表示H1ΔvipA有22 個(gè)差異基因,其中有18 個(gè)屬于T6SS 基因簇(圖2B)。已有研究表明,霍亂弧菌的鞘蛋白是高度保守的,其同源物已在所有T6SS結(jié)構(gòu)中發(fā)現(xiàn);并且VipA 和VipB 蛋白相結(jié)合是保證T6SS 活性的必要條件。在霍亂弧菌O1 V1552 中,vipA缺失會(huì)導(dǎo)致VipB 不能穩(wěn)定產(chǎn)生[9]。在同樣擁有T6SS 的土拉弗朗西斯菌(Francisella tularensis)中,vipA同源物基因的缺失,使vipB同源基因的轉(zhuǎn)錄水平下降兩倍,而vipB同源基因的缺失對(duì)vipA同源基因的轉(zhuǎn)錄水平并無(wú)影響[10]。以往研究表明vipA與vipB轉(zhuǎn)錄水平應(yīng)當(dāng)是相對(duì)穩(wěn)定的,而將vipA缺失vipB水平卻增高400 多倍(圖5),這種結(jié)果是從未報(bào)道過(guò)的。

    圖5 H1ΔvipA,H1ΔvipB和H1ΔclpV T6SS基因簇RNA-seq數(shù)據(jù)的qPCR驗(yàn)證Fig.5 Verification of RNA-seq data of genes of T6SS clusters by qPCR across H1ΔvipA,H1ΔvipB and H1ΔclpV

    H1ΔvipA,H1ΔvipB和H1ΔclpV共同擁有的差異基因有3 個(gè)(hapR,hapA,KW3_RS0115980)(圖2E)。在霍亂弧菌中,HapR是一種群體感應(yīng)通路主要的調(diào)節(jié)因子,可通過(guò)抑制霍亂弧菌多糖合成的調(diào)控基因vpsT及vpsR的表達(dá),來(lái)抑制霍亂弧菌生物膜形成和毒力基因的表達(dá);同時(shí)HapR 也可與T6SS附軸啟動(dòng)子結(jié)合,激活T6SS附軸的轉(zhuǎn)錄[11-13]。當(dāng)細(xì)胞密度較低時(shí),磷酸化的LuxO 激活Qrr sRNA 產(chǎn)生,激活A(yù)phA 的翻譯并抑制HapR 的翻譯,從而降低T6SS 附軸的轉(zhuǎn)錄水平;Qrr sRNA 也可通過(guò)與T6SS主基因簇的mRNA 堿基互補(bǔ)配對(duì)直接抑制T6SS 主軸轉(zhuǎn)錄[14]。當(dāng)細(xì)胞密度較高時(shí),LuxO 未磷酸化,Qrr sRNA 不轉(zhuǎn)錄,HapR 穩(wěn)定表達(dá),此時(shí)T6SS 主軸中的調(diào)控因子vasH 與HapR 一起促進(jìn)T6SS 附軸轉(zhuǎn)錄[11]。結(jié)合我們的實(shí)驗(yàn)結(jié)果,H1ΔvipA株中hapR是轉(zhuǎn)錄差異基因,并呈現(xiàn)明顯下降趨勢(shì);而hapR在H1ΔvipB和H1ΔclpV中呈明顯增高趨勢(shì)(圖2B-D)。這提示我們,hapR可能參與H1 株中vipA對(duì)T6SS 基因簇的轉(zhuǎn)錄調(diào)控,并與vipA轉(zhuǎn)錄趨勢(shì)一致。同時(shí),vipA可以調(diào)控T6SS 轉(zhuǎn)錄水平可能由于vipA與其他T6SS主軸基因共用一個(gè)啟動(dòng)子[6],并且vipA是T6SS主軸啟動(dòng)子后的第一個(gè)基因,可能作為反式作用因子結(jié)合阻遏蛋白,從而抑制后續(xù)T6SS基因簇轉(zhuǎn)錄。

    HapA 是hapA編碼的金屬蛋白酶,是一種霍亂弧菌O1 分泌的輔助毒素。其通過(guò)水解宿主黏蛋白、降低宿主細(xì)胞間的緊密連接、水解宿主的免疫蛋白等放方式加速疾病進(jìn)程[15]。在霍亂弧菌中,hapA的轉(zhuǎn)錄由HapR和RpoS共同調(diào)控,在營(yíng)養(yǎng)受限或細(xì)胞密度較高時(shí),群體感應(yīng)調(diào)控蛋白HapR 穩(wěn)定表達(dá),從而直接激活hapA的轉(zhuǎn)錄[16]。本實(shí)驗(yàn)中hapA和hapR表達(dá)趨勢(shì)一致并且與野生型有顯著差異性,即 在H1ΔvipA中hapA(log2FC=-3.89)和hapR(log2FC=-8.00)表達(dá)共降低,H1ΔvipB中hapA(log2FC=2.46)和hapR(log2FC=3.16)表達(dá)共增高,且H1ΔclpV中hapA(log2FC=2.27)和hapR(log2FC=3.78)表達(dá)共增高,提示hapA的轉(zhuǎn)錄變化可能由hapR變化引起。而三個(gè)敲除株中KW3_RS0115980轉(zhuǎn)錄水平均顯著上調(diào)(圖2B-D),但其基因編碼蛋白目前僅預(yù)測(cè)為假想蛋白,其具體功能需要進(jìn)一步揭示。

    ClpV 是由clpV編碼的ATP 酶,可通過(guò)水解VipA/VipB 鞘蛋白復(fù)合體,達(dá)到移除已收縮過(guò)的T6SS以啟動(dòng)新的T6SS分泌過(guò)程[17-18]。但將clpV敲除后,有多達(dá)86 個(gè)基因轉(zhuǎn)錄水平發(fā)生顯著改變,這些差異基因GO 富集到蛋白質(zhì)代謝相關(guān)通路(圖3C),KEGG 通路包括硫胺代謝通路(圖4C)。維生素B1又名硫胺素,是參與碳水化合物和氨基酸代謝的酶的輔因子,幾乎所有的生命體都需要硫胺素[19]。大多數(shù)細(xì)菌可以從頭合成硫胺素,當(dāng)缺乏硫胺素合成基因缺乏時(shí)則需要通過(guò)轉(zhuǎn)運(yùn)蛋白從外界攝取硫胺素[20]。目前細(xì)菌已經(jīng)進(jìn)化出許多硫胺素及前體轉(zhuǎn)運(yùn)蛋白,ABC 轉(zhuǎn)運(yùn)蛋白(ThiBPQ、ThiXYZ、ECFThiT和YkoEDC)、二級(jí)轉(zhuǎn)運(yùn)蛋白(NiaP)和促進(jìn)轉(zhuǎn)運(yùn)蛋白(PnuT)。在H1ΔclpV中thiC、thiF、thiH和thiS呈明顯轉(zhuǎn)錄增高趨勢(shì)(圖2D)。thiC在硫胺素合成中起著重要作用,其編碼蛋白ThiC 將硫胺素和嘌呤合成的分支點(diǎn)代謝物5-氨基咪唑核糖苷轉(zhuǎn)化為磷酸化的4-氨基-5-羥甲基-2-甲基嘧啶,這是合成有活性的焦磷酸硫胺素的關(guān)鍵[21]。有研究表示,在細(xì)菌、真菌和植物中硫胺素及其衍生物可以和RNA 結(jié)合以控制其基因的轉(zhuǎn)錄或翻譯[22]。可見(jiàn)硫胺素在轉(zhuǎn)錄調(diào)節(jié)方面有著重要作用,但至于在霍亂弧菌中是否也有基因轉(zhuǎn)錄調(diào)控的功能以及clpV缺失引起硫胺素合成相關(guān)基因轉(zhuǎn)錄水平升高,需要進(jìn)一步的實(shí)驗(yàn)探索。

    生物膜提供霍亂弧菌穩(wěn)定的生存環(huán)境,自然狀態(tài)下霍亂弧菌可以在浮游,單層生物膜中和成熟生物膜的生存狀態(tài)中切換,在不同的生存狀態(tài)下其體內(nèi)轉(zhuǎn)錄情況不同[23]。H1ΔvipA,H1ΔvipB和H1ΔclpV的KEGG 均富集到霍亂弧菌生物膜形成相關(guān)途徑(圖4A-C),其中H1ΔvipA的生物膜形成相關(guān)差異基因有hapA、hapR、KW3_RS0 107715;H1ΔvipB的生物膜形成相關(guān)差異基因有hapA、hapR、tcpS、KW3_RS0103100 和KW3_RS0 103115;H1ΔclpV的生物膜形成相關(guān)基因有hapA、KW3_RS0116000、hapR、KW3_RS0103100 和KW3_RS0103115 等?;魜y弧菌生物膜形成過(guò)程受多種途徑調(diào)控,其中包括群體感應(yīng)途徑[24]。已有研究表明霍亂弧菌C6706中hapR的表達(dá)時(shí)間與生物膜厚度,生物膜脫離率和霍亂弧菌定植率密切相關(guān)[25]。H1ΔvipA,H1ΔvipB和H1ΔclpV的KEGG 富集不僅富集到霍亂弧菌生物膜形成通路還都富集到群體感應(yīng)通路,推測(cè)H1株中群體感應(yīng)調(diào)控生物膜形成相關(guān)基因的轉(zhuǎn)錄水平變化。

    綜上所述,本研究發(fā)現(xiàn):vipA基因是一個(gè)存在與H1 菌株中具有多重功能的基因,一是編碼T6SS鞘結(jié)構(gòu)蛋白,二是參與T6SS 調(diào)控。同時(shí),研究表明hapR可能參與H1 菌株vipA對(duì)T6SS 的轉(zhuǎn)錄調(diào)控。結(jié)合前期研究結(jié)果,將存在于H1 菌株中的可能調(diào)控通路繪制簡(jiǎn)圖(圖6)。本研究為將霍亂弧菌新發(fā)菌株H1 的T6SS 作為防治靶標(biāo)奠定基礎(chǔ);也為霍亂弧菌T6SS的轉(zhuǎn)錄調(diào)控提供新的思路。

    圖6 霍亂弧菌H1株T6SS轉(zhuǎn)錄調(diào)控網(wǎng)絡(luò)Fig.6 T6SS transcription regulation network in Vibrio cholerae H1

    猜你喜歡
    硫胺素霍亂弧菌基因簇
    雛雞維生素B1缺乏癥的診治
    冬瓜高通量轉(zhuǎn)錄組測(cè)序及分析
    霍亂弧菌生物被膜發(fā)育與環(huán)境調(diào)控
    糖尿病酮癥酸中毒患者血清硫胺素水平改變及臨床意義
    血清凝集、基因測(cè)序聯(lián)合檢測(cè)群霍亂弧菌的應(yīng)用
    腸球菌萬(wàn)古霉素耐藥基因簇遺傳特性
    遺傳(2015年5期)2015-02-04 03:06:55
    海洋稀有放線菌 Salinispora arenicola CNP193 基因組新穎PKS 和NRPS基因簇的發(fā)掘
    反芻動(dòng)物硫胺素營(yíng)養(yǎng)研究進(jìn)展
    草食家畜(2014年6期)2014-03-27 01:21:16
    非O1群、非O139群霍亂弧菌致敗血癥1例報(bào)道
    三重實(shí)時(shí)熒光PCR檢測(cè)純培養(yǎng)物和環(huán)境水體中霍亂弧菌方法的建立及應(yīng)用*
    你懂的网址亚洲精品在线观看| 美女中出高潮动态图| 曰老女人黄片| 制服人妻中文乱码| 国产伦精品一区二区三区视频9| 成年美女黄网站色视频大全免费 | 一区二区三区精品91| 国产成人精品福利久久| 亚洲三级黄色毛片| 蜜桃久久精品国产亚洲av| 99热全是精品| 精品少妇久久久久久888优播| 777米奇影视久久| 少妇丰满av| 在线观看免费日韩欧美大片 | 中文字幕最新亚洲高清| www.色视频.com| 草草在线视频免费看| 少妇人妻 视频| 欧美bdsm另类| 大香蕉久久成人网| 黄片无遮挡物在线观看| .国产精品久久| 欧美变态另类bdsm刘玥| 亚洲人成网站在线播| 五月天丁香电影| 人人妻人人添人人爽欧美一区卜| 日本vs欧美在线观看视频| 久久99一区二区三区| 欧美人与善性xxx| 国产精品麻豆人妻色哟哟久久| 久久精品夜色国产| 久久99热这里只频精品6学生| 九九久久精品国产亚洲av麻豆| 国产精品 国内视频| 成年女人在线观看亚洲视频| 黄色视频在线播放观看不卡| 卡戴珊不雅视频在线播放| 亚洲av福利一区| 亚洲天堂av无毛| 制服丝袜香蕉在线| 欧美激情 高清一区二区三区| 狠狠精品人妻久久久久久综合| 美女cb高潮喷水在线观看| 街头女战士在线观看网站| 免费观看性生交大片5| 国产免费福利视频在线观看| 丰满饥渴人妻一区二区三| 一级片'在线观看视频| 亚洲精品视频女| 不卡视频在线观看欧美| 夜夜爽夜夜爽视频| 午夜精品国产一区二区电影| 下体分泌物呈黄色| 精品少妇内射三级| 国产视频首页在线观看| 日韩制服骚丝袜av| 在线观看免费日韩欧美大片 | 色视频在线一区二区三区| 啦啦啦在线观看免费高清www| 婷婷色综合大香蕉| 视频在线观看一区二区三区| 你懂的网址亚洲精品在线观看| 日本黄色日本黄色录像| 国产色婷婷99| 亚洲精品国产av蜜桃| 中文乱码字字幕精品一区二区三区| 精品少妇久久久久久888优播| 多毛熟女@视频| 十分钟在线观看高清视频www| 成人黄色视频免费在线看| 满18在线观看网站| 国产熟女午夜一区二区三区 | 精品少妇久久久久久888优播| 欧美亚洲 丝袜 人妻 在线| 国产在线视频一区二区| 国产在线一区二区三区精| 美女cb高潮喷水在线观看| 大码成人一级视频| 边亲边吃奶的免费视频| 中国美白少妇内射xxxbb| 亚洲精品一二三| 成人二区视频| 18禁在线无遮挡免费观看视频| 亚洲熟女精品中文字幕| 岛国毛片在线播放| 日韩视频在线欧美| 免费黄色在线免费观看| 国产综合精华液| 日韩伦理黄色片| 色网站视频免费| 免费黄频网站在线观看国产| 飞空精品影院首页| 另类精品久久| 晚上一个人看的免费电影| 妹子高潮喷水视频| 欧美丝袜亚洲另类| 在线观看www视频免费| 啦啦啦啦在线视频资源| 日本黄色片子视频| 人妻系列 视频| 韩国高清视频一区二区三区| 99热这里只有是精品在线观看| 亚洲成人一二三区av| 成人二区视频| 99热这里只有精品一区| 熟妇人妻不卡中文字幕| 蜜桃久久精品国产亚洲av| 高清毛片免费看| 亚洲精品乱码久久久v下载方式| 亚洲国产精品国产精品| 精品亚洲乱码少妇综合久久| 91成人精品电影| 又黄又爽又刺激的免费视频.| 国产高清三级在线| 丝袜喷水一区| 欧美日韩亚洲高清精品| 特大巨黑吊av在线直播| 欧美日韩亚洲高清精品| 亚洲精品乱码久久久v下载方式| 在线观看人妻少妇| 国产高清三级在线| 大又大粗又爽又黄少妇毛片口| 久久毛片免费看一区二区三区| 日韩熟女老妇一区二区性免费视频| 高清午夜精品一区二区三区| 一级,二级,三级黄色视频| 日韩av不卡免费在线播放| 亚洲精品亚洲一区二区| 欧美bdsm另类| 国产精品一国产av| 亚洲精品久久久久久婷婷小说| 国产精品 国内视频| 视频在线观看一区二区三区| 亚洲精品视频女| 精品人妻偷拍中文字幕| 男女啪啪激烈高潮av片| 国产黄色免费在线视频| 国产永久视频网站| 国产精品 国内视频| 亚洲精品日韩在线中文字幕| 精品人妻偷拍中文字幕| 国产免费一级a男人的天堂| 人妻夜夜爽99麻豆av| 国产成人精品久久久久久| 亚洲第一区二区三区不卡| 免费观看av网站的网址| 满18在线观看网站| 午夜福利网站1000一区二区三区| 五月开心婷婷网| 亚洲国产日韩一区二区| 岛国毛片在线播放| 大陆偷拍与自拍| 大话2 男鬼变身卡| 在线看a的网站| 边亲边吃奶的免费视频| 亚洲av欧美aⅴ国产| 久久久亚洲精品成人影院| 亚洲av欧美aⅴ国产| 亚洲精品一二三| 日韩成人伦理影院| 久久久久久久久久成人| 亚洲性久久影院| 国产探花极品一区二区| 日韩欧美一区视频在线观看| 一个人免费看片子| 少妇精品久久久久久久| 久久这里有精品视频免费| 少妇的逼水好多| 国产精品人妻久久久久久| 99国产精品免费福利视频| 大香蕉久久网| 在线天堂最新版资源| 亚洲第一区二区三区不卡| 18在线观看网站| 母亲3免费完整高清在线观看 | 亚洲精品久久久久久婷婷小说| 一区二区三区精品91| 高清av免费在线| 啦啦啦中文免费视频观看日本| 狂野欧美激情性xxxx在线观看| 成人影院久久| xxxhd国产人妻xxx| 久久精品国产自在天天线| 亚洲美女视频黄频| 国产成人精品婷婷| √禁漫天堂资源中文www| 美女xxoo啪啪120秒动态图| 久久久久久伊人网av| 亚洲美女视频黄频| 国产一区二区在线观看av| 精品国产一区二区久久| 欧美日韩视频精品一区| 日本av免费视频播放| 国产极品天堂在线| av播播在线观看一区| 久久久午夜欧美精品| 男女边摸边吃奶| 久久久久久久久久久久大奶| 亚洲经典国产精华液单| 亚洲精品美女久久av网站| 久久精品国产亚洲av涩爱| 亚洲在久久综合| 国产探花极品一区二区| 黄色配什么色好看| 岛国毛片在线播放| 国产成人精品在线电影| 亚洲欧美精品自产自拍| 国产熟女午夜一区二区三区 | 国产精品久久久久久精品古装| 少妇精品久久久久久久| 麻豆成人av视频| 九九在线视频观看精品| 国产黄频视频在线观看| 草草在线视频免费看| 色哟哟·www| av国产精品久久久久影院| 国内精品宾馆在线| 成人二区视频| 秋霞伦理黄片| 成年人免费黄色播放视频| 91久久精品电影网| 久久精品国产自在天天线| 一边亲一边摸免费视频| 一区二区三区免费毛片| 亚洲av不卡在线观看| 夜夜爽夜夜爽视频| 看免费成人av毛片| 国产成人精品久久久久久| 一区二区三区精品91| 亚洲国产av新网站| 狂野欧美激情性xxxx在线观看| 春色校园在线视频观看| 妹子高潮喷水视频| 国产一区有黄有色的免费视频| 精品少妇黑人巨大在线播放| 成人毛片a级毛片在线播放| 亚洲欧美成人综合另类久久久| 少妇的逼好多水| 亚洲国产精品一区三区| 色婷婷久久久亚洲欧美| 纯流量卡能插随身wifi吗| 色5月婷婷丁香| 国产av码专区亚洲av| 最黄视频免费看| 黑人高潮一二区| 一级二级三级毛片免费看| 国模一区二区三区四区视频| 国产精品国产三级国产专区5o| 精品亚洲成国产av| 亚洲性久久影院| 永久网站在线| 久久精品国产自在天天线| 久久久久视频综合| 国产综合精华液| 国产老妇伦熟女老妇高清| 三上悠亚av全集在线观看| 人妻 亚洲 视频| 久久女婷五月综合色啪小说| 汤姆久久久久久久影院中文字幕| 我的女老师完整版在线观看| 天堂8中文在线网| 日本欧美国产在线视频| 999精品在线视频| 亚洲丝袜综合中文字幕| 日韩熟女老妇一区二区性免费视频| 久久久久久久久久久久大奶| 久久精品久久久久久噜噜老黄| a级毛片免费高清观看在线播放| 国产免费一区二区三区四区乱码| 久久久久国产精品人妻一区二区| 成年av动漫网址| 80岁老熟妇乱子伦牲交| 全区人妻精品视频| 久久精品国产自在天天线| a级毛片免费高清观看在线播放| 熟女av电影| 岛国毛片在线播放| 久久久午夜欧美精品| 亚洲一级一片aⅴ在线观看| 91在线精品国自产拍蜜月| 亚洲欧美一区二区三区国产| 蜜桃久久精品国产亚洲av| 人妻一区二区av| 亚洲欧美一区二区三区黑人 | 水蜜桃什么品种好| 亚洲国产精品国产精品| 亚洲av日韩在线播放| √禁漫天堂资源中文www| 丝袜美足系列| 久久99热6这里只有精品| 青春草视频在线免费观看| videosex国产| 2022亚洲国产成人精品| 亚洲av福利一区| 国产精品久久久久久av不卡| 夜夜骑夜夜射夜夜干| 国产欧美日韩一区二区三区在线 | 国产伦精品一区二区三区视频9| 青青草视频在线视频观看| 欧美xxxx性猛交bbbb| 天天影视国产精品| 亚洲av二区三区四区| 男人操女人黄网站| 亚洲内射少妇av| 午夜福利影视在线免费观看| 亚洲国产欧美日韩在线播放| 熟女av电影| 男女边摸边吃奶| 国产成人精品婷婷| 午夜福利影视在线免费观看| 老熟女久久久| 五月开心婷婷网| 69精品国产乱码久久久| 99九九线精品视频在线观看视频| 熟女人妻精品中文字幕| 亚洲三级黄色毛片| 黄色配什么色好看| 日韩人妻高清精品专区| 成人国产av品久久久| 在现免费观看毛片| 在线亚洲精品国产二区图片欧美 | 国产伦精品一区二区三区视频9| 黑人猛操日本美女一级片| 一区二区三区乱码不卡18| 丝袜在线中文字幕| 午夜免费观看性视频| 久久国产精品大桥未久av| 亚洲精品,欧美精品| 国产成人精品婷婷| 午夜91福利影院| 久热久热在线精品观看| 久热这里只有精品99| 免费高清在线观看日韩| 日本av免费视频播放| 国产精品欧美亚洲77777| 亚洲精品第二区| 精品一品国产午夜福利视频| 男女边摸边吃奶| 男女国产视频网站| 一边摸一边做爽爽视频免费| 91精品一卡2卡3卡4卡| 国产黄色免费在线视频| av天堂久久9| 三级国产精品片| 亚洲精华国产精华液的使用体验| 亚洲一区二区三区欧美精品| 免费av中文字幕在线| 久久久久国产精品人妻一区二区| 久热这里只有精品99| 亚洲一区二区三区欧美精品| 亚洲国产精品一区三区| a级毛色黄片| 精品人妻熟女毛片av久久网站| 哪个播放器可以免费观看大片| av在线播放精品| 九九爱精品视频在线观看| 国产男女超爽视频在线观看| 成人国产麻豆网| av卡一久久| av在线老鸭窝| 国产女主播在线喷水免费视频网站| 在线观看三级黄色| 国产片特级美女逼逼视频| av专区在线播放| 亚洲av中文av极速乱| 一级毛片电影观看| 亚洲精品第二区| 在线精品无人区一区二区三| 亚洲天堂av无毛| 91精品伊人久久大香线蕉| 十八禁高潮呻吟视频| 国产免费福利视频在线观看| 又粗又硬又长又爽又黄的视频| 美女国产高潮福利片在线看| 久久久精品免费免费高清| 最新的欧美精品一区二区| 亚洲色图综合在线观看| 久热这里只有精品99| 九九久久精品国产亚洲av麻豆| 丝袜脚勾引网站| av有码第一页| 精品午夜福利在线看| 国产 精品1| 欧美激情 高清一区二区三区| 亚洲天堂av无毛| 99久久综合免费| 亚洲少妇的诱惑av| 夜夜看夜夜爽夜夜摸| 2018国产大陆天天弄谢| 亚洲五月色婷婷综合| 精品少妇黑人巨大在线播放| 国产国语露脸激情在线看| 丰满少妇做爰视频| 黄色欧美视频在线观看| 亚洲激情五月婷婷啪啪| 黑人欧美特级aaaaaa片| 国产欧美日韩一区二区三区在线 | 国产男女内射视频| 一边摸一边做爽爽视频免费| 在线精品无人区一区二区三| 国产国语露脸激情在线看| 一本大道久久a久久精品| h视频一区二区三区| 成人手机av| 国产成人a∨麻豆精品| 一级爰片在线观看| 亚洲精品乱码久久久久久按摩| 青青草视频在线视频观看| 国产深夜福利视频在线观看| 精品熟女少妇av免费看| 一级毛片电影观看| 亚洲精品中文字幕在线视频| 大话2 男鬼变身卡| 亚洲欧美日韩卡通动漫| 国产午夜精品一二区理论片| 高清黄色对白视频在线免费看| 成人影院久久| 成年av动漫网址| 三级国产精品欧美在线观看| 久久午夜综合久久蜜桃| 99热这里只有精品一区| 天堂中文最新版在线下载| 成人影院久久| 蜜桃久久精品国产亚洲av| 永久网站在线| 欧美人与性动交α欧美精品济南到 | 亚洲欧美一区二区三区黑人 | 两个人免费观看高清视频| 亚洲精品日韩av片在线观看| 99久久综合免费| 尾随美女入室| 69精品国产乱码久久久| 精品少妇内射三级| 国产高清三级在线| av网站免费在线观看视频| 国产精品一区二区在线不卡| 天堂俺去俺来也www色官网| 日韩,欧美,国产一区二区三区| 亚洲美女视频黄频| 免费观看的影片在线观看| 午夜免费观看性视频| 成人手机av| 大片免费播放器 马上看| 精品亚洲成国产av| 亚洲国产精品一区二区三区在线| 久久久精品94久久精品| 超碰97精品在线观看| 熟女电影av网| 视频中文字幕在线观看| 男人添女人高潮全过程视频| 26uuu在线亚洲综合色| 丰满少妇做爰视频| 国产精品国产av在线观看| 26uuu在线亚洲综合色| 丝袜脚勾引网站| av有码第一页| 男女边吃奶边做爰视频| 中文字幕亚洲精品专区| 男女免费视频国产| 色视频在线一区二区三区| 亚洲欧美色中文字幕在线| 亚洲熟女精品中文字幕| 在现免费观看毛片| 成年人免费黄色播放视频| 国产国拍精品亚洲av在线观看| 特大巨黑吊av在线直播| 边亲边吃奶的免费视频| 啦啦啦视频在线资源免费观看| 男女边吃奶边做爰视频| 欧美精品亚洲一区二区| 五月开心婷婷网| 一级毛片电影观看| 亚洲国产精品国产精品| 蜜臀久久99精品久久宅男| 夜夜骑夜夜射夜夜干| 欧美激情国产日韩精品一区| 日本猛色少妇xxxxx猛交久久| 国产日韩欧美在线精品| 国产欧美日韩综合在线一区二区| 成人二区视频| 国产黄色免费在线视频| 日韩亚洲欧美综合| 日韩熟女老妇一区二区性免费视频| 国产片内射在线| 国产一区有黄有色的免费视频| 人人妻人人澡人人看| 国产亚洲一区二区精品| 男的添女的下面高潮视频| 人妻人人澡人人爽人人| 免费看光身美女| 亚洲精品色激情综合| 国产精品一二三区在线看| 大香蕉久久网| 国产视频首页在线观看| 中文字幕久久专区| 最新中文字幕久久久久| 国产精品国产三级国产专区5o| 最近中文字幕2019免费版| 欧美性感艳星| 日韩精品免费视频一区二区三区 | 日韩成人伦理影院| 国产男女内射视频| 99久久人妻综合| 欧美日韩在线观看h| 黄色视频在线播放观看不卡| 色网站视频免费| 在线观看美女被高潮喷水网站| 亚洲欧美色中文字幕在线| 视频区图区小说| 妹子高潮喷水视频| 国产永久视频网站| 国产精品久久久久成人av| 国产一区二区三区综合在线观看 | 大片电影免费在线观看免费| 黑人巨大精品欧美一区二区蜜桃 | av在线观看视频网站免费| 女性被躁到高潮视频| 久久久国产精品麻豆| 久久人人爽av亚洲精品天堂| 日韩欧美一区视频在线观看| 一本—道久久a久久精品蜜桃钙片| 久久久久久久国产电影| 91国产中文字幕| 欧美xxxx性猛交bbbb| 寂寞人妻少妇视频99o| 日韩欧美精品免费久久| 九草在线视频观看| 亚洲综合色网址| 精品亚洲乱码少妇综合久久| 欧美人与性动交α欧美精品济南到 | 99久久精品国产国产毛片| 国产精品.久久久| 人人妻人人添人人爽欧美一区卜| 午夜免费观看性视频| av又黄又爽大尺度在线免费看| 高清在线视频一区二区三区| 亚洲精品日韩av片在线观看| 国产乱来视频区| 一级,二级,三级黄色视频| 国产老妇伦熟女老妇高清| 全区人妻精品视频| 国产成人a∨麻豆精品| 人妻制服诱惑在线中文字幕| 国产精品女同一区二区软件| 亚洲精华国产精华液的使用体验| 99九九线精品视频在线观看视频| 亚洲三级黄色毛片| 我的女老师完整版在线观看| 久久精品国产a三级三级三级| av福利片在线| 91精品国产九色| 国产精品无大码| 久久99热6这里只有精品| 亚洲av国产av综合av卡| 精品人妻熟女av久视频| 天天操日日干夜夜撸| 老司机影院成人| 寂寞人妻少妇视频99o| 夜夜骑夜夜射夜夜干| 99视频精品全部免费 在线| 王馨瑶露胸无遮挡在线观看| 国产精品.久久久| 久久 成人 亚洲| 国产在线免费精品| 亚洲精品av麻豆狂野| 免费大片黄手机在线观看| 久久久久久久精品精品| a级片在线免费高清观看视频| 欧美国产精品一级二级三级| 一区二区日韩欧美中文字幕 | 欧美xxⅹ黑人| av不卡在线播放| 制服丝袜香蕉在线| 老司机影院成人| 三级国产精品欧美在线观看| 免费日韩欧美在线观看| 寂寞人妻少妇视频99o| 色婷婷久久久亚洲欧美| 人妻夜夜爽99麻豆av| 亚洲国产色片| 国产免费一级a男人的天堂| av女优亚洲男人天堂| av网站免费在线观看视频| 亚洲欧美中文字幕日韩二区| 男的添女的下面高潮视频| 精品人妻熟女av久视频| 亚洲精品自拍成人| 少妇 在线观看| 亚洲色图 男人天堂 中文字幕 | 一级毛片我不卡| 亚洲av日韩在线播放| 国产 精品1| 男女边摸边吃奶| 青春草国产在线视频| 国产一区二区三区av在线| 婷婷成人精品国产| 久久精品国产自在天天线| 午夜精品国产一区二区电影| 少妇人妻 视频| videos熟女内射| 自拍欧美九色日韩亚洲蝌蚪91| 免费人成在线观看视频色| 97在线视频观看| 久久亚洲国产成人精品v| 人人澡人人妻人| 国产在线免费精品| 久久午夜福利片| 亚洲第一区二区三区不卡| 男人爽女人下面视频在线观看| 美女大奶头黄色视频| 免费黄色在线免费观看| 亚洲美女黄色视频免费看| 久久久久久久国产电影| 99热网站在线观看|