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

Testing the Population-Level Effects of Stress-Induced Susceptibility in the Ranavirus-Wood Frog System

Chronic exposure to physical, chemical, and biotic stressors can make animals more susceptible to infections. Such stress-induced susceptibility is widely expected to make disease-and epidemics-more likely and more severe. However, whether the impacts of stressors on individuals scale up to population-level outcomes is uncertain, both theoretically and empirically. We initiated ranavirus epidemics in replicate mesocosm populations of larval wood frogs (Lithobates sylvaticus) exposed to conditions known to impact their individual susceptibility to ranavirus infections: low and high salinity crossed factorially with ambient or elevated temperatures. Contrary to expectations, epidemics were no more likely or more severe in stressful conditions, although increased temperatures did speed their progression. We discuss several potential reasons why the effects of stressors did not scale up to epidemics, but our results suggest caution in assuming the individual-level effects of even well-studied stressors necessarily scale simply to population-level outcomes.

Leveraging machine learning to uncover multi-pathogen infection dynamics across co-distributed frog families

Background

Amphibians are experiencing substantial declines attributed to emerging pathogens. Efforts to understand what drives patterns of pathogen prevalence and differential responses among species are challenging because numerous factors related to the host, pathogen, and their shared environment can influence infection dynamics. Furthermore, sampling across broad taxonomic and geographic scales to evaluate these factors poses logistical challenges, and interpreting the roles of multiple potentially correlated variables is difficult with traditional statistical approaches. In this study, we leverage frozen tissues stored in natural history collections and machine learning techniques to characterize infection dynamics of three generalist pathogens known to cause mortality in frogs.

Methods

We selected 12 widespread and abundant focal taxa within three ecologically distinct, co-distributed host families (Bufonidae, Hylidae, and Ranidae) and sampled them across the eastern two-thirds of the United States of America. We screened and quantified infection loads via quantitative PCR for three major pathogens: the fungal pathogen Batrachochytrium dendrobatidis (Bd), double-stranded viruses in the lineage Ranavirus (Rv), and the alveolate parasite currently referred to as Amphibian Perkinsea (Pr). We then built balanced random forests (RF) models to predict infection status and intensity based on host taxonomy, age, sex, geography, and environmental variables and to assess relative variable importance across pathogens. Lastly, we used one-way analyses to determine directional relationships and significance of identified predictors.

Results

We found approximately 20% of individuals were infected with at least one pathogen (231 single infections and 25 coinfections). The most prevalent pathogen across all taxonomic groups was Bd (16.9%; 95% CI [14.9-19%]), followed by Rv (4.38%; 95% CI [3.35-5.7%]) and Pr (1.06%; 95% CI [0.618-1.82%]). The highest prevalence and intensity were found in the family Ranidae, which represented 74.3% of all infections, including the majority of Rv infection points, and had significantly higher Bd intensities compared to Bufonidae and Hylidae. Host species and environmental variables related to temperature were key predictors identified in RF models, with differences in importance among pathogens and host families. For Bd and Rv, infected individuals were associated with higher latitudes and cooler, more stable temperatures, while Pr showed trends in the opposite direction. We found no significant differences between sexes, but juvenile frogs had higher Rv prevalence and Bd infection intensity compared to adults. Overall, our study highlights the use of machine learning techniques and a broad sampling strategy for identifying important factors related to infection in multi-host, multi-pathogen systems.

Amphibian richness, rarity, threats, and conservation prospects across the U.S. National Park System

We assessed amphibian diversity, rarity, and threats across the National Park System (U.S.A.), which covers 3.5% of the country and 12% of federal lands. At least 230 of 354 (65%) amphibian species documented in the country occur on National Park Service lands. Of species in parks, 17% are at-risk globally and 20% are uncategorized, reflecting still-widespread data deficiencies. National parks in the Northwest and Northeast had the steepest species‒area relationships. Non-native crayfishes and amphibians occur within 50 km of 60% and 25% of parks, respectively, illustrating the broad threat of non-native predators. Projected mid-century (2040-2069) changes in climatic water deficit, based on 25 climate futures, produced an expected 34% increase in dryness across all national parks in the conterminous U.S.A. Our analyses highlight the extent and regional differences in current and future threats and reveal gaps in species protection, but also reveal opportunities for targeted expansion and active management.

