Breaking colinearity in the mouse HoxD complex

Teleost Hox code defines regional identities competent for the formation of dorsal and anal fins

The dorsal and anal fins can vary widely in position and length along the anterior-posterior axis in teleost fishes. However, the molecular mechanisms underlying the diversification of these fins remain unknown. Here, we used genetic approaches in zebrafish and medaka, in which the relative positions of the dorsal and anal fins are opposite, to demonstrate the crucial role of hox genes in the patterning of the teleost posterior body, including the dorsal and anal fins. By the CRISPR-Cas9-induced frameshift mutations and positional cloning of spontaneous dorsalfinless medaka, we show that various hox mutants exhibit the absence of dorsal or anal fins, or a stepwise posterior extension of these fins, with vertebral abnormalities. Our results indicate that multiple hox genes, primarily from hoxc-related clusters, encompass the regions responsible for the dorsal and anal fin formation along the anterior-posterior axis. These results further suggest that shifts in the anterior boundaries of hox expression which vary among fish species, lead to diversification in the position and size of the dorsal and anal fins, similar to how modulations in Hox expression can alter the number of anatomically distinct vertebrae in tetrapods. Furthermore, we show that hox genes responsible for dorsal fin formation are different between zebrafish and medaka. Our results suggest that a novel mechanism has occurred during teleost evolution, in which the gene network responsible for fin formation might have switched to the regulation downstream of other hox genes, leading to the remarkable diversity in the dorsal fin position.

Disappearance of Temporal Collinearity in Vertebrates and Its Eventual Reappearance

Simple Summary

In 1999 T. Kondo and D. Duboule performed excisions of posterior upstream DNA domains in mouse embryos and they observed that for an extended excision (including Evx gene) the Hox genes of the cluster were simultaneously expressed with the first Hoxd1 gene ‘as if’ Temporal Collinearity (TC) had disappeared. According to a Biophysical Model (BM) during Hox gene expression, Hox clusters behave similar toexpanding elastic springs. For the extended upstream DNA excision, BM predicts the TC disappearance and an experiment is proposed to test this BM prediction. In the chick limb bud C. Tickle et al. observed that the excision of the apical ectodermal ridge (AER) caused the inhibition of HoxA13 expression. However, the implantation of FGF soaked beads at the tip of the limb could surprisingly rescue HoxA13 expression after 24 h so that TC is restored.Brachyury transcription factor (TF) is essential in identifying the targets of this transcription and a chromatin immunoprecipitation microarray chip (ChIP-chip) was produced which can be inserted in the mouse embryonic cells. It is here proposed to insert this chip in the mutant cells where TC has disappeared and compare it to the limb bud case.Is TC restored? It is an important issue worth exploring.

Abstract

It was observed that a cluster of ordered genes (Hox1, Hox2, Hox3…) in the genome are activated in the ontogenetic units (1, 2, 3 …) of an embryo along the Anterior/Posterior axis following the same order of the Hox genes. This Spatial Collinearity (SC) is very strange since it correlates events of very different spatial dimensions. It was later observed in vertebrates, that, in the above ordering, first is Hox1expressed in ontogenetic unit 1, followed later by Hox2 in unit 2 and even later Hox3 in unit 3. This temporal collinearity (TC) is an enigma and even to-day is explored in depth. In 1999 T. Kondo and D. Duboule, after posterior upstream extended DNA excisions, concluded that the Hox cluster behaves ‘as if’ TC disappears. Here the consideration of TC really disappearing is taken face value and its repercussions are analyzed. Furthermore, an experiment is proposed to test TC disappearance. An outcome of this experiment could be the reappearance (partial or total) of TC.

A novel microdeletion upstream of HOXD13 in a Chinese family with synpolydactyly

Synpolydactyly (SPD) is a digital malformation with the typical clinical phenotype of the webbing of 3/4 fingers and/or 4/5 toes, and combined with polydactyly. In this study, we investigated a Chinese family with SPD and genetic analysis found that all of the affected individuals in the family carry a heterozygous 11,451 bp microdeletion at chr2:176933872-176945322 (GRCh37), which is located upstream of HOXD13 gene, the known disease gene for SPD1. All the affected individuals in the family carry the heterozygous deletion variant, and the variant co-segregated with SPD in the family. Thus, we speculate that the 11,451 bp microdeletion is the disease-causing variant in the family. To date, the microdeletion associating with SPD1 which we identified is the smallest deletion upstream of the HOXD13 gene and not altering the sequence of the HOXD13 gene.

Physical Laws Shape Up HOX Gene Collinearity

