Time-resolved dual RNA-seq of tissue uncovers Pseudomonas plecoglossicida key virulence genes in host-pathogen interaction with Epinephelus coioides

A novel small non-coding RNA 562 mediates the virulence of Pseudomonas plecoglossicida by regulating the expression of fliP, a key component of flagella T3SS

Pseudomonas plecoglossicida is a vital pathogen that poses a substantial risk to aquaculture. Small RNAs (sRNAs) are non-coding regulatory molecules capable of sensing environmental changes and modulating virulence-associated signaling pathways, such as the assembly of flagella. However, the relevant researches on P. plecoglossicida are an urgent need. Here, we report a novel sRNA, sRNA562, which has potential to regulate the post-transcriptional of fliP, a key component of the lateral flagellar type III secretion system. In this study, the effects of sRNA562 on the virulence of P. plecoglossicida and its role in regulating the pathogenic process were investigated through the use of a constructed sRNA562 deletion strain. The deletion of sRNA562 resulted in an up-regulation of fliP in P. plecoglossicida, and leading to increased swarming motility and enhanced the ability of biofilm formation, adhesion and chemotaxis. Subsequent artificial infection experiment demonstrated that the deletion of sRNA562 increased the virulence of P. plecoglossicida towards hybrid grouper, as evidenced by a reduction in survival rate, elevation of tissue bacterial load, and the exacerbation of histopathological damage. Further studies have found that the deletion of sRNA562 lead to an up-regulation of fliP expression during hybrid grouper infection, thereby enhancing bacterial swarming ability and ultimately heightening pathogenicity, leading to a dysregulated host response to infection, tissue damage and eventually death. Our work revealed a sRNA that exerts negative regulation on the expression of lateral flagella in P. plecoglossicida, thereby impacting its virulence. These findings provide a new perspective on the virulence regulation mechanism of P. plecoglossicida, contributing to a more comprehensive understanding in the field of pathogenicity research.

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

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

Vibrio splendidus Fur regulates virulence gene expression, swarming motility, and biofilm formation, affecting its pathogenicity in Apostichopus japonicus

Vibrio splendidus is an opportunistic pathogen that causes skin ulcer syndrome and results in huge losses to the Apostichopus japonicus breeding industry. Ferric uptake regulator (Fur) is a global transcription factor that affects varieties of virulence-related functions in pathogenic bacteria. However, the role of the V. splendidus fur (Vsfur) gene in the pathogenesis of V. splendidus remains unclear. Hence, we constructed a Vsfur knock-down mutant of the V. splendidus strain (MTVs) to investigate the role of the gene in the effect of biofilm, swarming motility, and virulence on A. japonicus. The result showed that the growth curves of the wild-type V. splendidus strain (WTVs) and MTVs were almost consistent. Compared with WTVs, the significant increases in the transcription of the virulence-related gene Vshppd mRNA were 3.54- and 7.33-fold in MTVs at the OD₆₀₀ of 1.0 and 1.5, respectively. Similarly, compared with WTVs, the significant increases in the transcription of Vsm mRNA were 2.10- and 15.92-fold in MTVs at the OD₆₀₀ of 1.0 and 1.5, respectively. On the contrary, the mRNA level of the flagellum assembly gene Vsflic was downregulated 0.56-fold in MTVs at the OD₆₀₀ of 1.0 compared with the WTVs. MTVs caused delayed disease onset time and reduced A. japonicus mortality. The median lethal doses of WTVs and MTVs were 9.116 × 10⁶ and 1.658 × 10¹¹ CFU·ml^-1, respectively. Compared with WTVs, the colonization abilities of MTVs to the muscle, intestine, tentacle, and coelomic fluid of A. japonicus were significantly reduced. Correspondingly, the swarming motility and biofilm formation in normal and iron-replete conditions were remarkably decreased compared with those of WTVs. Overall, these results demonstrate that Vsfur contributes to the pathogenesis of V. splendidus by regulating virulence-related gene expression and affecting its swarming and biofilm formation abilities.

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.

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.

Prevalence of Bacterial Coinfections with Vibrio harveyi in the Industrialized Flow-through Aquaculture Systems in Hainan Province: A Neglected High-Risk Lethal Causative Agent to Hybrid Grouper

