Next Article in Journal
Feline Foamy Virus is Highly Prevalent in Free-Ranging Puma concolor from Colorado, Florida and Southern California
Next Article in Special Issue
Susceptibility of Exopalaemon carinicauda to the Infection with Shrimp Hemocyte Iridescent Virus (SHIV 20141215), a Strain of Decapod Iridescent Virus 1 (DIV1)
Previous Article in Journal
Reassessment of Viroid RNA Cytosine Methylation Status at the Single Nucleotide Level
Previous Article in Special Issue
Description of a Natural Infection with Decapod Iridescent Virus 1 in Farmed Giant Freshwater Prawn, Macrobrachium rosenbergii
 
 
Search for Articles:
Title / Keyword
Author / Affiliation / Email
Journal
Article Type
 
 
Advanced Search
 
Section
Special Issue
Volume
Issue
Number
Page
 
Journals
Viruses
Volume 11
Issue 4
10.3390/v11040358
viruses-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
first_page
settings
Order Article Reprints
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Distribution and Phylogeny of Erythrocytic Necrosis Virus (ENV) in Salmon Suggests Marine Origin

by
Veronica A. Pagowski
1,†,
Gideon J. Mordecai
1,*,†,
Kristina M. Miller
2,*,
Angela D. Schulze
2,
Karia H. Kaukinen
2,
Tobi J. Ming
2,
Shaorong Li
2,
Amy K. Teffer
3,
Amy Tabata
2 and
Curtis A. Suttle
1,4,5,6,*
1
Department of Earth, Ocean and Atmospheric Sciences, University of British Columbia, 2207 Main Mall, Vancouver, BC V6T 1Z4, Canada
2
Pacific Biological Station, Fisheries and Oceans Canada, 3190 Hammond Bay Rd, Nanaimo, BC V9T 6N7, Canada
3
Biology Department, University of Victoria, 3800 Finnerty Rd, Victoria, BC V8P 5C2, Canada
4
Department of Microbiology and Immunology, University of British Columbia, 2207 Main Mall, Vancouver, BC V6T 1Z4, Canada
5
Department of Botany, University of British Columbia, 2207 Main Mall, Vancouver, BC V6T 1Z4, Canada
6
Institute for the Oceans and Fisheries, University of British Columbia, 2207 Main Mall, Vancouver, BC V6T 1Z4, Canada
*
Authors to whom correspondence should be addressed.
†
These authors contributed equally to this work.
Viruses 2019, 11(4), 358; https://doi.org/10.3390/v11040358
Submission received: 17 March 2019 / Revised: 12 April 2019 / Accepted: 16 April 2019 / Published: 18 April 2019
(This article belongs to the Special Issue Family Iridoviridae: Molecular and Ecological Studies of a Family Infecting Invertebrates and Ectothermic Vertebrates)
Browse Figures
Review Reports Versions Notes

Abstract

:
Viral erythrocytic necrosis (VEN) affects over 20 species of marine and anadromous fishes in the North Atlantic and North Pacific Oceans. However, the distribution and strain variation of its viral causative agent, erythrocytic necrosis virus (ENV), has not been well characterized within Pacific salmon. Here, metatranscriptomic sequencing of Chinook salmon revealed that ENV infecting salmon was closely related to ENV from Pacific herring, with inferred amino-acid sequences from Chinook salmon being 99% identical to those reported for herring. Sequence analysis also revealed 89 protein-encoding sequences attributed to ENV, greatly expanding the amount of genetic information available for this virus. High-throughput PCR of over 19,000 fish showed that ENV is widely distributed in the NE Pacific Ocean and was detected in 12 of 16 tested species, including in 27% of herring, 38% of anchovy, 17% of pollock, and 13% of sand lance. Despite frequent detection in marine fish, ENV prevalence was significantly lower in fish from freshwater (0.03%), as assessed with a generalized linear mixed effects model (p = 5.5 × 10−8). Thus, marine fish are likely a reservoir for the virus. High genetic similarity between ENV obtained from salmon and herring also suggests that transmission between these hosts is likely.

1. Introduction

Viral erythrocytic necrosis (VEN) is a disease associated with severe blood abnormalities in infected fish which has caused mass mortality in Pacific herring (Clupea pallasii) [1]. The disease is traditionally diagnosed by microscopic examination of stained blood smears for the presence of inclusion bodies within the cytoplasm of infected erythrocytes. Electron microscopy revealed that infected erythrocytes contained icosahedral virions, which were named erythrocytic necrosis virus (ENV) [2]. Although first described more than half a century ago, the virus is poorly characterized as attempts to propagate it in fish cell lines have been unsuccessful [3].
Based on the small amount of available ENV genomic sequence, the virus has been assigned to a new putative genus within the Iridoviridae (a family of double-stranded DNA viruses), comprising other erythrocytic viruses from ectothermic hosts [2,4]. Further genomic sequencing of ENV may provide a reference for studying genetic variation geographically, and among fish species. Currently, the only verified ENV sequences available in GenBank encode ATPase, the major capsid protein (MCP), DNA-dependent DNA polymerase, and DNA-dependent RNA polymerase.
In the NE Pacific Ocean, VEN has been described in Pacific herring [1] and in the marine phase of pink (Oncorhynchus gorbuscha), chum (Oncorhynchus keta), coho (Oncorhynchus kisutch), steelhead (Oncorhynchus mykiss), and Chinook (Oncorhynchus tshawytscha) salmon [5,6]. High geographic variability in VEN prevalence and disease susceptibility of chum, coho, sockeye, and Chinook salmon, as well as Pacific herring, suggests that ENV could help explain high year-to-year variability in the population dynamics of these keystone species throughout coastal regions of the NE Pacific Ocean [1,7,8]. In the laboratory, VEN has been induced after confining healthy salmon with diseased fish. Among salmonids, disease transmission has been demonstrated in pink and chum salmon [5,8,9]. Chinook, coho, and sockeye salmon appear to be more resistant to infection in challenge studies, as assessed by electron microscopy [5,9], despite natural epizootics in Chinook and coho salmon. To date, an explicit relationship between disease manifestation and viral load has not been described, and ENV isolation has not been accomplished in previous studies reporting disease transmission among species. There is substantial evidence that viruses in different genera within the family Iridoviridae cause different disease manifestation and severity [10,11,12,13,14,15]; thus, phylogenetic placement and genetic characterization of ENV may offer insight into factors contributing to disease outbreaks.
Considering the high prevalence of ENV in herring and salmon in the NE Pacific Ocean [1,16,17,18,19,20], there is a need to further characterize the virus in at-risk fish populations. Substantial herring stock declines in British Columbia since the 1970s, which have mainly been attributed to overfishing, are of concern as herring are a keystone species that underpin the coastal food web [21,22]. Similarly, Pacific salmon provide an important biological, economic, and cultural resource in British Columbia, and recent stock declines are likely to have large impacts on wildlife and fisheries in the region. The Pacific Salmon Commission, for example, reported a 60% reduction in Salish Sea Chinook abundance from 1984 to 2010 [23]. Furthermore, the Committee on the Status of Endangered Wildlife in Canada reports that half of the British Columbia Chinook salmon populations are endangered [24].
In the current study, we applied metatranscriptomic sequencing to examine the phylogenetic relationship of ENV in Chinook salmon and investigated the epidemiology and host tropism of ENV using high-throughput PCR. Our results show that ENV is widely prevalent across numerous fish species in the NE Pacific Ocean and that ENV in salmon and herring has high genetic similarity, suggesting that the virus may circulate among these species.

2. Materials and Methods

2.1. Fish Sampling

