Exploring the ecological drivers of bacteriophage diversity and functional viral potential in the skin of the axolotl Ambystoma altamirani

Bacteriophages play important roles in shaping microbial community dynamics across diverse environments. In the amphibian skin, most microbiome studies have focused on bacteria and their interactions with the fungus Batrachochytrium dendrobatidis (Bd), leaving other microbial components, including viruses, largely unexplored. Here, we present the first characterization of the viral community in the amphibian skin microbiome, focusing on ecological drivers of bacteriophage diversity and functional potential in the axolotl Ambystoma altamirani. Using public shotgun metagenomes, we found that the viral fraction was dominated by bacteriophages of the class Caudoviricetes. Bacteriophage diversity was significantly associated with local physicochemical parameters at the time of sampling, and showed a strong positive correlation with bacterial diversity, whereas no significant associations were detected with the presence of Bd. In addition, seasonality influenced the composition and properties of bacteria-bacteriophage co-abundance networks. Functional annotation of assembled bacteriophage sequences revealed a diverse functional potential, including putative auxiliary metabolic genes, superinfection exclusion, toxin-antitoxin, and virulence factors. Overall, these findings highlight the ecological relevance of bacteriophages in amphibian skin microbiomes and underscore the need for further studies on their role in the amphibian host's health.