Vibrio harveyi is one of the most serious bacterial pathogens to aquatic animals worldwide. Evidence is mounting that coinfections caused by multiple pathogens are common in nature and can alter the severity of diseases in marine animals. However, bacterial coinfections involving V. harveyi have received little attention in mariculture. In this study, the results of pathogen isolation indicated that bacterial coinfection was a common and overlooked risk for hybrid groupers (♀ Epinephelus polyphekadion × ♂ E. fuscoguttatus) reared in an industrialized flow-through pattern in Hainan Province. The artificial infection in hybrid groupers revealed that coinfections with V. harveyi strain GDH11385 (a serious lethal causative agent to groupers) and other isolated pathogens resulted in higher mortality (46.67%) than infection with strain GDH11385 alone (33.33%), whereas no mortality was observed in single infection with other pathogens. Furthermore, the intestine, liver and spleen of hybrid groupers are target organs for bacterial coinfections involving V. harveyi. Based on the infection patterns found in this study, we propose that V. harveyi may have a specific spatiotemporal expression pattern of virulence genes when infecting the host. Taken together, bacterial coinfection with V. harveyi is a neglected high-risk lethal causative agent to hybrid groupers in the industrialized flow-through aquaculture systems in Hainan Province.

Metabolomics insights into the interaction between Pseudomonas plecoglossicida and Epinephelus coioides

As a highly infectious epidemic in aquaculture, Pseudomonas plecoglossicida infection results in high mortality of teleosts and serious economic losses. Host–pathogen interactions shape the outcome of an infection, yet we still understand little about the molecular mechanism of these pathogen-mediated processes. Here, a P. plecoglossicida strain (NZBD9) and Epinephelus coioides were investigated as a model system to characterize pathogen-induced host metabolic remodeling over the course of infection. We present a non-targeted metabolomics profiling of E. coioides spleens from uninfected E. coioides and those infected with wild-type and clpV-RNA interference (RNAi) strains. The most significant changes of E. coioides upon infection were associated with amino acids, lysophospatidylcholines, and unsaturated fatty acids, involving disturbances in host nutritional utilization and immune responses. Dihydrosphingosine and fatty acid 16:2 were screened as potential biomarkers for assessing P. plecoglossicida infection. The silencing of the P. plecoglossicida clpV gene significantly recovered the lipid metabolism of infected E. coioides. This comprehensive metabolomics study provides novel insights into how P. plecoglossicida shape host metabolism to support their survival and replication and highlights the potential of the virulence gene clpV in the treatment of P. plecoglossicida infection in aquaculture.

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.

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

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

Vibrio splendidus virulence to Apostichopus japonicus is mediated by hppD through glutamate metabolism and flagellum assembly

Vibrio splendidus is the main opportunistic pathogen that causes skin ulcer syndrome in Apostichopus japonicus. hppDIn the present study, mutant V. splendidus with an in-frame deletion of hppDV.s. (MTVs) was constructed. The median lethal doses of wild-type V. splendidus (WTVs) and MTVs were 5.129 × 10⁶ and 2.606 × 10¹⁰ CFU mL⁻¹, respectively. RNA-Seq was performed using WTVs and MTVs cells at different growth stages to explore the mechanisms of the pathogenesis mediated by hppDV.s. Gene Ontology analysis showed that the expression levels of 105 genes involved in amino acid metabolism and protein binding were remarkably different between MTVs and WTVs. Kyoto Encyclopedia of Genes and Genomes analysis showed that the pathways of glutamate metabolism and flagellum assembly involved in biofilm formation and swarming motility were suppressed in MTVs. Correspondingly, the swarming motility, biofilm formation and colonisation of MTVs were remarkably decreased compared with those of WTVs. The results showed that 4-hppD catalyses tyrosine into fumarate, which could enhance glutamate metabolism and ATP production; promote flagellum assembly through the TCA cycle and lead to higher swarming, biofilm formation and colonisation abilities, to contribute to the pathogenesis of V. splendidus.

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.

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.

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

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

Full-Length Transcriptome: A Reliable Alternative for Single-Cell RNA-Seq Analysis in the Spleen of Teleost Without Reference Genome

