Next Article in Journal
Constraints on the Genetic and Antigenic Variability of Measles Virus
Previous Article in Journal
Validation of the Filovirus Plaque Assay for Use in Preclinical Studies
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

De Novo Transcriptome Analysis Shows That SAV-3 Infection Upregulates Pattern Recognition Receptors of the Endosomal Toll-Like and RIG-I-Like Receptor Signaling Pathways in Macrophage/Dendritic Like TO-Cells

Section of Aquatic Medicine and Nutrition, Department of Basic Sciences and Aquatic Medicine, Faculty of Veterinary Medicine and Biosciences, Norwegian University of Life Sciences, Ullevålsveien 72, P.O. Box 8146 Dep NO-0033 Oslo, Norway
*
Author to whom correspondence should be addressed.
Submission received: 15 January 2016 / Revised: 5 April 2016 / Accepted: 14 April 2016 / Published: 21 April 2016
(This article belongs to the Section Animal Viruses)

Abstract

:
A fundamental step in cellular defense mechanisms is the recognition of “danger signals” made of conserved pathogen associated molecular patterns (PAMPs) expressed by invading pathogens, by host cell germ line coded pattern recognition receptors (PRRs). In this study, we used RNA-seq and the Kyoto encyclopedia of genes and genomes (KEGG) to identify PRRs together with the network pathway of differentially expressed genes (DEGs) that recognize salmonid alphavirus subtype 3 (SAV-3) infection in macrophage/dendritic like TO-cells derived from Atlantic salmon (Salmo salar L) headkidney leukocytes. Our findings show that recognition of SAV-3 in TO-cells was restricted to endosomal Toll-like receptors (TLRs) 3 and 8 together with RIG-I-like receptors (RLRs) and not the nucleotide-binding oligomerization domain-like receptors NOD-like receptor (NLRs) genes. Among the RLRs, upregulated genes included the retinoic acid inducible gene I (RIG-I), melanoma differentiation association 5 (MDA5) and laboratory of genetics and physiology 2 (LGP2). The study points to possible involvement of the tripartite motif containing 25 (TRIM25) and mitochondrial antiviral signaling protein (MAVS) in modulating RIG-I signaling being the first report that links these genes to the RLR pathway in SAV-3 infection in TO-cells. Downstream signaling suggests that both the TLR and RLR pathways use interferon (IFN) regulatory factors (IRFs) 3 and 7 to produce IFN-a2. The validity of RNA-seq data generated in this study was confirmed by quantitative real time qRT-PCR showing that genes up- or downregulated by RNA-seq were also up- or downregulated by RT-PCR. Overall, this study shows that de novo transcriptome assembly identify key receptors of the TLR and RLR sensors engaged in host pathogen interaction at cellular level. We envisage that data presented here can open a road map for future intervention strategies in SAV infection of salmon.

1. Introduction

A crucial step in cellular defense mechanisms against viral infection is recognition of danger signals that initiate signaling pathways aimed at protecting host cells against pathogen invasion [1]. Apart from protecting host cells, recognition of microbial danger signals is a crucial step for targeted delivery of vaccine antigens into antigen presenting cells (APCs) as recently pointed out by Munang’andu and Evensen [2]. The major players in recognition of microbial invasion are pattern recognition receptors (PRRs) made of germ line coded receptors that recognize conserved microbial features called “pathogen associated molecular patterns” (PAMPs) [3]. In addition, PRRs also recognize endogenous host structures released after tissue damage called “damage associated molecular patterns” (DAMPs) [4]. The numbers of germ line coded PRRs is limited and as such PAMPs represent unique structures that are characteristic of several groups of pathogens.
Currently, there are different PRR families identified in vertebrates that serve as immune sensors of PAMPs and these include the Toll-like receptors (TLRs), retinoic acid-inducible gene I (RIG-I)-like receptors (RLRs), nucleotide oligomerization domain (NOD)-like receptors (NLRs), as well as the melanoma 2 (AIM2) like receptors (ALRs) and the cytoplasmic double stranded DNA sensors (CDSs) [5,6,7,8]. Members of the TLR family detect PAMPs from protozoa, bacteria, fungi and viruses and they can broadly be classified into TLRs found on cell surfaces and those found in endosomal compartments [1,9,10]. The NLRs with known functions mainly recognize bacteria while RLR are antiviral [11]. Thus far, 17 members of the TLRs have been identified in different fish species [12,13] while a genomic overview of NLRs found in fish was recently published by Laing et al. [14]. In addition, MDA5 and RIG-I receptors were recently cloned and characterized in salmonids [15]. However, there is little information regarding the signaling pathways induced by these PRRs in different fish species although several genes involved in the downstream signaling of PRR pathways have been cloned and characterized in different fish species [16,17]. One of the major drawbacks to elucidating the signaling pathways induced by different PRRs in fish is the general absence of tools such as knockout models that can be used to elucidate the functional roles of different genes expressed at different stages of the signaling pathways. While the search for signaling pathway analytical tools continues in fish, the emergence of RNA-seq has opened a new dimension in functional genomics in which a vast array of genes expressed in response to host–pathogen interaction can be sequenced at the same time thereby allowing for a global understanding of cellular responses induced by microbial invasion at transcript level [18]. As such, genome wide transcriptome data analysis can be used to identify networks of genes expressed in response to microbial invasions at the same time.
Hence, in the present study, we used a de novo assembly to generate a transcriptome of differentially expressed genes (DEGs) generated in response to salmonid alphavirus subtype 3 (SAV-3) infection in macrophages/dendritic like TO-cells derived from Atlantic salmon (Salmo salar L) headkidney leukocytes [19,20]. By using the Kyoto encyclopedia of genes and genomes (KEGG) pathway analysis, we wanted to find out the repertoire of genes linked to PRR pathways induced by SAV-3 infection in TO-cells. SAV-3 is the etiological agent for pancreas disease (PD) known to cause high economic losses in salmonids [21,22]. It is a member of the genus alphavirus in the family Togaviridae [23]. It contains a positive sense single stranded RNA (+ssRNA) genome with capped 5′ end and polyadenylated 3′ end that serves directly as messenger RNA (mRNA) for the translation of viral non-structural proteins upon entry and form the dsRNA intermediate during replication in infected cells [24,25]. As pointed out elsewhere [26,27], ssRNA is sensed by RIG-I and TLR-7/8, while dsRNA is sensed by TLR3, RIG-I and MDA5 [27,28,29,30]. The SAV-3 genome is subdivided into two open reading frames (ORFs). The first ORF encodes four non-structural proteins (nsPs) designated as nsP1-4 responsible for the transcription and replication of the viral RNA while the second ORF encode the structural proteins PE2-6K-E1 [31]. Based on the transcriptome analysis presented here, we demonstrate that the repertoire of PRR genes expressed in response to SAV-3 infection in TO-cells is comparable to the profile of genes linked to PRR signaling pathways induced by other alphavirus infections in mammalian cells. Further, we also show that pathway based analysis provides a contextual understanding of the biological relevance of DEGs expressed in a transcriptome. We envision that data presented here shall broaden our understanding of the cellular mechanisms used by fish cells to combat microbial invasion.

2. Materials and Methods

2.1. Cell Culture and Virus Infection

TO-cells derived from Atlantic salmon headkidney leukocytes characterized to possess macrophage/dendritic cell like properties [19,20], were propagated at 20 °C in HMEM (Eagle’s minimal essential medium (MEM) with Hanks’ balanced salt solution (BSS)) supplemented with l-glutamine, MEM nonessential amino acids, gentamicin sulphate, and 10% fetal bovine serum (FBS). When the cells were 80% confluent, one batch was inoculated with SAV-3 (Genebank accession JQ799139) [32] at multiplicity of infection (MOI) 1 while another batch was only exposed to the HMEM growth media. Thereafter, both the SAV-3 infected and non-infected TO-cells were incubated at 15 °C in HMEM maintenance media supplemented with 2% FBS. Cells from both the infected and non-infected groups were harvested after 48 h. Both the SAV-3 infected and non-infected cells were propagated in triplicates.

2.2. Total RNA Isolation

Extraction of total RNA from SAV-3 infected and non-infected TO-cells was carried out using the RNAeasy mini kit with on-column DNase treatment according to the manufacturers’ instructions (Qiagen, Hilden, Germany). The quality and concentration of RNA was analyzed using the ND1000 nanodrop (Nanodrop Technologies, Wilmington, NC, USA) and Agilent 2100 Bioanalyzer (Agilent Technologies, Santa Clara, CA, USA).

