Transcriptome analysis of the spleen provides insight into the immunoregulation of Mastacembelus armatus under Aeromonas veronii infection

Transcriptomics reveals that NAD(P)H affects the development of the Zig-zag eel (Mastacembelus armatus ♀) × Spiny eel (Sinobdella sinensis ♂) hybrid offspring leading to low hatching rates

Zig-zag eel (Mastacembelus armatus (2 n = 48)) and Spiny eel (Sinobdella sinensis (2 n = 48)) are two species of the Mastacembelidae family commonly found in southern China. Hybridization between the two has a very high deformity rate and a very low hatching rate. In order to investigate the reasons for this, the first hybridization between M. armatus and S. sinensis was carried out using artificial insemination, and the embryonic development of the hybrid offspring was examined using microphotography, and the malformations of the hybrid offspring were investigated by transcriptomics. The experiments showed that the average egg production was 4265.7 ± 322.94 (Mean ± SD), the average fertilization rate of hybrid offspring was 98.67 ± 0.58 % (Mean ± SD), the hatching rate was 12.06 ± 3.44 % (Mean ± SD), the deformity rate was 98.15 ± 3.21 % (Mean ± SD), and the embryonic development successively went through the five main stages of fertilized egg, egg cleavage, embryo formation, organogenesis, and exertion of membranes. Transcriptomics showed that the expression of NAD(P)H-related enzyme activity DEGs was increased, and many DEGs related to cell signaling molecule transmission and metabolic regulation are enriched in KEGG pathways, such as IL-17 signaling pathway, Osteoclast differentiation, TNF signaling pathway and MAPK signaling pathway. etc. The major types of DEGs corresponded to those coding for proteins. This study suggests that the high malformation rate in hybrid offspring may be caused by impaired synthesis of proteins during embryonic development.

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

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

Transcriptomic insights into the immune response of the intestine to Aeromonas veronii infection in northern snakehead (Channa argus)

Intestinal inflammation is a protective response that is implicated in bacterial enteritis triggered by gastrointestinal infection. The immune mechanisms elicited in teleost against the infection of Aeromonas veronii are largely unknown. In this study, we performed a de novo northern snakehead (Channa argus) transcriptome assembly using Illumina sequencing platform. On this basis we performed a comparative transcriptomic analysis of northern snakehead intestine from A. veronii-challenge and phosphate buffer solution (PBS)-challenge fish, and 2076 genes were up-regulated and 1598 genes were down-regulated in the intestines infected with A. veronii. The Gene Ontology (GO) enrichment analysis indicated that the differentially expressed genes (DEGs) were enriched to 27, 21 and 20 GO terms in biological process, cellular component, and molecular function, respectively. A Kyoto Encyclopedia of Genes and Genomes (KEGG) analysis showed that 420 DEGs were involved in 194 pathways. Moreover, 33 DEGs were selected for quantitative real-time PCR analysis to validate the RNA-seq data. The results reflected the consistency of the expression levels between qRT-PCR and RNA-seq data. In addition, a time-course analysis of the mRNA expression of 33 immune-related genes further indicated that the intestinal inflammation to A. veronii infection simultaneously regulated gene expression alterations. The present study provides transcriptome data of the teleost intestine, allowing us to understand the mechanisms of intestinal inflammation triggered by bacterial pathogens. DATA AVAILABILITY STATEMENT: All data supporting the findings of this study are available within the article and Supplementary files. The RNA-seq raw sequence data are available in NCBI short read archive (SRA) database under accession number PRJNA615958.

Characterization and expression analysis of the interferon regulatory factor (IRF) gene family in zig-zag eel (Mastacembelus armatus) against Aeromonas veronii infection

