c-Myb Binding Sites in Haematopoietic Chromatin Landscapes

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.

Transcription factor c-Myb: novel prognostic factor in osteosarcoma

The transcription factor c-Myb is an oncoprotein promoting cell proliferation and survival when aberrantly activated/expressed, thus contributing to malignant transformation. Overexpression of c-Myb has been found in leukemias, breast, colon and adenoid cystic carcinoma. Recent studies revealed its expression also in osteosarcoma cell lines and suggested its functional importance during bone development. However, the relevance of c-Myb in control of osteosarcoma progression remains unknown. A retrospective clinical study was carried out to assess a relationship between c-Myb expression in archival osteosarcoma tissues and prognosis in a cohort of high-grade osteosarcoma patients. In addition, MYB was depleted in metastatic osteosarcoma cell lines SAOS-2 LM5 and 143B and their growth, chemosensitivity, migration and metastatic activity were determined. Immunohistochemical analysis revealed that high c-Myb expression was significantly associated with poor overall survival in the cohort and metastatic progression in young patients. Increased level of c-Myb was detected in metastatic osteosarcoma cell lines and its depletion suppressed their growth, colony-forming capacity, migration and chemoresistance in vitro in a cell line-dependent manner. MYB knock-out resulted in reduced metastatic activity of both SAOS-2 LM5 and 143B cell lines in immunodeficient mice. Transcriptomic analysis revealed the c-Myb-driven functional programs enriched for genes involved in the regulation of cell growth, stress response, cell adhesion and cell differentiation/morphogenesis. Wnt signaling pathway was identified as c-Myb target in osteosarcoma cells. Taken together, we identified c-Myb as a negative prognostic factor in osteosarcoma and showed its involvement in the regulation of osteosarcoma cell growth, chemosensitivity, migration and metastatic activity.

A Dual-Site Inhibitor of CBP/p300 KIX is a Selective and Effective Modulator of Myb

The protein-protein interaction between the KIX motif of the transcriptional coactivator CBP/p300 and the transcriptional activator Myb is a high-value target due to its established role in certain acute myeloid leukemias (AML) and potential contributions to other cancers. However, the CBP/p300 KIX domain has multiple binding sites, several structural homologues, many binding partners, and substantial conformational plasticity, making it challenging to specifically target using small-molecule inhibitors. Here, we report a picomolar dual-site inhibitor (MybLL-tide) of the Myb-CBP/p300 KIX interaction. MybLL-tide has higher affinity for CBP/p300 KIX than any previously reported compounds while also possessing 5600-fold selectivity for the CBP/p300 KIX domain over other coactivator domains. MybLL-tide blocks the association of CBP and p300 with Myb in the context of the proteome, leading to inhibition of key Myb·KIX-dependent genes in AML cells. These results show that MybLL-tide is an effective, modifiable tool to selectively target the KIX domain and assess transcriptional effects in AML cells and potentially other cancers featuring aberrant Myb behavior. Additionally, the dual-site design has applicability to the other challenging coactivators that bear multiple binding surfaces.

Chromatin occupancy and target genes of the haematopoietic master transcription factor MYB

The transcription factor MYB is a master regulator in haematopoietic progenitor cells and a pioneer factor affecting differentiation and proliferation of these cells. Leukaemic transformation may be promoted by high MYB levels. Despite much accumulated molecular knowledge of MYB, we still lack a comprehensive understanding of its target genes and its chromatin action. In the present work, we performed a ChIP-seq analysis of MYB in K562 cells accompanied by detailed bioinformatics analyses. We found that MYB occupies both promoters and enhancers. Five clusters (C1-C5) were found when we classified MYB peaks according to epigenetic profiles. C1 was enriched for promoters and C2 dominated by enhancers. C2-linked genes were connected to hematopoietic specific functions and had GATA factor motifs as second in frequency. C1 had in addition to MYB-motifs a significant frequency of ETS-related motifs. Combining ChIP-seq data with RNA-seq data allowed us to identify direct MYB target genes. We also compared ChIP-seq data with digital genomic footprinting. MYB is occupying nearly a third of the super-enhancers in K562. Finally, we concluded that MYB cooperates with a subset of the other highly expressed TFs in this cell line, as expected for a master regulator.

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.

