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).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.
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.
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 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. 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.
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).
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.
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).
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.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.
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
- Kelly LM, Gilliland DG. Genetics of myeloid leukemias. Annu Rev Genomics Hum Genet. 2002; 3:179-98. PubMedhttps://doi.org/10.1146/annurev.genom.3.032802.115046Google Scholar
- Estey E, Dohner H. Acute myeloid leukaemia. Lancet. 2006; 368:1894-907. PubMedhttps://doi.org/10.1016/S0140-6736(06)69780-8Google Scholar
- Erickson P, Gao J, Chang KS, Look T, Whisenant E, Raimondi S. Identification of breakpoints in t(8;21) acute myelogenous leukemia and isolation of a fusion transcript, AML1/ETO, with similarity to Drosophila segmentation gene, runt. Blood. 1992; 80:1825-31. PubMedGoogle Scholar
- Miyoshi H, Kozu T, Shimizu K, Enomoto K, Maseki N, Kaneko Y. The t(8;21) translocation in acute myeloid leukemia results in production of an AML1-MTG8 fusion transcript. EMBO J. 1993; 12:2715-21. PubMedGoogle Scholar
- Licht JD. AML1 and the AML1-ETO fusion protein in the pathogenesis of t(8;21) AML. Oncogene. 2001; 20:5660-79. PubMedhttps://doi.org/10.1038/sj.onc.1204593Google Scholar
- Peterson LF, Zhang DE. The 8;21 translocation in leukemogenesis. Oncogene. 2004; 23:4255-62. PubMedhttps://doi.org/10.1038/sj.onc.1207727Google Scholar
- Liu P, Tarlé SA, Hajra A, Claxton DF, Marlton P, Freedman M. Fusion between transcription factor CBF beta/PEBP2 β and a myosin heavy chain in acute myeloid leukemia. Science. 1993; 261:1041-4. PubMedhttps://doi.org/10.1126/science.8351518Google Scholar
- Shigesada K, van de Sluis B, Liu PP. Mechanism of leukemogenesis by the inv(16) chimeric gene CBFB/PEBP2B-MHY11. Oncogene. 2004; 23:4297-307. PubMedhttps://doi.org/10.1038/sj.onc.1207748Google Scholar
- de Thé H, Chomienne C, Lanotte M, Degos L, Dejean A. The t(15;17) translocation of acute promyelocytic leukaemia fuses the retinoic acid receptor alpha gene to a novel transcribed locus. Nature. 1990; 347:558-61. PubMedhttps://doi.org/10.1038/347558a0Google Scholar
- de Thé H, Lavau C, Marchio A, Chomienne C, Degos L, Dejean A. The PML-RAR-α fusion mRNA generated by the t(15;17) translocation in acute promyelocytic leukemia encodes a functionally altered RAR. Cell. 1991; 66:675-84. PubMedhttps://doi.org/10.1016/0092-8674(91)90113-DGoogle Scholar
- Byrd JC, Mrózek K, Dodge RK, Carroll AJ, Edwards CG, Arthur DC. Pretreatment cytogenetic abnormalities are predictive of induction success, cumulative incidence of relapse, and overall survival in adult patients with de novo acute myeloid leukemia: results from Cancer and Leukemia Group B (CALGB 8461). Blood. 2002; 100:4325-36. PubMedhttps://doi.org/10.1182/blood-2002-03-0772Google Scholar
- Grimwade D, Walker H, Oliver F, Wheatley K, Harrison C, Harrison G. The importance of diagnostic cytogenetics on outcome in AML: analysis of 1,612 patients entered into the MRC AML 10 trial. The Medical Research Council Adult and Children’s Leukaemia Working Parties. Blood. 1998; 92:2322-33. PubMedGoogle Scholar
- Mrozek K, Marcucci G, Paschka P, Whitman SP, Bloomfield CD. Clinical relevance of mutations and gene-expression changes in adult acute myeloid leukemia with normal cytogenetics: are we ready for a prognostically prioritized molecular classification?. Blood. 2007; 109:431-48. PubMedhttps://doi.org/10.1182/blood-2006-06-001149Google Scholar
- Caligiuri MA, Schichman SA, Strout MP, Mrózek K, Baer MR, Frankel SR. Molecular rearrangement of the ALL-1 gene in acute myeloid leukemia without cytogenetic evidence of 11q23 chromosomal translocations. Cancer Res. 1994; 54:370-3. PubMedGoogle Scholar
- Gilliland DG, Griffin JD. The roles of FLT3 in hematopoiesis and leukemia. Blood. 2002; 100:1532-42. PubMedhttps://doi.org/10.1182/blood-2002-02-0492Google Scholar
- Reuter CW, Morgan MA, Bergmann L. Targeting the Ras signaling pathway: a rational, mechanism-based treatment for hematologic malignancies?. Blood. 2000; 96:1655-69. PubMedGoogle Scholar
- Pabst T, Mueller BU, Zhang P, Radomska HS, Narravula S, Schnittger S. Dominant-negative mutations of CEBPA, encoding CCAAT/enhancer binding protein-alpha (C/EBPα), in acute myeloid leukemia. Nat Genet. 2001; 27:263-70. PubMedhttps://doi.org/10.1038/85820Google Scholar
- Gombart AF, Hofmann WK, Kawano S, Takeuchi S, Krug U, Kwok SH. Mutations in the gene encoding the transcription factor CCAAT/enhancer binding protein alpha in myelodysplastic syndromes and acute myeloid leukemias. Blood. 2002; 99:1332-40. PubMedhttps://doi.org/10.1182/blood.V99.4.1332Google Scholar
- Falini B, Mecucci C, Tiacci E, Alcalay M, Rosati R, Pasqualucci L. Cytoplasmic nucleophosmin in acute myelogenous leukemia with a normal karyotype. N Engl J Med. 2005; 352:254-66. PubMedhttps://doi.org/10.1056/NEJMoa041974Google Scholar
- Falini B, Nicoletti I, Martelli MF, Mecucci C. Acute myeloid leukemia carrying cytoplasmic/mutated nucleophosmin (NPMc+ AML): biologic and clinical features. Blood. 2007; 109:874-85. PubMedhttps://doi.org/10.1182/blood-2006-07-012252Google Scholar
- Nannya Y, Sanada M, Nakazaki K, Hosoya N, Wang L, Hangaishi A. A robust algorithm for copy number detection using high-density oligonucleotide single nucleotide polymorphism genotyping arrays. Cancer Res. 2005; 65:6071-9. PubMedhttps://doi.org/10.1158/0008-5472.CAN-05-0465Google Scholar
- Engle LJ, Simpson CL, Landers JE. Using high-throughput SNP technologies to study cancer. Oncogene. 2006; 25:1594-601. PubMedhttps://doi.org/10.1038/sj.onc.1209368Google Scholar
- Yamamoto G, Nannya Y, Kato M, Sanada M, Levine RL, Kawamata N. Highly sensitive method for genomewide detection of allelic composition in nonpaired, primary tumor specimens by use of affymetrix single-nucleotide-polymorphism genotyping microarrays. Am J Hum Genet. 2007; 81:114-26. PubMedhttps://doi.org/10.1086/518809Google Scholar
- Pfeifer D, Pantic M, Skatulla I, Rawluk J, Kreutz C, Martens UM. Genome-wide analysis of DNA copy number changes and LOH in CLL using high-density SNP arrays. Blood. 2007; 109:1202-10. PubMedhttps://doi.org/10.1182/blood-2006-07-034256Google Scholar
- Lehmann S, Ogawa S, Raynaud SD, Sanada M, Nannya Y, Ticchioni M. Molecular allelokaryotyping of early stage untreated chronic lymphocytic leukemia. Cancer. 2008; 112:1296-305. PubMedhttps://doi.org/10.1002/cncr.23270Google Scholar
- Mullighan CG, Goorha S, Radtke I, Miller CB, Coustan-Smith E, Dalton JD. Genome-wide analysis of genetic alterations in acute lymphoblastic leukaemia. Nature. 2007; 446:758-64. PubMedhttps://doi.org/10.1038/nature05690Google Scholar
- Kawamata N, Ogawa S, Zimmermann M, Kato M, Sanada M, Hemminki K. Molecular allelokaryotyping of pediatric acute lymphoblastic leukemias by high resolution single nucleotide polymorphism oligonucleotide genomic microarray. Blood. 2008; 111:776-84. PubMedhttps://doi.org/10.1182/blood-2007-05-088310Google Scholar
- Flotho C, Steinemann D, Mullighan CG, Neale G, Mayer K, Kratz CP. Genome-wide single-nucleotide polymorphism analysis in juvenile myelomonocytic leukemia identifies uniparental disomy surrounding the NF1 locus in cases associated with neurofibromatosis but not in cases with mutant RAS or PTPN11. Oncogene. 2007; 26:5816-21. PubMedhttps://doi.org/10.1038/sj.onc.1210361Google Scholar
- Fitzgibbon J, Iqbal S, Davies A, O’shea D, Carlotti E, Chaplin T. Genome-wide detection of recurring sites of uniparental disomy in follicular and transformed follicular lymphoma. Leukemia. 2007; 21:1514-20. PubMedhttps://doi.org/10.1038/sj.leu.2404696Google Scholar
- Walker BA, Leone PE, Jenner MW, Li C, Gonzalez D, Johnson DC. Integration of global SNP-based mapping and expression arrays reveals key regions, mechanisms, and genes important in the pathogenesis of multiple myeloma. Blood. 2006; 108:1733-43. PubMedhttps://doi.org/10.1182/blood-2006-02-005496Google Scholar
- Raghavan M, Lillington DM, Skoulakis S, Debernardi S, Chaplin T, Foot NJ. Genome-wide single nucleotide polymorphism analysis reveals frequent partial uniparental disomy due to somatic recombination in acute myeloid leukemias. Cancer Res. 2005; 65:375-8. PubMedGoogle Scholar
- Fitzgibbon J, Smith LL, Raghavan M, Smith ML, Debernardi S, Skoulakis S. Association between acquired uniparental disomy and homozygous gene mutation in acute myeloid leukemias. Cancer Res. 2005; 65:9152-4. PubMedhttps://doi.org/10.1158/0008-5472.CAN-05-2017Google Scholar
- Tyybäkinoja A, Elonen E, Piippo K, Porkka K, Knuutila S. Oligonucleotide array-CGH reveals cryptic gene copy number alterations in karyotypically normal acute myeloid leukemia. Leukemia. 2007; 21:571-4. PubMedhttps://doi.org/10.1038/sj.leu.2404543Google Scholar
- Kiyoi H, Naoe T, Yokota S, Nakao M, Minami S, Kuriyama K. Internal tandem duplication of FLT3 associated with leukocytosis in acute promyelocytic leukemia: Leukemia Study Group of the Ministry of Health and Welfare (Kohseisho). Leukemia. 1997; 11:1447-52. PubMedhttps://doi.org/10.1038/sj.leu.2400756Google Scholar
- Schnittger S, Schoch C, Dugas M, Kern W, Staib P, Wuchter C. Analysis of FLT3 length mutations in 1003 patients with acute myeloid leukemia: correlation to cytogenetics, FAB subtype, and prognosis in the AMLCG study and usefulness as a marker for the detection of minimal residual disease. Blood. 2002; 100:59-66. PubMedhttps://doi.org/10.1182/blood.V100.1.59Google Scholar
- Schnittger S, Schoch C, Kern W, Mecucci C, Tschulik C, Martelli MF. Nucleophosmin gene mutations are predictors of favorable prognosis in acute myelogenous leukemia with a normal karyotype. Blood. 2005; 106:3733-9. PubMedhttps://doi.org/10.1182/blood-2005-06-2248Google Scholar
- Dicker F, Haferlach C, Kern W, Haferlach T, Schnittger S. Trisomy 13 is strongly associated with AML1/RUNX1 mutations and increased FLT3 expression in acute myeloid leukemia. Blood. 2007; 110:1308-16. PubMedhttps://doi.org/10.1182/blood-2007-02-072595Google Scholar
- Yin D, Xie D, Hofmann WK, Miller CW, Black KL, Koeffler HP. Methylation, expression, and mutation analysis of the cell cycle control genes in human brain tumors. Oncogene. 2002; 21:8372-8. PubMedhttps://doi.org/10.1038/sj.onc.1206031Google Scholar
- Bolstad BM, Irizarry RA, Astrand M, Speed TP. A comparison of normalization methods for high density oligonucleotide array data based on variance and bias. Bioinformatics. 2003; 19:185-93. PubMedhttps://doi.org/10.1093/bioinformatics/19.2.185Google Scholar
- Haferlach T, Kohlmann A, Schnittger S, Dugas M, Hiddemann W, Kern W. A global approach to the diagnosis of leukemia using gene expression profiling. Blood. 2005; 106:1189-98. PubMedhttps://doi.org/10.1182/blood-2004-12-4938Google Scholar
- Hemann MT, Bric A, Teruya-Feldstein J, Herbst A, Nilsson JA, Cordon-Cardo C. Evasion of the p53 tumour surveillance network by tumour-derived MYC mutants. Nature. 2005; 436:807-11. PubMedhttps://doi.org/10.1038/nature03845Google Scholar
- Dang CV, O’Donnell KA, Juopperi T. The great MYC escape in tumorigenesis. Cancer Cell. 2005; 8:177-8. PubMedhttps://doi.org/10.1016/j.ccr.2005.08.005Google Scholar
- Baxter EJ, Scott LM, Campbell PJ, East C, Fourouclas N, Swanton S. Acquired mutation of the tyrosine kinase JAK2 in human myeloproliferative disorders. Lancet. 2005; 365:1054-61. PubMedhttps://doi.org/10.1016/S0140-6736(05)71142-9Google Scholar
- Lee JW, Kim YG, Soung YH, Han KJ, Kim SY, Rhim HS. The JAK2 V617F mutation in de novo acute myelogenous leukemias. Oncogene. 2006; 25:1434-6. PubMedhttps://doi.org/10.1038/sj.onc.1209163Google Scholar
- Illmer T, Schaich M, Ehninger G, Thiede C, Tyrosine kinase mutations of JAK2 are rare events in AML but influence prognosis of patients with CBF-leukemias. Haematologica. 2007; 92:137-8. PubMedhttps://doi.org/10.3324/haematol.10489Google Scholar
- Side LE, Emanuel PD, Taylor B, Franklin J, Thompson P, Castleberry RP. Mutations of the NF1 gene in children with juvenile myelomonocytic leukemia without clinical evidence of neurofibromatosis, type 1. Blood. 1998; 92:267-72. PubMedGoogle Scholar
- Hernández JM, González MB, García JL, Ferro MT, Gutiérrez NC, Marynen P. Two cases of myeloid disorders and a t(8;12) (q12;p13). Haematologica. 2000; 85:31-4. PubMedGoogle Scholar
- Barjesteh van Waalwijk van Doorn-Khosrovani S, Spensberger D, de Knegt Y, Tang M, Löwenberg B, Delwel R. Somatic heterozygous mutations in ETV6 (TEL) and frequent absence of ETV6 protein in acute myeloid leukemia. Oncogene. 2005; 24:4129-37. PubMedGoogle Scholar
- R: A language and environment for statistical computing. R Foundation for Statistical Computing: Vienna, Austria; 2007. Google Scholar
- Serrano E, Carnicer MJ, Orantes V, Estivill C, Lasa A, Brunet S. Uniparental disomy may be associated with microsatellite instability in acute myeloid leukemia (AML) with a normal karyotype. Leuk Lymphoma. 2008; 49:1178-83. PubMedhttps://doi.org/10.1080/10428190802035941Google Scholar
- Gupta M, Raghavan M, Gale RE, Chelala C, Allen C, Molloy G. Novel regions of acquired uniparental disomy discovered in acute myeloid leukemia. Genes Chromosomes Cancer. 2008; 47:729-39. PubMedhttps://doi.org/10.1002/gcc.20573Google Scholar
- Raghavan M, Smith LL, Lillington DM, Chaplin T, Kakkas I, Molloy G. Segmental uniparental disomy is a commonly acquired genetic event in relapsed acute myeloid leukemia. Blood. 2008; 112:814-21. PubMedhttps://doi.org/10.1182/blood-2008-01-132431Google Scholar
- Gorletta TA, Gasparini P, D’Elios MM, Trubia M, Pelicci PG, Di Fiore PP. Frequent loss of heterozygosity without loss of genetic material in acute myeloid leukemia with a normal karyotype. Genes Chromosomes Cancer. 2005; 44:334-7. PubMedhttps://doi.org/10.1002/gcc.20234Google Scholar
- Tyybäkinoja A, Elonen E, Vauhkonen H, Saarela J, Knuutila S. Single nucleotide polymorphism microarray analysis of karyotypically normal acute myeloid leukemia reveals frequent copy number neutral loss of heterozygosity. Haematologica. 2008; 93:631-2. PubMedhttps://doi.org/10.3324/haematol.12232Google Scholar