A comprehensive reference transcriptome resource for the common house spider Parasteatoda tepidariorum

Evolutionary origin and distribution of leucine-rich repeat-containing G protein-coupled receptors in insects

Leucine-rich repeat-containing G protein-coupled receptors (LGRs) are crucial for animal growth and development. They were categorized into four types (A, B, C1, and C2) based on their sequence and domain structures. Despite the widespread distribution of LGRs across bilaterians, a thorough investigation of their distribution and evolutionary history remains elusive. Recent studies insect LGRs, especially the emergence of type C2 LGRs in various hemimetabolous insects, had prompted our study to address these problems. Initially, we traced the origins of LGRs by exploiting data from 99 species spanning 11 metazoan phyla, and discovered that type A and B LGRs originated from sponges, while type C LGRs originated from cnidarians. Subsequently, through comprehensive genomic and transcriptomic analyses across 565 species across 25 orders of insects, we found that both type A and C1 LGRs divided into two gene clusters. These clusters can be traced back to basal Insecta and an early ancestor of the Arthropoda, respectively. Furthermore, the absence of type B LGRs in wingless insects suggests a role in wing development, while the absence of type C2 LGRs in holometabolous insects hints at novel functions unrelated to insect metamorphosis. According to the origin of LGRs and the investigation of LGRs in insects, we speculated that type A and B LGRs appeared first among four types of LGRs, type A evolved into type C LGRs later, and type A and C1 LGRs independently duplicated during the evolutionary process. This study provides a more comprehensive view of the evolution of LGR genes than previously available, and sheds light on the evolutionary history and significance of LGRs in insect biology.

A taxon-restricted duplicate of Iroquois3 is required for patterning the spider waist

The chelicerate body plan is distinguished from other arthropod groups by its division of segments into 2 tagmata: the anterior prosoma ("cephalothorax") and the posterior opisthosoma ("abdomen"). Little is understood about the genetic mechanisms that establish the prosomal-opisthosomal (PO) boundary. To discover these mechanisms, we created high-quality genomic resources for the large-bodied spider Aphonopelma hentzi. We sequenced specific territories along the antero-posterior axis of developing embryos and applied differential gene expression analyses to identify putative regulators of regional identity. After bioinformatic screening for candidate genes that were consistently highly expressed in only 1 tagma (either the prosoma or the opisthosoma), we validated the function of highly ranked candidates in the tractable spider model Parasteatoda tepidariorum. Here, we show that an arthropod homolog of the Iroquois complex of homeobox genes is required for proper formation of the boundary between arachnid tagmata. The function of this homolog had not been previously characterized, because it was lost in the common ancestor of Pancrustacea, precluding its investigation in well-studied insect model organisms. Knockdown of the spider copy of this gene, which we designate as waist-less, in P. tepidariorum resulted in embryos with defects in the PO boundary, incurring discontinuous spider germ bands. We show that waist-less is required for proper specification of the segments that span the prosoma-opisthosoma boundary, which in adult spiders corresponds to the narrowed pedicel. Our results demonstrate the requirement of an ancient, taxon-restricted paralog for the establishment of the tagmatic boundary that defines Chelicerata.

De novo assembly of transcriptomes and differential gene expression analysis using short-read data from emerging model organisms - a brief guide

Many questions in biology benefit greatly from the use of a variety of model systems. High-throughput sequencing methods have been a triumph in the democratization of diverse model systems. They allow for the economical sequencing of an entire genome or transcriptome of interest, and with technical variations can even provide insight into genome organization and the expression and regulation of genes. The analysis and biological interpretation of such large datasets can present significant challenges that depend on the 'scientific status' of the model system. While high-quality genome and transcriptome references are readily available for well-established model systems, the establishment of such references for an emerging model system often requires extensive resources such as finances, expertise and computation capabilities. The de novo assembly of a transcriptome represents an excellent entry point for genetic and molecular studies in emerging model systems as it can efficiently assess gene content while also serving as a reference for differential gene expression studies. However, the process of de novo transcriptome assembly is non-trivial, and as a rule must be empirically optimized for every dataset. For the researcher working with an emerging model system, and with little to no experience with assembling and quantifying short-read data from the Illumina platform, these processes can be daunting. In this guide we outline the major challenges faced when establishing a reference transcriptome de novo and we provide advice on how to approach such an endeavor. We describe the major experimental and bioinformatic steps, provide some broad recommendations and cautions for the newcomer to de novo transcriptome assembly and differential gene expression analyses. Moreover, we provide an initial selection of tools that can assist in the journey from raw short-read data to assembled transcriptome and lists of differentially expressed genes.

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.

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.

FGF signalling is involved in cumulus migration in the common house spider Parasteatoda tepidariorum

Cell migration is a fundamental component during the development of most multicellular organisms. In the early spider embryo, the collective migration of signalling cells, known as the cumulus, is required to set the dorsoventral body axis. Here, we show that FGF signalling plays an important role during cumulus migration in the spider Parasteatoda tepidariorum. Spider embryos with reduced FGF signalling show reduced or absent cumulus migration and display dorsoventral patterning defects. Our study reveals that the transcription factor Ets4 regulates the expression of several FGF signalling components in the cumulus. In conjunction with a previous study, we show that the expression of fgf8 in the germ-disc is regulated via the Hedgehog signalling pathway. We also demonstrate that FGF signalling influences the BMP signalling pathway activity in the region around cumulus cells. Finally, we show that FGFR signalling might also influence cumulus migration in basally branching spiders and we propose that fgf8 might act as a chemo-attractant to guide cumulus cells towards the future dorsal pole of the spider embryo.

