Multi‐site generalised dissimilarity modelling: using zeta diversity to differentiate drivers of turnover in rare and widespread species

European river typologies fail to capture diatom, fish, and macrophyte community composition

Typology systems are frequently used in applied and fundamental ecology and are relevant for environmental monitoring and conservation. They aggregate ecosystems into discrete types based on biotic and abiotic variables, assuming that ecosystems of the same type are more alike than ecosystems of different types with regard to a specific property of interest. We evaluated whether this assumption is met by the Broad River Types (BRT), a recently proposed European river typology system, that classifies river segments based on abiotic variables, when it is used to group biological communities. We compiled data on the community composition of diatoms, fishes, and aquatic macrophytes throughout Europe and evaluated whether the composition is more similar in site groups with the same river type than in site groups of different river types using analysis of similarities, classification strength, typical species analysis, and the area under zeta diversity decline curves. We compared the performance of the BRT with those of four region-based typology systems, namely, Illies Freshwater Ecoregions, the Biogeographic Regions, the Freshwater Ecoregions of the World, and the Environmental Zones, as well as spatial autocorrelation (SA) classifications. All typology systems received low scores from most evaluation methods, relative to predefined thresholds and the SA classifications. The BRT often scored lowest of all typology systems. Within each typology system, community composition overlapped considerably between site groups defined by the types of the systems. The overlap tended to be the lowest for fishes and between Illies Freshwater Ecoregions. In conclusion, we found that existing broad-scale river typology systems fail to delineate site groups with distinct and compositionally homogeneous communities of diatoms, fishes, and macrophytes. A way to improve the fit between typology systems and biological communities might be to combine segment-based and region-based typology systems to simultaneously account for local environmental variation and historical distribution patterns, thus potentially improving the utility of broad-scale typology systems for freshwater biota.

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.

Zooplankton community turnover in response to a novel TiO2-coated nano-formulation of carbendazim and its constituents

Novel nanomaterial-based pesticide formulations are increasingly perceived as promising aids in the transition to more efficient agricultural production systems. The current understanding of potential unintended (eco)toxicological impacts of nano-formulated pesticides is scarce, in particular with regard to (non-target) aquatic organisms and ecosystems. The present study reports the results of a long-term freshwater mesocosm experiment which assessed responses of individual zooplankton taxa and communities to a novel TiO2-coated nano-formulation of the fungicide carbendazim. Population- and community trends were assessed and compared in response to the nano-formulation and its constituents applied individually (i.e. nano-sized TiO2, carbendazim) and in combination (i.e. nano-sized TiO2 & carbendazim). Minimal differences were observed between effects induced by the nano-formulation and its active ingredient (i.e. carbendazim) when applied at equivalent nominal test concentrations (4 μg L-1). Nano-sized TiO2 was found to affect zooplankton community trends when applied separately at environmentally realistic concentrations (20 μg L-1 nominal test concentration). However, when nano-sized TiO2 was applied in combination with carbendazim, nano-sized TiO2 was found not to alter effects on community trends induced by carbendazim. The findings of the current study provide an extensive and timely addition to the current body of work available on non-target impacts of nano-formulated pesticides.

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.

Zeta diversity patterns in metabarcoded lotic algal assemblages as a tool for bioassessment

Assessments of the ecological health of algal assemblages in streams typically focus on measures of their local diversity and classify individuals by morphotaxonomy. Such assemblages are often connected through various ecological processes, such as dispersal, and may be more accurately assessed as components of regional-, rather than local-scale assemblages. With recent declines in the costs of sequencing and computation, it has also become increasingly feasible to use metabarcoding to more accurately classify algal species and perform regional-scale bioassessments. Recently, zeta diversity has been explored as a novel method of constructing regional bioassessments for groups of streams. Here, we model the use of zeta diversity to investigate whether stream health can be determined by the landscape diversity of algal assemblages. We also compare the use of DNA metabarcoding and morphotaxonomy classifications in these zeta diversity-based bioassessments of regional stream health. From 96 stream samples in California, we used various orders of zeta diversity to construct models of biotic integrity for multiple assemblages of diatoms, as well as hybrid assemblages of diatoms in combination with soft-bodied algae, using taxonomy data generated with both DNA sequencing as well as traditional morphotaxonomic approaches. We compared our ability to evaluate the ecological health of streams with the performance of multiple algal indices of biological condition. Our zeta diversity-based models of regional biotic integrity were more strongly correlated with existing indices for algal assemblages classified using metabarcoding compared to morphotaxonomy. Metabarcoding for diatoms and hybrid algal assemblages involved rbcL and 18S V9 primers, respectively. Importantly, we also found that these algal assemblages, independent of the classification method, are more likely to be assembled under a process of niche differentiation rather than stochastically. Taken together, these results suggest the potential for zeta diversity patterns of algal assemblages classified using metabarcoding to inform stream bioassessments.

