Hox Expression Patterns Research Articles

The E3 ubiquitin ligase RNF220 maintains hindbrain Hox expression patterns through regulation of WDR5 stability.

The spatial and temporal linear expression of Hox genes establishes a regional Hox code, which is crucial for the antero-posterior (A-P) patterning, segmentation, and neuronal circuit development of the hindbrain. RNF220, an E3 ubiquitin ligase, is widely involved in neural development via targeting of multiple substrates. Here, we found that the expression of Hox genes in the pons was markedly up-regulated at the late developmental stage (post-embryonic day E15.5) in Rnf220-/- and Rnf220+/- mouse embryos. Single-nucleus RNA sequencing (RNA-seq) analysis revealed different Hox de-repression profiles in different groups of neurons, including the pontine nuclei (PN). The Hox pattern was disrupted and the neural circuits were affected in the PN of Rnf220+/- mice. We showed that this phenomenon was mediated by WDR5, a key component of the TrxG complex, which can be polyubiquitinated and degraded by RNF220. Intrauterine injection of WDR5 inhibitor (WDR5-IN-4) and genetic ablation of Wdr5 in Rnf220+/- mice largely recovered the de-repressed Hox expression pattern in the hindbrain. In P19 embryonal carcinoma cells, the retinoic acid-induced Hox expression was further stimulated by Rnf220 knockdown, which can also be rescued by Wdr5 knockdown. In short, our data suggest a new role of RNF220/WDR5 in Hox pattern maintenance and pons development in mice.

The E3 ubiquitin ligase RNF220 maintains hindbrain Hox expression patterns through regulation of WDR5 stability

The spatial and temporal linear expression of Hox genes establishes a regional Hox code, which is crucial for the antero-posterior (A-P) patterning, segmentation, and neuronal circuit development of the hindbrain. RNF220, an E3 ubiquitin ligase, is widely involved in neural development via targeting of multiple substrates. Here, we found that the expression of Hox genes in the pons was markedly up-regulated at the late developmental stage (post-embryonic day E15.5) in Rnf220-/- and Rnf220+/- mouse embryos. Single-nucleus RNA sequencing (RNA-seq) analysis revealed different Hox de-repression profiles in different groups of neurons, including the pontine nuclei (PN). The Hox pattern was disrupted and the neural circuits were affected in the PN of Rnf220+/- mice. We showed that this phenomenon was mediated by WDR5, a key component of the TrxG complex, which can be polyubiquitinated and degraded by RNF220. Intrauterine injection of WDR5 inhibitor (WDR5-IN-4) and genetic ablation of Wdr5 in Rnf220+/- mice largely recovered the de-repressed Hox expression pattern in the hindbrain. In P19 embryonal carcinoma cells, the retinoic acid-induced Hox expression was further stimulated by Rnf220 knockdown, which can also be rescued by Wdr5 knockdown. In short, our data suggest a new role of RNF220/WDR5 in Hox pattern maintenance and pons development in mice.

Plasticity Comparison of Two Stem Cell Sources with Different Hox Gene Expression Profiles in Response to Cobalt Chloride Treatment during Chondrogenic Differentiation.

The limited self-repair capacity of articular cartilage is a challenge for healing injuries. While mesenchymal stem/stromal cells (MSCs) are a promising approach for tissue regeneration, the criteria for selecting a suitable cell source remain undefined. To propose a molecular criterion, dental pulp stem cells (DPSCs) with a Hox-negative expression pattern and bone marrow mesenchymal stromal cells (BMSCs), which actively express Hox genes, were differentiated towards chondrocytes in 3D pellets, employing a two-step protocol. The MSCs' response to preconditioning by cobalt chloride (CoCl2), a hypoxia-mimicking agent, was explored in an assessment of the chondrogenic differentiation's efficiency using morphological, histochemical, immunohistochemical, and biochemical experiments. The preconditioned DPSC pellets exhibited significantly elevated levels of collagen II and glycosaminoglycans (GAGs) and reduced levels of the hypertrophic marker collagen X. No significant effect on GAGs production was observed in the preconditioned BMSC pellets, but collagen II and collagen X levels were elevated. While preconditioning did not modify the ALP specific activity in either cell type, it was notably lower in the DPSCs differentiated pellets compared to their BMSCs counterparts. These results could be interpreted as demonstrating the higher plasticity of DPSCs compared to BMSCs, suggesting the contribution of their unique molecular characteristics, including their negative Hox expression pattern, to promote a chondrogenic differentiation potential. Consequently, DPSCs could be considered compelling candidates for future cartilage cell therapy.

