- Open Access
IPNV with high and low virulence: host immune responses and viral mutations during infection
Virology Journal volume 8, Article number: 396 (2011)
Infectious pancreatic necrosis virus (IPNV) is an aquatic member of the Birnaviridae family that causes widespread disease in salmonids. IPNV is represented by multiple strains with markedly different virulence. Comparison of isolates reveals hyper variable regions (HVR), which are presumably associated with pathogenicity. However little is known about the rates and modes of sequence divergence and molecular mechanisms that determine virulence. Also how the host response may influence IPNV virulence is poorly described.
In this study we compared two field isolates of IPNV (NFH-Ar and NFH-El). The sequence changes, replication and mortality were assessed following experimental challenge of Atlantic salmon. Gene expression analyses with qPCR and microarray were applied to examine the immune responses in head kidney.
Significant differences in mortality were observed between the two isolates, and viral load in the pancreas at 13 days post infection (d p.i.) was more than 4 orders of magnitude greater for NFH-Ar in comparison with NFH-El. Sequence comparison of five viral genes from the IPNV isolates revealed different mutation rates and Ka/Ks ratios. A strong tendency towards non-synonymous mutations was found in the HRV of VP2 and in VP3. All mutations in VP5 produced precocious stop codons. Prior to the challenge, NFH-Ar and NFH-El possessed high and low virulence motifs in VP2, respectively. Nucleotide substitutions were noticed already during passage of viruses in CHSE-214 cells and their accumulation continued in the challenged fish. The sequence changes were notably directed towards low virulence. Co-ordinated activation of anti-viral genes with diverse functions (IFN-a1 and c, sensors - Rig-I, MDA-5, TLR8 and 9, signal transducers - Srk2, MyD88, effectors - Mx, galectin 9, galectin binding protein, antigen presentation - b2-microglobulin) was observed at 13 d p.i. (NFH-Ar) and 29 d p.i. (both isolates).
Mortality and expression levels of the immune genes were directly related to the rate of viral replication, which was in turn associated with sequences of viral genes. Rapid changes in the viral genome that dramatically reduced virus proliferation might indicate a higher susceptibility to protective mechanism employed by the host. Disease outbreak and mortality depend on a delicate balance between host defence, regulation of signalling cascades and virus genomic properties.
Despite vaccination programs, outbreaks of infectious pancreatic necrosis disease (IPN) are frequent in farmed salmon fry and post-smolts. Mortality rates observed in outbreaks vary considerably and have in part been ascribed to the inherited differences in susceptibility of the host [1–3]. Environmental stress [4–7] and the viral strains  also influence mortality. Atlantic salmon surviving an IPNV infection may become asymptomatic carriers of the virus for long periods [9, 10]. The production of virus may increase under stress, and carriers can shed the virus and infect surrounding fish .
The IPN virus (IPNV) is a bi-segmented double-stranded RNA (dsRNA) virus in the family Birnaviridae encoding 5 viral proteins. Segment B encodes the RNA-dependent RNA polymerase VP1. Segment A encodes a polyprotein which is cotranslationally cleaved by the viral encoded serine-lysine protease (VP4) releasing the proteins pVP2 and VP3 [12, 13]. pVP2 is further processed by host cell proteases to form the mature outer capsid protein VP2 , which is the most abundant virus protein and contains the antigenic regions responsible for induction of neutralizing antibodies in the host . VP3 is the inner structural protein, which bound to dsRNA constitutes the ribonucleoprotein core structure . Additionally VP3 is shown to bind VP1 and to self-associate strongly, indicating that it is a matrix protein . An alternative open reading frame (ORF) on Segment A encodes the small, arginine-rich, non-structural protein VP5. The biological function of IPNV VP5 remains to be determined.
The molecular basis of IPNV virulence and its interplay with host antiviral mechanisms are not fully understood. Sano and co-workers  were the first to suggest that the virulence of IPNV is associated with Segment A. Several studies using nucleotide sequence analyses have confirmed this and have shown that the VP2 residues 217 and 221 are the major determinant of virulence of IPNV serotype Sp strains. In addition, position 247 was seen as highly variable . Highly virulent isolates possess residues Thr217 and Ala221; moderate- to low-virulence strains have Pro217 and Ala221, while the strains containing Thr221 are almost avirulent, irrespective of the residue at position 217. IPNV isolates also differ in properties related to replication rate and the ability to cause persistent infections. These characteristics can be attributed to the same amino acids as those determining the virulence . Although some of the factors behind these mechanisms are known there are still many questions to be answered.
During viral infections the initial response of the immune system is the induction of type I interferons (IFN), which mediate antiviral and immunomodulatory activity. In Atlantic salmon three different subtypes of type I IFN have been identified: IFN-a, b and c . IFN-a1 and c are both expressed in head kidney and are induced by poly I:C . IFN-a1 has been shown to provide protection against IPNV in salmonid cells [21, 22]. The generation of anti-viral responses during infections requires a rapid viral sensing by pattern recognition receptors (PRRs). Toll-like receptors (TLRs) on the cell-surface or within endosomes recognize single-stranded RNA (ssRNA) and dsRNA , while the helicases RIG-I and MDA5 recognize ssRNA and dsRNA in the cytosol . Additionally, dsRNA are recognized by PKR . A number of PRRs have been identified in Atlantic salmon including RIG-I , MDA5 (GenBank: EG820831), PKR (GenBank: EF523422), TLR3 , TLR8 , TLR9  and TLR22 (GenBank: CAJ80696 and FM206383). The latter is a dsRNA-specific PRR found exclusively in lower vertebrates . Several studies have shown strong activation of immune genes upon challenge with highly virulent IPNV isolates [31, 32], and type I IFNs and the IFN-inducible Mx gene were among the most highly up-regulated genes. However, it remains to find, which PRRs are required for the induction of a systemic type I IFN response during IPNV infections.
Little is known about the relationship between the host responses and the virulence of different IPNV isolates. The latter can be associated with either down-regulation or excessive stimulation of innate immunity. Studies of IPNV infected cell-lines [33, 34] have shown inhibition of IFN signaling. In the current work we have assessed immune gene expression changes during an experimental challenge of salmon post-smolts with both a virulent and an avirulent IPNV field strain. In addition to quantitative real-time RT PCR (qPCR) analyses of selected genes we used cDNA microarray, which expanded the repertoire of genes.
The two IPNV isolates used in this study were originally collected from field outbreaks of IPN with 32% and 5% reported mortality, respectively. Initial characterization of the isolates in our lab showed that they both contained the high virulent motif Thr217 Ala221 Thr247 in VP2 defined by Santi et al . A previous bath challenge performed by our group included these isolates, and indicated that they differed substantially in virulence levels reflecting the field mortalities (unpublished). Sequencing of VP2 after this challenge confirmed the Thr217 Ala221 motif in both isolates. However, propagation of the isolates in cell-cultures before the challenge experiment reported here changed the isolate with the initial low mortality to contain a Thr217 Thr221 Thr247 motif (low virulent motif). Santi et al  suggested that the Ala-to-Thr substitution at position 221 in VP2 is a molecular determinant for the establishment of a persistent IPNV infection. Thus, we have analyzed the IPNV VP2 nucleotide composition in virus recovered from fish head kidney during the infection. Comparing sequences to the input strains could provide information about the rate and character of virus changes in vivo.
Two isolates of IPNV gave divergent mortality rates in a bath challenge
Before initiation of this challenge, head kidney from 10 fish was sampled and checked for the presence of IPNV by carrier tests in pools of 3 individuals. The lower detection limit was 0.88 virus particles g-1 tissue. No IPNV positive pools were found. In the bath challenge study, mortality caused by IPNV was initiated 12 days post infection (d p.i.) with the NFH-Ar isolate and there was a rapid increase in mortality until day 21. The experiment was terminated 30 days post IPNV challenge which was one week after the fish in both groups had ceased to die. At that point, fish infected with the NFH-Ar isolate showed far higher cumulative mortality (38.1%, average of two tanks) compared to the NFH-El isolate which gave only 1.6% cumulative mortality on average (Figure 1A). No mortality was observed among the control fish.
Virus quantification in pancreas and head kidney
The pancreatic tissue was used for quantification of virus as it is one of the primary target organs for IPNV replication. The viral titer was determined as infectious particles per gram tissue from 6 fish infected with each isolate at day 6, and from 8 fish at day 13 and 29. Viral titers increased for both isolates from day 6 to day 13, although for NFH-El (low-mortality groups) the replication was slower compared to NFH-Ar (high-mortality group). A further increase in titer from day 13 to day 29 was observed for NFH-El, whereas NFH-Ar showed decreased titers at 29 d p.i. (Figure 1B). The number of NFH-Ar infectious particles was 3.4 × 103/g in average at day 6, increased to 2.2 × 108 at day 13 and dropped to 8.5 × 105 by the third sampling at day 29. Whereas the number of infectious particles for the NFH-El isolate was slightly higher at the first sampling (average 1.6 × 104/g), only increasing to 3.3 × 104 at day 13 and further leaps to 6.4 × 107 at day 29. The individual variation was greater in the groups infected with NFH-El than the NFH-Ar infected groups.
The levels of IPNV VP2 RNA in the head kidney of NFH-Ar and NFH-El challenged fish were quantified by qPCR during the time-course. We calculated the relative quantity of VP2 by comparing the expression levels in the different samples to NFH-El sampled at 29 d.p.i. VP2 RNA was detected neither in the control fish nor in the two infected groups at 6 d p.i. The qPCR results for other time-points were overall in agreement with the viral titers determined in pancreas samples from the infected fish. The group infected with NFH-Ar had prevalence of 100% (7/7) at day 13 p.i. with an average 1000 fold higher VP2 RNA load to NFH-El 29 d p.i; while at day 29 p.i. the prevalence was 57% (4/7) and VP2 RNA levels were approximately 7 fold higher compared to NFH-El 29 d p.i. (Figure 1C, D). For NFH-El a decrease in IPNV VP2 RNA levels was observed from 13 d p.i. to 29 d p.i. in head kidney. At day 13 d p.i. the prevalence was 100% (8/8) for NFH-El (2 fold higher virus load than 29 d p.i.), while at day 29 it had decreased to 37.5% (3/8)
Sequence comparison of the viral genes
The sequences encoding the individual viral proteins were determined for each of the isolates to evaluate the rate and character of changes and to search for the differences associated with virulence. The virus isolates were sequenced at several stages of the experiment: initially after the field outbreaks, then after propagating the virus in the CHSE-214 cell-line, before the challenge and finally, in the infected fish. The rate of sequence changes was lowest in VP4 (7.4 sites/kb) being 2-2.7 times greater in other genes (Table 1). VP1 showed a higher rate of synonymous (silent) substitutions (Ka/Ks = 0.03). For VP2 and VP3 the non-synonymous (amino-acid replacement) nucleotide substitution rate was much higher (one order of magnitude greater in VP2 and VP3), and its distribution along the sequences strongly suggested presence of hyper variable regions (HVR). Importantly, the first of two HVRs of VP2 with the highest level of Ka/Ks (Figure 2) coincided with the region that is regarded as a determinant of virulence [8, 19]. All mutations in VP5 produced premature stop codons. The deduced amino acid sequences for VP2 and VP5 from the two isolates were aligned to each other and to the reference strain N1 before and after the challenge (Figure 3A, B, C and 3D). After sampling from the original field outbreaks (which were reported to give 32% mortality for the NFH-Ar and 5% mortality for the NFH-El isolate) the isolates were passaged once in CHSE-214 cells before sequence analyses. The VP2 sequences revealed a Thr217-Ala221-Thr247 motif in both isolates (origin) (Figure 3A, B and 3C). Originally, sequencing of VP5 NFH-Ar gave an ORF resulting in a 133 amino acid protein (15.2 kDa predicted), whereas VP5 from NFH-El had a premature stop-codon resulting in a very short ORF, encoding only 28 amino acid residues (3.3 kDa predicted) (Figure 3A and 3D). In the course of this study sequence changes in VP2 and VP5 accumulated continuously, during passage in culture and infection in the fish (Figure 3A). In order to obtain sufficient amount of virus to perform the bath challenge, the isolates were passaged additional times in CHSE-214 cells (2x for NFH-Ar and 3x and NFH-El). Sequencing of this material (challenge input) revealed a change into a Thr217-Thr221-Thr247 motif in VP2 for the NFH-El and a partial shift from Ala221 to Thr221 in NFH-Ar (Figure 3A, B and 3C). Additionally a partial shift from an Arg (R) in position 106 of VP5 to a stop-codon in this position was observed in NFH-Ar (Figure 3A and 3D). After the challenge viral RNA/cDNA from the head kidney of infected fish (n = 8) was sequenced and the VP2 motifs changed into Pro217-Ala221-Ala247 in both isolates during the infection process in the fish (Figure 3A, B and 3C). VP5 in NFH-Ar acquired a stop-codon in position 106, resulting in a 12.1 kDa protein, whereas NFH-El maintained it's short (3.3 kDa) VP5 (Figure 3A and 3D).
Microarray analyses revealed 211 differentially expressed genes, of which 81 and 130 genes responded to both isolates or to only NFH-Ar respectively; we did not find genes that were regulated exclusively in the NFH-El infected salmon. Both isolates affected genes within different functional groups. Equal down-regulation was observed in genes encoding motor proteins (myosins, tropomyosins, troponins), plasma proteins (plasminogen precursor, coagulation factor X, angiotensin I converting enzyme) and a suite of enzymes involved in xenobiotic metabolism and detoxification (Additional file 1). Similar changes were seen in immune related genes implicated in different processes (Table 2). It is worth noticing that genes involved in biosynthesis of eicosanoids, the inflammatory regulators of lipid origin and several complement factors are down-regulated. At the same time the virus responsive genes (VRG) were up-regulated only in fish infected with NFH-Ar. Some VRGs have well established roles in antiviral responses: antigen presentation (beta-2 microglobulin, MHCI, proteasome subunit), signal transduction downstream of IFN (Jak, STAT) and editing viral RNA (double-stranded RNA-specific adenosine deaminase - ADAR); several more genes are known as IFN-dependent. This group also includes genes with undetermined functions, which however have shown consistent induction in other virus infection studies that used the same microarray [35, 36] or the new oligonucleotide platform ; these are galectins and galectin binding protein, a tripartite motif protein and SRK2-like tyrosine kinase. Six immune genes were chosen for verification of microarray analyses and results produced with two independent methods were in a good agreement (Figure 4A). At 13 d p.i. significant up-regulation was observed only in salmon challenged with NFH-Ar. At 29 d p.i. expression changes were smaller with no difference between the isolates. The greatest expression changes were seen in SRK2.
The expression of selected immune genes not presented on the microarray was analyzed by qPCR (Figure 4B, C). The TLR8, 9, 22, RIG-I, MDA5 and PKR are implicated in viral detection [30, 38] and MyD88 is a TLR adaptor signaling molecule . Interestingly, all these nucleic acid recognition molecules were slightly down-regulated at 6 d p.i. in both infected groups (Figure 4B). At the later time-points, TLR22 did not show difference between the infected fish and controls. PKR was down-regulated in the infected fish both at day 6 and day 29 p.i., while a modest but significant (p < 0.05) increase in PKR levels could be detected at day 13 p.i for the NFH-El isolate. The remaining PRRs tested had expression profiles that were similar to those observed in other VRG (Figure 4A, C). TLR8 was up-regulated in NFH-Ar infected fish at both time-points which was significant (p < 0.05) for 13 d p.i, while expression change of TLR9 was significant (p < 0.05) at 29 d p.i. The cytoplasmic RNA sensors RIG-I and MDA5 showed 5-fold induction at 13 d p.i. in the NFH-Ar group which was significant (p < 0.05) for RIG-I. MDA5 also showed a large increase at 29 d p.i. in two individuals infected with NFH-El.
The antiviral genes Mx and IFN-γ (type II IFN) were the most highly induced transcripts in the study, and the NFH-Ar infected fish showed the greatest expression changes at 13 d p.i.; 280-fold for Mx and a significant (p < 0.05) 13-fold for IFN-γ (Figure 4C). Their expression levels decreased at 29 d p.i. to 13- and 3-fold up-regulation respectively (Figure 4C). In the NFH-El group Mx and IFN-γ showed highest expression at day 29 p.i. with 29- and 5-fold up-regulation respectively. As illustrated in Figure 5, exemplified by Mx and IFN-γ, there were variations in immune gene expression between individuals in both groups. Nevertheless, for most individuals the immune gene expression mirrored the virus levels, both the IPNV titers in pancreas and VP2 transcript levels in head kidney. Three type I IFN genes had different expression profiles during the challenge. While IFN-c did not respond to infection, IFN-a1 showed a moderate up-regulation both at day 13 (3.5-fold) and day 29 p.i. (2.5-fold) in NFH-Ar group, and at day 29 p.i. in the NFH-El group (5.9 fold). Expression of IFN-b in the head kidney was consistently low or undetected; however, an induction was observed at 13 d p.i. in three of eight NFH-Ar-infected fish (results not shown). Secreted IgM was slightly up-regulated at 29 d p.i for both isolates at similar levels (about 3-fold).
IPNV isolates have been reported to vary considerably in their virulence and pathogenicity for Atlantic salmon [18, 19]. Previous studies have indicated that a combination of IPNV virulence and host pathogen interactions determines the outcome of IPNV infections. Occurrence of multiple strains indicates rapid changes of the virus but little is known about the rate and character of mutations. It also remains undefined if susceptibility of salmon to IPNV could be associated with low or excessive immune responses. In this study two Norwegian field isolates of IPNV with marked difference in mortality were applied for experimental infection of Atlantic salmon. For the first time the sequence changes and immune responses to low and high virulence strain were compared within a single challenge test.
Previous studies have identified the outer capsid protein VP2 as the main determinant of IPNV virulence as it comprises all the neutralizing epitopes and cell attachment sites that determine host or cell specificity . The amino acid signature associated with virulence of different strains is also identified within the VP2 region [8, 19, 39]. The sequences of the two isolates in this study were determined before and after the challenge. Initially NFH-Ar had Thr217-Ala221-Thr247 while NFH-El had Thr217-Thr221-Thr247, implying high and low virulence of NFH-Ar and NFH-El respectively (Figure 3A). Thus the results from this study are in accordance with previously reports of the virulence motifs within VP2 . Rapid changes were seen during passage in cell culture and also during infection in the fish. After the challenge both isolates acquired Pro217-Ala221-Ala247 motif in VP2 associated with the moderate to low virulence . Sequence analyses revealed a remarkably high rate of non-synonymous substitutions in the HVR containing the virulence motif. Commonly, high Ka/Ks ratios point to the divergent selection meaning that sequence changes increase the rate of reproduction. However in this study we observed accumulation of mutations that correlated with a reduction in virus proliferation, thus being favorable for the fish. End of mortality in NFH-Ar group coincided with the loss of the high virulence motif. To explain this finding, virus modification as the host's defense mechanisms can be hypothesized. Virus editing is a rapidly expanding research area. At present, the best studied actors are adenosine deaminases (ADARs) that target regions of dsRNA, converting adenosine (A) to inosine (I) resulting in an A to guanosine (G) change after second strand synthesis . ADARs target mRNAs, transposable elements and RNA viruses' genomes. In mammals, several families of viruses show A to G mutations thought to be caused by ADARs . The induction of ADAR in the NFH-Ar infected fish observed in the microarray in this study might indicate a possible role for ADAR in editing IPNV. Most of the changes observed in our study, disregarding whether synonymous or non-synonymous, were from A to G or from T to C (or vice versa) associated with deamination (results not shown). However, whether the mutations in IPNV detected in this study are caused by salmon ADARs is an interesting question that needs to be addressed in future studies.
The two virus isolates used in this study showed greater discrepancies in the VP5 region than previously described isolates. The change of VP5 into a shorter protein in the NFH-Ar infected fish might imply a more benign virus since the longer 15.2 kDa VP5 protein was shown to have a potential antagonistic effect on the IFN response by inhibiting IFN-induced expression from the Mx promoter  and might thereby benefit the viral replication. However, the domains responsible for this function have not been mapped and can still be present in the shorter 12.1 kDa version. An early stop-codon located in NFH-El leaves this isolate with a severely truncated form of VP5, only 28 amino acids long, and to our knowledge such a mutation has not been reported in other surveys of IPNV field isolates. It has been suggested that VP5 has an anti-apoptotic function, which is probably not essential for the virulence or persistence of the virus [42, 43]. Functional significance of the observed changes in VP5 remained unclear. The amino acid sequences of VP1 ORF were identical between the two isolates throughout the study and VP4 were subject to very few mutations.
The mortality rates associated with the structural differences described above may be linked to the ability of virus to invade and replicate within the host cells and/or the scale and character of immune responses. Association between mortality and the virus titers was obvious (Figure 1). NFH-El characterized with low replication was avirulent while NFH-Ar infection was fatal at high rate of proliferation at 13 d p.i. To assess the immune responses, we used expression profiling with microarray and qPCR analyses of genes with well-established roles (IFNs, Mx and PRRs) and both approaches produced similar results. The expression levels of VRGs were apparently mirroring the viral titers and expression levels of the IPNV VP2. At 13 d p.i. VRG were induced in salmon infected with NFH-Ar but not with NFH-El consistently with the difference of virus titers. Up-regulation of VRG in NFH-E1 infected fish at 29 d p.i. was in line with the slight increase of virus titer and did not affect mortality. It is likely that the slower replication and concomitant slower spread of the NFH-El strain may allow time for a systemic induction of the host anti-viral system, including adaptive responses.
The IFN system is believed to have a crucial role in the first line of defense against virus infections, and in vitro studies have demonstrated that IPNV replication in cell cultures is efficiently inhibited by salmon IFN-a1 [21, 22]. Additionally, injection of synthetic IFN-inducers like CpG and poly I:C induce protection against IPNV in Atlantic salmon . In this study IFN-a1 was induced by both viral strains and major up-regulation was seen in the IFN-a1 dependent gene Mx. IFN-c, was not induced and even slightly down-regulated at 29 d p.i. for NFH-Ar. IFN-c is suggested to have a separate regulation from a and b and can be produced by a different cell population than IFN-a1 . IFN-b was not detected in this study or showed consistently low expression levels (data not shown). Despite the suggested role of type I IFN in restraining virus production, results in our lab has demonstrated that IFN-a1 does not completely inhibit IPNV growth but causes a delay in viral protein synthesis . Furthermore, our data suggest that IPNV-encoded proteins may be involved in weakening of IFN signaling . As a result high levels of viral proteins may impair the activity of IFN-induced genes, thus the higher replication rate of NFH-Ar compared to NFH-El may cause a more potent IFN-antagonizing effect of the NFH-Ar strain. However when interpreting the results from live pathogen challenges, it is important to keep in mind the complexity of such studies, where it is not straight forward whether an observed response is a strategy employed by the virus for its own benefit or a response by the host to control the virus.
Although innate immunity by itself represents a powerful system to combat viral invaders, many infections can only be cleared in combination with adaptive immunity. In this regard type I IFNs are known to promote the adaptive arm including both T cell mediated cellular responses and antibody production [45, 46]. It is likely that the reduction of virus titers by 29 d p.i. could be associated with the onset of adaptive immune responses. Unlike the genes implicated in the innate immunity, sIgM showed no expression changes at 13 d p.i. and was slightly induced at 29 d p.i.
Results of this study added knowledge to the understanding of the immune responses after IPNV infections. Expression of a panel of PRRs was assessed including several recently identified genes. TLR8, 9 RIG-I and MDA5 showed up-regulation and followed the same trend as other immune genes. Earlier we observed a modest increase of TLR8 and 9 expression during stimulation and infection [27, 28]. TLR22 is reported to recognize dsRNA in pufferfish and when over-expressed it induces type I IFN expression upon IPNV infection, which suggests a possible role for TLR22 in protection against IPNV . However, TLR22 was not induced by IPNV in this study. Unexpectedly, PKR, an IFN-inducible gene, was down-regulated at all time points except 13 d p.i. (Figure 4A). This was in contrast to Mx, another IFN-inducible gene and IFN-a1, which were up-regulated. Functional studies of salmon PKR have to our knowledge not been reported, however flounder PKR was up-regulated both in vitro and in vivo by a negative single stranded RNA virus (SMRV), which also induced Mx expression in vitro. PRRs, PKR and MyD88 were down-regulated at 6 d p.i. in both study groups. This could be explained by migration of leukocytes expressing the PRRs from the head-kidney into the bloodstream at early time-points of the infection.
IFN-γ was together with Mx the most highly induced immune gene in this study, and high levels of IFN-γ has been reported in other IPNV challenge experiments . IFN-γ is regarded as a typical Th1 cytokine which bridges the innate and adaptive immune responses. Fish IFN-γ share several functional properties with mammalian IFN-γ including macrophage activation [48–50] and rainbow trout IFN-γ is shown to signal through STAT1 [50, 51]. Recently, we have observed antiviral activity against IPNV by IFN-γ, although the effect was not as pronounced as described for IFN-a1 . Like in mammals, fish IFN-γ plays diverse roles in different facets of the immune system, and the increased levels of IFN-γ upon IPNV challenge detected here suggest anti-IPNV activity.
Our study aimed at investigation of different mortality associated with two IPNV isolates. High incidence of death of NFH-Ar infected salmon was in line with the initial presence of the virulence motif in VP2 and higher virus titers in the head kidney and pancreas compared to NHF-El. PPRs, IFNs and ISGs were induced to a higher magnitude by the high virulent strain NFH-Ar and correlated with the viral load in the head kidney. Subsequent decrease of virus titers and mortality could be due to the sequence changes in the hyper variable regions of IPNV together with the development of acquired immunity in fish.
Two IPNV isolates, referred to as NFH-Ar and NFH-El, collected from field outbreaks of IPN in 2004 were used in this study. Veterinarians had diagnosed IPN based on clinical observations and by use of an agglutination assay (Phadebact coating kit, Boule Diagnostics AB). The NFH-Ar isolate, obtained from an IPN outbreak in Frøya, Norway in 2004, was reported to give 32% mortality, while the NFH-El, obtained from an IPN outbreak in Alta, Norway in 2004, was reported to give 5% mortality. Head kidney samples collected during these outbreaks were sent to our lab and tissue homogenates were made and inoculated on CHSE-214 cells for propagation and sequencing. For NFH-Ar a second, and for NFH-El a third cell-culture passage of the virus strain was used in the present experiment. The input strains were sequenced before challenge.
IPNV challenge of Atlantic salmon
The challenge was carried out at Tromsø Aquaculture Research Station (Tromsø, Norway). Non-vaccinated Atlantic salmon, Salmo salar L., strain Aquagen IPNV sensitive (Aquagen, Kyrksæterøra, Norway), with an average size of 51 g, was used for the challenge two days after transfer to seawater. Bath challenge was performed as described by Johansen and Sommer  using an infectious dose of TCID50/ml = 5 × 105 of the field isolates NFH-Ar and NFH-El. Two parallel tanks supplied with 200 L of 10°C seawater were used for each isolate. Each tank contained 70 fish and a separate tank was used for 70 uninfected control fish. The fish were fed daily on commercial feed. The experiment was terminated 30 days after challenge.
Mortality was monitored throughout the experiment by counting dead fish. An IPNV rapid agglutination kit (Phadebact coating kit, Boule Diagnostics AB) was used on head kidney samples from dead fish to verify IPN as the cause of death. Percent cumulative mortality in infected fish was calculated compared to control fish using GraphPad Prism 4.0 tool for statistical analyses. Sampling was performed on surviving fish. Pancreatic tissue and head kidney were aseptically removed from 4 fish in each tank at 6, 13 and 29 d p.i. Pancreatic tissue was stored in L-15 medium on ice until frozen at -80°C, whereas head kidney samples were stored in RNAlater® (Ambion) at -80°C.
Quantification of IPNV by titration
Pancreatic tissue was homogenized using an Ultra Thurrax T25 basic crusher (IKA-WERKE). Homogenized tissue was diluted to 5% in Eagle's minimum essential medium (EMEM) (Gibco) supplemented with 100 μg/ml streptomycin, 60 μg/ml penicillin, 2 mM L-glutamine and 1% non-essential amino acids, before centrifuged at 15000 × g for 5 min at 4°C. Homogenates were inoculated onto CHSE-214 cells in 96 wells plates at a starting concentration of 0.5% (w/v). Individual viral titers were determined by end-point titration using 10-fold dilutions and 2 replicates and calculated by the TCID50 method .
RNA isolation and cDNA synthesis
For analyses of cellular and viral gene expression total RNA from head kidney sampled from 8 challenged and 4 control fish at each time point, was isolated using a combination of Trizol and PureLink RNA Kit (Invitrogen). RNA quality was assessed on an agarose gel and the quantity determined by NanoDrop 2000 spectrophotometry (Thermo Scientific). No RNA degradation was observed. After isolation, the RNA was DNase treated applying TURBO DNase (Ambion). cDNA was synthesized in a 25 μl reaction from 200 ng DNase treated total RNA primed with random hexamers (TaqMan Reverse Transcription Reagents kit, Applied Biosystems). The manufacturer's protocol was followed. For isolation of viral RNA after passage in cell culture QIAamp Viral RNA Mini Kit (Qiagen) was used according to the manufacturer's instructions.
Sequencing of virus isolates
Each virus isolate was sequenced from virus derived material under the following conditions; origin (virus homogenates from head kidney pooled from 2 fish in field outbreak passaged 1× in cell culture), input (virus homogenates from head kidney pooled from 2 fish in field outbreak passaged 2-3× in cell culture), during infection(head kidney from day 13 (NFH-Ar) or day 29 (NFH-El) infected fish, 8 pooled individuals from each group). Primers specific for each of the genes encoding the five IPNV proteins were used to amplify the individual genes (from head kidney cDNA pooled from 8 fish or cDNA obtained from head kidney virus homogenates from 2 fish passaged in cell culture) in a PCR reaction using Pfu Turbo Hotstart DNA polymerase (Stratagene). The PCR-products were sequenced in both directions using the same primers and BigDye3.1 chemistry and a 3100 gene analyzer. The sequences were aligned with BioEdit 7.0.5  and DnaSP V5  was used for evaluation of mutation rates and search for HVR.
Microarray analyses were performed at 13 d p.i. on head kidney samples of salmon infected with isolates NFH-Ar and NFH-El (5 individuals from each group, one microarray per individual). The salmonid fish microarray (SFA 2.0 or immunochip, Geo Omnibus GPL6154) contains 1800 unique clones printed each in 6 spot replicates. Pooled samples of uninfected salmon (equal amounts of RNA, n = 4) were used as a common reference. Test and reference RNA (10 μg) were labeled with the fluorescent dyes Cy5-dUTP and Cy3-dUTP respectively (Amersham Pharmacia), which were incorporated in cDNA using the SuperScript™ Indirect cDNA Labeling System (Invitrogen). Synthesis of cDNA was performed at 46°C for 3 hours in a 23 μl reaction volume, following RNA degradation with 2.5 M NaOH at 37°C for 15 min and alkaline neutralization with 2 M Hepes. Labeled cDNA was combined and purified with Microcon YM30 (Millipore). Microarray slides were pre-treated with 1% BSA fraction V, 5× SSC and 0.1% SDS for 30 min at 50°C and then washed with 2× SSC (3 min) followed by 0.2× SSC (3 min) at room temperature and hybridized over-night at 60°C in a cocktail containing 1.3× Denhardt's, 3× SSC, 0.3% SDS, 0.67 μg/μl polyadenylate and 1.4 μg/μl yeast tRNA. After hybridization slides were washed at room temperature in 0.5 × SSC and 0.1% SDS (15 min), 0.5 × SSC and 0.01% SDS (15 min), 0.06 × SSC (2 min) and 0.06 × SSC (1 min). Scanning was performed with GenePix 4100A and images were processed with GenePix 6.0 (Molecular Devices). The spots were filtered by criterionI (I-B)B/(S+S ≥ 0.6, where I and B are the mean signal and background intensities and SI S B are the standard deviations. Low quality spots were excluded from analyses and genes with less than 3 high quality spots on a slide were discarded. After subtraction of median background from median signal intensities, the expression ratios (ER) were calculated. Lowess normalization was performed first for the whole slide and next for twelve rows and four columns per slide. Statistical analyses were performed in two stages. First, technical accuracy was assessed by difference of log2-ER from zero in six spot replicates and genes with significant changes (Student's t-test, p < 0.05) in at least 3 of 5 individuals per group were selected. Next, analysis of biological replicates was carried out and differential expression was assessed by the mean fold change (> 1.5) and difference from control (one sample t-test, p < 0.05).
Quantitative real-time RT PCR (qPCR) analyses
PCR primers used with Sybr Green assay were designed using Vector NTI (Invitrogen) and synthesized by Invitrogen (table 3). The amplicon lengths set to be between 50 and 200 bases were checked on 1.5% agarose gel. TaqMan assays employing a hydrolysis probe were designed using Assays-by-design (Applied Biosystems) and synthesized by Applied Biosystems. PCR efficiency was calculated from tenfold serial dilutions of cDNA for each assay in triplicates. qPCR assays were conducted using 2× SYBR® Green Master Mix (Roche Diagnostics) in an optimized 12 μl reaction volume, using 1:10 diluted cDNA, with primer concentrations of 0.4-0.6 μM. PCR was performed in duplicate in 96-well optical plates on Light Cycler 480 (Roche Diagnostics) under the following conditions: 95°C for 5 min (pre-incubation), 95°C for 5 sec, 60°C for 15 sec, 72°C for 15 sec (amplification), 95°C for 5 sec, and 65°C for 1 min (melting curve). Forty cycles were performed. Assays employing a hydrolysis probe was conducted in a 20 μl reaction using 2.5 μl 1:10 fold diluted cDNA as template and 2× TaqMan Fast Universal PCR Master mix (Applied Biosystems). The expression of mRNA was measured in an ABI Prism 7500 FAST Cycler (Applied Biosystems) and the amplification profile was: 95°C for 20 sec followed by 40 cycles of 95°C for 3 sec and 60°C for 30 sec. PCR analyses of the SybrGreen assays for RIG-I, MDA5 and PKR at 13 and 29 d p.i. was also performed at the ABI Prism 7500 FAST Cycler. Relative expression of mRNA was evaluated by ΔΔCT, adjusted for PCR efficiency and normalized against the elongation factor (EF1AB). IPNV VP2 was undetected in control samples and at 6 d p.i., in order to calculate the ΔΔCT values the average of the NFH-El infected fish sampled at 29 d p.i was used as calibrator. Results were analyzed with ANOVA followed with Newman - Keuls test (p < 0.05). (See Table 1 for primers and probes used in the qPCR assays).
Guy DR, Bishop SC, Brotherstone S, Hamilton A, Roberts RJ, McAndrew BJ, Woolliams JA: Analysis of the incidence of infectious pancreatic necrosis mortality in pedigreed Atlantic salmon, Salmo salar L., populations. Journal of Fish Diseases 2006,29(11):637-647. 10.1111/j.1365-2761.2006.00758.x
Houston RD, Gheyas A, Hamilton A, Guy DR, Tinch AE, Taggart JB, McAndrew BJ, Haley CS, Bishop SC: Detection and Confirmation of a Major QTL Affecting Resistance to Infectious Pancreatic Necrosis (IPN) in Atlantic Salmon ( Salmo Salar ). Dev Biols 2008, 132: 199-204.
Ozaki A, Sakamoto T, Khoo S, Nakamura K, Coimbra MR, Akutsu T, Okamoto N: Quantitative trait loci (QTLs) associated with resistance/susceptibility to infectious pancreatic necrosis virus (IPNV) in rainbow trout ( Oncorhynchus mykiss ). Molecular Genetics and Genomics 2001,265(1):23-31. 10.1007/s004380000392
Jarp J, Gjerve AG, Olsen AB, Bruheim T: Risk factors for furunculosis, infectious pancreatic necrosis and mortality in post-smolt of Atlantic salmon, Salmo salar L. Journal of Fish Diseases 1995,18(1):67-78. 10.1111/j.1365-2761.1995.tb01267.x
Sundh H, Olsen R-E, Fridell F, Gadan K, Evensen Ø, Glette J, Taranger G-L, Myklebust R, Sundell K: The effect of hyperoxygenation and reduced flow in fresh water and subsequent infectious pancreatic necrosis virus challenge in sea water, on the intestinal barrier integrity in Atlantic salmon, Salmo salar L. Journal of Fish Diseases 2009,32(8):687-698. 10.1111/j.1365-2761.2009.01047.x
Fridell F, Gadan K, Sundh H, Taranger GL, Glette J, Olsen RE, Sundell K, Evensen Ø: Effect of hyperoxygenation and low water flow on the primary stress response and susceptibility of Atlantic salmon, Salmo salar L. to experimental challenge with IPN virus. Aquaculture 2007, 270:(1-4):23-35. 10.1016/j.aquaculture.2007.05.006
Taksdal T, Ramstad A, Stangeland K, Dannevig BH: Induction of infectious pancreatic necrosis (IPN) in covertly infected Atlantic salmon, Salmo salar L., post-smolts by stress exposure, by injection of IPN virus (IPNV) and by cohabitation. Journal of Fish Diseases 1998,21(3):193-204. 10.1111/j.1365-2761.1998.00092.x
Song H, Santi N, Evensen O, Vakharia VN: Molecular determinants of infectious pancreatic necrosis virus virulence and cell culture adaptation. J Virol 2005,79(16):10289-10299. 10.1128/JVI.79.16.10289-10299.2005
Munro ES, Gahlawat SK, Acosta F, Ellis AE: In infectious pancreatic necrosis virus carrier Atlantic salmon, Salmo salar L., post-smolts, almost all kidney macrophages ex vivo contain a low level of non-replicating virus. Journal of Fish Diseases 2006,29(1):43-48. 10.1111/j.1365-2761.2005.00680.x
Knott RM, Munro AL: The persistence of infectious pancreatic necrosis virus in Atlantic salmon. Vet Immunol Immunopathol 1986,12(1-4):359-364. 10.1016/0165-2427(86)90141-8
Saint-Jean SR, Borrego JJ, Perez-Prieto SI, Karl Maramorosch FAM, Aaron JS: Infectious Pancreatic Necrosis Virus: Biology, Pathogenesis, and Diagnostic Methods. In Advances in Virus Research. Volume 62. Academic Press; 2003:113-165.
Duncan R, Dobos P: The nucleotide sequence of infectious pancreatic necrosis virus (IPNV) dsRNA segment A reveals one large ORF encoding a precursor polyprotein. Nucleic Acids Res 1986,14(14):5934. 10.1093/nar/14.14.5934
Duncan R, Nagy E, Krell PJ, Dobos P: Synthesis of the infectious pancreatic necrosis virus polyprotein, detection of a virus-encoded protease, and fine structure mapping of genome segment A coding regions. J Virol 1987,61(12):3655-3664.
Magyar G, Dobos P: Expression of infectious pancreatic necrosis virus polyprotein and VP1 in insect cells and the detection of the polyprotein in purified virus. Virology 1994,198(2):437-445. 10.1006/viro.1994.1055
Heppell J, Tarrab E, Lecomte J, Berthiaume L, Arella M: Strain Variability and Localization of Important Epitopes on the Major Structural Protein (VP2) of Infectious Pancreatic Necrosis Virus1. Virology 1995,214(1):40-49. 10.1006/viro.1995.9956
Hjalmarsson A, Everitt E: Identification of IPNV-specified components released from productively infected RTG-2 cells following massive cytopathic effect. Arch Virol 1999,144(8):1487-1501. 10.1007/s007050050605
Pedersen T, Skjesol A, Jørgensen JB: VP3, a structural protein of infectious pancreatic necrosis virus, interacts with RNA-dependent RNA polymerase VP1 and with double-stranded RNA. Journal of Virology 2007,81(12):6652-6663. 10.1128/JVI.02831-06
Sano T, Tanaka K, Fukuzaki S: Immune response in adult trout against formalin-killed concentrated IPNV. Developments in Biological Standardization 1981, 49: 63-70.
Santi N, Vakharia VN, Evensen O: Identification of putative motifs involved in the virulence of infectious pancreatic necrosis virus. Virology 2004,322(1):31-40. 10.1016/j.virol.2003.12.016
Sun B, Robertsen B, Wang Z, Liu B: Identification of an Atlantic salmon IFN multigene cluster encoding three IFN subtypes with very different expression properties. Developmental & Comparative Immunology 2009,33(4):547-558. 10.1016/j.dci.2008.10.001
Ooi EL, Hirono I, Aoki T: Functional characterisation of the Japanese flounder, Paralichthys olivaceus , Mx promoter. Fish & Shellfish Immunology 2006,21(3):293-304. 10.1016/j.fsi.2005.12.004
Robertsen B, Bergan V, Rokenes T, Larsen R, Albuquerque A: Atlantic salmon interferon genes: cloning, sequence analysis, expression, and biological activity. J Interferon Cytokine Res 2003,23(10):601-612. 10.1089/107999003322485107
Alexopoulou L, Holt AC, Medzhitov R, Flavell RA: Recognition of double-stranded RNA and activation of NF-kappaB by Toll-like receptor 3. Nature 2001,413(6857):732-738. 10.1038/35099560
Kato H, Takeuchi O, Sato S, Yoneyama M, Yamamoto M, Matsui K, Uematsu S, Jung A, Kawai T, Ishii KJ, et al.: Differential roles of MDA5 and RIG-I helicases in the recognition of RNA viruses. Nature 2006,441(7089):101-105. 10.1038/nature04734
Meurs E, Chong K, Galabru J, Thomas NSB, Kerr IM, Williams BRG, Hovanessian AG: Molecular cloning and characterization of the human double-stranded RNA-activated protein kinase induced by interferon. Cell 1990,62(2):379-390. 10.1016/0092-8674(90)90374-N
Biacchesi S, LeBerre M, Lamoureux A, Louise Y, Lauret E, Boudinot P, Bremont M: Mitochondrial Antiviral Signaling Protein Plays a Major Role in Induction of the Fish Innate Immune Response against RNA and DNA Viruses. J Virol 2009,83(16):7815-7827. 10.1128/JVI.00404-09
Skjæveland I, Iliev DB, Strandskog G, Jørgensen JB: Identification and characterization of TLR8 and MyD88 homologs in Atlantic salmon ( Salmo salar ). Developmental & Comparative Immunology 2009,33(9):1011-1017. 10.1016/j.dci.2009.04.007
Skjæveland I, Iliev DB, Zou J, Jørgensen T, Jørgensen JB: A TLR9 homolog that is up-regulated by IFN-[gamma] in Atlantic salmon ( Salmo salar ). Developmental & Comparative Immunology 2008,32(6):603-607. 10.1016/j.dci.2007.10.011
Rebl A, Siegl E, Köllner B, Fischer U, Seyfert H-M: Characterization of twin toll-like receptors from rainbow trout ( Oncorhynchus mykiss ): Evolutionary relationship and induced expression by Aeromonas salmonicida salmonicida . Developmental & Comparative Immunology 2007,31(5):499-510. 10.1016/j.dci.2006.08.007
Matsuo A, Oshiumi H, Tsujita T, Mitani H, Kasai H, Yoshimizu M, Matsumoto M, Seya T: Teleost TLR22 Recognizes RNA Duplex to Induce IFN and Protect Cells from Birnaviruses. J Immunol 2008,181(5):3474-3485.
Ingerslev H-C, Rønneseth A, Pettersen EF, Wergeland HI: Differential Expression of Immune Genes in Atlantic Salmon ( Salmo salar L.) Challenged Intraperitoneally or by Cohabitation with IPNV. Scandinavian Journal of Immunology 2009,69(2):90-98. 10.1111/j.1365-3083.2008.02201.x
McBeath AJA, Snow M, Secombes CJ, Ellis AE, Collet B: Expression kinetics of interferon and interferon-induced genes in Atlantic salmon ( Salmo salar ) following infection with infectious pancreatic necrosis virus and infectious salmon anaemia virus. Fish & Shellfish Immunology 2007,22(3):230-241. 10.1016/j.fsi.2006.05.004
Collet B, Munro ES, Gahlawat S, Acosta F, Garcia J, Roemelt C, Zou J, Secombes CJ, Ellis AE: Infectious pancreatic necrosis virus suppresses type I interferon signalling in rainbow trout gonad cell line but not in Atlantic salmon macrophages. Fish & Shellfish Immunology 2007, 22:(1-2):44-56. 10.1016/j.fsi.2006.03.007
Skjesol A, Aamo T, Hegseth MN, Robertsen B, Jørgensen JB: The interplay between infectious pancreatic necrosis virus (IPNV) and the IFN system: IFN signaling is inhibited by IPNV infection. Virus Research 2009,143(1):53-60. 10.1016/j.virusres.2009.03.004
Schiøtz BL, Jørgensen SM, Rexroad C, Gjøen T, Krasnov A: Transcriptomic analysis of responses to infectious salmon anemia virus infection in macrophage-like cells. Virus Research 2008, 136:(1-2):65-74. 10.1016/j.virusres.2008.05.007
Jorgensen S, Afanasyev S, Krasnov A: Gene expression analyses in Atlantic salmon challenged with infectious salmon anemia virus reveal differences between individuals with early, intermediate and late mortality. BMC Genomics 2008,9(1):179. 10.1186/1471-2164-9-179
Krasnov A, Timmerhaus G, Afanasyev S, Jørgensen SM: Development and assessment of oligonucleotide microarrays for Atlantic salmon ( Salmo salar L.). Comparative Biochemistry and Physiology Part D: Genomics and Proteomics 2011,6(1):31-38. 10.1016/j.cbd.2010.04.006
Rebl A, Goldammer T, Seyfert H-M: Toll-like receptor signaling in bony fish. Veterinary Immunology and Immunopathology 2010, 134:(3-4):139-150. 10.1016/j.vetimm.2009.10.003
Blake S, Ma JY, Caporale DA, Jairath S, Nicholson BL: Phylogenetic relationships of aquatic birnaviruses based on deduced amino acid sequences of genome segment A cDNA. Dis Aquat Organ 2001,45(2):89-102.
Keegan LP, Gallo A, O'Connell MA: The many roles of an RNA editor. Nat Rev Genet 2001,2(11):869-878. 10.1038/35098584
Samuel CE: Adenosine deaminases acting on RNA (ADARs) are both antiviral and proviral. Virology 2011,411(2):180-193. 10.1016/j.virol.2010.12.004
Hong JR, Wu JL: Induction of apoptotic death in cells via Bad gene expression by infectious pancreatic necrosis virus infection. Cell Death Differ 2002,9(2):113-124. 10.1038/sj.cdd.4400933
Santi N, Song H, Vakharia VN, Evensen O: Infectious pancreatic necrosis virus VP5 is dispensable for virulence and persistence. J Virol 2005,79(14):9206-9216. 10.1128/JVI.79.14.9206-9216.2005
Jørgensen JB, Johansen LH, Steiro K, Johansen A: CpG DNA induces protective antiviral immune responses in Atlantic salmon ( Salmo salar L.). J Virol 2003,77(21):11471-11479. 10.1128/JVI.77.21.11471-11479.2003
Biron CA: Interferons alpha and beta as immune regulators--a new look. Immunity 2001,14(6):661-664. 10.1016/S1074-7613(01)00154-6
Prchal M, Pilz A, Simma O, Lingnau K, von Gabain A, Strobl B, Müller M, Decker T: Type I interferons as mediators of immune adjuvants for T- and B cell-dependent acquired immunity. Vaccine 2009, 27:(Supplement 6):G17-G20.
Zhu R, Zhang Y-B, Zhang Q-Y, Gui J-F: Functional Domains and the Antiviral Effect of the Double-Stranded RNA-Dependent Protein Kinase PKR from Paralichthys olivaceus . J Virol 2008,82(14):6889-6901. 10.1128/JVI.02385-07
Arts JAJ, Tijhaar EJ, Chadzinska M, Savelkoul HFJ, Verburg-van Kemenade BML: Functional analysis of carp interferon-gamma: Evolutionary conservation of classical phagocyte activation. Fish & Shellfish Immunology 2010,29(5):793-802. 10.1016/j.fsi.2010.07.010
Grayfer L, Garcia EG, Belosevic M: Comparison of macrophage antimicrobial responses induced by type II interferons of the goldfish ( Carassius auratus L.). J Biol Chem 2010,285(31):23537-23547. 10.1074/jbc.M109.096925
Zou J, Carrington A, Collet B, Dijkstra JM, Yoshiura Y, Bols N, Secombes C: Identification and bioactivities of IFN-gamma in rainbow trout Oncorhynchus mykiss : the first Th1-type cytokine characterized functionally in fish. J Immunol 2005,175(4):2484-2494.
Skjesol A, Hansen T, Shi CY, Thim HL, Jorgensen JB: Structural and functional studies of STAT1 from Atlantic salmon ( Salmo salar ). BMC Immunol 2010, 11: 17. 10.1186/1471-2172-11-17
Sun B, Skjaeveland I, Svingerud T, Zou J, Jorgensen J, Robertsen B: Antiviral Activity of Salmonid Interferon Gamma against Infectious Pancreatic Necrosis Virus and Salmonid Alphavirus and its Dependency on Type I Interferon. Journal of Virology 2011.
Johansen L-H, Sommer A-I: Multiplication of infectious pancreatic necrosis virus (IPNV) in head kidney and blood leucocytes isolated from Atlantic salmon, Salmo salar L. Journal of Fish Diseases 1995, 18: 147-156. 10.1111/j.1365-2761.1995.tb00272.x
Reed L, Muench H: A simple method for estimating fifty percent endpoints. American Journal of Hygiene 1938, 27: 493-497.
Hall TA: BioEdit: a user-friendly biological sequence alignment editor and analysis program for Windows 95/98/NT. Nucleic Acids Symposium Series 1999, 41: 95-98.
Rozas J: DNA Sequence Polymorphism Analysis Using DnaSP. Bioinformatics for DNA Sequence Analysis 2009, 537: 337-350. 10.1007/978-1-59745-251-9_17
This work was financed by The Research Council of Norway (Grant no. 183230/S10 and 185217)/S40). We appreciate the laboratory technical assistance from T. Aamo and H. L. Thim and the assistance in the IPNV challenge from R. I. Hansen and M. Marienborg, Tromsø Aquaculture Research Station.
The authors declare that they have no competing interests.
AS performed most of the sequencing, assisted with the challenge experiment and quantification of virus, in addition he drafted the manuscript together with IS. qPCR analyses were performed by IS and GT. ME and BNF performed the challenge experiment and some of the sequencing and virus quantification. The microarray analyzes were done by SMJ and GT. AK contributed to this work with gene expression and sequence data and drafting the manuscript. JBJ designed the experiments, analyzed data and contributed to finalizing the manuscript. All authors read and approved the final manuscript.
Astrid Skjesol, Ingrid Skjæveland contributed equally to this work.