Homeobox Gene Duplication and Divergence in Arachnids

An atlas of spider development at single-cell resolution provides new insights into arthropod embryogenesis

Spiders are a diverse order of chelicerates that diverged from other arthropods over 500 million years ago. Research on spider embryogenesis, particularly studies using the common house spider Parasteatoda tepidariorum, has made important contributions to understanding the evolution of animal development, including axis formation, segmentation, and patterning. However, we lack knowledge about the cells that build spider embryos, their gene expression profiles and fate. Single-cell transcriptomic analyses have been revolutionary in describing these complex landscapes of cellular genetics in a range of animals. Therefore, we carried out single-cell RNA sequencing of P. tepidariorum embryos at stages 7, 8 and 9, which encompass the establishment and patterning of the body plan, and initial differentiation of many tissues and organs. We identified 20 cell clusters, from 18.5 k cells, which were marked by many developmental toolkit genes, as well as a plethora of genes not previously investigated. We found differences in the cell cycle transcriptional signatures, suggestive of different proliferation dynamics, which related to distinctions between endodermal and some mesodermal clusters, compared with ectodermal clusters. We identified many Hox genes as markers of cell clusters, and Hox gene ohnologs were often present in different clusters. This provided additional evidence of sub- and/or neo-functionalisation of these important developmental genes after the whole genome duplication in an arachnopulmonate ancestor (spiders, scorpions, and related orders). We also examined the spatial expression of marker genes for each cluster to generate a comprehensive cell atlas of these embryonic stages. This revealed new insights into the cellular basis and genetic regulation of head patterning, hematopoiesis, limb development, gut development, and posterior segmentation. This atlas will serve as a platform for future analysis of spider cell specification and fate, and studying the evolution of these processes among animals at cellular resolution.

Insight into the adaptive role of arachnid genome-wide duplication through chromosome-level genome assembly of the Western black widow spider

Although spiders are one of the most diverse groups of arthropods, the genetic architecture of their evolutionary adaptations is largely unknown. Specifically, ancient genome-wide duplication occurring during arachnid evolution ~450 mya resulted in a vast assembly of gene families, yet the extent to which selection has shaped this variation is understudied. To aid in comparative genome sequence analyses, we provide a chromosome-level genome of the Western black widow spider (Latrodectus hesperus)-a focus due to its silk properties, venom applications, and as a model for urban adaptation. We used long-read and Hi-C sequencing data, combined with transcriptomes, to assemble 14 chromosomes in a 1.46 Gb genome, with 38,393 genes annotated, and a BUSCO score of 95.3%. Our analyses identified high repetitive gene content and heterozygosity, consistent with other spider genomes, which has led to challenges in genome characterization. Our comparative evolutionary analyses of eight genomes available for species within the Araneoidea group (orb weavers and their descendants) identified 1,827 single-copy orthologs. Of these, 155 exhibit significant positive selection primarily associated with developmental genes, and with traits linked to sensory perception. These results support the hypothesis that several traits unique to spiders emerged from the adaptive evolution of ohnologs-or retained ancestrally duplicated genes-from ancient genome-wide duplication. These comparative spider genome analyses can serve as a model to understand how positive selection continually shapes ancestral duplications in generating novel traits today within and between diverse taxonomic groups.

Single-cell RNA sequencing of mid-to-late stage spider embryos: new insights into spider development

BACKGROUND

The common house spider Parasteatoda tepidariorum represents an emerging new model organism of arthropod evolutionary and developmental (EvoDevo) studies. Recent technical advances have resulted in the first single-cell sequencing (SCS) data on this species allowing deeper insights to be gained into its early development, but mid-to-late stage embryos were not included in these pioneering studies.

RESULTS

Therefore, we performed SCS on mid-to-late stage embryos of Parasteatoda and characterized resulting cell clusters by means of in-silico analysis (comparison of key markers of each cluster with previously published information on these genes). In-silico prediction of the nature of each cluster was then tested/verified by means of additional in-situ hybridization experiments with additional markers of each cluster.

CONCLUSIONS

Our data show that SCS data reliably group cells with similar genetic fingerprints into more or less distinct clusters, and thus allows identification of developing cell types on a broader level, such as the distinction of ectodermal, mesodermal and endodermal cell lineages, as well as the identification of distinct developing tissues such as subtypes of nervous tissue cells, the developing heart, or the ventral sulcus (VS). In comparison with recent other SCS studies on the same species, our data represent later developmental stages, and thus provide insights into different stages of developing cell types and tissues such as differentiating neurons and the VS that are only present at these later stages.

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.

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

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

Evolution of the Spider Homeobox Gene Repertoire by Tandem and Whole Genome Duplication

Gene duplication generates new genetic material that can contribute to the evolution of gene regulatory networks and phenotypes. Duplicated genes can undergo subfunctionalization to partition ancestral functions and/or neofunctionalization to assume a new function. We previously found there had been a whole genome duplication (WGD) in an ancestor of arachnopulmonates, the lineage including spiders and scorpions but excluding other arachnids like mites, ticks, and harvestmen. This WGD was evidenced by many duplicated homeobox genes, including two Hox clusters, in spiders. However, it was unclear which homeobox paralogues originated by WGD versus smaller-scale events such as tandem duplications. Understanding this is a key to determining the contribution of the WGD to arachnopulmonate genome evolution. Here we characterized the distribution of duplicated homeobox genes across eight chromosome-level spider genomes. We found that most duplicated homeobox genes in spiders are consistent with an origin by WGD. We also found two copies of conserved homeobox gene clusters, including the Hox, NK, HRO, Irx, and SINE clusters, in all eight species. Consistently, we observed one copy of each cluster was degenerated in terms of gene content and organization while the other remained more intact. Focussing on the NK cluster, we found evidence for regulatory subfunctionalization between the duplicated NK genes in the spider Parasteatoda tepidariorum compared to their single-copy orthologues in the harvestman Phalangium opilio. Our study provides new insights into the relative contributions of multiple modes of duplication to the homeobox gene repertoire during the evolution of spiders and the function of NK genes.

