Integration of phylogenomics and molecular modeling reveals lineage-specific diversification of toxins in scorpions

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.

The remarkably enzyme-rich venom of the Big Bend Scorpion (Diplocentrus whitei)

Scorpion venoms have long been studied for their peptide discovery potential, with modern high-throughput venom-characterization techniques paving the way for the discovery of thousands of novel putative toxins. Research into these toxins has provided insight into the pathology and treatment of human diseases, even resulting in the development of one compound with Food and Drug Administration (FDA) approval. Although most of this research has focused on the toxins of scorpion species considered medically significant to humans, the venom of harmless scorpion species possess toxins that are homologous to those from medically significant species, indicating that harmless scorpion venoms may also serve as valuable sources of novel peptide variants. Furthermore, as harmless scorpions represent a vast majority of scorpion species diversity, and therefore venom toxin diversity, venoms from these species likely contain entirely new toxin classes. We sequenced the venom-gland transcriptome and venom proteome of two male Big Bend scorpions (Diplocentrus whitei), providing the first high-throughput venom characterization for a member of this genus. We identified a total of 82 toxins in the venom of D. whitei, 25 of which were identified in both the transcriptome and proteome, and 57 of which were only identified in the transcriptome. Furthermore, we identified a unique, enzyme-rich venom dominated by serine proteases and the first arylsulfatase B toxins identified in scorpions.

Let’s end taxonomic blank slates with molecular morphology

Many known evolutionary lineages have yet to be described formally due to a lack of traditional morphological characters. This is true for genetically distinctive groups within the amoeboid Placozoa animals, the protists in ponds, and the bacteria that cover nearly everything. These taxonomic tabula rasae, or blank slates, are problematic; without names, communication is hampered and other scientific progress is slowed. We suggest that the morphology of molecules be used to help alleviate this issue. Molecules, such as proteins, have structure. Proteins are even visualizable with X-ray crystallography, albeit more easily detected by and easier to work with using genomic sequencing. Given their structured nature, we believe they should not be considered as anything less than traditional morphology. Protein-coding gene content (presence/absence) can also be used easily with genomic sequences, and is a convenient binary character set. With molecular morphology, we believe that each taxonomic tabula rasa can be solved.

Evolution, expression patterns and distribution of novel ribbon worm predatory and defensive toxins

Ribbon worms are active predators that use an eversible proboscis to inject venom into their prey and defend themselves with toxic epidermal secretions. Previous work on nemertean venom has largely focused on just a few species and has not investigated the different predatory and defensive secretions in detail. Consequently, our understanding of the composition and evolution of ribbon worm venoms is still very limited. Here, we present a comparative study of nemertean venom combining RNA-seq differential gene expression analyses of venom-producing tissues, tandem mass spectrometry-based proteomics of toxic secretions, and mass spectrometry imaging of proboscis sections, to shed light onto the composition and evolution of predatory and defensive toxic secretions in Antarctonemertes valida. Our analyses reveal a wide diversity of putative defensive and predatory toxins with tissue-specific gene expression patterns and restricted distributions to the mucus and proboscis proteomes respectively, suggesting that ribbon worms produce distinct toxin cocktails for predation and defense. Our results also highlight the presence of numerous lineage-specific toxins, indicating that venom evolution is highly divergent across nemerteans, producing toxin cocktails that might be finely tuned to subdue different prey. Our data also suggest that the hoplonemertean proboscis is a highly specialized predatory organ that seems to be involved in a variety of biological functions besides predation, including secretion and sensory perception. Overall, our results advance our knowledge into the diversity and evolution of nemertean venoms and highlight the importance of combining different types of data to characterize toxin composition in understudied venomous organisms.

Burrowing into the forest: Phylogeny of the Asian forest scorpions (Scorpionidae: Heterometrinae) and the evolution of ecomorphotypes

Asian forest scorpions (Scorpionidae Latreille, 1802: Heterometrinae Simon, 1879) are distributed across South and Southeast Asia. All are fossorial, constructing burrows under stones or in open ground, in habitats differing in precipitation and vegetation cover, from rainforests and tropical deciduous forests to savanna and scrubland. The systematics of these scorpions has long been confused due to bad taxonomy and the absence of a phylogenetic framework. Although the monophyly of the group was previously confirmed as part of broader phylogenetic analyses based on exemplar species, the only quantitative analysis of species-level variation to date was based on overall similarity. This contribution presents the first species-level phylogenetic analysis of Asian Scorpionidae, based on 186 morphological characters and 4188 aligned base-pairs of DNA sequence data from two nuclear and three mitochondrial loci for 132 terminals including all 41 ingroup species and four outgroup species. Simultaneous analyses of the morphological and molecular datasets with parsimony, Maximum Likelihood and Bayesian Inference provided the framework for a revised classification presented elsewhere. In order to understand how adaptation following dispersal into new habitats has driven the morphological diversification of Asian forest scorpions, species were scored for 10 characters concerning morphology and burrow architecture, which contributed to an ensemble index of adaptation to habitat aridity. Species were classified into three ecomorphotypes based on the index, and ancestral state reconstruction of ecomorphotypes performed on the phylogeny. A pattern was recovered in which lineages and species occurring in different habitats on a continuum from wet (evergreen forest) to dry (savanna, scrubland) exhibited characters presumed to be adaptive and hence responsible for driving scorpion diversification.

