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

Mechanism of sturgeon intestinal inflammation induced by Yersinia ruckeri and the effect of florfenicol intervention

The mechanism by which Y. ruckeri infection induces enteritis in Chinese sturgeon remains unclear, and the efficacy of drug prevention and control measures is not only poor but also plagued with numerous issues. We conducted transcriptomic and 16 S rRNA sequencing analyses to examine the differences in the intestinal tract of hybrid sturgeon before and after Y. ruckeri infection and florfenicol intervention. Our findings revealed that Y. ruckeri induced the expression of multiple inflammatory factors, including il1β, il6, and various chemokines, as well as casp3, casp8, and multiple tumor necrosis factor family members, resulting in pathological injury to the body. Additionally, at the phylum level, the relative abundance of Firmicutes and Bacteroidota increased, while the abundance of Plesiomonas and Cetobacterium decreased at the genus level, altering the composition of the intestinal flora. Following florfenicol intervention, the expression of multiple apoptosis and inflammation-related genes was down-regulated, promoting tissue repair. However, the flora became further dysregulated, increasing the risk of infection. In conclusion, our analysis of the transcriptome and intestinal microbial composition demonstrated that Y. ruckeri induces intestinal pathological damage by triggering apoptosis and altering the composition of the intestinal microbiota. Florfenicol intervention can repair pathological damage, but it also exacerbates flora imbalance, leading to a higher risk of infection. These findings help elucidate the molecular mechanism of Y. ruckeri-induced enteritis in sturgeon and evaluate the therapeutic effect of drugs on intestinal inflammation in sturgeon.

C1QBP is a critical component in the immune response of large yellow croaker (Larimichthys crocea) against visceral white spot disease caused by Pseudomonas plecoglossicida

The large yellow croaker (Larimichthys crocea) stands as a cornerstone of mariculture in China due to its significant value. However, the threat of Pseudomonas plecoglossicida infection looms large, capable of triggering "visceral white spot disease" and subsequently inflicting severe economic ramifications. Through a prior genome-wide association analysis (GWAS) aimed at understanding the resistance of the large yellow croaker to this ailment, a pivotal player emerged: the complement component 1q binding protein, aptly named LcC1qbp. This protein assumes a crucial role in the activation of the complement system. This study delves deeper into the immune response by examining the expression patterns of LcC1QBP when confronted with P. plecoglossicida. The investigation into gene expression patterns reveals LcC1qbp's widespread presence, with its highest transcriptional abundance identified in the kidney tissues. Upon infection by P. plecoglossicida, the up-regulation of LcC1qbp in major immune organs manifests at both the transcriptional and translational levels. In the context of RNA interference, transcriptome analysis of C1qbp in HEK 293T cells uncovers 1327 differentially expressed genes (DEGs), featuring 41 significant immune genes. This includes pivotal components such as C1S and C3 of the complement system, along with IL11, IL12RB2, and Myd88, among others. The outcomes of enrichment analysis spotlight the prevalence of DEGs within key pathways like immune system development, myeloid leukocyte-mediated immunity, MAPK signaling, and other immune-related routes. By unveiling the immune response mechanisms of the large yellow croaker to P. plecoglossicida infection, this study bolsters our understanding. Furthermore, it lays the groundwork for pursuing effective strategies in both preventing and treating "visceral white spot disease" in the large yellow croaker.

Potential functions of the shared bacterial taxa in the citrus leaf midribs determine the symptoms of Huanglongbing

Instruction

Citrus is a globally important fruit tree whose microbiome plays a vital role in its growth, adaptability, and resistance to stress.

Methods

With the high throughput sequencing of 16S rRNA genes, this study focused on analyzing the bacterial community, especially in the leaf midribs, of healthy and Huanglongbing (HLB)-infected plants.

Results

We firstly identified the shared bacterial taxa in the midribs of both healthy and HLB-infected plants, and then analyzed their functions. Results showed that the shared bacterial taxa in midribs belonged to 62 genera, with approximately 1/3 of which modified in the infected samples. Furthermore, 366 metabolic pathways, 5851 proteins, and 1833 enzymes in the shared taxa were predicted. Among these, three metabolic pathways and one protein showed significant importance in HLB infection. With the random forest method, six genera were identified to be significantly important for HLB infection. Notably, four of these genera were also among the significantly different shared taxa. Further functional characterization of these four genera revealed that Pseudomonas and Erwinia likely contributed to plant defense against HLB, while Streptomyces might have implications for plant defense against HLB or the pathogenicity of Candidatus Liberibacter asiaticus (CLas).

Disccusion

Overall, our study highlights that the functions of the shared taxa in leaf midribs are distinguished between healthy and HLB-infected plants, and these microbiome-based findings can contribute to the management and protection of citrus crops against CLas.

Comparative transcriptomes reveal different tolerance mechanisms to Streptococcus agalactiae in hybrid tilapia, nile tilapia, and blue tilapia

Tilapia is one of the most economically important freshwater fish farmed in China. Streptococcosis outbreaks have been extensively documented in farmed tilapia species. Hybrid tilapia (Oreochromis niloticus ♀ × O. aureus ♂) exhibit greater disease resistance than Nile tilapia (O. niloticus) and blue tilapia (O. aureus). However, the molecular mechanism underlying the enhanced tolerance of hybrid tilapia is still poorly understood. In this study, comparative transcriptome analysis was performed to reveal the different tolerance mechanisms to Streptococcus agalactiae in the three tilapia lines. In total, 1982, 2355, and 2076 differentially expressed genes were identified at 48 h post-infection in hybrid tilapia, Nile tilapia, and blue tilapia, respectively. Functional enrichment analysis indicated that numerous metabolic and immune-related pathways were activated in all three tilapia lines. The differential expression of specific genes associated with phagosome, focal adhesion, cytokine-cytokine receptor interaction, and toll-like receptor signaling pathways contributed to the resistance of hybrid tilapia. Notably, immune response genes in hybrid tilapia, such as P38, TLR5, CXCR3, CXCL12, PSTPIP1, and TFR, were generally suppressed under normal conditions but selectively induced following pathogen challenge. These results expand our knowledge of the molecular mechanisms underlying S. agalactiae tolerance in hybrid tilapia and provide valuable insights for tilapia breeding programs.

