Metabarcoding for biodiversity inventory blind spots: A test case using the beetle fauna of an insular cloud forest.

Soils harbour a rich arthropod fauna, but many species are still not formally described (Linnaean shortfall) and the distribution of those already described is poorly understood (Wallacean shortfall). Metabarcoding holds much promise to fill this gap, however, nuclear copies of mitochondrial genes, and other artefacts lead to taxonomic inflation, which compromise the reliability of biodiversity inventories. Here, we explore the potential of a bioinformatic approach to jointly "denoise" and filter nonauthentic mitochondrial sequences from metabarcode reads to obtain reliable soil beetle inventories and address open questions in soil biodiversity research, such as the scale of dispersal constraints in different soil layers. We sampled cloud forest arthropod communities from 49 sites in the Anaga peninsula of Tenerife (Canary Islands). We performed whole organism community DNA (wocDNA) metabarcoding, and built a local reference database with COI barcode sequences of 310 species of Coleoptera for filtering reads and the identification of metabarcoded species. This resulted in reliable haplotype data after considerably reducing nuclear mitochondrial copies and other artefacts. Comparing our results with previous beetle inventories, we found: (i) new species records, potentially representing undescribed species; (ii) new distribution records, and (iii) validated phylogeographic structure when compared with traditional sequencing approaches. Analyses also revealed evidence for higher dispersal constraint within deeper soil beetle communities, compared to those closer to the surface. The combined power of barcoding and metabarcoding contribute to mitigate the important shortfalls associated with soil arthropod diversity data, and thus address unresolved questions for this vast biodiversity fraction.

Flightlessness in insects enhances diversification and determines assemblage structure across whole communities.

Dispersal limitation has been recurrently suggested to shape both macroecological patterns and microevolutionary processes within invertebrates. However, because of potential interactions among biological, environmental, temporal, and spatial variables, causal links among flight-related traits, diversification and spatial patterns of community assembly remain elusive. Integrating genetic variation within species across whole insect assemblages, within a simplified spatial and environmental framework, can be used to reduce the impact of these potentially confounding variables. Here, we used standardized sampling and mitochondrial DNA sequencing for a whole-community characterization of the beetle fauna inhabiting a singular forested habitat (laurel forest) within an oceanic archipelago setting (Canary Islands). The spatial structure of species assemblages together with species-level genetic diversity was compared at the archipelago and island scales for 104 winged and 110 wingless beetle lineages. We found that wingless beetle lineages have: (i) smaller range sizes at the archipelago scale, (ii) lower representation in younger island communities, (iii) stronger population genetic structure, and (iv) greater spatial structuring of species assemblages between and within islands. Our results reveal that dispersal limitation is a fundamental trait driving diversity patterns at multiple hierarchical levels by promoting spatial diversification and affecting the spatial configuration of entire assemblages at both island and archipelago scales.

Validated removal of nuclear pseudogenes and sequencing artefacts from mitochondrial metabarcode data.

