Deciduous forests form the dominant natural vegetation of Europe today, but were restricted to small refugia during Pleistocene cold stages, implying an evolutionary past shaped by recurrent range contractions and expansions. Cold-stage forest refugia were probably widespread in southern and central Europe, with the northwestern Balkan Peninsula being of particular importance. However, the actual number and location of deciduous forest refugia, as well as the connections between them, remain disputed. Here, we address the evolutionary dynamics of the deciduous forest understorey species Euphorbia carniolica as a proxy for past forest dynamics. To do so, we obtained genomic and morphometric data from populations representing the species' entire range, investigated phylogenetic position and intraspecific genetic variation, tested explicit demographic scenarios and applied species distribution models. Our data support two disjoint groups linked to separate refugia on the northwestern and central Balkan Peninsula. We find that genetic differentiation between groups started in the early Pleistocene via vicariance, suggesting a larger distribution in the past. Both refugia acted as sources for founder events to the southeastern Alps and the Carpathians; the latter were likely colonised before the last cold stage. In line with traditional views on the pre-Pleistocene origin of many southeastern European deciduous forest species, the origin of E. carniolica was dated to the late Pliocene. The fact that E. carniolica evolved at a time when a period of continuous forestation was ending in much of Eurasia provides an interesting biogeographical perspective on the past links between Eurasian deciduous forests and their biota.

Department of Botany and Biodiversity Research University of Vienna Vienna Austria

Department of Botany University of Innsbruck Innsbruck Austria

Faculty of Agricultural Environmental and Food Sciences Free University of Bozen Bolzano Bolzano Italy

Institute of Botany of the Czech Academy of Sciences Průhonice Czechia

Zobrazit více v PubMed

Abràmoff, M. , Magalhães, P. , & Ram, S. J. (2004). Image processing with ImageJ. Biophotonics International, 11, 36–42.

Asch, K. (2005). IGME 5000: 1: 5 million international geological map of Europe and adjacent areas [map]. BGR.

Bálint, M. , Ujvárosi, L. , Theissinger, K. , Lehrian, S. , Mészáros, N. , & Pauls, S. U. (2011). The Carpathians as a major diversity hotspot in Europe. In Zachos F. E. & Habel J. C. (Eds.), Biodiversity hotspots: Distribution and protection of conservation priority areas (pp. 189–205). Springer. 10.1007/978-3-642-20992-5_11 DOI

Bartha, L. , Sramkó, G. , Volkova, P. A. , Surina, B. , Ivanov, A. L. , & Banciu, H. L. (2015). Patterns of plastid DNA differentiation in DOI

Broennimann, O. , Fitzpatrick, M. C. , Pearman, P. B. , Petitpierre, B. , Pellissier, L. , Yoccoz, N. G. , Thuiller, W. , Fortin, M.‐J. , Randin, C. , Zimmermann, N. E. , Graham, C. H. , & Guisan, A. (2012). Measuring ecological niche overlap from occurrence and spatial environmental data. Global Ecology and Biogeography, 21(4), 481–497. 10.1111/j.1466-8238.2011.00698.x DOI

Burnham, K. P. , & Anderson, D. R. (Eds.). (2002). Model selection and multimodel inference (2nd ed.). Springer New York. 10.1007/b97636 DOI

Caković, D. , Cresti, L. , Stešević, D. , Schönswetter, P. , & Frajman, B. (2021). High genetic and morphological diversification of the DOI

Catchen, J. , Hohenlohe, P. , Bassham, S. , Amores, A. , & Cresko, W. (2013). Stacks: An analysis tool set for population genomics. Molecular Ecology, 22(11), 3124–3140. 10.1111/mec.12354 PubMed DOI PMC

Charles, K. L. , Bell, R. C. , Blackburn, D. C. , Burger, M. , Fujita, M. K. , Gvoždík, V. , Jongsma, G. F. M. , Kouete, M. T. , Leaché, A. D. , & Portik, D. M. (2018). Sky, sea, and forest islands: Diversification in the African leaf‐folding frog

Clement, M. , Posada, D. , & Crandall, K. A. (2000). TCS: A computer program to estimate gene genealogies. Molecular Ecology, 9(10), 1657–1659. 10.1046/j.1365-294x.2000.01020.x PubMed DOI

