Co-habiting amphibian species harbor unique skin bacterial communities in wild populations.

Although all plant and animal species harbor microbial symbionts, we know surprisingly little about the specificity of microbial communities to their hosts. Few studies have compared the microbiomes of different species of animals, and fewer still have examined animals in the wild. We sampled four pond habitats in Colorado, USA, where multiple amphibian species were present. In total, 32 amphibian individuals were sampled from three different species including northern leopard frogs (Lithobates pipiens), western chorus frogs (Pseudacris triseriata) and tiger salamanders (Ambystoma tigrinum). We compared the diversity and composition of the bacterial communities on the skin of the collected individuals via barcoded pyrosequencing of the 16S rRNA gene. Dominant bacterial phyla included Acidobacteria, Actinobacteria, Bacteriodetes, Cyanobacteria, Firmicutes and Proteobacteria. In total, we found members of 18 bacterial phyla, comparable to the taxonomic diversity typically found on human skin. Levels of bacterial diversity varied strongly across species: L. pipiens had the highest diversity; A. tigrinum the lowest. Host species was a highly significant predictor of bacterial community similarity, and co-habitation within the same pond was not significant, highlighting that the skin-associated bacterial communities do not simply reflect those bacterial communities found in their surrounding environments. Innate species differences thus appear to regulate the structure of skin bacterial communities on amphibians. In light of recent discoveries that some bacteria on amphibian skin have antifungal activity, our finding suggests that host-specific bacteria may have a role in the species-specific resistance to fungal pathogens.

Introduction

All species of plants and animals harbor assemblages of microbes and these symbionts often have vital roles in the health of their hosts. However, there have been few comprehensive studies assessing whether different species harbor distinct microbial symbionts. Work on mammals (Ley ), primates (Yildirim ) and plants (Redford ) suggests that there is some degree of host specificity of microbial symbionts across taxa. However, such studies have not yet been performed in wild populations, nor have these studies been designed to distinguish environmental factors from innate host factors. For different species in the wild, the question of host-specific microbial assemblages becomes difficult to address because species often occupy different habitat niches, and different environmental conditions may confound questions of host specificity. For several reasons, we have identified pond-dwelling larval amphibian communities as a model system to address whether different, co-habiting species harbor unique microbial communities. First, populations of larval amphibians observed before metamorphosis (for example, tadpoles), develop within a single pond and therefore have a known environmental origin. Second, ponds can harbor several species of co-habiting larval amphibians that occupy the same environmental niche during their development. Third, different ponds that harbor multiple species of larval amphibians can serve as replicates to examine patterns of host-associated microbial community assembly across different sites.

We focused specifically on skin-associated bacterial communities of amphibians as these communities likely provide the first line of defense against pathogens and, for co-habiting species, are exposed to the same environmental sources of microbes. One specific skin pathogen of global significance is the fungus Batrachochytrium dendrobatidis (Bd), which has a key role in amphibian declines worldwide (Skerratt ; Fisher ). Amphibian species, and some populations within species, vary in their susceptibility to disease caused by this fungal pathogen (Kilpatrick ). While there are many factors that may have a role in amphibian disease resistance including host genetics, immunology (Carey ), skin peptides (Rollins-Smith 2009) and environmental factors (Lips ), microbial communities living on the skin of amphibians may also have a crucial role in host defense mechanisms. Certain bacteria isolated from amphibian skin are known to produce fungicides and can inhibit growth of Bd (Harris ) and even reduce amphibian mortality under controlled experimental conditions (Harris ). However, it is not known if different species of wild amphibians can select unique skin microbial assemblages from exposure to the same microbial background. If specific microbial assemblages are involved in resistance to pathogens (such as Bd), and amphibian species vary in their susceptibility to disease, then we would expect that unique skin microbial composition across species or populations might provide a mechanism underlying this observed pattern of differences in susceptibility. Our goal in the present study was to use co-habiting, pond-dwelling larval amphibians to determine whether different species of amphibians harbor unique communities of bacteria on their skin. If they do, it suggests that species-specific bacterial communities have the potential to contribute to differences in the resistance of individual species to skin pathogens.