Effect of Agaricus bisporus Polysaccharides (ABPs) on anti-CCV immune response of channel catfish

Herein, the effects of Agaricus bisporus Polysaccharides (ABPs) on anti-channel catfish virus (CCV) infections to promote their application in channel catfish culture were explored. Transcriptome and metabolome analyses were conducted on the spleen of a CCV-infected channel catfish model fed with or without ABPs. CCV infections upregulated many immune and apoptosis-related genes, such as IL-6, IFN-α3, IFN-γ1, IL-26, Casp3, Casp8, and IL-10, and activated specific immunity mediated by B cells. However, after adding ABPs, the expression of inflammation-related genes decreased in CCV-infected channel catfish, and the inflammatory inhibitors NLRC3 were upregulated. Meanwhile, the expression of apoptosis-related genes was reduced, indicating that ABPs can more rapidly and strongly enhance the immunity of channel catfish to resist viral infection. Moreover, the metabonomic analysis showed that channel catfish had a high energy requirement during CCV infection, and ABPs could enhance the immune function of channel catfish. In conclusion, ABPs can enhance the antiviral ability of channel catfish by enhancing immune response and regulating inflammation. Thus, these findings provided new insights into the antiviral response effects of ABPs, which might support their application in aquaculture.

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.

miR-130a targets CiGadd45bb to modulate the inflammatory response to bacterial infection in Ctenopharyngodon idella kidney (CIK) cells

Septicemia is a systemic inflammatory response to bacterial infection that results in a hyper-inflammatory state, which could lead to septic shock and death in grass carp (Ctenopharyngodon idella). The aim of this study was to determine the underlying mechanism of microRNA (miR-130a) in bacteria-infected grass carp. Expression levels of miR-130a against Aeromonas hydrophila (A. hydrophila) infection in Ctenopharyngodon idella kidney cells (CIK) were analyzed. Luciferase reporter assay, quantitative reverse transcription-polymerase chain reaction were performed to explore whether Ctenopharyngodon idella growth arrest and DNA damage-inducible 45 (CiGadd45bb) was a target of miR-130a. MiR-130a mimic, inhibitor and miR-control were transfected to CIK respectively. After transfection, the expression levels of proinflammatory genes were determined. Here we show that CiGadd45bb as a target of miR-130a. We also confirmed that miR-130a levels were significantly higher after being stimulated for 4 h and lower after 12 h (P < 0.01). Overexpressing miR-130a strikingly inhibited p38, JNK, ERK and TNF-a genes (P < 0.01) and silencing miR-130a activated p38, JNK, ERK, TNF-a, IFN and IL-8 (P < 0.01). Our results provide a theoretical basis for studying the molecular mechanism underlying the regulation of inflammation by miR-130a in grass carp.

Dual RNA-seq of spleens extracted from channel catfish infected with Aeromonas veronii reveals novel insights into host-pathogen interactions

Interactions between host and pathogen are involving various dynamic changes in transcript expression and critical for understanding host immunity against infections and its associated pathogenesis. Herein, we established a model of channel catfish infected with Aeromonas veronii. The infected fish had prominent body surface bleeding, and the spleen showed hyperemia and swelling. Then, the spleen of channel catfish infected with A. veronii was analyzed by dual RNA sequencing (RNA-seq), and the transcriptome data were compared with uninfected channel catfish spleen or bacteria cultured in vitro. The transcript expression profile of pathogen-host interaction between A. veronii and channel catfish was successfully studied. During infection, the host was enriched for multiple immune-related signaling pathways, such as the Toll-like receptor signaling pathway, Cytokine-cytokine receptor interaction, and T cell receptor signaling pathway; and significantly upregulated for many innate immune-related genes, including IL-8. At the same time, we found that A. veronii mainly harmed the host spleen through hemolysin. Our current findings are of great significance in clarifying the pathogenesis of channel catfish induced by A. veronii and provide gene targets for developing preventive measures.

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.

Integrative time series of cellular, humoral and molecular response revealed immunotoxicity of bifenthrin to Chinese rare minnow (Gobiocypris rarus) following Pseudomonas fluorescens challenge

Bifenthrin is a common pesticide that is widespread in aquatic environments. Although it has been shown to be toxic to aquatic organisms, its immunotoxicity and mechanism are unclear. Herein, we reported the immunotoxicity of bifenthrin on adult Chinese rare minnow (Gobiocypris rarus) after 28 days of exposure to different concentrations of bifenthrin (0.1, 0.3, and 1.0 μg/L) and 36-h Pseudomonas fluorescens challenge. Bifenthrin inhibited the fish humoral immune response to bacteria by altering the lymphocyte and neutrophil ratios and decreasing the production of lysozyme, complement component 3, immunoglobulin M, and C-reactive protein, particularly were 1.0 μg/L. Bifenthrin caused intestinal damage and significantly reduced the volume of intestinal mucus at 12 and 36 hours postinjection (hpi) (p < 0.05). Moreover, 1.0 μg/L bifenthrin significantly increased the fish mortality and bacterial loads at 12 and 36 hpi (p < 0.05). RNA-seq analysis revealed several enriched genes involved in pathogen attachment and recognition, inflammatory responses, and complement system at the early-to-mid stage of infection (4-12 hpi). Overall, our results corroborated that bifenthrin induced immunotoxicity in Gobiocypris rarus, resulting in immune dysfunction of fish and increasing their sensitivity to bacterial infection and accelerating mortality. Moreover, 4-12 hpi was better than 36 hpi for analyzing immune responses against pathogen infection in fish exposed to bifenthrin.