Drivers of compositional turnover of narrow-ranged versus widespread naturalised woody plants in South Africa

Introduction

Alien trees and shrubs have become increasingly common invaders globally and have caused major negative impacts to ecosystems and society. Non-native woody plant species make up the majority of legislated invasive alien taxa in South Africa and contribute substantially to recorded negative impacts. It is of management interest to elucidate the macroecological processes that mediate the assembly of alien taxa, as this is expected to be associated with anthropogenic factors (e.g., human activity, introduction events, pathways of propagule dispersal mediated by humans) and bioclimatic factors (such as diurnal temperature range and precipitation gradients). These analyses require large species-occurrence datasets with comprehensive sampling across broad environmental conditions. Efforts of citizen scientists produce large numbers of occurrence records in a consistent manner which may be utilised for scientific investigations.

Methods

Research Grade occurrence data on naturalised plants of South Africa were extracted from the citizen scientist platform iNaturalist. Sampling bias was mitigated using statistical modelling of background points estimated from a Target Group of species which identifies well sampled communities. The drivers of assembly for alien plants at different range sizes were identified using multi-site generalised dissimilarity modelling (MS-GDM) of zeta diversity. The predicted compositional similarity between all cells was computed based on the subset of identified well sampled communities and using generalised dissimilarity modelling (GDM). From this, alien bioregions were identified using a k-means cluster analysis.

Results and Discussion

Bioclimatic factors significantly influenced community turnover in inland areas with large diurnal temperature ranges, and in areas with high precipitation. Communities separated by large geographical distances had significantly different compositions, indicating little contribution of long-range propagule movement by humans, and the presence of localised introduction hubs within the country which harbour unique species compositions. Analyses also showed a significant contribution of road density to turnover, which may be moderated by the habitat service provided by road verges. The same is true for natural dispersal via rivers in arid areas. The distribution of naturalised tree and shrub species is geographically clustered and forms six alien bioregions that are distinct from the South African biomes defined by native species distributionanalysis.

The species composition of local flea assemblages at a small scale in two South American regions is predominantly driven by niche-based mechanisms

We applied a step-down factor analysis (SDFA) and multi-site generalised dissimilarity modelling (MS-GDM) to local flea communities harboured by small mammals (i.e., collected at small sampling sites over a short time period) in two South American regions (Patagonia and the Northwestern Argentina) with the aim of understanding whether these communities were assembled via niche-based or dispersal-based processes. The SDFA allows us to determine whether clusters of flea assemblages across different types of climates, vegetation and soils can be distinguished (suggesting niche-based assembly). MS-GDM allows us to determine whether a substantial proportion of the variation in flea species turnover is explained by specific climate-associated, vegetation-associated and soil-associated variables (indicating niche-based assembly) or host turnover (indicating dispersal-based assembly). Mapping of assemblages on climate, vegetation and soil maps, according to their loadings on axis 1 or axis 2 of the SDFA, did not provide clear-cut results. Clusters of similar loadings could be recognized within some, but not other, climate, vegetation and soil types. However, MS-GDM demonstrated that the effect of environmental variables (especially air temperature) on flea compositional turnover was much stronger than that of host turnover, indicating the predominance of niche-based processes in local community assembly. A comparison of our results with those on the mechanisms that drive species assembly in regional communities allows us to conclude that local and regional communities result from the joint action of niche-based and dispersal-based processes, with the former more important at a smaller spatial scale and the latter at a larger spatial scale.

Applications of environmental DNA (eDNA) in agricultural systems: Current uses, limitations and future prospects

