Patterns and drivers of species richness and turnover of neo-endemic and palaeo-endemic vascular plants in a Mediterranean hotspot: the case of Crete, Greece

Tree communities and functional traits determine herbivore compositional turnover

There are many factors known to drive species turnover, although the mechanisms by which these operate are less clear. Based on comprehensive datasets from the largest tree diversity experiment worldwide (BEF-China), we used shared herbivore species (zeta diversity) and multi-site generalized dissimilarity modelling to investigate the patterns and determinants of species turnover of Lepidoptera herbivores among study plots across a gradient in tree species richness. We found that zeta diversity declined sharply with an increasing number of study plots, with complete changes in caterpillar species composition observed even at the fine spatial scale of our study. Plant community characteristics rather than abiotic factors were found to play key roles in driving caterpillar compositional turnover, although these effects varied with an increasing number of study plots considered, due to the varying contributions of rare and common species to compositional turnover. Our study reveals details of the impact of phylogeny- and trait-mediated processes of trees on herbivore compositional turnover, which has implications for forest management and conservation and shows potential avenues for maintenance of heterogeneity in herbivore communities.

Uncovering the determinants of biodiversity hotspots in China: Evidence from the drivers of multiple diversity metrics on insect assemblages and implications for conservation

Understanding large-scale biodiversity patterns and underlying mechanisms during the formation process is essential for guiding conservation efforts. However, previous studies on the identification and formation mechanism of diversity hotspots in China were often limited to a single (alpha) diversity metric, while multiple (beta or zeta) diversity has rarely been used for exploring drivers and conservation actions. Here, a comprehensive species distribution dataset consisting of representative families of three insect orders was compiled to explore biodiversity hotspots based on different algorithms. Furthermore, to assess the effects of environmental factors on hotspots, we fitted generalized additive mixed-effects models (GAMMs) for species richness, generalized dissimilarity models (GDMs) and multi-site generalized dissimilarity modeling (MS-GDM) for the total beta and zeta diversity. Our results showed that biodiversity hotspots were mainly concentrated in central and southern China, especially in mountainous areas with complex topography, which indicated the insects' affinity to montane environments. Further analyses based on multiple models showed that water-energy factors exerted the strongest explanatory power for the insect assemblage diversity in hotspots of both alpha and beta (or zeta) levels. Additionally, anthropogenic factors also exerted a significant effect on hotspots, and this effect was higher for beta diversity than for alpha diversity. Overall, our study elucidates a comprehensive analysis of the identification and underlying mechanism of biodiversity hotspots in China. Despite several limitations, we still believe that our findings can provide some new insights for conservation efforts in Chinese hotspots.

Environmental drivers of sympatric mammalian species compositional turnover in giant panda nature reserves: Implications for conservation

The charismatic giant panda (Ailuropoda melanoleuca) is an iconic species of wildlife conservation worldwide. As the most effective measure to protect giant pandas and their habitats, China has established 67 giant panda nature reserves (GPNR) during the last five decades, which also bring benefits to many sympatric medium- and large-bodied mammals (MLM). To better inform the planning of the GPNR network with the view of preserving regional MLM diversity, we investigated the zeta diversity (a novel index to measure species compositional turnover considering the contributions of both rare and common species) patterns (i.e. zeta decline and retention rate curve) of MLMs across 40 GPNRs. The effects of species' body mass and conservation status on the zeta diversity patterns were tested. Further, we applied the multi-site generalized dissimilarity modelling (MS-GDM) framework to explore the impacts of environmental and geographic distances on MLM turnover. The results indicated that there are a core set of 17 MLM species sympatric with the giant panda in the GPNRs. Species' body mass can affect the patterns of zeta decline and retention rate curves, and the number of large-bodied species shared by multiple GPNRs is higher than that of medium-bodied species across zeta orders. The MS-GDM revealed the important roles of difference in habitat heterogeneity and spatial distance between GPNRs in driving MLM turnover. Consequently, we advocate maintaining and increasing the diversity of (natural) habitats in GPNRs to protect giant panda's sympatric MLM diversity. The government should consider optimizing the GPNR network (e.g. incorporating multiple small GPNRs into one single large reserve) to capture the most turnover of MLMs, and the newly-established Giant Panda National Park is relevant to fulfilling this long-term goal.