Secreted peptidases contribute to virulence of fish pathogen Flavobacterium columnare

Flavobacterium columnare causes columnaris disease in freshwater fish in both natural and aquaculture settings. This disease is often lethal, especially when fish population density is high, and control options such as vaccines are limited. The type IX secretion system (T9SS) is required for F. columnare virulence, but secreted virulence factors have not been fully identified. Many T9SS-secreted proteins are predicted peptidases, and peptidases are common virulence factors of other pathogens. T9SS-deficient mutants, such as ΔgldN and ΔporV, exhibit strong defects in secreted proteolytic activity. The F. columnare genome has many peptidase-encoding genes that may be involved in nutrient acquisition and/or virulence. Mutants lacking individual peptidase-encoding genes, or lacking up to ten peptidase-encoding genes, were constructed and examined for extracellular proteolytic activity, for growth defects, and for virulence in zebrafish and rainbow trout. Most of the mutants retained virulence, but a mutant lacking 10 peptidases, and a mutant lacking the single peptidase TspA exhibited decreased virulence in rainbow trout fry, suggesting that peptidases contribute to F. columnare virulence.

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.

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.

α-MSH is partially involved in the immunomodulation of Nile tilapia (Oreochromis niloticus) antibacterial immunity

α-Melanocyte-stimulating hormone (α-MSH) is a well-studied neuropeptide controlling skin and hair color. Besides, numerous immunomodulation roles of α-MSH were recorded in humans and mice. However, the regulatory effects of α-MSH in teleost immunity haven't been well elucidated. In this study, several precursor molecules of α-MSH (POMCs) and its receptors (MCRs) in Nile tilapia (Oreochromis niloticus) were characterized, and their expression characteristics and specific functions on antibacterial immunity were determined. Overall, POMCs and MCRs were principally detected in the brain, skin, and liver, and were remarkably promoted post Streptococcus agalactiae infection. However, tiny POMCs and MCRs were observed in tilapia immune organs (head kidney and spleen) or lymphocytes, and no evident immunomodulation effect was detected in vitro. Moreover, the in vivo challenge experiments revealed that α-MSH protects tilapia from bacterial infection by regulating responses in the brain and intestine. This study lays theoretical data for a deeper comprehension of the immunomodulation mechanisms of teleost α-MSH and the evolutional process of the vertebrate melanocortin system.

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.

Involvement and characterization of NLRCs and pyroptosis-related genes in Nile tilapia (Oreochromis niloticus) immune response

Pyroptosis is an inflammatory and programmed cell death initiated by the formation of the inflammasome, which consists of NLR, ASC, and Caspase. Pyroptosis has received growing attention due to its association with innate immunity and various diseases. However, the involvement and induction of the NLRCs and pyroptosis-related genes in fish immunity remain poorly studied. In this study, several NLRCs and pyroptosis-related genes in Nile tilapia (Oreochromis niloticus) were identified and characterized. Their involvement in bacterial infection and expression profiles in Nile tilapia lymphocyte responses were also assessed. Overall, three NLRC members (NOD1, NOD2, and NLRC3) and five pyroptosis-related genes (ASC1, Caspase1, Gsdme, NLRP3, and NLRP14) in Nile tilapia were cloned and characterized. The transcript levels of these molecules were broadly distributed in various tissues with comparatively high expression in the gills, intestine, and spleen. Their transcripts were also induced during Streptococcus agalactiae or Aeromonas hydrophila infection. Moreover, they were primarily expressed in T cells, NCCs, and Mo/Mφ and showed antibacterial and partially antiviral responses. The present study lays a theoretical foundation for further investigation of the pyroptosis mechanisms in fish as well as the evolution of the antiviral roles of pyroptosis in vertebrates.

Enhanced Hemolytic Activity of Mesophilic Aeromonas salmonicida SRW-OG1 Is Brought about by Elevated Temperatures

Aeromonas salmonicida is a well-known cold-water pathogenic bacterium. Previously, we reported the first isolation of pathogenic A. salmonicida from diseased Epinephelus coioides, a kind of warm-water fish, and it was proved to be a putative mesophilic strain with potent pathogenicity to humans. In order to investigate the mechanisms underlying mesophilic growth ability and virulence, the transcriptome of A. salmonicida SRW-OG1 at 18, 28, and 37 °C was analyzed. The transcriptome of A. salmonicida SRW-OG1 at different temperatures showed a clear separation boundary, which might provide valuable information for the temperature adaptation and virulence regulation of A. salmonicida SRW-OG1. Interestingly, aerA and hlyA, the hemolytic genes encoding aerolysin and hemolysin, were found to be significantly up-regulated at 28 and 37 °C. Since aerolysin and hemolysin are the most well-known and -characterized virulence factors of pathogenic Aeromonas strains, the induction of aerA and hlyA was associated with the mesophilic virulence. Further study proved that the extracellular products (ECPs) purchased from A. salmonicida SRW-OG1 cultured at 28 and 37 °C showed elevated hemolytic activity and virulence than those at 18 °C. Moreover, the silence of aerA and hlyA led to significantly decreased hemolysis and virulence. Taken together, our results revealed that the mesophilic virulence of A. salmonicida SRW-OG1 might be due to the enhanced expression of aerA and hlyA induced by elevated temperatures.

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.