Global food production, food supply chains and food security are increasingly stressed by human population growth and loss of arable land, becoming more vulnerable to anthropogenic and environmental perturbations. Numerous mutualistic and antagonistic species are interconnected with the cultivation of crops and livestock and these can be challenging to identify on the large scales of food production systems. Accurate identifications to capture this diversity and rapid scalable monitoring are necessary to identify emerging threats (i.e. pests and pathogens), inform on ecosystem health (i.e. soil and pollinator diversity), and provide evidence for new management practices (i.e. fertiliser and pesticide applications). Increasingly, environmental DNA (eDNA) is providing rapid and accurate classifications for specific organisms and entire species assemblages in substrates ranging from soil to air. Here, we aim to discuss how eDNA is being used for monitoring of agricultural ecosystems, what current limitations exist, and how these could be managed to expand applications into the future. In a systematic review we identify that eDNA-based monitoring in food production systems accounts for only 4 % of all eDNA studies. We found that the majority of these eDNA studies target soil and plant substrates (60 %), predominantly to identify microbes and insects (60 %) and are biased towards Europe (42 %). While eDNA-based monitoring studies are uncommon in many of the world's food production systems, the trend is most pronounced in emerging economies often where food security is most at risk. We suggest that the biggest limitations to eDNA for agriculture are false negatives resulting from DNA degradation and assay biases, as well as incomplete databases and the interpretation of abundance data. These require in silico, in vitro, and in vivo approaches to carefully design, test and apply eDNA monitoring for reliable and accurate taxonomic identifications. We explore future opportunities for eDNA research which could further develop this useful tool for food production system monitoring in both emerging and developed economies, hopefully improving monitoring, and ultimately food security.

Spatial turnover of core and occasional bacterial taxa in the plastisphere from a plateau river, China

Plastic surfaces in the environment are a comparatively new niche for microbial colonization, also known as the "plastisphere". However, our understanding of the core and occasional bacterial taxa in the plastisphere is limited. Here, environmental plastic, water, and sediment samples were collected from 10 sites in a plateau river (Lhasa River, China) in September of 2019. The composition and spatial turnover of core and occasional bacterial taxa in the plastisphere were revealed via 16S rRNA gene sequencing and compared with water and sediment. The results indicated that deterministic processes dominated the habitat specialization that shaped the formation of core and occasional taxa in the plastisphere, water, and sediment of the Lhasa River because the decline in zeta diversity in the plastisphere, water, and sediment was more fitted to a power-law form rather than an exponential form. Proteobacteria (65.9%), Bacteroidetes (16.0%), and Cyanobacteria (11.7%) dominated the plastic core taxa. Core taxa rather than occasional taxa in the plastisphere had a lower (21.7%) proportion of OTUs and a higher (81.7%) proportion of average relative abundance than water and sediment, which were dominant in plastic bacterial communities. The spatial turnover of core and occasional bacterial taxa in the plastisphere was governed by abiotic as well as biotic factors. Specifically, the spatial turnover of core taxa in the plastisphere with high connectivity but low functional redundancy was easily affected by geographical distance, altitude, and heavy metals. Furthermore, strong drug resistance was found in the spatially persistent core taxa in the plastisphere. This study provides empirical support for the spatial turnover (species variation) and potential ecological mechanisms of bacterial communities in the plastisphere from river ecosystems.

Impact of marine protected areas on temporal stability of fish species diversity

Preserving biodiversity over time is a pressing challenge for conservation science. A key goal of marine protected areas (MPAs) is to maintain stability in species composition, via reduced turnover, to support ecosystem function. Yet, this stability is rarely measured directly under different levels of protection. Rather, evaluations of MPA efficacy generally consist of static measures of abundance, species richness, and biomass, and rare measures of turnover are limited to short-term studies involving pairwise (beta diversity) comparisons. Zeta diversity is a recently developed metric of turnover that allows for measurement of compositional similarity across multiple assemblages and thus provides more comprehensive estimates of turnover. We evaluated the effectiveness of MPAs at preserving fish zeta diversity across a network of marine reserves over 10 years in Batemans Marine Park, Australia. Snorkel transect surveys were conducted across multiple replicated and spatially interspersed sites to record fish species occurrence through time. Protection provided by MPAs conferred greater stability in fish species turnover. Marine protected areas had significantly shallower decline in zeta diversity compared with partially protected and unprotected areas. The retention of harvested species was four to six times greater in MPAs compared with partially protected and unprotected areas, and the stabilizing effects of protection were observable within 4 years of park implementation. Conversely, partial protection offered little to no improvement in stability, compared with unprotected areas. These findings support the efficacy of MPAs for preserving temporal fish diversity stability. The implementation of MPAs helps stabilize fish diversity and may, therefore, support biodiversity resilience under ongoing environmental change.

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.