High prevalence and concomitant infection of Ranavirus and Eustrongylides sp. in the invasive American Bullfrog in Brazil

American Bullfrogs, Aquarana catesbeiana, are invasive anuran species distributed worldwide. One of the adverse impacts that this species causes in native communities is as a reservoir host for pathogens and parasites. Here, we report the coinfection of two pathogenic organisms in A. catesbeiana: Ranavirus and the nematode Eustrongylides. Bullfrogs were collected in the wild in a pond close to the urban area of São Paulo, Brazil. The prevalence of both pathogens was high: 77% were infected with ranavirus with a mean viral load of 1010.3 viral copies, and 100% of the bullfrogs were infected by Eustrongylides sp. with a mean intensity of infection of 13.4 specimens per host. Four host specimens (31%) presented pathological signs that seemed to be related to the Eustrongylides sp. infection, such as internal organs adhered to each other due to high intensity and large size of the nematodes, ulcers, and raw flesh wounds caused by the nematode. The pathogenic and concomitant infections have potential zoonotic implications and raise concerns about human infection risks for Eustrongylides infection. Moreover, such infections may represent an additional level of threat to native communities through the potential shifts in patterns of parasite and pathogen transmission. Future research involving the native anuran community is essential to ascertain whether invasive bullfrogs are attenuating or exacerbating diseases such as ranavirosis and eustrongylidiosis.

Ranavirus and helminth parasite co-infection in invasive American bullfrogs in the Atlantic forest, Brazil

Emerging infectious diseases threaten amphibian species across the globe. In Brazil, the American bullfrog (Aquarana catesbeiana) is a highly invasive species that can potentially transmit parasites and pathogens to native amphibians. This is the first assessment of co-infection of Ranavirus and helminth macroparasites in invasive populations of bullfrogs in South America. We collected, measured, and euthanized 65 specimens of A. catesbeiana sampled from 9 sites across three states of Brazil in the Atlantic Forest biome. We collected and identified helminth macroparasites and sampled host liver tissue to test for the presence and load of Ranavirus with quantitative PCR. We documented patterns of prevalence, parasite load, and co-infection with generalized linear mixed models, generalized logistic regressions, and randomization tests. Most individual bullfrogs did not exhibit clinical signs of infection, but the overall Ranavirus prevalence was 27% (95% confidence interval, [CI 17-38]). Bullfrogs were infected with helminth macroparasites from 5 taxa. Co-infection of helminth macroparasites and Ranavirus was also common (21% CI [12-31]). Bullfrog size was positively correlated with total macroparasite abundance and richness, and the best-fitting model included a significant interaction between bullfrog size and Ranavirus infection status. We observed a negative correlation between Ranavirus viral load and nematode abundance (slope = -0.22, P = 0.03). Invasive bullfrogs (A. catesbeiana) in Brazil were frequently infected with both Ranavirus and helminth macroparasites, so adult bullfrogs could serve as reservoir hosts for both pathogens and parasites. However, many macroparasites collected were encysted and not developing. Coinfection patterns suggest a potential interaction between Ranavirus and macroparasites because helminth abundance increased with bullfrog size but was lower in Ranavirus infected individuals. Future studies of bullfrogs in the Atlantic Forest should investigate their potential role in pathogen and parasite transmission to native anurans.

Framework for multi-stressor physiological response evaluation in amphibian risk assessment and conservation

Controlled laboratory experiments are often performed on amphibians to establish causality between stressor presence and an adverse outcome. However, in the field, identification of lab-generated biomarkers from single stressors and the interactions of multiple impacts are difficult to discern in an ecological context. The ubiquity of some pesticides and anthropogenic contaminants results in potentially cryptic sublethal effects or synergistic effects among multiple stressors. Although biochemical pathways regulating physiological responses to toxic stressors are often well-conserved among vertebrates, different exposure regimes and life stage vulnerabilities can yield variable ecological risk among species. Here we examine stress-related biomarkers, highlight endpoints commonly linked to apical effects, and discuss differences in ontogeny and ecology that could limit interpretation of biomarkers across species. Further we identify promising field-based physiological measures indicative of potential impacts to health and development of amphibians that could be useful to anuran conservation. We outline the physiological responses to common stressors in the context of altered functional pathways, presenting useful stage-specific endpoints for anuran species, and discussing multi-stressor vulnerability in the larger framework of amphibian life history and ecology. This overview identifies points of physiological, ecological, and demographic vulnerability to provide context in evaluating the multiple stressors impacting amphibian populations worldwide for strategic conservation planning.

