Comparative Genomics of an Emerging Amphibian Virus

Age- and dose-dependent susceptibility of axolotls (Ambystoma mexicanum) by bath exposure to Ambystoma tigrinum virus (ATV)

Ranaviruses are large, dsDNA viruses that have significant ecological and economic impact on cold-blooded vertebrates. However, our understanding of the viral proteins and subsequent host immune response(s) that impact susceptibility to infection and disease is not clear. The ranavirus Ambystoma tigrinum virus (ATV), originally isolated from the Sonoran tiger salamander (Ambystoma mavortium stebbinsi), is highly pathogenic at low doses of ATV at all tiger salamander life stages and this model has been used to explore the host-pathogen interactions of ATV infection. However, inconsistencies in the availability of laboratory reared larval tiger salamanders required us to look at the well characterized axolotl (A. mexicanum) as a model for ATV infection. Data obtained from five infection experiments over different developmental timepoints suggest that axolotls are susceptible to ATV in an age- and dose-dependent manner. These data support the use of the ATV-axolotl model to further explore the host-pathogen interactions of ranavirus infections.

Interspecific introgression of MHC genes in Triturus newts: Evidence from multiple contact zones

The major histocompatibility complex (MHC) genes are central to the adaptive immune response in vertebrates. Selection generally maintains high MHC variation because the spectrum of recognized pathogens depends on MHC polymorphism. Novel alleles favoured by selection originate by interallelic recombination or de novo mutations but may also be acquired by introgression from related species. However, the extent and prevalence of MHC introgression remain an open question. In this study, we tested for MHC introgression in six hybrid zones formed by six Triturus newt species. We sequenced and genotyped the polymorphic second exons of the MHC class I and II genes and compared their interspecific similarity at various distances from the centre of the hybrid zone. We found evidence for introgression of both MHC classes in the majority of examined hybrid zones, with support for a more substantial class I introgression. Furthermore, the overall MHC allele sharing outside of hybrid zones was elevated between pairs of Triturus species with abutting ranges, regardless of the phylogenetic distance between them. No effect of past hybrid zone movement on MHC allele sharing was found. Finally, using previously published genome-wide data, we demonstrated that MHC introgression was more extensive than genome-wide introgression, supporting its adaptive potential. Our study thus provides evidence for the prevalence of MHC introgression across multiple Triturus hybrid zones, indicating that MHC introgression between divergent hybridizing species may be widespread and adaptive.

High prevalence of subclinical frog virus 3 infection in freshwater turtles of Ontario, Canada

Ranaviruses have been associated with chelonian mortality. In Canada, the first two cases of ranavirus were detected in turtles in 2018 in Ontario, although a subsequent survey of its prevalence failed to detect additional positive cases. To confirm the prevalence of ranavirus in turtles in Ontario, we used a more sensitive method to investigate if lower level persistent infection was present in the population. Here we report results via a combination of qPCR, PCR, Sanger sequencing and genome sequencing from turtles from across Ontario, with no clinical signs of illness. We found 2 positives with high viral load and 5 positives with low viral load. Histopathology found subtle histological changes. DNA sequences identified two types of frog virus 3 (FV3), and genome sequencing identified a ranavirus similar to wild-type FV3. Our results show that the virus has been present in Ontario's turtles as subclinical infections.

Genome analysis of Ranavirus frog virus 3 isolated from American Bullfrog (Lithobates catesbeianus) in South America

Ranaviruses (family Iridoviridae) cause important diseases in cold-blooded vertebrates. In addition, some occurrences indicate that, in this genus, the same virus can infect animals from different taxonomic groups. A strain isolated from a Ranavirus outbreak (2012) in the state of Sao Paulo, Brazil, had its genome sequenced and presented 99.26% and 36.85% identity with samples of Frog virus 3 (FV3) and Singapore grouper iridovirus (SGIV) ranaviruses, respectively. Eight potential recombination events among the analyzed sample and reference FV3 samples were identified, including a recombination with Bohle iridovirus (BIV) sample from Oceania. The analyzed sample presented several rearrangements compared to FV3 reference samples from North America and European continent. We report for the first time the complete genome of Ranavirus FV3 isolated from South America, these results contribute to a greater knowledge related to evolutionary events of potentially lethal infectious agent for cold-blooded animals.

Frog Virus 3 Genomes Reveal Prevalent Recombination between Ranavirus Lineages and Their Origins in Canada

