The house spider genome reveals an ancient whole-genome duplication during arachnid evolution

Evidence of positive selection on six spider developmental genes

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

Human follicular mites: Ectoparasites becoming symbionts

Most humans carry mites in the hair follicles of their skin for their entire lives. Follicular mites are the only metazoans tha continuously live on humans. We propose that Demodex folliculorum (Acari) represents a transitional stage from a host-injuring obligate parasite to an obligate symbiont. Here, we describe the profound impact of this transition on the genome and physiology of the mite. Genome sequencing revealed that the permanent host association of D. folliculorum led to an extensive genome reduction through relaxed selection and genetic drift, resulting in the smallest number of protein-coding genes yet identified among panarthropods. Confocal microscopy revealed that this gene loss coincided with an extreme reduction in the number of cells. Single uninucleate muscle cells are sufficient to operate each of the three segments that form each walking leg. While it has been assumed that the reduction of the cell number in parasites starts early in development, we identified a greater total number of cells in the last developmental stage (nymph) than in the terminal adult stage, suggesting that reduction starts at the adult or ultimate stage of development. This is the first evolutionary step in an arthropod species adopting a reductive, parasitic or endosymbiotic lifestyle. Somatic nuclei show underreplication at the diploid stage. Novel eye structures or photoreceptors as well as a unique human host melatonin-guided day/night rhythm are proposed for the first time. The loss of DNA repair genes coupled with extreme endogamy might have set this mite species on an evolutionary dead-end trajectory.

Chromosome-level genome assembly of the black widow spider Latrodectus elegans illuminates composition and evolution of venom and silk proteins

BACKGROUND

The black widow spider has both extraordinarily neurotoxic venom and three-dimensional cobwebs composed of diverse types of silk. However, a high-quality reference genome for the black widow spider was still unavailable, which hindered deep understanding and application of the valuable biomass.

FINDINGS

We assembled the Latrodectus elegans genome, including a genome size of 1.57 Gb with contig N50 of 4.34 Mb and scaffold N50 of 114.31 Mb. Hi-C scaffolding assigned 98.08% of the genome to 14 pseudo-chromosomes, and with BUSCO, completeness analysis revealed that 98.4% of the core eukaryotic genes were completely present in this genome. Annotation of this genome identified that repetitive sequences account for 506.09 Mb (32.30%) and 20,167 protein-coding genes, and specifically, we identified 55 toxin genes and 26 spidroins and provide preliminary analysis of their composition and evolution.

CONCLUSIONS

We present the first chromosome-level genome assembly of a black widow spider and provide substantial toxin and spidroin gene resources. These high-qualified genomic data add valuable resources from a representative spider group and contribute to deep exploration of spider genome evolution, especially in terms of the important issues on the diversification of venom and web-weaving pattern. The sequence data are also firsthand templates for further application of the spider biomass.

Genome of a giant isopod, Bathynomus jamesi, provides insights into body size evolution and adaptation to deep-sea environment

BACKGROUND

The deep-sea may be regarded as a hostile living environment, due to low temperature, high hydrostatic pressure, and limited food and light. Isopods, a species-rich group of crustaceans, are widely distributed across different environments including the deep sea and as such are a useful model for studying adaptation, migration, and speciation. Similar to other deep-sea organisms, giant isopods have larger body size than their shallow water relatives and have large stomachs and fat bodies presumably to store organic reserves. In order to shed light on the genetic basis of these large crustaceans adapting to the oligotrophic environment of deep-sea, the high-quality genome of a deep-sea giant isopod Bathynomus jamesi was sequenced and assembled.

RESULTS

B. jamesi has a large genome of 5.89 Gb, representing the largest sequenced crustacean genome to date. Its large genome size is mainly attributable to the remarkable proliferation of transposable elements (84%), which may enable high genome plasticity for adaptive evolution. Unlike its relatives with small body size, B. jamesi has expanded gene families related to pathways of thyroid and insulin hormone signaling that potentially contribute to its large body size. Transcriptomic analysis showed that some expanded gene families related to glycolysis and vesicular transport were specifically expressed in its digestive organs. In addition, comparative genomics and gene expression analyses in six tissues suggested that B. jamesi has inefficient lipid degradation, low basal metabolic rate, and bulk food storage, suggesting giant isopods adopt a more efficient mechanism of nutrient absorption, storage, and utilization to provide sustained energy supply for their large body size.

CONCLUSIONS

Taken together, the giant isopod genome may provide a valuable resource for understanding body size evolution and adaptation mechanisms of macrobenthic organisms to deep-sea environments.

Embryonic expression patterns of Wnt genes in the RTA-clade spider Cupiennius salei

Spiders represent widely used model organisms for chelicerate and even arthropod development and evolution. Wnt genes are important and evolutionary conserved factors that control and regulate numerous developmental processes. Recent studies comprehensively investigated the complement and expression of spider Wnt genes revealing conserved as well as diverged aspects of their expression and thus (likely) function among different groups of spiders representing Mygalomorphae (tarantulas), and both main groups of Araneae (true spiders) (Haplogynae/Synspermiata and Entelegynae). The allegedly most modern/derived group of entelegyne spiders is represented by the RTA-clade of which no comprehensive data on Wnt expression were available prior to this study. Here, we investigated the embryonic expression of all Wnt genes of the RTA-clade spider Cupiennius salei. We found that most of the Wnt expression patterns are conserved between Cupiennius and other spiders, especially more basally branching species. Surprisingly, most differences in Wnt gene expression are seen in the common model spider Parasteatoda tepidariorum (a non-RTA clade entelegyne species). These results show that data and conclusions drawn from research on one member of a group of animals (or any other organism) cannot necessarily be extrapolated to the group as a whole, and instead highlight the need for comprehensive taxon sampling.

Chromosome-level genome of the globe skimmer dragonfly (Pantala flavescens)

BACKGROUND

The globe skimmer dragonfly (Pantala flavescens) is a notable Odonata insect distributed in nature fields and farmlands worldwide, and it is commonly recognized as a natural enemy because it preys on agricultural pests and health pests. As one of the sister groups of winged insects, odonatan species are key to understanding the evolution of insect wings.

FINDINGS

