Transcriptome profiling provides gene resources for understanding gill immune responses in Japanese flounder (Paralichthys olivaceus) challenged with Edwardsiella tarda

Whole genome identification of olive flounder (Paralichthys olivaceus) cathepsin genes: Provides insights into its regulation on biotic and abiotic stresses response

Cathepsins are major lysosomal enzymes involved in essential physiological processes, including protein degradation, tissue differentiation, and innate or adaptive responses. Several kinds of cathepsins have been reported in teleost fishes, but no characterization have been performed for the inflammatory response of cathepsin family in olive flounder until now. In our current study, a total of 17 cathepsins in olive flounder were systematically identified and characterized. Phylogenetic analysis clearly indicated that the cathepsin genes was highly conserved. Analysis of structure and motifs exhibited high sequence similarity of cathepsin genes in olive flounder. Expression profiles of cathepsin genes in different tissues and developmental stages showed that cathepsins were temporally and spatially specific. RNA-seq analysis of bacteria and temperature stresses revealed that members of cathepsin were involved in inflammatory responses. Collectively, our findings would provide a further reference for understanding the molecular mechanisms of cathepsins in olive flounder.

Development of Nanopore sequencing-based full-length transcriptome database toward functional genome annotation of the Pacific oyster, Crassostrea gigas

The Pacific oyster (Crassostrea gigas) is a widely cultivated shellfish in the world, while its transcriptome diversity remains less unexplored due to the limitation of short reads. In this study, we used Oxford Nanopore sequencing to develop the full-length transcriptome database of C. gigas. We identified 77,920 full-length transcripts from 21,523 genes, and uncovered 9668 alternative splicing events and 87,468 alternative polyadenylation sites. Notably, a total of 16,721 novel transcripts were annotated in this work. Furthermore, integrative analysis of 25 publicly available RNA-seq datasets revealed the transcriptome diversity involved in post-transcriptional regulation in C. gigas. We further developed a Drupal based webserver, Cgtdb, which can be used for transcriptome visualization, sequence alignment, and functional genome annotation analyses. This work provides valuable resources and a useful tool for integrative analysis of various transcriptome datasets in C. gigas, which will serve as an essential reference for functional annotation of the oyster genome.

Investigation of Gene Networks in Three Components of Immune System Provides Novel Insights into Immune Response Mechanisms against Edwardsiella tarda Infection in Paralichthys olivaceus

As a quintessential marine teleost, Paralichthys olivaceus demonstrates vulnerability to a range of pathogens. Long-term infection with Edwardsiella tarda significantly inhibits fish growth and even induces death. Gills, blood, and kidneys, pivotal components of the immune system in teleosts, elicit vital regulatory roles in immune response processes including immune cell differentiation, diseased cell clearance, and other immunity-related mechanisms. This study entailed infecting P. olivaceus with E. tarda for 48 h and examining transcriptome data from the three components at 0, 8, and 48 h post-infection employing weighted gene co-expression network analysis (WGCNA) and protein-protein interaction (PPI) network analysis. Network analyses revealed a series of immune response processes after infection and identified multiple key modules and key, core, and hub genes including xpo1, src, tlr13, stat1, and mefv. By innovatively amalgamating WGCNA and PPI network methodologies, our investigation facilitated an in-depth examination of immune response mechanisms within three significant P. olivaceus components post-E. tarda infection. Our results provided valuable genetic resources for understanding immunity in P. olivaceus immune-related components and assisted us in further exploring the molecular mechanisms of E. tarda infection in teleosts.

Effects of heat stress on physiological parameters, biochemical parameters and expression of heat stress protein gene in Lateolabraxmaculatus

With global warming and the increasing frequency of extreme heat, we have to consider the heat tolerance of fish to sudden high temperatures. This study investigated the effects of high temperature (32 °C) on physiology and biochemistry and heat shock proteins (HSPs) gene of spotted sea bass (Lateolabrax maculatus). The spotted sea bass (14.7 ± 1.54 g) temporarily cultured at 26 °C was directly transferred to the 32 °C high temperature group, and the gill morphology, liver antioxidant activity, respiratory metabolism related enzyme activity and the expression of five HSP70 family members genes were measured at 3, 6, 9, 12, 24, 48, 72, and 96 h. The results showed that 32 °C had damage effect on gill tissue and antioxidant system, and the damage degree increased with high temperature. Respiratory rate and malondialdehyde increased gradually with the continuous heat stress. Superoxide dismutase and total antioxidant capacity increased briefly and then decreased continuously. Succinate dehydrogenase decreased to the lowest value at 24 h and then continued to increase. Lactate dehydrogenase decreased continuously; the expression of HSP70 increased rapidly and then decreased. These results indicated that the antioxidant system and HSP70 were activated under heat stress condition and provided protection to the body, but with the continuous high temperature, the protective effect was limited, and the fish body was irreversibly damaged. It is necessary to pay close attention to the temperature change in the production practice of spotted sea bass to reduce the influence caused by high temperature.

Identification of the matrix metalloproteinase (MMP) gene family in Japanese flounder (Paralichthys olivaceus): Involved in immune response regulation to temperature stress and Edwardsiella tarda infection

The Matrix metalloproteinase (MMP) gene family is responsible for regulating the degradation of Extra Cellular Matrix (ECM) proteins, which are important for physiological processes such as wound healing, tissue remodeling, and stress response. Although MMPs have been studied in many species, their role in immune response in Japanese flounder (Paralichthys olivaceus) is still not fully understood. This study conducted a comprehensive analysis of MMPs in flounder, including gene structures, evolutionary relationships, conserved domains, molecular evolution, and expression patterns. Analysis revealed that MMP genes could be grouped into 17 subfamilies and were evolutionarily conserved and functionally-constrained. Meanwhile, MMP genes were found to express in different embryonic and larval stages and might play the role of sentinel in healthy tissues. Furthermore, expression profiling showed that MMPs had diverse functions in environmental stress, with 60% (9/15) and 73% (11/15) of MMPs showing differential expression patterns under temperature stress and Edwardsiella tarda (E. tarda) infection, respectively. These findings provide a useful resource for understanding the immune functions of MMP genes in Japanese flounder.

Efficacy of bivalent vaccine against Aeromonas salmonicida and Edwardsiella tarda infections in turbot