Ranaviruses are pathogens associated with the decline of amphibian populations across much of their distribution. In North America, frog virus 3 (FV3) is a widely distributed pathogen with wild populations of amphibians harboring different lineages and putative recombinants between FV3 and common midwife toad virus (CMTV). These recombinants have higher pathogenicity, and CMTV-derived genes associated with virulence are reported in wild strains in Canada. However, while FV3 is linked to amphibian die-offs in North America, CMTVs have been reported only in commercial frog farms in North America. We sequenced complete genomes of 18 FV3 isolates from three amphibian species to characterize genetic diversity of the lineages in Canada and infer possible recombinant regions. The 18 FV3 isolates displayed different signals of recombination, varying from none to interspersed recombination with previously isolated CMTV-like viruses. In general, most recombination breakpoints were located within open reading frames (ORFs), generating new ORFs and proteins that were a mixture between FV3 and CMTV. A combined spatial and temporal phylogeny suggests the presence of the FV3 lineage in Canada is relatively contemporary (<100 years), corroborating the hypothesis that both CMTV- and FV3-like viruses spread to North America when the international commercial amphibian trade started. Our results highlight the importance of pathogen surveillance and viral dynamics using full genomes to more clearly understand the mechanisms of disease origin and spread.IMPORTANCE Amphibian populations are declining worldwide, and these declines have been linked to a number of anthropogenic factors, including disease. Among the pathogens associated with amphibian mortality, ranaviruses have caused massive die-offs across continents. In North America, frog virus 3 (FV3) is a widespread ranavirus that can infect wild and captive amphibians. In this study, we sequenced full FV3 genomes isolated from frogs in Canada. We report widespread recombination between FV3 and common midwife toad virus (CMTV). Phylogenies indicate a recent origin for FV3 in Canada, possibly as a result of international amphibian trade.

Pathogen Risk Analysis for Wild Amphibian Populations Following the First Report of a Ranavirus Outbreak in Farmed American Bullfrogs (Lithobates catesbeianus) from Northern Mexico

Ranaviruses are the second deadliest pathogens for amphibian populations throughout the world. Despite their wide distribution in America, these viruses have never been reported in Mexico, the country with the fifth highest amphibian diversity in the world. This paper is the first to address an outbreak of ranavirus in captive American bullfrogs (Lithobates catesbeianus) from Sinaloa, Mexico. The farm experienced high mortality in an undetermined number of juveniles and sub-adult bullfrogs. Affected animals displayed clinical signs and gross lesions such as lethargy, edema, skin ulcers, and hemorrhages consistent with ranavirus infection. The main microscopic lesions included mild renal tubular necrosis and moderate congestion in several organs. Immunohistochemical analyses revealed scant infected hepatocytes and renal tubular epithelial cells. Phylogenetic analysis of five partial ranavirus genes showed that the causative agent clustered within the Frog virus 3 clade. Risk assessment with the Pandora⁺ protocol demonstrated a high risk for the pathogen to affect amphibians from neighboring regions (overall Pandora risk score: 0.619). Given the risk of American bullfrogs escaping and spreading the disease to wild amphibians, efforts should focus on implementing effective containment strategies and surveillance programs for ranavirus at facilities undertaking intensive farming of amphibians.

Ranavirus genotypes in the Netherlands and their potential association with virulence in water frogs (Pelophylax spp.)

Ranaviruses are pathogenic viruses for poikilothermic vertebrates worldwide. The identification of a common midwife toad virus (CMTV) associated with massive die-offs in water frogs (Pelophylax spp.) in the Netherlands has increased awareness for emerging viruses in amphibians in the country. Complete genome sequencing of 13 ranavirus isolates collected from ten different sites in the period 2011–2016 revealed three CMTV groups present in distinct geographical areas in the Netherlands. Phylogenetic analysis showed that emerging viruses from the northern part of the Netherlands belonged to CMTV-NL group I. Group II and III viruses were derived from the animals located in the center-east and south of the country, and shared a more recent common ancestor to CMTV-amphibian associated ranaviruses reported in China, Italy, Denmark, and Switzerland. Field monitoring revealed differences in water frog host abundance at sites where distinct ranavirus groups occur; with ranavirus-associated deaths, host counts decreasing progressively, and few juveniles found in the north where CMTV-NL group I occurs but not in the south with CMTV-NL group III. Investigation of tandem repeats of coding genes gave no conclusive information about phylo-geographical clustering, while genetic analysis of the genomes revealed truncations in 17 genes across CMTV-NL groups II and III compared to group I. Further studies are needed to elucidate the contribution of these genes as well as environmental variables to explain the observed differences in host abundance.

Invertebrate Iridoviruses: A Glance over the Last Decade

Members of the family Iridoviridae (iridovirids) are large dsDNA viruses that infect both invertebrate and vertebrate ectotherms and whose symptoms range in severity from minor reductions in host fitness to systemic disease and large-scale mortality. Several characteristics have been useful for classifying iridoviruses; however, novel strains are continuously being discovered and, in many cases, reliable classification has been challenging. Further impeding classification, invertebrate iridoviruses (IIVs) can occasionally infect vertebrates; thus, host range is often not a useful criterion for classification. In this review, we discuss the current classification of iridovirids, focusing on genomic and structural features that distinguish vertebrate and invertebrate iridovirids and viral factors linked to host interactions in IIV6 (Invertebrate iridescent virus 6). In addition, we show for the first time how complete genome sequences of viral isolates can be leveraged to improve classification of new iridovirid isolates and resolve ambiguous relations. Improved classification of the iridoviruses may facilitate the identification of genus-specific virulence factors linked with diverse host phenotypes and host interactions.

