Dual RNA-Seq uncovers the function of an ABC transporter gene in the host-pathogen interaction between Epinephelus coioides and Pseudomonas plecoglossicida

The bacterial pathogen Pseudomonas plecoglossicida, its epidemiology, virulence factors, vaccine development, and host-pathogen interactions

OBJECTIVE

Pseudomoans plecoglossicida has been identified as a fish pathogen since 2000 and has caused serious infections in cultured Large Yellow Croakers Larimiththys crocea in coastal eastern China during recent years.

METHODS

Published literatures of this pathogen have been reviewed.

RESULT

Several strains with high genomic similarity have been isolated and identified; the bacteria induce natural infection at lower water temperatures (12.0-25.5°C) and induce numerous granulomas and nodules in the visceral organs of croakers. Researchers have investigated the epidemiology of P. plecoglossicida infection, identified major virulence factors, searched for pathogenic genes, analyzed host-pathogen interactions, and endeavored to develop efficient vaccines.

CONCLUSION

This paper provides an overview of these research advances to elucidate the virulence mechanisms of the pathogen and to promote vaccine development against infection.

Applying genetic technologies to combat infectious diseases in aquaculture

Disease and parasitism cause major welfare, environmental and economic concerns for global aquaculture. In this review, we examine the status and potential of technologies that exploit genetic variation in host resistance to tackle this problem. We argue that there is an urgent need to improve understanding of the genetic mechanisms involved, leading to the development of tools that can be applied to boost host resistance and reduce the disease burden. We draw on two pressing global disease problems as case studies-sea lice infestations in salmonids and white spot syndrome in shrimp. We review how the latest genetic technologies can be capitalised upon to determine the mechanisms underlying inter- and intra-species variation in pathogen/parasite resistance, and how the derived knowledge could be applied to boost disease resistance using selective breeding, gene editing and/or with targeted feed treatments and vaccines. Gene editing brings novel opportunities, but also implementation and dissemination challenges, and necessitates new protocols to integrate the technology into aquaculture breeding programmes. There is also an ongoing need to minimise risks of disease agents evolving to overcome genetic improvements to host resistance, and insights from epidemiological and evolutionary models of pathogen infestation in wild and cultured host populations are explored. Ethical issues around the different approaches for achieving genetic resistance are discussed. Application of genetic technologies and approaches has potential to improve fundamental knowledge of mechanisms affecting genetic resistance and provide effective pathways for implementation that could lead to more resistant aquaculture stocks, transforming global aquaculture.

Potential dsRNAs can be delivered to aquatic for defense pathogens

The use of antibiotics to facilitate resistance to pathogens in aquatic animals is a traditional method of pathogen control that is harmful to the environment and human health. RNAi is an emerging technology in which homologous small RNA molecules target specific genes for degradation, and it has already shown success in laboratory experiments. However, further research is needed before it can be applied in aquafarms. Many laboratories inject the dsRNA into aquatic animals for RNAi, which is obviously impractical and very time consuming in aquafarms. Therefore, to enable the use of RNAi on a large scale, the methods used to prepare dsRNA need to be continuously in order to be fast and efficient. At the same time, it is necessary to consider the issue of biological safety. This review summarizes the key harmful genes associated with aquatic pathogens (viruses, bacteria, and parasites) and provides potential targets for the preparation of dsRNA; it also lists some current examples where RNAi technology is used to control aquatic species, as well as how to deliver dsRNA to the target hydrobiont.

Metabolomics insights into the interaction between Pseudomonas plecoglossicida and Epinephelus coioides