Gene Expression Machine Learning Models Classify Pediatric AML Subtypes with High Performance

Pediatric acute myeloid leukemia (pAML) encompasses over 20 molecular subtypes driven by unique genetic alterations, including hallmark chromosomal rearrangements and less frequently, point mutations or tandem duplications. Collectively, many of these pAML molecular categories are enriched in pediatric populations and are not represented by current classification systems, including the recently updated WHO or ICC. While fusion detection from RNA-Seq-based approaches is robust, many fusion negative subtypes would need to be defined by expression-based approaches as mutation calling from RNA-Seq data is less developed. Nevertheless, this could be challenging for subtypes with similar transcriptional profiles, such as those with shared HOX expression patterns, including NPM1, NUP98r, UBTF, DEK::NUP214, and KMT2A-PTD. To aid in the appropriate molecular classification of pAML, which is crucial for prognosis, we developed and compared three gene expression-based classifiers. A total of 1707 pAML gene expression profiles were mapped and analyzed from three distinct sources (St. Jude = 659; TARGET = 168; AAML1031 = 880). Raw read count data was normalized and scaled to obtain a relative expression value in transcripts per million (TPM), which served as input for feature selection, model training, and testing. Ground truth labels for all 1707 samples were obtained through multi-omics analysis, including whole genome sequencing, to identify fusions and mutations. For validation purposes, the data was stratified by subtypes and split 70/30 into training (n=1187) and testing (n=520). Three machine learning models were selected: random forest, XGboost, and linear support vector machine (SVM). Each sample had gene expression TPM data for 60,754 transcripts, out of which 20,004 transcripts related to protein-coding genes were incorporated for feature selection. Feature selection was performed using a median absolute deviation (MAD) algorithm to select the 3000 transcripts with the highest variability. Top variable 3000 genes were selected to allow for adequate tuning of the number of predictors. Each model was independently trained using stratified cross-validation and Monte-Carlo search for hyperparameter tuning. The best model was selected based on the Matthews Correlation Coefficient (MCC). Each model was tested on the hold-out TPM set with z-score normalization. On the hold-out testing set (n=520), the linear SVM model outperformed the random forest and XGboost models on five performance metrics across all subtypes (sensitivity=0.9577; precision=0.9577; specificity=0.9978; F1=0.96; accuracy=0.9958). The random forest (sensitivity=0.9154; precision=0.9154; specificity=0.9955; F1=0.92; accuracy=0.9915) and XGboost models (sensitivity=0.9231; precision=0.9231; specificity=0.9960; F1=0.92; accuracy=0.9923) also performed well across all subtypes. Although feature selection was shared across all three models, performance within each subtype varied between models. The linear SVM model demonstrated strong performance overall, driven by high specificity in classifying the KMT2Ar subgroup (n=127) and equal sensitivity across the GLIS-rearranged (n=14), GATA1 (n=10), BCL11B (n=7), CBFB::MYH11 (n=60), CEBPA (n=33), and RUNX1::RUNX1T1 (n=78) subtypes. The primary difference in performance between models is the high false positive rate for KMT2Ar and NPM1 (n=55) in the random forest and XGboost models. A preliminary hypothesis for this might be due to the large representation of KMT2Ar and NPM1 in the training data (24.85% and 10.78%, respectively). Synthetic upsampling (SMOTE) for the training dataset (n=1369) counteracts bias towards the majority classes and increases performance for the random forest (sensitivity=0.9327; precision=0.9327; specificity=0.9965; F1=0.93; accuracy=0.9933), XGboost (sensitivity=0.9404; precision=0.9404; specificity=0.9969; F1=0.94; accuracy=0.9940) and linear SVM (sensitivity=0.9615; precision=0.9615; specificity=0.9980; F1=0.96; accuracy=0.9962) models. Conjointly, these models demonstrate the utility and effectiveness of a machine learning approach for classifying pAML samples from transcriptome sequencing data, which may have broad clinical and research utility, especially for fusion negative subtypes.

HOX genes in stem cells: Maintaining cellular identity and regulation of differentiation.

