Advertisement
Research Article
Transparent Process
Open Access

Dependence on Myb expression is attenuated in myeloid leukaemia with N-terminal CEBPA mutations

Giacomo Volpe  Correspondence email, Pierre Cauchy, David S Walton, Carl Ward, Daniel Blakemore, Rachael Bayley, Mary L Clarke, Luisa Schmidt, Claus Nerlov, Paloma Garcia, Stéphanie Dumon, Florian Grebien, Jon Frampton  Correspondence email
Published 15 March 2019. DOI: 10.26508/lsa.201800207

Abstract

Mutations at the N- or C-terminus of C/EBPα are frequent in acute myeloid leukaemia (AML) with normal karyotype. Here, we investigate the role of the transcription factor Myb in AMLs driven by different combinations of CEBPA mutations. Using knockdown of Myb in murine cell lines modelling the spectrum of CEBPA mutations, we show that the effect of reduced Myb depends on the mutational status of the two Cebpa alleles. Importantly, Myb knockdown fails to override the block in myeloid differentiation in cells with biallelic N-terminal C/EBPα mutations, demonstrating for the first time that the dependency on Myb is much lower in AML with this mutational profile. By comparing gene expression following Myb knockdown and chromatin immunoprecipitation sequencing data for the binding of C/EBPα isoforms, we provide evidence for a functional cooperation between C/EBPα and Myb in the maintenance of AML. This co-dependency breaks down when both alleles of CEBPA harbour N-terminal mutations, as a subset of C/EBPα-regulated genes only bind the short p30 C/EBPα isoform and, unlike other C/EBPα-regulated genes, do so without a requirement for Myb.

Introduction

Acute myeloid leukaemia (AML), one of the most common and deadliest forms of proliferative neoplasms, is established through a stepwise acquisition of genetic and epigenetic alterations that result in the malignant transformation of haematopoietic progenitor cells (Kelly & Gilliland, 2002; Moore, 2005). Often, AML arises through the collaboration between mutations affecting transcription factors (e.g., CEBPA, PU.1, and RUNX1) and signalling proteins (such as FLT3, RAS, and KIT) that lead to an aberrant proliferation capacity coupled with a disruption of terminal myeloid differentiation (Tenen, 2003; Rosenbauer & Tenen, 2007).

C/EBPα, a leucine zipper transcription factor with a known tumour suppressor function, has been demonstrated to play an important role in granulocytic development and in the maintenance of haematopoietic stem cell homeostasis (Porse et al, 2001, 2005; Zhang et al, 2004; Koschmieder et al, 2009; Welner et al, 2013; Ye et al, 2013). C/EBPα is translated as two major isoforms, namely a full-length 42-kD form (p42) and a truncated 30-kD protein (p30) that arises from a downstream translational initiation codon (Lin et al, 1993). Mutations in the CEBPA gene are frequently associated with leukaemia, being found in 8–14% of all de novo AML with normal karyotype (Nerlov, 2004; Leroy et al, 2005; Song et al, 2015) and typically involve both alleles. C/EBPα-mutant proteins are classified into two major groups: (i) C-terminal insertions or deletions within the basic region leucine zipper DNA-binding domain; and (ii) N-terminal mutations that lead to the complete ablation of p42 while retaining normal p30 function (Pabst et al, 2001; Leroy et al, 2005; Fasan et al, 2014). Most patients carrying CEBPA mutations harbour one allele with an N-terminal mutation and one with a C-terminal mutation, with homozygosity for N- or C-terminal mutations being less common (Gombart et al, 2002; Pabst & Mueller, 2007). Furthermore, several reports have demonstrated that biallelic mutations of CEBPA are associated with a favourable outcome, when not found in association with FLT3-activating mutations (Renneville et al, 2009; Dufour et al, 2010).

Efforts aimed at understanding how mutations or oncoproteins may cooperate in driving the leukaemogenesis have pointed to cooperation between C/EBPα and other transcription factors, such as RUNX1, MYB, and PU.1. We have previously demonstrated the functional cooperation of Myb and C/EBPα in the regulation of the Flt3 gene in both haematopoietic and leukaemia stem cells (Volpe et al, 2013, 2015). Our studies indicated that Myb and C/EBPα act cooperatively through their combined activity on promoter and intronic elements in the Flt3 gene (Volpe et al, 2013). Furthermore, we reported a strong linear correlation between expression of the two transcription factors and FLT3 RNA levels in human CN-AML, adding to an increasing body of evidence that points to MYB being a crucial component of leukaemia maintenance and oncogene addiction (Hess et al, 2006; Zuber et al, 2011; Clarke et al, 2017).

Our findings on the cooperation of Myb and C/EBPα in Flt3 gene regulation prompted us to investigate the global extent of this cooperation in leukaemia and to determine how manipulation of Myb expression might impact on the maintenance of C/EBPα-driven leukaemia. To address this, we performed genetic manipulation studies in murine haematopoietic progenitor cell lines harbouring either wild-type C/EBPα or the most frequently occurring combinations of biallelic CEBPA mutations, that is Nter/Nter or Nter/Cter to determine the biological and molecular consequences of reduced Myb activity on the leukaemia driven by those mutations. Here, we show that reducing Myb activity can override the differentiation barrier, although the dependency on Myb expression generally observed in leukaemia is minimal in the presence of CEBPA biallelic N-terminal mutations.

Materials and Methods

Cell lines

Cells were cultured in RPMI medium supplemented with 10% fetal bovine serum, 50U/ml penicillin, 50 μg/ml streptomycin, and 2 mM l-glutamine. The culture of FMH9 cells (Volpe et al, 2013) were supplemented with 50 ng/ml stem cell factor, 5 ng/ml GM-CSF, 5 ng/ml IL-3 (IL3), and 5 ng IL-6 (IL6), whereas KL cells (GV, JF, and FG, unpublished) and LL cells (Grebien et al, 2015; Schmidt et al, 2019) required 2 ng/ml IL3. Both KL and LL have been established by serial replating of E14.5 foetal liver cells obtained from mice homozygous for the Lp30 allele (Kirstetter et al, 2008) or carrying both Lp30 and K313KK alleles (Bereshchenko et al, 2009). Briefly, the KL and LL cell lines were obtained by performing six rounds of replating in M3434 semisolid medium (Stem Cell Technologies Inc) followed by an initial liquid culture for 4 wk in the presence of stem cell factor (50 ng/ml), IL3 (2 ng/ml), and IL6 (2 ng/ml). After this initial period, the cells were transferred into the culture medium described above. All cytokines were purchased from Peprotech EC.

Transfection experiments, cell viability, proliferation, apoptosis, and differentiation assays

