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.

Gene- and Species-Specific Hox mRNA Translation by Ribosome Expansion Segments.

Ribosomes have been suggested to directly control gene regulation, but regulatory roles for ribosomal RNA (rRNA) remain largely unexplored. Expansion segments (ESs) consist of multitudes of tentacle-like rRNA structures extending from the core ribosome in eukaryotes. ESs are remarkably variable in sequence and size across eukaryotic evolution with largely unknown functions. In characterizing ribosome binding to a regulatory element within a Homeobox (Hox) 5' UTR, we identify a modular stem-loop within this element that binds to a single ES, ES9S. Engineering chimeric, "humanized" yeast ribosomes for ES9S reveals that an evolutionary change in the sequence of ES9S endows species-specific binding of Hoxa9 mRNA to the ribosome. Genome editing to site-specifically disrupt the Hoxa9-ES9S interaction demonstrates the functional importance for such selective mRNA-rRNA binding in translation control. Together, these studies unravel unexpected gene regulation directly mediated by rRNA and how ribosome evolution drives translation of critical developmental regulators.

Giant lungfish genome elucidates the conquest of land by vertebrates.

Lungfishes belong to lobe-fined fish (Sarcopterygii) that, in the Devonian period, 'conquered' the land and ultimately gave rise to all land vertebrates, including humans1-3. Here we determine the chromosome-quality genome of the Australian lungfish (Neoceratodus forsteri), which is known to have the largest genome of any animal. The vast size of this genome, which is about 14× larger than that of humans, is attributable mostly to huge intergenic regions and introns with high repeat content (around 90%), the components of which resemble those of tetrapods (comprising mainly long interspersed nuclear elements) more than they do those of ray-finned fish. The lungfish genome continues to expand independently (its transposable elements are still active), through mechanisms different to those of the enormous genomes of salamanders. The 17 fully assembled lungfish macrochromosomes maintain synteny to other vertebrate chromosomes, and all microchromosomes maintain conserved ancient homology with the ancestral vertebrate karyotype. Our phylogenomic analyses confirm previous reports that lungfish occupy a key evolutionary position as the closest living relatives to tetrapods4,5, underscoring the importance of lungfish for understanding innovations associated with terrestrialization. Lungfish preadaptations to living on land include the gain of limb-like expression in developmental genes such as hoxc13 and sall1 in their lobed fins. Increased rates of evolution and the duplication of genes associated with obligate air-breathing, such as lung surfactants and the expansion of odorant receptor gene families (which encode proteins involved in detecting airborne odours), contribute to the tetrapod-like biology of lungfishes. These findings advance our understanding of this major transition during vertebrate evolution.

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.

'Building a perfect body': control of vertebrate organogenesis by PBX-dependent regulatory networks.

Pbx genes encode transcription factors that belong to the TALE (three-amino-acid loop extension) superclass of homeodomain proteins. We have witnessed a surge in information about the roles of this gene family as leading actors in the transcriptional control of development. PBX proteins represent a clear example of how transcription factors can regulate developmental processes by combinatorial properties, acting within multimeric complexes to implement activation or repression of transcription depending on their interaction partners. Here, we revisit long-emphasized functions of PBX transcription factors as cofactors for HOX proteins, major architects of the body plan. We further discuss new knowledge on roles of PBX proteins in different developmental contexts as upstream regulators of Hox genes-as factors that interact with non-HOX proteins and can work independently of HOX-as well as potential pioneer factors. Committed to building a perfect body, PBX proteins govern regulatory networks that direct essential morphogenetic processes and organogenesis in vertebrate development. Perturbations of PBX-dependent networks can cause human congenital disease and cancer.

Hox, homology, and parsimony: An organismal perspective.