Stem cells display a unique cell type within the body that has the capacity to self-renew and differentiate into specialized cell types. Compared to pluripotent stem cells, adult stem cells (ASC) such as mesenchymal stem cells (MSCs) and hematopoietic stem cells (HSCs) exhibit restricted differentiation capabilities that are limited to cell types typically found in the tissue of origin, which implicates that there must be a certain code or priming determined by the tissue of origin. HOX genes, a subset of homeobox genes encoding transcription factors that are generally repressed in undifferentiated pluripotent stem cells, emerged here as master regulators of cell identity and cell fate during embryogenesis, and in maintaining this positional identity throughout life as well as specifying various regional properties of respective tissues. Concurrently, intricate molecular circuits regulated by diverse stem cell-typical signaling pathways, balance stem cell maintenance, proliferation and differentiation. However, it still needs to be unraveled how stem cell-related signaling pathways establish and regulate ASC-specific HOX expression pattern with different temporal-spatial topography, known as the HOX code. This comprehensive review therefore summarizes the current knowledge of specific ASC-related HOX expression patterns and how these were integrated into stem cell-related signaling pathways. Understanding the mechanism of HOX gene regulation in stem cells may provide new ways to manipulate stem cell fate and function leading to improved and new approaches in the field of regenerative medicine.

Molecular Characterization of HOXA2 and HOXA3 Binding Properties.

The highly conserved HOX homeodomain (HD) transcription factors (TFs) establish the identity of different body parts along the antero–posterior axis of bilaterian animals. Segment diversification and the morphogenesis of different structures is achieved by generating precise patterns of HOX expression along the antero–posterior axis and by the ability of different HOX TFs to instruct unique and specific transcriptional programs. However, HOX binding properties in vitro, characterised by the recognition of similar AT-rich binding sequences, do not account for the ability of different HOX to instruct segment-specific transcriptional programs. To address this problem, we previously compared HOXA2 and HOXA3 binding in vivo. Here, we explore if sequence motif enrichments observed in vivo are explained by binding affinities in vitro. Unexpectedly, we found that the highest enriched motif in HOXA2 peaks was not recognised by HOXA2 in vitro, highlighting the importance of investigating HOX binding in its physiological context. We also report the ability of HOXA2 and HOXA3 to heterodimerise, which may have functional consequences for the HOX patterning function in vivo.

Disruption of CTCF Boundary Remodels Chromatin Domain and Reprograms HOX Gene Transcription in Acute Myeloid Leukemia

The organization of the mammalian genome into separate domains either facilitates or restricts interactions between regulatory elements, e.g. enhancers and promoters, for normal gene regulation. CTCF insulator constrains the temporal HOX gene expression patterns within the confined chromatin domain for normal development. However, dysregulation of HOX genes is also a common feature of acute myeloid leukemia (AML). The molecular mechanisms of aberrant HOX expression and associated AML pathogenesis remain unclear. Here, we found that CTCF binding site located between HOXA7 and HOXA9 genes (CBS7/9) acts as a chromatin boundary to establish and maintain AML specific aberrant HOX gene expression patterns and is essential for pathogenesis of myeloid malignancies.To unbiasedly and systematically investigate the mechanism underlying HOX gene regulation in AML, we employed a pooled CRISPR-Cas9 genetic knockout library screening to interrogate intergenic regulatory and chromatin boundary elements in all 4 HOX gene loci. Through this screening, we identified the CBS7/9 chromatin boundary located between HOXA7 and HOXA9 genes that is critical for maintaining aberrant expression of posterior HOXA genes including oncogene HOXA9 and HOXA10 in primary AML patients and cells. To further confirm the role of CBS7/9 in HOX gene regulation, two single-guide RNAs (sgRNAs) flanking 47 base pairs (bps) of the core CTCF motif located between HOXA7 and HOXA9 geneswere employed to guild Cas9 for deletion of the CBS7/9 boundary in 3 different subtypes of AML cells containing MLL rearrangements or NPM1C+ mutation. The defected CBS7/9 resulted in a decrease in CTCF binding to the caudal HOXA domain and an inhibition of posterior HOXA gene expression.To further delineate the role of CBS7/9 chromatin boundary in HOX gene organization and regulation, we then carried out genome wide ChIP-seq, ATAC-seq, 4C, and RNA-seq analysis to examine the changes in chromatin domain organization and corresponding gene expression patterns compared wild-type control and the CBS7/9 KO AML cells. Disruption of CBS7/9 boundary resulted in a spreading of repressive H3K27me3 domain into the caudal HOXA domain that subsequently impaired AML associated active TAD domain by altering enhancer/promoter chromatin accessibility and reducing ectopic interactions among the posterior HOXA genes. Consistent with the role of CBS7/9 boundary in HOX locus chromatin structure and organization, impairment of CBS7/9 boundary function reduced posterior HOXA gene expression and altered myeloid specific transcriptome profiles important for pathogenesis of myeloid malignancies suggesting that the CBS7/9 boundary mediated ectopic HOX expression patterns might be essential for pathogenesis of myeloid malignancies. To further test this possibility, the CBS 7/9 boundary was deleted in both MLL rearranged AML cells and primary patient samples by the CRISPR-Cas9 mediated genome editing and the edited cells were transplanted into NSG mice for the development of human AML in mice. Disruption of CBS chromatin boundary in the HOXA locus reduced human leukemic blast development and invasion as well as enhanced survival of both transplanted AML cell xenograft and patient derived xenograft (PDX) mouse models. Thus, our finding demonstrated that CTCF boundary not only constrains normal gene expression program for memory of cell identity, but also plays a role in maintaining oncogenic transcription program for leukemic transformation. Alteration of CTCF mediated topological boundaries in the oncogene loci provides therapeutic approaches for targeting oncogenic expression program and leukemogenesis. DisclosuresNo relevant conflicts of interest to declare.