Metabarcoding of Metazoa using mitochondrial genes may be confounded by both the accumulation of PCR and sequencing artefacts and the co-amplification of nuclear mitochondrial pseudogenes (NUMTs). The application of read abundance thresholds and denoising methods is efficient in reducing noise accompanying authentic mitochondrial amplicon sequence variants (ASVs). However, these procedures do not fully account for the complex nature of concomitant sequences and the highly variable DNA contribution of specimens in a metabarcoding sample. We propose, as a complement to denoising, the metabarcoding Multidimensional Abundance Threshold Evaluation (metaMATE) framework, a novel approach that allows comprehensive examination of multiple dimensions of abundance filtering and the evaluation of the prevalence of unwanted concomitant sequences in denoised metabarcoding datasets. metaMATE requires a denoised set of ASVs as input, and designates a subset of ASVs as being either authentic (mitochondrial DNA haplotypes) or nonauthentic ASVs (NUMTs and erroneous sequences) by comparison to external reference data and by analysing nucleotide substitution patterns. metaMATE (i) facilitates the application of read abundance filtering strategies, which are structured with regard to sequence library and phylogeny and applied for a range of increasing abundance threshold values, and (ii) evaluates their performance by quantifying the prevalence of nonauthentic ASVs and the collateral effects on the removal of authentic ASVs. The output from metaMATE facilitates decision-making about required filtering stringency and can be used to improve the reliability of intraspecific genetic information derived from metabarcode data. The framework is implemented in the metaMATE software (available at https://github.com/tjcreedy/metamate).

Long-term cloud forest response to climate warming revealed by insect speciation history.

Montane cloud forests are areas of high endemism, and are one of the more vulnerable terrestrial ecosystems to climate change. Thus, understanding how they both contribute to the generation of biodiversity, and will respond to ongoing climate change, are important and related challenges. The widely accepted model for montane cloud forest dynamics involves upslope forcing of their range limits with global climate warming. However, limited climate data provides some support for an alternative model, where range limits are forced downslope with climate warming. Testing between these two models is challenging, due to the inherent limitations of climate and pollen records. We overcome this with an alternative source of historical information, testing between competing model predictions using genomic data and demographic analyses for a species of beetle tightly associated to an oceanic island cloud forest. Results unequivocally support the alternative model: populations that were isolated at higher elevation peaks during the Last Glacial Maximum are now in contact and hybridizing at lower elevations. Our results suggest that genomic data are a rich source of information to further understand how montane cloud forest biodiversity originates, and how it is likely to be impacted by ongoing climate change.

The limited spatial scale of dispersal in soil arthropods revealed with whole-community haplotype-level metabarcoding.

Soil arthropod communities are highly diverse and critical for ecosystem functioning. However, our knowledge of spatial structure and the underlying processes of community assembly are scarce, hampered by limited empirical data on species diversity and turnover. We implement a high-throughput sequencing approach to generate comparative data for thousands of arthropods at three hierarchical levels: genetic, species and supra-specific lineages. A joint analysis of the spatial arrangement across these levels can reveal the predominant processes driving the variation in biological assemblages at the local scale. This multihierarchical approach was performed using haplotype-level COI metabarcoding of entire communities of mites, springtails and beetles from three Iberian mountain regions. Tens of thousands of specimens were extracted from deep and superficial soil layers and produced comparative phylogeographic data for >1,000 codistributed species and nearly 3,000 haplotypes. Local assemblage composition differed greatly between grasslands and forests and, within each habitat, showed strong spatial structure and high endemicity. Distance decay was high at all levels, even at the scale of a few kilometres or less. The local distance decay patterns were self-similar for the haplotypes and higher hierarchical entities, and this fractal structure was similar in all regions, suggesting that uniform processes of limited dispersal determine local-scale community assembly. Our results from whole-community metabarcoding provide insight into how dispersal limitations constrain mesofauna community structure within local spatial settings over evolutionary timescales. If generalized across wider areas, the high turnover and endemicity in the soil locally may indicate extremely high richness globally, challenging our current estimations of total arthropod diversity on Earth.

Building a global genomics observatory: Using GEOME (the Genomic Observatories Metadatabase) to expedite and improve deposition and retrieval of genetic data and metadata for biodiversity research.

Genetic data represent a relatively new frontier for our understanding of global biodiversity. Ideally, such data should include both organismal DNA-based genotypes and the ecological context where the organisms were sampled. Yet most tools and standards for data deposition focus exclusively either on genetic or ecological attributes. The Genomic Observatories Metadatabase (GEOME: geome-db.org) provides an intuitive solution for maintaining links between genetic data sets stored by the International Nucleotide Sequence Database Collaboration (INSDC) and their associated ecological metadata. GEOME facilitates the deposition of raw genetic data to INSDCs sequence read archive (SRA) while maintaining persistent links to standards-compliant ecological metadata held in the GEOME database. This approach facilitates findable, accessible, interoperable and reusable data archival practices. Moreover, GEOME enables data management solutions for large collaborative groups and expedites batch retrieval of genetic data from the SRA. The article that follows describes how GEOME can enable genuinely open data workflows for researchers in the field of molecular ecology.

Climate drives community-wide divergence within species over a limited spatial scale: evidence from an oceanic island.

Geographic isolation substantially contributes to species endemism on oceanic islands when speciation involves the colonisation of a new island. However, less is understood about the drivers of speciation within islands. What is lacking is a general understanding of the geographic scale of gene flow limitation within islands, and thus the spatial scale and drivers of geographical speciation within insular contexts. Using a community of beetle species, we show that when dispersal ability and climate tolerance are restricted, microclimatic variation over distances of only a few kilometres can maintain strong geographic isolation extending back several millions of years. Further to this, we demonstrate congruent diversification with gene flow across species, mediated by Quaternary climate oscillations that have facilitated a dynamic of isolation and secondary contact. The unprecedented scale of parallel species responses to a common environmental driver for evolutionary change has profound consequences for understanding past and future species responses to climate variation.

Dispersal limitation: Evolutionary origins and consequences in arthropods.

Niche and dispersal ability are key traits for explaining the geographical structuring of species into discrete populations, and its evolutionary significance. Beyond their individual effects, the interplay between species niche and its geographic limits, together with the evolutionary lability of dispersal ability, can underpin trait diversification and speciation when exposed to gradients of selection. In this issue of Molecular Ecology, two complementary papers demonstrate how evolutionary lability for dispersal ability linked to niche shift can drive such a model in a context that includes selection. Both papers investigate the evolution of dispersal limitation in arthropods across altitudinal gradients, but using taxa with contrasting ecologies. McCulloch et al. (2019) investigate the evolution of wing loss at higher altitudes in stoneflies, a taxon inhabiting freshwater systems. Suzuki et al. (2019) report a similar phenomenon, but involving wing reduction at higher altitudes in scorpionflies, a taxon associated with moist terrestrial habitats. Here, we compare and contrast the results of both studies to explore their broader implications for understanding diversification and speciation within arthropods.