Integration of RNAi and RNA-seq Reveals the Immune Responses of Epinephelus coioides to sigX Gene of Pseudomonas plecoglossicida

Genome-wide identification and characterization of large yellow croaker (Larimichthys crocea) suppressors of cytokine signaling (SOCS) in immune response to Pseudomonas plecoglossicida infection and acute hypoxia stress

The suppressor of cytokine signaling (SOCS) gene family is a group of genes involved in the negative regulation of cytokine signal transduction. The members of this family play a crucial role in regulating immune and inflammatory processes. However, comprehensive investigations of these genes have not yet been conducted in the economically significant fish large yellow croaker (Larimichthys crocea). In this study, a total of 13 SOCS genes (LcSOCS1a, LcSOCS1b, LcSOCS2, LcSOCS3a, LcSOCS3b, LcSOCS4, LcSOCS5a, LcSOCS5b, LcSOCS6, LcSOCS7a, LcSOCS7b, LcCISHa and LcCISHb) were identified and analyzed in L. crocea. The phylogenetic tree revealed a high conservation of SOCS genes in evolution, and the gene structure and motif analysis indicated a high similarity in the structure of LcSOCSs in the same subfamily. In addition, the expression patterns of LcSOCSs showed that LcSOCS1b was significantly down-regulated in all time under acute hypoxia stress, but it was markedly up-regulated throughout the entire process after P. plecoglossicida infection, revealing its different immune effects to two stresses. Besides, LcSOCS2a, LcSOCS6 and LcSOCS7a only participated in acute hypoxic stress, while LcSOCS5a was more sensitive to P. plecoglossicida infection. In summary, these results indicated that SOCS genes were involved in stress responses to both biological and non-biological stimuli, setting the foundation for deeper study on the functions of SOCS genes.

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.

Evaluation of immunogenicity and protective efficacy of the outer membrane porin F (OprF) against Pseudomonas plecoglossicida in large yellow croaker (Larimichthys crocea)

Large yellow croaker (Larimichthys crocea) farming dominates the marine aquaculture industry in China. However, the epidemic outbreaks of visceral white nodules disease (VWND), caused by bacterial pathogen Pseudomonas plecoglossicida, have emerged as a significant concern within the large yellow croaker industry. Although vaccination is considered to be an effective method for preventing and controlling P. plecoglossicida infection, there is currently no commercially available vaccine targeting this bacterium. In the present study, the outer membrane porin F (OprF) of P. plecoglossicida was characterized and revealed a high sequence similarity with that of other Pseudomonas species. The recombinant OprF protein (rOprF) produced in Escherichia coli was then evaluated for its immunogenicity and protective role against P. plecoglossicida in large yellow croaker. The rOprF was identified to have immunogenicity by Western blot using large yellow croaker anti-P. plecoglossicida sera. Additionally, the indirect immunofluorescence assay (IIFA) provided evidence indicating the surface exposure of OprF in P. plecoglossicida. Fish vaccinated twice via intraperitoneal (IP) injection with the purified rOprF combined with commercial adjuvant ISA 763A VG exhibited a relative percent survival (RPS) of 70.60% after challenge with virulent P. plecoglossicida strain through immersion. The administration of rOprF resulted in a notable increase in specific serum antibody levels and serum lysozyme activity compared to the control groups. The immune-related genes in the spleen and head kidney of rOprF-vaccinated fish were remarkably upregulated compared with the PBS-vaccinated sham group after the P. plecoglossicida challenge. In summary, the findings of this study suggest that rOprF exhibits considerable potential in inducing a robust immune response, making it a viable candidate for vaccination against P. plecoglossicida infection in large yellow croaker.

Differential Toxicity Responses between Hepatopancreas and Gills in Litopenaeus vannamei under Chronic Ammonia-N Exposure

Ammonia nitrogen is one of the main toxic substances in aquatic cultivation environments. Chronic exposure to excessive amounts of ammonia-N creates toxic consequences, retarding the growth of aquatic organisms. This study investigated the growth performance, morphological and physiological alterations, and transcriptome changes in the hepatopancreas and gills of white shrimp Litopenaeus vannamei. The results showed that there was no significant difference in the survival rate (p > 0.05), whereas growth performance was reduced significantly in the treated groups compared to the control groups (p < 0.05). Significant structural damage and vacuolation occurred in hepatopancreas and gill tissues in the treated groups. Superoxide dismutase (SOD) activity and Na+/K+-ATPase content were significantly increased by chronic ammonia-N exposure in the two tissue groups. In addition, catalase (CAT) activity and malondialdehyde (MDA) levels were significantly altered in the hepatopancreas groups (p < 0.05), whereas no differences were observed in the gill groups (p > 0.05). There were 890 and 1572 differentially expressed genes identified in the hepatopancreas (treated versus control groups) and gills (treated versus control groups), respectively, of L. vannamei under chronic ammonia-N exposure. Functional enrichment analysis revealed associations with oxidative stress, protein synthesis, lipid metabolism, and different serine proteases. The gills maintained cellular homeostasis mainly through high expression of cytoskeleton and transcription genes, whereas the hepatopancreas down-regulated related genes in the ribosome, proteasome, and spliceosome pathways. These genes and pathways are important in the biosynthesis and transformation of living organisms. In addition, both tissues maintained organismal growth primarily through lipid metabolism, which may serve as an effective strategy for ammonia-N resistance in L. vannamei. These results provided a new perspective in understanding the mechanisms of ammonia-N resistance in crustaceans.

