Investigating Differences across Host Species and Scales to Explain the Distribution of the Amphibian Pathogen Batrachochytrium dendrobatidis

Many pathogens infect more than one host species, and clarifying how these different hosts contribute to pathogen dynamics can facilitate the management of pathogens and can lend insight into the functioning of pathogens in ecosystems. In this study, we investigated a suite of native and non-native amphibian hosts of the pathogen Batrachochytrium dendrobatidis (Bd) across multiple scales to identify potential mechanisms that may drive infection patterns in the Colorado study system. Specifically, we aimed to determine if: 1) amphibian populations vary in Bd infection across the landscape, 2) amphibian community composition predicts infection (e.g., does the presence or abundance of any particular species influence infection in others?), 3) amphibian species vary in their ability to produce infectious zoospores in a laboratory infection, 4) heterogeneity in host ability observed in the laboratory scales to predict patterns of Bd prevalence in the landscape. We found that non-native North American bullfrogs (Lithobates catesbeianus) are widespread and have the highest prevalence of Bd infection relative to the other native species in the landscape. Additionally, infection in some native species appears to be related to the density of sympatric L. catesbeianus populations. At the smaller host scale, we found that L. catesbeianus produces more of the infective zoospore stage relative to some native species, but that this zoospore output does not scale to predict infection in sympatric wild populations of native species. Rather, landscape level infection relates most strongly to density of hosts at a wetland as well as abiotic factors. While non-native L. catesbeianus have high levels of Bd infection in the Colorado Front Range system, we also identified Bd infection in a number of native amphibian populations allopatric with L. catesbeianus, suggesting that multiple host species are important contributors to the dynamics of the Bd pathogen in this landscape.     

Host species composition influences infection severity among amphibians in the absence of spillover transmission

Wildlife epidemiological outcomes can depend strongly on the composition of an ecological community, particularly when multiple host species are affected by the same pathogen. However, the relationship between host species richness and disease risk can vary with community context and with the degree of spillover transmission that occurs among co‐occurring host species. We examined the degree to which host species composition influences infection by Batrachochytrium dendrobatidis (Bd), a widespread fungal pathogen associated with amphibian population declines around the world, and whether transmission occurs from one highly susceptible host species to other co‐occurring host species. By manipulating larval assemblages of three sympatric amphibian species in the laboratory, we characterized the relationship between host species richness and infection severity, whether infection mediates growth and survivorship differently across various combinations of host species, and whether Bd is transmitted from experimentally inoculated tadpoles to uninfected tadpoles. We found evidence of a dilution effect where Bd infection severity was dramatically reduced in the most susceptible of the three host species (Anaxyrus boreas). Infection also mediated survival and growth of all three host species such that the presence of multiple host species had both positive (e.g., infection reduction) and negative (e.g., mortality) effects on focal species. However, we found no evidence that Bd infection is transmitted by this species. While these results demonstrate that host species richness as well as species identity underpin infection dynamics in this system, dilution is not the product of reduced transmission via fewer infectious individuals of a susceptible host species. We discuss various mechanisms, including encounter reduction and antagonistic interactions such as competition and opportunistic cannibalism that may act in concert to mediate patterns of infection severity, growth, and mortality observed in multihost communities.     

Emerging disease of amphibians cured by elevated body temperature

The emerging infectious disease chytridiomycosis is thought to have contributed to many of the recent alarming declines in amphibian populations. Mortalities associated with these declines have often occurred during cooler seasons and at high elevations, suggesting that environmental temperature may be an important factor in disease emergence. We found that thermal environment affects the progress of the disease, and that housing frogs Litoria chloris at an environmental temperature of 37 degrees C for less than 16 h can clear them of the chytrid pathogen Batrachochytrium dendrobatidis. Our experiment demonstrated that elevated body temperatures similar to those experienced in behavioral fever and during normal thermoregulation can clear frogs of chytrid infection; therefore, variation in thermoregulatory opportunities and behaviors are likely to contribute to the differences in disease incidence observed among host species, populations, and regions. Although further refinement of the technique is needed to encompass various host species, appropriately applied thermal manipulations of amphibians and their enclosures may prove to be a safe and effective way of eliminating the fungal pathogen from captive amphibian populations and: preventing accidental spread of the pathogen when animals are translocated or released from captivity.     

