Cool temperatures reduce antifungal activity of symbiotic bacteria of threatened amphibians--implications for disease management and patterns of decline

Effects of temperature on the interaction between amphibian skin bacteria and Batrachochytrium dendrobatidis

Symbiotic relationships between animals and microbes are important for a range of functions, from digestion to protection from pathogens. However, the impact of temperature variation on these animal-microbe interactions remains poorly understood. Amphibians have experienced population declines and even extinctions on a global scale due to chytridiomycosis, a disease caused by chytrid fungi in the genus Batrachochytrium. Variation in susceptibility to this disease exists within and among host species. While the mechanisms generating differences in host susceptibility remain elusive, differences in immune system components, as well as variation in host and environmental temperatures, have been associated with this variation. The symbiotic cutaneous bacteria of amphibians are another potential cause for variation in susceptibility to chytridiomycosis, with some bacterial species producing antifungal metabolites that prevent the growth of Bd. The growth of both Bd and bacteria are affected by temperature, and thus we hypothesized that amphibian skin bacteria may be more effective at preventing Bd growth at certain temperatures. To test this, we collected bacteria from the skins of frogs, harvested the metabolites they produced when grown at three different temperatures, and then grew Bd in the presence of those metabolites under those same three temperatures in a three-by-three fully crossed design. We found that both the temperature at which cutaneous bacteria were grown (and metabolites produced) as well as the temperature at which Bd is grown can impact the ability of cutaneous bacteria to inhibit the growth of Bd. While some bacterial isolates showed the ability to inhibit Bd growth across multiple temperature treatments, no isolate was found to be inhibitive across all combinations of bacterial incubation or Bd challenge temperatures, suggesting that temperature affects both the metabolites produced and the effectiveness of those metabolites against the Bd pathogen. These findings move us closer to a mechanistic understanding of why chytridiomycosis outbreaks and related amphibian declines are often limited to certain climates and seasons.

Metabolites of Xenorhabdus bacteria are potent candidates for mitigating amphibian chytridiomycosis

Chytridiomycosis, caused by the chytrid fungus Batrachochytrium dendrobatidis (Bd), has caused extreme losses in amphibian biodiversity. Finding bacteria that produce metabolites with antifungal properties may turn out to be invaluable in the fight against this devastating disease. The entomopathogenic bacteria, Xenorhabdus szentirmaii and X. budapestensis produce secondary metabolites that are effective against a wide range of fungal plant pathogens. To assess whether they may also be effective against Bd, we extracted cell-free culture media (CFCM) from liquid cultures of X. szentirmaii and X. budapestensis and tested their ability to inhibit Bd growth in vitro. As a second step, using juvenile common toads (Bufo bufo) experimentally infected with Bd we also tested the in vivo antifungal efficacy of X. szentirmaii CFCM diluted to 2 and 10% (v/v), while also assessing possible malign side effects on amphibians. Results of the in vitro experiment documented highly effective growth inhibition by CFCMs of both Xenorhabdus species. The in vivo experiment showed that treatment with CFCM of X. szentirmaii applied at a dilution of 10% resulted in infection intensities reduced by ca. 73% compared to controls and to juvenile toads treated with CFCM applied at a dilution of 2%. At the same time, we detected no negative side effects of treatment with CFCM on toad survival and development. Our results clearly support the idea that metabolites of X. szentirmaii, and perhaps of several other Xenorhabdus species as well, may prove highly useful for the treatment of Bd infected amphibians.

Spatiotemporal and ontogenetic variation, microbial selection, and predicted Bd-inhibitory function in the skin-associated microbiome of a Rocky Mountain amphibian

Host-associated microbiomes play important roles in host health and pathogen defense. In amphibians, the skin-associated microbiota can contribute to innate immunity with potential implications for disease management. Few studies have examined season-long temporal variation in the amphibian skin-associated microbiome, and the interactions between bacteria and fungi on amphibian skin remain poorly understood. We characterize season-long temporal variation in the skin-associated microbiome of the western tiger salamander (Ambystoma mavortium) for both bacteria and fungi between sites and across salamander life stages. Two hundred seven skin-associated microbiome samples were collected from salamanders at two Rocky Mountain lakes throughout the summer and fall of 2018, and 127 additional microbiome samples were collected from lake water and lake substrate. We used 16S rRNA and ITS amplicon sequencing with Bayesian Dirichlet-multinomial regression to estimate the relative abundances of bacterial and fungal taxa, test for differential abundance, examine microbial selection, and derive alpha diversity. We predicted the ability of bacterial communities to inhibit the amphibian chytrid fungus Batrachochytrium dendrobatidis (Bd), a cutaneous fungal pathogen, using stochastic character mapping and a database of Bd-inhibitory bacterial isolates. For both bacteria and fungi, we observed variation in community composition through time, between sites, and with salamander age and life stage. We further found that temporal trends in community composition were specific to each combination of salamander age, life stage, and lake. We found salamander skin to be selective for microbes, with many taxa disproportionately represented relative to the environment. Salamander skin appeared to select for predicted Bd-inhibitory bacteria, and we found a negative relationship between the relative abundances of predicted Bd-inhibitory bacteria and Bd. We hope these findings will assist in the conservation of amphibian species threatened by chytridiomycosis and other emerging diseases.

Short-term continuous and pulse Pb exposure causes negative effects on skin histomorphological structure and bacterial composition of adult Pelophylax nigromaculatus

Lead (Pb) is a toxic heavy metal often present in the environment as a pulse in water. Traditional toxicity tests are usually carried out under conditions of continuous concentration, without considering the impact of pulse exposure on aquatic organisms. This study aimed to evaluate the effects of short-term continuous and pulse Pb exposures on the skin bacteria and histomorphological structure of Pelophylax nigromaculatus. Results showed that compared to the control (CON) and Pb continuous exposure group (CEPb), the Pb pulse exposure group (PEPb) showed the smallest size of granular glands, which would interfere with the permeability and secretory function of skin, making the individual more sensitive to external pollution. Lead exposure significantly changed the composition and diversity of skin bacteria. Compared to the CON and CEPb groups, the PEPb group showed a significant increase in the abundance of harmful bacteria (e.g., Bacteroidetes and Chryseobacterium) and a decrease in the abundance of beneficial bacteria (e.g., Pseudomonas). PICRUSt software showed that there were differences in the metabolic pathway of skin bacteria among the three groups (CON, CEPb, and PEPb). Overall, this study indicates that Pb pulse exposure can aggravate the toxicity of Pb for frog skin, providing a new framework for simulating short-term heavy metal exposure in the context of frog health.

Comparative Analysis of Skin Bacterial Diversity and Its Potential Antifungal Function Between Desert and Pine Forest Populations of Boreal Toads Anaxyrus boreas

The skin microbiome in amphibians has gained a lot of attention as some of its members play a protective role against pathogens such as the fungus Batrachochytrium dendrobatidis (Bd). The composition of skin bacterial communities has been suggested as one of the factors explaining differences in susceptibility to Bd among amphibian species and populations. The boreal toad Anaxyrus boreas is known to be susceptible to Bd, and severe population declines in its southeastern range have been documented. However, throughout A. boreas distribution, populations present differences in susceptibility to Bd infections which may be associated with differences in skin microbial diversity. This study compared the skin bacterial diversity and Bd infection levels of A. boreas in one desert population and one pine forest population from Baja California, Mexico. We found that desert and pine forest toad populations exhibit differences in skin bacterial community structure but show similar Bd infection levels. Using a predictive method, we found that the abundance of bacteria with potential Bd-inhibitory properties differed between uninfected and infected individuals but not between populations. Our data suggest that several bacteria in the skin community may be offering protection from Bd infections in these A. boreas populations. This study provides foundational evidence for future studies seeking to understand the skin-microbial variation among boreal toads' populations and its relation with Bd susceptibility.