A multi-site method to capture turnover in rare to common interactions in bipartite species networks

Ecological network structure is maintained by a generalist core of common species. However, rare species contribute substantially to both the species and functional diversity of networks. Capturing changes in species composition and interactions, measured as turnover, is central to understanding the contribution of rare and common species and their interactions. Due to a large contribution of rare interactions, the pairwise metrics used to quantify interaction turnover are, however, sensitive to compositional change in the interactions of, often rare, peripheral specialists rather than common generalists in the network. Here we expand on pairwise interaction turnover using a multi-site metric that enables quantifying turnover in rare to common interactions (in terms of occurrence of interactions). The metric further separates this turnover into interaction turnover due to species turnover and interaction rewiring. We demonstrate the application and value of this method using a host-parasitoid system sampled along gradients of environmental modification. In the study system, both the type and amount of habitat needed to maintain interaction composition depended on the properties of the interactions considered, that is, from rare to common. The analyses further revealed the potential of host switching to prevent or delay species loss, and thereby buffer the system from perturbation. Multi-site interaction turnover provides a comprehensive measure of network change that can, for example, detect ecological thresholds to habitat loss for rare to common interactions. Accurate description of turnover in common, in addition to rare, species and their interactions is particularly relevant for understanding how network structure and function can be maintained.

Rarity facets of biodiversity: Integrating Zeta diversity and Dark diversity to understand the nature of commonness and rarity

Measuring commonness and rarity is pivotal to ecology and conservation. Zeta diversity, the average number of species shared by multiple sets of assemblages, and Dark diversity, the number of species that could occur in an assemblage but are missing, have been recently proposed to capture two aspects of the commonness-rarity spectrum. Despite a shared focus on commonness and rarity, thus far, Zeta and Dark diversities have been assessed separately. Here, we review these two frameworks and suggest their integration into a unified paradigm of the "rarity facets of biodiversity." This can be achieved by partitioning Alpha and Beta diversities into five components (the Zeta, Eta, Theta, Iota, and Kappa rarity facets) defined based on the commonness and rarity of species. Each facet is assessed in traditional and multiassemblage fashions to bridge conceptual differences between Dark diversity and Zeta diversity. We discuss applications of the rarity facets including comparing the taxonomic, functional, and phylogenetic diversity of rare and common species, or measuring species' prevalence in different facets as a metric of species rarity. The rarity facets integrate two emergent paradigms in biodiversity science to better understand the ecology of commonness and rarity, an important endeavor in a time of widespread changes in biodiversity across the Earth.

Metabolic flexibility allows bacterial habitat generalists to become dominant in a frequently disturbed ecosystem

Ecological theory suggests that habitat disturbance differentially influences distributions of habitat generalist and specialist species. While well-established for macroorganisms, this theory has rarely been explored for microorganisms. Here we tested these principles in permeable (sandy) sediments, ecosystems with much spatiotemporal variation in resource availability and physicochemical conditions. Microbial community composition and function were profiled in intertidal and subtidal sediments using 16S rRNA gene amplicon sequencing and metagenomics, yielding 135 metagenome-assembled genomes. Community composition and metabolic traits modestly varied with sediment depth and sampling date. Several taxa were highly abundant and prevalent in all samples, including within the orders Woeseiales and Flavobacteriales, and classified as habitat generalists; genome reconstructions indicate these taxa are highly metabolically flexible facultative anaerobes and adapt to resource variability by using different electron donors and acceptors. In contrast, obligately anaerobic taxa such as sulfate reducers and candidate lineage MBNT15 were less abundant overall and only thrived in more stable deeper sediments. We substantiated these findings by measuring three metabolic processes in these sediments; whereas the habitat generalist-associated processes of sulfide oxidation and fermentation occurred rapidly at all depths, the specialist-associated process of sulfate reduction was restricted to deeper sediments. A manipulative experiment also confirmed habitat generalists outcompete specialist taxa during simulated habitat disturbance. Together, these findings show metabolically flexible habitat generalists become dominant in highly dynamic environments, whereas metabolically constrained specialists are restricted to narrower niches. Thus, an ecological theory describing distribution patterns for macroorganisms likely extends to microorganisms. Such findings have broad ecological and biogeochemical ramifications.

