Temporal Variation of the Skin Bacterial Community and Batrachochytrium dendrobatidis Infection in the Terrestrial Cryptic Frog Philoria loveridgei

Gene functions of the Ambystoma altamirani skin microbiome vary across space and time but potential antifungal genes are widespread and prevalent

Amphibian skin microbiomes can play a critical role in host survival against emerging diseases by protecting their host against pathogens. While a plethora of biotic and abiotic factors have been shown to influence the taxonomic diversity of amphibian skin microbiomes it remains unclear whether functional genomic diversity varies in response to temporal and environmental factors. Here we applied a metagenomic approach to evaluate whether seasonality, distinct elevations/sites, and pathogen presence influenced the functional genomic diversity of the A. altamirani skin microbiome. We obtained a gene catalogue of 92 107 nonredundant annotated genes and a set of 50 unique metagenome assembled genomes (MAGs). Our analysis showed that genes linked to general and potential antifungal traits significantly differed across seasons and sampling locations at different elevations. Moreover, we found that the functional genomic diversity of A. altamirani skin microbiome differed between B. dendrobatidis infected and not infected axolotls only during winter, suggesting an interaction between seasonality and pathogen infection. In addition, we identified the presence of genes and biosynthetic gene clusters (BGCs) linked to potential antifungal functions such as biofilm formation, quorum sensing, secretion systems, secondary metabolite biosynthesis, and chitin degradation. Interestingly genes linked to these potential antifungal traits were mainly identified in Burkholderiales and Chitinophagales MAGs. Overall, our results identified functional traits linked to potential antifungal functions in the A. altamirani skin microbiome regardless of variation in the functional diversity across seasons, elevations/sites, and pathogen presence. Our findings suggest that potential antifungal traits found in Burkholderiales and Chitinophagales taxa could be related to the capacity of A. altamirani to survive in the presence of Bd, although further experimental analyses are required to test this hypothesis.

Low-level pathogen infection and geographic location correlate with the skin microbiomes of Columbia spotted frogs (Rana luteiventris) in a montane landscape

The skin microbiome of amphibians can influence host susceptibility towards the fungal pathogen Batrachochytrium dendrobatidis (Bd), while simultaneously having the potential to be altered by Bd. Severe Bd infections are known to alter the amphibian skin microbiome; however, little is known about microbiome interactions in amphibians with low infection intensity. In addition to disease dynamics, environmental factors may influence the microbiome. To test for patterns in bacterial diversity based on pathogen infection and environmental factors, 399 Columbia spotted frogs (Rana luteiventris) were sampled throughout northern Idaho and northeastern Washington across two years. Bd prevalence and intensity were measured in 376 frogs, revealing a prevalence of 69%, but generally low infection intensity (Mean = 127 Bd zoospore equivalents among infected frogs). Skin bacterial communities were characterized in 92 frogs using 16S rRNA gene amplicon sequencing. Our results indicated correlations of decreasing Shannon diversity and evenness as infection intensity increased. Latitude was correlated with bacterial richness and Faith's Phylogenetic Diversity measures, indicating increased diversity in northern locations. Beta diversity (UniFrac) analyses revealed that skin microbiomes were distinct between infected and uninfected frogs, and infection intensity had a significant effect on microbiome composition. Site explained the majority of microbiome variation (weighted UniFrac: 57.5%), suggesting a combination of local habitat conditions explain variation, as only small proportions of variation could be explained by year, month, temperature, elevation, and latitude individually. Bacterial genera with potential for Bd-inhibitory properties were found with differential relative abundance in infected and uninfected frogs, with higher Stenotrophomonas and lower Pseudomonas relative abundance observed in infected frogs. Further study may indicate if Bd inhibition by members of the skin microbiome is an influence behind the low infection intensities observed and whether low Bd infection intensities are capable of altering skin microbiome composition.

Modeling the Distribution of the Chytrid Fungus Batrachochytrium dendrobatidis with Special Reference to Ukraine