Development and evolution of regionalization within the avian axial column

AbstractThe origin of birds from their terrestrial antecedents was accompanied by a wholesale transformation of their skeleton as they transitioned from a terrestrial to aerial realm. Part of this dramatic transformation is the reduction of separate vertebral elements into regional fusions to limit axial flexibility. This is partially mirrored within the development of the axial column, with regions of the axial column experiencing increasing morphological modularity and the loss of skeletal elements through vertebral fusions. Using a detailed description of the morphological development of the axial column in the model domestic chicken, Gallus gallus domesticus, we present a map of axial ossification based on discrete characters. Delays in ossification are found to occur in conjunction with the formation of fusions. Our study shows that the pattern and sequence of fusion and ossification during development may reflect the presence of independent modules as subsets within the typical regions of the avian axial column. Interestingly, few of these fusion modules correspond to the initial axial Hox expression patterns, suggesting another patterning mechanism is driving axial fusion regionalization. Additionally, two regions of fusion are discovered in the synsacrum. The anterior region of seven fused synsacrals may correspond to the non-ornithuran pygostylian synsacrum of the same number of vertebrae.

An atlas of anterior hox gene expression in the embryonic sea lamprey head: Hox-code evolution in vertebrates

In the hindbrain and the adjacent cranial neural crest (NC) cells of jawed vertebrates (gnathostomes), nested and segmentally-restricted domains of Hox gene expression provide a combinatorial Hox-code for specifying regional properties during head development. Extant jawless vertebrates, such as the sea lamprey (Petromyzon marinus), can provide insights into the evolution and diversification of this Hox-code in vertebrates. There is evidence for gnathostome-like spatial patterns of Hox expression in lamprey; however, the expression domains of the majority of lamprey hox genes from paralogy groups (PG) 1-4 are yet to be characterized, so it is unknown whether they are coupled to hindbrain segments (rhombomeres) and NC. In this study, we systematically describe the spatiotemporal expression of all 14 sea lamprey hox genes from PG1-PG4 in the developing hindbrain and pharynx to investigate the extent to which their expression conforms to the archetypal gnathostome hindbrain and pharyngeal hox-codes. We find many similarities in Hox expression between lamprey and gnathostome species, particularly in rhombomeric domains during hindbrain segmentation and in the cranial neural crest, enabling inference of aspects of Hox expression in the ancestral vertebrate embryonic head. These data are consistent with the idea that a Hox regulatory network underlying hindbrain segmentation is a pan vertebrate trait. We also reveal differences in hindbrain domains at later stages, as well as expression in the endostyle and in pharyngeal arch (PA) 1 mesoderm. Our analysis suggests that many Hox expression domains that are observed in extant gnathostomes were present in ancestral vertebrates but have been partitioned differently across Hox clusters in gnathostome and cyclostome lineages after duplication.

Long noncoding RNA HOTTIP cooperates with CCCTC-binding factor to coordinate HOXA gene expression

The spatiotemporal control of HOX gene expression is dependent on positional identity and often correlated to their genomic location within each loci. Maintenance of HOX expression patterns is under complex transcriptional and epigenetic regulation, which is not well understood. Here we demonstrate that HOTTIP, a lincRNA transcribed from the 5′ edge of the HOXA locus, physically associates with the CCCTC-binding factor (CTCF) that serves as an insulator by organizing HOXA cluster into disjoint domains, to cooperatively maintain the chromatin modifications of HOXA genes and thus coordinate the transcriptional activation of distal HOXA genes in human foreskin fibroblasts. Our results reveal the functional connection of HOTTIP and CTCF, and shed light on lincRNAs in gene activation and CTCF mediated chromatin organization.

