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.

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.

ACME dissociation: a versatile cell fixation-dissociation method for single-cell transcriptomics

Single-cell sequencing technologies are revolutionizing biology, but they are limited by the need to dissociate live samples. Here, we present ACME (ACetic-MEthanol), a dissociation approach for single-cell transcriptomics that simultaneously fixes cells. ACME-dissociated cells have high RNA integrity, can be cryopreserved multiple times, and are sortable and permeable. As a proof of principle, we provide single-cell transcriptomic data of different species, using both droplet-based and combinatorial barcoding single-cell methods. ACME uses affordable reagents, can be done in most laboratories and even in the field, and thus will accelerate our knowledge of cell types across the tree of life.

Duplication and expression of Sox genes in spiders

BACKGROUND

The Sox family of transcription factors is an important part of the genetic 'toolbox' of all metazoans examined to date and is known to play important developmental roles in vertebrates and insects. However, outside the commonly studied Drosophila model little is known about the repertoire of Sox family transcription factors in other arthropod species. Here we characterise the Sox family in two chelicerate species, the spiders Parasteatoda tepidariorum and Stegodyphus mimosarum, which have experienced a whole genome duplication (WGD) in their evolutionary history.

RESULTS

We find that virtually all of the duplicate Sox genes have been retained in these spiders after the WGD. Analysis of the expression of Sox genes in P. tepidariorum embryos suggests that it is likely that some of these genes have neofunctionalised after duplication. Our expression analysis also strengthens the view that an orthologue of vertebrate Group B1 genes, SoxNeuro, is implicated in the earliest events of CNS specification in both vertebrates and invertebrates. In addition, a gene in the Dichaete/Sox21b class is dynamically expressed in the spider segment addition zone, suggestive of an ancient regulatory mechanism controlling arthropod segmentation as recently suggested for flies and beetles. Together with the recent analysis of Sox gene expression in the embryos of other arthropods, our findings support the idea of conserved functions for some of these genes, including a potential role for SoxC and SoxD genes in CNS development and SoxF in limb development.

CONCLUSIONS

Our study provides a new chelicerate perspective to understanding the evolution and function of Sox genes and how the retention of duplicates of such important tool-box genes after WGD has contributed to different aspects of spider embryogenesis. Future characterisation of the function of these genes in spiders will help us to better understand the evolution of the regulation of important developmental processes in arthropods and other metazoans including neurogenesis and segmentation.

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

BACKGROUND

The duplication of genes can occur through various mechanisms and is thought to make a major contribution to the evolutionary diversification of organisms. There is increasing evidence for a large-scale duplication of genes in some chelicerate lineages including two rounds of whole genome duplication (WGD) in horseshoe crabs. To investigate this further, we sequenced and analyzed the genome of the common house spider Parasteatoda tepidariorum.

RESULTS

We found pervasive duplication of both coding and non-coding genes in this spider, including two clusters of Hox genes. Analysis of synteny conservation across the P. tepidariorum genome suggests that there has been an ancient WGD in spiders. Comparison with the genomes of other chelicerates, including that of the newly sequenced bark scorpion Centruroides sculpturatus, suggests that this event occurred in the common ancestor of spiders and scorpions, and is probably independent of the WGDs in horseshoe crabs. Furthermore, characterization of the sequence and expression of the Hox paralogs in P. tepidariorum suggests that many have been subject to neo-functionalization and/or sub-functionalization since their duplication.

CONCLUSIONS

Our results reveal that spiders and scorpions are likely the descendants of a polyploid ancestor that lived more than 450 MYA. Given the extensive morphological diversity and ecological adaptations found among these animals, rivaling those of vertebrates, our study of the ancient WGD event in Arachnopulmonata provides a new comparative platform to explore common and divergent evolutionary outcomes of polyploidization events across eukaryotes.

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

Background

The duplication of genes can occur through various mechanisms and is thought to make a major contribution to the evolutionary diversification of organisms. There is increasing evidence for a large-scale duplication of genes in some chelicerate lineages including two rounds of whole genome duplication (WGD) in horseshoe crabs. To investigate this further, we sequenced and analyzed the genome of the common house spider Parasteatoda tepidariorum.

Results

We found pervasive duplication of both coding and non-coding genes in this spider, including two clusters of Hox genes. Analysis of synteny conservation across the P. tepidariorum genome suggests that there has been an ancient WGD in spiders. Comparison with the genomes of other chelicerates, including that of the newly sequenced bark scorpion Centruroides sculpturatus, suggests that this event occurred in the common ancestor of spiders and scorpions, and is probably independent of the WGDs in horseshoe crabs. Furthermore, characterization of the sequence and expression of the Hox paralogs in P. tepidariorum suggests that many have been subject to neo-functionalization and/or sub-functionalization since their duplication.

Conclusions

Our results reveal that spiders and scorpions are likely the descendants of a polyploid ancestor that lived more than 450 MYA. Given the extensive morphological diversity and ecological adaptations found among these animals, rivaling those of vertebrates, our study of the ancient WGD event in Arachnopulmonata provides a new comparative platform to explore common and divergent evolutionary outcomes of polyploidization events across eukaryotes.

Electronic supplementary material

The online version of this article (doi:10.1186/s12915-017-0399-x) contains supplementary material, which is available to authorized users.

