MYB: A Key Transcription Factor in the Hematopoietic System Subject to Many Levels of Control.

MYB is a master regulator and pioneer factor highly expressed in hematopoietic progenitor cells (HPCs) where it contributes to the reprogramming processes operating during hematopoietic development. MYB plays a complex role being involved in several lineages of the hematopoietic system. At the molecular level, the MYB gene is subject to intricate regulation at many levels through several enhancer and promoter elements, through transcriptional elongation control, as well as post-transcriptional regulation. The protein is modulated by post-translational modifications (PTMs) such as SUMOylation restricting the expression of its downstream targets. Together with a range of interaction partners, cooperating transcription factors (TFs) and epigenetic regulators, MYB orchestrates a fine-tuned symphony of genes expressed during various stages of haematopoiesis. At the same time, the complex MYB system is vulnerable, being a target for unbalanced control and cancer development.

MYB as a Critical Transcription Factor and Potential Therapeutic Target in AML.

Myb was identified over four decades ago as the transforming component of acute leukemia viruses in chickens. Since then it has become increasingly apparent that dysregulated MYB activity characterizes many blood cancers, including acute myeloid leukemia, and that it represents the most "addictive" oncoprotein in many, if not all, such diseases. As a consequence of this tumor-specific dependency for MYB, it has become a major focus of efforts to develop specific antileukemia drugs. Much attention is being given to ways to interrupt the interaction between MYB and cooperating factors, in particular EP300/KAT3B and CBP/KAT3A. Aside from candidates identified through screening of small molecules, the most exciting prospect for novel drugs seems to be the design of peptide mimetics that interfere directly at the interface between MYB and its cofactors. Such peptides combine a high degree of target specificity with good efficacy including minimal effects on normal hematopoietic cells.

MYB regulates the SUMO protease SENP1 and its novel interaction partner UXT, modulating MYB target genes and the SUMO landscape.

SUMOylation is a post-translational modification frequently found on nuclear proteins, including transcription factors (TFs) and coactivators. By controlling the activity of several TFs, SUMOylation may have far-reaching effects. MYB is an example of a developmental TF subjected to SUMO-mediated regulation, through both SUMO conjugation and SUMO binding. How SUMO affects MYB target genes is unknown. Here, we explored the global effect of reduced SUMOylation of MYB on its downstream gene programs. RNA-Seq in K562 cells after MYB knockdown and rescue with mutants having an altered SUMO status revealed a number of differentially regulated genes and distinct gene ontology term enrichments. Clearly, the SUMO status of MYB both quantitatively and qualitatively affects its regulome. The transcriptome data further revealed that MYB upregulates the SUMO protease SENP1, a key enzyme that removes SUMO conjugation from SUMOylated proteins. Given this role of SENP1 in the MYB regulome, we expanded the analysis, mapped interaction partners of SENP1, and identified UXT as a novel player affecting the SUMO system by acting as a repressor of SENP1. MYB inhibits the expression of UXT suggesting that MYB is able not only to control a specific gene program directly but also indirectly by affecting the SUMO landscape through SENP1 and UXT. These findings suggest an autoactivation loop whereby MYB, through enhancing SENP1 and reducing UXT, is itself being activated by a reduced level of repressive SUMOylation. We propose that overexpressed MYB, seen in multiple cancers, may drive this autoactivation loop and contribute to oncogenic activation of MYB.

Pioneer transcription factors are associated with the modulation of DNA methylation patterns across cancers.

