Transcriptome profiling of spleen provides insights into the antiviral mechanism in Schizothorax prenanti after poly (I: C) challenge

Transcriptome-wide analyses of early immune responses in lumpfish leukocytes upon stimulation with poly(I:C)

Background

Both bacterial and viral diseases are a major threat to farmed fish. As the antiviral immune mechanisms in lumpfish (Cyclopterus lumpus L.) are poorly understood, lumpfish leukocytes were stimulated with poly(I:C), a synthetic analog of double stranded RNA, which mimic viral infections, and RNA sequencing was performed.

Methods

To address this gap, we stimulated lumpfish leukocytes with poly(I:C) for 6 and 24 hours and did RNA sequencing with three parallels per timepoint. Genome guided mapping was performed to define differentially expressed genes (DEGs).

Results

Immune genes were identified, and transcriptome-wide analyses of early immune responses showed that 376 and 2372 transcripts were significantly differentially expressed 6 and 24 hours post exposure (hpe) to poly(I:C), respectively. The most enriched GO terms when time had been accounted for, were immune system processes (GO:0002376) and immune response (GO:0006955). Analysis of DEGs showed that among the most highly upregulated genes were TLRs and genes belonging to the RIG-I signaling pathway, including LGP2, STING and MX, as well as IRF3 and IL12A. RIG-I was not identified, but in silico analyses showed that genes encoding proteins involved in pathogen recognition, cell signaling, and cytokines of the TLR and RIG-I signaling pathway are mostly conserved in lumpfish when compared to mammals and other teleost species.

Conclusions

Our analyses unravel the innate immune pathways playing a major role in antiviral defense in lumpfish. The information gathered can be used in comparative studies and lay the groundwork for future functional analyses of immune and pathogenicity mechanisms. Such knowledge is also necessary for the development of immunoprophylactic measures for lumpfish, which is extensively cultivated for use as cleaner fish in the aquaculture for removal of sea lice from Atlantic salmon (Salmo salar L.).

Spleen Transcriptome Profiling Reveals Divergent Immune Responses to LPS and Poly (I:C) Challenge in the Yellow Drum (Nibea albiflora)

The yellow drum (Nibea albiflora) is a marine teleost fish with strong disease resistance, yet the understanding of its immune response and key functional genes is fragmented. Here, RNA-Seq was used to investigate the regulation pathways and genes involved in the immune response to infection with lipopolysaccharide (LPS) and polyinosinic-polycytidylic acid (poly (I:C)) on the spleen of the yellow drum. There were fewer differentially expressed genes (DEGs) in the LPS-infected treatment group at either 6 or 48 h. Kyoto Encyclopedia of Genes and Genomes (KEGG) analysis showed that these DEGs were mainly significantly enriched in c5-branching dibasic acid metabolic and complement and coagulation cascades pathways. The yellow drum responded more strongly to poly (I:C) infection, with 185 and 521 DEGs obtained under 6 and 48 h treatments, respectively. These DEGs were significantly enriched in the Toll-like receptor signaling pathway, RIG-I-like receptor signaling pathway, Jak-STAT signaling pathway, NOD-like signaling pathway, and cytokine-cytokine receptor interaction. The key functional genes in these pathways played important roles in the immune response and maintenance of immune system homeostasis in the yellow drum. Weighted gene co-expression network analysis (WGCNA) revealed several important hub genes. Although the functions of some genes have not been confirmed, our study still provides significant information for further investigation of the immune system of the yellow drum.

Sequencing-based network analysis provides a core set of genes for understanding hemolymph immune response mechanisms against Poly I:C stimulation in Amphioctopus fangsiao