Hox genes are important regulators in animal development. They often show a mosaic of conserved (e.g., longitudinal axis patterning) and lineage-specific novel functions (e.g., development of skeletal, sensory, or locomotory systems). Despite extensive research over the past decades, it remains controversial at which node in the animal tree of life the Hox cluster evolved. Its presence already in the last common metazoan ancestor has been proposed, although the genomes of both putative earliest extant metazoan offshoots, the ctenophores and the poriferans, are devoid of Hox sequences. The lack of Hox genes in the supposedly "simple"-built poriferans and their low number in cnidarians and the basally branching bilaterians, the xenacoelomorphs, seems to support the classical notion that the number of Hox genes is correlated with the degree of animal complexity. However, the 4-fold increase of the Hox cluster in xiphosurans, a basally branching chelicerate clade, as well as the situation in some teleost fishes that show a multitude of Hox genes compared to, e.g., human, demonstrates, that there is no per se direct correlation between organismal complexity and Hox number. Traditional approaches have tried to base homology on the morphological level on shared expression profiles of individual genes, but recent data have shown that, in particular with respect to Hox and other regulatory genes, complex gene-gene interactions rather than expression signatures of individual genes alone are responsible for shaping morphological traits during ontogeny. Accordingly, for sound homology assessments and reconstructions of character evolution on organ system level, additional independent datasets (e.g., morphological, developmental) need to be included in any such analyses. If supported by solid data, proposed structural homology should be regarded as valid and not be rejected solely on the grounds of non-parsimonious distribution of the character over a given phylogenetic topology.

Hox genes are crucial regulators of periosteal stem cell identity.

Periosteal stem and progenitor cells (PSPCs) are major contributors to bone maintenance and repair. Deciphering the molecular mechanisms that regulate their function is crucial for the successful generation and application of future therapeutics. Here, we pinpoint Hox transcription factors as necessary and sufficient for periosteal stem cell function. Hox genes are transcriptionally enriched in periosteal stem cells and their overexpression in more committed progenitors drives reprogramming to a naïve, self-renewing stem cell-like state. Crucially, individual Hox family members are expressed in a location-specific manner and their stem cell-promoting activity is only observed when the Hox gene is matched to the anatomical origin of the PSPC, demonstrating a role for the embryonic Hox code in adult stem cells. Finally, we demonstrate that Hoxa10 overexpression partially restores the age-related decline in fracture repair. Together, our data highlight the importance of Hox genes as key regulators of PSPC identity in skeletal homeostasis and repair.

Hox genes in development and beyond.

Hox genes encode evolutionarily conserved transcription factors that are essential for the proper development of bilaterian organisms. Hox genes are unique because they are spatially and temporally regulated during development in a manner that is dictated by their tightly linked genomic organization. Although their genetic function during embryonic development has been interrogated, less is known about how these transcription factors regulate downstream genes to direct morphogenetic events. Moreover, the continued expression and function of Hox genes at postnatal and adult stages highlights crucial roles for these genes throughout the life of an organism. Here, we provide an overview of Hox genes, highlighting their evolutionary history, their unique genomic organization and how this impacts the regulation of their expression, what is known about their protein structure, and their deployment in development and beyond.

Origin and elaboration of a major evolutionary transition in individuality.

Obligate endosymbiosis, in which distantly related species integrate to form a single replicating individual, represents a major evolutionary transition in individuality1-3. Although such transitions are thought to increase biological complexity1,2,4-6, the evolutionary and developmental steps that lead to integration remain poorly understood. Here we show that obligate endosymbiosis between the bacteria Blochmannia and the hyperdiverse ant tribe Camponotini7-11 originated and also elaborated through radical alterations in embryonic development, as compared to other insects. The Hox genes Abdominal A (abdA) and Ultrabithorax (Ubx)-which, in arthropods, normally function to differentiate abdominal and thoracic segments after they form-were rewired to also regulate germline genes early in development. Consequently, the mRNAs and proteins of these Hox genes are expressed maternally and colocalize at a subcellular level with those of germline genes in the germplasm and three novel locations in the freshly laid egg. Blochmannia bacteria then selectively regulate these mRNAs and proteins to make each of these four locations functionally distinct, creating a system of coordinates in the embryo in which each location performs a different function to integrate Blochmannia into the Camponotini. Finally, we show that the capacity to localize mRNAs and proteins to new locations in the embryo evolved before obligate endosymbiosis and was subsequently co-opted by Blochmannia and Camponotini. This pre-existing molecular capacity converged with a pre-existing ecological mutualism12,13 to facilitate both the horizontal transfer10 and developmental integration of Blochmannia into Camponotini. Therefore, the convergence of pre-existing molecular capacities and ecological interactions-as well as the rewiring of highly conserved gene networks-may be a general feature that facilitates the origin and elaboration of major transitions in individuality.

