留言板

尊敬的读者、作者、审稿人, 关于本刊的投稿、审稿、编辑和出版的任何问题, 您可以本页添加留言。我们将尽快给您答复。谢谢您的支持!

姓名
邮箱
手机号码
标题
留言内容
验证码

Dual RNA-seq provides novel insight into the roles of dksA from Pseudomonas plecoglossicida in pathogen-host interactions with large yellow croakers (Larimichthys crocea)

Lu-Ying Wang Zi-Xu Liu Ling-Min Zhao Li-Xing Huang Ying-Xue Qin Yong-Quan Su Wei-Qiang Zheng Fan Wang Qing-Pi Yan

Lu-Ying Wang, Zi-Xu Liu, Ling-Min Zhao, Li-Xing Huang, Ying-Xue Qin, Yong-Quan Su, Wei-Qiang Zheng, Fan Wang, Qing-Pi Yan. Dual RNA-seq provides novel insight into the roles of dksA from Pseudomonas plecoglossicida in pathogen-host interactions with large yellow croakers (Larimichthys crocea). Zoological Research, 2020, 41(4): 410-422. doi: 10.24272/j.issn.2095-8137.2020.048
Citation: Lu-Ying Wang, Zi-Xu Liu, Ling-Min Zhao, Li-Xing Huang, Ying-Xue Qin, Yong-Quan Su, Wei-Qiang Zheng, Fan Wang, Qing-Pi Yan. Dual RNA-seq provides novel insight into the roles of dksA from Pseudomonas plecoglossicida in pathogen-host interactions with large yellow croakers (Larimichthys crocea). Zoological Research, 2020, 41(4): 410-422. doi: 10.24272/j.issn.2095-8137.2020.048

基于Dual RNA-seq研究dksA基因在变形假单胞菌与大黄鱼的病原体-宿主相互作用中的功能

doi: 10.24272/j.issn.2095-8137.2020.048

Dual RNA-seq provides novel insight into the roles of dksA from Pseudomonas plecoglossicida in pathogen-host interactions with large yellow croakers (Larimichthys crocea)