Functional characterization of Vibrio alginolyticus T3SS regulator ExsA and evaluation of its mutant as a live attenuated vaccine candidate in zebrafish (Danio rerio) model

Vibrio alginolyticus, a Gram-negative bacterium, is an opportunistic pathogen of both marine animals and humans, resulting in significant losses in the aquaculture industry. Type III secretion system (T3SS) is a crucial virulence mechanism of V. alginolyticus. In this study, the T3SS regulatory gene exsA, which was cloned from V. alginolyticus wild-type strain HY9901, is 861 bp encoding a protein of 286 amino acids. The ΔexsA was constructed by homologous recombination and Overlap-PCR. Although there was no difference in growth between HY9901 and ΔexsA, the ΔexsA exhibited significantly decreased extracellular protease activity and biofilm formation. Besides, the ΔexsA showed a weakened swarming phenotype and an ~100-fold decrease in virulence to zebrafish. Antibiotic susceptibility testing showed the HY9901ΔexsA was more sensitive to kanamycin, minocycline, tetracycline, gentamicin, doxycycline and neomycin. Compared to HY9901 there were 541 up-regulated genes and 663 down-regulated genes in ΔexsA, screened by transcriptome sequencing. qRT-PCR and β-galactosidase reporter assays were used to analyze the transcription levels of hop gene revealing that exsA gene could facilitate the expression of hop gene. Finally, Danio rerio, vaccinated with ΔexsA through intramuscular injection, induced a relative percent survival (RPS) value of 66.7% after challenging with HY9901 wild type strain. qRT-PCR assays showed that vaccination with ΔexsA increased the expression of immune-related genes, including GATA-1, IL6, IgM, and TNF-α in zebrafish. In summary, these results demonstrate the importance of exsA in V. alginolyticus and provide a basis for further investigations into the virulence and infection mechanism.

Function of the rpoD gene in Pseudomonas plecoglossicida pathogenicity and Epinephelus coioides immune response

Pseudomonas plecoglossicida is a Gram-negative pathogenic bacterium that causes visceral white spot disease in several marine fish species, resulting in high mortality and financial loss. Based on previous RNA sequencing (RNA-seq) results, rpoD gene expression is significantly up-regulated in P. plecoglossicida during infection, indicating that rpoD may contribute to bacterial pathogenicity. To investigate the role of this gene, five specific short hairpin RNAs (shRNAs) were designed and synthesized based on the rpoD gene sequence, with all five mutants exhibiting a significant decrease in rpoD gene expression in P. plecoglossicida. The mutant with the highest silencing efficiency (89.2%) was chosen for further study. Compared with the wild-type (WT) P. plecoglossicida strain NZBD9, silencing rpoD in the rpoD-RNA interference (RNAi) strain resulted in a significant decrease in growth, motility, chemotaxis, adhesion, and biofilm formation in P. plecoglossicida. Silencing of rpoD also resulted in a 25% increase in the survival rate, a one-day delay in the onset of death, and a significant decrease in the number of white spots on the spleen surface of infected orange-spotted groupers (Epinephelus coioides). In addition, rpoD expression and pathogen load were significantly lower in the spleens of E. coioides infected with the rpoD-RNAi strain than with the WT strain of P. plecoglossicida. We performed RNA-seq of E. coioides spleens infected with different P. plecoglossicida strains. Results showed that rpoD silencing in P. plecoglossicida led to a significant change in the infected spleen transcriptomes. In addition, comparative transcriptome analysis showed that silencing rpoD caused significant changes in complement and coagulation cascades and the IL-17 signaling pathway. Thus, this study revealed the effects of the rpoD gene on P. plecoglossicida pathogenicity and identified the main pathway involved in the immune response of E. coioides.

UBE2G1 Is a Critical Component of Immune Response to the Infection of Pseudomonas Plecoglossicida in Large Yellow Croaker (Larimichthys crocea)

The large yellow croaker (Larimichthys crocea) is one of the most economically valuable mariculture fish in China. Infection of Pseudomonas plecoglossicida can cause an outbreak of “internal organ white-spot disease”, which seriously affects the aquaculture of the large yellow croaker. Ubiquitylation is closely related to the post-translation modification of proteins and plays a vital role in many hosts’ immune defense pathways, while the E2-binding enzyme is a key factor in ubiquitination. Our previous genome-wide association study found that the ubiquitin-binding enzyme E2G1 (designed LcUbe2g1) was one of the candidate genes related to disease resistance in large yellow croaker. In this study, we analyzed the molecular characteristics, function, and immune mechanism of the LcUbe2g1. The full-length cDNA is 812 bp, with an open reading frame of 513 bp, encoding 170 amino acid residues. The results of the RT-qPCR and immunohistochemistry analysis revealed that its transcription and translation were significantly activated by the infection of P. plecoglossicida in large yellow croaker. Immunocytochemistry experiments verified the co-localization of LcUBE2G1 and the ubiquitin proteins in the head kidney cells of large yellow croaker. Through GST pull-down, we found that LcUBE2G1 interacted with NEDD8 to co-regulate the ubiquitination process. The above results indicate that LcUBE2G1 is essential in the regulation of ubiquitination against P. plecoglossicida infection in large yellow croaker, which lays a foundation for further study on the resistance mechanism of internal organ white-spot disease.

Methods to Evaluate Bacterial Motility and Its Role in Bacterial–Host Interactions

