Hybridization and introgression have played important roles in the history of various species, including lineage diversification and the evolution of adaptive traits. Hybridization can accelerate the development of reproductive isolation between diverging species, and thus valuable insight into the evolution of reproductive barrier formation may be gained by studying secondary contact zones. Hedgehogs of the genus Erinaceus, which are insectivores sensitive to changes in climate, are a pioneer model in Pleistocene phylogeography. The present study provides the first genome-wide SNP data regarding the Erinaceus hedgehogs species complex, offering a unique comparison of two secondary contact zones between Erinaceus europaeus and E. roumanicus. Results confirmed diversification of the genus during the Pleistocene period, and detected a new refugial lineage of E. roumanicus outside the Mediterranean basin, most likely in the Ponto-Caspian region. In the Central European zone, the level of hybridization was low, whereas in the Russian-Baltic zone, both species hybridise extensively. Asymmetrical gene flow from E. europaeus to E. roumanicus suggests that reproductive isolation varies according to the direction of the crosses in the hybrid zones. However, no loci with significantly different patterns of introgression were detected. Markedly different pre- and post-zygotic barriers, and thus diverse modes of species boundary maintenance in the two contact zones, likely exist. This pattern is probably a consequence of the different age and thus of the different stage of evolution of reproductive isolating mechanisms in each hybrid zone.

Center for Species Survival Smithsonian Conservation Biology Institute National Zoological Park Washington DC USA

Department of Animal Science and Food Processing Faculty of Tropical AgriSciences Czech University of Life Sciences Prague Prague Czech Republic

Department of Zoology Faculty of Science Charles University Prague Czech Republic

EMBL Genomics Core Facility Heidelberg Germany

Evolutionary Adaptive Genomics University of Potsdam Potsdam Germany

Institut Cavanilles de Biodiversitat i Biologia Evolutiva Valencia Spain

Lomonosov Moscow State University Moscow Russia

Zobrazit více v PubMed

Ai H, Fang X, Yang B, Huang Z, Chen H, Mao L, et al. Adaptation and possible ancient interspecies introgression in pigs identified by whole-genome sequencing. Nat Genet. 2015;47:217–225. PubMed

Alexander DH, Lange K. Enhancements to the ADMIXTURE algorithm for individual ancestry estimation. BMC Bioinforma. 2011;12:246. PubMed PMC

Alexander DH, Novembre J, Lange K. Fast model-based estimation of ancestry in unrelated individuals. Genome Res. 2009;19:1655–1664. PubMed PMC

Allen JRM, Hickler T, Singarayer JS, Sykes MT, Valdes PJ, Huntley B. Last glacial vegetation of Northern Eurasia. Quat Sci Rev. 2010;29:2604–2618.

Amori G (2016) Erinaceus europaeus. The IUCN Red List of Threatened Species: e.T29650A2791303

Amori G, Hutterer R, Kryštufek B, Yigit N, Mitsain G & Palomo LJ (2016a) Erinaceus concolor (errata version published in 2017) The IUCN Red List of Threatened Species: e.T40605A115174497

Amori G, Hutterer R, Kryštufek B, Yigit, N, Mitsain G & Palomo LJ (2016b) Erinaceus roumanicus (errata version published in 2017) The IUCN Red List of Threatened Species: e.T136344A115206348.

Anderson EC, Thompson EA. A model-based method for identifying species hybrids using multilocus genetic data. Genetics. 2002;160:1217–1229. PubMed PMC

Andrade P, Pinho C, i de Lanuza GP, Afonso S, Brejcha J, Rubin CJ, et al. Regulatory changes in pterin and carotenoid genes underlie balanced color polymorphisms in the wall lizard. Proc Natl Acad Sci USA. 2019;116:5633–5642. PubMed PMC

Babik W, Branicki W, Sandera M, Litvinchuk S, Borkin LJ, Irwin JT, et al. Mitochondrial phylogeography of the moor frog, Rana arvalis. Mol Ecol. 2004;13:1469–1480. PubMed