In recent years, more than one pathogenic organism has usually been isolated from diseased turbot Scophthalmus maximus, creating a pressing need for the development of combination vaccines to prevent fish diseases brought on simultaneously by various infections. In this study, the inactivated bivalent vaccine of Aeromonas salmonicida and Edwardsiella tarda was prepared by the formalin inactivation method. After challenge with A. salmonicida and E. tarda at 4 weeks post-vaccination in turbot, the relative percentage survival (RPS) of the inactivated bivalent vaccine was 77.1%. In addition, we assessed the effects of the inactivated bivalent vaccine and evaluated the immunological processes after immunization in a turbot model. Serum antibody titer and lysozyme activity of the vaccinated group were both upregulated and higher than that in control group after vaccination. The expression levels of genes (TLR2, IL-1β, CD4, MHCI, MHCⅡ) that related to antigen recognition, processing and presentation were also studied in the liver, spleen and kidney tissues of vaccinated turbot. All the detected genes in the vaccinated group had a significant upward trend, and most of them reached the maximum value at 3-4 weeks, which had significant differences from the control group, suggesting that antigen recognition, processing and presentation pathway was activated by the inactivated bivalent vaccine. Our study provides a basis for further application of the killed bivalent vaccine against A. salmonicida and E. tarda in turbot, making it good potential that can be applied in aquaculture.

Cecropin P1 antimicrobial peptide modulates differential expression of immune relevant genes in rainbow trout (Oncorhynchus mykiss) gill cell line, RTgill-W1

Accumulated evidence indicates that antimicrobial peptides modulate immune activities in fish. In this study, we profiled the differential expression patterns of representative immune relevant genes in an epithelial-like cell line of rainbow trout gill, RTgill-W1, in response to exposure of cecropin P1 antimicrobial peptide. RTgill-W1 cells were treated with synthetic cecropin P1 over time (0, 2, 4 and 24 h) with or without the present of lipopolysaccharide (LPS) or polyinosinic polycytidylic acid (PolyI:C). The relative abundances of each mRNA were measured by real-time quantitative PCR. The dose-response study revealed significant perturbations of mRNA levels of genes related to pro-inflammation, acute phase, surface proteins and transcription factors at 30 μM of cecropin P1. In addition, cecropin P1 altered the differential expression patterns that were induced by LPS or PolyI:C, at different time points in RTgill-W1. Overall, our results indicate that cecropin P1 exhibits pro-inflammation activity, modulate cell-cell interaction and cytokine signal transduction in rainbow trout gill cell, and may suggest a potential application of this peptide as an immune adjuvant for disease control in aquaculture.

Genome-wide identification of TNF receptor-associated factors in Japanese flounder (Paralichthys olivaceus) and functional analysis of resistance to temperature and Edwardsiella tarda stress

Tumor necrosis factor receptor-associated factors (TRAFs), as the signaling mediators of the tumor necrosis factor (TNFR) superfamily, toll-like receptors (TLR) and interleukin-1 receptor (IL-1R) superfamily, can activate downstream signal transduction pathways and play an important role in the body's immune process. In this study, six TRAF genes, namely PoTRAF2a, PoTRAF2b, PoTRAF3, PoTRAF4, PoTRAF6 and PoTRAF7, were identified and annotated in Japanese flounder by using bioinformatics methods. Phylogenetic analysis confirmed that TRAF genes can be divided into seven groups. Analysis of motif composition and gene structure demonstrated that all PoTRAF members were evolutionarily conserved. The expression patterns of PoTRAF genes were then further investigated in six different developmental stages and eleven tissues of healthy fish, and it was found that there were spatial and tissue specificities among the members. To investigate the immune response of Japanese flounder to abiotic and biotic stresses, we further analyzed the expression profile of PoTRAFs after temperature stress and pathogen challenge. The result showed that PoTRAF3 and PoTRAF4 were observably differentially expressed under temperature stress, indicating that they were involved in the immune response after temperature stress. The expression of PoTRAF2a, PoTRAF2b and PoTRAF4 was significantly different after E. tarda infection, suggesting that they might have antibacterial effects. These results would help to clarify the molecular roles of PoTRAF genes in the regulation of immune and inflammatory responses in Japanese flounder.

Focal adhesion signaling pathway involved in skin immune response of tongue sole Cynoglossus semilaevis to Vibrio vulnificus infection

Focal adhesion (FA) plays a key role in cell adhesion, migration and antibacterial immune, but it remained unclear in fish. In this study, half-smooth tongue sole Cynoglossus semilaevis were infected with Vibrio vulnificus, and then immune-related protein in the skin, especially for FA signaling pathway were screened and identified by iTRAQ analysis. Results showed that the differentially expressed proteins (DEPs) in skin immune response (eg., ITGA6, FN, COCH, AMBP, COL6A1, COL6A3, COL6A6, LAMB1, LAMC1, FLMNA) were firstly found in FA signaling pathway. Furthermore, the validation analysis of FA-related genes were basically consistent with the iTRAQ data at 36 hpi (r = 0.678, p < 0.01), and their spatio-temporal expressions were confirmed by qPCR analysis. The molecular characterization of vinculin of C. semilaevis was described. This study will provide a new perspective for understanding the molecular mechanism of FA signaling pathway in the skin immune response in marine fish.

Comprehensive analysis of immune parameters, mRNA and miRNA profiles, and immune genes expression in the gill of rainbow trout infected with infectious hematopoietic necrosis virus (IHNV)

Rainbow trout (Oncorhynchus mykiss) is a species of cold-water fish with important economic values, widely cultivated worldwide. However, the outbreak of infectious hematopoietic necrosis virus (IHNV) caused the large-scale death of rainbow trout and seriously restricted the development of the trout farming industry. In this study, the changes of immune parameters in different periods (6-, 12-, 24-, 48-, 72-, 96-, 120-, and 144 h post-infection (hpi)), transcriptome profiles of 48 hpi (T48G) compared to control (C48G), and key immune-related genes expression patterns were measured in rainbow trout gill following IHNV challenge through biochemical methods, RNA sequencing (RNA-seq), and quantitative real-time polymerase chain reaction (qRT-PCR). The results showed that alkaline phosphatase (AKP), acid phosphatase (ACP), total superoxide dismutase (T-SOD), alanine aminotransferase (ALT), and aspartate aminotransferase (AST) activities, as well as lysozyme (LZM) and malonaldehyde (MDA) content decreased and then increased during infection, and remained at a high level after 48 hpi (P < 0.05), whereas catalase (CAT) activity showed a significant peak at 48 hpi (P < 0.05). The mRNA and miRNA analysis identified 4343 differentially expressed genes (DEGs) and 11 differentially expressed miRNAs (DEMs), and numerous immune-related DEGs involved in the Toll-like receptor signaling pathway, apoptosis, DNA replication, p53 signaling, RIG-I-like receptor signaling pathway, and NOD-like receptor signaling pathway and expression were significantly up-regulated in T48Gm group, including tlr3, tlr7, tlr8, traf3, ifih1, trim25, dhx58, ddh58, hsp90a.1, nlrc3, nlrc5, socs3, irf3, irf7, casp7, mx1, and vig2. The integrated analysis identified several important miRNAs (ola-miR-27d-3p_R+5, gmo-miR-124-3-5p, ssa-miR-301a-5p_L+2, and ssa-miR-146d-3p) that targeted key immune-related DEGs. Expression analysis showed that tlr3, tlr7, traf3, ifih1, dhx58, hap90a.1, irf3, irf7, and mx1 genes increased and then decreased during infection, and peaked at 72 hpi (P < 0.05). However, trim25 expression peaked at 96 hpi (P < 0.05). This study contributes to understanding immune response of rainbow trout against IHNV infection, and provides new insights into the immune regulation mechanisms and disease resistance breeding studies.