In total, 19,652 fish comprising 16 species were sampled from freshwater and marine environments as previously described [17,25]. Briefly, 3228 freshwater samples were collected at hatcheries, through beach seining, or by smolt traps for wild salmon. Marine samples were obtained mostly by purse and beach seines, and by trawl. Typically, mixed tissue samples were dissected from fish using sterile procedures in the field [26] and frozen in RNAlater before nucleic acid extraction, although some fish were flash frozen in the field and dissected in the laboratory. Both procedures have been routinely used for examining individual fish for the presence of viral nucleic acids [16,18,26,27,28]. Samples were collected over the course of an 11-year period from 2007–2018 in a region spanning Alaska to Northern Washington (Figure S1) as part of a large pathogen-screening effort conducted by Fisheries and Oceans Canada.

2.2. Data Collection

The occurrence and abundance of ENV in fish was determined using the Fluidigm BioMark Platform at the Department of Fisheries and Oceans Canada [26]. Briefly, the platform provides an estimate of viral load based on copy number, as assessed by RT-qPCR. Copy numbers are calculated based on serial dilutions of artificial construct DNA standards. The calculated limit of detection (LOD) was applied to identify fish with amplifications above the 95% detection threshold [26]. ENV assays on this platform show 100% inclusivity (detection of all known strain variants of the targeted microbe) and 97.9% exclusivity (no detection of untargeted microbe species). Primers were originally obtained from a 100 bp ATPase-like protein partial gene sequence. Additional details on primer sequences as well as assay specificity and reliability are available in Miller et al. [26]. ENV monitoring was conducted alongside assays for 46 other infective agents on combined RNA and DNA extractions. In this study, ENV prevalence is reported as the proportion of fish with ENV detections, both with and without the LOD criteria applied. When not mentioned explicitly, prevalence values are reported with LOD criteria applied. Statistical analyses were done only on samples with LOD criteria applied.

2.3. Metatranscriptomic Sequencing and Bioinformatics

In order to isolate putative ENV sequences, several samples with high-load detections of this virus, as assessed using the Fluidigm BioMark, were selected for transcriptomic sequencing. In total, three aquaculture and one wild Chinook salmon were used to predict ENV proteins. Three additional fish, including an Atlantic salmon, a Chinook salmon, and a herring sample underwent an enrichment step for these predicted ENV proteins before bioinformatic analysis.
All fish except the wild Chinook salmon were sequenced on Illumina Next-Generation Sequencing (NGS) platforms using different RNA-seq protocols, depending on original targets and loads. To avoid DNA contamination from the host reducing the sequencing depth of the target virus, RNA sequencing was used to obtain transcriptomic sequences of the DNA virus. Each of these libraries was single tissue, either heart or spleen. Ribosomal RNA was removed from total RNA using the RiboMinus Invitrogen Eukaryote kit for RNA (Life Technologies, Carlsbad, CA, USA). The RNA-Seq library was prepared using the NEBNext Ultra RNA Library prep kit (New England BioLabs, Ipswich, MA, USA) with an average fragment size of 250-bp and was paired-end sequenced with 100-bp reads on the Illumina HiSeq analyzer (Illumina, San Diego, CA, USA).
For the wild Chinook sample, the ENV contigs were obtained using a similar RNAseq approach. This library was created using pooled tissue (gill, liver, heart, kidney and brain) and prepared with the ScriptSeq Complete Epidemiology NGS library kit (Illumina, San Diego, CA, USA). Briefly, ribosomal RNA was removed from Total RNA using the Epicentre ScriptSeq Complete Gold Kit (Epidemiology) (Illumina, San Diego, CA, USA) according to the manufacturer’s instructions. The ScriptSeq Index reverse primers were added to the cDNA during the final amplification step which involved 14 cycles. Finally, a paired-end 125 bp sequencing run was performed on the Illumina HiSeq System.
Finally, in order to enhance our sensitivity (i.e., NGS read depth and coverage) for ENV in the medium load (~1300 copies) farmed Atlantic salmon sample, and the Chinook and herring samples (~42,670 and 214,500 copies, respectively) we employed SureSelectXT enrichment technology (Agilent, Santa Clara, CA, USA). A custom set of RNA target enrichment probes (120 bp in length and staggered along the exome of interest) were designed for ENV as well as many other salmonid viruses assessed on our infectious agent monitoring platform. These sequences (497.266-kbp in total length) and subsequent bait oligonucleotides included all of the suspected ENV contigs previously assembled from high load samples. Baits which failed the SureSelect QA/QC parameters and/or significantly matched salmonid genes via BLAST searches were removed, leaving the final set of enrichment probes at 20,497. The mixed tissue (gill, atrium, ventricle, liver, pyloric caeca, spleen, head kidney, posterior kidney) RNA library, was prepared using the SureSelectXT low input (NGS) target workflow (Agilent, Santa Clara, CA, USA) with a SureSelectXT RNA Direct/XTHS modified protocol. Approximately 200 ng of total RNA was lyophilized and fed into the SureSelect Strand-Specific RNA library Prep kit (Agilent, Santa Clara, CA, USA) according to the manufacturer’s instructions. After the 2nd strand cDNA synthesis and end repair steps, the library prep was moved to the SureSelectXT low input reagent kit, starting with the end repair and A tailing step. Molecular barcoded adaptors were added using ligation and then amplified for 14 cycles according to manufacturer’s instructions to create a pre-capture RNAseq library. Samples were quantified with the Qubit dsDNA HS kit (Invitrogen, Carlsbad, CA, USA), qualified with the DNA12000 chips run on the Agilent 2100 Bioanalyzer (Agilent, Santa Clara, CA, USA), and pooled into a batch of 12 prior to hybridizing 1500 ng to the bait library. Hybridizations were incubated at 65 ℃, captured on streptavidin beads (Beckman Coulter, Brea, CA, USA) and washed at 70 ℃ according to manufacturer’s instructions before a post-capture amplification of 14 cycles. These final libraries were quantified with the Qubit dsDNA HS kit (Invitrogen, Carlsbad, CA, USA) and qualified with the DNA HS chips run on the Agilent 2100 Bioanalyzer (Agilent, Santa Clara, CA, USA). Finally, a paired-end 101 bp v2 300 kit sequencing run was performed on the Illumina Miseq (Illumina, San Diego, CA, USA), which included a 5% phiX spike in.
After adapter removal, Illumina MiSeq sequencing produced between 53.3 and 59.9 M reads for each Chinook salmon used to predict ENV amino-acid sequences (quality score >28). These reads were processed as detailed below. Adapters were removed using Trimmomatic [29], and the trimmed reads aligned with genome sequence from Atlantic Salmon [30] using the Burrows-Wheeler Aligner [31]. Unmapped sequences were extracted from the dataset using Samtools [32] and assembled into contiguous sequences (contigs) using SPAdes [33]. The translated contigs were queried against the non-redundant (NR) database in GenBank using DIAMOND [34]. Contigs with top hits to members of the family Iridoviridae were extracted in Microsoft Excel and the Qiime script filter_fasta.py [35] was adapted to retain these contigs as fasta files. GeneMark [36] was used to make protein predictions for putative ENV nucleotide sequences. Predicted proteins were subject to a BLAST search against the NR database [37], the lymphocystis disease virus genome, and a set of 47 conserved genes within Nucleo-Cytoplasmic Large DNA viruses (NCLDVs) [38]. Additionally, contigs of 500 bp in length or greater were extracted and translated into all six frames using Geneious version 9.1.8 [39]. A single, most likely translation frame was selected from the BLAST results, and those with BLAST results mapping to ENV were used to create phylogenies for ATPase, DNA-dependent DNA polymerase, DNA-dependent RNA polymerase, and the MCP. For each of these proteins, phylogenetic trees were mapped using available sequences from closely related iridoviruses using ClustalW for alignments and PhyML with Le Gascuel substitution model for tree creation [40,41]. For each phylogenetic tree, Spodoptera frugiperda ascovirus 1a was used as an outgroup. Assembled contigs have been submitted to GenBank under the accession numbers MK638669–MK638757. To increase confidence that sequences we attributed to ENV belong to this virus, we aligned sequences obtained from the farmed Atlantic salmon (632,935 reads), Chinook salmon (658,241 reads), and herring (800,139 reads) which underwent an enrichment step for ENV viral content to our putative ENV protein-encoding sequences using the using the Burrows-Wheeler Aligner [31].