Hox gene collinearity (HGC) is a multi-scalar property of many animal phyla particularly important in embryogenesis. It relates entities and events occurring in Hox clusters inside the chromosome DNA and in embryonic tissues. These two entities differ in linear size by more than four orders of magnitude. HGC is observed as spatial collinearity (SC), where the Hox genes are located in the order (Hox1, Hox2, Hox3 …) along the 3′ to 5′ direction of DNA in the genome and a corresponding sequence of ontogenetic units (E1, E2, E3, …) located along the Anterior—Posterior axis of the embryo. Expression of Hox1 occurs in E1, Hox2 in E2, Hox3 in E3, etc. Besides SC, a temporal collinearity (TC) has been also observed in many vertebrates. According to TC, first Hox1 is expressed in E1; later, Hox2 is expressed in E2, followed by Hox3 in E3, etc. Lately, doubt has been raised about whether TC really exists. A biophysical model (BM) was formulated and tested during the last 20 years. According to BM, physical forces are created which pull the Hox genes one after the other, driving them to a transcription factory domain where they are transcribed. The existing experimental data support this BM description. Symmetry is a physical–mathematical property of matter that was explored in depth by Noether who formulated a ground-breaking theory (NT) that applies to all sizes of matter. NT may be applied to biology in order to explain the origin of HGC in animals developing not only along the A/P axis, but also to animals with circular symmetry.

Extraction of peroxisome proliferator-activated receptor α agonist-induced lipid metabolism-related and unrelated genes in rat liver and analysis of their genomic location

Although peroxisome proliferator-activated receptor α (PPARα) agonists are obviously hepatocarcinogenic in rodents, they have been widely used for dyslipidemia and proven to be safe for clinical use without respect to the species difference. It is established that PPARα acts as a part of the transcription factor complex, but its precise mechanism is still unknown. Using the data of Toxicogenomics Database, reliable genes responsive to PPARα agonists, clofibrate, fenofibrate and WY-14,643, in rat liver, were extracted from both in vivo and in vitro data, and sorted by their fold increase. It was found that there were many genes responding to fibrates exclusively in vivo. Most of the in vivo specific genes appear to be unrelated to lipid metabolism and are not upregulated in the kidney. Fifty-seven genes directly related to cell proliferation were extracted from in vivo data, but they were not induced in vitro at all. Analysis of PPAR-responsive elements could not explain the observed difference in induction. To evaluate possible interaction between neighboring genes in gene expression, the correlation of the fold changes of neighboring genes for 22 drugs with various PPARα agonistic potencies were calculated for the genes showing more than 2.5 fold induction by 3 fibrates in vivo, and their genomic location was compared with that of the human orthologue. In the present study, many candidates of genes other than lipid metabolism were selected, and these could be good starting points to elucidate the mechanism of PPARα agonist-induced rodent-specific toxicity.

Hox Gene Collinearity May Be Related to Noether Theory on Symmetry and Its Linked Conserved Quantity

Hox Gene Collinearity (HGC) is a fundamental property that controls the development of many animal species, including vertebrates. In the Hox gene clusters, the genes are located in a sequential order Hox1, Hox2, Hox3, etc., along the 3’ to 5’ direction of the cluster in the chromosome. During Hox cluster activation, the Hox genes are expressed sequentially in the ontogenetic units D1, D2, D3, etc., along the anterior–posterior axis (A-P) of the early embryo. This collinearity, first observed by E.B. Lewis, is surprising because the spatial collinearity of these structures (Hox clusters and embryos) correlates entities that differ by about four orders of magnitude. Biomolecular mechanisms alone cannot explain such correlations. Long-range physical interactions, such as diffusion or electric attractions, should be involved. A biophysical model (BM) was formulated, which, in alignment with the biomolecular processes, successfully describes the existing vertebrate genetic engineering data. One hundred years ago, Emmy Noether made a fundamental discovery in mathematics and physics. She proved, rigorously, that a physical system obeying a symmetry law (e.g., rotations or self-similarity) is followed by a conserved physical quantity. It is argued here that HGC obeys a ‘primitive’ self-similarity symmetry. In this case, the associated primitive conserved quantity is the irreversibly increasing ‘ratchet’-like Hoxgene ordering where some genes may be missing. The genes of a vertebrate Hox clusterare located along a finite straight line. The same order follows the ontogenetic unitsof the vertebrate embryo. Therefore, HGC is a manifestation of a primitive Noether Theory (NT). NT may be applied to other than the vertebrate case, for instance, to animals with a circular topological symmetry. For example, the observed abnormal Hox gene ordering of the echinoderm Hox clusters may be reproduced by a double-strand break of the circular Hox gene ordering and its subsequent incorporation in the flanking chromosome.

Remembering where we are: Positional information in salamander limb regeneration

Fifty years ago, Lewis Wolpert defined an important question in developmental biology: how are cell fates determined by the positions of cells within a system? He proposed that cells retain positional values as if they lie within a coordinate system and that the interpretation of these values produces patterns in development. He referred to this concept as positional information. Though initially controversial, this concept of positional information has proven to be profoundly influential in developmental biology. One area in which the influence of Wolpert's theoretical work can be clearly demonstrated is the study of limb regeneration in salamanders. Here, we review the work in limb regeneration leading up to Wolpert defining the concept of positional information and how his theory has guided regeneration research over the subsequent 50 years.

Role of the Chromosome Architectural Factor SMCHD1 in X-Chromosome Inactivation, Gene Regulation, and Disease in Humans

