|
|
|||||||
Multiple Myeloma |
From the Molecular Cytogenetics Group, Centro Nacional de Investigaciones Oncológicas (CNIO), Madrid, Spain (CL, SA, JS, BF, DB, JCC); Department of Genetics, University of Navarra, Pamplona, Navarra, Spain (BS, MJC); Area of Cancer, Area of Cell Therapy and Haematology Service, University Clinic, University of Navarra and Foundation for Applied Medical Research, Pamplona, Navarra, Spain (FP)
Correspondence: Juan C. Cigudosa, PhD, Molecular Cytogenetics Group, Centro Nacional de Investigaciones Oncológicas (CNIO), Melchor Fernández Almagro 3, 28029 Madrid, Spain. E-mail address: jccigudosa{at}cnio.es
| ABSTRACT |
|---|
|
|
|---|
Design and Methods: For the first time, a high resolution array CGH with more than 40,000 probes, has been used to analyze 26 primary MM samples after the enrichment of CD138-positive plasma cells.
Results: This approach identified copy number imbalances in all cases. Bioinformatics strategies were optimized to perform data analysis allowing the segregation of hyperdiploid and non-hyperdiploid cases by array CGH. Additional analysis showed that structural chromosome rearrangements were more frequently seen in hyperdiploid cases. We also identified the same Xq21 duplication in nearly 20% of the cases, which originated through unbalanced chromosome translocations. High level amplifications and homozygous deletions were recurrently observed in our series and involved genes with meaningful function in cancer biology.
Interpretation and Conclusions: High resolution array CGH allowed us to identify copy number changes in 100% of the primary MM samples. We segregated different MM subgroups based on their genomic profiles which made it possible to identify homozygous deletions and amplifications of great genetic relevance in MM.
Key words: multiple myeloma, arrayCGH, hyperdiploid, amplification, homozygous deletion.
Multiple myeloma (MM) is an incurable malignant plasma cell (PC) neoplasia characterized by the accumulation of malignant plasma cells in the bone marrow. There is a wide variability in clinical features, responses to treatment, and survival times among patients.1 Chromosome aberrations, present in virtually all patients with MM2 include chromosome translocations involving the IG loci and/or copy number changes of partial or whole chromosomes.3 Two major genetic subtypes of MM have been defined: the hyperdiploid variant (H-MM), associated with multiple chromosome trisomies, and the non-hyperdiploid variant (NH-MM) with a high prevalence of IGH translocations.4 This classification seems to have relevance in prognosis since, recently, H-MM has been associated with a better survival and a similar response rate to treatment when compared to NH-MM.5
Chromosome-based comparative genomic hybridization (CGH) has shown that DNA copy number variations (CNV) in MM may affect almost all chromosomes.6–8 We recently found, by array-based CGH (array CGH) CNV that resulted in the overexpression of genes included in amplified chromosomal regions.9 Here, we characterized CNV in a panel of CD138-enriched primary MM samples with high density array CGH.
| Design and Methods |
|---|
|
|
|---|
|
Array CGH assays
Human Genome CGH 44k microarrays (Agilent Technologies, Palo Alto, CA, USA) were used in this study. This platform consists of approximately 44000 60-mer oligonucleotide probes that span the human genome with an average resolution below 50 Kb. The platform has gene focused coverage in order to ensure adequate coverage in most commonly studied genomic regions. Array CGH assays were performed as previously described11 with minor modifications. The DNA pool from ten helathy females was used as the reference DNA for all hybridizations. Slides were scanned using an Agilent 2565AA DNA Microarray Scanner (Agilent Technologies). Raw data can be freely accessed via our ftp server at www.cnio.es.
Images and data analysis
Microarray images were transformed to fluorescence intensities using Feature Extraction Software, version 8.1 (Agilent Technologies). Data were analyzed using Agilent CGH Analytics 3.2.25 software. Quality criteria were a derivative log ratio (DLR) spread lower than 0.3 log units and a signal to noise ratio for each channel greater than 30. DLR spread metrics estimates the log ratio noise by calculating the spread of log ratio differences between consecutive probes along all chromosomes. In order to establish gained and lost regions, a 500 Kb weighted moving-average window and a
-score threshold of ±2.5 were used.
score values were obtained for each probe and were categorized into gains and losses. Using this threshold, high level amplifications were defined as consecutive clones, spanning a region, which showed a log10 ratio higher than 0.3 (more than five copies) within a gained region and homozygous deletions as consecutive clones with a log10 ratio lower than –1 (no copies) and detected as losses by
–score. Those aberrant regions that did not reach these cut off values were considered low copy gains or single losses. For the subsequent analyses, regions of gains and losses were established. A region is defined as a group of at least three consecutive probes that are simultaneously and similarly changed.
Aberrations were categorized as numerical and structural changes. The latter type was divided into changes greater than 3Mb (>3Mb) and those less than 3Mb (<3Mb) regarding high resolution conventional CGH resolution (detailed data are provided in Supplementary Table 3).
Identification of H-MM and NH-MM
H-MM and NH-MM were easily classified in our series by array CGH. Based on the array CGH profile a copy number score (CNS) value was assigned to each chromosome. When there was not a whole chromosome gain or loss, a CNS value of 2 was assigned, when an entire chromosome gain occurred, a CNS of 3 was assigned and, finally, when a loss occurred the CNS was equal to 1. The sum of CNS values of each case was calculated and if the value was from 48 to 74 the case was classified as H-MM; when a different CNS was obtained, it was considered NH-MM. Ploidy status was further confirmed by a standard FISH assay based on centromeric probes (Vysis, Downers Grove, IL, USA).
Definition of smallest overlapping regions of imbalance (SORI)
Genomic aberration region data were used to obtain the SORI. Each SORI represented the most recurrent region of imbalance whose limits were defined by the minimum number of clones with the highest percentage of recurrence.12 SORI were defined independently for gains and losses. They were evaluated to define regions targeted by overlapping events in two or more samples. SORI limits were defined by the clones at the boundaries of the gained or lost region of the case showing the smallest unbalanced region. After SORI delimitation the presence or absence of belonging aberration was assigned to each sample (0=absence, 1=presence).
Bioinformatics and statistical analysis
GEPAS Suite software (http://gepas.bioinfo.cipf.es/) was used for unsupervised clustering, supervised analysis and gene annotation. A region was considered as statistically significant if the adjusted p value was below 0.05. For statistical comparison of the difference in the number of copy number changes in different groups, the U Mann-Whitney test was used with the SPSS v 13.0 package.
FISH and spectral karyotyping (SKY) analysis
FISH assays were carried out as described elsewhere.13 In order to identify whether there were any translocations involving the immunoglobulin (Ig) heavy and light chain loci, FISH assays were conducted with LSI IGH Dual Color Break Apart probe and IGL break apart probe.10 A noncommercial probe was used for the analysis of Xq21-qter duplication (Supplementary Table 3). Total Chromosome DNA X Probe (Qbiogene, Carlsbad, CA, USA) was used in combination with LSI IGH Dual Color Break Apart probe (Vysis), in order to study the involvement of IGH in Xq21 translocation. The probes used for the validation of amplification and homozygous deletions are described in supplementary Table 3.
|
| Results |
|---|
|
|
|---|
Clustering analysis identified two different H-MM subgroups
In order to obtain a reliable picture of common aberrations present in this series, we defined the SORI for gains and losses. Sixty-eight SORI were identified: 33 gains and 35 losses (Table 2). The most frequent gained SORI was 19p13.12-pter (65%) followed by 9q34, 19q13 and 19q13.2 (50%) and the most frequent losses were in 13q12-q21 (38%), 13q12.12 (35%), 8p12-p22 (27%) and 16p13.12 (23%). The median size was 1.5 Mb ranging from 0.026 Mb (a few genes) to 102 Mb (whole chromosomes). Deletion at 16p13.12, present in 23% of the cases, affected only ERCC4, a modifier gene reported to be involved in breast cancer.15 Unsupervised and supervised clustering analyses were performed using SORI as variables (Figure 1A). Unsupervised analysis was able to segregate H-MM from NH-MM, and after a supervised analysis, gained regions in chromosomes 15, 5, 19, 3, 9, 7 and 11q were found to be significantly (adjusted p-value <0.05) more frequent in the H-MM group. In the unsupervised analysis, our H-MM cases appeared to be distributed in two subgroups that statistically differed (adjusted p-value=0.03, Fishers exact test) according to the presence of gains in chromosome 7. Similarly, the unsupervised clustering analysis segregated two different subgroups within the NH-MM cases. One of the clusters, formed of only two cases, was characterized by the presence of a high number of CNV, such as gains in 8q24.21 (C-MYC) and absence of IGH translocations.
|
|
A novel duplication on Xq21.33 was identified in five out of the 26 primary cases and in the MM cell line L363. FISH studies allowed the mapping of the duplication breakpoints that clustered around a 10 Mb region in Xq21.33 (from nearly 86122579bp to 98277428 bp). SKY analysis (Figure 1) demonstrated that duplications were the unbalanced outcome of translocations with breakpoints at Xq21: in L363 involving chromosome 14q, t(X;14)(q21.33; q31.1); and in case #20 involving chromosome 21; t(X;21)(q21.33;q22.3). The breakpoint in chromosome 14 was revealed out by the array CGH profile and we confirmed that IGH was not involved by using a combination of a chromosome X painting probe (Qbiogene) and the LSI IGH Dual Color Break Apart probe (Vysis) in the same hybridization (Supplementary Figure 2).
| Discussion |
|---|
|
|
|---|
In this study, an array CGH platform was used for the analysis of copy number imbalances in CD138+ plasma cells. Although, as expected, most of the samples showed an unaltered karyotype at diagnosis, we found CNV in all of them. Observed frequencies of CNV, losses of chromosome 13q and duplications of chromosome 1q, and Ig genes translocations, were concordant with those in previous studies.6,7,9,12,17 Array CGH enabled efficient segregation of the series into H-MM and NH-MM cases; half of the cases were H-MM and half NH-MM. The presence of IGH translocations was associated with NH-MM; among the H-MM cases, only one showed IGH translocations and two IGL translocations. It has been suggested that translocations present in H-MM are secondary events.5,17
Array CGH provides a level of resolution that overcomes the limits of chromosome-based CGH which is, in the best conditions, near 3 Mb. Due to the high resolution of the platform (44K) we identified several aberrations smaller than this limit, distributed in both H-MM and NH-MM cases, which had not been previously described. We delimited several SORI, offering a reliable picture of the aberrations occurring in MM.12 Apart from segregating the two major genetic categories of MM, the unsupervised clustering based on the SORI aberration (Figure 1A) identified two H-MM subgroups, one of which was associated with chromosome 7 gains. Gains in 7q have been described in Burkitt's lymphoma,18 a combination in which they were associated with a poor outcome.
The approach used also allowed the identification of novel genomic markers and the definition of high level amplifications and homozygous deletions. Nine high level amplifications were identified. Some of them, such as 9q22.2, 19q13 and 16q22.2, had not been previously described in MM, but may harbor candidate oncogenes, as we have demonstrated in mantle cell lymphoma, in which we observed primary samples showing amplification and overexpression of SYK,19 or CCNE1, in 19q13, which has been shown to be amplified and overexpressed in endometrial and bladder neoplasia.20,21 Most of the homozygous deletions we found had not been previously described in MM. Reinforcing the putative pathogenic role of this finding, these regions also showed up as single copy losses in several additional samples of our series. As regards homozygous deletion in chromosome 13, ZNF198 seems to be a good candidate for the target of the deletion because, besides the deletion described in case #20, in case #15, the homozygous deletion covered only the ZNF198 gene. Furthermore, candidate genes such as PORIMIN, YAP1 and MMP8, which are located in the homozygous deletion that we found at 11q22.2, were reported to be deleted and down regulated in a previous study. Five cases showed duplication on Xq21.33. Although some CGH studies have reported data on Xq aberrations,22,23 we provide here the first description of the frequency and genetic origin of this marker in MM. In summary, our study represents the analysis of CNV in MM at the highest level of resolution ever reported. Our work includes some data confirming results obtained with other technological approaches, and provides a complete genomic characterization of the two genetic categories of MM. Additionally, we have demonstrated for the first time that extreme genomic aberrations, such as high level amplifications and homozygous deletions, are frequent events in MM. We believe that these data warrant further investigation to gain a deeper insight into the role of genomic instability in MM and its potential use for therapeutic purposes.
| Acknowledgements |
|---|
|
|
|---|
| Footnotes |
|---|
CL carried out all technical arrays, SKY and some FISH experiments. She contributed to the writing of the paper; BS performed FISH assays and contributed to the writing the paper; SA designed the study and contributed to the microarray analisis of the data; JS, BF and DB carried out microarray analyses and contributed to the analysis of the data; FP has sustantially contributed substantially with biological material and clinical data; MJC contributed with essential biological material, cytogenetic information and helped to design the experiments; JCC designed the work, collected all data and wrote the paper.
The authors reported no potential conflicts of interest.
Funding: CL and BS have PhD fellowships from Gobierno de Navarra. JS has a PhD fellowship from Ministerio de Educación y Ciencia. BF has a Marie Curie PhD Early Stage Research Training Fellowship. Work funded by Red Temática FIS G03/136 Mieloma Múltiple y otras Ganmapatías: de la Génesis a la Terapéutica and by Grant PI040555 from Fondo Investigaciones Sanitarias, Ministerio de Sanidad y Consumo to JCC.
Received for publication November 16, 2006. Accepted for publication March 15, 2007.
| References |
|---|
|
|
|---|
This article has been cited by other articles:
![]() |
B. I. Ferreira, J. F. Garcia, J. Suela, M. Mollejo, F. I. Camacho, A. Carro, S. Montes, M. A. Piris, and J. C. Cigudosa Comparative genome profiling across subtypes of low-grade B-cell lymphoma identifies type-specific and common aberrations that target genes with a role in B-cell neoplasia Haematologica, May 1, 2008; 93(5): 670 - 679. [Abstract] [Full Text] [PDF] |
||||
![]() |
A. Varshavsky From the Cover: Targeting the absence: Homozygous DNA deletions as immutable signposts for cancer therapy PNAS, September 18, 2007; 104(38): 14935 - 14940. [Abstract] [Full Text] [PDF] |
||||
| ||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||
| HOME | HELP | FEEDBACK | TABLE OF CONTENTS | ARCHIVE | SUBSCRIPTIONS |