Dynamic regulation of Nanog and stem cell-signaling pathways by Hoxa1 during early neuro-ectodermal differentiation of ES cells. Proc Natl Acad Sci U S A. 2017;114:5838-5845. [PMID: 28584089 DOI: 10.1073/pnas.1610612114]

Neural stem/progenitor cells from olfactory neuroepithelium collected by nasal brushing as a cell model reflecting molecular and cellular dysfunctions in schizophrenia

OBJECTIVES

Neural stem/progenitor cells derived from olfactory neuroepithelium (hereafter olfactory neural stem/progenitor cells, ONSPCs) are emerging as a potential tool in the exploration of psychiatric disorders. The present study intended to assess whether ONSPCs could help discern individuals with schizophrenia (SZ) from non-schizophrenic (NS) subjects by exploring specific cellular and molecular features.

METHODS

ONSPCs were collected from 19 in-patients diagnosed with SZ and 31 NS individuals and propagated in basal medium. Mitochondrial ATP production, expression of β-catenin and cell proliferation, which are described to be altered in SZ, were examined in freshly isolated or newly thawed ONSPCs after a few culture passages.

RESULTS

SZ-ONSPCs exhibited a lower mitochondrial ATP production and insensitivity to agents capable of positively or negatively affecting β-catenin expression with respect to NS-ONSPCs. As to proliferation, it declined in SZ-ONSPCs as the number of culture passages increased compared to a steady level of growth shown by NS-ONSPCs.

CONCLUSIONS

The ease and safety of sample collection as well as the differences observed between NS- and SZ-ONSPCs, may lay the groundwork for a new approach to obtain biological material from a large number of living individuals and gain a better understanding of the mechanisms underlying SZ pathophysiology.

Lola-I is a promoter pioneer factor that establishes de novo Pol II pausing during development

While the accessibility of enhancers is dynamically regulated during development, promoters tend to be constitutively accessible and poised for activation by paused Pol II. By studying Lola-I, a Drosophila zinc finger transcription factor, we show here that the promoter state can also be subject to developmental regulation independently of gene activation. Lola-I is ubiquitously expressed at the end of embryogenesis and causes its target promoters to become accessible and acquire paused Pol II throughout the embryo. This promoter transition is required but not sufficient for tissue-specific target gene activation. Lola-I mediates this function by depleting promoter nucleosomes, similar to the action of pioneer factors at enhancers. These results uncover a level of regulation for promoters that is normally found at enhancers and reveal a mechanism for the de novo establishment of paused Pol II at promoters.

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.

Ratiometric sensing of Pnt and Yan transcription factor levels confers ultrasensitivity to photoreceptor fate transitions in Drosophila

Cell state transitions are often triggered by large changes in the concentrations of transcription factors and therefore large differences in their stoichiometric ratios. Whether cells can elicit transitions using modest changes in the ratios of co-expressed factors is unclear. Here, we investigate how cells in the Drosophila eye resolve state transitions by quantifying the expression dynamics of the ETS transcription factors Pnt and Yan. Eye progenitor cells maintain a relatively constant ratio of Pnt/Yan protein, despite expressing both proteins with pulsatile dynamics. A rapid and sustained twofold increase in the Pnt/Yan ratio accompanies transitions to photoreceptor fates. Genetic perturbations that modestly disrupt the Pnt/Yan ratio produce fate transition defects consistent with the hypothesis that transitions are normally driven by a twofold shift in the ratio. A biophysical model based on cooperative Yan-DNA binding coupled with non-cooperative Pnt-DNA binding illustrates how twofold ratio changes could generate ultrasensitive changes in target gene transcription to drive fate transitions. Thus, coupling cell state transitions to the Pnt/Yan ratio sensitizes the system to modest fold-changes, conferring robustness and ultrasensitivity to the developmental program.

A conserved HOTAIRM1-HOXA1 regulatory axis contributes early to neuronal differentiation

HOTAIRM1 is unlike most long non-coding RNAs in that its sequence is highly conserved across mammals. Such evolutionary conservation points to it having a role in key cellular processes. We previously reported that HOTAIRM1 is required to curb premature activation of downstream HOXA genes in a cell model recapitulating their sequential induction during development. We found that it regulates 3' HOXA gene expression by a mechanism involving epigenetic and three-dimensional chromatin changes. Here we show that HOTAIRM1 participates in proper progression through the early stages of neuronal differentiation. We found that it can associate with the HOXA1 transcription factor and contributes to its downstream transcriptional program. Particularly, HOTAIRM1 affects the NANOG/POU5F1/SOX2 core pluripotency network maintaining an undifferentiated cell state. HOXA1 depletion similarly perturbed expression of these pluripotent factors, suggesting that HOTAIRM1 is a modulator of this transcription factor pathway. Also, given that binding of HOTAIRM1 to HOXA1 was observed in different cell types and species, our results point to this ribonucleoprotein complex as an integral part of a conserved HOTAIRM1-HOXA1 regulatory axis modulating the transition from a pluripotent to a differentiated neuronal state.

ATG7-mediated autophagy facilitates embryonic stem cell exit from naive pluripotency and marks commitment to differentiation

Macroautophagy/autophagy is a conserved cellular mechanism to degrade unneeded cytoplasmic proteins and organelles to recycle their components, and it is critical for embryonic stem cell (ESC) self-renewal and somatic cell reprogramming. Whereas autophagy is essential for early development of embryos, no information exists regarding its functions during the transition from naive-to-primed pluripotency. Here, by using an in vitro transition model of ESCs to epiblast-like cells (EpiLCs), we find that dynamic changes in ATG7-dependent autophagy are critical for the naive-to-primed transition, and are also necessary for germline specification. RNA-seq and ATAC-seq profiling reveal that NANOG acts as a barrier to prevent pluripotency transition, and autophagy-dependent NANOG degradation is important for dismantling the naive pluripotency expression program through decommissioning of naive-associated active enhancers. Mechanistically, we found that autophagy receptor protein SQSTM1/p62 translocated into the nucleus during the pluripotency transition period and is preferentially associated with K63 ubiquitinated NANOG for selective protein degradation. In vivo, loss of autophagy by ATG7 depletion disrupts peri-implantation development and causes increased chromatin association of NANOG, which affects neuronal differentiation by competitively binding to OTX2-specific neuroectodermal development-associated regions. Taken together, our findings reveal that autophagy-dependent degradation of NANOG plays a critical role in regulating exit from the naive state and marks distinct cell fate allocation during lineage specification.Abbreviations: 3-MA: 3-methyladenine; EpiLC: epiblast-like cell; ESC: embryonic stem cell; PGC: primordial germ cell.