Inhibition of Batrachochytrium dendrobatidis Infection by Skin Bacterial Communities in Wild Amphibian Populations

Skin-associated bacteria are known to inhibit infection by the fungal pathogen Batrachochytrium dendrobatidis (Bd) in amphibians. It has also been postulated that skin-associated bacterial community is related to Bd infection intensity. However, our understanding of host microbial dynamics and their importance in regulating Bd intensity is limited. We analyzed Bd infection and skin-associated bacteria from two amphibian species, the salamander Ambystoma rivulare and the frog Lithobates spectabilis that co-occurred in a tropical high-altitude site in central Mexico. Sixty-three percent of sampled salamander individuals and 80% of frog individuals tested positive for Bd. Overall, we registered 622 skin-associated bacterial genera, from which 73 are known to have Bd inhibitory effects. These inhibitory taxa represented a relative abundance of 50% in relation to total relative bacterial abundance. Our results indicated that, although sharing some bacterial taxa, bacterial community from the skin of both species was different in taxonomic composition and in relative abundance. Pseudomonas spp. and Stenotrophomonas spp. were among the five most abundant bacterial taxa of both species. Both bacterial taxa inhibit Bd infection. We detected that bacterial richness and relative abundance of inhibitory Bd bacteria were negatively related to intensity of Bd infection independent of species and seasons. Despite the high Bd prevalence in both host species, no dead or sick individuals were registered during field surveys. The relatively low levels of Bd load apparently do not compromise survival of host species. Therefore, our results suggested that individuals analyzed were able to survive and thrive under a dynamic relation with enzootic infections of Bd and their microbiota.

TEMPERATURE AS A DRIVER OF THE PATHOGENICITY AND VIRULENCE OF AMPHIBIAN CHYTRID FUNGUS BATRACHOCHYTRIUM DENDROBATIDIS: A SYSTEMATIC REVIEW

Chytridiomycosis, caused by the chytrid fungus Batrachochytrium dendrobatidis (Bd), is a leading cause of global amphibian declines. Severe infections with Bd can lead to cardiac arrest, and mass deaths during epidemics have been reported. Temperature, pH, salinity, and moisture are important determinants of the survival, growth, reproduction, and pathogenicity of Bd, as well as its effect on amphibian populations. Here, we synthesize current knowledge on the role of temperature as a driver of the pathogenicity and virulence of Bd to better understand the effects of temperature on amphibian defense mechanisms against infection. This review advises on research direction and management approaches to benefit amphibian populations affected by Bd. We conclude by offering guidelines for four levels of temperature monitoring in amphibian field studies to improve consistency between studies: regional climate, habitat, microhabitat, and amphibian host.

Exposure to a fungal pathogen increases the critical thermal minimum of two frog species

The ability of an organism to tolerate seasonal temperature changes, such as extremely cold temperatures during the winter, can be influenced by their pathogens. We tested how exposure to a virulent fungal pathogen, Batrachochytrium dendrobatidis (Bd), affected the critical thermal minimum (CTmin) of two frog species, Hyla versicolor (gray treefrog) and Lithobates palustris (pickerel frog). The CTmin is the minimum thermal performance point of an organism, which we estimated via righting response trials. For both frog species, we compared the righting response of Bd-exposed and Bd-unexposed individuals in either a constant (15ºC) environment or with decreasing temperatures (-1°C/2.5 min) starting from 15°C. The CTmin for both species was higher for Bd-exposed frogs than unexposed frogs, and the CTmin of H. versicolor was higher than L. palustris. We also found that Bd-exposed frogs of both species righted themselves significantly fewer times in both decreasing and constant temperature trials. Our findings show that pathogen exposure can reduce cold tolerance and limit the thermal performance range of hosts, which may lead to increased overwintering mortality.

Co-infecting pathogen lineages have additive effects on host bacterial communities

Amphibian skin bacteria may confer protection against the fungus Batrachochytrium dendrobatidis (Bd), but responses of skin bacteria to different Bd lineages are poorly understood. The global panzootic lineage (Bd-GPL) has caused amphibian declines and extinctions globally. However, other lineages are enzootic (Bd-Asia-2/Brazil). Increased contact rates between Bd-GPL and enzootic lineages via globalization pose unknown consequences for host-microbiome-pathogen dynamics. We conducted a laboratory experiment and used 16S rRNA amplicon-sequencing to assess: (i) whether two lineages (Bd-Asia-2/Brazil and Bd-GPL) and their recombinant, in single and mixed infections, differentially affect amphibian skin bacteria; (ii) and the changes associated with the transition to laboratory conditions. We determined no clear differences in bacterial diversity among Bd treatments, despite differences in infection intensity. However, we observed an additive effect of mixed infections on bacterial alpha diversity and a potentially antagonistic interaction between Bd genotypes. Additionally, observed changes in community composition suggest a higher ability of Bd-GPL to alter skin bacteria. Lastly, we observed a drastic reduction in bacterial diversity and a change in community structure in laboratory conditions. We provide evidence for complex interactions between Bd genotypes and amphibian skin bacteria during coinfections, and expand on the implications of experimental conditions in ecological studies.

Characterizing the culturable surface microbiomes of diverse marine animals

Biofilm-forming bacteria have the potential to contribute to the health, physiology, behavior and ecology of the host and serve as its first line of defense against adverse conditions in the environment. While metabarcoding and metagenomic information furthers our understanding of microbiome composition, fewer studies use cultured samples to study the diverse interactions among the host and its microbiome, as cultured representatives are often lacking. This study examines the surface microbiomes cultured from three shallow-water coral species and two whale species. These unique marine animals place strong selective pressures on their microbial symbionts and contain members under similar environmental and anthropogenic stress. We developed an intense cultivation procedure, utilizing a suite of culture conditions targeting a rich assortment of biofilm-forming microorganisms. We identified 592 microbial isolates contained within 15 bacterial orders representing 50 bacterial genera, and two fungal species. Culturable bacteria from coral and whale samples paralleled taxonomic groups identified in culture-independent surveys, including 29% of all bacterial genera identified in the Megaptera novaeangliae skin microbiome through culture-independent methods. This microbial repository provides raw material and biological input for more nuanced studies which can explore how members of the microbiome both shape their micro-niche and impact host fitness.

Winter is coming-Temperature affects immune defenses and susceptibility to Batrachochytrium salamandrivorans

Environmental temperature is a key factor driving various biological processes, including immune defenses and host-pathogen interactions. Here, we evaluated the effects of environmental temperature on the pathogenicity of the emerging fungal pathogen, Batrachochytrium salamandrivorans (Bsal), using controlled laboratory experiments, and measured components of host immune defense to identify regulating mechanisms. We found that adult and juvenile Notophthalmus viridescens died faster due to Bsal chytridiomycosis at 14°C than at 6 and 22°C. Pathogen replication rates, total available proteins on the skin, and microbiome composition likely drove these relationships. Temperature-dependent skin microbiome composition in our laboratory experiments matched seasonal trends in wild N. viridescens, adding validity to these results. We also found that hydrophobic peptide production after two months post-exposure to Bsal was reduced in infected animals compared to controls, perhaps due to peptide release earlier in infection or impaired granular gland function in diseased animals. Using our temperature-dependent susceptibility results, we performed a geographic analysis that revealed N. viridescens populations in the northeastern United States and southeastern Canada are at greatest risk for Bsal invasion, which shifted risk north compared to previous assessments. Our results indicate that environmental temperature will play a key role in the epidemiology of Bsal and provide evidence that temperature manipulations may be a viable disease management strategy.