Fish is considered as a supreme model for clarifying the evolution and regulatory mechanism of vertebrate immunity. However, the knowledge of distinct immune cell populations in fish is still limited, and further development of techniques advancing the identification of fish immune cell populations and their functions are required. Single cell RNA-seq (scRNA-seq) has provided a new approach for effective in-depth identification and characterization of cell subpopulations. Current approaches for scRNA-seq data analysis usually rely on comparison with a reference genome and hence are not suited for samples without any reference genome, which is currently very common in fish research. Here, we present an alternative, i.e. scRNA-seq data analysis with a full-length transcriptome as a reference, and evaluate this approach on samples from Epinephelus coioides-a teleost without any published genome. We show that it reconstructs well most of the present transcripts in the scRNA-seq data achieving a sensitivity equivalent to approaches relying on genome alignments of related species. Based on cell heterogeneity and known markers, we characterized four cell types: T cells, B cells, monocytes/macrophages (Mo/MΦ) and NCC (non-specific cytotoxic cells). Further analysis indicated the presence of two subsets of Mo/MΦ including M1 and M2 type, as well as four subsets in B cells, i.e. mature B cells, immature B cells, pre B cells and early-pre B cells. Our research will provide new clues for understanding biological characteristics, development and function of immune cell populations of teleost. Furthermore, our approach provides a reliable alternative for scRNA-seq data analysis in teleost for which no reference genome is currently available.

Transcriptome analysis reveals the molecular mechanisms of response to an emergent yellow-flower disease in green Chinese prickly ash (Zanthoxylum schinifolium)

Chinese prickly ash (Zanthoxylum) is extensively used as spice and traditional medicine in eastern Asian countries. Recently, an emergent yellow-flower disease (YFD) break out in green Chinese prickly ash (Zanthoxylum schinifolium, Qinghuajiao in Chinese) at Chongqing municipality, and then leads to a sharp reduction in the yield of Qinghuajiao, and thus results in great economic losses for farmers. To address the molecular response for the emergent YFD of Qinghuajiao, we analyzed the transcriptome of 12 samples including the leaves and inflorescences of asymptomatic and symptomatic plants from three different towns at Chongqing by high-throughput RNA-Seq technique. A total of 126,550 genes and 229,643 transcripts were obtained, and 21,054 unigenes were expressed in all 12 samples. There were 56 and 164 different expressed genes (DEGs) for the AL_vs_SL (asymptomatic leaf vs symptomatic leaf) and AF_vs_SF (asymptomatic flower vs symptomatic flower) groups, respectively. The results of KEGG analysis showed that the “phenylpropanoid biosynthesis” pathway that related to plant–pathogen interaction were found in AL_vs_SL and AF_vs_SF groups, and the “Plant–pathogen interaction” found in AF_vs_SF group, implying that this Qinghuajiao YFD might cause by plant pathogen. Interestingly, we detected 33 common unigenes for the 2 groups, and almost these unigenes were up-regulated in the symptomatic plants. Moreover, most of which were homologs to virus RNA, the components of viruses, implying that this YFD was related to virus. Our results provided a primary molecular basis for the prevention and treatment of YFD of Qinghuajiao trees.

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

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

Analysis of Gene Expression Profiles, Cytokines, and Bacterial Loads Relevant to Alcoholic Liver Disease Mice Infected With V. vulnificus

Patients with liver disease are susceptible to infection with Vibrio vulnificus (V. vulnificus), but the specific reasons remain elusive. Through RNA-seq, we found that when mice with alcoholic liver disease (ALD) were infected with V. vulnificus by gavage, compared with the Pair group, the small intestinal genes affecting intestinal permeability were upregulated; and the number of differentially expressed genes related to immune functions (e.g., such as cell chemotaxis, leukocyte differentiation, and neutrophil degranulation) decreased in the liver, spleen, and blood. Further analysis showed that the number of white blood cells decreased in the Pair group, whereas those in the ALD mice did not change significantly. Interestingly, the blood bacterial load in the ALD mice was about 100 times higher than that of the Pair group. After the ALD mice were infected with V. vulnificus, the concentrations of T cell proliferation-promoting cytokines (IL-2, IL-23) decreased. Therefore, unlike the Pair group, ALD mice had weaker immune responses, lower T cell proliferation-promoting cytokines, and higher bacterial loads post-infection, possibly increasing their susceptibility to V. vulnificus infection. These new findings we presented here may help to advance the current understanding of the reasons why patients with liver disease are susceptible to V. vulnificus infection and provides potential targets for further investigation in the context of treatment options for V. vulnificus sepsis in liver disease patient.

