amphibian skin‐associated microbiome across species, space and life history stages

Skin‐associated bacteria of amphibians are increasingly recognized for their role in defence against pathogens, yet we have little understanding of their basic ecology. Here, we use high‐throughput 16S rRNA gene sequencing to examine the host and environmental influences on the skin microbiota of the cohabiting amphibian species Anaxyrus boreas, Pseudacris regilla, Taricha torosa and Lithobates catesbeianus from the Central Valley in California. We also studied populations of Rana cascadae over a large geographic range in the Klamath Mountain range of Northern California, and across developmental stages within a single site. Dominant bacterial phylotypes on amphibian skin included taxa from Bacteroidetes, Gammaproteobacteria, Alphaproteobacteria, Firmicutes, Sphingobacteria and Actinobacteria. Amphibian species identity was the strongest predictor of microbial community composition. Secondarily, within a given amphibian species, wetland site explained significant variation. Amphibian‐associated microbiota differed systematically from microbial assemblages in their environments. Rana cascadae tadpoles have skin bacterial communities distinct from postmetamorphic conspecifics, indicating a strong developmental shift in the skin microbes following metamorphosis. Establishing patterns observed in the skin microbiota of wild amphibians and environmental factors that underlie them is necessary to understand skin symbiont community assembly, and ultimately, the role skin microbiota play in the extended host phenotype including disease resistance.