The amphioxusHox cluster: characterization, comparative genomics, and evolution

Chromosome-level genome assemblies of 2 hemichordates provide new insights into deuterostome origin and chromosome evolution

Deuterostomes are a monophyletic group of animals that includes Hemichordata, Echinodermata (together called Ambulacraria), and Chordata. The diversity of deuterostome body plans has made it challenging to reconstruct their ancestral condition and to decipher the genetic changes that drove the diversification of deuterostome lineages. Here, we generate chromosome-level genome assemblies of 2 hemichordate species, Ptychodera flava and Schizocardium californicum, and use comparative genomic approaches to infer the chromosomal architecture of the deuterostome common ancestor and delineate lineage-specific chromosomal modifications. We show that hemichordate chromosomes (1N = 23) exhibit remarkable chromosome-scale macrosynteny when compared to other deuterostomes and can be derived from 24 deuterostome ancestral linkage groups (ALGs). These deuterostome ALGs in turn match previously inferred bilaterian ALGs, consistent with a relatively short transition from the last common bilaterian ancestor to the origin of deuterostomes. Based on this deuterostome ALG complement, we deduced chromosomal rearrangement events that occurred in different lineages. For example, a fusion-with-mixing event produced an Ambulacraria-specific ALG that subsequently split into 2 chromosomes in extant hemichordates, while this homologous ALG further fused with another chromosome in sea urchins. Orthologous genes distributed in these rearranged chromosomes are enriched for functions in various developmental processes. We found that the deeply conserved Hox clusters are located in highly rearranged chromosomes and that maintenance of the clusters are likely due to lower densities of transposable elements within the clusters. We also provide evidence that the deuterostome-specific pharyngeal gene cluster was established via the combination of 3 pre-assembled microsyntenic blocks. We suggest that since chromosomal rearrangement events and formation of new gene clusters may change the regulatory controls of developmental genes, these events may have contributed to the evolution of diverse body plans among deuterostomes.

First Genome of Rock Lizard Darevskia valentini Involved in Formation of Several Parthenogenetic Species

The extant reptiles are one of the most diverse clades among terrestrial vertebrates and one of a few groups with instances of parthenogenesis. Due to the hybrid origin of parthenogenetic species, reference genomes of the parental species as well as of the parthenogenetic progeny are indispensable to explore the genetic foundations of parthenogenetic reproduction. Here, we report on the first genome assembly of rock lizard Darevskia valentini, a paternal species for several parthenogenetic lineages. The novel genome was used in the reconstruction of the comprehensive phylogeny of Squamata inferred independently from 7369 trees of single-copy orthologs and a supermatrix of 378 conserved proteins. We also investigated Hox clusters, the loci that are often regarded as playing an important role in the speciation of animal groups with drastically diverse morphology. We demonstrated that Hox clusters of D. valentini are invaded with transposons and contain the HoxC1 gene that has been considered to be lost in the amniote ancestor. This study provides confirmation for previous works and releases new genomic data that will contribute to future discoveries on the mechanisms of parthenogenesis as well as support comparative studies among reptiles.

The Degenerate Tale of Ascidian Tails

Ascidians are invertebrate chordates, with swimming chordate tadpole larvae that have distinct heads and tails. The head contains the small brain, sensory organs, including the ocellus (light) and otolith (gravity) and the presumptive endoderm, while the tail has a notochord surrounded by muscle cells and a dorsal nerve cord. One of the chordate features is a post-anal tail. Ascidian tadpoles are nonfeeding, and their tails are critical for larval locomotion. After hatching the larvae swim up toward light and are carried by the tide and ocean currents. When competent to settle, ascidian tadpole larvae swim down, away from light, to settle and metamorphose into a sessile adult. Tunicates are classified as chordates because of their chordate tadpole larvae; in contrast, the sessile adult has a U-shaped gut and very derived body plan, looking nothing like a chordate. There is one group of ascidians, the Molgulidae, where many species are known to have tailless larvae. The Swalla Lab has been studying the evolution of tailless ascidian larvae in this clade for over 30 years and has shown that tailless larvae have evolved independently several times in this clade. Comparison of the genomes of two closely related species, the tailed Molgula oculata and tailless Molgula occulta reveals much synteny, but there have been multiple insertions and deletions that have disrupted larval genes in the tailless species. Genomics and transcriptomics have previously shown that there are pseudogenes expressed in the tailless embryos, suggesting that the partial rescue of tailed features in their hybrid larvae is due to the expression of intact genes from the tailed parent. Yet surprisingly, we find that the notochord gene regulatory network is mostly intact in the tailless M. occulta, although the notochord does not converge and extend and remains as an aggregate of cells we call the "notoball." We expect that eventually many of the larval gene networks will become evolutionarily lost in tailless ascidians and the larval body plan abandoned, with eggs developing directly into an adult. Here we review the current evolutionary and developmental evidence on how the molgulids lost their tails.

