Dual RNA-seq reveals the effect of the flgM gene of Pseudomonas plecoglossicida on the immune response of Epinephelus coioides

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.

Prognostic prediction of lung adenocarcinoma by integrative analysis of RHOH expression and methylation

BACKGROUND AND OBJECTIVE

The development of epigenetics holds great promise for diagnosis and treatment of lung adenocarcinoma (LUAD). The purpose of this work was to analyze the correlation between Ras Homolog Gene Family Member H (RHOH) expression and methylation in patients with LUAD and its association with survival.

METHODS

Data related to gene expression, DNA methylation, and clinical features of LUAD from The Cancer Genome Atlas (TCGA) database were analyzed. A total of 50 patients were included in verification group. The methylation level of RHOH in verification group was detected by bisulfite amplicon sequencing.

RESULTS

The RHOH methylation level in TCGA cohort was significantly and negatively correlated with its expression level (Cor = -0.5, p = 2.687e-33). Patients with hypermethylation and low expression of RHOH had significantly worse prognosis than patients with hypomethylation and low expression of RHOH (TCGA: p = 0.004; validation cohort: p = 0.006, HR: 4.740, 95% CI: 1.567-14.340).

CONCLUSION

Our research revealed that RHOH may prove to be a new potential prognostic predictor for LUAD patients.

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.

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.

The contributions of fliG gene to the pathogenicity of Pseudomonas plecoglossicida and pathogen-host interactions with Epinephelus coioides

Pseudomonas plecoglossicida is a Gram-negative aerobic rod-shaped bacterium with polar flagella. It is the causative agent of visceral white spot disease in cultured fish, resulting in serious economic losses. In our previous study, RNA sequencing showed that the expression of the fliG gene in P. plecoglossicida is significantly up-regulated during infection of orange-spotted grouper (Epinephelus coioides). In this study, four P. plecoglossicida RNA interference (RNAi) mutants were successfully constructed by linking four short hairpin RNAs (shRNAs), which target different sites of the fliG gene, to pCM130/tac, respectively. The mRNA expression levels of the fliG gene in P. plecoglossicida were significantly decreased in four mutants. The shRNA-335 mutant (fliG-RNAi strain) showed the best silencing efficiency (88.2%) and was thus chosen for further analysis. Electron microscopy indicated that the flagella of the fliG-RNAi strain of P. plecoglossicida were shorter and finer than those of the wild type strain. The fliG-RNAi strain also showed significantly decreased mobility, chemotaxis, adhesion, and biofilm formation. Furthermore, compared with wild type strain infection, E. coioides infected with the fliG-RNAi strain exhibited a 0.5-d delay in the time of first death and 55% reduction in accumulated mortality, as well as milder splenic symptoms. RNAi of the fliG gene significantly affected the transcriptomes of both pathogen and host in the infected spleens of E. coioides. KEGG analysis revealed that the flagellar assembly pathway, bacterial chemotaxis pathway, and starch and sucrose metabolism pathway were significantly enriched in the pathogen at 3 days post infection (dpi). In contrast, the complement and coagulation cascade pathway and antigen processing and presentation pathway were significantly enriched in the host at 3 dpi. More immune-related pathways were enriched at 5 dpi and more differentially expressed genes were found in the complement and coagulation cascade and antigen processing and presentation pathways. Cytokine-cytokine receptor interaction, hematopoietic cell lineage, and IgA-producing intestinal immune network pathways were significantly enriched in the host at 5 dpi. These results indicate that fliG is an important virulence gene of P. plecoglossicida and contributes to the pathogenicity of P. plecoglossicida as well as pathogen-host interactions with E. coioides.

Transcriptome analysis revealed multiple immune processes and energy metabolism pathways involved in the defense response of the large yellow croaker Larimichthys crocea against Pseudomonas plecoglossicida

