Whole-genome analysis of piscine reovirus (PRV) shows PRV represents a new genus in family Reoviridae and its genome segment S1 sequences group it into two separate sub-genotypes
Virology Journal volume 10, Article number: 230 (2013)
Piscine reovirus (PRV) is a newly discovered fish reovirus of anadromous and marine fish ubiquitous among fish in Norwegian salmon farms, and likely the causative agent of heart and skeletal muscle inflammation (HSMI). HSMI is an increasingly economically significant disease in Atlantic salmon (Salmo salar) farms. The nucleotide sequence data available for PRV are limited, and there is no genetic information on this virus outside of Norway and none from wild fish.
RT-PCR amplification and sequencing were used to obtain the complete viral genome of PRV (10 segments) from western Canada and Chile. The genetic diversity among the PRV strains and their relationship to Norwegian PRV isolates were determined by phylogenetic analyses and sequence identity comparisons.
PRV is distantly related to members of the genera Orthoreovirus and Aquareovirus and an unambiguous new genus within the family Reoviridae. The Canadian and Norwegian PRV strains are most divergent in the segment S1 and S4 encoded proteins. Phylogenetic analysis of PRV S1 sequences, for which the largest number of complete sequences from different “isolates” is available, grouped Norwegian PRV strains into a single genotype, Genotype I, with sub-genotypes, Ia and Ib. The Canadian PRV strains matched sub-genotype Ia and Chilean PRV strains matched sub-genotype Ib.
PRV should be considered as a member of a new genus within the family Reoviridae with two major Norwegian sub-genotypes. The Canadian PRV diverged from Norwegian sub-genotype Ia around 2007 ± 1, whereas the Chilean PRV diverged from Norwegian sub-genotype Ib around 2008 ± 1.
The newly discovered piscine reovirus (PRV) belongs to the family Reoviridae, subfamily Spinareovirinae, probably in a new reovirus genus that is equally distant to the genera Orthoreovirus and Aquareovirus, although with 10 genome segments, PRV is like members of the genus Orthoreovirus and unlike the genus Aquareovirus with 11 segments. The Orthoreovirus genus can be divided into the fusogenic and non-fusogenic orthoreoviruses based on the ability of the fusogenic orthoreoviruses to induce cell-cell fusion during infection resulting in syncytium formation  by virtue of possession of a fusion-associated small transmembrane (FAST) protein . Whereas the non-fusogenic orthoreoviruses, Mammalian Orthoreovirus (MRV), are not clinically significant , the fusogenic orthoreoviruses Nelson Bay virus (NBV)  and Baboon Orthoreovirus (BRV)  that infect primates, Avian Orthoreovirus (ARV)  that infect birds, and Reptilian Orthoreovirus (RRV)  that infect reptiles, have been shown to cause significant and often fatal disease. Most recently, PRV has been shown to be more closely related with recognized orthoreoviruses than with recognized aquareoviruses, and does not encode a FAST protein and is therefore non-fusogenic .
PRV is associated with heart and skeletal muscle inflammation (HSMI) ; an emerging disease of marine-farmed Atlantic salmon , first recognized in 1999 in western Norway  and subsequently in Scotland . PRV has also been detected by real-time reverse transcription quantitative polymerase chain reaction (RT-qPCR) at a low prevalence in wild Atlantic salmon “S. salar”  and in certain marine fish species (Atlantic herring “Clupea harengus”, Capelin “Mallotus villosus”, Atlantic horse mackerel “Trachurus trachurus”, and Great silver smelt “Argentina silus”) along the coast of Norway . PRV was also detected in 3% of anadromous trout (sea-trout) “Salmo trutta” tested, but not in anadromous Arctic char “Salvelinus alpinus” . PRV is ubiquitous in Norwegian salmon farms [16, 17], but there is a significant increase in the viral load and tissue distribution during outbreaks of HSMI [2, 18]. The virus can be propagated in the GF-1 cell line , derived from the tissue of orange-spotted grouper, Epinephelus coioides, and cardiac and skeletal muscle pathology typical of HSMI can be reproduced in naïve Atlantic salmon by experimental inoculation with the supernatant from cell culture passaged PRV . Most recently, it has been reported that serum enzymes creatine kinase and lactate dehydrogenase, associated with cardiac injury in humans , are significantly correlated with HSMI histopathology in Atlantic salmon . Other reports doubt the pathogenicity of PRV, describing PRV as an opportunistic virus [23, 24] or non-pathogenic virus . The virus has been detected in marine-farmed Atlantic salmon in Chile [25, 26]. There is anecdotal evidence that it is also present in farmed Atlantic salmon and wild Pacific salmon in British Columbia-Canada , where 75% of 300 farm salmon reportedly tested positive for PRV  but no sequence information was reported.
The PRV genome comprises at least 10 dsRNA segments including three large (L), three medium (M), and four small (S) size-class RNA genome segments . To date only two “isolates”, both from marine-farmed Atlantic salmon from Norway, have been sequenced on all 10 genomic segments by high-throughput pyrosequencing of clinical samples: Reovirus sp. Salmo/GP-2010/NOR from HSMI , and CMS PRV from a CMS outbreak . However, only the Salmo/GP-2010/NOR sequences are accessible from the GenBank Database (GenBank accession numbers GU994013-GU994022). The coding assignments of genomic segments S1 and S4 initially reported to encode proteins with no identified homologs in orthoreoviruses and aquareoviruses , were recently shown to be reversed such that S1 encodes the major outer capsid protein, Outer clamp protein (σ3 and VP7 in MRV and aquareoviruses, respectively), and S4 encodes the virus attachment protein, Outer fiber protein (σ1 in MRV, which is absent in aquareoviruses) . A further complication is that several sequences of Norwegian PRV isolates were deposited in the GenBank database as S4 sequences [18, 28] but correspond to S1 sequences [2, this study], and sequences of the remaining 9 genomic segments for these virus isolates have not been reported.
The sequence data available for PRV strains are limited, with no genetic information on this virus outside of Norway, and none from wild fish despite the economical impact of HSMI on salmon aquaculture and the potential for transmission of PRV to wild salmon populations or from wild salmon to farmed salmon.
The primary goal of the present study was to determine the genetic diversity among PRV strains detected in tissue samples obtained from fish in western Canada, and in Chile, and their relationship to known Norwegian PRV sequences. We also attempted to sequence the complete genomes of three “isolates”, two Canadian and one Chilean to obtain more information about the taxonomic assignment of PRV.
Results and discussion
Amplification and sequencing of cDNA of genomic segments of PRV from fish samples
Piscine reovirus was readily detected by RT-qPCR during testing at the Atlantic Veterinary College laboratory in fish tissue samples from western Canada, and at the ETECMA diagnostic laboratory in fish tissue samples from Chile (data not shown). Consistent with observations elsewhere [2, 18], PRV is ubiquitous in marine-farmed salmon. PRV was consistently detected in gill tissue, identifying the gills as suitable target tissue and likely a primary transmission route for PRV. This is consistent with Orthoreovirus, which are commonly isolated from enteric and respiratory tract tissues .
Ten samples from western Canada with either low Ct values or unique case histories, host species, and sampling times listed in (Additional file 1: Table S1a) were selected for amplification and cloning of cDNA of viral genome segments. Four additional samples for which the 3′ portion of genome segment L1 had been PCR-amplified during the original testing for PRV (Table 1), were included in the analysis of PRV sequences.
Figure 1 shows the RT-PCR amplification and sequencing strategy used for the PRV genome segments, based on Canadian “isolate” 358. The new PRV nucleotide sequences (Additional file 2: Table S2) are available through the GenBank database . The complete PRV genome (10 segments) was amplified from 2 of 10 samples. Additional partial or full-length sequences were also obtained on PRV genomic segments L1 (12 samples), L2 (1 sample), L3 (2 samples), M1 (4 samples), M2 (2 samples), M3 (4 samples), S1 (4 samples), S2 (3 samples), S3 (1 sample), and S4 (1 sample). PRV sequences were obtained from four different western Canada fish species (Atlantic salmon “Salmo salar”, Cutthroat trout “Oncorhynchus clarkii”, Steelhead trout “Oncorhynchus mykiss”, and Chum salmon “Oncorhynchus keta”) (Table 1 and Additional file 2: Table S2). Failure to amplify transcripts from all PRV positive samples was attributed to variation in viral loads. It has been reported that fish are capable of reducing the viral load by the end of the production cycle . The differences in RT-PCR amplification could be due to differences between the PRV “isolates”. It is also possible that variations in transcription levels of different virus genes, and efficiency of PCR of the different targets contributed to the inability to amplify all 10 genome segment transcripts in some of the samples.
Among the Chilean PRV positive samples, 6 fish individually sampled from the same farm with low Ct values were selected for amplification of all 10 viral genome segments; sequences from one of these samples was used in the analysis. These were all fish kidney samples, which had significantly lower Ct values (Additional file 1: Table S1b) compared to the fish gill samples from Canada (Additional file 1: Table S1a). Three additional sequences on PRV genomic segments L3 (3 samples) and S1 (2 samples) were already available and were included in the analysis, for a total of 4 Chilean PRV “isolates” (Table 1 and Additional file 2: Table S2).
(Additional file 2: Table S2) shows the nucleotide and amino acid sequence identities of the new PRV isolates when compared to the single Norwegian PRV isolate, Salmo/GP-2010/NOR (GenBank accession numbers GU994013-GU994022). The largest nucleotide sequence differences between Canadian and Norwegian PRV strains are on segments M2 and S1 (96-97% sequence identity). However, at the amino acid sequence level, our analysis shows that Canadian and Norwegian PRV strains are most divergent in the S1 encoded proteins, the major outer capsid protein (Outer clamp protein) and the non-structural protein p13, and the S4 encoded virus attachment protein (Outer fiber) . The difference on the S4 protein is very interesting as it consists of a variable region of 18 residues at the C-terminus. This work is the first report of genomic analysis of PRV strains detected in tissue samples obtained from fish outside of Norway, extending the current geographical range of the characterized virus to both North and South America.
The PRV conserved terminal nucleotide sequences
Conserved terminal nucleotide sequences are useful for reovirus classification . Palacios et al. reported the complete genome sequence of the Norwegian PRV isolate Salmo/GP-2010/NOR including the conserved nucleotides at the 5′ end and the 3′ end of the genome (5′-GAUAAA/U------UCAUC-3). Table 2 compares these conserved terminal sequences to those of members of the Orthoreovirus and Aquareovirus genera. The conserved nucleotides 5′-GAUAAA/U were present at the 5′ ends in all the positive strands of each of the 10 genome segments of PRV, and are unique to PRV, whereas the 3′ conserved termini UCAUC-3′ are also conserved between PRV, and the Orthoreovirus and Aquareovirus genera (Table 2).
The PRV protein profile deduced from whole-genome sequence analysis
In the present study, the major open reading frames (ORFs) in the 10 PRV genomic segments, identified based on the first methionine of the ORF, vary in length from 315 codons in S4 to 1,290 codons in L2. The lengths of the non-coding regions ranged from 7 to 83 nucleotides at the 5′ end and from 44 to 89 at the 3′ end. The putative PRV gene products calculated from the nucleotide sequence data in this study are shown in Table 3. Only the S1 genome segment is bicistronic, encoding the Outer clamp protein and a nonstructural protein, p13, which is not a FAST protein . In this sense, PRV is similar to Mammalian orthoreovirus (MRV), which also does not have a FAST protein and is non-fusogenic. However, MRV differs from PRV in gene coding assignments: in MRV, Core RdRp (λ3) is encoded on segment L1; Core shell (λ1) on segment L3; Outer clamp (σ3) is encoded on S4; Outer fiber (σ1) and NS, other (σ1s) are encoded on S1 (Table 3). Probably the biggest difference is the switch in coding assignments of segments S1 and S4 [2, 10]. In most other orthoreoviruses, the Outer fiber protein is encoded on the same bi- or tricistronic S genome segment as the FAST protein and/or a poorly conserved nonstructural protein of unclear function .
Whole-genome sequence comparison to other members of family Reoviridae
The complete sequencing of Canadian PRV isolates 358, 371, and the Chilean PRV isolate CGA280-5 in the present study enabled us to elaborate the taxonomic grouping of the PRV isolates, and the phylogenetic relationships between Orthoreovirus, Aquareovirus, and PRV at the genome level. Nucleotide sequences of 13 selected members of family Reoviridae, belonging to Orthoreovirus, Aquareovirus, PRV, and the Bluetongue virus as the outgroup sequence used in the analyses are shown in (Additional file 3: Table S3). Because the equivalent PRV segment S4 gene (Outer fiber protein) is not present in Aquareovirus genus, we restricted our analysis to segments homologous to PRV L1, L2, L3, M1, M2, M3, S1, S2, and S3, using a segment to segment comparison approach as done by Palacios et al.. For the 13 isolates (Additional file 3: Table S3), we generated a phylogenetic tree for each of the 9 segments L1, L2, L3, M1, M2, M3, S1, S2, and S3. These trees are shown in Figures 2, 3, 4, 5, 6, 7, 8, 9 and 10. All the 9 trees show that PRV isolates cluster in a separate group. In 3 of the 9 trees (i.e., segments homologous to PRV L3, M1, and M2), isolates belonging to Orthoreovirus, Aquareovirus, and PRV grouped in clearly separate clusters, and would therefore clearly delineate PRV as a new genus separate from both Orthoreovirus and Aquareovirus. In 3 out of the 9 trees (i.e., segments homologous to PRV L1, M3, and S2), all the isolates inside the Orthoreovirus genus are in a separate group. In one tree (i.e., segment homologous to PRV L2), only the isolates inside the Aquareovirus genus are in a separate group. In 2 out of the 9 trees (i.e., segments homologous to PRV S1 and S3), none of the three groups of isolates (Orthoreovirus, Aquareovirus, and PRV) clearly clustered in a separate group although the PRV isolates formed a tight cluster. Thus while the distinction between Orthoreovirus genus and Aquareovirus genus is not consistent, the distinction between PRV isolates and these two genera is very consistent. These observations further argue for assigning PRV to a new genus within the family Reoviridae.
While use of concatenated sequences for phylogenetic analyses  for a segmented dsRNA virus may be of limited value because of issues of reassortment and how these may influence phylogenetic groupings, it was also attempted in this study as such a phylogenetic tree is more objective and comprehensively reflects the relationships among the different isolates in Reoviridae. Because the lengths of these concatemers were in the range of 18kbp to 23kbp, it was difficult to make sure they aligned well across the whole length. To improve the alignment quality, we first created an alignment for every homologous segment of the 13 isolates (Additional file 3: Table S3) separately, and then merged all the aligned nine segments together for each of the isolates. In this way, we were able to ensure that all alignments involved sequences of the homologous segments. The alignments were visually checked and a phylogenetic tree generated based on this concatemer alignment is shown in (Additional file 4: Figure S1a). The concatemer alignment was then examined for highly conserved regions, which would strongly support the assumption that the 13 isolates had a common ancestor and justify the current phylogenetic analysis. These highly conserved regions were found in segments homologous to PRV segments L1, L2, and L3. The computer software JALVIEW  was then used to extract these regions, which were then used to generate another phylogenetic tree (Additional file 4: Figure S1b). A comparison of Figures S1a and S1b revealed that while slight differences exist, both trees individually and in combination support three major groups: Aquareovirus genus, Orthoreovirus genus, and PRV isolates, i.e., they also support the classification of PRV as a member of a new genus within the family Reoviridae.
The distances inside the tree in Additional file 4: Figure S1a provide more insight about this proposal. The average distance between the isolates (MRV, BRV, NBV, ARV138) in the genus Orthoreovirus is 0.526. The average distance between an isolate of genus Orthoreovirus and an isolate of genus Aquareovirus is 0.619. The average distance between a PRV isolate and an isolate of genus Orthoreovirus is 0.588, which is much closer to 0.619 than to 0.526, allowing us to unambiguously conclude that PRV represents a new genus within the family Reoviridae.
Phylogenetic analysis and sequence diversity of PRV genomic segment S1
To determine the genetic relationship between the PRV strains from western Canada, Chile, and Norway, we compared the segment S1 sequences using phylogenetic analyses and percent nucleotide similarities. The sequences of seven of the Norwegian isolates were deposited in the GenBank database as S4 sequences  but correspond to S1 sequences of the PRV type strain, Norwegian isolate Reovirus sp. Salmo/GP-2011/NOR , and the PRV isolates in the present study. Figure 11 shows the phylogenetic tree generated with these sequences. It shows 4 Norwegian isolates 5433, 3817, 1921, 9326, are very close to the Canadian strains and another 4 Norwegian isolates 7243, 7030, GP-2010, 8286 are very close to Chilean strains. This indicates one PRV genotype in Norway, Genotype I, with two major sub-genotypes, which we designate Ia and Ib, with Canadian PRV strains in sub-genotype Ia and Chilean PRV strains in sub-genotype Ib, both with strong bootstrap support. The Canadian PRV strains form two subgroups, (167, 196, 358, 209, bootstrapping support 90.5%, and 163, 371, bootstrapping support 61.2%). The two PRV sub-genotypes are separated by a relatively long branch (Figure 11), suggesting that they have been evolving independently in Norway. Interestingly, the phylogenetic trees of the individual genome segment nucleotide sequences of Reoviridae (Figures 2, 4, 5, 6, 7 and 8, and 10) except for segments homologous to PRV L2 and S2 (Figures 3 and 9, respectively) also seem to support the existence of the two sub-genotypes of PRV.
Piscine reovirus segment S1 is bicistronic (Table 3), encoding the Outer clamp protein and a 124-aa protein, designated p13, that induces cytotoxicity . Thus, S1 may be relevant for virulence of PRV. Our sequence analysis of p13 protein showed the Chilean PRV strains had 100% amino acid sequence identity with the Norwegian strain Reovirus sp. Salmo/GP-2010/NOR, whereas the Canadian strains had ≤92.7% amino acid sequence identity with this PRV strain (Additional file 2: Table S2). In the S1 sequence phylogenetic trees (Figures 11 and 12), the Canadian PRV strains are most similar to the Norwegian PRV strains found in Atlantic salmon with HSMI outbreaks from the Lofoten Archipelago of Norway ; in contrast, the Chilean PRV strains are most similar to the strains found in Atlantic salmon farms without HSMI outbreaks near Trondheim, Norway  (Dr. Torstein Tengs, personal communication). These findings suggest the existence of PRV provides the potential for a HSMI outbreak, but other factors (including environment, stress, PRV/host contact types, PRV infection titre) determine whether a HSMI outbreak actually occurs. This requires further investigation.
Table 4 shows the percent sequence identities on segment S1 between the new Canadian and Chilean PRV strains and the Norwegian PRV isolates reported in the GenBank database. This analysis confirms the phylogenetic analysis (Figures 11 and 12) showing that all Canadian isolates belong to Norwegian sub-genotype Ia and Chilean isolates belong to Norwegian sub-genotype Ib. The Canadian and Norwegian PRV isolates of sub-genotype Ia showed nucleotide sequence identities ≥ 98.1% and amino acid sequence identities ≥ 98.2%. As noted in Table 4, the Chilean PRV S1 sequences 8857 and 337 have inserts relative to the other strains, which contributes to lower sequence identities with other strains. If these two Chilean isolates are excluded, then the Norwegian and Chilean PRV strains of sub-genotype Ib showed nucleotide sequence identities ≥ 98.1% and amino acid sequence identities ≥ 99.4%. The nucleotide and amino acid sequence identities between sub-genotypes Ia and I b strains were ≤ 96.9%. and ≤ 97.0%, respectively, (Table 4).
Divergence time estimation between Canadian and Chilean PRV and the Norwegian strains
Our analysis using BEAST simulation  shows the time when Canadian PRV isolates diverged from Norwegian PRV isolates was between 2006 and 2011; the time when Chilean PRV isolates diverged from Norwegian PRV isolates was between 2003 and 2010. These estimations were based on isolates, collection times (Additional file 1: Tables S1a and S1b) and all the information inside the phylogenetic tree (Figure 11). We also used the program BACKTRACK , which reads a phylogenetic tree with evolutionary distances and years of isolation for all the sequences and then generates a time interval for each inner node, to estimate the divergence times. With BACKTRACK, the Canadian isolates diverged from Norwegian isolates between 2007 and 2011; the Chilean isolates diverged from Norwegian isolates between 2004 and 2011. Both algorithms produced a wide range in the estimations. We tried to make the estimations more specific by using the knowledge we have about these sequences. For example, we believe the multiple insertion events of CGA8857 and CGA337 could be caused by some kind of environmental changes and the mutation rates during that period could be significantly higher than normal. Thus, we believe the most likely time when Canadian isolates diverged from Norwegian isolates was between 2006 and 2008, i.e., around 2007 ± 1; the most likely time when Chilean isolates diverged from Norwegian isolates was between 2007 and 2009, around 2008 ± 1. This evolutionary direction was confirmed with several outgroup sequences. The timeline for Canadian PRV is supported by observations that: 1) Heart lesions and HSMI type lesions were reported in British Columbia farmed Atlantic salmon beginning in 2008, in fish that had entered seawater in 2007. The pattern of inflammation in the heart was consistent with systemic immune stimulation; differentials include a bacterial or viral infection , and 2) a survey of juvenile Pink salmon “Oncorhynchus gorbuscha” in the Broughton Archipelago region of western Canada in April and May 2008 (200 samples in 44 pools) found no PRV when tested with a RT-qPCR assay targeting PRV L1 gene in 2010 . It is not known how the virus could have been transmitted from Norway to Canada since there have never been any authorized direct imports of Atlantic salmon eggs from Norway since 1985; recent imports have been from Washington State-USA (2001) and Iceland (2004–2009) . There is no information about the PRV situation in Washington State or Iceland. Horizontal spread and/or introduction of virus through wild fish migration are not reasonable routes of transmission. The distribution of PRV in Canada is uncertain since there is no national surveillance for it. In Chile, the presence of PRV in farmed Atlantic salmon reared in Chilean seawater was first detected in 2010 and published in a laboratory report in 2011 , which cited a high prevalence among the sites located in different areas of the Los Lagos region growing-up macro-zone. PRV could have been introduced to Chile through importation of Atlantic salmon eggs similarly to ISAV [37–39] albeit more recently than ISAV: most Atlantic salmon egg imports to Chile in 2008 were from Norway, and from 2009–2013 have been from Iceland . In Chile, HSMI will be recommended to be included on the List 3 of high-risk diseases , thus ensuring active surveillance for it in Chilean aquaculture. To better understand the molecular epidemiology of PRV, it will be necessary to know the situation of PRV in other salmon producing countries and in countries with wild salmonids. Not all countries have surveillance for HSMI; therefore, it may be advisable to consider PRV-HSMI as an emerging disease and initiate its surveillance.
The present work constitutes the first report of genomic analysis of PRV strains detected in tissue samples obtained from fish in western Canada, and in Chile, extending the geographical range of the characterized virus to Pacific shorelines of both North and South America. Our work suggests PRV entered both Chile and western Canada recently. We provide strong support for classification of PRV as a member of a new genus within the family Reoviridae. Our work groups PRV into one genotype, Genotype I, with two major sub-genotypes designated Ia and Ib, with Canadian PRV strains in sub-genotype Ia and Chilean PRV strains in sub-genotype Ib. Taken together, these findings raise awareness on PRV existence outside of Norway so that the aquaculture industry and wild fisheries managers worldwide can become proactive and curtail its international spread, as well as implement mitigation measures regionally and locally.
Fish samples and processing
All samples used in this study were submitted to the laboratory for testing for PRV and other aquatic animal viruses. Samples were either taken from fresh fish and put into microcentrifuge tubes containing RNAlater® (Ambion Inc., Foster City CA) preservative or were immediately placed in sterile whirlpak bags and shipped overnight by courier cold on ice to the testing laboratory; samples, which consisted of individual gill, heart, kidney or pooled gill and heart or kidney, heart and liver tissues, were either bagged individually or pooled (2–3 tissues per pool) for each fish. The fish tissue sample source is detailed in Additional file 1: Table S1. The samples from western Canada were either harvest samples of marine-farmed Atlantic salmon or wild fish samples from fish caught live under sport, scientific or First Nation licenses in regions where Atlantic salmon aquaculture sites were present (Additional file 1: Table S1a).
Each tissue (or pool of tissues) was weighed and macerated to a 10% suspension w/v in phosphate buffered saline (PBS) with 10x antibiotics. The specimen supernatant was used for RNA extraction. The samples from Chile were from marine-farmed Atlantic salmon after seawater transfer (Additional file 1: Table S1b) and were collected in RNAlater® preservative. Samples preserved in RNAlater® (Ambion Inc) were first washed three times with PBS and then homogenized as described above prior to total RNA extraction.
RT-PCR and nucleic acid sequencing
Total RNA was isolated using a modified total RNA extraction protocol with the RNeasy® mini Kit (QIAGEN). Briefly, total RNA was isolated from samples using 1.25 ml of TRIZOL Reagent (Invitrogen) and 375 μl of sample volume as previously described . The Viral RNA mini Kit (QIAGEN) was also utilized on selected samples following the manufacturer’s recommended protocol. In all cases, the extracted RNA was eluted in 20–50 μl of nuclease-free water, and RNA yields were quantified and purity analysed using the OD260/280 ratio and a NanoDrop ND-1000 spectrophotometer (Thermo Scientific). The eluted RNA was tested immediately following quantitation, or was stored frozen at −80°C until use.
RT-qPCR was run on the LightCycler 480 (Roche Applied Science), version 4.0. The crossing point (Cp) or threshold cycle (Ct) was determined by use of the maximum-second-derivative function on the LightCycler software release 1.5.0. The OneStep RT-PCR kit (QIAGEN) was employed for all RT-qPCR reactions according to the manufacturer’s specifications.
Sample RNA quality was based on RT-qPCR for elongation factor 1 alpha (ELF-1α) as internal control targeting either Atlantic salmon ELF-1α (GenBank accession number BT072490) or Chinook salmon ELF-1α (GenBank accession number FJ890356) carried out using Roche LightCycler® 480 RNA master Hydrolysis Probe kit (Roche Diagnostics). The following primers and probes were used.
For Atlantic salmon EF1α:
ASELF1a Forward – 5'- CGT GAC ATG AGG CAG ACA GT-3';
ASELF1a Reverse – 5'- CGG CCT TAA CAG CAG ACT TTG-3';
ASELF1a Probe – 5'-TGC TGT CGG TGT CAT CAA GGC T-3'; and
For Chinook salmon ELF1α:
CSELF1α Forward – 5′- GGT CAC CAC CTA CAT CAA GAA GA-3′;
CSELF1α Reverse – 5′- CCA ACC AGA GAT GGG CAC AAA G-3′;
CSELF1α Probe – 5′-TGG CTA CAA CCC TGC CAC TGT C-3′.
The probes were labelled at the 5′ end with 6-FAM and at the 3′ end with BHQ-1 quencher (Biosearch Technologies Inc.). The final concentrations of primers and probe in each case were 900 nM for each primer and 250 nM for the probe in a final volume of 25 μl. The following thermal cycling parameters were used: 1 cycle of RT for 3 min at 63°C followed by denaturation at 95°C for 3 s, and 45 cycles of denaturation at 95°C for 15 s, annealing and detection at 60°C for 1 min and extension at 72°C for 1 s. Ct values above 40 and no Ct values were defined as negative and these samples would be considered unfit for further testing if after re-extraction and repeated RT-qPCR yielded the same results.
The RT-qPCR assay for PRV used the primer-probe set sequences developed by Haugland et al.  targeting the PRV L1 gene. The primers were PRV-F– 5′-CCC CAT CCC TCA CAT ATG GAT A-3′and PRV-R– 5′-GGT GAA ATC ATC GCC AAC TCA-3′. The PRV probe, which was labelled at the 5′ end with 6-FAM and at the 3′ end with BHQ-1 quencher (Biosearch Technologies Inc.) was 5′-ATG TCC AGG TAT TTA CC-3′. The reaction conditions were the same as used by Palacios et al., but with 8 μl of template RNA. The following concentrations were used: 400 nM primer, 300 nM probe and 1.25 mM MgCl2. The following thermal cycling parameters were used: 1 cycle of RT for 30 min at 50°C followed by denaturation at 94°C for 15 min, and 45 cycles of denaturation at 94°C for 15 s, annealing at 54°C for 30 s and amplification and detection at 72°C for 15 s. Samples to be considered positive had Ct values up to 40 and with an exponential curve; Ct values between 40.1 and 45 were considered suspicious, and a sample was negative if there was no Ct value.
Because the laboratory did not have a PRV isolate from cell culture to use as positive control, samples with positive Ct values were further tested in classic RT-PCR targeting the 3′ portion of genome segment L1 with the following PCR primer pairs: PRV-L1 For1 – 5′-CAC TCA CCA ATG ACC CAA ATG C-3′; PRV-L1 Rev1 – 5′-TTG ACA GTC TGG CTA CTT CGG-3′ and/or PRV-L1 For2 – 5′-CTG AAC TGC TAG TTG AGG ATG G-3′; PRV-L1 Rev2 – 5′-GCC AAT CCA AAC AGA TTA GG-3′. These PCR primers and those used to amplify the 10 genomic segments of PRV, listed in (Additional file 5: Table S4) were designed based on the published PRV sequences . RT-PCR for the amplification of each viral genome segment was carried out by using the OneStep RT-PCR kit (QIAGEN). Briefly, the reaction mixture contained 1 μl of total RNA, 4 μl of 5X QIAGEN OneStep RT-PCR buffer, 0.8 μl of dNTPs, 0.5 μM (final concentration) of each primer pair, and 0.8 μl of QIAGEN OneStep RT-PCR enzyme mix in a final volume of 20 μl. Thermal cycling conditions were as follows: an initial cycle of 50°C for 40 min and 95°C for 10 min; 40 cycles of 95°C for 30 s, 54°C for 30 sec, 72°C for 70 s; and a final extension cycle of 72°C for 10 min. Amplified products were analyzed by electrophoresis on 1% agarose gel and purified using High Pure PCR Product Purification Kit (Roche). The PCR products were then either directly sequenced or they were cloned into the pCRII vector using a TOPO TA cloning kit (Invitrogen) in preparation for nucleotide sequencing. Plasmid DNA for sequencing was prepared as described before . DNA sequencing was performed as previously described  by ACGT Corporation (Toronto, Ontario, Canada). DNA Sequencing was done either directly on RT-PCR products or on plasmid DNA containing the cloned RT-PCR products obtained from reactions using total RNA from tissue samples.
Similarity analysis was performed using BLAST programs available via the National Center for Biotechnology Information  and the FASTA program package for personal computers . Analysis to identify putative ORFs and their predicted amino acid sequences and other protein characteristics was conducted using the Sequence Manipulation suite, version 2 .
The Canadian and Chilean PRV sequences used in the phylogenetic analyses are described in (Additional file 2: Table S2). All the Norwegian PRV sequences were obtained from GenBank . Sequences were processed using ClustalX 2.1 . The multiple sequence alignment was manually verified and adjusted to reach high quality alignment. The phylogenetic trees were generated when positions with gaps were excluded and corrections for multiple substitutions were used. Bootstrapping was performed for 1,000 times. In most cases, only the bootstrapping supports higher than 70% were noted. For some important branches, those bootstrapping values a little lower than 70% were also noted. To verify the evolution direction, outgroup sequences were used to determine the root of the phylogenetic trees.
Divergence time estimation in a rooted phylogenetic tree
BEAST v1.7.5  was used to estimate divergence time. To find the most suitable substitution model, we ran jModelTest 0.1.1  against the aligned sequences. The result shows K80 model  is the most suitable model. Based on this result, a similar model, HKY85 model , was chosen in the BEAST simulation. We believe the mutation rates among lineages could be different, and the uncorrelated relaxed molecular clock was chosen. Five million simulation steps were performed and enough effective sample sizes (ESSs) were generated. We also used program BACKTRACK , which reads a phylogenetic tree with evolutionary distances and years of isolation for all the sequences and then generates a time interval for each inner node, to estimate the divergence times.
Attoui H, Mertens PPC, Becnel J, Belaganahalli S, Bergoin M, Brussaard CP, Chappell JD, Ciarlet M, del Vas M, et al.: Family Reoviridae. In Virus Taxonomy: Ninth Report of the International Committee on Taxonomy of Viruses. Edited by: King AMQ, Adams MJ, Carstens EB, Lefkowitz EJ. San Diego, CA: Academic Press; 2011:541-637.
Palacios G, Lovoll M, Tengs T, Hornig M, Hutchison S, et al.: Heart and skeletal muscle inflammation of farmed salmon is associated with infection with a novel reovirus. PLoS One 2010, 5: e11487. 10.1371/journal.pone.0011487
Duncan R: Extensive sequence divergence and phylogenetic relationships between the fusogenic and nonfusogenic orthoreoviruses: a species proposal. Virology 1999, 260: 316. 10.1006/viro.1999.9832
Shmulevitz M, Duncan R: A new class of fusion-associated small trans-membrane (FAST) proteins encoded by the non-enveloped fusogenic reoviruses. EMBO J 2000, 19: 902-912. 10.1093/emboj/19.5.902
Sabin A: Reoviruses: a new group of respiratory and enteric viruses formerly classified as ECHO type 10 is described. Science 1959, 130: 1387-1389. 10.1126/science.130.3386.1387
Chua K, Crameri G, Hyatt A, Yu M, Tompang M, et al.: A previously unknown reovirus of bat origin is associated with an acute respiratory disease in humans. Proc Nat Acad Sci USA 2007, 104: 11424-11429. 10.1073/pnas.0701372104
Leland MM, Hubbard GB, Sentmore HT, Soike KF, Hilliard JK: Outbreak of orthoreovirus-induced meningoencephalomyelitis in baboons. Comp Med 2000, 50: 199-205.
van den Brand JM, Manvell R, Paul G, Kik MJ, Dorrestein GM: Reovirus infections associated with high mortality in psittaciformes in The Netherlands. Avian Pathol 2007, 36: 293-299. 10.1080/03079450701447309
Lamirande EW, Nichols DK, Owens JW, Gaskin JM, Jacobson ER: Isolation and experimental transmission of a reovirus pathogenic in ratsnakes ( Elaphe species). Virus Res 1999, 63: 135-141. 10.1016/S0168-1702(99)00067-2
Key T, Read J, Nibert ML, Duncan R: Piscine reovirus encodes a cytotoxic, non-fusogenic, integral membrane protein and previously unrecognized virion outer-capsid proteins. J Gen Virol 2013, 94: 1039-1050. 10.1099/vir.0.048637-0
Kibenge FSB, Godoy MG, Fast M, Workenhe S, Kibenge MJT: Countermeasures against viral diseases of farmed fish. Antiviral Res 2012, 95: 257-281. 10.1016/j.antiviral.2012.06.003
Kongtorp RT, Kjerstad A, Taksdal T, Guttvik A, Falk K: Heart and skeletal muscle inflammation in Atlantic salmon Salmo salar L: a new infectious disease. J Fish Dis 2004, 27: 351-358. 10.1111/j.1365-2761.2004.00549.x
Ferguson HW, Kongtorp RT, Taksdal T, Graham D, Falk K: An outbreak of disease resembling heart and skeletal muscle inflammation (HSMI) in Scottish farmed salmon ( Salmo salar L) with observations on myocardial regeneration. J Fish Dis 2005, 28: 119-123. 10.1111/j.1365-2761.2004.00602.x
Wiik-Nielsen CR, Løvoll M, Sandlund N, Faller R, Wiik-Nielsen J, Jensen BB: First detection of piscine reovirus (PRV) in marine fish species. Dis Aquat Org 2012, 97: 255-258. 10.3354/dao02425
Garseth ÅH, Fritsvold C, Opheim M, Skjerve E, Biering E: Piscine reovirus (PRV) in wild Atlantic salmon, Salmo salar L., and sea-trout, Salmo trutta L., in Norway. J Fish Dis 2012. 10.1111/j.1365-2761.2012.01450.x
National Veterinary Institute: Farmed Fish Health Report 2007. Edited by: Skjelstad H, Bornø G, Flesjå K, Hansen H, Nilsen H, Wasmuth MA, Hjeltnes B. Oslo: National Veterinary Institute; 2007:19.
Bornø G, Sviland C: The health situation in farmed salmonids 2010. In Farmed Fish Health Report 2010. Edited by: Bornø G, Sviland C, Hellberg H. Oslo: National Veterinary Institute; 2010:36.
Løvoll M, Alarcón M, Jensen BB, Taksdal T, Kristoffersen AB, Tengs T: Quantification of piscine reovirus (PRV) at different stages of Atlantic salmon Salmo salar production. Dis Aquat Org 2012, 99: 7-12. 10.3354/dao02451
Mikalsen AB, Haugland O, Rode M, Solbakk IT, Evensen O: Atlantic salmon reovirus infection causes a CD8 T Cell myocarditis in Atlantic salmon (Salmo salar L.). PLoS One 2012, 7: e37269. 10.1371/journal.pone.0037269
Chi SC, Hu WW, Lo BJ: Establishment and characterization of a continuous cell line (GF-1) derived from grouper, Epinephelus coioides (Hamilton): a cell line susceptible to grouper nervous necrosis virus (GNNV). J Fish Dis 1999, 22: 173-182. 10.1046/j.1365-2761.1999.00152.x
Kemp M, Donovan J, Higham H, Hooper J: Biochemical markers of myocardial injury. Br J Anaesth 2004, 93: 63-73. 10.1093/bja/aeh148
Yousaf MN, Powell MD: The effects of heart and skeletal muscle inflammation and cardiomyopathy syndrome on creatine kinase and lactate dehydrogenase levels in Atlantic salmon ( Salmo salar L.). Sci World J 2012, 741302: 9.
Løvoll M, Wiik-Nielsen J, Grove S, Wiik-Nielsen CR, Kristoffersen AB, Faller F, et al.: A novel totivirus and piscine reovirus (PRV) in Atlantic salmon ( Salmo salar ) with cardiomyopathy syndrome (CMS). Virol J 2010, 7: 309. 10.1186/1743-422X-7-309
Wiik-Nielsen J, Løvoll M, Fritsvold C, Kristoffersen AB, Haugland Ø, Hordvik I, Aamelfot M, Jirillo E, Koppang EO, Grove S: Characterization of myocardial lesions associated with cardiomyopathy syndrome in Atlantic salmon, Salmo salar L., using laser capture microdissection. J Fish Dis 2012, 35: 907-916. 10.1111/j.1365-2761.2012.01431.x
Anonymous 2012: Salmon mortality levels rise in H1. (Published August 17, 2012; Consulted August 20, 2012) http://fis.com/fis/worldnews/worldnews.asp?l=e&country=0&special=aquaculture&monthyear=&day=&id=54692&ndb=1&df=0
Bustos P, Rozas M, Bohle H, Iidefonso R, Sandoval A, Gaete A, Araya C, Grothusen H, Tapia E, Gallardo A, Rojas M, Lara M, Labra A, Gálvez C: Primer reporte de piscine reovirus en Salmon del Atlántico, Salmo salar, cultivado en Chile. ADL Diagnostic Chile Ltda. 2011, 7: 1-4.
Feinberg J: Science divided on fish virus found in Cultus Lake trout. http://www.theprogress.com/news/163480656.html
Løvoll M, Alarcón M, Jensen BB, Taksdal T, Kristoffersen AB, Tengs T: Quantification of piscine reovirus (PRV) at different stages of Atlantic salmon Salmo salar production. Dis Aquat Org 2012, 99: 7-12. supplementary material ) http://www.int-res.com/articles/suppl/d099p007_supp.pdf supplementary material ) 10.3354/dao02451
Benson DA, Cavanaugh M, Clark K, Karsch-Mizrachi I, Lipman DJ, Ostell J, Sayers EW: GenBank. Nucleic Acids Res 2013,41(D1):D36-D42. 10.1093/nar/gks1195
Yan X, Parent KN, Goodman RP, Tang J, Shou J, Nibert ML, Duncan R, Baker TS: Virion structure of baboon reovirus, a fusogenic orthoreovirus that lacks an adhesion fiber. J Virol 2011, 85: 7483-7495. 10.1128/JVI.00729-11
Waterhouse AM, Procter JB, Martin DMA, Clamp M, Barton GJ: Jalview Version 2-a multiple sequence alignment editor and analysis workbench. Bioinformatics 2009, 25: 1189-1191. 10.1093/bioinformatics/btp033
Drummond AJ, Suchard MA, Xie D, Rambaut A: Bayesian phylogenetics with BEAUti and the BEAST 1.7. Mol Biol Evol 2012, 29: 1969-1973. 10.1093/molbev/mss075
Kibenge FSB, Kibenge MJ, Wang Y, Qian B, Hariharan S, McGeachy S: Mapping of putative virulence motifs on infectious salmon anaemia virus surface glycoprotein genes. J Gen Virol 2007, 88: 3100-3111. 10.1099/vir.0.83097-0
The Honourable Bruce I: Cohen, Commissioner. Commission of Inquiry into the Decline of the Sockeye Salmon in the Fraser River, Final Report, October 2012, Exhibit #1549–310. Available at http://www.cohencommission.ca/en/Exhibits.php
Saksida SM, Marty GD, Jones SRM, Manchester HA, Diamond CL, Bidulka J, St-Hilaire S: Parasites and hepatic lesions among pink salmon, Oncorhynchus gorbuscha (Walbaum), during early seawater residence. J Fish Dis 2012, 35: 137-151. 10.1111/j.1365-2761.2011.01330.x
DFO: Public reporting on aquaculture in the Pacific Region – salmon egg imports. 2012. (Published 2012-03-07; Consulted 2012-12-27) http://www.pac.dfo-mpo.gc.ca/aquaculture/reporting-rapports/egg-oeuf-eng.htm
Vike S, Nylund S, Nylund A: ISA virus in Chile: evidence of vertical transmission. Arch Virol 2009, 154: 1-8. 10.1007/s00705-008-0251-2
Kibenge FSB, Godoy MG, Wang Y, Kibenge MJT, Gherardelli V, Mansilla S, Lisperger A, Jarpa M, Larroquete G, Avendaño F, Lara M, Gallardo A: Infectious salmon anaemia virus (ISAV) isolated from the ISA disease outbreaks in Chile diverged from ISAV isolates from Norway around 1996 and was disseminated around 2005, based on surface glycoprotein gene sequences. Virol J 2009, 6: 88. 10.1186/1743-422X-6-88
Plarre H, Nylund A, Karlsen M, Brevik Ø, Sæther PA: Evolution of infectious salmon anaemia virus (ISA virus). Arch Virol 2012, 157: 2309-2326. 10.1007/s00705-012-1438-0
Sernapesca: Estadística de importación de ovas por origen 2008 – 2013. http://www.sernapesca.cl/index.php?option=com_remository&Itemid=246&func=fileinfo&id=4696
Workenhe ST, Kibenge MJT, Iwamoto T, Kibenge FSB: Absolute quantitation of infectious salmon anaemia virus using different real-time reverse transcription PCR chemistries. J Virol Methods 2008, 154: 128-134. 10.1016/j.jviromet.2008.08.007
Haugland Ø, Mikalsen AB, Nilsen P, Lindmo K, Thu BJ, Eliassen TM, Roos N, Rode M, Evensen O: Cardiomyopathy syndrome of Atlantic salmon (Salmo salar L.) is caused by a double-stranded RNA virus of the Totiviridae family. J Virol 2011, 85: 5275-5286. 10.1128/JVI.02154-10
Kibenge FSB, Dybing JK, McKenna PK: Rapid procedure for large-scale isolation of plasmid DNA. Biotechniques 1991, 11: 65-67.
Kibenge FSB, Xu H, Kibenge MJT, Qian B, Joseph T: Characterization of gene expression on genomic segment 7 of infectious salmon anaemia virus. Virol J 2007, 4: 34. 10.1186/1743-422X-4-34
Altschul SF, Gish W, Miller W, Myers EW, Lipman DJ: Basic local alignment search tool. J Mol Biol 1990, 215: 403-410.
Pearson WR, Lipman DJ: Improved tools for biological sequence comparison. Proc Nat Acad Sci USA 1988, 85: 2444-2448. 10.1073/pnas.85.8.2444
Larkin MA, Blackshields G, Brown NP, Chenna R, McGettigan PA, McWilliam H, Valentin F, Wallace IM, Wilm A, Lopez R, Thompson JD, Gibson TJ, Higgins DG: Clustal W and Clustal X version 2.0. Bioinformatics 2007, 23: 2947-2948. 10.1093/bioinformatics/btm404
Posada D: jModelTest: Phylogenetic Model Averaging. Mol Biol Evol 2008, 25: 1253-1256. 10.1093/molbev/msn083
Kimura M: A simple method for estimating evolutionary rates of base substitutions through comparative studies of nucleotide sequences. J Mol Evol 1980, 16: 111-120. 10.1007/BF01731581
Hasegawa M, Kishino H, Yano T: Dating of human-ape splitting by a molecular clock of mitochondrial DNA. J Mol Evol 1985, 22: 160-174. 10.1007/BF02101694
Cheng L, Fang Q, Shah S, Atanasov IC, Zhou ZH: Subnanometer-resolution structures of the grass carp reovirus core and virion. J Mol Biol 2008, 382: 213-222. 10.1016/j.jmb.2008.06.075
This work was supported by the Virology Research Laboratory at the Atlantic Veterinary College, University of Prince Edward Island, Charlottetown, PE, Canada. Van City supported the work by Alexandra Morton. We thank Dr. Richard Routledge for submitting the Cutthroat trout samples for diagnostic testing.
The authors declare that they have no competing interests in this scientific work.
MJTK isolated total RNA from tissue samples, performed the RT-qPCR for ELF-1α and PRV and developed and performed the classic RT-PCR for 3′ portion of segment L1 of PRV, and helped to write the manuscript. TI designed the PRV PCR primers used to amplify transcripts of all the 10 genome segments, performed the classic RT-PCR and cloned all PCR products for sequencing and helped to write the manuscript. YW performed all the phylogenetic analyses and helped to write the manuscript. AM provided the Canadian samples for diagnostic testing and edited the manuscript. MGG coordinated the laboratory testing in Chile and edited the manuscript. FSBK coordinated all viral testing and DNA sequence analysis and helped to write the manuscript. All authors read and approved the final manuscript.
Molly JT Kibenge, Tokinori Iwamoto contributed equally to this work.
Electronic supplementary material
Additional file 1:Title: Additional results on the piscine reovirus (PRV) positive samples from Canada and Chile. Description: Two tables showing RT-qPCR and conventional RT-PCR results of fish tissue samples from Canada and Chile tested for PRV. (DOC 64 KB)
Additional file 2:List of new piscine reovirus (PRV) nucleotide sequences and their percent identity to Norwegian isolate Salmo/GP-2010/NOR .(DOC 94 KB)
Additional file 3:GenBank Accession numbers of genome segments of selected members of family Reoviridae used in phylogenetic comparison of nucleotide sequences of individual genome segments [].(DOC 40 KB)
Additional file 4:Title: Phylogenetic trees showing the relationships between isolates in family Reoviridae at the genome-level. Description: (Figure S1a) Concatenated sequences of nine homologous segments (segment L1, L2, L3, M1, M2, M3, S1, S2, S3) shared by piscine reovirus (PRV) and selected members of family Reoviridae, were used to generate a phylogenetic tree. (Figure S1b) Phylogeny of highly-conserved regions of concatemers in Figure S1a. (PDF 55 KB)
Additional file 5:Title: List of oligonucleotide primers used in amplification of piscine reovirus (PRV) genome segments. Description: Table listing oligonucleotide primers. (DOC 58 KB)
Authors’ original submitted files for images
Below are the links to the authors’ original submitted files for images.
About this article
Cite this article
Kibenge, M.J., Iwamoto, T., Wang, Y. et al. Whole-genome analysis of piscine reovirus (PRV) shows PRV represents a new genus in family Reoviridae and its genome segment S1 sequences group it into two separate sub-genotypes. Virol J 10, 230 (2013). https://doi.org/10.1186/1743-422X-10-230