The Oncogenic and Tumor Suppressive Long Non-Coding RNA–microRNA–Messenger RNA Regulatory Axes Identified by Analyzing Multiple Platform Omics Data from Cr(VI)-Transformed Cells and Their Implications in Lung Cancer

Chronic exposure to hexavalent chromium (Cr(VI)) causes lung cancer in humans, however, the underlying mechanism has not been well understood. Long non-coding RNAs (lncRNAs) and microRNAs (miRNAs) are commonly studied non-coding RNAs. miRNAs function mainly through interaction with the 3′-untranslated regions of messenger RNAs (mRNAs) to down-regulate gene expression. LncRNAs have been shown to function as competing endogenous RNAs (ceRNAs) to sponge miRNAs and regulate gene expression. It is now well accepted that lncRNAs and miRNAs could function as oncogenes or tumor suppressors. Dysregulations of lncRNAs and miRNAs have been shown to play important roles in cancer initiation, progression, and prognosis. To explore the mechanism of Cr(VI) lung carcinogenesis, we performed lncRNA, mRNA, and miRNA microarray analysis using total RNAs from our previously established chronic Cr(VI) exposure malignantly transformed and passage-matched control human bronchial epithelial BEAS-2B cells. Based on the differentially expressed lncRNAs, miRNAs, and mRNAs between the control (BEAS-2B-Control) and Cr(VI)-transformed (BEAS-Cr(VI)) cells and by using the lncRNA–miRNA interaction and miRNA target prediction algorithms, we identified three oncogenic (HOTAIRM1/miR-182-5p/ERO1A, GOLGA8B/miR-30d-5p/RUNX2, and PDCD6IPP2/miR-23a-3p/HOXA1) and three tumor suppressive (ANXA2P1/miR-20b-5p/FAM241A (C4orf32), MIR99AHG/miR-218-5p/GPM6A, and SH3RF3-AS1/miR-34a-5p/HECW2) lncRNA–miRNA–mRNA regulatory axes. Moreover, the relevance of these three oncogenic and three tumor suppressive lncRNA–miRNA–mRNA regulatory axes in lung cancer was explored by analyzing publicly available human lung cancer omics datasets. It was found that the identified three oncogenic lncRNA–miRNA–mRNA regulatory axes (HOTAIRM1/miR-182-5p/ERO1A, GOLGA8B/miR-30d-5p/RUNX2, and PDCD6IPP2/miR-23a-3p/HOXA1) and the three tumor suppressive lncRNA–miRNA–mRNA regulatory axes (ANXA2P1/miR-20b-5p/FAM241A (C4orf32), MIR99AHG/miR-218-5p/GPM6A, and SH3RF3-AS1/miR-34a-5p/HECW2) have significant diagnostic and prognosis prediction values in human lung cancer. In addition, our recent studies showed that Cr(VI)-transformed cells display cancer stem cell (CSC)-like properties. Further bioinformatics analysis identified the oncogenic lncRNA–miRNA–mRNA regulatory axes as the potential regulators of cancer stemness. In summary, our comprehensive analysis of multiple platform omics datasets obtained from Cr(VI)-transformed human bronchial epithelial cells identified several oncogenic and tumor suppressive lncRNA–miRNA–mRNA regulatory axes, which may play important roles in Cr(VI) carcinogenesis and lung cancer in general.

Pluripotency factors regulate the onset of Hox cluster activation in the early embryo

Pluripotent cells are a transient population of the mammalian embryo dependent on transcription factors, such as OCT4 and NANOG, which maintain pluripotency while suppressing lineage specification. However, these factors are also expressed during early phases of differentiation, and their role in the transition from pluripotency to lineage specification is largely unknown. We found that pluripotency factors play a dual role in regulating key lineage specifiers, initially repressing their expression and later being required for their proper activation. We show that Oct4 is necessary for activation of HoxB genes during differentiation of embryonic stem cells and in the embryo. In addition, we show that the HoxB cluster is coordinately regulated by OCT4 binding sites located at the 3' end of the cluster. Our results show that core pluripotency factors are not limited to maintaining the precommitted epiblast but are also necessary for the proper deployment of subsequent developmental programs.

HOXA1, a breast cancer oncogene

More than 25 years ago, the first literature records mentioned HOXA1 expression in human breast cancer. A few years later, HOXA1 was confirmed as a proper oncogene in mammary tissue. In the following two decades, molecular data about the mode of action of the HOXA1 protein, the factors contributing to activate and maintain HOXA1 gene expression and the identity of its target genes have accumulated and provide a wider view on the association of this transcription factor to breast oncogenesis. Large-scale transcriptomic data gathered from wide cohorts of patients further allowed refining the relationship between breast cancer type and HOXA1 expression. Several recent reports have reviewed the connection between cancer hallmarks and the biology of HOX genes in general. Here we take HOXA1 as a paradigm and propose an extensive overview of the molecular data centered on this oncoprotein, from what its expression modulators, to the interactors contributing to its oncogenic activities, and to the pathways and genes it controls. The data converge to an intricate picture that answers questions on the multi-modality of its oncogene activities, point towards better understanding of breast cancer aetiology and thereby provides an appraisal for treatment opportunities.

Transcriptional Regulation and Implications for Controlling Hox Gene Expression