Funds: This work was supported by the National Natural Science Foundation of China (31672694), Fujian Provincial Special Fund for Marine and Fishery Protection and Development (MCZ[2019]062), and Open Fund of Fujian Province Key Laboratory of Special Aquatic Formula Feed (2019-01)
More Information
  • 摘要: 变形假单胞菌(Pseudomonas plecoglossicida)是革兰氏阴性菌,呈杆状,有鞭毛,该菌在大黄鱼(Larimichthys crocea)养殖过程中引起内脏白点病并有非常高的死亡率,可造成严重的经济损失。经互作转录组数据分析和实时定量PCR验证,dksA基因在变形假单胞菌感染斜带石斑鱼(Epinephelus coioides)48小时后的表达量上调显著(log2FC=3.12 ,P<0.001)。利用RNAi技术显著降低了变形假单胞菌中基因dksA的表达量,最佳沉默效率达到96.23%。分别用变形假单胞菌野生株与效果最好的沉默菌株注射感染健康大黄鱼,观察大黄鱼的存活情况,并于注射后48小时取脾脏进行dual RNA-seq。与野生型毒株相比,感染沉默菌株的大黄鱼内脏白点病症状减轻,死亡时间延迟了2 天且死亡率降低20%。经dual RNA-seq 得知人工感染后病原和宿主的转录组发生了显著变化。dksA基因的沉默导致病原菌的细胞组分,鞭毛和核糖体组装相关基因大量下调,而fliC基因的显著上调可能是变形假单胞菌毒力维持的一种途径。GO和KEGG结果表明沉默菌株感染的大黄鱼的多个与免疫相关通路中炎性因子基因大量下调,在这些免疫途径中,受影响的基因与多个免疫应答过程相关。研究结果表明,dksA是变形假单胞菌的毒力基因,相比于野生株感染,沉默菌株感染的大黄鱼免疫反应较弱。
  • Figure  1.  Construction and characterization of dksA-RNAi strains of P. plecoglossicida

    A: Relative expression of gene dksA (in vivo/in vitro). B: dksA expression level in five mutant strains. C: Growth curve of two strains of P. plecoglossicida.

    Figure  2.  Infection of L. crocea by two strains of P. plecoglossicida

    A: Survival rate of L. crocea infected by P. plecoglossicida. B: Spleen appearance of infected L. crocea. C: Relative pathogen load of RNAi strain in spleen of L. crocea compared to that of wild-type strain. D: Expression levels of dksA gene in two strains of P. plecoglossicida in vitro (0 hpi) and in vivo (24, 48, 72 and 96 hpi).

    Figure  3.  DEG enrichment analysis of pathogen transcriptome data

    A: Volcano plot obtained from edgeR analysis of P. plecoglossicida transcriptome. B: Heat map of DEGs from transcripts of pathogen in host spleen (FDR<0.05, |log2FC|≥1), C: Five up-regulated and down-regulated DEGs were randomly selected for qRT-PCR verification.

    Figure  4.  GO and KEGG enrichment analysis of pathogen transcriptome

    A: GO enrichment analysis of pathogen DEGs. B: KEGG pathway enrichment analysis of pathogen DEGs. C: Heat map of DEGs for top 30 significantly enriched GO terms.

    Figure  5.  DEG enrichment analysis of L. crocea transcriptome data

    A: Volcano plot obtained from edgeR analysis of infected L. crocea RNA pools. B: Heat map analysis of spleen transcript DEGs in infected L. crocea (FDR<0.05, |log2FC|≥1), C: Five up-regulated and down-regulated DEGs of host were randomly selected for qRT-PCR verification.

    Figure  6.  GO and KEGG enrichment analysis of host transcriptome

    A: GO enrichment analysis of pathogen DEGs. B: Top 10 GO term-gene annotation enrichment analysis results. C: KEGG pathway enrichment analysis of host DEGs. D: DEGs of cytokine-cytokine receptor interaction (Ko04060), hematopoietic cell lineage (Ko04640). and Toll-like receptor signaling pathway (Ko04620).

    Figure  7.  Schematic overview of response of TNF signaling pathway of L. crocea to P. plecoglossicida infection

    Color indicates relative expression level of gene (RNAi strain/wild-type strain), red is up-regulated, blue is down-regulated. Darker color indicates greater change in gene expression.

  • [1] Alcamo E, Mizgerd JP, Horwitz BH, Bronson R, Beg AA, Scott M, Doerschuk CM, Hynes RO, Baltimore D. 2001. Targeted mutation of TNF receptor I rescues the RelA-deficient mouse and reveals a critical role for NF-κB in leukocyte recruitment. The Journal of Immunology, 167(3): 1592−1600. doi:  10.4049/jimmunol.167.3.1592
    [2] Anders S, Huber W. 2010. Differential expression analysis for sequence count data. Genome Biology, 11(10): R106. doi:  10.1186/gb-2010-11-10-r106
    [3] Chen JP, Pang W, Zhao ZW, Bi YH, Chen XW. 2019. Transcription profiles of skin and head kidney from goldfish suffering hemorrhagic septicemia with an emphasis on the TLR signaling pathway. Zoological Research, 40(4): 337−342. doi:  10.24272/j.issn.2095-8137.2019.028
    [4] Collette Y, Gilles A, Pontarotti P, Olive D. 2003. A co-evolution perspective of the TNFSF and TNFRSF families in the immune system. Trends in Immunology, 24(7): 387−394. doi:  10.1016/S1471-4906(03)00166-2
    [5] Crofts AA, Giovanetti SM, Rubin EJ, Poly FM, Gutierrez RL, Talaat KR, Porter CK, Riddle MS, DeNearing B, Brubaker J, Maciel M Jr, Alcala AN, Chakraborty S, Prouty MG, Savarino SJ, Davies BW, Trent MS. 2018. Enterotoxigenic E. coli virulence gene regulation in human infections. Proceedings of the National Academy of Sciences of the United States of America, 115(38): E8968−E8976. doi:  10.1073/pnas.1808982115
    [6] Dalebroux ZD, Yagi BF, Sahr T, Buchrieser C, Swanson MS. 2010. Distinct roles of ppGpp and DksA in Legionella pneumophila differentiation. Molecular Microbiology, 76(1): 200−219. doi:  10.1111/j.1365-2958.2010.07094.x
    [7] Delves PJ. 2020. Innate and adaptive systems of immunity. In: Rose NR, Mackay IR. The Autoimmune Diseases. 6th ed. San Diego: Academic Press, 45–61.
    [8] Gardell SJ, Craik CS, Clauser E, Goldsmith EJ, Stewart CB, Graf M, Rutter WJ. 1988. A novel rat carboxypeptidase, CPA2: characterization, molecular cloning, and evolutionary implications on substrate specificity in the carboxypeptidase gene family. Journal of Biological Chemistry, 263(33): 17828−17836.
    [9] Godahewa GI, Wickramaarachchi WDN, Whang I, Bathige SDNK, Lim BS, Choi CY, De Zoysa M, Noh JK, Lee J. 2014. Two carboxypeptidase counterparts from rock bream (Oplegnathus fasciatus): molecular characterization, genomic arrangement and immune responses upon pathogenic stresses. Veterinary Immunology and Immunopathology, 162(3-4): 180−191. doi:  10.1016/j.vetimm.2014.10.005
    [10] Guanzon DAV, Maningas MBB. 2018. Functional elucidation of LvToll 3 receptor from P. vannamei through RNA interference and its potential role in the shrimp antiviral response. Developmental & Comparative Immunology, 84: 172−180.
    [11] Guo LN, Huang LX, Su YQ, Qin YX, Zhao LM, Yan QP. 2018. secA, secD, secF, yajC, and yidC contribute to the adhesion regulation of Vibrio alginolyticus. Microbiologyopen, 7(2): e00551. doi:  10.1002/mbo3.551
    [12] He Y, Xu TT, Fossheim LE, Zhang XH. 2012. FliC, a flagellin protein, is essential for the growth and virulence of fish pathogen Edwardsiella tarda. PLoS One, 7(9): e45070. doi:  10.1371/journal.pone.0045070
    [13] Huang LX, Liu WJ, Jiang QL, Zuo YF, Su YQ, Zhao LM, Qin YX, Yan QP. 2018. Integration of Transcriptomic and Proteomic Approaches Reveals the Temperature-Dependent Virulence of Pseudomonas plecoglossicida. Frontiers in Cellular and Infection Microbiology, 8: 207. doi:  10.3389/fcimb.2018.00207
    [14] Huang LX, Zhao LM, Qi WL, Xu XJ, Zhang JN, Zhang JL, Yan QP. 2020. Temperature-specific expression of cspA1 contributes to activation of sigX during pathogenesis and intracellular survival in Pseudomonas plecoglossicida. Aquaculture, 518: 734861. doi:  10.1016/j.aquaculture.2019.734861
    [15] Huang LX, Zuo YF, Jiang QL, Su YQ, Qin YX, Xu XJ, Zhao LM, Yan QP. 2019. A metabolomic investigation into the temperature-dependent virulence of Pseudomonas plecoglossicida from large yellow croaker (Pseudosciaena crocea). Journal of Fish Diseases, 42(3): 431−446. doi:  10.1111/jfd.12957
    [16] Joosten LAB, Abdollahi-Roodsaz S, Dinarello CA, O'Neill L, Netea MG. 2016. Toll-like receptors and chronic inflammation in rheumatic diseases: new developments. Nature Reviews Rheumatology, 12(6): 344−357. doi:  10.1038/nrrheum.2016.61
    [17] Kamarthapu V, Epshtein V, Benjamin B, Proshkin S, Mironov A, Cashel M, Nudler E. 2016. ppGpp couples transcription to DNA repair in E. coli. Science, 352(6288): 993−996. doi:  10.1126/science.aad6945
    [18] Kasai K, Ishikawa T, Komata T, Fukuchi K, Chiba M, Nozaka H, Nakamura T, Sato T, Miura T. 2010. Novel L-amino acid oxidase with antibacterial activity against methicillin-resistant Staphylococcus aureus isolated from epidermal mucus of the flounder Platichthys stellatus. The FEBS Journal, 277(2): 453−465. doi:  10.1111/j.1742-4658.2009.07497.x
    [19] Kim D, Langmead B, Salzberg SL. 2015. HISAT: a fast spliced aligner with low memory requirements. Nature Methods, 12(4): 357−360. doi:  10.1038/nmeth.3317
    [20] Langmead B, Salzberg SL. 2012. Fast gapped-read alignment with Bowtie 2. Nature Methods, 9(4): 357−359. doi:  10.1038/nmeth.1923
    [21] Lemke JJ, Durfee T, Gourse RL. 2009. DksA and ppGpp directly regulate transcription of the Escherichia coli flagellar cascade. Molecular Microbiology, 74(6): 1368−1379. doi:  10.1111/j.1365-2958.2009.06939.x
    [22] Li CH, Xiong JB, Ding FF, Chen J. 2020. Immune and gut bacterial successions of large yellow croaker (Larimichthys crocea) during Pseudomonas plecoglossicida infection. Fish & Shellfish Immunology, 99: 176−183.
    [23] Liu WJ, Huang LX, Su YQ, QinYX, Zhao LM, Y an, Q P. 2017. Contributions of the oligopeptide permeases in multistep of vibrio alginolyticus pathogenesis. Microbiologyopen, 6(1): e00511.
    [24] Lu XJ, Chen J. 2019. Specific function and modulation of teleost monocytes/macrophages: polarization and phagocytosis. Zoological Research, 40(3): 146−150. doi:  10.24272/j.issn.2095-8137.2019.035
    [25] Luo G, Sun YJ, Huang LX, Su YQ, Zhao LM, Qin YX, Xu XJ, Yan QP. 2020. Time-resolved dual RNA-seq of tissue uncovers Pseudomonas plecoglossicida key virulence genes in host-pathogen interaction with Epinephelus coioides. Environmental Microbiology, 22(2): 677−693. doi:  10.1111/1462-2920.14884
    [26] Luo G, Xu XJ, Zhao LM, Qin YX, Huang LX, Su YQ, Yan QP. 2019b. clpV is a key virulence gene during in vivo Pseudomonas plecoglossicida infection. Journal of Fish Diseases, 42(7): 991−1000.
    [27] Luo G, Zhao LM, Xu XJ, Qin YX, Huang LX, Su YQ, Zheng WQ, Yan QP. 2019a. Integrated dual RNA-seq and dual iTRAQ of infected tissue reveals the functions of a diguanylate cyclase gene of Pseudomonas plecoglossicida in host-pathogen interactions with Epinephelus coioides. Fish & shellfish Immunology, 95: 481−490.
    [28] Mallik P, Paul BJ, Rutherford ST, Gourse RL, Osuna R. 2006. DksA is required for growth phase-dependent regulation, growth rate-dependent control, and stringent control of fis expression in Escherichia coli. Journal of Bacteriology, 188(16): 5775−5782. doi:  10.1128/JB.00276-06
    [29] Mandel M, Higa A. 1970. Calcium-dependent bacteriophage DNA infection. Journal of Molecular Biology, 53(1): 159−162. doi:  10.1016/0022-2836(70)90051-3
    [30] Mantovani A, Dinarello CA, Molgora M, Garlanda C. 2019. Interleukin-1 and Related Cytokines in the Regulation of Inflammation and Immunity. Immunity, 50(4): 778−795. doi:  10.1016/j.immuni.2019.03.012
    [31] Marshall JS, Warrington R, Watson W, Kim HL. 2018. An introduction to immunology and immunopathology. Allergy, Asthma & Clinical Immunology, 14(S2): 49.
    [32] Nishimori E, Kita-Tsukamoto K, Wakabayashi H. 2000. Pseudomonas plecoglossicida sp. nov., the causative agent of bacterial haemorrhagic ascites of ayu, Plecoglossus altivelis. International Journal of Systematic and Evolutionary Microbiology, 50(1): 83−89. doi:  10.1099/00207713-50-1-83
    [33] Nuss AM, Beckstette M, Pimenova M, Schmühl C, Opitz W, Pisano F, Heroven AK, Dersch P. 2017. Tissue dual RNA-seq allows fast discovery of infection-specific functions and riboregulators shaping host-pathogen transcriptomes. Proceedings of the National Academy of Sciences of the United States of America, 114(5): E791−E800. doi:  10.1073/pnas.1613405114
    [34] Palti Y. 2011. Toll-like receptors in bony fish: from genomics to function. Developmental & Comparative Immunology, 35(12): 1263−1272.
    [35] Paul BJ, Barker MM, Ross W, Schneider DA, Webb C, Foster JW, Gourse RL. 2004. DksA: a critical component of the transcription initiation machinery that potentiates the regulation of rRNA promoters by ppGpp and the initiating NTP. Cell, 118(3): 311−322. doi:  10.1016/j.cell.2004.07.009
    [36] Robinson MD, McCarthy DJ, Smyth GK. 2010. edgeR: a Bioconductor package for differential expression analysis of digital gene expression data. Bioinformatics, 26(1): 139−140. doi:  10.1093/bioinformatics/btp616
    [37] Rong DL, Wu QP, Xu MF, Zhang JM, Yu SB. 2017. Prevalence, virulence genes, antimicrobial susceptibility, and genetic diversity of Staphylococcus aureus from retail aquatic products in China. Frontiers in Microbiology, 8: 714. doi:  10.3389/fmicb.2017.00714
    [38] Saleh M, Kumar G, Abdel-Baki AA, Dkhil MA, El-Matbouli M, Al-Quraishy S. 2016. In vitro gene silencing of the fish microsporidian Heterosporis saurida by RNA interference. Nucleic Acid Therapeutics, 26(4): 250−256. doi:  10.1089/nat.2016.0613
    [39] Sharma A, Steichen AL, Jondle CN, Mishra BB, Sharma J. 2014. Protective role of Mincle in bacterial pneumonia by regulation of neutrophil mediated phagocytosis and extracellular trap formation. The Journal of Infectious Disease, 209(11): 1837−1846. doi:  10.1093/infdis/jit820
    [40] Sharma AK, Payne SM. 2006. Induction of expression of hfq by DksA is essential for Shigella flexneri virulence. Molecular Microbiology, 62(2): 469−479. doi:  10.1111/j.1365-2958.2006.05376.x
    [41] Stokes BA, Yadav S, Shokal U, Smith LC, Eleftherianos I. 2015. Bacterial and fungal pattern recognition receptors in homologous innate signaling pathways of insects and mammals. Frontiers in Microbiology, 6: 19.
    [42] Sun YJ, Luo G, Zhao LM, Huang LX, Qin YX, Su YQ, Yan QP. 2018. Integration of RNAi and RNA-seq Reveals the Immune Responses of Epinephelus coioides to sigX Gene of Pseudomonas plecoglossicida. Frontiers in Immunology, 9: 1624. doi:  10.3389/fimmu.2018.01624
    [43] Sun YJ, Nie P, Zhao LM, Huang LX, Qin YX, Xu XJ, Zhang JN, Yan QP. 2019a. Dual RNA-Seq unveils the role of the Pseudomonas plecoglossicida fliA gene in pathogen-host interaction with Larimichthys crocea. Microorganisms, 7(10): 443. doi:  10.3390/microorganisms7100443
    [44] Sun YJ, Zhuang ZX, Wang XR, Huang HB, Fu Q, Yan QP. 2019b. Dual RNA-seq reveals the effect of the flgM gene of Pseudomonas plecoglossicida on the immune response of Epinephelus coioides. Fish & Shellfish Immunology, 87: 515−523.
    [45] Takeda K, Akira S. 2015. Toll-like receptors. Current Protocols in Immunology, 109(14.12.1): 14.12.10.
    [46] Tan SY, Chowdhury S, Polak N, Gorrell MD, Weninger W. 2017. Fibroblast activation protein is dispensable in the anti-influenza immune response in mice. PLoS One, 12(2): e0171194. doi:  10.1371/journal.pone.0171194
    [47] Tang RQ, Luo G, Zhao LM, Huang LX, Qin YX, Xu XJ, Su YQ, Yan QP. 2019b. The effect of a LysR-type transcriptional regulator gene of Pseudomonas plecoglossicida on the immune responses of Epinephelus coioides. Fish & Shellfish Immunology, 89: 420−427.
    [48] Tang RQ, Zhao LM, Xu XJ, Huang LX, Qin YX, Su YQ, Yan QP. 2019a. Dual RNA-Seq uncovers the function of an ABC transporter gene in the host-pathogen interaction between Epinephelus coioides and Pseudomonas plecoglossicida. Fish & Shellfish Immunology, 92: 45−53.
    [49] Tang Y, Sun YJ, Zhao LM, Xu XJ, Huang LX, Qin YX, Su YQ, Yi GF, Yan QP. 2019c. Mechanistic insight into the roles of Pseudomonas plecoglossicida clpV gene in host-pathogen interactions with Larimichthys crocea by dual RNA-seq. Fish & Shellfish Immunology, 93: 344−353.
    [50] Tang Y, Xin G, Zhao LM, Huang LX, Qin YX, Su YQ, Zheng WQ, Wu B, Lin N, Yan QP. 2020. Novel insights into host-pathogen interactions of large yellow croakers (Larimichthys crocea) and pathogenic bacterium Pseudomonas plecoglossicida using time-resolved dual RNA-seq of infected spleens. Zoological Research, 41(3): 341−327.
    [51] Tao Z, Xu Y, Zhou SM, Qian D, Liu MH, Li WY, Xu WJ, Yan XJ. 2020. Acquisition of a type VI secretion system is critical for Pseudomonas plecoglossicida induced granulomas in fish internal organs. Aquaculture, 516: 734629. doi:  10.1016/j.aquaculture.2019.734629
    [52] Tao Z, Zhou T, Zhou SM, Wang GL. 2016. Temperature-regulated expression of type VI secretion systems in fish pathogen Pseudomonas plecoglossicida revealed by comparative secretome analysis. FEMS Microbiology Letters, 363(22): fnw261. doi:  10.1093/femsle/fnw261
    [53] Valenzuela-Miranda D, Gallardo-Escarate C. 2018. Dual RNA-seq uncovers metabolic amino acids dependency of the intracellular bacterium Piscirickettsia salmonis infecting Atlantic salmon. Frontiers in Microbiology, 9: 2877. doi:  10.3389/fmicb.2018.02877
    [54] Wang LY, Sun YJ, Zhao LM, Xu XJ, Huang LX, Qin YX, Su YQ, Zhang JN, Yan QP. 2019a. Dual RNA-seq uncovers the immune response of Larimichthys crocea to the secY gene of Pseudomonas plecoglossicida from the perspective of host-pathogen interactions. Fish & Shellfish Immunology, 93: 949−957.
    [55] Wang YQ, Wang XY, Ali F, Li ZQ, Fu YY, Yang XJ, Lin WX, Lin XM. 2019b. Comparative extracellular proteomics of Aeromonas hydrophila reveals iron-regulated secreted proteins as potential vaccine candidates. Frontiers in Immunology, 10: 256. doi:  10.3389/fimmu.2019.00256
    [56] Westermann AJ, Barquist L, Vogel J. 2017. Resolving host-pathogen interactions by dual RNA-seq. PLoS Pathogens, 13(2): e1006033. doi:  10.1371/journal.ppat.1006033
    [57] Westermann AJ, Förstner KU, Amman F, Barquist L, Chao YJ, Schulte LN, Müller L, Reinhardt R, Stadler PF, Vogel J. 2016. Dual RNA-seq unveils noncoding RNA functions in host-pathogen interactions. Nature, 529(7587): 496−501. doi:  10.1038/nature16547
    [58] Westermann AJ, Gorski SA, Vogel J. 2012. Dual RNA-seq of pathogen and host. Nature Reviews Microbiology, 10(9): 618−630. doi:  10.1038/nrmicro2852
    [59] Xie C, Mao XZ, Huang JJ, Ding Y, Wu JM, Dong S, Kong L, Gao G, Li CY, Wei LP. 2011. KOBAS 2.0: a web server for annotation and identification of enriched pathways and diseases. Nucleic Acids Research, 39(S2): W316−W322.
    [60] Yang J, Lu XJ, Chai FC, Chen J. 2016. Molecular characterization and functional analysis of a piscidin gene in large yellow croaker (Larimichthys crocea). Zoological Research, 37(6): 347−355.
    [61] Yang Y, Adeola AC, Xie HB, Zhang YP. 2018. Genomic and transcriptomic analyses reveal selection of genes for puberty in Bama Xiang pigs. Zoological Research, 39(6): 424−430. doi:  10.24272/j.issn.2095-8137.2018.068
    [62] Yao ZJ, Guo Z, Wang YQ, Li WX, Fu YY, Lin YX, Lin WX, Lin XM. 2019. Integrated succinylome and metabolome profiling reveals crucial role of S-ribosylhomocysteine lyase in quorum sensing and metabolism of Aeromonas hydrophila. Molecular & Cellular Proteomics, 18(2): 200−215.
    [63] Ye C, An X, Jiang YD, Ding BY, Shang F, Christiaens O, Taning CNT, Smagghe G, Niu JZ, Wang JJ. 2018. Induction of RNAi core machinery's gene expression by exogenous dsRNA and the effects of pre-exposure to dsRNA on the gene silencing efficiency in the pea aphid (Acyrthosiphon pisum). Frontiers in Physiology, 9: 1906.
    [64] Zhang BB, Luo G, Zhao LM, Huang LX, Qin YX, Su YQ, Yan QP. 2018a. Integration of RNAi and RNA-seq uncovers the immune responses of Epinephelus coioides to L321_RS19110 gene of Pseudomonas plecoglossicida. Fish & Shellfish Immunology, 81: 121−129.
    [65] Zhang BB, Zhuang ZX, Wang XR, Huang HB, Fu Q, Yan QP. 2019b. Dual RNA-Seq reveals the role of a transcriptional regulator gene in pathogen-host interactions between Pseudomonas plecoglossicida and Epinephelus coioides. Fish & Shellfish Immunology, 87: 778−787.
    [66] Zhang J, Wang Y, Guo H, Mao Z, Ge C. 2017. Identification and characterization of a phospholipase A1 activity type three secreted protein, PP_ExoU from Pseudomonas plecoglossicida NB2011, the causative agent of visceral granulomas disease in large yellow croaker (Larimichthys crocea). Journal of Fish Diseases, 40(6): 831−840. doi:  10.1111/jfd.12565
    [67] Zhang JT, Zhou SM, An SW, Chen L, Wang GL. 2014. Visceral granulomas in farmed large yellow croaker, Larimichthys crocea (Richardson), caused by a bacterial pathogen, Pseudomonas plecoglossicida. Journal of Fish Diseases, 37(2): 113−121. doi:  10.1111/jfd.12075
    [68] Zhang MM, Qin YX, Huang LX, Yan QP, Mao LL, Xu XJ, Wang SY, Zhang MM, Chen LW. 2019a. The role of sodA and sodB in Aeromonas hydrophila resisting oxidative damage to survive in fish macrophages and escape for further infection. Fish & Shellfish Immunology, 88: 489−495.
    [69] Zhang MM, Yan QP, Mao LL, Wang SY, Huang LX, Xu XJ, Qin YX. 2018b. KatG plays an important role in Aeromonas hydrophila survival in fish macrophages and escape for further infection. Gene, 672: 156−164. doi:  10.1016/j.gene.2018.06.029
    [70] Zuo YF, Zhao LM, Xu XJ, Zhang JN, Zhang JL, Yan QP, Huang LX. 2019. Mechanisms underlying the virulence regulation of new Vibrio alginolyticus ncRNA Vvrr1 with a comparative proteomic analysis. Emerging Microbes & Infections, 8(1): 1604−1618.
  • [1] Xi-Yao Ma, Tiao Ning, Adeniyi C. Adeola, Jie Li, Ali Esmailizadeh, Jacqueline K. Lichoti, Bernard R. Agwanda, Jainagul Isakova, Almaz A. Aldashev, Shi-Fang Wu, He-Qun Liu, Najmudinov Tojiddin Abdulloevich, Manilova Elena Afanasevna, Khudoidodov Behruz Ibrohimovich, Rahamon Akinyele Moshood Adedokun, Sunday Charles Olaogun, Oscar J. Sanke, Godwin F. Mangbon, Xi Chen, Wei-Kang Yang, Zhe Wang, Min-Sheng Peng, Sheila C. Ommeh, Yan Li, Ya-Ping Zhang.  Potential dual expansion of domesticated donkeys revealed by worldwide analysis on mitochondrial sequences, Zoological Research. doi: 10.24272/j.issn.2095-8137.2020.007
    [2] Yi Tang, Ge Xin, Ling-Min Zhao, Li-Xing Huang, Ying-Xue Qin, Yong-Quan Su, Wei-Qiang Zheng, Bin Wu, Nan Lin, Qing-Pi Yan.  Novel insights into host-pathogen interactions of large yellow croakers (Larimichthys crocea) and pathogenic bacterium Pseudomonas plecoglossicida using time-resolved dual RNA-seq of infected spleens, Zoological Research. doi: 10.24272/j.issn.2095-8137.2020.035
    [3] Chang-Zhang Feng, Can-Chao Yang, Wei Liang.  Nest sanitation facilitates egg recognition in the common tailorbird, a plaintive cuckoo host, Zoological Research. doi: 10.24272/j.issn.2095-8137.2019.054
    [4] Zheng-Bo Wang, Dong-Dong Qin, Xin-Tian Hu.  Engrafted newborn neurons could functionally integrate into the host neuronal network, Zoological Research. doi: 10.13918/j.issn.2095-8137.2017.005
    [5] Jing YANG, Xin-Jiang LU, Fang-Chao CHAI, Jiong CHEN.  Molecular characterization and functional analysis of a piscidin gene in large yellow croaker (Larimichthys crocea), Zoological Research. doi: 10.13918/j.issn.2095-8137.2016.6.347
    [6] Feng-Chao LANG, Xin LI, Olga VLADMIROVA, Zhuo-Ran LI, Gui-Jun CHEN, Yu XIAO, Li-Hong LI, Dan-Feng LU, Hong-Bo HAN, Ju-Min ZHOU.  Selective recruitment of host factors by HSV-1 replication centers, Zoological Research.
    [7] Jun WU, Yu-Hong SHI, Xue-Heng ZHANG, Chang-Hong LI, Ming-Yun LI, Jiong CHEN.  Molecular characterization of an IL-1β gene from the large yellow croaker (Larimichthys crocea) and its effect on fish defense against Vibrio alginolyticus infection, Zoological Research.
    [8] Guang-Liang ZHANG, Zhan-Hui TANG, Ti-Yu HONG, Jian YANG, Yu ZENG, Liang-Jing TAN, Qi LIU, Yi CHEN, Qi-Qi SHEN, Li-Biao ZHANG.  Host selection of ectoparasitic gamasid mites on Tylonycteris pachypus and T. robustula, Zoological Research. doi: 10.3724/SP.J.1141.2013.01021
    [9] Yanfeng ZHANG, Bing SU.  Peak identification for ChIP-seq data with no controls, Zoological Research. doi: 10.3724/SP.J.1141.2012.E05-06E121
    [10] CHANG Yu-mei, XU Wan-tu, CHI Bing-jie, GAO Guo-qiang, HAN Qi-xia, HE Wei, SUN Xiao-we.  Kinship Analysis of One Broodstock Population of Large Yellow Croaker Pseudosciaena crocea, Zoological Research. doi: 10.3724/SP.J.1141.2009.06620
    [11] CUI Shu-fang, QU Xiu-hua, HAN Shi-zhong, YU Chen-lin, SUN Wei, TANG Qiu.  Application of Dual-color Fluorescence Hybridization Chip Technique in the Monitoring of the Genetic Quality of Inbred Mice Strains, Zoological Research.
    [12] ZHANG Hui-Jie, WANG Hui-Jin, LI Jian-She, ZHANG Li-Ping, DONG Zhe-Sheng, XU Qi.  The Relationship between the Fitness of Host Plants to Liriomyza sativae and the Structure of Leaf Blades, Zoological Research.
    [13] YOU Yong-Long, LIN Dan-Jun, CHEN Lian-Yun.  Spermatogenesis of Teleosts,Pseudosciaena crocea, Zoological Research.
    [14] LIN Hao-ran.  The Interactions of Neuroendocrine Regulation on Reproduction and Growth in Fish, Zoological Research.
    [15] ZHANG Jian-xu, WANG Zu-wang, ZHANG Zhi-bin.  Odor Preference and Behavioral Interactions of Ratlike Hamsters (Cricetulus triton) During The Nonbreeding Season, Zoological Research.
    [16] ZHANG Jing, GU Bao-hong, PENG Shou-li, LIU Ci-quan.  Distributions of Triplet Codons in Messenger RNA Secondary Structures, Zoological Research.
    [17] DUAN Zhi-zhong, TU Guang-chou.  Properties of Phosphodiesterase of Agkistrodon Halys (Pallas) Venom (From Zhejiang Province) and its Actions on Both Rna and DNA, Zoological Research.
    [18] SHI Fu-tian, ZHANG Zo-run.  The Community of Intestinal Flagellates and its Significance in the Classification of Their Host Termites, Zoological Research.
    [19] LI Jing-yan, CHEN Yun-he, ZHOU Xing-liang, ZHENG Zi-xiu.  Immunoprotective Activity of Ribosomal Preparation From Poisonous Escherichia coli,the Pathogen of Piglet Dysentery, Zoological Research.
    [20] XU Yi-sheng, PHILIPPE Brulet.  Cytoplasmic Distribution of A Mouse Embryonal Carcinoma Cell Specific RNA,MS3 RNA, Zoological Research.
  • ZR-2020-034 Supporting Materials.doc
  • 加载中