Symbiont Transmission onto the Cell Surface of Early Oocytes in the Deep-Sea Clam Phreagena okutanii

Symbiotic associations with beneficial microorganisms endow a variety of host animals with adaptability to the environment. Stable transmission of symbionts across host generations is a key event in the maintenance of symbiotic associations through evolutionary time. However, our understanding of the mechanisms of symbiont transmission remains fragmentary. The deep-sea clam Phreagena okutanii harbors chemoautotrophic intracellular symbiotic bacteria in gill epithelial cells, and depends on these symbionts for nutrition. In this study, we focused on the association of these maternally transmitted symbionts with ovarian germ cells in juvenile female clams. First, we established a sex identification method for small P. okutanii individuals, and morphologically classified female germ cells observed in the ovary. Then, we investigated the association of the endosymbiotic bacteria with germ cells. We found that the symbionts were localized on the outer surface of the cell membrane of primary oocytes and not within the cluster of oogonia. Based on our findings, we discuss the processes and mechanisms of symbiont vertical transmission in P. okutanii.

Hox Genes Polymorphism Depicts Developmental Disruption of Common Sole Eggs

In sole aquaculture production, consistency in the quality of produced eggs throughout the year is unpredictable. Hox genes have a crucial role in controlling embryonic development and their genetic variation could alter the phenotype dramatically. In teleosts genome duplication led paralog hox genes to become diverged. Direct association of polymorphism in hoxa1a, hoxa2a & hoxa2b of Solea solea with egg viability indicates hoxa2b as a potential genetic marker. High Resolution Melt (HRM) analysis was carried out in 52 viable and 61 non-viable eggs collected at 54±6 hours post fertilization (hpf). Allelic and genotypic frequencies of polymorphism were analyzed and results illustrated a significantly increased risk for non-viability for minor alleles and their homozygous genotypes. Haplotype analysis demonstrated a significant recessive effect on the risk of non-viability, by increasing the odds of disrupting embryonic development up to three-fold. Phylogenetic analysis showed that the paralog genes hoxa2a and hoxa2b, are separated distinctly in two clades and presented a significant ω variation, revealing their diverged evolutionary rate.

The phylum Vertebrata: a case for zoological recognition

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

Hagfish and lamprey Hox genes reveal conservation of temporal colinearity in vertebrates

Hox genes exert fundamental roles for proper regional specification along the main rostro-caudal axis of animal embryos. They are generally expressed in restricted spatial domains according to their position in the cluster (spatial colinearity)-a feature that is conserved across bilaterians. In jawed vertebrates (gnathostomes), the position in the cluster also determines the onset of expression of Hox genes (a feature known as whole-cluster temporal colinearity (WTC)), while in invertebrates this phenomenon is displayed as a subcluster-level temporal colinearity. However, little is known about the expression profile of Hox genes in jawless vertebrates (cyclostomes); therefore, the evolutionary origin of WTC, as seen in gnathostomes, remains a mystery. Here, we show that Hox genes in cyclostomes are expressed according to WTC during development. We investigated the Hox repertoire and Hox gene expression profiles in three different species-a hagfish, a lamprey and a shark-encompassing the two major groups of vertebrates, and found that these are expressed following a whole-cluster, temporally staggered pattern, indicating that WTC has been conserved during the past 500 million years despite drastically different genome evolution and morphological outputs between jawless and jawed vertebrates.

Hox gene cluster of the ascidian, Halocynthia roretzi, reveals multiple ancient steps of cluster disintegration during ascidian evolution

BACKGROUND

Hox gene clusters with at least 13 paralog group (PG) members are common in vertebrate genomes and in that of amphioxus. Ascidians, which belong to the subphylum Tunicata (Urochordata), are phylogenetically positioned between vertebrates and amphioxus, and traditionally divided into two groups: the Pleurogona and the Enterogona. An enterogonan ascidian, Ciona intestinalis (Ci), possesses nine Hox genes localized on two chromosomes; thus, the Hox gene cluster is disintegrated. We investigated the Hox gene cluster of a pleurogonan ascidian, Halocynthia roretzi (Hr) to investigate whether Hox gene cluster disintegration is common among ascidians, and if so, how such disintegration occurred during ascidian or tunicate evolution.

RESULTS

Our phylogenetic analysis reveals that the Hr Hox gene complement comprises nine members, including one with a relatively divergent Hox homeodomain sequence. Eight of nine Hr Hox genes were orthologous to Ci-Hox1, 2, 3, 4, 5, 10, 12 and 13. Following the phylogenetic classification into 13 PGs, we designated Hr Hox genes as Hox1, 2, 3, 4, 5, 10, 11/12/13.a, 11/12/13.b and HoxX. To address the chromosomal arrangement of the nine Hox genes, we performed two-color chromosomal fluorescent in situ hybridization, which revealed that the nine Hox genes are localized on a single chromosome in Hr, distinct from their arrangement in Ci. We further examined the order of the nine Hox genes on the chromosome by chromosome/scaffold walking. This analysis suggested a gene order of Hox1, 11/12/13.b, 11/12/13.a, 10, 5, X, followed by either Hox4, 3, 2 or Hox2, 3, 4 on the chromosome. Based on the present results and those previously reported in Ci, we discuss the establishment of the Hox gene complement and disintegration of Hox gene clusters during the course of ascidian or tunicate evolution.

CONCLUSIONS

The Hox gene cluster and the genome must have experienced extensive reorganization during the course of evolution from the ancestral tunicate to Hr and Ci. Nevertheless, some features are shared in Hox gene components and gene arrangement on the chromosomes, suggesting that Hox gene cluster disintegration in ascidians involved early events common to tunicates as well as later ascidian lineage-specific events.

Segmental arithmetic: summing up the Hox gene regulatory network for hindbrain development in chordates

Organization and development of the early vertebrate hindbrain are controlled by a cascade of regulatory interactions that govern the process of segmentation and patterning along the anterior-posterior axis via Hox genes. These interactions can be assembled into a gene regulatory network that provides a framework to interpret experimental data, generate hypotheses, and identify gaps in our understanding of the progressive process of hindbrain segmentation. The network can be broadly separated into a series of interconnected programs that govern early signaling, segmental subdivision, secondary signaling, segmentation, and ultimately specification of segmental identity. Hox genes play crucial roles in multiple programs within this network. Furthermore, the network reveals properties and principles that are likely to be general to other complex developmental systems. Data from vertebrate and invertebrate chordate models are shedding light on the origin and diversification of the network. Comprehensive cis-regulatory analyses of vertebrate Hox gene regulation have enabled powerful cross-species gene regulatory comparisons. Such an approach in the sea lamprey has revealed that the network mediating segmental Hox expression was present in ancestral vertebrates and has been maintained across diverse vertebrate lineages. Invertebrate chordates lack hindbrain segmentation but exhibit conservation of some aspects of the network, such as a role for retinoic acid in establishing nested Hox expression domains. These comparisons lead to a model in which early vertebrates underwent an elaboration of the network between anterior-posterior patterning and Hox gene expression, leading to the gene-regulatory programs for segmental subdivision and rhombomeric segmentation. WIREs Dev Biol 2017, 6:e286. doi: 10.1002/wdev.286 For further resources related to this article, please visit the WIREs website.

Hox cluster characterization of Banna caecilian (Ichthyophis bannanicus) provides hints for slow evolution of its genome

BACKGROUND

Caecilians, with a discrete lifestyle, are the least explored group of amphibians. Though with distinct traits, many aspects of their biology are poorly investigated. Obtaining the caecilian genomic sequences will offer new perspectives and aid the fundamental studies in caecilian biology. The caecilian genomic sequences are also important and practical in the comparative genomics of amphibians. Currently, however, only sparse genomic sequences of caecilians are available. Hox genes, an old family of transcription factors playing central roles in the establishment of metazoan body plan. Understanding their structure and genomic organization may provide insights into the animal's genome, which is valuable for animals without a sequenced genome.

RESULTS

We sequenced and characterized the Hox clusters of Banna caecilian (Ichthyophis bannanicus) with a strategy combining long range PCR and genome walking. We obtained the majority of the four caecilian Hox clusters and identified 39 Hox genes, 5 microRNA genes and 1 pseudogene (ψHoxD12). There remained seven intergenic gaps we were unable to fill. From the obtained sequences, the caecilian Hox clusters contained less repetitive sequences and more conserved noncoding elements (CNEs) than the frog counterparts. We found that caecilian and coelacanth shared many more CNEs than frog and coelacanth did. Relative rate of sequence evolution showed that caecilian Hox genes evolved significantly more slowly than the other tetrapod species used in this study and were comparable to the slowly evolving coelacanth Hox genes. Phylogenetic tree of the four Hox clusters also revealed shorter branch length especially for the caecilian HoxA, HoxB and HoxD clusters. These features of the caecilian Hox clusters suggested a slowly evolving genome, which was supported by further analysis of a large orthologous protein dataset.

CONCLUSIONS

Our analyses greatly extended the knowledge about the caecilian Hox clusters from previous PCR surveys. From the obtained Hox sequences and the orthologous protein dataset, the caecilian Hox loci and its genome appear evolving comparatively slowly. As the basal lineage of amphibians and land vertebrate, this characteristic of the caecilian genome is valuable in the study concerning the genome biology and evolution of amphibians and early tetrapods.

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

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

A cDNA resource for gene expression studies of a hemichordate, Ptychodera flava

Recent investigations into the evolution of deuterostomes and the origin of chordates have paid considerable attention to hemichordates (acorn worms), as hemichordates and echinoderms are the closest chordate relatives. The present study prepared cDNA libraries from Ptychodera flava, to study expression and function of genes involved in development of the hemichordate body plan. Expressed sequence tag (EST) analyses of nine cDNA libraries yielded 18,832 cloned genes expressed in eggs, 18,739 in blastulae, 18,539 in gastrulae, 18,811 in larvae, 18,978 in juveniles, 11,802 in adult proboscis, 17,259 in stomochord, 11,886 in gills, and 11,580 in liver, respectively. A set of 34,159 uni-gene clones of P. flava was obtained. This cDNA resource will be valuable for studying temporal and spatial expression of acorn worm genes during development.

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

Background

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

Results

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

Conclusions

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

Evidence for at least six Hox clusters in the Japanese lamprey (Lethenteron japonicum)

Cyclostomes, comprising jawless vertebrates such as lampreys and hagfishes, are the sister group of living jawed vertebrates (gnathostomes) and hence an important group for understanding the origin and diversity of vertebrates. In vertebrates and other metazoans, Hox genes determine cell fate along the anteroposterior axis of embryos and are implicated in driving morphological diversity. Invertebrates contain a single Hox cluster (either intact or fragmented), whereas elephant shark, coelacanth, and tetrapods contain four Hox clusters owing to two rounds of whole-genome duplication ("1R" and "2R") during early vertebrate evolution. By contrast, most teleost fishes contain up to eight Hox clusters because of an additional "teleost-specific" genome duplication event. By sequencing bacterial artificial chromosome (BAC) clones and the whole genome, here we provide evidence for at least six Hox clusters in the Japanese lamprey (Lethenteron japonicum). This suggests that the lamprey lineage has experienced an additional genome duplication after 1R and 2R. The relative age of lamprey and human paralogs supports this hypothesis. Compared with gnathostome Hox clusters, lamprey Hox clusters are unusually large. Several conserved noncoding elements (CNEs) were predicted in the Hox clusters of lamprey, elephant shark, and human. Transgenic zebrafish assay indicated the potential of CNEs to function as enhancers. Interestingly, CNEs in individual lamprey Hox clusters are frequently conserved in multiple Hox clusters in elephant shark and human, implying a many-to-many orthology relationship between lamprey and gnathostome Hox clusters. Such a relationship suggests that the first two rounds of genome duplication may have occurred independently in the lamprey and gnathostome lineages.

Evolution of Hox gene clusters in deuterostomes

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

Broken colinearity of the amphioxus Hox cluster

Background

In most eumetazoans studied so far, Hox genes determine the identity of structures along the main body axis. They are usually linked in genomic clusters and, in the case of the vertebrate embryo, are expressed with spatial and temporal colinearity. Outside vertebrates, temporal colinearity has been reported in the cephalochordate amphioxus (the least derived living relative of the chordate ancestor) but only for anterior and central genes, namely Hox1 to Hox4 and Hox6. However, most of the Hox gene expression patterns in amphioxus have not been reported. To gain global insights into the evolution of Hox clusters in chordates, we investigated a more extended expression profile of amphioxus Hox genes.

Results

Here we report an extended expression profile of the European amphioxus Branchiostoma lanceolatum Hox genes and describe that all Hox genes, except Hox13, are expressed during development. Interestingly, we report the breaking of both spatial and temporal colinearity for at least Hox6 and Hox14, which thus have escaped from the classical Hox code concept. We show a previously unidentified Hox6 expression pattern and a faint expression for posterior Hox genes in structures such as the posterior mesoderm, notochord, and hindgut. Unexpectedly, we found that amphioxus Hox14 had the most divergent expression pattern. This gene is expressed in the anterior cerebral vesicle and pharyngeal endoderm. Amphioxus Hox14 expression represents the first report of Hox gene expression in the most anterior part of the central nervous system. Nevertheless, despite these divergent expression patterns, amphioxus Hox6 and Hox14 seem to be still regulated by retinoic acid.

Conclusions

Escape from colinearity by Hox genes is not unusual in either vertebrates or amphioxus and we suggest that those genes escaping from it are probably associated with the patterning of lineage-specific morphological traits, requiring the loss of those developmental constraints that kept them colinear.

Evolution of coding and non-coding genes in HOX clusters of a marsupial

Background

The HOX gene clusters are thought to be highly conserved amongst mammals and other vertebrates, but the long non-coding RNAs have only been studied in detail in human and mouse. The sequencing of the kangaroo genome provides an opportunity to use comparative analyses to compare the HOX clusters of a mammal with a distinct body plan to those of other mammals.

Results

Here we report a comparative analysis of HOX gene clusters between an Australian marsupial of the kangaroo family and the eutherians. There was a strikingly high level of conservation of HOX gene sequence and structure and non-protein coding genes including the microRNAs miR-196a, miR-196b, miR-10a and miR-10b and the long non-coding RNAs HOTAIR, HOTAIRM1 and HOXA11AS that play critical roles in regulating gene expression and controlling development. By microRNA deep sequencing and comparative genomic analyses, two conserved microRNAs (miR-10a and miR-10b) were identified and one new candidate microRNA with typical hairpin precursor structure that is expressed in both fibroblasts and testes was found. The prediction of microRNA target analysis showed that several known microRNA targets, such as miR-10, miR-414 and miR-464, were found in the tammar HOX clusters. In addition, several novel and putative miRNAs were identified that originated from elsewhere in the tammar genome and that target the tammar HOXB and HOXD clusters.

Conclusions

This study confirms that the emergence of known long non-coding RNAs in the HOX clusters clearly predate the marsupial-eutherian divergence 160 Ma ago. It also identified a new potentially functional microRNA as well as conserved miRNAs. These non-coding RNAs may participate in the regulation of HOX genes to influence the body plan of this marsupial.

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 anterior Hox regulatory elements among chordates

Background

The Hox family of transcription factors has a fundamental role in segmentation pathways and axial patterning of embryonic development and their clustered organization is linked with the regulatory mechanisms governing their coordinated expression along embryonic axes. Among chordates, of particular interest are the Hox paralogous genes in groups 1-4 since their expression is coupled to the control of regional identity in the anterior nervous system, where the highest structural diversity is observed.

Results

To investigate the degree of conservation in cis-regulatory components that form the basis of Hox expression in the anterior nervous system, we have used assays for transcriptional activity in ascidians and vertebrates to compare and contrast regulatory potential. We identified four regulatory sequences located near the CiHox1, CiHox2 and CiHox4 genes of the ascidian Ciona intestinalis which direct neural specific domains of expression. Using functional assays in Ciona and vertebrate embryos in combination with sequence analyses of enhancer fragments located in similar positions adjacent to Hox paralogy group genes, we compared the activity of these four Ciona cis-elements with a series of neural specific enhancers from the amphioxus Hox1-3 genes and from mouse Hox paralogous groups 1-4.

Conclusions

This analysis revealed that Kreisler and Krox20 dependent enhancers critical in segmental regulation of the hindbrain appear to be specific for the vertebrate lineage. In contrast, neural enhancers that function as Hox response elements through the action of Hox/Pbx binding motifs have been conserved during chordate evolution. The functional assays reveal that these Hox response cis-elements are recognized by the regulatory components of different and extant species. Together, our results indicate that during chordate evolution, cis-elements dependent upon Hox/Pbx regulatory complexes, are responsible for key aspects of segmental Hox expression in neural tissue and appeared with urochordates after cephalochordate divergence.

Evolutionary crossroads in developmental biology: amphioxus

The phylogenetic position of amphioxus, together with its relatively simple and evolutionarily conserved morphology and genome structure, has led to its use as a model for studies of vertebrate evolution. In particular, the recent development of technical approaches, as well as access to the complete amphioxus genome sequence, has provided the community with tools with which to study the invertebrate-chordate to vertebrate transition. Here, we present this animal model, discussing its life cycle, the model species studied and the experimental techniques that it is amenable to. We also summarize the major findings made using amphioxus that have informed us about the evolution of vertebrate traits.

Revisiting the origin of the vertebrate Hox14 by including its relict sarcopterygian members

Bilaterian Hox genes play pivotal roles in the specification of positional identities along the anteroposterior axis. Particularly in vertebrates, their regulation is tightly coordinated by tandem arrays of genes [paralogy groups (PGs)] in four gene clusters (HoxA-D). Traditionally, the uninterrupted Hox cluster (Hox1-14) of the invertebrate chordate amphioxus was regarded as an archetype of the vertebrate Hox clusters. In contrast to Hox1-13 that are globally regulated by the "Hox code" and are often phylogenetically conserved, vertebrate Hox14 members were only recently revealed to be present in an African lungfish, a coelacanth, chondrichthyans and a lamprey, and decoupled from the Hox code. In this study we performed a PCR-based search of Hox14 members from diverse vertebrates, and identified one in the Australian lungfish, Neoceratodus forsteri. Based on a molecular phylogenetic analysis, this gene was designated NfHoxA14. Our real-time RT-PCR suggested its hindgut-associated expression, previously observed also in cloudy catshark HoxD14 and lamprey Hox14α. It is likely that this altered expression scheme was established before the Hox cluster quadruplication, probably at the base of extant vertebrates. To investigate the origin of vertebrate Hox14, by including this sarcopterygian Hox14 member, we performed focused phylogenetic analyses on its relationship with other vertebrate posterior Hox PGs (Hox9-13) as well as amphioxus posterior Hox genes. Our results confirmed the hypotheses previously proposed by other studies that vertebrate Hox14 does not have any amphioxus ortholog, and that none of 1-to-1 pairs of vertebrate and amphioxus posterior Hox genes, based on their relative location in the clusters, is orthologous.

Transfection-based genomic readout for identifying rare transcriptional splice variants

Understanding the transcriptome, defined as the complete transcriptional component of the genome, is far more complex than originally considered. Even with the near fully resolved human and mouse genomes, for which extensive databases of transcribed sequence data (e.g., expressed sequence tags) are available, it is presently not possible to experimentally recover or computationally predict the full range of transcription products that derive from multiexon genes. Many genes are tightly regulated, which could include alternative processing of RNA, and lead to significant underrepresentation of many transcripts. A multitude of factors in addition to cell lineage- and developmental stage-specific expression as well as shortcomings in computational methods result in a less than complete understanding of transcriptional complexity. Here, we describe an approach to predict and evaluate a more complete repertoire of transcriptional products that derive from specific genetic loci with attention toward analysis of immune receptor genes. This approach is particularly useful in identifying gene products, including alternative splice forms, that originate from complex multigene families.

Evolution of conserved non-coding sequences within the vertebrate Hox clusters through the two-round whole genome duplications revealed by phylogenetic footprinting analysis

As a result of two-round whole genome duplications, four or more paralogous Hox clusters exist in vertebrate genomes. The paralogous genes in the Hox clusters show similar expression patterns, implying shared regulatory mechanisms for expression of these genes. Previous studies partly revealed the expression mechanisms of Hox genes. However, cis-regulatory elements that control these paralogous gene expression are still poorly understood. Toward solving this problem, the authors searched conserved non-coding sequences (CNSs), which are candidates of cis-regulatory elements. When comparing orthologous Hox clusters of 19 vertebrate species, 208 intergenic conserved regions were found. The authors then searched for CNSs that were conserved not only between orthologous clusters but also among the four paralogous Hox clusters. The authors found three regions that are conserved among all the four clusters and eight regions that are conserved between intergenic regions of two paralogous Hox clusters. In total, 28 CNSs were identified in the paralogous Hox clusters, and nine of them were newly found in this study. One of these novel regions bears a RARE motif. These CNSs are candidates for gene expression regulatory regions among paralogous Hox clusters. The authors also compared vertebrate CNSs with amphioxus CNSs within the Hox cluster, and found that two CNSs in the HoxA and HoxB clusters retain homology with amphioxus CNSs through the two-round whole genome duplications.

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.

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.

Conservation of gene linkage in dispersed vertebrate NK homeobox clusters

Nk homeobox genes are important regulators of many different developmental processes including muscle, heart, central nervous system and sensory organ development. They are thought to have arisen as part of the ANTP megacluster, which also gave rise to Hox and ParaHox genes, and at least some NK genes remain tightly linked in all animals examined so far. The protostome-deuterostome ancestor probably contained a cluster of nine Nk genes: (Msx)-(Nk4/tinman)-(Nk3/bagpipe)-(Lbx/ladybird)-(Tlx/c15)-(Nk7)-(Nk6/hgtx)-(Nk1/slouch)-(Nk5/Hmx). Of these genes, only NKX2.6-NKX3.1, LBX1-TLX1 and LBX2-TLX2 remain tightly linked in humans. However, it is currently unclear whether this is unique to the human genome as we do not know which of these Nk genes are clustered in other vertebrates. This makes it difficult to assess whether the remaining linkages are due to selective pressures or because chance rearrangements have "missed" certain genes. In this paper, we identify all of the paralogs of these ancestrally clustered NK genes in several distinct vertebrates. We demonstrate that tight linkages of Lbx1-Tlx1, Lbx2-Tlx2 and Nkx3.1-Nkx2.6 have been widely maintained in both the ray-finned and lobe-finned fish lineages. Moreover, the recently duplicated Hmx2-Hmx3 genes are also tightly linked. Finally, we show that Lbx1-Tlx1 and Hmx2-Hmx3 are flanked by highly conserved noncoding elements, suggesting that shared regulatory regions may have resulted in evolutionary pressure to maintain these linkages. Consistent with this, these pairs of genes have overlapping expression domains. In contrast, Lbx2-Tlx2 and Nkx3.1-Nkx2.6, which do not seem to be coexpressed, are also not associated with conserved noncoding sequences, suggesting that an alternative mechanism may be responsible for the continued clustering of these genes.

Evolution of developmental regulation in the vertebrate FgfD subfamily

Fibroblast growth factors (Fgfs) encode small signaling proteins that help regulate embryo patterning. Fgfs fall into seven families, including FgfD. Nonvertebrate chordates have a single FgfD gene; mammals have three (Fgf8, Fgf17, and Fgf18); and teleosts have six (fgf8a, fgf8b, fgf17, fgf18a, fgf18b, and fgf24). What are the evolutionary processes that led to the structural duplication and functional diversification of FgfD genes during vertebrate phylogeny? To study this question, we investigated conserved syntenies, patterns of gene expression, and the distribution of conserved noncoding elements (CNEs) in FgfD genes of stickleback and zebrafish, and compared them with data from cephalochordates, urochordates, and mammals. Genomic analysis suggests that Fgf8, Fgf17, Fgf18, and Fgf24 arose in two rounds of whole genome duplication at the base of the vertebrate radiation; that fgf8 and fgf18 duplications occurred at the base of the teleost radiation; and that Fgf24 is an ohnolog that was lost in the mammalian lineage. Expression analysis suggests that ancestral subfunctions partitioned between gene duplicates and points to the evolution of novel expression domains. Analysis of CNEs, at least some of which are candidate regulatory elements, suggests that ancestral CNEs partitioned between gene duplicates. These results help explain the evolutionary pathways by which the developmentally important family of FgfD molecules arose and the deduced principles that guided FgfD evolution are likely applicable to the evolution of developmental regulation in many vertebrate multigene families.

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.

Complete HOX cluster characterization of the coelacanth provides further evidence for slow evolution of its genome

The living coelacanth is a lobe-finned fish that represents an early evolutionary departure from the lineage that led to land vertebrates, and is of extreme interest scientifically. It has changed very little in appearance from fossilized coelacanths of the Cretaceous (150 to 65 million years ago), and is often referred to as a "living fossil." An important general question is whether long-term stasis in morphological evolution is associated with stasis in genome evolution. To this end we have used targeted genome sequencing for acquiring 1,612,752 bp of high quality finished sequence encompassing the four HOX clusters of the Indonesian coelacanth Latimeria menadoensis. Detailed analyses were carried out on genomic structure, gene and repeat contents, conserved noncoding regions, and relative rates of sequence evolution in both coding and noncoding tracts. Our results demonstrate conclusively that the coelacanth HOX clusters are evolving comparatively slowly and that this taxon should serve as a viable outgroup for interpretation of the genomes of tetrapod species.

Control of vertebrate Hox clusters by remote and global cis-acting regulatory sequences

Despite apparently shared structural organisation and functional roles, vertebrate Hox genes are controlled by regulatory mechanisms rather distinct from those of the prototypic Drosophila Antennapedia (ANT-C) and Bithorax (BX-C) Complexes. If individual regulatory modules have been shown to recapitulate specific Hox expression patterns, other experimental studies underscore that vertebrate Hox clusters are controlled in many of their functions in a global manner, through distinct mechanisms. We will discuss the different models that have been proposed to account for these global regulatory modes. In this context, the studies of the regulation of the HoxD complex during limb development highlighted the role of global regulatory elements and the different mechanisms associated to transform a structural organisation into distinct temporal and spatial expression domains. We will further discuss how these mechanisms may have benefited from the structure of the vertebrate homeotic clusters and reciprocally contribute to shape their evolution towards an increased level of organisation and compaction.

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.

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.

Atypical relaxation of structural constraints in Hox gene clusters of the green anole lizard

Hox genes control many aspects of embryonic development in metazoans. Previous analyses of this gene family revealed a surprising diversity in terms of gene number and organization between various animal species. In vertebrates, Hox genes are grouped into tightly organized clusters, claimed to be devoid of repetitive sequences. Here, we report the genomic organization of the four Hox loci present in the green anole lizard and show that they have massively accumulated retrotransposons, leading to gene clusters larger in size when compared to other vertebrates. In addition, similar repeats are present in many other development-related gene-containing regions, also thought to be refractory to such repetitive elements. Transposable elements are major sources of genetic variations, including alterations of gene expression, and hence this situation, so far unique among vertebrates, may have been associated with the evolution of the spectacular realm of morphological variations in the body plans of Squamata. Finally, sequence alignments highlight some divergent evolution in highly conserved DNA regions between vertebrate Hox clusters, which may coincide with the emergence of mammalian-specific features.

Computational prediction of amphioxus microRNA genes and their targets

Using a homology search based on genomic survey sequence analysis and microRNA (miRNA) secondary structure, a total of 51 miRNA genes encoding 30 amphioxus miRNAs were identified. These amphioxus miRNAs belong to 27 miRNA families. Many miRNA genes were found in multiple locations and in more than one genomic context. Analysis of amphioxus miRNA genes revealed nearly 50% miRNA genes identified positioning within introns of protein-coding genes. MicroRNA genes were also positioned intergenically between known protein-coding genes. MicroRNA genes were also clustered throughout the genome, indicating the potential for the cotranscription and coordinate expression of miRNAs located in a given cluster. Cross-species comparison indicates that of the 27 miRNA families, 17 families are shared by both protostomia and deuterostomia, and 2, 3 and 5 families are protostome-specific, deuterostome invertebrate-specific and chordate-specific, respectively. Computational predictions of amphioxus miRNA targets, taking into account the relationship between miRNA target genes and their host genes/neighboring genes, show that 49 target sites for 34 potential target genes were complementary to 19 miRNAs. For the intronic miRNA genes co-expressed with the host gene, they may not be involved in the regulation of the host gene expression. However, some intergenic miRNAs could participate in the modulation of their neighboring gene expression.

Genome duplication and the origin of the vertebrate skeleton

During vertebrate embryonic development, tissue patterning and differentiation are regulated by members of multigene families. Evolutionary expansion of these families is thought to have played a role in the evolution of anatomical complexity, including the origins of new cell and tissue types. A defining feature of vertebrates is an endoskeleton, the primary components of which are cartilage and bone. The molecular control of skeletal patterning has been the subject of intensive investigation for over two decades. More recently, comparative studies of organisms at key phylogenetic positions have highlighted the importance of gene duplication in the evolutionary diversification of connective tissues. Understanding the natural histories of gene families involved in skeletogenesis is therefore central to the issue of vertebrate skeletal evolution.

Unexpectedly large number of conserved noncoding regions within the ancestral chordate Hox cluster

The single amphioxus Hox cluster contains 15 genes and may well resemble the ancestral chordate Hox cluster. We have sequenced the Hox genomic complement of the European amphioxus Branchiostoma lanceolatum and compared it to the American species, Branchiostoma floridae, by phylogenetic footprinting to gain insights into the evolution of Hox gene regulation in chordates. We found that Hox intergenic regions are largely conserved between the two amphioxus species, especially in the case of genes located at the 3' of the cluster, a trend previously observed in vertebrates. We further compared the amphioxus Hox cluster with the human HoxA, HoxB, HoxC, and HoxD clusters, finding several conserved noncoding regions, both in intergenic and intronic regions. This suggests that the regulation of Hox genes is highly conserved across chordates, consistent with the similar Hox expression patterns in vertebrates and amphioxus.