Interacting Symbionts and Immunity in the Amphibian Skin Mucosome Predict Disease Risk and Probiotic Effectiveness

Pathogenesis is strongly dependent on microbial context, but development of probiotic therapies has neglected the impact of ecological interactions. Dynamics among microbial communities, host immune responses, and environmental conditions may alter the effect of probiotics in human and veterinary medicine, agriculture and aquaculture, and the proposed treatment of emerging wildlife and zoonotic diseases such as those occurring on amphibians or vectored by mosquitoes. Here we use a holistic measure of amphibian mucosal defenses to test the effects of probiotic treatments and to assess disease risk under different ecological contexts. We developed a non-invasive assay for antifungal function of the skin mucosal ecosystem (mucosome function) integrating host immune factors and the microbial community as an alternative to pathogen exposure experiments. From approximately 8500 amphibians sampled across Europe, we compared field infection prevalence with mucosome function against the emerging fungal pathogen Batrachochytrium dendrobatidis. Four species were tested with laboratory exposure experiments, and a highly susceptible species, Alytes obstetricans, was treated with a variety of temperature and microbial conditions to test the effects of probiotic therapies and environmental conditions on mucosome function. We found that antifungal function of the amphibian skin mucosome predicts the prevalence of infection with the fungal pathogen in natural populations, and is linked to survival in laboratory exposure experiments. When altered by probiotic therapy, the mucosome increased antifungal capacity, while previous exposure to the pathogen was suppressive. In culture, antifungal properties of probiotics depended strongly on immunological and environmental context including temperature, competition, and pathogen presence. Functional changes in microbiota with shifts in temperature provide an alternative mechanistic explanation for patterns of disease susceptibility related to climate beyond direct impact on host or pathogen. This nonlethal management tool can be used to optimize and quickly assess the relative benefits of probiotic therapies under different climatic, microbial, or host conditions.     

Composition of symbiotic bacteria predicts survival in Panamanian golden frogs infected with a lethal fungus

Symbiotic microbes can dramatically impact host health and fitness, and recent research in a diversity of systems suggests that different symbiont community structures may result in distinct outcomes for the host. In amphibians, some symbiotic skin bacteria produce metabolites that inhibit the growth of Batrachochytrium dendrobatidis (Bd), a cutaneous fungal pathogen that has caused many amphibian population declines and extinctions. Treatment with beneficial bacteria (probiotics) prevents Bd infection in some amphibian species and creates optimism for conservation of species that are highly susceptible to chytridiomycosis, the disease caused by Bd. In a laboratory experiment, we used Bd-inhibitory bacteria from Bd-tolerant Panamanian amphibians in a probiotic development trial with Panamanian golden frogs, Atelopus zeteki, a species currently surviving only in captive assurance colonies. Approximately 30% of infected golden frogs survived Bd exposure by either clearing infection or maintaining low Bd loads, but this was not associated with probiotic treatment. Survival was instead related to initial composition of the skin bacterial community and metabolites present on the skin. These results suggest a strong link between the structure of these symbiotic microbial communities and amphibian host health in the face of Bd exposure and also suggest a new approach for developing amphibian probiotics.     

Expression Profiling the Temperature-Dependent Amphibian Response to Infection by Batrachochytrium dendrobatidis

Amphibians are experiencing a panzootic of unprecedented proportions caused by the emergence of Batrachochytrium dendrobatidis (Bd). However, all species are not equally at risk of infection, and risk is further modified by environmental variables, specifically temperature. In order to understand how, and when, hosts mount a response to Bd we analysed infection dynamics and patterns of gene expression in the model amphibian species Silurana (Xenopus) tropicalis. Mathematical modelling of infection dynamics demonstrate the existence of a temperature-dependent protective response that is largely independent of the intrinsic growth-rate of Bd. Using temporal expression-profiling by microarrays and qRT-PCR, we characterise this response in the main amphibian lymphoid tissue, the spleen. We demonstrate that clearance of Bd at the host-optimal temperature is not clearly associated with an adaptive immune response, but rather is correlated with the induction of components of host innate immunity including the expression of genes that are associated with the production of the antimicrobial skin peptide preprocareulein (PPCP) as well as inflammatory responses. We find that adaptive immunity appears to be lacking at host-optimal temperatures. This suggests that either Bd does not stimulate, or suppresses, adaptive immunity, or that trade-offs exist between innate and adaptive limbs of the amphibian immune system. At cold temperatures, S. tropicalis loses the ability to mount a PPCP-based innate response, and instead manifests a more pronounced inflammatory reaction that is characterised by the production of proteases and higher pathogen burdens. This study demonstrates the temperature-dependency of the amphibian response to infection by Bd and indicates the influence that changing climates may exert on the ectothermic host response to pathogens.     

