Hox genes are not always Colinear

The functional evolution of collembolan Ubx on the regulation of abdominal appendage formation

Folsomia candida is a tiny soil-living arthropod belonging to the Collembola, which is an outgroup to Insecta. It resembles insects as having a pair of antennae and three pairs of thorax legs, while it also possesses three abdominal appendages: a ventral tube located in the first abdominal segment (A1), a retinaculum in A3, and a furca in A4. Collembolan Ubx and AbdA specify abdominal appendages, but they are unable to repress appendage marker gene Dll. The genetic basis of collembolan appendage formation and the mechanisms by which Ubx and AbdA regulate Dll transcription and appendage development remains unknown. In this study, we analysed the developmental transcriptomes of F. candida and identified candidate appendage formation genes, including Ubx (FcUbx). The expression data revealed the dominance of Dll over Ubx during the embryonic 3.5 and 4.5 days, suggesting that Ubx is deficient in suppressing Dll at early appendage formation stages. Furthermore, via electrophoretic mobility shift assays and dual luciferase assays, we found that the binding and repression capacity of FcUbx on Drosophila Dll resembles those of the longest isoform of Drosophila Ubx (DmUbx_Ib), while the regulatory mechanism of the C-terminus of FcUbx on Dll repression is similar to that of the crustacean Artemia franciscana Ubx (AfUbx), demonstrating that the function of collembolan Ubx is intermediate between that of Insecta and Crustacea. In summary, our study provides novel insights into collembolan appendage formation and sheds light on the functional evolution of Ubx. Additionally, we propose a model that collembolan Ubx regulates abdominal segments in a context-specific manner.

Comparative Hox genes expression within the dimorphic annelid Streblospio benedicti reveals patterning variation during development

Hox genes are transcriptional regulators that elicit cell positional identity along the anterior-posterior region of the body plan across different lineages of Metazoan. Comparison of Hox gene expression across distinct species reveals their evolutionary conservation; however, their gains and losses in different lineages can correlate with body plan modifications and morphological novelty. We compare the expression of 11 Hox genes found within Streblospio benedicti, a marine annelid that produces two types of offspring with distinct developmental and morphological features. For these two distinct larval types, we compare Hox gene expression through ontogeny using hybridization chain reaction (HCR) probes for in situ hybridization and RNA-seq data. We find that Hox gene expression patterning for both types is typically similar at equivalent developmental stages. However, some Hox genes have spatial or temporal differences between the larval types that are associated with morphological and life-history differences. This is the first comparison of developmental divergence in Hox gene expression within a single species and these changes reveal how body plan differences may arise in larval evolution.

40 years of the homeobox: mechanisms of Hox spatial-temporal collinearity in vertebrates

Animal body plans are established during embryonic development by the Hox genes. This patterning process relies on the differential expression of Hox genes along the head-to-tail axis. Hox spatial collinearity refers to the relationship between the organization of Hox genes in clusters and the differential Hox expression, whereby the relative order of the Hox genes within a cluster mirrors the spatial sequence of expression in the developing embryo. In vertebrates, the cluster organization is also associated with the timing of Hox activation, which harmonizes Hox expression with the progressive emergence of axial tissues. Thereby, in vertebrates, Hox temporal collinearity is intimately linked to Hox spatial collinearity. Understanding the mechanisms contributing to Hox temporal and spatial collinearity is thus key to the comprehension of vertebrate patterning. Here, we provide an overview of the main discoveries pertaining to the mechanisms of Hox spatial-temporal collinearity.

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.

Hox, homology, and parsimony: An organismal perspective

Hox genes are important regulators in animal development. They often show a mosaic of conserved (e.g., longitudinal axis patterning) and lineage-specific novel functions (e.g., development of skeletal, sensory, or locomotory systems). Despite extensive research over the past decades, it remains controversial at which node in the animal tree of life the Hox cluster evolved. Its presence already in the last common metazoan ancestor has been proposed, although the genomes of both putative earliest extant metazoan offshoots, the ctenophores and the poriferans, are devoid of Hox sequences. The lack of Hox genes in the supposedly "simple"-built poriferans and their low number in cnidarians and the basally branching bilaterians, the xenacoelomorphs, seems to support the classical notion that the number of Hox genes is correlated with the degree of animal complexity. However, the 4-fold increase of the Hox cluster in xiphosurans, a basally branching chelicerate clade, as well as the situation in some teleost fishes that show a multitude of Hox genes compared to, e.g., human, demonstrates, that there is no per se direct correlation between organismal complexity and Hox number. Traditional approaches have tried to base homology on the morphological level on shared expression profiles of individual genes, but recent data have shown that, in particular with respect to Hox and other regulatory genes, complex gene-gene interactions rather than expression signatures of individual genes alone are responsible for shaping morphological traits during ontogeny. Accordingly, for sound homology assessments and reconstructions of character evolution on organ system level, additional independent datasets (e.g., morphological, developmental) need to be included in any such analyses. If supported by solid data, proposed structural homology should be regarded as valid and not be rejected solely on the grounds of non-parsimonious distribution of the character over a given phylogenetic topology.

Comparative Hox genes expression within the dimorphic annelid Streblospio benedicti reveals patterning variation during development