Methylation of cytosines on DNA is a prominent modification associated with gene expression regulation. Aberrant DNA methylation patterns have recurrently been linked to dysregulation of the regulatory program in cancer cells. To shed light on the underlying molecular mechanism driving this process, we hypothesised that aberrant methylation patterns could be controlled by the binding of specific transcription factors (TFs) across cancer types. By combining DNA methylation arrays and gene expression data with TF binding sites (TFBSs), we explored the interplay between TF binding and DNA methylation in 19 cancer types. We performed emQTL (expression-methylation quantitative trait loci) analyses independently in each cancer type and identified 13 TFs whose expression levels are correlated with local DNA methylation patterns around their binding sites in at least 2 cancer types. The 13 TFs are mainly associated with local demethylation and are enriched for pioneer function, suggesting a specific role for these TFs in modulating chromatin structure and transcription in cancer patients. Furthermore, we confirmed that de novo methylation is precluded across cancers at CpGs lying in genomic regions enriched for TF binding signatures associated with SP1, CTCF, NRF1, GABPA, KLF9, and/or YY1. The modulation of DNA methylation associated with TF binding was observed at cis-regulatory regions controlling immune- and cancer-associated pathways, corroborating that the emQTL signals were derived from both cancer and tumor-infiltrating cells. As a case example, we experimentally confirmed that FOXA1 knock-down is associated with higher methylation in regions bound by FOXA1 in breast cancer MCF-7 cells. Finally, we reported physical interactions between FOXA1 with TET1 and TET2 both in an in vitro setup and in vivo at physiological levels in MCF-7 cells, adding further support for FOXA1 attracting TET1 and TET2 to induce local demethylation in cancer cells.

Dissecting the transactivation domain (tAD) of the transcription factor c-Myb to assess recent models of tAD function.

Transcription factors use a DNA-binding domain to localize their action and a transactivation domain (tAD) to stimulate activation of the associated gene. Recent work has renewed interest in how tADs activate genes, which remains poorly understood. Key features in the new models are exposure of short linear motifs (SLMs) and liquid-liquid phase separation (LLPS). Inspired by the new models for tAD function, we decided to revisit the tAD of the haematopoietic transcription factor c-Myb by performing a mutational analysis to see how these new models fit and potentially explain the tAD behaviour of this master regulator. We know that c-Myb has an acidic tAD, which contains a well-characterized SLM in the form of a LxxLL motif. By testing 12 alanine-scanning mutants and three mutants with major reorganization of its tAD in two mammalian reporter systems, we found a pattern of effects very close to what would be expected from the SLM-exposure model, with strong effects exerted by both acidic replacements and SLM mutation. When the same mutants were tested in a yeast system, the pattern of effects was dramatically different, with the SLM mutation exerting no effect, and tAD behaviour was much less affected by small alterations, as would be expected from a LLPS model. These observations are discussed in the light of the two new tAD models, and a two-step hypothesis for transactivation, combining both models, is proposed.

LXRα Regulates ChREBPα Transactivity in a Target Gene-Specific Manner through an Agonist-Modulated LBD-LID Interaction.

The cholesterol-sensing nuclear receptor liver X receptor (LXR) and the glucose-sensing transcription factor carbohydrate responsive element-binding protein (ChREBP) are central players in regulating glucose and lipid metabolism in the liver. More knowledge of their mechanistic interplay is needed to understand their role in pathological conditions like fatty liver disease and insulin resistance. In the current study, LXR and ChREBP co-occupancy was examined by analyzing ChIP-seq datasets from mice livers. LXR and ChREBP interaction was determined by Co-immunoprecipitation (CoIP) and their transactivity was assessed by real-time quantitative polymerase chain reaction (qPCR) of target genes and gene reporter assays. Chromatin binding capacity was determined by ChIP-qPCR assays. Our data show that LXRα and ChREBPα interact physically and show a high co-occupancy at regulatory regions in the mouse genome. LXRα co-activates ChREBPα and regulates ChREBP-specific target genes in vitro and in vivo. This co-activation is dependent on functional recognition elements for ChREBP but not for LXR, indicating that ChREBPα recruits LXRα to chromatin in trans. The two factors interact via their key activation domains; the low glucose inhibitory domain (LID) of ChREBPα and the ligand-binding domain (LBD) of LXRα. While unliganded LXRα co-activates ChREBPα, ligand-bound LXRα surprisingly represses ChREBPα activity on ChREBP-specific target genes. Mechanistically, this is due to a destabilized LXRα:ChREBPα interaction, leading to reduced ChREBP-binding to chromatin and restricted activation of glycolytic and lipogenic target genes. This ligand-driven molecular switch highlights an unappreciated role of LXRα in responding to nutritional cues that was overlooked due to LXR lipogenesis-promoting function.

