Phylogenomics, Diversification Dynamics, and Comparative Transcriptomics across the Spider Tree of Life

Mystery on the Bounty: The family-level status of Pacificana cockayni Hogg, 1904 (Araneae)

The family-level placement of the species Pacificana cockayniHogg, 1904 (Araneae, Miturgidae) has been ambiguous for over a century, with the monotypic genus Pacificana initially placed in Agelenidae, later transferred to Amaurobioidinae (Anyphaenidae), and presently in Miturgidae. A recent work describing the male and molecular data consisting of a single mitochondrial gene, cytochrome c oxidase subunit I, confirmed that the species is part of the marronoid clade; however, these data did not result in a conclusive family-level placement. Here, we use low-coverage whole genome sequencing (lcWGS) combined with data from the Sequence Read Archive to infer a phylogeny from ultraconserved elements (UCEs) and six legacy Sanger loci. Indications of potential family placements from prior work and the topologies from this study support a transfer of Pacificana Hogg, 1904 to Cycloctenidae Simon, 1898 (new family placement).

Comparative mitogenomics of Cheiracanthium species (Araneae: Cheiracanthiidae) with phylogenetic implication and evolutionary insights

The genus Cheiracanthium C. L. Koch, 1839 is the most species-rich genus of the family Cheiracanthiidae. Given the unavailability of information on the evolutionary biology and molecular taxonomy of this genus, here we sequenced nine mitochondrial genomes (mitogenomes) of Cheiracanthium species, four of which were fully annotated, and conducted comparative analyses with other well-characterized Araneae mitogenomes. We also provide phylogenetic insights on the genus Cheiracanthium. The circular mitogenomes of the Cheiracanthium contain 37 genes, including 13 protein-coding genes (PCGs), 22 transfer RNA genes (tRNAs), two ribosomal RNA genes (rRNAs) and one putative control region (CR). All genes show a high A+T bias, characterized by a negative AT skew and positive GC skew, along with numerous overlapped regions and intergenic spacers. Approximately half of the tRNAs lack TΨC and/or dihydrouracil (DHU) arm and are characterized with unpaired amino acid acceptor arms. Most PCGs used the standard ATN start codons and TAR termination codons. The mitochondrial gene order of Cheiracanthium differs significantly from the putative ancestral gene order (Limulus polyphemus). Our novel phylogenetic analyses infer Cheiracanthiidae to be the sister group of Salticidae in BI analysis, but as sister to the node with Miturgidae, Viridasiidae, Corinnidae, Selenopidae, Salticidae, and Philodromidae in ML analysis. We confirm that Cheiracanthium is paraphyletic, for the first time using molecular phylogenetic approaches, with the earliest divergence estimated at 67 Ma. Our findings enhance our understanding of Cheiracanthium taxonomy and evolution.

Reassessing the evolutionary relationships of tropical wandering spiders using phylogenomics: A UCE-based phylogeny of Ctenidae (Araneae) with the discovery of a new lycosoid family

Tropical wandering spiders (Ctenidae) are a diverse family of cursorial predators whose species richness peaks in the tropics. The phylogeny of Ctenidae has been examined using morphology and Sanger-based sequencing data, but these studies have been limited by taxon sampling and have often recovered low branch support for many intrafamilial phylogenetic relationships. Herein, we present the most extensive phylogenetic sampling of this family using genome-scale data, leveraging museum collections of all ctenid subfamilies from across the world. We obtained a well-resolved phylogeny of Ctenidae, with the majority of nodes showing maximal nodal support and topological congruence across different phylogenetic analyses. For the first time, we show with high support that Ancylometes is not within Ctenidae but is the sister lineage to all the remaining lycosoid families. Therefore, we propose Ancylometidae as a new family. We assess the phylogenetic position of Ctenidae within Lycosoidea using a variety of phylogenetic methods and tests, demonstrating that the previously proposed position of Ctenidae as the sister clade of Psechridae, based on phylotranscriptomic analyses, lacks phylogenetic support. As a new finding, this study shows that the subfamily Acantheinae, as currently delimited, is polyphyletic. Therefore, we erect the new ctenid subfamily Enoplocteninae to accommodate the Neotropical genera Enoploctenus, Chococtenus, and Phymatoctenus. Our phylogenomic results using UCE data resolve the position of several problematic genera (e.g., Califorctenus and Acantheis) and add support to other parts of the tree that received low support in the most recent Sanger-based phylogeny. We discuss some of the putative morphological synapomorphies of the main ctenid lineages within the phylogenetic framework provided by the molecular phylogenetic results of this study.

Molecular phylogenetics of nursery web spiders (Araneae: Pisauridae)

Pisauridae Simon, 1890, or nursery web spiders, are a large family with a worldwide distribution and very diverse life history strategies. Despite being named for their nursery webs, similar structures are built by some members of Ctenidae, Trechaleidae, and Oxyopidae. Pisauridae has no known morphological synapomorphies that circumscribe all members of the family, and delineation of subfamilies has been a longstanding issue. In addition, several recent molecular phylogenetic studies have called into question the monophyly of Pisauridae. Here, we infer a phylogeny of Pisauridae with nine genetic markers (12S, 16S, 18S, 28S, actin, COI, histone H3, ITS2, and NADH) combined with ultraconserved elements (UCEs) to test the monophyly of the group and determine intrafamilial relationships. Our study includes a total of 81 terminals (59 of them pisaurids). Our maximum likelihood and gene coalescence analyses strongly suggest that as currently circumscribed, the family Pisauridae is not monophyletic. To circumscribe Pisauridae as a monophyletic group, based on our molecular results, we propose the restoration of the family Dolomedidae Simon, 1876 (rank resurrected) to include the genera DolomedesLatreille, 1804 and BradystichusSimon, 1884. In addition, based on morphological evidence, we also place in Dolomedidae the genera MegadolomedesDavies & Raven, 1980; Tasmomedes Raven, 2018; Mangromedes Raven, 2018; OrnodolomedesRaven & Hebron, 2018; and CaledomedesRaven & Hebron, 2018. We provide a phylogenetic delimitation of the three subfamilies of Pisauridae: Halinae, Thaumasiinae, and Pisaurinae. Ten pisaurid genera are deemed incertae sedis in terms of their placement. In light of our hypothesized molecular phylogeny, we discuss the morphological characters and putative synapomorphies of Pisauridae and propose diagnostic characters for its subfamilies and for the family Dolomedidae, along with taxonomic notes about genera not included in our study.

Molecular Phylogenetic Relationships Based on Mitogenomes of Spider: Insights Into Evolution and Adaptation to Extreme Environments

In this study, we performed a comparative analysis based on a total of 255 spider mitogenomes and four outgroups, of which the mitogenomes of 39 species were assembled de novo, to explore the phylogenetic relationships and the adaptive evolution of mitogenomes. Results showed that Argyroneta aquatica had the longest mitochondrial length and the most pronounced codon preference to be UUA, followed by CCU. Codon usage frequencies were similar between families and codon usage in the mitogenome of spiders was mainly influenced by natural selection pressures rather than G/C mutation bias. Our phylogenetic topology clearly explained the evolutionary relationships among the spiders, and divergence time estimates indicated that the spiders originated in the early Devonian, and that the two clades of Mesothelae and Opisthothelae separated in the late Carboniferous. Ancestral range and trait reconstruction results supported the ancestral origin of spiders to the Devonian Nearctic realm, with the trapdoor being the original trait. Selection analysis detected positive selection signals in the ATP8 gene in Desis jiaxiangi. The ND5 gene is a convergent evolutionary gene between D. jiaxiangi and A. aquatica. Positive selection signals in the ATP8 gene and convergent selection sites in the ND5 gene may facilitate metabolic adaptation to the aquatic environment in two aquatic spiders. In conclusion, our analysis contributes to a better understanding of the taxonomic status, species diversity, mitochondrial characteristics, and environmental adaptations of these spiders.