Isolation and enzyme bioprospection of bacteria associated to Bruguiera cylindrica, a mangrove plant of North Sumatra, Indonesia

Mangrove-associated bacteria are of industrial interest due to their diverse and versatile enzyme properties. This study investigates the culturable bacteria from a wide range of habitat in a Bruguiera cylindrica mangrove ecosystem in North Sumatra. Screening of extracellular hydrolytic enzymes showed multiple potential traits in amylase, cellulase, chitinase, phosphatase, protease, and urease production by bacterial isolates. Molecular identification based on 16S rDNA region of a potential strain, Vibrio alginolyticus Jme3-20 is then reported as a newly proteolytic agent. The strain also showed a stable growth under salinity (NaCl) stress with considerable phosphate solubilization activities. Protease activity was enhanced by optimizing the 0.5 % (w/v) sucrose and soy peptone in the fermentation medium. SDS-PAGE and zymogram analysis showed the presence of a 35-kDa MW protease. Hence, our study revealed important insights into the bacterial diversity and activity in mangrove ecosystems, evidencing the importance of microbial exploration in this ecosystem.

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.

LncRNA-WAS and lncRNA-C8807 interact with miR-142a-3p to regulate the inflammatory response in grass carp

Septicemia of grass carp is a systemic inflammatory reaction caused by bacterial infection. More and more evidences show that long non-coding RNAs (lncRNAs) can participate in the regulation of inflammatory response. In the present study, lncRNA-WAS and lncRNA-C8807 were confirmed to be involved in the inflammatory response following infection with Aeromonas hydrophila. LncRNA-WAS and lncRNA-C8807 could interact with miR-142a-3p. LncRNA-WAS and lncRNA-C8807 interact with miR-142a-3p to effect pro-inflammatory genes and NF-κB pathway. Our results provide a theoretical basis for studying the molecular mechanism underlying the regulation of inflammation by lncRNA in grass carp.

CqsA/LuxS-HapR Quorum sensing circuit modulates type VI secretion system VflT6SS2 in Vibrio fluvialis

Vibrio fluvialis is an emerging enteric pathogen of increasing public health threat. Two quorum sensing (QS) systems, VfqI-VfqR and CqsA/LuxS-HapR, and two type VI secretion systems (T6SSs), VflT6SS1 and VflT6SS2, have been identified in V. fluvialis. Whether there exists any correlation between the two systems is unclear. In this study, we found that CqsA/LuxS-HapR circuit regulator LuxO represses while HapR activates VflT6SS2. The effect of LuxO is more pronounced at low cell density and is HapR-dependent. Deletion of hapR abolished Hcp expression and alleviated antibacterial virulence. However, these effects were rescued by HapR-expressing plasmid. Reporter fusion analyses showed that HapR is required for the promoter activities of VflT6SS2. Sequence inspection of the major cluster promoter revealed two potential Motif 1 HapR binding sites, and their bindings to HapR were confirmed by both electrophoretic mobility shift assay (EMSA) and DNase I footprinting assay. Meanwhile, two single Motif 2 sites were identified in tssD2_a (hcpA) and tssD2_b (hcpB) promoter regions of the orphan cluster which are less conserved and displayed lower affinities to HapR. Together, our study demonstrated that CqsA/LuxS-HapR QS manipulate VflT6SS2 in V. fluvialis, and this finding will enhance our understanding of possible crosstalk between T6SS and QS in microbes.

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.

Comparative genomics with a multidrug-resistant Klebsiella pneumoniae isolate reveals the panorama of unexplored diversity in Northeast Brazil

The emergence of community acquired infections increases the public health concern on K. pneumoniae and closely related bacteria among which antimicrobial resistance spreads. We report a multidrug-resistant K. pneumoniae isolate, B31, of a patient infected in the community and admitted to an intensive care unit in Northeast Brazil. Antimicrobial susceptibility and genome information were thoroughly investigated to characterize B31 in front of 172 sequenced strains of different countries. Assigned to the Sequence Type 15, which is globally spread, B31 presented extended spectrum beta-lactamase, tigecycline and ciprofloxacin resistance. Genome sequencing revealed most resistance genes being carried by plasmids with high dissemination potential. The absence of main virulence factors, like yersiniabactin and colibactin, apparently suggests a mild pathogenic strain which, on the contrary, persisted and caused severe infection in a previously healthy patient. The present study contributes to unveil the unclear genomic scenario of virulent and multidrug-resistant K. pneumoniae in Brazil.

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.