In 2010, a new skin-eating fungus, Batrachochytrium salamandrivorans (Bsal), was discovered killing salamanders in the Netherlands. Since then, the pathogen has spread to other European countries. Bsal is believed to be from Asia and is being translocated through the international trade of amphibians. To our knowledge, Bsal has not arrived to North America. As a proactive strategy for disease control, we evaluated how a range of environmental temperatures in North America could affect invasion risk of Bsal into a widely distributed salamander species, the eastern newt (Notophthalmus viridescens). Our results show that northeastern USA, southeastern Canada, and the higher elevations of the Appalachian Mountains have the greatest likelihood of Bsal invasion, when temperature-dependent susceptibility is included in risk analyses. Changes in eastern newt susceptibility to Bsal infection associated with temperature are likely an interaction between pathogen replication rate and host immune defenses, including changes in skin microbiome composition and the host’s ability to produce Bsal-killing proteins on the skin. Our study provides new insights into how latitude, elevation and season can impact the epidemiology of Bsal, and suggests that strategies that manipulate microclimate of newt habitats could be useful in managing Bsal outbreaks and that climate change will impact Bsal invasion probability.

Infection Status as the Basis for Habitat Choices in a Wild Amphibian

AbstractAnimals challenged with disease may select specific habitat conditions that help prevent or reduce infection. Whereas preinfection avoidance of habitats with a high risk of disease exposure has been documented in both captive and free-ranging animals, evidence of switching habitats after infection to support the clearing of the infection is limited to laboratory experiments. The extent to which wild animals proximately modify habitat choices in response to infection status thus remains unclear. We investigated preinfection behavioral avoidance and postinfection habitat switching using wild, radio-tracked boreal toads (Anaxyrus boreas boreas) in a population challenged with Batrachochytrium dendrobatidis (Bd), a pathogenic fungus responsible for a catastrophic panzootic affecting hundreds of amphibian species worldwide. Boreal toads did not preemptively avoid microhabitats with conditions conducive to Bd growth. Infected individuals, however, selected warmer, more open habitats, which were associated with elevated body temperature and the subsequent clearing of infection. Our results suggest that disease can comprise an important selective pressure on animal habitat and space use. Habitat selection models, therefore, may be greatly improved by including variables that quantify infection risk and/or the infection status of individuals through time.

Skin microbiome correlates with bioclimate and Batrachochytrium dendrobatidis infection intensity in Brazil's Atlantic Forest treefrogs

In Brazil’s Atlantic Forest (AF) biodiversity conservation is of key importance since the fungal pathogen Batrachochytrium dendrobatidis (Bd) has led to the rapid loss of amphibian populations here and worldwide. The impact of Bd on amphibians is determined by the host's immune system, of which the skin microbiome is a critical component. The richness and diversity of such cutaneous bacterial communities are known to be shaped by abiotic factors which thus may indirectly modulate host susceptibility to Bd. This study aimed to contribute to understanding the environment-host–pathogen interaction determining skin bacterial communities in 819 treefrogs (Anura: Hylidae and Phyllomedusidae) from 71 species sampled across the AF. We investigated whether abiotic factors influence the bacterial community richness and structure on the amphibian skin. We further tested for an association between skin bacterial community structure and Bd co-occurrence. Our data revealed that temperature, precipitation, and elevation consistently correlate with richness and diversity of the skin microbiome and also predict Bd infection status. Surprisingly, our data suggest a weak but significant positive correlation of Bd infection intensity and bacterial richness. We highlight the prospect of future experimental studies on the impact of changing environmental conditions associated with global change on environment-host–pathogen interactions in the AF.

Mechanisms underlying host persistence following amphibian disease emergence determine appropriate management strategies

Emerging infectious diseases have caused many species declines, changes in communities and even extinctions. There are also many species that persist following devastating declines due to disease. The broad mechanisms that enable host persistence following declines include evolution of resistance or tolerance, changes in immunity and behaviour, compensatory recruitment, pathogen attenuation, environmental refugia, density-dependent transmission and changes in community composition. Here we examine the case of chytridiomycosis, the most important wildlife disease of the past century. We review the full breadth of mechanisms allowing host persistence, and synthesise research on host, pathogen, environmental and community factors driving persistence following chytridiomycosis-related declines and overview the current evidence and the information required to support each mechanism. We found that for most species the mechanisms facilitating persistence have not been identified. We illustrate how the mechanisms that drive long-term host population dynamics determine the most effective conservation management strategies. Therefore, understanding mechanisms of host persistence is important because many species continue to be threatened by disease, some of which will require intervention. The conceptual framework we describe is broadly applicable to other novel disease systems.

Predicting Microbiome Function Across Space Is Confounded by Strain-Level Differences and Functional Redundancy Across Taxa

Variation in the microbiome among individual organisms may play a critical role in the relative susceptibility of those organisms to infection, disease, and death. However, predicting microbiome function is difficult because of spatial and temporal variation in microbial diversity, and taxonomic diversity is not predictive of microbiome functional diversity. Addressing this issue may be particularly important when addressing pandemic diseases, such as the global amphibian die-off associated with Bd. Some of the most important factors in probiotic development for disease treatment are whether bacteria with desired function can be found on native amphibians in the local environment. To address this issue, we isolated, sequenced, and assayed the cutaneous bacterial communities of Plethodon cinereus along a gradient of land use change. Our results suggest that cutaneous community composition, but not overall diversity, change with changes in land use, but this does not correspond to significant change in Bd-inhibitory function. We found that Bd-inhibition is a functionally redundant trait, but that level of inhibition varies over phylogenetic, spatial, and temporal scales. This research provides further evidence for the importance of continued examination of amphibian microbial communities across environmental gradients, including biotic and abiotic interactions, when considering disease dynamics.

Functional Redundancy of Batrachochytrium dendrobatidis Inhibition in Bacterial Communities Isolated from Lithobates clamitans Skin

The cutaneous microbial community can influence the health of amphibians exposed to Batrachochytrium dendrobatidis (Bd), a fungal pathogen that has contributed to recent amphibian declines. Resistance to Bd in amphibian populations is correlated with the presence of anti-Bd cutaneous microbes, which confer disease resistance by inhibiting Bd growth. I aimed to determine if green frogs (Lithobates clamitans), an abundant and widely distributed species in New Jersey, harbored bacteria that inhibit Bd and whether the presence and identity of these microbes varied among sites. I used in vitro challenge assays to determine if bacteria isolated from green frog skin could inhibit or enhance the growth of Bd. I found that green frogs at all sites harbored anti-Bd bacteria. However, there were differences in Bd inhibition capabilities among bacterial isolates identified as the same operational taxonomic unit (OTU), lending support to the idea that phylogenetic relatedness does not always predict Bd inhibition status. Additionally, anti-Bd bacterial richness did not vary by site, but the composition of anti-Bd bacterial taxa was distinct at each site. This suggests that there is functional redundancy of Bd inhibition across unique communities of anti-Bd symbionts found on frogs at different sites. These findings highlight the need to better elucidate the structure-function relationship of microbiomes and their role in disease resistance.

Temperature dependence of parasitic infection and gut bacterial communities in bumble bees

