Population Genomics of Adaptive Radiations: Exceptionally High Levels of Genetic Diversity and Recombination in an Endemic Spider From the Canary Islands.

The spider genus Dysdera has undergone a remarkable diversification in the oceanic archipelago of the Canary Islands, with ~60 endemic species having originated during the 20 million years since the origin of the archipelago. This evolutionary radiation has been accompanied by substantial dietary shifts, often characterised by phenotypic modifications encompassing morphological, metabolic and behavioural changes. Hence, these endemic spiders represent an excellent model for understanding the evolutionary drivers and to pinpoint the genomic determinants underlying adaptive radiations. Recently, we achieved the first chromosome-level genome assembly of one of the endemic species, D. silvatica, providing a high-quality reference sequence for evolutionary genomics studies. Here, we conducted a low coverage-based resequencing study of a natural population of D. silvatica from La Gomera island. Taking advantage of the new high-quality genome, we characterised genome-wide levels of nucleotide polymorphism, divergence and linkage disequilibrium, and inferred the demographic history of this population. We also performed comprehensive genome-wide scans for recent positive selection. Our findings uncovered exceptionally high levels of nucleotide diversity and recombination in this geographically restricted endemic species, indicative of large historical effective population sizes. We also identified several candidate genomic regions that are potentially under positive selection, highlighting relevant biological processes, such as vision and nitrogen extraction as potential adaptation targets. These processes may ultimately drive species diversification in this genus. This pioneering study of spiders that are endemic to an oceanic archipelago lays the groundwork for broader population genomics analyses aimed at understanding the genetic mechanisms driving adaptive radiation in island ecosystems.

Jumping through hoops: Structural rearrangements and accelerated mutation rates on Dendrodorididae (Mollusca: Nudibranchia) mitogenomes rumble their evolution.

The systematics of the family Dendrodorididae, with only three valid genera, is a challenge for integrative taxonomists. Its members lack hard structures for morphological comparisons and their mitochondrial and nuclear markers provide contradictory phylogenetic signals, making phylogenetic reconstructions difficult. This molecular discordance has been hypothesized to be the result of nuclear pseudogenes or exogenous contamination. However, these hypotheses have not been tested. Here, we assembled the first genome drafts of seven Dendrodorididae species to investigate the evolutionary processes of this family. Two of the mitogenomes displayed an identical structural rearrangement involving the translocation of three coding genes and five tRNAs, described for the first time in nudibranchs. In addition, we found particularly high dN and dN/dS values and multiple insertions and deletions on the mitochondrial genes of smooth Dendrodoris. In contrast, nuclear single-copy ortholog genes showed no such mutational differences. Models of protein structures from mitochondrial genes are conserved, suggesting conserved functionality. Phylogenies using mitogenomic and nuclear data showed that species with rearranged mitogenomes form a clade, although Dendrodorididae relationships remained unresolved. The present study provides novel evidence for accelerated mutation rates in the mitogenomes of Dendrodorididae, which presumably have implications on respiratory adaptation, and highlights the importance of using genomic data to unveil rare evolutionary processes, crucial for correctly inferring phylogenies.

The Non-Dereliction in Evolution: Trophic Specialisation Drives Convergence in the Radiation of Red Devil Spiders (Araneae: Dysderidae) in the Canary Islands.

Natural selection plays a key role in deterministic evolution, as clearly illustrated by the multiple cases of repeated evolution of ecomorphological characters observed in adaptive radiations. Unlike most spiders, Dysdera species display a high variability of cheliceral morphologies, which has been suggested to reflect different levels of specialization to feed on isopods. In this study, we integrate geometric morphometrics and experimental trials with a fully resolved phylogeny of the highly diverse endemic species from the Canary Islands to 1) quantitatively delimit the different cheliceral morphotypes present in the archipelago, 2) test their association with trophic specialization, as reported for continental species, 3) reconstruct the evolution of these ecomorphs throughout the diversification of the group, 4) test the hypothesis of convergent evolution of the different morphotypes, and 5) examine whether specialization constitutes a case of evolutionary irreversibility in this group. We show the existence of 9 cheliceral morphotypes and uncovered their significance for trophic ecology. Further, we demonstrate that similar ecomorphs evolved multiple times in the archipelago, providing a novel study system to explain how convergent evolution and irreversibility due to specialization may be combined to shape phenotypic diversification in adaptive radiations.

