Next Article in Journal
Triatomine Feeding Profiles and Trypanosoma cruzi Infection, Implications in Domestic and Sylvatic Transmission Cycles in Ecuador
Next Article in Special Issue
Challenges and Solutions to Viral Diseases of Finfish in Marine Aquaculture
Previous Article in Journal
A 2-Year Longitudinal Seroepidemiological Evaluation of Toxoplasma gondii Antibodies in a Cohort of Autochthonous Sheep from Central Portugal
Previous Article in Special Issue
Low Pathogenic Strain of Infectious Pancreatic Necrosis Virus (IPNV) Associated with Recent Outbreaks in Iranian Trout Farms
 
 
Search for Articles:
Title / Keyword
Author / Affiliation / Email
Journal
Article Type
 
 
Advanced Search
 
Section
Special Issue
Volume
Issue
Number
Page
 
Journals
Pathogens
Volume 10
Issue 1
10.3390/pathogens10010041
pathogens-logo
Submit to this Journal Review for this Journal Propose a Special Issue
► Article Menu

Article Menu

share Share announcement Help format_quote Cite question_answer Discuss in SciProfiles thumb_up
...
Endorse textsms
...
Comment
first_page
settings
Order Article Reprints
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Communication

Extensive Phylogenetic Analysis of Piscine Orthoreovirus Genomic Sequences Shows the Robustness of Subgenotype Classification

by
Marcos Godoy
1,2,3,
Daniel A. Medina
1,
Rudy Suarez
2,4,
Sandro Valenzuela
2,
Jaime Romero
5,
Molly Kibenge
6,
Yingwei Wang
7 and
Frederick Kibenge
6,*
1
Laboratorio de Biotecnología Aplicadas, Facultad de Medicina Veterinaria, Universidad San Sebastián, Sede De la Patagonia, Lago Panguipulli 1390, Puerto Montt 5480000, Chile
2
Centro de Investigaciones Biológicas Aplicadas (CIBA), Lago Panguipulli 1390, Puerto Montt 5480000, Chile
3
Doctorado en Acuicultura, Programa Cooperativo Universidad de Chile, Universidad CatÓlica del Norte, Coquimbo 1781421, Chile
4
Universidad CatÓlica del Norte, Coquimbo 1781421, Chile
5
Laboratorio de Biotecnología, Instituto de Nutrición y Tecnología de los Alimentos, Universidad de Chile, Santiago 7800024, Chile
6
Department of Pathology and Microbiology, Atlantic Veterinary College, University of Prince Edward Island, Charlottetown, PE C1A4P3, Canada
7
School of Mathematical and Computational Sciences, University of Prince Edward Island, Charlottetown, PE C1A4P3, Canada
*
Author to whom correspondence should be addressed.
Pathogens 2021, 10(1), 41; https://doi.org/10.3390/pathogens10010041
Submission received: 13 December 2020 / Revised: 4 January 2021 / Accepted: 5 January 2021 / Published: 7 January 2021
(This article belongs to the Special Issue Virulence Mechanisms, Detection and Control of Aquatic Animal Viruses)
Browse Figures
Versions Notes

Abstract

:
Piscine orthoreovirus (PRV) belongs to the family Reoviridae and has been described mainly in association with salmonid infections. The genome of PRV consists of about 23,600 bp, with 10 segments of double-stranded RNA, classified as small (S1 to S4), medium (M1, M2 and M3) and large (L1, L2 and L3); these range approximately from 1000 bp (segment S4) to 4000 bp (segment L1). How the genetic variation among PRV strains affects the virulence for salmonids is still poorly understood. The aim of this study was to describe the molecular phylogeny of PRV based on an extensive sequence analysis of the S1 and M2 segments of PRV available in the GenBank database to date (May 2020). The analysis was extended to include new PRV sequences for S1 and M2 segments. In addition, subgenotype classifications were assigned to previously published unclassified sequences. It was concluded that the phylogenetic trees are consistent with the original classification using the PRV genomic segment S1, which differentiates PRV into two major genotypes, I and II, and each of these into two subgenotypes, designated as Ia and Ib, and IIa and IIb, respectively. Moreover, some clusters of country- and host-specific PRV subgenotypes were observed in the subset of sequences used. This work strengthens the subgenotype classification of PRV based on the S1 segment and can be used to enhance research on the virulence of PRV.

1. Introduction

