Hox cluster duplications and the opportunity for evolutionary novelties

Hox genes play a key role in animal body plan development. These genes tend to occur in tightly linked clusters in the genome. Vertebrates and invertebrates differ in their Hox cluster number, with vertebrates having multiple clusters and invertebrates usually having only one. Recent evidence shows that vertebrate Hox clusters are structurally more constrained than invertebrate Hox clusters; they exclude transposable elements, do not undergo tandem duplications, and conserve their intergenic distances and gene order. These constraints are only relaxed after a cluster duplication. In contrast, invertebrate Hox clusters are structurally more plastic; tandem duplications are common, the linkage of Hox genes can change quickly, or they can lose their structural integrity completely. We propose that the constraints on vertebrate Hox cluster structure lead to an association between the retention of duplicated Hox clusters and adaptive radiations. After a duplication the constraints on Hox cluster structure are temporarily lifted, which opens a window of evolvability for the Hox clusters. If this window of evolvability coincides with an adaptive radiation, chances are that a modified Hox cluster becomes recruited in an evolutionary novelty and then both copies of duplicated Hox clusters are retained.     

Elephant shark (Callorhinchus milii) provides insights into the evolution of Hox gene clusters in gnathostomes

We have sequenced and analyzed Hox gene clusters from elephant shark, a holocephalian cartilaginous fish. Elephant shark possesses 4 Hox clusters with 45 Hox genes that include orthologs for a higher number of ancient gnathostome Hox genes than the 4 clusters in tetrapods and the supernumerary clusters in teleost fishes. Phylogenetic analysis of elephant shark Hox genes from 7 paralogous groups that contain all of the 4 members indicated an ((AB)(CD)) topology for the order of Hox cluster duplication, providing support for the 2R hypothesis (i.e., 2 rounds of whole-genome duplication during the early evolution of vertebrates). Comparisons of noncoding sequences of the elephant shark and human Hox clusters have identified a large number of conserved noncoding elements (CNEs), which represent putative cis-regulatory elements that may be involved in the regulation of Hox genes. Interestingly, in fugu more than 50% of these ancient CNEs have diverged beyond recognition in the duplicated (HoxA, HoxB, and HoxD) as well as the singleton (HoxC) Hox clusters. Furthermore, the b-paralogs of the duplicated fugu Hox clusters are virtually devoid of unique ancient CNEs. In contrast to fugu Hox clusters, elephant shark and human Hox clusters have lost fewer ancient CNEs. If these ancient CNEs are indeed enhancers directing tissue-specific expression of Hox genes, divergence of their sequences in vertebrate lineages might have led to altered expression patterns and presumably the functions of their associated Hox genes.     

Hox cluster genomics in the horn shark, Heterodontus francisci

Reconstructing the evolutionary history of Hox cluster origins will lead to insights into the developmental and evolutionary significance of Hox gene clusters in vertebrate phylogeny and to their role in the origins of various vertebrate body plans. We have isolated two Hox clusters from the horn shark, Heterodontus francisci. These have been sequenced and compared with one another and with other chordate Hox clusters. The results show that one of the horn shark clusters (HoxM) is orthologous to the mammalian HoxA cluster and shows a structural similarity to the amphioxus cluster, whereas the other shark cluster (HoxN) is orthologous to the mammalian HoxD cluster based on cluster organization and a comparison with noncoding and Hox gene-coding sequences. The persistence of an identifiable HoxA cluster over an 800-million-year divergence time demonstrates that the Hox gene clusters are highly integrated and structured genetic entities. The data presented herein identify many noncoding sequence motifs conserved over 800 million years that may function as genetic control motifs essential to the developmental process.     

Noncanonical role of Hox14 revealed by its expression patterns in lamprey and shark