In total, 5 × 106 FMH9, KL, or LL cells were electroporated with 300 mM of Myb siRNA (s70212, Ambion; Life Technologies) or a scrambled negative control siRNA (4390843 Silencer Select Negative Control #1; Life Technologies) using an Amaxa 4D-nucleofector with solution SF Cell Line (V4XC-2024; Lonza) and program EO-100 for FMH9 cells or solution P3 Primary Cells (V4XP-3024; Lonza) and program DS-120 for KL and LL cells.

After transfection, the cells were plated at a density of 106 cells/ml and viable cells counted and passaged at a ratio of 1:2 every 24 h for 4 consecutive days. Cell cycle analysis was performed by labelling transfected cells (48 h post-nucleofection) with 10 μM BrdU for 1 h. Cells were co-stained with 7-AAD (A9400-1mg; Sigma-Aldrich) and BrdU using the BrdU flow kit (8811-6600; BD Bioscience) according to the manufacturer’s instructions, as previously described (Bayley et al, 2018). Apoptosis analysis was performed using the Annexin V kit (eBioscience) as previously described (Volpe et al, 2015). The percentage of apoptotic cells was obtained by performing live cells gating. Proliferation analysis was performed using CellTrace carboxyfluorescein succinimidyl ester (CFSE) Cell Proliferation Kit (C34554; Thermofisher Scientific). Assessment of differentiation following Myb knockdown was achieved by flow cytometry/immunofluorescence staining of the cells with anti-CD11b PE-Cy7 (25-0112-81; eBioscience) anti-Gr-1 APC (14-5921-82; eBioscience), anti-CD135 PE (12-1351-81; eBioscience), and anti-CD117 PE-Cy5 (15-1171-82; eBioscience). Acquisition and analysis of flow cytometric data were performed using Cyan ADP with either Summit 4.4 software (Beckman Coulter) or FlowJo software (FlowJo, LLC).

Quantitative reverse transcriptase polymerase chain reaction (RT-PCR) analysis

106 cells from each line were harvested 24 h post-transfection. RNA was extracted using RNeasy Mini kit (QIAGEN), and first-strand cDNA synthesis was performed using standard protocols. Quantitative PCR reactions were performed using predesigned Taqman gene expression assays as previously described (Volpe et al, 2015).

Statistical analysis

Statistical significance was determined by performing t test for pairwise comparison, and the P-values are indicated where appropriate. Analysis of MYB expression in human patient array data presented in Fig 1A was performed using non-parametric Kruskal–Wallis test. All statistical analyses were performed using GraphPad Prism 7 (GraphPad Software Inc).

Figure 1.
Figure 1. Myb expression is required for the proliferation of CEBPA biallelic mutant cell lines.

(A) Scatter plot depicting the abundance of MYB transcript in subgroups of patients from the Verhaak et al (2009) dataset, characterised by the molecular abnormalities indicated in the graph. Statistical significance presented in this plot has been calculated using non-parametric the Kruskal–Wallis test. (B) Bar plot representing Myb mRNA quantification by quantitative RT-PCR in FMH9, KL, and LL cell lines, normalised against B2m house-keeping gene results. Statistical analysis was performed using t test (**P < 0.01 and *P < 0.05). (C) Quantitative RT-PCR of Myb transcript abundance in FMH9, KL, and LL cells 24 h post-transfection with Myb siRNA. Expression is normalised to B2m and standardised to the control samples. Error bars represent the SEM and numbers are plotted as mean ± SEM. Each plot is representative of six independent experiments (***P < 0.001 and *P < 0.05). (D) Bar plot indicating the cell viability and relative proliferation of FMH9, KL, and LL cells after Myb siRNA transfection relative to the corresponding negative control (***P < 0.001 and **P < 0.01, *P < 0.05). This plot represents an average of six independent experiments. (E) Flow cytometric analysis of cellular proliferation by CFSE incorporation in FMH9, KL, and LL cells after Myb siRNA transfection relative to the corresponding negative control. Continuous lines indicate cells transfected with the negative control siRNA, whereas dashed lines indicate siMyb-transfected cells. The statistical analysis was performed using t test on the geometric means of fluorescence intensity for each time point comparing siNEG versus siMyb-treated cells, as indicated by the colour-matched bar on top of every peak (**P < 0.01 and *P < 0.05). Each histogram is representative of six independent experiments. (F) Flow cytometric analysis of the cell cycle in FMH9, KL, and LL cells was performed by staining with 7-AAD at 48 h after Myb siRNA transfection and compared with the negative control. Percentages of cells in G0/G1 are indicated in each histogram. Each plot is indicative of six independent experiments. (G) Representative bar plot showing apoptosis analysis performed by Annexin V staining 72 h post Myb siRNA transfection in FMH9, KL, and LL cells (***P < 0.001 and **P < 0.01). Each bar plot represents an average of six independent experiments.

Bioinformatic analyses

Details of the methodologies used to analyse and compare RNA-seq, microarray, and ChIP-seq data, including Gene Ontology (GO) and Gene Set Enrichment Analysis (GSEA) comparisons, can be found in the Supplementary Information.

RNA-sequencing

For RNA-Seq, libraries were prepared using the Illumina TruSeq Stranded kit according to the manufacturer’s instructions. Sequencing was performed at the Institute of Medical Biochemistry, University of Veterinary Medicine, Vienna, Austria, and in Genomics Birmingham, University of Birmingham, Birmingham, UK, on Illumina HiSeq 2500 and NextSeq 500 sequencers, respectively.

Data availability

RNA-Seq data generated in this study are available at the Gene Expression Omnibus under series GSE119348.

Results

High MYB expression is associated with CEBPA mutations in AML

Previous studies have provided evidence for functional cooperation between C/EBPα and MYB in activating the expression of key genes for both haematopoietic and leukaemia stem cell functions in mouse and human (Volpe et al, 2013, 2015, 2017). Using publicly available AML patient profiling arrays (Verhaak et al, 2009) and focussing on specific molecular abnormalities found in cytogenetically normal AML (CN-AML), we observed that MYB expression levels are highest in those patients carrying biallelic CEBPA mutations (Fig 1A).

Given the lack of a suitable cellular system to investigate the relationship between MYB and C/EBPα in human AML, we used murine cell lines modelling the spectrum of CEBPA mutations. Bereshchenko and coworkers demonstrated that mutations in the C/EBPα protein efficiently drive leukaemia in vivo and that the combination of N- and C-terminal mutations were the most highly leukaemogenic, whereas biallelic C-terminal mutation resulted in the longest latency (Kirstetter et al, 2008; Bereshchenko et al, 2009). We generated cell lines by performing serial replating in semisolid medium of E14.5 foetal liver cells carrying knock-in modifications mimicking either the N-terminal LP30 (L allele) or the C-terminal K313KK (K allele) CEBPA mutations (Grebien et al, 2015). LL cells were used to represent AML harbouring biallelic N-terminal CEBPA mutations, whereas cells carrying one K mutation and one L mutation (KL cells) provided a model of leukaemia with both N- and C-terminal CEBPA mutations. AML with wild-type C/EBPα expression was modelled using a previously characterised myelomonocytic leukaemia cell line, namely, FMH9 (Volpe et al, 2013), which was established by ectopic expression of HoxA9 and Meis1 in bone marrow haematopoietic progenitor cells. The phenotype of these cells lines was investigated by determining the surface expression of several myeloid markers (CD11b, Gr1, Kit, and Flt3) (Fig S1A). Importantly, mRNA quantification indicated that both FMH9 and LL cells displayed similar Myb mRNA levels, whereas KL cells exhibited a higher level of Myb expression (Fig 1B), thus being in agreement with the observations obtained from the patient array data.

Figure S1.
Figure S1.

(A) Phenotypic characterisation of CD11b, Gr-1, Flt3, and Kit surface marker expression in FMH9, KL, and LL cells. (B) Related to Fig 1D. Bar plot depicting the siNEG/siMyb proliferation ratio comparison between FMH9 control and CEBPA-mutant cell lines. These data represent an average of six independent experiments. Statistical analysis was performed using t test (***P < 0.001, **P < 0.01). (C) Related to Fig 1E. Bar plots representing the average of CFSE geometric mean intensities (as indicated by the bar on top of every peak in Fig 1E) of FMH9, KL, and LL cells and their relative changes upon siMyb-mediated knockdown. Filled bars represent siNEG-treated cells, whereas empty bars indicated siMyb-treated cells. Geometric mean intensities have been calculated using FlowJo software and have been used to perform statistical analysis using the t test (**P < 0.01, *P < 0.05). These data represent an average of six independent experiments.

Manipulation of Myb expression does not reverse the differentiation block in cells carrying biallelic N-terminal CEBPA mutations

To investigate the requirement for Myb in the maintenance of CEBPA-driven leukaemia, we performed siRNA-mediated knockdown of Myb. Cells were transfected with siRNAs targeting either Myb or a scrambled negative control and were harvested after 24 h to determine the efficiency of knockdown. This analysis revealed a decrease in Myb transcripts by 60–80% in all cell lines (Fig 1C).

To determine the biological consequences of Myb knockdown cells were cultured for up to 96 h and cell numbers determined daily. Myb down-regulation induced growth retardation in the three cell lines, regardless of the CEBPA mutational status; albeit those cells showing a similar pattern, the growth defect observed in FMH9 cells was significantly more pronounced in comparison with the CEBPA mutant cells lines (Figs 1D and S1B). Flow cytometric analysis of CSFE dye dilution revealed that Myb knockdown induced a significant proliferation defect in FMH9 cells, whereas both KL and LL cells were unaffected (Figs 1E and S1C). The pattern of the proliferation defect observed in FMH9 cells 72 h post Myb down-regulation appeared to be bimodal, probably because of cells being induced to differentiate or due to the transient effect of the siRNA-mediated knockdown. By combined staining with 7-AAD and BRDU 48 h post-transfection, we observed that Myb knockdown led to a significant increase in the proportion of cells in the G0/G1 phase with a concomitant decrease in both S and G2/M phases in FMH9, whereas KL and LL cells showed no alteration in their capacity to progress through the cell cycle (Figs 1F and S2A, and B). This analysis also revealed an increase in the percentage of cells with less than 2n DNA content in CEBPA-mutant cell lines only, this being indicative of cells undergoing apoptosis/necrosis (Fig S2C). To confirm this observation, we performed Annexin V staining. This revealed a significant increase in the extent of apoptotic/necrotic cells in both KL and LL cells, whereas no increase was observed in FMH9 cells (Figs 1G and S3A).

Figure S2.
Figure S2. Related to Fig 1F and G.

(A) Phenotypic two-dimensional dot-plots representing the percentage of FMH9, KL, and LL cells in the G0/G1, S, and G2/M phases of the cell cycles and their changes in response to Myb siRNA-mediated down-regulation. The percentage of each population is indicated in every histogram. (B) Bar plots depicting the changes in the percentages of FMH9, KL, and LL cells in the G0/G1, S, and G2/M phases of the cell cycle as shown in the top panel. These data represent an average of six independent experiments. Statistical analysis was performed using t test (***P < 0.001, **P < 0.01). (C) Representative histograms showing cell cycle profile changes in response to Myb siRNA-mediated knockdown in FMH9, KL, and LL. The percentage of cells in G0/G1 and apoptotic/necrotic cells is indicated in each histogram. This panel is indicative of six independent experiments.

Figure S3.
Figure S3. Related to Fig 1G.

(A) Two-dimensional dot-plots representing the percentage of Annexin V+ cells in FMH9, KL, and LL cells in response to Myb siRNA-mediated down-regulation. Each dot plot is representative of six independent experiments. (B) Bar plots showing quantitative PCR analysis of Bcl2 and Bim expression in FMH9, KL, and LL cells in response to Myb knockdown in four independent experiments. Expression is normalised to β2m and standardised to the siNEG-treated samples for every cell line. Error bars represent the standard error or the mean. Statistical analysis was performed using t test (*P < 0.05).

Because it is accepted that homozygous CEBPA mutations lead to a block in myeloid lineage commitment, we investigated how Myb knockdown affects the differentiation capacity of cells in the presence of either wild-type or biallelic mutant C/EBPα. At 96 h, both FMH9 and KL cells exhibited a clear induction of myeloid differentiation as indicated by increased expression of Gr-1 and CD11b (Fig 2). However, this phenomenon was not observed in LL cells, which is intriguing because it suggests that leukaemia cells carrying biallelic N-terminal CEBPA mutations have a reduced dependency on Myb expression in respect to differentiation control.

Figure 2.
Figure 2. Suppression of Myb expression overrides myeloid differentiation block in FMH9, KL but not in LL cells.

Two-dimensional flow cytometric dot plot representing the analysis of CD11b and Gr-1 myeloid surface markers expression in FMH9, KL, and LL transfected with either Myb siRNA or the corresponding control. The percentage of double-positive cells is indicated in every plot. The box plots in the right panel shows an average of six independent experiments. Statistical significance was calculated using t test (***P < 0.001).

Molecular consequences of Myb manipulation in cell lines carrying either wild-type or mutant CEBPA

Myb was previously reported to suppress myeloid commitment and to promote self-renewal in haematopoietic progenitors and leukaemia cells (Lorenzo et al, 2011; Zhao et al, 2011). We wished to explore how the interplay between Myb and wild-type or mutant C/EBPα influences the transcriptome, so we performed RNA-seq analysis following Myb knockdown in the context of wild-type or mutant C/EBPα. In line with previous reports, inspection of the RNA-seq datasets for FMH9, KL, and LL cells showed the expected patterns for typical myeloid genes that are known to be Myb targets, such as Gfi1, Itgam/CD11b, Gr-1/Ly6d, and Ccnd2 (Fig S4A and B).

Figure S4.
Figure S4. Related to Fig 3.

(A) University of California at Santa Cruz (UCSC) genome browser screenshots of RNA-Seq performed in FMH9, KL, and LL cells with siMyb and siNEG treatment at differentially regulated genes. Profiles scaled to 1% GAPDH. (B) Absolute (top) and siNEG-relative (bottom) FPKM quantification of expression levels of genes are shown.

Analysis of the RNA-seq data by global correlation clustering of fragments per kilobase of transcript per million mapped reads (FPKM) in steady-state conditions (i.e., after control siRNA treatment only) revealed higher similarity between KL and LL cells, with FMH9 cells clustering on their own, in line with previous reports on biallelic CEBPA mutants (Wouters et al, 2009; Taskesen et al, 2011) (Fig 3A). Differential gene expression analysis comparing control and siMyb-transfected FMH9, KL, and LL cells revealed 790, 1217, and 40 genes being down-regulated, whereas 1364, 1668, and 329 genes were up-regulated (Fig S5A). This demonstrates that, whereas FMH9 and KL cells display a large number of genes responsive to Myb knockdown, LL cells exhibited only minor transcriptomic responses. Intersection of the gene expression changes showed generally a greater overlap of up-regulated genes than down-regulated genes (Fig 3B), whereas hierarchical clustering revealed a higher degree of similarity between gene expression changes in FMH9 and KL cells over changes observed in LL cells (Fig 3C). These findings confirm our hypothesis that the LL-mutant leukaemia phenotype is less dependent on Myb.

Figure 3.
Figure 3. Myb down-regulation causes concomitant differential regulation of leukaemia gene expression programmes in both FMH9 and KL but not LL cells.

(A) Spearman correlation clustering of steady-state, control scrambled negative siRNA-transfected FMH9, KL, and LL cells. (B) Venn diagram overlaps of differentially expressed genes in FMH9, KL, and LL cell lines following Myb knockdown. Left and right: significantly down- and up-regulated genes, respectively. (C) Hierarchical clustering of log2 fold changes resulting from siMyb treatment in FMH9, KL, and LL cells. (D) Average siMyb/siNEG log2 fold changes for leukaemia-relevant GO classes. (E) RT-qPCR gene expression analysis of differentiation, apoptosis, and cell cycle genes post control and siMyb transfection. Relative expression values are presented as ± SEM. Statistical analysis was performed using t test (***P < 0.001, **P < 0.01, and *P < 0.05). Each bar plot represents an average of six independent experiments.

Figure S5.
Figure S5. Related to Fig 3.

(A) Boxplot showing FPKM gene expression levels of genes differentially regulated in FMH9, KL, and LL cells following Myb knockdown. (B) Heat maps showing gene expression fold changes at individual genes from GO categories depicted in Fig 3D. The detailed list of genes from each pathway/group is provided in Table S2.

We next investigated whether Myb knockdown results in differential regulation of relevant gene ontologies. We examined negative regulation of G1 to S phase transition (Ashburner et al, 2000), leukaemia stem cell state (Gal et al, 2006), response to AML treatment (Bogni et al, 2006), and terminal myeloid differentiation, computing average log2 siMyb/control fold changes for each ontologies (Table S1). We found that FMH9 and KL cells, but not LL cells, exhibited down-regulation of the leukaemia stem cell gene expression programme following Myb knockdown (Figs 3D and S5B). Crucially, this analysis also showed that specifically for FMH9 and KL cells, siMyb treatment resulted in up-regulation of genes that are also up-regulated during the treatment of AML.

To confirm the observations obtained from the RNA-seq, we performed quantitative RT-PCR analysis of selected key genes. Significantly, Cebpa mRNA was down-regulated upon Myb knockdown in both FMH9 and KL cells, but not in LL cells. This suggests that the functional cross-regulation between Myb and C/EBPα could be lost in the presence of biallelic N-terminal CEBPA mutations. This analysis confirmed specific regulation by Myb of differentiation-related genes in FMH9 and KL cells but not in LL cells (Gfi1 and Sbno2), changes in the expression of genes leading to a negative impact on proliferation in FMH9 cells (Ccnd2 and Sema4d), and down- and up-regulation of genes related to increased apoptosis seen in KL and LL cells (Prune2 and Dusp1) (Figs 3E and S4B). Moreover, inspection of known Myb target genes revealed a significant repression of Bcl2 (Taylor et al, 1996; Salomoni et al, 1997), a known anti-apoptotic regulator, in KL cells; concomitantly, we observed the up-regulation of a pro-apoptotic gene normally anti-correlated with Myb in AML, namely, Bcl2l11 (Bim) (Jing et al, 2015), in LL cells. This is in agreement with the strong induction of apoptosis observed in these cells upon Myb knockdown (Fig S3A and B).

Overall, the transcriptome changes following Myb knockdown in cell lines with different CEBPA mutational status are consistent with the corresponding phenotypic changes and further demonstrate that cells with the biallelic LL C/EBPα configuration lack a major dependence on Myb.

Coincident binding of C/EBPα p42 and Myb correlates with gene repression by Myb, whereas genes positively controlled by Myb tend to bind C/EBPα p30

To shed light on the possible interplay between the mutational status of C/EBPα and the Myb-dependent regulation of C/EBPα target genes, we set out to investigate the chromatin binding properties and transcriptional effects of C/EBPα p30, C/EBPα p42, and K313KK-mutant C/EBPα isoforms. The N-terminal–mutant L allele (Kirstetter et al, 2008) leads to the expression of the p30 isoform only, whereas the C/EBPα K313KK–mutant allele gives rise to a C-terminal mutant that disables the DNA-binding domain, resulting in a block in differentiation (Bereshchenko et al, 2009). Because the binding dynamics of C/EBPα in double-mutant cells would be technically difficult to characterise, we used previously published chromatin immunoprecipitation sequencing (ChIP-Seq) data from a single isoform transfection model (Grebien et al, 2015). This study used HA-tagged Cebpa constructs transfected into the FDCP1 cell line (Bernard et al, 1991) and an immortalised IL3-dependent murine myeloid cell line that approximates to the wild-type C/EBPα leukaemia line FMH9 that we used for the Myb knockdown studies. This latter study concurrently provided gene expression microarray analysis of mock- and C/EBPα-transfected cells. As expected, there were significant increases in Cebpa transcript abundance following overexpression of C/EBPα p30, C/EBPα p42, and C/EBPα K313KK as compared with mock transfection (Fig S6A, compare with Fig S6B in FMH9, KL, and LL cells). Analysis of ChIP-seq data revealed 17,452, 22,873, and 68,432 peaks for C/EBPα p30-, C/EBPα p42-, and C/EBPα K313KK–transfected cells, respectively, which were mostly located in intergenic and intronic regions (Fig S6C and D). However, visual inspection of the C/EBPα K313KK dataset revealed low signal in K313KK peaks, hinting that those are not specific, consistent with the loss of binding due to the K313KK mutation in the DNA-binding domain (Grebien et al, 2015; Bereshchenko et al, 2009). We next characterised specific p30 and p42 peaks by ranking tag counts around merged summits by p30/p42 fold change and identified 3,585 p42-specific, 19,949 shared, and 4,421 p30-specific peaks (Fig 4A). Although we did not use the C/EBPα K313KK dataset as a direct base for comparison because of the deleterious effect of the K allele on DNA binding, and thus low signal to noise ratio, its binding pattern was mostly located in shared sites. Crucially, by retrieving tag counts for Myb ChIP-Seq datasets, we observed that Myb binding largely parallels that of p42 binding, both in FDCP1 and in an MLL-AF9/NrasG12D murine AML cell line (Roe et al, 2015), although there was an overlap with some regions that predominantly bind p30 (Fig 4A). GO analysis of the p42-specific peaks revealed signalling pathways involved in haematopoietic homeostasis and pro-apoptotic genes (Fig 4B). Conversely, p30-specific peaks were enriched in pluripotency genes, consistent with the leukaemia stem cell signature seen amongst the up-regulated genes following Myb knockdown in LL cells (Figs 3D and 4B, and S5B).

Figure S6.
Figure S6. Related to Fig 4.

(A) Cebpa transcript abundances in FDCP1 cells transfected with mock, p30, and p42 C/EBPα isoforms measured by RNA-Seq, relative to B2m. (B) Cebpa transcript abundances in FMH9, KL, and LL cells measured by quantitative RT-PCR, relative to B2m. This bar plot represents an average of six independent experiments. Statistical analysis was performed using t test (***P < 0.001, **P < 0.01, *P < 0.05). (C) UCSC genome browser screenshot of C/EBPα p30, p42, and K313KK ChIP-Seq datasets in the FDCP1 cell line at the Runx1 gene locus. (D) Stacked percentage bar plots showing genomic annotation of C/EBPα p30, p42, and K313KK ChIP-seq peaks in the FDCP1 cell line.

Figure 4.
Figure 4. p42 C/EBPα binding is linked with gene activation increased following Myb knockdown, whereas p30 C/EBPα binding correlates with gene repression independently of Myb.

(A) Heat maps sorted by C/EBPα p30/p42 tag count fold change of ChIP-seq signals for C/EBPα p42, C/EBPα p42 K313KK, and C/EBPα p30 isoforms in the FDCP1 cell line, as well as for Myb and C/EBPα p42 in the RN2 cell line. (B) GO analyses of p42 and p30 C/EBPα-specific peaks (left, right). (C) Gene set enrichment analyses of C/EBPα p30, C/EBPα p42, and C/EBPα p42-K313KK binding versus cognate-induced fold change (top left, top right, and bottom left, respectively). (A, D) Heat map showing siMyb/siNEG gene expression fold change in FMH9, KL, and LL cells sorted by C/EBPα p30/p42 ChIP-seq tag count fold change as in (A). (A, D, E) Box plots showing quantification of gene expression fold changes from (D) for the nearest genes from groups 1, 2, and 3 defined in (A). Means indicated.

To characterise the consequences of C/EBPα binding on gene expression, we performed GSEA using microarray datasets from FDCP1 cells expressing p30, p42, and K313KK, ranking by log2 fold change against mock transfection. We selected gene sets corresponding to the closest genes of the top 1,000 peaks for cognate ChIP-Seq datasets (p30, p42, and K313KK) in accordance with the constraints of GSEA. We observed significant correlations between p42-induced gene activation and p42 binding, as well as between p30-induced gene repression and p30 binding (Fig 4C). However, C/EBPα K313KK binding was not correlated with changes in gene expression.

We next asked how C/EBPα isoform binding correlates with genes whose expression is altered by Myb knockdown. We plotted gene expression fold changes caused by Myb knockdown in FMH9, KL, and LL cells against the ChIP-seq data ranking of p30/p42 binding. Two broad conclusions arise from this analysis: first, Myb-repressed genes largely bound p42 in the presence or absence of p30; second, genes that are positively regulated by Myb are more predominant amongst the group of genes that preferentially bind p30 (Fig 4D and E). Conversely, p30 binding was linked with down-regulation of gene expression upon Myb knockdown.

We next investigated the sequence content of p30- and p42-specific sites by performing motif discovery analysis in the sequences corresponding to these peaks. This revealed that C/EBP, AP-1, Ets, Myb, and Runx motifs were highly enriched in p42-specific peaks (Fig S7A). However, the C/EBP motif was not enriched in the p30-specific peaks. Instead, these peaks were enriched in CTCF, Nrf, and Ets motifs. To confirm these trends, we plotted motif matches amongst increasing p30/p42 fold change as above. Strikingly, C/EBP motifs seemed to be restricted to p42-specific and shared sites (Fig S7B). Myb, Runx, AP-1, and Ets motifs also followed this trend. Conversely, CTCF motifs were highly enriched in p30-specific sites and to some extent in shared sites corresponding to lower p42 binding. Nrf, Sp1, and CREB motifs also seemed to follow this trend, but not Elf motifs. To verify binding of cognate factors to these motifs as well as active and inactive transcriptional hallmarks, we made use of publicly available AML ChIP-Seq datasets (Yue et al, 2014; Roe et al, 2015). Our analysis revealed that C/EBPβ, an essential transcription factor in normal myeloid development (Goode et al, 2016), which is able to bind C/EBP motifs as well, also co-localised with p42 (Fig S7C). Furthermore, by retrieving tag counts for the transcription activation hallmark p300, we could also show that p42-specific and shared sites, but not p30-specific sites, correspond to putative enhancer sites, consistent with our GSEA analyses for this isoform (Figs 4C and S7D).

Figure S7.
Figure S7. Differential p42, p30 C/EBPα binding with CEBP, MYB motifs and co-localisation with active and repressive transcriptional hallmarks.

(A) Plots showing motif discovery results (x-axis, motif rank; y-axis, −logP) in p42-specific peaks (top) and p-30 specific peaks (bottom). (B) Heat maps showing matches of top p42-specific (top) and p30-specific motifs ranked by p30/p42 C/EBPα ChIP-Seq tag count fold change as in Fig 4A. (A, C) RN2 ChIP-Seq tag counts for transcription factors known to bind motifs from (A) as well p300 sorted as in (C) and Fig 4A. (D) Heat map showing MEL ChIP-Seq tag counts for transcription factors known to bind motifs from Fig S5A, as well as co-activators and co-repressors sorted by p30/p42 C/EBPα tag count fold change as in Fig 4A.

Discussion

In the present study, we have investigated the requirement for the transcription factor Myb in the maintenance of CN-AML driven by different combinations of CEBPA mutations in comparison with leukaemia characterised by the expression of wild-type C/EBPα. We show for the first time that the dependency on Myb is affected by the mutational status of C/EBPα. Compared with cells expressing wild-type C/EBPα, which show a proliferation and differentiation response to enforced reduction in Myb levels, leukaemia driven by biallelic CEBPA mutations exhibits distinct phenotypic responses that are reflected in changes in gene expression. Furthermore, leukaemia with biallelic Nter/Nter CEBPA mutations shows a reduced dependency on Myb, whereas Cter/Nter mutant–driven AML cells are as reliant on Myb as those expressing wild-type C/EBPα but exhibit a quite distinct pattern of phenotype and gene expression changes upon Myb knockdown.

Here, we show that knockdown of Myb in leukaemia cells harbouring either wild-type or Cter/Nter–mutant C/EBPα reverses the abnormal myeloid phenotype normally observed in leukaemia, whereas AML cells carrying biallelic Nter/Nter mutations exhibit persistence of the undifferentiated phenotype. These behaviours are reflected in distinct changes in gene expression. Hence, Myb knockdown in leukaemia cells with wild-type or Cter/Nter–mutant C/EBPα resulted in the loss of a leukaemia stem cell signature and the up-regulation of a gene expression pattern generally observed in patients that are responding to therapeutic treatment. In line with previous reports (Zhao et al, 2014), our analysis demonstrated Gfi1 and Sbno2, both of which encode transcription factors that can be related to the loss of differentiation block, to be positive and negative targets of Myb activity, respectively. The reduced dependency of Nter/Nter–mutant C/EBPα–driven AML cells on Myb is also paralleled by very little change in the number of genes affected by Myb knockdown, including no effect on Gfi1 or Sbno2. Interestingly, Myb targets the expression of the Cebpa gene itself in the wild-type and Cter/Nter C/EBPα contexts, suggesting positive feedback that is not seen in C/EBPα Nter/Nter cells. This result might provide a hint why LL leukaemia displays a different response to Myb manipulation. Analysis of cell cycle and apoptosis also revealed how knockdown of Myb can lead to quite a distinct phenotype; for instance, although no induction in apoptosis/necrosis was observed in C/EBPα wild-type cells, both mutant cell lines displayed a large increase in the percentage of Annexin V+ cells and the appearance of a sub2n population, which is indicative of cells undergoing cell death. However, it is possible that other biological pathways could be affected that lead to such a different response, with the cells perhaps being forced to engage in non-apoptotic fates, such as necroptosis, autophagic cell death, or pyroptosis (Tait et al, 2014); these possibilities remain to be further elucidated.

Considering the different involvement of p30 and p42 in the three leukaemia scenarios we have investigated, the fact that there are p42-only– and p30-only–bound genes indicates that a distinct response should be expected when Myb is reduced, especially when only p30 is present. Most p42- or p42+p30–binding genes also bind Myb, but interestingly, a significant number of p30-only target genes probably do not bind Myb. The genes that bound predominantly p42 fall into GO groups, including those associated with myeloid homeostasis and differentiation, consistent with the finding that Myb knockdown restores myeloid differentiation. Conversely, but in agreement with a persisting leukaemia phenotype not being affected by Myb manipulation, genes bound by p30 include pluripotency genes. This is in line with a previous report that mature cells can tap into stem cell regulatory networks when experiencing mutational hits in key differentiation factors (Soucie et al, 2016).

The vast majority of genes affected by Myb knockdown in the context of Nter/Nter biallelic–mutant C/EBPα are de-repressed and overlap with genes similarly affected in the context of Cter/Nter. That genes such as Dusp1 are not responsive to Myb changes in the context of wild-type C/EBPα presumably means that p42, being fully competent to dimerise and bind DNA is less dependent on Myb. Mapping these gene expression changes onto the profile of p42/p30 binding revealed that they can be expected normally to bind both p42 and p30. The large number of genes affected by Myb knockdown that distinguished the response of Cter/Nter from Nter/Nter are both down-regulated (1,188 genes) and up-regulated (1,423 genes), the majority of the former corresponding to genes that are preferentially bound by p30. In contrast, the up-regulated genes, normally repressed by Myb as seen in the Nter/Nter situation, tend to be genes that exhibit greater binding by p42 when it is present. Intriguingly, in the Nter/Nter situation, in which most of the genes affected by Myb knockdown reflect loss of Myb-dependent repression, these genes also fall into the category of preferential binding to p42. Because no p42 is present in the Nter/Nter leukaemia, this must mean that such genes can be regulated by Myb without a prerequisite for cooperation with C/EBPα p42 or alternatively that another C/EBP family protein such as C/EBPβ can act in place of C/EBPα. In this circumstance, the profile of C/EBPβ binding to genes in the context of myeloid leukaemia cells parallels that of C/EBPα p42. The consequence of this is that most of the genes bound by C/EBPα p42 that are Myb dependent would not exhibit such redundancy of C/EBP protein requirement.

By analysing the occurrence of transcription factor consensus binding motifs in the sequence content of the C/EBPα ChIP-seq peaks, we observed that the regulatory network maintaining the different leukaemia statuses involves different possible sets of transcription factor binding. p42-specific peaks are enriched in C/EBP, Ap1, Myb, Ets, and Runx motifs, whereas the p30-specific peaks contain mostly CTCF, Nrf, and Ets motifs. Importantly, we also observed that p42-specific peaks that are also bound by Myb are highly enriched in the transcription activation hallmark p300. This is to be expected as p300 has been demonstrated to be one of the most important cofactors for Myb and the inhibition of the Myb/p300 interaction is crucial for the maintenance of the leukaemia state (Pattabiraman et al, 2014). The same analysis performed for the p30-specific peaks revealed those sites may be enriched for CTCF in cells that do not express CEBPα p42 (Fig S7D).

In conclusion, we have shown that the nature of mutations in one transcription factor that drives leukaemia can dictate how another leukaemia-associated oncoprotein affects the maintenance of the leukaemia phenotype. The precise nature of the interactions between C/EBPα or its mutated variants and the Myb protein on specific genes that dictate the leukaemia phenotype remain to be elucidated, and it would be fruitful to assess in more detail the relevance of this interaction in human leukaemia patients harbouring those mutations. Added complexity in the involvement of C/EBP proteins in determining the Myb dependency of a given leukaemia most likely goes beyond solely the balance of the C/EBPα isoforms, especially given the possibility that co-expressed factors such as C/EBPβ might compete for binding to C/EBP motifs and heterodimerise with C/EBPα. Ultimately, our findings call for a larger study to determine how manipulating MYB would impact on the maintenance of both murine and human leukaemia driven by other genetic lesions.

Acknowledgements

The work was supported by a Bloodwise Programme Grant (12010) held by J Frampton, S Dumon, and P Garcia and through funding provided by the College of Medical and Dental Sciences of the University of Birmingham. This project has received funding from the European Research Council under the European Union’s Horizon 2020 research and innovation programme (grant agreement no. 636855/StG to F Grebien)

Authors Contributions

  • G Volpe: conceptualisation, data curation, investigation, and writing—original draft.

  • P Cauchy: data curation, formal analysis, and writing—original draft.

  • DS Walton: investigation.

  • C Ward: formal analysis.

  • D Blakemore: investigation.

  • R Bayley: investigation.

  • ML Clarke: investigation.

  • L Schmidt: investigation.

  • C Nerlov: resources.

  • P Garcia: funding acquisition and writing—review and editing.

  • S Dumon: conceptualisation, supervision, and funding acquisition.

  • F Grebien: conceptualisation, formal analysis, investigation, and writing—review and editing.

  • J Frampton: conceptualisation, supervision, funding acquisition, project administration, and writing—review and editing.

Conflict of Interest Statement

The authors declare no conflicts of interest.

  • Received October 7, 2018.
  • Revision received March 7, 2019.
  • Accepted March 7, 2019.
Creative Commons logoCreative Commons logohttps://creativecommons.org/licenses/by/4.0/

This article is available under a Creative Commons License (Attribution 4.0 International, as described at https://creativecommons.org/licenses/by/4.0/).

References

  1. ↵
    1. Ashburner M,
    2. Ball CA,
    3. Blake JA,
    4. Botstein D,
    5. Butler H,
    6. Cherry JM,
    7. Davis AP,
    8. Dolinski K,
    9. Dwight SS,
    10. Eppig JT, et al.
    (2000) Gene ontology: Tool for the unification of biology. The gene ontology consortium. Nat Genet 25: 25–29. doi:10.1038/75556
    OpenUrlCrossRefPubMed
  2. ↵
    1. Bayley R,
    2. Blakemore D,
    3. Cancian L,
    4. Dumon S,
    5. Volpe G,
    6. Ward C,
    7. Almaghrabi R,
    8. Gujar J,
    9. Reeve N,
    10. Raghavan M, et al.
    (2018) MYBL2 supports DNA double strand break repair in hematopoietic stem cells. Cancer Res 78: 5767–5779. doi:10.1158/0008-5472.CAN-18-0273
    OpenUrlAbstract/FREE Full Text
  3. ↵
    1. Bereshchenko O,
    2. Mancini E,
    3. Moore S,
    4. Bilbao D,
    5. Mansson R,
    6. Luc S,
    7. Grover A,
    8. Jacobsen SE,
    9. Bryder D,
    10. Nerlov C
    (2009) Hematopoietic stem cell expansion precedes the generation of committed myeloid leukemia-initiating cells in C/EBPalpha mutant AML. Cancer Cell 16: 390–400. doi:10.1016/j.ccr.2009.09.036
    OpenUrlCrossRefPubMed
  4. ↵
    1. Bernard O,
    2. Li M,
    3. Reid HH
    (1991) Expression of two different forms of fibroblast growth factor receptor 1 in different mouse tissues and cell lines. Proc Natl Acad Sci U S A 88: 7625–7629. doi:10.1073/pnas.88.17.7625
    OpenUrlAbstract/FREE Full Text
  5. ↵
    1. Bogni A,
    2. Cheng C,
    3. Liu W,
    4. Yang W,
    5. Pfeffer J,
    6. Mukatira S,
    7. French D,
    8. Downing JR,
    9. Pui CH,
    10. Relling MV
    (2006) Genome-wide approach to identify risk factors for therapy-related myeloid leukemia. Leukemia 20: 239–246. doi:10.1038/sj.leu.2404059
    OpenUrlCrossRefPubMed
  6. ↵
    1. Clarke M,
    2. Volpe G,
    3. Sheriff L,
    4. Walton D,
    5. Ward C,
    6. Wei W,
    7. Dumon S,
    8. Garcia P,
    9. Frampton J
    (2017) Transcriptional regulation of SPROUTY2 by MYB influences myeloid cell proliferation and stem cell properties by enhancing responsiveness to IL-3. Leukemia 31: 957–966. doi:10.1038/leu.2016.289
    OpenUrlCrossRef
  7. ↵
    1. Dufour A,
    2. Schneider F,
    3. Metzeler KH,
    4. Hoster E,
    5. Schneider S,
    6. Zellmeier E,
    7. Benthaus T,
    8. Sauerland MC,
    9. Berdel WE,
    10. Buchner T, et al.
    (2010) Acute myeloid leukemia with biallelic CEBPA gene mutations and normal karyotype represents a distinct genetic entity associated with a favorable clinical outcome. J Clin Oncol 28: 570–577. doi:10.1200/jco.2008.21.6010
    OpenUrlAbstract/FREE Full Text
  8. ↵
    1. Fasan A,
    2. Haferlach C,
    3. Alpermann T,
    4. Jeromin S,
    5. Grossmann V,
    6. Eder C,
    7. Weissmann S,
    8. Dicker F,
    9. Kohlmann A,
    10. Schindela S, et al.
    (2014) The role of different genetic subtypes of CEBPA mutated AML. Leukemia 28: 794–803. doi:10.1038/leu.2013.273
    OpenUrlCrossRefPubMed
  9. ↵
    1. Gal H,
    2. Amariglio N,
    3. Trakhtenbrot L,
    4. Jacob-Hirsh J,
    5. Margalit O,
    6. Avigdor A,
    7. Nagler A,
    8. Tavor S,
    9. Ein-Dor L,
    10. Lapidot T, et al.
    (2006) Gene expression profiles of AML derived stem cells; similarity to hematopoietic stem cells. Leukemia 20: 2147–2154. doi:10.1038/sj.leu.2404401
    OpenUrlCrossRefPubMed
  10. ↵
    1. Gombart AF,
    2. Hofmann WK,
    3. Kawano S,
    4. Takeuchi S,
    5. Krug U,
    6. Kwok SH,
    7. Larsen RJ,
    8. Asou H,
    9. Miller CW,
    10. Hoelzer D, et al.
    (2002) Mutations in the gene encoding the transcription factor CCAAT/enhancer binding protein alpha in myelodysplastic syndromes and acute myeloid leukemias. Blood 99: 1332–1340. doi:10.1182/blood.v99.4.1332
    OpenUrlAbstract/FREE Full Text
  11. ↵
    1. Goode DK,
    2. Obier N,
    3. Vijayabaskar MS,
    4. Lie ALM,
    5. Lilly AJ,
    6. Hannah R,
    7. Lichtinger M,
    8. Batta K,
    9. Florkowska M,
    10. Patel R, et al.
    (2016) Dynamic gene regulatory networks drive hematopoietic specification and differentiation. Dev Cell 36: 572–587. doi:10.1016/j.devcel.2016.01.024
    OpenUrlCrossRefPubMed
  12. ↵
    1. Grebien F,
    2. Vedadi M,
    3. Getlik M,
    4. Giambruno R,
    5. Grover A,
    6. Avellino R,
    7. Skucha A,
    8. Vittori S,
    9. Kuznetsova E,
    10. Smil D, et al.
    (2015) Pharmacological targeting of the Wdr5-MLL interaction in C/EBPalpha N-terminal leukemia. Nat Chem Biol 11: 571–578. doi:10.1038/nchembio.1859
    OpenUrlCrossRefPubMed
  13. ↵
    1. Hess JL,
    2. Bittner CB,
    3. Zeisig DT,
    4. Bach C,
    5. Fuchs U,
    6. Borkhardt A,
    7. Frampton J,
    8. Slany RK
    (2006) c-Myb is an essential downstream target for homeobox-mediated transformation of hematopoietic cells. Blood 108: 297–304. doi:10.1182/blood-2005-12-5014
    OpenUrlAbstract/FREE Full Text
  14. ↵
    1. Jing D,
    2. Bhadri VA,
    3. Beck D,
    4. Thoms JA,
    5. Yakob NA,
    6. Wong JW,
    7. Knezevic K,
    8. Pimanda JE,
    9. Lock RB
    (2015) Opposing regulation of BIM and BCL2 controls glucocorticoid-induced apoptosis of pediatric acute lymphoblastic leukemia cells. Blood 125: 273–283. doi:10.1182/blood-2014-05-576470
    OpenUrlAbstract/FREE Full Text
  15. ↵
    1. Kelly LM,
    2. Gilliland DG
    (2002) Genetics of myeloid leukemias. Annu Rev Genomics Hum Genet 3: 179–198. doi:10.1146/annurev.genom.3.032802.115046
    OpenUrlCrossRefPubMed
  16. ↵
    1. Kirstetter P,
    2. Schuster MB,
    3. Bereshchenko O,
    4. Moore S,
    5. Dvinge H,
    6. Kurz E,
    7. Theilgaard-Monch K,
    8. Mansson R,
    9. Pedersen TA,
    10. Pabst T, et al.
    (2008) Modeling of C/EBPalpha mutant acute myeloid leukemia reveals a common expression signature of committed myeloid leukemia-initiating cells. Cancer Cell 13: 299–310. doi:10.1016/j.ccr.2008.02.008
    OpenUrlCrossRefPubMed
  17. ↵
    1. Koschmieder S,
    2. Halmos B,
    3. Levantini E,
    4. Tenen DG
    (2009) Dysregulation of the C/EBPalpha differentiation pathway in human cancer. J Clin Oncol 27: 619–628. doi:10.1200/jco.2008.17.9812
    OpenUrlAbstract/FREE Full Text
  18. ↵
    1. Leroy H,
    2. Roumier C,
    3. Huyghe P,
    4. Biggio V,
    5. Fenaux P,
    6. Preudhomme C
    (2005) CEBPA point mutations in hematological malignancies. Leukemia 19: 329–334. doi:10.1038/sj.leu.2403614
    OpenUrlCrossRefPubMed
  19. ↵
    1. Lin FT,
    2. MacDougald OA,
    3. Diehl AM,
    4. Lane MD
    (1993) A 30-kDa alternative translation product of the CCAAT/enhancer binding protein alpha message: Transcriptional activator lacking antimitotic activity. Proc Natl Acad Sci U S A 90: 9606–9610. doi:10.1073/pnas.90.20.9606
    OpenUrlAbstract/FREE Full Text
  20. ↵
    1. Lorenzo PI,
    2. Brendeford EM,
    3. Gilfillan S,
    4. Gavrilov AA,
    5. Leedsak M,
    6. Razin SV,
    7. Eskeland R,
    8. Saether T,
    9. Gabrielsen OS
    (2011) Identification of c-Myb target genes in K562 cells reveals a role for c-Myb as a Master regulator. Genes Cancer 2: 805–817. doi:10.1177/1947601911428224
    OpenUrlCrossRefPubMed
  21. ↵
    1. Moore MA
    (2005) Converging pathways in leukemogenesis and stem cell self-renewal. Exp Hematol 33: 719–737. doi:10.1016/j.exphem.2005.04.011
    OpenUrlCrossRefPubMed
  22. ↵
    1. Nerlov C
    (2004) C/EBPalpha mutations in acute myeloid leukaemias. Nat Rev Cancer 4: 394–400. doi:10.1038/nrc1363
    OpenUrlCrossRefPubMed
  23. ↵
    1. Pabst T,
    2. Mueller BU
    (2007) Transcriptional dysregulation during myeloid transformation in AML. Oncogene 26: 6829–6837. doi:10.1038/sj.onc.1210765
    OpenUrlCrossRefPubMed
  24. ↵
    1. Pabst T,
    2. Mueller BU,
    3. Zhang P,
    4. Radomska HS,
    5. Narravula S,
    6. Schnittger S,
    7. Behre G,
    8. Hiddemann W,
    9. Tenen DG
    (2001) Dominant-negative mutations of CEBPA, encoding CCAAT/enhancer binding protein-alpha (C/EBPalpha), in acute myeloid leukemia. Nat Genet 27: 263–270. doi:10.1038/85820
    OpenUrlCrossRefPubMed
  25. ↵
    1. Pattabiraman DR,
    2. McGirr C,
    3. Shakhbazov K,
    4. Barbier V,
    5. Krishnan K,
    6. Mukhopadhyay P,
    7. Hawthorne P,
    8. Trezise A,
    9. Ding J,
    10. Grimmond SM, et al.
    (2014) Interaction of c-Myb with p300 is required for the induction of acute myeloid leukemia (AML) by human AML oncogenes. Blood 123: 2682–2690. doi:10.1182/blood-2012-02-413187
    OpenUrlAbstract/FREE Full Text
  26. ↵
    1. Porse BT,
    2. Bryder D,
    3. Theilgaard-Monch K,
    4. Hasemann MS,
    5. Anderson K,
    6. Damgaard I,
    7. Jacobsen SE,
    8. Nerlov C
    (2005) Loss of C/EBP alpha cell cycle control increases myeloid progenitor proliferation and transforms the neutrophil granulocyte lineage. J Exp Med 202: 85–96. doi:10.1084/jem.20050067
    OpenUrlAbstract/FREE Full Text
  27. ↵
    1. Porse BT,
    2. Pedersen TA,
    3. Xu X,
    4. Lindberg B,
    5. Wewer UM,
    6. Friis-Hansen L,
    7. Nerlov C
    (2001) E2F repression by C/EBPalpha is required for adipogenesis and granulopoiesis in vivo. Cell 107: 247–258. doi:10.1016/s0092-8674(01)00516-5
    OpenUrlCrossRefPubMed
  28. ↵
    1. Renneville A,
    2. Boissel N,
    3. Gachard N,
    4. Naguib D,
    5. Bastard C,
    6. de Botton S,
    7. Nibourel O,
    8. Pautas C,
    9. Reman O,
    10. Thomas X, et al.
    (2009) The favorable impact of CEBPA mutations in patients with acute myeloid leukemia is only observed in the absence of associated cytogenetic abnormalities and FLT3 internal duplication. Blood 113: 5090–5093. doi:10.1182/blood-2008-12-194704
    OpenUrlAbstract/FREE Full Text
  29. ↵
    1. Roe JS,
    2. Mercan F,
    3. Rivera K,
    4. Pappin DJ,
    5. Vakoc CR
    (2015) BET bromodomain inhibition suppresses the function of hematopoietic transcription factors in acute myeloid leukemia. Mol Cell 58: 1028–1039. doi:10.1016/j.molcel.2015.04.011
    OpenUrlCrossRefPubMed
  30. ↵
    1. Rosenbauer F,
    2. Tenen DG
    (2007) Transcription factors in myeloid development: Balancing differentiation with transformation. Nat Rev Immunol 7: 105–117. doi:10.1038/nri2024
    OpenUrlCrossRefPubMed
  31. ↵
    1. Salomoni P,
    2. Perrotti D,
    3. Martinez R,
    4. Franceschi C,
    5. Calabretta B
    (1997) Resistance to apoptosis in CTLL-2 cells constitutively expressing c-Myb is associated with induction of BCL-2 expression and Myb-dependent regulation of bcl-2 promoter activity. Proc Natl Acad Sci U S A 94: 3296–3301. doi:10.1073/pnas.94.7.3296
    OpenUrlAbstract/FREE Full Text
  32. ↵
    1. Schmidt L,
    2. Heyes E,
    3. Scheiblecker L,
    4. Eder T,
    5. Volpe G,
    6. Frampton J,
    7. Nerlov C,
    8. Valent P,
    9. Grembecka J,
    10. Grebien F
    (2019) CEBPA-mutated leukemia is sensitive to genetic and pharmacological targeting of the MLL1 complex. Leukemia. doi:10.1038/s41375-019-0382-3.
    OpenUrlCrossRef
  33. ↵
    1. Song G,
    2. Wang L,
    3. Bi K,
    4. Jiang G
    (2015) Regulation of the C/EBPalpha signaling pathway in acute myeloid leukemia (Review). Oncol Rep 33: 2099–2106. doi:10.3892/or.2015.3848
    OpenUrlCrossRef
  34. ↵
    1. Soucie EL,
    2. Weng Z,
    3. Geirsdottir L,
    4. Molawi K,
    5. Maurizio J,
    6. Fenouil R,
    7. Mossadegh-Keller N,
    8. Gimenez G,
    9. VanHille L,
    10. Beniazza M, et al.
    (2016) Lineage-specific enhancers activate self-renewal genes in macrophages and embryonic stem cells. Science 351: aad5510. doi:10.1126/science.aad5510
    OpenUrlAbstract/FREE Full Text
  35. ↵
    1. Tait SW,
    2. Ichim G,
    3. Green DR
    (2014) Die another way: Non-apoptotic mechanisms of cell death. J Cel Sci 127: 2135–2144. doi:10.1242/jcs.093575
    OpenUrlCrossRef
  36. ↵
    1. Taskesen E,
    2. Bullinger L,
    3. Corbacioglu A,
    4. Sanders MA,
    5. Erpelinck CA,
    6. Wouters BJ,
    7. van der Poel-van de Luytgaarde SC,
    8. Damm F,
    9. Krauter J,
    10. Ganser A, et al.
    (2011) Prognostic impact, concurrent genetic mutations, and gene expression features of AML with CEBPA mutations in a cohort of 1182 cytogenetically normal AML patients: Further evidence for CEBPA double mutant AML as a distinctive disease entity. Blood 117: 2469–2475. doi:10.1182/blood-2010-09-307280
    OpenUrlAbstract/FREE Full Text
  37. ↵
    1. Taylor D,
    2. Badiani P,
    3. Weston K
    (1996) A dominant interfering Myb mutant causes apoptosis in T cells. Genes Dev 10: 2732–2744. doi:10.1101/gad.10.21.2732
    OpenUrlAbstract/FREE Full Text
  38. ↵
    1. Tenen DG
    (2003) Disruption of differentiation in human cancer: AML shows the way. Nat Rev Cancer 3: 89–101.doi:10.1038/nrc989
    OpenUrlCrossRefPubMed
  39. ↵
    1. Verhaak RG,
    2. Wouters BJ,
    3. Erpelinck CA,
    4. Abbas S,
    5. Beverloo HB,
    6. Lugthart S,
    7. Lowenberg B,
    8. Delwel R,
    9. Valk PJ
    (2009) Prediction of molecular subtypes in acute myeloid leukemia based on gene expression profiling. Haematologica 94: 131–134. doi:10.3324/haematol.13299
    OpenUrlAbstract/FREE Full Text
  40. ↵
    1. Volpe G,
    2. Clarke M,
    3. Garcia P,
    4. Walton DS,
    5. Vegiopoulos A,
    6. Del Pozzo W,
    7. O'Neill LP,
    8. Frampton J,
    9. Dumon S
    (2015) Regulation of the Flt3 gene in haematopoietic stem and early progenitor cells. PLoS One 10: e0138257. doi:10.1371/journal.pone.0138257
    OpenUrlCrossRef
  41. ↵
    1. Volpe G,
    2. Walton DS,
    3. Del Pozzo W,
    4. Garcia P,
    5. Dasse E,
    6. O'Neill LP,
    7. Griffiths M,
    8. Frampton J,
    9. Dumon S
    (2013) C/EBPalpha and MYB regulate FLT3 expression in AML. Leukemia 27: 1487–1496. doi:10.1038/leu.2013.23
    OpenUrlCrossRef
  42. ↵
    1. Volpe G,
    2. Walton DS,
    3. Grainger DE,
    4. Ward C,
    5. Cauchy P,
    6. Blakemore D,
    7. Coleman DJL,
    8. Cockerill PN,
    9. Garcia P,
    10. Frampton J
    (2017) Prognostic significance of high GFI1 expression in AML of normal karyotype and its association with a FLT3-ITD signature. Sci Rep 7: 11148. doi:10.1038/s41598-017-11718-8
    OpenUrlCrossRef
  43. ↵
    1. Welner RS,
    2. Bararia D,
    3. Amabile G,
    4. Czibere A,
    5. Benoukraf T,
    6. Bach C,
    7. Wansa KD,
    8. Ye M,
    9. Zhang H,
    10. Iino T, et al.
    (2013) C/EBPalpha is required for development of dendritic cell progenitors. Blood 121: 4073–4081. doi:10.1182/blood-2012-10-463448
    OpenUrlAbstract/FREE Full Text
  44. ↵
    1. Wouters BJ,
    2. Lowenberg B,
    3. Erpelinck-Verschueren CA,
    4. van Putten WL,
    5. Valk PJ,
    6. Delwel R
    (2009) Double CEBPA mutations, but not single CEBPA mutations, define a subgroup of acute myeloid leukemia with a distinctive gene expression profile that is uniquely associated with a favorable outcome. Blood 113: 3088–3091. doi:10.1182/blood-2008-09-179895
    OpenUrlAbstract/FREE Full Text
  45. ↵
    1. Ye M,
    2. Zhang H,
    3. Amabile G,
    4. Yang H,
    5. Staber PB,
    6. Zhang P,
    7. Levantini E,
    8. Alberich-Jorda M,
    9. Zhang J,
    10. Kawasaki A, et al.
    (2013) C/EBPa controls acquisition and maintenance of adult haematopoietic stem cell quiescence. Nat Cell Biol 15: 385–394. doi:10.1038/ncb2698
    OpenUrlCrossRefPubMed
  46. ↵
    1. Yue F,
    2. Cheng Y,
    3. Breschi A,
    4. Vierstra J,
    5. Wu W,
    6. Ryba T,
    7. Sandstrom R,
    8. Ma Z,
    9. Davis C,
    10. Pope BD, et al.
    (2014) A comparative encyclopedia of DNA elements in the mouse genome. Nature 515: 355–364. doi:10.1038/nature13992
    OpenUrlCrossRefPubMed
  47. ↵
    1. Zhang P,
    2. Iwasaki-Arai J,
    3. Iwasaki H,
    4. Fenyus ML,
    5. Dayaram T,
    6. Owens BM,
    7. Shigematsu H,
    8. Levantini E,
    9. Huettner CS,
    10. Lekstrom-Himes JA, et al.
    (2004) Enhancement of hematopoietic stem cell repopulating capacity and self-renewal in the absence of the transcription factor C/EBP alpha. Immunity 21: 853–863. doi:10.1016/j.immuni.2004.11.006
    OpenUrlCrossRefPubMed
  48. ↵
    1. Zhao L,
    2. Glazov EA,
    3. Pattabiraman DR,
    4. Al-Owaidi F,
    5. Zhang P,
    6. Brown MA,
    7. Leo PJ,
    8. Gonda TJ
    (2011) Integrated genome-wide chromatin occupancy and expression analyses identify key myeloid pro-differentiation transcription factors repressed by Myb. Nucleic Acids Res 39: 4664–4679. doi:10.1093/nar/gkr024
    OpenUrlCrossRefPubMed
  49. ↵
    1. Zhao L,
    2. Ye P,
    3. Gonda TJ
    (2014) The MYB proto-oncogene suppresses monocytic differentiation of acute myeloid leukemia cells via transcriptional activation of its target gene GFI1. Oncogene 33: 4442–4449. doi:10.1038/onc.2013.419
    OpenUrlCrossRef
  50. ↵
    1. Zuber J,
    2. Rappaport AR,
    3. Luo W,
    4. Wang E,
    5. Chen C,
    6. Vaseva AV,
    7. Shi J,
    8. Weissmueller S,
    9. Fellmann C,
    10. Taylor MJ, et al.
    (2011) An integrated approach to dissecting oncogene addiction implicates a Myb-coordinated self-renewal program as essential for leukemia maintenance. Genes Dev 25: 1628–1640. doi:10.1101/gad.17269211
    OpenUrlAbstract/FREE Full Text
Back to top
Download PDF
Article Alerts
Email Article
Citation Tools
Myb dependence in CEBPA-mutant AML
Giacomo Volpe, Pierre Cauchy, David S Walton, Carl Ward, Daniel Blakemore, Rachael Bayley, Mary L Clarke, Luisa Schmidt, Claus Nerlov, Paloma Garcia, Stéphanie Dumon, Florian Grebien, Jon Frampton
Life Science Alliance Mar 2019, 2 (2) e201800207; DOI: 10.26508/lsa.201800207
del.icio.us logo Digg logo Reddit logo Twitter logo CiteULike logo Facebook logo Google logo Mendeley logo
Issue Cover

In this Issue

Volume 2, No. 2
April 2019
Advertisement

Jump to section

Subjects

Related Articles

  • No related articles found.

Cited By...

  • No citing articles found.

More in this TOC Section

Similar Articles