The large yellow croaker (Larimichthys crocea) aquaculture industry is suffering substantial financial losses caused by visceral white nodules disease resulting from Pseudomonas plecoglossicida infection. However, how L. crocea responds to P. plecoglossicida infection remains largely unknown. Here, we characterized the changes in the mRNA profile in the spleen of L. crocea upon P. plecoglossicida infection and explored the related defensive strategies. Sample clustering analysis and qRT-PCR indicated that P. plecoglossicida induced profound and reproducible transcriptome remodeling in the L. crocea spleen. Many innate immune-related genes, such as IL-17 signaling molecules, chemokines and chemokine receptors, complement components, TLR5 signaling molecules, and antimicrobial peptide hepcidins (Hamps), were upregulated by P. plecoglossicida and may play important roles in the L. crocea defense against P. plecoglossicida. The antibacterial activity of Hamp2-5 against P. plecoglossicida was further confirmed by using synthetic mature peptide of Hamp2-5. Additionally, significant enrichment of "Glycolysis/Gluconeogenesis", "Citrate cycle" and "Oxidative phosphorylation" pathways and a significant upregulation of all 6 rate-limiting enzyme genes (HK1, PFK, PKM, CS, IDH2, DLST) in the Glycolysis and Citrate cycle pathways in P. plecoglossicida-infected fish suggested that ATP synthesis may be accelerated to ensure energy supply in response to pathogenic infection. Altogether, our results not only identified the key immune-related genes and immune pathways that participated in the defense response of L. crocea against P. plecoglossicida, but also revealed a novel defensive strategy involving ATP synthesis in this species.

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.

The Effect of tonB Gene on the Virulence of Pseudomonas plecoglossicida and the Immune Response of Epinephelus coioides

Pseudomonas plecoglossicida is the causative agent of "visceral white spot disease" in cultured fish and has resulted in serious economic losses. tonB gene plays a crucial role in the uptake of nutrients from the outer membranes in Gram-negative bacteria. The previous results of our lab showed that the expression of tonB gene of P. plecoglossicida was significantly upregulated in the spleens of infected Epinephelus coioides. To explore the effect of tonB gene on the virulence of P. plecoglossicida and the immune response of E. coioides, tonB gene of P. plecoglossicida was knocked down by RNAi; and the differences between the wild-type strain and the tonB-RNAi strain of P. plecoglossicida were investigated. The results showed that all of the four mutants of P. plecoglossicida exhibited significant decreases in mRNA of tonB gene, and the best knockdown efficiency was 94.0%; the survival rate of E. coioides infected with the tonB-RNAi strain was 20% higher than of the counterpart infected with the wild strain of P. plecoglossicida. Meanwhile, the E. coioides infected with the tonB-RNAi strain of P. plecoglossicida carried less pathogens in the spleen and less white spots on the surface of the spleen; compared with the wild-type strain, the motility, chemotaxis, adhesion, and biofilm formation of the tonB-RNAi strain were significantly attenuated; the transcriptome data of E. coioides infected with the tonB-RNAi strain were different from the counterpart infected with the wild strain of P. plecoglossicida; the antigen processing and presentation pathway and the complement and coagulation cascade pathway were the most enriched immune pathways. The results indicated that tonB was a virulence gene of P. plecoglossicida; tonB gene was involved in the regulation of motility, chemotaxis, adhesion, and biofilm formation; tonB gene affected the immune response of E. coioides to P. plecoglossicida infection.

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.

An outbreak of visceral white nodules disease caused by Pseudomonas plecoglossicida at a water temperature of 12°C in cultured large yellow croaker (Larimichthys crocea) in China

Visceral white nodules disease (VWND) caused by Pseudomonas plecoglossicida is a common disease in cage-farmed large yellow croaker (Larimichthys crocea) in China. VWND usually occurred at water temperature of 16-19℃, resulting in high mortality in farmed large yellow croaker. Now, P. plecoglossicida as its pathogen has been considered nonpathogenic at 7-12℃. During February 2019, an infectious disease outbreak was observed in cage-farmed large yellow croaker at a water temperature of 12℃ in Ningde, China. This disease is characterized by white granulomatous lesions in internal organs of the diseased fish, which was similar with the symptoms of the VWND in large yellow croaker. Then, we isolated a bacterial strain named PQLYC4 from visceral lesions of the diseased fish. The experimental infection studies demonstrated that the strain PQLYC4 was the pathogen of the disease, which was further identified as P. plecoglossicida by the analysis of morphology, 16s rRNA gene homology and average nucleotide identity based on the whole genome sequence. Our results revealed that P. plecoglossicida strain PQLYC4 could cause the outbreak of the VWND at 12℃, a water temperature lower than that reported previously, thus providing new knowledges of prevalence and prevention of the VWND in large yellow croaker.

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.