Heart and skeletal muscle inflammation (HSMI) is an infectious disease caused by Piscine orthoreovirus (PRV) and was described for the first time in Norway [1,2,3]. PRV belongs to the genus Orthoreovirus [4], subfamily Spinareovirinae in the Reoviridae family [1].
PRV has been detected in numerous wild and farmed salmonid species in Norway [5,6,7,8,9], Canada [10,11,12,13,14], Chile [10,15], the United States [12,16], Japan [17], Scotland, Ireland, the Faroe Islands, Iceland, Germany, Sweden, Denmark, Italy, France [18,19,20,21,22,23] and the Czech Republic [24].
In Atlantic salmon (Salmo salar) and Coho salmon (Oncorhynchus kisutch) in Chile the clinical signs of affected fish are mainly characterized by systemic circulatory disturbance. Histopathologically, the fish present myocarditis and myositis of the red musculature [2,15,25,26]. Moreover, during PRV infection, the presence of intracytoplasmic inclusion bodies in erythrocytes has been described [26,27,28].
The genome of PRV consists of about 23,600 bp, which encode at least 11 proteins distributed across ten segments of double-stranded RNA, classified as small (S1, S2, S3 and S4), medium (M1, M2 and M3) and large (L1, L2 and L3); these range approximately from 1000 bp (segment S4) to 4000 bp (segment L1) [1,4,10]. Phylogenetic analysis, based mainly on the PRV S1 segment, differentiates the virus into two genotypes, I and II, and each of these into two major subgenotypes, designated as Ia and Ib, and IIa and IIb, respectively [10,15]. Alternatively, Garseth et al. [29] described four genogroups, I, II, III and IV. A complementary classification using the coding sequence of PRV segments described three subtypes of PRV, designated as PRV-1, PRV-2 and PRV-3 [17,19], with subgenotypes Ia and Ib making up the PRV-1 subtype and subgenotypes IIb and IIa corresponding to the PRV-2 and PRV-3 subtypes, respectively [30].
PRV subgenotypes Ia and Ib have been associated with HSMI in Atlantic salmon and HSMI-like disease in Rainbow trout (Oncorhynchus mykiss) and Coho salmon [1,3,7,13,15]. In Canada, the infection of Chinook salmon (O. tshawytscha) with PRV subgenotype Ia is associated with Jaundice disease [14,31] and in Japan, infection with the PRV-2 (subgenotype IIb) causes erythrocytic inclusion body syndrome (EIBS) [17]. No PRV subgenotype or subtype has been described to be host-specific or geographically exclusive, apart from PRV-Ib, which has not been found in Canada or the USA to date [10,32], and PRV-2, which to date has been reported only in Japan [17]. Consequently, the wide range of hosts, clinical signs and geographic distribution have turned PRV into an emerging virus in the salmon aquaculture industry worldwide.
The aim of this study was to revisit the molecular classification, based on phylogenetic analysis, of almost all PRV sequences available to date (May 2020) in the GenBank database, with a focus on the S1 and M2 genome segments, including unpublished new sequences generated for this work. Our results strengthen the subgenotype classification of PRV based on the S1 segment and allow the classification of previously unclassified, publicly available sequences into the PRV subgenotypes through phylogenetic dendrograms.

2. Results

Phylogenetic Analyses of PRV Genomic Sequences

We retrieved all the available sequences for the ten genomic segments of PRV listed under the Tax ID: txid1157337 in the National Center for Biotechnology Information (NCBI) GenBank database (Supplementary Table S1). New sequences for the S1 and M2 segments were added in preparation for this work, resulting in the obtainment of between 38 and 390 sequences with respect to the PRV segments used in the study (Table 1).
The sequences collected were used to perform an initial phylogenic analysis based on subgenotype classifications, obtaining a clear differentiation using PRV S1 and M2 segments. In accordance with previous observations, it was concluded that the use of these segments would enable the grouping of the PRV into subgenotypes [10,15]. We performed a restrictive analysis of the S1 and M2 segments, discarding sequences with less than 400 nucleotides for the S1 segment and less than 1000 nucleotides for the M2 segment, to cover the conserved coding sequences of σ3 and μ1 proteins. A total of 356 and 59 sequences, respectively, were retained in the end, including the 11 new sequences each for S1 and M2 segments reported in this work (Table 2). The circular dendrogram of the PRV S1 sequences is presented in Figure 1, showing the assignment of new and published unclassified sequences according to their positions in the dendrogram branches. A rectangular version of Figure 1 is provided to enable better exploration of the data (Supplementary Figure S1). The phylogenetic tree of PRV M2 sequences is presented in Figure 2, showing the assignment of new and published unclassified sequences according to their positions in the dendrogram branches. These phylogenetic analyses were able to group the PRV genogroups and allowed us to classify both new sequences (Table 2) and published unclassified sequences into the PRV subgenotypes (Supplementary Table S2). The remaining eight genomic segments had lower resolution and representativity and did not allow for a clear ordination using either subgenotype [10] or subtype [20] classifications (Supplementary Figure S2).
We then used 14 and 17 genomic sequences belonging to M2 and S1 segments, respectively, to graphically show the distances between the sequences according to their classification by subgenotype. Both genomic distances matrixes obtained show good fitness in differentiating between the PRV subgenotypes (Figure 3 and Figure 4). Then, we extended this analysis using 40 and 61 amino acidic sequences belonging to M2 and S1 segments that encode the proteins μ1 and σ3, respectively, including the new sequences produced in this work, classified by subgenotype, subtype, country and host. The distance matrix built using the amino acid sequences of the protein μ1 (M2 segment), presented in Supplementary Figure S3, shows less robustness in the heatmaps used to cluster PRV by subgenotype or subtype. In contrast, the protein σ3 heatmaps show better clustering using subgenotype and subtype classifications (Supplementary Figure S4). Thus, two clear clusters belonging to subgenotypes Ia and Ib were observed in the PRV-1 subtype; these served to strengthen this subtype classification (Figure 4 and Supplementary Figure S3). Country (Supplementary Figure S5) and host species (Supplementary Figure S6) clustering were observed when heatmaps were drawn using these metadata variables. We observed, mainly with the protein σ3 (S1 segment), that the Chile, Norway and Canada sequences shared low distance, suggesting high similarity. However, in the Chile sequences there were two marked groups, one sharing low distance sequences with Canada and Norway and the other clustering with a few Italy and Denmark sequences (Supplementary Figure S5). This agrees with previous reports of the presence of two subgenotypes in Chile [10,15]. Regarding the host species, the distance matrixes inside the Coho salmon sequences were found to have two clear clusters, thereby indicating that this host may carry two PRV subgenotypes. In addition, Coho salmon was found to share low distance clusters with Rainbow trout and a few sequences belonging to Atlantic salmon, suggesting that they harbor similar PRV subgenotypes. We found similar observations in the case of Atlantic salmon. First, we observed two clusters of low distance, suggesting that Atlantic salmon harbors two kinds of PRV sequences. Moreover, Atlantic salmon shared low distance clustering with Chinook salmon and a few sequences of Rainbow trout and Coho salmon (Supplementary Figure S6).

