Pathogen invasion history elucidates contemporary host pathogen dynamics

Chytrid infections exhibit historical spread and contemporary seasonality in a declining stream-breeding frog

Species with extensive geographical ranges pose special challenges to assessing drivers of wildlife disease, necessitating collaborative and large-scale analyses. The imperilled foothill yellow-legged frog (Rana boylii) inhabits a wide geographical range and variable conditions in rivers of California and Oregon (USA), and is considered threatened by the pathogen Batrachochytrium dendrobatidis (Bd). To assess drivers of Bd infections over time and space, we compiled over 2000 datapoints from R. boylii museum specimens (collected 1897-2005) and field samples (2005-2021) spanning 9° of latitude. We observed a south-to-north spread of Bd detections beginning in the 1940s and increase in prevalence from the 1940s to 1970s, coinciding with extirpation from southern latitudes. We detected eight high-prevalence geographical clusters through time that span the species' geographical range. Field-sampled male R. boylii exhibited the highest prevalence, and juveniles sampled in autumn exhibited the highest loads. Bd infection risk was highest in lower elevation rain-dominated watersheds, and with cool temperatures and low stream-flow conditions at the end of the dry season. Through a holistic assessment of relationships between infection risk, geographical context and time, we identify the locations and time periods where Bd mitigation and monitoring will be critical for conservation of this imperilled species.

Continent-wide recent emergence of a global pathogen in African amphibians

Introduction

Emerging infectious diseases are increasingly recognized as a global threat to wildlife. Pandemics in amphibians, caused by the fungal pathogen Batrachochytrium dendrobatidis (Bd), have resulted in biodiversity loss at a global scale. Genomic data suggest a complex evolutionary history of Bd lineages that vary in pathogenicity. Africa harbors a significant proportion of global amphibian biodiversity, and multiple Bd lineages are known to occur there; yet, despite the decline of many host species, there are currently no described Bd-epizootics. Here, we describe the historical and recent biogeographical spread of Bd and assess its risk to amphibians across the continent of Africa.

Methods

We provide a 165-year view of host-pathogen interactions by (i) employing a Bd assay to test 4,623 specimens (collected 1908–2013); (ii) compiling 12,297 published Bd records (collected 1852–2017); (iii) comparing the frequency of Bd-infected amphibians through time by both country and region; (iv) genotyping Bd lineages; (v) histologically identifying evidence of chytridiomycosis, and (vi) using a habitat suitability model to assess future Bd risk.

Results

We found a pattern of Bd emergence beginning largely at the turn of the century. From 1852–1999, we found low Bd prevalence (3.2% overall) and limited geographic spread, but after 2000 we documented a sharp increase in prevalence (18.7% overall), wider geographic spread, and multiple Bd lineages that may be responsible for emergence in different regions. We found that Bd risk to amphibians was highest in much of eastern, central, and western Africa.

Discussion

Our study documents a largely overlooked yet significant increase in a fungal pathogen that could pose a threat to amphibians across an entire continent. We emphasize the need to bridge historical and contemporary datasets to better describe and predict host-pathogen dynamics over larger temporal scales.

Host density has limited effects on pathogen invasion, disease-induced declines and within-host infection dynamics across a landscape of disease

1. Host density is hypothesized to be a major driver of variability in the responses and outcomes of wildlife populations following pathogen invasion. While the effects of host density on pathogen transmission have been extensively studied, these studies are dominated by theoretical analyses and small-scale experiments. This focus leads to an incomplete picture regarding how host density drives observed variability in disease outcomes in the field. 2. Here, we leveraged a dataset of hundreds of replicate amphibian populations that varied by orders of magnitude in host density. We used these data to test the effects of host density on three outcomes following the arrival of the amphibian-killing fungal pathogen Batrachochytrium dendrobatidis (Bd): the probability that Bd successfully invaded a host population and led to a pathogen outbreak, the magnitude of the host population-level decline following an outbreak and within-host infection dynamics that drive population-level outcomes in amphibian-pathogen systems. 3. Based on previous small-scale transmission experiments, we expected that populations with higher densities would be more likely to experience Bd outbreaks and would suffer larger proportional declines following outbreaks. To test these predictions, we developed and fitted a Hidden Markov Model that accounted for imperfectly observed disease outbreak states in the amphibian populations we surveyed. 4. Contrary to our predictions, we found minimal effects of host density on the probability of successful Bd invasion, the magnitude of population decline following Bd invasion and the dynamics of within-host infection intensity. Environmental conditions, such as summer temperature, winter severity and the presence of pathogen reservoirs, were more predictive of variability in disease outcomes. 5. Our results highlight the limitations of extrapolating findings from small-scale transmission experiments to observed disease trajectories in the field and provide strong evidence that variability in host density does not necessarily drive variability in host population responses following pathogen arrival. In an applied context, we show that feedbacks between host density and disease will not necessarily affect the success of reintroduction efforts in amphibian-Bd systems of conservation concern.