The SUMO protease SENP1 and the chromatin remodeler CHD3 interact and jointly affect chromatin accessibility and gene expression.

The small ubiquitin-like modifier (SUMO) post-translationally modifies lysine residues of transcription factors and co-regulators and thereby contributes to an important layer of control of the activities of these transcriptional regulators. Likewise, deSUMOylation of these factors by the sentrin-specific proteases (SENPs) also plays a role in gene regulation, but whether SENPs functionally interact with other regulatory factors that control gene expression is unclear. In the present work, we focused on SENP1, specifically, on its role in activation of gene expression investigated through analysis of the SENP1 interactome, which revealed that SENP1 physically interacts with the chromatin remodeler chromodomain helicase DNA-binding protein 3 (CHD3). Using several additional methods, including GST pulldown and co-immunoprecipitation assays, we validated and mapped this interaction, and using CRISPR-Cas9-generated CHD3- and SENP1-KO cells (in the haploid HAP1 cell line), we investigated whether these two proteins are functionally linked in regulating chromatin remodeling and gene expression. Genome-wide ATAC-Seq analysis of the CHD3- and SENP1-KO cells revealed a large degree of overlap in differential chromatin openness between these two mutant cell lines. Moreover, motif analysis and comparison with ChIP-Seq profiles in K562 cells pointed to an association of CHD3 and SENP1 with CCCTC-binding factor (CTCF) and SUMOylated chromatin-associated factors. Lastly, genome-wide RNA-Seq also indicated that these two proteins co-regulate the expression of several genes. We propose that the functional link between chromatin remodeling by CHD3 and deSUMOylation by SENP1 uncovered here provides another level of control of gene expression.

GSuite HyperBrowser: integrative analysis of dataset collections across the genome and epigenome.

Background: Recent large-scale undertakings such as ENCODE and Roadmap Epigenomics have generated experimental data mapped to the human reference genome (as genomic tracks) representing a variety of functional elements across a large number of cell types. Despite the high potential value of these publicly available data for a broad variety of investigations, little attention has been given to the analytical methodology necessary for their widespread utilisation. Findings: We here present a first principled treatment of the analysis of collections of genomic tracks. We have developed novel computational and statistical methodology to permit comparative and confirmatory analyses across multiple and disparate data sources. We delineate a set of generic questions that are useful across a broad range of investigations and discuss the implications of choosing different statistical measures and null models. Examples include contrasting analyses across different tissues or diseases. The methodology has been implemented in a comprehensive open-source software system, the GSuite HyperBrowser. To make the functionality accessible to biologists, and to facilitate reproducible analysis, we have also developed a web-based interface providing an expertly guided and customizable way of utilizing the methodology. With this system, many novel biological questions can flexibly be posed and rapidly answered. Conclusions: Through a combination of streamlined data acquisition, interoperable representation of dataset collections, and customizable statistical analysis with guided setup and interpretation, the GSuite HyperBrowser represents a first comprehensive solution for integrative analysis of track collections across the genome and epigenome. The software is available at: https://hyperbrowser.uio.no.

The adaptor protein ARA55 and the nuclear kinase HIPK1 assist c-Myb in recruiting p300 to chromatin.

LIM-domain proteins, containing multiple cysteine-rich zinc finger-like motifs, have been shown to play diverse roles in several cellular processes. A common theme is that they mediate important protein-protein interactions that are key to their function. Androgen receptor-associated protein 55 (ARA55) belongs to this family of bridging proteins containing four C-terminal LIM domains. It has a dual role with functions both at focal adhesions and in the nucleus, apparently shuttling between the two compartments. In the present work, we have expanded our understanding of its nuclear functions by showing that it interacts with three nuclear regulators not previously linked to ARA55. We first identified ARA55 as a novel interaction partner of the nuclear kinase HIPK1 and found that ARA55, like HIPK1, also interacts with the transcription factor c-Myb. In search of a function for these associations, we observed that the coactivator p300 not only binds to c-Myb, but to ARA55 as well. When combined, c-Myb, p300, HIPK1 and ARA55 caused strong synergistic activation of a chromatinized reporter gene. In parallel, all partners, including p300, were efficiently recruited to chromatin at the c-Myb-bound promoter. Consistent with this cooperation, we found that c-Myb and ARA55 share a common set of target genes in an osteosarcoma cellular context. We propose that ARA55 and HIPK1 assist c-Myb in recruiting the coactivator and acetyltransferase p300 to chromatin.

