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.

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.

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.

The Vitruvian spider: Segmenting and integrating over different body parts to describe ecophenotypic variation.

Understanding what drives the existing phenotypic variability has been a major topic of interest for biologists for generations. However, the study of the phenotype may not be straightforward. Indeed, organisms may be interpreted as composite objects, comprising different ecophenotypic traits, which are neither necessarily independent from each other nor do they respond to the same evolutionary pressures. For this reason, a deep biological understanding of the focal organism is essential for any morphological analysis. The spider genus Dysdera provides a particularly well-suited system for setting up protocols for morphological analyses that encompass a suit of morphological structures in any nonmodel system. This genus has undergone a remarkable diversification in the Canary Islands, where different species perform different ecological roles, exhibiting different levels of trophic specialization or troglomorphic adaptations, which translate into a remarkable interspecific morphological variability. Here, we seek to develop a broad guide, of which morphological characters must be considered, to study the effect of different ecological pressures in spiders and propose a general workflow that will be useful whenever researchers set out to investigate variation in the body plans of different organisms, with data sets comprising a set of morphological traits. We use geometric morphometric methods to quantify variation in different body structures, all of them with diverse phenotypic modifications in their chelicera, prosoma, and legs. We explore the effect of analyzing different combined landmark (LM) configurations of these characters and the degree of morphological integration that they exhibit. Our results suggest that different LM configurations of each of these body parts exhibit a higher degree of integration compared to LM configurations from different structures and that the analysis of each of these body parts captures different aspects of morphological variation, potentially related to different ecological factors.

For all audiences: Incorporating immature stages into standardised spider inventories has a major impact on the assessment of biodiversity patterns.

Although arthropods are the largest component of animal diversity, they are traditionally underrepresented in biological inventories and monitoring programmes. However, no biodiversity assessment can be considered informative without including them. Arthropod immature stages are often discarded during sorting, despite frequently representing more than half of the collected individuals. To date, little effort has been devoted to characterising the impact of discarding nonadult specimens on our diversity estimates. Here, we used a metabarcoding approach to analyse spiders from oak forests in the Iberian Peninsula, to assess (1) the contribution of juvenile stages to local diversity estimates, and (2) their effect on the diversity patterns (compositional differences) across assemblages. We further investigated the ability of metabarcoding to inform on abundance. We obtained 363 and 331 species as adults and juveniles, respectively. Including the species represented only by juveniles increased the species richness of the whole sampling in 35% with respect to those identified from adults. Differences in composition between assemblages were greatly reduced when immature stages were considered, especially across latitudes, possibly due to phenological differences. Moreover, our results revealed that metabarcoding data are to a certain extent quantitative, but some sort of taxonomic conversion factor may be necessary to provide accurate informative estimates. Although our findings do not question the relevance of the information provided by adult-based inventories, they also reveal that juveniles provide a novel and relevant layer of knowledge that, especially in areas with marked seasonality, may influence our interpretations, providing more accurate information from standardised biological inventories.

The chromosome-scale assembly of the Canary Islands endemic spider Dysdera silvatica (Arachnida, Araneae) sheds light on the origin and genome structure of chemoreceptor gene families in chelicerates.

Here, we present the chromosome-level genome assembly of Dysdera silvatica Schmidt, 1981, a nocturnal ground-dwelling spider endemic from the Canary Islands. The genus Dysdera has undergone a remarkable diversification in this archipelago mostly associated with shifts in the level of trophic specialization, becoming an excellent model to study the genomic drivers of adaptive radiations. The new assembly (1.37 Gb; scaffold N50 of 174.2 Mb), was performed using the chromosome conformation capture scaffolding technique, represents a continuity improvement of more than 4500 times with respect to the previous version. The seven largest scaffolds or pseudochromosomes, which cover 87% of the total assembly size, probably correspond with the seven chromosomes of the karyotype of this species, including a characteristic large X chromosome. To illustrate the value of this new resource we performed a comprehensive analysis of the two major arthropod chemoreceptor gene families (i.e., gustatory and ionotropic receptors). We identified 545 chemoreceptor sequences distributed across all pseudochromosomes, with a notable underrepresentation in the X chromosome. At least 54% of them localize in 83 genomic clusters with a significantly lower evolutionary distances between them than the average of the family, suggesting a recent origin of many of them. This chromosome-level assembly is the first high-quality genome representative of the Synspermiata clade, and just the third among spiders, representing a new valuable resource to gain insights into the structure and organization of chelicerate genomes, including the role that structural variants, repetitive elements and large gene families played in the extraordinary biology of spiders.

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.

