Unusual number and genomic organization of Hox genes in the tunicate Ciona intestinalis

Dissecting the evolutionary role of the Hox gene proboscipedia in Drosophila mouthpart diversification by full locus replacement

Hox genes determine positional codes along the head-to-tail axis. Here, we replaced the entire Drosophila melanogaster proboscipedia (pb) Hox locus, which controls the development of the proboscis and maxillary palps, with that from Drosophila mimica, a related species with highly modified mouthparts. The D. mimica replacement rescues most aspects of adult proboscis morphology; however, the shape and orientation of maxillary palps were modified, resembling D. mimica and closely related species. Expressing the D. mimica Pb protein in the D. melanogaster pattern fully rescued D. melanogaster morphology. However, the expression pattern directed by D. mimica pb cis-regulatory sequences differed from that of D. melanogaster pb in cells that produce altered maxillary structures, indicating that pb regulatory sequences can evolve in related species to alter morphology.

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.

Homeotic Genes: Clustering, Modularity, and Diversity

Hox genes code for transcription factors and are evolutionarily conserved. They regulate a plethora of downstream targets to define the anterior-posterior (AP) body axis of a developing bilaterian embryo. Early work suggested a possible role of clustering and ordering of Hox to regulate their expression in a spatially restricted manner along the AP axis. However, the recent availability of many genome assemblies for different organisms uncovered several examples that defy this constraint. With recent advancements in genomics, the current review discusses the arrangement of Hox in various organisms. Further, we revisit their discovery and regulation in Drosophila melanogaster. We also review their regulation in different arthropods and vertebrates, with a significant focus on Hox expression in the crustacean Parahyale hawaiensis. It is noteworthy that subtle changes in the levels of Hox gene expression can contribute to the development of novel features in an organism. We, therefore, delve into the distinct regulation of these genes during primary axis formation, segment identity, and extra-embryonic roles such as in the formation of hair follicles or misregulation leading to cancer. Toward the end of each section, we emphasize the possibilities of several experiments involving various organisms, owing to the advancements in the field of genomics and CRISPR-based genome engineering. Overall, we present a holistic view of the functioning of Hox in the animal world.

Towards reconstructing the dipteran demise of an ancient essential gene: E3 ubiquitin ligase Murine double minute

Genome studies have uncovered many examples of essential gene loss, raising the question of how ancient genes transition from essentiality to dispensability. We explored this process for the deeply conserved E3 ubiquitin ligase Murine double minute (Mdm), which is lacking in Drosophila despite the conservation of its main regulatory target, the cellular stress response gene p53. Conducting gene expression and knockdown experiments in the red flour beetle Tribolium castaneum, we found evidence that Mdm has remained essential in insects where it is present. Using bioinformatics approaches, we confirm the absence of the Mdm gene family in Drosophila, mapping its loss to the stem lineage of schizophoran Diptera and Pipunculidae (big-headed flies), about 95-85 million years ago. Intriguingly, this gene loss event was preceded by the de novo origin of the gene Companion of reaper (Corp), a novel p53 regulatory factor that is characterized by functional similarities to vertebrate Mdm2 despite lacking E3 ubiquitin ligase protein domains. Speaking against a 1:1 compensatory gene gain/loss scenario, however, we found that hoverflies (Syrphidae) and pointed-wing flies (Lonchopteridae) possess both Mdm and Corp. This implies that the two p53 regulators have been coexisting for ~ 150 million years in select dipteran clades and for at least 50 million years in the lineage to Schizophora and Pipunculidae. Given these extensive time spans of Mdm/Corp coexistence, we speculate that the loss of Mdm in the lineage to Drosophila involved further acquisitions of compensatory gene activities besides the emergence of Corp. Combined with the previously noted reduction of an ancestral P53 contact domain in the Mdm homologs of crustaceans and insects, we conclude that the loss of the ancient Mdm gene family in flies was the outcome of incremental functional regression over long macroevolutionary time scales.

The notochord gene regulatory network in chordate evolution: Conservation and divergence from Ciona to vertebrates