Motor neurons with limb-innervating character in the cervical spinal cord are sculpted by apoptosis based on the Hox code in chick embryo.

In the developing chick embryo, a certain population of motor neurons (MNs) in the non-limb-innervating cervical spinal cord undergoes apoptosis between embryonic days 4 and 5. However, the characteristics of these apoptotic MNs remain undefined. Here, by examining the spatiotemporal profiles of apoptosis and MN subtype marker expression in normal or apoptosis-inhibited chick embryos, we found that this apoptotic population is distinguishable by Foxp1 expression. When apoptosis was inhibited, the Foxp1+ MNs survived and showed characteristics of lateral motor column (LMC) neurons, which are of a limb-innervating subtype, suggesting that cervical Foxp1+ MNs are the rostral continuation of the LMC. Knockdown and misexpression of Foxp1 did not affect apoptosis progression, but revealed the role of Foxp1 in conferring LMC identity on the cervical MNs. Furthermore, ectopic expression of Hox genes that are normally expressed in the brachial region prevented apoptosis, and directed Foxp1+ MNs to LMC neurons at the cervical level. These results indicate that apoptosis in the cervical spinal cord plays a role in sculpting Foxp1+ MNs committed to LMC neurons, depending on the Hox expression pattern.

A Simple Predictive Enhancer Syntax for Hindbrain Patterning Is Conserved in Vertebrate Genomes.

BackgroundDetermining the function of regulatory elements is fundamental for our understanding of development, disease and evolution. However, the sequence features that mediate these functions are often unclear and the prediction of tissue-specific expression patterns from sequence alone is non-trivial. Previous functional studies have demonstrated a link between PBX-HOX and MEIS/PREP binding interactions and hindbrain enhancer activity, but the defining grammar of these sites, if any exists, has remained elusive.ResultsHere, we identify a shared sequence signature (syntax) within a heterogeneous set of conserved vertebrate hindbrain enhancers composed of spatially co-occurring PBX-HOX and MEIS/PREP transcription factor binding motifs. We use this syntax to accurately predict hindbrain enhancers in 89% of cases (67/75 predicted elements) from a set of conserved non-coding elements (CNEs). Furthermore, mutagenesis of the sites abolishes activity or generates ectopic expression, demonstrating their requirement for segmentally restricted enhancer activity in the hindbrain. We refine and use our syntax to predict over 3,000 hindbrain enhancers across the human genome. These sequences tend to be located near developmental transcription factors and are enriched in known hindbrain activating elements, demonstrating the predictive power of this simple model.ConclusionOur findings support the theory that hundreds of CNEs, and perhaps thousands of regions across the human genome, function to coordinate gene expression in the developing hindbrain. We speculate that deeply conserved sequences of this kind contributed to the co-option of new genes into the hindbrain gene regulatory network during early vertebrate evolution by linking patterns of hox expression to downstream genes involved in segmentation and patterning, and evolutionarily newer instances may have continued to contribute to lineage-specific elaboration of the hindbrain.

F9 EC 세포에서 레티노산에 의해 유도되는 Hoxc 유전자의 발현에 히스톤 메틸화가 미치는 영향

Department of Anatomy, Embryology Laboratory, Brain Korea 21 Pl us project for Medical Science, Yonsei University College of Medicine, Seoul 120-752, KoreaReceived December 1, 2014/Revised March 11, 2015/Accepted March 12, 2015Hox genes encode a highly conserved family of homeodomain-containing transcription factors control-ling vertebrate pattern formation along the anteroposterior body axis during embryogenesis. Retinoic acid (RA) is a key morphogen in embryogenesis and a critical regulator of both adult and embryonic cellular activity. Specifically, RA regulates Hox gene expression in mouse- or human-derived embry-onic carcinoma (EC) cells. Histone modification has been reported to play a pivotal role in the process of RA-induced gene expression and cell differentiation. As histone modification is thought to play an essential role in RA-induced Hox gene expression, we examined RA-induced initiation of collinear ex-pression of Hox genes and the corresponding histone modifications in F9 murine embryonic ter-atocarcinoma (EC) cells. Hox expression patterns and histone modifications were analyzed by semi- quantitative RT-PCR, RNA-sequencing, and chromatin immuno-precipitation (ChIP)-PCR analyses. The Hoxc4 gene (D0) was initiated earlier than the Hoxc5 to –c10 genes (D3) upon RA treatment (day 0 [D0], day 1 [D1], and day 3 [D3]). The Hox nonexpressing D0 sample had a strong repressive marker, H3K27me3, than the D1 and D3 samples. In the D1 and D3 samples, reduced enrichment of the H3K27me3 marker was observed in the whole cluster. The active H3K4me3 marker was closely asso-ciated with the collinear expression of Hoxc genes. Thus, the Hoxc4 gene (D1) and all Hoxc genes (D3) expressed H3K4me3 upon transcription activation. In conclusion, these data indicated that removing H3K27me3 and acquiring H3K4me3 regulated RA-induced Hoxc gene collinearity in F9 cells. Keywords :