2.4. Spatial Epidemiology Analysis

All statistical analyses and plots were performed in R version 3.5.2 [42]; scripts are available in File S1. For statistical analyses, ENV detection was categorized as positive or negative with LOD criteria applied and viral loads were quantified based on estimated viral copy numbers in host fish, following Miller et al. [26]. For map plotting only, ENV load was quantified by subdividing ENV copy number data into six categorical bins based on value with no LOD criteria (0, 0–5, 5–10, 10–100, 100–10,000, >10,000 viral copies). As sample sizes for other fish were small, statistical analyses were only conducted on herring and salmon. However, ENV prevalence and load was investigated for all sampled fish species. Heat maps were produced with an inverse-distance weighting function, using R packages “ggmap” and “gstat” [43,44], while plots were constructed with “ggplot2” [45]. Differences in load among species, age class, and habitat type were assessed using Kruskal-Wallis and post-hoc pairwise Dunn tests with Benjamini-Hochberg adjusted p values for multiple comparisons. Differences in ENV prevalence among categorical variables were assessed using Chi-squared tests of independence. A post-hoc Fisher’s exact test and Chi-squared tests with Bonferroni correction were conducted for pairwise comparisons to determine ENV prevalence differences among species and years, respectively, with “rcompanion” [46]. Spearman correlation was conducted to examine correlations among monthly prevalence between species. To account for residual variation in the data, a generalized linear mixed-effects model with Laplace approximation was implemented to examine differences in ENV prevalence between fresh and saltwater, with catch region, species, age class, season, population (hatchery or wild), and year considered as random effects with the R package “mlmRev” [47,48]. A similar model was used to assess differences in prevalence among age classes (smolt or adult) with habitat type used as an additional random effect instead of age class.

3. Results

3.1. Genetic Characterization

Metatranscriptomic sequencing of Chinook salmon with high ENV loads revealed high sequence identity between ENV sequences found in Chinook salmon and viral sequences from GenBank associated with Pacific herring. The available ENV sequences from Pacific herring (ATPase, DNA-dependent DNA polymerase, MCP, and DNA-dependent RNA polymerase) showed over 99% nucleic acid identity to ENV sequences from Chinook salmon in our study (Table 1). Sequences reported in Table 1 originate from heart tissue of one of the aquaculture Chinook salmon samples and from the wild Chinook salmon mixed-tissue specimen, both with a high load ENV detection (CT values of 10.7 and 13.4, respectively). The aquaculture Chinook salmon specimen had jaundice and several co-infections including Paranucleospora theridion, Piscine reovirus, and Renibacterium salmoninarum, as determined by the RT-qPCR assay. Phylogenies based on these genes (Figure 1) from ENV and its relatives (Table 2), showed that viral sequences from herring and salmon form a well-supported clade. A metatranscriptomic approach on a DNA virus only reveals virally expressed transcripts; thus, a full ENV genome could not be assembled. Within the metatranscriptomic sequences, there were 32 putative ENV proteins which consistently mapped to proteins of similar function from viruses in the family Iridoviridae (Table S2), as well as BLAST hits to 21 of the 47 core proteins in NCLDVs [38] (File S2). Putative ENV transcripts typically had between 20% and 60% nucleotide identity with other iridoviruses. These values are consistent with nucleotide identities reported between existing ENV sequences and other fish iridoviruses. In total 117 contigs with BLAST hits to iridoviruses were isolated from transcriptomic sequencing. Putative ENV sequences ranged in length from 152 to 4475 bp (File S3). From these, we predicted 89 protein-encoding sequences longer than 200-bp (Accession numbers MK638669-MK638757). When putative ENV protein-encoding sequences were aligned to reads obtained from the Atlantic salmon, Chinook salmon, and Pacific herring enriched for these contigs, we detected 35 of 117 putative ENV sequences. Results from this analysis are summarized in file S4.

3.2. Spatial Epidemiology

ENV is widely distributed in the NE Pacific Ocean and was detected in 12 of 16 tested species (Figure 2) throughout the sampling region. High viral loads were common within the Strait of Georgia, along the west coast of Vancouver Island, and in straits and channels throughout coastal northern British Columbia and southern coastal Alaska (Figure 3). Over 19,000 fish were tested for ENV, from 16 different species collected from marine and fresh waters spanning from Washington to Alaska. ENV prevalence was highest in anchovy and herring, occurring in over 27% of all sampled fish in these species and in over 37% of smolts. In herring, the proportion of fish in which ENV was detected was significantly higher than in any salmon species tested. We also report significantly lower ENV prevalence and load within salmon smolt (pprev = 1.2 × 10−17, pload = 1.8 × 10−4) and greater ENV prevalence among herring smolt (p = 4.5 × 10−6), when compared to respective adult prevalence. Among fish with detections of ENV, viral load was highest in herring and Chinook salmon and lowest in Atlantic salmon. ENV load also appeared low in chum and pink salmon; however, these differences were not statistically significant (Figure S2, Table S3). Significant differences in ENV load were found between herring and coho, Chinook, sockeye, and Atlantic salmon (p values are given in Tables S3 and S4).
Despite high ENV prevalence and load in coastal regions, ENV was rarely detected, and percent prevalence was statistically lower in all species sampled in freshwater (p < 2.2 × 10−16) (Figure 4). Furthermore, after applying a generalized linear mixed-effects model with Laplace approximation to account for residual variation in the data, categorical classification of habitat type as freshwater or saltwater still had a significant effect on ENV prevalence (p = 5.50 × 10−8). Given that 98% of the freshwater detections were below the LOD, these all represent very low-load detections, and are unlikely biologically relevant or may represent “false” positives. As such, all statistical tests are reported for samples with LOD criteria applied only.

3.3. Seasonal and Yearly Variation

In Atlantic salmon, ENV prevalence showed seasonal variation, with the lowest prevalence occurring in August, for both smolt and adult fish (Figure S3). A similar seasonal trend was observed in smolt sockeye and Chinook salmon. A Spearman correlation coefficient of 0.87 was found between monthly prevalence of ENV in sockeye and Atlantic salmon (p = 0.004). Despite very low ENV prevalence in Atlantic salmon (2%), farmed coho and Chinook salmon showed higher ENV prevalence (29% and 48%, respectively) compared to their wild counterparts (3% and 7%, respectively) (Figure S4). This difference was significant in Chinook salmon (p = 2.08 × 10−50). Among salmon, changes in ENV prevalence from year to year appear to occur synchronously (Figure S5) and ENV prevalence also varies significantly by year (p < 2.2 × 10−16, Figure S6). A 63% decrease from the average prevalence of 5% is observed after 2013.

4. Discussion