Teasing apart how symbiotic microbial communities are influenced by factors such as the environment, diet or innate host traits can now be addressed with the use of newer high-throughput survey techniques such as barcoded pyrosequencing of the 16S rRNA gene. These methods allow for observing the composition and relative abundance of microbes with greater accuracy than previous methods used to estimate bacterial community structure (Costello ). By applying these newer techniques in pond-dwelling amphibians, we hope to gain insight on the host-specific nature of skin microbes independent of environmental confounds as well as build a foundation to understand the potential for these microbes to have a role in host resistance against pathogens.

Materials and methods

Amphibian species and field collections

Four pond habitats in Colorado were identified in 2008 that each harbored at least two amphibian species that are common in Colorado and co-habit in ponds during the summer breeding season. Amphibian species sampled included Lithobates (=Rana) pipiens (northern leopard frog), Pseudacris triseriata (western chorus frog) and Ambystoma tigrinum (tiger salamander). Each pond was less than one hectare in surface area and located in a separate watershed. Site characteristics and the number of individual amphibians sampled per species per pond are described in Table 1. Amphibians were sampled for skin bacteria during August 2008 and all individuals from a given pond were sampled on the same day. Pre-metamorphic, larval amphibians were targeted because they develop over one season within a single pond, therefore, we can ensure that our cross-species comparisons only include individuals that shared the same environment and microbial inocula during their development. All three species named above are pond-breeders that lay eggs in early summer, which hatch larval swimming stages that develop over the course of the summer into metamorphic stages (Hammerson, 1999). This system allows us to control for confounding environmental conditions and isolate the effect of host species.

Table 1

Amphibian sampling scheme

Pond name	# Leopard frogs sampled	# Chorus frogs sampled	# Tiger salamanders sampled	Pond elevation (m)	Pond location in Colorado	Latitude longitude (UTM)
Monster bull	4	2	4	2357	Meeker	40.010359–107.625487
Chive	0	4	4	2598	Nederland	39.951877–105.463940
Horse	0	4	4	1702	Southwest Boulder	39.946662–105.241441
Eggleston3	3	4	0	1754	Southeast Boulder	39.927857–105.20637

All three species co-occurred at one site and other sites had two species co-occurring.

Larval amphibians were captured from the water using hand-held dipnets or a blocking net seine. Within a given species, larval individuals of a similar size were selected for skin microbe sampling. Larval amphibians were held temporarily (<10 min) in plastic containers with pond water before sampling. In total, 32 individual amphibians were sampled (7L. pipiens, 14 P. triseriata and 12 A. tigrinum; see Table 1) from four ponds. Fresh latex gloves were used for every individual amphibian handled. Before sampling, each individual was held and rinsed with 100 ml of sterile water three times to ensure that the skin sample primarily included skin-associated microbes rather than pond-associated material, including pond water, sediment and transient microbes (Culp ; Lauer ). Prior studies by Culp and Lauer demonstrate that the composition of bacteria obtained from amphibian skin versus rinse water mostly differ and represent unique assemblages, thus, we have some precedent to assume that most of the bacteria observed in this study are likely associated with the amphibians rather than transient bacteria from the environment. Immediately following rinsing, each amphibian was sampled using a sterile cotton-tipped swab brushed over the entire body of the amphibian for 30 s. The small size of larval amphibians (3–12 cm, depending on species) prevented us from swabbing a specific site on the body, so we opted to uniformly swab the entire body (head, mouth, tail, dorsal and ventral surfaces) in an effort to capture the total skin microbiome of each amphibian. All individuals were swabbed according to a standard protocol such that the swab was brushed in the direction from the head to the tail and the swab was turned slightly with each move to a new body surface (for example, dorsal stroke, rotate swab, then lateral body stroke) in order to capture each body surface on the swab without saturating the whole swab with one body part; this technique was manageable given the small size of the animals. Each swab was placed in a sterile vial and stored on dry ice. Within a few hours, each sample was then transferred to a −20 °C freezer for storage until DNA extraction. All animals were handled and released according to an approved Institutional Animal Care and Use Committee protocol.