The Impact of Whole Genome Duplication on the Evolution of the Arachnids

The proliferation of genomic resources for Chelicerata in the past 10 years has revealed that the evolution of chelicerate genomes is more dynamic than previously thought, with multiple waves of ancient whole genome duplications affecting separate lineages. Such duplication events are fascinating from the perspective of evolutionary history because the burst of new gene copies associated with genome duplications facilitates the acquisition of new gene functions (neofunctionalization), which may in turn lead to morphological novelties and spur net diversification. While neofunctionalization has been invoked in several contexts with respect to the success and diversity of spiders, the overall impact of whole genome duplications on chelicerate evolution and development remains imperfectly understood. The purpose of this review is to examine critically the role of whole genome duplication on the diversification of the extant arachnid orders, as well as assess functional datasets for evidence of subfunctionalization or neofunctionalization in chelicerates. This examination focuses on functional data from two focal model taxa: the spider Parasteatoda tepidariorum, which exhibits evidence for an ancient duplication, and the harvestman Phalangium opilio, which exhibits an unduplicated genome. I show that there is no evidence that taxa with genome duplications are more successful than taxa with unduplicated genomes. I contend that evidence for sub- or neofunctionalization of duplicated developmental patterning genes in spiders is indirect or fragmentary at present, despite the appeal of this postulate for explaining the success of groups like spiders. Available expression data suggest that the condition of duplicated Hox modules may have played a role in promoting body plan disparity in the posterior tagma of some orders, such as spiders and scorpions, but functional data substantiating this postulate are critically missing. Spatiotemporal dynamics of duplicated transcription factors in spiders may represent cases of developmental system drift, rather than neofunctionalization. Developmental system drift may represent an important, but overlooked, null hypothesis for studies of paralogs in chelicerate developmental biology. To distinguish between subfunctionalization, neofunctionalization, and developmental system drift, concomitant establishment of comparative functional datasets from taxa exhibiting the genome duplication, as well as those that lack the paralogy, is sorely needed.

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

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

Expression of netrin and its receptors uncoordinated-5 and frazzled in arthropods and onychophorans suggests conserved and diverged functions in neuronal pathfinding and synaptogenesis

BACKGROUND

Development of the nervous system and the correct connection of nerve cells require coordinated axonal pathfinding through an extracellular matrix. Outgrowing axons exhibit directional growth toward or away from external guidance cues such as Netrin. Guidance cues can be detected by growth cones that are located at the end of growing axons through membrane-bound receptors such as Uncoordianted-5 and Frazzled. Binding of Netrin causes reformation of the cytoskeleton and growth of the axon toward (or away from) the source of Netrin production.

RESULTS

Here, we investigate the embryonic mRNA expression patterns of netrin genes and their potential receptors, uncoordinated-5 and frazzled in arthropod species that cover all main branches of Arthropoda, that is, Pancrustacea, Myriapoda, and Chelicerata. We also studied the expression patterns in a closely related outgroup species, the onychophoran Euperipatoides kanangrensis, and provide data on expression profiles of these genes in larval tissues of the fly Drosophila melanogaster including the brain and the imaginal disks.

CONCLUSION

Our data reveal conserved and diverged aspects of neuronal guidance in Drosophila with respect to the other investigated species and suggest a conserved function in nervous system patterning of the developing appendages.

Diversity, duplication, and genomic organization of homeobox genes in Lepidoptera

Homeobox genes encode transcription factors with essential roles in patterning and cell fate in developing animal embryos. Many homeobox genes, including Hox and NK genes, are arranged in gene clusters, a feature likely related to transcriptional control. Sparse taxon sampling and fragmentary genome assemblies mean that little is known about the dynamics of homeobox gene evolution across Lepidoptera or about how changes in homeobox gene number and organization relate to diversity in this large order of insects. Here we analyze an extensive data set of high-quality genomes to characterize the number and organization of all homeobox genes in 123 species of Lepidoptera from 23 taxonomic families. We find most Lepidoptera have around 100 homeobox loci, including an unusual Hox gene cluster in which the lab gene is repositioned and the ro gene is next to pb A topologically associating domain spans much of the gene cluster, suggesting deep regulatory conservation of the Hox cluster arrangement in this insect order. Most Lepidoptera have four Shx genes, divergent zen-derived loci, but these loci underwent dramatic duplication in several lineages, with some moths having over 165 homeobox loci in the Hox gene cluster; this expansion is associated with local LINE element density. In contrast, the NK gene cluster content is more stable, although there are differences in organization compared with other insects, as well as major rearrangements within butterflies. Our analysis represents the first description of homeobox gene content across the order Lepidoptera, exemplifying the potential of newly generated genome assemblies for understanding genome and gene family evolution.

Hox genes in spiders: Their significance for development and evolution