The RNA-binding protein ELAV regulates Hox RNA processing, expression and function within the Drosophila nervous system.

The regulated head-to-tail expression of Hox genes provides a coordinate system for the activation of specific programmes of cell differentiation according to axial level. Recent work indicates that Hox expression can be regulated via RNA processing but the underlying mechanisms and biological significance of this form of regulation remain poorly understood. Here we explore these issues within the developing Drosophila central nervous system (CNS). We show that the pan-neural RNA-binding protein (RBP) ELAV (Hu antigen) regulates the RNA processing patterns of the Hox gene Ultrabithorax (Ubx) within the embryonic CNS. Using a combination of biochemical, genetic and imaging approaches we demonstrate that ELAV binds to discrete elements within Ubx RNAs and that its genetic removal reduces Ubx protein expression in the CNS leading to the respecification of cellular subroutines under Ubx control, thus defining for the first time a specific cellular role of ELAV within the developing CNS. Artificial provision of ELAV in glial cells (a cell type that lacks ELAV) promotes Ubx expression, suggesting that ELAV-dependent regulation might contribute to cell type-specific Hox expression patterns within the CNS. Finally, we note that expression of abdominal A and Abdominal B is reduced in elav mutant embryos, whereas other Hox genes (Antennapedia) are not affected. Based on these results and the evolutionary conservation of ELAV and Hox genes we propose that the modulation of Hox RNA processing by ELAV serves to adapt the morphogenesis of the CNS to axial level by regulating Hox expression and consequently activating local programmes of neural differentiation.

Comprehensive Analysis Of HOX Gene Expression and DNA Methylation From 189 Primary AMLs Demonstrates Canonical Patterns Associated With Hematopoietic Stem/Progenitors and Recurrent AML Mutations

▪HOX genes encode a family of homeodomain transcription factors with important roles in hematopoiesis. Expression of HOX genes is also a common feature of acute myeloid leukemia (AML), and functional studies have suggested that HOX-dependent pathways may contribute to leukemogenesis. Although HOX expression is known to correlate with specific AML mutations, the patterns of expression of all 39 HOX genes in primary AML samples, and their relationships with recurrent AML mutations, are incompletely understood. In addition, little is known about the influence of AML mutations on DNA methylation at the HOX loci, and the relationship between HOX gene expression and methylation in AML.In this study, we carried out a combined analysis of gene expression data from microarray and RNA-sequencing platforms and genome-wide DNA array-based methylation from 189 primary AML samples that have been previously characterized by either whole-genome or whole exome sequencing. We also measured expression and methylation using the same platforms from normal bone marrow subsets, including CD34+ cells, promyelocytes, monocytes, neutrophils and lymphocytes, and obtained expression data from CD34+ hematopoietic precursors generated from in vitro differentiation of human embryonic stem cells.Our analysis confirmed previous work on the general patterns of HOX expression in AML. The HOXA and HOXB genes showed variation both within each cluster and across the AMLs, although high level expression was restricted to a subset of these genes, including HOXA3, HOXA5, HOXA7, HOXA9, HOXA10, HOXB2-HOXB4, and HOXB6, as well as HOX cofactor MEIS1; HOXC and HOXD genes were minimally expressed in all of the samples. These observations were orthogonally validated by RNA-seq, and with a targeted Nanostring expression platform. Consistent with previous studies, MLL-positive AML samples (n=11) expressed only HOXA genes and MEIS1. AML samples with CBFB-MYH11 rearrangements (n=12) showed expression of only MEIS1, and HOXB2-HOXB4 at moderate levels; RUNX1-RUNX1T1 (n=7) and PML-RARA (n=19) samples did not detectably express any HOX genes. In AMLs with a normal karyotype (n=85), we observed two distinct patterns; one pattern displayed little or no HOX gene expression (7/85; 8%), and another displayed canonical expression of a specific subset of the HOXA and HOXB genes and MEIS1 (78/85; 92%) with similar relative HOX gene expression levels in all cases. Comparison of this pattern with normal bone marrow revealed the same HOX expression pattern in normal CD34+ cells; additional analysis showed that this pattern was confined to hematopoietic stem/progenitor cells, but was not seen in more mature cells, including other CD34+ subsets, promyelocytes, monocytes and neutrophils. We also measured HOX gene expression in CD34+ hematopoietic precursors generated from in vitro differentiation of human embryonic stem cells, which revealed expression of only MEIS1 and the canonical HOXB genes, suggesting that activation of these genes may represent the earliest events in the HOX pathway of hematopoietic development. Correlation of HOX expression with recurrent AML mutations by gene set enrichment analysis demonstrated a significant association with NPM1 (P<10-4) and DNMT3A (P<10-2) mutations, but not with other recurrent somatic mutations, including FLT3,IDH1/IDH2, and TET2. Methylation at the HOX loci demonstrated patterns that correlated with HOX expression, including hypomethylation at HOX promoters in samples with high expression. However, additional mutation-specific patterns were apparent. For example, NPM1-mutant AMLs demonstrated a distinct methylation pattern that included hypomethylation at the HOXB3 promoter, which was not shared with CBFB-MYH11 cases or other AMLs with HOXB3 expression.In summary, our comprehensive analysis demonstrates canonical expression and methylation patterns at the HOX loci in AML. These patterns correspond to specific recurrent AML mutations, and the dominant pattern in most normal karyotype AMLs mimics the signature of hematopoietic stem cells. This supports previous observations of developmental regulation of HOX genes in hematopoiesis, and implies that this normal stem cell signature is “captured” in the majority of AMLs with normal karyotype. In addition, distinct methylation patterns at HOX loci suggest that multiple regulatory mechanisms are involved in HOX expression in AML. Disclosures:No relevant conflicts of interest to declare.