As a highly infectious epidemic in aquaculture, Pseudomonas plecoglossicida infection results in high mortality of teleosts and serious economic losses. Host–pathogen interactions shape the outcome of an infection, yet we still understand little about the molecular mechanism of these pathogen-mediated processes. Here, a P. plecoglossicida strain (NZBD9) and Epinephelus coioides were investigated as a model system to characterize pathogen-induced host metabolic remodeling over the course of infection. We present a non-targeted metabolomics profiling of E. coioides spleens from uninfected E. coioides and those infected with wild-type and clpV-RNA interference (RNAi) strains. The most significant changes of E. coioides upon infection were associated with amino acids, lysophospatidylcholines, and unsaturated fatty acids, involving disturbances in host nutritional utilization and immune responses. Dihydrosphingosine and fatty acid 16:2 were screened as potential biomarkers for assessing P. plecoglossicida infection. The silencing of the P. plecoglossicida clpV gene significantly recovered the lipid metabolism of infected E. coioides. This comprehensive metabolomics study provides novel insights into how P. plecoglossicida shape host metabolism to support their survival and replication and highlights the potential of the virulence gene clpV in the treatment of P. plecoglossicida infection in aquaculture.

Transcriptomic analysis of Pseudomonas ogarae F113 reveals the antagonistic roles of AmrZ and FleQ during rhizosphere adaption

Rhizosphere colonization by bacteria involves molecular and cellular mechanisms, such as motility and chemotaxis, biofilm formation, metabolic versatility, or biosynthesis of secondary metabolites, among others. Nonetheless, there is limited knowledge concerning the main regulatory factors that drive the rhizosphere colonization process. Here we show the importance of the AmrZ and FleQ transcription factors for adaption in the plant growth-promoting rhizobacterium (PGPR) and rhizosphere colonization model Pseudomonas ogarae F113. RNA-Seq analyses of P. ogarae F113 grown in liquid cultures either in exponential and stationary growth phase, and rhizosphere conditions, revealed that rhizosphere is a key driver of global changes in gene expression in this bacterium. Regarding the genetic background, this work has revealed that a mutation in fleQ causes considerably more alterations in the gene expression profile of this bacterium than a mutation in amrZ under rhizosphere conditions. The functional analysis has revealed that in P. ogarae F113, the transcription factors AmrZ and FleQ regulate genes involved in diverse bacterial functions. Notably, in the rhizosphere, these transcription factors antagonistically regulate genes related to motility, biofilm formation, nitrogen, sulfur, and amino acid metabolism, transport, signalling, and secretion, especially the type VI secretion systems. These results define the regulon of two important bifunctional transcriptional regulators in pseudomonads during the process of rhizosphere colonization.

Immune responses and inorganic ion transport regulations of Epinephelus coioides in response to L321_RS13075 gene of Pseudomonas plecoglossicida

Pseudomonas plecoglossicida is a well-known pathogen of viscera granulomas disease in fish, which has led to severe economic losses. In our previous study, L321_RS13075 was predicted to be a key virulence gene of P. plecoglossicida during the host-pathogen interaction with Epinephelus coioides. To investigate the role of L321_RS13075 in the regulation of virulence in P. plecoglossicida, a L321_RS13075 knock-down strain was constructed. And a significant reduction in the ability of colonization, intracellular survival, motility, biofilm formation, and adhesion was detected in the L321_RS13075 knock-down strain. Compared with the wild-type strain, the silence of L321_RS13075 in P. plecoglossicida resulted in a significant change in the transcriptome of infected Epinephelus coioides (E. coioides). Results of COG and GO analysis on E. coioides showed that genes related to immune responses and inorganic ion transport were significantly affected by L321_RS13075 of P. plecoglossicida. Meanwhile, the interactions of the genes related to immune responses and inorganic ion transport were predicted, and the important hub genes were identified. Taken together, the results indicated that L321_RS13075 was a virulent gene of P. plecoglossicida, which significantly affected the immune responses and inorganic ion transport in E. coioides.

The contribution of exbB gene to pathogenicity of Pseudomonas plecoglossicida and its interactions with Epinephelus coioides

