A Good Compromise: Rapid and Robust Species Proxies for Inventorying Biodiversity Hotspots Using the Terebridae (Gastropoda: Conoidea)

Devising a reproducible approach for species delimitation of hyperdiverse groups is an ongoing challenge in evolutionary biology. Speciation processes combine modes of passive and adaptive trait divergence requiring an integrative taxonomy approach to accurately generate robust species hypotheses. However, in light of the rapid decline of diversity on Earth, complete integrative approaches may not be practical in certain species-rich environments. As an alternative, we applied a two-step strategy combining ABGD (Automated Barcode Gap Discovery) and Klee diagrams, to balance speed and accuracy in producing primary species hypotheses (PSHs). Specifically, an ABGD/Klee approach was used for species delimitation in the Terebridae, a neurotoxin-producing marine snail family included in the Conoidea. Delimitation of species boundaries is problematic in the Conoidea, as traditional taxonomic approaches are hampered by the high levels of variation, convergence and morphological plasticity of shell characters. We used ABGD to analyze gaps in the distribution of pairwise distances of 454 COI sequences attributed to 87 morphospecies and obtained 98 to 125 Primary Species Hypotheses (PSHs). The PSH partitions were subsequently visualized as a Klee diagram color map, allowing easy detection of the incongruences that were further evaluated individually with two other species delimitation models, General Mixed Yule Coalescent (GMYC) and Poisson Tree Processes (PTP). GMYC and PTP results confirmed the presence of 17 putative cryptic terebrid species in our dataset. The consensus of GMYC, PTP, and ABGD/Klee findings suggest the combination of ABGD and Klee diagrams is an effective approach for rapidly proposing primary species proxies in hyperdiverse groups and a reliable first step for macroscopic biodiversity assessment.     

Algorithmic single‐locus species delimitation: effects of sampling effort,variation and nonmonophyly in four methods and 1870 species of beetles

The vast number of undescribed species and the fast rate of biodiversity loss call for new approaches to speed up alpha taxonomy. A plethora of methods for delimiting species or operational taxonomic units (OTUs) based on sequence data have been published in recent years. We test the ability of four delimitation methods (BIN, ABGD, GMYC, PTP) to reproduce established species boundaries on a carefully curated DNA barcode data set of 1870 North European beetle species. We also explore how sampling effort, intraspecific variation, nearest neighbour divergence and nonmonophyly affect the OTU delimitations. All methods produced approximately 90% identity between species and OTUs. The effects of variation and sampling differed between methods. ABGD was sensitive to singleton sequences, while GMYC showed tendencies for oversplitting. The best fit between species and OTUs was achieved using simple rules to find consensus between discordant OTU delimitations. Using several approaches simultaneously allows the methods to compensate for each other's weaknesses. Barcode‐based OTU‐picking is an efficient way to delimit putative species from large data sets where the use of more sophisticated methods based on multilocus or genomic data is not feasible.     

Range‐constrained co‐occurrence simulation reveals little niche partitioning among rock‐dwelling Montenegrina land snails (Gastropoda: Clausiliidae)

Aim

Taxon co‐occurrence analysis is commonly used in ecology, but it has not been applied to range‐wide distribution data of partly allopatric taxa because existing methods cannot differentiate between distribution‐related effects and taxon interactions. Our first aim was to develop a taxon co‐occurrence analysis method that is also capable of taking into account the effect of species ranges and can handle faunistic records from museum databases or biodiversity inventories. Our second aim was to test the independence of taxon co‐occurrences of rock‐dwelling gastropods at different taxonomic levels, with a special focus on the Clausiliidae subfamily Alopiinae, and in particular the genus Montenegrina.

Location

Balkan Peninsula in south‐eastern Europe (46N–36N, 13.5E–28E).

Methods

