Open Access

CEBPA-regulated lncRNAs, new players in the study of acute myeloid leukemia

  • James M Hughes1,
  • Beatrice Salvatori2,
  • Federico M Giorgi2,
  • Irene Bozzoni1, 3, 4 and
  • Alessandro Fatica1Email author
Journal of Hematology & Oncology20147:69

DOI: 10.1186/s13045-014-0069-1

Received: 1 September 2014

Accepted: 16 September 2014

Published: 25 September 2014

Abstract

CCAAT/enhancer-binding protein-α (CEBPA) is a critical regulator of myeloid differentiation. Disruption of CEBPA function contributes to the development of acute myeloid leukemia (AML). CEBPA regulates a large number of protein coding genes of which several were shown to contribute to CEBPA function. In this study, we expand the analysis of CEBPA transcriptional targets to the newly identified class of long non-coding RNAs. We show that lncRNAs are a main component of the transcriptional program driven by C/EBPα and that many of these are also induced during granulocytic differentiation of AML cell lines supporting their relevance in proliferation arrest and differentiation.

Keywords

lncRNAs CEBPA AML

To the Editor

LncRNAs participate in multiple networks controlling cell differentiation and development [1], with their expression already associated with cancer and several disorders [2]. To what degree C/EBPα regulates the expression of lncRNAs is still largely unknown.

To investigate the effect of C/EBPα on the expression of lncRNAs we utilized the K562 AML cell line carrying a stable and Tet-on inducible CEBPA allele (Additional file 1 and Additional file 2: Figure S1). K562 cells lack endogenous C/EBPα and restoration of its expression induces proliferation arrest and granulocytic differentiation [3] (Additional file 1 and Additional file 2: Figure S1). Based on the expression of known C/EBPα transcriptional targets, we selected RNA extracted from 48 hours of induction (K562-C/EBPα) together with RNA extracted from control-induced cells (K562-CTR). Gene expression profiling was performed using the Agilent Whole Human Genome Oligo 8x60K v2 Microarrays from 4 biological replicates for each sample (Figure 1A). We identified 4605 mRNAs (2643 induced and 1962 repressed) and 930 lncRNAs (600 induced and 330 repressed) with significant differential expression (fold change ≥2 and p-value ≤ 0.05) between C/EBPα- and CTR- induced cells (Figure 1B and Additional file 1, Additional files 3 and 4: Tables S1 and S2). Appropriate expression patterns of many known coding transcriptional targets of C/EBPα confirmed the reliability of our gene expression analysis (Additional file 1 and Additional file 2: Figure S1). Gene set enrichment and Gene ontology analysis confirmed significant enrichment of known C/EBPα targets [4] (NES = 7.97, p = 1.65×10-15), coupled with down-regulation of cell cycle genes and upregulation of granulocytic differentiation pathways (Figure 2). Notably, we found the E2F1 motif to be negatively enriched in the promoters of C/EBPα repressed genes (NES = -7.18, p = 7.06×10-13), confirming the known role of C/EBPα in repressing E2F1 activity (Additional file 5: Figure S2). Expression of differentially induced lncRNAs was further validated by qRT-PCR in K562 -C/EBPα and -CTR cells (Figure 2A). When applicable, official lncRNA reference names were utilized. Otherwise, we refer to as lnc-CUs (lncRNA-C/EBPα-up-regulated) and lnc-DCs (lncRNA-C/EBPα down-regulated) for induced and repressed lncRNAs, respectively (Additional file 1 and Additional file 6: Table S3).
Figure 1

LncRNAs are a component of the C/EBPα-regulated transcriptional network in AML. (A) Experimental pipeline for identification of C/EBPα-regulated lncRNAs. See text and Additional file 1 for details (B) Heat maps of median centered log2 fold change values from microarray analysis of mRNAs and lncRNAs with significant differential expression (fold change ≥ 2 and adjusted p-value ≤ 0.05). Exact values are provided in Additional files 3 and 4: Table S1 and Table S2. Presence of CEBPB or CEBPD binding sites is indicated on the side of each heat maps.

Figure 2

Validation of selected lncRNAs. (A) Quantitative real time RT-PCR analysis of selected lncRNAs in K562-C/EBPα treated for 48 hours with doxycycline. Values were compared to K562-CTR cells treated for the same time and normalized with HPRT mRNA. The histograms represent the log2 fold change of the relative expression ± SEM from three replicates. TCONS names and chromosomal positions are provided as Additional file 6: Table S3. (B) Quantitative real time RT-PCR analysis of selected lncRNAs in NB4 cells treated for 72 hours with all-trans-retinoic acid (ATRA). Values were compared to untreated NB4 cells and normalized with HPRT mRNA. The histograms represent the log2 fold change of the relative expression ± SEM from three replicates.

In order to annotate the presence of putative C/EBPα binding sites in the promoter of differentially expressed genes, we used previously generated ChIP data sets for CEBPB (C/EBPα) and CEBPD (C/EBPα) in K562 cells [5], which exhibit identical DNA-binding specificities with C/EBPα [6]. We found several coding and non-coding differentially expressed genes bound by either CEBPB or CEBPD in their putative promoter region within a distance of -5 kb from the TSS (Figure 1B and Additional file 1, Additional file 7: Figure S3, Additional file 8: Table S4, Additional file 9: Table S5, Additional file 10: Table S6 and Additional file 11: Table S7).

