Increasing species sampling in chelicerate genomic-scale datasets provides support for monophyly of Acari and Arachnida

Major vault protein is part of an extracellular cement material in the Atlantic salmon louse (Lepeophtheirus salmonis)

Major vault protein (MVP) is the main component of the vault complex, which is a highly conserved ribonucleoprotein complex found in most eukaryotic organisms. MVP or vaults have previously been found to be overexpressed in multidrug-resistant cancer cells and implicated in various cellular processes such as cell signaling and innate immunity. The precise function of MVP is, however, poorly understood and its expression and probable function in lower eukaryotes are not well characterized. In this study, we report that the Atlantic salmon louse expresses three full-length MVP paralogues (LsMVP1-3). Furthermore, we extended our search and identified MVP orthologues in several other ecdysozoan species. LsMVPs were shown to be expressed in various tissues at both transcript and protein levels. In addition, evidence for LsMVP to assemble into vaults was demonstrated by performing differential centrifugation. LsMVP was found to be highly expressed in cement, an extracellular material produced by a pair of cement glands in the adult female salmon louse. Cement is important for the formation of egg strings that serve as protective coats for developing embryos. Our results imply a possible novel function of LsMVP as a secretory cement protein. LsMVP may play a role in structural or reproductive functions, although this has to be further investigated.

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.

Phylogenomic analysis and molecular identification of true fruit flies

The family Tephritidae in the order Diptera, known as true fruit flies, are agriculturally important insect pests. However, the phylogenetic relationships of true fruit flies, remain controversial. Moreover, rapid identification of important invasive true fruit flies is essential for plant quarantine but is still challenging. To this end, we sequenced the genome of 16 true fruit fly species at coverage of 47-228×. Together with the previously reported genomes of nine species, we reconstructed phylogenetic trees of the Tephritidae using benchmarking universal single-copy ortholog (BUSCO), ultraconserved element (UCE) and anchored hybrid enrichment (AHE) gene sets, respectively. The resulting trees of 50% taxon-occupancy dataset for each marker type were generally congruent at 88% nodes for both concatenation and coalescent analyses. At the subfamily level, both Dacinae and Trypetinae are monophyletic. At the species level, Bactrocera dorsalis is more closely related to Bactrocera latifrons than Bactrocera tryoni. This is inconsistent with previous conclusions based on mitochondrial genes but consistent with recent studies based on nuclear data. By analyzing these genome data, we screened ten pairs of species-specific primers for molecular identification of ten invasive fruit flies, which PCR validated. In summary, our work provides draft genome data of 16 true fruit fly species, addressing the long-standing taxonomic controversies and providing species-specific primers for molecular identification of invasive fruit flies.

Lower Ordovician synziphosurine reveals early euchelicerate diversity and evolution

Euchelicerata is a clade of arthropods comprising horseshoe crabs, scorpions, spiders, mites and ticks, as well as the extinct eurypterids (sea scorpions) and chasmataspidids. The understanding of the ground plans and relationships between these crown-group euchelicerates has benefited from the discovery of numerous fossils. However, little is known regarding the origin and early evolution of the euchelicerate body plan because the relationships between their Cambrian sister taxa and synziphosurines, a group of Silurian to Carboniferous stem euchelicerates with chelicerae and an unfused opisthosoma, remain poorly understood owing to the scarce fossil record of appendages. Here we describe a synziphosurine from the Lower Ordovician (ca. 478 Ma) Fezouata Shale of Morocco. This species possesses five biramous appendages with stenopodous exopods bearing setae in the prosoma and a fully expressed first tergite in the opisthosoma illuminating the ancestral anatomy of the group. Phylogenetic analyses recover this fossil as a member of the stem euchelicerate family Offacolidae, which is characterized by biramous prosomal appendages. Moreover, it also shares anatomical features with the Cambrian euarthropod Habelia optata, filling the anatomical gap between euchelicerates and Cambrian stem taxa, while also contributing to our understanding of the evolution of euchelicerate uniramous prosomal appendages and tagmosis.

Unveiling the Genetic Blueprint of a Desert Scorpion: A Chromosome-level Genome of Hadrurus arizonensis Provides the First Reference for Parvorder Iurida

Over 400 million years old, scorpions represent an ancient group of arachnids and one of the first animals to adapt to life on land. Presently, the lack of available genomes within scorpions hinders research on their evolution. This study leverages ultralong nanopore sequencing and Pore-C to generate the first chromosome-level assembly and annotation for the desert hairy scorpion, Hadrurus arizonensis. The assembled genome is 2.23 Gb in size with an N50 of 280 Mb. Pore-C scaffolding reoriented 99.6% of bases into nine chromosomes and BUSCO identified 998 (98.6%) complete arthropod single copy orthologs. Repetitive elements represent 54.69% of the assembled bases, including 872,874 (29.39%) LINE elements. A total of 18,996 protein-coding genes and 75,256 transcripts were predicted, and extracted protein sequences yielded a BUSCO score of 97.2%. This is the first genome assembled and annotated within the family Hadruridae, representing a crucial resource for closing gaps in genomic knowledge of scorpions, resolving arachnid phylogeny, and advancing studies in comparative and functional genomics.

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.

Vestigial organs alter fossil placements in an ancient group of terrestrial chelicerates

Vestigial organs provide a link between ancient and modern traits and therefore have great potential to resolve the phylogeny of contentious fossils that bear features not seen in extant species. Here we show that extant daddy-longlegs (Arachnida, Opiliones), a group once thought to possess only one pair of eyes, in fact additionally retain a pair of vestigial median eyes and a pair of vestigial lateral eyes. Neuroanatomical gene expression surveys of eye-patterning transcription factors, opsins, and other structural proteins in the daddy-longlegs Phalangium opilio show that the vestigial median and lateral eyes innervate regions of the brain positionally homologous to the median and lateral eye neuropils, respectively, of chelicerate groups like spiders and horseshoe crabs. Gene silencing of eyes absent shows that the vestigial eyes are under the control of the retinal determination gene network. Gene silencing of dachshund disrupts the lateral eyes, but not the median eyes, paralleling loss-of-function phenotypes in insect models. The existence of lateral eyes in extant daddy-longlegs bears upon the placement of the oldest harvestmen fossils, a putative stem group that possessed both a pair of median eyes and a pair of lateral eyes. Phylogenetic analysis of harvestman relationships with an updated understanding of lateral eye incidence resolved the four-eyed fossil group as a member of the extant daddy-longlegs suborder, which in turn resulted in older estimated ages of harvestman diversification. This work underscores that developmental vestiges in extant taxa can influence our understanding of character evolution, placement of fossils, and inference of divergence times.

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.

Diversity and Distribution of Mites (ACARI) Revealed by Contamination Survey in Public Genomic Databases

Acari (mites and ticks) are a biodiverse group of microarthropods within the Arachnida. Because of their diminutive size, mites are often overlooked. We hypothesized that mites, like other closely related microorganisms, could also contaminate public genomic database. Here, using a strategy based on DNA barcodes previously reported, we scanned contaminations related to mites (Acari, exclusive of Ixodida) in Genbank WGS/TSA database. In 22,114 assemblies (17,845 animal and 4269 plant projects), 1717 contigs in 681 assemblies (3.1%) were detected as mite contaminations. Additional taxonomic analysis showed the following: (1) most of the contaminants (1445/1717) were from the specimens of Magnoliopsida, Insecta and Pinopsida; (2) the contamination rates were higher in plant or TSA projects; (3) mite distribution among different classes of hosts varied considerably. Additional phylogenetic analysis of these contaminated contigs further revealed complicated mite-host associations. Overall, we conducted a first systemic survey and analysis of mite contaminations in public genomic database, and these DNA barcode related mite contigs will provide a valuable resource of information for understanding the diversity and phylogeny of mites.

Evolution and phylogenetic diversity of the aquaporin gene family in arachnids

