Ultra-deep sequencing of VHSV isolates contributes to understanding the role of viral quasispecies

Positive selection, genetic recombination, and intra-host evolution in novel equine coronavirus genomes and other members of the Embecovirus subgenus

There are several examples of coronaviruses in the Betacoronavirus subgenus Embecovirus that have jumped from an animal to the human host. Studying how evolutionary factors shape coronaviruses in non-human hosts may provide insight into the coronavirus host-switching potential. Equids, such as horses and donkeys, are susceptible to equine coronaviruses (ECoVs). With increased testing prevalence, several ECoV genome sequences have become available for molecular evolutionary analyses, especially those from the United States of America (USA). To date, no analyses have been performed to characterize evolution within coding regions of the ECoV genome. Here, we obtain and describe four new ECoV genome sequences from infected equines from across the USA presenting clinical symptoms of ECoV, and infer ECoV-specific and Embecovirus-wide patterns of molecular evolution. Within two of the four data sets analyzed, we find evidence of intra-host evolution within the nucleocapsid (N) gene, suggestive of quasispecies development. We also identify 12 putative genetic recombination events within the ECoV genome, 11 of which fall in ORF1ab. Finally, we infer and compare sites subject to positive selection on the ancestral branch of each major Embecovirus member clade. Specifically, for the two currently identified human coronavirus (HCoV) embecoviruses that have spilled from animals to humans (HCoV-OC43 and HCoV-HKU1), we find that there are 42 and 2 such sites, respectively, perhaps reflective of the more complex ancestral evolutionary history of HCoV-OC43, which involves several different animal hosts.IMPORTANCEThe Betacoronavirus subgenus Embecovirus contains coronaviruses that not only pose a health threat to animals and humans, but also have jumped from animal to human host. Equids, such as horses and donkeys are susceptible to equine coronavirus (ECoV) infections. No studies have systematically examined evolutionary patterns within ECoV genomes. Our study addresses this gap and provides insight into intra-host ECoV evolution from infected horses. Further, we identify and report natural selection pattern differences between two embecoviruses that have jumped from animals to humans [human coronavirus OC43 and HKU1 (HCoV-OC43 and HCoV-HKU1, respectively)], and hypothesize that the differences observed may be due to the different animal host(s) that each virus circulated in prior to its jump into humans. Finally, we contribute four novel, high-quality ECoV genomes to the scientific community.

Genomic and immunogenic changes of Piscine novirhabdovirus (Viral Hemorrhagic Septicemia Virus) over its evolutionary history in the Laurentian Great Lakes

A unique and highly virulent subgenogroup (-IVb) of Piscine novirhabdovirus, also known as Viral Hemorrhagic Septicemia Virus (VHSV), suddenly appeared in the Laurentian Great Lakes, causing large mortality outbreaks in 2005 and 2006, and affecting >32 freshwater fish species. Periods of apparent dormancy have punctuated smaller and more geographically-restricted outbreaks in 2007, 2008, and 2017. In this study, we conduct the largest whole genome sequencing analysis of VHSV-IVb to date, evaluating its evolutionary changes from 48 isolates in relation to immunogenicity in cell culture. Our investigation compares genomic and genetic variation, selection, and rates of sequence changes in VHSV-IVb, in relation to other VHSV genogroups (VHSV-I, VHSV-II, VHSV-III, and VHSV-IVa) and with other Novirhabdoviruses. Results show that the VHSV-IVb isolates we sequenced contain 253 SNPs (2.3% of the total 11,158 nucleotides) across their entire genomes, with 85 (33.6%) of them being non-synonymous. The most substitutions occurred in the non-coding region (NCDS; 4.3%), followed by the Nv- (3.8%), and M- (2.8%) genes. Proportionally more M-gene substitutions encoded amino acid changes (52.9%), followed by the Nv- (50.0%), G- (48.6%), N- (35.7%) and L- (23.1%) genes. Among VHSV genogroups and subgenogroups, VHSV-IVa from the northeastern Pacific Ocean has shown the fastest substitution rate (2.01x10^-3), followed by VHSV-IVb (6.64x10^-5) and by the VHSV-I, -II and-III genogroups from Europe (4.09x10^-5). A 2016 gizzard shad (Dorosoma cepedianum) from Lake Erie possessed the most divergent VHSV-IVb sequence. The in vitro immunogenicity analysis of that sample displayed reduced virulence (as did the other samples from 2016), in comparison to the original VHSV-IVb isolate (which had been traced back to 2003, as an origin date). The 2016 isolates that we tested induced milder impacts on fish host cell innate antiviral responses, suggesting altered phenotypic effects. In conclusion, our overall findings indicate that VHSV-IVb has undergone continued sequence change and a trend to lower virulence over its evolutionary history (2003 through present-day), which may facilitate its long-term persistence in fish host populations.