The amphibian's spleen as a source of biomarkers for ecotoxicity assessment: Historical review and trends

Amphibians are very sensitive to many environmental changes, so these animals are considered good bioindicator models for ecotoxicology. Given the importance of the amphibian spleen for hematopoietic and immune responses, this can be a key organ for the evaluation of biomarkers to monitor the health of individuals in nature or in captivity. In this systematic review, we searched databases and summarized the main findings concerning the amphibian spleen as a source of possible biomarkers applied in different scientific fields. The searches resulted in 83 articles published from 1923 to 2022, which applied the use of splenic samples to evaluate the effects of distinct stressors on amphibians. Articles were distributed in more than twenty countries, with USA, Europe, and Brazil, standing out among them. Publications focused mainly on anatomical and histomorphological characterization of the spleen, its physiology, and development. Recently, the use of splenic biomarkers in pathology and ecotoxicology began to grow but many gaps still need to be addressed in herpetological research. About 85 % of the splenic biomarkers showed responses to various stressors, which indicates that the spleen can provide numerous biomarkers to be used in many study fields. The limited amount of information on morphological description and splenic anatomy in amphibians may be a contributing factor to the underestimated use of splenic biomarkers in herpetological research around the world. We hope that this unprecedented review can instigate researchers to refine herpetological experimentation, using the spleen as a versatile and alternative source for biomarkers in ecotoxicology.

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.

Pathogenic selection promotes adaptive immune variations against serious bottlenecks in early invasions of bullfrogs

Adaptive genetic variations are key for understanding evolutionary processes influencing invasions. However, we have limited knowledge on how adaptive genetic diversity in invasive species responds to new pathogenic environments. Here, we compared variations in immune major histocompatibility complex (MHC) class-II β gene and neutral loci in relation to pathogenic chytrid fungus (Batrachochytrium dendrobatidis, Bd) infection across invasive and native populations of American bullfrog between China and United States (US). Chinese invasive populations show a 60% reduction in neutral cytb variations relative to US native populations, and there were similar MHC variation and functional diversity between them. One MHC allele private to China was under recent positive selection and associated with decreased Bd infection, partly explaining the lower Bd prevalence for Chinese populations than for native US populations. These results suggest that pathogen-mediated selection favors adaptive MHC variations and functional diversity maintenance against serious bottlenecks during the early invasions (within 15 generations) of bullfrogs.

Ranavirus infection does not reduce heat tolerance in a larval amphibian

Extreme heat events and emerging infectious diseases negatively impact wildlife populations, but the interacting effects of infection and host heat tolerance remain understudied. The few studies covering this subject have demonstrated that pathogens lower the heat tolerance of their hosts, which places infected hosts at a greater risk experiencing lethal heat stress. Here, we studied how ranavirus infection influenced heat tolerance in larval wood frogs (Lithobates sylvaticus). In line with similar studies, we predicted the elevated costs of ranavirus infection would lower heat tolerance, measured as critical thermal maximum (CT_max), compared to uninfected controls. Ranavirus infection did not reduce CT_max and there was a positive relationship between CT_max and viral loads. Our results demonstrate that ranavirus-infected wood frog larvae had no loss in heat tolerance compared to uninfected larvae, even at viral loads associated with high mortality rates, which contradicts the common pattern for other pathogenic infections in ectotherms. Larval anurans may prioritize maintenance of their CT_max when infected with ranavirus to promote selection of warmer temperatures during behavioral fever that can improve pathogen clearance. Our study represents the first to examine the effect of ranavirus infection on host heat tolerance, and because no decline in CT_max was observed, this suggests that infected hosts would not be under greater risk of heat stress.

A nested model with boosting and waning of immunity from Tilapia Lake Virus infection with distributed resistance to pathogens carrier-state