To study the roles of the exbB gene in Pseudomonas plecoglossicida during interactions with Epinephelus coioides, five short hairpin RNAs (shRNAs) were designed and synthesized to silence the exbB gene in P. plecoglossicida which resulted in significant reductions in exbB mRNA expression. The mutant with the best silencing efficiency (89.3%) was selected for further study. Silencing exbB in the exbB-RNA interference (RNAi) strain resulted in a 70% increase in the survival rate and a 3-day delay in the onset of infection in E. coioides. Silencing of the exbB gene also resulted in a significant decrease in the number of white spots on the spleen surface and in the spleen pathogen load. The results of dual RNA-seq showed that exbB silencing in P. plecoglossicida also resulted in a significant change in both the pathogen and host transcriptomes in the spleens of infected E. coioides. Comparative transcriptome analysis showed that silencing exbB caused significant changes in multiple signaling molecules and interaction- and immune system-related genes in E. coioides. Gene silencing also resulted in the differential expression of flagellar assembly and the bacterial secretion system in P. plecoglossicida during the infection period, and most of the DEGs were down-regulation. These host-pathogen interactions may make it easier for E. coioides to eliminate the exbB-RNAi strain of P. plecoglossicida, suggesting a significant decrease in the pathogenicity of this strain. These results indicated that exbB was a virulence gene of P. plecoglossicida which contributed a lot in the pathogen-host interactions with E. coioides.

Integration of RNA-seq and RNAi reveals the contribution of znuA gene to the pathogenicity of Pseudomonas plecoglossicida and to the immune response of Epinephelus coioides

Pseudomonas plecoglossicida is an important pathogen in aquaculture and causes serious economic losses. Our previous study indicated that znuA gene might play an important role in the pathogenicity of P. plecoglossicida. Five shRNAs were designed and synthesized to silence the znuA gene of P. plecoglossicida. Two of the five mutants of P. plecoglossicida exhibited significant reduction in the expression level of znuA mRNA with different efficiencies. The mutant with the highest silencing efficiency of 89.2% was chosen for further studies. Intrapleural injection of the znuA-RNAi strain at a dose of 10⁵  cfu/fish did not cause the death of Epinephelus coioides, and no significant signs were observed at the spleen surface of infected E. coioides, while the counterpart E. coioides infected by the same dose of wild-type strain of P. plecoglossicida all died in 5 days post-infection (dpi). The expression of znuA gene of znuA-RNAi strain in E. coioides was always lower than that in wild-type strain of P. plecoglossicida. The pathogen load in the early stage of infection was higher than that in the later stage of infection. Although the infection of the znuA-RNAi strain of P. plecoglossicida could induce the production of antibodies in E. coioides, it failed to produce a good immune protection against the infection of wild-type strain of P. plecoglossicida. Compared with the transcriptome data of E. coioides infected by the wild-type strain of P. plecoglossicida, the transcriptome data of E. coioides infected by the znuA-RNAi strain of P. plecoglossicida have altered significantly. Among them, KEGG enrichment analysis showed that the focal adhesion pathway was significantly enriched and exhibited the largest number of 302 DEMs (differentially expressed mRNAs). These results showed that the immune response of E. coioides to P. plecoglossicida infection was significantly affected by the RNAi of znuA gene.

Full-Length Transcriptome: A Reliable Alternative for Single-Cell RNA-Seq Analysis in the Spleen of Teleost Without Reference Genome

Fish is considered as a supreme model for clarifying the evolution and regulatory mechanism of vertebrate immunity. However, the knowledge of distinct immune cell populations in fish is still limited, and further development of techniques advancing the identification of fish immune cell populations and their functions are required. Single cell RNA-seq (scRNA-seq) has provided a new approach for effective in-depth identification and characterization of cell subpopulations. Current approaches for scRNA-seq data analysis usually rely on comparison with a reference genome and hence are not suited for samples without any reference genome, which is currently very common in fish research. Here, we present an alternative, i.e. scRNA-seq data analysis with a full-length transcriptome as a reference, and evaluate this approach on samples from Epinephelus coioides-a teleost without any published genome. We show that it reconstructs well most of the present transcripts in the scRNA-seq data achieving a sensitivity equivalent to approaches relying on genome alignments of related species. Based on cell heterogeneity and known markers, we characterized four cell types: T cells, B cells, monocytes/macrophages (Mo/MΦ) and NCC (non-specific cytotoxic cells). Further analysis indicated the presence of two subsets of Mo/MΦ including M1 and M2 type, as well as four subsets in B cells, i.e. mature B cells, immature B cells, pre B cells and early-pre B cells. Our research will provide new clues for understanding biological characteristics, development and function of immune cell populations of teleost. Furthermore, our approach provides a reliable alternative for scRNA-seq data analysis in teleost for which no reference genome is currently available.

