Transcriptome of Edwardsiella anguillarum in vivo and in vitro revealed two-component system, ABC transporter and flagellar assembly are three pathways pathogenic to European eel (Anguilla anguilla)

Comparative phenotype and transcriptome analysis revealed the role of ferric uptake regulator (Fur) in the virulence of Vibrio harveyi isolated from diseased American eel (Anguilla rostrata)

Vibrio harveyi is commonly found in salt and brackish water and is recognized as a serious bacterial pathogen in aquaculture worldwide. In this study, we cloned the ferric uptake regulator (fur) gene from V. harveyi wild-type strain HA_1, which was isolated from diseased American eels (Anguilla rostrata) and has a length of 450 bp, encoding 149 amino acids. Then, a mutant strain, HA_1-Δfur, was constructed through homologous recombination of a suicide plasmid (pCVD442). The HA_1-Δfur mutant exhibited weaker biofilm formation and swarming motility, and 18-fold decrease (5.5%) in virulence to the American eels; compared to the wild-type strain, the mutant strain showed time and diameter differences in growth and haemolysis, respectively. Additionally, the adhesion ability of the mutant strain was significantly decreased. Moreover, there were 15 different biochemical indicators observed between the two strains. Transcriptome analysis revealed that 875 genes were differentially expressed in the Δfur mutant, with 385 up-regulated and 490 down-regulated DEGs. GO and KEGG enrichment analysis revealed that, compared to the wild-type strain, the type II and type VI secretion systems (T2SS and T6SS), amino acid synthesis and transport and energy metabolism pathways were significantly down-regulated, but the ABC transporters and biosynthesis of siderophore group non-ribosomal peptides pathways were up-regulated in the Δfur strain. The qRT-PCR results further confirmed that DEGs responsible for amino acid transport and energy metabolism were positively regulated, but DEGs involved in iron acquisition were negatively regulated in the Δfur strain. These findings suggest that the virulence of the Δfur strain was significantly decreased, which is closely related to phenotype changing and gene transcript regulation.

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.

Bullfrogs (Lithobates catesbeianus) as a Potential Source of Foodborne Disease

In Mexico, bullfrogs (Lithobates catesbeianus) are produced as gourmet food. However, bullfrogs can be carriers of pathogens because the frogs' preferred living conditions occur in stagnant water. The present study aimed to identify bacteria that cause foodborne diseases or are associated with human diseases. For molecular identification, based on the sequential analysis by 16S rRNA or rpoD was conducted on all isolates obtained from bullfrog. A total of 91 bacterial isolates were obtained from bullfrogs; 14 genera and 23 species were identified, including Acinetobacter johnsonii 16.5%; Aeromonas media 14.3%; Aeromonas veronii 13.2%; Providencia rettgeri 7.7%; Citrobacter freundii 6.6%; Aeromonas caviae 4.4%; Aeromonas hydrophila and Elizabethkingia ursingii 3.3%; Pseudomonas stutzeri, Raoultella ornithinolytica, and Shewanella putrefaciens 2.2%; Acinetobacter guillouiae, Acinetobacter pseudolwoffii, Citrobacter portucalensis, Citrobacter werkmanii, Edwardsiella anguillarum, Klebsiella michiganensis, Kluyvera intermedia, Kocuria rosea, Myroides odoratimimus, Myroides odoratus, Proteus sp., and Proteus hauseri 1.1%. In this study, 49.4% of the isolates obtained cause foodborne disease, 19.8% are bacteria that play an important role in the spoilage of food, 5.5% of isolates have nosocomial significance, 13.2% of bacteria are considered to be pollutants of the ecosystem, and in the case of A. salmonicida and Edwardsiella anguillarum (12.1%) to have a negative impact on aquaculture. Acinetobacter pseudolwoffii and Citrobacter portucalensis have not been reported to cause disease. Lastly of these isolates, 97.8% (89/91) can cause disease by food consumption or by direct contact for immunocompromised persons. The presence of these bacteria in bullfrogs represents a significant problem for human health. There is evidence that these microorganisms are pathogenic and frogs may also be reservoirs.

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

