Ranavirus phylogenomics: Signatures of recombination and inversions among bullfrog ranaculture isolates

Hypersensitive response to interferon-stimulated gene (ISG)-inducing double-stranded RNA in American bullfrog tadpole fibroblasts

American bullfrogs are thought to be carriers of ranaviruses and contribute to their global spread via trade. Bullfrog tadpoles succumb to ranaviral infection's more severe and deadly effects than bullfrog adults. Presently, little is known about bullfrog tadpoles' innate antiviral immunity, possible due to the lack of available bullfrog tadpole cell lines. In this study, we describe a novel bullfrog tadpole fibroblast cell line named BullTad-leg. Its general cellular attributes, gene expression and function of class-A scavenger receptors (SR-As), and responses to poly IC (a synthetic dsRNA mimicking viral dsRNAs and a potent inducer of the interferon (IFN)-mediated antiviral responses) are investigated. Its abundant expression of vimentin corroborated with the cells' fibroblast morphology. BullTad-leg cells expressed transcripts of four SR-A members: SR-AI, SCARA3, SCARA4, and SCARA5, but transcripts of MARCO, the fifth SR-A member, were not detected. BullTad-leg cells expressed functional SR-As and could bind AcLDL. BullTad-leg cells exhibited cytotoxicity in response to poly IC treatment via SR-As. Additionally, very low doses of poly IC were able to induce dose-dependent expressions of ISGs including Mx, PKR, ISG20, and IFI35. This research sheds new light on the innate immune response, particularly SR-A biology and dsRNA responsiveness, in bullfrog tadpoles.

Fusing a TurboID tag with the Andrias davidianus ranavirus 2L reduced virus adsorption efficiency

Andrias davidianus ranavirus (ADRV) is a member of the genus ranavirus (family Iridoviridae). ADRV 2L is an envelope protein that could be essential in viral infection. In the present study, the function of ADRV 2L was investigated by fusion with the biotin ligase TurboID tag. A recombinant ADRV with a V5-TurboID tag fused in the N-terminal of 2L (ADRVT-2L) and a recombinant ADRV expressing V5-TurboID (ADRVT) were constructed, respectively. Infection of the recombinant viruses and wild-type ADRV (ADRVWT) in the Chinese giant salamander thymus cell line (GSTC) showed that ADRVT-2L had reduced cytopathic effect and lower virus titers than the other two viruses, indicating the fusion of a big tag affected ADRV infection. Analysis of the temporal expression profile showed that the expression of V5-TurboID-2L was delayed than wild-type 2L. However, electron microscopy found that the virion morphogenesis was not affected in ADRVT-2L-infected cells. Furthermore, the virus binding assay revealed that the adsorption efficiency of ADRVT-2L was considerably decreased compared to the other two viruses. Therefore, these data showed that linking the TurboID tag to ADRV 2L affected virus adsorption to the cell membrane, which suggested an important role of 2L in virus entry into cells.

Genome Characteristics of Two Ranavirus Isolates from Mandarin Fish and Largemouth Bass

Ranaviruses are promiscuous pathogens that threaten lower vertebrates globally. In the present study, two ranaviruses (SCRaV and MSRaV) were isolated from two fishes of the order Perciformes: mandarin fish (Siniperca chuatsi) and largemouth bass (Micropterus salmoides). The two ranaviruses both induced cytopathic effects in cultured cells from fish and amphibians and have the typical morphologic characteristics of ranaviruses. Complete genomes of the two ranaviruses were then sequenced and analyzed. Genomes of SCRaV and MSRaV have a length of 99, 405, and 99, 171 bp, respectively, and both contain 105 predicted open reading frames (ORFs). Eleven of the predicted proteins have differences between SCRaV and MSRaV, in which only one (79L) possessed a relatively large difference. A comparison of the sequenced six ranaviruses from the two fish species worldwide revealed that sequence identities of the six proteins (11R, 19R, 34L, 68L, 77L, and 103R) were related to the place where the virus was isolated. However, there were obvious differences in protein sequence identities between the two viruses and iridoviruses from other hosts, with more than half lower than 55%. Especially, 12 proteins of the two isolates had no homologs in viruses from other hosts. Phylogenetic analysis revealed that ranaviruses from the two fishes clustered in one clade. Further genome alignment showed five groups of genome arrangements of ranaviruses based on the locally collinear blocks, in which the ranaviruses, including SCRaV and MSRaV, constitute the fifth group. These results provide new information on the ranaviruses infecting fishes of Perciformes and also are useful for further research of functional genomics of the type of ranaviruses.

Genotypic Characterization of Infectious Spleen and Kidney Necrosis Virus (ISKNV) in Southeast Asian Aquaculture