The notochord is a structure required for support and patterning of all chordate embryos, from sea squirts to humans. An increasing amount of information on notochord development and on the molecular strategies that ensure its proper morphogenesis has been gleaned through studies in the sea squirt Ciona. This invertebrate chordate offers a fortunate combination of experimental advantages, ranging from translucent, fast-developing embryos to a compact genome and impressive biomolecular resources. These assets have enabled the rapid identification of numerous notochord genes and cis-regulatory regions, and provide a rather unique opportunity to reconstruct the gene regulatory network that controls the formation of this developmental and evolutionary chordate landmark. This chapter summarizes the morphogenetic milestones that punctuate notochord formation in Ciona, their molecular effectors, and the current knowledge of the gene regulatory network that ensures the accurate spatial and temporal orchestration of these processes.

De novo draft assembly of the Botrylloides leachii genome provides further insight into tunicate evolution

Tunicates are marine invertebrates that compose the closest phylogenetic group to the vertebrates. These chordates present a particularly diverse range of regenerative abilities and life-history strategies. Consequently, tunicates provide an extraordinary perspective into the emergence and diversity of these traits. Here we describe the genome sequencing, annotation and analysis of the Stolidobranchian Botrylloides leachii. We have produced a high-quality 159 Mb assembly, 82% of the predicted 194 Mb genome. Analysing genome size, gene number, repetitive elements, orthologs clustering and gene ontology terms show that B. leachii has a genomic architecture similar to that of most solitary tunicates, while other recently sequenced colonial ascidians have undergone genome expansion. In addition, ortholog clustering has identified groups of candidate genes for the study of colonialism and whole-body regeneration. By analysing the structure and composition of conserved gene linkages, we observed examples of cluster breaks and gene dispersions, suggesting that several lineage-specific genome rearrangements occurred during tunicate evolution. We also found lineage-specific gene gain and loss within conserved cell-signalling pathways. Such examples of genetic changes within conserved cell-signalling pathways commonly associated with regeneration and development that may underlie some of the diverse regenerative abilities observed in tunicates. Overall, these results provide a novel resource for the study of tunicates and of colonial ascidians.

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.

Comparative transcriptomics enlarges the toolkit of known developmental genes in mollusks

BACKGROUND

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

RESULTS

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

CONCLUSIONS

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

T-Box Genes and Developmental Gene Regulatory Networks in Ascidians

Ascidians are invertebrate chordates with a biphasic life cycle characterized by a dual body plan that displays simplified versions of chordate structures, such as a premetamorphic 40-cell notochord topped by a dorsal nerve cord and postmetamorphic pharyngeal slits. These relatively simple chordates are characterized by rapid development, compact genomes and ease of transgenesis, and thus provide the opportunity to rapidly characterize the genomic organization, developmental function, and transcriptional regulation of evolutionarily conserved gene families. This review summarizes the current knowledge on members of the T-box family of transcription factors in Ciona and other ascidians. In both chordate and nonchordate animals, these genes control a variety of morphogenetic processes, and their mutations are responsible for malformations and developmental defects in organisms ranging from flies to humans. In ascidians, T-box transcription factors are required for the formation and specialization of essential structures, including notochord, muscle, heart, and differentiated neurons. In recent years, the experimental advantages offered by ascidian embryos have allowed the rapid accumulation of a wealth of information on the molecular mechanisms that regulate the expression of T-box genes. These studies have also elucidated the strategies employed by these transcription factors to orchestrate the appropriate spatial and temporal deployment of the numerous target genes that they control.

Ciona intestinalis as a marine model system to study some key developmental genes targeted by the diatom-derived aldehyde decadienal

The anti-proliferative effects of diatoms, described for the first time in copepods, have also been demonstrated in benthic invertebrates such as polychaetes, sea urchins and tunicates. In these organisms PUAs (polyunsaturated aldehydes) induce the disruption of gametogenesis, gamete functionality, fertilization, embryonic mitosis, and larval fitness and competence. These inhibitory effects are due to the PUAs, produced by diatoms in response to physical damage as occurs during copepod grazing. The cell targets of these compounds remain largely unknown. Here we identify some of the genes targeted by the diatom PUA 2-trans-4-trans-decadienal (DD) using the tunicate Ciona intestinalis. The tools, techniques and genomic resources available for Ciona, as well as the suitability of Ciona embryos for medium-to high-throughput strategies, are key to their employment as model organisms in different fields, including the investigation of toxic agents that could interfere with developmental processes. We demonstrate that DD can induce developmental aberrations in Ciona larvae in a dose-dependent manner. Moreover, through a preliminary analysis, DD is shown to affect the expression level of genes involved in stress response and developmental processes.

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.