Exposure to microplastics impairs digestive performance, stimulates immune response and induces microbiota dysbiosis in the gut of juvenile guppy (Poecilia reticulata)

Microplastics (MPs) are widely distributing in aquatic environment. They are easily ingested by aquatic organisms and accumulate in digestive tract especially of intestine. To explore the potential effects of MPs on intestine, here we, using juvenile guppy (Poecilia reticulata) as experimental animal, investigated the response characteristics of digestion, immunity and gut microbiota. After exposure to 100 and 1000 μg/L concentrations of MPs (polystyrene; 32-40 μm diameters) for 28 days, we observed that MPs could exist in guppy gut and induce enlargement of goblet cells. Activities of digestive enzymes (trypsin, chymotrypsin, amylase and lipase) in guppy gut generally reduced. MPs stimulated the expression of immune cytokines (TNF-α, IFN-γ, TLR4 and IL-6). Through high throughput sequencing of 16S rRNA gene, decreases in diversity and evenness and changed composition of microbiota were found in guppy gut. PICRUSt analysis revealed that MPs might have effects on intestinal microbiota functions, such as inhibition of metabolism and repair pathway. Our findings suggested that MPs could retain in the gut of juvenile guppy, impair digestive performance, stimulate immune response and induce microbiota dysbiosis in guppy gut. The results obtained here provide new insights into the potential risks of MPs to aquatic animals.

Molecular characterization of complement 9 in Epinephelus coioides and differential expression analysis of classical complement genes following Vibrio alginolyticus challenge

Vibrio alginolyticus is posting an increasing threat to survival of grouper. Classical complement cascade can trigger initiation of immunity, while complement 9 (C9) is a major complement molecule involved in final step of membrane attack complex (MAC) formation. In this study, full-length EcC9 contained an ORF sequence of 1779 bp, encoding a polypeptide of 592 amino acids. A high-level expression of EcC9 mRNA was observed in liver. Following vibrio challenge, increased expression levels of EcC1q, EcBf/C2, EcC4, EcC6, EcC7 and EcC9 mRNA were detected in liver and kidney. These results implied that elevated expression level of classical complement pathway (CCP) and terminal complement components (TCCs) may assess toxicological effect of V. alginolyticus.

Diversity of virulence-associated genes in pathogenic Aeromonas hydrophila isolates and their in vivo modulation at varied water temperatures

Most environmental parameters have no consistent effect on the expression of bacterial genes responsible for their virulence. However, as fish are poikilothermic, the possibility of temperature variation having a pronounced effect on the expression of virulence-associated gene(s) of bacteria infecting the host needs to be investigated. In this study, the diversity of virulence genes in seven Aeromonas hydrophila isolates collected from diseased fish from different parts of India was characterized, and the effect of temperature variation on the extent of expression of their virulence was investigated. All bacterial isolates were screened for a total of nine bacterial virulent genes {aerolysin, hemolysin, cytoen, outer membrane protein TS (Omp TS), elastase, flagellin, lipase, β hemolysin and type 3 secretion system}, and the diversity in their presence or absence were marked at a particular in vitro condition. Three bacterial isolates (nos. 1, 7 and 2) were selected for further study, based on their ability to cause varied mortalities (20-100%) in Labeo rohita juveniles in intraperitoneal challenge study. Further, three isolates were injected intraperitoneally into L. rohita fingerlings at three different temperatures (i.e., 20, 28 and 37 °C) and at 6 h post-challenge, the kidney samples were collected to measure the levels of all nine bacterial virulence genes using semi-quantitative PCR. The maximum level of amplicons of virulence genes in all three A. hydrophila isolates was noticed at 28 °C as compared to 37 °C and 20 °C. It was also observed that haemolysin played a more prominent role in the expression of virulence, when compared to cytoen gene. Hence, it was concluded that water temperature does play a crucial role in governing virulence gene expression, and a temperature of 28 °C would be considered as suitable for looking into the pathogenicity of A. hydrophila for conducting any challenge study with this organism in tropical environment.

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.

