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.

Dual Functions of labial Resolve the Hox Logic of Chelicerate Head Segments.

Despite an abundance of gene expression surveys, comparatively little is known about Hox gene function in Chelicerata. Previous investigations of paralogs of labial (lab) and Deformed (Dfd) in a spider have shown that these play a role in tissue maintenance of the pedipalp segment (lab-1) and in patterning the first walking leg identity (Dfd-1), respectively. However, extrapolations of these data across chelicerates are hindered by the existence of duplicated Hox genes in arachnopulmonates (e.g., spiders and scorpions), which have resulted from an ancient whole genome duplication (WGD) event. Here, we investigated the function of the single-copy ortholog of lab in the harvestman Phalangium opilio, an exemplar of a lineage that was not subject to this WGD. Embryonic RNA interference against lab resulted in two classes of phenotypes: homeotic transformations of pedipalps to chelicerae, as well as reduction and fusion of the pedipalp and leg 1 segments. To test for combinatorial function, we performed a double knockdown of lab and Dfd, which resulted in a homeotic transformation of both pedipalps and the first walking legs into cheliceral identity, whereas the second walking leg is transformed into a pedipalpal identity. Taken together, these results elucidate a model for the Hox logic of head segments in Chelicerata. To substantiate the validity of this model, we performed expression surveys for lab and Dfd paralogs in scorpions and horseshoe crabs. We show that repetition of morphologically similar appendages is correlated with uniform expression levels of the Hox genes lab and Dfd, irrespective of the number of gene copies.

A plurality of morphological characters need not equate with phylogenetic accuracy: A rare genomic change refutes the placement of Solifugae and Pseudoscorpiones in Haplocnemata.

Recent advances in higher-level invertebrate phylogeny have leveraged shared features of genomic architecture to resolve contentious nodes across the tree of life. Yet, the interordinal relationships within Chelicerata have remained recalcitrant given competing topologies in recent molecular analyses. As such, relationships between topologically unstable orders remain supported primarily by morphological cladistic analyses. Solifugae, one such unstable chelicerate order, has long been thought to be the sister group of Pseudoscorpiones, forming the clade Haplocnemata, on the basis of eight putative morphological synapomorphies. The discovery, however, of a shared whole genome duplication placing Pseudoscorpiones in Arachnopulmonata provides the opportunity for a simple litmus test evaluating the validity of Haplocnemata. Here, we present the first developmental transcriptome of a solifuge (Titanopuga salinarum) and survey copy numbers of the homeobox genes for evidence of systemic duplication. We find that over 70% of the identified homeobox genes in T. salinarum are retained in a single copy, while representatives of the arachnopulmonates retain orthologs of those genes as two or more copies. Our results refute the placement of Solifugae in Haplocnemata. Subsequent reevaluation of putative interordinal morphological synapomorphies among chelicerates reveals a high incidence of homoplasy, reversals, and inaccurate coding within Haplocnemata and other small clades, as well as Arachnida more broadly, suggesting existing morphological character matrices are insufficient to resolve chelicerate phylogeny.

Under the hood: Phylogenomics of hooded tick spiders (Arachnida, Ricinulei) uncovers discordance between morphology and molecules.

Ricinulei or hooded tick-spiders are a cryptic and ancient group of arachnids. The order consists of around 100 highly endemic extant species restricted to the Afrotropics and the Neotropics along with 22 fossil species. Their antiquity and low vagility make them an excellent group with which to interrogate biogeographic questions. To date, only four molecular analyses have been conducted on the group and they failed to resolve the relationships of the main lineages and even recovering the non-monophyly of the three genera. These studies were limited to a few Sanger loci or phylogenomic analyses with at most seven ingroup samples. To increase phylogenetic resolution in this little-understood and poorly studied group, we present the most comprehensive phylogenomic study of Ricinulei to date leveraging the Arachnida ultra-conserved element probe set. With a data set of 473 loci across 96 ingroup samples, analyses resolved a monophyletic Neotropical clade consisting of four main lineages. Two of them correspond to the current genera Cryptocellus and Pseudocellus while topology testing revealed one lineage to likely be a phylogenetic reconstruction artefact. The fourth lineage, restricted to Northwestern, Andean South America, is consistent with the Cryptocellus magnus group, likely corresponding to the historical genus Heteroricinoides. Since we did not sample the type species for this old genus, we do not formally re-erect Heteroricinoides but our data suggest the need for a thorough morphological re-examination of Neotropical Ricinulei.

