Acute Myeloid Leukemia |
1 Division of Hematology and Oncology, Cedars-Sinai Medical Center, UCLA School of Medicine, Los Angeles, CA, USA
2 Department of Hematology and Oncology and
3 Department of Cell Therapy and Transplantation Medicine and the 21st Century COE Program, Graduate School of Medicine, University of Tokyo, Tokyo, Japan
4 Core Research for Evolutional Science and Technology, Japan Science and Technology Agency, Tokyo, Japan
5 Department of Medical Informatics and Biomathematics, University of Munster, Munster, Germany
6 MLL Munich Leukemia Laboratory, Munich, Germany
Correspondence: Tadayuki Akagi, Ph.D, Division of Hematology and Oncology, Cedars-Sinai Medical Center, UCLA School of Medicine 8700 Beverly Blvd, Los Angeles, CA90048, USA. E-mail:tadayuki{at}staff.kanazawa-u.ac.jp
|
|
|---|
Design and Methods: Thirty-eight normal karyotype acute myeloid leukemia/myelodysplastic syndrome samples were analyzed with high-density single-nucleotide polymorphism microarray using a new algorithm: allele-specific copy-number analysis using anonymous references (AsCNAR). Expression of mRNA in these samples was determined by mRNA microarray analysis.
Results: Eighteen samples (49%) showed either one or more genomic abnormalities including duplication, deletion and copy-number neutral loss of heterozygosity. Importantly, 12 patients (32%) had copy-number neutral loss of heterozygosity, causing duplication of either mutant FLT3 (2 cases), JAK2 (1 case) or AML1/RUNX1 (1 case); and each had loss of the normal allele. Nine patients (24%) had small copy-number changes (< 10 Mb) including deletions of NF1, ETV6/TEL, CDKN2A and CDKN2B. Interestingly, mRNA microarray analysis showed a relationship between chromosomal changes and mRNA expression levels: loss or gain of chromosomes led, respectively, to either a decrease or increase of mRNA expression of genes in the region.
Conclusions: This study suggests that at least one half of cases of normal karyotype acute myeloid leukemia/myelodysplastic syndrome have readily identifiable genomic abnormalities, as found by our analysis; the high frequency of copy-number neutral loss of heterozygosity is especially notable.
Key words: normal karyotype acute myeloid leukemia/myelodysplastic syndrome, SNP-chip, CNN-LOH.
|
|
|---|
Besides chromosomal abnormalities, the leukemic cells can have a variety of mutations involving individual genes. Activating mutations of the receptor tyrosine kinase, FMS-like tyrosine kinase 3 (FLT3) occur in about 30% AML patients; two major mutant forms occur: an internal tandem duplication (ITD) or a point mutation in the tyrosine kinase domain (TKD).15 Activating mutations at codon 12, 13 or 61 of either the NRAS or KRAS occur in 25% and 15% of AML patients, respectively.1,16 About 10–15% of AML samples have inactivating mutations of C/EBP
whose wild-type function is to enhance differentiation.17,18 Nucleophosmin1 (NPM1) is mutated in 50–60% of AML samples with normal karyotype.13,19 This protein has an important role in ribosome biogenesis, including nuclear export of ribosomal proteins. Mutant NPM1 has an aberrant nuclear export signal and remains localized in the cytoplasm.20
Single-nucleotide polymorphism microarray (SNP-chip) analysis is a new technique to examine the genome including any copy-number changes and loss of heterozygosity (LOH).21–23 Importantly, SNP-chip analysis can reveal cryptic abnormalities such as a small copy-number changes (< 10 Mb) or copy-number neutral loss of heterozygosity [CNN-LOH, also called uniparental disomy (UPD)] that cannot be detected by karyotype analysis. In addition, comparative genomic hybridization cannot detect CNN-LOH. SNP-chip analysis has been used in chronic lymphocytic leukemia,24,25 childhood acute lymphoblastic leukemia,26,27 juvenile myelomonocytic leukemia,28 follicular lymphoma,29 multiple myeloma,30 and AML.31,32,50–54
In the present study, we identified hidden abnormalities and novel disease-related genomic regions using 250 K SNP-chip analysis in samples from patients with normal karyotype AML/myelodysplastic syndrome (MDS). The use of CNAG (copy-number analysis for Affymetrix GeneChips) program21 and a new algorithm AsCNAR (allele-specific copy-number analysis using anonymous references)23 provided a highly sensitive technique to detect CNN-LOH, as well as, copy-number changes in AML/MDS genomes.
|
|
|---|
|
View this table: [in a new window] [Download PPT slide] |
Table 1. Baseline clinical characteristics of 38 cases of normal karyotype acute myeloid leukemia/myelodysplastic syndrome.
|
|
View this table: [in a new window] [Download PPT slide] |
Table 4. Chromosomal location of small copy-number changes.
|
Determination of SNP sequences, JAK2, FLT3, NPM1, and AML1/RUNX1 mutations, and other target genes in cases of CNN-LOH
To determine the SNP sequences, (SNP identities are rs7747259, rs1122637, rs9505293, rs6934027, rs280153 and rs191986) in case #38 chromosome 6p region, the genomic region of each SNP site was amplified by genomic polymerase chain reaction (PCR) using specific primers. For determination of JAK2 V617F mutation in case #20, genomic PCR was performed with specific primers. PCR products were purified and sequenced. The sequences of the primers are shown in Online Supplementary Tables S1 and S2. To determine the FLT3-ITD mutation, the PCR reaction was performed with specific primers, and the PCR products were separated on a 2.0% agarose gel stained with ethidium bromide as described previously.34,35 Mutations at exon 12 of the NPM1 gene were determined using a melting curve-based LightCycler assay (Roche Diagnostics, Mannheim, Germany).36 Denaturing high-performance liquid chromatography analysis was performed to determine the AML1/RUNX1 mutation in case #17 as described previously.37 Alterations of several tyrosine kinase genes including FGR (case #3 and #23), DDR1 (case #2 and 38), TYK2 (case #2), MATK (case #2), FER (case #8) and FGFR4 (case #8) were determined by either nucleotide sequencing of their exons and/or band-shifts of PCR products of exons after their electrophoresis and visualization on a gel (single strand conformation polymorphism), as described previously38 with minor modifications. The PCR reaction contained genomic DNA, 500 nM of each of the primers, 200 nM of each of the dNTP, 0.5 units of Taq DNA polymerase and 3 µCi [
-32P] dCTP in 20 µL PCR products were diluted 10-fold in the loading buffer (10 mM NaOH, 95% formamide, and 0.05% of both bromophenol blue and xylene cyanol). After denaturation at 94°C for 5 min, 2 mL of the samples were loaded onto a 6% non-denaturating polyacrylamide mutation detection enhancement gel (Bioproducts, Rockland, ME, USA) with 10% (v/v) glycerol and separated at 300 V for 20 h. The gel was dried and subjected to autoradiography.
Quantitative real-time polymerase chain reaction
Gene-dosages of chromosome 6p24.3 in case #38, and the MYC and CDKN2A genes in case #20 were determined by quantitative real-time PCR (iCycler, Bio-Rad, Hercules, CA, USA) using Sybr Green. To determine the relative gene dosage of each sample, the chromosome 2p21 region was measured as a control.27 The copy-number of the 2p21 region was normal, as determined by SNP-chip analysis, in these samples. The delta threshold cycle value (
Ct) was calculated from the given Ct value by the formula
Ct = (Ct sample – Ct control). The fold change was calculated as 2–
Ct. Primer sequences are shown in Online Supplementary Table S2.
Gene expression microarray analysis
Total RNA was isolated from AML/MDS cells and processed according to Affymetrix guidelines for analysis with HGU133 Plus 2.0 microarrays. Data were analyzed with R version 2.5.0 using Bioconductor version 2.0.49 Data were normalized using the robust multi-array average procedure.39 Since most regions that showed chromosomal abnormalities were not recurring, we were not able to compare individual genes across samples with statistical tests. To assess plausibility of large deletions and amplifications, we subtracted from each gene (in the respective region) mean expression of this gene in other cases: case #11 was compared with 37 normal karyotype AML/MDS cases; and cases #20, #4 and #5 were compared with other normal karyotype AML/MDS samples. We then calculated a mean expression difference for each region and considered a value below zero to be consistent with deletion and a value above zero to be consistent with amplification.
|
|
|---|
|
View this table: [in a new window] [Download PPT slide] |
Table 2. Copy-number changes in case #11 detected by SNP-chip analysis.
|
![]() View larger version (41K): [in a new window] [Download PPT slide] |
Figure 1. Genomic DNA of 37 acute myeloid leukemia samples with normal karyotype were subjected to SNP-chip analysis; genomic abnormalities are summarized. Pink, green and red bars/boxes indicate CNN-LOH, deletion and duplication, respectively. Nineteen patients (51%) showed no detectable genomic abnormalities (data not shown), whereas 18 patients (49%) had one or more genomic abnormalities. Deletion or duplication was found in nine patients (24%), and CNN-LOH occurred in 12 patients (32%). Chromosomal location, size and genes are shown in Tables 3 and 4.
|
|
View this table: [in a new window] [Download PPT slide] |
Table 3. Chromosomal regions identified as CNN-LOH.
|
Cases #2, #8, #12, #17 and #20 had CNN-LOH on 19p (13.41 Mb), 5q (103.77 Mb), 12q (96.23 Mb), 21q (29.54 Mb) and 9p (43.96 Mb), respectively. Although these regions of CNN-LOH occurred in only one case each, several interesting genes were found in the region, including INSR, TYK2, and MATK (case #2); APC, FER, FMS/FLT4, PDGFRB, ITK and FGFR4 (case #8), AML1/RUNX1 (case #17), and JAK2 and TEK (case #20) (Table 3).
Interestingly, cases #10 and #21 had a FLT3-ITD gene mutation (Table 3); case #17 had an AML1/RUNX1 frameshift caused by a deletion of cytosine at nucleotide 211 (Table 3). Sequencing of JAK2 in case #20 showed a homozygous canonical JAK2 mutation [V617F (GTC
TTC)] (Table 3). Each of these mutations occurred at a CNN-LOH. The data suggest that removal of a normal allele and duplication of the mutated allele is favored by the cancer cells.
Validation of copy number-neutral loss of heterozygosity
To validate CNN-LOH, we determined SNP sequences and gene-dosage in a CNN-LOH region using case #38 (Figure 2). If a chromosome has LOH, the nucleotide at the SNP site should not be heterozygous, but should be homozygous. We, therefore, examined six independent SNP sites in case #38 on the chromosome 6p region of CNN-LOH including rs7747259, rs1122637, rs9505293, rs6934027, rs280153 and rs191986. All six SNP sites showed only a single nucleotide; no SNP sites showed heterozygosity (Figure 2B). Each one of these sites is heterozygous in the general population at a frequency varying between 25% and 42% (Entrez SNP database, http://www.ncbi.nlm.nih.gov/sites/entrez?db=snp). These results strongly suggest that this region has LOH.
![]() View larger version (13K): [in a new window] [Download PPT slide] |
Figure 2. Validation of CNN-LOH (A) Region of CNN-LOH in chromosome 6 of case #38. Red dots represent SNP sites as probes and indicate total copy-number. The blue line represents an average of copy-number and shows gene dosage. Green bars represents heterozygous (hetero) SNP calls. Red and green lines show allele-specific copy-number (AsCN). Blue bars indicate LOH detected by heterozygous SNP calls. (B) Determination of SNP sequences in the 6p region. Six independent SNP sites were sequenced. All six SNP sites contained only a single nucleotide; no SNP site displayed heterozygosity. Results are consistent with CNN-LOH. (C) Determination of gene-dosage in the 6p region. Gene-dosage of 6p24.3 (CNN-LOH region) in case #38 is compared to that in normal genomic DNA using quantitative genomic real-time PCR. Levels of the gene-dosage were determined as a ratio between 6p24.3 and the reference genomic DNA, 2p21.
|
Chromosomal regions of copy-number change detected by SNP-chip analysis
Nine patients (24%) had small copy-number changes including deletions and/or duplications; the median size of the duplications and deletions was 0.3 Mb (range, 0.09–4.33 Mb) and 0.625 Mb (range, 0.11–5.87 Mb), respectively. As shown in Table 4, hemizygous deletions were found at 14q21.2 (0.3 Mb, case #2), 17q11.2 (2.7 Mb, case #4), 12p13.31 - p13.2 (2.91 Mb, case #5), 21q21.2 (0.44 Mb, case #7), 2q36.2 (0.41 Mb, case #9), 2p23.1 (0.56 Mb, case #13), 4q24 (1.08 Mb, case #20), 9p21.3 - p21.2 (5.87 Mb, case #20), 3p26.3 (0.69 Mb, case #26), and 8p23.2 (0.11 Mb, case #41). Cases #4, #5 and #20 had duplication at 1q43 (0.09 Mb), 18q21.2 (0.3 Mb), and 8q24.13 - q24.21 (4.33 Mb), respectively. These regions contain well-known oncogenes and tumor suppressor genes (Table 4). The tumor suppressor genes, NF1 and CDKN2A/CDKN2B, and the transcription factor, ETV6/TEL were deleted in cases #4, #20 and #5, respectively; and the oncogene MYC was duplicated in case #20.
Validation of copy-number changes
Next, we validated copy-number changes in cases #20 and #5 using different techniques. Case #20 had duplication at 8q24.13 - q24.21 (Figure 3A) and hemizygous deletion at 9p21.3 - p21.2 (Figure 3B); these regions contain the oncogene MYC and the tumor suppressor genes CDKN2A and CDKN2B, respectively. Relative gene-dosages of the MYC and CDKN2A genes were examined by quantitative genomic real-time PCR with the chromosome 2p21 region as a control. The level of the MYC gene was about 2-fold higher while the level of the CDKN2A gene was approximately 10-fold lower compared with normal genomic DNA (Figures 3C and D).
![]() View larger version (17K): [in a new window] [Download PPT slide] |
Figure 3. Validation of duplication and deletion: (A) Chromosome 8q24.13-q24.21 is duplicated. This region contains the oncogene MYC. (B) Chromosome 9p21.3 - p21.2 shows a deletion. The deleted region contains the tumor suppressor genes CDKN2A (p16/INK4A and p14/ARF) and CDKN2B ( p15/INK4B ). (C, D) Gene-dosages of the MYC gene (C) and the CDKN2A gene (D) region in case #20 are compared to normal genomic DNA by quantitative genomic real-time PCR. Levels of the gene-dosage are determined as a ratio between target gene and the reference genomic DNA, 2p21. (E) Case #5 had hemizygous deletion in chromosome 12p13.31-p13.2; this region contains the transcription factor ETV6/TEL gene. Physical localization and size are presented in Table 4. (F) FISH analysis of case #5 with probes for the ETV6/TEL region. Probes of fluorescein-labeled ETV6-downstream (normal region by SNP-chip analysis) and Texas-red-labeled ETV6-upstream (hemizygous deleted region by SNP-chip analysis) revealed one and two signals, respectively.
|
Relationship between genomic abnormalities and mutant genes within the region
In our normal karyotype AML/MDS samples, eight cases (21%) had FLT3-ITD and 14 cases (37%) had a NPM1 mutation (Table 1). We compared genomic abnormalities, and FLT3-ITD and NPM1 mutations (Online Supplementary Table S4). Both FLT3-ITD and NPM1 were mutated in two samples in group A (11%) and four cases in group B (22%). A single mutation of FLT3-ITD was found in one sample in group A (5%) and one case in group B (6%); a single mutation of NPM1 occurred in five samples in group A (26%) and three samples in group B (17%). These mutations were, therefore, dispersed between both groups A and B
Relationship between genomic abnormalities and gene expression
We compared genomic abnormalities and gene expression. mRNA microarray analysis was done on all samples.40 First, the level of mRNA expression in case #11 was compared with that in 37 normal karyotype AML samples. Affymetrix microarray analysis showed decreased average gene expression in the deleted regions and increased gene expression for regions with trisomy: the difference of average expression of genes located on deleted regions of chromosomes 5, 7, 17, as well as, trisomy 8, 21 and 22 were –0.21±0.01, –0.16 ±0.013, –0.27±0.018, +0.21±0.012, +0.22±0.022 and +0.15±0.013 (mean difference ± standard error), respectively (Figure 4A and data not shown).
![]() View larger version (40K): [in a new window] [Download PPT slide] |
Figure 4. mRNA microarray analyses of acute myeloid leukemia samples. (A) Relationship between genomic abnormality and gene expression in acute myeloid leukemia case #11. mRNA microarray analysis was performed on all samples, and expression levels of acute myeloid leukemia cells from case #11 were compared to those of 37 normal karyotype acute myeloid leukemia samples. Affymetrix microarray analysis showed decreased average gene expression in the deleted regions, and increased average gene expression for trisomy 8: the difference of mean expression of genes located in the deleted region of chromosomes 5 (upper, left), 7 (upper, right), 17 (lower, right) and trisomy 8 (lower, left) were –0.21±0.01, –0.16±0.013, –0.27±0.018, and +0.21±0.012 (mean difference ± standard error), respectively. (B) Expression levels in acute myeloid leukemia cells from cases #20, #4 and #5 were compared with those in 36 normal karyotype AML samples. The differences in mean expression of genes located on the deleted region of chromosome 9 in case #20 (right), chromosome 17 in case #4 (middle), and chromosome 12 in case #5 (left) were –0.15±0.07, –0.37±0.07, and –0.23±0.051 (mean difference ± standard error), respectively. Each spot (black and red, Y-axis) indicates one gene and reflects the difference between each case and the mean of the other cases. Red spots represent genes located on an aberrant chromosome. The X-axis shows the chromosomal location.
|
|
|
|---|
We also found small copy-number changes in some cases. Several features of case #20 are worthy of comment. The MYC gene was duplicated, and the CDKN2A (p16/INK4A and p14/ARF) and CDKN2B (p15/INK4B) genes were hemizygously deleted. Prominent expression of C-MYC protein is associated with stimulation of p14/ARF which inactivates MDM2, producing greater levels of p53 resulting in either apoptosis or slowing of cell growth which allows for DNA repair.41,42 However, when the p14/ARF gene is deleted, C-MYC has an unfettered ability to stimulate growth of the cells. Case #20 had this constellation of changes. Furthermore, this individual had a homozygous JAK2 mutation. JAK2 is mutated (codon 617, valine changed to phenylalanine) and constitutively active in nearly 100%, 50% and 30% of samples from patients with polycythemia vera, agnogenic myeloid metaplasia and essential thrombocythemia, respectively, as well as in 1–3% of AML cases.43–45 We do not know the prior history of this individual.
Some of the deleted genes are of particular interest; first, the tumor suppressor gene NF1 was deleted in case #4. Children with neurofibromatosis type-1 have inactivating mutations of the NF1 and an increased risk of developing juvenile myelomonocytic leukemia,46 and LOH at the NF1 gene locus occurs in this form of leukemia and other cancers. A recent study showed that three of 103 T-ALL (3%) samples and two of 71 AML samples with MLL rearrangements (3%) had deletion of the NF1 gene region; a mutation in the remaining NF1 allele was found in three samples, suggesting that NF1 inactivation might be involved in the development of leukemia. Second, concerning case #5 (deletion of ETV6/TEL), ETV6/TEL is a transcriptional repressor and is involved in various translocations associated with leukemia. About 30% of AML patients have loss of expression of the ETV6/TEL protein;47,48 mutations of ETV6/TEL were found in 2% of AML samples, and these mutants behaved in a dominant-negative fashion.48 Interestingly, previous array-comparative genome hybridization analysis of normal karyotype AML showed duplication of 8q24.13-q24.21 (including the MYC gene) and deletion of 12p12.3 (including the ETV6 gene);33 this constellation of alterations was also observed in our study.
Our microarray analysis showed that regions with copy-number loss or gain of chromosomal material were associated with either decreased or increased mRNA expression of genes in that same region, respectively, demonstrating the relationship between chromosomal status and gene expression. From an analysis perspective, we applied a descriptive approach and intended to assess plausibility of data. Some genes do indeed have higher expression values in deleted regions (Figure 4A, red points above zero) than in other cases, and some genes have lower values in trisomy (Figure 4A, red points below zero) than in other cases. However on average, expression in deleted regions is clearly lower than in non-deleted cases.
Because most regions are not recurring, we compared only one sample versus the rest (i.e. case #11 was compared with 37 normal karyotype AML/MDS cases; and cases #20, #4 and #5 were compared with other normal karyotype AML/MDS samples.) Various technical and biological sources of noise can confound the analysis. Overall, expression data appear to be consistent with chromosomal deletions and amplifications of the investigated regions. Further studies in larger cohorts of patients should enable prognostic stratification of patients in relation to their genomic changes and reveal new therapeutic targets.
# Current address: Department of Stem Cell Biology, Graduate School of Medical Science, Kanazawa University. 13-1 Takara-machi, Kanazawa, Ishikawa 920-8640, Japan ![]()
The online version of this article contains a supplementary appendix.
TA performed research, analyzed the data and wrote the paper; SO and MS performed SNP-chip analyses; GY and YN developed the CNAG; NK, AY, CWM and MD assisted in data analyses; SS, CH and TH provided AML samples, performed FISH analysis and aided in data analyses; HPK directed the overall study.
The authors declare no competing financial interests.
Funding: this work was supported by NIH grants as well as the Parker Hughes Fund. HPK is the holder of the Mark Goodson endowed Chair in Oncology Research and is a member of the Jonsson Cancer Center and the Molecular Biology Institute, UCLA. MD and TH are supported by the European Leukemia Network (funded by the 6th Framework Program of the European Community). The study is dedicated to the memory of David Golde, a mentor and friend.
Received for publication March 5, 2008. Revision received September 17, 2008. Accepted for publication October 6, 2008.
|
|
|---|
fusion mRNA generated by the t(15;17) translocation in acute promyelocytic leukemia encodes a functionally altered RAR. Cell 1991;66:675-84.[CrossRef][Web of Science][Medline]
), in acute myeloid leukemia. Nat Genet 2001;27:263-70.[CrossRef][Web of Science][Medline]This article has been cited by other articles:
![]() |
T. Akagi, L.-Y. Shih, S. Ogawa, J. Gerss, S. R. Moore, R. Schreck, N. Kawamata, D.-C. Liang, M. Sanada, Y. Nannya, et al. Single nucleotide polymorphism genomic arrays analysis of t(8;21) acute myeloid leukemia cells Haematologica, September 1, 2009; 94(9): 1301 - 1306. [Abstract] [Full Text] [PDF] |
||||
| ||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||