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

    麻杏石甘湯對流感病毒感染小鼠腸道菌群及趨化因子CCL5、CXCL10的影響

    2021-01-06 00:07:14盧芳國張香港陳純靜田維毅
    中草藥 2021年1期
    關(guān)鍵詞:奧司趨化因子菌門

    王 平,趙 澄,盧芳國*,吳 濤,張香港,陳純靜,肖 榮,寧 毅,魏 科,李 玲,蔡 琨,田維毅

    麻杏石甘湯對流感病毒感染小鼠腸道菌群及趨化因子CCL5、CXCL10的影響

    王 平1, 2,趙 澄1,盧芳國1*,吳 濤1,張香港1,陳純靜1,肖 榮1,寧 毅1,魏 科1,李 玲1,蔡 琨2,田維毅2

    1. 湖南中醫(yī)藥大學(xué),湖南 長沙 410208 2. 貴州中醫(yī)藥大學(xué),貴州 貴陽 550025

    探索麻杏石甘湯(Maxing Shigan Decoction)通過影響小鼠腸道菌群組成及趨化因子產(chǎn)生而防治流感病毒感染的潛在機(jī)制。通過鼻腔接種法建立A型流感病毒感染小鼠模型,經(jīng)ig給藥3 d和7 d后,觀察小鼠體質(zhì)量等一般情況;HE染色檢測結(jié)腸組織病理變化;免疫組化檢測結(jié)腸組織趨化因子CCL5、CXCL10的表達(dá);酶聯(lián)免疫法檢測結(jié)腸組織中CCL5、CXCL10含量,并取小鼠腸內(nèi)容物進(jìn)行16S rRNA基因V3~V4可變區(qū)測序、物種注釋及聚類,進(jìn)行Alpha多樣性和Beta多樣性分析,線性判別式分析 [linear discriminant analysis(LDA)effect size,LEfSe] 篩選組間差異物種,冗余分析及Spearman等級相關(guān)系數(shù)分析腸道菌群與CCL5、CXCL10的關(guān)聯(lián)性。給藥3 d后,模型組小鼠體質(zhì)量下降(<0.01),腸黏膜固有層中炎性細(xì)胞明顯浸潤且高表達(dá)CCL5和CXCL10(<0.01),結(jié)腸組織勻漿中CCL5和CXCL10含量增加(<0.01);與模型組比較,各藥物組小鼠體質(zhì)量明顯增加,炎性細(xì)胞浸潤減少,CCL5、CXCL10陽性表達(dá)細(xì)胞減少,腸組織勻漿中CCL5、CXCL10含量降低,以奧司他韋組和麻杏石甘湯組差異顯著(<0.05、0.01)。菌群多樣性分析和物種差異分析結(jié)果顯示,給藥3 d后,模型組小鼠腸道菌群組成與對照組及各藥物組存在明顯差異。模型組小鼠變形菌門、埃希菌屬相對豐度較對照組明顯增加(<0.05),厚壁菌門、乳桿菌屬、糞球菌屬相對豐度明顯降低(<0.05、0.01);與模型組比較,各藥物組變形菌門、埃希菌屬相對豐度降低(<0.05),奧司他韋組和麻杏石甘湯組厚壁菌門相對豐度顯著增加(<0.01),抗病毒顆粒組顫螺旋菌屬相對豐度增加(<0.05),麻杏石甘湯組乳桿菌屬、糞球菌屬相對豐度增加(<0.05、0.01)。給藥7 d后,模型組小鼠結(jié)腸組織病變減輕,趨化因子表達(dá)下降,變形菌門、埃希菌屬相對豐度較給藥3 d后明顯降低(<0.01);與模型組比較,各藥物組疣微菌門相對豐度明顯降低(<0.05、0.01),擬桿菌門相對豐度明顯升高(<0.05)。關(guān)聯(lián)性分析顯示,小鼠結(jié)腸組織中CCL5、CXCL10含量與埃希菌屬、克雷伯菌屬、梭菌屬、糞球菌屬等相對豐度顯著相關(guān)(<0.05、0.01)。A型流感病毒感染可引起小鼠腸道菌群結(jié)構(gòu)紊亂及免疫功能失衡,麻杏石甘湯通過調(diào)節(jié)腸道菌群結(jié)構(gòu)并影響趨化因子的產(chǎn)生,對流感病毒引起的腸道免疫損傷有一定的保護(hù)作用。

    A型流感病毒;麻杏石甘湯;腸道菌群;趨化因子;CCL5;CXCL10

    流行性感冒(以下簡稱流感)是由流感病毒(influenza virus)引起的一種急性呼吸道傳染病。A型流感病毒(influenza A virus,IAV)是引起人和動物感染的主要病原體。流感臨床表現(xiàn)多樣,常見的癥狀包括發(fā)熱、咳嗽、頭痛、乏力,并常伴有惡心、嘔吐、腹瀉、腹痛等消化道癥狀。有研究報道,流感患者胃腸道癥狀發(fā)生率為30.9%[1]。流感病毒感染致腸道菌群失調(diào)及黏膜局部免疫功能紊亂可能是流感容易出現(xiàn)消化道癥狀的重要原因[2]。

    腸道菌群從生命早期開始逐步有序地定植在消化道內(nèi),形成穩(wěn)定的腸道微生態(tài)系統(tǒng)。腸道菌群及其代謝產(chǎn)物與宿主的能量代謝、免疫穩(wěn)態(tài)以及健康狀況密切相關(guān)。越來越多的證據(jù)表明,腸道菌群種類和比例改變對流感病毒感染的結(jié)局可能產(chǎn)生積極影響或負(fù)面作用,這與腸道菌群及其代謝產(chǎn)物對宿主免疫功能的影響密不可分[3]。研究發(fā)現(xiàn),流感病毒感染后,腸道菌群結(jié)構(gòu)組成發(fā)生改變,并且腸上皮細(xì)胞可以釋放過量的趨化因子和促炎細(xì)胞因子,腸道菌群與黏膜局部炎癥因子之間相互作用,進(jìn)而影響流感的進(jìn)程[2,4]。Wang等[2]研究發(fā)現(xiàn),流感病毒感染可促進(jìn)肺源性CCR9+CD4+T細(xì)胞進(jìn)入腸道,分泌γ干擾素(interferon-γ,IFN-γ)干擾腸道微生物群穩(wěn)態(tài),引起腸道菌群失調(diào),而紊亂的腸道菌群又可促進(jìn)小腸上皮細(xì)胞分泌白細(xì)胞介素-15(interleukin-15,IL-15),誘導(dǎo)Th17細(xì)胞極化,促進(jìn)IL-17等細(xì)胞因子的產(chǎn)生,進(jìn)而加重腸組織損傷。Li等[4]研究發(fā)現(xiàn),流感病毒感染可促進(jìn)腸道中變形菌門腸桿菌科細(xì)菌生長,破壞黏膜屏障功能,誘導(dǎo)腸上皮細(xì)胞過度表達(dá)IL-22、IFN-α、IL-17A等促炎細(xì)胞因子,引起炎癥反應(yīng)和腸組織損傷。趨化因子CCL5(CC chemokine ligand 5)和CXCL10(C-X-C motif chemokine ligand 10)是引起流感免疫病理損傷的重要炎癥因子,能強(qiáng)烈趨化炎癥細(xì)胞向病灶部位聚集,與流感病毒感染引起的炎癥損傷程度呈正相關(guān)[5-6]。研究發(fā)現(xiàn)大腸埃希菌或肺炎克雷伯菌等可以促進(jìn)CCL5及CXCL10等趨化因子的表達(dá),從而加重炎性損傷程度[7-8]。因此,修復(fù)腸道微生態(tài)失衡,改善黏膜局部微環(huán)境對于防治流感有重要意義。源自漢代張仲景《傷寒論》的麻杏石甘湯(Maxing Shigan Decoction,MXSGD)是常用的防治流感的中藥經(jīng)方。本課題組前期研究發(fā)現(xiàn),該方具有干預(yù)病毒吸附、抑制病毒增殖、抑制趨化因子及炎癥介質(zhì)釋放等作用[9-11]。本研究以流感病毒感染BALB/c小鼠為模型,探討流感病毒對小鼠腸道菌群及腸黏膜局部趨化因子的影響,及MXSGD的干預(yù)作用,以期為腸道菌群作為流感防控新策略中的重要靶點(diǎn)提供依據(jù),并進(jìn)一步揭示MXSGD防治流感的可能機(jī)制,為拓展其臨床應(yīng)用提供理論支持和實(shí)驗(yàn)數(shù)據(jù)。

    1 材料

    1.1 動物

    6~8周BALB/c雄性小鼠60只,體質(zhì)量(18±2)g,購于湖南斯萊克景達(dá)實(shí)驗(yàn)動物有限公司 [許可證號SCXK(湘)2016-0002,動物批號43004700062561],于湖南中醫(yī)藥大學(xué)實(shí)驗(yàn)動物中心[許可證號SCXK(湘)2015-0003] 飼養(yǎng)。飼養(yǎng)條件:溫度(23±1)℃,濕度(65±5)%。所有實(shí)驗(yàn)符合湖南中醫(yī)藥大學(xué)動物實(shí)驗(yàn)倫理學(xué)規(guī)定。

    1.2 病毒株

    流感病毒小鼠肺適應(yīng)株(A型,IAV,A/PR/8/34),由湖南師范大學(xué)病毒研究室惠贈。10日齡雞胚尿囊腔接種培養(yǎng)傳代,血凝效價1∶640以上者用于實(shí)驗(yàn)。正式實(shí)驗(yàn)前根據(jù)文獻(xiàn)方法[12],測定小鼠病毒半數(shù)致死量(LD50)為1×10?3.77,每0.1 毫升病毒液用滅菌生理鹽水稀釋成50 LD50,置于冰袋中備用。

    1.3 藥物與試劑

    MXSGD(麻黃9 g、杏仁9 g、石膏18 g、炙甘草6 g):麻黃(批號1901029)、杏仁(批號2019032702)、石膏(碎,棉裹,批號1905250032)、炙甘草(批號1806004),購自湖南中醫(yī)藥大學(xué)附屬第一醫(yī)院門診藥房,經(jīng)湖南中醫(yī)藥大學(xué)第一附屬醫(yī)院藥學(xué)部戴冰教授鑒定為正品,符合《中國藥典》2020年版的要求;磷酸奧司他韋膠囊(規(guī)格75 mg/粒,批號M1050,意大利Roche S.p.A. 公司生產(chǎn),上海羅氏制藥有限公司分裝)購自中南大學(xué)湘雅二附院門診藥房,用蒸餾水充分溶解,制成質(zhì)量濃度為2.165 mg/mL的混懸液;抗病毒顆粒(規(guī)格9 g/袋,批號1906311,四川光大制藥有限公司)購自湖南中醫(yī)藥大學(xué)附屬第一醫(yī)院門診藥房,用蒸餾水充分溶解,制成質(zhì)量濃度為0.39 g/mL的溶液。以上實(shí)驗(yàn)用藥4 ℃避光保存?zhèn)溆?。兔抗小鼠RANTES抗體(批號ab189841)、CXCL10多克隆抗體(批號ab9938);驢抗兔IgG熒光二抗(批號ab175772)均購于Abcam公司;通用二步法試劑盒(批號PV-9000,北京中杉金橋生物技術(shù)有限公司);小鼠CXCL10/IP-10 ELISA試劑盒(批號F10933,上海西唐生物科技有限公司)、CCL5/ Rantes ELISA試劑盒(批號F11450,上海西唐生物科技有限公司)。

    1.4 主要儀器

    BSC-1300 IIA2生物安全柜(蘇州安泰空氣技術(shù)有限公司);SW-CJ-1FD超凈工作臺(蘇州安泰空氣技術(shù)有限公司);Axioscope 5光學(xué)掃描顯微鏡+Axiocam 503 color 顯微鏡攝像頭 [蔡司科技(蘇州)有限公司];BBY24M高通量組織細(xì)胞破碎儀(美國Next Advance公司);Spark多功能酶標(biāo)儀(瑞士Tecan公司);TP-200D電子分析天平(湘儀天平儀器設(shè)備有限公司);RMZ135精密輪轉(zhuǎn)切片機(jī)(德國LEICA公司);TGL 20M高速冷凍離心機(jī)(長沙湘智離心機(jī)儀器有限公司);IlluminaMiseq測序儀(美國Illumina公司)。

    2 方法

    2.1 動物分組

    實(shí)驗(yàn)設(shè)對照組、模型組、奧司他韋組(化學(xué)藥對照組)、抗病毒顆粒組(中藥對照組)和MXSGD組,每組12只小鼠。

    2.2 MXSGD水煎液的制備

    參照文獻(xiàn)方法[13-14],采用麻黃先煎方法制備MXSGD水煎液。按其組成比例分別稱量5劑藥的藥材量。先取麻黃45 g,加入藥材總量10倍體積的蒸餾水,武火煎煮,待沸騰后,調(diào)整為文火煎煮25 min,去沫,再加入石膏90 g、杏仁45 g、炙甘草30 g,繼續(xù)文火煎煮30 min,煎煮完畢后濾過;二煎加入7倍體積蒸餾水武火煮沸后再文火煎煮20 min,煎煮完畢后濾過,合并2次濾液,水浴濃縮至含生藥0.605 g/mL,4 ℃避光保存?zhèn)溆谩?/p>

    2.3 模型制備及給藥

    實(shí)驗(yàn)動物適應(yīng)性飼養(yǎng)2~3 d后稱質(zhì)量。除對照組外,各組小鼠按照本課題前期研究已建立的方法建立A型流感病毒感染模型[15]。用乙醚輕度麻醉小鼠后,每只小鼠鼻孔內(nèi)均勻滴入50 LD50流感病毒液0.05 mL,建立模型。對照組小鼠隔離飼養(yǎng)在同等條件下的房間,并按同樣方法鼻腔接種0.9%氯化鈉溶液0.05 mL。

    按臨床等效劑量(臨床等效劑量按動物每千克體質(zhì)量占人體表面積的比值計算[16])于感染后24 h開始ig給藥,每天1次,每次0.2 mL,連續(xù)給藥3 d和7 d。奧司他韋、抗病毒顆粒、MXSGD的給藥劑量分別為20.14 mg/(kg?d)、4.03 mL/(kg?d)、5.22 g/(kg?d)。對照組和模型組均ig等量生理鹽水。

    2.4 標(biāo)本采集

    分別于感染后第4天和第8天(給藥后第3天和第7天),禁水禁食8 h后,每組隨機(jī)選取6只小鼠,稱質(zhì)量后,摘眼球放血法處死小鼠,轉(zhuǎn)移至超凈工作臺內(nèi)解剖,用無菌鑷子將部分結(jié)腸段內(nèi)容物取下,置于EP管中,凍存于?80 ℃冰箱。收集一部分結(jié)腸組織置于4%多聚甲醛中固定;另一部分稱重后,液氮迅速冷凍保存?zhèn)溆谩?/p>

    2.5 結(jié)腸組織病理學(xué)觀察

    蘇木素-伊紅(HE)染色檢測結(jié)腸組織病理變化,用4%多聚甲醛固定結(jié)腸組織、石蠟包埋、切片(4~5 μm)、HE染色后,在光學(xué)顯微鏡下(×200)觀察病理變化。

    2.6 免疫組織化學(xué)檢測趨化因子CCL5、CXCL10在腸黏膜表達(dá)

    PV-9000通用二步法檢測腸組織中趨化因子CCL5、CXCL10蛋白表達(dá)水平。石蠟切片常規(guī)脫蠟,水化,抗原修復(fù),阻斷內(nèi)源性過氧化物酶,與CCL5、CXCL10一抗37 ℃孵育60 min,陰性對照組用PBS緩沖液代替一抗,PBS緩沖液沖洗后滴加反應(yīng)增強(qiáng)液,室溫孵育20 min,滴加增強(qiáng)酶標(biāo)山羊抗兔IgG聚合物,室溫孵育20 min,加DAB顯色劑處理5~8 min,蘇木素復(fù)染1 min,流水沖洗30 min,氨水返藍(lán),脫水、透明、封片,鏡下觀察結(jié)果。細(xì)胞膜和細(xì)胞質(zhì)棕褐色或者棕黃色著色為陽性細(xì)胞,以細(xì)胞不著色為陰性。用Axiocam 503 color采集圖像,選取結(jié)構(gòu)完整的1張切片,在高倍鏡下(×200),每張切片任選5個視野,用Image-proplus分析軟件進(jìn)行陽性顯色的平均積分光密度測定(陽性區(qū)域面積陽性強(qiáng)度),代表陽性細(xì)胞目的蛋白的表達(dá)量,然后求其平均值,作為該樣本的相對表達(dá)量。

    2.7 ELISA法檢測CCL5、CXCL10在腸組織勻漿中的含量

    從液氮中取出凍存結(jié)腸標(biāo)本,標(biāo)本融化后,剪碎,加入一定量預(yù)冷的PBS(pH 7.4),將標(biāo)本充分勻漿(制備成10%勻漿),離心20 min(4 ℃、2000 r/min),嚴(yán)格按ELISA試劑盒說明書進(jìn)行操作,運(yùn)用酶標(biāo)儀進(jìn)行檢測。

    2.8 16S rRNA基因V3~V4區(qū)測序

    采用Macherey-Nagel Kit試劑盒抽提腸內(nèi)容物DNA,通過Qubit Fluorometer對DNA濃度進(jìn)行定量,并通過瓊脂糖凝膠電泳對DNA完整性進(jìn)行評價,DNA質(zhì)量符合的樣本進(jìn)行文庫構(gòu)建。取30 ng DNA樣品及融合引物(341F:5’-ACTCCTACGGGA- GGCAGCAG-3’;806R:5’-GGACTACHVGGGTWT- CTAAT-3’)進(jìn)行PCR擴(kuò)增(16S rRNA基因V3~V4區(qū))。純化后溶于Elution Buffer,貼上標(biāo)簽,文庫構(gòu)建完成。采用IlluminaMiseqPE 300高通量測序平臺對檢測合格的文庫進(jìn)行雙末端(paired-end,PE)測序,獲得PE讀段(reads)。16S rRNA基因測序分析由武漢華大醫(yī)學(xué)檢驗(yàn)所有限公司高通量實(shí)驗(yàn)室完成。

    2.9 生物信息學(xué)分析

    下機(jī)數(shù)據(jù)經(jīng)過數(shù)據(jù)過濾,去除低質(zhì)量的reads,使用FLASH(Fast Length Adjustment of Short reads,v1.2.11)軟件,利用重疊關(guān)系將雙末端測序得到的成對reads拼接成原始序列(raw tags);raw tags經(jīng)進(jìn)一步去除嵌合體、短序列后得到優(yōu)質(zhì)序列(clean tags)。利用UPARSE軟件在97%相似度下將優(yōu)質(zhì)序列進(jìn)行聚類,獲得操作分類單元(operational taxonomic units,OTUs),通過RDP classifer軟件將OTU代表序列與數(shù)據(jù)庫Greengene比對,得到每個樣本的物種分類信息,并在各個水平(門,綱,目,科,屬,種)進(jìn)行物種注釋。基于聚類結(jié)果,進(jìn)行Alpha多樣性分析 [Chao1指數(shù)()、Ace指數(shù)(Ace)、香農(nóng)指數(shù)(Shannon index)、辛普森指數(shù)()]和Beta多樣性分析 [層級聚類樹分析(hierarchical clustering tree analysis)和主成分分析(principal component analysis,PCA)]。采用線性判別式分析 [linear discriminant analysis(LDA)effect size,LEfSe] 找到組間在豐度上有差異的物種,結(jié)果用LDA值分布柱狀圖及進(jìn)化分支圖表示,LDA值(log 10)>2的物種被認(rèn)為具有顯著性差異[17]。使用冗余分析(redundancy analysis,RDA)及Spearman系數(shù)分析腸道菌群與趨化因子水平的相關(guān)性。

    2.10 統(tǒng)計學(xué)方法處理

    3 結(jié)果

    3.1 一般情況觀察

    給藥3 d后,對照組小鼠食欲旺盛、外觀肥滿、被毛光澤、行動敏捷、呼吸均勻、大小便正常、體溫正常、體質(zhì)量無明顯變化。模型組食欲減退、被毛雜亂、豎毛明顯、弓背、扎堆、精神萎靡、倦怠懶動、呼吸短促、腹式呼吸明顯,小便正常,大便稀軟不成形,體溫降低、體質(zhì)量下降。與對照組比較,模型組小鼠體質(zhì)量差異有統(tǒng)計學(xué)意義(<0.01,圖1)。各藥物治療組小鼠呼吸漸平穩(wěn),被毛光澤度良好,大便濕潤成型,體質(zhì)量減輕后漸恢復(fù)。與模型組比較,各藥物治療組體質(zhì)量差異有統(tǒng)計學(xué)意義(<0.01,圖1)。給藥7 d后,對照組和各藥物組小鼠一般情況良好,模型組小鼠精神狀態(tài)好轉(zhuǎn),大便漸成形,體質(zhì)量增加,但與對照組和MXSGD組比較,差異仍有統(tǒng)計學(xué)意義(<0.05,圖1)。

    3.2 小鼠結(jié)腸組織病理變化

    光學(xué)顯微鏡下觀察各組小鼠結(jié)腸組織病理切片,結(jié)果見圖2。給藥3 d后,對照組結(jié)腸組織結(jié)構(gòu)完整,黏膜排列整齊,腺體無萎縮,黏膜下未見明顯炎性細(xì)胞浸潤。模型組結(jié)腸黏膜上皮水腫,部分腺體萎縮,黏膜固有層可見大量炎性細(xì)胞浸潤。經(jīng)藥物干預(yù)后,結(jié)腸組織損傷得到不同程度改善,固有層細(xì)胞浸潤明顯減輕。給藥7 d后,炎性細(xì)胞浸潤程度較給藥3 d后有所減輕,各藥物組結(jié)腸黏膜結(jié)構(gòu)完整,層次清晰,與對照組結(jié)構(gòu)相似。

    與對照組比較:*P<0.05 **P<0.01;與模型組比較:#P<0.05 ##P<0.01

    3.3 小鼠結(jié)腸組織CCL5、CXCL10蛋白表達(dá)

    由圖3、4可見,CCL5和CXCL10主要表達(dá)于胞漿或胞膜上。模型組結(jié)腸黏膜固有層中可見大量陽性表達(dá)細(xì)胞,對照組和各藥物組小鼠腸黏膜中也可見少量陽性表達(dá)細(xì)胞,尤其以CCL5表達(dá)更廣泛。由表1可知,給藥3 d后,模型組CCL5表達(dá)量較對照組顯著升高(<0.01),奧司他韋組、MXSGD組CCL5表達(dá)量較模型組顯著降低(<0.01);模型組CXCL10表達(dá)量較對照組顯著升高(<0.01),奧司他韋組、MXSGD組CXCL10表達(dá)量較模型組顯著降低(<0.01)。給藥7 d后,模型組CCL5表達(dá)量較對照組明顯升高(<0.05),MXSGD組CCL5表達(dá)量較模型組明顯降低(<0.05);各組CXCL10表達(dá)無顯著差異。由上述結(jié)果可知,正常情況下,小鼠結(jié)腸組織中少量表達(dá)CCL5和CXCL10。流感病毒感染可致結(jié)腸組織CCL5和CXCL10過量表達(dá)。經(jīng)藥物治療后,CCL5和CXCL10表達(dá)均有不同程度下降,與抗病毒顆粒比較,奧司他韋與MXSGD能更好地降低CCL5和CXCL10表達(dá)水平。給藥7 d后,模型組CCL5仍高于對照組(<0.05),MXSGD組CCL5表達(dá)低于模型組(<0.05)。

    圖2 各組小鼠結(jié)腸組織病理變化 (HE染色,×200)

    圖3 CCL5在各組小鼠結(jié)腸組織中的表達(dá) (PV-9000通用二步法,×200)

    圖4 CXCL10在各組小鼠結(jié)腸組織中的表達(dá) (PV-9000通用二步法,×200)

    表1 各組小鼠結(jié)腸組織CCL5、CXCL10平均光密度值(±s, n = 6)

    與對照組比較:*<0.05**<0.01;與模型組比較:#<0.05##<0.01,圖5同

    *< 0.05**< 0.01control group;#< 0.05##< 0.01model group, same as Figure 5

    3.4 小鼠結(jié)腸組織中CCL5、CXCL10的含量

    結(jié)果見圖5。給藥3 d后,模型組小鼠結(jié)腸組織中CCL5含量明顯高于對照組(<0.01);奧司他韋組、MXSGD組CCL5含量明顯低于模型組(<0.05);給藥7 d后,各組小鼠結(jié)腸組織中CCL5含量無顯著差異。給藥3 d后,模型組小鼠結(jié)腸組織中CXCL10含量明顯高于對照組(<0.01);奧司他韋組、抗病毒顆粒組、MXSGD組CXCL10含量明顯低于模型組(<0.05、0.01);給藥7 d后,各組小鼠結(jié)腸組織中CXCL10含量無顯著差異。

    圖5 各組小鼠結(jié)腸組織勻漿CCL5和CXCL10含量

    3.5 OTU聚類分析

    5組小鼠2個時間點(diǎn)共收集30個結(jié)腸內(nèi)容物樣品(每組每個時間點(diǎn)3個樣品),全部建庫合格,共獲得1 086 473條高質(zhì)量的tags,以97%相似度聚類,共獲得483個OTU,測序覆蓋深度(覆蓋指數(shù))0.998 8,平均每個樣品288個OTU,歸屬于9個門、50個屬。

    3.6 物種分布——微生物群落結(jié)構(gòu)分析

    3.6.1 門(Phylum)水平分布 給藥治療后,各組小鼠的腸道優(yōu)勢菌門構(gòu)成及相對豐度如圖6和表2所示。腸道菌群優(yōu)勢菌門以厚壁菌門(Firmicutes)和擬桿菌門(Bacteroidetes)為主。給藥3 d后,與對照組比較,模型組厚壁菌門相對豐度顯著降低(<0.01),變形菌門(Proteobacteria)相對豐度明顯增加(<0.05),厚壁菌門與擬桿菌門比值(F/B)降低,但差異不顯著;與模型組比較,奧司他韋組和MXSGD組厚壁菌門相對豐度顯著增加(<0.01),F(xiàn)/B值明顯增加(<0.05);各藥物組變形菌門相對豐度明顯降低(<0.05)。給藥7 d后,模型組疣微菌門(Verrucomicrobia)相對豐度較對照組明顯增加(<0.05);與模型組比較,奧司他韋組、抗病毒顆粒組、MXSGD組疣微菌門相對豐度明顯降低(<0.05、0.01),擬桿菌門相對豐度明顯升高(<0.05)。與給藥3 d后比較,給藥7 d后,模型組變形菌門相對豐度顯著下降(<0.01),F(xiàn)/B值明顯增加(<0.05)。

    3.6.2 屬水平分布 各組小鼠腸道菌群優(yōu)勢菌屬構(gòu)成及相對豐度如圖7及表3所示。給藥3 d后,正常菌群優(yōu)勢菌屬包括:(10.97%)、顫螺旋菌屬(6.76%)、艾克曼菌屬(6.27%)、(5.29%)、擬桿菌屬(2.62%)、乳桿菌屬(1.90%)、糞球菌屬(1.11%)、瘤胃球菌屬(3.46%)等;與對照組比較,模型組糞球菌屬、乳桿菌屬、普雷沃菌屬相對豐度明顯降低(<0.05),埃希菌屬()相對豐度明顯增加(<0.05);與模型組比較,奧司他韋組、抗病毒顆粒組和MXSGD組埃希菌屬相對豐度明顯降低(<0.05),抗病毒顆粒組顫螺旋菌屬相對豐度明顯增加(<0.05),MXSGD組糞球菌屬、乳酸桿菌屬相對豐度明顯增加(<0.05、0.01)。給藥7 d后,對照組優(yōu)勢菌屬包括:普雷沃菌屬(11.19%)、顫螺旋菌屬(8.01%)、乳桿菌屬(5.84%)、瘤胃球菌屬(2.00%)等;與對照組比較,模型組顫螺旋菌屬相對豐度明顯下降(<0.05);與給藥3 d后比較,給藥7 d后,模型組埃希菌屬相對豐度顯著下降(<0.01)。

    3.7 Alpha多樣性指數(shù)分析

    常用的Alpha多樣性指數(shù)包括Chao1指數(shù)、ACE指數(shù)、Shannon指數(shù)和Simpson指數(shù)。Chao1指數(shù)和ACE指數(shù)用來評價菌群的豐富度,數(shù)值越大,說明該菌群豐富度越高。Shannon指數(shù)和Simpson指數(shù)主要用來評價菌群多樣性,Shannon指數(shù)值越大,Simpson指數(shù)值越小,說明樣品的物種多樣性越高。給藥3 d和7 d后各組小鼠腸道菌群Alpha多樣性指數(shù)見表4。給藥3 d后,與對照組比較,模型組Chao1、ACE、Shannon指數(shù)降低,Simpson指數(shù)增高(<0.05、0.01);與模型組比較,抗病毒顆粒組Chao1、ACE指數(shù)升高,差異有統(tǒng)計學(xué)意義(<0.05)。奧司他韋組、抗病毒顆粒組、MXSGD組Shannon指數(shù)均顯著增高,Simpson指數(shù)均顯著降低(<0.01);給藥7 d后,與對照組比較,模型組Shannon指數(shù)明顯降低(<0.05);與模型組比較,抗病毒顆粒組Shannon指數(shù)明顯增高(<0.05)。上述結(jié)果提示,流感病毒感染4 d后(給藥3 d后),模型組小鼠腸道菌群豐富度和多樣性明顯降低,抗病毒顆粒能增加小鼠腸道菌群豐富度和多樣性;奧司他韋、MXSGD能增加小鼠腸道菌群多樣性;流感病毒感染8 d后(給藥7 d后),模型組小鼠腸道菌群多樣性仍低于對照組,抗病毒顆??梢栽黾有∈竽c道菌群多樣性。

    圖6 各組小鼠腸道菌群在門水平上分布

    表2 各組小鼠腸道菌群門水平相對豐度 (±s, n = 3)

    與對照組比較:*<0.05**<0.01;與模型組比較:#<0.05##<0.01;與給藥3 d后比較:^<0.05 ^^<0.01

    *< 0.05**< 0.01control group;#< 0.05##< 0.01model group; ^< 0.05 ^^< 0.01corresponding groups after 3 d of administration

    圖7 各組小鼠腸道菌群在屬水平分布

    表3 各組小鼠腸道菌群屬水平相對豐度(±s,n=3)

    與對照組比較:*<0.05**<0.01;與模型組比較:#<0.05##<0.01;與給藥治療3 d比較:^<0.05 ^^<0.01

    *< 0.05**< 0.01control group;#< 0.05##< 0.01model control group; ^< 0.05 ^^< 0.01corresponding groups after 3 d of administration

    3.8 Beta多樣性分析

    Beta多樣性分析主要對不同樣品/不同組間樣品微生物群落構(gòu)成進(jìn)行比較分析,常用PCA和層級聚類樹分析。PCA是常用的一種降維分析方法,通過將方差進(jìn)行分解,將多組數(shù)據(jù)的差異反映在二維坐標(biāo)圖上,圖中兩樣品點(diǎn)距離越近,表示它們的物種組成越相似,反之亦然。層級聚類樹可以清楚看出每個樣品的距離遠(yuǎn)近,根據(jù)樹枝的距離可以劃分出不同的分組。PCA結(jié)果(圖8-a)提示,給藥3 d后,模型組各樣品點(diǎn)與對照組樣品點(diǎn)距離較遠(yuǎn),表明兩組樣品之間細(xì)菌群落組成存在明顯差異,說明流感病毒感染可以影響正常小鼠腸道菌群的組成和結(jié)構(gòu)。藥物治療組各樣品點(diǎn)與對照組各樣品點(diǎn)距離接近,且明顯偏離于模型對照組各樣品點(diǎn),提示藥物組群落組成與對照組更相似。層級聚類樹分析結(jié)果(圖8-b)提示,模型組細(xì)菌群落結(jié)構(gòu)與其他各組有明顯差異。說明奧司他韋、抗病毒顆粒和MXSGD均可以調(diào)節(jié)流感病毒感染后小鼠腸道菌群組成結(jié)構(gòu)。給藥7 d后,各組PCA與層級聚類樹分析結(jié)果(圖8-c、d)提示,PCA結(jié)果顯示模型組樣品點(diǎn)仍然偏離于其余各組樣品點(diǎn),各藥物組樣品點(diǎn)與對照組樣品點(diǎn)間距離較給藥3 d后更接近。層級聚類樹分析結(jié)果顯示,MXSGD組與對照組腸道菌群組成更相似。

    表4 各組小鼠腸道菌群Alpha多樣性指數(shù)統(tǒng)計分析結(jié)果(±s,n=3)

    與對照組比較:*<0.05**<0.01;與模型組比較:#<0.05##<0.01

    *< 0.05**< 0.01control group;#< 0.05##< 0.01model group

    a、c-PCA b、d-層級聚類樹分析 A1~A3、F1~F3-對照組 B1~B3、G1~G3-模型組 C1~C3、H1~H3-奧司他韋組 D1~D3、I1~I(xiàn)3-抗病毒顆粒組 E1~E3、J1~J3-MXSGD組

    3.9 組間差異物種的篩選

    采用LEfSe篩選組間具有顯著差異的物種。LDA值超過設(shè)定值(常以>3.5為篩選標(biāo)準(zhǔn))的物種,即被認(rèn)為是組間具有統(tǒng)計學(xué)差異的生物標(biāo)記物(biomarker),并根據(jù)差異物種繪制物種進(jìn)化分支圖。在進(jìn)化分支圖中,由內(nèi)至外輻射的圓圈代表由界至屬(或種)的分類級別。在不同分類級別上的每一個小圓圈代表該水平下的一個分類,小圓圈直徑大小與相對豐度大小呈正比。無顯著差異的物種統(tǒng)一著色為黃色,差異物種生物標(biāo)記物跟隨組進(jìn)行著色。圖9-a~h顯示,給藥3 d后,模型組與其他各組比較具有顯著差異的優(yōu)勢物種為變形菌門、γ-變形桿菌綱(Gammaproteobacteria)、腸桿菌目(Enterobacteriales)、腸桿菌科(Enterobacteriaceae)、埃希菌屬、克雷伯菌屬、沙雷氏菌屬、變形桿菌屬,疣微菌門、疣微菌綱(Verrucomicrobiae)、脫鐵桿菌目(Deferribacterales)、脫鐵桿菌屬、多拉菌屬等;與模型組比較,對照組具有顯著優(yōu)勢的物種為厚壁菌門、梭菌綱(Clostridia)、瘤胃菌科(Ruminococcaceae)、瘤胃菌屬、顫螺旋菌屬、毛螺菌科(Lachnospiraceae)、糞球菌屬、紅蝽菌綱(Coriobacteriia)、紅蝽菌目(Coriobacteriales)、紅蝽菌科(Coriobacteriaceae)、羅斯氏菌屬等;奧司他韋組具有顯著優(yōu)勢的物種為厚壁菌門、梭菌綱、瘤胃球菌科、顫螺旋菌屬、乳桿菌目(Lactobacillales)等;抗病毒顆粒組具有顯著優(yōu)勢的物種為厚壁菌門、梭菌綱、瘤胃球菌科、顫螺旋菌屬;MXSGD組具有顯著優(yōu)勢的物種為普雷沃菌科(Prevotellaceae)、普雷沃菌屬、厚壁菌門、梭菌綱、瘤胃球菌科、顫螺旋菌屬、乳桿菌目、乳桿菌科(Lactobacillaceae)、乳桿菌屬、毛螺菌科、糞球菌屬、丁酸梭菌屬等。給藥7 d天后,各組間顯著差異物種不明顯(LDA值<3.5)。

    a-對照組和模型組LDA值分布柱狀圖 b-對照組和模型組進(jìn)化分支圖 c-模型組和奧司他韋組LDA值分布柱狀圖 d-模型組和奧司他韋組進(jìn)化分支圖 e-模型組和抗病毒顆粒組LDA值分布柱狀圖 f-模型組和抗病毒顆粒組進(jìn)化分支圖 g-模型組和MXSGD組LDA值分布柱狀圖 h-模型組和MXSGD組進(jìn)化分支圖

    a-histogram of LDA scores in normal control group and model control group b-cladogram in normal control group and model control group c- histograms of LDA scores in model control group and oseltamivir group d-cladogram in model control group and oseltamivir group e-histograms of LDA scores in model control group and antivirus granule group f-cladogram in model control group and antivirus granule group g-histograms of LDA scores in model control group and MXSGD group h-cladogram in model control group and MXSGD group

    圖9 各組小鼠腸道菌群LEfSE差異分析結(jié)果

    Fig. 9 LEfSE difference analysis in various groups

    3.10 腸道菌群與趨化因子關(guān)聯(lián)性分析

    通過RDA及Spearman等級相關(guān)系數(shù)分析各組小鼠結(jié)腸組織中趨化因子CCL5、CXCL10與腸道菌群的關(guān)聯(lián)性。結(jié)果(圖10)顯示,小鼠腸道菌群的分布與腸黏膜趨化因子CCL5和CXCL10含量的相關(guān)性具有顯著性差異。在門水平上:變形菌門、疣微菌門、脫鐵桿菌門、擬桿菌門與CCL5和CXCL10含量呈正相關(guān);厚壁菌門、放線菌門(Actinobacteria)與CCL5和CXCL10含量呈負(fù)相關(guān)。在屬水平上:埃希菌屬、梭菌屬、變形桿菌屬、克雷伯菌屬、沙雷氏菌屬、擬桿菌屬、屬與CCL5和CXCL10含量呈正相關(guān);乳桿菌屬、厭氧棒狀菌屬、糞球菌屬、顫螺旋菌屬、普雷沃菌屬、瘤胃球菌屬、丁酸梭菌屬、低嗜鹽細(xì)菌屬與CCL5和CXCL10含量呈負(fù)相關(guān)。其中,埃希菌屬、克雷伯菌屬、梭菌屬、變形菌屬、沙雷氏菌屬、糞球菌屬、厭氧棒狀菌屬與CCL5和CXCL10相關(guān)性顯著(<0.05、0.01);低嗜鹽細(xì)菌屬、擬桿菌屬、乳桿菌屬與CCL5相關(guān)性顯著(<0.05、0.01);顫螺旋菌屬與CXCL10相關(guān)性顯著(<0.05)。

    a-RDA分析(門水平) b-RDA分析(屬水平) c-Spearman相關(guān)熱圖 A1~A3-對照組 B1~B3-模型組 C1~C3-奧司他韋組 D1~D3-抗病毒顆粒組 E1~E3-MXSGD組

    4 討論

    MXSGD為漢代張仲景《傷寒論》名方,具有宣肺平喘、祛除寒濕、通腑瀉熱兼燥濕健脾之功效,是治療證屬熱毒襲肺和熱毒壅肺流感的基本方[18]。2019年12月以來爆發(fā)的新型冠狀病毒(severe acute respiratory syndrome coronavirus 2,SARS-CoV-2)感染肺炎(corona virus disease 2019,COVID-19)疫情中,MXSGD是國家及各地方衛(wèi)生管理部門發(fā)布的診療方案中被推薦使用頻次最高的中藥方劑之一[19],是國家衛(wèi)生健康委員會發(fā)布的系列《新型冠狀病毒的肺炎診療方案》中推薦在臨床治療期使用的“清肺排毒湯”的重要組成部分[20]。

    前期研究發(fā)現(xiàn),趨化因子可能是流感“肺病及腸”的物質(zhì)基礎(chǔ)之一[11]。流感病毒感染的單核/巨噬細(xì)胞可迅速產(chǎn)生CC家族趨化因子(如CCL5),以及CXC家族趨化因子(如CXCL10),這些趨化因子的急劇升高與流感病毒感染引起的炎癥病理損傷密切相關(guān)。CCL5也被命名為RANTES(regulated upon activation normal T cell expressed and secreted),由上皮細(xì)胞、CD8+T、CD4+T、單核巨噬細(xì)胞產(chǎn)生,對單核細(xì)胞、嗜堿性粒細(xì)胞和T細(xì)胞有強(qiáng)烈的趨化作用,能促進(jìn)白細(xì)胞向炎癥部位遷移和浸潤。CXCL-10又被稱作γ干擾素誘導(dǎo)蛋白(interferon-gamma-induced protein,IP-10),是一種屬于CXC趨化因子家族的小分子蛋白質(zhì),與流感病毒感染引起的肺部損傷的嚴(yán)重程度密切相關(guān),由IFN-γ、腫瘤壞死因子-α(tumor necrosis factor,TNF-α)、病毒RNA等誘導(dǎo)上皮細(xì)胞、內(nèi)皮細(xì)胞、單核巨噬細(xì)胞、成纖維細(xì)胞和淋巴細(xì)胞產(chǎn)生,可趨化T淋巴細(xì)胞、自然殺傷細(xì)胞和單核細(xì)胞聚集,促進(jìn)IFN-γ、IL-6、IL-8等多種促炎細(xì)胞因子分泌[21]。因此,抑制CCL5和CXCL10過表達(dá),對于控制流感病毒感染引起的炎癥反應(yīng)有重要意義。

    已有研究發(fā)現(xiàn),腸道菌群及其代謝產(chǎn)物可影響宿主局部或全身免疫功能。短鏈脂肪酸(short chain fatty acids,SCFAs)是腸道有益菌發(fā)酵膳食纖維的重要產(chǎn)物,主要包括乙酸、丙酸、丁酸等。已有研究發(fā)現(xiàn),厚壁菌門是丁酸鹽的重要來源,丁酸鹽是細(xì)胞的主要能量底物,厚壁菌門和擬桿菌門比值降低可以直接影響腸道菌群對膳食纖維代謝,使SCFAs濃度降低,厚壁菌門減少可誘發(fā)或加重局部炎癥反應(yīng)[22]。變形菌門是腸道的4個主要菌門(厚壁菌門、擬桿菌門、變形菌門、放線菌門)中隨時間變化最不穩(wěn)定的菌門,變形菌門豐度增加被視為反映腸道菌群紊亂的一個標(biāo)志及潛在的疾病診斷指標(biāo)[23]。變形菌門能夠有效利用炎癥反應(yīng)產(chǎn)生的硝酸鹽作為電子受體,進(jìn)行厭氧呼吸,因此在炎癥環(huán)境中,比依賴發(fā)酵生長的厚壁菌門和擬桿菌門更具增殖優(yōu)勢,其過度生長可進(jìn)一步加重局部炎癥反應(yīng)[23]。研究發(fā)現(xiàn),SCFAs可以誘導(dǎo)調(diào)節(jié)性T細(xì)胞(regulatory cells,Tregs)的產(chǎn)生和分化,從而促進(jìn)抑炎細(xì)胞因子IL-10產(chǎn)生,并且可以強(qiáng)烈減少促炎趨化因子CCL3、CCL4、CCL5、CXCL9、CXCL10、CXCL11的釋放,從而抑制白細(xì)胞遷徙,具有很強(qiáng)的抗炎作用[24]。普雷沃菌屬、瘤胃球菌屬、糞球菌屬、丁酸梭菌屬、乳酸桿菌屬均可增加腸道內(nèi)SCFAs含量,誘導(dǎo)小鼠腸道中Treg細(xì)胞的產(chǎn)生和分化,促進(jìn)IL-10的產(chǎn)生,發(fā)揮炎癥抑制作用[25-27]。顫螺旋菌屬在炎癥疾病中豐度降低,與健康呈正相關(guān),對人體健康具有潛在的重要性[28]。而過度生長的大腸埃希菌或肺炎克雷伯菌等可以促進(jìn)CCL5及CXCL10等趨化因子的表達(dá),從而加重炎性損傷程度[7-8]。

    研究發(fā)現(xiàn),流感病毒感染后模型組小鼠腸黏膜出現(xiàn)明顯病理變化,黏膜下層和固有層中見大量炎性細(xì)胞浸潤,炎性細(xì)胞胞漿和胞膜上有大量趨化因子CCL5和CXCL10表達(dá)。Alpha多樣性分析結(jié)果顯示模型組小鼠腸道菌群豐富度和多樣性明顯降低。Beta多樣性分析結(jié)果顯示模型組小鼠腸道菌群組成結(jié)構(gòu)與對照組和各藥物組明顯不同。病毒感染后第4天,與對照組比較,模型組厚壁菌門相對豐度明顯下降,厚壁菌門與擬桿菌門比值(F/B)明顯降低,變形菌門相對豐度顯著增加,埃希菌屬相對豐度增加,糞球菌屬、乳桿菌屬、普雷沃菌屬豐度下降。LEfSe分析結(jié)果顯示,克雷伯菌屬也是模型組的優(yōu)勢菌屬。經(jīng)藥物治療3 d后,各藥物組腸黏膜炎性細(xì)胞浸潤明顯減少,CCL5及CXCL10表達(dá)明顯降低;Alpha多樣性分析結(jié)果顯示各藥物組腸道菌群多樣性增加,抗病毒顆粒組菌群豐富度增加;Beta多樣性分析結(jié)果顯示,在各藥物組中,MXSGD組小鼠腸道菌群組成更接近對照組;與模型組比較,各藥物組變形菌門、埃希菌屬相對豐度均降低,奧司他韋組和MXSGD組厚壁菌門相對豐度及F/B值明顯增加,MXSGD組乳桿菌屬和糞球菌屬相對豐度顯著增加;此外,LEfSe分析的結(jié)果提示,MXSGD組的優(yōu)勢物種還包括瘤胃球菌科、顫螺旋菌屬、毛螺菌科、丁酸梭菌屬等。通過相關(guān)性分析進(jìn)一步發(fā)現(xiàn),埃希菌屬、克雷伯菌屬、變形桿菌屬等與CCL5和CXCL10含量呈顯著正相關(guān),糞球菌屬、厭氧棒狀菌等與CCL5和CXCL10含量呈顯著負(fù)相關(guān),乳桿菌屬與CCL5含量呈顯著負(fù)相關(guān),顫螺旋菌屬與CXCL10含量呈顯著負(fù)相關(guān)。CCL5和CXCL10表達(dá)下調(diào),也可能會影響抗病毒細(xì)胞因子IFN-γ產(chǎn)生,從而影響機(jī)體抗病毒感染免疫。但前期研究[29]發(fā)現(xiàn),MXSGD可以促進(jìn)細(xì)胞因子IL-2產(chǎn)生,IL-2可以刺激T細(xì)胞快速增殖分化產(chǎn)生IFN-γ。可見,MXSGD可以通過調(diào)節(jié)復(fù)雜的細(xì)胞因子網(wǎng)絡(luò),減少炎癥因子表達(dá)下調(diào)對機(jī)體抗感染免疫帶來的負(fù)面影響。

    綜上所述,流感病毒感染后,小鼠腸道菌群組成結(jié)構(gòu)明顯紊亂,且結(jié)腸組織中趨化因子CCL5和CXCL10的表達(dá)明顯增強(qiáng)。相關(guān)性分析結(jié)果顯示,腸道菌群與趨化因子CCL5和CXCL10表達(dá)顯著相關(guān)。奧司他韋、抗病毒顆粒和MXSGD均能夠不同程度修復(fù)流感病毒感染后紊亂的腸道菌群,并降低結(jié)腸組織中趨化因子CCL5和CXCL10的表達(dá)。與其他藥物組比較,MXSGD組腸道菌群結(jié)構(gòu)更接近對照組水平,MXSGD能更好促進(jìn)乳桿菌屬和糞球菌屬益生菌生長,且對流感病毒誘導(dǎo)的CCL5和CXCL10過表達(dá)均有一定的抑制作用。因此,腸道菌群可以作為MXSGD防治流感的一個潛在靶點(diǎn),MXSGD可能通過調(diào)節(jié)菌群結(jié)構(gòu)進(jìn)而下調(diào)趨化因子水平,緩解流感病毒感染所致免疫炎癥損傷,其具體作用機(jī)制值得進(jìn)一步深入研究。

    利益沖突 所有作者均聲明不存在利益沖突

    [1] Laetitia M, Charrel R N, Pierre-Emmanuel C,. Prevalence of gastrointestinal symptoms in patients with influenza, clinical significance, and pathophysiology of human influenza viruses in faecal samples: what do we know? [J]., 2015, 12: 215-224.

    [2] Wang J, Li F, Wei H,. Respiratory influenza virus infection induces intestinal immune injury via microbiota-mediated Th17 cell-dependent inflammation [J]., 2014, 211(12): 2397-2410.

    [3] 劉雨琪, 朱俊萍, 何秋水. 流感病毒與腸道菌群的相互作用對感染結(jié)局的影響 [J]. 醫(yī)學(xué)綜述, 2019, 25(14): 2746-2758.

    [4] Li H, Liu X, Chen F,. Avian influenza virus subtype H9N2 affects intestinal microbio-ta, barrier structure injury, and inflammatory intestinal disease in the chicken ileum [J]., 2018, 10(5): 270-284.

    [5] Tavares L P, Garcia C C, Gon?alves A,. ACKR2 contributes to pulmonary dysfunction by shaping CCL5: CCR5-dependent recruitment of lymphocytes during influenza A infection in mice [J]., 2020, 318(4): L655-L670.

    [6] 胡艮娣, 王曉泉, 陳祥, 等. 趨化因子IP-10在介導(dǎo)病毒感染性疾病中的研究進(jìn)展 [J]. 病毒學(xué)報, 2019, 35(4): 672-678.

    [7] Karlsson I, Hagman R, Guo Y,. Pathogenicand lipopolysaccharide enhance the expression of IL-8, CXCL5, and CXCL10 in canine endometrial stromal cells [J]., 2015, 84(1): 34-42.

    [8] 楊明, 麻雅婷, 何賞, 等. LIF和RANTES在單病原菌性血流感染小鼠模型中的表達(dá)及意義 [J]. 臨床檢驗(yàn)雜志, 2018, 36(1): 53-56.

    [9] 李玲, 魏科, 盧芳國, 等. 基于TLR4-MyD88-TRAF6信號通路的麻杏石甘湯抗A型流感病毒感染小鼠所致的病毒性肺損傷研究 [J]. 中草藥, 2017, 48(8): 1591-1596.

    [10] 張世鷹, 何谷良, 盧芳國, 等. 基于TLR7/8介導(dǎo)的IFN-α/β蛋白表達(dá)水平探討麻黃先煎之麻杏石甘湯抗流感病毒的機(jī)制 [J]. 中華中醫(yī)藥雜志, 2019, 34(3): 1188-1193.

    [11] 鄒莉, 何谷良, 盧芳國, 等. 麻杏石甘湯對流感病毒肺部感染模型肺、結(jié)腸組織MCP-1蛋白表達(dá)水平的影響 [J]. 中國實(shí)驗(yàn)方劑學(xué)雜志, 2018, 24(5): 100-106.

    [12] 郭元吉, 程小雯. 流行性感冒病毒及其試驗(yàn)技術(shù) [M]. 北京: 中國三峽出版社, 1997: 129-130.

    [13] 葛資宇, 童驕, 那婧婧, 等. 不同煎煮方法的麻杏石甘湯及其含藥血清對A型流感病毒神經(jīng)氨酸酶活性的影響 [J]. 中國中西醫(yī)結(jié)合雜志, 2016, 36(9): 1119-1123.

    [14] 饒毅, 張五萍, 魏惠珍, 等. 不同煎煮方法對麻杏石甘湯中成分變化研究 [J]. 中成藥, 2014, 36(2): 337-341.

    [15] 李玲, 盧芳國, 伍參榮, 等. 麻杏石甘湯對A型流感病毒感染小鼠免疫保護(hù)作用的研究 [J]. 中醫(yī)藥學(xué)報, 2010, 38(2): 25-28.

    [16] 賀石林, 王鍵, 王凈凈. 中醫(yī)科研設(shè)計與統(tǒng)計學(xué) [M]. 長沙: 湖南科學(xué)技術(shù)出版社2013: 48-49.

    [17] Segata N, Izard J, Waldron L,. Metagenomic biomarker discovery and explanation [J]., 2011, 12(6): R60.

    [18] 趙靜, 郭洪濤, 韓經(jīng)丹, 等. 中醫(yī)藥治療流行性感冒文本挖掘結(jié)果與診療方案的比較分析 [J]. 中醫(yī)雜志, 2014, 55(7): 612-616.

    [19] 王哲義, 懌澤, 曲稔棟, 等. 基于網(wǎng)絡(luò)藥理學(xué)的麻杏石甘湯治療新型冠狀病毒肺炎(COVID-19)作用機(jī)制探討 [J]. 中草藥, 2020, 51(8): 1996-2003.

    [20] 李修洋, 宋斌, 雷燁, 等.《新型冠狀病毒肺炎診療方案(試行第六版)》中醫(yī)診療方案解讀 [J]. 吉林中醫(yī)藥, 2020, 40(6): 701-708.

    [21] Wang W, Yang P, Zhong Y,. Monoclonal antibody against CXCL-10/IP-10 ameliorates influ enza A (H1N1) virus induced acute lung injury [J]., 2013, 23(4): 577-580.

    [22] Kho Z Y, Lal S K. The human gut microbiome-A potential controller of wellness and disease [J]., 2018, doi: 10.3389/fmicb.2018.01835.

    [23] Shin N R, Whon T W, Bae J W. Proteobacteria: Microbial signature of dysbiosis in gut microbiota [J]., 2015, 33(9): 496-503.

    [24] Nastasi C, Candela M, Bonefeld C M,. The effect of short-chain fatty acids on human monocyte-derived dendritic cells [J]., 2015, 5: 16148.

    [25] Koh A, De Vadder F, Kovatcheva-Datchary P,. From dietary fiber to host physiology: Short-chain fatty acids as key bacterial metabolites [J]., 2016, 165(6): 1332-1345.

    [26] Louis P, Hold G L, Flint H J. The gut microbiota, bacterial metabolites and colorectal cancer [J]., 2014, 12(10): 661-672.

    [27] 俞昊男, 劉志華. 腸道微生物與黏膜免疫研究的前沿進(jìn)展 [J]. 中國免疫學(xué)雜志, 2019, 35(16): 1921-1930.

    [28] Konikoff T, Gophna U. Oscillospira: A central, enigmatic component of the human gut microbiota [J]., 2016, 24(7): 523-524.

    [29] 李玲, 吳佳敏, 歐陽建軍, 等. 抗流感病毒性肺炎的有效中藥復(fù)方篩選及機(jī)制研究 [J]. 中國免疫學(xué)雜志, 2018, 34(8): 1168-1173.

    Effects of Maxing Shigan Decoction on intestinal flora and chemokines CCL5 and CXCL10 in mice infected with influenza virus

    WANG Ping1, 2, ZHAO Cheng1, LU Fang-guo1, WU Tao1, ZHANG Xiang-gang1, CHEN Chun-jing1, XIAO Rong1, NING Yi1, WEI Ke1, LI Ling1, CAI Kun2, TIAN Wei-yi2

    1. Hunan University of Chinese Medicine, Changsha 410208, China 2. Guizhou University of Chinese Medicine, Guiyang 550025, China

    To investigate the potential mechanism of Maxing Shigan Decoction (MXSGD, 麻杏石甘湯) by affecting the intestinal flora and chemokines in mice in order to prevent and treat the influenza virus infection.The infected mice model of influenza A virus was established by intranasal inoculation. After 3 and 7 d of gavage administration or saline, the state of mice was observed. Colon tissue was detected by HE staining. The expressions and contents of chemokines CCL5 and CXCL10 were detected by immunohistochemistry and enzyme-linked immunosorbent assay, respectively. The bacteria in colonic contents was sequenced, annotated and clustered by using the V3—V4 variable region of 16S rRNA. The Alpha diversity, Beta diversity, and the species difference among groups were analyzed by linear discriminant analysis (LDA) effect size, LEfSe. Redundancy analysis and Spearman’s rank correlation coefficient were used for correlation analysis between intestinal flora and CCL5, CXCL10.Compared with the normal group, the weight of mice in the model control group was decreased (< 0.01), the inflammatory cell infiltration in the lamina propria of intestinal mucosa was obvious and the expression of CCL5 and CXCL10 was upregulated (< 0.01), the levels of CCL5 and CXCL10 in the colon tissue were both increased (< 0.01) after 3 d treatment. Compared with the model control group, the weight of mice was increased in each treatment group. Inflammatory cell infiltration, CCL5 and CXCL10 positive cells, the contents of CCL5 and CXCL10 in intestinal mucosa were all decreased, especially in oseltamivir group and MXSGD group (< 0.01,< 0.05). There were significant differences in the composition of the gut microbiota among groups according to the analysis of flora diversity and the species differences. The composition of intestinal flora in the model control group was significantly different from the normal control group and each drug group. The relative abundances of Proteobacteria andin the model control group were significantly higher than the normal control group (< 0.05), while Firmicutes,andwere significantly lower (< 0.05). Compared with the model control group, the relative abundances of Proteobacteria andin each drug group were significantly lower (< 0.05), the relative abundance of Firmicutes in oseltamivir group and MXSGD group were increased significantly (< 0.01), the relative abundance ofin antiviral granule group was increased (< 0.05), and the relative abundance ofandin MXSGD group was increased significantly (< 0.01,< 0.05). After 7 d, pathological changes of colon tissue, the expression of chemokines and the abundance of Proteobacteria andspecies (< 0.01) in the model control group were all decreased. Compared with the model control group, the relative abundances of Verrucomicrobia in each drug group were significantly lower (< 0.01,< 0.05) and the relative abundance of Bacteroidetes were increased significantly (< 0.05). Correlation analysis showed that CCL5 and CXCL10 contents in the colon of mice were significantly correlated with the abundance of,,and(< 0.01,< 0.05).The infection of influenza A virus can make the intestinal flora structure disorganized and immune function in mice unbalanced. MXSGD may regulate the intestinal flora structure and then affect the production of chemokines, which has a certain protective effect on the intestinal immune damage caused by influenza virus.

    influenza A virus (IAV); Maxing Shigan Decoction; intestinal flora; chemokines; CCL5; CXCL10

    R285

    A

    0253 - 2670(2021)01 - 0160 - 16

    10.7501/j.issn.0253-2670.2021.01.020

    2020-04-20

    國家自然科學(xué)基金資助項目(81774126);國家自然科學(xué)基金資助項目(82074250);湖南省自然科學(xué)基金項目(2020JJ4063);湖南省教育廳創(chuàng)新平臺開放基金項目(17K067);湖南中醫(yī)藥大學(xué)中西醫(yī)結(jié)合一流學(xué)科開放基金(2018ZXYJH11);湖南中醫(yī)藥大學(xué)研究生創(chuàng)新課題立項項目(2019CX05);湖南中醫(yī)藥大學(xué)研究生創(chuàng)新課題立項項目(2019CX22);貴州省科技創(chuàng)新人才團(tuán)隊(黔科合平臺人才 [2020] 5010)

    王 平(1977—),女,貴州畢節(jié)人,副教授,湖南中醫(yī)藥大學(xué)在讀博士,研究方向?yàn)楦腥拘约膊〉闹嗅t(yī)藥防治研究。

    盧芳國,女,教授,博士研究生導(dǎo)師,研究方向?yàn)楦腥拘约膊〉闹嗅t(yī)藥防治研究。E-mail: lufangguo0731@163.com

    [責(zé)任編輯 潘明佳]

    猜你喜歡
    奧司趨化因子菌門
    持續(xù)吃奧司他韋可預(yù)防流感?
    科教新報(2025年5期)2025-01-14 00:00:00
    合理用藥—奧司他韋知多少
    藥物與人(2023年3期)2024-01-03 21:50:04
    野生樹鼩與人工飼養(yǎng)樹鼩消化道不同部位微生物組成的比較研究
    饑餓與重攝食對河蟹腸道菌群結(jié)構(gòu)的影響
    昆蟲體內(nèi)微生物多樣性的影響因素研究進(jìn)展
    妊娠期糖尿病腸道菌群變化及臨床價值分析
    趨化因子及其受體在腫瘤免疫中調(diào)節(jié)作用的新進(jìn)展
    不同劑量奧司他韋治療流行性感冒疑似患兒的臨床療效比較
    肝細(xì)胞癌患者血清趨化因子CXCR12和SA的表達(dá)及臨床意義
    奧司他韋治療小兒流行性感冒75例療效及對血清細(xì)胞因子的影響
    中文字幕免费在线视频6| 韩国高清视频一区二区三区| 熟妇人妻不卡中文字幕| 亚洲人成网站在线观看播放| 亚洲精品乱码久久久久久按摩| 日韩大片免费观看网站| av卡一久久| 国产精品国产三级国产专区5o| 日本三级黄在线观看| 色视频www国产| 久久99蜜桃精品久久| 亚洲第一区二区三区不卡| 最近中文字幕2019免费版| 亚洲经典国产精华液单| 久久女婷五月综合色啪小说 | 制服丝袜香蕉在线| 国产欧美另类精品又又久久亚洲欧美| 卡戴珊不雅视频在线播放| 久久久久性生活片| 99re6热这里在线精品视频| 神马国产精品三级电影在线观看| av一本久久久久| 五月天丁香电影| 久久影院123| 久久久久网色| a级一级毛片免费在线观看| 我的老师免费观看完整版| 亚洲精品乱久久久久久| 国产老妇伦熟女老妇高清| 免费黄频网站在线观看国产| 亚洲欧美精品自产自拍| 国产一区二区亚洲精品在线观看| 久久人人爽人人片av| 秋霞伦理黄片| 免费高清在线观看视频在线观看| 成人特级av手机在线观看| 色婷婷久久久亚洲欧美| 亚洲av成人精品一区久久| 国产av码专区亚洲av| 精品一区二区免费观看| 免费观看的影片在线观看| 乱码一卡2卡4卡精品| 亚洲成人av在线免费| 国产伦精品一区二区三区视频9| 69av精品久久久久久| 欧美成人精品欧美一级黄| 亚洲无线观看免费| 美女主播在线视频| 丝袜美腿在线中文| videos熟女内射| 精品国产一区二区三区久久久樱花 | 国产老妇女一区| 国产伦精品一区二区三区视频9| 欧美成人a在线观看| 熟妇人妻不卡中文字幕| 免费观看av网站的网址| 欧美另类一区| 欧美xxxx性猛交bbbb| 80岁老熟妇乱子伦牲交| 中文欧美无线码| 一级毛片久久久久久久久女| 欧美成人午夜免费资源| 日本免费在线观看一区| 2018国产大陆天天弄谢| 99热这里只有是精品在线观看| 国产精品不卡视频一区二区| 3wmmmm亚洲av在线观看| 97热精品久久久久久| 中文字幕人妻熟人妻熟丝袜美| 精品国产一区二区三区久久久樱花 | 国产午夜精品久久久久久一区二区三区| 听说在线观看完整版免费高清| 国产乱来视频区| 亚洲无线观看免费| 国产精品秋霞免费鲁丝片| 亚洲精品国产成人久久av| 国内少妇人妻偷人精品xxx网站| 国产毛片在线视频| 日本熟妇午夜| 色网站视频免费| 日韩av在线免费看完整版不卡| 国内揄拍国产精品人妻在线| 国产亚洲91精品色在线| 国产精品女同一区二区软件| 久久午夜福利片| 久久久久精品性色| 成人二区视频| 成人美女网站在线观看视频| 不卡视频在线观看欧美| 亚洲欧美中文字幕日韩二区| 国产精品久久久久久精品电影小说 | 国产 一区精品| 国产精品国产三级国产av玫瑰| 午夜亚洲福利在线播放| 麻豆成人午夜福利视频| 午夜爱爱视频在线播放| 久久99热这里只有精品18| 亚洲精品色激情综合| 99九九线精品视频在线观看视频| 美女内射精品一级片tv| 亚洲综合精品二区| 欧美一区二区亚洲| 国产av不卡久久| 亚洲国产av新网站| 女人久久www免费人成看片| 99热这里只有是精品在线观看| 久久99热这里只频精品6学生| 亚洲色图综合在线观看| 亚洲av.av天堂| 亚洲精品视频女| 久久久久久久久久久免费av| 内地一区二区视频在线| 18禁在线无遮挡免费观看视频| 久久人人爽人人爽人人片va| 只有这里有精品99| 亚洲成人av在线免费| 黄色一级大片看看| 交换朋友夫妻互换小说| 精品久久久久久久人妻蜜臀av| 高清视频免费观看一区二区| 国产永久视频网站| 精品久久久久久久久av| 男插女下体视频免费在线播放| av线在线观看网站| av在线亚洲专区| 男的添女的下面高潮视频| 国产成年人精品一区二区| 亚洲va在线va天堂va国产| 美女高潮的动态| 日韩在线高清观看一区二区三区| 两个人的视频大全免费| 国产成年人精品一区二区| 国精品久久久久久国模美| 国产黄色视频一区二区在线观看| 国产色婷婷99| 亚洲国产色片| 亚洲国产欧美人成| 国产精品嫩草影院av在线观看| 亚洲美女搞黄在线观看| 中文字幕亚洲精品专区| 久热这里只有精品99| 国产成人精品福利久久| 日韩一本色道免费dvd| 免费av毛片视频| 精品午夜福利在线看| 有码 亚洲区| 视频区图区小说| 国产午夜福利久久久久久| 中文资源天堂在线| 亚洲成人一二三区av| 能在线免费看毛片的网站| 国产一区二区亚洲精品在线观看| 亚洲av成人精品一二三区| 亚洲精品日韩在线中文字幕| 亚洲美女视频黄频| 舔av片在线| 成人国产av品久久久| 噜噜噜噜噜久久久久久91| 一级毛片 在线播放| 亚洲aⅴ乱码一区二区在线播放| 日本一二三区视频观看| 欧美成人精品欧美一级黄| 中文字幕久久专区| 欧美+日韩+精品| 日韩一本色道免费dvd| 亚洲av免费在线观看| 18禁裸乳无遮挡免费网站照片| 成年av动漫网址| 精品熟女少妇av免费看| 69人妻影院| 真实男女啪啪啪动态图| 91aial.com中文字幕在线观看| 91精品伊人久久大香线蕉| 国产精品麻豆人妻色哟哟久久| 亚洲一区二区三区欧美精品 | 高清午夜精品一区二区三区| 少妇裸体淫交视频免费看高清| 亚洲精品一区蜜桃| 街头女战士在线观看网站| 欧美日韩国产mv在线观看视频 | 婷婷色麻豆天堂久久| 在线观看免费高清a一片| 国产成人精品一,二区| 汤姆久久久久久久影院中文字幕| 一区二区三区乱码不卡18| 少妇的逼水好多| 26uuu在线亚洲综合色| 日本欧美国产在线视频| av专区在线播放| 国内精品美女久久久久久| 一级片'在线观看视频| 亚洲一级一片aⅴ在线观看| 欧美亚洲 丝袜 人妻 在线| 18禁裸乳无遮挡动漫免费视频 | 晚上一个人看的免费电影| 欧美高清成人免费视频www| 九九在线视频观看精品| 精品久久久久久电影网| 亚洲欧美日韩卡通动漫| 国产精品伦人一区二区| 中国美白少妇内射xxxbb| 久久精品久久久久久噜噜老黄| 国产亚洲精品久久久com| 国产伦理片在线播放av一区| 国产男女内射视频| 新久久久久国产一级毛片| 亚洲av国产av综合av卡| 在线 av 中文字幕| 国产淫语在线视频| 亚洲av日韩在线播放| 精品午夜福利在线看| 最近手机中文字幕大全| 在线观看一区二区三区激情| 三级男女做爰猛烈吃奶摸视频| 免费观看a级毛片全部| 午夜亚洲福利在线播放| 成人二区视频| 亚洲国产日韩一区二区| 亚洲欧美成人综合另类久久久| 简卡轻食公司| 在线a可以看的网站| 噜噜噜噜噜久久久久久91| 中文字幕人妻熟人妻熟丝袜美| 久久99精品国语久久久| 日本欧美国产在线视频| av天堂中文字幕网| 91在线精品国自产拍蜜月| 亚洲三级黄色毛片| 日本wwww免费看| 亚洲精品成人av观看孕妇| 免费av不卡在线播放| 成人漫画全彩无遮挡| 男女那种视频在线观看| 精品久久久久久久久亚洲| 亚洲国产精品成人久久小说| 插逼视频在线观看| 亚洲国产精品成人综合色| 午夜视频国产福利| 亚洲精品乱久久久久久| 亚洲av在线观看美女高潮| freevideosex欧美| 国产老妇伦熟女老妇高清| 国产黄片美女视频| 免费观看性生交大片5| 欧美潮喷喷水| 五月玫瑰六月丁香| 高清在线视频一区二区三区| 在线观看一区二区三区| 美女视频免费永久观看网站| 国产精品国产三级专区第一集| 美女cb高潮喷水在线观看| 内射极品少妇av片p| 九色成人免费人妻av| av线在线观看网站| 少妇丰满av| 午夜激情福利司机影院| 亚洲av欧美aⅴ国产| 国产成人a∨麻豆精品| 午夜免费男女啪啪视频观看| 久久久国产一区二区| 久久久久久久国产电影| 亚洲精品日韩av片在线观看| 天堂中文最新版在线下载 | 国产亚洲av片在线观看秒播厂| 欧美日韩在线观看h| 下体分泌物呈黄色| 我的老师免费观看完整版| 嫩草影院入口| 九草在线视频观看| 嫩草影院新地址| 国产欧美日韩精品一区二区| 国产精品国产三级国产av玫瑰| 久久久久久久大尺度免费视频| 日韩av不卡免费在线播放| 国语对白做爰xxxⅹ性视频网站| 在线精品无人区一区二区三 | 久久国内精品自在自线图片| 80岁老熟妇乱子伦牲交| 国产av码专区亚洲av| 新久久久久国产一级毛片| 久久精品人妻少妇| 国产精品精品国产色婷婷| 午夜爱爱视频在线播放| 国产一区二区三区综合在线观看 | 亚洲欧美日韩东京热| 久久久久久久久大av| 下体分泌物呈黄色| 久久这里有精品视频免费| 国产午夜精品久久久久久一区二区三区| 午夜福利视频精品| 国产日韩欧美亚洲二区| 高清午夜精品一区二区三区| 一边亲一边摸免费视频| 水蜜桃什么品种好| 国产毛片在线视频| 久久精品久久久久久噜噜老黄| 久久久久九九精品影院| 久久这里有精品视频免费| 深夜a级毛片| 日韩免费高清中文字幕av| 大香蕉久久网| 久久精品久久精品一区二区三区| 高清视频免费观看一区二区| 免费大片18禁| 国产伦理片在线播放av一区| 99九九线精品视频在线观看视频| 欧美日本视频| 久久精品综合一区二区三区| 成年女人在线观看亚洲视频 | 久热这里只有精品99| 国产一区亚洲一区在线观看| 亚洲av成人精品一区久久| 久久久成人免费电影| 美女视频免费永久观看网站| 哪个播放器可以免费观看大片| 麻豆乱淫一区二区| 国产亚洲5aaaaa淫片| 国产淫片久久久久久久久| 大陆偷拍与自拍| 国产成人免费无遮挡视频| 亚洲在久久综合| 亚洲av电影在线观看一区二区三区 | 日韩欧美精品v在线| 水蜜桃什么品种好| 久久精品人妻少妇| 亚洲欧美日韩另类电影网站 | 3wmmmm亚洲av在线观看| 国产欧美日韩一区二区三区在线 | av黄色大香蕉| av免费在线看不卡| 久久久a久久爽久久v久久| 亚洲综合色惰| 老师上课跳d突然被开到最大视频| 国内揄拍国产精品人妻在线| 日本三级黄在线观看| 久久精品久久久久久噜噜老黄| 久热久热在线精品观看| 26uuu在线亚洲综合色| 久久精品国产亚洲网站| 国产高清不卡午夜福利| 亚洲av.av天堂| 美女主播在线视频| 美女内射精品一级片tv| 涩涩av久久男人的天堂| 国产av码专区亚洲av| 听说在线观看完整版免费高清| 久久这里有精品视频免费| 午夜激情久久久久久久| 日韩欧美一区视频在线观看 | 天天躁夜夜躁狠狠久久av| 亚洲国产色片| 男女边摸边吃奶| 丝袜脚勾引网站| 国产高潮美女av| 成年人午夜在线观看视频| 亚洲欧美一区二区三区黑人 | 亚洲婷婷狠狠爱综合网| 欧美区成人在线视频| 六月丁香七月| 九九在线视频观看精品| 一本一本综合久久| 欧美精品一区二区大全| 乱系列少妇在线播放| 91久久精品电影网| 少妇人妻精品综合一区二区| 国产精品偷伦视频观看了| 麻豆乱淫一区二区| 久久国产乱子免费精品| 欧美极品一区二区三区四区| 久久99热这里只有精品18| 99久久人妻综合| 菩萨蛮人人尽说江南好唐韦庄| 亚洲人成网站在线播| 男女下面进入的视频免费午夜| 午夜精品国产一区二区电影 | 日韩免费高清中文字幕av| 又黄又爽又刺激的免费视频.| 黄色欧美视频在线观看| 亚洲国产av新网站| 中文字幕免费在线视频6| 永久网站在线| 午夜亚洲福利在线播放| 久久久久久久亚洲中文字幕| 国产一区二区在线观看日韩| 51国产日韩欧美| 欧美精品一区二区大全| 久久国内精品自在自线图片| 白带黄色成豆腐渣| 亚洲欧美日韩卡通动漫| 亚洲性久久影院| 在线观看国产h片| 国产午夜精品一二区理论片| 国产欧美另类精品又又久久亚洲欧美| 亚洲美女搞黄在线观看| 国产 一区精品| 三级国产精品片| 在线观看av片永久免费下载| 久久国内精品自在自线图片| 伊人久久精品亚洲午夜| 一级毛片久久久久久久久女| 一级爰片在线观看| 日韩精品有码人妻一区| 一区二区三区乱码不卡18| 黄色一级大片看看| 自拍欧美九色日韩亚洲蝌蚪91 | 在线观看av片永久免费下载| 日日摸夜夜添夜夜爱| 大又大粗又爽又黄少妇毛片口| 亚洲怡红院男人天堂| 国产69精品久久久久777片| 我的老师免费观看完整版| 久久韩国三级中文字幕| 特级一级黄色大片| 男女国产视频网站| 欧美+日韩+精品| 99久久精品热视频| 真实男女啪啪啪动态图| 国产一区二区亚洲精品在线观看| 伦精品一区二区三区| 久久6这里有精品| 免费看光身美女| 日日摸夜夜添夜夜爱| av天堂中文字幕网| 久久久久网色| 亚洲精品,欧美精品| 国产精品久久久久久久久免| 九九久久精品国产亚洲av麻豆| 97在线人人人人妻| 啦啦啦在线观看免费高清www| a级一级毛片免费在线观看| 久久影院123| 亚洲欧美日韩东京热| 亚洲高清免费不卡视频| 免费观看在线日韩| 久久精品夜色国产| 五月伊人婷婷丁香| 丝袜喷水一区| 成人特级av手机在线观看| 色5月婷婷丁香| 国产av国产精品国产| 国产高清不卡午夜福利| 简卡轻食公司| 亚洲熟女精品中文字幕| videos熟女内射| 国产淫片久久久久久久久| 亚洲精品aⅴ在线观看| 少妇熟女欧美另类| 日韩精品有码人妻一区| 国产免费又黄又爽又色| 我要看日韩黄色一级片| 伦理电影大哥的女人| 久久久久九九精品影院| 亚洲一区二区三区欧美精品 | 18禁裸乳无遮挡动漫免费视频 | 欧美激情在线99| 免费av不卡在线播放| 色婷婷久久久亚洲欧美| 国产精品偷伦视频观看了| 国产爽快片一区二区三区| 一本—道久久a久久精品蜜桃钙片 精品乱码久久久久久99久播 | 精品国产乱码久久久久久小说| 尤物成人国产欧美一区二区三区| 一级片'在线观看视频| 中文资源天堂在线| 大片免费播放器 马上看| 永久网站在线| 色视频www国产| 亚洲av成人精品一二三区| 欧美一级a爱片免费观看看| 狂野欧美激情性xxxx在线观看| 国产成人freesex在线| 最后的刺客免费高清国语| 交换朋友夫妻互换小说| 香蕉精品网在线| 一级毛片 在线播放| 人妻 亚洲 视频| 69人妻影院| 高清午夜精品一区二区三区| 99热全是精品| 精品一区二区免费观看| 国产在线男女| 国精品久久久久久国模美| 夜夜看夜夜爽夜夜摸| 亚洲欧美日韩另类电影网站 | 九色成人免费人妻av| 高清视频免费观看一区二区| 最后的刺客免费高清国语| 天天躁夜夜躁狠狠久久av| 欧美成人精品欧美一级黄| 国产精品国产三级国产av玫瑰| 小蜜桃在线观看免费完整版高清| 国产免费视频播放在线视频| 欧美精品一区二区大全| 中文字幕亚洲精品专区| 久久久国产一区二区| 欧美激情国产日韩精品一区| 婷婷色综合www| 汤姆久久久久久久影院中文字幕| 直男gayav资源| 一级片'在线观看视频| 亚洲在久久综合| 美女cb高潮喷水在线观看| 午夜福利在线在线| av线在线观看网站| 激情五月婷婷亚洲| 久久久久久久久久久免费av| 亚洲成人久久爱视频| 久久精品国产自在天天线| 成人综合一区亚洲| 久久精品久久精品一区二区三区| 亚洲av日韩在线播放| 免费看av在线观看网站| 国产欧美另类精品又又久久亚洲欧美| 国产成人一区二区在线| 免费观看的影片在线观看| 日本与韩国留学比较| 亚洲天堂国产精品一区在线| 在线观看国产h片| 在线a可以看的网站| 亚洲欧美日韩另类电影网站 | 久久久精品94久久精品| 亚洲经典国产精华液单| 51国产日韩欧美| 亚洲av成人精品一二三区| a级一级毛片免费在线观看| 丝袜美腿在线中文| 天天躁夜夜躁狠狠久久av| 内地一区二区视频在线| 丰满乱子伦码专区| av福利片在线观看| 日韩一区二区三区影片| 我的女老师完整版在线观看| 中文资源天堂在线| 国产亚洲av嫩草精品影院| av播播在线观看一区| av在线老鸭窝| 久久久久久久精品精品| 国产女主播在线喷水免费视频网站| 国产精品蜜桃在线观看| 日本猛色少妇xxxxx猛交久久| 大香蕉久久网| 最近2019中文字幕mv第一页| 少妇人妻 视频| 小蜜桃在线观看免费完整版高清| 国产免费一区二区三区四区乱码| 国产极品天堂在线| 乱码一卡2卡4卡精品| 久久综合国产亚洲精品| 日本-黄色视频高清免费观看| 亚洲精品中文字幕在线视频 | 日韩av不卡免费在线播放| 国产久久久一区二区三区| 国产视频首页在线观看| 全区人妻精品视频| 又爽又黄无遮挡网站| 亚洲自偷自拍三级| 亚洲欧美日韩东京热| 国产人妻一区二区三区在| 亚洲欧美成人综合另类久久久| 国产午夜精品一二区理论片| 一区二区三区乱码不卡18| 亚洲成人精品中文字幕电影| 精品一区二区三卡| av在线播放精品| 精品人妻偷拍中文字幕| 国产乱人偷精品视频| 偷拍熟女少妇极品色| 日本黄大片高清| 欧美极品一区二区三区四区| av一本久久久久| 在线天堂最新版资源| 嘟嘟电影网在线观看| 国产精品三级大全| 在线a可以看的网站| 啦啦啦啦在线视频资源| 啦啦啦在线观看免费高清www| 五月开心婷婷网| 99热全是精品| 亚洲精品色激情综合| 精品酒店卫生间| 韩国av在线不卡| 黑人高潮一二区| 亚洲最大成人中文| 日韩 亚洲 欧美在线| 菩萨蛮人人尽说江南好唐韦庄| 色5月婷婷丁香| 男人爽女人下面视频在线观看| 草草在线视频免费看| 中文资源天堂在线| 国产老妇伦熟女老妇高清| 国产精品秋霞免费鲁丝片| 亚洲无线观看免费| 久久久久网色| 亚洲伊人久久精品综合| 99re6热这里在线精品视频| 国产 精品1| av免费在线看不卡| 久久久久国产网址| 秋霞在线观看毛片| 精品午夜福利在线看| 一个人看的www免费观看视频| 国产亚洲一区二区精品| 亚洲av福利一区| 国产精品一区二区在线观看99| 韩国高清视频一区二区三区| 2018国产大陆天天弄谢| 欧美日韩综合久久久久久| 国产精品久久久久久久久免| 激情五月婷婷亚洲| 制服丝袜香蕉在线| 亚洲图色成人| 日韩成人av中文字幕在线观看| 午夜福利高清视频| 国产亚洲av片在线观看秒播厂| 亚洲欧美日韩卡通动漫|