Skip to main content
Advertisement
Browse Subject Areas
?

Click through the PLOS taxonomy to find articles in your field.

For more information about PLOS Subject Areas, click here.

  • Loading metrics

The effect of global warming on the Australian endemic orchid Cryptostylis leptochila and its pollinator

  • Marta Kolanowska ,

    Roles Conceptualization, Data curation, Formal analysis, Funding acquisition, Investigation, Methodology, Project administration, Resources, Software, Supervision, Validation, Visualization, Writing – original draft, Writing – review & editing

    martakolanowska@wp.pl

    Affiliation Faculty of Biology and Environmental Protection, Department of Geobotany and Plant Ecology, University of Lodz, Poland

  • Ewa Michalska

    Roles Data curation, Writing – review & editing

    Affiliation Faculty of Biology and Environmental Protection, Department of Geobotany and Plant Ecology, University of Lodz, Poland

Abstract

Ecological stability together with the suitability of abiotic conditions are crucial for long-term survival of any organism and the maintenance of biodiversity and self-sustainable ecosystems relies on species interactions. By influencing resource availability plants affect the composition of plant communities and ultimately ecosystem functioning. Plant-animal interactions are very complex and include a variety of exploitative and mutualistic relationships. One of the most important mutualistic interactions is that between plants and their pollinators. Coevolution generates clustered links between plants and their pollen vectors, but the pollination and reproductive success of plants is reduced by increase in the specialization of plant-animal interactions. One of the most specialized types of pollination is sexual deception, which occurs almost exclusively in Orchidaceae. In this form of mimicry, male insects are attracted to orchid flowers by chemical compounds that resemble insect female sex pheromones and pollinate the flowers during attempted copulations. These interactions are often species-specific with each species of orchid attracting only males of one or very few closely related species of insects. For sexually deceptive orchids the presence of a particular pollen vector is crucial for reproductive success and any reduction in pollinator availability constitutes a threat to the orchid. Because global warming is rapidly becoming the greatest threat to all organisms by re-shaping the geographical ranges of plants, animals and fungi, this paper focuses on predicting the effect of global warming on Cryptostylis leptochila, a terrestrial endemic in eastern Australia that is pollinated exclusively via pseudo copulation with Lissopimpla excelsa. As a species with a single pollinator this orchid is a perfect model for studies on the effect of global warming on plants and their pollen vectors. According to our predictions, global warming will cause a significant loss of suitable niches for C. leptochila. The potential range of this orchid will be 36%-75% smaller than currently and as a result the Eastern Highlands will become unsuitable for C. leptochila. On the other hand, some new niches will become available for this species in Tasmania. Simultaneously, climate change will result in a substantial expansion of niches suitable for the pollinator (44–82%). Currently ca. 71% of the geographical range of the orchid is also suitable for L. excelsa, therefore, almost 30% of the areas occupied by C. leptochila already lack the pollen vector. The predicted availability of the pollen vector increased under three of the climate change scenarios analysed. The predicted habitat loss is a serious threat to this orchid even with the potential colonization of Tasmania by this plant. In the reduced range of C. leptochila the pollen vector will also be present assuring fruit set in populations of this orchid. The genetic pool of the populations in New South Wales and Queensland will probably be lost.

Introduction

Species interactions are crucial for maintaining natural biodiversity and are essential for preserving self-sustainable ecosystems [1, 2]. By affecting resource availability and habitat structure plants affect plant communities and ultimately ecosystem functioning [3, 4]. Plant-plant interactions and their role in mediating the effect of environmental transformation are discussed by Brooker [5]. Facilitation among plants is considered to be a significant driver of ecosystem structure and functioning, especially in regions characterized by harsh environmental conditions [6]. On the other hand, competition for resources (water, nutrient, sunlight) not only influences the composition of communities it also has a role in shaping the distribution of species, as well as their evolution [7]. Plant-fungi interactions can also be beneficial or detrimental for plants. Fungal pathogens negatively affect plant physiology, whereas mutualistic fungi may enhance host defence against pathogens and symbiotic fungi can also improve plant nutrient uptake [8]. Plant-animal interactions are very complex and include a variety of biological relationships. Among these, plant-herbivore relations are examples of exploitation, whereas mutualistic connections include plants providing nesting sites and/or food rewards to ants, which often protect plants from herbivores or competing plants. Plant-pollinator and plant-seed disperser interactions are usually beneficial for both partners. The process of pollen transfer most often involves insects, birds or bats [9]. Noteworthy, coevolution generates clustered links between plants and their pollinators or seed dispersers [10], but the stability of pollination and reproductive success of plants is reduced by increased specialization of the plant-animal interaction [1113].

One of the most specialized types of pollination is sexual deception (so-called “pseudo copulation”) which occurs almost exclusively in Orchidaceae [14, 15]. In this Pouyannian mimicry plants lure their pollinators with fraudulent sex signals [16, 17] and the pollen transfer requires specific visual and chemical adaptations of the plant. Species utilizing sexual deception depend on a single or very few species of insects [1823].

Australia is a centre of diversity of plants pollinated via sexual deception, with 11 orchid genera luring insects with false sex signals [16]. These plants are found in a variety of habitats: woodlands (27%), shrubland (16%), forests with grassy or shrubby bushes (12%), damp spots near streams (11%) and exposed rocky outcrops (10%) [24]. Most species can grow in several habitats [16, 2531]. Unfortunately, numerous Australian orchids are listed as extinct, threatened or rare [32] and this catalogue includes numerous taxa pollinated via sexual deception, e.g. Caladenia (>80 species), Chiloglottis R. Br. (7 species), Cryptostylis R. Br. (4 species), Drakaea Lindl. (5 species) and Pterostylis R. Br. (> 40 species).

Cryptostylis leptochila F. Muell. ex Benth. the small tongue orchid is an Australian endemic. It is characterized by green leaves the under surface of which is reddish-purple. Plants produce 5–15 flowers with green, linear sepals and petals and a dark red, oblong-linear lip. The lip of this species is densely covered with glandular hairs and the callus forms a ridge, which is often discontinuous, flanked on either side by an irregular row of shiny black calli [33]. C. leptochila is only pollinated by Lissopimpla excelsa Costa [Pimplinae, Ichneumonidae], the so called orchid dupe wasp [34, 35]. The male wasp faces the top of the upright orchid lip and its abdomen and genital claspers probe the curved pocket at base of the lip, where the stigma and the pollinia are located [25]. This plant-insect relationship is unique in being the only known act of pollination that elicits the pollinator to ejaculate [36]. Moreover, all Australian representatives of Cryptostylis share Lissopimpla excelsa as a pollen vector [25].