Construction expression and immunogenicity of a novel trivalent outer membrane protein (OmpU-A-II) from three bacterial pathogens in Japanese eels (Anguilla japonica)

Vibrio vulnificus, Edwardsiella anguillarum and Aeromonas hydrophila are three common bacterial pathogens in cultivated eels. To protect farming eels from infection by these pathogens, a trivalent outer membrane protein (OMP) containing partial sequences of OmpU from V. vulnificus, OmpA from E. anguillarum and OmpII from A. hydrophila was expressed and purified; then, the OMP was used as a vaccine to immunize Japanese eels (Anguilla japonica). Whole-blood cell proliferation, antibody titres and complement and lysozyme activities were detected at different days post-immunization (dpi), and the relative per cent survival (RPS) was determined after eels were infected with V. vulnificus, E. anguillarum or A. hydrophila at 28 dpi. The results showed that the OMP significantly stimulates the antibody titres. At 14 days after the challenge (i.e. at 28 dpi), the RPS of OMP against V. vulnificus, E. anguillarum and A. hydrophila was 20%, 70% and 11.1%, respectively. The construction, expression and immunogenicity of a trivalent Omp were reported for the first time, and this study will provide a valuable reference for the development of fish multiplex vaccines.

Large-scale profiling of the proteome and dual transcriptome in Nile tilapia (Oreochromis niloticus) challenged with low- and high-virulence strains of Streptococcus agalactiae

Streptococcus agalactiae is a common pathogen in aquatic animals, especially tilapia, that hinders aquaculture development and leads to serious economic losses. Previously, a S. agalactiae strain named HN016 was identified from infected tilapia, and the attenuated strain YM001 was subsequently obtained by continuous passaging in Tryptic Soy Broth (TSB) medium. YM001 has been demonstrated as a safe vaccine for S. agalactiae infection in tilapia. To understand the molecular bases of the virulence of these two strains, we performed proteomic and transcriptomic analysis to reveal the protein and gene expression changes in the liver and intestine during the infection process. HN016 significantly decreased the contents of white blood cells (WBCs), neutrophils (NEUs), red blood cells (RBCs) and hematocrit (HCT) and increased the levels of total protein (TP), albumin (ALB) and globulin (GLO), while no such significant differences were observed when comparing the control with YM001. During the infection process, pathogenic peptidoglycan hydrolase, CSPA and membrane proteins were significantly differentially expressed between YM001 and HN016. Furthermore, both proteome and transcriptome data showed that the complement and coagulation cascades pathway and the antigen processing and presentation pathway were stimulated in the liver and intestine, respectively, by YM001 infection compared to HN016 infection. The interaction network analysis of key virulence genes from pathogens suggested that CSPA, as a key node, affects the expression of DOLPP1, MIPEP, PA2G4, OCIAD1, G3BP1 and CLIC5 with a positive correlation. The present evidence suggests that during the infection process, CSPA was the key genes contributing to low virulence in YM001.

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 Unveils the Role of the Pseudomonas plecoglossicida fliA Gene in Pathogen-Host Interaction with Larimichthys crocea

In the present study, Larimichthys crocea and Pseudomonas plecoglossicida were selected as a host-pathogen interaction model for teleosts and prokaryotic pathogens. Five shRNAs were designed and synthesized to silence the fliA gene, all of which resulted in pronounced reductions in fliA mRNA; the mutant strain with the best silencing efficiency of 92.16% was chosen for subsequent analysis. A significant decrease in motility, intracellular survival and escape was observed for the fliA-RNAi strain of P. plecoglossicida, whereby silencing of the fliA gene led to a 30% decrease in mortality and a four-day delay in the onset of infection in L. crocea. Moreover, silencing of P. plecoglossicida fliA resulted in a significant change in both the pathogen and host transcriptome in the spleens of infected L. crocea. Kyoto Encyclopedia of Genes and Genomes (KEGG) analysis of pathogen transcriptome data showed that silencing fliA resulted in downregulation of 18 flagellum-related genes; KEGG analysis of host transcriptome data revealed that infection with the fliA-RNAi strain caused upregulation of 47 and downregulation of 106 immune-related genes. These pathogen-host interactions might facilitate clearance of P. plecoglossicida by L. crocea, with a significant decrease in fliA-RNAi P. plecoglossicida strain virulence in L. crocea.

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.