Three Novel Spider Genomes Unveil Spidroin Diversification and Hox Cluster Architecture: Ryuthela nishihirai (Liphistiidae), Uloborus plumipes (Uloboridae) and Cheiracanthium punctorium (Cheiracanthiidae)

Spiders are a hyperdiverse taxon and among the most abundant predators in nearly all terrestrial habitats. Their success is often attributed to key developments in their evolution such as silk and venom production and major apomorphies such as a whole-genome duplication. Resolving deep relationships within the spider tree of life has been historically challenging, making it difficult to measure the relative importance of these novelties for spider evolution. Whole-genome data offer an essential resource in these efforts, but also for functional genomic studies. Here, we present de novo assemblies for three spider species: Ryuthela nishihirai (Liphistiidae), a representative of the ancient Mesothelae, the suborder that is sister to all other extant spiders; Uloborus plumipes (Uloboridae), a cribellate orbweaver whose phylogenetic placement is especially challenging; and Cheiracanthium punctorium (Cheiracanthiidae), which represents only the second family to be sequenced in the hyperdiverse Dionycha clade. These genomes fill critical gaps in the spider tree of life. Using these novel genomes along with 25 previously published ones, we examine the evolutionary history of spidroin gene and structural hox cluster diversity. Our assemblies provide critical genomic resources to facilitate deeper investigations into spider evolution. The near chromosome-level genome of the 'living fossil' R. nishihirai represents an especially important step forward, offering new insights into the origins of spider traits.

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.

Tracing the evolutionary and genetic footprints of atmospheric tillandsioids transition from land to air

Plant evolution is driven by key innovations of functional traits that enables their survivals in diverse ecological environments. However, plant adaptive evolution from land to atmospheric niches remains poorly understood. In this study, we use the epiphytic Tillandsioideae subfamily of Bromeliaceae as model plants to explore their origin, evolution and diversification. We provide a comprehensive phylogenetic tree based on nuclear transcriptomic sequences, indicating that core tillandsioids originated approximately 11.3 million years ago in the Andes. The geological uplift of the Andes drives the divergence of tillandsioids into tank-forming and atmospheric types. Our genomic and transcriptomic analyses reveal gene variations and losses associated with adaptive traits such as impounding tanks and absorptive trichomes. Furthermore, we uncover specific nitrogen-fixing bacterial communities in the phyllosphere of tillandsioids as potential source of nitrogen acquisition. Collectively, our study provides integrative multi-omics insights into the adaptive evolution of tillandsioids in response to elevated aerial habitats.

Systematics of cybaeid spiders endemic to the Japanese Archipelago, and their historical biogeographic implications (Araneae: Cybaeidae)

The epigean and subterranean spiders of the genus Cybaeus L. Koch, 1868 are distributed in the Holarctic, and are highly diversified in western North America and Japan. More than 100 species have been described from the Japanese Archipelago and several species assemblages have also been recognised among the Japanese Cybaeus based on their morphological similarities. However, their phylogenetic backbone remains unclear. Moreover, genus-level classification of several of the Japanese species appear to be questionable. In this study we performed extensive molecular analyses of the family Cybaeidae in East Asia, mainly the Japanese Archipelago, to clarify their fundamental phylogenetic relationships. We also conducted a divergence time estimation to provide insights into their historical biogeography and evolutionary history. Our divergence dating results indicate that the diversification of the major lineages of the East Asian cybaeid spiders might be related to the opening of the Sea of Japan c .20million years ago. On the basis of the morphological evaluation and obtained phylogenies, some East Asian species formerly placed in Cybaeus are transferred to the cybaeid genera Allocybaeina Bennett, 2020, Sincybaeus Wang & Zhang, 2022 and Cybaeina Chamberlin & Ivie, 1932, i.e. Allocybaeina petegarina (Yaginuma, 1972), comb. nov., Sincybaeus monticola (Kobayashi, 2006), comb. nov., Sincybaeus rarispinosus (Yaginuma, 1970), comb. nov., Sincybaeus yoshiakii (Yaginuma, 1968), comb. nov., and Cybaeina whanseunensis (Paik & Namkung, 1967), comb. nov. Our results clarify the genus-level richness of cybaeids in the Japanese Archipelago for the first time, and reveal the fine-scale phylogenetic relationships of Cybaeus species endemic to the Japanese Archipelago and adjacent regions. ZooBank: urn:lsid:zoobank.org:pub:AF2A3C0E-7F0F-4253-85BA-D995A075F00D.

Spider's Silk as a Potential Source of Antibiotics: An Integrative Review

Spiders produce webs, which are still a largely unexplored source of antibacterial compounds, although the reports of its application in the medical field. Therefore, this study aims to present an integrative review of the antibacterial activity of spider webs. The research was conducted using Google Scholar, Scielo, Web of Science, PubMed, ScienceDirect, Medline EBSCO, LILACS, and Embase. The inclusion criteria were original articles written in English that studied the antibiotic properties of the web or isolated compounds tested. The studies were compared according to the spider species studied, the type of web, treatment of the sample, type of antimicrobial test, and the results obtained. Nine hundred and seventy-three publications were found, and after applying the inclusion and exclusion criteria, sixteen articles were selected. Bacterial inhibition was found in seven studies against various species of bacteria such as Escherichia coli, Klebsiella pneumoniae, Staphylococcus aureus, Salmonella Typhi, Bacillus megaterium, Listeria monocytogenes, Acinetobacter baumannii, Streptococcus pneumoniae, Pasteurella multocida, and Bacillus subtilis. Additionally, there was no apparent relationship between the proximity of the spider species evaluated in the studies and the presence or absence of activity. Methodological problems detected may affected the reproducibility and reliability of the results in some studies, such as the lack of description of the web or microorganism strain, as well as the absence of adequate controls and treatments to sterilize the sample. Spider webs can be a valuable source of antibiotics; however, more studies are needed to confirm the real activity of the web or components involved.

Phylogenomic data exploration with increased sampling provides new insights into the higher-level relationships of butterflies and moths (Lepidoptera)

A robust and stable phylogenetic framework is a fundamental goal of evolutionary biology. As the third largest insect order in the world following Coleoptera and Diptera, Lepidoptera (butterflies and moths) play a central role in almost every terrestrial ecosystem as indicators of environmental change and serve as important models for biologists exploring questions related to ecology and evolutionary biology. However, for such a charismatic insect group, the higher-level phylogenetic relationships among its superfamilies are still poorly resolved. Compared to earlier phylogenomic studies, we increased taxon sampling among Lepidoptera (37 superfamilies and 68 families containing 263 taxa) and acquired a series of large amino-acid datasets from 69,680 to 400,330 for phylogenomic reconstructions. Using these datasets, we explored the effect of different taxon sampling with significant increases in the number of included genes on tree topology by considering a series of systematic errors using maximum-likelihood (ML) and Bayesian inference (BI) methods. Moreover, we also tested the effectiveness in topology robustness among the three ML-based models. The results showed that taxon sampling is an important determinant in tree robustness of accurate lepidopteran phylogenetic estimation. Long-branch attraction (LBA) caused by site-wise heterogeneity is a significant source of bias giving rise to unstable positions of ditrysian groups in phylogenomic reconstruction. Phylogenetic inference showed the most comprehensive framework to reveal the relationships among lepidopteran superfamilies, and presented some newly relationships with strong supports (Papilionoidea was sister to Gelechioidea and Immoidea was sister to Galacticoidea, respectively), but limited by taxon sampling, the relationships within the species-rich and relatively rapid radiation Ditrysia and especially Apoditrysia remain poorly resolved, which need to increase taxon sampling for further phylogenomic reconstruction. The present study demonstrates that taxon sampling is an important determinant for an accurate lepidopteran tree of life and provides some essential insights for future lepidopteran phylogenomic studies.

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.

