• <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)小鼠的灌胃給藥技巧
    欧美97在线视频| 成人毛片a级毛片在线播放| 国产av在哪里看| 国产精品一区二区三区四区久久| 大香蕉97超碰在线| 亚洲三级黄色毛片| 日韩人妻高清精品专区| 综合色av麻豆| 国产视频内射| 日本熟妇午夜| 国产不卡一卡二| 综合色丁香网| 国产一级毛片在线| 天天一区二区日本电影三级| 永久免费av网站大全| 免费不卡的大黄色大毛片视频在线观看 | 亚洲四区av| 日韩欧美国产在线观看| 欧美三级亚洲精品| 一级片'在线观看视频| 亚洲国产最新在线播放| av在线老鸭窝| 永久免费av网站大全| 久久久久久久久久黄片| 18禁在线无遮挡免费观看视频| 亚洲av电影不卡..在线观看| 亚洲av免费在线观看| 国产黄色免费在线视频| av专区在线播放| 91久久精品国产一区二区成人| av一本久久久久| 亚洲av福利一区| 在线观看美女被高潮喷水网站| 日韩三级伦理在线观看| 七月丁香在线播放| 亚洲av成人精品一二三区| av播播在线观看一区| 久久久午夜欧美精品| 日韩av不卡免费在线播放| 成人毛片60女人毛片免费| 身体一侧抽搐| 欧美精品一区二区大全| 国产精品一区二区性色av| 国产精品1区2区在线观看.| 九九在线视频观看精品| 又粗又硬又长又爽又黄的视频| 亚洲四区av| 国产高清三级在线| 欧美精品一区二区大全| av在线天堂中文字幕| 国产乱来视频区| 91精品一卡2卡3卡4卡| 亚洲av中文av极速乱| 国产成人福利小说| 国产精品一二三区在线看| 80岁老熟妇乱子伦牲交| 亚洲自拍偷在线| 久久综合国产亚洲精品| 欧美成人午夜免费资源| 亚洲欧美一区二区三区黑人 | 国产精品国产三级国产专区5o| 国产又色又爽无遮挡免| 亚洲成人精品中文字幕电影| 国产精品女同一区二区软件| 又爽又黄无遮挡网站| 搞女人的毛片| 最新中文字幕久久久久| 卡戴珊不雅视频在线播放| 亚洲av电影不卡..在线观看| 日本免费在线观看一区| 久久久成人免费电影| 亚洲国产av新网站| 亚洲人成网站在线播| 女的被弄到高潮叫床怎么办| 国产探花在线观看一区二区| 一个人看视频在线观看www免费| 身体一侧抽搐| 最近最新中文字幕免费大全7| 欧美性猛交╳xxx乱大交人| 成人午夜高清在线视频| 18禁在线无遮挡免费观看视频| 最近的中文字幕免费完整| 伊人久久国产一区二区| 欧美日本视频| 国产免费福利视频在线观看| 亚洲欧美日韩无卡精品| 精品久久久久久久久av| 日日摸夜夜添夜夜爱| 亚洲欧美清纯卡通| 成人av在线播放网站| av国产久精品久网站免费入址| 好男人视频免费观看在线| 国产片特级美女逼逼视频| 国产黄色视频一区二区在线观看| 免费黄色在线免费观看| 亚洲四区av| 最近2019中文字幕mv第一页| 久99久视频精品免费| 国产69精品久久久久777片| 美女高潮的动态| 中文字幕亚洲精品专区| 蜜桃亚洲精品一区二区三区| av在线亚洲专区| 哪个播放器可以免费观看大片| 日韩av在线大香蕉| 国产成人a区在线观看| 纵有疾风起免费观看全集完整版 | 久久久久精品性色| 97超碰精品成人国产| 国产一区二区三区av在线| 久久6这里有精品| 国产在线一区二区三区精| 亚洲av二区三区四区| 美女黄网站色视频| 韩国高清视频一区二区三区| 天天一区二区日本电影三级| 亚洲国产精品专区欧美| 午夜福利视频精品| 国产日韩欧美在线精品| 纵有疾风起免费观看全集完整版 | 国产淫片久久久久久久久| 久久这里有精品视频免费| 日日啪夜夜撸| 日本免费a在线| 成人亚洲精品av一区二区| 三级男女做爰猛烈吃奶摸视频| 一级毛片 在线播放| 国产片特级美女逼逼视频| 看黄色毛片网站| 久久国产乱子免费精品| 2021少妇久久久久久久久久久| 天堂网av新在线| 久久久亚洲精品成人影院| 听说在线观看完整版免费高清| 欧美 日韩 精品 国产| av福利片在线观看| 国产亚洲最大av| 91精品国产九色| 国产男女超爽视频在线观看| 一本—道久久a久久精品蜜桃钙片 精品乱码久久久久久99久播 | 国内揄拍国产精品人妻在线| 国产精品嫩草影院av在线观看| 国产午夜精品一二区理论片| 热99在线观看视频| 国产麻豆成人av免费视频| 欧美+日韩+精品| 草草在线视频免费看| 2021少妇久久久久久久久久久| 欧美变态另类bdsm刘玥| 亚洲欧美日韩无卡精品| 久久久色成人| 欧美人与善性xxx| 国产伦精品一区二区三区视频9| 精品久久久精品久久久| 精品久久久久久成人av| 18禁动态无遮挡网站| 天堂√8在线中文| 好男人视频免费观看在线| 久久精品国产亚洲网站| 可以在线观看毛片的网站| 亚洲av日韩在线播放| 男的添女的下面高潮视频| 男女边摸边吃奶| 欧美日韩精品成人综合77777| 少妇丰满av| 免费看美女性在线毛片视频| 简卡轻食公司| 亚洲不卡免费看| 久久精品国产亚洲av涩爱| 久久久午夜欧美精品| 菩萨蛮人人尽说江南好唐韦庄| 国国产精品蜜臀av免费| 一级毛片久久久久久久久女| 色综合站精品国产| av卡一久久| 亚洲高清免费不卡视频| 三级国产精品片| 国产成人精品福利久久| 亚洲最大成人av| 国产一级毛片在线| 人妻夜夜爽99麻豆av| 97人妻精品一区二区三区麻豆| 午夜福利视频精品| 80岁老熟妇乱子伦牲交| 美女黄网站色视频| 国产一级毛片在线| 男女啪啪激烈高潮av片| 春色校园在线视频观看| 亚洲欧美成人精品一区二区| 精品一区二区免费观看| 国产精品伦人一区二区| 天天躁夜夜躁狠狠久久av| 嫩草影院新地址| 国产精品麻豆人妻色哟哟久久 | 亚洲成人久久爱视频| 欧美激情久久久久久爽电影| 国产成人精品婷婷| 精品99又大又爽又粗少妇毛片| 两个人的视频大全免费| av国产久精品久网站免费入址| 成人午夜精彩视频在线观看| 免费电影在线观看免费观看| 特大巨黑吊av在线直播| 亚洲av中文av极速乱| 久久午夜福利片| 亚洲人成网站在线观看播放| 久久这里只有精品中国| 国产单亲对白刺激| 三级国产精品欧美在线观看| av国产久精品久网站免费入址| 中文字幕av在线有码专区| 草草在线视频免费看| 97精品久久久久久久久久精品| 亚洲av二区三区四区| 淫秽高清视频在线观看| 天堂网av新在线| 美女大奶头视频| 老司机影院成人| 看免费成人av毛片| 综合色av麻豆| 青春草视频在线免费观看| 国产一区二区在线观看日韩| 在现免费观看毛片| 成人欧美大片| 亚洲精品国产av成人精品| 欧美3d第一页| 久久久久久久久久久免费av| 国产成年人精品一区二区| 国产片特级美女逼逼视频| 日韩欧美一区视频在线观看 | xxx大片免费视频| 国产乱来视频区| 免费看光身美女| 高清日韩中文字幕在线| av女优亚洲男人天堂| 搞女人的毛片| 亚洲一级一片aⅴ在线观看| 亚洲成人久久爱视频| 国产精品99久久久久久久久| av免费在线看不卡| 男插女下体视频免费在线播放| 大又大粗又爽又黄少妇毛片口| 午夜日本视频在线| 一个人观看的视频www高清免费观看| 成人亚洲精品av一区二区| 欧美三级亚洲精品| 麻豆国产97在线/欧美| 国产乱来视频区| h日本视频在线播放| 国产男人的电影天堂91| 亚洲精品成人av观看孕妇| 亚洲av中文字字幕乱码综合| 2022亚洲国产成人精品| 搡老乐熟女国产| 日韩一区二区三区影片| 青春草视频在线免费观看| 男的添女的下面高潮视频| 久久精品综合一区二区三区| 午夜激情福利司机影院| 97超碰精品成人国产| 日本与韩国留学比较| 日本免费a在线| 啦啦啦啦在线视频资源| 中国美白少妇内射xxxbb| 黄片无遮挡物在线观看| 国产高清有码在线观看视频| 国产精品久久久久久久久免| 黄色一级大片看看| 久久99热这里只有精品18| 国产69精品久久久久777片| 精品一区二区三区人妻视频| 亚洲va在线va天堂va国产| 亚洲最大成人中文| 在线观看美女被高潮喷水网站| 午夜精品一区二区三区免费看| 日韩精品有码人妻一区| 欧美 日韩 精品 国产| 日韩精品青青久久久久久| 久久精品综合一区二区三区| 18禁动态无遮挡网站| 亚洲av中文av极速乱| 亚洲四区av| 久久久久久久久久久免费av| 久久97久久精品| 成人欧美大片| 国产av不卡久久| 尤物成人国产欧美一区二区三区| 少妇被粗大猛烈的视频| 国产精品伦人一区二区| 久久99蜜桃精品久久| 国产片特级美女逼逼视频| 成人高潮视频无遮挡免费网站| 精品一区在线观看国产| 精品国内亚洲2022精品成人| 亚洲国产欧美在线一区| 丝袜喷水一区| 在线免费观看的www视频| 欧美成人午夜免费资源| 亚洲欧美成人综合另类久久久| 啦啦啦啦在线视频资源| 亚洲精品,欧美精品| 久99久视频精品免费| 肉色欧美久久久久久久蜜桃 | 人妻系列 视频| 大话2 男鬼变身卡| 亚洲四区av| 大陆偷拍与自拍| 成人午夜高清在线视频| 日韩,欧美,国产一区二区三区| 精品一区在线观看国产| 免费人成在线观看视频色| 色吧在线观看| 十八禁国产超污无遮挡网站| 丰满人妻一区二区三区视频av| 国产久久久一区二区三区| 日本色播在线视频| 天堂俺去俺来也www色官网 | 亚洲欧美中文字幕日韩二区| 中文天堂在线官网| 能在线免费观看的黄片| 中文字幕av成人在线电影| 久久久久网色| 永久网站在线| 成年人午夜在线观看视频 | 综合色丁香网| 日韩在线高清观看一区二区三区| 五月伊人婷婷丁香| 国产精品熟女久久久久浪| 极品少妇高潮喷水抽搐| 久久精品国产亚洲av天美| 久久热精品热| 亚洲欧美清纯卡通| 亚洲人成网站高清观看| 日本-黄色视频高清免费观看| 欧美性猛交╳xxx乱大交人| 蜜桃久久精品国产亚洲av| 午夜福利高清视频| 午夜福利在线观看吧| 久久久欧美国产精品| 亚洲美女视频黄频| 麻豆成人午夜福利视频| 看黄色毛片网站| 九九久久精品国产亚洲av麻豆| 国产老妇伦熟女老妇高清| 成人高潮视频无遮挡免费网站| 国国产精品蜜臀av免费| 熟妇人妻不卡中文字幕| 国产免费一级a男人的天堂| 99九九线精品视频在线观看视频| av黄色大香蕉| 日韩伦理黄色片| 一级av片app| 中文欧美无线码| 日本黄色片子视频| 亚洲精品日韩av片在线观看| 欧美极品一区二区三区四区| 日本wwww免费看| 高清毛片免费看| 热99在线观看视频| 午夜久久久久精精品| 99久久精品热视频| 国产综合懂色| 亚洲精品色激情综合| 免费黄频网站在线观看国产| 亚洲av福利一区| 欧美一级a爱片免费观看看| 黑人高潮一二区| 亚洲在线观看片| 中文精品一卡2卡3卡4更新| 国产成人a区在线观看| 国产成年人精品一区二区| 国产精品麻豆人妻色哟哟久久 | 别揉我奶头 嗯啊视频| 国产精品av视频在线免费观看| 亚洲美女搞黄在线观看| 日韩制服骚丝袜av| 十八禁国产超污无遮挡网站| 一边亲一边摸免费视频| 国产伦精品一区二区三区四那| 久久久久久国产a免费观看| 久久久精品免费免费高清| 成人特级av手机在线观看| 日韩欧美一区视频在线观看 | 九草在线视频观看| 一区二区三区免费毛片| 欧美xxⅹ黑人| 在现免费观看毛片| 亚洲激情五月婷婷啪啪| 女人被狂操c到高潮| 亚洲成人精品中文字幕电影| 久久久久久久国产电影| 亚洲av.av天堂| 十八禁国产超污无遮挡网站| 欧美精品国产亚洲| 国产成人a区在线观看| 亚洲av免费在线观看| 夜夜爽夜夜爽视频| 综合色av麻豆| 伊人久久精品亚洲午夜| 大片免费播放器 马上看| 人体艺术视频欧美日本| av女优亚洲男人天堂| 日韩av不卡免费在线播放| 人妻系列 视频| 国产一区二区三区av在线| 晚上一个人看的免费电影| 精品国产三级普通话版| 免费黄频网站在线观看国产| 国产黄a三级三级三级人| 成人欧美大片| 免费大片18禁| 网址你懂的国产日韩在线| 精品久久久精品久久久| 1000部很黄的大片| 麻豆精品久久久久久蜜桃| 亚洲国产精品sss在线观看| 午夜精品国产一区二区电影 | 一级毛片电影观看| 特级一级黄色大片| 亚洲婷婷狠狠爱综合网| 日本黄大片高清| 国产美女午夜福利| 精品人妻视频免费看| 如何舔出高潮| 一级毛片 在线播放| 精品一区二区三卡| 亚洲av电影在线观看一区二区三区 | 亚洲最大成人中文| 能在线免费看毛片的网站| 免费观看在线日韩| 午夜激情欧美在线| eeuss影院久久| 尤物成人国产欧美一区二区三区| 久久久国产一区二区| 中文字幕av在线有码专区| 综合色丁香网| 老司机影院毛片| 美女cb高潮喷水在线观看| 高清午夜精品一区二区三区| 亚洲人与动物交配视频| 哪个播放器可以免费观看大片| 国产久久久一区二区三区| 亚洲av免费在线观看| 日韩欧美一区视频在线观看 | 亚洲精品中文字幕在线视频 | 亚洲怡红院男人天堂| 久久久久久久久久久丰满| 亚洲综合色惰| 午夜精品一区二区三区免费看| 又粗又硬又长又爽又黄的视频| 日韩视频在线欧美| 色综合站精品国产| 国产精品人妻久久久久久| 欧美性感艳星| 久久久久国产网址| 午夜福利在线观看吧| 久久精品久久精品一区二区三区| 岛国毛片在线播放| 性插视频无遮挡在线免费观看| 一个人看的www免费观看视频| 好男人视频免费观看在线| 国产av码专区亚洲av| 2021少妇久久久久久久久久久| 三级国产精品欧美在线观看| 三级毛片av免费| 国产永久视频网站| 天堂俺去俺来也www色官网 | 18禁裸乳无遮挡免费网站照片| 人妻系列 视频| 天天躁夜夜躁狠狠久久av| 国产麻豆成人av免费视频| 日产精品乱码卡一卡2卡三| 成人鲁丝片一二三区免费| 少妇人妻一区二区三区视频| 九九爱精品视频在线观看| 美女国产视频在线观看| 亚洲三级黄色毛片| 搡老乐熟女国产| 精品少妇黑人巨大在线播放| 久久韩国三级中文字幕| 欧美三级亚洲精品| 亚洲欧美一区二区三区国产| 中文欧美无线码| 欧美日韩亚洲高清精品| 七月丁香在线播放| 欧美变态另类bdsm刘玥| 久久韩国三级中文字幕| 男的添女的下面高潮视频| 国产一级毛片在线| 国产一区亚洲一区在线观看| 91狼人影院| 97热精品久久久久久| 在线观看av片永久免费下载| 亚洲精品国产av蜜桃| 亚洲国产欧美在线一区| 内射极品少妇av片p| 亚洲va在线va天堂va国产| 22中文网久久字幕| 欧美xxxx性猛交bbbb| 精品人妻偷拍中文字幕| 91精品国产九色| 97超视频在线观看视频| 精品一区二区三卡| 久久久久久国产a免费观看| 欧美高清成人免费视频www| 欧美日韩视频高清一区二区三区二| 精品不卡国产一区二区三区| 免费在线观看成人毛片| 亚洲综合色惰| 日韩在线高清观看一区二区三区| av国产久精品久网站免费入址| 精品人妻偷拍中文字幕| 午夜福利成人在线免费观看| 亚洲av成人av| 人妻制服诱惑在线中文字幕| 国产老妇伦熟女老妇高清| 国产av在哪里看| av黄色大香蕉| av一本久久久久| 亚洲av二区三区四区| 水蜜桃什么品种好| 亚洲av成人精品一二三区| 日韩一区二区视频免费看| 国精品久久久久久国模美| 国产免费一级a男人的天堂| 久热久热在线精品观看| 又粗又硬又长又爽又黄的视频| 十八禁网站网址无遮挡 | 午夜精品国产一区二区电影 | 免费观看性生交大片5| 国产成人精品久久久久久| 超碰97精品在线观看| 亚洲精品乱久久久久久| 小蜜桃在线观看免费完整版高清| 最近视频中文字幕2019在线8| 亚洲,欧美,日韩| 精品久久久久久久久亚洲| 日日干狠狠操夜夜爽| 中文资源天堂在线| 超碰av人人做人人爽久久| 伦精品一区二区三区| 亚洲欧美日韩卡通动漫| 大陆偷拍与自拍| 99久久精品国产国产毛片| 欧美激情久久久久久爽电影| 国产男女超爽视频在线观看| 国产精品久久久久久av不卡| 在线a可以看的网站| 国模一区二区三区四区视频| 大香蕉久久网| 97热精品久久久久久| 三级国产精品欧美在线观看| 能在线免费观看的黄片| 国产av在哪里看| 国产精品伦人一区二区| 久久这里只有精品中国| 天堂网av新在线| 欧美另类一区| 国产精品99久久久久久久久| 欧美zozozo另类| 免费看不卡的av| 在线a可以看的网站| 91精品伊人久久大香线蕉| 在线a可以看的网站| 最近手机中文字幕大全| 午夜福利网站1000一区二区三区| 有码 亚洲区| 亚洲av免费高清在线观看| 日本免费在线观看一区| 精品久久久精品久久久| 国产精品久久久久久久电影| 毛片女人毛片| 午夜福利高清视频| 亚洲一区高清亚洲精品| 男女下面进入的视频免费午夜| 尤物成人国产欧美一区二区三区| 亚州av有码| a级一级毛片免费在线观看| 国产精品精品国产色婷婷| 人妻制服诱惑在线中文字幕| 久久午夜福利片| 大话2 男鬼变身卡| 一级黄片播放器| 97人妻精品一区二区三区麻豆| 亚洲不卡免费看| 亚洲18禁久久av| 国产成人午夜福利电影在线观看| 国产黄色小视频在线观看| av天堂中文字幕网| 中文乱码字字幕精品一区二区三区 | 亚洲精品,欧美精品| 午夜精品国产一区二区电影 | 国产不卡一卡二| 美女被艹到高潮喷水动态| 一区二区三区免费毛片| 亚洲欧美日韩无卡精品| 久久这里有精品视频免费| 性色avwww在线观看| 亚洲成人中文字幕在线播放| 国产免费又黄又爽又色| 一个人观看的视频www高清免费观看| 国产极品天堂在线| 国产色婷婷99| av卡一久久| 亚洲自拍偷在线| 日韩欧美一区视频在线观看 | 国产精品1区2区在线观看.| 日本黄色片子视频| 2018国产大陆天天弄谢| 色网站视频免费| 欧美 日韩 精品 国产| 免费不卡的大黄色大毛片视频在线观看 | 插阴视频在线观看视频| 女人久久www免费人成看片| 成人无遮挡网站|