鸡白痢沙门氏菌感染对雏鸡脾脏miRNA表达谱的影响
Effects of Salmonella Pullorum Infection on miRNA Expression Profiles in Spleen of Chicks
-
Keywords:
- Salmonella Pullorum /
- SPF chicks /
- spleen /
- differentially expressed miRNA /
- immunity
-
鸡白痢沙门氏菌(Salmonella Pullorum)是禽类重要的病原体,可引起鸡白痢,死亡率高,对家禽业构成严重威胁[1-2]。鸡白痢在世界许多地区(中东、非洲、亚洲和中南美洲)流行,造成巨大的经济损失[3-4]。该病多发于雏鸡,其典型症状是排白色粪便,肛周绒毛粘有粪便。目前,暂无有效疫苗预防鸡白痢沙门氏菌,常用净化手段是淘汰阳性鸡。因此,探究鸡白痢沙门氏菌感染宿主的免疫机制对防治鸡白痢尤为重要。
microRNA (miRNA) 是在动物、植物和一些病毒中发现的一类非编码小RNA (长18~25 nt),可调控基因表达水平[5],在细胞增殖、分化和免疫等生物学过程中发挥控制作用[6]。miR-146a、miR-155和Let-7等miRNAs参与小鼠对沙门氏菌或脂多糖感染的免疫应答调节[7];gga-miR-125b-5p、gga-miR-34a-5p、gga-miR-1416-5p和gga-miR-1662可通过调节其靶基因在肠炎沙门氏菌感染中发挥重要作用[8];gga-miR-1306-5p通过抑制Toll相互作用蛋白调节肠炎沙门氏菌感染宿主的免疫反应,刺激炎症细胞因子的产生[9]。尽管与细菌感染相关的miRNA研究已取得不少成果,但关于miRNA在鸡白痢沙门氏菌感染宿主反应中的作用信息仍然有限。脾脏作为鸡体内最大的外周淋巴器官,对禽类免疫系统的功能比其他脊椎动物更重要[10]。鸡脾脏中有来自法氏囊的B细胞、来自胸腺的T细胞以及大量巨噬细胞,是重要的免疫器官;此外,脾脏还是鸡白痢沙门氏菌定植的靶器官。本研究采用Illumina Hiseq2500高通量测序技术和生物信息学分析方法,在全基因组范围内筛选鸡白痢沙门氏菌感染后SPF雏鸡脾脏差异表达miRNA,获取与感染相关的潜在miRNA,并进一步分析其可能参与调控的通路,以期为阐明miRNA在鸡白痢沙门氏菌感染过程的作用机制提供理论基础。
1. 材料与方法
1.1 试验菌株和动物
供试鸡白痢沙门氏菌CVCC1792购自中国兽医微生物菌种保藏管理中心;80只3日龄SPF雏鸡购自中国农业科学院哈尔滨兽医研究所国家禽类实验动物资源库。
1.2 试验方法
1.2.1 动物感染和样品采集
随机选取3日龄SPF雏鸡,感染组50只,对照组30只。以CVCC1792半数致死量饲喂感染组,即以含3.263×108 CFU的菌液0.2 mL饲喂感染组,对照组饲喂无菌 PBS 0.2 mL,2组隔离饲养。攻菌24 h后处死雏鸡,取感染组和对照组脾脏于−80 ℃保存,用于RNA提取。
1.2.2 miRNA文库构建及测序
感染组和对照组各选取9只雏鸡的脾脏,每3个脾脏混为1个样,采用Trizol提取RNA,并用1%琼脂糖凝胶电泳对RNA进行质检。按Small RNA样本制备试剂盒(TruSeq™ Small RNA Sample Preparation Small RNA Sample Prep Kits,illumina公司生产)说明书构建6个miRNA文库,即3个感染组文库(T1、T2和T3)及3个对照组文库(C1、C2和C3)。文库质检合格后送至联川生物,利用Illumina Hiseq 2500进行测序。
1.2.3 miRNA的鉴定与分析
运用联川生物开发的 ACGT101-miR软件 (LC Sciences,Houston,Texas,USA)去除高通量测序原始数据中带有接头污染、无法确定以及碱基比例大于10%的序列;再进行长度筛选,保留碱基长度为18~26 nt的序列,获得过滤后序列;将过滤后序列与各RNA数据库序列进行比对,如mRNA数据库、RFam数据库(包含rRNA、tRNA、snRNA和snoRNA等)以及Repbase数据库(重复序列数据库),过滤去除非miRNA序列,获得有效序列;将有效序列映射到鸡基因组和miRBase (v22.0)鸡前体miRNA上,获得已知miRNA。对于可以比对到鸡基因组、但比对不到鸡前体miRNA的序列,采用mfold软件进行二级结构预测,满足所有发夹结构所需原则,则获得新预测miRNA;最后采用Ensembl和miRBase数据库对已知和新预测miRNA进行序列注释。
1.2.4 差异表达miRNA的筛选及聚类分析
将鉴定到的已知miRNA与新预测的miRNA序列进行TPM (transcript per million)表达量统计,即每百万条序列的转录本,并对其进行表达量归一化处理[11],得到norm值。由于所测样本具有生物学重复,采用基于正态分布的方法计算P值;感染组和对照组差异表达的miRNA分析使用t检验,并以P<0.05为阈值。将各组差异miRNA集用于层次聚类分析,使用pheatmap绘制聚类热图。
1.2.5 差异表达miRNA靶基因的预测以及GO和KEGG富集分析
运用TargetScan (v5.0)和miRanda (v3.3a)进行差异表达的miRNA靶基因预测。筛选阈值为TargetScan score≥50且miRanda Energy<−10,获得预测靶基因; 使用 Gene Ontology 数据库(http://geneontology.org)和KEGG数据库(http://geneontology.org)进行GO功能富集分析和KEGG信号通路分析。
1.2.6 高通量测序结果的RT-qPCR验证
选取差异表达miRNA进行荧光定量验证,其引物序列见表1。按照miRNA 1st Strand cDNA Synthesis Kit (by stem-loop) 试剂盒说明书反转录总RNA,按照miRNA Universal SYBR qPCR Master Mix试剂盒说明书对反转录的cDNA在Bio-Rad CFX96上进行荧光定量,定量体系为20 µL,反应程序为:95 ℃ 5 min,1个循环;95 ℃ 10 s,60 ℃ 60 s,40个循环,使用仪器默认熔解程序采集熔解曲线。每个样品3个重复,以5S rRNA 为内参,用2−ΔΔCt 法计算miRNA的相对表达量。
表 1 miRNA RT-qPCR引物Table 1. Primers of miRNA RT-qPCRmiRNA 引物名称
primer name引物序列 (5′→3′)
primer sequencegga-let-7j-5p RTj5p GTCGTATCCAGTGCAGGGTCCGAGGTATTCGCACTGGATACGACAACTAT j5p-F GCGCGTGAGGTAGTAGGTTGT R1 AGTGCAGGGTCCGAGGTATT gga-miR-125b-3p RT125 GTCGTATCCAGTGCAGGGTCCGAGGTATTCGCACTGGATACGACAGGTCCCA 125-F CCGCGACAAGTCAGGCTC R1 AGTGCAGGGTCCGAGGTATT gga-miR-130b-3p RT130 GTCGTATCCAGTGCAGGGTCCGAGGTATTCGCACTGGATACGACACGCCC 130-F CGCGCAGTGCAATAATGAAA R1 AGTGCAGGGTCCGAGGTATT gga-let-7b RT7b GTCGTATCCAGTGCAGGGTCCGAGGTATTCGCACTGGATACGACAACCAC 7b-F GCGCGTGAGGTAGTAGGTTGT R1 AGTGCAGGGTCCGAGGTATT gga-miR-451 RT451 GTCGTATCCAGTGCAGGGTCCGAGGTATTCGCACTGGATACGACAAACTC 451-F CGCGAAACCGTTACCATTACT R1 AGTGCAGGGTCCGAGGTATT gga-miR-19b-3p RT19b GTCGTATCCAGTGCAGGGTCCGAGGTATTCGCACTGGATACGACTCAGTT 19b-F CAGCGTGTGTGCAAATCCAT R1 AGTGCAGGGTCCGAGGTATT gga-miR-153-3p RT153 GTCGTATCCAGTGCAGGGTCCGAGGTATTCGCACTGGATACGACTCACTT 153-F CCAGCGTGTTGCATAGTCAC R1 AGTGCAGGGTCCGAGGTATT gga-miR-21-3p RT21 CTCAACTGGTGTCGTGGAGTCGGCAATTCAGTTGAGGACAGC 21-F CGCGCAACAACAGTCGGTAG R2 ACTGGTGTCGTGGAGTCGGC gga-miR-147 RT147 CTCAACTGGTGTCGTGGAGTCGGCAATTCAGTTGAGGCAGAA 147-F GCGGTGTGCGGAAATGC R2 ACTGGTGTCGTGGAGTCGGC 5S rRNA 5S-F CCATACCACCCTGGAAACGC 5S-R TACTAACCGAGCCCGACCCT 1.2.7 数据统计与分析
数据统计分析采用R (v3.0.1)软件和GraphPad Prism 6进行;采用SPSS 26进行独立样本t检验,并进行显著性分析。
2. 结果与分析
2.1 测序数据质控及有效数据统计
RNA质检合格后构建感染组和对照组文库并上机测序;为保证数据的质量及可靠性,对原始数据进行质控,得到的有效数据见表2。
表 2 有效数据统计结果Table 2. Statistical results of valid data样品
samples原始序列
raw reads筛选序列
selected reads低质量序列
junk reads过滤后序列
clean readsRfam数据库序列
reads of Rfam database信使RNA
mRNA重复序列
repeats reads有效序列
valid reads有效序列占比/%
rate of valid readsC1 9380231 4988969 9889 4381373 844861 229212 321743 3373552 35.96 C2 11884024 2433385 12240 9438399 852781 158504 387259 8483589 71.39 C3 12181462 2560633 12531 9608298 962446 177670 390585 8525520 69.99 T1 22529885 12411658 10214 10108013 4452492 235137 2844813 5573614 24.74 T2 11378691 5185418 6500 6186773 2121851 148841 1354775 4006009 35.21 T3 11433648 5298273 23510 6111865 1255664 140035 575972 4790964 41.90 2.2 miRNA的鉴定和预测结果
鉴定获得817个miRNAs,包括725个已知miRNAs和92个新miRNAs,且miRNA在22 nt处的数量最多(图1),与Dicer酶切产物的长度分布基本一致,表明测序结果可靠。
2.3 差异表达miRNA的筛选与聚类
感染组和对照组共筛选到29个显著差异表达的miRNAs,其中24个为已知miRNAs。与对照组相比,感染组中有15个miRNAs的表达量显著上调,14个miRNAs的表达量显著下调。前10个差异显著的miRNAs为gga-miR-1a-3p、gga-miR-214、cgr-miR-27a-3p、gga-let-7j-5p、gga-miR-100-5p、gga-miR-125b-3p、gga-miR-130b-3p、gga-miR-451、gga-let-7d和gga-miR-147 (图2)。
由图3可知:对照组(C组)和感染组(T组)的样品聚成2簇,说明2组的miRNA表达谱不同。C组中,C2和C3聚为一类,说明两者的miRNA表达模式较相似;T组中,T2和T3聚为一类,说明两者的miRNA表达模式较相似。此外,差异表达的miRNA被分成2簇,第1簇miRNA在对照组中表达较高,包括gga-let-7b、gga-miR-18b-3p和gga-miR-100-5p等;第2簇在感染组中表达量较高,包括gga-miR-451、gga-miR-193-3p和gga-miR-147等。
2.4 差异表达miRNA的GO功能富集和KEGG通路富集分析
差异表达miRNA的靶基因共预测到10659个miRNA-mRNA互作靶点,其中gga-miR-214和gga-miR-21-3p预测靶基因最多,分别为802和606个。靶基因GO功能富集显示(图4):T细胞共刺激、凋亡过程负调控、T细胞活化、脂多糖介导的信号通路正调控和NIK/NF-κB信号的负调控等免疫相关生物学过程显著富集;显著富集的细胞组分主要有胞浆、膜的组成成分和内体等;显著富集的分子功能主要有GTP结合、GTP酶活性和跨膜转运蛋白结合等。
KEGG富集显示(图5):显著富集的免疫相关通路有FoxO信号通路、铁死亡、p53信号通路、凋亡、溶酶体、内吞作用和自噬—动物等。
2.5 差异表达miRNA的RT-qPCR验证结果
由图6可知:gga-miR-19b-3p、gga-miR-21-3p、gga-miR-147、gga-miR-153-3p和gga-miR-451的差异显著上调,gga-let-7b、gga-let7j-5p、gga-miR-125b-3p 和 gga-miR-130b-3p的差异显著下调。该结果与高通量测序结果一致,且差异表达倍数相近,说明高通量测序结果可靠且准确。
3. 讨论
本研究利用Illumina高通量测序平台检测鸡白痢沙门氏菌感染雏鸡脾脏中miRNA的表达谱,共获得817个miRNAs,包括725个已知miRNAs和92个新miRNAs。其中29个miRNAs的表达存在显著差异,是鸡白痢沙门氏菌感染潜在的候选miRNA;选取9个差异表达的miRNAs,利用RT-qPCR技术进行验证,结果显示9个差异表达miRNAs表达量趋势与高通量测序结果一致,说明Illumina测序数据具有覆盖率高、通量高和测序准确可靠等特性。
Toll样受体具有识别病原菌和激活固有免疫的作用。TLR4及其信号分子MyD88和 NF-кB介导革兰氏阴性菌脂多糖的刺激传导,是连接天然免疫和获得性免疫的桥梁[12]。CHEN等[13]研究显示:gga-miR-19b-3p过表达促进了新城疫病毒诱导的炎性细胞因子产生并抑制该病毒复制,通过靶向RNF11和ZMYND11激活NF-κB信号传导,增强炎性细胞因子产生。本研究中,感染鸡白痢沙门氏菌后脾脏中gga-miR-19b-3p显著上调,该miRNA是否也通过NF-κB通路抵抗鸡白痢沙门氏菌的感染还需后续研究证实。
本研究中差异表达的gga-miR-21-3p、gga-miR-130b-3p、 gga-miR-451 和 gga-miR-147 等也可参与调控 MAPK 和 NF-κB 等炎症相关通路,抵御细菌感染。在滑液囊支原体感染DF-1细胞模型中,gga-miR-21上调,靶向MAP3K1激活MAPK和NF-κB通路,促进TNF-α和IL-6等产生,并通过增加细胞周期和抑制细胞凋亡促进炎症细胞因子的产生和细胞增殖,以抵御支原体感染[14];在滑液囊支原体感染时,miR-130b-3p上调,激活PI3K/AKT/NF-κB途径,下调PTEN促进细胞增殖[15];在脂多糖处理的小鼠中,miR-130b-3p上调并靶向IRF1抑制经典活化巨噬细胞极化[16];gga-miR-451过表达并靶向YWHAZ降低了滑液囊支原体感染诱导的炎性细胞因子(如TNF-α、IL-1β和IL-6等)产生[17];miR-147已被证明与人类传染病、癌症、心血管疾病和神经退行性疾病等相关[18],在感染分枝杆菌的小鼠巨噬细胞中,miR-147-3p表达上调,抑制了IL-6和IL-10产生,并显著降低了分枝杆菌的存活率[19];核苷酸结合寡聚化结构域蛋白的激活可抑制小鼠肺内皮细胞中的miR-147-3p表达,导致TNF-α和IL-6等细胞因子上调,而小鼠TNF-α和IL-6 mRNA的3′UTR含有miR-147的预测结合位点[20],miR-147-3p表现为负调控炎症因子的表达。在GO富集分析中,脂多糖介导的信号通路正调控以及NIK/NF-κB信号的负调控等生物学过程也显著富集,说明gga-miR-19b-3p、gga-miR-21-3p、gga-miR-130b-3p、gga-miR-451和gga-miR-147可能参与调控MAPK和NF-κB等炎症相关通路,进而调控鸡白痢沙门氏菌的感染。
溶酶体是宿主防御系统的组成部分,含有大量酸性水解酶以抵抗细菌感染。通过与含有细菌的内体融合,溶酶体可经水解酶系统降解细菌[21]。然而,沙门氏菌已进化出逃避溶酶体降解的系统。沙门氏菌隐藏于含沙门氏菌液泡(Salmonella-containing vacuole,SCV)的特殊隔室中,从而阻止与溶酶体融合[22-23]。沙门氏菌的感染可减少宿主细胞溶酶体数量,而SCV不断分裂形成新的SCV,导致没有足够的溶酶体靶向不断形成的SCV,从而有助于沙门氏菌的胞内存活[17, 24]。鸡白痢沙门氏菌的持续感染机制主要基于巨噬细胞中的存活、复制和SCV分裂[25]。本研究中溶酶体通路在KEGG分析中显著富集,该通路预测到靶基因最多的前3个差异表达miRNAs为gga-miR-214、gga-miR-21-3p和gga-miR-19b-3p,它们可能是潜在参与调控溶酶体通路的miRNA。KEGG分析显示:差异表达miRNA预测靶基因显著富集于凋亡、内吞作用和自噬等重要免疫相关通路,其具体调控机制还需进一步研究证实。
4. 结论
成功鉴定了鸡白痢沙门氏菌感染雏鸡脾脏miRNA表达谱,获得29个鸡白痢沙门氏菌感染潜在相关miRNAs,它们可能通过调控NF-κB和溶酶体等免疫相关通路参与鸡白痢感染宿主引起的免疫反应。本研究为阐明miRNA在鸡白痢沙门氏菌感染中的调控机制提供了理论基础。
-
表 1 miRNA RT-qPCR引物
Table 1 Primers of miRNA RT-qPCR
miRNA 引物名称
primer name引物序列 (5′→3′)
primer sequencegga-let-7j-5p RTj5p GTCGTATCCAGTGCAGGGTCCGAGGTATTCGCACTGGATACGACAACTAT j5p-F GCGCGTGAGGTAGTAGGTTGT R1 AGTGCAGGGTCCGAGGTATT gga-miR-125b-3p RT125 GTCGTATCCAGTGCAGGGTCCGAGGTATTCGCACTGGATACGACAGGTCCCA 125-F CCGCGACAAGTCAGGCTC R1 AGTGCAGGGTCCGAGGTATT gga-miR-130b-3p RT130 GTCGTATCCAGTGCAGGGTCCGAGGTATTCGCACTGGATACGACACGCCC 130-F CGCGCAGTGCAATAATGAAA R1 AGTGCAGGGTCCGAGGTATT gga-let-7b RT7b GTCGTATCCAGTGCAGGGTCCGAGGTATTCGCACTGGATACGACAACCAC 7b-F GCGCGTGAGGTAGTAGGTTGT R1 AGTGCAGGGTCCGAGGTATT gga-miR-451 RT451 GTCGTATCCAGTGCAGGGTCCGAGGTATTCGCACTGGATACGACAAACTC 451-F CGCGAAACCGTTACCATTACT R1 AGTGCAGGGTCCGAGGTATT gga-miR-19b-3p RT19b GTCGTATCCAGTGCAGGGTCCGAGGTATTCGCACTGGATACGACTCAGTT 19b-F CAGCGTGTGTGCAAATCCAT R1 AGTGCAGGGTCCGAGGTATT gga-miR-153-3p RT153 GTCGTATCCAGTGCAGGGTCCGAGGTATTCGCACTGGATACGACTCACTT 153-F CCAGCGTGTTGCATAGTCAC R1 AGTGCAGGGTCCGAGGTATT gga-miR-21-3p RT21 CTCAACTGGTGTCGTGGAGTCGGCAATTCAGTTGAGGACAGC 21-F CGCGCAACAACAGTCGGTAG R2 ACTGGTGTCGTGGAGTCGGC gga-miR-147 RT147 CTCAACTGGTGTCGTGGAGTCGGCAATTCAGTTGAGGCAGAA 147-F GCGGTGTGCGGAAATGC R2 ACTGGTGTCGTGGAGTCGGC 5S rRNA 5S-F CCATACCACCCTGGAAACGC 5S-R TACTAACCGAGCCCGACCCT 表 2 有效数据统计结果
Table 2 Statistical results of valid data
样品
samples原始序列
raw reads筛选序列
selected reads低质量序列
junk reads过滤后序列
clean readsRfam数据库序列
reads of Rfam database信使RNA
mRNA重复序列
repeats reads有效序列
valid reads有效序列占比/%
rate of valid readsC1 9380231 4988969 9889 4381373 844861 229212 321743 3373552 35.96 C2 11884024 2433385 12240 9438399 852781 158504 387259 8483589 71.39 C3 12181462 2560633 12531 9608298 962446 177670 390585 8525520 69.99 T1 22529885 12411658 10214 10108013 4452492 235137 2844813 5573614 24.74 T2 11378691 5185418 6500 6186773 2121851 148841 1354775 4006009 35.21 T3 11433648 5298273 23510 6111865 1255664 140035 575972 4790964 41.90 -
[1] BARROW P A, FREITAS N O C. Pullorum disease and fowl typhoid-new thoughts on old diseases: a review[J]. Avian Pathology, 2011, 40(1): 1. DOI: 10.1080/03079457.2010.542575.
[2] SHIVAPRASAD H L. Fowl typhoid and pullorum disease[J]. Revue Scientifique et Technique, 2000, 19(2): 405. DOI: 10.20506/rst.19.2.1222.
[3] TADESSE S, ASHENAFI H, ASCHALEW Z. Seroprevalence study of newcastle disease in local chickens in central Ethiopia[J]. International Journal of Applied Research in Veterinary Medicine, 2005, 3(1): 25.
[4] TEFERI M, NEJASH A. Epidemiology and economic importance of Pullorum disease in poultry: a review[J]. Global Veterinaria, 2016, 17(3): 228. DOI: 10.5829/idosi.gv.2016.17.03.103123.
[5] BARTEL D P. MicroRNAs: genomics, biogenesis, mechanism, and function[J]. Cell, 2004, 116(2): 281. DOI: 10.1016/s0092-8674(04)00045-5.
[6] HE L, HANNON G J. MicroRNAs: small RNAs with a big role in gene regulation[J]. Nature Reviews Genetics, 2004, 5(7): 522. DOI: 10.1038/nrg1379.
[7] SCHNITGER A K, MACHOVA A, MUELLER R U, et al. Listeria monocytogenes infection in macrophages induces vacuolar-dependent host miRNA response[J]. PLoS One, 2011, 6(11): e27435. DOI: 10.1371/journal.pone.0027435.
[8] WU G X, QI Y K, LIU X Y, et al. Cecal microRNAome response to Salmonella enterica serovar enteritidis infection in white leghorn layer[J]. BMC Genomics, 2017, 18(1): 77. DOI: 10.1186/s12864-016-3413-8.
[9] SUN W W, LIU R R, LI P, et al. Chicken gga-miR-1306-5p targets Tollip and plays an important role in host response against Salmonella enteritidis infection[J]. Journal of Animal Science and Biotechnology, 2019, 10: 59. DOI: 10.1186/s40104-019-0365-2.
[10] SMITH K G, HUNT J L. On the use of spleen mass as a measure of avian immune system strength[J]. Oecologia, 2004, 138(1): 28. DOI: 10.1007/s00442-003-1409-y.
[11] LI X, SHAHID M Q, WU J W, et al. Comparative small RNA analysis of pollen development in autotetraploid and diploid rice[J]. International Journal of Molecular Sciences, 2016, 17(4): 499. DOI: 10.3390/ijms17040499.
[12] AKIRA S. Toll receptor families: structure and function[J]. Seminars in Immunology, 2004, 16: 1. DOI: 10.1016/j.smim.2003.10.001.
[13] CHEN Y, LIU W, XU H X, et al. gga-miR-19b-3p inhibits newcastle disease virus replication by suppressing inflammatory response via targeting RNF11 and ZMYND11[J]. Frontiers in Microbiology, 2019, 10: 2006. DOI: 10.3389/fmicb.2019.02006.
[14] ZHAO Y B, ZOU M Y, SUN Y F, et al. gga-miR-21 modulates Mycoplasma gallisepticum (HS strain)-induced inflammation via targeting MAP3K1 and activating MAPKs and NF-κB pathways[J]. Veterinary Microbiology, 2019, 237: 108407. DOI: 10.1016/j.vetmic.2019.108407.
[15] YUAN B, ZOU M Y, ZHAO Y B, et al. Up-regulation of miR-130b-3p activates the PTEN/PI3K/AKT/NF-κB pathway to defense against Mycoplasma gallisepticum (HS strain) infection of chicken[J]. International Journal of Molecular Sciences, 2018, 19(8): 2172. DOI: 10.3390/ijms19082172.
[16] GUO Q, ZHU X X, WEI R, et al. miR-130b-3p regulates M1 macrophage polarization via targeting IRF1[J]. Journal of Cellular Physiology, 2021, 236(3): 2008. DOI: 10.1002/jcp.29987.
[17] ZHAO Y B, ZHANG K, ZOU M Y, et al. gga-miR-451 negatively regulates Mycoplasma gallisepticum (HS strain)-induced inflammatory cytokine production via targeting YWHAZ[J]. International Journal of Molecular Sciences, 2018, 19(4): 1191. DOI: 10.3390/ijms19041191.
[18] LIN L, HU K B. MiR-147: functions and implications in inflammation and diseases[J]. Microrna, 2021, 10(2): 91. DOI: 10.2174/2211536610666210707113605.
[19] ZUO X S, WANG L, BAO Y Q, et al. The ESX-1 virulence factors downregulate miR-147-3p in Mycobacterium marinum-infected macrophages[J]. Infection and Immunity, 2020, 88(6): e00088. DOI: 10.1128/IAI.00088-20.
[20] VLACIL A K, VOLLMEISTER E, BERTRAMS W, et al. Identification of microRNAs involved in NOD-dependent induction of pro-inflammatory genes in pulmonary endothelial cells[J]. PLoS One, 2020, 15(4): e0228764. DOI: 10.1371/journal.pone.0228764.
[21] ESWARAPPA S M, NEGI V D, CHAKRABORTY S, et al. Division of the Salmonella-containing vacuole and depletion of acidic lysosomes in Salmonella-infected host cells are novel strategies of Salmonella enterica to avoid lysosomes[J]. Infection and Immunity, 2010, 78(1): 68. DOI: 10.1128/IAI.00668-09.
[22] HASHIM S, MUKHERJEE K, RAJE M, et al. Live Salmonella modulate expression of Rab proteins to persist in a specialized compartment and escape transport to lysosomes[J]. The Journal of Biological Chemistry, 2000, 275(21): 16281. DOI: 10.1074/jbc.275.21.16281.
[23] GARVIS S G, BEUZÓN C R, HOLDEN D W. A role for the PhoP/Q regulon in inhibition of fusion between lysosomes and Salmonella-containing vacuoles in macrophages[J]. Cellular Microbiology, 2001, 3(11): 731. DOI: 10.1046/j.1462-5822.2001.00153.x.
[24] LAHIRI A, ESWARAPPA S M, DAS P, et al. Altering the balance between pathogen containing vacuoles and lysosomes: a lesson from Salmonella[J]. Virulence, 2010, 1(4): 325. DOI: 10.4161/viru.1.4.12361.
[25] LI Q C, WANG X, XIA J, et al. Salmonella-containing vacuole development in avian cells and characteristic of cigR in Salmonella enterica serovar Pullorum replication within macrophages[J]. Veterinary Microbiology, 2018, 223: 65. DOI: 10.1016/j.vetmic.2018.07.013.