Nasal responses to elevated temperature and Francisella noatunensis infection in Atlantic cod (Gadus morhua)

We report the histological and transcriptomic changes in the olfactory organ of Atlantic cod exposed to Francisella noatunensis. Experimental infection was performed at either 12 °C or 17 °C. Infected fish presented the classic gross pathologies of francisellosis. Nasal morpho-phenotypic parameters were not significantly affected by elevated temperature and infection, except for the number of mucus cells in the 12 °C group seven weeks after the challenge. A higher number of genes were altered through time in the group reared at 17 °C. At termination, the nasal transcriptome of infected fish in both groups was similar to the control. When both infected groups were compared, 754 DEGs were identified, many of which were involved in signalling, defence, transmembrane and enzymatic processes. In conclusion, the study reveals that elevated temperature could trigger responses in the olfactory organ of Atlantic cod and shape the nasal response to F. noatunensis 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.

The potential regulatory role of the lncRNA-miRNA-mRNA axis in teleost fish

Research over the past two decades has confirmed that noncoding RNAs (ncRNAs), which are abundant in cells from yeast to vertebrates, are no longer "junk" transcripts but functional regulators that can mediate various cellular and physiological processes. The dysregulation of ncRNAs is closely related to the imbalance of cellular homeostasis and the occurrence and development of various diseases. In mammals, ncRNAs, such as long noncoding RNAs (lncRNAs) and microRNAs (miRNAs), have been shown to serve as biomarkers and intervention targets in growth, development, immunity, and disease progression. The regulatory functions of lncRNAs on gene expression are usually mediated by crosstalk with miRNAs. The most predominant mode of lncRNA-miRNA crosstalk is the lncRNA-miRNA-mRNA axis, in which lncRNAs act as competing endogenous RNAs (ceRNAs). Compared to mammals, little attention has been given to the role and mechanism of the lncRNA-miRNA-mRNA axis in teleost species. In this review, we provide current knowledge about the teleost lncRNA-miRNA-mRNA axis, focusing on its physiological and pathological regulation in growth and development, reproduction, skeletal muscle, immunity to bacterial and viral infections, and other stress-related immune responses. Herein, we also explored the potential application of the lncRNA-miRNA-mRNA axis in the aquaculture industry. These findings contribute to an enhanced understanding of ncRNA and ncRNA-ncRNA crosstalk in fish biology to improve aquaculture productivity, fish health and quality.

Long non‐coding RNA, a supreme post‐transcriptional immune regulator of bacterial or virus‐driven immune evolution in teleost

The global aquaculture boom, fuelled by a reduction in wild population and detection of novel viruses, has created a demanding market, hence, there is a pressing need to investigate the immune system of fish, further. As the most diverse community of vertebrates and a central contributor to the progressing global aquaculture market, teleost continues to draw vast scientific interest. Recent breakthroughs in multi‐omics technologies have provided a platform to understand the role of long non‐coding RNA (lncRNA) in the host immune system during infection. Emerging evidence shows that teleost lncRNA might have a regulatory role in immune responses, mostly through lncRNA–microRNA (miRNA) sponging. Teleost lncRNA shares a functionally active short sequence complement to target the miRNA which is conserved among the several fish species. Recent report suggests that rhabdovirus exploits a lncRNA in teleost and, to dodge the host immune mechanism and negatively regulate the immune system. This observation reveals the essentiality of lncRNA in pathogen‐driven immunity in teleost. Reports available on the function of teleost lncRNA are still in early stages and experimental verifications are a limiting factor. Unravelling the lncRNA‐mediated immune regulation in fishes could be used against the invading pathogens to strengthen the aquaculture production. This review elaborates on the experimentally identified and functionally characterized lncRNA and its regulatory role in the teleost immune response during infection and pathogen‐driven host immune evolution, which could eventually lead to achieving high standards in aquaculture productivity.

Integrative analyses of mRNA and microRNA expression profiles reveal the innate immune mechanism for the resistance to Vibrio parahaemolyticus infection in Epinephelus coioides

Vibrio parahaemolyticus, as one of the main pathogens of marine vibriosis, has brought huge losses to aquaculture. However, the interaction mechanism between V. parahaemolyticus and Epinephelus coioides remains unclear. Moreover, there is a lack of comprehensive multi-omics analysis of the immune response of grouper spleen to V. parahaemolyticus. Herein, E. coioides was artificially injected with V. parahaemolyticus, and it was found that the mortality was 16.7% in the early stage of infection, and accompanied by obvious histopathological lesions in the spleen. Furthermore, 1586 differentially expressed genes were screened by mRNA-seq. KEGG analysis showed that genes were significantly enriched in immune-related pathways, Acute-phase immune response, Apoptosis, Complement system and Cytokine-cytokine receptor interaction. As for miRNA-seq analysis, a total of 55 significantly different miRNAs were identified. Further functional annotation analysis indicated that the target genes of differentially expressed miRNAs were enriched in three important pathways (Phosphatidylinositol signaling system, Lysosome and Focal adhesions). Through mRNA-miRNA integrated analysis, 1427 significant miRNA–mRNA pairs were obtained and “p53 signaling pathway”, “Intestinal immune network for IgA production” were considered as two crucial pathways. Finally, miR-144-y, miR-497-x, novel-m0459-5p, miR-7133-y, miR-378-y, novel-m0440-5p and novel-m0084-3p may be as key miRNAs to regulate immune signaling pathways via the miRNA-mRNA interaction network. The above results suggest that the mRNA-miRNA integrated analysis not only sheds new light on the molecular mechanisms underlying the interaction between host and V. parahaemolyticus but also provides valuable and new insights into resistance to vibrio infection.

