Bimodal control of Hoxd gene transcription in the spinal cord defines two regulatory subclusters

Hox genes and patterning the vertebrate body

The diversity of vertebrate body plans is dizzying, yet stunning for the many things they have in common. Vertebrates have inhabited virtually every part of the earth from its coldest to warmest climates. They locomote by swimming, flying, walking, slithering, or climbing, or combinations of these behaviors. And they exist in many different sizes, from the smallest of frogs, fish and lizards to giraffes, elephants, and blue whales. Despite these differences, vertebrates follow a remarkably similar blueprint for the establishment of their body plan. Within the relatively small amount of time required to complete gastrulation, the process through which the three germ layers, ectoderm, mesoderm, and endoderm are created, the embryo also generates its body axis and is simultaneously patterned. For the length of this axis, the genes that distinguish the neck from the rib cage or the trunk from the sacrum are the Hox genes. In vertebrates, there was evolutionary pressure to maintain this set of genes in the organism. Over the past decades, much has been learned regarding the regulatory mechanisms that ensure the appropriate expression of these genes along the main body axes. Genetic functions continue to be explored though much has been learned. Much less has been discerned on the identity of co-factors used by Hox proteins for the specificity of transcriptional regulation or what downstream targets and pathways are critical for patterning events, though there are notable exceptions. Current work in the field is demonstrating that Hox genes continue to function in many organs long after directing early patterning events. It is hopeful continued research will shed light on remaining questions regarding mechanisms used by this important and conserved set of transcriptional regulators.

Shared retinoic acid responsive enhancers coordinately regulate nascent transcription of Hoxb coding and non-coding RNAs in the developing mouse neural tube

Signaling pathways regulate the patterns of Hox gene expression that underlie their functions in the specification of axial identity. Little is known about the properties of cis-regulatory elements and underlying transcriptional mechanisms that integrate graded signaling inputs to coordinately control Hox expression. Here, we optimized a single molecule fluorescent in situ hybridization (smFISH) technique with probes spanning introns to evaluate how three shared retinoic acid response element (RARE)-dependent enhancers in the Hoxb cluster regulate patterns of nascent transcription in vivo at the level of single cells in wild-type and mutant embryos. We predominately detect nascent transcription of only a single Hoxb gene in each cell, with no evidence for simultaneous co-transcriptional coupling of all or specific subsets of genes. Single and/or compound RARE mutations indicate that each enhancer differentially impacts global and local patterns of nascent transcription, suggesting that selectivity and competitive interactions between these enhancers is important to robustly maintain the proper levels and patterns of nascent Hoxb transcription. This implies that rapid and dynamic regulatory interactions potentiate transcription of genes through combined inputs from these enhancers in coordinating the retinoic acid response.

Sequential in cis mutagenesis in vivo reveals various functions for CTCF sites at the mouse HoxD cluster

Mammalian Hox gene clusters contain a range of CTCF binding sites. In addition to their importance in organizing a TAD border, which isolates the most posterior genes from the rest of the cluster, the positions and orientations of these sites suggest that CTCF may be instrumental in the selection of various subsets of contiguous genes, which are targets of distinct remote enhancers located in the flanking regulatory landscapes. We examined this possibility by producing an allelic series of cumulative in cis mutations in these sites, up to the abrogation of CTCF binding in the five sites located on one side of the TAD border. In the most impactful alleles, the global chromatin architecture of the locus was modified, yet not drastically, illustrating that CTCF sites located on one side of a strong TAD border are sufficient to organize at least part of this insulation. Spatial colinearity in the expression of these genes along the major body axis was nevertheless maintained, despite abnormal expression boundaries. In contrast, strong effects were scored in the selection of target genes responding to particular enhancers, leading to the misregulation of Hoxd genes in specific structures. Altogether, while most enhancer-promoter interactions can occur in the absence of this series of CTCF sites, the binding of CTCF in the Hox cluster is required to properly transform a rather unprecise process into a highly discriminative mechanism of interactions, which is translated into various patterns of transcription accompanied by the distinctive chromatin topology found at this locus. Our allelic series also allowed us to reveal the distinct functional contributions for CTCF sites within this Hox cluster, some acting as insulator elements, others being necessary to anchor or stabilize enhancer-promoter interactions, and some doing both, whereas they all together contribute to the formation of a TAD border. This variety of tasks may explain the amazing evolutionary conservation in the distribution of these sites among paralogous Hox clusters or between various vertebrates.