The current study expands on sequence available for ENV and demonstrates that the virus is widespread in salmon and herring and is often present at high load in these fish. The low genetic variation among viruses infecting salmon and herring has implications for potential host range and the taxonomic classification of these viruses.
Nucleotide variation of ENV from herring and Chinook salmon was low, with ENV sequences from Pacific herring [2,4] having >99% nucleotide identity with the sequence obtained from Chinook salmon in this study, indicating that both viruses belong within the same genus, and likely the same species. High sequence similarity between herring and Chinook salmon could be suggestive of viral spillover between hosts. Moreover, the phylogenetic analysis clearly places ENV within the Iridoviridae, but distinct from other viruses within the family.
After enrichment for putative ENV sequences in the herring, Chinook salmon, and Atlantic salmon samples, we found 35 of 117 predicted protein-encoding sequences, suggesting that not all ENV transcripts are present in every infected fish. This may depend on expression levels and infection stage. Indeed, variable gene expression during different stages of infection occurs commonly among iridoviruses [49,50,51,52]. Moreover, our estimates of prevalence are conservative, as the assay is sensitive to sequence variation, so related strains of ENV could be missed.
ENV was widely distributed geographically and among fish species in marine waters of western North America. Despite high ENV prevalence and load in coastal marine environments, viral load and prevalence were consistently low in freshwater environments. Of the 3622 fish analyzed from freshwater, ENV was only detected in one Chinook salmon specimen after LOD criteria were applied. Interestingly, there were few ENV detections in fish from the open waters north of Vancouver Island (Figure 3); whereas ENV was common in coastal areas, i.e., channels, inlets, and straits along coastal British Columbia and southern Alaska (Figure 3). We hypothesize that interspecies transmission may be more likely in these areas, where there are higher densities of salmon and other fish. Greater ENV prevalence in aquaculture Pacific salmon than in wild counterparts (Figure S4) supports the idea that transmission of the virus may be more common in coastal and high-density environments. There is some support for this hypothesis in the literature as previous studies have shown that chum salmon may contract VEN via waterborne exposure [8] and that increased stock density is a predictor of disease progression in fish infected with other iridoviruses [53,54]. We also report significantly lower ENV prevalence among salmon smolts and in adult herring when compared to respective adult and smolt age classes for these species, suggesting that susceptibility to ENV varies by host species and age class. ENV prevalence may not, however, directly correlate with disease manifestation mediated by the virus, as species with low ENV load and prevalence (<0.5%) in this study, including chum and pink salmon, are more susceptible to VEN than other salmon species [55,56]. Thus, salmon species in which ENV prevalence and load were the greatest (Chinook and coho salmon) may be able to sustain higher viral loads and display fewer clinical symptoms.
We hypothesize that salmon likely contract the virus from marine reservoirs, given the low detection of ENV in freshwater salmon, high sequence similarity between ENV in salmon and herring, and high viral prevalence in several species of marine fish including herring, anchovy, pollock, and sand lance. Furthermore, the prevalence of ENV in Pacific salmon is much higher than in Atlantic salmon, suggesting that farmed Atlantic salmon are not a significant reservoir for ENV transmission, despite high-density rearing in this species. Previous research has shown that overall infectious agent diversity and burden increases when sockeye salmon enter the ocean [16], suggesting that ENV could contribute to increased infection stress experienced by out-migrating salmon smolts.
Salmon migration, which varies among species and populations, may help explain yearly and seasonal variation in ENV prevalence. Peaks in ENV prevalence occur during spring and late fall for Chinook and sockeye salmon, with significant drops seen during peak salmon river runs in July and August. A similar pattern is observed in Atlantic salmon, which are stationary throughout the year and therefore may provide a useful sentinel to study seasonal ENV dynamics in wild salmon. Indeed, there is a significant correlation in monthly ENV prevalence between Atlantic and sockeye salmon. However, variation in migration patterns among species and stocks complicates the analysis of seasonal ENV prevalence changes. Migration routes and timing among Chinook salmon stocks, for example, is highly variable [57,58]. Furthermore, several sockeye salmon stocks were sampled during different parts of the year. Increases in ENV prevalence in winter may also arise as a result of herring migrations to coastal regions during this time [59]. Overall, monthly ENV prevalence dynamics in salmon were similar to those observed previously in herring [1], further substantiating the hypothesis that interactions with herring promote infection dynamics in wild salmon.
Previously, little was known about the epidemiology of ENV in salmon and marine fish, as most studies focused on ENV in herring. Hershberger et al. detected ENV in up to 67% of herring, with similar seasonal variation to that which we observed in salmon, with the greatest proportion of fish testing positive for ENV in summer months [1]. They also reported that ENV epizootics can arise and dissipate spontaneously in geographically isolated regions along the North Pacific coastline. Additionally, Teffer et al. investigated ENV prevalence in returning Chinook salmon and detected ENV in 16% of tagged males and 25% of females in the Chilliwack River [20]. Together, our research and these studies indicate that the virus is widely distributed on the west coast of British Columbia and Alaska and that salmon are likely infected once they enter the ocean, with herring or other marine fish likely acting as a reservoir for ENV.

Implications

Detection of ENV has not been conclusively linked to disease onset and further studies are required to characterize this relationship. VEN is a poorly characterized disease in Chinook and sockeye salmon, yet it is relatively common in wild at-risk populations of these species. In contrast, ENV was relatively rare in pink and chum salmon, even though these species are more susceptible to VEN than Chinook and sockeye salmon in challenge studies [5]. However, relatively few pink (n = 222) and chum (n = 191) salmon were sampled in our study. There were significant differences in viral loads among species, with lowest mean ENV loads occurring in Atlantic salmon. Low viral loads and prevalence could indicate higher virulence, which may lower the chance of transmission and detection compared to persistently infected fish [60]. Alternatively, species with large ranges in ENV copy number, such as herring and Chinook salmon, could carry a persistent infection which becomes virulent at higher loads. These species may transmit the virus to susceptible species such as pink and chum salmon. Future challenge studies which further characterize Chinook salmon infection may elucidate whether infection dynamics appear similar to those of herring, which have the most similar distribution of viral load.
ENV prevalence was lower in salmon smolts and adult herring. Similarly, Hershberger et al. reported more frequent VEN epizootics in juvenile herring, compared to adults [1]. Presumably, lower viral prevalence in smolt salmon arises because fewer smolt have been exposed to marine waters, where we propose the virus originates. Previous studies reported that osmoregulatory stress, such as transitions between saline and freshwater environments, could be implicated in herring mortality in fish infected with VEN [15,61]. If viral infection does, indeed, impact osmoregulatory capacity and adaptation, the relatively high prevalence and load of ENV detected in salmon soon after ocean entry could diminish their ability to properly acclimate to changes in salinity in their environment.
Numerous studies [10,11,12,13,14] have reported greater disease severity caused by iridoviruses that are closely related to ENV when the temperature increases. It has been suggested that below 20 °C, iridoviruses may remain dormant in teleost hosts [14]. Other members of Iridoviridae that infect fish and show high infection mortality typically occur in warmer climates, such as Southeast Asia and Australia. Similarly, VEN progression is most severe during the summer in Pacific salmon [5]. Changes in temperature were not investigated in this study, but seasonal and yearly variation in ENV prevalence suggests that environmental variables, such as temperature, may be important. A significant drop in ENV prevalence following 2013 coincides with a shift to a positive Pacific Decadal Oscillation Index and warming temperatures in the study region [62].
It is possible that disease progression intensifies at warmer temperatures, such that fewer fish harboring the virus survive. This interpretation is consistent with a decrease in the overall prevalence of salmon infectious agents in the region from 2012 to 2013 reported by Nekouei et al. [16]. Alternatively, fish infected with the virus may be weakened or more susceptible to other diseases at suboptimal temperatures. In the context of climate change, this is an interesting avenue of future research, and directly relevant to salmon populations, as evidence suggests that increasing coastal and oceanic temperatures can have significant and detrimental effects on salmon migration and spawning [63,64]. If VEN has a temperature-dependent onset similar to other diseases caused by iridoviruses, ENV-mediated mortality could further stress at-risk populations of salmon and herring in the NE Pacific Ocean.

5. Conclusions

This research demonstrates that ENV is highly prevalent in the NE Pacific Ocean. Low ENV prevalence in freshwater, high prevalence in marine fish, and seasonal variability corresponding to marine migrations of salmon and herring suggest that ENV originates from the marine environment. High prevalence in several marine fish species suggests that the virus is endemic and that these species are reservoirs of the virus. Moreover, the similarity between ENV sequences from Chinook salmon and those from Pacific herring indicates that transmission between these species is possible. Finally, we present 89 new protein-encoding sequences attributed to ENV in this study.

Supplementary Materials