Proteomics Analysis Reveals Bacterial Antibiotics Resistance Mechanism Mediated by ahslyA Against Enoxacin in Aeromonas hydrophila

Bacterial antibiotic resistance is a serious global problem; the underlying regulatory mechanisms are largely elusive. The earlier reports states that the vital role of transcriptional regulators (TRs) in bacterial antibiotic resistance. Therefore, we have investigated the role of TRs on enoxacin (ENX) resistance in Aeromonas hydrophila in this study. A label-free quantitative proteomics method was utilized to compare the protein profiles of the ahslyA knockout and wild-type A. hydrophila strains under ENX stress. Bioinformatics analysis showed that the deletion of ahslyA triggers the up-regulated expression of some vital antibiotic resistance proteins in A. hydrophila upon ENX stress and thereby reduce the pressure by preventing the activation of SOS repair system. Moreover, ahslyA directly or indirectly induced at least 11 TRs, which indicates a complicated regulatory network under ENX stress. We also deleted six selected genes in A. hydrophila that altered in proteomics data in order to evaluate their roles in ENX stress. Our results showed that genes such as AHA_0655, narQ, AHA_3721, AHA_2114, and AHA_1239 are regulated by ahslyA and may be involved in ENX resistance. Overall, our data demonstrated the important role of ahslyA in ENX resistance and provided novel insights into the effects of transcriptional regulation on antibiotic resistance in bacteria.

Transcriptome of Edwardsiella anguillarum in vivo and in vitro revealed two-component system, ABC transporter and flagellar assembly are three pathways pathogenic to European eel (Anguilla anguilla)

Edwardsiella anguillarum is one of the common bacterial pathogens for the cultivated eels in China. The aim of this study was to reveal the cause of E. anguillarum pathogenic to European eel (Anguilla anguilla) from the perspective of the transcriptome. In this study, we first prepared E. anguillarum cultured in vitro and analysed the whole transcriptome after extracting the total RNA. Then, eels were i.p injected with E. anguillarum, and total RNA were extracted from the liver of European eels 48 h after the infection. After sequencing the transcriptome, we obtained average 1.97 × 10⁸ clean reads cultured in vitro and 1.36 × 10⁵ clean reads located in vivo after annotating all reads into the genome of E. anguillarum. The whole transcriptome showed, compared to the E. anguillarum cultured in vitro, 503 significantly up and 657 significantly down-regulated different expressed genes (DEGs) were observed. KEGG analysis showed that 38 DEGs of Two-Component System, 41 DEGs of ABC transporter, and 10 DEGs flagellar assembly pathways were highly upregulated in E. anguillarum located in vivo. Then, we designed primers to analyse the up-regulated DEGs through qRT-PCR and confirmed some up-regulated DEGs. The results of this study provide important reference for the further study of pathogen-host interaction between E. anguillarum and European eel.

In vivo transcriptome analysis provides insights into host-dependent expression of virulence factors by Yersinia entomophaga MH96, during infection of Galleria mellonella

The function of microbes can be inferred from knowledge of genes specifically expressed in natural environments. Here, we report the in vivo transcriptome of the entomopathogenic bacterium Yersinia entomophaga MH96, captured during initial, septicemic, and pre-cadaveric stages of intrahemocoelic infection in Galleria mellonella. A total of 1285 genes were significantly upregulated by MH96 during infection; 829 genes responded to in vivo conditions during at least one stage of infection, 289 responded during two stages of infection, and 167 transcripts responded throughout all three stages of infection compared to in vitro conditions at equivalent cell densities. Genes upregulated during the earliest infection stage included components of the insecticidal toxin complex Yen-TC (chi1, chi2, and yenC1), genes for rearrangement hotspot element containing protein yenC3, cytolethal distending toxin cdtAB, and vegetative insecticidal toxin vip2. Genes more highly expressed throughout the infection cycle included the putative heat-stable enterotoxin yenT and three adhesins (usher-chaperone fimbria, filamentous hemagglutinin, and an AidA-like secreted adhesin). Clustering and functional enrichment of gene expression data also revealed expression of genes encoding type III and VI secretion system-associated effectors. Together these data provide insight into the pathobiology of MH96 and serve as an important resource supporting efforts to identify novel insecticidal agents.

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

