Large scale transgenic and cluster deletion analysis of the HoxD complex separate an ancestral regulatory module from evolutionary innovations

EVOLUTIONARY CO-OPTION OF AN ANCESTRAL CLOACAL REGULATORY LANDSCAPE DURING THE EMERGENCE OF DIGITS AND GENITALS

The transition from fins to limbs has been a rich source of discussion for more than a century. One open and important issue is understanding how the mechanisms that pattern digits arose during vertebrate evolution. In this context, the analysis of Hox gene expression and functions to infer evolutionary scenarios has been a productive approach to explain the changes in organ formation, particularly in limbs. In tetrapods, the transcription of Hoxd genes in developing digits depends on a well-characterized set of enhancers forming a large regulatory landscape1,2. This control system has a syntenic counterpart in zebrafish, even though they lack bona fide digits, suggestive of deep homology3 between distal fin and limb developmental mechanisms. We tested the global function of this landscape to assess ancestry and source of limb and fin variation. In contrast to results in mice, we show here that the deletion of the homologous control region in zebrafish has a limited effect on the transcription of hoxd genes during fin development. However, it fully abrogates hoxd expression within the developing cloaca, an ancestral structure related to the mammalian urogenital sinus. We show that similar to the limb, Hoxd gene function in the urogenital sinus of the mouse also depends on enhancers located in this same genomic domain. Thus, we conclude that the current regulation underlying Hoxd gene expression in distal limbs was co-opted in tetrapods from a preexisting cloacal program. The orthologous chromatin domain in fishes may illustrate a rudimentary or partial step in this evolutionary co-option.

Origin matters: Using a local reference genome improves measures in population genomics