Cicero Predicts cis-Regulatory DNA Interactions from Single-Cell Chromatin Accessibility Data.

Linking regulatory DNA elements to their target genes, which may be located hundreds of kilobases away, remains challenging. Here, we introduce Cicero, an algorithm that identifies co-accessible pairs of DNA elements using single-cell chromatin accessibility data and so connects regulatory elements to their putative target genes. We apply Cicero to investigate how dynamically accessible elements orchestrate gene regulation in differentiating myoblasts. Groups of Cicero-linked regulatory elements meet criteria of "chromatin hubs"-they are enriched for physical proximity, interact with a common set of transcription factors, and undergo coordinated changes in histone marks that are predictive of changes in gene expression. Pseudotemporal analysis revealed that most DNA elements remain in chromatin hubs throughout differentiation. A subset of elements bound by MYOD1 in myoblasts exhibit early opening in a PBX1- and MEIS1-dependent manner. Our strategy can be applied to dissect the architecture, sequence determinants, and mechanisms of cis-regulation on a genome-wide scale.

A region of the human HOXD cluster that confers polycomb-group responsiveness.

Polycomb group (PcG) proteins are essential for accurate axial body patterning during embryonic development. PcG-mediated repression is conserved in metazoans and is targeted in Drosophila by Polycomb response elements (PREs). However, targeting sequences in humans have not been described. While analyzing chromatin architecture in the context of human embryonic stem cell (hESC) differentiation, we discovered a 1.8kb region between HOXD11 and HOXD12 (D11.12) that is associated with PcG proteins, becomes nuclease hypersensitive, and then shows alteration in nuclease sensitivity as hESCs differentiate. The D11.12 element repressed luciferase expression from a reporter construct and full repression required a highly conserved region and YY1 binding sites. Furthermore, repression was dependent on the PcG proteins BMI1 and EED and a YY1-interacting partner, RYBP. We conclude that D11.12 is a Polycomb-dependent regulatory region with similarities to Drosophila PREs, indicating conservation in the mechanisms that target PcG function in mammals and flies.

Visualizing DNA folding and RNA in embryos at single-cell resolution.

The establishment of cell types during development requires precise interactions between genes and distal regulatory sequences. We have a limited understanding of how these interactions look in three dimensions, vary across cell types in complex tissue, and relate to transcription. Here we describe optical reconstruction of chromatin architecture (ORCA), a method that can trace the DNA path in single cells with nanoscale accuracy and genomic resolution reaching two kilobases. We used ORCA to study a Hox gene cluster in cryosectioned Drosophila embryos and labelled around 30 RNA species in parallel. We identified cell-type-specific physical borders between active and Polycomb-repressed DNA, and unexpected Polycomb-independent borders. Deletion of Polycomb-independent borders led to ectopic enhancer-promoter contacts, aberrant gene expression, and developmental defects. Together, these results illustrate an approach for high-resolution, single-cell DNA domain analysis in vivo, identify domain structures that change with cell identity, and show that border elements contribute to the formation of physical domains in Drosophila.

Nervous system-wide analysis of Hox regulation of terminal neuronal fate specification in Caenorhabditis elegans.