DNA extraction, bacterial 16S rRNA PCR amplification and pyrosequencing

Bacterial DNA was extracted from each swab using the PowerSoil DNA isolation kit (MoBio Laboratories, Carlsbad, CA, USA). For each of the 32 samples, the swab was loaded into a bead tube from the extraction kit and the remaining procedures were performed according to the manufacturer's protocol. Bacterial community composition was determined using barcoded pyrosequencing, following the procedure described in detail in Fierer , which includes both the PCR conditions and error-correcting barcoded primer sequences. The primer set used was the 27F/388R that targets the V2 region of the 16S small subunit ribosomal gene. No template and template from blank filters were included as negative controls throughout the entire process from DNA extraction to PCR amplification to check for contamination. PCR amplicons were pooled from each sample to equimolar ratio in a single tube, and sequenced at the University of South Carolina Environmental Genomics Core Facility on a 454 Life Sciences Sequence GS FLX instrument (Roche, Florence, SC, USA). We obtained an average of 1220 sequences per sample (range 780–1510), with read lengths averaging 258 bp in length. Previous work has shown that reads of this length over the V2 region are sufficient for accurate taxonomic identification to at least the family level and for community characterization (Liu ).

A portion of the extracted DNA was used to screen samples for Bd, the chytrid fungal pathogen. Frozen DNA was shipped to Pisces Molecular LLC (Boulder, CO, USA) for analysis using Bd-specific primers (Annis ), PCR amplification and gel electrophoresis to examine for the presence or absence of bands.

Sequence analysis

All sequence analyses were conducted using the Quantitative Insights Into Microbial Ecology (QIIME) pipeline (Caporaso ). In brief, the QIIME pipeline removed low-quality sequences (that is, those sequences <200 bp in length, had a quality score <25, contained ambiguous characters, having an unreadable barcode or did not contain the primer sequence). Remaining sequences were clustered using CD-HIT (Li and Godzik, 2006), using a minimum coverage of 99% and minimum sequence identity of 97%. A representative sequence was chosen for each phylotype by selecting the longest sequence with the highest number of hits to other sequences in the same phylotype. The representative sequences were aligned using PyNAST (Caporaso ) and taxonomy was assigned using the RDP classifier (Wang ). A phylogenetic tree of all aligned sequences was constructed with FastTree (Price ) and used in all downstream phylogenetic community comparisons. Community similarity (that is, β-diversity) analyses were performed using the taxon-based Bray–Curtis dissimilarity metric, and the phylogenetic UniFrac metric: the Bray–Curtis metric examines taxon overlap, whereas the UniFrac metric is a phylogenetic approach based on the fraction of shared branch length between two communities within a phylogenetic tree (Lozupone ). Phylogenetic trees were constructed using FastTree (Price ). All samples were rarified to 750 sequences per sample to remove sample heterogeneity, which impacts α- and β-diversity metrics (Lozupone ).

Statistical analyses

To address whether bacterial phylotype richness varied among species and/or among pond sites, we conducted an analysis of covariance. Phylotype richness per frog served as the response variable (n=32) and frog species and pond site served as two categorical covariables. This analysis was performed using the statistical software JMP (SAS version 8, Cary, NC, USA). The data were normally distributed as determined by a Shapiro–Wilk W goodness-of-fit test.

