Next Article in Journal
Could Bryophagous Beetles (Coleoptera: Byrrhidae) Help Us Understand Bryophyte Taxonomy? Preferences within the Hypnum cupressiforme Hedw. Species Complex
Next Article in Special Issue
Phytogeographical Analysis and Ecological Factors of the Distribution of Orchidaceae Taxa in the Western Carpathians (Local study)
Previous Article in Journal
Accumulation of Potentially Toxic Elements in Mosses Collected in the Republic of Moldova
Previous Article in Special Issue
Next Chapter in the Legend of Silphion: Preliminary Morphological, Chemical, Biological and Pharmacological Evaluations, Initial Conservation Studies, and Reassessment of the Regional Extinction Event
 
 
Search for Articles:
Title / Keyword
Author / Affiliation / Email
Journal
Article Type
 
 
Advanced Search
 
Section
Special Issue
Volume
Issue
Number
Page
 
Journals
Plants
Volume 10
Issue 3
10.3390/plants10030470
plants-logo
Submit to this Journal Review for this Journal Propose a Special Issue
► Article Menu

Article Menu

share Share announcement Help format_quote Cite question_answer Discuss in SciProfiles thumb_up
...
Endorse textsms
...
Comment
first_page
settings
Order Article Reprints
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

An Orchid in Retrograde: Climate-Driven Range Shift Patterns of Ophrys helenae in Greece

by
Martha Charitonidou
1,*,
Konstantinos Kougioumoutzis
2,* and
John M. Halley
1,*
1
Laboratory of Ecology, Department of Biological Applications and Technology, University of Ioannina, 45110 Ioannina, Greece
2
Department of Ecology and Systematics, Faculty of Biology, National and Kapodistrian University of Athens, 15784 Athens, Greece
*
Authors to whom correspondence should be addressed.
Plants 2021, 10(3), 470; https://doi.org/10.3390/plants10030470
Submission received: 31 January 2021 / Revised: 23 February 2021 / Accepted: 25 February 2021 / Published: 2 March 2021
(This article belongs to the Special Issue Taxonomy and Plant Conservation)
Browse Figures
Versions Notes

Abstract

:
Climate change is regarded as one of the most important threats to plants. Already species around the globe are showing considerable latitudinal and altitudinal shifts. Helen’s bee orchid (Ophrys helenae), a Balkan endemic with a distribution center in northwestern Greece, is reported to be expanding east and southwards. Since this southeastern movement goes against the usual expectations, we investigated via Species Distribution Modelling, whether this pattern is consistent with projections based on the species’ response to climate change. We predicted the species’ future distribution based on three different climate models in two climate scenarios. We also explored the species’ potential distribution during the Last Interglacial and the Last Glacial Maximum. O. helenae is projected to shift mainly southeast and experience considerable area changes. The species is expected to become extinct in the core of its current distribution, but to establish a strong presence in the mid- and high-altitude areas of the Central Peloponnese, a region that could have provided shelter in previous climatic extremes.

1. Introduction

The orchid family is among the earliest divergent angiosperms [1,2] and as such, have a nearly cosmopolitan distribution [3] and comprise more than 28,000 species [4], thus being one of the largest and most diverse plant families in the world. Known for their astonishing morphologies, the complexity of their life cycles, and their high levels of speciation [5,6,7], orchids have always been in the spotlight of scientists, being characterized as a great model for ecological studies to target conservation priorities [8]. Like other plants, orchids are expected to be responding to climate change (cf. [8,9,10,11,12,13]). However, for a number of reasons, such as their reproduction complexity and their dependence on other organisms during their life history, orchids are thought to be more vulnerable [13,14].
The last few decades, the trend of global warming is more evident, and the effects on living organisms can be clearly seen [15,16,17,18], with changes in the phenology and distribution of species, as well as possible extinctions attributed to the rapidly changing climate [15,16,19,20]. Many studies have focused on the range changes of species, demonstrating a poleward movement, primarily by the leading edge expanding mainly to the north. Also observed are changes along altitudinal gradients, with species expanding into higher elevations [17,19,21]. Parmesan and Yohe [16] by analyzing ca. 1700 taxa, found that species are facing significant range shifts, with an average poleward movement of 6.1 km per decade. Chen et al. [21] found almost a three times larger shift, with a median rate of a northward shift of 16.9 km per decade, and an uphill shift of 11.0 meters per decade. It is thus clear, that the expected pattern of range changes for plants is a poleward or uphill migration. Although the majority of observations follow this general trend, many have been inconclusive and there is a percentage of species (mainly forest plant taxa) that moved to the opposite direction, southwards and/or downhill, showing a wide variety of range shifts [19,21].
Orchids are no exception when it comes to climate change. Like all plants, Orchidaceae are affected by changes in climatic conditions worldwide, and are reported to follow the general poleward trends and uphill patterns of range shifts. In addition, orchids could face especially significant losses to their distributions, as well as population declines [22,23,24] because of the complexity of their life cycle and their symbiotic dependence on other organisms. In complex symbioses, the responses of all partners to climate change need to be congruous. Thus, any adaptive changes plants make that affect pollination services or the mycorrhizal symbiosis (e.g., [8,22,23]) must be consonant with changes made by their partners. Hence, orchids could face greater risks from global warming [14,25] than other families. This is even more alarming taking into consideration that ca. 40% of all plants are facing greater risk [26] and Orchidaceae are one of the families that are under-assessed regarding their extinction risk status [27].
Greece has been extensively floristically explored during the last two centuries [28,29]. Consequently, many studies exist dealing with the factors affecting the biogeographical and biodiversity patterns in Greece ([30] and references therein). Nevertheless, very few studies have taken into consideration the impacts of climate change on plant distribution patterns [31,32,33,34,35] and none has yet assessed if the projected future distribution shifts constitute a physiological response previously experienced by the given taxon in question. After all, incorporating past climate projections in ecological niche modelling may further our understanding of species responses to climate change [36], and could lead to the development of efficient conservation schemes [37].
Being a Mediterranean country with a high topographical and environmental heterogeneity, Greece stands out for its species-rich flora (> ca. 7,500 taxa, 20.4% endemic; [38,39]). One of the most notable aspects of the Greek floral diversity is its orchid species richness (193 orchid species and subspecies; [40]). Orchids can be found nearly everywhere in Greece; 85.4% of the country’s surface hosts orchid species, that are spread across different habitat types, from phrygana and thermophilous pine forests (e.g., several Anacamptis, Ophrys, Serapias taxa occur in these habitats), to fir forests and alpine grasslands (e.g., Cephalanthera, Dactylorhiza, Epipactis) [41]. Most of the orchid species that occur in Greece belong to the genus Ophrys (bee orchids, 91 taxa according to [40]), which has its distribution and speciation center in the eastern Mediterranean [40,42]. Ophrys is an ideal example of plant-pollinator coevolution, as the majority of Ophrys species use sexual deception as a pollination mechanism. By encompassing visual and olfactory signals, they deceive naïve male Hymenoptera to pseudocopulation, and through this procedure, they achieve pollination, and thus, reproduction [43]. Helen’s bee orchid (Ophrys helenae Renz) is the only known exception to this rule for Ophrys. This Balkan endemic relies not on sexual deception for pollination, but on shelter mimicry [43,44].
Recently, Tsiftsis and Tsiripidis [45] assessed the patterns of Greek orchid distributions in time and space, and found that the majority of them, has shown a range increase. Among others, Ophrys helenae seems to follow this pattern (see Table 3 in [45]). Described in 1928 on Corfu Island by Dr. Jany Renz [46], its known distribution until the 1980s was restricted to Epirus (NW Greece), and Corfu [47,48]. However, the last few decades, its distribution appears to be expanding, with populations being found in east and southern Greece [40,45]. This observation is intriguing, which seems to run counter to the patterns that are seen due to climate change. Parmesan and Hanley [17] have noted that it is these exceptions to rules have often been the main drivers of advance in ecological understanding.
The purpose of this paper is to investigate whether the conjecture, expressed in [40] for O. helenae, that “…its distribution seems to be undergoing an active expansion towards the east and southern Greece…”, is consistent with a more systematic analysis, and how climate change will affect the species’ range in the upcoming decades. By using Species Distribution Models (SDMs), we show that this initially puzzling behavior is in fact what we would expect. Moreover, we show that the southward movement of O. helenae materializes through an interplay of climatic change with other factors (e.g., topography, soil characteristics, bioclimatic variables, land-cover), and that the species has probably experienced this pattern of shift in the past.

2. Results