Genome-wide identification of olive flounder (Paralichthys olivaceus) SOCS genes: Involvement in immune response regulation to temperature stress and Edwardsiella tarda infection

The suppressors of cytokine signaling (SOCS) gene family participates in development and immunity through negative regulation of cytokine signaling pathways. Although the immune response of SOCS gene family members has been extensively characterized in teleost, no similar study has been reported in olive flounder yet. In our present study, a total of 13 SOCSs in olive flounder were identified and characterized systematically. By querying the SOCS sequences of ten teleost fish species, we found there were exactly more members of SOCSs in fish than mammals, which indicated that there were more duplication events occurred in fish than in higher vertebrates. Phylogenetic analysis clearly illuminated that SOCS genes were highly conserved. The analysis of gene structure and motif showed SOCS proteins of olive flounder shared a high level of sequence similarity strikingly. The expression profiles of tissues and developmental stages indicated that SOCS members had a kind of specificity in temporality and spatiality. RNA-Seq analysis of temperature stress and E. Tarda infection demonstrated SOCS members were involved in inflammatory response. In a word, our results would provide a further reference for understanding the mechanism of SOCS genes in olive flounder.

Comprehensive transcriptomics and proteomics analysis of Carassius auratus gills in response to Aeromonas hydrophila

As one of the mucosal barriers, fish gills represent the first line of defense against pathogen infection. However, the exact mechanism of gill mucosal immune response to bacterial infection still needs further investigation in fish. Here, to investigate pathological changes and molecular mechanisms of the mucosal immune response in the gills of crucian carp (Carassius auratus) challenged by Aeromonas hydrophila, the transcriptomics and proteomics were performed by using multi-omics analyses of RNA-seq coupled with iTRAQ techniques. The results demonstrated gill immune response were mostly related to the activation of complement and coagulation cascades, antigen processing and presentation, phagosome, NOD-like receptor (NLR) and nuclear factor κB (NFκB) signaling pathway. Selected 21 immune-related DEGs (ie., Clam, nfyal, snrpf, acin1b, psme, sf3b5, rbm8a, rbm25, prpf18, g3bp2, snrpd3l, tecrem-2, cfl-A, C7, lysC, ddx5, hsp90, α-2M, C9, C3 and slc4a1a) were verified for their immune roles in the A. hydrophila infection via using qRT-PCR assay. Meanwhile, some complement (C3, C7, C9, CFD, DF and FH) and antigen presenting (HSP90, MHC Ⅱ, CALR, CANX and PSME) proteins were significantly participated in the process of defense against infections in gill tissues, and protein-protein interaction (PPI) network displayed the immune signaling pathways and interactions among these DEPs. The correlation analysis indicated that the iTRAQ and qRT-PCR results was significantly correlated (Pearson's correlation coefficient = 0.70, p < 0.01). To our knowledge, the transcriptomics and proteomics of gills firstly identified by multi-omics analyses contribute to understanding on the molecular mechanisms of local mucosal immunity in cyprinid species.

Effects of Zinc Adaptation on Histological Morphology, Antioxidant Responses, and Expression of Immune-Related Genes of Grass Carp (Ctenopharyngodon idella)

The study was conducted to evaluate the effect of zinc adaptation on histological morphology and antioxidant and immune responses of grass carp(Ctenopharyngodon idella). A total of 180 young grass carp (20.0 ± 2.0 g) was equally distributed into 9 groups, and triplicate groups were subjected to 0 μg/L Zn²⁺ (control group), 200 μg/L Zn²⁺, and 300 μg/L Zn²⁺ solution for 42 days, respectively. The results indicated that the liver and gill have obvious pathological changes after long-term adaptation to zinc except the intestine; the zinc adaptation can positively influence intestinal morphology. The activities of GPX (glutathione peroxidase activity), SOD (superoxide dismutase), and CAT (Catalase) were significantly increased in zinc treatment groups (P < 0.05). The genes expression levels of CuZnSOD (copper zinc superoxide dismutase), CAT, Hsp70 (heat shock protein-70), IL-1b (interleukin-1-b), and TGF-β1 (transforming growth factor-β1) were upregulated in the gill and intestine of grass carp following waterborne adaptation to zinc solution for 42 days (P < 0.05). In conclusion, zinc adaptation has different effects on organs of grass carp and may reduce the inflammatory response of the body's gills and intestines by improving the body's antioxidant and anti-stress defense capabilities.

Integrative analysis of miRNA and mRNA expression associated with the immune response in the intestine of rainbow trout (Oncorhynchus mykiss) infected with infectious hematopoietic necrosis virus

Rainbow trout (Oncorhynchus mykiss), an economically important cold-water fish cultured worldwide, suffers from infectious hematopoietic necrosis virus (IHNV) infection, resulting in huge financial losses. In order to understand the immune response of rainbow trout during virus infection, we explored trout intestine transcriptome profiles following IHNV challenge, and identified 3355 differentially expressed genes (DEGs) and 80 differentially expressed miRNAs (DEMs). Transcriptome analysis revealed numerous DEGs involved in immune responses, such as toll-like receptor 3 (TLR3), toll-like receptor 7/8 (TLR7/8), tripartite motif-containing 25 (TRIM25), DExH-Box helicase 58 (DHX58), interferon-induced with helicase C domain 1 (IFIH1), interferon regulatory factor 3 (IRF3/7), signal transducer and activator of transcription 1 (STAT1) and heat shock protein 90-alpha 1 (HSP90A1). Integrated analysis identified five key miRNAs (miR-19-y, miR-181-z, miR-203-y, miR-143-z and miR-206-y) targeting at least two important immune genes (TRIM25, DHX58, STAT1, TLR7/8 and HSP90A1). Gene Ontology and Kyoto Encyclopedia of Genes and Genomes analyses showed that DEGs and target genes were significantly enriched in various immune-related terms including immune system process, binding, cell part and pathways of Toll-like receptor signalling, RIG-I-like receptor signalling, NOD-like receptor signalling, JAK-STAT signalling, PI3K-Akt signalling, NF-kappa B signalling, IL-17 signalling and AGE-RAGE signalling. In addition, protein-protein interaction networks (PPI) was used to display highly interactive DEG networks involving eight immune-related pathways. The expression trends of 12 DEGs and 10 DEMs were further verified by quantitative real-time PCR, which confirmed the reliability of the transcriptome sequencing results. This study expands our understanding of the immune response of rainbow trout infected with IHNV, and provides valuable resources for future studies on the immune molecular mechanism and disease resistance breeding.

Splenic protection network revealed by transcriptome analysis in inactivated vaccine-immunized flounder (Paralichthys olivaceus) against Edwardsiella tarda infection

The protective immune response produced by fish after vaccination is crucial for vaccine effectiveness. Our previous studies have shown inactivated vaccine against Edwardsiella tarda can induce immune response in flounder (Paralichthys olivaceus). To elucidate the protective immune response at the genetic level, in this study, flounder was immunized with inactivated E. tarda for 5 weeks, and then they were challenged with E. tarda. The spleen was dissected at 7th day post immunization, 1st and 7th day post challenge, respectively. Transcriptome analysis showed that average of 46 million clean reads were obtained per library, while percentage of clean reads being mapped to reference genome was more than 89% in all cases, which suggested good quality of samples. As for differentially expressed genes (DEGs) identification in inactivated E. tarda groups, at 7th day post immunization, 1422 DEGs were identified and significantly enriched in innate immune-related pathways, such as Phagosome, Cell adhesion molecules and NF-kappa B signaling pathway; At 1st post challenge day, 1210 DEGs were identified and enriched to Antigen processing and presentation and Cell adhesion molecules, indicating that the pathogen was rapidly recognized and delivered; At 7th post challenge day, 1929 DEGs were identified, belonged to Toll-like receptor signaling pathway, Antigen processing and presentation, Th1 and Th2 cell differentiation and Th17 cell differentiation. Compared to 7th post immunization day, 73 immune-associated DEGs were identified at 1st post challenge day. Protein-protein interaction networks analysis revealed 11 hub genes (TLR7, TLR3, CXCR4, IFIH1, TLR8 etc), associated with recognition of pathogens and activation of innate immunity; while for 7th post challenge day, 141 immune-associated DEGs were identified. 30 hub genes (IL6, STAT1, HSP90A.1, TLR7, IL12β etc) were associated with stimulation of lymphocyte differentiation and activation of cellular immunity. Ten immune-related genes were randomly selected for RT-qPCR validation at each time point. In conclusion, data revealed protection of flounder against E. tarda infection by inactivated vaccine is mediated via immediate recognition of pathogen and subsequently activation of cellular immunity. Results give new aspect for vaccine protection cascades, is good references for vaccine evaluation.

Gene expression patterns in Atlantic salmon (Salmo salar) with severe nephrocalcinosis

Nephrocalcinosis is a common disorder in farmed Atlantic salmon, but the consequences for the fish physiology are not well understood. We performed a transcriptome study in kidneys of Atlantic salmon (Salmo salar) smolts without and with severe chronic nephrocalcinosis (NC). The study revealed that numerous genes are differentially expressed in fish with NC compared with healthy salmon. The most evident changes in gene expression patterns in the NC group were a massive downregulation of metabolism and energy production, upregulation of signalling pathways important for tissue repair and function maintenance and upregulation of inflammatory responses. Overall, the extensive tissue damage and the gene regulation responses that affect salmon with severe nephrocalcinosis are highly likely to have dramatic consequences on fish survival.

Identification and functional analysis of dual-specificity phosphatases (DUSP) genes in Japanese flounder (Paralichthys olivaceus) against temperature and Edwardsiella tarda stress

Dual-specificity Phosphatases (DUSPs) are not only the key regulators of dephosphorylating and inactivating mitogen-activated protein kinases (MAPKs), but play a crucial role in the immune response. However, the role of DUSP genes in Japanese flounder (PoDUSPs) is still unclear. In this study, 28 DUSP genes in Japanese flounder were identified and classified based on the whole genome database. Phylogenetic analysis and protein structure analysis revealed that DUSPs had highly conserved domains in teleosts. Molecular evolution analysis indicated that the PoDUSP genes were conservative during evolution and were functional-constrained. Meanwhile, PoDUSP genes were found to express in different embryonic and larval stages which might play the role of sentinel in healthy organisms. Furthermore, PoDUSP genes' expression profiles after temperature stress and Edwardsiella tarda (E. tarda) infection were determined in Japanese flounder without precedent, and the results demonstrated that Podusp1, Podusp2 and Podusp16 were more respective to temperature variation whereas Podusp1 and Podusp6 were more respective to E. tarda infection. In summary, our results provide useful resources for understanding the immune responsibilities of DUSP genes in flatfish.

Temperature induces metabolic reprogramming in fish during bacterial infection

Water temperature elevation as a consequence of global warming results in increased incidence of bacterial disease, such as edwardsiellosis, in fish farming. Edwardsiellosis is caused by the bacterial pathogen Edwardsiella tarda and affects many farmed fish including flounder (Paralichthys olivaceus). Currently, the effect of temperature on the metabolic response of flounder to E. tarda infection is unclear. In this study, we found that compared to low temperature (15°C), high temperature (23°C) enhanced E. tarda dissemination in flounder tissues. To examine the impact of temperature on the metabolism of flounder induced by E. tarda, comparative metabolomics were performed, which identified a large number of metabolites responsive to E. tarda invasion and temperature alteration. During E. tarda infection, the metabolic profile induced by elevated temperature was mainly featured by extensively decreased amino acids and TCA intermediates such as succinate, a proven immune regulator. Further, 38 potential metabolite markers of temperature effect (MMTE) in association with bacterial infection were identified. When used as exogenous supplements, two of the MMTE, i.e., L-methionine and UDP-glucose, effectively upregulated the expression of pro-inflammatory cytokines and suppressed E. tarda infection in flounder leukocytes. Taken together, the results of this study indicate an important influence of temperature on the metabolism of flounder during bacterial infection, which eventually affects the survivability of the fish.

Genome-wide identification and analysis of scavenger receptors and their expression profiling in response to Edwardsiella tarda infection in Japanese flounder (Paralichthys olivaceus)

The scavenger receptors (SRs) gene family, as one of pattern recognition receptors, participates in the innate immune response in diverse lineages. However, the systematic identification, characteristics and functions of SRs family are lacking in teleost. Here, we identified all 19 SRs family members in Japanese flounder (Paralichthys olivaceus) based on the genome and transcriptome data. Phylogenetic and Ka/Ks analysis demonstrated that these SRs genes were divided into five classes and all exhibited pronounced purified selection pressures. Whole genome duplication event was found in colec12, scarb2, and lamp1. Gene structure, functional domain and motif distribution analyses indicated that SRs within the different subfamilies are severely conservative. SRs genes showed diverse expression patterns in the embryogenesis and unchanged tissues. The regulations of 14 SRs genes in blood, gill and kidney after E. tarda infection suggested their roles in innate immune response. Meanwhile, ten SRs genes were differentially expressed after E. tarda stimulation in macrophages in vitro. Then we proved that PoSCARA3 could suppress the activity of NF-κB and AP-1 in HEK 293T cells by dual-luciferase assays. In summary, this study provided valuable basis for further functional characterization and immune functions of SRs genes in P. olivaceus.

The Isolation and Full-Length Transcriptome Sequencing of a Novel Nidovirus and Response of Its Infection in Japanese Flounder (Paralichthys olivaceus)

A novel nidovirus, CSBV Bces-Po19, was isolated from the marine fish, Japanese flounder (Paralichthys olivaceus). The viral genome was 26,597 nucleotides long and shared 98.62% nucleotide identity with CSBV WHQSR4345. PacBio Sequel and Illumina sequencing were used to perform full-length transcriptome sequencing on CSBV Bces-Po19-sensitive (S) and -resistant (R) Japanese flounder. The results of negative staining revealed bacilliform and spherical virions. There were in total 1444 different genes between CSBV Bces-Po19 S and R groups, with 935 being up-regulated and 513 being down-regulated. Metabolism-, immune-, and RNA-related pathways were significantly enriched. Furthermore, CSBV Bces-Po19 infection induced alternative splicing (AS) events in Japanese flounder; the S group had a higher numbers of AS events (12,352) than the R group (11,452). The number of long non-coding RNA (lncRNA) in the S group, on the other hand, was significantly lower than in the R group. In addition to providing valuable information that sheds more light on CSBV Bces-Po19 infection, these research findings provide further clues for CSBV Bces-Po19 prevention and treatment.

Transcriptome profiling based on larvae at different time points after hatching provides a core set of gene resource for understanding the immune response mechanisms of the egg-protecting behavior against Vibrio anguillarum infection in Amphioctopus fangsiao

Mollusks have recently received increasing attention because of their unique immune systems. Mollusks such as Amphioctopus fangsiao are economically important cephalopods, and the effects of their egg-protecting behavior on the larval immune response are unclear. Meanwhile, little research has been done on the resistance response of cephalopod larvae infected with pathogenic bacteria such as Vibrio anguillarum. In this study, V. anguillarum was used to infect the primary hatching A. fangsiao larvae under different egg-protecting behaviors for 24 h, and a total of 7156 differentially expressed genes (DEGs) were identified at four time points after hatching based on transcriptome analysis. GO and KEGG enrichment analyses showed that multiple immune-related GO terms and KEGG signaling pathways were enriched. Protein-protein interaction networks (PPI networks) were used to search functional relationships between immune-related DEGs. Finally, 20 hub genes related to multiple gene functions or involved in multiple signaling pathways were identified, and their accuracy was verified using quantitative RT-PCR. PPI networks were first used to study the effects A. fangsiao larvae after infection with V. anguillarum under different egg-protecting behaviors. The results provide significant genetic resources for exploring invertebrate larval immune processes. The data lays a foundation for further study the immune response mechanisms for invertebrates after infection.

Interacting Effects of Sea Louse (Lepeophtheirus salmonis) Infection and Formalin-Killed Aeromonas salmonicida on Atlantic Salmon Skin Transcriptome

Lepeophtheirus salmonis (sea lice) and bacterial co-infection threatens wild and farmed Atlantic salmon performance and welfare. In the present study, pre-adult L. salmonis-infected and non-infected salmon were intraperitoneally injected with either formalin-killed Aeromonas salmonicida bacterin (ASAL) or phosphate-buffered saline (PBS). Dorsal skin samples from each injection/infection group (PBS/no lice, PBS/lice, ASAL/no lice, and ASAL/lice) were collected at 24 h post-injection and used for transcriptome profiling using a 44K salmonid microarray platform. Microarray results showed no clear inflammation gene expression signatures and revealed extensive gene repression effects by pre-adult lice (2,189 down and 345 up-regulated probes) in the PBS-injected salmon (PBS/lice vs. PBS/no lice), which involved basic cellular (e.g., RNA and protein metabolism) processes. Lice repressive effects were not observed within the group of ASAL-injected salmon (ASAL/lice vs. ASAL/no lice); on the contrary, the observed skin transcriptome changes –albeit of lesser magnitude (82 up and 1 down-regulated probes)– suggested the activation in key immune and wound healing processes (e.g., neutrophil degranulation, keratinocyte differentiation). The molecular skin response to ASAL was more intense in the lice-infected (ASAL/lice vs. PBS/lice; 272 up and 11 down-regulated probes) than in the non-infected fish (ASAL/no lice vs. PBS/no lice; 27 up-regulated probes). Regardless of lice infection, the skin’s response to ASAL was characterized by the putative activation of both antibacterial and wound healing pathways. The transcriptomic changes prompted by ASAL+lice co-stimulation (ASAL/lice vs. PBS/no lice; 1878 up and 3120 down-regulated probes) confirmed partial mitigation of lice repressive effects on fundamental cellular processes and the activation of pathways involved in innate (e.g., neutrophil degranulation) and adaptive immunity (e.g., antibody formation), as well as endothelial cell migration. The qPCR analyses evidenced immune-relevant genes co-stimulated by ASAL and lice in an additive (e.g., mbl2b, bcl6) and synergistic (e.g., hampa, il4r) manner. These results provided insight on the physiological response of the skin of L. salmonis-infected salmon 24 h after ASAL stimulation, which revealed immunostimulatory properties by the bacterin with potential applications in anti-lice treatments for aquaculture. As a simulated co-infection model, the present study also serves as a source of candidate gene biomarkers for sea lice and bacterial co-infection.

Assessing the Effect of Modified Clay on the Toxicity of Karenia mikimotoi Using Marine Medaka (Oryzias melastigma) as a Model Organism

Blooms of the toxic dinoflagellate Karenia mikimotoi could threaten the survival of marine life, and modified clay (MC) is considered a promising method for the control of harmful algal blooms. Here, using marine medaka as the model organism, the toxicity of K. mikimotoi before and after MC disposal was investigated. The results showed that only a certain density of intact K. mikimotoi cells could cause obvious damage to fish gills and lead to rapid death. A systematic analysis of morphology, physiology, and molecular biology parameters revealed that the fish gills exhibited structural damage, oxidative damage, osmotic regulation impairment, immune response activation, and signal transduction enhancement. MC can flocculate K. mikimotoi rapidly in water and reduce its toxicity by reducing the density of intact algae cells and hemolytic toxicity. The results indicate that MC is an effective and safe method for controlling K. mikimotoi blooms.

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.

Identification, evolution and expression analyses of mapk gene family in Japanese flounder (Paralichthys olivaceus) provide insight into its divergent functions on biotic and abiotic stresses response

Mitogen-activated protein kinase (MAPK) are a series of serine/threonine protein kinases showing evolutionary conservation, which can be activated by many stimulus signals and then transfer them from cell membrane to nucleus. MAPKs regulate a variety of biological processes, such as apoptosis, hormone signaling and immune response. In this study, 14 putative mapk genes in Japanese flounder were identified, and their basic physical and chemical properties were characterized. Phylogenetic analysis showed that mapk genes were divided into three main subfamilies, including ERK, JNK and the p38 MAPK. Selection pressure analysis revealed they were evolutionarily-constrained and undergone strong purifying selection. Gene structure and conserved protein motif comparison suggested high levels of conservation in members of mapk gene family. The expression patterns were further investigated in each embryonic and larval development stages and different tissues. In addition, RNA-seq analyses after bacteria and temperature stresses suggested mapk genes had different expression patterns. Three mapk genes showed significant differences in response to E. tarda challenge and five were induced significantly after temperature stress, indicating their potential functions. This systematic analysis provided valuable information for further understanding of the regulation mechanism of mapk gene family under different stresses in Japanese flounder.

M. fortuitum-induced CNS-pathology: Deciphering the role of canonical Wnt signaling, blood brain barrier components and cytokines

Molecular underpinning of mycobacteria-induced CNS-pathology is not well understood. In the present study, zebrafish were infected with Mycobacterium fortuitum and the prognosis of CNS-pathogenesis studied. We observed M. fortuitum triggers extensive brain-pathology. Evans blue extravasation demonstrated compromised blood-brain barrier (BBB) integrity. Further, decreased expression in tight-junction (TJ) and adherens junction complex (AJC) genes were noted in infected brain. Wnt-signaling has emerged as a major player in host-mycobacterial immunity but its involvement/role in brain-infection is not well studied. Sustained expression of wnt2, wnt3a, fzd5, lrp5/6 and β-catenin, with concordant decline in degradation complex components axin, gsk3β and β-catenin regulator capn2a were observed. The surge in ifng1 and tnfa expression preceding il10 and il4 suggested cytokine-interplay critical in M. fortuitum-induced brain-pathology. Therefore, we suggest adult zebrafish as a viable model for studying CNS-pathology and using the same, conclude that M. fortuitum infection is associated with repressed TJ-AJC gene expression and compromised BBB permeability. Our results implicate Wnt/β-catenin pathway in M. fortuitum-induced CNS-pathology wherein Th1-type signals facilitate bacterial clearance and Th2-type signals prevent the disease sequel.

RGD-binding integrins mediated phagocytosis involved in the entry of Edwardsiella tarda into mudskipper MO/MФ

The versatile fish pathogen Edwardsiella tarda is an intracellular pathogen with the ability to invade and replicate in host phagocytes. However, the mechanism mediating the uptake of E. tarda in fish monocytes/macrophages (MO/MΦ) is not yet understood. Generating mudskipper kidney-derived MO/MФ transcriptomic resources from mudskipper challenged by E. tarda is crucial for understanding the molecular mechanisms underlying the mudskipper invasion process. In the present study, a total of 1185 up-regulated and 885 down-regulated differentially expressed genes (DEGs) were identified using RNA-seq. Enrichment and pathway analysis of DEGs revealed the centrality of the phagosome and regulation of actin cytoskeleton pathways in pathogen entry. The progress of phagosome formation was observed by transmission electron microscopy. Eight conserved integrin (ITG) subunit genes, belonging to the phagocytic receptors, were found in the transcriptomic sequence data. Additionally, quantitative real-time PCR showed that the mRNA expressions of most ITG subunit genes were related to the different infection times of E. tarda and the different bacterial pathogens. Further assays demonstrated that phagocytosis of FITC-labeled E. tarda by mudskipper MO/MФ was significantly reduced by the tetrapeptide Asp-Gly-Arg-Ser (RGDS). In summary, phagocytosis is one of the entry pathways into mudskipper MO/MΦ, and RGD-binding ITGs are involved in the phagosome formation process.

The Dynamic Immune Response of Yellow Catfish (Pelteobagrus fulvidraco) Infected With Edwardsiella ictaluri Presenting the Inflammation Process

Edwardsiella ictaluri is a highly destructive pathogen in cultured yellow catfish, thus it was very necessary to study the immune response of yellow catfish against bacterial infection. In this study, RNA-Seq technology was used to study the immune response in two distinct tissues of yellow catfish at eight different time points (h) after E. ictaluri infection. The number of differentially expressed genes (DEGs) in the spleen and liver was low at 3 h and 6 h post-infection, respectively. Afterwards, the most number of DEGs in the spleen was detected at 72 h, while the number of DEGs in the liver maintained a high level from 24 h to 120 h. The GO and KEGG enrichment analyses of DEGs at different time points uncovered that cytokines were continuously transcribed at 6 h to 120 h; whereas the liver is the main organ that secretes the components of the complement system, and metabolic regulation was activated from 12 h to 120 h. Moreover, an overview of the inflammation response of yellow catfish was exhibited including pattern-recognition receptors, inflammatory cytokines, chemokines, complements, and inflammation-related signal pathways. The similar expression tendency of nine genes by qRT-PCR validated the accuracy of transcriptome analyses. The different transcriptomic profiles obtained from the spleen and liver will help to better understand the dynamic immune response of fish against bacterial infection, and will provide basic information for establishing effective measures to prevent and control diseases in fish.

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)