We introduced a taxon‐specific metric that characterizes the occurrence probability at a given location. This probability was calculated as a distance‐weighted mean of the taxon's presence and absence records at all sites. We applied corrections to account for the biases introduced by varying sampling intensity in our dataset. Then we used probabilistic null‐models to simulate taxon distributions under the null hypothesis of no taxon interactions and calculated pairwise and cumulated co‐occurrences. Independence of taxon occurrences was tested by comparing observed co‐occurrences to simulated values.

Results

We observed significantly fewer co‐occurrences among species and intra‐generic lineages of Montenegrina than expected under the assumption of no taxon interaction.

Main conclusions

Fewer than expected co‐occurrences among species and intra‐generic clades indicate that species divergence preceded niche partitioning. This suggests a primary role of non‐adaptive processes in the speciation of rock‐dwelling gastropods. The method can account for the effects of distributional constraints in range‐wide datasets, making it suitable for testing ecological, biogeographical, or evolutionary hypotheses where interactions of partly allopatric taxa are in question.     

Using DNA taxonomy to investigate the ecological determinants of plankton diversity: explaining the occurrence of Synchaeta spp. (Rotifera,Monogononta) in mountain lakes

1. The occurrence of unresolved complexes of cryptic species may hinder the identification of the main ecological drivers of biodiversity when different cryptic taxa have different ecological requirements. 2. We assessed factors influencing the occurrence of Synchaeta species (monogonont rotifers) in 17 waterbodies of the Trentino‐South Tyrol region in the Eastern Alps. To do so, we compared the results of using unresolved complexes of cryptic species, as is common practice in limnological studies based on morphological taxonomy, and having resolved cryptic complexes, made possible by DNA taxonomy. 3. To identify cryptic species, we used the generalised mixed Yule coalescent (GMYC) model. We investigated the relationship between the environment and the occurrence of Synchaeta spp. by multivariate ordination using two definitions of the units of diversity, namely (i) unresolved species complexes (morphospecies) and (ii) putative cryptic species (GMYC entities). Our expectation was that resolving complexes of cryptic species could provide more information than using morphospecies. 4. As expected, DNA taxonomy provided greater taxonomic resolution than morphological taxonomy. Further, environmental‐based multivariate ordination on cryptic species explained a significantly higher proportion of variance than that based on morphospecies. Occurrence of GMYC entities was related to total phosphorus (TP), whereas no relationship could be found between morphospecies and the environment. Moreover, different cryptic species within the same morphospecies showed different, and even opposite, preferences for TP. In addition, the wide geographical distribution of haplotypes and cryptic species indicated the absence of barriers to dispersal in Synchaeta.     

A unified model explains commonness and rarity on coral reefs

Abundance patterns in ecological communities have important implications for biodiversity maintenance and ecosystem functioning. However, ecological theory has been largely unsuccessful at capturing multiple macroecological abundance patterns simultaneously. Here, we propose a parsimonious model that unifies widespread ecological relationships involving local aggregation, species‐abundance distributions, and species associations, and we test this model against the metacommunity structure of reef‐building corals and coral reef fishes across the western and central Pacific. For both corals and fishes, the unified model simultaneously captures extremely well local species‐abundance distributions, interspecific variation in the strength of spatial aggregation, patterns of community similarity, species accumulation, and regional species richness, performing far better than alternative models also examined here and in previous work on coral reefs. Our approach contributes to the development of synthetic theory for large‐scale patterns of community structure in nature, and to addressing ongoing challenges in biodiversity conservation at macroecological scales.     

An assessment of species limits of the South American mouse genus Oligoryzomys (Rodentia,Cricetidae) using unilocus delimitation methods

