Comprehensive Relationship Analysis of the Long Noncoding RNAs (lncRNAs) and the Target mRNAs in Response to the Infection of Edwardsiella anguillarum in European eel (Anguilla anguilla) Inoculated with Freund's Adjuvant

Correlation Between Effector Gene Expression Targeted by lncRNAs in the Oomycete Fish Pathogen, Saprolegnia parasitica

Saprolegniasis caused by Saprolegnia parasitica leads to significant economic losses in the aquaculture industry worldwide. Effector proteins secreted by pathogens are key molecules involved in their pathogenicity and long non-coding lncRNAs (lncRNAs) act as regulators in these processes. However, little is known about the lncRNAs and effector proteins in S. parasitica. Here, we first identified 1027 lncRNAs during the developmental stages and infection process of S. parasitica. Compared with mRNAs, these lncRNAs had shorter sequences and exon lengths and lower expression levels. In addition, their sequence conservation among other oomycete species was also low. The S. parasitica lncRNAs were characterized according to developmental stage and infection time point. We also identified effector proteins using a computational pipeline. In total, 131 S. parasitica effector proteins were identified and classified into 34 families. The 47 genes encoding effector genes were neighbors of 39 lncRNAs, and there was a correlation between the transcription level of lncRNAs and their neighboring genes. Gain- and loss-of-function experiments revealed that lncRNA8375.2 promoted the expression of a neighboring effector gene, SpCAP. Our results provide new data on S. parasitica lncRNAs and effector proteins, and provide insights into the lncRNA-effector module involved in S. parasitica.

Differentially Expressed Genes and Alternative Splicing Analysis Revealed the Difference in Virulence to American Eels (Anguilla rostrata) Infected by Edwardsiella anguillarum and Aeromonas hydrophila

Edwardsiella anguillarum and Aeromonas hydrophila are two common bacterial pathogens affecting cultivated eels, and the differences in their virulence remain unclear. In this study, after two groups of American eels (Anguilla rostrata) were administered the LD50 dose of E. anguillarum and A. hydrophila, respectively, the histopathology of the liver, trunk kidney, and spleen, as well as transcriptomic RNA sequencing (RNA-seq) analysis of the spleen, was examined at three time points: pre-infection (Con group) and post-infection at 36 h (Ea_36 group, Ah_36 group) and 60 h (Ea_60 group, Ah_60 group). The results showed that the differences in pathological changes were characterized by severe hepatocyte edema at 36 h post-infection (hpi) and hepatocyte atrophy at 60 hpi in the livers of eels infected by A. hydrophila, in contrast to the severe atrophy of glomeruli in the trunk kidneys and numerous bacterial nodules in the spleens of eels infected by E. anguillarum. The RNA-seq results revealed 906 and 77 typical differentially expressed genes (DEGs) in eels infected with E. anguillarum and A. hydrophila, respectively, compared to the control eels. The DEGs between the infected and control groups were predominantly annotated in GO terms related to binding, catalytic activity, membrane part, cell part, and cellular process, as well as in KEGG pathways associated with human diseases and organismal systems. The GO enrichment analysis showed 83 and 146 differential GO terms, along with 32 and 78 differential KEGG pathways in two comparisons of Ea_36 vs Con versus Ah_36 vs Con and Ea_60 vs Con versus Ah_60 vs Con, respectively. Furthermore, the analysis of differential alternative splicing genes (DASs) showed 1244 and 1341 DASs out of 12,907 and 12,833 AS genes, respectively, in the comparisons of Ea_36 vs Ah_36 and Ea_60 vs Ah_60. These DASs were enriched in two common KEGG pathways: "NOD-like receptor signaling pathway" and "necroptosis" which shared 11 hub DASs. Finally, analysis of protein-protein interactions revealed that 91 of 412 cross DASs between Ea_36 vs Ah_36 and Ea_60 vs Ah_60 potentially play an essential role in the difference in virulence of E. anguillarum and A. hydrophila in American eels, with 12 encoded proteins being particularly notable. Together, this study is the first to report a comparative pathogenicity and RNA-seq analysis of E. anguillarum and A. hydrophila in American eels, shedding new light on our understanding of the differences in virulence as revealed by pathological changes, DEGs, and DASs, contributing to more effective control strategies to prevent outbreaks of bacterial infections.

Transcriptome analysis of host anti-Aeromonas hydrophila infection revealed the pathogenicity of A. hydrophila to American eels (Anguilla rostrata)

