Invasion from the cold past: extensive introgression of mountain hare (Lepus timidus) mitochondrial DNA into three other hare species in northern Iberia

Mitochondrial DNA introgression from Lepus timidus into Lepus granatensis and Lepus europaeus was recently reported in Iberia, although L. timidus presumably retreated from this region at the end of the last ice age. Here we assess the extent of this ancient mtDNA introgression by RFLP analysis of 695 specimens representing the three hare species present in Iberia. The introgressed L. timidus lineage was found in 23 of the 37 populations sampled. It is almost fixed in L. europaeus across its Iberian range in the Pyrenean foothills, and in L. granatensis, which occupies the rest of the peninsula, it is predominant in the north and gradually disappears further south. We also found it in Lepus castroviejoi, a species endemic to Cantabria. Multiple hybridizations and, potentially, a selective advantage for the L. timidus lineage can explain the remarkable taxonomic and geographical range of this mitochondrial introgression.     

The ubiquitous mountain hare mitochondria: multiple introgressive hybridization in hares, genus Lepus

Climatic oscillations during the glaciations forced dramatic changes in species distributions, such that some presently temperate regions were alternately occupied by temperate and arctic species. These species could have met and hybridized during climatic transitions. This phenomenon happened for three hare species present in Iberia (Lepus granatensis, Lepus europaeus and Lepus castroviejoi), which display high frequencies of mitochondrial DNA (mtDNA) from Lepus timidus, an arctic/boreal species presently extinct in Iberia. Here, we extend our previous geographical survey to determine whether the distribution of this mtDNA lineage extends beyond the northern half of the Iberian Peninsula, where it is found at high frequencies. We also review the taxonomy, distribution and molecular phylogeny of the genus Lepus. The phylogenetic inference reveals the presence of L. timidus-like mtDNA in several other hare species in Asia and North America, suggesting that the mitochondrial introgression observed in Iberia might be generalized. Comparison with the available nuclear gene phylogenies suggests that introgression could have happened repeatedly, possibly during different climatic transitions. We discuss demographic and adaptive scenarios that could account for the repetition in time and space of this spectacular phenomenon and suggest ways to improve our understanding of its determinants and consequences. Such high levels of introgressive hybridization should discourage attempts to revise hare taxonomy based solely on mtDNA.     

The diel activity pattern of mountain hare (Lepus timidus) on managed heather moorland in Scotland

The research presented in this paper provides an insight into the behavioral ecology of mountain hares on heather moorland in the Lammermuir Hills of southeast Scotland. We examine the seasonal and diel activity patterns using camera traps over a period of 12 months. The rate of camera detections was calculated for the different divisions of the 24‐hr cycle (daylight, dusk, night, and dawn). During autumn and winter (October–February), the activity pattern was crepuscular with greater activity at dusk than at dawn. Daylight activity was relatively low, and there was a regular pattern of small peaks of activity during the night. In spring and summer (March–September), peaks of crepuscular activity remained evident but daylight activity was much more prevalent than during autumn and winter, and night activity was lower. We discuss the problematic definition of twilight and present an explanation for seasonal changes in the pattern of diel activity that is linked to the reproductive cycle of the mountain hare.     

Diet of the Iberian hare (Lepus granatensis) in a mountain ecosystem

The diet of the Iberian hare (Lepus granatensis) was studied through microhistological pellet analysis in two areas from a mountain ecosystem in Central Portugal. Fecal pellets were collected monthly in 24 plots spatially distributed throughout the two study areas. For each period, a sample of 15 to 20 pellets was milled and 400 epidermal fragments were identified, by comparison with a reference collection. A wide range of plant species was observed in hare’s diet. Grasses represent the basis of the Iberian hare diet, with frequencies always higher than 50% in both study areas (annual average = 69.98%). Most of the 35 species of grasses assembled for the reference collection (91.43%) were identified in the pellets. Nevertheless, only six of these were consumed in proportions greater than 5%, being Anthoxanthum odoratum, Secale cereale and Agrostis spp. the species ingested in higher frequencies. The rate of grasses consumption reached 80.69% in winter but decreased in summer to around 55%. In this season, a concurrent rise in the ingestion of other plant groups, like herbs and shrubs, and of plant inflorescences was observed. This work provides the first results on the Iberian hare’s diet on mountain ecosystems, and suggests that the Iberian hare diet in a mountain ecosystem is similar to the observed in L. europaeus and L. timidus.     