Effect of Ferredoxin Receptor FusA on the Virulence Mechanism of Pseudomonas plecoglossicida

Pseudomonas plecoglossicida is an aerobic Gram-negative bacterium, which is the pathogen of “Visceral white spot disease” in large yellow croaker. P. plecoglossicida is a temperature-dependent bacterial pathogen in fish, which not only reduces the yield of large yellow croaker but also causes continuous transmission of the disease, seriously endangering the healthy development of fisheries. In this study, a mutant strain of fusA was constructed using homologous recombination technology. The results showed that knockout of P. plecoglossicida fusA significantly affected the ability of growth, adhesion, and biofilm formation. Temperature, pH, H₂O₂, heavy metals, and the iron-chelating agent were used to treat the wild type of P. plecoglossicida; the results showed that the expression of fusA was significantly reduced at 4°C, 12°C, and 37°C. The expression of fusA was significantly increased at pH 4 and 5. Cu²⁺ has a significant inducing effect on the expression of fusA, but Pb²⁺ has no obvious effect; the expression of fusA was significantly upregulated under different concentrations of H₂O₂. The expression of the fusA gene was significantly upregulated in the 0.5~4-μmol/l iron-chelating agent. The expression level of the fusA gene was significantly upregulated after the logarithmic phase. It was suggested that fusA included in the TBDR family not only was involved in the transport of ferredoxin but also played important roles in the pathogenicity and environment adaptation of P. plecoglossicida.

Operation of the Atypical Canonical Bone Morphogenetic Protein Signaling Pathway During Early Human Odontogenesis

Bone morphogenetic protein (BMP) signaling plays essential roles in the regulation of early tooth development. It is well acknowledged that extracellular BMP ligands bind to the type I and type II transmembrane serine/threonine kinase receptor complexes to trigger the BMP signaling pathway. Then, the receptor-activated Smad1/5/8 in cytoplasm binds to Smad4, the central mediator of the canonical BMP signaling pathway, to form transfer complexes for entering the nucleus and regulating target gene expression. However, a recent study revealed the functional operation of a novel BMP-mediated signaling pathway named the atypical BMP canonical signaling pathway in mouse developing tooth, which is Smad1/5/8 dependent but Smad4 independent. In this study, we investigated whether this atypical BMP canonical signaling is conserved in human odontogenesis. We showed that pSMAD1/5/8 is required for the expression of Msh homeobox 1 (MSX1), a well-defined BMP signaling target gene, in human dental mesenchyme, but the typical BMP canonical signaling is in fact not operating in the early human developing tooth, as evidenced by the absence of pSMAD1/5/8-SMAD4 complexes in the dental mesenchyme and translocation of pSMAD1/5/8, and the expression of MSX1 induced by BMP4 is mothers against decapentaplegic homolog 4 (SMAD4)-independent in human dental mesenchymal cells. Moreover, integrative analysis of RNA-Seq data sets comparing the transcriptome profiles of human dental mesenchymal cells with and without SMAD4 knockdown by siRNA displays unchanged expression profiles of pSMAD1/5/8 downstream target genes, further affirming the functional operation of the atypical canonical BMP signaling pathway in a SMAD1/5/8-dependent but SMAD4-independent manner in the dental mesenchyme during early odontogenesis in humans.

Combining of transcriptomic and proteomic data to mine immune-related genes and proteins in the liver of Cynoglossus semilaevis challenged with Vibrio anguillarum

The liver is a multi-functional organ including metabolism, substance synthesis, detoxification, and various immune functions, and its role in immunity has attracted more and more attention. However, research on the liver immune response of fish infected by pathogenic bacteria is currently lacking. In this study, the transcriptomics and proteomics of the liver of Cynoglossus semilaevis infected with Vibrio anguillarum were analyzed. A total of 1470 genes and 497 proteins were differentially expressed in the pairwise comparison of obvious symptoms of infection (HOSG), no obvious symptoms of infection (NOSG) and PBS treatment (CG). Gene ontology and KEGG enrichment pathways analysis showed that differentially expressed genes (DEGs) and differentially expressed proteins (DEPs) were mainly enriched in toll-like receptors (TLRs), complement and coagulation cascades, nucleotide oligomerization domain (NOD)-like receptors (NLRs), mitogen-activated protein kinase (MAPK) and phagosome signaling pathways, which suggested the combined action of the five pathways were significant to enhance the liver immune defense. The combination of transcriptomic and proteomic analysis showed that ITGβ1, C3, C5 and MRC1 were significantly up-regulated, which might play an important role in the liver immune response to the recognition of V. anguillarum, inflammatory response and phagocytosis. The transcriptome and proteome data we obtained provide information on some key genes and proteins for further study of the mechanism of liver immune response.

The Zinc Nutritional Immunity of Epinephelus coioides Contributes to the Importance of znuC During Pseudomonas plecoglossicida Infection

