Abstract
Müllerian mimicry was proposed to be an example of a coevolved mutualism promoted by population isolation in glacial refugia. This, however, has not been well supported in butterfly models. Here, we use genomic data to test this theory while examining the population genetics behind mimetic diversification in a pair of co-mimetic bumble bees, Bombus breviceps Smith and Bombus trifasciatus Smith. In both lineages, populations were structured by geography but not as much by color pattern, suggesting sharing of color alleles across regions of restricted gene flow and formation of mimicry complexes in the absence of genetic differentiation. Demographic analyses showed mismatches between historical effective population size changes and glacial cycles, and niche modeling revealed only mild habitat retraction during glaciation. Moreover, mimetic subpopulations of the same color form in the two lineages only in some cases exhibit similar population history and genetic divergence. Therefore, the current study supports a more complex history in this comimicry than a simple refugium-coevolution model.