2.3. Library Construction, Sequencing and Data Analysis for RNA-Seq

Library construction was carried out by pooling together triplicate samples obtained from total RNA extraction of SAV-3 infected and non-infected cells for RNA-Seq. Treatment of total RNA with DNase I to degrade any possible DNA contamination, enrichment using oligo(dT) magnetic beads, fragmentation into approximately 200 bp fragments, synthesis of first strand cDNA using random hexamer-primers followed by synthesis of the second strand together with end reparation coupled with 3′-end single nucleotide A (adenine) addition, ligation of sequence adaptors to the fragments and fragment enrichment by PCR amplification were also carried out as previously described in our studies [17]. Thereafter, quality check (QC step) was carried out using the Agilent 2100 Bioanaylzer and ABI StepOnePlus Real-Time PCR System (Bio-Rad) to qualify and quantify the sample library. Subsequently, library products were used for RNA-sequencing using Illumina HiSeqTM 2000, BGI-Hong Kong and clean reads were obtained after removal of adaptor sequences together with reads having >10% of unknown bases and reads with low quality bases (base with quality value ≤5) >50% in a read.

2.4. De Novo Assembly, Functional Annotation and Gene Ontology Classification

Once a library of clean reads was prepared, clean reads were then used for de novo transcriptome assembly using the Trinity software [33]. Thereafter, the assembled unigenes were annotated into different functional classifications after searching in different protein databases using the BlastX (version 2.2.23) alignment. The four public protein databases used include: (i) NCBI non-redundat (NR); (ii) Swiss-Prot; (iii) Kyoto Encyclopedia of Genes and Genomes (KEGG); and (iv) Cluster of Orthologous Groups (COG) at e-value < 0.00001. The direction of the identified unigenes was determined using the best alignments obtained from the four databases. In the case of conflicting results between different databases, the priority order: (i) NR; (ii) Swissprot; (iii) KEGG; and (iv) COG was used. BlastX data was used to extract the coding regions (CDS) from unigene sequences and translate them into peptide sequences. Unigenes not identified by BlastX were analyzed using ESTScan to predict their CDS and to decide their sequence direction while unigenes with NR annotation were further analyzed with Blast2go [34] to obtain their gene ontology (GO) annotations. The identified unigenes were classified according to GO functions using the Web Gene Ontology (WEGO) annotation software.

2.5. Identification of Differentially Expressed Genes

Mapped read counts for each gene generated from the functional annotation above were normalized for RNA length and total read counts in each lane using the reads per kilobase per million method (RPKM). As such, the RPKM method allowed for direct comparison of the number of transcripts between the SAV-3 infected and non-infected groups, which created the basis for identifying the differentially expressed genes (DEGs). We set the cutoff limit at 95% confidence interval for all RPKM values for each gene and used a rigorous algorithm to generate DEGs by comparing RPKM mapped reads from SAV-3 infected versus non-infected TO-cells. Only DEGs with a threshold of false discovery rate (FDR) <0.001 and an absolute value log2 ratio >1 were considered differentially expressed. Thereafter, all identified DEGs were mapped to GO annotations using the Blast2GO software [34] and were later assigned KEGG ortholog (KOs) identifiers for pathway analysis using the KEGG pathway analytical software using the zebrafish model.

2.6. Data Access

The RNA-sequencing data generated in this study have been deposited in the National Center for Biotechnology Information (NCBI) Gene Expression Omnibus (GEO) database accession number GSE64095 (www.ncbi.nih.gov/geo Accession number GSE64095) [35].

2.7. Validation of RNA-Seq Data and Virus Quantification

In order to confirm the validity of our RNA-seq data, 13 randomly selected DEGs shown to be up- or downregulated by RNA-seq were used for quantitative real-time PCR (qRT-PCR) analysis using the QuantiFast SYBR Green RT-PCR Kit (Qiagen) and the LightCycler 480 system (Roche). For each gene, the quantity of template, master mix final volume, reverse transcriptase, PCR initiation activation and cycles used per reaction were carried out as previously described [17]. Primer sequences used for RT-PCR are shown in Table 1. The specificity of each PCR product from each primer pair was confirmed by melting curve analysis and agarose gel analysis while the 2−∆∆Ct method was used to quantify the fold increase in gene expression levels relative to the control group. All quantifications were normalized using the β-actin endogenous gene, which has been shown to be a stable normalizer of different viral infections in Atlantic salmon in our studies [32,36,37]. For virus quantification, qRT-PCR was used to determine the quantity of virus in the SAV-3 infected and non-treated cells using the E2 SP expressed during virus replication using primer sequences used for E2 quantification are shown in Table 1 as previously described by Xu et al. [24].

3. Results

3.1. Gene Ontology Classification and KEGG Pathway Analysis

After filtration, a total of 20,115 unigenes were identified and assigned KOs identifiers. Analysis of DEGs using KOs resulted into 9315 genes being assigned to 252 pathways, which included the PRR pathways shown in Table 2. The significance of each pathway was set at p-value <0.05 while the cutoff for enrichment was set at Qvalue <0.50. The TLR pathway had the highest number of DEGs (112) followed by the NLR (89) and RLR pathways (79). Although the TLR pathway had a marginal significance (p = 0.058), it was more significant than the NLR pathway (p = 0.91), while the RLR pathway showed the highest significance (p = 0.024) of all the PRRs expressed in response to SAV-3 infection in TO-cells. As a result, the RLR pathway had the highest enrichment (Qvalue = 3.0117 × 10−1) followed by the TLR pathways (Qvalue = 4.8652 × 10−1) suggesting that both pathways played a pivotal role in the recognition of SAV-3 infection in TO-cells. On the contrary, the NLR pathway was not enriched (Qvalue = 1.0 × 100) indicating that NLRs had no capacity to recognize the invasion of SAV3 in TO-cells (Table 2).
Hence, in the next studies we focused on the RLRs and TLRs that were significantly enriched to identify the exact sensors within these PRRs that were able to recognize SAV3 infection in TO-cells. To identify the genes involved in downstream signaling after ligand binding of TLRs and RLRs to the SAV-3 PAMPs, we further analyzed the repertoire of DEGs expressed in each pathway as shown below.

3.2. Toll-Like Receptor Signaling Pathways

Table 3A shows the profile of genes upregulated in the TLR signaling pathway induced by SAV3 infection in TO-cells. Only the endosomal TLRs 3 and 8 were upregulated in which the fold increase for TLR8 was approximately threefold higher than TLR3. Among the IFN regulatory factors (IRFs), upregulation of IRF7 was at similar level as IRF3 (Table 3A). Other genes upregulated include IFN α/β receptor 1 (IFNAR1), IFN-a2, and the chemokines IFNγ induced protein 10 (IP-10) and IFN-inducible T-cell α chemoattractant (I-TAC). Figure 1 shows the KEGG network pathways in which only endosomal TLRs 3 and 8 were upregulated, as shown in Table 3A. Downstream signaling from the endosomal TLRs 3 and 8 show upregulation of IRF3 and IFR7 linked to upregulation of IFN-a2. In addition, Figure 1 also shows upregulation of the IFN-α/β receptor linked to upregulation of IP-10 and I-TAC.
To summarize the TLR pathways induced by SAV-3 infection in TO-cells, Figure 2 shows the TLR signaling pathway based on upregulated genes (Table 3A) excluding the downregulated genes (Table 3B).
Table 3B shows downregulated TLR pathway genes in TO-cells infected by SAV-3 in which the extracellular TLRs 1 and 2 were downregulated together with their downstream signaling genes belonging to the P13K-AKT and FADD-CASP8 pathways. In addition, Figure 1 shows that TLRs 4 and 5 were not expressed together with their downstream signaling genes like the translocating chain-associated membrane protein (TRAM), adaptor protein (TIRAP) and myeloid differentiation primary response protein (MyD88) in TO-cells infected by SAV-3. Put together, these data show that none of the extracellular TLRs were upregulated in response to SAV-3 infection in TO-cells.

3.3. RIG-I-Like Receptor Signaling Pathway