PIAS1 binds p300 and behaves as a coactivator or corepressor of the transcription factor c-Myb dependent on SUMO-status.

The PIAS proteins (Protein Inhibitor of Activated STATs) constitute a family of multifunctional nuclear proteins operating as SUMO E3 ligases and being involved in a multitude of interactions. They participate in a range of biological processes, also beyond their well-established role in the immune system and cytokine signalling. They act both as transcriptional corepressors and coactivators depending on the context. In the present work, we investigated mechanisms by which PIAS1 causes activation or repression of c-Myb dependent target genes. Analysis of global expression data shows that c-Myb and PIAS1 knockdowns affect a subset of common targets, but with a dual outcome consistent with a role of PIAS1 as either a corepressor or coactivator. Our mechanistic studies show that PIAS1 engages in a novel interaction with the acetyltransferase and coactivator p300. Interaction and ChIP analysis suggest a bridging function where PIAS1 enhances p300 recruitment to c-Myb-bound sites through interaction with both proteins. In addition, the E3 activity of PIAS1 enhances further its coactivation. Remarkably, the SUMO status of c-Myb had a decisive role, indicating a SUMO-dependent switch in the way PIAS1 affects c-Myb, either as a coactivator or corepressor. Removal of the two major SUMO-conjugation sites in c-Myb (2KR mutant), which enhances its activity significantly, turned PIAS1 into a corepressor. Also, p300 was less efficiently recruited to chromatin by c-Myb-2KR. We propose that PIAS1 acts as a "protein inhibitor of activated c-Myb" in the absence of SUMOylation while, in its presence, PIAS behaves as a "protein activator of repressed c-Myb".

In the loop: promoter-enhancer interactions and bioinformatics.

Enhancer-promoter regulation is a fundamental mechanism underlying differential transcriptional regulation. Spatial chromatin organization brings remote enhancers in contact with target promoters in cis to regulate gene expression. There is considerable evidence for promoter-enhancer interactions (PEIs). In the recent years, genome-wide analyses have identified signatures and mapped novel enhancers; however, being able to precisely identify their target gene(s) requires massive biological and bioinformatics efforts. In this review, we give a short overview of the chromatin landscape and transcriptional regulation. We discuss some key concepts and problems related to chromatin interaction detection technologies, and emerging knowledge from genome-wide chromatin interaction data sets. Then, we critically review different types of bioinformatics analysis methods and tools related to representation and visualization of PEI data, raw data processing and PEI prediction. Lastly, we provide specific examples of how PEIs have been used to elucidate a functional role of non-coding single-nucleotide polymorphisms. The topic is at the forefront of epigenetic research, and by highlighting some future bioinformatics challenges in the field, this review provides a comprehensive background for future PEI studies.

c-Myb Binding Sites in Haematopoietic Chromatin Landscapes.

Strict control of tissue-specific gene expression plays a pivotal role during lineage commitment. The transcription factor c-Myb has an essential role in adult haematopoiesis and functions as an oncogene when rearranged in human cancers. Here we have exploited digital genomic footprinting analysis to obtain a global picture of c-Myb occupancy in the genome of six different haematopoietic cell-types. We have biologically validated several c-Myb footprints using c-Myb knockdown data, reporter assays and DamID analysis. We show that our predicted conserved c-Myb footprints are highly dependent on the haematopoietic cell type, but that there is a group of gene targets common to all cell-types analysed. Furthermore, we find that c-Myb footprints co-localise with active histone mark H3K4me3 and are significantly enriched at exons. We analysed co-localisation of c-Myb footprints with 104 chromatin regulatory factors in K562 cells, and identified nine proteins that are enriched together with c-Myb footprints on genes positively regulated by c-Myb and one protein enriched on negatively regulated genes. Our data suggest that c-Myb is a transcription factor with multifaceted target regulation depending on cell type.