Bacterial motility is a widespread characteristic that can provide several advantages for the cell, allowing it to move towards more favorable conditions and enabling host-associated processes such as colonization. There are different bacterial motility types, and their expression is highly regulated by the environmental conditions. Because of this, methods for studying motility under realistic experimental conditions are required. A wide variety of approaches have been developed to study bacterial motility. Here, we present the most common techniques and recent advances and discuss their strengths as well as their limitations. We classify them as macroscopic or microscopic and highlight the advantages of three-dimensional imaging in microscopic approaches. Lastly, we discuss methods suited for studying motility in bacterial–host interactions, including the use of the zebrafish model.

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.

Transcriptomic analysis using dual RNA sequencing revealed a Pathogen-Host interaction after Edwardsiella anguillarum infection in European eel (Anguilla anguilla)

Many studies have explored differentially expressed genes (DEGs) between some pathogens and hosts, but no study has focused on the interaction of DEGs between Edwardsiella anguillarum (Ea) and Anguilla anguilla (Aa). In this study, we examined the interactions of DEGs during Ea infection and Aa anti-infection processes by dual RNA sequencing. Total RNA from in vitro and in vivo (Aa liver) Ea culture was extracted. Using high-throughput transcriptomics, significant DEGs that were expressed between Ea cultured in vitro versus in vivo and those in the liver of the infected group versus control group were identified. Protein-protein interactions between the pathogen and host were explored using Cytoscape according to the HPIDB 3.0 interaction transcription database. The results showed that the liver in the infection group presented with severe bleeding and a large number of thrombi in the hepatic vessels. We found 490 upregulated and 398 downregulated DEGs of Ea in vivo versus Ea cultured in vitro, and 2177 upregulated and 970 downregulated genes in the liver of the infected eels. Gene ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) analysis of the pathogen DEGs revealed that the upregulated genes were mainly enriched in migration, colonization, biofilm formation, and significantly enriched in ABC transport and quorum sensing; the downregulated genes were mainly involved in metabolism, information transduction, organelle formation, enzyme catalysis, molecular transport, and binding. GO of the host DEGs showed that metabolic process, catalytic activity, single organism metabolic process, small molecule binding, nucleotide binding, nucleotide phosphate binding, and anion binding were markedly enriched. Finally, we found that 79 Ea and 148 Aa proteins encoded by these DEGs were involved in an interaction network, and some pathogen (DegP, gcvP, infC, carB, rpoC, trpD, sthA, and FhuB) and host proteins (MANBA, STAT1, ETS2, ZEP1, TKT1, NMI and RBPMS) appear to play crucial roles in infection. Thus, determining the interaction networks revealed crucial molecular mechanisms underlying the process of pathogenic infection and host anti-infection.

Dual RNA-Seq of Trunk Kidneys Extracted From Channel Catfish Infected With Yersinia ruckeri Reveals Novel Insights Into Host-Pathogen Interactions

Host-pathogen intectarions are complex, involving large dynamic changes in gene expression through the process of infection. These interactions are essential for understanding anti-infective immunity as well as pathogenesis. In this study, the host-pathogen interaction was analyzed using a model of acute infection where channel catfish were infected with Yersinia ruckeri. The infected fish showed signs of body surface hyperemia as well as hyperemia and swelling in the trunk kidney. Double RNA sequencing was performed on trunk kidneys extracted from infected channel catfish and transcriptome data was compared with data from uninfected trunk kidneys. Results revealed that the host-pathogen interaction was dynamically regulated and that the host-pathogen transcriptome fluctuated during infection. More specifically, these data revealed that the expression levels of immune genes involved in Cytokine-cytokine receptor interactions, the NF-kappa B signaling pathway, the JAK-STAT signaling pathway, Toll-like receptor signaling and other immune-related pathways were significantly upregulated. Y. ruckeri mainly promote pathogenesis through the flagellum gene fliC in channel catfish. The weighted gene co-expression network analysis (WGCNA) R package was used to reveal that the infection of catfish is closely related to metabolic pathways. This study contributes to the understanding of the host-pathogen interaction between channel catfish and Y. ruckeri, more specifically how catfish respond to infection through a transcriptional perspective and how this infection leads to enteric red mouth disease (ERM) in these fish.

Genome-Wide Identification and Expression Analysis of Potential Antiviral Tripartite Motif Proteins (TRIMs) in Grass Carp (Ctenopharyngodon idella)

Tripartite motif proteins (TRIMs), especially B30.2 domain-containing TRIMs (TRIMs-B30.2), are increasingly well known for their antiviral immune functions in mammals, while antiviral TRIMs are far from being identified in teleosts. In the present study, we identified a total of 42 CiTRIMs from the genome of grass carp, Ctenopharyngodon idella, an important cultured teleost in China, based on hmmsearch and SMART analysis. Among these CiTRIMs, the gene loci of 37 CiTRIMs were located on different chromosomes and shared gene collinearities with homologous counterparts from human and zebrafish genomes. They possessed intact conserved RBCC or RB domain assemblies at their N-termini and eight different domains, including the B30.2 domain, at their C-termini. A total of 19 TRIMs-B30.2 were identified, and most of them were clustered into a large branch of CiTRIMs in the dendrogram. Tissue expression analysis showed that 42 CiTRIMs were universally expressed in various grass carp tissues. A total of 11 significantly differentially expressed CiTRIMs were found in two sets of grass carp transcriptomes during grass carp reovirus (GCRV) infection. Three of them, including Cibtr40, CiTRIM103 and CiTRIM109, which all belonged to TRIMs-B30.2, were associated with the type I interferon response during GCRV infection by weighted network co-expression and gene expression trend analyses, suggesting their involvement in antiviral immunity. These findings may offer useful information for understanding the structure, evolution, and function of TRIMs in teleosts and provide potential antiviral immune molecule markers for grass carp.

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.

Histopathology, antioxidant responses, transcriptome and gene expression analysis in triangle sail mussel Hyriopsis cumingii after bacterial infection

