Expression of arthropod distal limb-patterning genes in the onychophoran Euperipatoides kanangrensis

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

BACKGROUND

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

RESULTS

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

CONCLUSIONS

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

Molecular evidence for a single origin of ultrafiltration-based excretory organs

Excretion is an essential physiological process, carried out by all living organisms, regardless of their size or complexity.^1-3 Both protostomes (e.g., flies and flatworms) and deuterostomes (e.g., humans and sea urchins) possess specialized excretory organs serving that purpose. Those organs exhibit an astonishing diversity, ranging from units composed of just few distinct cells (e.g., protonephridia) to complex structures, built by millions of cells of multiple types with divergent morphology and function (e.g., vertebrate kidneys).^4,5 Although some molecular similarities between the development of kidneys of vertebrates and the regeneration of the protonephridia of flatworms have been reported,^6,7 the molecular underpinnings of the development of excretory organs have never been systematically studied in a comparative context.⁴ Here, we show that a set of transcription factors (eya, six1/2, pou3, sall, lhx1/5, and osr) and structural proteins (nephrin, kirre, and zo1) is expressed in the excretory organs of a phoronid, brachiopod, annelid, onychophoran, priapulid, and hemichordate that represent major protostome lineages and non-vertebrate deuterostomes. We demonstrate that the molecular similarity observed in the vertebrate kidney and flatworm protonephridia^6,7 is also seen in the developing excretory organs of those animals. Our results show that all types of ultrafiltration-based excretory organs are patterned by a conserved set of developmental genes, an observation that supports their homology. We propose that the last common ancestor of protostomes and deuterostomes already possessed an ultrafiltration-based organ that later gave rise to the vast diversity of extant excretory organs, including both proto- and metanephridia.

Loss of intermediate regions of perpendicular body axes contributed to miniaturization of tardigrades

Tardigrades have a miniaturized body plan. Miniaturization in tardigrades is associated with the loss of several organ systems and an intermediate region of their anteroposterior (AP) axis. However, how miniaturization has affected tardigrade legs is unclear. In arthropods and in onychophorans, the leg gap genes are expressed in regionalized proximodistal (PD) patterns in the legs. Functional studies indicate that these genes regulate growth in their respective expression domains and establish PD identities, partly through mutually antagonistic regulatory interactions. Here, we investigated the expression patterns of tardigrade orthologs of the leg gap genes. Rather than being restricted to a proximal leg region, as in arthropods and onychophorans, we detected coexpression of orthologues of homothorax and extradenticle broadly across the legs of the first three trunk segments in the tardigrade Hypsibius exemplaris. We could not identify a dachshund orthologue in tardigrade genomes, a gene that is expressed in an intermediate region of developing legs in arthropods and onychophorans, suggesting that this gene was lost in the tardigrade lineage. We detected Distal-less expression broadly across all developing leg buds in H. exemplaris embryos, unlike in arthropods and onychophorans, in which it exhibits a distally restricted expression domain. The broad expression patterns of the remaining leg gap genes in H. exemplaris legs may reflect the loss of dachshund and the accompanying loss of an intermediate region of the legs in the tardigrade lineage. We propose that the loss of intermediate regions of both the AP and PD body axes contributed to miniaturization of Tardigrada.

Expression of the zinc finger transcription factor Sp6-9 in the velvet worm Euperipatoides kanangrensis suggests a conserved role in appendage development in Panarthropoda

The Sp-family genes encode important transcription factors in animal development. Here we investigate the embryonic expression patterns of the complete set of Sp-genes in the velvet worm Euperipatoides kanangrensis (Onychophora), with a special focus on the Sp6-9 ortholog. In arthropods, Sp6-9, the ortholog of the Drosophila melanogaster D-Sp1 gene plays a conserved role in appendage development. Our data show that the expression of Sp6-9 during the development of the velvet worm is conserved, suggesting that the key function of the Sp6-9 gene dates back to at least the last common ancestor of arthropods and onychophorans and thus likely the last common ancestor of Panarthropoda.