The following are available online at https://www.mdpi.com/1999-4915/11/4/358/s1, Figure S1: Sampling effort and ENV detection map:, Table S1: List of taxa abbreviations used in tables and figures, Table S2: ENV BLAST summary, Figure S2: ENV copy number in mixed tissue samples by species, Table S3: Summary of p values from pairwise tests for differences in load (L) and prevalence (P) among species, Table S4: Summary of p values and sample size for tests reported in this study, Figure S3: Seasonal ENV prevalence by species, Figure S4: ENV prevalence by population, Figure S5: ENV prevalence by year and species, Figure S6: ENV prevalence by year, File S1: R Scripts, File S2: NCVOG Summary, File S3: ENV Contigs, File S4: Enrichment Summary.

Author Contributions

Conceptualization, G.J.M., K.M.M., and C.A.S.; Methodology, V.A.P., G.J.M., A.D.S., K.H.K., T.J.M., S.L., and A.K.T.; Formal analysis, V.A.P.; Investigation, V.A.P., G.J.M., A.D.S.; Resources, K.M.M. and C.A.S.; Data curation, A.T.; Writing—original draft preparation, V.A.P.; Writing—review and editing, G.J.M., C.A.S., and K.M.M.; Visualization, V.A.P.; Supervision, G.J.M., K.M.M., and C.A.S.; Project administration, G.J.M., C.A.S., and K.M.M.

Funding

Funding was provided through the Strategic Salmon Health Initiative co-funded by the Pacific Salmon Foundation, Genome British Columbia, Fisheries and Oceans Canada, and MITACS, as well as the NSERC summer 2018 Undergraduate Student Research Award.

Acknowledgments

