Abstract
Skin microbiomes in amphibians are complex systems that can be influenced by biotic and abiotic factors. In this study, we examined the effect of host species and environmental conditions on the skin bacterial and fungal microbiota of four obligate paedomorphic salamander species, commonly known as axolotls (Ambystoma andersoni, A. dumerilii, A. mexicanum, and A. taylori), all of them endemic to the Trans-Mexican Volcanic Belt. We found that despite their permanent aquatic lifestyle, these species present a host-specific skin microbiota that is distinct from aquatic communities. We identified skin-associated taxa that were unique to each host species and that differentiated axolotl species based on alpha and beta diversity metrics. Moreover, we identified a set of microbial taxa that were shared across hosts with high relative abundances across skin samples. Specifically, bacterial communities were dominated by Burkholderiales and Pseudomonadales bacterial orders and Capnodiales and Pleosporales fungal orders. Host species and environmental variables collectively explained more microbial composition variation in bacteria (R2 = 0.46) in comparison to fungi (R2 = 0.2). Our results contribute to a better understanding of the factors shaping the diversity and composition of skin microbial communities in Ambystoma. Additional studies are needed to disentangle the effects of specific host associated and environmental factors that could influence the skin microbiome of these endangered species.