Previously, the dual RNA-seq was carried out in a Pseudomonas plecoglossicida- Epinephelus coioides infection model to investigate the dynamics of pathogen-host interplay in vivo. ZnuC, a member of ZnuCBA Zn importer, was found transcriptionally up-regulated during infection. Thus, this study aimed to assess its role during the trade-off for Zn between host and P. plecoglossicida. ICP-MS analysis and fluorescent staining showed that Zn was withheld from serum and accumulated in the spleen, with increased Zn uptake in the Golgi apparatus of macrophages after infection. Additionally, growth assay, macrophage infection and animal infection after gene knockout / silencing revealed that znuC was necessary for growth in Zn-limiting conditions, colonization, intracellular viability, immune escape and virulence of P. plecoglossicida. Further analysis with dual RNA-seq revealed associations of host's Zn nutritional immunity genes with bacterial Zn assimilation genes. IL6 and ZIP4 played key roles in this network, and markedly affected znuB expression, intracellular viability and immune escape, as revealed by gene silencing. Moreover, EMSA and GFP reporter gene analysis showed that Fur sensed changes in Fe concentration to regulate znuCBA in P. plecoglossicida. Jointly, these findings suggest a trade-off for Zn between host and P. plecoglossicida, while ZnuC is important for P. plecoglossicida Zn acquisition.

RNA Interference in Aquaculture: A Small Tool for Big Potential

For decades, the tight regulatory functions of DNA and RNA have been the focus of extensive research with the goal of harnessing RNA molecules (e.g., microRNA and small interfering RNA) to control gene expression and to study biological functions. RNA interference (RNAi) has shown evidence of mediating gene expression, has been utilized to study functional genomics, and recently has potential in therapeutic agents. RNAi is a natural mechanism and a well-studied tool that can be used to silence specific genes. This method is also used in aquaculture as a research tool and to enhance immune responses. RNAi methods do have their limitations (e.g., immune triggering); efficient and easy-to-use RNAi methods for large-scale applications need further development. Despite these limitations, RNAi methods have been successfully used in aquaculture, in particular shrimp. This review discusses the uses of RNAi in aquaculture, such as immune- and production-related issues and the possible limitations that may hinder the application of RNAi in the aquaculture industry. Our challenge is to develop a highly potent in vivo RNAi delivery platform that could complete the desired action with minimal side effects and which can be applied on a large-scale with relatively little expense in the aquaculture industry.

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.

Integrated Analysis of circRNA-miRNA-mRNA Regulatory Networks in the Intestine of Sebastes schlegelii Following Edwardsiella tarda Challenge

Sebastes schlegelii, an important aquaculture species, has been widely cultured in East Asian countries. With the increase in the cultivation scale, various diseases have become major threats to the industry. Evidence has shown that non-coding RNAs (ncRNAs) have remarkable functions in the interactions between pathogens and their hosts. However, little is known about the mechanisms of circular RNAs (circRNAs) and coding RNAs in the process of preventing pathogen infection in the intestine in teleosts. In this study, we aimed to uncover the global landscape of mRNAs, circRNAs, and microRNAs (miRNAs) in response to Edwardsiella tarda infection at different time points (0, 2, 6, 12, and 24 h) and to construct regulatory networks for exploring the immune regulatory mechanism in the intestine of S. schlegelii. In total, 1,794 mRNAs, 87 circRNAs, and 79 miRNAs were differentially expressed. The differentially expressed RNAs were quantitatively validated using qRT-PCR. Kyoto Encyclopedia of Genes and Genomes (KEGG) enrichment analysis showed that most of the differentially expressed mRNA genes and the target genes of ncRNAs were related to immune signaling pathways, such as the NF-κB signal pathway, pathogen recognition receptors related to signaling pathways (Toll-like receptors and Nod-like receptors), and the chemokine signaling pathway. Based on these differentially expressed genes, 624 circRNA-miRNA pairs and 2,694 miRNA-mRNA pairs were predicted using the miRanda software. Integrated analyses generated 25 circRNA-miRNA-mRNA interaction networks. In a novel_circ_0004195/novel-530/IκB interaction network, novel_530 was upregulated, while its two targets, novel_circ_0004195 and IκB, were downregulated after E. tarda infection. In addition, two circRNA-miRNA-mRNA networks related to apoptosis (novel_circ_0003210/novel_152/apoptosis-stimulating of p53 protein 1) and interleukin (novel_circ_0001907/novel_127/interleukin-1 receptor type 2) were also identified in our study. We thus speculated that the downstream NF-κB signaling pathway, p53 signaling pathway, and apoptosis pathway might play vital roles in the immune response in the intestine of S. schlegelii. This study revealed a landscape of RNAs in the intestine of S. schlegelii during E. tarda infection and provided clues for further study on the immune mechanisms and signaling networks based on the circRNA-miRNA-mRNA axis in S. schlegelii.

Microglia exosomal miRNA-137 attenuates ischemic brain injury through targeting Notch1

Microglia are the resident immune cells in the central nervous system and play an essential role in brain homeostasis and neuroprotection in brain diseases. Exosomes are crucial in intercellular communication by transporting bioactive miRNAs. Thus, this study aimed to investigate the function of microglial exosome in the presence of ischemic injury and related mechanism. Oxygen-glucose deprivation (OGD)-treated neurons and transient middle cerebral artery occlusion (TMCAO)-treated mice were applied in this study. Western blotting, RT-PCR, RNA-seq, luciferase reporter assay, transmission electron microscope, nanoparticle tracking analysis, immunohistochemistry, TUNEL and LDH assays, and behavioral assay were applied in mechanistic and functional studies. The results demonstrated that exosomes derived from microglia in M2 phenotype (BV2-Exo) were internalized by neurons and attenuated neuronal apoptosis in response to ischemic injury in vitro and in vivo. BV2-Exo also decreased infarct volume and behavioral deficits in ischemic mice. Exosomal miRNA-137 was upregulated in BV2-Exo and participated in the partial neuroprotective effect of BV2-Exo. Furthermore, Notch1 was a directly targeting gene of exosomal miRNA-137. In conclusion, these results suggest that BV2-Exo alleviates ischemia-reperfusion brain injury through transporting exosomal miRNA-137. This study provides novel insight into microglial exosomes-based therapies for the treatment of ischemic brain injury.

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.

