Home-loving boreal hare mitochondria survived several invasions in Iberia: the relative roles of recurrent hybridisation and allele surfing.

Genetic introgression from a resident species into an invading close relative can result from repeated hybridisation along the invasion front and/or allele surfing on the expansion wave. Cases where the phenomenon is massive and systematic, such as for hares (genus Lepus) in Iberia, would be best explained by recurrent hybridisation but this is difficult to prove because the donor populations are generally extinct. In the Pyrenean foothills, Lepus europaeus presumably replaced Lepus granatensis recently and the present species border is parallel to the direction of invasion, so that populations of L. granatensis in the contact zone represent proxies of existing variation before the invasion. Among three pairs of populations sampled across this border, we find less differentiation of mitochondrial DNA (mtDNA) across than along it, as predicted under recurrent hybridisation at the invasion front. Using autosomal microsatellite loci and X- and Y-linked diagnostic loci, we show that admixture across the border is quasi-absent, making it unlikely that lack of interspecific mtDNA differentiation results from ongoing gene flow. Furthermore, we find that the local species ranges are climatically contrasted, making it also unlikely that ongoing ecology-driven movement of the contact account for mtDNA introgression. The lack of mtDNA differentiation across the boundary is mostly due to sharing of mtDNA from a boreal species currently extinct in Iberia (Lepus timidus) whose mitochondria have thus remained in place since the last deglaciation despite successive invasions by two other species. Home-loving mitochondria thus witness past species distribution rather than ongoing exchanges across stabilised contact zones.

Recurrent introgression of mitochondrial DNA among hares (Lepus spp.) revealed by species-tree inference and coalescent simulations.

Understanding recent speciation history requires merging phylogenetic and population genetics approaches, taking into account the persistence of ancestral polymorphism and possible introgression. The emergence of a clear phylogeny of hares (genus Lepus) has been hampered by poor genomic sampling and possible occurrence of mitochondrial DNA (mtDNA) introgression from the arctic/boreal Lepus timidus into several European temperate and possibly American boreal species. However, no formal test of introgression, taking also incomplete lineage sorting into account, has been done. Here, to clarify the yet poorly resolved species phylogeny of hares and test hypotheses of mtDNA introgression, we sequenced 14 nuclear DNA and 2 mtDNA fragments (8205 and 1113 bp, respectively) in 50 specimens from 11 hare species from Eurasia, North America, and Africa. By applying an isolation-with-migration model to the nuclear data on subsets of species, we find evidence for very limited gene flow from L. timidus into most temperate European species, and not into the American boreal ones. Using a multilocus coalescent-based method, we infer the species phylogeny, which we find highly incongruent with mtDNA phylogeny using parametric bootstrap. Simulations of mtDNA evolution under the speciation history inferred from nuclear genes did not support the hypothesis of mtDNA introgression from L. timidus into the American L. townsendii but did suggest introgression from L. timidus into 4 temperate European species. One such event likely resulted in the complete replacement of the aboriginal mtDNA of L. castroviejoi and of its sister species L. corsicanus. It is remarkable that mtDNA introgression in hares is frequent, extensive, and always from the same donor arctic species. We discuss possible explanations for the phenomenon in relation to the dynamics of range expansions and species replacements during the climatic oscillations of the Pleistocene.