Development and Infectious Disease in Hosts with Complex Life Cycles

Metamorphosis is often characterized by profound changes in morphology and physiology that can affect the dynamics of species interactions. For example, the interaction between a pathogen and its host may differ depending on the life stage of the host or pathogen. One pathogen that infects hosts with complex life cycles is the emerging fungal pathogen of amphibians, Batrachochytrium dendrobatidis (Bd). We sought to determine how conditions at the larval stage can affect variation in development and patterns of Bd infection across amphibian life stages. We used outdoor experimental mesocosms to simulate natural pond habitats and manipulated the presence of Bd, the larval density, and the number of host species in larvae of two co-occurring amphibian species (Rana cascadae and Pseudacris regilla). We found that infection differed between species throughout development; P. regilla consistently had higher infection severity compared to R. cascadae. Additionally, while up to 100% of larvae were infected, only 18.2% of R. cascadae and 81.5% of P. regilla were infected after metamorphosis. This indicates that amphibians have the ability to recover from Bd infection as they undergo metamorphosis. Higher larval densities in P. regilla led to a shorter larval period, and individuals with a shorter larval period had lower infection severity. This led to a trend where P. regilla larvae reared at high densities tended to have lower infection prevalence after metamorphosis. We also found that exposure to Bd increased larval mortality and prolonged the larval period in P. regilla, indicating that P. regilla are susceptible to the negative effects of Bd as larvae. This study demonstrates that host density, species composition, and pathogen exposure may all interact to influence development and infection in hosts with complex life cycles.     

Partitioning the net effect of host diversity on an emerging amphibian pathogen

The ‘dilution effect’ (DE) hypothesis predicts that diverse host communities will show reduced disease. The underlying causes of pathogen dilution are complex, because they involve non-additive (driven by host interactions and differential habitat use) and additive (controlled by host species composition) mechanisms. Here, we used measures of complementarity and selection traditionally employed in the field of biodiversity–ecosystem function (BEF) to quantify the net effect of host diversity on disease dynamics of the amphibian-killing fungus Batrachochytrium dendrobatidis (Bd). Complementarity occurs when average infection load in diverse host assemblages departs from that of each component species in uniform populations. Selection measures the disproportionate impact of a particular species in diverse assemblages compared with its performance in uniform populations, and therefore has strong additive and non-additive properties. We experimentally infected tropical amphibian species of varying life histories, in single- and multi-host treatments, and measured individual Bd infection loads. Host diversity reduced Bd infection in amphibians through a mechanism analogous to complementarity (sensu BEF), potentially by reducing shared habitat use and transmission among hosts. Additionally, the selection component indicated that one particular terrestrial species showed reduced infection loads in diverse assemblages at the expense of neighbouring aquatic hosts becoming heavily infected. By partitioning components of diversity, our findings underscore the importance of additive and non-additive mechanisms underlying the DE.     

After the epizootic: Host–pathogen dynamics in montane tropical amphibian communities with high prevalence of chytridiomycosis