The ensemble of small models (ESM) framework predictions was very good, with sufficient predictive power (TSS ≥ 0.95 for all algorithms and the ensemble prediction for both thinning procedures and climate databases, Figure S1). The intrathinning procedure and -climate database variation were statistically insignificant (Kruskal-Wallis ANOVA: H = 7.21, d.f. = 3, p = 0.06).
Precipitation of the wettest month (PWM) had the highest contribution among the response variables, followed by the potential evapotranspiration of the wettest quarter (PETWQ) for almost all thinning procedures and climate databases. The only exception was the CHELSA – Geographical Thinning combination (CHGEO), where the predicted mean value of soil organic carbon mass fraction at standard depth of 5 cm (SOC-5) was the second most important variable (Table 1).
The resulting habitat suitability maps (Figure 1 and Figures S2–S4) had high bioclimatic consistency for every combination of the thinning procedures and distribution areas (Figure 2). They were converted into binary maps, and then compared to the binary maps obtained for each Global Circulation Model (GCM), Representative Concentration Pathway (RCP) scenario, time-period, thinning procedure and climate database.
As the trends for the future potential distribution of Ophrys helenae were largely identical across all uncertainty sources, we selected to present the area-range change for one combination of climatic database and thinning procedure, the CHELSA – Geographical Thinning (CHGEO). Our results indicate that the response to climate is expected to contain a number of complexities; Ophrys helenae has experienced range expansion and contraction shifts due to periodic climate oscillations (Figure 3 and Figures S5–S7).
The extent of suitable areas attained was largest during the LIG under any thinning procedure/climate database (Figure 3 and Figures S5–S7). Climatically suitable areas for the species underwent slight contraction during the LGM compared to the LIG and moderate contraction in the current period, relative to the LGM. Furthermore, Ophrys helenae is expected to face considerable future area changes for all selected GCMs and RCPs, for the presented combination of climatic database and thinning procedure (CHGEO: −65.63% to 12.22% depending on the GCM and RCP combination; Table 2 and Figure 3) and all other combinations (Table S1 and Figures S5–S7). In all cases, there is a significant loss of area (Table 2 and Table S1), especially in NW Greece, near the current centre of its range (Figure 3 and Figures S5–S7). However, the species is also projected to expand its southern limits in the near future, shifting its distribution southwards to the Peloponnese, thus resembling its LIG extent of occurrence (Table 2 and Table S1, Figure 3 and Figures S5–S7). The centroids of all projected future distributions of Ophrys helenae appear to be lying to the south and southeast of the current distribution centroid, validating the southward shift of the species for future projections (Figure 4), with the sole exception of the HadGEM2 RCP 2.6 projection. The results are consistent through all combinations of climatic databases and occurrence data thinning procedures (Figures S8–S10).
We detected a statistically significant altitudinal shift between the different time-slices included in our analysis for all climate databases and thinning procedures (Table 3 and Table S2; Kruskal-Wallis ANOVA: H = 14345, d.f. = 8, p < 0.01). Ophrys helenae currently occurs in significantly lower altitudes compared to either its past (i.e., in the LIG and the LGM) or its future distribution, with the HadGEM2 RCP8.5 GCM/RCP combination presenting the highest mean altitude among all GCMs, RCPs and time-slices (mean altitude: 915 m a.s.l.; Table 3 and Table S2).

3. Discussion

Climate change is expected to continue and to intensify in the future, with a warming effect in a global scale. Plants have been affected by climatic changes in factors that play a major role in their life cycle (e.g., temperature, precipitation) [15,20], resulting in phenological alterations [18,49], as well as in range shifts. This anticipated warming is expected to drive plant populations to shift their distributions on both latitudinal and altitudinal gradients or even become extinct [50,51,52,53]. Species are moving poleward in both hemispheres, while on mountains the altitudinal ranges of species are expected to shift to greater elevation (cf. [16,17,21]). All of these patterns have been observed in numerous studies (e.g., [22,23,25,28,41,42]). However, there are always exceptions to the rule; for example, although most of the shifts observed have been northwards and uphill, there are few cases where plants followed in the opposite direction (see [19,21]). This has led to an increase in biotic homogenization all over the globe, since many range-restricted species are experiencing range contractions, while widespread and alien species are gaining ground, due to the intensifying effects of both climate and land-use change [54,55,56,57,58,59,60].
Orchids, like other plant species, are expected to be affected by the pressures of the changing climate [14]. A number of papers have been written about climate change and orchids, with the majority of them focusing on the physiological effects of changing climate, especially the relationships of orchids with their symbionts (mycorrhizal fungi and pollinators) [8,22,23,61]. The available literature about shifts of orchid distributions (e.g., [62,63,64,65]) has shown that the majority of studied species will be following the usual patterns of distribution change; decreasing of suitable habitat areas of orchids, and northwards or uphill shifts of distribution. Although orchid species that have been studied so far under the prism of climate change are following a distribution shift to higher latitudes and elevations, contrary migration is not an unusual phenomenon for plants and other organisms. Chen et al. [21] found that despite the significant movement northwards and uphill for the majority of species, there has been some movement in the opposite direction, with 22% of the studied taxa shifting towards the equator in latitude, and 25% moving downhill for elevation. Lenoir et al. [19] also found a significant minority going against the mainstream, with ca. 31% of the studied forest plants migrating downhill in W. Europe. The degree of fidelity to the mainstream trend will depend on a variety of functional traits such as generation time, seed mass [66] or, in the case of orchids especially, patterns of life history and fecundity [67].
We found that the expected movement of Helen’s bee orchid is anadromous based on its latitudinal shift; for most of the climatic models and scenarios that we used, the species is moving against the mainstream pattern, shifting its distribution southwards in the coming decades (Figure 4). On the other hand, the species is projected to move uphill in the future, as its mean altitudinal difference reaches ca. 300 m (Table 3). Thus, Ophrys helenae seems to be a cold-adapted species based on its response to climate change [16,19,21], as indicated by its distributional and altitudinal shifts from the LIG up to the current time-period, a trend observed for other range-restricted plant species as well [68]. Its slightly lower mean altitude currently observed is directly attributed to the extreme grazing pressure the species is facing in its core distributional area (i.e., Epirus in NW Greece) [69,70,71], that has probably forced many of its populations to extinction. There is a moderate distribution change variation among the GCMs and RCPs that we included in our analysis, but in all cases movement is projected to be mostly south-eastwards. This is in agreement with the pattern observed in [40], but in contrast to the expected poleward trend of most orchids noted in the review of [22]. A northwards range shift was observed by [72] for the lizard orchid (Himantoglossum hircinum) in the UK, as well as from [73] for the lady’s slipper orchid (Cypripedium calceolus) in Europe.
Ophrys helenae is predicted to experience moderate to high changes in range (Table 2, Figure 3 and Figure 4; see also Supplementary Material), retreating from its northern and western edges, while expanding mainly in the central Peloponnese. Nevertheless, the extreme changes associated with RCP 8.5 for all GCMs will be unequivocally negative for O. helenae distribution in Greece, since at least a 20% overall range contraction is predicted to occur. This aligns with the majority of studies for changes in distribution of orchid species. For instance, Ophrys argolica and O. delphinensis, two orchids occurring in Greece are projected to face a severe loss of their distribution areas, with the latter projected to become extinct [74]. The same holds true for all endemic Cretan orchids [32], which are expected to face severe mean losses of their distribution areas (min: 73.6% - Epipactis cretica; max: 99.3% - Ophrys omegaifera subsp. fleischmannii). In addition, according to [73], Cypripedium calceolus will face a major decline of its niche in the future (30–63% loss) under all studied scenarios, the same as for Epipactis helleborine, that is expected to undergo significant losses throughout its whole distribution range, despite its wide extent [62]. On the other hand, in the less pessimistic scenario (RCP 2.6), the overall range of Helen’s bee orchid shows a less contraction and may even expand (ca. 15%, CCSM4 model). Similar results were found by [72] for Himantoglossum hircinum, which is expected to gain from changing conditions due to climate change, and significantly widen its distribution in the UK northwards. These findings are also in line with those of [75] for Sardinian orchids, such as Anacamptis papilionacea var. papilionacea, Serapias parviflora, Ophrys bombyliflora and O. morisii, that are found to significantly expand their distributions in the future under more optimistic RCPs.
The patterns of change expected for O. helenae are complex and not easily contained within a simple description. Lenoir et al. [19] wrote that “…climate warming does not only affect species at their range boundaries, but its consequences ripple through the whole range of species.”, That is, climatic change will affect species throughout their range and not just at the boundaries. The simple model of a linearly shifting plant distribution in response to the climatic signal is simplified, just as the notion of plant distributions thin at their edges and aggregating towards an “abundant center” are also simplified [76]. This is clearly demonstrated here where the species is always projected to be extirpated in its current distribution center, and not just at its boundaries.
In this study, we explored Ophrys helenae distributions in the past, for two time slices, the Last Glacial Maximum (LGM) and the Last Interglacial (LIG), in order to investigate whether the projected future changes of O. helenae are following a recurring pattern. It is notable that a significant portion of the species’ niche is in the mid- to high-altitude areas of central Peloponnese. Thus, when comparing the projected distributions in the prehistoric past, we see that the expected gains in the future distributions, especially south in the central Peloponnese, are returning Helen’s bee orchid to its former “heartlands”. The species will reappear in areas where it should have been at the time of the Last Interglacial, when the climatic conditions were more or less similar to the ones expected in a warmer future [77,78].