Aquatic viruses can spread rapidly and widely in seawater for their high infective ability. Polyinosinic-polycytidylic acid (Poly I:C), a viral dsRNA analog, is an immunostimulant that has been proved to activate various immune responses of immune cells in invertebrate. Hemolymph is a critical site that host immune response in invertebrates, and its transcriptome information obtained from Amphioctopus fangsiao stimulated by Poly I:C is crucial for understanding the antiviral molecular mechanisms of this species. In this study, we analyzed gene expression data in A. fangsiao hemolymph tissue within 24 h under Poly I:C stimulation and found 1082 and 299 differentially expressed genes (DEGs) at 6 and 24 h, respectively. Union set (1,369) DEGs were selected for subsequent analyses. Gene Ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) functional enrichment analyses were carried out for identifying DEGs related to immunity. Several significant immune-related terms and pathways, such as toll-like receptor signaling pathways term, inflammatory response term, TNF signaling pathway, and chemokine signaling pathway were identified. A protein-protein interaction (PPI) network was constructed for examining the relationships among immune-related genes. Finally, 12 hub genes, including EGFR, ACTG1, MAP2K1, and other nine hub genes, were identified based on the KEGG enrichment analysis and PPI network. The quantitative RT-PCR (qRT-PCR) was used to verify the expression profile of 12 hub genes. This research provides a reference for solving the problem of high mortality of A. fangsiao and other mollusks and provides a reference for the future production of some disease-resistant A. fangsiao.

Expression profiling and inflammatory activation analysis of high-mobility group box 1 in Schizothorax prenanti

OBJECTIVE

High-mobility group box 1 (HMGB1) is a highly conserved nuclear protein and participates in the immune response to pathogens in bony fish. In this study, the structure and function of HMGB1 in the cyprinid fish Schizothorax prenanti (SpHMGB1) were investigated.

METHODS

The spatial structure of SpHMGB1 was predicted by CPHmodels. Quantitative reverse transcription PCR was used to detect the mRNA of SpHMGB1 in different tissues and Streptococcus agalactiae infection. The macrophage was treated with synthetic SpHMGB1-B box peptide to analyze the inflammatory activity.

RESULT

Structurally, SpHMGB1 had the conserved A box, B box, and acid tail compared with Zebrafish Danio rerio and mice Mus musculus. SpHMGB1 was universally expressed in various tissues, with the highest expression in the middle kidney. In vivo, SpHMGB1 was significantly induced in response to Streptococcus agalactiae infection in the blood and spleen. Synthetic SpHMGB1-B box peptide activated respiratory burst and up-regulated the messenger RNA expression of interleukin-1β, tumor necrosis factor α, interleukin-10, interferon regulatory factor 1, interferon regulatory factor 7, C-X-C motif chemokine ligand 11-1, C-X-C motif chemokine ligand 11-2, and toll-like receptor 4 in macrophages.

CONCLUSION

This study suggested that SpHMGB1 participated in the response to bacterial pathogens and that SpHMGB1-B box peptide played an important role in mediating the immune response of S. prenanti.

Fish Innate Immune Response to Viral Infection-An Overview of Five Major Antiviral Genes

Fish viral diseases represent a constant threat to aquaculture production. Thus, a better understanding of the cellular mechanisms involved in establishing an antiviral state associated with protection against virus replication and pathogenesis is paramount for a sustainable aquaculture industry. This review summarizes the current state of knowledge on five selected host innate immune-related genes in response to the most relevant viral pathogens in fish farming. Viruses have been classified as ssRNA, dsRNA, and dsDNA according to their genomes, in order to shed light on what those viruses may share in common and what response may be virus-specific, both in vitro (cell culture) as well as in vivo. Special emphasis has been put on trying to identify markers of resistance to viral pathogenesis. That is, those genes more often associated with protection against viral disease, a key issue bearing in mind potential applications into the aquaculture industry.

Chitosan nanoparticle formulation attenuates poly (I:C) induced innate immune responses against inactivated virus vaccine in Atlantic salmon (Salmo salar)