Baird NA, Etter PD, Atwood TS, Currey MC, Shiver AL, Lewis ZA, et al. Rapid SNP discovery and genetic mapping using sequenced RAD markers. PLoS One. 2008;3:3376. PubMed PMC

Bannikova AA, Lebedev VS, Abramov AV, Rozhnov VV. Contrasting evolutionary history of hedgehogs and gymnures (Mammalia: Erinaceomorpha) as inferred from a multigene study. Biol J Linn Soc. 2014;112:499–519.

Bannikova AA, Fedorova LV, Fedorov AN, Troitskiĭ AV, Grechko VV, Dolgov VA, et al. Comparison of repeat DNA sequences in the Erinaceidae mammal family by restriction analysis. Genetika. 1995;31:1498–1506. PubMed

Barton NH, Hewitt GM. Analysis of Hybrid Zones. Annu Rev Ecol Syst. 1985;16:113–148.

Bierne N, Welch J, Loire E, Bonhomme F, David P. The coupling hypothesis: why genome scans may fail to map local adaptation genes. Mol Ecol. 2011;20:2044–2072. PubMed

Bogdanov AS, Bannikova AA, Pirusskii YM, Formozov NA. The first genetic evidence of hybridization between West European and Northern white-breasted hedgehogs (Erinaceus europaeus and E. roumanicus) in Moscow region. Biol Bull. 2009;36:647–651. PubMed

Bouckaert R, Heled J, Kühnert D, Vaughan T, Wu CH, Xie D, et al. BEAST 2: A Software Platform for Bayesian Evolutionary Analysis. PLoS Comput Biol. 2014;10:e1003537. PubMed PMC

Cassola F (2016) Erinaceus amurensis (errata version published in 2017) The IUCN Red List of Threatened Species: e.T40604A115174360

Catchen J, Bassham S, Wilson T, Currey M, O’Brien C, Yeates Q, et al. The population structure and recent colonization history of Oregon threespine stickleback determined using restriction-site associated DNA-sequencing. Mol Ecol. 2013;22:2864–2883. PubMed PMC

Černá Bolfíková B, Eliášová K, Loudová M, Kryštufek B, Lymberakis P, Sándor AD, et al. Glacial allopatry vs. postglacial parapatry and peripatry: The case of hedgehogs. PeerJ. 2017;5:e3163. PubMed PMC

Cooper SJB, Ibrahim KM, Hewitt GM. Postglacial expansion and genome subdivision in the European grasshopper Chorthippus parallelus. Mol Ecol. 1995;4:49–60. PubMed

Corbet GB. The family Erinaceidae: a synthesis of its taxonomy, phylogeny, ecology and zoogeography. Mamm Rev. 1988;18:117–172.

Cordova CE, Harrison SP, Mudie PJ, Riehl S, Leroy SAG, Ortiz N. Pollen, plant macrofossil and charcoal records for palaeovegetation reconstruction in the Mediterranean-Black Sea Corridor since the Last Glacial Maximum. Quat Int. 2009;197:12–26.

Coyne J, Orr H. Speciation. Sunderland: Sinauer Associates; 2004.

Curto M, Winter S, Seiter A, Schmid L, Scheicher K, Barthel LMF, et al. Application of a SSR-GBS marker system on investigation of European hedgehog species and their hybrid zone dynamics. Ecol Evol. 2019;9:2814–2832. PubMed PMC

Danecek P, Auton A, Abecasis G, Albers CA, Banks E, DePristo MA, et al. The variant call format and VCFtools. Bioinformatics. 2011;27:2156–2158. PubMed PMC

Dobzhansky T. Genetics and the origin of species. New York: Columbia University Press; 1937.

Durand EY, Patterson N, Reich D, Slatkin M. Testing for ancient admixture between closely related populations. Mol Biol Evol. 2011;28:2239–2252. PubMed PMC

Filippucci MG, Lapini L. First data on the genetic differentiation between Erinaceus europaeus Linnaeus, 1758 and Erinaceus concolor Martin, 1838 in North-eastern Italy (Mammalia, Insectivora, Erinaceidae) Gortania. 1988;9:227–236.