Water flux across cells predominantly occurs through the pore formed by the aquaporin channels. Since water balance is one of the most important challenges to terrestrial animals, aquaporin evolution and diversity is known to play roles in animal terrestrialisation. Arachnids (Arthropoda: Chelicerata: Arachnida) are the second most diverse group and represent the pioneer land colonists in animals; however, there remains no thorough investigation on aquaporin evolution and diversity in this evolutionarily important lineage. Here we reported a phylogenetic study of aquaporin evolution and diversity using genomic data from 116 arachnid species covering almost all (15/16) extant orders. A previously unrecognised subfamily related to aquaporin-4 (i.e. Aqp4-like subfamily) via phylogenetic analysis was identified, suggesting certain underestimate of the arachnid aquaporin diversity in earlier studies probably due to limited taxonomic sampling. Further analysis indicates that this subfamily emerged deep within the life tree of arthropods. Gene tree of another Aqp4-like subfamily (PripL) shows an unexpected basal split between acariform mites (Acariformes) and other arachnids. A closer inspection demonstrated that the PripL evolved quickly and has been under differential selection pressure in acariform mites. Evidence is provided that the evolutionarily ancient Glp subfamily (i.e. aquaglyceroporin) is significantly expanded in terrestrial arachnids compared with their marine relatives. Finally, in spite of the phylogenetic diversity, there exists conservation of some exons in size, functional domain, and intron-insertion phase: an 81-bp and a 218-bp exon, respectively, in apq4-like and glp genes across Eumetazoa lineages including arachnids and human beings. Both exons encode the carboxyl-terminal NPA motif, implying the coding and splicing pressure during hundreds of million years of animal evolution. Hypotheses were tested to explore the possible link between these findings and arachnid terrestrialisation.

Genomic and transcriptomic analyses support a silk gland origin of spider venom glands

BACKGROUND

Spiders comprise a hyperdiverse lineage of predators with venom systems, yet the origin of functionally novel spider venom glands remains unclear. Previous studies have hypothesized that spider venom glands originated from salivary glands or evolved from silk-producing glands present in early chelicerates. However, there is insufficient molecular evidence to indicate similarity among them. Here, we provide comparative analyses of genome and transcriptome data from various lineages of spiders and other arthropods to advance our understanding of spider venom gland evolution.

RESULTS

We generated a chromosome-level genome assembly of a model spider species, the common house spider (Parasteatoda tepidariorum). Module preservation, GO semantic similarity, and differentially upregulated gene similarity analyses demonstrated a lower similarity in gene expressions between the venom glands and salivary glands compared to the silk glands, which questions the validity of the salivary gland origin hypothesis but unexpectedly prefers to support the ancestral silk gland origin hypothesis. The conserved core network in the venom and silk glands was mainly correlated with transcription regulation, protein modification, transport, and signal transduction pathways. At the genetic level, we found that many genes in the venom gland-specific transcription modules show positive selection and upregulated expressions, suggesting that genetic variation plays an important role in the evolution of venom glands.

CONCLUSIONS

This research implies the unique origin and evolutionary path of spider venom glands and provides a basis for understanding the diverse molecular characteristics of venom systems.

Molecular Basis for the Selectivity of the Succinate Dehydrogenase Inhibitor Cyflumetofen between Pest and Predatory Mites

Acaricides that act as inhibitors of the mitochondrial succinate dehydrogenase (SDHIs) provide excellent control of phytophagous mites but display limited toxicity to predatory mites and other beneficial organisms. However, the molecular mechanism of selectivity is not fully understood. Here, we first confirm that SDHI acaricides are over 10,000-fold more toxic to spider mites than predatory mites. Next, we show that differential penetration, pro-acaricide activation, or metabolism are most likely not the main reason for this selectivity. In contrast, the inhibition of AB-1 on the SDH target is approximately 200-fold more potent in spider mites compared to predatory mites, revealing strong target-site selectivity. Strikingly, a key motif associated with differential binding was identified and validated by gene editing in Drosophila. Our findings contribute to understanding the selectivity of SDHIs, which can be used for the rational design of selective acaricides in support of an integrated pest management.

Do mites eat and run? A systematic review of feeding and dispersal strategies

Dispersal is an important process affecting the survival of organisms and the structure and dynamics of communities and ecosystems in space and time. It is a multiphase phenomenon influenced by many internal and external factors. Dispersal syndromes can be complicated, but they are vital to our knowledge of the biology of any organism. We analysed dispersal ability in mites (Acariformes and Parasitiformes), a highly diverse group of wingless arthropods, taking into consideration various modes of dispersal, feeding strategies, body size and the number of articles published for each species. Based on 174 articles summarized for this study, it appears that mites are opportunistic when it comes to dispersal, regardless of their feeding habits, and are often able to adopt several different strategies as needs arise. Moreover, we find a significant positive relationship between the amount of research effort that was put into studying a given species and the number of modes of dispersal that were described. The most salient conclusion to be drawn from this positive correlation is that additional studies are needed, especially on a broader set of mite taxa, until the aforementioned correlation is no longer demonstrably significant.

Using PhyloSuite for molecular phylogeny and tree-based analyses

Phylogenetic analysis has entered the genomics (multilocus) era. For less experienced researchers, conquering the large number of software programs required for a multilocus-based phylogenetic reconstruction can be somewhat daunting and time-consuming. PhyloSuite, a software with a user-friendly GUI, was designed to make this process more accessible by integrating multiple software programs needed for multilocus and single-gene phylogenies and further streamlining the whole process. In this protocol, we aim to explain how to conduct each step of the phylogenetic pipeline and tree-based analyses in PhyloSuite. We also present a new version of PhyloSuite (v1.2.3), wherein we fixed some bugs, made some optimizations, and introduced some new functions, including a number of tree-based analyses, such as signal-to-noise calculation, saturation analysis, spurious species identification, and etc. The step-by-step protocol includes background information (i.e., what the step does), reasons (i.e., why do the step), and operations (i.e., how to do it). This protocol will help researchers quick-start their way through the multilocus phylogenetic analysis, especially those interested in conducting organelle-based analyses.

Re-evaluating and dating myriapod diversification with phylotranscriptomics under a regime of dense taxon sampling

Recent transcriptomic studies of myriapod phylogeny have been based on relatively small datasets with <40 myriapod terminals and variably supported or contradicted the traditional morphological groupings of Progoneata and Dignatha. Here we amassed a large dataset of 104 myriapod terminals, including multiple species for each of the four myriapod classes. Across the tree, most nodes are stable and well supported. Most analyses across a range of gene occupancy levels provide moderate to strong support for a deep split of Myriapoda into Symphyla + Pauropoda (=Edafopoda) and an uncontradicted grouping of Chilopoda + Diplopoda (=Pectinopoda nov.), as in other recent transcriptome-based analyses; no analysis recovers Progoneata or Dignatha as clades. As in all recent multi-locus and phylogenomic studies, chilopod interrelationships resolve with Craterostigmus excluded from Amalpighiata rather than uniting with other centipedes with maternal brood care in Phylactometria. Diplopod ordinal interrelationships are largely congruent with morphology-based classifications. Chilognathan clades that are not invariably advocated by morphologists include Glomerida + Glomeridesmida, such that the volvation-related characters of pill millipedes may be convergent, and Stemmiulida + Polydesmida more closely allied to Juliformia than to Callipodida + Chordeumatida. The latter relationship implies homoplasy in spinnerets and contradicts Nematophora. A time-tree with nodes calibrated by 25 myriapod and six outgroup fossil terminals recovers Cambrian-Ordovician divergences for the deepest splits in Myriapoda, Edafopoda and Pectinopoda, predating the terrestrial fossil record of myriapods as in other published chronograms, whereas age estimates within Chilopoda and Diplopoda overlap with or do not appreciably predate the calibration fossils. The grouping of Chilopoda and Diplopoda is recovered in all our analyses and is formalized as Pectinopoda nov., named for the shared presence of mandibular comb lamellae. New taxonomic proposals for Chilopoda based on uncontradicted clades are Tykhepoda nov. for the three blind families of Scolopendromorpha that share a "sieve-type" gizzard, and Taktikospina nov. for Scolopendromorpha to the exclusion of Mimopidae.

