Towards a Canary Islands barcode database for soil biodiversity: revealing cryptic and unrecorded mite species diversity within insular soils.

Soil arthropod diversity contributes to a high proportion of the total biodiversity on Earth. However, most soil arthropods are still undescribed, hindering our understanding of soil functioning and global biodiversity estimations. Inventorying soil arthropods using conventional taxonomical approaches is particularly difficult and costly due to the great species richness, abundance and local-scale heterogeneity of mesofauna communities and the poor taxonomic background knowledge of most lineages. To alleviate this situation, we have designed and implemented a molecular barcoding framework adapted to soil fauna. This pipeline includes different steps, starting with a morphology-based selection of specimens which are imaged. Then, DNA is extracted non-destructively. Both images and voucher specimens are used to assign a taxonomic identification, based on morphology that is further checked for consistency with molecular information. Using this procedure, we studied 239 specimens of mites from the Canary Islands including representatives of Mesostigmata, Sarcoptiformes and Trombidiformes, of which we recovered barcode sequences for 168 specimens that were morphologically identified to 49 species, with nine specimens that could only be identified at the genus or family levels. Multiple species delimitation analyses were run to compare molecular delimitations with morphological identifications, including ASAP, mlPTP, BINs and 3% and 8% genetic distance thresholds. Additionally, a species-level search was carried out at the Biodiversity Databank of the Canary Islands (BIOTA) to evaluate the number of species in our dataset that were not previously recorded in the archipelago. In parallel, a sequence-level search of our sequences was performed against BOLD Systems. Our results reveal that multiple morphologically identified species correspond to different molecular lineages, which points to significant levels of unknown cryptic diversity within the archipelago. In addition, we evidenced that multiple species in our dataset constituted new records for the Canary Islands fauna and that the information for these lineages within online genetic repositories is very incomplete. Our study represents the first systematic effort to catalogue the soil arthropod mesofauna of the Canary Islands and establishes the basis for the Canary Islands Soil Biodiversity barcode database. This resource will constitute a step forward in the knowledge of these arthropods in a region of special interest.

Dispersal ability and niche breadth influence interspecific variation in spider abundance and occupancy.

The relationship between species local abundance and their regional distribution (occupancy) is one of the most extensively recognized and investigated patterns in ecology. While exceptions exist, the generally held model is that locally abundant species also tend to be more widespread geographically. However, there is only a limited understanding of both the mechanisms driving this relationship, and their scale dependency. Here we use occupancy and abundance data for 123 species of spider from across the Canary Islands to understand how both dispersal ability and niche breadth might mediate variation among species for local abundance and occupancy. We test the predictions that (i) dispersal ability explains variation among species for both abundance and occupancy, and (ii) species with a higher degree of habitat specialization, reflecting more limited niche breadth, will have both higher occupancy and abundance. We find no evidence within habitat patches for an effect of dispersal ability on either local abundance or site occupancy, while across all patches species with higher dispersal ability tend to occupy more sites. Species largely restricted to laurel forests have higher abundance than species with broader niche breadth, but similar occupancy. The study revealed that dispersal ability and niche breadth were significant predictors of the abundance-occupancy relationship, highlighting the importance of both factors for understanding patterns of abundance and occupancy among spider species.

Inferring the ecological and evolutionary determinants of community genetic diversity.

Understanding the relative contributions of ecological and evolutionary processes to the structuring of ecological communities is needed to improve our ability to predict how communities may respond to future changes in an increasingly human-modified world. Metabarcoding methods make it possible to gather population genetic data for all species within a community, unlocking a new axis of data to potentially unveil the origins and maintenance of biodiversity at local scales. Here, we present a new eco-evolutionary simulation model for investigating community assembly dynamics using metabarcoding data. The model makes joint predictions of species abundance, genetic variation, trait distributions and phylogenetic relationships under a wide range of parameter settings (e.g. high speciation/low dispersal or vice versa) and across a range of community states, from pristine and unmodified to heavily disturbed. We first demonstrate that parameters governing metacommunity and local community processes leave detectable signatures in simulated biodiversity data axes. Next, using a simulation-based machine learning approach we show that neutral and non-neutral models are distinguishable and that reasonable estimates of several model parameters within the local community can be obtained using only community-scale genetic data, while phylogenetic information is required to estimate those describing metacommunity dynamics. Finally, we apply the model to soil microarthropod metabarcoding data from the Troodos mountains of Cyprus, where we find that communities in widespread forest habitats are structured by neutral processes, while high-elevation and isolated habitats act as an abiotic filter generating non-neutral community structure. We implement our model within the ibiogen R package, a package dedicated to the investigation of island, and more generally community-scale, biodiversity using community-scale genetic data.

