Spatial genetic structure and diversity of natural populations of Aesculus hippocastanum L. in Greece
Authors:
Łukasz Walas aff001; Petros Ganatsas aff002; Grzegorz Iszkuło aff001; Peter A. Thomas aff004; Monika Dering aff001
Authors place of work:
Institute of Dendrology, Polish Academy of Sciences, Parkowa, Kórnik, Poland
aff001; Aristotle University of Thessaloniki, School of Forestry and Natural Environment, Laboratory of Silviculture, Thessaloniki, Greece
aff002; Faculty of Biological Sciences, University of Zielona Góra, Prof. Z. Szafrana, Zielona Góra, Poland
aff003; School of Biological Sciences, Keele University, Staffordshire, United Kingdom
aff004; Harvard Forest, Harvard University, Petersham, MA, United States of America
aff005; Faculty of Forestry, Poznań University of Life Sciences, Wojska Polskiego, Poznań, Poland
aff006
Published in the journal:
PLoS ONE 14(12)
Category:
Research Article
doi:
https://doi.org/10.1371/journal.pone.0226225
Summary
Horse-chestnut (Aesculus hippocastanum L.) is an endemic and relict species from the Mediterranean biodiversity hotspot and a popular ornamental tree. Knowledge about the evolutionary history of this species remains scarce. Here, we ask what historical and ecological factors shaped the pattern of genetic diversity and differentiation of this species. We genotyped 717 individuals from nine natural populations using microsatellite markers. The influence of distance, topography and habitat variables on spatial genetic structure was tested within the approaches of isolation-by-distance and isolation-by-ecology. Species niche modeling was used to project the species theoretical range through time and space. The species showed high genetic diversity and moderate differentiation for which topography, progressive range contraction through the species’ history and long-term persistence in stable climatic refugia are likely responsible. A strong geographic component was revealed among five genetic clusters that are connected with very limited gene flow. The environmental variables were a significant factor in the spatial genetic structure. Modeling results indicated that future reduction of the species range may affect its survival. The possible impact of climate changes and high need of in situ conservation are discussed.
Keywords:
Population genetics – Species diversity – Seedlings – Climate change – Inbreeding – Paleogenetics – Gene flow – Paleoclimatology
Introduction
During the Cenozoic (66–23.03 Ma) period, Arcto-Tertiary flora covered a large part of the Northern Hemisphere and formed a unique type of forest ecosystem [1–3]. Species of Tertiary humid-temperate forests began to decline rapidly along with the aridization that appeared in the middle Miocene and were finally pushed out during cold periods of the Pleistocene into warm and humid refugia located in the Central Asia and the West Coast of the North America [4]. Many Tertiary species found safe shelter in the Mediterranean Basin, which is one of the important refugia for species of this ancient flora [5].
Accumulation of endemic taxa drives the formation of biodiversity hotspots and most conservation effort is expended on these endemic species due to their high extinction risk [6, 7]. The Mediterranean Basin is one of the biodiversity hotspots with an exceptional role in the preservation of unique species and genetic diversity [6, 8]. Currently, climate change is becoming a serious and the most urgent problem in this region. A wide spectrum of biological and socio-economic consequences are predicted that impose a challenge for the conservation of the biodiversity stored in the Mediterranean biodiversity hotspot [9–11]. Expected environmental changes may cause extinction of many endemic species which can persist only in a stable climate refugia [12]. Therefore, the protection of the Mediterranean relicts and endemics is an extremely urgent issue and produces a number of challenges.
Firstly, the persistence of Tertiary relicts depends on the conservation of the habitats to which they are frequently highly specialized [13, 14]. In the densely populated Mediterranean, increasing demand on land for agriculture and urban development means that protecting habitats for Tertiary relict and other narrowly distributed plants is difficult [15, 16]. Secondly, long-term survival of relict and endemic species depends on the genetic diversity which is a measure of the potential within a population to respond to natural selection [17, 18]. There are many studies concerning genetic diversity patterns and their causal factors in plants, especially trees, which are key elements of forest ecosystems and in the forest industry [19]. Less attention is generally paid to rare tree species with low economic importance. However, some of these can have high biodiversity impact. Study of relict and endemic species allows better understanding of the processes of extinction and survival, which is especially relevant in the context of expected climate-induced decline in worldwide populations [20–24].
Two basic sets of factors drive population genetic structure and divergence. One of them is the historical biogeography including historical climatic change and tectonic movements that have modified genetic connectivity among populations resulting in significant isolation-by-distance (IBD) [25, 26]. In the Mediterranean region, geological history has left strong and still readable signatures in the current spatial organization of genetic diversity and differentiation [27–30]. The second set of factors affecting the genetic structure are environmental variables including edaphic conditions and topographic heterogeneity that may drive adaptive divergence [31]. Habitat conditions are especially important to endemic and relict species as their populations tend to occur in a very narrow range of environmental conditions [13, 32–34]. Populations existing in isolated areas with dissimilar conditions, such as refugia, may undergo divergent selection that can lead to local adaptation and genetic divergence detectable at neutral variability [35]. This isolation-by-ecology (IBE), referring to significant negative correlation between genetic distance and ecological dissimilarity (edaphic conditions and topography), has recently become recognized as an almost equally important driver of intra-specific differentiation and speciation as the classic IBD model [31, 36, 37].
Horse-chestnut (Aesculus hippocastanum L.) is a Tertiary relict and naturally occurs only in the Balkan Peninsula [38]. During the Pliocene, Aesculus was broadly distributed in Europe and paleaobotanical records prove its presence in Africa and the Caucasus [39, 40]. Climate change in the Pleistocene pushed the species into mountain refugia in the Balkan Peninsula. The natural origin of horse-chestnut remained unknown until the end of 19th century, undoubtedly because of its dispersed and low-density occurrence in the high mountains [38, 41–43]. Paradoxically, the species has been widely cultivated as an ornamental tree in Europe since the 16th century [41].
As a paleoendemic of Tertiary origin, horse-chestnut has experienced a long and rich history during which it has experienced several major geological events that were crucial in the development of spatial genetic structure of tree species in the Mediterranean region [44]. Hence, we postulate that the historical factors might be major determinants of the species’ genetic diversity and differentiation patterns. Additionally, being an endemic species, horse-chestnut may show reduced genetic diversity [45, 46] which would be detrimental to retaining evolutionary potential [47, 48]. Accordingly, the International Union for Conservation of Nature has recommended urgent genetic research on natural stands of horse-chestnut [49].
In its favour, horse-chestnut is distributed across a large mountain system stretching c. 425 km along the Balkan Peninsula, implying the possibility of abrupt environmental gradients. Similar to other Tertiary relicts, the species requires humid conditions [13, 50, 51], and shows some habitat specialization [52]. According to Tsiroukis [53], 62% of natural populations are located in ravines with constantly flowing water and 11% with temporary water. The remaining small fraction is present in dry ravines or on arid and steep slopes (9% and 6%, respectively). Hence, apart from the historical factors, dissimilar ecological conditions might also potentially promote genetic differentiation in horse-chestnut through adaptive divergence, which relates to IBE [54, 55].
In this study we examined genetic diversity and differentiation within natural populations of horse-chestnut in Greece. We integrated genetic methods with ecological niche modeling to better understand the processes and factors that have governed species evolution since the Pleistocene. Such analyses could help to disentangle the roles of historical and ecological factors in shaping the species’ spatial genetic structure. We hypothesized that: (1) spatial genetic structure of horse-chestnut is mostly determined by historical factors, and the species natural populations are characterized by low genetic variability within and high genetic diversity among populations as expected for a palaeoendemic; (2) the range of horse-chestnut will be reduced in the future due to climate change.
Material and methods
Species description and population sampling
Horse-chestnut is native to mountainous regions of the Balkan Peninsula, particularly Greece, where it occurs at altitudes between 228 m and 1485 m a.s.l. [38]. The species is andromonoecious with the majority of flowers functionally male (c. 73%) [56]. Flowers, which are pollinated by insects, are an important source of nectar [57, 58]. However, horse-chestnut is ambophilous since some pollen can also be spread by wind [59]. Large, heavy seeds are dispersed mostly by gravity [60]. They may be secondarily dispersed by water since natural populations occur mostly along streams [53]. Seeds of the related Japanese horse-chestnut (Aesculus turbinata Blume) travel up to 14.5 m from the parent tree with a mean distance 12.2–44.7 m [61]. These distances are probably similar in horse-chestnut.
Due to fragmentary information of species distribution in Balkans, collection of the material in this study was limited to the core range of the species in Greece. Sites with the highest number of mature individuals were selected for investigation using data from Avtzis et al. [38]. Since most of the stands were characterized by small populations [38], we included nine of these in our study to ensure reliable statistical estimations based on microsatellite markers. Material was collected in June 2015. All field sampling was organized by the Laboratory of Silviculture of the Aristotle University of Thessaloniki, and was carried out in close collaboration with the local Forest Services. More specifically, field sampling was carried out after communication with the Directors of the local Forest Districts, and always arranging a local staff accompanies the research team for sampling procedure. Horse-chestnut is not a protected species in Greece and materials were not collected in protected areas, therefore special permissions were not required. In total, 717 mature individuals were sampled (Table 1, S1 Table). Additionally, to analyze how genetic diversity is transferred from adult trees to offspring, 191 one-year-old seedlings were collected in three populations that showed good regeneration (Karitsa I, Karitsa II and Perivoli). All sampled populations are distributed in three mountain ranges: Pindos Mts. (Ondria, Kalampaka, Dasos Nanitsa, Vaeni, Vathirrevma, Perivoli), Ossa Massif (Karitsa I and Karitsa II) and the more southerly located Parnassus Massif (Mariolata). Most of analyzed populations grew in the vicinity of streams; the exceptions were populations Ondria and Mariolata.
We noticed root suckers of Aesculus during the field study, but a full determination of the connections between individuals was not possible. We used a random pattern of sampling but estimated the level of clonality in each population using a clonal diversity parameter (genotypic richness) R = (G—1) / (N—1), in which G denotes the number of distinct genotypes scored in a population (genets) and N is the total number of trees sampled (ramets) [62]. The values of R range from 1, where each genotype is unique, to 0, when all sampled individuals have the same genotype (Table 1).
DNA isolation and genotyping
Genomic DNA was extracted from leaf tissue using protocols described by Dumolin et al. [63]. To determine the genetic variation and genetic diversity, 10 nuclear microsatellite loci (nSSRs) developed for A. turbinata were tested [64]. Preliminary analysis allowed final selection of eight polymorphic loci: AT3D6, AT5D2, AT6D2, AT6D8, AT6D11, AT6D12, AT7D1 and AT7D8. Two rejected primer pairs (AT6D17 and AT5D10) gave low quality products. PCRs were conducted in a final volume of 10 μL, containing 1 × reaction buffer with 2.5 mM MgCl2, 2 μM of dNTP mix, 0.5 U of VivaTaq polymerase (Novazym, Poznań, Poland), 0.4 μM of each starters and 100 ng of DNA. Reactions were conducted using the following thermal protocol: initial denaturation at 94°C/12 min., followed by 35 cycles of denaturation at 94°C/30 s, annealing at 52°C or 55°C (specific to locus) for 30 s, elongation at 72°C/60 s and final elongation at 72°C for 5 min. Products of amplification were analyzed using a 3130 Genetic Analyzer (Applied Biosystems, Foster City, California, USA) with internal size standard GeneScan LIZ-500 and genotypes were scored using GENEMAPPER vs. 4.0 (Applied Biosystems, Foster City, California, USA).
Genetic diversity and differentiation
All clones were removed before analysis of genetic data, identified using the “Find Clones” procedure in GENEALEX 6.4 [65]. This software was also used to estimate basic multilocus within-population diversity estimates such as mean number of alleles (A), effective number of alleles (Ae) and number of private alleles (Ap). INEST v. 2.0 [66] was used to calculate observed (Ho) and expected (He) heterozygosity while FSTAT v 2.9.3. [67] was used to compute allelic richness (Ar). A Bayesian approach implemented in INEST software was applied to estimate the inbreeding coefficient (FIS) including a ‘null alleles’ correction according to the individual inbreeding model (IIM). The estimation was run with 500,000 MCMC cycles with every 200th updated and a burn-in of 50,000. The Deviance Information Criterion (DIC) was used to compare the full model (‘nfb’, when FIS>0) with the random mating model (‘nb’, when FIS = 0) to assess the determinants of levels of homozygosity. The significance of the heterozygote deficiency in a sampled population was tested by U test [68] in GENEPOP [69], and p-values were obtained with the Markov chain algorithm using default settings. Frequency of the null alleles (Null) and Wright’s fixation index (FST) were estimated in FREENA with ENA (Excluding Null Alleles) correction [70]. This method was used to correct positive bias inducted by the presence of the null alleles. The hypothesis on the bottleneck effect in the demographic history of the horse-chestnut populations was tested with the M-ratio approach implemented in INEST.
Subdivision of the horse-chestnut gene pool was evaluated by a non-spatial Bayesian clustering model implemented in STRUCTURE 2.3.4. [71]. The procedure for each of the 10 independent runs included 104 of burn-in and 105 MCMC iterations with the maximum number of clusters set to K = 10. The STRUCTURE model assumed correlated allele frequencies within populations and allowed for mixed ancestry of individuals. To estimate the best-supported number of clusters (S1 Fig; see Supplemental Data with this article), Evanno's delta K method implemented in CLUMPAK [72] was used. Individuals that had an assignment for each cluster of less than 70% were recognized as admixed. In order to detect hidden population-substructure, additional STRUCTURE analysis for separate populations with K = 5 was performed with the same initial conditions. INSTRUCT software [73] was used to verify STRUCTURE results and to estimate selfing rate in each genetic cluster. Analyze was performed for K = 5 as two independent chains in mode 2 (infer populations structure and selfing rates) with Adaptive Independence Sampler, 106 MCMC iterations and 5×105 of burn-in.
In addition to the model-based approach, a discriminant analysis of principal components (DAPC) [74] was performed to infer homogenous genetic clusters. The main advantage of this method is that it works in the absence of any assumptions related to population genetic models and it is a quick computation. DAPC was conducted in the package ‘adegenet’ in R 3.4.3 [75, 76]. In contrast to STRUCTURE, this analytical method can detect hierarchical patterns in the spatial genetic structure and does not require a priori group definition [74]. DAPC is a multivariate two-stage procedure in which data are firstly transformed by principal component analysis (PCA) to remove correlation between variables, which are then submitted to discriminant analysis (DA). This part of the procedure aims at partitioning the total genetic variability into between-cluster and within-cluster components. The procedure optimizes the variance for the first component and minimizes the within-cluster component to obtain the best discriminative power and to define final clusters. To infer the optimal number of the genetic clusters function ‘find.cluster’ was used, which runs successive K-means clustering with increasing number of clusters K from 2 to 50 supported by Bayesian Information Criterion (BIC). Function ‘xvalDapc’ with 30 replications was used to perform cross-validation on varying numbers of principle components. DAPC was conducted with the function ‘dapc’. Fourteen PCs and an optimal number of clusters (10) were applied in estimations. Output was visualized using library ‘ggplot2’ in R [77]. The program POPULATIONS 1.2.30 was used to construct a phylogenetic tree, which shows the relationships among all genetic clusters determined by STRUCTURE [78]. A Cavalli-Sforza and Edwards (Dc) genetic distances was used to obtain a neighbor-joining (NJ) tree.
BAYESASS 3.0 (BA3) was used to infer recent migration rates among studied populations [79]. BA3 is based on the Bayesian MCMC method and estimates migration over the last few generations without assumptions on the Hardy-Weinberg equilibrium. It works best when migration is relatively low [80], which was expected in these isolated populations of horse-chestnut. To obtain a recommended acceptance rate within the range of 20–40% [81], MCMC mixing parameters were initially adjusted during the testing stage for inbreeding coefficient (0.30) and migration rate (0.75) which resulted in acceptance rates of 23% and 28%, respectively. Default settings of mixing parameters for allele frequencies were used. Afterwards, the final analysis was conducted with 107 MCMC iterations, 106 of burn-in and sapling every 1000th iteration. To check consistency, three independent runs with different random seeds were carried out. Finally, a script of Meirmans [80] applied in R [76] was used to calculate the Bayesian deviance for each run. The run with the lowest value of the Bayesian deviance was chosen and used for our interpretations.
Historical gene flow among populations was investigated with MIGRATE-n v. 3.2.6 [82, 83]. This tool computes mutation-scaled historical effective population size (Θ) and migration rates (M) using a coalescence approach. All runs used the Bayesian inference method and the Brownian motion mutation model with identical priors and parameter values, and identical mutation rates were assumed among loci. Testing runs were conducted to set uniform priors (minimum, maximum, delta) for Θ and M (0; 10000; 1000) for all subsequent runs. Each run finally consisted of 50,000 recorded steps at an increment of 100 steps, after a burn-in of 20,000 steps, and a static heating scheme (four chains set at 1, 1.5, 3, 105). We performed three independent runs with different initial seed numbers to verify consistency and used the Bezier approximation for the marginal likelihood to test which run best fit the data [83].
Climate niche modeling
The program MAXENT version 3.3.2 [84, 85] was used to build a climate niche model. This software uses a maximum entropy presence-only model to estimate probability of species occurrence. The dataset of 74 known localities (S1 Table) was collected from the available literature on horse-chestnut distribution [38, 42, 86–88]; each population was represented by a single point analyses. Bioclimatic variables with a resolution of 30 arc-sec were downloaded from the WorldClim database [89] but seven out of 19 variables were excluded because of strong correlation (S2 Table) evaluated using ENMTOOLS v1.3. [90]. Variables for the Mediterranean Basin and Europe were used. Analyses were performed for the present day as well as for two past periods (maximum glaciation c. 22 000 years ago and mid-Holocene c. 6 000 years ago) and for future conditions. Three scenarios of climate changes were used for the future projection: RCP 2.6 (Representative Concentration Pathway)—increase of average temperature by 1°C before the year 2065; RCP 4.5: + 1.4°C before the year 2065; and RCP 8.5: + 2°C before the year 2065 [91]. The Community Climate System Model (CCSM) of global climate was used [92]. Analyses were performed as a bootstrap with 100 replicates; in each replication 20% of the data were set aside as test points and a ‘random seed’ option was applied. The maximum iterations were set to 10,000, convergence threshold to 0.00001 and output was set to logistic. A Receiver Operating Characteristic (ROC) curve and Value of Area Under the Curve (AUC) were used to evaluate the model [93, 94]. Models included areas of Europe, North Africa and the Middle East. Results of analyses were visualized in QGIS 2.18.20 ‘Las Palmas’ [95].
Isolation-by-ecology and isolation-by-distance
A connectivity-resistance approach implemented in the software CIRCUITSCAPE v.4.0.5 [96] was applied to test the significance of the three possibly important factors in gene flow among horse-chestnut populations (distance, topographic complexity and habitat suitability) within the approach of IBE. CIRCUITSCAPE is based on the electrical circuit theory. It models the gene flow (connectivity) between pairs of populations as the resistance distance between nodes connected by resistors. The resistance surfaces are created by different landscape properties that may impede or modify the gene flow and the landscape itself is considered as an electric network. The genetic distance that is built between populations, and is attributable to landscape resistance, is termed isolation-by-resistance (IBR). The unique property of the CIRCUITSCAPE is that each added additional pathway and connection between nodes reduces the resistance. In other words, the conductance of the resistors reflects the intensity of gene flow between populations [97].
Three different resistance surfaces were developed as matrices. A matrix for isolation-by-distance model (IBD) was created using completely flat landscapes in the program QGIS 2.14.21 [95]; this software was also used for creating a relief raster to topographic complexity calculations. To test if the past or present environmental conditions were relevant for observed differentiation pattern (i.e. IBE), past and present layers of the niche suitability previously obtained with MAXENT were used [84, 85]. Matrices of habitat suitability were used as conductance surfaces. All matrices were tested against the matrices of the genetic distance (pairwise FST and FSTNA values) with a multiple matrix regression approach implemented in R (function ‘mantel.rtest’).
Results
Genetic diversity
Out of 717 individuals genotyped, 686 unique genotypes were identified (Table 1). In most populations distinct genotypes were detected (R = 1) or only a few individuals (1–3) were of vegetative origin (Table 1). The highest level of clonality (R = 0.82) was noted in Karitsa I, where 94 unique genotypes were detected among 114 trees sampled.
Screening the sampled individuals with eight nSSR gave 124 different alleles in total and number of alleles varied from 48 in Kalampaka to 88 in Karitsa I—mature (S3 Table). The number of alleles per locus varied from 1–3 (AT6D11) to 8–18 (AT6D12) and was similar to values for Aesculus turbinata [64]. The highest value of expected heterozygosity (He) was at locus AT7D1 in population Dasos Nanitsa (0.915) while the lowest (0.000) at AT6D11 in populations Kalampaka and Ondria, because this locus was monomorphic in these populations. The highest frequency of null alleles was observed in the AT6D2 locus in population Perivoli–mature (23.7%). The effective number of alleles was high at locus AT7D1 in Dasos Nanitsa and Vathirrevma (10.77 and 10.03, respectively), while in the Karitsa I population the highest values were at the AT6D12 (11.50) and AT7D8 (10.44) loci.
Estimators of the genetic diversity are presented in Table 2. The average number of alleles in a population (A) ranged from 6.00 in Kalampaka to 11.00 in Karitsa I. Average effective number of alleles (Ae) is 4.19, which is much smaller value than average number of alleles (8.73). The highest value of allelic richness (Ar) was observed in Vathirrevma (8.44), whereas the lowest was in Kalampaka (5.93). Observed heterozygosity (Ho) ranged from 0.442 (Mariolata) to 0.636 (Vathirrevma), while gene diversity (He) varied from 0.533 (Mariolata) to 0.746 (Vathirrevma). Significant positive FIS values, suggesting an excess of homozygotes, were observed for all analyzed populations. In Kalampaka, the inbreeding coefficient was highest, reaching 0.1520. According to DIC criterion, in the majority of the populations inbreeding was the likely reason for excess homozygotes. In three populations (Karitsa II–mature, Perivoli–mature and Perivoli–seedlings), null alleles are probably responsible for increased homozygosity level. Generally, null alleles were present in all populations with an average frequency of 6.3%. Private alleles were detected in almost all studied stands, and the highest number (4) was noted for Ondria, located at the northernmost part of the Pindos Mts. (Table 2).
Genetic diversity between mature parental populations and 1-year-old seedlings was very similar as indicated by He values, with the exception of Ar that appeared lower in seedlings (Table 2). However, allelic richness differences between mature and seedlings populations was statistically insignificant (Kruskall-Wallis test, P = 0.513). The allele frequency spectra showed that some alleles were absent in progeny. The average frequency of lost alleles in seedlings from Karitsa I was 4.48%, in Karitsa II 2.88%, and in Perivoli 2.89%. New alleles were also detected in all seedling populations. In Perivoli and Karitsa I the gain of new alleles was low (1.61% and 1.01%, respectively) but a surprisingly high percentage of new alleles, 6.08%, was noted in Karitsa II.
Spatial genetic structure
CIRCUITSCAPE indicated significant IBD (r = 0.49, P = 0.0034) for the set of analyzed populations. Statistically significant global FST values with and without the ENA correction attained similarly moderate values (FST = 0.109, FST = 0.114, respectively), but the pairwise estimations showed a wider range of variation (S4 Table). The highest value of FST with ENA correction was noted between Mariolata and Karitsa II (0.227), while the lowest was between the two closest stands, Karitsa I and Karitsa II (0.035). M-ratio test indicated a significant bottleneck effect in all studied populations except for Dasos Nanitsa (P = 0.125).
Analysis of recent migration rates among nine studied populations with BAYESASS showed that migration is very low (S2 Fig, S5 Table). The average proportion of migrants was 0.011, indicating that only 1.1% of individuals in a population are migrants. Higher migration was, however, detected between Karitsa I and Karitsa II (0.060) and between Vaeni and Kalampaka (0.042). Historical migration, estimated with MIGRATE-n, was more intensive (S3 Fig, S6 Table). Average numbers of migrants per generation (θM / 4) was 133.74. The highest migration was noted from Karitsa I to Karitsa II (1202.21), which can be readily explained by the short distance between these stands. Generally, the gene flow was asymmetrical and more intensive from west (Pindos) to east (Ossa) than in the opposite direction (S3 Fig, S6 Table). The exception to this generality may be the intensive gene flow from Karitsa I to Kalampaka (300.89). The lowest effective population size was in Vaeni (θ = 13.56) and the highest in Karitsa II (θ = 41.65).
The Bayesian analysis of spatial genetic structure made with STRUCTURE defined five homogenic genetic clusters among nine analyzed mature populations (Fig 1, S1 Fig, S7 Table). Three of the clusters consisted of a single population (Perivoli–cluster I, Vathirrevma–cluster II, Mariolata–cluster III), while cluster IV grouped populations from Karitsa I and Karitsa II. Cluster V was formed by the populations of Ondria, Kalampaka and Dasos Nanitsa (Fig 1). The Vaeni population exhibited a wide intermixing within all detected clusters but had a slightly higher percentage of membership to cluster V (34.81%) than to cluster II (30.12%). An interesting result is the presence of a significant genetic admixture in Karitsa I located in the Ossa Massif. Phylogenetic analysis made for five detected genetic clusters showed a division into east populations (Mariolata and Ossa Massif, clusters III and IV) and west populations (three clusters from Pindos Mts.: clusters I, II and V). The cluster from Perivoli turned out to be more similar to the south Pindos cluster than to north Pindos (S4 Fig). Eight out of nine analyzed mature populations showed a division into subpopulations (S5A Fig). In Mariolata two subpopulations were detected, in Vaeni three, in Dasos Nanitsa four; for Karitsa I, Vathirrevma and Perivoli the best K was 5. (S6 Fig). Only in Karitsa I was a substructure not detected (S5A Fig). STRUCTURE analyze for three seedling populations divided young individuals into two clusters, one for seedlings from Perivoli and a second from Karitsa I and II.
INSTRUCT results for K = 5 are very similar to STRUCTURE genetic groups (S5B Fig). Lowest selfing rate (0.095) was estimated for cluster IV (Ossa Massif) and highest (0.316) for cluster III (Mariolata). Values for other clusters ranged between 0.20 and 0.23 (Cluster I = 0.222, Cluster II = 0.226, Cluster V = 0.203).
DAPC revealed a hierarchical population genetic structure among the analyzed populations, with an optimal number of clusters of K = 10. Among those ten clusters, five main clusters can be readily distinguished, which corresponds to results obtained with STRUCTURE (Fig 2, S7 Fig). Two clusters are present only in Perivoli indicating a cryptic substructuring in this population (cluster 2 and 6), one is specific to the isolated population of Mariolata (cluster 3) and another (cluster 1) was found in populations from Ossa Massif (Karitsa I and Karitsa II). Clusters 5 and 7 were mainly in Vathirrevma, while remaining clusters (4, 8, 9, 10) were dispersed across all the analyzed populations, with some tendency to occur in populations from the Pindos Mts. DAPC conducted for seedling populations showed that young individuals from each stand formed a separate genetic cluster.
Isolation by ecology (IBE)
The genetic differentiation pattern was driven firstly by topography (r = 0.59, P = 0.0225), which partially proves IBE. Habitat suitability was not a significant factor of differentiation (r = 0.34, P = 0.010). The theoretical model of gene flow from CIRCUITSCAPE showed generally favorable conditions for gene flow in the core range, and pointed to impeded genetic connectivity with the most isolated population from Mariolata (S10 Fig). The genetic connectivity between the Pindos Mts. and Ossa Massif was also not affected by any physical barrier.
Climate niche modeling
MAXENT models based on our nine sampling locations and 65 records from the literature gave robust predictions on horse-chestnut’s past, present and future distribution (Fig 3, S8 Fig, S9 Fig). All models received AUC (Value of Area Under the Curve) greater than 0.995, indicating a very good fit of the model. The most important factors that limit the species distribution are mean temperature of the wettest quarter (relative contribution 29.0, Table 3), precipitation of the coldest quarter (21.1) and precipitation of the driest month (19.3). Quite important also are isothermality (13.9), which quantifies how large the day-to-night temperatures oscillations are relative to summer-winter temperatures, and precipitation seasonality (11.9).
The predictions of horse-chestnut’s current theoretical range mostly overlap our sampling locations from the Balkan Peninsula, especially the Pindos Mts. (suitability > 0.65). This indicates that our sampling for genetic analysis covered the core of the species potential range. There was some overprediction of the model in an eastern direction where the species does not currently occur. Specifically, the Pontic Mts. in Anatolia were highly supported as a possible area of the species distribution (Fig 3, S9A Fig).
During the period of the Last Glacial Maximum (LGM, c. 22,000 BP), the theoretical range of the species was likely much larger, and a westward expansion is evident (S8A Fig). Apart from the Balkan Peninsula, horse-chestnut could grow in wide areas of the northern coasts of the Mediterranean Sea and the Atlantic coast of Europe. The most suitable (>70%) were the territories by the Gulf of Lion, the Apennines, and the Atlantic coast in the northern part of the Iberian Peninsula (Cantabrian Mts.) and Brittany. In Anatolia, only a narrow stripe in the Pontic Mts. appears as a potential area of distribution but with only a low suitability (< 35%).
During the Mid-Holocene (c. 6000 BP), the theoretical range of horse-chestnut decreased, but it was still wider than the current range (S8B Fig). The major change was the elimination of the marginal areas delineated as less suitable during LGM (< 20%). The model also excluded areas by the Gulf of Lion, which were strongly predicted to have populations of horse-chestnut during the previous period. The main distribution during the Holocene were the Pindos Mts. on the Balkan Peninsula and the Cantabrian Mts. on the Iberian Peninsula. The western part of Anatolia appears as a suitable area (as high as 80%) despite its lack of suitability in the previous period.
In terms of future distribution, RCP 2.6 and RCP 4.5 climate changes scenarios are predicted not to have a significant impact on the central stands of horse-chestnut, although they will reduce the potential range (S9B Fig, S9C Fig). Extreme projected climate change in scenario RCP 8.5 may, however, pose a threat to the natural populations (S9D Fig). The models do not predict a change in bioclimatic conditions, which could result in a shift of the species range to the north.
Discussion
Level of diversity and its factors
Jiménez-Mejías et al. [44] stated that the level of genetic diversity within endemic species from the Mediterranean zone is related more with historic and environmental factors than with their biology. Because horse-chestnut is a paleoendemic and a relict with highly fragmented distribution we expected to find rather low genetic variability due to historical processes leading to range fragmentation and contraction [44]. Nevertheless, diversity assessment based on eight microsatellite loci was moderate and similar to other endemic woody species from the Mediterranean region [98–102]. However, inferences based on a small number of loci that do not represent the whole genome may be biased.
Being a Tertiary relict, horse-chestnut thrived for a long time under the relatively stable environment of climatic refugia. According to our modeling, its occurrence in Greece represents the core species distribution since at least the last glaciation, which could be an important factor in the diversity reported here [103]. Long-term survival keeps the old genetic diversity preserved and allows the new genes to accumulate [103, 104]. Recent reviews support the trend of moderate to high genetic variability among Mediterranean endemics [44, 105] and examples from other regions suggest that this may not be specific to the Mediterranean [106–111]. However, the factors underlying the retention of high genetic diversity in endemics are still rather elusive [51, 112, 113].
Over recent years, many studies have shown the high resilience of trees to negative effects of fragmentation due to their longevity, prolonged juvenile stage and high rate of gene flow. All of these are thought to retard the loss of genetic diversity even in endemic species [102, 107, 114, 115]. Hence, in the case of tree relicts, species biology can be as important as historical factors in its effect of genetic diversity. However, conclusions on high genetic variability based on mature populations may not reflect the true situation [116]. Accordingly, a closer look into differences in diversity levels between mature and seedling generations of horse-chestnut revealed some alarming symptoms of diversity leakage. Generally, two out of the three mature-offspring population pairs that were studied showed a suggested decrease in allelic richness and increased inbreeding, although differences were not statistically significant; and heterozygosity decreased in all offspring populations in comparison to mature ones (Table 2). Comparison of heterozygosity between analyzed loci show similar differences between mature and offspring populations (S3 Table); higher value of expected heterozygosity in Karitsa II–seedlings than in Karitsa II–mature probably results from a small number of analyzed individuals in mature population. Increased homozygosity in young ontogenetic stages is frequently noted, and inbreeding depression is expected to eliminate inbred individuals between seed and adult stages. However, the outflow of the allelic diversity noted in seedlings is worrying, even where the acquisition of new alleles was noted. Detection of new alleles means either that immigrants moved in, or that the new alleles were simply omitted in the mature population due to their low frequency and insufficient sampling.
In all populations of horse-chestnut investigated here (except for Dasos Nanitsa), a bottleneck was inferred. However, this estimation may be rather weak because it used only eight loci, where seven is considered a minimum in this type of analysis [117]. Considering that the range of the species has been reducing since at least last glacial cycle, a bottleneck effect ought to be assumed, although it is not a general rule among relict species [117–120]. Nevertheless, the diversity in studied populations was still maintained at a considerable level probably due to traits generally perceived as favouring high variability in tree species [19, 121, 122]. Longevity along with high fecundity would be the factors that might mitigate the effects of fragmentation on the loss of diversity [19]. High numbers of low-frequency alleles may suggest that populations experienced a less severe bottleneck and are now at the expansion phase.
Spatial genetic structure and its factors
In principle, becoming a paleoendemic involves the shrinking and disintegration of the once widespread ancestral range [5, 123]. Distribution in isolated stands may hamper the gene flow contributing to diversity loss due to genetic drift [48, 124]. Hence, considerable genetic differentiation in this group of plants is expected since population genetic theory predicts that allele frequencies diverge among small, fragmented and isolated populations [125, 126]. Indeed, the value of FST = 0.109 indicates a moderate level of genetic differentiation, comparable to several Mediterranean conifer species that do not have such a fragmented and narrow distribution as horse-chestnut [105, 127, 128]. Another Tertiary relict, Taxus baccata L. that has a much wider distribution was reported to attain a very large among-population divergence [129–131]. However, English yew suffers from strong limitations of gene flow, biparental inbreeding and almost no natural regeneration [132–135]. Poor regeneration and high seedling mortality is a common problem in Mediterranean relict species, for example Frangula alnus Mill. subsp. baetica (E. Rev. and Willk.) [136], Buxus balearica Lam. [137] and Olea europaea L. subsp. laperrinei (Batt. and Trabut) Cif. [138]. In comparison, the demographic situation in horse-chestnut populations seems to be better since abundant natural regeneration was reported in some populations [52].
Detection of significant IBD indicates that geographical separation among horse-chestnut populations is responsible for the differentiation detected. This means that historical factors, such as range contraction induced by climatic transformations, are important drivers of spatial genetic structure in horse-chestnut. Low admixture among populations from distinct geographic regions and generally a very low rate of recent gene flow both argue for limitations in current genetic connectivity among populations, which is a consequence of range disintegration that started even before the last glaciation. This fits with the highest values of FST in pairwise comparisons being found between the most distant and marginally located populations, Mariolata and Karitsa, and the remaining populations; the lowest differentiation was noted between Karitsa I and Karitsa II separated with only c. 6 km.
Both, DAPC and STUCTURE indicated strong divergence of marginal populations (Mariolata, Karitsa I and II) and the very distinct character at Perivoli from the northern Pindos compared to other populations from this part of the mountain range (Fig 2). Populations from southern Pindos were also different from the northern Pindos. Thus, DAPC suggests the existence of different gene pools for horse-chestnut in the northern and southern Pindos Mts. The fact that both southern Pindos populations (Vaeni and Vathirrevma) are located by the rivers of the Ionian Sea watershed, while the northern Pindos’ rivers flow to the Aegean Sea may contribute to the observed differentiation. Results from STRUCTURE were less exact, giving Vaeni an intermediate position between Vathirrevma and the northern Pindos populations (Fig 1).
Phylogeographic studies underline the major role of southern European mountains as the centers of refugial areas during the last glacial cycle [8]. The physiographic complexity of the mountain landscapes offered shelter during the climate harshness of the Pleistocene, especially for Tertiary relicts [4, 5, 29]. However, complex mountainous landscapes may also promote genetic divergence, especially in endemics [31, 44, 139]. Indeed, our IBE analysis indicated that out of the two ecological factors considered; only topography was an important driving factor in horse-chestnut differentiation. The relatively high FST value between Mariolata and the remaining populations (average value 0.195) particularly well exemplifies the issue. Mariolata is located at the southern edge of the species distribution in the Parnassus Range, which is in the southwestern spur of the Pindos Mts. This part, unlike to the main Pindos range which runs along a north-south axis, is more latitudinally-oriented and divided into many smaller mountain ranges [140]. Such spatial orientation of ridges may produce additional obstacles for gene flow [141, 142].
Contrary to our initial assumptions, habitat characteristics were not an important factor in differentiation. We expected to find at least a weak signal of adaptive divergence since some populations inhabit ecologically very dissimilar habitats. Additionally, limited gene flow would be a factor enhancing such divergence. One of the reasons that we failed to detect local adaptation may stem from the limited number of populations included in our investigations or too small an area covered (only the Greek part of the range); only two populations out of nine analyzed showed a distinct edaphic character, being located on steep and dry slopes. For T. baccata, Mayol et al. [131] confirmed a significant influence of environmental variables on the current pattern of genetic differentiation. In T. baccata, survival in spatially isolated refugia under divergent temperature regimes were significant determinants of current genetic differentiation and consequently lead to evolution of two evolutionary lineages adapted to different temperature ranges.
Clusters inferred by STRUCTURE and INSTRUCT were related to specific mountain ranges (Fig 1). Such structuring shows that the gene flow among regions, and even between adjacent populations, is too low in such a complex environment to counteract the effects of the accumulated differentiation (Fig 1). However, gene flow is also strictly related to mechanisms of pollen and seeds dispersal. Hence, species biology may be a limiting factor for gene flow. Pollen movement distance in closely related A. turbinata may reach on average of c. 180 m with maximum distance reported of over 700 m. Similar distance may be attained in horse-chestnut, and it may account for effective within-population gene flow resulting in high within-population diversity. However, this distance is too short to maintain among-population connectivity in so disjunct a range and such a complex environment. Generally, in animal-pollinated tree species, highly mobile pollinators such as birds can work effectively in fragmented landscapes but in the case of insect-pollinated trees, pollen movement is mostly local [143]. Pollen of tropical tree species is moved further by insects, but even here pollen dispersal over more than a few kilometers happens only occasionally [144].
Two populations from the easternmost margins of the species range, Karitsa I and II (Ossa Massif), are separated from the core range although genetic admixture from the northern Pindos was revealed by STRUCTURE. CIRCUITSCAPE. This indicates, at least theoretically, that there are no impenetrable barriers preventing gene exchange between the Pindos Mts. and the Ossa Massif (S10 Fig), and gene flow may be asymmetric—primarily in a west-east direction (S2 Fig, S3 Fig, S10 Fig). We conjecture that this may be related to the hydrochory of species linked to the river network as the rivers from the northern Pindos that run down to the Aegean Sea. There are no direct studies on seed dispersal by water in horse-chestnut, but its populations are restricted to the watercourses which suggests the possibility [53]. In Frangula alnus subsp. beatica and other Tertiary endemics growing by mountain rivers in southern Spain, winter flooding events were shown to play an important role in the secondary dispersal of seeds already spread by birds [136]. Seeds of horse-chestnut are primary barochoric but dispersal by rodents or other animals should be considered [145, 146]. Congeneric A. turbinata is dispersed for distance of up to 114.5 m [61]. Similar dispersal distances can be expected for horse-chestnut as the seed of both species are of comparable size. However, as with wind-mediated gene flow, animal-dispersion of seeds would tend to account for local gene flow while long-distance gene transfer would require the involvement of a more effective medium; here, the mountain rivers would suit perfectly. In their work, Robledo-Arnuncio et al. [141] showed that watersheds were involved in the phylogeographic pattern of Scots pine in the Northern Meseta of Spain.
Other explanations for the genetic admixture observed between the northern Pindos and Ossa Massif and for the overall moderate value of genetic differentiation include the historically larger distribution and more extensive gene exchange in the past. Strong support for this comes from the estimations of the historical migration rate (S3 Fig) and modeling potential distributions (Fig 3). Firstly, MIGRATE-n showed more intensive gene flow among studied populations in the past. Large and demographically stable populations in the past, connected via gene flow and longevity would synergistically buffer horse-chestnut from the negative effects of range fragmentation for a long time. Secondly, modeling of theoretical range clearly shows a possibility of wider distribution in the past. Accordingly, during LGM (ca. 22 ka BP) Aesculus sp. populations might have occupied extensive areas on the west coast of the Mediterranean Sea (Fig 3). Paleodistribution modeling predicted five main domains of species occurrence which were the Pindos Mts., Apennines, areas over the Gulf of Lion, Cantabrian Mts. and the Brittany Peninsula. Except for the two latter locations, fossil data support the presence of horse-chestnut in those areas during the Pleistocene [40, 147]. Mid-Holocene sites of horse-chestnut might even have existed in south-east Ukraine [148], which was also weakly supported by our distribution modeling (1% of suitability, S4 Fig). Based on these, we propose that Pleistocene refugial populations might have existed for horse-chestnut in the Gulf of Lion, Apennines and Pindos Mts., although we cannot distinguish the existence of such population in the Cantabrian Mts. or Brittany.
Lack of palaeobotanic data on the past occurrence of horse-chestnut in Anatolia makes its presence there tentative. MAXENT projected the existence of suitable habitats for horse-chestnut in the mid-Holocene despite it being absent during the last glaciation, which is rather puzzling. Considering its limited colonization ability due to heavy seeds and the mode of their dispersion, the Holocene colonization of these areas from the Balkan Peninsula across the Aegean Sea is highly unlikely. More mysterious is the introduction route of horse-chestnut from Turkey to Europe [41]. It is possible that against the prediction of niche modeling, horse-chestnut grew in Anatolia during the LGM in local, limited refugia from which it expanded in the mid-Holocene. Unfortunately, there are no fossil records from this period and the only late Pliocene pollen record was found in eastern Anatolia [149]. It is possible, that horse-chestnut was not native to Turkey when discovered by European botanists in the 16th century. It could have been introduced to Turkey from Greece or Albania in the past, for example at a time when the Balkan Peninsula was conquered by the Ottoman Empire. However, until new paleaobotanical data are collected, neither of these hypotheses can be verified.
Conservation remarks
Given the high fragmentation of the horse-chestnut range, there is a serious risk that the current level of divergence will be increasing since observed gene flow seems to be insufficient to maintain functional connectivity among the remnants dispersed in a complex mountainous landscape. Horse-chestnut is a habitat specialist, like other paleoendemic species [13], and may not be able to adapt to new conditions. Environmental change, human pressure and especially limitations of the species biology hinders it from rapid colonization and successful competition. Combination of these factors makes the future of this species uncertain.
Our results indicate increased homozygosity in natural populations, implying a risk of inbreeding depression and loss of diversity in the future (Table 2). In most cases, the homozygote excess was due to inbreeding, which may be accounted for by both selfing and biparental inbreeding. Although data on the mating system of horse-chestnut are missing, the conspecific A. turbinata shows a considerable level of selfing (8.3% self-pollinated seedlings) [150, 151]. Horse-chestnut is primarily entomophilous although the concentration of air-borne pollen and existence of fossil pollen records suggest that it may be partially anemophilous as well [39, 59]. The homozygosity excess could be also generated by a Wahlund effect. STRUCTURE revealed hidden population substructure in all (except one) mature populations that may suggest the existence small local breeding groups. The species reproductive biology in terms of limited capacity of pollen/seed dispersion could be the possible cause of the observed substructuring and thus, the homozygosity excess. However, this issue requires further studies.
Many species are currently at high risk from the negative impact of climate change. In the face of changing environmental conditions, endemic species may not be able to keep pace with these changes and therefore may be first to disappear [152]. For example, the range of the Chinese Tertiary relict Davidia involucrata Baill., is expected to be reduced to less than 30% of its current theoretical range [153]. In the Mediterranean region, the distribution of many plant species is limited by water availability [5, 154, 155] and drought is one of the most important factor that leads to the decline of forest in this region [156]. Future climate change in the Mediterranean may bring some risk to natural populations of horse-chestnut, mainly by the reduction of rainfall and increasing summer temperature [157]. According to future range modeling, a moderate scenario of climate change should not cause extinction of the populations in the Pindos Mts., but extreme changes (RCP 8.5 scenario) entails some risk for species survival.
The results obtained in this work are less optimistic about the future persistence of horse-chestnut than previous studies [52]. The most important contributing factors are also different. These result from using a more accurate model. However, it should be remembered that modeling provides only estimations and the differences between various models can be large, especially when considering future conditions [158, 159]. Unfortunately, current models of the future range of horse-chestnut do not show a large expansion of this species to the north to trace the optimum conditions. Thus, in situ conservation should be a priority and a conservation program should be immediately launched for Greek populations, which are the largest [49]. Unfortunately, as observed during this study, most of the natural populations are not protected and some of them are being destroyed by infrastructure expansion. Although the most valuable natural populations of the horse-chestnut are located in Greece, this species is not on the Red List of Greece and a reasonable conservation strategy is lacking [160]. So far horse-chestnut has been included in the national list of protected species of the Presidential Degree 67/1981, and some populations are within the NATURA 2000 network [38], but it is not enough. Much more needs to be done and surely can be done in order to protect this beautiful tree in its natural environment.
Supporting information
S1 Fig [tiff]
Optimal number of clusters based on Evanno's delta K.
S2 Fig [tiff]
Theoretical current gene flow between populations estimated with BAYESASS and visualized using QGIS.
S3 Fig [tiff]
Theoretical historical gene flow between populations estimated with MIGRATE-N and visualized using QGIS.
S4 Fig [tiff]
Neighbor-joining tree based on Cavalli-Sforza and Edwards genetic distances for clusters determined by STRUCTURE.
S5 Fig [tiff]
Results of STRUCTURE and INSTRUCT analyzes.
S6 Fig [tiff]
Best K for each STRUCTURE analyze for separate population, according to ln(Pr(X|K) values.
S7 Fig [tiff]
Ordination plot for the first two discriminant axes obtained DAPC analysis.
S8 Fig [tiff]
Theoretical range of in the past.
S9 Fig [tiff]
Current and future theoretical range of .
S10 Fig [tiff]
Theoretical gene flow in relation to environmental conditions determined using CIRCUITSCAPE.
S1 Table [2012]
Geographical coordinates and bioclimatic variables of 74 natural populations of horse-chestnut used in MAXENT analyze.
S2 Table [docx]
Bioclimatic variables from WorldClim database.
S3 Table [xlsx]
Parameters of genetic diversity of analyzed loci.
S4 Table [docx]
Matrix of genetic distance between populations.
S5 Table [docx]
Matrix of migration rate between populations as calculated by BAYEASS.
S6 Table [docx]
Matrix of migration (number of individuals per generation) between populations as calculated by MIGRATE-n.
S7 Table [docx]
Probability of belonging of each individual to each cluster according to STRUCTURE analysis.
Zdroje
1. Chaney RW. Tertiary centers and migration routes. Ecological Monographs 1947; 17: 139–148.
2. Wolfe JA. Tertiary climates and floristic relationships at high latitudes in the Northern Hemisphere. Palaeogeography, Palaeoclimatology, Palaeoecology 1980; 30: 313–323.
3. Mai DH. Palaeofloristic changes in Europe and the confirmation of the Arctotertiary-Palaeotropical geofloral concept. Review of Palaeobotany and Palynology 1991; 68: 29–36.
4. Milne RI, Abbott RJ. The origin and evolution of tertiary relict floras. Advances in Botanical Research 2002; 38: 281–314.
5. Thompson JD. Plant evolution in the Mediterranean. Oxford: Oxford University Press on Demand; 2005.
6. Myers N, Mittermeier RA, Mittermeier CG, Da Fonseca GA, Kent J. Biodiversity hotspots for conservation priorities. Nature 2000; 403: 853–858. doi: 10.1038/35002501 10706275
7. Mouillot D, Bellwood DR, Baraloto C, Chave J, Galzin R, Harmelin-Vivien M, Kulbicki M, Levergne S, Lavorel S, Mouguet N, Paine CET, Renaud J, Thuiller W. Rare Species Support Vulnerable Functions in High-Diversity Ecosystems. PLOS Biology 2003; 11: e1001569.
8. Feliner GN. Patterns and processes in plant phylogeography in the Mediterranean Basin. A review. Perspectives in Plant Ecology, Evolution and Systematics 2014; 16: 265–278.
9. Sarris D, Christodoulakis D, Körner C. 2007. Recent decline in precipitation and tree growth in the eastern Mediterranean. Global Change Biology 13: 1187–1200.
10. García-Ruiz JM, López-Moreno JI, Vicente-Serrano SM, Lasanta–Martínez T, Beguería S. 2011. Mediterranean water resources in a global change scenario. Earth-Science Reviews 2011; 105: 121–139.
11. Bruno D, Belmar O, Sánchez-Fernández D, Guareschi S, Millán A, Velasco J. Responses of Mediterranean aquatic and riparian communities to human pressures at different spatial scales. Ecological Indicators 2014; 45: 456–464.
12. Harrison S, Noss R. Endemism hotspots are linked to stable climatic refugia. Annals of Botany 2017; 119: 207–214. doi: 10.1093/aob/mcw248 28064195
13. Huang Y, Jacques FMB, Su T, Ferguson DK, Tang H, Chen W, Zhou Z. Distribution of Cenozoic plant relicts in China explained by drought in dry season. Scientific Reports 2015; 5: 14212. doi: 10.1038/srep14212 26369980
14. Molina-Venegas R, Aparicio A, Lavergne S, Arroyo J. Climatic and topographical correlates of plant palaeo-and neoendemism in a Mediterranean biodiversity hotspot. Annals of Botany 2016; 119: 229–238. doi: 10.1093/aob/mcw093 27288510
15. Vogiatzakis IN, Mannion AM, Griffiths GH. Mediterranean ecosystems: problems and tools for conservation. Progress in Physical Geography: Earth and Environment 2006; 30: 175–200.
16. Cox RL, Underwood EC. The Importance of Conserving Biodiversity Outside of Protected Areas in Mediterranean Ecosystems. PLOS ONE 2011; 6: e14508. doi: 10.1371/journal.pone.0014508 21249126
17. Blanco‐Pastor JL, Fernández‐Mazuecos M, Vargas P. Past and future demographic dynamics of alpine species: limited genetic consequences despite dramatic range contraction in a plant from the Spanish Sierra Nevada. Molecular Ecology 2013; 22: 4177–4195. doi: 10.1111/mec.12383 23844700
18. López-Pujol J, Martinell MC, Massó S, Blanché C, Sáez L. The ‘paradigm of extremes’: extremely low genetic diversity in an extremely narrow endemic species, Coristospermum huteri (Umbelliferae). Plant Systematics and Evolution 2013; 299: 439–446.
19. Petit RJ, Hampe A. Some Evolutionary Consequences of Being a Tree. Annual Review of Ecology, Evolution, and Systematics 2006; 37: 187–214.
20. Dirzo R, Young HS, Galetti M, Ceballos G, Isaac NJ, Collen B. Defaunation in the Anthropocene. Science 2014; 345: 401–406. doi: 10.1126/science.1251817 25061202
21. Geng Q, Yao Z, Yang J, He J, Wang D, Wang Z, Liu H. Effect of Yangtze River on population genetic structure of the relict plant Parrotia subaequalis in eastern China. Ecology and Evolution 2015; 5: 4617–4627. doi: 10.1002/ece3.1734 26668727
22. Maharramova EH, Safarov HM, Kozlowski G, Borsch T, Muller LA. Analysis of nuclear microsatellites reveals limited differentiation between Colchic and Hyrcanian populations of the wind-pollinated relict tree Zelkova carpinifolia (Ulmaceae). American Journal of Botany 2015; 102: 119–128. doi: 10.3732/ajb.1400370 25587154
23. McDowell NG, Allen CD. Darcy’s law predicts widespread forest mortality under climate warming. Nature Climate Change 2015; 5: 669.
24. Tamaki S, Isoda K, Takahashi M, Yamada H, Yamashita Y. Genetic structure and diversity in relation to the recently reduced population size of the rare conifer, Pseudotsuga japonica, endemic to Japan. Conservation Genetics 2018; 19: 1243–1255.
25. Haddad NM, Bowne DR, Cunningham A, Danielson BJ, Levey DJ, Sargent S, Spira T. Corridor Use by Diverse Taxa. Ecology 2003; 84: 609–615.
26. Townsend PA, Levey DJ. An Experimental Test of Whether Habitat Corridors Affect Pollen Transfer. Ecology 2005; 86: 466–475.
27. Thompson JD. Population differentiation in Mediterranean plants: insights into colonization history and the evolution and conservation of endemic species. Heredity 1999; 82: 229–236. doi: 10.1038/sj.hdy.6885040 10336696
28. Thompson JD, Lavergne S, Affre L, Gaudeul M, Debussche M. Ecological differentiation of Mediterranean endemic plants. Taxon 2005; 54: 967–976.
29. Médail F, Diadema K. Glacial refugia influence plant diversity patterns in the Mediterranean Basin. Journal of Biogeography 2009; 36: 1333–1345.
30. Hewitt GM. Mediterranean Peninsulas: The Evolution of Hotspots. In: Zachos FE, Habel JC, eds. Biodiversity Hotspots: Distribution and Protection of Conservation Priority Areas, Berlin, Heidelberg: Springer Berlin Heidelberg 2011; 123–147.
31. Noguerales V, Cordero PJ, Ortego J. Hierarchical genetic structure shaped by topography in a narrow-endemic montane grasshopper. BMC Evolutionary Biology 2016; 16: 96. doi: 10.1186/s12862-016-0663-7 27149952
32. Kruckeberg AR, Rabinowitz D. Biological aspects of endemism in higher plants. Annual Review of Ecology and Systematics 1985; 16: 447–479.
33. Lavergne S, Thompson JD, Garnier E, Debussche M. The biology and ecology of narrow endemic and widespread plants: a comparative study of trait variation in 20 congeneric pairs. Oikos 2004; 107: 505–518.
34. Hermant M, Prinzing A, Vernon P, Convey P, Hennion F. Endemic species have highly integrated phenotypes, environmental distributions and phenotype–environment relationships. Journal of Biogeography 2013; 40: 1583–1594.
35. Shafer ABA, Wolf JBW. Widespread evidence for incipient ecological speciation: a meta-analysis of isolation-by-ecology. Ecology Letters 2013; 16: 940–950. doi: 10.1111/ele.12120 23627762
36. Edelaar P, Alonso D, Lagerveld S, Senar JC, Björklund M. Population differentiation and restricted gene flow in Spanish crossbills: not isolation-by-distance but isolation-by-ecology. Journal of Evolutionary Biology 2012; 25: 417–430. doi: 10.1111/j.1420-9101.2011.02443.x 22239460
37. Wang IJ, Bradburd GS. Isolation by environment. Molecular Ecology 2014; 23: 5649–5662. doi: 10.1111/mec.12938 25256562
38. Avtzis ND, Avtzis DN, Vergos SG, Diamandis S. A contribution to the natural distribution of Aesculus hippocastanum (Hippocastanaceae) in Greece. Phytologia Balcanica 2007; 13: 183–187.
39. Maley J. Les changements climatiques de la fin du Tertiaire en Afrique: leur conséquence sur l’apparition du Sahara et de sa végétation. The Sahara and the Nile 1980; 63–86.
40. Mijarra JMP, Manzaneque FG, Morla C. Survival and long-term maintenance of tertiary trees in the Iberian Peninsula during the Pleistocene: first record of Aesculus L. (Hippocastanaceae) in Spain. Vegetation History and Archaeobotany 2008; 17: 351.
41. Lack HW. The Discovery and Rediscovery of the Horse Chestnut. Arnoldia 2002; 61: 15–19.
42. Peçi D, Mullaj A, Dervishi A. The natural distribution of horse chestnut (Aesculus hippocastanum L) in Albania. Journal of Institute Alb-Shkenca 2012; 5: 153–157. [in Albanian]
43. Prada D, Velloza TM, Toorop PE, Pritchard HW. Genetic population structure in horse chestnut (Aesculus hippocastanum L.): effects of human-mediated expansion across Europe. Plant Species Biology 2011; 26: 43–50.
44. Jiménez-Mejías P, Fernández-Mazuecos M, Amat ME, Vargas P. Narrow endemics in European mountains: high genetic diversity within the monospecific genus Pseudomisopates (Plantaginaceae) despite isolation since the late Pleistocene. Journal of Biogeography 2015; 42: 1455–1468.
45. Ellstrand NC, Elam DR. Population Genetic Consequences of Small Population Size: Implications for Plant Conservation. Annual Review of Ecology and Systematics 1993; 24: 217–242.
46. Jennings H, Wallin K, Brennan J, Valle AD, Guzman A, Hein D, Hunter S, Lewandowski A, Olson S, Parsons H, Scheidt S, Wang Z, Werra A, Kartzinel RY, Givnish TJ. Inbreeding, low genetic diversity, and spatial genetic structure in the endemic Hawaiian lobeliads Clermontia fauriei and Cyanea pilosa ssp. longipedunculata. Conservation Genetics 2016; 17: 497–502.
47. Ostfeld RS, Keesing F. Effects of Host Diversity on Infectious Disease. Annual Review of Ecology, Evolution, and Systematics 2012; 43: 157–182.
48. Frankham R, Bradshaw CJ, Brook BW. Genetics in conservation management: revised recommendations for the 50/500 rules, Red List criteria and population viability analyses. Biological Conservation 2014; 170: 56–63.
49. Allen D, Khela S. Aesculus hippocastanum (errata version published in 2018). The IUCN Red List of Threatened Species 2017. 2017. Available at: https://doi.org/e.T202914A122961065
50. Romo A, Iszkuło G, Seghir Taleb M, Walas Ł, Boratyński A. Taxus baccata in Morocco: a tree in regression in its southern extreme. Dendrobiology 2017; 78: 63–74.
51. Maharramova E, Huseynova I, Kolbaia S, Gruenstaeudl M, Borsch T, Muller LA. Phylogeography and population genetics of the riparian relict tree Pterocarya fraxinifolia (Juglandaceae) in the South Caucasus. Systematics and Biodiversity 2018; 16: 14–27.
52. Walas Ł, Dering M, Ganatsas P, Pietras M, Pers-Kamczyc E, Iszkuło G. The present status and potential distribution of relict populations of Aesculus hippocastanum L. in Greece and the diverse infestation by Cameraria ohridella Deschka & Dimić. Plant Biosystems 2018; 152: 1048–1058.
53. Tsiroukis A. Reproductive biology and ecology of horse chestnut (Aesculus hippocastanum L.) [in Greek]. 2008. Department of Botany, Faculty of Biology, National & Kapodistrian University of Athens.
54. Orsini L, Vanoverbeke J, Swillen I, Mergeay J, De Meester L. Drivers of population genetic differentiation in the wild: isolation by dispersal limitation, isolation by adaptation and isolation by colonization. Molecular Ecology 2013; 22: 5983–5999. doi: 10.1111/mec.12561 24128305
55. Lafontaine G, Napier JD, Petit RJ, Hu FS. Invoking adaptation to decipher the genetic legacy of past climate change. Ecology 2018; 99: 1530–1546. doi: 10.1002/ecy.2382 29729183
56. Weryszko-Chmielewska E, Chwil M. Structure of floral nectaries in Aesculus hippocastanum L. Acta Botanica Croatica 2017; 76: 41–48.
57. Free JB. The Flower Constancy of Honeybees. Journal of Animal Ecology 1963; 32: 119–131.
58. Kevan PG. How large bees, Bombus and Xylocopa (Apoidea Hymenoptera) forage on trees: optimality and patterns of movement in temperate and tropical climates. Ethology Ecology & Evolution 1990; 2: 233–242.
59. Weryszko-Chmielewska E, Tietze M, Michońska M. Ecological features of the flowers of Aesculus hippocastanum L. and characteristics of Aesculus L. pollen seasons under the conditions of central-eastern Poland. Acta Agrobotanica 2012; 65: 61–68.
60. Thomas PA, Alhamd O, Iszkuło G, Dering M, Mukassabi TA. Biological Flora of the British Isles: Aesculus hippocastanum. Journal of Ecology 2019; 107: 992–1030.
61. Hoshizaki K, Suzuki W, Nakashizuka T. Evaluation of secondary dispersal in a large-seeded tree Aesculus turbinata: a test of directed dispersal. Plant Ecology 1999; 144: 167–176.
62. Dorken ME, Eckert CG. Severely reduced sexual reproduction in northern populations of a clonal plant, Decodonverticillatus (Lythraceae). Journal of Ecology 2001; 89: 339–350
63. Dumolin S, Demesure B, Petit RJ. Inheritance of chloroplast and mitochondrial genomes in pedunculate oak investigated with an efficient PCR method. Theoretical and Applied Genetics 1995; 91: 1253–1256. doi: 10.1007/BF00220937 24170054
64. Minami E. Polymorphic microsatellite markers in Japanese horse chestnut Aesculus turbinata Blume. Molecular Ecology 1998; 7: 1616–1617. 9819913
65. Peakall R, Smouse PE. GENALEX 6: genetic analysis in Excel. Population genetic software for teaching and research. Molecular Ecology Notes 2006; 6: 288–295.
66. Chybicki, I. INEST 2.0 [Computer Software]. 2016.
67. Goudet J. FSTAT, a program to estimate and test gene diversity and fixation indices (version 2.9.3). 2001. Available at: http://www2.Unil.Ch/Popgen/Softwares/Fstat.Htm
68. Guo SW, Thompson EA. Performing the exact test of Hardy-Weinberg proportion for multiple alleles. Biometrics 1992; 361–372. 1637966
69. Rousset F. GENEPOP’007: a complete re-implementation of the genepop software for Windows and Linux. Molecular Ecology Resources 2008; 8: 103–106. doi: 10.1111/j.1471-8286.2007.01931.x 21585727
70. Chapuis M-P, Estoup A. Microsatellite Null Alleles and Estimation of Population Differentiation. Molecular Biology and Evolution 2007; 24: 621–631. doi: 10.1093/molbev/msl191 17150975
71. Pritchard JK, Stephens M, Donnelly P. Inference of population structure using multilocus genotype data. Genetics 2000; 155: 945–959. 10835412
72. Kopelman NM, Mayzel J, Jakobsson M, Rosenberg NA, Mayrose I. CLUMPAK a program for identifying clustering modes and packaging population structure inferences across K. Molecular Ecology Resources 2015; 15: 1179–1191. doi: 10.1111/1755-0998.12387 25684545
73. Gao H, Williamson S, Bustamante CD. A Markov chain Monte Carlo approach for joint inference of population structure and inbreeding rates from multilocus genotype data. Genetics 2007; 176: 1635–1651. doi: 10.1534/genetics.107.072371 17483417
74. Jombart T, Devillard S, Balloux F. Discriminant analysis of principal components: a new method for the analysis of genetically structured populations. BMC Genetics 2010; 11: 94. doi: 10.1186/1471-2156-11-94 20950446
75. Jombart T. adegenet: a R package for the multivariate analysis of genetic markers. Bioinformatics 2008; 24: 1403–1405. doi: 10.1093/bioinformatics/btn129 18397895
76. R Core Team. R: A language and environment for statistical computing. R Foundation for Statistical Computing, Vienna, Austria. 2013. Available at: http://www.R-project.org/
77. Wickham H. ggplot2: Elegant Graphics for Data Analysis. Springer-Verlag New York 2016.
78. Langella O. POPULATIONS, a free population genetics software. 2002. URL http://www.legs.cnrs-gif.fr.
79. Wilson GA, Rannala B. Bayesian inference of recent migration rates using multilocus genotypes. Genetics 2003; 163: 1177–1191. 12663554
80. Meirmans PG. Nonconvergence in Bayesian estimation of migration rates. Molecular Ecology Resources 2014; 14: 726–733. doi: 10.1111/1755-0998.12216 24373147
81. Rannala B, Zhu T, Yang Z. Tail paradox, partial identifiability, and influential priors in Bayesian branch length inference. Molecular Biology and Evolution 2011; 29: 325–335. doi: 10.1093/molbev/msr210 21890479
82. Beerli P, Felsenstein J. Maximum-likelihood estimation of migration rates and effective population numbers in two populations using a coalescent approach. Genetics 1999; 152: 763–773. 10353916
83. Beerli P, Palczewski M. Unified framework to evaluate panmixia and migration direction among multiple sampling locations. Genetics 2010; 185: 313–326. doi: 10.1534/genetics.109.112532 20176979
84. Phillips SJ, Dudík M, Schapire RE. A Maximum Entropy Approach to Species Distribution Modeling. In: Brodley C, eds. Proceedings of the Twenty-first International Conference on Machine Learning. New York, NY, USA: ACM, 2004; 83.
85. Elith J, Phillips SJ, Hastie T, Dudík M, Chee YE, Yates CJ. A statistical explanation of MaxEnt for ecologists. Diversity and Distributions 2011; 17: 43–57.
86. Valade R, Kenis M, Hernandez‐Lopez A, Augustin S, Mena NM, Magnoux E, Rougerie R, Lakatos F, Roques A, Lopez-Vaamonde C et al. Mitochondrial and microsatellite DNA markers reveal a Balkan origin for the highly invasive horse-chestnut leaf miner Cameraria ohridella (Lepidoptera, Gracillariidae). Molecular Ecology 2009; 18: 3458–3470. doi: 10.1111/j.1365-294X.2009.04290.x 19627490
87. Lees DC, Lack HW, Rougerie R, Hernandez-Lopez A, Raus T, Avtzis ND, Augustin S, Lopez-Vaamonde C. Tracking origins of invasive herbivores through herbaria and archival DNA: the case of the horse-chestnut leaf miner. Frontiers in Ecology and the Environment 2011; 9: 322–328.
88. Acevski J, Simovski B. Forest associations of the National Park Mavrovo in the Republic of Macedonia. In: Horodnic SA, Duduman M-L, Palaghianu C, eds. Proceedings of the International Conference Integrated Management of Environmental Resources-Suceava, November 2011. Suceava: Editura Universităţii "Ştefan cel Mare" 2012; 17–27.
89. Hijmans RJ, Cameron SE, Parra JL, Jones PG, Jarvis A. Very high resolution interpolated climate surfaces for global land areas. International Journal of Climatology 2005; 25: 1965–1978.
90. Kukwa M, Kolanowska M. Glacial refugia and the prediction of future habitat coverage of the South American lichen species Ochrolechia austroamericana. Scientific Reports 2016; 6: 38779. doi: 10.1038/srep38779 27929090
91. Collins M, Knutti R, Arblaster J, Dufresne J-L, Fichefet T, Friedlingstein P, Gao X, Gutowski WJ, Johns T, Krinner G, Shongwe M, Tebaldi C, Weaver AJ, Wehner MF, Allen MR, Andrews T, Beyerle U, Bitz CM, Bony S, Booth BBB. Long-term climate change: projections, commitments and irreversibility. In: Climate Change 2013: The Physical Science Basis. IPCC Working Group I Contribution to AR5. Eds. IPCC, Cambridge: Cambridge University Press. 2013.
92. Gent PR, Danabasoglu G, Donner LJ, Holland MM, Hunke EC, Jayne SR, Lawrence DM, Neale RB, Rasch PJ, Vertenstein M, Worley PH, Yang Z-L, Zhang M. The Community Climate System Model Version 4. Journal of Climate 2011; 24: 4973–4991.
93. Wang Z, Chang YI, Ying Z, Zhu L, Yang Y. A parsimonious threshold-independent protein feature selection method through the area under receiver operating characteristic curve. Bioinformatics 2007; 23: 2788–2794. doi: 10.1093/bioinformatics/btm442 17878205
94. Mas J-F, Soares Filho B, Pontius RG, Farfán Gutiérrez M, Rodrigues H. A Suite of Tools for ROC Analysis of Spatial Models. ISPRS International Journal of Geo-Information 2013; 2: 869–887.
95. QGIS Development Team. QGIS Geographic Information System. Open Source Geospatial Foundation Project. 2012. Available at: http://qgis.osgeo.org
96. McRae BH, Beier P. Circuit theory predicts gene flow in plant and animal populations. Proceedings of the National Academy of Sciences 2007; 104: 19885–19890.
97. McRae BH. Isolation by resistance. Evolution 2006; 60: 1551–1561. 17017056
98. Awad L, Fady B, Khater C, Roig A, Cheddadi R. Genetic structure and diversity of the endangered fir tree of Lebanon (Abies cilicica Carr.): implications for conservation. PLoS One 2014; 9: e90086. doi: 10.1371/journal.pone.0090086 24587219
99. Boratyński A, Wachowiak W, Dering M, Boratyńska K, Sękiewicz K, Sobierajska K, Jasińska AK, Klimko M, Montserrat JM, Romo A, Ok T, Didukh Y. The biogeography and genetic relationships of Juniperus oxycedrus and related taxa from the Mediterranean and Macaronesian regions. Botanical Journal of the Linnean Society 2014; 174: 637–653.
100. Dering M, Sękiewicz K, Boratyńska K, Litkowiec M, Iszkuło G, Romo A, Boratyński A. Genetic diversity and inter-specific relations of western Mediterranean relic Abies taxa as compared to the Iberian A. alba. Flora—Morphology, Distribution, Functional Ecology of Plants 2014; 209: 367–374.
101. Sękiewicz K, Dering M, Sękiewicz M, Boratyńska K, Iszkuło G, Litkowiec M., Ok T, Dagher-Kharrat MB, Boratyński A. Effect of geographic range discontinuity on species differentiation—East-Mediterranean Abies cilicica: a case study. Tree Genetics and Genomes 2015; 11: 810.
102. Eliades N-GH, Fady B, Gailing O, Leinemann L, Finkeldey R. Significant patterns of fine-scale spatial genetic structure in a narrow endemic wind-dispersed tree species, Cedrus brevifolia Henry. Tree Genetics & Genomes 2018; 14: 15.
103. Gómez A, Lunt DH. Refugia within refugia: patterns of phylogeographic concordance in the Iberian Peninsula. In: Weiss S, Ferrand N, eds. Phylogeography of southern European refugia. Springer, Dordrecht. 2007; 155–188.
104. Stewart JR, Lister AM, Barnes I, Dalén L. Refugia revisited: individualistic responses of species in space and time. Proceedings of the Royal Society of London B: Biological Sciences 2010; 277: 661–671.
105. Sękiewicz K, Dering M, Romo A, Dagher-Kharrat MB, Boratyńska K, Ok T, Boratyński A. Phylogenetic and biogeographic insights into long-lived Mediterranean Cupressus taxa with a schizo-endemic distribution and Tertiary origin. Botanical Journal of the Linnean Society 2018; 188: 190–212.
106. Médail F, Baumel A. Using phylogeography to define conservation priorities: The case of narrow endemic plants in the Mediterranean Basin hotspot. Biological Conservation 2018; 224: 258–266.
107. Torres-Díaz C, Ruiz E, González F, Fuentes G, Cavieres LA. Genetic Diversity in Nothofagus alessandrii (Fagaceae), an Endangered Endemic Tree Species of the Coastal Maulino Forest of Central Chile. Annals of Botany 2007; 100: 75–82. doi: 10.1093/aob/mcm073 17513870
108. Aleksić JM, Geburek T. Quaternary population dynamics of an endemic conifer, Picea omorika, and their conservation implications. Conservation Genetics 2014; 15: 87–107.
109. Dodd RS, DeSilva R. Long-term demographic decline and late glacial divergence in a Californian paleoendemic: Sequoiadendron giganteum (giant sequoia). Ecology and Evolution 2016; 6: 3342–3355. doi: 10.1002/ece3.2122 27252835
110. Millar MA, Byrne M, Coates DJ, Roberts JD. Contrasting diversity and demographic signals in sympatric narrow-range endemic shrubs of the south-west Western Australian semi-arid zone. Biological Journal of the Linnean Society 2016; 118: 315–329.
111. Hu Y, Dang M, Feng X, Woeste K, Zhao P. Genetic diversity and population structure in the narrow endemic Chinese walnut Juglans hopeiensis Hu: implications for conservation. Tree Genetics & Genomes 2017; 13: 91.
112. Jian H, Li S, Guo J, Li S, Wang Q, Yan H, Qiu, X, Zhang Y, Cai Z, Volis S, Tang K. High genetic diversity and differentiation of an extremely narrowly distributed and critically endangered decaploid rose (Rosa praelucens): implications for its conservation. Conservation Genetics 2018; 19: 761–776.
113. Lázaro-Nogal A, Matesanz S, García-Fernández A, Traveset A, Valladares F. Population size, center–periphery, and seed dispersers’ effects on the genetic diversity and population structure of the Mediterranean relict shrub Cneorum tricoccon. Ecology and Evolution 2017; 7: 7231–7242. doi: 10.1002/ece3.2940 28944013
114. Mitchell RJ, Karron JD, Holmquist KG, Bell JM. The influence of Mimulus ringens floral display size on pollinator visitation patterns. Functional Ecology 2004; 18: 116–124.
115. Hampe A, Petit RJ. Conserving biodiversity under climate change: the rear edge matters. Ecology Letters 2005; 8: 461–467. doi: 10.1111/j.1461-0248.2005.00739.x 21352449
116. Dering M, Rączka G, Szmyt J. Sex-specific pattern of spatial genetic structure in dioecious and clonal tree species, Populus alba L. Tree Genetics & Genomes 2016; 12: 70.
117. Garza JC, Williamson EG. Detection of reduction in population size using data from microsatellite loci. Molecular ecology 2001; 10: 305–318 doi: 10.1046/j.1365-294x.2001.01190.x 11298947
118. Lowe WH, Kovach RP, Allendorf FW. Population Genetics and Demography Unite Ecology and Evolution. Trends in Ecology & Evolution 2017; 32: 141–152.
119. Jump AS, Peñuelas J. Genetic effects of chronic habitat fragmentation in a wind-pollinated tree. Proceedings of the National Academy of Sciences 2006; 103: 8096–8100.
120. Karhu A, Vogl C, Moran GF, Bell JC, Savolainen O. Analysis of microsatellite variation in Pinus radiata reveals effects of genetic drift but no recent bottlenecks. Journal of Evolutionary Biology 2006; 19: 167–175. doi: 10.1111/j.1420-9101.2005.00982.x 16405588
121. Broquet T, Angelone S, Jaquiery J, Joly P, Lena J-P, Lengagne T, Plenet S, Luquet E, Perrin N. Genetic Bottlenecks Driven by Population Disconnection. Conservation Biology 2010; 24: 1596–1605. doi: 10.1111/j.1523-1739.2010.01556.x 20666803
122. Setoguchi H, Mitsui Y, Ikeda H, Nomura N, Tamura A. Genetic structure of the critically endangered plant Tricyrtis ishiiana (Convallariaceae) in relict populations of Japan. Conservation Genetics 2011; 12: 491–501.
123. Manchester SR, Chen Z-D, Lu A-M, Uemura K. Eastern Asian endemic seed plant genera and their paleogeographic history throughout the Northern Hemisphere. Journal of Systematics and Evolution 2009; 47: 1–42.
124. Young A, Boyle T, Brown T. The population genetic consequences of habitat fragmentation for plants. Trends in Ecology & Evolution 1996; 11: 413–418
125. Pluess AR, Stöcklin J. Genetic diversity and fitness in Scabiosa columbaria in the Swiss Jura in relation to population size. Conservation Genetics 2004; 5: 145–156.
126. Jackson ND, Fahrig L. Habitat amount, not habitat configuration, best predicts population genetic structure in fragmented landscapes. Landscape Ecology 2016; 31: 951–968.
127. Fady-Welterlen B. Is there really more biodiversity in Mediterranean forest ecosystems? Taxon 2005; 54: 905–910.
128. Piotti A, Leonarduzzi C, Postolache D, Bagnoli F, Spanu I, Brousseau L, Urbinati C, Leonardi S, Vendramin GG. Unexpected scenarios from Mediterranean refugial areas: disentangling complex demographic dynamics along the Apennine distribution of silver fir. Journal of Biogeography 2017; 44: 1547–1558.
129. Dubreuil M, Riba M, González‐Martínez SC, Vendramin GG, Sebastiani F, Mayol M. Genetic effects of chronic habitat fragmentation revisited: Strong genetic structure in a temperate tree, Taxus baccata (Taxaceae), with great dispersal capability. American Journal of Botany 2010; 97: 303–310. doi: 10.3732/ajb.0900148 21622391
130. González-Martínez SC, Dubreuil M, Riba M, Vendramin GG, Sebastiani F, Mayol M. Spatial genetic structure of Taxus baccata L. in the western Mediterranean Basin: Past and present limits to gene movement over a broad geographic scale. Molecular Phylogenetics and Evolution 2010; 55: 805–815. doi: 10.1016/j.ympev.2010.03.001 20211747
131. Mayol M, Riba M, González-Martínez SC, Bagnoli F, Beaulieu J-L, Berganzo E, Burgarella C, Dubreuil M, Krajmerová D, Paule L, Romšáková I, Vettori C, Vincenot L, Vendramin GG. 2015. Adapting through glacial cycles: insights from a long-lived tree (Taxus baccata). New Phytologist 208: 973–986. doi: 10.1111/nph.13496 26096330
132. García D, Ramón Obeso J. Facilitation by herbivore-mediated nurse plants in a threatened tree, Taxus baccata: local effects and landscape level consistency. Ecography 2003; 26: 739–750.
133. Devaney JL, Jansen MAK, Whelan PM. Spatial patterns of natural regeneration in stands of English yew (Taxus baccata L.); Negative neighbourhood effects. Forest Ecology and Management 2014; 321: 52–60.
134. Iszkuło G, Pers-Kamczyc E, Nalepka D, Rabska M, Walas Ł, Dering M. Postglacial migration dynamics helps to explain current scattered distribution of Taxus baccata. Dendrobiology 2016; 76: 81–89.
135. Bujoczek L, Bujoczek M. The dynamics of the Taxus baccata L. population and the factors affecting its regeneration in the Jasień Nature Reserve. Dendrobiology 2018; 80: 24–36.
136. Hampe A, Arroyo J. Recruitment and regeneration in populations of an endangered South Iberian Tertiary relict tree. Biological Conservation 2002; 107: 263–271.
137. Lázaro A, Traveset A, Castillo A. Spatial concordance at a regional scale in the regeneration process of a circum-Mediterranean relict (Buxus balearica): connecting seed dispersal to seedling establishment. Ecography 2006; 29: 683–696.
138. Baali-Cherif D, Besnard G. High genetic diversity and clonal growth in relict populations of Olea europaea subsp. laperrinei (Oleaceae) from Hoggar, Algeria. Annals of Botany 2005; 96: 823–830. doi: 10.1093/aob/mci232 16043438
139. Neophytou C, Konnert M, Fussi B. Western and eastern post-glacial migration pathways shape the genetic structure of sycamore maple (Acer pseudoplatanus L.) in Germany. Forest Ecology and Management 2019; 432: 83–93.
140. Skourlis K, Doutsos T. The Pindos Fold-and-thrust belt (Greece): inversion kinematics of a passive continental margin. International Journal of Earth Sciences, 2003; 92: 891–903.
141. Robledo-Arnuncio JJ, Collada C, Alia R, Gil L. Genetic structure of montane isolates of Pinus sylvestris L. in a Mediterranean refugial area. Journal of Biogeography 2005; 32: 595–605.
142. Ohsawa T, Ide Y. Global patterns of genetic variation in plant species along vertical and horizontal gradients on mountains. Global Ecology and Biogeography 2008; 17: 152–163.
143. Breed MF, Ottewell KM, Gardner MG, Marklund MH, Dormontt EE, Lowe AJ. Mating patterns and pollinator mobility are critical traits in forest fragmentation genetics. Heredity 2015; 115: 108. doi: 10.1038/hdy.2013.48 24002239
144. Ottewell K, Grey E, Castillo F, Karubian J. The pollen dispersal kernel and mating system of an insect-pollinated tropical palm, Oenocarpus bataua. Heredity 2012; 109: 332. doi: 10.1038/hdy.2012.40 22892637
145. Hoshizaki K. Rodent seed hoarding and regeneration of Aesculus turbinata: patterns, processes and implications. In: Tamura T., Sakio H. Ecology of Riparian Forests in Japan 2008; pp. 107–122. Springer, Tokyo.
146. Hoshizaki K, Suzuki W, Sasaki S. Impacts of secondary seed dispersal and herbivory on seedling survival in Aesculus turbinata. Journal of Vegetation Science 1997; 8: 735–742.
147. Magri D, Di Rita F, Aranbarri J, Fletcher W, González-Sampériz P. Quaternary disappearance of tree taxa from Southern Europe: Timing and trends. Quaternary Science Reviews 2017; 163: 23–55.
148. Gobet E, Schwörer Ch, van Leeuwen J, Wahab SA, Nielsen EH, Kotova N, Makhortykh S, Kiosak D, Tinner W. Vegetation shifts at the monumental Ukrainian site of Kamyana Mohyla during the Neolithisation period. In: Makhorytkh S., Capitani A., eds. Archaeology and paleoecology of the Ukrainian steppe. IA NAS, Kyiv, 2017; 51–65.
149. Van Zeist W, Woldring H. A postglacial pollen diagram from Lake Van in east Anatolia. Review of Palaeobotany and Palynology 1978; 26: 249–276
150. Feliner GN. Southern European glacial refugia: a tale of tales. Taxon 2011; 60: 365–372.
151. Isagi Y, Saito D, Kawaguchi H, Tateno R, Watanabe S. Effective pollen dispersal is enhanced by the genetic structure of an Aesculus turbinata population. Journal of Ecology 2007; 95: 983–990.
152. Casazza G, Giordani P, Benesperi R, Foggi B, Viciani D, Filigheddu R, Farris E, Bagella S, Pisanu S, Mariotti MG. Climate change hastens the urgency of conservation for range-restricted plant species in the central-northern Mediterranean region. Biological Conservation 2014; 179: 129–138.
153. Tang CQ, Dong Y-F, Herrando-Moraira S, Matsui T, Ohashi H, He L-Y, Nakao K. Tanaka N, Tomita M, Li X-S, Yan H-Z, Peng M-C, Hu J, Yang R-H, Li W-J, Yan K, Hou X, Zhang Z-Y, López-Pujol J. Potential effects of climate change on geographic distribution of the Tertiary relict tree species Davidia involucrata in China. Scientific Reports 2017; 7: 43822. doi: 10.1038/srep43822 28272437
154. Joffre R, Rambal S, Ratte JP. The dehesa system of southern Spain and Portugal as a natural ecosystem mimic. Agroforestry Systems 1999; 45: 57–79.
155. Pulido F, Valladares F, Calleja JA, Moreno G, González‐Bornay G. Tertiary relict trees in a Mediterranean climate: abiotic constraints on the persistence of Prunus lusitanica at the eroding edge of its range. Journal of Biogeography 2008; 35: 1425–1435.
156. Allen CD, Macalady AK, Chenchouni H, Bachelet D, McDowell N, Vennetier M, Kitzberger T, Rigling A, Breshears DD, Hogg EHT, Gonzalez P, Fensham R, Zhang Z, Castro J, Demidova N, Lim J-H, Allard G, Running SW, Semerci A, Cobb N. A global overview of drought and heat-induced tree mortality reveals emerging climate change risks for forests. Forest Ecology and Management 2010; 259: 660–684.
157. Giorgi F, Lionello P. Climate change projections for the Mediterranean region. Global and planetary change, 2008; 63: 90–104.
158. Smith DM, Scaife AA, Boer GJ, Caian M, Doblas-Reyes FJ, Guemas V, Hawkins E, Hazeleger W, Hermanson L, Ho CK, Ishii M, Kharin V, Kimoto M, Kirtman B, Lean J, Matei D, Merryfield WJ, Müller WA, Pohlmann H, Rosati A, Wouters B, Wyser K. Real-time multi-model decadal climate predictions. Climate Dynamics 2013; 41: 2875–2888.
159. Salazar E, Hammerling D, Wang X, Sansó B, Finley AO, Mearns LO. Observation-based blended projections from ensembles of regional climate models. Climatic Change 2016; 138: 55–69.
160. Phitos D, Constantinidis T, Kamari G. The Red Data Book of Rare and Threatened Plants of Greece. Hellenic Botanical Society, Patras. 2009.
Článek vyšel v časopise
PLOS One
2019 Číslo 12
- Řekněte mi, co jíte, a já vám řeknu, jaká je úroveň vaší kognice
- Původní český test ALBA pro odhalení mírných kognitivních deficitů získal certifikaci MZ ČR
- Duševní onemocnění se u teenagerů přenáší v rámci školní třídy
- Od kostiček k lidské kůži. Jak si postavit 3D biotiskárnu?
- Prof. Libor Grubhoffer: Nová varianta mpox je příležitostí pokročit v poznání původce varioly
Nejčtenější v tomto čísle
- Methylsulfonylmethane increases osteogenesis and regulates the mineralization of the matrix by transglutaminase 2 in SHED cells
- Oregano powder reduces Streptococcus and increases SCFA concentration in a mixed bacterial culture assay
- The characteristic of patulous eustachian tube patients diagnosed by the JOS diagnostic criteria
- Parametric CAD modeling for open source scientific hardware: Comparing OpenSCAD and FreeCAD Python scripts