Female mating rates and their fitness consequences in the common house spider Parasteatoda tepidariorum

Mating systems, with varying female mating rates occurring with the same partner (monandry) or with multiple mates (polyandry), can have far reaching consequences for population viability and the rate of gene flow. Here, we investigate the mating rates of the common house spider Parasteatoda tepidariorum (Theridiidae), an emerging model for genetic studies, with yet undescribed reproductive behavior. It is hypothesized that spiders belonging to this family have low re-mating rates. We paired females twice with the same male (monandry) or with different males (polyandry), and recorded behaviors, mating success and fitness resulting from single- and double-matings, either monandrous or polyandrous. Despite the study being explorative in nature, we predict successful matings to be more frequent during first encounters, to reduce female risk of remaining unmated. For re-mating to be adaptive, we expect higher fitness of double-mated females, and polyandrous females to experience highest mating success and fitness if reproductive gains are achieved by mating with multiple partners. We show that the majority of the females did not mate, and those that did mated only once, not necessarily on their first encounter. The likelihood of re-mating did not differ between monandrous and polyandrous encounters and female mating experience (mated once, twice monandrous, twice polyandrous) did not affect fitness, indicated by similar offspring production. Female twanging of the web leads to successful matings suggesting female behavioral receptivity. Cannibalism rates were low and mostly occurred pre-copulatory. We discuss how the species ecology, with potentially high mating costs for males and limited female receptivity, may shape a mating system with low mating rates.

Lineage-specific, fast-evolving GATA-like gene regulates zygotic gene activation to promote endoderm specification and pattern formation in the Theridiidae spider

Background

The process of early development varies across the species-rich phylum Arthropoda. Owing to the limited research strategies for dissecting lineage-specific processes of development in arthropods, little is known about the variations in early arthropod development at molecular resolution. The Theridiidae spider, Parasteatoda tepidariorum, has its genome sequenced and could potentially contribute to dissecting early embryonic processes.

Results

We present genome-wide identification of candidate genes that exhibit locally restricted expression in germ disc forming stage embryos of P. tepidariorum, based on comparative transcriptomes of isolated cells from different regions of the embryo. A subsequent pilot screen by parental RNA interference identifies three genes required for body axis formation. One of them is a GATA-like gene that has been fast evolving after duplication and divergence from a canonical GATA family gene. This gene is designated fuchi nashi (fuchi) after its knockdown phenotypes, where the cell movement toward the formation of a germ disc was reversed. fuchi expression occurs in cells outside a forming germ disc and persists in the endoderm. Transcriptome and chromatin accessibility analyses of fuchi pRNAi embryos suggest that early fuchi activity regulates chromatin state and zygotic gene activation to promote endoderm specification and pattern formation. We also show that there are many uncharacterized genes regulated by fuchi.

Conclusions

Our genome-based research using an arthropod phylogenetically distant from Drosophila identifies a lineage-specific, fast-evolving gene with key developmental roles in one of the earliest, genome-wide regulatory events, and allows for molecular exploration of the developmental variations in early arthropod embryos.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12915-022-01421-0.

The organizing role of Wnt signaling pathway during arthropod posterior growth

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

Expression of spidroin proteins in the silk glands of golden orb-weaver spiders

The expression of spidroins in the major ampullate, minor ampullate, flagelliform, and tubuliform silk glands of Trichonephila clavipes spiders was analyzed using proteomics analysis techniques. Spidroin peptides were identified and assigned to different gene products based on sequence concurrence when compared with the whole genome of the spider. It was found that only a relatively low proportion of the spidroin genes are expressed as proteins in any of the studied glands. In addition, the expression of spidroin genes in different glands presents a wide range of patterns, with some spidroins being found in a single gland exclusively, while others appear in the content of several glands. The combination of precise genomics, proteomics, microstructural, and mechanical data provides new insights both on the design principles of these materials and how these principles might be translated for the production of high-performance bioinspired artificial fibers.

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.

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

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

Identification and activity of monoamine oxidase in the orb-weaving spider Larinioides cornutus

Monoamine oxidase (MAO) is a mitochondrial membrane-bound enzyme that catalyzes the oxidative deamination of monoamines in a wide array of organisms. While the enzyme monoamine oxidase has been studied extensively in its role in moderating behavior in mammals, there is a paucity of research investigating this role in invertebrates, where the latter utilizes this enzyme in a major pathway to degrade monoamines. There is especially a dismal lack of information on how MAO influences activity in invertebrates, particularly in account of the circadian cycle. Previous studies revealed MAO degrades serotonin and norepinephrine in arachnids, but did not investigate other critically important compounds like octopamine. Larinioides cornutus is a species of orb-weaving spider that exhibits diel fluctuations in behavior, specifically levels of aggression. The monoamines octopamine and serotonin have been shown to influence aggressive behaviors in L. cornutus, thus this species was used to investigate if MAO is a potential site of regulation throughout the day. Not only did gene expression of MAO orthologs and MAO activity fluctuate at different times of day, but the enzymatic activity was substrate-specific producing a higher level of degradation of octopamine as compared to serotonin in vitro. This study further supports evidence that MAO has an active role in monoamine inactivation in invertebrates and provides a first look at how MAO ultimately may be regulating behavior in an invertebrate.

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.

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.

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.

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.

