Abstract
Variations within proteasome β (PSMB) genes, which encode the β subunits of the 20S proteasome, may affect proteasome function, assembly, and/or binding of proteasome inhibitors. To investigate the potential association between PSMB gene variants and treatment-emergent resistance to bortezomib and/or long-term outcomes, in the present study, PSMB gene sequence variation was characterized in tumor DNA samples from patients who participated in the phase 3 Assessment of Proteasome Inhibition for Extending Remissions (APEX) study of bortezomib versus high-dose dexamethasone for treatment of relapsed multiple myeloma. Twelve new PSMB variants were identified. No associations were found between PSMB single nucleotide polymorphism genotype frequency and clinical response to bortezomib or dexamethasone treatment or between PSMB single nucleotide polymorphism allelic frequency and pooled overall survival or time to progression. Although specific PSMB5 variants have been identified previously in preclinical models of bortezomib resistance, these variants were not detected in patient tumor samples collected after clinical relapse from bortezomib, which suggests that alternative mechanisms underlie bortezomib insensitivity. This study is registered at www.clinicaltrials.gov as NCT00048230.
Introduction
The 20S core of the 26S proteasome degrades polyubiquitinated intracellular proteins1 and is composed of 4 stacked rings,1-4 each with 7 α and 7 β subunits. Three constitutive proteasome β (PSMB) subunits, β5, β2, and β1 (encoded by the PSMB5, PSMB7, and PSMB6 genes, respectively),5 are responsible for chymotrypsin-like, trypsin-like, and post-glutamyl peptide hydrolyzing activities, respectively.6-8 On IFN-γ stimulation, constitutive β subunits are replaced by the IFN-inducible subunits β5i, β2i, and β1i (encoded by PSMB8, PSMB10, and PSMB9 genes, respectively)5 to form the immunoproteasome.6,9 Variations in PSMB subunits could potentially affect proteasome structure, assembly, function, and/or binding of proteasome inhibitors.
Bortezomib (trade name Velcade; Millennium Pharmaceuticals) is approved in the United States and Europe for the treatment of patients with multiple myeloma (MM),10,11 and in the United States for patients with relapsed mantle cell lymphoma.10 Bortezomib selectively binds to the β5 subunit, leading to full inhibition of ubiquitinated protein hydrolysis.7 In addition, bortezomib interacts with the β1 subunit8 and, when bound to the β5 subunit in the chymotryptic catalytic site, is in close proximity to the β6 subunit.12 Several studies have shown that PSMB5 variants can arise in vitro when tumor cell lines are cultured with bortezomib13-18 ; it remains unclear whether this mechanism is relevant for bortezomib resistance in the clinical setting.
The present study addressed whether variations in PSMB genes affect treatment-emergent resistance in bortezomib-treated MM patients or long-term outcome in MM patients. Sequence variation was characterized in coding regions of PSMB genes in pre- and posttreatment samples from patients who participated in the phase 3 Assessment of Proteasome Inhibition for Extending Remissions (APEX) trial of single-agent bortezomib versus high-dose dexamethasone (Dex) for the treatment of relapsed MM.19
Study design
Review boards at all participating institutions approved the (APEX) study,19 and BM aspirates were obtained from consenting patients in accordance with the Declaration of Helsinki during the APEX trial. Tumor cells were purified and frozen for nucleic acid isolation as described previously.20 Matching germline DNA samples were not collected. DNA samples were amplified using the QIAGEN REPLI-g whole genome amplification kit and used for PCR reactions with primers for coding regions of the PSMB1, PSMB5, PSMB6, PSMB8, PSMB9, and PSMB10 genes (supplemental Table 1, available on the Blood Web site; see the Supplemental Materials link at the top of the online article). Bidirectional DNA sequencing was performed to increase confidence in genetic variant identification. Sequence data were analyzed using Mutation Surveyor Version 2.61 (SoftGenetics) and Sequencher Version 4.8 (Gene Codes Corporation) software.
Allelic and genotype frequencies were typically compared with a weighted average of European population data from the National Center for Biotechnology Information single nucleotide polymorphism (SNP) database (dbSNP) using the Fisher exact test. P values were adjusted for multiplicity using the false discovery rate adjustment.21 The Fisher exact test and log-rank test were used to compare genotype frequencies with respect to clinical response and allelic frequencies with respect to pooled overall survival (OS) and time to progression (TTP), respectively. Full methodological details are provided in supplemental Methods.
Results and discussion
SNP frequency analysis
A total of 76 DNA samples were of adequate yield and quality for sequencing, including 47 (25 bortezomib-treated and 22 Dex-treated) pretreatment samples and 29 (16 bortezomib-treated and 13 Dex-treated) posttreatment samples. Paired pre- and posttreatment samples were available from 6 patients (3 bortezomib-treated and 3 Dex-treated). The dataset size limits formal statistical analyses of SNP associations with MM incidence or baseline characteristics; however, it is unique with respect to data regarding sensitivity to single-agent bortezomib and sampling before and after therapy. Allelic and genotype frequency of nonsynonymous SNPs in pre- and posttreatment MM samples did not differ significantly from population frequency data in dbSNP (Table 1 and supplemental Table 2), suggesting that nonsynonymous variants in PSMB are not specifically selected in MM. No unique nonsynonymous substitutions were observed in posttreatment samples. Further, recurrent variants (eg, in PSMB8) occurred at similar frequencies in pre- and posttreatment samples, suggesting that PSMB variation does not arise from selective pressure exerted during treatment and that these variants are more likely to represent naturally occurring germline SNPs. Supporting this, 3 SNPs detected in this study were also detected in a previous study of gene variants in MM,22 and 2 (rs2304974, rs2304975) were detected at similar frequencies in both studies.
Gene . | SNP ID . | Amino acids . | Alleles . | NCBI frequency . | Observed frequency . | FDR P . |
---|---|---|---|---|---|---|
Pretreatment samples | ||||||
PSMB1 | rs12717 | P11A | C/G | 0.61/0.39 | 0.58/0.42 | .813 |
PSMB5 | rs11543947 | R24C | C/T | NA* | 0.95/0.05 | NA |
PSMB6 | rs3169950 | A21A | G/A | 0.57/0.43 | 0.67/0.33 | .386 |
rs11552523 | T56T | T/C | 0.99/0.02† | 0.98/0.02 | > .999 | |
rs2304974 | P107A | C/G | 0.91/0.09† | 0.94/0.06 | .746 | |
rs2304975 | S115S | C/T | 0.64/0.36† | 0.77/0.23 | .081 | |
rs7468 | Y177Y | C/T | 0.91/0.09‡ | 0.93/0.07 | .740 | |
PSMB8 | rs17220206 | S70T | T/A | 0.97/0.03† | 1.00/0.00 | .573 |
PSMB9 | rs35100697 | G9E | G/A | 0.97/0.03† | 1.00/0.00 | > .999 |
rs241419 | V32I | G/A | 0.99/0.01 | 0.98/0.02 | .777 | |
rs17587 | R60H | G/A | 0.70/0.30 | 0.77/0.23 | .635 | |
rs20547 | A161A | A/G | 0.99/0.01 | 0.94/0.06 | .081 | |
rs17213861 | R173C | C/T | 0.98/0.02† | 1.00/0.00 | .740 | |
PSMB10 | rs20549 | L107L | T/C | 0.83/0.17 | 0.72/0.28 | .168 |
rs14178 | G151G | T/C | 0.95/0.05 | 0.87/0.13 | .057 | |
Posttreatment samples | ||||||
PSMB1 | rs12717 | P11A | C/G | 0.61/0.39 | 0.63/0.38 | > .999 |
PSMB5 | rs11543947 | R24C | C/T | NA* | 0.83/0.17 | NA |
PSMB6 | rs3169950 | A21A | G/A | 0.57/0.43 | 0.41/0.59 | .499 |
rs11552523 | T56T | T/C | 0.99/0.02† | 0.97/0.03 | .701 | |
rs2304974 | P107A | C/G | 0.91/0.09† | 0.97/0.03 | .701 | |
rs2304975 | S115S | C/T | 0.64/0.36† | 0.97/0.03 | < .001 | |
rs7468 | Y177Y | C/T | 0.91/0.09‡ | 0.96/0.04 | .701 | |
PSMB8 | rs17220206 | S70T | T/A | 0.97/0.03† | 1.00/0.00 | .701 |
PSMB9 | rs35100697 | G9E | G/A | 0.97/0.03† | 1.00/0.00 | > .999 |
rs241419 | V32I | G/A | 0.99/0.01 | 0.98/0.02 | .701 | |
rs17587 | R60H | G/A | 0.70/0.30 | 0.79/0.21 | .701 | |
rs20547 | A161A | A/G | 0.99/0.01 | 0.88/0.12 | .021 | |
rs17213861 | R173C | C/T | 0.98/0.02† | 1.00/0.00 | > .999 | |
PSMB10 | rs20549 | L107L | T/C | 0.83/0.17 | 0.72/0.28 | .499 |
rs14178 | G151G | T/C | 0.95/0.05 | 0.91/0.09 | .612 |
Gene . | SNP ID . | Amino acids . | Alleles . | NCBI frequency . | Observed frequency . | FDR P . |
---|---|---|---|---|---|---|
Pretreatment samples | ||||||
PSMB1 | rs12717 | P11A | C/G | 0.61/0.39 | 0.58/0.42 | .813 |
PSMB5 | rs11543947 | R24C | C/T | NA* | 0.95/0.05 | NA |
PSMB6 | rs3169950 | A21A | G/A | 0.57/0.43 | 0.67/0.33 | .386 |
rs11552523 | T56T | T/C | 0.99/0.02† | 0.98/0.02 | > .999 | |
rs2304974 | P107A | C/G | 0.91/0.09† | 0.94/0.06 | .746 | |
rs2304975 | S115S | C/T | 0.64/0.36† | 0.77/0.23 | .081 | |
rs7468 | Y177Y | C/T | 0.91/0.09‡ | 0.93/0.07 | .740 | |
PSMB8 | rs17220206 | S70T | T/A | 0.97/0.03† | 1.00/0.00 | .573 |
PSMB9 | rs35100697 | G9E | G/A | 0.97/0.03† | 1.00/0.00 | > .999 |
rs241419 | V32I | G/A | 0.99/0.01 | 0.98/0.02 | .777 | |
rs17587 | R60H | G/A | 0.70/0.30 | 0.77/0.23 | .635 | |
rs20547 | A161A | A/G | 0.99/0.01 | 0.94/0.06 | .081 | |
rs17213861 | R173C | C/T | 0.98/0.02† | 1.00/0.00 | .740 | |
PSMB10 | rs20549 | L107L | T/C | 0.83/0.17 | 0.72/0.28 | .168 |
rs14178 | G151G | T/C | 0.95/0.05 | 0.87/0.13 | .057 | |
Posttreatment samples | ||||||
PSMB1 | rs12717 | P11A | C/G | 0.61/0.39 | 0.63/0.38 | > .999 |
PSMB5 | rs11543947 | R24C | C/T | NA* | 0.83/0.17 | NA |
PSMB6 | rs3169950 | A21A | G/A | 0.57/0.43 | 0.41/0.59 | .499 |
rs11552523 | T56T | T/C | 0.99/0.02† | 0.97/0.03 | .701 | |
rs2304974 | P107A | C/G | 0.91/0.09† | 0.97/0.03 | .701 | |
rs2304975 | S115S | C/T | 0.64/0.36† | 0.97/0.03 | < .001 | |
rs7468 | Y177Y | C/T | 0.91/0.09‡ | 0.96/0.04 | .701 | |
PSMB8 | rs17220206 | S70T | T/A | 0.97/0.03† | 1.00/0.00 | .701 |
PSMB9 | rs35100697 | G9E | G/A | 0.97/0.03† | 1.00/0.00 | > .999 |
rs241419 | V32I | G/A | 0.99/0.01 | 0.98/0.02 | .701 | |
rs17587 | R60H | G/A | 0.70/0.30 | 0.79/0.21 | .701 | |
rs20547 | A161A | A/G | 0.99/0.01 | 0.88/0.12 | .021 | |
rs17213861 | R173C | C/T | 0.98/0.02† | 1.00/0.00 | > .999 | |
PSMB10 | rs20549 | L107L | T/C | 0.83/0.17 | 0.72/0.28 | .499 |
rs14178 | G151G | T/C | 0.95/0.05 | 0.91/0.09 | .612 |
All frequencies totaled 1.00, although some appear not to due to being rounded to the nearest hundredth unit.
FDR indicates false discovery rate; and NA, not applicable.
NCBI frequency data are inconclusive. Bold font indicates statistical significance (FDR P < .05).
NCBI European frequency data unavailable (listed frequency is weighted average from available NCBI population data).
NCBI European frequency data unavailable (different major/minor allele depending on population).
Novel variant identification
Twelve novel, low-frequency PSMB variants were identified that were not listed in dbSNP (Table 2); 10 were nonsynonymous variants. Two variants were located in PSMB5; one, a C/G substitution resulting in a S112R change in PSMB5 approximately 7 angstroms away from the bortezomib binding pocket, was observed in 1 pretreatment sample from a patient who achieved a partial response (PR) to bortezomib.
Gene . | Nucleotide change . | Amino acid change . | Samples in which variant identified, treatment (clinical response) . |
---|---|---|---|
PSMB1 | A/G | I50V | 1 pretreatment: bortezomib (PR) |
PSMB5 | C/G | S112R | 1 pretreatment: bortezomib (PR) |
C/T | A218V | 1 pretreatment: Dex (SD) | |
PSMB6 | C/T | A195V | 1 pretreatment: Dex (SD) |
G/A | A239T | 1 pretreatment: Dex (SD) | |
PSMB8 | G/A | G8R | 1 pretreatment: bortezomib (PD); 1 posttreatment: bortezomib (PR) |
T/C | A43A | 1 pretreatment: bortezomib (SD); 1 posttreatment: bortezomib (NC) | |
C/A | Q49K | 10 pretreatment: 6 bortezomib (2 PR, 3 SD, 1 PD), 4 Dex (1 MR, 3 SD); 3 posttreatment: 1 bortezomib (PR), 2 Dex (NC) | |
PSMB9 | G/A | A25T | 1 pretreatment: bortezomib (PR) |
PSMB10 | C/T | P33S | 1 pretreatment: bortezomib (PR) |
C/T | P154P | 1 posttreatment: bortezomib (PR) | |
G/A | G239S | 1 pretreatment: bortezomib (SD) |
Gene . | Nucleotide change . | Amino acid change . | Samples in which variant identified, treatment (clinical response) . |
---|---|---|---|
PSMB1 | A/G | I50V | 1 pretreatment: bortezomib (PR) |
PSMB5 | C/G | S112R | 1 pretreatment: bortezomib (PR) |
C/T | A218V | 1 pretreatment: Dex (SD) | |
PSMB6 | C/T | A195V | 1 pretreatment: Dex (SD) |
G/A | A239T | 1 pretreatment: Dex (SD) | |
PSMB8 | G/A | G8R | 1 pretreatment: bortezomib (PD); 1 posttreatment: bortezomib (PR) |
T/C | A43A | 1 pretreatment: bortezomib (SD); 1 posttreatment: bortezomib (NC) | |
C/A | Q49K | 10 pretreatment: 6 bortezomib (2 PR, 3 SD, 1 PD), 4 Dex (1 MR, 3 SD); 3 posttreatment: 1 bortezomib (PR), 2 Dex (NC) | |
PSMB9 | G/A | A25T | 1 pretreatment: bortezomib (PR) |
PSMB10 | C/T | P33S | 1 pretreatment: bortezomib (PR) |
C/T | P154P | 1 posttreatment: bortezomib (PR) | |
G/A | G239S | 1 pretreatment: bortezomib (SD) |
MR indicates minimal response; NC, no change; PD, progressive disease; SD, stable disease; and Dex, high-dose dexamethasone.
The PSMB5 A108T variant, which has been linked to bortezomib resistance in vitro,13,14,17,18 was not observed in any pre- or posttreatment samples. A previous study reported the absence of this variant in the germline DNA of MM patients, but tumor DNA was not characterized either before or after bortezomib therapy.23 In the present study, tumor samples collected after bortezomib treatment were from 10 patients who were relatively insensitive to bortezomib monotherapy (best response of minimal response, stable disease, or progressive disease) and from 6 patients who achieved a confirmed PR to bortezomib and subsequently relapsed on study before sample collection (supplemental Figure 1). In these cases, treatment-emergent resistance to single-agent bortezomib was independent of variants in the proteasome genes PSMB1, PSMB5, PSMB6, PSMB8, PSMB9, and PSMB10.
The most common new variant detected, a C/A substitution resulting in a Q49K change in PSMB8, was observed in 10 pretreatment and 3 posttreatment samples: 2 other PSMB8 variants were found, 1 each in 1 pre- and 1 posttreatment sample, and 9 other variants (in PSMB1, PSMB5, PSMB6, PSMB9, and PSMB10) were detected in 1 sample each. These variants may be rare SNPs or somatic changes that occurred in the context of MM.
Correlation with clinical outcomes
There were no associations between PSMB SNP genotype frequencies in a pooled dataset of pre- and posttreatment samples and subsequent patient response to bortezomib or Dex treatment (supplemental Table 3). However, 3 significant associations were identified between SNP allelic frequencies and OS or TTP (SNPs PSMB6 rs2304975, PSMB6 rs3169950, and PSMB9 rs241419; supplemental Table 4). As a caveat, 2 SNPs, PSMB9 rs241419 and PSMB6 rs2304975, were of low frequency, being found in only 3 and 5 patients, respectively. For a third SNP (PSMB6 rs3169950), patients with the A allele appeared to have shorter OS than those with the G allele, but not shorter TTP, as may have been expected. Acknowledging the limited dataset size, these findings require confirmation in larger populations.
The PSMB1 SNP rs12717 (C/G substitution resulting in a P11A change), which was reported recently to be associated with a relative progression-free survival benefit in relapsed follicular lymphoma patients treated with bortezomib-rituximab versus rituximab,24 was not associated with OS or TTP in this pooled APEX dataset. There was no statistically significant difference in TTP for rs12717 C/G heterozygotes with bortezomib (n = 12) versus Dex treatment (n = 11; P = .077). Because of the limited size of the current dataset, additional studies are required to further evaluate the PSMB1 rs12717 SNP in MM to enable analyses of SNP combinations that may contribute to drug sensitivity.
In summary, in the present study, no unique PSMB5 variants were detected in patient tumor samples collected after bortezomib treatment, including specimens from patients who were initially sensitive to bortezomib (confirmed PR) and then relapsed after prolonged therapy. These results suggest that the bortezomib insensitivity that develops in some initially responsive MM patients25 is not because of variants in PSMB5 or other catalytic proteasome subunits that form the bortezomib-binding pocket, and also indicate that alternative mechanisms underlie treatment-emergent bortezomib resistance in the clinical setting.
The online version of this article contains a data supplement.
The publication costs of this article were defrayed in part by page charge payment. Therefore, and solely to indicate this fact, this article is hereby marked “advertisement” in accordance with 18 USC section 1734.
Acknowledgments
The authors thank Emma Landers of FireKite for writing assistance, which was funded by Millennium Pharmaceuticals.
This work was supported by research funding from Millennium Pharmaceuticals and Janssen Global Services.
Authorship
Contribution: D.I.L., H.D., and G.M. designed the research; D.I.L., H.D., O.T., M.S., and G.M. performed the research; D.I.L., H.D., M.D.P, O.T., M.S., H.S., J.B., H.v.d.V., D.R., D.-L.E., W.L.T, G.M., and K.C.A. analyzed and interpreted the data; P.G.R., J.C., J.B., T.F., R.N., M.A., W.D., P.S., S.L., and K.C.A. contributed patients; D.I.L., H.D., O.T., M.S., H.S., P.G.R., J.C., J.B., T.F., R.N., M.A., W.D., P.S., S.L., G.M., and K.C.A. collected the data; M.D.P. and H.S. performed the statistical analysis; D.I.L., H.D., and G.M. wrote the initial draft of the manuscript; and all authors reviewed and critically revised the draft manuscript and approved the final manuscript.
Conflict-of-interest disclosure: D.I.L., H.D., M.D.P., O.T., M.S., H.S., D.-L.E., W.L.T, and G.M. are employees of Millennium Pharmaceuticals. D.-L.E. has ownership interests in Millennium/Takeda and Johnson & Johnson. P.G.R. has served on advisory boards for Millennium Pharmaceuticals, Celgene, Novartis, Bristol-Myers Squibb, and Janssen. J.C. and J.B. have received honoraria from Janssen and Celgene. J.B. has served on advisory boards for Janssen and Celgene and has received research funding from Celgene. T.F. has served on advisory boards/speakers' bureaus for Janssen. R.N. has served as a consultant and on advisory boards/speakers' bureaus for and received research funding from Millennium Pharmaceuticals, Celgene, and Onyx. M.A. has served on advisory boards/speakers' bureaus and as a consultant for Celgene, Millennium Pharmaceuticals, and Ortho Biotech Products and has received research support from Celgene and Millennium Pharmaceuticals. P.S. has served on advisory boards for Millennium Pharmaceuticals, Janssen, Celgene, Novartis, and Onyx and has received research support from Janssen and Celgene. S.L. has served as a consultant for Bristol-Myers Squibb, Celgene, Merck, Millennium Pharmaceuticals, Novartis, and Onyx. H.v.d.V. and D.R. are employees of and have ownership interest in Janssen Research & Development. K.C.A. has served on advisory boards for and received research funding from Bristol-Myers Squibb, Celgene, Merck, Millennium Pharmaceuticals, Novartis, and Onyx and has ownership interest in Acetylon. W.D. declares no competing financial interests.
Correspondence: David I. Lichter, Millennium Pharmaceuticals, 35 Landsdowne St, Cambridge, MA 02139; e-mail: david.lichter@MPI.com.
References
Author notes
D.I.L. and H.D. contributed equally to this work.
G.M. and K.C.A. share senior authorship.