Structural maintenance of chromosomes flexible hinge domain-containing 1 (SMCHD1) is an architectural factor critical for X-chromosome inactivation (XCI) and the repression of select autosomal gene clusters. In mice, homozygous nonsense mutations in Smchd1 cause female-specific embryonic lethality due to an XCI defect. However, although human mutations in SMCHD1 are associated with congenital arhinia and facioscapulohumeral muscular dystrophy type 2 (FSHD2), the diseases do not show a sex-specific bias, despite the essential nature of XCI in humans. To investigate whether there is a dosage imbalance for the sex chromosomes, we here analyze transcriptomic data from arhinia and FSHD2 patient blood and muscle cells. We find that X-linked dosage compensation is maintained in these patients. In mice, SMCHD1 controls not only protocadherin (Pcdh) gene clusters, but also Hox genes critical for craniofacial development. Ablating Smchd1 results in aberrant expression of these genes, coinciding with altered chromatin states and three-dimensional (3D) topological organization. In a subset of FSHD2 and arhinia patients, we also found dysregulation of clustered PCDH, but not HOX genes. Overall, our study demonstrates preservation of XCI in arhinia and FSHD2, and implicates SMCHD1 in the regulation of the 3D organization of select autosomal gene clusters.

HOXA cluster gene expression during osteoblast differentiation involves epigenetic control

The HOXA gene cluster is generally recognized as a pivotal mediator of positional identity in the skeletal system, expression of different orthologues conferring alternative locational phenotype of the vertebrate bone. Strikingly, however, the molecular mechanisms that regulate orthologue-specific expression of different HOXA cluster members in gestating osteoblasts remain largely obscure, but in analogy to the processes observed in acute lymphatic leukemia it is assumed that alternative methylation of HOXA promoter regions drives position specific expression patterns. In an effort to understand HOXA cluster gene expression in osteogenesis we characterize both expression and the epigenetic landscape of the HOXA gene cluster during in vitro osteoblast formation from mesenchymal precursors. We observe that osteoblast formation per se provokes strong upregulation of HOXA gene cluster expression, in particular of midcluster genes, and paradoxal downregulation of HOXA7 and HOXA10. These differences in expression appear related to promoter methylation. LnRNAs HOTAIR and HOTTIP, known to modulate HOXA expression, are also regulated by their promoter methylation processing, but do not correlate with HOXA cluster expression profile. We thus conclude that HOXA expression is profoundly regulated during osteoblast differentiation through canonical methylation-dependent mechanisms but not through the flanking lnRNAs.

Roles and regulation of histone methylation in animal development

Histone methylation can occur at various sites in histone proteins, primarily on lysine and arginine residues, and it can be governed by multiple positive and negative regulators, even at a single site, to either activate or repress transcription. It is now apparent that histone methylation is critical for almost all stages of development, and its proper regulation is essential for ensuring the coordinated expression of gene networks that govern pluripotency, body patterning and differentiation along appropriate lineages and organogenesis. Notably, developmental histone methylation is highly dynamic. Early embryonic systems display unique histone methylation patterns, prominently including the presence of bivalent (both gene-activating and gene-repressive) marks at lineage-specific genes that resolve to monovalent marks during differentiation, which ensures that appropriate genes are expressed in each tissue type. Studies of the effects of methylation on embryonic stem cell pluripotency and differentiation have helped to elucidate the developmental roles of histone methylation. It has been revealed that methylation and demethylation of both activating and repressive marks are essential for establishing embryonic and extra-embryonic lineages, for ensuring gene dosage compensation via genomic imprinting and for establishing body patterning via HOX gene regulation. Not surprisingly, aberrant methylation during embryogenesis can lead to defects in body patterning and in the development of specific organs. Human genetic disorders arising from mutations in histone methylation regulators have revealed their important roles in the developing skeletal and nervous systems, and they highlight the overlapping and unique roles of different patterns of methylation in ensuring proper development.

Homeotic transformation of tails into limbs in anurans

Anuran tadpoles can regenerate their tails after amputation. However, they occasionally form ectopic limbs instead of the lost tail part after vitamin A treatment. This is regarded as an example of a homeotic transformation. In this phenomenon, the developmental fate of the tail blastema is apparently altered from that of a tail to that of limbs, indicating a realignment of positional information in the blastema. Morphological observations and analyses of the development of skeletal elements during the process suggest that positional information in the blastema is rewritten from tail to trunk specification under the influence of vitamin A, resulting in limb formation. Despite the extensive information gained from morphological observations, a comprehensive understanding of this phenomenon also requires molecular data. We review previous studies related to anuran homeotic transformation. The findings of these studies provide a basis for evaluating major hypotheses and identifying molecular data that should be prioritized in future studies. Finally, we argue that positional information for the tail blastema changes to that for a part of the trunk, leading to homeotic transformations. To suggest this hypothesis, we present published data that favor the rewriting of positional information.

The three-dimensional genome: regulating gene expression during pluripotency and development