Oligoryzomys, as currently understood is formed by 25 living species, is the most diverse genus of the tribe Oryzomyini of the New World subfamily Sigmodontinae of cricetid rodents. Nonetheless, the species richness of Oligoryzomys seems to be an underestimate, given some species complex has been proposed in previous studies, at the time that large geographic areas remain to be sampled, and several taxonomic forms have not been assessed with contemporary approaches. In this study, we present a new assessment of the species diversity of Oligoryzomys based on multiple unilocus species delimitation methods (ABGD, BPP, PTP, GMYC and b GMYC), using 665 cytb gene sequences as evidence (532 gathered from Genbank and 133 obtained in this study). We sampled representatives of almost all currently known species of Oligoryzomys, at the time that extending the geographic coverage to the Central Andes, a large area that was largely unrepresented in previous studies. Phylogenetic relationships, based on a non‐redundant alignment, were inferred via maximum likelihood and Bayesian inference; an ultrametric tree, used in species delimitation analyses, was obtained using multiple secondary calibration points. Results of species delimitation methods are discussed at the light of previous knowledge (e.g., taxonomic history and geographic provenance of samples in relation to type localities) and the morphological assessments of some specimens. Results of the distinct delimitation methods are mostly congruent, being BPP and PTP the most sensible to estimate species delimitation, allowing us to suggest that Oligoryzomys is composed of 30 lineages of species level. Of these, 22 correspond to forms currently considered species; some of these include in their synonymy some forms currently considered valid species (e.g., yatesi would be a synonym of longicaudatus). The remaining eight lineages are candidate species that need to be further evaluated. This study, by advancing taxonomic hypothesis that should be further tested in future studies, constitutes a stepping‐stone for upcoming taxonomic and biogeographic studies centred on Oligoryzomys.     

Species delimitation and phylogeny in the genus Nasutitermes (Termitidae: Nasutitermitinae) in French Guiana

Species delimitation and identification can be arduous for taxa whose morphologic characters are easily confused, which can hamper global biodiversity assessments and pest species management. Exploratory methods of species delimitation that use DNA sequence as their primary information source to establish group membership and estimate putative species boundaries are useful approaches, complementary to traditional taxonomy. Termites of the genus Nasutitermes make interesting models for the application of such methods. They are dominant in Neotropical primary forests but also represent major agricultural and structural pests. Despite the prevalence, pivotal ecological role and economical impact of this group, the taxonomy of Nasutitermes species mainly depends on unreliable characters of soldier external morphology. Here, we generated robust species hypotheses for 79 Nasutitermes colonies sampled throughout French Guiana without any a priori knowledge of species affiliation. Sequence analysis of the mitochondrial cytochrome oxidase II gene was coupled with exploratory species‐delimitation tools, using the automatic barcode gap discovery method (ABGD) and a generalized mixed Yule‐coalescent model (GMYC) to propose primary species hypotheses (PSHs). PSHs were revaluated using phylogenetic analyses of two more loci (mitochondrial 16S rDNA and nuclear internal transcribed spacer 2) leading to 16 retained secondary species hypotheses (RSSH). Seven RSSHs, represented by 44/79 of the sampled colonies, were morphologically affiliated to species recognized as pests in the Neotropics, where they represent a real invasive pest potential in the context of growing ecosystem anthropization. Multigenic phylogenies based on combined alignments (1426–1784 bp) were also reconstructed to identify ancestral ecological niches and major‐pest lineages, revealing that Guyanese pest species do not form monophyletic groups.     

Influence of the geography of speciation on current patterns of coral reef fish biodiversity across the Indo‐Pacific