We present a high-quality reference genome of P. flavescens, which is the first chromosome-level genome in the Palaeoptera (Odonata and Ephemeroptera). The assembled genome size was 662 Mb, with a contig N50 of 16.2 Mb. Via Hi-C scaffolding, 648 Mb (97.9%) of contig sequences were clustered, ordered, and assembled into 12 large scaffolds, each corresponding to a natural chromosome. The X chromosome was identified by sequence coverage depth. The repetitive sequences and gene density of the X chromosome are similar to those of autosomal sequences, but the X chromosome shows a much lower degree of heterozygosity. Our analysis shows that the effective population size experienced 3 declining events, which may have been caused by climate change and environmental pollution.

CONCLUSIONS

The genome of P. flavescens provides more information on the biology and evolution of insects and will help for the use of this species in pest control.

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.

Phylogenetic analysis of forkhead transcription factors in the Panarthropoda

Fox genes encode transcription factors that contain a DNA binding domain, the forkhead domain, and are known from diverse animal species. The exact homology of the Fox genes of different species is debated and this makes inferences about the evolution of the Fox genes, and their duplications and losses difficult. We have performed phylogenetic analyses of the Fox gene complements of 32 panarthropod species. Our results confirm an ancestral complement of FoxA, FoxB, FoxC, FoxD, FoxF, FoxG, FoxJ1, FoxJ2/3, FoxK, FoxL1, FoxL2, FoxN1/4, FoxN2/3, FoxO, FoxP, and FoxQ2 in the Arthropoda, and additionally FoxH and FoxQ1 in the Panarthropoda (including tardigrades and onychophorans). We identify a novel Fox gene sub-family, that we designate as FoxT that includes two genes in Drosophila melanogaster, Circadianly Regulated Gene (Crg-1) and forkhead domain 3F (fd3F). In a very recent paper, the same new Fox gene sub-family was identified in insects (Lin et al. 2021). Our analysis confirms the presence of FoxT and shows that its members are present throughout Panarthropoda. We show that the hitherto unclassified gene CG32006 from the fly Drosophila melanogaster belongs to FoxJ1. We also detect gene losses: FoxE and FoxM were lost already in the panarthropod ancestor, whereas the loss of FoxH occurred in the arthropod ancestor. Finally, we find an ortholog of FoxQ1 in the bark scorpion Centruroides sculpturatus, confirmed not only by phylogenetic analysis, but also by forming an evolutionarily conserved gene cluster with FoxF, FoxC, and FoxL1. This suggests that FoxQ1 belongs to the ancestral Fox gene complement in panarthropods and also in chelicerates, but has been lost at the base of the mandibulate arthropods.

Genome size evolution in the diverse insect order Trichoptera

Background

Genome size is implicated in the form, function, and ecological success of a species. Two principally different mechanisms are proposed as major drivers of eukaryotic genome evolution and diversity: polyploidy (i.e., whole-genome duplication) or smaller duplication events and bursts in the activity of repetitive elements. Here, we generated de novo genome assemblies of 17 caddisflies covering all major lineages of Trichoptera. Using these and previously sequenced genomes, we use caddisflies as a model for understanding genome size evolution in diverse insect lineages.

Results

We detect a ∼14-fold variation in genome size across the order Trichoptera. We find strong evidence that repetitive element expansions, particularly those of transposable elements (TEs), are important drivers of large caddisfly genome sizes. Using an innovative method to examine TEs associated with universal single-copy orthologs (i.e., BUSCO genes), we find that TE expansions have a major impact on protein-coding gene regions, with TE-gene associations showing a linear relationship with increasing genome size. Intriguingly, we find that expanded genomes preferentially evolved in caddisfly clades with a higher ecological diversity (i.e., various feeding modes, diversification in variable, less stable environments).

Conclusion

Our findings provide a platform to test hypotheses about the potential evolutionary roles of TE activity and TE-gene associations, particularly in groups with high species, ecological, and functional diversities.

It takes Two: Discovery of Spider Pax2 Duplicates Indicates Prominent Role in Chelicerate Central Nervous System, Eye, as Well as External Sense Organ Precursor Formation and Diversification After Neo- and Subfunctionalization

Paired box genes are conserved across animals and encode transcription factors playing key roles in development, especially neurogenesis. Pax6 is a chief example for functional conservation required for eye development in most bilaterian lineages except chelicerates. Pax6 is ancestrally linked and was shown to have interchangeable functions with Pax2. Drosophila melanogaster Pax2 plays an important role in the development of sensory hairs across the whole body. In addition, it is required for the differentiation of compound eyes, making it a prime candidate to study the genetic basis of arthropod sense organ development and diversification, as well as the role of Pax genes in eye development. Interestingly, in previous studies identification of chelicerate Pax2 was either neglected or failed. Here we report the expression of two Pax2 orthologs in the common house spider Parasteatoda tepidariorum, a model organism for chelicerate development. The two Pax2 orthologs most likely arose as a consequence of a whole genome duplication in the last common ancestor of spiders and scorpions. Pax2.1 is expressed in the peripheral nervous system, including developing lateral eyes and external sensilla, as well as the ventral neuroectoderm of P. tepidariorum embryos. This not only hints at a conserved dual role of Pax2/5/8 orthologs in arthropod sense organ development but suggests that in chelicerates, Pax2 could have acquired the role usually played by Pax6. For the other paralog, Pt-Pax2.2, expression was detected in the brain, but not in the lateral eyes and the expression pattern associated with sensory hairs differs in timing, pattern, and strength. To achieve a broader phylogenetic sampling, we also studied the expression of both Pax2 genes in the haplogyne cellar spider Pholcus phalangioides. We found that the expression difference between paralogs is even more extreme in this species, since Pp-Pax2.2 shows an interesting expression pattern in the ventral neuroectoderm while the expression in the prosomal appendages is strictly mesodermal. This expression divergence indicates both sub- and neofunctionalization after Pax2 duplication in spiders and thus presents an opportunity to study the evolution of functional divergence after gene duplication and its impact on sense organ diversification.

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.

A masculinizing supergene underlies an exaggerated male reproductive morph in a spider

