Hox clusters and bilaterian phylogeny

Irreducible Complexity of Hox Gene: Path to the Canonical Function of the Hox Cluster

The evolution of major taxa is often associated with the emergence of new gene families. In all multicellular animals except sponges and comb jellies, the genomes contain Hox genes, which are crucial regulators of development. The canonical function of Hox genes involves colinear patterning of body parts in bilateral animals. This general function is implemented through complex, precisely coordinated mechanisms, not all of which are evolutionarily conserved and fully understood. We suggest that the emergence of this regulatory complexity was preceded by a stage of cooperation between more ancient morphogenetic programs or their individual elements. Footprints of these programs may be present in modern animals to execute non-canonical Hox functions. Non-canonical functions of Hox genes are involved in maintaining terminal nerve cell specificity, autophagy, oogenesis, pre-gastrulation embryogenesis, vertical signaling, and a number of general biological processes. These functions are realized by the basic properties of homeodomain protein and could have triggered the evolution of ParaHoxozoa and Nephrozoa subsequently. Some of these non-canonical Hox functions are discussed in our review.

Gene Loss may have Shaped the Cnidarian and Bilaterian Hox and ParaHox Complement

Hox and ParaHox transcription factors are important for specifying cell fates along the primary body axes during the development of most animals. Within Cnidaria, much of the research on Hox/ParaHox genes has focused on Anthozoa (anemones and corals) and Hydrozoa (hydroids) and has concentrated on the evolution and function of cnidarian Hox genes in relation to their bilaterian counterparts. Here we analyze together the full complement of Hox and ParaHox genes from species representing all four medusozoan classes (Staurozoa, Cubozoa, Hydrozoa, and Scyphozoa) and both anthozoan classes (Octocorallia and Hexacorallia). Our results show that Hox genes involved in patterning the directive axes of anthozoan polyps are absent in the stem leading to Medusozoa. For the first time, we show spatial and temporal expression patterns of Hox and ParaHox genes in the upside-down jellyfish Cassiopea xamachana (Scyphozoa), which are consistent with diversification of medusozoan Hox genes both from anthozoans and within medusozoa. Despite unprecedented taxon sampling, our phylogenetic analyses, like previous studies, are characterized by a lack of clear homology between most cnidarian and bilaterian Hox and Hox-related genes. Unlike previous studies, we propose the hypothesis that the cnidarian-bilaterian ancestor possessed a remarkably large Hox complement and that extensive loss of Hox genes was experienced by both cnidarian and bilaterian lineages.

Seeking Sense in the Hox Gene Cluster

The Hox gene cluster, responsible for patterning of the head-tail axis, is an ancestral feature of all bilaterally symmetrical animals (the Bilateria) that remains intact in a wide range of species. We can say that the Hox cluster evolved successfully only once since it is commonly the same in all groups, with labial-like genes at one end of the cluster expressed in the anterior embryo, and Abd-B-like genes at the other end of the cluster expressed posteriorly. This review attempts to make sense of the Hox gene cluster and to address the following questions. How did the Hox cluster form in the protostome-deuterostome last common ancestor, and why was this with a particular head-tail polarity? Why is gene clustering usually maintained? Why is there collinearity between the order of genes along the cluster and the positions of their expressions along the embryo? Why do the Hox gene expression domains overlap along the embryo? Why have vertebrates duplicated the Hox cluster? Why do Hox gene knockouts typically result in anterior homeotic transformations? How do animals adapt their Hox clusters to evolve new structural patterns along the head-tail axis?

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.

CRISPR/Cas9-Mediated Mutagenesis of Abdominal-A and Ultrabithorax in the Asian Corn Borer, Ostrinia furnacalis

Simple Summary

Homeotic genes encode transcription factors that coordinated the anatomical structure formation during the early embryonic development of organisms. In this study, we functionally characterized two homeotic genes, Abdominal-A (Abd-A) and Ultrabithorax (Ubx), in the Asian corn borer, Ostrinia furnacalis (a maize pest that has devastated the Asia-Pacific region) by using a CRISPR/Cas9 genome editing system. Our results show that the mutagenesis of OfAbd-A and OfUbx led to severe morphological defects in O. furnacalis, which included fused segments and segmental twist during the larval stage, and hollowed and incision-like segments during the pupal stage in OfAbd-A mutants, as well as defects in the wing-pad development in pupal and adult OfUbx mutants. Overall, knocking out Abd-A and Ubx in O. furnacalis resulted in the embryonic lethality to, and pleiotropic impact on, other homeotic genes. This study not only confirms the conserved body planning functions in OfAbd-A and OfUbx, but it also strengthens the control implications of these homeotic genes for lepidopteran pests.

Abstract

(1) Background: Abdominal-A (Abd-A) and Ultrabithorax (Ubx) are homeotic genes that determine the identity and morphology of the thorax and abdomen in insects. The Asian corn borer, Ostrinia furnacalis (Guenée) (Lepidoptera: Pyralidae), is a devastating maize pest throughout Asia, the Western Pacific, and Australia. Building on previous knowledge, we hypothesized that the knockout of Abd-A and Ubx would disrupt the abdominal body planning in O. furnacalis. (2) Methods: CRISPR/Cas9-targeted mutagenesis was employed to decipher the functions of these homeotic genes. (3) Results: Knockout insects demonstrated classical homeotic transformations. Specifically, the mutagenesis of OfAbd-A resulted in: (1) Fused segments and segmental twist during the larval stage; (2) Embryonic lethality; and (3) The pleiotropic upregulation of other homeotic genes, including Lab, Pd, Dfd, Antp, and Abd-B. The mutagenesis of OfUbx led to: (1) Severe defects in the wing pads, which limited the ability of the adults to fly and mate; (2) Female sterility; and (3) The pleiotropic upregulation of other homeotic genes, including Dfd, Abd-B, and Wnt1. (4) Conclusions: These combined results not only support our hypothesis, but they also strengthen the potential of using homeotic genes as molecular targets for the genetic control of this global insect pest.