Freund's complete adjuvant (FCA) and incomplete adjuvant (FIA), generally applied in subunit fishery vaccine, have not been explored on the molecular mechanism of the non-specific immune enhancement. As long noncoding RNAs (lncRNAs) play vital regulating roles in various biological activities, in this study, we examined the genome-wide expression of transcripts in the liver of European eel (Anguilla anguilla, Aa) inoculated with FCA and FIA (FCIA) to elucidate the regulators of lncRNAs in the process of Edwardsiella anguillarum (Ea) infection and Aa anti-Ea infection using strand-specific RNA-seq. After eels were challenged by Ea at 28 days post the first inoculation (dpi), compared to the control uninfected eels (Li group), the control infected eels (Con_Li group) showed severe bleeding, hepatocyte atrophy, and thrombi formed in the hepatic vessels of the liver, although eels inoculated with FCIA (FCIA_Li group) also formed slight thrombi in the hepatic vessels. Compared to the FCIA_Li group, there was about 10 times colony-forming unit (cfu) in the Con_Li group per 100 μg liver tissue, and the relative percent survival (RPS) of eels was 50% in FCIA_Li vs Con_Li. Using high-throughput transcriptomics, differential expressed genes (DEGs) and transcripts were identified and the results were verified using fluorescence real-time polymerase chain reaction (qRT-PCR). Interactions between the differential expressed lncRNAs (DE-lncRNAs) and the target DEGs were explored using Cytoscape according to their co-expression and co-location relationship. We found 13,499 lncRNAs (10,176 annotated and 3423 novel lncRNAs) between 3 comparisons of Con_Li vs Li, FCIA_Li vs Li, and FCIA_Li vs Con_Li, of which 111, 110, and 129 DE-lncRNAs were ascertained. Gene ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) analysis of DEGs targeted by DE-lncRNAs revealed these DEGs mainly involved in single-organism cellular process in BP, membrane in CC and binding in MF, and KEGG pathways showed that the target DEGs in co-expression and co-location enriched in cell adhesion molecules. Finally, 118 DE-lncRNAs target 1161 DEGs were involved in an interaction network of 8474 co-expression and 333 co-location-related links, of which 16 DE-lncRNAs play vital roles in anti-Ea infection. Taken together, the interaction networks revealed that DE-lncRNAs underlies the process of Ea infection and Aa anti-Ea infection.

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

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

Transcriptome RNA-seq revealed lncRNAs activated by Edwardsiella anguillarum post the immunization of OmpA protecting European eel (Anguilla anguilla) from being infected

Long noncoding RNAs (lncRNAs) play important roles in various biological activities as vital regulators. However, no study has focused on the lncRNA regulation of Outer membrane protein (OMP) immunization against aquatic bacterial infection. In this study, we examined the genome-wide expression of lncRNAs in the liver of European eel (Anguilla anguilla, Aa) administrated by a recombinant OmpA (rOmpA) from Edwardsiella anguillarum (Ea) to elucidate the functions of lncRNAs in the process of Ea infection and Aa anti-Ea infection using strand specific RNA-seq. Eels were challenged by Ea at 28 d post the immunization (dpi) of OmpA, and the result showed, compared to uninfected livers in the PBS group (Con group), the infected livers in the PBS group (Con_inf group) showed severe bleeding, hepatocyte atrophy and thrombi formed in the hepatic vessels; livers in the OmpA group (OmpA_inf) also formed slight thrombi in the hepatic vessels. The relative percent survival of eels in OmpA_inf vs Con_inf was 78.6%. Using high-throughput transcriptomics, we found 13405 lncRNAs in 3 compares of Con_inf vs Con, OmpA_inf vs Con and OmpA_inf vs Con_inf, of which 111, 129 and 158 DE-lncRNAs were ascertained. GO analysis of the DE-lncRNAs revealed the targeting DEGs were mainly involved in single-organism process, signaling, biological process and response to stimulus in BP, component of membrane in CC and binding in MF; KEGG pathways showed that the targeting DEGs in co-expression and co-location enriched in cell adhesion molecules. Finally, 54 DE-lncRNAs targeting 1675 DEGs were involved in an interaction network of 21692 co-expression and 483 co-location related links, of which 18 DE-lncRNAs appear to play crucial roles in anti-Ea infection. Thus, the interaction networks revealed crucial DE-lncRNAs underlying the process of Ea infection and Aa anti-Ea infection pre and post the immunization of OmpA.