Gut-Liver Immune Response and Gut Microbiota Profiling Reveal the Pathogenic Mechanisms of Vibrio harveyi in Pearl Gentian Grouper (Epinephelus lanceolatus × E. fuscoguttatus )

Vibrio harveyi causes vibriosis in nearly 70% of grouper (Epinephelus sp.), seriously limiting grouper culture. As well as directly inhibiting pathogens, the gut microbiota plays critical roles in immune homeostasis and provides essential health benefits to its host. However, there is still little information about the variations in the immune response to V. harveyi infection and the gut microbiota of grouper. To understand the virulence mechanism of V. harveyi in the pearl gentian grouper, we investigated the variations in the pathological changes, immune responses, and gut bacterial communities of pearl gentian grouper after exposure to differently virulent V. harveyi strains. Obvious histopathological changes were detected in heart, kidney, and liver. In particular, nodules appeared and huge numbers of V. harveyi cells colonized the liver at 12 h postinfection (hpi) with highly virulent V. harveyi. Although no V. harveyi was detected in the gut, the infection simultaneously induced a gut-liver immune response. In particular, the expression of 8 genes associated with cellular immune processes, including genes encoding inflammatory cytokines and receptors, and pattern recognition proteins, was markedly induced by V. harveyi infection, especially with the highly virulent V. harveyi strain. V. harveyi infection also induced significant changes in gut bacterial community, in which Vibrio and Photobacterium increased but Bradyrhizobium, Lactobacillus, Blautia, and Faecalibaculum decreased in the group infected with the highly virulent strain, with accounting for 82.01% dissimilarity. Correspondingly, four bacterial functions related to bacterial pathogenesis were increased by infection with highly virulent V. harveyi, whereas functions involving metabolism and genetic information processing were reduced. These findings indicate that V. harveyi colonizes the liver and induces a gut-liver immune response that substantially disrupts the composition of and interspecies interactions in the bacterial community in fish gut, thereby altering the gut-microbiota-mediated functions and inducing fish death.

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

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

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.

Integration of RNA-seq and RNAi provides a novel insight into the immune responses of Epinephelus coioides to the impB gene of Pseudomonas plecoglossicida

Pseudomonas plecoglossicida is a Gram-negative bacterium that causes visceral white spot disease in Epinephelus coioides and leads to severe aquatic economic losses. The RNA-seq results of a previous study showed that the expression of the impB gene in P. plecoglossicida was significantly upregulated during infection. Four shRNAs were designed and synthesized to silence the impB gene in P. plecoglossicida, and the maximum silencing efficiency was 95.2%. Intraperitoneal injection of the impB-RNAi strain of P. plecoglossicida did not cause E. coioides death, and the spleens of infected fish did not show significant clinical symptoms. Although the injection of the mutant strain increased the antibody titer in E. coioides serum, it could not effectively protect E. coioides against wild strain infection. Compared with E. coioides infected with the wild type strain, the RNA-seq results for E. coioides infected with the impB-RNAi strain differed greatly. The KEGG enrichment analysis showed that key genes of the chemokine signalling pathway of E. coioides were downregulated by the silencing of impB in P. plecoglossicida. Infection with the impB-RNAi strain of P. plecoglossicida through injection did not produce good immune protection against E. coioides. The present study provides a novel insight into the immune responses of E. coioides to the impB gene of P. plecoglossicida.

Virulence of the Pathogen Porphyromonas gingivalis Is Controlled by the CRISPR-Cas Protein Cas3

Porphyromonas gingivalis is a key pathogen of periodontitis, a polymicrobial disease characterized by a chronic inflammation that destroys the tissues supporting the teeth. Thus, understanding the virulence potential of P. gingivalis is essential to maintaining a healthy oral microbiome. In nonoral organisms, CRISPR-Cas systems have been shown to modulate a variety of microbial processes, including protection from exogenous nucleic acids, and, more recently, have been implicated in bacterial virulence. Previously, our clinical findings identified activation of the CRISPR-Cas system in patient samples at the transition to disease; however, the mechanism of contribution to disease remained unknown. The importance of the present study resides in that it is becoming increasingly clear that CRISPR-associated proteins have broader functions than initially thought and that those functions now include their role in the virulence of periodontal pathogens. Studying a P. gingivaliscas3 mutant, we demonstrate that at least one of the CRISPR-Cas systems is involved in the regulation of virulence during infection.