Interferon regulatory factors (IRFs) play an important role in innate and adaptive immune system. However, in teleosts, the data on IRFs is still scarce. Here, for the first time, we identified 11 members of IRFs from the zig-zag eel Mastacembelus armatus (MarIRF1-10). The deduced protein sequences are highly conserved among different fish species especially in DBD and IAD domain. Phylogenetic analysis indicated that MarIRFs preferentially grouped with fish species in Synbranchiformes or Perciformes. Expression analysis showed that MarIRFs were expressed in all nine tissues including spleen, gill, muscle and intestine. After infected by Aeromonas veronii, expression of MarIRF2, MaIRF4b and MaIRF5 were significantly upregulated in spleen, MarIRF1, MarIRF2 were significantly upregulated in kidney, but in liver, nearly all MarIRFs were downregulated. Taken together, this study first reported molecular characterization and expression patterns of 11 IRFs in the zig-zag eel. All these results will contribute a lot to better understanding the antibacterial mechanism of IRFs in teleosts.

Transcriptome analysis of the spleen provides insight into the immunoregulation of Scortum barcoo under Streptococcus agalactiae infection

Jade perch (Scortum barcoo) is a freshwater fish with substantial economic value, which has been widely cultivated all over the world. However, with the intensification and expansion of farming, several bacterial and viral diseases have occurred in jade perch. To understand the immune response of jade perch against Streptococcus agalactiae (Group B Streptococcus, GBS), we performed a histopathological examination and transcriptome sequencing of jade perch spleen after artificial bacterial infection. GBS infection can cause structural changes and even necrosis of the jade perch spleen, which may affect the survival of infected individuals. A total of 144,458 unigenes were obtained through de novo assembly of spleen transcriptome. Among them, 1821 unigenes were identified as DEGs, including 1415 up-regulated and 406 down-regulated unigenes in the infection group. Moreover, the analysis of GO and KEGG revealed that many GO terms and pathways were involved in the host immune response, such as Toll-like receptor signaling pathway, Cytokine-cytokine receptor interaction, and TNF signaling pathway. In addition, according to transcriptome data and qRT-PCR analysis, the expression levels of many cytokines that participate in the inflammatory response changed a lot after GBS infection. Overall, this transcriptomic analysis provided valuable information for studying the immune response of jade perch against bacterial infection.

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.

Acute septicemia and immune response of spotted sea bass (Lateolabrax maculatus) to Aeromonas veronii infection

A previous study confirmed that spotted sea bass (Lateolabrax maculatus), an economically important cultured species in East Asia, is a new host of Aeromonas veronii, which can cause acute death in these fish, but there is little in-depth understanding of this disease. In the present study, the virulence of 10 isolates of A. veronii derived from spotted sea bass was determined. It was found that the 18BJ181 isolate was a virulent strain and led to the fastest death of spotted sea bass. Death was determined to be within in 2-12 h, and resulted in abdominal effusion and varying degrees of hemorrhage in internal organs. Bacterial colonization analysis showed that the bacterial load in the spleen was highest, and was up to 3.1 × 10⁵ cfu g^-1. In addition, the bacteria proliferated massively in the blood and reached 2.4 × 10⁷ cfu mL^-1 at 12 h after 18BJ181 strain infection, which was also a typical feature of acute septicemia. Histopathology of the spleen revealed edema in interstitial tissue, degeneration, and necrosis in lymphoid tissue, and hemorrhage in the capillary network. Transcriptome analysis of the spleen showed that the expression level of HSP70, CCL19, and IL-1β was extremely significantly up-regulated at 8 h after infection (P < 0.01), and the expression of these genes was normal at 24 h. These results revealed that A. veronii infection could rapidly activate the chemokine signal pathway and stimulate the acute inflammatory response in the host. The bacterial colonization, pathological features, and gene expression patterns in immune pathways will help us to better understand acute septicemia in spotted sea bass caused by A. veronii.

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.

Comparative genomics analysis of strains from diverse sources reveals the evolutionary relationship of Aeromonas veronii