Machine learning approaches delimit cryptic taxa in a previously intractable species complex.

Cryptic species are not diagnosable via morphological criteria, but can be detected through analysis of DNA sequences. A number of methods have been developed for identifying species based on genetic data; however, these methods are prone to over-splitting taxa with extreme population structure, such as dispersal-limited organisms. Machine learning methodologies have the potential to overcome this challenge. Here, we apply such approaches, using a large dataset generated through hybrid target enrichment of ultraconserved elements (UCEs). Our study taxon is the Aoraki denticulata species complex, a lineage of extremely low-dispersal arachnids endemic to the South Island of Aotearoa New Zealand. This group of mite harvesters has been the subject of previous species delimitation studies using smaller datasets generated through Sanger sequencing and analytical approaches that rely on multispecies coalescent models and barcoding gap discovery. Those analyses yielded a number of putative cryptic species that seems unrealistic and extreme, based on what we know about species' geographic ranges and genetic diversity in non-cryptic mite harvesters. We find that machine learning approaches, on the other hand, identify cryptic species with geographic ranges that are similar to those seen in other morphologically diagnosable mite harvesters in Aotearoa New Zealand's South Island. We performed both unsupervised and supervised machine learning analyses, the latter with training data drawn either from animals broadly (vagile and non-vagile) or from a custom training dataset from dispersal-limited harvesters. We conclude that applying machine learning approaches to the analysis of UCE-derived genetic data is an effective method for delimiting species in complexes of low-vagility cryptic species, and that the incorporation of training data from biologically relevant analogues can be critically informative.

Phylogenomic data reveal three new families of poorly studied Solifugae (camel spiders).

The systematics of the arachnid order Solifugae have been an enigma, owing to challenges in interpreting morphology, a paucity of molecular phylogenetic studies sampling across the group, and a dearth of taxonomic attention for many lineages. Recent work has suggested that solifuge families largely exhibit contiguous distributions and reflect patterns of vicariance, with the exception of three families: Melanoblossidae, Daesiidae and Gylippidae. Morphological studies have cast doubt on their existing circumscriptions and the present composition of these taxa renders their distributions as disjunct. We leveraged ultraconserved elements (UCEs) to test the phylogenetic placement of three key lineages of Solifugae that cause these anomalous distributions: Dinorhax rostrumpsittaci (putative melanoblossid), Namibesia (putative daesiid), and Trichotoma (putative gylippid). Phylogenetic placement of these three genera based on UCEs rendered the families that harbor them as para- or polyphyletic, recovering instead relationships that better accord with a biogeographic history driven by vicariance. Toward a stable and phylogenetically informed classification of Solifugae, we establish three new families, Dinorhaxidae new rank, Namibesiidae new rank and Lipophagidae new rank.

The kinematics of amblypygid (Arachnida) pedipalps during predation: extreme elongation in raptorial appendages does not result in a proportionate increase in reach and closing speed.

The link between form and function is key to understanding the evolution of unique and/or extreme morphologies. Amblypygids, or whip spiders, are arachnids that often have highly elongated spined pedipalps. These limbs are used to strike at, and secure, prey before processing by the chelicerae. Amblypygi pedipalps are multifunctional, however, being used in courtship and contest, and vary greatly in form between species. Increased pedipalp length may improve performance during prey capture, but length could also be influenced by factors including territorial contest and sexual selection. Here, for the first time, we used high-speed videography and manual tracking to investigate kinematic differences in prey capture between amblypygid species. Across six morphologically diverse species, spanning four genera and two families, we created a total dataset of 86 trials (9-20 per species). Prey capture kinematics varied considerably between species, with differences being expressed in pedipalp joint angle ranges. In particular, maximum reach ratio did not remain constant with total pedipalp length, as geometric scaling would predict, but decreased with longer pedipalps. This suggests that taxa with the most elongated pedipalps do not deploy their potential length advantage to proportionally increase reach. Therefore, a simple mechanical explanation of increased reach does not sufficiently explain pedipalp elongation. We propose other factors to help explain this phenomenon, such as social interactions or sexual selection, which would produce an evolutionary trade-off in pedipalp length between prey capture performance and other behavioural and/or anatomical pressures.

Habituation to a predatory stimulus in a harvester (Arachnida, Opiliones).

