- Research
- Open access
- Published:
Single-cell sequencing suggests a conserved function of Hedgehog-signalling in spider eye development
EvoDevo volume 15, Article number: 11 (2024)
Abstract
Background
Spiders evolved different types of eyes, a pair of primary eyes that are usually forward pointing, and three pairs of secondary eyes that are typically situated more posterior and lateral on the spider’s head. The best understanding of arthropod eye development comes from the vinegar fly Drosophila melanogaster, the main arthropod model organism, that also evolved different types of eyes, the larval eyes and the ocelli and compound eyes of the imago. The gene regulatory networks that underlie eye development in this species are well investigated revealing a conserved core network, but also show several differences between the different types of eyes. Recent candidate gene approaches identified a number of conserved genes in arthropod eye development, but also revealed crucial differences including the apparent lack of some key factors in some groups of arthropods, including spiders.
Results
Here, we re-analysed our published scRNA sequencing data and found potential key regulators of spider eye development that were previously overlooked. Unlike earlier research on this topic, our new data suggest that Hedgehog (Hh)-signalling is involved in eye development in the spider Parasteatoda tepidariorum. By investigating embryonic gene expression in representatives of all main groups of spiders, we demonstrate that this involvement is conserved in spiders. Additionally, we identified genes that are expressed in the developing eyes of spiders, but that have not been studied in this context before.
Conclusion
Our data show that single-cell sequencing represents a powerful method to gain deeper insight into gene regulatory networks that underlie the development of lineage-specific organs such as the derived set of eyes in spiders. Overall, we gained deeper insight into spider eye development, as well as the evolution of arthropod visual system formation.
Background
The ability to process information provided by light is pervasive in animals as it represents an evolutionary advantage [52]. Starting with simple light-sensitive cells or organs, complex eyes evolved already more than half a billion years ago in the Cambrian (e.g. [3, 11, 24, 53, 87]) and may have contributed to the “Cambrian explosion” (e.g. [74]). Among arthropods, insects and spiders (Araneae) (Fig. 1A) have evolved the most complex eyes (e.g. [51], Morehouse 2020). Most spiders possess eight eyes, except for some groups that possess only six [51] or less, often due to a loss of some of the secondary eyes (e.g. [33]). The four pairs of eyes of (most) spiders are named after their position on the spider’s head, i.e. the anterior median (AM) eyes that represent the principal/primary eyes of spiders. In contrast, the anterior lateral (AL), the posterior lateral (PL), and the posterior median (PM) eyes represent the secondary eyes of spiders (Fig. 1B). There is accumulating evidence that the principal eyes of spiders are homologous with the ocelli of insects, and that the secondary eyes of spiders are homologous with the compound eyes of insects (reviewed in [63], Morehouse 2020). The three secondary pairs of eyes indeed likely evolved from compound eyes present in their aquatic ancestors albeit reducing the number of facets from about a couple of thousands to only three [63, 67].
Principal and secondary eyes of spiders both originate from the non-neurogenic ectoderm of the head lobes, but from two different regions: the principal eyes originate from the anterior edge of the head lobes, while the secondary eyes form at the lateral edge of the posterior head lobe close to the base of the chelicerae (e.g. [6, 77, 79]) (Fig. 1B). Later during development, the non-neurogenic ectoderm overgrows the neurogenic ectoderm translocating especially the principal eyes into a more anterior-facing position (Fig. 1B, C). The final position of the eyes, however, can vary in different groups of spiders (reviewed in [7, 63, 64, 86]. One of the main differences between principal and secondary eyes of spiders is that the retinas of the principal eyes are everted, while the retinas of the secondary eyes are inverted, suggesting different evolutionary origin of these two types of eyes [30]. The gene regulatory network underlying arthropod eye development appears to be widely conserved (e.g. [89]) and is thus far best studied in the vinegar fly Drosophila melanogaster (reviewed in e.g. [12, 18, 50, 92]). Interestingly, although a core gene regulatory network (GRN) appears to contribute to the development of all types of light-sensitive organs (eyes) in the fly, i.e. ocelli, compound eyes and Bolwig’s organs, there are still differences in gene regulation and interaction in the different types of eyes (reviewed in e.g. [18]). In a nutshell, Pax6 genes determine the optic fields in which the compound eyes and the ocelli develop, and they activate the retinal determination genes (RDGs) such as eyes absent (eya), and sine oculis (so/six1) ([18, 90]. Wnt-signalling inhibits eye development and Wnt gene expression thus restricts the expansion of the optic field (e.g. [54]). Downstream of the RDGs, and additionally activated by Hedgehog (Hh)-signalling, acts the proneural gene atonal (ato) in the developing eyes of Drosophila [46] that finally activates the expression of the photoreceptor specific gene glass (gl) [65].
In the common house spider Parasteatoda tepidariorum, Pax6 appears to be replaced by another Pax gene Pax2, at least in the secondary eyes [34], and Wnt-signalling likely represses the development of eyes [6]. Like in Drosophila, eya and so/six1 are expressed in all spider eyes ([20, 77, 79]), and also the RDG downstream factor ato is expressed in all eyes of Parasteatoda [6] (Fig. 1D). Data on gl expression in spiders are lacking, as well as data on Hh-signalling that would suggest a conserved function during the activation of retinal target genes. Instead, in Parasteatoda, gene expression studies suggest that not only Wnt-signalling but also Hh-signalling appear to restrict the eye fields and consequently that the eyes develop in regions of the head where neither Wnt- nor Hh-signalling occurs [6] (Fig. 1D).
We previously performed single-cell sequencing (SCS) in late-stage embryos of Parasteatoda revealing several cell clusters with a neurogenic fingerprint, including clusters representing the developing peripheral and sensory nervous system [60]. In the current study, we provide an improved analysis of this previous data using newer versions of analytical software and more restricted criteria for the definition of marker genes. The new analysis revealed a previously undetected cell cluster with an eye-specific fingerprint represented by numerous known Parasteatoda eye-marker genes (cf. e.g. [7, 79]). Interestingly, we found that the most specific genetic marker of this cluster is represented by a previously unstudied paralog of hh. In this paper we investigate the top markers of this cluster, including the previously unstudied hh paralog, and show that they are expressed in the developing eyes of Parasteatoda. For hh genes, we expanded our study to representatives of all main groups of spiders (Fig. 1A) demonstrating that Hh-signalling is involved in eye development in spiders as a whole (taking mRNA expression of the sole ligand of this pathway, hh (or its paralog hh2), as evidence for pathway activity [47]. Beyond that important finding, our data also identify additional genes that are expressed during spider eye development but that have not been studied in this context yet. Our data thus contribute to a deeper insight into eye development in spiders.
Methods
Research animals
The embryos of Parasteatoda come from our own culture in Uppsala. Embryos of Acanthoscurria and Ischnothele come from the culture of Matthias Pechmann (Cologne University). Fertilized females of Pholcus were collected in the basement of a house, 47559 Kranenburg, Germany, and now serve as the founders of a laboratory colony in Uppsala. Egg sacks were collected from females of Pardosa at the KFUM Uppsala Survival trainings facility, 75260 Uppsala, Sweden, during July and August 2021. Females of Pardosa were released back into the wild after carefully collecting their egg sacks. Embryos of the harvestman Phalangium opilio were obtained and treated as described in Janssen et al. [41].
Single-cell sequencing
Embryonic tissue dissociation, cell capture, cDNA library preparation, and single-cell RNA sequencing were performed at the Department of Developmental Biology and Gene Core facilities of the European Molecular Biology Laboratory (EMBL), Heidelberg, Germany. These processes as well as quality control of sequencing reads, and reference genome pre-processing are described in Medina-Jiménez et al. [60]. The raw sequencing reads used in the present study are the same as presented in Medina-Jiménez et al. [60]. We here present a de novo gene expression matrix generated after mapping the sequencing reads to the reference genome Ptep_3.0 (same version of the genome as in Medina-Jiménez et al. [60] using CellRanger version 7.0.1. instead of the previously used CellRanger version 6.0.2 (Supplementary File 1). Unlike the previous study, we did not apply any UMI number restrictions in the de novo analysis. The resulting gene expression matrix was loaded into Seurat ver. 4.3.0, an R package for analysis of single-cell RNA sequencing [25, 78]. Only cells with at least 250 expressed genes (min.features = 250), with no more than 5% of these genes being mitochondrial, and genes that are expressed in at least three cells (min.cells = 3) were kept for downstream analysis. Filtering of low quality, and non-informative cell barcodes was done by a back-and-forth strategy, in which cells were processed up to the clustering and marker selection steps, then clusters with markers that did not pass the pre-established threshold were removed, and the cells were processed again. Data normalization, variance stabilization, linear and non-linear dimensional reduction and clustering was done as in Medina-Jiménez et al. [60]. The Seurat function ‘FindAllMarkers’ was used to identify the differentially expressed genes (DEGs) to be used for cluster identification (Supplementary File 2). Parameters included selecting only the DEGs that are upregulated (only.pos = TRUE), and those that are expressed in at least 10% of the cells in a cluster (min.pct = 0.1). In addition, a return threshold of 1e−10 was established to avoid selecting markers that have a high adjusted p-value. The final subset is available in *.rds format (Supplementary File 3).
Gene orthology and gene ontology (GO) analysis
In order to identify Drosophila melanogaster orthologs of the cluster markers, we conducted reciprocal best hits (RBH) using the software MMSeqs2 [85]. For the gene ontology analysis of these orthologs, we applied the bioinformatics package clusterProfiler [97] using the function ‘enrichGO’ library and the Drosophila melanogaster database (or.Dm.eg.db) and focusing only on the analysis of GO: Biological Processes.
Phylogenetic analyses
Hedgehog (hh) genes were identified in the genome of the true spider Parasteatoda tepidariorum [80] and the transcriptomes of the true spiders Cupiennius salei [77], Pardosa amentata [26] and Pholcus phalangioides [39], the mygalomorph spiders Acanthoscurria geniculata [73] and Ischnothele caudata (TRINITY assembled transcriptome, embryonic stages 5–14 and the postembryonic stage), and the harvestman Phalangium opilio [81] performing reciprocal BLAST searches (tBLASTn) using the Hedgehog protein sequence from the vinegar fly Drosophila melanogaster as query. A phylogenetic analysis was performed using MrBayes [31] in the same way as previously described in Panara et al. [70]. Amino acid sequences of the complete coding regions of hh genes (as far as available) were aligned using T-Coffee with default parameters in MacVector v12.6.0, followed by manual editing (Supplementary Files 4–6). All phylogenetic analyses conducted in this paper were performed using a fixed WAG amino acid substitution model with gamma-distributed rate variation across sites (with four rate categories), unconstrained exponential prior probability distribution on branch lengths, and exponential prior for the gamma shape parameters for among-site rate variation was applied. Gene topology of the Hh-tree was analysed applying 1 million cycles for the Metropolis-Coupled Markov Chain Monte Carlo (MCMCMC) analysis (four chains; chain-heating temperature of 0.2). Markov chains were sampled every 200 cycles and default settings of 25% of samples were applied as burn-in. Clade support was calculated with posterior probabilities in MrBayes. The resulting tree was midpoint rooted.
We also performed a phylogenetic analysis with the identified CD36-family gene of Parasteatoda that we named CD36.1, additional closely related CD36-family genes of Parasteatoda, CD36-family genes from the zebrafish Danio rerio and the house mouse Mus musculus [98], and all previously identified Drosophila melanogaster CD36-family genes [66], using the same parameters as applied for the Hh-analysis (see above) but with 0.5 million cycles for the MCMCMC analysis (Supplementary Files 7 and 8). Gene identifiers are shown in the tree (Supplementary File 9).
The phylogenetic analysis of spider T-box genes including optomotor-blind (omb) genes has been published in Janssen and Budd [45].
In situ hybridization and data documentation
Fragments of all investigated genes were amplified by means of PCR using gene specific primers (Supplementary File 6). In all cases, a T7 RNA-polymerase promotor site (TAATACGACTCACTATAG) was added to the reverse primer to allow a direct use of the PCR amplicon as template for subsequent anti-sense RNA-probe synthesis [13]. Spider embryos and harvestman embryos were fixed for in situ hybridization as described in Prpic et al. [76] (for spiders) and Janssen et al. [41] (for the harvestmen) at room temperature on a moderately moving tumbling table for 16 to 24 h. All in situ hybridizations were performed as described in Janssen et al. [40]. All embryos were incubated in 1:10,000 Sybr-green in phosphate buffered saline pH 7.4 with 0.1% Tween-20 (PBST) for 30 min. Unincorporated Sybr-green was removed by washing with PBST (washing 3 times for 10 min on a tumbling table). Pictures of stained embryos were taken under a MZ-FLIII Leica dissection microscope equipped with a Leica DC490 digital camera and an external UV-light source. Linear adjustments were made using the image-processing software Adobe Photoshop 2022.
Results
Improved analysis of single-cell sequencing (SCS) data reveals new cell clusters and further specifies previously discovered cell clusters
Previously, we performed SCS on complete Parasteatoda embryos representing the developmental stages 10–12 (after [61]) in which we identified several cell clusters assigned to the developing central and peripheral nervous system [60]. Based on this pioneering study, we improved our data analysis that we publish in the current paper. This led to the discovery of new cell clusters (Fig. 2; Table 1) including new marker genes (Supplementary File 2), and the improvement of the specific genetic fingerprints of previously recovered cell clusters (cf. [60]). We based our cluster-predictions on the embryonic gene expression patterns of previously published marker genes (Table 1), newly provided expression data (Supplementary File 10), and a GO analysis (Supplementary Files 11–13). The new GO analysis provides information on the putative biological function of the cells in a given cell cluster, although this information has to be handled with care given the possibility of research biases such as above average interest in certain biological functions. This can lead to an above average of available comparative information on such a topic (and associated genes). And this can possibly lead to incorrect cluster characterization. A potential example is given by cluster C25 that according to our GO analysis would likely represent cells involved in cardiac development, although available gene expression data rather suggest another location of these cells (cf. Table 1). Whole mount in situ hybridization data that provide spatial expression data within the developing spider embryos are thus required to verify (or falsify) literature and GO analysis predictions.
In the new analysis, we identified the novel clusters C2 and C10 (both putatively representing the proximal ectoderm of developing appendages), C14 and C25 (both putatively representing the ventral ectoderm of the appendages), C19 (putatively representing a subset of the developing brain), C23 (putatively, a second cluster representing the developing stomodaeal area), C24 (putatively a second brain-specific cluster), C26 (putatively representing tissue of the pedipalps and/or the pedipalpal segment), C27 (putatively a second cluster representing the peripheral nervous system), C29 (putatively representing the book lungs and the tracheae), and C33 (putatively representing the immune system) (Table 1).
With the aid of the GO analysis, we also managed to further specify previously identified and new cell clusters. For example, clusters C0, C11, C15, C16, and C28 that were previously simply described as associated with the CNS and the midline/ventral sulcus can now be more accurately defined as also be involved in axon-guidance (Table 1; Supplementary Files 11–13). Clusters C1, C17, C19, C24, C27 appear to be generally involved in visual system development, although their marker genes are not necessarily expressed in the developing eyes themselves (examples provided in Supplementary File 10, see panels O, P, S).
We identified inter alia a small cell cluster (C32) represented by 39 cells and 25 differentially expressed marker genes (Fig. 2A; Supplementary Files 2 and 6). In our previous study, this cell cluster was named Cluster-XIX and described as “sensory nervous system of the head”-specific cell cluster (CNS-9) (Table 1); this cluster was represented by 79 cells and was defined by 42 marker genes [60]. Many of the marker genes of C32 have previously been studied in the context of spider eye development including six3.2, so1, eya, dac2, and Pax2.1 (e.g. [79]), suggesting that this cluster indeed specifically represents the developing eyes, rather than anterior sensory organs in general. Our present GO-analysis further substantiated this assumption (Supplementary File 11).We provide feature-plots of all 25 marker genes of cluster C32 (Supplementary File 14).
The most specific marker gene of C32 is hedgehog 2 (hh2), a paralog of the previously studied segment-polarity gene hedgehog (hh) (hereafter referred to as hh1). Intriguingly, this gene has not been addressed by earlier studies.
Paralogs of hedgehog
We identified three copies of hh in the transcriptomes of the entelegyne spiders Cupiennius and Pardosa, and the genome of the entelegyne spider Parasteatoda, two copies in the transcriptome of the synspermiatan spider Pholcus, and two copies in the transcriptome of the mygalomorph spider Ischnothele (Figs. 1A and 3). In another Mygalomorpha, Acanthoscurria, however, we only identified one copy of hh (Fig. 3). The latter likely represents a shortcoming of the sequenced embryonic transcriptome of Acanthoscurria rather than a gene loss, given that two copies exist in Ischnothele; a loss of hh2 in Acanthoscurria can, however, not be excluded.
Our phylogenetic analysis indicates that one spider hh paralog (Fig. 3, green background, referred to as the “hh1-group”) has been duplicated once more in RTA-class spiders represented by Cupiennius and Pardosa, but not so in the other groups of spiders (Figs. 1A and 3).
The second paralog of hh (Fig. 3, blue background, referred to as the “hh2-group”) duplicated in Parasteatoda, but not in the other investigated spiders. One of these copies underwent drastic changes losing the ‘Hedge’ domain and instead acquired two N-terminal EGF domains (Supplementary File 4) (see [10] for a review on hh and hh-related genes and their protein motifs).
In representatives of other groups of arthropods, an onychophoran, a tardigrade and a priapulid, we only identified one single copy of hh. Note that these sequences do not cluster according to animal phylogeny. The hh1 and hh2 groups of spiders, however, each form distinct clusters (Fig. 3). The source of this could be a whole genome duplication as suggested by some authors (e.g. [80, 82]), but in the case of hh, this would require more extensive analyses.
Expression of spider hedgehog (hh) genes
Expression of one hh paralog (i.e. hh1) has been described in detail in the spider Parasteatoda [2, 6, 15, 48, 69, 72].
All here identified hh1 genes in the investigated spiders are expressed in virtually the same pattern as in Parasteatoda, including the SPG-like expression in the form of transverse segmental stripes, expression in the appendages, dynamic expression in the SAZ, and expression in the stomodaeum and the pre-cheliceral region (Supplementary File 15). In Pardosa, that possesses two copies of hh1 (hh1a and hh1b), hh1a is not expressed in the SAZ during the process of segment addition (Supplementary File 15, panel B). In Acanthoscurria there are four small dots of expression anterior to the labrum in the neuronal ectoderm at late developmental stages; this expression, however, is not associated with the developing eyes (Supplementary File 15, panel O). hh1 is thus not expressed in the developing eyes of Parasteatoda [6] or any other investigated spider.
In all investigated spiders, Parasteatoda, Pardosa, Pholcus and Ischnothele, expression of hh2 genes is restricted to regions of the developing head (with the exception of the duplicated hh2-related gene of Parasteatoda). In Parasteatoda, hh2 is transiently expressed in the developing labrum (Fig. 4A). At around stage 10.2 (after [61]), expression of hh2 starts in the form of two dots laterally on either side of the head lobes in the non-neurogenic ectoderm in the place where the primordia of the secondary eyes are located (Fig. 4B). Later, this expression transforms into six distinct dots that represent the six secondary eyes (Fig. 4C). We did not detect any expression of hh2 in the principal eyes of Parasteatoda.
In Pardosa, we detect a similar pattern in the secondary eyes (Fig. 4D–F), but we also detect faint expression of hh2 in the principal eyes at later developmental stages (Fig. 4F, G).
In Pholcus, we find weak dot-like expression in the lateral non-neurogenic ectoderm that is likely associated with the developing secondary eyes (Fig. 4H–J). Moreover, we detect tiny dots of expression anterior in the non-neurogenic ectoderm where the primordia of the principal eyes are located (Fig. 4H, J). Later during development, we see strong expression in the principal and the secondary eyes (detection of the latter is, however, inhibited in two of the three secondary eyes due to the relatively early development of the cuticle that stains un-specifically and covers at this stage two of the three secondary eyes at either side of the head lobe) (Fig. 4J).
In Ischnothele, like in Parasteatoda, hh2 is expressed in the labrum (Fig. 4K–M). At late developmental stages, hh2 also is expressed in the primordia of the primary eyes (Fig. 4L, M) and weakly in the developing secondary eyes (Fig. 4L, M); the latter, however, is almost below detectable level at the stages for which we have embryos.
Although we did not detect expression of hh2 in the labrum in Pardosa and Pholcus, this could simply be due to the lack of embryos representing all developmental stages, including the stage during which hh2 is expressed in the labrum.
The duplicated and modified hh2-class gene found in Parasteatoda, referred to as Parasteatoda-hh2-related (Pt-hh2-r) is expressed in the dorsal field, and hence is likely involved in either midgut development (Supplementary File 16, panels A and B) (cf. [15]), yolk-metabolism, or the innate immune system (cf. [60]).
Expression of the Phalangium opilio hedgehog (hh) gene
We identified a single copy of hh in the genome and an embryonic transcriptome of the harvestman Phalangium [21, 81], and investigated expression of this gene in embryos ranging from stage 8 to 12 (staging after [22]). Phalangium hh is expressed in a typical segment-polarity gene pattern in the form of transverse segmental stripes (Fig. 5A–C, E), including a stripe in the pre-cheliceral region. Phalangium hh is also posteriorly expressed in the appendages (Fig. 5A, B) and in the stomodaeum (Fig. 5C-E). At later developmental stages, expression is seen in the developing eyes (Fig. 5F; cf. expression marked by an arrow with the position of the developing eyes [22, 23]).
Expression of additional eye-markers identified in our SCS experiment: SoxF1, CD36.1, optomotor-blind3, sevenless and engrailed-2
Sox genes have previously been studied in Parasteatoda, and expression of SoxF1 has been reported to be restricted to late developmental stages when it is expressed in the secondary eyes and the developing spinnerets ([5], their supplementary data). In addition to the previously reported pattern of SoxF1 we found that at late stages, expression is also in the principal eyes, not only the secondary eyes (Fig. 6A–C). Interestingly, we also discovered specific expression of SoxF2 in the anlagen of the principal and secondary eyes of Parasteatoda (Fig. 6D), suggesting that both genes are likely involved in eye development, and that the early and late expression of SoxF2 and SoxF1, respectively, could represent another case of sub-functionalization.
Parasteatoda CD36.1 is first expressed in a salt-and-pepper pattern in the central nervous system (Fig. 6F, G). Within the head, there are two domains of stronger expression in the neurogenic ectoderm close to the lateral furrow (Fig. 6A). Later, this tissue contributes to the developing brain that still strongly expresses CD36.1 (Fig. 6H, I). At this point during development there are eight distinct domains of expression that clearly correlate with the six secondary and the two principal eyes (Fig. 6H, I).
The Parasteatoda optomotor-blind3 (omb3) gene has also been investigated before, but expression has only been shown and described for earlier developmental stages prior to late eye-development ([55, 56]; aug3.g3790). Like other omb genes [36] (omb = omb1) (Supplementary File 16, panels C-K), also omb3 is expressed in the pre-cheliceral region and in dorsal tissue of prosomal and opisthosomal appendages (Fig. 6J, K). At later developmental stages, omb3 is expressed in a complex pattern in the developing appendages, and indeed also specifically in the developing secondary eyes, unlike its paralogs omb1, omb2 and omb4 (Fig. 6J, K and Supplementary File 16, panels C-K).
Previously published data on sevenless (sev) were not conclusive concerning expression in the developing eyes [60]. Here we show that sev is also expressed in the primordia of the secondary and principal eyes (Fig. 6L, M).
Likewise, previously published data on engrailed-2 (en2) did not cover later developmental stages [42]. Here we show that en2 is expressed in the six secondary eyes, but not in the two principal eyes, in late-stage embryos (Fig. 6N).
Expression of Parasteatoda epidermal growth factor receptor and glass
Epidermal growth factor receptor (Egfr) and glass (gl) play important roles during eye development in the fly Drosophila (reviewed in [18]). Although these two crucial factors were not among the enriched marker genes in the C32 cluster, we found the strongest expression of gl in a subset of C32 and Egfr was moderately expressed in all cell types (Fig. 7A). Therefore, we also investigated the embryonic expression to test whether they may play conserved roles in spider eye development as well.
We find that Egfr indeed is expressed ubiquitously, but is upregulated in certain tissues, including the primordia of the primary and secondary eyes (Fig. 7A-C). Interestingly, at later developmental stages, when the indentations of the secondary eyes have formed, Egfr appears to be specifically downregulated in these areas (Fig. 7C). Such downregulation is not seen in the area where the primary eyes form. Another interesting aspect of Egfr is its expression in the form of rings in the developing appendages that suggest a function in joint formation (Fig. 7D).
Expression of gl is restricted to later developmental stages and is weakly recognizable in the secondary eyes, but strongly in the primary eyes (Fig. 7A, E, F).
Discussion
The role of hh in spider and panarthropod eye development
The hh1 gene of spiders is expressed in the form of a transverse stripe in the pre-cheliceral region which is part of the head lobes. This domain, however, does not cover the primordia of the secondary eyes (Supplementary File 15). Baudouin-Gonzalez et al. [6] suggested that hh1 (hh in their paper) could aid Wnt-signalling in limiting the optical field and thus restrict eye growth/size. In any case, the expression of spider hh1 clearly shows that this gene cannot be involved in retina development and differentiation as it is the case in Drosophila (reviewed in e.g. [18]). In Drosophila, hh is activated by So (e.g. [74]), and subsequently, expression of ato is dependent on Hh-signalling (e.g. [14, 27]) making Hh-signalling a crucial component of eye development in the fly. The apparent lack of hh1 expression in the developing eyes of spiders raised the question if this crucial function of hh could have been gained somewhere in the lineage leading to Drosophila, or if it could have been lost in the lineage leading to spiders. Both scenarios would imply substantial differences in the eye-developmental gene regulatory network in either of these two groups of arthropods.
The finding of a second paralog of hh in spiders (i.e. hh2) answers this question in a very satisfying way, showing that Hh-signalling (as inferred by the expression of its ligand encoding gene, hh2) is indeed involved in eye development in both the fly and spiders. The relatively late onset of hh2 expression in the developing eyes is in line with a dependence on so as it also is the case in Drosophila (e.g. [71]; cf. [7]). The most likely scenario in spiders is that the two copies of hh have undergone sub-functionalization after a duplication that occurred in the lineage leading to Arachnopulmonata (i.e. scorpions, whip spiders, and spiders) (e.g. [82]), and indeed similar cases of sub-functionalization in spiders have been recorded frequently (e.g. [1, 41, 55, 80]).
Despite the fact that hh has been the subject of various studies in various groups of arthropods including other chelicerates, myriapods, crustaceans and insects beyond Drosophila (e.g. [16, 32, 35, 37, 62, 84]), data on the potential role of hh in eye development are scarce. In the scorpion Euscorpius flavicaudis, for example, expression has only been reported for one paralog of hh and during stages that reflect germ band extension [84]. The investigated hh gene, however, is expressed in the semi-lunar grooves at late germ band stages, and this is where the eyes of the scorpion will develop [58, 84]). In the myriapod Glomeris marginata, the single hh ortholog is expressed in exactly the place where the eyes will form ([35], their Fig. 6G; cf. [75] for expression of Pax6 as marker of Glomeris eye-primordia). In insects other than Drosophila, hh is at least expressed anteriorly in the head lobes during germ band formation and extension [4, 62], but again conclusive data on a potential function of hh in eye development in these species are missing. As we show in this paper, expression of hh in the harvestman Phalangium is very similar to that in other arthropods, including expression in the developing eyes. We thus conclude that this expression, including expression in the eyes, is ancestral and covered by the two paralogs in spiders.
In the outgroups, tardigrades and onychophorans, hh expression has only been studied in the closely related velvet worms Euperipatoides kanangrensis [38] and Euperipatoides rowelli [17], but with a focus on body axis segmentation. A closer look at the expression of hh in Euperipatoides kanangrensis, however, revealed no expression in the developing eyes, at least not prior to and including stage 20 [38]. Due to the lack of further comparative data, especially those from tardigrades, it thus remains unclear if the role of hh in eye development originated in the lineage leading to arthropods, the lineage leading to panarthropods, or even earlier, representing a case of deep homology (cf. data on hh function in vertebrates (e.g. [68, 83]).
Primary eyes of spiders with or without hh2 expression?
According to our data, expression of hh2 in the developing eyes is possibly not fully conserved in different groups of spiders, at least with respect to expression in the primary eyes (Fig. 4). The detected pattern, however, is not reflecting spider phylogeny [49] because we see expression in the primary eyes of the entelegyne RTA-class spider Pardosa, not in the more basally branching entelegyne spider Parasteatoda, but in the most basally branching synspermiatan spider Pholcus, and the mygalomorph spider Ischnothele (Fig. 4). This could mean independent recruitment of hh2 in the development of the primary eyes in the lineages leading to Ischnothele, Pholcus and Pardosa, or more likely a secondary loss in the lineage leading to Parasteatoda or failure to detect hh2 in Parasteatoda due to heterochrony (delayed onset of hh2 expression in the primary eyes). In this context, it may be interesting to mention that we cannot detect en2, the putative activator of hh2, in the primary eyes of Parasteatoda either [42] (Fig. 6N).
A very recent study reported differences in the expression of some RDGs in the different eye types of spiders suggesting that these differences may indeed be correlated with different eye types, specification/specialization of eye types, and consequently ecological adaptation of spiders [7]. Similar differences were previously reported for Pax6 genes in spiders (cf. [77, 79]). The differences seen in hh2 expression in the principal eyes of some spiders and the lack in others may thus contribute to the proposed differentiation, specialization and ecological adaptation reported previously.
Hypothesized ancient interactions in visual sense organ development
In one of his review papers on Drosophila eye development, Friedrich [18] suggested a core genetic program that mediates visual sense organ (eye) development, because the factors involved in this network represent key players of all visual sense organs in the fly, i.e. the compound eyes, the ocelli, and the Bolwig organs (present in the larvae). Here, Pax6 is seen as the direct activator of eya and so, dpp as an activator of eya, and wg (or Wnt-signalling in general) as an inhibitor of eya and so in optic field specification. eya and so, together with hh, then activate ato, and ato together with epidermal growth factor receptor (Egfr)-signalling activates gl, the final target of this GRN [18] (Fig. 8).
In this scenario, hh is thus a conserved ancient member of the GRN that regulates eye development. From top to bottom of this GRN and in the “eyes” of a spider, firstly Pax6 is not expressed in all developing eyes [7, 77, 79]. A recent study, however, revealed that another Pax gene, Pax2 could substitute for Pax6 [19, 34]). eya and so are both expressed in the primordia of all eyes in spiders [7, 79], and this is where Wnt-signalling mostly does not take place [6, 41]. Expression of dpp in or near the developing eyes, however, could not be detected raising the question if this component of the suggested GRN may not be part of the corresponding GRN in spiders [6]. The combined expression data from several groups of spiders presented in this paper, however, suggest that Hh-signalling, as represented by the expression of hh2 contributes to the development of all eyes (Fig. 4), and as predicted by Friedrich [18], also ato is expressed in all eyes of spiders [6, 7]. The expression of Egfr appears to contradict a specific function in spider eye development because it is downregulated in the region where the secondary eyes form. Interestingly, however, it has been reported previously that Egfr mRNA is downregulated in certain tissues of Drosophila as the result of elevated Egfr-signalling in these tissues [88]. The specific downregulation of Egfr mRNA in the secondary eye primordia (Fig. 7C) could thus be the result of preceding strong Egfr-signalling (Fig. 7B). The gene expression analysis is thus in line with a conserved role of Egfr-signalling and Ato in the regulation of the terminal gene gl which also is expressed in the primordia of all spider eyes (Fig. 7E, F). The combined data on spider eye development, including our new data on hh2, initially identified in our SCS data set, Egfr, and gl thus support the scenario that Friedrich [18] suggested about a conserved ancestral eye-developmental GRN (Fig. 8). Surprisingly, this network appears to be conserved in all eyes of spiders but has been modified in the eyes of flies. The eyes of spiders thus appear to represent less-derived optical organs, at least when compared to insects, making spider eyes interesting targets for further investigation and the understanding of eye evolution in arthropods and animals in general.
SoxF, a new player in spider eye development
SoxF genes are not expressed during panarthropod eye development [40], except for the previously reported expression of SoxF1 in the primordia of the developing secondary eyes of Parasteatoda [5]. Expression and function of Parasteatoda SoxF1 thus clearly represents a case of neo-functionalization and identifies SoxF1 as a new player in the concert of spider eye development. From the reported data, it appeared tempting to speculate that SoxF1 is a unique new factor of secondary eye development [5]. We, however, also detected faint expression of SoxF1 in the principal eyes of Parasteatoda (Fig. 6C). We conclude that SoxF1 is not a specific factor of secondary eye development, but rather a general new factor of spider eye development. The expression of SoxF2 (Fig. 6D) is in line with an earlier function of this gene in the primordia of both the secondary and principal eyes of the spider (cf. expression of so1 in Parasteatoda [79]), and thus SoxF orthologs could have been the target of temporal sub-functionalization in arachnopulmonates, or at least spiders, after the recruitment as eye-developmental factors. The remaining question is when SoxF has been recruited as an eye-developmental gene. Due to the lack of comparative data from other spiders, and indeed other chelicerates, this question needs to be addressed in the future.
CD36.1, a spider CD36-family gene that is expressed in the developing eyes
The CD36 gene family comprises four orthologs in mammals, six in the nematode worm Caenorhabditis elegans, two in an onychophoran, five in a tick (Chelicerata), eight in a myriapod, six to eight in crustaceans, and between 12 and 14 in insects [9, 29, 66, 93]. These genes fulfil a plethora of functions including cytoadhesion, carotenoid transport, and chemoreception. Many of these genes have not been investigated in Drosophila, other insects, or indeed any other (pan) arthropod species rendering comparative analysis difficult. At least in Drosophila, however, comprehensive gene expression studies and functional analysis have identified some specific functions of some members (summarized in [29]). Among them are two genes, NinaD and santa-maria, that are involved in the redistribution of carotenoids to the developing eyes (NinaD; [94]) and the biosynthesis of rhodopsin in neurons and glia cells outside the retina (santa-maria; [95]). The expression pattern of santa-maria in the ventral nerve cord, the brain, and the developing eyes is at least superficially similar to that of the Parasteatoda CD36.1 gene. This could suggest a conserved function of arthropod Group-2 CD36-family genes in arthropod nervous system and eye development. Addressing this question, however, would at least require a comprehensive analysis of arthropod CD36-family genes and accompanying gene expression analyses.
Conclusions
In this paper we show that single-cell sequencing (SCS) of developing embryos can reveal specific tissue- and structure-specific new and unexpected genes, and thus can provide insight into correlated gene regulatory networks (GRNs) beyond the classical candidate gene approach. SCS will thus undoubtedly help untangling complex GRNs in model and non-model organisms, especially with respect to differences between those new models and the well-studied established models such as the fly Drosophila, the mouse Mus and the nematode worm Caenorhabditis. Once identified, the evolutionary origin of new molecular players (genes) involved in the development of morphological structures can be tested by means of comparative studies as exemplified in this study for the hh2 paralog in different groups of spiders.
Data Availability
The datasets supporting the conclusions of this article are included within the article (and its additional files).
References
Aase-Remedios ME, Janssen R, Leite DJ, Sumner-Rooney L, McGregor AP. Evolution of the spider homeobox gene repertoire by tandem and whole genome duplication. Mol Biol Evol. 2023. https://doi.org/10.1093/molbev/msad239.
Akiyama-Oda Y, Oda H. Cell migration that orients the dorsoventral axis is coordinated with anteroposterior patterning mediated by Hedgehog signaling in the early spider embryo. Development. 2010;137(8):1263–73.
Arendt D, Wittbrodt J. Reconstructing the eyes of Urbilateria. Philos Trans R Soc Lond B Biol Sci. 2001;356(1414):1545–63.
Auman T, Chipman AD. Growth zone segmentation in the milkweed bug Oncopeltus fasciatus sheds light on the evolution of insect segmentation. BMC Evol Biol. 2018;18(1):178.
Baudouin-Gonzalez L, Schoenauer A, Harper A, Blakeley G, Seiter M, Arif S, Sumner-Rooney L, Russell S, Sharma PP, McGregor AP. The evolution of sox gene repertoires and regulation of segmentation in arachnids. Mol Biol Evol. 2021;38(8):3153–69.
Baudouin-Gonzalez L, Harper A, McGregor AP, Sumner-Rooney L. Regulation of eye determination and regionalization in the spider parasteatoda tepidariorum. Cells. 2022;11(4):631.
Baudouin-Gonzalez L, Schönauer A, Harper A, Arif S, Leite DJ, Steinhoff POM, Pechmann M, Telizhenko V, Pande A, Kosiol C, McGregor AP, Sumner-Rooney L. Development and patterning of a highly versatile visual system in spiders. biorxiv. 2023. https://doi.org/10.1101/2023.12.22.572789.
Bonatto Paese CL, Leite DJ, Schönauer A, McGregor AP, Russell S. Duplication and expression of Sox genes in spiders. BMC Evol Biol. 2018;18(1):205.
Boohar RT, Vandepas LE, Traylor-Knowles N, Browne WE. Phylogenetic and protein structure analyses provide insight into the evolution and diversification of the CD36 domain “Apex” among scavenger receptor class B proteins across Eukarya. Genome Biol Evol. 2023. https://doi.org/10.1093/gbe/evad218.
Bürglin TR. The Hedgehog protein family. Genome Biol. 2008;9(11):241.
Budd GE. The origin and evolution of the euarthropod labrum. Arthropod Struct Dev. 2021;62: 101048.
Casares F, Almudi I. Fast and furious 800. the retinal determination gene network in Drosophila. In: Castelli-Gair Hombría J, Bovolenta P, editors. Organogenetic gene networks: Genetic control of organ formation. Cham: Springer International Publishing; 2016. p. 95–124.
David R, Wedlich D. PCR-based RNA probes: a quick and sensitive method to improve whole mount embryo in situ hybridizations. Biotechniques. 2001;30(4):769–72.
Domínguez M. Dual role for Hedgehog in the regulation of the proneural gene atonal during ommatidia development. Development. 1999;126(11):2345–53.
Feitosa NM, Pechmann M, Schwager EE, Tobias-Santos V, McGregor AP, Damen WGM, Nunes da Fonseca R. Molecular control of gut formation in the spider Parasteatoda tepidariorum. Genesis. 2017. https://doi.org/10.1002/dvg.23033.
Farzana L, Brown SJ. Hedgehog signaling pathway function conserved in Tribolium segmentation. Dev Genes Evol. 2008;218(3–4):181–92.
Franke FA, Mayer G. Controversies surrounding segments and parasegments in onychophora: insights from the expression patterns of four “segment polarity genes” in the peripatopsid Euperipatoides rowelli. PLoS ONE. 2014;9(12): e114383.
Friedrich M. Ancient mechanisms of visual sense organ development based on comparison of the gene networks controlling larval eye, ocellus, and compound eye specification in Drosophila. Arthropod Struct Dev. 2006;35(4):357–78.
Friedrich M. Coming into clear sight at last: Ancestral and derived events during chelicerate visual system development. BioEssays. 2022;44(12): e2200163.
Gainett G, Ballesteros JA, Kanzler CR, Zehms JT, Zern JM, Aharon S, Gavish-Regev E, Sharma PP. Systemic paralogy and function of retinal determination network homologs in arachnids. BMC Genomics. 2020;21(1):811.
Gainett G, González VL, Ballesteros JA, Setton EVW, Baker CM, Barolo Gargiulo L, Santibáñez-López CE, Coddington JA, Sharma PP. The genome of a daddy-long-legs (Opiliones) illuminates the evolution of arachnid appendages. Proc Biol Sci. 2021;288(1956):20211168.
Gainett G, Crawford AR, Klementz BC, So C, Baker CM, Setton EVW, Sharma PP. Eggs to long-legs: embryonic staging of the harvestman Phalangium opilio (Opiliones), an emerging model arachnid. Front Zool. 2022;19(1):11.
Gainett G, Klementz BC, Blaszczyk P, Setton EVW, Murayama GP, Willemart R, Gavish-Regev E, Sharma PP. Vestigial organs alter fossil placements in an ancient group of terrestrial chelicerates. Curr Biol. 2024;34(6):1258–70.
Gehring WJ. Historical perspective on the development and evolution of eyes and photoreceptors. Int J Dev Biol. 2004;48(8–9):707–17.
Hao Y, Hao S, Andersen-Nissen E, Mauck WM 3rd, Zheng S, Butler A, Lee MJ, Wilk AJ, Darby C, Zager M, Hoffman P, Stoeckius M, Papalexi E, Mimitou EP, Jain J, Srivastava A, Stuart T, Fleming LM, Yeung B, Rogers AJ, McElrath JM, Blish CA, Gottardo R, Smibert P, Satija R. Integrated analysis of multimodal single-cell data. Cell. 2021;184(13):3573-3587.e29.
Harper A, Baudouin Gonzalez L, Schönauer A, Janssen R, Seiter M, Holzem M, Arif S, McGregor AP, Sumner-Rooney L. Widespread retention of ohnologs in key developmental gene families following whole-genome duplication in arachnopulmonates. G3 (Bethesda). 2021. https://doi.org/10.1093/g3journal/jkab299.
Heberlein U, Singh CM, Luk AY, Donohoe TJ. Growth and differentiation in the Drosophila eye coordinated by hedgehog. Nature. 1995;373(6516):709–11.
Heingård M, Turetzek N, Prpic NM, Janssen R. FoxB, a new and highly conserved key factor in arthropod dorsal-ventral (DV) limb patterning. EvoDevo. 2019;8(10):28.
Herboso L, Talamillo A, Pérez C, Barrio R. Expression of the scavenger receptor Class B type I (SR-BI) family in Drosophila melanogaster. Int J Dev Biol. 2011;55(6):603–11.
Homann H. Die Augen der Araneae: Anatomie, Ontogenie und Bedeutung für die Systematik (Chelicerata, Arachnida). Z Morphol Tiere. 1971;69:201–72.
Huelsenbeck JP, Ronquist F. MRBAYES: Bayesian inference of phylogenetic trees. Bioinformatics. 2001;17(8):754–5.
Hunnekuhl VS, Akam M. Formation and subdivision of the head field in the centipede Strigamia maritima, as revealed by the expression of head gap gene orthologues and hedgehog dynamics. EvoDevo. 2017;23(8):18.
Jäger P. Revision of the genus Sinopoda Jäger, 1999 in Laos with discovery of the first eyeless huntsman spider species (Sparassidae: Heteropodinae). Zootaxa. 2012;57:37–57.
Janeschik M, Schacht MI, Platten F, Turetzek N. 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. Front Ecol Evol. 2022;10: 810077.
Janssen R, Prpic NM, Damen WG. Gene expression suggests decoupled dorsal and ventral segmentation in the millipede Glomeris marginata (Myriapoda: Diplopoda). Dev Biol. 2004;268(1):89–104.
Janssen R, Feitosa NM, Damen WG, Prpic NM. The T-box genes H15 and optomotor-blind in the spiders Cupiennius salei, Tegenaria atrica and Achaearanea tepidariorum and the dorsoventral axis of arthropod appendages. Evol Dev. 2008;10(2):143–54.
Janssen R. Segment polarity gene expression in a myriapod reveals conserved and diverged aspects of early head patterning in arthropods. Dev Genes Evol. 2012;222(5):299–309.
Janssen R, Budd GE. Deciphering the onychophoran “segmentation gene cascade”: Gene expression reveals limited involvement of pair rule gene orthologs in segmentation, but a highly conserved segment polarity gene network. Dev Biol. 2013;382(1):224–34.
Janssen R, Schönauer A, Weber M, Turetzek N, Hogvall M, Goss GE, Patel NH, McGregor AP, Hilbrant M. The evolution and expression of panarthropod frizzled genes. Front Ecol Evol. 2015;3:96.
Janssen R, Andersson E, Betnér E, Bijl S, Fowler W, Höök L, Leyhr J, Mannelqvist A, Panara V, Smith K, Tiemann S. Embryonic expression patterns and phylogenetic analysis of panarthropod sox genes: insight into nervous system development, segmentation and gonadogenesis. BMC Evol Biol. 2018;18(1):88.
Janssen R, Pechmann M, Turetzek N. A chelicerate Wnt gene expression atlas: novel insights into the complexity of arthropod Wnt-patterning. EvoDevo. 2021;12(1):12.
Janssen R, Turetzek N, Pechmann M. Lack of evidence for conserved parasegmental grooves in arthropods. Dev Genes Evol. 2022;232(1):27–37.
Janssen R, Schomburg C, Prpic NM, Budd GE. A comprehensive study of arthropod and onychophoran Fox gene expression patterns. PLoS ONE. 2022;17(7): e0270790.
Janssen R, Budd GE. Expression of netrin and its receptors uncoordinated-5 and frazzled in arthropods and onychophorans suggests conserved and diverged functions in neuronal pathfinding and synaptogenesis. Dev Dyn. 2023;252(1):172–85.
Janssen R, Budd GE. New insights into mesoderm and endoderm development, and the nature of the onychophoran blastopore. Front Zool. 2024;21(1):2.
Jarman AP, Grell EH, Ackerman L, Jan LY, Jan YN. Atonal is the proneural gene for Drosophila photoreceptors. Nature. 1994;369(6479):398–400.
Jiménez-Jiménez C, Grobe K, Guerrero I. Hedgehog on the move: glypican-regulated transport and gradient formation in drosophila. Cells. 2024;13(5):418.
Kanayama M, Akiyama-Oda Y, Nishimura O, Tarui H, Agata K, Oda H. Travelling and splitting of a wave of hedgehog expression involved in spider-head segmentation. Nat Commun. 2011;11(2):500.
Kulkarni S, Wood HM, Hormiga G. Advances in the reconstruction of the spider tree of life: a roadmap for spider systematics and comparative studies. Cladistics. 2023;39(6):479–532.
Kumar JP. The molecular circuitry governing retinal determination. Biochim Biophys Acta. 2009;1789(4):306–14.
Land MF. The morphology and optics of spider eyes. In: Barth FG, editor. Neurobiology of Arachnids. Berlin: Springer; 1985. p. 53–78.
Land MF, Nilsson D-E. Animal eyes. Oxford, UK: Oxford Univ Press; 2012.
Lee MS, Jago JB, García-Bellido DC, Edgecombe GD, Gehling JG, Paterson JR. Modern optics in exceptionally preserved eyes of early Cambrian arthropods from Australia. Nature. 2011;474(7353):631–4.
Legent K, Treisman JE. Wingless signaling in Drosophila eye development. Methods Mol Biol. 2008;469:141–61.
Leite DJ, Baudouin-Gonzalez L, Iwasaki-Yokozawa S, Lozano-Fernandez J, Turetzek N, Akiyama-Oda Y, Prpic NM, Pisani D, Oda H, Sharma PP, McGregor AP. Homeobox gene duplication and divergence in arachnids. Mol Biol Evol. 2018;35(9):2240–53.
Leite DJ, Schönauer A, Blakeley G, Harper A, Garcia-Castro H, Baudouin-Gonzalez L, Wang R, Sarkis N, Nikola AG, Koka VSP, Kenny NJ, Turetzek N, Pechmann M, Solana J, McGregor AP. An atlas of spider development at single-cell resolution provides new insights into arthropod embryogenesis. EvoDevo. 2024;15(1):5.
Linne V, Eriksson BJ, Stollewerk A. Single-minded and the evolution of the ventral midline in arthropods. Dev Biol. 2012;364(1):66–76.
Loria SF, Prendini L. Homology of the lateral eyes of scorpiones: a six-ocellus model. PLoS ONE. 2014;9(12): e112913.
Medina-Jiménez BI, Budd GE, Janssen R. Panarthropod tiptop/teashirt and spalt orthologs and their potential role as “trunk”-selector genes. EvoDevo. 2021;12(1):7.
Medina-Jiménez BI, Budd GE, Janssen R. Single-cell RNA sequencing of mid-to-late stage spider embryos: new insights into spider development. BMC Genomics. 2024;25(1):150.
Mittmann B, Wolff C. Embryonic development and staging of the cobweb spider Parasteatoda tepidariorum C. L. Koch 1841 (syn.: Achaearanea tepidariorum; Araneomorphae Theridiidae;). Dev Genes Evol. 2012. https://doi.org/10.1007/s00427-012-0401-0.
Miyawaki K, Mito T, Sarashina I, Zhang H, Shinmyo Y, Ohuchi H, Noji S. Involvement of Wingless/Armadillo signaling in the posterior sequential segmentation in the cricket, Gryllus bimaculatus (Orthoptera), as revealed by RNAi analysis. Mech Dev. 2004;121(2):119–30.
Morehouse NI, Buschbeck EK, Zurek DB, Steck M, Porter ML. Molecular evolution of spider vision: new opportunities. Familiar Players Biol Bull. 2017;233(1):21–38.
Morehouse N. Spider vision. Curr Biol. 2020;30(17):R975–80.
Moses K, Ellis MC, Rubin GM. The glass gene encodes a zinc-finger protein required by Drosophila photoreceptor cells. Nature. 1989;340(6234):531–6.
Nichols Z, Vogt RG. The SNMP/CD36 gene family in diptera, hymenoptera and coleoptera: drosophila melanogaster, D. pseudoobscura, anopheles gambiae, aedes aegypti, apis mellifera, and tribolium castaneum. Insect Biochem Mol Biol. 2008;38(4):398–415.
Nilsson D-E. Optics and evolution of the compound eye. In: Stavenga DG, Hardie R, editors. Facets of Vision. Berlin: Springer; 1989. p. 30–73.
Neumann CJ, Nuesslein-Volhard C. Patterning of the zebrafish retina by a wave of sonic hedgehog activity. Science. 2000;289(5487):2137–9.
Oda H, Nishimura O, Hirao Y, Tarui H, Agata K, Akiyama-Oda Y. Progressive activation of Delta-Notch signaling from around the blastopore is required to set up a functional caudal lobe in the spider Achaearanea tepidariorum. Development. 2007;134(12):2195–205.
Panara V, Budd GE, Janssen R. Phylogenetic analysis and embryonic expression of panarthropod Dmrt genes. Front Zool. 2019;2(16):23.
Pauli T, Seimiya M, Blanco J, Gehring WJ. Identification of functional sine oculis motifs in the autoregulatory element of its own gene, in the eyeless enhancer and in the signalling gene hedgehog. Development. 2005;132(12):2771–82.
Pechmann M, McGregor AP, Schwager EE, Feitosa NM, Damen WG. Dynamic gene expression is required for anterior regionalization in a spider. Proc Natl Acad Sci U S A. 2009;106(5):1468–72.
Pechmann M. Embryonic development and secondary axis induction in the Brazilian white knee tarantula acanthoscurria geniculata, C. L. Koch, 1841 (Araneae; Mygalomorphae; Theraphosidae). Dev Genes Evol. 2020;230(2):75–94.
Plotnick RE, Dornbos SQ, Chen J. Information landscapes and sensory ecology of the Cambrian radiation. Paleobiology. 2010;36(2):303–17.
Prpic NM. Duplicated pax6 genes in glomeris marginata (Myriapoda: Diplopoda), an arthropod with simple lateral eyes. Zoology (Jena). 2005;108(1):47–53.
Prpic NM, Schoppmeier M, Damen WG. Collection and fixation of spider embryos. CSH Protoc. 2008. https://doi.org/10.1101/pdb.prot5067.
Samadi L, Schmid A, Eriksson BJ. Differential expression of retinal determination genes in the principal and secondary eyes of Cupiennius salei Keyserling (1877). EvoDevo. 2015;28(6):16.
Satija R, Farrell JA, Gennert D, Schier AF, Regev A. Spatial reconstruction of single-cell gene expression data. Nat Biotechnol. 2015;33(5):495–502.
Schomburg C, Turetzek N, Schacht MI, Schneider J, Kirfel P, Prpic NM, Posnien N. Molecular characterization and embryonic origin of the eyes in the common house spider Parasteatoda tepidariorum. EvoDevo. 2015;28(6):15.
Schwager EE, Sharma PP, Clarke T, Leite DJ, Wierschin T, Pechmann M, Akiyama-Oda Y, Esposito L, Bechsgaard J, Bilde T, Buffry AD, Chao H, Dinh H, Doddapaneni H, Dugan S, Eibner C, Extavour CG, Funch P, Garb J, Gonzalez LB, Gonzalez VL, Griffiths-Jones S, Han Y, Hayashi C, Hilbrant M, Hughes DST, Janssen R, Lee SL, Maeso I, Murali SC, Muzny DM, Nunes da Fonseca R, Paese CLB, Qu J, Ronshaugen M, Schomburg C, Schönauer A, Stollewerk A, Torres-Oliva M, Turetzek N, Vanthournout B, Werren JH, Wolff C, Worley KC, Bucher G, Gibbs RA, Coddington J, Oda H, Stanke M, Ayoub NA, Prpic NM, Flot JF, Posnien N, Richards S, McGregor AP. The house spider genome reveals an ancient whole-genome duplication during arachnid evolution. BMC Biol. 2017;15(1):62.
Sharma PP, Schwager EE, Extavour CG, Giribet G. Hox gene expression in the harvestman Phalangium opilio reveals divergent patterning of the chelicerate opisthosoma. Evol Dev. 2012;14(5):450–63.
Sharma PP, Kaluziak ST, Pérez-Porro AR, González VL, Hormiga G, Wheeler WC, Giribet G. Phylogenomic interrogation of arachnida reveals systemic conflicts in phylogenetic signal. Mol Biol Evol. 2014;31(11):2963–84.
Shkumatava A, Fischer S, Müller F, Strahle U, Neumann CJ. Sonic hedgehog, secreted by amacrine cells, acts as a short-range signal to direct differentiation and lamination in the zebrafish retina. Development. 2004;131(16):3849–58.
Simonnet F, Deutsch J, Quéinnec E. hedgehog is a segment polarity gene in a crustacean and a chelicerate. Dev Genes Evol. 2004;214(11):537–45.
Steinegger M, Söding J. MMseqs2 enables sensitive protein sequence searching for the analysis of massive data sets. Nat Biotechnol. 2017;35(11):1026–8.
Steinhoff POM, Harzsch S, Uhl G. Comparative neuroanatomy of the central nervous system in web-building and cursorial hunting spiders. J Comp Neurol. 2023. https://doi.org/10.1002/cne.25554.
Strausfeld NJ, Ma X, Edgecombe GD, Fortey RA, Land MF, Liu Y, Cong P, Hou X. Arthropod eyes: the early Cambrian fossil record and divergent evolution of visual systems. Arthropod Struct Dev. 2016;45(2):152–72.
Sturtevant MA, O’Neill JW, Bier E. Down-regulation of Drosophila Egf-r mRNA levels following hyperactivated receptor signaling. Development. 1994;120(9):2593–600.
Treisman JE. A conserved blueprint for the eye? BioEssays. 1999;21(10):843–50.
Treisman JE. Retinal differentiation in drosophila. Wiley Interdiscip Rev Dev Biol. 2013;2(4):545–57.
Turetzek N, Pechmann M, Schomburg C, Schneider J, Prpic NM. Neofunctionalization of a duplicate dachshund gene underlies the evolution of a novel leg segment in arachnids. Mol Biol Evol. 2016;33(1):109–21.
Tsachaki M, Sprecher SG. Genetic and developmental mechanisms underlying the formation of the Drosophila compound eye. Dev Dyn. 2012;241(1):40–56.
Vizueta J, Escuer P, Frías-López C, Guirao-Rico S, Hering L, Mayer G, Rozas J, Sánchez-Gracia A. Evolutionary history of major chemosensory gene families across panarthropoda. Mol Biol Evol. 2020;37(12):3601–15.
Voolstra O, Kiefer C, Hoehne M, Welsch R, Vogt K, von Lintig J. The drosophila class B scavenger receptor NinaD-I is a cell surface receptor mediating carotenoid transport for visual chromophore synthesis. Biochemistry. 2006;45(45):13429–37.
Wang T, Jiao Y, Montell C. Dissection of the pathway required for generation of vitamin A and for drosophila phototransduction. J Cell Biol. 2007;177(2):305–16.
Weller M, Tautz D. Prospero and snail expression during spider neurogenesis. Dev Genes Evol. 2003;213(11):554–66.
Wu T, Hu E, Xu S, Chen M, Guo P, Dai Z, Feng T, Zhou L, Tang W, Zhan L, Fu X, Liu S, Bo X, Yu G. clusterProfiler 4.0: a universal enrichment tool for interpreting omics data. Innovation. 2021;2(3):100141.
Zhang M, Xu Y, Li L, Wei S, Zhang S, Liu Z. Identification, evolution and expression of a CD36 homolog in the basal chordate amphioxus Branchiostoma japonicum. Fish Shellfish Immunol. 2013;34(2):546–55.
Zhang N (née Turetzek). Developmental studies of appendage patterning and formation in spiders. PhD thesis. 2016. Georg-August University Göttingen, Germany.
Funding
Open access funding provided by Uppsala University. Funding was provided to RJ and GEB by the Swedish Natural Science Council (VR) (grant no. 2022-03522), and the Marie Skłodowska-Curie Action (MSCA), Innovative Training Network (ITN), H20202-MSCAITN-2017 “EvoCELL” (grant no. 766053). This work was supported in part by the Deutsche Forschungsgemeinschaft priority program “The genomic basis of evolutionary innovations” (SPP2349; Project No. 503325252 awarded to MP [PE 2075/4-1] and NP [PO 1648/6-1]).
Author information
Authors and Affiliations
Contributions
RJ, GEB and BIM-J designed the study. BIM-J performed the single-cell-sequencing study. BIM-J and RJ performed in situ hybridization experiments for all spiders except Ischnothele caudata. RJ assembled the final figures. RJ and BIM-J wrote the first draft of the manuscript. MP provided bioinformatic data and embryos of Acanthoscurria geniculata. MP and NP provided all data concerning Ischnothele caudata, such as bioinformatic data, in situ hybridizations and figure panels. All authors discussed the data and contributed to the writing of the final version of the manuscript.
Corresponding author
Ethics declarations
Competing interests
The authors declare no competing interests.
Additional information
Publisher's Note
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
Supplementary Information
13227_2024_230_MOESM9_ESM.tif
Additional file 9: Phylogenetic tree of CD36-family genes. Bayesian analysis using MrBayes applying 0.5 million cycles for the Metropolis-Coupled Markov Chain Monte Carlo (MCMCMC). The tree is midpoint rooted. Node labels represent posterior possibilities. The scale bar represents 0.3 amino acid substitutions per site. The identified CD36.1 gene is marked with a red asterisk. Note that this gene appears to possess a second paralog in the spider (XP021004437.1). The sequence of Parasteatoda CD36.1. Is most similar to that of Drosophila Group2A and B genes but does not represent a one-to-one ortholog of either of these six Drosophila CD36-family genes (see text for further information). Species abbreviations: Dm, Drosophila melanogaster; Dr, Danio rerio; Mm, Mus musculus; Pt, Parasteatoda tepidariorum
13227_2024_230_MOESM10_ESM.tif
Additional file 10: Additional gene expression data. Expression of unc0466 (A-D), opa (E-H), lim1 (I, J), lim2 (K, L), unc4232 (M, N), scro (O, P), unc3221 (Q), unc7000 (R), gsc (S), mmr1 (T, U), mmr2 (V), and FMRFa-R (W). In all panels, anterior is to the left, except panel W (dissected leg, dorsal to the left; distal end pointing downwards). Arrows in panel W point to expression inside the leg. Abbreviations: aSp, anterior spinneret; bl, book lung; ch, chelicera; h, heart; hl, head lobes; L, leg; pSp, posterior spinneret; s, stomodaeum.
13227_2024_230_MOESM15_ESM.tif
Additional file 15: Expression of spider hh1 genes. In all panels, anterior is to the left. Panels A, C, D, F, G, H, L, M, N, O and R represent lateral views. Panels B, E, I, J, P and Q represent lateral views (in these panels, dorsal is up). Panel K represents a dorsal view. Note the segment-polarity gene-like expression of all hh1-group genes, expression along the appendages, expression in the segment addition zone (saz), and expression in the anlagen of the stomodaeum (asterisks in all panels). The arrow in panel O points to four dots of expression in the neuronal ectoderm that thus are not associated with the developing eyes. In all panels, full circles mark expression in the pre-cheliceral region, but note that this expression is not at the place where the eyes will form. Panels B´, H´, I´, L´, O´ and P´ represent Sybr-green staining of corresponding embryos. Developmental stages are indicated after Mittmann and Wolff [61]; note that developmental stages in different species of spiders have been defined by comparable morphological landmarks such as the overall shape of the embryo and the length of the appendages. Abbreviations: ch, cheliceral segment; L, leg-bearing segment; pp, pedipalp-bearing segment; saz, segment addition zone.
13227_2024_230_MOESM16_ESM.tif
Additional file 16: Expression of Parasteatoda hh2-r (panels A and B), omb2 (panels C-E), omb3 (panels F-I), and omb4 (panels J and K). In all panels, anterior is to the left. The asterisks in panels C, D, and F-H mark expression in the developing labrum. The primordia of the secondary eyes are encircled in one half of the head lobes in panels D, E, J, and K. Note that there is no expression of omb2 and omb4. The dot-like expression seen in panels A and B are in the yolk and the dorsal field, but not in the embryo proper. Abbreviations: hl, head lobe; saz, segment addition zone.
Rights and permissions
Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http://creativecommons.org/licenses/by/4.0/.
About this article
Cite this article
Medina‑Jiménez, B.I., Budd, G.E., Pechmann, M. et al. Single-cell sequencing suggests a conserved function of Hedgehog-signalling in spider eye development. EvoDevo 15, 11 (2024). https://doi.org/10.1186/s13227-024-00230-6
Received:
Accepted:
Published:
DOI: https://doi.org/10.1186/s13227-024-00230-6