Hox genes are known for their role in the specification of typical body plan features in animals. Evolutionary changes in Hox gene function are believed to be involved in the emergence of the diverse body plans we observe in animals today. Spiders share many body plan features with other arthropods, but also have numerous unique traits of their own. Studies of spider Hox genes have already provided insights into evolutionarily conserved and derived features of the spider body plan and their genetic basis. However, many aspects of Hox gene biology have been insufficiently studied in spiders so far. In this review, we highlight previous comparative studies of Hox genes in spiders and their significance for our understanding of the evolution of the spider body plan. We also identify aspects of Hox gene biology that need to be studied in greater detail. Many spider Hox genes have not been investigated beyond their mRNA expression patterns, and the role of Hox genes with apparently plesiomorphic or dual functions, like ftz and Hox3 is still unclear. Spiders have a duplicated Hox gene cluster, but possible sub- or neofunctionalisation of duplicates have not yet been studied systematically. Future research should therefore focus on these issues, in addition to the role of Polycomb and trithorax-mediated regulation, the identification of regulatory regions, cofactors or spider-specific target genes, and the significance of non-coding RNAs transcribed from within the Hox cluster and even from the antisense strand of particular Hox genes.

The lower Cambrian lobopodian Cardiodictyon resolves the origin of euarthropod brains

For more than a century, the origin and evolution of the arthropod head and brain have eluded a unifying rationale reconciling divergent morphologies and phylogenetic relationships. Here, clarification is provided by the fossilized nervous system of the lower Cambrian lobopodian Cardiodictyon catenulum, which reveals an unsegmented head and brain comprising three cephalic domains, distinct from the metameric ventral nervous system serving its appendicular trunk. Each domain aligns with one of three components of the foregut and with a pair of head appendages. Morphological correspondences with stem group arthropods and alignments of homologous gene expression patterns with those of extant panarthropods demonstrate that cephalic domains of C. catenulum predate the evolution of the euarthropod head yet correspond to neuromeres defining brains of living chelicerates and mandibulates.

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.

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.

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.

Emergence of distinct syntenic density regimes is associated with early metazoan genomic transitions

Background

Animal genomes are strikingly conserved in terms of local gene order (microsynteny). While some of these microsyntenies have been shown to be coregulated or to form gene regulatory blocks, the diversity of their genomic and regulatory properties across the metazoan tree of life remains largely unknown.

Results

Our comparative analyses of 49 animal genomes reveal that the largest gains of synteny occurred in the last common ancestor of bilaterians and cnidarians and in that of bilaterians. Depending on their node of emergence, we further show that novel syntenic blocks are characterized by distinct functional compositions (Gene Ontology terms enrichment) and gene density properties, such as high, average and low gene density regimes. This is particularly pronounced among bilaterian novel microsyntenies, most of which fall into high gene density regime associated with higher gene coexpression levels. Conversely, a majority of vertebrate novel microsyntenies display a low gene density regime associated with lower gene coexpression levels.

Conclusions