Table 4A shows RLR pathway genes that were upregulated in response to SAV-3 infection in TO-cells of which LPG2 had the highest expression followed by RIG-1 and MDA5, respectively. Downstream signaling showed upregulation of the endoplasmic reticulum mediator of IRF3 activation (MITA) and mitochondria IFNβ promoter stimulator I (IPS-I), which is also known as the mitochondrial antiviral-signaling protein (MAVS). Among the genes that regulate the expression of RIG-I, the tripartite motif-containing protein 25 (TRIM25) was upregulated when Ubiquitin carboxyl-terminal hydrolase CYLD was downregulated. Other upregulated genes included IRF3, IRF7, IP-10 and IFN-a2. Figure 3 shows genes differentially expressed for the RLR pathway induced by SAV-3 infection in TO-cells in which RIG-I and MDA5 were upregulated. In addition, Figure 3 shows that TRIM25 was linked to upregulation of RIG-I when CYLD was downregulated. Further, Figure 3 also shows that downstream signaling of RIG-I and MDA5 converge on the interferon-beta promoter stimulator 1 (IPS-I), which is linked to MITA found in the endoplasmic reticulum while downstream signaling from IPS-1 via IRF3 and IRF7 culminated in upregulation of IFN-a2. To summarize the RLR signaling pathway induced by SAV-3 infection in TO-cells, Figure 2 shows a pathway based on upregulated RLR genes (Table 4A) excluding the downregulated genes. Table 4B shows the RLR genes that were downregulated in TO-cells during SAV-3 infection. Consistent with the TLR pathways (Table 3B and Figure 1), genes involving the FADD-CASP8 signaling pathways were downregulated together with genes that signal via the TNF-receptor associated factor 2 (TRAF2) and 6 (TRAF6) pathways.

3.4. Quantitative Real-Time Polymerase Chain Reaction Test and Virus Quantification

Figure 4 shows randomly selected genes for qRT-PCR analysis of the TLR and RLR pathway genes detected by RNA-seq. In the situation where duplicated copies of the selected genes existed [13,38], only genes that had the highest significant (p-value) of expression were used for qRT-PCR analysis. Among the TLR genes, TLRs 3 and 8 were upregulated, while, among the RLR genes, RIG-I, MDA5, LPG2 and IPS-1 were also upregulated. In addition, IRF3, IRF7 and IP-10 reported from both TLR and RLR pathways were also upregulated. Genes that were downregulated included the toll-interacting protein (TOLLIP), TRAF6, p38b1 mitogen activated protein kinase (p38) and ras-related protein rac1 (Rac1). Overall, Figure 4 shows genes that genes that were up- or downregulated by RNA-seq were also up- or downregulated by qRT-PCR thereby confirming the validity of our RNA-seq data. In terms of virus quantification, SAV-3 was only detected in the infected cells at mean Cp value 20 (n = 3) of the E2 structural protein detected by qRT-PCR. E2 gene transcription was not detected by qRT-PCR in the non-infected cells.

4. Discussion

4.1. Transcriptome Signaling Pathway Analysis

Gene ontology (GO) and KEGG are among the most commonly used databases for functional annotation. While GO terms are mostly used for the annotation of individual genes [39], KEGG pathways are widely used for annotations in which genes can be grouped into network maps [39,40], which provides a functional understanding of genes that work together in a pathway. KEGG pathway maps can be plotted into biological pathways of model and none-model organisms [40,41]. Based on statistical analyses, pathways with significant enrichment can be determined to ensure that only pathways having relevant biological implications are used in the analysis of transcriptome data. Thus, pathways identified as the most enriched can be useful in identifying relevant genes activated in response to stimuli while genes that rank high in a pathway could serve as potential candidates for testing the validity of the pathway using functional studies such as the use of knock-out-models. Hence, in the case of duplicated genes, isoforms having the highest significance (p = value) of differential expression were used for the validation of RNA-seq data by qRT-PCR because they were considered to have the highest impact on influencing the outcome of the PRR network pathways induced by SAV-3 infection in TO-cells. In this study, the KEGG pathway analysis showed significant enrichment of the endosomal TLRs and RLRs and not the NLRs in TO-cells infected by SAV3. As such, further analysis paved the way to deciphering the pathway network of genes involved in TLR and RLR signaling expressed in response to SAV-3 infection in TO-cells as shown below. Detection of increased levels of the E2 structural protein in TO-cells 48 h post infection by qRT-PCR consolidates the fact that the TLR and RLR genes differentially expressed in this study were induced by SAV-3 infection in TO-cells. Hence, the pathway analysis carried out in this study suggests that SAV-3 infection could be using the TLR and RLR pathways to produce high levels of type I IFNs in TO-cells. Overall this study shows that pathway based analyses improves the analytical power to identify the most important genes expressed in response to stimuli in a de novo assembled transcriptome. Further, the study also shows that the use of pathway based analysis to decipher molecular networks of genes expressed in a transcriptome provides a contextual understanding of biological processes induced by microbial invasion unlike the tedious work of trying to identify the biological functions of individual genes expressed in a transcriptome using a non-systematic approach.

4.2. Toll-Like Receptor Signaling Genes

The TLR family is made of extracellular and intracellular receptors able to recognize PAMPs from different pathogens [42,43]. In the present study, only TLRs 3 and 8 were upregulated suggesting that recognition of SAV-3 in TO-cells could be by the endosomal TLRs and not the cell surface TLRs. Studies in mammals have shown that TLR8 is only expressed by phagocytic cells such as macrophages and dendritic cells (DCs) [44] while TLR3 has been shown to be a dsRNA innate immune receptor primarily expressed by macrophages and DCs [28]. Put together these observations support observations made by Pettersen et al. [19] who showed that TO-cells derived from Atlantic salmon leukocytes possess macrophage/dendritic cell like properties and, hence, their ability to express both TLRs 3 and 8 in response to SAV-3 infection further firms up this notion. In fish, TLR3 is expressed in high levels in different organs inclusive of mucosal and lymphoid organs suggesting that it could play an important role in sensing pathogens at portals of entry in mucosal organs as well as recognizing pathogens that reach the lymphoid organs after entering the systemic environment [45,46,47]. On the contrary, TLR8 is mainly expressed in lymphoid organs [17,38,48] suggesting that it plays a vital role in sensing viruses that get to lymphoid organs via APCs. Studies in humans have shown that TLR8 binds to viral ssRNA [26] suggesting that it could be using similar mechanisms to bind to the ssRNA genome of SAV-3 in TO-cells. On the other hand, TLR3 has been shown to be specific for the recognition of viral dsRNA [28]. This is supported by observations made by Weber et at [49] who showed that viral ssRNA produce intermediate dsRNA which is recognized by TLR3 and given that alphaviruses replicases form dsRNA intermediates as pointed by Smerdou et al. [50], it is likely that the sensing of SAV-3 infection by TLR3 in TO-cells could be by recognition of the intermediate dsRNA produced in its replicative form. However, there is need for detailed investigations to support these observations given that such information has proved to be useful in the targeting of alphavirus replicons in APCs for the optimization of vaccine performance as shown in the case of mammalian alphaviruses [51,52,53]. In terms of downstream signaling, upregulation of TLR3 and TLR8 was linked to upregulation of IRF3 and IRF7 culminating in upregulation of IFN-a2 which is in line with observations seen in higher vertebrates that TLR3 and TLR8 produce type I IFNs via IRF3 and IRF7 [54,55]. In addition, the study also shows that activation of IFN-α/β receptors was linked to upregulation of STAT1 via the JAK/STAT pathway resulting in upregulation of IP-10 and I-TAC, which are chemoattractants for T-cell responses in virus infected cells [56]. All in all, the repertoire of TLR genes expressed by TO-cells in response to SAV-3 infection conforms to genes expressed by the endosomal TLRs signaling pathways expressed in mammalian cells suggesting that fish macrophages and DC-like cells could be using similar mechanisms to those used by mammalian macrophages and DCs to combat intracellular microbial invasion [9,54,57,58]. However, there is need for detailed studies using knockout systems to demonstrate the functional mechanisms of individual genes expressed in the TLR signaling pathway shown in this study.

4.3. RIG-I-Like Receptor Signaling Genes