Different AML cell lines are widely used to study the block of differentiation in AML because they can be differentiated in mature and functional myeloid cells by treatment with specific agents. Thus, we analysed the expression of selected lncRNAs in NB4 cells, which are able to undergo granulocytic differentiation by treatment with all-trans retinoic acid (ATRA) [7]. Notably, the majority of validated C/EBPα-induced lncRNAs in K562 are also significantly upregulated by ATRA in NB4 (21 out of 26), suggesting that they may play a role in the differentiation process (Figure 2B). Nevertheless, upon validated lncRNAs repressed by C/EBPα treatment in K562, 6 out of 8 showed opposite trend while 2 were not significantly expressed in NB4 (data not shown). This behaviour still remains to be explained and extended to the study of more lncRNAs in NB4 cell line: we speculate it could be due to silencing of diverse cellular settings between K562 and NB4 cell lines.

In summary, this study shows that lncRNAs are a main component of the transcriptional program driven by C/EBPα. We identified more than 900 lncRNAs regulated by C/EBPα in K562. We confirmed that the majority of these are also induced during granulocytic differentiation of AML cell lines supporting their relevance in proliferation arrest and differentiation. How many of the lncRNAs identified in this study are directly involved in regulating differentiation programmes of AML is an interesting question that warrants further investigations.

Moreover, regardless of function, this work indicates that changes in lncRNAs expression might also have diagnostic applications in AML with CEBPA mutations.

Authors’ contribution

JMH and BS conducted experiments, participated in research design and interpretation of data. FMG performed data analysis. IB participated in research design and provided financial support. AF designed research, wrote the manuscript, and provided financial support. All authors read and approved the final manuscript.

Additional files

Declarations

Acknowledgments

The authors would like to thank Prof K. Nerlov for CEBPA plasmid, Dr A. Rosa and Dr A. Brivanlou for the ePiggyBac inducible transposon system, M. Arceci and M. Marchioni for technical assistance. This work was supported by FP7-PEOPLE-2011-ITN Project HemID (289611), Italian Epigenomics Flagship Project (EPIGEN) and “Research Projects of National Interest” (PRIN).

Authors’ Affiliations

(1)
Department of Biology and Biotechnology C. Darwin, Sapienza University of Rome
(2)
Department of Systems Biology, Herbert Irving Comprehensive Cancer Center, Columbia University Medical Center
(3)
Center for Life Nano Science@Sapienza, Istituto Italiano di Tecnologia
(4)
Institute Pasteur Fondazione Cenci-Bolognetti, Sapienza University of Rome

References

  1. Fatica A, Bozzoni I: Long non-coding RNAs: new players in cell differentiation and development. Nat Rev Genet. 2014, 15: 7-21. 10.1038/nrg3606.View ArticlePubMedGoogle Scholar
  2. Batista PJ, Chang HY: Long noncoding RNAs: cellular address codes in development and disease. Cell. 2013, 152: 1298-1307. 10.1016/j.cell.2013.02.012.PubMed CentralView ArticlePubMedGoogle Scholar
  3. Tavor S, Park DJ, Gery S, Vuong PT, Gombart AF, Koeffler HP: Restoration of C/EBPalpha expression in a BCR-ABL + cell line induces terminal granulocytic differentiation. J Biol Chem. 2003, 278: 52651-52659. 10.1074/jbc.M307077200.View ArticlePubMedGoogle Scholar
  4. Subramanian A, Tamayo P, Mootha VK, Mukherjee S, Ebert BL, Gillette MA, Paulovich A, Pomeroy SL, Golub TR, Lander ES, Mesirov JP: Gene set enrichment analysis: A knowledge-based approach for interpreting genome-wide expression profiles. PNAS. 2005, 102: 15545-15550. 10.1073/pnas.0506580102.PubMed CentralView ArticlePubMedGoogle Scholar
  5. Johnson DS, Mortazavi A, Myers RM, Wold B: Genome-wide mapping of in vivo protein-DNA interactions. Science. 2007, 316 (5830): 1497-1502. 10.1126/science.1141319.View ArticlePubMedGoogle Scholar
  6. Tsukada J, Yoshida Y, Kominato Y, Auron PE: The CCAAT/enhancer (C/EBP) family of basic-leucine zipper (bZIP) transcription factors is a multifaceted highly-regulated system for gene regulation. Cytokine. 2011, 54: 6-19. 10.1016/j.cyto.2010.12.019.View ArticlePubMedGoogle Scholar
  7. Roussel MJ, Lanotte M: Maturation sensitive and resistant t(15;17) NB4 cell lines as tools for APL physiopathology: nomenclature of cells and repertory of their known genetic alterations and phenotypes. Oncogene. 2001, 20: 7287-7291. 10.1038/sj.onc.1204863.View ArticlePubMedGoogle Scholar

Copyright

© Hughes et al.; licensee BioMed Central Ltd. 2014

This article is published under license to BioMed Central Ltd. This is an Open Access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/4.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly credited. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated.

Advertisement