Hox genes play key roles in axial patterning and regulating the regional identity of cells and tissues in a wide variety of animals from invertebrates to vertebrates. Nested domains of Hox expression generate a combinatorial code that provides a molecular framework for specifying the properties of tissues along the A–P axis. Hence, it is important to understand the regulatory mechanisms that coordinately control the precise patterns of the transcription of clustered Hox genes required for their roles in development. New insights are emerging about the dynamics and molecular mechanisms governing transcriptional regulation, and there is interest in understanding how these may play a role in contributing to the regulation of the expression of the clustered Hox genes. In this review, we summarize some of the recent findings, ideas and emerging mechanisms underlying the regulation of transcription in general and consider how they may be relevant to understanding the transcriptional regulation of Hox genes.

Enhancer architecture-dependent multilayered transcriptional regulation orchestrates RA signaling-induced early lineage differentiation of ESCs

Signaling pathway-driven target gene transcription is critical for fate determination of embryonic stem cells (ESCs), but enhancer-dependent transcriptional regulation in these processes remains poorly understood. Here, we report enhancer architecture-dependent multilayered transcriptional regulation at the Halr1–Hoxa1 locus that orchestrates retinoic acid (RA) signaling-induced early lineage differentiation of ESCs. We show that both homeobox A1 (Hoxa1) and Hoxa adjacent long non-coding RNA 1 (Halr1) are identified as direct downstream targets of RA signaling and regulated by RARA/RXRA via RA response elements (RAREs). Chromosome conformation capture-based screens indicate that RA signaling promotes enhancer interactions essential for Hoxa1 and Halr1 expression and mesendoderm differentiation of ESCs. Furthermore, the results also show that HOXA1 promotes expression of Halr1 through binding to enhancer; conversely, loss of Halr1 enhances interaction between Hoxa1 chromatin and four distal enhancers but weakens interaction with chromatin inside the HoxA cluster, leading to RA signaling-induced Hoxa1 overactivation and enhanced endoderm differentiation. These findings reveal complex transcriptional regulation involving synergistic regulation by enhancers, transcription factors and lncRNA. This work provides new insight into intrinsic molecular mechanisms underlying ESC fate determination during RA signaling-induced early differentiation.

Inability to switch from ARID1A-BAF to ARID1B-BAF impairs exit from pluripotency and commitment towards neural crest formation in ARID1B-related neurodevelopmental disorders

Subunit switches in the BAF chromatin remodeler are essential during development. ARID1B and its paralog ARID1A encode for mutually exclusive BAF subunits. De novo ARID1B haploinsufficient mutations cause neurodevelopmental disorders, including Coffin-Siris syndrome, which is characterized by neurological and craniofacial features. Here, we leveraged ARID1B^+/− Coffin-Siris patient-derived iPSCs and modeled cranial neural crest cell (CNCC) formation. We discovered that ARID1B is active only during the first stage of this process, coinciding with neuroectoderm specification, where it is part of a lineage-specific BAF configuration (ARID1B-BAF). ARID1B-BAF regulates exit from pluripotency and lineage commitment by attenuating thousands of enhancers and genes of the NANOG and SOX2 networks. In iPSCs, these enhancers are maintained active by ARID1A-containing BAF. At the onset of differentiation, cells transition from ARID1A- to ARID1B-BAF, eliciting attenuation of the NANOG/SOX2 networks and triggering pluripotency exit. Coffin-Siris patient cells fail to perform the ARID1A/ARID1B switch, and maintain ARID1A-BAF at the pluripotency enhancers throughout all stages of CNCC formation. This leads to persistent NANOG/SOX2 activity which impairs CNCC formation. Despite showing the typical neural crest signature (TFAP2A/SOX9-positive), ARID1B-haploinsufficient CNCCs are also aberrantly NANOG-positive. These findings suggest a connection between ARID1B mutations, neuroectoderm specification and a pathogenic mechanism for Coffin-Siris syndrome.

Mutations in the ARID1B subunit of the BAF chromatin remodeling complex are associated with the neurodevelopmental Coffin-Siris syndrome. Here the authors reveal that there is a transition from ARID1A-containing complexes to ARID1B during cranial neural crest cell differentiation that is impaired in Coffin-Siris patient-derived cells, which is important for exit from pluripotency.

Oleamide, a Sleep-Inducing Supplement, Upregulates Doublecortin in Hippocampal Progenitor Cells via PPARα

BACKGROUND

Doublecortin (DCX), a microtubule associated protein, has emerged as a central biomarker of hippocampal neurogenesis. However, molecular mechanisms by which DCX is regulated are poorly understood.

OBJECTIVE

Since sleep is involved with the acquisition of memory and oleamide or 9-Octadecenamide (OCT) is a sleep-inducing supplement in human, we examined whether OCT could upregulate DCX in hippocampal progenitor cells (HPCs).

METHODS

We employed real-time PCR, western blot, immunostaining, chromatin immunoprecipitation, lentiviral transduction in HPCs, and the calcium influx assay.

RESULTS

OCT directly upregulated the transcription of Dcx in HPCs via activation of peroxisome proliferator-activated receptor α (PPARα), a lipid-lowering transcription factor. We observed that, HPCs of Ppara-null mice displayed significant impairment in DCX expression and neuronal differentiation as compared to that of wild-type mice. Interestingly, treatment with OCT stimulated the differentiation process of HPCs in wild-type, but not Ppara-null mice. Reconstruction of PPARα in mouse Ppara-null HPCs restored the expression of DCX, which was further stimulated with OCT treatment. In contrast, a dominant-negative mutant of PPARα significantly attenuated the stimulatory effect of OCT on DCX expression and suppressed neuronal differentiation of human neural progenitor cells. Furthermore, RNA microarray, STRING, chromatin immunoprecipitation, site-directed mutagenesis, and promoter reporter assay have identified DCX as a new target of PPARα.

CONCLUSION

These results indicate that OCT, a sleep supplement, directly controls the expression of DCX and suggest that OCT may be repurposed for stimulating the hippocampal neurogenesis.

Pharmacological Manipulation of Wnt/β-Catenin Signaling Pathway in Human Neural Precursor Cells Alters Their Differentiation Potential and Neuronal Yield