Allometry and ecology shape eye size evolution in spiders

Eye size affects many aspects of visual function, but eyes are costly to grow and maintain. The allometry of eyes can provide insight into this trade-off, but this has mainly been explored in species that have two eyes of equal size. By contrast, animals possessing larger visual systems can exhibit variable eye sizes within individuals. Spiders have up to four pairs of eyes whose sizes vary dramatically, but their ontogenetic, static, and evolutionary allometry has not yet been studied in a comparative context. We report variable dynamics in eye size across 1,098 individuals in 39 species and 8 families, indicating selective pressures and constraints driving the evolution of different eye pairs and lineages. Supplementing our sampling with a recently published phylogenetically comprehensive dataset, we confirmed these findings across more than 400 species; found that ecological factors such as visual hunting, web building, and circadian activity correlate with eye diameter; and identified significant allometric shifts across spider phylogeny using an unbiased approach, many of which coincide with visual hunting strategies. The modular nature of the spider visual system provides additional degrees of freedom and is apparent in the strong correlations between maximum/minimum investment and interocular variance and three key ecological factors. Our analyses suggest an antagonistic relationship between the anterior and posterior eye pairs. These findings shed light on the relationship between spider visual systems and their diverse ecologies and how spiders exploit their modular visual systems to balance selective pressures and optical and energetic constraints.

Scoutknife: A naïve, whole genome informed phylogenetic robusticity metric

Background: The phylogenetic bootstrap, first proposed by Felsenstein in 1985, is a critically important statistical method in assessing the robusticity of phylogenetic datasets. Core to its concept was the use of pseudo sampling - assessing the data by generating new replicates derived from the initial dataset that was used to generate the phylogeny. In this way, phylogenetic support metrics could overcome the lack of perfect, infinite data. With infinite data, however, it is possible to sample smaller replicates directly from the data to obtain both the phylogeny and its statistical robusticity in the same analysis. Due to the growth of whole genome sequencing, the depth and breadth of our datasets have greatly expanded and are set to only expand further. With genome-scale datasets comprising thousands of genes, we can now obtain a proxy for infinite data. Accordingly, we can potentially abandon the notion of pseudo sampling and instead randomly sample small subsets of genes from the thousands of genes in our analyses. Methods: We introduce Scoutknife, a jackknife-style subsampling implementation that generates 100 datasets by randomly sampling a small number of genes from an initial large-gene dataset to jointly establish both a phylogenetic hypothesis and assess its robusticity. We assess its effectiveness by using 18 previously published datasets and 100 simulation studies. Results: We show that Scoutknife is conservative and informative as to conflicts and incongruence across the whole genome, without the need for subsampling based on traditional model selection criteria. Conclusions: Scoutknife reliably achieves comparable results to selecting the best genes on both real and simulation datasets, while being resistant to the potential biases caused by selecting for model fit. As the amount of genome data grows, it becomes an even more exciting option to assess the robusticity of phylogenetic hypotheses.

Insight into the adaptive role of arachnid genome-wide duplication through chromosome-level genome assembly of the Western black widow spider

Although spiders are one of the most diverse groups of arthropods, the genetic architecture of their evolutionary adaptations is largely unknown. Specifically, ancient genome-wide duplication occurring during arachnid evolution ~450 mya resulted in a vast assembly of gene families, yet the extent to which selection has shaped this variation is understudied. To aid in comparative genome sequence analyses, we provide a chromosome-level genome of the Western black widow spider (Latrodectus hesperus)-a focus due to its silk properties, venom applications, and as a model for urban adaptation. We used long-read and Hi-C sequencing data, combined with transcriptomes, to assemble 14 chromosomes in a 1.46 Gb genome, with 38,393 genes annotated, and a BUSCO score of 95.3%. Our analyses identified high repetitive gene content and heterozygosity, consistent with other spider genomes, which has led to challenges in genome characterization. Our comparative evolutionary analyses of eight genomes available for species within the Araneoidea group (orb weavers and their descendants) identified 1,827 single-copy orthologs. Of these, 155 exhibit significant positive selection primarily associated with developmental genes, and with traits linked to sensory perception. These results support the hypothesis that several traits unique to spiders emerged from the adaptive evolution of ohnologs-or retained ancestrally duplicated genes-from ancient genome-wide duplication. These comparative spider genome analyses can serve as a model to understand how positive selection continually shapes ancestral duplications in generating novel traits today within and between diverse taxonomic groups.

Metazoa-level USCOs as markers in species delimitation and classification

Metazoa-level universal single-copy orthologs (mzl-USCOs) are universally applicable markers for DNA taxonomy in animals that can replace or supplement single-gene barcodes. Previously, mzl-USCOs from target enrichment data were shown to reliably distinguish species. Here, we tested whether USCOs are an evenly distributed, representative sample of a given metazoan genome and therefore able to cope with past hybridization events and incomplete lineage sorting. This is relevant for coalescent-based species delimitation approaches, which critically depend on the assumption that the investigated loci do not exhibit autocorrelation due to physical linkage. Based on 239 chromosome-level assembled genomes, we confirmed that mzl-USCOs are genetically unlinked for practical purposes and a representative sample of a genome in terms of reciprocal distances between USCOs on a chromosome and of distribution across chromosomes. We tested the suitability of mzl-USCOs extracted from genomes for species delimitation and phylogeny in four case studies: Anopheles mosquitos, Drosophila fruit flies, Heliconius butterflies and Darwin's finches. In almost all instances, USCOs allowed delineating species and yielded phylogenies that corresponded to those generated from whole genome data. Our phylogenetic analyses demonstrate that USCOs may complement single-gene DNA barcodes and provide more accurate taxonomic inferences. Combining USCOs from sources that used different versions of ortholog reference libraries to infer marker orthology may be challenging and, at times, impact taxonomic conclusions. However, we expect this problem to become less severe as the rapidly growing number of reference genomes provides a better representation of the number and diversity of organismal lineages.

Facilitating taxonomy and phylogenetics: An informative and cost-effective protocol integrating long amplicon PCRs and third-generation sequencing

Phylogenetic inference has become a standard technique in integrative taxonomy and systematics, as well as in biogeography and ecology. DNA barcodes are often used for phylogenetic inference, despite being strongly limited due to their low number of informative sites. Also, because current DNA barcodes are based on a fraction of a single, fast-evolving gene, they are highly unsuitable for resolving deeper phylogenetic relationships due to saturation. In recent years, methods that analyse hundreds and thousands of loci at once have improved the resolution of the Tree of Life, but these methods require resources, experience and molecular laboratories that most taxonomists do not have. This paper introduces a PCR-based protocol that produces long amplicons of both slow- and fast-evolving unlinked mitochondrial and nuclear gene regions, which can be sequenced by the affordable and portable ONT MinION platform with low infrastructure or funding requirements. As a proof of concept, we inferred a phylogeny of a sample of 63 spider species from 20 families using our proposed protocol. The results were overall consistent with the results from approaches based on hundreds and thousands of loci, while requiring just a fraction of the cost and labour of such approaches, making our protocol accessible to taxonomists worldwide.

The polyketide to fatty acid transition in the evolution of animal lipid metabolism