The role of speciation processes in shaping current biodiversity patterns represents a major scientific question for ecologists and biogeographers. Hence, numerous methods have been developed to determine the geography of speciation based on co‐occurrence between sister‐species. Most of these methods rely on the correlation between divergence time and several metrics based on the geographic ranges of sister‐taxa (i.e. overlap, asymmetry). The relationship between divergence time and these metrics has scarcely been examined in a spatial context beyond regression curves. Mapping this relationship across spatial grids, however, may unravel how speciation processes have shaped current biodiversity patterns through space and time. This can be particularly relevant for coral reef fishes of the Indo‐Pacific since the origin of the exceptional concentration of biodiversity in the Indo‐Australian Archipelago (IAA) has been actively debated, with several alternative hypotheses involving species diversification and dispersal. We reconstructed the phylogenetic relationships between three species‐rich families of coral reef fish (Chaetodontidae, Labridae, Pomacentridae) and calculated co‐occurrence metrics between closely related lineages of those families. We demonstrated that repeated biogeographic processes can be identified in present‐day species distribution by projecting co‐occurrence metrics between related lineages in a geographical context. Our study also evidence that sister‐species do not co‐occur randomly across the Indo‐Pacific, but tend to overlap their range within the IAA. We identified the imprint of two important biogeographic processes that caused this pattern in 48% of the sister‐taxa considered: speciation events within the IAA and repeated divergence between the Indian and Pacific Ocean, with subsequent secondary contact in the IAA.     

Spatial heterogeneity in the Mediterranean Biodiversity Hotspot affects barcoding accuracy of its freshwater fishes

Incomplete knowledge of biodiversity remains a stumbling block for conservation planning and even occurs within globally important Biodiversity Hotspots (BH). Although technical advances have boosted the power of molecular biodiversity assessments, the link between DNA sequences and species and the analytics to discriminate entities remain crucial. Here, we present an analysis of the first DNA barcode library for the freshwater fish fauna of the Mediterranean BH (526 spp.), with virtually complete species coverage (498 spp., 98% extant species). In order to build an identification system supporting conservation, we compared species determination by taxonomists to multiple clustering analyses of DNA barcodes for 3165 specimens. The congruence of barcode clusters with morphological determination was strongly dependent on the method of cluster delineation, but was highest with the general mixed Yule‐coalescent (GMYC) model‐based approach (83% of all species recovered as GMYC entity). Overall, genetic morphological discontinuities suggest the existence of up to 64 previously unrecognized candidate species. We found reduced identification accuracy when using the entire DNA‐barcode database, compared with analyses on databases for individual river catchments. This scale effect has important implications for barcoding assessments and suggests that fairly simple identification pipelines provide sufficient resolution in local applications. We calculated Evolutionarily Distinct and Globally Endangered scores in order to identify candidate species for conservation priority and argue that the evolutionary content of barcode data can be used to detect priority species for future IUCN assessments. We show that large‐scale barcoding inventories of complex biotas are feasible and contribute directly to the evaluation of conservation priorities.     

The community and ecosystem consequences of intraspecific diversity: a meta‐analysis

Understanding the relationships between biodiversity and ecosystem functioning has major implications. Biodiversity–ecosystem functioning relationships are generally investigated at the interspecific level, although intraspecific diversity (i.e. within‐species diversity) is increasingly perceived as an important ecological facet of biodiversity. Here, we provide a quantitative and integrative synthesis testing, across diverse plant and animal species, whether intraspecific diversity is a major driver of community dynamics and ecosystem functioning. We specifically tested (i) whether the number of genotypes/phenotypes (i.e. intraspecific richness) or the specific identity of genotypes/phenotypes (i.e. intraspecific variation) in populations modulate the structure of communities and the functioning of ecosystems, (ii) whether the ecological effects of intraspecific richness and variation are strong in magnitude, and (iii) whether these effects vary among taxonomic groups and ecological responses. We found a non‐linear relationship between intraspecific richness and community and ecosystem dynamics that follows a saturating curve shape, as observed for biodiversity–function relationships measured at the interspecific level. Importantly, intraspecific richness modulated ecological dynamics with a magnitude that was equal to that previously reported for interspecific richness. Our results further confirm, based on a database containing more than 50 species, that intraspecific variation also has substantial effects on ecological dynamics. We demonstrated that the effects of intraspecific variation are twice as high as expected by chance, and that they might have been underestimated previously. Finally, we found that the ecological effects of intraspecific variation are not homogeneous and are actually stronger when intraspecific variation is manipulated in primary producers than in consumer species, and when they are measured at the ecosystem rather than at the community level. Overall, we demonstrated that the two facets of intraspecific diversity (richness and variation) can both strongly affect community and ecosystem dynamics, which reveals the pivotal role of within‐species biodiversity for understanding ecological dynamics.     

