Seasonal dynamics and potential drivers of ranavirus epidemics in wood frog populations

Environmental DNA concentrations vary greatly across productive and degradative conditions, with implications for the precision of population estimates

Population size is an important metric to inform the conservation and management of species. For aquatic species, environmental DNA (eDNA) concentration has been suggested for non-invasively estimating population size. However, many biotic and abiotic factors simultaneously influence the production and degradation of eDNA which can alter the relationship between population size and eDNA concentration. We investigated the influence of temperature, salinity, and ranavirus infection on eDNA concentrations using tadpole mesocosms. Using linear regression models, we tested the influence of each experimental treatment on eDNA concentrations at three time points before and during epidemics. Prior to infection, elevated temperatures lowered eDNA concentrations, indicating that degradation was the driving force influencing eDNA concentrations. During early epidemics, no treatments strongly influenced eDNA concentrations and in late epidemics, productive forces dominated as ranavirus intensity and dead organisms increased eDNA concentrations. Finally, population size was only an important predictor of eDNA concentration in late epidemics and we observed high levels of variation between samples of replicate mesocosms. We demonstrate the complexities of several interacting factors influencing productive and degradative forces, variation in influences on eDNA concentration over short time spans, and examine the limitations of estimating population sizes from eDNA with precision in semi-natural conditions.

Ontogeny drives shifts in skin bacterial communities in facultatively paedomorphic salamanders

Microbiomes are major determinants of host growth, development and survival. In amphibians, host-associated bacteria in the skin can inhibit pathogen infection, but many processes can influence the structure and composition of the community. Here we quantified the shifts in skin-associated bacteria across developmental stages in the striped newt (Notophthalmus perstriatus), a threatened salamander species with a complex life history and vulnerable to infection by the amphibian chytrid fungus Batrachochytrium dendrobatidis and ranavirus. Our analyses show that pre-metamorphic larval and paedomorphic stages share similar bacterial compositions, and that the changes in the microbiome coincided with physiological restructuring during metamorphosis. Newts undergoing metamorphosis exhibited microbiome compositions that were intermediate between paedomorphic and post-metamorphic stages, further supporting the idea that metamorphosis is a major driver of host-associated microbes in amphibians. We did not find support for infection-related disruption of the microbiome, though infection replicates were small for each respective life stage.

The wood frog (Rana sylvatica): An emerging comparative model for anuran immunity and host-ranavirus interactions

The wood frog (Rana sylvatica) is widely distributed across North America and is the only amphibian found north of the Arctic Circle due to its remarkable ability to tolerate whole-body freezing. Recent mass mortalities attributable to Ranavirus spp. (family Iridoviridae) in wild juvenile wood frogs, coupled with the apparent high susceptibility of wood frogs to experimental infection with frog virus 3 (FV3), the type species of the Ranavirus genus, or FV3-like isolates underscore the serious threat ranaviruses poses to wood frog populations. Despite the ecological relevance and unique life history of wood frogs, our understanding of the wood frog immune system and antiviral response to ranaviral infections is in its infancy. Here we aim to (1) synthesize the limited knowledge of wood frog immune defences, (2) review recent progress in establishing the wood frog as a study system for ranavirus infection, and (3) highlight the future use of wood frogs as a model anuran to provide insight into the evolution of anuran immune systems and antiviral responses.

Widespread amphibian Perkinsea infections associated with Ranidae hosts, cooler months and Ranavirus co-infection

Amphibians suffer from large-scale population declines globally, and emerging infectious diseases contribute heavily to these declines. Amphibian Perkinsea (Pr) is a worldwide anuran pathogen associated with mass mortality events, yet little is known about its epidemiological patterns, especially in comparison to the body of literature on amphibian chytridiomycosis and ranavirosis. Here, we establish Pr infection patterns in natural anuran populations and identify important covariates including climate, host attributes and co-infection with Ranavirus (Rv). We used quantitative (q)PCR to determine the presence and intensity of Pr and Rv across 1234 individuals sampled throughout central Florida in 2017-2019. We then implemented random forest ensemble learning models to predict infection with both pathogens based on physiological and environmental characteristics. Perkinsea infected 32% of all sampled anurans, and Pr prevalence was significantly elevated in Ranidae frogs, cooler months, metamorphosed individuals and frogs co-infected with Rv, while Pr intensity was significantly higher in ranid frogs and individuals collected dead. Ranavirus prevalence was 17% overall and was significantly higher in Ranidae frogs, metamorphosed individuals, locations with higher average temperatures, and individuals co-infected with Pr. Perkinsea prevalence was significantly higher than Rv prevalence across months, regions, life stages and species. Among locations, Pr prevalence was negatively associated with crayfish prevalence and positively associated with relative abundance of microhylids, but Rv prevalence did not associate with any tested co-variates. Co-infections were significantly more common than single infections for both pathogens, and we propose that Pr infections may propel Rv infections because seasonal Rv infection peaks followed Pr infection peaks and random forest models found Pr intensity was a leading factor explaining Rv infections. Our study elucidates epidemiological patterns of Pr in Florida and suggests that Pr may be under-recognized as a cause of anuran declines, especially in the context of pathogen co-infection.