To determine the differences in bacterial community composition across individuals (β-diversity), we used two different distance metrics: a Bray–Curtis similarity index (a taxonomic metric) and the UniFrac algorithm, which provides a measure of phylogenetic distance between communities from individual samples (Lozupone ). β-Diversity patterns were visualized using a MDS (non-metric multidimensional scaling) ordination approach (Clarke and Warwick, 2001). To test for significant differences in community composition, we used a two-way nested ANOSIM (analysis of similarity) to examine whether host amphibian species and/or pond site were significant predictors of variability in bacterial communities across individual amphibian samples (PRIMER-E). We ran two separate analyses, one with Bray–Curtis similarity as the response variable and one with unweighted UniFrac distances as the response variable.

Lastly, to visualize the relative abundances of the most common bacterial groups across amphibian species per pond, we constructed a heat map figure. By removing all sequences that comprised ⩽5% of the total number of classified sequences across all samples, 19 phylotypes remained. By reducing the data set to the 19 most abundant phylotypes, visualization of the compositional differences was tractable and issues related to sequencing error are avoided (Quince ; Huse ; Kunin ).

Results

From the rarefied 750 sequences per sample (n=32), we found between 15 and 345 unique phylotypes per sample representing 18 different bacterial phyla. Diversity results thus apply only to the dominant members of the community, and it is possible that sequences scored as absent are in fact present at lower abundance; however, the communities are substantially different from one another when sampled at this depth. The dominant phyla across all samples were Acidobacteria, Actinobacteria, Bacteriodetes, Cyanobacteria, Firmicutes and Proteobacteria. A list of all phylotypes and their taxonomic affiliation is provided in Supplementary Information (Supplementary Table S1).

Phylotype richness varied between the amphibian species (Figure 1). Northern leopard frogs had the highest bacterial phylotype richness relative to the first 750 sequences sampled (mean=287 phylotypes per sample), western chorus frogs had a moderate richness (mean=185) and tiger salamanders had the lowest (mean=66). These phylotype richness patterns varied significantly among species (P<0.0001, r2=0.60) but were not significantly different among ponds (P=0.47), as determined by an analysis of covariance. We cannot effectively estimate absolute levels of diversity because sequencing errors may lead to an inflation of phylotype numbers (Quince ) and because deeper sequencing depths would undoubtedly increase the number of phylotypes identified. However, since all samples were compared at an equivalent sequencing depth, these results do suggest that there are strong differences in relative levels of bacterial diversity across the amphibian species sampled here.

Figure 1

Phylotype richness per amphibian species. The number of unique phylotypes per individual amphibian, estimated from a sampling depth of 750 sequences per sample. Each symbol represents the estimate of unique phylotypes or operational taxonomic units (OTUs) per amphibian individuals at a given site. Different symbols correspond to different pond locations.

There were strong differences in the composition of the skin-associated bacterial communities found on the different amphibian species, a pattern evident in the ordination plot (Figure 2). These patterns are confirmed by the ANOSIM analyses; with both the taxon-based Bray–Curtis distance and the unweighted UniFrac phylogenetic distance metrics, host species identity was a highly significant predictor of skin bacterial community uniqueness (P<0.001; Table 2). Host species explained a large proportion of the variance in the skin bacterial communities (86% using Bray–Curtis and 91% using unweighted UniFrac). Different pond sites were not significant predictors of skin bacterial communities on amphibians (Table 2), although variation among sites within species does exist, as evident from Figure 2.

Figure 2

Ordination plot. β-Diversity patterns were visualized using a MDS ordination approach with skin-associated bacterial community differences represented as Bray–Curtis distances. The 2D stress (or measure of distortion) in the MDS configuration was relatively low (0.13), so the 2D distance between points in the ordination plot is a good representation of the degree of similarity between each sample's bacterial community. Each point represents the skin bacterial community of an individual amphibian; color indicates species and shape indicates pond location.

Table 2

Results of the two-way nested analysis of similarity (ANOSIM)

	Bray–Curtis similarity		Unweighted UniFrac distance
	Global R	Significance	Global R	Significance