Hox genes are transcriptional regulators that elicit cell positional identity along the anterior-posterior region of the body plan across different lineages of Metazoan. Comparison of Hox gene expression across distinct species reveals their evolutionary conservation, however their gains and losses in different lineages can correlate with body plan modifications and morphological novelty. We compare the expression of eleven Hox genes found within Streblospio benedicti, a marine annelid that produces two types of offspring with distinct developmental and morphological features. For these two distinct larval types, we compare Hox gene expression through ontogeny using HCR (hybridization chain reaction) probes for in-situ hybridization and RNA-seq data. We find that Hox gene expression patterning for both types is typically similar at equivalent developmental stages. However, some Hox genes have spatial or temporal differences between the larval types that are associated with morphological and life-history differences. This is the first comparison of developmental divergence in Hox genes expression within a single species and these changes reveal how body plan differences may arise in larval evolution.

Emerging questions on the mechanisms and dynamics of 3D genome evolution in spiralians

Information on how 3D genome topology emerged in animal evolution, how stable it is during development, its role in the evolution of phenotypic novelties and how exactly it affects gene expression is highly debated. So far, data to address these questions are lacking with the exception of a few key model species. Several gene regulatory mechanisms have been proposed, including scenarios where genome topology has little to no impact on gene expression, and vice versa. The ancient and diverse clade of spiralians may provide a crucial testing ground for such mechanisms. Sprialians have followed distinct evolutionary trajectories, with some clades experiencing genome expansions and/or large-scale genome rearrangements, and others undergoing genome contraction, substantially impacting their size and organisation. These changes have been associated with many phenotypic innovations in this clade. In this review, we describe how emerging genome topology data, along with functional tools, allow for testing these scenarios and discuss their predicted outcomes.

CRISPR/Cas9 and FLP-FRT mediated regulatory dissection of the BX-C of Drosophila melanogaster

The homeotic genes or Hox define the anterior-posterior (AP) body axis formation in bilaterians and are often present on the chromosome in an order collinear to their function across the AP axis. However, there are many cases wherein the Hox are not collinear, but their expression pattern is conserved across the AP axis. The expression pattern of Hox is attributed to the cis-regulatory modules (CRMs) consisting of enhancers, initiators, or repressor elements that regulate the genes in a segment-specific manner. In the Drosophila melanogaster Hox complex, the bithorax complex (BX-C) and even the CRMs are organized in an order that is collinear to their function in the thoracic and abdominal segments. In the present study, the regulatorily inert regions were targeted using CRISPR/Cas9 to generate a series of transgenic lines with the insertion of FRT sequences. These FRT lines are repurposed to shuffle the CRMs associated with Abd-B to generate modular deletion, duplication, or inversion of multiple CRMs. The rearrangements yielded entirely novel phenotypes in the fly suggesting the requirement of such complex manipulations to address the significance of higher order arrangement of the CRMs. The functional map and the transgenic flies generated in this study are important resources to decipher the collective ability of multiple regulatory elements in the eukaryotic genome to function as complex modules.

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.

Annotation of Hox cluster and Hox cofactor genes in the Asian citrus psyllid, Diaphorina citri, reveals novel features

Hox genes and their cofactors are essential developmental genes specifying regional identity in animals. Hox genes have a conserved arrangement in clusters in the same order in which they specify identity along the anterior-posterior axis. A few insect species have breaks in the cluster, but these are exceptions. We annotated the 10 Hox genes of the Asian citrus psyllid Diaphorina citri, and found a split in its Hox cluster between the Deformed and Sex combs reduced genes - the first time a break at this position has been observed in an insect Hox cluster. We also annotated D. citri orthologs of the Hox cofactor genes homothorax, PKNOX and extradenticle and found an additional copy of extradenticle in D. citri that appears to be a retrogene. Expression data and sequence conservation suggest that the extradenticle retrogene may have retained the original extradenticle function and allowed divergence of the parental extradenticle gene.

Phylogenetic analysis of forkhead transcription factors in the Panarthropoda

Fox genes encode transcription factors that contain a DNA binding domain, the forkhead domain, and are known from diverse animal species. The exact homology of the Fox genes of different species is debated and this makes inferences about the evolution of the Fox genes, and their duplications and losses difficult. We have performed phylogenetic analyses of the Fox gene complements of 32 panarthropod species. Our results confirm an ancestral complement of FoxA, FoxB, FoxC, FoxD, FoxF, FoxG, FoxJ1, FoxJ2/3, FoxK, FoxL1, FoxL2, FoxN1/4, FoxN2/3, FoxO, FoxP, and FoxQ2 in the Arthropoda, and additionally FoxH and FoxQ1 in the Panarthropoda (including tardigrades and onychophorans). We identify a novel Fox gene sub-family, that we designate as FoxT that includes two genes in Drosophila melanogaster, Circadianly Regulated Gene (Crg-1) and forkhead domain 3F (fd3F). In a very recent paper, the same new Fox gene sub-family was identified in insects (Lin et al. 2021). Our analysis confirms the presence of FoxT and shows that its members are present throughout Panarthropoda. We show that the hitherto unclassified gene CG32006 from the fly Drosophila melanogaster belongs to FoxJ1. We also detect gene losses: FoxE and FoxM were lost already in the panarthropod ancestor, whereas the loss of FoxH occurred in the arthropod ancestor. Finally, we find an ortholog of FoxQ1 in the bark scorpion Centruroides sculpturatus, confirmed not only by phylogenetic analysis, but also by forming an evolutionarily conserved gene cluster with FoxF, FoxC, and FoxL1. This suggests that FoxQ1 belongs to the ancestral Fox gene complement in panarthropods and also in chelicerates, but has been lost at the base of the mandibulate arthropods.