Candidate gene screen for potential interaction partners and regulatory targets of the Hox gene labial in the spider Parasteatoda tepidariorum

The Hox gene labial (lab) governs the formation of the tritocerebral head segment in insects and spiders. However, the morphology that results from lab action is very different in the two groups. In insects, the tritocerebral segment (intercalary segment) is reduced and lacks appendages, whereas in spiders the corresponding segment (pedipalpal segment) is a proper segment including a pair of appendages (pedipalps). It is likely that this difference between lab action in insects and spiders is mediated by regulatory targets or interacting partners of lab. However, only a few such genes are known in insects and none in spiders. We have conducted a candidate gene screen in the spider Parasteatoda tepidariorum using as candidates Drosophila melanogaster genes known to (potentially) interact with lab or to be expressed in the intercalary segment. We have studied 75 P. tepidariorum genes (including previously published and duplicated genes). Only 3 of these (proboscipedia-A (pb-A) and two paralogs of extradenticle (exd)) showed differential expression between leg and pedipalp. The low success rate points to a weakness of the candidate gene approach when it is applied to lineage specific organs. The spider pedipalp has no counterpart in insects, and therefore relying on insect data apparently cannot identify larger numbers of factors implicated in its specification and formation. We argue that in these cases a de novo approach to gene discovery might be superior to the candidate gene approach.

Ovarian Transcriptomic Analyses in the Urban Human Health Pest, the Western Black Widow Spider

Due to their abundance and ability to invade diverse environments, many arthropods have become pests of economic and health concern, especially in urban areas. Transcriptomic analyses of arthropod ovaries have provided insight into life history variation and fecundity, yet there are few studies in spiders despite their diversity within arthropods. Here, we generated a de novo ovarian transcriptome from 10 individuals of the western black widow spider (Latrodectus hesperus), a human health pest of high abundance in urban areas, to conduct comparative ovarian transcriptomic analyses. Biological processes enriched for metabolism—specifically purine, and thiamine metabolic pathways linked to oocyte development—were significantly abundant in L. hesperus. Functional and pathway annotations revealed overlap among diverse arachnid ovarian transcriptomes for highly-conserved genes and those linked to fecundity, such as oocyte maturation in vitellogenin and vitelline membrane outer layer proteins, hormones, and hormone receptors required for ovary development, and regulation of fertility-related genes. Comparative studies across arachnids are greatly needed to understand the evolutionary similarities of the spider ovary, and here, the identification of ovarian proteins in L. hesperus provides potential for understanding how increased fecundity is linked to the success of this urban pest.

Expression profile of genes encoding allatoregulatory neuropeptides in females of the spider Parasteatoda tepidariorum (Araneae, Theridiidae)

Allatoregulatory neuropeptides are multifunctional proteins that take part in the synthesis and secretion of juvenile hormones. In insects, allatostatins are inhibitors of juvenile hormone biosynthesis in the corpora allata while allatotropins, act as stimulators. By quantitative real-time PCR, we analyzed the gene expression of allatostatin A (PtASTA), allatostatin B (PtASTB), allatostatin C (PtASTC), allatotropin (PtAT) and their receptors (PtASTA-R, PtASTB-R, PtASTC-R, PtAT-R) in various tissues in different age groups of female spiders. In the presented manuscript, the presence of allatostatin A, allatostatin C, and allatotropin are reported in females of the spider P. tepidariorum. The obtained results indicated substantial differences in gene expression levels for allatoregulatory neuropeptides and their receptors in the different tissues. Additionally, the gene expression levels also varied depending on the female age. Strong expression was observed coinciding with sexual maturation in the neuroendocrine and nervous system, and to a lower extent in the digestive tissues and ovaries. Reverse trends were observed for the expression of genes encoding the receptors of these neuropeptides. In conclusion, our study is the first hint that the site of synthesis and secretion is fulfilled by similar structures as observed in other arthropods. In addition, the results of the analysis of spider physiology give evidence that the general functions like regulation of the juvenile hormone synthesis, regulation of the digestive tract and ovaries action, control of vitellogenesis process by the neuropeptides seem to be conserved among arthropods and are the milestone to future functional studies.

Niemann-Pick proteins type C2 are identified as olfactory related genes of Pardosa pseudoannulata by transcriptome and expression profile analysis