Between amphibian species	0.868	<0.001	0.914	<0.001
Between pond sites	0.089	0.391	−0.178	0.714

We examined the differences among bacterial communities from different amphibian species nested within different pond sites using two different community similarity measurements.

The dominant phylotypes across all amphibian samples and their relative abundances are presented in Figure 3, with further details provided in Supplementary Table S1. The most dominant phylotype across all three amphibian species was classified as belonging to the genus Curvibacter (Betaproteobacteria: Comamonadaceae), which accounted for 12–97% of the sequences from individual samples, depending on the species and pond location (Figure 3). A Bacteroidetes phylotype belonging to the group Sphingobacteriales was dominant on tiger salamanders from one pond, comprising 51% of the community, but otherwise, the remaining phylotypes observed on salamanders comprised ⩽6% of the bacterial sequences. A species of Ideonella (Betaproteobacteria) was abundant on chorus frogs from one pond, comprising 18% of the community and another phylotype belonging to the Actinobacteria (Actinomycetales) was also abundant on chorus frogs at two of the pond sites. Leopard frogs had many phylotypes present at lower abundances. For instance, several phylotypes comprised between 3% and 8% of the leopard frog skin communities, and these phylotypes belonged to a range of different phyla and subphyla.

Figure 3

Heat map. A color-scale heat map demonstrates the relative abundance of the 19 most common bacterial phylotypes across species and ponds; red cells indicate higher proportional abundance and blue cells indicate lower proportional abundance. The relative abundance of each phylotype was converted to a proportion of the total number of classified sequences. The total number of unclassified sequences summed together for each amphibian sample comprised on average only 5.7% (±s.d.=0.06) of the sequences per sample and were removed from the table. In order to condense the bacterial groups and present the most common phylotypes in the heat map, we further removed all phylotypes that comprised ⩽5% of the total number of classified sequences across all samples. The proportion of sequences that are represented in the heat map for each species in each pond site varied between 55–98%. Nineteen unique phylotypes remained and the average proportion was calculated per amphibian species per pond.

One individual western chorus frog (P. triseriata) from Chive Pond was found to be positive for Bd, the chytrid fungal skin pathogen. Due to a lack of other positive individuals in this data set, we cannot draw any conclusions about the possible interactions of chytrid infection with the skin bacterial communities.

Discussion

Our understanding of the microbes associated with plants and animals, including humans, is in a phase of rapid discovery owing, in part, to the development of high-throughput sequence-based approaches that allow us to document the structure and composition of numerous complex bacterial communities simultaneously. With improved community-level data, we can explore questions regarding how microbiomes compare across individuals, species or environments, yet there is a paucity of studies comparing microbiomes across animal species. Ley compared gut microbiomes of 60 different mammal species and discovered that both host phylogeny and diet appear to influence the bacterial communities. Roughly half of the samples they observed were from animals in the wild, but from populations that did not overlap in their geographic range, and the other half were from captive zoo animals. Yildirim compared the gut microbiome of three different primate species from wild troops whose home ranges overlapped. While the authors found highly significant differences in the gut bacteria among the primate species, it is difficult to completely disentangle the role of social interactions (for example, grooming, mating, aggressive interactions) and close physical contact. It is possible that the similarity of the gut communities of primate individuals is, in part, due to closeness and physical contact within a troop rather than phylogenetic host differences. Thus, many wild animal species may not directly overlap in their habitat use or level of social contact with other species, and these factors may confound our understanding of species-specific differences of their microbiomes.