Aeromonas veronii (A. veronii, AV) strains are emerging zoonotic and aquatic pathogens, yet we know very little about their genomics. This study aims to utilize comparative genomics to investigate the intraspecific genetic diversity, differences in virulence factors and evolutionary mechanisms of A. veronii strains from diverse sources and to fundamentally demonstrate their pathogenic mechanisms. We conducted comparative genomics analysis of 39 A. veronii strains from different sources and found that 1993 core genes are shared by these strains and that these shared core genes may be necessary to maintain the basic characteristics of A. veronii. Additionally, phylogenetic relationship analysis based on these shared genes revealed that a distant relationship between the AMC34 strain and the other 38 strains but that, the genetic relationship among the 38 strains is relatively close, indicating that AMC34 may not belong to A. veronii. Furthermore, analysis of shared core genes and average nucleotide identity (ANI) values showed no obvious correlation with the location of A. veronii isolation and genetic relationship. Our research indicates the evolutionary mechanism of A. veronii from different sources and provides new insights for a deeper understanding of its pathogenic mechanism.

Map of dimorphic switching‑related signaling pathways in Sporothrix schenckii based on its transcriptome

Sporothrix schenckii (S. schenckii) induces sporotrichosis, which has gained attention in recent years due to its worldwide prevalence. The dimorphic switching process is essential for the pathogenesis of S. schenckii. Previously, overexpression of several signal transduction genes, including SsDRK1 and SsSte20, was observed during the mycelium‑to‑yeast transition; these were necessary for asexual development, yeast‑phase cell formation, cell wall integrity and melanin synthesis. However, the mechanisms of the signaling pathways during dimorphic switching of S. schenckii remain unclear. In the present study, transcriptome sequencing of the 48‑h induced yeast forms and mycelium of S. schenckii was performed. In total, 24,904,510 high‑quality clean reads were obtained from mycelium samples and 22,814,406 from 48‑h induced yeast form samples. Following assembly, 31,779 unigene sequences were obtained with 52.98% GC content (The proportion of guanine G and cytosine C to all bases in nucleic acid). The results demonstrated that 12,217 genes, including genes involved in signal transduction and chitin synthesis, were expressed differentially between the two stages. According to these results, a map of the signaling pathways, including two‑component and heterotrimeric G‑protein signaling systems, Ras and MAPK cascades associated with the dimorphic switch, was drawn. Taken together, the transcriptome data and analysis performed in the present study lay the foundation for further research into the molecular mechanisms controlling the dimorphic switch of S. schenckii and support the development of anti‑S. schenckii strategies targeting genes associated with signaling pathways.

Novel insights into the immune regulatory effects of Megalobrama amblycephala intelectin on the phagocytosis and killing activity of macrophages

Previous studies have found that the expression level of Megalobrama amblycephala intelectin (MaINTL) increased significantly post Aeromonas hydrophila infection, and recombinant MaINTL (rMaINTL) protein could activate macrophages and enhance the phagocytosis and killing activity of macrophages. In order to reveal the immune regulatory mechanisms of MaINTL, primary M. amblycephala macrophages were treated with endotoxin-removed rMaINTL and GST-tag proteins, then total RNA were extracted and used for comparative Digital Gene Expression Profiling (DGE). 1247 differentially expressed genes were identified by comparing rMaINTL and GST-tag treated macrophage groups, including 482 up-regulated unigenes and 765 down-regulated unigenes. In addition, eleven randomly selected differentially expressed genes were verified by qRT-PCR, and most of them shared the similar expression patterns as that of DGE results. GO enrichment revealed that the differentially expressed genes were mainly concentrated in the membrane part and cytoskeleton of cellular component, the binding and signal transducer activity of molecular function, the cellular process, regulation of biological process, signaling and localization of biological process, most of which might related with the phagocytosis and killing activity of macrophages. KEGG analysis revealed the activation and involvement of differentially expressed genes in immune related pathways, such as Tumor necrosis factor (TNF) signaling pathway, Interleukin 17 (IL-17) signaling pathway, Toll-like receptor signaling pathway, and NOD like receptor signaling pathway, etc. In these pathways, TNF-ɑ, Activator protein-1 (AP-1), Myeloid differentiation primary response protein MyD88 (MyD88), NF-kappa-B inhibitor alpha (ikBɑ) and other key signaling factors were significantly up-regulated. These results will be helpful to clarify the immune regulatory mechanisms of fish intelectin on macrophages, thus providing a theoretical basis for the prevention and control of fish bacterial diseases.

Transcriptome analysis of the Macrobrachium nipponense hepatopancreas provides insights into immunoregulation under Aeromonas veronii infection