Phenotypic regionalization of the vertebral column in the thorny skate Amblyraja radiata: Stability and variation

Regionalization of the vertebral column occurred early during vertebrate evolution and has been extensively investigated in mammals. However, less data are available on vertebral regions of crown gnathostomes. This is particularly true for batoids (skates, sawfishes, guitarfishes, and rays) whose vertebral column has long been considered to be composed of the same two regions as in teleost fishes despite the presence of a synarcual. However, the numerous vertebral units in chondrichthyans may display a more complex regionalization pattern than previously assumed and the intraspecific variation of such pattern deserves a thorough investigation. In this study, we use micro-computed tomography (µCT) scans of vertebral columns of a growth series of thorny skates Amblyraja radiata to provide the first fine-scale morphological description of vertebral units in a batoids species. We further investigate axial regionalization using a replicable clustering analysis on presence/absence of vertebral elements to decipher the regionalization of the vertebral column of A. radiata. We identify four vertebral regions in this species. The two anteriormost regions, named synarcual and thoracic, may undergo strong developmental or functional constraints because they display stable patterns of shapes and numbers of vertebral units across all growth stages. The third region, named hemal transitional, is characterized by high inter-individual morphological variation and displays a transition between the monospondylous (one centrum per somite) to diplospondylous (two centra per somite) conditions. The posteriormost region, named caudal, is subdivided into three sub-regions with shapes changing gradually along the anteroposterior axis. These regionalized patterns are discussed in light of ecological habits of skates.

Improved Understanding of the Role of Gene and Genome Duplications in Chordate Evolution With New Genome and Transcriptome Sequences

Comparative approaches to understanding chordate genomes have uncovered a significant role for gene duplications, including whole genome duplications (WGDs), giving rise to and expanding gene families. In developmental biology, gene families created and expanded by both tandem and WGDs are paramount. These genes, often involved in transcription and signalling, are candidates for underpinning major evolutionary transitions because they are particularly prone to retention and subfunctionalisation, neofunctionalisation, or specialisation following duplication. Under the subfunctionalisation model, duplication lays the foundation for the diversification of paralogues, especially in the context of gene regulation. Tandemly duplicated paralogues reside in the same regulatory environment, which may constrain them and result in a gene cluster with closely linked but subtly different expression patterns and functions. Ohnologues (WGD paralogues) often diversify by partitioning their expression domains between retained paralogues, amidst the many changes in the genome during rediploidisation, including chromosomal rearrangements and extensive gene losses. The patterns of these retentions and losses are still not fully understood, nor is the full extent of the impact of gene duplication on chordate evolution. The growing number of sequencing projects, genomic resources, transcriptomics, and improvements to genome assemblies for diverse chordates from non-model and under-sampled lineages like the coelacanth, as well as key lineages, such as amphioxus and lamprey, has allowed more informative comparisons within developmental gene families as well as revealing the extent of conserved synteny across whole genomes. This influx of data provides the tools necessary for phylogenetically informed comparative genomics, which will bring us closer to understanding the evolution of chordate body plan diversity and the changes underpinning the origin and diversification of vertebrates.

Inferring Tunicate Relationships and the Evolution of the Tunicate Hox Cluster with the Genome of Corella inflata

Tunicates, the closest living relatives of vertebrates, have served as a foundational model of early embryonic development for decades. Comparative studies of tunicate phylogeny and genome evolution provide a critical framework for analyzing chordate diversification and the emergence of vertebrates. Toward this goal, we sequenced the genome of Corella inflata (Ascidiacea, Phlebobranchia), so named for the capacity to brood self-fertilized embryos in a modified, "inflated" atrial chamber. Combining the new genome sequence for Co. inflata with publicly available tunicate data, we estimated a tunicate species phylogeny, reconstructed the ancestral Hox gene cluster at important nodes in the tunicate tree, and compared patterns of gene loss between Co. inflata and Ciona robusta, the prevailing tunicate model species. Our maximum-likelihood and Bayesian trees estimated from a concatenated 210-gene matrix were largely concordant and showed that Aplousobranchia was nested within a paraphyletic Phlebobranchia. We demonstrated that this relationship is not an artifact due to compositional heterogeneity, as had been suggested by previous studies. In addition, within Thaliacea, we recovered Doliolida as sister to the clade containing Salpida and Pyrosomatida. The Co. inflata genome provides increased resolution of the ancestral Hox clusters of key tunicate nodes, therefore expanding our understanding of the evolution of this cluster and its potential impact on tunicate morphological diversity. Our analyses of other gene families revealed that several cardiovascular associated genes (e.g., BMP10, SCL2A12, and PDE2a) absent from Ci. robusta, are present in Co. inflata. Taken together, our results help clarify tunicate relationships and the genomic content of key ancestral nodes within this phylogeny, providing critical insights into tunicate evolution.

Hox gene expression during development of the phoronid Phoronopsis harmeri

Background

Phoronida is a small group of marine worm-like suspension feeders, which together with brachiopods and bryozoans form the clade Lophophorata. Although their development is well studied on the morphological level, data regarding gene expression during this process are scarce and restricted to the analysis of relatively few transcription factors. Here, we present a description of the expression patterns of Hox genes during the embryonic and larval development of the phoronid Phoronopsis harmeri.

Results