HoxD transcription factors define monosynaptic sensory-motor specificity in the developing spinal cord

The specificity of monosynaptic connections between proprioceptive sensory neurons and their recipient spinal motor neurons depends on multiple factors, including motor neuron positioning and dendrite morphology, axon projection patterns of proprioceptive sensory neurons in the spinal cord, and the ligand-receptor molecules involved in cell-to-cell recognition. However, with few exceptions, the transcription factors engaged in this process are poorly characterized. Here, we show that members of the HoxD family of transcription factors play a crucial role in the specificity of monosynaptic sensory-motor connections. Mice lacking Hoxd9, Hoxd10 and Hoxd11 exhibit defects in locomotion but have no obvious defects in motor neuron positioning or dendrite morphology through the medio-lateral and rostro-caudal axes. However, we found that quadriceps motor neurons in these mice show aberrant axon development and receive inappropriate inputs from proprioceptive sensory axons innervating the obturator muscle. These genetic studies demonstrate that the HoxD transcription factors play an integral role in the synaptic specificity of monosynaptic sensory-motor connections in the developing spinal cord.

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.

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

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

Deletion of 2.7 kb near HOXD3 in an Arabian horse with occipitoatlantoaxial malformation

In the horse, the term occipitoatlantoaxial malformation (OAAM) is used to describe a developmental defect in which the first cervical vertebra (atlas) resembles the base of the skull (occiput) and the second cervical vertebra (axis) resembles the atlas. Affected individuals demonstrate an abnormal posture and varying degrees of ataxia. The homeobox (HOX) gene cluster is involved in the development of both the axial and appendicular skeleton. Hoxd3-null mice demonstrate a strikingly similar phenotype to Arabian foals with OAAM. Whole-genome sequencing was performed in an OAAM-affected horse (OAAM1) and seven unaffected Arabian horses. Visual inspection of the raw reads within the region of HOXD3 identified a 2.7-kb deletion located 4.4 kb downstream of the end of HOXD4 and 8.2 kb upstream of the start of HOXD3. A genotyping assay revealed that both parents of OAAM1 were heterozygous for the deletion. Additional genotyping identified two of 162 heterozygote Arabians, and the deletion was not present in 371 horses of other breeds. Comparative genomics studies have revealed that this region is highly conserved across species and that the entire genomic region between Hoxd4 and Hoxd3 is transcribed in mice. Two additional Arabian foals diagnosed with OAAM (OAAM 2 and 3) were genotyped and did not have the 2.7-kb deletion. Closer examination of the phenotype in these cases revealed notable variation. OAAM3 also had facial malformations and a patent ductus arteriosus, and the actual malformation at the craniocervical junction differed. Genetic heterogeneity may exist across the HOXD locus in Arabian foals with OAAM.

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.

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.

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.

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.

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.

Genetic and functional modularity of Hox activities in the specification of limb-innervating motor neurons

A critical step in the assembly of the neural circuits that control tetrapod locomotion is the specification of the lateral motor column (LMC), a diverse motor neuron population targeting limb musculature. Hox6 paralog group genes have been implicated as key determinants of LMC fate at forelimb levels of the spinal cord, through their ability to promote expression of the LMC-restricted genes Foxp1 and Raldh2 and to suppress thoracic fates through exclusion of Hoxc9. The specific roles and mechanisms of Hox6 gene function in LMC neurons, however, are not known. We show that Hox6 genes are critical for diverse facets of LMC identity and define motifs required for their in vivo specificities. Although Hox6 genes are necessary for generating the appropriate number of LMC neurons, they are not absolutely required for the induction of forelimb LMC molecular determinants. In the absence of Hox6 activity, LMC identity appears to be preserved through a diverse array of Hox5–Hox8 paralogs, which are sufficient to reprogram thoracic motor neurons to an LMC fate. In contrast to the apparently permissive Hox inputs to early LMC gene programs, individual Hox genes, such as Hoxc6, have specific roles in promoting motor neuron pool diversity within the LMC. Dissection of motifs required for Hox in vivo specificities reveals that either cross-repressive interactions or cooperativity with Pbx cofactors are sufficient to induce LMC identity, with the N-terminus capable of promoting columnar, but not pool, identity when transferred to a heterologous homeodomain. These results indicate that Hox proteins orchestrate diverse aspects of cell fate specification through both the convergent regulation of gene programs regulated by many paralogs and also more restricted actions encoded through specificity determinants in the N-terminus.