A targeted gene phylogenetic framework to investigate diversification in the highly diverse yet geographically restricted red devil spiders (Araneae, Dysderidae).

The family Dysderidae is a highly diverse group of nocturnal ground-dwelling and active-hunter spiders. Dysderids are mostly restricted to the Western Palearctic, and particularly rich and abundant around the Mediterranean region. Interestingly, the distribution of species richness among its 24 genera and three subfamilies is highly biased-80% of its 644 documented species belong to just two genera, Dysdera (326) and Harpactea (211). Dysderidae provides an excellent study case for evolutionary and ecological research. It includes cases of trophic specialization, which are uncommon among spiders, and exhibit other remarkable biological (e.g. holocentric chromosomes), behavioural (e.g. cryptic female choice), evolutionary (e.g. adaptive radiation) and ecological features (e.g. recurrent colonization of the subterranean environment). The lack of a quantitative hypothesis on its phylogenetic structure has hampered its potential as a testing ground for evolutionary, biogeographical and ecological hypotheses. Here, we present the results of a target, multi-locus phylogenetic analysis, using mitochondrial (cox1, 16s and 12s) and nuclear genes (h3, 28s and 18s), of the most exhaustive taxonomic sample within Dysderidae (104 spp.) to date and across related families (Synspermiata) (83 spp.). We estimate divergence times using a combination of fossil and biogeographic node calibrations and use this timeline to identify shifts in diversification rates. Our results support the monophyly of the Dysderidae subfamilies Rhodinae and Dysderinae but reject Harpacteinae as currently defined. Moreover, the clades recovered within Harpacteinae do not support its current taxonomy. The origin of the family most likely post-dated the break-up of Pangea, and cave colonization may be older than previously considered. After correcting for the taxonomic artefacts, we identified a significant shift in diversification rates at the base of the genus Dysdera. Although the unique coexistence of specialist and generalist diets within the lineage could be suggested as the potential driver for the rate acceleration, further quantitative analyses would be necessary to test this hypothesis.

Monophyly, Taxon Sampling, and the Nature of Ranks in the Classification of Orb-Weaving Spiders (Araneae: Araneoidea).

We address some of the taxonomic and classification changes proposed by Kuntner et al. (2019) in a comparative study on the evolution of sexual size dimorphism in nephiline spiders. Their proposal to recircumscribe araneids and to rank the subfamily Nephilinae as a family is fundamentally flawed as it renders the family Araneidae paraphyletic. We discuss the importance of monophyly, outgroup selection, and taxon sampling, the subjectivity of ranks, and the implications of the age of origin criterion to assign categorical ranks in biological classifications. We explore the outcome of applying the approach of Kuntner et al. (2019) to the classification of spiders with emphasis on the ecribellate orb-weavers (Araneoidea) using a recently published dated phylogeny. We discuss the implications of including the putative sister group of Nephilinae (the sexually dimorphic genus Paraplectanoides) and the putative sister group of Araneidae (the miniature, monomorphic family Theridiosomatidae). We propose continuation of the phylogenetic classification put forth by Dimitrov et al. (2017), and we formally rank Nephilinae and Phonognathinae as subfamilies of Araneidae. Our classification better reflects the understanding of the phylogenetic placement and evolutionary history of nephilines and phonognathines while maintaining the diagnosability of Nephilinae. It also fulfills the fundamental requirement that taxa must be monophyletic, and thus avoids the paraphyly of Araneidae implied by Kuntner et al. (2019).

Repeated colonization, adaptive radiation and convergent evolution in the sheet-weaving spiders (Linyphiidae) of the south Pacific Archipelago of Juan Fernandez.

We report on the colonization and diversification of linyphiid spiders in the Pacific oceanic archipelago of Juan Fernandez. About 50 spider species occur naturally in these islands, most of them endemic and about half of them are linyphiids. Linyphiidae includes no fewer than 15 species of Laminacauda and three of Neomaso (with several additional undescribed species in the latter genus), all of them single island endemics. There are three additional linyphiid endemic genera, two monotypic and one, Juanfernandezia, with two species. Unlike the rather uniform somatic morphology and small ground sheet webs of the continental Laminacauda and Neomaso species, the Juan Fernandez endemics exhibit morphological features and life history traits that are very rare or unknown in any other linyphiids. A multi-locus phylogenetic analysis confirms at least five independent Juan Fernandez colonizations of Linyphiidae, two within the same genus, and three of which underwent subsequent local diversification. Different calibrations suggest alternative colonization timelines, some at odds with island ages, but all agree on similar diversification timings of the endemic lineages. Rare phenotypic traits (e.g. gigantism, massive chelicerae or elongated legs) evolved multiple times independently within the islands. Based on the remarkable levels of eco-phenotypic differentiation in locally diversified species showing densely packed distributions, we propose that Laminacauda, and probably Neomaso, constitute a case of adaptive radiation.

