Frequent genomic abnormalities in acute myeloid leukemia/myelodysplastic syndrome with normal karyotype
Tadayuki Akagi, Seishi Ogawa, Martin Dugas, Norihiko Kawamata, Go Yamamoto, Yasuhito Nannya, Masashi Sanada, Carl W. Miller, Amanda Yung, Susanne Schnittger, Torsten Haferlach, Claudia Haferlach, H. Phillip Koeffler

Author Affiliations

  1. Tadayuki Akagi1,#,
  2. Seishi Ogawa2,3,4,
  3. Martin Dugas5,
  4. Norihiko Kawamata1,
  5. Go Yamamoto2,
  6. Yasuhito Nannya2,
  7. Masashi Sanada3,4,
  8. Carl W. Miller1,
  9. Amanda Yung1,
  10. Susanne Schnittger6,
  11. Torsten Haferlach6,
  12. Claudia Haferlach6 and
  13. H. Phillip Koeffler1
  1. 1Division of Hematology and Oncology, Cedars-Sinai Medical Center, UCLA School of Medicine, Los Angeles, CA, USA
  2. 2Department of Hematology and Oncology and
  3. 3Department of Cell Therapy and Transplantation Medicine and the 21st Century COE Program, Graduate School of Medicine, University of Tokyo, Tokyo, Japan
  4. 4Core Research for Evolutional Science and Technology, Japan Science and Technology Agency, Tokyo, Japan
  5. 5Department of Medical Informatics and Biomathematics, University of Munster, Munster, Germany
  6. 6MLL Munich Leukemia Laboratory, Munich, Germany
  1. 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
View Abstract

Abstract

Background Acute myeloid leukemia is a clonal hematopoietic malignant disease; about 45–50% of cases do not have detectable chromosomal abnormalities. Here, we identified hidden genomic alterations and novel disease-related regions in normal karyotype acute myeloid leukemia/myelodysplastic syndrome samples.

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.

Introduction

Acute myeloid leukemia (AML) is a clonal malignant hematopoietic disease characterized by a block in differentiation, resulting in accumulation of immature myeloid cells.1,2 Karyotypic analyses have revealed several frequent chromosomal translocations producing fusion genes associated with AML. The t(8;21)(q22;q22) translocation is one of these abnormal karyotypes, and this translocation produces AML1-ETO fusion products.3,4 The AML1-ETO blocks hematopoietic differentiation and enhances self-renewal of human and murine hematopoietic stem cells.5,6 The fusion product apparently binds to AML1 target genes and represses their transcription.5,6 The inv(16)(p13q22) or t(16;16)(p13;q22) produces the leukemogenic CBFB-MYH11 fusion gene which blocks differentiation of hematopoietic stem cells by inhibiting the function of AML1.7,8 Acute promyelocytic leukemia cells usually have t(15;17)(q22;q11-21) producing PML-RARA fusion products which also behave as a transcriptional repressor.9,10 Other frequent translocations include t(9;11), t(6;11), inv(3)/t(3;3) and t(6;9).11 Trisomy 8, 11, 13, 21 and 22, and deletion of chromosome 5/5q, 7/7q, 17/17p and 20/20q also occur moderately frequently.2,11,12 About 45–50% of AML patients have no detectable chromosomal abnormalities.13,14 In general, these individuals with a normal karyotype in their leukemic cells show an intermediate prognosis.13,14

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).2123 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,5054

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.

Design and Methods

Patients’ samples

Samples from 30 anonymized patients with normal karyotype AML and 8 anonymized patients with normal karyotype MDS (age, 33–88 years; median, 62 years) were examined. These samples were isolated from bone marrow at diagnosis. The patients’ age, gender, diagnosis, white blood cell count (WBC), karyotype and additional mutations of FLT3 and NPM1 are summarized in Table 1. This study was approved by Cedars-Sinai Medical Center (IRB number 4485).

High-density SNP-chip analysis