Bacterial disease outbreaks in filter feeder bivalve Hyriopsis cumingii as water contamination become more frequent in the water ecosystem, especially in intensive aquaculture habitats. To characterize host-pathogen interactions between H. cumingii and bacterial infection, we investigated the effects of Stenotrophomonas maltophilia HOP3 and Aeromonas veronii GL1 on the antioxidant response, tissue invasion and transcriptome expression of H. cumingii by infectivity trials. We showed that bacterial infections resulted in tubular necrosis of the hepatopancreas and induced the acute immune response in H. cumingii. The transcriptomic study identified a total of 5957 differentially expressed genes (DEGs) after A. veronii challenge. These DEGs were implicated in 302 KEGG pathways, notably in Apoptosis, Phagosome and Lysosome. The results showed that the relative expressions of all six immune-related DEGs were effectively stimulated with A. veronii, accompanied by tissue differences. Overall, these findings will contribute to an analysis of the immune response of H. cumingii to bacterial infection at the transcriptomic level and its genomic resource for research.

Genetics Responses to Hypoxia and Reoxygenation Stress in Larimichthys crocea Revealed via Transcriptome Analysis and Weighted Gene Co-Expression Network

Simple Summary

Hypoxia, which occurs frequently in aquaculture, can cause serious harm to all aspects of the growth, reproduction and metabolism of cultured fish. Due to the intolerance of Larimichthys crocea to hypoxia, Larimichthys crocea often floats head or even dies under hypoxic environment. However, the molecular mechanism of hypoxia tolerance in Larimichthys crocea has not been fully described. Therefore, the aim of this study was to explore the hub regulatory genes under hypoxic stress environment by transcriptome analysis of three key tissues (liver, blood and gill) in Larimichthys crocea. We identified a number of important genes that exercise different regulatory functions. Overall, this study will provide important clues to the molecular mechanisms of hypoxia tolerance in Larimichthys crocea.

Abstract

The large yellow croaker (Larimichthys crocea) is an important marine economic fish in China; however, its intolerance to hypoxia causes widespread mortality. To understand the molecular mechanisms underlying hypoxia tolerance in L. crocea, the transcriptome gene expression profiling of three different tissues (blood, gills, and liver) of L. crocea exposed to hypoxia and reoxygenation stress were performed. In parallel, the gene relationships were investigated based on weighted gene co-expression network analysis (WGCNA). Accordingly, the Gene Ontology and Kyoto Encyclopedia of Genes and Genomes enrichment analysis showed that several pathways (e.g., energy metabolism, signal transduction, oxygen transport, and osmotic regulation) may be involved in the response of L. crocea to hypoxia and reoxygenation stress. In addition, also, four key modules (darkorange, magenta, saddlebrown, and darkolivegreen) that were highly relevant to the samples were identified by WGCNA. Furthermore, some hub genes within the association module, including RPS16, EDRF1, KCNK5, SNAT2, PFKL, GSK-3β, and PIK3CD, were found. This is the first study to report the co-expression patterns of a gene network after hypoxia stress in marine fish. The results provide new clues for further research on the molecular mechanisms underlying hypoxia tolerance in L. crocea.

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.

Exploring genetic resistance to infectious salmon anaemia virus in Atlantic salmon by genome-wide association and RNA sequencing

Background

Infectious Salmonid Anaemia Virus (ISAV) causes a notifiable disease that poses a large threat for Atlantic salmon (Salmo salar) aquaculture worldwide. There is no fully effective treatment or vaccine, and therefore selective breeding to increase resistance to ISAV is a promising avenue for disease prevention. Genomic selection and potentially genome editing can be applied to enhance host resistance, and these approaches benefit from improved knowledge of the genetic and functional basis of the target trait. The aim of this study was to characterise the genetic architecture of resistance to ISAV in a commercial Atlantic salmon population and study its underlying functional genomic basis using RNA Sequencing.

Results

A total of 2833 Atlantic salmon parr belonging to 194 families were exposed to ISAV in a cohabitation challenge in which cumulative mortality reached 63% over 55 days. A total of 1353 animals were genotyped using a 55 K SNP array, and the estimate of heritability for the trait of binary survival was 0.13–0.33 (pedigree-genomic). A genome-wide association analysis confirmed that resistance to ISAV was a polygenic trait, albeit a genomic region in chromosome Ssa13 was significantly associated with resistance and explained 3% of the genetic variance. RNA sequencing of the heart of 16 infected (7 and 14 days post infection) and 8 control fish highlighted 4927 and 2437 differentially expressed genes at 7 and 14 days post infection respectively. The complement and coagulation pathway was down-regulated in infected fish, while several metabolic pathways were up-regulated. The interferon pathway showed little evidence of up-regulation at 7 days post infection but was mildly activated at 14 days, suggesting a potential crosstalk between host and virus. Comparison of the transcriptomic response of fish with high and low breeding values for resistance highlighted TRIM25 as being up-regulated in resistant fish.

Conclusions

ISAV resistance shows moderate heritability with a polygenic architecture, but a significant QTL was detected on chromosome 13. A mild up-regulation of the interferon pathway characterises the response to the virus in heart samples from this population of Atlantic salmon, and candidate genes showing differential expression between samples with high and low breeding values for resistance were identified.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12864-021-07671-6.

Early immune response in large yellow croaker (Larimichthys crocea) after immunization with oral vaccine