Converging on the orb: denser taxon sampling elucidates spider phylogeny and new analytical methods support repeated evolution of the orb web.

High throughput sequencing and phylogenomic analyses focusing on relationships among spiders have both reinforced and upturned long-standing hypotheses. Likewise, the evolution of spider webs-perhaps their most emblematic attribute-is being understood in new ways. With a matrix including 272 spider species and close arachnid relatives, we analyze and evaluate the relationships among these lineages using a variety of orthology assessment methods, occupancy thresholds, tree inference methods and support metrics. Our analyses include families not previously sampled in transcriptomic analyses, such as Symphytognathidae, the only araneoid family absent in such prior works. We find support for the major established spider lineages, including Mygalomorphae, Araneomorphae, Synspermiata, Palpimanoidea, Araneoidea and the Retrolateral Tibial Apophysis Clade, as well as the uloborids, deinopids, oecobiids and hersiliids Grade. Resulting trees are evaluated using bootstrapping, Shimodaira-Hasegawa approximate likelihood ratio test, local posterior probabilities and concordance factors. Using structured Markov models to assess the evolution of spider webs while accounting for hierarchically nested traits, we find multiple convergent occurrences of the orb web across the spider tree-of-life. Overall, we provide the most comprehensive spider tree-of-life to date using transcriptomic data and use new methods to explore controversial issues of web evolution, including the origins and multiple losses of the orb web.

The Grass was Greener: Repeated Evolution of Specialized Morphologies and Habitat Shifts in Ghost Spiders Following Grassland Expansion in South America.

While grasslands, one of Earth's major biomes, are known for their close evolutionary ties with ungulate grazers, these habitats are also paramount to the origins and diversification of other animals. Within the primarily South American spider subfamily Amaurobioidinae (Anyphaenidae), several species are found living in the continent's grasslands, with some displaying putative morphological adaptations to dwelling unnoticed in the grass blades. Herein, a dated molecular phylogeny provides the backbone for analyses revealing the ecological and morphological processes behind these spiders' grassland adaptations. The multiple switches from Patagonian forests to open habitats coincide with the expansion of South America's grasslands during the Miocene, while the specialized morphology of several grass-dwelling spiders originated at least three independent times and is best described as the result of different selective regimes operating on macroevolutionary timescales. Although grass-adapted lineages evolved towards different peaks in adaptive landscape, they all share one characteristic: an anterior narrowing of the prosoma allowing spiders to extend the first two pairs of legs, thus maintaining a slender resting posture in the grass blade. By combining phylogenetic, morphological, and biogeographic perspectives we disentangle multiple factors determining the evolution of a clade of terrestrial invertebrate predators alongside their biomes.

Arm-less mitochondrial tRNAs conserved for over 30 millions of years in spiders.

BACKGROUND: In recent years, Next Generation Sequencing (NGS) has accelerated the generation of full mitogenomes, providing abundant material for studying different aspects of molecular evolution. Some mitogenomes have been observed to harbor atypical sequences with bizarre secondary structures, which origins and significance could only be fully understood in an evolutionary framework. RESULTS: Here we report and analyze the mitochondrial sequences and gene arrangements of six closely related spiders in the sister genera Parachtes and Harpactocrates, which belong to the nocturnal, ground dwelling family Dysderidae. Species of both genera have compacted mitogenomes with many overlapping genes and strikingly reduced tRNAs that are among the shortest described within metazoans. Thanks to the conservation of the gene order and the nucleotide identity across close relatives, we were able to predict the secondary structures even on arm-less tRNAs, which would be otherwise unattainable for a single species. They exhibit aberrant secondary structures with the lack of either DHU or TΨC arms and many miss-pairings in the acceptor arm but this degeneracy trend goes even further since at least four tRNAs are arm-less in the six spider species studied. CONCLUSIONS: The conservation of at least four arm-less tRNA genes in two sister spider genera for about 30 myr suggest that these genes are still encoding fully functional tRNAs though they may be post-transcriptionally edited to be fully functional as previously described in other species. We suggest that the presence of overlapping and truncated tRNA genes may be related and explains why spider mitogenomes are smaller than those of other invertebrates.