The canonical Wnt/β-catenin pathway is a master-regulator of cell fate during embryonic and adult neurogenesis and is therefore a major pharmacological target in basic and clinical research. Chemical manipulation of Wnt signaling during in vitro neuronal differentiation of stem cells can alter both the quantity and the quality of the derived neurons. Accordingly, the use of Wnt activators and blockers has become an integral part of differentiation protocols applied to stem cells in recent years. Here, we investigated the effects of the glycogen synthase kinase-3β inhibitor CHIR99021, which upregulates β-catenin agonizing Wnt; and the tankyrase-1/2 inhibitor XAV939, which downregulates β-catenin antagonizing Wnt. Both drugs and their potential neurogenic and anti-neurogenic effects were studied using stable lines human neural precursor cells (hNPCs), derived from embryonic stem cells, which can be induced to generate mature neurons by chemically-defined conditions. We found that Wnt-agonism by CHIR99021 promotes induction of neural differentiation, while also reducing cell proliferation and survival. This effect was not synergistic with those of pro-neural growth factors during long-term neuronal differentiation. Conversely, antagonism of Wnt by XAV939 consistently prevented neuronal progression of hNPCs. We show here how these two drugs can be used to manipulate cell fate and how self-renewing hNPCs can be used as reliable human in vitro drug-screening platforms.

Reporter gene systems for the identification and characterization of cancer stem cells

Cancer stem cells (CSCs) are tumor cells that share functional characteristics with normal and embryonic stem cells. CSCs have increased tumor-initiating capacity and metastatic potential and lower sensitivity to chemo- and radiotherapy, with important roles in tumor progression and the response to therapy. Thus, a current goal of cancer research is to eliminate CSCs, necessitating an adequate phenotypic and functional characterization of CSCs. Strategies have been developed to identify, enrich, and track CSCs, many of which distinguish CSCs by evaluating the expression of surface markers, the initiation of specific signaling pathways, and the activation of master transcription factors that control stemness in normal cells. We review and discuss the use of reporter gene systems for identifying CSCs. Reporters that are under the control of aldehyde dehydrogenase 1A1, CD133, Notch, Nanog homeobox, Sex-determining region Y-box 2, and POU class 5 homeobox can be used to identify CSCs in many tumor types, track cells in real time, and screen for drugs. Thus, reporter gene systems, in combination with in vitro and in vivo functional assays, can assess changes in the CSCs pool. We present relevant examples of these systems in the evaluation of experimental CSCs-targeting therapeutics, demonstrating their value in CSCs research.

The Hox protein conundrum: The "specifics" of DNA binding for Hox proteins and their partners

Homeotic genes (Hox genes) are homeodomain-transcription factors involved in conferring segmental identity along the anterior-posterior body axis. Molecular characterization of HOX protein function raises some interesting questions regarding the source of the binding specificity of the HOX proteins. How do HOX proteins regulate common and unique target specificity across space and time? This review attempts to summarize and interpret findings in this area, largely focused on results from in vitro and in vivo studies in Drosophila and mouse systems. Recent studies related to HOX protein binding specificity compel us to reconsider some of our current models for transcription factor-DNA interactions. It is crucial to study transcription factor binding by incorporating components of more complex, multi-protein interactions in concert with small changes in binding motifs that can significantly impact DNA binding specificity and subsequent alterations in gene expression. To incorporate the multiple elements that can determine HOX protein binding specificity, we propose a more integrative Cooperative Binding model.

Androgen and glucocorticoid receptor direct distinct transcriptional programs by receptor-specific and shared DNA binding sites

The glucocorticoid (GR) and androgen (AR) receptors execute unique functions in vivo, yet have nearly identical DNA binding specificities. To identify mechanisms that facilitate functional diversification among these transcription factor paralogs, we studied them in an equivalent cellular context. Analysis of chromatin and sequence suggest that divergent binding, and corresponding gene regulation, are driven by different abilities of AR and GR to interact with relatively inaccessible chromatin. Divergent genomic binding patterns can also be the result of subtle differences in DNA binding preference between AR and GR. Furthermore, the sequence composition of large regions (>10 kb) surrounding selectively occupied binding sites differs significantly, indicating a role for the sequence environment in guiding AR and GR to distinct binding sites. The comparison of binding sites that are shared shows that the specificity paradox can also be resolved by differences in the events that occur downstream of receptor binding. Specifically, shared binding sites display receptor-specific enhancer activity, cofactor recruitment and changes in histone modifications. Genomic deletion of shared binding sites demonstrates their contribution to directing receptor-specific gene regulation. Together, these data suggest that differences in genomic occupancy as well as divergence in the events that occur downstream of receptor binding direct functional diversification among transcription factor paralogs.

CTCF-binding element regulates ESC differentiation via orchestrating long-range chromatin interaction between enhancers and HoxA

Proper expression of Homeobox A cluster genes (HoxA) is essential for embryonic stem cell (ESC) differentiation and individual development. However, mechanisms controlling precise spatiotemporal expression of HoxA during early ESC differentiation remain poorly understood. Herein, we identified a functional CTCF-binding element (CBE⁺⁴⁷) closest to the 3'-end of HoxA within the same topologically associated domain (TAD) in ESC. CRISPR-Cas9-mediated deletion of CBE⁺⁴⁷ significantly upregulated HoxA expression and enhanced early ESC differentiation induced by retinoic acid (RA) relative to wild-type cells. Mechanistic analysis by chromosome conformation capture assay (Capture-C) revealed that CBE⁺⁴⁷ deletion decreased interactions between adjacent enhancers, enabling formation of a relatively loose enhancer-enhancer interaction complex (EEIC), which overall increased interactions between that EEIC and central regions of HoxA chromatin. These findings indicate that CBE⁺⁴⁷ organizes chromatin interactions between its adjacent enhancers and HoxA. Furthermore, deletion of those adjacent enhancers synergistically inhibited HoxA activation, suggesting that these enhancers serve as an EEIC required for RA-induced HoxA activation. Collectively, these results provide new insight into RA-induced HoxA expression during early ESC differentiation, also highlight precise regulatory roles of the CTCF-binding element in orchestrating high-order chromatin structure.

Genome-Wide Binding Analyses of HOXB1 Revealed a Novel DNA Binding Motif Associated with Gene Repression