Many vaccine formulations, in particular vaccines based on inactivated virus, needs adjuvants to boost immunogenicity. In aquaculture, mineral and plant oil are used as adjuvant in commercial vaccines, and the advent of oil-adjuvanted vaccines was crucial to aquaculture development. Nevertheless, some of these approved vaccines display suboptimal performance in the field compared to experimental conditions. Therefore, there is a need to improve adjuvants and delivery methods for fish vaccines against viruses. We used RNA sequencing of Atlantic salmon head kidney to analyse the difference in gene expression 24 h after injection of different experimental vaccine formulations. We compared five different formulations in addition to a PBS control: inactivated virus alone (group V), soluble poly (I:C) (group P), nanoparticles containing poly (I:C) (group N), soluble poly (I:C) + inactivated virus (group PV) and finally nanoparticles containing poly (I:C) + inactivated virus (group NV). Our results showed poly (I:C)'s ability as adjuvant and its capacity influence innate immune genes expression in Atlantic salmon. Soluble poly (I:C) upregulated multiple immune related genes and was more effective compared to poly (I:C) formulated into chitosan nanoparticles (more than 10 fold increase in differentially expressed genes, DEGs). However, inclusion of inactivated ISA virus in the nanoparticle vaccine, increased the number of DEGs fivefold suggesting a synergistic effect of adjuvant and antigen. Our results indicate that the way poly (I:C) is formulated and the presence of antigen is important for the magnitude of the innate immune response in Atlantic salmon.

Characterization of antiviral immune response induced by poly(I:C) in macrophages of farmed large yellow croaker (Larimichthys crocea)

Fish tend to rely more on their innate immunity to executing defense against viral infection by inducing antiviral gene production. However, the expression pattern and underlying mechanism of fish antiviral responses have yet to be fully defined. In the present study, an in vitro viral infection model was established by exposing head kidney-derived macrophages of large yellow croaker to virus analog, poly(I:C). Transcriptome analysis indicated that poly(I:C) appeared to induce potent antiviral activity featuring dominant interferon a3 (IFNa3) expression through activation of toll-like receptors (TLRs)/TIR-domain-containing adapter-inducing interferon-β (TRIF) and retinoic acid-inducible gene I-like receptors (RLRs)/mitochondrial antiviral signaling protein (MAVS) pathways. Inhibition of nuclear factor κB (NF-κB) and stimulator of interferon genes (STING)/interferon regulatory factor 3 (IRF3) pathways diminished the expression of IFNa3. Mechanistically, transcription factors including p65 and IRF3 could promote expression of IRF3, and activated IRF3 alone further increased the transcriptional activity of IFNa3. We also characterized the promoter of IFNa3 with direct IRF3 binding site which was sufficient to render the transcription of IFNa3. This effect was attenuated after deletion or mutation of the IRF3 binding sites. Taken together, our findings illustrate the distinct transcriptional profiling of fish macrophages triggered by poly(I:C). Also, this work provides new insights into the molecular mechanism underpinning coordinated activation of pathogen recognition and signaling transduction in the antiviral responses of non-model fish species.

Transcriptome-based evaluation and validation of suitable housekeeping gene for quantification real-time PCR under specific experiment condition in teleost fishes

Quantification real-time PCR (qRT-PCR) is a common method in analysis of gene expression, but the stable reference genes for the normalization analysis have not been appreciated before identifying expression pattern of genes in teleost fishes. In this study, we selected eight candidate reference genes (18S, Actin, EF-1α, 40S, B2M, TUBA, UBCE and GAPDH) basing on transcriptome analysis and the traditional housekeeping genes, and analyzed the stability of the reference genes in spleen, head kidney and head kidney leukocytes (HKL) after pathogen challenge in Schizothorax prenanti (S. prenanti). Three common programs (geNorm, NormFinder and Bestkeeper) were used to evaluate the stability of the candidate reference genes. Two reference genes, Actin and EF-1α presented higher stability, while 18S and GAPDH were the lower stable genes, both in in vitro and in vivo. An important immune gene, toll-like receptor 22a (TLR22a), was selected to validate the stability of the proposed reference genes (Actin and EF-1α) across different experiment treatments. The results reveal that Actin and EF-1α are quite suitable reference genes for the normalization analysis. Otherwise, using the most stable gene Actin to validate the reliable of transcriptome data showed the high correlation between the fold change of transcriptome data and qRT-PCR data. In conclusion, our study not only acquired the suitable reference gene for the qRT-PCR assay under specific experiment condition, but also provided a comprehensive method to evaluate and validate the reference gene based on transcriptome analysis in teleost fishes.