Fitak RR (2019) optM: an R package to optimize the number of migration edges using threshold models. v0.1.3. J Hered. https://cran.r-project.org/web/packages/OptM/index.html

García-Vázquez A, Pinto Llona AC, Grandal-d’Anglade A. Post-glacial colonization of Western Europe brown bears from a cryptic Atlantic refugium out of the Iberian Peninsula. Hist Biol. 2019;31:618–630.

Garrison E, Marth G (2012) Haplotype-based variant detection from short-read sequencing. arXiv:1207.3907. http://arxiv.org/abs/1207.3907

Gompert Z, Buerkle CA. Bayesian estimation of genomic clines. Mol Ecol. 2011;20:2111–2127. PubMed

Gompert Z, Buerkle CA. Bgc: Software for Bayesian estimation of genomic clines. Mol Ecol Resour. 2012;12:1168–1176. PubMed

Green RE, Krause J, Briggs AW, Maricic T, Stenzel U, Kircher M, et al. A draft sequence of the neandertal genome. Science. 2010;328:710–722. PubMed PMC

Harrison RG, Larson EL. Hybridization, introgression, and the nature of species boundaries. J Hered. 2014;105:795–809. PubMed

He K, Chen J-H, Gould GC, Yamaguchi N, Ai H-S, Wang Y-X, et al. An estimation of erinaceidae phylogeny: a combined analysis approach. PLoS One. 2012;7:e39304. PubMed PMC

Hewitt GM. Hybrid zones-natural laboratories for evolutionary studies. Trends Ecol Evol. 1988;3:158–167. PubMed

Hewitt GM. Post-glacial re-colonization of European biota. Biol J Linn Soc. 1999;68:87–112.

Jackson DB. The breeding biology of introduced hedgehogs (Erinaceus europaeus) on a Scottish Island: Lessons for population control and bird conservation. J Zool. 2006;268:303–314.

Johannesson K, Le Moan A, Perini S, André C. A Darwinian Laboratory of Multiple Contact Zones. Trends Ecol Evol. 2020;35:1021–1036. PubMed

Jombart T, Ahmed I. Adegenet 1.3-1: New tools for the analysis of genome-wide SNP data. Bioinformatics. 2011;27:3070–3071. PubMed PMC

Korneliussen TS, Albrechtsen A, Nielsen R. ANGSD: analysis of next generation sequencing data. BMC Bioinforma. 2014;15:356. doi: 10.1186/s12859-014-0356-4. PubMed DOI PMC

Kotlík P, Deffontaine V, Mascheretti S, Zima J, Michaux JR, Searle JB. A northern glacial refugium for bank voles (Clethrionomys glareolus) Proc Natl Acad Sci USA. 2006;103:14860–14864. PubMed PMC

Kratochvil J. Zur Frage der Verbreitung des Igels (Erinaceus) in der ČSSR. Zool List. 1966;15:291–304.

Kratochvíl J. Zur Kenntnis der Igel der Gattung Erinaceus in der ČSSR (Insectivora, Mamm.) Zool List. 1975;24:297–312.

Kumar V, Lammers F, Bidon T, Pfenninger M, Kolter L, Nilsson MA, et al. The evolutionary history of bears is characterized by gene flow across species. Sci Rep. 2017;7:46487. PubMed PMC

Langmead B, Trapnell C, Pop M, Salzberg SL. Ultrafast and memory-efficient alignment of short DNA sequences to the human genome. Genome Biol. 2009;10:R25. PubMed PMC

Li G, Davis BW, Eizirik E, Murphy WJ. Phylogenomic evidence for ancient hybridization in the genomes of living cats (Felidae) Genome Res. 2016;26:1–11. PubMed PMC

Li H, Handsaker B, Wysoker A, Fennell T, Ruan J, Homer N, et al. The Sequence Alignment/Map format and SAMtools. Bioinformatics. 2009;25:2078–2079. PubMed PMC