Animals synthesize simple lipids using a distinct fatty acid synthase (FAS) related to the type I polyketide synthase (PKS) enzymes that produce complex specialized metabolites. The evolutionary origin of the animal FAS and its relationship to the diversity of PKSs remain unclear despite the critical role of lipid synthesis in cellular metabolism. Recently, an animal FAS-like PKS (AFPK) was identified in sacoglossan molluscs. Here, we explore the phylogenetic distribution of AFPKs and other PKS and FAS enzymes across the tree of life. We found AFPKs widely distributed in arthropods and molluscs (>6300 newly described AFPK sequences). The AFPKs form a clade with the animal FAS, providing an evolutionary link bridging the type I PKSs and the animal FAS. We found molluscan AFPK diversification correlated with shell loss, suggesting AFPKs provide a chemical defense. Arthropods have few or no PKSs, but our results indicate AFPKs contributed to their ecological and evolutionary success by facilitating branched hydrocarbon and pheromone biosynthesis. Although animal metabolism is well studied, surprising new metabolic enzyme classes such as AFPKs await discovery.

Burma Terrane Amber Fauna Shows Connections to Gondwana and Transported Gondwanan Lineages to the Northern Hemisphere (Araneae: Palpimanoidea)

Burmese amber is a significant source of fossils that documents the mid-Cretaceous biota. This deposit was formed around 99 Ma on the Burma Terrane, which broke away from Gondwana and later collided with Asia, although the timing is disputed. Palpimanoidea is a dispersal-limited group that was a dominant element of the Mesozoic spider fauna, and has an extensive fossil record, particularly from Burmese amber. Using morphological and molecular data, evolutionary relationships of living and fossil Palpimanoidea are examined. Divergence dating with fossils as terminal tips shows timing of diversification is contemporaneous with continental breakup.Ancestral range estimations show widespread ancestral ranges that divide into lineages that inherit different Pangean fragments, consistent with vicariance. Our results suggest that the Burmese amber fauna has ties to Gondwana due to a historical connection in the Early Cretaceous, and that the Burma Terrane facilitated biotic exchange by transporting lineages from Gondwana into the Holarctic in the Cretaceous.

Evolution of the Spider Homeobox Gene Repertoire by Tandem and Whole Genome Duplication

Gene duplication generates new genetic material that can contribute to the evolution of gene regulatory networks and phenotypes. Duplicated genes can undergo subfunctionalization to partition ancestral functions and/or neofunctionalization to assume a new function. We previously found there had been a whole genome duplication (WGD) in an ancestor of arachnopulmonates, the lineage including spiders and scorpions but excluding other arachnids like mites, ticks, and harvestmen. This WGD was evidenced by many duplicated homeobox genes, including two Hox clusters, in spiders. However, it was unclear which homeobox paralogues originated by WGD versus smaller-scale events such as tandem duplications. Understanding this is a key to determining the contribution of the WGD to arachnopulmonate genome evolution. Here we characterized the distribution of duplicated homeobox genes across eight chromosome-level spider genomes. We found that most duplicated homeobox genes in spiders are consistent with an origin by WGD. We also found two copies of conserved homeobox gene clusters, including the Hox, NK, HRO, Irx, and SINE clusters, in all eight species. Consistently, we observed one copy of each cluster was degenerated in terms of gene content and organization while the other remained more intact. Focussing on the NK cluster, we found evidence for regulatory subfunctionalization between the duplicated NK genes in the spider Parasteatoda tepidariorum compared to their single-copy orthologues in the harvestman Phalangium opilio. Our study provides new insights into the relative contributions of multiple modes of duplication to the homeobox gene repertoire during the evolution of spiders and the function of NK genes.

Single-cell transcriptomics reveals the brain evolution of web-building spiders

Spiders are renowned for their efficient capture of flying insects using intricate aerial webs. How the spider nervous systems evolved to cope with this specialized hunting strategy and various environmental clues in an aerial space remains unknown. Here we report a brain-cell atlas of >30,000 single-cell transcriptomes from a web-building spider (Hylyphantes graminicola). Our analysis revealed the preservation of ancestral neuron types in spiders, including the potential coexistence of noradrenergic and octopaminergic neurons, and many peptidergic neuronal types that are lost in insects. By comparing the genome of two newly sequenced plesiomorphic burrowing spiders with three aerial web-building spiders, we found that the positively selected genes in the ancestral branch of web-building spiders were preferentially expressed (42%) in the brain, especially in the three mushroom body-like neuronal types. By gene enrichment analysis and RNAi experiments, these genes were suggested to be involved in the learning and memory pathway and may influence the spiders' web-building and hunting behaviour. Our results provide key sources for understanding the evolution of behaviour in spiders and reveal how molecular evolution drives neuron innovation and the diversification of associated complex behaviours.

Advances in the reconstruction of the spider tree of life: A roadmap for spider systematics and comparative studies

In the last decade and a half, advances in genetic sequencing technologies have revolutionized systematics, transforming the field from studying morphological characters or a few genetic markers, to genomic datasets in the phylogenomic era. A plethora of molecular phylogenetic studies on many taxonomic groups have come about, converging on, or refuting prevailing morphology or legacy-marker-based hypotheses about evolutionary affinities. Spider systematics has been no exception to this transformation and the inter-relationships of several groups have now been studied using genomic data. About 51 500 extant spider species have been described, all with a conservative body plan, but innumerable morphological and behavioural peculiarities. Inferring the spider tree of life using morphological data has been a challenging task. Molecular data have corroborated many hypotheses of higher-level relationships, but also resulted in new groups that refute previous hypotheses. In this review, we discuss recent advances in the reconstruction of the spider tree of life and highlight areas where additional effort is needed with potential solutions. We base this review on the most comprehensive spider phylogeny to date, representing 131 of the 132 spider families. To achieve this sampling, we combined six Sanger-based markers with newly generated and publicly available genome-scale datasets. We find that some inferred relationships between major lineages of spiders (such as Austrochiloidea, Palpimanoidea and Synspermiata) are robust across different classes of data. However, several new hypotheses have emerged with different classes of molecular data. We identify and discuss the robust and controversial hypotheses and compile this blueprint to design future studies targeting systematic revisions of these problematic groups. We offer an evolutionary framework to explore comparative questions such as evolution of venoms, silk, webs, morphological traits and reproductive strategies.

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.

Phylogenomics illuminates the evolution of orb webs, respiratory systems and the biogeographic history of the world's smallest orb-weaving spiders (Araneae, Araneoidea, Symphytognathoids)

The miniature orb weaving spiders (symphytognathoids) are a group of small spiders (< 2 mm), including the smallest adult spider Patu digua (0.37 mm in body length), that have been classified into five families. The species of one of its constituent lineages, the family Anapidae, build a remarkable diversity of webs (ranging from orbs to sheet webs and irregular tangles) and even include a webless kleptoparasitic species. Anapids are also exceptional because of the extraordinary diversity of their respiratory systems. The phylogenetic relationships of symphytognathoid families have been recalcitrant with different classes of data, such as, monophyletic with morphology and its concatenation with Sanger-based six markers, paraphyletic (including a paraphyletic Anapidae) with solely Sanger-based six markers, and polyphyletic with transcriptomes. In this study, we capitalized on a large taxonomic sampling of symphytognathoids, focusing on Anapidae, and using de novo sequenced ultraconserved elements (UCEs) combined with UCEs recovered from available transcriptomes and genomes. We evaluated the conflicting relationships using a variety of support metrics and topology tests. We found support for the phylogenetic hypothesis proposed using morphology to obtain the "symphytognathoids'' clade, Anterior Tracheal System (ANTS) Clade and monophyly of the family Anapidae. Anapidae can be divided into three major lineages, the Vichitra Clade (including Teutoniella, Holarchaea, Sofanapis and Acrobleps), the subfamily Micropholcommatinae and the Orb-weaving anapids (Owa) Clade. Biogeographic analyses reconstructed a hypothesis of multiple long-distance transoceanic dispersal events, potentially influenced by the Antarctic Circumpolar Current and West Wind Drift. In symphytognathoids, the ancestral anterior tracheal system transformed to book lungs four times and reduced book lungs five times. The posterior tracheal system was lost six times. The orb web structure was lost four times independently and transformed into sheet web once.