Hox genes encode evolutionarily conserved transcription factors that specify regional identities along the anterior-posterior (A-P) axis. Although some Hox genes are known to regulate the differentiation of certain neurons, to what extent Hox genes are involved in the terminal specification of the entire nervous system is unclear. Here, we systematically mapped the expression of all six Hox genes in C. elegans nervous system and found Hox expression in 97 (32%) of the 302 neurons in adult hermaphrodites. Our results are generally consistent with previous high-throughput expression mapping and single-cell transcriptomic studies. Detailed analysis of the fate markers for these neurons revealed that Hox genes regulate the differentiation of 29 (25%) of the 118 classes of C. elegans neurons. Hox genes not only regulate the specification of terminal neuronal fates through multiple mechanisms but also control subtype diversification along the A-P axis. The widespread involvement of Hox genes in neuronal differentiation indicates their roles in establishing complex nervous systems.

A novel post-developmental role of the Hox genes underlies normal adult behavior.

The molecular mechanisms underlying the stability of mature neurons and neural circuits are poorly understood. Here we explore this problem and discover that the Hox genes are a component of the genetic program that maintains normal neural function in adult Drosophila. We show that post-developmental downregulation of the Hox gene Ultrabithorax (Ubx) in adult neurons leads to substantial anomalies in flight. Mapping the cellular basis of these effects reveals that Ubx is required within a subset of dopaminergic neurons, and cell circuitry analyses in combination with optogenetics allow us to link these dopaminergic neurons to flight control. Functional imaging experiments show that Ubx is necessary for normal dopaminergic activity, and neuron-specific RNA-sequencing defines two previously uncharacterized ion channel-encoding genes as potential mediators of Ubx behavioral roles. Our study thus reveals a novel role of the Hox system in controlling adult behavior and neural function. Based on the broad evolutionary conservation of the Hox system across distantly related animal phyla, we predict that the Hox genes might play neurophysiological roles in adult forms of other species, including humans.

The anterior Hox gene ceh-13 and elt-1/GATA activate the posterior Hox genes nob-1 and php-3 to specify posterior lineages in the C. elegans embryo.

Hox transcription factors play a conserved role in specifying positional identity during animal development, with posterior Hox genes typically repressing the expression of more anterior Hox genes. Here, we dissect the regulation of the posterior Hox genes nob-1 and php-3 in the nematode C. elegans. We show that nob-1 and php-3 are co-expressed in gastrulation-stage embryos in cells that previously expressed the anterior Hox gene ceh-13. This expression is controlled by several partially redundant transcriptional enhancers. These enhancers act in a ceh-13-dependant manner, providing a striking example of an anterior Hox gene positively regulating a posterior Hox gene. Several other regulators also act positively through nob-1/php-3 enhancers, including elt-1/GATA, ceh-20/ceh-40/Pbx, unc-62/Meis, pop-1/TCF, ceh-36/Otx, and unc-30/Pitx. We identified defects in both cell position and cell division patterns in ceh-13 and nob-1;php-3 mutants, suggesting that these factors regulate lineage identity in addition to positional identity. Together, our results highlight the complexity and flexibility of Hox gene regulation and function and the ability of developmental transcription factors to regulate different targets in different stages of development.

Widespread employment of conserved C. elegans homeobox genes in neuronal identity specification.

Homeobox genes are prominent regulators of neuronal identity, but the extent to which their function has been probed in animal nervous systems remains limited. In the nematode Caenorhabditis elegans, each individual neuron class is defined by the expression of unique combinations of homeobox genes, prompting the question of whether each neuron class indeed requires a homeobox gene for its proper identity specification. We present here progress in addressing this question by extending previous mutant analysis of homeobox gene family members and describing multiple examples of homeobox gene function in different parts of the C. elegans nervous system. To probe homeobox function, we make use of a number of reporter gene tools, including a novel multicolor reporter transgene, NeuroPAL, which permits simultaneous monitoring of the execution of multiple differentiation programs throughout the entire nervous system. Using these tools, we add to the previous characterization of homeobox gene function by identifying neuronal differentiation defects for 14 homeobox genes in 24 distinct neuron classes that are mostly unrelated by location, function and lineage history. 12 of these 24 neuron classes had no homeobox gene function ascribed to them before, while in the other 12 neuron classes, we extend the combinatorial code of transcription factors required for specifying terminal differentiation programs. Furthermore, we demonstrate that in a particular lineage, homeotic identity transformations occur upon loss of a homeobox gene and we show that these transformations are the result of changes in homeobox codes. Combining the present with past analyses, 113 of the 118 neuron classes of C. elegans are now known to require a homeobox gene for proper execution of terminal differentiation programs. Such broad deployment indicates that homeobox function in neuronal identity specification may be an ancestral feature of animal nervous systems.