Linck E, Battey CJ. Minor allele frequency thresholds strongly affect population structure inference with genomic data sets. Mol Ecol Resour. 2019;19:639–647. PubMed

Lischer HEL, Excoffier L. PGDSpider: an automated data conversion tool for connecting population genetics and genomics programs. Bioinformatics. 2012;28:298–299. PubMed

Liu S, Lorenzen ED, Fumagalli M, Li B, Harris K, Xiong Z, et al. Population genomics reveal recent speciation and rapid evolutionary adaptation in polar bears. Cell. 2014;157:785–794. PubMed PMC

Malinsky M, Challis RJ, Tyers AM, Schiffels S, Terai Y, Ngatunga BP, et al. Genomic islands of speciation separate cichlid ecomorphs in an East African crater lake. Science. 2015;350:1493–1498. PubMed PMC

Mallet J. Hybridization as an invasion of the genome. Trends Ecol Evol. 2005;20:229–237. PubMed

Martin SH, Davey JW, Jiggins C. Evaluating the Use of ABBA-BABA Statistics to Locate Introgressed Loci. Mol Biol Evol. 2015;32:244–257. PubMed PMC

Martin SH, Dasmahapatra KK, Nadeau NJ, Salazar C, Walters JR, Simpson F, et al. Genome-wide evidence for speciation with gene flow in Heliconius butterflies. Genome Res. 2013;23:1817–1828. PubMed PMC

Nadeau NJ, Ruiz M, Salazar P, Counterman B, Medina JA, Ortiz-Zuazaga H, et al. Population genomics of parallel hybrid zones in the mimetic butterflies, H. melpomene and H. erato. Genome Res. 2014;24:1316–1333. PubMed PMC

Patterson N, Moorjani P, Luo Y, Mallick S, Rohland N, Zhan Y, et al. Ancient admixture in human history. Genetics. 2012;192:1065–1093. PubMed PMC

Payseur BA, Rieseberg LH. A genomic perspective on hybridization and speciation. Mol Ecol. 2016;25:2337–2360. PubMed PMC

Pickrell JK, Pritchard JK. Inference of population splits and mixtures from genome-wide allele frequency data. PLoS Genet. 2012;8:e1002967. PubMed PMC

Poduschka W, Poduschka C. Kreuzungsversuche an mitteleuropäischen Igeln. Säugetierkundliche Mitteilungen. 1983;31:1–12.

Posada D. jModelTest: phylogenetic model averaging. Mol Biol Evol. 2008;25:1253–1256. PubMed

Purcell S, Neale B, Todd-Brown K, Thomas L, Ferreira MAR, Bender D, et al. PLINK: A tool set for whole-genome association and population-based linkage analyses. Am J Hum Genet. 2007;81:559–575. PubMed PMC

R Core Team (2013) R: A language and environment for statistical computing. R Foundation for Statistical Computing, Vienna, Austria. http://www.R-project.org/

Rambaut A (2018) FigTree v1.4.4. http://tree.bio.ed.ac.uk/software-/figtree/

Rambaut A, Drummond A (2007) Tracer v1.4. http://tree.bio.ed.ac.uk/software/tracer/

Reeve NJ, Morris PA. Mating strategy in the hedgehog (Erinaceus europaeus) J Zool. 1986;210:613–614.

Reeve N (1994) Hedgehogs. T. & A.D. Poyser: London

Reich D, Green RE, Kircher M, Krause J, Patterson N, Durand EY, et al. Genetic history of an archaic hominin group from Denisova cave in Siberia. Nature. 2010;468:1053–1060. PubMed PMC

Remington CL. Suture-Zones of Hybrid Interaction Between Recently Joined Biotas. In: Dobzhansky T, Hecht MK, Steere CW, editors. Evolutionary Biology. New York: Plenum Press; 1968. pp. 321–428.

Rockman MV. The QTN program and the alleles that matter for evolution: All that’s gold does not glitter. Evolution (NY) 2012;66:1–17. PubMed PMC