图(7)
计量
  • 文章访问数:  867
  • HTML全文浏览量:  262
  • PDF下载量:  38
  • 被引次数: 0
出版历程
  • 收稿日期:  2020-02-20
  • 录用日期:  2020-04-21
  • 网络出版日期:  2020-05-13
  • 刊出日期:  2020-07-01

Dual RNA-seq provides novel insight into the roles of dksA from Pseudomonas plecoglossicida in pathogen-host interactions with large yellow croakers (Larimichthys crocea)

doi: 10.24272/j.issn.2095-8137.2020.048
    基金项目:  This work was supported by the National Natural Science Foundation of China (31672694), Fujian Provincial Special Fund for Marine and Fishery Protection and Development (MCZ[2019]062), and Open Fund of Fujian Province Key Laboratory of Special Aquatic Formula Feed (2019-01)
    通讯作者: E-mail: yanqp@jmu.edu.cn

摘要: 变形假单胞菌(Pseudomonas plecoglossicida)是革兰氏阴性菌,呈杆状,有鞭毛,该菌在大黄鱼(Larimichthys crocea)养殖过程中引起内脏白点病并有非常高的死亡率,可造成严重的经济损失。经互作转录组数据分析和实时定量PCR验证,dksA基因在变形假单胞菌感染斜带石斑鱼(Epinephelus coioides)48小时后的表达量上调显著(log2FC=3.12 ,P<0.001)。利用RNAi技术显著降低了变形假单胞菌中基因dksA的表达量,最佳沉默效率达到96.23%。分别用变形假单胞菌野生株与效果最好的沉默菌株注射感染健康大黄鱼,观察大黄鱼的存活情况,并于注射后48小时取脾脏进行dual RNA-seq。与野生型毒株相比,感染沉默菌株的大黄鱼内脏白点病症状减轻,死亡时间延迟了2 天且死亡率降低20%。经dual RNA-seq 得知人工感染后病原和宿主的转录组发生了显著变化。dksA基因的沉默导致病原菌的细胞组分,鞭毛和核糖体组装相关基因大量下调,而fliC基因的显著上调可能是变形假单胞菌毒力维持的一种途径。GO和KEGG结果表明沉默菌株感染的大黄鱼的多个与免疫相关通路中炎性因子基因大量下调,在这些免疫途径中,受影响的基因与多个免疫应答过程相关。研究结果表明,dksA是变形假单胞菌的毒力基因,相比于野生株感染,沉默菌株感染的大黄鱼免疫反应较弱。