Diet- and salinity-induced modifications of the gut microbiota are associated with differential physiological responses to ranavirus infection in Rana sylvatica

Greater knowledge of how host-microbiome interactions vary with anthropogenic environmental change and influence pathogenic infections is needed to better understand stress-mediated disease outcomes. We investigated how increasing salinization in freshwaters (e.g. due to road de-icing salt runoff) and associated increases in growth of nutritional algae influenced gut bacterial assembly, host physiology and responses to ranavirus exposure in larval wood frogs (Rana sylvatica). Elevating salinity and supplementing a basic larval diet with algae increased larval growth and also increased ranavirus loads. However, larvae given algae did not exhibit elevated kidney corticosterone levels, accelerated development or weight loss post-infection, whereas larvae fed a basic diet did. Thus, algal supplementation reversed a potentially maladaptive stress response to infection observed in prior studies in this system. Algae supplementation also reduced gut bacterial diversity. Notably, we observed higher relative abundances of Firmicutes in treatments with algae-a pattern consistent with increased growth and fat deposition in mammals-that may contribute to the diminished stress responses to infection via regulation of host metabolism and endocrine function. Our study informs mechanistic hypotheses about the role of microbiome mediation of host responses to infection that can be tested in future experiments in this host-pathogen system. This article is part of the theme issue 'Amphibian immunity: stress, disease and ecoimmunology'.

Environmental DNA/RNA for pathogen and parasite detection, surveillance, and ecology

Detection of pathogens, parasites, and other symbionts in environmental samples via eDNA/eRNA (collectively eNA) is an increasingly important source of information about their occurrence and activity. There is great potential for using such detections as a proxy for infection of host organisms in connected habitats, for pathogen monitoring and surveillance, and for early warning systems for disease. However, many factors require consideration, and appropriate methods developed and verified, in order that eNA detections can be reliably interpreted and adopted for surveillance and assessment of disease risk, and potentially inclusion in international standards, such as the World Organisation for Animal Health guidelines. Disease manifestation results from host-symbiont-environment interactions between hosts, demanding a multifactorial approach to interpretation of eNA signals.

Pathogenesis of Bohle iridovirus infection in Krefft's freshwater turtle hatchlings (Emydura macquarii krefftii)

Ranaviruses have been detected in over 12 families of reptiles including many genera of turtles, tortoises, and terrapins, but the pathogenesis of these infections is still poorly understood. Krefft's river turtle hatchlings (N = 36; Emydura macquarii krefftii) were inoculated intramuscularly with Bohle iridovirus (BIV, Ranavirus, isolate) or saline, and euthanized at 9 timepoints (3 infected and 1 control per timepoint) over a 24-day period. Samples of lung, liver, kidney, and spleen were collected for quantitative polymerase chain reaction (PCR); internal organs, skin, and oral cavity samples were fixed for histopathological examination. The earliest lesions, at 8 days postinoculation (dpi), were lymphocytic inflammation of the skin and fibrinoid necrosis of regional vessels at the site of inoculation, and mild ulcerative necrosis with lymphocytic and heterophilic inflammation in the oral, nasal, and tongue mucosae. Fibrinonecrotic foci with heterophilic inflammation were detected in spleen and gonads at 16 dpi. Multifocal hepatic necrosis, heterophilic inflammation, and occasional basophilic intracytoplasmic inclusion bodies were observed at 20 dpi, along with ulcerative lymphocytic and heterophilic tracheitis and bronchitis. Tracheitis, bronchitis, and rare bone marrow necrosis were present at 24 dpi. Of the viscera tested for ranaviral DNA by PCR, the liver and spleen had the highest viral loads throughout infection, and thus appeared to be major targets of viral replication. Testing of whole blood by qPCR was the most-effective ante-mortem method for detecting ranaviral infection compared with oral swabs. This study represents the first time-dependent pathogenesis study of a ranaviral infection in turtles.