High temperatures (e.g., fever) and gut microbiota can both influence host resistance to infection. However, effects of temperature-driven changes in gut microbiota on resistance to parasites remain unexplored. We examined the temperature dependence of infection and gut bacterial communities in bumble bees infected with the trypanosomatid parasite Crithidia bombi. Infection intensity decreased by over 80% between 21 and 37°C. Temperatures of peak infection were lower than predicted based on parasite growth in vitro, consistent with mismatches in thermal performance curves of hosts, parasites and gut symbionts. Gut bacterial community size and composition exhibited slight but significant, non-linear, and taxon-specific responses to temperature. Abundance of total gut bacteria and of Orbaceae, both negatively correlated with infection in previous studies, were positively correlated with infection here. Prevalence of the bee pathogen-containing family Enterobacteriaceae declined with temperature, suggesting that high temperature may confer protection against diverse gut pathogens. Our results indicate that resistance to infection reflects not only the temperature dependence of host and parasite performance, but also temperature-dependent activity of gut bacteria. The thermal ecology of gut parasite-symbiont interactions may be broadly relevant to infectious disease, both in ectothermic organisms that inhabit changing climates, and in endotherms that exhibit fever-based immunity.

The Skin Microbiota of Eleutherodactylus Frogs: Effects of Host Ecology, Phylogeny, and Local Environment

Amphibian skin microbiota has a potential protective role against diseases. However, the effects of environmental and host factors on symbiotic bacterial communities are not well understood. Caribbean frogs in the genus Eleutherodactylus represent a case of congeneric species that differ in ecological specialization by the process of adaptive radiation. For a small clade of Eleutherodactylus from Puerto Rico, we investigated the role of local environments, host species, and microhabitat in the composition of their skin microbiome. The potential congruence between microbial communities in hosts that are most closely related phylogenetically was also addressed. We hypothesized that the skin microbiota of Eleutherodactylus frogs would be mostly associated to microhabitat use, but also differ according to locality, and to a lesser extent to host species. To test this hypothesis, we swabbed the skin of a total of 98 adult individuals of seven Eleutherodactylus species distributed in two nearby localities in Puerto Rico, and sequenced the V4 region of the 16S rRNA gene. Results showed that locality had the greatest effect on determining skin bacterial communities of amphibian hosts, but this effect was stronger on the composition (based on presence/absence) than on its structure (based on sequence abundance). The most ecologically distinct host, E. cooki, and the generalist E. coqui presented, respectively, the most dissimilar and similar microbiota compared to other hosts. Host phylogeny showed a weak influence on skin microbiota. Results suggest that both local environment and ecological specialization are structuring the skin bacterial community in these Eleutherodactylus species, but that characteristics intrinsic to species may also render unique hosts the ability to maintain distinct microbiotas.

Reservoir frogs: seasonality of Batrachochytrium dendrobatidis infection in robber frogs in Dominica and Montserrat

Emerging infectious diseases are an increasingly important threat to wildlife conservation, with amphibian chytridiomycosis, caused by Batrachochytrium dendrobatidis, the disease most commonly associated with species declines and extinctions. However, some amphibians can be infected with B. dendrobatidis in the absence of disease and can act as reservoirs of the pathogen. We surveyed robber frogs (Eleutherodactylus spp.), potential B. dendrobatidis reservoir species, at three sites on Montserrat, 2011-2013, and on Dominica in 2014, to identify seasonal patterns in B. dendrobatidis infection prevalence and load (B. dendrobatidis genomic equivalents). On Montserrat there was significant seasonality in B. dendrobatidis prevalence and B. dendrobatidis load, both of which were correlated with temperature but not rainfall. B. dendrobatidis prevalence reached 35% in the cooler, drier months but was repeatedly undetectable during the warmer, wetter months. Also, B. dendrobatidis prevalence significantly decreased from 53.2% when the pathogen emerged on Montserrat in 2009 to a maximum 34.8% by 2011, after which it remained stable. On Dominica, where B. dendrobatidis emerged seven years prior to Montserrat, the same seasonal pattern was recorded but at lower prevalence, possibly indicating long-term decline. Understanding the dynamics of disease threats such as chytridiomycosis is key to planning conservation measures. For example, reintroductions of chytridiomycosis-threatened species could be timed to coincide with periods of low B. dendrobatidis infection risk, increasing potential for reintroduction success.

Fungal disease and temperature alter skin microbiome structure in an experimental salamander system

Pathogens compete with host microbiomes for space and resources. Their shared environment impacts pathogen-microbiome-host interactions, which can lead to variation in disease outcome. The skin microbiome of red-backed salamanders (Plethodon cinereus) can reduce infection by the pathogen Batrachochytrium dendrobatidis (Bd) at moderate infection loads, with high species richness and high abundance of competitors as putative mechanisms. However, it is unclear if the skin microbiome can reduce epizootic Bd loads across temperatures. We conducted a laboratory experiment to quantify skin microbiome and host responses (P. cinereus: n = 87) to Bd at mimicked epizootic loads across temperatures (13, 17 and 21°C). We quantified skin microbiomes using 16S rRNA gene metabarcoding and identified operational taxonomic units (OTUs) taxonomically similar to culturable bacteria known to kill Bd (anti-Bd OTUs). Prior to pathogen exposure, temperature changed the microbiome (OTU richness decreased by 12% and the abundance of anti-Bd OTUs increased by 18% per degree increase in temperature), but these changes were not predictive of disease outcome. After exposure, Bd changed the microbiome (OTU richness decreased by 0.1% and the abundance of anti-Bd OTUs increased by 0.2% per 1% increase in Bd load) and caused high host mortality across temperatures (35/45: 78%). Temperature indirectly impacted microbiome change and mortality through its direct effect on pathogen load. We did not find support for the microbiome impacting Bd load or host survival. Our research reveals complex host, pathogen, microbiome and environmental interactions to demonstrate that during epizootic events the microbiome will be unlikely to reduce pathogen invasion, even for putatively Bd-resistant species.

Ecological patterns in the skin microbiota of frogs from tropical Australia

The microbiota of frog skin can play an important role in protecting against diseases and parasites. The frog skin microbial community represents a complex mix of microbes that are promoted by the chemical environment of the frog skin and influenced by the animal's immediate past environment. The microbial communities of six species of frogs sampled from the campus of Charles Darwin University (CDU) were more similar within species than between species. The microbiota of the introduced cane toad (Rhinella marina) was most dissimilar among the species. Pairwise comparisons showed that the microbial communities of each species were different, except for the terrestrial Litoria nasuta and the arboreal L. rothii. The microbial communities of the six species were not related to ecological habit (arboreal or terrestrial), and neither was the alpha diversity of the microbes. The core microbes (defined as being on ≥90% of individuals of a species or group) were significantly different among all species, although 89 microbial operational taxonomic units (OTUs) were core microbes for all six species at CDU. Two species, Rhinella marina and Litoria rothii, were sampled at additional sites approximately 10 and 30 km from CDU. The microbial communities and the core OTU composition were different among the sites, but there were nevertheless 194 (R. marina) and 181 (L. rothii) core OTUs present at all three sites. Thus, the core microbiota varied with respect to geographic range and sample size.

Temperature-mediated inhibition of a bumblebee parasite by an intestinal symbiont

