- Research article
- Open Access
A phyloclimatic study of Cyclamen
© Yesson and Culham; licensee BioMed Central Ltd. 2006
- Received: 27 February 2006
- Accepted: 20 September 2006
- Published: 20 September 2006
The impact of global climate change on plant distribution, speciation and extinction is of current concern. Examining species climatic preferences via bioclimatic niche modelling is a key tool to study this impact. There is an established link between bioclimatic niche models and phylogenetic diversification. A next step is to examine future distribution predictions from a phylogenetic perspective. We present such a study using Cyclamen (Myrsinaceae), a group which demonstrates morphological and phenological adaptations to its seasonal Mediterranean-type climate. How will the predicted climate change affect future distribution of this popular genus of garden plants?
We demonstrate phylogenetic structure for some climatic characteristics, and show that most Cyclamen have distinct climatic niches, with the exception of several wide-ranging, geographically expansive, species. We reconstruct climate preferences for hypothetical ancestral Cyclamen. The ancestral Cyclamen lineage has a preference for the seasonal Mediterranean climate characteristic of dry summers and wet winters.
Future bioclimatic niches, based on BIOCLIM and Maxent models, are examined with reference to a future climate scenario for the 2050s. Over the next 50 years we predict a northward shift in the area of climatic suitability, with many areas of current distribution becoming climatically unsuitable. The area of climatic suitability for every Cyclamen species is predicted to decrease. For many species, there may be no areas with a suitable climate regardless of dispersal ability, these species are considered to be at high risk of extinction. This risk is examined from a phylogenetic perspective.
Examining bioclimatic niches from a phylogenetic perspective permits novel interpretations of these models. In particular, reconstruction of ancestral niches can provide testable hypothesis about the historical development of lineages. In the future we can expect a northwards shift in climatic suitability for the genus Cyclamen. If this proves to be the case then dispersal is the best chance of survival, which seems highly unlikely for ant-dispersed Cyclamen. Human-assisted establishment of Cyclamen species well outside their native ranges offers hope and could provide the only means of dispersal to potentially suitable future environments. Even without human intervention the phylogenetic perspective demonstrates that major lineages could survive climate change even if many species are lost.
- Daily Precipitation
- Extinction Risk
- Future Climate Scenario
- Area Loss
- Maxent Model
The prospect of global climate change has directed interest towards investigating the impact of the environment on floral and faunal distribution, speciation and extinction [1, 2]. One way to investigate species response to climate is through examination of climatic preferences by constructing bioclimatic niche models (these are also known as species distribution models or environmental niche models) [3–6]. These methods establish preferences of a given species, based on its known distribution, and provide a model of the climate parameters correlating with this. One of the earliest and simplest methods is BIOCLIM, which uses the minimum and maximum (or 95th percentiles) of observed values for each climate parameter to define the environmental niche [3, 7]. BIOCLIM's models are more conducive to interpretation than some more complicated methodologies , although many comparisons demonstrate that more complex algorithms such as Maxent can have greater predictive value under most conditions [4, 9]. Once built, the models can be used in conjunction with different climate scenarios and timeframes to estimate past [10–13], present [9, 14, 15] and future [2, 6, 16–18] distributions.
There is an established link between bioclimatic niche models and phylogenetic diversification. Peterson et al.  suggest that bioclimatic envelopes are heritable and are conserved across evolutionary time. Martinez-Meyer et al.  demonstrated this using bioclimatic niche models of Passerina birds to successfully predict the distribution of sister species. This is further supported by a wider link between climate and phylogenetic diversification . Many researchers are now examining species' climatic preferences across phylogenetic trees [10, 20–23]. Of these studies, those concerned with distributions have focussed on present or past distributions [10, 20–23]. Yet bioclimatic niche models have also been used to predict future distributions, and their impact on extinction risk [2, 6, 24]. A clear next step is to examine future distribution predictions from a phylogenetic perspective.
Distribution data were gathered for Cyclamen from herbarium specimens and distribution maps and these data were used to develop climate profiles for each species. The difference between direct observations and distribution maps was tested and shown to give similar results. BIOCLIM  bioclimatic niche models were produced using modelled present day climate data from the Climate Research Unit . These models and underlying climatic parameters were examined from a phylogenetic perspective on the recent, complete species-level phylogeny of Compton et al. , and a reconstruction of the ancestral Cyclamen bioclimatic niche model was performed.
Each species' bioclimatic niche model was examined within a future climate scenario from the Intergovernmental Panel on Climate Change for the 2050s . Future areas of climatic suitability were compared with present day suitable areas and actual distribution. Many species show no future area of climatic suitability and these are considered to be at high risk of extinction due to climate change. This extinction risk measure is examined from a phylogenetic perspective. For full details of methods see methods section below.
The models built via the distribution (area) maps of Grey-Wilson  provide a good alternative to those built from direct (point) observations, they also provide the opportunity to study models for all 21 species of Cyclamen in a consistent fashion. If there is any bias, it is likely to be due to the distribution maps encompassing an area greater than the true range, which will have the affect of widening the bioclimatic envelope . However, the high agreement of point and area data suggest this effect is minimal, and all subsequent analyses are performed using models built from distribution maps.
Cyclamen climate preferences.
Dev Tmp (0.36)*
Pcp Wet Mnth (0.68)
Tmp Warm Mnth (0.44)
Mean Tmp (0.36)*
Dev Pcp (0.60)
Num Mnth Tmp >0 (0.37)*
Tmp Cool Mnth (0.34)*
Pcp Warm Mnth (0.26)*
Tmp Range (0.37)*
Pcp Dry Mnth (0.35)*
Min Tmp Cool Mnth (0.34)*
Pcp Cool Mnth (0.61)
Mean Daily Pcp (0.52)
Hi Tmp Warm Mnth (0.41)
C. creticum (Dörfl.) Hildebr.
C. repandum Sm.
C. balearicum Willk.
C. rohlfsianum Aschers.
C. graecum Link
C. persicum Mill.
C. somalense Thulin & Warfa
C. hederifolium Ait.
C. africanum Boiss. & Reut.
C. purpurascens Mill.
C. colchicum Alboff
C. parviflorum Pobedimova
C. pseudibericum Hildebr.
C. cyprium Kotschy
C. libanoticum Hildebr.
C. cilicium Boiss. & Heldr.
C. mirabile Hildebr.
C. intaminatum (Meikle) C. Grey-Wilson
C. trochopteranthum Schwarz
C. coum Mill.
C. elegans Boiss. & Buhse
The substantial difference in winter minimum temperatures (-8 to +21°C) is indicative of the contrast between both inland vs. coastal distributions and latitudinal/altitudinal variation. C. somalense is substantially nearer the equator than any other species and experiences warm summers and winters, whilst C. purpurascens, the most northerly species, tolerates freezing winter temperatures for at least four months of the year. C. creticum, an island species with its climate controlled by proximity to the sea, experiences one of the lowest annual temperature ranges with a mean summer temperature of +25°C compared with +11°C in the coldest winter month. C. colchicum, a species growing in continental Asia, experiences one of the highest annual temperature ranges with a mean summer temperature of +18°C compared with -8°C in the coldest month.
Ancestral reconstruction of climate preference.
Pcp Wet Mnth
Tmp Warm Mnth
Num Mnth Tmp >0
Tmp Cool Mnth
Pcp Warm Mnth
Pcp Dry Mnth
Min Tmp Cool Mnth
Pcp Cool Mnth
Mean Daily Pcp
Hi Tmp Warm Mnth
We can examine which of the extant species of Cyclamen is closest in climatic preference to our model reconstruction. It is not possible to use the full interpredictivity measure to determine which of the present day Cyclamen is most like the ancestral lineage. However, in the area selected by the model you can find the species C. hederifolium, C. repandum, C. coum, C. graecum and C. intaminatum. No extant species has a distribution similar to this area (Kappa < 0.02 for all species). However, a crude measure of model similarity based on the sum of squares of the differences of means for each climate layer suggests C. balearicum is overall the most similar to our ancestral reconstruction.
Similarity of bioclimatic niche models within selected groups of Cyclamen.
Group (Node No)
Subgen. Corticata (Node 15)
Subgen. Cyclamen (-)
Subgen. Gyrophoebe (-)
Subgen. Persicum (Node 8)
Subgen. Psilanthum (Node 2)
Random groups of 8 species
100 × 8
Random groups of 5 species
100 × 5
Random groups of 3 species
100 × 3
Future area projections
Occupied and modelled range sizes and extinction risk estimates for Cyclamen. Range sizes based on Grey-Wilson's maps .
Ancestral model reconstruction
Many Cyclamen share a preference for seasonal Mediterranean climate characteristics. The ancestral Cyclamen was probably well suited to this environment, and may have had a similar climatic preference to the extant species C. balearicum. A question remains on the timing of the origin of Cyclamen and how this compares with estimates of the emergence of the Mediterranean climate zone.
Although, for some ancestral model reconstructions, no present-day area was within the suitable envelope, this does not invalidate the model, as past climates could exhibit climate types not seen in the present day . However, failure to predict any suitable areas of the relevant palaeoclimate could demonstrate the invalidity of the models. Yesson and Culham  have demonstrated the plausibility of this general approach to estimates for ancestral areas based on reconstruction of climate preferences for ancestral lineages in other taxa.
There is not a phylogenetic pattern of inherited range size for Cyclamen. Several sister pairs of species are seen to have drastically different ranges, all but two species pairs on the phylogeny differ in their total area of distribution by more than 90%. This fits the pattern of schizo-endemic distribution of C. balearicum and C. repandum found by Thompson et al. , whereby a single wide ranging species is assumed to be progenitor to a narrow ranging isolated sister species. This is supported by the lack of overlap in bioclimatic niche models. This suggests a combination of geographic, edaphic and climatic constraints are limiting the maximum area of distribution of some species more than others.
Whilst many of the constituent climatic variables demonstrate phylogenetic conservancy, the resulting models which combine the individual variables into bioclimatic envelopes show little overlap between adjacent nodes. It may be the case that measures for comparing these envelopes are inadequate when the majority of species are climatically isolated from one-another. The only point of comparison of bioclimatic niche models that is currently used is to compare the areas that they select when overlaid into present day data [19, 39]. This is also true for the comparison measures used in model validation such as Cohen's Kappa . This is not a direct measure of model similarity. Furthermore such measures may 'hide' real differences for climate types not evident in the present day.
Some studies have found patterns of risk associated with a particular type of climate [18, 41]. Thuiller et al.  found species from Mediterranean climates to be at lower risk of extinction than species from other climates, whereas Thuiller et al.  found species in Mediterranean regions to be the most sensitive to the changing climate, with species in mountainous regions under great threat but other species having the potential to expand. For Cyclamen there is no particular climate-type which is more at risk than another. There is little overlap between any of the bioclimatic niche models of the high-risk Cyclamen species, many of these species occur in differing climates and geographies. For example the highest risk group (BIOCLIM models) includes the Mediterranean-type climate of C. rohlfsianum; the continental climate of the mountain dwelling C. colchicum; the island endemic C. cyprium; and the African C. somalense which is barely 10° north of the equator. However, for the Maxent models the wide-ranging, more generalist species are consistently assigned a high risk.
Measures of risk
These broad classifications of risk can hide some noteworthy detail. For example, although C. balearicum is designated Endangered (>50% area loss), it is noted that the BIOCLIM future area selection gives no suitable area for this species within the Balearics, and less than 5% of its present day French distribution is considered climatically suitable, indeed most of the predicted suitable climate occurs in Northern France. It is interesting to note that the Balearics are still suitable for other Cyclamen species. It is reassuring to note that the French populations of C. balearicum have a far higher genetic diversity than the Balearic populations, but the latter demonstrate the higher ecological diversity .
There is a link between current range size and extinction risk due to climate change. For the BIOCLIM models every species with a current distribution of under 35,000 km2 is considered the highest risk with the exception of C. pseudibericum. In contrast no species with a range size above 35,000 km2 loses all area. Range size and proportion of area lost is strongly negatively correlated for the BIOCLIM (Pearson correlation = -0.816, P < 0.001). In contrast the Maxent models show weak positive correlation of proportion loss and range size (Pearson correlation = 0.332, P = 0.142). The changes in range size show a distinct scattering across the cladograms whether modelled with BIOCLIM or Maxent (fig 12).
There are of course other factors than climate limiting species distributions. At the fine scale factors such as soil, aspect, and surrounding vegetation will have a role in determining the establishment of individuals given an otherwise suitable climate. Such fine-scale limiting factors are of secondary importance when modelling species distribution at a global or regional scale . These criticisms of species distribution modelling and the extinction risk measures based on them are discussed at length in other places [6, 19, 41–47] and we refer the reader to these references.
The BIOCLIM models are restrictive in their definitions of the edge of a suitable climate, predictions for all areas are essentially a binary presence/absence and do not allow any climatic parameter to be outside the observed range . For example, the BIOCLIM models predict no Cyclamen in the UK, but we know that three species are naturalised there. Alternative modelling algorithms such as Maxent have the advantage of producing probability measures rather than binary presence/absence, which provide more realistic edges to ranges . However, for the UK naturalised Cyclamen, the Maxent models also fail to select any of their naturalised distribution within the selected threshold.
Both modelling algorithms predict large range loss for all Cyclamen species, and the predicted future diversity maps are broadly similar (fig 11). However, predictions for individual species differ dramatically. The Maxent models show the highest area loss, yet the BIOCLIM models suggest many more species at the highest extinction risk. This kind of algorithm dependent difference has been observed in other studies [44, 49].
For the BIOCLIM predictions, it is evident that range size is a key factor in determining extinction risk due to climate change. Those species at high risk tend to have smaller distributions and smaller climatically suitable areas. This is in general agreement with others who have found a relationship between niche breadth and projected range loss [18, 41], (it is trivial to note that for Cyclamen the widest-ranging species have the broadest niches). However, the relationship between range size and range loss reported by Thuiller et al.  is a bell shaped curve with the largest ranges suffering substantial losses, although it is possible that no Cyclamen would fall into the category of truly wide-ranging. In contrast the Maxent models demonstrate the opposite pattern (albeit weakly), with widest-ranging species losing proportionally more area. It should be noted, however, that bioclimatic niche models may be less accurate for wide-ranging generalist species [4, 50].
Thuiller et al.  considered selecting the best algorithm for each species in turn based on validation scores. Such a policy, in this case, would combine the highest area loss of both methods producing high overall area loss and high numbers of species with 100% projected area loss. However it is to be noted that present day validation may not reflect model validity for the future [46, 51].
The main advantage of BIOCLIM in the context of Phyloclimatic modelling is its independent treatment of climatic variables. As far as we are aware, this is the only method that permits independent optimisation of environmental characteristics on a phylogeny following established phylogenetic reconstruction methods. These optimised variables can then be directly converted into a BIOCLIM bioclimatic niche models. Reconstruction of ancestral states on a phylogeny requires a distance metric to directly compare the output models of different species. There are many problems associated with using the more complex (and often more precise) models currently being developed. For example, many discard non-informative input parameters, meaning models are not directly comparable amongst species as different species models will be built with different subsets of variables. Other models are based on the amalgamation of multiple model outputs, again making direct comparison of models difficult.
Examining climate preference envelopes from a phylogenetic perspective brings new insight into this field. Ancestral climate envelope reconstruction can provide testable hypotheses about the development of lineages. In the case of Cyclamen, the ancestral lineage most probably developed in a Mediterranean type climate similar to parts of present day Greece and Turkey.
Climate parameters often show phylogenetic conservancy within Cyclamen. Despite this, the resulting bioclimatic envelopes show startling contrast between many sister species. The consequence is that while individual species are often at high risk, extinction is not predicted for most major lineages. The pattern of contrasting niches may not always be the case and it is likely that in groups demonstrating clearer patterns of niche conservatism it may be possible to assign risk of area loss due to a changing climate by phylogenetic proximity. This approach may be useful where data for full scale niche modelling is not available.
Many Cyclamen species are at high risk from the changing climate, and none will remain unaffected. Over the next fifty years we can expect a northwards shift in climatic suitability for the genus Cyclamen. Many species face the prospect of their local climate changing so much that their current distribution will be outside their current observed climate tolerance. If this proves to be the case then survival requires either adapting over a very short timeframe (implying exaptive changes or potentially suicidal rates of natural selection) or dispersing to new areas. The long distance migration option seems highly unlikely for ant-dispersed Cyclamen. This is worsened by the presence of many geographic barriers, the largest being the Mediterranean sea. Although all Cyclamen are listed in CITES, which protects them from trade exploitation from the wild, this does not per-se prevent their extinction though habitat loss by climate change. There are, however, still reasons for hope, the garden use of several Cyclamen species in more northern areas of Europe has facilitated the establishment of some Cyclamen species well outside the climatic zone of their native ranges. This suggests at least some climatic exaptation and the presence of a nucleus of plants in an area that will become climatically optimal for Cyclamen in the future. This provides a novel means of dispersal to, and colonization of, potentially suitable future environments.
The 21 species recognised by Compton et al.  were used for this study. Their phylogenetic study of Cyclamen is used to define relationships within Cyclamen for this paper. Their study produced phylogenies based upon nuclear (ITS) and plastid (trn L-F) DNA, as well as morphological data. The phylogenetic trees have strong topological similarity, differing mainly in degree of resolution. One of the two maximum parsimony trees based on the combined ITS + trn L-F dataset is used throughout this analysis. The alternative topology differs only in the position of C. graecum relative to C. rohlfsianum within series Persicum, our chosen topology agrees with findings of Clennett  with regards to this clade.
Locality data were collected for each of the 21 Cyclamen species. These data came from four sources. Firstly from Cyclamen Society collecting trips, which has made many trips over the period 1987–2004 covering Greece and the Greek islands, Israel, Sardinia and Turkey, providing 2,315 observations for 12 species . Secondly, data from the Global Biodiversity Information Facility , provides 232 observations for 6 species from 6 institutions covering Austria, Croatia, Germany, Greece and the Greek islands, Israel, and Italy (a full list of institutions in provided in the appendix [see Additional file 1]). Thirdly, data from specimen labels at the University of Reading Herbarium (RNG) containing 61 observations for 12 species covering the Balearics, Algeria, Austria, the Caucasus, Corsica, Croatia, Cyprus, Greece and the Greek islands, Israel, Italy and Turkey.
Points were excluded if they fell in an ocean or sea. Exclusions were minimal, comprising fewer than 3% of the total observations, primarily being records from near coastlines with poor geographic resolution (fig 3 shows an example of poor geographic resolution). Species with fewer than 10 observations were excluded from analyses, following Graham et al. . Data from the UK National Biodiversity Network for C. hederifolium, C. repandum and C. coum were excluded from the analysis as these represent records outside the native range of these species . Exclusions reduced the number of analysable Cyclamen down to 13 species.
The fourth source for distribution data was the detailed distribution maps of Grey-Wilson . This alternative approach allowed analysis of all Cyclamen species. These maps were digitally captured and pseudo distribution points were produced by taking the centroid of every overlapping quarter degree square. Quarter degrees were chosen as this is the resolution of the climate data used in this study. This produced distribution data for all 21 Cyclamen species. Due to its restricted distribution C. libanoticum has only 5 pseudo-localities following this methodology, and is the only species that falls below the 10 observations recommended by Graham et al. , however reasonable models were still produced for this species, and this species is included in all analyses.
Cyclamen range size was estimated by summing the area of each quarter degree grid square that overlapped with the digitised distribution maps.
Present day observed climate data were taken from an observed climatology dataset (known as CRU CL1) from the Climate Research Unit [32, 54]. Future data were taken from the Hadley centre general circulation model (HadCM3) from the IPCC website . The period 2040–2069 hereafter referred to as 2050s was chosen as the future time-frame. This allows sufficient time for the effects of climate change to be detectable, whilst minimising the overall timeframe to reduce uncertainties [55, 56] There are many different future climate scenarios based on different predictions of socio-economic factors, we have chosen the A2 scenario representing the mid-range of climate change severity . In order to produce future climate datasets which are directly comparable with the observed data, we applied modelled changes in climate to the present day observed data. The modelled changes were calculated by subtracting the 2050s prediction from directly comparable modelled present day climate following Peterson et al.  and IPCC recommendations .
Both the present and future datasets were processed into biologically meaningful climate profile parameters . Standard deviation of mean temperature, mean daily precipitation in wettest month, mean temperature in warmest month, mean temperature, standard deviation of mean precipitation, number of months with minimum temperature above freezing, mean temperature in coolest month, mean daily precipitation in warmest month, annual temperature range, mean daily precipitation in driest month, lowest temperature in coolest month, mean daily precipitation in coolest month, mean daily precipitation, highest temperature in warmest month. These are similar to Busby  and were processed using the "Climate Data Processor" plug-in developed by Tim Sutton as part of the Quantum GIS project .
The mean, standard deviation, minimum and maximum values were collected for all climatic parameters, for all species. Each value was independently optimised as a continuous character across the chosen phylogeny, using square change parsimony optimisation implemented in Mesquite , following the methodology of Graham et al. . A table showing each of these optimised values for every node on the tree is provided in the appendix [see Additional file 2]. These values define a BIOCLIM bioclimatic niche model  for each node on the tree. Models were created and manipulated using the BIOCLIM algorithm of the openModeller software package version 0.3.4 .
Additional models were constructed using the maximum entropy method (Maxent) for modelling species geographic distributions . This new method has been shown to be very effective for species distribution modelling . Models were built using the same data as the BIOCLIM models using Maxent version 2.3.0 . Thresholds were calculated for models to maximise the Kappa statistic. These thresholds were used to determine presence/absence for present day and future projections.
The quantitative convergence index (QVI)  was calculated for each optimised climate parameter. This measure is analogous to 1 minus the retention index for discrete characters. A randomisation test of this value was performed . A random shuffle of the terminal node values was optimised on the chosen topology, this was repeated 100 times for each character. The observed QVI was compared with the distribution of the randomisation replicates to test if the observed value is different from a random placement of data on the fixed phylogeny. If the observed QVI was outside the 95% confidence interval of the randomisation replicates then it was considered significantly different from random .
Only the variables demonstrating phylogenetic conservancy, with reference to the QVI, were used in the construction of ancestral bioclimatic niche models. This reduced the number of variables used in this stage of the modelling from 14 to 8.
Similarity of bioclimatic niche models was assessed using interpredictivity calculations [19, 39]. This involves the overlaying of locality records of one species (Species A) into the predicted area of another (Species B). Similarity is measured as the percentage of points of A falling within the prediction area of the model of B, and vice versa . Note this is not a symmetric measure, for example, if the envelope for A completely encompasses B then we could have B⊂A but A⊄B, therefore both directions of similarity were performed. Similarity was also assessed using Cohen's Kappa statistic [40, 64]. Mean similarity for a group of species was assessed by averaging similarity values of each pair of species within the group.
Summary 'hotspot' maps were produced by creating binary presence/absence maps for each species based on their bioclimatic niche model output, then summing the number of species predicted in each quarter degree cell. This was accomplished using software accessible through the openModeller project .
Present day bioclimatic niche models were examined within the 2050s climate scenario to produce a future distribution estimate based on climate suitability. Species with zero area predicted as climatically suitable are selected as high extinction risk . When predicting future distributions and extinction risk it is necessary to take dispersal rates into account [2, 66]. Plant migration rates have been estimated at 20–40 km over a 100 year timescale , more recently even lower rates have been suggested , it is highly unlikely that the ant-dispersed Cyclamen will exceed these dispersal rates. Our analysis uses a fifty year timescale and a geographic resolution of approx 25 × 25 km squares, so we employ a model permitting dispersal only within cells immediately adjacent to the current distribution. Therefore species which have future predictions which do not overlap with this the slightly expanded present distribution are also considered at risk, though this is regarded as lower risk than the case where there is no suitable area. Risk categories are assigned according to IUCN classifications based on area loss over 505 years [16, 69], namely: Extinction 100% loss, Critical >80%, Endangered >50%, Vulnerable >30%
We wish to thank: BBSRC for grant 45/BEP17792 for Biodiversity World; The University of Reading for CY's studentship; The Cyclamen Society for financial support and access to extensive unpublished data; The Biodiversity World team; Chloe Lee for populating a database of Cyclamen Society collection details; Sue Rumsey for digitising Reading herbarium specimen labels; Ronald Rutherford for validating the identity of the specimens; and two anonymous reviewers for their helpful comments; Peter Brewer for assistance with Maxent models.
- Davies TJ, Savolainen V, Chase MW, Moat J, Barraclough TG: Environmental energy and evolutionary rates in flowering plants. Proc R Soc Lond Ser B-Biol Sci. 2004, 271: 2195-2200.View ArticleGoogle Scholar
- Thomas CD, Cameron A, Green RE, Bakkenes M, Beaumont LJ, Collingham YC, Erasmus BFN, de Siqueira MF, Grainger A, Hannah L, Hughes L, Huntley B, van Jaarsveld AS, Midgley GF, Miles L, Ortega-Huerta MA, Peterson AT, Phillips OL, Williams SE: Extinction risk from climate change. Nature. 2004, 427: 145-148. 10.1038/nature02121.View ArticlePubMedGoogle Scholar
- Nix HA: A Biogeographic Analysis of Australian Elapid Snakes. Australian Flora and Fauna Series Number 7: Atlas of Elapid snakes of Australia. Edited by: Longmore, R. 1986, , Australian Government Publishing Service, Canberra., 4-15.Google Scholar
- Elith J, Graham CH, Anderson RP, Dudik M, Ferrier S, Guisan A, Hijmans RJ, Huettmann F, Leathwick JR, Lehmann A, Li J, Lohmann LG, Loiselle BA, Manion G, Moritz C, Nakamura M, Nakazawa Y, Overton JM, Peterson AT, Phillips SJ, Richardson K, Scachetti-Pereira R, Schapire RE, Soberon J, Williams S, Wisz MS, Zimmermann NE: Novel methods improve prediction of species' distributions from occurrence data. Ecography. 2006, 29: 129-151. 10.1111/j.2006.0906-7590.04596.x.View ArticleGoogle Scholar
- Guisan A, Thuiller W: Predicting species distribution: offering more than simple habitat models. Ecol Lett. 2005, 8: 993-1009. 10.1111/j.1461-0248.2005.00792.x.View ArticleGoogle Scholar
- Peterson AT, Tian H, Martínez-Meyer E, Soberon J, Sanchez-Cordero V, Huntley B: Modeling Distributional Shifts of Individual Species and Biomes. Climate change and biodiversity. Edited by: Lovejoy TE and Hannah LJ. 2005, New Haven Conn., Yale University Press, 211-229.Google Scholar
- Busby JR: BIOCLIM - A Bioclimatic Analysis and Prediction System. Nature Conservation: cost effective biological surveys and data analysis. Edited by: Margules CR and Austin, M.P. 1991, Canberra, CSIROGoogle Scholar
- Stockman AK, Beamer DA, Bond JE: An evaluation of a GARP model as an approach to predicting the spatial distribution of non-vagile invertebrate species. Divers Distrib. 2006, 12: 81-89. 10.1111/j.1366-9516.2006.00225.x.View ArticleGoogle Scholar
- Phillips SJ, Anderson RP, Schapire RE: Maximum entropy modeling of species geographic distributions. Ecol Model. 2006, 190: 231-259. 10.1016/j.ecolmodel.2005.03.026.View ArticleGoogle Scholar
- Hugall A, Moritz C, Moussalli A, Stanisic J: Reconciling paleodistribution models and comparative phylogeography in the Wet Tropics rainforest land snail Gnarosophia bellendenkerensis (Brazier 1875). Proc Natl Acad Sci U S A. 2002, 99: 6112-6117. 10.1073/pnas.092538699.PubMed CentralView ArticlePubMedGoogle Scholar
- Martinez-Meyer E, Townsend Peterson A, Hargrove WW: Ecological niches as stable distributional constraints on mammal species, with implications for Pleistocene extinctions and climate change projections for biodiversity. Glob Ecol Biogeogr. 2004, 13: 305-314. 10.1111/j.1466-822X.2004.00107.x.View ArticleGoogle Scholar
- Peterson AT, Soberon J, Sanchez-Cordero V: Conservatism of ecological niches in evolutionary time. Science. 1999, 285: 1265-1267. 10.1126/science.285.5431.1265.View ArticlePubMedGoogle Scholar
- Bonaccorso E, Koch I, Peterson AT: Pleistocene fragmentation of Amazon species’ ranges. Divers Distrib. 2006, 12: 157-164. 10.1111/j.1366-9516.2005.00212.x.View ArticleGoogle Scholar
- Robertson MP, Caithness N, Villet MH: A PCA-based modelling technique for predicting environmental suitability for organisms from presence records. Divers Distrib. 2001, 7: 1-2. 10.1046/j.1472-4642.2001.00094.x.View ArticleGoogle Scholar
- Stockwell DRB: Improving ecological niche models by data mining large environmental datasets for surrogate models. Ecol Model. 2006, 192: 188-196. 10.1016/j.ecolmodel.2005.05.029.View ArticleGoogle Scholar
- Thuiller W, Broennimann O, Hughes G, Alkemade JRM, Midgley GF, Corsi F: Vulnerability of African mammals to anthropogenic climate change under conservative land transformation assumptions. Glob Change Biol. 2006, 12: 424-440. 10.1111/j.1365-2486.2006.01115.x.View ArticleGoogle Scholar
- Broennimann O, Thuiller W, Hughes G, Midgley GF, Alkemade JMR, Guisan A: Do geographic distribution, niche property and life form explain plants' vulnerability to global change?. Global Change Biology. 2006, 12: 1079-1093. 10.1111/j.1365-2486.2006.01157.x.View ArticleGoogle Scholar
- Thuiller W, Lavorel S, Araújo MB, Sykes MT, Prentice IC: Climate change threats to plant diversity in Europe. Proc Natl Acad Sci U S A. 2005, 102: 8245-8250. 10.1073/pnas.0409902102.PubMed CentralView ArticlePubMedGoogle Scholar
- Martinez-Meyer E, Peterson AT, Navarro-Siguenza AG: Evolution of seasonal ecological niches in the Passerina buntings (Aves : Cardinalidae). Proc R Soc Lond Ser B-Biol Sci. 2004, 271: 1151-1157. 10.1098/rspb.2003.2564.View ArticleGoogle Scholar
- Graham CH, Ron SR, Santos JC, Schneider CJ, Moritz C: Integrating phylogenetics and environmental niche models to explore speciation mechanisms in dendrobatid frogs. Evolution. 2004, 58: 1781-1793. 10.1554/03-274.View ArticlePubMedGoogle Scholar
- Hardy CR, Linder HP: Intraspecific Variability and Timing in Ancestral Ecology Reconstruction: A Test Case from the Cape Flora. Syst Biol. 2005, 54: 299–316-View ArticlePubMedGoogle Scholar
- Hoffmann MH: Evolution of the realized climatic niche in the genus Arabidopsis (Brassicaceae). Evolution. 2005, 59: 1425-1436. 10.1554/05-033.PubMedGoogle Scholar
- Yesson C, Culham A: Phyloclimatic Modelling: Combining Phylogenetics and Bioclimatic Modelling. Systematic Biology. 2006, 55 (5): 1-18.View ArticleGoogle Scholar
- Peterson AT, Ortega-Huerta MA, Bartley J, Sanchez-Cordero V, Soberon J, Buddemeier RH, Stockwell DRB: Future projections for Mexican faunas under global climate change scenarios. Nature. 2002, 416: 626-629. 10.1038/416626a.View ArticlePubMedGoogle Scholar
- Compton JA, Clennett JCB, Culham A: Nomenclature in the dock. Overclassification leads to instability: a case study in the horticulturally important genus Cyclamen (Myrsinaceae). Bot J Lin Soc. 2004, 146: 339-349. 10.1111/j.1095-8339.2004.00322.x.View ArticleGoogle Scholar
- Grey-Wilson C: Cyclamen: a guide for gardeners, horticulturalists and botanists. 2003, London, Batsford, 224-NewGoogle Scholar
- DeBussche M, Garnier E, Thompson JD: Exploring the causes of variation in phenology and morphology in Mediterranean geophytes: a genus-wide study of Cyclamen. Bot J Linnean Soc. 2004, 145: 469-484. 10.1111/j.1095-8339.2004.00298.x.View ArticleGoogle Scholar
- Website: The Cyclamen Society. [http://www.cyclamen.org/]
- Thulin M, Warfa AM: Cyclamen (Primulaceae) in Tropical Africa. Plant Syst Evol. 1989, 166: 249-252. 10.1007/BF00935953.View ArticleGoogle Scholar
- Stace CA: New flora of the British Isles. 1997, Cambridge, Cambridge University Press, xxvii, 1130-2ndGoogle Scholar
- Affre L, Thompson JD, Debussche M: Genetic structure of continental and island populations of the Mediterranean endemic Cyclamen balearicum (Primulaceae). Am J Bot. 1997, 84: 437-451. 10.2307/2446019.View ArticlePubMedGoogle Scholar
- Website: The Climate Research Unit (Data archive). [http://www.cru.uea.ac.uk/cru/data/]
- Website: The Intergovernmental Panel on Climate Change. [http://www.ipcc.ch/]
- Peterson AT: Predicting the geography of species' invasions via ecological niche modeling. Q Rev Biol. 2003, 78: 419-433. 10.1086/378926.View ArticlePubMedGoogle Scholar
- Peterson AT, Martínez-Meyer E, Gonzalez-Salazar C: Reconstructing the Pleistocene geography of the Aphelocoma jays (Corvidae). Divers Distrib. 2004, 10: 237-246. 10.1111/j.1366-9516.2004.00097.x.View ArticleGoogle Scholar
- Thompson JD, Lavergne S, Affre L, Gaudeul M, Debussche M: Ecological differentiation of Mediterranean endemic plants. Taxon. 2005, 54: 967-976.View ArticleGoogle Scholar
- Faith DP: Conservation Evaluation and Phylogenetic Diversity. Biol Conserv. 1992, 61: 1-10. 10.1016/0006-3207(92)91201-3.View ArticleGoogle Scholar
- Pillon Y, Fay MF, Shipunov AB, Chase MW: Species diversity versus phylogenetic diversity: A practical study in the taxonomically difficult genus Dactylorhiza (Orchidaceae). Biol Conserv. 2006, 129: 4-13. 10.1016/j.biocon.2005.06.036.View ArticleGoogle Scholar
- Peterson AT, Vieglais DA: Predicting species invasions using ecological niche modeling: New approaches from bioinformatics attack a pressing problem. Bioscience. 2001, 51: 363-371. 10.1641/0006-3568(2001)051[0363:PSIUEN]2.0.CO;2.View ArticleGoogle Scholar
- Fielding AH, Bell JF: A review of methods for the assessment of prediction errors in conservation presence/absence models. Environ Conserv. 1997, 24: 38-49. 10.1017/S0376892997000088.View ArticleGoogle Scholar
- Thuiller W, Lavorel S, Araújo MB: Niche properties and geographical extent as predictors of species sensitivity to climate change. Glob Ecol Biogeogr. 2005, 14: 347-357. 10.1111/j.1466-822X.2005.00162.x.View ArticleGoogle Scholar
- Pearson RG, Dawson TP: Predicting the impacts of climate change on the distribution of species: are bioclimate envelope models useful?. Glob Ecol Biogeogr. 2003, 12: 361-371. 10.1046/j.1466-822X.2003.00042.x.View ArticleGoogle Scholar
- Pearson RG, Dawson TP: Bioclimate envelope models: what they detect and what they hide - response to Hampe (2004). Glob Ecol Biogeogr. 2004, 13: 471-473. 10.1111/j.1466-822X.2004.00112.x.View ArticleGoogle Scholar
- Pearson RG, Thuiller W, Araújo MB, Martínez-Meyer E, Brotons L, McClean C, Miles L, Segurado P, Dawson TP, Lees DC: Model-based uncertainty in species range prediction. J Biogeogr. 2006, 33: 1704-1711.View ArticleGoogle Scholar
- Hampe A: Bioclimate envelope models: what they detect and what they hide. Glob Ecol Biogeogr. 2004, 13: 469-471. 10.1111/j.1466-822X.2004.00090.x.View ArticleGoogle Scholar
- Araujo MB, Pearson RG, Thuiller W, Erhard M: Validation of species-climate impact models under climate change. Glob Change Biol. 2005, 11: 1504-1513. 10.1111/j.1365-2486.2005.01000.x.View ArticleGoogle Scholar
- Beaumont LJ, Hughes L, Poulsen M: Predicting species distributions: use of climatic parameters in BIOCLIM and its impact on predictions of species' current and future distributions. Ecol Model. 2005, 186: 250-269. 10.1016/j.ecolmodel.2005.01.030.View ArticleGoogle Scholar
- Guisan A, Zimmermann NE: Predictive habitat distribution models in ecology. Ecol Model. 2000, 135: 147-186. 10.1016/S0304-3800(00)00354-9.View ArticleGoogle Scholar
- Thuiller W: BIOMOD - optimizing predictions of species distributions and projecting potential future shifts under global change. Glob Change Biol. 2003, 9: 1353-1362. 10.1046/j.1365-2486.2003.00666.x.View ArticleGoogle Scholar
- Bomhard B, Richardson DM, Donaldson JS, Hughes GO, Midgley GF, Raimondo DC, Rebelo AG, Rouget M, Thuiller W: Potential impacts of future land use and climate change on the Red List status of the Proteaceae in the Cape Floristic Region, South Africa. Glob Change Biol. 2005, 11: 1452-1468. 10.1111/j.1365-2486.2005.00997.x.View ArticleGoogle Scholar
- Thuiller W: Patterns and uncertainties of species' range shifts under climate change. Glob Change Biol. 2004, 10: 2020-2027. 10.1111/j.1365-2486.2004.00859.x.View ArticleGoogle Scholar
- Clennett JCB: An analysis and revision of Cyclamen L. with emphasis on subgenus Gyrophoebe O. Schwarz. Bot J Linnean Soc. 2002, 138: 473-481. 10.1046/j.1095-8339.2002.00034.x.View ArticleGoogle Scholar
- Website: The Global Biodiversity Information Facility. [http://www.gbif.org/]
- New M, Hulme M, Jones P: Representing twentieth-century space-time climate variability. Part I: Development of a 1961-90 mean monthly terrestrial climatology. J Clim. 1999, 12: 829-856. 10.1175/1520-0442(1999)012<0829:RTCSTC>2.0.CO;2.View ArticleGoogle Scholar
- Houghton JT, Intergovernmental Panel on Climate Change. Working Group I.: Climate change 2001 : the scientific basis. 2001, Cambridge, U.K. ; New York, Published for the Intergovernmental Panel of Climate Change by Cambridge University Press, x, 881-Google Scholar
- McCarthy JJ, Intergovernmental Panel on Climate Change. Working Group II.: Climate change 2001 : impacts, adaptation, and vulnerability : contribution of Working Group II to the third assessment report of the Intergovernmental Panel on Climate Change. 2001, Cambridge, UK ; New York, Cambridge University Press, x, 1032-Google Scholar
- Wilby RL, Charles SP, Zorita E, Timbal B, Whetton P, Mearnes LO: Guidelines for use of climate scenarios developed from statistical downscaling methods. [http://ipcc-ddc.cru.uea.ac.uk/guidelines/dgm_no2_v1_09_2004.pdf]
- Website: Quantum GIS - An open source geographic information system. [http://www.qgis.org/]
- Maddison WP, Maddison DR: Mesquite: A modular system for evolutionary analysis. Version 1.05. http://mesquiteproject.org. 2004, [http://mesquiteproject.org]Google Scholar
- Website: openModeller - An open source spatial distribution modelling library. [http://openmodeller.sourceforge.net/]
- Website: Maxent software and datasets. [http://www.cs.princeton.edu/~schapire/maxent/]
- Ackerly DD, Donoghue MJ: Leaf size, sapling allometry, and Corner's rules: Phylogeny and correlated evolution in maples (Acer). Am Nat. 1998, 152: 767-791. 10.1086/286208.View ArticlePubMedGoogle Scholar
- Prinzing A, Durka W, Klotz S, Brandl R: The niche of higher plants: evidence for phylogenetic conservatism. Proc R Soc Lond Ser B-Biol Sci. 2001, 268: 2383-2389. 10.1098/rspb.2001.1801.View ArticleGoogle Scholar
- Segurado P, Araújo MB: An evaluation of methods for modelling species distributions. J Biogeogr. 2004, 31: 1555-1568. 10.1111/j.1365-2699.2004.01076.x.View ArticleGoogle Scholar
- Buckley LB, Roughgarden J: Biodiversity conservation - Effects of changes in climate and land use. Nature. 2004, 430: U2-U2. 10.1038/nature02717.View ArticleGoogle Scholar
- Pearson RG: Climate change and the migration capacity of species. Trends Ecol Evol. 2006, 21: 111-113. 10.1016/j.tree.2005.11.022.View ArticlePubMedGoogle Scholar
- Davis MB, Shaw RG: Range shifts and adaptive responses to Quaternary climate change. Science. 2001, 292: 673-679. 10.1126/science.292.5517.673.View ArticlePubMedGoogle Scholar
- McLachlan JS, Clark JS, Manos PS: Molecular indicators of tree migration capacity under rapid climate change. Ecology. 2005, 86: 2088-2098.View ArticleGoogle Scholar
- IUCN. Species Survival Commission.: IUCN red list categories and criteria. 2001, Gland ; Cambridge, IUCN, 30-3.1.Google Scholar
This article is published under license to BioMed Central Ltd. This is an Open Access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/2.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.