Edwardsiella anguillarum is one of the common bacterial pathogens for the cultivated eels in China. The aim of this study was to reveal the cause of E. anguillarum pathogenic to European eel (Anguilla anguilla) from the perspective of the transcriptome. In this study, we first prepared E. anguillarum cultured in vitro and analysed the whole transcriptome after extracting the total RNA. Then, eels were i.p injected with E. anguillarum, and total RNA were extracted from the liver of European eels 48 h after the infection. After sequencing the transcriptome, we obtained average 1.97 × 10⁸ clean reads cultured in vitro and 1.36 × 10⁵ clean reads located in vivo after annotating all reads into the genome of E. anguillarum. The whole transcriptome showed, compared to the E. anguillarum cultured in vitro, 503 significantly up and 657 significantly down-regulated different expressed genes (DEGs) were observed. KEGG analysis showed that 38 DEGs of Two-Component System, 41 DEGs of ABC transporter, and 10 DEGs flagellar assembly pathways were highly upregulated in E. anguillarum located in vivo. Then, we designed primers to analyse the up-regulated DEGs through qRT-PCR and confirmed some up-regulated DEGs. The results of this study provide important reference for the further study of pathogen-host interaction between E. anguillarum and European eel.

The Effect of Acute Erythromycin Exposure on the Swimming Ability of Zebrafish (Danio rerio) and Medaka (Oryzias latipes)

