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.

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.

Species delimitation and mitogenome phylogenetics in the subterranean genus Pseudoniphargus (Crustacea: Amphipoda).

The amphi-Atlantic distributions exhibited by many thalassoid stygobiont (obligate subterranean) crustaceans have been explained by fragmentation by plate tectonics of ancestral shallow water marine populations. The amphipod stygobiont genus Pseudoniphargus is distributed across the Mediterranean region but also in the North Atlantic archipelagos of Bermuda, Azores, Madeira and the Canaries. We used species delimitation methods and mitogenome phylogenetic analyses to clarify the species diversity and evolutionary relationships within the genus and timing their diversification. Analyses included samples from the Iberian Peninsula, northern Morocco, the Balearic, Canarian, Azores and Madeira archipelagoes plus Bermuda. In most instances, morphological and molecular-based species delimitation analyses yielded consistent results. Notwithstanding, in a few cases either incipient speciation with no involvement of detectable morphological divergence or species crypticism were the most plausible explanations for the disagreement found between morphological and molecular species delimitations. Phylogenetic analyses based on a robust calibrated mitochondrial tree suggested that Pseudoniphargus lineages have a younger age than for other thalassoid amphipods displaying a disjunct distribution embracing both sides of the Atlantic Ocean. A major split within the family was estimated to occur at the Paleocene, when a lineage from Northern Iberian Peninsula diverged from the rest of pseudoniphargids. Species diversification in the peri-Mediterranean area was deduced to occur in early Miocene to Tortonian times, while in the Atlantic islands it started in the Pliocene. Our results show that the current distribution pattern of Pseudoniphargus resulted from a complex admix of relatively ancient vicariance events and several episodes of long- distance dispersal.

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.

The imprint of geologic history on within-island diversification of woodlouse-hunter spiders (Araneae, Dysderidae) in the Canary Islands.

Geological processes and ecological adaptation are major drivers of diversification on oceanic islands. Although diversification in these islands is often interpreted as resulting from dispersal or island hopping rather than vicariance, this may not be the case in islands with complex geological histories. The island of Tenerife, in the Canary Islands, emerged in the late Miocene as 3 precursor islands that were subsequently connected and reisolated by volcanic cycles. The spider Dysdera verneaui is endemic to the island of Tenerife, where it is widely distributed throughout most island habitats, providing an excellent model to investigate the role of physical barriers and ecological adaptation in shaping within-island diversity. Here, we present evidence that the phylogeographic patterns of this species trace back to the independent emergence of the protoislands. Molecular markers (mitochondrial genes cox1, 16S, and nad1 and the nuclear genes ITS-2 and 28S) analyzed from 100 specimens (including a thorough sampling of D. verneaui populations and additional outgroups) identify 2 distinct evolutionary lineages that correspond to 2 precursor islands, each with diagnostic genital characters indicative of separate species status. Episodic introgression events between these 2 main evolutionary lineages explain the observed incongruence between mitochondrial and nuclear markers, probably as a result of the homogenization of their ITS-2 sequence types. The most widespread lineage exhibits a complex population structure, which is compatible with either secondary contact, following connection of deeply divergent lineages, or alternatively, a back colonization from 1 precursor island to another.

On hypogean Roncocreagris (Arachnida: Pseudoscorpiones: Neobisiidae) from Portugal, with descriptions of three new species.

Three new hypogean species of the Iberian genus Roncocreagris Mahnert, 1974 are described from mainland Portugal: R. borgesi sp. nov. and R. gepesi sp. nov. from caves in the Sicó massif, and R. occidentalis sp. nov. from caves in the Montejunto and Cesaredas karst plateau. This brings to nine the number of known hypogean species of the mostly Iberian genus Roncocreagris: five from Portugal and four from Spain. Ecological comments and new localities for some of the previously known species are also included.

Temperature variation in caves and its significance for subterranean ecosystems.

Climate change affects all ecosystems, but subterranean ecosystems are repeatedly neglected from political and public agendas. Cave habitats are home to unknown and endangered species, with low trait variability and intrinsic vulnerability to recover from human-induced disturbances. We studied the annual variability and cyclicity of temperatures in caves vis-à-vis surface in different climatic areas. We hypothesize that cave temperatures follow the average temperature pattern at the surface for each location with a slight delay in the signal, but we found three different thermal patterns occurring in caves: (1) high positive correlation and a similar thermal pattern to the surface, (2) low correlation and a slight thermal delay of the signal from the surface, and (3) high negative correlation with an extreme delay from the surface. We found daily thermal cycles in some caves, which may potentially control the circadian rhythms of cave organisms. Our results show that caves had lower thermal amplitude than the surface, and that thermal averages within caves approximately correspond to the to the annual average of surface temperature. Caves buffer external temperature and act as refugia for biota in extreme climatic events. Likewise, temperature increases at surface will lead to increment in caves, threatening subterranean biota and ecosystem services.