The potential sensing molecules and signal cascades for protecting teleost fishes against lipopolysaccharide

Lipopolysaccharide (LPS) is a classical pathogen-associated molecular pattern that can trigger strong inflammatory response mainly by TLR4-mediated signaling pathway in mammals, but the molecular mechanism of anti-LPS immunity is unclear in teleost fishes. In this study, we analyzed the gene expression features based on transcriptome analysis in Schizothorax prenanti (S. prenanti), after stimulation with two sources of LPS from Aeromonas hydrophila and Escherichia coli (Ah. LPS and Ecoli. LPS). 921 different expression genes (DEGs) after Ah. LPS stimulation and 975 DEGs after Ecoli.LPS stimulation were acquired, but only 706 and 750 DEGs were successfully annotated into the databases, respectively. Both of two groups of DGEs were significantly enriched into immune-related pathways by KEGG enrichment analysis, such as "Toll-like receptor signaling pathway", "Cytokine-cytokine receptor interaction" and "JAK-STAT signaling pathway". The annotated DEGs from Ah. LPS and Ecoli. LPS stimulation shared 470 DEGs, including 88 immune-related DEGs (IRGs) identified mainly by KEGG enrichment to immune-related signaling pathways. Among the shared IRGs, four pattern-recognition genes (TLR5, TLR25, PTX3 and C1q) were induced with high expression foldchange, and IFN-γ and relative genes also showed higher expression levels than control. Meanwhile, inflammatory signals were highlighted by upregulating the expression of inflammatory cytokines (IL-1β, IL-10 and IL-8). Moreover, some non-shared IRGs (including TLR2 and TLR4) were identified, suggesting that different sources of LPS own different potentials for the induction of immune gene expression. In conclusion, TLR5, TLR25, PTX3 and C1q may function as the sensing molecules to catch the invasion signal of LPS. The anti-LPS immune response may be involved into TLR25/TLR5-mediated inflammatory signals that regulate subsequently the activation of PTX3/C1q-modulated complement pathway upon the induction of PTX3 expression, and the crosstalk between IFN-γ and TLR signaling pathways in teleost fishes. This study will contribute to further explore the molecular mechanism of LPS-induced immunity in teleost fishes.

Arginine supplementation improves growth, antioxidant capacity, immunity and disease resistance of juvenile Chinese mitten crab, Eriocheir sinensis

To investigate the effects of arginine (Arg) on the growth, antioxidant capacity, immunity and disease resistance of juvenile Chinese mitten crab, three diets containing Arg levels at 1.72% (control), 2.73% and 3.72% were formulated and fed to Chinese mitten crab (0.22 ± 0.03 g) for eight weeks. The weight gain, ecdysterone and growth hormone in the serum, relative expression of insulin-like growth factor 2 in the hepatopancreas significantly increased in crabs fed the 2.73% and 3.72% Arg diets. The protein and lipid contents significantly increased in crabs fed the 3.72% Arg diet. The feed conversion ratios in crabs fed the diets with Arg additions were lower than in the control. Arg supplementation also enhanced the antioxidative capacity by increasing the activities of superoxide dismutase, catalase and the relative expression of Kelch-like ECH-associated protein 1 gene in the hepatopancreas, which subsequently decreased malondialdehyde content in the hepatopancreas. Besides, Arg also decreased nitric oxide content in the serum and the activity of nitric oxide synthetase in the hepatopancreas. The relative mRNA levels of crustin, relish, lysozyme and cryptocyanin genes were significantly upregulated by Arg supplementation. The activities of acid phosphatase and alkaline phosphatase in the serum significantly increased in crabs fed the 3.72% Arg diet than those in the control. Similarly, the relative mRNA levels of crustin, cryptocyanin and proPO genes were significantly upregulated in crabs fed the 2.73% Arg diet after lipopolysaccharide challenge, and in crabs fed the 3.72% Arg diet after the Poly (I:C) challenge. The crabs fed the 2.73% and 3.72% Arg diets had higher survival rate after bacterial infection than those fed the control diet. This study indicates that the addition of Arg to the diet at 2.7-3.7% can improve the growth, survival, antioxidant capacity, immunity and disease resistance in juvenile Chinese mitten crab.

