Community assembly and metaphylogeography of soil biodiversity: Insights from haplotype-level community DNA metabarcoding within an oceanic island.

Most of our understanding of island diversity comes from the study of aboveground systems, while the patterns and processes of diversification and community assembly for belowground biotas remain poorly understood. Here, we take advantage of a relatively young and dynamic oceanic island to advance our understanding of ecoevolutionary processes driving community assembly within soil mesofauna. Using whole organism community DNA (wocDNA) metabarcoding and the recently developed metaMATE pipeline, we have generated spatially explicit and reliable haplotype-level DNA sequence data for soil mesofaunal assemblages sampled across the four main habitats within the island of Tenerife. Community ecological and metaphylogeographic analyses have been performed at multiple levels of genetic similarity, from haplotypes to species and supraspecific groupings. Broadly consistent patterns of local-scale species richness across different insular habitats have been found, whereas local insular richness is lower than in continental settings. Our results reveal an important role for niche conservatism as a driver of insular community assembly of soil mesofauna, with only limited evidence for habitat shifts promoting diversification. Furthermore, support is found for a fundamental role of habitat in the assembly of soil mesofauna, where habitat specialism is mainly due to colonization and the establishment of preadapted species. Hierarchical patterns of distance decay at the community level and metaphylogeographical analyses support a pattern of geographic structuring over limited spatial scales, from the level of haplotypes through to species and lineages, as expected for taxa with strong dispersal limitations. Our results demonstrate the potential for wocDNA metabarcoding to advance our understanding of biodiversity.

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).

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.

MtDNA metagenomics reveals large-scale invasion of belowground arthropod communities by introduced species.

Using a series of standardized sampling plots within forest ecosystems in remote oceanic islands, we reveal fundamental differences between the structuring of aboveground and belowground arthropod biodiversity that are likely due to large-scale species introductions by humans. Species of beetle and spider were sampled almost exclusively from single islands, while soil-dwelling Collembola exhibited more than tenfold higher species sharing among islands. Comparison of Collembola mitochondrial metagenomic data to a database of more than 80 000 Collembola barcode sequences revealed almost 30% of sampled island species are genetically identical, or near identical, to individuals sampled from often very distant geographic regions of the world. Patterns of mtDNA relatedness among Collembola implicate human-mediated species introductions, with minimum estimates for the proportion of introduced species on the sampled islands ranging from 45% to 88%. Our results call for more attention to soil mesofauna to understand the global extent and ecological consequences of species introductions.