We identified sequences of eight Hox genes in the transcriptome of Ph. harmeri and determined their expression pattern during embryonic and larval development using whole mount in situ hybridization. We found that none of the Hox genes is expressed during embryonic development. Instead their expression is initiated in the later developmental stages, when the larval body is already formed. In the investigated initial larval stages the Hox genes are expressed in the non-collinear manner in the posterior body of the larvae: in the telotroch and the structures that represent rudiments of the adult worm. Additionally, we found that certain head-specific transcription factors are expressed in the oral hood, apical organ, preoral coelom, digestive system and developing larval tentacles, anterior to the Hox-expressing territories.

Conclusions

The lack of Hox gene expression during early development of Ph. harmeri indicates that the larval body develops without positional information from the Hox patterning system. Such phenomenon might be a consequence of the evolutionary intercalation of the larval form into an ancestral life cycle of phoronids. The observed Hox gene expression can also be a consequence of the actinotrocha representing a “head larva”, which is composed of the most anterior body region that is devoid of Hox gene expression. Such interpretation is further supported by the expression of head-specific transcription factors. This implies that the Hox patterning system is used for the positional information of the trunk rudiments and is, therefore, delayed to the later larval stages. We propose that a new body form was intercalated to the phoronid life cycle by precocious development of the anterior structures or by delayed development of the trunk rudiment in the ancestral phoronid larva.

De novo transcription of multiple Hox cluster genes takes place simultaneously in early Xenopus tropicalis embryos

hox genes are found as clusters in the genome in most bilaterians. The order of genes in the cluster is supposed to be correlated with the site of expression along the anterior-posterior body axis and the timing of expression during development, and these correlations are called spatial and temporal collinearity, respectively. Here we studied the expression dynamics of all hox genes of the diploid species Xenopus tropicalis in four Hox clusters (A–D) by analyzing high-temporal-resolution RNA-seq databases and the results showed that temporal collinearity is not supported, which is consistent with our previous data from allotetraploid Xenopuslaevis. Because the temporal collinearity hypothesis implicitly assumes the collinear order of gene activation, not mRNA accumulation, we determined for the first time the timing of when new transcripts of hox genes are produced, by detecting pre-spliced RNA in whole embryos with reverse transcription and quantitative PCR (RT-qPCR) for all hoxa genes as well as several selected hoxb, hoxc and hoxd genes. Our analyses showed that, coinciding with the RNA-seq results, hoxa genes started to be transcribed in a non-sequential order, and found that multiple genes start expression almost simultaneously or more posterior genes could be expressed earlier than anterior ones. This tendency was also found in hoxb and hoxc genes. These results suggest that temporal collinearity of hox genes is not held during early development of Xenopus.

Summary: qPCR analysis for de novo transcription of hox genes suggest that temporal collinearity is not held for all hox genes during early development of Xenopus tropicalis.

Vertebrate hox temporal collinearity: does it exist and what is it's function?

Hox temporal collinearity (TC) is a mysterious feature of embryogenesis. This article is opportune because of a recent challenge to TC’s existence This challenge is examined and the evidence that TC does exist is presented. Its function is discussed. Temporal collinearity is thought to be important because it lays the basis for Hox spatial collinearity and the vertebrate A-P axial pattern. The time-space translation mechanism whereby this occurs is examined.

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.

Staggered Hox expression is more widespread among molluscs than previously appreciated

Hox genes are expressed along the anterior-posterior body axis in a colinear fashion in the majority of bilaterians. Contrary to polyplacophorans, a group of aculiferan molluscs with conserved ancestral molluscan features, gastropods and cephalopods deviate from this pattern by expressing Hox genes in distinct morphological structures and not in a staggered fashion. Among conchiferans, scaphopods exhibit many similarities with gastropods, cephalopods and bivalves, however, the molecular developmental underpinnings of these similar traits remain unknown. We investigated Hox gene expression in developmental stages of the scaphopod Antalis entalis to elucidate whether these genes are involved in patterning morphological traits shared by their kin conchiferans. Scaphopod Hox genes are predominantly expressed in the foot and mantle but also in the central nervous system. Surprisingly, the scaphopod mid-stage trochophore exhibits a near-to staggered expression of all nine Hox genes identified. Temporal colinearity was not found and early-stage and late-stage trochophores, as well as postmetamorphic individuals, do not show any apparent traces of staggered expression. In these stages, Hox genes are expressed in distinct morphological structures such as the cerebral and pedal ganglia and in the shell field of early-stage trochophores. Interestingly, a re-evaluation of previously published data on early-stage cephalopod embryos and of the gastropod pre-torsional veliger shows that these developmental stages exhibit traces of staggered Hox expression. Considering our results and all gene expression and genomic data available for molluscs as well as other bilaterians, we suggest a last common molluscan ancestor with colinear Hox expression in predominantly ectodermal tissues along the anterior-posterior axis. Subsequently, certain Hox genes have been co-opted into the patterning process of distinct structures (apical organ or prototroch) in conchiferans.

Two Tier Hox Collinearity Mediates Vertebrate Axial Patterning