Essay the (unusual) heuristic value of Hox gene clusters; a matter of time?

Ever since their first report in 1984, Antennapedia-type homeobox (Hox) genes have been involved in such a series of interesting observations, in particular due to their conserved clustered organization between vertebrates and arthropods, that one may legitimately wonder about the origin of this heuristic value. In this essay, I first consider different examples where Hox gene clusters have been instrumental in providing conceptual advances, taken from various fields of research and mostly involving vertebrate embryos. These examples touch upon our understanding of genomic evolution, the revisiting of 19th century views on the relationships between development and evolution and the building of a new framework to understand long-range and pleiotropic gene regulation during development. I then discuss whether the high value of the Hox gene family, when considered as an epistemic object, is related to its clustered structure (and the absence thereof in some animal species) and, if so, what is it in such particular genetic oddities that made them so generous in providing the scientific community with interesting information.

The genetic factors of bilaterian evolution

The Cambrian explosion was a unique animal radiation ~540 million years ago that produced the full range of body plans across bilaterians. The genetic mechanisms underlying these events are unknown, leaving a fundamental question in evolutionary biology unanswered. Using large-scale comparative genomics and advanced orthology evaluation techniques, we identified 157 bilaterian-specific genes. They include the entire Nodal pathway, a key regulator of mesoderm development and left-right axis specification; components for nervous system development, including a suite of G-protein-coupled receptors that control physiology and behaviour, the Robo-Slit midline repulsion system, and the neurotrophin signalling system; a high number of zinc finger transcription factors; and novel factors that previously escaped attention. Contradicting the current view, our study reveals that genes with bilaterian origin are robustly associated with key features in extant bilaterians, suggesting a causal relationship.

Dorsoventral decoupling of Hox gene expression underpins the diversification of molluscs

In contrast to the Hox genes in arthropods and vertebrates, those in molluscs show diverse expression patterns with differences reported among lineages. Here, we investigate 2 phylogenetically distant molluscs, a gastropod and a polyplacophoran, and show that the Hox expression in both species can be divided into 2 categories. The Hox expression in the ventral ectoderm generally shows a canonical staggered pattern comparable to the patterns of other bilaterians and likely contributes to ventral patterning, such as neurogenesis. The other category of Hox expression on the dorsal side is strongly correlated with shell formation and exhibits lineage-specific characteristics in each class of mollusc. This generalized model of decoupled dorsoventral Hox expression is compatible with known Hox expression data from other molluscan lineages and may represent a key characteristic of molluscan Hox expression. These results support the concept of widespread staggered Hox expression in Mollusca and reveal aspects that may be related to the evolutionary diversification of molluscs. We propose that dorsoventral decoupling of Hox expression allowed lineage-specific dorsal and ventral patterning, which may have facilitated the evolution of diverse body plans in different molluscan lineages.

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.

Clustered brachiopod Hox genes are not expressed collinearly and are associated with lophotrochozoan novelties

Temporal collinearity is often considered the main force preserving Hox gene clusters in animal genomes. Studies that combine genomic and gene expression data are scarce, however, particularly in invertebrates like the Lophotrochozoa. As a result, the temporal collinearity hypothesis is currently built on poorly supported foundations. Here we characterize the complement, cluster, and expression of Hox genes in two brachiopod species, Terebratalia transversa and Novocrania anomalaT. transversa has a split cluster with 10 genes (lab, pb, Hox3, Dfd, Scr, Lox5, Antp, Lox4, Post2, and Post1), whereas N. anomala has 9 genes (apparently missing Post1). Our in situ hybridization, real-time quantitative PCR, and stage-specific transcriptomic analyses show that brachiopod Hox genes are neither strictly temporally nor spatially collinear; only pb (in T. transversa), Hox3 (in both brachiopods), and Dfd (in both brachiopods) show staggered mesodermal expression. Thus, our findings support the idea that temporal collinearity might contribute to keeping Hox genes clustered. Remarkably, expression of the Hox genes in both brachiopod species demonstrates cooption of Hox genes in the chaetae and shell fields, two major lophotrochozoan morphological novelties. The shared and specific expression of Hox genes, together with Arx, Zic, and Notch pathway components in chaetae and shell fields in brachiopods, mollusks, and annelids provide molecular evidence supporting the conservation of the molecular basis for these lophotrochozoan hallmarks.

Comparative transcriptomics enlarges the toolkit of known developmental genes in mollusks

BACKGROUND

Mollusks display a striking morphological disparity, including, among others, worm-like animals (the aplacophorans), snails and slugs, bivalves, and cephalopods. This phenotypic diversity renders them ideal for studies into animal evolution. Despite being one of the most species-rich phyla, molecular and in silico studies concerning specific key developmental gene families are still scarce, thus hampering deeper insights into the molecular machinery that governs the development and evolution of the various molluscan class-level taxa.

RESULTS

Next-generation sequencing was used to retrieve transcriptomes of representatives of seven out of the eight recent class-level taxa of mollusks. Similarity searches, phylogenetic inferences, and a detailed manual curation were used to identify and confirm the orthology of numerous molluscan Hox and ParaHox genes, which resulted in a comprehensive catalog that highlights the evolution of these genes in Mollusca and other metazoans. The identification of a specific molluscan motif in the Hox paralog group 5 and a lophotrochozoan ParaHox motif in the Gsx gene is described. Functional analyses using KEGG and GO tools enabled a detailed description of key developmental genes expressed in important pathways such as Hedgehog, Wnt, and Notch during development of the respective species. The KEGG analysis revealed Wnt8, Wnt11, and Wnt16 as Wnt genes hitherto not reported for mollusks, thereby enlarging the known Wnt complement of the phylum. In addition, novel Hedgehog (Hh)-related genes were identified in the gastropod Lottia cf. kogamogai, demonstrating a more complex gene content in this species than in other mollusks.

CONCLUSIONS