Twenty years, eight legs, one concept: describing spider biodiversity in Zootaxa (Arachnida: Araneae).

Zootaxa published more than a thousand papers on Araneae from 2002 to the present, including descriptions of 3,833 new spider species and 177 new genera. Here we summarise the key contributions of Zootaxa to our current knowledge of global spider diversity. We provide a historical account of the researchers that have actively participated as editors, and recognize the more than 1,000 reviewers without whom none of this would have been possible. We conduct a simple analysis of the contributions by authors and geographic region, which allows us to uncover some of the underlying trends in current spider taxonomy. In addition, we examine some of the milestones in twenty years of spider systematic research in Zootaxa. Finally, we discuss future prospects of spider taxonomy and the role that Zootaxa and its younger sister journal Megataxa will play in it. We would like to dedicate this contribution to the memory of Norman I. Platnick, a crucial figure in the advancement of spider systematics.

A new species of Cerroneuroterus Melika Pujade-Villar from the Eastern Palaearctic (Hymenoptera, Cynipidae, Cynipini).

A new gall wasp species, Cerroneuroterus yukawamasudai Pujade-Villar Melika sp. nov., is described from China. The asexual generation induces leaf lenticular galls, and the sexual generation induces catkin galls on Q. acutissima and Q. variabilis. Sexual females and males are not described yet. Data on the diagnosis, distribution, and biology of the new species are provided, including a key to sexual and asexual females of all described Cerroneuroterus species. Previous misidentifications of host galls recorded in Japan are commented on. Molecular analyses were performed to assess the validity of the genus Cerroneuroterus and the affiliation of the new species.

Current status of Andricus mairei (Kieffer), with synonymization of two species from China (Hymenoptera: Cynipidae).

The original description of Parandricus mairei Kieffer, 1906 included a misinterpretation of a relevant forewing trait. The species was subsequently transferred to the genus Andricus, despite presenting two very peculiar morphological characters, namely a simple tarsal claw and the absence of irradiating carinae on the lower face. Similarly, the original descriptions of Andricus deqingis Wang, Gui, Chen, 2013 and A. flavus Pujade-Villar, Wang, Guo Chen, 2014 included some relevant mistakes. Here, we present the results of a molecular analysis that reveals that individuals of the three species are genetically very similar; supporting the proposal that A. mairei is a senior synonym of the other two species, A. deqingis n. syn. and A. flavus n. syn. In addition, our results indicate that Parandricus renders Andricus paraphyletic, which supports that Parandricus is a junior synonym of Andricus. We re-describe and illustrate the relevant characters of A. mairei and provide additional comments on the characters erroneously interpreted in the former descriptions and give notes about the biology and intraspecific variability.

A New Genus of Oak Gall Wasp, Striatoandricus Pujade-Villar (Hymenoptera: Cynipidae: Cynipini) from America with Descriptions of Two New Mexican Species.

A new genus of cynipid oak gall wasp, Striatoandricus Pujade-Villar (Hymenoptera: Cynipidae: Cynipini), is described. Striatoandricus gen. nov. includes four previously described species, Andricus nievesaldreyi n. comb., A. georgei n. comb., A. maesi n. comb., and A. barriosi n. comb., which induce pubescent leaves or twig galls on Quercus belonging to Quercus section. Two new species from México are also described: S. cuixarti Pujade-Villar n. sp. and S. sanchezi Pujade-Villar n. sp. in Quercus section. Descriptions of the genus and diagnostic characters, including DNA sequence data, are presented. This new genus is supported by both morphological and molecular data.

A database of functional traits for spiders from native forests of the Iberian Peninsula and Macaronesia.