Comparative Mitogenomics of Jumping Spiders with First Complete Mitochondrial Genomes of Euophryini (Araneae: Salticidae)

Salticidae is the most species-rich family of spiders with diverse morphology, ecology and behavior. However, the characteristics of the mitogenomes within this group are poorly understood with relatively few well-characterized complete mitochondrial genomes. In this study, we provide completely annotated mitogenomes for Corythalia opima and Parabathippus shelfordi, which represent the first complete mitogenomes of the tribe Euophryini of Salticidae. The features and characteristics of the mitochondrial genomes are elucidated for Salticidae by thoroughly comparing the known well-characterized mitogenomes. The gene rearrangement between trnL2 and trnN was found in two jumping spider species, Corythalia opima and Heliophanus lineiventris Simon, 1868. Additionally, the rearrangement of nad1 to between trnE and trnF found in Asemonea sichuanensis Song & Chai, 1992 is the first protein-coding gene rearrangement in Salticidae, which may have an important phylogenetic implication for the family. Tandem repeats of various copy numbers and lengths were discovered in three jumping spider species. The codon usage analyses showed that the evolution of codon usage bias in salticid mitogenomes was affected by both selection and mutational pressure, but selection may have played a more important role. The phylogenetic analyses provided insight into the taxonomy of Colopsus longipalpis (Żabka, 1985). The data presented in this study will improve our understanding of the evolution of mitochondrial genomes within Salticidae.

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

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

Evolution and comparative morphology of raptorial feet in spiders

Spiders are among the most diverse animals, which developed different morphological and behavioral traits for capturing prey. We studied the anatomy and functionality of the rare and apomorphic raptorial spider feet using 3D reconstruction modeling, among other imaging techniques. The evolutionary reconstruction of the raptorial feet (tarsus plus pretarsus) features using a composite tree of spiders, indicating that similar traits emerged three times independently in Trogloraptoridae, Gradungulinae, and Doryonychus raptor (Tetragnathidae). The characteristics defining the raptorial feet are an interlocked complex merging of the base of the elongated prolateral claw with the pretarsal sclerotized ring, with the former clasping against the tarsus. Raptorial feet even flex over robust raptorial macrosetae forming a reduced tarsal version of a catching basket to encase prey during hunting. Our results show that Celaeniini (Araneidae) and Heterogriffus berlandi (Thomisidae), taxa previously compared with raptorial spiders, lack the raptorial feet key characteristics and the tarsal-catching basket. We make predictions about the possible behavior of the abovementioned taxa that will need to be tested by observing living specimens. We conclude that multiple morphological tarsal and pretarsal micro-structures define the raptorial foot functional unit and recommend a comprehensive evaluation before assigning this configuration to any spider taxa.

Assessing the impact of fire on spiders through a global comparative analysis

In many regions fire regimes are changing due to anthropogenic factors. Understanding the responses of species to fire can help to develop predictive models and inform fire management decisions. Spiders are a diverse and ubiquitous group and can offer important insights into the impacts of fire on invertebrates and whether these depend on environmental factors, phylogenetic history or functional traits. We conducted phylogenetic comparative analyses of data from studies investigating the impacts of fire on spiders. We investigated whether fire affects spider abundance or presence and whether ecologically relevant traits or site-specific factors influence species' responses to fire. Although difficult to make broad generalizations about the impacts of fire due to variation in site- and fire-specific factors, we find evidence that short fire intervals may be a threat to some spiders, and that fire affects abundance and species compositions in forests relative to other vegetation types. Orb and sheet web weavers were also more likely to be absent after fire than ambush hunters, ground hunters and other hunters suggesting functional traits may affect responses. Finally, we show that analyses of published data can be used to detect broad-scale patterns and provide an alternative to traditional meta-analytical approaches.

A novel probe set for the phylogenomics and evolution of RTA spiders

Spiders are important models for evolutionary studies of web building, sexual selection and adaptive radiation. The recent development of probes for UCE (ultra-conserved element)-based phylogenomic studies has shed light on the phylogeny and evolution of spiders. However, the two available UCE probe sets for spider phylogenomics (Spider and Arachnida probe sets) have relatively low capture efficiency within spiders, and are not optimized for the retrolateral tibial apophysis (RTA) clade, a hyperdiverse lineage that is key to understanding the evolution and diversification of spiders. In this study, we sequenced 15 genomes of species in the RTA clade, and using eight reference genomes, we developed a new UCE probe set (41 845 probes targeting 3802 loci, labelled as the RTA probe set). The performance of the RTA probes in resolving the phylogeny of the RTA clade was compared with the Spider and Arachnida probes through an in-silico test on 19 genomes. We also tested the new probe set empirically on 28 spider species of major spider lineages. The results showed that the RTA probes recovered twice and four times as many loci as the other two probe sets, and the phylogeny from the RTA UCEs provided higher support for certain relationships. This newly developed UCE probe set shows higher capture efficiency empirically and is particularly advantageous for phylogenomic and evolutionary studies of RTA clade and jumping spiders.

Geometric regularity in webs of non-orb-weaving spiders

Geometric regularity of spider webs has been intensively studied in orb-weaving spiders, although it is not exclusive of orb weavers. Here, we document the geometrically regular, repetitive elements in the webs of the non-orb-weaving groups Leptonetidae and Telemidae for the first time. Similar to orb weavers, we found areas with regularly spaced parallel lines in the webs of Calileptoneta helferi, Sulcia sp., and cf. Pinelema sp. Furthermore, we provide a detailed account of the regular webs of Ochyrocera (Ochyroceratidae). The sections of the web with regularly disposed parallel lines are built as U-shaped modules reminiscent of orb webs. It has been suggested that the regularly spaced parallel lines in the webs of Ochyroceratidae and Psilodercidae may be produced in a single sweep of their posterior lateral spinnerets, which have regularly spaced aciniform gland spigots, perhaps involving expansion of the spinnerets. To test this hypothesis, we compared the spacing between parallel lines with the spacing between spigots, searched for expansible membranes in the spinnerets, and examined the junctions of regularly spaced lines. The distance between parallel lines was 10-20 times the distance between spigots, and we found no expansible membranes, and the intersection of parallel lines are cemented, which opposes the single sweep hypothesis. Furthermore, we found cues of viscid silk in the parallel lines of the psilodercid Althepus and broadened piriform gland spigots that may be responsible of its production. Finally, we evaluated the presence or absence of geometrically regular web elements across the spider tree of life. We found reports of regular webs in 31 spider families, including 20 families that are not orb weavers and hypothesize that the two basic aspects of regularity (parallel lines spaced at regular intervals, and radial lines spaced at regular angles) probably appeared many times in the evolution of spiders.

The deep-rooted origin of disulfide-rich spider venom toxins