Diamond Blackfan anemia is mediated by hyperactive Nemo-like kinase

Diamond Blackfan Anemia (DBA) is a congenital bone marrow failure syndrome associated with ribosomal gene mutations that lead to ribosomal insufficiency. DBA is characterized by anemia, congenital anomalies, and cancer predisposition. Treatment for DBA is associated with significant morbidity. Here, we report the identification of Nemo-like kinase (NLK) as a potential target for DBA therapy. To identify new DBA targets, we screen for small molecules that increase erythroid expansion in mouse models of DBA. This screen identified a compound that inhibits NLK. Chemical and genetic inhibition of NLK increases erythroid expansion in mouse and human progenitors, including bone marrow cells from DBA patients. In DBA models and patient samples, aberrant NLK activation is initiated at the Megakaryocyte/Erythroid Progenitor (MEP) stage of differentiation and is not observed in non-erythroid hematopoietic lineages or healthy erythroblasts. We propose that NLK mediates aberrant erythropoiesis in DBA and is a potential target for therapy.

Diamond Blackfan Anemia (DBA) is a congenital bone marrow failure syndrome that is associated with anemia. Here, the authors examine the role of Nemo-like kinase (NLK) in erythroid cells in the pathogenesis of DBA and as a potential target for therapy.

c-Myb regulates tumorigenic potential of embryonal rhabdomyosarcoma cells

Rhabdomyosarcomas (RMS) are a heterogeneous group of mesodermal tumors, the most common sub-types are embryonal (eRMS) and alveolar (aRMS) rhabdomyosarcoma. Immunohistochemical analysis revealed c-Myb expression in both eRMS and aRMS. c-Myb has been reported to be often associated with malignant human cancers. We therefore investigated the c-Myb role in RMS using cellular models of RMS. Specific suppression of c-Myb by a lentiviral vector expressing doxycycline (Dox)-inducible c-Myb shRNA inhibited proliferation, colony formation, and migration of the eRMS cell line (RD), but not of the aRMS cell line (RH30). Upon c-Myb knockdown in eRMS cells, cells accumulated in G0/G1 phase, the invasive behaviour of cells was repressed, and elevated levels of myosin heavy chain, marker of muscle differentiation, was detected. Next, we used an RD-based xenograft model to investigate the role of c-Myb in eRMS tumorigenesis in vivo. We found that Dox administration did not result in efficient suppression of c-Myb in growing tumors. However, when c-Myb-deficient RD cells were implanted into SCID mice, we observed inefficient tumor grafting and attenuation of tumor growth during the initial stages of tumor expansion. The presented study suggests that c-Myb could be a therapeutic target in embryonal rhabdomyosarcoma assuming that its expression is ablated.

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.

MYB - A regulatory factor in hematopoiesis

MYB is a transcription factor which was identified in birds as a viral oncogene (v-MYB). Its cellular counterpart was subsequently isolated as c-MYB which has three functional domains - DNA binding domain, transactivation domain and negative regulatory domain. c-MYB is essential for survival, and deletion of both alleles of the gene results in embryonic death. It is highly expressed in hematopoietic cells, thymus and neural tissue, and required for T and B lymphocyte development and erythroid maturation. Additionally, aberrant MYB expression has been found in numerous solid cancer cells and human leukemia. Recent studies have also implicated c-MYB in the regulation of expression of fetal hemoglobin which is highly beneficial to the β-hemoglobinopathies (beta thalassemia and sickle cell disease). These findings suggest that MYB could be a potential therapeutic target in leukemia, and possibly also a target for therapeutic increase of fetal hemoglobin in the β-hemoglobinopathies.

MYB induces the expression of the oncogenic corepressor SKI in acute myeloid leukemia