Amphibians are the most threatened group of vertebrates. While habitat loss poses the greatest threat to amphibians, a spreading fungal disease caused by Batrachochytrium dendrobatidis Longcore, Pessier & D.K. Nichols 1999 (Bd) is seriously affecting an increasing number of species. Although Bd is widely prevalent, there are identifiable heterogeneities in the pathogen's distribution that are linked to environmental parameters. Our objective was to identify conditions that affect the geographic distribution of this pathogen using species distribution models (SDMs) with a special focus on Eastern Europe. SDMs can help identify hotspots for future outbreaks of Bd but perhaps more importantly identify locations that may be environmental refuges ("coldspots") from infection. In general, climate is considered a major factor driving amphibian disease dynamics, but temperature in particular has received increased attention. Here, 42 environmental raster layers containing data on climate, soil, and human impact were used. The mean annual temperature range (or 'continentality') was found to have the strongest constraint on the geographic distribution of this pathogen. The modeling allowed to distinguish presumable locations that may be environmental refuges from infection and set up a framework to guide future search (sampling) of chytridiomycosis in Eastern Europe.

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.

The skin microbiota of the axolotl Ambystoma altamirani is highly influenced by metamorphosis and seasonality but not by pathogen infection

BACKGROUND

Microbiomes have been increasingly recognized as major contributors to host health and survival. In amphibians, bacterial members of the skin microbiota protect their hosts by inhibiting the growth of the fungal pathogen Batrachochytrium dendrobatidis (Bd). Even though several studies describe the influence of biotic and abiotic factors over the skin microbiota, it remains unclear how these symbiotic bacterial communities vary across time and development. This is particularly relevant for species that undergo metamorphosis as it has been shown that host physiology and ecology drastically influence diversity of the skin microbiome.

RESULTS

We found that the skin bacterial communities of the axolotl A. altamirani are largely influenced by the metamorphic status of the host and by seasonal variation of abiotic factors such as temperature, pH, dissolved oxygen and conductivity. Despite high Bd prevalence in these samples, the bacterial diversity of the skin microbiota did not differ between infected and non-infected axolotls, although relative abundance of particular bacteria were correlated with Bd infection intensity.

CONCLUSIONS

Our work shows that metamorphosis is a crucial process that shapes skin bacterial communities and that axolotls under different developmental stages respond differently to environmental seasonal variations. Moreover, this study greatly contributes to a better understanding of the factors that shape amphibian skin microbiota, especially in a largely underexplored group like axolotls (Mexican Ambystoma species).

Temporal variation in skin microbiota of cohabitating amphibians

Temporal changes and transmission patterns in host-associated microbial communities have important implications for host health. The diversity of amphibian skin microbial communities is associated with disease outcome in amphibians exposed to the fungal pathogen Batrachochytrium dendrobatidis (Bd). To successfully develop conservation strategies against Bd, we need a comprehensive understanding of how skin microbes are maintained and transmitted over time within populations. We used 16S rRNA sequence analysis to compare Epipedobates anthonyi frogs housed with one conspecific to frogs housed singly at four time points over the course of 1 year. We found that both α and β diversity of frog skin bacterial communities changed significantly over the course of the experiment. Specifically, we found that bacterial communities of cohabitating frogs became more similar over time. We also observed that some bacterial taxa were differentially abundant between frogs housed singly and frogs housed with a conspecific. These results suggest that conspecific contact may play a role in mediating amphibian skin microbial diversity and that turnover of skin microbial communities can occur across time. Our findings provide rationale for future studies exploring horizontal transmission as a potential mechanism of host-associated microbial maintenance in amphibians.

Experimental test of microbiome protection across pathogen doses reveals importance of resident microbiome composition

The commensal microbes inhabiting a host tissue can interact with invading pathogens and host physiology in ways that alter pathogen growth and disease manifestation. Prior work in house finches (Haemorhous mexicanus) found that resident ocular microbiomes were protective against conjunctival infection and disease caused by a relatively high dose of Mycoplasma gallisepticum. Here, we used wild-caught house finches to experimentally examine whether protective effects of the resident ocular microbiome vary with the dose of invading pathogen. We hypothesized that commensal protection would be strongest at low M. gallisepticum inoculation doses because the resident microbiome would be less disrupted by invading pathogen. Our five M. gallisepticum dose treatments were fully factorial with an antibiotic treatment to perturb resident microbes just prior to M. gallisepticum inoculation. Unexpectedly, we found no indication of protective effects of the resident microbiome at any pathogen inoculation dose, which was inconsistent with the prior work. The ocular bacterial communities at the beginning of our experiment differed significantly from those previously reported in local wild-caught house finches, likely causing this discrepancy. These variable results underscore that microbiome-based protection in natural systems can be context dependent, and natural variation in community composition may alter the function of resident microbiomes in free-living animals.