Environment is associated with chytrid infection and skin microbiome richness on an amphibian rich island (Taiwan)

Growing evidence suggests that the origins of the panzootic amphibian pathogens Batrachochytrium dendrobatidis (Bd) and Batrachochytrium salamandrivorans (Bsal) are in Asia. In Taiwan, an island hotspot of high amphibian diversity, no amphibian mass mortality events linked to Bd or Bsal have been reported. We conducted a multi-year study across this subtropical island, sampling 2517 individuals from 30 species at 34 field sites, between 2010 and 2017, and including 171 museum samples collected between 1981 and 2009. We analyzed the skin microbiome of 153 samples (6 species) from 2017 in order to assess any association between the amphibian skin microbiome and the probability of infection amongst different host species. We did not detect Bsal in our samples, but found widespread infection by Bd across central and northern Taiwan, both taxonomically and spatially. Museum samples show that Bd has been present in Taiwan since at least 1990. Host species, geography (elevation), climatic conditions and microbial richness were all associated with the prevalence of infection. Host life-history traits, skin microbiome composition and phylogeny were associated with lower prevalence of infection for high altitude species. Overall, we observed low prevalence and burden of infection in host populations, suggesting that Bd is enzootic in Taiwan where it causes subclinical infections. While amphibian species in Taiwan are currently threatened by habitat loss, our study indicates that Bd is in an endemic equilibrium with the populations and species we investigated. However, ongoing surveillance of the infection is warranted, as changing environmental conditions may disturb the currently stable equilibrium.

Social Behavior, Community Composition, Pathogen Strain, and Host Symbionts Influence Fungal Disease Dynamics in Salamanders

The emerging fungal pathogen, Batrachochytrium dendrobatidis (Bd), which can cause a fatal disease called chytridiomycosis, is implicated in the collapse of hundreds of host amphibian species. We describe chytridiomycosis dynamics in two co-occurring terrestrial salamander species, the Santa Lucia Mountains slender salamander, Batrachoseps luciae, and the arboreal salamander, Aneides lugubris. We (1) conduct a retrospective Bd-infection survey of specimens collected over the last century, (2) estimate present-day Bd infections in wild populations, (3) use generalized linear models (GLM) to identify biotic and abiotic correlates of infection risk, (4) investigate susceptibility of hosts exposed to Bd in laboratory trials, and (5) examine the ability of host skin bacteria to inhibit Bd in culture. Our historical survey of 2,866 specimens revealed that for most of the early 20th century (~1920–1969), Bd was not detected in either species. By the 1990s the proportion of infected specimens was 29 and 17% (B. luciae and A. lugubris, respectively), and in the 2010s it was 10 and 17%. This was similar to the number of infected samples from contemporary populations (2014–2015) at 10 and 18%. We found that both hosts experience signs of chytridiomycosis and suffered high Bd-caused mortality (88 and 71% for B. luciae and A. lugubris, respectively). Our GLM revealed that Bd-infection probability was positively correlated with intraspecific group size and proximity to heterospecifics but not to abiotic factors such as precipitation, minimum temperature, maximum temperature, mean temperature, and elevation, or to the size of the hosts. Finally, we found that both host species contain symbiotic skin-bacteria that inhibit growth of Bd in laboratory trials. Our results provide new evidence consistent with other studies showing a relatively recent Bd invasion of amphibian host populations in western North America and suggest that the spread of the pathogen may be enabled both through conspecific and heterospecific host interactions. Our results suggest that wildlife disease studies should assess host-pathogen dynamics that consider the interactions and effects of multiple hosts, as well as the historical context of pathogen invasion, establishment, and epizootic to enzootic transitions to better understand and predict disease dynamics.

Tracking, Synthesizing, and Sharing Global Batrachochytrium Data at AmphibianDisease.org

