| HOME | HELP | FEEDBACK | SUBSCRIPTIONS | ARCHIVE | SEARCH | TABLE OF CONTENTS |
Molecular Biology, Pathobiology, and Genetics |
1 Hematopathology Unit, 2 Hematology and Pathology Departments, and 3 Genomics Unit, Hospital Clínic, IDIBAPS, University of Barcelona; and 4 Hematology Department, Hospital de la Santa Creu i Sant Pau, Barcelona, Spain
Requests for reprints: Jordi Esteve, Hematology Department, Hospital Clínic, Villarroel 170, 08036 Barcelona, Spain. Phone: 34-93-227-54-28; Fax: 34-93-227-54 84; E-mail: jesteve{at}clinic.ub.es.
| Abstract |
|---|
|
|
|---|
, RUNX1-RUNX1T1, and CBFß-MYH11 rearrangements. The relative expression of 46 genes, selected according to their differential expression in the high-density array study, was analyzed by low-density arrays in an additional series of 40 patients, which included 7 MYST3-CREBBP AML cases. Thus, genes such as prolactin (PRL) and proto-oncogene RET were confirmed to be specifically overexpressed in MYST3-CREBBP samples whereas genes such as CCND2, STAT5A, and STAT5B were differentially underexpressed in this AML category. Interestingly, MYST3-CREBBP AML exhibited a characteristic pattern of HOX expression, with up-regulation of HOXA9, HOXA10, and cofactor MEIS1 and marked down-regulation of other homeobox genes. This profile, with overexpression of FLT3, HOXA9, MEIS1, AKR7A2, CHD3, and APBA2, partially resembles that of AML with MLL rearrangement. In summary, this study shows the distinctive gene expression profile of MYST3-CREBBP AML, with overexpression of RET and PRL and a specific pattern of HOX gene expression. (Cancer Res 2006; 66(14): 6947-54) | Introduction |
|---|
|
|
|---|
, RUNX1-RUNX1T1 (AML1-ETO), and CBFß-MYH11, found in
25% of de novo AML cases, constitute abnormal transcriptional factors causing a disturbed program of myeloid differentiation. Moreover, each of these translocations defines a specific leukemia subtype associated with a favorable prognosis (1, 2). In this setting, translocation t(8;16)(p11;p13) is an infrequent recurrent chromosomal abnormality found in both de novo and therapy-related AML cases after treatment with topoisomerase II inhibitors (36). These patients present with specific clinical and biological features, such as a blast population with a myelomonocytic stage of differentiation, frequent extramedullary involvement, severe coagulation disorder, and a poor outcome (3, 5, 6). At the molecular level, translocation t(8;16) fuses MYST histone acetyltransferase (monocytic leukemia)-3 (MYST3; formerly named MOZ) and CREB binding protein (Rubinstein-Taybi syndrome; CREBBP, or CBP) genes, both encoding proteins with histone acetyltransferase activity (4, 79). MYST3 has been shown to modulate gene transcription through activation of the transcription factor complex RUNX1 (8, 10). Moreover, the protein complex MYST3-RUNX1 has been found to increase during normal monocytic differentiation. In its turn, CREBBP protein also regulates transcription by means of histone acetyltransferase activity and by binding to several proteins with key cell cycle functions, such as p53 and nuclear factor
B (8, 9). Therefore, an inhibition of RUNX1-mediated transcription by MYST3-CREBBP fusion protein has been hypothesized to be the main mechanism of leukemogenesis in this AML variety (10). However, the precise pathways disrupted by this chimerical protein are mostly unknown. Analysis of gene expression profile might contribute to refine the classification of AML based on biological grounds and to assign the prognostic risk of a given subtype more accurately (1114). Furthermore, genomic analysis of AML might provide a deeper insight into the underlying disease mechanisms.
In this study, we have examined the gene expression profile of AML with the MYST3-CREBBP fusion gene to determine the specific signature of this leukemia compared with other well-defined AML subtypes and to define possible molecular pathways involved in the pathogenesis of this leukemia.
| Materials and Methods |
|---|
|
|
|---|
; n = 3); AML with t(8;21)(q22;q22) (RUNX1-RUNX1T1; n = 3); AML with inv(16)/t(16;16) (CBFß-MYH11; n = 3); AML with t(9;11)(p22;q23) (MLLT3-MLL; n = 1); acute monocytic leukemias (n = 8); and two cases of AML with multilineage dysplasia. In 10 of these 23 cases, an internal tandem duplication of the FLT3 gene was detected whereas mutations of nucleophosmin (NPM) gene were found in 6 cases (Supplementary Table A). Two cases with MYST3-CREBBP rearrangement and one monocytic-differentiated AML (cases 1, 3, and 18) followed exposure to topoisomerase II inhibitors. Cases 1 and 18 also presented amplification of MLL gene.
To confirm the findings of the previous global gene expression profile study, a second subset of 40 AML patients (subset B, Supplementary Table B) was studied using TaqMan low-density arrays (Applied Biosytems, Foster City, CA; see below). This subset of patients included the three MYST3-CREBBP AML cases previously studied by high-density array and four additional MYST3-CREBBP samples with no appropriate material for the genome-wide assay. In addition, an independent set of 33 AML samples was included in this study. These samples corresponded to AML with well-characterized rearrangements (PML-RAR
, n = 3; RUNX1-RUNX1T1, n = 3; CBFß-MYH11, n = 3; MLL-rearranged AML, n = 9) and normal karyotype AML (n = 15). Mutations of FLT3-ITD, NPM, and MLL abnormalities are detailed in Supplementary Table B. None of the MYST3-CREBBP samples analyzed harbored either FLT3-ITD or NPM mutations. Four patients (cases 2, 3, 7, and 18) presented with a therapy-related AML. The clinical and laboratory data of the seven MYST3-CREBBP AML cases were previously reported (16). All peripheral blood and bone marrow samples were obtained with informed consent according to the guidelines of the Ethical Committee of the participating Institutions.
RNA extraction and cDNA synthesis. Total RNA was isolated from peripheral blood and bone marrow by standard methods (17). In one case, we obtained RNA from a paraffin-embedded tissue using a phenol-extraction method. In all cases, integrity of RNA was examined with Agilent 2100 Bioanalyzer (Agilent Technologies, Palo Alto, CA). One microgram of RNA was reverse transcribed to cDNA using random primers with the High Capacity cDNA Archive Kit (Applied Biosystems).
Molecular analyses. The detection of transcript type I of MYST3-CREBBP rearrangement was analyzed by reverse transcription-PCR (RT-PCR) with primers MOZ3558F and CPB335R, as previously described (16). The presence of PML-RAR
, RUNX1-RUNX1T1, and CBFß-MYH11 rearrangement was analyzed by RT-PCR following published conditions (18). The analysis of FLT3-ITD was done using primers and conditions previously published (19). The presence of mutations in exon 12 of NPM gene was studied by amplification of genomic DNA or cDNA with a fluorescently labeled forward primer and subsequent analysis of the PCR product in an Automatic sequencer (Abi Prism 310) using the Genescan software as described (20). MLL rearrangement was studied by Southern blot analysis using probe B859 after digestion with HindIII and BamHI enzymes (21). Cases without available DNA, with a noninformative karyotype, or showing abnormalities at 11q23 region by conventional cytogenetics were studied by fluorescence in situ hybridization with the LSI MLL Dual Color Probe (Vysis) as described (22). Partial tandem duplication was specifically analyzed by RT-PCR as previously described (23).
High-density array study: RNA purification, labeling, and hybridization. In this part of the study, 23 AML samples with confirmed high quality RNA by Agilent 2100 Bioanalyzer analysis and sufficient amount of RNA (5 µg) were included. Amplified biotinylated complementary RNA was produced with an in vitro transcription labeling reaction and was subsequently hybridized to Affymetrix HU133A oligonucleotide arrays following the Affymetrix protocol for high-density arrays (details are provided in supplementary material). Scans were carried out on Agilent G2500A GeneArray scanner (Agilent Technologies, Waldbronn, Germany) and the fluorescence intensities of scanned arrays were analyzed with the Affymetrix GeneChip software.
Statistical analysis. Affymetrix Microarray Suite Software version 5.0 (MAS5.0) was used for the quantification of the expression level of target genes. HU133A microarray raw expression intensities were scaled to a target intensity of 200 units. To exclude genes with minimal variation across samples, only genes with a mean (SD) of normalized values between 0.7 and 10 were filtered. Thereafter, we selected those genes with an expression level of
20 in
25% of samples. The 2,683 resulting genes were studied by means of unsupervised two-dimensional cluster analysis with dChip v1.3 software using the default clustering algorithm defined as 1 r, where r is the Pearson correlation coefficient between standardized expression values (make 0 and SD 1) and the centroid linkage method. To identify those genes with significant differences in their expression level among different AML categories, we used a random-variance F test as described in BRB ArrayTools software (BRB ArrayTools developed by Dr. Richard Simon and Amy Peng Lam)5 with all probe sets, but assigning an arbitrary value of 10 to genes with an expression level below 10 units. A significance level of 0.001 was chosen to reduce the number of false positive results. For each one of the differentially expressed genes, a ratio between the mean expression value in MYST3-CREBBP samples and the mean expression in every AML category was calculated. Additionally, a t test comparing the gene expression in MYST3-CREBBP and that in the remaining samples as a whole, with a significance level of 0.001, was used as another method of detection of differentially expressed genes in MYST3-CREBBP samples.
Quantitative real-time RT-PCR. A selection of 46 genes was subsequently studied by real-time PCR using TaqMan Low Density Arrays (Applied Biosystems) in an additional series of 40 AML patients (subset B, Supplementary Table B). These genes were selected on the basis of their differential expression in MYST3-CREBBP according to high-density array analysis or their oncogenic potential in leukemia. The list of 46 genes is provided in Supplementary Table C. Briefly, cDNA was obtained from these cases, loaded onto the low-density arrays, and amplified using standard conditions in an Abi Prism 7900HT Sequence Detection System (Applied Biosystems). All samples were tested in duplicate and the average value between replicates was taken as the specific level of expression of a given gene. To quantify the relative expression of each gene, the Ct values were normalized for endogenous reference (
Ct = Cttarget Ctß2-microglobulin) and compared with a calibrator using the 
Ct method. As calibrator, the average Ct value of each gene in all samples grouped together was taken. Comparison of the relative expression of the 46 genes in MYST3-CREBBP AML with that in the remaining AML samples was done using a t test comparing two groups, with a significance level of 0.05. In addition, an ANOVA test was also used to compare the relative expression of these genes after defining different AML categories (significance level, 0.05).
| Results |
|---|
|
|
|---|
(cluster 4), and, on the other hand, AML with multilineage dysplasia (cluster 5)]. In contrast, the remaining samples (no. 13-21), defined by their monocytic differentiation and the absence of the above-mentioned fusion proteins, were distributed among the different clusters of the array forming a heterogeneous group. One of the samples harboring a CBFß-MYH11 rearrangement, with a minimally differentiated phenotype (M0 subtype), segregated from two other cases with the same molecular alteration but presenting with a myelomonocytic phenotype.
|
As an additional method to study possible shared patterns of gene expression between MYST3-CREBBP and other AML, we selected those genes with high or low expression in MYST3-CREBBP leukemias and a similar pattern of expression in one of the other AML categories (Supplementary Tables E and F; Supplementary Fig. B). This analysis revealed that MYST3-CREBBP leukemias had 33 genes overexpressed in common with AML of monocytic lineage (i.e., AKR7A2, CHD3, and AK2) and 19 with PML-RAR
, but only a minority of genes in common with other AML subtypes.
Finally, when MYST3-CREBBP cases were compared with the other samples as a group using a t test, 237 genes (53 overexpressed and 184 underexpressed) showed a differential expression (Supplementary Tables G and H). Using this analysis, a high expression of the aforementioned gene RET could be observed. Additionally, an overexpression of genes with a role in transcription (SATB1), genes involved in apoptosis (OPTN), and CEBP
, a crucial gene in the myeloid differentiation process, was observed.
Low-density Array Analysis
Supervised analysis. Forty-six genes were selected according to their differential expression in MYST3-CREBBP samples in the high-density array and/or their relevant oncogenic role in leukemia. The relative expression of these genes was analyzed in 40 AML patients (subset B, Supplementary Table B). To define the group of genes characteristic of MYST3-CREBBP subtype, a t test with two groups was used to compare the gene expression in MYST3-CREBBP samples and the remaining AML patients. Twenty-two genes were significantly overexpressed in MYST3-CREBBP samples whereas three genes were underexpressed in this group of leukemias (Table 1
). Subsequently, an ANOVA test was used to compare the gene relative expression among the seven AML categories defined by their underlying molecular abnormality (A, MYST3-CREBBP; B, PML-RAR
; C, RUNX1-RUNX1T1; D, CBFß-MYH11; E, MLL-translocated samples; F, MLL partial tandem duplication; and G, miscellaneous). The ANOVA test yielded five different gene expression profiles. First (profile 1, Fig. 2A
), a group of nine genes (PRL, C20orf103, RET, GGA2, ICSBP1, ITGA7, DAP, IRAK1 and PPARG) were found to be highly expressed in MYST3-CREBBP AML with absent or low expression in the remaining samples. Genes such as FLT3, HOXA9, MEIS1, AKR7A2, CHD3, and APBA2 composed a second pattern of expression (profile 2, Fig. 2B), with high expression level in both MYST3-CREBBP and AML with monocytic phenotype. In the third profile observed (profile 3, Fig. 2C), some genes had a high expression level in MYST3-CREBBP AML and were variably expressed in other AML categories. This group included SATB1, LMO2, CEBP
, SURF1, HIST1H2A, DDB2, and LGALS3 genes.
|
|
|
Expression of Homeobox Genes
Given the high expression of several homeobox (HOX) genes observed in MYST3-CREBBP, we did an unsupervised analysis on the whole series of patients focused on the relative level of expression of the HOX genes present in the array. Several patterns of HOX genes expression were seen in different AML categories (Fig. 4A
). Thus, overexpression of 5'-HOX genes (HOXA9, HOXB9, and HOXA10) was characteristically found in all MYST3-CREBBP samples and most cases of monocytic AML. Other HOX genes (HOXA2-A7 and HOXB2-B7) showed a high expression in the group of monocytic leukemias whereas their expression was low in MYST3-CREBBP AML. In contrast, HOX genes were underexpressed in good-risk cytogenetic AML.
|
| Discussion |
|---|
|
|
|---|
The combination of several methods of comparative analysis allowed the identification of groups of genes with a differential expression in distinct AML categories. First, a subset of genes seemed to be highly characteristic of MYST3-CREBBP AML, being up-regulated in these cases and showing a low or absent expression in the remaining AML categories. Thus, genes such as PRL, C20orf103, RET, GGA2, ICSBP1, ITGA7, DAP, IRAK1, and PPARG were overexpressed almost exclusively in MYST3-CREBBP cases. Among those genes, PRL and RET have been occasionally reported to be involved in leukemogenesis (2433). In this regard, an increased expression of prolactin protein in blast populations has been observed in anecdotal cases of monocytic-lineage leukemia as well as in the eosinophilic cell line Eol-1 (2529). In a recent study, prolactin expression in Eol-1 cells was shown to involve different signaling pathways induced by cyclic AMP whereas inhibition of Jak-STAT5 pathway resulted in up-regulation of prolactin (28). Interestingly, STAT5A and STAT5B genes were found to be significantly underexpressed in our MYST3-CREBBP AML samples, suggesting a negative regulatory effect between prolactin and STAT5 proteins. On its turn, RET gene is a proto-oncogene that encodes a tyrosine kinase receptor expressed during normal myelomonocytic differentiation and, accordingly, has been found to be predominantly expressed in AML of monocytic phenotype (3032). Moreover, RET mRNA levels are typically low among immature CD34+ hemopoietic progenitors whereas overexpression of this protein has been associated with coexpression of adhesion molecules such as CD56 (30). These observations resemble the characteristic phenotype of MYST3-CREBBP AML, defined by common CD34 negativity, high expression of monocytic antigens, and frequent coexpression of CD56 and NG2 (16). Of note, RET has recently been reported as one of the most characteristic genes in one of the 16 AML clusters defined according to its genomic profile in a recent study by Valk et al. (14). Although most of the cases forming that cluster were monocytic-differentiated leukemias (14), modulation of RET does not seem to be merely attributable to a monocytic differentiation process because, in the present study, the expression level of RET was significantly higher in MYST3-CREBBP cases than in other monocytic-lineage leukemia samples. Somatic and germ-line RET mutations leading to gene activation are responsible for several human diseases, including multiple endocrine neoplasia types 2A and 2B and papillary thyroid carcinomas (24). Nevertheless, mutations of RET as a mechanism of overexpression were discarded in the present MYST3-CREBBP series, in accordance with a previous study analyzing diverse AML subtypes (33).
One of the most striking findings of this study was the similarities observed between MYST3-CREBBP and MLL-rearranged leukemias. MYST3-CREBBP cases presented high levels of homeobox genes (HOXA9 and HOXA10), their cofactor MEIS1, and the receptor with tyrosine-kinase activity FLT3, all of them typically up-regulated in MLL leukemias (3436). HOX genes are transcription factors required for a proper hematopoietic development and constitute downstream targets of MLL protein (37). In this regard, MLL-rearranged leukemias are typically characterized by an impaired pattern of HOX expression, and HOXA9, HOXA10, and the cofactor MEIS1 are up-regulated in virtually all lymphoid, myeloid, and biphenotypic lineage MLL-rearranged leukemia subtypes (38). In addition to the unexpected HOX overexpression in MYST3-CREBBP cases, the results of the present study confirmed previous findings on the pattern of HOX expression in different AML subtypes (3941). Thus, whereas an overall low expression of HOX genes was observed in cases of AML with favorable cytogenetics, high levels of the members of HOX family HOXA9 and HOXA10 were seen in MYST3-CREBBP and MLL-rearranged leukemias. In contrast, other HOX genes analyzed, such as HOXA3, HOXA7, and HOXB5, were expressed at different levels in MLL-translocated cases, MLL with partial tandem duplication, and other monocytic leukemias, but were not expressed in MYST3-CREBBP AML cases.
Similarities between MYST3-CREBBP and MLL leukemias comprised other genes. Thus, overexpression of AKR7A2, PBX3, NICAL, and IRAK1B genes, observed in MYST3-CREBBP subtype, was also found in AML with MLL rearrangement, as reported in a recent study by Kohlmann et al. (36). Moreover, a coincidental expression of several genes (RET, C20orf103, GGA2, GAGED2, AKR7A2, and AK2) was observed between our MYST3-CREBBP samples and the above-mentioned cluster no. 16 of the study of Valk et al. (14). Of note, 5 of the 11 patients included in this cluster had 11q23 abnormalities.
As a common mechanism of leukemogenesis, chimeric protein PML-RAR
in acute promyelocytic leukemia and RUNX1-RUNX1T1 and CBFß-MYH11, associated to core binding factor leukemias, respectively, induce a constitutive transcriptional repression, leading to a blockage of the normal myeloid differentiation program (42). This contrasts with the presumed function of fusion gene products derived from the rearrangement of MLL with different partners, thought to produce a constitutive transcriptional activation through a gain-of-function mechanism resulting in an inappropriate expression of target genes such as HOX (42). In this context, the gene expression signature of MYST3-CREBBP rearrangement obtained in this study seemed to be similar to that of MLL-rearranged leukemia and markedly different from those of CBF-AML and acute promyelocytic leukemia. Therefore, the leukemogenic effect of MYST3-CREBBP fusion gene could rely on a deregulated modulation of downstream targets, resembling that of MLL chimeras, probably due to impaired histone acetyltransferase activity of the proteins involved in this translocation. Moreover, the adverse prognosis classically associated with this entity differs from that of CBF-AML and acute promyelocytic leukemia and is similar to that of MLL-rearranged AML.
With regard to the underexpressed genes, the CREBBP expression was decreased in MYST3-CREBBP samples, although not reaching a significant difference, thus suggesting a negative regulation of the chimerical protein MYST3-CREBBP over native CREBBP transcript. Additionally, down-regulation of WT1 was observed in MYST3-CREBBP, as well as in RUNX1-RUNX1T1, as compared with the other AML categories. WT1 gene has been shown to be overexpressed in >90% of AML cases and it has recently been proposed as a molecular marker for minimal residual disease studies in AML (43). However, the low levels of WT1 in MYST3-CREBBP leukemias observed in our study would hamper this strategy in the follow-up of the minimal residual disease in this group of leukemia.
In summary, the double strategy followed, based on the selection of a group of genes according to their differential expression in the high-density array assay and further analysis by real-time PCR in an additional set of patients, allowed the assessment of gene profile in a group of seven MYST3-CREBBP patients. Of note, a distinctive gene expression signature of MYST3-CREBBP leukemias was observed, characterized by the overexpression of homeobox genes HOXA9, HOXA10, and their cofactor MEIS1; the up-regulation of the oncogenes RET and PRL; and the decreased expression of genes such as CCND2, STAT5, and WT1. This profile harbors some similarities with that of MLL-rearranged leukemias, thus suggesting a partially common leukemogenic pathway.
| Acknowledgments |
|---|
The costs of publication of this article were defrayed in part by the payment of page charges. This article must therefore be hereby marked advertisement in accordance with 18 U.S.C. Section 1734 solely to indicate this fact.
We thank all the staff of the Hemopathology Unit and Montse Sánchez from the Genomics Unit for their excellent work.
| Footnotes |
|---|
5 http://linus.nci.nih.gov/BRB-ArrayTools.html. ![]()
Received 12/22/05. Revised 4/28/06. Accepted 5/16/06.
| References |
|---|
|
|
|---|
This article has been cited by other articles:
![]() |
K. I. Mills, A. Kohlmann, P. M. Williams, L. Wieczorek, W.-m. Liu, R. Li, W. Wei, D. T. Bowen, H. Loeffler, J. M. Hernandez, et al. Microarray-based classifiers and prognosis models identify subgroups with distinct clinical outcomes and high risk of AML transformation of myelodysplastic syndrome Blood, July 30, 2009; 114(5): 1063 - 1072. [Abstract] [Full Text] [PDF] |
||||
![]() |
U. Bacher, A. Kohlmann, and T. Haferlach Perspectives of gene expression profiling for diagnosis and therapy in haematological malignancies Brief Funct Genomic Proteomic, May 27, 2009; (2009) elp011v1. [Abstract] [Full Text] [PDF] |
||||
![]() |
B. Xu, D. Geerts, K. Qian, H. Zhang, and G. Zhu Myeloid ecotropic viral integration site 1 (MEIS) 1 involvement in embryonic implantation Hum. Reprod., June 1, 2008; 23(6): 1394 - 1406. [Abstract] [Full Text] [PDF] |
||||
![]() |
X.-L. Liu, J.-Y. Yuan, J.-W. Zhang, X.-H. Zhang, and R.-X. Wang Differential gene expression in human hematopoietic stem cells specified toward erythroid, megakaryocytic, and granulocytic lineage J. Leukoc. Biol., October 1, 2007; 82(4): 986 - 1002. [Abstract] [Full Text] [PDF] |
||||
![]() |
E. Morgado, S. Albouhair, and C. Lavau Flt3 is dispensable to the Hoxa9/Meis1 leukemogenic cooperation Blood, May 1, 2007; 109(9): 4020 - 4022. [Abstract] [Full Text] [PDF] |
||||
| ||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||
| HOME | HELP | FEEDBACK | SUBSCRIPTIONS | ARCHIVE | SEARCH | TABLE OF CONTENTS |
| Cancer Research | Clinical Cancer Research |
| Cancer Epidemiology Biomarkers & Prevention | Molecular Cancer Therapeutics |
| Molecular Cancer Research | Cancer Prevention Research |
| Cancer Prevention Journals Portal | Cancer Reviews Online |
| Annual Meeting Education Book | Meeting Abstracts Online |