Infectious spleen and kidney necrosis virus (ISKNV) is a species within the genus Megalocytivirus (family Iridoviridae), which causes high mortality disease in many freshwater and marine fish species. ISKNV was first reported in Asia and is an emerging threat to aquaculture with increasing global distribution, in part due to its presence in ornamental fish with clinical and subclinical infections. The species ISKNV includes three genotypes: red seabream iridovirus (RSIV), turbot reddish body iridovirus (TRBIV), and ISKNV. There is an increasing overlap in the recognized range of susceptible fish hosts and the geographic distribution of these distinct genotypes. To better understand the disease caused by ISKNV, a nucleic acid hybridization capture enrichment was used prior to sequencing to characterize whole genomes from archived clinical specimens of aquaculture and ornamental fish from Southeast Asia (n = 16). The method was suitable for tissue samples containing 2.50 × 104–4.58 × 109 ISKNV genome copies mg−1. Genome sequences determined using the hybridization capture method were identical to those obtained directly from tissues when there was sufficient viral DNA to sequence without enrichment (n = 2). ISKNV genomes from diverse locations, environments, and hosts had very high similarity and matched established genotype classifications (14 ISKNV genotype Clade 1 genomes with >98.81% nucleotide similarity). Conversely, two different genotypes were obtained at the same time and location (RSIV and ISKNV from grouper, Indonesia with 92.44% nucleotide similarity). Gene-by-gene analysis with representative ISKNV genomes identified 59 core genes within the species (>95% amino acid identity). The 14 Clade 1 ISKNV genomes in this study had 100% aa identity for 92–105 of 122 predicted genes. Despite high overall sequence similarity, phylogenetic analyses using single nucleotide polymorphisms differentiated isolates from different host species, country of origin, and time of collection. Whole genome studies of ISKNV and other megalocytiviruses enable genomic epidemiology and will provide information to enhance disease control in aquaculture.

Empirical evidence for effects of invasive American Bullfrogs on occurrence of native amphibians and emerging pathogens

Invasive species and emerging infectious diseases are two of the greatest threats to biodiversity. American Bullfrogs (Rana [Lithobates] catesbeiana), which have been introduced to many parts of the world, are often linked with declines in native amphibians via predation and the spread of emerging pathogens such as amphibian chytrid fungus (Batrachochytrium dendrobatidis [Bd]) and ranaviruses. Although many studies have investigated the potential role of bullfrogs in the decline of native amphibians, analyses that account for shared habitat affinities and imperfect detection have found limited support for clear effects. Similarly, the role of bullfrogs in shaping the patch-level distribution of pathogens is unclear. We used eDNA methods to sample 233 sites in the southwestern USA and Sonora, Mexico (2016-2018) to estimate how the presence of bullfrogs affects the occurrence of four native amphibians, Bd, and ranaviruses. Based on two-species, dominant-subordinate occupancy models fitted in a Bayesian context, federally threatened Chiricahua Leopard Frogs (Rana chiricahuensis) and Western Tiger Salamanders (Ambystoma mavortium) were eight times (32% vs. 4%) and two times (36% vs. 18%), respectively, less likely to occur at sites where bullfrogs occurred. Evidence for the negative effects of bullfrogs on Lowland Leopard Frogs (Rana yavapaiensis) and Northern Leopard Frogs (Rana pipiens) was less clear, possibly because of smaller numbers of sites where these native species still occurred and because bullfrogs often occur at lower densities in streams, the primary habitat for Lowland Leopard Frogs. At the community level, Bd was most likely to occur where bullfrogs co-occurred with native amphibians, which could increase the risk to native species. Ranaviruses were estimated to occur at 33% of bullfrog-only sites, 10% of sites where bullfrogs and native amphibians co-occurred, and only 3% of sites where only native amphibians occurred. Of the 85 sites where we did not detect any of the five target amphibian species, we also did not detect Bd or ranaviruses; this suggests other hosts do not drive the distribution of these pathogens in our study area. Our results provide landscape-scale evidence that bullfrogs reduce the occurrence of native amphibians and increase the occurrence of pathogens, information that can clarify risks and aid the prioritization of conservation actions.

Infectious disease threats to amphibians in Greece: new localities positive for Batrachochytrium dendrobatidis