Under the volcano: phylogeography and evolution of the cave-dwelling Palmorchestia hypogaea (Amphipoda, Crustacea) at La Palma (Canary Islands).

BACKGROUND: The amphipod crustacean Palmorchestia hypogaea occurs only in La Palma (Canary Islands) and is one of the few terrestrial amphipods in the world that have adapted to a strictly troglobitic life in volcanic cave habitats. A surface-dwelling closely related species (Palmorchestia epigaea) lives in the humid laurel forest on the same island. Previous studies have suggested that an ancestral littoral Orchestia species colonized the humid forests of La Palma and that subsequent drought episodes in the Canaries reduced the distribution of P. epigaea favouring the colonization of lava tubes through an adaptive shift. This was followed by dispersal via the hypogean crevicular system. RESULTS: P. hypogaea and P. epigaea did not form reciprocally monophyletic mitochondrial DNA clades. They showed geographically highly structured and genetically divergent populations with current gene flow limited to geographically close surface locations. Coalescence times using Bayesian estimations assuming a non-correlated relaxed clock with a normal prior distribution of the age of La Palma, together with the lack of association of habitat type with ancestral and recent haplotypes, suggest that their adaptation to cave life is relatively ancient. CONCLUSION: The data gathered here provide evidence for multiple invasions of the volcanic cave systems that have acted as refuges. A re-evaluation of the taxonomic status of the extant species of Palmorchestia is needed, as the division of the two species by habitat and ecology is unnatural. The information obtained here, and that from previous studies on hypogean fauna, shows the importance of factors such as the uncoupling of morphological and genetic evolution, the role of climatic change and regressive evolution as key processes in leading to subterranean biodiversity.

A combined field survey and molecular identification protocol for comparing forest arthropod biodiversity across spatial scales.

Obtaining fundamental biodiversity metrics such as alpha, beta and gamma diversity for arthropods is often complicated by a lack of prior taxonomic information and/or taxonomic expertise, which can result in unreliable morphologically based estimates. We provide a set of standardized ecological and molecular sampling protocols that can be employed by researchers whose taxonomic skills may be limited, and where there may be a lack of robust a priori information regarding the regional pool of species. These protocols combine mass sampling of arthropods, classification of samples into parataxonomic units (PUs) and selective sampling of individuals for mtDNA sequencing to infer biological species. We sampled ten lowland rainforest plots located on the volcanic oceanic island of Réunion (Mascarene archipelago) for spiders, a group with limited taxonomic and distributional data for this region. We classified adults and juveniles into PUs and then demonstrated the reconciliation of these units with presumed biological species using mtDNA sequence data, ecological data and distributional data. Because our species assignment protocol is not reliant upon prior taxonomic information, or taxonomic expertise, it minimizes the problem of the Linnean shortfall to yield diversity estimates that can be directly compared across independent studies. Field sampling can be extended to other arthropod groups and habitats by adapting our field sampling protocol accordingly.

A geographical distribution database of the genus Dysdera in the Canary Islands (Araneae, Dysderidae).

The ground-dweller spider genus Dysdera shows very high species richness on the oceanic archipelago of the Canary Islands, providing one of the most outstanding examples of island radiation among spiders, only paralleled by Tetragnatha spiders on the Hawaiian archipelago. A georeferenced database of the 48 Dysdera species occurring in the Canary Islands was assembled to facilitate ongoing and future research on this remarkable lineage. All species are endemic to the archipelago except for the cosmopolitan Dysdera crocata. The dataset consists of 794 distributional records documented from 1971 to 2015, each locality being represented only once per species. Distribution maps are provided for each species, along with basic diversity and distribution information. The database and geographical maps included in this article stand for the most updated, accurate and complete information on the distribution of the spider genus Dysdera in the Canary Islands.

Diversification of the forest beetle genus Tarphius on the Canary Islands, and the evolutionary origins of island endemics.