Hox genes are arranged in uninterrupted clusters in vertebrate genomes, and the nested patterns of their expression define spatial identities in multiple embryonic tissues. The ancestral Hox cluster of vertebrates has long been thought to consist of, maximally, 13 Hox genes. However, recently, Hox14 genes were discovered in three chordate lineages, the coelacanth, cartilaginous fishes, and amphioxus, but their expression patterns have not yet been analyzed. We isolated Hox14 cDNAs from the Japanese lamprey and cloudy catshark. These genes were not expressed in the central nervous systems, somites, or fin buds/folds but were expressed in a restricted cell population surrounding the hindgut. The lack of Hox14 expression in most of the embryonic axial elements, where nested Hox expressions define spatial identities, suggests a decoupling of Hox14 genes' regulation from the ancestral regulatory mechanism. The relaxation of preexisting constraint for collinear expression may have permitted the secondary losses of this Hox member in the tetrapod and teleost lineages.     

Molecular evolution of the HoxA cluster in the three major gnathostome lineages

The duplication of Hox clusters and their maintenance in a lineage has a prominent but little understood role in chordate evolution. Here we examined how Hox cluster duplication may influence changes in cluster architecture and patterns of noncoding sequence evolution. We sequenced the entire duplicated HoxAa and HoxAb clusters of zebrafish (Danio rerio) and extended the 5' (posterior) part of the HoxM (HoxA-like) cluster of horn shark (Heterodontus francisci) containing the hoxa11 and hoxa13 orthologs as well as intergenic and flanking noncoding sequences. The duplicated HoxA clusters in zebrafish each house considerably fewer genes and are dramatically shorter than the single HoxA clusters of human and horn shark. We compared the intergenic sequences of the HoxA clusters of human, horn shark, zebrafish (Aa, Ab), and striped bass and found extensive conservation of noncoding sequence motifs, i.e., phylogenetic footprints, between the human and horn shark, representing two of the three gnathostome lineages. These are putative cis-regulatory elements that may play a role in the regulation of the ancestral HoxA cluster. In contrast, homologous regions of the duplicated HoxAa and HoxAb clusters of zebrafish and the HoxA cluster of striped bass revealed a striking loss of conservation of these putative cis-regulatory sequences in the 3' (anterior) segment of the cluster, where zebrafish only retains single representatives of group 1, 3, 4, and 5 (HoxAa) and group 2 (HoxAb) genes and in the 5' part of the clusters, where zebrafish retains two copies of the group 13, 11, and 9 genes, i.e., AbdB-like genes. In analyzing patterns of cis-sequence evolution in the 5' part of the clusters, we explicitly looked for evidence of complementary loss of conserved noncoding sequences, as predicted by the duplication-degeneration-complementation model in which genetic redundancy after gene duplication is resolved because of the fixation of complementary degenerative mutations. Our data did not yield evidence supporting this prediction. We conclude that changes in the pattern of cis-sequence conservation after Hox cluster duplication are more consistent with being the outcome of adaptive modification rather than passive mechanisms that erode redundancy created by the duplication event. These results support the view that genome duplications may provide a mechanism whereby master control genes undergo radical modifications conducive to major alterations in body plan. Such genomic revolutions may contribute significantly to the evolutionary process.     

Class 3 Hox genes in insects and the origin of zen.

We have cloned, from a beetle and a locust, genes that are homologous to the class 3 Hox genes of vertebrates. Outside the homeobox they share sequence motifs with the Drosophila zerknüllt (zen) and z2 genes, and like zen, are expressed only in extraembryonic membranes. We conclude that the zen genes of Drosophila derive from a Hox class 3 sequence that formed part of the common ancestral Hox cluster, but that in insects this (Hox) gene has lost its role in patterning the anterio-posterior axis of the embryo, and acquired a new function. In the lineage leading to Drosophila, the zen genes have diverged particularly rapidly.     

Highly conserved syntenic blocks at the vertebrate Hox loci and conserved regulatory elements within and outside Hox gene clusters

