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.

Changing from summer-brown to winter-white pelage or plumage is a crucial adaptation to seasonal snow in more than 20 mammal and bird species. Many of these species maintain nonwhite winter morphs, locally adapted to less snowy conditions, which may have evolved independently. Mountain hares (Lepus timidus) from Fennoscandia were introduced into the Faroe Islands in 1855. While they were initially winter-white, within ∼65 y all Faroese hares became winter-gray, a morph that occurs in the source population at low frequency. The documented population history makes this a valuable model for understanding the genetic basis and evolution of the seasonal trait polymorphism. Through whole-genome scans of differentiation and single-nucleotide polymorphism (SNP) genotyping, we associated winter coat color polymorphism to the genomic region of the pigmentation gene Agouti, previously linked to introgression-driven winter coat color variation in the snowshoe hare (Lepus americanus). Lower Agouti expression in the skin of winter-gray individuals during the autumn molt suggests that regulatory changes may underlie the color polymorphism. Variation in the associated genomic region shows signatures of a selective sweep in the Faroese population, suggesting that positive selection drove the fixation of the variant after the introduction. Whole-genome analyses of several hare species revealed that the winter-gray variant originated through introgression from a noncolor changing species, in keeping with the history of ancient hybridization between the species. Our findings show the recurrent role of introgression in generating winter coat color variation by repeatedly recruiting the regulatory region of Agouti to modulate seasonal coat color change.

Introgressive hybridization is an important and widespread evolutionary process, but the relative roles of neutral demography and natural selection in promoting massive introgression are difficult to assess and an important matter of debate. Hares from the Iberian Peninsula provide an appropriate system to study this question. In its northern range, the Iberian hare, Lepus granatensis, shows a northwards gradient of increasing mitochondrial DNA (mtDNA) introgression from the arctic/boreal L. timidus, which it presumably replaced after the last glacial maximum. Here, we asked whether a south-north expansion wave of L. granatensis into L. timidus territory could underlie mtDNA introgression, and whether nuclear genes interacting with mitochondria (“mitonuc” genes) were affected. We extended previous RNA-sequencing and produced a comprehensive annotated transcriptome assembly for L. granatensis. We then genotyped 100 discovered nuclear SNPs in 317 specimens spanning the species range. The distribution of allele frequencies across populations suggests a northwards range expansion, particularly in the region of mtDNA introgression. We found no correlation between variants at 39 mitonuc genes and mtDNA introgression frequency. Whether the nuclear and mitochondrial genomes coevolved will need a thorough investigation of the hundreds of mitonuc genes, but range expansion and species replacement likely promoted massive mtDNA introgression.

Hybridization may often be an important source of adaptive variation, but the extent and long-term impacts of introgression have seldom been evaluated in the phylogenetic context of a radiation. Hares (Lepus) represent a widespread mammalian radiation of 32 extant species characterized by striking ecological adaptations and recurrent admixture. To understand the relevance of introgressive hybridization during the diversification of Lepus, we analyzed whole exome sequences (61.7 Mb) from 15 species of hares (1- 4 individuals per species), spanning the global distribution of the genus, and two outgroups. We used a coalescent framework to infer species relationships and divergence times, despite extensive genealogical discordance. We found high levels of allele sharing among species and show that this reflects extensive incomplete lineage sorting and temporally layered hybridization. Our results revealed recurrent introgression at all stages along the Lepus radiation, including recent gene flow between extant species since the last glacial maximum, but also pervasive ancient introgression occurring since near the origin of the hare lineages. We show that ancient hybridization between northern hemisphere species has resulted in shared variation of potential adaptive relevance to highly seasonal environments, including genes involved in circadian rhythm regulation, pigmentation, and thermoregulation. Our results illustrate how the genetic legacy of ancestral hybridization may persist across a radiation, leaving a long-lasting signature of shared genetic variation that may contribute to adaptation.