Competition between organisms is often mediated by environmental factors, including temperature. In animal intestines, nonpathogenic symbionts compete physically and chemically against pathogens, with consequences for host infection. We used metabolic theory-based models to characterize differential responses to temperature of a bacterial symbiont and a co-occurring trypanosomatid parasite of bumblebees, which regulate body temperature during flight and incubation. We hypothesized that inhibition of parasites by bacterial symbionts would increase with temperature, due to symbionts having higher optimal growth temperatures than parasites. We found that a temperature increase over the range measured in bumblebee colonies would favour symbionts over parasites. As predicted by our hypothesis, symbionts reduced the optimal growth temperature for parasites, both in direct competition and when parasites were exposed to symbiont spent medium. Inhibitory effects of the symbiont increased with temperature, reflecting accelerated growth and acid production by symbionts. Our results indicate that high temperatures, whether due to host endothermy or environmental factors, can enhance the inhibitory effects of symbionts on parasites. Temperature-modulated manipulation of microbiota could be one explanation for fever- and heat-induced reductions of infection in animals, with consequences for diseases of medical and conservation concern.

pH-mediated inhibition of a bumble bee parasite by an intestinal symbiont

Gut symbionts can augment resistance to pathogens by stimulating host-immune responses, competing for space and nutrients, or producing antimicrobial metabolites. Gut microbiota of social bees, which pollinate many crops and wildflowers, protect hosts against diverse infections and might counteract pathogen-related bee declines. Bumble bee gut microbiota, and specifically abundance of Lactobacillus 'Firm-5' bacteria, can enhance resistance to the trypanosomatid parasite Crithidia bombi. However, the mechanism underlying this effect remains unknown. We hypothesized that the Firm-5 bacterium Lactobacillus bombicola, which produces lactic acid, inhibits C. bombi via pH-mediated effects. Consistent with our hypothesis, L. bombicola spent medium inhibited C. bombi growth via reduction in pH that was both necessary and sufficient for inhibition. Inhibition of all parasite strains occurred within the pH range documented in honey bees, though sensitivity to acidity varied among strains. Spent medium was slightly more potent than HCl, d- and l-lactic acids for a given pH, suggesting that other metabolites also contribute to inhibition. Results implicate symbiont-mediated reduction in gut pH as a key determinant of trypanosomatid infection in bees. Future investigation into in vivo effects of gut microbiota on pH and infection intensity would test the relevance of these findings for bees threatened by trypanosomatids.

Disentangling causes of seasonal infection prevalence patterns: tropical tadpoles and chytridiomycosis as a model system

Identifying the factors that affect pathogen prevalence is critical to understanding the effects of wildlife diseases. We aimed to examine drivers of seasonal changes in the prevalence of infection by the amphibian chytrid fungus Batrachochytrium dendrobatidis in tadpoles. Because tadpoles may be important reservoirs for this disease, examining them will aid in understanding how chytridiomycosis affects entire amphibian populations. We hypothesized that temperature is a strong driver of prevalence of Bd in tadpoles, and the accumulation of infection as tadpoles become larger and older also drives prevalence in this system. We studied Litoria rheocola, a tropical rainforest stream frog with seasonal recruitment of annual tadpoles, and surveyed 6 streams in northeastern Queensland, Australia. Comparisons among models relating infection status to stream type, season, their interaction, tadpole age, and water temperature showed that age explained a large portion of the variance in infection status. Across sites and seasons, larger, older tadpoles had increased mean probabilities of infection, indicating that a large component of the variation among individuals was related to age, and thus to cumulative infection risk. Our results indicate that in systems with annual tadpoles, seasonal changes in infection prevalence may be strongly affected by seasonal patterns of tadpole growth and development in addition to stream type, season, and water temperature. These effects may then influence prevalence of infection in terrestrial individuals in species that have relatively frequent contact with water. This reinforces the need to integrate studies of the drivers of pathogen prevalence across all host life history stages.

Increased Numbers of Culturable Inhibitory Bacterial Taxa May Mitigate the Effects of Batrachochytrium dendrobatidis in Australian Wet Tropics Frogs

Symbiotic bacterial communities resident on amphibian skin can benefit their hosts. For example, antibiotic production by community members can control the pathogen Batrachochytrium dendrobatidis (Bd) and it is possible for these community members to be used as probiotics to reduce infection levels. In the early 1990s, the emergence of Bd caused declines and disappearances of frogs in the Australian Wet Tropics; the severity of its effects varied among species and sites. Some species have since recolonized despite enzootic Bd within their populations. This variation in history among species and sites provided an opportunity to investigate the role of anti-fungal cutaneous bacteria in protecting frogs against Bd infection. We collected cutaneous swab samples from three species of frogs at two upland and two lowland sites in the Wet Tropics, and used in vitro challenge assays to identify culturable Bd-inhibitory bacterial isolates for further analysis. We sequenced DNA from cultured inhibitory isolates to identify taxa, resulting in the classification of 16 Bd-inhibitory OTUs, and determined whether inhibitory taxa were associated with frog species, site, or intensity of infection. We present preliminary results showing that the upper limit of Bd infection intensity was negatively correlated with number of inhibitory OTUs present per frog indicating that increased numbers of Bd-inhibiting taxa may play a role in reducing the intensity of Bd infections, facilitating frog coexistence with enzootic Bd. One upland site had a significantly lower prevalence of Bd infection, a significantly higher proportion of frogs with one or more culturable Bd-inhibitory OTUs, a greater number of inhibitory bacterial genera present per frog, and statistically significant clustering of individual frogs with similar Bd-inhibitory signatures when compared to all other sites. This suggests that Bd-inhibitory taxa are likely to be particularly important to frogs at this site and may have played a role in their ability to recolonize following population declines. Our findings suggest that the use of multi-taxon Bd-inhibitory probiotics to support at-risk amphibian populations may be more effective than single-taxon alternatives.

Fragile coexistence of a global chytrid pathogen with amphibian populations is mediated by environment and demography

Unravelling the multiple interacting drivers of host-pathogen coexistence is crucial in understanding how an apparently stable state of endemism may shift towards an epidemic and lead to biodiversity loss. Here, we investigate the apparent coexistence of the global amphibian pathogen Batrachochytrium dendrobatidis (Bd) with Bombina variegata populations in The Netherlands over a 7-year period. We used a multi-season mark-recapture dataset and assessed potential drivers of coexistence (individual condition, environmental mediation and demographic compensation) at the individual and population levels. We show that even in a situation with a clear cost incurred by endemic Bd, population sizes remain largely stable. Current environmental conditions and an over-dispersed pathogen load probably stabilize disease dynamics, but as higher temperatures increase infection probability, changing environmental conditions, for example a climate-change-driven rise in temperature, could unbalance the current fragile host-pathogen equilibrium. Understanding the proximate mechanisms of such environmental mediation and of site-specific differences in infection dynamics can provide vital information for mitigation actions.

Temperature-Dependent Effects of Cutaneous Bacteria on a Frog's Tolerance of Fungal Infection

Defense against pathogens is one of many benefits that bacteria provide to animal hosts. A clearer understanding of how changes in the environment affect the interactions between animals and their microbial benefactors is needed in order to predict the impact and dynamics of emerging animal diseases. Due to its dramatic effects on the physiology of animals and their pathogens, temperature may be a key variable modulating the level of protection that beneficial bacteria provide to their animal hosts. Here we investigate how temperature and the makeup of the skin microbial community affect the susceptibility of amphibian hosts to infection by Batrachochytrium dendrobatidis (Bd), one of two fungal pathogens known to cause the disease chytridiomycosis. To do this, we manipulated the skin bacterial communities of susceptible hosts, northern cricket frogs (Acris crepitans), prior to exposing these animals to Bd under two different ecologically relevant temperatures. Our manipulations included one treatment where antibiotics were used to reduce the skin bacterial community, one where the bacterial community was augmented with the antifungal bacterium, Stenotrophomonas maltophilia, and one in which the frog’s skin bacterial community was left intact. We predicted that frogs with reduced skin bacterial communities would be more susceptible (i.e., less resistant to and/or tolerant of Bd infection), and frogs with skin bacterial communities augmented with the known antifungal bacterium would be less susceptible to Bd infection and chytridiomycosis. However, we also predicted that this interaction would be temperature dependent. We found a strong effect of temperature but not of skin microbial treatment on the probability and intensity of infection in Bd-exposed frogs. Whether temperature affected survival; however, it differed among our skin microbial treatment groups, with animals having more S. maltophilia on their skin surviving longer at 14 but not at 26°C. Our results suggest that temperature was the predominant factor influencing Bd’s ability to colonize the host (i.e., resistance) but that the composition of the cutaneous bacterial community was important in modulating the host’s ability to survive (i.e., tolerate) a heavy Bd infection.