Knowledge of the diverse DNA binding specificities of transcription factors is important for understanding their specific regulatory functions in animal development and evolution. We have examined the genome-wide binding properties of the mouse HOXB1 protein in embryonic stem cells differentiated into neural fates. Unexpectedly, only a small number of HOXB1 bound regions (7%) correlate with binding of the known HOX cofactors PBX and MEIS. In contrast, 22% of the HOXB1 binding peaks display co-occupancy with the transcriptional repressor REST. Analyses revealed that co-binding of HOXB1 with PBX correlates with active histone marks and high levels of expression, while co-occupancy with REST correlates with repressive histone marks and repression of the target genes. Analysis of HOXB1 bound regions uncovered enrichment of a novel 15 base pair HOXB1 binding motif HB1RE (HOXB1 response element). In vitro template binding assays showed that HOXB1, PBX1, and MEIS can bind to this motif. In vivo, this motif is sufficient for direct expression of a reporter gene and over-expression of HOXB1 selectively represses this activity. Our analyses suggest that HOXB1 has evolved an association with REST in gene regulation and the novel HB1RE motif contributes to HOXB1 function in part through a repressive role in gene expression.

Multiple roles of HOX proteins in Metastasis: Let me count the ways

Knowledge of the role of HOX proteins in cancer has been steadily accumulating in the last 25 years. They are encoded by 39 HOX genes arranged in 4 distinct clusters, and have unique and redundant function in all types of cancers. Many HOX genes behave as oncogenic transcriptional factors regulating multiple pathways that are critical to malignant progression in a variety of tumors. Some HOX proteins have dual roles that are tumor-site specific, displaying both oncogenic and tumor suppressor function. The focus of this review is on how HOX proteins contribute to growth or suppression of metastasis. The review will cover HOX protein function in the critical aspects of epithelial-mesenchymal transition, in cancer stem cell sustenance and in therapy resistance, manifested as distant metastasis. The emerging role of adiposity in both initiation and progression of metastasis is described. Defining the role of HOX genes in the metastatic process has identified candidates for targeted cancer therapies that may combat the metastatic process. We will discuss potential therapeutic opportunities, particularly in pathways influenced by HOX proteins.

Differential abilities to engage inaccessible chromatin diversify vertebrate Hox binding patterns

Although Hox genes encode for conserved transcription factors (TFs), they are further divided into anterior, central and posterior groups based on their DNA-binding domain similarity. The posterior Hox group expanded in the deuterostome clade and patterns caudal and distal structures. We aimed to address how similar Hox TFs diverge to induce different positional identities. We studied Hox TF DNA-binding and regulatory activity during an in vitro motor neuron differentiation system that recapitulates embryonic development. We found diversity in the genomic binding profiles of different Hox TFs, even among the posterior group paralogs that share similar DNA-binding domains. These differences in genomic binding were explained by differing abilities to bind to previously inaccessible sites. For example, the posterior group HOXC9 had a greater ability to bind occluded sites than the posterior HOXC10, producing different binding patterns and driving differential gene expression programs. From these results, we propose that the differential abilities of posterior Hox TFs to bind to previously inaccessible chromatin drive patterning diversification.This article has an associated 'The people behind the papers' interview.

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.

Galectin-13/placental protein 13: redox-active disulfides as switches for regulating structure, function and cellular distribution

Galectin-13 (Gal-13) plays numerous roles in regulating the relationship between maternal and fetal tissues. Low expression levels or mutations of the lectin can result in pre-eclampsia. The previous crystal structure and gel filtration data show that Gal-13 dimerizes via formation of two disulfide bonds formed by Cys136 and Cys138. In the present study, we mutated them to serine (C136S, C138S and C136S/C138S), crystalized the variants and solved their crystal structures. All variants crystallized as monomers. In the C136S structure, Cys138 formed a disulfide bond with Cys19, indicating that Cys19 is important for regulation of reversible disulfide bond formation in this lectin. Hemagglutination assays demonstrated that all variants are inactive at inducing erythrocyte agglutination, even though gel filtration profiles indicate that C136S and C138S could still form dimers, suggesting that these dimers do not exhibit the same activity as wild-type (WT) Gal-13. In HeLa cells, the three variants were found to be distributed the same as with WT Gal-13. However, a Gal-13 variant (delT221) truncated at T221 could not be transported into the nucleus, possibly explaining why women having this variant get pre-eclampsia. Considering the normally high concentration of glutathione in cells, WT Gal-13 should exist mostly as a monomer in cytoplasm, consistent with the monomeric variant C136S/C138S, which has a similar ability to interact with HOXA1 as WT Gal-13.

lncRNA HotairM1 Depletion Promotes Self-Renewal of Cancer Stem Cells through HOXA1-Nanog Regulation Loop

In cancer cells, a gain of stemness may have profound implications for tumor initiation, aggressiveness, and clinical outcome. However, the molecular mechanisms underlying the self-renewal maintenance of cancer stem-like cells (CSCs) remain elusive. Here, based on analysis of transcriptome sequencing, we identified a long noncoding RNA (lncRNA) named HotairM1, which is weakly expressed in human colorectal carcinoma and uveal melanoma, and a much lower expression in corresponding CSCs. Our results showed that HotairM1 depletion could promote CSC self-renewal and tumor propagation. Mechanistically, HotairM1 recruit EZH2 and SUZ12 to the promoter of its target gene HOXA1, leading to histone H3K27 trimethylation and epigenetic silencing of HOXA1. The silence of HOXA1 subsequently induces the H3K27 acetylation at the enhancer site of Nanog gene to upregulate its expression. The enrichment of Nanog could further inhibit HOXA1 expression, forming a reciprocal regulation loop augmenting the stemness maintaining effect. In summary, our results revealed a lncRNA-based regulatory loop that sustains self-renewal of CSCs, which highlights the critical role of HotairM1 in CSC development through the HOXA1-Nanog signaling loop.

SOX2Mediates Carbon Nanotube-Induced Fibrogenesis and Fibroblast Stem Cell Acquisition