Chance and predictability in evolution: The genomic basis of convergent dietary specializations in an adaptive radiation.

The coexistence of multiple eco-phenotypes in independently assembled communities makes island adaptive radiations the ideal framework to test convergence and parallelism in evolution. In the radiation of the spider genus Dysdera in the Canary Islands, species diversification occurs concomitant with repeated events of trophic specialization. These dietary shifts, to feed primarily on woodlice, are accompanied by modifications in morphology (mostly in the mouthparts), behaviour and nutritional physiology. To gain insight into the molecular basis of this adaptive radiation, we performed a comprehensive comparative transcriptome analysis of five Canary Island Dysdera endemics representing two evolutionary and geographically independent events of dietary specialization. After controlling for the potential confounding effects of hemiplasy, our differential gene expression and selective constraint analyses identified a number of genetic changes that could be associated with the repeated adaptations to specialized diet of woodlice, including some related to heavy metal detoxification and homeostasis, the metabolism of some important nutrients and venom toxins. Our results shed light on the genomic basis of an extraordinary case of dietary shift convergence associated with species diversification. We uncovered putative molecular substrates of convergent evolutionary changes at different hierarchical levels, including specific genes, genes with equivalent functions and even particular amino acid positions. This study improves our knowledge of rapid adaptive radiations and provides new insights into the predictability of evolution.

DOMINO: development of informative molecular markers for phylogenetic and genome-wide population genetic studies in non-model organisms.

MOTIVATION: The development of molecular markers is one of the most important challenges in phylogenetic and genome wide population genetics studies, especially in studies with non-model organisms. A highly promising approach for obtaining suitable markers is the utilization of genomic partitioning strategies for the simultaneous discovery and genotyping of a large number of markers. Unfortunately, not all markers obtained from these strategies provide enough information for solving multiple evolutionary questions at a reasonable taxonomic resolution. RESULTS: We have developed Development Of Molecular markers In Non-model Organisms (DOMINO), a bioinformatics tool for informative marker development from both next generation sequencing (NGS) data and pre-computed sequence alignments. The application implements popular NGS tools with new utilities in a highly versatile pipeline specifically designed to discover or select personalized markers at different levels of taxonomic resolution. These markers can be directly used to study the taxa surveyed for their design, utilized for further downstream PCR amplification in a broader set taxonomic scope, or exploited as suitable templates to bait design for target DNA enrichment techniques. We conducted an exhaustive evaluation of the performance of DOMINO via computer simulations and illustrate its utility to find informative markers in an empirical dataset. AVAILABILITY AND IMPLEMENTATION: DOMINO is freely available from www.ub.edu/softevol/domino CONTACT: elsanchez@ub.edu or jrozas@ub.eduSupplementary information: Supplementary data are available at Bioinformatics online.

The life and adventures of an eight-legged castaway: Colonization and diversification of Philisca ghost spiders on Robinson Crusoe Island (Araneae, Anyphaenidae).

Oceanic archipelagoes, by their young origin and isolation, provide privileged settings to study the origin and diversification of species. Here, we study the anyphaenid spider genus Philisca, endemic to the Valdivian temperate rainforest, which includes species living both on the mainland as well as on the Robison Crusoe Island in the Juan Fernández archipelago. Anyphaenids, as many spiders, are potentially good colonizers due their ability for ballooning, an airborne dispersal mediated by strands of silk that are caught in the wind. We use a molecular approach to estimate both the phylogenetic relationships and the timeframe of species diversification of Philisca, with the aim to infer its evolutionary history. We further estimate the rates of speciation on both the insular and continental Philisca species and score the microhabitat used by each species and their sizes as a proxy to evaluate ecological niche diversification within the island. Most analyses support the monophyly of Philisca, with the exclusion of Philisca tripunctata. Our results reveal colonization from a single lineage that postdated the origin of the island, followed by rapid (∼2Ma) diversification. The ancestral microhabitat was most likely leaf-dwelling but we identify two independent microhabitat shifts. Our data provides evidence that Philisca has undergone an adaptive radiation on the Robison Crusoe Island.