The amphibian fungal disease chytridiomycosis, caused by Batrachochytrium dendrobatidis (Bd), poses a great threat to global amphibian biodiversity. In Peruvian cloud forests of the Kosñipata Valley of Manu National Park where chytrid infection is highly prevalent, we have monitored species-rich amphibian communities since 1996. An epizootic of chytridiomycosis is thought to have caused the disappearance of 35% of species richness in the early 2000s. We investigated the post-epizootic Bd prevalence and infection intensity within the remnant amphibian community from 2008 to 2015, and modeled Bd dynamics as a function of species, season, reproductive mode, life stage, and elevation. Prevalence was higher in 2012–2015 than in 2008–2009, but overall prevalence has remained fairly constant (~50%) post-epizootic. We also found that while prevalence decreased with elevation during the wet season, it generally increased with elevation during the dry season, potentially due to seasonal changes in temperature and precipitation. In aquatic habitats, Bd is likely maintained through a single, stream-breeding, putative reservoir species (which survived epizootics, in contrast to other aquatic-breeding species). The now-dominant terrestrial-breeding species allow Bd to persist and spread in terrestrial habitats, possibly through individual dispersal into naïve areas. We conclude that Bd prevalence in the Kosñipata Valley has stabilized over time, suggesting that Bd is now enzootic. Long-term monitoring of host infection is important because temporal changes in prevalence and infection intensity can cause changes in host species richness and abundance, which in turn may alter the trajectory of host–pathogen dynamics.     

Interaction between breeding habitat and elevation affects prevalence but not infection intensity of Batrachochytrium dendrobatidis in Brazilian anuran assemblages

Chytridiomycosis, an infectious disease of amphibians, is caused by the fungus Batrachochytrium dendrobatidis (Bd) and has been linked to declining amphibian populations worldwide. The susceptibility of amphibians to chytridiomycosis-induced population declines is potentially influenced by many factors, including environmental characteristics, differences among host species and the growth of the pathogen itself. We investigated the effects of elevation and breeding habitat on Bd prevalence and individual infection intensity (zoospore loads) in 3 anuran assemblages of the Atlantic Coastal Forest of Brazil. Bd infection intensity was strongly influenced by elevation and breeding habitat, but we found no evidence of an interaction between those 2 variables in explaining the number of zoospores sampled from individual frogs. In contrast, Bd infection odds were predicted by elevation and by an interaction between elevation and breeding habitat, such that frogs had a higher probability of Bd infection in lotic habitats at low elevations. Our results indicate that Bd persists across a wide variety of habitats and elevations in the Atlantic Coastal Forest. Prevalence and infection intensity of Bd are highest at high elevations where overall environmental conditions for Bd are most favorable. In addition, at low elevations amphibian host habitat choice is also an important determinant of infection. Our study highlights the need to investigate interacting variables of host ecology and the environment simultaneously.     

Infection risk decreases with increasing mismatch in host and pathogen environmental tolerances

The fungal pathogen Batrachochytrium dendrobatidis (Bd) has caused the greatest known wildlife pandemic, infecting over 500 amphibian species. It remains unclear why some host species decline from disease‐related mortality whereas others persist. We introduce a conceptual model that predicts that infection risk in ectotherms will decrease as the difference between host and pathogen environmental tolerances (i.e. tolerance mismatch) increases. We test this prediction using both local‐scale data from Costa Rica and global analyses of over 11 000 Bd infection assays. We find that infection prevalence decreases with increasing thermal tolerance mismatch and with increasing host tolerance of habitat modification. The relationship between environmental tolerance mismatches and Bd infection prevalence is generalisable across multiple amphibian families and spatial scales, and the magnitude of the tolerance mismatch effect depends on environmental context. These findings may help explain patterns of amphibian declines driven by a global wildlife pandemic.     

Disease hotspots or hot species? Infection dynamics in multi‐host metacommunities controlled by species identity,not source location

Pathogen persistence in host communities is influenced by processes operating at the individual host to landscape‐level scale, but isolating the relative contributions of these processes is challenging. We developed theory to partition the influence of host species, habitat patches and landscape connectivity on pathogen persistence within metacommunities of hosts and pathogens. We used this framework to quantify the contributions of host species composition and habitat patch identity on the persistence of an amphibian pathogen across the landscape. By sampling over 11 000 hosts of six amphibian species, we found that a single host species could maintain the pathogen in 91% of observed metacommunities. Moreover, this dominant maintenance species contributed, on average, twice as much to landscape‐level pathogen persistence compared to the most influential source patch in a metacommunity. Our analysis demonstrates substantial inequality in how species and patches contribute to pathogen persistence, with important implications for targeted disease management.     

Risk of Exposure to Eastern Equine Encephalomyelitis Virus Increases with the Density of Northern Cardinals