We would like to thank Christoph Deeg for feedback on phylogenetic analyses and insightful suggestions.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Hershberger, P.K.; Elder, N.E.; Grady, C.A.; Gregg, J.L.; Pacheco, C.A.; Greene, C.; Rice, C.; Meyers, T.R. Prevalence of Viral Erythrocytic Necrosis in Pacific Herring and Epizootics in Skagit Bay, Puget Sound, Washington. J. Aquat. Anim. Health 2009, 21, 1–7. [Google Scholar] [CrossRef] [PubMed]
  2. Emmenegger, E.J.; Glenn, J.A.; Winton, J.R.; Batts, W.N.; Gregg, J.L.; Hershberger, P.K. Molecular Identification of Erythrocytic Necrosis Virus (ENV) from the Blood of Pacific Herring (Clupea pallasii). Vet. Microbiol. 2014, 174, 16–26. [Google Scholar] [CrossRef] [PubMed]
  3. Winton, J.R.; Hershberger, P.K. Viral Erythrocytic Necrosis: Chapter 2.2.7. In FHS Blue Book: Suggested Procedures for the Detection and Identification of Certain Finfish and Shellfish Pathogens; American Fisheries Society: Bethesda, MA, USA, 2014. [Google Scholar]
  4. Purcell, M.K.; Pearman-Gillman, S.; Thompson, R.L.; Gregg, J.L.; Hart, L.M.; Winton, J.R.; Emmenegger, E.J.; Hershberger, P.K. Identification of the Major Capsid Protein of Erythrocytic Necrosis Virus (ENV) and Development of Quantitative Real-Time PCR Assays for Quantification of ENV DNA. J. Vet. Diagn. Invest. 2016, 28, 382–391. [Google Scholar] [CrossRef] [PubMed]
  5. Evelyn, T.P.T.; Traxler, G.S. Viral Erythrocytic Necrosis: Natural Occurrence in Pacific Salmon and Experimental Transmission. J. Fish. Res. Board Can. 1978, 35, 903–907. [Google Scholar] [CrossRef]
  6. Rohovec, J.S.; Amandi, A. Incidence of Viral Erythrocytic Necrosis among Hatchery Reared Salmonids of Oregon. Fish Pathol. 1981, 15, 135–141. [Google Scholar] [CrossRef] [Green Version]
  7. Haney, D.C.; Hursh, D.A.; Mix, M.C.; Winton, J.R. Physiological and Hematological Changes in Chum Salmon Artificially Infected with Erythrocytic Necrosis Virus. J. Aquat. Anim. Health 1992, 4, 48–57. [Google Scholar] [CrossRef]
  8. MacMillian, J.R.; Mulcahy, D. Artificial Transmission to and Susceptibility of Puget Sound Fish to Viral Erythrocytic Necrosis (VEN). J. Fish. Res. Board Can. 1979, 36, 1097–1101. [Google Scholar] [CrossRef]
  9. Eaton, W.D. Artificial Transmission of Erythrocytic Necrosis Virus (ENV) from Pacific Herring in Alaska to Chum, Sockeye, and Pink Salmon. J. Appl. Ichthyol. 1990, 6, 136–141. [Google Scholar] [CrossRef]
  10. International Committee on Taxonomy of Viruses; King, A.M.Q.; International Union of Microbiological Societies; Virology Division. Virus Taxonomy: Classification and Nomenclature of Viruses: Ninth Report of the International Committee on Taxonomy of Viruses; Academic Press: London, UK, 2012; ISBN 978-0-12-384684-6. [Google Scholar]
  11. Huang, S.M.; Tu, C.; Tseng, C.H.; Huang, C.C.; Chou, C.-C.; Kuo, H.-C.; Chang, S.-K. Genetic Analysis of Fish Iridoviruses Isolated in Taiwan during 2001–2009. Arch. Virol. 2011, 156, 1505–1515. [Google Scholar] [CrossRef]
  12. Wang, C.S.; Shih, H.H.; Ku, C.C.; Chen, S.N. Studies on Epizootic Iridovirus Infection Among Red Sea Bream, Pagrus major (Temminck & Schlegel), Cultured in Taiwan. J. Fish Dis. 2003, 26, 127–133. [Google Scholar] [PubMed]
  13. Liu, H.I.; Lin, Y.C.; Chiou, P.P.; Chou, H.-Y. Temperature-Dependent Pathogenicity of Grouper Iridovirus of Taiwan (TGIV). J. Mar. Sci. Technol. 2016, 24, 637–644. [Google Scholar]
  14. Wang, Y.Q.; Lü, L.; Weng, S.P.; Huang, J.N.; Chan, S.-M.; He, J.G. Molecular Epidemiology and Phylogenetic Analysis of a Marine Fish Infectious Spleen and Kidney Necrosis Virus-like (ISKNV-like) Virus. Arch. Virol. 2007, 152, 763–773. [Google Scholar] [CrossRef] [PubMed]
  15. Meyers, T.R.; Hauck, A.K.; Blankenbeckler, W.D.; Minicucci, T. First report of Viral Erythrocytic Necrosis in Alaska, USA, Associated with Epizootic Mortality in Pacific Herring, Clupea harengus pallasi (Valenciennes). J. Fish Dis. 1986, 9, 479–491. [Google Scholar] [CrossRef]
  16. Nekouei, O.; Vanderstichel, R.; Ming, T.; Kaukinen, K.H.; Thakur, K.; Tabata, A.; Laurin, E.; Tucker, S.; Beacham, T.D.; Miller, K.M. Detection and Assessment of the Distribution of Infectious Agents in Juvenile Fraser River Sockeye Salmon, Canada, in 2012 and 2013. Front. Microbiol. 2018, 9. [Google Scholar] [CrossRef] [PubMed]
  17. Tucker, S.; Li, S.; Kaukinen, K.H.; Patterson, D.A.; Miller, K.M. Distinct Seasonal Infectious Agent Profiles in Life-history Variants of Juvenile Fraser River Chinook Salmon: An Application of High-Throughput Genomic Screening. PLoS ONE 2018, 13, e0195472. [Google Scholar] [CrossRef] [PubMed]
  18. Thakur, K.K.; Vanderstichel, R.; Li, S.; Laurin, E.; Tucker, S.; Neville, C.; Tabata, A.; Miller, K.M. A Comparison of Infectious Agents Between Hatchery-Enhanced and Wild Out-Migrating Juvenile Chinook Salmon (Oncorhynchus tshawytscha) from Cowichan River, British Columbia. FACETS 2018. [Google Scholar] [CrossRef]
  19. Laurin, E.; Jaramillo, D.; Vanderstichel, R.; Ferguson, H.; Kaukinen, K.H.; Schulze, A.D.; Keith, I.R.; Gardner, I.A.; Miller, K.M. Histopathological and Novel High-Throughput Molecular Monitoring Data from Farmed Salmon (Salmo salar and Oncorhynchus spp.) in British Columbia, Canada, from 2011–2013. Aquaculture 2019, 499, 220–234. [Google Scholar] [CrossRef]
  20. Teffer, A.K.; Bass, A.L.; Miller, K.M.; Patterson, D.A.; Juanes, F.; Hinch, S.G. Infections, Fisheries Capture, Temperature, and Host Responses: Multistressor Influences on Survival and Behaviour of Adult Chinook Salmon. Can. J. Fish. Aquat. Sci. 2018, 75, 2069–2083. [Google Scholar] [CrossRef]
  21. DFO. Stock Assessment Report on Pacific Herring in British Columbia. Canadian Science Advisory Secretariat Science Advisory Report 2009/059. Available online: http://www.dfo-mpo.gc.ca/Library/338390.pdf (accessed on 17 April 2019).
  22. Varpe, Ø.; Fiksen, Ø.; Slotte, A. Meta-Ecosystems and Biological Energy Transport from Ocean to Coast: The Ecological Importance of Herring Migration. Oecologia 2005, 146, 443. [Google Scholar] [CrossRef] [PubMed]
  23. Pacific Salmon Commission. 2010 Annual Report of Catches and Escapements; Report TCChinook (11)-2; Pacific Salmon Commission: Vancouver, BC, Canada, 2011. [Google Scholar]
  24. Government of Canada. Summary of COSEWIC Wildlife Species Assessments. November 2018. Available online: https://www.canada.ca/en/environment-climate-change/services/committee-status-endangered-wildlife/assessments/wildlife-species-assessment-summary-nov-2018.html (accessed on 28 January 2019).
  25. Miller, K.M.; Teffer, A.; Tucker, S.; Li, S.; Schulze, A.D.; Trudel, M.; Juanes, F.; Tabata, A.; Kaukinen, K.H.; Ginther, N.G.; et al. Infectious Disease, Shifting Climates, and Opportunistic Predators: Cumulative Factors Potentially Impacting Wild Salmon Declines. Evol. Appl. 2014, 7, 812–855. [Google Scholar] [CrossRef] [PubMed]
  26. Miller, K.M.; Gardner, I.A.; Vanderstichel, R.; Burnley, T.; Schulze, A.D.; Li, S.; Tabata, A.; Kaukinen, K.H.; Ming, T.J.; Ginther, N.G. Report on the Performance Evaluation of the Fluidigm BioMark Platform for High Throughput Microbe Monitoring in Salmon; Fisheries and Oceans Canada: Ottawa, ON, Canada, 2016.
  27. 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]
  28. Bass, A.L.; Hinch, S.G.; Teffer, A.K.; Patterson, D.A.; Miller, K.M. A Survey of Microparasites Present in Adult Migrating Chinook Salmon (Oncorhynchus tshawytscha) in South-Western British Columbia Determined by High-Throughput Quantitative Polymerase Chain Reaction. J. Fish Dis. 2017, 40, 453–477. [Google Scholar] [CrossRef]
  29. Bolger, A.M.; Lohse, M.; Usadel, B. Trimmomatic: A Flexible Trimmer for Illumina Sequence Data. Bioinforma. Oxf. Engl. 2014, 30, 2114–2120. [Google Scholar] [CrossRef] [PubMed]
  30. Davidson, W.S.; Koop, B.F.; Jones, S.J.; Iturra, P.; Vidal, R.; Maass, A.; Jonassen, I.; Lien, S.; Omholt, S.W. Sequencing the Genome of the Atlantic Salmon (Salmo salar). Genome Biol. 2010, 11, 403. [Google Scholar] [PubMed]
  31. Li, H.; Durbin, R. Fast and Accurate Short Read Alignment with Burrows-Wheeler Transform. Bioinform. Oxf. Engl. 2009, 25, 1754–1760. [Google Scholar] [CrossRef] [PubMed]
  32. Li, H.; Handsaker, B.; Wysoker, A.; Fennell, T.; Ruan, J.; Homer, N.; Marth, G.; Abecasis, G.; Durbin, R. 1000 Genome Project Data Processing Subgroup the Sequence Alignment/Map Format and SAMtools. Bioinforma. Oxf. Engl. 2009, 25, 2078–2079. [Google Scholar] [CrossRef] [PubMed]
  33. Bankevich, A.; Nurk, S.; Antipov, D.; Gurevich, A.A.; Dvorkin, M.; Kulikov, A.S.; Lesin, V.M.; Nikolenko, S.I.; Pham, S.; Prjibelski, A.D.; et al. SPAdes: A New Genome Assembly Algorithm and its Applications to Single-Cell Sequencing. J. Comput. Biol. J. Comput. Mol. Cell Biol. 2012, 19, 455–477. [Google Scholar] [CrossRef] [PubMed]
  34. Buchfink, B.; Xie, C.; Huson, D.H. Fast and Sensitive Protein Alignment using DIAMOND. Nat. Methods 2015, 12, 59–60. [Google Scholar] [CrossRef] [PubMed]
  35. Caporaso, J.G.; Kuczynski, J.; Stombaugh, J.; Bittinger, K.; Bushman, F.D.; Costello, E.K.; Fierer, N.; Peña, A.G.; Goodrich, J.K.; Gordon, J.I.; et al. QIIME Allows Analysis of High-Throughput Community Sequencing Data. Nat. Methods 2010, 7, 335–336. [Google Scholar] [CrossRef]
  36. Borodovsky, M.; Mills, R.; Besemer, J.; Lomsadze, A. Prokaryotic Gene Prediction Using GeneMark and GeneMark.hmm. Curr. Protoc. Bioinform. 2003, 4, 4.5.1–4.5.16. [Google Scholar] [CrossRef] [PubMed]
  37. Altschul, S.F.; Gish, W.; Miller, W.; Myers, E.W.; Lipman, D.J. Basic Local Alignment Search Tool. J. Mol. Biol. 1990, 215, 403–410. [Google Scholar] [CrossRef]
  38. Yutin, N.; Wolf, Y.I.; Raoult, D.; Koonin, E.V. Eukaryotic Large Nucleo-Cytoplasmic DNA Viruses: Clusters of Orthologous Genes and Reconstruction of Viral Genome Evolution. Virol. J. 2009, 6, 223. [Google Scholar] [CrossRef] [PubMed]
  39. Kearse, M.; Moir, R.; Wilson, A.; Stones-Havas, S.; Cheung, M.; Sturrock, S.; Buxton, S.; Cooper, A.; Markowitz, S.; Duran, C.; et al. Geneious Basic: An Integrated and Extendable Desktop Software Platform for the Organization and Analysis of Sequence Data. Bioinformatics 2012, 28, 1647–1649. [Google Scholar] [CrossRef] [PubMed]
  40. Guindon, S.; Dufayard, J.-F.; Lefort, V.; Anisimova, M.; Hordijk, W.; Gascuel, O. New Algorithms and Methods to Estimate Maximum-Likelihood Phylogenies: Assessing the Performance of PhyML 3.0. Syst. Biol. 2010, 59, 307–321. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  41. Thompson, J.D.; Higgins, D.G.; Gibson, T.J. CLUSTAL W: Improving the Sensitivity of Progressive Multiple Sequence Alignment Through Sequence Weighting, Position-Specific Gap Penalties and Weight Matrix Choice. Nucleic Acids Res. 1994, 22, 4673–4680. [Google Scholar] [CrossRef] [PubMed]
  42. R Core Team R: A Language and Environment for Statistical Computing. Ed. Found. Stat. Comput. Vienna Austria 2013.
  43. Pebesma, E.J. Multivariable Geostatistics in S: The gstat Package. Comput. Geosci. 2004, 30, 683–691. [Google Scholar] [CrossRef]
  44. Heuvelink, G.B.M. Spatio-Temporal Interpolation Using gstat by Benedikt Gräler, Edzer Pebesma and Gerard Heuvelink. R J. 2016, 8, 204–218. [Google Scholar]
  45. Kahle, D.; Wickham, H. ggmap: Spatial Visualization with ggplot2. R J. 2013, 5, 144. [Google Scholar] [CrossRef]
  46. Mangiafico, S. rcompanion: Functions to Support Extension Education Program Evaluation [R Statistical Package]. 2016. Available online: https://cran.r-project.org/web/packages/rcompanion/ (accessed on 17 April 2018).
  47. Bates, D.; Mächler, M.; Bolker, B.; Walker, S. Fitting Linear Mixed-Effects Models Using lme4. J. Stat. Softw. 2015, 1. [Google Scholar] [CrossRef]
  48. Hothorn, T.; Bretz, F.; Westfall, P. Simultaneous Inference in General Parametric Models. Biom. J. 2008, 50, 346–363. [Google Scholar] [CrossRef] [Green Version]
  49. Teng, Y.; Hou, Z.; Gong, J.; Liu, H.; Xie, X.; Zhang, L.; Chen, X.; Qin, Q.W. Whole-Genome Transcriptional Profiles of a Novel Marine Fish Iridovirus, Singapore Grouper Iridovirus (SGIV) in Virus-Infected Grouper Spleen Cell Cultures and in Orange-Spotted Grouper, Epinephulus coioides. Virology 2008, 377, 39–48. [Google Scholar] [CrossRef] [PubMed]
  50. Chen, L.M.; Wang, F.; Song, W.; Hew, C.L. Temporal and Differential Gene Expression of Singapore Grouper Iridovirus. J. Gen. Virol. 2006, 87, 2907–2915. [Google Scholar] [CrossRef] [PubMed]
  51. İnce, İ.A.; Özcan, O.; Ilter-Akulke, A.Z.; Scully, E.D.; Özgen, A. Invertebrate Iridoviruses: A Glance Over the Last Decade. Viruses 2018, 10, 161. [Google Scholar] [CrossRef]
  52. Sample, R.; Bryan, L.; Long, S.; Majji, S.; Hoskins, G.; Sinning, A.; Olivier, J.; Chinchar, V.G. Inhibition of Iridovirus Protein Synthesis and Virus Replication by Antisense Morpholino Oligonucleotides Targeted to the Major Capsid Protein, the 18 kDa Immediate-Early Protein, and a Viral Homolog of RNA Polymerase II. Virology 2007, 358, 311–320. [Google Scholar] [CrossRef] [PubMed]
  53. Kurobe, T.; MacConnell, E.; Hudson, C.; McDowell, T.S.; Mardones, F.O.; Hedrick, R.P. Iridovirus Infections among Missouri River Sturgeon: Initial Characterization, Transmission, and Evidence for Establishment of a Carrier State. J. Aquat. Anim. Health 2011, 23, 9–18. [Google Scholar] [CrossRef]
  54. Georgiadis, M.P.; Hedrick, R.P.; Johnson, W.O.; Yun, S.; Gardner, I.A. Risk Factors for Outbreaks of Disease Attributable to White Sturgeon Iridovirus and White Sturgeon Herpesvirus-2 at a Commercial Sturgeon Farm. Am. J. Vet. Res. 2000, 61, 1232–1240. [Google Scholar] [CrossRef]
  55. Traxler, G.; Bell, G. Pathogens Associated with Impounded Pacific Herring Clupea harengus pallasi, with Emphasis on Viral Erythrocytic Necrosis (VEN) and Atypical Aeromonas salmonicida. Dis. Aquat. Organ. 1988, 5, 93–100. [Google Scholar] [CrossRef]
  56. Hershberger, P.; Hart, A.; Gregg, J.; Elder, N.; Winton, J. Dynamics of Viral Hemorrhagic Septicemia, Viral Erythrocytic Necrosis and Ichthyophoniasis in Confined Juvenile Pacific Herring Clupea pallasii. Dis. Aquat. Organ. 2006, 70, 201–208. [Google Scholar] [CrossRef]
  57. Weitkamp, L.A. Marine Distributions of Chinook Salmon from the West Coast of North America Determined by Coded Wire Tag Recoveries. Trans. Am. Fish. Soc. 2010, 139, 147–170. [Google Scholar] [CrossRef]
  58. Myers, K.W.; Klovach, N.V.; Gritsenko, O.F.; Urawa, S.; Royer, T.C. Stock-Specific Distributions of Asian and North American Salmon in the Open Ocean, Interannual Changes, and Oceanographic Conditions; NPAFC: Vancouver, BC, Canada, 2007. [Google Scholar]
  59. Hourston, A.S.; Haegele, C.W. Herring on Canada’s Pacific Coast; Department of Fisheries and Oceans: Longueuil, QC, USA, 1980.
  60. Schroeder, D.C.; Martin, S.J. Deformed wing virus. Virulence 2012, 3, 589–591. [Google Scholar] [CrossRef] [PubMed]
  61. Bell, Q.R.; Traxler, G.S. First Record of Viral Erythrocytic Necrosis and Ceratomyxa Shasta Noble, 1950 (Myxozoa: Myxosporea) in Feral Pink Salmon (Oncorhynchus gorbuscha Walbaum). J. Wildl. Dis. 1985, 21, 169–171. [Google Scholar] [CrossRef] [PubMed]
  62. Peterson, W.; Morgan, C.; Casillas, E.; Peterson, J.; Fisher, J.; W Ferguson, J. Ocean Ecosystem Indicators of Salmon Marine Survival in the Northern California Current. 2019. Available online: https://www.fwspubs.org/doi/suppl/10.3996/042010-JFWM-009/suppl_file/10.3996_042010-jfwm-009.s6.pdf (accessed on 17 April 2019).
  63. Goniea, T.; Bjornn, T.C.; Peery, C.A.; Bennett, D.H.; Stuehrenberg, L.C. Behavioral Thermoregulation and Slowed Migration by Adult Fall Chinook Salmon in Response to High Columbia River Water Temperatures. Trans. Am. Fish. Soc. 2006, 135, 408–419. [Google Scholar] [CrossRef] [Green Version]
  64. Beacham, T.D.; Murray, C.B. Effect of Female Size, Egg Size, and Water Temperature on Developmental Biology of Chum Salmon (Oncorhynchus keta) from the Nitinat River, British Columbia. Can. J. Fish. Aquat. Sci. 1985, 42, 1755–1765. [Google Scholar] [CrossRef]