Mechanistic insight into the roles of Pseudomonas plecoglossicida clpV gene in host-pathogen interactions with Larimichthys crocea by dual RNA-seq

Large yellow croaker (Larimichthys crocea) is an economical important farmed fish in China. "Visceral White Spot Disease" caused by Pseudomonas plecoglossicida is a disease with a high mortality rate in cage-cultured L. crocea in recent years and resulted in heavy economy lossess. The dual RNA-seq results of previous study showed that the expression of clpV gene in P. plecoglossicida was significantly up-regulated during infection. RNAi significantly reduced the expression of clpV in P. plecoglossicida with maximum silencing efficiency of 96.1%. Compared with the wild type strain, infection of clpV-RNAi strain resulted in a delayed onset time and a 25% reduction in mortality of L. crocea, as well as lessening the symptoms of the spleen. The results of dual RNA-seq of L. crocea infected by clpV-RNAi strain of P. plecoglossicida changed considerably, compared with the counterpart infected with the wild strain. The KEGG enrichment analysis showed that Cytokine-cytokine receptor interaction, Toll-like receptor signaling pathway, C-type lectin receptor signaling pathway and MAPK signaling pathway of L. crocea were most affected by the silence of clpV in P. plecoglossicida. RNAi of clpV resulted in the downregulation of genes in flagella assembly pathway and a weaker immune response of host.

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

As an important pathogen in aquaculture, Pseudomonas plecoglossicida has caused heavy losses. The expression of an ABC transporter gene-L321_23611 of P. plecoglossicida at 18 °C was found significant higher than those at 28 °C by RNA-seq and qRT-PCR. RNAi significantly reduced the content of L321_23611 mRNA in P. plecoglossicida with a maximal decrease of 89.2%. Compared with the wild type strain, the infection of L321_23611-RNAi strain resulted in the reduction in mortality and the onset time delay of a kind of marine teleosts, Epinephelus coioides. The results of dual RNA-seq showed that the RNAi of L321_23611 resulted in a significant change in both pathogen and host transcriptome in the spleens of infected E. coioides. The result of GO and KEGG analysis from dual RNA-seq data showed both host genes of chemokine signaling pathway, coagulation and complement system, hematopoietic cell lineage pathway as well as hemoglobin complex GO term and pathogenic genes of bacterial-type flagellum-dependent cell mortality GO term and flagellar assembly, biosynthesis of amino acids and lysine biosynthesis systems pathways were mainly affected by L321_23611 gene of P. plecoglossicida. The results indicated that: 1. ABC transporter gene-L321_23611 was a virulent gene of P. plecoglossicida. 2. Both the activation of the host immune pathways and depression of pathogenic virulence-related pathways facilitated E. coioides to remove L321_23611-RNAi strain than the wild type strain of P. plecoglossicida.

clpV is a key virulence gene during in vivo Pseudomonas plecoglossicida infection

Interaction between bacterial pathogen and aquatic animal host is exceedingly complex, which involves large dynamic changes in gene expression during different stages of the disease. However, research on identifying key virulence genes based on the dynamics of gene expression changes of a one-sided bacterial pathogen in tissue has not been reported so far across different stages of infectious disease. The clpV for the T6SS of Pseudomonas plecoglossicida was identified for a candidate for key virulence gene based on dynamic changes of gene expression. For the Epinephelus coioides infected using clpV-RNAi strain, no deaths were observed up to 20 dpi. The spleens, kidneys and livers of all the E. coioides that received clpV-RNAi strain failed to develop visible nodules at 5-8 dpi, with the swelling gradually disappearing. The burdens of clpV-RNAi strain in the spleen and blood were greatly reduced at most of the time points after injection, and the burdens of clpV-RNAi strain in the head kidneys and trunk kidneys also had a sharp reduction from 72 to 120 hpi. This paper provides a new insight into the discovery of key virulence genes of pathogens in infected tissue systems.