Although the significance of antiviral effects of TLRs 3 and 8 in macrophages and DCs is indisputable, the key viral sensors for other cell types for intracellular recognition of infection are RLRs [59]. In the case of ssRNA viruses, the major PAMP recognized by RLRs is the 5′-triphosphate (ppp-) RNA [60,61,62]. Single stranded 5′-ppp-RNAs that lack 2’-O-methylation of the 5′-cap, but bear a 5′-ppp-RNA, are specifically from viruses, which serve as a molecular signature for distinguishing self from non-self mRNAs [63,64]. This PAMP is recognized by RIG-I and it is expressed by several viruses including alphaviruses [56] suggesting that the sensing of SAV-3 by RIG-I in this study could be through the same PAMP used to bind to RIG-I in the cytosol for mammalian viruses [60,61,62]. Although RIG-I signaling is dependent on 5′-ppp-RNA binding, it requires ubiquitination by TRIM25 [65], oligomerization by MITA [60] and IPS-I multimerization on the mitochondrion-associated endoplasmic reticulum [66,67], which could account for the expression of these genes in this study. Apart from the 5-ppp-RNA PAMP, the presence of viral dsRNA in host cells is recognized as a non-self-entity given that vertebrate cells do not encode the RNA-dependent-RNA-polymerase (RdRp) encoded in RNA viral genomes [27]. Unlike RIG-I that detects the 5′-ppp-RNA PAMP, MDA5 functions as a sensor for dsRNA in the cytosol [29,30,68]. As pointed out by Nikonov et al. [27], the +ssRNA of alphaviruses serves as an mRNA that is transcribed to form a dsRNA in the cytosol during replication. Given that both RIG-I and MDA5 bind to dsRNA, it is likely that the dsRNA formed during virus replication serves as a ligand that binds to RIG-I and MDA5. Hence, upregulation of MDA5 in this study could be linked to detection of dsRNA generated from SAV-3 replication in TO-cells, which conforms to observations made for other alphaviruses [27]. Expression of MDA5 and RIG-I has been shown to increase following infection by dsRNA viruses such as infectious pancreatic necrosis virus (IPNV) in Atlantic salmon [69] suggesting that these PRRs could play an important role in virus recognition by APCs and lead to their activation to enhance antigen uptake, processing and presention for activation of the adaptive immune system [70]. On the other hand, LGP2 has been shown to potentiate the function of RIG-I while blocking the function of MDA5 in mammals [71,72]. The KEGG pathway analysis used in this study shows that the RIG-I and MDA5 pathways converge on the IPS-I adaptor, which is in line with observations seen in mammals [73,74]. Similar to observations made for TLR3 and TLR8 signaling pathways above, the RLR signaling pathway analysis carried out in this study shows that downstream signaling via IPS-1 culminate in production of type I IFNs using the IRF3 and IRF7 signaling pathways which is in line with observations seen in higher vertebrates [27,52,73]. Overall, the repertoire of genes clustered in the RLR signaling pathway generated in this study conforms to genes induced by other alphaviruses in higher vertebrate dendritic cells [75] suggesting that SAV-3 infection in TO-cells uses similar mechanisms to produce type I IFNs and anti-inflammatory cytokines used to combat alphavirus infection in higher vertebrates.

5. Conclusions

In this study, we have shown that the repertoire of genes linked to PRRs induced by SAV-3 infection in TO-cells is comparable to genes induced by other alphaviruses in Mammalia [66]. Among the TLRs, only endosomal TLRs 3 and 8 were upregulated while RIG-I, MDA5 and LPG2 were upregulated among the RLRs suggesting that these PRRs are essential for the sensing of SAV-3 infection in TO-cells. Both TLR and RLR signaling pathways were linked to upregulation of IRF3 and IRF7 culminating in upregulation of IFN-a2. This study links the expression TRIM25 to the RIG-I signaling pathway being the first report that points to the possible involvement of this gene in the recognition of SAV-3 infection in TO-cells. Finally, upregulation of IFN-a2 observed from the TLR and RLR pathways suggest that SAV-3 functions through these pathways for potent induction of IFN-a2. It is important to point out that data generated here is based on transcriptome analysis of different genes expressed in response to SAV-3 infection in TO-cells, there is need for further studies to consolidate these findings using functional studies such as gene knockout systems to elucidate the functional activities of all the genes expressed in this study. Overall, the study shows that a pathway-based approach improves the analytical power of transcriptome data analysis and that it provides a contextual approach to understanding the biological relevance of DEGs induced by microbial invasion.

Acknowledgments

This study was funded by the Research Council of Norway, project “The Atlantic salmon genome sequence as a tool for precision breeding” project number 226275, and the Targeted disease prophylaxis in European Fish farming project, grant agreement 311993. We thank the Beijing Genomics Institute (BGI) for assistance with some bioinformatics analysis.

Author Contributions

Cheng Xu: carried out laboratory experiments and data analysis; Hetron Mweemba Munang’andu: carried out data analysis and preparation of the manuscript; and Øystein Evensen: mobilizing of resources, data analysis, preparation of the manuscript and overall supervision of the project. All authors read and approved submission of the manuscript for publication.

Conflicts of Interest

Authors declare no conflict of interest.