Land use, but not distance, drives fungal beta diversity

Fungi are one of the most diverse taxonomic groups on the planet, but much of their diversity and community organization remains unknown, especially at local scales. Indeed, a consensus on how fungal communities change across spatial or temporal gradients-beta diversity-remains nascent. Here, we use a data set of plant-associated fungal communities (leaf, root, and soil) across multiple land uses from a New Zealand-wide study to look at fungal community turnover at small spatial scales (<1 km). Using hierarchical Bayesian beta regressions and Hill-number-based diversity profiles, we show that fungal communities are often markedly dissimilar at even small distances, regardless of land use. Moreover, diversity profile plots indicate that leaf, root, and soil-associated communities show different patterns in the dominance or rarity of dissimilar species. Leaf-associated communities differed from site to site in their low-abundance species, whereas root-associated communities differed between sites in the dominant species; soil-associated communities were intermediate. Land-use differences were largely driven by the lower turnover between high-productivity grassland sites. Further, we discuss the implications and benefits of using diversity profile plots of turnover to draw inferences into the mechanisms of how communities are structured across spatial gradients.

Measuring change in biological communities: multivariate analysis approaches for temporal datasets with low sample size

Effective and robust ways to describe, quantify, analyse, and test for change in the structure of biological communities over time are essential if ecological research is to contribute substantively towards understanding and managing responses to ongoing environmental changes. Structural changes reflect population dynamics, changes in biomass and relative abundances of taxa, and colonisation and extinction events observed in samples collected through time. Most previous studies of temporal changes in the multivariate datasets that characterise biological communities are based on short time series that are not amenable to data-hungry methods such as multivariate generalised linear models. Here, we present a roadmap for the analysis of temporal change in short-time-series, multivariate, ecological datasets. We discuss appropriate methods and important considerations for using them such as sample size, assumptions, and statistical power. We illustrate these methods with four case-studies analysed using the R data analysis environment.

Compositional turnover in ecto- and endoparasite assemblages of an African bat, Miniopterus natalensis (Chiroptera, Miniopteridae): effects of hierarchical scale and host sex

We studied the compositional turnover in infracommunities and component communities of ecto- and endoparasites infesting a bat, Miniopterus natalensis (Chiroptera, Miniopteridae), across seven sampling sites using the zeta diversity metric (measuring similarity between multiple communities) and calculating zeta decline and retention rate (both scales) and zeta decay (component communities). We asked whether the patterns of zeta diversity differ between (a) infracommunities and component communities; (b) ecto- and endoparasites and (c) subsets of communities infecting male and female bats. The pattern of compositional turnover differed between infracommunities and component communities in endoparasites only. The shape of zeta decline for infracommunities indicated that there were approximately equal probabilities of ecto- and endoparasitic species to occur on/in any bat individual within a site. The shape of zeta decline for component communities suggested the stochasticity of ectoparasite turnover, whereas the turnover of endoparasites was driven by niche-based processes. Compositional turnover in component communities of ectoparasites was more spatially dependent than that of endoparasites. Spatial independence of compositional turnover in endoparasites was due to subcommunities harboured by female bats. We conclude that the patterns of compositional turnover in infracommunities were similar in ecto- and endoparasites, whereas the patterns of turnover in component communities differed between these groups.

Soil Bacterial Communities Exhibit Strong Biogeographic Patterns at Fine Taxonomic Resolution