Bacterial community dynamics on bats and the implications for pathogen resistance

The bats skin microbiota plays an important role in reducing pathogen infection, including the deadly fungal pathogen Pseudogymnoascus destructans, the causative agent of white-nose syndrome. However, the dynamic of skin bacterial communities response to environmental perturbations remains poorly described. We characterized skin bacterial community over time and space in Rhinolophus ferrumequinum, a species with high resistance to the infection with P. destructans. We collected environmental covariate data to determine what factors influenced changes in community structure. We observed significant temporal and spatial shifts in the skin bacterial community, which was mainly associated with variation in operational taxonomic units. The skin bacterial community differed by the environmental microbial reservoirs and was most influenced by host body condition, bat roosting temperature and geographic distance between sites, but was not influenced by pathogen infection. Furthermore, the skin microbiota was enriched in particular taxa with antifungal abilities, such as Enterococcus, Burkholderia, Flavobacterium, Pseudomonas, Corynebacterium and Rhodococcus. And specific strains of Pseudomonas, Corynebacterium and Rhodococcus even inhibited P. destructans growth. Our findings provide new insights in characterizing the variation in bacterial communities can inform us about the processes of driving community assembly and predict the host's ability to resist or survive pathogen infection.

Longitudinal patterns in the skin microbiome of wild, individually marked frogs from the Sierra Nevada, California

The amphibian skin microbiome has been the focus of numerous studies because of the protective effects that some bacteria provide against the pathogenic fungus Batrachochytrium dendrobatidis, which has caused a global panzootic among amphibians. However, the mechanisms driving community structure and function in the amphibian skin microbiome are still poorly understood, and longitudinal analyses of the skin microbiome have not yet been conducted in wild populations. In this study, we investigate longitudinal patterns in the skin microbiome of 19 individually marked adult frogs from two wild populations of the endangered Sierra Nevada yellow-legged frog (Rana sierrae), sampled over the course of 2 years. We found that individuals with low bacterial diversity (dominated by order Burkhorderiales) had significantly more stable bacterial communities than those with higher diversity. Amplicon sequence variants (ASVs) with high relative abundance were significantly less transient than those with low relative abundance, and ASVs with intermediate-level relative abundances experienced the greatest volatility over time. Based on these results, we suggest that efforts to develop probiotic treatments to combat B. dendrobatidis should focus on bacteria that are found at high relative abundances in some members of a population, as these strains are more likely to persist and remain stable in the long term.

Assembly of the amphibian microbiome is influenced by the effects of land-use change on environmental reservoirs

A growing focus in microbial ecology is understanding how beneficial microbiome function is created and maintained through various assembly mechanisms. This study explores the role of both the environment and disease in regulating the composition of microbial species in the soil and on amphibian hosts. We compared the microbial communities of Plethodon cinereus salamanders along a land-use gradient in the New York metropolitan area and paired these with associated soil cores. Additionally, we characterized the diversity of bacterial and fungal symbionts that putatively inhibit the pathogenic fungus Batrachochytrium dendrobatidis. We predicted that variation in skin microbial community composition would correlate with changes seen in the soil which functions as the regional species pool. We found that salamanders and soil share many microbial taxa but that these two communities exhibit differences in the relative abundances of the bacterial phyla Acidobacteria, Actinobacteria, and Proteobacteria and the fungal phyla Ascomycota and genus Basidiobolus. Microbial community composition varies with changes in land-use associated factors creating site-specific compositions. By employing a quantitative, null-based assembly model, we identified that dispersal limitation, variable selection, and drift guide assembly of microbes onto their skin, creating high dissimilarity between individuals with likely consequences in disease preventative function.

