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

Dual RNA-Seq Unveils Candidate Key Virulence Genes of Vibrio harveyi at the Early Stage of Infection in Hybrid Grouper (♀ Epinephelus polyphekadion × ♂ E. fuscoguttatus)

Vibrio harveyi is a major bacterial pathogen that causes disease in aquaculture animals worldwide. Although V. harveyi consistently harbors a range of traditional virulence genes, it remains unclear which specific genes are crucial for virulence at different infection stages. Dual RNA-seq is a cutting-edge RNA sequencing technology that is ideal for investigating the gene expression patterns of pathogens within the host, which is highly effective in identifying key virulence genes. In previous artificial infection experiments, we have identified the liver of hybrid grouper (♀ Epinephelus polyphekadion × ♂ E. fuscoguttatus) as the main target organ for pathogenic V. harveyi GDH11385 during the initial infection phase. To further explore the key virulence factors of V. harveyi at the early stage of infection, the liver of the hybrid grouper infected with strain GDH11385 was analyzed here by dual RNA-seq. The transcriptome data were compared with that of in vitro cultured bacteria. The results showed that 326 and 1140 DEGs (differentially expressed genes) were significantly up- and down-regulated, respectively, at 4 h post-infection (hpi). Further pathway enrichment analyses revealed that these up-regulated DEGs in vivo were mainly enriched in siderophore biosynthesis and transport, type VI secretion system (T6SS), flagellar assembly, glycolysis/gluconeogenesis, and ribosome. Notably, all genes involved in the metabolism and utilization of vibrioferrin (a carboxylate class of siderophore produced by Vibrio), and most of the genes within one of three T6SSs, were significantly up-regulated in vivo. This indicates that siderophore-dependent iron competition and T6SS-mediated delivery of virulence factors are vital for the successful colonization of V. harveyi at the early stage of infection. This study provides more precise clues to reveal the virulence mechanism of V. harveyi during the initial phase of infection.

Molecular characterization, expression pattern and immunologic function of CD82a in large yellow croaker (Larimichthys crocea)

Visceral white spot disease (VWND) caused by Pseudomonas plecoglossicida poses a major threat to the sustainable development of large yellow croaker (Larimichthys crocea) aquaculture. Genome-wide association analysis (GWAS) and RNA-seq research indicated that LcCD82a play an important role in resistance to visceral white spot disease in L. crocea, but the molecular mechanism of LcCD82a response to P. plecoglossicida infection is still unclear. In this study, we cloned and validated the Open Reading Frame (ORF) sequence of LcCD82a and explored the expression profile of LcCD82a in various tissues of L.crocea. In addition, two different transcript variants (LcCD82a-L and LcCD82a-S) of LcCD82a were identified that exhibit alternative splicing patterns after P. plecoglossicida infection, which may be closely related to the immune regulation during pathogenetic process of VWND. In order to explore the function of LcCD82a, we purified the recombinant protein of LcCD82a-L and LcCD82a-S. The bacterial agglutination and apoptosis function analysis showed that LcCD82a may involve in extracellular bacterial recognition, agglutination, and at the same time participate in the process of antigen presentation and induction of cell apoptosis. Collectively, our studies demonstrate that LcCD82a plays a crucial role in regulating apoptosis and antimicrobial immunity.

Pseudomonas plecoglossicida fliP gene affects the immune response of Epinephelus fuscoguttatus ♀×Epinephelus lanceolatus ♂ to infection

Pseudomonas plecoglossicida is a pathogen that causes visceral white spot disease in a variety of teleosts. The protein encoded by fliP gene is involved in the assembly of bacterial flagella, which plays a vital role in bacterial pathogenicity. However, the roles of the fliP gene on the host immune response remain unclear. Here, we compared the pathogenicity of fliP gene-deleted (ΔfliP) strain, fliP gene-complemented (C-ΔfliP) strain and wild-type (NZBD9) strain of P. plecoglossicida to hybrid grouper (Epinephelus fuscoguttatus ♀ × E. lanceolatus ♂), and explored the impacts of fliP gene on the immune response of hybrid grouper to P. plecoglossicida infection by using RNA-seq. In this study, the grouper in the ΔfliP strain-infected group had a 30% higher survival rate than those in the NZBD9 strain-infected group. In addition, the deletion of fliP gene decreased bacterial load in the spleen, intestine, liver as well as head kidney of hybrid grouper and the tissues damage were weakened. Moreover, the infection of hybrid grouper spleen by the ΔfliP strain induced 1,189 differential expression genes compared with the counterpart infected by NZBD9 strain. KEGG enrichment analysis showed that 9 immune-related pathways, 5 signal transduction pathways, and 3 signaling molecules and interaction pathways were significantly enriched. qRT-PCR analysis revealed that the ΔfliP strain mainly up-regulated the expression of inflammation related genes (IL-6, IL-12, IL-1β, IL-10, CXCL8, CXCL10) and immune regulation related genes (TLR2, P65, MyD88, P85, AKT), but down-regulated the expression of cell death related genes (FoxO1, Bim, PLK2 and LDHA) during infection. Based on the above results, fliP gene contributed to the pathogenicity of P. plecoglossicida to hybrid grouper (E. fuscoguttatus ♀ × E. lanceolatus ♂), deletion of fliP gene promoted the inflammation and immune response of hybrid grouper to P. plecoglossicida infection, which accelerating host clearance of pathogen and reducing tissue damages.