Extraembryonic tissue in chelicerates: a review and outlook

The formation of extraembryonic membranes (EEMs) contributes to the proper development of many animals. In arthropods, the formation and function of EEMs have been studied best in insects. Regarding the development of extraembryonic tissue in chelicerates (spiders and relatives), most information is available for spiders (Araneae). Especially two populations of cells have been considered to represent EEMs in spiders. The first of these potential EEMs develops shortly after egg deposition, opposite to a radially symmetrical germ disc that forms in one hemisphere of the egg and encloses the yolk. The second tissue, which has been described as being extraembryonic is the so-called dorsal field, which is required to cover the dorsal part of the developing spider germ rudiment before proper dorsal closure. In this review, we summarize the current knowledge regarding the formation of potential extraembryonic structures in the Chelicerata. We describe the early embryogenesis of spiders and other chelicerates, with a special focus on the formation of the potential extraembryonic tissues.

This article is part of the theme issue ‘Extraembryonic tissues: exploring concepts, definitions and functions across the animal kingdom’.

One-way ticket to the blue: A large-scale, dated phylogeny revealed asymmetric land-to-water transitions in acariform mites (Acari: Acariformes)

Acariform mites are an ancient and megadiverse lineage that may have experienced a complex pattern of invasions into terrestrial and aquatic habitats. These among-realm transitions may relate to periods of turmoil in Earth's history or be simply results of uneven biodiversity patterns across habitats. Here, we inferred a dated, representative acariform phylogeny (five genes, 9,200 bp aligned, 367 terminals belonging to 150 ingroup plus 15 outgroup families, 23 fossil calibration points) which was used to infer transitions between marine/freshwater/terrestrial habitats. We detected four unambiguous transitions from terrestrial to freshwater habitats (Hydrozetes, Naiadacarus, Fusohericia, Afronothrus, Homocaligus); one from freshwater to marine (Pontarachnidae), and four from marine to brackish or freshwater transitions (all among Halacaridae: Acarothrix; Halacarellus petiti; Copidognathus sp.; clade Limnohalacarus + Soldanellonyx + Porohalacarus + Porolohmannella). One transition to the sea was inferred ambiguously with respect to the ancestor being either terrestrial or freshwater (Hyadesiidae), and another must be most carefully examined by adding potential related taxa (Selenoribatidae + Fortuyniidae). Finally, we inferred a single, remarkable transition from aquatic to terrestrial habitats involving early evolution of the large and ecologically diverse lineage: the ancestor of the Halacaridae + Parasitengona clade was probably freshwater given our dataset, thus making terrestrial Parasitengona secondarily terrestrial. Overall, our results suggested a strong asymmetry in environmental transitions: the majority occurred from terrestrial to aquatic habitats. This asymmetry is probably linked to mites' biological properties and uneven biodiversity patterns across habitats rather than Earth's geological history. Since the land holds more acariform diversity than water habitats, a shift from the former is more likely than from the latter. We inferred the following relationships: alicid endeostigmatid + eriophyoid (Alycidae, (Nanorchestidae, (Nematalycidae, Eriophyoidea))) being sister group to the remaining Acariformes: (proteonematalycid Endeostigmata, alicorhagiid Endeostigmata, Trombidiformes, Oribatida (including Astigmata)). Trombidiform relationships had several novel rearrangements: (i) traditional Eupodina lacked support for the inclusion of Bdelloidea; (ii) Teneriffidae, traditionally placed among Anystina, was consistently recovered in a clade including Heterostigmata in Eleutherengona; (iii) several lineages, such as Adamystidae, Paratydeidae, Caeculidae and Erythracaridae, were recovered in a large clade along other Anystina and Eleutherengona, suggesting single origins of several fundamental character states, such as the reduction of the cheliceral fixed digit and development of the palpal thumb-claw complex.

The dynamic changes of genes revealed that persistently overexpressed genes drive the evolution of cyflumetofen resistance in Tetranychus cinnabarinus

Changes in gene expression are associated with the evolution of pesticide resistance in arthropods. In this study, transcriptome sequencing was performed in 3 different resistance levels (low, L; medium, M; and high, H) of cyflumetofen-resistant strain (YN-CyR). A total of 1 685 genes, including 97 detoxification enzyme genes, were upregulated in all 3 stages, of which 192 genes, including 11 detoxification enzyme genes, showed a continuous increase in expression level with resistance development (L to H). RNA interference experiments showed that overexpression of 7 genes (CYP392A1, TcGSTd05, CCE06, CYP389A1, TcGSTz01, CCE59, and CYP389C2) is involved in the development of cyflumetofen resistance in Tetranychus cinnabarinus. The recombinant CYP392A1 can effectively metabolize cyflumetofen, while CCE06 can bind and sequester cyflumetofen in vitro. We compared 2 methods for rapid screening of resistance molecular markers, including short-term induction and 1-time high-dose selection. Two detoxification enzyme genes were upregulated in the field susceptible strain (YN-S) by induction with 20% lethal concentration (LC₂₀ ) of cyflumetofen. However, 16 detoxification enzyme genes were upregulated by 1-time selection with LC₈₀ of cyflumetofen. Interestingly, the 16 genes were overexpressed in all 3 resistance stages. These results indicated that 1 685 genes that were upregulated at the L stage constituted the basis of cyflumetofen resistance, of which 192 genes in which upregulation continued to increase were the main driving force for the development of resistance. Moreover, the 1-time high-dose selection is an efficient way to rapidly obtain the resistance-related genes that can aid in the development of resistance markers and resistance management in mites.

Was There a Cambrian Explosion on Land? The Case of Arthropod Terrestrialization

Arthropods, the most diverse form of macroscopic life in the history of the Earth, originated in the sea. Since the early Cambrian, at least ~518 million years ago, these animals have dominated the oceans of the world. By the Silurian-Devonian, the fossil record attests to arthropods becoming the first animals to colonize land, However, a growing body of molecular dating and palaeontological evidence suggests that the three major terrestrial arthropod groups (myriapods, hexapods, and arachnids), as well as vascular plants, may have invaded land as early as the Cambrian-Ordovician. These dates precede the oldest fossil evidence of those groups and suggest an unrecorded continental "Cambrian explosion" a hundred million years prior to the formation of early complex terrestrial ecosystems in the Silurian-Devonian. We review the palaeontological, phylogenomic, and molecular clock evidence pertaining to the proposed Cambrian terrestrialization of the arthropods. We argue that despite the challenges posed by incomplete preservation and the scarcity of early Palaeozoic terrestrial deposits, the discrepancy between molecular clock estimates and the fossil record is narrower than is often claimed. We discuss strategies for closing the gap between molecular clock estimates and fossil data in the evolution of early ecosystems on land.

The organizing role of Wnt signaling pathway during arthropod posterior growth

Wnt signaling pathways are recognized for having major roles in tissue patterning and cell proliferation. In the last years, remarkable progress has been made in elucidating the molecular and cellular mechanisms that underlie sequential segmentation and axial elongation in various arthropods, and the canonical Wnt pathway has emerged as an essential factor in these processes. Here we review, with a comparative perspective, the current evidence concerning the participation of this pathway during posterior growth, its degree of conservation among the different subphyla within Arthropoda and its relationship with the rest of the gene regulatory network involved. Furthermore, we discuss how this signaling pathway could regulate segmentation to establish this repetitive pattern and, at the same time, probably modulate different cellular processes precisely coupled to axial elongation. Based on the information collected, we suggest that this pathway plays an organizing role in the formation of the body segments through the regulation of the dynamic expression of segmentation genes, via controlling the caudal gene, at the posterior region of the embryo/larva, that is necessary for the correct sequential formation of body segments in most arthropods and possibly in their common segmented ancestor. On the other hand, there is insufficient evidence to link this pathway to axial elongation by controlling its main cellular processes, such as convergent extension and cell proliferation. However, conclusions are premature until more studies incorporating diverse arthropods are carried out.