4. Materials and Methods

4.1. Study Species

Helen’s bee Orchid (Ophrys helenae Renz – sensu [38]) is a Balkan endemic orchid, with its distribution centre lying in northwestern Greece (Figure 5), where it is locally common, while it is also occurring in Albania. A perennial species, with an average height of 26.5 cm, O. helenae is one of the most easily recognisable bee orchids, due to its large, cherry-red labellum that lacks a speculum design. In addition, it is the only Ophrys species that enlists shelter mimicry as its pollination strategy, rather than sexual deception – the dominant one for bee orchids [43,44]. It can be found in bloom from late March through April and May, in full sun and/or semishade sites among shrubs and forest openings. Recent confirmed observations including populations from East Attica and the Peloponnese reveal an active species’ distribution, with an eastward and southward expansion [40].

4.2. Species Occurrence Data

Ophrys helenae occurrence data (721 occurrences) were obtained from the database of the Orchid Flora of Greece project ([40] and pers. comm. with Asst. Prof. Spyros Tsiftsis).
We environmentally and geographically thinned the species’ occurrences to avoid pseudoreplication and reduce sampling bias [79,80]. The geographical data cleaning and organizing procedure followed the procedure outlined in [81], using functions from the ‘biogeo’ 1.0 [81] and ‘spThin’ 0.1.0 [80] R packages. The environmental filtering procedure can improve model performance [82] and was based on the representative and uncorrelated environmental variables occurring in the study area (see environmental data below) following [82]. Finally, we evaluated whether any geographical sampling bias existed in our species occurrence data by comparing the statistical distance distribution observed in our dataset to a simulated distribution expected under random sampling via the ‘sampbias’ 1.0.4 [83] R package.

4.3. Environmental Data

Current, future and past (Last Glacial Maximum (LGM) and Last Interglacial (LIG)) climatic data were obtained from the WorldClim [84] and the CHELSA [85] databases at a 30 sec resolution (except for LGM data: 2.5 arc-min resolution). We were thus able to assess the bioclimatic consistency and congruence of our models [86] and estimate climate-database uncertainty [87].
We obtained sixteen additional climatic variables at the same resolution via the ‘envirem’ 2.2 [88] R package based on the bioclimatic data from WorldClim and CHELSA for all time-slices. We selected three Global Circulation Models (GCMs) based on [89] and two different Intergovernmental Panel on Climate Change (IPCC) scenarios from the Representative Concentration Pathways (RCP) family: RCP2.6 (mild scenario) and RCP8.5 (severe scenario). We extracted soil variables from the SoilGrids 250 m database [90]. We extracted elevation data from the CGIAR-CSI data-portal [91], which were then aggregated and resampled using functions from the ‘raster’ 3.3.13 R package [92] in order to match the resolution of the other environmental variables. We estimated five topographical variables (slope, aspect, heat load index, topographic position index and terrain ruggedness index) using functions from the ‘raster’ 3.3.13 [92] and ‘spatialEco’ 1.2-0 R packages [93]. We also created a raster layer with the calcareous substrates in Greece, using functions from the ‘sf’ 0.9.6 [94] and the ‘fasterize’ 1.0.3 R package [95], based on the Geological Map of Greece [96].
From this initial set of 50 predictors, only twelve were not highly correlated (Spearman rank correlation < 0.7 and Variance Inflation Factor < 5 – [97]). Multi-collinearity assessment was performed using the ‘vifcor’ function from the ‘usdm’ 1.1.18 [98] R package.

4.4. Species Distribution Models

4.4.1. Model Parameterization and Evaluation

We modelled the realized climatic niche of Ophrys helenae by combining the available occurrence data with current environmental predictors with the ‘biomod2’ 3.4.6 [99] and ‘ecospat’ 3.1 [100] R packages. We used three different modelling algorithms for our study species: Random Forest (RF), Classification Tree Analysis (CTA) and Artificial Neural Networks (ANN) in an ensemble modelling scheme, to reduce model algorithm uncertainty [87,101]. We generated pseudoabsences following the recommendations of [102] at a minimum distance of 42.9 km from presence locations, which equals the median autocorrelation of the non-collinear environmental variables, using the ‘blockCV’ 2.1.1 [103] R package. We followed the ensemble of small models (ESM) framework [104,105,106], since the occurrence/predictors ratio was lower than 20 [79]. ESMs were subsequently calibrated by fitting bivariate models, which were then averaged into an ensemble model using weights based on model performances. For all models, prevalence was equal to 0.5. We used the True Skill Statistic (TSS; [107]) to evaluate the models’ predictive performance based on a repeated (10 times) 80-20 split-sampling approach. We used null model significance testing [108] to evaluate the performance of our model, which outperformed the null expectation at p < 0.001.

4.4.2. Model Projections

We used the calibrated models to predict the suitable area for Ophrys helenae under current, future and past conditions [101]. The contribution of each model to the ensemble forecast was weighted according to its TSS score. Models with a TSS score < 0.8 (poorly calibrated models) were excluded from building projections. The final models used for spatial projections were calibrated using 100% of the data. Models were binarised using the value maximizing the TSS score as the threshold for distinguishing presence and absence predictions. As a conservative approach, the suitability of all cells having non-zero values in the clamping mask was set to zero [79]. Finally, we applied a mask representing urban and suburban areas to eliminate any cells that are unsuitable regardless of the environmental conditions.

4.4.3. Area Range Change

Using functions from the ‘biomod2’ 3.4.6 [99] R package, we were able to investigate if Ophrys helenae will experience range reduction or increase under future and past conditions.

4.5. Bioclimatic Congruence and Consistency

We followed the framework of [86] in order to construct the bioclimatic congruence and consistency maps for Ophrys helenae for every time-period that was available in both climate databases.

4.6. Distribution Changes in Latitudinal and Altitudinal Gradient

We estimated the distribution centroids of all cases, in order to test if and how the distribution of Ophrys helenae may have shifted in all time-slices (present, past and future), and combinations of climate databases (CHELSA – WorldClim) and occurrence data thinning procedures (Geographical – Environmental), using the ‘st_centroid’ function of the ‘sf’ 0.9.6 [94] R-package. All generated points were plotted in maps via the QGIS v.3.14.16 ‘Pi’ [109]. Also, in order to test the species’ altitudinal shift, we compared the altitude the species is predicted to appear in all time-slices via Kruskal-Wallis tests.

Supplementary Materials