flgC gene is involved in the virulence regulation of Pseudomonas plecoglossicida and affects the immune response of Epinephelus coioides

As a pathogen of cultured teleosts, Pseudomonas plecoglossicida has caused significant economic losses. flgC plays an important role in encoding flagellar basal-body rod proteins. Our previous studies revealed the high expression of P. plecoglossicida flgC in infected Epinephelus coioides. To explore the role of flgC in the virulence of P. plecoglossicida and the immune response of E. coioides to the infection of P. plecoglossicida, flgC gene of P. plecoglossicida was knocked down by RNA interference (RNAi). The results showed that the flgC gene in all four mutants of P. plecoglossicida was significantly knocked down, and the mutant with the best knockdown efficiency of 94.3% was selected for subsequent studies. Compared with the NZBD9 strain of P. plecoglossicida, the flgC-RNAi strain showed a significantly decrease in chemotaxis, motility, adhesion, and biofilm formation. Furthermore, compared with the E. coioides infected with the NZBD9 strain, the infection of flgC-RNAi strain resulted in the infected E. coioides a 1.5-day delay in the time of first death and an 80% increase in survival rate, far fewer white nodules upon the spleen surfaces, and lower pathogen load in the spleens. RNAi of flgC significantly influenced the metabolome and transcriptome of the spleen in infected E. coioides. KEGG enrichment analysis exhibited that the Toll-like receptor signaling pathway was the most enriched immune pathway; the most significantly enriched metabolic pathways were associated with Linoleic acid metabolism, Choline metabolism in cancer, and Glycerophospholipid metabolism. Further combined analysis of transcriptome and metabolome indicated significant correlations among pantothenate and CoA biosynthesis, beta-alanine metabolism, lysosome metabolites, and related genes. These results suggested that flgC was a pathogenic gene of P. plecoglossicida; flgC was associated with the regulation of chemotaxis, motility, biofilm formation, and adhesion; flgC influenced the immune response of E. coioides to the infection of P. plecoglossicida.

Transcriptomic and metabolomic insights into the role of the flgK gene in the pathogenicity of Pseudomonas plecoglossicida to orange-spotted grouper (Epinephelus coioides)

Pseudomonas plecoglossicida is the pathogen responsible for visceral white spot disease in large yellow croaker (Larimichthys crocea) and orange-spotted grouper (Epinephelus coioides). Previously, RNA sequencing showed that P. plecoglossicida flgK gene expression was significantly up-regulated in orange-spotted grouper spleens during infection. To explore the role of flgK in P. plecoglossicida pathogenicity, RNA interference (RNAi) was performed to silence the P. plecoglossicida flgK gene, and the mutant (flgK-RNAi strain) with the best silencing efficiency (89.40%) was chosen for further study. Results showed that flgK gene silencing significantly attenuated P. plecoglossicida motility, adhesion, and biofilm formation. Compared to those fish infected with the wild-type strain of P. plecoglossicida, orange-spotted grouper infected with the flgK-RNAi strain showed a 55% increase in the survival rate and a one-day delay in time of first death, with fewer pathogens in the spleen and fewer white spots on the spleen surface. RNAi of flgK significantly affected the transcriptome and metabolome of the spleen in infected orange-spotted grouper. Kyoto Encyclopedia of Genes and Genomes (KEGG) enrichment analysis showed that the C-type lectin receptor signaling pathway was the most significantly changed immune-related pathway and the mitogen-activated protein kinase (MAPK) signaling pathway was related to multiple immune-related pathways. Furthermore, arginine biosynthesis and glycerophospholipid metabolism were the most significantly changed metabolism-related pathways. These findings suggest that flgK is a virulence gene of P. plecoglossicida. Furthermore, flgK appears to be involved in the regulation of motility, adhesion, and biofilm formation in P. plecoglossicida, as well as in the regulation of inflammatory and immune responses of orange-spotted grouper to P. plecoglossicida infection.

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.

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.

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.

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.

Staphylococcus aureus Synergized with Candida albicans to Increase the Pathogenesis and Drug Resistance in Cutaneous Abscess and Peritonitis Murine Models

The mixed species of Staphylococcus aureus and Candida albicans can cause infections on skin, mucosa or bloodstream; however, mechanisms of their cross-kingdom interactions related to pathogenesis and drug resistance are still not clear. Here an increase of S. aureus proliferation and biofilm formation was observed in S. aureus and C. albicans dual-species culture, and the synergistic pathogenic effect was then confirmed in both local (cutaneous abscess) and systemic infection (peritonitis) murine models. According to the transcriptome analysis of the dual-species culture, virulence factors of S. aureus were significantly upregulated. Surprisingly, the beta-lactams and vancomycin-resistant genes in S. aureus as well as azole-resistant genes in C. albicans were also significantly increased. The synergistic effects on drug resistance to both antibacterial and antifungal agents were further proved both in vitro and in cutaneous abscess and peritonitis murine models treated by methicillin, vancomycin and fluconazole. The synergistic interactions between S. aureus and C. albicans on pathogenesis and drug resistance highlight the importance of targeting the microbial interactions in polyspecies-associated infections.