SET1A/COMPASS and shadow enhancers in the regulation of homeotic gene expression.

The homeotic (Hox) genes are highly conserved in metazoans, where they are required for various processes in development, and misregulation of their expression is associated with human cancer. In the developing embryo, Hox genes are activated sequentially in time and space according to their genomic position within Hox gene clusters. Accumulating evidence implicates both enhancer elements and noncoding RNAs in controlling this spatiotemporal expression of Hox genes, but disentangling their relative contributions is challenging. Here, we identify two cis-regulatory elements (E1 and E2) functioning as shadow enhancers to regulate the early expression of the HoxA genes. Simultaneous deletion of these shadow enhancers in embryonic stem cells leads to impaired activation of HoxA genes upon differentiation, while knockdown of a long noncoding RNA overlapping E1 has no detectable effect on their expression. Although MLL/COMPASS (complex of proteins associated with Set1) family of histone methyltransferases is known to activate transcription of Hox genes in other contexts, we found that individual inactivation of the MLL1-4/COMPASS family members has little effect on early Hox gene activation. Instead, we demonstrate that SET1A/COMPASS is required for full transcriptional activation of multiple Hox genes but functions independently of the E1 and E2 cis-regulatory elements. Our results reveal multiple regulatory layers for Hox genes to fine-tune transcriptional programs essential for development.

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.

Multi-axial self-organization properties of mouse embryonic stem cells into gastruloids.

The emergence of multiple axes is an essential element in the establishment of the mammalian body plan. This process takes place shortly after implantation of the embryo within the uterus and relies on the activity of gene regulatory networks that coordinate transcription in space and time. Whereas genetic approaches have revealed important aspects of these processes1, a mechanistic understanding is hampered by the poor experimental accessibility of early post-implantation stages. Here we show that small aggregates of mouse embryonic stem cells (ESCs), when stimulated to undergo gastrulation-like events and elongation in vitro, can organize a post-occipital pattern of neural, mesodermal and endodermal derivatives that mimic embryonic spatial and temporal gene expression. The establishment of the three major body axes in these 'gastruloids'2,3 suggests that the mechanisms involved are interdependent. Specifically, gastruloids display the hallmarks of axial gene regulatory systems as exemplified by the implementation of collinear Hox transcriptional patterns along an extending antero-posterior axis. These results reveal an unanticipated self-organizing capacity of aggregated ESCs and suggest that gastruloids could be used as a complementary system to study early developmental events in the mammalian embryo.

hox gene expression predicts tetrapod-like axial regionalization in the skate, Leucoraja erinacea.

The axial skeleton of tetrapods is organized into distinct anteroposterior regions of the vertebral column (cervical, trunk, sacral, and caudal), and transitions between these regions are determined by colinear anterior expression boundaries of Hox5/6, -9, -10, and -11 paralogy group genes within embryonic paraxial mesoderm. Fishes, conversely, exhibit little in the way of discrete axial regionalization, and this has led to scenarios of an origin of Hox-mediated axial skeletal complexity with the evolutionary transition to land in tetrapods. Here, combining geometric morphometric analysis of vertebral column morphology with cell lineage tracing of hox gene expression boundaries in developing embryos, we recover evidence of at least five distinct regions in the vertebral skeleton of a cartilaginous fish, the little skate (Leucoraja erinacea). We find that skate embryos exhibit tetrapod-like anteroposterior nesting of hox gene expression in their paraxial mesoderm, and we show that anterior expression boundaries of hox5/6, hox9, hox10, and hox11 paralogy group genes predict regional transitions in the differentiated skate axial skeleton. Our findings suggest that hox-based axial skeletal regionalization did not originate with tetrapods but rather has a much deeper evolutionary history than was previously appreciated.