Genomic DNA was isolated from AML/MDS cells, and the DNA was subjected to GeneChip Human mapping 250 K array NspI microarray (SNP-chip, Affymetrix, Santa Clara, CA, USA) as described previously.21,23 Hybridization, washing and signal detection were performed on GeneChip Fluidics Station 400 and GeneChip scanner 3000 according to the manufacturer’s protocols (Affymetrix). Microarray data were analyzed for determination of both total and allelic-specific copy-number using the CNAG program as previously described21,23 with minor modifications; the status of copy-numbers as well as CNN-LOH at each SNP was inferred using the algorithms based on hidden Markov models.21,23 GNAGraph software was used for clustering of AML samples with regards to their copy-number changes, as well as CNN-LOH.27 Size, position and location of genes were identified with UCSC Genome Browser http://genome.ucsc.edu. Copy-number changes, including duplication and deletion, were identified by allele-specific CNAG software.23,27 These copy-number changes include copy-number variant and physiological deletion at the immunoglobulin and T-cell receptor genes. Copy-number variants as described previously at http://projects.tcag.ca/variation and physiological deletions were eliminated manually, and other regions detected by allele-specific CNAG software are listed on Table 4.

Fluorescence in situ hybridization analysis

Bone marrow samples from AML patients were used for interphase fluorescence in situ hybridization (FISH) analysis. The FISH studies were performed using the following probes: D5S721 (5p15.2), D5S23 (5p15.2), D7Z1 (centromere of chromosome 7), ABL (9q34.12), EGR1 (5q31.2), D7S486 (7q31), TP53 (17p13.1), D8Z2 (centromere of chromosome 8), AML1 (21q22.12) and BCR (22q11.23) (ABBOTT/VYSIS, Des Plaines, IL, USA). Probes for the 12p13 region [fluorescein-labeled ETV6-downstream region (483 kb-length) and Texas-red-labeled ETV6-upstream region (264 kb-length)] were used for FISH analysis in case #5. The ETV6 probes were obtained from ABBOTT/VYSIS.

Table 1.

Baseline clinical characteristics of 38 cases of normal karyotype acute myeloid leukemia/myelodysplastic syndrome.

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.

Results

Proof of principal