The following are available online at https://www.mdpi.com/2223-7747/10/3/470/s1, Figure S1. True Skill Statistic (TSS) values for every source of uncertainty included in the present study. THGEO: geographical thinning. THENV: environmental thinning. HUALPHA: distribution area estimated with the alpha-hull method. CH: CHELSA climate database. WC: WorldClim climate database. Panels A, B, and C correspond to the thinning procedure, climate database, and all the aforementioned uncertainty sources combined, respectively. Figure S2. Habitat suitability map for Ophrys helenae, for the combination of CHELSA climate database and environmental thinning procedure. Orange border delineates the species’ distributional area based on the Atlas of Greek Orchids. Figure S3. Habitat suitability map for Ophrys helenae, for the combination of WorldClim climate database and geographical thinning procedure. Orange border delineates the species’ distributional area based on the Atlas of Greek Orchids. Figure S4. Habitat suitability map for Ophrys helenae, for the combination of WorldClim climate database and environmental thinning procedure. Orange border delineates the species’ distributional area based on the Atlas of Greek Orchids. Figure S5. Past, present and future potential distribution maps for Ophrys helenae, based on the environmentally thinned occurrence data and climate data from CHELSA (CHENV) that show the transition from one time-slice to another. In the central panel, green coloring indicates the cells the species currently potentially occupies, while the orange lines delineate the species’ distributional area. Left-hand panel: past potential distribution maps for showing the transition from the Last Interglacial (ca. 120–140 Kya) to the Last Glacial Maximum (ca. 20 Kya) and from the Last Glacial Maximum (ca. 20 Kya) to the present time-period. Red, blue and grey coloring depict cells that show loss, gain, or remained stable, respectively, compared to current the previous conditions. Future potential distribution maps represent combinations of three Global Circulation Models (BCC, CCSM4, and HadGEM2-ES) and two Representative Concentration Pathways (RCP 2.6 and RCP 8.5) showing the transition from the present time-period to each respective GCM and RCP combination. The red grid cells in the right-hand panel indicate that the species is currently present at these areas but will not be in the future. The grey grid cells in the right-hand panel indicate that the species is not currently present at these areas and it will not be in the future. The blue grid cells in the right-hand panel indicate that the species currently occupies and will continue to occupy these areas in the future. Figure S6. Past, present and future potential distribution maps for Ophrys helenae, based on the geographically thinned occurrence data and climate data from WorldClim (WCGEO) that show the transition from one time-slice to another. In the central panel, green coloring indicates the cells the species currently potentially occupies, while the orange lines delineate the species’ distributional area. Left-hand panel: past potential distribution maps for showing the transition from the Last Interglacial (ca. 120–140 Kya) to the Last Glacial Maximum (ca. 20 Kya) and from the Last Glacial Maximum (ca. 20 Kya) to the present time-period. Red, blue and grey coloring depict cells that show loss, gain, or remained stable, respectively, compared to current the previous conditions. Future potential distribution maps represent combinations of three Global Circulation Models (BCC, CCSM4, and HadGEM2-ES) and two Representative Concentration Pathways (RCP 2.6 and RCP 8.5) showing the transition from the present time-period to each respective GCM and RCP combination. The red grid cells in the right-hand panel indicate that the species is currently present at these areas but will not be in the future. The grey grid cells in the right-hand panel indicate that the species is not currently present at these areas and it will not be in the future. The blue grid cells in the right-hand panel indicate that the species currently occupies and will continue to occupy these areas in the future. Figure S7. Past, present and future potential distribution maps for Ophrys helenae, based on the environmentally thinned occurrence data and climate data from WorldClim (WCENV) that show the transition from one time-slice to another. In the central panel, green coloring indicates the cells the species currently potentially occupies, while the orange lines delineate the species’ distributional area. Left-hand panel: past potential distribution maps for showing the transition from the Last Interglacial (ca. 120–140 Kya) to the Last Glacial Maximum (ca. 20 Kya) and from the Last Glacial Maximum (ca. 20 Kya) to the present time-period. Red, blue and grey coloring depict cells that show loss, gain, or remained stable, respectively, compared to current the previous conditions. Future potential distribution maps represent combinations of three Global Circulation Models (BCC, CCSM4, and HadGEM2-ES) and two Representative Concentration Pathways (RCP 2.6 and RCP 8.5) showing the transition from the present time-period to each respective GCM and RCP combination. The red grid cells in the right-hand panel indicate that the species is currently present at these areas but will not be in the future. The grey grid cells in the right-hand panel indicate that the species is not currently present at these areas and it will not be in the future. The blue grid cells in the right-hand panel indicate that the species currently occupies and will continue to occupy these areas in the future. Figure S8. Distribution centroids for past, present and future projected distributions of Ophrys helenae in Greece, for the combination of CHELSA climatic database, and environmental thinning procedure for occurrence data. Figure S9. Distribution centroids for past, present and future projected distributions of Ophrys helenae in Greece, for the combination of WorldClim climatic database, and geographical thinning procedure for occurrence data. Figure S10. Distribution centroids for past, present and future projected distributions of Ophrys helenae in Greece, for the combination of WorldClim climatic database, and environmental thinning procedure for occurrence data. Table S1. Percentages of range loss and gain, an overall range change for past and future projections of Ophrys helenae distribution. For the past, Last Interglacial and Last Glacial Maximum (LIG and LGM, respectively) are presented, while for the future, all selected Global Circulation Models (GCMs) and Representative Concentration Pathways (RCPs). BC: BCC-CSM-1, CC: CCSM4, HE: HadGEM2-ES. The presented values are for all climatic databases, and occurrence data thinning procedures. Table S2. Mean altitude for past, present and future projections of Ophrys helenae distribution. For the past, Last Interglacial and Last Glacial Maximum (LIG and LGM, respectively) are presented, while for the future, all selected Global Circulation Models (GCMs) and Representative Concentration Pathways (RCPs). BC: BCC-CSM-1, CC: CCSM4, HE: HadGEM2-ES. The presented values are for all climatic databases, and occurrence data thinning procedures.

Author Contributions

Conceptualization, M.C., K.K.; investigation, M.C. and K.K.; methodology, K.K. and M.C.; formal analysis, M.C. and K.K.; resources M.C.; supervision, J.M.H.; writing—original draft preparation, M.C., K.K.; writing—review and editing, M.C., K.K. and J.M.H.; visualization, M.C. and K.K. All authors have read and agreed to the published version of the manuscript.

Funding

This study is part of first author’s doctoral research, funded by the General Secretariat for Research and Technology (GSRT) and the Hellenic Foundation for Research and Innovation (HFRI) (Scholarship Code: 2359).

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

Not applicable.

Acknowledgments

We would like to thank Spyros Tsiftsis for providing the occurrences of Ophrys helenae from the database of the Orchid Flora of Greece project.

Conflicts of Interest

The authors declare no conflict of interest. The funders had no role in the design of the study, in the collection, analyses, or interpretation of data, in the writing of the manuscript, or in the decision to publish the results.