The use of de novo transcriptome assembly and well-designed in silico protocols proved to be a robust approach for surveying and mining large sequence data in a wide range of non-model mollusks. The data presented herein constitute only a small fraction of the information retrieved from the analysed molluscan transcriptomes, which can be promptly employed in the identification of novel genes and gene families, phylogenetic inferences, and other studies using molecular tools. As such, our study provides an important framework for understanding some of the underlying molecular mechanisms involved in molluscan body plan diversification and hints towards functions of key developmental genes in molluscan morphogenesis.

HOX gene complement and expression in the planarian Schmidtea mediterranea

Background

Freshwater planarians are well known for their regenerative abilities. Less well known is how planarians maintain spatial patterning in long-lived adult animals or how they re-pattern tissues during regeneration. HOX genes are good candidates to regulate planarian spatial patterning, yet the full complement or genomic clustering of planarian HOX genes has not yet been described, primarily because only a few have been detectable by in situ hybridization, and none have given morphological phenotypes when knocked down by RNAi.

Results

Because the planarian Schmidteamediterranea (S. mediterranea) is unsegmented, appendage less, and morphologically simple, it has been proposed that it may have a simplified HOX gene complement. Here, we argue against this hypothesis and show that S. mediterranea has a total of 13 HOX genes, which represent homologs to all major axial categories, and can be detected by whole-mount in situ hybridization using a highly sensitive method. In addition, we show that planarian HOX genes do not cluster in the genome, yet 5/13 have retained aspects of axially restricted expression. Finally, we confirm HOX gene axial expression by RNA deep-sequencing 6 anterior–posterior “zones” of the animal, which we provide as a dataset to the community to discover other axially restricted transcripts.

Conclusions

Freshwater planarians have an unappreciated HOX gene complexity, with all major axial categories represented. However, we conclude based on adult expression patterns that planarians have a derived body plan and their asexual lifestyle may have allowed for large changes in HOX expression from the last common ancestor between arthropods, flatworms, and vertebrates. Using our in situ method and axial zone RNAseq data, it should be possible to further understand the pathways that pattern the anterior–posterior axis of adult planarians.

Electronic supplementary material

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

An Overview of Hox Genes in Lophotrochozoa: Evolution and Functionality

Hox genes are regulators of animal embryonic development. Changes in the number and sequence of Hox genes as well as in their expression patterns have been related to the evolution of the body plan. Lophotrochozoa is a clade of Protostomia characterized by several phyla which show a wide morphological diversity. Despite that the works summarized in this review emphasize the fragmentary nature of the data available regarding the presence and expression of Hox genes, they also offer interesting insight into the evolution of the Hox cluster and the role played by Hox genes in several phyla. However, the number of genes involved in the cluster of the lophotrochozoan ancestor is still a question of debate. The data presented here suggest that at least nine genes were present while two other genes, Lox4 and Post-2, may either have been present in the ancestor or may have arisen as a result of duplication in the Brachiopoda-Mollusca-Annelida lineage. Spatial and temporal collinearity is a feature of Hox gene expression which was probably present in the ancestor of deuterostomes and protostomes. However, in Lophotrochozoa, it has been detected in only a few species belonging to Annelida and Mollusca.

The TALE face of Hox proteins in animal evolution

Hox genes are major regulators of embryonic development. One of their most conserved functions is to coordinate the formation of specific body structures along the anterior-posterior (AP) axis in Bilateria. This architectural role was at the basis of several morphological innovations across bilaterian evolution. In this review, we traced the origin of the Hox patterning system by considering the partnership with PBC and Meis proteins. PBC and Meis belong to the TALE-class of homeodomain-containing transcription factors and act as generic cofactors of Hox proteins for AP axis patterning in Bilateria. Recent data indicate that Hox proteins acquired the ability to interact with their TALE partners in the last common ancestor of Bilateria and Cnidaria. These interactions relied initially on a short peptide motif called hexapeptide (HX), which is present in Hox and non-Hox protein families. Remarkably, Hox proteins can also recruit the TALE cofactors by using specific PBC Interaction Motifs (SPIMs). We describe how a functional Hox/TALE patterning system emerged in eumetazoans through the acquisition of SPIMs. We anticipate that interaction flexibility could be found in other patterning systems, being at the heart of the astonishing morphological diversity observed in the animal kingdom.

Unexpected co-linearity of Hox gene expression in an aculiferan mollusk

BACKGROUND

Mollusca is an extremely diverse animal phylum that includes the aculiferans (worm-like aplacophorans and eight-shelled polyplacophorans) and their sister group, the conchiferans, comprising monoplacophorans, bivalves (clams, mussels), gastropods (snails, slugs), scaphopods (tusk shells) and cephalopods (squids, octopuses). Studies on mollusks have revealed an overall number of 11 Hox genes in seven out of eight molluscan "class"-level taxa, but expression data of key developmental regulators such as homeotic genes are only available for three gastropod and two cephalopod species. These show that Hox genes are involved in the formation of specific features including shell, foot, funnel or tentacles and not in antero-posterior body plan patterning as in most other bilaterian animals. The role of Hox genes in non-conchiferan (i.e., aculiferan) mollusks remains entirely unknown.

RESULTS

Here we present the first data on the expression of seven Hox genes in apolyplacophoran mollusk, Acanthochitona crinita. In A. crinita the Hox genes Acr-Hox1-5, Hox7 and Post2 are expressed in a co-linear pattern along the antero-posterior axis, but not in molluscan-specific features such as the shell or the foot. The expression pattern is restricted to the post-trochal region and the transcripts are present in ecto-, endo- and mesodermal cell layers. Contrary to the situation in gastropods and cephalopods, we did neither find Hox gene expression in distinct neural subsets of A. crinita, nor in its developing shell plates.

CONCLUSIONS

Our analysis and comparison with other lophotrochozoans indicate that the basal role of Hox genes is in antero-posterior axis patterning in mollusks, similar to the vast majority of bilaterian animals, and that this role has been conserved in polyplacophorans, while co-option into patterning of evolutionary novelties emerged either at the base of Conchifera or independently in gastropods and cephalopods. These morphological innovations most likely contributed to the evolutionary success of its representatives, as exemplified by, e.g., the wide ecological range and species richness of gastropods.