Several studies have investigated habituation in a defensive context, but few have addressed responses to dangerous stimuli. In such cases, animals should not habituate since this could cost their lives. Here we have stimulated individuals of the harvester Mischonyx squalidus with a predatory stimulus (squeezing with tweezers) in repeated trials within and between days, and measured the occurrence and magnitude of nipping, a defensive behavior. Contrary to our expectations, they did habituate to this stimulus. The probability and magnitude of response declined over trials during each of three days of testing in a typical habituation pattern. During the trials we also observed other defensive behaviors. We discuss our results mainly considering alternative defensive responses. Our data show that we lack information on (1) the role played by the ambiguity of stimuli, (2) the role played by subsequent stimuli and (3) the importance of the array of defensive behaviors of a species in understanding habituation. Although ubiquitous across animals and therefore expected, habituation is described for the first time in the order Opiliones.

Functional traits of predators and decomposer prey determine context dependency in trophic control over ecosystems.

Research Highlight: Piccoli, G. C. d. O., Antiqueira, P. A. P., Srivastava, D. S., & Romero, G. Q. (2024). Trophic cascades within and across ecosystems: The role of anti-predatory defences, predator type and detritus quality. Journal of Animal Ecology, 00, 1-14. https://doi.org/10.1111/1365-2656.14063. Ecosystem functioning is controlled by the interplay between bottom-up supply of limiting nutrients and top-down animal feedback effects. However, the degree of animal versus nutrient control is context-dependent. A key challenge lies in characterizing this context dependency which is hypothesized to depend on differences in animal functional traits. Reporting on an important experiment, Piccoli et al. (2014) evaluate how interactions among functionally different predators and decomposer prey create context dependency in top-down control of a model system-tropical bromeliad tank ecosystems. Bromeliad plants hold water in their tanks supporting microcosm ecosystems containing terrestrial and aquatic insect larvae and arachnids. The ecosystems are supported by nutrients in plant litter that rains down from forest canopies into the tanks. Nutrients are released after litter is decomposed by a functionally diverse community of larval insect decomposers that differ in feeding mode and antipredator defence strategy. This decomposer community is preyed upon by an exclusively narrowly ranging aquatic insect larval predator and widely ranging spider predator that crosses between the aquatic and surrounding terrestrial ecosystems. Experimental manipulation of the animal community to test for the degree of control by predators mediated by the functionally diverse prey community included four treatments: (i) a control with the detritivores composing different function groups but without predators, (ii) the cross-ecosystem spider predator added, (iii) the purely aquatic damselfly larvae predator added and (iv) both predator types added to capture their interacting effect on ecosystem function (decomposition, nutrient release, and plant growth). Notably, the study resolved the causal pathways and strengths of direct and indirect control using structural equation modelling. These findings reveal how context dependency arises due to different capacities of the predators alone and together to overcome prey defences and control their abundances, with attendant cascading effects that diminished as well as enhanced decomposition and nutrient release to support bromeliad plant production. The study reveals that predators have a decided, albeit qualitatively and quantitatively different, hand in shaping the degree of bottom-up control through feedback effect on the release of limiting nutrients. This ground-breaking study provides a way forward in understanding the mechanisms determining context dependency in the control over ecosystem functioning.

Comprehensive Species Sampling and Sophisticated Algorithmic Approaches Refute the Monophyly of Arachnida.

Deciphering the evolutionary relationships of Chelicerata (arachnids, horseshoe crabs, and allied taxa) has proven notoriously difficult, due to their ancient rapid radiation and the incidence of elevated evolutionary rates in several lineages. Although conflicting hypotheses prevail in morphological and molecular data sets alike, the monophyly of Arachnida is nearly universally accepted, despite historical lack of support in molecular data sets. Some phylotranscriptomic analyses have recovered arachnid monophyly, but these did not sample all living orders, whereas analyses including all orders have failed to recover Arachnida. To understand this conflict, we assembled a data set of 506 high-quality genomes and transcriptomes, sampling all living orders of Chelicerata with high occupancy and rigorous approaches to orthology inference. Our analyses consistently recovered the nested placement of horseshoe crabs within a paraphyletic Arachnida. This result was insensitive to variation in evolutionary rates of genes, complexity of the substitution models, and alternative algorithmic approaches to species tree inference. Investigation of sources of systematic bias showed that genes and sites that recover arachnid monophyly are enriched in noise and exhibit low information content. To test the impact of morphological data, we generated a 514-taxon morphological data matrix of extant and fossil Chelicerata, analyzed in tandem with the molecular matrix. Combined analyses recovered the clade Merostomata (the marine orders Xiphosura, Eurypterida, and Chasmataspidida), but merostomates appeared nested within Arachnida. Our results suggest that morphological convergence resulting from adaptations to life in terrestrial habitats has driven the historical perception of arachnid monophyly, paralleling the history of numerous other invertebrate terrestrial groups.

