Differences in sensitivity to the fungal pathogen Batrachochytrium dendrobatidis among amphibian populations

Host life stage susceptibility to the chytrid fungus in a Neotropical torrent frog

Pathogen-host systems become complex when they involve life histories with multiple stages. Understanding these complexities is particularly important for investigating the infection dynamics of the amphibian pathogen Batrachochytrium dendrobatidis (Bd). Here, we investigate whether Bd infection susceptibility differs between host life stages and determine the influence of environmental factors on Bd infection rates across remnant populations of a Neotropical torrent frog. We found that Bd infection probability varies between tadpoles and adults in Hylodes phyllodes, with tadpoles exhibiting a higher likelihood of infection. Tadpoles are tolerant to Bd, acting as zoospore reservoirs, potentially aiding in the pathogen's persistence in the environment and infecting other susceptible hosts. Topographic complexity, species richness, the human footprint, precipitation seasonality and diurnal temperature variations influenced Bd infection rates. Conservation strategies should encompass both host life stages, monitoring from larvae to adults, while also evaluating threats synergistically, such as the human footprint, to effectively predict and mitigate the impact of Bd on susceptible populations.

Differential recruitment drives pathogen-mediated competition between species in an amphibian chytridiomycosis system

Pathogens that infect multiple host species have an increased capacity to cause extinctions through parasite-mediated apparent competition. Given unprecedented and continuing losses of biodiversity due to Batrachochytrium dendrobatidis (Bd), the causative fungus of the amphibian skin disease chytridiomycosis, a robust understanding of the mechanisms driving cross-species infection dynamics is essential. Here, we used stage-structured, susceptible-infected compartmental models to explore drivers of Bd-mediated apparent competition between two sympatric amphibians, the critically endangered Litoria spenceri and the non-threatened Litoria lesueurii. We additionally simulated the impact of plausible L. spenceri conservation management interventions on competitive outcomes between these two species. Despite being more susceptible to disease than its competitor, a high relative rate of recruitment allowed the non-threatened L. lesueurii to reach substantially higher densities than L. spenceri in our baseline models, applying a strong absolute force of infection on L. spenceri as an amplifying host. However, simulated management interventions which bolstered L. spenceri recruitment (i.e., captive breeding and release, removal of predatory non-native trout) spurred strong recoveries of L. spenceri while simultaneously (1) increasing the force of Bd infection in the environment and (2) reducing L. lesueurii population density. At high and moderate elevations, combined captive breeding/release and non-native trout removal were sufficient to make L. spenceri the most abundant species. Overall, our results demonstrate the importance of recruitment in moderating pathogen dynamics of multi-host amphibian chytridiomycosis systems. While infection-based parameters are undoubtedly important in Bd management, modifying relative rates of recruitment can substantially alter pathogen-mediated competition between species of an amphibian community.

Naïve and induced tolerance of 15 amphibian populations to three commonly applied insecticides

Human impacts on ecological communities are pervasive and species must either move or adapt to changing environmental conditions. For environments polluted by contaminants, researchers have found hundreds of target pest species evolving increased tolerance, but we have substantially fewer cases of evolved tolerance in non-target species. When species do evolve increased tolerance, inducible tolerance can provide immediate protection and favor the evolution of increased tolerance over generations via genetic assimilation. Using a model larval amphibian (wood frogs, Rana sylvatica), we examined the tolerance of 15 populations from western Pennsylvania and eastern New York (USA), when first exposed to no pesticide or sublethal concentrations and subsequently exposed to lethal concentrations of three common insecticides (carbaryl, chlorpyrifos, and diazinon). We found high variation in naïve tolerance among the populations for all three insecticides. We also discovered that nearly half of the populations exhibited inducible tolerance, though the degree of inducible tolerance (magnitude of tolerance plasticity; MoTP) varied. We observed a cross-tolerance pattern of the populations between chlorpyrifos and diazinon, but no pattern of similar MoTP among the pesticides. With populations combined from two regions, increased tolerance was not associated with proximity to agricultural fields, but there were correlations between proximity to agriculture and MoTP. Collectively, these results suggests that amphibian populations possess a wide range of naïve tolerance to common pesticides, with many also being able to rapidly induce increased tolerance. Future research should examine inducible tolerance in a wide variety of other taxa and contaminants to determine the ubiquity of these responses to anthropogenic factors.

Batrachochytrium dendrobatidis strain affects transcriptomic response in liver but not skin in latitudinal populations of the common toad (Bufo bufo)

Batrachochytrium dendrobatidis (Bd) is a fungal pathogen that has decimated amphibian populations worldwide for several decades. We examined the changes in gene expression in response to Bd infection in two populations of the common toad, Bufo bufo, in a laboratory experiment. We collected B. bufo eggs in southern and northern Sweden, and infected the laboratory-raised metamorphs with two strains of the global panzoonotic lineage Bd-GPL. Differential expression analysis showed significant differences between infected and control individuals in both liver and skin. The skin samples showed no discernible differences in gene expression between the two strains used, while liver samples were differentiated by strain, with one of the strains eliciting no immune response from infected toads. Immune system genes were overexpressed in skin samples from surviving infected individuals, while in liver samples the pattern was more diffuse. Splitting samples by population revealed a stronger immune response in northern individuals. Differences in transcriptional regulation between populations are particularly relevant to study in Swedish amphibians, which may have experienced varying exposure to Bd. Earlier exposure to this pathogen and subsequent adaptation or selection pressure may contribute to the survival of some populations over others, while standing genetic diversity in different populations may also affect the infection outcome.

Body size mediates latitudinal population differences in the response to chytrid fungus infection in two amphibians

Factors behind intraspecific variation in sensitivity to pathogens remain poorly understood. We investigated how geographical origin in two North European amphibians affects tolerance to infection by the chytrid fungus Batrachochytrium dendrobatidis (Bd), a generalist pathogen which has caused amphibian population declines worldwide. We exposed newly metamorphosed individuals of moor frog Rana arvalis and common toad Bufo bufo from two latitudinal regions to two different BdGPL strains. We measured survival and growth as infections may cause sub-lethal effects in fitness components even in the absence of mortality. Infection loads were higher in B. bufo than in R. arvalis, and smaller individuals had generally higher infection loads. B. bufo had high mortality in response to Bd infection, whereas there was little mortality in R. arvalis. Bd-mediated mortality was size-dependent and high-latitude individuals were smaller leading to high mortality in the northern B. bufo. Bd exposure led to sub-lethal effects in terms of reduced growth suggesting that individuals surviving the infection may have reduced fitness mediated by smaller body size. In both host species, the Swedish Bd strain caused stronger sublethal effects than the British strain. We suggest that high-latitude populations can be more vulnerable to chytrids than those from lower latitudes and discuss the possible mechanisms how body size and host geographical origin contribute to the present results.

Comparative Analysis of Skin Bacterial Diversity and Its Potential Antifungal Function Between Desert and Pine Forest Populations of Boreal Toads Anaxyrus boreas

The skin microbiome in amphibians has gained a lot of attention as some of its members play a protective role against pathogens such as the fungus Batrachochytrium dendrobatidis (Bd). The composition of skin bacterial communities has been suggested as one of the factors explaining differences in susceptibility to Bd among amphibian species and populations. The boreal toad Anaxyrus boreas is known to be susceptible to Bd, and severe population declines in its southeastern range have been documented. However, throughout A. boreas distribution, populations present differences in susceptibility to Bd infections which may be associated with differences in skin microbial diversity. This study compared the skin bacterial diversity and Bd infection levels of A. boreas in one desert population and one pine forest population from Baja California, Mexico. We found that desert and pine forest toad populations exhibit differences in skin bacterial community structure but show similar Bd infection levels. Using a predictive method, we found that the abundance of bacteria with potential Bd-inhibitory properties differed between uninfected and infected individuals but not between populations. Our data suggest that several bacteria in the skin community may be offering protection from Bd infections in these A. boreas populations. This study provides foundational evidence for future studies seeking to understand the skin-microbial variation among boreal toads' populations and its relation with Bd susceptibility.

Chloride Toxicity to Native Freshwater Species in Natural and Reconstituted Prairie Pothole Waters

Oil and gas extraction in the Prairie Pothole Region (PPR) of the northern USA has resulted in elevated chloride concentrations in ground and surface water due to widespread contamination with highly saline produced water, or brine. The toxicity of chloride is poorly understood in the high hardness waters characteristic of the region. We evaluated the toxicity of chloride to two endemic species, Daphnia magna (water flea) and Lemna gibba (duckweed), exposed in field-collected waters (hardness ~ 3000 mg/L as CaCO₃) and reconstituted waters (hardness 370 mg/L as CaCO₃) intended to mimic PPR background waters. We also investigated the role of chloride in the toxicity of water reconstituted to mimic legacy brine-contaminated wetlands, using two populations of native Pseudacris maculata (Boreal Chorus Frog). Chloride toxicity was similar in field-collected and reconstituted waters for both D. magna (LC50s 3070-3788 mg Cl^-1/L) and L. gibba (IC50s 2441-2887). Although hardness can ameliorate chloride toxicity at low to high hardness, we did not observe additional protection as hardness increased from 370 to ~ 3000 mg/L. In P. maculata exposures, chloride did not fully explain toxicity. Chloride sensitivity also differed between populations, with mortality at 2000 mg Cl^-/L in one population but not the other, and population-specific growth responses. Overall, these results (1) document toxicity to native species at chloride concentrations occurring in the PPR, (2) indicate that very high hardness in the region's waters may not provide additional protection against chloride and (3) highlight challenges of brine investigations, including whether surrogate study populations are representative of local populations.

Challenging a host-pathogen paradigm: Susceptibility to chytridiomycosis is decoupled from genetic erosion

The putatively positive association between host genetic diversity and the ability to defend against pathogens has long attracted the attention of evolutionary biologists. Chytridiomycosis, a disease caused by the chytrid fungus Batrachochytrium dendrobatidis (Bd), has emerged in recent decades as a cause of dramatic declines and extinctions across the amphibian clade. Bd susceptibility can vary widely across populations of the same species, but the relationship between standing genetic diversity and susceptibility has remained notably underexplored so far. Here, we focus on a putatively Bd-naive system of two mainland and two island populations of the common toad (Bufo bufo) at the edge of the species' range and use controlled infection experiments and dd-RAD sequencing of >10 000 SNPs across 95 individuals to characterize the role of host population identity, genetic variation and individual body mass in mediating host response to the pathogen. We found strong genetic differentiation between populations and marked variation in their susceptibility to Bd. This variation was not, however, governed by isolation-mediated genetic erosion, and individual heterozygosity was even found to be negatively correlated with survival. Individual survival during infection experiments was strongly positively related to body mass, which itself was unrelated to population of origin or heterozygosity. Our findings underscore the general importance of context-dependency when assessing the role of host genetic variation for the ability of defence against pathogens.

Direct and Latent Effects of Pathogen Exposure Across Native and Invasive Amphibian Life Stages

Emerging infectious diseases are one of the multiple factors contributing to the current "biodiversity crisis". As part of the worldwide biodiversity crisis, amphibian populations are declining globally. Chytridiomycosis, an emerging infectious disease, caused by the fungal pathogen Batrachochytrium dendrobatidis (Bd), is a major cause of amphibian population declines. This fungus primarily affects keratinized structures in larval, juvenile, and adult amphibians as well as heart function. However, we know little about how Bd can impact embryos as well as potential latent effects of Bd exposure over ontogeny. Using two different Bd strains and multiple exposure times, we examined the effects of Bd exposure in Pacific chorus frog (Pseudacris regilla), Western toad (Anaxyrus boreas) and American bullfrog (Lithobates catesbeianus) life stages. Using a factorial experimental design, embryos of these three species were exposed to Bd at early and late embryonic stages, with some individuals re-exposed after hatching. Embryonic Bd exposure resulted in differential survival as a function of host species, Bd strain and timing of exposure. P. regilla experienced embryonic mortality when exposed during later developmental stages to one Bd strain. There were no differences across the treatments in embryonic mortality of A. boreas and embryonic mortality of L. catesbeianus occurred in all Bd exposure treatments. We detected latent effects in A. boreas and L. catesbeianus larvae, as mortality increased when individuals had been exposed to any of the Bd strains during the embryonic stage. We also detected direct effects on larval mortality in all three anuran species as a function of Bd strain, and when individuals were double exposed (late in the embryonic stage and again as larvae). Our results suggest that exposure to Bd can directly affect embryo survival and has direct and latent effects on larvae survival of both native and invasive species. However, these impacts were highly context dependent, with timing of exposure and Bd strain influencing the severity of the effects.

Variability in environmental persistence but not per capita transmission rates of the amphibian chytrid fungus leads to differences in host infection prevalence

Heterogeneities in infections among host populations may arise through differences in environmental conditions through two mechanisms. First, environmental conditions may alter host exposure to pathogens via effects on survival. Second, environmental conditions may alter host susceptibility, making infection more or less likely if contact between a host and pathogen occurs. Further, host susceptibility might be altered through acquired resistance, which hosts can develop, in some systems, through exposure to dead or decaying pathogens and their metabolites. Environmental conditions may alter the rates of pathogen decomposition, influencing the likelihood of hosts developing acquired resistance. The present study primarily tests how environmental context influences the relative contributions of pathogen survival and per capita transmission on host infection prevalence using the amphibian chytrid fungus (Batrachochytrium dendrobatidis; Bd) as a model system. Secondarily, we evaluate how environmental context influences the decomposition of Bd because previous studies have shown that dead Bd and its metabolites can illicit acquired resistance in hosts. We conducted Bd survival and infection experiments and then fit models to discern how Bd mortality, decomposition and per capita transmission rates vary among water sources [e.g. artificial spring water (ASW) or water from three ponds]. We found that infection prevalence differed among water sources, which was driven by differences in mortality rates of Bd, rather than differences in per capita transmission rates. Bd mortality rates varied among pond water treatments and were lower in ASW compared to pond water. These results suggest that variation in Bd infection dynamics could be a function of environmental factors in waterbodies that result in differences in exposure of hosts to live Bd. In contrast to the persistence of live Bd, we found that the rates of decomposition of dead Bd did not vary among water sources, which may suggest that exposure of hosts to dead Bd or its metabolites might not commonly vary among nearby sites. Ultimately, a mechanistic understanding of the environmental dependence of free-living pathogens could lead to a deeper understanding of the patterns of outbreak heterogeneity, which could inform surveillance and management strategies.

Effects of invasive larval bullfrogs (Rana catesbeiana) on disease transmission, growth and survival in the larvae of native amphibians

The mechanisms by which invasive species negatively affect native species include competition, predation, and the introduction of novel pathogens. Moreover, if an invasive species is a competent disease reservoir, it may facilitate the long-term maintenance and spread of pathogens in ecological assemblages and drive the extinction of less tolerant or less resistant species. Disease-driven loss of biodiversity is exemplified by the amphibian–chytrid fungus system. The disease chytridiomycosis is caused by the aquatic chytrid fungus Batrachochytrium dendrobatidis (Bd) in anurans and is associated with worldwide amphibian population declines and extinctions. For amphibian species that metamorphose and leave infected aquatic habitats, the mechanisms by which Bd persists over winter in these habitats remains a critical open question. A leading hypothesis is that American bullfrogs (Rana catesbeiana), a worldwide invasive species, are tolerant to Bd and serve as a reservoir host for Bd during winter months and subsequently infect native species that return to breed in spring. Using outdoor mesocosms, we experimentally examined if two strains of Bd could overwinter in aquatic systems, in the presence or absence of bullfrog tadpoles, and if overwintered Bd could be transmitted to tadpoles of two spring-breeding species: Pacific treefrogs (Pseudacris regilla) and Cascades frogs (Rana cascadae). We found that only 4 of 448 total animals (one bullfrog and three spring breeders) tested positive for Bd after overwintering. Moreover, two of the three infected spring breeders emerged from tanks that contained overwintered Bd but in the absence of infected bullfrogs. This suggests that Bd can persist over winter without bullfrogs as a reservoir host. We found no effect of Bd strain on bullfrog survival after overwintering. For Pacific treefrogs, Bd exposure did not significantly affect mass at or time to metamorphosis while exposure to bullfrogs reduced survival. For Cascades frogs, we found an interactive effect of Bd strain and bullfrog presence on time to metamorphosis, but no main or interactive effects on their survival or mass at metamorphosis. In short, bullfrog tadpoles rarely retained and transmitted Bd infection in our experiment and we found limited evidence that Bd successfully overwinters in the absence of bullfrog tadpoles and infects spring-breeding amphibians.

Guidance for Developing Amphibian Population Models for Ecological Risk Assessment

Despite widespread acceptance of the utility of population modeling and advocacy of this approach for a more ecologically relevant perspective, it is not routinely incorporated in ecological risk assessments (ERA). A systematic framework for situation-specific model development is one of the major challenges to broadly adopting population models in ERA. As risk assessors confront the multitude of species and chemicals requiring evaluation, an adaptable stepwise guide for model parameterization would facilitate this process. Additional guidance on interpretation of model output and evaluating uncertainty would further contribute to establishing consensus on good modeling practices. We build on previous work that created a framework and decision guide for developing population models for ERA by focusing on data types, model structure, and extrinsic stressors relevant to anuran amphibians. Anurans have a unique life cycle with varying habitat requirements and high phenotypic plasticity. These species belong to the amphibian class, which is facing global population decline in large part due to anthropogenic stressors, including chemicals. We synthesize information from databases and literature relevant to amphibian risks to identify traits that influence exposure likelihood, inherent sensitivity, population vulnerability, and environmental constraints. We link these concerns with relevant population modeling methods and structure in order to evaluate pesticide effects with appropriate scale and parameterization. A standardized population modeling approach, with additional guidance for anuran ERA, offers an example method for quantifying population risks and evaluating long-term impacts of chemical stressors to populations. Integr Environ Assess Manag 2020;16:223-233. © 2019 SETAC.

Relationships between glucocorticoids and infection with Batrachochytrium dendrobatidis in three amphibian species

It is often hypothesized that organisms exposed to environmental change may experience physiological stress, which could reduce individual quality and make them more susceptible to disease. Amphibians are amongst the most threatened taxa, particularly in the context of disease, but relatively few studies explore links between stress and disease in amphibian species. Here, we use the fungal pathogen Batrachochytrium dendrobatidis (Bd) and amphibians as an example to explore relationships between disease and glucocorticoids (GCs), metabolic hormones that comprise one important component of the stress response. While previous work is limited, it has largely identified positive relationships between GCs and Bd-infection. However, the causality remains unclear and few studies have integrated both baseline (GC release that is related to standard, physiological functioning) and stress-induced (GC release in response to an acute stressor) measures of GCs. Here, we examine salivary corticosterone before and after exposure to a stressor, in both field and captive settings. We present results for Bd-infected and uninfected individuals of three amphibian species with differential susceptibilities to this pathogen (Rana catesbeiana, R. clamitans, and R. sylvatica). We hypothesized that prior to stress, baseline GCs would be higher in Bd-infected animals, particularly in more Bd-susceptible species. We also expected that after exposure to a stressor, stress-induced GCs would be lower in Bd-infected animals. These species exhibited significant interspecific differences in baseline and stress induced corticosterone, though other variables like sex, body size, and day of year were usually not predictive of corticosterone. In contrast to most previous work, we found no relationships between Bd and corticosterone for two species (R. catesbeiana and R. clamitans), and in the least Bd-tolerant species (R. sylvatica) animals exhibited context-dependent differences in relationships between Bd infection and corticosterone: Bd-positive R. sylvatica had significantly lower baseline and stress-induced corticosterone, with this pattern being stronger in the field than in captivity. These results were surprising, as past work in other species has more often found elevated GCs in Bd-positive animals, a pattern that aligns with well-documented relationships between chronically high GCs, reduced individual quality, and immunosuppression. This work highlights the potential relevance of GCs to disease susceptibility in the context of amphibian declines, while underscoring the importance of characterizing these relationships in diverse contexts.

Host age alters amphibian susceptibility to Batrachochytrium dendrobatidis, an emerging infectious fungal pathogen

Parasites and pathogens are often aggregated in a minority of susceptible hosts within a population, with a majority of individuals harboring low infection intensities. However, determining the relative importance of host traits to explain this heterogeneity is a challenge. One ecologically important pathogen is Batrachochytrium dendrobatidis (Bd), which causes the disease chytridiomycosis and has been associated with many amphibian population declines worldwide. For many hosts, post-metamorphic stages are generally more susceptible than the larval stage. Yet, examination of the effects of Bd infection at different ages within a life stage, has received little attention. This study investigated the hypothesis that recently-post-metamorphic frogs were more sensitive to chytridiomycosis than older frogs, and that sensitivity to Bd infection decreased as frogs aged. We examined this relationship with Pacific treefrogs (Pseudacris regilla) and red legged frogs (Rana aurora). Age had a strong effect on susceptibility to infection, infection intensity, and survival-but not in the directions we had predicted. In both host species, an increase in age was associated with frogs becoming more susceptible to Bd infection, harboring larger infection intensities, and greater risk of mortality. This suggests that the timing of Bd exposure may influence amphibian population dynamics.

Effects of Emerging Infectious Diseases on Amphibians: A Review of Experimental Studies

Numerous factors are contributing to the loss of biodiversity. These include complex effects of multiple abiotic and biotic stressors that may drive population losses. These losses are especially illustrated by amphibians, whose populations are declining worldwide. The causes of amphibian population declines are multifaceted and context-dependent. One major factor affecting amphibian populations is emerging infectious disease. Several pathogens and their associated diseases are especially significant contributors to amphibian population declines. These include the fungi Batrachochytrium dendrobatidis and B. salamandrivorans, and ranaviruses. In this review, we assess the effects of these three pathogens on amphibian hosts as found through experimental studies. Such studies offer valuable insights to the causal factors underpinning broad patterns reported through observational studies. We summarize key findings from experimental studies in the laboratory, in mesocosms, and from the field. We also summarize experiments that explore the interactive effects of these pathogens with other contributors of amphibian population declines. Though well-designed experimental studies are critical for understanding the impacts of disease, inconsistencies in experimental methodologies limit our ability to form comparisons and conclusions. Studies of the three pathogens we focus on show that host susceptibility varies with such factors as species, host age, life history stage, population and biotic (e.g., presence of competitors, predators) and abiotic conditions (e.g., temperature, presence of contaminants), as well as the strain and dose of the pathogen, to which hosts are exposed. Our findings suggest the importance of implementing standard protocols and reporting for experimental studies of amphibian disease.

Host-pathogen metapopulation dynamics suggest high elevation refugia for boreal toads

Emerging infectious diseases are an increasingly common threat to wildlife. Chytridiomycosis, caused by the fungal pathogen Batrachochytrium dendrobatidis (Bd), is an emerging infectious disease that has been linked to amphibian declines around the world. Few studies exist that explore amphibian-Bd dynamics at the landscape scale, limiting our ability to identify which factors are associated with variation in population susceptibility and to develop effective in situ disease management. Declines of boreal toads (Anaxyrus boreas boreas) in the southern Rocky Mountains are largely attributed to chytridiomycosis but variation exists in local extinction of boreal toads across this metapopulation. Using a large-scale historic data set, we explored several potential factors influencing disease dynamics in the boreal toad-Bd system: geographic isolation of populations, amphibian community richness, elevational differences, and habitat permanence. We found evidence that boreal toad extinction risk was lowest at high elevations where temperatures may be suboptimal for Bd growth and where small boreal toad populations may be below the threshold needed for efficient pathogen transmission. In addition, boreal toads were more likely to recolonize high elevation sites after local extinction, again suggesting that high elevations may provide refuge from disease for boreal toads. We illustrate a modeling framework that will be useful to natural resource managers striving to make decisions in amphibian-Bd systems. Our data suggest that in the southern Rocky Mountains high elevation sites should be prioritized for conservation initiatives like reintroductions.

Gene expression differs in susceptible and resistant amphibians exposed to Batrachochytrium dendrobatidis

Chytridiomycosis, the disease caused by the fungal pathogen Batrachochytrium dendrobatidis (Bd), has devastated global amphibian biodiversity. Nevertheless, some hosts avoid disease after Bd exposure even as others experience near-complete extirpation. It remains unclear whether the amphibian adaptive immune system plays a role in Bd defence. Here, we describe gene expression in two host species-one susceptible to chytridiomycosis and one resistant-following exposure to two Bd isolates that differ in virulence. Susceptible wood frogs (Rana sylvatica) had high infection loads and mortality when exposed to the more virulent Bd isolate but lower infection loads and no fatal disease when exposed to the less virulent isolate. Resistant American bullfrogs (R. catesbeiana) had high survival across treatments and rapidly cleared Bd infection or avoided infection entirely. We found widespread upregulation of adaptive immune genes and downregulation of important metabolic and cellular maintenance components in wood frogs after Bd exposure, whereas American bullfrogs showed little gene expression change and no evidence of an adaptive immune response. Wood frog responses suggest that adaptive immune defences may be ineffective against virulent Bd isolates that can cause rapid physiological dysfunction. By contrast, American bullfrogs exhibited robust resistance to Bd that is likely attributable, at least in part, to their continued upkeep of metabolic and skin integrity pathways as well as greater antimicrobial peptide expression compared to wood frogs, regardless of exposure. Greater understanding of these defences will ultimately help conservationists manage chytridiomycosis.

Effects of amphibian phylogeny, climate and human impact on the occurrence of the amphibian-killing chytrid fungus

Chytridiomycosis, due to the fungus Batrachochytrium dendrobatidis (Bd), has been associated with the alarming decline and extinction crisis of amphibians worldwide. Because conservation programs are implemented locally, it is essential to understand how the complex interactions among host species, climate and human activities contribute to Bd occurrence at regional scales. Using weighted phylogenetic regressions and model selection, we investigated geographic patterns of Bd occurrence along a latitudinal gradient of 1500 km within a biodiversity hot spot in Chile (1845 individuals sampled from 253 sites and representing 24 species), and its association with climatic, socio-demographic and economic variables. Analyses show that Bd prevalence decreases with latitude although it has increased by almost 10% between 2008 and 2013, possibly reflecting an ongoing spread of Bd following the introduction of Xenopus laevis. Occurrence of Bd was higher in regions with high gross domestic product (particularly near developed centers) and with a high variability in rainfall regimes, whereas models including other bioclimatic or geographic variables, including temperature, exhibited substantially lower fit and virtually no support based on Akaike weights. In addition, Bd prevalence exhibited a strong phylogenetic signal, with five species having high numbers of infected individuals and higher prevalence than the average of 13.3% across all species. Taken together, our results highlight that Bd in Chile might still be spreading south, facilitated by a subset of species that seem to play an important epidemiological role maintaining this pathogen in the communities, in combination with climatic and human factors affecting the availability and quality of amphibian breeding sites. This information may be employed to design conservation strategies and mitigate the impacts of Bd in the biodiversity hot spot of southern Chile, and similar studies may prove useful to disentangle the role of different factors contributing to the emergence and spread of this catastrophic disease.

Host niche may determine disease-driven extinction risk

The fungal pathogen Batrachochytrium dendrobatidis (Bd) drives declines and extinctions in amphibian communities. However, not all regions and species are equally affected. Here, we show that association with amphibian aquatic habitat types (bromeliad phytotelmata versus stream) across Central America results in the odds of being threatened by Bd being five times higher in stream microhabitats. This differential threat of Bd was supported in our study by a significantly lower prevalence of Bd in bromeliad-associated amphibian species compared to riparian species in Honduran cloud forests. Evidence that the bromeliad environment is less favorable for Bd transmission is exemplified by significantly less suitable physicochemical conditions and higher abundance of Bd-ingesting micro-eukaryotes present in bromeliad water. These factors may inhibit aquatic Bd zoospore survival and the development of an environmental reservoir of the pathogen. Bromeliad phytotelmata thus may act as environmental refuges from Bd, which contribute to protecting associated amphibian communities against chytridiomycosis-driven amphibian declines that threaten the nearby riparian communities.

Evolved pesticide tolerance influences susceptibility to parasites in amphibians

Because ecosystems throughout the globe are contaminated with pesticides, there is a need to understand how natural populations cope with pesticides and the implications for ecological interactions. From an evolutionary perspective, there is evidence that pesticide tolerance can be achieved via two mechanisms: selection for constitutive tolerance over multiple generations or by inducing tolerance within a single generation via phenotypic plasticity. While both mechanisms can allow organisms to persist in contaminated environments, they might result in different performance trade-offs including population susceptibility to parasites. We have identified 15 wood frog populations that exist along a gradient from close to agriculture and high, constitutive pesticide tolerance to far from agriculture and inducible pesticide tolerance. Using these populations, we investigated the relationship between evolutionary responses to the common insecticide carbaryl and host susceptibility to the trematode Echinoparyphium lineage 3 and ranavirus using laboratory exposure assays. For Echinoparyphium, we discovered that wood frog populations living closer to agriculture with high, constitutive tolerance experienced lower loads than populations living far from agriculture with inducible pesticide tolerance. For ranavirus, we found no relationship between the mechanism of evolved pesticide tolerance and survival, but populations living closer to agriculture with high, constitutive tolerance experienced higher viral loads than populations far from agriculture with inducible tolerance. Land use and mechanisms of evolved pesticide tolerance were associated with susceptibility to parasites, but the direction of the relationship is dependent on the type of parasite, underscoring the complexity between land use and disease outcomes. Collectively, our results demonstrate that evolved pesticide tolerance can indirectly influence host-parasite interactions and underscores the importance of including evolutionary processes in ecotoxicological studies.