Emerging infectious diseases have been especially devastating to amphibians, the most endangered class of vertebrates. For amphibians, the greatest disease threat is chytridiomycosis, caused by one of two chytridiomycete fungal pathogens Batrachochytrium dendrobatidis (Bd) and Batrachochytrium salamandrivorans (Bsal). Research over the last two decades has shown that susceptibility to this disease varies greatly with respect to a suite of host and pathogen factors such as phylogeny, geography (including abiotic factors), host community composition, and historical exposure to pathogens; yet, despite a growing body of research, a comprehensive understanding of global chytridiomycosis incidence remains elusive. In a large collaborative effort, Bd-Maps was launched in 2007 to increase multidisciplinary investigations and understanding using compiled global Bd occurrence data (Bsal was not discovered until 2013). As its database functions aged and became unsustainable, we sought to address critical needs utilizing new technologies to meet the challenges of aggregating data to facilitate research on both Bd and Bsal. Here, we introduce an advanced central online repository to archive, aggregate, and share Bd and Bsal data collected from around the world. The Amphibian Disease Portal (https://amphibiandisease.org) addresses several critical community needs while also helping to build basic biological knowledge of chytridiomycosis. This portal could be useful for other amphibian diseases and could also be replicated for uses with other wildlife diseases. We show how the Amphibian Disease Portal provides: (1) a new repository for the legacy Bd-Maps data; (2) a repository for sample-level data to archive datasets and host published data with permanent DOIs; (3) a flexible framework to adapt to advances in field, laboratory, and informatics technologies; and (4) a global aggregation of Bd and Bsal infection data to enable and accelerate research and conservation. The new framework for this project is built using biodiversity informatics best practices and metadata standards to ensure scientific reproducibility and linkages across other biological and biodiversity repositories.

Longitudinal patterns in the skin microbiome of wild, individually marked frogs from the Sierra Nevada, California

The amphibian skin microbiome has been the focus of numerous studies because of the protective effects that some bacteria provide against the pathogenic fungus Batrachochytrium dendrobatidis, which has caused a global panzootic among amphibians. However, the mechanisms driving community structure and function in the amphibian skin microbiome are still poorly understood, and longitudinal analyses of the skin microbiome have not yet been conducted in wild populations. In this study, we investigate longitudinal patterns in the skin microbiome of 19 individually marked adult frogs from two wild populations of the endangered Sierra Nevada yellow-legged frog (Rana sierrae), sampled over the course of 2 years. We found that individuals with low bacterial diversity (dominated by order Burkhorderiales) had significantly more stable bacterial communities than those with higher diversity. Amplicon sequence variants (ASVs) with high relative abundance were significantly less transient than those with low relative abundance, and ASVs with intermediate-level relative abundances experienced the greatest volatility over time. Based on these results, we suggest that efforts to develop probiotic treatments to combat B. dendrobatidis should focus on bacteria that are found at high relative abundances in some members of a population, as these strains are more likely to persist and remain stable in the long term.

Divergent regional evolutionary histories of a devastating global amphibian pathogen

Emerging infectious diseases are a pressing threat to global biological diversity. Increased incidence and severity of novel pathogens underscores the need for methodological advances to understand pathogen emergence and spread. Here, we use genetic epidemiology to test, and challenge, key hypotheses about a devastating zoonotic disease impacting amphibians globally. Using an amplicon-based sequencing method and non-invasive samples we retrospectively explore the history of the fungal pathogen Batrachochytrium dendrobatidis (Bd) in two emblematic amphibian systems: the Sierra Nevada of California and Central Panama. The hypothesis in both regions is the hypervirulent Global Panzootic Lineage of Bd (BdGPL) was recently introduced and spread rapidly in a wave-like pattern. Our data challenge this hypothesis by demonstrating similar epizootic signatures can have radically different underlying evolutionary histories. In Central Panama, our genetic data confirm a recent and rapid pathogen spread. However, BdGPL in the Sierra Nevada has remarkable spatial structuring, high genetic diversity and a relatively older history inferred from time-dated phylogenies. Thus, this deadly pathogen lineage may have a longer history in some regions than assumed, providing insights into its origin and spread. Overall, our results highlight the importance of integrating observed wildlife die-offs with genetic data to more accurately reconstruct pathogen outbreaks.

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.

Differences in Fungal Disease Dynamics in Co-occurring Terrestrial and Aquatic Amphibians

The fungal pathogen, Batrachochytrium dendrobatidis (Bd), has devastated biodiversity and ecosystem health and is implicated as a driver of mass amphibian extinctions. This 100-year study investigates which environmental factors contribute to Bd prevalence in a fully terrestrial species, and determines whether infection patterns differ between a fully terrestrial amphibian and more aquatic host species. We performed a historical survey to quantify Bd prevalence in 1127 Batrachoseps gregarius museum specimens collected from 1920 to 2000, and recent data from 16 contemporary (live-caught) B. gregarius populations from the southwestern slopes of the Sierra Nevada mountains in California, USA. We compared these results to Bd detection rates in 1395 historical and 1033 contemporary specimens from 10 species of anurans and 427 historical Taricha salamander specimens collected throughout the Sierra Nevada mountains. Our results indicate that Bd dynamics in the entirely terrestrial species, B. gregarius, differ from aquatic species in the same region in terms of both seasonal patterns of Bd abundance and in the possible timing of Bd epizootics.