Inhibitory activities of short linear motifs underlie Hox interactome specificity in vivo

Hox proteins are well-established developmental regulators that coordinate cell fate and morphogenesis throughout embryogenesis. In contrast, our knowledge of their specific molecular modes of action is limited to the interaction with few cofactors. Here, we show that Hox proteins are able to interact with a wide range of transcription factors in the live Drosophila embryo. In this context, specificity relies on a versatile usage of conserved short linear motifs (SLiMs), which, surprisingly, often restrains the interaction potential of Hox proteins. This novel buffering activity of SLiMs was observed in different tissues and found in Hox proteins from cnidarian to mouse species. Although these interactions remain to be analysed in the context of endogenous Hox regulatory activities, our observations challenge the traditional role assigned to SLiMs and provide an alternative concept to explain how Hox interactome specificity could be achieved during the embryonic development.

DOI:http://dx.doi.org/10.7554/eLife.06034.001

In all animals, it is important that cells are correctly organised into tissues and organs. This organisation starts in the embryo, and cells are instructed to perform different roles depending on their position within the body.

A family of proteins called the Hox proteins coordinates the organisation of the cells in the animal embryo by binding to and controlling the expression of specific genes. To properly control their target genes, Hox proteins need to interact with other proteins called transcription factors that can also bind to the genes. However, only a few of these transcription factors have been identified so far, and it is not clear how Hox proteins are able to interact with them.

Here, Baëza, Viala, Heim et al. identified several more transcription factors that can bind to the Hox proteins in fruit fly embryos. The experiments show that Hox proteins are able to bind to many transcription factors that are very different from each other. Baëza, Viala, Heim et al. also show that two short sections within the Hox proteins known as short linear motifs are important for controlling these interactions. A fly Hox protein that was missing these motifs was able to interact with new transcription factors. This inhibitory role was found in Hox proteins from mice and sea anemones, suggesting that these motifs may play the same role in all animals.

Baëza, Viala, Heim et al.'s findings challenge the traditional view of the role of the short linear motifs in interactions between proteins. Also, the findings provide an alternative explanation for how the Hox proteins are only able to interact with particular transcription factors in animal embryos. The next step will be to find out whether the inhibitory role of short linear motifs could more generally apply to many other protein families.

DOI:http://dx.doi.org/10.7554/eLife.06034.002

Hox genes pattern the anterior-posterior axis of the juvenile but not the larva in a maximally indirect developing invertebrate, Micrura alaskensis (Nemertea)

BACKGROUND

The pilidium larva is a novel body plan that arose within a single clade in the phylum Nemertea - the Pilidiophora. While the sister clade of the Pilidiophora and the basal nemerteans develop directly, pilidiophorans have a long-lived planktotrophic larva with a body plan distinctly different from that of the juvenile. Uniquely, the pilidiophoran juvenile develops inside the larva from several discrete rudiments. The orientation of the juvenile with respect to the larval body varies within the Pilidiophora, which suggests that the larval and juvenile anteroposterior (AP) axes are patterned differently. In order to gain insight into the evolutionary origins of the pilidium larva and the mechanisms underlying this implied axial uncoupling, we examined the expression of the Hox genes during development of the pilidiophoran Micrura alaskensis.

RESULTS

We identified sequences of nine Hox genes and the ParaHox gene caudal through a combination of transcriptome analysis and molecular cloning, and determined their expression pattern during development using in situ hybridization in whole-mounted larvae. We found that Hox genes are first expressed long after the pilidium is fully formed and functional. The Hox genes are expressed in apparently overlapping domains along the AP axis of the developing juvenile in a subset of the rudiments that give rise to the juvenile trunk. Hox genes are not expressed in the larval body at any stage of development.

CONCLUSIONS

While the Hox genes pattern the juvenile pilidiophoran, the pilidial body, which appears to be an evolutionary novelty, must be patterned by some mechanism other than the Hox genes. Although the pilidiophoran juvenile develops from separate rudiments with no obvious relationship to the embryonic formation of the larva, the Hox genes appear to exhibit canonical expression along the juvenile AP axis. This suggests that the Hox patterning system can maintain conserved function even when widely decoupled from early polarity established in the egg.

Tracking context-specific transcription factors regulating hox activity

BACKGROUND

Hox proteins are key developmental regulators involved in almost every embryonic tissue for specifying cell fates along longitudinal axes or during organ formation. It is thought that the panoply of Hox activities relies on interactions with tissue-, stage-, and/or cell-specific transcription factors. High-throughput approaches in yeast or cell culture systems have shown that Hox proteins bind to various types of nuclear and cytoplasmic components, illustrating their remarkable potential to influence many different cell regulatory processes. However, these approaches failed to identify a relevant number of context-specific transcriptional partners, suggesting that these interactions are hard to uncover in non-physiological conditions. Here we discuss this problematic.

RESULTS

In this review, we present intrinsic Hox molecular signatures that are probably involved in multiple (yet specific) interactions with transcriptional partners. We also recapitulate the current knowledge on Hox cofactors, highlighting the difficulty to tracking context-specific cofactors through traditional large-scale approaches.

CONCLUSION

We propose experimental approaches that will allow a better characterisation of interaction networks underlying Hox contextual activities in the next future.

Hox proteins display a common and ancestral ability to diversify their interaction mode with the PBC class cofactors

Hox protein function during development and evolution relies on conserved multiple interaction modes with cofactors of the PBC and Meis families.

Hox transcription factors control a number of developmental processes with the help of the PBC class proteins. In vitro analyses have established that the formation of Hox/PBC complexes relies on a short conserved Hox protein motif called the hexapeptide (HX). This paradigm is at the basis of the vast majority of experimental approaches dedicated to the study of Hox protein function. Here we questioned the unique and general use of the HX for PBC recruitment by using the Bimolecular Fluorescence Complementation (BiFC) assay. This method allows analyzing Hox-PBC interactions in vivo and at a genome-wide scale. We found that the HX is dispensable for PBC recruitment in the majority of investigated Drosophila and mouse Hox proteins. We showed that HX-independent interaction modes are uncovered by the presence of Meis class cofactors, a property which was also observed with Hox proteins of the cnidarian sea anemone Nematostella vectensis. Finally, we revealed that paralog-specific motifs convey major PBC-recruiting functions in Drosophila Hox proteins. Altogether, our results highlight that flexibility in Hox-PBC interactions is an ancestral and evolutionary conserved character, which has strong implications for the understanding of Hox protein functions during normal development and pathologic processes.