Past connections with the mainland structure patterns of insular species richness in a continental-shelf archipelago (Aegean Sea, Greece)

Recent research in island biogeography has highlighted the important role of late Quaternary sea-level fluctuations in shaping biogeographic patterns in insular systems but focused on oceanic systems. Through this study, we aim investigate how late Quaternary sea-level fluctuations shaped species richness patterns in continental-shelf island systems. Focusing on the Aegean archipelago, we first compiled maps of the area's geography using published data, under three sea-level stands: (a) current; (b) median sea-level over the last nine glacial-interglacial cycles (MSL); and (c) Last Glacial Maximum (LGM). We gathered taxon-island occurrences for multiple chorotypes of angiosperms, butterflies, centipedes, and reptiles. We investigated the impact of present-day and past geographic settings on chorological groups by analyzing island species-area relationships (ISARs) and using generalized linear mixed models (GLMMs) selection based on multiple metrics of goodness of fit. Our results confirm that the Aegean's geography has changed dramatically since the LGM, whereas the MSL only modestly differs from the present configuration. Apart for centipedes, paleogeographic changes affected both native and endemic species diversity through altering connections between land-bridge islands and the mainland. On land-bridge islands, we detected over-representation of native species and under-representation of endemics. Unlike oceanic islands, sea-level-driven increase of isolation and area contraction did not strongly shape patterns of species richness. Furthermore, the LGM configurations rather than the MSL configuration shaped patterns of endemic species richness. This suggests that even short episodes of increased connectivity with continental populations are sufficient to counteract the genetic differentiation of insular populations. On the other hand, the over-representation of native nonendemic species on land-bridge islands reflected MSL rather than LGM mainland connections. Our study shows that in terms of processes affecting species richness patterns, continental archipelagos differ fundamentally from oceanic systems because episodic connections with the mainland have profound effects on the biota of land-bridge islands.

Plant Endemism Centres and Biodiversity Hotspots in Greece

Biodiversity hotspots (BH) cover a small fraction of the Earth's surface, yet host numerous endemics. Human-induced biodiversity loss has been increasing worldwide, despite attempts to halt the extinction crisis. There is thus an urgent need to efficiently allocate the available conservation funds in an optimised conservation prioritization scheme. Identifying BH and endemism centres (EC) is therefore a valuable tool in conservation prioritization and planning. Even though Greece is one of the most plant species-rich European countries, few studies have dealt with the identification of BH or EC and none has ever incorporated phylogenetic information or extended to the national scale. Consequently, we are unaware of the extent that Special Areas of Conservation (SAC) of the Natura 2000 network efficiently protect Greek plant diversity. Here, we located for the first time at a national scale and in a phylogenetic framework, the areas serving as BH and EC, and assessed the effectiveness of the Greek SAC in safeguarding them. BH and EC are mainly located near mountainous areas, and in areas supposedly floristically impoverished, such as the central Aegean islands. A critical re-assessment of the Greek SAC might be needed to minimize the extinction risk of the Greek endemics, by focusing the conservation efforts also on the BH and EC that fall outside the established Greek SAC.

Spatial Phylogenetics, Biogeographical Patterns and Conservation Implications of the Endemic Flora of Crete (Aegean, Greece) under Climate Change Scenarios

Human-induced biodiversity loss has been accelerating since the industrial revolution. The climate change impacts will severely alter the biodiversity and biogeographical patterns at all scales, leading to biotic homogenization. Due to underfunding, a climate smart, conservation-prioritization scheme is needed to optimize species protection. Spatial phylogenetics enable the identification of endemism centers and provide valuable insights regarding the eco-evolutionary and conservation value, as well as the biogeographical origin of a given area. Many studies exist regarding the conservation prioritization of mainland areas, yet none has assessed how climate change might alter the biodiversity and biogeographical patterns of an island biodiversity hotspot. Thus, we conducted a phylogenetically informed, conservation prioritization study dealing with the effects of climate change on Crete’s plant diversity and biogeographical patterns. Using several macroecological analyses, we identified the current and future endemism centers and assessed the impact of climate change on the biogeographical patterns in Crete. The highlands of Cretan mountains have served as both diversity cradles and museums, due to their stable climate and high topographical heterogeneity, providing important ecosystem services. Historical processes seem to have driven diversification and endemic species distribution in Crete. Due to the changing climate and the subsequent biotic homogenization, Crete’s unique bioregionalization, which strongly reminiscent the spatial configuration of the Pliocene/Pleistocene Cretan paleo-islands, will drastically change. The emergence of the ‘Anthropocene’ era calls for the prioritization of biodiversity-rich areas, serving as mixed-endemism centers, with high overlaps among protected areas and climatic refugia.

