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

    瑞士乳桿菌調(diào)控小鼠腸道菌群變化規(guī)律的研究

    2018-01-08 02:47:25臧凱麗崔文靜馬新穎趙林森閻亞麗陳慶森
    食品科學(xué) 2018年1期
    關(guān)鍵詞:梭菌灌胃菌群

    臧凱麗,賈 彥,崔文靜,馬新穎,王 泳,趙林森,趙 培,葉 雷,閻亞麗,*,陳慶森,*

    (1.天津商業(yè)大學(xué)生物技術(shù)與食品科學(xué)學(xué)院,天津市食品生物技術(shù)重點(diǎn)實(shí)驗(yàn)室,天津 300134;2.河北一然生物科技有限公司,河北 石家莊 050899)

    瑞士乳桿菌調(diào)控小鼠腸道菌群變化規(guī)律的研究

    臧凱麗1,賈 彥1,崔文靜2,馬新穎2,王 泳1,趙林森2,趙 培1,葉 雷1,閻亞麗1,*,陳慶森1,*

    (1.天津商業(yè)大學(xué)生物技術(shù)與食品科學(xué)學(xué)院,天津市食品生物技術(shù)重點(diǎn)實(shí)驗(yàn)室,天津 300134;2.河北一然生物科技有限公司,河北 石家莊 050899)

    為探究瑞士乳桿菌(Lactobacillus helveticus)TS206對小鼠腸道菌群的調(diào)控規(guī)律,更好地闡述益生菌維持腸道菌群穩(wěn)態(tài)的機(jī)制,本研究以雄性BALB/c小鼠為實(shí)驗(yàn)對象,將其分為L. helveticus TS206組、對照組、空白組。L. helveticus TS206組和對照組分別每天灌胃L. helveticus TS206菌懸液和等量的生理鹽水,連續(xù)灌胃7 周,空白組不灌胃。每周采集小鼠新鮮糞便樣品,通過試劑盒提取糞便細(xì)菌總DNA,利用Ion torrent個(gè)人化操作基因組測序平臺測序技術(shù)對16S rRNA基因的V6區(qū)進(jìn)行高通量測序,最后通過生物信息學(xué)和多變量統(tǒng)計(jì)學(xué)方法對測序數(shù)據(jù)進(jìn)行分析。結(jié)果顯示:所有測序序列在97%相似水平劃分得到1 617 個(gè)操作分類單位(operational taxonomic units,OTUs),被劃分為8 個(gè)門,硬壁菌門(Firmicutes)和擬桿菌門(Bacteroidetes)為各組小鼠中的優(yōu)勢菌門,占總序列數(shù)的97.49%,且硬壁菌門豐度最高,超過70%;優(yōu)勢菌科主要為S24-7、毛螺菌科(Lachnospiraceae)、瘤胃球菌科(Ruminococcaceae)、理研菌科(Rikenellaceae)和擬桿菌科(Bacteroidaceae),且與對照組相比L. helveticus TS206組腸桿菌科細(xì)菌數(shù)量較低并呈下降趨勢。經(jīng)LEfSe分析有53 個(gè)關(guān)鍵OTUs與2 組小鼠腸道糞便菌群結(jié)構(gòu)顯著相關(guān),其中27 個(gè)OTUs在對照組中富集,分屬于腸桿菌科、瘤胃球菌科、毛螺菌科和梭菌目;26 個(gè)OTUs在L. helveticus TS206組中富集,分屬于梭菌目、毛螺菌科和擬桿菌屬,其中與對照組相比,梭菌目豐度水平高。經(jīng)主成分分析,兩組樣品完全分離,菌群的整體結(jié)構(gòu)存在差異,且富集的OTUs之間呈負(fù)相關(guān)關(guān)系。綜上所述,在實(shí)驗(yàn)階段,灌胃雄性BALB/c小鼠一定劑量的L. helveticus TS206后,經(jīng)擴(kuò)增子測序分析獲得的初步結(jié)論為L. helveticus TS206能改變腸道菌群整體結(jié)構(gòu)的功效,表現(xiàn)出了抑制腸道有害微生物生長的作用,并通過促進(jìn)部分有益菌的增殖來維持腸道菌群結(jié)構(gòu)的穩(wěn)態(tài)。

    瑞士乳桿菌TS206;腸道菌群;Ion torrent個(gè)人化操作基因組測序技術(shù);生物信息學(xué);多變量統(tǒng)計(jì)學(xué)分析

    Lederberg等[1]提出了人是一種由自身細(xì)胞和腸道微生物組成的“超級生物體”。人體腸道菌群編碼的基因組可被視為“人的第二基因組”,人的基因組與人體腸道微生物宏基因組相互協(xié)調(diào)、和諧一致,共同參與人體的營養(yǎng)攝入、代謝和免疫的調(diào)控過程。到目前為止,眾多的研究表明腸道菌群的結(jié)構(gòu)和功能與機(jī)體的健康狀況息息相關(guān),腸道菌群影響宿主的代謝顯型[2]、營養(yǎng)物質(zhì)的產(chǎn)生和吸收[3]以及調(diào)整和改善機(jī)體的免疫系統(tǒng)[4]。腸道內(nèi)菌群的任何動態(tài)平衡的打破都會影響胃腸道內(nèi)的正?;顒?,造成腸道內(nèi)的疾病,如炎性腸病、腸應(yīng)激綜合癥、過敏性反應(yīng)等[5],特別是與飲食不當(dāng)造成的“現(xiàn)代文明病”如肥胖、糖尿病、心腦血管疾病、結(jié)腸癌等發(fā)生發(fā)展具有密切的關(guān)系[6-9]。因此,腸道菌群結(jié)構(gòu)可以忠實(shí)而精細(xì)地反映人體的健康狀況,而益生菌可以通過改變腸道內(nèi)微生物結(jié)構(gòu),維持微生物區(qū)系的平衡來調(diào)控機(jī)體的代謝功能和腸黏膜屏障功能,這對腸道相關(guān)疾病的預(yù)防和治療具有重要意義[10]。目前,國內(nèi)外常用的益生菌菌種主要集中在功能特性比較清晰的乳桿菌和雙歧桿菌屬,如嗜酸乳桿菌(Lactobacillus acidophilus)、瑞士乳桿菌(Lactobacillus helveticus)、干酪乳桿菌(Lactobacillus casei)、植物乳桿菌(Lactobacillus plantarum)以及雙歧桿菌(Bif i dobacterium)等[11-12]。

    近年來,高通量測序技術(shù)是宏基因組學(xué)研究中應(yīng)用最廣泛的測序技術(shù)[13],主流的高通量測序技術(shù)主要包括454焦磷酸測序平臺、美國Illumina公司的Hiseq和Miseq測序平臺以及美國Life Technologies公司的Ion torrent個(gè)人化操作基因組測序儀(personal genome machine,PGM)和Ion proton測序平臺。基于16S rRNA基因可變區(qū)高通量測序的技術(shù)(擴(kuò)增子測序)被廣泛用于腸道微生物群落結(jié)構(gòu)分析中,為解析腸道菌群結(jié)構(gòu)的變化與人體健康的關(guān)系提供了很好的技術(shù)平臺。

    到目前為止,就L. helveticus而言,國內(nèi)外的研究主要集中在其發(fā)酵產(chǎn)生的生物活性物質(zhì)(比如血管緊張素轉(zhuǎn)化酶抑制肽等)對機(jī)體健康的促進(jìn)作用。有關(guān)L. helveticus調(diào)控小鼠腸道菌群為主題的國內(nèi)外文獻(xiàn)及研究報(bào)道較少。王友湘等[14]采用選擇性培養(yǎng)基探討了L. helveticus TS206干預(yù)對小鼠腸道菌群結(jié)構(gòu)的影響,發(fā)現(xiàn)小鼠腸道內(nèi)乳酸菌、雙歧桿菌、腸球菌等有益菌的數(shù)量顯著增加,提示L. helveticus TS206能夠促進(jìn)小鼠腸道內(nèi)有益菌的增殖,能夠維持小鼠腸道微生態(tài)系統(tǒng)的平衡。而Frece等[15]通過傳統(tǒng)的平板計(jì)數(shù)法檢測小鼠糞便中相關(guān)細(xì)菌的數(shù)量,發(fā)現(xiàn)灌胃L. helveticus 8 d后小鼠腸道中乳酸桿菌數(shù)量增加而腸桿菌以及梭菌屬細(xì)菌減少,這也與王友湘等[14]的研究結(jié)論相吻合。Taverniti等[16]指出,一些體外研究表明L. helveticus表現(xiàn)出許多常見的益生菌特性,如腸胃中生存能力、定植上皮細(xì)胞和抵抗病原體的能力,體內(nèi)小鼠模型研究表明L. helveticus可以預(yù)防胃腸道感染,增強(qiáng)其對病原體的抵抗力,調(diào)節(jié)宿主免疫反應(yīng),并影響腸道菌群的組成。

    鑒于針對L. helveticus的這類研究方法較少,所以本研究參考Singh等[17]給小鼠灌胃Bifidobacterium bifidum MIMBb75,利用定量聚合酶鏈?zhǔn)椒磻?yīng)(polymerase chain reaction,PCR)技術(shù)檢測小鼠腸道菌群變化的方法,以L. helveticus TS206為研究對象,利用擴(kuò)增子測序技術(shù)對長期L. helveticus TS206干預(yù)小鼠的腸道菌群結(jié)構(gòu)進(jìn)行分析,探究L. helveticus TS206調(diào)控小鼠腸道菌群結(jié)構(gòu)變化的規(guī)律并鑒定出與其干預(yù)顯著相關(guān)的關(guān)鍵菌屬。本研究將豐富L. helveticus TS206的益生功能,為該菌作為功能性食品改善人體健康提供可靠的理論依據(jù),還為以腸道菌群為靶點(diǎn)通過益生菌的攝入來預(yù)防各類慢性疾病帶來了全新的思路和方法。

    1 材料與方法

    1.1 材料與試劑

    L. helveticus TS206,天津商業(yè)大學(xué)實(shí)驗(yàn)室保藏。

    110 只6 周齡BALB/c雄性小鼠,SPF級,購自中國人民解放軍軍事醫(yī)學(xué)科學(xué)院實(shí)驗(yàn)動物中心,許可證號SCXK(京)2009-0017。小鼠飼料為普通繁育飼料,SPF級,購于中國人民解放軍軍事醫(yī)學(xué)科學(xué)院實(shí)驗(yàn)動物中心。

    MRS液體培養(yǎng)基 北京奧博星生物技術(shù)有限責(zé)任公司;乳清粉液體培養(yǎng)基(質(zhì)量分?jǐn)?shù)10%的乳清粉液體培養(yǎng)基)、腸道菌計(jì)數(shù)瓊脂(violet red bile dextrose agar,VRBDA)培養(yǎng)基 青島海博生物有限公司;QIAamp Fast DNA Stool Mini Kit 德國QIAGEN企業(yè);Pfu DNA聚合酶、Ion Plus Fragment Library Kit、Ion Xpress? Barcode Adapters Kit、Ion PGM?Template OT2 200 Kit、Ion PGM? Sequencing 200 Kit v2美國賽默飛世爾科技公司;MiniBEST Agarose Gel DNA Extraction Kit 寶生物工程(大連)有限公司;TAE緩沖液(50×) 上海生物工程有限公司;High Sensitivity DNA Kit 美國安捷倫公司; Qubit?dsDNA HS Assay Kit 美國Invitrogen公司;Agencourt?AMPure?XP Kit美國貝克曼公司。

    1.2 儀器與設(shè)備

    900 SERIES超低溫冰箱 美國賽默飛世爾科技公司;SW-CJ-2DD單人雙面凈化工作臺 蘇州凈化設(shè)備有限公司;3K18高速冷凍離心機(jī) 美國Sigma公司;DYY-2C電泳儀 北京六一儀器廠;Ion torrent個(gè)人化操作基因組測序平臺、INS1005527模板制備系統(tǒng)(Ion One TouchTM2)、8441-21模板富集系統(tǒng)(Ion OneTouch? ES)美國賽默飛世爾科技公司;331-1生物分析儀 美國安捷倫公司;N8050200熒光定量PCR儀 美國ABI公司;Qubit 2.0核酸蛋白定量儀、12321D磁力架 美國Invitrogen公司;II-3DNA紫外-可見分光光度計(jì) 英國柏諾公司。

    1.3 方法

    1.3.1 L. helveticus TS206的培養(yǎng)和富集

    將保藏的L. helveticus TS206于乳清粉液體培養(yǎng)基40 ℃活化2 代。將活化后的L. helveticus TS206液按體積分?jǐn)?shù)1%的接種量接種在MRS液體培養(yǎng)基中,40 ℃靜置培養(yǎng)18 h至對數(shù)期,5 000 r/min 4 ℃離心10 min收集菌體細(xì)胞,無菌生理鹽水洗滌2 次,最后用10 mL無菌生理鹽水懸浮使菌體濃度為3×108~4×108CFU/mL,以備后續(xù)動物實(shí)驗(yàn)灌胃使用。

    1.3.2 動物實(shí)驗(yàn)

    1.3.2.1 動物分組與灌胃

    110 只SPF級雄性BALB/c小鼠適應(yīng)性喂養(yǎng)1 周后隨機(jī)分為L. helveticus TS206組、對照組和空白組,其中L. helveticus TS206組45 只,對照組45 只,空白組20 只;L. helveticus TS206組小鼠每天灌胃0.2 mL L. helveticus TS206菌懸液,對照組小鼠每天灌胃等量的生理鹽水,連續(xù)灌胃7 周,灌胃結(jié)束后再繼續(xù)自由飲食飲水飼養(yǎng)1 周,空白組不灌胃;實(shí)驗(yàn)期間各組小鼠均飼喂普通繁育飼料,自由進(jìn)食和飲水,同時(shí)對小鼠的毛發(fā)色澤、飲食狀況、活動程度進(jìn)行記錄。

    1.3.2.2 小鼠糞便采集

    分別在實(shí)驗(yàn)開始后第0、1、2、3、4、5、6、7、8周早上8 點(diǎn),采用逼迫法采集L. helveticus TS206組和對照組小鼠新鮮糞便,在第2、4、6、8周采集空白組小鼠糞便,并立即放入冰盒保存后帶回實(shí)驗(yàn)室進(jìn)行后續(xù)實(shí)驗(yàn)。

    1.3.3 Ion torrent PGM測序技術(shù)探究L. helveticus TS206對小鼠腸道菌群結(jié)構(gòu)的影響

    1.3.3.1 小鼠糞便中細(xì)菌基因組DNA的提取

    迅速稱取-80 ℃凍存的1.3.2.2節(jié)所采集小鼠糞便樣品0.2 g,按照QIAamp Fast DNA Stool Mini Kit的使用說明提取細(xì)菌基因組DNA。用DNA紫外-可見分光光度計(jì)進(jìn)行核酸質(zhì)量濃度檢測,并用0.8%瓊脂糖凝膠電泳檢測基因組DNA的完整性。

    1.3.3.2 腸道細(xì)菌16S rRNA V6區(qū)的PCR擴(kuò)增

    利用細(xì)菌通用引物擴(kuò)增1 6 S r R N A V 6區(qū)片段,上游引物為以下兩引物的等量混合:9 6 7 F(5’-CAACGCGAAGAACCTTACC-3’)和967F(5’-ATACGCGAGGAACCTTACC-3’),下游引物為1046R(5’-CGACARCCATGCASCACCT-3’)。PCR反應(yīng)體系(50 μL):5 μL緩沖液(10×)(無Mg2+)、3 μL 25 mmol/L MgSO4、1 μL 10 mmol/L的dNTP、上下游引物(10 μmol/L)各0.5 μL,50 ng腸道細(xì)菌基因組模板DNA,最后加入0.4 μL Pfu DNA 聚合酶(2.5 U/μL),加無菌水至50 μL。為避免非特異性擴(kuò)增采用降落PCR的方法進(jìn)行擴(kuò)增,反應(yīng)程序:95 ℃預(yù)變性3 min;95 ℃變性1 min;退火溫度從65~55 ℃,每個(gè)循環(huán)降0.5 ℃,55 ℃ 15 個(gè)循環(huán),退火30 s;72 ℃延伸1 min,一共35 個(gè)循環(huán);最后72 ℃延伸7 min。PCR產(chǎn)物用3%的瓊脂糖凝膠電泳檢測并用MiniBEST Agarose Gel DNA Extraction Kit純化。純化后的擴(kuò)增產(chǎn)物用Qubit?dsDNA HS Assay Kit進(jìn)行PCR產(chǎn)物的定量。

    1.3.3.3 Ion torrent測序

    先用Ion Plus Fragment Library Kit和Ion Xpress?Barcode Adapters Kit進(jìn)行測序文庫的構(gòu)建。然后進(jìn)行測序模板的制備,Ion OneTouch? ES儀器自動對攜帶有測序模板的Ion Sphere? Particles進(jìn)行富集。最后在Ion torrent個(gè)人化操作基因組測序平臺上進(jìn)行上機(jī)測序。

    1.3.3.4 測序數(shù)據(jù)生物信息學(xué)和多變量統(tǒng)計(jì)學(xué)分析

    原始數(shù)據(jù)利用Fast QC軟件進(jìn)行質(zhì)控,并利用NG Stoolkits過濾掉低質(zhì)量序列,利用Usearch方法對所得的高質(zhì)量序列,根據(jù)樣品條碼進(jìn)行樣品的分選,在97%的相似水平下劃分分類操作單元(operational taxonomic units,OTUs),生成以序列數(shù)代表的每個(gè)樣本中每個(gè)OTUs豐度的OTU Table(biom文件)。然后按照微生物生態(tài)學(xué)定量分析(quantitative insights into microbial ecology,QIIME)[21]的流程進(jìn)行物種分類地位的確定,并進(jìn)行稀疏曲線的繪制,計(jì)算觀測到的物種(observed species)指數(shù)和Shannon指數(shù)。最后,計(jì)算結(jié)果用R語言

    3.1.3 軟件進(jìn)行作圖?;谏鲜龅玫降腛TU Table應(yīng)用基于線性判別分析(linear discriminant analysis,LDA)的LEfSe(linear discriminant analysis effect size)算法篩選組間差異OTUs。根據(jù)這些差異OTUs的豐度利用R語言3.1.3軟件進(jìn)行豐度聚類熱圖的繪制和斯皮爾曼相關(guān)性系數(shù)的計(jì)算,并利用Cytoscape 3.2.1軟件繪制互作網(wǎng)絡(luò)圖。

    1.4 數(shù)據(jù)統(tǒng)計(jì)分析

    采用了R語言軟件的wilcox檢驗(yàn)對數(shù)據(jù)進(jìn)行了顯著性分析,P<0.05表示差異顯著。

    2 結(jié)果與分析

    2.1 L. helveticus TS206干預(yù)后小鼠腸道微生物整體特征

    2.1.1 各樣品中序列和OTUs數(shù)目概況

    22 個(gè)樣本中除去2 個(gè)誤差較大的樣本,余下的20 個(gè)樣本共獲得3 001 460 條高質(zhì)量序列,每個(gè)樣品的序列分布情況如表1所示。全部序列的長度均在80~120 bp的范圍內(nèi),由于細(xì)菌16S rRNA基因的V6區(qū)平均長度為100 bp,所以本次實(shí)驗(yàn)的測序長度基本能覆蓋V6區(qū)。

    表1 各樣品中序列和OTUs(97%相似水平)分布Table 1 Numbers of reads and OTUs (97% similarity level) in the library of each sample

    在本實(shí)驗(yàn)中,根據(jù)16S rRNA基因序列相似度與細(xì)菌分類地位之間的對應(yīng)關(guān)系,在97%的序列相似性水平上經(jīng)過Usearch的方法劃分OTUs。每個(gè)樣本的OTUs分布如表1所示,20 個(gè)樣本一共獲得了1 617 個(gè)OTUs。

    2.1.2 L. helveticus TS206干預(yù)后小鼠腸道微生物整體結(jié)構(gòu)的改變

    通過97%相似水平劃分得到1 617 個(gè)OTUs,OTUs的代表序列與Greengene 1308軟件的數(shù)據(jù)庫進(jìn)行比對注釋得到每個(gè)OTU的物種注釋情況,最后統(tǒng)計(jì)各個(gè)樣本中各分類水平相對豐度,結(jié)果如圖1所示。

    圖1 L. helveticus TS206對小鼠腸道菌群門水平上相對豐度的影響Fig. 1 Effect of L. helveticus TS206 on relative abundance of dominant bacterial phyla in the intestine of mice

    圖1顯示所有OTUs被劃分為8 個(gè)門,其中硬壁菌門(Firmicutes)和擬桿菌門(Bacteroidetes)為各組小鼠糞便中的優(yōu)勢菌門,占總序列數(shù)的97.49%。其中硬壁菌門的含量最豐富,約占總序列數(shù)的70.66%,含有1 147 個(gè)OTUs,其中99.3%的序列屬于Clostridia綱,含有1 135 個(gè)OTUs。而擬桿菌門則占總序列數(shù)的26.83%,含有313 個(gè)OTUs全部屬于Bacteroidales目,其中54.9%的序列是S24-7科細(xì)菌,分屬177 個(gè)OTUs;其次是擬桿菌科(Bacteroidaceae)占17.8%,分屬37 個(gè)OTUs;理研菌科(Rikenellaceae)含有36 個(gè)OTUs,占11%;其他序列分屬于普雷沃氏菌科(Prevotellaceae)、Odoribacteraceae科、Paraprevotellaceae科和紫單孢菌科(Porphyromonadaceae),含有63 個(gè)OTUs。其余的如放線菌門(Actinobacteria)3 個(gè)OTUs、藍(lán)細(xì)菌門(Cyanobacteria)9 個(gè)OTUs、脫鐵桿菌門(Deferribacteres)1 個(gè)OTU、變形菌門(Proteobacteria)17 個(gè)OTUs、無壁菌門(Tenericutes)59 個(gè)OTUs和疣微菌門(Verrucomicrobia)3 個(gè)OTUs,它們的相對豐度較低,在各樣品中所占比例基本都不足1%,另外,還有0.5%左右的無法確定具體分類地位的細(xì)菌。

    從圖1可以看到,L. helveticus TS206灌胃組小鼠腸道擬桿菌門(Bacteroidetes)細(xì)菌的相對豐度隨著灌胃L. helveticus TS206時(shí)間延長有降低的趨勢,在灌胃第6周擬桿菌門細(xì)菌相對豐度降到最低,由最初的27.5%降到17.8%;對照組小鼠腸道擬桿菌門細(xì)菌豐度雖然在第2、3周(B2、B3)有所波動,但實(shí)驗(yàn)期間基本維持在同一水平。伴隨著擬桿菌門細(xì)菌相對豐度的降低,L. helveticus TS206灌胃組小鼠腸道硬壁菌門細(xì)菌相對豐度相應(yīng)升高,由起初的71.4%增加到第6周的80.7%。

    圖2 L. helveticus TS206對小鼠腸道菌群科水平上相對豐度的影響Fig. 2 Effect of L. helveticus TS206 on relative abundance of major bacterial families in the intestine of mice

    各組小鼠腸道菌群在科水平上主要物種相對豐度的變化如圖2所示。在科的水平上豐度占優(yōu)勢地位的主要為S24-7、毛螺菌科(Lachnospiraceae)、瘤胃球菌科(Ruminococcaceae)、擬桿菌科(Bacteroidaceae)和理研菌科(Rikenellaceae)的細(xì)菌。經(jīng)過干預(yù)后,對瘤胃球菌科細(xì)菌作用持續(xù)性較好,對動物的生理影響也比較持久。而毛螺菌科細(xì)菌在灌胃L. helveticus TS206后相對豐度波動較大,但是總體上都高于對照組。同時(shí),在灌胃L. helveticus TS206后有顯著變化的理研菌科從起初的2.9%降低到第7周(A7)的0.7%,也顯著低于對照組(4.1%);S24-7科細(xì)菌從最初的27.5%降低到第6周(A6)的11.3%,也低于對照組的平均水平(17.4%)。綜合在門水平上的比較結(jié)果,可以得知灌胃L. helveticus TS206后小鼠腸道菌群擬桿菌門細(xì)菌豐度的降低可能是由于理研菌科和S24-7科細(xì)菌的降低引起的。

    另外,對各組小鼠腸道中豐度較低的腸桿菌科(Enterobacteriaceae)細(xì)菌的相對豐度也進(jìn)行了比較,發(fā)現(xiàn)雖然它們在各組小鼠腸道中的總體相對豐度較低,且在各個(gè)樣本中的個(gè)體差異也比較大,但是總體而言在L. helveticus TS206組小鼠腸道中腸桿菌科細(xì)菌的相對豐度要比對照組低一個(gè)數(shù)量級。這一結(jié)果通過傳統(tǒng)的平板計(jì)數(shù)的方法得到了驗(yàn)證,結(jié)果如圖3所示。在整個(gè)實(shí)驗(yàn)過程中,L. helveticus TS206組小鼠腸道中腸菌科細(xì)菌數(shù)量具有降低趨勢,到灌胃第4周(A4)降到最低,相比于對照組,隨后一直處于較低水平;而對照組小鼠腸桿菌科細(xì)菌則是隨著實(shí)驗(yàn)時(shí)間的延長而增加。雖然在L. helveticus TS206灌胃前兩周時(shí),對照組小鼠腸桿菌數(shù)量都低于L. helveticus TS206組,差異極顯著(P<0.01),但是從第3周開始一直到實(shí)驗(yàn)結(jié)束,L. helveticus TS206灌胃組小鼠腸道內(nèi)腸桿菌數(shù)量一直顯著地低于對照組(P<0.05)。

    圖3 瑞士乳桿菌對小鼠腸道中腸桿菌科細(xì)菌數(shù)量的影響Fig. 3 Effect of L. helveticus TS206 on the number of intestinal Enterobacteriaceae in mice

    2.2 L. helveticus TS206對小鼠腸道微生物α多樣性的影響

    為了進(jìn)一步研究L. helveticus TS206對小鼠腸道微生物整體結(jié)構(gòu)的影響,從各樣本腸道菌群的豐富度和均勻度的層面研究小鼠腸道微生物的多樣性。以各樣品在97%相似水平上劃分的OTUs數(shù)據(jù)為對象,通過QIIME分析了各樣本的α多樣性,首先通過稀釋曲線和反映樣本物種豐富度的observed species指數(shù)來研究兩組樣本物種多樣性的差異,其次通過反應(yīng)微生物多樣性和均勻度的Shannon指數(shù)來研究兩組樣本的多樣性和均勻度的差異,發(fā)現(xiàn)L. helveticus TS206對小鼠腸道微生物α多樣性造成了一定程度的影響。

    圖4 對照組和益生菌組各樣本稀釋曲線Fig. 4 Rarefaction curves of OUT diversity for control and probiotic groups

    圖5 L. helveticus TS206對腸道微生物α多樣性指數(shù)的影響Fig. 5 Effect of L. helveticus TS206 on α diversity index of gut microbiota

    圖4結(jié)果表明各組樣本的OTUs數(shù)目隨著測序深度的增加基本達(dá)到飽和,說明當(dāng)前測序深度足以發(fā)現(xiàn)各樣本生境中的大部分物種。而2 組樣本的稀釋曲線交錯(cuò),說明2 組小鼠腸道物種多樣性并沒有大的差異;這與圖5a兩組小鼠腸道微生物observed species指數(shù)差異分析箱圖的結(jié)果一致,雖然差異分析箱圖顯示對照組要略高于L. helveticus TS206組,但差異不顯著(P>0.05)。

    圖5b中兩組小鼠腸道微生物Shannon指數(shù)的差異分析也出現(xiàn)同樣結(jié)果,在灌胃L. helveticus TS206后α多樣性指數(shù)與對照組相比較低,但是差異不顯著(P>0.05),以上說明L. helveticus TS206對小鼠腸道微生物整體多樣性沒有顯著影響(P>0.05)。

    2.3 與2 組小鼠腸道糞便菌群結(jié)構(gòu)差異顯著相關(guān)的腸道關(guān)鍵OTUs

    上述腸道微生物整體結(jié)構(gòu)的分析發(fā)現(xiàn),2 組小鼠腸道微生物具有一定程度的差異。采用LEfSe差異分析算法篩選出L. helveticus TS206組與對照組間具有顯著差異的OTUs,結(jié)果如圖6所示。

    圖6 與L. helveticus TS206灌胃相關(guān)的腸道關(guān)鍵OTUsFig. 6 Key OTUs of gut microbiota responding to L. helveticus TS206

    通過LEfSe分析尋找到53 個(gè)OTUs是造成L. helveticus TS206組和對照組小鼠腸道菌群結(jié)構(gòu)差異的關(guān)鍵物種。相對于對照組,27 個(gè)OTUs在L. helveticus TS206組中豐度較低,26 個(gè)OTUs在L. helveticus TS206組豐度富集(圖6a)。

    從圖6b中可以發(fā)現(xiàn),共有3 個(gè)OTUs(OTU1 488、OTU1 588、OTU1 598)屬于腸桿菌科和1 個(gè)屬于瘤胃球菌科的OTUs均在對照組中富集;同時(shí),在這53 個(gè)差異OTUs中,有15 個(gè)OTUs屬于梭菌目(Clostridiales),而7 個(gè)屬于梭菌目的OTUs以及部分屬于擬桿菌屬的OTUs在L. helveticus TS206組豐度更高。在53 個(gè)差異OTUs中一共有10 個(gè)OTUs屬于毛螺菌科,其中有5 個(gè)(OTU1 248、OTU801、OTU771、OTU1 355、OTU 1 295)在L. helveticus TS206組富集,另外5 個(gè)(OTU468、OTU462、OTU1 238、OTU952、OTU659)在對照組中富集,這些均說明即使屬于同一分類水平的物種在腸道中可能會發(fā)揮不同的功能。

    2.4 與L. helveticus TS206相關(guān)的腸道關(guān)鍵物種對各組樣本腸道微生物整體結(jié)構(gòu)差異的影響

    為了探究通過LEfSe分析尋找到53 個(gè)關(guān)鍵OTUs對2 組樣本的整體分布情況,基于53 個(gè)關(guān)鍵OTUs進(jìn)行主成分分析,結(jié)果如圖7所示。

    圖7 關(guān)鍵OTUs對2 組樣本菌群結(jié)構(gòu)分布的影響Fig. 7 Effect of key OTUs on gut microbiota structures

    通過主成分分析可以發(fā)現(xiàn),基于53 個(gè)關(guān)鍵OTUs的相對豐度進(jìn)行降維分析,L. helveticus TS206組所有樣本聚為一類,對照組所有樣本也聚為一類,并且2 組在圖7中差異顯著(P<0.05),這從另一層面上說明這53 個(gè)差異顯著OTUs可能是造成L. helveticus TS206組與對照組小鼠腸道菌群結(jié)構(gòu)差異的關(guān)鍵物種。

    2.5 關(guān)鍵OTUs在2 組樣本間互作網(wǎng)絡(luò)的特征

    為了進(jìn)一步研究53 個(gè)差異OTUs在2 組樣本間的互作關(guān)系的特征,根據(jù)53 個(gè)OTUs在2 組樣本中的豐度信息計(jì)算斯皮爾曼相關(guān)性系數(shù),根據(jù)各個(gè)OTUs間的相關(guān)系數(shù)大小繪制互作網(wǎng)絡(luò)圖,如圖8所示。在對照組富集的和在L. helveticus TS206組富集的OTUs之間互作關(guān)系緊密,并且在2 組富集的OTUs之間均呈負(fù)相關(guān)關(guān)系,這也充分說明這些差異OTUs在兩組樣本的腸道中可能存在競爭抑制的關(guān)系。并且,只在對照組富集的屬于腸桿菌科的OTUs與L. helveticus TS206組多個(gè)OTUs具有負(fù)相關(guān)關(guān)系,最為明顯的是與屬于梭菌目的OTU41有負(fù)相關(guān)關(guān)系,這也提示該梭菌目的OTUs可能是L. helveticus TS206組中關(guān)鍵的OTUs,并能抵制腸桿菌科物種的增殖作用。另外,只在對照組中富集的屬于瘤胃球菌科的OTUs也與L. helveticus TS206組多個(gè)OTUs呈現(xiàn)負(fù)相關(guān)關(guān)系;此外,在L. helveticus TS206組富集的豐度較高的擬桿菌屬的OTUs也與對照組多個(gè)OTUs呈現(xiàn)負(fù)相關(guān)關(guān)系。這些在互作網(wǎng)絡(luò)圖中呈現(xiàn)了豐富連接關(guān)系的OTUs節(jié)點(diǎn)可能是與L. helveticus TS206組和對照組小鼠腸道微生物結(jié)構(gòu)差異有關(guān)的關(guān)鍵OTUs,有待于更深入的探討。

    圖8 關(guān)鍵OTUs間互作關(guān)系網(wǎng)絡(luò)圖Fig. 8 Interaction network among key OTUs

    3 討 論

    人體腸道有種類繁多的微生物,而這些微生物與機(jī)體一直處于相互作用的動態(tài)平衡中,一旦這種平衡失調(diào)就會造成腸道疾病甚至是系統(tǒng)性的慢性疾病。在腸道微生物與宿主相互作用的過程中,宿主自身的基因型、年齡和免疫系統(tǒng)以及飲食結(jié)構(gòu)和營養(yǎng)的攝入等因素共同影響著腸道菌群的形成及多樣性組成[18-20]。食品的攝入一方面會改變腸道菌群的結(jié)構(gòu),另一方面食品中難以被人體吸收的成分通過腸道菌群的作用影響食品的營養(yǎng)。而益生菌作為當(dāng)前人們?nèi)粘I钪谐S玫氖称?,能夠有效地維持機(jī)體與腸道菌群處于相互作用的動態(tài)平衡中。

    本研究發(fā)現(xiàn)各組小鼠腸道菌群在門水平上的組成基本無差異,都以擬桿菌門和硬壁菌門為主,相對豐度超過95%,且硬壁菌門豐度最高,超過70%。有其他研究發(fā)現(xiàn),在正常人和小鼠腸道中擬桿菌門和硬壁菌門是兩大優(yōu)勢菌門,幾乎占腸道所有細(xì)菌的90%以上[21-22],并且硬壁菌門是最為優(yōu)勢的一類菌,大多數(shù)屬于梭菌綱[23],其次是豐度較低的放線菌門和變形菌門[24-25]。在本實(shí)驗(yàn)中,L. helveticus TS206組小鼠腸道菌群隨著灌胃時(shí)間的延長,擬桿菌門豐度有降低的趨勢,而硬壁菌門有升高的趨勢。

    屬于硬壁菌門的梭菌目的球形梭菌和柔嫩梭菌亞群中,大多數(shù)細(xì)菌是產(chǎn)丁酸鹽細(xì)菌[26],主要通過丁酸激酶或通過丁酰CoA、乙酰CoA轉(zhuǎn)移酶的作用產(chǎn)生丁酸[27]。丁酸是腸上皮細(xì)胞能量的主要來源,而且在調(diào)節(jié)炎癥性反應(yīng)、細(xì)胞增殖和凋亡及抗結(jié)直腸癌方面也有重要作用[28-29],例如Faecalibacterium prausnitzii就是典型的丁酸產(chǎn)生菌,具有明顯的抗炎作用[30]。本研究發(fā)現(xiàn)各組小鼠腸道中豐度最高的是梭菌目細(xì)菌,通過LEfSe分析發(fā)現(xiàn)梭菌目是L. helveticus TS206造成小鼠腸道菌群結(jié)構(gòu)差異的關(guān)鍵菌屬,且在L. helveticus TS206組小鼠腸道中豐度更高,這表明L. helveticus TS206可能會通過增加腸道中梭菌目中丁酸產(chǎn)生菌的數(shù)量來維持腸道菌群結(jié)構(gòu)的穩(wěn)定和腸道的健康。另外,也有研究顯示梭菌在機(jī)體腸道中可以通過分泌短鏈脂肪酸來誘導(dǎo)和調(diào)節(jié)調(diào)節(jié)性T細(xì)胞的增殖和分化來塑造腸道內(nèi)穩(wěn)定的免疫系統(tǒng)[31-32],在本實(shí)驗(yàn)中,L. helveticus TS206誘導(dǎo)的腸道梭菌目細(xì)菌的增加可能會對穩(wěn)定腸道免疫平衡有重要作用。但是由于測序片段較短、分類地位不夠明確,有待進(jìn)一步研究這些梭菌目細(xì)菌的具體分類地位。

    此外,本實(shí)驗(yàn)還發(fā)現(xiàn)在L. helveticus TS206灌胃后,豐度顯著增加的菌屬中大部分是對腸道健康有益的菌屬,而在對照組中一些對機(jī)體健康不利的細(xì)菌豐度較高。例如,由LEfSe分析得知L. helveticus TS206組中跟丁酸產(chǎn)生相關(guān)的梭菌目細(xì)菌具有更高的豐度,而瘤胃菌科只在對照組中富集。另一方面,在對照組中小鼠腸道中腸桿菌科細(xì)菌豐度顯著高于L. helveticus TS206組,這個(gè)結(jié)果也通過傳統(tǒng)的平板培養(yǎng)計(jì)數(shù)實(shí)驗(yàn)得到驗(yàn)證。腸桿菌科細(xì)菌是常見的內(nèi)毒素產(chǎn)生菌,能夠誘導(dǎo)機(jī)體產(chǎn)生炎癥性反應(yīng),是典型的條件致病菌[21]。王友湘[14]、Frece[15]等在先前的研究中也發(fā)現(xiàn)L. helveticus灌胃能夠顯著增加腸道內(nèi)乳酸菌和雙歧桿菌的數(shù)量,抑制腸桿菌的增殖,這與本研究得到的結(jié)論相似。

    在本研究中,還鑒定出與L. helveticus TS206灌胃顯著相關(guān)的其他菌屬。目前對這些菌屬的研究較少,而且對這些菌屬的具體作用也不明確。例如,分屬于毛螺菌科和瘤胃菌科的OTUs與2 組小鼠腸道菌群差異都顯著相關(guān),而且在L. helveticus TS206組小鼠腸道中豐度顯著增加的毛螺菌科OTUs數(shù)更多。有研究顯示,瘤胃菌科細(xì)菌在結(jié)直腸癌高風(fēng)險(xiǎn)人群腸道內(nèi)要高于低風(fēng)險(xiǎn)人群[33]。另一項(xiàng)研究也描述了毛螺菌科細(xì)菌與結(jié)直腸癌的相關(guān)性,研究者在評估結(jié)直腸癌與腸道菌群的相關(guān)性是發(fā)現(xiàn)腺瘤患者的一些OTUs相對豐度較高,包括與瘤胃菌科、假單胞菌屬和紫單胞菌科相關(guān)的OTUs,而與擬桿菌屬、毛螺菌科、梭菌目及其他梭菌屬相關(guān)的OTUs則相對豐度較低[34]。但是,瘤胃菌科和毛螺菌科細(xì)菌都是宿主腸道中常見菌屬,因此在本研究中2 組小鼠腸道中都有與之相關(guān)的OTUs的相對豐度增高,它們在腸道中的具體功能還應(yīng)進(jìn)一步通過功能分析來驗(yàn)證。

    4 結(jié) 論

    本研究通過給健康小鼠長期灌胃L. helveticus TS206,利用Ion torrent PGM測序技術(shù)和多變量統(tǒng)計(jì)學(xué)方法分析得知L. helveticus TS206組與對照組小鼠腸道菌群結(jié)構(gòu)具有顯著差異。根據(jù)測序的數(shù)據(jù)分析,得到了以下幾點(diǎn)提示:從門的水平上分析,各組小鼠腸道糞便中的主要菌群主要由硬壁菌門和擬桿菌門組成,同時(shí)L. helveticus TS206能夠一定程度上增加腸道微生物中硬壁菌門的相對豐度,并伴隨著擬桿菌門相對豐度的降低;從科水平上分析,豐度占優(yōu)勢地位的細(xì)菌主要為S24-7、毛螺菌科、瘤胃球菌科、理研菌科和擬桿菌科;此外,結(jié)合門水平分析結(jié)果,L. helveticus TS206組擬桿菌門細(xì)菌豐度的降低是由理研菌科和S24-7科豐度降低引起的;同時(shí),與對照組相比,L. helveticus TS206組腸桿菌科細(xì)菌數(shù)量具有降低趨勢。另外,經(jīng)LEfSe分析篩選出L. helveticus TS206組與對照組間具有顯著差異的53 個(gè)關(guān)鍵OTUs,其中腸桿菌科和瘤胃菌科在對照組中富集,梭菌目在L. helveticus TS206組中具有更高的豐度;經(jīng)主成分分析證實(shí)這53 個(gè)關(guān)鍵OTUs可能是造成L. helveticus TS206組與對照組小鼠腸道菌群結(jié)構(gòu)差異的關(guān)鍵物種;且通過互作網(wǎng)絡(luò)圖可知,2 組富集的OTUs之間均呈負(fù)相關(guān)關(guān)系,可見這些差異OTUs在2 組樣本的腸道中可能存在競爭抑制的關(guān)系。因此,研究經(jīng)生物信息學(xué)分析,初步顯示L. helveticus TS206可以通過抑制腸道有害微生物的生長、促進(jìn)有益菌的增殖來維持腸道菌群結(jié)構(gòu)的穩(wěn)態(tài)。然而,研究尚需進(jìn)一步的定量PCR等驗(yàn)證工作,以充分證實(shí)L. helveticus TS206對腸道菌群的調(diào)控能力。

    [1] LEDERBERG J. Infectious history[J]. Science, 2000, 288: 287-293.DOI:10.1126/science.288.5464.287.

    [2] 李旻. 人體腸道菌群結(jié)構(gòu)與宿主代謝的相關(guān)性研究[D]. 上海: 上海交通大學(xué), 2009: 15-18.

    [3] GOODMAN A L, MCNULTY N P, ZHAO Y, et al. Identifying genetic determinants needed to establish a human gut symbiontin its habitat[J]. Cell Host and Microbe, 2009, 6(3): 279-289. DOI:10.1016/j.chom.2009.08.003.

    [4] LEE Y K, MAZMANIAN S K. Has the microbiota played a critical role in the evolution of the adaptive immune system?[J]. Science,2010, 330: 1768-1773. DOI:10.1126/science.1195568.

    [5] ZHU Y, MICHELLE T, JOBIN C, et al. Gut microbiota and probiotics in colon tumorigenesis[J]. Cancer Letters, 2011, 309(2): 119-127.DOI:10.1128/mSphere.00001-15.

    [6] TURNBAUGH P J, HAMADY M, YATSUNENKO T, et al. A core gut microbiome in obese and lean twins[J]. Nature, 2008, 457: 480-484. DOI:10.1038/nature07540.

    [7] B?CKHED F, MANCHESTER J K, SEMENKOVICH C F, et al.Mechanisms underlying the resistance to diet-induced obesity in germfree mice[J]. Proceedings of the National Academy of Sciences, 2007,104(3): 979-984. DOI:10.1073/pnas.0605374104.

    [8] WEN L, LEY R E, VOLCHKOV P Y, et al. Innate immunity and intestinal microbiota in the development of type 1 diabetes[J]. Nature,2008, 455: 1109-1113. DOI:10.1038/nature07336.

    [9] QIN J J, LI Y R, CAI Z M, et al. A metagenome-wide association study of gut microbiota in type 2 diabetes[J]. Nature, 2012, 490: 55-60.

    [10] VéTIZOU M, PITT J M, DAILLèRE R, et al. Anticancer immunotherapy by CTLA-4 blockade relies on the gut microbiota[J].Science, 2015, 350: 1079-1084. DOI:10.1007/s00109-016-1401-8.

    [11] SALMINEN S J, GUEIMONDE M, ISOLAURI E. Probiotics that modify disease risk[J]. Nutrition, 2005, 135(5): 1294-1298.

    [12] 張和平. 我國益生乳酸菌研究與產(chǎn)業(yè)化現(xiàn)狀及發(fā)展對策[J]. 生物產(chǎn)業(yè)技術(shù), 2009(6): 53-55. DOI:10.3969/j.issn.1674-0319.2009.06.005.

    [13] 葉雷, 閆亞麗, 陳慶森, 等. 高通量測序技術(shù)在腸道微生物宏基因組學(xué)研究中的應(yīng)用[J]. 中國食品學(xué)報(bào), 2016, 16(7): 216-223.DOI:10.16429/j.1009-7848.2016.07.029.

    [14] 王友湘, 陳慶森. 瑞士乳桿菌對小鼠腸道微生物區(qū)系的影響[J]. 食品科學(xué), 2008, 29(9): 542-546. DOI:10.3321/j.issn:1002-6630.2008.09.129.

    [15] FRECE J, KOS B, SVETEC I K, et al. Synbiotic effect of Lactobacillus helveticus M92 and prebiotics on the intestinal microflora and immune system of mice[J]. Dairy Research, 2009,76(1): 98-104. DOI:10.1017/S0022029908003737.

    [16] TAVERNITI V, GUGLIELMETTI S. Health-promoting properties of Lactobacillus helveticus[J]. Frontiers in Microbiology, 2012, 3: 392.DOI:10.3389/fmicb.2012.00392.

    [17] SINGH N, ARIOLI S, WANG A, et al. Impact of Bifidobacterium bifidum MIMBb75 on mouse intestinal microorganisms[J]. FEMS Microbiology Ecology, 2013, 85(2): 369-375. DOI:10.1111/1574-6941.12124.

    [18] LEY R E, LOZUPONE C A, HAMADY M, et al. Worlds within worlds: evolution of the vertebrate gut microbiota[J]. Nature Reviews Microbiology, 2008, 6(10): 776-788. DOI:10.1371/journal.pgen.1001314.

    [19] ZHANG C H, ZHANG M H, PANG X Y, et al. Structural resilience of the gut microbiota in adult mice under high-fat dietary perturbations[J].ISME Journal, 2012, 6(10): 1848-1857.

    [20] TURNBAUGH P J, BACKHED F, FULTON L, et al. Diet-induced obesity is linked to marked but reversible alterations in the mouse distal gut microbiome[J]. Cell Host and Microbe, 2008, 3(4): 213-223.DOI:10.1016/j.chom.2008.02.015.

    [21] ECKBURG P B, BIK E M, BERNSTEIN C N, et al. Diversity of the human intestinal microbial flora[J]. Science, 2005, 308: 1635-1638.DOI:10.1126/science.1110591.

    [22] TURNBAUGH P J, LEY R E, MAHOWALD M A, et al. An obesityassociated gut microbiome with increased capacity for energy harvest[J]. Nature, 2006, 444: 1027-1031. DOI:10.1038/nature05414.

    [23] HAYASHI H, SAKAMOTO M, BENNO Y. Phylogenetic analysis of the human gut microbiota using 16S rDNA clone libraries and strictly anaerobic culture-based methods[J]. Microbiology and Immunology,2002, 46(8): 535-548.

    [24] FEI N, ZHAO L P. An opportunistic pathogen isolated from the gut of an obese human causes obesity in germfree mice[J]. ISME Journal,2013, 7(4): 880-884. DOI:10.1038/ismej.2012.153.

    [25] SHKOPOROV A N, KHOKHLOVA E V, KULAGINA E V, et al.Application of several molecular techniques to study numerically predominant Bifidobacterium spp. and Bacteroidales order strains in the feces of healthy children[J]. Bioscience, Biotechnology and Biochemistry, 2008, 72(3): 742-748. DOI:10.1271/bbb.70628.

    [26] PRYDE S E, DUNCAN S H, HOLD G L, et al. The microbiology of butyrate formation in the human colon[J]. FEMS Microbiology Letters, 2002, 217(2): 133-139. DOI:10.1016/j.cam.2009.02.103.

    [27] FUKUDA S, TOH H, HASE K, et al. Bif i dobacteria can protect from enteropathogenic infection through production of acetate[J]. Nature,2011, 469: 543-547. DOI:10.1038/nature09646.

    [28] VANHOUTVIN S A, TROOST F J, HAMER H M, et al. Butyrateinduced transcriptional changes in human colonic mucosa[J]. PLoS ONE, 2009, 4(8): 6759-6766. DOI:10.1371/journal.pone.0006759.

    [29] FUNG K Y, COSGROVE L, LOCKETT T, et al. A review of the potential mechanisms for the lowering of colorectal oncogenesis by butyrate[J]. British Journal of Nutrition, 2012, 108(5): 820-831.DOI:10.1017/S0007114512001948.

    [30] SOKOL H, PIGNEUR B, WATTERLOT L, et al. Faecalibacterium prausnitzii is an anti-inf l ammatory commensal bacterium identif i ed by gut microbiota analysis of Crohn disease patients[J]. Proceedings of the National Academy of Sciences, 2008, 105(43): 16731-16736.

    [31] ATARASHI K, TANOUE T, SHIMA T, et al. Induction of colonic regulatory T cells by indigenous Clostridium species[J]. Science, 2011,331: 337-341. DOI:10.1053/j.gastro.2011.07.013.

    [32] FURUSAWA Y, OBATA Y, FUKUDA S, et al. Commensal microbederived butyrate induces the differentiation of colonic regulatory T cells[J]. Nature, 2013, 506: 446-450. DOI:10.1038/nature12721.

    [33] MOORE W E, MOORE L H. Intestinal floras of populations that have a high risk of colon cancer[J]. Applied and Environmental Microbiology, 1995, 61(9): 3202-3207.

    [34] ZACKULAR J P, ROGERS M A, RUFFIN M T, et al. The human gut microbiome as a screening tool for colorectal cancer[J]. Cancer Prevention Research, 2014, 7(11): 1112-1121. DOI:10.1158/1940-6207.CAPR-14-0129.

    Modulation of Probiotic Lactobacillus helveticus on Gut Microbiota in Mice

    ZANG Kaili1, JIA Yan1, CUI Wenjing2, MA Xinying2, WANG Yong1, ZHAO Linsen2, ZHAO Pei1, YE Lei1, YAN Yali1,*, CHEN Qingsen1,*
    (1. Tianjin Key Laboratory of Food Biotechnology, College of Biotechnology and Food Science, Tianjin University of Commerce,Tianjin 300134, China; 2. Hebei Inatural Biological Technical Company, Shijiazhuang 050899, China)

    Objective: To explore the role of the probiotic strain Lactobacillus helveticus TS206 in regulating the intestinal microf l ora in mice and to unravel the mechanisms by which probiotics can maintain intestinal fl ora homeostasis. Methods:Male BALB/c mice were divided into probiotic, control and blank groups. The probiotic and control groups were respectively gavaged with L. helveticus suspension and an equal volume of normal saline once a day for 7 consecutive weeks, respectively. The blank group did not receive gavage. Stool samples were collected every week for extracting genomic DNA by means of a commercial kit. The V6 variable region of the 16S rRNA gene was subjected to highthroughput sequencing by using the Ion torrent personal genome machine system and the obtained data were analyzed using bioinformatic tools and multivariable statistical analysis. Results: All sequence reads were delineated into 1 617 operational taxonomic units (OTUs) at a 97% similarity level, belonging to 8 major bacterial phyla, with Firmicutes and Bacteroidetes being the dominant ones in all groups, accounting for 97.49% of the total number of sequences, amongwhich Firmicutes showed the highest abundance of more than 70%. The strain S24-7, Lachnospiraceae, Ruminococcaceae,Rikenellaceae, and Bacteroidaceae were the dominant bacterial families, and the Enterobacteriaceae population in the probiotic group was lower than in the control group and fell during the administration period. Totally 53 key OTUs were signif i cantly associated with the bacterial community structure in the intestine of the two groups of mice as determined by linear discriminant analysis effect size (LEfSe) analysis, of which, 27 OTUs were enriched in the control group, belonging to Enterobacteriaceae, Ruminococcaceae, Lachnospiraceae and Clostridium, and 26 OTUs in the L. helveticus TS206 group,belonging to Clostridium, whose abundance level was higher than in the control group, Lachnospiraceae and Bacteroides.Principal component analysis showed that the fecal microbial communities of the probiotic group were completely separated from those of the control group. Overall there were signif i cant differences in the structure of intestinal fl ora between the two groups, and the OTUs in the groups showed a negative correlation. Conclusion: L. helveticus TS206 can change the structure of intestinal fl ora, inhibiting the growth of harmful microbes, promoting the proliferation of benef i cial bacteria proliferation and consequently maintaining the structure of the intestinal fl ora at a steady state.

    Lactobacillus helveticus TS206; gut microbiota; Ion torrent personal genome machine sequencing technology;bioinformatics; multivariate statistical analysis

    10.7506/spkx1002-6630-201801024

    TS201.3

    A

    1002-6630(2018)01-0156-09

    臧凱麗, 賈彥, 崔文靜, 等. 瑞士乳桿菌調(diào)控小鼠腸道菌群變化規(guī)律的研究[J]. 食品科學(xué), 2018, 39(1): 156-164.

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

    ZANG Kaili, JIA Yan, CUI Wenjing, et al. Modulation of probiotic Lactobacillus helveticus on gut microbiota in mice[J]. Food Science, 2018, 39(1): 156-164. (in Chinese with English abstract)

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

    2017-04-10

    國家自然科學(xué)基金面上項(xiàng)目(31071522)

    臧凱麗(1993—),女,碩士研究生,研究方向?yàn)樯锘钚晕镔|(zhì)與腸道健康。E-mail:1623249305@qq.com

    *通信作者簡介:閻亞麗(1962—)女,副教授,碩士,研究方向?yàn)榘l(fā)酵生物技術(shù)、微生物學(xué)。E-mail:yyali@tjcu.edu.cn

    陳慶森(1957—),男,教授,碩士,研究方向?yàn)榘l(fā)酵生物技術(shù)、食源性生物活性物質(zhì)與腸道健康。E-mail:chqsen@tjcu.edu.cn

    猜你喜歡
    梭菌灌胃菌群
    “云雀”還是“貓頭鷹”可能取決于腸道菌群
    中老年保健(2022年2期)2022-08-24 03:20:50
    丁酸梭菌的篩選、鑒定及生物學(xué)功能分析
    中國飼料(2021年17期)2021-11-02 08:15:10
    “水土不服”和腸道菌群
    科學(xué)(2020年4期)2020-11-26 08:27:06
    小鼠、大鼠灌胃注意事項(xiàng)
    豬肉毒梭菌中毒癥的發(fā)生和診療方法
    來曲唑灌胃后EM大鼠病灶體積及COX-2 mRNA、survivin蛋白表達(dá)變化
    肉牛剩余采食量與瘤胃微生物菌群關(guān)系
    西藏牦牛肉毒梭菌中毒病的防治
    西藏科技(2015年10期)2015-09-26 12:10:24
    咽部菌群在呼吸道感染治療中的臨床應(yīng)用
    實(shí)驗(yàn)小鼠的灌胃給藥技巧
    亚洲天堂av无毛| 欧美日韩视频精品一区| 国产不卡av网站在线观看| 欧美激情久久久久久爽电影 | 亚洲国产欧美一区二区综合| h视频一区二区三区| 老司机亚洲免费影院| 亚洲avbb在线观看| 一边摸一边抽搐一进一小说 | 日韩成人在线观看一区二区三区| 天天躁日日躁夜夜躁夜夜| 一进一出抽搐动态| 女性生殖器流出的白浆| 国产亚洲精品第一综合不卡| 我要看黄色一级片免费的| 久久久国产一区二区| 99riav亚洲国产免费| 欧美黑人欧美精品刺激| 丝袜在线中文字幕| 中文字幕av电影在线播放| bbb黄色大片| 丝袜美腿诱惑在线| 亚洲国产看品久久| 麻豆乱淫一区二区| 男男h啪啪无遮挡| 黄色a级毛片大全视频| 精品久久久久久久毛片微露脸| 高潮久久久久久久久久久不卡| 国产一区有黄有色的免费视频| 久久人妻av系列| 日韩一卡2卡3卡4卡2021年| 久久久国产成人免费| 一边摸一边抽搐一进一小说 | 国产精品久久久人人做人人爽| 老司机靠b影院| 99国产精品一区二区三区| 国产单亲对白刺激| 日本av手机在线免费观看| 欧美乱妇无乱码| 亚洲av电影在线进入| 欧美在线黄色| 12—13女人毛片做爰片一| 国产精品熟女久久久久浪| 在线播放国产精品三级| 欧美乱妇无乱码| 国产有黄有色有爽视频| 热99久久久久精品小说推荐| 午夜视频精品福利| 王馨瑶露胸无遮挡在线观看| 一级片'在线观看视频| a级毛片在线看网站| a级片在线免费高清观看视频| 日韩中文字幕欧美一区二区| 制服人妻中文乱码| 麻豆国产av国片精品| 精品久久久久久久毛片微露脸| 欧美精品高潮呻吟av久久| 女同久久另类99精品国产91| 淫妇啪啪啪对白视频| 国产欧美日韩综合在线一区二区| 久久中文字幕一级| 久久久精品免费免费高清| 美女视频免费永久观看网站| 搡老熟女国产l中国老女人| 丰满人妻熟妇乱又伦精品不卡| 精品熟女少妇八av免费久了| 满18在线观看网站| 国产一区二区三区综合在线观看| videos熟女内射| 亚洲专区国产一区二区| 国产成人av激情在线播放| 亚洲成国产人片在线观看| 久久亚洲真实| 日本黄色日本黄色录像| 欧美av亚洲av综合av国产av| 国产精品99久久99久久久不卡| 成人影院久久| 亚洲一卡2卡3卡4卡5卡精品中文| 俄罗斯特黄特色一大片| 一进一出好大好爽视频| 成人国产一区最新在线观看| 在线 av 中文字幕| 国产精品电影一区二区三区 | 香蕉国产在线看| 最近最新中文字幕大全免费视频| 亚洲熟妇熟女久久| 久久精品国产亚洲av香蕉五月 | av免费在线观看网站| 免费在线观看日本一区| 亚洲第一欧美日韩一区二区三区 | 日本wwww免费看| 久久天躁狠狠躁夜夜2o2o| 亚洲av日韩精品久久久久久密| av天堂在线播放| 国产精品久久久久久精品电影小说| 精品一区二区三区av网在线观看 | 亚洲色图综合在线观看| aaaaa片日本免费| 亚洲国产欧美一区二区综合| 日韩熟女老妇一区二区性免费视频| 国产主播在线观看一区二区| 国产欧美日韩一区二区三区在线| 麻豆成人av在线观看| netflix在线观看网站| 91大片在线观看| 一级片'在线观看视频| 亚洲国产av影院在线观看| 桃红色精品国产亚洲av| 国产男女超爽视频在线观看| 久久毛片免费看一区二区三区| 国产一区二区三区视频了| netflix在线观看网站| 久久精品亚洲av国产电影网| 99精国产麻豆久久婷婷| 极品少妇高潮喷水抽搐| 欧美 亚洲 国产 日韩一| 亚洲人成电影免费在线| a在线观看视频网站| 久久亚洲真实| 天堂中文最新版在线下载| 久久久国产精品麻豆| 免费av中文字幕在线| 久久久精品94久久精品| 亚洲成国产人片在线观看| 欧美日韩一级在线毛片| av超薄肉色丝袜交足视频| 色94色欧美一区二区| a级片在线免费高清观看视频| 免费少妇av软件| 欧美 日韩 精品 国产| 欧美 亚洲 国产 日韩一| 怎么达到女性高潮| 国产亚洲午夜精品一区二区久久| 69精品国产乱码久久久| 香蕉丝袜av| 亚洲熟妇熟女久久| 1024香蕉在线观看| 免费少妇av软件| 亚洲精品在线观看二区| 999精品在线视频| 亚洲午夜精品一区,二区,三区| 在线观看舔阴道视频| 久久99热这里只频精品6学生| 伦理电影免费视频| 国产精品免费大片| 久久久久视频综合| 性少妇av在线| 亚洲国产av影院在线观看| 亚洲一卡2卡3卡4卡5卡精品中文| 久久精品亚洲av国产电影网| 午夜福利在线观看吧| 久久人妻福利社区极品人妻图片| 高清视频免费观看一区二区| 日本欧美视频一区| 蜜桃国产av成人99| 手机成人av网站| av电影中文网址| 亚洲中文日韩欧美视频| 国产精品电影一区二区三区 | 国产精品美女特级片免费视频播放器 | 午夜福利免费观看在线| 国产不卡av网站在线观看| 丁香欧美五月| 人人妻,人人澡人人爽秒播| 久久婷婷成人综合色麻豆| 国产视频一区二区在线看| 一级黄色大片毛片| 老司机在亚洲福利影院| 激情在线观看视频在线高清 | 国产一区二区 视频在线| 熟女少妇亚洲综合色aaa.| 女人精品久久久久毛片| 日韩有码中文字幕| 制服诱惑二区| 1024视频免费在线观看| 男女免费视频国产| 国产免费福利视频在线观看| 天堂中文最新版在线下载| 香蕉丝袜av| 狠狠狠狠99中文字幕| 免费女性裸体啪啪无遮挡网站| 老司机影院毛片| 国内毛片毛片毛片毛片毛片| 视频区欧美日本亚洲| 一二三四社区在线视频社区8| 日韩欧美一区二区三区在线观看 | 精品一区二区三区av网在线观看 | 一本综合久久免费| www.精华液| 人妻久久中文字幕网| 人妻久久中文字幕网| 成人三级做爰电影| 一本色道久久久久久精品综合| 日本一区二区免费在线视频| 亚洲一卡2卡3卡4卡5卡精品中文| 香蕉国产在线看| 波多野结衣一区麻豆| 如日韩欧美国产精品一区二区三区| 91国产中文字幕| 日本av免费视频播放| 最黄视频免费看| 一区二区av电影网| 极品人妻少妇av视频| 国产国语露脸激情在线看| 国产极品粉嫩免费观看在线| 亚洲一卡2卡3卡4卡5卡精品中文| 成年人黄色毛片网站| 成人手机av| 久久性视频一级片| 制服诱惑二区| 欧美黄色淫秽网站| 精品欧美一区二区三区在线| 自拍欧美九色日韩亚洲蝌蚪91| 黑人猛操日本美女一级片| 正在播放国产对白刺激| 丰满迷人的少妇在线观看| 亚洲精品国产色婷婷电影| 欧美大码av| 丰满迷人的少妇在线观看| 搡老乐熟女国产| 国产人伦9x9x在线观看| 丰满迷人的少妇在线观看| 日本av手机在线免费观看| 黄色毛片三级朝国网站| 美国免费a级毛片| 午夜激情av网站| 亚洲精品中文字幕在线视频| 97在线人人人人妻| 午夜福利一区二区在线看| 国产精品免费一区二区三区在线 | 久久精品亚洲精品国产色婷小说| 免费日韩欧美在线观看| 在线观看免费视频日本深夜| 亚洲成人免费电影在线观看| 飞空精品影院首页| 亚洲伊人色综图| 女人高潮潮喷娇喘18禁视频| 国产在线一区二区三区精| 极品少妇高潮喷水抽搐| 久久国产精品男人的天堂亚洲| 少妇被粗大的猛进出69影院| 一二三四社区在线视频社区8| 成在线人永久免费视频| 69av精品久久久久久 | 成人永久免费在线观看视频 | 久久久久精品人妻al黑| 国产男女超爽视频在线观看| 国产福利在线免费观看视频| 亚洲av国产av综合av卡| 国产极品粉嫩免费观看在线| 精品欧美一区二区三区在线| 精品一区二区三区av网在线观看 | 又紧又爽又黄一区二区| 国产成人欧美| 男女午夜视频在线观看| 国产日韩欧美在线精品| 久久精品国产亚洲av香蕉五月 | 男女无遮挡免费网站观看| 国产精品国产av在线观看| av有码第一页| 久久久久久久精品吃奶| 香蕉国产在线看| 极品人妻少妇av视频| 人妻久久中文字幕网| 搡老乐熟女国产| 99精品在免费线老司机午夜| 久久精品国产综合久久久| 五月天丁香电影| 黑人巨大精品欧美一区二区mp4| 亚洲 国产 在线| 久热这里只有精品99| 在线天堂中文资源库| 免费观看av网站的网址| 久久天躁狠狠躁夜夜2o2o| 一区在线观看完整版| 天天躁夜夜躁狠狠躁躁| 国产一卡二卡三卡精品| 亚洲人成电影观看| 十分钟在线观看高清视频www| 国产精品久久电影中文字幕 | 欧美激情高清一区二区三区| 日韩人妻精品一区2区三区| 国产av又大| 久久久久久亚洲精品国产蜜桃av| 国产精品久久电影中文字幕 | 国产淫语在线视频| 91麻豆精品激情在线观看国产 | 国产福利在线免费观看视频| 亚洲精华国产精华精| 免费观看人在逋| 国产av又大| 国产精品免费大片| 国产99久久九九免费精品| 岛国在线观看网站| 建设人人有责人人尽责人人享有的| 乱人伦中国视频| 午夜福利乱码中文字幕| 十八禁网站网址无遮挡| 肉色欧美久久久久久久蜜桃| 国产精品成人在线| 国产单亲对白刺激| 久久99热这里只频精品6学生| 亚洲一区中文字幕在线| 亚洲国产欧美日韩在线播放| 国产精品麻豆人妻色哟哟久久| 日韩欧美三级三区| 精品国产一区二区三区久久久樱花| 一边摸一边抽搐一进一小说 | 男女边摸边吃奶| 精品一区二区三区四区五区乱码| 亚洲欧洲精品一区二区精品久久久| 中文字幕高清在线视频| 欧美日韩精品网址| 老熟妇仑乱视频hdxx| 精品一品国产午夜福利视频| 国产免费现黄频在线看| 超色免费av| 色综合欧美亚洲国产小说| 日日夜夜操网爽| 在线观看66精品国产| 国产成人系列免费观看| 美女高潮到喷水免费观看| 国产精品 国内视频| 亚洲国产av新网站| 最近最新中文字幕大全电影3 | 在线 av 中文字幕| 久久久久久久大尺度免费视频| 成年动漫av网址| 国产男女超爽视频在线观看| 国产av又大| 一级毛片电影观看| 成年人午夜在线观看视频| 国产男靠女视频免费网站| 捣出白浆h1v1| 久久久久久久久久久久大奶| 老司机午夜十八禁免费视频| 久久久久精品国产欧美久久久| 王馨瑶露胸无遮挡在线观看| 在线观看www视频免费| 日韩视频在线欧美| 国产一区二区三区在线臀色熟女 | 在线观看免费高清a一片| 亚洲欧美日韩另类电影网站| av天堂在线播放| 亚洲精品美女久久久久99蜜臀| 日韩人妻精品一区2区三区| 日韩大片免费观看网站| 男男h啪啪无遮挡| 婷婷成人精品国产| 午夜福利一区二区在线看| 香蕉丝袜av| 日韩有码中文字幕| 成人影院久久| av不卡在线播放| 三上悠亚av全集在线观看| 午夜福利,免费看| 男女边摸边吃奶| 国产精品免费大片| 国产精品熟女久久久久浪| 人人妻,人人澡人人爽秒播| 91精品三级在线观看| 亚洲欧美一区二区三区黑人| 国产成人一区二区三区免费视频网站| 久久国产精品人妻蜜桃| 亚洲午夜精品一区,二区,三区| 日韩成人在线观看一区二区三区| 女人久久www免费人成看片| 国产精品成人在线| 国产精品欧美亚洲77777| 欧美乱码精品一区二区三区| 午夜精品国产一区二区电影| 黄片小视频在线播放| 日韩大码丰满熟妇| www.熟女人妻精品国产| 97人妻天天添夜夜摸| 在线观看一区二区三区激情| 天天添夜夜摸| bbb黄色大片| 看免费av毛片| 天堂俺去俺来也www色官网| 热99久久久久精品小说推荐| 日韩人妻精品一区2区三区| 2018国产大陆天天弄谢| 热99国产精品久久久久久7| 丝袜喷水一区| 亚洲九九香蕉| 热99re8久久精品国产| 国产av国产精品国产| 在线观看www视频免费| 18禁国产床啪视频网站| 久久精品成人免费网站| svipshipincom国产片| 黄色视频在线播放观看不卡| 天天躁狠狠躁夜夜躁狠狠躁| 久久人妻福利社区极品人妻图片| av片东京热男人的天堂| 中文字幕高清在线视频| 日韩免费av在线播放| 日韩三级视频一区二区三区| 亚洲天堂av无毛| 精品一区二区三区av网在线观看 | 国产亚洲精品第一综合不卡| 露出奶头的视频| 视频区图区小说| 日日夜夜操网爽| 欧美+亚洲+日韩+国产| 亚洲人成77777在线视频| 日本撒尿小便嘘嘘汇集6| 欧美日韩国产mv在线观看视频| 精品亚洲成国产av| 91成年电影在线观看| 国产精品美女特级片免费视频播放器 | 最新美女视频免费是黄的| 黄网站色视频无遮挡免费观看| 人妻一区二区av| 亚洲色图 男人天堂 中文字幕| 高潮久久久久久久久久久不卡| 91精品三级在线观看| 免费看十八禁软件| 中文字幕最新亚洲高清| 亚洲精品在线美女| 国产欧美亚洲国产| 国产一区二区在线观看av| 亚洲专区字幕在线| 丝袜喷水一区| 69av精品久久久久久 | 韩国精品一区二区三区| 午夜福利免费观看在线| av片东京热男人的天堂| 热99国产精品久久久久久7| h视频一区二区三区| 汤姆久久久久久久影院中文字幕| 9色porny在线观看| 亚洲国产成人一精品久久久| 国产一区二区激情短视频| 国产精品秋霞免费鲁丝片| 日韩中文字幕视频在线看片| 99国产精品99久久久久| 十八禁高潮呻吟视频| 日本a在线网址| 久久ye,这里只有精品| 妹子高潮喷水视频| 19禁男女啪啪无遮挡网站| 精品一区二区三区四区五区乱码| 午夜福利影视在线免费观看| 国产不卡一卡二| 男女边摸边吃奶| 精品国产乱子伦一区二区三区| 18禁观看日本| 欧美激情高清一区二区三区| 成人三级做爰电影| 国产亚洲av高清不卡| 19禁男女啪啪无遮挡网站| e午夜精品久久久久久久| 下体分泌物呈黄色| 久久av网站| 午夜精品久久久久久毛片777| 中文字幕高清在线视频| 欧美精品高潮呻吟av久久| a级毛片在线看网站| 久久99热这里只频精品6学生| 亚洲精品在线观看二区| 国产福利在线免费观看视频| 搡老乐熟女国产| 蜜桃国产av成人99| 国产精品自产拍在线观看55亚洲 | 久久午夜亚洲精品久久| 王馨瑶露胸无遮挡在线观看| 国产淫语在线视频| 国产黄频视频在线观看| 99国产精品99久久久久| tube8黄色片| 久久精品亚洲av国产电影网| 在线观看66精品国产| 99国产精品一区二区蜜桃av | 欧美成人免费av一区二区三区 | 十分钟在线观看高清视频www| 五月天丁香电影| 免费在线观看完整版高清| 亚洲人成伊人成综合网2020| 水蜜桃什么品种好| 日本黄色日本黄色录像| 亚洲精品乱久久久久久| 一级黄色大片毛片| 久久久国产精品麻豆| 国产男女内射视频| 老司机靠b影院| 亚洲国产精品一区二区三区在线| 国产日韩一区二区三区精品不卡| 999久久久国产精品视频| 久久精品国产亚洲av高清一级| 久久人人爽av亚洲精品天堂| 在线观看人妻少妇| 91国产中文字幕| 精品人妻1区二区| 一夜夜www| 正在播放国产对白刺激| 欧美+亚洲+日韩+国产| av超薄肉色丝袜交足视频| 国产成人免费无遮挡视频| 国产淫语在线视频| 久9热在线精品视频| www日本在线高清视频| 91精品三级在线观看| 国产精品亚洲av一区麻豆| 丰满饥渴人妻一区二区三| 精品国产一区二区三区久久久樱花| 搡老乐熟女国产| avwww免费| 国产高清videossex| 国产精品免费一区二区三区在线 | 欧美日韩福利视频一区二区| 免费在线观看影片大全网站| 在线观看免费日韩欧美大片| 黄网站色视频无遮挡免费观看| 国产黄频视频在线观看| 国产高清视频在线播放一区| 欧美日韩一级在线毛片| 久久久久久久大尺度免费视频| 免费人妻精品一区二区三区视频| 午夜激情av网站| 日韩大码丰满熟妇| 国产精品久久久久久精品古装| 一区二区三区精品91| 成人特级黄色片久久久久久久 | 国产一区二区 视频在线| 日韩欧美三级三区| 久久人妻av系列| 男女床上黄色一级片免费看| 91老司机精品| 十八禁网站网址无遮挡| 啦啦啦视频在线资源免费观看| 亚洲精品美女久久av网站| 中国美女看黄片| 久久久国产一区二区| 80岁老熟妇乱子伦牲交| 久久久精品94久久精品| 他把我摸到了高潮在线观看 | 日本五十路高清| aaaaa片日本免费| 精品人妻1区二区| 夜夜夜夜夜久久久久| 国产xxxxx性猛交| av片东京热男人的天堂| 欧美亚洲日本最大视频资源| 欧美人与性动交α欧美精品济南到| 建设人人有责人人尽责人人享有的| 免费人妻精品一区二区三区视频| 日本黄色视频三级网站网址 | 国产成人一区二区三区免费视频网站| 日韩免费av在线播放| 精品福利观看| 午夜福利在线观看吧| 国产精品亚洲一级av第二区| 桃红色精品国产亚洲av| 在线观看www视频免费| 制服诱惑二区| 不卡一级毛片| 国产熟女午夜一区二区三区| 久久久久久久国产电影| 久久久精品国产亚洲av高清涩受| 搡老岳熟女国产| 人人妻人人添人人爽欧美一区卜| 免费女性裸体啪啪无遮挡网站| 日本黄色视频三级网站网址 | 两性午夜刺激爽爽歪歪视频在线观看 | 久久久精品免费免费高清| 久久国产亚洲av麻豆专区| 国产精品久久久久久人妻精品电影 | 天天操日日干夜夜撸| 女人被躁到高潮嗷嗷叫费观| 欧美中文综合在线视频| 老司机影院毛片| 成人特级黄色片久久久久久久 | 在线亚洲精品国产二区图片欧美| 午夜成年电影在线免费观看| 久久精品人人爽人人爽视色| 日韩一区二区三区影片| 久久久久久亚洲精品国产蜜桃av| 99re6热这里在线精品视频| 不卡av一区二区三区| 中文字幕人妻熟女乱码| 考比视频在线观看| 久久精品国产综合久久久| 久久九九热精品免费| 国产日韩欧美在线精品| 免费一级毛片在线播放高清视频 | 女人被躁到高潮嗷嗷叫费观| 亚洲欧洲日产国产| 欧美激情极品国产一区二区三区| 亚洲专区国产一区二区| 黄色片一级片一级黄色片| 91麻豆精品激情在线观看国产 | 怎么达到女性高潮| 久久久久久久久久久久大奶| 国产男靠女视频免费网站| 黄色怎么调成土黄色| 久久久久久久久久久久大奶| 男女午夜视频在线观看| 一级毛片女人18水好多| 一级,二级,三级黄色视频| 91成人精品电影| 亚洲色图综合在线观看| 久久精品成人免费网站| 无遮挡黄片免费观看| 一级毛片精品| 欧美亚洲 丝袜 人妻 在线| 欧美在线黄色| 亚洲av片天天在线观看| 日本av免费视频播放| 人妻 亚洲 视频| 高清欧美精品videossex| 女警被强在线播放| 人妻一区二区av| 亚洲色图av天堂|