Hox proteins are key transcriptional regulators of animal development, famously helping to determine identity along the anterior-posterior body axis. Although their evolution and developmental roles are well established, the molecular mechanisms underlying their specific functions remain poorly characterized. The current dominant view is that interaction with different members of the PBC family of transcription factors confers specific DNA-binding properties on different Hox proteins. However, this idea conflicts with in vitro evidence that a short “hexapeptide” (HX) motif shared by most Hox proteins is solely responsible for generic PBC recruitment. Here we have used the BiFC (bimolecular fluorescence complementation) method to address the global importance of the HX motif for Hox-PBC interactions in living cells and living animals including fruit flies and chick embryos. We observe that most interactions between Hox and PBC proteins do not depend on HX, and that alternative protein motifs are widely used for PBC recruitment in vivo. We also show that DNA binding by a second family of cofactors, the Meis proteins, unmasks these alternative interaction modes and that this property is conserved not only across Bilateria, but also in the basal animal phylum Cnidaria. Taken together, our results demonstrate that Hox-PBC partnership relies on multiple interaction modes, which can be influenced by additional transcriptional partners. We propose that this ancestral feature has been essential for ensuring Hox functional plasticity during development and evolution.

The origin of patterning systems in bilateria-insights from the Hox and ParaHox genes in Acoelomorpha

Hox and ParaHox genes constitute two families of developmental regulators that pattern the Anterior–Posterior body axis in all bilaterians. The members of these two groups of genes are usually arranged in genomic clusters and work in a coordinated fashion, both in space and in time. While the mechanistic aspects of their action are relatively well known, it is still unclear how these systems evolved. For instance, we still need a proper model of how the Hox and ParaHox clusters were assembled over time. This problem is due to the shortage of information on gene complements for many taxa (mainly basal metazoans) and the lack of a consensus phylogenetic model of animal relationships to which we can relate our new findings. Recently, several studies have shown that the Acoelomorpha most probably represent the first offshoot of the Bilateria. This finding has prompted us, and others, to study the Hox and ParaHox complements in these animals, as well as their activity during development. In this review, we analyze how the current knowledge of Hox and ParaHox genes in the Acoelomorpha is shaping our view of bilaterian evolution.

Evolution of invertebrate deuterostomes and Hox/ParaHox genes

Transcription factors encoded by Antennapedia-class homeobox genes play crucial roles in controlling development of animals, and are often found clustered in animal genomes. The Hox and ParaHox gene clusters have been regarded as evolutionary sisters and evolved from a putative common ancestral gene complex, the ProtoHox cluster, prior to the divergence of the Cnidaria and Bilateria (bilaterally symmetrical animals). The Deuterostomia is a monophyletic group of animals that belongs to the Bilateria, and a sister group to the Protostomia. The deuterostomes include the vertebrates (to which we belong), invertebrate chordates, hemichordates, echinoderms and possibly xenoturbellids, as well as acoelomorphs. The studies of Hox and ParaHox genes provide insights into the origin and subsequent evolution of the bilaterian animals. Recently, it becomes apparent that among the Hox and ParaHox genes, there are significant variations in organization on the chromosome, expression pattern, and function. In this review, focusing on invertebrate deuterostomes, I first summarize recent findings about Hox and ParaHox genes. Next, citing unsolved issues, I try to provide clues that might allow us to reconstruct the common ancestor of deuterostomes, as well as understand the roles of Hox and ParaHox genes in the development and evolution of deuterostomes.

Selection of distinct Hox-Extradenticle interaction modes fine-tunes Hox protein activity

Hox genes encode transcription factors widely used for diversifying animal body plans in development and evolution. To achieve functional specificity, Hox proteins associate with PBC class proteins, Pre-B cell leukemia homeobox (Pbx) in vertebrates, and Extradenticle (Exd) in Drosophila, and were thought to use a unique hexapeptide-dependent generic mode of interaction. Recent findings, however, revealed the existence of an alternative, UbdA-dependent paralog-specific interaction mode providing diversity in Hox-PBC interactions. In this study, we investigated the basis for the selection of one of these two Hox-PBC interaction modes. Using naturally occurring variations and mutations in the Drosophila Ultrabithorax protein, we found that the linker region, a short domain separating the hexapeptide from the homeodomain, promotes an interaction mediated by the UbdA domain in a context-dependent manner. While using a UbdA-dependent interaction for the repression of the limb-promoting gene Distalless, interaction with Exd during segment-identity specification still relies on the hexapeptide motif. We further show that distinctly assembled Hox-PBC complexes display subtle but distinct repressive activities. These findings identify Hox-PBC interaction as a template for subtle regulation of Hox protein activity that may have played a major role in the diversification of Hox protein function in development and evolution.

Expression of Hox genes during the larval development of the snail, Gibbula varia (L.)-further evidence of non-colinearity in molluscs

