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.

Stronger Geographic Limitations Shape a Rapid Turnover and Potentially Highly Connected Network of Core Bacteria on Microplastics.

Core microbiota is shared microbial taxa within the same habitat, which is important for understanding the stable and consistent components of the complex microbial assembly. However, information on the microplastic core bacteria from the river ecosystems is poor. Here, we investigated the composition and function of microplastic core bacteria from the Three Gorges Reservoir area along the approximate 662 km of the Yangtze River via full-length 16S rRNA gene sequencing, compared with those in water, sediment, and soil. The results showed that the spatial turnover of bacterial communities in four habitats supported deterministic processes dominated by niche differentiation, which shaped their core bacteria. The composition and function of microplastic core bacteria were significantly different from those in the other three habitats. Rhodobacteraceae was the main component of microplastic core bacteria, while the main component of core bacteria in water, sediment, and soil were Burkholderiaceae (21.90%), Burkholderiaceae (5.01%), Nitrosomonadaceae (4.61%), respectively. Furthermore, microplastic core bacteria had stronger geographic limitations along the Yangtze River in the Three Gorges Reservoir area. Stronger geographic limitations shaped the rapid community turnover and a potentially more connected network for the microplastic core bacteria than water, sediment, and soil. More importantly, microplastic core bacteria had strong potential functions of drug resistance and could cause risks to ecosystems and human health. Microplastic core bacteria were mainly influenced by sediment core bacteria, although the bacteria colonizing on microplastics could be from all the contact environments and original sources. These findings provide important insights into the composition, function, and association of microplastic core bacteria with their surrounding environment.

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.

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.

Impactos de las Áreas Protegidas Marinas sobre la Estabilidad Temporal de la Diversidad de Especies de Peces Resumen A medida que avanza el tiempo, la conservación de la biodiversidad es un reto apremiante para las ciencias de la conservación. Un objetivo importante de las áreas marinas protegidas (AMP) es mantener la estabilidad de la composición de especies, por medio de rotaciones reducidas, para así ayudar a la función del ecosistema. Sin embargo, esta estabilidad casi no se mide directamente bajo diferentes niveles de protección. En su lugar, las evaluaciones de eficiencia de las AMP generalmente consisten en medidas estáticas de abundancia, riqueza de especies y biomasa, y las pocas medidas de la rotación están limitadas a los estudios a corto plazo que involucran comparaciones por pares (diversidad beta). La diversidad zeta es una medida recientemente desarrollada de la rotación, la cual permite la medición de las similitudes en la composición en múltiples ensamblajes, proporcionando así estimaciones más completas de la rotación. Evaluamos la efectividad que tienen las AMP en la conservación de la diversidad zeta de los peces en una red de reservas marinas durante diez años en el Parque Marino Bateman, Australia. Se realizaron censos en transecto con snorkel en varios sitios replicados e intercalados espacialmente para registrar la presencia de especies de peces a lo largo del tiempo. La protección proporcionada por las AMP otorgó una mayor estabilidad en la rotación de especies de peces. Las áreas marinas protegidas tuvieron una declinación significativamente más baja de la diversidad zeta que las áreas parcialmente protegidas o desprotegidas. La retención de especies pescadas fue 4-6 veces mayor en las AMP que en las áreas desprotegidas o parcialmente protegidas, y los efectos estabilizadores de la protección fueron observables a partir de cuatro años de la implementación del parque. De manera opuesta, la protección parcial ofreció poca o ninguna estabilidad, comparada con las áreas desprotegidas. Estos descubrimientos respaldan la eficiencia que tienen las AMP en la conservación de la estabilidad temporal de la diversidad de especies de peces. La implementación de las AMP ayuda a estabilizar la diversidad de peces y por lo tanto puede fomentar la resiliencia de la biodiversidad frente al cambio ambiental en curso.

Landscape analyses using eDNA metabarcoding and Earth observation predict community biodiversity in California.