Rounding up the usual suspects: a standard target-gene approach for resolving the interfamilial phylogenetic relationships of ecribellate orb-weaving spiders with a new family-rank classification (Araneae, Araneoidea).

We test the limits of the spider superfamily Araneoidea and reconstruct its interfamilial relationships using standard molecular markers. The taxon sample (363 terminals) comprises for the first time representatives of all araneoid families, including the first molecular data of the family Synaphridae. We use the resulting phylogenetic framework to study web evolution in araneoids. Araneoidea is monophyletic and sister to Nicodamoidea rank. n. Orbiculariae are not monophyletic and also include the RTA clade, Oecobiidae and Hersiliidae. Deinopoidea is paraphyletic with respect to a lineage that includes the RTA clade, Hersiliidae and Oecobiidae. The cribellate orb-weaving family Uloboridae is monophyletic and is sister group to a lineage that includes the RTA Clade, Hersiliidae and Oecobiidae. The monophyly of most Araneoidea families is well supported, with a few exceptions. Anapidae includes holarchaeids but the family remains diphyletic even if Holarchaea is considered an anapid. The orb-web is ancient, having evolved by the early Jurassic; a single origin of the orb with multiple "losses" is implied by our analyses. By the late Jurassic, the orb-web had already been transformed into different architectures, but the ancestors of the RTA clade probably built orb-webs. We also find further support for a single origin of the cribellum and multiple independent losses. The following taxonomic and nomenclatural changes are proposed: the cribellate and ecribellate nicodamids are grouped in the superfamily Nicodamoidea rank n. (Megadictynidae rank res. and Nicodamidae stat. n.). Araneoidea includes 17 families with the following changes: Araneidae is re-circumscribed to include nephilines, Nephilinae rank res., Arkyidae rank n., Physoglenidae rank n., Synotaxidae is limited to the genus Synotaxus, Pararchaeidae is a junior synonym of Malkaridae (syn. n.), Holarchaeidae of Anapidae (syn. n.) and Sinopimoidae of Linyphiidae (syn. n.).

The spider tree of life: phylogeny of Araneae based on target-gene analyses from an extensive taxon sampling.