In the early 2000s, numerous cases of European amphibian population declines and mass die-offs started to emerge. Investigating those events led to the discovery that wild European amphibians were confronted with grave disease threats caused by introduced pathogens, namely the amphibian and the salamander chytrid fungi Batrachochytrium dendrobatidis (Bd) and B. salamandrivorans (Bsal) and ranaviruses. In Greece, Bd was previously documented among wild amphibian populations in 2 different locations and 3 different species. However, no disease-related mass declines or mortality events have been reported. In this work, we build upon previous findings with new, subsequently obtained data, resulting in a 225-sample dataset of 14 species from 17 different locations throughout Greece, in order to examine the occurrence status of all 3 pathogens responsible for emerging infectious diseases in European amphibians. No positive samples for Bsal or ranavirus were recorded in any location. We confirmed the presence of Bd in 4 more localities and in 4 more species, including 1 urodelan (Macedonian crested newt Triturus macedonicus) and 1 introduced anuran (American bullfrog Lithobates catesbeianus). All insular localities were negative for Bd, except for Crete, where Bd was identified in 2 different locations. Again, no mass declines or die-offs were recorded in any Bd-positive area or elsewhere. However, given the persistence of Bd across Greece over the past ~20 yr, monitoring efforts should continue, and ideally be further expanded.

Climate warming triggers the emergence of native viruses in Iberian amphibians

The number of epizootics in amphibian populations caused by viruses of the genus Ranavirus is increasing worldwide. Yet, causes for pathogen emergence are poorly understood. Here, we confirmed that the Common midwife toad virus (CMTV) and Frog virus 3 (FV3) are responsible for mass mortalities in Iberia since the late 1980s. Our results illustrate the Iberian Peninsula as a diversity hotspot for the highly virulent CMTV. Although this pattern of diversity in Europe is consistent with spread by natural dispersal, the exact origin of the emergence of CMTV remains uncertain. Nevertheless, our data allow hypothesizing that the Iberian Peninsula might harbor the ancestral population of CMTVs that could have spread into the rest of Europe. In addition, we found that climate warming could be triggering the CMTV outbreaks, supporting its endemic status in the Iberian Peninsula.

Replication and transcription machinery for ranaviruses: components, correlation, and functional architecture

BACKGROUND

Ranaviruses (family Iridoviridae) are promiscuous pathogens that can infect across species barriers in poikilotherms and can replicate in amphibian and fish cells and even in cultured mammalian cells. However, as nucleocytoplasmic large DNA viruses (NCLDVs), their replication and transcription mechanisms remain largely unknown. Here, we screened and uncovered the replication and transcription machinery of two ranaviruses, Andrias davidianus ranavirus (ADRV) and Rana grylio virus (RGV), by a combination of methods, including the isolation of proteins on nascent DNA, recombinant virus-based affinity, and NanoLuc complementation assay.

RESULTS

The ranavirus replication and transcription machinery was deeply dissected and identified as a complicated apparatus containing at least 30 viral and 6 host proteins. The viral proteins ADRV-47L/RGV-63R (DNA polymerase, vDPOL), ADRV-23L/RGV-91R (proliferating cell nuclear antigen, vPCNA), ADRV-85L/RGV-27R (single-stranded DNA binding protein, vSSB), ADRV-88L/RGV-24R (vhelicase/primase), etc., constitute the core replisome. Specifically, the core of the transcription complex, the viral RNA polymerase, contain the host RNAPII subunits Rpb3, Rpb6, and Rpb11, which was a first report in NCLDVs. Furthermore, correlations and interactions among these factors in the machinery were described. Significantly, the replisome core protein vDPOL (ADRV-47L) can interact with numerous viral and host proteins and could act as a linker and regulation center in viral DNA replication and transcription. Thus, these results depicted an architecture for ranavirus replication and transcription.

CONCLUSIONS

Up to 36 components from ranavirus and their host were found to form viral replisomes and transcription complexes using a series of precise methods, which further constructed an architecture for ranavirus replication and transcription in which vDPOL was a key central factor and various components correlated and cooperated. Therefore, it provides a cornerstone for further understanding the mechanisms of the replication and transcription of ranaviruses which can ensure the efficient production of progeny virus and adaptation to cross-species infection.

Genomic Sequencing of Ranavirus Isolates from a Three-Spined Stickleback (Gasterosteus aculeatus) and a Red-Legged Frog (Rana aurora)

Two ranavirus isolates were recovered during a wildlife disease investigation in California in 1996. Preliminary testing at the time of sample collection indicated that the two isolates were identical. Phylogenetic analysis of the full genomes of these two isolates revealed that they are a single strain of frog virus 3.

Patterns of infection, origins, and transmission of ranaviruses among the ectothermic vertebrates of Asia

Ranaviral infections, a malady of ectothermic vertebrates, are becoming frequent, severe, and widespread, causing mortality among both wild and cultured species, raising odds of species extinctions and economic losses. This increase in infection is possibly due to the broad host range of ranaviruses and the transmission of these pathogens through regional and international trade in Asia, where outbreaks have been increasingly reported over the past decade. Here, we focus attention on the origins, means of transmission, and patterns of spread of this infection within the region. Infections have been recorded in both cultured and wild populations in at least nine countries/administrative regions, together with mass die-offs in some regions. Despite the imminent seriousness of the disease in Asia, surveillance efforts are still incipient. Some of the viral strains within Asia may transmit across host-taxon barriers, posing a significant risk to native species. Factors such as rising temperatures due to global climate change seem to exacerbate ranaviral activity, as most known outbreaks have been recorded during summer; however, data are still inadequate to verify this pattern for Asia. Import risk analysis, using protocols such as Pandora+, pre-border pathogen screening, and effective biosecurity measures, can be used to mitigate introduction of ranaviruses to uninfected areas and curb transmission within Asia. Comprehensive surveillance using molecular diagnostic tools for ranavirus species and variants will help in understanding the prevalence and disease burden in the region. This is an important step toward conserving native biodiversity and safeguarding the aquaculture industry.