The key differentially expressed genes and proteins related to immune response in the spleen of pufferfish (Takifugu obscurus) infected by Aeromonas hydrophila

The immune mechanism elicited in pufferfish (Takifugu obscurus) against the invasion of Aeromonas hydrophila is still poorly understood. We examined the spleen of pufferfish at the transcriptome and proteome levels by using Illumina-seq and TMT coupled mass spectrometry after 12 h infection by A. hydrophila, respectively. A total of 2,339 genes (1,512 up-regulated and 827 down-regulated) and 537 (237 up-regulated and 300 down-regulated) proteins were identified. GO and KEGG analyses revealed that the responses to stimulus were the main biological processes, intestinal immune network for IgT production and calcium signaling pathway. Fourteen genes (8 up-regulated and 6 down-regulated) and proteins (5 up-regulated and 9 down-regulated) involved immune responses or signal transduction were validated by qRT-PCR and parallel reaction monitoring to confirm the reliability of the transcriptomic and proteomic analyses, respectively. Moreover, qRT-PCR and flow cytometry were used to detect dynamics of the genes in calcium signaling pathway and changes of concentration of cytoplasm Ca²⁺ in spleen cells within a 72 h challenge. This study provides the findings regarding immune response, especially intestinal immune network for IgT production pathway and calcium signaling pathway at the molecular, protein and cellular in pufferfish after infection by A. hydrophila. These results would provide a new insight and molecular targets into the response to pathogenic infection in pufferfish.

Transcriptome analysis of differential expressed genes in hepatopancreas of Procambarus clarkii challenged with peptidoglycan

Procambarus clarkii is one of the most economically important species in Chinese aquaculture, and is widely cultured. Infection of P. clarkii populations with bacterial pathogens causes high mortality and great economic loss, therefore disease control is of significant economic importance. P. clarkii is a model system for studying immune responses in invertebrates, and its immune system consists solely of the innate response. In the present study, we examined gene expression related to immune function in P. clarkii in response to pathogen challenge. The transcriptome of hepatopancreas tissue from P. clarkii challenged with peptidoclycan (PGN) was analyzed and compared to control specimens. After assembly and annotation, 48,661 unigenes were identified with an average length of 671.54 bp. A total of 2533 differentially expressed genes (DEGs) were obtained, including 765 significantly up-regulated unigenes and 1757 significantly down-regulated unigenes. Gene ontology (GO) analysis demonstrated 19 biological process subcategories, 16 cellular component subcategories, and 17 molecular function subcategories that were enriched among these DEGs. Enrichment analysis using the Kyoto Encyclopedia of Genes and Genomes (KEGG) database revealed enrichment among immune responses pathways. Taken together, this study not only enriches the existing P. clarkii transcriptome database, but also elucidates immune responses of crayfish that are activated in response to PGN challenge.

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 analysis of spleen reveals the signal transduction of toll-like receptors after Aeromonas hydrophila infection in Schizothorax prenanti

Schizothorax prenanti (S. prenanti), an important species of economical fish in Southwest China, is susceptible to Aeromonas hydrophila (Ah). To understand the immune response to Ah, the transcriptome profiling of spleen of S. prenanti was analyzed after Ah infection. A total of 6, 213 different expression genes (DEGs) were obtained, including 3, 066 up-regulated DEGs and 3, 147 down-regulated DEGs. These DEGs were annotated by KEGG and GO databases, so that the immune-related DEGs (IRDs) can be identified and classified. Then, the interesting IRDs were screened to build heat map, and the reliability of the transcriptome data was validated by qPCR. In order to clarify the mechanism of signal transduction in the anti-bacterial immunity, the signaling pathway initiated by TLRs was predicted. In this pathway, TLR25 and TLR5 mediate the NF-κB and AP-1 signals via MyD88-dependent pathway. Meanwhile, the type I IFN (IFNα/β) induced by IRF1 and IRF3/7 may play an important role in the anti-bacterial immunity. In conclusion, this study preliminarily provides insights into the mechanism of signal transduction after Ah infection in S. prenanti, which contributes to exploring the complex anti-bacterial immunity.