Our study provides first evidence for evolutionary transitions between different modes of microsyntenic block regulation that coincide with key events of metazoan evolution. Moreover, the microsyntenic profiling strategy and interactive online application (Syntenic Density Browser, available at: http://synteny.csb.univie.ac.at/) we present here can be used to explore regulatory properties of microsyntenic blocks and predict their coexpression in a wide-range of animal genomes.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12864-022-08304-2.

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.

Regulation of Eye Determination and Regionalization in the Spider Parasteatoda tepidariorum

Animal visual systems are enormously diverse, but their development appears to be controlled by a set of conserved retinal determination genes (RDGs). Spiders are particular masters of visual system innovation, and offer an excellent opportunity to study the evolution of animal eyes. Several RDGs have been identified in spider eye primordia, but their interactions and regulation remain unclear. From our knowledge of RDG network regulation in Drosophila melanogaster, we hypothesize that orthologs of Pax6, eyegone, Wnt genes, hh, dpp, and atonal could play important roles in controlling eye development in spiders. We analyzed the expression of these genes in developing embryos of the spider Parasteatodatepidariorum, both independently and in relation to the eye primordia, marked using probes for the RDG sine oculis. Our results support conserved roles for Wnt genes in restricting the size and position of the eye field, as well as for atonal initiating photoreceptor differentiation. However, we found no strong evidence for an upstream role of Pax6 in eye development, despite its label as a master regulator of animal eye development; nor do eyg, hh or dpp compensate for the absence of Pax6. Conversely, our results indicate that hh may work with Wnt signaling to restrict eye growth, a role similar to that of Sonichedgehog (Shh) in vertebrates.

De novo transcriptome sequencing of the northern fowl mite, Ornithonyssus sylviarum, shed light on parasitiform poultry mites evolution and its chemoreceptor repertoires

The northern fowl mite (NFM), Ornithonyssus sylviarum, and the poultry red mite (PRM), Dermanyssus gallinae, are the most serious pests of poultry, both of which have an expanding global prevalence. Research on NFM has been constrained by a lack of genomic and transcriptomic data. Here, we report and analyze the first global transcriptome data across all mite live stages and sexes. A total of 28,999 unigenes were assembled, of which 19,750 (68.10%) were annotated using seven functional databases. The biological function of these unigenes was classified using the GO, KOG, and KEGG databases. To gain insight into the chemosensory receptor-based system of parasitiform mites, we furthermore assessed the gene repertoire of gustatory receptors (GRs) and ionotropic receptors (IRs), both of which encode putative ligand-gated ion channel proteins. While these receptors are well characterized in insect model species, our understanding of chemosensory detection in mites and ticks is in its infancy. To address this paucity of data, we identified 9 IR/iGluRs and 2 GRs genes by analyzing transcriptome data in the NFM, while 9 GRs and 41 IR/iGluRs genes were annotated in the PRM genome. Taken together, the transcriptomic and genomic annotation of these two species provide a valuable reference for studies of parasitiform mites and also help to understand how chemosensory gene family expansion/contraction events may have been reshaped by an obligate parasitic lifestyle compared with their free-living closest relatives. Future studies should include additional species to validate this observation and functional characterization of the identified proteins as a step forward in identifying tools for controlling these poultry pests.

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.

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.

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

Horizontally acquired cysteine synthase genes undergo functional divergence in lepidopteran herbivores

Horizontal gene transfer (HGT) plays an important role in evolutionary processes as organisms adapt to their environments, and now cases of gene duplication after HGT in eukaryotes are emerging at an increasing rate. However, the fate and roles of the duplicated genes over time in eukaryotes remain unclear. Here we conducted a comprehensive analysis of the evolution of cysteine synthase (CYS) in lepidopteran insects. Our results indicate that HGT-derived CYS genes are widespread and have undergone duplication following horizontal transfer in many lepidopteran insects. Moreover, lepidopteran CYS proteins not only have β-cyanoalanine synthase activity but also possess cysteine synthase activity that is involved in sulfur amino acid biosynthesis. Duplicated CYS genes show marked divergence in gene expression patterns and enzymatic properties, suggesting that they probably have undergone subfunctionalization and/or neofunctionalization in Lepidoptera. The gene transfer of CYS genes and subsequent duplication appears to have facilitated the adaptation of lepidopteran insects to different diets and promoted their ecological diversification. Overall, this study provides valuable insights into the ecological and evolutionary contributions of CYS in lepidopteran insects.

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.

Horizontal Gene Transfer and Gene Duplication of β-Fructofuranosidase Confer Lepidopteran Insects Metabolic Benefits

Horizontal gene transfer (HGT) is a potentially critical source of material for ecological adaptation and the evolution of novel genetic traits. However, reports on posttransfer duplication in organism genomes are lacking, and the evolutionary advantages conferred on the recipient are generally poorly understood. Sucrase plays an important role in insect physiological growth and development. Here, we performed a comprehensive analysis of the evolution of insect β-fructofuranosidase transferred from bacteria via HGT. We found that posttransfer duplications of β-fructofuranosidase were widespread in Lepidoptera and sporadic occurrences of β-fructofuranosidase were found in Coleoptera and Hymenoptera. β-fructofuranosidase genes often undergo modifications, such as gene duplication, differential gene loss, and changes in mutation rates. Lepidopteran β-fructofuranosidase gene (SUC) clusters showed marked divergence in gene expression patterns and enzymatic properties in Bombyx mori (moth) and Papilio xuthus (butterfly). We generated SUC1 mutations in B. mori using CRISPR/Cas9 to thoroughly examine the physiological function of SUC. BmSUC1 mutant larvae were viable but displayed delayed growth and reduced sucrase activities that included susceptibility to the sugar mimic alkaloid found in high concentrations in mulberry. BmSUC1 served as a critical sucrase and supported metabolic homeostasis in the larval midgut and silk gland, suggesting that gene transfer of β-fructofuranosidase enhanced the digestive and metabolic adaptation of lepidopteran insects. These findings highlight not only the universal function of β-fructofuranosidase with a link to the maintenance of carbohydrate metabolism but also an underexplored function in the silk gland. This study expands our knowledge of posttransfer duplication and subsequent functional diversification in the adaptive evolution and lineage-specific adaptation of organisms.

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.

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.

Arachnid monophyly: Morphological, palaeontological and molecular support for a single terrestrialization within Chelicerata

The majority of extant arachnids are terrestrial, but other chelicerates are generally aquatic, including horseshoe crabs, sea spiders, and the extinct eurypterids. It is necessary to determine whether arachnids are exclusively descended from a single common ancestor (monophyly), because only that relationship is compatible with one land colonisation in chelicerate evolutionary history. Some studies have cast doubt on arachnid monophyly and recast the origins of their terrestrialization. These include some phylogenomic analyses placing horseshoe crabs within Arachnida, and from aquatic Palaeozoic stem-group scorpions. Here, we evaluate the possibility of arachnid monophyly by considering morphology, fossils and molecules holistically. We argue arachnid monophyly obviates the need to posit reacquisition/retention of aquatic characters such as gnathobasic feeding and book gills without trabeculae from terrestrial ancestors in horseshoe crabs, and that the scorpion total-group contains few aquatic taxa. We built a matrix composed of 200 slowly-evolving genes and re-analysed two published molecular datasets. We retrieved arachnid monophyly where other studies did not - highlighting the difficulty of resolving chelicerate relationships from current molecular data. As such, we consider arachnid monophyly the best-supported hypothesis. Finally, we inferred that arachnids terrestrialized during the Cambrian-Ordovician using the slow-evolving molecular matrix, in agreement with recent analyses.

Hedgehog signaling controls segmentation dynamics and diversity via msx1 in a spider embryo

Hedgehog (Hh) signaling plays fundamental roles in animal body patterning. Understanding its mechanistic complexity requires simple tractable systems that can be used for these studies. In the early spider embryo, Hh signaling mediates the formation of overall anterior-posterior polarity, yet it remains unclear what mechanisms link the initial Hh signaling activity with body axis segmentation, in which distinct periodic stripe-forming dynamics occur depending on the body region. We performed genome-wide searches for genes that transcriptionally respond to altered states of Hh signaling. Characterization of genes negatively regulated by Hh signaling suggested that msx1, encoding a conserved transcription factor, functions as a key segmentation gene. Knockdown of msx1 prevented all dynamic processes causing spatial repetition of stripes, including temporally repetitive oscillations and bi-splitting, and temporally nonrepetitive tri-splitting. Thus, Hh signaling controls segmentation dynamics and diversity via msx1 These genome-wide data from an invertebrate illuminate novel mechanistic features of Hh-based patterning.

Genomic resources and toolkits for developmental study of whip spiders (Amblypygi) provide insights into arachnid genome evolution and antenniform leg patterning

BACKGROUND

The resurgence of interest in the comparative developmental study of chelicerates has led to important insights, such as the discovery of a genome duplication shared by spiders and scorpions, inferred to have occurred in the most recent common ancestor of Arachnopulmonata (a clade comprising the five arachnid orders that bear book lungs). Nonetheless, several arachnid groups remain understudied in the context of development and genomics, such as the order Amblypygi (whip spiders). The phylogenetic position of Amblypygi in Arachnopulmonata posits them as an interesting group to test the incidence of the proposed genome duplication in the common ancestor of Arachnopulmonata, as well as the degree of retention of duplicates over 450 Myr. Moreover, whip spiders have their first pair of walking legs elongated and modified into sensory appendages (a convergence with the antennae of mandibulates), but the genetic patterning of these antenniform legs has never been investigated.

RESULTS

We established genomic resources and protocols for cultivation of embryos and gene expression assays by in situ hybridization to study the development of the whip spider Phrynus marginemaculatus. Using embryonic transcriptomes from three species of Amblypygi, we show that the ancestral whip spider exhibited duplications of all ten Hox genes. We deploy these resources to show that paralogs of the leg gap genes dachshund and homothorax retain arachnopulmonate-specific expression patterns in P. marginemaculatus. We characterize the expression of leg gap genes Distal-less, dachshund-1/2 and homothorax-1/2 in the embryonic antenniform leg and other appendages, and provide evidence that allometry, and by extension the antenniform leg fate, is specified early in embryogenesis.

CONCLUSION

This study is the first step in establishing P. marginemaculatus as a chelicerate model for modern evolutionary developmental study, and provides the first resources sampling whip spiders for comparative genomics. Our results suggest that Amblypygi share a genome duplication with spiders and scorpions, and set up a framework to study the genetic specification of antenniform legs. Future efforts to study whip spider development must emphasize the development of tools for functional experiments in P. marginemaculatus.

Temperature fluctuations during embryonic development implicated in a naturally occurring instance of abnormal spinnerets in the spider Australomimetus maculosus (Araneae, Mimetidae)

We record developmental abnormalities of the spinnerets in a field-collected adult male specimen of Australomimetus maculosus. These include (1) a supernumerary right posterior lateral spinneret (PLS), (2) ectopic piriform silk gland spigots and tartipores on the left PLS that are normally restricted to anterior lateral spinnerets (ALSs), and (3) what appear to be ectopic ALS sensilla on the left posterior median spinneret (PMS). Published results of teratological experiments and climate data for the collection site indicate that fluctuating sub- and supra-optimal temperatures during embryogenesis may have been responsible for these anomalies. This specimen thus supports the view that spinneret abnormalities, among other aberrations, may be induced when embryos of entelegyne spiders are exposed to fluctuations between high and low temperatures, whether in the laboratory or, as here, in nature. To our knowledge, the ectopic structures seen on the left PLS and left PMS have not been observed previously. Their locations are consistent with a hypothesis by which only the lateral portion of the araneomorph ALS is serially homologous to the PLS, while the remainder of the ALS, along with the colulus/cribellum, is homologous to the PMS.

The gene-rich genome of the scallop Pecten maximus

BACKGROUND

The king scallop, Pecten maximus, is distributed in shallow waters along the Atlantic coast of Europe. It forms the basis of a valuable commercial fishery and plays a key role in coastal ecosystems and food webs. Like other filter feeding bivalves it can accumulate potent phytotoxins, to which it has evolved some immunity. The molecular origins of this immunity are of interest to evolutionary biologists, pharmaceutical companies, and fisheries management.

FINDINGS

Here we report the genome assembly of this species, conducted as part of the Wellcome Sanger 25 Genomes Project. This genome was assembled from PacBio reads and scaffolded with 10X Chromium and Hi-C data. Its 3,983 scaffolds have an N50 of 44.8 Mb (longest scaffold 60.1 Mb), with 92% of the assembly sequence contained in 19 scaffolds, corresponding to the 19 chromosomes found in this species. The total assembly spans 918.3 Mb and is the best-scaffolded marine bivalve genome published to date, exhibiting 95.5% recovery of the metazoan BUSCO set. Gene annotation resulted in 67,741 gene models. Analysis of gene content revealed large numbers of gene duplicates, as previously seen in bivalves, with little gene loss, in comparison with the sequenced genomes of other marine bivalve species.

CONCLUSIONS

The genome assembly of P. maximus and its annotated gene set provide a high-quality platform for studies on such disparate topics as shell biomineralization, pigmentation, vision, and resistance to algal toxins. As a result of our findings we highlight the sodium channel gene Nav1, known to confer resistance to saxitoxin and tetrodotoxin, as a candidate for further studies investigating immunity to domoic acid.

The common house spider Parasteatoda tepidariorum

The common house spider Parasteatoda tepidariorum, belonging to the Chelicerata in the phylum Arthropoda, has emerged as an experimental system for studying mechanisms of development from an evolutionary standpoint. In this article, we review the distinct characteristics of P. tepidariorum, the major research questions relevant to this organism, and the available key methods and resources. P. tepidariorum has a relatively short lifecycle and, once mated, periodically lays eggs. The morphogenetic field of the P. tepidariorum embryo is cellular from an early stage and exhibits stepwise symmetry-breaking events and stripe-forming processes that are associated with body axes formation and segmentation, respectively, before reaching the arthropod phylotypic stage. Self-regulatory capabilities of the embryonic field are a prominent feature in P. tepidariorum. The mechanisms and logic underlying the evolvability of heritable patterning systems at the phylum level could be one of the major avenues of research investigated using this animal. The sequenced genome reveals whole genome duplication (WGD) within chelicerates, which offers an invertebrate platform for investigating the potential roles of WGD in animal diversification and evolution. The development and evolution of lineage-specific organs, including the book lungs and the union of spinnerets and silk glands, are attractive subjects of study. Studies using P. tepidariorum can benefit from the use of parental RNA interference, microinjection applications (including cell labeling and embryonic RNA interference), multicolor fluorescence in situ hybridization, and laser ablation as well as rich genomic and transcriptomic resources. These techniques enable functional gene discoveries and the uncovering of cellular and molecular insights.

A Cambrian-Ordovician Terrestrialization of Arachnids

Understanding the temporal context of terrestrialization in chelicerates depends on whether terrestrial groups, the traditional Arachnida, have a single origin and whether or not horseshoe crabs are primitively or secondarily marine. Molecular dating on a phylogenomic tree that recovers arachnid monophyly, constrained by 27 rigorously vetted fossil calibrations, estimates that Arachnida originated during the Cambrian or Ordovician. After the common ancestor colonized the land, the main lineages appear to have rapidly radiated in the Cambrian-Ordovician boundary interval, coinciding with high rates of molecular evolution. The highest rates of arachnid diversification are detected between the Permian and Early Cretaceous. A pattern of ancient divergence estimates for terrestrial arthropod groups in the Cambrian while the oldest fossils are Silurian (seen in both myriapods and arachnids) is mirrored in the molecular and fossil records of land plants. We suggest the discrepancy between molecular and fossil evidence for terrestrialization is likely driven by the extreme sparseness of terrestrial sediments in the rock record before the late Silurian.

six3 acts upstream of foxQ2 in labrum and neural development in the spider Parasteatoda tepidariorum

Anterior patterning in animals is based on a gene regulatory network, which comprises highly conserved transcription factors like six3, pax6 and otx. More recently, foxQ2 was found to be an ancestral component of this network but its regulatory interactions showed evolutionary differences. In most animals, foxQ2 is a downstream target of six3 and knockdown leads to mild or no epidermal phenotypes. In contrast, in the red flour beetle Tribolium castaneum, foxQ2 gained a more prominent role in patterning leading to strong epidermal and brain phenotypes and being required for six3 expression. However, it has remained unclear which of these novel aspects were insect or arthropod specific. Here, we study expression and RNAi phenotype of the single foxQ2 ortholog of the spider Parasteatoda tepidariorum. We find early anterior expression similar to the one of insects. Further, we show an epidermal phenotype in the labrum similar to the insect phenotype. However, our data indicate that foxQ2 is positioned downstream of six3 like in other animals but unlike insects. Hence, the epidermal and neural pattering function of foxQ2 is ancestral for arthropods while the upstream role of foxQ2 may have evolved in the lineage leading to the insects.

Experimental duplication of bilaterian body axes in spider embryos: Holm's organizer and self-regulation of embryonic fields

Bilaterally symmetric body plans of vertebrates and arthropods are defined by a single set of two orthogonal axes, the anterior-posterior (or head-tail) and dorsal-ventral axes. In vertebrates, and especially amphibians, complete or partial doubling of the bilaterian body axes can be induced by two different types of embryological manipulations: transplantation of an organizer region or bi-sectioning of an embryo. Such axis doubling relies on the ability of embryonic fields to flexibly respond to the situation and self-regulate toward forming a whole body. This phenomenon has facilitated experimental efforts to investigate the mechanisms of vertebrate body axes formation. However, few studies have addressed the self-regulatory capabilities of embryonic fields associated with body axes formation in non-vertebrate bilaterians. The pioneer spider embryologist Åke Holm reported twinning of spider embryos induced by both types of embryological manipulations in 1952; yet, his experiments have not been replicated by other investigators, and access to spider or non-vertebrate twins has been limited. In this review, we provide a historical background on twinning experiments in spiders, and an overview of current twinning approaches in familiar spider species and related molecular studies. Moreover, we discuss the benefits of the spider model system for a deeper understanding of the ancestral mechanisms of body axes formation in arthropods, as well as in bilaterians.

Transcriptional changes in Toxoplasma gondii in response to treatment with monensin

Background

Infection with the apicomplexan protozoan parasite T. gondii can cause severe and potentially fatal cerebral and ocular disease, especially in immunocompromised individuals. The anticoccidial ionophore drug monensin has been shown to have anti-Toxoplasma gondii properties. However, the comprehensive molecular mechanisms that underlie the effect of monensin on T. gondii are still largely unknown. We hypothesized that analysis of T. gondii transcriptional changes induced by monensin treatment can reveal new aspects of the mechanism of action of monensin against T. gondii.

Methods

Porcine kidney (PK)-15 cells were infected with tachyzoites of T. gondii RH strain. Three hours post-infection, PK-15 cells were treated with 0.1 μM monensin, while control cells were treated with medium only. PK-15 cells containing intracellular tachyzoites were harvested at 6 and 24 h post-treatment, and the transcriptomic profiles of T. gondii-infected PK-15 cells were examined using high-throughput RNA sequencing (RNA-seq). Quantitative real-time PCR was used to verify the expression of 15 differentially expressed genes (DEGs) identified by RNA-seq analysis.

Results

A total of 4868 downregulated genes and three upregulated genes were identified in monensin-treated T. gondii, indicating that most of T. gondii genes were suppressed by monensin. Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway enrichment analysis of T. gondii DEGs showed that T. gondii metabolic and cellular pathways were significantly downregulated. Spliceosome, ribosome, and protein processing in endoplasmic reticulum were the top three most significantly enriched pathways out of the 30 highly enriched pathways detected in T. gondii. This result suggests that monensin, via down-regulation of protein biosynthesis in T. gondii, can limit the parasite growth and proliferation.

Conclusions

Our findings provide a comprehensive insight into T. gondii genes and pathways with altered expression following monensin treatment. These data can be further explored to achieve better understanding of the specific mechanism of action of monensin against T. gondii.

The forkhead box containing transcription factor FoxB is a potential component of dorsal-ventral body axis formation in the spider Parasteatoda tepidariorum

In the spider, determination of the dorsal-ventral body (DV) axis depends on the interplay of the dorsal morphogen encoding gene decapentaplegic (Dpp) and its antagonist, short gastrulation (sog), a gene that is involved in the correct establishment of ventral tissues. Recent work demonstrated that the forkhead domain encoding gene FoxB is involved in dorsal-ventral axis formation in spider limbs. Here, Dpp likely acts as a dorsal morphogen, and FoxB is likely in control of ventral tissues as RNAi-mediated knockdown of FoxB causes dorsalization of the limbs. In this study, we present phenotypes of FoxB knockdown that demonstrate a function in the establishment of the DV body axis. Knockdown of FoxB function leads to embryos with partially duplicated median germ bands (Duplicitas media) that are possibly the result of ectopic activation of Dpp signalling. Another class of phenotypes is characterized by unnaturally slim (dorsal-ventrally compressed) germ bands in which ventral tissue is either not formed, or is specified incorrectly, likely a result of Dpp over-activity. These results suggest that FoxB functions as an antagonist of Dpp signalling during body axis patterning, similarly as it is the case in limb development. FoxB thus represents a general player in the establishment of dorsal-ventral structures during spider ontogeny.

Major Revisions in Arthropod Phylogeny Through Improved Supermatrix, With Support for Two Possible Waves of Land Invasion by Chelicerates

Deep phylogeny involving arthropod lineages is difficult to recover because the erosion of phylogenetic signals over time leads to unreliable multiple sequence alignment (MSA) and subsequent phylogenetic reconstruction. One way to alleviate the problem is to assemble a large number of gene sequences to compensate for the weakness in each individual gene. Such an approach has led to many robustly supported but contradictory phylogenies. A close examination shows that the supermatrix approach often suffers from two shortcomings. The first is that MSA is rarely checked for reliability and, as will be illustrated, can be poor. The second is that, to alleviate the problem of homoplasy at the third codon position of protein-coding genes due to convergent evolution of nucleotide frequencies, phylogeneticists may remove or degenerate the third codon position but may do it improperly and introduce new biases. We performed extensive reanalysis of one of such "big data" sets to highlight these two problems, and demonstrated the power and benefits of correcting or alleviating these problems. Our results support a new group with Xiphosura and Arachnopulmonata (Tetrapulmonata + Scorpiones) as sister taxa. This favors a new hypothesis in which the ancestor of Xiphosura and the extinct Eurypterida (sea scorpions, of which many later forms lived in brackish or freshwater) returned to the sea after the initial chelicerate invasion of land. Our phylogeny is supported even with the original data but processed with a new "principled" codon degeneration. We also show that removing the 1673 codon sites with both AGN and UCN codons (encoding serine) in our alignment can partially reconcile discrepancies between nucleotide-based and AA-based tree, partly because two sequences, one with AGN and the other with UCN, would be identical at the amino acid level but quite different at the nucleotide level.

Developmental gene expression as a phylogenetic data class: support for the monophyly of Arachnopulmonata

Despite application of genome-scale datasets, the phylogenetic placement of scorpions within arachnids remains contentious between two different phylogenetic data classes. Paleontologists continue to recover scorpions in a basally branching position, partly owing to their morphological similarity to extinct marine orders like Eurypterida (sea scorpions). Phylogenomic datasets consistently recover scorpions in a derived position, as the sister group of Tetrapulmonata (a clade of arachnids that includes spiders). To adjudicate between these hypotheses using a rare genomic change (RGC), we leveraged the recent discovery of ancient paralogy in spiders and scorpions to assess phylogenetic placement. We identified homologs of four transcription factors required for appendage patterning (dachshund, homothorax, extradenticle, and optomotor blind) in arthropods that are known to be duplicated in spiders. Using genomic resources for a spider, a scorpion, and a harvestman, we conducted gene tree analyses and assayed expression patterns of scorpion gene duplicates. Here we show that scorpions, like spiders, retain two copies of all four transcription factors, whereas arachnid orders like mites and harvestmen bear a single copy. A survey of embryonic expression patterns of the scorpion paralogs closely matches those of their spider counterparts, with one paralog consistently retaining the putatively ancestral pattern found in the harvestman, as well as the mite, and/or other outgroups. These data comprise a rare genomic change in chelicerate phylogeny supporting the inference of a distal placement of scorpions. Beyond demonstrating the diagnostic power of developmental genetic data as a phylogenetic data class, a derived placement of scorpions within the arachnids, together with an array of stem-group Paleozoic scorpions that occupied marine habitats, effectively rules out a scenario of a single colonization of terrestrial habitat within Chelicerata, even in tree topologies contrived to recover the monophyly of Arachnida.

Duplication and diversification of trehalase confers evolutionary advantages on lepidopteran insects

Gene duplication provides a major source of new genes for evolutionary novelty and ecological adaptation. However, the maintenance of duplicated genes and their relevance to adaptive evolution has long been debated. Insect trehalase (Treh) plays key roles in energy metabolism, growth, and stress recovery. Here, we show that the duplication of Treh in Lepidoptera (butterflies and moths) is linked with their adaptation to various environmental stresses. Generally, two Treh genes are present in insects: Treh1 and Treh2. We report three distinct forms of Treh in lepidopteran insects, where Treh1 was duplicated into two gene clusters (Treh1a and Treh1b). These gene clusters differ in gene expression patterns, enzymatic properties, and subcellular localizations, suggesting that the enzymes probably underwent sub- and/or neofunctionalization in the lepidopteran insects. Interestingly, selective pressure analysis provided significant evidence of positive selection on duplicate Treh1b gene in lepidopteran insect lineages. Most positively selected sites were located in the alpha-helical region, and several sites were close to the trehalose binding and catalytic sites. Subcellular adaptation of duplicate Treh1b driven by positive selection appears to have occurred as a result of selected changes in specific sequences, allowing for rapid reprogramming of duplicated Treh during evolution. Our results suggest that gene duplication of Treh and subsequent functional diversification could increase the survival rate of lepidopteran insects through various regulations of intracellular trehalose levels, facilitating their adaptation to diverse habitats. This study provides evidence regarding the mechanism by which gene family expansion can contribute to species adaptation through gene duplication and subsequent functional diversification.

The Dual α-Amidation System in Scorpion Venom Glands

Many peptides in scorpion venoms are amidated at their C-termini. This post-translational modification is paramount for the correct biological function of ion channel toxins and antimicrobial peptides, among others. The discovery of canonical amidation sequences in transcriptome-derived scorpion proproteins suggests that a conserved enzymatic α-amidation system must be responsible for this modification of scorpion peptides. A transcriptomic approach was employed to identify sequences putatively encoding enzymes of the α-amidation pathway. A dual enzymatic α-amidation system was found, consisting of the membrane-anchored, bifunctional, peptidylglycine α-amidating monooxygenase (PAM) and its paralogs, soluble monofunctional peptidylglycine α-hydroxylating monooxygenase (PHMm) and peptidyl-α-hydroxyglycine α-amidating lyase (PALm). Independent genes encode these three enzymes. Amino acid residues responsible for ion coordination and enzymatic activity are conserved in these sequences, suggesting that the enzymes are functional. Potential endoproteolytic recognition sites for proprotein convertases in the PAM sequence indicate that PAM-derived soluble isoforms may also be expressed. Sequences potentially encoding proprotein convertases (PC1 and PC2), carboxypeptidase E (CPE), and other enzymes of the α-amidation pathway, were also found, confirming the presence of this pathway in scorpions.

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

Chelicerates are a diverse group of arthropods, represented by such forms as predatory spiders and scorpions, parasitic ticks, humic detritivores, and marine sea spiders (pycnogonids) and horseshoe crabs. Conflicting phylogenetic relationships have been proposed for chelicerates based on both morphological and molecular data, the latter usually not recovering arachnids as a clade and instead finding horseshoe crabs nested inside terrestrial Arachnida. Here, using genomic-scale datasets and analyses optimised for countering systematic error, we find strong support for monophyletic Acari (ticks and mites), which when considered as a single group represent the most biodiverse chelicerate lineage. In addition, our analysis recovers marine forms (sea spiders and horseshoe crabs) as the successive sister groups of a monophyletic lineage of terrestrial arachnids, suggesting a single colonisation of land within Chelicerata and the absence of wholly secondarily marine arachnid orders.

A Critical Appraisal of the Placement of Xiphosura (Chelicerata) with Account of Known Sources of Phylogenetic Error

Horseshoe crabs (Xiphosura) are traditionally regarded as sister group to the clade of terrestrial chelicerates (Arachnida). This hypothesis has been challenged by recent phylogenomic analyses, but the non-monophyly of Arachnida has consistently been disregarded as artifactual. We re-evaluated the placement of Xiphosura among chelicerates using the most complete phylogenetic data set to date, expanding outgroup sampling, and including data from whole genome sequencing projects. In spite of uncertainty in the placement of some arachnid clades, all analyses show Xiphosura consistently nested within Arachnida as the sister group to Ricinulei (hooded tick spiders). It is apparent that the radiation of arachnids is an old one and occurred over a brief period of time, resulting in several consecutive short internodes, and thus is a potential case for the confounding effects of incomplete lineage sorting (ILS). We simulated coalescent gene trees to explore the effects of increasing levels of ILS on the placement of horseshoe crabs. In addition, common sources of systematic error were evaluated, as well as the effects of fast-evolving partitions and the dynamics of problematic long branch orders. Our results indicated that the placement of horseshoe crabs cannot be explained by missing data, compositional biases, saturation, or ILS. Interrogation of the phylogenetic signal showed that the majority of loci favor the derived placement of Xiphosura over a monophyletic Arachnida. Our analyses support the inference that horseshoe crabs represent a group of aquatic arachnids, comparable to aquatic mites, breaking a long-standing paradigm in chelicerate evolution and altering previous interpretations of the ancestral transition to the terrestrial habitat. Future studies testing chelicerate relationships should approach the task with a sampling strategy where the monophyly of Arachnida is not held as the premise.