Expression of NK genes that are not part of the NK cluster in the onychophoran Euperipatoides rowelli (Peripatopsidae)

Background

NK genes are a group of homeobox transcription factors that are involved in various molecular pathways across bilaterians. They are typically divided into two subgroups, the NK cluster (NKC) and NK-linked genes (NKL). While the NKC genes have been studied in various bilaterians, corresponding data of many NKL genes are missing to date. To further investigate the ancestral roles of NK family genes, we analyzed the expression patterns of NKL genes in the onychophoran Euperipatoides rowelli.

Results

The NKL gene complement of E. rowelli comprises eight genes, including BarH, Bari, Emx, Hhex, Nedx, NK2.1, vax and NK2.2, of which only NK2.2 was studied previously. Our data for the remaining seven NKL genes revealed expression in different structures associated with the developing nervous system in embryos of E. rowelli. While NK2.1 and vax are expressed in distinct medial regions of the developing protocerebrum early in development, BarH, Bari, Emx, Hhex and Nedx are expressed in late developmental stages, after all major structures of the nervous system have been established. Furthermore, BarH and Nedx are expressed in distinct mesodermal domains in the developing limbs.

Conclusions

Comparison of our expression data to those of other bilaterians revealed similar patterns of NK2.1, vax, BarH and Emx in various aspects of neural development, such as the formation of anterior neurosecretory cells mediated by a conserved molecular mechanism including NK2.1 and vax, and the development of the central and peripheral nervous system involving BarH and Emx. A conserved role in neural development has also been reported from NK2.2, suggesting that the NKL genes might have been primarily involved in neural development in the last common ancestor of bilaterians or at least nephrozoans (all bilaterians excluding xenacoelomorphs). The lack of comparative data for many of the remaining NKL genes, including Bari, Hhex and Nedx currently hampers further evolutionary conclusions. Hence, future studies should focus on the expression of these genes in other bilaterians, which would provide a basis for comparative studies and might help to better understand the role of NK genes in the diversification of bilaterians.

Electronic supplementary material

The online version of this article (10.1186/s12861-019-0185-9) contains supplementary material, which is available to authorized users.

Expression of NK cluster genes in the onychophoran Euperipatoides rowelli: implications for the evolution of NK family genes in nephrozoans

Background

Understanding the evolution and development of morphological traits of the last common bilaterian ancestor is a major goal of the evo-devo discipline. The reconstruction of this "urbilaterian" is mainly based on comparative studies of common molecular patterning mechanisms in recent model organisms. The NK homeobox genes are key players in many of these molecular pathways, including processes regulating mesoderm, heart and neural development. Shared features seen in the expression patterns of NK genes have been used to determine the ancestral bilaterian characters. However, the commonly used model organisms provide only a limited view on the evolution of these molecular pathways. To further investigate the ancestral roles of NK cluster genes, we analyzed their expression patterns in the onychophoran Euperipatoides rowelli.

Results

We identified nine transcripts of NK cluster genes in E. rowelli, including single copies of NK1, NK3, NK4, NK5, Msx, Lbx and Tlx, and two copies of NK6. All of these genes except for NK6.1 and NK6.2 are expressed in different mesodermal organs and tissues in embryos of E. rowelli, including the anlagen of somatic musculature and the heart. Furthermore, we found distinct expression patterns of NK3, NK5, NK6, Lbx and Msx in the developing nervous system. The same holds true for the NKL gene NK2.2, which does not belong to the NK cluster but is a related gene playing a role in neural patterning. Surprisingly, NK1, Msx and Lbx are additionally expressed in a segment polarity-like pattern early in development-a feature that has been otherwise reported only from annelids.

Conclusion