A preliminary study on the seasonal body temperature rhythms of the captive mountain hare (Lepus timidus)

The mountain hare (Lepus timidus) is a year-round active herbivore adapted to survive the boreal winter. Captive mountain hares (N = 4) were implanted with intraabdominal thermosensitive loggers to record their core body temperature (T_b) for a year and during food deprivation (8–48 h) in summer and winter. The average T_b was 38.7 ± 0.01 °C in summer and 38.3 ± 0.01 °C in winter. The yearly T_b correlated positively with the ambient temperature. The 24-h T_b was the highest from late scotophase to early photophase in summer and winter and the lowest during middle-late photophase in summer or during early-middle scotophase in winter. The range of the 24-h oscillations in T_b increased in three animals in winter. Food deprivation did not induce hypothermia in summer or winter. These preliminary data suggest that the mountain hare can spare a modest amount of energy with the wintertime reduction in T_b.     

The occurrence of mountain hare mitochondrial DNA in wild brown hares

If interspecific hybrids are fertile and backcross to either parental species, transmission of mitochondrial DNA over the species barrier can occur. To investigate if such transmission has occurred between the brown hare Lepus europeus Pall and the mountain hare L. timidus L. in Scandinavia, an analysis of genetic variation in mitochondrial DNA from 36 hares, collected from 15 localities, was performed. Sequence divergence of mtDNA between species was estimated at 8 ± 1% (SD). Intraspecific mtDNA sequence divergence varied between 0.09 and 0.38% in brown hares and 0.10 and 1.44% in mountain hares. In six out of 18 brown hares examined, two different haplotypes of mountain hare origin were detected, demonstrating a transmission of mtDNA haplotypes from mountain hares to brown hares. The results indicate that interspecific hybridization between the two species occurs in wild populations.     

The effect of establishing native woodland on habitat selection and ranging of moorland mountain hares (Lepus timidus), a flexible forager

Mountain hares Lepus timidus L. typify species that occupy a broad geographic range and have flexible foraging and nutritional strategies. Such species may show a range of responses to habitat modification. This study aimed to provide a basis for prediction of the impact of mountain hares on woodland establishment, and of woodland establishment on mountain hare distribution. The selection of and the extent of incorporation of new woodland into the home range of mountain hares was investigated in an area where Scots pine Pinus sylvestris L. woodland was establishing within their usual habitat in Britain, upland heather moorland. Seasonal home, day and night-range sizes of radio-tracked mountain hares were determined using the multinuclear probability polygon technique and analysed using residual maximum likelihood (REML). Habitat selection was analysed using compositional analysis. Three main habitat types were available to hares: heather moorland with trees, heather moorland and grassland-mire. Mean home-range size of mountain hares in summer was 10.3 ha and in winter 9.6 ha. There were no significant seasonal or sex differences in home-range size. Females selected grassland-mire habitat in summer and showed no strong selection for any habitat in winter. Males selected heather moorland in both summer and winter. Heather moorland with trees was not selected preferentially by mountain hares of either sex in summer or winter. The absence of selection for areas of newly establishing-Scots pine woodland suggests that any browsing damage to trees by hares is most likely to be a function of the local abundance of mountain hares, rather than a result of active preference of hares for the modified habitat.     

Phylogeography of the brown hare (Lepus europaeus) in Europe: a legacy of south‐eastern Mediterranean refugia?

