Abstract
Wolbachia is the most frequent bacterial endosymbiont of arthropods and nematodes. Although it is mostly vertically transmitted, from parent to offspring through the egg cytoplasm, horizontal transfer of Wolbachia is thought to be common over evolutionary timescales. However, the relative frequency of each transmission mechanism has not been studied systematically in closely related species. Additionally, while Wolbachia is generally regarded as a reproductive manipulator, it is unclear how frequently the symbiont induces such effects. In this study, we investigated the presence, phenotypes and phylogenetic relationships among Wolbachia strains in whole genome sequence data for 18 European butterfly sister‐species pairs. We find that sister‐species share Wolbachia strains more often than random species pairs and that the probability of strain sharing is higher for younger pairs of host species, especially those with greater range overlap. We also find that split times between Wolbachia strains that infect the same sister‐species pair generally pre‐date host divergence, ruling out co‐divergence in favour of horizontal transfer. However, some strains are younger than the mitochondrial split times of their hosts, so introgressive transfer cannot be ruled out in some cases. In addition, all newly assembled Wolbachia genomes contained putative homologues of genes associated with cytoplasmic incompatibility and male killing. This supports the potential for reproductive manipulation in Wolbachia strains infecting European butterflies, which until now was only inferred from mitochondrial diversity patterns. Our results show that horizontal and introgressive transfer of Wolbachia are frequent even between recently speciated host taxa, suggesting the symbiont's turnover rate is higher than had been inferred previously from surveys of distantly related hosts.