The flightless beetle genus Tarphius Erichson (Coleoptera: Colydiidae) is a distinctive element of the beetle fauna of the Canary Islands with 29 species distributed across the five western islands. The majority of Tarphius species are rare and intimately associated with the monteverde forest and only two species occur on more than one island. In this study we investigate the phylogeography of the Canary Island Tarphius, and their relationship to Tarphius from the more northerly archipelagos of Madeira and the Azores using maximum parsimony and Bayesian inference analysis of mitochondrial cytochrome oxidase I and II sequence data. We use geological datings for the Canary Islands, Azores, and Madeira to calibrate specific nodes of the tree for the estimation of divergence times using a penalized likelihood method. Data suggest that the Canary Island species assemblage is of some antiquity, however, much of this species diversity is relatively recent in origin. The phylogenetic relationships of species inhabiting the younger islands of El Hierro and La Palma indicate that colonization events between islands have probably been a significant factor in the evolutionary history of the Canary Island species assemblage. A comparison of molecular phylogenetic studies of arthropods on the Canary Islands suggests that, in the evolution of the arthropod species community of an island, the origin of endemic species is initially the result of colonizing lineages differentiating from their source populations. However, as an island matures a greater proportion of endemic species originate from intra-island speciation.

Mitogenomic phylogenetic analysis supports continental-scale vicariance in subterranean thalassoid crustaceans.

Many continental subterranean water crustaceans ("stygobionts") display extreme disjunct distributions, where different species in the same genus are isolated on continents or islands separated by broad oceanic expanses. Despite their freshwater habitat, most of these taxa appear to be most closely related to typical marine groups ("thalassoid" origin). Among the hadzioids-thalassoid amphipods including the stygobiont families Hadziidae, Pseudoniphargidae, and Metacrangonyctidae-several genera are restricted to inland groundwaters ranging from the Caribbean region to the Mediterranean and Middle East, including interspersed oceanic islands. This distribution might have arisen from Tethyan vicariance triggered by the sequential occlusion of the former Tethys Sea, a vast circumtropical ocean existing from the Middle Jurassic up to 20 million years ago (mya). Previous studies have been based on morphological analyses or limited DNA sequence data, making it difficult to test this hypothesis. We used complete mitochondrial protein-coding gene sequences, mainly obtained by next-generation sequencing methods and a nuclear ribosomal gene to resolve the phylogeny and to establish a time frame for diversification of the family Metacrangonyctidae (Amphipoda). The results were consistent with the plate tectonics vicariance hypothesis, with major diversifications occurring between 96 and 83 mya.

Evolution and diversification of the forest and hypogean ground-beetle genus Trechus in the Canary Islands.

The beetle genus Trechus (Carabidae) is represented in the Macaronesian Islands by 43 endemic species. The Canary Islands have 16 endemic species, with two adapted to hypogean life. Phylogenetic relationships among 177 individuals of 38 Canarian, Madeiran, Azorean and continental Trechus species were examined using mitochondrial DNA and nuclear internal transcribed spacer 2 (ITS2) sequence data. Results show two main lineages in the Canaries: one comprising two sister groups with species from the laurel forest of La Gomera and Tenerife, and the other containing the single species from Gran Canaria and a species complex in the four western islands including two troglobites. Calibrations were applied to a linearized tree using a relaxed molecular clock method to estimate the major evolutionary divergence times of the Canarian Trechus species. Although the species assemblage in this archipelago is relatively ancient (7-8 million years), much of the species diversity is recent. Transition to the hypogean environment is more consistent with the "adaptive shift" rather than with the competing "climatic relict" hypothesis.

Living on the edge: demographic and phylogeographical patterns in the woodlouse-hunter spider Dysdera lancerotensis Simon, 1907 on the eastern volcanic ridge of the Canary Islands.

The Eastern Canary Islands are the emerged tips of a continuous volcanic ridge running parallel to the northeastern African coast, originated by episodic volcanic eruptions that can be traced back to the Miocene and that, following a major period of quiescence and erosion, continued from the Pliocene to the present day. The islands have been periodically connected by eustatic sea-level changes resulting from Pleistocene glacial cycles. The ground-dwelling spider Dysdera lancerotensis Simon, 1907 occurs along the entire ridge, except on recent barren lavas and sand dunes, and is therefore an ideal model for studying the effect of episodic geological processes on terrestrial organisms. Nested clade and population genetic analyses using 39 haplotypes from 605 base pairs of mitochondrial DNA cytochrome c oxidase I sequence data, along with phylogenetic analyses including two additional mitochondrial genes, uncover complex phylogeographical and demographic patterns. Our results indicate that D. lancerotensis colonized the ridge from north to south, in contrast to what had been expected given the SSW-NNE trend of volcanism and to what had been reported for other terrestrial arthropods. The occurrence of several episodes of extinction, recolonization and expansion are hypothesized for this species, and areas that act as refugia during volcanic cycles are identified. Relaxed molecular clock methods reveal divergence times between main haplotype lineages that suggest an older origin of the northern islets than anticipated based on geological evidence. This study supports the key role of volcanism in shaping the distribution of terrestrial organisms on oceanic islands and generates phylogeographical predictions that warrant further research into other terrestrial endemisms of this fascinating region.