A two tier mechanism mediates Hox collinearity. Besides the familiar collinear chromatin modification within each Hox cluster (nanocollinearity), there is also a macrocollinearity tier. Individual Hox clusters and individual cells are coordinated and synchronized to generate multiscale (macro and nano) collinearity in the early vertebrate embryo. Macro-collinearity is mediated by three non-cell autonomous Hox–Hox interactions. These mediate temporal collinearity in early NOM (non-organizer mesoderm), time space translation where temporal collinearity is translated to spatial collinearity along the early embryo’s main body axis and neural transformation, where Hox expression is copied monospecifically from NOM mesoderm to overlying neurectoderm in the late gastrula. Unlike nanocollinearity, which is Hox cluster restricted, axial macrocollinearity extends into the head and EAD domains, thus covering the whole embryonic anterior-posterior (A-P) axis. EAD: extreme anterior domain, the only A-P axial domain anterior to the head. The whole time space translation mechanism interacts with A-P signaling pathways via “decision points,” separating different domains on the axis.

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.

Mechanisms of Specificity for Hox Factor Activity

Metazoans encode clusters of paralogous Hox genes that are critical for proper development of the body plan. However, there are a number of unresolved issues regarding how paralogous Hox factors achieve specificity to control distinct cell fates. First, how do Hox paralogs, which have very similar DNA binding preferences in vitro, drive different transcriptional programs in vivo? Second, the number of potential Hox binding sites within the genome is vast compared to the number of sites bound. Hence, what determines where in the genome Hox factors bind? Third, what determines whether a Hox factor will activate or repress a specific target gene? Here, we review the current evidence that is beginning to shed light onto these questions. In particular, we highlight how cooperative interactions with other transcription factors (especially PBC and HMP proteins) and the sequences of cis-regulatory modules provide a basis for the mechanisms of Hox specificity. We conclude by integrating a number of the concepts described throughout the review in a case study of a highly interrogated Drosophila cis-regulatory module named “The Distal-less Conserved Regulatory Element” (DCRE).

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.

TCF/Lef regulates the Gsx ParaHox gene in central nervous system development in chordates

BACKGROUND

The ParaHox genes play an integral role in the anterior-posterior (A-P) patterning of the nervous system and gut of most animals. The ParaHox cluster is an ideal system in which to study the evolution and regulation of developmental genes and gene clusters, as it displays similar regulatory phenomena to its sister cluster, the Hox cluster, but offers a much simpler system with only three genes.

RESULTS

Using Ciona intestinalis transgenics, we isolated a regulatory element upstream of Branchiostoma floridae Gsx that drives expression within the central nervous system of Ciona embryos. The minimal amphioxus enhancer region required to drive CNS expression has been identified, along with surrounding sequence that increases the efficiency of reporter expression throughout the Ciona CNS. TCF/Lef binding sites were identified and mutagenized and found to be required to drive the CNS expression. Also, individual contributions of TCF/Lef sites varied across the regulatory region, revealing a partial division of function across the Bf-Gsx-Up regulatory element. Finally, when all TCF/Lef binding sites are mutated CNS expression is not only abolished, but a latent repressive function is also unmasked.

CONCLUSIONS

We have identified a B. floridae Gsx upstream regulatory element that drives CNS expression within transgenic Ciona intestinalis, and have shown that this CNS expression is dependent upon TCF/Lef binding sites. We examine the evolutionary and developmental implications of these results, and discuss the possibility of TCF/Lef not only as a regulator of chordate Gsx, but as a deeply conserved regulatory factor controlling all three ParaHox genes across the Metazoa.

The origin of the Hox/ParaHox genes, the Ghost Locus hypothesis and the complexity of the first animal

A key aim in evolutionary biology is to deduce ancestral states to better understand the evolutionary origins of clades of interest and the diversification process(es) that has/have elaborated them. These ancestral deductions can hit difficulties when undetected loss events are misinterpreted as ancestral absences. With the ever-increasing amounts of animal genomic sequence data, we are gaining a much clearer view of the preponderance of differential gene losses across animal lineages. This has become particularly clear with recent progress in our understanding of the origins of the Hox/ParaHox developmental control genes relative to the earliest branching lineages of the animal kingdom: the sponges (Porifera), comb jellies (Ctenophora) and placozoans (Placozoa). These reassessments of the diversity and complexity of developmental control genes in the earliest animal ancestors need to go hand-in-hand with complementary advances in comparative morphology, phylogenetics and palaeontology to clarify our understanding of the complexity of the last common ancestor of all animals. The field is currently undergoing a shift from the traditional consensus of a sponge-like animal ancestor from which morphological and molecular elaboration subsequently evolved, to a scenario of a more complex animal ancestor, with subsequent losses and simplifications in various lineages.

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.

Time space translation: a hox mechanism for vertebrate a-p patterning

The vertebrate A-P axis is a time axis. The head is made first and more and more posterior levels are made at later and later stages. This is different to the situation in most other animals, for example, in Drosophila. Central to this timing is Hox temporal collinearity (see below). This occurs rarely in the animal kingdom but is characteristic of vertebrates and is used to generate the primary axial Hox pattern using time space translation and to integrate successive derived patterns (see below). This is thus a different situation than in Drosophila, where the primary pattern guiding Hox spatial collinearity is generated externally, by the gap and segmentation genes.

Coral comparative genomics reveal expanded Hox cluster in the cnidarian-bilaterian ancestor

The key developmental role of the Hox cluster of genes was established prior to the last common ancestor of protostomes and deuterostomes and the subsequent evolution of this cluster has played a major role in the morphological diversity exhibited in extant bilaterians. Despite 20 years of research into cnidarian Hox genes, the nature of the cnidarian-bilaterian ancestral Hox cluster remains unclear. In an attempt to further elucidate this critical phylogenetic node, we have characterized the Hox cluster of the recently sequenced Acropora digitifera genome. The A. digitifera genome contains two anterior Hox genes (PG1 and PG2) linked to an Eve homeobox gene and an Anthox1A gene, which is thought to be either a posterior or posterior/central Hox gene. These data show that the Hox cluster of the cnidarian-bilaterian ancestor was more extensive than previously thought. The results are congruent with the existence of an ancient set of constraints on the Hox cluster and reinforce the importance of incorporating a wide range of animal species to reconstruct critical ancestral nodes.