References

  1. Arancibia, S.A.; Beltran, C.J.; Aguirre, I.M.; Silva, P.; Peralta, A.L.; Malinarich, F.; Hermoso, M.A. Toll-like receptors are key participants in innate immune responses. Biol. Res. 2007, 40, 97–112. [Google Scholar] [CrossRef] [PubMed]
  2. Munang’andu, H.M.; Evensen, Ø. A review of intra- and extracellular antigen delivery systems for virus vaccines of finfish. J. Immunol. Res. 2015, 2015. [Google Scholar] [CrossRef] [PubMed]
  3. Medzhitov, R.; PrestonHurlburt, P.; Janeway, C.A. A human homologue of the Drosophila Toll protein signals activation of adaptive immunity. Nature 1997, 388, 394–397. [Google Scholar] [PubMed]
  4. Matzinger, P. The danger model: A renewed sense of self. Science 2002, 296, 301–305. [Google Scholar] [CrossRef] [PubMed]
  5. Brunette, R.L.; Young, J.M.; Whitley, D.G.; Brodsky, I.E.; Malik, H.S.; Stetson, D.B. Extensive evolutionary and functional diversity among mammalian AIM2-like receptors. J. Exp. Med. 2012, 209, 1969–1983. [Google Scholar] [CrossRef] [PubMed]
  6. Germain, R.N. An innately interesting decade of research in immunology. Nat. Med. 2004, 10, 1307–1320. [Google Scholar] [CrossRef] [PubMed]
  7. Hornung, V.; Ablasser, A.; Charrel-Dennis, M.; Bauernfeind, F.; Horvath, G.; Caffrey, D.R.; Latz, E.; Fitzgerald, K.A. AIM2 recognizes cytosolic dsDNA and forms a caspase-1-activating inflammasome with ASC. Nature 2009, 458, 514–518. [Google Scholar] [CrossRef] [PubMed]
  8. Medzhitov, R.; Shevach, E.M.; Trinchieri, G.; Mellor, A.L.; Munn, D.H.; Gordon, S.; Libby, P.; Hansson, G.K.; Shortman, K.; Dong, C.; et al. Highlights of 10 years of immunology in Nature Reviews Immunology. Nat. Rev. Immunol. 2011, 11, 693–702. [Google Scholar] [CrossRef] [PubMed]
  9. Akira, S.; Takeda, K. Toll-like receptor signalling. Nat. Rev. Immunol. 2004, 4, 499–511. [Google Scholar] [CrossRef] [PubMed]
  10. Chtarbanova, S.; Imler, J.L. Microbial sensing by Toll receptors: A historical perspective. Arterioscler. Thromb. Vasc. Biol. 2011, 31, 1734–1738. [Google Scholar] [CrossRef] [PubMed]
  11. Creagh, E.M.; O’Neill, L.A.J. TLRs, NLRs and RLRs: A trinity of pathogen sensors that co-operate in innate immunity. Trends Immunol. 2006, 27, 352–357. [Google Scholar] [CrossRef] [PubMed]
  12. Palti, Y. Toll-like receptors in bony fish: From genomics to function. Dev. Comp. Immunol. 2011, 35, 1263–1272. [Google Scholar] [CrossRef] [PubMed]
  13. Rebl, A.; Goldammer, T.; Seyfert, H.M. Toll-like receptor signaling in bony fish. Vet. Immunol. Immunopathol. 2010, 134, 139–150. [Google Scholar] [CrossRef] [PubMed]
  14. Laing, K.J.; Purcell, M.K.; Winton, J.R.; Hansen, J.D. A genomic view of the NOD-like receptor family in teleost fish: Identification of a novel NLR subfamily in zebrafish. Bmc Evol. Biol. 2008, 8. [Google Scholar] [CrossRef] [PubMed]
  15. Chang, M.X.; Collet, B.; Nie, P.; Lester, K.; Campbell, S.; Secombes, C.J.; Zou, J. Expression and functional characterization of the RIG-I-Like receptors MDA5 and LGP2 in rainbow trout (Oncorhynchus mykiss). J. Virol. 2011, 85, 8403–8412. [Google Scholar] [CrossRef] [PubMed]
  16. Iliev, D.B.; Sobhkhez, M.; Fremmerlid, K.; Jorgensen, J.B. MyD88 interacts with interferon regulatory factor (IRF) 3 and IRF7 in Atlantic Salmon (Salmo salar) transgenic SsMyD88 modulates the IRF-Induced type I interferon response and accumulates in aggresomes. J. Biol. Chem. 2011, 286, 42715–42724. [Google Scholar] [CrossRef] [PubMed]
  17. Skjaeveland, I.; Iliev, D.B.; Strandskog, G.; Jorgensen, J.B. Identification and characterization of TLR8 and MyD88 homologs in Atlantic salmon (Salmo salar). Dev. Comp. Immunol. 2009, 33, 1011–1017. [Google Scholar] [CrossRef] [PubMed]
  18. Xu, C.; Evensen, O.; Munang’andu, H.M. De novo assembly and transcriptome analysis of Atlantic salmon macrophage/dendritic-like TO cells following type I IFN treatment and Salmonid alphavirus subtype-3 infection. BMC Genom. 2015, 16. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  19. Pettersen, E.F.; Ingerslev, H.C.; Stavang, V.; Egenberg, M.; Wergeland, H.I. A highly phagocytic cell line TO from Atlantic salmon is CD83 positive and M-CSFR negative, indicating a dendritic-like cell type. Fish Shellfish Immunol. 2008, 25, 809–819. [Google Scholar] [CrossRef] [PubMed]
  20. Wergeland, H.I.; Jakobsen, R.A. A salmonid cell line (TO) for production of infectious salmon anaemia virus (ISAV). Dis. Aquat. Org. 2001, 44, 183–190. [Google Scholar] [CrossRef] [PubMed]
  21. Munro, A.L.S.; Ellis, A.E.; Mcvicar, A.H.; Mclay, H.A.; Needham, E.A. An exocrine pancreas disease of farmed atlantic salmon in Scotland. Helgol. Meeresunters. 1984, 37, 571–586. [Google Scholar]
  22. Murphy, T.M.; Rodger, H.D.; Drinan, E.M.; Gannon, F.; Kruse, P.; Korting, W. The sequential pathology of pancreas disease in atlantic salmon farms in Ireland. J. Fish Dis. 1992, 15, 401–408. [Google Scholar] [CrossRef]
  23. Weston, J.H.; Welsh, M.D.; McLoughlin, M.F.; Todd, D. Salmon pancreas disease virus, an alphavirus infecting farmed Atlantic salmon, Salmo salar L. Virology 1999, 256, 188–195. [Google Scholar] [CrossRef] [PubMed]
  24. Rheme, C.; Ehrengruber, M.U.; Grandgirard, D. Alphaviral cytotoxicity and its implication in vector development. Exp. Physiol. 2005, 90, 45–52. [Google Scholar] [CrossRef] [PubMed]
  25. Strauss, J.H.; Strauss, E.G. The alphaviruses: Gene expression, replication, and evolution. Microbiol. Rev. 1994, 58, 491–562. [Google Scholar] [PubMed]
  26. Bernard, M.A.; Han, X.B.; Inderbitzin, S.; Agbim, I.; Zhao, H.; Koziel, H.; Tachado, S.D. HIV-derived ssRNA binds to TLR8 to induce inflammation-driven macrophage foam cell formation. PLoS ONE 2014, 9. [Google Scholar] [CrossRef] [PubMed]
  27. Nikonov, A.; Molder, T.; Sikut, R.; Kiiver, K.; Mannik, A.; Toots, U.; Lulla, A.; Lulla, V.; Utt, A.; Merits, A.; et al. RIG-I and MDA-5 detection of viral RNA-dependent RNA polymerase activity restricts positive-strand RNA virus replication. PLoS Pathog. 2013, 9. [Google Scholar] [CrossRef] [PubMed]
  28. Alexopoulou, L.; Holt, A.C.; Medzhitov, R.; Flavell, R.A. Recognition of double-stranded RNA and activation of NF-κB by Toll-like receptor 3. Nature 2001, 413, 732–738. [Google Scholar] [CrossRef] [PubMed]
  29. Feng, Q.; Hato, S.V.; Langereis, M.A.; Zoll, J.; Virgen-Slane, R.; Peisley, A.; Hur, S.; Semler, B.L.; van Rij, R.P.; van Kuppeveld, F.J.M. MDA5 detects the double-stranded RNA replicative form in picornavirus-infected cells. Cell Rep. 2012, 2, 1187–1196. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  30. Triantafilou, K.; Vakakis, E.; Kar, S.; Richer, E.; Evans, G.L.; Triantafilou, M. Visualisation of direct interaction of MDA5 and the dsRNA replicative intermediate form of positive strand RNA viruses. J. Cell Sci. 2012, 125 Pt 20, 4761–4769. [Google Scholar] [CrossRef] [PubMed]
  31. Knipe, D.M.; Howley, P.M. Fundamental Virology, 4th ed.; Lippincott Williams & Wilkins: New York, NY, USA, 2001. [Google Scholar]
  32. Xu, C.; Guo, T.C.; Mutoloki, S.; Haugland, O.; Marjara, I.S.; Evensen, O. α interferon and not γ interferon inhibits salmonid alphavirus subtype 3 replication In Vitro. J. Virol. 2010, 84, 8903–8912. [Google Scholar] [CrossRef] [PubMed]
  33. Hua, W.P.; Zhang, Y.; Song, J.; Zhao, L.J.; Wang, Z.Z. De novo transcriptome sequencing in Salvia miltiorrhiza to identify genes involved in the biosynthesis of active ingredients. Genomics 2011, 98, 272–279. [Google Scholar]
  34. Conesa, M.; Götz, S.; García-Gómez, J.M.; Terol, J.; Talón, M.; Robles, M. Blast2GO: a universal tool for annotation, visualization and analysis in functional genomics research. Bioinformatics 2005, 21, 3674–3676. [Google Scholar]
  35. GEO Accession viewer – NCBI. Available online: http://0-www-ncbi-nlm-nih-gov.brum.beds.ac.uk/geo/query/acc.cgi?acc=GSE64095 (accessed on 18 April 2016).
  36. Munang’andu, H.M.; Fredriksen, B.N.; Mutoloki, S.; Dalmo, R.A.; Evensen, O. The kinetics of CD4+ and CD8+ T-cell gene expression correlate with protection in Atlantic salmon (Salmo salar L) vaccinated against infectious pancreatic necrosis. Vaccine 2013, 31, 1956–1963. [Google Scholar] [CrossRef] [PubMed]
  37. Xu, C.; Guo, T.C.; Mutoloki, S.; Haugland, O.; Evensen, O. Gene expression studies of host response to Salmonid alphavirus subtype 3 experimental infections in Atlantic salmon. Vet. Res. 2012, 43, 78. [Google Scholar] [CrossRef] [PubMed]
  38. Lee, P.T.; Zou, J.; Holland, J.W.; Martin, S.A.; Kanellos, T.; Secombes, C.J. Identification and characterization of TLR7, TLR8a2, TLR8b1 and TLR8b2 genes in Atlantic salmon (Salmo salar). Dev. Comp. Immunol. 2013, 41, 295–305. [Google Scholar] [CrossRef] [PubMed]
  39. Lomax, J. Get ready to GO! A biologist’s guide to the Gene Ontology. Brief. Bioinform. 2005, 6, 298–304. [Google Scholar] [CrossRef] [PubMed]
  40. Li, C.; Liu, D.R.; Li, G.G.; Wang, H.H.; Li, X.W.; Zhang, W.; Wu, Y.L.; Chen, L. CD97 promotes gastric cancer cell proliferation and invasion through exosome-mediated MAPK signaling pathway. World J. Gastroenterol. 2015, 21, 6215–6228. [Google Scholar] [CrossRef] [PubMed]
  41. Liu, X.; Wang, J.; Sun, G. Identification of key genes and pathways in renal cell carcinoma through expression profiling data. Kidney Blood Press Res. 2015, 40, 288–297. [Google Scholar] [CrossRef] [PubMed]
  42. Kawai, T.; Akira, S. TLR signaling. Cell Death Differen. 2006, 13, 816–825. [Google Scholar] [CrossRef] [PubMed]
  43. Kawai, T.; Akira, S. The role of pattern-recognition receptors in innate immunity: Update on Toll-like receptors. Nat. Immunol. 2010, 11, 373–384. [Google Scholar] [CrossRef] [PubMed]
  44. Schroder, M.; Bowie, A.G. TLR3 in antiviral immunity: Key player or bystander? Trends Immunol. 2005, 26, 462–468. [Google Scholar] [CrossRef] [PubMed]
  45. Rodriguez, M.F.; Wiens, G.D.; Purcell, M.K.; Palti, Y. Characterization of Toll-like receptor 3 gene in rainbow trout (Oncorhynchus mykiss). Immunogenetics 2005, 57, 510–519. [Google Scholar] [CrossRef] [PubMed]
  46. Abid, A.; Davies, S.J.; Waines, P.; Emery, M.; Castex, M.; Gioacchini, G.; Carnevali, O.; Bickerdike, R.; Romero, J.; Merrifield, D.L. Dietary synbiotic application modulates Atlantic salmon (Salmo salar) intestinal microbial communities and intestinal immunity. Fish Shellfish Immunol. 2013, 35, 1948–1956. [Google Scholar] [CrossRef] [PubMed]
  47. Munang’andu, H.M.; Mutoloki, S.; Evensen, O. A review of the immunological mechanisms following mucosal vaccination of finfish. Front. Immunol. 2015, 6, 427. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  48. Palti, Y.; Rodriguez, M.F.; Gahr, S.A.; Purcell, M.K.; Rexroad, C.E., III; Wiens, G.D. Identification, characterization and genetic mapping of TLR1 loci in rainbow trout (Oncorhynchus mykiss). Fish Shellfish Immunol. 2010, 28, 918–926. [Google Scholar] [CrossRef] [PubMed]
  49. Weber, F.; Wagner, V.; Rasmussen, S.B.; Hartmann, R.; Paludan, S.R. Double-stranded RNA is produced by positive-strand RNA viruses and DNA viruses but not in detectable amounts by negative-strand RNA viruses. J. Virol. 2006, 80, 5059–5064. [Google Scholar] [CrossRef] [PubMed]
  50. Smerdou, C.; Liljestrom, P. Non-viral amplification systems for gene transfer: Vectors based on alphaviruses. Curr. Opin. Mol. Ther. 1999, 1, 244–251. [Google Scholar] [PubMed]
  51. Knudsen, M.L.; Johansson, D.X.; Kostic, L.; Nordstron, E.K.L.; Tegerstedt, K.; Pasetto, A.; Applequist, S.E.; Ljungberg, K.; Sirard, J.C.; Liljestrom, P. The adjuvant activity of alphavirus replicons is enhanced by incorporating the microbial molecule flagellin into the replicon. PLoS ONE 2013, 8, e65964. [Google Scholar] [CrossRef] [PubMed]
  52. Rayner, J.O.; Dryga, S.A.; Kamrud, K.I. Alphavirus vectors and vaccination. Rev. Med. Virol. 2002, 12, 279–296. [Google Scholar] [CrossRef] [PubMed]
  53. Vander Veen, R.L.; Harris, D.L.; Kamrud, K.I. Alphavirus replicon vaccines. Anim. Health Res. Rev. 2012, 13, 1–9. [Google Scholar] [CrossRef] [PubMed]
  54. Cervantes, J.L.; Weinerman, B.; Basole, C.; Salazar, J.C. TLR8: The forgotten relative revindicated. Cell. Mol. Immunol. 2012, 9, 434–438. [Google Scholar] [CrossRef] [PubMed]
  55. Jensen, S.; Thomsen, A.R. Sensing of RNA viruses: A review of innate immune receptors involved in recognizing RNA virus invasion. J. Virol. 2012, 86, 2900–2910. [Google Scholar] [CrossRef] [PubMed]
  56. Renn, C.N.; Sanchez, D.J.; Ochoa, M.T.; Legaspi, A.J.; Oh, C.K.; Liu, P.T.; Krutzik, S.R.; Sieling, P.A.; Cheng, G.H.; Modlin, R.L. TLR activation of langerhans cell-like dendritic cells triggers an antiviral immune response. J. Immunol. 2006, 177, 298–305. [Google Scholar] [CrossRef] [PubMed]
  57. Blasius, A.L.; Beutler, B. Intracellular Toll-like receptors. Immunity 2010, 32, 305–315. [Google Scholar] [CrossRef] [PubMed]
  58. Takeda, K.; Akira, S. TLR signaling pathways. Semin. Immunol. 2004, 16, 3–9. [Google Scholar] [CrossRef] [PubMed]
  59. Kato, H.; Sato, S.; Yoneyama, M.; Yamamoto, M.; Uematsu, S.; Matsui, K.; Tsujimura, T.; Takeda, K.; Fujita, T.; Takeuchi, O.; et al. Cell type-specific involvement of RIG-I in antiviral response. Immunity 2005, 23, 19–28. [Google Scholar] [CrossRef] [PubMed]
  60. Gack, M.U. Mechanisms of RIG-I-Like receptor activation and manipulation by viral pathogens. J. Virol. 2014, 88, 5213–5216. [Google Scholar] [CrossRef] [PubMed]
  61. Leung, D.W.; Amarasinghe, G.K. Structural insights into RNA recognition and activation of RIG-I-like receptors. Curr. Opin. Struct. Biol. 2012, 22, 297–303. [Google Scholar] [CrossRef] [PubMed]
  62. Ramos, H.J.; Gale, M. RIG-I like receptors and their signaling crosstalk in the regulation of antiviral immunity. Curr. Opin. Virol. 2011, 1, 167–176. [Google Scholar] [CrossRef] [PubMed]
  63. Abbas, Y.M.; Pichmair, A.; Gorna, M.W.; Superti-Furga, G.; Nagar, B. Structural basis for viral 5′-ppp-RNA recognition by human IFIT proteins. Nature 2013, 494, 60–64. [Google Scholar] [CrossRef] [PubMed]
  64. Pichlmair, A.; Lassnig, C.; Eberle, C.A.; Gorna, M.W.; Baumann, C.L.; Burkard, T.R.; Burckstummer, T.; Stefanovic, A.; Krieger, S.; Bennett, K.L.; et al. IFIT1 is an antiviral protein that recognizes 5′-triphosphate RNA. Nat. Immunol. 2011, 12, 624–630. [Google Scholar] [CrossRef] [PubMed]
  65. Gack, M.U.; Shin, Y.C.; Joo, C.H.; Urano, T.; Liang, C.; Sun, L.J.; Takeuchi, O.; Akira, S.; Chen, Z.J.; Inoue, S.S.; et al. TRIM25 RING-finger E3 ubiquitin ligase is essential for RIG-I-mediated antiviral activity. Nature 2007, 446, 916–919. [Google Scholar] [CrossRef] [PubMed]
  66. Horner, S.M.; Liu, H.M.; Park, H.S.; Briley, J.; Gale, M. Mitochondrial-associated ER membranes form MAVS-anchored innate immune synapses that are targeted by hepatitis C virus. Cytokine 2011, 56, 65. [Google Scholar] [CrossRef]
  67. Hou, F.J.; Sun, L.J.; Zheng, H.; Skaug, B.; Jiang, Q.X.; Chen, Z.J.J. MAVS forms functional prion-like aggregates to activate and propagate antiviral innate immune response. Cell 2011, 146, 448–461. [Google Scholar] [CrossRef] [PubMed]
  68. Peisley, A.; Jo, M.H.; Lin, C.; Wu, B.; Orme-Johnson, M.; Walz, T.; Hohng, S.; Hur, S. Kinetic mechanism for viral dsRNA length discrimination by MDA5 filaments. Proc. Natl. Acad. Sci. USA 2012, 109, E3340–E3349. [Google Scholar] [CrossRef] [PubMed]
  69. Skjesol, A.; Skjaeveland, I.; Elnaes, M.; Timmerhaus, G.; Fredriksen, B.N.; Jorgensen, S.M.; Krasnov, A.; Jorgensen, J.B. IPNV with high and low virulence: Host immune responses and viral mutations during infection. Virol. J. 2011, 8, 396. [Google Scholar] [CrossRef] [PubMed]
  70. Munang’andu, H.M.; Fredriksen, B.N.; Mutoloki, S.; Dalmo, R.A.; Evensen, O. Antigen dose and humoral immune response correspond with protection for inactivated infectious pancreatic necrosis virus vaccines in Atlantic salmon (Salmo salar L). Vet. Res. 2013, 44, 7. [Google Scholar] [CrossRef] [PubMed]
  71. Childs, K.S.; Randall, R.E.; Goodbourn, S. LGP2 plays a critical role in sensitizing mda-5 to activation by double-stranded RNA. PLoS ONE 2013, 8. [Google Scholar] [CrossRef] [PubMed]
  72. Rodriguez, K.R.; Bruns, A.M.; Horvath, C.M. MDA5 and LGP2: Accomplices and antagonists of antiviral signal transduction. J. Virol. 2014, 88, 8194–8200. [Google Scholar] [CrossRef] [PubMed]
  73. Kawai, T.; Takahashi, K.; Sato, S.; Coban, C.; Kumar, H.; Kato, H.; Ishii, K.J.; Takeuchi, O.; Akira, S. IPS-1, an adaptor triggering RIG-I- and Mda5-mediated type I interferon induction. Nat. Immunol. 2005, 6, 981–988. [Google Scholar] [CrossRef] [PubMed]
  74. Loo, Y.M.; Fornek, J.; Crochet, N.; Bajwa, G.; Perwitasari, O.; Martinez-Sobrido, L.; Akira, S.; Gill, M.A.; Garcia-Sastre, A.; Katze, M.G.; et al. Distinct RIG-I and MDA5 signaling by RNA viruses in innate immunity. J. Virol. 2008, 82, 335–345. [Google Scholar] [CrossRef] [PubMed]
  75. Hidmark, A.S.; McInerney, G.M.; Nordstrom, E.K.L.; Douagi, I.; Werner, K.M.; Liljestrom, P.; Hedestam, G.B.K. Early α/β interferon production by myeloid dendritic cells in response to UV-inactivated virus requires viral entry and interferon regulatory factor 3 but not MyD88. J. Virol. 2005, 79, 10376–10385. [Google Scholar] [CrossRef] [PubMed]