References

  1. Givnish, T.J.; Spalink, D.; Ames, M.; Lyon, S.P.; Hunter, S.J.; Zuluaga, A.; Doucette, A.; Caro, G.G.; McDaniel, J.; Clements, M.A.; et al. Orchid historical biogeography, diversification, Antarctica and the paradox of orchid dispersal. J. Biogeogr. 2016, 43, 1905–1916. [Google Scholar] [CrossRef]
  2. Poinar, G., Jr.; Rasmussen, F.N. Orchids from the past, with a new species in Baltic amber. Bot. J. Linn. Soc. 2017, 183, 327–333. [Google Scholar] [CrossRef]
  3. Chase, M.W. Classification of Orchidaceae in the Age of DNA data. Curtis’s Bot. Mag. Bot. Mag. 2005, 22, 2–7. [Google Scholar] [CrossRef]
  4. WCSP World Checklist of Selected Plant Families. Available online: http://wcsp.science.kew.org/ (accessed on 28 February 2020).
  5. Darwin, C. On the Various Contrivances by Which British and Foreign Orchids Are Fertilised by Insects: And on the Good Effect of Intercrossing; John Murray: London, UK, 1862; ISBN 9780511910197. [Google Scholar]
  6. Shefferson, R.P.; Jacquemyn, H.; Kull, T.; Hutchings, M.J. The demography of terrestrial orchids: Life history, population dynamics and conservation. Bot. J. Linn. Soc. 2020, 192, 315–332. [Google Scholar] [CrossRef]
  7. Bateman, R.M. Circumscribing species in the European orchid flora: Multiple datasets interpreted in the context of speciation mechanisms. Ber. aus den Arb. Heim. Orchid. Beih. 2012, 8, 160–212. [Google Scholar]
  8. Swarts, N.D.; Dixon, K.W. Terrestrial orchid conservation in the age of extinction. Ann. Bot. 2009, 104, 543–556. [Google Scholar] [CrossRef] [Green Version]
  9. Wraith, J.; Pickering, C. Quantifying anthropogenic threats to orchids using the IUCN Red List. Ambio 2018, 47, 307–317. [Google Scholar] [CrossRef] [Green Version]
  10. Wraith, J.; Norman, P.; Pickering, C. Orchid conservation and research: An analysis of gaps and priorities for globally Red Listed species. Ambio 2020, 49, 1601–1611. [Google Scholar] [CrossRef]
  11. Liu, H.; Feng, C.L.; Luo, Y.B.; Chen, B.S.; Wang, Z.S.; Gu, H.Y. Potential challenges of climate change to orchid conservation in a Wild Orchid Hotspot in Southwestern China. Bot. Rev. 2010, 76, 174–192. [Google Scholar] [CrossRef]
  12. Hágsater, E.; Dumont, V. Orchids: Status Survey and Conservation Action Plan; IUCN: Gland, Switzerland, 1996; ISBN 2831703255. [Google Scholar]
  13. Seaton, P.; Kendon, J.P.; Pritchard, H.W.; Murti Puspitaningtyas, D.; Marks, T.R. Orchid conservation: The next ten years. Lankesteriana 2013, 13, 93–101. [Google Scholar] [CrossRef] [Green Version]
  14. Fay, M.F. Orchid conservation: How can we meet the challenges in the twenty-first century? Bot. Stud. 2018, 59, 16. [Google Scholar] [CrossRef] [Green Version]
  15. Walther, G.-R.; Post, E.; Convey, P.; Menzel, A.; Parmesan, C.; Beebee, T.J.C.; Fromentin, J.-M.; Hoegh-Guldberg, O.; Bairlein, F. Ecological responses to recent climate change. Nature 2002, 416, 389–395. [Google Scholar] [CrossRef]
  16. Parmesan, C.; Yohe, G. A globally coherent fingerprint of climate change. Nature 2003, 421, 37–42. [Google Scholar] [CrossRef]
  17. Parmesan, C.; Hanley, M.E. Plants and climate change: Complexities and surprises. Ann. Bot. 2015, 116, 849–864. [Google Scholar] [CrossRef] [PubMed]
  18. Easterling, D.R.; Meehl, G.A.; Parmesan, C.; Changnon, S.A.; Karl, T.R.; Mearns, L.O. Climate Extremes: Observations, Modeling, and Impacts. Science 2000, 289, 2068–2074. [Google Scholar] [CrossRef] [Green Version]
  19. Lenoir, J.; Gégout, J.C.; Marquet, P.A.; De Ruffray, P.; Brisse, H. A significant upward shift in plant species optimum elevation during the 20th century. Science 2008, 320, 1768–1771. [Google Scholar] [CrossRef]
  20. Parmesan, C. Ecological and evolutionary responses to recent climate change. Annu. Rev. Ecol. Evol. Syst. 2006, 37, 637–669. [Google Scholar] [CrossRef] [Green Version]
  21. Chen, I.C.; Hill, J.K.; Ohlemüller, R.; Roy, D.B.; Thomas, C.D. Rapid range shifts of species associated with high levels of climate warming. Science 2011, 333, 1024–1026. [Google Scholar] [CrossRef]
  22. Seaton, P.T.; Hu, H.; Perner, H.; Pritchard, H.W. Ex Situ Conservation of Orchids in a Warming World. Bot. Rev. 2010, 76, 193–203. [Google Scholar] [CrossRef]
  23. Kull, T.; Selgis, U.; Peciña, M.V.; Metsare, M.; Ilves, A.; Tali, K.; Sepp, K.; Kull, K.; Shefferson, R.P. Factors influencing IUCN threat levels to orchids across Europe on the basis of national red lists. Ecol. Evol. 2016, 6, 6245–6265. [Google Scholar] [CrossRef] [PubMed]
  24. Fay, M.F.; Pailler, T.; Dixon, K.W. Orchid conservation: Making the links. Ann. Bot. 2015, 116, 377–379. [Google Scholar] [CrossRef] [Green Version]
  25. Djordjević, V.; Tsiftsis, S. The Role of Ecological Factors in Distribution and Abundance of Terrestrial Orchids. In Orchids Phytochemistry, Biology and Horticulture; Mérillon, J.-M., Kodja, H., Eds.; Springer Nature: Cham, Switzerland, 2020; pp. 1–71. ISBN 9783030112578. [Google Scholar]
  26. Nic Lughadha, E.; Bachman, S.P.; Leão, T.C.C.; Forest, F.; Halley, J.M.; Moat, J.; Acedo, C.; Bacon, K.L.; Brewer, R.F.A.; Gâteblé, G.; et al. Extinction risk and threats to plants and fungi. Plants People Planet 2020, 2, 389–408. [Google Scholar] [CrossRef]
  27. Bachman, S.P.; Nic Lughadha, E.M.; Rivers, M.C. Quantifying progress toward a conservation assessment for all plants. Conserv. Biol. 2018, 32, 516–524. [Google Scholar] [CrossRef] [Green Version]
  28. Cribb, P.; Lack, H.W.; Mabberley, D.J. The Flora Graeca Story. Sibthorp, Bauer, and Hawkins in the Levant. Kew Bull. 1999. [Google Scholar] [CrossRef]
  29. Strid, A. The botanical exploration of Greece. Plant Syst. Evol. 2020, 306, 1–23. [Google Scholar] [CrossRef]
  30. Kougioumoutzis, K.; Kokkoris, I.P.; Panitsa, M.; Kallimanis, A.; Strid, A.; Dimopoulos, P. Plant Endemism Centres and Biodiversity Hotspots in Greece. Biology (Basel) 2021, 10, 72. [Google Scholar] [CrossRef]
  31. Kougioumoutzis, K.; Kokkoris, I.P.; Panitsa, M.; Trigas, P.; Strid, A.; Dimopoulos, P. Plant Diversity Patterns and Conservation Implications under Climate-Change Scenarios in the Mediterranean: The Case of Crete (Aegean, Greece). Diversity 2020, 12, 270. [Google Scholar] [CrossRef]
  32. Kougioumoutzis, K.; Kokkoris, I.P.; Panitsa, M.; Trigas, P.; Strid, A.; Dimopoulos, P. Spatial phylogenetics, biogeographical patterns and conservation implications of the endemic flora of Crete (Aegean, Greece) under climate change scenarios. Biology (Basel) 2020, 9, 199. [Google Scholar] [CrossRef] [PubMed]
  33. Fassou, G.; Kougioumoutzis, K.; Iatrou, G.; Trigas, P.; Papasotiropoulos, V. Genetic diversity and range dynamics of Helleborus odorus subsp. cyclophyllus under different climate change scenarios. Forests 2020, 11, 620. [Google Scholar] [CrossRef]
  34. Stathi, E.; Kougioumoutzis, K.; Abraham, E.M.; Trigas, P.; Ganopoulos, I.; Avramidou, E.V.; Tani, E. Population genetic variability and distribution of the endangered Greek endemic Cicer graecum under climate change scenarios. AoB Plants 2020, 12, plaa007. [Google Scholar] [CrossRef]
  35. Tsiftsis, S.; Djordjević, V. Modelling sexually deceptive orchid species distributions under future climates: The importance of plant–pollinator interactions. Sci. Rep. 2020, 10, 10623. [Google Scholar] [CrossRef]
  36. Provan, J.; Bennett, K.D. Phylogeographic insights into cryptic glacial refugia. Trends Ecol. Evol. 2008, 23, 564–571. [Google Scholar] [CrossRef]
  37. Bai, Y.; Wei, X.; Li, X. Distributional dynamics of a vulnerable species in response to past and future climate change: A window for conservation prospects. PeerJ 2018, 2018, e4287. [Google Scholar] [CrossRef]
  38. Dimopoulos, P.; Raus, T.; Bergmeier, E.; Constantinidis, T.; Iatrou, G.; Kokkini, S.; Strid, A.; Tzanoudakis, D. Vascular Plants of Greece: An Annotated Checklist; Botanischer Garten und Botanischers Museum Berlin-Dahlem: Berlin, Germany; Hellenic Botanical Society: Athens, Greece, 2013. [Google Scholar]
  39. Dimopoulos, P.; Raus, T.; Bergmeier, E.; Constantinidis, T.; Iatrou, G.; Kokkini, S.; Strid, A.; Tzanoudakis, D. Vascular plants of Greece: An annotated checklist. Supplement. Willdenowia 2016, 46, 301–347. [Google Scholar] [CrossRef] [Green Version]
  40. Tsiftsis, S.; Antonopoulos, Z. Atlas of the Greek Orchids; Mediterraneo Editions Rethymnon: Crete, Greece, 2017. [Google Scholar]
  41. Tsiftsis, S.; Štípková, Z.; Kindlmann, P. Role of way of life, latitude, elevation and climate on the richness and distribution of orchid species. Biodivers. Conserv. 2019, 28, 75–96. [Google Scholar] [CrossRef]
  42. Breitkopf, H.; Onstein, R.E.; Cafasso, D.; Schlüter, P.M.; Cozzolino, S. Multiple shifts to different pollinators fuelled rapid diversification in sexually deceptive Ophrys orchids. New Phytol. 2015, 207, 377–389. [Google Scholar] [CrossRef]
  43. Paulus, H.; Gack, C. Pollinators as prepollinating isolation factors: Evolution and speciation in {IOphrys} (Orchidaceae). Isr. J. Bot. 1990, 39, 43–79. [Google Scholar]
  44. Paulus, H.F. Deceived males—Pollination biology of the Mediterranean orchid genus Ophrys (Orchidaceae). J. Eur. Orchid. 2006, 38, 303–351. [Google Scholar]
  45. Tsiftsis, S.; Tsiripidis, I. Temporal and spatial patterns of orchid species distribution in Greece: Implications for conservation. Biodivers. Conserv. 2020, 29, 3461–3489. [Google Scholar] [CrossRef]
  46. Renz, J. Zur Kenntnis der griechischen Orchideen. Repert. Specierum Nov. Regni Veg. 1928, 25, 225–270. [Google Scholar] [CrossRef]
  47. Danesch, O.; Danesch, E. Ophrys helenae: Die Geschichte einer Orchidee. Kosmos 1975, 75, 185–188. [Google Scholar]
  48. Danesch, O.; Danesch, E. Ophrys helenae (Orchidaceae), a Neglected Species of the Balkan Peninsula. Plant Syst. Evol. 1977, 127, 11–22. [Google Scholar] [CrossRef]
  49. Menzel, A.; Sparks, T.H.; Estrella, N.; Koch, E.; Aaasa, A.; Ahas, R.; Alm-Kübler, K.; Bissolli, P.; Braslavská, O.; Briede, A.; et al. European phenological response to climate change matches the warming pattern. Glob. Chang. Biol. 2006, 12, 1969–1976. [Google Scholar] [CrossRef]
  50. Urban, M.C. Accelerating extinction risk from climate change. Science 2015, 348, 571–573. [Google Scholar] [CrossRef] [Green Version]
  51. Thuiller, W.; Lavorel, S.; Araújo, M.B.; Sykes, M.T.; Prentice, I.C.; Thuiller, W.; Lavorel, S.; Arau, M.B. Climate change threats to plant diversity in Europe. Proc. Natl. Acad. Sci. USA 2005, 102, 8245–8250. [Google Scholar] [CrossRef] [Green Version]
  52. Bellard, C.; Bertelsmeier, C.; Leadley, P.; Thuiller, W.; Courchamp, F. Impacts of climate change on the future of biodiversity. Ecol. Lett. 2012, 15, 365–377. [Google Scholar] [CrossRef] [Green Version]
  53. Calinger, K.M.; Queenborough, S.; Curtis, P.S. Herbarium specimens reveal the footprint of climate change on flowering trends across north-central North America. Ecol. Lett. 2013, 16, 1037–1044. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  54. Enquist, B.J.; Feng, X.; Boyle, B.; Maitner, B.; Newman, E.A.; Jørgensen, P.M.; Roehrdanz, P.R.; Thiers, B.M.; Burger, J.R.; Corlett, R.T.; et al. The commonness of rarity: Global and future distribution of rarity across land plants. Sci. Adv. 2019, 5, eaaz0414. [Google Scholar] [CrossRef] [Green Version]
  55. Newbold, T. Future effects of climate and land-use change on terrestrial vertebrate community diversity under different scenarios. Proc. R. Soc. B Biol. Sci. 2018, 285, 20180792. [Google Scholar] [CrossRef] [PubMed]
  56. Newbold, T.; Hudson, L.N.; Contu, S.; Hill, S.L.L.; Beck, J.; Liu, Y.; Meyer, C.; Phillips, H.R.P.; Scharlemann, J.P.W.; Purvis, A. Widespread winners and narrow-ranged losers: Land use homogenizes biodiversity in local assemblages worldwide. PLoS Biol. 2018, 16, e2006841. [Google Scholar] [CrossRef] [PubMed]
  57. Newbold, T.; Oppenheimer, P.; Etard, A.; Williams, J.J. Tropical and Mediterranean biodiversity is disproportionately sensitive to land-use and climate change. Nat. Ecol. Evol. 2020. [Google Scholar] [CrossRef] [PubMed]
  58. Powers, R.P.; Jetz, W. Global habitat loss and extinction risk of terrestrial vertebrates under future land-use-change scenarios. Nat. Clim. Chang. 2019, 9, 323–329. [Google Scholar] [CrossRef]
  59. Pyšek, P.; Hulme, P.E.; Simberloff, D.; Bacher, S.; Blackburn, T.M.; Carlton, J.T.; Dawson, W.; Essl, F.; Foxcroft, L.C.; Genovesi, P.; et al. Scientists’ warning on invasive alien species. Biol. Rev. 2020, 95, 1511–1534. [Google Scholar] [CrossRef]
  60. Seebens, H.; Blackburn, T.M.; Dyer, E.E.; Genovesi, P.; Hulme, P.E.; Jeschke, J.M.; Pagad, S.; Pyšek, P.; Winter, M.; Arianoutsou, M.; et al. No saturation in the accumulation of alien species worldwide. Nat. Commun. 2017, 8, 14435. [Google Scholar] [CrossRef] [PubMed]
  61. Hutchings, M.J.; Robbirt, K.M.; Roberts, D.L.; Davy, A.J. Vulnerability of a specialized pollination mechanism to climate change revealed by a 356-year analysis. Bot. J. Linn. Soc. 2018, 186, 498–509. [Google Scholar] [CrossRef]
  62. Kolanowska, M. Niche Conservatism and the Future Potential Range of Epipactis helleborine (Orchidaceae). PLoS ONE 2013, 8, e77352. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  63. Kolanowska, M.; Kras, M.; Lipińska, M.; Mystkowska, K.; Szlachetko, D.L.; Naczk, A.M. Global warming not so harmful for all plants-response of holomycotrophic orchid species for the future climate change. Sci. Rep. 2017, 7, 12704. [Google Scholar] [CrossRef] [Green Version]
  64. Reina-Rodríguez, G.A.; Rubiano Mejía, J.E.; Castro Llanos, F.A.; Soriano, I. Orchids distribution and bioclimatic niches as a strategy to climate change in areas of tropical dry forest in Colombia. Lankesteriana 2017, 17, 17–47. [Google Scholar] [CrossRef] [Green Version]
  65. Vogt-Schilb, H.; Munoz, F.; Richard, F.; Schatz, B. Recent declines and range changes of orchids in Western Europe (France, Belgium and Luxembourg). Biol. Conserv. 2015, 190, 133–141. [Google Scholar] [CrossRef]
  66. Veresoglou, S.D.; Halley, J.M. Seed mass predicts migration lag of European trees. Ann. For. Sci. 2018, 75, 86. [Google Scholar] [CrossRef] [Green Version]
  67. Charitonidou, M.; Halley, J.M. What goes up must come down—Why high fecundity orchids challenge conservation beliefs. Biol. Conserv. 2020, 252, 108835. [Google Scholar] [CrossRef]
  68. You, J.; Qin, X.; Ranjitkar, S.; Lougheed, S.C.; Wang, M.; Zhou, W.; Ouyang, D.; Zhou, Y.; Xu, J.; Zhang, W.; et al. Response to climate change of montane herbaceous plants in the genus Rhodiola predicted by ecological niche modelling. Sci. Rep. 2018, 8, 5879. [Google Scholar] [CrossRef] [PubMed]
  69. Tzedakis, P.C. Vegetation variability in Greece during the last interglacial. Neth. J. Geosci. 2000, 79, 355–367. [Google Scholar] [CrossRef] [Green Version]
  70. Charitonidou, M.; Stara, K.; Kougioumoutzis, K.; Halley, J.M. Implications of salep collection for the conservation of the Elder-flowered orchid (Dactylorhiza sambucina) in Epirus, Greece. J. Biol. Res. Thessalon. 2019, 26, 1–13. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  71. Gerasimidis, A.; Panajiotidis, S.; Fotiadis, G.; Korakis, G. Review of the Quaternary vegetation history of Epirus (NW Greece). Phytol. Balc. 2009, 15, 29–37. [Google Scholar]
  72. van der Meer, S.; Jacquemyn, H.; Carey, P.D.; Jongejans, E. Recent range expansion of a terrestrial orchid corresponds with climate-driven variation in its population dynamics. Oecologia 2016, 181, 435–448. [Google Scholar] [CrossRef] [PubMed]
  73. Kolanowska, M.; Jakubska-Busse, A. Is the lady’s-slipper orchid (Cypripedium calceolus) likely to shortly become extinct in Europe?—Insights based on ecological niche modelling. PLoS ONE 2020, 15, e0228420. [Google Scholar] [CrossRef]
  74. Djordjević, V.; Tsiftsis, S.; Lakušić, D.; Jovanović, S.; Stevanović, V. Orchid species richness and composition in relation to vegetation types. Wulfenia 2020, 27, 183–210. [Google Scholar]
  75. Ongaro, S.; Martellos, S.; Bacaro, G.; De Agostini, A.; Cogoni, A.; Cortis, P. Distributional pattern of sardinian orchids under a climate change scenario. Community Ecol. 2018, 19, 223–232. [Google Scholar] [CrossRef]
  76. Tzortzaki, A.E.; Vokou, D.; Halley, J.M. Campanula lingulata populations on Mt. Olympus, Greece: Where’s the “abundant centre”? J. Biol. Res. 2017, 24, 1–13. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  77. Cowie, J. Climate Change Biological and Human Aspects; Cambridge University Press: Cambridge, UK, 2007; ISBN 9780521696197. [Google Scholar]
  78. Qu, Y.; Luo, X.; Zhang, R.; Song, G.; Zou, F.; Lei, F. Lineage diversification and historical demography of a montane bird Garrulax elliotii—Implications for the Pleistocene evolutionary history of the eastern Himalayas. BMC Evol. Biol. 2011, 11, 174. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  79. Elith, J.; Kearney, M.; Phillips, S. The art of modelling range-shifting species. Methods Ecol. Evol. 2010, 1, 330–342. [Google Scholar] [CrossRef]
  80. Aiello-Lammens, M.E.; Boria, R.A.; Radosavljevic, A.; Vilela, B.; Anderson, R.P. spThin: An R package for spatial thinning of species occurrence records for use in ecological niche models. Ecography 2015, 38, 541–545. [Google Scholar] [CrossRef]
  81. Robertson, M.P.; Visser, V.; Hui, C. Biogeo: An R package for assessing and improving data quality of occurrence record datasets. Ecography 2016, 39, 394–401. [Google Scholar] [CrossRef] [Green Version]
  82. Varela, S.; Anderson, R.P.; García-Valdés, R.; Fernández-González, F. Environmental filters reduce the effects of sampling bias and improve predictions of ecological niche models. Ecography 2014, 37, 1084–1091. [Google Scholar] [CrossRef]
  83. Zizka, A.; Antonelli, A.; Silvestro, D. sampbias, a method for quantifying geographic sampling biases in species distribution data. Ecography 2021, 44, 25–32. [Google Scholar] [CrossRef]
  84. Hijmans, R.J.; Cameron, S.E.; Parra, J.L.; Jones, P.G.; Jarvis, A. Very high resolution interpolated climate surfaces for global land areas. Int. J. Climatol. 2005, 25, 1965–1978. [Google Scholar] [CrossRef]
  85. Karger, D.N.; Conrad, O.; Böhner, J.; Kawohl, T.; Kreft, H.; Soria-Auza, R.W.; Zimmermann, N.E.; Linder, H.P.; Kessler, M. Climatologies at high resolution for the earth’s land surface areas. Sci. Data 2017, 4, 170122. [Google Scholar] [CrossRef] [Green Version]
  86. Morales-Barbero, J.; Vega-Álvarez, J. Input matters matter: Bioclimatic consistency to map more reliable species distribution models. Methods Ecol. Evol. 2019, 10, 212–224. [Google Scholar] [CrossRef]
  87. Araújo, M.B.; Anderson, R.P.; Barbosa, A.M.; Beale, C.M.; Dormann, C.F.; Early, R.; Garcia, R.A.; Guisan, A.; Maiorano, L.; Naimi, B.; et al. Standards for distribution models in biodiversity assessments. Sci. Adv. 2019, 5, eaat4858. [Google Scholar] [CrossRef] [Green Version]
  88. Title, P.O.; Bemmels, J.B. ENVIREM: An expanded set of bioclimatic and topographic variables increases flexibility and improves performance of ecological niche modeling. Ecography 2018, 41, 291–307. [Google Scholar] [CrossRef] [Green Version]
  89. McSweeney, C.F.; Jones, R.G.; Lee, R.W.; Rowell, D.P. Selecting CMIP5 GCMs for downscaling over multiple regions. Clim. Dyn. 2015, 44, 3237–3260. [Google Scholar] [CrossRef] [Green Version]
  90. Hengl, T.; Mendes de Jesus, J.; Heuvelink, G.B.M.; Ruiperez Gonzalez, M.; Kilibarda, M.; Blagotić, A.; Shangguan, W.; Wright, M.N.; Geng, X.; Bauer-Marschallinger, B.; et al. SoilGrids250m: Global Gridded Soil Information Based on Machine Learning. PLoS ONE 2017, 12, e0169748. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  91. Jarvis, A.; Reuter, H.I.; Nelson, A.; Guevara, E. Hole-Filled SRTM for the Globe Version 4: Data Grid; GIAR Consortium for Spatial Information: Washington, DC, USA, 2008; Available online: http://srtm.csi.cgiar.org/ (accessed on 10 February 2021).
  92. Hijmans, R.J. Package ‘raster’—Geographic Data Analysis and Modeling. CRAN Repos. 2019. Available online: https://rspatial.org/raster (accessed on 10 February 2021).
  93. Evans, J.S. ‘spatialEco’ R Package Version 1.2-0. 2019. Available online: https://github.com/jeffreyevans/spatialEco (accessed on 10 February 2021).
  94. Pebesma, E. Simple features for R: Standardized support for spatial vector data. R J. 2018, 10, 439. [Google Scholar] [CrossRef] [Green Version]
  95. Ross, N. ‘fasterize’: Fast Polygon to Raster Conversion. R Package version 1.0.3. 2020. Available online: https://github.com/ecohealthalliance/fasterize (accessed on 10 February 2021).
  96. Bornovas, I.; Rondogianni-Tsiambaou, T. Geological Map of Greece. Scale 1:500,000, 2nd ed.; Institute of Geology and Mineral Exploration, Division of General Geology and Economic Geology: Athens, Greece, 1983. [Google Scholar]
  97. Dormann, C.F.; Elith, J.; Bacher, S.; Buchmann, C.; Carl, G.; Carré, G.; Marquéz, J.R.G.; Gruber, B.; Lafourcade, B.; Leitão, P.J.; et al. Collinearity: A review of methods to deal with it and a simulation study evaluating their performance. Ecography 2013, 36, 27–46. [Google Scholar] [CrossRef]
  98. Naimi, B.; Hamm, N.A.S.; Groen, T.A.; Skidmore, A.K.; Toxopeus, A.G. Where is positional uncertainty a problem for species distribution modelling? Ecography 2014, 37, 191–203. [Google Scholar] [CrossRef]
  99. Thuiller, W.; Lafourcade, B.; Engler, R.; Araújo, M.B. BIOMOD—A platform for ensemble forecasting of species distributions. Ecography 2009, 32, 369–373. [Google Scholar] [CrossRef]
  100. Di Cola, V.; Broennimann, O.; Petitpierre, B.; Breiner, F.T.; D’Amen, M.; Randin, C.; Engler, R.; Pottier, J.; Pio, D.; Dubuis, A.; et al. ecospat: An R package to support spatial analyses and modeling of species niches and distributions. Ecography 2017, 40, 774–787. [Google Scholar] [CrossRef]
  101. Araújo, M.B.; New, M. Ensemble forecasting of species distributions. Trends Ecol. Evol. 2007, 22, 42–47. [Google Scholar] [CrossRef]
  102. Barbet-Massin, M.; Jiguet, F.; Albert, C.H.; Thuiller, W. Selecting pseudo-absences for species distribution models: How, where and how many? Methods Ecol. Evol. 2012, 3, 327–338. [Google Scholar] [CrossRef]
  103. Valavi, R.; Elith, J.; Lahoz-Monfort, J.J.; Guillera-Arroita, G. blockCV: An r package for generating spatially or environmentally separated folds for k-fold cross-validation of species distribution models. Methods Ecol. Evol. 2019, 10, 225–232. [Google Scholar] [CrossRef] [Green Version]
  104. Breiner, F.T.; Nobis, M.P.; Bergamini, A.; Guisan, A. Optimizing ensembles of small models for predicting the distribution of species with few occurrences. Methods Ecol. Evol. 2018, 9, 802–808. [Google Scholar] [CrossRef] [Green Version]
  105. Breiner, F.T.; Guisan, A.; Bergamini, A.; Nobis, M.P. Overcoming limitations of modelling rare species by using ensembles of small models. Methods Ecol. Evol. 2015, 6, 1210–1218. [Google Scholar] [CrossRef]
  106. Breiner, F.T.; Guisan, A.; Nobis, M.P.; Bergamini, A. Including environmental niche information to improve IUCN Red List assessments. Divers. Distrib. 2017, 23, 484–495. [Google Scholar] [CrossRef] [Green Version]
  107. Allouche, O.; Tsoar, A.; Kadmon, R. Assessing the accuracy of species distribution models: Prevalence, kappa and the true skill statistic (TSS). J. Appl. Ecol. 2006, 43, 1223–1232. [Google Scholar] [CrossRef]
  108. Raes, N.; Ter Steege, H. A null-model for significance testing of presence-only species distribution models. Ecography 2007, 30, 727–736. [Google Scholar] [CrossRef]
  109. QGIS Development Team. QGIS Geographic Information System. Open Source Geospatial Foundation. 2019. Available online: http://www.qgis.org (accessed on 10 February 2021).