3. Discussion

We extended the geographical range of the already characterized PRV isolates by compiling all the available sequences to date (May 2020) associated with the NCBI Taxonomy ID: 1157337 of PRV. Moreover, we added new sequences for the S1 and M2 segments of 11 PRV isolates from Chile, where salmonids are non-native species and where Coho salmon, Rainbow trout and Atlantic salmon are farmed in the same regions, which contributes to the understanding of the phylogenetics of PRV. This study supports the original classification using the PRV genomic segment S1, which differentiates PRV into two major genotypes, I and II, and each of them into two subgenotypes designated as Ia and Ib, and IIa and IIb, respectively [10,15]. Through massive phylogenetic dendrograms with all available sequences, we provided additional evidence that strengthens this classification using new PRV S1 and M2 segment sequences. These dendrograms allowed us to classify, by subgenotype, already published sequences of S1 and M2 segments that were previously unclassified in the GenBank database. This assignment showed that a significant number of the publicly available sequences belong to the PRV-1 subtype (subgenotypes Ia and Ib) and, to a lesser extent, to the PRV-3 subtype (subgenotype IIa). However, the lack of sequences of PRV-2 (subgenotype IIb) points to the need to perform further analysis to understand the contribution of this variant to the worldwide PRV phylogeny. Furthermore, this kind of study should be complemented by using methodologies based on whole genome sequencing.
The distance analyses of the PRV genomic M2 and S1 segments and the amino acid sequences of proteins σ3 (S1 segment) and μ1 (M2 segment) show clear clustering between sequences belonging to the same subgenotype. This contrasts with the PRV-1 subtype classification [20], which shows, in protein σ3 alignments, two marked clusters belonging to subgenotypes Ia and Ib, thereby strengthening the subgenotype classification. Although no country- or host-specific PRV subgenotype has been described, some clusters were observed in the subset of sequences used. These latter observations indicate the need to perform a deep meta-analysis to elucidate whether relationships between virulence and subgenotypes can be described. In fact, the PRV subgenotypes Ia and Ib (PRV-1) have been associated with HSMI in Atlantic salmon, this being a trending topic in aquaculture surveillance. In farmed Atlantic salmon, PRV-Ib has been shown to be the cause of HSMI [3]—whereas experimental studies have shown PRV-Ia to be of low virulence, although it is associated with moderate heart lesions in sockeye salmon [31]—and in Pacific Canada the associated disease in farmed Atlantic salmon was referred to as HSMI-like disease [33]. The higher virulence of PRV-Ib was linked to the evolution of PRV in Norwegian salmonid aquaculture possibly resulting from gene reassortment involving PRV segments S1 and M2 [34]. This prompted the development of a real-time reverse transcription polymerase chain reaction (RT-qPCR) assay specific for PRV-Ib associated with HSMI [35]. However, the study linking PRV-Ib to higher virulence [34] was largely speculative because, although eight of the PRV segment S1 sequences used in that study were of PRV-Ia from HSMI in BC Canada [13,14], the so-called “HSMI Clade” included only PRV-Ib isolates. Moreover, some of the PRV-Ib isolates have not been associated with HSMI [6,36]. In fact, an experimental study comparing the virulence of three PRV-1b isolates and three PRV-1a isolates in farmed Atlantic salmon showed no differences in PRV RNA in blood cells among the different isolates when measured by RT-qPCR from one to ten weeks post-challenge. However, two of the PRV-1b isolates resulted in more severe heart histopathology scores than the other PRV-1b isolate and the three PRV-1a isolates, although even in these PRV isolates with less severe histopathological lesions, a few individual fish were observed to have lesion scores high enough to qualify as HSMI [37]. It therefore follows that the variation in virulence found in PRV-1b isolates would be expected to also occur among PRV-1a isolates. The identification of the sequence relationships and the subgenotype classification performed in the present study can be used to enhance research on the virulence of PRV.

4. Materials and Methods

4.1. Extensive Phylogeny of PRV Segments

Sequences in FAST-All (FASTA) format for all PRV segments were retrieved from the NCBI GenBank public database [38]. Sequences from 2007 to the present date (May 2020) were selected by means of the rentrez package [39] under R environment v3.5.2 [40], utilizing the NCBI Taxonomy ID: txid1157337 assigned to PRV. To increase the accuracy of our analysis, multiple sequence alignments (MSAs) [41] were performed using the T-Coffee software [42]. Additionally, IQ-TREE, software to calculate trees using the maximum likelihood algorithm [43], was used to find the best model substitution automatically through its ModelFinder feature [44]. Trees were built using the ggtree [45] and ggplot2 packages [46], both based on the R statistical language [40] and scripted utilizing the integrated development environment R-Studio [47]. To graphically show the distance between σ3 and μ1 amino acidic sequences, the nucleotide sequences were translated to amino acidic sequences using Blastx [48]. Heatmaps were drawn using the distance matrix calculated by the phangorn R library [49] from the alignments obtained by the maximum likelihood algorithm using the Clustal Omega software, including the model substitution parameter [50]. The resulting distance information was postprocessed and formatted using the library reshape2 [51] and represented using the pheatmap package [52] under the R statistical language.