As the skin microbiome is in more constant contact with the external environment than the gut microbiome, we might expect this to obscure patterns of phylogenetic host specificity. For example, we may expect the skin microbiome to consist of many exogenous bacteria originating from the environment (Belden and Harris, 2007), rather than bacteria that are host specific. Although the human skin microbiome and its spatial variability has been reasonably well studied (Costello ; Grice ), patterns of microbial assemblages living on the skin of other animals is largely unknown. Without fur or feathers, amphibians provide an excellent model system for studying skin-associated microbial communities. Pre-metamorphic amphibians in pond habitats provide a rare example of multispecies communities that co-exist temporally and spatially in aquatic habitats. We can be confident that all larval individuals present in a given pond during a breeding season originated there and are exposed to the same microbial inocula over the course of their larval development. If there are no differences in skin bacterial communities across species that co-habit in a pond, but there are differences across pond sites, then we can argue that the composition of skin-associated bacterial communities is most strongly influenced by the pond environment. Conversely, if we observe that the skin communities are unique to different amphibian species and that those patterns are consistent across multiple sites, then we can assert that the skin bacterial communities are demonstrating evidence for host specificity. Our results indicate strong support for the latter case, skin bacteria of amphibians appear to be strongly species specific across different sites, even when species co-exist in the same pond (Table 2).

The number of unique skin bacterial phyla varied between 10 and 18 phyla depending upon amphibian species. Using similar techniques, Costello found the human skin microbiome is comprised of 18 different bacterial phyla. Therefore, the phylum-level bacterial richness found on the skin of amphibians appears high and on par with that of human skin. Previous studies have reported a diversity of bacteria from amphibian skin using different methods. Using a PCR-based fingerprinting technique (denaturing gradient gel electrophoresis) of 16S rRNA fragments, Lauer reported up to 19 unique bacteria from eastern red-backed salamanders (Plethodon cinereus), representing four different phyla. Likewise, Woodhams used a culture-based approach and found 40 unique bacterial isolates, representing three different phyla, from 70 frog (Rana muscosa) individuals (mean of eight isolates per frog individual). As in other environments, culture-based approaches likely capture only a small fraction of the bacterial diversity present on amphibian skin (Pace, 1997). While culture-based approaches allow one to determine the specific attributes of individual bacterial taxa, such as their ability to limit fungal growth (Lauer , 2008), combining molecular with culture-based approaches is clearly important for understanding the role of bacterial communities on their amphibian hosts.

All of the phyla reported from previous studies of bacteria on amphibian skin that employed different methods, including Firmicutes, Actinobacteria, Proteobacteria and Bacteroidetes (Culp ; Lauer , 2008; Woodhams ; Lam ), were also found in the present study and comprised some of the more abundant groups observed. Given that the barcoded pyrosequencing method allowed us to observe 14 additional phyla, we can examine community-level patterns and composition to a far greater degree. Closer examination of the composition of the bacterial communities revealed that, for several amphibians sampled in this study, their skin bacterial communities were dominated by a few phylotypes (Figure 3). For example, on tiger salamanders and chorus frogs from four different pond sites, between 42–97% of the bacterial community was dominated by members of the genus Curvibacter, a Betaproteobacteria in the family Comamonadaceae. This was by far the most abundant phylotype across all species in all sites. This genus and its close relatives are commonly found in a wide variety of natural aquatic systems (Shaw ; Mueller-Spitz ; Hutalle-Schmelzer ) and they appear to be primarily freshwater organisms with a wide range of temperature tolerance (9–40 °C) and pH tolerance (5.5–8.5) for growth (Ding and Yokota, 2004). Due to the lack of comparable studies, we do not know if members of this genus are commonly associated with other species of amphibians found at other locations. Another dominant phylotype was identified as being related to Ideonella (Betaproteobacteria) and was particularly abundant on western chorus frogs. The sequence matches that of a symbiont found on the epithelium of Hydra, a small freshwater cnidarian (Fraune ). A few of the other more abundant phylotypes found on the amphibians also match sequences of bacterial strains that are previously known as symbionts living with hosts. These include a strain of Sphingomonas (Alphaproteobacteria) reported to be a symbiont with plant hopper insects (Tang ), a strain of Hydrogenophilus (Betaproteobacteria) reported from the human oral microbiome (Dewhirst ) and a strain of Caulobacter (Alphaproteobacteria) reported from the human gut (Frank ).