Seasonal changes and the unexpected impact of environmental disturbance on skin bacteria of individual amphibians in a natural habitat

Amphibians host diverse skin bacteria that have a role in pathogen defense, but these skin communities could change over time and impact this function. Here, we monitored individual Eastern red-spotted newts (Notophthalmus viridescens; N = 17) for 2 years in a field pond enclosure and assessed the effects of season and disturbance on skin bacterial community dynamics. We created disturbances by adding additional pond substrate to the enclosure at two timepoints. We planned to sample the skin bacterial community and metabolite profiles of each newt every 6 weeks; we ultimately sampled eight individuals at least six times. We used 16S rRNA gene amplicon sequencing to characterize the bacterial communities and HPLC-MS for metabolite profiling. We found that disturbance had a dramatic effect on skin bacterial communities and metabolite profiles, while season had an effect only using select metrics. There were seven core bacterial taxa (97% OTUs) that were found on all newts in all seasons, pre- and post-disturbance. Lastly, there was a correlation between bacterial and metabolite profiles post-disturbance, which was not observed pre-disturbance. This longitudinal study suggests that environmental disturbances can have lasting effects on skin bacterial communities that overwhelm seasonal changes, although the core bacteria remain relatively consistent over 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.

Relationships among host microbiota, parasite resistance or tolerance, and host fitness

Article impact statement: Parasite-induced shifts in host microbiota that lead to parasite resistance or tolerance may have unintended consequences for host fitness.

Altitudinal Effects on Innate Immune Response of a Subterranean Rodent

Immune defense is costly to maintain and deploy, and the optimal investment into immune defense depends on risk of infection. Altitude is a natural environmental factor that is predicted to affect parasite abundance, with lower parasite abundance predicted at higher altitudes due to stronger environmental stressors, which reduce parasite transmission. Using high and low altitude populations of the Turkish blind mole-rat (TBMR) Nannospalax xanthodon, we tested for effects of altitude on constitutive innate immune defense. Field studies were performed with 32 wild animals in 2017 and 2018 from two low- and one high-altitude localities in the Central Taurus Mountains, at respective altitudes of 1010 m, 1115 m, and 2900 m above sea level. We first compared innate standing immune defense as measured by the bacteria-killing ability of blood serum. We then measured corticosterone stress hormone levels, as stressful conditions may affect immune response. Finally, we compared prevalence and intensity of gastrointestinal parasites of field-captured TBMR. We found that the bacteria-killing ability of serum is greater in the mole-rat samples from high altitude. There was no significant difference in stress (corticosterone) levels between altitude categories. Coccidian prevalence and abundance were significantly higher in 2017 than 2018 samples, but there was no significant difference in prevalence, abundance, or intensity between altitudes, or between sexes. Small sample sizes may have reduced power to detect true differences; nevertheless, this study provides support that greater standing innate immunity in high altitude animals may reflect greater investment into constitutive defense.

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.

Development and Comparison of Complementary Methods to Study Potential Skin and Inhalational Exposure to Pathogens During Personal Protective Equipment Doffing

BACKGROUND

Fluorescent tracers are often used with ultraviolet lights to visibly identify healthcare worker self-contamination after doffing of personal protective equipment (PPE). This method has drawbacks, as it cannot detect pathogen-sized contaminants nor airborne contamination in subjects' breathing zones.

METHODS

A contamination detection/quantification method was developed using 2-µm polystyrene latex spheres (PSLs) to investigate skin contamination (via swabbing) and potential inhalational exposure (via breathing zone air sampler). Porcine skin coupons were used to estimate the PSL swabbing recovery efficiency and limit of detection (LOD). A pilot study with 5 participants compared skin contamination levels detected via the PSL vs fluorescent tracer methods, while the air sampler quantified potential inhalational exposure to PSLs during doffing.

RESULTS

Average PSL skin swab recovery efficiency was 40% ± 29% (LOD = 1 PSL/4 cm2 of skin). In the pilot study, all subjects had PSL and fluorescent tracer skin contamination. Two subjects had simultaneously located contamination of both types on a wrist and hand. However, for all other subjects, the PSL method enabled detection of skin contamination that was not detectable by the fluorescent tracer method. Hands/wrists were more commonly contaminated than areas of the head/face (57% vs 23% of swabs with PSL detection, respectively). One subject had PSLs detected by the breathing zone air sampler.