Host-pathogen interactions are highly complex, involving large dynamic changes in gene expression during infection. These interactions are fundamental to understanding anti-infection immunity of hosts, as well as the pathogenesis of pathogens. For bacterial pathogens interacting with animal hosts, time-resolved dual RNA-seq of infected tissue is difficult to perform due to low pathogen load in infected tissue. In this study, an acute infection model of Larimichthys crocea infected by Pseudomonas plecoglossicida was established. The spleens of infected fish exhibited typical symptoms, with a maximum bacterial load at two days post-injection (dpi). Time-resolved dual RNA-seq of infected spleens was successfully applied to study host-pathogen interactions between L. crocea and P. plecoglossicida. The spleens of infected L. crocea were subjected to dual RNA-seq, and transcriptome data were compared with those of noninfected spleens or in vitro cultured bacteria. Results showed that pathogen-host interactions were highly dynamically regulated, with corresponding fluctuations in host and pathogen transcriptomes during infection. The expression levels of many immunogenes involved in cytokine-cytokine receptor, Toll-like receptor signaling, and other immune-related pathways were significantly up-regulated during the infection period. Furthermore, metabolic processes and the use of oxygen in L. crocea were strongly affected by P. plecoglossicida infection. The WGCNA results showed that the metabolic process was strongly related to the entire immune process. For P. plecoglossicida, the expression levels of motility-related genes and flagellum assembly-related genes were significantly up-regulated. The results of this study may help to elucidate the interactions between L. crocea and P. plecoglossicida.

Integration of RNA-seq and RNAi provides a novel insight into the immune responses of Epinephelus coioides to the impB gene of Pseudomonas plecoglossicida

Pseudomonas plecoglossicida is a Gram-negative bacterium that causes visceral white spot disease in Epinephelus coioides and leads to severe aquatic economic losses. The RNA-seq results of a previous study showed that the expression of the impB gene in P. plecoglossicida was significantly upregulated during infection. Four shRNAs were designed and synthesized to silence the impB gene in P. plecoglossicida, and the maximum silencing efficiency was 95.2%. Intraperitoneal injection of the impB-RNAi strain of P. plecoglossicida did not cause E. coioides death, and the spleens of infected fish did not show significant clinical symptoms. Although the injection of the mutant strain increased the antibody titer in E. coioides serum, it could not effectively protect E. coioides against wild strain infection. Compared with E. coioides infected with the wild type strain, the RNA-seq results for E. coioides infected with the impB-RNAi strain differed greatly. The KEGG enrichment analysis showed that key genes of the chemokine signalling pathway of E. coioides were downregulated by the silencing of impB in P. plecoglossicida. Infection with the impB-RNAi strain of P. plecoglossicida through injection did not produce good immune protection against E. coioides. The present study provides a novel insight into the immune responses of E. coioides to the impB gene of P. plecoglossicida.

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

Pseudomonas plecoglossicida is a rod-shaped, gram-negative bacterium with flagella. It causes visceral white spot disease and high mortality in Larimichthys crocea during culture, resulting in serious economic loss. Analysis of transcriptome and quantitative real-time polymerase chain reaction (PCR) data showed that dksA gene expression was significantly up-regulated after 48 h of infection with Epinephelus coioides (log 2FC=3.12, P<0.001). RNAi of five shRNAs significantly reduced the expression of dksA in P. plecoglossicida, and the optimal silencing efficiency was 96.23%. Compared with wild-type strains, the symptoms of visceral white spot disease in L. crocea infected with RNAi strains were reduced, with time of death delayed by 48 h and mortality reduced by 25%. The dksA silencing led to a substantial down-regulation in cellular component-, flagellum-, and ribosome assembly-related genes in P. plecoglossicida, and the significant up-regulation of fliC may be a way in which virulence is maintained in P. plecoglossicida. The GO and KEGG results showed that RNAi strain infection in L. crocea led to the down-regulation of inflammatory factor genes in immune-related pathways, which were associated with multiple immune response processes. Results also showed that dksA was a virulence gene in P. plecoglossicida. Compared with the wild-type strains, RNAi strain infection induced a weaker immune response in L. crocea.