Aeromonas hydrophila is a commonly pathogenic bacterium in cultivated eels, but its pathogenicity to American eel (Anguilla rostrata) and the molecular mechanism of host anti-A. hydrophila infection remains uncertain. In this study, LD50 of A. hydrophila to American eels was determined and bacterial load in the liver and kidney of eels was assessed post 2.56 doses of LD50 of A. hydrophila infection. The results showed that the LD50 of A. hydrophila to American eels was determined to be 3.9 × 105 cfu/g body weight (7.8 × 106 cfu/fish), and the bacterial load peaked at 36 h post the infection (hpi) in the liver. Then, the histopathology was highlighted by congestion in splenic blood vessels, atrophied glomeruli, and necrotic hepatocytes. Additionally, the results of qRT-PCR revealed that 18 host immune-related genes showed significantly up or downregulated post-infection compare to that of pre-infection. Finally, results of the RNA-seq revealed 10 hub DEGs and 7 encoded proteins play essential role to the anti-A. hydrophila infection in American eels. Pathogenicity of A. hydrophila to American eels and RNA-seq of host anti-A. hydrophila infection were firstly reported in this study, shedding new light on our understanding of the A. hydrophila pathogenesis and the host immune response to the A. hydrophila infection strategies in gene transcript.

Transcriptome Analysis of Host Anti-Vibrio harveyi Infection Revealed the Pathogenicity of V. harveyi to American Eel (Anguilla rostrata)

Vibrio harveyi, a recently discovered pathogenic bacterium isolated from American eels (Anguilla rostrata), poses uncertainties regarding its pathogenesis in American eel and the molecular mechanisms underlying host defense against V. harveyi infection. This study aimed to determine the LD50 of V. harveyi in American eel and assess the bacterial load in the liver, spleen, and kidney post-infection with the LD50 dose. The results showed that the LD50 of V. harveyi via intraperitoneal injection in American eels over a 14d period was determined to be 1.24 × 103 cfu/g body weight (6.2 × 104 cfu/fish). The peak bacterial load occurred at 36 h post-infection (hpi) in all three organs examined. Histopathology analysis revealed hepatic vein congestion and thrombi, tubular vacuolar degeneration, and splenic bleeding. Moreover, quantitative reverse transcription polymerase chain reaction (qRT-PCR) results indicated significant up or downregulation of 18 host immune- or anti-infection-related genes post 12 to 60 hpi following the infection. Additionally, RNA sequencing (RNA-seq) unveiled 7 hub differentially expressed genes (DEGs) and 11 encoded proteins play crucial roles in the anti-V. harveyi response in American eels. This study firstly represents the comprehensive report on the pathogenicity of V. harveyi to American eels and RNA-seq of host's response to V. harveyi infection. These findings provide valuable insights into V. harveyi pathogenesis and the strategies employed by the host's immune system at the transcriptomic level to combat V. harveyi infection.

RNA-seq analysis revealed the pathogenicity of Vibrio vulnificus to American eel (Anguilla rostrata) and the strategy of host anti-V. vulnificus infection

Vibrio vulnificus is a commonly pathogenic bacterium in cultivated eels, but its pathogenicity to American eel (Anguilla rostrata) and the molecular mechanism of host anti-V. vulnificus infection remains uncertain. In this study, American eels were infected with different dose of V. vulnificus to determine the LD50. Then, bacterial load in the liver and kidney histopathology were assessed post the LD50 of V. vulnificus infection. Additionally, gene expressions of 18 immune related genes in the liver, spleen and kidney were detected. Furthermore, transcriptome sequencing and enrichment of differentially expressed genes (DEGs) were analyzed in the eel spleens between pre-infection (Con_0), post-36 h (Vv_36), and post-60 h (Vv_60) infection. The results showed that LD50 of V. vulnificus to American eels was determined to be 5.0 × 105 cfu/g body weight, and the bacterial load peaked at 24 and 12 h post the infection (hpi) in the kidney and liver, respectively. The histopathology was highlighted by necrotic hepatocytes and splenic cells, congestion blood vessels in liver and spleen, atrophied glomeruli and vacuolization of renal tubular epithelial cells. The results of RT-PCR revealed that 18 host immune-related genes showed significantly up or downregulated expression post-infection compare to that of pre-infection. Finally, results of the RNA-seq revealed 16 DEGs play essential role to the immunosuppression in American eels, and the protein-protein interactions shed light on the widespread upregulation GEGs related to metabolism and immune response maintained the host cell homeostasis post the V. vulnificus infection, shedding new light on our understanding of the V. vulnificus pathogenesis towards understudied American eel and the host anti-V. vulnificus infection strategies in gene transcript.

Pathogenicity of Edwardsiella anguillarum to American eels (Anguilla rostrata) and RNA-seq analysis of host immune response to the E. anguillarum infection

Edwardsiella anguillarum is a commonly pathogenic bacterium in cultivated eels, but its pathogenicity to American eel (Anguilla rostrata) and the molecular mechanism of host anti-E. anguillarum infection remains uncertain. In this study, LD50 of E. anguillarum to American eels was determined and bacterial load in the liver and kidney of eels was assessed post the LD50 of E. anguillarum infection. The results showed that LD50 of E. anguillarum to American eels was determined to be 2.5 × 105 cfu/g body weight, and the bacterial load peaked at 36 and 72 h post the infection (hpi) in the kidney and liver, respectively. Then, the histopathology was highlighted by congestion in splenic blood vessels, atrophied glomeruli, and necrotic hepatocytes, as well as ultrastructural pathology in the kidney were charactered by acute nephritis, showing necrosis of the renal tubular epithelial cells, glomerular capillaries dilate, mitochondria swelling and ribosomes separate from the endoplasmic reticulum. Furthermore, the results of qRT-PCR revealed that 12 host immune-related genes showed significantly up or downregulated post-infection compare to that of pre-infection. Finally, results of the RNA-seq revealed 6 hub DEGs play essential role to the anti-E. anguillarum infection in American eels. Pathogenicity of E. anguillarum to American eels and hub genes related host anti- E. anguillarum infection were firstly reported in this study, shedding new light on our understanding of the E. anguillarum pathogenesis and the host immune response to the E. anguillarum infection strategies in gene transcript.

