Fungal infection alters the selection,dispersal and drift processes structuring the amphibian skin microbiome

Symbiotic microbial communities are important for host health, but the processes shaping these communities are poorly understood. Understanding how community assembly processes jointly affect microbial community composition is limited because inflexible community models rely on rejecting dispersal and drift before considering selection. We developed a flexible community assembly model based on neutral theory to ask: How do dispersal, drift and selection concurrently affect the microbiome across environmental gradients? We applied this approach to examine how a fungal pathogen affected the assembly processes structuring the amphibian skin microbiome. We found that the rejection of neutrality for the amphibian microbiome across a fungal gradient was not strictly due to selection processes, but was also a result of species‐specific changes in dispersal and drift. Our modelling framework brings the qualitative recognition that niche and neutral processes jointly structure microbiomes into quantitative focus, allowing for improved predictions of microbial community turnover across environmental gradients.     

Incapacitating effects of fungal coinfection in a novel pathogen system

As globalization lowers geographic barriers to movement, coinfection with novel and enzootic pathogens is increasingly likely. Novel and enzootic pathogens can interact synergistically or antagonistically, leading to increased or decreased disease severity. Here we examine host immune responses to coinfection with two closely related fungal pathogens: Batrachochytrium dendrobatidis (Bd) and Batrachochytrium salamandrivorans (Bsal). Both pathogens have had detrimental effects on amphibian populations, with Bd now largely enzootic, while Bsal is currently spreading and causing epizootics. Recent experimental work revealed that newts coinfected with Bd and Bsal had significantly higher mortality than those infected with either pathogen alone. Here we characterize host immunogenomic responses to chytrid coinfection relative to single infection. Across several classes of immune genes including pattern recognition receptors, cytokines, and MHC, coinfected host gene expression was weakly upregulated or comparable to that seen in single Bd infection, but significantly decreased when compared to Bsal infection. Combined with strong complement pathway downregulation and keratin upregulation, these results indicate that coinfection with Bd and Bsal compromises immune responses active against Bsal alone. As Bsal continues to invade naïve habitats where Bd is enzootic, coinfection will be increasingly common. If other Bd‐susceptible species in the region have similar responses, interactions between the two pathogens could cause severe population and community‐level declines.     

Daphnia predation on the amphibian chytrid fungus and its impacts on disease risk in tadpoles

Direct predation upon parasites has the potential to reduce infection in host populations. For example, the fungal parasite of amphibians, Batrachochytrium dendrobatidis (Bd), is commonly transmitted through a free‐swimming zoospore stage that may be vulnerable to predation. Potential predators of Bd include freshwater zooplankton that graze on organisms in the water column. We tested the ability of two species of freshwater crustacean (Daphnia magna and D. dentifera) to consume Bd and to reduce Bd density in water and infection in tadpoles. In a series of laboratory experiments, we allowed Daphnia to graze in water containing Bd while manipulating Daphnia densities, Daphnia species identity, grazing periods and concentrations of suspended algae (Ankistrodesmus falcatus). We then exposed tadpoles to the grazed water. We found that high densities of D. magna reduced the amount of Bd detected in water, leading to a reduction in the proportion of tadpoles that became infected. Daphnia dentifera, a smaller species of Daphnia, also reduced Bd in water samples, but did not have an effect on tadpole infection. We also found that algae affected Bd in complex ways. When Daphnia were absent, less Bd was detected in water and tadpole samples when concentrations of algae were higher, indicating a direct negative effect of algae on Bd. When Daphnia were present, however, the amount of Bd detected in water samples showed the opposite trend, with less Bd when densities of algae were lower. Our results indicate that Daphnia can reduce Bd levels in water and infection in tadpoles, but these effects vary with species, algal concentration, and Daphnia density. Therefore, the ability of predators to consume parasites and reduce infection is likely to vary depending on ecological context.     

Diverse genotypes of the amphibian‐killing fungus produce distinct phenotypes through plastic responses to temperature

