Symbiotic bacteria on amphibian skin can inhibit growth of the fungus Batrachochytrium dendrobatidis (Bd) that has caused dramatic population declines and extinctions of amphibians in the Neotropics. It remains unclear how the amphibians' skin microbiota is influenced by environmental bacterial reservoirs, host-associated factors such as susceptibility to pathogens, and pathogen presence in tropical amphibians. We sampled skin bacteria from five co-occurring frog species that differ in Bd susceptibility at one Bd-naive site, and sampled one of the non-susceptible species from Bd-endemic and Bd-naive sites in Panama. We hypothesized that skin bacterial communities (1) would be distinct from the surrounding environment regardless of the host habitat, (2) would differ between Bd susceptible and non-susceptible species and (3) would differ on hosts in Bd-naive and Bd-endemic sites. We found that skin bacterial communities were enriched in bacterial taxa that had low relative abundances in the environment. Non-susceptible species had very similar skin bacterial communities that were enriched in particular taxa such as the genera Pseudomonas and Acinetobacter. Bacterial communities of Craugastor fitzingeri in Bd-endemic sites were less diverse than in the naive site, and differences in community structure across sites were explained by changes in relative abundance of specific bacterial taxa. Our results indicate that skin microbial structure was associated with host susceptibility to Bd and might be associated to the history of Bd presence at different sites.

Vertebrates, including amphibians, host diverse symbiotic microbes that contribute to host disease resistance. Globally, and especially in montane tropical systems, many amphibian species are threatened by a chytrid fungus, Batrachochytrium dendrobatidis (Bd), that causes a lethal skin disease. Bd therefore may be a strong selective agent on the diversity and function of the microbial communities inhabiting amphibian skin. In Panamá, amphibian population declines and the spread of Bd have been tracked. In 2012, we completed a field survey in Panamá to examine frog skin microbiota in the context of Bd infection. We focused on three frog species and collected two skin swabs per frog from a total of 136 frogs across four sites that varied from west to east in the time since Bd arrival. One swab was used to assess bacterial community structure using 16S rRNA amplicon sequencing and to determine Bd infection status, and one was used to assess metabolite diversity, as the bacterial production of anti-fungal metabolites is an important disease resistance function. The skin microbiota of the three Panamanian frog species differed in OTU (operational taxonomic unit, ~bacterial species) community composition and metabolite profiles, although the pattern was less strong for the metabolites. Comparisons between frog skin bacterial communities from Panamá and the US suggest broad similarities at the phylum level, but key differences at lower taxonomic levels. In our field survey in Panamá, across all four sites, only 35 individuals (~26%) were Bd infected. There was no clustering of OTUs or metabolite profiles based on Bd infection status and no clear pattern of west-east changes in OTUs or metabolite profiles across the four sites. Overall, our field survey data suggest that different bacterial communities might be producing broadly similar sets of metabolites across frog hosts and sites. Community structure and function may not be as tightly coupled in these skin symbiont microbial systems as it is in many macro-systems.

Animal-associated microbiomes are integral to host health, yet key biotic and abiotic factors that shape host-associated microbial communities at the global scale remain poorly understood. We investigated global patterns in amphibian skin bacterial communities, incorporating samples from 2,349 individuals representing 205 amphibian species across a broad biogeographic range. We analyzed how biotic and abiotic factors correlate with skin microbial communities using multiple statistical approaches. Global amphibian skin bacterial richness was consistently correlated with temperature-associated factors. We found more diverse skin microbiomes in environments with colder winters and less stable thermal conditions, compared to environments with warm winters and less annual temperature variation. We used bioinformatically predicted bacterial growth rates, dormancy genes, and antibiotic synthesis genes, as well as inferred bacterial thermal growth optima to propose mechanistic hypotheses that may explain the observed patterns. We conclude that temporal and spatial characteristics of the host's macro-environment mediate microbial diversity.

Chytridiomycosis is a globally emerging disease of amphibians and the leading cause of population declines and extirpations at species-diverse montane sites in Central America. We continued long-term monitoring efforts for the presence of the fungal pathogen Batrachochytrium dendrobatidis (Bd) and for amphibian populations at two sites in western Panama, and we began monitoring at three new sites to the east.Population declines associated with chytridiomycosis emergence were detected at Altos de Campana National Park. We also detected Bd in three species east of the Panama Canal at Soberania National Park, and prevalence data suggests that Bd may be enzootic in the lowlands of the park. However, no infected frogs were found further east at Torti (prevalence <7.5% with 95% confidence). Our results suggest that Panama's diverse and not fully described amphibian communities east of the canal are at risk. Precise predictions of future disease emergence events are not possible until factors underlying disease emergence, such as dispersal, are understood.

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.

Host-associated microbial communities are ubiquitous among animals, and serve important functions. For example, the bacterial skin microbiome of amphibians can play a role in preventing or reducing infection by the amphibian chytrid fungus, Batrachochytrium dendrobatidis. Evidence suggests that environmental bacteria likely serve as a source pool for at least some of the members of the amphibian skin bacterial community, underscoring the potential for local environmental changes to disrupt microbial community source pools that could be critical to the health of host organisms. However, few studies have assessed variation in the amphibian skin microbiome along clear environmental gradients, and so we know relatively little about how local environmental conditions influence microbiome diversity. We sampled the skin bacterial communities of Coqui frogs, Eleutherodactylus coqui (N = 77), along an elevational gradient in eastern Puerto Rico (0–875 m), with transects in two land use types: intact forest (N = 4 sites) and disturbed (N = 3 sites) forest. We found that alpha diversity (as assessed by Shannon, Simpson, and Phylogenetic Diversity indices) varied across sites, but this variation was not correlated with elevation or land use. Beta diversity (community structure), on the other hand, varied with site, elevation and land use, primarily due to changes in the relative abundance of certain bacterial OTUs (∼species) within these communities. Importantly, although microbiome diversity varied, E. coqui maintained a common core microbiota across all sites. Thus, our findings suggest that environmental conditions can influence the composition of the skin microbiome of terrestrial amphibians, but that some aspects of the microbiome remain consistent despite environmental variation.

The fungal disease chytridiomycosis, caused by Batrachochytrium dendrobatidis (Bd), is one of the main causes of amphibian population declines and extinctions all over the world. In the Neotropics, this fungal disease has caused catastrophic declines in the highlands as it has spread throughout Central America down to Panamá. In this study, we determined the prevalence and intensity of Bd infection in three species of frogs in one highland and four lowland tropical forests, including two lowland regions in eastern Panamá in which the pathogen had not been detected previously. Bd was present in all the sites sampled with a prevalence ranging from 15–34%, similar to other Neotropical lowland sites. The intensity of Bd infection on individual frogs was low, ranging from average values of 0.11–24 zoospore equivalents per site. Our work indicates that Bd is present in anuran communities in lowland Panamá, including the Darién province, and that the intensity of the infection may vary among species from different habitats and with different life histories. The population-level consequences of Bd infection in amphibian communities from the lowlands remain to be determined. Detailed studies of amphibian species from the lowlands will be essential to determine the reason why these species are persisting despite the presence of the pathogen.

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.