Analysis of Alternative Splicing and Long Noncoding RNAs After the Edwardsiella anguillarum Infected the Immunized European Eels (Anguilla anguilla) Revealed the Role of Outer Membrane Protein A in OmpA Subunit Vaccine

Edwardsiella anguillarum is a bacterium that commonly infects cultivated eels. Outer membrane protein A (OmpA) emulsified with Freund's adjuvant has been shown to be an effective fishery vaccine against this pathogen. However, the specific roles of OmpA in the vaccine have not been fully explored. In this study, we performed RNA-seq in the liver of a European eel (Anguilla anguilla) after challenge with E. anguillarum in eels previously immunized with an OmpA subunit vaccine. Our aim was to elucidate the differentially alternative splicing (DAS) and differentially expressed long noncoding RNAs (DE-lncRNAs) using a genome-wide transcriptome. The results showed after that at 28 days post-immunization, eels challenged with E. anguillarum (Con_inf) exhibited severe pathological changes in the liver. In contrast, the OmpA infused eels (OmpA_inf group) showed infiltrated lymphocytes, while Freund's adjuvant-inoculated eels (FCIA_inf group) showed edema of hepatocytes and blood coagulation. The relative percent survival (RPS) was 77.7% and 44.4% for OmpA_inf and FCIA_inf compared to the Con_inf group. We identified 37 DE-lncRNAs and 293 DAS genes between OmpA_inf and FCIA_inf. Interactions between DAS gene-expressed proteins indicated that 66 expressed proteins formed 20 networks. Additionally, 33 DE-lncRNAs interacted with 194 target genes formed 246 and 41 networks in co-expression and co-location. Taken together, our findings demonstrate that the OmpA subunit vaccine elicits a higher RPS and provides novel insights into the role of OmpA through DAS genes and DE-lncRNAs perspective. These results are significant for the development of fishery subunit vaccines.

Inoculation of Freund's adjuvant in European eel (Anguilla anguilla) revealed key KEGG pathways and DEGs of host anti-Edwardsiella anguillarum infection

Freund's complete (FCA) and incomplete adjuvants (FIA), generally applied in subunit fishery vaccine, have not been explored on the molecular mechanism of the nonspecific immune enhancement. In this study, we examined the RNA-seq in the spleen of European eel (Anguilla anguilla) inoculated with FCA and FIA (FCIA group) to elucidate the key KEGG pathways and differential expressed genes (DEGs) in the process of Edwardsiella anguillarum infection and A. anguilla anti-E. anguillarum infection using genome-wide transcriptome. After eels were challenged by E. anguillarum at 28 d post the first inoculation (dpi), compared to the control uninfected eels (Con group), the control infected eels (Con_inf group) showed severe pathological changes in the liver, kidney and spleen, although infected eels post the inoculation of FCIA (FCIA_inf group) also formed slight bleeding. Compared to the FCIA_inf group, there was more than 10 times colony forming unit (cfu) in the Con_inf group per 100 μg spleen, kidney or blood, and the relative percent survival (RPS) of eels was 44.4% in FCIA_inf vs Con_inf. Compared to the Con group, the SOD activity in the FCIA group increased significantly in the liver and spleen. Using high-throughput transcriptomics, DEGs were identified and 29 genes were verified using fluorescence real-time polymerase chain reaction (qRT-PCR). The result of DEGs clustering showed 9 samples in 3 groups of Con, FCIA and FCIA_inf were similar, contrast to distinct differences of 3 samples in the Con_inf group. We found 3795 up and 3548 down regulated DEGs in the compare of FCIA_inf vs Con_inf, of which 5 enriched KEGG pathways of "Lysosome", "Autophagy", "Apoptosis", "C-type lectin receptor signaling" and "Insulin signaling" were ascertained, and 26 of 30 top GO terms in the compare were significantly enriched. Finally, protein-protein interactions between the DEGs of the 5 KEGG pathways and other DEGs were explored using Cytoscape 3.9.1. The compare of FCIA_inf vs Con_inf showed 110 DEGs from the 5 pathways and 718 DEGs from other pathways formed total of 9747° in a network, of which 9 hub DEGs play vital roles in anti-infection or apoptosis. Together, the interaction networks revealed that 9 DEGs involved in the 5 pathways underlies the key process of A. anguilla anti-E. anguillarum infection or host cell apoptosis.