Phenotypes are the target of selection and affect the ability of organisms to persist in variable environments. Phenotypes can be influenced directly by genes and/or by phenotypic plasticity. The amphibian‐killing fungus Batrachochytrium dendrobatidis (Bd) has a global distribution, unusually broad host range, and high genetic diversity. Phenotypic plasticity may be an important process that allows this pathogen to infect hundreds of species in diverse environments. We quantified phenotypic variation of nine Bd genotypes from two Bd lineages (Global Pandemic Lineage [GPL] and Brazil) and a hybrid (GPL‐Brazil) grown at three temperatures (12, 18 and 24°C). We measured five functional traits including two morphological traits (zoospore and zoosporangium sizes) and three life history traits (carrying capacity, time to fastest growth and exponential growth rate) in a phylogenetic framework. Temperature caused highly plastic responses within each genotype, with all Bd genotypes showing phenotypic plasticity in at least three traits. Among genotypes, Bd generally showed the same direction of plastic response to temperature: larger zoosporangia, higher carrying capacity, longer time to fastest growth and slower exponential growth at lower temperatures. The exception was zoospore size, which was highly variable. Our findings indicate that Bd genotypes have evolved novel phenotypes through plastic responses to temperature over very short timescales. High phenotypic variability likely extends to other traits and may facilitate the large host range and rapid spread of Bd.     

Impacts of thermal mismatches on chytrid fungus Batrachochytrium dendrobatidis prevalence are moderated by life stage,body size,elevation and latitude

Global climate change is increasing the frequency of unpredictable weather conditions; however, it remains unclear how species‐level and geographic factors, including body size and latitude, moderate impacts of unusually warm or cool temperatures on disease. Because larger and lower‐latitude hosts generally have slower acclimation times than smaller and higher‐latitude hosts, we hypothesised that their disease susceptibility increases under ‘thermal mismatches’ or differences between baseline climate and the temperature during surveying for disease. Here, we examined how thermal mismatches interact with body size, life stage, habitat, latitude, elevation, phylogeny and International Union for Conservation of Nature (IUCN) conservation status to predict infection prevalence of the chytrid fungus Batrachochytrium dendrobatidis (Bd) in a global analysis of 32 291 amphibian hosts. As hypothesised, we found that the susceptibility of larger hosts and hosts from lower latitudes to Bd was influenced by thermal mismatches. Furthermore, hosts of conservation concern were more susceptible than others following thermal mismatches, suggesting that thermal mismatches might have contributed to recent amphibian declines.     

Unlocking the story in the swab: A new genotyping assay for the amphibian chytrid fungus Batrachochytrium dendrobatidis

One of the most devastating emerging pathogens of wildlife is the chytrid fungus, Batrachochytrium dendrobatidis (Bd), which affects hundreds of amphibian species around the world. Genomic data from pure Bd cultures have advanced our understanding of Bd phylogenetics, genomic architecture and mechanisms of virulence. However, pure cultures are laborious to obtain and whole‐genome sequencing is comparatively expensive, so relatively few isolates have been genetically characterized. Thus, we still know little about the genetic diversity of Bd in natural systems. The most common noninvasive method of sampling Bd from natural populations is to swab amphibian skin. Hundreds of thousands of swabs have been collected from amphibians around the world, but Bd DNA collected via swabs is often low in quality and/or quantity. In this study, we developed a custom Bd genotyping assay using the Fluidigm Access Array platform to amplify 192 carefully selected regions of the Bd genome. We obtained robust sequence data for pure Bd cultures and field‐collected skin swabs. This new assay has the power to accurately discriminate among the major Bd clades, recovering the basic tree topology previously revealed using whole‐genome data. Additionally, we established a critical value for initial Bd load for swab samples (150 Bd genomic equivalents) above which our assay performs well. By leveraging advances in microfluidic multiplex PCR technology and the globally distributed resource of amphibian swab samples, noninvasive skin swabs can now be used to address critical spatial and temporal questions about Bd and its effects on declining amphibian populations.     

The thermal mismatch hypothesis explains host susceptibility to an emerging infectious disease