Ranavirus Amplification in Low-Diversity Amphibian Communities

In an era where emerging infectious diseases are a serious threat to biodiversity, epidemiological patterns need to be identified, particularly the complex mechanisms driving the dynamics of multi-host pathogens in natural communities. Many amphibian species have faced unprecedented population declines associated with diseases. Yet, specific processes shaping host-pathogen relationships within and among communities for amphibian pathogens such as ranaviruses (RV) remain poorly understood. To address this gap, we conducted a comprehensive study of RV in low-diversity amphibian communities in north-western Canada to assess the effects of biotic factors (species identity, species richness, abundance) and abiotic factors (conductivity, pH) on the pathogen prevalence and viral loads. Across 2 years and 18 sites, with communities of up to three hosts (wood frog, Rana sylvatica; boreal chorus frog, Pseudacris maculata; Canadian toad, Anaxyrus hemiophrys), we observed that RV prevalence nearly doubled with each additional species in a community, suggesting an amplification effect in aquatic, as well as terrestrial life-history stages. Infection intensity among infected wood frogs and boreal chorus frogs also significantly increased with an increase in species richness. Interestingly, we did not observe any effects of host abundance or abiotic factors, highlighting the importance of including host identity and species richness when investigating multi-host pathogens. Ultimately, only such a comprehensive approach can improve our understanding of complex and often highly context-dependent host-pathogen interactions.

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.

Host-multiparasite interactions in amphibians: a review

Parasites, including viruses, bacteria, fungi, protists, helminths, and arthropods, are ubiquitous in the animal kingdom. Consequently, hosts are frequently infected with more than one parasite species simultaneously. The assessment of such co-infections is of fundamental importance for disease ecology, but relevant studies involving non-domesticated animals have remained scarce. Many amphibians are in decline, and they generally have a highly diverse parasitic fauna. Here we review the literature reporting on field surveys, veterinary case studies, and laboratory experiments on co-infections in amphibians, and we summarize what is known about within-host interactions among parasites, which environmental and intrinsic factors influence the outcomes of these interactions, and what effects co-infections have on hosts. The available literature is piecemeal, and patterns are highly diverse, so that identifying general trends that would fit most host–multiparasite systems in amphibians is difficult. Several examples of additive, antagonistic, neutral, and synergistic effects among different parasites are known, but whether members of some higher taxa usually outcompete and override the effects of others remains unclear. The arrival order of different parasites and the time lag between exposures appear in many cases to fundamentally shape competition and disease progression. The first parasite to arrive can gain a marked reproductive advantage or induce cross-reaction immunity, but by disrupting the skin and associated defences (i.e., skin secretions, skin microbiome) and by immunosuppression, it can also pave the way for subsequent infections. Although there are exceptions, detrimental effects to the host are generally aggravated with increasing numbers of co-infecting parasite species. Finally, because amphibians are ectothermic animals, temperature appears to be the most critical environmental factor that affects co-infections, partly via its influence on amphibian immune function, partly due to its direct effect on the survival and growth of parasites. Besides their importance for our understanding of ecological patterns and processes, detailed knowledge about co-infections is also crucial for the design and implementation of effective wildlife disease management, so that studies concentrating on the identified gaps in our understanding represent rewarding research avenues.

Preparatory immunity: Seasonality of mucosal skin defences and Batrachochytrium infections in Southern leopard frogs

