Latent developmental potential to form limb-like skeletal structures in zebrafish.

Changes in appendage structure underlie key transitions in vertebrate evolution. Addition of skeletal elements along the proximal-distal axis facilitated critical transformations, including the fin-to-limb transition that permitted generation of diverse modes of locomotion. Here, we identify zebrafish mutants that form supernumerary long bones in their pectoral fins. These new bones integrate into musculature, form joints, and articulate with neighboring elements. This phenotype is caused by activating mutations in previously unrecognized regulators of appendage patterning, vav2 and waslb, that function in a common pathway. This pathway is required for appendage development across vertebrates, and loss of Wasl in mice causes defects similar to those seen in murine Hox mutants. Concordantly, formation of supernumerary bones requires Hox11 function, and mutations in the vav2/wasl pathway drive enhanced expression of hoxa11b, indicating developmental homology with the forearm. Our findings reveal a latent, limb-like pattern ability in fins that is activated by simple genetic perturbation.

Members of an array of zinc-finger proteins specify distinct Hox chromatin boundaries.

Partitioning of repressive from actively transcribed chromatin in mammalian cells fosters cell-type-specific gene expression patterns. While this partitioning is reconstructed during differentiation, the chromatin occupancy of the key insulator, CCCTC-binding factor (CTCF), is unchanged at the developmentally important Hox clusters. Thus, dynamic changes in chromatin boundaries must entail other activities. Given its requirement for chromatin loop formation, we examined cohesin-based chromatin occupancy without known insulators, CTCF and Myc-associated zinc-finger protein (MAZ), and identified a family of zinc-finger proteins (ZNFs), some of which exhibit tissue-specific expression. Two such ZNFs foster chromatin boundaries at the Hox clusters that are distinct from each other and from MAZ. PATZ1 was critical to the thoracolumbar boundary in differentiating motor neurons and mouse skeleton, while ZNF263 contributed to cervicothoracic boundaries. We propose that these insulating activities act with cohesin, alone or combinatorially, with or without CTCF, to implement precise positional identity and cell fate during development.

Dissecting CTCF site function in a tense HoxD locus.

In this issue of Genes & Development, Amândio and colleagues (pp. 1490-1509) dissect the function of a cluster of several CTCF binding sites in the HoxD cluster by iterative deletions in mice. They found additive functions for some, and intriguingly found that some sites function as insulators, while others function as anchors for enhancer-promoter interactions. These functions vary depending on developmental context. The work provides new insights into the roles played by CTCF in regulating developmental patterns and 3D chromatin organization.

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.

A six-amino-acid motif is a major determinant in functional evolution of HOX1 proteins.

Gene duplication and divergence is a major driver in the emergence of evolutionary novelties. How variations in amino acid sequences lead to loss of ancestral activity and functional diversification of proteins is poorly understood. We used cross-species functional analysis of Drosophila Labial and its mouse HOX1 orthologs (HOXA1, HOXB1, and HOXD1) as a paradigm to address this issue. Mouse HOX1 proteins display low (30%) sequence similarity with Drosophila Labial. However, substituting endogenous Labial with the mouse proteins revealed that HOXA1 has retained essential ancestral functions of Labial, while HOXB1 and HOXD1 have diverged. Genome-wide analysis demonstrated similar DNA-binding patterns of HOXA1 and Labial in mouse cells, while HOXB1 binds to distinct targets. Compared with HOXB1, HOXA1 shows an enrichment in co-occupancy with PBX proteins on target sites and exists in the same complex with PBX on chromatin. Functional analysis of HOXA1-HOXB1 chimeric proteins uncovered a novel six-amino-acid C-terminal motif (CTM) flanking the homeodomain that serves as a major determinant of ancestral activity. In vitro DNA-binding experiments and structural prediction show that CTM provides an important domain for interaction of HOXA1 proteins with PBX. Our findings show that small changes outside of highly conserved DNA-binding regions can lead to profound changes in protein function.

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.

Hox11-expressing interstitial cells contribute to adult skeletal muscle at homeostasis.

Interstitial stromal cells play critical roles in muscle development, regeneration and repair and we have previously reported that Hoxa11 and Hoxd11 are expressed in the interstitial cells of muscles attached to the zeugopod, and are crucial for the proper embryonic patterning of these muscles. Hoxa11eGFP expression continues in a subset of muscle interstitial cells through adult stages. The induction of Hoxa11-CreERT2-mediated lineage reporting (Hoxa11iTom) at adult stages in mouse results in lineage induction only in the interstitial cells. However, Hoxa11iTom+ cells progressively contribute to muscle fibers at subsequent stages. The contribution to myofibers exceeds parallel Pax7-CreERT2-mediated lineage labeling. Nuclear-specific lineage labeling demonstrates that Hoxa11-expressing interstitial cells contribute nuclear contents to myofibers. Crucially, at no point after Hoxa11iTom induction are satellite cells lineage labeled. When examined in vitro, isolated Hoxa11iTom+ interstitial cells are not capable of forming myotubes, but Hoxa11iTom+ cells can contribute to differentiating myotubes, supporting Hox-expressing interstitial cells as a new population of muscle progenitors, but not stem cells. This work adds to a small but growing body of evidence that supports a satellite cell-independent source of muscle tissue in vivo.

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.

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.