Hox genes in vertebrates are clustered, and the organization of the clusters has been highly conserved during evolution. The conservation of Hox clusters has been attributed to enhancers located within and outside the Hox clusters that are essential for the coordinated "temporal" and "spatial" expression patterns of Hox genes in developing embryos. To identify evolutionarily conserved regulatory elements within and outside the Hox clusters, we obtained contiguous sequences for the conserved syntenic blocks from the seven Hox loci in fugu and carried out a systematic search for conserved noncoding sequences (CNS) in the human, mouse, and fugu Hox loci. Our analysis has uncovered unusually large conserved syntenic blocks at the HoxA and HoxD loci. The conserved syntenic blocks at the human and mouse HoxA and HoxD loci span 5.4 Mb and 4 Mb and contain 21 and 19 genes, respectively. The corresponding regions in fugu are 16- and 12-fold smaller. A large number of CNS was identified within the Hox clusters and outside the Hox clusters spread over large regions. The CNS include previously characterized enhancers and overlap with the 5' global control regions of HoxA and HoxD clusters. Most of the CNS are likely to be control regions involved in the regulation of Hox and other genes in these loci. We propose that the regulatory elements spread across large regions on either side of Hox clusters are a major evolutionary constraint that has maintained the exceptionally long syntenic blocks at the HoxA and HoxD loci.     

From the Cover: Role of a polymorphism in a Hox/Pax-responsive enhancer in the evolution of the vertebrate spine

Patterning of the vertebrate skeleton requires the coordinated activity of Hox genes. In particular, Hox10 proteins are essential to set the transition from thoracic to lumbar vertebrae because of their rib-repressing activity. In snakes, however, the thoracic region extends well into Hox10-expressing areas of the embryo, suggesting that these proteins are unable to block rib formation. Here, we show that this is not a result of the loss of rib-repressing properties by the snake proteins, but rather to a single base pair change in a Hox/Paired box (Pax)-responsive enhancer, which prevents the binding of Hox proteins. This polymorphism is also found in Paenungulata, such as elephants and manatees, which have extended rib cages. In vivo, this modified enhancer failed to respond to Hox10 activity, supporting its role in the extension of rib cages. In contrast, the enhancer could still interact with Hoxb6 and Pax3 to promote rib formation. These results suggest that a polymorphism in the Hox/Pax-responsive enhancer may have played a role in the evolution of the vertebrate spine by differently modulating its response to rib-suppressing and rib-promoting Hox proteins.     

A long-range regulatory element of Hoxc8 identified by using the pClasper vector.

Hox genes are located in highly conserved clusters. The significance of this organization is unclear, but one possibility is that regulatory regions for individual genes are dispersed throughout the cluster and shared with other Hox genes. This hypothesis is supported by studies on several Hox genes in which even large genomic regions immediately surrounding the gene fail to direct the complete expression pattern in transgenic mice. In particular, previous studies have identified proximal regulatory regions that are primarily responsible for early phases of mouse Hoxc8 expression. To locate additional regulatory regions governing expression during the later periods of development, a yeast homologous recombination-based strategy utilizing the pClasper vector was employed. Using homologous recombination into pClasper, we cloned a 27-kb region around the Hoxc8 gene from a yeast artificial chromosome. A reporter gene was introduced into the coding region of the isolated gene by homologous recombination in yeast. This large fragment recapitulates critical aspects of Hoxc8 expression in transgenic mice. We show that the regulatory elements that maintain the anterior boundaries of expression in the neural tube and paraxial mesoderm are located between 11 and 19 kb downstream of the gene.     

Two steps in the evolution of Antennapedia-class vertebrate homeobox genes

Antennapedia-class vertebrate homeobox genes have been classified with regard to their chromosomal locations and nucleotide sequence similarities within the 183-base-pair homeobox domain. The results of these comparisons support the view that in mammals and most likely the vertebrates, four clusters of homeobox genes exist that were created by duplications of an entire primordial gene cluster. We present evidence that this primordial cluster arose by local gene duplications of homeoboxes that were present before the divergence of arthropods and chordates. Sequence analyses indicate that the expansion of the primordial gene cluster complex was accompanied by diversification, whereas conservation predominated after the duplications of entire homeobox gene clusters.     

Dispersal of NK homeobox gene clusters in amphioxus and humans

The Drosophila melanogaster genome has six physically clustered NK-related homeobox genes in just 180 kb. Here we show that the NK homeobox gene cluster was an ancient feature of bilaterian animal genomes, but has been secondarily split in chordate ancestry. The NK homeobox gene clusters of amphioxus and vertebrates are each split and dispersed at two equivalent intergenic positions. From the ancestral NK gene cluster, only the Tlx-Lbx and NK3-NK4 linkages have been retained in chordates. This evolutionary pattern is in marked contrast to the Hox and ParaHox gene clusters, which are compact in amphioxus and vertebrates, but have been disrupted in Drosophila.     

Elephant shark sequence reveals unique insights into the evolutionary history of vertebrate genes: A comparative analysis of the protocadherin cluster

Cartilaginous fishes are the oldest living phylogenetic group of jawed vertebrates. Here, we demonstrate the value of cartilaginous fish sequences in reconstructing the evolutionary history of vertebrate genomes by sequencing the protocadherin cluster in the relatively small genome (910 Mb) of the elephant shark (Callorhinchus milii). Human and coelacanth contain a single protocadherin cluster with 53 and 49 genes, respectively, that are organized in three subclusters, Pcdhalpha, Pcdhbeta, and Pcdhgamma, whereas the duplicated protocadherin clusters in fugu and zebrafish contain >77 and 107 genes, respectively, that are organized in Pcdhalpha and Pcdhgamma subclusters. By contrast, the elephant shark contains a single protocadherin cluster with 47 genes organized in four subclusters (Pcdhdelta, Pcdhepsilon, Pcdhmu, and Pcdhnu). By comparison with elephant shark sequences, we discovered a Pcdhdelta subcluster in teleost fishes, coelacanth, Xenopus, and chicken. Our results suggest that the protocadherin cluster in the ancestral jawed vertebrate contained more subclusters than modern vertebrates, and the evolution of the protocadherin cluster is characterized by lineage-specific differential loss of entire subclusters of genes. In contrast to teleost fish and mammalian protocadherin genes that have undergone gene conversion events, elephant shark protocadherin genes have experienced very little gene conversion. The syntenic block of genes in the elephant shark protocadherin locus is well conserved in human but disrupted in fugu. Thus, the elephant shark genome appears to be less prone to rearrangements compared with teleost fish genomes. The small and "stable" genome of the elephant shark is a valuable reference for understanding the evolution of vertebrate genomes.     

Breakup of a homeobox cluster after genome duplication in teleosts

Several families of homeobox genes are arranged in genomic clusters in metazoan genomes, including the Hox, ParaHox, NK, Rhox, and Iroquois gene clusters. The selective pressures responsible for maintenance of these gene clusters are poorly understood. The ParaHox gene cluster is evolutionarily conserved between amphioxus and human but is fragmented in teleost fishes. We show that two basal ray-finned fish, Polypterus and Amia, each possess an intact ParaHox cluster; this implies that the selective pressure maintaining clustering was lost after whole-genome duplication in teleosts. Cluster breakup is because of gene loss, not transposition or inversion, and the total number of ParaHox genes is the same in teleosts, human, mouse, and frog. We propose that this homeobox gene cluster is held together in chordates by the existence of interdigitated control regions that could be separated after locus duplication in the teleost fish.     

Ciona intestinalis Hox gene cluster: Its dispersed structure and residual colinear expression in development

Ascidians, belonging to the subphylum Urochordata, the earliest branch from the lineage to the vertebrates, exhibit a prototypical morphogenesis of chordates in the larval development, although they subsequently metamorphose into adults with a unique body structure. Recent draft genome analysis of the ascidian Ciona intestinalis has identified 9 Hox genes, which, however, have been located on five scaffolds. Similarly, expression patterns of Ciona Hox genes are largely unknown, although some data have been available for a few Hox member genes. Thus, the cluster structure and colinearity of Hox genes are still an enigma in C. intestinalis. To address these issues, we used fluorescence in situ hybridization and whole-mount in situ hybridization techniques and examined the genomic organization and spatiotemporal expression of all Hox as well as extended Hox member genes (Evx and Mox) of C. intestinalis. We found that seven of nine Ciona Hox genes are located on a single chromosome with some ordering exceptions, and the other genes, including Evx and Mox, are located on three other chromosomes. Some Ciona Hox genes, if not all, exhibited spatially coordinated expression within the larval central nervous system and the gut of the juvenile. In light of these observations, we suggest that the cluster organization and colinearity of the Hox genes are under dispersing and disintegrating conditions in C. intestinalis.