Accurately predicting the impacts of climate change on wildlife health requires a deeper understanding of seasonal rhythms in host-pathogen interactions. The amphibian pathogen, Batrachochytrium dendrobatidis (Bd), exhibits seasonality in incidence; however, the role that biological rhythms in host defences play in defining this pattern remains largely unknown. The aim of this study was to examine whether host immune and microbiome defences against Bd correspond with infection risk and seasonal fluctuations in temperature and humidity. Over the course of a year, five populations of Southern leopard frogs (Rana [Lithobates] sphenocephala) in Tennessee, United States, were surveyed for host immunity, microbiome and pathogen dynamics. Frogs were swabbed for pathogen load and skin bacterial diversity and stimulated to release stored antimicrobial peptides (AMPs). Secretions were analysed to estimate total hydrophobic peptide concentrations, presence of known AMPs and effectiveness of Bd growth inhibition in vitro. The diversity and proportion of bacterial reads with a 99% match to sequences of isolates known to inhibit Bd growth in vitro were used as an estimate of predicted anti-Bd function of the skin microbiome. Batrachochytrium dendrobatidis dynamics followed the expected seasonal fluctuations-peaks in cooler months-which coincided with when host mucosal defences were most potent against Bd. Specifically, the concentration and expression of stored AMPs cycled synchronously with Bd dynamics. Although microbiome changes followed more linear trends over time, the proportion of bacteria that can function to inhibit Bd growth was greatest when risk of Bd infection was highest. We interpret the increase in peptide storage in the fall and the shift to a more anti-Bd microbiome over winter as a preparatory response for subsequent infection risk during the colder periods when AMP synthesis and bacterial growth is slow and pathogen pressure from this cool-adapted fungus is high. Given that a decrease in stored AMP concentrations as temperatures warm in spring likely means greater secretion rates, the subsequent decrease in prevalence suggests seasonality of Bd in this host may be in part regulated by annual immune rhythms, and dominated by the effects of temperature.

A non-invasive method to assess environmental contamination with avian pathogens: beak and feather disease virus (BFDV) detection in nest boxes

Indirect transmission of pathogens can pose major risks to wildlife, yet the presence and persistence of wildlife pathogens in the environment has been little studied. Beak and feather disease virus (BFDV) is of global conservation concern: it can infect all members of the Psittaciformes, one of the most threatened bird orders, with infection often being lethal. Indirect transmission of BFDV through contaminated nest hollows has been proposed as a major infection source. However, data on whether and for how long nest sites in the wild remain contaminated have been absent. We determined the BFDV status of birds (parents and nestlings) for 82 nests of Crimson Rosellas, Platycercus elegans and Eastern Rosellas, Platycercus eximius. In 11 of these nests (13.4%, 95% confidence interval 6.9-22.7), we found an infected parent or nestling. Using nest swabs, we then compared BFDV presence at three points in time (before, during and after breeding) in three groups of nest boxes. These were nest boxes occupied by infected birds, and two control groups (nest boxes occupied by uninfected birds, and unoccupied nest boxes). Detection of BFDV on nest swabs was strongly associated with the infection status of parents in each nest box and with the timing of breeding. During breeding, boxes occupied by BFDV-positive birds were significantly more likely to have BFDV-positive nest swabs than boxes occupied by BFDV-negative birds; nest swabs tested BFDV-positive in 80% (28.4-99.5) of nests with parental antigen excretion, 66.7% (9.4-99.2) of nests occupied by parents with BFDV-positive cloacal swabs and 66.7% (22.3-95.7) of nests occupied by parents with BFDV-positive blood. 0% (0-52.2) of nests with BFDV-positive nestlings had BFDV-positive nest swabs. Across all boxes occupied by BFDV-positive birds (parents or nestlings), no nest swabs were BFDV-positive before breeding, 36.4% (95% CI 10.9-69.2) were positive during breeding and 9.1% (0.2-41.3) remained positive after breeding. BFDV was present on nest swabs for up to 3.7 months. Our study provides novel insights into the potential role of nest cavities and other fomites in indirect transmission of BFDV, and possibly other pathogens, and offers a non-invasive method for surveillance of pathogens in wild bird populations.

Salinity stress increases the severity of ranavirus epidemics in amphibian populations

The stress-induced susceptibility hypothesis, which predicts chronic stress weakens immune defences, was proposed to explain increasing infectious disease-related mass mortality and population declines. Previous work characterized wetland salinization as a chronic stressor to larval amphibian populations. Thus, we combined field observations with experimental exposures quantifying epidemiological parameters to test the role of salinity stress in the occurrence of ranavirus-associated mass mortality events. Despite ubiquitous pathogen presence (94%), populations exposed to salt runoff had slightly more frequent ranavirus related mass mortality events, more lethal infections, and 117-times greater pathogen environmental DNA. Experimental exposure to chronic elevated salinity (0.8-1.6 g l^-1 Cl^-) reduced tolerance to infection, causing greater mortality at lower doses. We found a strong negative relationship between splenocyte proliferation and corticosterone in ranavirus-infected larvae at a moderate elevation of salinity, supporting glucocorticoid-medicated immunosuppression, but not at high salinity. Salinity alone reduced proliferation further at similar corticosterone levels and infection intensities. Finally, larvae raised in elevated salinity had 10 times more intense infections and shed five times as much virus with similar viral decay rates, suggesting increased transmission. Our findings illustrate how a small change in habitat quality leads to more lethal infections and potentially greater transmission efficiency, increasing the severity of ranavirus epidemics.