The End of a 60-year Riddle: Identification and Genomic Characterization of an Iridovirus, the Causative Agent of White Fat Cell Disease in Zooplankton

The planktonic freshwater crustacean of the genus Daphnia are a model system for biomedical research and, in particular, invertebrate-parasite interactions. Up until now, no virus has been characterized for this system. Here we report the discovery of an iridovirus as the causative agent of White Fat Cell Disease (WFCD) in Daphnia WFCD is a highly virulent disease of Daphnia that can easily be cultured under laboratory conditions. Although it has been studied from sites across Eurasia for more than 60 years, its causative agent had not been described, nor had an iridovirus been connected to WFCD before now. Here we find that an iridovirus-the Daphnia iridescent virus 1 (DIV-1)-is the causative agent of WFCD. DIV-1 has a genome sequence of about 288 kbp, with 39% G+C content and encodes 367 predicted open reading frames. DIV-1 clusters together with other invertebrate iridoviruses but has by far the largest genome among all sequenced iridoviruses. Comparative genomics reveal that DIV-1 has apparently recently lost a substantial number of unique genes but has also gained genes by horizontal gene transfer from its crustacean host. DIV-1 represents the first invertebrate iridovirus that encodes proteins to purportedly cap RNA, and it contains unique genes for a DnaJ-like protein, a membrane glycoprotein and protein of the immunoglobulin superfamily, which may mediate host-pathogen interactions and pathogenicity. Our findings end a 60-year search for the causative agent of WFCD and add to our knowledge of iridovirus genomics and invertebrate-virus interactions.

Complete Genome Sequence of Frog virus 3, Isolated from a Strawberry Poison Frog (Oophaga pumilio) Imported from Nicaragua into the Netherlands

Frog virus 3 was isolated from a strawberry poison frog (Oophaga pumilio) imported from Nicaragua via Germany to the Netherlands, and its complete genome sequence was determined. Frog virus 3 isolate Op/2015/Netherlands/UU3150324001 is 107,183 bp long and has a nucleotide similarity of 98.26% to the reference Frog virus 3 isolate.

From fish to frogs and beyond: Impact and host range of emergent ranaviruses

Ranaviruses are pathogens of ectothermic vertebrates, including amphibians. We reviewed patterns of host range and virulence of ranaviruses in the context of virus genotype and postulate that patterns reflect significant variation in the historical and current host range of three groups of Ranavirus: FV3-like, CMTV-like and ATV-like ranaviruses. Our synthesis supports previous hypotheses about host range and jumps: FV3s are amphibian specialists, while ATVs are predominantly fish specialists that switched once to caudate amphibians. The most recent common ancestor of CMTV-like ranaviruses and FV3-like forms appears to have infected amphibians but CMTV-like ranaviruses may circulate in both amphibian and fish communities independently. While these hypotheses are speculative, we hope that ongoing efforts to describe ranavirus genetics, increased surveillance of host species and targeted experimental assays of susceptibility to infection and/or disease will facilitate better tests of the importance of hypothetical evolutionary drivers of ranavirus virulence and host range.

Molecular epidemiology of Epizootic haematopoietic necrosis virus (EHNV)

Low genetic diversity of Epizootic haematopoietic necrosis virus (EHNV) was determined for the complete genome of 16 isolates spanning the natural range of hosts, geography and time since the first outbreaks of disease. Genomes ranged from 125,591-127,487 nucleotides with 97.47% pairwise identity and 106-109 genes. All isolates shared 101 core genes with 121 potential genes predicted within the pan-genome of this collection. There was high conservation within 90,181 nucleotides of the core genes with isolates separated by average genetic distance of 3.43 × 10^-4 substitutions per site. Evolutionary analysis of the core genome strongly supported historical epidemiological evidence of iatrogenic spread of EHNV to naïve hosts and establishment of endemic status in discrete ecological niches. There was no evidence of structural genome reorganization, however, the complement of non-core genes and variation in repeat elements enabled fine scale molecular epidemiological investigation of this unpredictable pathogen of fish.

Comparative Genomics of Amphibian-like Ranaviruses, Nucleocytoplasmic Large DNA Viruses of Poikilotherms

Recent research on genome evolution of large DNA viruses has highlighted a number of incredibly dynamic processes that can facilitate rapid adaptation. The genomes of amphibian-like ranaviruses – double-stranded DNA viruses infecting amphibians, reptiles, and fish (family Iridoviridae) – were examined to assess variation in genome content and evolutionary processes. The viruses studied were closely related, but their genome content varied considerably, with 29 genes identified that were not present in all of the major clades. Twenty-one genes had evidence of recombination, while a virus isolated from a captive reptile appeared to be a mosaic of two divergent parents. Positive selection was also found to be acting on more than a quarter of Ranavirus genes and was found most frequently in the Spanish common midwife toad virus, which has had a severe impact on amphibian host communities. Efforts to resolve the root of this group by inclusion of an outgroup were inconclusive, but a set of core genes were identified, which recovered a well-supported species tree.