Evolutionarily conserved inhibitory uORFs sensitize Hox mRNA translation to start codon selection stringency.

Translation start site selection in eukaryotes is influenced by context nucleotides flanking the AUG codon and by levels of the eukaryotic translation initiation factors eIF1 and eIF5. In a search of mammalian genes, we identified five homeobox (Hox) gene paralogs initiated by AUG codons in conserved suboptimal context as well as 13 Hox genes that contain evolutionarily conserved upstream open reading frames (uORFs) that initiate at AUG codons in poor sequence context. An analysis of published cap analysis of gene expression sequencing (CAGE-seq) data and generated CAGE-seq data for messenger RNAs (mRNAs) from mouse somites revealed that the 5' leaders of Hox mRNAs of interest contain conserved uORFs, are generally much shorter than reported, and lack previously proposed internal ribosome entry site elements. We show that the conserved uORFs inhibit Hox reporter expression and that altering the stringency of start codon selection by overexpressing eIF1 or eIF5 modulates the expression of Hox reporters. We also show that modifying ribosome homeostasis by depleting a large ribosomal subunit protein or treating cells with sublethal concentrations of puromycin leads to lower stringency of start codon selection. Thus, altering global translation can confer gene-specific effects through altered start codon selection stringency.

Embryonic timing, axial stem cells, chromatin dynamics, and the Hox clock.

Collinear regulation of Hox genes in space and time has been an outstanding question ever since the initial work of Ed Lewis in 1978. Here we discuss recent advances in our understanding of this phenomenon in relation to novel concepts associated with large-scale regulation and chromatin structure during the development of both axial and limb patterns. We further discuss how this sequential transcriptional activation marks embryonic stem cell-like axial progenitors in mammals and, consequently, how a temporal genetic system is further translated into spatial coordinates via the fate of these progenitors. In this context, we argue the benefit and necessity of implementing this unique mechanism as well as the difficulty in evolving an alternative strategy to deliver this critical positional information.

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.

Somite development and regionalisation of the vertebral axial skeleton.

A critical stage in the development of all vertebrate embryos is the generation of the body plan and its subsequent patterning and regionalisation along the main anterior-posterior axis. This includes the formation of the vertebral axial skeleton. Its organisation begins during early embryonic development with the periodic formation of paired blocks of mesoderm tissue called somites. Here, we review axial patterning of somites, with a focus on studies using amniote model systems - avian and mouse. We summarise the molecular and cellular mechanisms that generate paraxial mesoderm and review how the different anatomical regions of the vertebral column acquire their specific identity and thus shape the body plan. We also discuss the generation of organoids and embryo-like structures from embryonic stem cells, which provide insights regarding axis formation and promise to be useful for disease modelling.

Identification of in vivo Hox13-binding sites reveals an essential locus controlling zebrafish brachyury expression.

During early embryogenesis, the vertebrate embryo extends from anterior to posterior because of the progressive addition of cells from a posteriorly localized neuromesodermal progenitor (NMp) population. An autoregulatory loop between Wnt and Brachyury/Tbxt is required for NMps to retain mesodermal potential and, hence, normal axis development. We recently showed that Hox13 genes help to support body axis formation and to maintain the autoregulatory loop, although the direct Hox13 target genes were unknown. Here, using a new method for identifying in vivo transcription factor-binding sites, we identified more than 500 potential Hox13 target genes in zebrafish. Importantly, we found two highly conserved Hox13-binding elements far from the tbxta transcription start site that also contain a conserved Tcf7/Lef1 (Wnt response) site. We show that the proximal of the two elements is sufficient to confer somitogenesis-stage expression to a tbxta promoter that, on its own, only drives NMp expression during gastrulation. Importantly, elimination of this proximal element produces shortened embryos due to aberrant formation of the most posterior somites. Our study provides a potential direct connection between Hox13 and regulation of the Wnt/Brachyury loop.

Dynamic extrinsic pacing of the HOX clock in human axial progenitors controls motor neuron subtype specification.