4.2. Nucleic Acid Sequencing of New Isolates

Viral sequencing was performed following the methods previously described to sequence genomic PRV segments [10,15]. Briefly, the RNA was isolated using a modified extraction protocol based on equilibrated phenol and stabilized chloroform, isoamyl alcohol (25:24:1) (PanReac Applichem, Chicago, IL, USA), and repurified using the E.Z.N.A.® Total RNA Kit I ( Omega Bio-tek, Inc., Norcross, GA, USA). The RNA was tested by RT-qPCR to detect the presence of PRV, using ELF-1α as internal control, following the protocol described by Kibenge et al. [10]. Only samples with a cycle threshold (Ct) below 25 were utilized to perform the DNA sequencing of segments S1 and M2 using the primers described in Kibenge et al. [10]. Sequences were deposited in the NCBI GenBank database under the accession numbers listed in Table 2.

Supplementary Materials

The following are available online at https://www.mdpi.com/2076-0817/10/1/41/s1, Figure S1: Rectangular version of phylogenetic distances from the alignment of S1 sequences above 400 nucleotides; Figure S2: Phylogenetic trees of segments L1, L2, L3, M1, M3, S2, S3 and S4; Figure S3: Phylogenetic distances from the alignment of 40 μ1 amino acid sequences grouped by genotype and represented as heatmaps; Figure S4: Phylogenetic distances from the alignment of 61 σ3 amino acid sequences grouped by genotype and represented as heatmaps; Figure S5: Phylogenetic distances from the alignment of μ1 sequences in Figure S4, grouped by country and represented as heatmaps; Figure S6: Phylogenetic distances from the alignment of σ3 sequences in Figure S4, grouped by host and represented as heatmaps; Table S1: GenBank Accession Numbers and metadata for PRV sequences of all ten genomic segments listed under the Tax ID: txid1157337 in the NCBI GenBank database to date (May 2020). Accession numbers and metadata of all sequences retrieved from NCBI GenBank used in this study; Table S2: GenBank accession numbers and metadata of PRV segments S1 and M2 sequences listed under the Tax ID: txid1157337 in the NCBI GenBank database to date (May 2020) used in this study.

Author Contributions

Conceptualization, D.A.M. and M.G.; methodology, D.A.M., M.G., S.V., R.S., J.R., M.K., F.K., Y.W.; software, S.V.; validation, Y.W., F.K. and M.K.; formal analysis, S.V.; investigation, M.G., R.S.; resources, M.G.; data curation, D.A.M.; writing—original draft preparation, D.A.M., M.G.; writing—review and editing, D.A.M., M.G., S.V., R.S., J.R., M.K., F.K., Y.W.; visualization, D.A.M.; supervision, M.G.; project administration, M.G.; funding acquisition, M.G. All authors have read and agreed to the published version of the manuscript.

Funding

This research was funded by the Centro de Investigaciones Biológicas Aplicadas (CIBA) and supported by the Facultad de Medicina Veterinaria, Laboratorio de Biotecnología, Universidad San Sebastián, Sede de La Patagonia.

Data Availability Statement

The data presented in this study are available in the Supplementary Materials.

Conflicts of Interest

The authors declare no conflict of interest. The funders had no role in the design of the study; in the collection, analyses, or interpretation of data; in the writing of the manuscript, or in the decision to publish the results.

