Gene Diversification of an Emerging Pathogen: A Decade of Mutation in a Novel Fish Viral Hemorrhagic Septicemia (VHS) Substrain since Its First Appearance in the Laurentian Great Lakes

Transcriptome Profiling of Oncorhynchus mykiss Infected with Low or Highly Pathogenic Viral Hemorrhagic Septicemia Virus (VHSV)

The rainbow trout (Oncorhynchus mykiss) is the most important produced species in freshwater within the European Union, usually reared in intensive farming systems. This species is highly susceptible to viral hemorrhagic septicemia (VHS), a severe systemic disease widespread globally throughout the world. Viral hemorrhagic septicemia virus (VHSV) is the etiological agent and, recently, three classes of VHSV virulence (high, moderate, and low) have been proposed based on the mortality rates, which are strictly dependent on the viral strain. The molecular mechanisms that regulate VHSV virulence and the stimulated gene responses in the host during infection are not completely unveiled. While some preliminary transcriptomic studies have been reported in other fish species, to date there are no publications on rainbow trout. Herein, we report the first time-course RNA sequencing analysis on rainbow trout juveniles experimentally infected with high and low VHSV pathogenic Italian strains. Transcriptome analysis was performed on head kidney samples collected at different time points (1, 2, and 5 days post infection). A large set of notable genes were found to be differentially expressed (DEGs) in all the challenged groups (e.s. trim63a, acod1, cox-2, skia, hipk1, cx35.4, ins, mtnr1a, tlr3, tlr7, mda5, lgp2). Moreover, the number of DEGs progressively increased especially during time with a greater amount found in the group infected with the high VHSV virulent strain. The gene ontology (GO) enrichment analysis highlighted that functions related to inflammation were modulated in rainbow trout during the first days of VHSV infection, regardless of the pathogenicity of the strain. While some functions showed slight differences in enrichments between the two infected groups, others appeared more exclusively modulated in the group challenged with the highly pathogenic strain.

Population genetic analyses unveiled genetic stratification and differential natural selection signatures across the G-gene of viral hemorrhagic septicemia virus

Introduction: Viral hemorrhagic septicemia virus (VHSV) is the most lethal pathogen in aquaculture, infecting more than 140 fish species in marine, estuarine, and freshwater environments. Viral hemorrhagic septicemia virus is an enveloped RNA virus that belongs to the family Rhabdoviridae and the genus Novirhabdovirus. The current study is designed to infer the worldwide Viral hemorrhagic septicemia virus isolates' genetic diversity and evolutionary dynamics based on G-gene sequences. Methods: The complete G-gene sequences of viral hemorrhagic septicemia virus were retrieved from the public repositories with known timing and geography details. Pairwise statistical analysis was performed using Arlequin. The Bayesian model-based approach implemented in STRUCTURE software was used to investigate the population genetic structure, and the phylogenetic tree was constructed using MEGA X and IQ-TREE. The natural selection analysis was assessed using different statistical approaches, including IFEL, MEME, and SLAC. Results and Discussion: The global Viral hemorrhagic septicemia virus samples are stratified into five genetically distinct subpopulations. The STRUCTURE analysis unveiled spatial clustering of genotype Ia into two distinct clusters at K = 3. However, at K = 5, the genotype Ia samples, deposited from Denmark, showed temporal distribution into two groups. The analyses unveiled that the genotype Ia samples stratified into subpopulations possibly based on spatiotemporal distribution. Several viral hemorrhagic septicemia virus samples are characterized as genetically admixed or recombinant. In addition, differential or subpopulation cluster-specific natural selection signatures were identified across the G-gene codon sites among the viral hemorrhagic septicemia virus isolates. Evidence of low recombination events elucidates that genetic mutations and positive selection events have possibly driven the observed genetic stratification of viral hemorrhagic septicemia virus samples.

Current status and advances of fish vaccines in Malaysia