Bacteria have been inferred to exhibit relatively weak biogeographic patterns. To what extent such findings reflect true biological phenomena or methodological artifacts remains unclear. Here, we addressed this question by analyzing the turnover of soil bacterial communities from three data sets. We applied three methodological innovations: (i) design of a hierarchical sampling scheme to disentangle environmental from spatial factors driving turnover; (ii) resolution of 16S rRNA gene amplicon sequence variants to enable higher-resolution community profiling; and (iii) application of the new metric zeta diversity to analyze multisite turnover and drivers. At fine taxonomic resolution, rapid compositional turnover was observed across multiple spatial scales. Turnover was overwhelmingly driven by deterministic processes and influenced by the rare biosphere. The communities also exhibited strong distance decay patterns and taxon-area relationships, with z values within the interquartile range reported for macroorganisms. These biogeographical patterns were weakened upon applying two standard approaches to process community sequencing data: clustering sequences at 97% identity threshold and/or filtering the rare biosphere (sequences lower than 0.05% relative abundance). Comparable findings were made across local, regional, and global data sets and when using shotgun metagenomic markers. Altogether, these findings suggest that bacteria exhibit strong biogeographic patterns, but these signals can be obscured by methodological limitations. We advocate various innovations, including using zeta diversity, to advance the study of microbial biogeography.IMPORTANCE It is commonly thought that bacterial distributions show lower spatial variation than for multicellular organisms. In this article, we present evidence that these inferences are artifacts caused by methodological limitations. Through leveraging innovations in sampling design, sequence processing, and diversity analysis, we provide multifaceted evidence that bacterial communities in fact exhibit strong distribution patterns. This is driven by selection due to factors such as local soil characteristics. Altogether, these findings suggest that the processes underpinning diversity patterns are more unified across all domains of life than previously thought, which has broad implications for the understanding and management of soil biodiversity.

Spatial and temporal turnover of parasite species and parasite-host interactions: a case study with fleas and gamasid mites parasitic on small mammals

We studied patterns of ectoparasite species turnover and pairwise ectoparasite-host interactions across space and time in fleas and mites harboured by small mammals using a novel metric, zeta diversity (similarity between multiple communities). We asked whether the zeta diversity of parasites and their interactions with hosts follow a similar spatial or temporal trend. We found substantial differences in some (zeta decline and retention rate) but not in other (zeta decay) spatial patterns of zeta diversity between species and interactions, whereas the differences between the patterns of the temporal species versus interaction zeta diversity occurred to a much lesser extent. In particular, the parametric form of zeta decline suggested that the distribution of ectoparasite species across localities is driven mainly by niche-based processes, whereas the spatial distribution of flea-host and mite-host interactions is predominantly stochastic. We also found much stronger variation in the number of shared species and interactions over space than over time. Parasite community composition, in terms of species, appeared to be much more temporally stable than that in terms of parasite-host interactions. The parametric form of temporal zeta decline indicated that both parasite communities and parasite-host networks are assembled over time via niche-based processes.

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.

Drivers of compositional turnover are related to species' commonness in flea assemblages from four biogeographic realms: zeta diversity and multi-site generalised dissimilarity modelling

We investigated drivers of species turnover in fleas parasitic on small mammals in four biogeographic realms using novel methodology (zeta diversity, and Multi-Site Generalised Dissimilarity Modelling). We asked whether (i) flea turnover was better explained by host turnover or environmental variables; (ii) different factors drive the turnover of rare and widespread fleas; (iii) the factors affecting the turnover of rare or widespread fleas differ between realms; and (iv) environmental variables drive flea turnover directly or via their effects on hosts. Dissimilarity in host species composition was the most important factor affecting flea turnover in all realms. In the Afrotropics, the Nearctic, and the Neotropics, this was true mainly for rare species, whereas the zeta diversity of the Palearctic hosts exerted a strong effect on the turnover of both rare and widespread fleas. Dissimilarity in temperature contributed the most to the turnover of rare fleas in the Neotropics and the Palearctic, whereas the turnover of widespread species in these realms was strongly affected by dissimilarity in precipitation. In the Nearctic, dissimilarity in precipitation or temperature mostly affected the turnover of rare fleas or common species, respectively. In the Afrotropics, dissimilarity in the Normalised Difference Vegetation Index and temperature affected the turnover of all species, independently of their level of commonness, while dissimilarity in rainfall was important for the turnover of rare fleas. The responses of flea assemblages to environmental factors represented a combination of direct responses and responses mediated via effects on host turnover. We conclude that host turnover is a more important factor than environmental dissimilarity in its effect on flea species turnover. However, the relative effects of host composition and environment, as well as those of temperature, precipitation and the amount of vegetation, on flea turnover differ (i) between realms and (ii) between rare and common fleas.

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.

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

Background