In many species, individuals can develop into strikingly different morphs, which are determined by a simple Mendelian locus. How selection shapes loci that control complex phenotypic differences remains poorly understood. In the spider Oedothorax gibbosus, males develop either into a 'hunched' morph with conspicuous head structures or as a fast-developing 'flat' morph with a female-like appearance. We show that the hunched-determining allele contains a unique genomic fragment of approximately 3 megabases that is absent in the flat-determining allele. This fragment comprises dozens of genes that duplicated from genes found at the same as well as different chromosomes. All functional duplicates, including a duplicate of the key sexual differentiation regulatory gene doublesex, show male-specific expression, which illustrates their integrated role as a masculinizing supergene. Our findings demonstrate how extensive indel polymorphisms and duplications of regulatory genes may contribute to the evolution of co-adapted gene clusters, sex-limited reproductive morphs and the enigmatic evolution of exaggerated sexual traits in general.

Lack of evidence for conserved parasegmental grooves in arthropods

In the arthropod model species Drosophila melanogaster, a dipteran fly, segmentation of the anterior–posterior body axis is under control of a hierarchic gene cascade. Segmental boundaries that form morphological grooves are established posteriorly within the segmental expression domain of the segment-polarity gene (SPG) engrailed (en). More important for the development of the fly, however, are the parasegmental boundaries that are established at the interface of en expressing cells and anteriorly adjacent wingless (wg) expressing cells. In Drosophila, both segmental and transient parasegmental grooves form. The latter are positioned anterior to the expression of en. Although the function of the SPGs in establishing and maintaining segmental and parasegmental boundaries is highly conserved among arthropods, parasegmental grooves have only been reported for Drosophila, and a spider (Cupiennius salei). Here, we present new data on en expression, and re-evaluate published data, from four distantly related spiders, including Cupiennius, and a distantly related chelicerate, the harvestman Phalangium opilio. Gene expression analysis of en genes in these animals does not corroborate the presence of parasegmental grooves. Consequently, our data question the general presence of parasegmental grooves in arthropods.

MicroRNAs as Indicators into the Causes and Consequences of Whole-Genome Duplication Events

Whole-genome duplications (WGDs) have long been considered the causal mechanism underlying dramatic increases to morphological complexity due to the neo-functionalization of paralogs generated during these events. Nonetheless, an alternative hypothesis suggests that behind the retention of most paralogs is not neo-functionalization, but instead the degree of the inter-connectivity of the intended gene product, as well as the mode of the WGD itself. Here, we explore both the causes and consequences of WGD by examining the distribution, expression, and molecular evolution of microRNAs (miRNAs) in both gnathostome vertebrates as well as chelicerate arthropods. We find that although the number of miRNA paralogs tracks the number of WGDs experienced within the lineage, few of these paralogs experienced changes to the seed sequence, and thus are functionally equivalent relative to their mRNA targets. Nonetheless, in gnathostomes, although the retention of paralogs following the 1R autotetraploidization event is similar across the two subgenomes, the paralogs generated by the gnathostome 2R allotetraploidization event are retained in higher numbers on one subgenome relative to the second, with the miRNAs found on the preferred subgenome showing both higher expression of mature miRNA transcripts and slower molecular evolution of the precursor miRNA sequences. Importantly, WGDs do not result in the creation of miRNA novelty, nor do WGDs correlate to increases in complexity. Instead, it is the number of miRNA seed sequences in the genome itself that not only better correlate to instances in complexification, but also mechanistically explain why complexity increases when new miRNA families are established.

Widespread retention of ohnologs in key developmental gene families following whole-genome duplication in arachnopulmonates

Whole-genome duplications (WGDs) have occurred multiple times during animal evolution, including in lineages leading to vertebrates, teleosts, horseshoe crabs, and arachnopulmonates. These dramatic events initially produce a wealth of new genetic material, generally followed by extensive gene loss. It appears, however, that developmental genes such as homeobox genes, signaling pathway components and microRNAs are frequently retained as duplicates (so-called ohnologs) following WGD. These not only provide the best evidence for WGD, but an opportunity to study its evolutionary consequences. Although these genes are well studied in the context of vertebrate WGD, similar comparisons across the extant arachnopulmonate orders are patchy. We sequenced embryonic transcriptomes from two spider species and two amblypygid species and surveyed three important gene families, Hox, Wnt, and frizzled, across these and 12 existing transcriptomic and genomic resources for chelicerates. We report extensive retention of putative ohnologs, further supporting the ancestral arachnopulmonate WGD. We also found evidence of consistent evolutionary trajectories in Hox and Wnt gene repertoires across three of the six arachnopulmonate orders, with interorder variation in the retention of specific paralogs. We identified variation between major clades in spiders and are better able to reconstruct the chronology of gene duplications and losses in spiders, amblypygids, and scorpions. These insights shed light on the evolution of the developmental toolkit in arachnopulmonates, highlight the importance of the comparative approach within lineages, and provide substantial new transcriptomic data for future study.

The Tetragnatha kauaiensis Genome Sheds Light on the Origins of Genomic Novelty in Spiders

Spiders (Araneae) have a diverse spectrum of morphologies, behaviors, and physiologies. Attempts to understand the genomic-basis of this diversity are often hindered by their large, heterozygous, and AT-rich genomes with high repeat content resulting in highly fragmented, poor-quality assemblies. As a result, the key attributes of spider genomes, including gene family evolution, repeat content, and gene function, remain poorly understood. Here, we used Illumina and Dovetail Chicago technologies to sequence the genome of the long-jawed spider Tetragnatha kauaiensis, producing an assembly distributed along 3,925 scaffolds with an N50 of ∼2 Mb. Using comparative genomics tools, we explore genome evolution across available spider assemblies. Our findings suggest that the previously reported and vast genome size variation in spiders is linked to the different representation and number of transposable elements. Using statistical tools to uncover gene-family level evolution, we find expansions associated with the sensory perception of taste, immunity, and metabolism. In addition, we report strikingly different histories of chemosensory, venom, and silk gene families, with the first two evolving much earlier, affected by the ancestral whole genome duplication in Arachnopulmonata (∼450 Ma) and exhibiting higher numbers. Together, our findings reveal that spider genomes are highly variable and that genomic novelty may have been driven by the burst of an ancient whole genome duplication, followed by gene family and transposable element expansion.

Insights into the Karyotype Evolution of Charinidae, the Early-Diverging Clade of Whip Spiders (Arachnida: Amblypygi)