Mesoporous silica nanoparticles (MSNs) have been used in the field of biomedicine as antigen carriers and adjuvants for protective antigens. In the present study, an oral nanovaccine against Vibrio alginolyticus was prepared employing MSNs as carriers. The uptake of the dihydrolipoamide dehydrogenase (DLDH) antigens in the intestine of large yellow croaker was evaluated using an immunohistochemistry assay. Additionally, the effects of the nanovaccine on the early immune response in large yellow croaker were investigated via oral vaccination. The presence of the antigens was detected in the mucosa and lamina propria of the foregut, midgut, and hindgut of large yellow croaker at 3 h following oral immunization. The expression levels of cytokines (i.e., lysozyme, IFN-γ, IFITM, TNF-α, IL-1β, IL-2, IL-4, IL-10, and IL-13) in the intestine, spleen, and head kidney tissues of large yellow croaker before and after the immune challenge were determined via RT-qPCR assay. The obtained results revealed that the expression levels of lysozyme, IFN-γ, IFITM, TNF-α, IL-1β, IL-2, IL-4, IL-10, and IL-13 in the intestine and head kidney of the vaccinated large yellow croaker, as well as the expression of lysozyme, IL-1β, and IL-10 in the spleen, exhibited time-dependent oscillation regulation patterns. Notably, the nanovaccine immunization could induce early (6 h) and high expression of IFN-γ in the spleen and kidney tissues after the bacterial infection. The current study supplements the available data on the early immune response to fish nanovaccines. It also provides a valuable theoretical basis for the future development of large yellow croaker oral vaccines.

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

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

Coordination of primary metabolism and virulence factors expression mediates the virulence of Vibrio parahaemolyticus towards cultured shrimp (Penaeus vannamei)

AIMS

Acute hepatopancreatic necrosis disease (AHPND) caused by Vibrio parahaemolyticus has emerged as a severe bacterial disease of cultured shrimp. To identify the key virulence factors, two AHPND-causing V. parahaemolyticus (Vp_AHPND ) strains (123 and 137) and two non-Vp_AHPND strains (HZ56 and ATCC 17082) were selected.

METHODS AND RESULTS

Challenge tests showed that the four strains exhibited different virulence towards shrimp with cumulative mortalities at 48 h postinfection (hpi) ranging from 10 to 92%. The expression of pirAB^VP in strain 123 and 137 was not significantly different. Genomic analysis revealed that the two Vp_AHPND strains contain a plasmid with the PirAB^VP toxins (pirAB^VP ) flanked by the insertion sequence (ISVal1) that has been identified in various locations of chromosomes in Vp_AHPND strains. The two Vp_AHPND strains possessed almost identical virulence factors, while ISVal1 disrupted three genes related to flagellar motility in strain 137. Phenotype assay showed that strain 123 possessed the highest growth rate and swimming motility, followed by strain 137, suggesting that the disruption of essential genes mediated by ISVal1 significantly affected the virulence level. Transcriptome analysis of two Vp_AHPND strains (123 and 137) further suggested that virulence genes related to the capsule, flagella and primary metabolism were highly expressed in strain 123.

CONCLUSIONS

Here for the first time, it is demonstrated that the virulence of Vp_AHPND is not only determined by the expression of pirAB^VP , but also is mediated by ISVal1 which affects the genes involved in flagellar motility and primary metabolism.

SIGNIFICANCE AND IMPACT OF THE STUDY

The genomic and transcriptomic analysis of Vp_AHPND strains provides valuable information on the virulence factors affecting the pathogenicity of Vp_AHPND.

Effects of oligochitosan on the growth, immune responses and gut microbes of tilapia (Oreochromis niloticus)

The immunomodulatory effects of oligochitosan have been demonstrated in several fish. However, the underlying mechanisms are not well characterized. The profound interplay between gut microbes and aquaculture has received much scientific attention but understanding the alternations of microbes populating in gut of tilapia (Oreochromis niloticus) fed with oligochitosan remains enigmatic. In this study, the effects of oligochitosan on the growth, immune responses and gut microbes of tilapia were investigated. The feeding trial was conducted in triplicates with the control diet supplemented with oligochitosan at different concentrations (0, 100, 200, 400 or 800 mg/kg). Following a six-week feeding trial, body weights of the fish supplemented with 200 mg/kg and 400 mg/kg oligochitosan were significantly higher than that of the control group. To address the immune responses stimulated by oligochitosan, by the quantitative real time PCR (qRT-PCR), the mRNA expression levels of CSF, IL-1β, IgM, TLR2 and TLR3 genes from head kidney were all significantly up-regulated in the 400 mg/kg group compared to the control. To characterize the gut microbes, bacterial samples were collected from the foregut, midgut, and hindgut, respectively and were subjected to high-throughput sequencing of 16S rDNA. The results showed that significantly lower abundance of Fusobacterium was detected in the hindgut of 400 mg/kg group compared to the control. Additionally, beta-diversity revealed that both gut habitat and oligochitosan had effects on the gut bacterial assembly. To further elucidate the mechanism underlying the effects of oligochitosan on bacterial assembly, the results showed that difference dosages of dietary oligochitosan could alter the specific metabolic pathways and functions of the discriminatory bacterial taxa, resulting in the different bacterial assemblies. To test the antibacterial ability of tilapia fed with oligochitosan, when the tilapias were challenged with Aeromonas hydrophila, the mortality of groups fed with dietary oligochitosan was significantly lower than that of the control. Taken together, appropriate dietary oligochitosan could improve growth, immune responses and alter the bacterial flora in the intestine of tilapia, so as to play a role in fighting against the bacterial infection.

Skin proteome profiling of tongue sole (Cynoglossus semilaevis) challenged with Vibrio vulnificus