Acute myeloid leukemia (AML) arises through clonal expansion of transformed myeloid progenitor cells. The SKI proto-oncogene is highly upregulated in different solid tumors and leukemic cells, but little is known about its transcriptional regulation during leukemogenesis. MYB is an important hematopoietic transcription factor involved in proliferation as well as differentiation and upregulated in most human acute leukemias. Here, we find that MYB protein binds within the regulatory region of the SKI gene in AML cells. Reporter gene assays using MYB binding sites present in the SKI gene locus show MYB-dependent transcriptional activation. SiRNA-mediated depletion of MYB in leukemic cell lines reveals that MYB is crucial for SKI gene expression. Consistently, we observed a positive correlation of MYB and SKI expression in leukemic cell lines and in samples of AML patients. Moreover, MYB and SKI both were downregulated by treatment with histone deacetylase inhibitors. Strikingly, differentiation of AML cells induced by depletion of MYB is attenuated by overexpression of SKI. Our findings identify SKI as a novel MYB target gene, relevant for the MYB-induced differentiation block in leukemic cells.

Transcription factor c-Myb inhibits breast cancer lung metastasis by suppression of tumor cell seeding

Metastasis accounts for most of cancer-related deaths. Paracrine signaling between tumor cells and the stroma induces changes in the tumor microenvironment required for metastasis. Transcription factor c-Myb was associated with breast cancer (BC) progression but its role in metastasis remains unclear. Here we show that increased c-Myb expression in BC cells inhibits spontaneous lung metastasis through impaired tumor cell extravasation. On contrary, BC cells with increased lung metastatic capacity exhibited low c-Myb levels. We identified a specific inflammatory signature, including Ccl2 chemokine, that was expressed in lung metastatic cells but was suppressed in tumor cells with higher c-Myb levels. Tumor cell-derived Ccl2 expression facilitated lung metastasis and rescued trans-endothelial migration of c-Myb overexpressing cells. Clinical data show that the identified inflammatory signature, together with a MYB expression, predicts lung metastasis relapse in BC patients. These results demonstrate that the c-Myb-regulated transcriptional program in BCs results in a blunted inflammatory response and consequently suppresses lung metastasis.

Effector Regulatory T Cell Differentiation and Immune Homeostasis Depend on the Transcription Factor Myb

FoxP3-expressing regulatory T (Treg) cells are essential for maintaining immune homeostasis. Activated Treg cells undergo further differentiation into an effector state that highly expresses genes critical for Treg cell function, although how this process is coordinated on a transcriptional level is poorly understood. Here, we demonstrate that mice lacking the transcription factor Myb in Treg cells succumbed to a multi-organ inflammatory disease. Myb was specifically expressed in, and required for the differentiation of, thymus-derived effector Treg cells. The combination of transcriptome and genomic footprint analyses revealed that Myb directly regulated a large proportion of the gene expression specific to effector Treg cells, identifying Myb as a critical component of the gene regulatory network controlling effector Treg cell differentiation and function.

Dual Mechanism of Rag Gene Repression by c-Myb during Pre-B Cell Proliferation

Developing B lymphocytes undergo clonal expansion following successful immunoglobulin heavy chain gene rearrangement. During this proliferative burst, expression of the Rag genes is transiently repressed to prevent the generation of double-stranded DNA (dsDNA) breaks in cycling large pre-B cells. The Rag genes are then reexpressed in small, resting pre-B cells for immunoglobulin light chain gene rearrangement. We previously identified c-Myb as a repressor of Rag transcription during clonal expansion using Abelson murine leukemia virus-transformed B cells. Nevertheless, the molecular mechanisms by which c-Myb achieved precise spatiotemporal repression of Rag expression remained obscure. Here, we identify two mechanisms by which c-Myb represses Rag transcription. First, c-Myb negatively regulates the expression of the Rag activator Foxo1, an activity dependent on M303 in c-Myb's transactivation domain, and likely the recruitment of corepressors to the Foxo1 locus by c-Myb. Second, c-Myb represses Rag transcription directly by occupying the Erag enhancer and antagonizing Foxo1 binding to a consensus forkhead site in this cis-regulatory element that we show is crucial for Rag expression in Abelson pre-B cell lines. This work provides important mechanistic insight into how spatiotemporal expression of the Rag genes is tightly controlled during B lymphocyte development to prevent mistimed dsDNA breaks and their deleterious consequences.

Plasma Triglyceride Levels May Be Modulated by Gene Expression of IQCJ, NXPH1, PHF17 and MYB in Humans