Whip spiders (Amblypygi) represent an ancient order of tetrapulmonate arachnids with a low diversity. Their cytogenetic data are confined to only a few reports. Here, we analyzed the family Charinidae, a lineage almost at the base of the amblypygids, providing an insight into the ancestral traits and basic trajectories of amblypygid karyotype evolution. We performed Giemsa staining, selected banding techniques, and detected 18S ribosomal DNA and telomeric repeats by fluorescence in situ hybridization in four Charinus and five Sarax species. Both genera exhibit a wide range of diploid chromosome numbers (2n = 42-76 and 22-74 for Charinus and Sarax, respectively). The 2n reduction was accompanied by an increase of proportion of biarmed elements. We further revealed a single NOR site (probably an ancestral condition for charinids), the presence of a (TTAGG)n telomeric motif localized mostly at the chromosome ends, and an absence of heteromorphic sex chromosomes. Our data collectively suggest a high pace of karyotype repatterning in amblypygids, with probably a high ancestral 2n and its subsequent gradual reduction by fusions, and the action of pericentric inversions, similarly to what has been proposed for neoamblypygids. The possible contribution of fissions to charinid karyotype repatterning, however, cannot be fully ruled out.

A chelicerate Wnt gene expression atlas: novel insights into the complexity of arthropod Wnt-patterning

The Wnt genes represent a large family of secreted glycoprotein ligands that date back to early animal evolution. Multiple duplication events generated a set of 13 Wnt families of which 12 are preserved in protostomes. Embryonic Wnt expression patterns (Wnt-patterning) are complex, representing the plentitude of functions these genes play during development. Here, we comprehensively investigated the embryonic expression patterns of Wnt genes from three species of spiders covering both main groups of true spiders, Haplogynae and Entelegynae, a mygalomorph species (tarantula), as well as a distantly related chelicerate outgroup species, the harvestman Phalangium opilio. All spiders possess the same ten classes of Wnt genes, but retained partially different sets of duplicated Wnt genes after whole genome duplication, some of which representing impressive examples of sub- and neo-functionalization. The harvestman, however, possesses a more complete set of 11 Wnt genes but with no duplicates. Our comprehensive data-analysis suggests a high degree of complexity and evolutionary flexibility of Wnt-patterning likely providing a firm network of mutational protection. We discuss the new data on Wnt gene expression in terms of their potential function in segmentation, posterior elongation, and appendage development and critically review previous research on these topics. We conclude that earlier research may have suffered from the absence of comprehensive gene expression data leading to partial misconceptions about the roles of Wnt genes in development and evolution.

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.

The biology and evolution of spider venoms

Spiders are diverse, predatory arthropods that have inhabited Earth for around 400 million years. They are well known for their complex venom systems that are used to overpower their prey. Spider venoms contain many proteins and peptides with highly specific and potent activities suitable for biomedical or agrochemical applications, but the key role of venoms as an evolutionary innovation is often overlooked, even though this has enabled spiders to emerge as one of the most successful animal lineages. In this review, we discuss these neglected biological aspects of spider venoms. We focus on the morphology of spider venom systems, their major components, biochemical and chemical plasticity, as well as ecological and evolutionary trends. We argue that the effectiveness of spider venoms is due to their unprecedented complexity, with diverse components working synergistically to increase the overall potency. The analysis of spider venoms is difficult to standardize because they are dynamic systems, fine-tuned and modified by factors such as sex, life-history stage and biological role. Finally, we summarize the mechanisms that drive spider venom evolution and highlight the need for genome-based studies to reconstruct the evolutionary history and physiological networks of spider venom compounds with more certainty.

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.

The chromosome-scale assembly of the Canary Islands endemic spider Dysdera silvatica (Arachnida, Araneae) sheds light on the origin and genome structure of chemoreceptor gene families in chelicerates

Here, we present the chromosome-level genome assembly of Dysdera silvatica Schmidt, 1981, a nocturnal ground-dwelling spider endemic from the Canary Islands. The genus Dysdera has undergone a remarkable diversification in this archipelago mostly associated with shifts in the level of trophic specialization, becoming an excellent model to study the genomic drivers of adaptive radiations. The new assembly (1.37 Gb; scaffold N50 of 174.2 Mb), was performed using the chromosome conformation capture scaffolding technique, represents a continuity improvement of more than 4500 times with respect to the previous version. The seven largest scaffolds or pseudochromosomes, which cover 87% of the total assembly size, probably correspond with the seven chromosomes of the karyotype of this species, including a characteristic large X chromosome. To illustrate the value of this new resource we performed a comprehensive analysis of the two major arthropod chemoreceptor gene families (i.e., gustatory and ionotropic receptors). We identified 545 chemoreceptor sequences distributed across all pseudochromosomes, with a notable underrepresentation in the X chromosome. At least 54% of them localize in 83 genomic clusters with a significantly lower evolutionary distances between them than the average of the family, suggesting a recent origin of many of them. This chromosome-level assembly is the first high-quality genome representative of the Synspermiata clade, and just the third among spiders, representing a new valuable resource to gain insights into the structure and organization of chelicerate genomes, including the role that structural variants, repetitive elements and large gene families played in the extraordinary biology of spiders.

Early, nonlethal ploidy and genome size quantification using confocal microscopy in zebrafish embryos

Ploidy transitions through whole genome duplication have shaped evolution by allowing the sub- and neo-functionalization of redundant copies of highly conserved genes to express novel traits. The nuclear:cytoplasmic (n:c) ratio is maintained in polyploid vertebrates resulting in larger cells, but body size is maintained by a concomitant reduction in cell number. Ploidy can be manipulated easily in most teleosts, and the zebrafish, already well established as a model system for biomedical research, is therefore an excellent system in which to study the effects of increased cell size and reduced cell numbers in polyploids on development and physiology. Here we describe a novel technique using confocal microscopy to measure genome size and determine ploidy non-lethally at 48 h post-fertilization (hpf) in transgenic zebrafish expressing fluorescent histones. Volumetric analysis of myofiber nuclei using open-source software can reliably distinguish diploids and triploids from a mixed-ploidy pool of embryos for subsequent experimentation. We present an example of this by comparing heart rate between confirmed diploid and triploid embryos at 54 hpf.