Cotto, O. , Wessely, J. , Georges, D. , Klonner, G. , Schmid, M. , Dullinger, S. , Thuiller, W. , & Guillaume, F. (2017). A dynamic eco‐evolutionary model predicts slow response of alpine plants to climate warming. Nature Communications, 8, 15399. 10.1038/ncomms15399 PubMed DOI PMC

Danecek, P. , Auton, A. , Abecasis, G. , Albers, C. A. , Banks, E. , DePristo, M. A. , Handsaker, R. E. , Lunter, G. , Marth, G. T. , Sherry, S. T. , McVean, G. , & Durbin, R. (2011). The variant call format and VCFtools. Bioinformatics, 27(15), 2156–2158. 10.1093/bioinformatics/btr330 PubMed DOI PMC

Drummond, A. J. , Suchard, M. A. , Xie, D. , & Rambaut, A. (2012). Bayesian phylogenetics with BEAUti and the beast 1.7. Molecular Biology and Evolution, 29(8), 1969–1973. 10.1093/molbev/mss075 PubMed DOI PMC

Đurović, S. Z. , Temunović, M. , Niketić, M. , Tomović, G. , Schönswetter, P. , & Frajman, B. (2021). Impact of quaternary climatic oscillations on phylogeographic patterns of three habitat‐segregated DOI

Dvořáková, H. , Fér, T. , & Marhold, K. (2010). Phylogeographic pattern of the European forest grass species DOI

Ehlers, J. , Gibbard, P. L. , & Hughes, P. D. (2011). Quaternary glaciations – Extent and chronology: A closer look (Vol. 15). Elsevier.

Ellenberg, H. , & Leuschner, C. (2010). Vegetation Mitteleuropas mit den Alpen in ökologischer, dynamischer und historischer Sicht (6th ed.). UTB.

Evanno, G. , Regnaut, S. , & Goudet, J. (2005). Detecting the number of clusters of individuals using the software STRUCTURE: A simulation study. Molecular Ecology, 14, 2611–2620. 10.1111/j.1365-294X.2005.02553.x PubMed DOI

Frajman, B. , Graniszewska, M. , & Schönswetter, P. (2016). Evolutionary patterns and morphological diversification within the European members of the

Frajman, B. , & Schönswetter, P. (2011). Giants and dwarfs: Molecular phylogenies reveal multiple origins of annual spurges within PubMed DOI

Gernhard, T. (2008). The conditioned reconstructed process. Journal of Theoretical Biology, 253(4), 769–778. 10.1016/j.jtbi.2008.04.005 PubMed DOI