Species co‐occurrence networks show reptile community reorganization under agricultural transformation

Agricultural transformation represents one of the greatest threats to biodiversity, causing degradation and loss of habitat, leading to changes in the richness and composition of communities. These changes in richness and composition may, in turn, lead to altered species co‐occurrence, but our knowledge of this remains limited. We used a novel co‐occurrence network approach to examine the impact of agricultural transformation on reptile community structure within two large (> 172 000 km²; 224 sites) agricultural regions in southeastern Australia. We contrasted assemblages from sites surrounded by intact and modified landscapes and tested four key hypotheses that agricultural transformation leads to (H1) declines in species richness, (H2) altered assemblages, (H3) declines in overall co‐occurrence, and (H4) complex restructuring of pairwise associations. We found that modified landscapes differed in composition but not richness compared with intact sites. Modified landscapes were also characterized by differences in co‐occurrence network structure; with species sharing fewer sites with each other (reduced co‐occurrence connectance), fewer highly‐connected species (truncation of the frequency distribution of co‐occurrence degree) and increased modularity of co‐occurrence networks. Critically, overall loss of co‐occurrence was underpinned by complex changes to the number and distribution of pair‐wise co‐occurrence links, with 41–44% of species also gaining associations with other species. Change in co‐occurrence was not correlated with changes in occupancy, nor by functional trait membership, allowing a novel classification of species susceptibility to agricultural transformation. Our study reveals the value of using co‐occurrence analysis to uncover impacts of agricultural transformation that may be masked in conventional studies of species richness and community composition.     

Phylogeny and co‐occurrence of mammal species on Southeast Asian islands

Aim Islands have often been used as model systems in community ecology. The incorporation of information on phylogenetic relatedness of species in studies of island assemblage structure is still uncommon, but could provide valuable insights into the processes of island community assembly. We propose six models of island community assembly that make different predictions about the associations between co‐occurrences of species pairs on islands, phylogenetic relatedness and ecological similarity. We then test these models using data on mammals of Southeast Asian islands. Location Two hundred and forty islands of the Sundaland region of Southeast Asia. Methods We quantified the co‐occurrence of species pairs on islands, and identified pairs that co‐occur more frequently (positive co‐occurrence) or less frequently (negative co‐occurrence) than expected under null models. We then examined the distributions of these significantly deviating pairs with respect to phylogenetic relatedness and ecological differentiation, and compared these patterns with those predicted by the six community assembly models. We used permutation regression to test whether co‐occurrence patterns are predicted by relatedness, body size difference or difference in diet quality. Separate co‐occurrence matrices were analysed in this way for seven mammal families and four smaller subsets of the islands of Sundaland. Results In many matrices, average numbers of negative co‐occurrences were higher than expected under null models. This is consistent with assemblage structuring by competition, but may also result from low geographic overlap of species pairs, which contributes to negative co‐occurrences at the archipelago‐wide level. Distributions of species pairs within plots of phylogenetic distance × ecological differentiation were consistent with competition, habitat filtering or within‐island speciation models, depending on the taxon. Regressions indicated that co‐occurrence was more likely among closely related species pairs within the Viverridae and Sciuridae, but in most matrices phylogenetic distance was unrelated to co‐occurrence. Main conclusions Simple deterministic models linking co‐occurrence with phylogeny and ecology are a useful framework for interpreting distributions and assemblage structure of island species. However, island assemblages in Sundaland have probably been shaped by a complex idiosyncratic set of interacting ecological and evolutionary processes, limiting the predictive power of such models.     

From null community to non-randomly structured actual plant assemblages: parsimony analysis of species co-occurrences