Overall, several bacterial types that are related to those found in freshwater systems (for example, Nevskia; Pladdies ) are abundant on amphibian skin, as well as several that are related to known symbionts of other organisms. While we used sterile water rinsing of amphibians specifically to remove transient microbes before swab sampling (sensu Culp ), it is possible that not all of the bacteria obtained from the swab are actual epibiotic symbionts of the amphibians. However, since we found that amphibian species inhabiting the same pond harbored distinct bacterial communities, it is likely that most of the taxa were epibiotic and not simply transient microbes from the pond environment. Furthermore, we do not know whether transient microbes potentially interact with the epibiotic community on amphibians or simply exist temporarily in proximity to the skin microbiome. The mechanism by which amphibian species acquire unique skin communities is not known. Sources of microbial inocula may include vertical transmission (from parent to offspring), horizontal transmission (between individuals) and a variety of environmental sources (water, sediment, algae, invertebrates, organic matter, etc.). While some terrestrial species of amphibians exhibit parental care of egg nests and young (Banning ), the aquatic species studied in this system do not display parental care; thus, vertical transmission seems unlikely in this system. However, Fraune provide evidence that freshwater cnidarian (Hydra) embryos obtain maternally derived peptides that drive the assembly of the epithelial bacterial community, indicating that it may be possible for organisms without extended parental care to receive vertical transmission of peptides that select for the type of community that resides on the offspring. In any case, the assembly of microbial communities that live on amphibians is an intriguing area that can be examined with both field observations and experimental approaches.

While Curvibacter is clearly the most abundant bacteria observed on amphibians in this data set, the presence and relative abundance of other phylotypes largely account for the species-specific patterns observed (Figure 3). Pond sites were not significant in determining the differences among individual amphibians (Table 2); however, we can see that there is some variation across sites; for example, chorus frogs from four different sites have variation in their communities (Figure 3). Three of the Colorado sites occur on the Front Range and one is on the Western Slope (∼200 miles away), yet several dominant phylotypes are still found in both regions (for example, Curvibacter, a Sphingobacteriales phylotype, a Bacteroidaceae phylotype and Sphingomonas; see Figure 3). It is possible to study larger-scale regional variation (for example, continental scale) for some wider ranging species such as leopard frogs, to see if dominant groups shift at different spatial scales; however, these patterns are completely unknown.

Skin-associated bacteria have the potential to influence susceptibility to infectious diseases that affect amphibians. Several studies have cultured bacteria from various amphibians and demonstrated that certain isolates can inhibit the growth of fungi, including Bd (Harris , 2009b; Lauer , 2008; Woodhams ). Harris used an anti-Bd isolate (Janthinobacterium lividum) in a bioaugmentation experiment and demonstrated success in reducing mortality rates due to Bd for susceptible amphibians in the laboratory. The same bacterial strain and other anti-Bd bacteria were found to be more abundant naturally on endangered mountain yellow-legged frog (Rana sierrae) populations that have co-existed with Bd relative to populations that have been extirpated due to Bd (Woodhams ; Lam ). Just as some populations of amphibians differ with respect to their tolerance of Bd, different species of amphibians are known to vary in their ability to tolerate and resist infection with Bd (Garner ; Kilpatrick ). Our finding that amphibians have species-specific bacterial communities on their skin further suggests that these skin bacterial communities have the potential to influence species differences in susceptibility to Bd. Amphibians that harbor more antifungal bacterial species are likely to be better at tolerating or resisting infection with Bd. Broadly, how microbiotas influence the susceptibility to infectious disease is a next frontier with many potential applications relevant to human medicine (Stecher ) and wildlife conservation (Harris ; Lam ).