DNA Barcoding and Phylogeny of Acari Species Based on ITS and COI Markers

Acari harbor numerous minute species of agricultural economic importance, mainly Tetranychidae and Phytoseiidae. Great efforts have been established by means of recovering morphological, molecular, and phylogenetic traits for species identification. Traditional identification still relies on external diagnostic characters, which are limited and usually exhibit large phenotypic plasticity within the species, rendering them useless for species delimitation and identification. We decided to increase the number of sequences of the Acari mitochondrial COI (Cytochrome C oxidase I) marker and ITS nuclear ribosomal DNA region for species identification in Tetranychidae and Phytoseiidae. The molecular data allow us to establish species boundaries and phylogenetic relationships among several clades of Acari, mainly Tetranychidae and Phytoseiidae. Sequence comparisons between complete COI and the Acari mitochondrial COI, ITS1-5,8S-ITS2, and ITS2 among all Acari sequences have demonstrated that the selected regions, even small, gave enough informative positions for both species’ identification and phylogenetic studies. Analyses of both DNA regions have unveiled their use as species identification characters, with special emphasis on Acari mitochondrial COI for Tetranychidae and Phytoseiidae species in comparison with the Folmer fragment, which has been universally used as a barcode marker. We demonstrated that the Acari mitochondrial COI region is also a suitable marker to establish a barcode dataset for Acari identification. Our phylogenetic analyses are congruent with other recent works, showing that Acari is a monophyletic group, of which Astigmata, Ixodida, Mesostigmata, Oribatida, and Prostigmata are also monophyletic.

Comparative Genomics Reveals Insights into the Divergent Evolution of Astigmatic Mites and Household Pest Adaptations

Highly diversified astigmatic mites comprise many medically important human household pests such as house dust mites causing ∼1–2% of all allergic diseases globally; however, their evolutionary origin and diverse lifestyles including reversible parasitism have not been illustrated at the genomic level, which hampers allergy prevention and our exploration of these household pests. Using six high-quality assembled and annotated genomes, this study not only refuted the monophyly of mites and ticks, but also thoroughly explored the divergence of Acariformes and the diversification of astigmatic mites. In monophyletic Acariformes, Prostigmata known as notorious plant pests first evolved, and then rapidly evolving Astigmata diverged from soil oribatid mites. Within astigmatic mites, a wide range of gene families rapidly expanded via tandem gene duplications, including ionotropic glutamate receptors, triacylglycerol lipases, serine proteases and UDP glucuronosyltransferases. Gene diversification after tandem duplications provides many genetic resources for adaptation to sensing environmental signals, digestion, and detoxification in rapidly changing household environments. Many gene decay events only occurred in the skin-burrowing parasitic mite Sarcoptes scabiei. Throughout the evolution of Acariformes, massive horizontal gene transfer events occurred in gene families such as UDP glucuronosyltransferases and several important fungal cell wall lytic enzymes, which enable detoxification and digestive functions and provide perfect drug targets for pest control. This comparative study sheds light on the divergent evolution and quick adaptation to human household environments of astigmatic mites and provides insights into the genetic adaptations and even control of human household pests.

A microCT-based atlas of the central nervous system and midgut in sea spiders (Pycnogonida) sheds first light on evolutionary trends at the family level

Background

Pycnogonida (sea spiders) is the sister group of all other extant chelicerates (spiders, scorpions and relatives) and thus represents an important taxon to inform early chelicerate evolution. Notably, phylogenetic analyses have challenged traditional hypotheses on the relationships of the major pycnogonid lineages (families), indicating external morphological traits previously used to deduce inter-familial affinities to be highly homoplastic. This erodes some of the support for phylogenetic information content in external morphology and calls for the study of additional data classes to test and underpin in-group relationships advocated in molecular analyses. In this regard, pycnogonid internal anatomy remains largely unexplored and taxon coverage in the studies available is limited.

Results

Based on micro-computed X-ray tomography and 3D reconstruction, we created a comprehensive atlas of in-situ representations of the central nervous system and midgut layout in all pycnogonid families. Beyond that, immunolabeling for tubulin and synapsin was used to reveal selected details of ganglionic architecture. The ventral nerve cord consistently features an array of separate ganglia, but some lineages exhibit extended composite ganglia, due to neuromere fusion. Further, inter-ganglionic distances and ganglion positions relative to segment borders vary, with an anterior shift in several families. Intersegmental nerves target longitudinal muscles and are lacking if the latter are reduced. Across families, the midgut displays linear leg diverticula. In Pycnogonidae, however, complex multi-branching diverticula occur, which may be evolutionarily correlated with a reduction of the heart.

Conclusions

Several gross neuroanatomical features are linked to external morphology, including intersegmental nerve reduction in concert with trunk segment fusion, or antero-posterior ganglion shifts in partial correlation to trunk elongation/compaction. Mapping on a recent phylogenomic phylogeny shows disjunct distributions of these traits. Other characters show no such dependency and help to underpin closer affinities in sub-branches of the pycnogonid tree, as exemplified by the tripartite subesophageal ganglion of Pycnogonidae and Rhynchothoracidae. Building on this gross anatomical atlas, future studies should now aim to leverage the full potential of neuroanatomy for phylogenetic interrogation by deciphering pycnogonid nervous system architecture in more detail, given that pioneering work on neuron subsets revealed complex character sets with unequivocal homologies across some families.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12983-022-00459-8.

Highly diversified mitochondrial genomes provide new evidence for interordinal relationships in the Arachnida

Arachnida is an exceptionally diverse class in the Arthropoda, consisting of 20 orders and playing crucial roles in the terrestrial ecosystems. However, their interordinal relationships have been debated for over a century. Rearranged or highly rearranged mitochondrial genomes (mitogenomes) were consistently found in this class, but their various extent in different lineages and efficiency for resolving arachnid phylogenies are unclear. Here, we reconstructed phylogenetic trees using mitogenome sequences of 290 arachnid species to decipher interordinal relationships as well as diversification through time. Our results recovered monophyly of ten orders (i.e. Amblypygi, Araneae, Ixodida, Mesostigmata, Opiliones, Pseudoscorpiones, Ricinulei, Sarcoptiformes, Scorpiones and Solifugae), while rejecting monophyly of the Trombidiformes due to the unstable position of the Eriophyoidea. The monophyly of Acari (subclass) was rejected, possibly due to the long-branch attraction of the Pseudoscorpiones. The monophyly of Arachnida was further rejected because the Xiphosura nested within arachnid orders with unstable positions. Mitogenomes that are highly rearranged in mites but less rearranged or conserved in the remaining lineages point to their exceptional diversification in mite orders; however, shared derived mitochondrial (mt) gene clusters were found within superfamilies rather than interorders, confusing phylogenetic signals in arachnid interordinal relationships. Molecular dating results show that arachnid orders have ancient origins, ranging from the Ordovician to the Carboniferous, yet have significantly diversified since the Cretaceous in orders Araneae, Mesostigmata, Sarcoptiformes, and Trombidiformes. By summarizing previously resolved key positions of some orders, we propose a plausible arachnid tree of life. Our results underline a more precise framework for interordinal phylogeny in the Arachnida and provide new insights into their ancient evolution.

Eggs to long-legs: embryonic staging of the harvestman Phalangium opilio (Opiliones), an emerging model arachnid

Background