The role of gene duplication and unconstrained selective pressures in the melanopsin gene family evolution and vertebrate circadian rhythm regulation

Melanopsin is a photosensitive cell protein involved in regulating circadian rhythms and other non-visual responses to light. The melanopsin gene family is represented by two paralogs, OPN4x and OPN4m, which originated through gene duplication early in the emergence of vertebrates. Here we studied the melanopsin gene family using an integrated gene/protein evolutionary approach, which revealed that the rhabdomeric urbilaterian ancestor had the same amino acid patterns (DRY motif and the Y and E conterions) as extant vertebrate species, suggesting that the mechanism for light detection and regulation is similar to rhabdomeric rhodopsins. Both OPN4m and OPN4x paralogs are found in vertebrate genomic paralogons, suggesting that they diverged following this duplication event about 600 million years ago, when the complex eye emerged in the vertebrate ancestor. Melanopsins generally evolved under negative selection (ω = 0.171) with some minor episodes of positive selection (proportion of sites = 25%) and functional divergence (θ(I) = 0.349 and θ(II) = 0.126). The OPN4m and OPN4x melanopsin paralogs show evidence of spectral divergence at sites likely involved in melanopsin light absorbance (200F, 273S and 276A). Also, following the teleost lineage-specific whole genome duplication (3R) that prompted the teleost fish radiation, type I divergence (θ(I) = 0.181) and positive selection (affecting 11% of sites) contributed to amino acid variability that we related with the photo-activation stability of melanopsin. The melanopsin intracellular regions had unexpectedly high variability in their coupling specificity of G-proteins and we propose that Gq/11 and Gi/o are the two G-proteins most-likely to mediate the melanopsin phototransduction pathway. The selection signatures were mainly observed on retinal-related sites and the third and second intracellular loops, demonstrating the physiological plasticity of the melanopsin protein group. Our results provide new insights on the phototransduction process and additional tools for disentangling and understanding the links between melanopsin gene evolution and the specializations observed in vertebrates, especially in teleost fish.

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 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.

Evolution of the insulin-like growth factor binding protein (IGFBP) family

The evolution of the IGF binding protein (IGFBP) gene family has been difficult to resolve. Both chromosomal and serial duplications have been suggested as mechanisms for the expansion of this gene family. We have identified and annotated IGFBP sequences from a wide selection of vertebrate species as well as Branchiostoma floridae and Ciona intestinalis. By combining detailed sequence analysis with sequence-based phylogenies and chromosome information, we arrive at the following scenario: the ancestral chordate IGFBP gene underwent a local gene duplication, resulting in a gene pair adjacent to a HOX cluster. Subsequently, the gene family expanded in the two basal vertebrate tetraploidization (2R) resulting in the six IGFBP types that are presently found in placental mammals. The teleost fish ancestor underwent a third tetraploidization (3R) that further expanded the IGFBP repertoire. The five sequenced teleost fish genomes retain 9-11 of IGFBP genes. This scenario is supported by the phylogenies of three adjacent gene families in the HOX gene regions, namely the epidermal growth factor receptors (EGFR) and the Ikaros and distal-less (DLX) transcription factors. Our sequence comparisons show that several important structural components in the IGFBPs are ancestral vertebrate features that have been maintained in all orthologs, for instance the integrin interaction motif Arg-Gly-Asp in IGFBP-2. In contrast, the Arg-Gly-Asp motif in IGFBP-1 has arisen independently in mammals. The large degree of retention of IGFBP genes after the ancient expansion of the gene family strongly suggests that each gene evolved distinct and important functions early in vertebrate evolution.

Regulation of the Caenorhabditis elegans posterior Hox gene egl-5 by microRNA and the polycomb-like gene sop-2