Out of India, thrice: diversification of Asian forest scorpions reveals three colonizations of Southeast Asia

The 'Out of India' hypothesis is often invoked to explain patterns of distribution among Southeast Asian taxa. According to this hypothesis, Southeast Asian taxa originated in Gondwana, diverged from their Gondwanan relatives when the Indian subcontinent rifted from Gondwana in the Late Jurassic, and colonized Southeast Asia when it collided with Eurasia in the early Cenozoic. A growing body of evidence suggests these events were far more complex than previously understood, however. The first quantitative reconstruction of the biogeography of Asian forest scorpions (Scorpionidae Latreille, 1802: Heterometrinae Simon, 1879) is presented here. Divergence time estimation, ancestral range estimation, and diversification analyses are used to determine the origins, dispersal and diversification patterns of these scorpions, providing a timeline for their biogeographical history that can be summarized into four major events. (1) Heterometrinae diverged from other Scorpionidae on the African continent after the Indian subcontinent became separated in the Cretaceous. (2) Environmental stresses during the Cretaceous-Tertiary (KT) mass extinction caused range contraction, restricting one clade of Heterometrinae to refugia in southern India (the Western Ghats) and Sri Lanka (the Central Highlands). (3) Heterometrinae dispersed to Southeast Asia three times during India's collision with Eurasia, the first dispersal event occurring as the Indian subcontinent brushed up against the western side of Sumatra, and the other two events occurring as India moved closer to Eurasia. (4) Indian Heterometrinae, confined to southern India and Sri Lanka during the KT mass extinction, recolonized the Deccan Plateau and northern India, diversifying into new, more arid habitats after environmental conditions stabilized. These hypotheses, which are congruent with the geological literature and biogeographical analyses of other taxa from South and Southeast Asia, contribute to an improved understanding of the dispersal and diversification patterns of taxa in this biodiverse and geologically complex region.

Another one bites the dust: taxonomic sampling of a key genus in phylogenomic datasets reveals more non-monophyletic groups in traditional scorpion classification

Historically, morphological characters have been used to support the monophyly, composition, and phylogenetic relationships of scorpion families. Although recent phylogenomic analyses have recovered most of these traditional higher-level relationships as non-monophyletic, certain key taxa have yet to be sampled using a phylogenomic approach. Salient among these is the monotypic genus Caraboctonus Pocock,1893, the type species of the family Caraboctonidae Kraepelin, 1905. Here, we examined the putative monophyly and phylogenetic placement of this family, sampling the library of C. keyserlingi Pocock, 1893 using high throughput transcriptomic sequencing. Our phylogenomic analyses recovered Caraboctonidae as polyphyletic due to the distant placement of the genera Caraboctonus and Hadrurus Thorell, 1876. Caraboctonus was stably recovered as the sister-group of the monotypic family Superstitioniidae Stahnke, 1940, whereas Hadrurus formed an unstable relationship with Uroctonus Thorell, 1876and Belisarius Simon, 1879. Four-cluster likelihood mapping revealed that the instability inherent to the placement of Hadrurus, Uroctonus and Belisarius was attributable to significant gene tree conflict in the internodes corresponding to their divergences. To redress the polyphyly of Caraboctonidae, the following systematic actions have been taken: (1) the family Caraboctonidae has been delimited to consist of 23 species in the genera Caraboctonus and Hadruroides Pocock, 1893; (2) Caraboctonidae, previously included in the superfamily Iuroidea Thorell, 1876 or as incertae sedis, is transferred to the superfamily Caraboctonoidea (new rank); (3) the superfamily Hadruroidea (new rank) is established and the status of Hadrurinae Stahnke, 1973 is elevated to family (Hadruridae new status) including 9 species in the genera Hadrurus and Hoffmannihadrurus Fet & Soleglad, 2004 and (4) we treat Uroctonus and Belisarius as insertae sedis with respect to superfamilial placement. Our systematic actions engender the monophyly of both Iuroidea and Caraboctonidae. Future phylogenomic investigations should target similar taxon-poor and understudied lineages of potential phylogenetic significance, which are anticipated to reveal additional non-monophyletic groups.

