Abstract
Deletions on chromosome 13 are thought to be one of the most important prognostic features in Multiple Myeloma (MM). The biology underlying this is, however, uncertain. Chromosome 13 abnormalities have been evaluated conventionally by FISH using probes for 13q14, covering the retinoblastoma gene (RB1) region. Typically, for recurrent regions of loss of heterozygosity (LOH) it is possible to map a minimally deleted region within which an important gene may be located. This should be the case with 13q−, or alternatively there may be linkage with another genetic lesion, which could be contributing to the poor prognosis. Following the implementation of high-density single-nucleotide polymorphism (SNP) array, it is now possible to genotype the whole human genome with a mapping resolution of less than 50 Kb. Thus, the SNP array approach offers an opportunity to analyze both copy number abnormalities and LOH simultaneously. The aim of this study was to determine the numerical alterations, LOH and changes in the gene expression profile of the chromosome 13 in MM, and its possible association with other genetic events. For this purpose, we analyzed 17 patients included on the Myeloma IX trial with deletion of 13q14 compared with 22 samples without deletion, using Affymetrix 50K SNP arrays and Affymetrix U133 Plus 2 expression array. IGH translocations and 13q deletion were determined by FISH. dChipSNP and WGSA programs were used to analyze the data. With respect to 13q14, there was 100% correlation between FISH and SNP array results. 16 out of 17 cases with deletion of the RB1 gene by FISH analysis showed loss of 13q arm by SNP array, demonstrating that loss of the whole chromosome 13 is responsible for 13q deletions found in MM in >90% of cases, with only one case showing a defined region of deletion of chromosome 13 (13q14.11–13q21.2). Using gene expression arrays we could not define a specific pattern characteristic of expression loss in genes at 13q. Lower RB1 expression levels were not only restricted to cases with del(13). However, samples containing IGH translocations (t(11;14) and t(4;14)) without del(13) showed up to 4 times more RB1 expression, suggesting that MM evolution in cases containing IGH translocations is independent of RB1 expression. Interestingly, the hyperdiploid cases with and without del(13) expressed similar level of RB1. We also investigated whether other key cell cycle regulatory genes were associated with del(13); in particular, 4 cases showed 9p21 LOH by SNP array and no different gene expression levels, which suggest that LOH does not seem to be a mechanism of lost of expression of CDKN2A, CDKN2B and p14/ARF. We could not find any significant correlation with del(13) and expression of cell cycle regulatory genes, apart from 8/17 samples with del(13) that had low expression of p53 gene, including 6 t(4;14) cases and 2 t(11;14) cases. Also, 2 cases without monosomy 13 (1 with t(4;14) and 1 with t(11;14)), showed low p53 expression levels. However, SNP array data did not show any deletion at 17p in 38 cases, with the exception of a case with monosomy 13 and t(11;14) in which SNP array data showed loss at 17pter-17q21.2 and FISH detected p53 deletion. Further investigation between the association of p53 and del(13) are ongoing and maybe useful in defining the biology of this poor subgroup of patients.
Integrative DNA copy number detection and genotyping from sequencing and array-based platforms
