Gene expression profiling for molecular classification of multiple myeloma in newly diagnosed patients

Blood ◽  
2010 ◽  
Vol 116 (14) ◽  
pp. 2543-2553 ◽  
Author(s):  
Annemiek Broyl ◽  
Dirk Hose ◽  
Henk Lokhorst ◽  
Yvonne de Knegt ◽  
Justine Peeters ◽  
...  

Abstract To identify molecularly defined subgroups in multiple myeloma, gene expression profiling was performed on purified CD138+ plasma cells of 320 newly diagnosed myeloma patients included in the Dutch-Belgian/German HOVON-65/GMMG-HD4 trial. Hierarchical clustering identified 10 subgroups; 6 corresponded to clusters described in the University of Arkansas for Medical Science (UAMS) classification, CD-1 (n = 13, 4.1%), CD-2 (n = 34, 1.6%), MF (n = 32, 1.0%), MS (n = 33, 1.3%), proliferation-associated genes (n = 15, 4.7%), and hyperdiploid (n = 77, 24.1%). Moreover, the UAMS low percentage of bone disease cluster was identified as a subcluster of the MF cluster (n = 15, 4.7%). One subgroup (n = 39, 12.2%) showed a myeloid signature. Three novel subgroups were defined, including a subgroup of 37 patients (11.6%) characterized by high expression of genes involved in the nuclear factor kappa light-chain-enhancer of activated B cells pathway, which include TNFAIP3 and CD40. Another subgroup of 22 patients (6.9%) was characterized by distinct overexpression of cancer testis antigens without overexpression of proliferation genes. The third novel cluster of 9 patients (2.8%) showed up-regulation of protein tyrosine phosphatases PRL-3 and PTPRZ1 as well as SOCS3. To conclude, in addition to 7 clusters described in the UAMS classification, we identified 3 novel subsets of multiple myeloma that may represent unique diagnostic entities.

Blood ◽  
2007 ◽  
Vol 110 (11) ◽  
pp. 1494-1494
Author(s):  
Abderrahman Abdelkefi ◽  
John de Vos ◽  
Said Assou ◽  
Tarek Ben Othman ◽  
Jean-Francois Rossi ◽  
...  

Abstract Background: Thalidomide which represents an effective treatment strategy for relapsed/refractory multiple myeloma, actually represents a standard of care also for newly diagnosed multiple myeloma patients. Methods: In the present study, we adopted a gene expression profiling (GEP) strategy in an attempt to predict response (> 50% reduction in serum M protein) to primary therapy with thalidomide-dexamethasone for newly diagnosed multiple myeloma. Plasma cells (CD138+) were purified from bone marrow aspirates from 17 patients at diagnosis, before initiation of treatment with thalidomide-dexamethasone. GEP was performed using the Affymetrix U133 Plus_2 microarray platform. The Affymetrix output (CEL files) was imported into Genespring 7.3 (Agilent technologies) microarray analysis software, where data files were normalized across chips using GCRMA and to the 50th percentile, followed by per gene normalization to median. Criteria of response were those established by Bladè et al. Results: After sufficient follow-up, responders (n=9) and nonresponders (n=8) were identified, and gene expression differences in baselines samples were examined. Of the 11000 genes surveyed, Wilcoxon rank sum test identified 149 genes that distinguished response from non response. A multivariate step-wise discriminant analysis (MSDA) revealed that 14 of the 149 genes could be used in a response predictor model (see table). Of interest, the gene list encompasses WXSC1, known to be involved in the chromosomal translocation t(4;14) (p16.3;q32.3) in multiple myeloma. Conclusion: These results could be the first step to adopt microfluidic cards, in an attempt to select at diagnosis patients who will respond favourably to a particular treatment strategy. List of 14 genes able to predict response to primary therapy with thalidomide-dexamethasone for newly diagnosed multiple myeloma. Gene ID Gene Name Chromosomal location 212771_at C10orf38 10p13 229874_x_at LOC400741 1p36.13 219690_at U2AF1L4 19q13.12 202207_at ARL7 2q37.1 243819_at GNG2 14q21 203753_at TCF4 18q21.1 235400_at FCRLM1 1q23.3 211474_s_at SERPINB6 6p25 226785_at ATP11C Xq27.1 215440_s_at BEXL1 Xq22.1–q22.3 209054_s_at WXSC1 4p16.3 227168_at FLJ25967 22p12.1 213355_at ST3GAL6 3q12.1 223218_s_at NFKBIZ 3p12–q12