Figure 1. The KEGG pathway analysis for the Toll like receptor (TLR) signaling pathway differentially expressed genes (DEGs) expressed in response to SAV-3 infection in TO-cells. Red squares show unpregulated genes while green squares represent downregulated genes. Square having both red and green represent a mixed expression of upregulated (red) and downregulated (green) unigenes for the gene represented. Black squares show that the genes represented were not expressed in TO-cells.
Figure 1. The KEGG pathway analysis for the Toll like receptor (TLR) signaling pathway differentially expressed genes (DEGs) expressed in response to SAV-3 infection in TO-cells. Red squares show unpregulated genes while green squares represent downregulated genes. Square having both red and green represent a mixed expression of upregulated (red) and downregulated (green) unigenes for the gene represented. Black squares show that the genes represented were not expressed in TO-cells.
Viruses 08 00114 g001
Figure 2. A summary of the Toll like receptors (TLR) and RIG-I like receptor (RLR) pathway genes based on upregulated genes shown in Table 3A and Table 4A excluding the downregulated genes shown in Table 3B and Table 4B. Pathway A shows the RLR signaling pathway, while Pathway B shows the TLR signaling pathway.
Figure 2. A summary of the Toll like receptors (TLR) and RIG-I like receptor (RLR) pathway genes based on upregulated genes shown in Table 3A and Table 4A excluding the downregulated genes shown in Table 3B and Table 4B. Pathway A shows the RLR signaling pathway, while Pathway B shows the TLR signaling pathway.
Viruses 08 00114 g002
Figure 3. The KEGG pathway analysis of the RIG-I-like receptor DEGs expressed in response to SAV-3 infection TO-cells. Red squares show upregulated genes while green squares represent downregulated genes. Square having both red and green represent a mixed expression of upregulated (red) and downregulated (green) unigenes for the gene represented. Black squares show that the genes represented were not expressed in TO-cells.
Figure 3. The KEGG pathway analysis of the RIG-I-like receptor DEGs expressed in response to SAV-3 infection TO-cells. Red squares show upregulated genes while green squares represent downregulated genes. Square having both red and green represent a mixed expression of upregulated (red) and downregulated (green) unigenes for the gene represented. Black squares show that the genes represented were not expressed in TO-cells.
Viruses 08 00114 g003
Figure 4. Shows the Pearson’s correlation coefficency (r2 = 0.813, p = 0.0023) of RNA-seq versus qRT-PCR. Dots show number of RNA-Seq versus qRT-PCR pairs analyzed while the dotted lines show the 95% Confidence Interval of the correlation coefficiency line.
Figure 4. Shows the Pearson’s correlation coefficency (r2 = 0.813, p = 0.0023) of RNA-seq versus qRT-PCR. Dots show number of RNA-Seq versus qRT-PCR pairs analyzed while the dotted lines show the 95% Confidence Interval of the correlation coefficiency line.
Viruses 08 00114 g004
Table 1. Primers used for quantitative real time PCR.
Table 1. Primers used for quantitative real time PCR.
Primer NameSequenceGeneBank Accession No.
SAV-3 E2-FCAGTGAAATTCGATAAGAAGTGCAAEF675594
SAV-3 E2-RTGGGAGTCGCTGGTAAAGGT
β-Actin-FCCAGTCCTGCTCACTGAGGCAF012125
β-Actin-RGGTCTCAAACATGATCTGGGTCA
IP-10-FTGCCAGAACATGGAGATCATEF619047
IP-10-RTTTACTGCACACTCCTTTGGTT
TLR3-FTTTGATGAGTCTCCGCCAACTCCAKP231342
TLR3-RAATCTGCGAGGGACACAAAGGTCT
TLR8-FACAAGAAAGAATGCCTCAATGTCANM_001161693
TLR8-RCACCCAGTCTGACACCAACA
IRF3-FTGGACCAATCAGGAGCGAACFJ517643
IRF3-RAGCCCACGCCTTGAAAATAA
IRF7-FGAGGAGTGGGCAGAGAACTANM_001171850
IRF7-RTTCTGGGAGACTGGCTGGG
STAT1-FCGGGCCCTGTCACTGTTCGQ325309
STAT1-RGGCATACAGGGCTGTCTCT
RIG-I-FGACGGTCAGCAGGGTGTACTNM_001163699
RIG-I-RCCCGTGTCCTAACGAACAGT
MDA5-FAGAGCCCGTCCAAAGTGAAGTNM_001195179
MDA5-RGTTCAGCATAGTCAAAGGCAGGTA
LGP2-FGTGGCAGGCAATGGGGAATGFN396358
LGP2-RCCTCCAGTGTAATAGCGTATCAATCC
TOLLIP-FACCATTAGCACCCAACGAGBT045489
TOLLIP-RTGGGAGTAATACGCAGGAAG
RAC1-FGACAGGAAGACTACGACAGACNM_001160673
RAC1-RTCAAAGGAGGCAGGACTCAC
TRAF6-FACAGACTGTCCAAAGGCTC
TRAF6-RTCATTGCGCTGCATCATC
P38-FTCCACGCCAAGAGAACCTACNM_001123715
P38-RACATCATTGAACTCCTCCAGAC
Table 2. KEGG analysis of the pattern recognition receptors induced by SAV-3 infection in TO-cells.
Table 2. KEGG analysis of the pattern recognition receptors induced by SAV-3 infection in TO-cells.
ParametersToll Like ReceptorRIG-I-Like ReceptorNOD Like Receptor
Pathway IDKo04620Ko04622ko04621
Pathway significance0.0580.0240.9101
Pathway enrichment4.865157 × 10−13.011733 × 10−11.00000 × 100
Total KO genes201152011520115
All genes with pathway annotation931593159315
All genes in each pathway216144212
DEGs1127989
Table 3A. Toll like receptor pathway genes upregulated in response to SAV-3 infection in TO-cells.
Table 3A. Toll like receptor pathway genes upregulated in response to SAV-3 infection in TO-cells.
Gene NameAbbr.NCBIUnigKORegLog2 ratiop-Value
Toll like receptor 3TLR3|DAA64469.1|Unig9113K05401Up2.61407.4290 × 10−71
Toll like receptor 8TLR8|NP_001155165.1|Unig2363K10170Up4,04623.4201 × 10−5
Signal transducer and activator of transcription 1STAT1|NP_001134757.1|CL2066.2K11220Up6.272133.0554 × 10−68
Interferon regulatory factor 3IRF3|ACL68544.1|Unig4271K05411Up3.36445.3137 × 10−135
Interferon regulatory factor 7IRF7|NP_001165321.1|Unig10251K09447Up3.19701,13523 × 10−22
Interferon αIFN-a2|NP_001117042.1|Unig5589K05414Up7.6042
Interferon α receptor 1IFNAR1|NP_001268239.1|Unig34816K05130Up1.86408.9299 × 10−66
IFNγ induced protein 10IP-10|ACI69209.1|Unig8163K12671Up7.52332.2267 × 10−112
IFN-inducible T-cell α chemoattractantI-TAC|NM_0011412293.1|Unig1740K12762Up9,55672
Table 3B. Toll like receptor pathway genes downregulated during SAV-3 infection in TO-cells.
Table 3B. Toll like receptor pathway genes downregulated during SAV-3 infection in TO-cells.
Gene NameAbbr.NCBIUnigeneKORegLog2 ratiop-Value
Receptor interacting serine/threonine protein kinase 1RIP1|NP_001036815.1|Unig17924K02861Down−1.30562.0987 × 10−8
Caspase 8CASP8|XP_001335163|CL4461.1K04398Down−1.18002.4337 × 10−142
Toll like receptor 1TLR1|ACV92064.1|Unig41380K05398Down−4.25143.3298 × 10−5
Toll like receptor 2TLR2|CCK73195.1|Unig9045K10159Down−165891.7562 × 10−41
Transcription factor AP-1AP-1|XP_004369047.1|CL3191.1K04448Down−2.20253.9706 × 10−9
Extracellular signal-regulated kinaseERK|BAD23843.1|Unig24550K04371Down−1.80161.5449 × 10−5
NF-kappa-B inhibitor αNFκBα|ACI67986.1|CL8473.1K04735Down−1.39231.6765 × 10−13
TANK-binding kinase 1TBK1|JF241943.1|Unig5544K05410Down−1.26191.0212 × 10−124
TNF receptor associated factor 6TRAF6Unig40008K03175Down−3.15833.552 × 10−4
Interleukin 8IL-8|NP_001134182.1|Unig7278K10030Down−1.73686.5401 × 10−93
Kinase 1-binding protein 1TAB1|XP_002662286.2|Unig1972K04403Down−1.56147.0199 × 10−66
Kinase 1-binding protein 2TAB2|XP_003971436.1|CL4395K04404Down−2.079259.2399 × 10−28
Phosphatidylinositol-4,5-bisphosphate 3-kinasePI3K|XP_003455769.1|CL120K02649Down−1,92955,7289 × 10−39
RAC-α serine/threonine-protein kinase (AkT)AkT|ACH70834.1|CL5806K04456Down−1.434782.9777 × 10−19
Mitogen-activated protein kinase kinase 6 MKK6|AAV52830|Unig80K04433Down−1.65907.5296 × 10−169
Mitogen-activated protein kinase kinase 4MKK4|ACI33552.1|CL292.2KO4430Down−1.803971.9570 × 10−19
p38b1 mitogen activated protein kinasep38|EF123660.1|Unig10574K04441Down−1.61427.29346 × 10−9
Table 4A. RIG-I-like receptor pathway genes upregulated in response to SAV-3 infection in TO-cells.
Table 4A. RIG-I-like receptor pathway genes upregulated in response to SAV-3 infection in TO-cells.
Gene NameAbbr.NCBIUnigeneKORegLog2ratiop-Value
Retinoic acid-inducible gene-IRIG-I|NP_001157171.1|Unig7848K12646Up3.93773.171 × 10−102
Melanoma differentiation associated gene 5MDA5|NP_001182108.1|Unig6816K12647Up1.77881.8065 × 10−172
Laboratory of genetics and physiology 2LPG2 |NP_001133649.1|CL8555K12649Up5.51280
Interferon promoter stimulating protein 1IPS-1|NP_001161824.1|Uni12389K12648Up1.55883.0551 × 10−18
Tripartite motif-containing protein 25TRIM25|ACN11060.1|CL8518.2K10652UP4.20801.8646 × 10−5
IFNγ induced protein 10IP-10|ACI69209.1|Unig8163K12671Up7.52332.2267 × 10−112
OptineurinOptn|NP_001133761.1|CL4866K07210Up1.98530
Interferon regulatory factor 3IRF3|ACL68544.1|Unig4271K05411Up3.36445.3137 × 10−135
Interferon regulatory factor 7IRF7|NP_001165321.1|Unig8533K09447Up11.0551.1354 × 10−22
Interferon a2IFN-a2|NP_001117042.1|Unig5589K05414Up7.60420
Table 4B. RIG-I-like receptor pathway genes downregulated during SAV-3 infection in TO-cells.
Table 4B. RIG-I-like receptor pathway genes downregulated during SAV-3 infection in TO-cells.
Gene NameAbbr.NCBIUnigeneKORegLog2ratiop-Value
NLR family member X1NLRX1|AFY26970.1|Unig11078K12653Down−1.65001.5422 × 10−58
Autophagy protein 5Atg5|ACN11274.1|Unig5014K08339Down−1.37487.961 × 10−10
Interleukin 8IL-8|ABA86669.1|Unig7278K10030Down−1.73686.5401 × 10−93
Ubiquitin carboxyl-terminal hydrolase CYLDCyLD|XP_004068277.1|Unig18460K08601Down−1.98264.4252 × 10−11
Suppressor of IKK-epsilonSIKE|ACI33887.1|Unig21218K12656Down−1.53442.5131 × 10−7
TNF receptor type 1-associated death domain TRADD|Q1M161.1|CL8304K03171Down−1.66651.8945 × 10−20
TGF-β-activated kinase 1TAK1 |AAT07829.1|CL7020.1K04427Down−1.16433.2707 × 10−14
Nuclear factor κ-B NFκβ|HM771267|CL8473.1K04735Down−1.39231.6765 × 10−13
TANK-binding kinase 1 TBK1|JF241943.1|Unig5544K05410Down−1.26191.0212 × 10−124
TNF receptor-associated factor 2TRAF2|NP_001167255.1|Unig16762K03173Down−1.14451.0105 × 10−4
NF-κ-B inhibitorIκB|CAC85086.1|Unig31637K02581Down−2.83691.9042 × 10−4
Caspase 8CASP8|XP_001335163|CL4461.1K04398Down−1.18002.4337 × 10−142
Caspase 10CASP10|CAE51933.1|CL7349K04400Down−1.22522.8250 × 10−23
Receptor-interacting threonine-protein kinase 1RIPK1|NP_001036815.1|Unig17924K02861Down−1.30562.0987 × 10−8