This paper proposes and analyzes an immune-structured population model of tilapia subject to Tilapia Lake Virus (TiLV) disease. The model incorporates within-host dynamics, used to describe the interaction between the pathogen, the immune system and the waning of immunity. Individuals infected with a low dose acquire a low immunity level and those infected with a high dose acquire a high level of immunity. Since individuals' immune status plays an important role in the spread of infectious diseases at the population level, the within-host dynamics are connected to the between-host dynamics in the population. We define an explicit formula for the reproductive number [Formula: see text] and show that the disease-free equilibrium is locally asymptotically stable when [Formula: see text], while it is unstable when [Formula: see text]. Furthermore, we prove that an endemic equilibrium exists. We also study the influence of the initial distribution of host resistance on the spread of the disease, and find that hosts' initial resistance plays a crucial role in the disease dynamics. This suggests that the genetic selection aiming to improve hosts' initial resistance to TiLV could help fight the disease. The results also point out the crucial role played by the inoculum size. We find that the higher the initial inoculum size, the faster the dynamics of infection. Moreover, if the initial inoculum size is below a certain threshold, it may not result in an outbreak at the between-host level. Finally, the model shows that there is a strong negative correlation between heterogeneity and the probability of pathogen invasion.

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.

First evidence of ranavirus in native and invasive amphibians in Colombia

Ranaviruses can cause mass mortality events in amphibians, thereby becoming a threat to populations that are already facing dramatic declines. Ranaviruses affect all life stages and persist in multiple amphibian hosts. The detrimental effects of ranavirus infections to amphibian populations have already been observed in the UK and in North America. In Central and South America, the virus has been reported in several countries, but the presence of the genus Ranavirus (Rv) in Colombia is unknown. To help fill this knowledge gap, we surveyed for Rv in 60 species of frogs (including one invasive species) in Colombia. We also tested for co-infection with Batrachochytrium dendrobatidis (Bd) in a subset of individuals. For Rv, we sampled 274 vouchered liver tissue samples collected between 2014 and 2019 from 41 localities covering lowlands to mountaintop páramo habitat across the country. Using quantitative polymerase chain reaction (qPCR) and end-point PCR, we detected Rv in 14 individuals from 8 localities, representing 6 species, including 5 native frogs of the genera Osornophryne, Pristimantis and Leptodactylus, and the invasive American bullfrog Rana catesbeiana. Bd was detected in 7 of 140 individuals, with 1 co-infection of Rv and Bd in an R. catesbeiana specimen collected in 2018. This constitutes the first report of ranavirus in Colombia and should set off alarms about this new emerging threat to amphibian populations in the country. Our findings provide some preliminary clues about how and when Rv may have spread and contribute to understanding how the pathogen is distributed globally.

Non-Lethal Detection of Frog Virus 3-Like (RUK13) and Common Midwife Toad Virus-Like (PDE18) Ranaviruses in Two UK-Native Amphibian Species

Ranaviruses have been involved in amphibian mass mortality events worldwide. Effective screening to control this pathogen is essential; however, current sampling methods are unsuitable for the detection of subclinical infections. Non-lethal screening is needed to prevent both further spread of ranavirus and losses of at-risk species. To assess non-lethal sampling methods, we conducted two experiments: bath exposing common frogs to RUK13 ranavirus at three concentrations, and exposing common toads to RUK13 or PDE18. Non-lethal sampling included buccal, digit, body and tank swabs, along with toe clips and stool taken across three time-points post-exposure. The presence/load of ranavirus was examined using quantitative PCR in 11 different tissues obtained from the same euthanised animals (incl. liver, gastro-intestinal tract and kidney). Buccal swab screening had the highest virus detection rate in both species (62% frogs; 71% toads) and produced consistently high virus levels compared to other non-lethal assays. The buccal swab was effective across multiple stages of infection and differing infection intensities, though low levels of infection were more difficult to detect. Buccal swab assays competed with, and even outperformed, lethal sampling in frogs and toads, respectively. Successful virus detection in the absence of clinical signs was observed (33% frogs; 50% toads); we found no difference in detectability for RUK13 and PDE18. Our results suggest that buccal swabbing could replace lethal sampling for screening and be introduced as standard practice for ranavirus surveillance.