Erythromycin is a widely used antibiotic, and erythromycin contamination may pose a threat to aquatic organisms. However, little is known about the adverse effects of erythromycin on swimming ability. To quantify erythromycin-induced damage to fish swimming ability, Oryzias latipes and Danio rerio were acutely exposed to erythromycin. The swimming ability of the experimental fish was measured after exposure to varying doses of erythromycin (2 µg/L, 20 µg/L, 200 µg/L, and 2 mg/L) for 96 h. Burst speed (U_burst) and critical swimming speed (U_crit) of experimental fish significantly decreased. In addition, gene expression analysis of O. latipes and D. rerio under erythromycin treatment (2 mg/L) showed that the expression of genes related to energy metabolism in the muscle was significantly reduced in both species of fish. However, the gene expression pattern in the head of the two species was differentially impacted; D. rerio showed endocrine disruption, while phototransduction was impacted in O. latipes. The results of our study may be used as a reference to control erythromycin pollution in natural rivers.

Gene network analysis reveals a core set of genes involved in the immune response of Japanese flounder (Paralichthys olivaceus) against Vibrio anguillarum infection

Japanese flounder (Paralichthys olivaceus) is one of the most economically important marine fish cultured in north Asia. Vibrio anguillarum is a severe bacterial pathogen to Japanese flounder and many other aquaculture species. In order to understand the immune response of flounder during bacterial infection, we systematically examined the transcriptome profiles of flounder spleen at three time points after V. anguillarum challenge. More than one billion high quality reads were obtained, approximately 80.70% of which were successfully mapped to the reference genome of flounder. A total of 6060, 4688 and 4235 differentially expressed genes (DEGs) were captured at 6, 12 and 24-h post-infection, respectively. The DEGs exhibited dynamic changes in expression and were assigned into four different profiles based on expression trend. GO and KEGG analysis showed that the DEGs were enriched in various immune-related terms, including response to stimulation, immune system and pathways of cytokine-cytokine receptor interaction, Jak-STAT signaling and Toll-like receptor signaling. Furthermore, a network of highly interactive DEGs involved in 11 immune-related pathways was detected by utilizing the weighted co-expressing network analysis (WGCNA). Accordingly, 26 hub genes were discovered that constituted an elaborate immune regulatory network and functioned mainly in pathogen recognition, antigen processing, and molecular signaling. The results of this study provided the first systematical transcriptome profile of flounder in association with V. anguillarum infection and can serve as a valuable resource of target genes for future studies on the molecular mechanisms underlying the immune defense of flounder against bacterial infection.