In Caenorhabditis elegans, the domains of Hox gene expression are controlled by the novel global regulatory gene sop-2. We identified a region located 3' of the Hox gene egl-5 that promotes ectopic expression of an egl-5 reporter gene in a sop-2 mutant. SOP-2 could directly block positive regulatory factors acting in this region, or it could block their expression. We identified three possible miRNA binding sites within the egl-5 3' untranslated region (UTR). Cognate microRNAs are expressed in relevant tissues and can block egl-5 expression when expressed from a transgene. Mutation of the putative binding sites in the egl-5 3'UTR resulted in a modest degree of misexpression of a minimal egl-5 reporter gene, suggesting that microRNAs may contribute to the tight restriction of egl-5 expression to particular cell lineages.

Evolutionary conservation of vertebrate notochord genes in the ascidian Ciona intestinalis

To reconstruct a minimum complement of notochord genes evolutionarily conserved across chordates, we scanned the Ciona intestinalis genome using the sequences of 182 genes reported to be expressed in the notochord of different vertebrates and identified 139 candidate notochord genes. For 66 of these Ciona genes expression data were already available, hence we analyzed the expression of the remaining 73 genes and found notochord expression for 20. The predicted products of the newly identified notochord genes range from the transcription factors Ci-XBPa and Ci-miER1 to extracellular matrix proteins. We examined the expression of the newly identified notochord genes in embryos ectopically expressing Ciona Brachyury (Ci-Bra) and in embryos expressing a repressor form of this transcription factor in the notochord, and we found that while a subset of the genes examined are clearly responsive to Ci-Bra, other genes are not affected by alterations in its levels. We provide a first description of notochord genes that are not evidently influenced by the ectopic expression of Ci-Bra and we propose alternative regulatory mechanisms that might control their transcription.

MicroRNAs in the Hox network: an apparent link to posterior prevalence

Homeobox (Hox) transcription factors confer anterior-posterior (AP) axial coordinates to vertebrate embryos. Hox genes are found in clusters that also contain genes for microRNAs (miRNAs). Our analysis of predicted miRNA targets indicates that Hox cluster-embedded miRNAs preferentially target Hox mRNAs. Moreover, the presumed Hox target genes are predominantly situated on the 3' side of each Hox miRNA locus. These results suggest that Hox miRNAs help repress more anterior programmes, thereby reinforcing posterior prevalence, which is the hierarchical dominance of posterior over anterior Hox gene function that is observed in bilaterians. In this way, miRNA-mediated regulation seems to recapitulate interactions at other levels of gene expression, some more ancestral, within a network under stabilizing selection.

Phylogenetic and chromosomal analyses of multiple gene families syntenic with vertebrate Hox clusters

Background

Ever since the theory about two rounds of genome duplication (2R) in the vertebrate lineage was proposed, the Hox gene clusters have served as the prime example of quadruplicate paralogy in mammalian genomes. In teleost fishes, the observation of additional Hox clusters absent in other vertebrate lineages suggested a third tetraploidization (3R). Because the Hox clusters occupy a quite limited part of each chromosome, and are special in having position-dependent regulation within the multi-gene cluster, studies of syntenic gene families are needed to determine the extent of the duplicated chromosome segments. We have analyzed in detail 14 gene families that are syntenic with the Hox clusters to see if their phylogenies are compatible with the Hox duplications and the 2R/3R scenario. Our starting point was the gene family for the NPY family of peptides located near the Hox clusters in the pufferfish Takifugu rubripes, the zebrafish Danio rerio, and human.

Results

Seven of the gene families have members on at least three of the human Hox chromosomes and two families are present on all four. Using both neighbor-joining and quartet-puzzling maximum likelihood methods we found that 13 families have a phylogeny that supports duplications coinciding with the Hox cluster duplications. One additional family also has a topology consistent with 2R but due to lack of urochordate or cephalocordate sequences the time window when these duplications could have occurred is wider. All but two gene families also show teleost-specific duplicates.

Conclusion

Based on this analysis we conclude that the Hox cluster duplications involved a large number of adjacent gene families, supporting expansion of these families in the 2R, as well as in the teleost 3R tetraploidization. The gene duplicates presumably provided raw material in early vertebrate evolution for neofunctionalization and subfunctionalization.

PCR survey of Xenoturbella bocki Hox genes

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

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.

Genome regulation by polycomb and trithorax proteins

Polycomb group (PcG) and trithorax group (trxG) proteins are critical regulators of numerous developmental genes. To silence or activate gene expression, respectively, PcG and trxG proteins bind to specific regions of DNA and direct the posttranslational modification of histones. Recent work suggests that PcG proteins regulate the nuclear organization of their target genes and that PcG-mediated gene silencing involves noncoding RNAs and the RNAi machinery.