Share and Cite

MDPI and ACS Style

Xu, C.; Evensen, Ø.; Munang’andu, H.M. De Novo Transcriptome Analysis Shows That SAV-3 Infection Upregulates Pattern Recognition Receptors of the Endosomal Toll-Like and RIG-I-Like Receptor Signaling Pathways in Macrophage/Dendritic Like TO-Cells. Viruses 2016, 8, 114. https://0-doi-org.brum.beds.ac.uk/10.3390/v8040114

AMA Style

Xu C, Evensen Ø, Munang’andu HM. De Novo Transcriptome Analysis Shows That SAV-3 Infection Upregulates Pattern Recognition Receptors of the Endosomal Toll-Like and RIG-I-Like Receptor Signaling Pathways in Macrophage/Dendritic Like TO-Cells. Viruses. 2016; 8(4):114. https://0-doi-org.brum.beds.ac.uk/10.3390/v8040114

Chicago/Turabian Style

Xu, Cheng, Øystein Evensen, and Hetron Mweemba Munang’andu. 2016. "De Novo Transcriptome Analysis Shows That SAV-3 Infection Upregulates Pattern Recognition Receptors of the Endosomal Toll-Like and RIG-I-Like Receptor Signaling Pathways in Macrophage/Dendritic Like TO-Cells" Viruses 8, no. 4: 114. https://0-doi-org.brum.beds.ac.uk/10.3390/v8040114

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop