Experimental evolution alters the rate and temporal pattern of population growth in Batrachochytrium dendrobatidis, a lethal fungal pathogen of amphibians

Coinfection with chytrid genotypes drives divergent infection dynamics reflecting regional distribution patterns

By altering the abundance, diversity, and distribution of species-and their pathogens-globalization may inadvertently select for more virulent pathogens. In Brazil's Atlantic Forest, a hotspot of amphibian biodiversity, the global amphibian trade has facilitated the co-occurrence of previously isolated enzootic and panzootic lineages of the pathogenic amphibian-chytrid (Batrachochytrium dendrobatidis, 'Bd') and generated new virulent recombinant genotypes ('hybrids'). Epidemiological data indicate that amphibian declines are most severe in hybrid zones, suggesting that coinfections are causing more severe infections or selecting for higher virulence. We investigated how coinfections involving these genotypes shapes virulence and transmission. Overall, coinfection favored the more virulent and competitively superior panzootic genotype, despite dampening its transmission potential and overall virulence. However, for the least virulent and least competitive genotype, coinfection increased both overall virulence and transmission. Thus, by integrating experimental and epidemiological data, our results provide mechanistic insight into how globalization can select for, and propel, the emergence of introduced hypervirulent lineages, such as the globally distributed panzootic lineage of Bd.

Predicted Batrachochytrium dendrobatidis infection sites in Guyana, Suriname, and French Guiana using the species distribution model maxent

The fungal pathogen Batrachochytrium dendrobatidis (Bd) which causes that amphibian disease chytridiomycosis is expanding its worldwide range from an Asian origin, infecting amphibians in a growing number of countries. Modelling the potential range of this amphibian pathogen using environmental variables and presence data could advance our understanding of at-risk areas and species in locations with limited surveillance to date. We used a species distribution model to assess Bd habitat suitability in the three Guiana's (Guyana, Suriname, and French Guiana) in South America. The model output showed that all three countries have substantial areas where Bd could grow and proliferate, and maximum temperature of the warmest month was the top predictor of suitable Bd habitat, inversely correlated with modeled Bd occurrence. Predicted Bd infection areas in Guyana and French Guiana were large and localized whereas possible sites in Suriname were more scattered throughout the country. The areas projected as potential suitable in Suriname were mostly high elevation regions. These results could help inform efficiencies for development of a proactive monitoring program that could alert managers of novel Bd outbreaks for focused mitigation actions to forestall the spread of this amphibian disease.

Chytridiomycosis Outbreak in a Chilean Giant Frog (Calyptocephalella gayi) Captive Breeding Program: Genomic Characterization and Pathological Findings

Emerging infectious diseases in wildlife are increasingly associated with animal mortality and species declines, but their source and genetic characterization often remains elusive. Amphibian chytridiomycosis, caused by the fungus Batrachochytrium dendrobatidis (Bd), has been associated with catastrophic and well-documented amphibian population declines and extinctions at the global scale. We used histology and whole-genome sequencing to describe the lesions caused by, and the genetic variability of, two Bd isolates obtained from a mass mortality event in a captive population of the threatened Chilean giant frog (Calyptocephalella gayi). This was the first time an association between Bd and high mortality had been detected in this charismatic and declining frog species. Pathological examinations revealed that 30 dead metamorphosed frogs presented agnathia or brachygnathia, a condition that is reported for the first time in association with chytridiomycosis. Phylogenomic analyses revealed that Bd isolates (PA1 and PA2) from captive C. gayi group with other Bd isolates (AVS2, AVS4, and AVS7) forming a single highly supported Chilean Bd clade within the global panzootic lineage of Bd (BdGPL). These findings are important to inform the strengthening of biosecurity measures to prevent the impacts of chytridiomycosis in captive breeding programs elsewhere.

Early presence of Batrachochytrium dendrobatidis in Mexico with a contemporary dominance of the global panzootic lineage

Chytridiomycosis, caused by the fungus Batrachochytrium dendrobatidis (Bd), is a devastating infectious disease of amphibians. Retrospective studies using museum vouchers and genetic samples supported the hypothesis that Bd colonized Mexico from North America and then continued to spread into Central and South America, where it led to dramatic losses in tropical amphibian biodiversity (the epizootic wave hypothesis). While these studies suggest that Bd has been in Mexico since the 1970s, information regarding the historical and contemporary occurrence of different pathogen genetic lineages across the country is limited. In the current study, we investigated the historical and contemporary patterns of Bd in Mexico. We combined the swabbing of historical museum vouchers and sampling of wild amphibians with a custom Bd genotyping assay to assess the presence, prevalence, and genetic diversity of Bd over time in Mexico. We found Bd-positive museum specimens from the late 1800s, far earlier than previous records and well before recent amphibian declines. With Bd genotypes from samples collected between 1975-2019, we observed a contemporary dominance of the global panzootic lineage in Mexico and report four genetic subpopulations and potential for admixture among these populations. The observed genetic variation did not have a clear geographic signature or provide clear support for the epizootic wave hypothesis. These results provide a framework for testing new questions regarding Bd invasions and their temporal relationship to observed amphibian declines in the Americas.

Mechanisms underlying host persistence following amphibian disease emergence determine appropriate management strategies

Emerging infectious diseases have caused many species declines, changes in communities and even extinctions. There are also many species that persist following devastating declines due to disease. The broad mechanisms that enable host persistence following declines include evolution of resistance or tolerance, changes in immunity and behaviour, compensatory recruitment, pathogen attenuation, environmental refugia, density-dependent transmission and changes in community composition. Here we examine the case of chytridiomycosis, the most important wildlife disease of the past century. We review the full breadth of mechanisms allowing host persistence, and synthesise research on host, pathogen, environmental and community factors driving persistence following chytridiomycosis-related declines and overview the current evidence and the information required to support each mechanism. We found that for most species the mechanisms facilitating persistence have not been identified. We illustrate how the mechanisms that drive long-term host population dynamics determine the most effective conservation management strategies. Therefore, understanding mechanisms of host persistence is important because many species continue to be threatened by disease, some of which will require intervention. The conceptual framework we describe is broadly applicable to other novel disease systems.

Experimental methodologies can affect pathogenicity of Batrachochytrium salamandrivorans infections

Controlled experiments are one approach to understanding the pathogenicity of etiologic agents to susceptible hosts. The recently discovered fungal pathogen, Batrachochytrium salamandrivorans (Bsal), has resulted in a surge of experimental investigations because of its potential to impact global salamander biodiversity. However, variation in experimental methodologies could thwart knowledge advancement by introducing confounding factors that make comparisons difficult among studies. Thus, our objective was to evaluate if variation in experimental methods changed inferences made on the pathogenicity of Bsal. We tested whether passage duration of Bsal culture, exposure method of the host to Bsal (water bath vs. skin inoculation), Bsal culturing method (liquid vs. plated), host husbandry conditions (aquatic vs. terrestrial), and skin swabbing frequency influenced diseased-induced mortality in a susceptible host species, the eastern newt (Notophthalmus viridescens). We found that disease-induced mortality was faster for eastern newts when exposed to a low passage isolate, when newts were housed in terrestrial environments, and if exposure to zoospores occurred via water bath. We did not detect differences in disease-induced mortality between culturing methods or swabbing frequencies. Our results illustrate the need to standardize methods among Bsal experiments. We provide suggestions for future Bsal experiments in the context of hypothesis testing and discuss the ecological implications of our results.

Isolation and maintenance of Batrachochytrium salamandrivorans cultures

Discovered in 2013, the chytrid fungus Batrachochytrium salamandrivorans (Bsal) is an emerging amphibian pathogen that causes ulcerative skin lesions and multifocal erosion. A closely related pathogen, B. dendrobatidis (Bd), has devastated amphibian populations worldwide, suggesting that Bsal poses a significant threat to global salamander biodiversity. To expedite research into this emerging threat, we seek to standardize protocols across the field so that results of laboratory studies are reproducible and comparable. We have collated data and experience from multiple labs to standardize culturing practices of Bsal. Here we outline common culture practices including a medium for standardized Bsal growth, standard culturing protocols, and a method for isolating Bsal from infected tissue.

Chytrid fungi and global amphibian declines

Discovering that chytrid fungi cause chytridiomycosis in amphibians represented a paradigm shift in our understanding of how emerging infectious diseases contribute to global patterns of biodiversity loss. In this Review we describe how the use of multidisciplinary biological approaches has been essential to pinpointing the origins of amphibian-parasitizing chytrid fungi, including Batrachochytrium dendrobatidis and Batrachochytrium salamandrivorans, as well as to timing their emergence, tracking their cycles of expansion and identifying the core mechanisms that underpin their pathogenicity. We discuss the development of the experimental methods and bioinformatics toolkits that have provided a fuller understanding of batrachochytrid biology and informed policy and control measures.

Eco-evolutionary rescue promotes host-pathogen coexistence

Emerging infectious pathogens are responsible for some of the most severe host mass mortality events in wild populations. Yet, effective pathogen control strategies are notoriously difficult to identify, in part because quantifying and forecasting pathogen spread and disease dynamics is challenging. Following an outbreak, hosts must cope with the presence of the pathogen, leading to host-pathogen coexistence or extirpation. Despite decades of research, little is known about host-pathogen coexistence post-outbreak when low host abundances and cryptic species make these interactions difficult to study. Using a novel disease-structured N-mixture model, we evaluate empirical support for three host-pathogen coexistence hypotheses (source-sink, eco-evolutionary rescue, and spatial variation in pathogen transmission) in a Neotropical amphibian community decimated by Batrachochytrium dendrobatidis (Bd) in 2004. During 2010-2014, we surveyed amphibians in Parque Nacional G. D. Omar Torríjos Herrera, Coclé Province, El Copé, Panama. We found that the primary driver of host-pathogen coexistence was eco-evolutionary rescue, as evidenced by similar amphibian survival and recruitment rates between infected and uninfected hosts. Average apparent monthly survival rates of uninfected and infected hosts were both close to 96%, and the expected number of uninfected and infected hosts recruited (via immigration/reproduction) was less than one host per disease state per 20-m site. The secondary driver of host-pathogen coexistence was spatial variation in pathogen transmission as we found that transmission was highest in areas of low abundance but there was no support for the source-sink hypothesis. Our results indicate that changes in the host community (i.e., through genetic or species composition) can reduce the impacts of emerging infectious disease post-outbreak. Our disease-structured N-mixture model represents a valuable advancement for conservation managers trying to understand underlying host-pathogen interactions and provides new opportunities to study disease dynamics in remnant host populations decimated by virulent pathogens.

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.

Distribution modeling and lineage diversity of the chytrid fungus Batrachochytrium dendrobatidis (Bd) in a central African amphibian hotspot

The amphibian disease chytridiomycosis in amphibians is caused by the chytrid fungus Batrachochytrium dendrobatidis (Bd) and has resulted in dramatic declines and extinctions of amphibian populations worldwide. A hypervirulent, globally-dispersed pandemic lineage (Bd-GPL) is thought to be largely responsible for population declines and extinctions, although numerous endemic lineages have also been found. Recent reports of amphibian declines have been linked to the emergence of Bd in Cameroon, a major hotspot of African amphibian diversity. However, it is not known whether Bd-GPL or other lineages have been found in this region. This study therefore aims to examine Bd lineage diversity in the region and predict the distribution of this pathogen under current and future climate conditions using data from this study and from historical records. Almost 15% (52/360) of individuals tested positive for Bd using a standard quantitative PCR diagnostic. Infected amphibians were found at all eight sites sampled in this study. Species distribution models generated in BIOMOD2 indicate that areas with highest predicted environmental suitability occur in the Cameroon highlands and several protected areas throughout the country. These areas of high environmental suitability for Bd are projected to shift or decrease in size under future climate change. However, montane regions with high amphibian diversity are predicted to remain highly suitable. Phylogenetic analysis of the ITS sequences obtained from a set of positive Bd samples indicate that most fall within the Bd-GPL lineage while the remainder group with isolates from either Brazil or South Korea. Although more in depth phylogenetic analyses are needed, identification of Bd-GPL lineages in areas of high amphibian diversity emphasizes the need to continue to monitor for Bd and develop appropriate conservation strategies to prevent its further spread.

Genomic epidemiology of the emerging pathogen Batrachochytrium dendrobatidis from native and invasive amphibian species in Chile

Emerging fungal diseases represent a threat to food security, animal and human health worldwide. Amphibian chytridiomycosis, caused by the fungus Batrachochytrium dendrobatidis (Bd), has been associated with catastrophic and well-documented amphibian population declines and extinctions. For the first time, Bd was cultured from native and non-native wild amphibians in Chile. Phylogenomic analyses revealed that Chilean isolates AVS2, AVS4 and AVS7 group within the global panzootic lineage of Bd (BdGPL) in a single highly supported clade that includes a genotype previously isolated from the United Kingdom. Our results extend the known distribution of BdGPL in South America and suggest a single and relatively recent introduction of BdGPL into the country, providing additional support to the role of anthropogenic activity in the global spread of this panzootic lineage.

Resistance, tolerance and environmental transmission dynamics determine host extinction risk in a load-dependent amphibian disease

While disease-induced extinction is generally considered rare, a number of recently emerging infectious diseases with load-dependent pathology have led to extinction in wildlife populations. Transmission is a critical factor affecting disease-induced extinction, but the relative importance of transmission compared to load-dependent host resistance and tolerance is currently unknown. Using a combination of models and experiments on an amphibian species suffering extirpations from the fungal pathogen Batrachochytrium dendrobatidis (Bd), we show that while transmission from an environmental Bd reservoir increased the ability of Bd to invade an amphibian population and the extinction risk of that population, Bd-induced extinction dynamics were far more sensitive to host resistance and tolerance than to Bd transmission. We demonstrate that this is a general result for load-dependent pathogens, where non-linear resistance and tolerance functions can interact such that small changes in these functions lead to drastic changes in extinction dynamics.

Batrachochytrium dendrobatidis and the Decline and Survival of the Relict Leopard Frog

Epizootic disease caused by the fungal pathogen Batrachochytrium dendrobatidis (Bd) is a major driver of amphibian declines, yet many amphibians declined before the pathogen was described. The Relict Leopard Frog, Rana onca (=Lithobates onca), was nearly extinct, with the exception of populations within a few geothermal springs. Growth of Bd, however, is limited by high water temperature, and geothermal springs may have provided refuge during outbreaks of chytridiomycosis. We conducted field surveys and laboratory experiments to assess the susceptibility of R. onca to Bd. In the field, we found Bd at one of the two areas where remnant populations of R. onca still occur, but not in the other. In the laboratory, we infected juvenile frogs from these two areas with two hypervirulent Bd isolates associated with declines in other ranid species. In our experiments, these Bd isolates did not affect survivorship of R. onca and most infections (64%) were cleared by the end of the experiments. We propose that R. onca either has inherent resistance to Bd or has recently evolved such resistance. These results may be important for conservation efforts aimed at establishing new populations of R. onca across a landscape where Bd exists. Resistance, however, varies among life stages, and we also did not assess Bd from the local environment. We caution that the resistance we observed for young frogs under laboratory conditions may not translate to the situation for R. onca in the wild.

Diversity in growth patterns among strains of the lethal fungal pathogen Batrachochytrium dendrobatidis across extended thermal optima

The thermal sensitivities of organisms regulate a wide range of ecological interactions, including host–parasite dynamics. The effect of temperature on disease ecology can be remarkably complex in disease systems where the hosts are ectothermic and where thermal conditions constrain pathogen reproductive rates. Amphibian chytridiomycosis, caused by the pathogen Batrachochytrium dendrobatidis (Bd), is a lethal fungal disease that is influenced by temperature. However, recent temperature studies have produced contradictory findings, suggesting that our current understanding of thermal effects on Bd may be incomplete. We investigated how temperature affects three different Bd strains to evaluate diversity in thermal responses. We quantified growth across the entire thermal range of Bd, and beyond the known thermal limits (T_max and T_min). Our results show that all Bd strains remained viable and grew following 24 h freeze (−12 °C) and heat shock (28 °C) treatments. Additionally, we found that two Bd strains had higher logistic growth rates (r) and carrying capacities (K) at the upper and lower extremities of the temperature range, and especially in low temperature conditions (2–3 °C). In contrast, a third strain exhibited relatively lower growth rates and carrying capacities at these same thermal extremes. Overall, our results suggest that there is considerable variation among Bd strains in thermal tolerance, and they establish a new thermal sensitivity profile for Bd. More generally, our findings point toward important questions concerning the mechanisms that dictate fungal thermal tolerances and temperature-dependent pathogenesis in other fungal disease systems.

First in Vivo Batrachochytrium dendrobatidis Transcriptomes Reveal Mechanisms of Host Exploitation, Host-Specific Gene Expression, and Expressed Genotype Shifts

For generalist pathogens, host species represent distinct selective environments, providing unique challenges for resource acquisition and defense from host immunity, potentially resulting in host-dependent differences in pathogen fitness. Gene expression modulation should be advantageous, responding optimally to a given host and mitigating the costs of generalism. Batrachochytrium dendrobatidis (Bd), a fungal pathogen of amphibians, shows variability in pathogenicity among isolates, and within-strain virulence changes rapidly during serial passages through artificial culture. For the first time, we characterize the transcriptomic profile of Bd in vivo, using laser-capture microdissection. Comparison of Bd transcriptomes (strain JEL423) in culture and in two hosts (Atelopus zeteki and Hylomantis lemur), reveals >2000 differentially expressed genes that likely include key Bd defense and host exploitation mechanisms. Variation in Bd transcriptomes from different amphibian hosts demonstrates shifts in pathogen resource allocation. Furthermore, expressed genotype variant frequencies of Bd populations differ between culture and amphibian skin, and among host species, revealing potential mechanisms underlying rapid changes in virulence and the possibility that amphibian community composition shapes Bd evolutionary trajectories. Our results provide new insights into how changes in gene expression and infecting population genotypes can be key to the success of a generalist fungal pathogen.

Insights From Genomics Into Spatial and Temporal Variation in Batrachochytrium dendrobatidis

Advances in genetics and genomics have provided new tools for the study of emerging infectious diseases. Researchers can now move quickly from simple hypotheses to complex explanations for pathogen origin, spread, and mechanisms of virulence. Here we focus on the application of genomics to understanding the biology of the fungal pathogen Batrachochytrium dendrobatidis (Bd), a novel and deadly pathogen of amphibians. We provide a brief history of the system, then focus on key insights into Bd variation garnered from genomics approaches, and finally, highlight new frontiers for future discoveries. Genomic tools have revealed unexpected complexity and variation in the Bd system suggesting that the history and biology of emerging pathogens may not be as simple as they initially seem.

Disentangling host, pathogen, and environmental determinants of a recently emerged wildlife disease: lessons from the first 15 years of amphibian chytridiomycosis research

The amphibian fungal disease chytridiomycosis, which affects species across all continents, recently emerged as one of the greatest threats to biodiversity. Yet, many aspects of the basic biology and epidemiology of the pathogen, Batrachochytrium dendrobatidis (Bd), are still unknown, such as when and from where did Bd emerge and what is its true ecological niche? Here, we review the ecology and evolution of Bd in the Americas and highlight controversies that make this disease so enigmatic. We explore factors associated with variance in severity of epizootics focusing on the disease triangle of host susceptibility, pathogen virulence, and environment. Reevaluating the causes of the panzootic is timely given the wealth of data on Bd prevalence across hosts and communities and the recent discoveries suggesting co-evolutionary potential of hosts and Bd. We generate a new species distribution model for Bd in the Americas based on over 30,000 records and suggest a novel future research agenda. Instead of focusing on pathogen "hot spots," we need to identify pathogen "cold spots" so that we can better understand what limits the pathogen's distribution. Finally, we introduce the concept of "the Ghost of Epizootics Past" to discuss expected patterns in postepizootic host communities.