Correlation of Prenatal and Postnatal Diagnosis in Umbilical-Portal-Systemic Venous Shunts.

INTRODUCTION: Umbilical-portal-systemic venous shunts (UPSVS) are rare anomalies in the development of the fetal venous system. There are several postnatal and prenatal classifications of hepatic venous anomalies but the link between them is missing. We aimed to review the prenatal to postnatal diagnosis correlation in UPSVS at our center. METHODS: It is a retrospective study of patients diagnosed with UPSVS between 2019 and 2021 at our institution. Demographic, obstetric, genetic, and neonatal data were reviewed with special focus on prenatal and postnatal ultrasounds. RESULTS: A total of seven patients were diagnosed with UPSVS at a median of 24 (20-34) weeks of gestational age. All patients were male and 62% were Caucasian. None of the patients had chromosomopathies or cardiac anomalies. One patient had renal ectopia, another one had a single umbilical artery, and a third one suffered from intrauterine growth retardation. An umbilico-systemic shunt (USS) was found in two patients and a ductus venosus-systemic shunt (DVSS) in the rest. Patients with USS were diagnosed postnatally with intrahepatic portosystemic shunts. One of the DVSS patients was transferred to another hospital and the other four had normal postnatal cardiac ultrasounds, with normal abdominal ultrasounds in two patients and lack of postnatally abdominal control in the other two. All babies were found to be doing well at a median follow-up of 1 month (0-24). CONCLUSION: There is a knowledge gap in the natural history of UPSVS between fetal and neonatal life. Building bridges between prenatal and postnatal research is mandatory in order to understand these rare anomalies.

Community metabarcoding reveals the relative role of environmental filtering and spatial processes in metacommunity dynamics of soil microarthropods across a mosaic of montane forests.

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.

Collective and harmonized high throughput barcoding of insular arthropod biodiversity: Toward a Genomic Observatories Network for islands.

Current understanding of ecological and evolutionary processes underlying island biodiversity is heavily shaped by empirical data from plants and birds, although arthropods comprise the overwhelming majority of known animal species, and as such can provide key insights into processes governing biodiversity. Novel high throughput sequencing (HTS) approaches are now emerging as powerful tools to overcome limitations in the availability of arthropod biodiversity data, and hence provide insights into these processes. Here, we explored how these tools might be most effectively exploited for comprehensive and comparable inventory and monitoring of insular arthropod biodiversity. We first reviewed the strengths, limitations and potential synergies among existing approaches of high throughput barcode sequencing. We considered how this could be complemented with deep learning approaches applied to image analysis to study arthropod biodiversity. We then explored how these approaches could be implemented within the framework of an island Genomic Observatories Network (iGON) for the advancement of fundamental and applied understanding of island biodiversity. To this end, we identified seven island biology themes at the interface of ecology, evolution and conservation biology, within which collective and harmonized efforts in HTS arthropod inventory could yield significant advances in island biodiversity research.

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.

Toward global integration of biodiversity big data: a harmonized metabarcode data generation module for terrestrial arthropods.