Spider venoms are a complex concoction of enzymes, polyamines, inorganic salts, and disulfide-rich peptides (DRPs). Although DRPs are widely distributed and abundant, their bevolutionary origin has remained elusive. This knowledge gap stems from the extensive molecular divergence of DRPs and a lack of sequence and structural data from diverse lineages. By evaluating DRPs under a comprehensive phylogenetic, structural and evolutionary framework, we have not only identified 78 novel spider toxin superfamilies but also provided the first evidence for their common origin. We trace the origin of these toxin superfamilies to a primordial knot - which we name 'Adi Shakti', after the creator of the Universe according to Hindu mythology - 375 MYA in the common ancestor of Araneomorphae and Mygalomorphae. As the lineages under evaluation constitute nearly 60% of extant spiders, our findings provide fascinating insights into the early evolution and diversification of the spider venom arsenal. Reliance on a single molecular toxin scaffold by nearly all spiders is in complete contrast to most other venomous animals that have recruited into their venoms diverse toxins with independent origins. By comparatively evaluating the molecular evolutionary histories of araneomorph and mygalomorph spider venom toxins, we highlight their contrasting evolutionary diversification rates. Our results also suggest that venom deployment (e.g. prey capture or self-defense) influences evolutionary diversification of DRP toxin superfamilies.

Molecular phylogeny of the tropical wandering spiders (Araneae, Ctenidae) and the evolution of eye conformation in the RTA clade

Tropical wandering spiders (Ctenidae) are a diverse group of cursorial predators with its greatest species richness in the tropics. Traditionally, Ctenidae are diagnosed based on the presence of eight eyes arranged in three rows (a 2-4-2 pattern). We present a molecular phylogeny of Ctenidae, including for the first time representatives of all of its subfamilies. The molecular phylogeny was inferred using five nuclear (histone H3, 28S, 18S, Actin and ITS-2) and four mitochondrial (NADH, COI, 12S and 16S) markers. The final matrix includes 259 terminals, 103 of which belong to Ctenidae and represent 28 of the current 49 described genera. We estimated divergence times by including fossils as calibration points and biogeographic events, and used the phylogenetic hypothesis obtained to reconstruct the evolution of the eye conformation in the retrolateral tibial apophysis (RTA) clade. Ctenidae and its main lineages originated during the Paleocene-Eocene and have diversified in the tropics since then. However, in some analyses Ctenidae was recovered as polyphyletic as the genus Ancylometes Bertkau, 1880 was placed as sister to Oxyopidae. Except for Acantheinae, in which the type genus Acantheis Thorell, 1891 is placed inside Cteninae, the four recognized subfamilies of Ctenidae are monophyletic in most analyses. The ancestral reconstruction of the ocular conformation in the retrolateral tibial apophysis clade suggests that the ocular pattern of Ctenidae has evolved convergently seven times and that it has originated from ocular conformations of two rows of four eyes (4-4) and the ocular pattern of lycosids (4-2-2). We also synonymize the monotypic genus Parabatinga Polotov & Brescovit, 2009 with Centroctenus Mello-Leitão, 1929. We discuss some of the putative morphological synapomorphies of the main ctenid lineages within the phylogenetic framework offered by the molecular phylogenetic results of the study.

Transcriptome analysis of the spider Phonotimpus pennimani reveals novel toxin transcripts

Background

Phonotimpus pennimani (Araneae, Phrurolithidae) is a small-sized (3-5 mm) spider endemic to the Tacaná volcano in Chiapas, Mexico, where it is found in soil litter of cloud forests and coffee plantations. Its venom composition has so far not been investigated, partly because it is not a species of medical significance. However, it does have an important impact on the arthropod populations of its natural habitat.

Methods

Specimens were collected in Southeastern Mexico (Chiapas) and identified taxonomically by morphological characteristics. A partial sequence from the mitochondrial gene coxI was amplified. Sequencing on the Illumina platform of a transcriptome library constructed from 12 adult specimens revealed 25 toxin or toxin-like genes. Transcripts were validated (RT-qPCR) by assessing the differential expression of the toxin-like PpenTox1 transcript and normalising with housekeeping genes.

Results

Analysis of the coxI-gene revealed a similarity to other species of the family Phrurolithidae. Transcriptome analysis also revealed similarity with venom components of species from the families Ctenidae, Lycosidae, and Sicariidae. Expression of the toxin-like PpenTox1 gene was different for each developmental stage (juvenile or adult) and also for both sexes (female or male). Additionally, a partial sequence was obtained for the toxin-like PpenTox1 from DNA.

Conclusion

Data from the amplification of the mitochondrial coxI gene confirmed that P. pennimani belongs to the family Phrurolithidae. New genes and transcripts coding for venom components were identified.

Chromosome-level genome and the identification of sex chromosomes in Uloborus diversus

The orb web is a remarkable example of animal architecture that is observed in families of spiders that diverged over 200 million years ago. While several genomes exist for araneid orb-weavers, none exist for other orb-weaving families, hampering efforts to investigate the genetic basis of this complex behavior. Here we present a chromosome-level genome assembly for the cribellate orb-weaving spider Uloborus diversus. The assembly reinforces evidence of an ancient arachnid genome duplication and identifies complete open reading frames for every class of spidroin gene, which encode the proteins that are the key structural components of spider silks. We identified the 2 X chromosomes for U. diversus and identify candidate sex-determining loci. This chromosome-level assembly will be a valuable resource for evolutionary research into the origins of orb-weaving, spidroin evolution, chromosomal rearrangement, and chromosomal sex determination in spiders.

The evolutionary history of cribellate orb-weaver capture thread spidroins

Background

Spiders have evolved two types of sticky capture threads: one with wet adhesive spun by ecribellate orb-weavers and another with dry adhesive spun by cribellate spiders. The evolutionary history of cribellate capture threads is especially poorly understood. Here, we use genomic approaches to catalog the spider-specific silk gene family (spidroins) for the cribellate orb-weaver Uloborus diversus.

Results

We show that the cribellar spidroin, which forms the puffy fibrils of cribellate threads, has three distinct repeat units, one of which is conserved across cribellate taxa separated by ~ 250 Mya. We also propose candidates for a new silk type, paracribellar spidroins, which connect the puffy fibrils to pseudoflagelliform support lines. Moreover, we describe the complete repeat architecture for the pseudoflagelliform spidroin (Pflag), which contributes to extensibility of pseudoflagelliform axial fibers.

Conclusions

Our finding that Pflag is closely related to Flag, supports homology of the support lines of cribellate and ecribellate capture threads. It further suggests an evolutionary phase following gene duplication, in which both Flag and Pflag were incorporated into the axial lines, with subsequent loss of Flag in uloborids, and increase in expression of Flag in ecribellate orb-weavers, explaining the distinct mechanical properties of the axial lines of these two groups.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12862-022-02042-5.

Biogeography and eye size evolution of the ogre-faced spiders

Net-casting spiders (Deinopidae) comprise a charismatic family with an enigmatic evolutionary history. There are 67 described species of deinopids, placed among three genera, Deinopis, Menneus, and Asianopis, that are distributed globally throughout the tropics and subtropics. Deinopis and Asianopis, the ogre-faced spiders, are best known for their giant light-capturing posterior median eyes (PME), whereas Menneus does not have enlarged PMEs. Molecular phylogenetic studies have revealed discordance between morphology and molecular data. We employed a character-rich ultra-conserved element (UCE) dataset and a taxon-rich cytochrome-oxidase I (COI) dataset to reconstruct a genus-level phylogeny of Deinopidae, aiming to investigate the group's historical biogeography, and examine PME size evolution. Although the phylogenetic results support the monophyly of Menneus and the single reduction of PME size in deinopids, these data also show that Deinopis is not monophyletic. Consequently, we formally transfer 24 Deinopis species to Asianopis; the transfers comprise all of the African, Australian, South Pacific, and a subset of Central American and Mexican species. Following the divergence of Eastern and Western deinopids in the Cretaceous, Deinopis/Asianopis dispersed from Africa, through Asia and into Australia with its biogeographic history reflecting separation of Western Gondwana as well as long-distance dispersal events.