The precise expression of genes in time and space during embryogenesis is largely influenced by communication between enhancers and promoters, which is propagated and governed by the physical proximity of these elements in the nucleus. Here, we review how chromatin domains organize the genome by guiding enhancers to their target genes thereby preventing non-specific interactions with other neighboring regions. We also discuss the dynamics of chromatin interactions between enhancers and promoters, as well as the consequent changes in gene expression, that occur in pluripotent cells and during development. Finally, we evaluate how genomic changes such as deletions, inversions and duplications affect 3D chromatin configuration overall and lead to ectopic enhancer-promoter contacts, and thus gene misexpression, which can contribute to abnormal development and disease.

Physical Forces May Cause the HoxD Gene Cluster Elongation

Hox gene collinearity was discovered be Edward B. Lewis in 1978. It consists of the Hox1, Hox2, Hox3 ordering of the Hox genes in the chromosome from the telomeric to the centromeric side of the chromosome. Surprisingly, the spatial activation of the Hox genes in the ontogenetic units of the embryo follows the same ordering along the anterior-posterior embryonic axis. The chromosome microscale differs from the embryo macroscale by 3 to 4 orders of magnitude. The traditional biomolecular mechanisms are not adequate to comprise phenomena at so divergent spatial domains. A Biophysical Model of physical forces was proposed which can bridge the intermediate space and explain the results of genetic engineering experiments. Recent progress in constructing instruments and achieving high resolution imaging (e.g., 3D DNA FISH, STORM etc.) enable the assessment of the geometric structure of the chromatin during the different phases of Hox gene activation. It is found that the mouse HoxD gene cluster is elongated up to 5-6 times during Hox gene transcription. These unexpected findings agree with the BM predictions. It is now possible to measure several physical quantities inside the nucleus during Hox gene activation. New experiments are proposed to test further this model.

MicroRNA filters Hox temporal transcription noise to confer boundary formation in the spinal cord

The initial rostrocaudal patterning of the neural tube leads to differential expression of Hox genes that contribute to the specification of motor neuron (MN) subtype identity. Although several 3' Hox mRNAs are expressed in progenitors in a noisy manner, these Hox proteins are not expressed in the progenitors and only become detectable in postmitotic MNs. MicroRNA biogenesis impairment leads to precocious expression and propagates the noise of Hoxa5 at the protein level, resulting in an imprecise Hoxa5-Hoxc8 boundary. Here we uncover, using in silico simulation, two feed-forward Hox-miRNA loops accounting for the precocious and noisy Hoxa5 expression, as well as an ill-defined boundary phenotype in Dicer mutants. Finally, we identify mir-27 as a major regulator coordinating the temporal delay and spatial boundary of Hox protein expression. Our results provide a novel trans Hox-miRNA circuit filtering transcription noise and controlling the timing of protein expression to confer robust individual MN identity.

Exploring CTCF and cohesin related chromatin architecture at HOXA gene cluster in primary human fibroblasts

Spatial expression patterns of homeobox (HOX) genes delineate positional identity of primary fibroblasts from different topographic sites. The molecular mechanism underlying the establishing or maintaining of HOX gene expression pattern remains an attractive developmental issue to be addressed. Our previous work suggested a critical role of CTCF/cohesin-mediated higher- order chromatin structure in RA-induced HOXA activation in human teratocarcinoma NT2/D1 cells. This study investigated the recruitment of CTCF and cohesin, and the higher-order chromatin structure of the HOXA locus in fetal lung and adult foreskin fibroblasts, which display complementary HOXA gene expression patterns. Chromatin contacts between the CTCF-binding sites were observed with lower frequency in human foreskin fibroblasts. This observation is consistent with the lower level of cohesin recruitment and 5' HOXA gene expression in the same cells. We also showed that CTCF-binding site A56 (CBSA56) related chromatin structures exhibit the most notable changes in between the two types of cell, and hence may stand for one of the key CTCF-binding sites for cell-type specific chromatin structure organization. Together, these results imply that CTCF/cohesin coordinates HOXA cluster higher-order chromatin structure and expression during development, and provide insight into the relationship between cell-type specific chromatin organization and the spatial collinearity.

A discrete transition zone organizes the topological and regulatory autonomy of the adjacent tfap2c and bmp7 genes

Despite the well-documented role of remote enhancers in controlling developmental gene expression, the mechanisms that allocate enhancers to genes are poorly characterized. Here, we investigate the cis-regulatory organization of the locus containing the Tfap2c and Bmp7 genes in vivo, using a series of engineered chromosomal rearrangements. While these genes lie adjacent to one another, we demonstrate that they are independently regulated by distinct sets of enhancers, which in turn define non-overlapping regulatory domains. Chromosome conformation capture experiments reveal a corresponding partition of the locus in two distinct structural entities, demarcated by a discrete transition zone. The impact of engineered chromosomal rearrangements on the topology of the locus and the resultant gene expression changes indicate that this transition zone functionally organizes the structural partition of the locus, thereby defining enhancer-target gene allocation. This partition is, however, not absolute: we show that it allows competing interactions across it that may be non-productive for the competing gene, but modulate expression of the competed one. Altogether, these data highlight the prime role of the topological organization of the genome in long-distance regulation of gene expression.