Parasites typically have broader thermal limits than hosts, so large performance gaps between pathogens and their cold‐ and warm‐adapted hosts should occur at relatively warm and cold temperatures, respectively. We tested this thermal mismatch hypothesis by quantifying the temperature‐dependent susceptibility of cold‐ and warm‐adapted amphibian species to the fungal pathogen Batrachochytrium dendrobatidis (Bd) using laboratory experiments and field prevalence estimates from 15 410 individuals in 598 populations. In both the laboratory and field, we found that the greatest susceptibility of cold‐ and warm‐adapted hosts occurred at relatively warm and cool temperatures, respectively, providing support for the thermal mismatch hypothesis. Our results suggest that as climate change shifts hosts away from their optimal temperatures, the probability of increased host susceptibility to infectious disease might increase, but the effect will depend on the host species and the direction of the climate shift. Our findings help explain the tremendous variation in species responses to Bd across climates and spatial, temporal and species‐level variation in disease outbreaks associated with extreme weather events that are becoming more common with climate change.     

Underestimated ranges and overlooked refuges from amphibian chytridiomycosis

Aim

Accurately documenting and predicting declines and shifts in species’ distributions is fundamental for implementing effective conservation strategies and directing future research; species distribution models (SDM) have become a powerful tool for such work. Nevertheless, much of the data used to create these models are opportunistic and often violate some of their basic assumptions. We use amphibian declines and extinctions linked to the fungus Batrachochytrium dendrobatidis (Bd) to examine how sampling biases in data collection can affect what we know of this disease and its effect on amphibians in the wild.

Location

Queensland, Australia.

Methods

We developed a distribution model for Bd incorporating known locality records for Bd and a subset of climatic variables that should correctly characterize its distribution. We tested this (original) model with additional surveys, recorded new Bd observations in novel environments and reran the distribution model. We then investigated the difference between the original and new models, and used frog abundance and infection status data from two of these new localities to look at the susceptibility of the torrent frog Litoria nannotis to chytridiomycosis.

Results

While largely correct, the original SDM underestimated the distribution of Bd; sampling in ‘unsuitable’ drier environments discovered abundant populations of susceptible frogs with pathogen prevalences of up to 100%. The validation surveys further uncovered a new population of the frog Litoria lorica coexisting with the pathogen; this species was previously believed to be an extinct rain forest endemic.

Main conclusion

Our results indicate that SDMs constructed using opportunistically collected data can be biased if species are not at equilibrium with their environment or because environmental gradients have not been adequately sampled. For disease ecology, the better estimations of pathogen distribution may lead to the discovery of new populations persisting at the edge of their range, which has important implications for the conservation of species threatened by chytridiomycosis.

     

Amphibian‐killing chytrid in Brazil comprises both locally endemic and globally expanding populations

Chytridiomycosis, caused by the fungus Batrachochytrium dendrobatidis (Bd), is the emerging infectious disease implicated in recent population declines and extinctions of amphibian species worldwide. Bd strains from regions of disease‐associated amphibian decline to date have all belonged to a single, hypervirulent clonal genotype (Bd‐GPL). However, earlier studies in the Atlantic Forest of southeastern Brazil detected a novel, putatively enzootic lineage (Bd‐Brazil), and indicated hybridization between Bd‐GPL and Bd‐Brazil. Here, we characterize the spatial distribution and population history of these sympatric lineages in the Brazilian Atlantic Forest. To investigate the genetic structure of Bd in this region, we collected and genotyped Bd strains along a 2400‐km transect of the Atlantic Forest. Bd‐Brazil genotypes were restricted to a narrow geographic range in the southern Atlantic Forest, while Bd‐GPL strains were widespread and largely geographically unstructured. Bd population genetics in this region support the hypothesis that the recently discovered Brazilian lineage is enzootic in the Atlantic Forest of Brazil and that Bd‐GPL is a more recently expanded invasive. We collected additional hybrid isolates that demonstrate the recurrence of hybridization between panzootic and enzootic lineages, thereby confirming the existence of a hybrid zone in the Serra da Graciosa mountain range of Paraná State. Our field observations suggest that Bd‐GPL may be more infective towards native Brazilian amphibians, and potentially more effective at dispersing across a fragmented landscape. We also provide further evidence of pathogen translocations mediated by the Brazilian ranaculture industry with implications for regulations and policies on global amphibian trade.     

Infection with Batrachochytrium dendrobatidis is common in tropical lowland habitats: Implications for amphibian conservation