Teleost-specific TLR25 identified from Schizothorax prenanti may recognize bacterial/viral components and activate NF-κB and type I IFNs signaling pathways

TLR25 is a new member of TLR1 family that is only identified in teleosts, but its function in immune response is still unclear. In current study, the coding sequence (CDS) of TLR25 was cloned from Schizothorax prenanti (named spTLR25), and spTLR25 is 2454 bp in length and coding a protein of 817 aa. The spTLR25 contains a signal peptide, twenty leucine-rich repeat (LRR) domains, a LRR C-terminal (LRRCT) motif, a transmembrane region and a Toll/IL-1 receptor (TIR) domain. Phylogenetic analysis indicates that spTLR25 has the closest relationship with Cyprinus carpio (C. carpio) TLR25-2. The 3D structure of spTLR25 exhibits 5 α-helices and 3 β-sheets in the TIR domain, and 8 α-helices and 6 β-sheets in the LRR domains. The spTLR25 is mainly expressed in immune-related tissues and peripheral blood leukocytes (PBL). Furthermore, the expression levels of spTLR25 were upregulated in spleen, head kidney and liver while S. prenanti was challenged with LPS or Aeromonas hydrophila (Ah), and the upregulation was also detected in head kidney leukocytes (HKL) after LPS and Poly (I:C) stimulation. The luciferase reporter assay demonstrated that NF-κB and type I IFNs signaling pathways can be activated by spTLR25, and this process may involve in the cascade amplification of TLR25-MyD88 signaling. In addition, the co-localization analysis showed that spTLR25 localizes to intracellular region. Taken together, our results reveal that teleost-specific TLR25 may be a multifunctional receptor for recognizing both LPS and Poly (I:C) and may activate NF-κB and type I IFNs signaling pathways.

Transcriptome analysis provides insights into the immune responsive pathways and genes in the head kidney of tiger grouper (Epinephelus fuscoguttatus) fed with Spatholobus suberectus, Phellodendron amurense, or Eclipta prostrata

The tiger grouper, Epinephelus fuscoguttatus, is an economically important fish in Southeast Asia but has been plagued by several diseases. Spatholobus suberectus (S), Phellodendron amurense (P), and Eclipta prostrate (E) are three commonly used Chinese medicinal herbs. Although previous pharmacological and clinical studies indicated that S, P, and E possess a variety of beneficial functions in mammals, little is known about their functions in farmed fish and the underlying molecular mechanism of their actions. Challenge tests in this study showed that after 14 days of diet supplement, all these herbs could effectively enhance the disease resistance of E. fuscoguttatus against Vibrio harveyi. However, the non-specific immune parameters of the herb-supplemented groups were not significantly different from the control group. To further explore the molecular mechanism of herbal immune-regulating effects on E. fuscoguttatus, transcriptome sequencing and RNA-Seq technique were applied on E. fuscoguttatus kidney. De novo transcriptome assembly of E. fuscoguttatus kidney yield 80,014 unigenes, among which, 44,901 (56.12%) were annotated with at least one of the public databases (Nr, Nt, Swiss-Prot, KEGG, COG, GO). Among these, 22,738, 11,700 and 27,457 unigenes were assigned to 57, 25 and 258 categories of GO, COG and KEGG databases, respectively. Using Solexa/Illumina's DGE platform, a total of 231, 186 and 144 putative differentially expressed genes (DEGs) were detected in P, E and S group compared with the control group. GO analysis indicated that in P and E, down-regulated DEGs were dominant in almost every GO term; whereas in S, up-regulated DEGs were more dominant. KEGG pathway analysis revealed that putative DEGs in all three herb groups were obviously enriched in the pathways related to infective diseases and immune system. We also identified a number of immune relative genes and pathways (TLR5, IL8 and MAPK pathway, for instance) associated with P, E and S's regulatory effects on E. fuscoguttatus. This study will enrich the E. fuscoguttatus transcriptome database, contribute to a better understanding of the molecular mechanisms associated with the immunoregulatory activities of Chinese medicinal herbs on teleost and provide valuable information on the prevention of grouper Vibrio diseases using Chinese medicinal herbs.