Aim We analysed the population genetics of the brown hare (Lepus europaeus) in order to test the hypothesis that this species migrated into central Europe from a number of late glacial refugia, including some in Asia Minor. Location Thirty‐three localities in Greece, Bulgaria, Italy, Croatia, Serbia, Poland, Switzerland, Austria, France, Germany, the Netherlands, Spain, the United Kingdom, Turkey and Israel. Methods In total, 926 brown hares were analysed for mitochondrial DNA (mtDNA) variation by restriction fragment length polymorphism (RFLP) performed on polymerase chain reaction‐amplified products spanning cytochrome b (cyt b)/control region (CR), cytochrome oxidase I (COI) and 12S–16S rRNA. In addition, sequence analysis of the mtDNA CR‐I region was performed on 69 individuals, and the data were compared with 137 mtDNA CR‐I sequences retrieved from GenBank. Results The 112 haplotypes detected were partitioned into five phylogeographically well‐defined major haplogroups, namely the ‘south‐eastern European type haplogroup’ (SEEh), ‘Anatolian/Middle Eastern type haplogroup’ (AMh), ‘European type haplogroup, subgroup A’ (EUh‐A), ‘European type haplogroup, subgroup B’ (EUh‐B) and ‘Intermediate haplogroup’ (INTERh). Sequence data retrieved from GenBank were consistent with the haplogroups determined in this study. In Bulgaria and north‐eastern Greece numerous haplotypes of all five haplogroups were present, forming a large overlap zone. Main conclusions The mtDNA results allow us to infer post‐glacial colonization of large parts of Europe from a late glacial/early Holocene source population in the central or south‐central Balkans. The presence of Anatolian/Middle Eastern haplotypes in the large overlap zone in Bulgaria and north‐eastern Greece reveals gene flow from Anatolia to Europe across the late Pleistocene Bosporus land‐bridge. Although various restocking operations could be partly responsible for the presence of unexpected haplotypes in certain areas, we nevertheless trace a strong phylogeographic signal throughout all regions under study. Throughout Europe, mtDNA results indicate that brown hares are not separated into discernable phyletic groups.     

The genomic legacy from the extinct Lepus timidus to the three hare species of Iberia: contrast between mtDNA, sex chromosomes and autosomes

Extensive interspecific genetic introgression is often reported, and appraising its genomic impact can serve to determine whether it results from selection on specific loci or from demographic processes affecting the whole genome. The three species of hares present in the Iberian Peninsula harbour high frequencies of mitochondrial DNA (mtDNA) from Lepus timidus, an arctic/boreal species now extinct in the region. This could result from the invasive replacement of L. timidus by the temperate species during deglaciation but should then have left traces in the nuclear genome. We typed single nucleotide polymorphisms (SNPs) discovered by sequencing 10 autosomal loci, two X-linked and one Y-linked in species-wide samples of the four taxa. Based on lineage-diagnostic SNPs, we detected no trace of L. timidus sex chromosomes in Iberia. From the frequencies of inferred haplotypes, autosomal introgression into L. granatensis appeared mostly sporadic but always widespread instead of restricted to the north as mtDNA. Autosomal introgression into Iberian L. europaeus , inhabiting the Pyrenean foothills, was hardly detectable, despite quasi-fixation of L. timidus mtDNA. L. castroviejoi , endemic to the Cantabrian Mountains and fixed for L. timidus mtDNA, showed little traces of autosomal introgression. The absence of sex-chromosome introgression presumably resulted from X-linked hybrid male unfitness. The contrasting patterns between the autosomes and mtDNA could reflect general gender asymmetric processes such as frequency-dependent female assortative mating, lower mtDNA migration and higher male dispersal, but adaptive mtDNA introgression cannot be dismissed. Additionally, we document reciprocal introgression between L. europaeus and both L. granatensis in Iberia and L. timidus outside Iberia.     

Steppes, savannahs, forests and phytodiversity reservoirs during the Pleistocene in the Iberian Peninsula

A palaeobotanical analysis of the Pleistocene floras and vegetation in the Iberian Peninsula shows the existence of patched landscapes with Pinus woodlands, deciduous and mixed forests, parklands (savannah-like), shrublands, steppes and grasslands. Extinctions of Arctotertiary woody taxa are recorded during the Early and Middle Pleistocene, but glacial refugia facilitated the survival of a number of temperate, Mediterranean and Ibero-North African woody angiosperms. The responses of Iberian vegetation to climatic changes during the Pleistocene have been spatially and temporarily complex, including rapid changes of vegetation in parallel to orbital and suborbital variability, and situations of multi-centennial resilience or accommodation to climatic changes. Regional characteristics emerged as soon as for the Middle Pleistocene, if not earlier: Ericaceae in the Atlantic coast indicating wetter climate, thermo-mediterranean elements in the south as currently, and broad-leaf trees in the northeastern. Overall, steppe landscapes and open Pinus woodlands prevailed over many continental regions during the cold spells of the Late Pleistocene. The maintenance of a high phytodiversity during the glacials was linked to several refuge zones in the coastal shelves of the Mediterranean and intramountainous valleys. Northern Iberia, especially on coastal areas, was also patched with populations of tree species, and this is not only documented by palaeobotanical data (pollen and charcoal) but also postulated by phylogeographical models.     