BoDV-1 infection induces neuroinflammation by activating the TLR4/MyD88/IRF5 signaling pathway, leading to learning and memory impairment in rats

Borna disease virus (BoDV-1) can infect the hippocampus and limbic lobes of newborn rodents, causing cognitive deficits and abnormal behavior. Studies have found that neuroinflammation caused by viral infection in early life can affect brain development and impair learning and memory function, revealing the important role of neuroinflammation in cognitive impairment caused by viral infection. However, there is no research to explore the pathogenic mechanism of BoDV-1 in cognition from the direction of neuroinflammation. We established a BoDV-1 infection model in rats, and tested the learning and memory impairment by Morris water maze (MWM) experiment. RNAseq was introduced to detect changes in the gene expression profile of BoDV-1 infection, focusing on inflammation factors and related signaling pathways. BoDV-1 infection impairs the learning and memory of Sprague-Dawley rats in the MWM test and increases the expression of inflammatory cytokines in the hippocampus. RNAseq analysis found 986 differentially expressed genes (DEGs), of which 845 genes were upregulated and 141 genes were downregulated, and 28 genes were found to be enriched in the toll-like receptor (TLR) pathway. The expression of TLR4, MyD88, and IRF5 in the hippocampus was significantly changed in the BoDV-1 group. Our results indicate that BoDV-1 infection stimulates TLR4/MyD88/IRF5 pathway activation, causing the release of downstream inflammatory factors, which leads to neuroinflammation in rats. Neuroinflammation may play a significant role in learning and memory impairment caused by BoDV-1 infection.

A novel antiviral coumarin derivative as a potential agent against WSSV infection in shrimp seedling culture

White spot syndrome virus (WSSV), a double-stranded DNA virus that infects crustaceans, is the most serious viral pathogen affecting shrimp farming worldwide. To reduce the economic losses caused by WSSV, we screened a novel coumarin derivative from a small molecule drug library, N-(4-((4-(((2-oxo-2H-chromen-7-yl)oxy)methyl)-1H-1,2,3-triazol-1-yl)sulfonyl)phenyl)acetamide (N2905), to evaluate its anti-WSSV effects in vivo. We determined that compound N2905, up to a concentration of 20 mg/L, significantly decreased the number of WSSV copies in Litopenaeus vannamei post-larvae, with a maximum inhibitory rate of > 90 %, and increased the survival rate of WSSV-infected post-larvae. Pre-treatment and post-treatment assays indicated that N2905 could treat, but not prevent, WSSV infections. When WSSV was preincubated with N2905 for 1-4 h, the incidence of viral infections was significantly reduced and survival time of post-larvae extended to 120 h. A stability study of N2905 provided a reference for its practical use. Considering the antiviral stability of N2905 in culture water within 2 d, continuous N2905 exchange was performed, showing a significant decrease in viral load at 120 h post-infection (hpi) and a 55 % increase in survival of WSSV-infected post-larvae. Overall, our study demonstrated the potential of N2905 as an antiviral agent.

Identification and expression analysis of Langerhans cells marker Langerin/CD207 in grasscarp, Ctenopharyngodon idella

Langerhans cells (LCs) play an essential role in the initiation of immune response and maintenance of immune tolerance. However, the function and the molecular markers of grass carp LCs remains unclear. The grass carp LCs were firstly identified by immunofluorescence (IF) using a commercial anti-human Langerin/CD207 polyclonal antibody (pAb) and transmissionelectronmicroscope (TEM) technology in this study. After that, a cDNA sequence that homology with human and mouse CD207 gene was obtained by the bBLASTn program in NCBI. The open reading frame (ORF) of the grass carp CD207 gene contains 903 bp encoding 300 amino acids which consisted of a transmembrane domain, a coiled-coil domain and a CLECT domain. Furthermore, the result of quantitative real-time PCR (qRT-PCR) indicated that this gene was expressed in all tested tissues, and mainly expressed in immune organs such as the gill, trunk kidney, head kidney, spleen and skin. To explore the role of CD207 gene in the immune responses induced by bacteria, an immersed infection model of grass carp with Flavobacterium columnare was constructed, and the optimal infection dose was determined to be 1.0 × 10⁸ CFU/mL. Moreover, the qRT-PCR results indicated that the expression levels of CD207 gene were significantly upregulated at 6 h, 12 h, 1 d, 3 d and 7 d in the spleen, and significantly downregulated at 5 d in the head kidney, at 12 h and 5 d in the gill, and at all time points in the skin after F. columnare infection. This result suggested that the grass carp CD207 gene may play an important role in antigen processing and presentation. Our results in this study suggested that CD207 gene is also existed in teleosts, and this study provided a molecular basis to analyzed the biological function of grass carp CD207 gene and the critical roles of LCs in the immune responses induced by bacterial infections.