RNA-seq liver transcriptome analysis reveals an activated MHC-I pathway and an inhibited MHC-II pathway at the early stage of vaccine immunization in zebrafish
© Yang et al.; licensee BioMed Central Ltd. 2012
Received: 10 February 2012
Accepted: 17 July 2012
Published: 17 July 2012
Zebrafish (Danio rerio) is a prominent vertebrate model of human development and pathogenic disease and has recently been utilized to study teleost immune responses to infectious agents threatening the aquaculture industry. In this work, to clarify the host immune mechanisms underlying the protective effects of a putative vaccine and improve its immunogenicity in the future efforts, high-throughput RNA sequencing technology was used to investigate the immunization-related gene expression patterns of zebrafish immunized with Edwardsiella tarda live attenuated vaccine.
Average reads of 18.13 million and 14.27 million were obtained from livers of zebrafish immunized with phosphate buffered saline (mock) and E. tarda vaccine (WED), respectively. The reads were annotated with the Ensembl zebrafish database before differential expressed genes sequencing (DESeq) comparative analysis, which identified 4565 significantly differentially expressed genes (2186 up-regulated and 2379 down-regulated in WED; p<0.05). Among those, functional classifications were found in the Gene Ontology database for 3891 and in the Kyoto Encyclopedia of Genes and Genomes database for 3467. Several pathways involved in acute phase response, complement activation, immune/defense response, and antigen processing and presentation were remarkably affected at the early stage of WED immunization. Further qPCR analysis confirmed that the genes encoding the factors involved in major histocompatibility complex (MHC)-I processing pathway were up-regulated, while those involved in MHC-II pathway were down-regulated.
These data provided insights into the molecular mechanisms underlying zebrafish immune response to WED immunization and might aid future studies to develop a highly immunogenic vaccine against gram-negative bacteria in teleosts.
Zebrafish (Danio rerio) is a well-established vertebrate model of human development and hematopoiesis . However, as the aquaculture industry grows to meet the needs of an expanding global population, there has been a rapid increase in the research of teleost physiology and immunology. The zebrafish model system has proven to be a useful tool for studying infectious diseases that are natural threats to fish species of important human food sources, such as rock bream (Oplegnathus fasciatus) and turbot (Psetta maxima) . Several studies of the zebrafish immune system have provided initial insights into host-microbe interactions with both pathogens and commensals [3, 4], and the molecular processes mediating clearance of pathogenic infections [5–7]. Not surprisingly, many features of the zebrafish defense responses following pathogen invasion resembled those of other mammals, including humans [7–9].
The adult zebrafish has been used in a few studies to investigate new vaccines against specific pathogenic species. For example, Novoa et al. studied the efficacy of a vaccine derived by reverse genetics against viral hemorrhagic septicemia virus in fish by using the zebrafish as a model system, and revealed that the vaccine was protective even at low temperatures. In addition, Cui et al. used zebrafish to study an attenuated live Mycobacterium marinum vaccine aimed at reducing mycobacteriosis in freshwater and marine fish. Finally, Xiao et al. developed an edwardsiellosis zebrafish model to screen attenuated live Edwardsiella tarda vaccine candidates in order to identify those most highly effective for subsequent development for industry use. However, so far, no work involves the immune-related pathways underlying the zebrafish response to vaccination.
In order to design a novel and effective vaccine, it is essential to gain a comprehensive understanding of the immune responses elicited in host upon vaccination. To date, most of the studies of the teleost immune system have focused on head kidney or/and spleen . However, the vertebrate liver has recently been recognized as an essential immune organ [14–16], accommodating a variety of cell types [14, 17], including those primarily involved in immune activities. Since the liver receives blood from both the systemic circulation and the intestine, it is exposed to a wide array of antigens. Therefore, its immune-related cellular components can manifest a broad range of immune reactions [14, 18]. For example, the liver lymphocyte population includes both innate immune cells (such as the macrophages and natural killer cells) and adaptive immune cells (such as conventional T lymphocytes [14, 19] that recognize and respond to antigenic peptides presented by the major histocompatibility complex (MHC)-I or -II [18, 20]). As such, different infectious pathogens would be expected to induce distinctive profiles of immune responses in the liver , which might be manipulated to create specific and effective therapeutic strategies.
Several methods exist by which to determine the comprehensive transcriptomic profile of a pathogen-specific immune response, including microarray and quantitative real-time PCR [22, 23]. However, the high-throughput RNA sequencing (RNA-seq) technology offers several advantages over the other profiling applications. Not only is RNA-seq independent on predefined probes, which facilitates the discovery of new transcript variants, but the sequence platform also produces low background noise, which allow for distinction between closely homologous genes and detection of weakly expressed transcripts . In addition, concurrent advances in the bioinformatic algorithms used to analyze the RNA-seq data have allowed for better interpretation of the whole transcriptomic profile and provided further insights into complex molecular processes. The RNA-seq approach has already been successfully applied to several infectious disease models of zebrafish [25–27], including zebrafish embryo infected with Salmonella, and adult zebrafish and embryos infected by Mycobacteria[25, 28]. In addition, other fish species infection models have been subjected to RNA-seq analysis, including large yellow croaker (Larimichthys crocea) infected by Aeromonas hydrophila and Japanese seabass (Lateolabrax japonicus) infected by Vibrio Harveyi, but the overall immune-related transcription profiles have differed among species [25–30]. No reports exist in the literature of RNA-seq technology used to analyze the changes in an infected fish transcriptome profile induced upon vaccine treatment.
Edwardsiellosis, caused by the gram-negative Edwardsiella tarda, is currently one of the most economically disastrous infectious diseases affecting the global aquaculture industry . E. tarda displays polymorphic phenotypes and has a broad range of hosts from aquatic invertebrates to higher vertebrates, including birds, reptiles, mammals, and even humans . In developing a putative live attenuated vaccine against edwardsiellosis, Xiao et al. constructed an E. tarda mutant (WED) with low residual virulence. Although the mutant was capable of inducing robust protection in zebrafish and turbot, the antibody titers detected in sera were relatively low. By thoroughly understanding the immune-mechanism of zebrafish induced by the putative live attenuated vaccine, a more immunogenic vaccine may be able to be generated. To this end, we performed a comparative gene transcription analysis of livers from mock-immunized and WED-immunized zebrafish using RNA-seq technology to investigate their differential transcriptsomic profile. Furthermore, 12 genes associated with MHC antigen processing were analyzed by qPCR and the results revealed an activated MHC-I pathway and an inhibited MHC-II pathway during the early stage of vaccine immunization. It was prompted that WED conferred a robust protection in zebrafish by eliciting an effective cell immunity via the MHC-I pathway.
RNA-seq of liver transcriptome
To better understand the early stage immune response of zebrafish immunized with WED, six Solexa cDNA libraries were constructed from the livers of mock-immunized and WED-immunized zebrafish (Additional file 1). Biological replicates were pooled to make representative samples for deep sequencing analysis. Across the two groups of triplicate data, after normalization of the generated 95 bp PE raw reads, 15,683,828, 13,040,780 and 25,660,654 reads were obtained from C1-C3, and 16,306,312, 15,589,848 and 10,906,906 reads from V1-V3, respectively. To assess the quality of sequencing, the reads were mapped to the zebrafish reference genome. From the reads of each group, successful mapping occurred for 10,823,266 (C1), 9,584,828 (C2), 18,321,987 (C3), 12,209,418 (V1), 11,675,593 (V2) and 8,605,104 (V3) reads. However, 4,860,562, 3,455,952 and 7,338,667 unmapped reads were generated from C1-C3, while 4,096,894, 3,914,255 and 2,301,802 unmapped reads were found in V1-V3, respectively; we plan to conduct de novo analysis of these unmapped reads to generate a better reference of immune-relevant genes in zebrafish.
Analysis of differential expression among WED- and mock-immunized zebrafish liver
Statistically significant KEGG classifications of differentially expressed genes
Protein processing in endoplasmic reticulum
Arginine and proline metabolism
Primary bile acid biosynthesis
Fructose and mannose metabolism
Glycine, serine and threonine metabolism
In general, based on the results from GO analysis (by BiNGO plugin) and KEGG pathway analysis, the up- and down-regulated genes that were highly related to immune response of fish after WED immunization, significantly grouped into acute phase response (APR), complement activation, immune/defense response, and antigen processing and presentation pathway.
The acute phase response is conserved in zebrafish liver following WED immunization
Significantly differentially expressed genes in acute phase response
Transferrin receptor (p90, CD71)
Complement component 1, q subcomponent
Complement component 2
Complement factor B
Serum amyloid P component
Complement component C9
Complement factor H, like 2
Transferrin receptor 2
Fibrinogen alpha chain
Complement component c3a
PREDICTED: complement C3-H2-like
Plasminogen activator inhibitor type 1
CCAAT/enhancer binding protein, beta
Fibrinogen, gamma polypeptide
Coagulation factor II (Thrombin) receptor
ATPase, Cu++ transporting
Toll-like receptor 5
Complement component 7
Leukocyte cell-derived chemotaxin 2
Serum amyloid protein A
Traditionally, complement has been considered as a supportive first line of defense against microbial intruders . In WED-immunized zebrafish liver, three isoforms of complement C3, and the complement C4, C2, C7 and C9 were remarkably up-regulated. The C1 inhibitor and the C1q were up-regulated and down-regulated, respectively (Table 2).
Hepatocytes, which account for 80% of the liver mass , are the primary site of synthesis for all the genes involved in ion-binding and transporting [33, 35, 37]. In the RNA-seq data, more than 20 differentially expressed genes involved in ion-binding and transporting were strongly induced in zebrafish liver upon WED immunization (Table 2). These included haptoglobin, hemopexin, ceruloplasmin, transferrin receptor 2, ATPase, and Cu2+ transporting alpha polypeptide. Intelectin , which is involved in iron homeostasis, binding and transport, was one of the most up-regulated genes (129-fold) in the ion-binding and transporting category; However, the functions of intelectin in the contexts of normal iron metabolism and disease defense in zebrafish need to be further clarified. Members of the transferrin and ferritin families were significantly affected to result in obvious up- and down-regulation in zebrafish liver by WED immunization (Table 2). Leukocyte cell–derived chemotaxin 2 (LECT2) , originally named for its possible neutrophil chemotactic activity in vitro, was strongly induced by 222.8-fold in WED-immunized group, but its function in zebrafish remains unknown.
WED immunization induces defense responses and signaling transduction pathways
Immunity is a complex process of tightly controlled signals that involve a broad array of receptors, cytokines, enzymes, signal transducers, transcription factors, and other functional proteins. In our study, WED immunization increased dramatically the expression of cytokine genes related to the Jak-STAT, MAPK, TGF-β, apoptosis and VEGF signaling pathways. Therefore, WED-induced gene expression in zebrafish liver might facilitate protection against E. tarda by activating these pathways. Similar results were obtained in a previous study of large yellow croaker spleen during A. hydrophila infection . Since the majority of the differentially expressed genes in these signaling pathways were up-regulated in our study, it is possible that the WED immunization of zebrafish is capable of triggering a vigorous adaptive immune response.
WED immunization induces the antigen processing and presentation pathway
Differentially expressed genes related to the antigen processing and presentation pathways
MHC class II, transactivator
Beta chain, MHC class II
MHC class II DAB gene
Novel MHC II beta chain protein
Lysosomal membrane glycoprotein 2
MHC, class II invariant chain, CD74 molecule,
Cathepsin S, b.1
CD68 antigen variant
CD151 antigen, like
Proteasome activator complex subunit 3
Integral membrane protein 1(STT3)
Proteasome activator PA28 subunit 2
Proteasome maturation protein
TAP binding protein-like(TAPBPL)
Proteasome activator PA28 subunit 1
MHC class I ZE protein
T-cell immunomodulatory protein
Heat shock protein 4a
Proteasome (prosome, macropain) subunit, alpha
ER-resident chaperone calreticulin
Proteasome (prosome, macropain) subunit, beta
Proteasome (prosome, macropain) inhibitor
Proteasome (prosome, macropain)
Heat shock protein HSP 90 kDa alpha 2
Proteasome (macropain) 26 S subunit
Heat shock cognate 70 kDa protein
GATA-binding protein 2a
Cathepsin L, 1 a
DnaJ (Hsp40) homolog, subfamily B
Up-regulated genes with established roles in immune responses comprised another large functional category, indicating that active immune-surveillance, immune signaling, and immune cell activation were triggered in the WED-immunized zebrafish liver, like the MHC-I ZE protein (by 3.09-fold). However, the MHC class II DAB, MHC class II beta chain, MHC class II invariant chain (CD74 molecule), MHC class II transactivator (CIITA), and cathepsin S were down-regulated at this stage (Table 3).
qPCR analysis of MHC processing pathways
At present, molecular studies on the immune response to pathogens in fish models are mainly focused on infectious disease pathogenesis. RNA-seq and microarray-based transcriptome profiling studies have revealed that the teleosts are useful in vivo models for identifying host determinants of responses to bacterial infection [27–30]. Furthermore, the RNA-seq approach has already been successfully applied to several infectious disease models of zebrafish [25–27]. However, none have applied the RNA-seq technology to elucidate the immune-related pathways underlying the zebrafish response to vaccination for more effective vaccine evaluation. In this work, in order to gain comprehensive insight into the immunogenetics of zebrafish following immunization with the putative E. tarda live attenuated vaccine, a high-throughput deep sequencing-by-synthesis technology was used to investigate the immunization-related gene expression patterns. DESeq analysis identified 4565 significantly differentially expressed genes in the zebrafish liver following WED immunization. GO and KEGG analysis revealed that the genes involved in the ER protein processing as well as the phagosome and antigen processing and presentation pathways are regulated at the early stage following WED immunization (Table 1 and Figure 2). Significantly, two class MHC pathways were found to be reversely regulated upon immunization, and the MHC class I pathway was activated and the MHC class II pathway was inhibited (Figures 45 and 6). Both the RNA-seq results and qPCR data from our study of zebrafish liver during the early stage after WED immunization indicated that activation of the MHC-I processing pathway in teleosts could elicit cellular immune responses for protection.
Once bacterial vaccines are administrated into the animal host, they are often internalized by phagocytes via different entry mechanism. However, the subsequent issues involved in microbial sensing and antigen processing are not well defined. In the conventional paradigm, MHC class II molecules present antigenic fragments acquired by the endocytic route to the immune system for recognition and activation of CD4+ T cells . MHC class I molecules, on the other hand, are restricted to surveying the cytosol for endogenous antigen from intracellular pathogens (such as bacteria, parasites, and viruses), tumors, or self-proteins, which are degraded into proteasomal products and then presented on MHC class I molecules to CD8+ T cells, thus exersting an irreplaceable role on cellular-mediated immuno-protection toward intracellular pathogens [43, 44]. E. tarda is believed to be an intracellular pathogen that can survive and replicate within large phagosomes in macrophages . Since WED is an attenuated strain from wild type E. tarda, it could be assumed that WED bacteria possess the ability to survive in phagosomes of APC cells and the internalized bacteria are recognized as endogenous or exogenous antigen which would be presented or cross-presented by the MHC-I pathway, and finally evoking a CD8+ CTL-mediated response to achieve immune protection.
In MHC-I antigen processing pathway, antigenic peptides are degraded in the cytoplasm by proteasome, then translocated into the ER and loaded onto MHC-I molecules with the help of several protein components. PA28, as an important proteasome activator, is a heterohexameric ring that binds to one or both ends of the 20 S proteasome [43, 46]. Upon binding, it increases the catalytic activity of all three of the proteasome active sites, leads to changes in substrate cleavage, thereby generating more MHC class I-presented peptides [46, 47]. Khan et al. reported that constitutive proteasomes were replaced with immune-proteasomes in mice livers starting at two days after Listeria monocytogenes infection. Immuno-proteasomes support the generation of MHC class I epitopes and shape immune-dominance hierarchies of CD8+ T cells . In mice, this switch is marked by the up-regulation of proteasome activator PA28 subunits, which alter the fragmentation of polypeptides through the proteasome and are inducible by IFN-γ . The study of immune responses to E. ictaluri infection in blue catfish liver demonstrated that both the PA28α and PA28β were up-regulated . In the study described herein, the genes encoding PA28 subunit 1, PA28 subunit 2 and PA28 subunit 3 were all up-regulated in zebrafish liver, which suggested a shift toward MHC class I antigen processing occurred at the early stage after WED immunization.
Heat shock proteins (HSP) are a type of highly conserved and ubiquitously expressed proteins that play an essential role as molecular chaperones in protein folding and transport within the cell  and possess the ability to stimulate MHC class I antigen processing . HSP/peptide complexes are taken up by APC via specific receptors, whose signaling leads to MHC-I presentation of HSP-associated peptides and the induction of specific CD8+ cytotoxic T cells . The antigenic peptides chaperoned by HSPs are known to be more efficient, by orders of magnitude, than the free peptides for presentation by MHC-I [49, 50]. In our work, three heat shock proteins (heat shock cognate 70 kDa protein, heat shock protein 4a and heat shock protein 90 kDa alpha 2) were found to be up-regulated following WED immunization, and the activated HSPs suggested that the internalized WED bacteria were processed and loaded onto MHC class I molecules, ultimately initiating initiate the CTLs.
As cited above, MHC class I molecules present antigenic peptides on cell surface for recognition by CD8+ T cells . Like other glycoproteins, the folding and assembly of MHC class I molecules require interactions with a number of chaperone molecules in the ER, some of which are specific to MHC class I molecules . Among the known ER chaperones, endoplasmin (grp94) possesses the ability to bind peptides suitable for assembly on to MHC class I molecules together with calreticulin . Calreticulin and calnexin are specialized ER lectin-binding chaperones to bind transiently to newly-synthesized glycoproteins, but the calreticulin has been suggested as unique to interactions with the HSP/grp94 complex, which leads to recruitment of ER protein 57 . The interaction between calnexin and MHC class I molecules is believed to stabilize the class I heavy chain and help it to associate with the β2m component [51, 53]. In this work, the three ER chaperons, calreticulin, calnexin and endoplasmin (grp94), were all found to be induced in WED-immunized zebrafish liver, providing further evidence that an active MHC class I processing pathway was stimulated by WED immunization. In addition, TAP binding protein, another molecule involved in MHC class I antigen loading [44, 49, 51, 53], and MHC class I complex ZE protein were also up-regulated in WED-immunized zebrafish liver, strongly suggesting a vigorous activation of the MHC-I processing pathway.
The MHC antigen processing-associated genes from zebrafish have been extensively characterized. However, little is known about their expression patterns in zebrafish following vaccine immunization. Recently, the coordinated up-regulation of MHC class I-related components including MHC class I alpha chain, β2m, calreticulin, endoplasmin, PA28α and PA28β were reported in large yellow croaker following poly I:C injection  and in catfish following an intracellular bacterial infection . In this work, the RNA-seq data were given to show a coordinated down-regulation of several MHC class II antigen processing and presentation components, including the MHC-II DAB, MHC-II beta chain, MHC-II invariant chain (CD74), MHC class II transactivator (CIITA), cathepsin B and lysosomal membrane glycoprotein 2 (lamp2). This complex process is illustrated in Figure 4 and the differentially expressed genes are listed in Table 3. Furthermore, qPCR data confirmed the co-inhibition of lamp2, MHC-II dab, CD74, and CIITA in zebrafish liver and spleen (Figure 6). In previous researches, a remarkable inhibition of MHC-II expression and antigen presentation was ever reported in some pathogen infection models, including Brucella abortus, and Mycobacterium tuberculosis[56–58]. For pathogens, an ability to impair the antigen processing and presentation of host has been proposed to facilitate chronic infection by decreasing T cell responses to microbial antigens. For vaccines, however, the underlying significance of suppression of the MHC-II expression and antigen presentation remains unknown.
In conclusion, in this work, zebrafish was used as a model to investigate the host immune mechanisms underlying the protective effects of the E. tarda live attenuated vaccine. RNA-seq data revealed that the coordinate up-regulation of MHC-I processing pathways and down-regulation of MHC-II-associated pathways occurred at the early stage of vaccine immunization, providing insights into the molecular mechanisms of immune protection. The successful application of RNA-seq technology in the vaccine-zebrafish interaction model in this work established a new experimental platform for investigating the vaccine-specific host immune responses in a comprehensive and sensitive manner. Future studies using this approach will likely provide further significant insights into the detailed mechanisms of teleost immunity that will benefit the aquaculture industry, both from economic and human food source perspectives.
Fish and immunization
Healthy zebrafish (Danio rerio), weighing 0.3 ± 0.1 g and about 6 months of age, were obtained from the animal center at the East China University of Science and Technology (Shanghai, China) and maintained at 22 ± 2°C in a zebrafish cultivation system with a photo-period of 12:12 h (light : dark). Aquaria were supplied with flow-through dechlorinated and continuously aerated water at a rate of approximately 2×10-4 min-1. After at least one week of acclimatization, they were randomly divided into six treatment groups (70 fish per group) including three immunized groups (V1-V3) and three control groups (C1-C3), and the fish in each group were cultured in a separate tank. The fish in V1-V3 groups were intramuscularly (i.m.) injected with 1×105 CFU.fish-1 of WED bacteria in 5 μl phosphate buffered saline (PBS), as previously described , and the fish in C1-C3 groups were i.m. injected with 5 μl PBS alone. After two days of immunization, 20 fish from each of the three WED-immunized and three mock-immunized groups were sacrificed under anesthesia to obtain liver samples, and subsequently stored at −80°C until RNA extraction for RNA-seq analysis. Meanwhile, 10 fish from each group were sacrificed under anesthesia at days 1, 2, 3, and 5 post-immunization to obtain liver and spleen tissue samples, and subsequently stored at −80°C until RNA extraction for real-time qPCR analysis. All the zebrafish were handled in compliance with the local animal welfare regulations and maintained according to standard protocols (http://ZFIN.org). The immunization experiment was approved by the animal center at the East China University of Science and Technology (Shanghai, China).
Library preparation and sequencing
Total RNA was extracted from each tissue sample using the TRIzol reagent (Invitrogen, USA) according to the manufacturer’s instructions. To remove residual genomic DNA, the RNA samples were incubated with 10 units of DNA-free DNAse I (Ambion, USA) for 30 min at 37°C. The quality and quantity of the purified RNA were determined by measuring the absorbance at 260 nm/280 nm (A260/A280) using a Nanodrop ND-1000 spectrophotometer (LabTech, USA). RNA integrity was further verified by electrophoresis through a 1.5% (w/v) agarose gel.
Poly (A) mRNA was isolated from the total RNA samples with oligo (dT) magnetic beads (Invitrogen). The purified mRNA was fragmented by the RNA fragmentation kit (Ambion) and applied as template for first-strand cDNA synthesis using random hexamer primers and reverse transcriptase (Invitrogen). The second-strand cDNA was synthesized using RNase H (Invitrogen) and DNA polymerase I (New England Biolabs, USA). The Illumina Genomic DNA Sample Prep kit (Illumina, USA) was used to generate 120 bp paired-end (PE) cDNA libraries by following the manufacturer’s protocol. The libraries were loaded onto flow cell channels for sequencing on the Illumina HiSeq 2000 instrument by the Chinese National Human Genome Center (Shanghai, China). A total of six paired-end cDNA libraries of zebrafish livers were constructed for each of the test groups of WED-immunized and mock-immunized fish. Triplicate biological replicates were performed for each group. Raw data (tag sequences) were deposited in the NCBI database under submission number SRA048658.2.
The Illumina HiSeq 2000 system-generated 120 bp raw PE reads were first processed by the FASTX-Toolkit to remove the reads with sequencing adaptors and of low quality (phred quality <5). Then, the Burrows-Wheeler Aligner’s Smith-Waterman Alignment (BWA) program was used to align the remaining reads to the reference zebrafish mRNA from the Ensembl database . The transcription level of each gene was deduced by determining the total number of reads mapped to each gene using Picard tools (http://picard.sourceforge.net/). Differentially expressed genes were identified by the DESeq package in R software , using two-fold change (log2 (fold-change) ≥ 1 or ≤ -1) and p-value<0.05 (cut-off at 5% false discovery rate (FDR)) as the threshold. After data normalization by the p-value and FDR calculation, the resulting expression intensity values were analyzed by the MA plot-based method, as described by Wang et al..
Functional analysis of differentially expressed genes
The Database for Annotation, Visualization and Integrated Discovery (DAVID, v6.7)  was used to investigate functional enrichment for over- and under-expressed genes by more than two-fold in the WED-immunized group relative to the mock-immunized group. Gene functional enrichment was performed using the default parameters in DAVID to obtain an adjusted p-value <0.05 for the test gene group versus the zebrafish gene ontology (GO) annotation set. The fold-enrichment cut-off suggested for DAVID functional annotation is 1.5 . In addition, the significantly up-regulated genes from the differentially expressed genes dataset were further analyzed by investigating the corresponding GO biological processes. Furthermore, GO analysis of genes transcribed at different levels was also performed using the Biological Networks Gene Ontology (BiNGO) tool, which is based on the Cytoscape software (http://www.cytoscape.org). The hypergeometric test with Benjamini & Hochberg False Discovery Rate (FDR) was performed using the default parameters to obtain an adjusted p-value (<0.05) between the test gene group and the merged non-redundancy zebrafish (Danio rerio) and mouse (Mus musculus) GO annotation set. Finally, the web-based Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway analysis program run by the KEGG Automatic Annotation Server (KAAS) (http://www.genome.jp/tools/kaas/) was used to obtain functional annotation of genes by performing basic local alignment search tool (BLAST)-mediated comparisons against the manually-curated KEGG GENES database . We merged the most current KEGG GENES entries for Danio rerio and Mus musculus to generate a reference dataset and used the bi-directional best hit (BBH) information method to further analyze the significantly differentially expressed genes to gain insights into the related biological pathways.
To verify the differential expression detected by sequencing, qPCR was performed using the ABI Prism 7500 Detection System (Applied Biosystems, USA) with SYBR Green (Roche, USA) as the fluorescent detection dye, according to the manufacturer’s protocol. First-strand cDNA was synthesized from 1 μg of total mRNA, as described above, and applied as a template for qPCR with gene-specific primers (Additional file 4). Primers were designed using Primer Express 3 software. To determine the PCR efficiency, we first generated a standard curve by amplifying ten-fold serial dilutions of cDNA using primers to both the gene of interest and an internal control (β-actin), and all primers were optimized until PCR efficiency values fell in 1.80-2.15. The qPCR thermal cycling conditions for all reactions were 95°C for 15 min, followed by 40 cycles of 95°C for 5 s, 60°C for 20 s, and 72°C for 20 s. All qPCR reactions were performed for three biological replicates, and the data for each sample were expressed relative to the expression levels of β-actin by using the 2-ΔΔCT method . Independent-sample t-test in the SPSS software (Version 11.5, SPSS Inc.) was used to determine statistical significance. Significant differences were considered at p<0.01.
This work was supported by grants from the National High Technology Research and Development Program of China (No. 2008AA092501), the earmarked fund for China Agro-industry Research System (No. CARS-50), and the Science Foundation of Ministry of Education of China (No. Ncet-09-0344).
- Felix E, Graham JL: Zebrafish as a model for vertebrate hematopoiesis. Curr Opin Pharmacol. 2010, 10: 563-570. 10.1016/j.coph.2010.05.004.View ArticleGoogle Scholar
- Lieschke GJ, Currie PD: Animal models of human disease: zebrafish swim into view. Nat Rev Genet. 2007, 8: 353-367. 10.1038/nrg2091.View ArticlePubMedGoogle Scholar
- Meeker N, Trede N: Immunology and zebrafish: spawning new models of human disease. Dev Comp Immunol. 2008, 32: 745-757. 10.1016/j.dci.2007.11.011.View ArticlePubMedGoogle Scholar
- Sullivan C, Kim CH: Zebrafish as a model for infectious disease and immune function. Fish Shellfish Immunol. 2008, 25: 341-350. 10.1016/j.fsi.2008.05.005.View ArticlePubMedGoogle Scholar
- Chao CC, Hsu PC, Jen CF, Chen IH, Wang CH, Chan HC, Tsai PW, Tung KC, Wang CH, Lan CY, Chuang YJ: Zebrafish as a model host for candida albicans infection. Infect Immun. 2010, 78: 2512-2521. 10.1128/IAI.01293-09.PubMed CentralView ArticlePubMedGoogle Scholar
- Vojtech LN, George ES, Carla C, Ostland V, Hansen JD: Host immune response and acute disease in a zebrafish model of francisella pathogenesis. Infect Immun. 2009, 77: 914-925. 10.1128/IAI.01201-08.PubMed CentralView ArticlePubMedGoogle Scholar
- Van der Sar AM, Appelmelk BJ, Vandenbroucke-Grauls CMJE, Bitter W: A star with stripes: zebrafish as an infection model. Trends Microbiol. 2004, 12: 451-457. 10.1016/j.tim.2004.08.001.View ArticlePubMedGoogle Scholar
- Kanther M, Rawls JF: Host-microbe interactions in the developing zebrafish. Curr Opin Immunol. 2010, 22: 10-19. 10.1016/j.coi.2010.01.006.PubMed CentralView ArticlePubMedGoogle Scholar
- Tobin DM, May RC, Wheeler RT: Zebrafish: A see-through host and a fluorescent toolbox to probe host–pathogen interaction. PLoS Pathog. 2012, 8 (1): e1002349-10.1371/journal.ppat.1002349.PubMed CentralView ArticlePubMedGoogle Scholar
- Novoa B, Romero A, Mulero V, Rodríguez I, Fernández I, Figueras A: Zebrafish (Danio rerio) as a model for the study of vaccination against viral haemorrhagic septicemia virus (VHSV). Vaccine. 2006, 24: 5806-5816. 10.1016/j.vaccine.2006.05.015.View ArticlePubMedGoogle Scholar
- Cui Z, Samuel-Shaker D, Watral V, Kent ML: Attenuated mycobacterium marinum protects zebrafish against mycobacteriosis. J Fish Dis. 2010, 33: 371-375. 10.1111/j.1365-2761.2009.01115.x.PubMed CentralView ArticlePubMedGoogle Scholar
- Xiao JF, Chen T, Wang QY, Liu Q, Wang X, Lv YZ, Wu HZ, Zhang YX: Search for live attenuated vaccine candidate against edwardsiellosis by mutating virulence-related genes of fish pathogen Edwardsiella tarda. Lett Appl Microbiol. 2011, 53: 430-437. 10.1111/j.1472-765X.2011.03126.x.View ArticlePubMedGoogle Scholar
- Xu XP, Zhang LC, Weng SP, Huang ZJ, Lu J, Lan DM, Zhong XJ, Yu XQ, Xu AL, He JG: A zebrafish (Danio rerio) model of infectious spleen and kidney necrosis virus (ISKNV) infection. Virology. 2008, 376: 1-12. 10.1016/j.virol.2007.12.026.View ArticlePubMedGoogle Scholar
- Crispe IN: The liver as a lymphoid organ. Annu Rev Immunol. 2009, 27: 147-163. 10.1146/annurev.immunol.021908.132629.View ArticlePubMedGoogle Scholar
- Hofmann J, Greter M, Du Pasquier L, Becher B: B-cells need a proper house, whereas T-cells are happy in a cave: the dependence of lymphocytes on secondary lymphoid tissues during evolution. Trends Immunol. 2010, 31: 144-153. 10.1016/j.it.2010.01.003.View ArticlePubMedGoogle Scholar
- Dancygier H: The liver as an immune organ. Clin Hepatol. 2010, 11: 141-152.Google Scholar
- Wittamer V, Bertrand JY, Gutschow PW, Traver D: Characterization of the mononuclear phagocyte system in zebrafish. Blood. 2011, 117: 7126-7135. 10.1182/blood-2010-11-321448.View ArticlePubMedGoogle Scholar
- Plitas G, Burt BM, Stableford JA, Nguyen HM, Welles AP, DeMatteo RP: Dendritic cells are required for effective cross-presentation in the murine liver. Hepatology. 2008, 47: 1343-1351.View ArticlePubMedGoogle Scholar
- Warren A, Le Couteur DG, Fraser R, Bowen DG, McCaughan GW, Bertolino P: T lymphocytes interact with hepatocytes through fenestrations in murine liver sinusoidal endothelial cells. Hepatology. 2006, 44: 1182-1190. 10.1002/hep.21378.View ArticlePubMedGoogle Scholar
- Ebrahimkhani MR, Mohar I, Crispe IN: Cross-presentation of antigen by diverse subsets of liver cells. Hepatology. 2011, 54: 1379-1387. 10.1002/hep.24508.PubMed CentralView ArticlePubMedGoogle Scholar
- Gao B, Jeong W, Tian ZG: Liver: an organ with predominant innate immunity. Hepatology. 2008, 47: 729-736.View ArticlePubMedGoogle Scholar
- Morozova O, Hirst M, Marra MA: Applications of new sequencing technologies for transcriptome analysis. Annu Rev Genomics Hum Genet. 2009, 10: 135-151. 10.1146/annurev-genom-082908-145957.View ArticlePubMedGoogle Scholar
- Han X, Wu X, Chung WY, Li T, Nekrutenko A, Altman NS, Chen G, Ma H: Transcriptome of embryonic and neonatal mouse cortex by highthroughput RNA sequencing. Proc Natl Acad Sci USA. 2009, 106: 12741-12746. 10.1073/pnas.0902417106.PubMed CentralView ArticlePubMedGoogle Scholar
- Tang F, Barbacioru C, Wang Y, Nordman E, Lee C, Xu N, Wang X, Bodeau J, Tuch BB, Siddiqui A, Lao K, Surani MA: mRNA-Seq whole-transcriptome analysis of a single cell. Nat Methods. 2009, 6: 377-382. 10.1038/nmeth.1315.View ArticlePubMedGoogle Scholar
- Hegedus Z, Zakrzewskaa A, Ágostonb VC, Ordas A, Rácz P, Mink M, Spaink HP, Meijer AH: Deep sequencing of the zebrafish transcriptome response to mycobacterium infection. Mol Immunol. 2009, 46: 2918-2930. 10.1016/j.molimm.2009.07.002.View ArticlePubMedGoogle Scholar
- Ordas A, Hegedus Z, Henkel CV, Stockhammer OW, Butler D, Jansen HJ, Racz P, Mink M, Spaink HP, Meijer AH: Deep sequencing of the innate immune transcriptomic response of zebrafish embryos to Salmonella infection. Fish Shellfish Immunol. 2011, 31: 716-724. 10.1016/j.fsi.2010.08.022.View ArticlePubMedGoogle Scholar
- Liselotte V, Hong J, Per U, Outi H, Juha K: The zebrafish transcriptome during early development. BMC Dev Biol. 2011, 11: 30-10.1186/1471-213X-11-30.View ArticleGoogle Scholar
- Van der Sar AH, Spaink HP, Zakrzewska A, Bitter W, Meijer AH: Specificity of the zebrafish host transcriptome response to acute and chronic mycobacterial infection and the role of innate and adaptive immune components. Mol Immunol. 2009, 46: 2317-2332. 10.1016/j.molimm.2009.03.024.View ArticlePubMedGoogle Scholar
- Mu YN, Ding F, Cui P, Ao JG, Hu SN, Chen XH: Transcriptome and expression profiling analysis revealed changes of multiple signaling pathways involved in immunity in the large yellow croaker during Aeromonas hydrophila infection. BMC Genomics. 2010, 11: 506-10.1186/1471-2164-11-506.PubMed CentralView ArticlePubMedGoogle Scholar
- Xiang LX, He D, Dong WR, Zhang YW, Shao JZ: Deep sequencing-based transcriptome profiling analysis of bacteria-challenged Lateolabrax japonicus reveals insight into the immune-relevant genes in marine fish. BMC Genomics. 2010, 11: 472-10.1186/1471-2164-11-472.PubMed CentralView ArticlePubMedGoogle Scholar
- Mohanty BR, Sahoo PK: Edwardsiellosis in fish: a brief review. J Biosci. 2007, 32: 1331-1344. 10.1007/s12038-007-0143-8.View ArticlePubMedGoogle Scholar
- Wang LK, Feng ZX, Wang X, Wang XW, Zhang XG: DEGseq: an R package for identifying differentially expressed genes from RNA-seq data. Bioinformatics. 2010, 26: 136-138. 10.1093/bioinformatics/btp612.View ArticlePubMedGoogle Scholar
- Gerwick L, Corley-Smith G, Bayne CJ: Gene transcript changes in individual rainbow trout livers following an inflammatory stimulus. Fish Shellfish Immunol. 2007, 22: 157-171. 10.1016/j.fsi.2006.04.003.View ArticlePubMedGoogle Scholar
- Peatman E, Terhune J, Baoprasertkul P, Xu P, Nandi S, Wang S, Somridhivej B, Kucuktas H, Li P, Dunham R, Liu Z: Microarray analysis of gene expression in the blue catfish liver reveals early activation of the MHC class I pathway after infection with Edwardsiella ictaluri. Mol Immunol. 2008, 45: 553-566. 10.1016/j.molimm.2007.05.012.View ArticlePubMedGoogle Scholar
- Lin B, Chen SW, Cao Z, Lin YQ, Mo DZ, Zhang HB, Gu JD, Dong ML, Liu ZH, Xu AL: Acute phase response in zebrafish upon Aeromonas salmonicida and Staphylococcus aureus infection: Striking similarities and obvious differences with mammals. Mol Immunol. 2007, 44: 295-301. 10.1016/j.molimm.2006.03.001.View ArticlePubMedGoogle Scholar
- Nakao M, Tsujikura M, Ichiki S, Vo TK, Somamoto T: The complement system in teleost fish: progress of post-homolog-hunting researches. Dev Comp Immunol. 2011, 35: 1296-1308. 10.1016/j.dci.2011.03.003.View ArticlePubMedGoogle Scholar
- John BV, Sharifa L: The binding and transport of alternative metals by transferring. Biochimica et Biophysica Acta. 2012, 1820: 362-378. 10.1016/j.bbagen.2011.07.003.View ArticleGoogle Scholar
- Stockhammer OW, Zakrzewska A, Hegeduˆs Z, Spaink HP, Meijer AH: Transcriptome profiling and functional analyses of the zebrafish embryonic innate immune response to Salmonella infection. J Immunol. 2009, 182: 5641-5653. 10.4049/jimmunol.0900082.View ArticlePubMedGoogle Scholar
- Sepulcre MP, Francisca AP, Azucena LM, Roca FJ, Meseguer J, Cayuela ML, Mulero V: Evolution of lipopolysaccharide (LPS) recognition and signaling: fish TLR4 does not recognize LPS and negatively regulates NF-κB activation. J Immunol. 2009, 182: 1836-1845. 10.4049/jimmunol.0801755.View ArticlePubMedGoogle Scholar
- Seki E, Brenner DA: Toll-like receptors and adaptor molecules in liver disease: update. Hepatology. 2008, 48: 322-335. 10.1002/hep.22306.View ArticlePubMedGoogle Scholar
- Vyas JM, Van der Veen AG, Ploegh HL: The known unknowns of antigen processing and presentation. Nat Rev Immunol. 2008, 8: 607-618. 10.1038/nri2368.PubMed CentralView ArticlePubMedGoogle Scholar
- Van den Hoorn T, Paul P, Jongsma ML, Neefjes J: Routes to manipulate MHC class II antigen presentation. Curr Opin Immunol. 2011, 23: 88-95. 10.1016/j.coi.2010.11.002.View ArticlePubMedGoogle Scholar
- Sijts EJ, Kloetzel PM: The role of the proteasome in the generation of MHC class I ligands and immune responses. Cell Mol Life Sci. 2011, 68: 1491-1502. 10.1007/s00018-011-0657-y.PubMed CentralView ArticlePubMedGoogle Scholar
- Ramachandra L, Simmons D, Harding CV: MHC molecules and microbial antigen processing in phagosomes. Curr Opin Immunol. 2009, 21: 98-104. 10.1016/j.coi.2009.01.001.PubMed CentralView ArticlePubMedGoogle Scholar
- Okuda J, Arikawa Y, Takeuchi Y: Intracellular replication of Edwardsiella tarda in murine macrophage is dependent on the type III secretion system and induces an up-regulation of anti-apoptotic NF-kB target genes protecting the macrophage from staurosporine-induced apoptosis. Microb Pathog. 2006, 41: 226-240. 10.1016/j.micpath.2006.08.002.View ArticlePubMedGoogle Scholar
- Rock KL, Farfán-Arribas DJ, Shen L: Proteases in MHC class I presentation and cross-presentation. J Immunol. 2010, 184: 9-15. 10.4049/jimmunol.0903399.PubMed CentralView ArticlePubMedGoogle Scholar
- Yamano T, Sugahara H, Mizukami S, Murata S, Chiba T, Tanaka K, Yui K, Udono H: Allele-selective effect of PA28 in MHC class I antigen processing. J Immunol. 2008, 181: 1655-1664.View ArticlePubMedGoogle Scholar
- Khan S, van den Broek M, Schwarz K, de Giuli R, Diener PA, Groettrup M: Immunoproteasomes largely replace constitutive proteasomes during an antiviral and antibacterial immune response in the liver. J Immunol. 2001, 167: 6859-6868.View ArticlePubMedGoogle Scholar
- Hartl FU: Molecular chaperones in cellular protein folding. Nature. 1996, 381: 571-579. 10.1038/381571a0.View ArticlePubMedGoogle Scholar
- Osterloh A, Breloer M: Heat shock proteins: linking danger and pathogen recognition. Med Microbiol Immunol. 2008, 197: 1-8.View ArticlePubMedGoogle Scholar
- Nicchitta CV, Reed RC: The immunological properties of endoplasmic reticulum chaperones: a conflict of interest?. Essays Biochem. 2000, 36: 15-25.View ArticlePubMedGoogle Scholar
- Ireland BS, Brockmeier U, Howe CM, Elliott T, Williams DB: Lectin-deficient calreticulin retains full functionality as a chaperone for class I histocompatibility molecules. Mol Biol Cell. 2008, 19: 2413-2423. 10.1091/mbc.E07-10-1055.PubMed CentralView ArticlePubMedGoogle Scholar
- Crotzer VL, Blum JS: Autophagy and its role in MHC-mediated antigen presentation. J Immunol. 2009, 182: 3335-3341. 10.4049/jimmunol.0803458.PubMed CentralView ArticlePubMedGoogle Scholar
- Liu G, Zheng W, Chen X: Molecular cloning of proteasome activator PA28-beta subunit of large yellow croaker (Pseudosciana crocea) and its coordinated up-regulation with MHC class I alpha-chain and beta(2)- microglobulin in poly I:C-treated fish. Mol Immunol. 2007, 44: 1190-1197. 10.1016/j.molimm.2006.06.024.View ArticlePubMedGoogle Scholar
- Barrionuevo P, Cassataro J, Delpino MV: Brucella abortus inhibits major histocompatibility complex class II expression and antigen processing through interleukin-6 secretion via toll-like receptor 2. Infect Immun. 2008, 76: 250-262. 10.1128/IAI.00949-07.PubMed CentralView ArticlePubMedGoogle Scholar
- Torres M, Ramachandra L, Rojas RE, Bobadilla K, Thomas J, Canaday DH, Harding CV, Boom WH: Role of phagosomes and major histocompatibility complex class II (MHC-II) compartment in MHC-II antigen processing of Mycobacterium tuberculosis in human macrophages. Infect Immun. 2006, 74: 1621-1630. 10.1128/IAI.74.3.1621-1630.2006.PubMed CentralView ArticlePubMedGoogle Scholar
- Reuter MA, Pecora ND, Harding CV, Canaday DH, McDonald D: Mycobacterium tuberculosis promotes HIV trans-infection and suppresses major histocompatibility complex class II antigen processing by dendritic cells. J Virol. 2010, 84: 8549-8560. 10.1128/JVI.02303-09.PubMed CentralView ArticlePubMedGoogle Scholar
- Noss EH, Harding CV, Boom WH: Mycobacterium tuberculosis inhibits MHC Class II antigen processing in murine bone marrow macrophages. Cell Immunol. 2000, 201: 63-74. 10.1006/cimm.2000.1633.View ArticlePubMedGoogle Scholar
- Li H, Durbin R: Fast and accurate long-read alignment with Burrows-Wheeler transform. Bioinformatics. 2010, 26: 589-595. 10.1093/bioinformatics/btp698.PubMed CentralView ArticlePubMedGoogle Scholar
- Anders S, Huber W: Differential expression analysis for sequence count data. Genome Biol. 2010, 11: R106-10.1186/gb-2010-11-10-r106.PubMed CentralView ArticlePubMedGoogle Scholar
- Huang DW, Sherman BT, Lempicki RA: Systematic and integrative analysis of large gene lists using DAVID bioinformatics resources. Nat Protocols. 2009, 4: 44-57.View ArticleGoogle Scholar
- Moriya Y, Itoh M, Okuda S, Yoshizawa AC, Kanehisa M: KAAS: an automatic genome annotation and pathway reconstruction server. Nucl Acids Res. 2007, 35: 182-185. 10.1093/nar/gkm321.View ArticleGoogle Scholar
- Thomas DS, Livak KJ: Analyzing real-time PCR data by the comparative CT method. Nat Protocols. 2008, 3: 1101-1108. 10.1038/nprot.2008.73.View ArticleGoogle Scholar
This article is published under license to BioMed Central Ltd. This is an Open Access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/2.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.