The comparative embryology of Chelicerata has greatly advanced in recent years with the integration of classical studies and genetics, prominently spearheaded by developmental genetic works in spiders. Nonetheless, the understanding of the evolution of development and polarization of embryological characters in Chelicerata is presently limited, as few non-spider species have been well studied. A promising focal species for chelicerate evo-devo is the daddy-long-legs (harvestman) Phalangium opilio, a member of the order Opiliones. Phalangium opilio, breeds prolifically and is easily accessible in many parts of the world, as well as tractable in a laboratory setting. Resources for this species include developmental transcriptomes, a draft genome, and protocols for RNA interference, but a modern staging system is critically missing for this emerging model system.

Results

We present a staging system of P. opilio embryogenesis that spans the most important morphogenetic events with respect to segment formation, appendage elongation and head development. Using time-lapse imaging, confocal microscopy, colorimetric in situ hybridization, and immunohistochemistry, we tracked the development of synchronous clutches from egg laying to adulthood. We describe key events in segmentation, myogenesis, neurogenesis, and germ cell formation.

Conclusion

Considering the phylogenetic position of Opiliones and the unduplicated condition of its genome (in contrast to groups like spiders and scorpions), this species is poised to serve as a linchpin for comparative studies in arthropod development and genome evolution. The staging system presented herein provides a valuable reference for P. opilio that we anticipate being useful to the arthropod evo-devo community, with the goal of revitalizing research in the comparative development of non-spider arachnids.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12983-022-00454-z.

Selectivity and molecular stress responses to classical and botanical acaricides in the predatory mite Phytoseiulus persimilis Athias-Henriot (Acari: Phytoseiidae)

BACKGROUND

Acaricide application remains an integral component of integrated pest management (IPM) for the two-spotted spider mite Tetranychus urticae. Species and strains of phytoseiid predatory mites vary significantly in their response to acaricides. For the success of IPM, it is imperative to identify the determinants of selectivity and molecular stress responses of acaricides in predatory mites.

RESULTS

The three classical acaricides bifenazate, cyflumetofen, and fenbutatin oxide did not affect the survival and fecundity of Phytoseiulus persimilis regardless of the route of exposure. Selectivity of the orange oil and terpenoid blend-based botanical acaricides was low via a combination of direct exposure, acaricide-laced diet, and residual exposure but improved when limiting exposure only to diet. To gain insights into the molecular stress responses, the transcriptome of P. persimilis was assembled. Subsequent gene expression analysis of predatory mites orally exposed to fenbutatin oxide and orange oil yielded only a limited xenobiotic stress response. In contrast, P. persimilis exhibited target-site resistance mutations, including I260M in SdhB, I1017M in CHS1, and kdr and super-kdr in VGSC. Extending the screen using available Phytoseiidae sequences uncovered I136T, S141F in cytb, G119S in AChE, and A2083V in ACC, well-known target-sites of acaricides.

CONCLUSION

Selectivity of the tested botanical acaricides to P. persimilis was low but could be enhanced by restricting exposure to a single route. Differential gene expression analysis did not show a robust induced stress response after sublethal exposure. In contrast, this study uncovered target-site mutations that may help to explain the physiological selectivity of several classical acaricides to phytoseiid predators.

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.

Circular RNA, circ1-3p, is Involved in Cyflumetofen Resistance by Acting as a Competitive RNA against miR-1-3p in Tetranychus cinnabarinus

As a newly recognized type of noncoding RNA, circular RNA can mediate a variety of physiological changes in mammals by regulating the post-transcriptional expression level of genes. However, the function of circRNA in the evolution of pesticide resistance in arthropods is still unknown. In this study, 2546 circRNAs were identified in Tetranychus cinnabarinus by transcriptome sequencing. The differentially expressed gene analysis indicated that 44 circRNAs were overexpressed in a cyflumetofen-resistant strain, of which a circRNA (named circ1-3p) was found to contain the response elements of miR-1-3p, an miRNA that is involved in cyflumetofen resistance by targeting TcGSTm04. The circular structure of circ1-3p was further determined using a divergent primer. The results of different molecular assays in vitro and in vivo showed that circ1-3p can compete with TcGSTm04 in miR-1-3p binding. The colocalization of circ1-3p and miR-1-3p was found using fluorescence in situ hybridization, suggesting that circ1-3p can directly sponge miR-1-3p in T. cinnabarinus. In addition, silencing the expression of circ1-3p resulted in the upregulation of miR-1-3p and the downregulation of TcGSTm04 as well as a significant increase in the sensitivity of T. cinnabarinus to cyflumetofen. All these pieces of evidence indicates that overexpressed circ1-3p promotes the expression of TcGSTm04 through sponging miR-1-3p, thereby involving in the evolution of cyflumetofen resistance in T. cinnabarinus.

The visual pathway in sea spiders (Pycnogonida) displays a simple serial layout with similarities to the median eye pathway in horseshoe crabs

BACKGROUND

Phylogenomic studies over the past two decades have consolidated the major branches of the arthropod tree of life. However, especially within the Chelicerata (spiders, scorpions, and kin), interrelationships of the constituent taxa remain controversial. While sea spiders (Pycnogonida) are firmly established as sister group of all other extant representatives (Euchelicerata), euchelicerate phylogeny itself is still contested. One key issue concerns the marine horseshoe crabs (Xiphosura), which recent studies recover either as sister group of terrestrial Arachnida or nested within the latter, with significant impact on postulated terrestrialization scenarios and long-standing paradigms of ancestral chelicerate traits. In potential support of a nested placement, previous neuroanatomical studies highlighted similarities in the visual pathway of xiphosurans and some arachnopulmonates (scorpions, whip scorpions, whip spiders). However, contradictory descriptions of the pycnogonid visual system hamper outgroup comparison and thus character polarization.

RESULTS

To advance the understanding of the pycnogonid brain and its sense organs with the aim of elucidating chelicerate visual system evolution, a wide range of families were studied using a combination of micro-computed X-ray tomography, histology, dye tracing, and immunolabeling of tubulin, the neuropil marker synapsin, and several neuroactive substances (including histamine, serotonin, tyrosine hydroxylase, and orcokinin). Contrary to previous descriptions, the visual system displays a serial layout with only one first-order visual neuropil connected to a bilayered arcuate body by catecholaminergic interneurons. Fluorescent dye tracing reveals a previously reported second visual neuropil as the target of axons from the lateral sense organ instead of the eyes.

CONCLUSIONS

Ground pattern reconstruction reveals remarkable neuroanatomical stasis in the pycnogonid visual system since the Ordovician or even earlier. Its conserved layout exhibits similarities to the median eye pathway in euchelicerates, especially in xiphosurans, with which pycnogonids share two median eye pairs that differentiate consecutively during development and target one visual neuropil upstream of the arcuate body. Given multiple losses of median and/or lateral eyes in chelicerates, and the tightly linked reduction of visual processing centers, interconnections between median and lateral visual neuropils in xiphosurans and arachnopulmonates are critically discussed, representing a plausible ancestral condition of taxa that have retained both eye types.

Comparative anatomy of the rostrosoma of Solifugae, Pseudoscorpiones and Acari

We compare the microscopic anatomy of the mouthparts of representative species of Solifugae, Pseudoscorpiones and Parasitiformes (Acari). Specifically, we focus on the epistome, the labrum, the lateral lips (= endites of the pedipalpal coxae) and the musculature of the pharyngeal suction pump. We provide evidence that the labrum is reduced in Solifugae, but present and functional in Pseudoscorpiones and Acari. The epistome constitutes the entire dorsal face of the rostrosoma in Solifugae, but is internalized into the prosoma in Pseudoscorpiones. In Acari, the epistome shows an ancestral morphology, probably close to the ground pattern of chelicerates. The lateral lips of Solifugae contribute to the ventral face of the rostrosoma and the two lips of the mouth opening. In Solifugae, the ventral rostrosoma also includes a sclerite that might derive from a tritosternum. In Pseudoscorpiones, the lateral lips remain independent of the rostrosoma, they interlock ventral to the rostrosoma forming a perioral space. Here, the rostrosoma has an unpaired ventral lip of unresolved morphological origin, which is, however, clearly distinct from the lateral lips of Solifugae. The pharyngeal suction pump differs in all three clades in attachment, number of muscles and origin of muscles. We interpret the data as evidence for independent, parallel evolution of elements of the ground pattern of the (eu)chelicerate mouth parts. Based on the morphological elements of a common euchelicerate ground plan, the rostrosoma evolved independently in the three clades. We reject earlier hypotheses that consider the rostrosoma a character to support a phylogenetic relationship of the three clades.

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 remains 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 lineages 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.