Building divergent body plans with similar genetic pathways

Deuterostome animals exhibit widely divergent body plans. Echinoderms have either radial or bilateral symmetry, hemichordates include bilateral enteropneust worms and colonial pterobranchs, and chordates possess a defined dorsal-ventral axis imposed on their anterior-posterior axis. Tunicates are chordates only as larvae, following metamorphosis the adults acquire a body plan unique for the deuterostomes. This paper examines larval and adult body plans in the deuterostomes and discusses two distinct ways of evolving divergent body plans. First, echinoderms and hemichordates have similar feeding larvae, but build a new adult body within or around their larvae. In hemichordates and many direct-developing echinoderms, the adult is built onto the larva, with the larval axes becoming the adult axes and the larval mouth becoming the adult mouth. In contrast, indirect-developing echinoderms undergo radical metamorphosis where adult axes are not the same as larval axes. A second way of evolving a divergent body plan is to become colonial, as seen in hemichordates and tunicates. Early embryonic development and gastrulation are similar in all deuterostomes, but, in chordates, the anterior-posterior axis is established at right angles to the animal-vegetal axis, in contrast to hemichordates and indirect-developing echinoderms. Hox gene sequences and anterior-posterior expression patterns illuminate deuterostome phylogenetic relationships and the evolution of unique adult body plans within monophyletic groups. Many genes that are considered vertebrate 'mesodermal' genes, such as nodal and brachyury T, are likely to ancestrally have been involved in the formation of the mouth and anus, and later were evolutionarily co-opted into mesoderm during vertebrate development.

Hox genes are not always Colinear

The deuterostomes are the clade of animals for which we have the most detailed understanding of Hox cluster organisation. With the Hox cluster of amphioxus (Branchiostoma floridae) we have the best prototypical, least derived Hox cluster for the group, whilst the urochordates present us with some of the most highly derived and disintegrated clusters. Combined with the detailed mechanistic understanding of vertebrate Hox regulation, the deuterostomes provide much of the most useful data for understanding Hox cluster evolution. Considering both the prototypical and derived deuterostome Hox clusters leads us to hypothesize that Temporal Colinearity is the main constraining force on Hox cluster organisation, but until we have a much deeper understanding of the mechanistic basis for this phenomenon, and know how widespread across the Bilateria the mechanism(s) is/are, then we cannot know how the Hox cluster of the last common bilaterian operated and what have been the major evolutionary forces operating upon the Hox gene cluster.

Cataloging transcription factor and major signaling molecule genes for functional genomic studies in Ciona intestinalis

The ascidian Ciona intestinalis provides an excellent experimental system for functional genomic studies because (1) its genome has been sequenced, (2) the transcription factor genes and genes for major signal transduction molecules have been extensively screened and annotated on a genome-wide scale using the molecular phylogenetical method, and (3) their embryonic expression profiles have been almost completely determined. However, the entire genetic structure, including the 5' and 3' untranslated regions and the protein-coding regions, of most gene models used in these prior studies is not always supported by cDNA evidence, and thus, these gene models are potentially imprecise. To facilitate functional genomic studies based on precise gene structures, our present study determined 406 cDNA sequences for 357 transcription factor genes and 112 cDNA sequences for 107 signal transduction molecule genes, greatly improving the previous gene models and revealing transcript variants for 44 genes. Considering these data alongside those of previously characterized genes deposited in the DNA Data Bank of Japan/European Molecular Biology Laboratory/GENBANK databases, 95.6% of the catalogued transcription factor genes (373/390) and 98.3% of the catalogued signal transduction molecule genes (117/119) have now been verified by cDNA sequences. Thus, the present study greatly improves the resources available for functional genomic studies in C. intestinalis.

Ciona intestinalis: chordate development made simple