Coordinated motor behaviors—as complex as playing a musical instrument or as simple as walking—rely on the ability of motor neurons within the spinal cord to navigate towards and establish specific connections with muscles in the limbs. The establishment of connections between motor neurons and limb muscles is mediated through the actions of genes encoding Hox proteins, a large family of transcription factors conserved amongst all metazoans. However, the specific requirements for Hox genes in motor neuron specification and patterns of muscle connectivity are poorly understood. We have found that members of the Hox6 gene paralog group (Hoxa6, Hoxc6, and Hoxb6) contribute to diverse aspects of motor neuron subtype differentiation. Hox6 gene activity is required during two critical phases of motor neuron development: first as motor axons select a trajectory toward the forelimb and second as they choose specific muscles to innervate. At the molecular level, these two functions are encoded by distinct peptide domains within Hox proteins. This work indicates that Hox proteins execute their critical functions in motor neurons through intrinsic modules that confer distinct specificities and that these activities are central in the genetic network required for motor neuron differentiation.

Polycomb repressive complex 1 activities determine the columnar organization of motor neurons

Polycomb repressive complexes (PRCs) establish and maintain gene repression through chromatin modifications, but their specific roles in cell fate determination events are poorly understood. Here we show an essential role for the PRC1 component Bmi1 in motor neuron (MN) subtype differentiation through dose-dependent effects on Hox gene expression. While Bmi1 is dispensable for generating MNs as a class, it has an essential role in specifying and determining the position of Hox-dependent MN columnar and pool subtypes. These actions are mediated through limiting anterior Hox expression boundaries, functions deployed in post-mitotic MNs, temporally downstream from morphogen gradients. Within the HoxC gene cluster, we found a progressive depletion of PRC-associated marks from rostral to caudal levels of the spinal cord, corresponding to major demarcations of MN subtypes. Selective ablation of Bmi1 elicits a derepression of more posterior Hox genes, leading to a switch in MN fates. Unexpectedly, Hox patterns and MN fates appear to be sensitive to absolute PRC1 activity levels; while reducing Bmi1 switches forelimb lateral motor column (LMC) MNs to a thoracic preganglionic (PGC) identity, elevating Bmi1 expression at thoracic levels converts PGC to LMC MNs. These results suggest that graded PRC1 activities are essential in determining MN topographic organization.

RNA-Seq defines novel genes, RNA processing patterns and enhancer maps for the early stages of nephrogenesis: Hox supergenes

During kidney development the cap mesenchyme progenitor cells both self renew and differentiate into nephrons. The balance between renewal and differentiation determines the final nephron count, which is of considerable medical importance. An important goal is to create a precise genetic definition of the early differentiation of cap mesenchyme progenitors. We used RNA-Seq to transcriptional profile the cap mesenchyme progenitors and their first epithelial derivative, the renal vesicles. The results provide a global view of the changing gene expression program during this key period, defining expression levels for all transcription factors, growth factors, and receptors. The RNA-Seq was performed using two different biochemistries, with one examining only polyadenylated RNA and the other total RNA. This allowed the analysis of noncanonical transcripts, which for many genes were more abundant than standard exonic RNAs. Since a large fraction of enhancers are now known to be transcribed the results also provide global maps of potential enhancers. Further, the RNA-Seq data defined hundreds of novel splice patterns and large numbers of new genes. Particularly striking was the extensive sense/antisense transcription and changing RNA processing complexities of the Hox clusters.