Virulence variation among strains of the emerging infectious fungus Batrachochytrium dendrobatidis (Bd) in multiple amphibian host species

Emerging infectious diseases have been documented in numerous plant and animal populations. The infectious disease amphibian chytridiomycosis, caused by the fungus Batrachochytrium dendrobatidis (Bd), is associated with global amphibian population declines. While much Bd-amphibian research has centered on response variation in hosts, a paucity of information exists on how variation in the pathogen, such as strain differences, affects infection dynamics. To examine how different Bd strains may differentially impact multiple hosts, we conducted laboratory experiments to measure 2 infection outcomes, viz. host survival and pathogen load, in 3 amphibian host species (Pacific treefrog, western toad, and Cascades frog) after exposure to 3 different Bd strains (an additional fourth Bd strain was tested in toads only). Our results confirm that the infection response differs among host species. Western toads experienced significant mortality, but Pacific treefrogs and Cascades frogs did not. Interestingly, our experiment also captured strain-dependent virulence variation but only in 1 host species, the western toad. Increased mortality was observed in 2 of the 4 Bd strains tested in this host species. Toads were also the only host species found to have variable pathogen load dependent on strain type; individuals exposed to the Panama strain harbored significantly higher loads compared to all other strains. These findings underscore the dynamic nature of Bd infection, showing that virulence can vary contingent on host and strain type. We highlight the importance of both host- and pathogen-dependent factors in determining overall infection virulence and show the need for in vivo testing to fully assess pathogenicity.

Complex interactive effects of water mold, herbicide, and the fungus Batrachochytrium dendrobatidis on Pacific treefrog Hyliola regilla hosts

Infectious diseases pose a serious threat to global biodiversity. However, their ecological impacts are not independent of environmental conditions. For example, the pathogenic fungus Batrachochytrium dendrobatidis (Bd), which has contributed to population declines and extinctions in many amphibian species, interacts with several environmental factors to influence its hosts, but potential interactions with other pathogens and environmental contaminants are understudied. We examined the combined effects of Bd, a water mold (Achlya sp.), and the herbicide Roundup® Regular (hereafter, Roundup®) on larval Pacific treefrog Hyliola regilla hosts. We employed a 2 wk, fully factorial laboratory experiment with 3 ecologically realistic levels (0, 1, and 2 mg l-1 of active ingredient) of field-formulated Roundup®, 2 Achlya treatments (present and absent), and 2 Bd treatments (present and absent). Our results were consistent with sublethal interactive effects involving all 3 experimental factors. When Roundup® was absent, the proportion of Bd-exposed larvae infected with Bd was elevated in the presence of Achlya, consistent with Achlya acting as a synergistic cofactor that facilitated the establishment of Bd infection. However, this Achlya effect became nonsignificant at 1 mg l-1 of the active ingredient of Roundup® and disappeared at the highest Roundup® concentration. In addition, Roundup® decreased Bd loads among Bd-exposed larvae. Our study suggests complex interactive effects of a water mold and a contaminant on Bd infection in amphibian hosts. Achlya and Roundup® were both correlated with altered patterns of Bd infection, but in different ways, and Roundup® appeared to remove the influence of Achlya on Bd.

Effect of Simultaneous Amphibian Exposure to Pesticides and an Emerging Fungal Pathogen, Batrachochytrium dendrobatidis

Amphibian declines have been linked to numerous factors, including pesticide use and the fungal pathogen Batrachochytrium dendrobatidis (Bd). Moreover, research has suggested a link between amphibian sensitivity to Bd and pesticide exposure. We simultaneously exposed postmetamorphic American toads (Anaxyrus americanus), western toads (A. boreas), spring peepers (Pseudacris crucifer), Pacific treefrogs (P. regilla), leopard frogs (Lithobates pipiens), and Cascades frogs (Rana cascadae) to a factorial combination of two pathogen treatments (Bd⁺, Bd^-) and four pesticide treatments (control, ethanol vehicle, herbicide mixture, and insecticide mixture) for 14 d to quantify survival and infection load. We found no interactive effects of pesticides and Bd on anuran survival and no effects of pesticides on infection load. Mortality following Bd exposure increased in spring peepers and American toads and was dependent upon snout-vent length in western toads, American toads, and Pacific treefrogs. Previous studies reported effects of early sublethal pesticide exposure on amphibian Bd sensitivity and infection load at later life stages, but we found simultaneous exposure to sublethal pesticide concentrations and Bd had no such effect on postmetamorphic juvenile anurans. Future research investigating complex interactions between pesticides and Bd should employ a variety of pesticide formulations and Bd strains and follow the effects of exposure throughout ontogeny.