Draft genome sequence of scale drop disease virus (SDDV) retrieved from metagenomic investigation of infected barramundi, Lates calcarifer (Bloch, 1790)

Scale drop disease virus (SDDV) is a novel viral pathogen considered to be distributed in farmed barramundi (Lates calcarifer) in South-East Asia. Despite the severity of the disease, only limited genomic information related to SDDV is available. In this study, samples of SDDV-infected fish collected in 2019 were used. The microbiome of brain tissue was investigated using Illumina HiSeq DNA sequencing. Taxonomic analysis showed that SDDV was the main pathogen contained in the affected barramundi. De novo metagenome assembly recovered the SDDV genome, named isolate TH2019, 131 kb in length, and comprised of 135 ORFs. Comparison between this genome and the Singaporean SDDV reference genome revealed that the nucleotide identity within the aligned region was 99.97%. Missense, frameshift, insertion and deletion mutations were identified in 26 ORFs. Deletion of four deduced amino acid sequence in ORF_030L, identical to the SDDV isolate previously identified in Thailand, would be a potential biomarker for future strain classification. Interestingly, the genome of SDDV TH2019 harboured a unique 7,695-bp-long genomic region containing six hypothetical protein-encoded genes. Collectively, this study demonstrated that the SDDV genome can be sequenced directly, although with limited coverage depth, using metagenomic analysis of barramundi sample with severe infection.

Evolutionary trajectory of fish Piscine novirhabdovirus (=Viral Hemorrhagic Septicemia Virus) across its Laurentian Great Lakes history: Spatial and temporal diversification

Piscine novirhabdovirus = Viral Hemorrhagic Septicemia Virus (VHSV) first appeared in the Laurentian Great Lakes with large outbreaks from 2005 to 2006, as a new and novel RNA rhabdovirus subgenogroup (IVb) that killed >30 fish species. Interlude periods punctuated smaller more localized outbreaks in 2007, 2010, and 2017, although some fishes tested positive in the intervals. There have not been reports of outbreaks or positives from 2018, 2019, or 2020. Here, we employ a combined population genetics and phylogenetic approach to evaluate spatial and temporal evolutionary trajectory on its G-gene sequence variation, in comparison with whole-genome sequences (11,083 bp) from a subset of 44 individual isolates (including 40 newly sequenced ones). Our results show that IVb (N = 184 individual fish isolates) diversified into 36 G-gene haplotypes from 2003 to 2017, stemming from two originals ("a" and "b"). G-gene haplotypes "a" and "b" differed by just one synonymous single-nucleotide polymorphism (SNP) substitution, remained the most abundant until 2011, then disappeared. Group "a" descendants (14 haplotypes) remained most prevalent in the Upper and Central Great Lakes, with eight (51%) having nonsynonymous substitutions. Group "b" descendants primarily have occurred in the Lower Great Lakes, including 22 haplotypes, of which 15 (68%) contained nonsynonymous changes. Evolutionary patterns of the whole-genome sequences (which had 34 haplotypes among 44 isolates) appear congruent with those from the G-gene. Virus populations significantly diverged among the Upper, Central, and Lower Great Lakes, diversifying over time. Spatial divergence was apparent in the overall patterns of nucleotide substitutions, while amino acid changes increased temporally. VHSV-IVb thus significantly differentiated across its less than two decades in the Great Lakes, accompanied by declining outbreaks and virulence. Continuing diversification likely allowed the virus to persist at low levels in resident fish populations, and may facilitate its potential for further and future spread to new habitats and nonacclimated hosts.

Sequencing of animal viruses: quality data assurance for NGS bioinformatics

BACKGROUND

Next generation sequencing (NGS) is becoming widely used among diagnostics and research laboratories, and nowadays it is applied to a variety of disciplines, including veterinary virology. The NGS workflow comprises several steps, namely sample processing, library preparation, sequencing and primary/secondary/tertiary bioinformatics (BI) analyses. The latter is constituted by a complex process extremely difficult to standardize, due to the variety of tools and metrics available. Thus, it is of the utmost importance to assess the comparability of results obtained through different methods and in different laboratories. To achieve this goal, we have organized a proficiency test focused on the bioinformatics components for the generation of complete genome sequences of salmonid rhabdoviruses.

METHODS

Three partners, that performed virus sequencing using different commercial library preparation kits and NGS platforms, gathered together and shared with each other 75 raw datasets which were analyzed separately by the participants to produce a consensus sequence according to their own bioinformatics pipeline. Results were then compared to highlight discrepancies, and a subset of inconsistencies were investigated more in detail.

RESULTS