Immune and gut bacterial successions of large yellow croaker (Larimichthys crocea) during Pseudomonas plecoglossicida infection

Large yellow croaker (Larimichthys crocea, LYC) aquaculture is being threatened by intensive infectious diseases. Relevant studies have focused on LYC immune responses to infection. By contrast, little is known how and to what extent the gut microbiota responds to infection. Here, we explored the interactions between LYC immune responses and gut bacterial communities during Pseudomonas plecoglossicida infection. P. plecoglossicida successfully colonized into LYC gut microbiota, resulting in an increasing mortality rate. Relative gene expressions of pro-inflammatory cytokines (TNF-α1, TNF-α2 and IL-1β) and anti-inflammatory cytokine (IL-10) were consistently and significantly induced by P. plecoglossicida infection, whereas non-specific immune enzymes activities were only enhanced at the early infection stages. P. plecoglossicida infection caused an irreversible disruption in the gut microbiota, of which infection and hours post infection constrained 16.2% and 5.6% variations, respectively. In addition, top 18 discriminatory taxa that were responsible for the difference between treatments were identified, whose abundances were significantly associated with the immune activities of LYC. Using a structural equation modeling (SEM), we found that gut bacterial communities were primarily governed by the conjointly direct (-0.33) and indirect (0) effects of infection, which subsequently affect host immune responses. Our results suggest that an irreversible dysbiosis in gut microbiota could be the causality of increasing mortality. To our knowledge, this is the first study to provide an integrated overview among pathogen infection, immune response and gut microbiota of LYC.

Effects of bacterial haemorrhagic septicemia on the immune response of Leiocassis longirostris by RNA-Seq and microRNA-Seq

Leiocassis longirostris is a common fish variety that is widely cultivated in China, during the breeding process however, it is highly susceptible to bacterial haemorrhagic septicemia, which can cause great economic loss for farmers. To understand the immune responses of L. longirostris to Aeromonas hydrophila infection, Illumina sequencing was employed to identify changes in the mRNA and miRNA in spleen tissue. In this study, a total of 92.16 and 95.61 million (M) high-quality transcriptome reads were generated from the control group (CG) and experimental group (EG) spleen samples, respectively, and 207 up-regulated and 185 down-regulated genes were identified. These genes were enriched in 29 GO terms and 30 KEGG pathways (P ≤ 0.05), including cytokine-cytokine receptor interaction and complement and coagulation cascades, with 17 up-regulated genes and 12 down-regulated genes related to immune responses in the EG relative to the CG. Based on the zebrafish genome, miRNA-seq identified a total of 343 miRNAs, of which 15 were up-regulated and 10 were down-regulated (fold-change ≥2 or ≤0.5 and P ≤ 0.05). Target gene prediction and KEGG enrichment analysis revealed that all of the target genes were concentrated in 13 pathways associated with immune response, including the mTOR signaling pathway and the TGF-beta signaling pathway. The expression patterns of 8 differentially expressed genes and 4 miRNAs involved in immune response were validated by quantitative real-time RT-PCR. These results have provided valuable insights into the molecular mechanisms underlying the immune response of L. longirostris to bacterial haemorrhagic septicemia.

Time-resolved dual RNA-seq of tissue uncovers Pseudomonas plecoglossicida key virulence genes in host-pathogen interaction with Epinephelus coioides