Thanks to their transparent and rapidly developing mosaic embryos, ascidians (or sea squirts) have been a model system for embryological studies for over a century. Recently, ascidians have entered the postgenomic era, with the sequencing of the Ciona intestinalis genome and the accumulation of molecular resources that rival those available for fruit flies and mice. One strength of ascidians as a model system is their close similarity to vertebrates. Literature reporting molecular homologies between vertebrate and ascidian tissues has flourished over the past 15 years, since the first ascidian genes were cloned. However, it should not be forgotten that ascidians diverged from the lineage leading to vertebrates over 500 million years ago. Here, we review the main similarities and differences so far identified, at the molecular level, between ascidian and vertebrate tissues and discuss the evolution of the compact ascidian genome.

Consequences of the Evolution of the GABAA Receptor Gene Family

1. This paper reviews the evolution of the family of genes present in mammals and other vertebrates that encode gamma-aminobutyric acid (GABA) type A (GABA(A)) receptors, which are the major inhibitory neurotransmitter receptors in the central nervous system (CNS). In mammals, 16 different polypeptides (alpha1-alpha6, beta1-beta3, gamma1-gamma3, delta, epsilon, pi, and theta) have been identified, using recombinant DNA techniques, each of which is encoded by a distinct gene. The products of these genes assemble in diverse combinations to form a variety of receptor subtypes that have different sensitivities to a number of clinically relevant compounds, such as the benzodiazepines (BZs). 2. Based on a number of chromosomal mapping techniques, the majority of the GABA(A) receptor genes have been localized, in man, in four clusters on chromosomes 4, 5, 15, and the X. Furthermore, the genes that are present within these clusters have a conserved transcriptional orientation. It has, therefore, been proposed that the clusters arose largely as a consequence of two whole-genome doublings that occurred during chordate evolution, and that the ancestral cluster contained an "alpha-like," a "beta-like," and a "gamma-like" subunit gene. 3. Our laboratory has identified two additional GABA(A) receptor polypeptides (the beta4 and gamma4 subunits) in a number of vertebrate species; these do not appear to be present in mammals. We discuss here the relationship of the corresponding genes to other GABA(A) receptor genes, and conclude that their products are orthologous to the mammalian theta and epsilon subunits, respectively. 4. The GABA(A) receptor has a number of binding sites for compounds such as BZs, barbiturates, neurosteroids, and certain volatile anaesthetics. However, the only site at which endogenous compounds are thought to be active is the steroid site; this binds steroids such as certain metabolites of progesterone and deoxycorticosterone, which are synthesized in the periphery and CNS. Since the in vivo functional relevance, if any, of binding sites for other classes of compounds (such as the BZs) is unknown, the significance of differences in primary sequence, between different receptor subunits, is uncertain. We suggest that a possibly more important consequence of gene duplication is that it permitted greater flexibility in the level, pattern and regulation of expression of GABA(A) receptor genes.

Remodelling of the homeobox gene complement in the tunicate Oikopleura dioica

Homeodomain transcription factors are involved in many developmental processes and have been intensely studied in a few model organisms, such as mouse, Drosophila and Caenorhabditis elegans. Homeobox genes fall into 10 classes (ANTP, PRD, POU, LIM, TALE, SIX, Cut, ZFH, HNF1, Prox) and 89 different families/groups, all of which are present in vertebrates. Additional groups may be uncovered by further genome annotation, particularly of complex vertebrate genomes. Eight of these groups have been found only in vertebrates, but not in the genome of the tunicate Ciona intestinalis. The other 81 groups of homeobox gene that have been detected in vertebrates so far probably appeared during the early evolution of bilaterians or earlier, as they are also present outside the chordates. How the homeobox genes evolved during and after the main radiation of the bilaterians remains poorly understood, as only a few animal genomes have been sequenced completely. However, drastic changes have occurred at least in the lineage of C. elegans , such as loss of several Hox genes and Hox cluster fragmentation . Here we report considerable alterations of the homeobox gene complement in the tunicate lineage.

A saturation screen for cis-acting regulatory DNA in the Hox genes of Ciona intestinalis

A screen for the systematic identification of cis-regulatory elements within large (>100 kb) genomic domains containing Hox genes was performed by using the basal chordate Ciona intestinalis. Randomly generated DNA fragments from bacterial artificial chromosomes containing two clusters of Hox genes were inserted into a vector upstream of a minimal promoter and lacZ reporter gene. A total of 222 resultant fusion genes were separately electroporated into fertilized eggs, and their regulatory activities were monitored in larvae. In sum, 21 separable cis-regulatory elements were found. These include eight Hox linked domains that drive expression in nested anterior-posterior domains of ectodermally derived tissues. In addition to vertebrate-like CNS regulation, the discovery of cis-regulatory domains that drive epidermal transcription suggests that C. intestinalis has arthropod-like Hox patterning in the epidermis.

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.

