Unusual gene order and organization of the sea urchin hox cluster

Evolutionary plasticity in nematode Hox gene complements and genomic loci arrangement

Hox genes are central to metazoan body plan formation, patterning and evolution, playing a critical role in cell fate decisions early in embryonic development in invertebrates and vertebrates. While the archetypical Hox gene cluster consists of members of nine ortholog groups (HOX1-HOX9), arrayed in close linkage in the order in which they have their anterior-posterior patterning effects, nematode Hox gene sets do not fit this model. The Caenorhabditis elegans Hox gene set is not clustered and contains only six Hox genes from four of the ancestral groups. The pattern observed in C. elegans is not typical of the phylum, and variation in orthologue set presence and absence and in genomic organisation has been reported. Recent advances in genome sequencing have resulted in the availability of many novel genome assemblies in Nematoda, especially from taxonomic groups that had not been analysed previously. Here, we explored Hox gene complements in high-quality genomes of 80 species from all major clades of Nematoda to understand the evolution of this key set of body pattern genes and especially to probe the origins of the "dispersed" cluster observed in C. elegans. We also included the recently available high-quality genomes of some Nematomorpha as an outgroup. We find that nematodes can have Hox genes from up to six orthology groups. While nematode Hox "clusters" are often interrupted by unrelated genes we identify species in which the cluster is intact and not dispersed.

The brittle star genome illuminates the genetic basis of animal appendage regeneration

Species within nearly all extant animal lineages are capable of regenerating body parts. However, it remains unclear whether the gene expression programme controlling regeneration is evolutionarily conserved. Brittle stars are a species-rich class of echinoderms with outstanding regenerative abilities, but investigations into the genetic bases of regeneration in this group have been hindered by the limited genomic resources. Here we report a chromosome-scale genome assembly for the brittle star Amphiura filiformis. We show that the brittle star genome is the most rearranged among echinoderms sequenced so far, featuring a reorganized Hox cluster reminiscent of the rearrangements observed in sea urchins. In addition, we performed an extensive profiling of gene expression during brittle star adult arm regeneration and identified sequential waves of gene expression governing wound healing, proliferation and differentiation. We conducted comparative transcriptomic analyses with other invertebrate and vertebrate models for appendage regeneration and uncovered hundreds of genes with conserved expression dynamics, particularly during the proliferative phase of regeneration. Our findings emphasize the crucial importance of echinoderms to detect long-range expression conservation between vertebrates and classical invertebrate regeneration model systems.

Sea cucumbers: an emerging system in evo-devo

A challenge for evolutionary developmental (evo-devo) biology is to expand the breadth of research organisms used to investigate how animal diversity has evolved through changes in embryonic development. New experimental systems should couple a relevant phylogenetic position with available molecular tools and genomic resources. As a phylum of the sister group to chordates, echinoderms extensively contributed to our knowledge of embryonic patterning, organ development and cell-type evolution. Echinoderms display a variety of larval forms with diverse shapes, making them a suitable group to compare the evolution of embryonic developmental strategies. However, because of the laboratory accessibility and the already available techniques, most studies focus on sea urchins and sea stars mainly. As a comparative approach, the field would benefit from including information on other members of this group, like the sea cucumbers (holothuroids), for which little is known on the molecular basis of their development. Here, we review the spawning and culture methods, the available morphological and molecular information, and the current state of genomic and transcriptomic resources on sea cucumbers. With the goal of making this system accessible to the broader community, we discuss how sea cucumber embryos and larvae can be a powerful system to address the open questions in evo-devo, including understanding the origins of bilaterian structures.

Comparative Analysis of Bivalve and Sea Urchin Genetics and Development: Investigating the Dichotomy in Bilateria

This comprehensive review presents a comparative analysis of early embryogenesis in Protostomia and Deuterostomia, the first of which exhibit a mosaic pattern of development, where cells are fated deterministically, while Deuterostomia display a regulatory pattern of development, where the fate of cells is indeterminate. Despite these fundamental differences, there are common transcriptional mechanisms that underline their evolutionary linkages, particularly in the field of functional genomics. By elucidating both conserved and unique regulatory strategies, this review provides essential insights into the comparative embryology and developmental dynamics of these groups. The objective of this review is to clarify the shared and distinctive characteristics of transcriptional regulatory mechanisms. This will contribute to the extensive areas of functional genomics, evolutionary biology and developmental biology, and possibly lay the foundation for future research and discussion on this seminal topic.

Spiralian genomics and the evolution of animal genome architecture

Recent developments in sequencing technologies have greatly improved our knowledge of phylogenetic relationships and genomic architectures throughout the tree of life. Spiralia, a diverse clade within Protostomia, is essential for understanding the evolutionary history of parasitism, gene conversion, nervous systems and animal body plans. In this review, we focus on the current hypotheses of spiralian phylogeny and investigate the impact of long-read sequencing on the quality of genome assemblies. We examine chromosome-level assemblies to highlight key genomic features that have driven spiralian evolution, including karyotype, synteny and the Hox gene organization. In addition, we show how chromosome rearrangement has influenced spiralian genomic structures. Although spiralian genomes have undergone substantial changes, they exhibit both conserved and lineage-specific features. We recommend increasing sequencing efforts and expanding functional genomics research to deepen insights into spiralian biology.

Analysis of the P. lividus sea urchin genome highlights contrasting trends of genomic and regulatory evolution in deuterostomes