Our results indicate that the NK cluster genes were involved in mesoderm and neural development in the last common ancestor of bilaterians or at least nephrozoans (i.e., bilaterians to the exclusion of xenacoelomorphs). By comparing our data from an onychophoran to those from other bilaterians, we critically review the hypothesis of a complex "urbilaterian" with a segmented body, a pulsatile organ or heart, and a condensed mediolaterally patterned nerve cord.

Gene Regulatory Networks, Homology, and the Early Panarthropod Fossil Record

The arthropod body plan is widely believed to have derived from an ancestral form resembling Cambrian-aged fossil lobopodians, and interpretations of morphological and molecular data have long favored this hypothesis. It is possible, however, that appendages and other morphologies observed in extinct and living panarthropods evolved independently. The key to distinguishing between morphological homology and homoplasy lies in the study of developmental gene regulatory networks (GRNs), and specifically, in determining the unique genetic circuits that construct characters. In this study, I discuss character identity and panarthropod appendage evolution within a developmental GRN framework, with a specific focus on potential limb character identity networks ("ChINs"). I summarize recent molecular studies, and argue that current data do not rule out the possibility of independent panarthropod limb evolution. The link between character identity and GRN architecture has broad implications for homology assessment, and this genetic framework offers alternative approaches to fossil character coding, phylogenetic analyses, and future research into the origin of the arthropod body plan.

Expression and function of spineless orthologs correlate with distal deutocerebral appendage morphology across Arthropoda

The deutocerebral (second) head segment is putatively homologous across Arthropoda, in spite of remarkable disparity of form and function of deutocerebral appendages. In Mandibulata this segment bears a pair of sensory antennae, whereas in Chelicerata the same segment bears a pair of feeding appendages called chelicerae. Part of the evidence for the homology of deutocerebral appendages is the conserved function of homothorax (hth), which has been shown to specify antennal or cheliceral fate in the absence of Hox signaling, in both mandibulate and chelicerate exemplars. However, the genetic basis for the morphological disparity of antenna and chelicera is not understood. To test whether downstream targets of hth have diverged in a lineage-specific manner, we examined the evolution of the function and expression of spineless (ss), which in two holometabolous insects is known to act as a hth target and distal antennal determinant. Toward expanding phylogenetic representation of gene expression data, here we show that strong expression of ss is observed in developing antennae of a hemimetabolous insect, a centipede, and an amphipod crustacean. By contrast, ss orthologs are not expressed throughout the cheliceral limb buds of spiders or harvestmen during developmental stages when appendage fate is specified. RNA interference-mediated knockdown of ss in Oncopeltus fasciatus, which bears a simple plesiomorphic antenna, resulted in homeotic distal antenna-to-leg transformation, comparable to data from holometabolous insect counterparts. Knockdown of hth in Oncopeltus fasciatus abrogated ss expression, suggesting conservation of upstream regulation. These data suggest that ss may be a flagellar (distal antennal) determinant more broadly, and that this function was acquired at the base of Mandibulata.

Gene expression reveals evidence for EGFR-dependent proximal-distal limb patterning in a myriapod

Evolution of segmented limbs is one of the key innovations of Arthropoda, allowing development of functionally specific specialized head and trunk appendages, a major factor behind their unmatched evolutionary success. Proximodistal limb patterning is controlled by two regulatory networks in the vinegar fly Drosophila melanogaster, and other insects. The first is represented by the function of the morphogens Wingless (Wg) and Decapentaplegic (Dpp); the second by the EGFR-signaling cascade. While the role of Wg and Dpp has been studied in a wide range of arthropods representing all main branches, that is, Pancrustacea (= Hexapoda + Crustacea), Myriapoda and Chelicerata, investigation of the potential role of EGFR-signaling is restricted to insects (Hexapoda). Gene expression analysis of Egfr, its potential ligands, and putative downstream factors in the pill millipede Glomeris marginata (Myriapoda: Diplopoda), reveals that-in at least mandibulate arthropods-EGFR-signaling is likely a conserved regulatory mechanism in proximodistal limb patterning.