High Value-Added Application of Natural Products in Crop Protection: Semisynthesis and Acaricidal Activity of Limonoid-Type Derivatives and Investigation of Their Biocompatible O/W Nanoemulsions as Agronanopesticide Candidates

The increasingly serious resistance of Tetranychus cinnabarinus Boisduval to a wide range of insecticides/acaricides poses a major challenge to their control. The citrus processing industry generates a huge quantity of various wastes that contain many limonoids. To effectively utilize these byproducts and discover more potent green acaricidal molecules as sustainable alternatives for traditional resistant pesticides, various limonoid-type derivatives (halogenated/seven-membered lactam derivatives of obacunone and halogenated/oxime esters/oxime ethers/seven-membered lactam derivatives of limonin) were synthesized based on a diversity-oriented synthetic strategy. The key steric configurations of 10 derivatives were further confirmed by X-ray crystallography. Compound 9m, which displayed greater than 9.7-fold potent acaricidal activity of limonin, was of preeminence. In addition, some interesting structure-activity relationships were observed. Moreover, a biocompatible O/W nanoemulsion delivery system was used to prepare the limonin-based agronanoacaricide, which exhibited pronounced control efficiency against T. cinnabarinus Boisduval in the greenhouse. This systematic investigation will provide valuable information and guidance for future value-added applications of novel eco-friendly natural product-based nanopesticides.

What Is an “Arachnid”? Consensus, Consilience, and Confirmation Bias in the Phylogenetics of Chelicerata

The basal phylogeny of Chelicerata is one of the opaquest parts of the animal Tree of Life, defying resolution despite application of thousands of loci and millions of sites. At the forefront of the debate over chelicerate relationships is the monophyly of Arachnida, which has been refuted by most analyses of molecular sequence data. A number of phylogenomic datasets have suggested that Xiphosura (horseshoe crabs) are derived arachnids, refuting the traditional understanding of arachnid monophyly. This result is regarded as controversial, not least by paleontologists and morphologists, due to the widespread perception that arachnid monophyly is unambiguously supported by morphological data. Moreover, some molecular datasets have been able to recover arachnid monophyly, galvanizing the belief that any result that challenges arachnid monophyly is artefactual. Here, we explore the problems of distinguishing phylogenetic signal from noise through a series of in silico experiments, focusing on datasets that have recently supported arachnid monophyly. We assess the claim that filtering by saturation rate is a valid criterion for recovering Arachnida. We demonstrate that neither saturation rate, nor the ability to assemble a molecular phylogenetic dataset supporting a given outcome with maximal nodal support, is a guarantor of phylogenetic accuracy. Separately, we review empirical morphological phylogenetic datasets to examine characters supporting Arachnida and the downstream implication of a single colonization of terrestrial habitats. We show that morphological support of arachnid monophyly is contingent upon a small number of ambiguous or incorrectly coded characters, most of these tautologically linked to adaptation to terrestrial habitats.

Small body size of pseudoscorpions and a distinct architecture of the ovary: A step to miniaturization?

Chelicerata, the second largest subphylum of Arthropoda, includes invertebrates with a wide range of body size. Pseudoscorpions are among small or miniature chelicerates which exhibit several morphological, anatomical, and developmental features related to miniaturization, e.g., replacement of book lungs by tracheae, unpaired gonads, and matrotrophic development of the embryos outside the female body, in the brood sac. In this paper, we show the ovary structure of two pseudoscorpion species, Cheiridium museorum and Apocheiridium ferum (Cheiridiidae). Both cheiridiids are one of the smallest pseudoscorpions. The results of our observations conducted in light, transmission electron, and confocal microscopy demonstrate that the ovary of C. museorum and A. ferum, displays a significant structural difference that is unusual for chelicerates. The difference concerns the spatially restricted position of the germarium. We show that such ovary architecture results in a significantly reduced number of growing oocytes and in consequence a reduced number of deposited eggs. A centrally located germarium implies also a modified pattern of ovary development during oocyte growth due to long distance migration of the germline and the accompanying somatic cells. Herein, we postulate that such an ovary structure is related to the pseudoscorpion's small body size and it is a step towards miniaturization in the smaller pseudoscorpions species.

The evolution of insect biodiversity

Insects comprise over half of all described animal species. Together with the Protura (coneheads), Collembola (springtails) and Diplura (two-pronged bristletails), insects form the Hexapoda, a terrestrial arthropod lineage characterised by possessing six legs. Exponential growth of genome-scale data for the hexapods has substantially altered our understanding of the origin and evolution of insect biodiversity. Phylogenomics has provided a new framework for reconstructing insect evolutionary history, resolving their position among the arthropods and some long-standing internal controversies such as the placement of the termites, twisted-winged insects, lice and fleas. However, despite the greatly increased size of phylogenomic datasets, contentious relationships among key insect clades remain unresolved. Further advances in insect phylogeny cannot rely on increased depth and breadth of genome and taxon sequencing. Improved modelling of the substitution process is fundamental to countering tree-reconstruction artefacts, while gene content, modelling of duplications and deletions, and comparative morphology all provide complementary lines of evidence to test hypotheses emerging from the analysis of sequence data. Finally, the integration of molecular and morphological data is key to the incorporation of fossil species within insect phylogeny. The emerging integrated framework of insect evolution will help explain the origins of insect megadiversity in terms of the evolution of their body plan, species diversity and ecology. Future studies of insect phylogeny should build upon an experimental, hypothesis-driven approach where the robustness of hypotheses generated is tested against increasingly realistic evolutionary models as well as complementary sources of phylogenetic evidence.

Phylogenomic Subsampling and the Search for Phylogenetically Reliable Loci

Phylogenomic subsampling is a procedure by which small sets of loci are selected from large genome-scale data sets and used for phylogenetic inference. This step is often motivated by either computational limitations associated with the use of complex inference methods or as a means of testing the robustness of phylogenetic results by discarding loci that are deemed potentially misleading. Although many alternative methods of phylogenomic subsampling have been proposed, little effort has gone into comparing their behavior across different data sets. Here, I calculate multiple gene properties for a range of phylogenomic data sets spanning animal, fungal, and plant clades, uncovering a remarkable predictability in their patterns of covariance. I also show how these patterns provide a means for ordering loci by both their rate of evolution and their relative phylogenetic usefulness. This method of retrieving phylogenetically useful loci is found to be among the top performing when compared with alternative subsampling protocols. Relatively common approaches such as minimizing potential sources of systematic bias or increasing the clock-likeness of the data are found to fare worse than selecting loci at random. Likewise, the general utility of rate-based subsampling is found to be limited: loci evolving at both low and high rates are among the least effective, and even those evolving at optimal rates can still widely differ in usefulness. This study shows that many common subsampling approaches introduce unintended effects in off-target gene properties and proposes an alternative multivariate method that simultaneously optimizes phylogenetic signal while controlling for known sources of bias.

The genome of a daddy-long-legs (Opiliones) illuminates the evolution of arachnid appendages

