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

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

OBJECTIVE

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

METHODS

Published literatures of this pathogen have been reviewed.

RESULT

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

CONCLUSION

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

Transcriptomic analysis of Vibrio alginolyticus challenged by Rhizoma coptidis reveals mechanisms of virulence genes

Rhizoma coptidis, a Chinese herbal medicine widely used to treat various bacterial infections, has the potential to develop antibiotic substitutes to overcome the drug resistance of Vibrio alginolyticus. To study the inhibitory effect of R. coptidis on V. alginolyticus, we sequenced the transcriptomes of three groups of samples of wild-type V. alginolyticus (CK) and V. alginolyticus, which were stressed by 5 mg/mL R. coptidis for 2 h (RC_2 h) and 4 h (RC_4 h). CK was compared with RC_2 h and RC_4 h, respectively, and a total of 1565 differentially expressed genes (DEGs) (988 up-regulated and 577 down-regulated) and 1737 DEGs (1152 up-regulated and 585 down-regulated) were identified. Comparing RC_2 h with RC_4 h, 156 DEGs (114 up-regulated and 42 down-regulated) were identified. The ability of biofilm formation and motility of V. alginolyticus altered upon with different concentrations of R. coptidis. Interestingly, relative expression patterns of virulence genes appeared statistically significantly varied, upon different concentrations of R. coptidis extract. DEGs were annotated to the Gene Ontology (GO) database for function enrichment analysis and the Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway analysis, the results showed that the main enriched pathways, was those related to the virulence of V. alginolyticus. This study provides a new perspective for understanding the complex pathogenic mechanism of V. alginolyticus. R. coptidis could potnetially be used as alternative or complimnetary to antibiotics to treat infections after further research.

Genetic analysis of Vibrio alginolyticus challenged by Fructus schisandrae reveals the mechanism of virulence genes

Due to the abusive use of antibiotics, bacterial resistance has become a global problem and poses severe threats to aquaculture. The drug-resistant diseases caused by Vibrio alginolyticus have caused significant economic losses to cultured marine fish. Fructus schisandrae is used to treat inflammatory diseases in China and Japan. There have been no reports of bacterial molecular mechanisms associated with F. schisandrae stress. In this study, the inhibiting effect of F. schisandrae on the growth of V. alginolyticus was detected to understand response mechanisms at the molecular level. The antibacterial tests were analyzed via next-generation deep sequencing technology (RNA sequencing, RNA-seq). Wild V. alginolyticus (CK) was compared with V. alginolyticus, F. schisandrae incubated for 2 h, and V. alginolyticus, F. schisandrae incubated for 4 h. Our results revealed that there were 582 genes (236 upregulated and 346 downregulated) and 1068 genes (376 upregulated and 692 downregulated), respectively. Differentially expressed genes (DEGs) were involved in the following functional categories: metabolic process, single-organism process, catalytic activity, cellular process, binding, membrane, cell part, cell, and localization. FS_2 h was compared with FS_4 h, and 21 genes (14 upregulated and 7 downregulated) were obtained. The RNA-seq results were validated by detecting the expression levels of 13 genes using quantitative real-time polymerase chain reaction (qRT-PCR). The qRT-PCR results matched those of the sequencing, which reinforced the reliability of the RNA-seq. The results revealed the transcriptional response of V. alginolyticus to F. schisandrae, which will provide new ideas for studying V. alginolyticus' complex virulence molecular mechanism and the possibility of developing Schisandra to prevent and treat drug-resistant diseases.

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.

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.

Pneumolysin promotes host cell necroptosis and bacterial competence during pneumococcal meningitis as shown by whole-animal dual RNA-seq

Pneumolysin is a major virulence factor of Streptococcus pneumoniae that plays a key role in interaction with the host during invasive disease. How pneumolysin influences these dynamics between host and pathogen interaction during early phase of central nervous system infection in pneumococcal meningitis remains unclear. Using a whole-animal in vivo dual RNA sequencing (RNA-seq) approach, we identify pneumolysin-specific transcriptional responses in both S. pneumoniae and zebrafish (Danio rerio) during early pneumococcal meningitis. By functional enrichment analysis, we identify host pathways known to be activated by pneumolysin and discover the importance of necroptosis for host survival. Inhibition of this pathway using the drug GSK'872 increases host mortality during pneumococcal meningitis. On the pathogen's side, we show that pneumolysin-dependent competence activation is crucial for intra-host replication and virulence. Altogether, this study provides new insights into pneumolysin-specific transcriptional responses and identifies key pathways involved in pneumococcal meningitis.

Potential dsRNAs can be delivered to aquatic for defense pathogens

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

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.

Dual RNAseq highlights the kinetics of skin microbiome and fish host responsiveness to bacterial infection

Background

Tenacibaculum maritimum is a fish pathogen known for causing serious damage to a broad range of wild and farmed marine fish populations worldwide. The recently sequenced genome of T. maritimum strain NCIMB 2154^T provided unprecedented information on the possible molecular mechanisms involved in the virulence of this species. However, little is known about the dynamic of infection in vivo, and information is lacking on both the intrinsic host response (gene expression) and its associated microbiota. Here, we applied complementary omic approaches, including dual RNAseq and 16S rRNA gene metabarcoding sequencing using Nanopore and short-read Illumina technologies to unravel the host–pathogen interplay in an experimental infection system using the tropical fish Platax orbicularis as model.