In arthropods, the large majority of studies on olfaction have been mainly focused on insects, whereas little on Arachnida, even though olfaction is very important in arachnid behavior. Pardosa pseudoannulata is one of the most common wandering spiders in rice fields, as the important natural enemy against a range of pests. However, little is known about the potential chemosensory proteins involved in olfactory behavior of these spiders. Niemann-Pick proteins type C2 (NPC2) as a new class of binding and transport proteins for semiochemicals in arthropods especially ticks and mites has received more attention in recent years. In this study, six NPC2s namely PpseNPC1-6 were newly identified in the appendages of P. pseudoannulata based on transcriptome data. A phylogenetic analysis indicated that all of P. pseudoannulata NPC2s were clustered together forming one clade with high posterior probability values. In addition, the sequences shared the same subclade with the NPC2 sequences of ticks and scorpion. The motif-patterns indicated that PpseNPC2-5 had the common pattern with the two-spotted spider mite Tetranychus urticae and the ant Trachymyrmex cornetzi. Furthermore, quantitative real-time PCR (qPCR) measurements were conducted to evaluate the expression profile of these genes in various tissues of P. pseudoannulata. It was found that most NPC2s (PpseNPC2-1, PpseNPC2-2, PpseNPC2-5 and PpseNPC2-6) were highly expressed in adult pedipalps and chelicerae. Owing to the functional olfactory organs in Chelicerata of pedipalps, our results supported a putative role of NPC2s as new odorant carriers in P. pseudoannulata.

A proteotranscriptomic study of silk-producing glands from the orb-weaving spiders

A proteotranscriptomic approach provides a biochemical basis for understanding the intricate spinning process and complex structural features of spider silk proteins.

The evolution of selector gene function: Expression dynamics and regulatory interactions of tiptop/teashirt across Arthropoda

The transcription factors spineless (ss) and tiptop/teashirt (tio/tsh) have been shown to be selectors of distal appendage identity in an insect, but it is unknown how they regulate one another. Here, we examined the regulatory relationships between these two determinants in the milkweed bug Oncopeltus faciatus, using maternal RNA interference (RNAi). We show that Ofas-ss RNAi embryos bear distally transformed antennal buds with heterogeneous Ofas-tio/tsh expression domains comparable to wild type legs. In the reciprocal experiment, Ofas-tio/tsh RNAi embryos bear distally transformed walking limb buds with ectopic expression of Ofas-ss in the distal leg primordia. These data suggest that Ofas-ss is required for the maintenance of Ofas-tio/tsh expression in the distal antenna, whereas Ofas-tio/tsh represses Ofas-ss in the leg primordia. To assess whether expression boundaries of tio/tsh are associated with the trunk region more generally, we surveyed the expression of one myriapod and two chelicerate tio/tsh homologs. Our expression survey suggests that tio/tsh could play a role in specifying distal appendage identity across Arthropoda, but Hox regulation of tio/tsh homologs has been evolutionarily labile.

An overview of Phoneutria nigriventer spider venom using combined transcriptomic and proteomic approaches

Phoneutria nigriventer is one of the largest existing true spiders and one of the few considered medically relevant. Its venom contains several neurotoxic peptides that act on different ion channels and chemical receptors of vertebrates and invertebrates. Some of these venom toxins have been shown as promising models for pharmaceutical or biotechnological use. However, the large diversity and the predominance of low molecular weight toxins in this venom have hampered the identification and deep investigation of the less abundant toxins and the proteins with high molecular weight. Here, we combined conventional and next-generation cDNA sequencing with Multidimensional Protein Identification Technology (MudPIT), to obtain an in-depth panorama of the composition of P. nigriventer spider venom. The results from these three approaches showed that cysteine-rich peptide toxins are the most abundant components in this venom and most of them contain the Inhibitor Cysteine Knot (ICK) structural motif. Ninety-eight sequences corresponding to cysteine-rich peptide toxins were identified by the three methodologies and many of them were considered as putative novel toxins, due to the low similarity to previously described toxins. Furthermore, using next-generation sequencing we identified families of several other classes of toxins, including CAPs (Cysteine Rich Secretory Protein-CRiSP, antigen 5 and Pathogenesis-Related 1-PR-1), serine proteinases, TCTPs (translationally controlled tumor proteins), proteinase inhibitors, metalloproteinases and hyaluronidases, which have been poorly described for this venom. This study provides an overview of the molecular diversity of P. nigriventer venom, revealing several novel components and providing a better basis to understand its toxicity and pharmacological activities.

Expression and function of the zinc finger transcription factor Sp6-9 in the spider Parasteatoda tepidariorum

Zinc finger transcription factors of the Sp6-9 group are evolutionarily conserved in all metazoans and have important functions in, e.g., limb formation and heart development. The function of Sp6-9-related genes has been studied in a number of vertebrates and invertebrates, but data from chelicerates (spiders and allies) was lacking so far. We have isolated the ortholog of Sp6-9 from the common house spider Parasteatoda tepidariorum and the cellar spider Pholcus phalangioides. We show that the Sp6-9 gene in these spider species is expressed in the developing appendages thus suggesting a conserved role in limb formation. Indeed, RNAi with Sp6-9 in P. tepidariorum leads not only to strong limb defects, but also to the loss of body segments and head defects in more strongly affected animals. Together with a new expression domain in the early embryo, these data suggest that Sp6-9 has a dual role P. tepidariorum. The early role in head and body segment formation is not known from other arthropods, but the role in limb formation is evolutionarily highly conserved.

Evolution of Chemosensory Gene Families in Arthropods: Insight from the First Inclusive Comparative Transcriptome Analysis across Spider Appendages