The Wnt and Delta-Notch signalling pathways interact to direct pair-rule gene expression via caudal during segment addition in the spider Parasteatoda tepidariorum

In short germ arthropods, posterior segments are added sequentially from a growth zone or segment addition zone (SAZ) during embryogenesis. Studies in spiders such as the common house spider, Parasteatoda tepidariorum, have provided insights into the gene regulatory network (GRN) that underlies the development of the SAZ, and revealed the involvement of two important signalling pathways. It was shown that Wnt8 maintains a pool of undifferentiated cells in the SAZ, but this ligand is also required for dynamic Delta (Dl) expression associated with the formation of new segments. However, it remains unclear how these pathways interact during SAZ formation and subsequently regulate segment addition. Here we show that Delta-Notch signalling is required for Wnt8 expression in posterior SAZ cells, but represses the expression of this Wnt gene in anterior SAZ cells. We also found that these two signalling pathways are required for the expression of the spider orthologues of the segmentation genes even-skipped (eve) and runt-1 (run-1), at least in part via the transcription factor encoded by caudal (cad). Moreover, it appears that dynamic expression of eve in this spider does not require a feedback loop with run-1, as is found in the pair-rule circuit of the beetle Tribolium. Taken together, our results suggest that the development of posterior segments in Parasteatoda is directed by dynamic interactions between Wnt8 and Delta-Notch signalling that are read out by cad, which is necessary but not sufficient to regulate the expression of the pair-rule genes eve and run-1. Our study therefore provides new insights towards better understanding the evolution and developmental regulation of segmentation in other arthropods including insects.

Analysis of the Wnt gene repertoire in an onychophoran provides new insights into the evolution of segmentation

BACKGROUND

The Onychophora are a probable sister group to Arthropoda, one of the most intensively studied animal phyla from a developmental perspective. Pioneering work on the fruit fly Drosophila melanogaster and subsequent investigation of other arthropods has revealed important roles for Wnt genes during many developmental processes in these animals.

RESULTS

We screened the embryonic transcriptome of the onychophoran Euperipatoides kanangrensis and found that at least 11 Wnt genes are expressed during embryogenesis. These genes represent 11 of the 13 known subfamilies of Wnt genes.

CONCLUSIONS

Many onychophoran Wnt genes are expressed in segment polarity gene-like patterns, suggesting a general role for these ligands during segment regionalization, as has been described in arthropods. During early stages of development, Wnt2, Wnt4, and Wnt5 are expressed in broad multiple segment-wide domains that are reminiscent of arthropod gap and Hox gene expression patterns, which suggests an early instructive role for Wnt genes during E. kanangrensis segmentation.

Evolutionary conservation of the eumetazoan gene regulatory landscape

Despite considerable differences in morphology and complexity of body plans among animals, a great part of the gene set is shared among Bilateria and their basally branching sister group, the Cnidaria. This suggests that the common ancestor of eumetazoans already had a highly complex gene repertoire. At present it is therefore unclear how morphological diversification is encoded in the genome. Here we address the possibility that differences in gene regulation could contribute to the large morphological divergence between cnidarians and bilaterians. To this end, we generated the first genome-wide map of gene regulatory elements in a nonbilaterian animal, the sea anemone Nematostella vectensis. Using chromatin immunoprecipitation followed by deep sequencing of five chromatin modifications and a transcriptional cofactor, we identified over 5000 enhancers in the Nematostella genome and could validate 75% of the tested enhancers in vivo. We found that in Nematostella, but not in yeast, enhancers are characterized by the same combination of histone modifications as in bilaterians, and these enhancers preferentially target developmental regulatory genes. Surprisingly, the distribution and abundance of gene regulatory elements relative to these genes are shared between Nematostella and bilaterian model organisms. Our results suggest that complex gene regulation originated at least 600 million yr ago, predating the common ancestor of eumetazoans.

Evolution of eye morphology and rhodopsin expression in the Drosophila melanogaster species subgroup

A striking diversity of compound eye size and shape has evolved among insects. The number of ommatidia and their size are major determinants of the visual sensitivity and acuity of the compound eye. Each ommatidium is composed of eight photoreceptor cells that facilitate the discrimination of different colours via the expression of various light sensitive Rhodopsin proteins. It follows that variation in eye size, shape, and opsin composition is likely to directly influence vision. We analyzed variation in these three traits in D. melanogaster, D. simulans and D. mauritiana. We show that D. mauritiana generally has larger eyes than its sibling species, which is due to a combination of larger ommatidia and more ommatidia. In addition, intra- and inter-specific differences in eye size among D. simulans and D. melanogaster strains are mainly caused by variation in ommatidia number. By applying a geometric morphometrics approach to assess whether the formation of larger eyes influences other parts of the head capsule, we found that an increase in eye size is associated with a reduction in the adjacent face cuticle. Our shape analysis also demonstrates that D. mauritiana eyes are specifically enlarged in the dorsal region. Intriguingly, this dorsal enlargement is associated with enhanced expression of rhodopsin 3 in D. mauritiana. In summary, our data suggests that the morphology and functional properties of the compound eyes vary considerably within and among these closely related Drosophila species and may be part of coordinated morphological changes affecting the head capsule.