English Abstract

王路英, 刘子旭, 赵玲敏, 黄力行, 覃映雪, 苏永全, 郑炜强, 王凡, 鄢庆枇. 基于Dual RNA-seq研究dksA基因在变形假单胞菌与大黄鱼的病原体-宿主相互作用中的功能[J]. 动物学研究, 2020, 41(4): 410-422. doi: 10.24272/j.issn.2095-8137.2020.048
引用本文: 王路英, 刘子旭, 赵玲敏, 黄力行, 覃映雪, 苏永全, 郑炜强, 王凡, 鄢庆枇. 基于Dual RNA-seq研究dksA基因在变形假单胞菌与大黄鱼的病原体-宿主相互作用中的功能[J]. 动物学研究, 2020, 41(4): 410-422. doi: 10.24272/j.issn.2095-8137.2020.048
Lu-Ying Wang, Zi-Xu Liu, Ling-Min Zhao, Li-Xing Huang, Ying-Xue Qin, Yong-Quan Su, Wei-Qiang Zheng, Fan Wang, Qing-Pi Yan. Dual RNA-seq provides novel insight into the roles of dksA from Pseudomonas plecoglossicida in pathogen-host interactions with large yellow croakers (Larimichthys crocea). Zoological Research, 2020, 41(4): 410-422. doi: 10.24272/j.issn.2095-8137.2020.048
Citation: Lu-Ying Wang, Zi-Xu Liu, Ling-Min Zhao, Li-Xing Huang, Ying-Xue Qin, Yong-Quan Su, Wei-Qiang Zheng, Fan Wang, Qing-Pi Yan. Dual RNA-seq provides novel insight into the roles of dksA from Pseudomonas plecoglossicida in pathogen-host interactions with large yellow croakers (Larimichthys crocea). Zoological Research, 2020, 41(4): 410-422. doi: 10.24272/j.issn.2095-8137.2020.048
    • Infection is an exceedingly complex process involving strong interactions between pathogen and host (Luo et al., 2020). Both pathogen and host must mobilize all available resources to win this life and death battle, and many changes during infection are reflected in their respective transcripts (Zhang et al., 2018a). Therefore, simultaneous detection of transcriptome profiles during infection can provide insight into the pathogenic mechanisms and host immune responses (Tang et al., 2019a). For a long time, due to technical limitations, studies on infection have focused on either the host or the pathogen (Sun et al., 2018). The advancement of dual RNA-seq, which can simultaneously detect both pathogen and host transcriptomes, has provided a powerful and advantageous tool for studying various infection models and pathogen-host interactions (Westermann et al., 2012, 2016, 2017; Valenzuela-Miranda & Gallardo-Escarate, 2018). Dual RNA-seq combined with RNAi has also been used to study the role of virulence genes in host-pathogen interactions (Nuss et al., 2017; Wang et al., 2019a). In addition, dual RNA-seq and dual iTRAQ have been applied to explore gene functions at the multi-omics level (Luo et al., 2019a).

      Large yellow croakers (Larimichthys crocea), an economically important marine fish in China, are widely cultured in the Fujian and Zhejiang provinces (Yang et al., 2016). The most important factor threatening L. crocea culture is the frequent occurrence of infection epidemics (Tang et al., 2019c). Visceral white spot disease, one of the most destructive diseases in L. crocea, is caused by Pseudomonas plecoglossicida, a short rod-shaped gram-negative bacterium (Huang et al., 2018; Zhang et al., 2014). Pseudomonas plecoglossicida was first isolated from ayu fish (Plecoglossus altivelis) suffering from bacterial hemorrhagic ascites (Nishimori et al., 2000). In view of its potential to cause great harm to aquaculture, the pathogenic mechanism of P. plecoglossicida has attracted considerable attention (Huang et al., 2019; Tao et al., 2016) and several virulence genes have been identified (Tao et al., 2020; Zhang et al., 2017). Earlier studies reported that P. plecoglossicida commonly circulates in orange-spotted groupers (Epinephelus coioides) and can easily adapt to and proliferate in the splenic environment, resulting in significantly higher pathogen loads in spleens (>200 times) than in other tissues (Luo et al., 2020). In addition, P. plecoglossicida infection can cause irreversible disruption of gut microbiota, resulting in increasing mortality (Li et al., 2020). Furthermore, CspA1 is known to contribute to P. plecoglossicida virulence in a temperature-specific manner via regulation of sigX expression (Huang et al., 2020).

      In previous studies from our laboratory, we found that the dksA gene from P. plecoglossicida was highly expressed during host infection (data deposited in GenBank SRA database under accession numbers SRP114910 and SRP115064). The dksA gene encodes an RNA polymerase-binding transcription factor (Kamarthapu et al., 2016). The gene has a variety of important functions, such as regulating rRNA promoter activity (Paul et al., 2004), HFQ gene expression and virulence of Shigella flexneri (Sharma & Payne, 2006), growth rate of Escherichia coli (Mallik et al., 2006), and assembly of flagella (Dalebroux et al., 2010). Thus, based on our previous laboratory research and reported literature, dksA may play a role in the pathogenicity of P. plecoglossicida. To date, however, there are no reports on the function of dksA in pathogen-host interactions.

      To explore the roles of dksA in host-pathogen interactions between P. plecoglossicida and L. crocea, a dksA-silenced strain of P. plecoglossicida was constructed by RNAi technology and differences in virulence between wild-type and RNAi strains were analysed. The spleens of L. crocea infected by the wild-type or dksA-RNAi strains of P. plecoglossicida were subjected to dual RNA-seq. The present study provides novel insight into host-pathogen interactions between P. plecoglossicida and L. crocea.

    • The highly pathogenic wild-type strain of P. plecoglossicida (NZBD9) was isolated from the spleen of a diseased large yellow croaker (Huang et al., 2018). Pseudomonas plecoglossicida was shake-cultured (220 r/min) in Luria Bertani (LB) broth at 18 or 28 °C. In addition, Escherichia coli DH5α was obtained from the Beijing Tiangen Company (China) and cultured in LB broth at 37 °C and 220 r/min.

    • The RNAi strains were constructed according to Sun et al. (2018). Five short hairpin RNA sequences targeting the dksA gene were designed using the RNAi website (http://rnaidesigner.thermofisher.com/rnaiexpress/setOption.do?designOption=shrna&pid=708587103220684543), and then synthesized by Shanghai Generay Biotech Co., Ltd. (China) (Supplementary Table S1). Each oligonucleotide was annealed and ligated to the pCM130/tac vector linearized with the restriction enzymes NsiI and BsrGI (New England Biolabs, USA) using T4 DNA ligase (New England Biolabs, USA) (Guo et al., 2018). The preparation of E. coli DH5α competent cells was performed by the CaCl2 method (Mandel & Higa, 1970). The recombinant pCM130/tac vector was transformed into competent E. coli DH5α by heat shock, and then extracted for electroporating into P. plecoglossicida competent cells (Tang et al., 2019b). Finally, the expression of dksA in five dksA-RNAi strains of P. plecoglossicida was verified by quantitative real-time polymerase chain reaction (qRT-PCR).

    • All fish experiments were carried out strictly following the “Guide for the Care and Use of Laboratory Animals” established by the National Institutes of Health. All animal protocols were approved by the Jimei University Animal Ethics Committee (Acceptance No. JMULAC201159).

      In total, 200 healthy L. crocea (average weight ~50 g) were obtained from Ningde (Fujian, China) and acclimatized to laboratory conditions for one week. To detect the pathogenicity of P. plecoglossicida in L. crocea, 60 fish were randomly divided into three groups. Each acclimatized fish was intrapleurally injected with 104 colony forming units per gram fish (cfu/g) of the wild-type or RNAi strain of P. plecoglossicida, respectively. As the negative control, other L. crocea individuals were intrapleurally injected with 0.1 mL of phosphate-buffered saline (PBS). The water temperature throughout the experiment was maintained at 18±1 °C. The status of injected fish was recorded twice a day. For spleen sampling, six spleens from wild-type or RNAi strain-infected L. crocea were randomly sampled at 48 h post infection (hpi) for dual RNA-seq, with two spleens mixed as one sample. In addition, six spleens from L. crocea infected by the wild-type or RNAi strain were randomly sampled at 6, 12, 24, 48, 72, and 96 hpi for pathogen load and dksA expression assays.

    • DNA purification of spleen samples was performed according to the instructions provided with the EasyPure Marine Animal Genomic DNA Kit (TransGen Biotech, China). The extracted genomic DNA was stored at –20 °C until use.

      Total RNA was extracted using an Eastep® Super Total RNA Extraction Kit (Shanghai Promega Biological Products, Ltd., China). The quality of total RNA was checked by agarose gel electrophoresis. cDNA was synthesized by TransScript All-in-One First-Strand cDNA Synthesis SuperMix for qPCR (One-Step gDNA Removal) (TransGen Biotech, China) (Liu et al., 2017). The synthesized cDNA was used as a new sample template for qRT-PCR and then stored at –20 °C until use.

    • qRT-PCR was performed using a QuantStudio 6 Flex Real-Time PCR System (Life Technologies, USA). All primer sequences were designed using Primer Premier 5.0 (Supplementary Table S2). Bacterial gene expression was normalized using 16S rDNA, L. crocea gene expression was normalized using β-actin, and the relative level of gene expression was calculated using the 2-ΔΔCt method. To assess the pathogen load of P. plecoglossicida in the infected spleens, the gyrB gene copy number was used to assess the number of P. plecoglossicida.

    • Sequencing experiments were carried out using an Illumina TruseqTM RNA Sample Prep Kit (Illumina, USA). Total RNA was extracted from tissue samples using TRIzol® reagent, with concentration and purity measured using a Nanodrop 2000, RNA integrity detected by agarose gel electrophoresis, and RIN values determined by an Agilent 2100. A single database requires a total RNA of 1 μg, concentration of ≥50 ng/μL, and OD260/280 of between 1.8 and 2.2. Eukaryotic and prokaryotic mRNAs were simultaneously obtained by removing rRNA, and fragmentation was carried out by adding a fragmentation buffer. After reverse-synthesizing the cDNA, an end repair mix was added to make it blunt-ended, and then poly(A) was added to the 3' end for ligation of the Y-shaped linker. Sequencing was performed on an Illumina HiSeq4000 sequencing platform from Majorbio Biotech Co., Ltd. (China). The RNA sequence data were deposited in the GenBank SRA database under accession Nos. PRJNA607373 and SRP176599. The basis of transcriptome data analysis is high-quality sequencing. Analysis showed that the distribution of the A/T/G/C base content was uniform (Supplementary Figure S2), the base mass distribution of the sequence data met the requirements of subsequent analysis (Supplementary Table S3), and the repeatability between the three repeated samples was good (Supplementary Figure S1).

    • Trimming and quality control of raw Illumina reads were performed using Sickle (https://github.com/najoshi/sickle) and SeqPrep (https:// github.com/jstjohn/SeqPrep) with default settings. For RNA-seq, clean data were mapped to the genome of P. plecoglossicida strain NB2011 (NCBI RefSeq accession number: NZ_ASJX00000000.1) using Bowtie2 (Langmead & Salzberg, 2012). Clean data were mapped to the genome of L. crocea (NCBI RefSeq accession number: GCF_000972845.2) using Hisat2 (Kim et al., 2015).

    • Differential expression was determined using edgeR, which performs differential expression calculations based on mRNA read count data and a negative binomial distribution model (Anders & Huber, 2010; Robinson et al., 2010). The screening criteria for significant DEMs were: FDR<0.05 and |log2FC|>=1. The DEMs were then subjected to enrichment analysis by hypergeometric distribution testing using Goatools (https://github.com/tanghaibao/goatools) and KOBAS (http://kobas.cbi.pku.edu. cn/home.do) (Xie et al., 2011).

    • All data are expressed as means±standard deviation (SD) from at least three sets of independent experiments. Data analysis was performed using SPSS 18.0 (SPSS Inc., USA), and one-way analysis of variance with Dunnett’s test was used. P-values of <0.05 were considered statistically significant.

    • Based on qRT-PCR, the expression level of dksA in the P. plecoglossicida-infected spleens at 48 hpi was six times higher than that in the in vitro culture, consistent with the transcriptome analysis results (Figure 1A). The expression level of dksA was significantly down-regulated in the five mutant strains. The mRNA levels in the shRNA-31, shRNA-49, shRNA-81, shRNA-87, and shRNA-249 mutant strains were only 32.19%, 30.80%, 61.18%, 3.77%, and 9.80% that of the wild-type strain, respectively (Figure 1B). The shRNA-87 mutant strain (hereinafter referred to as the dksA-RNAi or RNAi strain) exhibited the lowest dksA RNA level and was thus selected for subsequent analysis. The growth rates of the dksA-RNAi strain and wild-type strain of P. plecoglossicida were determined, although no significant differences between the two strains were observed (Figure 1C).

      Figure 1.  Construction and characterization of dksA-RNAi strains of P. plecoglossicida

    • Compared with the wild-type strain, the dksA-RNAi strain of P. plecoglossicida exhibited a significant decrease in virulence, as observed by the 25% increase in the survival rate of infected L. crocea and 48 h delay in first death. No deaths were recorded in the negative control group of L. crocea injected with PBS (Figure 2A). Furthermore, there were significant differences in spleen appearance between the two groups of L. crocea infected by the dksA-RNAi or wild-type strains of P. plecoglossicida. The spleens of L. crocea infected by the wild-type strain showed a large number of typical white nodules on the surface at 60 hpi, whereas the spleens of L. crocea infected by the dksA-RNAi strain displayed only a small number of white spots on the surface (Figure 2B).

      Figure 2.  Infection of L. crocea by two strains of P. plecoglossicida

      Throughout the infection process, the pathogen load in the spleens of L. crocea infected with the dksA-RNAi strain was always lower than that in L. crocea infected with the wild-type strain of P. plecoglossicida, and showed a tendency to increase gradually with the increase in infection time (Figure 2C). In vivo expression levels of dksA in the dksA-RNAi and wild-type strains of P. plecoglossicida were high throughout the infection process, and the expression levels in the dksA-RNAi strain were always lower than that in the wild-type strain (Figure 2D).

    • We used edgeR software to calculate gene expression levels. The screening criteria for significant differentially expressed genes (DEGs) were FDR<0.05 and |log2FC|>=1. From the constructed volcano map, a total of 4 988 P. plecoglossicida mRNAs were obtained from the transcriptome of the L. crocea spleens infected with the dksA-RNAi strain. Compared with the wild-type strain, we identified 145 differentially expressed P. plecoglossicida mRNAs in the dksA-RNAi strain-infected spleens, 24 of which were up-regulated and 121 of which were down-regulated (Figure 3A). The heat map (Figure 3B) shows pathogenic genes whose differential expression in the three samples exhibited >2-fold change by transcriptome sequencing. The most variable up-regulated gene was L321_RS17380 (log2FC=5.07), and the gene with the largest down-regulated fold-change was flgC (log2FC=–15.02) (Figure 3B). Five up-regulated and down-regulated DEGs of the pathogen were randomly selected for qRT-PCR detection. The qRT-PCR results were consistent with the transcriptome sequencing results (Figure 3C).

      Figure 3.  DEG enrichment analysis of pathogen transcriptome data

    • Gene Ontology (GO) enrichment analysis of P. plecoglossicida was performed based on Goatools using Fisher’s exact test. Results showed that there were 292 enriched terms, including 64 significantly enriched terms related to cellular component (12), biological process (45), and molecular function (7). The top 10 terms were organelle part, cell part, cellular nitrogen compound biosynthetic process, bacterial-type flagellum-dependent cell motility, structural molecule activity, cilium or flagellum-dependent cell motility, archaeal or bacterial-type flagellum-dependent cell motility, cell motility, movement of cell or subcellular component, and cellular process (Figure 4A). KOBAS was applied for KEGG pathway enrichment analysis, with 69 terms found to be enriched, including three significantly enriched terms (i.e., flagellar assembly, photosynthesis, and ribosome) (Figure 4B). In total, DEGS were enriched in 30 terms (P<0.01), as displayed in detail in Figure 4C. Compared with the wild-type strain, most DEGs were down-regulated, except for L321-RS02785, L321-RS02790, L321-RS14705, L321-RS17525, L321-RS12450, and L321-RS11375, which were up-regulated and primarily enriched in single-organism cellular process (GO:0044763) and cellular process (GO:0009987) (Figure 4C).

      Figure 4.  GO and KEGG enrichment analysis of pathogen transcriptome

    • The software used for gene expression calculation and the screening criteria for significant DEGs were the same as those used for the pathogen. From the constructed volcanic map, 27 520 L. crocea mRNAs were obtained from the spleen transcriptome of L. crocea infected with the dksA-RNAi strain. Compared with the wild-type-infected L. crocea, we identified 970 differentially expressed P. plecoglossicida mRNAs in the dksA-RNAi strain-infected spleens, 366 of which were up-regulated and 604 of which were down-regulated (Figure 5A). The heat map (Figure 5B) shows the expression levels of the top 50 up-regulated and down-regulated genes in the three samples of each group. The greatest fold-change in the up-regulated genes was for LOC104939330 (log2FC=8.51), and greatest fold-change in the down-regulated genes was for CPA2 (log2FC=–5.4). Five up-regulated and down-regulated DEGs of the host were randomly selected for qRT-PCR detection. The qRT-PCR results of these genes were consistent with the transcriptome sequencing results (Figure 5C).

      Figure 5.  DEG enrichment analysis of L. crocea transcriptome data

    • GO enrichment analysis of DEGs in L. crocea identified 292 enriched terms, including 21 significantly enriched terms related to cellular component (1), biological process (9), and molecular function (11) (Figure 6A). The top 10 significantly enriched terms for the DEGs are displayed in detail in Figure 6B. Based on z-score, more down-regulated genes were enriched in endopeptidase activity (GO: 0004175), with peptidase activity acting on L-amino acid peptides (GO: 0070011) showing the highest significance (Figure 6B). KEGG analysis identified 265 enriched KEGG pathways, including four significantly enriched immune-related pathways (i.e., Toll-like receptor signaling pathway, tumor necrosis factor (TNF) signaling pathway, hematopoietic cell lineage, and cytokine-cytokine receptor interaction) (Figure 6C). According to the heat map, compared with the wild-type strain-infected spleens, 48 DEGs (39 down-regulated, nine up-regulated) were enriched in the cytokine-cytokine receptor interaction (ko04060) pathway. We identified 14 DEGs (11 down-regulated, three up-regulated) enriched in the hematopoietic cell lineage (ko04640) pathway. In addition, 17 DEGs (16 down-regulated genes, one up-regulated) were enriched in the Toll-like receptor signaling pathway (ko04620) (Figure 8D).

      Figure 6.  GO and KEGG enrichment analysis of host transcriptome

      According to the mapped TNF signaling pathway, many genes in the spleen infected with the silenced strain showed significant changes in pathways compared to the spleen infected with the wild-type strain, with 23 genes found to be significantly down-regulated. Intracellular signaling (negative), transcription factors, leukocyte recruitment, cell adhesion, and other related gene expression levels changed, with major changes in inflammatory cytokine genes IL1b (–5.22), IL6 (–4.02) and synthesis of inflammatory mediator gene ptgs2 (–4.52) (Figure 7).

      Figure 7.  Schematic overview of response of TNF signaling pathway of L. crocea to P. plecoglossicida infection

    • The ability of pathogens to infect hosts is mainly regulated by virulence genes (Crofts et al., 2018; Yao et al., 2019). In recent years, dozens of virulence genes of aquatic pathogens have been identified (Rong et al., 2017; Zhang et al., 2018b), including several virulence genes of P. plecoglossicida (Wang et al., 2019b). Virulence genes can affect host immune responses (Sun et al., 2018) and are involved in host-pathogen interactions (Sun et al., 2019a; Tang et al., 2020). To date, however, no studies have reported on the effects of the dksA gene on host immune response.

      RNA interference technology can specifically reduce the expression of genes and has been widely used to explore gene functions (Zhang et al., 2019a). RNA interference technology results in different silencing efficiency of different aquatic pathogens (Guanzon & Maningas, 2018; Saleh et al., 2016), and different shRNAs have different silencing efficiency for the same gene (Sun et al., 2019b). To achieve good silencing, it is necessary to design several different shRNAs for a gene. In the present study, five shRNAs were designed to silence the dksA gene by RNAi. Among the five RNAi strains, dksA-shRNA-87 had the highest silencing efficiency (96.23%), which is more efficient than the silencing of most genes of aquatic pathogens (Ye et al., 2018; Zuo et al., 2019). The stability of gene silencing is crucial for the study of gene function. In the present study, the expression of dksA in the host was higher than that in vitro, indicating that the gene may be related to the pathogenicity of P. plecoglossicida. Moreover, the dksA gene in the dksA-RNAi strain of P. plecoglossicida was persistently silenced during the infection process and its relative expression was always lower than that of the wild-type strain. These results indicate that the RNAi technique was reliable in this study, thereby laying the foundation for subsequent research.

      In this study, L. crocea was artificially infected with P. plecoglossicida. After infection with the dksA-silenced strain, death time was delayed by 48 h and the mortality rate was reduced by 25% compared with L. crocea infected with wild-type P. plecoglossicida. After 60 h of infection, the spleens of L. crocea showed typical symptoms of visceral white spot disease, with fewer symptoms observed in the RNAi strain. RNAi of the dksA gene resulted in a decrease in the pathogenicity of P. plecoglossicida to L. crocea, suggesting that dksA may contribute to the virulence of P. plecoglossicida. Several genes have been verified to be associated with P. plecoglossicida virulence, the silencing of which results in a reduction in mortality of experimental fish (Tang et al., 2019a). The silencing of some genes also results in a dramatic decline in the virulence of pathogens, and therefore shows potential in the development of attenuated vaccines (Luo et al., 2020).

      Transcriptome sequencing is a powerful tool for studying gene function (Yang et al., 2018), and has been applied in host-pathogen interaction studies (Zhang et al., 2019b). In this study, the transcriptome of the silenced strain changed significantly during the infection process, and a total of 4 988 P. plecoglossicida mRNAs were detected in the samples. Compared with the wild-type strain, we identified 145 differentially expressed mRNAs in the dksA-RNAi strain-infected spleens, 24 of which were up-regulated and 121 of which were down-regulated. Significantly altered genes were analyzed for GO and KEGG enrichment. As a result, cellular component-, flagellum-, and ribosome assembly-related genes of P. plecoglossicida were down-regulated after dksA gene silencing, and the up-regulated gene L321-RS14705 (fliC) (log2FC=3.24) was enriched in seven functional GO terms. Research has shown that DksA and ppGpp in Escherichia coli inhibit the expression of the flagellar cascade during the stationary phase and following starvation, thus affecting flagella and ribosome assembly (Lemke et al., 2009). The loss of the fliC gene in Edwardsiella tarda can damage bacterial growth, reduce motility, decrease biofilm formation, and decrease secretion of virulence-related proteins involved in the type III secretion system (TTSS) (He et al., 2012). From this perspective, the silencing of the dksA gene in P. plecoglossicida inhibited the expression of a large number of genes related to flagella and ribosome assembly, and the significant up-regulation of the fliC gene may be a way to maintain virulence in P. plecoglossicida.

      Infected spleens were chosen for dual RNA-seq because they are an important immune organ (Chen et al., 2019). Dual RNA-seq can synchronously detect transcriptome changes in both host and pathogen (Luo et al., 2019a). During infection of L. crocea, the transcriptome changed significantly. Compared with the wild-type infected group, the most significantly down-regulated gene in the dksA-RNAi strain-infected group was CPA2. CPA2 is a member of the carboxypeptidase gene family. In bacteria, carboxypeptidases play a key role in the immune response to viral infections (Gardell et al., 1988; Godahewa et al., 2014). Based on enrichment analysis, most significantly enriched GO terms were related to peptidase activity, with the greatest impact on endopeptidase activity (GO: 0004175) and peptidase activity acting on L-amino acid peptides (GO: 0070011). Many peptidases are related to the immune response. Endopeptidases and L-amino acid peptides have immune-related functions in organisms, e.g., the role of endopeptidases in the immune response against influenza in mice (Tan et al., 2017) and L-amino acid antibacterial activity in the mucus layer of flounder Platichthys stellatus (Kasai et al., 2010). Some up-regulated genes in the pathogenic bacteria found here have not been reported previously in relation to pathogenicity, and further studies are needed to determine whether they are new virulence genes.

      KEGG enrichment analysis of the transcriptome data of L. crocea identified four significantly changed immune-related pathways, i.e., Toll-like receptor signaling pathway, TNF signaling pathway, hematopoietic cell lineage, and cytokine-cytokine receptor interaction. Hematopoietic cell lineage plays an important role in the immune response (Delves, 2020). Through hematopoietic cell lineage, hematopoietic stem cells differentiate into different blood cells, including T cells, natural killer (NK) cells, basophils, macrophages, and B cells, in response to various stimuli (Lu & Chen, 2019; Marshall et al., 2018). The cytokine-cytokine receptor interaction pathway is mainly involved in neutrophil infiltration during host immune response (Mantovani et al., 2019). Cytokines are soluble proteins secreted by donor cells in response to stimuli and transported to target cells through the circulatory system (Sharma et al., 2014). Studies have shown that Epinephelus coioides infected with L321_RS19110 gene-silenced P. plecoglossicida strains can significantly affect cytokine-cytokine receptor interactions (Zhang et al., 2018a). Toll-like receptors play an important role in a host’s ability to recognize pathogens and generate an immune response, which they regulate by promoting inflammatory cytokines (Palti, 2011; Takeda & Akira, 2015). The Toll-like receptor signaling pathway is also involved in the immune response of L. crocea to the secY gene of P. plecoglossicida (Wang et al., 2019a). TNF is a proinflammatory cytokine, which mediates inflammatory responses and regulates immune functions, with abnormal TNF signal transduction also related to inflammatory diseases (Joosten et al., 2016). To ensure body health, an organism will clear senescent and diseased cells by inducing apoptosis through the tumor necrosis factor superfamily (TNFSF) ligand (Collette et al., 2003). Studies have shown that interactions between tumor necrosis factor receptor type 1 (TNFR1) and nuclear factor kappa-B (NF-κB) are essential for maintenance of the TNFR1 pathway activity and activation of inflammatory cytokines that induce leukocyte recruitment (Alcamo et al., 2001). In addition, the TNFR1 pathway can also prepare antibodies for further pathogen clearance (Stokes et al., 2015). In the current study, significantly enriched inflammatory cytokines and cellular regulatory genes were down-regulated in the TNF pathway, similar to results reported in our previous study (Tang et al., 2019c). To summarize, the down-regulation of pro-inflammatory genes involved in the four significantly changed immune-related pathways may be due to the virulence of the RNAi strain being weaker than that of the wild-type strain. Thus, the immune response of the large yellow croakers infected with the RNAi strain was weaker than those infected with the wild-type strain.

    • In conclusion, dksA is a virulence gene of P. plecoglossicida. The silencing of dksA resulted in the down-regulation of cell component-, flagellum-, and ribosome assembly-related genes of P. plecoglossicida, thereby reducing the virulence of P. plecoglossicida. Through analysis of transcriptome data, we found that the fliC gene of the RNAi strain was significantly up-regulated in the course of infection, which may be a way in which to maintain the virulence of P. plecoglossicida. In addition, compared with those infected with the wild-type strain, the immune response of L. crocea infected by the RNAi strain was weaker. Most down-regulated GO terms in L. crocea infected with the RNAi strain were related to peptidase activity. KEGG enrichment analysis showed that genes related to inflammatory factors in four immune-related pathways were down-regulated in L. crocea infected with the RNAi strain. Therefore, L. crocea appears to be more resistant to infection by RNAi strains.

    • Supplementary data to this article can be found online.

    • The authors declare that they have no competing interests.

    • Q.P.Y., L.X.H., and Y.Q.S. designed the study. Q.P.Y. supervised the analyses. L.Y.W, L.M.Z., and W.Q.Z. performed the fish experiments. Z.X.L. cultivated the bacteria. L.Y.W. extracted DNA and RNA. L.Y.W., Z.X.L., L.X.H., Y.X.Q., and F.W. performed bioinformatics analysis. L.Y.W. wrote the manuscript with input from other authors. Q.P.Y., L.X.H., and Y.Q.S. revised the manuscript. All authors read and approved the final version of the manuscript.

参考文献 (70)
补充材料:
ZR-2020-034 Supporting Materials.doc

目录

    /

    返回文章
    返回