The Geography of Sexual Conflict: A Synthetic Review.

AbstractSexual conflict is a mechanism of selection driven by the divergent fitness interests between females and males. This disagreement can be great enough to promote antagonistic/defensive traits and behaviors. Although the existence of sexual conflict has been identified in many species, less research has explored the conditions that initially promote sexual conflict in animal mating systems. In previous work in Opiliones, we observed that morphological traits associated with sexual conflict occurred only in species from northern localities. We hypothesized that by shortening and compartmentalizing time periods optimal for reproduction, seasonality represents a geographic condition sufficient to promote sexual conflict. We conducted a systematic review of the literature on reproductive traits and behaviors. Using standardized criteria, we reviewed publications to identify whether subjects occurred in a temperate (high-seasonality) or tropical (low-seasonality) biome. After identifying and adjusting for a publication bias toward temperate research, we identified no significant difference in the strength of sexual conflict between temperate and tropical study systems. A comparison between the distribution of taxa studied in sexual conflict articles and articles focused on general biodiversity indicates that species with conflict-based mating systems more accurately represent the distribution of terrestrial animal species. These findings contribute to ongoing efforts to characterize the origins of sexual conflict as well as life history traits that covary with sexual conflict.

Too early for the ferry: The biogeographic history of the Assamiidae of southeast Asia (Chelicerata: Opiliones, Laniatores).

Opiliones (harvestmen) have come to be regarded as an abundant source of model groups for study of historical biogeography, due to their ancient age, poor dispersal capability, and high fidelity to biogeographic terranes. One of the least understood harvestman groups is the Paleotropical Assamiidae, one of the more diverse families of Opiliones. Due to a labyrinthine taxonomy, poorly established generic and subfamilial boundaries, and the lack of taxonomic keys for the group, few efforts have been undertaken to decipher relationships within this arachnid lineage. Neither the monophyly of the family, nor its exact placement in the harvestman phylogeny, have been established. Here, we assessed the internal phylogeny of Assamiidae using a ten-locus Sanger dataset, sampling key lineages putatively ascribed to this family for five of the ten markers. Our analyses recovered Assamiidae as a monophyletic group, in a clade with the primarily Afrotropical Pyramidopidae and the southeast Asian Beloniscidae. Internal relationships of assamiids disfavored the systematic validity of subfamilies, with biogeography reflecting much better phylogenetic structure than the existing higher-level taxonomy. To assess whether the Asian assamiids came to occupy Indo-Pacific terranes via rafting on the Indian subcontinent, we performed divergence dating to infer the age of the family. Our results show that Indo-Pacific clades are ancient, originating well before the Cretaceous and therefore predate a vicariant mechanism commonly encountered for Paleotropical taxa.

Stabilized Morphological Evolution of Spiders Despite Mosaic Changes in Foraging Ecology.

A prominent question in animal research is how the evolution of morphology and ecology interacts in the generation of phenotypic diversity. Spiders are some of the most abundant arthropod predators in terrestrial ecosystems and exhibit a diversity of foraging styles. It remains unclear how spider body size and proportions relate to foraging style, and if the use of webs as prey capture devices correlates with changes in body characteristics. Here, we present the most extensive data set to date of morphometric and ecological traits in spiders. We used this data set to estimate the change in spider body sizes and shapes over deep time and to test if and how spider phenotypes are correlated with their behavioral ecology. We found that phylogenetic variation of most traits best fitted an Ornstein-Uhlenbeck model, which is a model of stabilizing selection. A prominent exception was body length, whose evolutionary dynamics were best explained with a Brownian Motion (free trait diffusion) model. This was most expressed in the araneoid clade (ecribellate orb-weaving spiders and allies) that showed bimodal trends toward either miniaturization or gigantism. Only few traits differed significantly between ecological guilds, most prominently leg length and thickness, and although a multivariate framework found general differences in traits among ecological guilds, it was not possible to unequivocally associate a set of morphometric traits with the relative ecological mode. Long, thin legs have often evolved with aerial webs and a hanging (suspended) locomotion style, but this trend is not general. Eye size and fang length did not differ between ecological guilds, rejecting the hypothesis that webs reduce the need for visual cue recognition and prey immobilization. For the inference of the ecology of species with unknown behaviors, we propose not to use morphometric traits, but rather consult (micro-)morphological characters, such as the presence of certain podal structures. These results suggest that, in contrast to insects, the evolution of body proportions in spiders is unusually stabilized and ecological adaptations are dominantly realized by behavioral traits and extended phenotypes in this group of predators. This work demonstrates the power of combining recent advances in phylogenomics with trait-based approaches to better understand global functional diversity patterns through space and time. [Animal architecture; Arachnida; Araneae; extended phenotype; functional traits; macroevolution; stabilizing selection.].