Numerous species of amphibians declined in Central America during the 1980s and 1990s. These declines mostly affected highland stream amphibians and have been primarily linked to chytridiomycosis, a deadly disease caused by the chytrid fungus Batrachochytrium dendrobatidis (Bd). Since then, the majority of field studies on Bd in the Tropics have been conducted in midland and highland environments (>800 m) mainly because the environmental conditions of mountain ranges match the range of ideal abiotic conditions for Bd in the laboratory. This unbalanced sampling has led researchers to largely overlook host–pathogen dynamics in lowlands, where other amphibian species declined during the same period. We conducted a survey testing for Bd in 47 species (n = 348) in four lowland sites in Costa Rica to identify local host–pathogen dynamics and to describe the abiotic environment of these sites. We detected Bd in three sampling sites and 70% of the surveyed species. We found evidence that lowland study sites exhibit enzootic dynamics with low infection intensity and moderate to high prevalence (55% overall prevalence). Additionally, we found evidence that every study site represents an independent climatic zone, where local climatic differences may explain variations in Bd disease dynamics. We recommend more detection surveys across lowlands and other sites that have been historically considered unsuitable for Bd occurrence. These data can be used to identify sites for potential disease outbreaks and amphibian rediscoveries.     

Species distribution models predict the geographic expansion of an enzootic amphibian pathogen

Globally, numerous amphibian species have declined due to the introduction of the chytrid fungus Batrachochytrium dendrobatidis (Bd). However, the understanding of the spatiotemporal dynamics remains incomplete. Therefore, estimating the current geographic distribution of Bd is urgently needed, especially in countries like Costa Rica, where susceptible species are still recovering from Bd‐driven declines. We conducted model tuning and spatial analysis to compare the habitat suitability for epizootic and enzootic Bd in Costa Rica and to identify data‐deficient regions, opportunistic sampling, and Bd hotspots. Our dataset combined two methods of detection (histology and PCR methods) for a total of 451 Bd‐positive records from 34 localities. We found that the distribution of enzootic Bd in Costa Rica increased 60% since previous estimates in the early 2000s and extended to highlands and dry lowlands that were considered unsuitable for Bd. We also found that Bd is common across protected lands (80%) and within the herpetological provinces containing the highest amphibian richness and endemism in Costa Rica. Opportunistic sampling of Bd has focused on sites where epizootics occurred with the strongest intensity, leading to deficient or absent sampling across the Talamanca Range, the Nicoya Peninsula, and the northern lowlands. Our results showed that PCR increased the power of Bd detection in lowlands and favored the identification of Bd hotspots across the Caribbean side of Costa Rica. Our results add to the understanding of disease spread during enzootics and can be used to identify new hotspots for disease to mitigate future outbreaks of this pathogen. Abstract in Spanish is available with online material     

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.     

Elevated temperatures alter an ant‐aphid mutualism

Ant‐hemipteran mutualisms are keystone interactions that can be variously affected by warming: these mutualisms can be strengthened or weakened, or the species can transition to new mutualist partners. We examined the effects of elevated temperatures on an ant‐aphid mutualism in the subalpine zone of the Rocky Mountains in Colorado, USA. In this system, inflorescences of the host plant, Ligusticum porteri Coult. & Rose (Apiaceae), are colonized by the ant‐tended aphid Aphis asclepiadis Fitch or less frequently by the non‐ant tended aphid Cavariella aegopodii (Scopoli) (both Hemiptera: Aphididae). Using an 8‐year observational study, we tested for two key mechanisms by which ant‐hemipteran mutualisms may be altered by climate change: shifts in species identity and phenological mismatch. Whereas the aphid species colonizing the host plant is not changing in response to year‐to‐year variation in temperature, we found evidence that a phenological mismatch between ants and aphids could occur. In warmer years, colonization of host plant inflorescences by ants is decreased, whereas for A. asclepiadis aphids, host plant colonization is mostly responsive to date of snowmelt. We also experimentally established A. asclepiadis colonies on replicate host plants at ambient and elevated temperatures. Ant abundance did not differ between aphid colonies at ambient vs. elevated temperatures, but ants were less likely to engage in tending behaviors on aphid colonies at elevated temperatures. Sugar composition of aphid honeydew was also altered by experimental warming. Despite reduced tending by ants, aphid colonies at elevated temperatures had fewer intraguild predators. Altogether, our results suggest that higher temperatures may disrupt this ant‐aphid mutualism through both phenological mismatch and by altering benefits exchanged in the interaction.     

