Transcriptome analysis of Drosophila S2 cells infected by Listeria innocua
-
摘要:目的
英诺克李斯特菌Listeria innocua为李斯特菌属Listeria的一种非致病菌，与致病菌单增李斯特菌L. monocytogenes具有共同的毒力祖先，本文旨在研究英诺克李斯特菌侵染后宿主基因表达的变化情况，为宿主调控、防御单增李斯特菌的危害提供参考。
方法用英诺克李斯特菌侵染果蝇S2细胞，利用转录组测序检测该细菌侵染后果蝇S2细胞的基因表达变化，采用qPCR验证英诺克李斯特菌侵染后果蝇S2细胞中差异表达的基因。
结果英诺克李斯特侵染3 h后，宿主果蝇S2细胞中Toll/Imd信号通路发生显著上调，吞噬体和霍乱弧菌Vibrio cholerae感染信号通路显著下调；抗菌肽基因DmDef(DmDefensin)、DmDro(DmDrosomycin)、DmDpt A(DmDptericin A)、DmDpt B(DmDptericin B)、DmMtk(DmMetchnikowin)、DmCec A2(DmCecropin A2)、DmAtt A(DmAttacin A)、DmAtt B(DmAttacin B)、DmAtt D(DmAttacin D)、DmCec B(DmCecropin B)均被显著地诱导表达，其中，上调幅度最大的基因是DmDef，上调了9.805倍。qPCR验证结果显示，DmMtk、DmAtt A、DmDro和DmDef基因表达分别上调了8.180、7.533、7.204和4.569倍。
结论英诺克李斯特菌侵染果蝇S2细胞后，变化最显著的基因类群为抗菌肽。本研究全面调查了英诺克李斯特菌侵染后宿主基因表达的变化情况，为揭示非致病菌侵染后宿主细胞的反应、细菌与宿主的互作规律提供了参考。
Abstract:ObjectiveListeria innocua is a non-pathogenic bacterium from the Listeria genus, which harbors the virulence factors evolved from the same ancestor with the pathogenic bacterium L. monocytogenes. This study aims to investigate the transcriptional variations of host cells after L. innocuainfection
, and provide a basis for host regulation and prevention of damage from L. monocytogenes. MethodWe used L. innocua to infect Drosophila melanogaster S2 cells and analyzed the change of gene expression in Drosophila S2 cells by transcriptome sequencing. The differentially expressed genes in Drosophila S2 cells infected by L. innocua was verified by qPCR.
ResultThe Toll/Imd signaling pathways were significantly upregulated in Drosophila S2 cells after L. innocua infection for three hours, while the phagosome and Vibrio cholerae infection signaling pathways were significantly downregulated. The antimicrobial peptide genes including DmDef (DmDefensin), DmDro (DmDrosomycin), DmDpt A (DmDptericin A), DmDpt B (DmDptericin B), DmMtk (DmMetchnikowin), DmCec A2 (DmCecropin A2), DmAtt A (DmAttacin A), DmAtt B (DmAttacin B), DmAtt D (DmAttacin D), and DmCec B (DmCecropin B) were significantly induced after L. innocua infection. Among them, the most upregulated gene was DmDef with 9.805 fold change. The qPCR verification results showed that the expressions of DmMtk, DmAtt A, DmDro and DmDef genes were upregulated by 8.180, 7.533, 7.204 and 4.569 fold.
ConclusionAfter L. innocua infection, the genes with the most significant change in Drosophila S2 cells are antimicrobial peptide genes. This study offers a comprehensive investigation of gene expression changes in Drosophila S2 cells after L. innocua infection, and provide a reference for revealing the response of host cells evoked by non-pathogenic bacteria as well as interaction studies between bacterial pathogens and hosts.
-
李斯特菌属Listeria是革兰阳性菌，隶属厚壁菌门Firmicutes芽孢杆菌纲Bacilli芽孢杆菌目Bacillales，该属包括17个种。狭义上的李斯特菌属包括6个种，即单增李斯特菌L. monocytogenes、绵羊李斯特菌L. ivanovii、威尔斯李斯特菌L. welshimeri、默氏李斯特菌L. murrayi、西尔李斯特菌L. seeligeri和英诺克李斯特菌L. innocua，公认的致病菌只有2种，即单增李斯特菌和绵羊李斯特菌[1-2]。其中单增李斯特菌致病性最强，可导致人与动物患李斯特菌病[3]。绵羊李斯特菌又称伊氏李斯特菌，仅能引起动物李斯特菌病[4]，但也有引起人患李斯特菌病的个别案例[2]。此外，非常罕见英诺克李斯特菌和西尔李斯特菌引起人李斯特病的报道，且二者不具有溶血活性，因此它们一般被认为是非致病菌。其中，英诺克李斯特菌常伴随单增李斯特菌出现，并被视为后者的指示菌[5]。李斯特菌属成员的中间宿主范围较广，常常生存于蔬菜、乳类和肉类制品中；适应性强，在2~42 ℃均能正常生长[6-7]。李斯特菌病的潜伏期一般小于24 h[8]，临床表现为脑膜炎、菌血症、眼内炎、颈淋巴结炎、胆囊炎、骨髓炎及关节炎等，易感人群为老人、孕妇、新生儿和免疫缺陷患者[9]。在动物中，由绵羊李斯特菌引起的李斯特菌病可表现为流产、败血症等，但不引起脑膜炎[10]。鉴于李斯特菌有危害性大、宿主谱广和适应性强的特点，很多国家都已经采取相关措施防控食品中李斯特菌的污染，一些发达国家如美国等已将李斯特菌病定为法定传染病，但目前我国仅在部分省市对李斯特菌病开展相关检测工作[11]。
目前对李斯特菌的研究多集中在单增李斯特菌上，对其他李斯特菌的研究较少。已知的单增李斯特菌毒力因子主要有溶血素基因hly编码的李斯特溶血素O(Listerriolysin O，LLO)、plcA编码的磷脂酰肌醇特异性磷脂酶C(Phospholipase C，PI-PLC)、plcB基因编码的磷脂酰胆碱磷脂酶(Phosphatidylcholine-specific phospholipase C，PC-PLC)、SmcL基因编码的神经磷脂酶C(Nerve phoshpolipase C，NPI-PLC)、肌动蛋白A (Actin A，Act A)、内化素(Internalins，Inl)和细胞壁水解酶(p60)[12]。其中LLO是重要的毒力因子，相对分子质量为60 000，是一种孔成形毒素，与细胞膜表面分子结合后可以直接插入细胞膜上，是李斯特菌逃逸初级吞噬和双层膜吞噬体所必须的毒力因子[13-14]。此外，LLO还可以调节巨噬细胞的细胞因子表达，从而诱导树突细胞凋亡。当李斯特菌借助细胞内吞侵入宿主细胞后被包裹在吞噬体中，随后李斯特菌分泌LLO和PI-PLC裂解吞噬体膜，从而逃逸进入宿主细胞浆中[15]。但此时，李斯特菌仍存在被宿主细胞自噬所清除的危险，而Act A在李斯特菌从宿主胞质中逃逸的过程中发挥着重要的作用[16]。Act A是一种表面蛋白，当李斯特菌逃出吞噬体进入细胞浆中时，Act A可以模仿Wiskott-Aldrich综合征蛋白(Wiskott-Aldrich syndrome protein，WASP)机制，调用宿主细胞自身的肌动蛋白，这些肌动蛋白在李斯特菌一端形成1条彗星状的尾巴，赋予李斯特菌在胞质中运动的能力，并向邻近的细胞迁移，实现细胞间侵染[17]。Act A还可以与宿主细胞中的Ena/VASP和Arp2/3复合体相结合，以此逃避宿主细胞的自噬识别系统[18]。
研究表明，英诺克李斯特菌与单增李斯特菌的关系密切，两者可能具有共同的毒力祖先，在进化过程中英诺克李斯特菌丢失了编码主要毒力因子的基因[19]。英诺克李斯特菌在基因组上与单增李斯特菌非常相似，例如两者都有300多个编码运输蛋白的基因，包括80多个碳转运蛋白，两者都含有糖酵解所必需的酶[20]。典型的英诺克李斯特菌已在进化过程中缺失了LIPI-1毒力岛，但Johnson等[21]和Moura等[22]研究发现某些天然非典型英诺克李斯特菌株却含有完整的LIPI-1毒力岛，这类菌株具有溶血活性，可以主动穿过肠上皮，传播到肝脏和脾脏，但危害程度远低于单增李斯特菌。由于英诺克李斯特菌与单增李斯特菌的亲缘关系最为接近，所以在一些食品的研究中，英诺克李斯特菌常常作为单增李斯特菌的替代菌种[23]。例如，Sorrentino等[24]发现在低pH条件下3−苯基乳酸可以抑制英诺克李斯特菌的生长。另外Tremonte等[25]发现，英诺克李斯特菌在面对酸胁迫时会显著上调普遍胁迫蛋白(Universal stress protein，USP)基因的表达，表明在抵御酸胁迫时USP基因发挥着至关重要的作用。
目前，李斯特菌引起的相关并发症尚无特效药，也鲜有其侵染后寄主基因组水平的详细反应情况的研究报道。由于典型的英诺克李斯特菌被认为是非致病菌[26]，人们对宿主与其互作的机制关注较少。鉴于英诺克李斯特菌与单增李斯特菌的亲缘关系最近，我们如果能够深入挖掘其侵染后宿主基因表达的变化情况，将为宿主调控、防御单增李斯特菌的危害提供对比与参考。利用英诺克李斯特菌可以侵染果蝇Drosophila melanogaster S2细胞的特性，我们在英诺克李斯特菌侵染果蝇S2细胞3 h后，对果蝇S2细胞进行了转录组测序。
1. 材料与方法
1.1 细胞与菌株的培养
英诺克李斯特菌(ATCC，33090，由上海海洋大学提供)用李斯特菌培养基(HB4160，海博生物技术有限公司)活化和培养。果蝇S2细胞(由中国科学院上海植物生理生态研究所提供)在含有体积分数为5%胎牛血清(FBS500-S，AusGeneX，澳大利亚)的Insect SIM SF培养基(MSF1-1，信和生物公司)中培养。
1.2 英诺克李斯特菌侵染果蝇S2细胞
果蝇S2细胞铺板培养18 h后，用50 nmol/L水溶性胆固醇(134428，BIOBERRY)孵育30 min。
试验组处理：胆固醇孵育30 min后加入英诺克李斯特菌悬浮液(约20个细菌)在28 ℃条件下孵育1.5 h。然后用pH 7.4的PBS缓冲液洗涤果蝇S2细胞，在含100 µmol/L CuSO4和10 µg/mL庆大霉素的Insect SIM SF培养基基中继续培养果蝇S2细胞3 h。对照组处理：不加英诺克李斯特菌悬浮液，其他操作同试验组。培养结束后，收集试验组和对照组果蝇S2细胞进行后续试验。
1.3 RNA的抽提及转录组测序
用TRIzol试剂(9108-1，Takara，日本)抽提感染英诺克李斯特菌后的果蝇S2细胞中的总RNA。使用2100 Bioanalyze(Agilent Technologies，美国)定量总RNA模板。通过凝胶电泳确定每个样品中总RNA的完整性。根据浓度和质量值，分别选取试验组与对照组3个最好的RNA样品用于文库制备和Illumina Novaseq 6000测序。原始配对的末端读数修整后由SeqPrep( https://github.com/jstjohn/SeqPrep)和Sickle( https://github.com/najoshi/sickle)使用默认参数进行质量控制。然后使用TopHat软件( http://tophat.cbcb.umd.edu/，版本2.1.1)以定向模式将原始的读数与参考基因组比对。根据每百万条映射片段中每千个碱基的转录本所包含的片段(Fragments per kilobase of transcript per million fragments mapped，FPKM)计算每个转录本的表达水平。使用RSEM软件( http://deweylab.biostat.wisc.edu/rsem/)定量基因丰度。R统计软件包EdgeR( http://www.bioconductor.org/packages/2.12/bioc/html/edgeR.html)用于基因差异表达分析。此外，使用Goatools( https://github.com/tanghaibao/Goatools)和KOBAS 2.1.1 ( http://kobas.cbi.pku.edu.cn/download.php)进行GO和KEGG分析。对照组和试验组样品均设置3个生物学重复。
1.4 实时定量PCR检测
对差异基因进行qPCR验证：用TRIzol试剂分别提取对照组与试验组果蝇S2细胞的总RNA，使用PrimeScriptTM RT试剂盒(RR047A，Takara，日本)反转录生成cDNA，进行qPCR检测。以果蝇的rp49基因作为内参基因。对照组和试验组样品均设置3个生物学重复。
利用qPCR对英诺克李斯特菌的增殖情况进行检测：英诺克李斯特菌侵染果蝇S2细胞3 h后，收集果蝇S2细胞，用基因组提取试剂盒(D3129-01，Magen)提取总基因组，以此为模板，以果蝇rp49基因作为内参基因，用英诺克李斯特菌的凋亡抑制蛋白(Inhibitor of apoptosis protein，IAP)基因特异性引物进行qPCR检测英诺克李斯特菌的增殖情况[27]。
1.5 数据分析
采用Microsoft Office Excel2010软件绘图。试验组和对照组的差异显著性通过t检验进行统计分析。
按以下公式计算每个基因的差异表达倍数：
$$ {\text{差异倍数}}{\rm{ = }}\frac{{\text{试验组基因表达量}}}{{\text{对照组基因表达量}}}{\text{。}} $$ 2. 结果与分析
2.1 英诺克李斯特菌在果蝇S2细胞内的增殖情况
利用qPCR对英诺克李斯特菌的增殖情况进行检测，结果表明，英诺克李斯特菌的IAP基因拷贝数剧烈上升，果蝇S2细胞中侵染了大量的英诺克李斯特菌(图1)。因此，我们在英诺克李斯特菌侵染果蝇S2细胞3 h后收集果蝇S2细胞提取总RNA用于转录组测序以及后续的qPCR验证。
![]()
图 1 英诺克李斯特菌侵染果蝇S2细胞3 h后英诺克李斯特菌增殖情况图中数据为平均值±标准差；“***”表示差异达到0.001的显著水平(t检验)Figure 1. Proliferation of Listeria innocua after its infection in Drosophila S2 cells for three hoursDatum in the figure is $\overline x $ ± SD; “***” indicates significant difference at the level of 0.001 (t test)2.2 果蝇S2细胞转录组测序数据质控与评估
对照组和试验组分别平均产生了59 125 670和60 779 885条原始数据，错误率均在0.025 6%以下(表1)。果蝇S2细胞被英诺克李斯特菌侵染前基因组的GC质量分数在50%左右，被侵染3 h后基因组GC质量分数变化不大(表1)。对照组和试验组的3次生物学重复之间的相关性系数均在0.994 9以上，表明各组之间重复性良好，可进行下一步分析。
表 1 果蝇S2细胞转录组测序数据统计表Table 1. Statistics of transcriptome sequencing data for Drosophila S2 cells分组
Grouping样本编号
Sample ID原始序列读数
No. of raw reads有效序列读数
No. of clean reads错误率/%
Error rateGC质量分数/%
GC content对照组
Control group1 59 639 278 59 073 172 0.025 6 51.16 2 61 599 522 61 100 672 0.024 3 51.54 3 56 138 212 55 717 464 0.024 2 51.67 试验组
Test group4 64 092 782 63 648 024 0.024 0 51.43 5 63 143 600 62 708 658 0.024 1 51.56 6 55 103 274 54 732 148 0.024 3 51.40 2.3 果蝇S2细胞基因表达差异
英诺克李斯特菌侵染宿主果蝇S2细胞3 h后，韦恩分析表明对照组与试验组共同表达的基因有7 808个，仅在试验组表达的基因有98个，仅在对照组表达的基因有187个。
与对照组相比，试验组共有54个基因呈现显著性差异(上/下调倍数≥2)，其中21个基因上调，33个基因下调，变化幅度最大的是DmCR46081，上调了232倍，其次是DmDef(DmDefensin)、DmDro(DmDrosomycin)、DmDpt A(DmDptericin A)和DmDpt B(DmDptericin B)等抗菌肽基因。DmCR46081基因编码rRNA前体加工蛋白，提示果蝇S2细胞被英诺克李斯特菌侵染后某些基因的转录和翻译水平大幅度提升(表2)。
表 2 试验组与对照组果蝇S2细胞基因表达的差异倍数Table 2. Fold change in gene expression of Drosophila S2 cells between test group and control group基因
Gene基因描述
Gene description差异倍数
Fold change基因
Gene基因描述
Gene description差异倍数
Fold changeDmCR46081 CR46081 232.000 DmMuc30E Mucin 30E 0.473 DmCG34330 Uncharacterized protein, transcript variant A 14.649 DmWhamy WHAMM and JMY related, transcript variant C 0.464 DmDef Defensin 9.805 DmGgt-1 Gamma-glutamyl transpeptidase, transcript variant B 0.463 DmDro Drosocin, transcript variant A 7.629 DmCG15528 Uncharacterized protein 0.454 DmDpt A Diptericin A 7.447 DmCG9279 Uncharacterized protein, transcript variant C 0.447 DmDpt B Diptericin B 5.825 DmCG7848 Uncharacterized protein 0.432 DmEdin Elevated during infection 5.272 DmMuc68Ca Mucin 68Ca 0.423 DmMtk Metchnikowin 4.888 DmCG8219 Uncharacterized protein 0.414 DmCec A2 Cecropin A2 4.174 DmCG1246 Uncharacterized protein, transcript variant E 0.403 DmAtt D Attacin D 3.878 DmTsp42Ep Tetraspanin 42Ep, transcript variant A 0.390 DmAtt A Attacin A 3.820 DmAOX2 Aldehyde oxidase 2 0.384 DmAtt B Attacin B, transcript variant A 3.558 DmCG4409 Uncharacterized protein 0.384 DmPGRP-SD Peptidoglycan recognition protein SD 3.380 DmSmydA-7 SET and MYND domain containing, arthropod-
specific, member 7, transcript variant A0.346 DmPGRP-SA Peptidoglycan recognition protein SA,
Transcript variant B3.282 DmCR32636 CR32636 0.322 DmCec B Cecropin B 3.209 DmCd Cardinal 0.314 DmCG31274 Uncharacterized protein 3.131 DmVha100-4 Vacuolar H+ ATPase 100 kD subunit 4 0.303 DmTsf1 Transferrin 1, transcript variant A 2.963 DmCR45215 CR45215 0.294 DmPsf1 Psf1 2.436 DmCG17140 Uncharacterized protein, transcript variant B 0.280 DmCG43175 Uncharacterized protein, transcript variant A 2.419 DmGr59f Gustatory receptor 59f 0.278 DmObp99a Odorant-binding protein 99a, transcript variant B 2.151 DmCG34124 Uncharacterized protein 0.276 DmSp212 Serine-peptidase 212 2.080 DmCtr1B Copper transporter 1B 0.251 DmNdl Nudel 0.495 DmCG13992 Uncharacterized protein, transcript variant C 0.250 DmCep290 Cep290, transcript variant A 0.483 DmCG10514 Uncharacterized protein 0.237 DmCyp4d14 Cyp4d14, transcript variant A 0.481 DmVha68-3 Vacuolar H+ ATPase 68 kD subunit 3 0.235 DmCc2d2a Coiled-coil and C2 domain containing 2A 0.478 DmDhd Deadhead, transcript variant A 0.231 DmCG34056 Uncharacterized protein 0.475 DmCG12105 Uncharacterized protein, transcript variant B 0.230 DmMtg Mind the gap, transcript variant G 0.473 DmCG42326 Uncharacterized protein, transcript variant D 0.149 2.4 果蝇S2细胞差异表达基因的功能注释分析
GO注释分析显示，21个上调的基因共涉及10个GO条目，可分为生物进程、细胞组分和分子功能3大类，其中富集基因数目最多的GO条目分别是胞外区(9个)、应激反应(5个)、代谢进程(4个)、免疫系统进程(4个)、多细胞生物进程(3个)、催化活性(3个)、结合(3个)、细胞进程(1个)、细胞成分组织或生物合成(1个)、胞外区组件(1个)(表3)。
表 3 英诺克李斯特菌侵染后果蝇S2细胞上调基因GO注释Table 3. GO annotation of upregulated genes in Drosophila S2 cells after Listeria innocua infection分类 Category GO条目 GO term 所含基因 Contained gene 生物进程
Biological process应激反应 Response to stimulus DmDef, DmPGRP-SD, DmPGRP-SA, DmCec A2, DmCec B 细胞成分组织或生物合成
Cellular component organization or biogenesisDmCR46081 免疫系统进程 Immune system process DmPGRP-SD, DmPGRP-SA, DmCec A2, DmCec B 代谢进程 Metabolic process DmSp212, DmCR46081, DmPGRP-SD, DmPGRP-SA 细胞进程 Cellular process DmCR46081 多细胞生物进程 Multi-organism process DmPGRP-SD, DmCec A2, DmCec B 细胞组分
Cellular component胞外区 Extracellular region DmDpt B, DmDef, DmTsf1, DmAtt B, DmAtt A, DmPGRP-SD, DmAtt D, DmCec A2, DmCec B 胞外区组件 Extracellular region part DmDef 分子功能
Molecular function结合 Binding DmObp99a, DmPGRP-SD, DmPGRP-SA 催化活性 Catalytic activity DmSp212, DmPGRP-SD, DmPGRP-SA 下调的33个基因共涉及22个GO条目，其中富集基因数最多的GO条目分别是催化活性(11个)、代谢进程(10个)、细胞膜(10个)、单组织进程(8个)、细胞膜组件(8个)、结合(8个)、细胞进程(7个)、定位(6个)、细胞(6个)、细胞组件(6个)、细胞器(4个)、细胞器组件(3个)、大分子复合物(3个)、转录活性(3个)、排毒(1个)、生物进程调控(1个)、多细胞组织进程(1个)、生物调控(1个)、应激反应(1个)、胞外区(1个)、抗氧化活性(1个)、电子转运活性(1个)(表4)。
表 4 英诺克李斯特菌侵染后果蝇S2细胞下调基因GO注释Table 4. GO annotation of downregulated genes in Drosophila S2 cells after Listeria innocua infection分类 Category GO条目 GO term 所含基因 Contained gene 生物进程
Biological process多细胞组织进程 Multicellular organismal process DmGr59f 生物进程调控 Regulation of biological process DmDhd 排毒 Detoxification DmCd 应激反应 Response to stimulus DmCd 生物调控 Biological regulation DmDhd 定位 Localization DmVha68-3, DmCG9279, DmCG17140, DmCG8219, DmCtr1B, DmVha100-4 单组织进程 Single-organism process DmVha68-3, DmCG9279, DmCG7848, DmAOX2, DmCd, DmDhd, DmVha100-4, DmCyp4d14 代谢进程 Metabolic process DmCG15528, DmVha68-3, DmNdl, DmCG7848, DmAOX2, DmMtg, DmCd, DmDhd, DmGgt-1, DmCyp4d14 细胞进程 Cellular process DmCG15528, DmVha68-3, DmCG9279, DmCG7848, DmCd, DmDhd, DmGgt-1 细胞组分
Cellular component细胞 Cell DmCG9279, DmCG17140, DmCG8219, DmGr59f, DmDhd, DmVha100-4 细胞器 Organelle DmCG9279, DmCG17140, DmDhd, DmVha100-4 细胞器组件 Organelle part DmCG9279, DmCG17140, DmVha100-4 细胞组件 Cell part DmCG9279, DmCG17140, DmCG8219, DmGr59f, DmDhd, DmVha100-4 细胞膜 Membrane DmVha68-3, DmNdl, DmCG34056, DmCG1246, DmCG17140,
DmTsp42Ep, DmCd, DmGr59f, DmCtr1B, DmVha100-4细胞膜组件 Membrane part DmVha68-3, DmNdl, DmCG1246, DmTsp42Ep, DmCd, DmGr59f, DmCtr1B, DmVha100-4 胞外区 Extracellular region DmMtg 大分子复合物 Macromolecular complex DmVha68-3, DmCG9279, DmVha100-4 分子功能
Molecular function电子转运活性 Electron carrier activity DmAOX2 转录活性 Transporter activity DmVha68-3, DmCtr1B, DmVha100-4 结合 Binding DmVha68-3, DmCG9279, DmCG7848, DmAOX2, DmMtg, DmCd, DmCG8219, DmCyp4d14 催化活性 Catalytic activity DmCG15528, DmVha68-3, DmNdl, DmCG7848, DmAOX2, DmCG34056, DmCd, DmDhd, DmGgt-1, DmCG10514, DmCyp4d14 抗氧化活性 Antioxidant activity DmCd 结果显示，上调基因多集中在胞外区，下调基因多集中在代谢进程、细胞膜和催化活性，说明被英诺克李斯特菌侵染后果蝇S2细胞的某些代谢反应下降。
2.5 果蝇S2细胞差异表达基因的富集分析
2.5.1 上调基因KEGG富集
KEGG富集分析表明，显著上调的21个基因只注释到1个KEGG 通路，即Toll/Imd通路(P<0.001)，包括4个果蝇抗菌肽基因DmDpt A、DmDef、DmCec A2(Cecropin A2)和DmCec B(Cecropin B)。
2.5.2 下调基因KEGG富集
KEGG富集分析结果(表5)表明，下调的33个基因注释到23个KEGG 通路。
表 5 英诺克李斯特菌侵染后果蝇S2细胞下调基因KEGG富集表Table 5. KEGG enrichment of downregulated genes in Drosophila S2 cells after Listeria innocua infection通路 Pathway P 所含基因 Contained gene 霍乱弧菌感染 Vibrio cholerae infection 0.033 DmVha68-3, DmVha100-4 类风湿关节炎 Rheumatoid arthritis 0.034 DmVha68-3, DmVha100-4 幽门螺杆菌感染中的上皮细胞信号传导
Epithelial cell signaling in Helicobacter pylori infection0.036 DmVha68-3, DmVha100-4 集合管酸分泌 Collecting duct acid secretion 0.037 DmVha68-3, DmVha100-4 突触小泡循环 Synaptic vesicle cycle 0.042 DmVha68-3, DmVha100-4 咖啡因代谢 Caffeine metabolism 0.043 DmAOX2 牛磺酸和牛磺酸的代谢 Taurine and hypotaurine metabolism 0.049 DmGgt-1 吞噬体 Phagosome 0.049 DmVha68-3, DmVha100-4 花生四烯酸代谢 Arachidonic acid metabolism 0.072 DmGgt-1 氧化磷酸化 Oxidative phosphorylation 0.078 DmVha68-3, DmVha100-4 矿物质吸收 Mineral absorption 0.082 DmCtr1B 粘蛋白型 O−聚糖的生物合成 Mucin type O-glycan biosynthesis 0.089 DmCG34056 药物代谢−其他酶 Drug metabolism - other enzymes 0.155 DmAOX2 NOD样受体信号通路 NOD-like receptor signaling pathway 0.158 DmDhd 丙酮酸代谢 Pyruvate metabolism 0.163 DmCG7848 铂耐药 Platinum drug resistance 0.182 DmCtr1B 过氧化物酶体 Peroxisome 0.197 DmAOX2 结核 Tuberculosis 0.198 DmVha100-4 谷胱甘肽代谢 Glutathione metabolism 0.204 DmGgt-1 流体剪切应力与动脉粥样硬化 Fluid shear stress and atherosclerosis 0.209 DmDhd mTOR信号通路 mTOR signaling pathway 0.216 DmVha68-3 溶酶体 Lysosome 0.247 DmVha100-4 嘌呤代谢 Purine metabolism 0.262 DmAOX2 2.6 免疫相关通路中的基因分析
英诺克李斯特菌侵染果蝇S2细胞后的转录组测序数据经GO注释和KEGG富集分析显示，显著上调的基因主要与免疫相关。从24个KEGG通路中筛选出10个与免疫相关的KEGG 通路，分别是Toll/Imd通路、类风湿关节炎、吞噬体、霍乱弧菌感染、幽门螺杆菌感染中的上皮细胞信号传导、NOD样受体信号通路、铂耐药、结核、流体剪切应力与动脉粥样硬化、mTOR信号通路，共涉及到8个基因，其中上调幅度最大的是DmDef，上调了9.805倍，下调幅度最大的是DmDhd，下调了4.329倍。
2.7 果蝇S2细胞转录组中抗菌肽的qPCR验证
英诺克李斯特菌侵染果蝇S2细胞后，变化最为显著的基因类群为抗菌肽。转录组分析显示，英诺克李斯特菌侵染果蝇S2细胞3 h后共10个抗菌肽基因显著上调，包括DmDef、DmDro、DmDpt A、DmDpt B、DmCec A2、DmCec B、DmMtk(DmMetchnikowin)、DmAtt D(DmAttacin D)、DmAtt A(DmAttacin A)、DmAtt B(DmAttacin B)，其中上调幅度最大的是DmDef(表2)。
我们选取9个抗菌肽基因，利用qPCR进行了验证，其中上调幅度最大的是DmMtk，上调8.180倍，DmAtt A 上调7.533倍，DmDro 上调7.204倍， DmDef 上调4.569倍(图2)，与转录组测序结果的趋势基本一致，说明英诺克李斯特菌侵染果蝇S2细胞后，抗菌肽基因受到显著的诱导表达。
![]()
图 2 免疫相关基因表达的qPCR验证结果Figure 2. qPCR verification results for expression of immune related genes3. 讨论与结论
已知的果蝇免疫信号通路有3条，分别是Toll通路、Imd通路和JAK/STAT通路。果蝇虽然没有获得性免疫，但是依然能对不同病原菌的入侵做出不同的先天免疫反应，革兰阳性菌和真菌所引起的免疫反应主要由Toll通路介导，而革兰阴性菌和细胞壁含DAP型肽聚糖的阳性菌引起的免疫反应主要由Imd通路介导[28-29]。Toll通路或Imd通路被激活后可诱导下游抗菌肽的表达[30]。已知的果蝇抗菌肽主要有7种，分别是Attacin、Diptericin、Drosocin、Defensin、Cecropin、Metchnikowin和Drosomycin。抗菌肽的表达由Toll和Imd通路调控，若Toll通路和Imd通路缺失，将不能诱导抗菌肽的表达[31-33]。
英诺克李斯特菌属于革兰阳性菌，其细胞壁肽聚糖成分主要为DAP型肽聚糖[34]，DAP型肽聚糖可以激活果蝇的Imd通路。已有研究表明，果蝇的PGRP-SD作为胞外受体不参与Toll通路，而是参与Imd通路的激活[35]。本研究发现DmPGRP-SD出现显著上调，提示英诺克李斯特菌侵染果蝇S2细胞后Imd通路可能被激活。抗菌肽Attacin、Cecropin和Diptericin的表达主要受Imd通路的调控[36-37]，本研究结果显示果蝇的抗菌肽基因DmAtt A、DmAtt B、DmAtt D 和 DmCec A2 和 DmCec B均出现显著上调，显示了Imd通路的激活。而抗菌肽Drosomycin、Defensin和Metchnikowin的表达主要受Toll通路调控[38]，本研究结果中显示果蝇的抗菌肽基因DmDro、DmDef和DmMtk也出现了显著上调，提示英诺克李斯特菌可能同时引起了Imd和Toll通路的激活，但具体是Toll通路还是Imd通路起主要作用，则需要进一步的研究。抗菌肽被认为具有广谱性，但Hanson等[39]研究发现抗菌肽针对某些病原微生物具有特异性，Imd通路调控的抗菌肽对单增李斯特菌并不起作用[34]。Hanson等[39]研究还发现果蝇的spz基因缺失后对英诺克李斯特菌表现出极高的敏感性，提示在抵御英诺克李斯特菌的过程中，Toll通路调控的抗菌肽起重要作用。鉴于上述研究结果和单增李斯特菌与英诺克李斯特菌相近的亲缘关系，猜测宿主在抵御英诺克李斯特菌和单增李斯特菌侵染的过程中Toll通路调控的抗菌肽可能起到关键作用，但Toll通路的激活机制需要进一步深入研究。
英诺克李斯特菌作为一种胞内寄生菌，通过内吞作用侵入宿主细胞后会被包裹在吞噬体中，英诺克李斯特菌在宿主细胞中的增殖需要它从宿主的吞噬体中逃逸至细胞质中，我们注意到液泡型H+离子腺苷三磷酸酶(Vacuolar H+-ATPase，V-ATPase)基因DmVha68-3和DmVha100-4均出现了显著下调，KEGG富集结果显示它们富集到吞噬体等通路中。DmVha68-3编码果蝇V-ATPase-A亚基，DmVha100-4编码V-ATPase-a亚基，为V-ATPase复合体形成的关键亚基[40]，V-ATPase广泛存在于细胞器膜上，如内体、溶酶体及吞噬体膜，可以将H+泵入这些细胞器中以维持这些细胞器的相对酸性内环境，促进细胞器对内吞或是细胞自噬包裹的底物或是入侵的微生物进行降解[41-42]。单增李斯特菌的毒力因子LLO正常发挥作用时需要酸性环境，在pH 5.6时LLO的溶血活性最高，当pH超过7.0时LLO的活性基本丧失[43]。典型的英诺李斯特菌缺失毒力岛，而其侵染果蝇S2细胞后上述2个DmV-ATPase基因的下调是否跟毒力岛缺失有关目前尚不清楚，值得我们今后深入思考。此外，是否可以通过调控宿主细胞的V-ATPase活性从而抑制单增李斯特菌毒力因子LLO的活性对其进行防控，也是有益的思考。
-
![]()
图 1 英诺克李斯特菌侵染果蝇S2细胞3 h后英诺克李斯特菌增殖情况
图中数据为平均值±标准差；“***”表示差异达到0.001的显著水平(t检验)
Figure 1. Proliferation of Listeria innocua after its infection in Drosophila S2 cells for three hours
Datum in the figure is $\overline x $ ± SD; “***” indicates significant difference at the level of 0.001 (t test)
![]()
图 2 免疫相关基因表达的qPCR验证结果
Figure 2. qPCR verification results for expression of immune related genes
表 1 果蝇S2细胞转录组测序数据统计表
Table 1 Statistics of transcriptome sequencing data for Drosophila S2 cells
分组
Grouping样本编号
Sample ID原始序列读数
No. of raw reads有效序列读数
No. of clean reads错误率/%
Error rateGC质量分数/%
GC content对照组
Control group1 59 639 278 59 073 172 0.025 6 51.16 2 61 599 522 61 100 672 0.024 3 51.54 3 56 138 212 55 717 464 0.024 2 51.67 试验组
Test group4 64 092 782 63 648 024 0.024 0 51.43 5 63 143 600 62 708 658 0.024 1 51.56 6 55 103 274 54 732 148 0.024 3 51.40 
下载: 导出CSV
表 2 试验组与对照组果蝇S2细胞基因表达的差异倍数
Table 2 Fold change in gene expression of Drosophila S2 cells between test group and control group
基因
Gene基因描述
Gene description差异倍数
Fold change基因
Gene基因描述
Gene description差异倍数
Fold changeDmCR46081 CR46081 232.000 DmMuc30E Mucin 30E 0.473 DmCG34330 Uncharacterized protein, transcript variant A 14.649 DmWhamy WHAMM and JMY related, transcript variant C 0.464 DmDef Defensin 9.805 DmGgt-1 Gamma-glutamyl transpeptidase, transcript variant B 0.463 DmDro Drosocin, transcript variant A 7.629 DmCG15528 Uncharacterized protein 0.454 DmDpt A Diptericin A 7.447 DmCG9279 Uncharacterized protein, transcript variant C 0.447 DmDpt B Diptericin B 5.825 DmCG7848 Uncharacterized protein 0.432 DmEdin Elevated during infection 5.272 DmMuc68Ca Mucin 68Ca 0.423 DmMtk Metchnikowin 4.888 DmCG8219 Uncharacterized protein 0.414 DmCec A2 Cecropin A2 4.174 DmCG1246 Uncharacterized protein, transcript variant E 0.403 DmAtt D Attacin D 3.878 DmTsp42Ep Tetraspanin 42Ep, transcript variant A 0.390 DmAtt A Attacin A 3.820 DmAOX2 Aldehyde oxidase 2 0.384 DmAtt B Attacin B, transcript variant A 3.558 DmCG4409 Uncharacterized protein 0.384 DmPGRP-SD Peptidoglycan recognition protein SD 3.380 DmSmydA-7 SET and MYND domain containing, arthropod-
specific, member 7, transcript variant A0.346 DmPGRP-SA Peptidoglycan recognition protein SA,
Transcript variant B3.282 DmCR32636 CR32636 0.322 DmCec B Cecropin B 3.209 DmCd Cardinal 0.314 DmCG31274 Uncharacterized protein 3.131 DmVha100-4 Vacuolar H+ ATPase 100 kD subunit 4 0.303 DmTsf1 Transferrin 1, transcript variant A 2.963 DmCR45215 CR45215 0.294 DmPsf1 Psf1 2.436 DmCG17140 Uncharacterized protein, transcript variant B 0.280 DmCG43175 Uncharacterized protein, transcript variant A 2.419 DmGr59f Gustatory receptor 59f 0.278 DmObp99a Odorant-binding protein 99a, transcript variant B 2.151 DmCG34124 Uncharacterized protein 0.276 DmSp212 Serine-peptidase 212 2.080 DmCtr1B Copper transporter 1B 0.251 DmNdl Nudel 0.495 DmCG13992 Uncharacterized protein, transcript variant C 0.250 DmCep290 Cep290, transcript variant A 0.483 DmCG10514 Uncharacterized protein 0.237 DmCyp4d14 Cyp4d14, transcript variant A 0.481 DmVha68-3 Vacuolar H+ ATPase 68 kD subunit 3 0.235 DmCc2d2a Coiled-coil and C2 domain containing 2A 0.478 DmDhd Deadhead, transcript variant A 0.231 DmCG34056 Uncharacterized protein 0.475 DmCG12105 Uncharacterized protein, transcript variant B 0.230 DmMtg Mind the gap, transcript variant G 0.473 DmCG42326 Uncharacterized protein, transcript variant D 0.149
下载: 导出CSV
表 3 英诺克李斯特菌侵染后果蝇S2细胞上调基因GO注释
Table 3 GO annotation of upregulated genes in Drosophila S2 cells after Listeria innocua infection
分类 Category GO条目 GO term 所含基因 Contained gene 生物进程
Biological process应激反应 Response to stimulus DmDef, DmPGRP-SD, DmPGRP-SA, DmCec A2, DmCec B 细胞成分组织或生物合成
Cellular component organization or biogenesisDmCR46081 免疫系统进程 Immune system process DmPGRP-SD, DmPGRP-SA, DmCec A2, DmCec B 代谢进程 Metabolic process DmSp212, DmCR46081, DmPGRP-SD, DmPGRP-SA 细胞进程 Cellular process DmCR46081 多细胞生物进程 Multi-organism process DmPGRP-SD, DmCec A2, DmCec B 细胞组分
Cellular component胞外区 Extracellular region DmDpt B, DmDef, DmTsf1, DmAtt B, DmAtt A, DmPGRP-SD, DmAtt D, DmCec A2, DmCec B 胞外区组件 Extracellular region part DmDef 分子功能
Molecular function结合 Binding DmObp99a, DmPGRP-SD, DmPGRP-SA 催化活性 Catalytic activity DmSp212, DmPGRP-SD, DmPGRP-SA
下载: 导出CSV
表 4 英诺克李斯特菌侵染后果蝇S2细胞下调基因GO注释
Table 4 GO annotation of downregulated genes in Drosophila S2 cells after Listeria innocua infection
分类 Category GO条目 GO term 所含基因 Contained gene 生物进程
Biological process多细胞组织进程 Multicellular organismal process DmGr59f 生物进程调控 Regulation of biological process DmDhd 排毒 Detoxification DmCd 应激反应 Response to stimulus DmCd 生物调控 Biological regulation DmDhd 定位 Localization DmVha68-3, DmCG9279, DmCG17140, DmCG8219, DmCtr1B, DmVha100-4 单组织进程 Single-organism process DmVha68-3, DmCG9279, DmCG7848, DmAOX2, DmCd, DmDhd, DmVha100-4, DmCyp4d14 代谢进程 Metabolic process DmCG15528, DmVha68-3, DmNdl, DmCG7848, DmAOX2, DmMtg, DmCd, DmDhd, DmGgt-1, DmCyp4d14 细胞进程 Cellular process DmCG15528, DmVha68-3, DmCG9279, DmCG7848, DmCd, DmDhd, DmGgt-1 细胞组分
Cellular component细胞 Cell DmCG9279, DmCG17140, DmCG8219, DmGr59f, DmDhd, DmVha100-4 细胞器 Organelle DmCG9279, DmCG17140, DmDhd, DmVha100-4 细胞器组件 Organelle part DmCG9279, DmCG17140, DmVha100-4 细胞组件 Cell part DmCG9279, DmCG17140, DmCG8219, DmGr59f, DmDhd, DmVha100-4 细胞膜 Membrane DmVha68-3, DmNdl, DmCG34056, DmCG1246, DmCG17140,
DmTsp42Ep, DmCd, DmGr59f, DmCtr1B, DmVha100-4细胞膜组件 Membrane part DmVha68-3, DmNdl, DmCG1246, DmTsp42Ep, DmCd, DmGr59f, DmCtr1B, DmVha100-4 胞外区 Extracellular region DmMtg 大分子复合物 Macromolecular complex DmVha68-3, DmCG9279, DmVha100-4 分子功能
Molecular function电子转运活性 Electron carrier activity DmAOX2 转录活性 Transporter activity DmVha68-3, DmCtr1B, DmVha100-4 结合 Binding DmVha68-3, DmCG9279, DmCG7848, DmAOX2, DmMtg, DmCd, DmCG8219, DmCyp4d14 催化活性 Catalytic activity DmCG15528, DmVha68-3, DmNdl, DmCG7848, DmAOX2, DmCG34056, DmCd, DmDhd, DmGgt-1, DmCG10514, DmCyp4d14 抗氧化活性 Antioxidant activity DmCd
下载: 导出CSV
表 5 英诺克李斯特菌侵染后果蝇S2细胞下调基因KEGG富集表
Table 5 KEGG enrichment of downregulated genes in Drosophila S2 cells after Listeria innocua infection
通路 Pathway P 所含基因 Contained gene 霍乱弧菌感染 Vibrio cholerae infection 0.033 DmVha68-3, DmVha100-4 类风湿关节炎 Rheumatoid arthritis 0.034 DmVha68-3, DmVha100-4 幽门螺杆菌感染中的上皮细胞信号传导
Epithelial cell signaling in Helicobacter pylori infection0.036 DmVha68-3, DmVha100-4 集合管酸分泌 Collecting duct acid secretion 0.037 DmVha68-3, DmVha100-4 突触小泡循环 Synaptic vesicle cycle 0.042 DmVha68-3, DmVha100-4 咖啡因代谢 Caffeine metabolism 0.043 DmAOX2 牛磺酸和牛磺酸的代谢 Taurine and hypotaurine metabolism 0.049 DmGgt-1 吞噬体 Phagosome 0.049 DmVha68-3, DmVha100-4 花生四烯酸代谢 Arachidonic acid metabolism 0.072 DmGgt-1 氧化磷酸化 Oxidative phosphorylation 0.078 DmVha68-3, DmVha100-4 矿物质吸收 Mineral absorption 0.082 DmCtr1B 粘蛋白型 O−聚糖的生物合成 Mucin type O-glycan biosynthesis 0.089 DmCG34056 药物代谢−其他酶 Drug metabolism - other enzymes 0.155 DmAOX2 NOD样受体信号通路 NOD-like receptor signaling pathway 0.158 DmDhd 丙酮酸代谢 Pyruvate metabolism 0.163 DmCG7848 铂耐药 Platinum drug resistance 0.182 DmCtr1B 过氧化物酶体 Peroxisome 0.197 DmAOX2 结核 Tuberculosis 0.198 DmVha100-4 谷胱甘肽代谢 Glutathione metabolism 0.204 DmGgt-1 流体剪切应力与动脉粥样硬化 Fluid shear stress and atherosclerosis 0.209 DmDhd mTOR信号通路 mTOR signaling pathway 0.216 DmVha68-3 溶酶体 Lysosome 0.247 DmVha100-4 嘌呤代谢 Purine metabolism 0.262 DmAOX2
下载: 导出CSV
-
[1] 李翠云. 单核细胞李斯特菌研究近况[J]. 中国热带医学, 2010, 10(1): 120-122. [2] 蔡雪薛. 单核细胞增生李斯特菌毒力相关基因的筛选鉴定及其功能研究[D]. 扬州: 扬州大学, 2017. [3] VÁZQUEZ-BOLAND J A, DOMÍNGUEZ-BERNAL G, GONZÁLEZ-ZORN B, et al. Pathogenicity islands and virulence evolution in Listeria[J]. Microbes Infect, 2001, 3(7): 571-584. doi: 10.1016/S1286-4579(01)01413-7
[4] GAILLARD J L, BERCHE P, SANSONETTI P. Transposon mutagenesis as a tool to study the role of hemolysin in the virulence of Listeria monocytogenes[J]. Infect Immun, 1986, 52(1): 50-55. doi: 10.1128/IAI.52.1.50-55.1986
[5] 陈健舜, 江玲丽, 方维焕. 李斯特菌毒力因子及其进化[J]. 微生物学报, 2007(4): 738-742. doi: 10.3321/j.issn:0001-6209.2007.04.035 [6] ORSI R H, WIEDMANN M. Characteristics and distribution of Listeriaspp., including Listeria species newly described since 2009[J]. Appl Microbiol Biotechnol, 2016, 100(12): 5273-5287. doi: 10.1007/s00253-016-7552-2
[7] JUNTTILA J R, NIEMELÄ S I, HIRN J. Minimum growth temperatures of Listeria monocytogenes and non‐haemolytic Listeria[J]. J Appl Bacteriol, 1988, 65(4): 321-327. doi: 10.1111/j.1365-2672.1988.tb01898.x
[8] 王旭. 基于比较基因组学和转录组学技术揭示单增李斯特菌毒力因子的研究[D]. 上海: 上海海洋大学, 2016. [9] 杜秀嫒. 李斯特菌致病性的研究及部分基因的功能分析[D]. 兰州: 甘肃农业大学, 2018. [10] SCALLAN E, HOEKSTRA R M, ANGULO F J, et al. Foodborne illness acquired in the United States: Major pathogens[J]. Emerg Infect Dis, 2011, 17(1): 7-15. doi: 10.3201/eid1701.P11101
[11] LI W, BAI L, FU P, et al. The epidemiology of Listeria monocytogenes in China[J]. Foodborne Pathog Dis, 2018, 15(8): 459-466. doi: 10.1089/fpd.2017.2409
[12] 郎需龙, 王兴龙. 李斯特菌毒力因子及其调控机制的研究进展[J]. 中国畜牧兽医, 2011, 38(3): 195-198. [13] ALVAREZ-DOMINGUEZ C, ROBERTS R, STAHL P D. Internalized Listeria monocytogenes modulates intracellular trafficking and delays maturation of the phagosome[J]. J Cell Sci, 1997, 110(Pt 6): 731-743.
[14] MARQUIS H, GOLDFINE H, PORTNOY D A. Proteolytic pathways of activation and degradation of a bacterial phospholipase C during intracellular infection by Listeria monocytogenes[J]. J Cell Biol, 1997, 137(6): 1381-1392. doi: 10.1083/jcb.137.6.1381
[15] PIZARRO-CERDÁ J, COSSART P. Subversion of cellular functions by Listeria monocytogenes[J]. J Pathol, 2006, 208(2): 215-223. doi: 10.1002/path.1888
[16] KOCKS C, GOUIN E, TABOURET M, et al. L. monocytogenes-induced actin assembly requires the actA gene product, a surface protein[J]. Cell, 1992, 68(3): 521-531.
[17] COSSART P, TOLEDO-ARANA A. Listeria monocytogenes, a unique model in infection biology: An overview[J]. Microbes Infect, 2008, 10(9): 1041-1050. doi: 10.1016/j.micinf.2008.07.043
[18] CHONG R, SWISS R, BRIONES G, et al. Regulatory mimicry in Listeria monocytogenes actin-based motility[J]. Cell Host Microbe, 2009, 6(3): 268-278. doi: 10.1016/j.chom.2009.08.006
[19] DEN BAKKER H C, CUMMINGS C A, FERREIRA V, et al. Comparative genomics of the bacterial genus Listeria: Genome evolution is characterized by limited gene acquisition and limited gene loss[J]. BMC Genomics, 2010, 11(1): 688. doi: 10.1186/1471-2164-11-688
[20] MILILLO S R, FRIEDLY E C, SALDIVAR J C, et al. A Review of the ecology, genomics, and stress response of Listeria innocua and Listeria monocytogenes[J]. Crit Rev Food Sci Nutr, 2012, 52(8): 712-725. doi: 10.1080/10408398.2010.507909
[21] JOHNSON J, JINNEMAN K, STELMA G, et al. Natural atypical Listeria innocua strains with Listeria monocytogenes pathogenicity island 1 genes[J]. Appl Environ Microbiol, 2004, 7(7): 4256-4266.
[22] MOURA A, DISSON O, LAVINA M, et al. Atypical hemolytic Listeria innocua isolates are virulent, albeit less than Listeria monocytogenes[J]. Infect Immun, 2019, 87(4): e00758-18.
[23] NYHAN L, JOHNSON N, BEGLEY M, et al. Comparison of predicted and impedance determined growth of Listeria innocua in complex food matrices[J]. Food Microbiol, 2020, 87: 103381. doi: 10.1016/j.fm.2019.103381
[24] SORRENTINO E, TREMONTE P, SUCCI M, et al. Detection of antilisterial activity of 3-phenyllactic acid using Listeria innocua as a model[J]. Front Microbiol, 2018, 9: 1373.
[25] TREMONTE P, SUCCI M, COPPOLA R, et al. Homology-based modeling of universal stress protein from Listeria innocua up-regulated under acid stress conditions[J]. Front Microbiol, 2016, 7: 1998.
[26] YILDIRIM Z, YERLIKAYA S, ÖNCÜL N, et al. Inhibitory effect of lactococcin BZ against Listeria innocua and indigenous microbiota of fresh beef[J]. Food Technol Biotechnol, 2016, 54(3): 317-323.
[27] LIU D. Identification, subtyping and virulence determination of Listeria monocytogenes, an important foodborne pathogen[J]. J Med Microbiol, 2006, 55(6): 645-659. doi: 10.1099/jmm.0.46495-0
[28] 刘甜, 罗开珺. 果蝇Toll和Imd信号通路中的功能结构域[J]. 环境昆虫学报, 2011, 33(3): 388-395. doi: 10.3969/j.issn.1674-0858.2011.03.015 [29] STUART L M, EZEKOWITZ R A. Phagocytosis and comparative innate immunity: Learning on the fly[J]. Nat Rev Immunol, 2008, 8(2): 131-141. doi: 10.1038/nri2240
[30] TANJI T, HU X, WEBER A N R, et al. Toll and Imd pathways synergistically activate an innate immune response in Drosophila melanogaster[J]. Mol Cell Biol, 2007, 27(12): 4578-4588. doi: 10.1128/MCB.01814-06
[31] VODOVAR N, VINALS M, LIEHL P, et al. Drosophila host defense after oral infection by an entomopathogenic Pseudomonas species[J]. Proc Natl Acad Sci USA, 2005, 102(32): 11414-11419. doi: 10.1073/pnas.0502240102
[32] TZOU P, REICHHART J M, LEMAITRE B. Constitutive expression of a single antimicrobial peptide can restore wild-type resistance to infection in immunodeficient Drosophila mutants[J]. Proc Natl Acad Sci USA, 2002, 99(4): 2152-2157. doi: 10.1073/pnas.042411999
[33] GOTTAR M, GOBERT V, MICHEL T, et al. The Drosophila immune response against Gram-negative bacteria is mediated by a peptidoglycan recognition protein[J]. Nature, 2002, 416(6881): 640-644.
[34] 陈康康, 吕志强. 昆虫肽聚糖识别蛋白研究进展[J]. 昆虫学报, 2014, 57(8): 969-978. [35] IATSENKO I, KONDO S, MENGIN-LECREULX D, et al. PGRP-SD, an extracellular pattern-recognition receptor, enhances peptidoglycan-mediated activation of the Drosophila Imd pathway[J]. Immunity, 2016, 45(5): 1013-1023. doi: 10.1016/j.immuni.2016.10.029
[36] LEMAITRE B, REICHHART J M, HOFFMANN J A. Drosophila host defense: Differential induction of antimicrobial peptide genes after infection by various classes of microorganisms[J]. Proc Natl Acad Sci USA, 1997, 94(26): 14614-14619. doi: 10.1073/pnas.94.26.14614
[37] IMLER J, BULET P. Antimicrobial peptides in Drosophila: Structures, activities and gene regulation[J]. Chem Immunol Allergy, 2005, 86: 86648.
[38] 程廷才, 王根洪, 李娟, 等. 昆虫抗菌肽基因表达调控机理[J]. 蚕学通讯, 2004, 24(3): 20-26. doi: 10.3969/j.issn.1006-0561.2004.03.008 [39] HANSON M A, DOSTÁLOVÁ A, CERONI C, et al. Synergy and remarkable specificity of antimicrobial peptides in vivo using a systematic knockout approach[J]. eLife, 2019, 8: 44341.
[40] ALLAN A K, DU J, DAVIES S A, et al. Genome-wide survey of V-ATPase genes in Drosophila reveals a conserved renal phenotype for lethal alleles[J]. Physiol Genomics, 2005, 22(2): 128-138. doi: 10.1152/physiolgenomics.00233.2004
[41] 游海燕, 邓云, 覃文新. V-ATPases的功能及其抑制剂研究进展[J]. 生命科学, 2009, 21(4): 499-503. [42] HOLLIDAY L S. Vacuolar H+-ATPase: An essential multitasking enzyme in physiology and pathophysiology[J]. New J Sci, 2014, 2014: 1-21.
[43] 谭炳乾, 何启盖, 肖军, 等. 单核细胞增多性李斯特菌hlyA基因序列及溶血素活性测定[J]. 中国兽医学报, 2009, 29(2): 161-165. -
期刊类型引用(2)
1. 邓健浩，刘力，舒婵，吴坤，吴坤钟，姚兵华，李康，田铃，黄志君. 单增李斯特菌、格氏李斯特菌或威尔斯李斯特菌侵染果蝇S2细胞后的转录分析. 华南农业大学学报. 2025(01): 12-24 . 
本站查看
2. 徐佳卉，刘华伟，孙小桐，唐章晨，黄敏，赵萍. 家蚕核转录因子NF-κB对不同细胞系中免疫相关基因表达的调控分析. 生物工程学报. 2024(10): 3649-3665 . 百度学术
其他类型引用(0)
计量
- 文章访问数: 1827
- HTML全文浏览量: 15
- PDF下载量: 2423
- 被引次数: 2