To identify hidden abnormalities in AML/MDS with a normal karyotype, 37 samples were analyzed by 250K SNP-chip microarray. One additional case (case #11) had only 13 metaphases and chromosomal abnormalities were not detected on karyotypic analysis; this sample did, however, have numerous genetic abnormalities identified by SNP-chip including hemizygous deletions of 3p25.1-p24.3 (2.29 Mb), 3p24.2-p24.1 (3.96 Mb), 3p23-q12.1 (66.55 Mb), 5q11.2-q-terminal (124.89 Mb), 7q11.23-q36.1 (76.04 Mb), 7q36.2 (0.78 Mb), 11q23.3-q-terminal (18.24 Mb), 17p-terminal-q11.1 (22.48 Mb), and 17q11.2-q12 (4.42 Mb); duplications of 3p24.3 (2.14 Mb), 5p15.31 (1.83 Mb), and 5p14.3-q11.2 (35.53 Mb); and trisomy of chromosomes 8, 21 and 22 (Table 2). To confirm these SNP-chip results, we performed extensive FISH analysis. The number of signals for probes D5S721 (5p15.2), D5S23 (5p15.2), D7Z1 (centromere of chromosome 7) and ABL (9q34.12) was two, and SNP-chip analysis also showed normal copy number (2n) consistent with the SNP-chip data. The EGR1 (5q31.2), D7S486 (7q31) and TP53 (17p13.1) probes revealed one signal; and these regions also showed hemizygous deletion (1n) by SNP-chip analysis. D8Z2 (centromere of chromosome 8), AML1 (21q22.12) and BCR (22q11.23) probes showed three or four signals, and SNP-chip analysis also indicated trisomy (3n) of these chromosomes. Chromosome 9 was normal by both SNP-chip and FISH analyses. As summarized in Online Supplementary Table S3, the results of SNP-chip and FISH analyses were completely congruent. Taken together, these results suggest that SNP-chip analysis reflected the genomic changes.

SNP-chip analysis in 37 normal karyotype acute myeloid leukemia/myelodysplastic syndrome samples

SNP-chip analysis of samples from 37 patients with normal karyotype AML/MDS revealed several genomic copy-number changes, as well as CNN-LOH. Nineteen patients (51%) had a normal genome by SNP-chip analysis (group A). In contrast, 18 patients (49%) had one or more genomic abnormalities (group B) (Figure 1). Deletions and/or duplications were found in nine patients (24%). Twelve patients (32%) had CNN-LOH. In group B, 14 cases (78% of the 18 samples) had only one genomic change; one case (6%) had two genomic abnormalities (case #5); two cases (11%) had three changes (case #2 and #4) and one case (5%) had four genomic alterations (case #20).

We also compared the relationship between the genomic changes and the French-American British classification of the 15 AML and 3 MDS samples in group B. In the AML samples, 11 cases had CNN-LOH, three cases had a duplication and seven cases had a deletion. The one AML M1 sample (case #10) had CNN-LOH; and the two AML M5b samples (cases #3 and #37) had CNN-LOH in their chromosomes. In the four AML M4 samples, cases #38, #21 and #2 had CNN-LOH, and cases #2 and #7 had a small deletion. In eight AML M2 samples, five (cases #4, #8, #17, #20, and #23) had CNN-LOH, three (cases #4, #5 and #20) had a duplication, and five (cases # 4, #5, #9, #20 and #26) had a deletion. In the three MDS samples, one sample (case #12) had CNN-LOH, and two samples (cases #13 and #41) had a deletion (Figure 1, Tables 3 and 4). Taken together, these results show that the patients who were categorized as having normal karyotype AML/MDS had easily recognizable deletions, duplications and/or CNN-LOH of their genome.

Table 2.

Copy-number changes in case #11 detected by SNP-chip analysis.

Chromosomal region and candidate genes in CNN-LOH detected by SNP-chip analysis

Previous studies demonstrated CNN-LOH in AML samples at a frequency of 15–20%.31,32,50,51,53,54 Our analysis with AsCNAR (allele-specific copy-number analysis using anonymous references) revealed CNN-LOH in 32% of the AML/MDS samples with a normal karyotype; the median size of the CNN-LOH was 30.91 Mb (range, 11.76 Mb–103.77 Mb). We found some cases with a recurrent region of CNN-LOH. Cases #3 and #23 had CNN-LOH on 1p, and the common region of CNN-LOH (30.85 Mb) included the tyrosine kinase genes (FGR, EPHA2 and EPHB2) and an imprinted tumor suppressor gene TP73 (Table 3). Cases #2 and #38 had CNN-LOH on 6p and the common region of CNN-LOH (30.97 Mb) contained the tyrosine kinase gene DDR1 (Table 3). Cases #4 and #37 had CNN-LOH on 8q and the common region of CNN-LOH (11.76 Mb) contained the tyrosine kinase gene PTK (Table 3). CNN-LOH of the whole region of 13q was found in cases #10 and #21; this region contains the FLT3, FLT1, BRCA2 and RB1 genes (Table 3).

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.

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.

Next, we determined gene-dosage of the region to exclude the possibility of hemizygous deletion. The gene-dosage of 6p24.3 in case #38 was compared to that of normal genomic DNA using quantitative genomic real-time PCR by comparing the ratio between 6p24.3 and the reference genomic DNA, 2p21. As shown in Figure 2C, the amount of DNA at this site for case #38 was almost the same as that for normal genomic DNA, indicating that this region is not deleted. Taken together, our sequence data and gene dosage study validated the results of our SNP-chip analysis, clearly showing CNN-LOH at 6p24.3.

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

Chromosome 12p13.31 - p13.2 was deleted in case #5; this region contains the transcription factor ETV6/TEL (Figure 3E). FISH analysis with a probe of fluorescein-labeled ETV6-downstream (normal copy-number region) revealed two signals and a probe of Texas-red-labeled ETV6-upstream (hemizygous deleted region) revealed one signal (Figure 3F), validating the observations from SNP-chip analysis.

Table 3.

Chromosomal regions identified as CNN-LOH.

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.

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

Table 4.

Chromosomal location of small copy-number changes.

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

Next, we examined the relationship between small copy-number changes and mRNA expression levels in the region. For this analysis, we chose deleted regions on chromosome 9 in case #20 (Figure 3B), chromosome 17 in case #4 (Table 4) and chromosome 12 in case #5 (Figure 3E). The differences in mean expression of genes located in deleted regions of chromosomes 9 (case #20), 17 (case #4), and 12 (case #5) were −0.15±0.07, −0.37±0.07, and −0.23±0.051 (mean difference ± standard error), respectively (Figure 4B). These results showed that large and small copy-number changes led to alterations of mRNA expression. In addition, the difference in mean expression of genes located in the CNN-LOH regions of each sample was comparable to that in normal copy-number samples, suggesting that CNN-LOH does not contribute to aberrant levels of gene expression (data not shown).

Discussion

Our genome-wide SNP-chip analysis of normal karyotype AML/MDS showed that 49% of these samples had one or more genomic abnormalities including deletions, duplications and CNN-LOH. Previous studies demonstrated that CNN-LOH occurs in AML samples at a frequency of 15–20%.31,32,50,51,53,54 Of interest, about 40% of cases of relapsed of AML had CNN-LOH.52 In our analysis, 32% of samples had CNN-LOH, and these regions of CNN-LOH contain several tyrosine kinase and tumor suppressor genes that may be candidate target genes in normal karyotype AML/MDS. In fact, the FLT3-ITD (13q12.2), JAK2 V617F (9p24.1) and deletion of a cytosine at nucleotide 211 of AML1/RUNX1 (21q22.12) occurred in areas of CNN-LOH resulting in duplication of these mutant genes and loss of the normal allele. A prior paradigm was that CNN-LOH marked the location of a mutated tumor suppressor gene, but it is clear that CNN-LOH can also be a signpost of an activated (mutated) oncogene. Of note, several CNN-LOH, including a region on chromosome 1p (cases #3 and #23), 6p (cases #2 and #38), 8q (cases #4 and 37) and 13q (cases #10 and #21), occurred in more than one sample. In addition, CNN-LOH of these regions, as well as several other unique CNN-LOH regions in our cohort, were also found in other studies.50,51,53 Although these alterations are not frequent, shared regions of CNN-LOH clearly highlight their importance. These findings prompted us to screen genes located in CNN-LOH regions. We focused on tyrosine kinase genes including FGR (cases #3 and #23), DDR1 (cases #2 and #38), TYK2 (case #2), MATK (case #2), FER (case #8) and FGFR4 (case #8), and either determined their exon nucleotide sequences or looked for single strand conformation polymorphism band-shifts of PCR products of the exons. However, these genes did not have detectable mutations (data not shown). Nevertheless, we believe that these CNN-LOH, as well as deletions and duplications, are acquired somatic mutations. We examined these regions for known copy-number polymorphisms (web site, http://projects.tcag.ca/variation) and found none. Also previously, we compared SNP-chip data between matched samples of acute promyelocytic leukemia and normal genomic DNA from the same individual (Akagi et al., unpublished data) and found that CNN-LOH occurred only in the leukemia samples but not in the corresponding germline DNA. Furthermore, SNP-chip analysis easily detected a deletion on chromosome 3 (0.69 Mb) in case #28 in the AML sample which was not present in the remission bone marrow sample from the same individual (Online Supplementary Figure S1). Taken together, these findings suggest that the alterations detected by SNP-chip analysis are somatic mutations.

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.4345 We do not know the prior history of this individual.

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.

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.

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.

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.

Acknowledgments

We thank members of our laboratory for helpful discussions.

Footnotes

  • CH and HPK are co-last authors.

  • # 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.

  • Authorship and Disclosures

    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 March 5, 2008.
  • Revision received September 17, 2008.
  • Accepted October 6, 2008.

References

  1. 1.
  2. 2.
  3. 3.
  4. 4.
  5. 5.
  6. 6.
  7. 7.
  8. 8.
  9. 9.
  10. 10.
  11. 11.
  12. 12.
  13. 13.
  14. 14.
  15. 15.
  16. 16.
  17. 17.
  18. 18.
  19. 19.
  20. 20.
  21. 21.
  22. 22.
  23. 23.
  24. 24.
  25. 25.
  26. 26.
  27. 27.
  28. 28.
  29. 29.
  30. 30.
  31. 31.
  32. 32.
  33. 33.
  34. 34.
  35. 35.
  36. 36.
  37. 37.
  38. 38.
  39. 39.
  40. 40.
  41. 41.
  42. 42.
  43. 43.
  44. 44.
  45. 45.
  46. 46.
  47. 47.
  48. 48.
  49. 49.
  50. 50.
  51. 51.
  52. 52.
  53. 53.
  54. 54.
View Abstract