Temperature-dependent infection of freshwater turtle hatchlings, Emydura macquarii krefftii, inoculated with a ranavirus isolate (Bohle iridovirus, Iridoviridae)

Fish, amphibians, and reptiles exhibit temperature-dependent ranaviral disease. We performed an experimental infection at four different environmental temperatures (16, 22, 28, and 34 °C) to investigate the effect of temperature on ranaviral infection in Krefft’s turtle ( Emydura macquarii krefftii). Infection rates and viral loads were determined by quantitative polymerase chain reaction to detect ranaviral DNA in liver samples at 21 d postexposure. The rate of infection differed across the temperature treatment groups. Infection rates were 44%, 90%, 60%, and 10% for the 16, 22, 28, and 34 °C temperature groups, respectively. Highest viral load was observed in the 28 °C temperature group, and there was a statistically significant difference in viral load between the 16 and 28 °C temperature groups ( p = 0.027). Based on the results of this study, the temperature of maximal infection rate for ranaviral infection in Krefft’s river turtles is estimated to be 23.2 °C (SD = 4.5). The findings of this study can inform management decisions in terms of disease control and treatment and form a platform for modelling disease outbreaks.

Detecting Amphibians in Agricultural Landscapes Using Environmental DNA Reveals the Importance of Wetland Condition

Amphibians are declining worldwide, in part because of large-scale degradation of habitat from agriculture and pervasive pathogens. Yet a common North American amphibian, the wood frog (Lithobates sylvaticus), ranges widely and persists in agricultural landscapes. Conventional survey techniques rely on visual encounters and dip-netting efforts, but detectability limits the ability to test for the effects of environmental variables on amphibian habitat suitability. We used environmental DNA to determine the presence of wood frogs and an amphibian pathogen (ranavirus) in Prairie Pothole wetlands and investigated the effects of 32 water quality, wetland habitat, and landscape-level variables on frog presence at sites representing different degrees of agricultural intensity. Several wetland variables influenced wood frog presence, the most influential being those associated with wetland productivity (i.e., nutrients), vegetation buffer width, and proportion of the surrounding landscape that is comprised of other water bodies. Wood frog presence was positively associated with higher dissolved phosphorus (>0.4 mg/L), moderate dissolved nitrogen (0.1-0.2 mg/L), lower chlorophyll a (≤15 µg/L), wider vegetation buffers (≥10 m), and more water on the landscape (≥0.25). These results highlight the effects of environmental factors at multiple scales on the presence of amphibians in this highly modified landscape-namely the importance of maintaining wetland water quality, vegetation buffers, and surrounding habitat heterogeneity. Environ Toxicol Chem 2019;38:2750-2763. © 2019 SETAC.

Effects of historic and projected climate change on the range and impacts of an emerging wildlife disease

The global trend of increasing environmental temperatures is often predicted to result in more severe disease epidemics. However, unambiguous evidence that temperature is a driver of epidemics is largely lacking, because it is demanding to demonstrate its role among the complex interactions between hosts, pathogens, and their shared environment. Here, we apply a three-pronged approach to understand the effects of temperature on ranavirus epidemics in UK common frogs, combining in vitro, in vivo, and field studies. Each approach suggests that higher temperatures drive increasing severity of epidemics. In wild populations, ranavirosis incidents were more frequent and more severe at higher temperatures, and their frequency increased through a period of historic warming in the 1990s. Laboratory experiments using cell culture and whole animal models showed that higher temperature increased ranavirus propagation, disease incidence, and mortality rate. These results, combined with climate projections, predict severe ranavirosis outbreaks will occur over wider areas and an extended season, possibly affecting larval recruitment. Since ranaviruses affect a variety of ectothermic hosts (amphibians, reptiles, and fish), wider ecological damage could occur. Our three complementary lines of evidence present a clear case for direct environmental modulation of these epidemics and suggest management options to protect species from disease.

Multi-tool diagnosis of an outbreak of ranavirosis in amphibian tadpoles in the Canadian boreal forest