Projecting the Global Distribution of the Emerging Amphibian Fungal Pathogen, Batrachochytrium dendrobatidis, Based on IPCC Climate Futures

Projected changes in climate conditions are emerging as significant risk factors to numerous species, affecting habitat conditions and community interactions. Projections suggest species range shifts in response to climate change modifying environmental suitability and is supported by observational evidence. Both pathogens and their hosts can shift ranges with climate change. We consider how climate change may influence the distribution of the emerging infectious amphibian chytrid fungus, Batrachochytrium dendrobatidis (Bd), a pathogen associated with worldwide amphibian population losses. Using an expanded global Bd database and a novel modeling approach, we examined a broad set of climate metrics to model the Bd-climate niche globally and regionally, then project how climate change may influence Bd distributions. Previous research showed that Bd distribution is dependent on climatic variables, in particular temperature. We trained a machine-learning model (random forest) with the most comprehensive global compilation of Bd sampling records (~5,000 site-level records, mid-2014 summary), including 13 climatic variables. We projected future Bd environmental suitability under IPCC scenarios. The learning model was trained with combined worldwide data (non-region specific) and also separately per region (region-specific). One goal of our study was to estimate of how Bd spatial risks may change under climate change based on the best available data. Our models supported differences in Bd-climate relationships among geographic regions. We projected that Bd ranges will shift into higher latitudes and altitudes due to increased environmental suitability in those regions under predicted climate change. Specifically, our model showed a broad expansion of areas environmentally suitable for establishment of Bd on amphibian hosts in the temperate zones of the Northern Hemisphere. Our projections are useful for the development of monitoring designs in these areas, especially for sensitive species and those vulnerable to multiple threats.

Leaf Litter Inhibits Growth of an Amphibian Fungal Pathogen

Past studies have found a heterogeneous distribution of the amphibian chytrid fungal pathogen, Batrachochytrium dendrobatidis (Bd). Recent studies have accounted for some of this heterogeneity through a positive association between canopy cover and Bd abundance, which is attributed to the cooling effect of canopy cover. We questioned whether leaf litter inputs that are also associated with canopy cover might also alter Bd growth. Leaf litter inputs exhibit tremendous interspecific chemical variation, and we hypothesized that Bd growth varies with leachate chemistry. We also hypothesized that Bd uses leaf litter as a growth substrate. To test these hypotheses, we conducted laboratory trials in which we exposed cultures of Bd to leachate of 12 temperate leaf litter species at varying dilutions. Using a subset of those 12 litter species, we also exposed Bd to pre-leached litter substrate. We found that exposure to litter leachate and substrate reduced Bd spore and sporangia densities, although there was substantial variation among treatments. In particular, Bd densities were inversely correlated with concentrations of phenolic acids. We conducted a field survey of phenolic concentrations in natural wetlands which verified that the leachate concentrations in our lab study are ecologically relevant. Our study reinforces prior indications that positive associations between canopy cover and Bd abundance are likely mediated by water temperature effects, but this phenomenon might be counteracted by changes in aquatic chemistry from leaf litter inputs.

Population-specific toxicity of six insecticides to the trematode Echinoparyphium sp

The ubiquitous use of pesticides has increased concerns over their direct and indirect effects on disease dynamics. While studies examining the effects of pesticides on host-parasite interactions have largely focused on how pesticides influence the host, few studies have considered the effects of pesticides on parasites. We investigated the toxicity of six common insecticides at six environmentally-relevant concentrations to cercariae of the trematode Echinoparyphium from two populations. All six insecticides reduced the survival of cercariae (overall difference between mortality in control vs pesticide exposure = 86·2 ± 8·7%) but not in a predictable dose-dependent manner. These results suggest that Echinoparyphium are sensitive to a broad range of insecticides commonly used in the USA. The lack of a clear dose-dependent response in Echinoparyphium highlights the potential limitations of toxicity assays in predicting pesticide toxicity to parasites. Finally, population-level variation in cercarial susceptibility to pesticides underscores the importance of accounting for population variation as overlooking this variation can limit our ability to predict toxicity in nature. Collectively, this work demonstrates that consideration of pesticide toxicity to parasites is important to understanding how pesticides ultimately shape disease dynamics in nature.

Life history linked to immune investment in developing amphibians

The broad diversity of amphibian developmental strategies has been shaped, in part, by pathogen pressure, yet trade-offs between the rate of larval development and immune investment remain poorly understood. The expression of antimicrobial peptides (AMPs) in skin secretions is a crucial defense against emerging amphibian pathogens and can also indirectly affect host defense by influencing the composition of skin microbiota. We examined the constitutive or induced expression of AMPs in 17 species at multiple life-history stages. We found that AMP defenses in tadpoles of species with short larval periods (fast pace of life) were reduced in comparison with species that overwinter as tadpoles and grow to a large size. A complete set of defensive peptides emerged soon after metamorphosis. These findings support the hypothesis that species with a slow pace of life invest energy in AMP production to resist potential pathogens encountered during the long larval period, whereas species with a fast pace of life trade this investment in defense for more rapid growth and development.