Identification of Multigene Biomarker for Shrimp White Feces Syndrome by Full-Length Transcriptome Sequencing

The pacific white shrimp, Litopenaeus vannamei, with the largest shrimp industry production in the world, is currently threatened by a severe disease, white feces syndrome (WFS), which cause devastating losses globally, while its causal agents remain largely unknown. Herein, compared to the Control shrimp by metagenomic analysis, we firstly investigated that the altered functions of intestinal microbial community in WFS shrimp were the enrichment of bacterial chemotaxis and flagellar assembly pathways, hinting at a potential role of pathogenic bacteria for growth and development, which might be related to WFS occurrence. Single-molecule real-time (SMRT) sequencing was to further identify the gene structure and gene regulation for more clues in WFS aetiology. Totally 50,049 high quality transcripts were obtained, capturing 39,995 previously mapped and 10,054 newly detected transcripts, which were annotated to 30,554 genes. A total of 158 differentially expressed genes (DEGs) were characterized in WFS shrimp. These DEGs were strongly associated with various immune related genes that regulated the expression of multiple antimicrobial peptides (e.g., antilipopolysaccharide factors, penaeidins, and crustin), which were further experimentally validated using quantitative PCR on transcript level. Collectively, multigene biomarkers were identified to be closely associated with WFS, especially those functional alterations in microbial community and the upregulated immune related gene with antibacterial activities. Our finding not only inspired our cogitation on WFS aetiology from both microbial and host immune response perspectives with combined metagenomic and full-length transcriptome sequencing, but also provided valuable information for enhancing shrimp aquaculture.