Rostro-caudal patterning of vertebrates depends on the temporally progressive activation of HOX genes within axial stem cells that fuel axial embryo elongation. Whether the pace of sequential activation of HOX genes, the 'HOX clock', is controlled by intrinsic chromatin-based timing mechanisms or by temporal changes in extrinsic cues remains unclear. Here, we studied HOX clock pacing in human pluripotent stem cell-derived axial progenitors differentiating into diverse spinal cord motor neuron subtypes. We show that the progressive activation of caudal HOX genes is controlled by a dynamic increase in FGF signaling. Blocking the FGF pathway stalled induction of HOX genes, while a precocious increase of FGF, alone or with GDF11 ligand, accelerated the HOX clock. Cells differentiated under accelerated HOX induction generated appropriate posterior motor neuron subtypes found along the human embryonic spinal cord. The pacing of the HOX clock is thus dynamically regulated by exposure to secreted cues. Its manipulation by extrinsic factors provides synchronized access to multiple human neuronal subtypes of distinct rostro-caudal identities for basic and translational applications.This article has an associated 'The people behind the papers' interview.

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.

Discovery of seven hox genes in zebrafish thrombopoiesis.

Thrombopoiesis is the production of platelets from megakaryocytes in the bone marrow of mammals. In fish, thrombopoiesis involves the formation of thrombocytes without megakaryocyte-like precursors but derived from erythrocyte thrombocyte bi-functional precursor cells. One unique feature of thrombocyte differentiation involves the maturation of young thrombocytes in circulation. In this study, we investigated the role of hox genes in zebrafish thrombopoiesis to model platelet production. We selected hoxa10b, hoxb2a, hoxc5a, hoxd3a, and hoxc11b from thrombocyte RNA expression data, and checked whether they are expressed in young or mature thrombocytes. We found hoxa10b, hoxb2a, hoxc5a, and hoxd3a were expressed in both young and mature thrombocytes and hoxc11b was expressed in only young thrombocytes. We then performed knockdowns of these 5 hox genes and found hoxc11b knockdown resulted in thrombocytosis and the rest showed thrombocytopenia. To identify hox genes that could have been missed by the above datasets, we performed knockdowns 47 hox genes in the zebrafish genome and found hoxa9a, and hoxb1a knockdowns resulted in thrombocytopenia and they were expressed in both young and mature thrombocytes. In conclusion, our comprehensive knockdown study identified Hoxa10b, Hoxb2a, Hoxc5a, Hoxd3a, Hoxa9a, and Hoxb1a, as positive regulators and Hoxc11b, as a negative regulator for thrombocyte development.

Distinct roles of the Hox genes Ultrabithorax and abdominal-A in scorpionfly embryonic proleg development.

The abdominal appendages of larval insects have a complex evolutionary history of gain and loss, but the regulatory mechanisms underlying the abdominal appendage development remain largely unclear. Here, we investigated the embryogenesis of abdominal prolegs in the scorpionfly Panorpa liui Hua (Mecoptera: Panorpidae) using in situ hybridization and parental RNA interference. The results show that RNAi-mediated knockdown of Ultrabithorax (Ubx) led to a homeotic transformation of the first abdominal segment (A1) into the third thoracic segment (T3) and changed the distributions of the downstream target Distal-less (Dll) expression but did not affect the expression levels of Dll. Knockdown of abdominal-A (abd-A) resulted in malformed segments, abnormal prolegs and disrupted Dll expression. The results demonstrate that the gene Ubx maintains an ancestral role of modulating A1 appendage fate without preventing Dll initiation, and a secondary adaptation of abd-A evolves the ability to specify abdominal segments and proleg identity. We conclude that changes in abdominal Hox gene expression and their target genes regulate abdominal appendage morphology during the evolutionary course of holometabolous larvae.

Irreducible Complexity of Hox Gene: Path to the Canonical Function of the Hox Cluster.

The evolution of major taxa is often associated with the emergence of new gene families. In all multicellular animals except sponges and comb jellies, the genomes contain Hox genes, which are crucial regulators of development. The canonical function of Hox genes involves colinear patterning of body parts in bilateral animals. This general function is implemented through complex, precisely coordinated mechanisms, not all of which are evolutionarily conserved and fully understood. We suggest that the emergence of this regulatory complexity was preceded by a stage of cooperation between more ancient morphogenetic programs or their individual elements. Footprints of these programs may be present in modern animals to execute non-canonical Hox functions. Non-canonical functions of Hox genes are involved in maintaining terminal nerve cell specificity, autophagy, oogenesis, pre-gastrulation embryogenesis, vertical signaling, and a number of general biological processes. These functions are realized by the basic properties of homeodomain protein and could have triggered the evolution of ParaHoxozoa and Nephrozoa subsequently. Some of these non-canonical Hox functions are discussed in our review.

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.