Feline calicivirus strain 2280 p30 antagonizes type I interferon-mediated antiviral innate immunity through directly degrading IFNAR1 mRNA

Feline calicivirus (FCV) belongs to the Caliciviridae, which comprises small RNA viruses of both medical and veterinary importance. Once infection has occurred, FCV can persist in the cat population, but the molecular mechanism of how it escapes the innate immune response is still unknown. In this study, we found FCV strain 2280 to be relatively resistant to treatment with IFN-β. FCV 2280 infection inhibited IFN-induced activation of the ISRE (Interferon-stimulated response element) promoter and transcription of ISGs (Interferon-stimulated genes). The mechanistic analysis showed that the expression of IFNAR1, but not IFNAR2, was markedly reduced in FCV 2280-infected cells by inducing the degradation of IFNAR1 mRNA, which inhibited the phosphorylation of downstream adaptors. Further, overexpression of the FCV 2280 nonstructural protein p30, but not p30 of the attenuated strain F9, downregulated the expression of IFNAR1 mRNA. His-p30 fusion proteins were produced in Escherichia coli and purified, and an in vitro digestion assay was performed. The results showed that 2280 His-p30 could directly degrade IFNAR1 RNA but not IFNAR2 RNA. Moreover, the 5’UTR of IFNAR1 mRNA renders it directly susceptible to cleavage by 2280 p30. Next, we constructed two chimeric viruses: rFCV 2280-F9 p30 and rFCV F9-2280 p30. Compared to infection with the parental virus, rFCV 2280-F9 p30 infection displayed attenuated activities in reducing the level of IFNAR1 and inhibiting the phosphorylation of STAT1 and STAT2, whereas rFCV F9-2280 p30 displayed enhanced activities. Animal experiments showed that the virulence of rFCV 2280-F9 p30 infection was attenuated but that the virulence of rFCV F9-2280 p30 was increased compared to that of the parental viruses. Collectively, these data show that FCV 2280 p30 could directly and selectively degrade IFNAR1 mRNA, thus blocking the type I interferon-induced activation of the JAK-STAT signalling pathway, which may contribute to the pathogenesis of FCV infection.

Vaccination against FCV has been available for many years and has effectively reduced the incidence of clinical disease. However, vaccines cannot prevent infection, and vaccinated cats can still become persistently infected by FCV, suggesting that FCV has evolved several strategies for counteracting various components of the innate and adaptive immune systems. Here, we show that FCV strain 2280 is resistant to the antiviral effect of IFN. The molecular mechanism by which this occurs is that FCV 2280 infection blocks the JAK-STAT pathway through promoting the degradation of IFNAR1 mRNA by FCV p30 protein. An in vitro degradation assay demonstrated that 2280 p30, but not p30 of the vaccine strain F9, could directly and selectively decay IFNAR1 RNA. The exchange of p30 between 2280 and F9 strains using a reverse genetic system also showed that 2280 p30 is a key factor that contributes to the resistance to IFN and enhances virulence. Our findings reveal a new mechanism evolved by FCV to circumvent the host antiviral response.

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.

Characterization of a highly lethal barramundi (Lates calcarifer) model of Pseudomonas plecoglossicida infection

Pseudomonas plecoglossicida is a highly lethal causative agent associated with severe economic losses in aquaculture industry. P. plecoglossicida has been documented as a highly alarming pathogen in a wide variety of freshwater cultured fish including ayu (Plecoglossus altivelis), rainbow trout (Oncorhynchus mykiss) and pejerrey (Odontesthes bonariensis), and marine cultured fish such as large yellow croaker (Larimichthys crocea) and orange-spotted grouper (Epinephelus coioides) etc. Fish infected with P. plecoglossicida usually exhibited various symptoms, including lethargy, inappetence, disorientation, abdominal swelling with severe ascites and numerous white spots covered on the surface of spleen tissue. In present study, barramundi, zebrafish, spotted seabass and mandarinfish were investigated as potential hosts of P. plecoglossicida. Among them, barramundi was confirmed the most sensitive host fish species for P. plecoglossicida infection. Dynamic histopathology revealed that P. plecoglossicida caused various histopathological effects to barramundi: a) spleen: granulomas appeared at 2 days post infection (dpi) and matured at 4 dpi; b) liver: steatosis at 1 dpi and fat necrosis over time, and damaged the most compared to spleens and metanephros; c) metanephros: Bowman capsule space became larger and glomerulus shrank were even collapsed at 1 dpi; d) ascites: either bacterium or melanin were wrapped in cells from ascites. All these results indicated that P. plecoglossicida could cause systemic diseases with typical clinical sighs to barramundi and would be an alarming pathogen to barramundi industry.

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.

Characterization of a novel lncRNA (SETD3-OT) in turbot (Scophthalmus maximus L.)