Comparative transcriptome analysis of the gills of Procambarus clarkii provides novel insights into the immune-related mechanism of copper stress tolerance

The red-swamp crayfish (Procambarus clarkii) is the most important economic shrimp species in China, and is an important model crustacean organism in many fields of research. In crustaceans, gills interface directly with the ambient environment and thus play a vital role in the toxicology. In the context of increasing environmental heavy metal pollution, the relationship between copper (Cu²⁺) stress and the immune response of P. clarkii has recently received considerable attention. However, impact of Cu²⁺ on the crayfish immune system is still not fully understood. In this study, we used Illumina sequencing technology to perform a transcriptome analysis of the gills of P. clarkii after 24 h of Cu²⁺ treatment. A total of 37,226,812 unigenes were assembled, and 1943 unigenes were significantly differentially expressed between the control and Cu²⁺ treatment groups. Functional categorization of differentially expressed genes (DEGs) revealed that genes related to antioxidant activity, detoxication, metabolic processes, biosynthetic processes, and immune system processes were differentially regulated during Cu²⁺ stress. In addition, DEGs in the immune system were classified as being related to the MAPK signaling pathway, purine metabolism, Toll and Imd signaling pathway, PI3K-Akt signaling pathway and Hippo signaling pathway. Five genes (CuZnSOD, CAT, IDH1, PHYH and DECR2) were significantly up-regulated in the peroxisome pathway, which plays an important role in reacting to oxidative stress. Importantly, qRT-PCR validation of the results for seven genes chosen at random (NDK, ATP6L, ATP5C1, RPS14, RPL22e, CTSF and HSP90A) confirmed the Illumina sequencing results. This study provides a valuable starting point for further studies to elucidate the molecular basis of the immune system's response to Cu²⁺ stress in crayfish.

Integrated transcriptome analysis of immunological responses in the pearl sac of the triangle sail mussel (Hyriopsis cumingii) after mantle implantation

For pearl culture of bivalve Hyriopsis cumingii, implantation of the sabio may cause nucleus discharge and increased host death rates. We performed a transcriptome analysis of the pearl sac of H. cumingii for 30 days after mantle implantation; 293863 unigenes were obtained, and 27176 unigenes were identified using nr, nt, KO, Swiss-Prot, Pfam, GO, and KOG databases. We detected 4878 differentially expressed genes (DEGs) through pairwise comparisons. We speculated that the physical condition of the recipient mussels returned to normal in about one month; the period was divided into six vital phases (0, 2 h-6 h, 12 h-24 h, 48 h to 7 days, 14 days and 30 days) on the basis of the overall similarities in DEGs. We compared the DEGs between time points and identified key immune-related genes. Our findings provide information on the immunological reactions induced by implantation in pearl mussels.

Edwardsiella tarda-induced miR-7a functions as a suppressor in PI3K/AKT/GSK3β signaling pathway by targeting insulin receptor substrate-2 (IRS2a and IRS2b) in Paralichthys olivaceus