Genome-Wide Identification of Soybean ABC Transporters Relate to Aluminum Toxicity

ATP-binding cassette (ABC) transporter proteins are a gene super-family in plants and play vital roles in growth, development, and response to abiotic and biotic stresses. The ABC transporters have been identified in crop plants such as rice and buckwheat, but little is known about them in soybean. Soybean is an important oil crop and is one of the five major crops in the world. In this study, 255 ABC genes that putatively encode ABC transporters were identified from soybean through bioinformatics and then categorized into eight subfamilies, including 7 ABCAs, 52 ABCBs, 48 ABCCs, 5 ABCDs, 1 ABCEs, 10 ABCFs, 111 ABCGs, and 21 ABCIs. Their phylogenetic relationships, gene structure, and gene expression profiles were characterized. Segmental duplication was the main reason for the expansion of the GmABC genes. Ka/Ks analysis suggested that intense purifying selection was accompanied by the evolution of GmABC genes. The genome-wide collinearity of soybean with other species showed that GmABCs were relatively conserved and that collinear ABCs between species may have originated from the same ancestor. Gene expression analysis of GmABCs revealed the distinct expression pattern in different tissues and diverse developmental stages. The candidate genes GmABCB23, GmABCB25, GmABCB48, GmABCB52, GmABCI1, GmABCI5, and GmABCI13 were responsive to Al toxicity. This work on the GmABC gene family provides useful information for future studies on ABC transporters in soybean and potential targets for the cultivation of new germplasm resources of aluminum-tolerant soybean.

Implications of a cheliceral axial duplication in Tetragnatha versicolor (Araneae: Tetragnathidae) for arachnid deuterocerebral appendage development

The homology of the arachnid chelicera with respect to other head appendages in Panarthropoda has long been debated. Gene expression data and the re-interpretation of early transitional fossils have supported the homology of the deutocerebrum and its associated appendages, implying a homology between primary antennae (mandibulates), chelicerae (euchelicerates), and chelifores (sea spiders). Nevertheless, comparatively little is known about the mechanistic basis of proximo-distal (PD) axis induction in chelicerates, much less the basis for cheliceral fate specification. Here, we describe a new cheliceral teratology in the spider Tetragnatha versicolor Walckenaer, 1841, which consists on a duplication of the PD axis of the left chelicera associated with a terminal secondary schistomely on the fang of the lower axis. This duplication offers clues as to potential shared mechanisms of PD axis formation in the chelicera. We review the state of knowledge on PD axis induction mechanisms in arthropods and identify elements of gene regulatory networks that are key for future functional experiments of appendage development in non-insect model systems. Such investigations would allow a better understanding of PD axis induction of modified and poorly studied arthropod limbs (e.g., chelicerae, chelifores, and ovigers).

Panarthropod tiptop/teashirt and spalt orthologs and their potential role as "trunk"-selector genes

BACKGROUND

In the vinegar fly Drosophila melanogaster, the homeodomain containing transcription factor Teashirt (Tsh) appears to specify trunk identity in concert with the function of the Hox genes. While in Drosophila there is a second gene closely related to tsh, called tiptop (tio), in other arthropods species only one copy exists (called tio/tsh). The expression of tsh and tio/tsh, respectively, is surprisingly similar among arthropods suggesting that its function as trunk selector gene may be conserved. Other research, for example on the beetle Tribolium castaneum, questions even conservation of Tsh function among insects. The zinc-finger transcription factor Spalt (Sal) is involved in the regulation of Drosophila tsh, but this regulatory interaction does not appear to be conserved in Tribolium either. Whether the function and interaction of tsh and sal as potential trunk-specifiers, however, is conserved is still unclear because comparative studies on sal expression (except for Tribolium) are lacking, and functional data are (if at all existing) restricted to Insecta.

RESULTS

Here, we provide additional data on arthropod tsh expression, show the first data on onychophoran tio/tsh expression, and provide a comprehensive investigation on sal expression patterns in arthropods and an onychophoran.

CONCLUSIONS

Our data support the idea that tio/tsh genes are involved in the development of "trunk" segments by regulating limb development. Our data suggest further that the function of Sal is indeed unlikely to be conserved in trunk vs head development like in Drosophila, but early expression of sal is in line with a potential homeotic function, at least in Arthropoda.

Lineage-level divergence of copepod glycerol transporters and the emergence of isoform-specific trafficking regulation

Transmembrane conductance of small uncharged solutes such as glycerol typically occurs through aquaglyceroporins (Glps), which are commonly encoded by multiple genes in metazoan organisms. To date, however, little is known concerning the evolution of Glps in Crustacea or what forces might underly such apparent gene redundancy. Here, we show that Glp evolution in Crustacea is highly divergent, ranging from single copy genes in species of pedunculate barnacles, tadpole shrimps, isopods, amphipods and decapods to up to 10 copies in diplostracan water fleas although with monophyletic origins in each lineage. By contrast the evolution of Glps in Copepoda appears to be polyphyletic, with surprisingly high rates of gene duplication occurring in a genera- and species-specific manner. Based upon functional experiments on the Glps from a parasitic copepod (Lepeophtheirus salmonis), we show that such lineage-level gene duplication and splice variation is coupled with a high rate of neofunctionalization. In the case of L. salmonis, splice variation of a given gene resulted in tissue- or sex-specific expression of the channels, with each variant evolving unique sites for protein kinase C (PKC)- or protein kinase A (PKA)-regulation of intracellular membrane trafficking. The combined data sets thus reveal that mutations favouring a high fidelity control of intracellular trafficking regulation can be a selection force for the evolution and retention of multiple Glps in copepods.

The genome of the venomous snail Lautoconus ventricosus sheds light on the origin of conotoxin diversity

BACKGROUND

Venoms are deadly weapons to subdue prey or deter predators that have evolved independently in many animal lineages. The genomes of venomous animals are essential to understand the evolutionary mechanisms involved in the origin and diversification of venoms.

RESULTS