Hox transcription factors, a subfamily of homeobox genes, are expressed in distinct, often overlapping domains along the anterior-posterior body axis of animal embryos. Here, we report the sequence and expression pattern of Hox2, Hox3, Hox4, Hox5, Lox5, Hox7, Lox4, and Lox2 in different larval stages during the encapsulated development of the marine gastropod Gibbula varia. Our results show that all Gva-Hox genes are expressed in ectoderm-derived cells. Hox2, Hox3, Hox4, Hox5, and Hox7 are expressed in overlapping patterns in the pedal, pleural, oesophageal, and visceral ganglia, supporting the ancestral role of Hox genes in the neurogenesis processes in bilaterians. Gva-Hox1, Gva-Post2, and Gva-Post1 genes are involved in shell morphogenesis and have apparently lost their role in neurogangliogenesis. Lox5, Lox4, and Lox2 are expressed in different cells of the apical organ during the earlier larval stage (trochophore) and the cerebral ganglia during later larval stages (veliger). These results support the hypothesis that apical organ neurosensory cells contribute to the formation of cerebral ganglia commissures during metamorphosis. Gva-Hox7 and Gva-Lox4 are additionally expressed in the prototroch of the trochophore and in the velar area of the veliger larvae. This contradicts with the expression of these genes in the annelids, where most of Hox genes are expressed in the posttrochal area and are involved in segmental determination. Therefore, expression of Hox genes may serve as an example of co-option and plasticity of gene function during evolution of gastropods.

Are the deuterostome posterior Hox genes a fast-evolving class?

There has been a great deal of interest in analysing the molecular evolution of the Hox cluster using both bioinformatic and experimental approaches. The posterior Hox genes have been of particular interest to both groups of biologists for a number of reasons: they appear to be associated with the evolution of a number of morphological novelties; the protostomes appear to be have lost a highly-conserved and functionally important amino acid motif (the hexapeptide motif) from their posterior Hox genes; and deuterostome posterior Hox genes seem to be evolving more quickly than all other Hox genes. In this chapter I will discuss the last of these points. The idea that Deuterostome posterior Hox genes were evolving more quickly than other Hox genes was first suggested by David Ferrier and colleagues. In this chapter, I start by introducing the posterior Hox genes--their distribution among the animal phyla and the likely sequence of duplications that led to this distribution. I then introduce the idea of 'deuterostome posterior flexibility' and examine this hypothesis in light of more recent phylogenetic and genomic work on the Hox cluster. Finally, I discuss some new approaches that could be used to test directly for differential rates of evolution among Hox genes and to assess what might underlie these differences.

Evolution of Hox complexes

Recent years have seen a plethora of ideas and hypotheses, and lots of debate, about the origin and evolution of the Hox gene cluster. Here I will attempt to summarize these hypotheses, identify their strengths and weaknesses and highlight the types of new data that may lead to further resolution of the competing ideas. The major theme is that Hox genes originated very early in animal evolution and extensive independent duplications occurred in major lineages. Duplications however have not been the only route to change in the composition and structure of the Hox cluster, as extensive gene losses have occurred as well. Indeed it is gene loss that is one of the main obstacles in our understanding of the origin and evolution of Hox clusters. Matters should be improved with wider taxon sampling along with a clearer understanding of how duplicated genes evolve.

A conserved cluster of three PRD-class homeobox genes (homeobrain, rx and orthopedia) in the Cnidaria and Protostomia

Background

Homeobox genes are a superclass of transcription factors with diverse developmental regulatory functions, which are found in plants, fungi and animals. In animals, several Antennapedia (ANTP)-class homeobox genes reside in extremely ancient gene clusters (for example, the Hox, ParaHox, and NKL clusters) and the evolution of these clusters has been implicated in the morphological diversification of animal bodyplans. By contrast, similarly ancient gene clusters have not been reported among the other classes of homeobox genes (that is, the LIM, POU, PRD and SIX classes).

Results

Using a combination of in silico queries and phylogenetic analyses, we found that a cluster of three PRD-class homeobox genes (Homeobrain (hbn), Rax (rx) and Orthopedia (otp)) is present in cnidarians, insects and mollusks (a partial cluster comprising hbn and rx is present in the placozoan Trichoplax adhaerens). We failed to identify this 'HRO' cluster in deuterostomes; in fact, the Homeobrain gene appears to be missing from the chordate genomes we examined, although it is present in hemichordates and echinoderms. To illuminate the ancestral organization and function of this ancient cluster, we mapped the constituent genes against the assembled genome of a model cnidarian, the sea anemone Nematostella vectensis, and characterized their spatiotemporal expression using in situ hybridization. In N. vectensis, these genes reside in a span of 33 kb with the same gene order as previously reported in insects. Comparisons of genomic sequences and expressed sequence tags revealed the presence of alternative transcripts of Nv-otp and two highly unusual protein-coding polymorphisms in the terminal helix of the Nv-rx homeodomain. A population genetic survey revealed the Rx polymorphisms to be widespread in natural populations. During larval development, all three genes are expressed in the ectoderm, in non-overlapping territories along the oral-aboral axis, with distinct temporal expression.

Conclusion

We report the first evidence for a PRD-class homeobox cluster that appears to have been conserved since the time of the cnidarian-bilaterian ancestor, and possibly even earlier, given the presence of a partial cluster in the placozoan Trichoplax. Very similar clusters comprising these three genes exist in Nematostella and diverse protostomes. Interestingly, in chordates, one member of the ancestral cluster (homeobrain) has apparently been lost, and there is no linkage between rx and orthopedia in any of the vertebrates. In Nematostella, the spatial expression of these three genes along the body column is not colinear with their physical order in the cluster but the temporal expression is, therefore, using the terminology that has been applied to the Hox cluster genes, the HRO cluster would appear to exhibit temporal but not spatial colinearity. It remains to be seen whether the mechanisms responsible for the evolutionary conservation of the HRO cluster are the same mechanisms responsible for cohesion of the Hox cluster and other ANTP-class homeobox clusters that have been widely conserved throughout animal evolution.

A non-tree-based comprehensive study of metazoan Hox and ParaHox genes prompts new insights into their origin and evolution

BACKGROUND

Hox and the closely-related ParaHox genes, which emerged prior to the divergence between cnidarians and bilaterians, are the most well-known members of the ancient genetic toolkit that controls embryonic development across all metazoans. Fundamental questions relative to their origin and evolutionary relationships remain however unresolved. We investigate here the evolution of metazoan Hox and ParaHox genes using the HoxPred program that allows the identification of Hox genes without the need of phylogenetic tree reconstructions.

RESULTS