Pseudo copulation between an orchid and specific species of insect needs further study because global warming [37, 38] constitutes a serious threat to global biodiversity [3951] and ecological interactions [5257]. Many studies indicate that frequent, prolonged droughts adversely affect the pollination and seed dispersal of plants, which are crucial for reproduction [58, 59]. Orchids are particularly at risk [6068] because global warming adversely affects the specific symbiotic associations of these plants with fungi and their complex relationships with pollinators. Any reduction in the availability of pollen vectors are especially dangerous for sexually deceptive plants whose long-term persistence [fruit production] depends on particular species of insects.

Undeniably, suitable environmental conditions together with the stability of a series of interactions critical for plant growth, development and reproductive success are essential for orchid survival [66]. The aim of this study was to evaluate the availability of the specific pollinator of C. leptochila in various scenarios of climate change using ecological niche modelling.

Methods

List of localities

The databases of the localities of Cryptostylis leptochila and Lissopimpla excelsa were compiled based on information in public facilities accessed through Atlas of Living Australia [69], which includes data from Shire of Kalamunda Biodiversity Inventor, Australia’s Virtual Herbarium, Australian National Insect Collection, BowerBird, Tasmanian Museum and Art Gallery Invertebrate Collection, Questagame, Queen Victoria Museum and Art Gallery and South Australian Museum, as well as iNaturalist.org and Global Biodiversity Information Facility [70, 71].

From the total of 1442 localities for C. leptochila and 856 for L. excelsa in the repositories only records that were georeferenced with a minimum of 1 km precision were selected (S1 Annex).

To reduce sampling bias spatial thinning was done using SDMtoolbox 2.3 for ArcGIS [72, 73]. The location data were spatially filtered at 5 km2, 10 km2, 15 km2, 20 km2 and 25 km2, in areas of high, medium-high, medium, medium-low and low habitat heterogeneity, respectively. This method of graduated filtering was used to maximize the number of spatially independent localities. The final database included 210 records for C. leptochila and 161 for L. excelsa (S2 Annex).

Climatic niche modelling

Species distribution and ecological niche modelling are methods commonly used to predict the potential distribution of suitable habitats for various organisms [59, 7478]. The modelling of the current and future distribution of the species studied was done using the maximum entropy method implemented in MaxEnt version 3.3.2 [7981], which is based on presence-only observations. Bioclimatic variables in 30 arc-seconds of interpolated climate surface downloaded from WorldClim v. 2.1 were used for the modelling [82]. Five of 19 variables were removed from the analyses due to high correlation between them (>0.9) as indicated by Pearsons’ correlation coefficient computed using SDMtoolbox 2.3 for ArcGIS [72, 73] in order to avoid problems associated with auto-correlation. The final list of bioclimatic variables used in the analyses is provided in S1 Table. Because some previous studies [83] suggest that modelling for a small area is more reliable than calculating habitat suitability at a global scale, the area of the analysis was restricted to 5.454–44.838°S and 111.079° -156.921°E.

Predictions of the future extent of the climatic niches of C. leptochila and its pollinator in 2080–2100 were made using climate projections developed by the CNRM/CERFACS modelling group for Coupled Model Intercomparison Project (CNRM-CM6-1) for four Shared Socio-economic Pathways (SSPs): 1–2.6, 2–4.5, 3–7.0 and 5–8.5 [8486]. The layers in 2.5 arc-minutes were re-scaled to fit bioclimatic variables. SSPs are trajectories adopted by the Intergovernmental Panel on Climate Change (IPCC), which offer a broader view of a “business as usual” world without future climate policy, with global warming in 2100 ranging from a low of 3.1°C to a high of 5.1°C above pre-industrial levels.