A genome-wide association study (GWAS) by our group identified loci associated with the plasma triglyceride (TG) response to ω-3 fatty acid (FA) supplementation in IQCJ, NXPH1, PHF17 and MYB. Our aim is to investigate potential mechanisms underlying the associations between single nucleotide polymorphisms (SNPs) in the four genes and TG levels following ω-3 FA supplementation. 208 subjects received 3 g/day of ω-3 FA (1.9–2.2 g of EPA and 1.1 g of docosahexaenoic acid (DHA)) for six weeks. Plasma TG were measured before and after the intervention. 67 SNPs were selected to increase the density of markers near GWAS hits. Genome-wide expression and methylation analyses were conducted on respectively 30 and 35 participants’ blood sample together with in silico analyses. Two SNPs of IQCJ showed different affinities to splice sites depending on alleles. Expression levels were influenced by genotype for one SNP in NXPH1 and one in MYB. Associations between 12 tagged SNPs of IQCJ, 26 of NXPH1, seven of PHF17 and four of MYB and gene-specific CpG site methylation levels were found. The response of plasma TG to ω-3 FA supplementation may be modulated by the effect of DNA methylation on expression levels of genes revealed by GWAS.

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".

MYB-QKI rearrangements in angiocentric glioma drive tumorigenicity through a tripartite mechanism

Angiocentric gliomas are pediatric low-grade gliomas (PLGGs) without known recurrent genetic drivers. We performed genomic analysis of new and published data from 249 PLGGs including 19 Angiocentric Gliomas. We identified MYB-QKI fusions as a specific and single candidate driver event in Angiocentric Gliomas. In vitro and in vivo functional studies show MYB-QKI rearrangements promote tumorigenesis through three mechanisms: MYB activation by truncation, enhancer translocation driving aberrant MYB-QKI expression, and hemizygous loss of the tumor suppressor QKI. This represents the first example of a single driver rearrangement simultaneously transforming cells via three genetic and epigenetic mechanisms in a tumor.

Transcription factor motif quality assessment requires systematic comparative analysis

Transcription factor (TF) binding site prediction remains a challenge in gene regulatory research due to degeneracy and potential variability in binding sites in the genome. Dozens of algorithms designed to learn binding models (motifs) have generated many motifs available in research papers with a subset making it to databases like JASPAR, UniPROBE and Transfac. The presence of many versions of motifs from the various databases for a single TF and the lack of a standardized assessment technique makes it difficult for biologists to make an appropriate choice of binding model and for algorithm developers to benchmark, test and improve on their models. In this study, we review and evaluate the approaches in use, highlight differences and demonstrate the difficulty of defining a standardized motif assessment approach. We review scoring functions, motif length, test data and the type of performance metrics used in prior studies as some of the factors that influence the outcome of a motif assessment. We show that the scoring functions and statistics used in motif assessment influence ranking of motifs in a TF-specific manner. We also show that TF binding specificity can vary by source of genomic binding data. We also demonstrate that information content of a motif is not in isolation a measure of motif quality but is influenced by TF binding behaviour. We conclude that there is a need for an easy-to-use tool that presents all available evidence for a comparative analysis.

Transcription factor motif quality assessment requires systematic comparative analysis

Transcription factor (TF) binding site prediction remains a challenge in gene regulatory research due to degeneracy and potential variability in binding sites in the genome. Dozens of algorithms designed to learn binding models (motifs) have generated many motifs available in research papers with a subset making it to databases like JASPAR, UniPROBE and Transfac. The presence of many versions of motifs from the various databases for a single TF and the lack of a standardized assessment technique makes it difficult for biologists to make an appropriate choice of binding model and for algorithm developers to benchmark, test and improve on their models. In this study, we review and evaluate the approaches in use, highlight differences and demonstrate the difficulty of defining a standardized motif assessment approach. We review scoring functions, motif length, test data and the type of performance metrics used in prior studies as some of the factors that influence the outcome of a motif assessment. We show that the scoring functions and statistics used in motif assessment influence ranking of motifs in a TF-specific manner. We also show that TF binding specificity can vary by source of genomic binding data. Finally, we demonstrate that information content of a motif is not in isolation a measure of motif quality but is influenced by TF binding behaviour. We conclude that there is a need for an easy-to-use tool that presents all available evidence for a comparative analysis.