Unlike hexapods and vertebrates, in chelicerates, knowledge of the specific molecules involved in chemoreception comes exclusively from the comparative analysis of genome sequences. Indeed, the genomes of mites, ticks and spiders contain several genes encoding homologs of some insect membrane receptors and small soluble chemosensory proteins. Here, we conducted for the first time a comprehensive comparative RNA-Seq analysis across different body structures of a chelicerate: the nocturnal wandering hunter spider Dysdera silvatica Schmidt 1981. Specifically, we obtained the complete transcriptome of this species as well as the specific expression profile in the first pair of legs and the palps, which are thought to be the specific olfactory appendages in spiders, and in the remaining legs, which also have hairs that have been morphologically identified as chemosensory. We identified several ionotropic (Ir) and gustatory (Gr) receptor family members exclusively or differentially expressed across transcriptomes, some exhibiting a distinctive pattern in the putative olfactory appendages. Furthermore, these IRs were the only known olfactory receptors identified in such structures. These results, integrated with an extensive phylogenetic analysis across arthropods, uncover a specialization of the chemosensory gene repertoire across the body of D. silvatica and suggest that some IRs likely mediate olfactory signaling in chelicerates. Noticeably, we detected the expression of a gene family distantly related to insect odorant-binding proteins (OBPs), suggesting that this gene family is more ancient than previously believed, as well as the expression of an uncharacterized gene family encoding small globular secreted proteins, which appears to be a good chemosensory gene family candidate.

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.

Rapid diversification of homothorax expression patterns after gene duplication in spiders

Background

Gene duplications provide genetic material for the evolution of new morphological and physiological features. One copy can preserve the original gene functions while the second copy may evolve new functions (neofunctionalisation). Gene duplications may thus provide new genes involved in evolutionary novelties.

Results

We have studied the duplicated homeobox gene homothorax (hth) in the spider species Parasteatoda tepidariorum and Pholcus phalangioides and have compared these data with previously published data from additional spider species. We show that the expression pattern of hth1 is highly conserved among spiders, consistent with the notion that this gene copy preserves the original hth functions. By contrast, hth2 has a markedly different expression profile especially in the prosomal appendages. The pattern in the pedipalps and legs consists of several segmental rings, suggesting a possible role of hth2 in limb joint development. Intriguingly, however, the hth2 pattern is much less conserved between the species than hth1 and shows a species specific pattern in each species investigated so far.

Conclusions

We hypothesise that the hth2 gene has gained a new patterning function after gene duplication, but has then undergone a second phase of diversification of its new role in the spider clade. The evolution of hth2 may thus provide an interesting example for a duplicated gene that has not only contributed to genetic diversity through neofunctionalisation, but beyond that has been able to escape evolutionary conservation after neofunctionalisation thus forming the basis for further genetic diversification.

Electronic supplementary material

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

Molecular control of gut formation in the spider Parasteatoda tepidariorum

The development of a digestive system is an essential feature of bilaterians. Studies of the molecular control of gut formation in arthropods have been studied in detail in the fruit fly Drosophila melanogaster. However, little is known in other arthropods, especially in noninsect arthropods. To better understand the evolution of arthropod alimentary system, we investigate the molecular control of gut development in the spider Parasteatoda tepidariorum (Pt), the primary chelicerate model species for developmental studies. Orthologs of the ectodermal genes Pt-wingless (Pt-wg) and Pt-hedgehog (Pt-hh), of the endodermal genes, Pt-serpent (Pt-srp) and Pt-hepatocyte-nuclear factor-4 (Pt-hnf4) and of the mesodermal gene Pt-twist (Pt-twi) are expressed in the same germ layers during spider gut development as in D. melanogaster. Thus, our expression data suggest that the downstream molecular components involved in gut development in arthropods are conserved. However, Pt-forkhead (Pt-fkh) expression and function in spiders is considerably different from its D. melanogaster ortholog. Pt-fkh is expressed before gastrulation in a cell population that gives rise to endodermal and mesodermal precursors, suggesting a possible role for this factor in specification of both germ layers. To test this hypothesis, we knocked down Pt-fkh via RNA interference. Pt-fkh RNAi embryos not only fail to develop a proper gut, but also lack the mesodermal Pt-twi expressing cells. Thus, in spiders Pt-fkh specifies endodermal and mesodermal germ layers. We discuss the implications of these findings for the evolution and development of gut formation in Ecdysozoans.

Comparative transcriptomics of Entelegyne spiders (Araneae, Entelegynae), with emphasis on molecular evolution of orphan genes

Next-generation sequencing technology is rapidly transforming the landscape of evolutionary biology, and has become a cost-effective and efficient means of collecting exome information for non-model organisms. Due to their taxonomic diversity, production of interesting venom and silk proteins, and the relative scarcity of existing genomic resources, spiders in particular are excellent targets for next-generation sequencing (NGS) methods. In this study, the transcriptomes of six entelegyne spider species from three genera (Cicurina travisae, C. vibora, Habronattus signatus, H. ustulatus, Nesticus bishopi, and N. cooperi) were sequenced and de novo assembled. Each assembly was assessed for quality and completeness and functionally annotated using gene ontology information. Approximately 100 transcripts with evidence of homology to venom proteins were discovered. After identifying more than 3,000 putatively orthologous genes across all six taxa, we used comparative analyses to identify 24 instances of positively selected genes. In addition, between ~ 550 and 1,100 unique orphan genes were found in each genus. These unique, uncharacterized genes exhibited elevated rates of amino acid substitution, potentially consistent with lineage-specific adaptive evolution. The data generated for this study represent a valuable resource for future phylogenetic and molecular evolutionary research, and our results provide new insight into the forces driving genome evolution in taxa that span the root of entelegyne spider phylogeny.