Investigation of mortalities in isolated wild amphibian populations presents diagnostic difficulties that can hinder reaching a definitive diagnosis for the cause of death. Disease can only be diagnosed when pathogen presence (e.g. detection by PCR) is linked to tissue lesions (histopathology) in the host. We report a 2-site outbreak of ranavirosis in wild anuran tadpoles in the boreal forest of Wood Buffalo National Park, Canada, diagnosed by histologic and molecular techniques. Mortalities occurred in wood frog Rana sylvatica tadpoles and boreal chorus frog Pseudacris maculata tadpoles. Lack of mortality in sympatric Canadian toad Bufo (Anaxyrus) hemiophrys tadpoles suggested lower disease susceptibility in this species. In the former 2 species, ranavirosis was diagnosed based on consistent histopathology, immunohistochemistry (IHC), in situ hybridization (ISH), and quantitative PCR results. The most common histopathologic lesion present in wood and boreal chorus frog tadpoles was necrosis of the skin, oral mucosa, renal tubular epithelium, renal hematopoietic tissue, and branchial epithelium. Mild hepatic and pancreatic necrosis and rare intracytoplasmic inclusion bodies in hepatocytes were less common. Skeletal and connective tissues in budding limbs often had multifocal to coalescing necrosis and were intensely positive for ranavirus, with IHC staining even in areas where no obvious necrosis could be observed. Abundant IHC and ISH staining in actively growing tissues support a link between disease emergence and amphibian developmental stage. Our findings provide a definitive diagnosis of ranavirosis in free-living amphibians and highlight the effectiveness of multi-tool approaches to mortality investigation and elucidation of pathogenesis of ranavirosis in wild amphibians.

eDNA Increases the Detectability of Ranavirus Infection in an Alpine Amphibian Population

The early detection and identification of pathogenic microorganisms is essential in order to deploy appropriate mitigation measures. Viruses in the Iridoviridae family, such as those in the Ranavirus genus, can infect amphibian species without resulting in mortality or clinical signs, and they can also infect other hosts than amphibian species. Diagnostic techniques allowing the detection of the pathogen outside the period of host die-off would thus be of particular use. In this study, we tested a method using environmental DNA (eDNA) on a population of common frogs (Rana temporaria) known to be affected by a Ranavirus in the southern Alps in France. In six sampling sessions between June and September (the species' activity period), we collected tissue samples from dead and live frogs (adults and tadpoles), as well as insects (aquatic and terrestrial), sediment, and water. At the beginning of the breeding season in June, one adult was found dead; at the end of July, a mass mortality of tadpoles was observed. The viral DNA was detected in both adults and tadpoles (dead or alive) and in water samples, but it was not detected in insects or sediment. In live frog specimens, the virus was detected from June to September and in water samples from August to September. Dead tadpoles that tested positive for Ranavirus were observed only on one date (at the end of July). Our results indicate that eDNA can be an effective alternative to tissue/specimen sampling and can detect Ranavirus presence outside die-offs. Another advantage is that the collection of water samples can be performed by most field technicians. This study confirms that the use of eDNA can increase the performance and accuracy of wildlife health status monitoring and thus contribute to more effective surveillance programs.

Evaluating observer bias and seasonal detection rates in amphibian pathogen eDNA collections by citizen scientists

We trained volunteers from conservation organizations to collect environmental DNA (eDNA) from 21 ponds with amphibian communities that had a history of Batrachochytrium dendrobatidis (Bd) and ranavirus (Rv) infections. Volunteers were given sampling kits to filter pond water and preserve eDNA on filter paper, as were the principal investigators (PIs), who made independent collections within 48 h of volunteer collections. Using multi-scale occupancy modeling, we found no evidence to suggest the observer who collected the water sample (volunteer or PI) influenced either the probability of capturing eDNA on a filter or the probability of detecting extracted eDNA in a quantitative PCR (qPCR) reaction. The cumulative detection probability of Bd eDNA at a pond decreased from May through July 2017 because there was a decrease in the probability of detecting eDNA in qPCR reactions. In contrast, cumulative detection probability increased from May to July for Rv due to a higher probability of capturing eDNA on filters later in the year. Our models estimate that both pathogens could be detected with 95% confidence in as few as 5 water samples taken in June or July tested with either 4 or 3 qPCR reactions, respectively. Our eDNA protocols appeared to detect pathogens with 95% confidence using considerably fewer samples than protocols which typically recommend sampling ≥30 individual animals. In addition, eDNA sampling could reduce some biosecurity concerns, jurisdictional and institutional permitting, and stress to biota at ponds.