The oriental river prawn Macrobrachium nipponense is a commercially important freshwater shrimp that is widely farmed in China. Aeromonas veronii is a conditional pathogen of farmed shrimp, which has caused huge economic losses to the industry. Therefore, there is urgency to study the host-pathogen interactions between M. nipponense and A. veronii to screen individuals with antimicrobial resistance. In this study, we examined the hepatopancreas of moribund M. nipponense infected with A. veronii and healthy individuals at both the histopathological and transcriptomic levels. We showed that A. veronii infection resulted in tubular necrosis of the M. nipponense hepatopancreas. Such changes likely affect assimilation, storage, and excretion by the hepatopancreas, which could ultimately affect the survival and growth of infected individuals. Among the 61,345 unigenes obtained through RNA sequencing and de novo transcriptome assembly, 232 were differentially expressed between the two groups. KEGG and GO analyses revealed that these differentially expressed genes were implicated in pathways, including PPAR, PI3K/AKT, and AMPK signaling. The results of this study will contribute to an analysis of the immune response of M. nipponense to A. veronii infection at the transcriptomic level. Furthermore, the RNA-seq data generated here provide an important genomic resource for research on M. nipponense in the absence of a reference genome.

Transcriptome analysis reveals seven key immune pathways of Japanese flounder (Paralichthys olivaceus) involved in megalocytivirus infection

Megalocytivirus is a serious viral pathogen to many farmed fish including Japanese flounder (Paralichthys olivaceus). In this study, in order to systematically identify host immune genes induced by megalocytivirus infection, we examined the transcription profiles of flounder infected by megalocytivirus for 2, 6, and 8 days. Compared with uninfected fish, virus-infected fish exhibited 1242 differentially expressed genes (DEGs), with 225, 275, and 877 DEGs occurring at 2, 6, and 8 days post infection, respectively. Of these DEGs, 728 were upregulated and 659 were downregulated. The majority of DEGs were time-specific and formed four distinct expression profiles well correlated with the time of infection. The DEGs were classified into diverse Gene Ontology (GO) functional terms and enriched in 27 Kyoto Encyclopedia of Genes and Genomes (KEGG) pathways, approximately one third of which were related to immunity. Weighted co-expression network analysis (WGCNA) was used to identify 16 key immune DEGs belonging to seven immune pathways (RIG-I-like receptor signaling pathway, JAK-STAT signaling pathway, TLR signaling pathway, cytokine-cytokine receptor interaction, phagosome, apoptosis, and p53 signaling pathway). These pathways interacted extensively and formed complicated networks. This study provided a global picture of megalocytivirus-induced gene expression profiles of flounder at the transcriptome level and uncovered a set of key immune genes and pathways closely linked to megalocytivirus infection. These results provided a set of targets for future delineation of the key factors implicated in the anti-megalocytivirus immunity of flounder.

Transcriptome analysis and histopathological observations of Geloina erosa gills upon Cr(VI) exposure

The heavy metal contamination like Cr(VI) has been increased by human activities and that threats the ecosystem health of mangrove areas. Bioindicator is an emerging tool in the environmental contamination assessment. The objective of this study was to investigate the Geloina erosa response mechanisms and sensitivities of several biomarkers in the Cr(VI) exposure and identify the G. erosa capability of being used as heavy metals bioindicator. In this study, G. erosa was exposed to 100 μmol·L^-1 Cr(VI) for 48 h. After transcriptome sequencing, a total of 134,817 unigenes were obtained, including 12,555 up-regulated and 18,829 down-regulated differentially expressed genes and were validated through quantitative real-time PCR. In addition, a total of 12,185 SSRs and 1,428,214 candidate SNPs were identified from all the G. erosa transcriptome libraries. Histopathology of the gill indicated the Cr(VI) exposure induced damage of the organ leading to its immunization, detoxification or apoptosis reactions. Among eight genes of the selected biomarkers, Calm, HSP70, CYP450, ATG5, TLR2, MYD88 and CASP8 were up-regulated, while TLR4 was down-regulated in response to the Cr(VI) exposure.