Certain nanosized particles like carbon nanotubes (CNTs) are known to induce pulmonary fibrosis, but the underlying mechanisms are unclear, and efforts to prevent this disease are lacking. Fibroblast-associated stem cells (FSCs) have been suggested as a critical driver of fibrosis induced by CNTs by serving as a renewable source of extracellular matrix-producing cells; however, a detailed understanding of this process remains obscure. Here, we demonstrated that single-walled CNTs induced FSC acquisition and fibrogenic responses in primary human lung fibroblasts. This was indicated by increased expression of stem cell markers (e.g., CD44 and ABCG2) and fibrogenic markers (e.g., collagen and α-SMA) in CNT-exposed cells. These cells also showed increased sphere formation, anoikis resistance, and aldehyde dehydrogenase (ALDH) activities, which are characteristics of stem cells. Mechanistic studies revealed sex-determining region Y-box 2 (SOX2), a self-renewal associated transcription factor, as a key driver of FSC acquisition and fibrogenesis. Upregulation and colocalization of SOX2 and COL1 were found in the fibrotic lung tissues of CNT-exposed mice via oropharyngeal aspiration after 56 days. The knockdown of SOX2 by gene silencing abrogated the fibrogenic and FSC-inducing effects of CNTs. Chromatin immunoprecipitation assays identified SOX2-binding sites on COL1A1 and COL1A2, indicating SOX2 as a transcription factor in collagen synthesis. SOX2 was also found to play a critical role in TGF-β-induced fibrogenesis through its collagen- and FSC-inducing effects. Since many nanomaterials are known to induce TGF-β, our findings that SOX2 regulate FSCs and fibrogenesis may have broad implications on the fibrogenic mechanisms and treatment strategies of various nanomaterial-induced fibrotic disorders.

Accumulation of kynurenine elevates oxidative stress and alters microRNA profile in human bone marrow stromal cells

Kynurenine, a metabolite of tryptophan breakdown, has been shown to increase with age, and plays a vital role in a number of age-related pathophysiological changes, including bone loss. Accumulation of kynurenine in bone marrow stromal cells (BMSCs) has been associated with a decrease in cell proliferation and differentiation, though the exact mechanism by which kynurenine mediates these changes is poorly understood. MiRNAs have been shown to regulate BMSC function, and accumulation of kynurenine may alter the miRNA expression profile of BMSCs. The aim of this study was to identify differentially expressed miRNAs in human BMSCs in response to treatment with kynurenine, and correlate miRNAs function in BMSCs biology through bioinformatics analysis. Human BMSCs were cultured and treated with and without kynurenine, and subsequent miRNA isolation was performed. MiRNA array was performed to identify differentially expressed miRNA. Microarray analysis identified 50 up-regulated, and 36 down-regulated miRNAs in kynurenine-treated BMSC cultures. Differentially expressed miRNA included miR-1281, miR-330-3p, let-7f-5p, and miR-493-5p, which are important for BMSC proliferation and differentiation. KEGG analysis found up-regulated miRNA targeting glutathione metabolism, a pathway critical for removing oxidative species. Our data support that the kynurenine dependent degenerative effect is partially due to changes in the miRNA profile of BMSCs.

A distal enhancer maintaining Hoxa1 expression orchestrates retinoic acid-induced early ESCs differentiation

Retinoic acid (RA) induces rapid differentiation of embryonic stem cells (ESCs), partly by activating expression of the transcription factor Hoxa1, which regulates downstream target genes that promote ESCs differentiation. However, mechanisms of RA-induced Hoxa1 expression and ESCs early differentiation remain largely unknown. Here, we identify a distal enhancer interacting with the Hoxa1 locus through a long-range chromatin loop. Enhancer deletion significantly inhibited expression of RA-induced Hoxa1 and endoderm master control genes such as Gata4 and Gata6. Transcriptome analysis revealed that RA-induced early ESCs differentiation was blocked in Hoxa1 enhancer knockout cells, suggesting a requirement for the enhancer. Restoration of Hoxa1 expression partly rescued expression levels of ∼40% of genes whose expression changed following enhancer deletion, and ∼18% of promoters of those rescued genes were directly bound by Hoxa1. Our data show that a distal enhancer maintains Hoxa1 expression through long-range chromatin loop and that Hoxa1 directly regulates downstream target genes expression and then orchestrates RA-induced early differentiation of ESCs. This discovery reveals mechanisms of a novel enhancer regulating RA-induced Hoxa genes expression and early ESCs differentiation.

Nanog regulates Pou3f1 expression at the exit from pluripotency during gastrulation

Pluripotency is regulated by a network of transcription factors that maintain early embryonic cells in an undifferentiated state while allowing them to proliferate. NANOG is a critical factor for maintaining pluripotency and its role in primordial germ cell differentiation has been well described. However, Nanog is expressed during gastrulation across all the posterior epiblast, and only later in development is its expression restricted to primordial germ cells. In this work, we unveiled a previously unknown mechanism by which Nanog specifically represses genes involved in anterior epiblast lineage. Analysis of transcriptional data from both embryonic stem cells and gastrulating mouse embryos revealed Pou3f1 expression to be negatively correlated with that of Nanog during the early stages of differentiation. We have functionally demonstrated Pou3f1 to be a direct target of NANOG by using a dual transgene system for the controlled expression of Nanog Use of Nanog null ES cells further demonstrated a role for Nanog in repressing a subset of anterior neural genes. Deletion of a NANOG binding site (BS) located nine kilobases downstream of the transcription start site of Pou3f1 revealed this BS to have a specific role in the regionalization of the expression of this gene in the embryo. Our results indicate an active role of Nanog inhibiting neural regulatory networks by repressing Pou3f1 at the onset of gastrulation.This article has an associated First Person interview with the joint first authors of the paper.

Epigenomic signatures underpin the axonal regenerative ability of dorsal root ganglia sensory neurons