Here, we report the chromosome-level genome of the venomous Mediterranean cone snail, Lautoconus ventricosus (Caenogastropoda: Conidae). The total size of the assembly is 3.59 Gb; it has high contiguity (N50 = 93.53 Mb) and 86.6 Mb of the genome assembled into the 35 largest scaffolds or pseudochromosomes. On the basis of venom gland transcriptomes, we annotated 262 complete genes encoding conotoxin precursors, hormones, and other venom-related proteins. These genes were scattered in the different pseudochromosomes and located within repetitive regions. The genes encoding conotoxin precursors were normally structured into 3 exons, which did not necessarily coincide with the 3 structural domains of the corresponding proteins. Additionally, we found evidence in the L. ventricosus genome for a past whole-genome duplication event by means of conserved gene synteny with the Pomacea canaliculata genome, the only one available at the chromosome level within Caenogastropoda. The whole-genome duplication event was further confirmed by the presence of a duplicated hox gene cluster. Key genes for gastropod biology including those encoding proteins related to development, shell formation, and sex were located in the genome.

CONCLUSIONS

The new high-quality L. ventricosus genome should become a reference for assembling and analyzing new gastropod genomes and will contribute to future evolutionary genomic studies among venomous animals.

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.

Evolutionary pattern of karyotypes and meiosis in pholcid spiders (Araneae: Pholcidae): implications for reconstructing chromosome evolution of araneomorph spiders

BACKGROUND

Despite progress in genomic analysis of spiders, their chromosome evolution is not satisfactorily understood. Most information on spider chromosomes concerns the most diversified clade, entelegyne araneomorphs. Other clades are far less studied. Our study focused on haplogyne araneomorphs, which are remarkable for their unusual sex chromosome systems and for the co-evolution of sex chromosomes and nucleolus organizer regions (NORs); some haplogynes exhibit holokinetic chromosomes. To trace the karyotype evolution of haplogynes on the family level, we analysed the number and morphology of chromosomes, sex chromosomes, NORs, and meiosis in pholcids, which are among the most diverse haplogyne families. The evolution of spider NORs is largely unknown.

RESULTS

Our study is based on an extensive set of species representing all major pholcid clades. Pholcids exhibit a low 2n and predominance of biarmed chromosomes, which are typical haplogyne features. Sex chromosomes and NOR patterns of pholcids are diversified. We revealed six sex chromosome systems in pholcids (X0, XY, X1X20, X1X2X30, X1X2Y, and X1X2X3X4Y). The number of NOR loci ranges from one to nine. In some clades, NORs are also found on sex chromosomes.

CONCLUSIONS

The evolution of cytogenetic characters was largely derived from character mapping on a recently published molecular phylogeny of the family. Based on an extensive set of species and mapping of their characters, numerous conclusions regarding the karyotype evolution of pholcids and spiders can be drawn. Our results suggest frequent autosome-autosome and autosome-sex chromosome rearrangements during pholcid evolution. Such events have previously been attributed to the reproductive isolation of species. The peculiar X1X2Y system is probably ancestral for haplogynes. Chromosomes of the X1X2Y system differ considerably in their pattern of evolution. In some pholcid clades, the X1X2Y system has transformed into the X1X20 or XY systems, and subsequently into the X0 system. The X1X2X30 system of Smeringopus pallidus probably arose from the X1X20 system by an X chromosome fission. The X1X2X3X4Y system of Kambiwa probably evolved from the X1X2Y system by integration of a chromosome pair. Nucleolus organizer regions have frequently expanded on sex chromosomes, most probably by ectopic recombination. Our data suggest the involvement of sex chromosome-linked NORs in achiasmatic pairing.

Oscillating waves of Fox, Cyclin and CDK gene expression indicate unique spatiotemporal control of cell cycling during nervous system development in onychophorans

Forkhead box (Fox) genes code for a class of transcription factors with many different fundamental functions in animal development including cell cycle control. Other important factors of cell cycle control are Cyclins and Cyclin-dependent kinases (CDKs). Here we report on the oscillating expression of three Fox genes, FoxM, FoxN14 (jumeaux) and FoxN23 (Checkpoint suppressor like-1), Cyclins and CDKs in an onychophoran, a representative of a relatively small group of animals that are closely related to the arthropods. Expression of these genes is in the form of several waves that start as dot-like domains in the center of each segment and then transform into concentric rings that run towards the periphery of the segments. This oscillating gene expression, however, occurs exclusively along the anterior-posterior body axis in the tissue ventral to the base of the appendages, a region where the central nervous system and the enigmatic ventral and preventral organs of the onychophoran develop. We suggest that the oscillating gene expression and the resulting waves of expression we report are likely correlated with cell cycle control during the development of the onychophoran nervous system. This intriguing patterning appears to be unique for onychophorans as it is not found in any of the arthropods we also investigated in this study, and is likely correlated with the slow embryonic development of onychophorans compared to arthropods.

Striking parallels between dorsoventral patterning in Drosophila and Gryllus reveal a complex evolutionary history behind a model gene regulatory network

Dorsoventral pattering relies on Toll and BMP signalling in all insects studied so far, with variations in the relative contributions of both pathways. Drosophila and the beetle Tribolium share extensive dependence on Toll, while representatives of more distantly related lineages like the wasp Nasonia and bug Oncopeltus rely more strongly on BMP signalling. Here, we show that in the cricket Gryllus bimaculatus, an evolutionarily distant outgroup, Toll has, like in Drosophila, a direct patterning role for the ventral half of the embryo. In addition, Toll polarises BMP signalling, although this does not involve the conserved BMP inhibitor Sog/Chordin. Finally, Toll activation relies on ovarian patterning mechanisms with striking similarity to Drosophila. Our data suggest two surprising hypotheses: (1) that Toll's patterning function in Gryllus and Drosophila is the result of convergent evolution or (2) a Drosophila-like system arose early in insect evolution and was extensively altered in multiple independent lineages.

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

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

A chromosome-level genome of the spider Trichonephila antipodiana reveals the genetic basis of its polyphagy and evidence of an ancient whole-genome duplication event

Background

The spider Trichonephila antipodiana (Araneidae), commonly known as the batik golden web spider, preys on arthropods with body sizes ranging from ∼2 mm in length to insects larger than itself (>20‒50 mm), indicating its polyphagy and strong dietary detoxification abilities. Although it has been reported that an ancient whole-genome duplication event occurred in spiders, lack of a high-quality genome has limited characterization of this event.

Results