Molecular Confirmation of Ranavirus Infection in Amphibians From Chad, Africa

Ranaviruses are DNA viruses (Family Iridoviridae; Subfamily Alphairidovirinae) and ranaviral disease is considered an emerging infectious disease of ectothermic vertebrates. Ranavirus infection can have varying pathological effects on infected amphibians, reptiles, and fish, most notably causing significant mortality events and population declines. Despite having a broad global range with reports from six continents, only a single incidental finding in Xenopus longipes from mainland Africa (Cameroon) is known and lacks molecular confirmation. Thus, there is a considerable knowledge gap concerning ranaviruses in Africa. We opportunistically obtained tissue samples from 160 amphibians representing five genera (Hoplobatrachus, Hylarana, Ptychadena, Pyxicephalus, and Xenopus) and two turtles (Pelomedusa sp.) from Chad, Africa. Samples were tested for ranavirus infection using a conventional PCR assay targeting the major capsid protein (MCP). A total of 25/160 (16%) frogs tested positive including 15/87 (17%) Hoplobatrachus occipitalis, 10/58 (17%) Ptychadena spp., 0/3 Pyxicephalus spp., 0/9 Xenopus spp., and 0/3 Hylarana spp. One of two turtles tested positive. Partial MCP gene sequences indicated all samples were >98% similar to several frog virus 3 (FV3)-like sequences. Additional gene targets (DNA polymerase [DNApol], ribonucleotide reductase alpha [RNR- α], ribonucleotide reductase beta subunit [RNR- β]) were sequenced to provide further detailed classification of the virus. Sequences of individual gene targets indicate that the ranavirus detected in frogs in Chad is most similar to tiger frog virus (TFV), a FV3-like virus previously isolated from diseased amphibians cultured in China and Thailand. Full genome sequencing of one sample indicates that the Chad frog virus (CFV) is a well-supported sister group to the TFVs previously determined from Asia. This work represents the first molecular confirmation of ranaviruses from Africa and is a first step in comparing ranavirus phylogeography on a local and global scale.

Genetic diversity of the carp oedema virus in France

The koi sleepy disease of carp caused by the carp oedema virus (CEV) was observed on farms and in ponds in France since the 2010s. Samples of CEV collected in France over a period of eight years were characterized at the molecular level by sequencing the partial p4a gene. All the sequences, except one, fell into two well-defined genogroups. Sequences obtained from CEV detected in common carp generally clustered in genogroup I and sequences from CEV detected in the koi were assigned to genogroup II. A particular sample was different to the others and represented a putative new genogroup possibly arose from a recombination event between a genogroup II sequence and one from an unknown genogroup. Compared with sequences from CEV of other countries, most of the French sequences exhibited high degree of DNA identities with those published previously, indicating identical sources of viruses. The sequence diversity suggests multiple introductions of the viruses in France. Among the French sequences, two genogroup-specific molecular markers were identified. One was an insertion/deletion identified within a microsatellite and other was a group of single nucleotide polymorphisms. CEV seems to generate genetic diversity via diverse mechanisms: substitutions, indels and recombination events.

Genomic sequencing of a frog virus 3 strain from cultured American bullfrogs (Lithobates catesbeianus) in Brazil

Frog virus 3 (FV3) was detected in cultured bullfrogs in Southeast Brazil. Phylodynamic analysis revealed recombination events in this strain that were nearly identical to those detected in North American and Brazilian FV3 strains. These data suggest that international trade of live bullfrogs has spread recombinant strains of FV3.

FATAL RANAVIRUS INFECTION IN A GROUP OF ZOO-HOUSED MELLER'S CHAMELEONS (TRIOCEROS MELLERI)

