Understanding biogeographical patterns in the western Balkan Peninsula using environmental niche modelling and geostatistics in polymorphic Edraianthus tenuifolius

The Balkan Peninsula represents one of the three southern European glacial refugia where biodiversity persisted throughout the climatically unstable Quaternary. This study considered the ‘refugia within refugia’ model, which assumes the environmental heterogeneity over time and space in larger refugia. To better understand patterns shaped during the Quaternary climatic oscillations, the present and last glacial maximum (LGM) environmental conditions and current morphological variability of Edraianthus tenuifolius, an endemic plant of the western Balkans with a well-known therphical structure, were considered. Potential present and LGM distributions were studied through environmental niche modelling using 161 data of occurrences and six bioclimatic variables, hindcasting the model to LGM conditions using three different global circulation models. To test the geographical variability of the reproductive region, 41 characters of 667 inflorescences from 35 populations within the distribution range were measured. Geographical patterns, using geostatistics together with univariate and multivariate statistical approaches, were analysed. The environmental niche model suggested the current potential distribution in correspondence to its known occurrences. The hindcast to LGM conditions suggested two separate areas of environmental suitability, one in the present-day northern Adriatic coast of Croatia (Istrian Peninsula, Kvarner) and another on the present-day south-eastern Adriatic coast (southern Dalmatia, Montenegro and northern Albania). Morphological variability showed a similar pattern, where southern populations separated from northern populations according to a major split in the central part of its distribution range (central Dalmatia). On other levels, stronger barriers were predicted to separate northern populations from the eastern Istrian Peninsula and the Kvarner area. The results suggested congruent biogeographical patterns to the already known phylogeographical structure. Both environmental niche modelling and the geographical variability of morphological characters suggested spatial partitioning, indicating the potential presence of two separate refugia during the LGM.