Sea urchins are emblematic models in developmental biology and display several characteristics that set them apart from other deuterostomes. To uncover the genomic cues that may underlie these specificities, we generated a chromosome-scale genome assembly for the sea urchin Paracentrotus lividus and an extensive gene expression and epigenetic profiles of its embryonic development. We found that, unlike vertebrates, sea urchins retained ancestral chromosomal linkages but underwent very fast intrachromosomal gene order mixing. We identified a burst of gene duplication in the echinoid lineage and showed that some of these expanded genes have been recruited in novel structures (water vascular system, Aristotle's lantern, and skeletogenic micromere lineage). Finally, we identified gene-regulatory modules conserved between sea urchins and chordates. Our results suggest that gene-regulatory networks controlling development can be conserved despite extensive gene order rearrangement.

A single-cell RNA-seq analysis of early larval cell-types of the starfish, Patiria pectinifera: Insights into evolution of the chordate body plan

Ambulacrarians (echinoderms and hemichordates) are a sister group to chordates; thus, their larval cell-types may provide clues about evolution of chordate body plans. Although most genic information accumulated to date pertains to sea urchin embryogenesis, starfish embryogenesis represents a more ancestral mode than that of sea urchins. We performed single-cell RNA-seq analysis of cell-types from gastrulae and bipinnarial larvae of the starfish, Patiria pectinifera, and categorized them into 22 clusters, each of which is composed of cells with specific, shared profiles of development-relevant gene expression. Oral and aboral ectoderm, apical plate, hindgut or archenteron, midgut or intestine, pharynx, endomesoderm, stomodeum, and mesenchyme of the gastrulae, and neurons, ciliary bands, enterocoel and muscle of larvae were characterized by expression profiles of at least two relevant transcription factor genes and signaling molecular genes. Expression of Hox2, Hox7, Hox9/10, and Hox11/13b was detected in cells of clusters that form the larval enterocoel. By comparing homologous gene expression profiles in chordate embryos, we discuss and propose how the chordate body plan evolved from a deuterostome ancestor, from which the echinoderm body plan also evolved.

Worms and gills, plates and spines: the evolutionary origins and incredible disparity of deuterostomes revealed by fossils, genes, and development

Deuterostomes are the major division of animal life which includes sea stars, acorn worms, and humans, among a wide variety of ecologically and morphologically disparate taxa. However, their early evolution is poorly understood, due in part to their disparity, which makes identifying commonalities difficult, as well as their relatively poor early fossil record. Here, we review the available morphological, palaeontological, developmental, and molecular data to establish a framework for exploring the origins of this important and enigmatic group. Recent fossil discoveries strongly support a vermiform ancestor to the group Hemichordata, and a fusiform active swimmer as ancestor to Chordata. The diverse and anatomically bewildering variety of forms among the early echinoderms show evidence of both bilateral and radial symmetry. We consider four characteristics most critical for understanding the form and function of the last common ancestor to Deuterostomia: Hox gene expression patterns, larval morphology, the capacity for biomineralization, and the morphology of the pharyngeal region. We posit a deuterostome last common ancestor with a similar antero-posterior gene regulatory system to that found in modern acorn worms and cephalochordates, a simple planktonic larval form, which was later elaborated in the ambulacrarian lineage, the ability to secrete calcium minerals in a limited fashion, and a pharyngeal respiratory region composed of simple pores. This animal was likely to be motile in adult form, as opposed to the sessile origins that have been historically suggested. Recent debates regarding deuterostome monophyly as well as the wide array of deuterostome-affiliated problematica further suggest the possibility that those features were not only present in the last common ancestor of Deuterostomia, but potentially in the ur-bilaterian. The morphology and development of the early deuterostomes, therefore, underpin some of the most significant questions in the study of metazoan evolution.

A Chromosome-level Genome Assembly of the Highly Heterozygous Sea Urchin Echinometra sp. EZ Reveals Adaptation in the Regulatory Regions of Stress Response Genes

Echinometra is the most widespread genus of sea urchin and has been the focus of a wide range of studies in ecology, speciation, and reproduction. However, available genetic data for this genus are generally limited to a few select loci. Here, we present a chromosome-level genome assembly based on 10x Genomics, PacBio, and Hi-C sequencing for Echinometra sp. EZ from the Persian/Arabian Gulf. The genome is assembled into 210 scaffolds totaling 817.8 Mb with an N50 of 39.5 Mb. From this assembly, we determined that the E. sp. EZ genome consists of 2n = 42 chromosomes. BUSCO analysis showed that 95.3% of BUSCO genes were complete. Ab initio and transcript-informed gene modeling and annotation identified 29,405 genes, including a conserved Hox cluster. E. sp. EZ can be found in high-temperature and high-salinity environments, and we therefore compared E. sp. EZ gene families and transcription factors associated with environmental stress response ("defensome") with other echinoid species with similar high-quality genomic resources. While the number of defensome genes was broadly similar for all species, we identified strong signatures of positive selection in E. sp. EZ noncoding elements near genes involved in environmental response pathways as well as losses of transcription factors important for environmental response. These data provide key insights into the biology of E. sp. EZ as well as the diversification of Echinometra more widely and will serve as a useful tool for the community to explore questions in this taxonomic group and beyond.

Echiuran Hox genes provide new insights into the correspondence between Hox subcluster organization and collinearity pattern

In many bilaterians, Hox genes are generally clustered along the chromosomes and expressed in spatial and temporal order. In vertebrates, the expression of Hox genes follows a whole-cluster spatio-temporal collinearity (WSTC) pattern, whereas in some invertebrates the expression of Hox genes exhibits a subcluster-level spatio-temporal collinearity pattern. In bilaterians, the diversity of collinearity patterns and the cause of collinearity differences in Hox gene expression remain poorly understood. Here, we investigate genomic organization and expression pattern of Hox genes in the echiuran worm Urechis unicinctus (Annelida, Echiura). Urechis unicinctus has a split cluster with four subclusters divided by non-Hox genes: first subcluster (Hox1 and Hox2), second subcluster (Hox3), third subcluster (Hox4, Hox5, Lox5, Antp and Lox4), fourth subcluster (Lox2 and Post2). The expression of U. unicinctus Hox genes shows a subcluster-based whole-cluster spatio-temporal collinearity (S-WSTC) pattern: the anterior-most genes in each subcluster are activated in a spatially and temporally colinear manner (reminiscent of WSTC), with the subsequent genes in each subcluster then being very similar to their respective anterior-most subcluster gene. Combining genomic organization and expression profiles of Hox genes in different invertebrate lineages, we propose that the spatio-temporal collinearity of invertebrate Hox is subcluster-based.

Human follicular mites: Ectoparasites becoming symbionts

Most humans carry mites in the hair follicles of their skin for their entire lives. Follicular mites are the only metazoans tha continuously live on humans. We propose that Demodex folliculorum (Acari) represents a transitional stage from a host-injuring obligate parasite to an obligate symbiont. Here, we describe the profound impact of this transition on the genome and physiology of the mite. Genome sequencing revealed that the permanent host association of D. folliculorum led to an extensive genome reduction through relaxed selection and genetic drift, resulting in the smallest number of protein-coding genes yet identified among panarthropods. Confocal microscopy revealed that this gene loss coincided with an extreme reduction in the number of cells. Single uninucleate muscle cells are sufficient to operate each of the three segments that form each walking leg. While it has been assumed that the reduction of the cell number in parasites starts early in development, we identified a greater total number of cells in the last developmental stage (nymph) than in the terminal adult stage, suggesting that reduction starts at the adult or ultimate stage of development. This is the first evolutionary step in an arthropod species adopting a reductive, parasitic or endosymbiotic lifestyle. Somatic nuclei show underreplication at the diploid stage. Novel eye structures or photoreceptors as well as a unique human host melatonin-guided day/night rhythm are proposed for the first time. The loss of DNA repair genes coupled with extreme endogamy might have set this mite species on an evolutionary dead-end trajectory.

Chromosomal-Level Genome Assembly of the Sea Urchin Lytechinus variegatus Substantially Improves Functional Genomic Analyses

Lytechinus variegatus is a camarodont sea urchin found widely throughout the western Atlantic Ocean in a variety of shallow-water marine habitats. Its distribution, abundance, and amenability to developmental perturbation make it a popular model for ecologists and developmental biologists. Here, we present a chromosomal-level genome assembly of L. variegatus generated from a combination of PacBio long reads, 10× Genomics sequencing, and HiC chromatin interaction sequencing. We show L. variegatus has 19 chromosomes with an assembly size of 870.4 Mb. The contiguity and completeness of this assembly are reflected by a scaffold length N50 of 45.5 Mb and BUSCO completeness score of 95.5%. Ab initio and transcript-informed gene modeling and annotation identified 27,232 genes with an average gene length of 12.6 kb, comprising an estimated 39.5% of the genome. Repetitive regions, on the other hand, make up 45.4% of the genome. Physical mapping of well-studied developmental genes onto each chromosome reveals nonrandom spatial distribution of distinct genes and gene families, which provides insight into how certain gene families may have evolved and are transcriptionally regulated in this species. Lastly, aligning RNA-seq and ATAC-seq data onto this assembly demonstrates the value of highly contiguous, complete genome assemblies for functional genomics analyses that is unattainable with fragmented, incomplete assemblies. This genome will be of great value to the scientific community as a resource for genome evolution, developmental, and ecological studies of this species and the Echinodermata.

Physical Laws Shape Up HOX Gene Collinearity

Hox gene collinearity (HGC) is a multi-scalar property of many animal phyla particularly important in embryogenesis. It relates entities and events occurring in Hox clusters inside the chromosome DNA and in embryonic tissues. These two entities differ in linear size by more than four orders of magnitude. HGC is observed as spatial collinearity (SC), where the Hox genes are located in the order (Hox1, Hox2, Hox3 …) along the 3′ to 5′ direction of DNA in the genome and a corresponding sequence of ontogenetic units (E1, E2, E3, …) located along the Anterior—Posterior axis of the embryo. Expression of Hox1 occurs in E1, Hox2 in E2, Hox3 in E3, etc. Besides SC, a temporal collinearity (TC) has been also observed in many vertebrates. According to TC, first Hox1 is expressed in E1; later, Hox2 is expressed in E2, followed by Hox3 in E3, etc. Lately, doubt has been raised about whether TC really exists. A biophysical model (BM) was formulated and tested during the last 20 years. According to BM, physical forces are created which pull the Hox genes one after the other, driving them to a transcription factory domain where they are transcribed. The existing experimental data support this BM description. Symmetry is a physical–mathematical property of matter that was explored in depth by Noether who formulated a ground-breaking theory (NT) that applies to all sizes of matter. NT may be applied to biology in order to explain the origin of HGC in animals developing not only along the A/P axis, but also to animals with circular symmetry.

Transcriptomic analysis of sea star development through metamorphosis to the highly derived pentameral body plan with a focus on neural transcription factors

The Echinodermata is characterized by a secondarily evolved pentameral body plan. While the evolutionary origin of this body plan has been the subject of debate, the molecular mechanisms underlying its development are poorly understood. We assembled a de novo developmental transcriptome from the embryo through metamorphosis in the sea star Parvulastra exigua. We use the asteroid model as it represents the basal-type echinoderm body architecture. Global variation in gene expression distinguished the gastrula profile and showed that metamorphic and juvenile stages were more similar to each other than to the pre-metamorphic stages, pointing to the marked changes that occur during metamorphosis. Differential expression and gene ontology (GO) analyses revealed dynamic changes in gene expression throughout development and the transition to pentamery. Many GO terms enriched during late metamorphosis were related to neurogenesis and signalling. Neural transcription factor genes exhibited clusters with distinct expression patterns. A suite of these genes was up-regulated during metamorphosis (e.g. Pax6, Eya, Hey, NeuroD, FoxD, Mbx, and Otp). In situ hybridization showed expression of neural genes in the CNS and sensory structures. Our results provide a foundation to understand the metamorphic transition in echinoderms and the genes involved in development and evolution of pentamery.

Multiple origins of feeding head larvae by the Early Cambrian

In many animals the head develops early, most of the body axis later. A larva composed mostly of the developing front end therefore can attain mobility and feeding earlier in development. Fossils, functional morphology, and inferred homologies indicate that feeding head larvae existed by the Early Cambrian in members of three major clades of animals: ecdysozoans, lophotrochozoans, and deuterostomes. Some of these early larval feeding mechanisms were also those of juveniles and adults (the lophophore of brachiopod larvae and possibly the ciliary band of the dipleurula of hemichordates and echinoderms); some were derived from structures that previously had other functions (appendages of the nauplius). Trochophores that swim with a preoral band of cilia, the prototroch, originated before divergence of annelids and molluscs, but evidence of larval growth and thus a prototrochal role in feeding is lacking for molluscs until the Ordovician. Feeding larvae that definitely originated much later, as in insects, teleost fish, and amphibians, develop all or nearly all of what will become the adult body axis before they begin feeding. On present evidence, head larvae, including feeding head larvae, evolved multiple times early in the evolution of bilaterian animals and never since.

Genomic insights of body plan transitions from bilateral to pentameral symmetry in Echinoderms

Echinoderms are an exceptional group of bilaterians that develop pentameral adult symmetry from a bilaterally symmetric larva. However, the genetic basis in evolution and development of this unique transformation remains to be clarified. Here we report newly sequenced genomes, developmental transcriptomes, and proteomes of diverse echinoderms including the green sea urchin (L. variegatus), a sea cucumber (A. japonicus), and with particular emphasis on a sister group of the earliest-diverged echinoderms, the feather star (A. japonica). We learned that the last common ancestor of echinoderms retained a well-organized Hox cluster reminiscent of the hemichordate, and had gene sets involved in endoskeleton development. Further, unlike in other animal groups, the most conserved developmental stages were not at the body plan establishing phase, and genes normally involved in bilaterality appear to function in pentameric axis development. These results enhance our understanding of the divergence of protostomes and deuterostomes almost 500 Mya.

Li et al. investigate the evolution and genetic basis of the adult pentameral body plan in echinoderms using genomic, transcriptomic, and proteomic data. They determine that the last common ancestor of echinoderms contained an organized Hox cluster and endoskeleton genes, and suggest that cooption of bilateral development genes was involved in evolution of the pentameric body plan.

Hox Gene Collinearity May Be Related to Noether Theory on Symmetry and Its Linked Conserved Quantity

Hox Gene Collinearity (HGC) is a fundamental property that controls the development of many animal species, including vertebrates. In the Hox gene clusters, the genes are located in a sequential order Hox1, Hox2, Hox3, etc., along the 3’ to 5’ direction of the cluster in the chromosome. During Hox cluster activation, the Hox genes are expressed sequentially in the ontogenetic units D1, D2, D3, etc., along the anterior–posterior axis (A-P) of the early embryo. This collinearity, first observed by E.B. Lewis, is surprising because the spatial collinearity of these structures (Hox clusters and embryos) correlates entities that differ by about four orders of magnitude. Biomolecular mechanisms alone cannot explain such correlations. Long-range physical interactions, such as diffusion or electric attractions, should be involved. A biophysical model (BM) was formulated, which, in alignment with the biomolecular processes, successfully describes the existing vertebrate genetic engineering data. One hundred years ago, Emmy Noether made a fundamental discovery in mathematics and physics. She proved, rigorously, that a physical system obeying a symmetry law (e.g., rotations or self-similarity) is followed by a conserved physical quantity. It is argued here that HGC obeys a ‘primitive’ self-similarity symmetry. In this case, the associated primitive conserved quantity is the irreversibly increasing ‘ratchet’-like Hoxgene ordering where some genes may be missing. The genes of a vertebrate Hox clusterare located along a finite straight line. The same order follows the ontogenetic unitsof the vertebrate embryo. Therefore, HGC is a manifestation of a primitive Noether Theory (NT). NT may be applied to other than the vertebrate case, for instance, to animals with a circular topological symmetry. For example, the observed abnormal Hox gene ordering of the echinoderm Hox clusters may be reproduced by a double-strand break of the circular Hox gene ordering and its subsequent incorporation in the flanking chromosome.

Bacterial artificial chromosomes as recombinant reporter constructs to investigate gene expression and regulation in echinoderms

Genome sequences contain all the necessary information-both coding and regulatory sequences-to construct an organism. The developmental process translates this genomic information into a three-dimensional form. One interpretation of this translation process can be described using gene regulatory network (GRN) models, which are maps of interactions among regulatory gene products in time and space. As high throughput investigations reveal increasing complexity within these GRNs, it becomes apparent that efficient methods are required to test the necessity and sufficiency of regulatory interactions. One of the most complete GRNs for early development has been described in the purple sea urchin, Strongylocentrotus purpuratus. This work has been facilitated by two resources: a well-annotated genome sequence and transgenes generated in bacterial artificial chromosome (BAC) constructs. BAC libraries played a central role in assembling the S. purpuratus genome sequence and continue to serve as platforms for generating reporter constructs for use in expression and regulatory analyses. Optically transparent echinoderm larvae are highly amenable to transgenic approaches and are therefore particularly well suited for experiments that rely on BAC-based reporter transgenes. Here, we discuss the experimental utility of BAC constructs in the context of understanding developmental processes in echinoderm embryos and larvae.

Hox gene expression in postmetamorphic juveniles of the brachiopod Terebratalia transversa

BACKGROUND

Hox genes encode a family of homeodomain containing transcription factors that are clustered together on chromosomes of many Bilateria. Some bilaterian lineages express these genes during embryogenesis in spatial and/or temporal order according to their arrangement in the cluster, a phenomenon referred to as collinearity. Expression of Hox genes is well studied during embryonic and larval development of numerous species; however, relatively few studies focus on the comparison of pre- and postmetamorphic expression of Hox genes in animals with biphasic life cycle. Recently, the expression of Hox genes was described for embryos and larvae of Terebratalia transversa, a rhynchonelliformean brachiopod, which possesses distinct metamorphosis from planktonic larvae to sessile juveniles. During premetamorphic development, T. transversa does not exhibit spatial collinearity and several of its Hox genes are recruited for the morphogenesis of novel structures. In our study, we determined the expression of Hox genes in postmetamorphic juveniles of T. transversa in order to examine metamorphosis-related changes of expression patterns and to test whether Hox genes are expressed in the spatially collinear way in the postmetamorphic juveniles.

RESULTS

Hox genes are expressed in a spatially non-collinear manner in juveniles, generally showing similar patterns as ones observed in competent larvae: genes labial and post1 are expressed in chaetae-related structures, sex combs reduced in the shell-forming epithelium, whereas lox5 and lox4 in dorso-posterior epidermis. After metamorphosis, expression of genes proboscipedia, hox3, deformed and antennapedia becomes restricted to, respectively, shell musculature, prospective hinge rudiments and pedicle musculature and epidermis.

CONCLUSIONS

All developmental stages of T. transversa, including postmetamorphic juveniles, exhibit a spatial non-collinear Hox genes expression with only minor changes observed between pre- and postmetamorphic stages. Our results are concordant with morphological observation that metamorphosis in rhynchonelliformean brachiopods, despite being rapid, is rather gradual. The most drastic changes in Hox gene expression patterns observed during metamorphosis could be explained by the inversion of the mantle lobe, which relocates some of the more posterior larval structures into the anterior edge of the juveniles. Co-option of Hox genes for the morphogenesis of novel structures is even more pronounced in postmetamorphic brachiopods when compared to larvae.

Two more Posterior Hox genes and Hox cluster dispersal in echinoderms

BACKGROUND

Hox genes are key elements in patterning animal development. They are renowned for their, often, clustered organisation in the genome, with supposed mechanistic links between the organisation of the genes and their expression. The widespread distribution and comparable functions of Hox genes across the animals has led to them being a major study system for comparing the molecular bases for construction and divergence of animal morphologies. Echinoderms (including sea urchins, sea stars, sea cucumbers, feather stars and brittle stars) possess one of the most unusual body plans in the animal kingdom with pronounced pentameral symmetry in the adults. Consequently, much interest has focused on their development, evolution and the role of the Hox genes in these processes. In this context, the organisation of echinoderm Hox gene clusters is distinctive. Within the classificatory system of Duboule, echinoderms constitute one of the clearest examples of Disorganized (D) clusters (i.e. intact clusters but with a gene order or orientation rearranged relative to the ancestral state).

RESULTS

Here we describe two Hox genes (Hox11/13d and e) that have been overlooked in most previous work and have not been considered in reconstructions of echinoderm Hox complements and cluster organisation. The two genes are related to Posterior Hox genes and are present in all classes of echinoderm. Importantly, they do not reside in the Hox cluster of any species for which genomic linkage data is available.

CONCLUSION

Incorporating the two neglected Posterior Hox genes into assessments of echinoderm Hox gene complements and organisation shows that these animals in fact have Split (S) Hox clusters rather than simply Disorganized (D) clusters within the Duboule classification scheme. This then has implications for how these genes are likely regulated, with them no longer covered by any potential long-range Hox cluster-wide, or multigenic sub-cluster, regulatory mechanisms.

The phylum Vertebrata: a case for zoological recognition

The group Vertebrata is currently placed as a subphylum in the phylum Chordata, together with two other subphyla, Cephalochordata (lancelets) and Urochordata (ascidians). The past three decades, have seen extraordinary advances in zoological taxonomy and the time is now ripe for reassessing whether the subphylum position is truly appropriate for vertebrates, particularly in light of recent advances in molecular phylogeny, comparative genomics, and evolutionary developmental biology. Four lines of current research are discussed here. First, molecular phylogeny has demonstrated that Deuterostomia comprises Ambulacraria (Echinodermata and Hemichordata) and Chordata (Cephalochordata, Urochordata, and Vertebrata), each clade being recognized as a mutually comparable phylum. Second, comparative genomic studies show that vertebrates alone have experienced two rounds of whole-genome duplication, which makes the composition of their gene family unique. Third, comparative gene-expression profiling of vertebrate embryos favors an hourglass pattern of development, the most conserved stage of which is recognized as a phylotypic period characterized by the establishment of a body plan definitively associated with a phylum. This mid-embryonic conservation is supported robustly in vertebrates, but only weakly in chordates. Fourth, certain complex patterns of body plan formation (especially of the head, pharynx, and somites) are recognized throughout the vertebrates, but not in any other animal groups. For these reasons, we suggest that it is more appropriate to recognize vertebrates as an independent phylum, not as a subphylum of the phylum Chordata.

The sea cucumber genome provides insights into morphological evolution and visceral regeneration

Apart from sharing common ancestry with chordates, sea cucumbers exhibit a unique morphology and exceptional regenerative capacity. Here we present the complete genome sequence of an economically important sea cucumber, A. japonicus, generated using Illumina and PacBio platforms, to achieve an assembly of approximately 805 Mb (contig N50 of 190 Kb and scaffold N50 of 486 Kb), with 30,350 protein-coding genes and high continuity. We used this resource to explore key genetic mechanisms behind the unique biological characters of sea cucumbers. Phylogenetic and comparative genomic analyses revealed the presence of marker genes associated with notochord and gill slits, suggesting that these chordate features were present in ancestral echinoderms. The unique shape and weak mineralization of the sea cucumber adult body were also preliminarily explained by the contraction of biomineralization genes. Genome, transcriptome, and proteome analyses of organ regrowth after induced evisceration provided insight into the molecular underpinnings of visceral regeneration, including a specific tandem-duplicated prostatic secretory protein of 94 amino acids (PSP94)-like gene family and a significantly expanded fibrinogen-related protein (FREP) gene family. This high-quality genome resource will provide a useful framework for future research into biological processes and evolution in deuterostomes, including remarkable regenerative abilities that could have medical applications. Moreover, the multiomics data will be of prime value for commercial sea cucumber breeding programs.

Echinoderms, ubiquitous in the marine environment, are important from evolutionary, ecological, and socioeconomic perspectives. Together with chordates and hemichordates, they form the deuterostome clade, making them a crucial node in the study of chordate ancestry. Within echinoderms, class Holothuroidea is unique; its members (the sea cucumbers) display remarkable regenerative abilities and play key roles as sediment bioturbators and symbiotic hosts, and many are prized in the seafood and pharmaceutical industries. The sea cucumber genome therefore has the potential to significantly contribute to our understanding of important evolutionary and biological processes and help enhance aquaculture programs. Here we present a high-quality genome sequence for the economically important species Apostichopus japonicus. Through comparative analyses, we identified 763 echinoderm-specific gene families enriched in genes encoding membrane proteins, ion channels, and signal transduction proteins. Marker genes associated with the notochord and gill slits were also found, providing valuable insight into the origin of chordates. The reduced number and low expression levels of biomineralization genes reflect the skeletal degeneration seen in sea cucumbers. Importantly, 2 gene families appeared to be expanded in A. japonicus and may play crucial roles in its heightened regenerative potential. Together, findings from the sea cucumber genome provide important and novel insights into echinoderm and deuterostome biology.

Topological Domains, Metagenes, and the Emergence of Pleiotropic Regulations at Hox Loci

Hox gene clusters of jaw vertebrates display a tight genomic organization, which has no equivalent in any other bilateria genomes sequenced thus far. It was previously argued that such a topological consolidation toward a condensed, metagenic structure was due to the accumulation of long-range regulations flanking Hox loci, a phenomenon made possible by the successive genome duplications that occurred at the root of the vertebrate lineage, similar to gene neofunctionalization but applied to a coordinated multigenic system. Here, we propose that the emergence of such large vertebrate regulatory landscapes containing a range of global enhancers was greatly facilitated by the presence of topologically associating domains (TADs). These chromatin domains, mostly constitutive, may have been used as genomic niches where novel regulations could evolve due to both the preexistence of a structural backbone poised to integrate novel regulatory inputs, and a highly adaptive transcriptional readout. We propose a scenario for the coevolution of such TADs and the emergence of pleiotropy at ancestral vertebrate Hox loci.

Coping with living in the soil: the genome of the parthenogenetic springtail Folsomia candida

Background

Folsomia candida is a model in soil biology, belonging to the family of Isotomidae, subclass Collembola. It reproduces parthenogenetically in the presence of Wolbachia, and exhibits remarkable physiological adaptations to stress. To better understand these features and adaptations to life in the soil, we studied its genome in the context of its parthenogenetic lifestyle.

Results

We applied Pacific Bioscience sequencing and assembly to generate a reference genome for F. candida of 221.7 Mbp, comprising only 162 scaffolds. The complete genome of its endosymbiont Wolbachia, was also assembled and turned out to be the largest strain identified so far. Substantial gene family expansions and lineage-specific gene clusters were linked to stress response. A large number of genes (809) were acquired by horizontal gene transfer. A substantial fraction of these genes are involved in lignocellulose degradation. Also, the presence of genes involved in antibiotic biosynthesis was confirmed. Intra-genomic rearrangements of collinear gene clusters were observed, of which 11 were organized as palindromes. The Hox gene cluster of F. candida showed major rearrangements compared to arthropod consensus cluster, resulting in a disorganized cluster.

Conclusions

The expansion of stress response gene families suggests that stress defense was important to facilitate colonization of soils. The large number of HGT genes related to lignocellulose degradation could be beneficial to unlock carbohydrate sources in soil, especially those contained in decaying plant and fungal organic matter. Intra- as well as inter-scaffold duplications of gene clusters may be a consequence of its parthenogenetic lifestyle. This high quality genome will be instrumental for evolutionary biologists investigating deep phylogenetic lineages among arthropods and will provide the basis for a more mechanistic understanding in soil ecology and ecotoxicology.

Electronic supplementary material

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

Rotiferan Hox genes give new insights into the evolution of metazoan bodyplans

The phylum Rotifera consists of minuscule, nonsegmented animals with a unique body plan and an unresolved phylogenetic position. The presence of pharyngeal articulated jaws supports an inclusion in Gnathifera nested in the Spiralia. Comparison of Hox genes, involved in animal body plan patterning, can be used to infer phylogenetic relationships. Here, we report the expression of five Hox genes during embryogenesis of the rotifer Brachionus manjavacas and show how these genes define different functional components of the nervous system and not the usual bilaterian staggered expression along the anteroposterior axis. Sequence analysis revealed that the lox5-parapeptide, a key signature in lophotrochozoan and platyhelminthean Hox6/lox5 genes, is absent and replaced by different signatures in Rotifera and Chaetognatha, and that the MedPost gene, until now unique to Chaetognatha, is also present in rotifers. Collectively, our results support an inclusion of chaetognaths in gnathiferans and Gnathifera as sister group to the remaining spiralians.

Rotifers are microscopic animals with an unusual, nonsegmented body plan consisting of a head, trunk and foot. Here, Fröbius and Funch investigate the role of Hox genes—which are widely used in animal body plan patterning—in rotifer embryogenesis and find non-canonical expression in the nervous system.

Echinoderm development and evolution in the post-genomic era

The highly recognizable animals within the phylum Echinodermata encompass an enormous disparity of adult and larval body plans. The extensive knowledge of sea urchin development has culminated in the description of the exquisitely detailed gene regulatory network (GRN) that governs the specification of various embryonic territories. This information provides a unique opportunity for comparative studies in other echinoderm taxa to understand the evolution and developmental mechanisms underlying body plan change. This review focuses on recent work that has utilized new genomic resources and systems-level experiments to address questions of evolutionary developmental biology. In particular, we synthesize the results of several recent studies from various echinoderm classes that have explored the development and evolution of the larval skeleton, which is a major feature that distinguishes the two predominant larval subtypes within the Phylum. We specifically examine the ways in which GRNs can evolve, either through cis regulatory and/or protein-level changes in transcription factors. We also examine recent work comparing evolution across shorter time scales that occur within and between species of sea urchin, and highlight the kinds of questions that can be addressed by these comparisons. The advent of new genomic and transcriptomic datasets in additional species from all classes of echinoderm will continue to empower the use of these taxa for evolutionary developmental studies.

Hox Gene Collinearity: From A-P Patterning to Radially Symmetric Animals

Hox gene collinearity relates the gene order of the Hox cluster in the chromosome (telomeric to centromeric end) with the serial activation of these genes in the ontogenetic units along the Anterior-Posterior embryonic axis. Although this collinearity property is well respected in bilaterians (e.g. vertebrates), it is violated in other animals. The A-P axis is established in the early embryo of the sea urchin. Subsequently, rotational symmetry is superimposed when the vestibula larva is formed. In analogy to the linear A-P case, it is here hypothesized that the circular topology of the ontogenetic modules is associated to the architectural restructuring of the Hox loci where the two discrete ends of the Hox cluster approach each other so that an almost circular DNA contour is created. In the evolutionary process the circular mode undergoes double strand breaks and the generated cluster ends are attached to the open ends of the flanking chromosome. This event may lead to a novel gene ordering associated with an evolutionary innovation. For example, the loss of Hox4 is followed by the formation of a shorter gene circular arrangement. The opening of this contour at the missing Hox4 location and its connection to the chromosomal flanking ends leads to a new diversification namely the creation of the unusual gene order of the sea urchin Hox cluster.

Low coverage sequencing of three echinoderm genomes: the brittle star Ophionereis fasciata, the sea star Patiriella regularis, and the sea cucumber Australostichopus mollis

BACKGROUND

There are five major extant groups of Echinodermata: Crinoidea (feather stars and sea lillies), Ophiuroidea (brittle stars and basket stars), Asteroidea (sea stars), Echinoidea (sea urchins, sea biscuits, and sand dollars), and Holothuroidea (sea cucumbers). These animals are known for their pentaradial symmetry as adults, unique water vascular system, mutable collagenous tissues, and endoskeletons of high magnesium calcite. To our knowledge, the only echinoderm species with a genome sequence available to date is Strongylocentrotus pupuratus (Echinoidea). The availability of additional echinoderm genome sequences is crucial for understanding the biology of these animals.

FINDINGS

Here we present assembled draft genomes of the brittle star Ophionereis fasciata, the sea star Patiriella regularis, and the sea cucumber Australostichopus mollis from Illumina sequence data with coverages of 12.5x, 22.5x, and 21.4x, respectively.

CONCLUSIONS

These data provide a resource for mining gene superfamilies, identifying non-coding RNAs, confirming gene losses, and designing experimental constructs. They will be important comparative resources for future genomic studies in echinoderms.

Algebraic Dynamic Programming over general data structures

Background

Dynamic programming algorithms provide exact solutions to many problems in computational biology, such as sequence alignment, RNA folding, hidden Markov models (HMMs), and scoring of phylogenetic trees. Structurally analogous algorithms compute optimal solutions, evaluate score distributions, and perform stochastic sampling. This is explained in the theory of Algebraic Dynamic Programming (ADP) by a strict separation of state space traversal (usually represented by a context free grammar), scoring (encoded as an algebra), and choice rule. A key ingredient in this theory is the use of yield parsers that operate on the ordered input data structure, usually strings or ordered trees. The computation of ensemble properties, such as a posteriori probabilities of HMMs or partition functions in RNA folding, requires the combination of two distinct, but intimately related algorithms, known as the inside and the outside recursion. Only the inside recursions are covered by the classical ADP theory.

Results

The ideas of ADP are generalized to a much wider scope of data structures by relaxing the concept of parsing. This allows us to formalize the conceptual complementarity of inside and outside variables in a natural way. We demonstrate that outside recursions are generically derivable from inside decomposition schemes. In addition to rephrasing the well-known algorithms for HMMs, pairwise sequence alignment, and RNA folding we show how the TSP and the shortest Hamiltonian path problem can be implemented efficiently in the extended ADP framework. As a showcase application we investigate the ancient evolution of HOX gene clusters in terms of shortest Hamiltonian paths.

Conclusions

The generalized ADP framework presented here greatly facilitates the development and implementation of dynamic programming algorithms for a wide spectrum of applications.

Hemichordate genomes and deuterostome origins

Acorn worms, also known as enteropneust (literally, 'gut-breathing') hemichordates, are marine invertebrates that share features with echinoderms and chordates. Together, these three phyla comprise the deuterostomes. Here we report the draft genome sequences of two acorn worms, Saccoglossus kowalevskii and Ptychodera flava. By comparing them with diverse bilaterian genomes, we identify shared traits that were probably inherited from the last common deuterostome ancestor, and then explore evolutionary trajectories leading from this ancestor to hemichordates, echinoderms and chordates. The hemichordate genomes exhibit extensive conserved synteny with amphioxus and other bilaterians, and deeply conserved non-coding sequences that are candidates for conserved gene-regulatory elements. Notably, hemichordates possess a deuterostome-specific genomic cluster of four ordered transcription factor genes, the expression of which is associated with the development of pharyngeal 'gill' slits, the foremost morphological innovation of early deuterostomes, and is probably central to their filter-feeding lifestyle. Comparative analysis reveals numerous deuterostome-specific gene novelties, including genes found in deuterostomes and marine microbes, but not other animals. The putative functions of these genes can be linked to physiological, metabolic and developmental specializations of the filter-feeding ancestor.

Patterning of anteroposterior body axis displayed in the expression of Hox genes in sea cucumber Apostichopus japonicus

The presence of an anteroposterior body axis is a fundamental feature of bilateria. Within this group, echinoderms have secondarily evolved pentameral symmetric body plans. Although all echinoderms present bilaterally symmetric larval stages, they dramatically rearrange their body axis and develop a pentaradial body plan during metamorphosis. Therefore, the location of their anteroposterior body axis in adult forms remains a contentious issue. Unlike other echinoderms, sea cucumbers present an obvious anteroposterior axis not rearranged during metamorphosis, thus representing an interesting group to study their anteroposterior axis patterning. Hox genes are known to play a broadly conserved role in anteroposterior axis patterning in deuterostomes. Here, we report the expression patterns of Hox genes from early development to pentactula stage in sea cucumber. In early larval stages, five Hox genes (AjHox1, AjHox7, AjHox8, AjHox11/13a, and AjHox11/13b) were expressed sequentially along the archenteron, suggesting that the role of anteroposterior patterning of the Hox genes is conserved in bilateral larvae of echinoderms. In doliolaria and pentactula stages, eight Hox genes (AjHox1, AjHox5, AjHox7, AjHox8, AjHox9/10, AjHox11/13a, AjHox11/13b, and AjHox11/13c) were expressed sequentially along the digestive tract, following a similar expression pattern to that found in the visceral mesoderm of other bilateria. Unlike other echinoderms, pentameral expression patterns of AjHox genes were not observed in sea cucumber. Altogether, we concluded that AjHox genes are involved in the patterning of the digestive tract in both larvae and metamorphosis of sea cucumbers. In addition, the anteroposterior axis in sea cucumbers might be patterned like that of other bilateria.

Do echinoderm genomes measure up?

Echinoderm genome sequences are a corpus of useful information about a clade of animals that serve as research models in fields ranging from marine ecology to cell and developmental biology. Genomic information from echinoids has contributed to insights into the gene interactions that drive the developmental process at the molecular level. Such insights often rely heavily on genomic information and the kinds of questions that can be asked thus depend on the quality of the sequence information. Here we describe the history of echinoderm genomic sequence assembly and present details about the quality of the data obtained. All of the sequence information discussed here is posted on the echinoderm information web system, Echinobase.org.

Genome-wide analysis of homeobox genes from Mesobuthus martensii reveals Hox gene duplication in scorpions

Homeobox genes belong to a large gene group, which encodes the famous DNA-binding homeodomain that plays a key role in development and cellular differentiation during embryogenesis in animals. Here, one hundred forty-nine homeobox genes were identified from the Asian scorpion, Mesobuthus martensii (Chelicerata: Arachnida: Scorpiones: Buthidae) based on our newly assembled genome sequence with approximately 248 × coverage. The identified homeobox genes were categorized into eight classes including 82 families: 67 ANTP class genes, 33 PRD genes, 11 LIM genes, five POU genes, six SINE genes, 14 TALE genes, five CUT genes, two ZF genes and six unclassified genes. Transcriptome data confirmed that more than half of the genes were expressed in adults. The homeobox gene diversity of the eight classes is similar to the previously analyzed Mandibulata arthropods. Interestingly, it is hypothesized that the scorpion M. martensii may have two Hox clusters. The first complete genome-wide analysis of homeobox genes in Chelicerata not only reveals the repertoire of scorpion, arachnid and chelicerate homeobox genes, but also shows some insights into the evolution of arthropod homeobox genes.

Phylogenomic analysis of echinoderm class relationships supports Asterozoa

While some aspects of the phylogeny of the five living echinoderm classes are clear, the position of the ophiuroids (brittlestars) relative to asteroids (starfish), echinoids (sea urchins) and holothurians (sea cucumbers) is controversial. Ophiuroids have a pluteus-type larva in common with echinoids giving some support to an ophiuroid/echinoid/holothurian clade named Cryptosyringida. Most molecular phylogenetic studies, however, support an ophiuroid/asteroid clade (Asterozoa) implying either convergent evolution of the pluteus or reversals to an auricularia-type larva in asteroids and holothurians. A recent study of 10 genes from four of the five echinoderm classes used 'phylogenetic signal dissection' to separate alignment positions into subsets of (i) suboptimal, heterogeneously evolving sites (invariant plus rapidly changing) and (ii) the remaining optimal, homogeneously evolving sites. Along with most previous molecular phylogenetic studies, their set of heterogeneous sites, expected to be more prone to systematic error, support Asterozoa. The homogeneous sites, in contrast, support an ophiuroid/echinoid grouping, consistent with the cryptosyringid clade, leading them to posit homology of the ophiopluteus and echinopluteus. Our new dataset comprises 219 genes from all echinoderm classes; analyses using probabilistic Bayesian phylogenetic methods strongly support Asterozoa. The most reliable, slowly evolving quartile of genes also gives highest support for Asterozoa; this support diminishes in second and third quartiles and the fastest changing quartile places the ophiuroids close to the root. Using phylogenetic signal dissection, we find heterogenous sites support an unlikely grouping of Ophiuroidea + Holothuria while homogeneous sites again strongly support Asterozoa. Our large and taxonomically complete dataset finds no support for the cryptosyringid hypothesis; in showing strong support for the Asterozoa, our preferred topology leaves the question of homology of pluteus larvae open.

On a possible evolutionary link of the stomochord of hemichordates to pharyngeal organs of chordates

As a group closely related to chordates, hemichordate acorn worms are in a key phylogenic position for addressing hypotheses of chordate origins. The stomochord of acorn worms is an anterior outgrowth of the pharynx endoderm into the proboscis. In 1886 Bateson proposed homology of this organ to the chordate notochord, crowning this animal group "hemichordates." Although this proposal has been debated for over a century, the question still remains unresolved. Here we review recent progress related to this question. First, the developmental mode of the stomochord completely differs from that of the notochord. Second, comparison of expression profiles of genes including Brachyury, a key regulator of notochord formation in chordates, does not support the stomochord/notochord homology. Third, FoxE that is expressed in the stomochord-forming region in acorn worm juveniles is expressed in the club-shaped gland and in the endostyle of amphioxus, in the endostyle of ascidians, and in the thyroid gland of vertebrates. Based on these findings, together with the anterior endodermal location of the stomochord, we propose that the stomochord has evolutionary relatedness to chordate organs deriving from the anterior pharynx rather than to the notochord.

Genomic organization of Hox and ParaHox clusters in the echinoderm, Acanthaster planci

The organization of echinoderm Hox clusters is of interest due to the role that Hox genes play in deuterostome development and body plan organization, and the unique gene order of the Hox complex in the sea urchin Strongylocentrotus purpuratus, which has been linked to the unique development of the axial region. Here, it has been reported that the Hox and ParaHox clusters of Acanthaster planci, a corallivorous starfish found in the Pacific and Indian oceans, generally resembles the chordate and hemichordate clusters. The A. planci Hox cluster shared with sea urchins the loss of one of the medial Hox genes, even-skipped (Evx) at the anterior of the cluster, as well as organization of the posterior Hox genes.

Evolutionary genomics of fast evolving tunicates

Tunicates have been extensively studied because of their crucial phylogenetic location (the closest living relatives of vertebrates) and particular developmental plan. Recent genome efforts have disclosed that tunicates are also remarkable in their genome organization and molecular evolutionary patterns. Here, we review these latter aspects, comparing the similarities and specificities of two model species of the group: Oikopleura dioica and Ciona intestinalis. These species exhibit great genome plasticity and Oikopleura in particular has undergone a process of extreme genome reduction and compaction that can be explained in part by gene loss, but is mostly due to other mechanisms such as shortening of intergenic distances and introns, and scarcity of mobile elements. In Ciona, genome reorganization was less severe being more similar to the other chordates in several aspects. Rates and patterns of molecular evolution are also peculiar in tunicates, being Ciona about 50% faster than vertebrates and Oikopleura three times faster. In fact, the latter species is considered as the fastest evolving metazoan recorded so far. Two processes of increase in evolutionary rates have taken place in tunicates. One of them is more extreme, and basically restricted to genes encoding regulatory proteins (transcription regulators, chromatin remodeling proteins, and metabolic regulators), and the other one is less pronounced but affects the whole genome. Very likely adaptive evolution has played a very significant role in the first, whereas the functional and/or evolutionary causes of the second are less clear and the evidence is not conclusive. The evidences supporting the incidence of increased mutation and less efficient negative selection are presented and discussed.

How Hox genes can shed light on the place of echinoderms among the deuterostomes

Background

The Hox gene cluster ranks among the greatest of biological discoveries of the past 30 years. Morphogenetic patterning genes are remarkable for the systems they regulate during major ontogenetic events, and for their expressions of molecular, temporal, and spatial colinearity. Recent descriptions of exceptions to these colinearities are suggesting deep phylogenetic signal that can be used to explore origins of entire deuterostome phyla. Among the most enigmatic of these deuterostomes in terms of unique body patterning are the echinoderms. However, there remains no overall synthesis of the correlation between this signal and the variations observable in the presence/absence and expression patterns of Hox genes.

Results

Recent data from Hox cluster analyses shed light on how the bizarre shift from bilateral larvae to radial adults during echinoderm ontogeny can be accomplished by equally radical modifications within the Hox cluster. In order to explore this more fully, a compilation of observations on the genetic patterns among deuterostomes is integrated with the body patterning trajectories seen across the deuterostome clade.

Conclusions

Synthesis of available data helps to explain morphogenesis along the anterior/posterior axis of echinoderms, delineating the origins and fate of that axis during ontogeny. From this, it is easy to distinguish between ‘seriality’ along echinoderm rays and true A/P axis phenomena such as colinearity within the somatocoels, and the ontogenetic outcomes of the unique translocation and inversion of the anterior Hox class found within the Echinodermata. An up-to-date summary and integration of the disparate lines of research so far produced on the relationship between Hox genes and pattern formation for all deuterostomes allows for development of a phylogeny and scenario for the evolution of deuterostomes in general, and the Echinodermata in particular.

Oral-aboral identity displayed in the expression of HpHox3 and HpHox11/13 in the adult rudiment of the sea urchin Holopneustes purpurescens

Hox genes are noted for their roles in specifying axial identity in bilateral forms. In the radial echinoderms, the axis whose identity Hox genes might specify remains unclear. From the expression of Hox genes in the development of the sea urchin Holopneustes purpurescens reported here and that reported previously, we clarify the axis that might be specified by Hox genes in echinoderms. The expression of HpHox11/13 here is described at three developmental stages. The expression is around the rim of the blastopore in gastrulae, in the archenteron wall and adjacent mesoderm in early vestibula larvae, and in a patch of mesoderm close to the archenteron wall in later vestibula larvae. The retained expression of HpHox11/13 in the patch of mesoderm in the later vestibula larvae is, we suggest, indicative of a posterior or an aboral growth zone. The expression of HpHox3 at the echinoid-rudiment stage, in contrast, is in oral mesoderm beneath the epineural folds, concentrated in sites where the first three adult spines form. With the expression of HpHox5 and HpHox11/13 reported previously, the expressions here support the role of Hox genes in specifying oral-aboral identity in echinoderms. How such specification and a posterior growth zone add support to a concept of the structural homology between echinoderms and chordates is discussed.

Homeobox genes expressed during echinoderm arm regeneration

Regeneration in echinoderms has proved to be more amenable to study in the laboratory than the more classical vertebrate models, since the smaller genome size and the absence of multiple orthologs for different genes in echinoderms simplify the analysis of gene function during regeneration. In order to understand the role of homeobox-containing genes during arm regeneration in echinoderms, we isolated the complement of genes belonging to the Hox class that are expressed during this process in two major echinoderm groups: asteroids (Echinaster sepositus and Asterias rubens) and ophiuroids (Amphiura filiformis), both of which show an extraordinary capacity for regeneration. By exploiting the sequence conservation of the homeobox, putative orthologs of several Hox genes belonging to the anterior, medial, and posterior groups were isolated. We also report the isolation of a few Hox-like genes expressed in the same systems.

Evolution of Hox gene clusters in deuterostomes

Hox genes, with their similar roles in animals as evolutionarily distant as humans and flies, have fascinated biologists since their discovery nearly 30 years ago. During the last two decades, reports on Hox genes from a still growing number of eumetazoan species have increased our knowledge on the Hox gene contents of a wide range of animal groups. In this review, we summarize the current Hox inventory among deuterostomes, not only in the well-known teleosts and tetrapods, but also in the earlier vertebrate and invertebrate groups. We draw an updated picture of the ancestral repertoires of the different lineages, a sort of “genome Hox bar-code” for most clades. This scenario allows us to infer differential gene or cluster losses and gains that occurred during deuterostome evolution, which might be causally linked to the morphological changes that led to these widely diverse animal taxa. Finally, we focus on the challenging family of posterior Hox genes, which probably originated through independent tandem duplication events at the origin of each of the ambulacrarian, cephalochordate and vertebrate/urochordate lineages.

Intact cluster and chordate-like expression of ParaHox genes in a sea star

Background

The ParaHox genes are thought to be major players in patterning the gut of several bilaterian taxa. Though this is a fundamental role that these transcription factors play, their activities are not limited to the endoderm and extend to both ectodermal and mesodermal tissues. Three genes compose the ParaHox group: Gsx, Xlox and Cdx. In some taxa (mostly chordates but to some degree also in protostomes) the three genes are arranged into a genomic cluster, in a similar fashion to what has been shown for the better-known Hox genes. Sea urchins possess the full complement of ParaHox genes but they are all dispersed throughout the genome, an arrangement that, perhaps, represented the primitive condition for all echinoderms. In order to understand the evolutionary history of this group of genes we cloned and characterized all ParaHox genes, studied their expression patterns and identified their genomic loci in a member of an earlier branching group of echinoderms, the asteroid Patiria miniata.

Results

We identified the three ParaHox orthologs in the genome of P. miniata. While one of them, PmGsx is provided as maternal message, with no zygotic activation afterwards, the other two, PmLox and PmCdx are expressed during embryogenesis, within restricted domains of both endoderm and ectoderm. Screening of a Patiria bacterial artificial chromosome (BAC) library led to the identification of a clone containing the three genes. The transcriptional directions of PmGsx and PmLox are opposed to that of the PmCdx gene within the cluster.

Conclusions

The identification of P. miniata ParaHox genes has revealed the fact that these genes are clustered in the genome, in contrast to what has been reported for echinoids. Since the presence of an intact cluster, or at least a partial cluster, has been reported in chordates and polychaetes respectively, it becomes clear that within echinoderms, sea urchins have modified the original bilaterian arrangement. Moreover, the sea star ParaHox domains of expression show chordate-like features not found in the sea urchin, confirming that the dynamics of gene expression for the respective genes and their putative regulatory interactions have clearly changed over evolutionary time within the echinoid lineage.