Blood ◽  
2004 ◽  
Vol 104 (11) ◽  
pp. 73-73 ◽  
Author(s):  
Dirk Hose ◽  
Jean-Francois Rossi ◽  
Carina Ittrich ◽  
John deVos ◽  
Axel Benner ◽  
...  

Abstract AIM was to establish a new molecular classification of Multiple Myeloma (MM) based on changes in global gene expression attributable to cytogenetic aberrations detected by interphase FISH (iFISH) in order to (i) predict event free survival (EFS) and (ii) investigate differentially expressed genes as basis for a group specific and risk adapted therapy. PATIENTS AND METHODS. Bone marrow aspirates of 105 newly diagnosed MM-patients (65 trial (TG) / 40 independent validation group (VG)) and 7 normal donors (ND) were CD138-purified by magnetic activated cell sorting. RNA was in-vitro transcribed and hybridised to Affymetrix HG U133 A+B GeneChip (TG) and HG U133 2.0 plus arrays (VG). CCND1 and CCND2 expression was verified by real time RT-PCR. iFISH was performed on purified MM-cells using probes for chromosomes 11q23, 11q13, 13q14, 17p13 and the IgH-translocations t(4;14) and t(11;14). Expression data were normalised (Bioconductor package gcrma) and nearest shrunken centroids (NSC) applied to calculate and cross validate a predictor on 40 patients of the TG with a comprehensive iFISH panel available combined with CCND overexpression. Differentially expressed genes were identified using empirical Bayes statistics for pairwise comparison. RESULTS. Overexpression of a D-type cyclin (D1 or D2) was found in 61/65 patients with MM compared to ND. CCND3 overexpression only appeared concomitantly with CCND2 overexpression. Four groups could be distinguished: (1.1) CCND1 (11q13) overexpression and trisomy 11q13, (1.2) CCND1 overexpression and translocations involving 11q13 i.e. t(11;14), (2.1) CCND2 overexpression without 11q13+, t(11;14), t(4;14), (2.2) CCND2 overexpression with t(4;14) and FGFR3 upregulation. A predictor of 6 to 566 genes correctly classifies all 40 patients of the TG (estimated cross validated error rate 0%). An independent VG of 40 patients was used. Genes with highest scores in NSC are: (1.1) CCND1, ribosomal proteins (e.g. RPL 28, 29), GPX1, CCRL2, (1.2) CCND1, TGIF, and NCAM (non-overexpression), (2.1) CCND2, (2.2) FGFR3, WHSC1, CCND2, IRTA2, SELL, and MAGED4. Distribution of clinical parameters (i.e. β2M, Durie Salmon stages, ISS) was not significantly different between the groups. The distribution of del(13)(q14q14) was (1.1) 31.5%, (1.2) 37.5%, (2.1) 37.5% and (2.2) 100%. (p<0.01). I.e. HGF, DKK1, VCAM, CD163 are differentially expressed between all 4 groups and ND (adjusted p<0.001). The groups defined by the predictor show a significantly different EFS after autologous stem cell transplantation according to the GMMG-HD3 protocol (median: (1.1) 18 / (1.2) not reached (no event) / (2.1) 22 / (2.2) 6 months; log-rank-test: p=0.004). CONCLUSION. CCND1 or CCND2 overexpression is nearly ubiquitous in MM and attributable to defined cytogenetic aberrations. Gene expression and iFISH allow a molecular classification of MM which can be predicted by gene expression profiling alone. Groups in the classification show a distinctive pattern in gene expression as well as a different EFS interpretable as risk stratification and indicator of therapeutic targets.


Blood ◽  
2006 ◽  
Vol 108 (11) ◽  
pp. 114-114
Author(s):  
Guido Tricot ◽  
Fenghuang Zhan ◽  
Bart Barlogie ◽  
Yongsheng Huang ◽  
Jeffrey Sawyer ◽  
...  