Culture Media and Individual Hosts Affect the Recovery of Culturable Bacterial Diversity from Amphibian Skin

One current challenge in microbial ecology is elucidating the functional roles of the large diversity of free-living and host-associated bacteria identified by culture-independent molecular methods. Importantly, the characterization of this immense bacterial diversity will likely require merging data from culture-independent approaches with work on bacterial isolates in culture. Amphibian skin bacterial communities have become a recent focus of work in host-associated microbial systems due to the potential role of these skin bacteria in host defense against the pathogenic fungus Batrachochytrium dendrobatidis (Bd), which is associated with global amphibian population declines and extinctions. As there is evidence that some skin bacteria may inhibit growth of Bd and prevent infection in some cases, there is interest in using these bacteria as probiotic therapy for conservation of at-risk amphibians. In this study, we used skin swabs from American toads (Anaxyrus americanus) to: (1) assess the diversity and community structure of culturable amphibian skin bacteria grown on high and low nutrient culture media, (2) determine which culture media recover the highest proportion of the total skin bacterial community of individual toads relative to culture-independent data, and (3) assess whether the plated communities from the distinct media types vary in their ability to inhibit Bd growth in in-vitro assays. Overall, we found that culture media with low nutrient concentrations facilitated the growth of more diverse bacterial taxa and grew distinct communities relative to media with higher nutrient concentrations. Use of low nutrient media also resulted in culturing proportionally more of the bacterial diversity on individual toads relative to the overall community defined using culture-independent methods. However, while there were differences in diversity among media types, the variation among individual hosts was greater than variation among media types, suggesting that swabbing more individuals in a population is the best way to maximize culture collections, regardless of media type. Lastly, the function of the plated communities against Bd did not vary across culture media type or between high and low nutrient media. These results inform current efforts for developing a probiotic-based approach for amphibian conservation and help to ensure that culture collections are capturing the majority of the important diversity in these systems.

Inhibition of Fungal Pathogens across Genotypes and Temperatures by Amphibian Skin Bacteria

Symbiotic bacteria may dampen the impacts of infectious diseases on hosts by inhibiting pathogen growth. However, our understanding of the generality of pathogen inhibition by different bacterial taxa across pathogen genotypes and environmental conditions is limited. Bacterial inhibitory properties are of particular interest for the amphibian-killing fungal pathogens (Batrachochytrium dendrobatidis and Batrachochytrium salamandrivorans), for which probiotic applications as conservation strategies have been proposed. We quantified the inhibition strength of five putatively B. dendrobatidis-inhibitory bacteria isolated from woodland salamander skin against six Batrachochytrium genotypes at two temperatures (12 and 18°C). We selected six genotypes from across the Batrachochytrium phylogeny: B. salamandrivorans, B. dendrobatidis-Brazil and four genotypes of the B. dendrobatidis Global Panzootic Lineage (GPL1: JEL647, JEL404; GPL2: SRS810, JEL423). We performed 96-well plate challenge assays in a full factorial design. We detected a Batrachochytrium genotype by temperature interaction on bacterial inhibition score for all bacteria, indicating that bacteria vary in ability to inhibit Batrachochytrium depending on pathogen genotype and temperature. Acinetobacter rhizosphaerae moderately inhibited B. salamandrivorans at both temperatures (μ = 46–53%), but not any B. dendrobatidis genotypes. Chryseobacterium sp. inhibited three Batrachochytrium genotypes at both temperatures (μ = 5–71%). Pseudomonas sp. strain 1 inhibited all Batrachochytrium genotypes at 12°C and four Batrachochytrium genotypes at 18°C (μ = 5–100%). Pseudomonas sp. strain 2 and Stenotrophomonas sp. moderately to strongly inhibited all six Batrachochytrium genotypes at both temperatures (μ = 57–100%). All bacteria consistently inhibited B. salamandrivorans. Using cluster analysis of inhibition scores, we found that more closely related Batrachochytrium genotypes grouped together, suggesting that bacterial inhibition strength may be predictable based on Batrachochytrium relatedness. We conclude that bacterial inhibition capabilities change among bacterial strains, Batrachochytrium genotypes and temperatures. A comprehensive understanding of bacterial inhibitory function, across pathogen genotypes and temperatures, is needed to better predict the role of bacterial symbionts in amphibian disease ecology. For targeted conservation applications, we recommend using bacterial strains identified as strongly inhibitory as they are most likely to produce broad-spectrum antimicrobial agents at a range of temperatures.

Variation in Metabolite Profiles of Amphibian Skin Bacterial Communities Across Elevations in the Neotropics

Both the structure and function of host-associated microbial communities are potentially impacted by environmental conditions, just as the outcomes of many free-living species interactions are context-dependent. Many amphibian populations have declined around the globe due to the fungal skin pathogen, Batrachochytrium dendrobatidis (Bd), but enivronmental conditions may influence disease dynamics. For instance, in Panamá, the most severe Bd outbreaks have occurred at high elevation sites. Some amphibian species harbor bacterial skin communities that can inhibit the growth of Bd, and therefore, there is interest in understanding whether environmental context could also alter these host-associated microbial communities in a way that might ultimately impact Bd dynamics. In a field survey in Panamá, we assessed skin bacterial communities (16S rRNA amplicon sequencing) and metabolite profiles (HPLC-UV/Vis) of Silverstoneia flotator from three high- and three low-elevation populations representing a range of environmental conditions. Across elevations, frogs had similar skin bacterial communities, although one lowland site appeared to differ. Interestingly, we found that bacterial richness decreased from west to east, coinciding with the direction of Bd spread through Panamá. Moreover, metabolite profiles suggested potential functional variation among frog populations and between elevations. While the frogs have similar bacterial community structure, the local environment might shape the metabolite profiles. Ultimately, host-associated community structure and function could be dependent on environmental conditions, which could ultimately influence host disease susceptibility across sites.

Antifungal Bacteria on Woodland Salamander Skin Exhibit High Taxonomic Diversity and Geographic Variability

