- Benchmark Genetics Norway AS, Bergen, Norway
Infectious pancreatic necrosis (IPN) is an important viral disease of salmonids that can affect fish during various life cycles. In Atlantic salmon, selecting for genetically resistant fish against IPN has been one of the most highly praised success stories in the history of fish breeding. During the late 2000s, the findings that resistance against this disease has a significant genetic component, which is mainly controlled by variations in a single gene, have helped to reduce the IPN outbreaks to a great extent. In this paper, we present the identification of a new variant of the IPN virus from a field outbreak in Western Norway that had caused mortality, even in genetically resistant salmon. We recovered and assembled the full-length genome of this virus, following the deep-sequencing of the head-kidney transcriptome. The comparative sequence analysis revealed that for the critical amino acid motifs, previously found to be associated with the degree of virulence, the newly identified variant is similar to the virus’s avirulent form. However, we detected a set of deduced amino acid residues, particularly in the hypervariable domain of the VP2, that collectively are unique to this variant compared to all other reference sequences assessed in this study. We suggest that these mutations have likely equipped the virus with the capacity to escape the host defence mechanism more efficiently, even in the genetically deemed IPN resistant fish.
Introduction
Infectious pancreatic necrosis (IPN) is one of the leading viral diseases of the farmed Atlantic salmon (Salmo salar L.). The disease was first reported in 1941, following an outbreak in brook trout (Salvelinus fontinalis L.) in Canada (M’Gonigle, 1941). Further, the characterisation of the causative agent, the IPN virus (IPNV), was performed in 1960 (Wolf et al., 1960). In Norway, the virus was first isolated in 1975 from freshwater rainbow trout (Oncorhynchus mykiss L.) (Hastein and Krogsrud, 1976) and was designated as a notifiable disease from 1991 till 2008 (Sommerset et al., 2020). The virus can infect Atlantic salmon during all of its developmental stages, but the fish are especially susceptible as fry, during start-feeding, and post-smolts, soon after transfer to seawater (Roberts and Pearson, 2005; Sommerset et al., 2020). As suggested by the name, the pancreas is the primary target tissue of IPNV. However, the liver has also been shown to be one of the key organs affected by this virus (Ellis et al., 2010; Munang’Andu et al., 2013). In addition to significant economic losses, this disease is a concern for animal welfare, as survivors following infection can continue to infect naïve fish (Roberts and Pearson, 2005).
IPNV is a double-stranded RNA virus that belongs to the Birnaviridae family. A characteristic feature of the birnaviruses is their possession of a bi-segmented genome (i.e., A and B segments), contained within a non-enveloped, icosahedral capsid (Dobos and Roberts, 1983; Duncan and Dobos, 1986; Duncan et al., 1987). The smaller genomic segment, the B segment, is about 2.5 kb in size and consists of a single open reading frame (ORF), encoding an RNA-dependent, RNA polymerase VP1 (approximately 90 kDa) (Duncan and Dobos, 1986; Duncan et al., 1987). The A-segment is approximately 3 kb, and it contains two partially overlapping ORFs (Duncan and Dobos, 1986). The first ORF encodes VP5, a small, cysteine-rich, non-structural protein of about 17 kDa. A precursor polyprotein (NH2-preVP2-NS VP4 protease-VP3-COOH) of approximately 106 kDa constitutes the second and the larger ORF. The encoded VP4 protease cleavages the polyprotein into its two main components, the preVP2 and VP3. The preVP2 further matures to become VP2 and forms the outer capsid protein, containing the neutralizing epitopes and sites that facilitate cell attachment (Heppell et al., 1995b; Dopazo, 2020). The VP3 is the inner capsid protein (Duncan and Dobos, 1986; Duncan et al., 1987) but also, in association with the VP1, it seems to be involved in viral packaging and replication (Tacken et al., 2002).
Based on cross-neutralizing tests, the aquatic birnaviruses are broadly divided into two main serogroups, A and B (Hill and Way, 1995). The serogroup A comprises of nine different serotypes, A1-A9. Most of the isolates from the United States fall within the A1 serotype (West Buxton; WB), the A2-A5 include mainly the European and Chilean isolates (A2 (Spjarup; Sp), A3 (Abild; Ab), A4 (Hecht, He) and A5 (Te)), and A6-A9 are predominantly variants from Canada (A6 (Canada 1; C1), A7 (Canada 2; C2), A8 (Canada 3; C3) and A9 (Jasper)) (Hill and Way, 1995; Blake et al., 2001; Mutoloki and Evensen, 2011). The serogroup B consists of only one serotype (B1 (Tellinavirus; TV1)) that is non-pathogenic in fish. Based on the analysis of variations in the nucleotide and the deduced amino acids of the VP2, Blake et al. (2001) later proposed that aquatic birnaviruses constitute six major genogroups. These groups correspond relatively well with the geographical origins and serological classifications previously suggested. According to this classification, WB and the Jasper strains (A1 and A9) constitute genogroup 1, the A3 serotype forms genogroup 2, Te and C1 (A5 and A6) make genogroup 3, the two Canadian strains, C2 and C3 (A7 and A8) form genogroup 4, all European and Chilean isolates of serotype A2 cluster into genogroup 5 and the He strain (A4) constitutes genogroup 6 (Blake et al., 2001). A seventh genogroup was also suggested based on isolates recovered from Japanese aquatic sources (Nishizawa et al., 2005).
Despite extensive vaccination efforts, disinfection and targeted breeding programs, IPN outbreaks remains a concern in farmed Atlantic salmon. Related mortalities can be substantial and, in part, can be a function of the genetic makeup of the host (Guy et al., 2006; Houston et al., 2008; Moen et al., 2009), management and environmental factors (Jarp et al., 1995; Sundh et al., 2009), and the strain and specific genetic variations carried by the virus (Song et al., 2005; Skjesol et al., 2011). There are different examples of isolates belonging to the same serotype to cause different degrees of mortality in the host or losing their virulence after serial passages in cell culture (Dorson et al., 1978; Bruslind and Reno, 2000; Santi et al., 2004; Shivappa et al., 2004; Song et al., 2005). Indeed, many studies have investigated the molecular basis of IPNV virulence in Atlantic salmon (e.g., Santi et al., 2004, 2005; Song et al., 2005; Skjesol et al., 2011; Mutoloki et al., 2016). Although the genetic details of such processes are not yet fully understood, most studies speculated, and chiefly investigated, the genetic variations in segment A concerning the degree of virulence in IPNV. In particular, specific amino acid motifs in the VP2 of the polyprotein have been suggested to be most prevalent and in association with the virulent forms of the pathogen (Bruslind and Reno, 2000; Munang’andu et al., 2016; Mutoloki et al., 2016; Santi et al., 2004; Shivappa et al., 2004; Song et al., 2005). Variations in the amino acid residues in positions 217, 221 and 247 of the VP2 have been of particular focus (Dopazo, 2020; Julin et al., 2013; Munang’andu et al., 2016; Mutoloki et al., 2016; Santi et al., 2004; Song et al., 2005). Further, assessment of the recombinant viral strains has indicated that a threonine (T) residue rather than proline (P) in position 217 of the VP2 might be the most crucial amino acid that can result in the loss or gain of virulence (e.g., Song et al., 2005). Nonetheless, there has also been evidence from field outbreaks and challenge experiments that IPNV variants with a P substitute in this position can lead to high mortalities (Bain et al., 2008; Ruane et al., 2009, 2015; Tapia et al., 2015; Maj-Paluch et al., 2020).
Despite concerns and economic losses due to infection by IPNV, increased resistance against this virus through selective breeding is among the most remarkable success stories in the history of aquaculture and salmon breeding in particular. The initial attempt to estimate the genetic parameters controlling this trait, which was based on mortality in the field and controlled challenge testing, had shown that resistance to IPNV has a significant and heritable genetic component (h2 of 0.17–0.62) (Wetten et al., 2007; Kjøglum et al., 2008). However, the breakthrough came with the finding that resistance to IPNV is mainly controlled by variations in a single quantitative trait locus (QTL) on chromosome 26 (Houston et al., 2008, Houston et al., 2010; Moen et al., 2009). Subsequent work identified mutations in the epithelial cadherin gene (cdh1) as the likely causative genetic variation for IPNV resistance in Atlantic salmon (Moen et al., 2015). Following the inclusion of this QTL into the breeding programs in Norway, the industry witnessed a sharp decline in the number of IPN outbreaks throughout the country, dropping from 223 to only 19 reported cases from 2009 to 2019 (Sommerset et al., 2020).
This paper reports the identification and whole-genome sequence assembly and analysis of a new IPNV variant. The isolates were recovered from a field outbreak in the West of Norway following reported mortalities, due to an IPN infection, among QTL homozygous or heterozygous IPN resistant fish. We describe the amino acid motifs that distinguish these isolates from other phylogenetically related variants and suggest that these mutations are likely to play an essential role in this variant pathogenicity.
Materials and Methods
Field Outbreak and Sample Collection
In May of 2019, there was an IPN field outbreak in a Western region of Norway. The affected fish were Atlantic salmon of the SalmoBreed strain, delivered as unvaccinated, IPN QTL, all favorable homozygous or heterozygous for the QTL locus (Houston et al., 2008; Moen et al., 2009). Fish were at the post-smolt stage with an approximate average weight of about 700–800 gr at the time of the outbreak. The reported mortality reached about 10%, a high percentage for a field outbreak among QTL carrying fish. The IPNV infection was identified following an inspection by the veterinary authorities. The diagnosis was later validated by detection of the VP2 segment of the viral genome using polymerase chain reaction (qRT-PCR) of the infected head-kidney and through histopathological and immunohistochemical examinations of the liver and pancreatic tissues (Figures 1A,B), all performed at the Fish Vet Group Norway (http://fishvetgroup.no) (Supplementary Table S1). Tissues were fixed in a 4% formalin solution (4% formalin, 0.08 M sodium phosphate, pH 7.0), processed using Thermo Scientific Excelsior® and embedded in Histowax with the aid of Tissue—Tek®, TEC 5 (Sakura) embedding system. The embedded tissues were then sectioned at 1.5–2 µm using a Leica RM 2255 Microtome and stained with Hematoxylin-Eosin (HE). The stained sections were scanned in an Aperio ScanScope AT Turbo slide scanner and read using Aperio Image Scope (Leica Biosystems). The scans were examined and scored in a blind fashion, i.e., without information about the animal’s history. A clip from the adipose fin was also sampled and stored in 95% ethanol for DNA extraction and genotyping of ten fish, five moribund and five dead. We genotyped these fish for the genetic markers used to assign individuals as QTL carriers (Houston et al., 2008; Moen et al., 2009, 2015) and confirmed that all animals carry IPN QTL.
![www.frontiersin.org](https://www.frontiersin.org/files/Articles/635185/fgene-12-635185-HTML-r1/image_m/fgene-12-635185-g001.jpg)
FIGURE 1. (A) Histology section of the pancreatic tissue infected with the infectious pancreatic necrosis virus (IPNV). Arrows indicate acute necrosis of pancreatic acinar cells. (B) Immunohistochemical staining of the pancreas tissue using antibodies specific for IPNV. Red colours show cells infected with this virus (arrow).
Transcriptome Sequencing and Viral Genome Assembly
The head-kidney was collected from six infected fish (http://fishvetgroup.no/pcr-2/#1452163253406-d81f6dce-bdf9), all at a moribund state, and stored in RNALater® (Ambion). Total RNA extraction was performed using the RNeasy® Plus mini kit (Qiagen). Nanodrop ND-1000 Spectrophotometer (NanoDrop Technologies) was used to assess the concentration and purity of the extracted genetic materials. The Norwegian Sequencing Centre performed library construction and sequencing of the total RNA on an Illumina HiSeq 4000 platform as paired-end (PE) 150 bp reads.
We first cleaned the sequence data by removing low-quality nucleotides and sequencing adapters using Trimmomatic (v.0.36) (Bolger et al., 2014). The ribosomal RNA (rRNA) was then identified through a similarity search against the SILVA rRNA database (Quast et al., 2013) and excluded from subsequent analysis. The host-specific transcriptome was detected by aligning the remaining reads against the Atlantic salmon reference genome assembly (ICSASG_v2) (Lien et al., 2016) using TopHat (v.2.0.13) (Trapnell et al., 2009, Trapnell et al., 2012). The un-aligned reads were then blasted against the IPNV reference assemblies (ASM397166v1, ASM397170v1, ASM397174v1, ASM397168v1, ASM397172V1 and ViralMultiSegProj15024) by setting “task” as “megablast,” “word size” to 7 and “expectation value threshold” to 1.00 e−07. All the reads that aligned to the IPNV reference sequences were extracted, pooled into a single sequence datafile per animal, and then fed into Trinity (v.2.11.0) (Grabherr et al., 2011) with the default parameter settings to construct assemblies. We refer to the assembled viral genomes as BG_x_y, where x denotes the animal’s id (1–6) and y refers to the viral genome’s specific segment (e.g., VP1 or polyprotein).
Comparative Genomics and Phylogenetic Analysis
We compared the assembled viral genomes to one another (Supplementary Tables S2, S3), and against various reference isolates, with different degrees of virulence (Dopazo, 2020; Julin et al., 2013; Munang’andu et al., 2016; Mutoloki et al., 2016; Santi et al., 2004; Song et al., 2005) (Table 1) or strains representing different genogroups, as previously suggested (Blake et al., 2001) (Table 2). Further, we assessed the deduced amino acids of these isolates to any IPNV sequence in the NCBI database with a complete polyprotein or VP1 information (Data Sheets S1–S3). The sequences included 90 isolates with the complete coding sequence for the polyprotein and 120 for the VP1. Both the nucleotide and amino acid alignments were performed using ClustalX (Larkin et al., 2007) implemented in Unipro UGENE (v36.0) (Okonechnikov et al., 2012) or Muscle (Edgar, 2004) implemented in Seqotron (v1.0.1) (Fourment and Holmes, 2016). Unipro UGENE and Mega X (v10.1.7) (Kumar et al., 2018) were also used to assess sequence similarity, calculate pair-wise distances for nucleotide and amino acid data, and construct phylogenetic associations using the neighbour-joining (NJ) method. To assess the confidence of node assignments, we performed 1,000 bootstrap replications.
![www.frontiersin.org](https://www.frontiersin.org/files/Articles/635185/fgene-12-635185-HTML-r1/image_m/fgene-12-635185-t001.jpg)
TABLE 1. Amino acid sequences in segment A of the newly identified IPNV variant compared with previously reported key residues, suggested to be in association with the viral virulence (Santi et al., 2004).
![www.frontiersin.org](https://www.frontiersin.org/files/Articles/635185/fgene-12-635185-HTML-r1/image_m/fgene-12-635185-t002.jpg)
TABLE 2. Amino acid residues and their positions in the polyprotein region of the IPNV isolates reported in this study that are either exclusively unique or occur with a very low frequency in other strains.
Results and Discussion
Transcriptome Sequencing and Genome Assembly
This paper describes the genomic features of IPNV isolates responsible for a field outbreak in an Atlantic salmon population. The post-mortem examination identified IPN infection as the chief cause of mortality among all fish investigated. The histopathological and immunohistochemical examinations (Figures 1A,B; Supplementary Table S1) and the viral genetic material detected through qRT-PCR amplification (Supplementary Table S1) further confirmed earlier diagnosis. We found that the genotypes of all animals examined in this study were consistent with an IPN resistant animal (Houston et al., 2008; Moen et al., 2009; Moen et al., 2015). Further, assessment of the transcriptome sequence data confirmed that all fish carried at least one favorable copy of the causal mutation in the cdh1 gene, conferring resistance against IPNV (Moen et al., 2015). The causal mutation, which is in the form of a single nucleotide polymorphic (SNP) variation, is located in the protein-coding sequence of cdh1 (ssa26:15192533; C/T) and can change Proline (P) to Serine (S), with the former amino acid found to be associated with the IPNV resistance (Moen et al., 2015). It is expected that in a heterozygote fish, the favorable allele exhibit a close to complete dominant effect in response to the virus (Moen et al., 2015). Among the animals investigated in this study, three fish were homozygous for the favorable allele (C/C), and the other three were heterozygous (C/T). The extracted RNA was sequenced to high depth, with an average of 60 million PE reads per animal. Following the rRNA’s exclusion, more than 45 million PE reads per fish remained for the assemblage of the host’s viral genome and profiling gene expression. Together, the two segments of the IPNV genome consist of about 5,800 nucleotides. On average, we obtained 56-fold coverage of the viral genome from each animal in our sequence data.
Comparative Sequence Analysis of the Polyprotein
Comparative genomic analysis between the six assembled genomes revealed at least two different, closely linked variants, caused the field outbreak in the population investigated (Figures 2A,B; Supplementary Table S2A, S2B). The nucleotide sequences of the polyprotein are identical in four of the assemblies (BG_1_polyprot to BG_4_polyprot) and contain nine polymorphic sites compared to the other two assembled sequences (i.e., BG_5_polyprot and BG_6_polyprot). Two of the polymorphisms are nonsynonymous, changing the amino acid residues in positions 47 (aspartate (D) to glutamate (E)) and 717 (lysine (K) to glutamine (Q)).
![www.frontiersin.org](https://www.frontiersin.org/files/Articles/635185/fgene-12-635185-HTML-r1/image_m/fgene-12-635185-g002.jpg)
FIGURE 2. Neighbour-Joining tree, showing the phylogenetic relationship between the (A) Polyprotein and (B) VP5’s deduced amino acid data from different IPNV reference strains and the assemblies reported in this study. The accession ids of the sequences are presented in Table 2. The associated genogroups are in parentheses. Except for the Sp 975/99, which was recovered from the Shetland Isles, Scotland (Smail et al., 2006), all the Sp serotypes are based on different field outbreaks found in the Northwest of Norway (Santi et al., 2004). BG 1 to 6 represent isolates identified in the current work. The confidence of association between sequences was estimated using bootstrap testing (1,000 replicates) and is shown next to the branches. The evolutionary distances were computed using the Poisson correction method (Zuckerkandl and Pauling, 1965) and are in the units of the number of amino acid substitutions per site. The rate variation among sites was modeled with a gamma distribution (shape parameter = 1). All ambiguous positions were removed in a pairwise deletion fashion.
In a study aiming to trace the IPNV infection signatures from hatcheries to the sea (Kristoffersen et al., 2018), the authors have recovered a partial fragment of the polyprotein (2,772 bp) from an isolate (accession number: MH562009), with very similar nucleotide (∼99.6%) and amino acid (∼99.8–100%) sequences to the ones reported here (Supplementary Figures S1A, S1B). Recently, in controlled infection challenge testing, the mortality induced by this isolate among full-sib families was compared with a standard, virulent form of the virus (AY379740.1). While the infection with the latter variant resulted in approximately 30% mortality, mainly among the non-QTL fish, the former isolate caused more than 75% mortality, treating the QTL and non-QTL fish alike (Ørpetveit et al. in preparation; https://www.vetinst.no/nyheter/endringer-i-ipn-viruset-gjor-fisken-mer-utsatt-for-sykdom, https://www.fishfarmingexpert.com/article/changes-in-ipn-virus-make-salmon-more-susceptible, https://salmobreed.no/articles/benchmark-genetics-intensiverer-forskningen-pa-nytt-ipn-virus). Comparison between the nucleotide sequences of this isolate with our assembled genomes shows divergence ranging from 0.003 to 0.004 (Supplementary Figure S1A). However, this isolate has deduced amino acid sequences identical to the translated sequences from our first four assemblies (Supplementary Figure S1B), indicating that this variant is likely to be the ancestral form and the other isolates (i.e., BG_5_VP2 and BG_6_VP2) are more recently derived. One of the main characteristics of RNA viruses is their high mutation rates, mainly attributed to their lack of an effective proofreading activity. Therefore, finding different sequence variants from the same outbreak might not be surprising, but it indeed is important, as it will allow one to follow the trajectory of pathogen evolution through time and space.
Comparative assessment of the polyprotein nucleotide (data not shown) as well as the amino acid sequence data with multiple reference strains (Figure 2A) shows our assemblies are most closely related to the Sp serotype (genogroup 5). However, they also constitute a distinct group within this serotype compared to Sp isolates assessed in this study (Figure 2A; Supplementary Figure S2A). Several investigations have previously reported variations in some amino acid residues, particularly in the VP2 of the polyprotein, to generally be associated with a higher degree of the IPNV virulence (Bruslind and Reno, 2000; Munang’andu et al., 2016; Mutoloki et al., 2016; Santi et al., 2004; Shivappa et al., 2004; Song et al., 2005). Interestingly, the amino acid motifs found in the six assembled genomes in this study are mainly consistent with the virus’s avirulent form (Table 1) (Santi et al., 2004; Shivappa et al., 2004; Song et al., 2005). For instance, for the most critical positions associated with the virulence, i.e., 217 and 221 of the VP2, the amino acid residues of our assembled genomes are proline (P) and threonine (T) (Table 1; Supplementary File S1), while it is expected that in the virulent form, the corresponding residues to be T and A, respectively (Santi et al., 2004; Song et al., 2005; Mutoloki et al., 2016). Similar patterns have also been reported from the isolates recovered from outbreaks, with the amino acid signature comparable to the non-virulent forms in the key positions (Bain et al., 2008; Maj-Paluch et al., 2020; Ruane et al., 2009, Ruane et al., 2015; Tapia et al., 2015). On the other hand, we identified ten amino acid residues in the assembled genomes’ polyprotein that are either exclusively or with a disproportionate frequency in the reported isolates (Table 2; Supplementary File S1). Nine of these sites are within the VP2, and one is in the VP4. In the VP2, these residues are all close to one another, encompassing the hypervariable region (Blake et al., 2001; Heppell et al., 1995b). Specifically, these residues are in positions 245, 248, 252, 255, 257, 278, 282, 285, 321 (Table 2; Supplementary File S1). This data provides additional support for the previously suggested notion that the amino acid motifs in the VP2 capsid protein hypervariable region might play an essential role in the degree of the virulence displayed between different IPNV isolates (Bruslind and Reno, 2000; Skjesol et al., 2011). It is also likely that this specific constellation of the amino acid residues reported in these isolates had helped the virus to escape the host defence barrier, despite fish being regarded as genetically resistant to IPNV. In the VP4, the deduced amino acid residue, unique to the assembled isolates (Table 2; Supplementary File S1), is located in position 585. The amino acid substitution is a glycine (G) residue, compared to the other sequences investigated in this study, where they carry an aspartate (D) (Table 2).
While this paper was under review, another article also reported similar variations from isolates recovered from outbreaks in Scotland (Benkaroun et al., 2021). Taken together, these results confirm the virulence of this variant and provide some indication of the geographical regions it has spread.
Comparative Sequence Analysis of VP5
All assemblies contain the VP5, the partial overlapping, non-structural protein, of segment A. Similar to the polyprotein, the phylogenetic analysis of the VP5 grouped the assemblies into two distinct groups and further clustered them with the isolates from the Sp serotype (Figure 2B). The two assembled groups differ in only one amino acid residue at position 90, with the first four assemblies carrying an arginine (R) while the two other assemblies have a glycine (G) (Supplementary File S3). Sequence analysis of the deduced amino acids of VP5 did not reveal any unique residue in the assembled isolates compared to all other sequences. However, the amino acid in position 98 differs between the Sp isolates (Supplementary File S3). At this position, all the investigated Sp isolates carry an R, while in the assemblies, the residue is tryptophan (W). It has been shown that the VP5’s protein is generally produced during the initial stages of replication (Heppell et al., 1995a), and at least in the Ab serotype, it seems that this protein to be involved in anti-apoptotic functions (Hong and Wu, 2002). However, this product’s importance in viral growth or virulence has so far remained ambiguous (e.g., Julin et al., 2013; Dopazo 2020).
Comparative Sequence Analysis of VP1
The phylogenetic analysis of VP1 clustered the assemblies into three distinct clades (Supplementary Figure S2B) with the amino acid variations in positions 119 (threonine (T) to methionine (M)) and 701 (T to alanine (A) and serine (S)) (Supplementary File S2; Supplementary Table S3A). The nucleotide sequence comparisons show a greater divergence rate in the VP1 of the assemblies than those of the polyprotein, suggesting that the B segment of the genome has a higher rate of nucleotide mutation (Supplementary Tables S1, S2). Comparing the deduced amino acid sequences of the assembled genomes to the VP1 section of 120 IPNV sequences further revealed only a single, almost unique residue in our isolates (Supplementary File S2). The amino acid residue in position 656 of the assembled isolates is glutamate (E), while in the majority of other strains, the amino acid is either aspartate (D) or alanine (A). The only other isolates, carrying an E in this position, are two sequences recovered from Atlantic salmon (KY548520; collection date: 2001) and Brown trout (Salmo trutta L.) (KY548519; collection date: 1989) in Finland following outbreaks (Holopainen et al., 2017). However, it should be noted that so far, the role of VP1 in determining the degree of IPNV virulence has remained ambiguous. For example, in a study conducted by Song et al. 2005, through construction and comparative assessment of chimeric IPNV strains, the authors suggested that VP1 is not involved in this pathogen’s virulence. On the other hand, the data reported by Shivappa et al. 2004 suggests that the amino acids on the B segment of the genome in combination with specific residues in VP2 might help modulate the magnitude of virulence. Work in the infectious bursal disease virus (IBDV) has also shown that variations in the amino acid residues of the VP1 can change the kinetics of viral replication and influence the virus’s virulence (Liu and Vakharia, 2004). VP1 interacts with VP3 to form VP1-VP3 complexes, and as such, can affect various aspects of the virus biology, including its replication efficiency (Tacken et al., 2002; Pedersen et al., 2007).
Concluding Remarks
In light of the findings reported in this study and similar works, a key point worth noting is the importance of understanding the genetic interactions between the host and the pathogen in determining an infection’s outcome. As a general trend, research studies tend to mainly focus on either the host or the pathogen’s genetic polymorphisms while investigating the genetic basis of disease-related phenotypic variations. While the merits of such a targeted approach cannot be disputed, it is becoming increasingly important to be aware of the detailed genetic and genomic variations and interactions between the host and the microbe. In the case of monogenic infectious diseases, it is relatively trivial to speculate whether a new form of a pathogen interacts in fundamentally different ways with the host’s defence mechanism. For instance, considering that variations in a single gene mainly determine resistance to IPNV in Atlantic salmon, one might expect that the mutations acquired by the variant reported in this paper might have offered an adaptive advantage, assisting the pathogen to establish successful infection, irrespective of the variations in the cdh1 gene. This, for example, can be examined by challenging and comparative transcriptome analysis of the QTL and non-QTL full-sib animals following infection with any of these isolates and the common variant in a controlled setting. Of course, pinpointing the causative mutations and understanding and validating the molecular basis of such adaptive variations is a substantial undertaking by itself. The task of identifying and untangling the dynamics of how genotypic variations in the host and the pathogen influence one another and affect the outcome of infection becomes more challenging while investigating diseases with polygenetic nature. This point becomes even more critical, considering our rapidly changing environments that can facilitate the transmission dynamics and the geographical spread of the pathogens (Reid et al., 2019). Fortunately, advancements in technologies such as sequencing, genotyping, computational analysis and molecular biology can now provide us with many necessary tools if we are set to understand the dynamics of the interactions between the host and the microbe at their detailed molecular levels.
Data Availability Statement
All datasets generated for this study are included in the article (Supplementary Material) and have also been deposited in GenBank under the accession numbers MW496366, MW496367, MW496368, MW496369, MW496370, MW496371, MW496372, MW496373, MW496374, MW496375, MW496376 and MW496377.
Ethics Statement
Ethical review and approval was not required for the animal study as the data is based on a naturally occurring disease outbreak in a field.
Author Contributions
GM, SJ and BH arranged for tissue sampling, shipment and all laboratory and sequencing procedures. HM performed analyses of the sequence data. BH and HM drafted the article. All authors contributed to the development of the paper, interpretation of the results and improvement of the article.
Conflict of Interest
Author BH, SJ, GM, HM were employed by company Benchmark Genetics Norway AS.
Publisher’s Note
All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors, and the reviewers. Any product that may be evaluated in this article, or claim that may be made by its manufacturer, is not guaranteed or endorsed by the publisher.
Acknowledgments
We would like to thank the Fish Vet Group Norway, especially Dr Marianne Kraugerud and Simon Rey, for the laboratory work and the Norwegian Sequencing Centre for excellent work and guidance.
Supplementary Material
The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fgene.2021.635185/full#supplementary-material
References
Bain, N., Gregory, A., and Raynard, R. S. (2008). Genetic analysis of infectious pancreatic necrosis virus from Scotland. J. Fish. Dis. 31, 37–47. doi:10.1111/j.1365-2761.2007.00864.x
Benkaroun, J., Muir, K. F., Allshire, R., Tamer, C., and Weidmann, M. (2021). Isolation of a New Infectious Pancreatic Necrosis Virus (IPNV) Variant from a Fish Farm in Scotland. Viruses, 13. doi:10.3390/v13030385
Blake, S., Ma, J.-Y., Caporale, D., Jairath, S., and Nicholson, B. (2001). Phylogenetic relationships of aquatic birnaviruses based on deduced amino acid sequences of genome segment A cDNA. Dis. Aquat. Org. 45, 89–102. doi:10.3354/dao045089
Bolger, A. M., Lohse, M., and Usadel, B. (2014). Trimmomatic: A flexible trimmer for Illumina sequence data. Bioinformatics 30, 2114–2120. doi:10.1093/bioinformatics/btu170
Bruslind, L., and Reno, P. (2000). Virulence comparison of three Buhl-subtype isolates of infectious pancreatic necrosis virus in brook trout fry. J. Aquat. Anim. Hlth. 12, 301–315. doi:10.1577/1548-8667(2000)012<0301:vcotbs>2.0.co;2
Dobos, P., and Roberts, T. E. (1983). The molecular biology of infectious pancreatic necrosis virus: A review. Can. J. Microbiol. 29, 377–384. doi:10.1139/m83-062
Dopazo, C. P. (2020). The infectious pancreatic necrosis virus (IPNV) and its virulence determinants: What is known and what should be known. Pathogens 9, 94. doi:10.3390/pathogens9020094
Dorson, M., Castric, J., and Torchy, C. (1978). Infectious pancreatic necrosis virus of salmonids: biological and antigenic features of a pathogenic strain and of a non-pathogenic variant selected in RTG-2 cells. J. Fish. Dis. 1, 309–320. doi:10.1111/j.1365-2761.1978.tb00035.x
Duncan, R., and Dobos, P. (1986). The nudeotide sequence of infectious pancreatic necrosis virus (IPNV) dsRNA segment a reveals one large orf encoding a precursor polyprotein. Nucl. Acids Res. 14, 5934. doi:10.1093/nar/14.14.5934
Duncan, R., Nagy, E., Krell, P. J., and Dobos, P. (1987). 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. 61, 3655–3664. doi:10.1128/jvi.61.12.3655-3664.1987
Edgar, R. C. (2004). MUSCLE: Multiple sequence alignment with high accuracy and high throughput. Nucleic Acids Res. 32, 1792–1797. doi:10.1093/nar/gkh340
Ellis, A. E., Cavaco, A., Petrie, A., Lockhart, K., Snow, M., and Collet, B. (2010). Histology, immunocytochemistry and qRT-PCR analysis of Atlantic salmon, Salmo salar L., post-smolts following infection with infectious pancreatic necrosis virus (IPNV). J. Fish. Dis. 33, 803–818. doi:10.1111/j.1365-2761.2010.01174.x
Fourment, M., and Holmes, E. C. (2016). Seqotron: A user-friendly sequence editor for Mac OS X. BMC Res. Notes 9, 106. doi:10.1186/s13104-016-1927-4
Grabherr, M. G., Haas, B. J., Yassour, M., Levin, J. Z., Thompson, D. A., Amit, I., et al. (2011). Full-length transcriptome assembly from RNA-Seq data without a reference genome. Nat. Biotechnol. 29, 644–652. doi:10.1038/nbt.1883
Guy, D. R., Bishop, S. C., Brotherstone, S., Hamilton, A., Roberts, R. J., McAndrew, B. J., et al. (2006). Analysis of the incidence of infectious pancreatic necrosis mortality in pedigreed Atlantic salmon, Salmo salar L., populations. J. Fish. Dis. 29, 637–647. doi:10.1111/j.1365-2761.2006.00758.x
Håstein, T., and Krogsrud, J. (1976). Infectious Pancreatic Necrosis First Isolation of Virus from Fish in Norway. Acta Vet. Scand. 17, 109–111. doi:10.1186/bf03547948
Heppell, J., Tarrab, E., Berthiaume, L., Lecomte, J., and Arella, M. (1995a). Characterization of the small open reading frame on genome segment A of infectious pancreatic necrosis virus. J. Gen. Virol. 76, 2091–2096. doi:10.1099/0022-1317-76-8-2091
Heppell, J., Tarrab, E., Lecomte, J., Berthiaume, L., and Arella, M. (1995b). Strain Variability and Localization of Important Epitopes on the Major Structural Protein (VP2) of Infectious Pancreatic Necrosis Virus1. Virology 214, 40–49. doi:10.1006/viro.1995.9956
Hill, B. J., and Way, K. (1995). Serological classification of infectious pancreatic necrosis (IPN) virus and other aquatic birnaviruses. Annu. Rev. Fish Dis. 5, 55–77. doi:10.1016/0959-8030(95)00011-9
Holopainen, R., Eriksson-Kallio, A. M., and Gadd, T. (2017). Molecular characterisation of infectious pancreatic necrosis viruses isolated from farmed fish in Finland. Arch. Virol. 162, 3459–3471. doi:10.1007/s00705-017-3525-8
Hong, J.-R., and Wu, J.-L. (2002). Induction of apoptotic death in cells via Bad gene expression by infectious pancreatic necrosis virus infection. Cell Death Differ 9, 113–124. doi:10.1038/sj.cdd.4400933
Houston, R. D., Haley, C. S., Hamilton, A., Guy, D. R., Mota-Velasco, J. C., Gheyas, A. A., et al. (2010). The susceptibility of Atlantic salmon fry to freshwater infectious pancreatic necrosis is largely explained by a major QTL. Heredity 105, 318–327. doi:10.1038/hdy.2009.171
Houston, R. D., Haley, C. S., Hamilton, A., Guy, D. R., Tinch, A. E., Taggart, J. B., et al. (2008). Major quantitative trait loci affect resistance to infectious pancreatic necrosis in Atlantic salmon (Salmo salar). Genetics 178, 1109–1115. doi:10.1534/genetics.107.082974
Jarp, J., Gjevre, A. G., Olsen, A. B., and Bruheim, T. (1995). Risk factors for furunculosis, infectious pancreatic necrosis and mortality in post-smolt of Atlantic salmon, Salmo solar L. J. Fish. Dis. 18, 67–78. doi:10.1111/j.1365-2761.1995.tb01267.x
Julin, K., Mennen, S., and Sommer, A.-I. (2013). Study of virulence in field isolates of infectious pancreatic necrosis virus obtained from the northern part of Norway. J. Fish. Dis. 36, 89–102. doi:10.1111/j.1365-2761.2012.01423.x
Kjøglum, S., Henryon, M., Aasmundstad, T., and Korsgaard, I. (2008). Selective breeding can increase resistance of Atlantic salmon to furunculosis, infectious salmon anaemia and infectious pancreatic necrosis. Aquac. Res. 39, 498–505.
Kristoffersen, A. B., Devold, M., Aspehaug, V., Gjelstenli, O., Breck, O., and Bang Jensen, B. (2018). Molecular tracing confirms that infection with infectious pancreatic necrosis virus follows the smolt from hatchery to grow-out farm. J. Fish. Dis. 41, 1601–1607. doi:10.1111/jfd.12844
Kumar, S., Stecher, G., Li, M., Knyaz, C., and Tamura, K. (2018). MEGA X: Molecular evolutionary genetics analysis across computing platforms. Mol. Biol. Evol. 35, 1547–1549. doi:10.1093/molbev/msy096
Larkin, M. A., Blackshields, G., Brown, N. P., Chenna, R., Mcgettigan, P. A., McWilliam, H., et al. (2007). Clustal W and Clustal X version 2.0. Bioinformatics 23, 2947–2948. doi:10.1093/bioinformatics/btm404
Lien, S., Koop, B. F., Sandve, S. R., Miller, J. R., Kent, M. P., Nome, T., et al. (2016). The Atlantic salmon genome provides insights into rediploidization. Nature 533, 200–205. doi:10.1038/nature17164
Liu, M., and Vakharia, V. N. (2004). VP1 protein of infectious bursal disease virus modulates the virulence in vivo. Virology 330, 62–73. doi:10.1016/j.virol.2004.09.009
Maj-Paluch, J., Matras, M., Borzym, E., Stachnik, M., and Reichert, M. (2020). Phylogenetic characterization of Polish isolates of infectious pancreatic necrosis virus in salmonid fish. J. Fish. Dis. 43, 1443–1451. doi:10.1111/jfd.13249
M’Gonigle, R. H. (1941). Acute Catarrhal Enteritis of Salmonid Fingerlings. Trans. Am. Fish. Soc. 70, 297–303.
Moen, T., Baranski, M., Sonesson, A. K., and Kjøglum, S. (2009). Confirmation and fine-mapping of a major QTL for resistance to infectious pancreatic necrosis in Atlantic salmon (Salmo salar): Population-level associations between markers and trait. BMC Genomics 10, 368. doi:10.1186/1471-2164-10-368
Moen, T., Torgersen, J., Santi, N., Davidson, W. S., Baranski, M., Ødegård, J., et al. (2015). Epithelial cadherin determines resistance to infectious pancreatic necrosis virus in Atlantic salmon. Genetics 200, 1313–1326. doi:10.1534/genetics.115.175406
Munang'andu, H. M., Fredriksen, B. N., Mutoloki, S., Dalmo, R. A., and Evensen, Ø. (2013). Antigen dose and humoral immune response correspond with protection for inactivated infectious pancreatic necrosis virus vaccines in Atlantic salmon (Salmo salar L). Vet. Res. 44, 7. doi:10.1186/1297-9716-44-7
Munang’andu, H. M., Santi, N., Fredriksen, B. N., Løkling, K. E., and Evensen, Ø. (2016). A systematic approach towards optimizing a cohabitation challenge model for infectious pancreatic necrosis virus in Atlantic salmon (Salmo salar L). PLoS One 11, e0148467.
Mutoloki, S., and Evensen, Ø. (2011). Sequence similarities of the capsid gene of Chilean and European isolates of infectious pancreatic necrosis virus point towards a common origin. J. Gen. Virol. 92, 1721–1726. doi:10.1099/vir.0.030270-0
Mutoloki, S., Jøssund, T. B., Ritchie, G., Munang'andu, H. M., and Evensen, Ø. (2016). Infectious pancreatic necrosis virus causing clinical and subclinical infections in atlantic salmon have different genetic fingerprints. Front. Microbiol. 7, 1393. doi:10.3389/fmicb.2016.01393
Nishizawa, T., Kinoshita, S., and Yoshimizu, M. (2005). An approach for genogrouping of Japanese isolates of aquabirnaviruses in a new genogroup, VII, based on the VP2/NS junction region. J. Gen. Virol. 86, 1973–1978. doi:10.1099/vir.0.80438-0
Okonechnikov, K., Golosova, O., Fursov, M., Varlamov, A., Vaskin, Y., Efremov, I., et al. (2012). Unipro UGENE: A unified bioinformatics toolkit. Bioinformatics 28, 1166–1167. doi:10.1093/bioinformatics/bts091
Pedersen, T., Skjesol, A., and Jørgensen, J. B. (2007). VP3, a Structural Protein of Infectious Pancreatic Necrosis Virus, Interacts with RNA-dependent RNA Polymerase VP1 and with Double-Stranded RNA. J. Virol. 81, 6652–6663. doi:10.1128/jvi.02831-06
Quast, C., Pruesse, E., Yilmaz, P., Gerken, J., Schweer, T., Yarza, P., et al. (2013). The SILVA ribosomal RNA gene database project: Improved data processing and web-based tools. Nucleic Acids Res. 41, D590–D596. doi:10.1093/nar/gks1219
Reid, G., Gurney-Smith, H., Marcogliese, D., Knowler, D., Benfey, T., Garber, A., et al. (2019). Climate change and aquaculture: Considering biological response and resources. Aquacult. Environ. Interact. 11, 569–602. doi:10.3354/aei00332
Roberts, R. J., and Pearson, M. D. (2005). Infectious pancreatic necrosis in Atlantic salmon,Salmo salarL. J. Fish. Dis. 28, 383–390. doi:10.1111/j.1365-2761.2005.00642.x
Ruane, N. M., McCarthy, L. J., Swords, D., and Henshilwood, K. (2009). Molecular differentiation of infectious pancreatic necrosis virus isolates from farmed and wild salmonids in Ireland. J. Fish. Dis. 32, 979–987. doi:10.1111/j.1365-2761.2009.01080.x
Ruane, N. M., McCleary, S. J., McCarthy, L. J., and Henshilwood, K. (2015). Phylogenetic analysis of infectious pancreatic necrosis virus in Ireland reveals the spread of a virulent genogroup 5 subtype previously associated with imports. Arch. Virol. 160, 817–824. doi:10.1007/s00705-014-2307-9
Santi, N., Song, H., Vakharia, V. N., and Evensen, Ø. (2005). Infectious Pancreatic Necrosis Virus VP5 Is dispensable for virulence and Persistence. J. Virol. 79, 9206–9216. doi:10.1128/jvi.79.14.9206-9216.2005
Santi, N., Vakharia, V. N., and Evensen, Ø. (2004). Identification of putative motifs involved in the virulence of infectious pancreatic necrosis virus. Virology 322, 31–40. doi:10.1016/j.virol.2003.12.016
Shivappa, R., Song, H., Yao, K., Aas-Eng, A., Evensen, Ø., and Vakharia, V. (2004). Molecular characterization of Sp serotype strains of infectious pancreatic necrosis virus exhibiting differences in virulence. Dis. Aquat. Org. 61, 23–32. doi:10.3354/dao061023
Skjesol, A., Skjæveland, I., Elnæs, M., Timmerhaus, G., Fredriksen, B. N., Jørgensen, S. M., et al. (2011). IPNV with high and low virulence: Host immune responses and viral mutations during infection. Virol. J. 8, 396. doi:10.1186/1743-422x-8-396
Smail, D. A., Bain, N., Bruno, D. W., King, J. A., Thompson, F., Pendrey, D. J., et al. (2006). Infectious pancreatic necrosis virus in Atlantic salmon, Salmo salar L., post-smolts in the Shetland Isles, Scotland: Virus identification, histopathology, immunohistochemistry and genetic comparison with Scottish mainland isolates. J. Fish. Dis. 29, 31–41. doi:10.1111/j.1365-2761.2005.00678.x
Sommerset, I., Walde, C. S., Bang Jensen, B., Bornø, B., Haukaas, A., and Brun, E. (2020). Fiskehelserapporten 2019. Norway.
Song, H., Santi, N., Evensen, Ø., and Vakharia, V. N. (2005). Molecular determinants of Infectious Pancreatic Necrosis Virus virulence and cell culture adaptation. J. Virol. 79, 10289–10299. doi:10.1128/jvi.79.16.10289-10299.2005
Sundh, H., Olsen, R.-E., Fridell, F., Gadan, K., Evensen, O., Glette, J., et al. (2009). 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 salarL. J. Fish. Dis. 32, 687–698. doi:10.1111/j.1365-2761.2009.01047.x
Tacken, M. G., Peeters, B. P., Thomas, A. A., Rottier, P. J., and Boot, H. J. (2002). Infectious Bursal Disease Virus Capsid Protein VP3 Interacts both with VP1, the RNA-dependent RNA Polymerase, and with Viral Double-Stranded RNA. J. Virol. 76, 11301–11311. doi:10.1128/jvi.76.22.11301-11311.2002
Tapia, D., Eissler, Y., Torres, P., Jorquera, E., Espinoza, J., and Kuznar, J. (2015). Detection and phylogenetic analysis of infectious pancreatic necrosis virus in Chile. Dis. Aquat. Org. 116, 173–184. doi:10.3354/dao02912
Trapnell, C., Pachter, L., and Salzberg, S. L. (2009). TopHat: Discovering splice junctions with RNA-Seq. Bioinformatics 25, 1105–1111. doi:10.1093/bioinformatics/btp120
Trapnell, C., Roberts, A., Goff, L., Pertea, G., Kim, D., Kelley, D. R., et al. (2012). Differential gene and transcript expression analysis of RNA-seq experiments with TopHat and Cufflinks. Nat. Protoc. 7, 562–578. doi:10.1038/nprot.2012.016
Wetten, M., Aasmundstad, T., Kjøglum, S., and Storset, A. (2007). Genetic analysis of resistance to infectious pancreatic necrosis in Atlantic salmon (Salmo salar L.). Aquaculture 272, 111–117. doi:10.1016/j.aquaculture.2007.08.046
Wolf, K., Snieszko, S. F., Dunbar, C. E., and Pyle, E. (1960). Virus Nature of Infectious Pancreatic Necrosis in Trout. Proc. Soc. Exp. Biol. Med. 104, 25743. doi:10.3181/00379727-104-25743
Keywords: RNA sequencing, genome assembly, quantitative trait locus (QTL), Atlantic salmon, infectious pancreatic necrosis virus (IPNV)
Citation: Hillestad B, Johannessen S, Melingen GO and Moghadam HK (2021) Identification of a New Infectious Pancreatic Necrosis Virus (IPNV) Variant in Atlantic Salmon (Salmo salar L.) that can Cause High Mortality Even in Genetically Resistant Fish. Front. Genet. 12:635185. doi: 10.3389/fgene.2021.635185
Received: 16 December 2020; Accepted: 18 October 2021;
Published: 26 November 2021.
Edited by:
Jiuzhou Song, University of Maryland, United StatesReviewed by:
Oystein Evensen, Norwegian University of Life Sciences, NorwayYoanna Eissler, Universidad de Valparaíso, Chile
Ruben Avendaño-Herrera, Andres Bello University, Chile
Copyright © 2021 Hillestad, Johannessen, Melingen and Moghadam. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.
*Correspondence: Hooman K. Moghadam, hooman.moghadam@bmkgenetics.com