Exploring species richness and turnover patterns and their drivers can provide new insights into underlying mechanisms shaping community assembly, with significant implications for biodiversity conservation. Here, we explored diversity patterns of non-endemic, neo-endemic and palaeo-endemic vascular plants in Crete, Greece, a Mediterranean hotspot of plant richness and endemism. We evaluated the relationship between α-diversity and environmental (bioclimatic variables, topography), and anthropogenic variables by Generalized Additive Models, after accounting for spatial autocorrelation. Then, we quantified turnover using the novel concept of zeta diversity (the number of shared species by multiple sites), a framework which allows to explore the full spectrum of compositional turnover, the contribution of rare and widespread species to observed patterns and the underlying processes shaping them. Finally, we explored the abiotic and biotic effects, i.e. how well one category of species (non-endemics, palaeo-endemics, neo-endemics) predicts the patterns of the other categories, on zeta diversity by multi-site Generalized Dissimilarity Modelling.

Results

We found a strong correlation between neo-endemic and palaeo-endemic α-diversity, with climate, topography, and human impact driving species richness. Zeta diversity analysis revealed a sharper decrease of shared palaeo-endemic species, followed by neo-endemics, and then by non-endemics with the number of sites considered to estimate compositional turnover. Perhaps, the narrow distributions of palaeo-endemics as relict species and often habitat specialists, thus persisting locally, and of neo-endemics that may have not reached yet their potential geographical range, resulted in the observed zeta diversity decline pattern. Deterministic processes controlled species turnover of rare non-endemic and neo-endemic species, while deterministic and stochastic processes contributed similarly to palaeo-endemic turnover. However, stochasticity dominates in the case of widespread species in all occasions. The environmental and anthropogenic variables were poor predictors of compositional turnover, especially of widespread species. However, the non-endemic species composition was correlated to rare palaeo-endemics and neo-endemics, highlighting the importance of biotic effects in driving turnover patterns.

Conclusions

It seems that centers of neo-endemism of vascular plants coincide with centers of palaeo-endemism in Crete, but species richness and species turnover are shaped by different drivers.

Measuring continuous compositional change using decline and decay in zeta diversity

Incidence, or compositional, matrices are generated for a broad range of research applications in biology. Zeta diversity provides a common currency and conceptual framework that links incidence-based metrics with multiple patterns of interest in biology, ecology, and biodiversity science. It quantifies the variation in species (or OTU) composition of multiple assemblages (or cases) in space or time, to capture the contribution of the full suite of narrow, intermediate, and wide-ranging species to biotic heterogeneity. Here we provide a conceptual framework for the application and interpretation of patterns of continuous change in compositional diversity using zeta diversity. This includes consideration of the survey design context, and the multiple ways in which zeta diversity decline and decay can be used to examine and test turnover in the identity of elements across space and time. We introduce the zeta ratio-based retention rate curve to quantify rates of compositional change. We illustrate these applications using 11 empirical data sets from a broad range of taxa, scales, and levels of biological organization-from DNA molecules and microbes to communities and interaction networks-including one of the original data sets used to express compositional change and distance decay in ecology. We show (1) how different sample selection schemes used during the calculation of compositional change are appropriate for different data types and questions, (2) how higher orders of zeta may in some cases better detect shifts and transitions, and (3) the relative roles of rare vs. common species in driving patterns of compositional change. By exploring the application of zeta diversity decline and decay, including the retention rate, across this broad range of contexts, we demonstrate its application for understanding continuous turnover in biological systems.

Using alpha, beta, and zeta diversity in describing the health of stream-based benthic macroinvertebrate communities

Ecological monitoring of streams has frequently focused on measures describing the taxonomic, and sometimes functional, α diversity of benthic macroinvertebrates (BMIs) within a single sampled community. However, as many ecological processes effectively link BMI stream communities there is a need to describe groups of communities using measures of regional diversity. Here we demonstrate a role for incorporating both a traditional pairwise measure of community turnover, β diversity, in assessing community health as well as ζ diversity, a more generalized framework for describing similarity between multiple communities. Using 4,395 samples of BMI stream communities in California, we constructed a model using measures of α, β, and ζ diversity, which accounted for 71.7% of among-watershed variation in the mean health of communities, as described by the California Streams Condition Index (CSCI). We also investigated the use of ζ diversity in assessing models of stochastic vs. niche assembly across communities of BMIs within watersheds, with the niche assembly model found to be the likelier of the two.