Identification of c-Myb Target Genes in K562 Cells Reveals a Role for c-Myb as a Master Regulator.

The c-Myb transcription factor is an important regulator of hematopoietic cell development. c-Myb is expressed in immature hematopoietic cells and plays a direct role in lineage fate selection, cell cycle progression, and differentiation of myeloid as well as B- and T-lymphoid progenitor cells. As a DNA-binding transcription factor, c-Myb regulates specific gene programs through activation of target genes. Still, our understanding of these programs is incomplete. Here, we report a set of novel c-Myb target genes, identified using a combined approach: specific c-Myb knockdown by 2 different siRNAs and subsequent global expression profiling, combined with the confirmation of direct binding of c-Myb to the target promoters by ChIP assays. The combination of these 2 approaches, as well as additional validation such as cloning and testing the promoters in reporter assays, confirmed that MYADM, LMO2, GATA2, STAT5A, and IKZF1 are target genes of c-Myb. Additional studies, using chromosome conformation capture, demonstrated that c-Myb target genes may directly interact with each other, indicating that these genes may be coordinately regulated. Of the 5 novel target genes identified, 3 are transcription factors, and one is a transcriptional co-regulator, supporting a role of c-Myb as a master regulator controlling the expression of other transcriptional regulators in the hematopoietic system.

A SUMO-regulated activation function controls synergy of c-Myb through a repressor-activator switch leading to differential p300 recruitment.

Synergy between transcription factors operating together on complex promoters is a key aspect of gene activation. The ability of specific factors to synergize is restricted by sumoylation (synergy control, SC). Focusing on the haematopoietic transcription factor c-Myb, we found evidence for a strong SC linked to SUMO-conjugation in its negative regulatory domain (NRD), while AMV v-Myb has escaped this control. Mechanistic studies revealed a SUMO-dependent switch in the function of NRD. When NRD is sumoylated, the activity of c-Myb is reduced. When sumoylation is abolished, NRD switches into being activating, providing the factor with a second activation function (AF). Thus, c-Myb harbours two AFs, one that is constitutively active and one in the NRD being SUMO-regulated (SRAF). This double AF augments c-Myb synergy at compound natural promoters. A similar SUMO-dependent switch was observed in the regulatory domains of Sp3 and p53. We show that the change in synergy behaviour correlates with a SUMO-dependent differential recruitment of p300 and a corresponding local change in histone H3 and H4 acetylation. We therefore propose a general model for SUMO-mediated SC, where SUMO controls synergy by determining the number and strength of AFs associated with a promoter leading to differential chromatin signatures.

HIPK1 interacts with c-Myb and modulates its activity through phosphorylation.

The transcription factor v-Myb is a potent inducer of myeloid leukaemias, and its cellular homologue c-Myb plays a crucial role in the regulation of haematopoiesis. In a yeast two-hybrid (Y2H) screening we identified the nuclear kinase HIPK1 as an interaction partner for human c-Myb. The interaction involves a C-terminal region of HIPK1, while a bipartite interaction surface was identified in c-Myb, including regions in its N-terminal DNA-binding domain as well as in its C-terminal region. HIPK1 and c-Myb co-localize in distinct nuclear foci upon co-transfection. c-Myb appears to be phosphorylated by HIPK1 in its negative regulatory domain as supported by both in vivo and in vitro data. A functional assay revealed that HIPK1 repressed the ability of c-Myb to activate a chromatin embedded target gene, mim-1, in haematopoetic cells. Our findings point to a novel link between an important kinase and a key regulator of haematopoiesis.

SUMO modification regulates the transcriptional activity of FLASH.

FLASH is a huge multifunctional nuclear protein that has been linked to apoptotic signalling, transcriptional control and Cajal body function. To gain further insight into the functions of the FLASH protein, we performed a yeast two-hybrid screening with FLASH as bait and identified the SUMO-conjugating enzyme Ubc9 as an interaction partner. The main interaction surface for Ubc9 was found in the C-terminal part of FLASH, which is also a target for sumoylation. We identified K1813 as the major sumoylation site in FLASH, being enhanced by the SUMO E3 ligases Pc2 and PIASy. Disruption of this SUMO-conjugation site did not change the speckled subnuclear localization of FLASH, but it caused a reduction in FLASH activity as measured in a Gal4-tethering assay. Interestingly, the SUMO-specific protease SENP1 activated FLASH in the same assay. Overall, our results point to a complex involvement of sumoylation in modulating the function of FLASH.