Genome sequencing enables answering fundamental questions about the genetic basis of adaptation, population structure and epigenetic mechanisms. Yet, we usually need a suitable reference genome for mapping population-level resequencing data. In some model systems, multiple reference genomes are available, giving the challenging task of determining which reference genome best suits the data. Here, we compared the use of two different reference genomes for the three-spined stickleback (Gasterosteus aculeatus), one novel genome derived from a European gynogenetic individual and the published reference genome of a North American individual. Specifically, we investigated the impact of using a local reference versus one generated from a distinct lineage on several common population genomics analyses. Through mapping genome resequencing data of 60 sticklebacks from across Europe and North America, we demonstrate that genetic distance among samples and the reference genomes impacts downstream analyses. Using a local reference genome increased mapping efficiency and genotyping accuracy, effectively retaining more and better data. Despite comparable distributions of the metrics generated across the genome using SNP data (i.e. π, Tajima's D and FST ), window-based statistics using different references resulted in different outlier genes and enriched gene functions. A marker-based analysis of DNA methylation distributions had a comparably high overlap in outlier genes and functions, yet with distinct differences depending on the reference genome. Overall, our results highlight how using a local reference genome decreases reference bias to increase confidence in downstream analyses of the data. Such results have significant implications in all reference-genome-based population genomic analyses.

A brain-wide analysis maps structural evolution to distinct anatomical module

The vertebrate brain is highly conserved topologically, but less is known about neuroanatomical variation between individual brain regions. Neuroanatomical variation at the regional level is hypothesized to provide functional expansion, building upon ancestral anatomy needed for basic functions. Classically, animal models used to study evolution have lacked tools for detailed anatomical analysis that are widely used in zebrafish and mice, presenting a barrier to studying brain evolution at fine scales. In this study, we sought to investigate the evolution of brain anatomy using a single species of fish consisting of divergent surface and cave morphs, that permits functional genetic testing of regional volume and shape across the entire brain. We generated a high-resolution brain atlas for the blind Mexican cavefish Astyanax mexicanus and coupled the atlas with automated computational tools to directly assess variability in brain region shape and volume across all populations. We measured the volume and shape of every grossly defined neuroanatomical region of the brain and assessed correlations between anatomical regions in surface fish, cavefish, and surface × cave F2 hybrids, whose phenotypes span the range of surface to cave. We find that dorsal regions of the brain are contracted, while ventral regions have expanded, with F2 hybrid data providing support for developmental constraint along the dorsal-ventral axis. Furthermore, these dorsal-ventral relationships in anatomical variation show similar patterns for both volume and shape, suggesting that the anatomical evolution captured by these two parameters could be driven by similar developmental mechanisms. Together, these data demonstrate that A. mexicanus is a powerful system for functionally determining basic principles of brain evolution and will permit testing how genes influence early patterning events to drive brain-wide anatomical evolution.

Synthetic regulatory reconstitution reveals principles of mammalian Hox cluster regulation

Precise Hox gene expression is crucial for embryonic patterning. Intra-Hox transcription factor binding and distal enhancer elements have emerged as the major regulatory modules controlling Hox gene expression. However, quantifying their relative contributions has remained elusive. Here, we introduce "synthetic regulatory reconstitution," a conceptual framework for studying gene regulation, and apply it to the HoxA cluster. We synthesized and delivered variant rat HoxA clusters (130 to 170 kilobases) to an ectopic location in the mouse genome. We found that a minimal HoxA cluster recapitulated correct patterns of chromatin remodeling and transcription in response to patterning signals, whereas the addition of distal enhancers was needed for full transcriptional output. Synthetic regulatory reconstitution could provide a generalizable strategy for deciphering the regulatory logic of gene expression in complex genomes.

Developmental and evolutionary comparative analysis of a regulatory landscape in mouse and chicken

Modifications in gene regulation are driving forces in the evolution of organisms. Part of these changes involve cis-regulatory elements (CREs), which contact their target genes through higher-order chromatin structures. However, how such architectures and variations in CREs contribute to transcriptional evolvability remains elusive. We use Hoxd genes as a paradigm for the emergence of regulatory innovations, as many relevant enhancers are located in a regulatory landscape highly conserved in amniotes. Here, we analysed their regulation in murine vibrissae and chicken feather primordia, two skin appendages expressing different Hoxd gene subsets, and compared the regulation of these genes in these appendages with that in the elongation of the posterior trunk. In the two former structures, distinct subsets of Hoxd genes are contacted by different lineage-specific enhancers, probably as a result of using an ancestral chromatin topology as an evolutionary playground, whereas the gene regulation that occurs in the mouse and chicken embryonic trunk partially relies on conserved CREs. A high proportion of these non-coding sequences active in the trunk have functionally diverged between species, suggesting that transcriptional robustness is maintained, despite considerable divergence in enhancer sequences.

Summary: Analyses of the relationships between chromatin architecture and regulatory activities at the HoxD locus show that ancestral transcription patterns can be maintained while new regulations evolve.

Diversification and Functional Evolution of HOX Proteins

Gene duplication and divergence is a major contributor to the generation of morphological diversity and the emergence of novel features in vertebrates during evolution. The availability of sequenced genomes has facilitated our understanding of the evolution of genes and regulatory elements. However, progress in understanding conservation and divergence in the function of proteins has been slow and mainly assessed by comparing protein sequences in combination with in vitro analyses. These approaches help to classify proteins into different families and sub-families, such as distinct types of transcription factors, but how protein function varies within a gene family is less well understood. Some studies have explored the functional evolution of closely related proteins and important insights have begun to emerge. In this review, we will provide a general overview of gene duplication and functional divergence and then focus on the functional evolution of HOX proteins to illustrate evolutionary changes underlying diversification and their role in animal evolution.

A CYC-RAD-DIV-DRIF interaction likely pre-dates the origin of floral monosymmetry in Lamiales

BACKGROUND

An outstanding question in evolutionary biology is how genetic interactions defining novel traits evolve. They may evolve either by de novo assembly of previously non-interacting genes or by en bloc co-option of interactions from other functions. We tested these hypotheses in the context of a novel phenotype-Lamiales flower monosymmetry-defined by a developmental program that relies on regulatory interaction among CYCLOIDEA, RADIALIS, DIVARICATA, and DRIF gene products. In Antirrhinum majus (snapdragon), representing Lamiales, we tested whether components of this program likely function beyond their previously known role in petal and stamen development. In Solanum lycopersicum (tomato), representing Solanales which diverged from Lamiales before the origin of Lamiales floral monosymmetry, we additionally tested for regulatory interactions in this program.

RESULTS

We found that RADIALIS, DIVARICATA, and DRIF are expressed in snapdragon ovaries and developing fruit, similar to their homologs during tomato fruit development. In addition, we found that a tomato CYCLOIDEA ortholog positively regulates a tomato RADIALIS ortholog.

CONCLUSION

Our results provide preliminary support to the hypothesis that the developmental program defining floral monosymmetry in Lamiales was co-opted en bloc from a function in carpel development. This expands our understanding of novel trait evolution facilitated by co-option of existing regulatory interactions.

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.

An atlas of seven zebrafish hox cluster mutants provides insights into sub/neofunctionalization of vertebrate Hox clusters

Vertebrate Hox clusters are comprised of multiple Hox genes that control morphology and developmental timing along multiple body axes. Although results of genetic analyses using Hox-knockout mice have been accumulating, genetic studies in other vertebrates have not been sufficient for functional comparisons of vertebrate Hox genes. In this study, we isolated all of the seven hox cluster loss-of-function alleles in zebrafish using the CRISPR-Cas9 system. Comprehensive analysis of the embryonic phenotype and X-ray micro-computed tomography scan analysis of adult fish revealed several species-specific functional contributions of homologous Hox clusters along the appendicular axis, whereas important shared general principles were also confirmed, as exemplified by serial anterior vertebral transformations along the main body axis, observed in fish for the first time. Our results provide insights into discrete sub/neofunctionalization of vertebrate Hox clusters after quadruplication of the ancient Hox cluster. This set of seven complete hox cluster loss-of-function alleles provide a formidable resource for future developmental genetic analysis of the Hox patterning system in zebrafish.

Dbx2 regulation in limbs suggests interTAD sharing of enhancers

BACKGROUND

During tetrapod limb development, the HOXA13 and HOXD13 transcription factors are critical for the emergence and organization of the autopod, the most distal aspect where digits will develop. Since previous work had suggested that the Dbx2 gene is a target of these factors, we set up to analyze in detail this potential regulatory interaction.

RESULTS

We show that HOX13 proteins bind to mammalian-specific sequences at the vicinity of the Dbx2 locus that have enhancer activity in developing digits. However, the functional inactivation of the DBX2 protein did not elicit any particular phenotype related to Hox genes inactivation in digits, suggesting either redundant or compensatory mechanisms. We report that the neighboring Nell2 and Ano6 genes are also expressed in distal limb buds and are in part controlled by the same Dbx2 enhancers despite being localized into two different topologically associating domains (TADs) flanking the Dbx2 locus.

CONCLUSIONS

We conclude that Hoxa13 and Hoxd genes cooperatively activate Dbx2 expression in developing digits through binding to mammalian specific regulatory sequences in the Dbx2 neighborhood. Furthermore, these enhancers can overcome TAD boundaries in either direction to co-regulate a set of genes located in distinct chromatin domains.

HOX gene cluster (de)regulation in brain: from neurodevelopment to malignant glial tumours

HOX genes encode a family of evolutionarily conserved homeodomain transcription factors that are crucial both during development and adult life. In humans, 39 HOX genes are arranged in four clusters (HOXA, B, C, and D) in chromosomes 7, 17, 12, and 2, respectively. During embryonic development, particular epigenetic states accompany their expression along the anterior-posterior body axis. This tightly regulated temporal-spatial expression pattern reflects their relative chromosomal localization, and is critical for normal embryonic brain development when HOX genes are mainly expressed in the hindbrain and mostly absent in the forebrain region. Epigenetic marks, mostly polycomb-associated, are dynamically regulated at HOX loci and regulatory regions to ensure the finely tuned HOX activation and repression, highlighting a crucial epigenetic plasticity necessary for homeostatic development. HOX genes are essentially absent in healthy adult brain, whereas they are detected in malignant brain tumours, namely gliomas, where HOX genes display critical roles by regulating several hallmarks of cancer. Here, we review the major mechanisms involved in HOX genes (de)regulation in the brain, from embryonic to adult stages, in physiological and oncologic conditions. We focus particularly on the emerging causes of HOX gene deregulation in glioma, as well as on their functional and clinical implications.

Distal Limb Patterning Requires Modulation of cis-Regulatory Activities by HOX13

The combinatorial expression of Hox genes along the body axes is a major determinant of cell fate and plays a pivotal role in generating the animal body plan. Loss of HOXA13 and HOXD13 transcription factors (HOX13) leads to digit agenesis in mice, but how HOX13 proteins regulate transcriptional outcomes and confer identity to the distal-most limb cells has remained elusive. Here, we report on the genome-wide profiling of HOXA13 and HOXD13 in vivo binding and changes of the transcriptome and chromatin state in the transition from the early to the late-distal limb developmental program, as well as in Hoxa13^−/−; Hoxd13^−/−limbs. Our results show that proper termination of the early limb transcriptional program and activation of the late-distal limb program are coordinated by the dual action of HOX13 on cis-regulatory modules.

The HoxD cluster is a dynamic and resilient TAD boundary controlling the segregation of antagonistic regulatory landscapes

The mammalian HoxD cluster lies between two topologically associating domains (TADs) matching distinct enhancer-rich regulatory landscapes. During limb development, the telomeric TAD controls the early transcription of Hoxd genes in forearm cells, whereas the centromeric TAD subsequently regulates more posterior Hoxd genes in digit cells. Therefore, the TAD boundary prevents the terminal Hoxd13 gene from responding to forearm enhancers, thereby allowing proper limb patterning. To assess the nature and function of this CTCF-rich DNA region in embryos, we compared chromatin interaction profiles between proximal and distal limb bud cells isolated from mutant stocks where various parts of this boundary region were removed. The resulting progressive release in boundary effect triggered inter-TAD contacts, favored by the activity of the newly accessed enhancers. However, the boundary was highly resilient, and only a 400-kb deletion, including the whole-gene cluster, was eventually able to merge the neighboring TADs into a single structure. In this unified TAD, both proximal and distal limb enhancers nevertheless continued to work independently over a targeted transgenic reporter construct. We propose that the whole HoxD cluster is a dynamic TAD border and that the exact boundary position varies depending on both the transcriptional status and the developmental context.

Control of growth and gut maturation by HoxD genes and the associated lncRNA Haglr

During embryonic development, Hox genes participate in the building of a functional digestive system in metazoans, and genetic conditions involving these genes lead to important, sometimes lethal, growth retardation. Recently, this phenotype was obtained after deletion of Haglr, the Hoxd antisense growth-associated long noncoding RNA (lncRNA) located between Hoxd1 and Hoxd3 In this study, we have analyzed the function of Hoxd genes in delayed growth trajectories by looking at several nested targeted deficiencies of the mouse HoxD cluster. Mutant pups were severely stunted during the suckling period, but many recovered after weaning. After comparing seven distinct HoxD alleles, including CRISPR/Cas9 deletions involving Haglr, we identified Hoxd3 as the critical component for the gut to maintain milk-digestive competence. This essential function could be abrogated by the dominant-negative effect of HOXD10 as shown by a genetic rescue approach, thus further illustrating the importance of posterior prevalence in Hox gene function. A role for the lncRNA Haglr in the control of postnatal growth could not be corroborated.

The mitogen-activated protein kinase pathway in melanoma part I - Activation and primary resistance mechanisms to BRAF inhibition

Mitogen-activated protein kinase (MAPK) pathway has an important role in normal cells and can be activated under physiological conditions. MAPK pathway activation is a fundamental step in several intracellular processes requiring a sequential phosphorylation of the different pathway components. In normal cells, when MAPK pathway activation occurs, it leads to cell growth and differentiation. In order to prevent persistent MAPK pathway activation, physiological upstream negative feedback also takes place. In cells harbouring BRAFV600 mutations, the process leading to MAPK pathway activation is different, and the negative physiological feedback does not exist thus leading to permanent MAPK pathway activation, which ultimately can lead to uncontrolled proliferation. Targeted therapy with rapidly accelerated fibrosarcoma - B (BRAF) and/or mitogen-activated extracellular signal-regulated kinase kinase (MEK) inhibitors is indicated in patients with metastatic melanoma harboring BRAFV600 mutations. However, several different resistance mechanisms to this therapy were identified. In this review, we focus on primary or intrinsic resistance mechanisms to BRAF and MEK inhibition. In this setting, although a BRAF mutation is identified, there is no response to treatment with either BRAF or MEK inhibitor.

Reorganisation of Hoxd regulatory landscapes during the evolution of a snake-like body plan

Within land vertebrate species, snakes display extreme variations in their body plan, characterized by the absence of limbs and an elongated morphology. Such a particular interpretation of the basic vertebrate body architecture has often been associated with changes in the function or regulation of Hox genes. Here, we use an interspecies comparative approach to investigate different regulatory aspects at the snake HoxD locus. We report that, unlike in other vertebrates, snake mesoderm-specific enhancers are mostly located within the HoxD cluster itself rather than outside. In addition, despite both the absence of limbs and an altered Hoxd gene regulation in external genitalia, the limb-associated bimodal HoxD chromatin structure is maintained at the snake locus. Finally, we show that snake and mouse orthologous enhancer sequences can display distinct expression specificities. These results show that vertebrate morphological evolution likely involved extensive reorganisation at Hox loci, yet within a generally conserved regulatory framework.

Gene regulation at a distance: From remote enhancers to 3D regulatory ensembles

Large-scale identification of elements associated with gene expression revealed that many of them are located extremely far from gene transcriptional start sites. We review here the growing evidence that show that distal cis-acting elements provide key instructions to genes, as genetic variation affecting them is growingly identified as an importance source of phenotypic diversity and disease. We discuss the different mechanisms that allow these elements to exert their regulatory functions, in a robust and specific manner, despite the large genomic distances separating them from their target genes. We particularly focus on the role of the structural organization of the genome in guiding such regulatory interactions.

A genome-wide association study identifies a genomic region for the polycerate phenotype in sheep (Ovis aries)

Horns are a cranial appendage found exclusively in Bovidae, and play important roles in accessing resources and mates. In sheep (Ovies aries), horns vary from polled to six-horned, and human have been selecting polled animals in farming and breeding. Here, we conducted a genome-wide association study on 24 two-horned versus 22 four-horned phenotypes in a native Chinese breed of Sishui Fur sheep. Together with linkage disequilibrium (LD) analyses and haplotype-based association tests, we identified a genomic region comprising 132.0–133.1 Mb on chromosome 2 that contained the top 10 SNPs (including 4 significant SNPs) and 5 most significant haplotypes associated with the polycerate phenotype. In humans and mice, this genomic region contains the HOXD gene cluster and adjacent functional genes EVX2 and KIAA1715, which have a close association with the formation of limbs and genital buds. Our results provide new insights into the genetic basis underlying variable numbers of horns and represent a new resource for use in sheep genetics and breeding.

Independent regulation of vertebral number and vertebral identity by microRNA-196 paralogs

The Hox genes play a central role in patterning the embryonic anterior-to-posterior axis. An important function of Hox activity in vertebrates is the specification of different vertebral morphologies, with an additional role in axis elongation emerging. The miR-196 family of microRNAs (miRNAs) are predicted to extensively target Hox 3' UTRs, although the full extent to which miR-196 regulates Hox expression dynamics and influences mammalian development remains to be elucidated. Here we used an extensive allelic series of mouse knockouts to show that the miR-196 family of miRNAs is essential both for properly patterning vertebral identity at different axial levels and for modulating the total number of vertebrae. All three miR-196 paralogs, 196a1, 196a2, and 196b, act redundantly to pattern the midthoracic region, whereas 196a2 and 196b have an additive role in controlling the number of rib-bearing vertebra and positioning of the sacrum. Independent of this, 196a1, 196a2, and 196b act redundantly to constrain total vertebral number. Loss of miR-196 leads to a collective up-regulation of numerous trunk Hox target genes with a concomitant delay in activation of caudal Hox genes, which are proposed to signal the end of axis extension. Additionally, we identified altered molecular signatures associated with the Wnt, Fgf, and Notch/segmentation pathways and demonstrate that miR-196 has the potential to regulate Wnt activity by multiple mechanisms. By feeding into, and thereby integrating, multiple genetic networks controlling vertebral number and identity, miR-196 is a critical player defining axial formulae.

Epigenetic mechanisms and boundaries in the regulation of mammalian Hox clusters

Hox gene expression imparts segment identity to body structures along the anterior-posterior axis and is tightly governed by higher order chromatin mechanisms. Chromatin regulatory features of the homeotic complex are best defined in Drosophila melanogaster, where multiple cis-regulatory elements have been identified that ensure collinear Hox gene expression patterns in accordance with their genomic organization. Recent studies focused on delineating the epigenetic features of the vertebrate Hox clusters have helped reveal their dynamic chromatin organization and its impact on gene expression. Enrichment for the 'activating' H3K4me3 and 'repressive' H3K27me3 histone modifications is a particularly strong read-out for transcriptional status and correlates well with the evidence for chromatin loop domain structures and stage specific topological changes at these loci. However, it is not clear how such distinct domains are imposed and regulated independent of each other. Comparative analysis of the chromatin structure and organization of the homeotic gene clusters in fly and mammals is increasingly revealing the functional conservation of chromatin mediated mechanisms. Here we discuss the case for interspersed boundary elements existing within mammalian Hox clusters along with their possible roles and mechanisms of action. Recent studies suggest a role for factors other than the well characterized vertebrate boundary factor CTCF, such as the GAGA binding factor (GAF), in maintaining chromatin domains at the Hox loci. We also present data demonstrating how such regulatory elements may be involved in organizing higher order structure and demarcating active domains of gene expression at the mammalian Hox clusters.

Overcoming resistance to BRAF inhibition in BRAF-mutated metastatic melanoma

Almost 50% of metastatic melanoma patients harbor a BRAF(V600) mutation and the introduction of BRAF inhibitors has improved their treatment options. BRAF inhibitors vemurafenib and dabrafenib achieved improved overall survival over chemotherapy and have been approved for the treatment of BRAF-mutated metastatic melanoma. However, most patients develop mechanisms of acquired resistance and about 15% of them do not achieve tumor regression at all, due to intrinsic resistance to therapy. Moreover, early adaptive responses limit the initial efficacy of BRAF inhibition, leading mostly to incomplete responses that may favor the selection of a sub-population of resistant clones and the acquisition of alterations that cause tumor regrowth and progressive disease. The purpose of this paper is to review the mechanisms of resistance to therapy with BRAF inhibitors and to discuss the strategies to overcome them based on pre-clinical and clinical evidences.

Structure, function and evolution of topologically associating domains (TADs) at HOX loci

Hox genes encode transcription factors necessary for patterning the major developing anterior to posterior embryonic axis. In addition, during vertebrate evolution, various subsets of this gene family were co-opted along with the emergence of novel body structures, such as the limbs or the external genitalia. The morphogenesis of these axial structures thus relies in part upon the precisely controlled transcription of specific Hox genes, a mechanism involving multiple long-range enhancers. Recently, it was reported that such regulatory mechanisms were largely shared between different developing tissues, though with some specificities, suggesting the recruitment of ancestral regulatory modalities from one tissue to another. The analysis of chromatin architectures at HoxD and HoxA loci revealed the existence of two flanking topologically associating domains (TADs), precisely encompassing the adjacent regulatory landscapes. Here, we discuss the function of these TADs in the control of Hox gene regulation and we speculate about their capacity to serve as structural frameworks for the emergence of novel enhancers. In this view, TADs may have been used as genomic niches to evolve pleiotropic regulations found at many developmental loci.

Hox genes control vertebrate body elongation by collinear Wnt repression

In vertebrates, the total number of vertebrae is precisely defined. Vertebrae derive from embryonic somites that are continuously produced posteriorly from the presomitic mesoderm (PSM) during body formation. We show that in the chicken embryo, activation of posterior Hox genes (paralogs 9–13) in the tail-bud correlates with the slowing down of axis elongation. Our data indicate that a subset of progressively more posterior Hox genes, which are collinearly activated in vertebral precursors, repress Wnt activity with increasing strength. This leads to a graded repression of the Brachyury/T transcription factor, reducing mesoderm ingression and slowing down the elongation process. Due to the continuation of somite formation, this mechanism leads to the progressive reduction of PSM size. This ultimately brings the retinoic acid (RA)-producing segmented region in close vicinity to the tail bud, potentially accounting for the termination of segmentation and axis elongation.

DOI:http://dx.doi.org/10.7554/eLife.04379.001

In humans and other vertebrates, the number of bones (vertebrae) in the spine is determined early in development. The vertebrae form from blocks of tissue called somites that make segments along the body axis—a virtual line running from the head to the tail-end—of the embryo. The somites form as the embryo increases in length, with new somites forming periodically at the back near the embryo's tail-end.

A family of genes called the Hox genes are involved in controlling the formation of the somites. However, it is not known whether they directly control the number of somites that form, or whether they control the length of the body of the embryo.

Denans et al. studied the Hox genes in chicken embryos. The experiments suggest that the activation of some of the Hox genes in a structure called the tail-bud, which is found at the tail-end of the embryo, slow down the elongation of the body. The Hox genes achieve this by repressing the activity of a signaling pathway called Wnt so that Wnt activity in the tail-bud progressively decreases as the embryo develops.

The elongation of the body stops when the levels of a molecule called retinoic acid increase in the tail-bud, which causes the loss of the stem cells that are needed to make the somites. Denans et al.'s findings suggest that Hox genes influence the timing of the halt in elongation, which in turn is important for determining the total number of somites that form. Understanding how Hox genes control the formation of the cells that will make up the somites and influence Wnt signaling is a major challenge for the future.

DOI:http://dx.doi.org/10.7554/eLife.04379.002

YAP regulates the expression of Hoxa1 and Hoxc13 in mouse and human oral and skin epithelial tissues

Yes-associated protein (YAP) is a Hippo signaling transcriptional coactivator that plays pivotal roles in stem cell proliferation, organ size control, and tumor development. The downstream targets of YAP have been shown to be highly context dependent. In this study, we used the embryonic mouse tooth germ as a tool to search for the downstream targets of YAP in ectoderm-derived tissues. Yap deficiency in the dental epithelium resulted in a small tooth germ with reduced epithelial cell proliferation. We compared the gene expression profiles of embryonic day 14.5 (E14.5) Yap conditional knockout and YAP transgenic mouse tooth germs using transcriptome sequencing (RNA-Seq) and further confirmed the differentially expressed genes using real-time PCR and in situ hybridization. We found that YAP regulates the expression of Hoxa1 and Hoxc13 in oral and dental epithelial tissues as well as in the epidermis of skin during embryonic and adult stages. Sphere formation assay suggested that Hoxa1 and Hoxc13 are functionally involved in YAP-regulated epithelial progenitor cell proliferation, and chromatin immunoprecipitation (ChIP) assay implies that YAP may regulate Hoxa1 and Hoxc13 expression through TEAD transcription factors. These results provide mechanistic insights into abnormal YAP activities in mice and humans.

Genetic interactions between Shox2 and Hox genes during the regional growth and development of the mouse limb

The growth and development of the vertebrate limb relies on homeobox genes of the Hox and Shox families, with their independent mutation often giving dose-dependent effects. Here we investigate whether Shox2 and Hox genes function together during mouse limb development by modulating their relative dosage and examining the limb for nonadditive effects on growth. Using double mRNA fluorescence in situ hybridization (FISH) in single embryos, we first show that Shox2 and Hox genes have associated spatial expression dynamics, with Shox2 expression restricted to the proximal limb along with Hoxd9 and Hoxa11 expression, juxtaposing the distal expression of Hoxa13 and Hoxd13. By generating mice with all possible dosage combinations of mutant Shox2 alleles and HoxA/D cluster deletions, we then show that their coordinated proximal limb expression is critical to generate normally proportioned limb segments. These epistatic interactions tune limb length, where Shox2 underexpression enhances, and Shox2 overexpression suppresses, Hox-mutant phenotypes. Disruption of either Shox2 or Hox genes leads to a similar reduction in Runx2 expression in the developing humerus, suggesting their concerted action drives cartilage maturation during normal development. While we furthermore provide evidence that Hox gene function influences Shox2 expression, this regulation is limited in extent and is unlikely on its own to be a major explanation for their genetic interaction. Given the similar effect of human SHOX mutations on regional limb growth, Shox and Hox genes may generally function as genetic interaction partners during the growth and development of the proximal vertebrate limb.

Temporal dynamics and developmental memory of 3D chromatin architecture at Hox gene loci

Hox genes are essential regulators of embryonic development. Their step-wise transcriptional activation follows their genomic topology and the various states of activation are subsequently memorized into domains of progressively overlapping gene products. We have analyzed the 3D chromatin organization of Hox clusters during their early activation in vivo, using high-resolution circular chromosome conformation capture. Initially, Hox clusters are organized as single chromatin compartments containing all genes and bivalent chromatin marks. Transcriptional activation is associated with a dynamic bi-modal 3D organization, whereby the genes switch autonomously from an inactive to an active compartment. These local 3D dynamics occur within a framework of constitutive interactions within the surrounding Topological Associated Domains, indicating that this regulation process is mostly cluster intrinsic. The step-wise progression in time is fixed at various body levels and thus can account for the chromatin architectures previously described at a later stage for different anterior to posterior levels.DOI: http://dx.doi.org/10.7554/eLife.02557.001.

The genetics of murine Hox loci: TAMERE, STRING, and PANTHERE to engineer chromosome variants

Following their duplications at the base of the vertebrate clade, Hox gene clusters underwent remarkable sub- and neo-functionalization events. Many of these evolutionary innovations can be associated with changes in the transcriptional regulation of their genes, where an intricate relationship between the structure of the gene cluster and the architecture of the surrounding genomic landscape is at play. Here, we report on a portfolio of in vivo genome engineering strategies in mice, which have been used to probe and decipher the genetic and molecular underpinnings of the complex regulatory mechanisms implemented at these loci.

Clustering of tissue-specific sub-TADs accompanies the regulation of HoxA genes in developing limbs

HoxA genes exhibit central roles during development and causal mutations have been found in several human syndromes including limb malformation. Despite their importance, information on how these genes are regulated is lacking. Here, we report on the first identification of bona fide transcriptional enhancers controlling HoxA genes in developing limbs and show that these enhancers are grouped into distinct topological domains at the sub-megabase scale (sub-TADs). We provide evidence that target genes and regulatory elements physically interact with each other through contacts between sub-TADs rather than by the formation of discreet “DNA loops”. Interestingly, there is no obvious relationship between the functional domains of the enhancers within the limb and how they are partitioned among the topological domains, suggesting that sub-TAD formation does not rely on enhancer activity. Moreover, we show that suppressing the transcriptional activity of enhancers does not abrogate their contacts with HoxA genes. Based on these data, we propose a model whereby chromatin architecture defines the functional landscapes of enhancers. From an evolutionary standpoint, our data points to the convergent evolution of HoxA and HoxD regulation in the fin-to-limb transition, one of the major morphological innovations in vertebrates.

Hox genes encode transcription factors with crucial roles during development. These genes are grouped in four different clusters names HoxA, B, C, and D. Mutations in genes of the HoxA and D clusters have been found in several human syndromes, affecting in some cases limb development. Despite their essential role and contrary to the genes of the HoxD cluster, little is known about how the HoxA genes are regulated. Here, we identified a large set of regulatory elements controlling HoxA genes during limb development. By studying spatial chromatin organization at the HoxA region, we found that the regulatory elements are spatially clustered regardless of their activity. Clustering of enhancers define tissue-specific chromatin domains that interact specifically with each other and with active genes in the limb. Our findings give support to the emerging concept that chromatin architecture defines the functional properties of genomes. Additionally, our study suggests a common constraint of the chromatin topology in the evolution of HoxA and HoxD regulation in the emergence of the hand/foot, which is one of the major morphological innovations in vertebrates.

The genetic landscape of clinical resistance to RAF inhibition in metastatic melanoma

Most patients with BRAF(V600)-mutant metastatic melanoma develop resistance to selective RAF kinase inhibitors. The spectrum of clinical genetic resistance mechanisms to RAF inhibitors and options for salvage therapy are incompletely understood. We performed whole-exome sequencing on formalin-fixed, paraffin-embedded tumors from 45 patients with BRAF(V600)-mutant metastatic melanoma who received vemurafenib or dabrafenib monotherapy. Genetic alterations in known or putative RAF inhibitor resistance genes were observed in 23 of 45 patients (51%). Besides previously characterized alterations, we discovered a "long tail" of new mitogen-activated protein kinase (MAPK) pathway alterations (MAP2K2, MITF) that confer RAF inhibitor resistance. In three cases, multiple resistance gene alterations were observed within the same tumor biopsy. Overall, RAF inhibitor therapy leads to diverse clinical genetic resistance mechanisms, mostly involving MAPK pathway reactivation. Novel therapeutic combinations may be needed to achieve durable clinical control of BRAF(V600)-mutant melanoma. Integrating clinical genomics with preclinical screens may model subsequent resistance studies.

Targeted disruption of Hotair leads to homeotic transformation and gene derepression

Long noncoding RNAs (lncRNAs) are thought to be prevalent regulators of gene expression, but the consequences of lncRNA inactivation in vivo are mostly unknown. Here, we show that targeted deletion of mouse Hotair lncRNA leads to derepression of hundreds of genes, resulting in homeotic transformation of the spine and malformation of metacarpal-carpal bones. RNA sequencing and conditional inactivation reveal an ongoing requirement of Hotair to repress HoxD genes and several imprinted loci such as Dlk1-Meg3 and Igf2-H19 without affecting imprinting choice. Hotair binds to both Polycomb repressive complex 2, which methylates histone H3 at lysine 27 (H3K27), and Lsd1 complex, which demethylates histone H3 at lysine 4 (H3K4) in vivo. Hotair inactivation causes H3K4me3 gain and, to a lesser extent, H3K27me3 loss at target genes. These results reveal the function and mechanisms of Hotair lncRNA in enforcing a silent chromatin state at Hox and additional genes.

Hox genes regulation in vertebrates

Hox genes encode transcription factors defining cellular identities along the major and secondary body axes. Their coordinated expression in both space and time is critical for embryonic patterning. Accordingly, Hox genes transcription is tightly controlled at multiple levels, and involves an intricate combination of local and long-range cis-regulatory elements. Recent studies revealed that in addition to transcription factors, dynamic patterns of histone marks and higher-order chromatin structure are important determinants of Hox gene regulation. Furthermore, the emerging picture suggests an involvement of various species of non-coding RNA in targeting activating and repressive complexes to Hox clusters. I review these recent developments and discuss their relevance to the control of Hox gene expression in vivo, as well as to our understanding of transcriptional regulatory mechanisms.

Transgene- and locus-dependent imprinting reveals allele-specific chromosome conformations

When positioned into the integrin α-6 gene, an Hoxd9lacZ reporter transgene displayed parental imprinting in mouse embryos. While the expression from the paternal allele was comparable with patterns seen for the same transgene when present at the neighboring HoxD locus, almost no signal was scored at this integration site when the transgene was inherited from the mother, although the Itga6 locus itself is not imprinted. The transgene exhibited maternal allele-specific DNA hypermethylation acquired during oogenesis, and its expression silencing was reversible on passage through the male germ line. Histone modifications also corresponded to profiles described at known imprinted loci. Chromosome conformation analyses revealed distinct chromatin microarchitectures, with a more compact structure characterizing the maternally inherited repressed allele. Such genetic analyses of well-characterized transgene insertions associated with a de novo-induced parental imprint may help us understand the molecular determinants of imprinting.

Decoupling the function of Hox and Shh in developing limb reveals multiple inputs of Hox genes on limb growth

Limb development relies on an exquisite coordination between growth and patterning, but the underlying mechanisms remain elusive. Anterior-posterior and proximal-distal specification initiates in early limb bud concomitantly with the proliferative expansion of limb cells. Previous studies have shown that limb bud growth initially relies on fibroblast growth factors (FGFs) produced in the apical ectodermal ridge (AER-FGFs), the maintenance of which relies on a positive-feedback loop involving sonic hedgehog (Shh) and the BMP antagonist gremlin 1 (Grem1). The positive cross-regulation between Shh and the HoxA and HoxD clustered genes identified an indirect effect of Hox genes on the maintenance of AER-FGFs but the respective function of Shh and Hox genes in this process remains unknown. Here, by uncoupling Hox and Shh function, we show that HoxA and HoxD genes are required for proper AER-FGFs expression, independently of their function in controlling Shh expression. In addition, we provide evidence that the Hox-dependent control of AER-FGF expression is achieved through the regulation of key mesenchymal signals, namely Grem1 and Fgf10, ensuring proper epithelial-mesenchymal interactions. Notably, HoxA and HoxD genes contribute to both the initial activation of Grem1 and the subsequent anterior expansion of its expression domain. We propose that the intricate interactions between Hox genes and the FGF and Shh signaling pathways act as a molecular network that ensures proper limb bud growth and patterning, probably contributing to the coordination of these two processes.

Chromatin organization and global regulation of Hox gene clusters

During development, a properly coordinated expression of Hox genes, within their different genomic clusters is critical for patterning the body plans of many animals with a bilateral symmetry. The fascinating correspondence between the topological organization of Hox clusters and their transcriptional activation in space and time has served as a paradigm for understanding the relationships between genome structure and function. Here, we review some recent observations, which revealed highly dynamic changes in the structure of chromatin at Hox clusters, in parallel with their activation during embryonic development. We discuss the relevance of these findings for our understanding of large-scale gene regulation.

Duplications of hox gene clusters and the emergence of vertebrates

The vertebrate body plan is characterized by an increased complexity relative to that of all other chordates and large-scale gene amplifications have been associated with key morphological innovations leading to their remarkable evolutionary success. Here, we use compound full Hox clusters deletions to investigate how Hox genes duplications may have contributed to the emergence of vertebrate-specific innovations. We show that the combined deletion of HoxA and HoxB leads to an atavistic heart phenotype, suggesting that the ancestral HoxA/B cluster was co-opted to help in diversifying the complex organ in vertebrates. Other phenotypic effects observed seem to illustrate the resurgence of ancestral (plesiomorphic) features. This indicates that the duplications of Hox clusters were associated with the recruitment or formation of novel cis-regulatory controls, which were key to the evolution of many vertebrate features and hence to the evolutionary radiation of this group.

Delineating a conserved genetic cassette promoting outgrowth of body appendages

The acquisition of the external genitalia allowed mammals to cope with terrestrial-specific reproductive needs for internal fertilization, and thus it represents one of the most fundamental steps in evolution towards a life on land. How genitalia evolved remains obscure, and the key to understanding this process may lie in the developmental genetics that underpins the early establishment of the genital primordium, the genital tubercle (GT). Development of the GT is similar to that of the limb, which requires precise regulation from a distal signaling epithelium. However, whether outgrowth of the GT and limbs is mediated by common instructive signals remains unknown. In this study, we used comprehensive genetic approaches to interrogate the signaling cascade involved in GT formation in comparison with limb formation. We demonstrate that the FGF ligand responsible for GT development is FGF8 expressed in the cloacal endoderm. We further showed that forced Fgf8 expression can rescue limb and GT reduction in embryos deficient in WNT signaling activity. Our studies show that the regulation of Fgf8 by the canonical WNT signaling pathway is mediated in part by the transcription factor SP8. Sp8 mutants elicit appendage defects mirroring WNT and FGF mutants, and abolishing Sp8 attenuates ectopic appendage development caused by a gain-of-function β-catenin mutation. These observations indicate that a conserved WNT-SP8-FGF8 genetic cassette is employed by both appendages for promoting outgrowth, and suggest a deep homology shared by the limb and external genitalia.

Mammalian limbs and external genitalia are body appendages specialized for locomotion and internal fertilization, respectively. Despite their marked anatomical and functional differences, development of the limb and external genitalia appears to involve similar genetic controls, and some have suggested that regulatory mechanisms common to both might be evolutionarily linked. One essential aspect for appendage development is the establishment and maintenance of a separated proximodistal developmental axis apart from the main body axis, which is often instructed by a distal signaling epithelium. Herein, we adopted comprehensive mouse genetic approaches to investigate regulatory mechanisms underlying the distal signaling center in the limb and the GT, and uncovered a conserved genetic cassette that is utilized by both paired and unpaired appendages to establish a distal signaling center in the epithelium that mediates subsequent proximodistal outgrowth. Our results further suggested that the evolution of the external genital organ involved co-option of the same genetic program underpinning limb development.

Hoxa11 and Hoxd11 regulate chondrocyte differentiation upstream of Runx2 and Shox2 in mice

During limb development, posterior Hox genes of the Hoxa- and Hoxd cluster provide positional information along the limb axis. Here we report a new function for Hoxa11 and Hoxd11 in regulating the early steps of chondrocyte differentiation. We analyzed forelimbs of Hoxa11^−/−;d11^−/− and Ulnaless mice, which are characterized by specifically shortened zeugopods. By detailed morphological and molecular analyses, we show that loss of Hoxa11 and Hoxd11 in the ulna of both mutants leads to an arrest of chondrocyte differentiation at a step before the separation into round and columnar cells takes place. Furthermore, we demonstrate that Hoxa11 and Hoxd11 act upstream of Runx2 and Shox2, two key regulators of chondrocyte differentiation. We hypothesize that Runx2 activates Shox2 in early chondrocytes, which at later stages induces Runx2 expression to regulate hypertrophic differentiation. These results give insight into mechanisms by which positional information might be translated into a specific bone pattern.

Loss of abdominal muscle in Pitx2 mutants associated with altered axial specification of lateral plate mesoderm

Sequence specific transcription factors (SSTFs) combinatorially define cell types during development by forming recursively linked network kernels. Pitx2 expression begins during gastrulation, together with Hox genes, and becomes localized to the abdominal lateral plate mesoderm (LPM) before the onset of myogenesis in somites. The somatopleure of Pitx2 null embryos begins to grow abnormally outward before muscle regulatory factors (MRFs) or Pitx2 begin expression in the dermomyotome/myotome. Abdominal somites become deformed and stunted as they elongate into the mutant body wall, but maintain normal MRF expression domains. Subsequent loss of abdominal muscles is therefore not due to defects in specification, determination, or commitment of the myogenic lineage. Microarray analysis was used to identify SSTF families whose expression levels change in E10.5 interlimb body wall biopsies. All Hox9-11 paralogs had lower RNA levels in mutants, whereas genes expressed selectively in the hypaxial dermomyotome/myotome and sclerotome had higher RNA levels in mutants. In situ hybridization analyses indicate that Hox gene expression was reduced in parts of the LPM and intermediate mesoderm of mutants. Chromatin occupancy studies conducted on E10.5 interlimb body wall biopsies showed that Pitx2 protein occupied chromatin sites containing conserved bicoid core motifs in the vicinity of Hox 9-11 and MRF genes. Taken together, the data indicate that Pitx2 protein in LPM cells acts, presumably in combination with other SSTFs, to repress gene expression, that are normally expressed in physically adjoining cell types. Pitx2 thereby prevents cells in the interlimb LPM from adopting the stable network kernels that define sclerotomal, dermomyotomal, or myotomal mesenchymal cell types. This mechanism may be viewed either as lineage restriction or specification.

Heterochronically early decline of Hox expression prior to cartilage formation in the avian hindlimb zeugopod

The fibula, a zeugopod bone in the hindlimb, exhibits various morphologies in tetrapod species. The fibula in some species has a similar length with the other zeugopod element, the tibia, while other species have obvious differences in the sizes of the two elements. In the avian hindlimb, for example, the fibula is extremely short, thin, and truncated. Basic morphology of the fibula is established during development, and cartilage primordium of the bone emerges in a certain region defined by a distinct combination of expression of Hox genes (Hox code). In order to elucidate how the different morphologies are produced from a region that is defined as the fixed Hox code, we examined spatial and temporal patterns of Hoxd11/Hoxd12 expression in the developing limb bud, which defines the region from which the fibula emerges, in comparison with the sites of precartilaginous mesenchymal condensations representing regions for cartilage formation among chick, mouse, and gecko embryos. We found that in the chick hindlimb, expression of Hoxd11/Hoxd12 decreased and disappeared from the presumptive zeugopod region before cartilage formation. This heterochronically early decline of expression of Hox genes is strongly correlated with the peculiar trait of the fibula in the avian hindlimb, since in the other species examined, expression of those genes continued after the onset of cartilage formation. This is morphological phenotype-related because the early disappearance was not seen in the chick forelimb. Our results suggest that temporal change of the Hox code governs diversification in morphology of homologous structures among related species.

Patterns of cis regulatory variation in diverse human populations

The genetic basis of gene expression variation has long been studied with the aim to understand the landscape of regulatory variants, but also more recently to assist in the interpretation and elucidation of disease signals. To date, many studies have looked in specific tissues and population-based samples, but there has been limited assessment of the degree of inter-population variability in regulatory variation. We analyzed genome-wide gene expression in lymphoblastoid cell lines from a total of 726 individuals from 8 global populations from the HapMap3 project and correlated gene expression levels with HapMap3 SNPs located in cis to the genes. We describe the influence of ancestry on gene expression levels within and between these diverse human populations and uncover a non-negligible impact on global patterns of gene expression. We further dissect the specific functional pathways differentiated between populations. We also identify 5,691 expression quantitative trait loci (eQTLs) after controlling for both non-genetic factors and population admixture and observe that half of the cis-eQTLs are replicated in one or more of the populations. We highlight patterns of eQTL-sharing between populations, which are partially determined by population genetic relatedness, and discover significant sharing of eQTL effects between Asians, European-admixed, and African subpopulations. Specifically, we observe that both the effect size and the direction of effect for eQTLs are highly conserved across populations. We observe an increasing proximity of eQTLs toward the transcription start site as sharing of eQTLs among populations increases, highlighting that variants close to TSS have stronger effects and therefore are more likely to be detected across a wider panel of populations. Together these results offer a unique picture and resource of the degree of differentiation among human populations in functional regulatory variation and provide an estimate for the transferability of complex trait variants across populations.

A function for all posterior Hoxd genes during digit development?

BACKGROUND

Four posterior Hoxd genes, from Hoxd13 to Hoxd10, are collectively regulated during the development of tetrapod digits. Besides the well-documented role of Hoxd13, the function of the neighboring genes has been difficult to evaluate due to the close genetic linkage and potential regulatory interferences. We used a combination of five small nested deletions in cis, involving from two to four consecutive genes of the Hoxd13 to Hoxd9 loci, in mice, to evaluate their combined functional importance.

RESULTS

We show that deletions leading to a gain of function of Hoxd13, via regulatory re-allocation, generate abnormal phenotypes, in agreement with the dominant negative role of this gene. We also show that Hoxd10, Hoxd11, and Hoxd12 all seem to play a genuine role in digit development, though less compelling than that of Hoxd13. In contrast, the nearby Hoxd9 contributed no measurable function in digits.

CONCLUSIONS

We conclude that a slight and transient deregulation of Hoxd13 expression can readily affect the relative lengths of limb segments and that all posterior Hoxd genes likely contribute to the final limb morphology. We discuss the difficulty to clearly assess the functional share of individual genes within such a gene family, where closely located neighbors, coding for homologous proteins, are regulated by a unique circuitry and all contribute to shape the distal parts of our appendages.

A genetic approach to the transcriptional regulation of Hox gene clusters

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

Recruitment of 5' Hoxa genes in the allantois is essential for proper extra-embryonic function in placental mammals

The Hox gene family is well known for its functions in establishing morphological diversity along the anterior-posterior axis of developing embryos. In mammals, one of these genes, Hoxa13, is crucial for embryonic survival, as its function is required for the proper expansion of the fetal vasculature in the placenta. Thus, it appears that the developmental strategy specific to placental mammals is linked, at least in part, to the recruitment of Hoxa13 function in developing extra-embryonic tissues. Yet, the mechanism underlying this extra-embryonic recruitment is unknown. Here, we provide evidence that this functional novelty is not exclusive to Hoxa13 but is shared with its neighboring Hoxa11 and Hoxa10 genes. We show that the extra-embryonic function of these three Hoxa genes stems from their specific expression in the allantois, an extra-embryonic hallmark of amniote vertebrates. Interestingly, Hoxa10-13 expression in the allantois is conserved in chick embryos, which are non-placental amniotes, suggesting that the extra-embryonic recruitment of Hoxa10, Hoxa11 and Hoxa13 most likely arose in amniotes, i.e. prior to the emergence of placental mammals. Finally, using a series of targeted recombination and transgenic assays, we provide evidence that the regulatory mechanism underlying Hoxa expression in the allantois is extremely complex and relies on several cis-regulatory sequences.