Diverse bacteria inhabit amphibian skin; some of those bacteria inhibit growth of the fungal pathogen Batrachochytrium dendrobatidis Yet there has been no systematic survey of anti-B. dendrobatidis bacteria across localities, species, and elevations. This is important given geographic and taxonomic variations in amphibian susceptibility to B. dendrobatidis Our collection sites were at locations within the Appalachian Mountains where previous sampling had indicated low B. dendrobatidis prevalence. We determined the numbers and identities of anti-B. dendrobatidis bacteria on 61 Plethodon salamanders (37 P. cinereus, 15 P. glutinosus, 9 P. cylindraceus) via culturing methods and 16S rRNA gene sequencing. We sampled co-occurring species at three localities and sampled P. cinereus along an elevational gradient (700 to 1,000 meters above sea level [masl]) at one locality. We identified 50 anti-B. dendrobatidis bacterial operational taxonomic units (OTUs) and found that the degree of B. dendrobatidis inhibition was not correlated with relatedness. Five anti-B. dendrobatidis bacterial strains occurred on multiple amphibian species at multiple localities, but none were shared among all species and localities. The prevalence of anti-B. dendrobatidis bacteria was higher at Shenandoah National Park (NP), VA, with 96% (25/26) of salamanders hosting at least one anti-B. dendrobatidis bacterial species compared to 50% (7/14) at Catoctin Mountain Park (MP), MD, and 38% (8/21) at Mt. Rogers National Recreation Area (NRA), VA. At the individual level, salamanders at Shenandoah NP had more anti-B. dendrobatidis bacteria per individual (μ = 3.3) than those at Catoctin MP (μ = 0.8) and at Mt. Rogers NRA (μ = 0.4). All salamanders tested negative for B. dendrobatidis Anti-B. dendrobatidis bacterial species are diverse in central Appalachian Plethodon salamanders, and their distribution varied geographically. The antifungal bacterial species that we identified may play a protective role for these salamanders.IMPORTANCE Amphibians harbor skin bacteria that can kill an amphibian fungal pathogen, Batrachochytrium dendrobatidis Some amphibians die from B. dendrobatidis infection, whereas others do not. The bacteria that can kill B. dendrobatidis, called anti-B. dendrobatidis bacteria, are thought to influence the B. dendrobatidis infection outcome for the amphibian. Yet how anti-B. dendrobatidis bacterial species vary among amphibian species and populations is unknown. We determined the distribution of anti-B. dendrobatidis bacterial species among three salamander species (n = 61) sampled at three localities. We identified 50 unique anti-B. dendrobatidis bacterial species and found that all of the tested salamanders were negative for B. dendrobatidis Five anti-B. dendrobatidis bacterial species were commonly detected, suggesting a stable, functional association with these salamanders. The number of anti-B. dendrobatidis bacteria per individual varied among localities but not among co-occurring salamander species, demonstrating that environment is more influential than host factors in structuring the anti-B. dendrobatidis bacterial community. These anti-B. dendrobatidis bacteria may serve a protective function for their salamander hosts.

Amphibian skin microbiota exhibits temporal variation in community structure but stability of predicted Bd-inhibitory function

Host-associated microbiomes are increasingly recognized to contribute to host disease resistance; the temporal dynamics of their community structure and function, however, are poorly understood. We investigated the cutaneous bacterial communities of three newt species, Ichthyosaura alpestris, Lissotriton vulgaris and Triturus cristatus, at approximately weekly intervals for 3 months using 16S ribosomal RNA amplicon sequencing. We hypothesized cutaneous microbiota would vary across time, and that such variation would be linked to changes in predicted fungal-inhibitory function. We observed significant temporal variation within the aquatic phase, and also between aquatic and terrestrial phase newts. By keeping T. cristatus in mesocosms, we demonstrated that structural changes occurred similarly across individuals, highlighting the non-stochastic nature of the bacterial community succession. Temporal changes were mainly associated with fluctuations in relative abundance rather than full turnover of bacterial operational taxonomic units (OTUs). Newt skin microbe fluctuations were not correlated with that of pond microbiota; however, a portion of community variation was explained by environmental temperature. Using a database of amphibian skin bacteria that inhibit the pathogen Batrachochytrium dendrobatidis (Bd), we found that the proportion of reads associated with 'potentially' Bd-inhibitory OTUs did not vary temporally for two of three newt species, suggesting that protective function may be maintained despite temporal variation in community structure.

Dietary Carotenoid Supplementation Enhances the Cutaneous Bacterial Communities of the Critically Endangered Southern Corroboree Frog (Pseudophryne corroboree)

The rapid spread of infectious disease has resulted in the decline of animal populations globally. Amphibians support a diversity of microbial symbionts on their skin surface that help to inhibit pathogen colonisation and reduce disease susceptibility and virulence. These cutaneous microbial communities represent an important component of amphibian immune defence, however, very little is known about the environmental factors that influence the cutaneous microbiome. Here, we characterise the cutaneous bacterial communities of a captive colony of the critically endangered Australian southern corroboree frog, Pseudophyrne corroboree, and examine the effect of dietary carotenoid supplementation on bacterial abundance, species richness and community composition. Individuals receiving a carotenoid-supplemented diet exhibited significantly higher bacterial abundance and species richness as well as an altered bacterial community composition compared to individuals that did not receive dietary carotenoids. Our findings suggest that dietary carotenoid supplementation enhances the cutaneous bacteria community of the southern corroboree frog and regulates the presence of bacteria species within the cutaneous microbiome. Our study is the second to demonstrate that carotenoid supplementation can improve amphibian cutaneous bacterial community dynamics, drawing attention to the possibility that dietary manipulation may assist with the ex situ management of endangered species and improve resilience to lethal pathogens such as Batrachochytrium dendrobatidis (Bd).

Environmental fluctuations and host skin bacteria shift survival advantage between frogs and their fungal pathogen

Fluctuating environments can modulate host-pathogen interactions by providing a temporary advantage to one of the interacting organisms. However, we know very little about how environmental conditions facilitate beneficial interactions between hosts and their microbial communities, resulting in individual persistence with a particular pathogen. Here, we experimentally infected Eleutherodactylus coqui frogs with the fungal pathogen Batrachochytrium dendrobatidis (Bd) under environmental conditions known to confer the survival advantage to the host during the warm-wet season, or alternatively to the pathogen during the cool-dry season. We used 16S rRNA amplicon sequencing to quantify changes in bacterial richness and phylogenetic diversity, and identified operational taxonomic units (OTUs) that became overrepresented or suppressed as a consequence of Bd infection. During the warm-wet season, frogs limited Bd infections, recruited putatively beneficial bacteria and returned to pre-infection levels of richness and phylogenetic diversity. In contrast, during the cool-dry season, Bd infections kept increasing through time, and bacterial diversity remained constant. Our findings confirm that infection outcome not only depends on abiotic factors, but also on biotic interactions between hosts and their associated bacterial communities.

Diversity of cutaneous microbiome of Pelophylax perezi populations inhabiting different environments

Potential synergetic interaction between chemicals, climate change and the emergence of opportunistic diseases is of utmost concern within the amphibian decline scenario. Understand the structure and dynamic of this microbiome and how environmental stressors act on this community is a priority. The present study aimed to: i) characterize the skin microbiome of Pelophylax perezi frog by looking for variations between populations in reference and under stress conditions (one metal contaminated and another with salinity fluctuations) and ii) evaluate the tolerance of skin-isolated bacteria to chemical contamination. Skin swabs were collected from frog populations inhabiting three reference sites (LB, VA and AM), one metal-rich contaminated (TP) and one brackish lentic aquatic system (SL). The frogs' skin microbiome was characterized by culture independent method (DGGE) and by cultivation methods. DGGE showed a characteristic profile in frogs from TP population. Results of recovered communities revealed low morphotype diversity and density (Colony Forming Units per frog) on individuals from TP population, comparatively to the other sampled populations. Isolated bacteria identified based on 16S rRNA gene sequence belong mainly to the classes Actinobacteria and Alphaproteobacteria. Ecotoxicological assays with acid-metal contaminated effluent (ETP) showed that the percentage of tolerant strains was higher in frogs from TP population than in SL and LB populations. In conclusion, results suggest that: a) environmental specific characteristics, as the presence of chemical contaminants, influence the composition of amphibian microbiome, comparing with sites without such contaminants; b) there are differences in microbiome composition between populations; and c) bacteria historically exposed to effluent may evolve tolerance to this kind of contamination.

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.