De novo transcriptome assembly and analysis of differential gene expression following lipopolysaccharide challenge in Pelteobagrus fulvidraco

The yellow catfish, Pelteobagrus fulvidraco, has been recognized as an important freshwater aquaculture species in Eastern and Southeast Asia. To gain a better understanding of the immune response in P. fulvidraco, we analyzed its transcriptome following stimulation with lipopolysaccharide (LPS). Phosphate buffer saline (PBS) was used as control. Following assembly and annotation, 72,152 unigenes with an average length of 1090 bp were identified. A total of 370 differentially expressed genes (DEGs) in the P. fulvidraco were observed at 12 h post LPS treatment, including 197 up-regulated genes and 173 down-regulated genes. Clusters of Orthologous Groups of proteins (KOG/COG) annotation demonstrated that a total of 18,819 unigenes classified into 26 categories. Gene ontology (GO) analysis revealed 20 biological process subcategories, 7 cellular component subcategories and 20 molecular function subcategories. The Kyoto Encyclopedia of Genes and Genomes (KEGG) enrichment analysis identified immune responses pathways. Quantitative reverse transcription polymerase chain reaction measured the expression of 18 genes involved in the immune response. CXCL2-like chemokine (CXCL2), goose-type lysozyme (LYZ G), and cathepsin K (CTSK) were significantly up-regulated. This study enriches the P. fulvidraco transcriptome database and provides insight into the immune response of P. fulvidraco against infection.

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.

Exploring gene expression changes in the amphioxus gill after poly(I:C) challenge using digital expression profiling

Amphioxus, a cephalochordate, is a key model animal for studying the evolution of vertebrate immunity. Recently, studies have revealed that microRNA (miRNA) expression profiles change significantly in the amphioxus gill after immune stimulation, but it remains largely unknown how gene expression responds to immune stress. Elucidating gene expression changes in the amphioxus gill will provide a deeper understanding of the evolution of gill immunity in vertebrates. Here, we used high-throughput RNA sequencing technology (RNA-seq) to conduct tag-based digital gene expression profiling (DGE) analyses of the gills of control Branchiostoma belcheri and of those exposed to the viral mimic, poly(I:C) (pIC). Six libraries were created for the control and treatment groups including three biological replicates per group. A total of 1999 differently expressed genes (DEGs) were obtained, with 571 and 1428 DEGs showing up- or down-regulation, respectively, in the treatment group. Enrichment analysis of gene ontology (GO) terms and pathways revealed that the DEGs were primarily related to immune and defense response, apoptosis, human disease, cancer, protein metabolism, enzyme activity, and regulatory processes. In addition, eight DEGs were randomly selected to validate the RNA-seq data using real-time quantitative PCR (qRT-PCR), and the results confirmed the accuracy of the RNA-seq approach. Next, we screened eight key responding genes to examine the dynamic changes in expression levels at different time points in more detail. The results indicated that expressions of TRADD, MARCH, RNF31, NF-κb, CYP450, TNFRSF6B, IFI and LECT1 were induced to participate in the antiviral response against pIC. This study provides a valuable resource for understanding the role of the amphioxus gill in antiviral immunity and the evolution of gill immunity in vertebrates.

Genome-wide gene expression analysis in the amphioxus, Branchiostoma belcheri after poly (I: C) challenge using strand-specific RNA-seq