Ecosystems globally are under threat from ongoing anthropogenic environmental change. Effective conservation management requires more thorough biodiversity surveys that can reveal system-level patterns and that can be applied rapidly across space and time. Using modern ecological models and community science, we integrate environmental DNA and Earth observations to produce a time snapshot of regional biodiversity patterns and provide multi-scalar community-level characterization. We collected 278 samples in spring 2017 from coastal, shrub, and lowland forest sites in California, a complex ecosystem and biodiversity hotspot. We recovered 16,118 taxonomic entries from eDNA analyses and compiled associated traditional observations and environmental data to assess how well they predicted alpha, beta, and zeta diversity. We found that local habitat classification was diagnostic of community composition and distinct communities and organisms in different kingdoms are predicted by different environmental variables. Nonetheless, gradient forest models of 915 families recovered by eDNA analysis and using BIOCLIM variables, Sentinel-2 satellite data, human impact, and topographical features as predictors, explained 35% of the variance in community turnover. Elevation, sand percentage, and photosynthetic activities (NDVI32) were the top predictors. In addition to this signal of environmental filtering, we found a positive relationship between environmentally predicted families and their numbers of biotic interactions, suggesting environmental change could have a disproportionate effect on community networks. Together, these analyses show that coupling eDNA with environmental predictors including remote sensing data has capacity to test proposed Essential Biodiversity Variables and create new landscape biodiversity baselines that span the tree of life.

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.

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.

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.

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.

Drivers of species turnover vary with species commonness for native and alien plants with different residence times.

Communities comprising alien species with different residence times are natural experiments allowing the assessment of drivers of community assembly over time. Stochastic processes (such as dispersal and fluctuating environments) should be the dominant factors structuring communities of exotic species with short residence times. In contrast, communities should become more similar, or systematically diverge, if they contain exotics with increasing resident times, due to the increasing importance of deterministic processes (such as environmental filtering). We use zeta diversity (the number of species shared by multiple assemblages) to explore the relationship between the turnover of native species and two categories of alien species with different residence times (archaeophytes [introduced between 4000 BC and 1500 AD] and neophytes [introduced after 1500 AD]) in a network of nature reserves in central Europe. By considering multiple assemblages simultaneously, zeta diversity allows us to determine the contribution of rare and widespread species to turnover. Specifically, we explore the relative effects of assembly processes representing isolation by distance, environmental filtering, and environmental stochasticity (fluctuating environments) on zeta diversity using Multi-Site Generalized Dissimilarity Modelling (MS-GDM). Four clusters of results emerged. First, stochastic processes for structuring plant assemblages decreased in importance with increasing residence time. Environmental stochasticity only affected species composition for neophytes, offering possibilities to predict the spread debt of recent invasions. Second, native species turnover was well explained by environmental filtering and isolation by distance, although these factors did not explain the turnover of archaeophytes and neophytes. Third, native and alien species compositions were only correlated for rare species, whereas turnover in widespread alien species was surprisingly unrelated to the composition of widespread native species. Site-specific approaches would therefore be more appropriate for the monitoring and management of rare alien species, whereas species-specific approaches would suit widespread species. Finally, the size difference of nature reserves influences not only native species richness, but also their richness-independent turnover. A network of reserves must therefore be designed and managed using a variety of approaches to enhance native diversity, while controlling alien species with different residence times and degrees of commonness.

Diatom endemism and taxonomic turnover: Assessment in high-altitude alpine lakes covering a large geographical range.

Diatoms are widely used as ecological indicators and show various degrees of endemism. Many studies that support the idea of endemic species integrate several climate zones, a variety of ecosystem types, and often focus on a global scale. Here, we investigated whether endemism could be detected when considering a homogeneous type of ecosystem in a single climate zone. We sampled stone biofilms at 40-50 cm depth in high-altitude lakes in the Alpine climate zone. A total of 149 samples were obtained from the French and Georgian mountains, two areas separated by ∼3000 km. Using Amplicon Sequence Variants derived from DNA metabarcoding, we assessed taxonomic turnover and Zeta-diversity (a measure of endemism). We ran haplotype networks and phylogenetic tests to measure geographical signal in the phylogenies of dominant taxa. The French and Georgian communities shared 51 % of species. Species that were not shared across both regions were mostly rare, and often not characteristic of lakes but of neighboring habitats instead. In contrast, at the sub-species level, 87 % of the genotypes showed restricted distributions. Whereas endemism was the rule at sub-species level, most species were shared across both French and Georgian lakes, suggesting that geographic barriers strongly limited dispersal at the sub-species level but not species level. Dominant species hosted higher levels of sub-specific diversity than rare species. In contrast to global-scale studies, we did not find any significant geographical structuring in the phylogeny of the investigated species. This could indicate ongoing dispersal at a frequency fast enough to prevent allopatric divergence, yet slow enough to prevent sharing most haplotypes between France and Georgia. These results have implications for biomonitoring: depending on the taxonomic level chosen, robust generic tools (species level) or tools dedicated to a region able to discriminate fine pressures differences (sub-species level) may be developed.