Testing phylogeographic predictions on an active volcanic island: Brachyderes rugatus (Coleoptera: Curculionidae) on La Palma (Canary Islands).

Volcanic islands with well-characterized geological histories can provide ideal templates for generating and testing phylogeographic predictions. Many studies have sought to utilize these to investigate patterns of colonization and speciation within groups of closely related species across a number of islands. Here we focus attention within a single volcanic island with a well-characterized geological history to develop and test phylogeographic predictions. We develop phylogeographic predictions within the island of La Palma of the Canary Islands and test these using 69 haplotypes from 570 base pairs of mitochondrial DNA cytochrome oxidase II sequence data for 138 individuals of Brachyderes rugatus rugatus, a local endemic subspecies of curculionid beetle occurring throughout the island in the forests of Pinus canariensis. Although geological data do provide some explanatory power for the phylogeographic patterns found, our network-based analyses reveal a more complicated phylogeographic history than initial predictions generated from data on the geological history of the island. Reciprocal illumination of geological and phylogeographic history is also demonstrated with previous geological speculation gaining phylogeographic corroboration from our analyses.

Genetic structure, phylogeography and demography of two ground-beetle species endemic to the Tenerife laurel forest (Canary Islands).

The volcanic island of Tenerife (Canary archipelago) was formerly covered at 600-1200 m above sea level on most of its northern side by a cloud forest holding much of the endemic insect fauna. In the most significant surviving patches of this laurel forest at the eastern and western tips of the island occur two forest-specialist, closely related species of Eutrichopus (Coleoptera, Carabidae); here we present data on mitochondrial DNA variation among populations of these species. In total, 116 individuals from 16 localities were sampled and a 638 bp fragment of the cytochrome oxidase subunit II gene was sequenced, obtaining evidence for two distinct evolutionary lineages, in accordance with morphological and biogeographical data. Volcanic events at approximately 0.7 Ma might be responsible for vicariance and the fragmentation of the geographical range of an ancestral species, causing the establishment of two matrilineal lineages. Using nested clade and historical demography analyses we infer past cycles of demographic bottlenecks followed by population expansion, mostly in agreement with the geological time scale of volcanic events. Recent trends, however, refer to fragmentation of the cloud forest due to human intervention.

Phylogeography of the endangered darkling beetle species of Pimelia endemic to Gran Canaria (Canary Islands).

Phylogenetic and geographical nested clade analysis (NCA) methods were applied to mitochondrial DNA sequences of Pimelia darkling beetles (Coleoptera, Tenebrionidae) endemic to Gran Canaria, an island in the Canary archipelago. The three species P. granulicollis, P. estevezi and P. sparsa occur on the island, the latter with three recognized subspecies. Another species, P. fernandezlopezi (endemic to the island of La Gomera) is a close relative of P. granulicollis based on partial Cytochrome Oxidase I mtDNA sequences obtained in a previous study. Some of these beetles are endangered, so phylogeographical structure within species and populations can help to define conservation priorities. A total of about 700 bp of Cytochrome Oxidase II were examined in 18 populations and up to 75 individuals excluding outgroups. Among them, 22 haplotypes were exclusive to P. granulicollis and P. estevezi and 31 were from P. sparsa. Phylogenetic analysis points to the paraphyly of Gran Canarian Pimelia, as the La Gomera P. fernandezlopezi haplotypes are included in them, and reciprocal monophyly of two species groups: one constituted by P. granulicollis, P. estevezi and P. fernandezlopezi (subgenus Aphanaspis), and the other by P. sparsa'sensu lato'. The two species groups show a remarkably high mtDNA divergence. Within P. sparsa, different analyses all reveal a common result, i.e. conflict between current subspecific taxonomic designations and evolutionary units, while P. estevezi and P. fernandezlopezi are very close to P. granulicollis measured at the mtDNA level. Geographical NCA identifies several cases of nonrandom associations between haplotypes and geography that may be caused by allopatric fragmentation of populations with some cases of restriction of gene flow or range expansion. Analyses of molecular variance and geographical NCA allow definition of evolutionary units for conservation purposes in both species-groups and suggest scenarios in which vicariance caused by geological history of the island may have shaped the pattern of the mitochondrial genetic diversity of these beetles.