To study the effect of Edwardsiella tarda infection on miRNAs expression profile in Japanese flounder, fish were injected intraperitoneally with E. tarda. The miRNAs involved in regulating immune responses were analyzed by high-throughput sequencing. A total of 164 mature miRNAs were identified, of which 17 miRNAs were differentially expressed (DE miRNAs) after E. tarda infection, indicating that they were immune-related miRNAs. To further examine the relationship between the miRNAs and their predicted target mRNAs, a total of 22 predicted target mRNAs, mainly related to endocytic signaling pathway, NF-κB signaling pathway, and p53 signaling pathway, were detected with miRNA mimics in HEK-293T cells by dual-luciferase reporter experiments. Finally, we confirmed that insulin receptor substrate-2 (IRS2a and IRS2b) were regulated by miR-7a. And the target sites of the 3' untranslated region (UTR) of IRS2a and IRS2b were verified by dual-luciferase reporter experiments. Furthermore, we found that the E. tarda and LPS significantly increased host miR-7a expression. In vivo and in vitro studies revealed that IRS2-mediated PI3K/AKT/GSK3β signaling pathway was suppressed. Taken together, these results implied that miR-7a might be a key regulator of PI3K/AKT/GSK3β signaling pathway via suppressing the IRS2a and IRS2b genes.

Impaired immune function and structural integrity in the gills of common carp (Cyprinus carpio L.) caused by chlorpyrifos exposure: Through oxidative stress and apoptosis

As one of the mucosal lymphatic tissues, the gill is an important immune organ in fish. Water environmental pollutants enter fish body through the gill. Therefore, the gill is the initial site where pollutants produce toxic effects in water. Chlorpyrifos (CPF), a broad-spectrum organophosphate insecticide, is widely used for agricultural pests and causes river pollution. In the present study, we investigated histopathological effect, oxidative stress indexes (SOD, GSH, T-AOC, and MDA), and apoptosis-related genes (P53, PUMA, Bax, Bcl-2, Apaf-1, Caspase-9, and Caspase-3) in the gills of common carp exposed to CPF. The results indicated that CPF exposure decreased SOD, T-AOC, and GSH; increased MDA; decreased Bcl-2 mRNA expression; and increased P53, PUMA, Bax, Apaf-1, Caspase-9, and Caspase-3 mRNA expressions in common carp gills. Our results proved that CPF exposure caused oxidative stress and apoptosis in common carp gills; CPF exposure destroyed the structural integrity and affected the immune function through oxidative stress and apoptosis in common carp gills. These will provide evidence for the toxic effects of water environmental pollutants on immune function and structural integrity in fish gills.

Comparative Study of Immune Reaction Against Bacterial Infection From Transcriptome Analysis

Transcriptome analysis is a powerful tool that enables a deep understanding of complicated physiological pathways, including immune responses. RNA sequencing (RNA-Seq)-based transcriptome analysis and various bioinformatics tools have also been used to study non-model animals, including aquaculture species for which reference genomes are not available. Rapid developments in these techniques have not only accelerated investigations into the process of pathogenic infection and defense strategies in fish, but also used to identify immunity-related genes in fish. These findings will contribute to fish immunotherapy for the prevention and treatment of bacterial infections through the design of more specific and effective immune stimulants, adjuvants, and vaccines. Until now, there has been little information regarding the universality and diversity of immune reactions against pathogenic infection in fish. Therefore, one of the aims of this paper is to introduce the RNA-Seq technique for examination of immune responses in pathogen-infected fish. This review also aims to highlight comparative studies of immune responses against bacteria, based on our previous findings in largemouth bass (Micropterus salmoides) against Nocardia seriolae, gray mullet (Mugil cephalus) against Lactococcus garvieae, orange-spotted grouper (Epinephelus coioides) against Vibrio harveyi, and koi carp (Cyprinus carpio) against Aeromonas sobria, using RNA-seq techniques. We demonstrated that only 39 differentially expressed genes (DEGs) were present in all species. However, the number of specific DEGs in each species was relatively higher than that of common DEGs; 493 DEGs in largemouth bass against N. seriolae, 819 DEGs in mullets against L. garvieae, 909 in groupers against V. harveyi, and 1471 in carps against A. sobria. The DEGs in different fish species were also representative of specific immune-related pathways. The results of this study will enhance our understanding of the immune responses of fish, and will aid in the development of effective vaccines, therapies, and disease-resistant strains.

Transcriptome profiling reveals key roles of phagosome and NOD-like receptor pathway in spotting diseased Strongylocentrotus intermedius

Spotting disease is a common disease in the process of aquaculture and restocking of the sea urchin Strongylocentrotus intermedius and leads to mass mortality. To characterize the molecular processes and candidate genes related to spotting disease in S. intermedius, we conducted next-generation sequencing to assess the key genes/pathways in spotting diseased sea urchin (DUG) compared to healthy ones (HUG). A total of 321.1 million clean reads were obtained and assembled into 93,877 Unigenes with an N50 of 1185 bp, in which 86.48% of them matched to the genome sequence of the sea urchin S. purpuratus and 27,456 Unigenes mapped to Nr database. Salmon expression analysis revealed 1557 significantly differently expressed genes (DEGs) between DUG and HUG. These DEGs were enriched into 151 KEGG pathways including a core set of immune correlated pathways notably in phagosome and NOD-like receptor signaling. DUG displayed an obvious downregulation in these immune pathways. The expression patterns of six DEGs were confirmed by RT-qPCR, and the expressions were consistent with the results of RNA-seq. Furthermore, 15,990 SSRs were identified and a total of 235,249 and 295,567 candidate SNPs were identified from DUG and HUG, respectively. All these results provided basic information for our understanding of spotting disease outbreak in sea urchin.

Transcriptome Analysis Based on RNA-Seq in Understanding Pathogenic Mechanisms of Diseases and the Immune System of Fish: A Comprehensive Review

In recent years, with the advent of next-generation sequencing along with the development of various bioinformatics tools, RNA sequencing (RNA-Seq)-based transcriptome analysis has become much more affordable in the field of biological research. This technique has even opened up avenues to explore the transcriptome of non-model organisms for which a reference genome is not available. This has made fish health researchers march towards this technology to understand pathogenic processes and immune reactions in fish during the event of infection. Recent studies using this technology have altered and updated the previous understanding of many diseases in fish. RNA-Seq has been employed in the understanding of fish pathogens like bacteria, virus, parasites, and oomycetes. Also, it has been helpful in unraveling the immune mechanisms in fish. Additionally, RNA-Seq technology has made its way for future works, such as genetic linkage mapping, quantitative trait analysis, disease-resistant strain or broodstock selection, and the development of effective vaccines and therapies. Until now, there are no reviews that comprehensively summarize the studies which made use of RNA-Seq to explore the mechanisms of infection of pathogens and the defense strategies of fish hosts. This review aims to summarize the contemporary understanding and findings with regard to infectious pathogens and the immune system of fish that have been achieved through RNA-Seq technology.