Because projections of future changes in the distribution and coverage of suitable niches of the orchid studied is based exclusively on bioclimatic data, the importance of other factors such as soil properties and land cover were evaluated by creating present-time models using all available variables. The models created for present time based on bioclimatic data only and based on combined bioclimatic, soil and land cover data were then compared with the known occurrences of C. leptochila. The data on the distribution of soil classes, soil pH, clay content, sand content, nitrogen content and soil organic carbon content were obtained from Global Soil Information [87] (http://www.soilgrids.org) with a 250 m2 resolution and upscaled to fit the resolution and extent of the bioclimatic variables. The Global Land Cover by National Mapping Organizations ver. 3 (GLCNMO) was the source of information for the land cover analyses. The complete list of variables used in the present time analyses is presented in S1 Table.

In all analyses the maximum number of iterations was set to 10000 and convergence threshold to 0.00001. The neutral (= 1) regularization multiplier value and auto features were used. All samples were added to the background. The “random seed" option provided a random test partition and background subset for each run and 30% of the samples were used as test points. The run was performed as a bootstrap with 100 replicates. The output was set to logistic. In addition, the “fade by clamping” function in MaxEnt precluded extrapolations outside the environmental range of the training data [88]. All analyses of GIS data were carried out using ArcGis 10.6 (Esri, Redlands, CA, USA). The evaluation of the created models was done using the area under the curve (AUC) [89, 90] and True Skill Statistic (TSS) [91, 92].

SDMtoolbox 2.3 for ArcGIS [72, 73] was used to visualize changes in the distribution of suitable niches of the orchid studied and its pollinator due to global warming. To compare the distribution model created for current climatic conditions with future predictions all SDMs were converted into binary rasters and projected using the Goode homolosine as a projection. The presence threshold was estimated based on the values for grids in which the species studied occur in models developed using current data. About 80% of known localities for C. leptochila are located in grids with values > 0.4 and 78% of those for L. excelsa in grids with values > 0.3. These two thresholds were used to create binary rasters [93]. To calculate the availability of the pollinator for the orchid studied the overlap in the binary models of both organisms was calculated.

Results

Model evaluation and limiting factors

The scores of the values of both, TSS and AUC tests were high with 0.780–0.976 and 0.965–0.990, respectively (S2 Table), which indicates the high accuracy of the models.

The crucial factor limiting the distribution of C. leptochila in models based on bioclimatic data only was precipitation in the driest month (bio14–67.3%), much less important was the maximum temperature in the warmest month (bio5–10.6%) and the mean diurnal range (bio2–10.5%). For L. excelsa essential variables were the precipitation in coldest quarter (bio19–37.4%) and precipitation in the driest month (bio14–25.7%). A somewhat smaller contribution was recorded for the mean diurnal range (bio2–12.0%). The effect of environmental variables on the MaxEnt prediction is presented in Fig 1. Models created for the orchid studied based on all available variables and produced using only bioclimatic data are very similar, with a slight overestimation of the potential range predicted by the model including all variables (Fig 2). Also, the jackknife test of variable importance indicated that bioclimatic data are crucial for constructing models (Fig 1).

thumbnail
Fig 1.

A. Effect of environmental variables on the MaxEnt prediction based only on bioclimatic data. The curves show how the predicted probability of occurrence changes as each environmental variable is varied, keeping all other environmental variables at their average value. The curves show the mean response of the 100 replicate Maxent runs (red) and the mean +/- one standard deviation (blue). B. Results of the jack-knife test of variable importance for the model created for C. leptochila using all of the variables. Values shown are averages for replicate runs.

https://doi.org/10.1371/journal.pone.0280922.g001

thumbnail
Fig 2. Current potential distribution of Cryptostylis leptochila (A-C) and Lissopimpla excelsa (D).

A—model based on all variables and all georeferenced records of C. leptochila, B—map showing the localities of C. leptochila used in ENM studies, C-D—model based only on bioclimatic data with all georeferenced records for the species studied. Maps generated in ArcGIS based on niche modelling (MaxEnt) results.

https://doi.org/10.1371/journal.pone.0280922.g002

Current and potential future distribution of orchid and its pollinator

The models of current distribution of niches suitable for C. leptochila and L. excelsa are consistent with the known geographical ranges of these species (Fig 2). Based on these models global warming will result in a significant loss of suitable niches for C. leptochila (Figs 3 and 4; Table 1), but it will be favourable for L. excelsa (Figs 5 and 6, Table 1). The potential range of the orchid will be 36.064%(SSP1-2.6)-75.970% (SSP5-8.5) smaller than it is currently. In the SSP1-2.6 scenario habitat loss is expected to occur along the Eastern Highlands in New South Wales, but habitats will become available in north-east Tasmania, the Furneaux Group and some localities near Mt. Kościuszko. A similar situation is predicted by SSP2-4.5 and SSP3-7.0 models, but the habitat loss in the Eastern Highlands will be much more severe and will occur in both New South Wales and Queensland. In the SSP3-7.0 scenario, however, numerous new niches will become available for this species in southern- and northern- eastern Tasmania. In the SSP5-8.5 scenario almost all suitable niches in the Eastern Highlands will cease to exist, except for the Tasman Sea coast near the Croajingolong National Park.

thumbnail
Fig 3. Predicted future distribution of suitable niches for C. leptochila based on the SSP1-2.6 (A), SSP2-4.5 (B), SSP3-7.0 (C) and SSP5-8.5 (D) climate change scenarios.

Maps generated in ArcGIS based on niche modelling (MaxEnt) results.

https://doi.org/10.1371/journal.pone.0280922.g003

thumbnail
Fig 4. Predicted changes in the distribution of suitable niches for Cryptostylis leptochila based on the SSP1-2.6 (A), SSP2-4.5 (B), SSP3-7.0 (C) and SSP5-8.5 (D) climate change scenarios.

Maps generated in ArcGIS based on niche modelling (MaxEnt) results.

https://doi.org/10.1371/journal.pone.0280922.g004

thumbnail
Fig 5. Predicted future distribution of suitable niches for Lissopimpla excelsa based on the SSP1-2.6 (A), SSP2-4.5 (B), SSP3-7.0 (C) and SSP5-8.5 (D) climate change scenarios.

Maps generated in ArcGIS based on niche modelling (MaxEnt) results.

https://doi.org/10.1371/journal.pone.0280922.g005

thumbnail
Fig 6. Predicted changes in the distribution of suitable niches for Lissopimpla excelsa based on SSP1-2.6 (A), SSP2-4.5 (B), SSP3-7.0 (C) and SSP5-8.5 (D) climate change scenarios.

Maps generated in ArcGIS based on niche modelling (MaxEnt) results.

https://doi.org/10.1371/journal.pone.0280922.g006

thumbnail
Table 1. Changes in the coverage (km2) of suitable niches for C. leptochila and L. excelsa.

https://doi.org/10.1371/journal.pone.0280922.t001

On the other hand, global warming will result in a substantial increase in the coverage of niches suitable for L. excelsa. Based on SSP3-7.0, which is the most advantageous scenario for this species, the potential coverage will be 82.109% greater than currently. Based on the worst-case scenario (SSP1-2.6), the coverage of suitable niches for this species will increase to 44.103%. The habitat loss for L. excelsa predicted by all the models due to global warming is not significant and will occur near the Byfield State Forest. The increase in the coverage of suitable niches is expected to occur along the edges of the current potential range and in central and southern Tasmania, south part of Spencer Gulf and around Cape Leeuwin.

Currently most of the populations of this orchid are likely to be pollinated as 71.43% of this orchid’s geographical range is also suitable for L. excelsa (Fig 7). The availability of the pollen vector is predicted to increase in three climate change scenarios: SSP1-2.6 (83.33%), SSP3-7.0 (72.73%) and SSP5-8.5 (90%). In SSP2-4.5 part of the range of the orchid will overlap with the potential range of the insect.

thumbnail
Fig 7. Overlap of suitable niches for Cryptostylis leptochila and Lissopimpla excelsa currently (A) and based on the SSP1-2.6 (B), SSP2-4.5 (C), SSP3-7.0 (D) and SSP5-8.5 (E) climate change scenarios.

Maps generated in ArcGIS based on niche modelling (MaxEnt) results.

https://doi.org/10.1371/journal.pone.0280922.g007

Discussion

The ecology of specialized pollination is very complex. In this study the possible effect of climate change on the distribution of a sexually deceptive orchid and its pollinator were explored, which revealed that the reproductive success of C. leptochila depends on many ecological interactions. The pollinator of C. leptochila is a parasitoid that is dependent on the availability of a host insects that feed on particular plants. In addition, germination of the orchid is dependent on symbiosis with mycorrhizal fungi, which are also often dependent on the occurrence of other species of plants. Here we discuss the results of our analysis and limitations of this study, which requires a more detailed insight into the ecological interactions between plants, insects and fungi.

Areas with low and very high values for precipitation in the driest month are unsuitable for C. leptochila and since global warming will significantly change the water regimes in the world it is not surprising that numerous currently suitable niches for this orchid will change and probably become unsuitable for this species. The two other variables important for the occurrence of this species are related to temperature. Cryptostylis leptochila prefers areas where the maximum temperature in the warmest month is ca. 25°C and the mean diurnal range is ca. 5°C. The predicted climate changes [94] indicate that Australia will be hotter, with longer and more intense heatwaves. As a result of global warming the Eastern Highlands will become unsuitable for C. leptochila, however some new niches will become available for this species in Tasmania. The predicted loss in suitable niches of 36–76% is a serious threat to the orchid even with its potential expansion in Tasmania. The genetic pool of the populations in New South Wales and Queensland will probably be lost.

In contrast to previous studies on the future of orchids pollinated via pseudo copulation with a single species of insect [37] this research predicts that the orchid studied will face significant habitat loss, which will be accompanied by an expansion in the geographical range of its pollen vector. There is several studies on the European sexually deceptive genus Ophrys L. which is currently the most diverse taxon comprising of orchids pollinated via pseudo copulation [95]. Hutchings et al. [96] as well as Robbirt et al. [97] studied only one of five known pollinators of O. sphegodes Mill. They provide proof of disruption of the phenology and future interactions between this orchid and a bee (Andrena nigroaenea Kirby) and the potential effect of climate change on another three Hymeoptera (Andrena nitida Müller, Andrena ovatula Kirby, Andrena paucisquama Noskiewicz) and one Diptera (Myopa tessallatipennis Motschulsky) that are pollen vectors for the early spider-orchid [98101]. Tsiftis & Djordjević [102] include the presence of the exclusive pollinator of Greek Ophrys argolica Fleischm. ex Vierh. and O. delphinensis O. Danesch & E. Danesch as a variable in their models of future occurrence of these orchids. This approach, however, is misleading as the vegetative growth of bee-orchid is not dependent on the existence of a pollen vector. The methodology used by Tsiftis & Djordjević [102] is more applicable to studies on parasitic plants, which are entirely dependent on the presence of their hosts.

Our predictions indicate Tasmania as a future potentially suitable area for C. leptichila and L. excelsa. The orchid has the potential to migrate into this region as it is already present in mainland Australia as well as in Furneaux Group islands. The long-distance dispersal [> 150 km] is not uncommon in terrestrial orchids [103, 104]. Orchid dupe wasp is already present in New Zealand [105] so its migration to Tasmania is also probable.

In the limited number of niches that will still be available for C. leptichila the pollen vector will be present assuring fruit set in orchid populations. There is only one long-term study on a species of orchid pollinated via sexual deception [96, 97] and on one of at least five known pollinators of this taxon. This study indicates that ongoing climate warming, especially warmer springs, is likely to increase the frequency of years in which the plant experiences low pollination success. According to information in public databases over the last 70 years the flowering time of C. leptochila has not changed despite climate change. While few occurrences of L. excelsa were reported before 1980, based on the public databases records, the pollinator is available during the whole flowering period of this orchid (S3 Table). On the other hand, based on museum records Brunton Martin et al. [106] demonstrate that male orchid dupe wasps can be less abundant between May and September so that the synchronization of the flowering of the orchid and male insect activity could be threatened by climate change modifying their phenology. Interestingly, the same researchers [106] also show that populations of L. excelsa are male biased in the areas where C. leptochila occurs allowing the persistence of the pollen vector despite sperm wastage during pseudo copulation. Though global warming can potentially cause modifications in insect behaviour, the daily activity of L. excelsa on C. leptochila flowers currently does not seem to be limited by any particular climatic factor like humidity or sunlight [107, 108]. Male L. excelsa remain active over a broad range of temperatures [109] (20.1–36.8°) and based on the maps of future climatic conditions [82] the temperatures in localities where C. leptochila occurs will not exceed that tolerated by these insects.

This study explored the distribution of an orchid and its pollinator, but not the long-term survival of the orchid dupe wasp, which is also threatened by the effects of global warming. Females of L. excelsa parasitize large caterpillars such as budworms and armyworms [110] and the stability of orchid dupe wasp population depends on the availability of these larvae. For Lepidoptera one of the major threats associated with climate change are temporal mismatches with availability of food sources [111]. Global warming affects the development of insect herbivores more than their hosts leading to desynchronization of flowering time and insect activity period [112]. Increased temperatures trigger also asynchrony by supporting faster development of caterpillars than their parasitoids [111, 113]. The complex relationship of L. excelsa with the species it parasitizes could not be included in this study because the list of Lepidoptera species parasitized by orchid dupe wasp is incomplete. Also, the list of such Lepidoptera would have to be compiled and modelled to estimate their chances of survival when faced by changes in climate.

Obviously, pollinator availability is not the only ecological factor that could change due to global warming becoming a threat to C. leptochila. Orchids have numerous symbiotic relationships with endophytic fungi [114, 115] and bacteria [116]. The interaction with a mycorrhizal fungus, which provides essential nutrients for seedling development is vital for orchids [117] and the seedling recruitment success is strongly associated with the availability of a suitable fungal strain. While some orchids are mycorrhizal generalists [118], other are characterized by specific relationships with fungal genera [119]. The nature of the symbiotic relationships of C. leptochila with fungi and bacteria are poorly known and as a consequence this aspect of orchid ecology could not be evaluated in this study. In numerous mycological studies there is no species-level identification of symbiotic fungi only taxonomic units (OTUs) [120]. Without proper naming of orchid symbionts, a dataset of the distributions of particular fungi could not be created and included in the ENM analyses.

Abiotic factors are also likely to affect the distributions of orchids as germination of their seeds is dependent on soil moisture, pH and organic content [121]. Also, extreme events like heat waves, hurricanes, etc, are expected to occur more frequently as a result of climate change [122], which could substantially affect orchid populations. According to Brown [123] apart from the factors listed above, inappropriate fire regimes, groundwater extraction, spread of exotic plants and grazing should be considered as potential threats to the survival of Australian orchids. Unfortunately, because it is not possible to predict how the variables listed above will change in various scenarios of climate change we could not incorporate these factors in our analyses.

To summarize, the analyses revealed that climate change will adversely affect sexually deceptive orchids, but at the same time the pollinator of C. leptochila will benefit from global warming. In the worst-case scenario this orchid will lose 75% of its currently suitable niches, but even in its very limited geographical range the availability of the pollen vector will be higher than currently. The ecological interactions and climatic variables that could not be modelled using available data can be potential additional threats to the orchid studied.

Supporting information

S1 Annex. Records of Cryptostylis leptochila and Lissopimpla excelsa georeferenced with a minimum precision of 1 km.

https://doi.org/10.1371/journal.pone.0280922.s001

(XLSX)

S2 Annex. Records of Cryptostylis leptochila and Lissopimpla excelsa included in the ENM.

https://doi.org/10.1371/journal.pone.0280922.s002

(XLSX)

S1 Table. List of variables used in the modelling.

https://doi.org/10.1371/journal.pone.0280922.s003

(DOCX)

S2 Table. Results of the statistical evaluation of the models.

https://doi.org/10.1371/journal.pone.0280922.s004

(DOCX)

S3 Table. Calendar of the occurrence the flowering of C. leptochila and activity of L. excelsa in four periods of time (1951–1980, 1981–2000, 2001–2010, 2011–2020) compiled based on public databases.

https://doi.org/10.1371/journal.pone.0280922.s005

(DOCX)

Acknowledgments

We are grateful to prof. Anthony Dixon for the language correction of the manuscript. We would like to thank anonymous Reviewers for their valuable comments which helped us to improve the paper.

References

  1. 1. Bairey E, Kelsic ED, Kishony R. High-order species interactions shape ecosystem diversity. Nature Communications. 2016;7(1):12285. pmid:27481625
  2. 2. Allesina S, Tang S. Stability criteria for complex ecosystems. Nature. 2012;483(7388):205–8. pmid:22343894
  3. 3. Martorell C, Freckleton RP. Testing the roles of competition, facilitation and stochasticity on community structure in a species-rich assemblage. Journal of Ecology. 2014;102(1):74–85.
  4. 4. Grime JP. Benefits of plant diversity to ecosystems: immediate, filter and founder effects. Journal of Ecology. 1998;86(6):902–10.
  5. 5. Brooker RW. Plant–plant interactions and environmental change. New Phytologist. 2006;171(2):271–84. pmid:16866935
  6. 6. Brooker RW, Maestre FT, Callaway RM, Lortie CL, Cavieres LA, Kunstler G, et al. Facilitation in plant communities: the past, the present, and the future. Journal of Ecology. 2008;96(1):18–34.
  7. 7. Craine JM, Dybzinski R. Mechanisms of plant competition for nutrients, water and light. Functional Ecology. 2013;27(4):833–40.
  8. 8. Zeilinger S, Gupta VK, Dahms TE, Silva RN, Singh HB, Upadhyay RS, et al. Friends or foes? Emerging insights from fungal interactions with plants. FEMS Microbiol Rev. 2016;40(2):182–207. pmid:26591004
  9. 9. Phatlane WM. Pollination in Plants. London: IntechOpen; 2018.
  10. 10. Waser NM, Ollerton J. Plant-Pollinator Interactions: From Specialization to Generalization. Chicago: University of Chicago Press; 2005.
  11. 11. Wolf LL, Stiles FG. Adaptations for the ’Fail-safe’ Pollination of Specialized Ornithophilous Flowers. The American Midland Naturalist. 1989;121(1):1–10.
  12. 12. Johnson SD, Steiner KE. Generalization versus specialization in plant pollination systems. Trends in Ecology & Evolution. 2000;15(4):140–3. pmid:10717682
  13. 13. Armbruster WS. The specialization continuum in pollination systems: diversity of concepts and implications for ecology, evolution and conservation. Functional Ecology. 2017;31(1):88–100.
  14. 14. Scopece G, Cozzolino S, Johnson S, Schiestl F. Pollination Efficiency and the Evolution of Specialized Deceptive Pollination Systems. American Naturalist. 2010;175(1):98–105. pmid:19909087
  15. 15. Ellis AG, Johnson SD. Floral mimicry enhances pollen export: the evolution of pollination by sexual deceit outside of the orchidaceae. Am Nat. 2010;176(5):E143–51. pmid:20843263
  16. 16. Gaskett AC. Orchid pollination by sexual deception: pollinator perspectives. Biol Rev Camb Philos Soc. 2011;86(1):33–75. pmid:20377574
  17. 17. Vereecken NJ, Schiestl FP. The evolution of imperfect floral mimicry. Proc Natl Acad Sci U S A. 2008;105(21):7484–8. pmid:18508972
  18. 18. Ames O, Ames B. Pollination of orchids through pseudocopulation. Botanical Museum Leaflets, Harvard University. 1937;5(1):1–xix.
  19. 19. Schiestl FP. Pollination: Sexual Mimicry Abounds. Current Biology. 2010;20(23):R1020–R2. pmid:21145017
  20. 20. Vereecken N, Wilson C, Hotling S, Schulz S, Banketov S, Mardulyn P. Pre-adaptations and the evolution of pollination by sexual deception: Cope’s rule of specialization revisited. Proceedings of the Royal Society B-Biological Sciences. 2012;279(1748):4786–94. pmid:23055065
  21. 21. Martel C, Francke W, Ayasse M. The chemical and visual bases of the pollination of the Neotropical sexually deceptive orchid Telipogon peruvianus. New Phytol. 2019;223(4):1989–2001. pmid:31074029
  22. 22. Ayasse M, Stokl J, Francke W. Chemical ecology and pollinator-driven speciation in sexually deceptive orchids. Phytochemistry. 2011;72(13):1667–77. pmid:21497864
  23. 23. Peakall R, Whitehead MR. Floral odour chemistry defines species boundaries and underpins strong reproductive isolation in sexually deceptive orchids. Ann Bot. 2014;113(2):341–55. pmid:24052555
  24. 24. Rutgrink A. Sexual deception in flowering plants (PhD Thesis). Amsterdam: VU University Amsterdam; 2016.
  25. 25. Gaskett A. Floral shape mimicry and variation in sexually deceptive orchids with a shared pollinator. Biological Journal of the Linnean Society. 2012;106(3):469–81.
  26. 26. Boyden TC. The pollination biology of Calypso bulbosa var. Americana (Orchidaceae): Initial deception of bumblebee visitors. Oecologia. 1982;55(2):178–84. pmid:28311231
  27. 27. Schiestl F. On the success of a swindle: pollination by deception in orchids. Naturwissenschaften. 2005;92(6):255–64. pmid:15931514
  28. 28. Johnson S, Steiner K, Kaiser R. Deceptive pollination in two subspecies of Disa spathulata (Orchidaceae) differing in morphology and floral fragrance. Plant Systematics and Evolution. 2005;255(1–2):87–98.
  29. 29. Stokl J, Paulus H, Dafni A, Schulz C, Francke W, Ayasse M. Pollinator attracting odour signals in sexually deceptive orchids of the Ophrys fusca group. Plant Systematics and Evolution. 2005;254(1–2):105–20.
  30. 30. Cozzolino S, Scopece G, Lussu M, Cortis P, Schiestl FP. Do floral and ecogeographic isolation allow the co-occurrence of two ecotypes of. Ecol Evol. 2021;11(15):9917–31.
  31. 31. Cohen C, Liltved WR, Colville JF, Shuttleworth A, Weissflog J, Svatoš A, et al. Sexual deception of a beetle pollinator through floral mimicry. Curr Biol. 2021;31(9):1962–9.e6. pmid:33770493
  32. 32. Backhouse GN. Are our orchids safe down under? A national assessment of threatened orchids in Australia. Lankesteriana. 2007;7:28–43.
  33. 33. Gaskett AC, Herberstein ME. Colour mimicry and sexual deception by Tongue orchids (Cryptostylis). Naturwissenschaften. 2010;97(1):97–102. pmid:19798479
  34. 34. Coleman E. Pollination of an Australian orchid by the male ichneumonid Lissopimpla semipunctata, Kirby. Transactions of the Entomological Society of London. 1928;2:533–9.
  35. 35. Coleman E. Further observations on the pseudocopulation of the male Lissopimpla semipunctata Kirby (Hymenoptera parasitica) with the Australian orchid Cryptostylis leptochila. Proceedings of the Royal Entomological Society of London. 1938;13:82–4.
  36. 36. Brunton-Martin AL, Gaskett AC, Kokko H. Resilience of haplodiploids to being exploited by sexually deceptive plants. Oikos. 2021;130(11):2053–63.
  37. 37. Kolanowska M, Michalska E, Konowalik K. The impact of global warming on the niches and pollinator availability of sexually deceptive orchid with a single pollen vector. Sci Total Environ. 2021;795:148850. pmid:34246141
  38. 38. Tsiftsis S, Djordjevic V. Modelling sexually deceptive orchid species distributions under future climates: the importance of plant-pollinator interactions. Sci Rep. 2020;10(1):10623. pmid:32606363
  39. 39. Alabia I, Molinos J, Saitoh S, Hirata T, Hirawake T, Mueter F. Multiple facets of marine biodiversity in the Pacific Arctic under future climate. Science of the Total Environment. 2020;744. pmid:32721679
  40. 40. Arribas P, Abellan P, Velasco J, Bilton D, Millan A, Sanchez-Fernandez D. Evaluating drivers of vulnerability to climate change: a guide for insect conservation strategies. Global Change Biology. 2012;18(7):2135–46.
  41. 41. Bateman B, Taylor L, Wilsey C, Wu J, LeBaron G, Langham G. Risk to North American birds from climate change-related threats. Conservation Science and Practice. 2020;2(8).
  42. 42. Berdugo M, Quant J, Wason J, Dovciak M. Latitudinal patterns and environmental drivers of moss layer cover in extratropical forests. Global Ecology and Biogeography. 2018;27(10):1213–24.
  43. 43. Botkin D, Saxe H, Araujo M, Betts R, Bradshaw R, Cedhagen T, et al. Forecasting the effects of global warming on biodiversity. Bioscience. 2007;57(3):227–36.
  44. 44. Brodie J, Post E, Laurance WF. Climate change and tropical biodiversity: a new focus. Trends Ecol Evol. 2012;27(3):145–50. pmid:21975172
  45. 45. Descombes P, Wisz M, Leprieur F, Parravicini V, Heine C, Olsen S, et al. Forecasted coral reef decline in marine biodiversity hotspots under climate change. Global Change Biology. 2015;21(7):2479–87. pmid:25611594
  46. 46. Fordham D, Wigley T, Brook B. Multi-model climate projections for biodiversity risk assessments. Ecological Applications. 2011;21(8):3317–31.
  47. 47. Geppert C, Perazza G, Wilson RJ, Bertolli A, Prosser F, Melchiori G, et al. Consistent population declines but idiosyncratic range shifts in Alpine orchids under global change. Nat Commun. 2020;11(1):5835. pmid:33203870
  48. 48. Kasemsap K, Gezerman A, Corbacioglu B, Gurjar B, Gezerman A, Corbacioglu B, et al. Global Warming and Climate Change: Challenges and Impacts. Effective Solutions To Pollution Mitigation For Public Welfare. 2018:44–68.
  49. 49. McLaughlin J, Hellmann J, Boggs C, Ehrlich P. Climate change hastens population extinctions. Proceedings of the National Academy of Sciences of the United States of America. 2002;99(9):6070–4. pmid:11972020
  50. 50. Schwartz M, Iverson L, Prasad A, Matthews S, O’Connor R. Predicting extinctions as a result of climate change. Ecology. 2006;87(7):1611–5. pmid:16922312
  51. 51. Villen-Perez S, Heikkinen J, Salemaa M, Makipaa R. Global warming will affect the maximum potential abundance of boreal plant species. Ecography. 2020;43(6):801–11.
  52. 52. Buckley L, Angert A, LaDeau S, Ostfeld R. Get real: putting models of climate change and species interactions in practice. Climate Change and Species Interactions: Ways Forward. 2013;1297:126–38. pmid:23796014
  53. 53. Spence A, Tingley M. The challenge of novel abiotic conditions for species undergoing climate-induced range shifts. Ecography. 2020;43(11):1571–90.
  54. 54. Williams MA, Holland CV, Donohue I. Warming can alter host behavior in a similar manner to infection with behavior-manipulating parasites. Oecologia. 2020;194(1–2):65–74. pmid:32876762
  55. 55. Gårdmark A, Huss M. Individual variation and interactions explain food web responses to global warming. Philos Trans R Soc Lond B Biol Sci. 2020;375(1814):20190449. pmid:33131431
  56. 56. Graciá E. Biotic interactions matter in phylogeography research: Integrative analysis of demographic, genetic and distribution data to account for them. Mol Ecol. 2020;29(23):4503–5. pmid:33070388
  57. 57. Kissling W, Dormann C, Groeneveld J, Hickler T, Kuhn I, McInerny G, et al. Towards novel approaches to modelling biotic interactions in multispecies assemblages at large spatial extents. Journal of Biogeography. 2012;39(12):2163–78.
  58. 58. Aguirre-Gutierrez J, Kissling W, Biesmeijer J, WallisDeVries M, Reemer M, Carvalheiro L. Historical changes in the importance of climate and land use as determinants of Dutch pollinator distributions. Journal of Biogeography. 2017;44(3):696–707.
  59. 59. Kolanowska M. The future of a montane orchid species and the impact of climate change on the distribution of its pollinators and magnet species. Global Ecology and Conservation. 2021;32.
  60. 60. Fay M. Orchid conservation: how can we meet the challenges in the twenty-first century? Botanical Studies. 2018;59. pmid:29872972
  61. 61. Fay M, Pailler T, Dixon K. Orchid conservation: making the links. Annals of Botany. 2015;116(3):377–9. pmid:26311710
  62. 62. Gale S, Fischer G, Cribb P, Fay M. Orchid conservation: bridging the gap between science and practice. Botanical Journal of the Linnean Society. 2018;186(4):425–34.
  63. 63. Kaye TN, Bahm MA, Thorpe AS, Gray EC, Pfingsten I, Waddell C. Population extinctions driven by climate change, population size, and time since observation may make rare species databases inaccurate. PLoS One. 2019;14(10):e0210378. pmid:31622344
  64. 64. Liu Q, Chen J, Corlett R, Fan X, Yu D, Yang H, et al. Orchid conservation in the biodiversity hotspot of southwestern China. Conservation Biology. 2015;29(6):1563–72. pmid:26372504
  65. 65. McInerny GJ, Roberts DL, Davy AJ, Cribb PJ. Significance of sighting rate in inferring extinction and threat. Conserv Biol. 2006;20(2):562–7. pmid:16903116
  66. 66. Swarts ND, Dixon KW. Terrestrial orchid conservation in the age of extinction. Ann Bot. 2009;104(3):543–56. pmid:19218582
  67. 67. Zhou Z, Shi R, Zhang Y, Xing X, Jin X. Orchid conservation in China from 2000 to 2020: Achievements and perspectives. Plant Divers. 2021;43(5):343–9. pmid:34816060
  68. 68. Wraith J, Pickering C. Quantifying anthropogenic threats to orchids using the IUCN Red List. Ambio. 2018;47(3):307–17. pmid:29043611
  69. 69. Atlas of Living Australia. Atlas of Living Australia occurrence download. 2020.
  70. 70. GBIF. Occurrence dataset for L. excelsa. 2020.
  71. 71. GBIF. Occurrence data for C. leptochila. 2020.
  72. 72. Brown J. SDMtoolbox: a python-based GIS toolkit for landscape genetic, biogeographic and species distribution model analyses. Methods in Ecology and Evolution. 2014;5(7):694–700.
  73. 73. Brown JL, Bennett JR, French CM. SDMtoolbox 2.0: the next generation Python-based GIS toolkit for landscape genetic, biogeographic and species distribution model analyses. PeerJ. 2017;5:e4095. pmid:29230356
  74. 74. Al-Qaddi N, Vessella F, Stephan J, Al-Eisawi D, Schirone B. Current and future suitability areas of kermes oak (Quercus coccifera L.) in the Levant under climate change. Regional Environmental Change. 2017;17(1):143–56.
  75. 75. Bergamin R, Debastiani V, Joner D, Lemes P, Guimaraes T, Loyola R, et al. Loss of suitable climatic areas for Araucaria forests over time. Plant Ecology & Diversity. 2019;12(2):115–26.
  76. 76. Faroni-Perez L. Climate and environmental changes driving idiosyncratic shifts in the distribution of tropical and temperate worm reefs. Journal of the Marine Biological Association of the United Kingdom. 2017;97(5):1023–35.
  77. 77. Hernandez-Lambrano R, Gonzalez-Moreno P, Sanchez-Agudo J. Towards the top: niche expansion of Taraxacum officinale and Ulex europaeus in mountain regions of South America. Austral Ecology. 2017;42(5):577–89.
  78. 78. Yan Y, Wang Y, Feng C, Wan P, Chang K. Potential distributional changes of invasive crop pest species associated with global climate change. Applied Geography. 2017;82:83–92.
  79. 79. Elith J, Phillips S, Hastie T, Dudik M, Chee Y, Yates C. A statistical explanation of MaxEnt for ecologists. Diversity and Distributions. 2011;17(1):43–57.
  80. 80. Phillips S, Anderson R, Schapire R. Maximum entropy modeling of species geographic distributions. Ecological Modelling. 2006;190(3–4):231–59.
  81. 81. Phillips S, Dudik M. Modeling of species distributions with Maxent: new extensions and a comprehensive evaluation. Ecography. 2008;31(2):161–75.
  82. 82. Fick S, Hijmans R. WorldClim 2: new 1-km spatial resolution climate surfaces for global land areas. International Journal of Climatology. 2017;37(12):4302–15.
  83. 83. Barve N, Barve V, Jimenez-Valverde A, Lira-Noriega A, Maher S, Peterson A, et al. The crucial role of the accessible area in ecological niche modeling and species distribution modeling. Ecological Modelling. 2011;222(11):1810–9.
  84. 84. McGee R, Williams S, Poulton R, Moffitt T. A longitudinal study of cannabis use and mental health from adolescence to early adulthood. Addiction. 2000;95(4):491–503. pmid:10829326
  85. 85. Meinshausen M, Nicholls Z, Lewis J, Gidden M, Vogel E, Freund M, et al. The shared socio-economic pathway (SSP) greenhouse gas concentrations and their extensions to 2500. Geoscientific Model Development. 2020;13(8):3571–605.
  86. 86. Li J, Chen X, Kurban A, Van de Voorde T, De Maeyer P, Zhang C. Coupled SSPs-RCPs scenarios to project the future dynamic variations of water-soil-carbon-biodiversity services in Central Asia. Ecological Indicators. 2021;129.
  87. 87. Poggio L, de Sousa LM, Batjes NH, Heuvelink GBM, Kempen B, Ribeiro E, et al. SoilGrids 2.0: producing soil information for the globe with quantified spatial uncertainty. Soil. 2021;7:217–40.
  88. 88. Owens H, Campbell L, Dornak L, Saupe E, Barve N, Soberon J, et al. Constraints on interpretation of ecological niche models by limited environmental ranges on calibration areas. Ecological Modelling. 2013;263:10–8.
  89. 89. Gill N, Sangermano F. Africanized honeybee habitat suitability: a comparison between models for southern Utah and southern California. Applied Geography. 2016;76:14–21.
  90. 90. Veloz S. Spatially autocorrelated sampling falsely inflates measures of accuracy for presence-only niche models. Journal of Biogeography. 2009;36(12):2290–9.
  91. 91. Allouche O, Steinitz O, Rotem D, Rosenfeld A, Kadmon R. Incorporating distance constraints into species distribution models. Journal of Applied Ecology. 2008;45(2):599–609.
  92. 92. Tsoar A, Allouche O, Steinitz O, Rotem D, Kadmon R. A comparative evaluation of presence-only methods for modelling species distribution. Diversity and Distributions. 2007;13(4):397–405.
  93. 93. Kolanowska M, Rewicz A, Baranow P. Ecological niche modeling of the pantropical orchid Polystachya concreta (Orchidaceae) and its response to climate change. Sci Rep. 2020;10(1):14801. pmid:32908206
  94. 94. Herold N, Ekström M, Kala J, Goldie J, Evans JP. Australian climate extremes in the 21st century according to a regional climate model ensemble: Implications for health and agriculture. Weather and Climate Extremes. 2018;20:54–68.
  95. 95. Baguette M, Bertrand JAM, Stevens VM, Schatz B. Why are there so many bee-orchid species? Adaptive radiation by intra-specific competition for mnesic pollinators. Biol Rev Camb Philos Soc. 2020;95(6):1630–63. pmid:32954662
  96. 96. Hutchings MJ, Robbirt KM, Roberts DL, Davy AJ. Vulnerability of a specialized pollination mechanism to climate change revealed by a 356-year analysis. Botanical Journal of the Linnean Society. 2018;186(4):498–509.
  97. 97. Robbirt KM, Roberts DL, Hutchings MJ, Davy AJ. Potential disruption of pollination in a sexually deceptive orchid by climatic change. Curr Biol. 2014;24(23):2845–9. pmid:25454786
  98. 98. Weyland H. Bestäubungsbiologische Beobachtungen an Ophrys sphegodes Miller und Ophrys tommasinii Visiani und einigen anderen Ophrys-Arten in Istrien und Griechenland (Peloponnes). Berichte Aus Den Arbeitskreisen Heimische Orchideen. 2013;30(2):160–84.
  99. 99. Lara Ruiz J. Polinizadores y visitantes de Ophrys L. en la Península Ibérica e Islas Baleares. Micobotánica-Jaén. 2010;5(3):http://www.micobotanicajaen.com/Revista/Articulos/JLaraR/Polinizadores/Ophrys.html.
  100. 100. Paulus HF. Zur Bestäubungsbiologie von Serapias lingua und einiger Ophrys-Arten in Kroatien (Orchidaceae und Insecta, Apoidea). Journal Europäischer Orchideen. 2014;46(3/4):501–55.
  101. 101. Wilcox Y. Confirmation de la pollinisation d’un Ophrys par la mouche Myopa tessellatipennis Motschulsky. L’Orchidophile. 2012;43(3):181–3.
  102. 102. Tsiftsis S, Djordjević V. Modelling sexually deceptive orchid species distributions under future climates: the importance of plant-pollinator interactions. Sci Rep. 2020;10(1):10623. pmid:32606363
  103. 103. Vanden Broeck A, Van Landuyt W, Cox K, De Bruyn L, Gyselings R, Oostermeijer G, et al. High levels of effective long-distance dispersal may blur ecotypic divergence in a rare terrestrial orchid. BMC Ecology. 2014;14(1):20.
  104. 104. Hedren M, Olofsson S, Paun O. Orchid colonization: multiple parallel dispersal events and mosaic genetic structure in Dactylorhiza majalis ssp. lapponica on the Baltic island of Gotland. Annals of Botany. 2018;122(6):1019–31. pmid:29955767
  105. 105. Parrot AW. New Zealand Ichneumonidae II. Transactions and Proceedings of the Royal Society of New Zealand. 1952;80:155–7.
  106. 106. Brunton Martin AL, Gaskett AC, O’Hanlon JC. Museum records indicate male bias in pollinators of sexually deceptive orchids. Naturwissenschaften. 2021;108(4):25. pmid:34091791
  107. 107. Weinstein A, Davis B, Menz M, Dixon K, Phillips R. Behaviour of sexually deceived ichneumonid wasps and its implications for pollination in Cryptostylis (Orchidaceae). Biological Journal of the Linnean Society. 2016;119(2):283–98.
  108. 108. Alyssa M. Weinstein A, Davis BJ, Menz MHM, Dixon KW, Phillips RD. Behaviour of sexually deceived ichneumonid wasps and its implications for pollination in Cryptostylis (Orchidaceae). Biological Journal of the Linnean Society. 2016(119):283–98.
  109. 109. Tomlinson S, Phillips RD. Metabolic rate, evaporative water loss and field activity in response to temperature in an ichneumonid wasp. Journal of Zoology. 2012;287(2):81–90.
  110. 110. Baker GH, Tann CR. Long-term changes and host plant differences in the incidence of parasitoids attacking Helicoverpa spp. (Lepidoptera: Noctuidae) in agricultural landscapes in eastern Australia. Austral Entomology. 2020;59(2):386–97.
  111. 111. Hill GM, Kawahara AY, Daniels JC, Bateman CC, Scheffers BR. Climate change effects on animal ecology: butterflies and moths as a case study. Biological Reviews. 2021;96(5):2113–26. pmid:34056827
  112. 112. Pelini SL, Prior KM, Parker DJ, Dzurisin JDK, Lindroth RL, Hellmann JJ. Chapter 11—Climate Change and Temporal and Spatial Mismatches in Insect Communities. In: Letcher TM, editor. Climate Change. Amsterdam: Elsevier; 2009. p. 215–31.
  113. 113. Van Nouhuys S, Lei G. Parasitoid-Host Metapopulation Dynamics: The Causes and Consequences of Phenological Asynchrony. Journal of Animal Ecology. 2004;73(3):526–35.
  114. 114. Ma X, Kang J, Nontachaiyapoom S, Wen T, Hyde K. Non-mycorrhizal endophytic fungi from orchids. Current Science. 2015;109(1):72–87.
  115. 115. Sarsaiya S, Shi J, Chen J. A comprehensive review on fungal endophytes and its dynamics on Orchidaceae plants: current research, challenges, and future possibilities. Bioengineered. 2019;10(1):316–34. pmid:31347943
  116. 116. Kaur J, Sharma J. Orchid Root Associated Bacteria: Linchpins or Accessories? Front Plant Sci. 2021;12:661966. pmid:34249034
  117. 117. Rasmussen HN, Rasmussen FN. Orchid Mycorrhiza: Implications of a Mycophagous Life Style. Oikos. 2009;118(3):334–45.
  118. 118. Waud M, Brys R, Van Landuyt W, Lievens B, Jacquemyn H. Mycorrhizal specificity does not limit the distribution of an endangered orchid species. Mol Ecol. 2017;26(6):1687–701. pmid:28100022
  119. 119. Li T, Wu S, Yang W, Selosse M-A, Gao J. How Mycorrhizal Associations Influence Orchid Distribution and Population Dynamics. Frontiers in Plant Science. 2021;12. pmid:34025695
  120. 120. Oktalira FT, Whitehead MR, Linde CC. Mycorrhizal specificity in widespread and narrow-range distributed Caladenia orchid species. Fungal Ecology. 2019;42:100869.
  121. 121. McCormick MK, Jacquemyn H. What constrains the distribution of orchid populations? New Phytologist. 2014;202(2):392–400.
  122. 122. Seneviratne SI, Nicholls N, Easterling D, Goodess CM, Kanae S, Kossin J, et al. Changes in Climate Extremes and their Impacts on the Natural Physical Environment. In: Field CB, Barros V, Stocker TF, Qin D, Dokken DJ, Ebi KL, et al., editors. Managing the Risks of Extreme Events and Disasters to Advance Climate Change Adaptation. Cambridge, UK, and New York, NY, USA: Cambridge University Pre; 2012. p. 109–230.
  123. 123. Brown AP, Thomson-Dans C, Marchant N. Western Australia’s Threatened Flora. Perth: Department of Conservation and Land Management; 1998.