The specificity of enhancer-gene interactions is fundamental to the execution of gene regulatory programs underpinning embryonic development and cell differentiation. However, our understanding of the mechanisms conferring specificity to enhancers and target gene interactions is limited. In this study, we characterize the cis-regulatory organization of a large genomic locus consisting of two developmental genes, Tfap2c and Bmp7. We show that this locus is structurally partitioned into two distinct domains by the constitutive action of a discrete transition zone located between the two genes. This separation restricts selectively the functional action of enhancers to the genes present within the same domain. Interestingly, the effects of this region as a boundary are relative, as it allows some competing interactions to take place across domains. We show that these interactions modulate the functional output of a brain enhancer on its primary target gene resulting in the spatial restriction of its expression domain. These results support a functional link between topological chromatin domains and allocation of enhancers to genes. They further show that a precise adjustment of chromatin interaction levels fine-tunes gene regulation by long-range enhancers.

CTCF controls HOXA cluster silencing and mediates PRC2-repressive higher-order chromatin structure in NT2/D1 cells

HOX cluster genes are activated sequentially in their positional order along the chromosome during vertebrate development. This phenomenon, known as temporal colinearity, depends on transcriptional silencing of 5' HOX genes. Chromatin looping was recently identified as a conserved feature of silent HOX clusters, with CCCTC-binding factor (CTCF) binding sites located at the loop bases. However, the potential contribution of CTCF to HOX cluster silencing and the underlying mechanism have not been established. Here, we demonstrate that the HOXA locus is organized by CTCF into chromatin loops and that CTCF depletion causes significantly enhanced activation of HOXA3 to -A7, -A9 to -A11, and -A13 in response to retinoic acid, with the highest effect observed for HOXA9. Our subsequent analyses revealed that CTCF facilitates the stabilization of Polycomb repressive complex 2 (PRC2) and trimethylated lysine 27 of histone H3 (H3K27me3) at the human HOXA locus. Our results reveal that CTCF functions as a controller of HOXA cluster silencing and mediates PRC2-repressive higher-order chromatin structure.

Variant PRC1 complex-dependent H2A ubiquitylation drives PRC2 recruitment and polycomb domain formation

Chromatin modifying activities inherent to polycomb repressive complexes PRC1 and PRC2 play an essential role in gene regulation, cellular differentiation, and development. However, the mechanisms by which these complexes recognize their target sites and function together to form repressive chromatin domains remain poorly understood. Recruitment of PRC1 to target sites has been proposed to occur through a hierarchical process, dependent on prior nucleation of PRC2 and placement of H3K27me3. Here, using a de novo targeting assay in mouse embryonic stem cells we unexpectedly discover that PRC1-dependent H2AK119ub1 leads to recruitment of PRC2 and H3K27me3 to effectively initiate a polycomb domain. This activity is restricted to variant PRC1 complexes, and genetic ablation experiments reveal that targeting of the variant PCGF1/PRC1 complex by KDM2B to CpG islands is required for normal polycomb domain formation and mouse development. These observations provide a surprising PRC1-dependent logic for PRC2 occupancy at target sites in vivo.

•

Variant PRC1 complex-dependent H2AK119ub1 leads to binding of PRC2 and H3K27me3

•

Canonical PRC1 complexes fail to efficiently deposit H2AK119ub1 and recruit PRC2

•

A variant KDM2B/PCGF1/PRC1 complex is required for polycomb domain formation at CGIs

•

Failure to target KDM2B/PCGF1/PRC1 causes polycomb phenotypes and lethality in mice

Homeotic gene regulation: a paradigm for epigenetic mechanisms underlying organismal development

The organization of eukaryotic genome into chromatin within the nucleus eventually dictates the cell type specific expression pattern of genes. This higher order of chromatin organization is established during development and dynamically maintained throughout the life span. Developmental mechanisms are conserved in bilaterians and hence they have body plan in common, which is achieved by regulatory networks controlling cell type specific gene expression. Homeotic genes are conserved in metazoans and are crucial for animal development as they specify cell type identity along the anterior-posterior body axis. Hox genes are the best studied in the context of epigenetic regulation that has led to significant understanding of the organismal development. Epigenome specific regulation is brought about by conserved chromatin modulating factors like PcG/trxG proteins during development and differentiation. Here we discuss the conserved epigenetic mechanisms relevant to homeotic gene regulation in metazoans.

The regulation of Hox gene expression during animal development

Hox genes encode a family of transcriptional regulators that elicit distinct developmental programmes along the head-to-tail axis of animals. The specific regional functions of individual Hox genes largely reflect their restricted expression patterns, the disruption of which can lead to developmental defects and disease. Here, we examine the spectrum of molecular mechanisms controlling Hox gene expression in model vertebrates and invertebrates and find that a diverse range of mechanisms, including nuclear dynamics, RNA processing, microRNA and translational regulation, all concur to control Hox gene outputs. We propose that this complex multi-tiered regulation might contribute to the robustness of Hox expression during development.

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.

Impact of copy number variations (CNVs) on long-range gene regulation at the HoxD locus

Copy number variations are genomic structural variants that are frequently associated with human diseases. Among these copy number variations, duplications of DNA segments are often assumed to lead to dosage effects by increasing the copy number of either genes or their regulatory elements. We produced a series of large targeted duplications within a conserved gene desert upstream of the murine HoxD locus. This DNA region, syntenic to human 2q31-32, contains a range of regulatory elements required for Hoxd gene transcription, and it is often disrupted and/or reorganized in human genetic conditions collectively known as the 2q31 syndrome. Unexpectedly, one such duplication led to a transcriptional down-regulation in developing digits by impairing physical interactions between the target genes and their upstream regulatory elements, thus phenocopying the effect obtained when these enhancer sequences are deleted. These results illustrate the detrimental consequences of interrupting highly conserved regulatory landscapes and reveal a mechanism where genomic duplications lead to partial loss of function of nearby located genes.

Comparison of models for the collinearity of hox genes in the developmental axes of vertebrates

Hox gene clusters are very frequent in many animal genomes and their role in development is pivotal. Particularly in vertebrates, intensive efforts have established several properties of Hox clusters. The collinearity of Hox gene expressions (spatial, temporal and quantitative) is a common feature of the vertebrates. During the last decade, genetic engineering experiments have revealed some important facets of collinearity during limb and trunk development in mice. Two models have been proposed to explain all these properties. On one hand the ‘two-phases model’ makes use of the molecular regulatory mechanisms acting on the Hox genes. On the other hand, the’biophysical model’ is based on the signals transduced inside the cell nucleus and the generation of forces which apply on the cluster and lead to a coordinated activation of Hox genes. The two models differ fundamentally and a critical and detailed comparison is presented. Furthermore, experiments are proposed for which the two models provide divergent predictions. The outcome of these experiments will help to decide which of the two models is valid (if any).

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.

A regulatory archipelago controls Hox genes transcription in digits

The evolution of digits was an essential step in the success of tetrapods. Among the key players, Hoxd genes are coordinately regulated in developing digits, where they help organize growth and patterns. We identified the distal regulatory sites associated with these genes by probing the three-dimensional architecture of this regulatory unit in developing limbs. This approach, combined with in vivo deletions of distinct regulatory regions, revealed that the active part of the gene cluster contacts several enhancer-like sequences. These elements are dispersed throughout the nearby gene desert, and each contributes either quantitatively or qualitatively to Hox gene transcription in presumptive digits. We propose that this genetic system, which we call a "regulatory archipelago," provides an inherent flexibility that may partly underlie the diversity in number and morphology of digits across tetrapods, as well as their resilience to drastic variations.

Evolution of networks for body plan patterning; interplay of modularity, robustness and evolvability

A major goal of evolutionary developmental biology (evo-devo) is to understand how multicellular body plans of increasing complexity have evolved, and how the corresponding developmental programs are genetically encoded. It has been repeatedly argued that key to the evolution of increased body plan complexity is the modularity of the underlying developmental gene regulatory networks (GRNs). This modularity is considered essential for network robustness and evolvability. In our opinion, these ideas, appealing as they may sound, have not been sufficiently tested. Here we use computer simulations to study the evolution of GRNs' underlying body plan patterning. We select for body plan segmentation and differentiation, as these are considered to be major innovations in metazoan evolution. To allow modular networks to evolve, we independently select for segmentation and differentiation. We study both the occurrence and relation of robustness, evolvability and modularity of evolved networks. Interestingly, we observed two distinct evolutionary strategies to evolve a segmented, differentiated body plan. In the first strategy, first segments and then differentiation domains evolve (SF strategy). In the second scenario segments and domains evolve simultaneously (SS strategy). We demonstrate that under indirect selection for robustness the SF strategy becomes dominant. In addition, as a byproduct of this larger robustness, the SF strategy is also more evolvable. Finally, using a combined functional and architectural approach, we determine network modularity. We find that while SS networks generate segments and domains in an integrated manner, SF networks use largely independent modules to produce segments and domains. Surprisingly, we find that widely used, purely architectural methods for determining network modularity completely fail to establish this higher modularity of SF networks. Finally, we observe that, as a free side effect of evolving segmentation and differentiation in combination, we obtained in-silico developmental mechanisms resembling mechanisms used in vertebrate development.

Concerted involvement of Cdx/Hox genes and Wnt signaling in morphogenesis of the caudal neural tube and cloacal derivatives from the posterior growth zone