For a variety of infectious diseases, the richness of the community of potential host species has emerged as an important factor in pathogen transmission, whereby a higher richness of host species is associated with a lowered disease risk. The proposed mechanism driving this pattern is an increased likelihood in species-rich communities that infectious individuals will encounter dead-end hosts. Mosquito-borne pathogen systems potentially are exceptions to such “dilution effects” because mosquitoes vary their rates of use of vertebrate host species as bloodmeal sources relative to host availabilities. Such preferences may violate basic assumptions underlying the hypothesis of a dilution effect in pathogen systems. Here, we describe development of a model to predict exposure risk of sentinel chickens to eastern equine encephalitis virus (EEEV) in Walton County, Florida between 2009 and 2010 using avian species richness as well as densities of individual host species potentially important to EEEV transmission as candidate predictor variables. We found the highest support for the model that included the density of northern cardinals, a highly preferred host of mosquito vectors of EEEV, as a predictor variable. The highest-ranking model also included Culiseta melanura abundance as a predictor variable. These results suggest that mosquito preferences for vertebrate hosts influence pathogen transmission.     

The relationship between the emergence of Batrachochytrium dendrobatidis,the international trade in amphibians and introduced amphibian species

Chytridiomycosis is an emerging infectious disease of amphibians caused by the chytrid Batrachochytrium dendrobatidis. The disease has been associated with global amphibian declines and species extinctions, however the principle drivers that underly the emergence of chytridiomycosis remain unclear. Current evidence suggests that the world trade in amphibians is implicated in the emergence of chytridiomycosis. Here, we review the evidence that the amphibian trade is driving the emergence of chytridiomycosis by (1) spreading infected animals worldwide, (2) introducing non-native infected animals into naïve populations and (3) amplifying infection of amphibians by co-housing, followed by untreated discharge of infectious zoospores into water supplies. We conclude that the evidence that the amphibian trade is contributing to the spread of Batrachochytrium dendrobatidis is strong, and that specific actions are necessary to prevent the introduction of the pathogen into thus-far uninfected areas. Specifically, we recommend the development of national risk-abatement plans, focused on firstly preventing introduction of Bd into disease free areas, and secondly, decreasing the impact of the disease on populations that are currently infected.     

External Reinfection of a Fungal Pathogen Does not Contribute to Pathogen Growth

Chytridiomycosis is an emerging infectious disease of amphibians caused by the fungal pathogen Batrachochytrium dendrobatidis (Bd), which has led to devastating declines in amphibian populations worldwide. Current theory predicts that Bd infections are maintained through both reproduction on the host’s skin and reinfection from sources outside of the host. To investigate the importance of external reinfection on pathogen burden, we infected captive-bred individuals of the highly susceptible Panamanian Golden Frog, Atelopus glyphus, and wild-caught glass frogs, Espadarana prosoblepon, with Bd. We housed the animals in one of three treatments: individually, in heterospecific pairs, and in conspecific pairs. For 8 weeks, we measured the Bd load and shedding rate of all frogs. We found that Atelopus had high rates of increase in both Bd load and shedding rate, but pathogen growth rates did not differ among treatments. The infection intensity of Espadarana co-housed with Atelopus was indistinguishable from those housed singly and those in conspecific pairs, despite being exposed to a large external source of Bd zoospores. Our results indicate that Bd load in both species is driven by pathogen replication within an individual, with reinfection from outside the host contributing little to the amplification of host fungal load.     

Whether the Weather Drives Patterns of Endemic Amphibian Chytridiomycosis: A Pathogen Proliferation Approach