House spider genome uncovers evolutionary shifts in the diversity and expression of black widow venom proteins associated with extreme toxicity

BACKGROUND

Black widow spiders are infamous for their neurotoxic venom, which can cause extreme and long-lasting pain. This unusual venom is dominated by latrotoxins and latrodectins, two protein families virtually unknown outside of the black widow genus Latrodectus, that are difficult to study given the paucity of spider genomes. Using tissue-, sex- and stage-specific expression data, we analyzed the recently sequenced genome of the house spider (Parasteatoda tepidariorum), a close relative of black widows, to investigate latrotoxin and latrodectin diversity, expression and evolution.

RESULTS

We discovered at least 47 latrotoxin genes in the house spider genome, many of which are tandem-arrayed. Latrotoxins vary extensively in predicted structural domains and expression, implying their significant functional diversification. Phylogenetic analyses show latrotoxins have substantially duplicated after the Latrodectus/Parasteatoda split and that they are also related to proteins found in endosymbiotic bacteria. Latrodectin genes are less numerous than latrotoxins, but analyses show their recruitment for venom function from neuropeptide hormone genes following duplication, inversion and domain truncation. While latrodectins and other peptides are highly expressed in house spider and black widow venom glands, latrotoxins account for a far smaller percentage of house spider venom gland expression.

CONCLUSIONS

The house spider genome sequence provides novel insights into the evolution of venom toxins once considered unique to black widows. Our results greatly expand the size of the latrotoxin gene family, reinforce its narrow phylogenetic distribution, and provide additional evidence for the lateral transfer of latrotoxins between spiders and bacterial endosymbionts. Moreover, we strengthen the evidence for the evolution of latrodectin venom genes from the ecdysozoan Ion Transport Peptide (ITP)/Crustacean Hyperglycemic Hormone (CHH) neuropeptide superfamily. The lower expression of latrotoxins in house spiders relative to black widows, along with the absence of a vertebrate-targeting α-latrotoxin gene in the house spider genome, may account for the extreme potency of black widow venom.

Combining independent de novo assemblies optimizes the coding transcriptome for nonconventional model eukaryotic organisms

Background

Next-generation sequencing (NGS) technologies are arguably the most revolutionary technical development to join the list of tools available to molecular biologists since PCR. For researchers working with nonconventional model organisms one major problem with the currently dominant NGS platform (Illumina) stems from the obligatory fragmentation of nucleic acid material that occurs prior to sequencing during library preparation. This step creates a significant bioinformatic challenge for accurate de novo assembly of novel transcriptome data. This challenge becomes apparent when a variety of modern assembly tools (of which there is no shortage) are applied to the same raw NGS dataset. With the same assembly parameters these tools can generate markedly different assembly outputs.

Results

In this study we present an approach that generates an optimized consensus de novo assembly of eukaryotic coding transcriptomes. This approach does not represent a new assembler, rather it combines the outputs of a variety of established assembly packages, and removes redundancy via a series of clustering steps. We test and validate our approach using Illumina datasets from six phylogenetically diverse eukaryotes (three metazoans, two plants and a yeast) and two simulated datasets derived from metazoan reference genome annotations. All of these datasets were assembled using three currently popular assembly packages (CLC, Trinity and IDBA-tran). In addition, we experimentally demonstrate that transcripts unique to one particular assembly package are likely to be bioinformatic artefacts. For all eight datasets our pipeline generates more concise transcriptomes that in fact possess more unique annotatable protein domains than any of the three individual assemblers we employed. Another measure of assembly completeness (using the purpose built BUSCO databases) also confirmed that our approach yields more information.

Conclusions

Our approach yields coding transcriptome assemblies that are more likely to be closer to biological reality than any of the three individual assembly packages we investigated. This approach (freely available as a simple perl script) will be of use to researchers working with species for which there is little or no reference data against which the assembly of a transcriptome can be performed.

Electronic supplementary material

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

Neuropeptide Evolution: Chelicerate Neurohormone and Neuropeptide Genes may reflect one or more whole genome duplications

Four genomes and two transcriptomes from six Chelicerate species were analyzed for the presence of neuropeptide and neurohormone precursors and their GPCRs. The genome from the spider Stegodyphus mimosarum yielded 87 neuropeptide precursors and 101 neuropeptide GPCRs. High neuropeptide transcripts were also found in the trancriptomes of three other spiders, Latrodectus hesperus, Parasteatoda tepidariorum and Acanthoscurria geniculata. For the scorpion Mesobuthus martensii the numbers are 79 and 74 respectively. The very small genome of the house dust mite, Dermatophagoides farinae, on the other hand contains a much smaller number of such genes. A few new putative Arthropod neuropeptide genes were discovered. Thus, both spiders and the scorpion have an achatin gene and in spiders there are two different genes encoding myosuppressin-like peptides while spiders also have two genes encoding novel LGamides. Another finding is the presence of trissin in spiders and scorpions, while neuropeptide genes that seem to be orthologs of Lottia LFRYamide and Platynereis CCRFamide were also found. Such genes were also found in various insect species, but seem to be lacking from the Holometabola. The Chelicerate neuropeptide and neuropeptide GPCR genes often have paralogs. As the large majority of these are probably not due to local gene duplications, is not impossible that they reflect the effects of one or more ancient whole genome duplications.