A group of five juvenile Meller's chameleons (Trioceros melleri) experienced 100% mortality over a period of 1 mo due to ranavirus infection. The index case was found dead without premonitory signs. The three subsequent cases presented with nonspecific clinical signs (lethargy, decreased appetite, ocular discharge) and were ultimately euthanatized. The final case died after initially presenting with skin lesions. Postmortem examination revealed thin body condition in all five animals and mild coelomic effusion and petechiae affecting the tongue and kidneys of one animal. Microscopically, all animals had multifocal necrosis of the spleen, liver, and kidney; four of five animals had necrosis of the nasal cavity; and two of five had necrosis of adrenal tissue, bone marrow, and skin. Numerous basophilic intracytoplasmic inclusions were present in the liver of all animals and nasal mucosa of three of the five animals. Consensus polymerase chain reaction for herpesvirus and adenovirus were negative, whereas ranavirus quantitative polymerase chain reaction was positive. Virus isolation followed by whole genome sequencing and Bayesian phylogenetic analysis classified the isolates as a strain of frog virus 3 (FV3) most closely related to an FV3 isolate responsible for a previous outbreak in the zoo's eastern box turtle (Terrapene carolina carolina) group. This case series documents the first known occurrence of ranavirus-associated disease in chameleons and demonstrates the potential for interspecies transmission between chelonian and squamate reptiles.

Global patterns of ranavirus detections

Ranaviruses are emerging pathogens of poikilothermic vertebrates. In 2015 the Global Ranavirus Reporting System (GRRS) was established as a centralized, open access, online database for reports of the presence (and absence) of ranavirus around the globe. The GRRS has multiple data layers (e.g., location, date, host(s) species, and methods of detection) of use to those studying the epidemiology, ecology, and evolution of this group of viruses. Here we summarize the temporal, spatial, diagnostic, and host-taxonomic patterns of ranavirus reports in the GRRS. The number, distribution, and host diversity of ranavirus reports have increased dramatically since the mid 1990s, presumably in response to increased interest in ranaviruses and the conservation of their hosts, and also the availability of molecular diagnostics. Yet there are clear geographic and taxonomic biases among the reports. We encourage ranavirus researchers to add their studies to the portal because such collation can provide collaborative opportunities and unique insights to our developing knowledge of this pathogen and the emerging infectious disease that it causes.

Genome Sequence of a Ranavirus Isolated from a Red-Eared Slider (Trachemys scripta elegans) in Poland

The red-eared slider (RES) ranavirus (RESRV) was isolated from a free-ranging RES turtle that died with evidence of respiratory disease. The RESRV genome sequence (106,878 bp) was determined, and phylogenetic analysis revealed that it is a common midwife toad virus (CMTV) strain. This study is the first report of CMTV in RES.

First genetic characterization of sturgeon mimiviruses in Ukraine

A group of pathogenic nucleocytoplasmic large DNA viruses (NCLDVs) related to the Mimiviridae family infect farmed sturgeons across Europe, causing mild-to-severe losses. One of these viruses, Acipenser iridovirus-European (AcIV-E), was identified in six sturgeon species. During the 2018-2019 period, nine sick Siberian (A. baerii) and Russian (A. gueldenstaedtii) sturgeons were sampled in Ukrainian farms and tested for the presence of AcIV-E using real-time PCR. The presence of AcIV-E was confirmed in some samples. High-resolution melting (HRM) assay and Sanger sequencing demonstrated the presence in three farms of two alleles of the major capsid protein (MCP) gene, called var1 and var2. Five samples carried both var1 and var2 at varying ratios, and the sixth sample was infected with only var1. These results constitute the first detection of AcIV-E in Ukraine and the first detection of a sample carrying only var1. The full-length sequences of the MCP genes confirmed the existence of two genetic lineages of AcIV-E, tentatively named V1 and V2, each displaying multiple substitutions in the MCP gene. Some of the MCP sequences showed a genetic relationship to both V1 and V2 lineages, depending on the fragment examined. Most likely, these sequences resulted from recombination events.

Systems Biology Approaches for Therapeutics Development Against COVID-19

Understanding the systems biology approaches for promoting the development of new therapeutic drugs is attaining importance nowadays. The threat of COVID-19 outbreak needs to be vanished for global welfare, and every section of research is focusing on it. There is an opportunity for finding new, quick, and accurate tools for developing treatment options, including the vaccine against COVID-19. The review at this moment covers various aspects of pathogenesis and host factors for exploring the virus target and developing suitable therapeutic solutions through systems biology tools. Furthermore, this review also covers the extensive details of multiomics tools i.e., transcriptomics, proteomics, genomics, lipidomics, immunomics, and in silico computational modeling aiming towards the study of host-virus interactions in search of therapeutic targets against the COVID-19.

Selected Emerging Infectious Diseases of Amphibians

This article updates the understanding of two extirpation-driving infectious diseases, Batrachochytrium dendrobatidis and Batrachochytrium salamandrivorans, and Ranavirus. Experimental studies and dynamic, multifactorial population modeling have outlined the epidemiology and future population impacts of B dendrobatidis, B salamandrivorans, and Ranavirus. New genomic findings on divergent fungal and viral pathogens can help optimize control and disease management strategies. Although there have been major advances in knowledge of amphibian pathogens, controlled studies are needed to guide population recovery to elucidate and evaluate transmission routes for several pathogens, examine environmental control, and validate new diagnostic tools to confirm the presence of disease.