References

  1. Palacios, G.; Løvoll, M.; Tengs, T.; Hornig, M.; Hutchison, S.; Hui, J.; Kongtorp, R.-T.; Savji, N.; Bussetti, A.V.; Solovyov, A.; et al. Heart and Skeletal Muscle Inflammation of Farmed Salmon Is Associated with Infection with a Novel Reovirus. PLoS ONE 2010, 5, e11487. [Google Scholar] [CrossRef]
  2. Kongtorp, R.T.; 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. [Google Scholar] [CrossRef]
  3. Wessel, Ø.; Braaen, S.; Alarcón, M.; Haatveit, H.; Roos, N.; Markussen, T.; Tengs, T.; Dahle, M.K.; Rimstad, E. Infection with purified Piscine orthoreovirus demonstrates a causal relationship with heart and skeletal muscle inflammation in Atlantic salmon. PLoS ONE 2017, 12, e0183781. [Google Scholar] [CrossRef] [Green Version]
  4. Markussen, T.; Dahle, M.K.; Tengs, T.; Løvoll, M.T.; Finstad, Ø.W.; Wiik-Nielsen, C.R.; Grove, S.; Lauksund, R.S.; Robertsen, B.; Rimstad, E. Sequence Analysis of the Genome of Piscine Orthoreovirus (PRV) Associated with Heart and Skeletal Muscle Inflammation (HSMI) in Atlantic Salmon (Salmo salar). PLoS ONE 2013, 8, e70075. [Google Scholar] [CrossRef]
  5. Wiik-Nielsen, C.; Løvoll, M.; Sandlund, N.; Faller, R.; Wiik-Nielsen, J.; Jensen, B.B. First detection of piscine reovirus (PRV) in marine fish species. Dis. Aquat. Org. 2012, 97, 255–258. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  6. 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, 36, 483–493. [Google Scholar] [CrossRef] [PubMed]
  7. Olsen, A.B.; Hjortaas, M.; Tengs, T.; Hellberg, H.; Johansen, R. First Description of a New Disease in Rainbow Trout (Oncorhynchus mykiss (Walbaum)) Similar to Heart and Skeletal Muscle Inflammation (HSMI) and Detection of a Gene Sequence Related to Piscine Orthoreovirus (PRV). PLoS ONE 2015, 10, e0131638. [Google Scholar] [CrossRef] [PubMed]
  8. Madhun, A.S.; Isachsen, C.H.; Omdal, L.M.; Einen, A.C.B.; Bjørn, P.A.; Nilsen, R.; Karlsbakk, E. Occurrence of salmonid alphavirus (SAV) and piscine orthoreovirus (PRV) infections in wild sea trout Salmo trutta in Norway. Dis. Aquat. Org. 2016, 120, 109–113. [Google Scholar] [CrossRef] [Green Version]
  9. Madhun, A.S.; Isachsen, C.H.; Omdal, L.M.; Einen, A.C.B.; Maehle, S.; Wennevik, V.; Niemelä, E.; Svåsand, T.; Karlsbakk, E. Prevalence of piscine orthoreovirus and salmonid alphavirus in sea-caught returning adult Atlantic salmon (Salmo salar L.) in northern Norway. J. Fish Dis. 2018, 41, 797–803. [Google Scholar] [CrossRef] [Green Version]
  10. Kibenge, M.J.T.; Iwamoto, T.; Wang, Y.; Morton, A.; Godoy, M.; Kibenge, F.S.B. 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. 2013, 10, 230. [Google Scholar] [CrossRef] [Green Version]
  11. Marty, G.D.; Morrison, D.B.; Bidulka, J.; Joseph, T.; Siah, A. Piscine reovirus in wild and farmed salmonids in British Columbia, Canada: 1974–2013. J. Fish Dis. 2015, 38, 713–728. [Google Scholar] [CrossRef]
  12. Siah, A.; Morrison, D.B.; Fringuelli, E.; Savage, P.S.; Richmond, Z.; Johns, R.; Purcell, M.K.; Johnson, S.C.; Saksida, S.M. Piscine Reovirus: Genomic and Molecular Phylogenetic Analysis from Farmed and Wild Salmonids Collected on the Canada/US Pacific Coast. PLoS ONE 2015, 10, e0141475. [Google Scholar] [CrossRef] [PubMed]
  13. Di Cicco, E.; Ferguson, H.W.; Schulze, A.D.; Kaukinen, K.H.; Li, S.; Vanderstichel, R.; Wessel, Ø.; Rimstad, E.; Gardner, I.A.; Hammell, K.L.; et al. Heart and skeletal muscle inflammation (HSMI) disease diagnosed on a British Columbia salmon farm through a longitudinal farm study. PLoS ONE 2017, 12, e0171471. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  14. Di Cicco, E.; Ferguson, H.W.; Kaukinen, K.H.; Schulze, A.D.; Li, S.; Tabata, A.; Günther, O.P.; Mordecai, G.; Suttle, C.A.; Miller, K.M. The same strain of Piscine orthoreovirus (PRV-1) is involved in the development of different, but related, diseases in Atlantic and Pacific Salmon in British Columbia. Facets 2018, 3, 599–641. [Google Scholar] [CrossRef] [Green Version]
  15. Godoy, M.G.; Kibenge, M.J.T.; Wang, Y.; Suarez, R.; Leiva, C.; Vallejos, F.; Kibenge, F.S.B. First description of clinical presentation of piscine orthoreovirus (PRV) infections in salmonid aquaculture in Chile and identification of a second genotype (Genotype II) of PRV. Virol. J. 2016, 13, 98. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  16. Purcell, M.K.; Powers, R.L.; Evered, J.; Kerwin, J.; Meyers, T.R.; Stewart, B.; Winton, J.R. Molecular testing of adult Pacific salmon and trout (Oncorhynchus spp.) for several RNA viruses demonstrates widespread distribution of piscine orthoreo-virus (PRV) in Alaska and Washington. J. Fish Dis. 2018, 41, 347–355. [Google Scholar] [CrossRef] [PubMed]
  17. Takano, T.; Nawata, A.; Sakai, T.; Matsuyama, T.; Ito, T.; Kurita, J.; Terashima, S.; Yasuike, M.; Nakamura, Y.; Fujiwara, A.; et al. Full-Genome Sequencing and Confirmation of the Causative Agent of Erythrocytic Inclusion Body Syndrome in Coho Salmon Identifies a New Type of Piscine Orthoreovirus. PLoS ONE 2016, 11, e0165424. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  18. Gunnarsdóttir, H.M. Screening for Viruses in Atlantic Salmon in Iceland. ISAV, PMCV and PRV. 2017. Available online: https://skemman.is/handle/1946/27629 (accessed on 6 January 2021).
  19. Kuehn, R.; Stoeckle, B.C.; Young, M.; Popp, L.; Taeubert, J.-E.; Pfaffl, M.W.; Geist, J. Identification of a piscine reovirus-related pathogen in proliferative darkening syndrome (PDS) infected brown trout (Salmo trutta fario) using a next-generation technology detection pipeline. PLoS ONE 2018, 13, e0206164. [Google Scholar] [CrossRef] [Green Version]
  20. Dhamotharan, K.; Vendramin, N.; Markussen, T.; Wessel, Ø.; Cuenca, A.; Nyman, I.B.; Olsen, A.B.; Tengs, T.; Dahle, M.K.; Rimstad, E. Molecular and Antigenic Characterization of Piscine orthoreovirus (PRV) from Rainbow Trout (Oncorhynchus mykiss). Viruses 2018, 10, 170. [Google Scholar] [CrossRef] [Green Version]
  21. Vendramin, N.; Cuenca, A.; Sørensen, J.; Alencar, A.L.F.; Christiansen, D.H.; Jacobsen, J.A.; Axen, C.; Lieffrig, F.; Ruane, N.; Martin, P.; et al. Presence and genetic variability of Piscine orthoreovirus genotype 1 (PRV-1) in wild salmonids in Northern Europe and North Atlantic Ocean. J. Fish Dis. 2019, 42, 1107–1118. [Google Scholar] [CrossRef]
  22. Bigarré, L.; Boitard, P.M.; Labrut, S.; Jamin, M. Point sur la Situation Epidemiologique du Syndrome HSMI Chez les Salmonides en France. Bulletin Epidémiologique. Santé Animale-Alimentation. 2018. Available online: https://be.anses.fr/sites/default/files/O-026_2018-06-22_HSMI_bigarre-VF.pdf (accessed on 6 January 2021).
  23. Adamek, M.; Hellmann, J.; Flamm, A.; Teitge, F.; Vendramin, N.; Fey, D.; Riße, K.; Blakey, F.; Rimstad, E.; Steinhagen, D. Detection of piscine orthoreoviruses (PRV-1 and PRV-3) in Atlantic salmon and rainbow trout farmed in Germany. Transbound. Emerg. Dis. 2019, 66, 14–21. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  24. Pojezdal, Ľ.; Adamek, M.; Syrová, E.; Steinhagen, D.; Minářová, H.; Papežíková, I.; Seidlova, V.; Reschová, S.; Palíková, M. Health Surveillance of Wild Brown Trout (Salmo trutta fario) in the Czech Republic Revealed a Coexistence of Proliferative Kidney Disease and Piscine Orthoreovirus-3 Infection. Pathogens 2020, 9, 604. [Google Scholar] [CrossRef] [PubMed]
  25. Kongtorp, R.T.; Taksdal, T.; Lyngøy, A. Pathology of heart and skeletal muscle inflammation (HSMI) in farmed Atlantic salmon Salmo salar. Dis. Aquat. Org. 2004, 59, 217–224. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  26. Ferguson, H.W.; Kongtorp, R.T.; Taksdal, T.; Graham, D.; Falk, K. An outbreak of disease resembling heart and skeletal muscle inflammation in Scottish farmed salmon, Salmo salar L., with observations on myocardial regeneration. J. Fish Dis. 2005, 28, 119–123. [Google Scholar] [CrossRef] [PubMed]
  27. Finstad, Ø.W.; Dahle, M.K.; Lindholm, T.H.; Nyman, I.B.; Løvoll, M.; Wallace, C.; Olsen, C.M.; Storset, A.K.; Rimstad, E. Piscine orthoreovirus (PRV) infects Atlantic salmon erythrocytes. Vet. Res. 2014, 45, 1–13. [Google Scholar] [CrossRef] [Green Version]
  28. Wessel, Ø.; Olsen, C.M.; Rimstad, E.; Dahle, M.K. Piscine orthoreovirus (PRV) replicates in Atlantic salmon (Salmo salar L.) erythrocytes ex vivo. Vet. Res. 2015, 46, 26. [Google Scholar] [CrossRef] [Green Version]
  29. Garseth, Å.H.; Ekrem, T.; Biering, E. Phylogenetic Evidence of Long Distance Dispersal and Transmission of Piscine Reovirus (PRV) between Farmed and Wild Atlantic Salmon. PLoS ONE 2013, 8, e82202. [Google Scholar] [CrossRef]
  30. Kibenge, F.S.B. Emerging viruses in aquaculture. Curr. Opin. Virol. 2019, 34, 97–103. [Google Scholar] [CrossRef]
  31. Garver, K.A.; Marty, G.D.; Cockburn, S.N.; Richard, J.; Hawley, L.M.; Müller, A.; Thompson, R.L.; Purcell, M.K.; Saksida, S. Piscine reovirus, but not Jaundice Syndrome, was transmissible to Chinook Salmon, Oncorhynchus tshawytscha (Walbaum), Sockeye Salmon, Oncorhynchus nerka (Walbaum), and Atlantic Salmon, Salmo salar L. J. Fish Dis. 2016, 39, 117–128. [Google Scholar] [CrossRef]
  32. Siah, A.; Breyta, R.B.; Warheit, K.I.; Gagne, N.; Purcell, M.K.; Morrison, D.; Powell, J.F.F.; Johnson, S.C. Genomes reveal genetic diversity of Piscine orthoreovirus in farmed and free-ranging salmonids from Canada and USA. Virus Evol. 2020, 6, veaa054. [Google Scholar] [CrossRef]
  33. Polinski, M.P.; Marty, G.D.; Snyman, H.N.; Garver, K.A. Piscine othoreovirus demonstrates high infectivity but low virulence in Atlantic salmon of Pacific Canada. Sci. Rep. 2019, 9, 3297. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  34. Dhamotharan, K.; Tengs, T.; Wessel, Ø.; Braaen, S.; Nyman, I.B.; Hansen, E.F.; Christiansen, D.H.; Dahle, M.K.; Rimstad, E.; Markussen, T. Evolution of the Piscine orthoreovirus Genome Linked to Emergence of Heart and Skeletal Muscle Inflammation in Farmed Atlantic Salmon (Salmo salar). Viruses 2019, 11, 465. [Google Scholar] [CrossRef] [Green Version]
  35. Siah, A.; Knutsen, E.; Richmond, Z.; Mills, M.; Frisch, K.; Powell, J.F.F.; Brevik, Ø.; Duesund, H. Real-time RT-qPCR assay to detect sequences in the Piscine orthoreovirus-1 genome segment S1 associated with heart and skeletal muscle inflammation in Atlantic salmon. J. Fish Dis. 2020, 43, 955–962. [Google Scholar] [CrossRef] [PubMed]
  36. Løvoll, M.; Alarcón, M.; Jensen, B.B.; Taksdal, T.; Kristoffersen, A.B.; Tengs, T. Quantification of piscine reovirus (PRV) at different stages of Atlantic salmon Salmo salar production. Dis. Aquat. Org. 2012, 99, 7–12. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  37. Wessel, Ø.; Hansen, E.F.; Dahle, M.K.; Alarcon, M.; Vatne, N.A.; Nyman, I.B.; Soleim, K.B.; Dhamotharan, K.; Timmerhaus, G.; Markussen, T.; et al. Piscine Orthoreovirus-1 Isolates Differ in Their Ability to Induce Heart and Skeletal Muscle Inflammation in Atlantic Salmon (Salmo salar). Pathogens 2020, 9, 1050. [Google Scholar] [CrossRef]
  38. Clark, K.; Karsch-Mizrachi, I.; Lipman, D.J.; Ostell, J.; Sayers, E.W. GenBank. Nucleic Acids Res. 2016, 44, D67–D72. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  39. Winter, D.J. Rentrez: An R package for the NCBI eUtils. API R J. 2017, 9, 520–526. [Google Scholar] [CrossRef]
  40. R Core Team. R: A Language and Environment for Statistical Computing; Foundation for Statistical Computing: Vienna, Austria, 2013. [Google Scholar]
  41. Wheeler, T.J.; Kececioglu, J.D. Multiple alignment by aligning alignments. Bioinformatics 2007, 23, i559–i568. [Google Scholar] [CrossRef] [PubMed]
  42. Notredame, C.; Higgins, D.G.; Heringa, J. T-coffee: A novel method for fast and accurate multiple sequence alignment. J. Mol. Biol. 2000, 302, 205–217. [Google Scholar] [CrossRef] [Green Version]
  43. Nguyen, L.-T.; Schmidt, H.A.; Von Haeseler, A.; Minh, B.Q. IQ-TREE: A Fast and Effective Stochastic Algorithm for Estimating Maximum-Likelihood Phylogenies. Mol. Biol. Evol. 2015, 32, 268–274. [Google Scholar] [CrossRef]
  44. Kalyaanamoorthy, S.; Minh, B.Q.; Wong, T.K.F.; Von Haeseler, A.; Jermiin, L.S. ModelFinder: Fast model selection for accurate phylogenetic estimates. Nat. Methods 2017, 14, 587–589. [Google Scholar] [CrossRef] [Green Version]
  45. Yu, G. Using ggtree to Visualize Data on Tree-Like Structures. Curr. Protoc. Bioinform. 2020, 69, e96. [Google Scholar] [CrossRef] [PubMed]
  46. Ginestet, C.E. ggplot2: Elegant Graphics for Data Analysis. J. R. Stat. Soc. Ser. A Stat. Soc. 2011, 174, 245–246. [Google Scholar] [CrossRef]
  47. Booth, D.S.; Szmidt-Middleton, H.; King, N.; Westbrook, M.J.; Young, S.L.; Kuo, A.; Abedin, M.; Chapman, J.; Fairclough, S.R.; Hellsten, U.; et al. RStudio: Integrated Development for R. Nature 2018. [Google Scholar] [CrossRef]
  48. McGinnis, S.; Madden, T.L. BLAST: At the core of a powerful and diverse set of sequence analysis tools. Nucleic Acids Res. 2004, 32, W20–W25. [Google Scholar] [CrossRef] [PubMed]
  49. Schliep, K.P. phangorn: Phylogenetic analysis in R. Bioinformatics 2010, 27, 592–593. [Google Scholar] [CrossRef] [Green Version]
  50. Larkin, M.A.; Blackshields, G.; Brown, N.P.; Chenna, R.; Mcgettigan, P.A.; McWilliam, H.; Valentin, F.; Wallace, I.M.; Wilm, A.; Lopez, R.; et al. Clustal W and Clustal X version 2.0. Bioinformatics 2007, 23, 2947–2948. [Google Scholar] [CrossRef] [Green Version]
  51. Wickham, H. Reshaping Data with the reshape Package. J. Stat. Softw. 2007, 21, 1–20. [Google Scholar] [CrossRef]
  52. Kolde, R. Pheatmap: Pretty Heatmaps; R Package (Version 1.0.8). 2015. Available online: https://cran.r-project.org/web/packages/pheatmap/index.html (accessed on 6 January 2019).