Bacterial pathogen-host interactions are highly dynamic, regulated processes that have been primarily investigated using in vitro assays. The dynamics of bacterial pathogen-host interplay in vivo are poorly understood. Using time-resolved dual RNA-seq in a Pseudomonas plecoglossicida-Epinephelus coioides infection model, we observed that bacterial genes encoding classical virulence factors and host genes involved in immune regulation were dynamically expressed during infection. Using network inferencing, we were able to predict interspecies regulatory networks linking bacterial virulence genes to host immune genes. Together with gene co-expression network analysis of the pathogen, secY was predicted to be a key virulence gene for P. plecoglossicida pathogenicity in the host, fliN was predicted to be a less important virulence gene. The results of bioinformatics prediction were confirmed by animal infection experiments. Our work provides the first paradigm to study dynamic alterations of bacterial pathogen and host interactions based on the elucidation of time-resolved interactive transcriptomes in vivo, and may be developed into a novel and universal method for revealing the true complexity of the bacterial infection process.

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

The interactions between host and pathogen is exceedingly complex, which involves alterations at multiple molecular layers. However, research to simultaneously monitor the alterations of transcriptome and proteome between a bacterial pathogen and aquatic animal host through integrated dual RNA-seq and dual iTRAQ of tissue during infection is currently lacking. The important role of a diguanylate cyclase gene (L321_RS15240) in pathogenicity of Pseudomonas plecoglossicida against Epinephelus coioides was suggested by previous dual RNA-seq of our lab. Then L321_RS15240-RNAi strains of P. plecoglossicida were constructed with pCM130/tac, and the mutant with the best silencing effect was selected for follow-up study. The RNAi of L321_RS15240 resulted in a significant decrease in bacterial virulence of P. plecoglossicida. The E. coioides spleens infected by wild type strain or L321_RS15240-RNAi strain of P. plecoglossicida were subjected to dual RNA-seq and dual iTRAQ, respectively. The results showed that: RNAi of L321_RS15240 led to 1)alterations of host transcriptome associated with complement and coagulation cascades, ribosome, arginine and proline metabolism, and oxidative phosphorylation; 2)high expression of host proteins which related to phagosome and metabolism responses (metabolism of glutathione, amino sugar and nucleotide sugar); 3)the highly differentially expression of host lncRNAs and miRNAs. The differentially expressed proteins and mRNAs of pathogen were different after infection, but the functions of these proteins and mRNAs were mainly related to metabolism and virulence. This study provides a new insight to comprehensively understand the gene functions of pathogens and hosts at multiple molecular layers during in vivo infection.

Dual RNA-seq uncovers the immune response of Larimichthys crocea to the secY gene of Pseudomonas plecoglossicida from the perspective of host-pathogen interactions

Pseudomonas plecoglossicida is a Gram-negative aerobic bacterium that causes high mortality and serious economic losses in some commercial marine fish. Expression of secY was found to be significantly upregulated at 18 °C compared to 28 °C by RNA-seq and qRT-PCR. All five tested recombinant vectors (pCM130/tac + shRNA) significantly reduced secY mRNA levels in P. plecoglossicida. The recombinant vector encoding shRNA-1165 exhibited the best gene-silencing efficiency, 82.4% and was used to create an RNAi strain for further studies. Compared with the wildtype strain, infections of Larimichthys crocea with the RNAi strain resulted in a 2-day delay in onset time and a 35% reduction in mortality, as well as the alleviation of spleen symptoms. The spleens of L. crocea infected by the wild type or RNAi strain of P. plecoglossicida were subjected to dual RNA-seq at 2 dpi. Compared with the wildtype strain, infection of P. plecoglossicida with the RNAi strain resulted in significant changes in the transcriptomes of both host and pathogen. KEGG analysis showed that the complement and coagulation cascade and the Toll-like receptor signalling pathway were the most enriched host pathways. In the pathogen, genes of the "Sec secretion system" were significantly downregulated. This downregulation of "Sec secretion system" genes hindered the secretion of bacterial proteins and reduced the virulence of P. plecoglossicida. Thus, it was easier for L. crocea to clear the RNAi strain of P. plecoglossicida, and the immune response was similarly reduced. The results indicated that secY was a virulence gene of P. plecoglossicida and played roles in the host-pathogen interactions of L. crocea and P. plecoglossicida.