Links between prey assemblages and poison frog toxins: a landscape ecology approach to assess how biotic interactions affect species phenotypes

Ecological studies of species pairs showed that biotic interactions promote phenotypic change and eco-evolutionary feedbacks. However, it is unclear how phenotypes respond to synergistic interactions with multiple taxa. We investigate whether interactions with multiple prey species explain spatially structured variation in the skin toxins of the Neotropical poison frog Oophaga pumilio. Specifically, we assess how dissimilarity (i.e., beta diversity) of alkaloid-bearing arthropod prey assemblages (68 ant species) and evolutionary divergence between frog populations (from a neutral genetic marker) contribute to frog poison dissimilarity (toxin profiles composed of 230 different lipophilic alkaloids sampled from 934 frogs at 46 sites). We find that models that incorporate spatial turnover in the composition of ant assemblages explain part of the frog alkaloid variation, and we infer unique alkaloid combinations across the range of O. pumilio. Moreover, we find that alkaloid variation increases weakly with the evolutionary divergence between frog populations. The analytical framework proposed here can be extended to other multi-trophic systems, coevolutionary mosaics, microbial assemblages, and ecosystem services.