Aim A striking case of mitochondrial DNA (mtDNA) introgression occurs in the Iberian Peninsula: the mtDNA of Lepus timidus, a species that inhabited this territory during the Last Glacial Maximum (LGM), is almost fixed in the Iberian range of Lepus europaeus, is very frequent in the northern half of Iberia within Lepus granatensis, and is fixed in Lepus castroviejoi. Here, we used a transdisciplinary approach to understand the ecological drivers of the distribution of this introgressed mtDNA in the Iberian hare species. Location Europe, with an emphasis on the Iberian Peninsula. Methods First, we modelled the climatic niche of L. timidus to predict the favourability for the species during the LGM. Second, we gathered data regarding L. timidus mtDNA introgression for 1137 individuals (139 populations) from the Iberian Peninsula. We tested whether the climatic favourability for L. timidus during the LGM reflects the frequency of the introgressed mtDNA in the carrier hare populations. Results The climatic favourability for the LGM is positively related to the presence of individuals with L. timidus mtDNA, and discriminates better than by chance between individuals with and without this lineage. At the population level, the favourability values relate to the observed frequencies of mtDNA introgression and discriminate between populations with and without L. timidus mtDNA. Main conclusions The results suggest that the geographical distribution of L. timidus mtDNA introgression is linked to L. timidus distribution in Iberia during the LGM. Our study helps to improve our understanding of the role of climate in the reticulated evolutionary history of hares, and is a first step towards understanding the ecological impact of mtDNA introgression on Iberian hares.

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 ecologydriven 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.

Color molts from summer brown to winter white coats have evolved in several species to maintain camouflage year‐round in environments with seasonal snow. Despite the eco‐evolutionary relevance of this key phenological adaptation, its molecular regulation has only recently begun to be addressed. Here, we analyze skin transcription changes during the autumn molt of the mountain hare (Lepus timidus) and integrate the results with an established model of gene regulation across the spring molt of the closely related snowshoe hare (L. americanus). We quantified differences in gene expression among three stages of molt progression—“brown” (early molt), “intermediate,” and “white” (late molt). We found 632 differentially expressed genes, with a major pulse of expression early in the molt, followed by a milder one in late molt. The functional makeup of differentially expressed genes anchored the sampled molt stages to the developmental timeline of the hair growth cycle, associating anagen to early molt and the transition to catagen to late molt. The progression of color change was characterized by differential expression of genes involved in pigmentation, circadian, and behavioral regulation. We found significant overlap between differentially expressed genes across the seasonal molts of mountain and snowshoe hares, particularly at molt onset, suggesting conservatism of gene regulation across species and seasons. However, some discrepancies suggest seasonal differences in melanocyte differentiation and the integration of nutritional cues. Our established regulatory model of seasonal coat color molt provides an important mechanistic context to study the functional architecture and evolution of this crucial seasonal adaptation.

Hares (genus Lepus) provide clear examples of repeated and often massive introgressive hybridization and striking local adaptations. Genomic studies on this group have so far relied on comparisons to the European rabbit (Oryctolagus cuniculus) reference genome. Here, we report the first de novo draft reference genome for a hare species, the mountain hare (Lepus timidus), and evaluate the efficacy of whole-genome re-sequencing analyses using the new reference versus using the rabbit reference genome. The genome was assembled using the ALLPATHS-LG protocol with a combination of overlapping pair and mate-pair Illumina sequencing (77x coverage). The assembly contained 32,294 scaffolds with a total length of 2.7 Gb and a scaffold N50 of 3.4 Mb. Re-scaffolding based on the rabbit reference reduced the total number of scaffolds to 4,205 with a scaffold N50 of 194 Mb. A correspondence was found between 22 of these hare scaffolds and the rabbit chromosomes, based on gene content and direct alignment. We annotated 24,578 protein coding genes by combining ab-initio predictions, homology search, and transcriptome data, of which 683 were solely derived from hare-specific transcriptome data. The hare reference genome is therefore a new resource to discover and investigate hare-specific variation. Similar estimates of heterozygosity and inferred demographic history profiles were obtained when mapping hare whole-genome re-sequencing data to the new hare draft genome or to alternative references based on the rabbit genome. Our results validate previous reference-based strategies and suggest that the chromosome-scale hare draft genome should enable chromosome-wide analyses and genome scans on hares.