Vibrio vulnificus is a major pathogen of cultured Cynoglossus semilaevis and results in skin ulceration and haemorrhage, but the proteomic mechanism of skin immunity against V. vulnificus remains unclear. In this study, we investigated the histopathology and skin immune response in C. semilaevis with V. vulnificus infection at the protein levels, the differential proteomic profiling of its skin was examined by using iTRAQ and LC-MS/MS analyses. A total of 951 proteins were identified in skin, in which 134 and 102 DEPs were screened at 12 and 36 hpi, respectively. Selected eleven immune-related DEPs (pvβ, Hsp71, MLC1, F2, α2ML, HCII, C3, C5, C8β, C9 and CD59) were verified for their immune roles in the V. vulnificus infection via using qRT-PCR assay. KEGG enrichment analysis revealed that most of the identified immune proteins were significantly associated with complement and coagulation cascades, antigen processing and presentation, salivary secretion and phagosome pathways. To our knowledge, this study is the first to describe the proteome response of C. semilaevis skin against V. vulnificus infection. The outcome of this study contributed to provide a new perspective for understanding the molecular mechanism of local skin mucosal immunity, and facilitating the development of novel mucosal vaccination strategies in fish.

Virulence factors impair epithelial junctions during bacterial infection

Epithelial cells are typically connected through different types of cell junctions that are localized from the apical membrane to the basal surface. In this way, epithelium cells form the first barrier against pathogenic microorganisms and prevent their entry into internal organs and the circulatory system. Recent studies demonstrate that bacterial pathogens disrupt epithelial cell junctions through targeting junctional proteins by secreted virulence factors. In this review, we discuss the diverse strategies used by common bacterial pathogens, including Pseudomonas aeruginosa, Helicobacter pylori, and enteropathogenic Escherichia coli, to disrupt epithelial cell junctions during infection. We also discuss the potential of targeting the pathogenic mechanisms in the treatment of pathogen-associated diseases.

Time-resolved RNA-seq provided a new understanding of intestinal immune response of European eel (Anguilla anguilla) following infection with Aeromonas hydrophila

No studies systematically examined the intestinal immune response for yellow stage of European eel (Anguilla anguilla) with Aeromonas hydrophila infection by time-resolved RNA-seq. Here, we examined transcriptional profiles of the intestines at three-time points following infection with A. hydrophila. Intraperitoneal injections caused mortalities within 48 h post-injection (hpi), with the survival rate 87.5% at 24 hpi and 83.9% at 48 hpi. The result from KEGG pathway enrichment analysis showed that the immune related "cytosolic DNA-sensing pathway" was significantly enriched at the first and second time points (6 hpi and 18 hpi), with the up-regulated expression of irf3, il1b, tnfaip3, cxcl8a, ap1-2, c-fos, polr3d, polr3g and polr3k both at 6 hpi and 18 hpi, but not at the third time point (36 hpi). According to the KEGG annotation, 326 immune and inflammation-related DEGs were found. The co-expression network of those 326 DEGs revealed the existence of three modules, and tlr1 was found to be in the center of the biggest module which contained massive DEGs from "signal transduction" and "transport and catabolism". The c3 isoforms showed different expression pattern among the three time points, indicating a unique activation of complement systems at 18 hpi. Furthermore, two cathelicidins (aaCATH_1 and aaCATH_2) were highly up-regulated at the first two time points, and the bacterial growth inhibition assay revealed their antibacterial properties against A. hydrophila. Our data indicated the important roles of cytosolic DNA-sensing pathway, as well as transcripts including tlr1, c3, polr and cathelicidins in the intestine of A. anguilla in response to A. hydrophila infection. The present study will provide leads for functional studies of host-pathogen interactions.

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.

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

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

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

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

Snakehead vesiculovirus (SHVV) infection alters striped snakehead (Ophicephalus striatus) cells (SSN-1) glutamine metabolism and apoptosis pathways

Snakehead vesiculovirus (SHVV) causes enormous economic losses in snakehead fish (Ophicephalus striatus) culture. Understanding replication mechanisms of virus is considerable significance in preventing and treating viral disease. In our previous studies, we have reported that glutamine starvation could significant inhibit the replication of SHVV. Furthermore, we also showed that SHVV infection could cause apoptosis of striped snakehead fish cells (SSN-1). However, the underlying mechanisms remain enigmatic. To decipher the relationships among the viral infection, glutamine starvation and apoptosis, SSN-1 cells transcriptomic profilings of SSN-1 cells infected with or without SHVV under glutamine deprived condition were analyzed. RNA-seq was used to identify differentially expressed genes (DEGs). Our data revealed that 1215 up-regulated and 226 down-regulated genes at 24 h post-infection were involved in MAPK, apoptosis, RIG-1-like and toll-like receptors pathways and glutamine metabolism. Subsequently, DEGs of glutamine metabolism and apoptosis pathways were selected to validate the sequencing data by quantitative real-time PCR (qRT-PCR). The expression patterns of both transcriptomic data and qRT-PCR were consistent. We observed that lack of glutamine alone could cause mild cellular apoptosis. However, lack of glutamine together with SHVV infection could synergistically enhance cellular apoptosis. When the cells were cultured in complete medium with glutamine, overexpression of glutaminase (GLS), an essential enzyme for glutamine metabolism, could significantly enhance the SHVV replication. While, SHVV replication was decreased in cells when GLS was knocked down by specific siRNA, indicating that glutamine metabolism was essential for viral replication. Furthermore, the expression level of caspase-3 and Bax was significantly decreased in SHVV infected cells with GLS overexpression. By contrast, they were significantly increased in SHVV infected cells with GLS silence by SiRNA, indicating that SHVV infection activated the Bax and caspase-3 pathways to induce apoptosis independent of glutamine. Our results reveal that SHVV replication and starvation of glutamine could synergistically promote the cellular apoptosis, which will pave a new way for developing strategies against the vial infection.