Using "Omics" and Integrated Multi-Omics Approaches to Guide Probiotic Selection to Mitigate Chytridiomycosis and Other Emerging Infectious Diseases

Emerging infectious diseases in wildlife are responsible for massive population declines. In amphibians, chytridiomycosis caused by Batrachochytrium dendrobatidis, Bd, has severely affected many amphibian populations and species around the world. One promising management strategy is probiotic bioaugmentation of antifungal bacteria on amphibian skin. In vivo experimental trials using bioaugmentation strategies have had mixed results, and therefore a more informed strategy is needed to select successful probiotic candidates. Metagenomic, transcriptomic, and metabolomic methods, colloquially called "omics," are approaches that can better inform probiotic selection and optimize selection protocols. The integration of multiple omic data using bioinformatic and statistical tools and in silico models that link bacterial community structure with bacterial defensive function can allow the identification of species involved in pathogen inhibition. We recommend using 16S rRNA gene amplicon sequencing and methods such as indicator species analysis, the Kolmogorov-Smirnov Measure, and co-occurrence networks to identify bacteria that are associated with pathogen resistance in field surveys and experimental trials. In addition to 16S amplicon sequencing, we recommend approaches that give insight into symbiont function such as shotgun metagenomics, metatranscriptomics, or metabolomics to maximize the probability of finding effective probiotic candidates, which can then be isolated in culture and tested in persistence and clinical trials. An effective mitigation strategy to ameliorate chytridiomycosis and other emerging infectious diseases is necessary; the advancement of omic methods and the integration of multiple omic data provide a promising avenue toward conservation of imperiled species.

Microbiome Variation Across Amphibian Skin Regions: Implications for Chytridiomycosis Mitigation Efforts

Cutaneous bacteria may play an important role in the resistance of amphibians to the pathogenic fungus Batrachochytrium dendrobatidis (Bd). Microbial communities resident on hosts' skin show topographical diversity mapping to skin features, as demonstrated by studies of the human microbiome. We examined skin microbiomes of wild and captive fire-bellied toads (Bombina orientalis) for differences across their body surface. We found that bacterial communities differed between ventral and dorsal skin. Wild toads showed slightly higher bacterial richness and diversity in the dorsal compared to the ventral region. On the other hand, captive toads hosted a higher richness and diversity of bacteria on their ventral than their dorsal skin. Microbial community composition and relative abundance of major bacterial taxonomic groups also differed between ventral and dorsal skin in all populations. Furthermore, microbiome diversity patterns varied as a function of their Bd infection status in wild toads. Bacterial richness and diversity was greater, and microbial community structure more complex, in wild than captive toads. The results suggest that bacterial community structure is influenced by microhabitats associated with skin regions. These local communities may be differentially modified when interacting with environmental bacteria and Bd. A better understanding of microbiome variation across skin regions will be needed to assess how the skin microbiota affects the abilities of amphibian hosts to resist Bd infection, especially in captive breeding programs.

Natural disturbance reduces disease risk in endangered rainforest frog populations

Natural disturbances can drive disease dynamics in animal populations by altering the microclimates experienced by hosts and their pathogens. Many pathogens are highly sensitive to temperature and moisture, and therefore small changes in habitat structure can alter the microclimate in ways that increase or decrease infection prevalence and intensity in host populations. Here we show that a reduction of rainforest canopy cover caused by a severe tropical cyclone decreased the risk of endangered rainforest frogs (Litoria rheocola) becoming infected by a fungal pathogen (Batrachochytrium dendrobatidis). Reductions in canopy cover increased the temperatures and rates of evaporative water loss in frog microhabitats, which reduced B. dendrobatidis infection risk in frogs by an average of 11–28% in cyclone-damaged areas, relative to unaffected areas. Natural disturbances to the rainforest canopy can therefore provide an immediate benefit to frogs by altering the microclimate in ways that reduce infection risk. This could increase host survival and reduce the probability of epidemic disease outbreaks. For amphibian populations under immediate threat from this pathogen, targeted manipulation of canopy cover could increase the availability of warmer, drier microclimates and therefore tip the balance from host extinction to coexistence.

Phylogenetic distribution of symbiotic bacteria from Panamanian amphibians that inhibit growth of the lethal fungal pathogen Batrachochytrium dendrobatidis

The introduction of next-generation sequencing has allowed for greater understanding of community composition of symbiotic microbial communities. However, determining the function of individual members of these microbial communities still largely relies on culture-based methods. Here, we present results on the phylogenetic distribution of a defensive functional trait of cultured symbiotic bacteria associated with amphibians. Amphibians are host to a diverse community of cutaneous bacteria and some of these bacteria protect their host from the lethal fungal pathogen Batrachochytrium dendrobatidis (Bd) by secreting antifungal metabolites. We cultured over 450 bacterial isolates from the skins of Panamanian amphibian species and tested their interactions with Bd using an in vitro challenge assay. For a subset of isolates, we also completed coculture experiments and found that culturing isolates with Bd had no effect on inhibitory properties of the bacteria, but it significantly decreased metabolite secretion. In challenge assays, approximately 75% of the bacterial isolates inhibited Bd to some extent and these inhibitory isolates were widely distributed among all bacterial phyla. Although there was no clear phylogenetic signal of inhibition, three genera, Stenotrophomonas, Aeromonas and Pseudomonas, had a high proportion of inhibitory isolates (100%, 77% and 73%, respectively). Overall, our results demonstrate that antifungal properties are phylogenetically widespread in symbiotic microbial communities of Panamanian amphibians and that some functional redundancy for fungal inhibition occurs in these communities. We hope that these findings contribute to the discovery and development of probiotics for amphibians that can mitigate the threat of chytridiomycosis.

Host-specific thermal profiles affect fitness of a widespread pathogen

Host behavior can interact with environmental context to influence outcomes of pathogen exposure and the impact of disease on species and populations. Determining whether the thermal behaviors of individual species influence susceptibility to disease can help enhance our ability to explain and predict how and when disease outbreaks are likely to occur. The widespread disease chytridiomycosis (caused by the fungal pathogen Batrachochytrium dendrobatidis, Bd) often has species-specific impacts on amphibian communities; some host species are asymptomatic, whereas others experience mass mortalities and population extirpation. We determined whether the average natural thermal regimes experienced by sympatric frog species in nature, in and of themselves, can account for differences in vulnerability to disease. We did this by growing Bd under temperatures mimicking those experienced by frogs in the wild. At low and high elevations, the rainforest frogs Litoria nannotis, L. rheocola, and L. serrata maintained mean thermal regimes within the optimal range for pathogen growth (15-25°C). Thermal regimes for L. serrata, which has recovered from Bd-related declines, resulted in slower pathogen growth than the cooler and less variable thermal regimes for the other two species, which have experienced more long-lasting declines. For L. rheocola and L. serrata, pathogen growth was faster in thermal regimes corresponding to high elevations than in those corresponding to low elevations, where temperatures were warmer. For L. nannotis, which prefers moist and thermally stable microenvironments, pathogen growth was fastest for low-elevation thermal regimes. All of the thermal regimes we tested resulted in pathogen growth rates equivalent to, or significantly faster than, rates expected from constant-temperature experiments. The effects of host body temperature on Bd can explain many of the broad ecological patterns of population declines in our focal species, via direct effects on pathogen fitness. Understanding the functional response of pathogens to conditions experienced by the host is important for determining the ecological drivers of disease outbreaks.