We present a chromosome‐level T. antipodiana genome constructed on the basis of PacBio and Hi-C sequencing. The assembled genome is 2.29 Gb in size with a scaffold N50 of 172.89 Mb. Hi‐C scaffolding assigned 98.5% of the bases to 13 pseudo-chromosomes, and BUSCO completeness analysis revealed that the assembly included 94.8% of the complete arthropod universal single-copy orthologs (n = 1,066). Repetitive elements account for 59.21% of the genome. We predicted 19,001 protein-coding genes, of which 96.78% were supported by transcriptome-based evidence and 96.32% matched protein records in the UniProt database. The genome also shows substantial expansions in several detoxification-associated gene families, including cytochrome P450 mono-oxygenases, carboxyl/cholinesterases, glutathione-S-transferases, and ATP-binding cassette transporters, reflecting the possible genomic basis of polyphagy. Further analysis of the T. antipodiana genome architecture reveals an ancient whole-genome duplication event, based on 2 lines of evidence: (i) large-scale duplications from inter-chromosome synteny analysis and (ii) duplicated clusters of Hox genes.

Conclusions

The high-quality T. antipodiana genome represents a valuable resource for spider research and provides insights into this species’ adaptation to the environment.

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.

A conserved role for arrow in posterior axis patterning across Arthropoda

Segmentation is a key characteristic of Arthropoda that is linked to the evolutionary success of this lineage. It has previously been shown in both vertebrates and short germ insects that posterior segmentation requires canonical Wnt (cWnt) signaling, which maintains the expression of Caudal and the posterior growth zone; disruption of cWnt signaling incurs posterior truncations in these lineages due to the loss of the tail bud. However, comparable datasets for Wnt signaling are limited outside of holometabolous insects, due to incomparable phenotypic spectra and inefficacy of gene misexpression methods in certain model species. We applied RNA interference (RNAi) against the Wnt co-receptor arrow (arr), a key member of the cWnt signaling pathway in holometabolous insects and vertebrates, to examine posterior axis elongation of the cobweb spider Parasteatoda tepidariorum (short germ embryogenesis; one Wnt8 homolog), the cricket Gryllus bimaculatus (intermediate germ; one Wnt8 homolog), and the milkweed bug Oncopeltus fasciatus (short germ; two Wnt8 homologs). Knockdown of arr in insects resulted in posterior truncations affecting the gnathos through the abdomen in O. fasciatus, whereas posterior truncations only affected the T3 segment through the abdomen in G. bimaculatus. Spider embryos with disrupted arr expression exhibited defects along the entire axis, including segmentation defects throughout the germband. RNA-Seq-based differential gene expression analysis of severe Ptep-arr loss-of-function phenotypes at two developmental stages was used to confirm that knockdown of Ptep-arr results in systemic disruption of the Wnt pathway. Intriguingly, we found that knockdown of arr did not abrogate Wnt8 expression in any of the three species, with cad expression additionally retained in severe loss-of-function phenotypes in the cricket and the spider. Together with data from a holometabolous insect, our results suggest that cWnt signaling is not required for maintenance of Wnt8 expression across Arthropoda. These outcomes underscore the diagnostic power of differential gene expression analyses in characterizing catastrophic phenotypes in emerging model species.

Characterization of ecdysteroid biosynthesis in the pond wolf spider, Pardosa pseudoannulata

Ecdysteroids, as the key growth hormones, regulate moulting, metamorphosis and reproduction in arthropods. Ecdysteroid biosynthesis is catalysed by a series of cytochrome P450 monooxygenases (CYP450s) encoded by Halloween genes, including spook (spo), phantom (phm), disembodied (dib), shadow (sad) and shade (shd). The ecdysteroid biosynthesis in insects is clear with 20-hydroxyecdysone (20E) as the main ecdysteroid. However, the information on the major ecdysteroids in arachnids is limited. In this study, Halloween genes spo, dib, sad and shd, but not phm, were identified in the pond wolf spider, Pardosa pseudoannulata. Phylogenetic analysis grouped arachnid and insect Halloween gene products into two CYP450 clades, the CYP2 clan (spo and phm) and the mitochondrial clan (dib, sad, and shd). In P. pseudoannulata, the temporal expression profile of the four Halloween genes in concurrence with spiderling moulting with steady increase in the course of the 2nd instar followed by a rapid dropdown once moulting was completed. Spatially, the four Halloween genes were highly expressed in spiderling abdomen and in the ovaries of female adults. In parallel, ponasterone A (PA), but not 20E, was detected by LC-MS/MS analysis in P. pseudoannulata, and it was demonstrated as a functional ecdysteroid in the spider by accelerating of moulting with PA addition. The present study revealed the different ecdysteroid biosynthesis pathways in spiders and insects.

Phylogenomic Resolution of Sea Spider Diversification through Integration of Multiple Data Classes

Despite significant advances in invertebrate phylogenomics over the past decade, the higher-level phylogeny of Pycnogonida (sea spiders) remains elusive. Due to the inaccessibility of some small-bodied lineages, few phylogenetic studies have sampled all sea spider families. Previous efforts based on a handful of genes have yielded unstable tree topologies. Here, we inferred the relationships of 89 sea spider species using targeted capture of the mitochondrial genome, 56 conserved exons, 101 ultraconserved elements, and 3 nuclear ribosomal genes. We inferred molecular divergence times by integrating morphological data for fossil species to calibrate 15 nodes in the arthropod tree of life. This integration of data classes resolved the basal topology of sea spiders with high support. The enigmatic family Austrodecidae was resolved as the sister group to the remaining Pycnogonida and the small-bodied family Rhynchothoracidae as the sister group of the robust-bodied family Pycnogonidae. Molecular divergence time estimation recovered a basal divergence of crown group sea spiders in the Ordovician. Comparison of diversification dynamics with other marine invertebrate taxa that originated in the Paleozoic suggests that sea spiders and some crustacean groups exhibit resilience to mass extinction episodes, relative to mollusk and echinoderm lineages.

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.

Spider Diversification Through Space and Time

