- Research article
- Open Access
RNA-Seq analysis of differentially expressed genes relevant to innate and adaptive immunity in cecropin P1 transgenic rainbow trout (Oncorhynchus mykiss)
BMC Genomics volume 19, Article number: 760 (2018)
In the past years, our laboratory successfully generated transgenic rainbow trout bearing cecropin P1 transgene. These fish exhibited resistant characteristic to infection by Aeromonas salmonicida, Infectious Hematopoietic Necrosis Virus (IHNV) and Ceratomyxa shasta (a parasitic pathogen). Previously, treating rainbow trout macrophage cells (RTS-11) with cecropin B, pleurocidin and CF17, respectively, resulted in elevated expression of two pro-inflammatory genes, e.g. cyclooxygenase-2 (cox-2) and interleukin-1β (il-1β). In addition, a profiling of global gene expression by 44 k salmonid microarray analysis was conducted, and the results showed that immune relevant processes have been perturbed in cecopin P1 transgenic rainbow trout. Therefore, we hypothesized that cecropin P1 may not only eliminate pathogens directly, but also modulate the host immune systems, leading to increased resistance against pathogen infections. To confirm this hypothesis, we performed de novo mRNA deep sequencing (RNA-Seq) to analyze the transcriptomic expression profiles in three immune competent tissues of cecropin P1 transgenic rainbow trout.
De novo sequencing of mRNA of the rainbow trout spleen, liver and kidney tissues were conducted by second-generation Illumina system, followed by Trinity assembly. Tissue specific unigenes were obtained, and annotated according to the Gene Ontology (GO) and the Nucleotide Basic Local Alignment Search Tool (BLAST). Over 2000 differentially expressed genes (DEGs) were determined by normalized ratio of Reads Per Kilobase of transcript per million mapped reads (RPKM) among the transgenic and non-transgenic fish in a tissue specific manner, and there were 82 DEGs in common among the three tissues. In addition, the enrichment analysis according to Gene Ontology Biological Process (GO:BP), and Kyoto Encyclopedia of Genes and Genomes (KEGG) based pathway analysis associated with innate/adaptive immunity of fish were also performed to illustrate the altered immune-related functions in each tissue.
According to the RNA-Seq data, the correlations between alteration of gene expression profiles and the functional perturbations of the host immune processes were revealed. In comparison with the results of cDNA microarray analysis conducted by Lo et al., the overall results supported our hypothesis that the gene product of cecropin P1 transgene may not only directly eliminate pathogens, but also modulate the host immune system. Results of this study present valuable genetic information for Oncorhynchus mykiss, and will benefit future studies on the immunology of this fish species.
Outbreak of diseases caused by pathogen infection is one of the most serious bottlenecks in the aquaculture industry worldwide. Traditional strategies including vaccination, treatment of diseased fish with antibiotics, and artificial selection of disease resistant fish strains by traditional approach have been used to control fish diseases. However, deficiencies, such as high economic costs of vaccination, selection of antibiotic-resistant pathogens, lacking effective cure of viral infection, and low degree of overall protection of fish population associated with these traditional strategies have limited their feasibilities. Therefore, effective means of controlling fish diseases in the aquaculture industry are awaiting to be developed, and transgenic fish technology may provide a solution to this problem. Mediated via the transgenic fish technology, directly modifying the unwanted genetic traits that cause the vulnerability of fish to infection by pathogens, or introducing specific genes that may confer resistance to pathogen infection into the fish genome might achieve the purpose of protecting fish from invasion by pathogens (see review by Chen et al., 2014) .
Cecropin B, which was first identified in cecropia moth, Hyalophora cecropia, is one of the antimicrobial peptide (AMP) family member proteins playing an essential role in the innate immunity of insects . According to Shai , the zwitterionic AMPs, 31 to 37 amino acid residues in length and α-helical cationic amphipathic peptides, can be integrated into the cellular membrane of pathogens and resulted in the formation of pores on the cellular membrane leading to ultimate lysis of pathogen cells . Cecropin P1, identified from nematode inhabiting in the porcine small intestine, was found to be more potent against Gram-negative bacteria than Gram-positive bacteria , and was thus chosen as transgenic target gene in our laboratory for production of transgenic rainbow trout.
In the past years, transgenic medaka  and rainbow trout  that harbor cecropin B or cecropin P1 transgene were generated by microinjecting the transgene into fertilized eggs or by electroporating the transgene into sperm. Through repeated challenge studies, families of heterozygous cecropin P1 transgenic rainbow trout in the second and third generations were shown to be resistant to infection by Aeromonas salmonicida, infectious hematopoietic necrosis virus (IHNV) and Ceratomyxa shasta (a common parasite infecting rainbow trout) [7, 8]. In addition, Chiou et al. demonstrated that in vitro treatment of rainbow trout RTS11 cells (a trout macrophage cell line) with cecropin B and its synthetic analogue (CF-17) led to up-regulation of interleukin-1β (il-1β), which is a pro-inflammatory cytokine, and cyclooxygenase (cox-2), an enzyme essential for the inflammatory modulation . Although AMPs can eliminate pathogens directly, Chiou et al.  and Lai  further suggested that AMPs might also initiate an immune-modulatory process in the host immune system. To support this perception, a preliminary study on the analysis of gene expression profiles of immune related genes in two families of cecropin P1 transgenic rainbow trout was conducted by Lo et al.  in our laboratory via cDNA microarray analysis on a 44 k custom made salmonid chip ; many tissue-specific differentially expressed genes (DEGs), namely in the spleen, liver, and kidney, were determined. Organ specific functional perturbations of the host innate/adaptive immune pathways were identified including phagocytosis, lysosomal processing, complement activation, antigen presenting, and leukocyte migration. Furthermore, disturbance of biological processes, which may contribute indirectly to host immunity, were also determined such as lipid metabolic process, cellular focal adhesion, and extracellular matrix (ECM)-organization . By combining these facts, we hypothesized that transgenic cecropin P1 may not only eliminate pathogens directly, but also modulate the host’s innate and adaptive immunity. However, cDNA microarray conducted by Lo et al. has several limitations, e.g. (1) hybridization based florescence detection limits the measurement accuracy of the expression levels, especially the transcripts existing in very high or low copy numbers ; (2) only 44 k genes were analyzed and didn’t provide a global view of the total transcriptome; (3) the gene chip was customized according to Salmo salar genome and considered not thoroughly reflecting the trout genomic information. To overcome these limitations, and to extend the study for more solid evidence, a cross-platform confirmation is highly desirable.
The recent advances of the next generation sequencing technology (NGS), i.e., deep RNA sequencing (RNA-Seq), has exerted tremendous impact on studies on “de novo” construction of the transcriptome without a reference genome, evaluation of nucleotide variations, and evaluation of methylation patterns of genes etc. [15,16,17]. This technology has advantages over the cDNA microarray analysis in the following aspects: firstly, high levels of data reproducibility leading to reduction of technical replications for the experiments; secondly, allowing easy identification and quantification of the expression of isoforms and unknown transcripts; thirdly, increasing the popularity of high throughput sequencing technologies resulted in a significant reduction of the cost the RNA-Seq experiments. To confirm the mentioned hypothesis in the current study, de novo mRNA deep sequencing by the Illumina second generation system was carried out in three different immune relevant tissues (namely the spleen, liver, and kidney) of cecropin P1 transgenic rainbow trout, and followed by Trinity assembly of the data . By sorting of Reads Per Kilobase of transcript per million mapped reads (RPKM), the DEGs were determined. From the GeneCodis enrichment analysis [19, 20], we discovered functional alterations of biological processes with identifiers of Gene Ontology: Biological Process (GO:BP) and Kyoto Encyclopedia of Genes and Genomes (KEGG) databases. Finally, by establishing KEGG-based pathways analyses, we further revealed that immune relevant processes in the spleen and kidney, and energy metabolism relevant processes in the liver were significantly perturbed. These results strongly support our hypothesis, and will benefit future studies on the genetics of fish immunology. Furthermore, these perturbed DEGs may also serve as biological markers for artificial selection and breeding of disease resistant aquaculture important fish species.
De novo sequencing, reads assembly, and annotation of reference database from non-transgenic fish
Pools of RNA samples from three tissues (i.e., the spleen, liver, and kidney) of non-transgenic rainbow trout were used to generate the reference genome library database. By removing reads of adaptor, reads of unknown nucleotides greater than 5% and reads with low quality bases (base quality ≦ 10), 103,048,352 paired-end clean reads from total nucleotides of 9,274,351,680 were obtained. After de novo assembly by the Trinity program , 257,490 contigs with the average length of 285 bp were collected, and a total of 141,850 unigenes with the average length of 523 bp were determined (Table 1). Subsequently, the unigenes were annotated by BLAST-p according to six databases, namely nucleotide collection NR/NT, Swiss-Prot, KEGG, GO, and Clusters of Orthologous (COG). The annotated unigenes were than used as reference genome library for further transgenic transcriptomic sequence mapping and assembly.
Sequencing and assembly of cecropin P1 transgenic rainbow trout
Messenger RNA isolated from three tissues (the spleen, liver, and kidney) of cecropin P1 transgenic rainbow trout were individually sequenced by Illumina HiSeq 2000. By eliminating reads of adaptors, reads of unknown nucleotides larger than 10% and reads with low quality bases (base quality ≦ 5) greater than 50%, clean reads were obtained in a tissues specific manner. Afterwards, the clean reads were mapped to reference database generated from non-transgenic control by the aligning program of short oligonucleotide analysis package: SOAPaligner/SOAP2  with no more than two mismatches allowed to generate tissues specific unigenes. A sum of 131,671 unigenes was determined, and the overall good qualities of reads and aligning results were summarized in the Additional file 1.
Determination of differential expression profiles
By calculating the ratio of RPKM of unigenes  among transgenic and non-transgenic fish, the differentially expressed gene profiles were determined with a threshold of two folds been set. The genes with RPKM ratio greater than two folds (transgenic/non-transgenic > 2 or < 0.5) were considered as significantly up- or down-regulated, respectively, and were defined as differentially expressed genes (DEGs). As shown in Fig. 1, there were 1873/1402 genes in the spleen, 431/622 genes in the liver, and 916/1323 genes in the kidney were determined as significantly up−/down-regulated. After evaluating RPKM ratio in base two logarithm and eliminating redundancy by probes annotated with identical geneID, a total of 3275, 2239 and 1053 DEGs were identified in the spleen, kidney and liver, respectively. Additionally, by sorting RPKM ratio, a sum of 82 DEGs are in common among three tissues, 375 DEGs between the spleen and kidney, 86 DEGs between the spleen and liver, and 108 DEGs between the kidney and liver were identified (Fig. 2a). Moreover, these 82 DEGs belong to immune related biological functions, showing consistent or inverse expression patterns among the three tissues (Fig. 2c), where color gradient denoting relative fold changes of RPKM ratio. To confirm the gene expression level obtained by RNA-Seq, selected genes covering a wide range of expression ratios were determined by real-time quantitative reverse transcription poly chain reaction (RT-qPCR) analysis. By comparing the two independent platforms, a general agreement with an acceptable degree of linear correlation (R2 = 0.81, Fig. 2b) was observed, which indicating the expression dataset was solid.
Enrichment analysis of biology terms of GO and KEGG pathway by GeneCodis
To uncover the thematic association of altered gene expression patterns in cecropin P1 transgenic rainbow trout, enrichment analysis by GeneCodis was conducted. GeneCodis was developed as an over-representation analysis (ORA) approach [19, 20] and could be performed in either singular or modular enriching biological terms via different databases. The DEGs of each tissue and the total genes from reference genome were inputted to inquire GO biological process, GOSlim process and KEGG pathways. In Fig. 3, we reported the resulting tag clouds of modular GeneCodis enrichment analyses. The modular analyses, tissue specifically, indicated the intersections of sets of terms of biological processes among all and each of the three different annotation databases. Therefore, it revealed the most conspicuous terms with distinct profiles of enrichment. Afterwards, the statistical significance of terms in each tissue were ranked by hypergeometric and chi-square methods to discover the highly perturbed biological processes in each tissue studied. As shown in Table 2, the top 5 most perturbed biological functions in the spleen of the transgenic fish are cell adhesion, negative regulation of cell proliferation, innate and total immune response, and cytokine-mediated signaling pathways. In the kidney, functional biological processes found greatly altered include peroxisome proliferator-activated receptor (PPAR) signaling, cell adhesion, hematopoietic cell lineage, regulation of immune response, and inflammatory response. In the liver, fewer of the biological terms were enriched such as lipid metabolic process, immune response, and chemokine signaling pathways. As expected, energy efficiency correlated processes were also impacted in the liver of transgenic fish including PPAR signaling, starch and sucrose metabolism, and cholesterol metabolic and homeostasis. To discover a better broad-spectrum of tissue-specific direct and indirect effects on immune relevant DEGs, singular enrichment results were also analyzed (Additional file 2) since singular enrichment analyses enlighten the union of all sets of biological terms in each annotation database. The ECM-receptor interaction and hematopoietic cell lineage were altered in the kidney (Additional file 2, KEGG terms). Finally, starch and sucrose metabolism, chemokine signaling, and cytokine-cytokine receptor interaction were influenced in the liver as well (KEGG terms in the Additional file 2).
Analysis of functional perturbed immune relevant pathways
The alteration of biological processes in the spleen, liver and kidney of the transgenic rainbow trout were studied by GeneCodis enrichment analysis. To characterize genes contributing to immune relevant pathways that were affected by cecropin P1 transgene product, and to determine the signaling pathways that has been impacted, the pre-defined unigenes and DEGs from three tissues were subjected to analysis against KEGG pathways in order to assign their functions within the tissue-specific biological processes. By inputting a sum of 37,632 unigenes in all three tissues (the spleen, liver and kidney), 1286 DEGs in the spleen, 465 DEGs in the liver, and 936 DEGs in the kidney into the KEGG pathway database, 239, 212 and 235 KEGG pathways were mapped in the spleen, liver and kidney, respectively. By sorting the ratios among DEGs and all unigenes with pathway annotation (p-value < 0.05), followed by combining to the tissue-specific functions, the significantly perturbed immune relevant pathways were summarized in Table 3. Via ranking the ratio in the order of statistical significance (Table 3, I/R ratio), it reveals that out of 735 total annotated unigenes, a sum of 51 DEGs associate to leukocyte trans-endothelial migration, 45/562 to cytokine interaction, 32/493 to Janus kinase (JAK) to signal transducer and activator of transcription (STAT) signaling, and 40/913 to chemokine signaling were impacted in the transgenic spleen. Similarly, several DEGs in the pathways of the complement & coagulation cascade, the Toll-like receptor (TLR), the antigen processing & presentation and Fc ε RI signaling & Fc γ RI-mediated phagocytosis were greatly altered in the spleen. In the kidney, DEGs in the leukocyte trans-endothelial migration and hematopoietic cell lineage pathways were impacted (Table 3). By analyzing liver-associated pathways, DEGs in the pathways of phagosomal activity, fatty acid biosynthesis, and complement & coagulation cascade were significantly perturbed (Table 3). Interestingly, highly disturbed PPAR signaling pathways were discovered in both liver and kidney, and this result may indicate a major alteration in regulating energy efficiency occurred in transgenic rainbow trout.
Positively charged and amphipathic AMPs, e.g. cecropin P1, have been characterized for their activities of cytotoxic elimination of bacteria and viruses as their primary function in host innate immunity . As reviewed by Hilchie et al. , AMPs were also found to have multifaceted immunomodulatory effects in a variety of different hosts. The immunomodulatory activities of these AMPs were found to exert differential expression of cytokines and chemokines, differentiation of leukocytes, and elevation of damage repairing and so on in animal models [11, 24, 25]. Recently, transgenic rainbow trout bearing cecropin P1 transgene was produced in our laboratory, and the elevation of resistance to bacterial, viral and parasitic infections were also demonstrated through repeated challenge studies . Therefore, it is conceivable to use these transgenic fish as experimental animals to address the question whether cecropin P1 transgene product can modulate host immune system or other genetic traits that may lead to increasing resistant to microbial infections in the host. In a preliminary cDNA microarray study conducted Lo et al. , they had unveiled that several functional alterations on the expression of immune related genes have taken place in the spleen, liver, and kidney of the transgenic fish. However, due to the inherited disadvantages associated with the cDNA microarray study, as mentioned in the Introduction, a confirmation study via a second experimental approach is required in order to provide a more solid evidence to support our hypothesis, and thus RNA-Seq analysis was adopted in this study .
Spleen: The major innate/adaptive immune relevant organ
The major functions of the spleen in fish had been well reviewed by scientists [26,27,28,29,30]. Very much alike to other higher vertebrates, the spleen of teleost fish is merged in a complex of splenic ellipsoids, melanomacrophages centers (MMCs) and lymphoid tissue , this complex filters and traps blood cells, and processes antigens . Other functions of fish spleen are also known to include immune memories mediated by presenting lymphocytic antigen , phagocytosis of macrophages stimulated by complement receptors to C3 signaling , and naïve CD4+ T cells differentiating into T helper cells induced by major histocompatibility complex (MHC) . Since spleen plays critical roles in fish immunity, we believe that one of the major functional perturbations in cecropin P1 transgenic rainbow trout may occur in the spleen. To address this hypothesis, immune relevant pathways were generated by a custom made visualization tool from Shin et al. , and mapped with RPKM ratios of pre-defined transgenic spleen DEGs. As excepted, significant alterations in the cytokine/chemokine mediated JAK-STAT signaling pathway (Fig. 4) was observed in the spleen of transgenic rainbow trout. In addition, the heat-map analysis revealed a significant enhancement of cytokine and chemokine signaling via ligands and receptors binding, and followed by up-regulation of down-stream JAK kinases and STAT transcription factors. In addition, the STATs target genes that correlated to different biological processes, namely immune response (GO: 0006955), response to viruses (GO: 0009615), positive regulation of immune response (GO: 0050778, GO: 0002684), and JAK-related cytokine signaling (PMID: 27044867), were also shown partially elevated in Fig. 4. Such results may indicate that the splenoid lymphocytes may produce more chemokines and cytokines, and these ligands may provide positive feedback looping regulation of the proliferation and differentiation of lymphocytes in the spleen, and resulting in an overall enhancement of immune response in transgenic rainbow trout. Similar RNA-Seq results of tilapia spleen were reported by Zhu et al.  that the Streptococcus iniae challenged fish were shown significantly perturbed of expression levels of C-C/C-X-C motif chemokines, C-X-C chemokines receptors, interleukins, interferons, and MHC class I/II antigens at its acute phase after infection. Additionally, the transcriptomic analysis of spleen of Atlantic salmon conducted by Dettleff et al. , demonstrated that the fish altered their expression profiles of interleukins/interleukin receptors (il10, il10rb, and il13ra2), interferon/interferon receptors (ifng and ifngr1), in responding of viral infection of infectious salmon anemia virus (ISAV). In our transgenic rainbow trout, the Toll-like receptor signaling was also highly promoted in the spleen of transgenic fish by establishing up-regulated genes of Toll-like receptors (tlr1, 2, 3, 7 and 8), and the important adaptor protein (myd88). The expression levels of the down-stream kinases (e.g. tbk1, p38 and pi3k), transcription factors (irf3, irf7 and stat1), and targeting genes (cd40, cxcl10 and cxcl11) were all elevated (Additional file 3: Figure S4). The stimulated Toll-like receptor signaling may result in an enhancing chemotaxis of immune cells and increasing pro-inflammatory effects via Toll-like to NF-κB pathway. Moreover, complement components of c1q, c4, c3ar1, c5ar1, c6 and c7 were found up-regulated in the complement and coagulation pathway of the transgenic spleen (Additional file 3: Figure S3), and may indicate an increasing expression of complement complex in the transgenic fish. Finally, the surface integrins (itga4 and itgb1) of leukocytes, essential matrix metalloproteinases (mmp2 and mmp9), docking actins (actn1 and actb1), and recognition receptor: cxcr4 were increased their expression levels in leukocyte trans-endothelial migration pathway (Additional file 4: Figure S6). These findings may suggest that leukocytes might have higher recognition and adhesion abilities targeting onto epitheliums. The overall observed perturbation in the spleen by the cecropin P1 transgene product strongly support our hypothesis that enhancement of the spleen function leads to acceleration of innate/adaptive immune responses in cecropin P1 transgenic rainbow trout.
Kidney: The major hematopoiesis relevant organ
While in the absence of bone marrow, the anterior kidney of teleost fish plays an equivalent role of hematopoiesis as bone marrow in higher vertebrates. For instance, B-lymphocytes are produced in fish kidney, as well as myeloid lineage cells (e.g. monocytes, macrophages, and granulocytes) [34, 35]. In the current study, three important functional perturbations were observed in the kidney of the transgenic rainbow trout, and the spraying RPKM ratios of DEGs to custom-made pathways for detail evaluation were performed as well. For the kidney, the gene expression profiles of (1) hematopoietic cell lineage, (2) leukocyte trans-endothelial migration, and (3) PPAR signaling pathways were obviously altered. In the hematopoietic cell lineage, the expression of many lymphocytic associated surface antigen genes were down regulated. For example, down regulation of the expression cd2/cd3 (pro-T cell) and cd8 (T cell) genes was observed. In addition, the expression of cd9/cd22 (pre/pro-B cell), cd37 (immature B cell), igm/igd (B cell), cd11b/cd13 (pro-monocyte) and cd115 (monocyte) genes was down-regulated as well. By combining the finding of cDNA microarray analysis conducted by Lo et al.  and the results of the current study, a suppressed translational machinery is suggested and further reduced lymphocytic proliferation and differentiation may be inferred in transgenic kidney. However, enhancement of maturation and activation of lymphocytes were observed in transgenic spleen (see discussed in the spleen section), so the results from the kidney could be counter balanced by the spleen effects, and to a homeostasis of lymphocytic function been maintained in the transgenic rainbow trout. In contrast to spleen, kidney was demonstrated a suppressed leukocyte trans-endothelial migration in transgenic rainbow trout. Namely, surface integrins (itga4, itgb7), recognition antigens (cd11a and cd11b), and docking actins (actg1 and actn1) were down-regulated in the transgenic kidney (Additional file 4: Figure S7). In short, the lacking of surface antigens of lymphocytes and leukocytes may be the result from decreasing production of blood cells in the early stages of hematopoietic lineage in the transgenic kidney. Our results were supported by a cross-species transcriptomic analysis conducted by Xiang et al. . RNA samples from head kidney and spleen of Vibrio harveyi challenged Japanese sea bass (Lateolabrax japonicus) were sequenced, and functional perturbations of antigen presenting and processing (mhc1/2, hsp70, hsp90, and calnexin), pattern recognition (collectin12 and itga3) and regulators of hematopoiesis (gfi1 and gfi1b) were reported by Xiang et al. . Finally, the PPAR signaling pathway was significantly promoted in the transgenic kidney, and the overall effects of enhancing PPAR signaling will be discussed in the liver section.
Liver: The major energy efficiency relevant organ
As reviewed by Jenne et al. , liver was identified as one of the most important immune relevant organs in mammals . Alike to mammals, fish liver also plays critical roles in their immunity . Additionally, liver is responsible for many energy metabolism processes as well, including metabolic protein synthesis, degradation, and fatty acids (FA) biosynthesis . According to Bransden et al. , salmonid immunity was highly correlated to the alteration of their FA compositions. In addition, interrupting fish energy efficiency via starvation  or treating fish with additional fatty acids  could change the innate and adaptive immune responses in Atlantic salmon. In the study by Lo et al. , DEGs contributing to energy metabolism pathways were observed in the transgenic liver, and they suggested that it might indirectly affect fish immunity. Here, our RNA-Seq results confirmed the findings of Lo et al.  that pathways of PPAR signaling and FA biosynthesis were greatly perturbed in the transgenic liver. For PPAR signaling (Additional file 3: Figure S5), FA transport protein (slc27a1) and binding protein (fabp1) were both up-regulated. Moreover, the expression levels of target genes that are responsible for lipid transportation (pltp and apoa1), FA oxidation (acadm), and adipocyte differentiation (acdc and angptl4) were elevated as well. An additional supportive finding was provided by other scientists in the RNA-Seq transcriptomic analysis of the liver of Edwardsiella tarda vaccine immunized zebrafish . Yang et al.  reported differentially expression levels of energy metabolism associated genes such as apolipoprotein (apoa4), plasminogen (plg), fibrinogen (fgg) and serum amyloids (apcs and saa) in the liver of zebrafish (Danio rerio) at its acute phase after vaccine immunization. Very interestingly, the PPAR signaling was also perturbed in the transgenic kidney (Additional file 4: Figure S9): slc27a1, fabp1, acdc, angptl4 and acadm were consistently up-regulated in both the liver and kidney. Although some of the DEGs involving in the PPAR signaling are down-regulated, the up-regulated DEGs suggest an overall elevation of lipid metabolism in both the transgenic liver and kidney, and, therefore, an enhancing lipid energy metabolism may be expected in the transgenic fish. As revealed by Yang et al. , the vaccine immunized zebrafish altered the expressions of lysozyme (lygl1) and lysosomal membrane protein 2 (lamp2) in the liver. Here, we discovered the promotion of phagosomal and lysosomal activity in the liver of cecropin P1 transgenic rainbow trout. The collectin (colec11) and C-lectin receptors (cd206 and cd209) were stimulated as well as mhc1, and genes correlated to lysosomal activity were also perturbed such as rilp, lamp, and ctsl (Additional file 4: Figure S8). These results may indicate an increasing endoplasmic reticulum mediated phagocytosis, which play an essential role in lysosomal elimination of pathogens, in the transgenic liver.
In this study, we confirmed the findings of cDNA microarray analysis conducted by Lo et al.  in a cross-platform approach, RNA-Seq, and more solid evidence were obtained to support our hypothesis. The cecropin P1 transgene not only directly eliminate bacterial, virial and parasitic pathogens, but also modulate the host innate and adaptive immunity in rainbow trout. Many functional perturbations found in the spleen and kidney might contribute to directly increasing fish immune responses; some others in the liver and kidney might provide indirect effects via altering the energy metabolisms. The genetic data shown in this study provides evidence, expands pool of knowledge about fish immunity, and may benefit further studies in marine species. Finally, the identified DEGs committing to achieve immunomodulation may further serve as biological markers for artificial selection and breeding of disease resistant fish species in aquaculture.
Three immune competent tissues, namely spleen, liver and kidney, were harvested from two individual fish of two transgenic families respectively and one non-transgenic family (one year old) maintained in the Salmon Disease Laboratory at Oregon State University (OSU protocol # 4282). Prior to tissue collection, fish were euthanized by treatment with MS222 (Tricaine mesylate, Sigma-Aldrich, St. Louis, Mo) following the specification in OSU protocol #4282. For each tissue, one RNA sample was prepared from a pooled tissue of two individual fish of the control population (non-transgenic fish) and two RNA samples each was prepared from a pooled tissue of two individual fish of two independent transgenic families. Total RNA samples were prepared by the TRIzol extraction method following the protocol provided by the manufacturer (Invitrogen, Carlsbad, CA). By treating all RNA samples with RNase free DNase-I (M610A, Promega, Madison, WI), the genomic DNA contamination was removed. The quality of RNA samples were assessed by an Agilent Bioanalyzer (Agilent Technologies, Santa Clara, CA), and concentrations by a NanoDrop spectrophotometer (NanoDrop Technologies, Wilmington, DE).
Transcriptomic sequencing and reads assembly
Illumina second generation sequencing was performed by Beijing Genomic Institute (BGI) as a commercial service. RNA samples of the non-transgenic rainbow trout were used to generate reference database. Poly(A)+-RNA was isolated from the total RNA by magnetic beads conjugated with oligo(dT). The poly(A)+-RNA was fragmented to about 200 bases, and subjected to first strand cDNA synthesis using random hexamer as primers. After addition of buffer, dNTPs, RNase H and DNA polymerase I, the second strand cDNA was synthesized. The double-strand cDNA was purified, and modified by 5′-end phosphorylation and 3′-end addition of single nucleotide adenine. After ligation to the sequencing adaptors, the suitable cDNA fragments were selected as templates for the PCR amplification, and the products were sequenced by Illumina HiSeq™ 2000. After sequencing, clean reads were collected by removal reads of adaptors, reads containing > 5% unknown nucleotides and reads with low quality bases (base quality ≦ 10). Subsequently, the transcriptomic de novo assembly was carried out with the short reads assembling program, Trinity  to generate unigenes. After generating reference database from non-transgenic rainbow trout, tissue specific RNA samples from two transgenic families were sequenced by the same procedure. By removing the noise, reads of adaptor, reads of unknown bases > 10%, and reads of low quality bases (base quality ≦ 5) greater than 50%, clean reads were produced and mapped to reference database with SOAPaligner/SOAP2 . During the mapping, no more than two mismatches were allowed in the alignment.
Unigene annotation and functional classification
Homology searches of the assembled unigenes were performed against public protein databases of Nr, Swiss-Prot, KEGG, and COG by blastx (the e-value < 10− 5). In addition, functional classification of unigenes was also conducted via Blast2GO program  for non-redundant GO annotation, and web gene ontology annotation plot (WEGO) program  for GO classification.
Differential expression of Unigenes
The gene expression level was calculated by parameter of RPKM . RPKM is defined as RPKM(Χ) = 106*C/(NL*10− 3), where RPKM(Χ) stands for the expression of gene Χ, C equals to the number of reads that specifically aligned to gene Χ, N is the total number of reads that aligned to all genes, and L stands for the number of bases of gene Χ. After determining RPKM of unigenes, the base two logarithm of RPKM ratio of transgenic over non-transgenic group was calculated. For those unigenes with RPKM equals to zero were excluded, and the unigenes with RPKM ratio greater than 2 folds, i.e., log2(transgenic/non-transgenic > 1 or < − 1), were defined as differentially expressed genes (DEGs). If there was more than one unigene relevant to a gene of interest, the longest unigene was used to determine the expression level and coverage. For screening the significance of gene expression, strict algorithm was customized by BGI, as a commercial service, based on the general method of Audic et al. . A false detection rate (FDR) was pre-set to a number no larger than 0.01, following the procedure of Benjanini-Hochberg  which was used to determine the threshold of P-value in multiple tests of differential expression genes.
Reverse transcription (RT) quantitative real-time PCR (RT-qPCR)
Superscript III reverse transcriptase (18080–044, Life Technology) and oligo-(dT)18 were used to reverse transcribe two μg of DNase-treated RNA samples into first strand cDNA by following conditions provided by the manufacturer. The resulting cDNA products were diluted with DNase-free water into final volume of 100 μL. The qPCR assay was conducted in a 96 wells plate in the C1000 thermal cycler/CFX96 Real-Time PCR Detection System (Bio-Rad, Hercules, CA). The reaction mixture, 20 μL, contains 1 μL template, 0.5 μM gene specific primers (Additional file 5), 1x SsoFast EvaGreen Supermix (172–5201, Bio-Rad), and 0.01 μM fluorescein (170–8780, Bio-Rad). The amplification program consisted of initial denaturing for 2 min at 98 °C, followed by 40 cycles of 5 s at 98 °C and 30 s at 59 °C for annealing and synthesis. For quality control of the amplified products, a melting curve, 65 °C to 95 °C with 0.5 °C increments every 5 s, was performed after each amplification. Subsequently, the threshold cycle (Ct) values were collected by CFX manager (Bio-Rad) and analyzed by standard normalizer of 2-ΔΔCt to determine the differential expression (transgenic vs non-transgenic). For the assay of each gene, the efficiencies of reverse transcription and PCR were determined to be within the limit of reproducibility.
In the enrichment analysis by GeneCodis , both pre-defined DEGs of each tissue from the transgenic fish and the total genes of the non-transgenic fish library, serving as a reference gene list, were inputted to inquire GO: Biological Process, GOSlim Process and KEGG pathways. The default statistical parameters (minimum number of genes: 3, statistical test: hypergeometric and chi-square; p-value correction = false detection rate: FDR) were applied for both singular and modular enrichment analyses.
For the software versions, databases and parameters, utilized in this study, were summarized in Additional file 6.
Beijing Genomic Institute
Nucleotide Basic Local Alignment Search Tool
Clusters of orthologous
- cox2 :
Differentially expressed gene
False detection rate
Gene Ontology Biological Process
Infectious hematopoietic necrosis virus
Infectious Hematopoietic Necrosis Virus
- il1β :
Infectious salmon anemia virus
Kyoto Encyclopedia of Genes and Genomes
Major histocompatibility complex
Next generation sequencing
Peroxisome proliferator-activated receptor
Quantitative reverse transcription poly chain reaction
Reads per kilobase of transcript per million mapped reads
Short oligonucleotide analysis package
Signal transducer and activator of transcription
Web gene ontology annotation plot
Chen TT, Lin CM, Chen MJ, Lo JH, Chiou PP, Gong HY, Wu JL, Chen MHC, Yarish C. In: Kim S-K, editor.Transgenic technology in marine organisms. Springer handbook of marine biotechnology. Berlin, Heidelberg: Springer; 2015. p. 387–412.
Boman HG, Nilsson-Faye I, Paul K, Rasmuson T Jr. Insect immunity. I. Characteristics of an inducible cell-free antibacterial reaction in hemolymph of Samia cynthia pupae. Infect Immun. 1974;10(1):136–45.
Shai Y. Mechanism of the binding, insertion and destabilization of phospholipid bilayer membranes by alpha-helical antimicrobial and cell non-selective membrane-lytic peptides. Biochim Biophys Acta. 1999;1462(1–2):55–70.
Boman HG. Peptide antibiotics and their role in innate immunity. Annu Rev Immunol. 1995;13:61–92.
Zhang G, Ross CR, Blecha F. Porcine antimicrobial peptides: new prospects for ancient molecules of host defense. Vet Res. 2000;31(3):277–96.
Sarmasik A, Warr G, Chen TT. Production of transgenic medaka with increased resistance to bacterial pathogens. Mar Biotechnol (NY). 2002;4(3):310–22.
Chiou PP, Chen MJ, Lin CM, Khoo J, Larson J, Holt R, Leong JA, Thorgarrd G, Chen TT. Production of homozygous transgenic rainbow trout with enhanced disease resistance. Mar Biotechnol (NY). 2014;16(3):299–308.
Chen TT, Chen MJ: Transgenic Fish Technology: Principle and Application. In “Encyclopedia of Marine Biotechnology” (ed by S.K. Kim), Wiley Publisher. 2017 (in press).
Chiou PP, Lin CM, Perez L, Chen TT. Effect of cecropin B and a synthetic analogue on propagation of fish viruses in vitro. Mar Biotechnol (NY). 2002;4(3):294–302.
Chiou PP, Khoo J, Bols NC, Douglas S, Chen TT. Effects of linear cationic alpha-helical antimicrobial peptides on immune-relevant genes in trout macrophages. Dev Comp Immunol. 2006;30(9):797–806.
Lai Y, Gallo RL. AMPed up immunity: how antimicrobial peptides have multiple roles in immune defense. Trends Immunol. 2009;30(3):131–41.
Lo JH, Lin CM, Chen MJ, Chen TT. Altered gene expression patterns of innate and adaptive immunity pathways in transgenic rainbow trout harboring Cecropin P1 transgene. BMC Genomics. 2014;15:887–900.
Jantzen SG, Sanderson DS, von Schalburg KR, Yasuike M, Marass F, Koop BF. A 44K microarray dataset of the changing transcriptome in developing Atlantic salmon (Salmo salar L.). BMC Res Notes. 2011;4:88.
Zhao S, Fung-Leung WP, Bittner A, Karen N, Liu X. Comparison of RNA-seq and microarray in transcriptomic profiling of activated T cells. PLoS One. 2014;9(1):e78644.
Costa-Silva J, Domingues D, Lopes FM: RNA-Seq differential expression analysis: an extended review and software tool. PLoS One 2017, 12: e0190152. http//doi.org/10.137/journaal.pone.0190152.
Mortazavi A, Williams BA, McCue K, Schaeffer L, Wold B. Mapping and quantifying transcriptomes by RNA-Seq. Nat Methods. 2008;5:621–8.
Cones A, Madrigal P, Tarazona S, Gomez-Cabrero D, Cervera A, McPherson A, Szczesniak MW, Gaffiney DJ, Elo LL, Zhang X. A survey of best practice for RNA-Seq data analysis. Genome Biol. 2016;17:13–31.
Grabherr MG, Haas BJ, Yassour M, Levin JZ, Thompson DA, Amit I, Adiconis X, Fan L, Raychowdhury R, Zeng Q, Chen Z, Mauceli E, Hacohen N, Gnirke A, Rhind N, di Palma F, Birren BW, Nusbaum C, Lindblad-Toh K, Friedman N, Regev A. Full-length transcriptome assembly from RNA-seq data without a reference genome. Nat Biotechnol. 2011;29:644–52.
Carmona-Saez P, Chagoyen M, Tirado F, Carazo JM, Pascual-Montano A. GENECODIS: a web-based tool for finding significant concurrent annotations in gene lists. Genome Biol. 2007;8(1):R3.
Tabas-Madrid D, Nogales-Cadenas R, Pascual-Montano A. GeneCodis3: a non-redundant and modular enrichment analysis tool for functional genomics. Nucleic Acids Res. 2012. https://doi.org/10.1093/nar/gks402.
Li R, Yu C, Li Y, Lam TW, Yiu SM, Kristiansen K, Wang J. SOAP2: an improved ultrafast tool for short read alignment. Bioinformatics. 2009;25(15):1966–7.
Mortazavi A, Williams BA, McCue K, Schaeffer L, Wold B. Mapping and quantifying mammalian transcriptomes by RNA-seq. Nat Meth. 2008;5(7):621–8.
Brogden KA, Ackermann M, McCray PB Jr, Tack BF. Antimicrobial peptides in animals and their role in host defences. Int J Antimicrob Agents. 2003;22(5):465–78.
Hilchie AL, Wuerth K, Hancock RE. Immune modulation by multifaceted cationic host defense (antimicrobial) peptides. Nat Chem Biol. 2013;9(12):761–8.
Scott MG, Dullaghan E, Mookherjee N, Glavas N, Waldbrook M, Thompson A, Wang A, Lee K, Doria S, Hamill P, Yu JJ, Li Y, Donini O, Guarna MM, Finlay BB, North JR, Hancock RE. An anti-infection peptide that selectively modulates the innate immune response. Nat Biotechnol. 2007;25(4):465–72.
Uribe C, Folch H, Enriquez R, Moran G. Innate and adaptive immunity in teleost fish: a review. Vet Med. 2011;56(10):486–503.
Zapata AG, Chiba A, Varas A. Cells and tissues of the immune system of fish. San Diego: The Fish Immune System Organism, Pathogen and Environment Academic Press; 1996. p. 1–62.
Espenes A, Press CM, Dannevig BH, Landsverk T. Immune-complex trapping in the splenic ellipsoids of rainbow trout (Oncorhynchus mykiss). Cell Tissue Res. 1995;282(1):41–8.
Johnson E, Smith P. Attachment and phagocytosis by salmon macrophages of agarose beads coated with human C3b and C3bi. Dev Comp Immunol. 1984;8(3):623–30.
Vallejo AN, Miller NW, William Clem L. Antigen processing and presentation in teleost immune responses. Annu Rev Fish Dis. 1992;2(0):73–89.
Zhao Y, Hoang TH, Joshi P, Hong SH, Giardina C, Shin DG. A route-based pathway analysis framework integrating mutation information and gene expression data. Methods. 2017;124:3–12.
Zhu J, Li C, Ao Q, Tan Y, Luo Y, Guo Y, Lan G, Jiang H, Gan X. Transcriptomic profiling revealed the signatures of acute immune response in tilapia (Oreochromis niloticus) following Streptococcus iniae challenge. Fish Shellfish Immunol. 2015;46(2):346–53.
Dettleff P, Moen T, Santi N, Martinez V. Transcriptomic analysis of spleen infected with infectious salmon anemia virus reveals distinct pattern of viral replication on resistant and susceptible Atlantic salmon (Salmo salar). Fish Shellfish Immunol. 2017;61:187–93.
Orkin SH, Zon LI. Hematopoiesis: an evolving paradigm for stem cell biology. Cell. 2008;132(4):631–44.
Zwollo P, Cole S, Bromage E, Kaattari S. B cell heterogeneity in the teleost kidney: evidence for a maturation gradient from anterior to posterior kidney. J Immunol. 2005;174(11):6608–16.
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–83.
Jenne CN, Kubes P. Immune surveillance by the liver. Nat Immunol. 2013;14(10):996–1006.
Bruslé J, Anadon GG. The structure and function of fish liver. Fish mMorphology: Science Publishers; 1996. p. 77–93.
Bransden MP, Carter CG, Nichols PD. Replacement of fish oil with sunflower oil in feeds for Atlantic salmon (Salmo salar L.): effect on growth performance, tissue fatty acid composition and disease resistance. Comp Biochem Physiol B Biochem Mol Biol. 2003;135(4):611–25.
Martin SA, Douglas A, Houlihan DF, Secombes CJ. Starvation alters the liver transcriptome of the innate immune response in Atlantic salmon (Salmo salar). BMC Genomics. 2010;11:418.
Martinez-Rubio L, Morais S, Evensen O, Wadsworth S, Vecino JG, Ruohonen K, Bell JG, Tocher DR. Effect of functional feeds on fatty acid and eicosanoid metabolism in liver and head kidney of Atlantic salmon (Salmo salar L.) with experimentally induced heart and skeletal muscle inflammation. Fish Shellfish Immunol. 2013;34(6):1533–45.
Yang D, Liu Q, Yang M, Wu H, Wang Q, Xiao J, Zhang Y. 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. BMC Genomics. 2012;13:319–35.
Conesa A, Gotz S, Garcia-Gomez JM, Terol J, Talon M, Robles M. Blast2GO: a universal tool for annotation, visualization and analysis in functional genomics research. Bioinformatics. 2005;21(18):3674–6.
Ye J, Fang L, Zhang Y, Chen J, Zhang Z, Wang J, Li S, Li R, Bolund L, Wang J. WEGO: a web tool for plotting GO annotations. Nucleic Acids Res. 2006;34:W293–7.
Audic S, Claverie JM. The significance of digital gene expression profiles. Genome Res. 1997;7(10):986–95.
Benjamini Y, Yekutieli D. The control of the false discovery rate in multiple testing under dependency. Ann Stat. 2001;29(4):1165–88.
We thank members at the SDL of OSU for assistance in maintaining transgenic rainbow trout. We also thank Dr. D.G. Shin and his group at CSE of UCONN for helping develop the visualization tool for pathway analysis.
This research was supported by a grant from the U.S. Department of Agriculture (CONTR 58–1930–0-009) to T.T.C.
Availability of data and materials
The datasets generated and analyzed during the current research are available in GEO database as accession: GSE111373. The raw data (fastq format), and the aligning library (fasta format) can be inquired from: https://www.ncbi.nlm.nih.gov/geo/query/acc.cgi?acc=GSE111373
Ethics approval and consent to participate
Both transgenic and non-transgenic rainbow trout (one year old) were maintained in the Salmon Disease Laboratory at Oregon State University by following OSU protocol # 4282.
Consent for publication
The authors declare that there is no conflict of interest that would prejudice the impartiality of this scientific research.
Supplemental figures for pie charts of qualities of clean reads and results of sequencing saturation analyses to show quality of aligning, and supplemental table to show overall mapping results (to reference genes) per each sample. (PDF 206 kb)
Singular enrichment analyses among three tissues, grouped by both GOSlim biological process and KEGG terms. (PDF 121 kb)
Supplement figures for pathways generated from custom-made visualization tool. (PDF 743 kb)
Supplement figures from KEGG database and sprayed by the DEGs expression profiles. (PDF 240 kb)
Primers for RT-qPCR analysis. (XLSX 11 kb)
Summarizing of software versions, databases and parameters for Bioinformatics analyses. (XLSX 9 kb)
About this article
Cite this article
Han, YC., Lin, CM. & Chen, T.T. RNA-Seq analysis of differentially expressed genes relevant to innate and adaptive immunity in cecropin P1 transgenic rainbow trout (Oncorhynchus mykiss). BMC Genomics 19, 760 (2018). https://doi.org/10.1186/s12864-018-5141-8
- Oncorhynchus mykiss
- Disease resistant transgenic fish
- Antimicrobial peptide
- Cecropin P1