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.

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.

Taxonomic Sampling and Rare Genomic Changes Overcome Long-Branch Attraction in the Phylogenetic Placement of Pseudoscorpions.

Long-branch attraction is a systematic artifact that results in erroneous groupings of fast-evolving taxa. The combination of short, deep internodes in tandem with long-branch attraction artifacts has produced empirically intractable parts of the Tree of Life. One such group is the arthropod subphylum Chelicerata, whose backbone phylogeny has remained unstable despite improvements in phylogenetic methods and genome-scale data sets. Pseudoscorpion placement is particularly variable across data sets and analytical frameworks, with this group either clustering with other long-branch orders or with Arachnopulmonata (scorpions and tetrapulmonates). To surmount long-branch attraction, we investigated the effect of taxonomic sampling via sequential deletion of basally branching pseudoscorpion superfamilies, as well as varying gene occupancy thresholds in supermatrices. We show that concatenated supermatrices and coalescent-based summary species tree approaches support a sister group relationship of pseudoscorpions and scorpions, when more of the basally branching taxa are sampled. Matrix completeness had demonstrably less influence on tree topology. As an external arbiter of phylogenetic placement, we leveraged the recent discovery of an ancient genome duplication in the common ancestor of Arachnopulmonata as a litmus test for competing hypotheses of pseudoscorpion relationships. We generated a high-quality developmental transcriptome and the first genome for pseudoscorpions to assess the incidence of arachnopulmonate-specific duplications (e.g., homeobox genes and miRNAs). Our results support the inclusion of pseudoscorpions in Arachnopulmonata (new definition), as the sister group of scorpions. Panscorpiones (new name) is proposed for the clade uniting Scorpiones and Pseudoscorpiones.

Transcriptomic Analysis of Pseudoscorpion Venom Reveals a Unique Cocktail Dominated by Enzymes and Protease Inhibitors.

Transcriptomic and genomic analyses have illuminated the diversity of venoms in three of the four venomous arachnid orders (scorpions, spiders, and ticks). To date, no venom gland transcriptome analysis has been available for pseudoscorpions, the fourth venomous arachnid lineage. To redress this gap, we sequenced an mRNA library generated from the venom glands of the species Synsphyronus apimelus (Garypidae). High-throughput sequencing by the Illumina protocol, followed by de novo assembly, resulted in a total of 238,331 transcripts. From those, we annotated 131 transcripts, which code for putative peptides/proteins with similar sequences to previously reported venom components available from different arachnid species in protein databases. Transcripts putatively coding for enzymes showed the richest diversity, followed by other venom components such as peptidase inhibitors, cysteine-rich peptides, and thyroglobulin 1-like peptides. Only 11 transcripts were found that code for putatively low molecular mass spider toxins. This study constitutes the first report of the diversity of components within pseudoscorpion venom.