We show that HoxPred provides an efficient and accurate classification of Hox and ParaHox genes in their respective homology groups, including Hox paralogous groups (PGs). We analyzed more than 10,000 sequences from 310 metazoan species, from 6 genome projects and the complete UniProtKB database. The HoxPred program and all results arranged in the Datab'Hox database are freely available at http://cege.vub.ac.be/hoxpred/. Results for the genome-scale studies are coherent with previous studies, and also brings knowledge on the Hox repertoire and clusters for newly-sequenced species. The unprecedented scale of this study and the use of a non-tree-based approach allows unresolved key questions about Hox and ParaHox genes evolution to be addressed.

CONCLUSIONS

Our analysis suggests that the presence of a single type of Posterior Hox genes (PG9-like) is ancestral to bilaterians, and that new Posterior PGs would have arisen in deuterostomes through independent gene duplications. Four types of Central genes would also be ancestral to bilaterians, with two of them, PG6- and PG7-like that gave rise, in protostomes, to the UbdA- and ftz/Antp/Lox5-type genes, respectively. A fifth type of Central genes (PG8) would have emerged in the vertebrate lineage. Our results also suggest the presence of Anterior (PG1 and PG3), Central and Posterior Hox genes in the cnidarians, supporting an ancestral four-gene Hox cluster. In addition, our data support the relationship of the bilaterian ParaHox genes Gsx and Xlox with PG3, and Cdx with the Central genes. Our study therefore indicates three possible models for the origin of Hox and ParaHox in early metazoans, a two-gene (Anterior/PG3--Central/Posterior), a three-gene (Anterior/PG1, Anterior/PG3 and Central/Posterior), or a four-gene (Anterior/PG1--Anterior/PG3--Central--Posterior) ProtoHox cluster.

Regulation of Hox activity: insights from protein motifs

Deciphering the molecular bases of animal body plan construction is a central question in developmental and evolutionary biology. Genome analyses of a number of metazoans indicate that widely conserved regulatory molecules underlie the amazing diversity of animal body plans, suggesting that these molecules are reiteratively used for multiple purposes. Hox proteins constitute a good example of such molecules and provide the framework to address the mechanisms underlying transcriptional specificity and diversity in development and evolution. Here we examine the current knowledge of the molecular bases of Hox-mediated transcriptional control, focusing on how this control is encoded within protein sequences and structures. The survey suggests that the homeodomain is part of an extended multifunctional unit coordinating DNA binding and activity regulation and highlights the need for further advances in our understanding of Hox protein activity.

Early origin of the bilaterian developmental toolkit

Whole-genome sequences from the choanoflagellate Monosiga brevicollis, the placozoan Trichoplax adhaerens and the cnidarian Nematostella vectensis have confirmed results from comparative evolutionary developmental studies that much of the developmental toolkit once thought to be characteristic of bilaterians appeared much earlier in the evolution of animals. The diversity of transcription factors and signalling pathway genes in animals with a limited number of cell types and a restricted developmental repertoire is puzzling, particularly in light of claims that such highly conserved elements among bilaterians provide evidence of a morphologically complex protostome-deuterostome ancestor. Here, I explore the early origination of elements of what became the bilaterian toolkit, and suggest that placozoans and cnidarians represent a depauperate residue of a once more diverse assemblage of early animals, some of which may be represented in the Ediacaran fauna (c. 585-542 Myr ago).

Bilaterian phylogeny: a broad sampling of 13 nuclear genes provides a new Lophotrochozoa phylogeny and supports a paraphyletic basal acoelomorpha

During the past decade, great progress has been made in clarifying the relationships among bilaterian animals. Studies based on a limited number of markers established new hypotheses such as the existence of three superclades (Deuterostomia, Ecdysozoa, and Lophotrochozoa) but left major questions unresolved. The data sets used to the present either bear few characters for many taxa (i.e., the ribosomal genes) or present many characters but lack many phyla (such as recent phylogenomic approaches) failing to provide definitive answers for all the regions of the bilaterian tree. We performed phylogenetic analyses using a molecular matrix with a high number of characters and bilaterian phyla. This data set is built from 13 genes (8,880 bp) belonging to 90 taxa from 27 bilaterian phyla. Probabilistic analyses robustly support the three superclades, the monophyly of Chordata, a spiralian clade including Brachiozoa, the basal position of a paraphyletic Acoelomorpha, and point to an ecdysozoan affiliation for Chaetognatha. This new phylogeny not only agrees with most classical molecular results but also provides new insights into the relationships between lophotrochozoans and challenges the results obtained using high-throughput strategies, highlighting the problems associated with the current trend to increase gene number rather than taxa.

Classification of sequence signatures: a guide to Hox protein function

Hox proteins are part of the conserved superfamily of homeodomain-containing transcription factors and play fundamental roles in shaping animal body plans in development and evolution. However, molecular mechanisms underlying their diverse and specific biological functions remain largely enigmatic. Here, we have analyzed Hox sequences from the main evolutionary branches of the Bilateria group. We have found that four classes of Hox protein signatures exist, which together provide sufficient support to explain how different Hox proteins differ in their control and function. The homeodomain and its surrounding sequences accumulate nearly all signatures, constituting an extended module where most of the information distinguishing Hox proteins is concentrated. Only a small fraction of these signatures has been investigated at the functional level, but these show that approaches relying on Hox protein alterations still have a large potential for deciphering molecular mechanisms of Hox differential control.

Retinoic acid and Wnt/beta-catenin have complementary roles in anterior/posterior patterning embryos of the basal chordate amphioxus