Expression-Linked Patterns of Codon Usage, Amino Acid Frequency, and Protein Length in the Basally Branching Arthropod Parasteatoda tepidariorum

Spiders belong to the Chelicerata, the most basally branching arthropod subphylum. The common house spider, Parasteatoda tepidariorum, is an emerging model and provides a valuable system to address key questions in molecular evolution in an arthropod system that is distinct from traditionally studied insects. Here, we provide evidence suggesting that codon usage, amino acid frequency, and protein lengths are each influenced by expression-mediated selection in P. tepidariorum. First, highly expressed genes exhibited preferential usage of T3 codons in this spider, suggestive of selection. Second, genes with elevated transcription favored amino acids with low or intermediate size/complexity (S/C) scores (glycine and alanine) and disfavored those with large S/C scores (such as cysteine), consistent with the minimization of biosynthesis costs of abundant proteins. Third, we observed a negative correlation between expression level and coding sequence length. Together, we conclude that protein-coding genes exhibit signals of expression-related selection in this emerging, noninsect, arthropod model.

Pervasive microRNA Duplication in Chelicerates: Insights from the Embryonic microRNA Repertoire of the Spider Parasteatoda tepidariorum

MicroRNAs are small (∼22 nt) noncoding RNAs that repress translation and therefore regulate the production of proteins from specific target mRNAs. microRNAs have been found to function in diverse aspects of gene regulation within animal development and many other processes. Among invertebrates, both conserved and novel, lineage specific, microRNAs have been extensively studied predominantly in holometabolous insects such as Drosophila melanogaster However little is known about microRNA repertoires in other arthropod lineages such as the chelicerates. To understand the evolution of microRNAs in this poorly sampled subphylum, we characterized the microRNA repertoire expressed during embryogenesis of the common house spider Parasteatoda tepidariorum We identified a total of 148 microRNAs in P. tepidariorum representing 66 families. Approximately half of these microRNA families are conserved in other metazoans, while the remainder are specific to this spider. Of the 35 conserved microRNAs families 15 had at least two copies in the P. tepidariorum genome. A BLAST-based approach revealed a similar pattern of duplication in other spiders and a scorpion, but not among other chelicerates and arthropods, with the exception of a horseshoe crab. Among the duplicated microRNAs we found examples of lineage-specific tandem duplications, and the duplication of entire microRNA clusters in three spiders, a scorpion, and in a horseshoe crab. Furthermore, we found that paralogs of many P. tepidariorum microRNA families exhibit arm switching, which suggests that duplication was often followed by sub- or neofunctionalization. Our work shows that understanding the evolution of microRNAs in the chelicerates has great potential to provide insights into the process of microRNA duplication and divergence and the evolution of animal development.

Neuropeptide evolution: Chelicerate neurohormone and neuropeptide genes may reflect one or more whole genome duplications

Four genomes and two transcriptomes from six Chelicerate species were analyzed for the presence of neuropeptide and neurohormone precursors and their GPCRs. The genome from the spider Stegodyphus mimosarum yielded 87 neuropeptide precursors and 120 neuropeptide GPCRs. Many neuropeptide transcripts were also found in the transcriptomes of three other spiders, Latrodectus hesperus, Parasteatoda tepidariorum and Acanthoscurria geniculata. For the scorpion Mesobuthus martensii the numbers are 79 and 93 respectively. The very small genome of the house dust mite, Dermatophagoides farinae, on the other hand contains a much smaller number of such genes. A few new putative Arthropod neuropeptide genes were discovered. Thus, both spiders and the scorpion have an achatin gene and in spiders there are two different genes encoding myosuppressin-like peptides while spiders also have two genes encoding novel LGamides. Another finding is the presence of trissin in spiders and scorpions, while neuropeptide genes that seem to be orthologs of Lottia LFRYamide and Platynereis CCRFamide were also found. Such genes were also found in various insect species, but seem to be lacking from the Holometabola. The Chelicerate neuropeptide and neuropeptide GPCR genes often have paralogs. As the large majority of these are probably not due to local gene duplications, is plausible that they reflect the effects of one or more ancient whole genome duplications.

The evolution of insect germline specification strategies

The establishment of the germline is essential for sexually reproducing organisms. In animals, there are two major strategies to specify the germline: maternal provision and zygotic induction. The molecular basis of the maternal provision mode has been well characterized in several model organisms (fly, frog, fish, and nematode), while that of the zygotic induction mode has mainly been studied in mammalian models such as the mouse. Shifts in germline determination modes occur unexpectedly frequently and many such shifts have occurred several times among insects. Given their general tractability and rapidly increasing genomic and genetic tools applicable to many species, the insects present a uniquely powerful model system for understanding major transitions in reproductive strategies, and developmental processes in general.

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.

Evaluation of the impact of RNA preservation methods of spiders for de novo transcriptome assembly