Genetic and codon usage bias analyses of major capsid protein gene in Ranavirus

The Ranavirus (one genus of Iridovidae family) is an emerging pathogen that infects fish, amphibian, and reptiles, and causes great economical loss and ecological threat to farmed and wild animals globally. The major capsid protein (MCP) has been used as genetic typing marker and as target to design vaccines. Herein, the codon usage pattern of 73 MCP genes of Ranavirus and Lymphocystivirus are studied by calculating effective number of codons (ENC), relative synonymous codon usage (RSCU), codon adaptation index (CAI), and relative codon deoptimization index (RCDI), and similarity index (SiD). The Ranavirus are confirmed to be classified into five groups by using phylogenetic analysis, and varied nucleotide compositions and hierarchical cluster analysis based on RSCU. The results revealed different codon usage patterns among Lymphocystivirus and five groups of Ranavirus. Ranavirus had six over-represented codons ended with G/C nucleotide, while Lymphocystivirus had six over-represented codons ended with A/T nucleotide. A comparative analysis of parameters that define virus and host relatedness in terms of codon usage were analyzed indicated that Amphibian-like ranaviruses (ALRVs) seem to possess lower ENC values and higher CAIs in contrast to other ranaviruses isolated from fishes, and two groups (FV3-like and CMTV-like group) of them had received higher selection pressure from their hosts as having higher relative codon deoptimization index (RCDI) and similarity index (SiD). The correspondence analysis (COA) and Spearman's rank correlation analyses revealed that nucleotide compositions, relative dinucleotide frequency, mutation pressure, and natural translational selection shape the codon usage pattern in MCP genes and the ENC-GC_3S and neutrality plots indicated that the natural selection is the predominant factor. These results contribute to understanding the evolution of Ranavirus and their adaptions to their hosts.

Depauperate major histocompatibility complex variation in the endangered reticulated flatwoods salamander (Ambystoma bishopi)

Reticulated flatwoods salamander (Ambystoma bishopi) populations began decreasing dramatically in the 1900s. Contemporary populations are small, isolated, and may be susceptible to inbreeding and reduced adaptive potential because of low genetic variation. Genetic variation at immune genes is especially important as it influences disease susceptibility and adaptation to emerging infectious pathogens, a central conservation concern for declining amphibians. We collected samples from across the extant range of this salamander to examine genetic variation at major histocompatibility complex (MHC) class Iα and IIβ exons as well as the mitochondrial control region. We screened tail or toe tissue for ranavirus, a pathogen associated with amphibian declines worldwide. Overall, we found low MHC variation when compared to other amphibian species and did not detect ranavirus at any site. MHC class Iα sequencing revealed only three alleles with a nucleotide diversity of 0.001, while MHC class IIβ had five alleles with a with nucleotide diversity of 0.004. However, unique variation still exists across this species' range with private alleles at three sites. Unlike MHC diversity, mitochondrial variation was comparable to levels estimated for other amphibians with nine haplotypes observed, including one haplotype shared across all sites. We hypothesize that a combination of a historic disease outbreak and a population bottleneck may have contributed to low MHC diversity while maintaining higher levels of mitochondrial DNA variation. Ultimately, MHC data indicated that the reticulated flatwoods salamander may be at an elevated risk from infectious diseases due to low levels of immunogenetic variation necessary to combat novel pathogens.

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.

Amphibian ranaviruses in Europe: important directions for future research

Ranaviruses are an emerging group of pathogens capable of infecting all cold-blooded vertebrates. In Europe, ranaviruses pose a particularly potent threat to wild amphibian populations. Since the 1980s research on amphibian-infecting ranaviruses in Europe has been growing. The wide distribution of amphibian populations in Europe, the ease with which many are monitored, and the tractable nature of counterpart ex situ experimental systems have provided researchers with a unique opportunity to study many aspects of host–ranavirus interactions in the wild. These characteristics of European amphibian populations will also enable researchers to lead the way as the field of host–ranavirus interactions progresses. In this review, we provide a summary of the current key knowledge regarding amphibian infecting ranaviruses throughout Europe. We then outline important areas of further research and suggest practical ways each could be pursued. We address the study of potential interactions between the amphibian microbiome and ranaviruses, how pollution may exacerbate ranaviral disease either as direct stressors of amphibians or indirect modification of the amphibian microbiome. Finally, we discuss the need for continued surveillance of ranaviral emergence in the face of climate change.

Phylogenomic characterization of ranaviruses isolated from cultured fish and amphibians in Thailand