Fish diseases have a significant negative influence on the Malaysian aquaculture industry. Since the 1980s, the sector has grown in size, which has resulted in a rise in the prevalence of infectious outbreaks affecting both freshwater and marine cultured fish species. Demand for commercially available fish vaccinations is predicted to increase as infectious disease outbreaks continue to occur. In Malaysia, aquaculture vaccine research and development (R&D) are still in its infancy, with most efforts concentrating on producing vaccines against bacterial infections, most notably streptococcosis, vibriosis, and motile Aeromonas septicemia. Despite several attempts, no homegrown vaccine has been effectively introduced into the manufacturing pipeline to date. At the moment, only three imported aquatic vaccines have received full permission, a far cry from the 314 and 60 vaccines licensed in the poultry and porcine industries, respectively. This review will describe recent findings regarding the development of aquaculture vaccines for certain fish species and diseases in Malaysia. In our opinion, R&D on fish vaccines is critical to the aquaculture industry's viability.

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.

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.

Comparative effects of Novirhabdovirus genes on modulating constitutive transcription and innate antiviral responses, in different teleost host cell types

BACKGROUND

Infectious hematopoietic necrosis virus (IHNV) and viral hemorrhagic septicemia virus (VHSV) are highly contagious, pathogenic Novirhabdoviruses affecting fish and are thusly notifiable diseases with the World Organization for Animal Health. This study assessed the relative capacities of IHNV and VHSV genes to modulate host general transcription and explores the abilities of specific IHNV genes to interfere with the interferon pathway in heterogenous teleost cell-lines.

METHODS

Optimized protocols allowed for efficient transient transfections in EPC, BF-2, RTG-2 and RTgill-W1 cell lines of plasmids encoding IHNV (M genogroup) and VHSV (-IVb genotype) genes, including N, P, M, G and NV. Their impact on general cellular transcription was measured 48 hours post transfection (hpt) with luciferase constructs driven by a modified β-Actin promoter (pCAG). Their modulation of the innate antiviral immune response was characterized 72 hpt, using luciferase constructs measuring rainbow trout Type I IFN or MX-1 promoter augmentation, upon MAVS co-transfection.

RESULTS

M was generally confirmed as the strongest constitutive transcriptional suppressor while IHNV P, but not VHSV P, augmented constitutive transcription in fibroblastic cell types. Cell-specific effects were observed for viral G gene, with VHSV G exhibiting suppression of basal transcription in EPC and BF-2 but not in trout cells; while IHNV G was stimulatory in RTG-2, but inhibitory in RTgill-W1. NV consistently stimulated constitutive transcription, with higher augmentation patterns seen in fibroblastic compared to epithelial cells, and for IHNV NV compared to VHSV NV. The innate antiviral immune response, focusing on the IFN pathway, was silenced by IHNV M in all cell lines tested. IHNV N showed a dose-dependent suppression of type I IFN, but with minor effects on MX-1. IHNV P and G played minor IFN-inhibitory roles, consistent and dose-dependent only for G in rainbow trout cells. IHNV NV mediated a consistent stimulatory effect on either Type I IFN or MX-1, but much less pronounced in RTgill-W1.

CONCLUSIONS

This study extends our understanding of Novirhabdoviruses-host interaction, showing differential innate immune responses in heterogenous cell types. Viral regulators of innate immune signaling are identified, either as dose-dependent suppressors (such as M and N) or stimulators (mainly NV), indicating novel targets for the design of more efficient vaccination strategies.

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.

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.

Identification of the functional regions of the viral haemorrhagic septicaemia virus (VHSV) NV protein: Variants that improve function