Identical genomic organization of two hemichordate hox clusters

Genomic comparisons of chordates, hemichordates, and echinoderms can inform hypotheses for the evolution of these strikingly different phyla from the last common deuterostome ancestor. Because hox genes play pivotal developmental roles in bilaterian animals, we analyzed the Hox complexes of two hemichordate genomes. We find that Saccoglossus kowalevskii and Ptychodera flava both possess 12-gene clusters, with mir10 between hox4 and hox5, in 550 kb and 452 kb intervals, respectively. Genes hox1-hox9/10 of the clusters are in the same genomic order and transcriptional orientation as their orthologs in chordates, with hox1 at the 3' end of the cluster. At the 5' end, each cluster contains three posterior genes specific to Ambulacraria (the hemichordate-echinoderm clade), two forming an inverted terminal pair. In contrast, the echinoderm Strongylocentrotus purpuratus contains a 588 kb cluster of 11 orthologs of the hemichordate genes, ordered differently, plausibly reflecting rearrangements of an ancestral hemichordate-like ambulacrarian cluster. Hox clusters of vertebrates and the basal chordate amphioxus have similar organization to the hemichordate cluster, but with different posterior genes. These results provide genomic evidence for a well-ordered complex in the deuterostome ancestor for the hox1-hox9/10 region, with the number and kind of posterior genes still to be elucidated.

Global posterior prevalence is unique to vertebrates: a dance to the music of time?

We reach the conclusion that posterior prevalence, a collinear property considered important for Hox complex function, is so far unique, in a global form, to vertebrates. Why is this? We suspect this is because posterior prevalence is explicitly connected to the vertebrate form of Hox temporal collinearity, which is central to axial patterning.

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.

A genetic approach to the transcriptional regulation of Hox gene clusters

The evolution of vertebrate genomes was accompanied by an astounding increase in the complexity of their regulatory modalities. Genetic redundancy resulting from large-scale genome duplications at the base of the chordate tree was repeatedly exploited by the functional redeployment of paralogous genes via innovations in their regulatory circuits. As a paradigm of such regulatory evolution, we have extensively studied those control mechanisms at work in-cis over vertebrate Hox gene clusters. Here, we review the portfolio of genetic strategies that have been developed to tackle the intricate relationship between genomic topography and the transcriptional activities in this gene family, and we describe some of the mechanistic insights we gained by using the HoxD cluster as an example. We discuss the high heuristic value of this system in our general understanding of how changes in transcriptional regulation can diversify gene function and thereby fuel morphological 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.

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.

Evolution of Hox gene clusters in gnathostomes: insights from a survey of a shark (Scyliorhinus canicula) transcriptome

It is now well established that there were four Hox gene clusters in the genome of the last common ancestor of extant gnathostomes. To better understand the evolution of the organization and expression of these genomic regions, we have studied the Hox gene clusters of a shark (Scyliorhinus canicula). We sequenced 225,580 expressed sequence tags from several embryonic cDNA libraries. Blast searches identified corresponding transcripts to almost all the HoxA, HoxB, and HoxD cluster genes. No HoxC transcript was identified, suggesting that this cluster is absent or highly degenerate. Using Hox gene sequences as probes, we selected and sequenced seven clones from a bacterial artificial chromosome library covering the complete region of the three gene clusters. Mapping of cDNAs to these genomic sequences showed extensive alternative splicing and untranslated exon sharing between neighboring Hox genes. Homologous noncoding exons could not be identified in transcripts from other species using sequence similarity. However, by comparing conserved noncoding sequences upstream of these exons in different species, we were able to identify homology between some exons. Some alternative splicing variants are probably very ancient and were already coded for by the ancestral Hox gene cluster. We also identified several transcripts that do not code for Hox proteins, are probably not translated, and all but one are in the reverse orientation to the Hox genes. This survey of the transcriptome of the Hox gene clusters of a shark shows that the high complexity observed in mammals is a gnathostome ancestral feature.

Improving Hox protein classification across the major model organisms

The family of Hox-proteins has been a major focus of research for over 30 years. Hox-proteins are crucial to the correct development of bilateral organisms, however, some uncertainty remains as to which Hox-proteins are functionally equivalent across different species. Initial classification of Hox-proteins was based on phylogenetic analysis of the 60 amino acid homeodomain. This approach was successful in classifying Hox-proteins with differing homeodomains, but the relationships of Hox-proteins with nearly identical homeodomains, yet distinct biological functions, could not be resolved. Correspondingly, these ‘problematic’ proteins were classified into one large unresolved group. Other classifications used the relative location of the Hox-protein coding genes on the chromosome (synteny) to further resolve this group. Although widely used, this synteny-based classification is inconsistent with experimental evidence from functional equivalence studies. These inconsistencies led us to re-examine and derive a new classification for the Hox-protein family using all Hox-protein sequences available in the GenBank non-redundant protein database (NCBI-nr). We compare the use of the homeodomain, the homeodomain with conserved flanking regions (the YPWM and linker region), and full length Hox-protein sequences as a basis for classification of Hox-proteins. In contrast to previous attempts, our approach is able to resolve the relationships for the ‘problematic’ as well as ABD-B-like Hox-proteins. We highlight differences to previous classifications and clarify the relationships of Hox-proteins across the five major model organisms, Caenorhabditis elegans, Drosophila melanogaster, Branchiostoma floridae, Mus musculus and Danio rerio. Comparative and functional analysis of Hox-proteins, two fields crucial to understanding the development of bilateral organisms, have been hampered by difficulties in predicting functionally equivalent Hox-proteins across species. Our classification scheme offers a higher-resolution classification that is in accordance with phylogenetic as well as experimental data and, thereby, provides a novel basis for experiments, such as comparative and functional analyses of Hox-proteins.

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.