Ranaviruses are emerging pathogens associated with worldwide epizootics in farmed and wild ectothermic vertebrates. In this study, we determined the full genomes of eight ranaviruses isolated from marbled sleeper goby ( Oxyeleotris marmorata), goldfish ( Carassius auratus), guppy ( Poecilia reticulata), tiger frog ( Hoplobatrachus tigerinus), Asian grass frog ( Fejervarya limnocharis), and East Asian bullfrog ( H. rugulosus) cultured or imported into Thailand. These ranaviral isolates induced the same cytopathic effects (i.e., progression of coalescing round plaques) in epithelioma papulosum cyprini (EPC) cell cultures. Transmission electron microscopy of infected EPC cells revealed cytoplasmic viral particles with ultrastructural features typical for ranaviruses. Pairwise genetic comparisons of the complete major capsid protein coding sequences from the Thai ranaviruses displayed the highest identity (99.8%–100%) to a ranavirus (tiger frog virus; TFV) isolated from diseased tiger frogs cultured in China, a slightly lower identity (99.3%–99.4%) to a ranavirus (Wamena virus; WV) isolated from diseased green tree pythons ( Morelia viridis) illegally exported from Papua New Guinea, and a lower identity to 35 other ranaviruses (93.7%–98.6%). Phylogenomic analyses supported the eight Thai ranaviruses, Chinese TFV, and WV as a subclade within a larger frog virus 3 clade. Our findings confirm the spread of TFV among cultured fish and amphibians in Asia and likely in reptiles in Oceania. Biosecurity measures are needed to ensure TFV does not continue to spread throughout Southeast Asia and to other parts of the world via international trade.

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.

Ranavirus infection dynamics and shedding in American bullfrogs: consequences for spread and detection in trade

American bullfrogs Lithobates catesbeianus are thought to be important in the global spread of ranaviruses-often lethal viruses of cold-blooded vertebrates-because they are commonly farmed, dominate international trade, and may be 'carriers' of ranavirus infections. However, whether American bullfrogs are easily infected and maintain long-lasting ranavirus infections, or are refractory to or rapidly clear infections, remains unknown. We tracked the dynamics of ranavirus in American bullfrogs through time and with temperature in multiple types of samples and also screened shipments from commercial suppliers to determine whether we could detect subclinical infections. Collectively, we found that tadpoles and juveniles were commonly infected at moderate doses, and while some died, others controlled and appeared to clear their infections. Some individuals, however, harbored subclinical infections for up to 49 d, suggesting that American bullfrogs may be important carriers. Indeed, tadpoles and metamorphosed frogs from 2 of 5 commercial suppliers harbored subclinicial infections. Juveniles at warmer temperatures had less intense but still persistent infections. Because diagnostic performance was strongly related to infection intensity, non-lethal samples (i.e. tail or toe clips, swabs, and environmental DNA) had only a moderate chance of detecting subclinical infections. Even internal tissues may fail to detect subclinical infections. However, viral shedding was correlated with the intensity of infection, so while subclinically infected tadpoles shed virus for 35-49 d, the low levels might lead to little transmission. We suggest that a quantitative focus on virus dynamics within hosts can provide a more nuanced view of ranavirus infections and the risk presented by American bullfrogs in trade.

Characterization of a Novel Megalocytivirus Isolated from European Chub (Squalius cephalus)

A novel virus from moribund European chub (Squalius cephalus) was isolated on epithelioma papulosum cyprini (EPC) cells. Transmission electron microscopic examination revealed abundant non-enveloped, hexagonal virus particles in the cytoplasm of infected EPC cells consistent with an iridovirus. Illumina MiSeq sequence data enabled the assembly and annotation of the full genome (128,216 bp encoding 108 open reading frames) of the suspected iridovirus. Maximum Likelihood phylogenetic analyses based on 25 iridovirus core genes supported the European chub iridovirus (ECIV) as being the sister species to the recently-discovered scale drop disease virus (SDDV), which together form the most basal megalocytivirus clade. Genetic analyses of the ECIV major capsid protein and ATPase genes revealed the greatest nucleotide identity to members of the genus Megalocytivirus including SDDV. These data support ECIV as a novel member within the genus Megalocytivirus. Experimental challenge studies are needed to fulfill River’s postulates and determine whether ECIV induces the pathognomonic microscopic lesions (i.e., megalocytes with basophilic cytoplasmic inclusions) observed in megalocytivirus infections.

Low intraspecific variation of Frog virus 3 with evidence for novel FV3-like isolates in central and northwestern Canada