Hox, ParaHox, ProtoHox: facts and guesses

The Hox gene cluster has captivated the imagination of evolutionary and developmental biologists worldwide. In this review, the origin of the Hox and ParaHox gene clusters by duplication of a ProtoHox gene cluster, and the changes in their gene numbers in major Metazoan Transitions are reviewed critically. Re-evaluation of existing data and recent findings in Cnidarians, Acoels, and critical stages of vertebrate evolution suggest alternative scenarios for the origin, structure, and changes in Hox gene numbers in relevant events of Metazoan evolution. I discuss opposing views and propose that (i) the ProtoHox cluster had only two genes, and not four as commonly believed: a corollary is that the origin of Bilaterians was coincident with the invention of new Hox and ParaHox gene classes, which may have facilitated such a transition; (ii) the ProtoHox cluster duplication was a cis duplication event, rather than a trans duplication event, as previously suggested, and (iii) the ancestral vertebrate cluster possessed 14 Hox genes, and not the 13 generally assumed. These hypotheses could be verified or refuted in the near future, but they may help critical discussion of the evolution of the Hox/ParaHox family in the metazoan kingdom.

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

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

Hox cluster disintegration with persistent anteroposterior order of expression in Oikopleura dioica

Tunicate embryos and larvae have small cell numbers and simple anatomical features in comparison with other chordates, including vertebrates. Although they branch near the base of chordate phylogenetic trees, their degree of divergence from the common chordate ancestor remains difficult to evaluate. Here we show that the tunicate Oikopleura dioica has a complement of nine Hox genes in which all central genes are lacking but a full vertebrate-like set of posterior genes is present. In contrast to all bilaterians studied so far, Hox genes are not clustered in the Oikopleura genome. Their expression occurs mostly in the tail, with some tissue preference, and a strong partition of expression domains in the nerve cord, in the notochord and in the muscle. In each tissue of the tail, the anteroposterior order of Hox gene expression evokes spatial collinearity, with several alterations. We propose a relationship between the Hox cluster breakdown, the separation of Hox expression domains, and a transition to a determinative mode of development.

Exclusion of repetitive DNA elements from gnathostomeHox clusters

Despite their homology and analogous function, the Hox gene clusters of vertebrates and invertebrates are subject to different constraints on their structural organization. This is demonstrated by a drastically different distribution of repetitive DNA elements in the Hox cluster regions. While gnathostomes have a strong tendency to exclude repetitive DNA elements from the inside of their Hox clusters, no such trend can be detected in the Hox gene clusters of protostomes. Repeats "invade" the gnathostome Hox clusters from the 5' and 3' ends while the core of the clusters remains virtually free of repetitive DNA. This invasion appears to be correlated with relaxed constraints associated with gene loss after cluster duplications.

Chromosomal mapping of ANTP class homeobox genes in amphioxus: piecing together ancestral genomes

Homeobox genes encode DNA-binding proteins, many of which are implicated in the control of embryonic development. Evolutionarily, most homeobox genes fall into two related clades: the ANTP and the PRD classes. Some genes in ANTP class, notably Hox, ParaHox, and NK genes, have an intriguing arrangement into physical clusters. To investigate the evolutionary history of these gene clusters, we examined homeobox gene chromosomal locations in the cephalochordate amphioxus, Branchiostoma floridae. We deduce that 22 amphioxus ANTP class homeobox genes localize in just three chromosomes. One contains the Hox cluster plus AmphiEn, AmphiMnx, and AmphiDll. The ParaHox cluster resides in another chromosome, whereas a third chromosome contains the NK type homeobox genes, including AmphiMsx and AmphiTlx. By comparative analysis we infer that clustering of ANTP class homeobox genes evolved just once, during a series of extensive cis-duplication events of genes early in animal evolution. A trans-duplication event occurred later to yield the Hox and ParaHox gene clusters on different chromosomes. The results obtained have implications for understanding the origin of homeobox gene clustering, the diversification of the ANTP class of homeobox genes, and the evolution of animal genomes.