Axonal injury results in regenerative success or failure, depending on whether the axon lies in the peripheral or the CNS, respectively. The present study addresses whether epigenetic signatures in dorsal root ganglia discriminate between regenerative and non-regenerative axonal injury. Chromatin immunoprecipitation for the histone 3 (H3) post-translational modifications H3K9ac, H3K27ac and H3K27me3; an assay for transposase-accessible chromatin; and RNA sequencing were performed in dorsal root ganglia after sciatic nerve or dorsal column axotomy. Distinct histone acetylation and chromatin accessibility signatures correlated with gene expression after peripheral, but not central, axonal injury. DNA-footprinting analyses revealed new transcriptional regulators associated with regenerative ability. Machine-learning algorithms inferred the direction of most of the gene expression changes. Neuronal conditional deletion of the chromatin remodeler CCCTC-binding factor impaired nerve regeneration, implicating chromatin organization in the regenerative competence. Altogether, the present study offers the first epigenomic map providing insight into the transcriptional response to injury and the differential regenerative ability of sensory neurons.

The pluripotency transcription factor Nanog represses glutathione reductase gene expression in mouse embryonic stem cells

OBJECTIVE

Redox homeostasis maintenance is essential to bring about cellular functions. Particularly, embryonic stem cells (ESCs) have high fidelity mechanisms for DNA repair, high activity of different antioxidant enzymes and low levels of oxidative stress. Although the expression and activity of antioxidant enzymes are reduced throughout the differentiation, the knowledge about the transcriptional regulation of genes involved in defense against oxidative stress is yet restricted. Since glutathione is a central component of a complex system involved in preserving cellular redox status, we aimed to study whether the expression of the glutathione reductase (Gsr) gene, which encodes an essential enzyme for cellular redox homeostasis, is modulated by the transcription factors critical for self-renewal and pluripotency of ESCs.

RESULTS

We found that Gsr gene is expressed in ESCs during the pluripotent state and it was upregulated when these cells were induced to differentiate, concomitantly with Nanog decreased expression. Moreover, we found an increase in Gsr mRNA levels when Nanog was downregulated by a specific shRNA targeting this transcription factor in ESCs. Our results suggest that Nanog represses Gsr gene expression in ESCs, evidencing a role of this crucial pluripotency transcription factor in preservation of redox homeostasis in stem cells.

Autocrine TGF-β1/miR-200s/miR-221/DNMT3B regulatory loop maintains CAF status to fuel breast cancer cell proliferation

Cancer-associated fibroblasts (CAFs) remain active even in the absence of cancer cells. However, the molecular mechanism underlying the sustained active status of CAFs is largely unrevealed. We found that in CAFs, DNMT3B was not only a target of miR-200b, miR-200c and miR-221, but was able to induce DNA methylation of miR-200s promoters. DNMT3B eventually reached a stably high level by the counteracting effect of decreasing miR-200b/c and increasing miR-221 in normal fibroblasts (NFs) with long-term exogenous TGF-β1 treatment, and DNMT3B further led to a low level of miR-200s which established CAF activation. Meanwhile, miR-200s/miR-221/DNMT3B signaling sustained autocrine TGF-β1 maintaining active CAF status. Destruction of the autocrine TGF-β1/miR-200s/miR-221/DNMT3B signaling led to demethylation of miR-200s promoters and further restored the NF phenotypes. Moreover, we confirmed that TCF12, the target of miR-141, stimulated c-Myc/Cyclin D1 axis in breast cancer cells to promote cancer growth by enhancing CXCL12 of CAFs. The current study reveals that the TGF-β1/miR-200s/miR-221/DNMT3B regulatory loop is responsible for the maintenance of CAFs status and is also necessary for CAF function in promoting malignance of breast cancer, which provides a potential target for CAF-driven therapeutic strategy.

Hox genes: Downstream "effectors" of retinoic acid signaling in vertebrate embryogenesis

One of the major regulatory challenges of animal development is to precisely coordinate in space and time the formation, specification, and patterning of cells that underlie elaboration of the basic body plan. How does the vertebrate plan for the nervous and hematopoietic systems, heart, limbs, digestive, and reproductive organs derive from seemingly similar population of cells? These systems are initially established and patterned along the anteroposterior axis (AP) by opposing signaling gradients that lead to the activation of gene regulatory networks involved in axial specification, including the Hox genes. The retinoid signaling pathway is one of the key signaling gradients coupled to the establishment of axial patterning. The nested domains of Hox gene expression, which provide a combinatorial code for axial patterning, arise in part through a differential response to retinoic acid (RA) diffusing from anabolic centers established within the embryo during development. Hence, Hox genes are important direct effectors of retinoid signaling in embryogenesis. This review focuses on describing current knowledge on the complex mechanisms and regulatory processes, which govern the response of Hox genes to RA in several tissue contexts including the nervous system during vertebrate development.

A Hox-TALE regulatory circuit for neural crest patterning is conserved across vertebrates

In jawed vertebrates (gnathostomes), Hox genes play an important role in patterning head and jaw formation, but mechanisms coupling Hox genes to neural crest (NC) are unknown. Here we use cross-species regulatory comparisons between gnathostomes and lamprey, a jawless extant vertebrate, to investigate conserved ancestral mechanisms regulating Hox2 genes in NC. Gnathostome Hoxa2 and Hoxb2 NC enhancers mediate equivalent NC expression in lamprey and gnathostomes, revealing ancient conservation of Hox upstream regulatory components in NC. In characterizing a lamprey hoxα2 NC/hindbrain enhancer, we identify essential Meis, Pbx, and Hox binding sites that are functionally conserved within Hoxa2/Hoxb2 NC enhancers. This suggests that the lamprey hoxα2 enhancer retains ancestral activity and that Hoxa2/Hoxb2 NC enhancers are ancient paralogues, which diverged in hindbrain and NC activities. This identifies an ancestral mechanism for Hox2 NC regulation involving a Hox-TALE regulatory circuit, potentiated by inputs from Meis and Pbx proteins and Hox auto-/cross-regulatory interactions.

Silver Nanoparticles Compromise Female Embryonic Stem Cell Differentiation through Disturbing X Chromosome Inactivation

