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

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.

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.

INFECTION DYNAMICS OF BATRACHOCHYTRIUM DENDROBATIDIS IN TWO FROG SPECIES INHABITING QUITO'S METROPOLITAN GUANGÜILTAGUA PARK, ECUADOR

Batrachochytrium dendrobatidis (Bd) infection is one of the principal causes of amphibian declines worldwide. The presence of Bd has been determined in Gastrotheca riobambae tadpoles that inhabit ponds in Quito's Metropolitan Guangüiltagua Park, Ecuador. This study sought to determine whether these tadpoles are infected and to determine the presence of chytridiomycosis in another frog species, Pristimantis unistrigatus, which also inhabits the park and has different reproductive biology and distinct behavioral habits. We used end-point and real-time PCR techniques to detect and quantify Bd infection. At 1 yr, samples were taken from the skin of P. unistrigatus using swabs and were also taken from the mouthparts of G. riobambae tadpoles. It was found that the two species were infected with a Bd prevalence of 39% (53/135) in G. riobambae tadpoles and 15% (57/382) in P. unistrigatus frogs. The two types of samples (tissue and swabs) from mouthparts showed differences in the zoospores per microliter loads (x̄=1,376.7±3,450.2 vs. x̄=285.0±652.3). Moreover, a correlation (r2=0.621) was discovered between the monthly mean maximum temperature of the pond with disease prevalence in G. riobambae tadpoles. Infection levels in the P. unistrigatus population varied significantly over time, and distance to the pond was a determinant factor for infection intensity.

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.

Batrachochytrium dendrobatidis infection in amphibians predates first known epizootic in Costa Rica

Emerging infectious diseases are a growing threat to biodiversity worldwide. Outbreaks of the infectious disease chytridiomycosis, caused by the fungal pathogen Batrachochytrium dendrobatidis (Bd), are implicated in the decline and extinction of numerous amphibian species. In Costa Rica, a major decline event occurred in 1987, more than two decades before this pathogen was discovered. The loss of many species in Costa Rica is assumed to be due to Bd-epizootics, but there are few studies that provide data from amphibians in the time leading up to the proposed epizootics. In this study, we provide new data on Bd infection rates of amphibians collected throughout Costa Rica, in the decades prior to the epizootics. We used a quantitative PCR assay to test for Bd presence in 1016 anuran museum specimens collected throughout Costa Rica. The earliest specimen that tested positive for Bd was collected in 1964. Across all time periods, we found an overall infection rate (defined as the proportion of Bd-positive individuals) of 4%. The number of infected individuals remained relatively low across all species tested and the range of Bd-positive specimens was shown to be geographically constrained up until the 1980s; when epizootics are hypothesized to have occurred. After that time, infection rate increased three-fold, and the range of specimens tested positive for Bd increased, with Bd-positive specimens collected across the entire country. Our results suggest that Bd dynamics in Costa Rica are more complicated than previously thought. The discovery of Bd's presence in the country preceding massive declines leads to a number of different hypotheses: 1) Bd invaded Costa Rica earlier than previously known, and spread more slowly than previously reported; 2) Bd invaded multiple times and faded out; 3) an endemic Bd lineage existed; 4) an earlier Bd lineage evolved into the current Bd lineage or hybridized with an invasive lineage; or 5) an earlier Bd lineage went extinct and a new invasion event occurred causing epizootics. To help visualize areas where future studies should take place, we provide a Bd habitat suitability model trained with local data. Studies that provide information on genetic lineages of Bd are needed to determine the most plausible spatial-temporal, host-pathogen dynamics that could best explain the epizootics resulting in amphibian declines in Costa Rica and throughout Central America.

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.

Non-lethal isolation of the fungal pathogen Batrachochytrium dendrobatidis (Bd) from amphibians

The ability to isolate and purify pathogens is important for the study of infectious disease. A protocol for isolating Batrachochytrium dendrobatidis (Bd), a lethal pathogen of amphibians, has been available for over a decade, but the method relies on sacrificing infected animals. We validated a non-lethal protocol for Bd isolation that uses biopsy punches from toe webbing to collect skin samples from live amphibians in remote field locations. We successfully isolated Bd from the Cascades frog Rana cascadae and found a positive association between Bd infection and probability of Bd growth in culture. Recapture rates of sampled animals suggest that our isolation protocol did not affect frog survival. The ability to collect isolates from live animals will facilitate investigations of the biology of Bd and enhance amphibian conservation efforts.

Pathogen invasion and non-epizootic dynamics in Pacific newts in California over the last century

