Email updates

Keep up to date with the latest news and content from BMC Ecology and BioMed Central.

Open Access Research article

Introgression of mountain hare (Lepus timidus) mitochondrial DNA into wild brown hares (Lepus europaeus) in Denmark

Tina Fredsted1*, Trine Wincentz2 and Palle Villesen3

Author Affiliations

1 Department of Ecology and Genetics, Institute of Biological Sciences, University of Aarhus, Ny Munkegade, building 1540, DK-8000 Aarhus C, Denmark

2 National Environmental Research Institute, Dept. of Wildlife Ecology and Biodiversity, Grenaavej 14, DK-8410 Ronde, Denmark

3 BiRC – Bioinformatics Research Center, University of Aarhus, H.-Guldbergs Gade 10, Building 1090, DK-8000 Aarhus C, Denmark

For all author emails, please log on.

BMC Ecology 2006, 6:17  doi:10.1186/1472-6785-6-17

The electronic version of this article is the complete one and can be found online at: http://www.biomedcentral.com/1472-6785/6/17


Received:12 June 2006
Accepted:15 November 2006
Published:15 November 2006

© 2006 Fredsted et al; licensee BioMed Central Ltd.

This is an Open Access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/2.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

Abstract

Background

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.

Results

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

Conclusion

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.

Background

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 [1]. 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 [2]. Additional factors like predators (especially foxes), birds of prey and traffic may have influenced the population locally [3]. 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 [4].

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 [5]. During the late nineteenth century the brown hare was introduced in southern Sweden as a game animal [6]. 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 [7]. 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 [8].

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 [12] and introgression of foreign mtDNA is likely to occur in several hare species from Asia [13]. 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 [14]. 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 [15]. The F1 hybrids are morphological intermediates between the two species [16].

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.

Results

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

thumbnailFigure 1. Location and frequency of the Lepus timidus haplotype in Denmark. Lepus timidus was found in five (large map) of the eight (small map) sampling locations.

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.

thumbnailFigure 2. Bayesian phylogeny of 370 bp from the mitochondrial D-loop of 20 Lepus timidus and 28 Lepus europaeus haplotypes from Europe, including the 19 Danish Lepus europaeus haplotypes, and the divergent Danish haplotype. The numbers show nodal support, given as posterior probabilities (for selected nodes, posterior probabilities from three independent runs are shown).

Discussion

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 [17]). 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 [5], and c) that the two species are able to hybridize [9], 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 [18]. 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 [19]. As approximately 15% of Swedish L. europaeus specimens from sympatric areas carried transmitted L. timidus mtDNA [10], 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 [10]. In opposition to this, due to its observed high frequency in brown hares (93%) Melo-Ferreira et al. [20] 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. [17] 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.

Conclusion

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.

Methods

Sample collection

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.

Sequencing

DNA was extracted from hare muscle tissue using a slight modification of the Chelex protocol [21]. 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' [22]. 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

Sequences were edited using Bioedit version 7.0.0 [23]. Identical haplotypes among the 385 sequences were found using DNAcollapser [24].

Phylogenetic analyses

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 [25] and cropped to the shortest sequence (370 bp). Phylogenies and nodal support were estimated using MrBayes, version 3.0b4 [26] 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 [29]. 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.

Table 1. Sampling localities and accession numbers of D-loop sequences used in the phylogenetic analysis (figure 2).

Authors' contributions

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.

Acknowledgements

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.

References

  1. Flux JEC, Angermann R: The hares and the jackrabbits. In Rabbits, Hares and Pikas Status Survey and Conservation Action Plan. Edited by Chapman JA, Flux JEC. Gland, Switzerland , IUCN; 1990:61-94. OpenURL

  2. 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. Publisher Full Text OpenURL

  3. 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. PubMed Abstract | BioMed Central Full Text | PubMed Central Full Text OpenURL

  4. The Danish Game Bag Record http://vildtudbytte.dmu.dk/ ;

  5. Mitchell-Jones AJ: The atlas of European mammals. In T & AD Poyser natural history. London , T & AD Poyser; 1999:xi, 484 p.. OpenURL

  6. 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. OpenURL

  7. 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. Publisher Full Text OpenURL

  8. 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. OpenURL

  9. Thulin CG, Jaarola M, Tegelstrom H: The occurrence of mountain hare mitochondrial DNA in wild brown hares.

    Molecular Ecology 1997, 6(5):463-467. PubMed Abstract | Publisher Full Text OpenURL

  10. 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. OpenURL

  11. Thulin CG, Tegelström H, Fredga K: Haplotype diversity of mountain hare mtDNA among native mountain hares and introduced brown hares in Scandinavia. In ANNALES ZOOLOGICI FENNICI. Volume 40. THE FINNISH ZOOLOGICAL PUBLISHING BOARD; 2003::45-52. OpenURL

  12. 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. PubMed Abstract | Publisher Full Text OpenURL

  13. 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 OpenURL

  14. 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. Publisher Full Text OpenURL

  15. Gustavson I, Sundt CO: Anwendung von künstlicher Befruchtung bei der Hybridisierung von zwei Hasenarten.

    Zeitschrift für Jagdwissenschaft 1965, 11:155-158. Publisher Full Text OpenURL

  16. Notini G: Om harens biologi.

    Svenska Jägareförbundets Medelande 1941, 4:192. OpenURL

  17. 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. OpenURL

  18. Hansen K, Jørgensen A, Henriksen P: Haren. i naturen, opdræt, sygdomme. Frederiksberg , Jordbrugsforlaget; 1990.

  19. Kleist MH: Haren. In Natur og museum, 34 årgang, nr 1 (marts 1995). Århus , Naturhistorisk Museum; 1995. OpenURL

  20. 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. PubMed Abstract | Publisher Full Text OpenURL

  21. 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. PubMed Abstract OpenURL

  22. 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. PubMed Abstract | Publisher Full Text OpenURL

  23. 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. OpenURL

  24. DNAcollapser, part of the fasta toolbox (FaBox) [http://www.daimi.au.dk/~biopv/php/fabox/] webcite

  25. Edgar RC: MUSCLE: multiple sequence alignment with high accuracy and high throughput.

    Nucleic Acids Res 2004, 32(5):1792-1797. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  26. Huelsenbeck JP, Ronquist F: MRBAYES: Bayesian inference of phylogenetic trees.

    Bioinformatics 2001, 17(8):754-755. PubMed Abstract | Publisher Full Text OpenURL

  27. Lewis PO: Phylogenetic systematics turns over a new leaf.

    Trends In Ecology And Evolution 2001, 16(1):30-37. Publisher Full Text OpenURL

  28. Shoemaker JS, Painter IS, Weir BS: Bayesian statistics in genetics: a guide for the uninitiated.

    Trends Genet 1999, 15(9):354-358. PubMed Abstract | Publisher Full Text OpenURL

  29. Huelsenbeck JP, Imennov NS: Geographic origin of human mitochondrial DNA: accommodating phylogenetic uncertainty and model comparison.

    Syst Biol 2002, 51(1):155-165. PubMed Abstract | Publisher Full Text OpenURL