The CRISPR (clustered regularly interspaced short palindromic repeat)-Cas system is a unique genomic entity that provides prokaryotic cells with adaptive and heritable immunity. Initial studies identified CRISPRs as central elements used by bacteria to protect against foreign nucleic acids; however, emerging evidence points to CRISPR involvement in bacterial virulence. The present study aimed to identify the participation of one CRISPR-Cas protein, Cas3, in the virulence of the oral pathogen Porphyromonas gingivalis, an organism highly associated with periodontitis. Our results show that compared to the wild type, a mutant with a deletion of the Cas3 gene, an essential nuclease part of the class 1 type I CRISPR-Cas system, increased the virulence of P. gingivalis. In vitro infection modeling revealed only mildly enhanced production of proinflammatory cytokines by THP-1 cells when infected with the mutant strain. Dual transcriptome sequencing (RNA-seq) analysis of infected THP-1 cells showed an increase in expression of genes associated with pathogenesis in response to Δcas3 mutant infection, with the target of Cas3 activities in neutrophil chemotaxis and gene silencing. The importance of cas3 in controlling virulence was corroborated in a Galleria mellonella infection model, where the presence of the Δcas3 mutant resulted in a statistically significant increase in mortality of G. mellonella. A time-series analysis of transcription patterning during infection showed that G. mellonella elicited very different immune responses to the wild-type and the Δcas3 mutant strains and revealed a rearrangement of association in coexpression networks. Together, these observations show for the first time that Cas3 plays a significant role in regulating the virulence of P. gingivalis.

IMPORTANCEPorphyromonas gingivalis is a key pathogen of periodontitis, a polymicrobial disease characterized by a chronic inflammation that destroys the tissues supporting the teeth. Thus, understanding the virulence potential of P. gingivalis is essential to maintaining a healthy oral microbiome. In nonoral organisms, CRISPR-Cas systems have been shown to modulate a variety of microbial processes, including protection from exogenous nucleic acids, and, more recently, have been implicated in bacterial virulence. Previously, our clinical findings identified activation of the CRISPR-Cas system in patient samples at the transition to disease; however, the mechanism of contribution to disease remained unknown. The importance of the present study resides in that it is becoming increasingly clear that CRISPR-associated proteins have broader functions than initially thought and that those functions now include their role in the virulence of periodontal pathogens. Studying a P. gingivaliscas3 mutant, we demonstrate that at least one of the CRISPR-Cas systems is involved in the regulation of virulence during infection.

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.

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.

MicroRNA negatively regulates NF-κB-mediated immune responses by targeting NOD1 in the teleost fish Miichthys miiuy

Inflammation is a self-protection mechanism that can be triggered when innate immune cells detect infection. Eradication of pathogen infection requires appropriate immune and inflammatory responses, but excessive inflammatory responses can cause uncontrolled inflammation, autoimmune diseases, or pathogen dissemination. Mounting evidence has shown that microRNAs (miRNAs) in mammals act as important and versatile regulators of innate immunity and inflammation. However, miRNA-mediated regulation networks are largely unknown in inflammatory responses in lower vertebrates. Here miR-144 and miR-217 are identified as negative regulators in teleost inflammatory responses. We find that Vibrio harveyi and lipopolysaccharide (LPS) treatment significantly upregulate the expression of fish miR-144 and miR-217. Upregulated miR-144 and miR-217 suppress LPS-induced inflammatory cytokine expression by targeting nucleotide-binding oligomerization domain-containing protein 1 (NOD1), thereby avoiding excessive inflammatory responses. In addition, miR-144 and miR-217 regulate inflammatory responses through NOD1-induced nuclear factor kappa (NF-kB) signaling pathways. These findings demonstrate that miR-144 and miR-217 play regulatory roles in inflammatory responses by modulating the NOD1-induced NF-κB signaling pathway.

Cold adaptation in drylands: transcriptomic insights into cold-stressed Nostoc flagelliforme and characterization of a hypothetical gene with cold and nitrogen stress tolerance

