BackgroundUnderstanding the underlying processes shaping spatial patterns of genetic structure in free-ranging organisms is a central topic in evolutionary biology. Here, we aim to disentangle the relative importance of neutral (i.e. genetic drift) and local adaptation (i.e. ecological divergence) processes in the evolution of spatial genetic structure of the Morales grasshopper (Chorthippus saulcyi moralesi), a narrow-endemic taxon restricted to the Central Pyrenees. More specifically, we analysed range-wide patterns of genetic structure and tested whether they were shaped by geography (isolation-by-distance, IBD), topographic complexity and present and past habitat suitability models (isolation-by-resistance, IBR), and environmental dissimilarity (isolation-by-environment, IBE).ResultsDifferent clustering analyses revealed a deep genetic structure that was best explained by IBR based on topographic complexity. Our analyses did not reveal a significant role of IBE, a fact that may be due to low environmental variation among populations and/or consequence of other ecological factors not considered in this study are involved in local adaptation processes. IBR scenarios informed by current and past climate distribution models did not show either a significant impact on genetic differentiation after controlling for the effects of topographic complexity, which may indicate that they are not capturing well microhabitat structure in the present or the genetic signal left by dispersal routes defined by habitat corridors in the past.ConclusionsOverall, these results indicate that spatial patterns of genetic variation in our study system are primarily explained by neutral divergence and migration-drift equilibrium due to limited dispersal across abrupt reliefs, whereas environmental variation or spatial heterogeneity in habitat suitability associated with the complex topography of the region had no significant effect on genetic discontinuities after controlling for geography. Our study highlights the importance of considering a comprehensive suite of potential isolating mechanisms and analytical approaches in order to get robust inferences on the processes promoting genetic divergence of natural populations.Electronic supplementary materialThe online version of this article (doi:10.1186/s12862-016-0663-7) contains supplementary material, which is available to authorized users.

Anthropogenic habitat fragmentation has altered the distribution and population sizes in many organisms worldwide. For this reason, understanding the demographic and genetic consequences of this process is necessary to predict the fate of populations and establish management practices aimed to ensure their viability. In this study, we analyse whether the spatial configuration of remnant semi-natural habitat patches within a chronically fragmented landscape has shaped the patterns of genetic diversity and structure in the habitat-specialist esparto grasshopper (Ramburiella hispanica). In particular, we predict that agricultural lands constitute barriers to gene flow and hypothesize that fragmentation has restricted interpopulation dispersal and reduced local levels of genetic diversity. Our results confirmed the expectation that isolation and habitat fragmentation have reduced the genetic diversity of local populations. Landscape genetic analyses based on circuit theory showed that agricultural land offers ∽1000 times more resistance to gene flow than semi-natural habitats, indicating that patterns of dispersal are constrained by the spatial configuration of remnant patches of suitable habitat. Overall, this study shows that semi-natural habitat patches act as corridors for interpopulation gene flow and should be preserved due to the disproportionately large ecological function that they provide considering their insignificant area within these human-modified landscapes.

Although resolving phylogenetic relationships and establishing species limits are primary goals of systematics, these tasks remain challenging at both conceptual and analytical levels. Here, we integrated genomic and phenotypic data and employed a comprehensive suite of coalescent-based analyses to develop and evaluate competing phylogenetic and species delimitation hypotheses in a recent evolutionary radiation of grasshoppers (Chorthippus binotatus group) composed of two species and eight putative subspecies. To resolve the evolutionary relationships within this complex, we first evaluated alternative phylogenetic hypotheses arising from multiple schemes of genomic data processing and contrasted genetic-based inferences with different sources of phenotypic information. Second, we examined the importance of number of loci, demographic priors, number and kind of phenotypic characters and sex-based trait variation for developing alternative species delimitation hypotheses. The best-supported topology was largely compatible with phenotypic data and showed the presence of two clades corresponding to the nominative species groups, one including three well-resolved lineages and the other comprising a four-lineage polytomy and a well-differentiated sister taxon. Integrative species delimitation analyses indicated that the number of employed loci had little impact on the obtained inferences but revealed the higher power provided by an increasing number of phenotypic characters and the usefulness of assessing their phylogenetic information content and differences between sexes in among-taxa trait variation. Overall, our study highlights the importance of integrating multiple sources of information to test competing phylogenetic hypotheses and elucidate the evolutionary history of species complexes representing early stages of divergence where conflicting inferences are more prone to appear.

Understanding the factors promoting species formation is a major task in evolutionary research. Here, we employ an integrative approach to study the evolutionary history of the Californian scrub white oak species complex (genus Quercus). To infer the relative importance of geographical isolation and ecological divergence in driving the speciation process, we (i) analysed inter- and intraspecific patterns of genetic differentiation and employed an approximate Bayesian computation (ABC) framework to evaluate different plausible scenarios of species divergence. In a second step, we (ii) linked the inferred divergence pathways with current and past species distribution models (SDMs) and (iii) tested for niche differentiation and phylogenetic niche conservatism across taxa. ABC analyses showed that the most plausible scenario is the one considering the divergence of two main lineages followed by a more recent pulse of speciation. Genotypic data in conjunction with SDMs and niche differentiation analyses support that different factors (geography vs. environment) and modes of speciation (parapatry, allopatry and maybe sympatry) have played a role in the divergence process within this complex. We found no significant relationship between genetic differentiation and niche overlap, which probably reflects niche lability and/or that multiple factors, have contributed to speciation. Our study shows that different mechanisms can drive divergence even among closely related taxa representing early stages of species formation and exemplifies the importance of adopting integrative approaches to get a better understanding of the speciation process.