Non-virion (NV) protein is essential for an efficient replication increasing the pathogenicity of the Salmonid novirhabdovirus (formerly IHNV), Piscine novirhabdovirus (formerly VHSV), and Hirame novirhabdovirus (HIRV). The interferon system, apoptosis, and other immune-related genes are modulated by NV to finally induce a deficient antiviral state in the cell. However, little is known about the VHSV NV regions involved in function and location. Here, eight different NV 07.71 fragments and eleven NV 07.71 mutants derived from the region between the two first α-helices have been studied in order to establish the mx and il8 transcript levels in ZF4 cells and the subcellular location. As a result, we determined that the N-terminal part of NV preserves the same ability as the wild-type (wt) NV in mx/il8 modulation and it also shares the subcellular location. Among NV mutants, some induced mx upregulation (N34A, C35A, D38A, and S40A) but maintained the il8 levels stable when compared to wt-NV in ZF4. Four NV mutants (D28A, N31A, L33A, and F37A) were not affected by the mutation and showed mx and il8 transcript levels similar to wt-NV. Surprisingly, mutants D36A, R39A, and D41A induced a stronger downregulation of both mx and il8 transcript levels than wt-NV, suggesting that a more stable structure and an improved interaction with ligands could be achieved through these mutations. Amino acids at positions 36 and 39 are conserved among known VHSV NV proteins whereas at position 41 two different amino acids have been described. To date, no natural NV proteins with alanine at positions 36, 39, and 41 have been found. In addition, wt-NV, all NV mutants, and one N-terminal NV fragment were located at cytoplasm with a characteristic pattern, which might support that cytoplasm is the site for interaction with candidate ligands such as PPM1Bb. Taken together, the data presented in this work indicated that NV function relies on the first part of the molecule and is dependent on tertiary structure rather than on the linear one. This study could lead to a better knowledge of VHSV escape from fish antiviral mechanisms as well as to future studies on immune targets.

Role of Viral Hemorrhagic Septicemia Virus Matrix (M) Protein in Suppressing Host Transcription

Viral hemorrhagic septicemia virus (VHSV) is a pathogenic fish rhabdovirus found in discrete locales throughout the Northern Hemisphere. VHSV infection of fish cells leads to upregulation of the host's virus detection response, but the virus quickly suppresses interferon (IFN) production and antiviral gene expression. By systematically screening each of the six VHSV structural and nonstructural genes, we identified matrix protein (M) as the virus' most potent antihost protein. Only M of VHSV genotype IV sublineage b (VHSV-IVb) suppressed mitochondrial antiviral signaling protein (MAVS) and type I IFN-induced gene expression in a dose-dependent manner. M also suppressed the constitutively active simian virus 40 (SV40) promoter and globally decreased cellular RNA levels. Chromatin immunoprecipitation (ChIP) studies illustrated that M inhibited RNA polymerase II (RNAP II) recruitment to gene promoters and decreased RNAP II C-terminal domain (CTD) Ser2 phosphorylation during VHSV infection. However, transcription directed by RNAP I to III was suppressed by M. To identify regions of functional importance, M proteins from a variety of VHSV strains were tested in cell-based transcriptional inhibition assays. M of a particular VHSV-Ia strain, F1, was significantly less potent than IVb M at inhibiting SV40/luciferase (Luc) expression yet differed by just 4 amino acids. Mutation of D62 to alanine alone, or in combination with an E181-to-alanine mutation (D62A E181A), dramatically reduced the ability of IVb M to suppress host transcription. Introducing either M D62A or D62A E181A mutations into VHSV-IVb via reverse genetics resulted in viruses that replicated efficiently but exhibited less cytotoxicity and reduced antitranscriptional activities, implicating M as a primary regulator of cytopathicity and host transcriptional suppression.IMPORTANCE Viruses must suppress host antiviral responses to replicate and spread between hosts. In these studies, we identified the matrix protein of the deadly fish novirhabdovirus VHSV as a critical mediator of host suppression during infection. Our studies indicated that M alone could block cellular gene expression at very low expression levels. We identified several subtle mutations in M that were less potent at suppressing host transcription. When these mutations were engineered back into recombinant viruses, the resulting viruses replicated well but elicited less toxicity in infected cells and activated host innate immune responses more robustly. These data demonstrated that VHSV M plays an important role in mediating both virus-induced cell toxicity and viral replication. Our data suggest that its roles in these two processes can be separated to design effective attenuated viruses for vaccine candidates.

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.