In total, we observed 526 discrepancies, of which 39.5% were located at genome termini, 14.1% at intergenic regions and 46.4% at coding regions. Among these, 10 SNPs and 99 indels caused changes in the protein products. Overall reproducibility was 99.94%. Based on the analysis of a subset of inconsistencies investigated more in-depth, manual curation appeared the most critical step affecting sequence comparability, suggesting that the harmonization of this phase is crucial to obtain comparable results. The analysis of a calibrator sample allowed assessing BI accuracy, being 99.983%.

CONCLUSIONS

We demonstrated the applicability and the usefulness of BI proficiency testing to assure the quality of NGS data, and recommend a wider implementation of such exercises to guarantee sequence data uniformity among different virology laboratories.

Genome and Phylogenetic Analysis of Infectious Hematopoietic Necrosis Virus Strain SNU1 Isolated in Korea

Infectious hematopoietic necrosis virus (IHNV), one of the most important pathogenic fish viruses, affects trout fisheries and causes considerable economic losses. Currently, in Korea, more studies on IHNV infection are being reported. However, relatively less data is available on Korean isolates than on those from other countries. Few studies have focused on gene sequence analyses of IHNV glycoprotein (G) gene and almost none have focused on other gene fragments. Therefore, considering the dearth of adequate phylogenetic and genomic studies on Korean IHNV strains because of the lack of data, our study aimed to provide sufficient relevant data by sequencing the complete genome of the IHNV strain SNU1, which was recently isolated from a Korean rainbow trout farm. Moreover, we focused on expanding the perspectives on the phylogenesis of IHNV isolates from Korea and other Asian countries. IHNV was isolated from pooled hematopoietic tissue samples using Epithelioma papulosum cyprinid (EPC) cells, and phylogenetic analysis and genome study were conducted using complete G, N, and nonvirion (NV) gene sequences. Our main achievements were the development of a phylogenetic analytical method based on the NV gene and complete genome sequence analysis of the IHNV strain SNU1, which was compared with other Asian isolate sequences.

Investigation of round goby viral haemorrhagic septicaemia outbreak in New York

Eleven viral haemorrhagic septicaemia virus (VHSV) genotype IVb isolates were sequenced, and their genetic variation explored to determine the source of a VHS outbreak on the eastern shore of Cayuga Lake. An active fish kill of round gobies (Neogobius melanostomus, Pallas) was intensively sampled at King Ferry, NY and nearby Long Point State Park in May 2017. Gross lesions observed on 67 moribund round gobies and two rock bass (Ambloplites rupestris, Rafinesque) included moderately haemorrhagic internal organs and erythematous areas on the head, flank, and fins. RT-qPCR tests for VHSV were positive for all 69 fish. Viral isolation on epithelioma papulosum cyprinid cells showed cytopathic effect characteristic of VHSV for six round goby samples from King Ferry. The complete nucleotide sequence of the VHSV IVb genomes of five Cayuga Lake round goby isolates were derived on an Illumina platform along with 2017 VHSV IVb isolates from round gobies collected from the following: Lake Erie near Dunkirk, NY; the St. Lawrence River near Clayton and Cape Vincent, NY; and Lake St. Lawrence near Massena, NY. The phylogenetic tree created from these aligned sequences and four other complete VHSV IVb genomes shows Cayuga Lake isolates are closely related to the Lake Erie isolates.

DNA vaccination for finfish aquaculture

In fish, DNA vaccines have been shown to give very high protection in experimental facilities against a number of viral diseases, particularly diseases caused by rhabdoviruses. However, their efficacy in generating protection against other families of fish viral pathogens is less clear. One DNA vaccine is currently in use commercially in fish farms in Canada and the commercialisation of another was authorised in Europe in 2017. The mechanism of action of DNA vaccines, including the role of the innate immune responses induced shortly after DNA vaccination in the activation of the adaptive immunity providing longer term specific protection, is still not fully understood. In Europe the procedure for the commercialisation of a veterinary DNA vaccine requires the resolution of certain concerns particularly about safety for the host vaccinated fish, the consumer and the environment. Relating to consumer acceptance and particularly environmental safety, a key question is whether a DNA vaccinated fish is considered a Genetically Modified Organism (GMO). In the present opinion paper these key aspects relating to the mechanisms of action, and to the development and the use of DNA vaccines in farmed fish are reviewed and discussed.

Complete sequences of 4 viral hemorrhagic septicemia virus IVb isolates and their virulence in northern pike fry