We present a phylogenetic analysis of spiders using a dataset of 932 spider species, representing 115 families (only the family Synaphridae is unrepresented), 700 known genera, and additional representatives of 26 unidentified or undescribed genera. Eleven genera of the orders Amblypygi, Palpigradi, Schizomida and Uropygi are included as outgroups. The dataset includes six markers from the mitochondrial (12S, 16S, COI) and nuclear (histone H3, 18S, 28S) genomes, and was analysed by multiple methods, including constrained analyses using a highly supported backbone tree from transcriptomic data. We recover most of the higher-level structure of the spider tree with good support, including Mesothelae, Opisthothelae, Mygalomorphae and Araneomorphae. Several of our analyses recover Hypochilidae and Filistatidae as sister groups, as suggested by previous transcriptomic analyses. The Synspermiata are robustly supported, and the families Trogloraptoridae and Caponiidae are found as sister to the Dysderoidea. Our results support the Lost Tracheae clade, including Pholcidae, Tetrablemmidae, Diguetidae, Plectreuridae and the family Pacullidae (restored status) separate from Tetrablemmidae. The Scytodoidea include Ochyroceratidae along with Sicariidae, Scytodidae, Drymusidae and Periegopidae; our results are inconclusive about the separation of these last two families. We did not recover monophyletic Austrochiloidea and Leptonetidae, but our data suggest that both groups are more closely related to the Cylindrical Gland Spigot clade rather than to Synspermiata. Palpimanoidea is not recovered by our analyses, but also not strongly contradicted. We find support for Entelegynae and Oecobioidea (Oecobiidae plus Hersiliidae), and ambiguous placement of cribellate orb-weavers, compatible with their non-monophyly. Nicodamoidea (Nicodamidae plus Megadictynidae) and Araneoidea composition and relationships are consistent with recent analyses. We did not obtain resolution for the titanoecoids (Titanoecidae and Phyxelididae), but the Retrolateral Tibial Apophysis clade is well supported. Penestomidae, and probably Homalonychidae, are part of Zodarioidea, although the latter family was set apart by recent transcriptomic analyses. Our data support a large group that we call the marronoid clade (including the families Amaurobiidae, Desidae, Dictynidae, Hahniidae, Stiphidiidae, Agelenidae and Toxopidae). The circumscription of most marronoid families is redefined here. Amaurobiidae include the Amaurobiinae and provisionally Macrobuninae. We transfer Malenellinae (Malenella, from Anyphaenidae), Chummidae (Chumma) (new syn.) and Tasmarubriinae (Tasmarubrius, Tasmabrochus and Teeatta, from Amphinectidae) to Macrobuninae. Cybaeidae are redefined to include Calymmaria, Cryphoeca, Ethobuella and Willisius (transferred from Hahniidae), and Blabomma and Yorima (transferred from Dictynidae). Cycloctenidae are redefined to include Orepukia (transferred from Agelenidae) and Pakeha and Paravoca (transferred from Amaurobiidae). Desidae are redefined to include five subfamilies: Amphinectinae, with Amphinecta, Mamoea, Maniho, Paramamoea and Rangitata (transferred from Amphinectidae); Ischaleinae, with Bakala and Manjala (transferred from Amaurobiidae) and Ischalea (transferred from Stiphidiidae); Metaltellinae, with Austmusia, Buyina, Calacadia, Cunnawarra, Jalkaraburra, Keera, Magua, Metaltella, Penaoola and Quemusia; Porteriinae (new rank), with Baiami, Cambridgea, Corasoides and Nanocambridgea (transferred from Stiphidiidae); and Desinae, with Desis, and provisionally Poaka (transferred from Amaurobiidae) and Barahna (transferred from Stiphidiidae). Argyroneta is transferred from Cybaeidae to Dictynidae. Cicurina is transferred from Dictynidae to Hahniidae. The genera Neoramia (from Agelenidae) and Aorangia, Marplesia and Neolana (from Amphinectidae) are transferred to Stiphidiidae. The family Toxopidae (restored status) includes two subfamilies: Myroinae, with Gasparia, Gohia, Hulua, Neomyro, Myro, Ommatauxesis and Otagoa (transferred from Desidae); and Toxopinae, with Midgee and Jamara, formerly Midgeeinae, new syn. (transferred from Amaurobiidae) and Hapona, Laestrygones, Lamina, Toxops and Toxopsoides (transferred from Desidae). We obtain a monophyletic Oval Calamistrum clade and Dionycha; Sparassidae, however, are not dionychans, but probably the sister group of those two clades. The composition of the Oval Calamistrum clade is confirmed (including Zoropsidae, Udubidae, Ctenidae, Oxyopidae, Senoculidae, Pisauridae, Trechaleidae, Lycosidae, Psechridae and Thomisidae), affirming previous findings on the uncertain relationships of the "ctenids" Ancylometes and Cupiennius, although a core group of Ctenidae are well supported. Our data were ambiguous as to the monophyly of Oxyopidae. In Dionycha, we found a first split of core Prodidomidae, excluding the Australian Molycriinae, which fall distantly from core prodidomids, among gnaphosoids. The rest of the dionychans form two main groups, Dionycha part A and part B. The former includes much of the Oblique Median Tapetum clade (Trochanteriidae, Gnaphosidae, Gallieniellidae, Phrurolithidae, Trachelidae, Gnaphosidae, Ammoxenidae, Lamponidae and the Molycriinae), and also Anyphaenidae and Clubionidae. Orthobula is transferred from Phrurolithidae to Trachelidae. Our data did not allow for complete resolution for the gnaphosoid families. Dionycha part B includes the families Salticidae, Eutichuridae, Miturgidae, Philodromidae, Viridasiidae, Selenopidae, Corinnidae and Xenoctenidae (new fam., including Xenoctenus, Paravulsor and Odo, transferred from Miturgidae, as well as Incasoctenus from Ctenidae). We confirm the inclusion of Zora (formerly Zoridae) within Miturgidae.

Imprints of multiple glacial refugia in the Pyrenees revealed by phylogeography and palaeodistribution modelling of an endemic spider.