With advances in high-throughput sequencing technologies, de novo transcriptome sequencing and assembly has become a cost-effective method to obtain comprehensive genetic information of a species of interest, especially in nonmodel species with large genomes such as spiders. However, high-quality RNA is essential for successful sequencing, and sample preservation conditions require careful consideration for the effective storage of field-collected samples. To this end, we report a streamlined feasibility study of various storage conditions and their effects on de novo transcriptome assembly results. The storage parameters considered include temperatures ranging from room temperature to -80°C; preservatives, including ethanol, RNAlater, TRIzol and RNAlater-ICE; and sample submersion states. As a result, intact RNA was extracted and assembly was successful when samples were preserved at low temperatures regardless of the type of preservative used. The assemblies as well as the gene expression profiles were shown to be robust to RNA degradation, when 30 million 150-bp paired-end reads are obtained. The parameters for sample storage, RNA extraction, library preparation, sequencing and in silico assembly considered in this work provide a guideline for the study of field-collected samples of spiders.

Neofunctionalization of a Duplicate dachshund Gene Underlies the Evolution of a Novel Leg Segment in Arachnids

The acquisition of a novel function, or neofunctionalization, protects duplicated genes from redundancy and subsequent loss, and is a major force that drives adaptive evolution. Neofunctionalization has been inferred for many duplicated genes based on differences in regulation between the parental gene and its duplicate. However, only few studies actually link the new function of a duplicated gene to a novel morphological or physiological character of the organism. Here we show that the duplication of dachshund (dac) in arachnids (spiders and allies) is linked with the evolution of a novel leg segment, the patella. We have studied dac genes in two distantly related spider species, the entelegyne spider Parasteatoda tepidariorum and the haplogyne spider Pholcus phalangioides. Both species possess two paralogous dac genes that duplicated before the split between entelegyne and haplogyne spiders. In contrast to the evolutionarily highly conserved dac1, its duplicate dac2 is strongly expressed in the patella leg segment during embryogenesis in both species. Using parental RNA interference in P. tepidariorum we show that dac2 is required for the development of the patella segment. If dac2 function is impaired, then the patella is fused with the tibia into a single leg segment. Thus, removing the function of dac2 experimentally reverts P. tepidariorum leg morphology into a stage before the duplication of dac and the evolution of the patella segment. Our results indicate that the origin of the patella is the result of the duplication and subsequent neofunctionalization of dac in the arachnid lineage.

Comparative analysis of tissue-specific transcriptomes in the funnel-web spider Macrothele calpeiana (Araneae, Hexathelidae)

The funnel-web spider Macrothele calpeiana is a charismatic Mygalomorph with a great interest in basic, applied and translational research. Nevertheless, current scarcity of genomic and transcriptomic data of this species clearly limits the research in this non-model organism. To overcome this limitation, we launched the first tissue-specific enriched RNA-seq analysis in this species using a subtractive hybridization approach, with two main objectives, to characterize the specific transcriptome of the putative chemosensory appendages (palps and first pair of legs), and to provide a new set of DNA markers for further phylogenetic studies. We have characterized the set of transcripts specifically expressed in putative chemosensory tissues of this species, much of them showing features shared by chemosensory system genes. Among specific candidates, we have identified some members of the iGluR and NPC2 families. Moreover, we have demonstrated the utility of these newly generated data as molecular markers by inferring the phylogenetic position M. calpeina in the phylogenetic tree of Mygalomorphs. Our results provide novel resources for researchers interested in spider molecular biology and systematics, which can help to expand our knowledge on the evolutionary processes underlying fundamental biological questions, as species invasion or biodiversity origin and maintenance.

Molecular characterization and embryonic origin of the eyes in the common house spider Parasteatoda tepidariorum

Background

Two visual systems are present in most arthropod groups: median and lateral eyes. Most of our current knowledge about the developmental and molecular mechanisms involved in eye formation in arthropods comes from research in the model system Drosophila melanogaster. Here, a core set of retinal determination genes, namely, sine-oculis (so), eyes absent (eya), dachshund (dac), and the two pax6 orthologues eyeless (ey) and twin of eyeless (toy) govern early retinal development. By contrast, not much is known about the development of the up-to-eight eyes present in spiders. Therefore, we analyzed the embryonic expression of core retinal determination genes in the common house spider Parasteatoda tepidariorum.

Results

We show that the anlagen of the median and lateral eyes in P. tepidariorum originate from different regions of the non-neurogenic ectoderm in the embryonic head. The median eyes are specified as two individual anlagen in an anterior median position in the developing head and subsequently move to their final position following extensive morphogenetic movements of the non-neurogenic ectoderm. The lateral eyes develop from a more lateral position. Intriguingly, they are specified as a unique field of cells that splits into the three individual lateral eyes during late embryonic development. Using gene expression analyses, we identified a unique combination of determination gene expression in the anlagen of the lateral and median eyes, respectively.

Conclusions

This study of retinal determination genes in the common house spider P. tepidariorum represents the first comprehensive analysis of the well-known retinal determination genes in arthropods outside insects. The development of the individual lateral eyes via the subdivision of one single eye primordium might be the vestige of a larger composite eye anlage, and thus supports the notion that the composite eye is the plesiomorphic state of the lateral eyes in arthropods. The molecular distinction of the two visual systems is similar to the one described for compound eyes and ocelli in Drosophila, suggesting that a unique core determination network for median and lateral eyes, respectively, might have been in place already in the last common ancestor of spiders and insects.

Electronic supplementary material

The online version of this article (doi:10.1186/s13227-015-0011-9) contains supplementary material, which is available to authorized users.