Null community is a spatio‐temporal abstraction of an initial regional species pool from which local species pools and actual community assemblages are organized. Any process that causes joint responses of species with similar susceptibilities affects community assembly. Through time, sequential assembly processes change the composition of a species pool in a way analogous to the one in which evolutionary processes promote character changes from an ancestor to current species. The segregation of species occurrences in an actual community suggests that assembly processes non‐randomly structured the observed community assemblages. However, going backwards to imply the causes of a particular arrangement of species is a non‐trivial challenge. I merge these premises with the philosophical and methodological foundations of cladistics. I propound parsimony analysis of species co‐occurrences as an outstanding means of devising operational hypotheses about the assembly of any non‐randomly structured set of actual community assemblages related to a common species pool. To explore this approach, I used field data gathered in a suite of 10 wetland assemblages. First, I tested independence of 101 plant species occurrences by a null model. As significant non‐random species co‐occurrence was detected, I applied a parsimony analysis taking the species occurrences as attributes, the assemblages as terminal units, and a putative null community constituted by all the present local species as the root of the assembly suite. The analysis produced four most parsimonious trees of assembly relationships. These trees maximize the number of similarities among community assemblages that can be explained by the sole fact of sharing a common regional species pool. One most parsimonious spatio‐temporal arrangement of species occurrence changes was reconstructed on one of the trees. I interpret this reconstruction in terms of assembly events, species exclusions and recruitments, showing the potentialities of this analysis to formulate operational hypotheses about community organization.     

Co‐occurrence is not evidence of ecological interactions

There is a rich amount of information in co‐occurrence (presence–absence) data that could be used to understand community assembly. This proposition first envisioned by Forbes (1907) and then Diamond (1975) prompted the development of numerous modelling approaches (e.g. null model analysis, co‐occurrence networks and, more recently, joint species distribution models). Both theory and experimental evidence support the idea that ecological interactions may affect co‐occurrence, but it remains unclear to what extent the signal of interaction can be captured in observational data. It is now time to step back from the statistical developments and critically assess whether co‐occurrence data are really a proxy for ecological interactions. In this paper, we present a series of arguments based on probability, sampling, food web and coexistence theories supporting that significant spatial associations between species (or lack thereof) is a poor proxy for ecological interactions. We discuss appropriate interpretations of co‐occurrence, along with potential avenues to extract as much information as possible from such data.     

Independent and interactive effects of two facilitators on their habitat‐providing host plant,Spartina alterniflora

The role of habitat‐providing species in facilitating associated species abundance and diversity is recognized as a key structuring force in many ecosystems. Reciprocal facilitation by associates, often involving multiple species, can be important for the maintenance of the host species. As with other multi‐species interactions (e.g. multiple predator effects), non‐additive relationships may be common among these associates, yet relatively few studies have examined potential interactions among multiple facilitator species. We combined field surveys and a mesocosm experiment to examine the independent and interactive effects of two co‐occurring facilitator species, ribbed mussels Geukensia demissa and fiddler crabs Uca pugilator, on their host salt marsh plant species, cordgrass Spartina alterniflora. We also experimentally examined how these relationships varied across different host plant genotypes. Overall, facilitator effects increased with increasing facilitator density. There was a significant interaction between mussel and fiddler crab presence, indicating that the effects of each species on cordgrass were dependent on the presence of the other facilitator species. In addition, there were strong interactions among mussels, fiddler crabs, and plant genotype, with greater variation in the performance of individual genotypes when fiddler crabs were absent. Our work reinforces the importance of considering multiple responses when assessing the functional redundancy of co‐occurring facilitators, as species are seldom completely redundant across the range of services they provide. It also highlights that the strength and direction of species interactions can vary due to genetic variation within the interacting species.     

Measuring size and composition of species pools: a comparison of dark diversity estimates