Environmental stressors, especially low temperature, are very common on the earth's dryland systems. Terrestrial cyanobacteria have evolved with cold adaptability in addition to extreme dryness and high irradiation resistance. The dryland soil surface-dwelling species, Nostoc flagelliforme, serves as a potential model organism to gain insights into cyanobacterial cold adaptation. In this study, we performed transcriptomic analysis of N. flagelliforme samples in response to low temperature. The results revealed that the biological processes, such as terpenoid biosynthetic process, oxidoreductase activity, carbohydrate metabolism, biosynthesis of secondary metabolites, lipid and nitrogen metabolism, were significantly and dynamically changed during the cold stress. It was noteworthy that the transcription of the denitrification pathway for ammonia accumulation was enhanced, implying an importance for nitrogen utilization in stress resistance. In addition, characterization of a cold-responsive hypothetical gene csrnf1 found that it could greatly improve the cold-resistant performance of cells when it was heterologously expressed in transgenic Nostoc sp. PCC 7120. It was also found that csrnf1 transgenic strain exhibited resistance to nitrogen-deficient environmental stress. Considering that dryland cyanobacteria have to cope with low temperature on infertile soils, this study would enrich our understanding on the importance of multifunction of the genes for environmental cold adaptation in drylands.

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.

Ayu (Plecoglossus altivelis) CD46 isoforms protect the cells from autologous complement attack

CD46 is an important immune regulatory receptor with multiple functions. However, studies on the function of teleost CD46, especially the different CD46 isoforms are limited. In this study, we identified three membrane cofactor protein (MCP, CD46) gene isoforms from ayu (Plecoglossus altivelis) and tentatively named as PaCD46 isoforms. PaCD46 isoforms were generated by alternative splicing and all consisted of four conserved short consensus repeats (SCRs), and the variable serine-threonine-proline-rich domain, transmembrane hydrophobic domain, and cytoplasmic tail. Phylogenetic analysis showed that the isoforms clustered together with other fish CD46 and then with higher animal CD46. Western blotting analysis of peripheral blood mononuclear cells (PBMC) revealed three bands, all of which had much larger molecular weights than the theoretical values of the three PaCD46 isoforms. Moreover, three PaCD46 isoforms were individually expressed on HEK293 cells, and Western blotting showed the similar band profile to that of PBMC. The recombinant extracellular domain of the PaCD46 isoforms, obtained by expression in Pichia pastoris, significantly reduced hemolysis activity of ayu sera. Furthermore, each of the three PaCD46 isoforms respectively protected the HEK293 cells expressing the isoform. The isoforms were also identified for their protection of autologous PBMC from complement activation. These results provided the first evidence that PaCD46 isoforms may be complement regulatory proteins to prevent complement-induced damage to self-tissue.

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.

Immune and gut bacterial successions of large yellow croaker (Larimichthys crocea) during Pseudomonas plecoglossicida infection

Large yellow croaker (Larimichthys crocea, LYC) aquaculture is being threatened by intensive infectious diseases. Relevant studies have focused on LYC immune responses to infection. By contrast, little is known how and to what extent the gut microbiota responds to infection. Here, we explored the interactions between LYC immune responses and gut bacterial communities during Pseudomonas plecoglossicida infection. P. plecoglossicida successfully colonized into LYC gut microbiota, resulting in an increasing mortality rate. Relative gene expressions of pro-inflammatory cytokines (TNF-α1, TNF-α2 and IL-1β) and anti-inflammatory cytokine (IL-10) were consistently and significantly induced by P. plecoglossicida infection, whereas non-specific immune enzymes activities were only enhanced at the early infection stages. P. plecoglossicida infection caused an irreversible disruption in the gut microbiota, of which infection and hours post infection constrained 16.2% and 5.6% variations, respectively. In addition, top 18 discriminatory taxa that were responsible for the difference between treatments were identified, whose abundances were significantly associated with the immune activities of LYC. Using a structural equation modeling (SEM), we found that gut bacterial communities were primarily governed by the conjointly direct (-0.33) and indirect (0) effects of infection, which subsequently affect host immune responses. Our results suggest that an irreversible dysbiosis in gut microbiota could be the causality of increasing mortality. To our knowledge, this is the first study to provide an integrated overview among pathogen infection, immune response and gut microbiota of LYC.

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

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