Pathogens 10 00041 g001 550
Figure 1. Circular dendrogram of S1 PRV genomic segment. The phylogenetic analysis was performed using sequences with at least 400 nucleotides. New and unknown subgenotype sequences were assigned according to their positions in the dendrogram branches.
Figure 1. Circular dendrogram of S1 PRV genomic segment. The phylogenetic analysis was performed using sequences with at least 400 nucleotides. New and unknown subgenotype sequences were assigned according to their positions in the dendrogram branches.
Pathogens 10 00041 g001
Pathogens 10 00041 g002 550
Figure 2. Phylogenetic tree of PRV M2 genomic segment. The dendrogram was built utilizing sequences with at least 1000 nucleotides. New and unknown subgenotype sequences were assigned according to their positions in the dendrogram branches.
Figure 2. Phylogenetic tree of PRV M2 genomic segment. The dendrogram was built utilizing sequences with at least 1000 nucleotides. New and unknown subgenotype sequences were assigned according to their positions in the dendrogram branches.
Pathogens 10 00041 g002
Pathogens 10 00041 g003 550
Figure 3. Phylogenetic distance heatmap for 14 PRV M2 genomic sequences. The colors represent the subgenotype and the distance matrix value.
Figure 3. Phylogenetic distance heatmap for 14 PRV M2 genomic sequences. The colors represent the subgenotype and the distance matrix value.
Pathogens 10 00041 g003
Pathogens 10 00041 g004 550
Figure 4. Phylogenetic distance heatmap obtained from 17 PRV S1 genomic sequences. The colors represent the subgenotype and the distance matrix value.3. Discussion.
Figure 4. Phylogenetic distance heatmap obtained from 17 PRV S1 genomic sequences. The colors represent the subgenotype and the distance matrix value.3. Discussion.
Pathogens 10 00041 g004
Table
Table 1. Number of sequences for each Piscine orthoreovirus (PRV) genomic segment used in the study.
Table 1. Number of sequences for each Piscine orthoreovirus (PRV) genomic segment used in the study.
PRV Genomic Segment
L1L2L3M1M2M3S1S2S3S4
Number of sequences53384041614339010142100
Table
Table 2. New sequences of PRV segments S1 and M2 produced in this study.
Table 2. New sequences of PRV segments S1 and M2 produced in this study.
GenBank Accession NumberPRV SegmentPRV SubgenotypeCollection DateFish Species 1Location in ChileFarming Conditions
MT598071M2Ib06-04-2018O. kisutchChiloé (Calen)Seawater
MT598072M2Ia12-04-2018S. salarChiloé (Quellón)Seawater
MT598073M2Ia26-04-2018O. kisutchChiloé (Quellón)Seawater
MT598068M2Ia12-06-2018O. kisutchChiloé (Quellón)Seawater
MT598074M2Ia19-06-2018O. kisutchChiloé (Caucahue)Seawater
MT598069M2Ia21-06-2018O. kisutchChiloé (Detico)Seawater
MT598075M2Ib27-03-2018O. kisutchPargua (P. Montt) 2Freshwater
MT598076M2Ib05-12-2018S. salarChiloé (Voigue)Seawater
MT598077M2Ib19-12-2019S. salarLencarFeshwater
MT598078M2Ib19-12-2019S. salarLencaFreshwater
MT598070M2Ib04-12-2018S. salarChiloé (Quicavi)Seawater
MT598079S1Ib06-04-2018O. kisutchChiloé (Calen)Seawater
MT598080S1Ia12-04-2018S. salarChiloé (Quellón)Seawater
MT598081S1Ia26-04-2018O. kisutchChiloé (Quellón)Seawater
MT598082S1Ia12-06-2018O. kisutchChiloé (Quellón)Seawater
MT598083S1Ia19-06-2018O. kisutchChiloé (Caucahue)Seawater
MT598084S1Ia21-06-2018O. kisutchChiloé (Detico)Seawater
MT598085S1Ib27-03-2018O. kisutchPargua (P. Montt) 2Freshwater
MT598086S1Ib04-12-2018S. salarChiloé (Quicavi)Seawater
MT598087S1Ib05-12-2018S. salarChiloé (Voigue)Seawater
MT598088S1Ib17-06-2019S. salarPargua (P. Montt) 2Seawater
MT598089S1Ia07-08-2019O. kisutchPuerto CisneSeawater
1 O. kisutch (Oncorhynchus kisutch); S. salar (Salmo salar); 2 P. Montt, Puerto Montt.
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Godoy, M.; Medina, D.A.; Suarez, R.; Valenzuela, S.; Romero, J.; Kibenge, M.; Wang, Y.; Kibenge, F. Extensive Phylogenetic Analysis of Piscine Orthoreovirus Genomic Sequences Shows the Robustness of Subgenotype Classification. Pathogens 2021, 10, 41. https://doi.org/10.3390/pathogens10010041

AMA Style

Godoy M, Medina DA, Suarez R, Valenzuela S, Romero J, Kibenge M, Wang Y, Kibenge F. Extensive Phylogenetic Analysis of Piscine Orthoreovirus Genomic Sequences Shows the Robustness of Subgenotype Classification. Pathogens. 2021; 10(1):41. https://doi.org/10.3390/pathogens10010041

Chicago/Turabian Style

Godoy, Marcos, Daniel A. Medina, Rudy Suarez, Sandro Valenzuela, Jaime Romero, Molly Kibenge, Yingwei Wang, and Frederick Kibenge. 2021. "Extensive Phylogenetic Analysis of Piscine Orthoreovirus Genomic Sequences Shows the Robustness of Subgenotype Classification" Pathogens 10, no. 1: 41. https://doi.org/10.3390/pathogens10010041

APA Style

Godoy, M., Medina, D. A., Suarez, R., Valenzuela, S., Romero, J., Kibenge, M., Wang, Y., & Kibenge, F. (2021). Extensive Phylogenetic Analysis of Piscine Orthoreovirus Genomic Sequences Shows the Robustness of Subgenotype Classification. Pathogens, 10(1), 41. https://doi.org/10.3390/pathogens10010041

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

Article Metrics

Back to TopTop