Decrease in Cdx dosage in an allelic series of mouse Cdx mutants leads to progressively more severe posterior vertebral defects. These defects are corrected by posterior gain of function of the Wnt effector Lef1. Precocious expression of Hox paralogous 13 genes also induces vertebral axis truncation by antagonizing Cdx function. We report here that the phenotypic similarity also applies to patterning of the caudal neural tube and uro-rectal tracts in Cdx and Wnt3a mutants, and in embryos precociously expressing Hox13 genes. Cdx2 inactivation after placentation leads to posterior defects, including incomplete uro-rectal septation. Compound mutants carrying one active Cdx2 allele in the Cdx4-null background (Cdx2/4), transgenic embryos precociously expressing Hox13 genes and a novel Wnt3a hypomorph mutant all manifest a comparable phenotype with similar uro-rectal defects. Phenotype and transcriptome analysis in early Cdx mutants, genetic rescue experiments and gene expression studies lead us to propose that Cdx transcription factors act via Wnt signaling during the laying down of uro-rectal mesoderm, and that they are operative in an early phase of these events, at the site of tissue progenitors in the posterior growth zone of the embryo. Cdx and Wnt mutations and premature Hox13 expression also cause similar neural dysmorphology, including ectopic neural structures that sometimes lead to neural tube splitting at caudal axial levels. These findings involve the Cdx genes, canonical Wnt signaling and the temporal control of posterior Hox gene expression in posterior morphogenesis in the different embryonic germ layers. They shed a new light on the etiology of the caudal dysplasia or caudal regression range of human congenital defects.

A regulatory 'landscape effect' over the HoxD cluster

Faithful expression of Hox genes in both time and space is essential for proper patterning of the primary body axis. Transgenic approaches in vertebrates have suggested that this collinear activation process is regulated in a largely gene cluster-autonomous manner. In contrast, more recently co-opted expression specificities, required in other embryonic structures, depend upon long-range enhancer sequences acting from outside the gene clusters. This regulatory dichotomy was recently questioned, since gene activation along the trunk seems to be partially regulated by signals located outside of the cluster. We investigated these alternative regulatory strategies by engineering a large inversion that precisely separates the murine HoxD complex from its centromeric neighborhood. Mutant animals displayed posterior transformations along with subtle deregulations of Hoxd genes, indicating an impact of the centromeric landscape on the fine-tuning of Hoxd gene expression. Proximal limbs were also affected, suggesting that this 'landscape effect' is generic and impacts upon regulatory mechanisms of various qualities and evolutionary origins.

Multiple promoters and alternative splicing: Hoxa5 transcriptional complexity in the mouse embryo

Background

The genomic organization of Hox clusters is fundamental for the precise spatio-temporal regulation and the function of each Hox gene, and hence for correct embryo patterning. Multiple overlapping transcriptional units exist at the Hoxa5 locus reflecting the complexity of Hox clustering: a major form of 1.8 kb corresponding to the two characterized exons of the gene and polyadenylated RNA species of 5.0, 9.5 and 11.0 kb. This transcriptional intricacy raises the question of the involvement of the larger transcripts in Hox function and regulation.

Methodology/Principal Findings

We have undertaken the molecular characterization of the Hoxa5 larger transcripts. They initiate from two highly conserved distal promoters, one corresponding to the putative Hoxa6 promoter, and a second located nearby Hoxa7. Alternative splicing is also involved in the generation of the different transcripts. No functional polyadenylation sequence was found at the Hoxa6 locus and all larger transcripts use the polyadenylation site of the Hoxa5 gene. Some larger transcripts are potential Hoxa6/Hoxa5 bicistronic units. However, even though all transcripts could produce the genuine 270 a.a. HOXA5 protein, only the 1.8 kb form is translated into the protein, indicative of its essential role in Hoxa5 gene function. The Hoxa6 mutation disrupts the larger transcripts without major phenotypic impact on axial specification in their expression domain. However, Hoxa5-like skeletal anomalies are observed in Hoxa6 mutants and these defects can be explained by the loss of expression of the 1.8 kb transcript. Our data raise the possibility that the larger transcripts may be involved in Hoxa5 gene regulation.

Significance

Our observation that the Hoxa5 larger transcripts possess a developmentally-regulated expression combined to the increasing sum of data on the role of long noncoding RNAs in transcriptional regulation suggest that the Hoxa5 larger transcripts may participate in the control of Hox gene expression.

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

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

Uncoupling time and space in the collinear regulation of Hox genes

During development of the vertebrate body axis, Hox genes are transcribed sequentially, in both time and space, following their relative positions within their genomic clusters. Analyses of animal genomes support the idea that Hox gene clustering is essential for coordinating the various times of gene activations. However, the eventual collinear ordering of the gene specific transcript domains in space does not always require genomic clustering. We analyzed these complex regulatory relationships by using mutant alleles at the mouse HoxD locus, including one that splits the cluster into two pieces. We show that both positive and negative regulatory influences, located on either side of the cluster, control an early phase of collinear expression in the trunk. Interestingly, this early phase does not systematically impact upon the subsequent expression patterns along the main body axis, indicating that the mechanism underlying temporal collinearity is distinct from those acting during the second phase. We discuss the potential functions and evolutionary origins of these mechanisms, as well as their relationship with similar processes at work during limb development.

Hox genes encode proteins that control embryonic development along the head-to-tail axis. These genes are clustered in one site on the chromosome and their respective positions within the cluster determine their time and place of activation. Here, by using a large set of targeted mutations disturbing the integrity of the gene cluster, we show that the spatial organization of expression domains does not directly depend upon the timing of activation as was previously suggested. This uncoupling between space and time in the regulation of these Hox genes coincides with the existence of two major phases of regulation. The first is time-dependent and involves global regulatory influences, located outside the gene cluster, whereas the second relies upon more local regulatory elements, likely interspersed between the genes, inside the cluster. These results provide the bases for future analyses of collinear mechanisms and indicate that different types of collinearities are not necessarily related, neither in function, nor in their evolutionary histories.

Brpf1, a subunit of the MOZ histone acetyl transferase complex, maintains expression of anterior and posterior Hox genes for proper patterning of craniofacial and caudal skeletons

The epigenetic mechanism involving chromatin modification plays a critical role in the maintenance of the expression of Hox genes. Here, we characterize a mutant of the medaka fish, named biaxial symmetries (bis), in which brpf1, a subunit of the MOZ histone acetyl transferase (HAT) complex, is mutated. The bis mutant displayed patterning defects both in the anterior-posterior axis of the craniofacial skeleton and the dorsal-ventral axis of the caudal one. In the anterior region, the bis mutant exhibited craniofacial cartilage homeosis. The expression of Hox genes was decreased in the pharyngeal arches, suggesting that the pharyngeal segmental identities were altered in the bis mutant. In the posterior region, the bis mutant exhibited abnormal patterning of the caudal skeleton, which ectopically formed at the dorsal side of the caudal fin. The expression of Zic genes was decreased at the posterior region, suggesting that the dorsal-ventral axis formation of the posterior trunk was disrupted in the bis mutant. We also found that the MOZ-deficient mice exhibited an abnormal patterning of their craniofacial and cervical skeletons and a decrease of Hox transcripts. We propose a common role of the MOZ HAT complex in vertebrates, a complex which is required for the proper patterning for skeletal development.

Establishment of Hox vertebral identities in the embryonic spine precursors

The vertebrate spine exhibits two striking characteristics. The first one is the periodic arrangement of its elements-the vertebrae-along the anteroposterior axis. This segmented organization is the result of somitogenesis, which takes place during organogenesis. The segmentation machinery involves a molecular oscillator-the segmentation clock-which delivers a periodic signal controlling somite production. During embryonic axis elongation, this signal is displaced posteriorly by a system of traveling signaling gradients-the wavefront-which depends on the Wnt, FGF, and retinoic acid pathways. The other characteristic feature of the spine is the subdivision of groups of vertebrae into anatomical domains, such as the cervical, thoracic, lumbar, sacral, and caudal regions. This axial regionalization is controlled by a set of transcription factors called Hox genes. Hox genes exhibit nested expression domains in the somites which reflect their linear arrangement along the chromosomes-a property termed colinearity. The colinear disposition of Hox genes expression domains provides a blueprint for the regionalization of the future vertebral territories of the spine. In amniotes, Hox genes are activated in the somite precursors of the epiblast in a temporal colinear sequence and they were proposed to control their progressive ingression into the nascent paraxial mesoderm. Consequently, the positioning of the expression domains of Hox genes along the anteroposterior axis is largely controlled by the timing of Hox activation during gastrulation. Positioning of the somitic Hox domains is subsequently refined through a crosstalk with the segmentation machinery in the presomitic mesoderm. In this review, we focus on our current understanding of the embryonic mechanisms that establish vertebral identities during vertebrate development.

Distinct roles and regulations for HoxD genes in metanephric kidney development

Hox genes encode homeodomain-containing proteins that control embryonic development in multiple contexts. Up to 30 Hox genes, distributed among all four clusters, are expressed during mammalian kidney morphogenesis, but functional redundancy between them has made a detailed functional account difficult to achieve. We have investigated the role of the HoxD cluster through comparative molecular embryological analysis of a set of mouse strains carrying targeted genomic rearrangements such as deletions, duplications, and inversions. This analysis allowed us to uncover and genetically dissect the complex role of the HoxD cluster. Regulation of metanephric mesenchyme-ureteric bud interactions and maintenance of structural integrity of tubular epithelia are differentially controlled by some Hoxd genes during renal development, consistent with their specific expression profiles. We also provide evidence for a kidney-specific form of colinearity that underlies the differential expression of two distinct sets of genes located on both sides and overlapping at the Hoxd9 locus. These insights further our knowledge of the genetic control of kidney morphogenesis and may contribute to understanding certain congenital kidney malformations, including polycystic kidney disease and renal hypoplasia.

Hox genes encode proteins that control embryonic development along the head-to-tail axis and in multiple organs. Here, we show that several members of this gene family are necessary for the normal development of the mammalian kidneys. These genes are clustered in one site on the chromosome and their respective positions within the group determine which component of the kidneys they will contribute to. Using a large collection of engineered mutations in this system, we show that these genes are required both for the growth of the kidneys and for their proper organization, such that mutations in some genes reduce the size of the organs, whereas mutations in others induce polycystic kidneys. Our set of genetic rearrangements also allowed us to localize the position of regulatory sequences, which control the expression of these genes during kidney development.