Indirect development, transdifferentiation and the macroregulatory evolution of metazoans

It is proposed here that a biphasic life cycle with partial dedifferentiation of intermediate juvenile or larval stages represents the mainstream developmental mode of metazoans. Developmental plasticity of differentiated cells is considered the essential characteristic of indirect development, rather than the exclusive development of the adult from 'set-aside' cells. Many differentiated larval cells of indirect developers resume proliferation, partially dedifferentiate and contribute to adult tissues. Transcriptional pluripotency of differentiated states has premetazoan origins and seems to be facilitated by histone variant H2A.Z. Developmental plasticity of differentiated states also facilitates the evolution of polyphenism. Uncertainty remains about whether the most recent common ancestor of protostomes and deuterostomes was a direct or an indirect developer, and how the feeding larvae of bilaterians are related to non-feeding larvae of sponges and cnidarians. Feeding ciliated larvae of bilaterians form their primary gut opening by invagination, which seems related to invagination in cnidarians. Formation of the secondary gut opening proceeds by protostomy or deuterostomy, and gene usage suggests serial homology of the mouth and anus. Indirect developers do not use the Hox vector to build their ciliated larvae, but the Hox vector is associated with the construction of the reproductive portion of the animal during feeding-dependent posterior growth. It is further proposed that the original function of the Hox cluster was in gonad formation rather than in anteroposterior diversification.

Cis-regulatory characterization of sequence conservation surrounding the Hox4 genes

Hox genes are key regulators of anterior-posterior axis patterning and have a major role in hindbrain development. The zebrafish Hox4 paralogs have strong overlapping activities in hindbrain rhombomeres 7 and 8, in the spinal cord and in the pharyngeal arches. With the aim to predict enhancers that act on the hoxa4a, hoxb4a, hoxc4a and hoxd4a genes, we used sequence conservation around the Hox4 genes to analyze all fish:human conserved non-coding sequences by reporter assays in stable zebrafish transgenesis. Thirty-four elements were functionally tested in GFP reporter gene constructs and more than 100 F1 lines were analyzed to establish a correlation between sequence conservation and cis-regulatory function, constituting a catalog of Hox4 CNEs. Sixteen tissue-specific enhancers could be identified. Multiple alignments of the CNEs revealed paralogous cis-regulatory sequences, however, the CNE sequence similarities were found not to correlate with tissue specificity. To identify ancestral enhancers that direct Hox4 gene activity, genome sequence alignments of mammals, teleosts, horn shark and the cephalochordate amphioxus, which is the most basal extant chordate possessing a single prototypical Hox cluster, were performed. Three elements were identified and two of them exhibited regulatory activity in transgenic zebrafish, however revealing no specificity. Our data show that the approach to identify cis-regulatory sequences by genome sequence alignments and subsequent testing in zebrafish transgenesis can be used to define enhancers within the Hox clusters and that these have significantly diverged in their function during evolution.

Clustered Fox genes in lophotrochozoans and the evolution of the bilaterian Fox gene cluster

FoxC, FoxF, FoxL1 and FoxQ1 genes have been shown to be clustered in some animal genomes, with mesendodermal expression hypothesised as a selective force maintaining cluster integrity. Hypotheses are, however, constrained by a lack of data from the Lophotrochozoa. Here we characterise members of the FoxC, FoxF, FoxL1 and FoxQ1 families from the annelid Capitella teleta and the molluscs Lottia gigantea and Patella vulgata. We cloned FoxC, FoxF, FoxL1 and FoxQ1 genes from C. teleta, and FoxC, FoxF and FoxL1 genes from P. vulgata, and established their expression during development. We also examined their genomic organisation in C. teleta and L. gigantea, and investigated local syntenic relationships. Our results show mesodermal and anterior gut expression is a common feature of these genes in lophotrochozoans. In L. gigantea FoxC, FoxF and FoxL1 are closely linked, while in C. teleta Ct-foxC and Ct-foxL1 are closely linked, with Ct-foxF and Ct-foxQ1 on different scaffolds. Adjacent to these genes there is limited evidence of local synteny. This demonstrates conservation of genomic organisation and expression of these genes can be traced in all three bilaterian Superphyla. These data are evaluated against competing theories for the long-term maintenance of gene clusters.

Coordinated spatial and temporal expression of Hox genes during embryogenesis in the acoel Convolutriloba longifissura

Background

Hox genes are critical for patterning the bilaterian anterior-posterior axis. The evolution of their clustered genomic arrangement and ancestral function has been debated since their discovery. As acoels appear to represent the sister group to the remaining Bilateria (Nephrozoa), investigating Hox gene expression will provide an insight into the ancestral features of the Hox genes in metazoan evolution.

Results