The gene expression associated with immune response to bacteria/bacterial mimic has been extensively analyzed in amphioxus, but remains largely unknown about how gene are involved in the immune response to viral invasion at expression level. Here, we analyze the rRNA-depleted transcriptomes of Branchiostoma belcheri using strand-specific RNA-seq in response to the viral mimic, poly (I:C) (pIC). A total of 5,317 differentially expressed genes were detected at treatment group by comparing with control. The gene with the most significant expression changes (top 15) after pIC challenge and 7 immune-related categories involving 58 differently expressed genes were scrutinized. By functional enrichment analysis of differently expressed genes, gene ontology terms involving response to stress and stimulus, apoptosis, catabolic and metabolic processes and enzyme activity were overrepresented, and several pathways related to immune signaling, immune response, cancer, apoptosis, viral disease, metabolism were activated after pIC injection. A positive correlation between the qRT-PCR and strand-specific RNA-seq data confirmed the accuracy of the RNA-seq results. Additionally, the expression of genes encoding NLRC5, CASP1, CASP6, CYP450, CAT, and MDA5 were induced in B. belcheri under pIC challenge. Our experiments provide insight into the immune response of amphioxus to pIC and valuable gene expression information for studying the evolution of antiviral immunity in vertebrates.

Genomic and transcriptomic approaches to study immunology in cyprinids: What is next?

Accelerated by the introduction of Next-Generation Sequencing (NGS), a number of genomes of cyprinid fish species have been drafted, leading to a highly valuable collective resource of comparative genome information on cyprinids (Cyprinidae). In addition, NGS-based transcriptome analyses of different developmental stages, organs, or cell types, increasingly contribute to the understanding of complex physiological processes, including immune responses. Cyprinids are a highly interesting family because they comprise one of the most-diversified families of teleosts and because of their variation in ploidy level, with diploid, triploid, tetraploid, hexaploid and sometimes even octoploid species. The wealth of data obtained from NGS technologies provides both challenges and opportunities for immunological research, which will be discussed here. Correct interpretation of ploidy effects on immune responses requires knowledge of the degree of functional divergence between duplicated genes, which can differ even between closely-related cyprinid fish species. We summarize NGS-based progress in analysing immune responses and discuss the importance of respecting the presence of (multiple) duplicated gene sequences when performing transcriptome analyses for detailed understanding of complex physiological processes. Progressively, advances in NGS technology are providing workable methods to further elucidate the implications of gene duplication events and functional divergence of duplicates genes and proteins involved in immune responses in cyprinids. We conclude with discussing how future applications of NGS technologies and analysis methods could enhance immunological research and understanding.

Transcriptome-wide analysis of immune-responsive microRNAs against poly (I:C) challenge in Branchiostoma belcheri by deep sequencing and bioinformatics

Amphioxus is a key experimental animal for studying the evolution of vertebrate immune system. However, we still do not know about the roles of microRNAs (miRNAs) under viral stress in amphioxus. In this study, we sequenced six small RNA libraries (three biological replicates were included in the treatments challenged by the viral mimic, poly (I:C) (pIC) and control groups, respectively) from Branchiostoma belcheri. A total of 151 known miRNAs, 197 new miRNAs (named novel_mir, including nine conserved miRNAs) were identified by deep sequencing from the six libraries. We primarily focused on differentially expressed miRNAs (DEMs) after pIC challenge. Next, we screened a total of 77 DEMs, including 27 down- and 50 up-regulated DEMs in response to pIC challenge. Furthermore, we used real-time quantitative PCR (qRT-PCR) to verify the expression levels of 10 randomly selected DEMs. Target genes likely regulated by DEMs were predicted, and functional enrichment analyses of these targets were performed using bioinformatics approach. MiRNA targets of DEMs are primarily involved in immune response, diseases, cancer and regulation process, and could be largely linked to 14 immune-related signaling pathways, including NF-kappa B, NOD-like receptor, RIG-I-like receptor and endocytosis. The present study for the first time explores key regulatory roles of miRNAs in the innate antiviral immune response in amphioxus, and will provide insight into the molecular basis of antiviral immunity and evolution of immune-related miRNAs.