Frog virus 3 (FV3) and FV3-like ranaviruses can infect a variety of cold-blooded aquatic species and present a primary threat to amphibians across the globe. Previous studies of FV3-like viruses have largely investigated higher-level phylogenetic distinctions of these pathogens via portions of the conserved major capsid protein (MCP), and the putative virulence gene vIF-2α. Few studies, however, have investigated the spatial distribution of FV3 variants at the population level3-data that can be used to further understand the spatial epidemiology of this disease. In this study, we sequenced the MCP and vIF-2α of 127 FV3-positive amphibians sampled from Canadian water bodies in Ontario, northeastern Alberta, and southern Northwest Territories to explore whether intraspecific genetic variation exists within FV3. There was a lack of variation at the 2 markers across these regions, suggesting that there is a lack of FV3 sequence diversity in Canada, which may hint at a single source of infection that has spread. However, an undocumented variant termed Wood Buffalo ranavirus (WBRV) was detected in samples from 3 sites in Alberta and Northwest Territories that clustered within the FV3-like lineage with 99.3% sequence homology for MCP. For vIF-2α, all sequences were the expected truncated variant except for 6 samples in Ontario. These latter sequences were suggestive of recombination with common midwife toad virus (CMTV). The lack of variation suggests that higher-resolution genome analyses will be required to further explore the spatial spread and intraspecific variation of the disease.

Ranaviruses and reptiles

Ranaviruses can infect many vertebrate classes including fish, amphibians and reptiles, but for the most part, research has been focused on non-reptilian hosts, amphibians in particular. More recently, reports of ranaviral infections of reptiles are increasing with over 12 families of reptiles currently susceptible to ranaviral infection. Reptiles are infected by ranaviruses that are genetically similar to, or the same as, the viruses that infect amphibians and fish; however, physiological and ecological differences result in differences in study designs. Although ranaviral disease in reptiles is often influenced by host species, viral strain and environmental differences, general trends in pathogenesis are emerging. More experimental studies using a variety of reptile species, life stages and routes of transmission are required to unravel the complexity of wild ranavirus transmission. Further, our understanding of the reptilian immune response to ranaviral infection is still lacking, although the considerable amount of work conducted in amphibians will serve as a useful guide for future studies in reptiles.

Co-exposure to multiple ranavirus types enhances viral infectivity and replication in a larval amphibian system

Multiple pathogens commonly co-occur in animal populations, yet few studies demonstrate how co-exposure of individual hosts scales up to affect transmission. Although viruses in the genus Ranavirus are globally widespread, and multiple virus species or strains likely co-occur in nature, no studies have examined how co-exposure affects infection dynamics in larval amphibians. We exposed individual northern red-legged frog Rana aurora larvae to 2 species of ranavirus, namely Ambystoma tigrinum virus (ATV), frog virus 3 (FV3), or an FV3-like strain isolated from a frog-culturing facility in Georgia, USA (RCV-Z2). We compared single-virus to pairwise co-exposures while experimentally accounting for dosage. Co-exposure to ATV and FV3-like strains resulted in almost twice as many infected individuals compared to single-virus exposures, suggesting an effect of co-exposure on viral infectivity. The viral load in infected individuals exposed to ATV and FV3 was also higher than the single-dose FV3 treatment, suggesting an effect of co-exposure on viral replication. In a follow-up experiment, we examined how the co-occurrence of ATV and FV3 affected epizootics in mesocosm populations of larval western chorus frogs Pseudacris triseriata. Although ATV did not generally establish within host populations (<4% prevalence), when ATV and FV3 were both present, this co-exposure resulted in a larger epizootic of FV3. Our results emphasize the importance of multi-pathogen interactions in epizootic dynamics and have management implications for natural and commercial amphibian populations.

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.

Poor biosecurity could lead to disease outbreaks in animal populations

Human-mediated disease outbreaks due to poor biosecurity practices when processing animals in wild populations have been suspected. We tested whether not changing nitrile gloves between processing wood frog (Lithobates sylvaticus) tadpoles and co-housing individuals increased pathogen transmission and subsequent diseased-induced mortality caused by the emerging pathogen, ranavirus. We found that not changing gloves between processing infected and uninfected tadpoles resulted in transmission of ranavirus and increased the risk of mortality of uninfected tadpoles by 30X. Co-housing tadpoles for only 15 minutes with 10% of individuals infected resulted in ranavirus transmission and 50% mortality of uninfected tadpoles. More extreme mortality was observed when the co-housing infection prevalence was >10%. Our results illustrate that human-induced disease outbreaks due to poor biosecurity practices are possible in wild animal populations.

Genomic Sequencing of Ranaviruses Isolated from Edible Frogs (Pelophylax esculentus)

Ranaviruses were isolated from wild edible frogs (Pelophylax esculentus) during epizootics in Denmark and Italy. Phylogenomic analyses revealed that these isolates are closely related and belong to a clade of ranaviruses that includes the Andrias davidianus ranavirus (ADRV), common midwife toad ranavirus (CMTV), Testudo hermanni ranavirus (THRV), and pike-perch iridovirus (PPIV).