Precipitation Constrains Amphibian Chytrid Fungus Infection Rates in a Terrestrial Frog Assemblage in Jamaica,West Indies

We model Batrachochytrium dendrobatidis (Bd) infection rates in Jamaican frogs—one of the most threatened amphibian fauna in the world. The majority of species we surveyed were terrestrial direct‐developing frogs or frogs that breed in tank bromeliads, rather than those that use permanent water bodies to breed. Thus, we were able to investigate the climatic correlates of Bd infection in a frog assemblage that does not rely on permanent water bodies. We sampled frogs for Bd across all of the major habitat types on the island, used machine learning algorithms to identify climatic variables that are correlated with infection rates, and extrapolated infection rates across the island. We compared the effectiveness of the machine learning algorithms for species distribution modeling in the context of our study, and found that infection rate rose quickly with precipitation in the driest month. Infection rates also increased with mean temperature in the warmest quarter until 22 °C, and remained relatively level thereafter. Both of these results are in accordance with previous studies of the physiology of Bd. Based on our environmental results, we suggest that frogs occupying high‐precipitation habitats with cool rainy‐season temperatures, though zcurrently experiencing low frequencies of infection, may experience an increase in infection rates as global warming increases temperatures in their habitat.     

Mitigating amphibian chytridiomycosis with bioaugmentation: characteristics of effective probiotics and strategies for their selection and use

Probiotic therapy through bioaugmentation is a feasible disease mitigation strategy based on growing evidence that microbes contribute to host defences of plants and animals. Amphibians are currently threatened by the rapid global spread of the pathogen, Batrachochytrium dendrobatidis (Bd), which causes the disease chytridiomycosis. Bioaugmentation of locally occurring protective bacteria on amphibians has mitigated this disease effectively in laboratory trials and one recent field trial. Areas still naïve to Bd provide an opportunity for conservationists to proactively implement probiotic strategies to prevent further amphibian declines. In areas where Bd is endemic, bioaugmentation can facilitate repatriation of susceptible amphibians currently maintained in assurance colonies. Here, we synthesise the current research in amphibian microbial ecology and bioaugmentation to identify characteristics of effective probiotics in relation to their interactions with Bd, their host, other resident microbes and the environment. To target at‐risk species and amphibian communities, we develop sampling strategies and filtering protocols that result in probiotics that inhibit Bd under ecologically relevant conditions and persist on susceptible amphibians. This filtering tool can be used proactively to guide amphibian disease mitigation and can be extended to other taxa threatened by emerging infectious diseases.     

Batesian mimics influence the evolution of conspicuousness in an aposematic salamander

Conspicuousness, or having high contrast relative to the surrounding background, is a common feature of unpalatable species. Several hypotheses have been proposed to explain the occurrence of conspicuousness, and while most involve the role of conspicuousness as a direct signal of unpalatability to potential predators, one hypothesis suggests that exaggerated conspicuousness may evolve in unpalatable species to reduce predator confusion with palatable species (potential Batesian mimics). This hypothesis of antagonistic coevolution between palatable and unpalatable species hinges on the ‘cost of conspicuousness’, in which conspicuousness increases the likelihood of predation more in palatable species than in unpalatable species. Under this mimicry scenario, four patterns are expected: (i) mimics will more closely resemble local models than models from other localities, (ii) there will be a positive relationship between mimic and model conspicuousness, (iii) models will be more conspicuous in the presence of mimics, and (iv) when models and mimics differ in conspicuousness, mimics will be less conspicuous than models. We tested these predictions in the salamander mimicry system involving Notophthalmus viridescens (model) and one colour morph of Plethodon cinereus (mimic). All predictions were supported, indicating that selection for Batesian mimicry not only influences the evolution of mimics, but also the evolution of the models they resemble. These findings indicate that mimicry plays a large role in the evolution of model warning signals, particularly influencing the evolution of conspicuousness.     

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.