Chelicerate arthropods exhibit dynamic genome evolution, with ancient whole-genome duplication (WGD) events affecting several orders. Yet, genomes remain unavailable for a number of poorly studied orders, such as Opiliones (daddy-long-legs), which has hindered comparative study. We assembled the first harvestman draft genome for the species Phalangium opilio, which bears elongate, prehensile appendages, made possible by numerous distal articles called tarsomeres. Here, we show that the genome of P. opilio exhibits a single Hox cluster and no evidence of WGD. To investigate the developmental genetic basis for the quintessential trait of this group-the elongate legs-we interrogated the function of the Hox genes Deformed (Dfd) and Sex combs reduced (Scr), and a homologue of Epidermal growth factor receptor (Egfr). Knockdown of Dfd incurred homeotic transformation of two pairs of legs into pedipalps, with dramatic shortening of leg segments in the longest leg pair, whereas homeosis in L3 is only achieved upon double Dfd + Scr knockdown. Knockdown of Egfr incurred shortened appendages and the loss of tarsomeres. The similarity of Egfr loss-of-function phenotypic spectra in insects and this arachnid suggest that repeated cooption of EGFR signalling underlies the independent gains of supernumerary tarsomeres across the arthropod tree of life.

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.

Fine structure of the epicuticular secretion coat and associated glands of Pedipalpi and Palpigradi (Arachnida)

Pedipalpi Latreille, 1810 is a poorly studied clade of arachnids comprising the whip spiders (Amblypygi Thorell, 1883), short-tailed whip scorpions (Schizomida Petrunkevitch, 1945) and whip scorpions (Thelyphonida Cambridge, 1872). It has recently been shown that whip spiders coat their exoskeleton with a solid cement layer (cerotegument) that forms elaborate microstructures and turns the cuticle into a super-hydrophobic state. The amblypygid cerotegument provides taxonomic information due to its fine structural diversity, but its presence and variation in the sister groups was previously unknown. The present contribution reports the surface structure of the cuticle in species of Palpigradi, Thelyphonida, and Schizomida to determine if these taxa possess a solid epicuticular secretion coat. Scanning electron microscopy revealed that in addition to Amblypygi only species of Thelyphonida possess solid epicuticular secretion layers. Unlike in Amblypygi, in the Thelyphonida this layer does not usually form microstructures and is less rigidly attached to the underlying cuticle. A species of Typopeltis Pocock, 1894, which exhibited globular structures analogous to the amblypygid cerotegument, was an exception. Glandular structures associated with cement secretions in Amblypygi and Thelyphonida were considered homologous due to similar structure. Solid epicuticular secretion coats were absent from Schizomida, which is interpreted as a secondary loss despite the presence of slit-like glandular openings that appear to produce such epicuticular secretions. The micro-whip scorpion order Palpigradi Thorell, 1900 exhibited markedly different cuticular surface structures and lacked solid epicuticular secretions, consistent with the hypothesis that this order is not closely related to Pedipalpi. These results enhance the knowledge of the small, enigmatic orders of Arachnida.

Six-state Amino Acid Recoding is not an Effective Strategy to Offset Compositional Heterogeneity and Saturation in Phylogenetic Analyses

Six-state amino acid recoding strategies are commonly applied to combat the effects of compositional heterogeneity and substitution saturation in phylogenetic analyses. While these methods have been endorsed from a theoretical perspective, their performance has never been extensively tested. Here, we test the effectiveness of six-state recoding approaches by comparing the performance of analyses on recoded and non-recoded data sets that have been simulated under gradients of compositional heterogeneity or saturation. In our simulation analyses, non-recoding approaches consistently outperform six-state recoding approaches. Our results suggest that six-state recoding strategies are not effective in the face of high saturation. Furthermore, while recoding strategies do buffer the effects of compositional heterogeneity, the loss of information that accompanies six-state recoding outweighs its benefits. In addition, we evaluate recoding schemes with 9, 12, 15, and 18 states and show that these consistently outperform six-state recoding. Our analyses of other recoding schemes suggest that under conditions of very high compositional heterogeneity, it may be advantageous to apply recoding using more than six states, but we caution that applying any recoding should include sufficient justification. Our results have important implications for the more than 90 published papers that have incorporated six-state recoding, many of which have significant bearing on relationships across the tree of life. [Compositional heterogeneity; Dayhoff 6-state recoding; S&R 6-state recoding; six-state amino acid recoding; substitution saturation.]

Evidence for sponges as sister to all other animals from partitioned phylogenomics with mixture models and recoding

Resolving the relationships between the major lineages in the animal tree of life is necessary to understand the origin and evolution of key animal traits. Sponges, characterized by their simple body plan, were traditionally considered the sister group of all other animal lineages, implying a gradual increase in animal complexity from unicellularity to complex multicellularity. However, the availability of genomic data has sparked tremendous controversy as some phylogenomic studies support comb jellies taking this position, requiring secondary loss or independent origins of complex traits. Here we show that incorporating site-heterogeneous mixture models and recoding into partitioned phylogenomics alleviates systematic errors that hamper commonly-applied phylogenetic models. Testing on real datasets, we show a great improvement in model-fit that attenuates branching artefacts induced by systematic error. We reanalyse key datasets and show that partitioned phylogenomics does not support comb jellies as sister to other animals at either the supermatrix or partition-specific level.

Interrogating Genomic-Scale Data to Resolve Recalcitrant Nodes in the Spider Tree of Life

Genome-scale data sets are converging on robust, stable phylogenetic hypotheses for many lineages; however, some nodes have shown disagreement across classes of data. We use spiders (Araneae) as a system to identify the causes of incongruence in phylogenetic signal between three classes of data: exons (as in phylotranscriptomics), noncoding regions (included in ultraconserved elements [UCE] analyses), and a combination of both (as in UCE analyses). Gene orthologs, coded as amino acids and nucleotides (with and without third codon positions), were generated by querying published transcriptomes for UCEs, recovering 1,931 UCE loci (codingUCEs). We expected that congeners represented in the codingUCE and UCEs data would form clades in the presence of phylogenetic signal. Noncoding regions derived from UCE sequences were recovered to test the stability of relationships. Phylogenetic relationships resulting from all analyses were largely congruent. All nucleotide data sets from transcriptomes, UCEs, or a combination of both recovered similar topologies in contrast with results from transcriptomes analyzed as amino acids. Most relationships inferred from low-occupancy data sets, containing several hundreds of loci, were congruent across Araneae, as opposed to high occupancy data matrices with fewer loci, which showed more variation. Furthermore, we found that low-occupancy data sets analyzed as nucleotides (as is typical of UCE data sets) can result in more congruent relationships than high occupancy data sets analyzed as amino acids (as in phylotranscriptomics). Thus, omitting data, through amino acid translation or via retention of only high occupancy loci, may have a deleterious effect in phylogenetic reconstruction.

Transcriptome analysis revealed that multiple genes were related to the cyflumetofen resistance of Tetranychus cinnabarinus (Boisduval)

Metabolic resistance is one of the main causes of acaricide resistance. Many previous studies focused on the function of specific genes in insecticides/acaricides resistance. However, during the development of resistance, the overall dynamic of expression levels of detoxification enzyme genes in mites is still unclear. Tetranychus cinnabarinus, a major agricultural pest, which is notorious for developing resistance to acaricides rapidly. In this study, a field susceptible strain (YS) was continuously selected for 16, 25 and 32 generations, and developed to low resistance (7.83-fold, L), medium resistance (17.23-fold, M) and high resistance (86.05-fold, H), respectively. Transcriptome sequencing was performed in YS, L, M and H strains. Overall, compared with YS strain, the number of differential expression genes increased slightly with the development of cyflumetofen-resistance. As for detoxification genes, the median of fold change of up-regulated P450、CCE and GST genes was higher than those of all up-regulated genes in three resistance level, but only the number and the median of fold change of up-regulated P450 genes was increased slightly with the development of resistance. In addition, synergism experiments also proved that P450 and GST genes were the major contributors to the metabolic resistance of cyflumetofen of T. cinnabarinus. These results showed that the resistance of T. cinnabarinus to cyflumetofen was related to many resistant genes, among which P450 genes could play crucial roles in cyflumefen resistance.