Lineage-specific, fast-evolving GATA-like gene regulates zygotic gene activation to promote endoderm specification and pattern formation in the Theridiidae spider

Background

The process of early development varies across the species-rich phylum Arthropoda. Owing to the limited research strategies for dissecting lineage-specific processes of development in arthropods, little is known about the variations in early arthropod development at molecular resolution. The Theridiidae spider, Parasteatoda tepidariorum, has its genome sequenced and could potentially contribute to dissecting early embryonic processes.

Results

We present genome-wide identification of candidate genes that exhibit locally restricted expression in germ disc forming stage embryos of P. tepidariorum, based on comparative transcriptomes of isolated cells from different regions of the embryo. A subsequent pilot screen by parental RNA interference identifies three genes required for body axis formation. One of them is a GATA-like gene that has been fast evolving after duplication and divergence from a canonical GATA family gene. This gene is designated fuchi nashi (fuchi) after its knockdown phenotypes, where the cell movement toward the formation of a germ disc was reversed. fuchi expression occurs in cells outside a forming germ disc and persists in the endoderm. Transcriptome and chromatin accessibility analyses of fuchi pRNAi embryos suggest that early fuchi activity regulates chromatin state and zygotic gene activation to promote endoderm specification and pattern formation. We also show that there are many uncharacterized genes regulated by fuchi.

Conclusions

Our genome-based research using an arthropod phylogenetically distant from Drosophila identifies a lineage-specific, fast-evolving gene with key developmental roles in one of the earliest, genome-wide regulatory events, and allows for molecular exploration of the developmental variations in early arthropod embryos.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12915-022-01421-0.

Phylogeny and biogeography of the ancient spider family Filistatidae (Araneae) is consistent both with long-distance dispersal and vicariance following continental drift

Filistatids, the crevice weavers, are an ancient family of cribellate spiders without extant close relatives. As one of the first lineages of araneomorph spiders, they present a complicated mixture of primitive and derived characters that make them a key taxon to elucidate the phylogeny of spiders, as well as the evolution of phenotypic characters in this group. Their moderate diversity (187 species in 19 genera) is distributed mainly in arid and semi-arid subtropical zones of all continents, except Antarctica. The objective of this paper is to generate a comprehensive phylogenetic hypothesis for this family to advance the understanding of its morphological evolution and biogeography, as well as lay the basis for a natural classification scheme. By studying the morphology using optical and electronic microscopy techniques, we produced a matrix of 302 morphological characters coded for a sample of 103 species of filistatids chosen to represent the phylogenetic diversity of the family. In addition, we included sequences of four molecular markers (COI, 16S, H3 and 28S; 3787 aligned positions) of 70 filistatid species. The analysis of the data (morphological, molecular, and combined) consistently indicates the separation of the Filistatidae into two subfamilies, Prithinae and Filistatinae, in addition to supporting several groups of genera: Filistata, Zaitunia and an undescribed genus from Madagascar; Sahastata and Kukulcania; all Prithinae except Filistatinella and Microfilistata; Antilloides and Filistatoides; a large Old World group including Pritha, Tricalamus, Afrofilistata, Labahitha, Yardiella, Wandella and putative new genera; and a South American group formed by Lihuelistata, Pikelinia and Misionella. Pholcoides is transferred to Filistatinae and Microfilistata is transferred to Prithinae, and each represents the sister group to the remaining genera of its own subfamily. Most genera are valid, although Pikelinia is paraphyletic with respect to Misionella, so we consider the two genera as synonyms and propose a few new generic combinations. Considering the new phylogenetic hypothesis, we discuss the evolution of some morphological character systems and the biogeography of the family. The ages of divergence between clades were estimated using a total-evidence tip-dating approach by including fossils of Filistatidae and early spider clades; this approach resulted in younger age estimates than those obtained with traditional node-dating. Filistatidae is an ancient family that started diversifying in the Mesozoic and most genera date to the Cretaceous. Clades displaying transcontinental distributions were most likely affected by continental drift, but at least one clade shows unequivocal signs of transoceanic long-distance dispersal.

Evolutionary morphology of sperm in pholcid spiders (Pholcidae, Synspermiata)

Background

Pholcidae represent one of the largest and most diverse spider families and have been subject to various studies regarding behavior and reproductive biology. In contrast to the solid knowledge on phylogeny and general reproductive morphology, the primary male reproductive system is strongly understudied, as it has been addressed only for few species. Those studies however suggested a high diversity of sperm and seminal secretions across the family. To address this disparity and reconstruct the evolution of sperm traits, we investigate the primary male reproductive system of pholcid spiders by means of light, X-ray, and transmission electron microscopy using a comprehensive taxon sampling with 46 species from 33 genera, representing all five subfamilies.

Results

Our data show a high disparity of sperm morphology and seminal secretions within pholcids. We document several sperm characters that are unique for pholcids, such as a helical band (Pholcinae) or a lamellate posterior centriolar adjunct material (Modisiminae). Character mapping revealed several putative synapomorphies for individual taxa. With regard to sperm transfer forms, we found that synspermia occur only in the subfamily Ninetinae, whereas the other subfamilies have cleistospermia. In several species with cleistospermia, we demonstrate that spermatids remain fused until late stages of spermiogenesis before ultimately separating shortly before the coiling process. Additionally, we explored the previously hypothesized correlation between sperm size and minimum diameter of the spermophor in the male palpal organ. We show that synspermia differ strongly in size whereas cleistospermia are rather uniform, but neither transfer form is positively correlated with the diameter of the spermophor.

Conclusions

Our data revealed a dynamic evolution of sperm characters, with convergences across all subfamilies and a high level of homoplasy. The present diversity can be related to subfamily level and allows for assignments of specific subtypes of spermatozoa. Our observations support the idea that Ninetinae are an ancestral clade within Pholcidae that have retained synspermia and that synspermia represent the ancestral sperm transfer form of Pholcidae.

Comparative mitogenomic analyses provide evolutionary insights into the retrolateral tibial apophysis clade (Araneae: Entelegynae)

The retrolateral tibial apophysis (RTA) clade is the largest spider lineage within Araneae. To better understand the diversity and evolution, we newly determined mitogenomes of ten RTA species from six families and performed a comparative mitogenomics analysis by combining them with 40 sequenced RTA mitogenomes available on GenBank. The ten mitogenomes encoded 37 typical mitochondrial genes and included a large non-coding region (putative control region). Nucleotide composition and codon usage were well conserved within the RTA clade, whereas diversity in sequence length and structural features was observed in control region. A reversal of strand asymmetry in nucleotide composition, i.e., negative AT-skews and positive GC-skews, was observed in each RTA species, likely resulting from mitochondrial gene rearrangements. All protein-coding genes were evolving under purifying selection, except for atp8 whose Ka/Ks was larger than 1, possibly due to positive selection or selection relaxation. Both mutation pressure and natural selection might contribute to codon usage bias of 13 protein-coding genes in the RTA lineage. Phylogenetic analyses based on mitogenomic data recovered a family-level phylogeny within the RTA; {[(Oval calamistrum clade, Dionycha), Marronoid clade], Sparassidae}. This study characterized RTA mitogenomes and provided some new insights into the phylogeny and evolution of the RTA clade.

Improving the Annotation of the Venom Gland Transcriptome of Pamphobeteus verdolaga, Prospecting Novel Bioactive Peptides