Vertebrate GAGA factor associated insulator elements demarcate homeotic genes in the HOX clusters

BackgroundHox genes impart segment identity to body structures along the anterior-posterior axis and are crucial for the proper development of all organisms. Multiple regulatory elements, best defined in Drosophila melanogaster, ensure that Hox expression patterns follow the spatial and temporal colinearity reflected in their tight genomic organization. However, the precise mechanisms that regulate colinear patterns of Hox gene expression remain unclear, especially in higher vertebrates where it is not fully determined how the distinct activation domains of the tightly clustered Hox genes are defined independently of each other. Here, we report the identification of a large number of novel cis-elements at mammalian Hox clusters that can help in regulating their precise expression pattern.ResultsWe have identified DNA elements at all four murine Hox clusters that show poor association with histone H3 in chromatin immunoprecipitation (ChIP)-chip tiling arrays. The majority of these elements lie in the intergenic regions segregating adjacent Hox genes; we demonstrate that they possess efficient enhancer-blocking activity in mammalian cells. Further, we find that these histone-free intergenic regions bear GA repeat motifs and associate with the vertebrate homolog of the GAGA binding boundary factor. This suggests that they can act as GAGA factor-dependent chromatin boundaries that create independent domains, insulating each Hox gene from the influence of neighboring regulatory elements.ConclusionsOur results reveal a large number of potential regulatory elements throughout the murine Hox clusters. We further demarcate the precise location of several novel cis-elements bearing chromatin boundary activity that appear to segregate successive Hox genes. This reflects a pattern reminiscent of the organization of homeotic genes in Drosophila, where such regulatory elements have been characterized. Our findings thus provide new insights into the regulatory processes and evolutionarily conserved epigenetic mechanisms that control homeotic gene expression.

Hoxc Gene Collinear Expression and Epigenetic Modifications Established during Embryogenesis Are Maintained until after Birth

The Hox genes, which are organized into clusters on different chromosomes, are key regulators of embryonic anterior-posterior (A-P) body pattern formation and are expressed at specific times and in specific positions in developing vertebrate embryos. Previously, we have shown that histone methylation patterns are closely correlated with collinear Hox gene expression patterns along the A-P axis of E14.5 mouse embryos. Since histone modification is thought to play a crucial mechanistic role in the highly coordinated pattern of collinear Hox gene expression, we examined the maintenance of the spatial collinear expression pattern of Hoxc genes and the corresponding histone modifications during embryogenesis and in early postnatal mice. Hox expression patterns and histone modifications were analyzed by semi-quantitative RT-PCR and chromatin immunoprecipitation (ChIP)-PCR analyses, respectively. The spatiotemporal expression patterns of Hoxc genes in a cluster were maintained until the early postnatal stage (from E8.5 through P5). Examination of histone modifications in E14.5 and P5 tissues revealed that level of H3K27me3 is only a weak correlation with collinear Hoxc gene expression in the trunk regions although diminished in general, however the enrichment of H3K4me3 is strongly correlated with the gene expression in both stages. In summary, the initial spatiotemporal collinear expression pattern of Hoxc genes and epigenetic modifications are maintained after birth, likely contributing to the establishment of the gene expression code for position in the anatomic body axis throughout the entire life of the organism.