The widespread use of silver nanoparticles (AgNPs) has raised substantial health risks to human beings. Despite a wealth of progress on toxicity studies, the understanding of the adverse effects on fetuses, embryos, and early stage cells is still rather limited, particularly under low-dose exposure settings. Moreover, nearly all previous studies ascribed AgNP-induced toxic effects to oxidative stress. Differently, we here unearthed a mechanism, namely, interruption of X chromosome inactivation (XCI) in female mouse embryonic stem cells (mESCs). Albeit with no observable cytotoxicity, significant differentiation retardation was found in female mESCs upon low-dose AgNP exposure. Mechanistic investigations uncovered expedited inactivation for the inactive X chromosome (Xi) and attenuated maintenance of the active X chromosome (Xa) state during mESC differentiation upon the challenge of low-dose AgNPs, indicative of disordered XCI. Thereby, a few X-linked genes (which are closely involved in orchestrating ESC differentiation) were found to be repressed, partially attributable to reinforced enrichment of histone modification ( e. g., histone 3 lysine 27 trimethylation, H3K27me3) on their promoter regions, as the result of disordered XCI. In stark contrast to female mESCs, no impairment of differentiation was observed in male mESCs under low-dose AgNP exposure. All considered, our data unearthed that AgNPs at low concentrations compromised the differentiation program of female mESCs through disturbing XCI. Thus, this work would provide a model for the type of studies necessary to advance the understandings on AgNP-induced developmental toxicity.

A workflow for simplified analysis of ATAC-cap-seq data in R

Background

Assay for Transposase-Accessible Chromatin (ATAC)-cap-seq is a high-throughput sequencing method that combines ATAC-seq with targeted nucleic acid enrichment of precipitated DNA fragments. There are increased analytical difficulties arising from working with a set of regions of interest that may be small in number and biologically dependent. Common statistical pipelines for RNA sequencing might be assumed to apply but can give misleading results on ATAC-cap-seq data. A tool is needed to allow a nonspecialist user to quickly and easily summarize data and apply sensible and effective normalization and analysis.

Results

We developed atacR to allow a user to easily analyze their ATAC enrichment experiment. It provides comprehensive summary functions and diagnostic plots for studying enriched tag abundance. Application of between-sample normalization is made straightforward. Functions for normalizing based on user-defined control regions, whole library size, and regions selected from the least variable regions in a dataset are provided. Three methods for detecting differential abundance of tags from enriched methods are provided, including bootstrap t, Bayes factor, and a wrapped version of the standard exact test in the edgeR package. We compared the precision, recall, and F-score of each detection method on resampled datasets at varying replicate, significance threshold, and genes changed and found that the Bayes factor method had the greatest overall detection power, though edgeR was slightly stronger in simulations with lower numbers of genes changed.

Conclusions

Our package allows a nonspecialist user to easily and effectively apply methods appropriate to the analysis of ATAC-cap-seq in a reproducible manner. The package is implemented in pure R and is fully interoperable with common workflows in Bioconductor.

Coupling the roles of Hox genes to regulatory networks patterning cranial neural crest

The neural crest is a transient population of cells that forms within the developing central nervous system and migrates away to generate a wide range of derivatives throughout the body during vertebrate embryogenesis. These cells are of evolutionary and clinical interest, constituting a key defining trait in the evolution of vertebrates and alterations in their development are implicated in a high proportion of birth defects and craniofacial abnormalities. In the hindbrain and the adjacent cranial neural crest cells (cNCCs), nested domains of Hox gene expression provide a combinatorial'Hox-code' for specifying regional properties in the developing head. Hox genes have been shown to play important roles at multiple stages in cNCC development, including specification, migration, and differentiation. However, relatively little is known about the underlying gene-regulatory mechanisms involved, both upstream and downstream of Hox genes. Furthermore, it is still an open question as to how the genes of the neural crest GRN are linked to Hox-dependent pathways. In this review, we describe Hox gene expression, function and regulation in cNCCs with a view to integrating these genes within the emerging gene regulatory network for cNCC development. We highlight early roles for Hox1 genes in cNCC specification, proposing that this may be achieved, in part, by regulation of the balance between pluripotency and differentiation in precursor cells within the neuro-epithelium. We then describe what is known about the regulation of Hox gene expression in cNCCs and discuss this from the perspective of early vertebrate evolution.

The O-GlcNAc transferase OGT interacts with and post-translationally modifies the transcription factor HOXA1

HOXA1 belongs to the HOX family of transcription factors which are key regulators of animal development. Little is known about the molecular pathways controlling HOXA1. Recent data from our group revealed distinct partner proteins interacting with HOXA1. Among them, OGT is an O-linked N-acetylglucosamine (O-GlcNAc) transferase modifying a variety of proteins involved in different cellular processes including transcription. Here, we confirm OGT as a HOXA1 interactor, we characterise which domains of HOXA1 and OGT are required for the interaction, and we provide evidence that OGT post-translationally modifies HOXA1. Mass spectrometry experiments indeed reveal that HOXA1 can be phosphorylated on the AGGTVGSPQYIHHSY peptide and that upon OGT expression, the phosphate adduct is replaced by an O-GlcNAc group.

PRDM14, a putative histone methyl-transferase, interacts with and decreases the stability and activity of the HOXA1 transcription factor

Understanding how the activity of transcription factors like HOX proteins is regulated remains a widely open question. In a recent screen for proteins interacting with HOXA1, we identified a PRDM protein family member, PRDM14, which is known to be transiently co-expressed with HOXA1 in epiblast cells before their specification towards somatic versus germ cell fate. Here, we confirm PRDM14 is an interactor of HOXA1 and we identify the homeodomain of HOXA1 as well as the PR domain and Zinc fingers of PRDM14 to be required for the interaction. An 11-His repeat of HOXA1 previously highlighted to contribute to HOXA1-mediated protein-protein interactions is also involved. At a functional level, we provide evidence that HOXA1 displays an unexpectedly long half-life and demonstrate that PRDM14 can reduce the stability and affect the transcriptional activity of HOXA1.

HOXA1 and TALE proteins display cross-regulatory interactions and form a combinatorial binding code on HOXA1 targets

Hoxa1 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 Hoxa1 and TALE 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.