Ronquist F, Huelsenbeck JP. MrBayes 3: Bayesian phylogenetic inference under mixed models. Bioinformatics. 2003;19:1572–1574. PubMed

roumanicus in Central Europe. Heredity 108:248–255 PubMed PMC

Santucci F, Emerson BC, Hewitt GM. Mitochondrial DNA phylogeography of European hedgehogs. Mol Ecol. 1998;7:1163–1172. PubMed

Schmitt T, Varga Z. Extra-Mediterranean refugia: The rule and not the exception? Front Zool. 2012;9:22. PubMed PMC

Schwarz G. Estimating the dimension of a model. Ann Stat. 1978;6:461–464.

Seddon JM, Santucci F, Reeve NJ, Hewitt GM. DNA footprints of European hedgehogs, Erinaceus europaeus and E. concolor: Pleistocene refugia, postglacial expansion and colonization routes. Mol Ecol. 2001;10:2187–2198. PubMed

Smadja CM, Butlin RK. A framework for comparing processes of speciation in the presence of gene flow. Mol Ecol. 2011;20:5123–5140. PubMed

Sommer RS. When east met west: the sub-fossil footprints of the west European hedgehog and the northern white-breasted hedgehog during the Late Quaternary in Europe. J Zool. 2007;273:82–89.

Stephens PA, Sutherland WJ. Consequences of the Allee effect for behaviour, ecology and conservation. Trends Ecol Evol. 1999;14:401–405. PubMed

Sukumaran J, Holder MT. DendroPy: A Python library for phylogenetic computing. Bioinformatics. 2010;26:1569–1571. PubMed

Szarmach SJ, Brelsford A, Witt CC, Toews DPL. Comparing divergence landscapes from reduced-representation and whole genome resequencing in the yellow-rumped warbler (Setophaga coronata) species complex. Mol Ecol. 2021;30:5994–6005. PubMed

Taberlet P, Bouvet J. Mitochondrial DNA polymorphism, phylogeography, and conservation genetics of the brown bear Ursus arctos in Europe. Proc R Soc Lond Ser B Biol Sci. 1994;255:195–200. PubMed

Taylor SA, Larson EL. Insights from genomes into the evolutionary importance and prevalence of hybridization in nature. Nat Ecol Evol. 2019;3:170–177. PubMed

Varga Z. Extra-mediterranean refugia, post-glacial vegetation history and area dynamics in Eastern Central Europe. In: Habel J, Assmann T, editors. Relict Species: Phylogeography and Conservation Biology. Berlin, Heidelberg: Springer-Verlag; 2010. pp. 57–87.

Vijay N, Bossu CM, Poelstra JW, Weissensteiner MH, Suh A, Kryukov AP, et al. Evolution of heterogeneous genome differentiation across multiple contact zones in a crow species complex. Nat Commun. 2016;7:13195–13205. PubMed PMC

Weir BS, Cockerham CC. Estimating F-statistics for the analysis of population structure. Evolution (NY) 1984;38:1358. PubMed

Wu DD, Ding XD, Wang S, Wójcik JM, Zhang Y, Tokarska M, et al. Pervasive introgression facilitated domestication and adaptation in the Bos species complex. Nat Ecol Evol. 2018;2:1139–1145. PubMed

Zheng Y, Janke A. Gene flow analysis method, the D-statistic, is robust in a wide parameter space. BMC Bioinforma. 2018;19:10. PubMed PMC

Zieliński P, Dudek K, Arntzen JW, Palomar G, Niedzicka M, Fijarczyk A, et al. Differential introgression across newt hybrid zones: evidence from replicated transects. Mol Ecol. 2019;28:4811–4824. PubMed

Zolotareva KI, Belokon MM, Belokon YS, Rutovskaya MV, Hlyap LA, Starykov VP, et al. Genetic diversity and structure of the hedgehogs Erinaceus europaeus and Erinaceus roumanicus: evidence for ongoing hybridization in Eastern Europe. Biol J Linn Soc. 2021;132:174–195.