Abstract The International Staging System (ISS), based on B2-microglobulin and albumin levels at the time of diagnosis, has now generally been adopted as a new prognostic classification system for multiple myeloma (MM). While readily and widely applicable, ISS does not account for genetic disease features, such as metaphase (CA) and interphase fluorescence in situ hybridization (FISH) cytogenetic abnormalities, which when examined in the context of standard prognostic variables, confer higher hazards of relapse and disease-related death. Recently, gene expression profiling (GEP) uncovered the major prognostic significance for outcome of high expression of CKS1B, a gene mapping to an amplicon at chromosome 1q21. We have performed a comprehensive study of CA, FISH, GEP and ISS staging in 351 newly diagnosed MM patients, treated uniformly on Total Therapy 2. We have analyzed outcome based on a combination of high CKS1B by GEP and CA. GEP-based t(11;14) was prognostically favorable, irrespective of expression of CKS1B and, therefore, was removed from the group of patients with high CKS1B expression. After this adjustment, with the combination of both CA and high CKS1B (approximately 10% of all patients) conferred a very poor outcome with only 24% and 40% of such patients being event-free and/surviving at 3 years, compared with 72% and 84% for the others (p values : &lt;.0001). Such patients fared poorly, irrespective of their ISS stage. Similar prognostic information could be gained by combining CA with FISH-defined amplification of 1q21 and t(11;14). Because of their major prognostic impact, all newly diagnosed patients should be tested for these genetic markers. Novel treatment modalities are justified in the small subgroup of such poor prognosis patients, since they derive only a minor benefit from advances in MM therapy. CKS1B Q4 + CA (with no CCND1) vs. Others CKS1B Q4 + CA (with no CCND1) vs. Others


Blood ◽  
2016 ◽  
Vol 128 (22) ◽  
pp. 3264-3264 ◽  
Author(s):  
Ryan K Van Laar ◽  
Ivan Borrelo ◽  
David Jabalayan ◽  
Ruben Niesvizky ◽  
Aga Zielinski ◽  
...  

Abstract Background: There is a global consensus that multiple myeloma patients with high-risk disease require additional monitoring and therapy compared to low/standard risk patients in order to maximize their chances of survival. Current diagnostic guidelines recommend FISH-based assessment of chromosomal aberrations to determine risk status (i.e. t(14;20), t(14;16), t(4;14) and/or Del17p), however, studies show FISH for MM may have a 20-30% QNS rate and is up to 15% discordant between laboratories, even when starting from isolated plasma cells. In this study we demonstrate that MyPRS gene expression profiling reproduces the key high risk translocations for MM risk stratification, in addition to having other significant advantages. Methods: Reproducibility studies show that MyPRS results are less than 1% discordant starting from isolated plasma cells and return successful results in up to 95% of cases. 270 MM patients from Johns Hopkins University (MD) and Weill Cornell Medicine (NY) had both FISH and MyPRS gene expression profiling performed between 2012 and 2016 using standard and previously published methodology, respectively. Results: Retrospective review of the matched FISH and MyPRS results showed: 25/28 (89%) patients wish FISH-identified t(4;14) were classified as MMSET (MS) subtype. 10/10 (100%) patients with t(14;16) or t(14;20) were classified as MAF-like (MF) subtype 62/67 (93%) patients with t(11;14) were assigned to the Cyclin D (1 or 2) subtype. Patients with FISH hyperdiploidy status were classified as the Hyperdiploid (HY) subtype or had multiple gains detected by the separate MyPRS Virtual Karyotype (VK) algorithm, included in MyPRS. TP53del was seen in patients with multiple molecular subtypes, predominantly Proliferation (PR) and MMSET (MS). Assessment of TP53 function by gene expression is a more clinically relevant prognostic marker than TP53del, as dysregulation of the tumor suppressor is affected by mutations as well as deletions. Analysis of the TP53 expression in the 39 patients with delTP53 showed a statistically significant difference, compared to patients without this deletion (P<0.0001). Conclusion: Gene expression profiling is a superior and more reliable method for determining an individual patients' prognostic risk status. The molecular subtypes of MM, as reported by Signal Genetics MyPRS assay, are driven by large-scale changes in gene expression caused by or closely associated with chromosomal changes, including translocations. Physicians who are managing myeloma patients and wishing to base their assessment of risk on R-ISS or mSMART Guidelines may obtain the required data points from either FISH or MyPRS, with the latter offering lower QNS rates, higher reproducibility, assessment of a larger number of cells and a substantially lower price point ($5,480 vs. $1,912; 2016 CMS data). A larger cohort study is now underway to further validate these observations. Figure GEP-based TP53 expression in patients with and without Del17p. P<0.0001 Figure. GEP-based TP53 expression in patients with and without Del17p. P<0.0001 Disclosures Van Laar: Signal Genetics, Inc.: Employment. Borrelo:Sidney Kimmel Cancer Institute: Employment. Jabalayan:Weill Cornell Medical Center: Employment. Niesvizky:Celgene: Consultancy, Research Funding, Speakers Bureau; Takeda: Consultancy, Research Funding, Speakers Bureau; Onyx: Consultancy, Research Funding, Speakers Bureau. Zielinski:Signal Genetics, Inc.: Employment. Leigh:Signal Genetics, Inc.: Employment. Brown:Signal Genetics, Inc.: Employment. Bender:Signal Genetics, Inc.: Employment.