Ecological theory and biodiversity conservation have traditionally relied on the number of species recorded at a site, but it is agreed that site richness represents only a portion of the species that can inhabit particular ecological conditions, that is, the habitat‐specific species pool. Knowledge of the species pool at different sites enables meaningful comparisons of biodiversity and provides insights into processes of biodiversity formation. Empirical studies, however, are limited due to conceptual and methodological difficulties in determining both the size and composition of the absent part of species pools, the so‐called dark diversity. We used >50,000 vegetation plots from 18 types of habitats throughout the Czech Republic, most of which served as a training dataset and 1083 as a subset of test sites. These data were used to compare predicted results from three quantitative methods with those of previously published expert estimates based on species habitat preferences: (1) species co‐occurrence based on Beals' smoothing approach; (2) species ecological requirements, with envelopes around community mean Ellenberg values; and (3) species distribution models, using species environmental niches modeled by Biomod software. Dark diversity estimates were compared at both plot and habitat levels, and each method was applied in different configurations. While there were some differences in the results obtained by different methods, particularly at the plot level, there was a clear convergence, especially at the habitat level. The better convergence at the habitat level reflects less variation in local environmental conditions, whereas variation at the plot level is an effect of each particular method. The co‐occurrence agreed closest the expert estimate, followed by the method based on species ecological requirements. We conclude that several analytical methods can estimate species pools of given habitats. However, the strengths and weaknesses of different methods need attention, especially when dark diversity is estimated at the plot level.     

Using null model analysis of species co-occurrences to deconstruct biodiversity patterns and select indicator species

Aim Using total species richness to characterize biodiversity may mask multiple response patterns of species. We propose a null model analysis of species co‐occurrence‐based classification to identify sets of species that may have similar (within‐groups) and distinct (between groups) response patterns to their environment. The classification should also provide an explicit framework for selecting indicator species with characteristic co‐occurrence patterns to predict overall species richness. Location Côte‐Nord, Québec, Canada. Methods We combined null‐model of species co‐occurrence and cluster analysis to identify species groups within diverse assemblages of ground‐dwelling and flying beetles of stands in a boreal forest mosaic; we then examined their co‐occurrence and response patterns to habitat characteristics. Best subset regressions were used to select indicator species of richness within each group, from which indicators of total species richness were selected. Results The identified species groups appeared to display contrasting co‐occurrence and response patterns to at least one of the stand‐level habitat characteristics. Among flying beetles, for example, richness increased with stand‐level heterogeneity for two groups and decreased for two other groups, but the relationship was non‐significant for the total richness. We identified 28 indicator species that explained > 80% (validated by bootstrap analysis) of the variation in total species richness. Predictive performance of indicators was higher than when their co‐occurrence were reshuffled, even under a highly constrained null model, indicating that co‐occurrence patterns contributed to their predictive performance. Main conclusions Co‐occurrence‐based classification appears as a promising and effective tool for deconstructing biodiversity into species groups which reflect their ecological commonalities and differences, thus reducing the risk of making faulty inferences about the causes underlying overall diversity patterns. The method provides an explicit framework for selecting indicator species representing different species groups that may reflect the multiple responses of species co‐occurring with them. Indicator species can be effective for predicting overall species richness.     

Ephemeral ecological speciation and the latitudinal biodiversity gradient

The richness of biodiversity in the tropics compared to high‐latitude parts of the world forms one of the most globally conspicuous patterns in biology, and yet few hypotheses aim to explain this phenomenon in terms of explicit microevolutionary mechanisms of speciation and extinction. We link population genetic processes of selection and adaptation to speciation and extinction by way of their interaction with environmental factors to drive global scale macroecological patterns. High‐latitude regions are both cradle and grave with respect to species diversification. In particular, we point to a conceptual equivalence of “environmental harshness” and “hard selection” as eco‐evolutionary drivers of local adaptation and ecological speciation. By describing how ecological speciation likely occurs more readily at high latitudes, with such nascent species especially prone to extinction by fusion, we derive the ephemeral ecological speciation hypothesis as an integrative mechanistic explanation for latitudinal gradients in species turnover and the net accumulation of biodiversity.