BACKGROUND: There is an increasing demand for databases including species trait information for biodiversity and community ecology studies. The existence of trait databases is useful for comparative studies within taxa or geographical regions, but there is low availability of databases for certain organisms. Here we present an open access functional trait database for spiders from Macaronesia and the Iberian Peninsula, recording several morphological and ecological traits related to the species life histories, microhabitat and trophic preferences. NEW INFORMATION: We present a database that includes 12 biological traits for 506 spider species present in natural forests of the Iberian Peninsula (Spain) and three Macaronesian archipelagoes (Azores, Madeira and Canary Islands). The functional trait database consists of two sections:individual-level data for six morphological traits (total body size, prosoma length, prosoma width, prosoma height, tibia I length and fang length), based on direct measurements of 2844 specimens of all spider species; andspecies-level aggregate data for 12 traits (same 6 morphological traits as in the previous section plus dispersal ability, vertical stratification, circadian activity, foraging strategy, trophic specialization and colonization status), based on either the average of the direct measurements or bibliographic searches.This functional trait database will serve as a data standard for currently ongoing analyses that require trait and functional diversity statistics.

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

The draft genome sequence of the spider Dysdera silvatica (Araneae, Dysderidae): A valuable resource for functional and evolutionary genomic studies in chelicerates.

BACKGROUND: We present the draft genome sequence of Dysdera silvatica, a nocturnal ground-dwelling spider from a genus that has undergone a remarkable adaptive radiation in the Canary Islands. RESULTS: The draft assembly was obtained using short (Illumina) and long (PaciBio and Nanopore) sequencing reads. Our de novo assembly (1.36 Gb), which represents 80% of the genome size estimated by flow cytometry (1.7 Gb), is constituted by a high fraction of interspersed repetitive elements (53.8%). The assembly completeness, using BUSCO and core eukaryotic genes, ranges from 90% to 96%. Functional annotations based on both ab initio and evidence-based information (including D. silvatica RNA sequencing) yielded a total of 48,619 protein-coding sequences, of which 36,398 (74.9%) have the molecular hallmark of known protein domains, or sequence similarity with Swiss-Prot sequences. The D. silvatica assembly is the first representative of the superfamily Dysderoidea, and just the second available genome of Synspermiata, one of the major evolutionary lineages of the "true spiders" (Araneomorphae). CONCLUSIONS: Dysderoids, which are known for their numerous instances of adaptation to underground environments, include some of the few examples of trophic specialization within spiders and are excellent models for the study of cryptic female choice. This resource will be therefore useful as a starting point to study fundamental evolutionary and functional questions, including the molecular bases of the adaptation to extreme environments and ecological shifts, as well of the origin and evolution of relevant spider traits, such as the venom and silk.

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.

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.

A DNA barcode-assisted annotated checklist of the spider (Arachnida, Araneae) communities associated to white oak woodlands in Spanish National Parks.

BACKGROUND: A large scale semi-quantitative biodiversity assessment was conducted in white oak woodlands in areas included in the Spanish Network of National Parks, as part of a project aimed at revealing biogeographic patterns and identify biodiversity drivers. The semi-quantitative COBRA sampling protocol was conducted in sixteen 1-ha plots across six national parks using a nested design. All adult specimens were identified to species level based on morphology. Uncertain delimitations and identifications due to either limited information of diagnostic characters or conflicting taxonomy were further investigated using DNA barcode information. NEW INFORMATION: We identified 376 species belonging to 190 genera in 39 families, from the 8,521 adults found amongst the 20,539 collected specimens. Faunistic results include the discovery of 7 new species to the Iberian Peninsula, 3 new species to Spain and 11 putative new species to science. As largely expected by environmental features, the southern parks showed a higher proportion of Iberian and Mediterranean species than the northern parks, where the Palearctic elements were largely dominant. The analysis of approximately 3,200 DNA barcodes generated in the present study, corroborated and provided finer resolution to the morphologically based delimitation and identification of specimens in some taxonomically challenging families. Specifically, molecular data confirmed putative new species with diagnosable morphology, identified overlooked lineages that may constitute new species, confirmed assignment of specimens of unknown sexes to species and identified cases of misidentifications and phenotypic polymorphisms.