We describe the expression of anterior, central and posterior class Hox genes and the ParaHox ortholog Cdx in the acoel Convolutriloba longifissura. Expression of all three Hox genes begins contemporaneously after gastrulation and then resolves into staggered domains along the anterior-posterior axis, suggesting that the spatial coordination of Hox gene expression was present in the bilaterian ancestor. After early surface ectodermal expression, the anterior and central class genes are expressed in small domains of putative neural precursor cells co-expressing ClSoxB1, suggesting an evolutionary early function of Hox genes in patterning parts of the nervous system. In contrast, the expression of the posterior Hox gene is found in all three germ layers in a much broader posterior region of the embryo.

Conclusion

Our results suggest that the ancestral set of Hox genes was involved in the anterior-posterior patterning of the nervous system of the last common bilaterian ancestor and were later co-opted for patterning in diverse tissues in the bilaterian radiation. The lack of temporal colinearity of Hox expression in acoels may be due to a loss of genomic clustering in this clade or, alternatively, temporal colinearity may have arisen in conjunction with the expansion of the Hox cluster in the Nephrozoa.

Features of the ancestral bilaterian inferred from Platynereis dumerilii ParaHox genes

Background

The ParaHox gene cluster is the evolutionary sister to the Hox cluster. Whilst the role of the Hox cluster in patterning the anterior-posterior axis of bilaterian animals is well established, and the organisation of vertebrate Hox clusters is intimately linked to gene regulation, much less is known about the more recently discovered ParaHox cluster. ParaHox gene clustering, and its relationship to expression, has only been described in deuterostomes. Conventional protostome models (Drosophila melanogaster and Caenorhabditis elegans) are secondarily derived with respect to ParaHox genes, suffering gene loss and cluster break-up.

Results

We provide the first evidence for ParaHox gene clustering from a less-derived protostome animal, the annelid Platynereis dumerilii. Clustering of these genes is thus not a sole preserve of the deuterostome lineage within Bilateria. This protostome ParaHox cluster is not entirely intact however, with Pdu-Cdx being on the opposite end of the same chromosome arm from Pdu-Gsx and Pdu-Xlox. From the genomic sequence around the P. dumerilii ParaHox genes the neighbouring genes are identified, compared with other taxa, and the ancestral arrangement deduced.

Conclusion

We relate the organisation of the ParaHox genes to their expression, and from comparisons with other taxa hypothesise that a relatively complex pattern of ParaHox gene expression existed in the protostome-deuterostome ancestor, which was secondarily simplified along several invertebrate lineages. Detailed comparisons of the gene content around the ParaHox genes enables the reconstruction of the genome surrounding the ParaHox cluster of the protostome-deuterostome ancestor, which existed over 550 million years ago.

Uncoupling time and space in the collinear regulation of Hox genes

During development of the vertebrate body axis, Hox genes are transcribed sequentially, in both time and space, following their relative positions within their genomic clusters. Analyses of animal genomes support the idea that Hox gene clustering is essential for coordinating the various times of gene activations. However, the eventual collinear ordering of the gene specific transcript domains in space does not always require genomic clustering. We analyzed these complex regulatory relationships by using mutant alleles at the mouse HoxD locus, including one that splits the cluster into two pieces. We show that both positive and negative regulatory influences, located on either side of the cluster, control an early phase of collinear expression in the trunk. Interestingly, this early phase does not systematically impact upon the subsequent expression patterns along the main body axis, indicating that the mechanism underlying temporal collinearity is distinct from those acting during the second phase. We discuss the potential functions and evolutionary origins of these mechanisms, as well as their relationship with similar processes at work during limb development.

Hox genes encode proteins that control embryonic development along the head-to-tail axis. These genes are clustered in one site on the chromosome and their respective positions within the cluster determine their time and place of activation. Here, by using a large set of targeted mutations disturbing the integrity of the gene cluster, we show that the spatial organization of expression domains does not directly depend upon the timing of activation as was previously suggested. This uncoupling between space and time in the regulation of these Hox genes coincides with the existence of two major phases of regulation. The first is time-dependent and involves global regulatory influences, located outside the gene cluster, whereas the second relies upon more local regulatory elements, likely interspersed between the genes, inside the cluster. These results provide the bases for future analyses of collinear mechanisms and indicate that different types of collinearities are not necessarily related, neither in function, nor in their evolutionary histories.

Hox cluster duplication in the basal teleost Hiodon alosoides (Osteoglossomorpha)

Large-scale—even genome-wide—duplications have repeatedly been invoked as an explanation for major radiations. Teleosts, the most species-rich vertebrate clade, underwent a “fish-specific genome duplication” (FSGD) that is shared by most ray-finned fish lineages. We investigate here the Hox complement of the goldeye (Hiodon alosoides), a representative of Osteoglossomorpha, the most basal teleostean clade. An extensive PCR survey reveals that goldeye has at least eight Hox clusters, indicating a duplicated genome compared to basal actinopterygians. The possession of duplicated Hox clusters is uncoupled to species richness. The Hox system of the goldeye is substantially different from that of other teleost lineages, having retained several duplicates of Hox genes for which crown teleosts have lost at least one copy. A detailed analysis of the PCR fragments as well as full length sequences of two HoxA13 paralogs, and HoxA10 and HoxC4 genes places the duplication event close in time to the divergence of Osteoglossomorpha and crown teleosts. The data are consistent with—but do not conclusively prove—that Osteoglossomorpha shares the FSGD.