Hadrurid Scorpion Toxins: Evolutionary Conservation and Selective Pressures

Scorpion toxins are thought to have originated from ancestral housekeeping genes that underwent diversification and neofunctionalization, as a result of positive selection. Our understanding of the evolutionary origin of these peptides is hindered by the patchiness of existing taxonomic sampling. While recent studies have shown phylogenetic inertia in some scorpion toxins at higher systematic levels, evolutionary dynamics of toxins among closely related taxa remain unexplored. In this study, we used new and previously published transcriptomic resources to assess evolutionary relationships of closely related scorpions from the family Hadruridae and their toxins. In addition, we surveyed the incidence of scorpine-like peptides (SLP, a type of potassium channel toxin), which were previously known from 21 scorpion species. We demonstrate that scorpine-like peptides exhibit gene duplications. Our molecular analyses demonstrate that only eight sites of two SLP copies found in scorpions are evolving under positive selection, with more sites evolving under negative selection, in contrast to previous findings. These results show evolutionary conservation in toxin diversity at shallow taxonomic scale.

Dissecting Toxicity: The Venom Gland Transcriptome and the Venom Proteome of the Highly Venomous Scorpion Centruroides limpidus (Karsch, 1879)

Venom glands and soluble venom from the Mexican scorpion Centruroides limpidus (Karsch, 1879) were used for transcriptomic and proteomic analyses, respectively. An RNA-seq was performed by high-throughput sequencing with the Illumina platform. Approximately 80 million reads were obtained and assembled into 198,662 putative transcripts, of which 11,058 were annotated by similarity to sequences from available databases. A total of 192 venom-related sequences were identified, including Na⁺ and K⁺ channel-acting toxins, enzymes, host defense peptides, and other venom components. The most diverse transcripts were those potentially coding for ion channel-acting toxins, mainly those active on Na⁺ channels (NaScTx). Sequences corresponding to β- scorpion toxins active of K⁺ channels (KScTx) and λ-KScTx are here reported for the first time for a scorpion of the genus Centruroides. Mass fingerprint corroborated that NaScTx are the most abundant components in this venom. Liquid chromatography coupled to mass spectometry (LC-MS/MS) allowed the identification of 46 peptides matching sequences encoded in the transcriptome, confirming their expression in the venom. This study corroborates that, in the venom of toxic buthid scorpions, the more abundant and diverse components are ion channel-acting toxins, mainly NaScTx, while they lack the HDP diversity previously demonstrated for the non-buthid scorpions. The highly abundant and diverse antareases explain the pancreatitis observed after envenomation by this species.

Phylogenomics facilitates stable scorpion systematics: Reassessing the relationships of Vaejovidae and a new higher-level classification of Scorpiones (Arachnida)

The Neartic family Vaejovidae (Scorpiones: Chactoidea) has long been treated as a diverse and systematically cohesive group of scorpions, but its monophyly and relationship to other scorpion families have historically been questioned. Morphological data have supported its monophyly and a variety of phylogenetic placements within the superfamily Chactoidea. Recent phylogenomic analyses have instead recovered vaejovids as polyphyletic (albeit with minimal taxonomic sampling) and Chactoidea as paraphyletic. Here, we reexamined the monophyly and phylogenetic placement of the family Vaejovidae, sampling 17 new vaejovid libraries using high throughput transcriptomic sequencing. Our phylogenomic analyses revealed a previous misplacement of Smeringurus mesaensis. Regardless, we recovered Vaejovidae as diphyletic due to the placement of the enigmatic genus Uroctonus. The remaining vaejovids formed a clade that was strongly supported as the sister group of the superfamily Scorpionoidea, a placement insensitive to matrix completeness or concatenation vs. species tree approaches to inferring the tree topology. Chactoidea was invariably recovered as a paraphyletic group due to the nested placement of Scorpionoidea. As first steps to resolving the paraphyly of Chactoidea, we take the following systematic actions: (1) we establish the superfamily Superstitionoidea (new superfamily) to accommodate Superstitioniidae; (2) we restore Vaejovoidea (status revalidated) as a valid superfamily that excludes Uroctonus; and (3) we treat the families Caraboctonidae, Troglotayosicidae, and the subfamily Uroctoninae as incertae sedis with respect to superfamilial placement. Our systematic actions thus establish the monophyly of the presently redefined Chactoidea and Vaejovoidea.