Viruses 11 00358 g001 550
Figure 1. Substitution-rate optimized maximum likelihood phylogenetic trees of putative ENV DNA-dependent DNA polymerase (A), DNA-dependent RNA polymerase (B), major capsid protein (C), and ATPase (D) sequences (SEQ#) mapped with related iridoviruses. Alignments were based on nucleotide sequences for the major capsid protein and amino acid sequences for all other trees. Branch labels indicate bootstrap support values for 100 re-samplings and the scale bar indicates substitution rate. GenBank reference sequences and contig IDs are listed below (Table 2), with colors indicating genera groups. Putative ENV sequence lengths are listed in Table 1.
Figure 1. Substitution-rate optimized maximum likelihood phylogenetic trees of putative ENV DNA-dependent DNA polymerase (A), DNA-dependent RNA polymerase (B), major capsid protein (C), and ATPase (D) sequences (SEQ#) mapped with related iridoviruses. Alignments were based on nucleotide sequences for the major capsid protein and amino acid sequences for all other trees. Branch labels indicate bootstrap support values for 100 re-samplings and the scale bar indicates substitution rate. GenBank reference sequences and contig IDs are listed below (Table 2), with colors indicating genera groups. Putative ENV sequence lengths are listed in Table 1.
Viruses 11 00358 g001
Viruses 11 00358 g002 550
Figure 2. Apparent ENV prevalence in mixed tissue smolt samples by species. Error bars indicate 95% confidence intervals. Blue bars indicate proportions with limit of detection (LOD) criteria applied and coral bars indicate proportions without LOD criteria applied. Values indicate sample sizes for each species.
Figure 2. Apparent ENV prevalence in mixed tissue smolt samples by species. Error bars indicate 95% confidence intervals. Blue bars indicate proportions with limit of detection (LOD) criteria applied and coral bars indicate proportions without LOD criteria applied. Values indicate sample sizes for each species.
Viruses 11 00358 g002
Viruses 11 00358 g003 550
Figure 3. (A) Full sampling area heat map of the mean calculated ENV copy number, based on interpolated values. Copy number values are binned into 6 numerical categories (see color legend), (B) shows Vancouver Island inset map. Adults and smolts of all species are shown and LOD criteria are not applied.
Figure 3. (A) Full sampling area heat map of the mean calculated ENV copy number, based on interpolated values. Copy number values are binned into 6 numerical categories (see color legend), (B) shows Vancouver Island inset map. Adults and smolts of all species are shown and LOD criteria are not applied.
Viruses 11 00358 g003
Viruses 11 00358 g004 550
Figure 4. Saline (SW) and freshwater (FW) ENV prevalence (A) and load (B) among salmon species and herring. Samples with LOD criteria applied are in blue and samples without LOD criteria applied are in coral. Printed values indicate sample sizes and error bars indicate 95% confidence intervals. Only smolts are shown.
Figure 4. Saline (SW) and freshwater (FW) ENV prevalence (A) and load (B) among salmon species and herring. Samples with LOD criteria applied are in blue and samples without LOD criteria applied are in coral. Printed values indicate sample sizes and error bars indicate 95% confidence intervals. Only smolts are shown.
Viruses 11 00358 g004
Table
Table 1. Basic Local Alignment Search Tool (BLAST) summary of erythrocytic necrosis virus (ENV) sequences from Chinook salmon compared with available herring derived ENV sequences on GenBank. The origin of each sequence is indicated as aquaculture (Aq.) or wild Chinook salmon. Sequences used in phylogenies are in bold.
Table 1. Basic Local Alignment Search Tool (BLAST) summary of erythrocytic necrosis virus (ENV) sequences from Chinook salmon compared with available herring derived ENV sequences on GenBank. The origin of each sequence is indicated as aquaculture (Aq.) or wild Chinook salmon. Sequences used in phylogenies are in bold.
Reference SequenceSequence ID% Identity
(Nucleotide)		% Identity
(Amino Acid)		Sequence Alignment
Length (AA)
ATPase
KJ730210.1 (partial)		80 (Aq.)99.9100849
46 (Wild)99.698.774
57 (Wild)99.598.567
54 (Wild)10010077
DNA-dependent DNA polymerase
KJ756347.1 (partial)		73 (Wild)99.998.359
79 (Aq.)99.899.7897
19 (Wild)99.599.5213
MCP
KT211480.1 (partial, Puget Sound)		86 (Aq.)99.61001436
75 (Wild)100100175
DNA-dependent RNA polymerase
KJ756346.1 (partial)		78 (Aq.)99.996.8972
Table
Table 2. GenBank reference sequences and viral species used to create phylogenies, with colors indicating genera groups corresponding to phylogenetic trees in Figure 1.
Table 2. GenBank reference sequences and viral species used to create phylogenies, with colors indicating genera groups corresponding to phylogenetic trees in Figure 1.
AbbreviationAccession NumbersTaxon
DNA-Dependent DNA PolymeraseDNA-Dependent RNA PolymeraseMCPATPase
NAAAC54632.1YP_762407.1NC_008361.1NC_008361.1Spodoptera ascovirus
ENVAIQ77732.1AIQ77731.1KT211480.1AIN76233.1Erythrocytic necrosis
Megalocytivirus
RSIVBAA28669.1BAK14252.1AB109371.1BAK14298.1Red Sea Bream iridovirus
ISKNVCAZ73994.1AF370008.1NP_612331.1Infectious spleen and kidney necrosis virus
TRBIVADE34365.1ADE34378.1AY590687.2ADE34443.1Turbot reddish body iridovirus
RBIVAAT71848.1Rock Bream iridovirus
Ranavirus
FV3NC_005946.1ASH99239.1AHM26101.1Frog virus 3
EHNVACO25234.1YP_009182042.1AY187045.1YP_009182084.1Epizootic haematopoietic necrosis virus
GIVAY666015.1KX284838.1Grouper iridovirus
ECVAMZ05024.1YP_006347705.1European catfish virus
ATVALN36639.1YP_003852.1Ambystoma tigrinum virus
LMRVKM516719.1Lacerta monticola ranavirus
LMBVKU507317.1Largemouth bass ranavirus
Lymphocystivirus
LCDV(Sa)YP_009342128.1AY823414.1Lymphocystis disease virus (various strains)
LCDV(China)YP_073534.1NC_005902.1YP_073585.1
LCDV1AAX54510.1
Iridovirus
IIV6NC_003038.1AAK82288.1NP_149647.1Invertebrate iridescent virus 6
IIVNC_023615.1Invertebrate iridescent virus
AVIVNC_024451.1NC_024451.1Armadillium vulgare iridescent virus
Aedes taen.NC_008187.1Aedes taeniorhynchus iridescent virus
Chloriridovirus
IIV3YP_654692.1YP_654581.1YP_654693.1Invertebrate iridescent virus 3
RMIVCAC84133.1Regular mosquito iridescent virus

Share and Cite

MDPI and ACS Style

Pagowski, V.A.; Mordecai, G.J.; Miller, K.M.; Schulze, A.D.; Kaukinen, K.H.; Ming, T.J.; Li, S.; Teffer, A.K.; Tabata, A.; Suttle, C.A. Distribution and Phylogeny of Erythrocytic Necrosis Virus (ENV) in Salmon Suggests Marine Origin. Viruses 2019, 11, 358. https://doi.org/10.3390/v11040358

AMA Style

Pagowski VA, Mordecai GJ, Miller KM, Schulze AD, Kaukinen KH, Ming TJ, Li S, Teffer AK, Tabata A, Suttle CA. Distribution and Phylogeny of Erythrocytic Necrosis Virus (ENV) in Salmon Suggests Marine Origin. Viruses. 2019; 11(4):358. https://doi.org/10.3390/v11040358

Chicago/Turabian Style

Pagowski, Veronica A., Gideon J. Mordecai, Kristina M. Miller, Angela D. Schulze, Karia H. Kaukinen, Tobi J. Ming, Shaorong Li, Amy K. Teffer, Amy Tabata, and Curtis A. Suttle. 2019. "Distribution and Phylogeny of Erythrocytic Necrosis Virus (ENV) in Salmon Suggests Marine Origin" Viruses 11, no. 4: 358. https://doi.org/10.3390/v11040358

APA Style

Pagowski, V. A., Mordecai, G. J., Miller, K. M., Schulze, A. D., Kaukinen, K. H., Ming, T. J., Li, S., Teffer, A. K., Tabata, A., & Suttle, C. A. (2019). Distribution and Phylogeny of Erythrocytic Necrosis Virus (ENV) in Salmon Suggests Marine Origin. Viruses, 11(4), 358. https://doi.org/10.3390/v11040358

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