No abstract

The clustered Hox genes, which are highly conserved across metazoans, encode homeodomain-containing transcription factors that provide a blueprint for segmental identity along the body axis. Recent studies have underscored that in addition to encoding Hox genes, the homeotic clusters contain key noncoding RNA genes that play a central role in development. In this study, we have taken advantage of genome-wide approaches to provide a detailed analysis of retinoic acid (RA)-induced transcriptional and epigenetic changes within the homeotic clusters of mouse embryonic stem cells. Although there is a general colinear response, our analyses suggest a lack of strict colinearity for several genes in the HoxA and HoxB clusters. We have identified transcribed novel noncoding RNAs (ncRNAs) and their cis-regulatory elements that function in response to RA and demonstrated that the expression of these ncRNAs from both strands represent some of the most rapidly induced transcripts in ES cells. Finally, we have provided dynamic analyses of chromatin modifications for the coding and noncoding genes expressed upon activation and suggest that active transcription can occur in the presence of chromatin modifications and machineries associated with repressed transcription state over the clusters. Overall, our data provide a resource for a better understanding of the dynamic nature of the coding and noncoding transcripts and their associated chromatin marks in the regulation of homeotic gene transcription during development.

Homeobox a1 (Hoxa1) is one of the most rapidly induced genes in ES cell differentiation and it is the earliest expressed Hox gene in the mouse embryo. In this study, we used genomic approaches to identify Hoxa1-bound regions during early stages of ES cell differentiation into the neuro-ectoderm. Within 2 h of retinoic acid treatment, Hoxa1 is rapidly recruited to target sites that are associated with genes involved in regulation of pluripotency, and these genes display early changes in expression. The pattern of occupancy of Hoxa1 is dynamic and changes over time. At 12 h of differentiation, many sites bound at 2 h are lost and a new cohort of bound regions appears. At both time points the genome-wide mapping reveals that there is significant co-occupancy of Nanog (Nanog homeobox) and Hoxa1 on many common target sites, and these are linked to genes in the pluripotential regulatory network. In addition to shared target genes, Hoxa1 binds to regulatory regions of Nanog, and conversely Nanog binds to a 3′ enhancer of Hoxa1. This finding provides evidence for direct cross-regulatory feedback between Hoxa1 and Nanog through a mechanism of mutual repression. Hoxa1 also binds to regulatory regions of Sox2 (sex-determining region Y box 2), Esrrb (estrogen-related receptor beta), and Myc, which underscores its key input into core components of the pluripotential regulatory network. We propose a model whereby direct inputs of Nanog and Hoxa1 on shared targets and mutual repression between Hoxa1 and the core pluripotency network provides a molecular mechanism that modulates the fine balance between the alternate states of pluripotency and differentiation.

Members of the Spalt gene family encode putative transcription factors characterized by seven to nine C2H2 zinc finger motifs. Four genes have been identified in mice-Spalt1 to Spalt4 (Sall1 to Sall4). Spalt homologues are widely expressed in neural and mesodermal tissues during early embryogenesis. Sall3 is normally expressed in mice from embryonic day 7 (E7) in the neural ectoderm and primitive streak and subsequently in the brain, peripheral nerves, spinal cord, limb buds, palate, heart, and otic vesicles. We have generated a targeted disruption of Sall3 in mice. Homozygous mutant animals die on the first postnatal day and fail to feed. Examination of the oral structures of these animals revealed that abnormalities were present in the palate and epiglottis from E16.5. In E10.5 embryos, deficiencies in cranial nerves that normally innervate oral structures, particularly the glossopharyngeal nerve (IX), were observed. These studies indicate that Sall3 is required for the development of nerves that are derived from the hindbrain and for the formation of adjacent branchial arch derivatives.

We have identified a new member of the spalt-like gene family in mice, msal-3. We compared the expression patterns of msal-3 and msal-1 during development and show that they exhibit overlapping yet exclusive patterns of expression in the developing forebrain, diencephalon, midbrain/hindbrain boundary and spinal cord. Both genes are expressed from E7 in opposite gradients in primitive streak mesoderm. Subsequently their transcripts are localized to regions of mesenchyme/epithelial interaction in the palate, heart, limbs, anal and urogenital region.

has diverse functional roles in differentiation and development. We identify and characterize properties of regions bound by HOXA1 on a genome-wide basis in differentiating mouse ES cells. HOXA1-bound regions are enriched for clusters of consensus binding motifs for HOX, PBX, and MEIS, and many display co-occupancy of PBX and MEIS. PBX and MEIS are members of the TALE family and genome-wide analysis of multiple TALE members (PBX, MEIS, TGIF, PREP1, and PREP2) shows that nearly all HOXA1 targets display occupancy of one or more TALE members. The combinatorial binding patterns of TALE proteins define distinct classes of HOXA1 targets, which may create functional diversity. Transgenic reporter assays in zebrafish confirm enhancer activities for many HOXA1-bound regions and the importance of HOX-PBX and TGIF motifs for their regulation. Proteomic analyses show that HOXA1 physically interacts on chromatin with PBX, MEIS, and PREP family members, but not with TGIF, suggesting that TGIF may have an independent input into HOXA1-bound regions. Therefore, TALE proteins appear to represent a wide repertoire of HOX cofactors, which may coregulate enhancers through distinct mechanisms. We also discover extensive auto- and cross-regulatory interactions among the and genes, indicating that the specificity of HOXA1 during development may be regulated though a complex cross-regulatory network of HOXA1 and TALE proteins. This study provides new insight into a regulatory network involving combinatorial interactions between HOXA1 and TALE 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.

Hoxa1 has important functional roles in neural crest specification, hindbrain patterning and heart and ear development, yet the enhancers and genes that are targeted by Hoxa1 are largely unknown. In this study, we performed a comprehensive analysis of Hoxa1 target genes using genome-wide Hoxa1 binding data in mouse ES cells differentiated with retinoic acid (RA) into neural fates in combination with differential gene expression analysis in Hoxa1 gain- and loss-of-function mouse and zebrafish embryos. Our analyses reveal that Hoxa1-bound regions show epigenetic marks of enhancers, occupancy of Hox cofactors and differential expression of nearby genes, suggesting that these regions are enriched for enhancers. In support of this, 80 of them mapped to regions with known reporter activity in transgenic mouse embryos based on the Vista enhancer database. Two additional enhancers in Dok5 and Wls1 were shown to mediate neural expression in developing mouse and zebrafish. Overall, our analysis of the putative target genes indicate that Hoxa1 has input to components of major signaling pathways, including Wnt, TGF-β, Hedgehog and Hippo, and frequently does so by targeting multiple components of a pathway such as secreted inhibitors, ligands, receptors and down-stream components. We also identified genes implicated in heart and ear development, neural crest migration and neuronal patterning and differentiation, which may underlie major Hoxa1 mutant phenotypes. Finally, we found evidence for a high degree of evolutionary conservation of many binding regions and downstream targets of Hoxa1 between mouse and zebrafish. Our genome-wide analyses in ES cells suggests that we have enriched for in vivo relevant target genes and pathways associated with functional roles of Hoxa1 in mouse development.