Metazoan metabarcoding is emerging as an essential strategy for inventorying biodiversity, with diverse projects currently generating massive quantities of community-level data. The potential for integrating across such data sets offers new opportunities to better understand biodiversity and how it might respond to global change. However, large-scale syntheses may be compromised if metabarcoding workflows differ from each other. There are ongoing efforts to improve standardization for the reporting of inventory data. However, harmonization at the stage of generating metabarcode data has yet to be addressed. A modular framework for harmonized data generation offers a pathway to navigate the complex structure of terrestrial metazoan biodiversity. Here, through our collective expertise as practitioners, method developers, and researchers leading metabarcoding initiatives to inventory terrestrial biodiversity, we seek to initiate a harmonized framework for metabarcode data generation, with a terrestrial arthropod module. We develop an initial set of submodules covering the 5 main steps of metabarcode data generation: (i) sample acquisition; (ii) sample processing; (iii) DNA extraction; (iv) polymerase chain reaction amplification, library preparation, and sequencing; and (v) DNA sequence and metadata deposition, providing a backbone for a terrestrial arthropod module. To achieve this, we (i) identified key points for harmonization, (ii) reviewed the current state of the art, and (iii) distilled existing knowledge within submodules, thus promoting best practice by providing guidelines and recommendations to reduce the universe of methodological options. We advocate the adoption and further development of the terrestrial arthropod module. We further encourage the development of modules for other biodiversity fractions as an essential step toward large-scale biodiversity synthesis through harmonization.

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.

Dispersal ability and its consequences for population genetic differentiation and diversification.

Dispersal ability is known to influence geographical structuring of genetic variation within species, with a direct relationship between low vagility and population genetic structure, which can potentially give rise to allopatric speciation. However, our general understanding of the relationship between dispersal ability, population differentiation and lineage diversification is limited. To address this issue, we sampled mitochondrial DNA variation within lineages of beetles and spiders across the Canary Islands to explore the relationships between dispersal ability, differentiation within lineages and diversification. We found positive relationships between population genetic structure and diversification for both beetles and spiders. Comparisons between dispersive and non-dispersive lineages revealed significant differences for both lineage differentiation and diversification. For both taxa, non-dispersive lineages had stronger population genetic structure. Genus-level endemic species richness and proxies for diversification rate within genera were higher in non-dispersive taxa for both beetles and spiders. Comparisons of average and maximum node divergences within genera suggest that species turnover may be higher in non-dispersive genera. Our results reveal a model where dispersal limitation may shape the diversity of lineages across evolutionary timescales by positively influencing intraspecific and species diversity, moderated by higher extinction rates compared to more dispersive lineages.

Coming of age for COI metabarcoding of whole organism community DNA: Towards bioinformatic harmonisation.

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.

Efficacy and Safety of Taulinoplasty Compared with the Minimally Invasive Repair of Pectus Excavatum Approach to Correct Pectus Excavatum.

Background: Minimally invasive repair of pectus excavatum (MIRPE) technique is the current most used surgical method for pectus excavatum (PE) correction. To avoid MIRPE-required invasion of the mediastinum or pleural cavity, we developed taulinoplasty approach as an alternative option for the reduction of the sunken chest. Materials and Methods: This is a single-center unblind nonrandomized prospective pilot study, involving 26 MIRPE and 23 taulinoplasty patients. The primary safety endpoint was the presence of surgical complications, and the primary efficacy endpoint was Haller index measured postoperatively in taulinoplasty patients. Secondary endpoints were intensive care unit (ICU) and hospital length of stay (LOS), duration of surgical procedure, and postoperative pain management. Chi-square, Mann-Whitney, or Student "t" tests were used for comparison as appropriate. Results: There were no difference in median age (14.5 and 14 years), Haller index (4.63 and 4.17), or comorbidities. MIRPE and taulinoplasty procedures showed no difference regarding major or minor surgical complications. The efficacy endpoint of taulinoplasty was measured 5.5 (1.5-12) months after the procedure, exhibiting a significant reduction in the Haller index [4.17 (3.7-4.7) at baseline and 3.7 (2.9-4.1) postoperatively, P = .03]. Taulinoplasty required less surgical time (60.4 ± 15.5 versus 70.7 ± 15.9 minutes, P < .01); shorter ICU and hospital LOS (P < .01, respectively); and required fewer days of peridural, intravenous, and oral analgesia (P < .01, respectively). Conclusions: In this pilot study, taulinoplasty technique was as safe as the MIRPE approach and effective for the correction of PE, reducing surgical time, ICU and hospital LOS, and achieving better postoperative pain control.

Dispersal limitations and long-term persistence drive differentiation from haplotypes to communities within a tropical sky-island: Evidence from community metabarcoding.