Two dominant forms of multisite similarity decline - Their origins and interpretation.

The number of species shared by two or more sites is a fundamental measure of spatial variation in species composition. As more sites are included in the comparison of species composition, the average number of species shared across them declines, with a rate increasingly dependent on only the most widespread species. In over 80% of empirical communities, models of decline in shared species across multiple sites (multisite similarity decline) follow one of two distinct forms. An exponential form is assumed to reflect stochastic assembly and a power law form niche-based sorting, yet these explanations are largely untested, and little is known of how the two forms arise in nature. Using simulations, we first show that the distribution of the most widespread species largely differentiates the two forms, with the power law increasingly favored where such species occupy more than ~75% of sites. We reasoned the less cosmopolitan distribution of widespread species within exponential communities would manifest as differences in community biodiversity properties, specifically more aggregated within-species distributions, less even relative abundance distributions, and weaker between-species spatial associations. We tested and largely confirmed these relationships using 80 empirical datasets, suggesting that the form of multisite similarity decline offers a basis to predict how landscape-scale loss or gain of widespread species is reflected in different local-scale community structures. Such understanding could, for example, be used to predict changes in local-scale competitive interactions following shifts in widespread species' distributions. We propose multiple explanations for the origin of exponential decline, including high among-site abiotic variation, sampling of highly specialized (narrow niche width) taxa, and strong dispersal limitation. We recommend these are evaluated as alternative hypotheses to stochastic assembly.

Rare, common, alien and native species follow different rules in an understory plant community.

Biological invasions are a leading threat to biodiversity globally. Increasingly, ecosystems experience multiple introductions, which can have significant effects on patterns of diversity. The way these communities assemble will depend partly on whether rare and common alien species respond to environmental predictors in the same manner as rare and common native species, but this is not well understood. To examine this question across four national parks in south-eastern Australia, we sampled the understory plant community of eucalypt-dominated dry forest subject to multiple plant introductions. The drivers of diversity and turnover in alien and native species of contrasting frequency of occurrence (low, intermediate, and high) were each tested individually. We found alien species diversity and turnover were both strongly associated with abiotic conditions (e.g., soil pH), while distance had little influence because of the greater extent of occurrence and more homogeneous composition of common aliens. In contrast, native species diversity was not associated with abiotic conditions and their turnover was as strongly influenced by distance as by abiotic conditions. In both alien and native species, however, the most important predictors of turnover changed with frequency of occurrence. Although local coexistence appears to be facilitated by life history trade-offs, species richness of aliens and natives was negatively correlated and native species might face greater competition in areas with more neutral soils (e.g., pH > ~5.5) where alien richness and relative frequency were both highest. We conclude that diversity and turnover in the generally more widespread alien species are mainly driven by species sorting along an environmental gradient associated with pH and nutrient availability, whereas turnover of native species is driven by more neutral processes associated with dispersal limitation. We show alien and native plant species respond to different environmental factors, as do rare and common species within each component.

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.

Termite Assemblage and Damage on Tree Trunks in Fast-Growing Teak Plantations of Different Age: A Case Study in West Java, Indonesia.

We conducted surveys of termite assemblages and tree trunks damaged by termites in teakwood (Tectona grandis L.f.) plantations. The surveys were conducted in five-, six-, and nine-year-old plantations. We used a standardized belt-transect to collect termites and build tree inventories. Data of collected termites at the genus and functional-group levels and termites' diversity between plantations were compared with their attack rate. The results showed that four genera of soil recycler termites belonging to groups IIf and III were present across the plantations. Distribution analysis suggested that termite communities might develop from a stochastic distribution to a nonrandom co-occurrence distribution over time. Diversity analysis showed an increased nestedness-resultant diversity contribution to the total dissimilarity over time. Observed attacks on tree trunks were superficial and limited to the outer bark, with group IIf as the main contributor. Furthermore, the level of damage done by termites to tree trunks was positively correlated with increases in the group IIf occupancy area and overgrown understory vegetation. Plantation management by maintaining an adequate understory might suppress termite attacks on fast-growing teakwood, although in the case of our study, termite attacks are inevitable when termites from group IIf were already present.

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.

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.

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.

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.