Invasion of the Fungal Pathogen Batrachochytrium dendrobatidis on California Islands

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.

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.

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.

Pathogen invasion history elucidates contemporary host pathogen dynamics

Amphibians, the most threatened group of vertebrates, are seen as indicators of the sixth mass extinction on earth. Thousands of species are threatened with extinction and many have been affected by an emerging infectious disease, chytridiomycosis, caused by the fungal pathogen, Batrachochytrium dendrobatidis (Bd). However, amphibians exhibit different responses to the pathogen, such as survival and population persistence with infection, or mortality of individuals and complete population collapse after pathogen invasion. Multiple factors can affect host pathogen dynamics, yet few studies have provided a temporal view that encompasses both the epizootic phase (i.e. pathogen invasion and host collapse), and the transition to a more stable co-existence (i.e. recovery of infected host populations). In the Sierra Nevada mountains of California, USA, conspecific populations of frogs currently exhibit dramatically different host/ Bd-pathogen dynamics. To provide a temporal context by which present day dynamics may be better understood, we use a Bd qPCR assay to test 1165 amphibian specimens collected between 1900 and 2005. Our historical analyses reveal a pattern of pathogen invasion and eventual spread across the Sierra Nevada over the last century. Although we found a small number of Bd-infections prior to 1970, these showed no sign of spread or increase in infection prevalence over multiple decades. After the late 1970s, when mass die offs were first noted, our data show Bd as much more prevalent and more spatially spread out, suggesting epizootic spread. However, across the ~400km2 area, we found no evidence of a wave-like pattern, but instead discovered multiple, nearly-simultaneous invasions within regions. We found that Bd invaded and spread in the central Sierra Nevada (Yosemite National Park area) about four decades before it invaded and spread in the southern Sierra Nevada (Sequoia and Kings Canyon National Parks area), and suggest that the temporal pattern of pathogen invasion may help explain divergent contemporary host pathogen dynamics.

Introduced bullfrog facilitates pathogen invasion in the western United States

Batrachochytrium dendrobatidis (Bd), a causal agent of the amphibian fungal skin disease chytridiomycosis, has been implicated in the decline and extinction of over 200 species worldwide since the 1970s. Despite almost two decades of research, the history of Bd and its global spread is not well understood. However, the spread of the Global Panzootic Lineage of Bd (Bd-GPL), the lineage associated with amphibian die-offs, has been linked with the American bullfrog (Rana [Aqurana] catesbeiana) and global trade. Interestingly, R. catesbeiana is native to the eastern U.S., where no Bd-related declines have been observed despite Bd’s presence since the late 1800s. In contrast Bd has been found to have emerged in California and Mexico in the 1960s and 1970s, after which epizootics (i.e., epidemics in wildlife) ensued. We hypothesize that Bd-GPL spread from the eastern U.S. with the introduction of R. catesbeiana into the western US, resulting in epizootics and declines of native host species. Using museum records, we investigated the historical relationship between R. catesbeiana and Bd invasion in the western US and found that R. catesbeiana arrived in the same year or prior to Bd in most western watersheds that had data for both species, suggesting that Bd-GPL may have originated in the eastern US and R. catesbeiana may have facilitated Bd invasion in the western US. To predict areas with greatest suitability for Bd, we created a suitability model by integrating habitat suitability and host availability. When we incorporated invasion history with high Bd suitability, we found that watersheds with non-native R. catesbeiana in the mountain ranges of the West Coast have the highest disease risk. These findings shed light on the invasion history and disease dynamics of Bd in North America. Targeted historical surveys using archived specimens in natural history collections and present-day field surveys along with more localized, community-level studies, monitoring, and surveillance are needed to further test this hypothesis and grow our understanding of the disease ecology and host-pathogen dynamics of Bd.

Recent Emergence of a Chytrid Fungal Pathogen in California Cascades Frogs (Rana cascadae)

The pathogenic fungus Batrachochytrium dendrobatidis (Bd) has been associated with global amphibian declines, but it is often difficult to discern the relative importance of Bd as a causal agent in declines that have already occurred. Retrospective analyses of museum specimens have allowed researchers to associate the timing of Bd arrival with the timing of past amphibian declines. Cascades frogs (Rana cascadae) have experienced dramatic declines in northern California, but it is not clear whether the onset of these declines corresponds to the arrival of Bd. We used quantitative real-time PCR assays of samples collected from museum specimens to determine historical Bd prevalence in the northern California range of Cascades frogs. We detected Bd in 13 of 364 (3.5%) Cascades frog specimens collected between 1907 and 2003, with the first positive result from 1978. A Bayesian analysis suggested that Bd arrived in the region between 1973 and 1978, which corresponds well with the first observations of declines in the 1980s.