LncRNAs have been demonstrated to play pivotal roles in virous biological processes, especially the gene expression regulation, including transcriptional regulation, posttranscriptional control and epigenetic processes. However, most of the current studies of lncRNAs are still limited in mammalian species, the investigations of functional roles of lncRNAs in teleost species are still lacking. In current study, we identified a novel lncRNA (SETD3-OT) in turbot, with 2,504 bp full-length obtained by 5' and 3' RACE, located in turbot chromosome 17, ranged from 20,933,835 to 20,936,302 bp. In addition, 8 neighboring genes of SETD3-OT were identified within 100 kbp in genome location. From the annotation of the neighboring adjacent genes, SETD3-OT might involve in regulation of cell apoptosis and cycle, the immune cell development, and the immune response against infection, and its expression pattern is similar to majority of the neighboring genes following Aeromonas salmonicida challenge. Intriguingly, SETD3-OT showed significant high expression levels in mucosal surfaces (intestine, gill and skin), and was dramatically down-regulated in these mucosal tissues following Vibrio anguillarum challenge, especially in gill and skin. In addition, SETD3-OT was distributed in nucleus, it might regulate the neighboring genes in cis or in trans. Taken together, our results provide insights for lncRNA in fish innate immunity, further studies should be conduct to explore the detailed molecular mechanism of the gene regulation between SETD3-OT and its neighboring genes.

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.

Immunization of a novel bivalent outer membrane protein simultaneously resisting Aeromonas hydrophila, Edwardsiella anguillarum and Vibrio vulnificus infection in European eels (Angullia angullia)

In cultivated European eels, Aeromonas hydrophila, Edwardsiella anguillarum and Vibrio vulnificus are three important bacterial pathogens. In this study, European eels (Anguilla anguilla) were immunized by the bivalent expression products of the outer membrane protein (Omp) gene from A. hydrophila (OmpⅡ) and E. anguillarum (OmpA), and the effects of the bivalent protein (rOmpⅡ-A) on the immune function of the European eel were detected. Three hundred eels were divided average into three groups of PBS, adjuvant and rOmp. Eels of three goups were injected intraperitoneal with 0.2 mL of PBS (0.01 mol/L, pH7.4), PBS + F (PBS mixed equal volume of freund's uncomplete adjuvant) or rOmpⅡ-A (1 mg mL^-1 rOmpⅡ-A mixed equal volume of freund's uncomplete adjuvant). Four immune-related genes expression, proliferation of whole blood cells, serum and skin mucus antibody titer, superoxide dismutase (SOD) activity and the relative percent of survival (RPS) were studied at different days (or hours) post the immunization. The results showed that the igm, lysC, mhc2 and sod gene in the liver, spleen, kidney and intestine tract were significant increased in the Omp group; On the 28 day post the immunization (dpi), blood cell proliferation was increased in the Omp group, and on the 14, 21, 28 and 42 dpi, antibody titers in serum and mucus of the Omp group were significantly higher than that of the PBS and adjuvant group, regardless of coating with bacteria or Omp antigen. The SOD activity of Omp group increased significantly in liver, kidney, skin mucus and serum from 14 to 42 dpi, especially in serum. Eels chanllenged by A. hydrophila, E. anguillarum and V. vulnificus in the bivalent Omp group showed the RPS were 83.33%, 55.56% and 44.44%, respectively. The results of this study showed that immunization of the bivalent Omp could effectively improve the immune function of European eels, and produced effectively protection to A. hydrophila and E. anguillarum infection. Simultaneously, the bivalent Omp also produced distinct cross-protection to the eels challenged by V. vulnificus.

Aryl hydrocarbon receptor is required for immune response in Epinephelus coioides and Danio rerio infected by Pseudomonas plecoglossicida

Aryl hydrocarbon receptor (AhR), a ligand-dependent transcriptional factor that responds to environmental chemicals, has been recently found to be closely associated with immune response in mammals. Pseudomonas plecoglossicida (P. plecoglossicida) is a temperature-dependent bacterial pathogen of visceral white spot disease in fish. Using dual RNA-seq, we previously evaluated the expression levels of ahr1a, ahr1b, ahr2 and cyp1a in the spleen of Epinephelus coioides at different time points after infection with P. plecoglossicida. In the present study, the expression levels of ahr1a, ahr1b, ahr2 and cyp1a in different organs of E. coioides and Danio rerio showed similar trends after being infected by P. plecoglossicida. It also was noted that liver, intestine, spleen, and heart were the most obviously affected organs, and ahr2 particularly showed a dramatically increase in the spleen. Subsequently, macrophages of E. coioides were isolated, and then infected by P. plecoglossicida, followed by quantitative reverse transcription polymerase chain reaction (qRT-PCR) assay, which revealed that the expression level of ahr1a in macrophages was significantly down-regulated, while expression levels of ahr1b, ahr2 and cyp1a were noticeably up-regulated. Eventually, it was noted that ahr1b and ahr2 were knocked-down in macrophages, and intracellular survival rate and immune escape rate of P. plecoglossicida were markedly improved. Taken together, ahr1a, ahr1b, ahr2 and cyp1a participate in the immune response to P. plecoglossicida in different organs of fish, while ahr1b and ahr2 may play pivotal roles in the immune response of spleen and macrophages.

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.

The effect of a LysR-type transcriptional regulator gene of Pseudomonas plecoglossicida on the immune responses of Epinephelus coioides

