Bryant, Barb

Abstract Introduction: We investigated tumor gene expression profiling (PGx) to identify predictive classifiers of response and elucidate the mechanism of efficacy of proteasome inhibition in multiple myeloma (MM) by bortezomib, a first-in-class proteasome inhibitor. Analysis of 44 samples from the SUMMIT phase 2 trial identified predictive classifiers of response in patients (pts) with refractory MM (ASH 2002, Abstr 1519). We now report on the PGx analysis of the SUMMIT-derived classifier in an independent sample set from the APEX phase 3 trial, a multicenter randomized study of 669 pts that demonstrated the superiority of bortezomib to high-dose dex for relapsed MM following 1 to 3 prior therapies (Richardson NEJM 2005). As well, we report on a new classifier of response built by combining samples from SUMMIT and MLN34101-040 (040) (a study of bortezomib in pts who failed dex in APEX or were not eligible for APEX) and tested in the APEX samples. Methods: 93 participating clinical centers collected pretreatment bone marrow aspirates and enriched tumor cells on site before centralized RNA extraction and microarray analysis. Approximately 85% of pts consented to whole-genome analysis. Quality control for RNA integrity and quantity, microarray data quality, and tumor cell purity resulted in 129 patient samples from APEX evaluable for PGx; 67 were from the bortezomib arm and 62 from the dex arm. 48 samples were evaluable from 040. Patient characteristics of this subset were similar to those of the overall study population. Predictive classifiers were identified and tested via separate training and test datasets. Results: Consistent with the SUMMIT analysis, the APEX and 040 datasets exhibit biologic features expected of MM, including the overexpression of translocation target genes such as cyclin D1, MAF, and MMSET/FGFR3. Classifiers built using only the SUMMIT samples predict the outcome of pts in the APEX bortezomib arm with an overall accuracy of 51% (sensitivity [Sn] = 49%, specificity [Sp] = 55%; P = .81) and in the dex arm with an overall accuracy of 56% (Sn = 64%, Sp = 50%; P = .31). If the classifier was built using data from both the SUMMIT and 040 trials (87 total samples), the prediction of response in the APEX bortezomib arm was statistically significant (66% accuracy, Sn = 89%, Sp = 42%; P = .006). The same classifier only predicts the response of dex-treated pts with an overall accuracy of 46% (Sn = 72%, Sp = 28%; P = .88), suggesting specificity for bortezomib. Many of the overexpressed genes associated with bortezomib insensitivity relate to protein synthesis (EIF4B, RPL5) and RNA binding (SFRS10, CPSF6). Future Directions: This work highlights practical and scientific issues that inform future efforts to implement PGx in clinical trials. Predictive accuracy may be influenced by the clinical end point, limited sample size, and/or technical considerations associated with distinct trials and genomic data. The influence of study differences (relating to pts’ time from diagnosis, number of prior therapies, molecular subtypes, or hybridization characteristics) on the ability to build classifiers and identify biologic pathways related to response are being examined; these and potential insights on drug mechanisms will be presented.

Abstract Introduction: Clinical development of bortezomib in multiple myeloma (MM) included pharmacogenomic gene expression profiling (GEP) of pre-treatment tumor biopsies in order to evaluate the feasibility of developing predictive classifiers in prospective clinical trials and to better define the biological pathways related to drug sensitivity. By partitioning the samples into train (n=101) and test (n=80) datasets, statistically significant classifiers of response and survival were identified and then validated [Mulligan G, et al: Blood. 2007; 109:3177–88]. Here we describe specific genes and biological pathways associated with bortezomib efficacy in relapsed MM. Methods: Using tumor GEP data from 188 patients treated with single-agent bortezomib we identified genes associated with overall survival (OS), as well as genes associated with response (R) or progressive disease (PD) after bortezomib therapy. These clinical endpoints were also analyzed with biological pathway algorithms including GSEA [Subramanian A, et al: Proc Natl Acad Sci2005; 102:15545–50] Ingenuity Pathway Analysis and Gene Ontology/GATHER. These results were then compared to GEP data from a patient subset (n=76) treated with single-agent dexamethasone and also cell lines selected for differential sensitivity to bortezomib. Results: A comparison of R and PD highlighted numerous differentially expressed genes; approximately 200 were over-expressed in R and 500 were over-expressed in PD (t-test <0.01). R-associated genes include signaling molecules like TRADD, CFLAR and TANK, while PD-associated genes included oncogenes (NRAS), apoptotic regulators (DAP3) and cancer antigens (CTAG). Pathway tools helped to integrate and extend these observations, consistently identifying NFkB activation and overexpression of adhesion pathways as R-associated. Overexpression of protein synthesis and mitochondrial pathways were PD-associated. The overexpression of pathways such as NFkB, adhesion and cytokines appear linked specifically with bortezomib and not dexamethasone sensitivity. These approaches were also applied to patient survival. In contrast to analyses of single genes, the pathway tools highlighted links between response and survival, including adhesion and inflammatory cytokine pathway expression. Conclusions: This study highlights the diversity of biological pathways expressed in MM and the utility of clinical genomics to associate pathways with activity of single-agent bortezomib.

Abstract Multiple myeloma is a disease of the antibody producing plasma cells in the bone marrow. It has been characterized by frequent genomic alterations including gains and losses of chromosomes. Some of these alterations result in poor clinical prognosis. Bortezomib, a first in class proteasome inhibitor, is approved for the treatment of relapsed multiple myeloma. To better understand how the genome wide changes observed in myeloma relate to prognosis and treatment response, we have used single nucleotide polymorphism (SNP) array technology to assess genotypes and DNA copy number at up to 100,000 genomic loci. DNA from 147 bone marrow tumor biopsies collected in multi-center phase II and III clinical trials of relapsed multiple myeloma patients prior to treatment with Bortezomib were hybridized on SNP arrays to assess genomic aberrations. Commonly seen genomic alterations in myeloma were observed, including deletions of chromosome 13 (37%), 1p (25%), 6q (23%), and amplifications on 1q (35%). We found that 35% of this relapsed myeloma population exhibited a hyperdiploid genome characterized by gains of chromosomes 3, 5, 7, 9, 11, 15, and 19. Other notable genomic deletions included 17p (22%), 8p (27%) and 10p (16%). A correlation of response to bortezomib therapy in patients with differential gains and losses at the chromosomal level is ongoing. Consistent with previous reports, response to bortezomib was observed in patients with poor prognostic chr 13 deletions. A more in-depth analysis of specific loci associated with clinical outcome will be discussed. In summary, this data will help further define genomic aberrations in multiple myeloma and their association with treatment outcome.