Revisiting a selection of target genes for the hematopoietic transcription factor c-Myb using chromatin immunoprecipitation and c-Myb knockdown.

The transcription factor c-Myb is an important regulator of hematopoiesis required for proper development of most blood cell lineages in vertebrates. An increasing number of target genes for c-Myb are being published, although with little or no overlap between the lists of genes reported. This raises the question of which criteria a bona fide c-Myb-target gene should satisfy. In the present paper, we have analyzed a set of previously reported target genes using chromatin immunoprecipitation (ChIP) and siRNA-mediated knockdown. Among the seven well-studied c-Myb target genes that we analyzed by ChIP, only ADA, c-MYC and MAT2A seemed to be occupied by c-Myb under our experimental settings in the Myb-positive cell lines Jurkat and HL60. After siRNA-mediated knockdown of c-Myb expression, the expression levels of two out of three ChIP positive Myb target genes, ADA and c-MYC, were strongly affected. These results clearly demonstrate the importance of combining different methods for target gene validation and suggest that a combination of ChIP and c-Myb knockdown may represent a powerful approach to identify a core collection of c-Myb target genes.

The chromatin remodeling factor Mi-2alpha acts as a novel co-activator for human c-Myb.

The c-Myb protein belongs to a group of early hematopoietic transcription factors that are important for progenitor generation and proliferation. These factors have been hypothesized to participate in establishing chromatin patterns specific for hematopoietic genes. In a two-hybrid screening we identified the chromatin remodeling factor Mi-2alpha as an interaction partner for human c-Myb. The main interacting domains were mapped to the N-terminal region of Mi-2alpha and the DNA-binding domain of c-Myb. Surprisingly, functional analysis revealed that Mi-2alpha, previously studied as a subunit in the NuRD co-repressor complex, enhanced c-Myb-dependent reporter activation. Consistently, knock-down of endogenous Mi-2alpha in c-Myb-expressing K562 cells, led to down-regulation of the c-Myb target genes NMU and ADA. When wild-type and helicase-dead Mi-2alpha were compared, the Myb-Mi-2alpha co-activation appeared to be independent of the ATPase/DNA helicase activity of Mi-2alpha. The rationale for the unexpected co-activator function seems to lie in a dual function of Mi-2alpha, by which this factor is able to repress transcription in a helicase-dependent and activate in a helicase-independent fashion, as revealed by Gal4-tethering experiments. Interestingly, desumoylation of c-Myb potentiated the Myb-Mi-2alpha transactivational co-operation, as did co-transfection with p300.

c-Myb associates with PML in nuclear bodies in hematopoietic cells.

The c-Myb transcription factor plays a central role in the regulation of cell growth and differentiation of hematopoietic cells. Being the product of a proto-oncogene, one would expect c-Myb function to be modulated by signal transduction pathways, but our knowledge on such regulation of c-Myb is rather limited. Recently, we and others showed that c-Myb is subjected to sumoylation and that this posttranslational modification has considerable effect on c-Myb's activity. Interestingly, many proteins subjected to SUMO-1 conjugation associate with the promyelocytic leukemia (PML) protein and localize to PML nuclear bodies (PML NBs). Although the precise molecular function of PML NBs still remains to be defined, they seem to play a role in regulation of gene expression and are linked to specific cellular signaling. We show here that c-Myb localizes to PML NBs and that c-Myb interacts with PML as judged by immunofluorescence microcopy and co-immunoprecipitation experiments. Enforced expression of PML IV was shown to enhance c-Myb-dependent reporter activation. Our results imply a role for PML and possibly other components of PML NBs in regulating c-Myb's activity. This novel link between c-Myb and PML, two gene products being implicated in leukemic disorders, suggests that previously unknown mechanisms for regulating c-Myb's activity involving PML may exist.