Introduction {#s1}
============

Transcriptome analysis is used to study principal pathways of development, cellular fate, physiology, activity, and disease progression. RNA sequencing (RNA-Seq) is a modern technology for transcriptome profiling that uses next-generation sequencing (NGS). Advancements in bioinformatics has significantly supported RNA-Seq technology to accelerate the knowledge on transcriptomes ([@B1]).

In the aquaculture field, there is wide utility for RNA-Seq in various applications, such as understanding the development of embryo and larvae, toxicology, environmental stress, effect of dietary conditions, and discovery of novel transcripts ([@B2]--[@B4]). In addition, RNA-Seq has been used in many studies of fish immunology ([@B5], [@B6]). Pathogenic infection is a major concern for maintaining economic sustainability in natural and farmed fish; it results in high mortality and economical loss in aquaculture. An innate immunity is a front line of host defense, producing effectors that directly to the pathogen and attack it. An adaptive immune system is also present in teleost, including humoral and cellular mechanisms. To reduce disease outbreaks, it is essential to understand the immune mechanisms in fish during pathogenic infections. This knowledge will support the development of effective vaccines and adjuvants against pathogens. However, it has not been reported that the universality and diversity of immune reactions against pathogenic infection in fish.

In this paper, we firstly introduce the RNA-Seq technique and current knowledge for investigations of immune responses in pathogen-infected fish. This review also aims to highlight comparative studies of fish immune responses against bacteria based on our previous studies that we demonstrated the transcriptome of bacteria infected fish, 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*.

Advantage of RNA-seq Analysis in Fish Aquaculture {#s2}
=================================================

Since the 2000s, hybridization-based microarray has been used to examine fish immunology in aquacultures; some early examples included Japanese flounder (*Paralichthys olivaceus*) ([@B7]), rainbow trout (*Oncorhynchus mykiss*) ([@B8]), and Atlantic salmon (*Salmo salar*) ([@B9]). Although species-specific probes should be designed, microarray technology could provide us with high throughput gene expression data. For transcriptome analysis, series analysis of gene expression (SAGE) and cap analysis gene expression (CAGE) have also been utilized. SAGE and CAGE, which are tag-based technology, are more precise; however, the number of genes that can be analyzed at one time is lower as compared with that of microarrays. Recently, the number of reports that uses RNA-Seq in aquaculture studies has rapidly increased. For the complete detail of RNA-Seq methodology, please refer the nicely reviews focused on aquaculture field ([@B3], [@B6]). The advantage of RNA-Seq is that it could determine expression levels of low-level transcripts as well as each splice variant isoforms. Current focus in aquaculture fish research is to examine organisms that do not process reference gene sequences; unigenes are obtained via *de novo* assembly using Trinity or similar programs without requiring reference gene sequences ([@B10]). RNA-Seq could provide novel transcript sequences, thereby expanding our current list of known transcripts in fish.

Annotation, Enrichment Analysis, and Pathway Analysis Using *de novo* Assembly Data {#s3}
===================================================================================

Although RNA-Seq technology could be applied to non-model animals, there are problems associated with the functional annotation and enrichment analysis of transcripts data. Transcripts sequences following assembly are usually searched via several databases, such as NCBI nucleotide sequences (NT), NCBI non-redundant protein (NR), Clusters of Orthologous Groups (COGs) ([@B11]), Kyoto Encyclopedia of Genes and Genomes (KEGG) ([@B12], [@B13]), gene ontology (GO) ([@B14]), and InterPro annotation ([@B15]). In our previous study on orange-spotted grouper (*Epinephelus coioides*), a total of 79,128 unigenes were identified and aligned with each database; 58,926 (74.47%) in NT, 43,576 (55.07%) in NR, 14,750 (18.64%) in COG, 34711 (43.87%) in KEGG, and 4232 (5.35%) in GO ([@B16]). The number of genes aligned with existing genes in the database was the highest in NT, while gene alignment was the lowest in GO. Differences in the number of aligned genes between the databases were similar to those found in other aquaculture studies ([@B17]--[@B19]). Although the GO database could provide enrichment analysis and pathway analysis, due to low gene alignment, results may be limited to generalized conclusions. There are programs that can convert gene IDs from one database to those of another database. For example, DAVID and ID Converter Systems are able to change gene IDs from NT to that of GO. However, these systems are not very useful in genes of aquaculture species. Currently, the KEGG database shows a relatively high number of aligned genes, which allows for enrichment analysis of aquaculture species. To obtain generalized conclusions using transcriptome data, it is essential that systems are developed to describe non-model organisms. There are databases now under construction that contain transcriptome information of aquaculture species. In the European common carp (*Cyprinus carpio*), a wide range of data on tissue-specific gene expression and translation ([@B20]) has been presented. These datasets will allow us to investigate immune responses in aquaculture species via transcriptome analysis.

Immune Responses Against Pathogens in Teleost Using RNA-seq {#s4}
===========================================================

Although fish are constantly exposed to risk of microorganism pathogens, fish could keep in shape to act immune mechanisms against pathogens. In a first line of immune response, fish are protected by non-specific humoral factors including growth inhibiting substances (transferrin and antiproteases), lysins (lysozyme, C-reactive protein, and bactericidal peptides), and making a link with non-specific phagocyte responses. Second, fish produce antibody constitutes for a specific humoral defense inhibiting bacterial adherence and invasion of non-phagocytic host cells and counteracting toxins from bacterial ([@B21]). After developing the technology of molecular biology, the immune-related gene functions and responses against pathogens have been one of the major topic in fish immunology field. To investigate the expression pattern of immune-associated genes, real-time PCR is usually performed. However, this method is expensive and not recommended for a genome-wide survey of gene expression. As RNA-Seq could provide us with quantitative data on transcript expression levels, this technique has been commonly used to identify genes that respond to pathogenic conditions during exogenous challenge. [Table 1](#T1){ref-type="table"} lists transcriptome analysis studies that examined immune regulations in teleost. In this review, we will also focus on studies that investigated immune responses to pathogen or their mimic molecules, using RNA-Seq analysis.

###### 

List of Transcriptome report in immune systems of teleost.

  ----------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------
  **Pathogen or stimulator**      **Fish**                        **Target tissue or cell**            **Responding genes and pathways**                                                                                                                                                                **References**
  ------------------------------- ------------------------------- ------------------------------------ ------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------ ----------------
  **GRAM-NEGATIVE BACTERIA**                                                                                                                                                                                                                                                                            

  *Aeromonas hydrophila*          Blunt snout bream\              Mixture of RNA sample form tissues   Toll-like receptor signaling, complement cascade.                                                                                                                                                ([@B22])
                                  (*Megalobrama amblycephala*)                                                                                                                                                                                                                                          

  *Aeromonas hydrophila*          Darkbarbel catfish\             Spleen                               Complement and coagulation cascades, pattern recognition pathways, Nk cell mediated cytotoxicity, Fcγ receptor-mediated phagocytosis, B cell and T cell receptor signaling pathway.              ([@B23])
                                  (*Pelteobagrus vachellii*)                                                                                                                                                                                                                                            

  *Aeromonas hydrophila*          Golden mahseer\                 Liver                                Th1 and 2 cell differentiation, pathogen recognition, complement cascade pathways.                                                                                                               ([@B24])
                                  (*Tor putitora*)                                                                                                                                                                                                                                                      

  *Aeromonas hydrophila*          Grass carp\                     Spleen                               Antigen processing and presentation, Phagocytosis, Complement system, Cytokines                                                                                                                  ([@B25])
                                  (*Ctenopharygodon idella*)                                                                                                                                                                                                                                            

  *Aeromonas hydrophila*          Grass carp\                     Spleen                               Complement and coagulation cascades pathway.                                                                                                                                                     ([@B26])
                                  (*Ctenopharyngodon idella*)                                                                                                                                                                                                                                           

  *Aeromonas salmonicida*         Rainbow trout\                  Spleen                               Hematopoietic cell lineage and complement and coagulation cascades.                                                                                                                              ([@B27])
                                  (*Oncorhynchus mykiss*)                                                                                                                                                                                                                                               

  *Aeromonas sobria*              Koi Carp\                       Muscle, spleen                       Toll-like receptor signaling and complement cascade pathway.                                                                                                                                     ([@B19])
                                  (*Cyprinus carpio*)                                                                                                                                                                                                                                                   

  *Edwardsiella ictaluri*         Yellow catfish\                 Spleen                               DEGs linked to several innate and adaptive immune related pathways.                                                                                                                              ([@B28])
                                  (*Pelteobagrus fulvidraco*)                                                                                                                                                                                                                                           

  *Edwardsiella tarda*            Japanese flounder\              Blood                                A total of 30 immune related DEGs were extracted based on protein--protein interaction networks.                                                                                                 ([@B29])
                                  (*Paralichthys olivaceus*)                                                                                                                                                                                                                                            

  *Edwardsiella tarda*            Japanese flounder\              Gill                                 Based on protein--protein interaction networks, 24 genes were found as key regulators in immune responses.                                                                                       ([@B30])
                                  (*Paralichthys olivaceus*)                                                                                                                                                                                                                                            

  *Edwardsiella tarda*            Japanese flounder\              Kidney                               PI3K family, antigen processing and presentation.                                                                                                                                                ([@B31])
                                  (*Paralichthys olivaceus*)                                                                                                                                                                                                                                            

  *Edwardsiella tarda*\           Zebrafish\                      Liver                                Complement activation, antigen processing and presentation pathway.                                                                                                                              ([@B32])
  (mutant with low virulence)     (*Danio rerio*)                                                                                                                                                                                                                                                       

  *Escherichia coli*\             Black rockcod\                  Liver                                MHC class I and II antigen-processing moleules.                                                                                                                                                  ([@B33])
  (Heat-killed)                   (*Notothenia coriiceps*)                                                                                                                                                                                                                                              

  *Flavobacterium columnare*      Channel catfish\                Gill                                 Innate immune genes\                                                                                                                                                                             ([@B34])
                                  (*Ictalurus punctatus*)                                              (iNOS2b, lysozyme C, IL-8, and TNF-α) and mucosal immuen factors\                                                                                                                                
                                                                                                       (CD103 and IL-17)                                                                                                                                                                                

  *Flavobacterium columnare*      Mandarin fish\                  Head kidney                          Upregulated cytokine genes\                                                                                                                                                                      ([@B35])
                                  (*Siniperca chuatsi*)                                                (IL-1β, IL-8, CXCL14, CCK3, CCL25, Angpt1, Thpo, MMP-9)                                                                                                                                          

  *Flavobacterium columnare*      Topmouth culter\                Head kidney                          Phagosome pathway, complement activation                                                                                                                                                         ([@B36])
                                  *(Culter alburnus)*                                                                                                                                                                                                                                                   

  *Piscirickettsia salmonis*      Atlantic salmon\                Brain, head kidney, spleen           Endocytosis, Phagosome, Heme metabolism                                                                                                                                                          ([@B37])
                                  (*Salmo salar*)                                                                                                                                                                                                                                                       

  *Pseudomonas fluorescens*       Whitefish\                      Embryo                               Ion binding, aminoacyl-tRNA-biosynthesis, complement cascade, MHC class I and II, TNF-α, T-cell differentiation.                                                                                 ([@B38])
                                  (*Coregnous palaea*)                                                                                                                                                                                                                                                  

  *Vibrio alginolyticus*          Orange-spotted grouper\         Whole body of larvae                 Complementation cascades, phagosome activity, antigen processing, antigen presentation pathway.                                                                                                  ([@B39])
                                  (*Epinephelus coioides*)                                                                                                                                                                                                                                              

  *Vibrio alginolyticus*          Giant grouper\                  Whole body of larvae                 Inductions of TLR5, IL-1β, IL-8, hepcidin genes                                                                                                                                                  ([@B40])
                                  (*Epinephelus lanceolatus*)                                                                                                                                                                                                                                           

  *Vibrio anguillarum*            Half-smooth tongue sole\        Mixture of RNA sample form tissues   Complement cascades, antigen processing and presentation, toll-like receptor signaling, NOD-like receptor signaling                                                                              ([@B41])
                                  (*Cynoglossus semilaevis*)                                                                                                                                                                                                                                            

  *Vibrio anguillarum*            Japanese seabass\               Head kidney, liver, spleen           Toll-like receptor signaling, NOD-like receptor signaling, RIG-I-like receptor signaling, T cell receptor signaling                                                                              ([@B42])
                                  (*Lateolabrax japonicas*)                                                                                                                                                                                                                                             

  *Vibrio anguillarum*            Miiuy croaker\                  Spleen                               Phagosome pathway, NF-κB signaling pathway, hematopoietic cell lineage, cytokine--cytokine receptor interaction.                                                                                 ([@B43])
                                  (*Miichthys miiuy*)                                                                                                                                                                                                                                                   

  *Vibrio anguillarum*            Southern flounder\              Liver                                Upregulated hematopoiesis related genes, β-hemoglobin and erythropoiten.                                                                                                                         ([@B44])
                                  (*Paralichthys lethostigma*)                                                                                                                                                                                                                                          

  *Vibrio anguillarum*            Turbot\                         Intestine                            Pathogen attachment and recognition, antioxidant/apoptosis, mucus barrier modification, immune activation/inflammation                                                                           ([@B45])
                                  (*Scophthalmus maximus*)                                                                                                                                                                                                                                              

  *Vibrio harveyi*                Orange-spotted grouper\         Head kidney, spleen                  NOD-like receptor signaling, toll-like receptor signaling, NF-κB signaling, Jak-STAT signaling                                                                                                   ([@B16])
                                  (*Epinephelus coioides*)                                                                                                                                                                                                                                              

  *Vibrio parahaemolyticus*       Brown-marbled grouper\          Spleen                               Cell killing, antioxidant activity                                                                                                                                                               ([@B46])
                                  (*Epinephelus fuscoguttatus)*                                                                                                                                                                                                                                         

  *Vibrio parahaemolyticus*       Zebrafish\                      Whole body of larvae                 Complement cascades, chemokine, TNF signaling, NF-κB signaling, JAK-STAT signaling                                                                                                               ([@B47])
                                  (*Danio rerio*)                                                                                                                                                                                                                                                       

  *Yersinia ruckeri*              Amur sturgeon\                  Spleen                               Antigen processing and presentation, complement cascades, T cell receptor signaling, B and T cell receptor signaling.                                                                            ([@B48])
                                  (*Acipenser schrenckii*)                                                                                                                                                                                                                                              

  **GRAM-POSITIVE BACTERIA**                                                                                                                                                                                                                                                                            

  *Lactococcus garvieae*          Gray mullet\                    Head kidney, spleen                  Complement and coagulation cascades, toll-like receptor signaling, antigen processing and presentation                                                                                           ([@B18])
                                  (*Mugil cephalus*)                                                                                                                                                                                                                                                    

  *Nocardia seriolae*             Largemouth bass\                Spleen                               Cytokines, toll-like receptor signaling, T cell receptor signaling, NF-κB signaling, JAK-STAT signaling                                                                                          ([@B17])
                                  (*Micropterus salmoides*)                                                                                                                                                                                                                                             

  *Streptococcus agalactiae*      Hybrid tilapia\                 Whole body of larvae                 Toll-like receptor signaling, leukocyte transendothelial migration, TNF signaling, PI3K-Akt signaling, Jak-STAT signaling MAPK signaling                                                         ([@B49])
                                  (*Oreochromis* spp.)                                                                                                                                                                                                                                                  

  *Streptococcus agalactiae*      Mozambique tilapia\             Spleen                               Toll-like receptor signaling, chemokine signaling, antigen processing and presentation, NF-κB signaling, TNF signaling, cytokines                                                                ([@B50])
                                  (*Oreochromis mossambicus*)                                                                                                                                                                                                                                           

  *Streptococcus agalactiae*      Nile tilapia\                   Spleen                               Pathogen attachment and recognition, antioxidant/apoptosis, cytoskeletal dynamics regulation, immune activation                                                                                  ([@B51])
                                  (*Oreochromis niloticus*)                                                                                                                                                                                                                                             

  *Streptococcus dysgalactiae*    Soiny mullet\                   Spleen                               Complemen cascades, toll-like receptor signaling, chemokine signaling pathway, B snd T cell receptor signaling, antigen processing and presentation, natural killer cell mediated cytotoxicity   ([@B52])
                                  (*Liza haematocheila*)                                                                                                                                                                                                                                                

  **VIRUS**                                                                                                                                                                                                                                                                                             

  Carassius auratus herpesvirus   Gibel carp\                     Head kidney                          Chemokine signaling, toll-like receptor signaling, type I IFN, IL-6 etc.                                                                                                                         ([@B53])
                                  (*Carassius gibelio*)                                                                                                                                                                                                                                                 

  Grass carp reovirus             Grass carp\                     Gill, intestine, liver, spleen       Complemen cascades, antigen presentation, Nk cell mediated cytotoxicity, proteasome, lysosome, peroxisome, phagosome                                                                             ([@B54])
                                  (*Ctenopharyngodon idella*)                                                                                                                                                                                                                                           

  Grass carp reovirus             Grass carp\                     Spleen                               Complement cascades, hematopoietic cell lineage, phagosome, cytokine-cytokine receptor interaction                                                                                               ([@B55])
                                  (*Ctenopharyngodon idella*)                                                                                                                                                                                                                                           

  Grass carp reovirus             Grass carp\                     Kidney                               Complement and coagulation cascades, Glycolysis/gluconeogenesis                                                                                                                                  ([@B56])
                                  (*Ctenopharyngodon idellus*)                                                                                                                                                                                                                                          

  Grass carp reovirus             Grass Carp\                     Kidney cell line                     Focal adhesion, ECM--receptor interactions in early stage of infection, phagosome, lysosome in the midle stage                                                                                   ([@B57])
                                  (*Ctenopharyngodon idella*)                                                                                                                                                                                                                                           

  Lymphocystis disease virus      Japanese flounder\              Gill                                 Cell cycle, DNA replication, proteasome, p53 signaling, TNF signaling                                                                                                                            ([@B58])
                                  (*Paralichthys olivaceus*)                                                                                                                                                                                                                                            

  Nervous necrosis virus          Grouper\                        Kidney cell line                     Endoplasmic reticulum stress response                                                                                                                                                            ([@B59])
                                  (*Epinephelus spp*)                                                                                                                                                                                                                                                   

  Nervous necrosis virus          Asian seabass\                  Epithelial cell line                 Upregulation of pro-inflammatory cytokines, type I interferon, chemokines                                                                                                                        ([@B60])
                                  (*Lates calcarifer*)                                                                                                                                                                                                                                                  

  Salmon anemia virus             Atlantic salmon\                Spleen                               Protein translation; protein processing in endoplasmic reticulum, proteasome, phagosome, lysosome, antigen processing and presentation                                                           ([@B61])
                                  (*Salmo salar*)                                                                                                                                                                                                                                                       

  Spring viraemia of carp virus   Fathead minnow\                 Epithelioma papulosum cyprinid\      Oxidative stress, apoptosis, cytoskeleton, interferon system                                                                                                                                     ([@B62])
                                  (*Pimephales promelas*)         (EPC) cell line                                                                                                                                                                                                                       

  Spring viraemia of carp virus   Zebrafish\                      Brain, spleen                        Influenza A pathway, herpes simplex infection pathway, tuberculosis and toxoplasmosis pathways                                                                                                   ([@B63])
                                  (*Danio rerio*)                                                                                                                                                                                                                                                       

  **MIMIC MOLECULES**                                                                                                                                                                                                                                                                                   

  Lipopolysaccharide              Zebrafish\                      Whole body of larvae                 Chemokines and G protein-coupled receptor signaling.\                                                                                                                                            ([@B64])
                                  (*Danio rerio*)                                                      (mild dose treatment). Down-modulated proinflammatory genes after sublethal dose treatment.                                                                                                      

  Lipopolysaccharide              Yellow catfish\                 Liver                                Up-regulated expressions of CXCL2-like chemokine, goose-type lysozyme, and cathepsin K                                                                                                           ([@B65])
                                  (*Pelteobagrus fulvidraco*)                                                                                                                                                                                                                                           

  Poly(I:C)                       Large yellow croaker\           Spleen                               Toll-like receptor signaling, RIG-I-like receptors signaling, JAK-STAT signaling, and T-cell receptor signaling                                                                                  ([@B66])
                                  (*Pseudosciaena crocea*)                                                                                                                                                                                                                                              

  Poly(I:C)                       Miiuy croaker\                  Spleen                               Cytokine--cytokine receptor interaction, complement cascades, NF-κB signaling, toll-like receptor signaling                                                                                      ([@B67])
                                  (*Miichthys miiuy*)                                                                                                                                                                                                                                                   

  Poly(I:C)                       Black rockcod\                  Liver                                Upregulation of TNFα and TNF2 expressions                                                                                                                                                        ([@B33])
                                  (*Notothenia coriiceps*)                                                                                                                                                                                                                                              

  Poly(I:C)                       *Schizothorax prenanti*         Spleen                               MDA5 and JAK mediated signaling pathways                                                                                                                                                         ([@B68])
  ----------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------

*Aeromonus hydrophila* is a Gram-negative bacterium, and causes a wide spectrum of diseases in vertebrates ([@B69], [@B70]). It is a major pathogen in aquaculture farms, and leads to high mortalities and economic losses worldwide ([@B71], [@B72]). In blunt snout bream (*Megalobrama amblycephala*), RNA-Seq analysis was conducted with RNA from several tissues, and 238 differentially expressed unigenes were identified in infected fish ([@B22]). In grass carp (*Ctenopharygodon idella*), 2121 DEGs were identified in spleens of *A. hydrophila* (6 hpi)-infected fish, some of which were involved in phagocytosis, the complement system, and cytokine production ([@B25]). Using transcriptome analysis, another study showed that *A. hydrophila* infected grass carp exhibited 2992 DEGs in the spleen, which were associated with the complement and coagulation cascades ([@B26]). In golden mahseer (*Tor putitora*), DEGs in *A. hydrophila*-infected livers were mainly associated with Th1/2 cell differentiation pathways, as well as in pathogen recognition and complement system ([@B24]).

*Flavobacterium columnare* is a Gram-negative bacterium, and causes columnaris in freshwater fish ([@B73]). This disease induces pathological changes, and damages epidermal tissues, gills, and the skin ([@B74]). In channel catfish (*Ictalurus punctatus*), the transcript profile of *F. columnare-*infected gills was examined using RNA-Seq to investigate differences in susceptibility to *F. columnare* ([@B34]). In resistant fish, the expression level of innate immune-associated genes (iNOS2b, lysozyme C, IL-8, and TNFα) was found to be elevated. In susceptible fish, the expression of secreted mucin forms, mucosal immune factors (CD103 and IL-17a), and rhamnose-binding lectin ([@B34]) was upregulated. The transcriptomic profiles of *F. columnare*-infected and non-infected mandarin fish (*Siniperca chuatsi*) have been reported using the head kidney *F. columnare*-infected and non-infected group ([@B35]). The results indicated that 1019 genes were differentially expressed between the two groups, of which 27 were immune-related ([@B35]). A similar study using the head kidney *F. columnare*-infected topmouth culter (*Culter alburnus*) ([@B36]) was also conducted. A total of 4037 DEGs (1217 upregulated and 2820 downregulated genes) were identified, and were found to be involved in phagosome formation, carbohydrate metabolism, amino acid metabolism, and lipid metabolism ([@B36]).

*Streptococcus agalactiae*, a Gram-positive round bacterium, is a harmful aquaculture pathogen that leads to enormous economic losses in various teleost ([@B75]--[@B78]). Transcriptome analysis of hybrid tilapia (*Oreochromis* spp.) after *S. agalactiae* infection was conducted, and results indicated that DEGs are mainly involved in immune-related pathways, especially Toll-like receptor signaling and leukocyte transendothelial migration ([@B49]). Moreover, time-course expression profile of genes suggested that induction of the NADPH oxidase complex and piscidin is mediated by Toll-like receptor pathways ([@B49]). Another research group conducted RNA-Seq analysis in tilapia (*Oreochromis niloticus*) spleens following *S. agalactiae* infections ([@B51]). A total of 2822 DEGs were detected, many of which were involved in pathogen attachment and recognition, antioxidant/apoptosis, cytoskeletal rearrangement, and immune activation ([@B51]). Wang et al. ([@B50]) focused on the relation between temperature and bacterial infection. They showed that temperature influences mRNA profiles of the spleen in tilapia during *S. agalactiae* infections. In addition, it was suggested that DEGs are involved in immune responses and oxygen related metabolisms ([@B50]).

*Vibrio alginolyticus* is a halophilic Gram-negative bacterium that causes septicemias, ulcers, exophthalmia, and corneal opaqueness in marine fish worldwide ([@B79], [@B80]). Transcriptome analysis in larvae of orange-spotted grouper (*Epinephelus coioides*) revealed that the expression of genes involved in the complement pathway and antimicrobial peptides is enhanced upon *V. alginolyticus* infection ([@B39]). In addition, transcriptome profiles of giant grouper (*Epinephelus lanceolatus*) larvae infected with *Vibrio alginolyticus* suggested that TLR5 signaling induces secretion of several cytokines (IL-1β and IL-8) ([@B40]).

Diversity of Immune Responses Among Species and Pathogens {#s5}
=========================================================

In the previous section, we introduced various RNA-seq analyses conducted in fish with bacterial infections. We have also previously published **four** research papers that conducted transcriptome analysis on infected fish, namely largemouth bass (*Micropterus salmoides*) against *Nocardia seriolae* ([@B17]), gray mullet (*Mugil cephalus*) against *Lactococcus garvieae* ([@B18]), orange-spotted grouper (*Epinephelus coioides*) against *Vibrio harveyi* ([@B16]), and koi carp (*Cyprinus carpio*) against *Aeromonas sobria* ([@B19]). Based on the transcriptome data from these reports, we gained a deeper understanding of immune responses to bacterial infections. However, there is little information regarding the universality and diversity of immune reactions of fish against pathogenic infections. Here, we investigated specific genes and pathways that are involved in each bacterial infection in various fish species. In this study, we used DEGs (transcripts from spleen at 1 dpi with log2 \> 1 or \< −1 between infected and control group) with KEGG-annotations. We **first** identified overlapping and specific genes that were up- or down- regulated in each species. Venn diagrams ([Figure 1](#F1){ref-type="fig"}) showed that only 39 DEGs (25 up-regulated and 14 down regulated) were involved in all species. The number of specific DEGs in each species was relatively higher than that of common DEGs; 493 DEGs (167 up-regulated and 326 down regulated) were found in largemouth bass against *N. seriolae*, 819 DEGs (291 up-regulated and 528 down regulated) were found in mullets against *L. garvieae*, 909 DEGs (601 up-regulated and 308 down regulated) were found in groupers (*Epinephelus coioides*) against *V. harveyi*, and 1471 DEGs (1,001 up-regulated and 470 down regulated) were found in carps against *A. sobria* ([Figure 1](#F1){ref-type="fig"}).

![Venn diagrams showing overlaps of up and down regulated genes among each fish with bacterial challenge. The numbers indicate up (red arrow) and down (blue arrow) regulated genes in each categories.](fimmu-10-00153-g0001){#F1}

Of the common DEGs, we found several immune-related genes that were upregulated, including C4 (complement component 4), CCL19 (C-C motif chemokine 19), and SOCS1 (suppressor of cytokine signaling 1) ([Table S1](#SM4){ref-type="supplementary-material"}). The complement system is an important innate immune system that functions to detect pathogenic infections in both vertebrates and invertebrates. C4 is an important part of the classical and lectin pathways, which form enzymes C3 and C5 convertases ([@B81], [@B82]). CCL19, a CC chemokine that is expressed in lymphoid organs, manages the migration of antigen presenting cells and lymphocytes ([@B83]). In teleost, there are also various reports that investigated potential chemokine genes and their chemotactic activity ([@B84]--[@B86]). SOCS1 is a regulator of JAK/STAT signaling, and is induced by type I interferon (IFN) and IFN-γ via binding and blocking of JAK2 activation ([@B87]). It has been reported that SOCS1 acts as an inhibitor of IFN-mediated signaling in Atlantic salmon (*Salmo salar*) ([@B88]). From other reports of the transcriptome analysis ([Table 1](#T1){ref-type="table"}), complement system, JAK/STAT signaling and chemokine systems are also commonly appeared in responding pathways to bacterial infections. Therefore, it is suggested that these genes contribute to early immune responses following bacterial infections (within 24 h).

*Nocardia seriolae* is a filamentous Gram-positive bacterium that causes nocardiosis with high mortality in many fish species in Japan, Taiwan and Japan. The infected fish showed a lethal granulomatous disease of the skin, muscle, spleen, kidney, and liver tissues ([@B89]). Unlike other bacterial species from our previous studies, *N. seriolae* is an intercellular bacteria. To determine specific DEGs elicited by *N. seriolae* infections, we performed functional enrichment of the KEGG pathway for specific up-regulated genes in largemouth bass. As shown in [Table S2](#SM5){ref-type="supplementary-material"}, specific upregulated genes were assigned to 11 KEGG pathways; based on the enrichment analysis. From the enrichment analysis, Notch signaling pathway was focused and illustrated using expression levels of RNA-seq data from all four fish species ([Figure 2A](#F2){ref-type="fig"}). Results indicated that Notch1 and HES1 (hairy and enhancer of split 1) were specifically upregulated in largemouth bass against *N. seriolae*. The Notch and HES1 axis present in hematopoietic cells and stroma of the thymus plays an important role in T cell development ([@B90], [@B91]). Additionally, in the "cytokine-cytokine receptor interaction" pathway, IL12RB1 (interleukin 12 receptor β-1) and IL12RB2 (interleukin 12 receptor β2) in largemouth bass against *N. seriolae* were upregulated. IL12B is a ligand of IL12RBs, and is highly, but not specifically expressed, in largemouth bass ([Figure 2B](#F2){ref-type="fig"}). IL-12, a heterodimetric cytokine consisting of p35 and p40 subunits, is a key regulator of T helper 1 development (Th1), which promotes cellular immunity against intracellular pathogens. In Amberjack (*Seriola dumerili*), administrated recombinant IL-12 and formalin-killed *N. seriolae* showed the higher survival rates after challenged with *N. seriolae*, compared to vehicle and FKC only groups ([@B92]). These pathways promote immune reactions against *N. seriola*e during early stages of the infection, and are candidates for infection prevention and adjuvants in fish.

![Pathway map of Notch signaling **(A)** and Th1 differentiation **(B)** in KEGG. In each gene boxes, the gene expression levels are shown in 4 fish (C, Carp; G, Grouper; L, Largemouth bass; M, Mullet) spleen 1 day after infection with *A. sobria, V. harveyi, N. seriolae*, and *L. garvieae*, respectively, when compared to the control group. The lower expression levels of genes are shown in green, and the higher expression levels of genes are shown in red. Undetected genes are shown by white coloring (see color legend in figure).](fimmu-10-00153-g0002){#F2}

*Aeromonas sobria* is a Gram-negative, motile, rod-shaped bacterium that has been isolated from many diseased fish ([@B93]--[@B95]). In the spleen and head kidney of disease fish, necrotized spleen cells and hemorrhagic pulps were observed ([@B93]). From the extracted data of specific up-regulated genes (1001 genes) in koi carp against *A sobria* ([Figure 1](#F1){ref-type="fig"}), we performed functional enrichment of the KEGG pathway for specific up-regulated genes. As shown in [Table S3](#SM6){ref-type="supplementary-material"}, specific upregulated genes of koi carp against *A. sobria* are associated with 45 KEGG pathways. As shown in [Figure S1](#SM1){ref-type="supplementary-material"}, regulation of the actin cytoskeleton was activated during *A. sobria* infection. In addition, CXCL12 (C-X-C motif chemokine 12) and CXCR4 (C-X-C motif chemokine receptor 4) in were also up-regulated in koi carps during *A sobria* infections ([Figure S1](#SM1){ref-type="supplementary-material"}). The CXCL12-CXCR4 axis modulates various immune functions, such as induction of hematopoiesis and accumulation of immune cells in inflamed tissues ([@B96]). Therefore, the CXCL12-CXCR4 axis may function in reorganization of hematopoiesis in injured tissues during *A. sobria* infections.

*Vibrio harveyi* is one of the major photogenes of a luminescent Gram-negative bacterium, which impacts to wide range of aquaculture species ([@B97]--[@B99]). The 601 specific upregulated genes in orange-spotted grouper against *Vibrio harveyi* ([Figure 1](#F1){ref-type="fig"}) were assigned to 8 KEGG pathways ([Table S4](#SM7){ref-type="supplementary-material"}). We focused on the ErbB signaling pathway, and found that expression of TGFα (transforming growth factor α) and its receptor, ERBB1 (epidermal growth factor receptor), were upregulated ([Figure S2](#SM2){ref-type="supplementary-material"}). Previous studies have shown that TGFα promotes the expression and activity of TLR5 and TLR9 in skin keratinocytes ([@B100]). In our previous study, expressions of TRL5 and its downstream genes in the spleen were found to be enhanced 2 days following *V. harveyi* infections ([@B16]). While the immunological function of TGFα in the spleen of fish is unclear, we hypothesize that TGFα is a key regulator for prevention of *V. harveyi* infection in fish.

*Lactococcus garvieae* is a Gram-positive, facultative anaerobic, non-motile bacterium, and affects freshwater and marine cultured fish species worldwide ([@B101], [@B102]). Functional enrichment analysis of the KEGG pathway was performed to determine specific upregulated genes (291 genes) in gray mullets against *L*. *garvieae* ([Figure 1](#F1){ref-type="fig"}). Specific upregulated genes were mapped to 10 KEGG pathways during *L*. *garvieae* infection in gray mullets ([Table S5](#SM8){ref-type="supplementary-material"}). Results indicated that the IL-17 signaling pathway is clearly enhanced during the infection, as illustrated by [Figure S3](#SM3){ref-type="supplementary-material"}. IL-17 is composed of six ligands (IL-17A to F), and plays critical roles in inflammatory responses and host defenses during invasion by extracellular pathogens. Binding of IL-17s to its five perspective receptors (IL-17Rs; IL-17RA to E), induces inflammatory and immune responses ([@B103]). In *vitro* study, exogenous IL-17A induced bacterial clearance in F. tularensis LVS- live vaccine strain infected cells ([@B104]). Additionally, in mice model, it has been reported that *in vivo* administration of IL-17A moderately delays time of death from lethal infection of *Francisella tularensis* live vaccine strain ([@B105]). While we did not detect expression of IL-17 ligands in this study, we found that expressions of IL-17RB, IL17RC, and IL-17RE were up-regulated in gray mullets infected with *L*. *garvieae*. There are studies that aimed to identify and characterize IL-17 and IL17Rs in fish ([@B106], [@B107]). However, functional differentiations of teleost IL-17s and these receptors remain elusive. Our findings on the expression pattern of IL17Rs will provide useful models that can be used to investigate immune functions of IL17s in teleost.

*A. sobria* and *V. harveyi* are classified to Gram-negative bacteria. Therefore, we approached to find the immune-related genes and pathways using commonly DEG (379 upregulated genes and 91 downregulated genes) of koi carp against *A. sobria* and orange-spotted grouper against *V. harveyi* ([Figure 1](#F1){ref-type="fig"}). However, any immune-related pathways were not assigned by KEGG enrichments analysis. TLR4, which is the pathogen recognized receptor for the Gram-negative bacteria specific lipopolysaccharide, is not highly expression in the spleen of koi carp against *A. sobria* and orange-spotted grouper against *V. harveyi*. While, we could find the up-regulated immune-related gene of these two species, such as pattern recognition receptors (TLR6 and TLR5), cytokines and chemokines (CSF3 and CCL21), lysosome related genes (LYPLA3 and SLC11A1) and caspase recruitment domain-containing protein (Card) 9 ([Table S6](#SM9){ref-type="supplementary-material"}). *N. seriolae* and *L. garvieae* are classified to Gram-positive bacteria. Therefore, we investigated to commonly immune-related genes and pathways in Gram-positive bacteria using commonly DEGs (32 upregulated genes and 63 downregulated genes) of largemouth bass against *N. seriolae* and gray mullet against *L. garvieae* ([Figure 1](#F1){ref-type="fig"}). Up-regulated immune-related genes of these two species were identified, such as IL-6, TNF Receptor Superfamily Member 11b (TNFRSF11B), interferon regulatory factors (IRF4 and IRF8) and CD83 ([Table S7](#SM10){ref-type="supplementary-material"}). Although it is unclear the pathways to induce these up-regulated genes, these immune related genes may become the marker and immune factors in Gram-negative or positive bacterial infection.

Conclusion and Future Perspectives {#s6}
==================================

In this review, we first introduced applications of the RNA-Seq technology in aquaculture studies. The RNA-Seq technology has allowed us to identify many novel genes, and to investigate the expression patterns at various conditions in non-model teleost. Therefore, findings based on this technology have accelerated research in the aquaculture field. Additionally, high throughput quantification by RNA-seq could be used to identify pathogen, and to evaluate the efficacy of vaccines and adjuvants against pathogen *in vivo*. We also summarized current knowledge on immune responses to pathogenic challenges via RNA-Seq in teleost. In this study, we could identify the specific pathway in each fish against bacteria species, Noth1 signaling and IL-12 signaling pathway in largemouth bass against *N. seriolae*, CXCL12 and CXCR4 signaling pathway in koi carp against *A. sobria*, TGFα signaling pathway in orange-spotted grouper against *V. harveyi*, and IL-17 signaling pathway in gray mullet against *L. garvieae*. These types of studies are increasing, and have enormously aided in our understanding of pathogenic strategies and immune defense systems in aquaculture fish. However, there remains certain limitations of RNA-Seq analysis in aquaculture species. Additionally, the RNA-seq technology could be used to expand existing datasets on splicing variants in mRNA and SNP. Currently, differences in immune response against different pathogens are not well-described. In this study, we attempted to investigate both species-specific and common immune related genes that are up-regulated during bacterial infections based on our previous RNA-seq data. Secondary use of RNA-seq datasets may be essential for preparation of future RNA-seq studies in aquaculture species, which can further deepen our understanding of specific immune functions against pathogens. In aquaculture field, these deep and particular understanding of immune response against each pathogens will provide us to more accurate diagnosis of disease and develop a more effective vaccine and adjuvant of each pathogens.

Author Contributions {#s7}
====================

SM analyzed transcriptome data and wrote the paper. P-CW and S-CC reviewed the paper.

Conflict of Interest Statement
------------------------------

The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.

We thank Genomics Bioscience Technology Co. Ltd. (Taipei, Taiwan) for assistance with transcriptome analysis.

**Funding.** This research was funded by National Science Council, Taiwan, grant number NSC 104-2622-B-020-002-CC1 and MOST 107-2313-B-020−012 -MY3.

Supplementary Material {#s9}
======================

The Supplementary Material for this article can be found online at: <https://www.frontiersin.org/articles/10.3389/fimmu.2019.00153/full#supplementary-material>

###### 

Pathway map of Regulation of actin cytoskeleton in KEGG. In each gene boxes, the gene expression levels are shown in 4 fish (C, Carp; G, Grouper; L, Largemouth bass; M, Mullet) spleen 1 day after infection with *A. sobria, V. harveyi, N. seriolae*, and *L. garvieae*, respectively, when compared to the control group. The lower expression levels of genes are shown in green, and the higher expression levels of genes are shown in red. Undetected genes are shown by white coloring (see color legend in figure).

###### 

Click here for additional data file.

###### 

Pathway map of ErbB signaling pathway in KEGG. In each gene boxes, the gene expression levels are shown in 4 fish (C, Carp; G, Grouper; L, Largemouth bass; M, Mullet) spleen 1 day after infection with *A. sobria, V. harveyi, N. seriolae*, and *L. garvieae*, respectively, when compared to the control group. The lower expression levels of genes are shown in green, and the higher expression levels of genes are shown in red. Undetected genes are shown by white coloring (see color legend in figure).

###### 

Click here for additional data file.

###### 

Pathway map of Il-17 signaling pathway in KEGG. In each gene boxes, the gene expression levels are shown in 4 fish (C, Carp; G, Grouper; L, Largemouth bass; M, Mullet) spleen 1 day after infection with *A. sobria, V. harveyi, N. seriolae*, and *L. garvieae*, respectively, when compared to the control group. The lower expression levels of genes are shown in green, and the higher expression levels of genes are shown in red. Undetected genes are shown by white coloring (see color legend in figure).

###### 

Click here for additional data file.

###### 

Expression levels of commonly up or down regulated genes among 4 fish infected with bacteria. Each expression data indicated fold-change (log2) of bacterial infection/control groups.

###### 

Click here for additional data file.

###### 

Functional enrichment KEGG pathway for specific up-regulated genes of largemouth bass against *N. seriolae*.

###### 

Click here for additional data file.

###### 

Functional enrichment KEGG pathway for specific up-regulated genes of koi carp against *A. sobria*.

###### 

Click here for additional data file.

###### 

Functional enrichment KEGG pathway for specific up-regulated genes of orange-spotted grouper against *V. harveyi*.

###### 

Click here for additional data file.

###### 

Functional enrichment KEGG pathway for specific up-regulated genes of gray mullet against *Lactococcus garvieae*.

###### 

Click here for additional data file.

###### 

Expression levels of commonly up regulated immune-related genes among koi carp against *A. sobria* and orange-spotted grouper against *V. harveyi*. Each expression data indicated fold-change (log2) of bacterial infection/control groups.

###### 

Click here for additional data file.

###### 

Expression levels of commonly up regulated immune-related genes among largemouth bass against *N. seriolae* and gray mullet against *L. garvieae*. Each expression data indicated fold-change (log2) of bacterial infection/control groups.

###### 

Click here for additional data file.

[^1]: Edited by: Hetron Mweemba Munang\'andu, Norwegian University of Life Sciences, Norway

[^2]: Reviewed by: Ming-Wei Lu, National Taiwan Ocean University, Taiwan; Syarul Nataqain Baharum, National University of Malaysia, Malaysia

[^3]: This article was submitted to Comparative Immunology, a section of the journal Frontiers in Immunology