Comparative biology of spatial navigation in three arachnid orders (Amblypygi, Araneae, and Scorpiones).

From both comparative biology and translational research perspectives, there is escalating interest in understanding how animals navigate their environments. Considerable work is being directed towards understanding the sensory transduction and neural processing of environmental stimuli that guide animals to, for example, food and shelter. While much has been learned about the spatial orientation behavior, sensory cues, and neurophysiology of champion navigators such as bees and ants, many other, often overlooked animal species possess extraordinary sensory and spatial capabilities that can broaden our understanding of the behavioral and neural mechanisms of animal navigation. For example, arachnids are predators that often return to retreats after hunting excursions. Many of these arachnid central-place foragers are large and highly conducive to scientific investigation. In this review we highlight research on three orders within the Class Arachnida: Amblypygi (whip spiders), Araneae (spiders), and Scorpiones (scorpions). For each, we describe (I) their natural history and spatial navigation, (II) how they sense the world, (III) what information they use to navigate, and (IV) how they process information for navigation. We discuss similarities and differences among the groups and highlight potential avenues for future research.

Congruence between ultraconserved element-based matrices and phylotranscriptomic datasets in the scorpion Tree of Life.

Scorpions are ancient and historically renowned for their potent venom. Traditionally, the systematics of this group of arthropods was supported by morphological characters, until recent phylogenomic analyses (using RNAseq data) revealed most of the higher-level taxa to be non-monophyletic. While these phylogenomic hypotheses are stable for almost all lineages, some nodes have been hard to resolve due to minimal taxonomic sampling (e.g. family Chactidae). In the same line, it has been shown that some nodes in the Arachnid Tree of Life show disagreement between hypotheses generated using transcritptomes and other genomic sources such as the ultraconserved elements (UCEs). Here, we compared the phylogenetic signal of transcriptomes vs. UCEs by retrieving UCEs from new and previously published scorpion transcriptomes and genomes, and reconstructed phylogenies using both datasets independently. We reexamined the monophyly and phylogenetic placement of Chactidae, sampling an additional chactid species using both datasets. Our results showed that both sets of genome-scale datasets recovered highly similar topologies, with Chactidae rendered paraphyletic owing to the placement of Nullibrotheas allenii. As a first step toward redressing the systematics of Chactidae, we establish the family Anuroctonidae (new family) to accommodate the genus Anuroctonus.

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.

Retrospective analysis of histologic lesions in captive arachnids.

Invertebrates, including arachnids, are a common taxon in zoological collections. Invertebrate medicine and pathology are emerging subspecialties, but there is limited reference material or published resources describing histologic lesions in arachnids. Histopathology of 26 captive arachnids (20 spiders and 6 scorpions) from institutional collections was reviewed. Most animals were found dead with limited clinical signs. Tissues evaluated included body wall (cuticle and epidermis), skeletal muscle, book lungs, digestive tract (pharynx, esophagus, sucking stomach, midgut tube, midgut diverticula, and stercoral pocket), central and peripheral nervous system, heart, hemolymph vessels and sinuses, Malpighian tubules, coxal glands, and gonads. Inflammation was frequent (24/26, 92%), and seen in multiple organs (18/24, 75%) with the midgut diverticulum most commonly affected (14/24, 58%) followed by the book lungs (13/24 arachnids, 54%), and body wall (8/24 arachnids, 33%). Inflammation comprised hemocyte accumulation, hemocytic coagula, melanization, and nodulation. Infectious agents, including bacteria (11/26, 42%), fungi (10/26, 38%), and parasites (2/26, 8%), were seen within inflammatory aggregates. Coinfection with multiple infectious agents was common (6/24, 25%). No etiologic agent was identified in 7/24 (29%) cases with inflammatory lesions. Lesions suggestive of decreased nutritional status or increased metabolic rate included midgut diverticula atrophy in 11/26 (42%) animals and skeletal muscle atrophy in 6/26 (23%) animals. Atrophic lesions were seen in combination with infection (8/11, 73%), pregnancy (2/11, 18%), male sex (2/11, 18%), or without other lesions (1/11, 9%). Other suspected contributors to death included dysecdysis-associated trauma (2/26, 8%) and uterine intussusception (1/26, 4%). No animals had neoplasia.