Mediterranean mountain ranges harbour highly endemic biota in islandlike habitats. Their topographic diversity offered the opportunity for mountain species to persist in refugial areas during episodes of major climatic change. We investigate the role of Quaternary climatic oscillations in shaping the demographic history and distribution ranges in the spider Harpactocrates ravastellus, endemic to the Pyrenees. Gene trees and multispecies coalescent analyses on mitochondrial and nuclear DNA sequences unveiled two distinct lineages with a hybrid zone around the northwestern area of the Catalan Pyrenees. The lineages were further supported by morphological differences. Climatic niche-based species distribution models (SDMs) identified two lowland refugia at the western and eastern extremes of the mountain range, which would suggest secondary contact following postglacial expansion of populations from both refugia. Neutrality test and approximate Bayesian computation (ABC) analyses indicated that several local populations underwent severe bottlenecks followed by population expansions, which in combination with the deep population differentiation provided evidence for population survival during glacial periods in microrefugia across the mountain range, in addition to the main Atlantic and Mediterranean (western and eastern) refugia. This study sheds light on the complexities of Quaternary climatic oscillations in building up genetic diversity and local endemicity in the southern Europe mountain ranges.

Extant primitively segmented spiders have recently diversified from an ancient lineage.

Living fossils are lineages that have retained plesiomorphic traits through long time periods. It is expected that such lineages have both originated and diversified long ago. Such expectations have recently been challenged in some textbook examples of living fossils, notably in extant cycads and coelacanths. Using a phylogenetic approach, we tested the patterns of the origin and diversification of liphistiid spiders, a clade of spiders considered to be living fossils due to their retention of arachnid plesiomorphies and their exclusive grouping in Mesothelae, an ancient clade sister to all modern spiders. Facilitated by original sampling throughout their Asian range, we here provide the phylogenetic framework necessary for reconstructing liphistiid biogeographic history. All phylogenetic analyses support the monophyly of Liphistiidae and of eight genera. As the fossil evidence supports a Carboniferous Euramerican origin of Mesothelae, our dating analyses postulate a long eastward over-land dispersal towards the Asian origin of Liphistiidae during the Palaeogene (39-58 Ma). Contrary to expectations, diversification within extant liphistiid genera is relatively recent, in the Neogene and Late Palaeogene (4-24 Ma). While no over-water dispersal events are needed to explain their evolutionary history, the history of liphistiid spiders has the potential to play prominently in vicariant biogeographic studies.

Winter is coming: Miocene and Quaternary climatic shifts shaped the diversification of Western-Mediterranean Harpactocrates (Araneae, Dysderidae) spiders.

Past climatic shifts have played a major role in generating and shaping biodiversity. Quaternary glacial cycles are the better known examples of dramatic climatic changes endured by ecosystems in temperate regions. Although still a matter of debate, some authors suggest that glaciations promoted speciation. Here we investigate the effect of past climatic changes on the diversification of the ground-dwelling spider genus Harpactocrates, distributed across the major mountain ranges of the western Mediterranean. Concatenated and species-tree analyses of multiple mitochondrial and nuclear loci, combined with the use of fossil and biogeographic calibration points, reveal a Miocene origin of most nominal species, but also unravel several cryptic lineages tracing back to the Pleistocene. We hypothesize that the Miocene Climatic Transition triggered major extinction events in the genus but also promoted its subsequent diversification. Under this scenario, the Iberian mountains acted as an island-like system, providing shelter to Harpactocrates lineages during the climate shifts and favouring isolation between mountain ranges. Quaternary glacial cycles contributed further to the diversification of the group by isolating lineages in peripheral refugia within mountain ranges. In addition, we recovered some unique biogeographic patterns, such as the colonization of the Alps and the Apennines from the Iberian Peninsula.

Taxonomic accuracy and complementarity between bulk and eDNA metabarcoding provides an alternative to morphology for biological assessment of freshwater macroinvertebrates.