As an important pathogen in aquaculture, Pseudomonas plecoglossicida has caused heavy losses. It was determined with RNA-seq that the expression of a LysR-type transcriptional regulator gene (L321_20267) of P. plecoglossicida at 18 °C was significantly higher than that at 28 °C, which was verified by quantitative real-time PCR (qRT-PCR). RNAi significantly reduced the content of L321_20267 mRNA in P. plecoglossicida, with a maximal decrease of 90.63%. Compared with the wild-type strain, infection with the L321_20267-RNAi strain resulted in a 50% reduction in mortality and an onset time delay of Epinephelus coioides, as well as alleviation of the symptoms in E. coioides spleens. Compared with the wild-type strain of P. plecoglossicida, the L321_20267-RNAi strain resulted in a significant change in the spleen transcriptome of infected E. coioides. The results of GO and KEGG analysis showed that genes of serine hydrolase activity, the antigen processing and presentation pathway, the B cell receptor signalling pathway and the chemokine signalling pathway were most affected by the L321_20267 gene of P. plecoglossicida. Meanwhile, the immune genes were related to different numbers of miRNAs and lncRNAs, and some miRNAs were related to more than one gene. The results indicated that 1. L321_20267 is a virulence gene of P. plecoglossicida; 2. the upregulation of the immune pathways facilitated E. coioides to remove the L321_20267-RNAi strain compared with the wild-type strain of P. plecoglossicida; and 3. the immune genes were regulated by miRNA and lncRNA in a complex manner.

Protection against Pseudomonas plecoglossicida in Epinephelus coioides immunized with a cspA1-knock-down live attenuated vaccine

Pseudomonas plecoglossicida is well-known as the cause of viscera granulomas disease in fish. In this study, a cspA1 knock-down strain was constructed and tested in Epinephelus coioides to observe the changes in virulence and evaluate its potential as an attenuated live vaccine. The results showed that the cspA1 knock-down strain caused a significant reduction in the ability of biofilm formation, motility, adhesion and virulence. E. coioides vaccinated with cspA1 knock-down strain were more tolerant of the infection by wild-type P. plecoglossicida. The relative percent survival value of E. coioides vaccinated with cspA1 knock-down strain reached 80% after challenging with wild-type P. plecoglossicida. In the meanwhile, the expression level of genes associated with immunity, including IL-1β, IgM, MHC-I and MHC-II, was up-regulated after vaccination, indicating that the cspA1 knock-down strain can induce effective and durable immune response in E. coioides and it may be an effective attenuated live vaccine candidate for the prevention of infections by P. plecoglossicida.

Dual RNA-Seq Unveils Pseudomonas plecoglossicida htpG Gene Functions During Host-Pathogen Interactions With Epinephelus coioides

Pseudomonas plecoglossicida is a temperature-dependent opportunistic pathogen which is associated with a variety of diseases in fish. During the development of "white nodules" disease, the expression of htpG in P. plecoglossicida was found to be significantly up-regulated at its virulent temperature of 18°C. The infection of htpG-RNAi strain resulted in the onset time delay, reduction in mortality and infection symptoms in spleen of Epinephelus coioides, and affected the bacterial tissue colonization. In order to reveal the effect of htpG silencing of P. plecoglossicida on the virulence regulation in P. plecoglossicida and immune response in E. coioides, dual RNA-seq was performed and a pathogen-host integration network was constructed. Our results showed that infection induced the expression of host genes related to immune response, but attenuated the expression of bacterial virulence genes. Novel integration was found between host immune genes and bacterial virulence genes, while IL6, IL1R2, IL1B, and TLR5 played key roles in the network. Further analysis with GeneMANIA indicated that flgD and rplF might play key roles during the htpG-dependent virulence regulation, which was in accordance with the reduced biofilm production, motility and virulence in htpG-RNAi strain. Meanwhile, IL6 and IL1B were found to play key roles during the defense against P. plecoglossicida, while CELA2, TRY, CPA1, CPA2, and CPB1 were important targets for P. plecoglossicida attacking to the host.

Phenotypic characterization, virulence, and immunogenicity of Pseudomonas plecoglossicida rpoE knock-down strain

Pseudomonas plecoglossicida, a temperature dependent bacterial pathogen in fish, expresses rpoE gene that is sensitive to temperature and probably critical for pathogen virulence and disease development. In this study, the rpoE silence strain rpoE-RNAi-1 was constructed by gene knock-down. The rpoE-RNAi-1 displayed significant changes in biofilm formation, swarming motility, adhesion and virulence. Meanwhile, vaccination of grouper with rpoE-RNAi-1 led to a relative percent survival (RPS) value of 85% after challenged with the wild-type P. plecoglossicida. qRT-PCR assays showed that vaccination with rpoE-RNAi-1 enhanced the expression of immune-related genes, including MHC-I, MHC-II, IgM, and IL-1β, indicating that it was able to induce humoral and cell-mediated immune response in grouper. These results validated the possibility of rpoE as a potential target for constructing P. plecoglossicida live attenuated vaccine.

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

Pseudomonas plecoglossicida is an important and highly pathogenic bacterium for aquaculture and causes serious losses. The expression level of flgM was found to be significantly upregulated post-infection compared with in vitro results, which was confirmed by quantitative real-time PCR. RNAi significantly reduced the expression level of flgM mRNA of P. plecoglossicida. Compared with infection with the wild-type strain, infection with the flgM-RNAi strain resulted in a delay in death and a 75% reduction in the mortality of Epinephelus coioides, followed by alleviation of the symptoms in E. coioides spleen. Moreover, compared with infection with the wild-type strain, infection with the flgM-RNAi strain of P. plecoglossicida resulted in a significant change in the transcriptome of the spleens of infected E. coioides and P. plecoglossicida. KEGG analysis for E. coioides showed that genes of 17 immune pathways were most affected by flgM-RNAi of P. plecoglossicida. Among them, the expression of mhc2, zap70, rhoh, tlr2, ca79a, hcst and cd32 in E. coioides spleen was predicted to be negatively related to flgM in P. plecoglossicida but positively related to genes involved in communication, metabolism and motility.