Emerging infectious disease is a growing threat to global biodiversity. The infectious disease chytridiomycosis, caused by the fungal pathogen Batrachochytrium dendrobatidis (Bd) has led to the decline and extinction of hundreds of amphibian species. Severe Bd-caused epizootics have been documented in North, Central and South America—with many of the research focused on anurans. California, where Bd-related epizootics and amphibian declines have been reported, has some of the highest diversity of salamanders. After more than a decade since the first known epizootic in California, little is known about Bd disease dynamics in salamanders. Pacific newts (Genus: Taricha) are ideal study species because of their abundance, wide geographic range, occurrence in both aquatic and terrestrial habitats, and how little is known about Bd infection dynamics for this group. We conducted a retrospective study to determine the relationship between Pacific newts and the fungal pathogen. We tested 1895 specimens collected between 1889–2009 and found no evidence of Bd-infected Pacific newts until the late 1940’s. Although we estimate that Bd emerged in this genus and rapidly spread geographically throughout California, we did not find evidence for epizootic dynamics. Bd infection prevalence and intensity, two measures commonly used to estimate dynamics, remained consistently low over time; suggesting Pacific newts may not be highly susceptible. Also, we found the timing of first Bd emergence in Pacific newts predate Bd emergence in other California salamander species. In addition, we found several environmental and anthropogenic factors correlated with Bd prevalence which may help explain Bd disease dynamics in the genus Taricha. Pacific newts may be a reservoir species that signal pathogen invasion into California salamanders, though further studies are needed.

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.

Rapid extirpation of a North American frog coincides with an increase in fungal pathogen prevalence: Historical analysis and implications for reintroduction

As extinctions continue across the globe, conservation biologists are turning to species reintroduction programs as one optimistic tool for addressing the biodiversity crisis. For repatriation to become a viable strategy, fundamental prerequisites include determining the causes of declines and assessing whether the causes persist in the environment. Invasive species-especially pathogens-are an increasingly significant factor contributing to biodiversity loss. We hypothesized that Batrachochytrium dendrobatidis (Bd), the causative agent of the deadly amphibian disease chytridiomycosis, was important in the rapid (<10 years) localized extirpation of a North American frog (Rana boylii) and that Bd remains widespread among extant amphibians in the region of extirpation. We used an interdisciplinary approach, combining interviews with herpetological experts, analysis of archived field notes and museum specimen collections, and field sampling of the extant amphibian assemblage to examine (1) historical relative abundance of R. boylii; (2) potential causes of R. boylii declines; and (3) historical and contemporary prevalence of Bd. We found that R. boylii were relatively abundant prior to their rapid extirpation, and an increase in Bd prevalence coincided with R. boylii declines during a time of rapid change in the region, wherein backcountry recreation, urban development, and the amphibian pet trade were all on the rise. In addition, extreme flooding during the winter of 1969 coincided with localized extirpations in R. boylii populations observed by interview respondents. We conclude that Bd likely played an important role in the rapid extirpation of R. boylii from southern California and that multiple natural and anthropogenic factors may have worked in concert to make this possible in a relatively short period of time. This study emphasizes the importance of recognizing historical ecological contexts in making future management and reintroduction decisions.

Greater Species Richness of Bacterial Skin Symbionts Better Suppresses the Amphibian Fungal Pathogen Batrachochytrium Dendrobatidis

The symbiotic microbes that grow in and on many organisms can play important roles in protecting their hosts from pathogen infection. While species diversity has been shown to influence community function in many other natural systems, the question of how species diversity of host-associated symbiotic microbes contributes to pathogen resistance is just beginning to be explored. Understanding diversity effects on pathogen resistance could be particularly helpful in combating the fungal pathogen Batrachochytrium dendrobatidis (Bd) which has caused dramatic population declines in many amphibian species and is a major concern for amphibian conservation. Our study investigates the ability of host-associated bacteria to inhibit the proliferation of Bd when grown in experimentally assembled biofilm communities that differ in species number and composition. Six bacterial species isolated from the skin of Cascades frogs (Rana cascadae) were used to assemble bacterial biofilm communities containing 1, 2, 3, or all 6 bacterial species. Biofilm communities were grown with Bd for 7 days following inoculation. More speciose bacterial communities reduced Bd abundance more effectively. This relationship between bacterial species richness and Bd suppression appeared to be driven by dominance effects-the bacterial species that were most effective at inhibiting Bd dominated multi-species communities-and complementarity: multi-species communities inhibited Bd growth more than monocultures of constituent species. These results underscore the notion that pathogen resistance is an emergent property of microbial communities, a consideration that should be taken into account when designing probiotic treatments to reduce the impacts of infectious disease.