Plant Diversity Patterns and Conservation Implications under Climate-Change Scenarios in the Mediterranean: The Case of Crete (Aegean, Greece)

Climate change poses a great challenge for biodiversity conservation. Several studies exist regarding climate change’s impacts on European plants, yet none has investigated how climate change will affect the extinction risk of the entire endemic flora of an island biodiversity hotspot, with intense human disturbance. Our aim is to assess climate change’s impacts on the biodiversity patterns of the endemic plants of Crete (S Aegean) and provide a case-study upon which a climate-smart conservation planning strategy might be set. We employed a variety of macroecological analyses and estimated the current and future biodiversity, conservation and extinction hotspots in Crete. We evaluated the effectiveness of climatic refugia and the Natura 2000 network of protected areas (PAs) for protecting the most vulnerable species and identified the taxa of conservation priority based on the Evolutionary Distinct and Globally Endangered (EDGE) index. The results revealed that high altitude areas of Cretan mountains constitute biodiversity hotspots and areas of high conservation and evolutionary value. Due to the “escalator to extinction” phenomenon, these areas are projected to become diversity “death-zones” and should thus be prioritised. Conservation efforts should be targeted at areas with overlaps among PAs and climatic refugia, characterised by high diversity and EDGE scores. This conservation-prioritisation planning will allow the preservation of evolutionary heritage, trait diversity and future ecosystem services for human well-being and acts as a pilot for similar regions worldwide.

Patterns of zeta diversity in ectoparasite communities harboured by small mammals at three hierarchical scales: taxon-invariance and scale-dependence

We studied compositional turnover in communities of fleas and mites harboured by small mammals using zeta diversity metric (similarity between multiple communities) and asked whether the patterns of zeta diversity decline with an increase in the number of communities differ between taxa and hierarchical scales [infracommunities (parasite assemblages on individual hosts), component communities (parasite assemblages harboured by host populations), and compound communities (all parasite species in a locality)]. The average number of shared species declined with an increasing number of communities (zeta order). It attained zero at higher orders in infracommunities of both taxa with the shape of the zeta decline being best fitted by the negative exponential function, and the retention rate curves being modal. In contrast, zeta diversity values for compound communities of mites and fleas did not attain zero at higher zeta orders, and the form of the zeta decline was best fitted by the power-law function, whereas the retention rate curves were asymptotic. In component communities, the form of zeta decline was best fitted by either exponential or power-law function in dependence of whether communities were considered within a host across localities or across hosts within a locality and whether ubiquitous species were taken into account. Our main conclusions are that (a) the rules governing compositional turnover in parasite communities for the lowest and the highest hierarchical scales are taxon-invariant but scale-dependent and (b) species composition of infracommunities is mainly driven by stochastic assembly processed, whereas that of compound communities is mainly driven by niche-based processes.

Habitat Islands on the Aegean Islands (Greece): Elevational Gradient of Chasmophytic Diversity, Endemism, Phytogeographical Patterns and need for Monitoring and Conservation

The Aegean archipelago, characterized as a natural laboratory for research concerning plant species diversity and phytogeography has a complex geological and paleogeographical history that varies among its phytogeographical areas. A different combination of factors of variable intensity and duration time drives patterns of its impressive plant species richness and endemism. Cliffs, a conspicuous feature of the Aegean landscape, consist of biologically closed communities that serve as refugia for obligate chasmophytes, the majority of which are Greek or Aegean endemics, and for this reason, they are also considered as habitat islands on the Aegean islands. A synoptic analysis is presented concerning chasmophytic plant diversity focusing on endemic obligate chasmophytes. Phytogeographical patterns of obligate chasmophytes, and especially the endemic ones as well as their elevational range and distribution and zeta diversity, are analyzed and discussed in the frame of climatic change, mentioning that the most threatened endemic obligate chasmophytes are those specialized in high elevation areas, and focusing on the need for monitoring and conservation.