Investigation of round goby viral haemorrhagic septicaemia outbreak in New York

Salinity tolerance of Round Goby: Informing invasion potential in North American coastal watersheds

Since being introduced into the Laurentian Great Lakes in the 1990s, round goby (Neogobius melanostomus) has spread rapidly, reaching the Hudson River Estuary, NY in 2021. To address the expansion potential into saline environments from this North American coastal invasion front, we experimentally assessed the salinity tolerance of adult round gobies. Water temperatures vary widely in temperate aquatic ecosystems, and our study is the first to investigate the effect of temperature on the salinity tolerance of round goby, conducting tolerance trials at three temperatures: a preferred temperature of 20°C, 26°C reflective of summer conditions, and 5°C reflective of winter conditions. Adult round gobies were subjected to weekly salinity increases of 3 parts per thousand (ppt), concluding at 33ppt. Study specimens were monitored for stress cues (behavior changes and color changes), and mortality. We found significant salinity tolerance differences dependent on water temperature, with the highest tolerance at 5°C and the lowest tolerance at 26°C. By 30ppt, survival was 87% at 5°C and only 7% at 26°C. Based on mortality results, round goby expansion may occur year-round into brackish portions (<21ppt) of the Hudson River Estuary as far south as the Harlem River. This would open access to portions of Long Island Sound, potentially rendering other coastal watersheds to be at risk of invasion. However, temperature-dependent salinity tolerance findings suggest round goby expansion potential into high salinity habitats may be seasonally dependent, with expansion opportunities occurring in colder months and expansion barriers occurring in warmer months. To assess longer term survival and body condition, another experiment maintained round gobies at four sustained salinities (≤ 21ppt) for ten weeks at 20°C. Growth and hepatosomatic index at 21ppt were significantly lower (p<0.001) than at 1, 9, and 15ppt, indicating sustained exposure to higher salinities may affect energy stores, potentially limiting establishment potential.

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.