The pandemic amphibian disease chytridiomycosis often exhibits strong seasonality in both prevalence and disease-associated mortality once it becomes endemic. One hypothesis that could explain this temporal pattern is that simple weather-driven pathogen proliferation (population growth) is a major driver of chytridiomycosis disease dynamics. Despite various elaborations of this hypothesis in the literature for explaining amphibian declines (e.g., the chytrid thermal-optimum hypothesis) it has not been formally tested on infection patterns in the wild. In this study we developed a simple process-based model to simulate the growth of the pathogen Batrachochytrium dendrobatidis (Bd) under varying weather conditions to provide an a priori test of a weather-linked pathogen proliferation hypothesis for endemic chytridiomycosis. We found strong support for several predictions of the proliferation hypothesis when applied to our model species, Litoria pearsoniana, sampled across multiple sites and years: the weather-driven simulations of pathogen growth potential (represented as a growth index in the 30 days prior to sampling; GI₃₀) were positively related to both the prevalence and intensity of Bd infections, which were themselves strongly and positively correlated. In addition, a machine-learning classifier achieved ∼72% success in classifying positive qPCR results when utilising just three informative predictors 1) GI₃₀, 2) frog body size and 3) rain on the day of sampling. Hence, while intrinsic traits of the individuals sampled (species, size, sex) and nuisance sampling variables (rainfall when sampling) influenced infection patterns obtained when sampling via qPCR, our results also strongly suggest that weather-linked pathogen proliferation plays a key role in the infection dynamics of endemic chytridiomycosis in our study system. Predictive applications of the model include surveillance design, outbreak preparedness and response, climate change scenario modelling and the interpretation of historical patterns of amphibian decline.     

Assessing host extinction risk following exposure to Batrachochytrium dendrobatidis

Wildlife diseases are increasingly recognized as a major threat to biodiversity. Chytridiomycosis is an emerging infectious disease of amphibians caused by the fungus Batrachochytrium dendrobatidis (Bd). Using a mathematical model and simulations, we study its effects on a generic riparian host population with a tadpole and adult life stage. An analytical expression for the basic reproduction quotient, Q_o, of the pathogen is derived. By sampling the entire relevant parameter space, we perform a statistical assessment of the importance of all considered parameters in determining the risk of host extinction, upon exposure to Bd. We find that Q_o not only gives a condition for the initial invasion of the fungus, but is in fact the best predictor for host extinction. We also show that the role of tadpoles, which in some species tolerate infections, is ambivalent. While tolerant tadpoles may provide a reservoir for the fungus, thus facilitating its persistence or even amplifying its outbreaks, they can also act as a rescue buffer for a stressed host population. Our results have important implications for amphibian conservation efforts.     

Host persistence or extinction from emerging infectious disease: insights from white-nose syndrome in endemic and invading regions

Predicting species'' fates following the introduction of a novel pathogen is a significant and growing problem in conservation. Comparing disease dynamics between introduced and endemic regions can offer insight into which naive hosts will persist or go extinct, with disease acting as a filter on host communities. We examined four hypothesized mechanisms for host–pathogen persistence by comparing host infection patterns and environmental reservoirs for Pseudogymnoascus destructans (the causative agent of white-nose syndrome) in Asia, an endemic region, and North America, where the pathogen has recently invaded. Although colony sizes of bats and hibernacula temperatures were very similar, both infection prevalence and fungal loads were much lower on bats and in the environment in Asia than North America. These results indicate that transmission intensity and pathogen growth are lower in Asia, likely due to higher host resistance to pathogen growth in this endemic region, and not due to host tolerance, lower transmission due to smaller populations, or lower environmentally driven pathogen growth rate. Disease filtering also appears to be favouring initially resistant species in North America. More broadly, determining the mechanisms allowing species persistence in endemic regions can help identify species at greater risk of extinction in introduced regions, and determine the consequences for disease dynamics and host–pathogen coevolution.     

Thermoregulatory behaviour affects prevalence of chytrid fungal infection in a wild population of Panamanian golden frogs

Predicting how climate change will affect disease dynamics requires an understanding of how the environment affects host–pathogen interactions. For amphibians, global declines and extinctions have been linked to a pathogenic chytrid fungus, Batrachochytrium dendrobatidis. Using a combination of body temperature measurements and disease assays conducted before and after the arrival of B. dendrobatidis, this study tested the hypothesis that body temperature affects the prevalence of infection in a wild population of Panamanian golden frogs (Atelopus zeteki). The timing of first detection of the fungus was consistent with that of a wave of epidemic infections spreading south and eastward through Central America. During the epidemic, many golden frogs modified their thermoregulatory behaviour, raising body temperatures above their normal set point. Odds of infection decreased with increasing body temperature, demonstrating that even slight environmental or behavioural changes have the potential to affect an individual''s vulnerability to infection. The thermal dependency of the relationship between B. dendrobatidis and its amphibian hosts demonstrates how the progression of an epidemic can be influenced by complex interactions between host and pathogen phenotypes and the environments in which they are found.