- Research article
- Open Access
Introgression of mountain hare (Lepus timidus) mitochondrial DNA into wild brown hares (Lepus europaeus) in Denmark
BMC Ecologyvolume 6, Article number: 17 (2006)
In Europe the mountain hare (Lepus timidus) exists in Great Britain, Norway, Sweden, Finland, parts of the Alps and in Eastern Europe, but not in Denmark. Interspecific hybridization has been demonstrated between native Swedish mountain hares and introduced brown hares (Lepus europaeus). During the data collection in a study concerning Danish brown hares we identified 16 hares with a single very divergent haplotype.
Phylogenetic analysis shows that the divergent Danish haplotype is most closely related to the Swedish mountain hare. The frequency of Lepus timidus mtDNA haplotype in the Eastern Danish hare populations is estimated to 6%.
In contrast to what is known, the Danish hare populations are not pure L. europaeus populations but include introgressed brown hares with Swedish L. timidus mtDNA. The most probable explanation of this is natural migration or translocation of introgressed brown hares from Sweden. The impurity of hare populations has implications for conservation and population genetics.
The brown hare (Lepus europaeus Pallas, 1778) is widely distributed throughout Europe, up to 60°N, in Asia Minor and probably south to Israel. It is a popular game animal and its range has expanded to the east both by natural dispersion and by translocations to central and far-east Siberia . Worldwide it has also been imported to many other countries, outside its natural distribution range, e.g. South and North America, Australia and New Zealand.
The brown hare is common all over Denmark and is the only Lepus species that exists in the country. Since 1960 the hare population has declined drastically in Europe and in Denmark. The primary cause of the decline is habitat changes in relation to agricultural intensification, whereas effects of climatic changes and predator abundance have increased by loss of year-round access to high-quality food and cover . Additional factors like predators (especially foxes), birds of prey and traffic may have influenced the population locally . In Denmark this is reflected by a large drop in annual number of hares shot from more than 400,000 before 1960 to 67,600 in 2004/2005 .
The general range of the mountain hare (Lepus timidus Linnaeus, 1758) covers much of the Palaearctic. In Europe it exists in Great Britain, Norway, Sweden, Finland, parts of the Alps and in Eastern Europe but it does not exist in Denmark . During the late nineteenth century the brown hare was introduced in southern Sweden as a game animal . In Scandinavia the mountain hare populations have been retreating northwards since the introduction of the brown hare, presumably as a result of a gradual competitive exclusion by the latter . Mountain hares supposedly colonized Scandinavia through repeated immigration waves over two different post-glacial colonization routes; one from the south and one from the northeast .
Given the sympatric distribution of several hare species, e.g. in Sweden, Italy, and Spain, a number of studies focused on the possibility of interspecific hybridization. The first evidence of interspecific hybridization in hares concerns L. timidus and L. europaeus in Sweden [9–11] where unidirectional introgression of the native L. timidus mtDNA occurs in the introduced L. europaeus. However, introgression of mitochondrial DNA into other hare species has been demonstrated, e.g. the introgression of L. timidus mtDNA into L. granatensis, L. europaeus and L. castroviejoi in northern Iberia  and introgression of foreign mtDNA is likely to occur in several hare species from Asia . These findings have lead to the conclusion that the genetic integrity of many European species largely depends on differences in behaviour and ecology which, at best, offer semi-permeable isolation . In captivity, mountain hare females mate with brown hare males and produce viable offspring, whereas the reverse crosses do not happen spontaneously but can be performed successfully with insemination . The F1 hybrids are morphological intermediates between the two species .
During the data analysis in a study concerning Danish brown hares we identified 16 hares with a very divergent haplotype. Here we show that these 16 individuals from the Danish wild brown hare populations carry mountain hare mtDNA haplotypes.
During alignment of 385 hare sequences, 16 sequences (6 males, 10 females) were clearly different from the rest and collapsed into a single haplotype. The 16 individuals were obtained from five of the eight populations sampled, and these five populations are all located on islands in the eastern part of Denmark (Figure 1). The frequency of sampled individuals with L. timidus mtDNA in those five populations averaged 6% (a crude estimate across all eight populations is 4.16%, Figure 1).
A BLAST search against GenBank with the divergent haplotype mtDNA returned L. timidus as the best hit (99% similarity), but not a single L. europaeus sequence in the first 50 hits. The subsequent phylogenetic analysis (Figure 2) revealed that the divergent Danish haplotype is grouped with L. timidus samples from across Europe (100% posterior probability in all three independent runs) and most closely related to the Swedish L. timidus haplotypes (100% posterior probability in all three independent runs, green clade in figure 2). The other Danish haplotypes included in the phylogenetic analysis, the 19 L. europaeus haplotypes, are grouped with the L. europaeus clade as expected, and they are most closely related to German haplotypes (100% posterior probability from three independent runs, red clade in figure 2). Hence, the Danish haplotypes clearly originate from two different species.
The results show the presence of mountain hare mtDNA in Danish brown hare populations. From mtDNA alone we cannot decide whether these 16 individuals are true mountain hares or introgressed brown hares. This would require further analysis using nuclear gene sequences (e.g. transferrin) (there are no differential diagnostic microsatellite loci for the two species ). However, given a) the frequency of the L. timidus haplotype where it is present (6% in the five populations), b) that L. timidus is not reported to be a native breeding species in Denmark , and c) that the two species are able to hybridize , these individuals most likely represent introgressed brown hares.
The most likely origins of introgression are; 1) natural migration of mountain hares from Sweden followed by introgression, or 2) natural migration or translocation of introgressed brown hares from Sweden. However, since the documented distribution of mountain hares in Sweden does not include southern Sweden [9, 10], natural migration of mountain hares is less probable than introduction of introgressed brown hares (carrying mountain hare mtDNA) from Sweden.
The fact that The Baltic Sea occasionally freezes over in very severe winters facilitates natural migration of introgressed brown hares from southern Sweden across the sea barrier to the Eastern Danish islands (Figure 1). This is also supported by the results from the genetic analyses (Figure 2).
During the mid 1980'ies a network of hare breeding facilities was established in Denmark, and it was initiated with imported hares primarily from Italy, France, Hungary and Sweden . No combined record is kept of the amount and origin of hares imported to Denmark. In 1993, when farming hares were prohibited in Denmark, there were 100 hare-farms in Denmark, with an annual export of up to 5,500 hares . As approximately 15% of Swedish L. europaeus specimens from sympatric areas carried transmitted L. timidus mtDNA , it is possible that phenotypic brown hare specimens imported from Sweden for breeding purposes have carried the L. timidus mtDNA and later escaped to the natural hare population.
The detection of only one haplotype in the 16 individuals shows low mtDNA diversity in the total population of hybrids in Denmark. The geographical distribution of the 16 intogressed brown hares (Figure 1) and the lack of haplotype diversity support a common origin in southern Sweden.
The implications of introgressed brown hares in the Danish hare populations are currently unclear. It has been suggested that hares with an alien mtDNA have a lowered fitness as a result of a functional incompatibility between the cytoplasmatic mitochondrial genes and the cell nucleus . In opposition to this, due to its observed high frequency in brown hares (93%) Melo-Ferreira et al.  suggested, that the ancient L. timidus mtDNA observed in the Iberian Peninsular might have some selective advantage depending on the nuclear background.
Thulin et al.  raised the question as to whether any 'pure' population of brown hares exists anywhere. Our findings give further support to this statement. This has important consequences for conservation and population genetics, e.g. problematic definition of management units, unclear fitness effects, genetically mixed populations and the inclusion of nuclear markers will be necessary for future hare studies.
Contrary to what is known, the Danish hare populations are not pure L. europaeus populations but include introgressed brown hares with Swedish L. timidus mtDNA. The most probable explanation is natural migration or translocation of introgressed brown hares from Sweden. The impurity of hare populations has implications for conservation and population genetics.
During 2003–2006 The National Environmental Research Institute of Denmark collected tissue samples from shot, wild hares at eight locations in Denmark (Figure 1). The hares were autopsied and several organs and muscle tissue samples were taken out for later analyses.
DNA was extracted from hare muscle tissue using a slight modification of the Chelex protocol . The tissue was placed into 1.5 mL Eppendorf tubes containing 200 μl of 20% Chelex resin solution (BIO-Rad, Hercules, CA, USA). Tubes were vortexed briefly, boiled at 100°C for 20 min in a heating block, centrifuged 3 min at 13000 rpm and stored at minus 20°C.
A 494 bp fragment of the mtDNA D-loop (control region) was amplified via PCR, using the mammalian control region primers L15997 5'-CACCATTAGCACCCAAAGCT-3' located in the tRNA gene and H16498 5'-CCTGAAGTAGGAACCAGATG-3' . PCR was performed in a total volume of 10 μl and contained 1 μl buffer (1.5 mM MgCl2), 1.6 μl dNTP (1.25 mM of A, C, G, T, respectively), 0.5 μl of each primer (10 pmol/μl), 0.1 μl Taq polymerase (Amersham Pharmacia Biotech), topped up with distilled water to 9 μl and 1 μl template DNA (40–60 ng/μl) was added. Cycling conditions were 94°C for 3 min, and 35 cycles of 94°C for 30 s, 50°C for 40 s and 72°C for 60 s and a final extension step of 72°C for 7 min. Sequencing was conducted under BigDye™ terminator cycling conditions, the reacted products purified using ethanol precipitation and run using an Automatic Sequencer ABI 3730xl. Both strands were sequenced in all samples.
Haplotype editing and collapsing
Nine L. europaeus and 20 L. timidus sequences were downloaded from GenBank. Sequences with geographical information were selected to represent haplotypes from all over Europe and Russia. Together with the 19 Danish L. europaeus and the divergent haplotype, these 49 sequences (table 1) were aligned using Muscle  and cropped to the shortest sequence (370 bp). Phylogenies and nodal support were estimated using MrBayes, version 3.0b4  under a Bayesian framework [27, 28] using a general time reversible substitution model (GTR + γ) and the coalescent branch length model. Bayesian analysis was initiated with random starting trees, run for 5 × 106 generations, and the Markov chain was sampled every 1000 generations. Model parameters were estimated directly from the data and three independent replicates were conducted to avoid entrapment in local optima . The initial 1,250 trees were discarded as "burn-in" and the remaining 3,750 trees were used to estimate nodal support as posterior probabilities.
Flux JEC, Angermann R: The hares and the jackrabbits. Rabbits, Hares and Pikas Status Survey and Conservation Action Plan. Edited by: Chapman JA, Flux JEC. 1990, Gland, Switzerland , IUCN, 61-94.
Smith RK, Vaughan Jennings N, Harris S: A quantitative analysis of the abundance and demography of European hares Lepus europaeus in relation to habitat type, intensity of agriculture and climate. Mammal Review. 2005, 35 (1): 1-24. 10.1111/j.1365-2907.2005.00057.x.
Schmidt NM, Asferg T, Forchhammer MC: Long-term patterns in European brown hare population dynamics in Denmark: effects of agriculture, predation and climate. BMC Ecol. 2004, 4: 15-10.1186/1472-6785-4-15.
The Danish Game Bag Record. http://vildtudbytte.dmu.dk/
Mitchell-Jones AJ: The atlas of European mammals. T & AD Poyser natural history. 1999, London , T & AD Poyser, xi, 484 p.-
Lönnberg E: On hybrids between Lepus timidus and Lepus europaeus (Pallas 1778). from southern Sweden. Proceedings of Zoological Society of London. 1905, 1: 278-287.
Thulin CG: The distribution of mountain hares Lepus timidus in Europe: a challenge from brown hares L. europaeus?. Mammal Review. 2003, 33 (1): 29-42. 10.1046/j.1365-2907.2003.00008.x.
Bergengren A: On genetics, evolution and history of distribution of the heath-hare, a distinct population of the Arctic Hare, Lepus timidus (Linneaus 1758). Swedish Wildlife. 1969, 6: 381-460.
Thulin CG, Jaarola M, Tegelstrom H: The occurrence of mountain hare mitochondrial DNA in wild brown hares. Molecular Ecology. 1997, 6 (5): 463-467. 10.1046/j.1365-294X.1997.t01-1-00199.x.
Thulin CG, Tegelstrom H: Biased geographical distribution of mitochondrial DNA that passed the species barrier from mountain hares to brown hares (genus Lepus): an effect of genetic incompatibility and mating behaviour?. Journal of Zoology, London,. 2002, 258: 299-306.
Thulin CG, Tegelström H, Fredga K: Haplotype diversity of mountain hare mtDNA among native mountain hares and introduced brown hares in Scandinavia. ANNALES ZOOLOGICI FENNICI. 2003, THE FINNISH ZOOLOGICAL PUBLISHING BOARD, 40 (1): 45-52.
Alves PC, Ferrand N, Suchentrunk F, Harris DJ: Ancient introgression of Lepus timidus mtDNA into L. granatensis and L. europaeus in the Iberian Peninsula. Mol Phylogenet Evol. 2003, 27 (1): 70-80. 10.1016/S1055-7903(02)00417-7.
Ben Slimen H, Suchentrunk F, Shahin AB, Ben Ammar Elgaaied A: Phylogenetic analysis of mtCR-1 sequences of Tunisian and Egyptian hares (Lepus sp. or spp., Lagomorpha) with different coat colours. Mammalian Biology - Zeitschrift fur Saugetierkunde. In Press, Corrected Proof:
Robinson TJ, Matthee CA: Phylogeny and evolutionary origins of the Leporidae: a review of cytogenetics, molecular analyses and a supermatrix analysis. Mammal Review. 2005, 35 (3-4): 231-247. 10.1111/j.1365-2907.2005.00073.x.
Gustavson I, Sundt CO: Anwendung von künstlicher Befruchtung bei der Hybridisierung von zwei Hasenarten. Zeitschrift für Jagdwissenschaft. 1965, 11: 155-158. 10.1007/BF01964745.
Notini G: Om harens biologi. Svenska Jägareförbundets Medelande. 1941, 4: 192-
Thulin CG, Stone J, Tegelström H, Walker CW: Species assignment and hybrid identification among Scandinavian hares Lepus europaeus and L. timidus. Wildlife Biology. 2006, 12: 29-38.
Hansen K, Jørgensen A, Henriksen P: Haren. i naturen, opdræt, sygdomme. 1990, Frederiksberg , Jordbrugsforlaget
Kleist MH: Haren. Natur og museum, 34 årgang, nr 1 (marts 1995). 1995, Århus , Naturhistorisk Museum
Melo-Ferreira J, Boursot P, Suchentrunk F, Ferrand N, Alves PC: Invasion from the cold past: extensive introgression of mountain hare (Lepus timidus) mitochondrial DNA into three other hare species in northern Iberia. Molecular Ecology. 2005, 14 (8): 2459-2464. 10.1111/j.1365-294X.2005.02599.x.
Walsh PS, Metzger DA, Higuchi R: Chelex 100 as a medium for simple extraction of DNA for PCR-based typing from forensic material. Biotechniques. 1991, 10 (4): 506-513.
Gerloff U, Hartung B, Fruth B, Hohmann G, Tautz D: Intracommunity relationships, dispersal pattern and paternity success in a wild living community of Bonobos (Pan paniscus) determined from DNA analysis of faecal samples. Proc Biol Sci. 1999, 266 (1424): 1189-1195. 10.1098/rspb.1999.0762.
Hall TA: BioEdit: a user-friendly biological sequence alignment editor and analysis program for Windows 95/98/NT. Nucleic Acids Symposium Series. 1999, 41: 95-98.
DNAcollapser, part of the fasta toolbox (FaBox). [http://www.daimi.au.dk/~biopv/php/fabox/]
Edgar RC: MUSCLE: multiple sequence alignment with high accuracy and high throughput. Nucleic Acids Res. 2004, 32 (5): 1792-1797. 10.1093/nar/gkh340.
Huelsenbeck JP, Ronquist F: MRBAYES: Bayesian inference of phylogenetic trees. Bioinformatics. 2001, 17 (8): 754-755. 10.1093/bioinformatics/17.8.754.
Lewis PO: Phylogenetic systematics turns over a new leaf. Trends In Ecology And Evolution. 2001, 16 (1): 30-37. 10.1016/S0169-5347(00)02025-5.
Shoemaker JS, Painter IS, Weir BS: Bayesian statistics in genetics: a guide for the uninitiated. Trends Genet. 1999, 15 (9): 354-358. 10.1016/S0168-9525(99)01751-5.
Huelsenbeck JP, Imennov NS: Geographic origin of human mitochondrial DNA: accommodating phylogenetic uncertainty and model comparison. Syst Biol. 2002, 51 (1): 155-165. 10.1080/106351502753475934.
We would like to thank the people involved in sample collection. This study was financed jointly by WWF Biodiversity Fund, Denmark (TF), the Augustinus Foundation, Denmark (TF) and National Environmental Research Insitute (NERI), Denmark. TF was further supported by a PhD grant from the Faculty of Science at the University of Aarhus. PV was supported by a grant from the Danish Cancer Society. TW was supported by a PhD grant from the National Environmental Research Institute, Denmark. We thank Cino Pertoldi and Liselotte W. Andersen for helpful comments on the manuscript.
TF carried out all the lab-work, edited and analyzed the sequence data and drafted the manuscript. PV carried out the phylogenetic analysis, prepared the figures and drafted the manuscript. TW organized the samples and drafted the manuscript. All authors have read and approved the final manuscript.