CONCLUSIONS

This study provides a well-characterized method that can be used to quantitate levels of skin and inhalational contact with simulant pathogen particles. The PSL method serves as a complement to the fluorescent tracer method to study PPE doffing self-contamination.

Moving Beyond the Host: Unraveling the Skin Microbiome of Endangered Costa Rican Amphibians

Some neotropical amphibians, including a few species in Costa Rica, were presumed to be "extinct" after dramatic population declines in the late 1980s but have been rediscovered in isolated populations. Such populations seem to have evolved a resistance/tolerance to Batrachochytrium dendrobatidis (Bd), a fungal pathogen that causes a deadly skin disease and is considered one of the main drivers of worldwide amphibian declines. The skin microbiome is an important component of the host's innate immune system and is associated with Bd-resistance. However, the way that the bacterial diversity of the skin microbiome confers protection against Bd in surviving species remains unclear. We studied variation in the skin microbiome and the prevalence of putatively anti-Bd bacterial taxa in four co-habiting species in the highlands of the Juan Castro Blanco National Park in Costa Rica using 16S rRNA amplicon sequencing. Lithobates vibicarius, Craugastor escoces, and Isthmohyla rivularis have recently been rediscovered, whereas Isthmohyla pseudopuma has suffered population fluctuations but has never disappeared. To investigate the life stage at which the protective skin microbiome is shaped and when shifts occur in the diversity of putatively anti-Bd bacteria, we studied the skin microbiome of tadpoles, juveniles and adults of L. vibicarius. We show that the skin bacterial composition of sympatric species and hosts with distinct Bd-infection statuses differs at the phyla, family, and genus level. We detected 94 amplicon sequence variants (ASVs) with putative anti-Bd activity pertaining to distinct bacterial taxa, e.g., Pseudomonas spp., Acinetobacter johnsonii, and Stenotrophomonas maltophilia. Bd-uninfected L. vibicarius harbored 79% more putatively anti-Bd ASVs than Bd-infected individuals. Although microbiome composition and structure differed across life stages, the diversity of putative anti-Bd bacteria was similar between pre- and post-metamorphic stages of L. vibicarius. Despite low sample size, our results support the idea that the skin microbiome is dynamic and protects against ongoing Bd presence in endangered species persisting after their presumed extinction. Our study serves as a baseline to understand the microbial patterns in species of high conservation value. Identification of microbial signatures linked to variation in disease susceptibility might, therefore, inform mitigation strategies for combating the global decline of amphibians.

Skin bacterial communities of neotropical treefrogs vary with local environmental conditions at the time of sampling

The amphibian skin microbiome has been the focus of recent studies aiming to better understand the role of these microbial symbionts in host defense against disease. However, host-associated microbial communities are complex and dynamic, and changes in their composition and structure can influence their function. Understanding temporal variation of bacterial communities on amphibian skin is critical for establishing baselines from which to improve the development of mitigation techniques based on probiotic therapy and provides long-term host protection in a changing environment. Here, we investigated whether microbial communities on amphibian skin change over time at a single site. To examine this, we collected skin swabs from two pond-breeding species of treefrogs, Agalychnis callidryas and Dendropsophus ebraccatus, over 4 years at a single lowland tropical pond in Panamá. Relative abundance of operational taxonomic units (OTUs) based on 16S rRNA gene amplicon sequencing was used to determine bacterial community diversity on the skin of both treefrog species. We found significant variation in bacterial community structure across long and short-term time scales. Skin bacterial communities differed across years on both species and between seasons and sampling days only in D. ebraccatus. Importantly, bacterial community structures across days were as variable as year level comparisons. The differences in bacterial community were driven primarily by differences in relative abundance of key OTUs and explained by rainfall at the time of sampling. These findings suggest that skin-associated microbiomes are highly variable across time, and that for tropical lowland sites, rainfall is a good predictor of variability. However, more research is necessary to elucidate the significance of temporal variation in bacterial skin communities and their maintenance for amphibian conservation efforts.

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.