Changing topographic Hox expression in blood vessels results in regionally distinct vessel wall remodeling

SummaryThe distinct topographic Hox expression patterns observed in vascular smooth muscle cells (VSMCs) of the adult cardiovascular system suggest that these transcriptional regulators are critical for maintaining region-specific physiological properties of blood vessels. To test this proposition, we expanded the vascular Hoxc11 expression domain normally restricted to the lower limbs by utilizing an innovative integrated tetracycline regulatory system and Transgelin promoter elements to induce Hoxc11 expression universally in VSMCs of transgenic mice. Ectopic Hoxc11 expression in carotid arteries, aortic arch and descending aorta resulted in drastic vessel wall remodeling involving elastic laminae fragmentation, medial smooth muscle cell loss, and intimal lesion formation. None of these alterations were observed upon induction of Hoxc11 transgene expression in the femoral artery, i.e. the natural Hoxc11 activity domain, although this vessel was greatly enlarged, comparable to the topographically restricted vascular changes seen in Hoxc11−/− mice. To begin defining Hoxc11-controlled pathways of vascular remodeling, we performed immunolabeling studies in conjunction with co-transfection and chromatin immunoprecipitation (ChIP) assays using mouse vascular smooth muscle (MOVAS) cells. The results suggest direct transcriptional control of two members of the matrix metalloproteinase (Mmp) family, including Mmp2 and Mmp9 that are known as key players in the inception and progression of vascular remodeling events. In summary, the severe vascular abnormalities resulting from the induced dysregulated expression of a Hox gene with regional vascular patterning functions suggest that proper Hox function and regulation is critical for maintaining vascular functional integrity.

Hox and Pbx Factors Control Retinoic Acid Synthesis during Hindbrain Segmentation

In vertebrate embryos, retinoic acid (RA) synthesized in the mesoderm by Raldh2 emanates to the hindbrain neuroepithelium, where it induces anteroposterior (AP)-restricted Hox expression patterns and rhombomere segmentation. However, how appropriate spatiotemporal RA activity is generated in the hindbrain is poorly understood. By analyzing Pbx1/Pbx2 and Hoxa1/Pbx1 null mice, we found that Raldh2 is itself under the transcriptional control of these factors and that the resulting RA-deficient phenotypes can be partially rescued by exogenous RA. Hoxa1-Pbx1/2-Meis2 directly binds a specific regulatory element that is required to maintain normal Raldh2 expression levels in vivo. Mesoderm-specific Xhoxa1 and Xpbx1b knockdowns in Xenopus embryos also result in Xraldh2 downregulation and hindbrain defects similar to mouse mutants, demonstrating conservation of this Hox-Pbx-dependent regulatory pathway. These findings reveal a feed-forward mechanism linking Hox-Pbx-dependent RA synthesis during early axial patterning with the establishment of spatially restricted Hox-Pbx activity in the developing hindbrain.

XMeis3 Is Necessary for Mesodermal Hox Gene Expression and Function

Hox transcription factors provide positional information during patterning of the anteroposterior axis. Hox transcription factors can co-operatively bind with PBC-class co-factors, enhancing specificity and affinity for their appropriate binding sites. The nuclear localisation of these co-factors is regulated by the Meis-class of homeodomain proteins. During development of the zebrafish hindbrain, Meis3 has previously been shown to synergise with Hoxb1 in the autoregulation of Hoxb1. In Xenopus XMeis3 posteriorises the embryo upon ectopic expression. Recently, an early temporally collinear expression sequence of Hox genes was detected in Xenopus gastrula mesoderm (see intro. P3). There is evidence that this sequence sets up the embryo's later axial Hox expression pattern by time-space translation. We investigated whether XMeis3 is involved in regulation of this early mesodermal Hox gene expression. Here, we present evidence that XMeis3 is necessary for expression of Hoxd1, Hoxb4 and Hoxc6 in mesoderm during gastrulation. In addition, we show that XMeis3 function is necessary for the progression of gastrulation. Finally, we present evidence for synergy between XMeis3 and Hoxd1 in Hoxd1 autoregulation in mesoderm during gastrulation.