Abstract
Many organisms synthesize secondary metabolites against natural enemies. However, to which environmental factors the production of these metabolites is adjusted to is poorly investigated in animals, especially so in vertebrates. Bufadienolides are steroidal compounds that are present in a wide range of plants and animals and, if present in large quantities, can provide protection against natural enemies, such as pathogens. In a correlative study involving 16 natural populations we investigated how variation in bufadienolide content of larval common toads (Bufo bufo) is associated with the bacterial community structure of their aquatic environment. We also evaluated pond size, macrovegetation cover, and the abundance of predators, conspecifics and other larval amphibians. We measured toxin content of tadpoles using HPLC-MS and determined the number of bufadienolide compounds (NBC) and the total quantity of bufadienolides (TBQ). AICc-based model selection revealed strong relationships of NBC and TBQ with bacterial community structure of the aquatic habitat as well as with the presence of conspecific tadpoles. The observed relationships may have arisen due to adaptation to local bacterial communities, phenotypic plasticity, differential biotransformation of toxin compounds by different bacterial communities, or a combination of these processes. Bacterial groups that contribute to among-population variation in toxin content remain to be pinpointed, but our study suggesting that toxin production may be influenced by the bacterial community of the environment represents an important step towards understanding the ecological and evolutionary processes leading to microbiota-mediated variation in skin toxin profiles of aquatic vertebrates.