Blood ◽  
2008 ◽  
Vol 112 (11) ◽  
pp. 1705-1705
Author(s):  
Ricky D Edmondson ◽  
Sheeno P Thyparambil ◽  
Veronica MacLeod ◽  
Bart Barlogie ◽  
John D. Shaughnessy

Abstract Although melphalan-based autologous stem cell transplantation has improved prognosis for patients diagnosed with Multiple Myeloma, survival varies from a few months to more than 15 years with an individual’s risk not accurately predicted with standard prognostic variables. Correlating genome-wide mRNA expression profiles in purified myeloma cells with outcome, we recently showed that that the differential expression of 70 genes could identify patients at high risk for early disease related death [1]. The utility of a high throughput proteomics platform in the analysis of clinical samples has great potential but as of yet none have been firmly established. Herein, we describe the use of such a platform and its utility in stratifying patients with Multiple Myeloma in terms of high and low risk disease. Preliminary analysis indicates that the proteomics data can separate the patients into risk groups, although the proteins responsible for the assignment are not identical to the 70 genes identified in the gene expression profiling experiments. In addition to the proteomic analysis of plasma cells enriched using anti-CD138 immunomagnetic beads from mononuclear cell fractions of bone marrow aspirates from newly diagnosed myeloma patients; we have performed (in triplicate) LCMS profiling on plasma cells from 30 patients isolated prior to and 48 hours after a single test-dose application of bortezomib at 1.0mg/m2. An aliquot of 100,000 plasma cells was enzymatically digested with trypsin and a fraction (~5,000 cells) analyzed using our proteomics platform (an Eksigent nanoHPLC coupled to a ThermoElectron LTQ-Orbitrap with data analyzed using the Elucidator software package from Rosetta Biosoftware). The correlation of the proteomic profiles to gene expression profiles and clinical parameters will be presented. The analysis of proteins that were observed to change (p&lt;0.01) in abundance after the single agent dose of the proteasome inhibitor bortezomib yielded an unanticipated finding; the abundance of 30 proteins associated with the proteasome were observed to increase in a subset of patients. The majority of the patients with the increased levels of proteasome related proteins are predicted by GEP to have high risk disease. The proteomic data will be discussed in terms of its utility in the identification of activated pathways as well as in the development of a prognostic indicator as was achieved using gene expression profiling.


Blood ◽  
2013 ◽  
Vol 122 (21) ◽  
pp. 3123-3123
Author(s):  
Bart Barlogie ◽  
Emily Hansen ◽  
Sarah Waheed ◽  
Jameel Muzaffar ◽  
Monica Grazziutti ◽  
...  