Gömöry, D. , Paule, L. , Shvadchak, I. , Flaviu, P. , Sulkowska, M. , Hynek, V. , & Longauer, R. (2003). Spatial patterns of the genetic differentiation in European beech (

Govaerts, R. , Frodin, D. G. , & Radcliffe‐Smith, A. (2000). World checklist and bibliography of Euphorbiaceae (and Pandaceae) (Vol. 2). Royal Botanic Gardens, Kew.

Gubler, M. , Henne, P. D. , Schwörer, C. , Boltshauser‐Kaltenrieder, P. , Lotter, A. F. , Brönnimann, S. , & Tinner, W. (2018). Microclimatic gradients provide evidence for a glacial refugium for temperate trees in a sheltered hilly landscape of northern Italy. Journal of Biogeography, 45(11), 2564–2575. 10.1111/jbi.13426 DOI

Gutenkunst, R. N. , Hernandez, R. D. , Williamson, S. H. , & Bustamante, C. D. (2009). Inferring the joint demographic history of multiple populations from multidimensional SNP frequency data. PLoS Genetics, 5(10), e1000695. PubMed PMC

Hewitt, G. M. (1999). Post‐glacial re‐colonization of European biota. Biological Journal of the Linnean Society, 68(1–2), 87–112. 10.1111/j.1095-8312.1999.tb01160.x DOI

Horn, J. W. , Xi, Z. , Riina, R. , Peirson, J. A. , Yang, Y. , Dorsey, B. L. , Berry, P. E. , Davis, C. C. , & Wurdack, K. J. (2014). Evolutionary bursts in PubMed DOI

Horvat, I. , Glavač, I. , & Ellenberg, H. (1974). Vegetation Südosteuropas. Gustav Fischer.

Huang, X.‐C. , German, D. A. , & Koch, M. A. (2020). Temporal patterns of diversification in Brassicaceae demonstrate decoupling of rate shifts and mesopolyploidization events. Annals of Botany, 125(1), 29–47. 10.1093/aob/mcz123 PubMed DOI PMC

Hurka, H. , Friesen, N. , Bernhardt, K.‐G. , Neuffer, B. , Smirnov, S. , Shmakov, A. , & Blattner, F. (2019). The Eurasian steppe belt: Status quo, origin and evolutionary history. Turczaninowia, 22(3), 5–71. 10.14258/turczaninowia.22.3.1 DOI

Janes, J. K. , Miller, J. M. , Dupuis, J. R. , Malenfant, R. M. , Gorrell, J. C. , Cullingham, C. I. , & Andrew, R. L. (2017). The K = 2 conundrum. Molecular Ecology, 26(14), 3594–3602. 10.1111/mec.14187 PubMed DOI

Jiménez‐Alfaro, B. , Girardello, M. , Chytrý, M. , Svenning, J.‐C. , Willner, W. , Gégout, J.‐C. , Agrillo, E. , Campos, J. A. , Jandt, U. , Kącki, Z. , Šilc, U. , Slezák, M. , Tichý, L. , Tsiripidis, I. , Turtureanu, P. D. , Ujházyová, M. , & Wohlgemuth, T. (2018). History and environment shape species pools and community diversity in European beech forests. Nature Ecology & Evolution, 2, 483–490. 10.1038/s41559-017-0462-6 PubMed DOI

Jombart, T. , & Ahmed, I. (2011). Adegenet 1.3‐1: New tools for the analysis of genome‐wide SNP data. Bioinformatics, 27(21), 3070–3071. 10.1093/bioinformatics/btr521 PubMed DOI PMC

Kaltenrieder, P. , Belis, C. A. , Hofstetter, S. , Ammann, B. , Ravazzi, C. , & Tinner, W. (2009). Environmental and climatic conditions at a potential glacial refugial site of tree species near the southern alpine glaciers. New insights from multiproxy sedimentary studies at Lago della Costa (Euganean Hills, northeastern Italy). Quaternary Science Reviews, 28(25), 2647–2662. 10.1016/j.quascirev.2009.05.025 DOI

Kearse, M. , Moir, R. , Wilson, A. , Stones‐Havas, S. , Cheung, M. , Sturrock, S. , Buxton, S. , Cooper, A. , Markowitz, S. , Duran, C. , Thierer, T. , Ashton, B. , Meintjes, P. , & Drummond, A. (2012). Geneious basic: An integrated and extendable desktop software platform for the organization and analysis of sequence data. Bioinformatics, 28(12), 1647–1649. 10.1093/bioinformatics/bts199 PubMed DOI PMC

Leroy, S. A. G. (2007). Progress in palynology of the Gelasian–Calabrian Stages in Europe: Ten messages. Revue de Micropaleontologie, 50(4), 293–308. 10.1016/j.revmic.2006.08.001 DOI

Leroy, S. A. G. , Arpe, K. , & Mikolajewicz, U. (2011). Vegetation context and climatic limits of the early Pleistocene hominin dispersal in Europe. Quaternary Science Reviews, 30(11), 1448–1463. 10.1016/j.quascirev.2010.01.017 DOI

Lewis, P. O. (2001). A likelihood approach to estimating phylogeny from discrete morphological character data. Systematic Biology, 50(6), 913–925. 10.1080/106351501753462876 PubMed DOI

Magri, D. , Vendramin, G. G. , Comps, B. , Dupanloup, I. , Geburek, T. , Gömöry, D. , Latałowa, M. , Litt, T. , Paule, L. , Roure, J. M. , Tantau, I. , van der Knaap, W. O. , Petit, R. J. , & de Beaulieu, J.‐L. (2006). A new scenario for the quaternary history of European beech populations: Palaeobotanical evidence and genetic consequences. New Phytologist, 171(1), 199–221. 10.1111/j.1469-8137.2006.01740.x PubMed DOI

Magyari, E. (2002). Holocene biogeography of

Marinček, L. (1995). Contribution to demarcation and phytogeographic division of the Illyrian floral province, based on vegetation and flora. Gortania, 16, 99–124.

Meusel, H. , Jäger, E. , & Rauschert, S. (1978). Vergleichende Chorologie der zentraleuropäischen Flora II: Vol. II. Gustav Fischer Verlag.

Milne, R. I. , & Abbott, R. J. (2002). The origin and evolution of tertiary relict floras. In Advances in botanical research (Vol. 38, pp. 281–314). Academic Press. 10.1016/S0065-2296(02)38033-9 DOI

Mitka, J. , Baba, W. , & Szczepanek, K. (2014). Putative forest glacial refugia in the Western and eastern Carpathians. Modern Phytomorphology, 5, 85–92. 10.5281/zenodo.161009 DOI

Mráz, P. , & Ronikier, M. (2016). Biogeography of the Carpathians: Evolutionary and spatial facets of biodiversity. Biological Journal of the Linnean Society, 119(3), 528–559. 10.1111/bij.12918 DOI

Müller, K. (2005). SeqState: Primer design and sequence statistics for phylogenetic DNA datasets. Applied Bioinformatics, 4(1), 65–69. 10.2165/00822942-200504010-00008 PubMed DOI

Nikolić, T. (2021). Rasprostranjenost

Nylander, J. A. A. (2004). MrAIC.Pl. program distributed by the author. Evolutionary Biology Centre, Uppsala University.

Pahlevani, A. H. , & Frajman, B. (2023). Widespread, but less than assumed: Populations of DOI

Paris, J. R. , Stevens, J. R. , & Catchen, J. M. (2017). Lost in parameter space: A road map for stacks. Methods in Ecology and Evolution, 8(10), 1360–1373. 10.1111/2041-210X.12775 DOI

Parks, C. R. , & Wendel, J. F. (1990). Molecular divergence between Asian and north American species of DOI

Pattengale, N. D. , Alipour, M. , Bininda‐Emonds, O. R. P. , Moret, B. M. E. , & Stamatakis, A. (2010). How many bootstrap replicates are necessary? Journal of Computational Biology, 17(3), 337–354. 10.1089/cmb.2009.0179 PubMed DOI

Paun, O. , Turner, B. , Trucchi, E. , Munzinger, J. , Chase, M. W. , & Samuel, R. (2016). Processes driving the adaptive radiation of a tropical tree ( PubMed DOI PMC

Pignatti, S. , Guarino, R. , & la Rosa, M. (2019). Flora d'Italia (Vol. 2, 2nd ed.). Edagricole.

Polatschek, A. (1971). Die Verwandtschaftsgruppe um

Poldini, L. (2002). Nuovo Atlante corologico delle piante vascolari nel Friuli Venezia Giulia (2nd ed.). Regione autonoma Friuli Venezia Giulia, Azienda Parchi e Foreste Regionali; Università degli studi di Trieste, Dipartimento di Biologia.

Portik, D. M. , Leaché, A. D. , Rivera, D. , Barej, M. F. , Burger, M. , Hirschfeld, M. , Rödel, M.‐O. , Blackburn, D. C. , & Fujita, M. K. (2017). Evaluating mechanisms of diversification in a Guineo‐Congolian tropical forest frog using demographic model selection. Molecular Ecology, 26(19), 5245–5263. 10.1111/mec.14266 PubMed DOI

Pritchard, J. K. , Stephens, M. , & Donnelly, P. (2000). Inference of population structure using multilocus genotype data. Genetics, 155(2), 945–959. 10.1111/j.1471-8286.2007.01758.x PubMed DOI PMC

Prokhanov, Y. I. (1949). Genus 856.

R Core Team . (2022). R: A language and environment for statistical computing . R Foundation for Statistical Computing. https://www.R‐project.org/

Rambaut, A. (2014).

Rambaut, A. , Drummond, A. J. , Xie, D. , Baele, G. , & Suchard, M. A. (2018). Posterior summarization in Bayesian phylogenetics using tracer 1.7. Systematic Biology, 67(5), 901–904. 10.1093/sysbio/syy032 PubMed DOI PMC

Razgour, O. , Forester, B. , Taggart, J. B. , Bekaert, M. , Juste, J. , Ibáñez, C. , & Manel, S. (2019). Considering adaptive genetic variation in climate change vulnerability assessment reduces species range loss projections. Proceedings of the National Academy of Sciences of the United States of America, 116, 10418–10423. PubMed PMC

Rešetnik, I. , Frajman, B. , & Schönswetter, P. (2016). Heteroploid PubMed DOI

Riina, R. , Peirson, J. A. , Geltman, D. V. , Molero, J. , Frajman, B. , Pahlevani, A. , Barres, L. , Morawetz, J. J. , Salmaki, Y. , Zarre, S. , Kryukov, A. , Bruyns, P. V. , & Berry, P. E. (2013). A worldwide molecular phylogeny and classification of the leafy spurges, DOI

Ronikier, M. (2011). Biogeography of high‐mountain plants in the Carpathians: An emerging phylogeographical perspective. Taxon, 60(2), 373–389. 10.1002/tax.602008 DOI

Ronquist, F. , Teslenko, M. , van der Mark, P. , Ayres, D. L. , Darling, A. , Höhna, S. , Larget, B. , Liu, L. , Suchard, M. A. , & Huelsenbeck, J. P. (2012). MrBayes 3.2: Efficient Bayesian phylogenetic inference and model choice across a large model space. Systematic Biology, 61(3), 539–542. 10.1093/sysbio/sys029 PubMed DOI PMC

Ru, Y. , Mandáková, T. M. , Lysak, M. A. , & Koch, M. A. (2022). The evolutionary history of PubMed DOI PMC

Salzmann, U. , Haywood, A. M. , Lunt, D. J. , Valdes, P. J. , & Hill, D. J. (2008). A new global biome reconstruction and data‐model comparison for the middle Pliocene. Global Ecology and Biogeography, 17(3), 432–447. 10.1111/j.1466-8238.2008.00381.x DOI

Shaw, J. , Lickey, E. B. , Schilling, E. E. , & Small, R. L. (2007). Comparison of whole chloroplast genome sequences to choose noncoding regions for phylogenetic studies in angiosperms: The tortoise and the hare III. American Journal of Botany, 94(3), 275–288. 10.3732/ajb.94.3.275 PubMed DOI

Simmons, M. P. , & Ochoterena, H. (2000). Gaps as characters in sequence‐based phylogenetic analyses. Systematic Biology, 49(2), 369–381. 10.1093/sysbio/49.2.369 PubMed DOI

Slovák, M. , Kučera, J. , Turis, P. , & Zozomová‐Lihová, J. (2012). Multiple glacial refugia and postglacial colonization routes inferred for a woodland geophyte, DOI

Španiel, S. , & Rešetnik, I. (2022). Plant phylogeography of the Balkan Peninsula: Spatiotemporal patterns and processes. Plant Systematics and Evolution, 308(5), 38. 10.1007/s00606-022-01831-1 DOI

Stamatakis, A. (2014). RAxML version 8: A tool for phylogenetic analysis and post‐analysis of large phylogenies. Bioinformatics, 30(9), 1312–1313. 10.1093/bioinformatics/btu033 PubMed DOI PMC

Stark, J. R. , & Fridley, J. D. (2022). Microclimate‐based species distribution models in complex forested terrain indicate widespread cryptic refugia under climate change. Global Ecology and Biogeography, 31(3), 562–575. 10.1111/geb.13447 DOI

Stevanoski, I. , Kuzmanović, N. , Dolenc Koce, J. , Schönswetter, P. , & Frajman, B. (2020). Disentangling relationships between the amphi‐Adriatic DOI

Stojilkovič, V. , Záveská, E. , & Frajman, B. (2022). From Western Asia to the Mediterranean Basin: Diversification of the widespread PubMed DOI PMC

Sun, Y. , Fung, K.‐P. , Leung, P.‐C. , & Shaw, P.‐C. (2005). A phylogenetic analysis of PubMed DOI

Svenning, J.‐C. , Normand, S. , & Kageyama, M. (2008). Glacial refugia of temperate trees in Europe: Insights from species distribution modelling. Journal of Ecology, 96(6), 1117–1127. 10.1111/j.1365-2745.2008.01422.x DOI

Swofford, D. (2002). PAUP*. Phylogenetic analysis using parsimony (*and other Methods). Version 4.0b10. In DOI

Szabó, P. , Kovács, J. , Kocsis, L. , Vennemann, T. , Domingo, L. , Újvári, G. , Halmai, Á. , Pirkhoffer, E. , & Codrea, V. (2022). Pliocene—Early Pleistocene continental climate and vegetation in Europe based on stable isotope compositions of mammal tooth enamel. Quaternary Science Reviews, 288, 107572. 10.1016/j.quascirev.2022.107572 DOI

Tel‐zur, N. , Abbo, S. , Myslabodski, D. , & Mizrahi, Y. (1999). Modified CTAB procedure for DNA isolation from epiphytic cacti of the genera DOI

Thuiller, W. , Lafourcade, B. , Engler, R. , & Araújo, M. B. (2009). BIOMOD – A platform for ensemble forecasting of species distributions. Ecography, 32(3), 369–373. 10.1111/j.1600-0587.2008.05742.x DOI

Tiffney, B. H. , & Manchester, S. R. (2001). The use of geological and paleontological evidence in evaluating plant phylogeographic hypotheses in the Northern Hemisphere Tertiary. International Journal of Plant Sciences, 162(S6), S3–S17. 10.1086/323880 DOI

Török, K. , Podani, J. , & Borhidi, A. (1989). Numerical revision of the Fagion illyricum alliance. Vegetatio, 81(1), 169–180. 10.1007/BF00045522 DOI

Trinajstić, I. (1992). A contribution to the phytogeographical classification of the Illyrian floral element. Acta Botanica Croatica, 51, 135–142.

Turrill, W. B. (1929). The plant life of the Balkan peninsula. A phytogeographical study. Clarendon Press.

Urbaniak, J. , Kwiatkowski, P. , & Ronikier, M. (2018). Postglacial history and current population genetic diversity of a central‐European forest plant DOI

Warren, D. L. , Glor, R. E. , & Turelli, M. (2008). Environmental niche equivalency versus conservatism: Quantitative approaches to niche evolution. Evolution, 62(11), 2868–2883. 10.1111/j.1558-5646.2008.00482.x PubMed DOI

Wen, J. (2001). Evolution of Eastern Asian–Eastern North American biogeographic disjunctions: A few additional issues. International Journal of Plant Sciences, 162(S6), S117–S122. 10.1086/322940 DOI

Wickham, H. (2016). Ggplot2 (2nd ed.). Springer.

Williamson, P. G. (1987). Selection or constraint?: A proposal on the mechanism for stasis *. In Rates of evolution (p. 334). Routledge. 10.4324/9780429293849-6 DOI

Willner, W. , Di Pietro, R. , & Bergmeier, E. (2009). Phytogeographical evidence for post‐glacial dispersal limitation of European beech forest species. Ecography, 32(6), 1011–1018. 10.1111/j.1600-0587.2009.05957.x DOI

Willner, W. , Wessely, J. , Gattringer, A. , Moser, D. , Záveská, E. , Dullinger, S. , Schönswetter, P. , & Hülber, K. (2023). Postglacial range formation of temperate forest understorey herbs – Insights from a spatiotemporally explicit modelling approach. Global Ecology and Biogeography, 32, 1046–1058. PubMed PMC

Yesson, C. , Toomey, N. H. , & Culham, A. (2009). DOI

Záveská, E. , Kirschner, P. , Frajman, B. , Wessely, J. , Willner, W. , Gattringer, A. , Hülber, K. , Lazić, D. , Dobeš, C. , & Schönswetter, P. (2021). Evidence for glacial refugia of the forest understorey species PubMed DOI PMC

Zimmermann, W. , Hegi, G. , & Beger, H. (1924). Euphorbiaceae. In Illustrierte Flora von Mittel‐Europa: Vol. V (pp. 113–188). Lehmanns Verlag.

Spatiotemporal Diversification of Forest Understorey Species Reveals the Existence of Multiple Pleistocene Forest Refugia in Central Europe

Post-glacial range formation of temperate forest understorey herbs - Insights from a spatio-temporally explicit modelling approach