A role for Wnt/beta-catenin signaling in axial patterning has been demonstrated in animals as basal as cnidarians, while roles in axial patterning for retinoic acid (RA) probably evolved in the deuterostomes and may be chordate-specific. In vertebrates, these two pathways interact both directly and indirectly. To investigate the evolutionary origins of interactions between these two pathways, we manipulated Wnt/beta-catenin and RA signaling in the basal chordate amphioxus during the gastrula stage, which is the RA-sensitive period for anterior/posterior (A/P) patterning. The results show that Wnt/beta-catenin and RA signaling have distinctly different roles in patterning the A/P axis of the amphioxus gastrula. Wnt/beta-catenin specifies the identity of the ends of the embryo (high Wnt = posterior; low Wnt = anterior) but not intervening positions. Thus, upregulation of Wnt/beta-catenin signaling induces ectopic expression of posterior markers at the anterior tip of the embryo. In contrast, RA specifies position along the A/P axis, but not the identity of the ends of the embryo-increased RA signaling strongly affects the domains of Hox expression along the A/P axis but has little or no effect on the expression of either anterior or posterior markers. Although the two pathways may both influence such things as specification of neuronal identity, interactions between them in A/P patterning appear to be minimal.

Are Hox genes ancestrally involved in axial patterning? Evidence from the hydrozoan Clytia hemisphaerica (Cnidaria)

BACKGROUND

The early evolution and diversification of Hox-related genes in eumetazoans has been the subject of conflicting hypotheses concerning the evolutionary conservation of their role in axial patterning and the pre-bilaterian origin of the Hox and ParaHox clusters. The diversification of Hox/ParaHox genes clearly predates the origin of bilaterians. However, the existence of a "Hox code" predating the cnidarian-bilaterian ancestor and supporting the deep homology of axes is more controversial. This assumption was mainly based on the interpretation of Hox expression data from the sea anemone, but growing evidence from other cnidarian taxa puts into question this hypothesis.

METHODOLOGY/PRINCIPAL FINDINGS

Hox, ParaHox and Hox-related genes have been investigated here by phylogenetic analysis and in situ hybridisation in Clytia hemisphaerica, an hydrozoan species with medusa and polyp stages alternating in the life cycle. Our phylogenetic analyses do not support an origin of ParaHox and Hox genes by duplication of an ancestral ProtoHox cluster, and reveal a diversification of the cnidarian HOX9-14 genes into three groups called A, B, C. Among the 7 examined genes, only those belonging to the HOX9-14 and the CDX groups exhibit a restricted expression along the oral-aboral axis during development and in the planula larva, while the others are expressed in very specialised areas at the medusa stage.

CONCLUSIONS/SIGNIFICANCE

Cross species comparison reveals a strong variability of gene expression along the oral-aboral axis and during the life cycle among cnidarian lineages. The most parsimonious interpretation is that the Hox code, collinearity and conservative role along the antero-posterior axis are bilaterian innovations.

Lophotrochozoa internal phylogeny: new insights from an up-to-date analysis of nuclear ribosomal genes

Resolving the relationships among animal phyla is a key biological problem that remains to be solved. Morphology is unable to determine the relationships among most phyla and although molecular data have unveiled a new evolutionary scenario, they have their own limitations. Nuclear ribosomal genes (18S and 28S rDNA) have been used effectively for many years. However, they are considered of limited use for resolving deep divergences such as the origin of the bilaterians, due to certain drawbacks such as the long-branch attraction (LBA) problem. Here, we attempt to overcome these pitfalls by combining several methods suggested in previous studies and routinely used in contemporary standard phylogenetic analyses but that have not yet been applied to any bilaterian phylogeny based on these genes. The methods used include maximum likelihood and Bayesian inference, the application of models with rate heterogeneity across sites, wide taxon sampling and compartmentalized analyses for each problematic clade. The results obtained show that the combination of the above-mentioned methodologies minimizes the LBA effect, and a new Lophotrochozoa phylogeny emerges. Also, the Acoela and Nemertodermatida are confirmed with maximum support as the first branching bilaterians. Ribosomal RNA genes are thus a reliable source for the study of deep divergences in the metazoan tree, provided that the data are treated carefully.

Genomic organization and expression demonstrate spatial and temporal Hox gene colinearity in the lophotrochozoan Capitella sp. I

Hox genes define regional identities along the anterior–posterior axis in many animals. In a number of species, Hox genes are clustered in the genome, and the relative order of genes corresponds with position of expression in the body. Previous Hox gene studies in lophotrochozoans have reported expression for only a subset of the Hox gene complement and/or lack detailed genomic organization information, limiting interpretations of spatial and temporal colinearity in this diverse animal clade. We studied expression and genomic organization of the single Hox gene complement in the segmented polychaete annelid Capitella sp. I. Total genome searches identified 11 Hox genes in Capitella, representing 11 distinct paralog groups thought to represent the ancestral lophotrochozoan complement. At least 8 of the 11 Capitella Hox genes are genomically linked in a single cluster, have the same transcriptional orientation, and lack interspersed non-Hox genes. Studying their expression by situ hybridization, we find that the 11 Capitella Hox genes generally exhibit spatial and temporal colinearity. With the exception of CapI-Post1, Capitella Hox genes are all expressed in broad ectodermal domains during larval development, consistent with providing positional information along the anterior–posterior axis. The anterior genes CapI-lab, CapI-pb, and CapI-Hox3 initiate expression prior to the appearance of segments, while more posterior genes appear at or soon after segments appear. Many of the Capitella Hox genes have either an anterior or posterior expression boundary coinciding with the thoracic–abdomen transition, a major body tagma boundary. Following metamorphosis, several expression patterns change, including appearance of distinct posterior boundaries and restriction to the central nervous system. Capitella Hox genes have maintained a clustered organization, are expressed in the canonical anterior–posterior order found in other metazoans, and exhibit spatial and temporal colinearity, reflecting Hox gene characteristics that likely existed in the protostome–deuterostome ancestor.

PCR survey of Xenoturbella bocki Hox genes

Xenoturbella bocki has recently been identified as one of the most basal deuterostomes, although an even more basal phylogenetic position cannot be ruled out. Here we report on a polymerase chain reaction survey of partial Hox homeobox sequences of X. bocki. Surprisingly, we did not find evidence for more than five Hox genes, one clear labial/PG1 ortholog, one posterior gene most similar to the PG9/10 genes of Ambulacraria, and three central group genes whose precise assignment to a specific paralog group remains open. We furthermore report on a re-evaluation of the available published evidence of Hox genes in other basal deuterostomes.