Plants 10 00470 g001 550
Figure 1. Habitat suitability map for Ophrys helenae, for the combination of CHELSA climate database and geographical thinning procedure. Orange border delineates the species’ distributional area based on the Atlas of Greek Orchids [40].
Figure 1. Habitat suitability map for Ophrys helenae, for the combination of CHELSA climate database and geographical thinning procedure. Orange border delineates the species’ distributional area based on the Atlas of Greek Orchids [40].
Plants 10 00470 g001
Plants 10 00470 g002 550
Figure 2. Bioclimatic consistency maps for both thinning procedures (geographical and environmental) of Ophrys helenae. From top-left to bottom-right: geographically-thinned, based on WorldClim; geographically-thinned, based on CHELSA; Environmentally-thinned, based on WorldClim; environmentally-thinned, based on CHELSA.
Figure 2. Bioclimatic consistency maps for both thinning procedures (geographical and environmental) of Ophrys helenae. From top-left to bottom-right: geographically-thinned, based on WorldClim; geographically-thinned, based on CHELSA; Environmentally-thinned, based on WorldClim; environmentally-thinned, based on CHELSA.
Plants 10 00470 g002
Plants 10 00470 g003 550
Figure 3. Past, present and future potential distribution maps for Ophrys helenae, based on the geographically thinned occurrence data and climate data from CHELSA (CHGEO) that show the transition from one time-slice to another. In the central panel, green coloring indicates the cells the species currently potentially occupies, while the orange lines delineate the species’ distributional area. Left-hand panel: past potential distribution maps for showing the transition from the Last Interglacial (ca. 120–140 Kya) to the Last Glacial Maximum Maximum (ca. 20 Kya) and from the Last Glacial Maximum (ca. 20 Kya) to the present time-period. Red, blue and grey coloring depict cells that show loss, gain, or remained stable, respectively, compared to current the previous conditions. Future potential distribution maps represent combinations of three Global Circulation Models (BCC, CCSM4, and HadGEM2-ES) and two Representative Concentration Pathways (RCP 2.6 and RCP 8.5) showing the transition from the present time-period to each respective GCM and RCP combination. The red grid cells in the right-hand panel indicate that the species is currently present at these areas but will not be in the future. The grey grid cells in the right-hand panel indicate that the species is not currently present at these areas and it will not be in the future. The blue grid cells in the right-hand panel indicate that the species currently occupies and will continue to occupy these areas in the future.
Figure 3. Past, present and future potential distribution maps for Ophrys helenae, based on the geographically thinned occurrence data and climate data from CHELSA (CHGEO) that show the transition from one time-slice to another. In the central panel, green coloring indicates the cells the species currently potentially occupies, while the orange lines delineate the species’ distributional area. Left-hand panel: past potential distribution maps for showing the transition from the Last Interglacial (ca. 120–140 Kya) to the Last Glacial Maximum Maximum (ca. 20 Kya) and from the Last Glacial Maximum (ca. 20 Kya) to the present time-period. Red, blue and grey coloring depict cells that show loss, gain, or remained stable, respectively, compared to current the previous conditions. Future potential distribution maps represent combinations of three Global Circulation Models (BCC, CCSM4, and HadGEM2-ES) and two Representative Concentration Pathways (RCP 2.6 and RCP 8.5) showing the transition from the present time-period to each respective GCM and RCP combination. The red grid cells in the right-hand panel indicate that the species is currently present at these areas but will not be in the future. The grey grid cells in the right-hand panel indicate that the species is not currently present at these areas and it will not be in the future. The blue grid cells in the right-hand panel indicate that the species currently occupies and will continue to occupy these areas in the future.
Plants 10 00470 g003
Plants 10 00470 g004 550
Figure 4. Distribution centroids for current and future projected distributions of Ophrys helenae in Greece, based on the combination of CHELSA climatic database and geographical thinning procedure (CHGEO) for occurrence data.
Figure 4. Distribution centroids for current and future projected distributions of Ophrys helenae in Greece, based on the combination of CHELSA climatic database and geographical thinning procedure (CHGEO) for occurrence data.
Plants 10 00470 g004
Plants 10 00470 g005 550
Figure 5. Occurrences of Helen’s bee orchid in Greece, according to the database of the Orchid Flora of Greece project (Atlas of Greek Orchids [40]). Picture of Ophrys helenae individuals blooming in Ioannina, NW Greece, April 2019 (photo credits: Kalliopi Stara).
Figure 5. Occurrences of Helen’s bee orchid in Greece, according to the database of the Orchid Flora of Greece project (Atlas of Greek Orchids [40]). Picture of Ophrys helenae individuals blooming in Ioannina, NW Greece, April 2019 (photo credits: Kalliopi Stara).
Plants 10 00470 g005
Table
Table 1. The two most important variables for each combination of climatic database and thinning procedure. CHGEO: CHELSA—Geographical Thinning, CHENV: CHELSA—Environmental Thinning, WCGEO: WorldClim—Geographical Thinning, WCENV: WorldClim—Environmental Thinning, PWM: precipitation of wettest month; PETWQ: mean monthly potential evapotraspiration of wettest quarter; SOC-5: predicted mean value of soil organic carbon mass fraction at standard depth of 5 cm.
Table 1. The two most important variables for each combination of climatic database and thinning procedure. CHGEO: CHELSA—Geographical Thinning, CHENV: CHELSA—Environmental Thinning, WCGEO: WorldClim—Geographical Thinning, WCENV: WorldClim—Environmental Thinning, PWM: precipitation of wettest month; PETWQ: mean monthly potential evapotraspiration of wettest quarter; SOC-5: predicted mean value of soil organic carbon mass fraction at standard depth of 5 cm.
Database/ThinningVariable CodeVariable Importance
CHGEOPWM0.874
SOC-50.849
CHENVPWM0.899
PETWQ0.843
WCGEOPWM0.929
PETWQ0.855
WCENVPWM0.908
PETWQ0.852
Table
Table 2. Percentages of range loss and gain, a overall range change for past and future projections of Ophrys helenae distribution. For the past, Last Interglacial and Last Glacial Maximum (LIG and LGM, respectively) are presented, while for the future, all selected Global Circulation Models (GCMs) and Representative Concentration Pathways (RCPs) are presented (BC: BCC-CSM-1, CC: CCSM4, HE: HadGEM2-ES). The presented values are for the CHELSA climatic database, and for the geographical thinning procedure (CHGEO).
Table 2. Percentages of range loss and gain, a overall range change for past and future projections of Ophrys helenae distribution. For the past, Last Interglacial and Last Glacial Maximum (LIG and LGM, respectively) are presented, while for the future, all selected Global Circulation Models (GCMs) and Representative Concentration Pathways (RCPs) are presented (BC: BCC-CSM-1, CC: CCSM4, HE: HadGEM2-ES). The presented values are for the CHELSA climatic database, and for the geographical thinning procedure (CHGEO).
Database/ThinningTime SliceTransitionGCMRange Loss (%)Range Gain (%)Range Change (%)
CHGEOPastLIG to LGM 4.480.67−3.81
LGM to Present 29.882.14−27.74
FuturePresent to RCP 2.6BC10.4322.6512.22
CC4.6919.3814.68
HE35.2116.97−18.24
Present to RCP 8.5BC40.0419.96−20.08
CC47.2714.37−32.91
HE72.076.43−65.63
Table
Table 3. Mean altitude for past, present and future projections of Ophrys helenae distribution. For the past, Last Interglacial and Last Glacial Maximum (LIG and LGM, respectively) are presented, while for the future, all selected Global Circulation Models (GCMs) and Representative Concentration Pathways (RCPs). BC: BCC-CSM-1, CC: CCSM4, HE: HadGEM2-ES. The presented values are for the CHELSA climatic database and geographical thinning procedure (CHGEO).
Table 3. Mean altitude for past, present and future projections of Ophrys helenae distribution. For the past, Last Interglacial and Last Glacial Maximum (LIG and LGM, respectively) are presented, while for the future, all selected Global Circulation Models (GCMs) and Representative Concentration Pathways (RCPs). BC: BCC-CSM-1, CC: CCSM4, HE: HadGEM2-ES. The presented values are for the CHELSA climatic database and geographical thinning procedure (CHGEO).
Database/ThinningTime SliceRCPPeriod/GCMMean Altitude (m)
CHGEOPast LIG671
LGM679
Current Present637
FutureRCP 2.6BC653
CC634
HE866
RCP 8.5BC712
CC860
HE915
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Charitonidou, M.; Kougioumoutzis, K.; Halley, J.M. An Orchid in Retrograde: Climate-Driven Range Shift Patterns of Ophrys helenae in Greece. Plants 2021, 10, 470. https://doi.org/10.3390/plants10030470

AMA Style

Charitonidou M, Kougioumoutzis K, Halley JM. An Orchid in Retrograde: Climate-Driven Range Shift Patterns of Ophrys helenae in Greece. Plants. 2021; 10(3):470. https://doi.org/10.3390/plants10030470

Chicago/Turabian Style

Charitonidou, Martha, Konstantinos Kougioumoutzis, and John M. Halley. 2021. "An Orchid in Retrograde: Climate-Driven Range Shift Patterns of Ophrys helenae in Greece" Plants 10, no. 3: 470. https://doi.org/10.3390/plants10030470

APA Style

Charitonidou, M., Kougioumoutzis, K., & Halley, J. M. (2021). An Orchid in Retrograde: Climate-Driven Range Shift Patterns of Ophrys helenae in Greece. Plants, 10(3), 470. https://doi.org/10.3390/plants10030470

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop