Hox cluster genomics in the horn shark, Heterodontus francisci

HOX-Gene Cluster Organization and Genome Duplications in Fishes and Mammals: Transcript Variant Distribution along the Anterior–Posterior Axis

Hox genes play a crucial role in morphogenesis, especially in anterior–posterior body axis patterning. The organization of Hox clusters in vertebrates is a result of several genome duplications: two rounds of duplication in the ancestors of all vertebrates and a third round that was specific for teleost fishes. Teleostei cluster structure has been significantly modified in the evolutionary processes by Hox gene losses and co-options, while mammals show no such tendency. In mammals, the Hox gene number in a single cluster is stable and generally large, and the numbers are similar to those in the Chondrichthyes. Hox gene alternative splicing activity slightly differs between fishes and mammals. Fishes and mammals have differences in their known alternative splicing activity for Hox gene distribution along the anterior–posterior body axis. The analyzed fish groups—the Coelacanthiformes, Chondrichthyes, and Teleostei—all have higher known alternative mRNA numbers from the anterior and posterior regions, whereas mammals have a more uniform Hox transcript distribution along this axis. In fishes, most Hox transcripts produce functioning proteins, whereas mammals have significantly more known transcripts that do not produce functioning proteins.

Shark and ray genomics for disentangling their morphological diversity and vertebrate evolution

Developmental studies of sharks and rays (elasmobranchs) have provided much insight into the process of morphological evolution of vertebrates. Although those studies are supposedly fueled by large-scale molecular sequencing information, whole-genome sequences of sharks and rays were made available only recently. One compelling difficulty of elasmobranch developmental biology is the low accessibility to embryonic study materials and their slow development. Another limiting factor is the relatively large size of their genomes. Moreover, their large body sizes restrict sustainable captive breeding, while their high body fluid osmolarity prevents reproducible cell culturing for in vitro experimentation, which has also limited our knowledge of their chromosomal organization for validation of genome sequencing products. This article focuses on egg-laying elasmobranch species used in developmental biology and provides an overview of the characteristics of the shark and ray genomes revealed to date. Developmental studies performed on a gene-by-gene basis are also reviewed from a whole-genome perspective. Among the popular regulatory genes studied in developmental biology, I scrutinize shark homologs of Wnt genes that highlight vanishing repertoires in many other vertebrate lineages, as well as Hox genes that underwent an unexpected modification unique to the elasmobranch lineage. These topics are discussed together with insights into the reconstruction of developmental programs in the common ancestor of vertebrates and its subsequent evolutionary trajectories that mark the features that are unique to, and those characterizing the diversity among, cartilaginous fishes.

Deciphering The Potential Role of Hox Genes in Pancreatic Cancer

The Hox gene family plays an important role in organogenesis and animal development. Currently, 39 Hox genes that are clustered in four chromosome regions have been identified in humans. Emerging evidence suggests that Hox genes are involved in the development of the pancreas. However, the expression of Hox genes in pancreatic tumor tissues has been investigated in only a few studies. In addition, whether specific Hox genes can promote or suppress cancer metastasis is not clear. In this article, we first review the recent progress in studies on the role of Hox genes in pancreatic cancer. By comparing the expression profiles of pancreatic cancer cells isolated from genetically engineered mice established in our laboratory with three different proliferative and metastatic abilities, we identified novel Hox genes that exhibited tumor-promoting activity in pancreatic cancer. Finally, a potential oncogenic mechanism of the Hox genes was hypothesized.

An atlas of anterior hox gene expression in the embryonic sea lamprey head: Hox-code evolution in vertebrates

In the hindbrain and the adjacent cranial neural crest (NC) cells of jawed vertebrates (gnathostomes), nested and segmentally-restricted domains of Hox gene expression provide a combinatorial Hox-code for specifying regional properties during head development. Extant jawless vertebrates, such as the sea lamprey (Petromyzon marinus), can provide insights into the evolution and diversification of this Hox-code in vertebrates. There is evidence for gnathostome-like spatial patterns of Hox expression in lamprey; however, the expression domains of the majority of lamprey hox genes from paralogy groups (PG) 1-4 are yet to be characterized, so it is unknown whether they are coupled to hindbrain segments (rhombomeres) and NC. In this study, we systematically describe the spatiotemporal expression of all 14 sea lamprey hox genes from PG1-PG4 in the developing hindbrain and pharynx to investigate the extent to which their expression conforms to the archetypal gnathostome hindbrain and pharyngeal hox-codes. We find many similarities in Hox expression between lamprey and gnathostome species, particularly in rhombomeric domains during hindbrain segmentation and in the cranial neural crest, enabling inference of aspects of Hox expression in the ancestral vertebrate embryonic head. These data are consistent with the idea that a Hox regulatory network underlying hindbrain segmentation is a pan vertebrate trait. We also reveal differences in hindbrain domains at later stages, as well as expression in the endostyle and in pharyngeal arch (PA) 1 mesoderm. Our analysis suggests that many Hox expression domains that are observed in extant gnathostomes were present in ancestral vertebrates but have been partitioned differently across Hox clusters in gnathostome and cyclostome lineages after duplication.

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.

Functional and Comparative Genomics of Hoxa2 Gene cis-Regulatory Elements: Evidence for Evolutionary Modification of Ancestral Core Element Activity

Hoxa2 is an evolutionarily conserved developmental regulatory gene that functions to specify rhombomere (r) and pharyngeal arch (PA) identities throughout the Osteichthyes. Japanese medaka (Oryzias latipes) hoxa2a, like orthologous Hoxa2 genes from other osteichthyans, is expressed during embryogenesis in r2–7 and PA2-7, whereas the paralogous medaka pseudogene, ψhoxa2b, is expressed in noncanonical Hoxa2 domains, including the pectoral fin buds. To understand the evolution of cis-regulatory element (CRE) control of gene expression, we conducted eGFP reporter gene expression studies with extensive functional mapping of several conserved CREs upstream of medaka hoxa2a and ψhoxa2b in transient and stable-line transgenic medaka embryos. The CREs tested were previously shown to contribute to directing mouse Hoxa2 gene expression in r3, r5, and PA2-4. Our results reveal the presence of sequence elements embedded in the medaka hoxa2a and ψhoxa2b upstream enhancer regions (UERs) that mediate expression in r4 and the PAs (hoxa2a r4/CNCC element) or in r3–7 and the PAs ψhoxa2b r3–7/CNCC element), respectively. Further, these elements were shown to be highly conserved among osteichthyans, which suggests that the r4 specifying element embedded in the UER of Hoxa2 is a deeply rooted rhombomere specifying element and the activity of this element has been modified by the evolution of flanking sequences that redirect its activity to alternative developmental compartments.

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

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

SkateBase, an elasmobranch genome project and collection of molecular resources for chondrichthyan fishes

Chondrichthyan fishes are a diverse class of gnathostomes that provide a valuable perspective on fundamental characteristics shared by all jawed and limbed vertebrates. Studies of phylogeny, species diversity, population structure, conservation, and physiology are accelerated by genomic, transcriptomic and protein sequence data. These data are widely available for many sarcopterygii (coelacanth, lungfish and tetrapods) and actinoptergii (ray-finned fish including teleosts) taxa, but limited for chondrichthyan fishes.  In this study, we summarize available data for chondrichthyes and describe resources for one of the largest projects to characterize one of these fish, Leucoraja erinacea, the little skate.  SkateBase ( http://skatebase.org) serves as the skate genome project portal linking data, research tools, and teaching resources.

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.

Hox gene clusters of early vertebrates: do they serve as reliable markers for genome evolution?

Hox genes, responsible for regional specification along the anteroposterior axis in embryogenesis, are found as clusters in most eumetazoan genomes sequenced to date. Invertebrates possess a single Hox gene cluster with some exceptions of secondary cluster breakages, while osteichthyans (bony vertebrates) have multiple Hox clusters. In tetrapods, four Hox clusters, derived from the so-called two-round whole genome duplications (2R-WGDs), are observed. Overall, the number of Hox gene clusters has been regarded as a reliable marker of ploidy levels in animal genomes. In fact, this scheme also fits the situations in teleost fishes that experienced an additional WGD. In this review, I focus on cyclostomes and cartilaginous fishes as lineages that would fill the gap between invertebrates and osteichthyans. A recent study highlighted a possible loss of the HoxC cluster in the galeomorph shark lineage, while other aspects of cartilaginous fish Hox clusters usually mark their conserved nature. In contrast, existing resources suggest that the cyclostomes exhibit a different mode of Hox cluster organization. For this group of species, whose genomes could have differently responded to the 2R-WGDs from jawed vertebrates, therefore the number of Hox clusters may not serve as a good indicator of their ploidy level.

Body plan convergence in the evolution of skates and rays (Chondrichthyes: Batoidea)

Skates, rays and allies (Batoidea) comprise more than half of the species diversity and much of the morphological disparity among chondrichthyan fishes, the sister group to all other jawed vertebrates. While batoids are morphologically well characterized and have an excellent fossil record, there is currently no consensus on the interrelationships of family-level taxa. Here we construct a resolved, robust and time-calibrated batoid phylogeny using mitochondrial genomes, nuclear genes, and fossils, sampling densely across taxa. Data partitioning schemes, biases in the sequence data, and the relative informativeness of each fossil are explored. The molecular phylogeny is largely congruent with morphology crownward in the tree, but the branching orders of major batoid groups are mostly novel. Body plan convergence appears to be widespread in batoids. A depressed, rounded pectoral disk supported to the snout tip by fin radials, common to skates and stingrays, is indicated to have been derived independently by each group, while the long, spiny rostrum of sawfishes similarly appears to be convergent with that of sawsharks, which are not batoids. The major extant batoid lineages are inferred to have arisen relatively rapidly from the Late Triassic into the Jurassic, with long stems followed by subsequent radiations in each group around the Cretaceous/Tertiary boundary. The fossil record indicates that batoids were affected with disproportionate severity by the end-Cretaceous extinction event.

Hox clusters of the bichir (Actinopterygii, Polypterus senegalus) highlight unique patterns of sequence evolution in gnathostome phylogeny

Teleost fishes have extra Hox gene clusters owing to shared or lineage-specific genome duplication events in rayfinned fish (actinopterygian) phylogeny. Hence, extrapolating between genome function of teleosts and human or even between different fish species is difficult. We have sequenced and analyzed Hox gene clusters of the Senegal bichir (Polypterus senegalus), an extant representative of the most basal actinopterygian lineage. Bichir possesses four Hox gene clusters (A, B, C, D); phylogenetic analysis supports their orthology to the four Hox gene clusters of the gnathostome ancestor. We have generated a comprehensive database of conserved Hox noncoding sequences that include cartilaginous, lobe-finned, and ray-finned fishes (bichir and teleosts). Our analysis identified putative and known Hox cis-regulatory sequences with differing depths of conservation in Gnathostoma. We found that although bichir possesses four Hox gene clusters, its pattern of conservation of noncoding sequences is mosaic between outgroups, such as human, coelacanth, and shark, with four Hox gene clusters and teleosts, such as zebrafish and pufferfish, with seven or eight Hox gene clusters. Notably, bichir Hox gene clusters have been invaded by DNA transposons and this trend is further exemplified in teleosts, suggesting an as yet unrecognized mechanism of genome evolution that may explain Hox cluster plasticity in actinopterygians. Taken together, our results suggest that actinopterygian Hox gene clusters experienced a reduction in selective constraints that surprisingly predates the teleost-specific genome duplication.

A general scenario of Hox gene inventory variation among major sarcopterygian lineages

BACKGROUND

Hox genes are known to play a key role in shaping the body plan of metazoans. Evolutionary dynamics of these genes is therefore essential in explaining patterns of evolutionary diversity. Among extant sarcopterygians comprising both lobe-finned fishes and tetrapods, our knowledge of the Hox genes and clusters has largely been restricted in several model organisms such as frogs, birds and mammals. Some evolutionary gaps still exist, especially for those groups with derived body morphology or occupying key positions on the tree of life, hindering our understanding of how Hox gene inventory varied along the sarcopterygian lineage.

RESULTS

We determined the Hox gene inventory for six sarcopterygian groups: lungfishes, caecilians, salamanders, snakes, turtles and crocodiles by comprehensive PCR survey and genome walking. Variable Hox genes in each of the six sarcopterygian group representatives, compared to the human Hox gene inventory, were further validated for their presence/absence by PCR survey in a number of related species representing a broad evolutionary coverage of the group. Turtles, crocodiles, birds and placental mammals possess the same 39 Hox genes. HoxD12 is absent in snakes, amphibians and probably lungfishes. HoxB13 is lost in frogs and caecilians. Lobe-finned fishes, amphibians and squamate reptiles possess HoxC3. HoxC1 is only present in caecilians and lobe-finned fishes. Similar to coelacanths, lungfishes also possess HoxA14, which is only found in lobe-finned fishes to date. Our Hox gene variation data favor the lungfish-tetrapod, turtle-archosaur and frog-salamander relationships and imply that the loss of HoxD12 is not directly related to digit reduction.

CONCLUSIONS

Our newly determined Hox inventory data provide a more complete scenario for evolutionary dynamics of Hox genes along the sarcopterygian lineage. Limbless, worm-like caecilians and snakes possess similar Hox gene inventories to animals with less derived body morphology, suggesting changes to their body morphology are likely due to other modifications rather than changes to Hox gene numbers. Furthermore, our results provide basis for future sequencing of the entire Hox clusters of these animals.

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.

Genome enablement of the notothenioidei: genome size estimates from 11 species and BAC libraries from 2 representative taxa

The perciform suborder Notothenoidei provides a compelling opportunity to study the adaptive radiation of a marine species flock in the cold Southern Ocean surrounding Antarctica. To enable genome-level studies of these psychrophilic fishes, we estimated the sizes of the genomes of 11 Antarctic species and generated high-quality BAC libraries for 2, the notothen Notothenia coriiceps and the icefish Chaenocephalus aceratus. Our results indicate that evolution of phylogenetically derived notothenioid families, [e.g., the icefishes (Channichthyidae)], was accompanied by genome expansion. Species (n=6) of the basal family Nototheniidae had C values that ranged between 0.98 and 1.20 pg, whereas those of the icefishes, the notothenioid crown group, were 1.66-1.83 pg (n=4 species). The BAC libraries VMRC-19 (N. coriiceps) and VMRC-21 (C. aceratus) comprised 12X and 10X coverage of the respective genomes and had average insert sizes of 138 and 168 kb. Greater than 60% of paired BAC ends sampled from each library ( approximately 0.1% of each genome) contained repetitive sequences, and the repetitive element landscapes of the 2 genomes (13.4% of the N. coriiceps genome and 14.5% for C. aceratus) were similar. The representation and depth of coverage of the libraries were verified by identification of multiple Hox gene contigs: six discrete Hox clusters were found in N. coriiceps and at least five Hox clusters were found in C. aceratus. Given the unusual anatomical and physiological adaptations of the notothenioids, the availability of these BAC libraries sets the stage for expanded analysis of the psychrophilic mode of life.

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.

The evolution of gnathostome development: Insight from chondrichthyan embryology

Chondrichthyans (cartilaginous fishes) represent one of the two lineages of gnathostomes, the other being the osteicthyans (bony fishes). Classical studies on chondrichthyan embryology have strongly impacted our views of vertebrate body plan evolution, while recent studies highlight oviparous chondrichthyans as emerging vertebrate model systems that are amenable to experimental embryological manipulation. Here, we review three particular areas of interest in the field of chondrichthyan developmental biology-gastrulation, neural development, and appendage patterning-and we discuss recent findings within a broader chondrichthyan-osteichthyan comparative framework. In some cases, comparative studies of chondrichthyan and osteichthyan development reveal conserved patterns of gene expression in common developmental contexts. Studies of chondrichthyan gastrulation reveal conserved patterns of developmental gene expression, despite highly divergent modes of mesendoderm internalization, while molecular characterization of chondrichthyan neurogenic placodes indicates a conservation of placode transcription factor expression across gnathostome phylogeny. In other cases, comparative studies of chondrichthyan and osteichthyan development yield evidence of shared patterning mechanisms functioning in different developmental contexts. This is exemplified by studies on the development of chondrichthyan appendages-paired fins, median fins, and gill rays. These have demonstrated that a retinoic acid-responsive Shh-expressing signaling center functions to pattern the endoskeleton of gnathostome paired fins and chondrichthyan gill rays, while expression patterns of Tbx18 and HoxD family members are shared by gnathostome paired fins and chondrichthyan median fins. These findings fuel novel hypotheses of developmental genetic homology, and demonstrate how comparative studies of gnathostome development can provide insight into the evolutionary processes that underlie morphological diversity.

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.

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.

Comparative genomics of chondrichthyan Hoxa clusters

BACKGROUND

The chondrichthyan or cartilaginous fish (chimeras, sharks, skates and rays) occupy an important phylogenetic position as the sister group to all other jawed vertebrates and as an early lineage to diverge from the vertebrate lineage following two whole genome duplication events in vertebrate evolution. There have been few comparative genomic analyses incorporating data from chondrichthyan fish and none comparing genomic information from within the group. We have sequenced the complete Hoxa cluster of the Little Skate (Leucoraja erinacea) and compared to the published Hoxa cluster of the Horn Shark (Heterodontus francisci) and to available data from the Elephant Shark (Callorhinchus milii) genome project.

RESULTS

A BAC clone containing the full Little Skate Hoxa cluster was fully sequenced and assembled. Analyses of coding sequences and conserved non-coding elements reveal a strikingly high level of conservation across the cartilaginous fish, with twenty ultraconserved elements (100%,100 bp) found between Skate and Horn Shark, compared to three between human and marsupials. We have also identified novel potential non-coding RNAs in the Skate BAC clone, some of which are conserved to other species.

CONCLUSION

We find that the Little Skate Hoxa cluster is remarkably similar to the previously published Horn Shark Hoxa cluster with respect to sequence identity, gene size and intergenic distance despite over 180 million years of separation between the two lineages. We suggest that the genomes of cartilaginous fish are more highly conserved than those of tetrapods or teleost fish and so are more likely to have retained ancestral non-coding elements. While useful for isolating homologous DNA, this complicates bioinformatic approaches to identify chondrichthyan-specific non-coding DNA elements.

Molecular organization of 5S rDNA in sharks of the genus Rhizoprionodon: insights into the evolutionary dynamics of 5S rDNA in vertebrate genomes

In this study, we attempted a molecular characterization of the 5S rDNA in two closely related species of carcharhiniform sharks, Rhizoprionodon lalandii and Rhizoprionodon porosus, as well as a further comparative analysis of available data on lampreys, several fish groups and other vertebrates. Our data show that Rhizoprionodon sharks carry two 5S rDNA classes in their genomes: a short repeat class (termed class I) composed of approximately 185 bp repeats, and a large repeat class (termed class II) arrayed in approximately 465 bp units. These classes were differentiated by several base substitutions in the 5S coding region and by completely distinct non-transcribed spacers (NTS). In class II, both species showed a similar composition for both the gene coding region and the NTS region. In contrast, class I varied extensively both within and between the two shark species. A comparative analysis of 5S rRNA gene sequences of elasmobranchs and other vertebrates showed that class I is closely related to the bony fishes, whereas the class II gene formed a separate cartilaginous clade. The presence of two variant classes of 5S rDNA in sharks likely maintains the tendency for dual ribosomal classes observed in other fish species. The present data regarding the 5S rDNA organization provide insights into the dynamics and evolution of this multigene family in the fish genome, and they may also be useful in clarifying aspects of vertebrate genome evolution.

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.

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.

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.

Amniote phylogenomics: testing evolutionary hypotheses with BAC library scanning and targeted clone analysis of large-scale DNA sequences from reptiles

Phylogenomics research integrating established principles of systematic biology and taking advantage of the wealth of DNA sequences being generated by genome science holds promise for answering long-standing evolutionary questions with orders of magnitude more primary data than in the past. Although it is unrealistic to expect whole-genome initiatives to proceed rapidly for commercially unimportant species such as reptiles, practical approaches utilizing genomic libraries of large-insert clones pave the way for a phylogenomics of species that are nevertheless essential for testing evolutionary hypotheses within a phylogenetic framework. This chapter reviews the case for adopting genome-enabled approaches to evolutionary studies and outlines a program for using bacterial artificial chromosome (BAC) libraries or plasmid libraries as a basis for completing "genome scans" of reptiles. We have used BACs to close a critical gap in the genome database for Reptilia, the sister group of mammals, and present the methodological approaches taken to achieve this as a guideline for designing similar comparative studies. In addition, we provide a detailed step-by-step protocol for BAC-library screening and shotgun sequencing of specific clones containing target genes of evolutionary interest. Taken together, the genome scanning and shotgun sequencing techniques offer complementary diagnostic potential and can substantially increase the scale and power of analyses aimed at testing evolutionary hypotheses for nonmodel species.

Retracted: Gene duplication and functional evolution of Hox genes in fishes

With their power to shape animal morphology, few genes have captured the imagination of biologists as much as the evolutionarily conserved members of the Hox clusters. Hox genes encode transcription factors that play a key role in specifying the body plan in metazoans and are therefore essential in explaining patterns of evolutionary diversity. While each Hox cluster contains the same genes among the different mammalian species, this does not happen in ray-finned fish, in which both the number and organization of Hox genes and even Hox clusters are variable. Teleost fishes provide the first unambiguous support for ancient whole-genome duplication (third round) in an animal lineage. The number of genes differs in each cluster as a result of increased freedom to mutate after duplication. This has also allowed them to diverge and to adopt novel developmental roles. In this review, the authors have firstly focused on broadly outlining the duplication of Hoxgenes in fishes and discussing how comparative genomics is elucidating the molecular changes associated with the evolution of Hox genes expression and developmental function in the teleost fishes.Additional related research aspects, such as imaging of roles of microRNAs, chromatin regulation and evolutionary findings are also discussed.

An ultraconserved Hox-Pbx responsive element resides in the coding sequence of Hoxa2 and is active in rhombomere 4

The Hoxa2 gene has a fundamental role in vertebrate craniofacial and hindbrain patterning. Segmental control of Hoxa2 expression is crucial to its function and several studies have highlighted transcriptional regulatory elements governing its activity in distinct rhombomeres. Here, we identify a putative Hox–Pbx responsive cis-regulatory sequence, which resides in the coding sequence of Hoxa2 and is an important component of Hoxa2 regulation in rhombomere (r) 4. By using cell transfection and chromatin immunoprecipitation (ChIP) assays, we show that this regulatory sequence is responsive to paralogue group 1 and 2 Hox proteins and to their Pbx co-factors. Importantly, we also show that the Hox–Pbx element cooperates with a previously reported Hoxa2 r4 intronic enhancer and that its integrity is required to drive specific reporter gene expression in r4 upon electroporation in the chick embryo hindbrain. Thus, both intronic as well as exonic regulatory sequences are involved in Hoxa2 segmental regulation in the developing r4. Finally, we found that the Hox–Pbx exonic element is embedded in a larger 205-bp long ultraconserved genomic element (UCE) shared by all vertebrate genomes. In this respect, our data further support the idea that extreme conservation of UCE sequences may be the result of multiple superposed functional and evolutionary constraints.

Biphasic Hoxd gene expression in shark paired fins reveals an ancient origin of the distal limb domain

The evolutionary transition of fins to limbs involved development of a new suite of distal skeletal structures, the digits. During tetrapod limb development, genes at the 5' end of the HoxD cluster are expressed in two spatiotemporally distinct phases. In the first phase, Hoxd9-13 are activated sequentially and form nested domains along the anteroposterior axis of the limb. This initial phase patterns the limb from its proximal limit to the middle of the forearm. Later in development, a second wave of transcription results in 5' HoxD gene expression along the distal end of the limb bud, which regulates formation of digits. Studies of zebrafish fins showed that the second phase of Hox expression does not occur, leading to the idea that the origin of digits was driven by addition of the distal Hox expression domain in the earliest tetrapods. Here we test this hypothesis by investigating Hoxd gene expression during paired fin development in the shark Scyliorhinus canicula, a member of the most basal lineage of jawed vertebrates. We report that at early stages, 5'Hoxd genes are expressed in anteroposteriorly nested patterns, consistent with the initial wave of Hoxd transcription in teleost and tetrapod paired appendages. Unexpectedly, a second phase of expression occurs at later stages of shark fin development, in which Hoxd12 and Hoxd13 are re-expressed along the distal margin of the fin buds. This second phase is similar to that observed in tetrapod limbs. The results indicate that a second, distal phase of Hoxd gene expression is not uniquely associated with tetrapod digit development, but is more likely a plesiomorphic condition present the common ancestor of chondrichthyans and osteichthyans. We propose that a temporal extension, rather than de novo activation, of Hoxd expression in the distal part of the fin may have led to the evolution of digits.

Survey sequencing and comparative analysis of the elephant shark (Callorhinchus milii) genome

Owing to their phylogenetic position, cartilaginous fishes (sharks, rays, skates, and chimaeras) provide a critical reference for our understanding of vertebrate genome evolution. The relatively small genome of the elephant shark, Callorhinchus milii, a chimaera, makes it an attractive model cartilaginous fish genome for whole-genome sequencing and comparative analysis. Here, the authors describe survey sequencing (1.4x coverage) and comparative analysis of the elephant shark genome, one of the first cartilaginous fish genomes to be sequenced to this depth. Repetitive sequences, represented mainly by a novel family of short interspersed element-like and long interspersed element-like sequences, account for about 28% of the elephant shark genome. Fragments of approximately 15,000 elephant shark genes reveal specific examples of genes that have been lost differentially during the evolution of tetrapod and teleost fish lineages. Interestingly, the degree of conserved synteny and conserved sequences between the human and elephant shark genomes are higher than that between human and teleost fish genomes. Elephant shark contains putative four Hox clusters indicating that, unlike teleost fish genomes, the elephant shark genome has not experienced an additional whole-genome duplication. These findings underscore the importance of the elephant shark as a critical reference vertebrate genome for comparative analysis of the human and other vertebrate genomes. This study also demonstrates that a survey-sequencing approach can be applied productively for comparative analysis of distantly related vertebrate genomes.

Evolutionary constraint on Otx2 neuroectoderm enhancers-deep conservation from skate to mouse and unique divergence in teleost

Otx2 is a paired type homeobox gene that plays essential roles in each step and site of head development in vertebrates. In the mouse, Otx2 expression in the anterior neuroectoderm is regulated primarily by two distinct enhancers: anterior neuroectoderm (AN) and forebrain/midbrain (FM) enhancers at 92 kb and 75 kb 5'of the Otx2 locus, respectively. The AN enhancer has activity in the entire anterior neuroectoderm at headfold and early somite stages, whereas the FM enhancer is subsequently active in the future caudal forebrain and midbrain ectoderm. In tetrapods, both AN and FM enhancers are conserved, whereas the AN region is missing in teleosts, despite overt Otx2 expression in the anterior neuroectoderm. Here, we show that zebrafish and fugu FM regions drive expression not only in the forebrain and midbrain but also in the anterior neuroectoderm at headfold stage. The analysis of coelacanth and skate genomic Otx2 orthologues suggests that the utilization of the two enhancers, AN and FM, is an ancestral condition. In contrast, the AN enhancer has been specifically lost in the teleost lineage with a compensatory establishment of AN activity within the FM enhancer. Furthermore, the AN activity in the fish FM enhancer was established by recruiting upstream factors different from those that direct the tetrapod AN enhancer, yet zebrafish FM enhancer is active in both mouse and zebrafish anterior neuroectoderm at the headfold stage.

Comparative genomics in vertebrate evolution and development

The vast quantities of publicly available DNA sequencing data and genome resources are enabling biologists to investigate age-old problems in biology that were not addressable previously. In this review, we discuss how comparative genomics is practiced and how the data can be used to make biological inferences with respect to vertebrate evolution and development. Examples are taken from the well-known HOX clusters, which are always a high-priority target for genomic analyses due to their inferred role in the evolution of metazoans. In addition, we briefly discuss the application of genomic approaches to problems in comparative endocrinology.

Construction of a nurse shark (Ginglymostoma cirratum) bacterial artificial chromosome (BAC) library and a preliminary genome survey

BACKGROUND

Sharks are members of the taxonomic class Chondrichthyes, the oldest living jawed vertebrates. Genomic studies of this group, in comparison to representative species in other vertebrate taxa, will allow us to theorize about the fundamental genetic, developmental, and functional characteristics in the common ancestor of all jawed vertebrates.

AIMS

In order to obtain mapping and sequencing data for comparative genomics, we constructed a bacterial artificial chromosome (BAC) library for the nurse shark, Ginglymostoma cirratum.

RESULTS

The BAC library consists of 313,344 clones with an average insert size of 144 kb, covering ~4.5 x 1010 bp and thus providing an 11-fold coverage of the haploid genome. BAC end sequence analyses revealed, in addition to LINEs and SINEs commonly found in other animal and plant genomes, two new groups of nurse shark-specific repetitive elements, NSRE1 and NSRE2 that seem to be major components of the nurse shark genome. Screening the library with single-copy or multi-copy gene probes showed 6-28 primary positive clones per probe of which 50-90% were true positives, demonstrating that the BAC library is representative of the different regions of the nurse shark genome. Furthermore, some BAC clones contained multiple genes, making physical mapping feasible.

CONCLUSION

We have constructed a deep-coverage, high-quality, large insert, and publicly available BAC library for a cartilaginous fish. It will be very useful to the scientific community interested in shark genomic structure, comparative genomics, and functional studies. We found two new groups of repetitive elements specific to the nurse shark genome, which may contribute to the architecture and evolution of the nurse shark genome.

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.

Molecular evolution of duplicated ray finned fish HoxA clusters: increased synonymous substitution rate and asymmetrical co-divergence of coding and non-coding sequences

In this study the molecular evolution of duplicated HoxA genes in zebrafish and fugu has been investigated. All 18 duplicated HoxA genes studied have a higher non-synonymous substitution rate than the corresponding genes in either bichir or paddlefish, where these genes are not duplicated. The higher rate of evolution is not due solely to a higher non-synonymous-to-synonymous rate ratio but to an increase in both the non-synonymous as well as the synonymous substitution rate. The synonymous rate increase can be explained by a change in base composition, codon usage, or mutation rate. We found no changes in nucleotide composition or codon bias. Thus, we suggest that the HoxA genes may experience an increased mutation rate following cluster duplication. In the non-Hox nuclear gene RAG1 only an increase in non-synonymous substitutions could be detected, suggesting that the increased mutation rate is specific to duplicated Hox clusters and might be related to the structural instability of Hox clusters following duplication. The divergence among paralog genes tends to be asymmetric, with one paralog diverging faster than the other. In fugu, all b-paralogs diverge faster than the a-paralogs, while in zebrafish Hoxa-13a diverges faster. This asymmetry corresponds to the asymmetry in the divergence rate of conserved non-coding sequences, i.e., putative cis-regulatory elements. These results suggest that the 5' HoxA genes in the same cluster belong to a co-evolutionary unit in which genes have a tendency to diverge together.

An evolutionary constraint: strongly disfavored class of change in DNA sequence during divergence of cis-regulatory modules

The DNA of functional cis-regulatory modules displays extensive sequence conservation in comparisons of genomes from modestly distant species. Patches of sequence that are several hundred base pairs in length within these modules are often seen to be 80-95% identical, although the flanking sequence cannot even be aligned. However, it is unlikely that base pairs located between the transcription factor target sites of cis-regulatory modules have sequence-dependent function, and the mechanism that constrains evolutionary change within cis-regulatory modules is incompletely understood. We chose five functionally characterized cis-regulatory modules from the Strongylocentrotus purpuratus (sea urchin) genome and obtained orthologous regulatory and flanking sequences from a bacterial artificial chromosome genome library of a congener, Strongylocentrotus franciscanus. As expected, single-nucleotide substitutions and small indels occur freely at many positions within the regulatory modules of these two species, as they do outside the regulatory modules. However, large indels (>20 bp) are statistically almost absent within the regulatory modules, although they are common in flanking intergenic or intronic sequence. The result helps to explain the patterns of evolutionary sequence divergence characteristic of cis-regulatory DNA.

Hoxd13 expression in the developing limbs of the short-tailed fruit bat, Carollia perspicillata

Bat forelimbs are highly specialized for sustained flight, providing a unique model to explore the genetic programs that regulate vertebrate limb diversity. Hoxd9-13 genes are important regulators of stylopodium, zeugopodium, and autopodium development and thus evolutionary changes in their expression profiles and biochemical activities may contribute to divergent limb morphologies in vertebrates. We have isolated the genomic region that includes Hoxd12 and Hoxd13 from Carollia perspicillata, the short-tailed fruit bat. The bat Hoxd13 gene encodes a protein that shares 95% identity with human and mouse HOXD13. The expression pattern of bat Hoxd13 mRNA during limb development was compared with that of mouse. In bat and mouse hindlimbs, the expression patterns of Hoxd13 are relatively similar. However, although the forelimb Hoxd13 expression patterns in both organisms during early limb bud stages are similar, at later stages they diverge; the anterior expression boundary of bat Hoxd13 is posterior-shifted relative to the mouse. These findings, compared with the Hoxd13 expression profiles of other vertebrates, suggest that divergent Hoxd13 expression patterns may contribute to limb morphological variation.

Comparative cis-regulatory analyses identify new elements of the mouse Hoxc8 early enhancer

The Hoxc8 early enhancer is a 200 bp region that controls the early phase of Hoxc8 expression during mouse embryonic development. This enhancer defines the domain of Hoxc8 expression in the neural tube and mesoderm of the posterior regions of the developing embryo. Five distinct cis-acting elements, A-E, were previously shown to govern early phase Hoxc8 expression. Significant divergence between mammalian and fish Hoxc8 early enhancer sequences and activities suggested additional cis-acting elements. Here we describe four additional cis-acting elements (F-I) within the 200 bp Hoxc8 early enhancer region identified by comparative regulatory analysis and transgene-mutation studies. These elements affect posterior neural tube and mesoderm expression of the reporter gene, either singly or in combination. Surprisingly, these new elements are missing from the zebrafish and Fugu Hoxc8 early enhancer sequences. Considering that fish enhancers direct robust reporter expression in transgenic mouse embryos, it is tempting to postulate that fish and mammalian Hoxc8 early enhancers utilize different sets of elements to direct Hoxc8 early expression. These observations reveal a remarkable plasticity in the Hoxc8 early enhancer, suggesting different modes of initiation and establishment of Hoxc8 expression in different species. We postulate that extensive restructuring and remodeling of Hox cis-regulatory regions occurring in different taxa lead to relatively different Hox expression patterns, which in turn may act as a driving force in generating diverse axial morphologies.

The Dlx gene complement of the leopard shark, Triakis semifasciata, resembles that of mammals: implications for genomic and morphological evolution of jawed vertebrates

Extensive gene duplication is thought to have occurred in the vertebrate lineage after it diverged from cephalochordates and before the divergence of lobe- and ray-finned fishes, but the exact timing remains obscure. This timing was investigated by analysis of the Dlx gene family of a representative cartilaginous fish, the leopard shark, Triakis semifasciata. Dlx genes encode homeodomain transcription factors and are arranged in mammals as three convergently transcribed bigene clusters. Six Dlx genes were cloned from Triakis and shown to be orthologous to single mammalian Dlx genes. At least four of these are arranged in bigene clusters. Phylogenetic analyses of Dlx genes were used to propose an evolutionary scenario in which two genome duplications led to four Dlx bigene clusters in a common ancestor of jawed vertebrates, one of which was lost prior to the diversification of the group. Dlx genes are known to be involved in jaw development, and changes in Dlx gene number are mapped to the same branch of the vertebrate tree as the origin of jaws.

Evolution of microRNAs located withinHox gene clusters

MicroRNAs (miRNAs) form an abundant class of non-coding RNA genes that have an important function in post-transcriptional gene regulation and in particular modulate the expression of developmentally important transcription factors including Hox genes. Two families of microRNAs are genomically located in intergenic regions in the Hox clusters of vertebrates. Here we describe their evolution in detail. We show that the micro RNAs closely follow the patterns of protein evolution in the Hox clusters, which is characterized by cluster duplications followed by differential gene loss.

Organization of Hox genes in ascidians: Present, past, and future

Hox genes have been regarded to play a central role in anterior-posterior patterning of the animal body. Variations of Hox genes among animal species in the number, order on a chromosome, and the developmental expression pattern may reflect an evolutionary history. Therefore, it is definitely necessary to characterize Hox genes of wide variety of animal species, especially the species occupying key positions in the animal phylogeny. Ascidians, belonging to the subphylum Urochordata, are one of the sister groups of vertebrates in the phylum Chordata. Recent studies have shown that nine Hox genes of Ciona intestinalis, an ascidian species, are present on two chromosomes in the genome. In this review, we discuss the present state of Hox genes in ascidians, focusing on their novel chromosomal organization and expression pattern with unique features and how the novel organization has evolved in relation to the unique body plan of ascidians.

The shark HoxN cluster is homologous to the human HoxD cluster

The statistical analysis of phylogenetic footprints in the two known horn shark Hox clusters and the four mammalian clusters shows that the shark HoxN cluster is HoxD-like. This finding implies that the most recent common ancestor of jawed vertebrates had at least four Hox clusters, including those which are orthologous to the four mammalian Hox clusters.

BAC libraries and comparative genomics of aquatic chordate species

The bacterial artificial chromosome (BAC) system is useful for creating a representation of the genomes of target species. The system is advantageous in that it can accommodate exogenous inserts that are very large (>100 kilobases, kb), thereby allowing entire eukaryotic genes (including flanking regulatory regions) to be encompassed in a single clone. The interest in BACs has recently been spawned by vast improvements in high throughput genomic sequencing such that comparisons of orthologous regions from different genomes (comparative genomics) are being routinely investigated, and comprise a significant component, of all major sequencing centers. In this review, we discuss the general principles of BAC cloning, the resources that are currently available, and some of the applications of the technology. It is not intended to be an exhaustive treatise; rather our goal is to provide a primer of the BAC technology in order to make readers aware of these resources and how they may utilize them in their own research programs.

The duplication of the Hox gene clusters in teleost fishes

Higher teleost fishes, including zebrafish and fugu, have duplicated their Hox genes relative to the gene inventory of other gnathostome lineages. The most widely accepted theory contends that the duplicate Hox clusters orginated synchronously during a single genome duplication event in the early history of ray-finned fishes. In this contribution we collect and re-evaluate all publicly available sequence information. In particular, we show that the short Hox gene fragments from published PCR surveys of the killifish Fundulus heteroclitus, the medaka Oryzias latipes and the goldfish Carassius auratus can be used to determine with little ambiguity not only their paralog group but also their membership in a particular cluster.Together with a survey of the genomic sequence data from the pufferfish Tetraodon nigroviridis we show that at least percomorpha, and possibly all eutelosts, share a system of 7 or 8 orthologous Hox gene clusters. There is little doubt about the orthology of the two teleost duplicates of the HoxA and HoxB clusters. A careful analysis of both the coding sequence of Hox genes and of conserved non-coding sequences provides additional support for the "duplication early" hypothesis that the Hox clusters in teleosts are derived from eight ancestral clusters by means of subsequent gene loss; the data remain ambiguous, however, in particular for the HoxC clusters.Assuming the "duplication early" hypothesis we use the new evidence on the Hox gene complements to determine the phylogenetic positions of gene-loss events in the wake of the cluster duplication. Surprisingly, we find that the resolution of redundancy seems to be a slow process that is still ongoing. A few suggestions on which additional sequence data would be most informative for resolving the history of the teleostean Hox genes are discussed.

Comparative genomic analysis of vertebrate Hox3 and Hox4 genes

We used a comparative genomic approach to identify putative cis-acting regulatory sequences of the zebrafish hoxb3a and hoxb4a genes. We aligned genomic sequences spanning the clustered Hoxb1 to Hoxb5 genes from pufferfish, mice, and humans with the zebrafish hoxba and hoxbb cluster sequences. We identified multiple blocks of conserved sequences in non-coding regions within and surrounding the Hoxb3/b4 gene locus; a subset of these blocks are conserved in the zebrafish hoxbb cluster, despite loss of hoxb3/b4 genes. Overall, we find that the architecture of the Hoxb3/b4 loci and of the conserved sequence elements is very similar in teleosts and mammals. Our analyses also revealed two alternative transcripts of the zebrafish hoxb3a gene and an exon sequence unusually located 10 kb upstream of adjacent hoxb4a; an equivalent murine Hoxb3 exon has not yet been confirmed. We show that many of the Hoxb3/b4 conserved non-coding sequences correlate with functional neural enhancers previously described in the mouse. Further, within the conserved non-coding sequences we have identified binding sites for transcription factors, including Kreisler/Valentino, Krox20, Hox, and Pbx, some of which had not been previously described for the mouse. Finally, we demonstrate that the regulatory sequences of zebrafish hoxa3a are divergent with respect to the mouse ortholog Hoxa3, or the paralog hoxb3a. Despite limited conservation of regulatory sequences, zebrafish hoxa3a and hoxb3a genes share very similar expression profiles.