Results

We showed that the infection of the host is characterised by an enhancement of functions associated with antibiotic and glucans catabolism functions but a reduction of sulfate assimilation process in T. maritimum. The fish host concurrently displays a large panel of immune effectors, notably involving innate response and triggering acute inflammatory response. In addition, our results suggest that fish activate an adaptive immune response visible through the stimulation of T-helper cells, Th17, with congruent reduction of Th2 and T-regulatory cells. Fish were, however, largely sensitive to infection, and less than 25% survived after 96 hpi. These surviving fish showed no evidence of stress (cortisol levels) or significant difference in microbiome diversity compared with controls at the same sampling time. The presence of T. maritimum in resistant fish skin and the total absence of any skin lesions suggest that these fish did not escape contact with the pathogen, but rather that some mechanisms prevented pathogens entry. In resistant individuals, we detected up-regulation of specific immune-related genes differentiating resistant individuals from controls at 96 hpi, which suggests a possible genomic basis of resistance, although no genetic variation in coding regions was found.

Conclusion

Here we focus in detail on the interplay between common fish pathogens and host immune response during experimental infection. We further highlight key actors of defence response, pathogenicity and possible genomic bases of fish resistance to T. maritimum.

Supplementary Information

The online version contains supplementary material available at 10.1186/s42523-021-00097-1.

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.

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.

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.

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

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

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.

Transcriptome profiling of avian pathogenic Escherichia coli and the mouse microvascular endothelial cell line bEnd.3 during interaction

Background

Avian pathogenic Escherichia coli (APEC), an important extraintestinal pathogenic E. coli, causes colibacillosis, an acute and mostly systemic disease involving multiple organ lesions such as meningitis. Meningitis-causing APEC can invade the host central nervous system by crossing the blood–brain barrier (BBB), which is a critical step in the development of meningitis. However, the bacteria-host interaction mechanism in this process remains unclear.

Methods

In this study, we examined E. coli and bEnd.3 cells transcriptomes during infection and mock infection to investigate the global transcriptional changes in both organisms using RNA sequencing approach.

Results

When APEC infected the bEnd.3 cells, several significant changes in the expression of genes related to cell junctional complexes, extracellular matrix degradation, actin cytoskeleton rearrangement, immune activation and the inflammatory response in bEnd.3 cells were observed as compared to the mock infection group. Thus, the immune activation of bEnd.3 cells indicated that APEC infection activated host defenses. Furthermore, APEC may exploit cell junction degradation to invade the BBB. In addition, amino acid metabolism and energy metabolism related genes were downregulated and the protein export pathway related genes were upregulated in APEC cultured with bEnd.3 cells, compared to that in control. Thus, APEC may encounter starvation and express virulence factors during incubation with bEnd.3 cells.

Conclusion

This study provides a comprehensive overview of transcriptomic changes that occur during APEC infection of bEnd.3 cells, and offers insights into the bacterial invasion strategies and the subsequent host defense mechanism.

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

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

Investigating the Role of FlhF Identifies Novel Interactions With Genes Involved in Flagellar Synthesis in Campylobacter jejuni

FlhF is a key protein required for complete flagellar synthesis, and its deletion results in the complete absence of a flagella and thus motility in Campylobacter jejuni. However, the specific mechanism still remains unknown. In this study, RNA-Seq, EMSAs, ChIP-qPCR and β-Galactosidase assays were performed to elucidate the novel interactions between FlhF and genes involved in flagellar synthesis. Results showed that FlhF has an overall influence on the transcription of flagellar genes with an flhF mutant displaying down-regulation of most flagellar related genes. FlhF can directly bind to the flgI promoter to regulate its expression, which has significant expression change in an flhF mutant. The possible binding site of FlhF to the flgI promoter was explored by continuously narrowing the flgI promoter region and performing further point mutations. Meanwhile, FlhF can directly bind to the promoters of rpoD, flgS, and fliA encoding early flagellin regulators, thereby directly or indirectly regulating the synthesis of class I, II, and III flagellar genes, respectively. Collectively, this study demonstrates that FlhF may directly regulate the transcription of flagellar genes by binding to their promoters as a transcriptional regulator, which will be helpful in understanding the mechanism of FlhF in flagellar biosynthetic and bacterial flagellation in general.

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.

An SNP-Based Genetic Map and QTL Mapping for Growth Traits in the Red-Spotted Grouper (Epinephelus akaara)

The red-spotted grouper (Epinephelus akaara) is one of the most commercially important aquatic species in China. However, its seedstock has low larval survival rates, and its stability is confronted with the danger of overexploitation. In this study, a high-density genetic map was constructed using 3435 single nucleotide polymorphisms (SNPs) from 142 first generation (F₁) full-sib offspring and two parents of a red-spotted grouper population. The total genetic length of the map was 2300.12 cM with an average intermarker distance of 0.67 cM. Seventeen genome-wide significant quantitative trait loci (QTLs) for growth-related traits were detected on 24 linkage groups, including 5 QTLs for full length, 7 QTLs for body length, and 5 QTLs for body weight. The contribution values of explained phenotypic variance ranged from 10.7% to 12.9%. Moreover, 13 potential candidate genes for growth-related traits were identified. Collectively, these findings will be useful for conducting marker-assisted selection of the red-spotted grouper in future studies.

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.