Dual RNA-Seq reveals the role of a transcriptional regulator gene in pathogen-host interactions between Pseudomonas plecoglossicida and Epinephelus coioides

Pseudomonas plecoglossicida is a highly pathogenic bacterium for maricultured fish and causes serious losses. A transcriptional regulator gene RK21_RS10315 was found up-regulated during the whole infection process, which was confirmed by qRT-PCR. Five shRNA were designed to silence RK21_RS10315 gene, and the gene expression was reduced up to 96.1%. Compared with the counterpart infected with wild type strain, the infection of RK21_RS10315-RNAi strain resulted in the death time delay, and 90% reduction in mortality of Epinephelus coioides, as well as the alleviation in the symptoms of E. coioides spleen. Moreover, compared with the fish infected with wild type strain, the infection of RK21_RS10315-RNAi strain of P. plecoglossicida resulted in a significant change both in transcriptome of spleen of infected E. coioides and P. plecoglossicida. The KEGG analysis showed that genes of 16 immune pathways in E. coioides were affected by the silence of RK21_RS10315 of P. plecoglossicida. Among them, intestinal immune network for IgA production pathway and leukocyte transendothelial migration pathway were more prominent than other pathways. 19 euk-DEMs in these immune pathways had varying degrees of correlation with 19 pro-DEMs, and the expression of ipxA, grpE, yhbJ, truD and suhB from 19 pro-DEMs were predicted more related to RK21_RS10315 in P. plecoglossicida.

A metabolomic investigation into the temperature-dependent virulence of Pseudomonas plecoglossicida from large yellow croaker (Pseudosciaena crocea)

Pseudomonas plecoglossicida is associated with multiple fish diseases, and temperature is one of the most important environmental factors related to its outbreak. To elucidate the influence of temperature variation on the pathogen, the global metabolomics of P. plecoglossicida (NZBD9) were analysed at the virulent (18°C) and avirulent (12°C and 28°C) temperatures. The result showed that the levels of Phosphoric acid, Tyrosine, Spermidine and Sucrose were significantly reduced，while Itaconic acid, Glucaric acid and Isomaltose were increased in P. plecoglossicida at 18°C. These metabolic adjustments assist P. plecoglossicida to survive in adverse environments, proliferate in the host, colonize and resist host immune clearance during the initial steps of infection. The results suggested that L321_03626 and L321_18122 genes played a key role in the regulation of these metabolic adaptions and thus regulated P. plecoglossicida virulence at virulent temperature, which was proved by further gene silencing and artificial infection. The present study, for the first time, determines the P. plecoglossicida metabolomic responses to temperature variation, which is helpful to explore its pathogenic mechanism and provides reference for disease control.

The regulation of oxidative phosphorylation pathway on Vibrio alginolyticus adhesion under adversities

Vibrio alginolyticus is one of the most important pathogens in mariculture and leading to heavy losses. After treatment with Cu²⁺, Pb²⁺, and low pH, the expression of oxidative phosphorylation pathway genes, including coxA, coxB, coxC, ccoN, ccoO, and ccoQ, was found commonly downregulated by RNA‐seq as well as quantitative real‐time PCR. RNAi significantly reduced the expression of coxA, coxB, coxC, ccoN, ccoO, and ccoQ in V. alginolyticus. Compared with the wild‐type strain, the adhesion abilities of RNAi strains of V. alginolyticus were significantly impaired, as well as their cytochrome C oxidase activity. ccoQ appeared to be more important in the regulation of bacterial adhesion in these target genes, while ccoO was relatively weak in the regulation of the adhesion. Meanwhile, the changes of temperature, salinity, pH, and starvation affected coxA, coxB, coxC, ccoN, ccoO, and ccoQ expression remarkably. These findings indicated that: the oxidative phosphorylation pathway is a critical regulator of adhesion in V. alginolyticus; coxA, coxB, coxC, ccoN, ccoO, and ccoQ regulate the bacterial adhesion in response to environmental changes such as temperature, salinity, pH, and starvation.

Mechanisms underlying the virulence regulation of new Vibrio alginolyticus ncRNA Vvrr1 with a comparative proteomic analysis

The incidence of Vibrio alginolyticus infections has increased in recent years due to the influence of climate change and rising sea temperature. Vibrio virulence regulatory RNA 1 (Vvrr1) is a newly found noncoding RNA (ncRNA) predicted to be closely related to the adhesion ability of V. alginolyticus based on the previous RNA-seq. In this study, the target genes of Vvrr1 were fully screened and verified by constructing Vvrr1-overexpressing strains and using the proteome sequencing technology. Pyruvate kinase I (pykF) gene was predicted to be a chief target gene of Vvrr1 involved in virulence regulation. The adhesion ability, biofilm formation and virulence were significantly reduced in the Vvrr1-overexpressing and the pykF-silenced strain compared with the wild strains. Similar to the overexpression of Vvrr1, the silencing of pykF also reduced the expression level of virulence genes, such as ndk, eno, sdhB, glpF, and cysH. Meanwhile, by constructing the "pykF-GFP" fusion expression plasmid and using the GFP reporter gene analysis in Escherichia coli, the fluorescence intensity of the strain containing Vvrr1 whole ncRNA sequence vector was found to be significantly weakened. These indicated that Vvrr1 participated in the virulence regulation mechanism of V. alginolyticus by interacting with the virulence gene pykF.