Horseshoe crab genomes reveal the evolution of genes and microRNAs after three rounds of whole genome duplication

Whole genome duplication (WGD) has occurred in relatively few sexually reproducing invertebrates. Consequently, the WGD that occurred in the common ancestor of horseshoe crabs ~135 million years ago provides a rare opportunity to decipher the evolutionary consequences of a duplicated invertebrate genome. Here, we present a high-quality genome assembly for the mangrove horseshoe crab Carcinoscorpius rotundicauda (1.7 Gb, N50 = 90.2 Mb, with 89.8% sequences anchored to 16 pseudomolecules, 2n = 32), and a resequenced genome of the tri-spine horseshoe crab Tachypleus tridentatus (1.7 Gb, N50 = 109.7 Mb). Analyses of gene families, microRNAs, and synteny show that horseshoe crabs have undergone three rounds (3R) of WGD. Comparison of C. rotundicauda and T. tridentatus genomes from populations from several geographic locations further elucidates the diverse fates of both coding and noncoding genes. Together, the present study represents a cornerstone for improving our understanding of invertebrate WGD events on the evolutionary fates of genes and microRNAs, at both the individual and population level. We also provide improved genomic resources for horseshoe crabs, of applied value for breeding programs and conservation of this fascinating and unusual invertebrate lineage.

Genomic Determination of Reproductive Mode in Facultatively Parthenogenetic Opiliones

Sexual reproduction may pose myriad short-term costs to females. Despite these costs, sexual reproduction is near ubiquitous. Facultative parthenogenesis is theorized to mitigate some of the costs of sex, as individuals can participate in occasional sex to limit costs while obtaining many benefits. However, most theoretical models assume sexual reproduction is fixed following mating, with no possibility of clutches of mixed reproductive ontogeny. Therefore, we asked: if coercive males are present at high frequency in a population of facultative parthenogens, will their clutches be solely sexually produced, or will there be evidence of sexually and asexually-produced offspring? How will their offspring production compare to conspecifics in low-frequency male populations? We addressed our questions by collecting females and egg clutches of the facultatively parthenogenetic Opiliones species Leiobunum manubriatum and L. globosum. In L. manubriatum, females from populations with few males were not significantly more fecund than females from populations with higher male relative frequency, despite the potential release of the former from sexual conflict. We used 3 genotyping methods along with a custom set of DNA capture probes to reveal that offspring of L. manubriatum from these high male populations were primarily produced via asexual reproduction. This is surprising because sex ratios in these southern populations approach equality, increasing the probability for females to encounter mates and produce offspring sexually. We additionally found evidence for reproductive polymorphisms within populations. Rapid and accurate SNP genotyping data will continue to allow us to address broader evolutionary questions regarding the role of facultative reproductive modes in the maintenance of sex.

Phylogenomic Analysis of Ultraconserved Elements Resolves the Evolutionary and Biogeographic History of Segmented Trapdoor Spiders

The segmented trapdoor spiders (Liphistiidae) are the sole surviving family of the suborder Mesothelae, which forms the sister lineage to all other living spiders. Liphistiids have retained a number of plesiomorphic traits and their present-day distribution is limited to East and Southeast Asia. Studying this group has the potential to shed light on the deep evolutionary history of spiders, but the phylogeny and divergence times of the family have not been resolved with confidence. We performed phylogenomic and molecular dating analyses of 2,765 ultraconserved element loci from 185 liphistiid taxa. Our analyses show that the crown group of Liphistiidae appeared in the mid-Cretaceous at 102 Ma (95% credibility interval 92-113 Ma), but it was not until the Neogene that much of the diversification within the family occurred in mainland Southeast and East Asia. This diversification was coincident with tectonic events such as the extension of the East Asian continental margin, as well as geological upheavals in Indochina induced by the collision between India and Asia. Our study highlights the important role of major tectonic events in shaping the evolutionary history, present-day diversity, and geographical distribution of mesothele and liphistiid spiders.

Systemic paralogy and function of retinal determination network homologs in arachnids

BACKGROUND

Arachnids are important components of cave ecosystems and display many examples of troglomorphisms, such as blindness, depigmentation, and elongate appendages. Little is known about how the eyes of arachnids are specified genetically, let alone the mechanisms for eye reduction and loss in troglomorphic arachnids. Additionally, duplication of Retinal Determination Gene Network (RDGN) homologs in spiders has convoluted functional inferences extrapolated from single-copy homologs in pancrustacean models.

RESULTS

We investigated a sister species pair of Israeli cave whip spiders, Charinus ioanniticus and C. israelensis (Arachnopulmonata, Amblypygi), of which one species has reduced eyes. We generated embryonic transcriptomes for both Amblypygi species, and discovered that several RDGN homologs exhibit duplications. We show that duplication of RDGN homologs is systemic across arachnopulmonates (arachnid orders that bear book lungs), rather than being a spider-specific phenomenon. A differential gene expression (DGE) analysis comparing the expression of RDGN genes in field-collected embryos of both species identified candidate RDGN genes involved in the formation and reduction of eyes in whip spiders. To ground bioinformatic inference of expression patterns with functional experiments, we interrogated the function of three candidate RDGN genes identified from DGE using RNAi in the spider Parasteatoda tepidariorum. We provide functional evidence that one of these paralogs, sine oculis/Six1 A (soA), is necessary for the development of all arachnid eye types.

CONCLUSIONS

Our work establishes a foundation to investigate the genetics of troglomorphic adaptations in cave arachnids, and links differential gene expression to an arthropod eye phenotype for the first time outside of Pancrustacea. Our results support the conservation of at least one RDGN component across Arthropoda and provide a framework for identifying the role of gene duplications in generating arachnid eye diversity.

A genome-wide screening for RNAi pathway proteins in Acari

BACKGROUND

RNA interference (RNAi) is a highly conserved, sequence-specific gene silencing mechanism present in Eukaryotes. Three RNAi pathways are known, namely micro-RNA (miRNA), piwi-interacting RNA (piRNA) and short interfering RNA (siRNA). However, little knowledge exists about the proteins involved in these pathways in Acari. Moreover, variable successes has been obtained in gene knockdown via siRNA pathway in their functional genomics and management. We hypothesized that the clue may be in the variability of the composition and the efficacy of siRNA machinery among Acari.

RESULTS

Both comparative genomic analyses and domain annotation suggest that all the analyzed species have homologs of putative core proteins that mediate cleaving of targeted genes via the three RNAi pathways. We identified putative homologs of Caenorhabditis elegans RNA-dependent RNA polymerase (RdRP) protein in all species though no secondary Argonaute homologs that operate with this protein in siRNA amplification mechanism were found, suggesting that the siRNA amplification mechanism present in Acari may be distinct from that described in C. elegans. Moreover, the genomes of these species do not encode homologs of C. elegans systemic RNAi defective-1 (Sid-1) protein that mediate silencing of the mRNA target throughout the treated organisms suggesting that the phenomena of systemic RNAi that has been reported in some Acari species probably occur through a different mechanism. However, homologs of putative RNAi spreading defective-3 (Rsd-3) protein and scavenger receptors namely Eater and SR-CI that mediate endocytosis cellular update of dsRNA in C. elegans and Drosophila melanogaster were found in Acari genomes. This result suggests that cellular dsRNA uptake in Acari is endocytosis-dependent. Detailed phylogenetic analyses of core RNAi pathway proteins in the studied species revealed that their evolution is compatible with the proposed monophyletic evolution of this group.

CONCLUSIONS

Our analyses have revealed the potential activity of all three pathways in Acari. Still, much experimental work remains to be done to confirm the mechanisms behind these pathways in particular those that govern systemic/parental RNAi and siRNA amplification in Acari. Disclosure of these mechanisms will facilitate the development of new and specific management tools for the harmful species and enrichment of the beneficial species.