Disentangling the relative role of environmental filtering and spatial processes in driving metacommunity structure across mountainous regions remains challenging, as the way we quantify spatial connectivity in topographically and environmentally heterogeneous landscapes can influence our perception of which process predominates. More empirical data sets are required to account for taxon‐ and context‐dependency, but relevant research in understudied areas is often compromised by the taxonomic impediment. Here we used haplotype‐level community DNA metabarcoding, enabled by stringent filtering of amplicon sequence variants (ASVs), to characterize metacommunity structure of soil microarthropod assemblages across a mosaic of five forest habitats on the Troodos mountain range in Cyprus. We found similar β diversity patterns at ASV and species (OTU, operational taxonomic unit) levels, which pointed to a primary role of habitat filtering resulting in the existence of largely distinct metacommunities linked to different forest types. Within‐habitat turnover was correlated to topoclimatic heterogeneity, again emphasizing the role of environmental filtering. However, when integrating landscape matrix information for the highly fragmented Quercus alnifolia habitat, we also detected a major role of spatial isolation determined by patch connectivity, indicating that stochastic and niche‐based processes synergistically govern community assembly. Alpha diversity patterns varied between ASV and OTU levels, with OTU richness decreasing with elevation and ASV richness following a longitudinal gradient, potentially reflecting a decline of genetic diversity eastwards due to historical pressures. Our study demonstrates the utility of haplotype‐level community metabarcoding for characterizing metacommunity structure of complex assemblages and improving our understanding of biodiversity dynamics across mountainous landscapes worldwide.

Metabarcoding of DNA extracted from community samples of whole organisms (whole organism community DNA, wocDNA) is increasingly being applied to terrestrial, marine and freshwater metazoan communities to provide rapid, accurate and high resolution data for novel molecular ecology research. The growth of this field has been accompanied by considerable development that builds on microbial metabarcoding methods to develop appropriate and efficient sampling and laboratory protocols for whole organism metazoan communities. However, considerably less attention has focused on ensuring bioinformatic methods are adapted and applied comprehensively in wocDNA metabarcoding. In this study we examined over 600 papers and identified 111 studies that performed COI metabarcoding of wocDNA. We then systematically reviewed the bioinformatic methods employed by these papers to identify the state‐of‐the‐art. Our results show that the increasing use of wocDNA COI metabarcoding for metazoan diversity is characterised by a clear absence of bioinformatic harmonisation, and the temporal trends show little change in this situation. The reviewed literature showed (i) high heterogeneity across pipelines, tasks and tools used, (ii) limited or no adaptation of bioinformatic procedures to the nature of the COI fragment, and (iii) a worrying underreporting of tasks, software and parameters. Based upon these findings we propose a set of recommendations that we think the metabarcoding community should consider to ensure that bioinformatic methods are appropriate, comprehensive and comparable. We believe that adhering to these recommendations will improve the long‐term integrative potential of wocDNA COI metabarcoding for biodiversity science.

High‐throughput sequencing (HTS) is increasingly being used for the characterization and monitoring of biodiversity. If applied in a structured way, across broad geographical scales, it offers the potential for a much deeper understanding of global biodiversity through the integration of massive quantities of molecular inventory data generated independently at local, regional and global scales. The universality, reliability and efficiency of HTS data can potentially facilitate the seamless linking of data among species assemblages from different sites, at different hierarchical levels of diversity, for any taxonomic group and regardless of prior taxonomic knowledge. However, collective international efforts are required to optimally exploit the potential of site‐based HTS data for global integration and synthesis, efforts that at present are limited to the microbial domain. To contribute to the development of an analogous strategy for the nonmicrobial terrestrial domain, an international symposium entitled “Next Generation Biodiversity Monitoring” was held in November 2019 in Nicosia (Cyprus). The symposium brought together evolutionary geneticists, ecologists and biodiversity scientists involved in diverse regional and global initiatives using HTS as a core tool for biodiversity assessment. In this review, we summarize the consensus that emerged from the 3‐day symposium. We converged on the opinion that an effective terrestrial Genomic Observatories network for global biodiversity integration and synthesis should be spatially led and strategically united under the umbrella of the metabarcoding approach. Subsequently, we outline an HTS‐based strategy to collectively build an integrative framework for site‐based biodiversity data generation.

The integration of genetic information with ecological and phenotypic data constitutes an effective approach to gain insight into the mechanisms determining interpopulation variability and the evolutionary processes underlying local adaptation and incipient speciation. Here, we use the Pyrenean Morales grasshopper (Chorthippus saulcyi moralesi) as study system to (i) analyse the relative role of genetic drift and selection in range-wide patterns of phenotypic differentiation and (ii) identify the potential selective agents (environment, elevation) responsible for variation. We also test the hypothesis that (iii) the development of dispersal-related traits is associated with different parameters related to population persistence/turnover, including habitat suitability stability over the last 120 000 years, distance to the species distribution core and population genetic variability. Our results indicate that selection shaped phenotypic differentiation across all the studied morphological traits (body size, forewing length and shape). Subsequent analyses revealed that among-population differentiation in forewing length was significantly explained by a temperature gradient, suggesting an adaptive response to thermoregulation or flight performance under contrasting temperature regimes. We found support for our hypothesis predicting a positive association between the distance to the species distribution core and the development of dispersal-related morphology, which suggests an increased dispersal capability in populations located at range edges that, in turn, exhibit lower levels of genetic variability. Overall, our results indicate that range-wide patterns of phenotypic variation are partially explained by adaptation in response to local environmental conditions and differences in habitat persistence between core and peripheral populations.