Spiders (Araneae) make up a remarkably diverse lineage of predators that have successfully colonized most terrestrial ecosystems. All spiders produce silk, and many species use it to build capture webs with an extraordinary diversity of forms. Spider diversity is distributed in a highly uneven fashion across lineages. This strong imbalance in species richness has led to several causal hypotheses, such as codiversification with insects, key innovations in silk structure and web architecture, and loss of foraging webs. Recent advances in spider phylogenetics have allowed testing of some of these hypotheses, but results are often contradictory, highlighting the need to consider additional drivers of spider diversification. The spatial and historical patterns of diversity and diversification remain contentious. Comparative analyses of spider diversification will advance only if we continue to make progress with studies of species diversity, distribution, and phenotypic traits, together with finer-scale phylogenies and genomic data.

Chromosome-level reference genome of the European wasp spider Argiope bruennichi: a resource for studies on range expansion and evolutionary adaptation

BACKGROUND

Argiope bruennichi, the European wasp spider, has been investigated intensively as a focal species for studies on sexual selection, chemical communication, and the dynamics of rapid range expansion at a behavioral and genetic level. However, the lack of a reference genome has limited insights into the genetic basis for these phenomena. Therefore, we assembled a high-quality chromosome-level reference genome of the European wasp spider as a tool for more in-depth future studies.

FINDINGS

We generated, de novo, a 1.67 Gb genome assembly of A. bruennichi using 21.8× Pacific Biosciences sequencing, polished with 19.8× Illumina paired-end sequencing data, and proximity ligation (Hi-C)-based scaffolding. This resulted in an N50 scaffold size of 124 Mb and an N50 contig size of 288 kb. We found 98.4% of the genome to be contained in 13 scaffolds, fitting the expected number of chromosomes (n = 13). Analyses showed the presence of 91.1% of complete arthropod BUSCOs, indicating a high-quality assembly.

CONCLUSIONS

We present the first chromosome-level genome assembly in the order Araneae. With this genomic resource, we open the door for more precise and informative studies on evolution and adaptation not only in A. bruennichi but also in arachnids overall, shedding light on questions such as the genomic architecture of traits, whole-genome duplication, and the genomic mechanisms behind silk and venom evolution.

Insights into how development and life-history dynamics shape the evolution of venom

Venomous animals are a striking example of the convergent evolution of a complex trait. These animals have independently evolved an apparatus that synthesizes, stores, and secretes a mixture of toxic compounds to the target animal through the infliction of a wound. Among these distantly related animals, some can modulate and compartmentalize functionally distinct venoms related to predation and defense. A process to separate distinct venoms can occur within and across complex life cycles as well as more streamlined ontogenies, depending on their life-history requirements. Moreover, the morphological and cellular complexity of the venom apparatus likely facilitates the functional diversity of venom deployed within a given life stage. Intersexual variation of venoms has also evolved further contributing to the massive diversity of toxic compounds characterized in these animals. These changes in the biochemical phenotype of venom can directly affect the fitness of these animals, having important implications in their diet, behavior, and mating biology. In this review, we explore the current literature that is unraveling the temporal dynamics of the venom system that are required by these animals to meet their ecological functions. These recent findings have important consequences in understanding the evolution and development of a convergent complex trait and its organismal and ecological implications.

Evolution of MicroRNA Biogenesis Genes in the Sterlet (Acipenser ruthenus) and Other Polyploid Vertebrates

MicroRNAs play a crucial role in eukaryotic gene regulation. For a long time, only little was known about microRNA-based gene regulatory mechanisms in polyploid animal genomes due to difficulties of polyploid genome assembly. However, in recent years, several polyploid genomes of fish, amphibian, and even invertebrate species have been sequenced and assembled. Here we investigated several key microRNA-associated genes in the recently sequenced sterlet (Acipenser ruthenus) genome, whose lineage has undergone a whole genome duplication around 180 MYA. We show that two paralogs of drosha, dgcr8, xpo1, and xpo5 as well as most ago genes have been retained after the acipenserid-specific whole genome duplication, while ago1 and ago3 genes have lost one paralog. While most diploid vertebrates possess only a single copy of dicer1, we strikingly found four paralogs of this gene in the sterlet genome, derived from a tandem segmental duplication that occurred prior to the last whole genome duplication. ago1,3,4 and exportins1,5 look to be prone to additional segment duplications producing up to four-five paralog copies in ray-finned fishes. We demonstrate for the first time exon microsatellite amplification in the acipenserid drosha2 gene, resulting in a highly variable protein product, which may indicate sub- or neofunctionalization. Paralogous copies of most microRNA metabolism genes exhibit different expression profiles in various tissues and remain functional despite the rediploidization process. Subfunctionalization of microRNA processing gene paralogs may be beneficial for different pathways of microRNA metabolism. Genetic variability of microRNA processing genes may represent a substrate for natural selection, and, by increasing genetic plasticity, could facilitate adaptations to changing environments.

Diversity and evolution of the P450 family in arthropods

The P450 family (CYP genes) of arthropods encodes diverse enzymes involved in the metabolism of foreign compounds and in essential endocrine or ecophysiological functions. The P450 sequences (CYPome) from 40 arthropod species were manually curated, including 31 complete CYPomes, and a maximum likelihood phylogeny of nearly 3000 sequences is presented. Arthropod CYPomes are assembled from members of six CYP clans of variable size, the CYP2, CYP3, CYP4 and mitochondrial clans, as well as the CYP20 and CYP16 clans that are not found in Neoptera. CYPome sizes vary from two dozen genes in some parasitic species to over 200 in species as diverse as collembolans or ticks. CYPomes are comprised of few CYP families with many genes and many CYP families with few genes, and this distribution is the result of dynamic birth and death processes. Lineage-specific expansions or blooms are found throughout the phylogeny and often result in genomic clusters that appear to form a reservoir of catalytic diversity maintained as heritable units. Among the many P450s with physiological functions, six CYP families are involved in ecdysteroid metabolism. However, five so-called Halloween genes are not universally represented and do not constitute the unique pathway of ecdysteroid biosynthesis. The diversity of arthropod CYPomes has only partially been uncovered to date and many P450s with physiological functions regulating the synthesis and degradation of endogenous signal molecules (including ecdysteroids) and semiochemicals (including pheromones and defense chemicals) remain to be discovered. Sequence diversity of arthropod P450s is extreme, and P450 sequences lacking the universally conserved Cys ligand to the heme have evolved several times. A better understanding of P450 evolution is needed to discern the relative contributions of stochastic processes and adaptive processes in shaping the size and diversity of CYPomes.

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.