Determining biological status of freshwater ecosystems is critical for ensuring ecosystem health and maintaining associated services to such ecosystems. Freshwater macroinvertebrates respond predictably to environmental disturbances and are widely used in biomonitoring programs. However, many freshwater species are difficult to capture and sort from debris or substrate and morphological identification is challenging, especially larval stages, damaged specimens, or hyperdiverse groups such as Diptera. The advent of high throughput sequencing technologies has enhanced DNA barcoding tools to automatise species identification for whole communities, as metabarcoding is increasingly used to monitor biodiversity. However, recent comparisons have revealed little congruence between morphological and molecular-based identifications. Using broad range universal primers for DNA barcode marker cox1, we compare community composition captured between morphological and molecular-based approaches from different sources - tissue-based (bulk benthic and bulk drift samples) and environmental DNA (eDNA, filtered water) metabarcoding - for samples collected along a gradient of anthropogenic disturbances. For comparability, metabarcoding taxonomic assignments were filtered by taxa included in the standardised national biological metric IBMWP. At the family level, bulk benthic metabarcoding showed the highest congruence with morphology, and the most abundant taxa were captured by all techniques. Richness captured by morphology and bulk benthic metabarcoding decreased along the gradient, whereas richness recorded by eDNA remained constant and increased downstream when sequencing bulk drift. Estimates of biological metrics were higher using molecular than morphological identification. At species level, diversity captured by bulk benthic samples were higher than the other techniques. Importantly, bulk benthic and eDNA metabarcoding captured different and complementary portions of the community - benthic versus water column, respectively - and their combined use is recommended. While bulk benthic metabarcoding can likely replace morphology using similar benthic biological indices, water eDNA will require new metrics because this technique sequences a different portion of the community.

Optimal inventorying and monitoring of taxonomic, phylogenetic and functional diversity.

Comparable data is essential to understand biodiversity patterns. While assemblage or community inventorying requires comprehensive sampling, monitoring focuses on as few components as possible to detect changes. Quantifying species, their evolutionary history, and the way they interact requires studying changes in taxonomic (TD), phylogenetic (PD) and functional diversity (FD). Here we propose a method for the optimization of sampling protocols for inventorying and monitoring assemblages or communities across these three diversity dimensions taking sampling costs into account. We used Iberian spiders and Amazonian bats as two case-studies. The optimal combination of methods for inventorying and monitoring required optimizing the accumulation curve of α-diversity and minimizing the difference between sampled and estimated ß-diversity (bias), respectively. For Iberian spiders, the optimal combination for TD, PD and FD allowed sampling at least 50% of estimated diversity with 24 person-hours of fieldwork. The optimal combination of six person-hours allowed reaching a bias below 8% for all dimensions. For Amazonian bats, surveying all the 12 sites with mist-nets and 0 or 1 acoustic recorders was the optimal combination for almost all diversity types, resulting in >89% of the diversity and <10% bias with roughly a third of the cost. Only for phylogenetic α-diversity, the best solution was less clear and involved surveying both with mist nets and acoustic recorders. The widespread use of optimized and standardized sampling protocols and regular repetition in time will radically improve global inventory and monitoring of biodiversity. We strongly advocate for the global adoption of sampling protocols for both inventory and monitoring of taxonomic, phylogenetic and functional diversity.

From biomarkers to community composition: Negative effects of UV/chlorine-treated reclaimed urban wastewater on freshwater biota.

The use of urban wastewater reclaimed water has recently increased across the globe to restore stream environmental flows and mitigate the effects of water scarcity. Reclaimed water is disinfected using different treatments, but their effects into the receiving rivers are little studied. Physiological bioassays and biomarkers can detect sub-lethal effects on target species, but do not provide information on changes in community structure. In contrast, official monitoring programs use community structure information but often at coarse taxonomic resolution level that may fail to detect species level impacts. Here, we combined commonly used biomonitoring approaches from organism physiology to community species composition to scan a broad range of effects of disinfection of reclaimed water by UV-light only and both UV/chlorine on the biota. We (1) performed bioassays in one laboratory species (water flea Daphnia magna) and measured biomarkers in two wild species (caddisfly Hydropsyche exocellata and the barbel Luciobarbus graellsii), (2) calculated standard indices of biotic quality (IBQ) for diatoms, benthic macroinvertebrates, and fishes, and (3) analysed community species composition of eukaryotes determined by Cytochrome Oxidase C subunit I (cox1) metabarcoding. Only the UV/chlorine treatment caused significant changes in feeding rates of D. magna and reduced antioxidant defenses, increased anaerobic metabolism and altered the levels of lipid peroxidiation in H. exocellata. However, inputs of reclaimed water were significantly associated with a greater proportion of circulating neutrophils and LG-PAS cells in L. graellsii. Despite IBQ did not discriminate between the two water treatments, metabarcoding data detected community composition changes upon exposure to UV/chlorine reclaimed water. Overall, despite the effects of UV/chlorine-treated water were transient, our study suggests that UV-light treated is less harmful for freshwater biota than UV/chlorine-treated reclaimed water, but those effects depend of the organizational level.