The Evolution of Sox Gene Repertoires and Regulation of Segmentation in Arachnids.

The Sox family of transcription factors regulates many processes during metazoan development, including stem cell maintenance and nervous system specification. Characterizing the repertoires and roles of these genes can therefore provide important insights into animal evolution and development. We further characterized the Sox repertoires of several arachnid species with and without an ancestral whole-genome duplication and compared their expression between the spider Parasteatoda tepidariorum and the harvestman Phalangium opilio. We found that most Sox families have been retained as ohnologs after whole-genome duplication and evidence for potential subfunctionalization and/or neofunctionalization events. Our results also suggest that Sox21b-1 likely regulated segmentation ancestrally in arachnids, playing a similar role to the closely related SoxB gene, Dichaete, in insects. We previously showed that Sox21b-1 is required for the simultaneous formation of prosomal segments and sequential addition of opisthosomal segments in P. tepidariorum. We studied the expression and function of Sox21b-1 further in this spider and found that although this gene regulates the generation of both prosomal and opisthosomal segments, it plays different roles in the formation of these tagmata reflecting their contrasting modes of segmentation and deployment of gene regulatory networks with different architectures.

Evidence of positive selection on six spider developmental genes.

Spiders constitute more than 49,000 described species distributed all over the world, and all ecological environments. Their order, Araneae, is defined by a set of characteristics with no parallel among their arachnid counterparts (e.g., spinnerets, silk glands, chelicerae that inoculate venom, among others). Changes in developmental pathways often underlie the evolution of morphological synapomorphies, and as such spiders are a promising model to study the role of developmental genes in the origin of evolutionary novelties. With that in mind, we investigated changes in the evolutionary regime of a set of six developmental genes, using spiders as our model. The genes were mainly chosen for their roles in spinneret ontogeny, yet they are pleiotropic, and it is likely that the origins of other unique morphological phenotypes are also linked to changes in their sequences. Our results indicate no great differences in the selective pressures on those genes when comparing spiders to other arachnids, but a few site-specific positive selection evidence were found in the Araneae lineage. These findings lead us to new insights on spider evolution that are to be further tested.

How many long branch orders occur in Chelicerata? Opposing effects of Palpigradi and Opilioacariformes on phylogenetic stability.

Excepting a handful of nodes, phylogenetic relationships between chelicerate orders remain poorly resolved, due to both the incidence of long branch attraction artifacts and the limited sampling of key lineages. It has recently been shown that increasing representation of basal nodes plays an outsized role in resolving the higher-level placement of long-branch chelicerate orders. Two lineages have been consistently undersampled in chelicerate phylogeny. First, sampling of the miniaturized order Palpigradi has been restricted to a fragmentary transcriptome of a single species. Second, sampling of Opilioacariformes, a rarely encountered and key group of Parasitiformes, has been restricted to a single exemplar. These two lineages exhibit dissimilar properties with respect to branch length; Opilioacariformes shows relatively low evolutionary rate compared to other Parasitiformes, whereas Palpigradi possibly acts as another long-branch order (an effect that may be conflated with the degree of missing data). To assess these properties and their effects on tree stability, we constructed a phylogenomic dataset of Chelicerata wherein both lineages were sampled with three terminals, increasing the representation of these taxa per locus. We examined the effect of subsampling phylogenomic matrices using (1) taxon occupancy, (2) evolutionary rate, and (3) a principal components-based approach. We further explored the impact of taxon deletion experiments that mitigate the effect of long branches. Here, we show that Palpigradi constitutes a fourth long-branch chelicerate order (together with Acariformes, Parasitiformes, and Pseudoscorpiones), which further destabilizes the chelicerate backbone topology. By contrast, the slow-evolving Opilioacariformes were consistently recovered within Parasitiformes, with certain subsampling practices recovering their placement as the sister group to the remaining Parasitiformes. Whereas the inclusion of Opilioacariformes always resulted in the non-monophyly of Acari with support, deletion of Opilioacariformes from datasets consistently incurred the monophyly of Acari except in matrices constructed on the basis of evolutionary rate. Our results strongly suggest that Acari is an artifact of long- branch attraction.