Four viral hemorrhagic septicemia virus (VHSV) genotype IVb isolates were sequenced, their genetic variation explored, and comparative virulence assayed with experimental infections of northern pike Esox lucius fry. In addition to the type strain MI03, the complete 11183 bp genome of the first round goby Neogobius melanostomus isolate from the St. Lawrence River, and the 2013 and 2014 isolates from gizzard shad Dorosoma cepedianum die-offs in Irondequoit Bay, Lake Ontario and Dunkirk Harbor, Lake Erie were all deep sequenced on an Illumina platform. Mutations documented in the 11 yr since the MI03 index case from Lake St. Clair muskellunge Esox masquinongy showed 87 polymorphisms among the 4 isolates. Twenty-six mutations were non-synonymous and located at 18 different positions within the matrix protein, glycoprotein, non-virion protein, and RNA polymerase genes. The same 4 isolates were used to infect northern pike fry by a single 1 h bath exposure. Cumulative percent mortality varied from 42.5 to 62.5%. VHSV was detected in 57% (41/72) of the survivors at the end of the 21-d trial, suggesting that the virus was not rapidly cleared. Lesions were observed in many of the moribund and dead northern pike, such as hemorrhaging in the skin and fins, as well as hydrocephalus. Mean viral load measured from the trunk and visceral tissues of MI03-infected pike was significantly higher than the quantities detected in fish infected with the most recent isolates of genotype IVb, but there were no differences in cumulative mortality observed.

Contribution of Next-Generation Sequencing to Aquatic and Fish Virology

The recent technological advances in nucleic acid sequencing, called next-generation sequencing (NGS), have revolutionized the field of genomics and have also influenced viral research. Aquatic viruses, and especially those infecting fish, have also greatly benefited from NGS technologies, which provide a huge amount of molecular information at a low cost in a relatively short period of time. Here, we review the use of the current high-throughput sequencing platforms with a special focus on the associated challenges (regarding sample preparation and bioinformatics) in their applications to the field of aquatic virology, especially for: (i) discovering novel viruses that may be associated with fish mortalities, (ii) elucidating the mechanisms of pathogenesis, and finally (iii) studying the molecular epidemiology of these pathogens.

The Promise of Whole Genome Pathogen Sequencing for the Molecular Epidemiology of Emerging Aquaculture Pathogens

Aquaculture is the fastest growing food-producing sector, and the sustainability of this industry is critical both for global food security and economic welfare. The management of infectious disease represents a key challenge. Here, we discuss the opportunities afforded by whole genome sequencing of bacterial and viral pathogens of aquaculture to mitigate disease emergence and spread. We outline, by way of comparison, how sequencing technology is transforming the molecular epidemiology of pathogens of public health importance, emphasizing the importance of community-oriented databases and analysis tools.

A bioinformatics approach reveals seven nearly-complete RNA-virus genomes in bivalve RNA-seq data

Viral metagenomics (viromics) can provide a great contribution in expanding the knowledge of viruses and the relationship with their hosts. Viromic studies on marine organisms are still at a very early stage and only little efforts have been spent in the identification of viruses associated to marine invertebrates to date, leaving the complexity of marine viromes associated to bivalve hosts almost completely unexplored. However, the potential use of viromic approaches in the management of viral diseases affecting aquacultured species has been recently evidenced by the flourishing of studies on the Ostreid herpesvirus type-1, which has been associated with bivalve mortality events. Herein we discuss an effective pipeline to retrieve and reconstruct nearly complete and previously unreported viral genomes from existing host RNA-seq data. As a case study, we report the identification of seven RNA-virus genomes within the frame of a highly diversified viral community that characterizes both Crassostrea gigas and Mytilus galloprovincialis samples collected from the lagoon of Goro (Italy).

Can VHS Virus Bypass the Protective Immunity Induced by DNA Vaccination in Rainbow Trout?

DNA vaccines encoding viral glycoproteins have been very successful for induction of protective immunity against diseases caused by rhabdoviruses in cultured fish species. However, the vaccine concept is based on a single viral gene and since RNA viruses are known to possess high variability and adaptation capacity, this work aimed at evaluating whether viral haemorrhagic septicaemia virus (VHSV), an RNA virus and member of Rhabdoviridae family, was able to evade the protective immune response induced by the DNA vaccination of rainbow trout. The experiments comprised repeated passages of a highly pathogenic VHSV isolate in a fish cell line in the presence of neutralizing fish serum (in vitro approach), and in rainbow trout immunized with the VHS DNA vaccine (in vivo approach). For the in vitro approach, the virus collected from the last passage (passaged virus) was as sensitive as the parental virus to serum neutralization, suggesting that the passaging did not promote the selection of virus populations able to bypass the neutralization by serum antibodies. Also, in the in vivo approach, where virus was passaged several times in vaccinated fish, no increased virulence nor increased persistence in vaccinated fish was observed in comparison with the parental virus. However, some of the vaccinated fish did get infected and could transmit the infection to naïve cohabitant fish. The results demonstrated that the DNA vaccine induced a robust protection, but also that the immunity was non-sterile. It is consequently important not to consider vaccinated fish as virus free in veterinary terms.