Prevalence of Ranavirus Infection in Three Anuran Species across South Korea

To cope with amphibian die-offs caused by ranavirus, it is important to know the underlying ranavirus prevalence in a region. We studied the ranavirus prevalence in tadpoles of two native and one introduced anuran species inhabiting agricultural and surrounding areas at 49 locations across eight provinces of South Korea by applying qPCR. The local ranavirus prevalence and the individual infection rates at infected locations were 32.6% and 16.1%, respectively, for Dryophytes japonicus (Japanese tree frog); 25.6% and 26.1% for Pelophylax nigromaculatus (Black-spotted pond frog); and 30.5% and 50.0% for Lithobates catesbeianus (American bullfrog). The individual infection rate of L. catesbeianus was significantly greater than that of D. japonicus. The individual infection rate of P. nigromaculatus was related to the site-specific precipitation and air temperature. The individual infection rate gradually increased from Gosner development stage 39, and intermittent infection was confirmed in the early and middle developmental stages. Our results show that ranavirus is widespread among wild amphibians living in agricultural areas of South Korea, and mass die-offs by ranavirus could occur at any time.

Obligations of Researchers and Managers to Respect Wetlands: Practical Solutions to Minimizing Field Monitoring Impacts

Research and field monitoring can disturb wetland integrity. Adoption of ethical field practices is needed to limit monitoring induced stressors such as trampling, non-native seed and invertebrate dispersal, and disease and fungal spread. We identify a linear pathway of deterioration highlighting stressors that can progress to cumulative impacts, consequences, and losses at the site scale. The first step to minimize disturbance is to assess and classify the current ecosystem quality. We present a tiered framework for wetland classification and link preventative measures to the wetland tier. Preventative measures are recommended at various intensities respective to the wetland tier, with higher tiered wetlands requiring more intense preventative measures. In addition, preventative measures vary by time of implementation (before, during, and after the wetland visit) to mitigate impacts at various temporal scales. The framework is designed to increase transparency of field monitoring impacts and to promote the adoption of preventative measures. Implementing preventative measures can build accountability and foster a greater appreciation for our roles as researchers and managers in protecting wetlands.

First record of ranavirus (Ranavirus sp.) in Siberia, Russia

Ranaviruses are a group of double-strand DNA viruses that infect fish, amphibians and reptiles. These viruses are responsible for mass fish and amphibian mortality events worldwide, both in the wild and at fish and amphibian farms. The number of detected epizootics has grown significantly in recent years. In Eastern Europe and Northern Asia, including Russia, very few ranavirus monitoring studies have been conducted, in contrast with Western Europe and America. In the present work, we used a qPCR assay to survey for the first time the amphibian populations of West Siberia (Russia) for the presence of ranaviruses. In total, we studied 252 tissue samples from six amphibian species, collected across West Siberia from the south to the Arctic regions. We report a single infected sample: a common toad (Bufo bufo) captured near Tyumen city. The phylogenetic analysis showed that the detected virus strain belongs to the CMTV lineage. This is only the second observation of Ranavirus in Russia.

Pooled samples and eDNA-based detection can facilitate the "clean trade" of aquatic animals

The regional and international trade of live animals facilitates the movement, spillover, and emergence of zoonotic and epizootic pathogens around the world. Detecting pathogens in trade is critical for preventing their continued movement and introduction, but screening a sufficient fraction to ensure rare infections are detected is simply infeasible for many taxa and settings because of the vast numbers of animals involved—hundreds of millions of live animals are imported into the U.S.A. alone every year. Batch processing pools of individual samples or using environmental DNA (eDNA)—the genetic material shed into an organism’s environment—collected from whole consignments of animals may substantially reduce the time and cost associated with pathogen surveillance. Both approaches, however, lack a framework with which to determine sampling requirements and interpret results. Here I present formulae for pooled individual samples (e.g,. swabs) and eDNA samples collected from finite populations and discuss key assumptions and considerations for their use with a focus on detecting Batrachochytrium salamandrivorans, an emerging pathogen that threatens global salamander diversity. While empirical validation is key, these formulae illustrate the potential for eDNA-based detection in particular to reduce sample sizes and help bring clean trade into reach for a greater number of taxa, places, and contexts.