Abstract Intra-tumoral heterogeneity (ITH) is increasingly viewed as the Achilles heel of treatment failure in malignant disease including multiple myeloma (MM). Most MM patients harbor focal lesions (FL) that are recognized on MRI long before bone destruction is detectable by conventional X-ray examination. Serial MRI examinations show that eventually 60% of patients will achieve resolution of FL (MRI-CR). However, this will lag behind the onset of a clinical CR by 18 to 24 months, thus attesting to the biological differences between FL and diffuse MM growth patterns. Consequently, we performed concurrent gene expression profiling (GEP) analyses of plasma cells (PC) from both random bone marrow (RBM) via iliac crest and FL. Our primary aims were to first compare the molecular profiles of FL vs. RBM, second to determine if ITH existed (as defined molecular subgroup and risk), and finally to investigate if the bone marrow micro-environment (ME) contained a biologically interesting signature. A total of 176 patients were available for this study with a breakdown of: TT3 (n=23), TT4 for low-risk (n=131) and TT5 for high-risk MM (n=22). Regarding the molecular analyses of PCs, GEP-based risk (GEP-70, GEP-5) and molecular subgroup correspondence were examined for commonalties and differences between RBM and FL. A “filtering” approach for ME genes was also developed and bone marrow biopsy (BMBx) GEP data derived from this method is under analysis. PC risk correspondence between FL and RBM was 86% for GEP70 and 88% for the GEP5 model. Additionally, 82% had a molecular subgroup concordance, however, they did differ among subgroups (p=0.020) by Fisher's Exact Test. A lower concordance was noted in the CD2, LB, and PR subgroups (67%, 69%, 73%, respectively). GEP70 and GEP5 risk concordance between RBM and FL samples by molecular subgroup was also examined. The overall correlation coefficients were 0.619 (GEP70) and 0.597 (GEP5). The best correspondence was noted for CD1, MF and PR subgroups especially for the GEP5 model. HY, LB and MS showed intermediate correlations, while CD2 fared worst with values of only 0.322 for GEP70 and 0.267 for GEP5 model. Figure 1 portrays these data in more detail for the GEP70 and GEP5 models. Good correlations were noted between RBM and FL based risk scores in case of molecular subgroup concordance (left panels) in both GEP5 and GEP70 risk models, whereas considerable scatter existed in case of subgroup discordance (right panels). The clinical implications in TT4 regarding RBM and FL derived risk and molecular subgroup information, viewed in the context of standard prognostic baseline variables are portrayed in Table 1. High B2M levels at both cut-points imparted inferior OS and PFS as did low hemoglobin. Although present in 42% of patients, cytogenetic abnormalities (CA) did not affect outcomes. FL-based GEP5-defined high-risk designation conferred poor OS and PFS. B2M>5.5mg/L and FL-derived GEP5 high-risk MM, pertaining to 29% and 11% of patients, survived the multivariate model for both OS and PFS. Next, in examining PC-GEP differences among RBM and FL sites, 199 gene probes were identified with a false discovery rate (FDR) of 1x10-6. Additionally, 55 of the 199 belong to four molecular networks of inter related genes associated with: lipid metabolism, cellular movement, growth and proliferation, and cell-to-cell interactions. Multivariate analysis identified the GEP5 high risk designation of focal lesion PCs to be significantly prognostic with a HR=3.73 (p=0.023).Table 1Cox regression analysis of variables linked to overall and progression-free survival in TT4.Overall SurvivalProgression-Free SurvivalVariablen/N (%)HR (95% CI)P-valueHR (95% CI)P-valueMultivariateB2M > 5.5 mg/L38/130 (29%)3.71 (1.49, 9.22)0.0053.84 (1.58, 9.31)0.003FL GEP5 High Risk14/130 (11%)3.68 (1.19, 11.41)0.0243.73 (1.20, 11.62)0.023HR- Hazard Ratio, 95% CI- 95% Confidence Interval, P-value from Wald Chi-Square Test in Cox RegressionNS2- Multivariate results not statistically significant at 0.05 level. All univariate p-values reported regardless of significance.Multivariate model uses stepwise selection with entry level 0.1 and variable remains if meets the 0.05 level.A multivariate p-value greater than 0.05 indicates variable forced into model with significant variables chosen using stepwise selection. Disclosures: No relevant conflicts of interest to declare.


Sign in / Sign up

Export Citation Format

Share Document