Abstract
Despite its potential for providing a deeper understanding of both evolutionary processes and epidemic dynamics, the reciprocal interrelationship between mating system evolution and sexually transmitted diseases remains largely unexplored. Here we developed an agent-based model simulating the evolution of two different female mating strategies (monandry versus polyandry) under the spread of a hypothetical sexually transmitted disease implying reproductive costs. Our results strongly support the existence of feedback loops between the pathogen's transmissibility and the evolution of mating strategies. Importantly, we found several unexpected, non-linear emerging behaviours of the system, as well as tipping points which were undocumented under the conditions we considered here: i) medium-high/high probabilities of disease transmission per sexual contact reveal switches between disease-free and endemic outcomes, ii) counterintuitively, the disappearance of the pathogen might be a good indicator and predictor of the imminent extinction of the polyandrous genotype/phenotype/strategy from the population, and iii) probabilities of transmission above a medium-low threshold can offset the spread of polyandry even when this behaviour entails pronounced benefits. More broadly, our results illustrate that research into reciprocal influences between the dynamics of disease spread and sexual behaviour can provide valuable insights into disease transmission and the evolution of reproductive strategies, as well as into the sensitivity of mating systems' evolution to small variations in ecological and disease contexts.