Neutral theory proposes that dispersal stochasticity is one of the main drivers of local diversity. Haplotypes-level genetic variation can now be efficiently sampled from across whole communities, thus making it possible to test neutral predictions from the genetic to species-level diversity, and higher. However, empirical data is still limited, with the few studies to date coming from temperate latitudes. Here, we focus on a tropical mountain within the Transmexican Volcanic Belt to evaluate spatially fine-scale patterns of arthropod community assembly to understand the role of dispersal limitation and landscape features as drivers of diversity. We sampled whole-communities of arthropods for eight orders at a spatial scale ranging from 50 m to 19 km, using whole community metabarcoding. We explored multiple hierarchical levels, from individual haplotypes to lineages at 0.5, 1.5, 3, 5, and 7.5% similarity thresholds, to evaluate patterns of richness, turnover, and distance decay of similarity with isolation-by-distance and isolation-by-resistance (costs to dispersal given by landscape features) approaches. Our results showed that distance and altitude influence distance decay of similarity at all hierarchical levels. This holds for arthropod groups of contrasting dispersal abilities, but with different strength depending on the spatial scale. Our results support a model where local-scale differentiation mediated by dispersal constraints, combined with long-term persistence of lineages, is an important driver of diversity within tropical sky islands.

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

Connecting high-throughput biodiversity inventories: Opportunities for a site-based genomic framework for global integration and synthesis.

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.

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.

Irreversible habitat specialization does not constrain diversification in hypersaline water beetles.

Specialization to extreme environments is often considered an evolutionary dead-end, leading to irreversible adaptations and reduced evolvability. There is, however, mixed evidence for this macroevolutionary pattern, and limited data from speciose lineages. Here, we tested the effect of habitat specialization to hypersaline waters in the diversification rates of aquatic beetles of the genus Ochthebius (Coleoptera, Hydraenidae), using a molecular phylogeny with more than 50% of the 546 recognized species, including representatives of all but one of the nine recognized subgenera and 17 species groups. Phylogenies were built combining mitochondrial and nuclear genes, with the addition of 42 mitochondrial genomes. Using Bayesian methods of character reconstruction, we show that hypersaline tolerance is an irreversible ecological specialization that arose multiple times. Two lineages of Ochthebius experienced a significant increase in diversification rates, one of them inhabiting hypersaline waters, but there was no overall correlation with habitat or any significant decrease in diversification rates despite the irreversibility of hypersaline tolerance. Our study tested for the first time the impact of hypersaline habitat specialization on diversification rates, finding no support for it to be an evolutionary dead-end. On the contrary, multiple and ancient lineages fully adapted to these extreme osmotic conditions have persisted and diversified over a long evolutionary timescale.

Mitochondrial Metagenomics Reveals the Ancient Origin and Phylodiversity of Soil Mites and Provides a Phylogeny of the Acari.

High-throughput DNA methods hold great promise for phylogenetic analysis of lineages that are difficult to study with conventional molecular and morphological approaches. The mites (Acari), and in particular the highly diverse soil-dwelling lineages, are among the least known branches of the metazoan Tree-of-Life. We extracted numerous minute mites from soils in an area of mixed forest and grassland in southern Iberia. Selected specimens representing the full morphological diversity were shotgun sequenced in bulk, followed by genome assembly of short reads from the mixture, which produced >100 mitochondrial genomes representing diverse acarine lineages. Phylogenetic analyses in combination with taxonomically limited mitogenomes available publicly resulted in plausible trees defining basal relationships of the Acari. Several critical nodes were supported by ancestral-state reconstructions of mitochondrial gene rearrangements. Molecular calibration placed the minimum age for the common ancestor of the superorder Acariformes, which includes most soil-dwelling mites, to the Cambrian-Ordovician (likely within 455-552 Ma), whereas the origin of the superorder Parasitiformes was placed later in the Carboniferous-Permian. Most family-level taxa within the Acariformes were dated to the Jurassic and Triassic. The ancient origin of Acariformes and the early diversification of major extant lineages linked to the soil are consistent with a pioneering role for mites in building the earliest terrestrial ecosystems.