The diet of the mountain hare (Lepus timidus hibernicus) on coastal grassland

Across most of their range in Europe, mountain hares are usually restricted to upland areas with poor food quality. In these areas they generally feed on browse species such as heather or twigs and barks of trees. On lowland areas in Europe, with better food quality, the mountain hare is replaced by the brown hare ( Lepus europaeus ) which feeds predominantly on greasses. This khas led some authors to conclude that mountain hares are primarily adapted for browsing. In the absence of brown hares in Ireland, mountain hares are found on a wide variety of habitats including grassland. On grassland, their diet consists almost exclusively of grasses, up to 94% of their annual diet, which is more than has been reported for brown hares on similar habitat. Based on this evidence, and other work, it is proposed that the mountain hare in primarily a grazing animal and competitive exclusion by brown hares may underlie much of their present distribution in Europe.     

Contagious mucocutaneous dermatitis of the mountain hare (Lepus timidus): pathology and cause

Contagious mucocutaneous dermatitis is a frequently encountered disease of mountain hares (Lepidus timidus) in Finland. We describe the histopathologic changes and propose an etiologic cause for this disorder. Fifty-three cases collected during 1982-2000 were examined histologically. Transmission electron microscopy was performed in one case. In fully developed lesions, keratinocytes in epidermis and follicular infundibula were swollen and contained large eosinophilic intracytoplasmic inclusion bodies with marked reticular and ballooning degeneration. In later stages, there was marked necrosis and ulceration with severe pyogranulomatous and suppurative inflammation. At this stage, no viral inclusions were detectable, but secondary Staphylococcus warnerii infection was present in most cases. In late lesions, there was dermal fibrosis with epidermal hyperplasia. No spiral-shaped bacteria suggesting treponematosis were detected at any stage. Ultrastructurally, swollen epidermal and follicle infundibular cells contained round intracytoplasmic inclusion bodies with a myriad of virions typical of poxvirus with a biconcave nucleocapsid core, two lateral bodies, and a clearly discernible outer lipoprotein capsule. The findings suggest that contagious mucocutaneous dermatitis in mountain hares is a viral disease caused by a poxvirus. The disease is often complicated by secondary bacterial infection, most commonly S. warneri.     

Mitochondrial DNA phylogeography of Lissotriton boscai (Caudata, Salamandridae): evidence for old, multiple refugia in an Iberian endemic

In Europe, southern peninsulas served as refugia during cold periods in the Pleistocene, acting both as centres of origin of endemisms and as sources from which formerly glaciated areas were recolonized during interglacial periods. Previous studies have revealed that within the main refugial areas, intraspecific lineages often survived in allopatric refugia. We analysed two mitochondrial markers (nad4, control region, approximately 1.4 kb) in 103 individuals representing the entire distribution of Lissotriton boscai, a newt endemic to the western Iberian Peninsula. We inferred the evolutionary history of the species through phylogenetic, phylogeographic and historical demographic analyses. The results revealed unexpected, deep levels of geographically structured genetic variability. We identified two main evolutionary lineages, each containing three well-supported clades. The first historical split involved populations from central-southwestern coastal Portugal and the ancestor of all the remaining populations around 5.8 million years ago. Both lineages were subsequently fragmented into different population groups between 2.5 and 1.2 million years ago. According to nested clade analysis, at lower hierarchical levels the patterns suggest restricted gene flow with isolation by distance, whereas at higher levels the clades exhibit signatures of contiguous range expansion. Bayesian Skyline Plots show recent bottlenecks, followed by demographic expansions in all lineages. The significant genetic structure found is consistent with long-term survival of populations in allopatric refugia, supporting the 'refugia-within-refugia' scenario for southern European peninsulas. The comparison of our results with other co-distributed species highlights the generality of this hypothesis for the Iberian herpetofauna and suggests that Mediterranean refuges had more relevance for the composition and distribution of present biodiversity patterns than currently acknowledged. We briefly discuss the taxonomic and conservation implications of our results.