Spider venoms constitute a trove of novel peptides with biotechnological interest. Paucity of next-generation-sequencing (NGS) data generation has led to a description of less than 1% of these peptides. Increasing evidence supports the underestimation of the assembled genes a single transcriptome assembler can predict. Here, the transcriptome of the venom gland of the spider Pamphobeteus verdolaga was re-assembled, using three free access algorithms, Trinity, SOAPdenovo-Trans, and SPAdes, to obtain a more complete annotation. Assembler's performance was evaluated by contig number, N50, read representation on the assembly, and BUSCO's terms retrieval against the arthropod dataset. Out of all the assembled sequences with all software, 39.26% were common between the three assemblers, and 27.88% were uniquely assembled by Trinity, while 27.65% were uniquely assembled by SPAdes. The non-redundant merging of all three assemblies' output permitted the annotation of 9232 sequences, which was 23% more when compared to each software and 28% more when compared to the previous P. verdolaga annotation; moreover, the description of 65 novel theraphotoxins was possible. In the generation of data for non-model organisms, as well as in the search for novel peptides with biotechnological interest, it is highly recommended to employ at least two different transcriptome assemblers.

Using all gene families vastly expands data available for phylogenomic inference

Traditionally, single-copy orthologs have been the gold standard in phylogenomics. Most phylogenomic studies identify putative single-copy orthologs using clustering approaches and retain families with a single sequence per species. This limits the amount of data available by excluding larger families. Recent advances have suggested several ways to include data from larger families. For instance, tree-based decomposition methods facilitate the extraction of orthologs from large families. Additionally, several methods for species tree inference are robust to the inclusion of paralogs and could use all of the data from larger families. Here, we explore the effects of using all families for phylogenetic inference by examining relationships among 26 primate species in detail and by analyzing five additional data sets. We compare single-copy families, orthologs extracted using tree-based decomposition approaches, and all families with all data. We explore several species tree inference methods, finding that identical trees are returned across nearly all subsets of the data and methods for primates. The relationships among Platyrrhini remain contentious; however, the species tree inference method matters more than the subset of data used. Using data from larger gene families drastically increases the number of genes available and leads to consistent estimates of branch lengths, nodal certainty and concordance, and inferences of introgression in primates. For the other data sets, topological inferences are consistent whether single-copy families or orthologs extracted using decomposition approaches are analyzed. Using larger gene families is a promising approach to include more data in phylogenomics without sacrificing accuracy, at least when high-quality genomes are available.

Patterns in schizomid flagellum shape from elliptical Fourier analysis

The arachnid order Schizomida is a relatively understudied group of soil-dwelling predators found on all continents except Antarctica. While efforts to understand their biology are growing, there is still much to know about them. A curious aspect of their morphology is the male flagellum, a sexually dimorphic, tail-like structure which differs in shape across the order and functions in their courtship rituals. The flagellar shape is important for taxonomic classification, yet few efforts have been made to examine shape diversity across the group. Using elliptical Fourier analysis, a type of geometric morphometrics based on shape outline, we quantified shape differences across a combined nearly 550 outlines in the dorsal and lateral views, categorizing them based on genus, family, biogeographic realm, and habitat, with special emphasis on Caribbean and Cuban fauna. We tested for allometric relationships, differences in disparity based on locations and sizes in morphospace among these categories, and for clusters of shapes in morphospace. We found multiple differences in all categories despite apparent overlaps in morphospace, evolutionary allometry, and evidence for discrete clusters in some flagellum shapes. This study can serve as a foundation for further study on the evolution, diversification, and taxonomic utility of the male flagellum.

Phylogenomics of paleoendemic lampshade spiders (Araneae, Hypochilidae, Hypochilus), with the description of a new species from montane California

Hypochilus is a relictual lineage of Nearctic spiders distributed disjunctly across the United States in three montane regions (California, southern Rocky Mountains, southern Appalachia). Phylogenetic resolution of species relationships in Hypochilus has been challenging, and conserved morphology coupled with extreme genetic divergence has led to uncertain species limits in some complexes. Here, Hypochilus interspecies relationships have been reconstructed and cryptic speciation more critically evaluated using a combination of ultraconserved elements, mitochondrial CO1 by-catch, and morphology. Phylogenomic data strongly support the monophyly of regional clades and support a ((California, Appalachia), southern Rocky Mountains) topology. In Appalachia, five species are resolved as four lineages (H.thorelli Marx, 1888 and H.coylei Platnick, 1987 are clearly sister taxa), but the interrelationships of these four lineages remain unresolved. The Appalachian species H.pococki Platnick, 1987 is recovered as monophyletic but is highly genetically structured at the nuclear level. While algorithmic analyses of nuclear data indicate many species (e.g., all H.pococki populations as species), male morphology instead reveals striking stasis. Within the California clade, nuclear and mitochondrial lineages of H.petrunkevitchi Gertsch, 1958 correspond directly to drainage basins of the southern Sierra Nevada, with H.bernardino Catley, 1994 nested within H.petrunkevitchi and sister to the southernmost basin populations. Combining nuclear, mitochondrial, geographical, and morphological evidence a new species from the Tule River and Cedar Creek drainages is described, Hypochilusxomotesp. nov. We also emphasize the conservation issues that face several microendemic, habitat-specialized species in this remarkable genus.

Limitations of phylogenomic data can drive inferred speciation rate shifts

Biodiversity analyses of phylogenomic timetrees have produced many high-profile examples of shifts in the rate of speciation across the tree of life. Temporally correlated events in ecology, climate, and biogeography are frequently invoked to explain these rate shifts. In a re-examination of 15 genomic timetrees and 25 major published studies of the pattern of speciation through time, we observed an unexpected correlation between the timing of reported rate shifts and the information content of sequence alignments. Here, we show that the paucity of sequence variation and insufficient species sampling in phylogenomic datasets are the likely drivers of many inferred speciation rate shifts, rather than the proposed biological explanations. Therefore, data limitations can produce predictable but spurious signals of rate shifts even when speciation rates may be similar across taxa and time. Our results suggest that the reliable detection of speciation rate shifts requires the acquisition and assembly of long phylogenomic alignments with near-complete species sampling and accurate estimates of species richness for the clades of study.

Mitogenomics and hidden-trait models reveal the role of phoresy and host shifts in the diversification of parasitoid blister beetles (Coleoptera: Meloidae)

Changes in life history traits are often considered speciation triggers and can have dramatic effects on the evolutionary history of a lineage. Here, we examine the consequences of changes in two life history traits, host‐type and phoresy, in the hypermetamorphic blister beetles, Meloidae. Subfamilies Nemognathinae and Meloinae exhibit a complex life cycle involving multiple metamorphoses and parasitoidism. Most genera and tribes are bee‐parasitoids, and include phoretic or nonphoretic species, while two tribes feed on grasshopper eggs. These different life strategies are coupled with striking differences in species richness among clades. We generated a mitogenomic phylogeny for Nemognathinae and Meloinae, confirming the monophyly of these two clades, and used the dated phylogeny to explore the association between diversification rates and changes in host specificity and phoresy, using state‐dependent speciation and extinction (SSE) models that include the effect of hidden traits. To account for the low taxon sampling, we implemented a phylogenetic‐taxonomic approach based on birth‐death simulations, and used a Bayesian framework to integrate parameter and phylogenetic uncertainty. Results show that the ancestral hypermetamorphic Meloidae was a nonphoretic bee‐parasitoid, and that transitions towards a phoretic bee‐parasitoid and grasshopper parasitoidism occurred multiple times. Nonphoretic bee‐parasitoid lineages exhibit significantly higher relative extinction and lower diversification rates than phoretic bee‐and grasshopper‐parasitoids, but no significant differences were found between the latter two strategies. This suggests that Orthopteran host shifts and phoresy contributed jointly to the evolutionary success of the parasitoid meloidae. We also demonstrate that SSE models can be used to identify hidden traits coevolving with the focal trait in driving a lineage's diversification dynamics.