• Direct analysis of the HLA-presented peptidome identifies a distinct antigenic signature in MM.

  • T-cell responses for these antigens are detectable exclusively in MM patients and can be induced in vitro in response-naive patients.

Direct analysis of HLA-presented antigens by mass spectrometry provides a comprehensive view on the antigenic landscape of different tissues/malignancies and enables the identification of novel, pathophysiologically relevant T-cell epitopes. Here, we present a systematic and comparative study of the HLA class I and II presented, nonmutant antigenome of multiple myeloma (MM). Quantification of HLA surface expression revealed elevated HLA molecule counts on malignant plasma cells compared with normal B cells, excluding relevant HLA downregulation in MM. Analyzing the presentation of established myeloma-associated T-cell antigens on the HLA ligandome level, we found a substantial proportion of antigens to be only infrequently presented on primary myelomas or to display suboptimal degrees of myeloma specificity. However, unsupervised analysis of our extensive HLA ligand data set delineated a panel of 58 highly specific myeloma-associated antigens (including multiple myeloma SET domain containing protein) which are characterized by frequent and exclusive presentation on myeloma samples. Functional characterization of these target antigens revealed peptide-specific, preexisting CD8+ T-cell responses exclusively in myeloma patients, which is indicative of pathophysiological relevance. Furthermore, in vitro priming experiments revealed that peptide-specific T-cell responses can be induced in response-naive myeloma patients. Together, our results serve to guide antigen selection for T-cell–based immunotherapy of MM.

Antigen-specific immunotherapy holds the potential to induce clinically effective anticancer T-cell responses1,2  and might be harnessed to guide and increase the specificity of cancer immunotherapy in future combination trials.3  To this end, the exact knowledge of tumor-associated/specific T-cell epitopes is crucial. After decades of research into overexpressed tumor antigens, more recently the focus has shifted to the patient-individualized identification of mutation-derived neoantigens.4,5  The encouraging findings of these new studies6-8  have led to neoepitopes being viewed as the dominant targets of anticancer immune responses.9-11 

However, analyzing the antigenome of hematologic malignancies, we have recently demonstrated that nonmutated antigens are relevant targets of spontaneous antileukemia T-cell responses.12,13  The strategy implemented in these studies differentially maps the naturally presented HLA ligandomes of hematologic cells in health and disease by mass spectrometry and was found to efficiently identify relevant tumor-associated antigens (TAAs).

Here, we translated this approach to multiple myeloma (MM), a low-grade B-cell lymphoma, characterized by the proliferation of malignant plasma cells in the bone marrow.14  Despite recent advances in treatment, including high-dose chemotherapy followed by autologous stem cell transplantation, novel immunomodulatory drugs, and proteasome inhibitors, MM remains largely incurable.15,16  This is mostly due to the persistence of minimal residual disease (MRD), which leads to high relapse rates.17,18  So far, the only established immunotherapeutic approach for MM is allogenic stem cell transplantation, which is associated with a high morbidity and mortality and remains an option for only a fraction of patients.19-21  Antigen-specific T-cell–based immunotherapy,22,23  especially in the constellation of MRD characterized by favorable effector-to-target ratios, might present an effective, low side-effect option.24 

An array of myeloma-associated T-cell antigens has been described in previous studies.25-35  Most of these antigens were identified based on gene expression analysis and reverse immunology. Some of these antigens (WT1,36,37  RHAMM,38,39  hTERT,40  and Survivin40,41 ) have already found their way into clinical trials, showing promising results in terms of induction of specific T-cell responses as well as clinical responses in single patients. However, broad clinical effectiveness has not yet been achieved. These previous studies were restricted to very few HLA allotypes and single antigens/epitopes,42  limiting both the population of patients eligible for this therapeutic approach and the spectrum of inducible tumor-specific T-cell responses. Of note, recent studies demonstrated lacking degrees of tumor association for several of these tumor antigens, both on the transcriptome level43  and importantly also on the level of HLA-restricted presentation.12,13 

By analyzing the antigenic landscape of MM directly on the HLA ligand level, we here provide new insights on antigenic distribution/specificity and identify a panel of novel myeloma-associated epitopes suited for antigen-specific immunotherapy.

Patients, blood, and bone marrow samples

Bone marrow mononuclear cells (BMNCs) and peripheral blood mononuclear cells (PBMCs) from MM patients at the time of diagnosis or at relapse before therapy, as well as PBMCs, BMNCs and granulocytes of healthy volunteers (HVs), were isolated by density gradient centrifugation (Biocoll; Biochrom GmbH) and erythrocyte lysis (EL buffer; Qiagen). Informed consent was obtained in accordance with the Declaration of Helsinki protocol. The study was performed according to the guidelines of the local ethics committee (142/2013BO2). Patient characteristics are provided in Table 1.44  HLA typing was carried out by the Department of Hematology and Oncology, University of Tübingen, Tübingen, Germany.

Table 1

Patient characteristics

UPNSexAge (y)ISS StageDurie and SalmonLight chain restrictionCytogenetic riskPrevious therapyHLA typingExperiment
62 3B λ High Yes A*26, A*30, B*15, B*42 Q; L 
55 1A Non Standard No A*01, A*24, B*18, B*08 
69 3B κ Standard No A*02, A*01, B*08, B*37 
72 3A κ Standard No A*02, A*33, B*15, B*18 
54 3A κ Low No A*03, A*26, B*40, B*55 
60 3B λ High No A*02, A*24, B*07, B*27 
39 3A κ Standard No A*02, A*03, B*07, B*35 
73 3B λ Unknown No A*24, A*25, B*39, B*40 
52 1B κ Unknown No A*02, B*07, B*44 
10 47 3A κ Standard No A*03, A*33, B*07 
11 72 3A κ High Yes A*02, B*27, B*44 
12 68 3A κ Standard No n.d. 
13 74 3A κ Standard No n.d. 
14 50 3A λ Standard No A*01, A*68, B*08, B*35 
15 62 3A κ Standard Yes n.d. 
16 74 2B κ Standard Yes n.d. 
17 68 3A κ Unknown Yes A*03, A*33, B*14, B*18 
18 73 1A κ Standard Yes A*02, A*32, B*40, B*44 
19 70 1A κ Standard No n.d. 
20 53 3A κ Standard No n.d. 
21 59 3A λ Standard Yes n.d. 
22 57 1A κ Standard Yes n.d. 
23 65 1A λ Standard Yes n.d. 
24 76 3A λ Standard Yes n.d. 
25 74 2A κ Standard Yes A*01, A*02, B*08, B*44 
26 67 3A λ Standard Yes n.d. 
27 50 1A λ Unkwnon Yes A*01, A*03, B*35, B*40 
28 70 3A λ Unknown No n.d. 
29 75 3A κ Standard Yes n.d. 
UPNSexAge (y)ISS StageDurie and SalmonLight chain restrictionCytogenetic riskPrevious therapyHLA typingExperiment
62 3B λ High Yes A*26, A*30, B*15, B*42 Q; L 
55 1A Non Standard No A*01, A*24, B*18, B*08 
69 3B κ Standard No A*02, A*01, B*08, B*37 
72 3A κ Standard No A*02, A*33, B*15, B*18 
54 3A κ Low No A*03, A*26, B*40, B*55 
60 3B λ High No A*02, A*24, B*07, B*27 
39 3A κ Standard No A*02, A*03, B*07, B*35 
73 3B λ Unknown No A*24, A*25, B*39, B*40 
52 1B κ Unknown No A*02, B*07, B*44 
10 47 3A κ Standard No A*03, A*33, B*07 
11 72 3A κ High Yes A*02, B*27, B*44 
12 68 3A κ Standard No n.d. 
13 74 3A κ Standard No n.d. 
14 50 3A λ Standard No A*01, A*68, B*08, B*35 
15 62 3A κ Standard Yes n.d. 
16 74 2B κ Standard Yes n.d. 
17 68 3A κ Unknown Yes A*03, A*33, B*14, B*18 
18 73 1A κ Standard Yes A*02, A*32, B*40, B*44 
19 70 1A κ Standard No n.d. 
20 53 3A κ Standard No n.d. 
21 59 3A λ Standard Yes n.d. 
22 57 1A κ Standard Yes n.d. 
23 65 1A λ Standard Yes n.d. 
24 76 3A λ Standard Yes n.d. 
25 74 2A κ Standard Yes A*01, A*02, B*08, B*44 
26 67 3A λ Standard Yes n.d. 
27 50 1A λ Unkwnon Yes A*01, A*03, B*35, B*40 
28 70 3A λ Unknown No n.d. 
29 75 3A κ Standard Yes n.d. 

Cytogenetic risk according to the IMWG risk stratification.44 

F, female; IMWG, International Myeloma Working Group; ISS, International Staging System; L, HLA ligand isolation; M, male; n.d., not determined; Q, HLA quantification; UPN, uniform patient number; y, years.

Myeloma cell lines

For HLA ligandome analysis, myeloma cell lines (MCLs; U266, RPMI 8226, JJN3, LP-1, MM.1S) were cultured in the recommended cell media (RPMI 1640, Gibco; Iscove modified Dulbecco medium [IMDM], Lonza) supplemented with 10%/20% fetal calf serum, 100 IU/L penicillin, 100 mg/L streptomycin, and 2 mmol/L glutamine at 37°C and 5% CO2. The MCLs RPMI 8226, JJN3, MM.1S and LP-1 were obtained from the Department of Hematology and Oncology, University of Tübingen, Tübingen, Germany.

Quantification of HLA surface expression

HLA surface expression on MM patient and HV bone marrow cells including CD38+CD138+ myeloma cells/plasma cells, CD19+CD20+ B cells, CD3+ T cells, and CD34+CD38 hematopoietic progenitor cells (HPCs) were analyzed using the QIFIKIT bead–based quantitative flow cytometric assay (Dako) according to the manufacturer’s instructions as described previously.12  In brief, samples were stained with the pan-HLA class I–specific monoclonal antibody (mAb) W6/32, HLA-DR–specific mAb L243 (produced in-house) or immunoglobulin G isotype control (BioLegend), respectively. Surface marker staining was carried out with directly labeled CD138, anti-κ, anti-λ, CD19, CD20 (BioLegend) and CD38, CD3, and CD34 (BD) antibodies. 7-aminoactinomycin D (BioLegend) was added as viability marker immediately before flow cytometric analysis on a LSR Fortessa (BD Biosciences).

Isolation of HLA ligands from primary samples and MCLs

HLA class I and II molecules were isolated using standard immunoaffinity purification as described45  using the pan-HLA class I–specific mAb W6/32, the pan-HLA class II–specific mAb Tü39, and the HLA-DR–specific mAb L243 (produced in-house).

Analysis of HLA ligands by liquid chromatography–tandem mass spectrometry

HLA ligand extracts were analyzed in 5 technical replicates as described previously.13  In brief, peptide samples were separated by nanoflow high-performance liquid chromatography (RSLCnano; Thermo Fisher Scientific) using a 50 μm × 25 cm PepMap rapid separation liquid chromatography column (Thermo Fisher Scientific) and a gradient ranging from 2.4% to 32.0% acetonitrile over the course of 90 minutes. Eluting peptides were analyzed in an online-coupled LTQ Orbitrap XL mass spectrometer (Thermo Fisher Scientific) using a top 5 collision-induced dissociation fragmentation method.

Database search and spectral annotation

Data processing was performed as described previously.13  In brief, the Mascot search engine (Mascot 2.2.04; Matrix Science) was implemented to search the human proteome as comprised in the Swiss-Prot database (20 279 reviewed protein sequences, September 2013) without enzymatic restriction. Potential mutated HLA ligands were searched implementing a database containing the human proteome concatenated with proteins containing single amino acid variants listed in the COSMIC database (http://cancer.sanger.ac.uk/cosmic/). Only recurrent single amino acid variants described in 2 or more samples of hematologic origin were included. Oxidized methionine was allowed as a dynamic modification. The false discovery rate (FDR) was estimated using the Percolator algorithm46  and set to 5%. Peptide lengths were limited to 8 to 12 aa for HLA class I and 12 to 25 aa for HLA class II. Protein inference was disabled, allowing for multiple protein annotations of peptides. HLA annotation was performed using SYFPEITHI47  or an extended in-house database. Experimental validation of peptide identifications and HLA annotations was performed by mass spectrometric and functional characterization of synthetic peptides for a subset of peptides.

Peptide and HLA peptide monomer synthesis

The automated peptide synthesizer EPS221 (Abimed) was used to synthesize peptides using the 9-fluorenylmethyl-oxycarbonyl/tert-butyl (Fmoc/tBu) strategy.48  Synthetic peptides were used for validation of liquid chromatography–tandem mass spectrometry (MS) identifications as well as for functional experiments. Biotinylated recombinant HLA molecules and fluorescent HLA-peptide tetramers were produced as described previously.49 

Amplification of peptide-specific T cells

PBMCs from MM patients and HVs were cultured as described previously.12,13  In brief, for CD8+ T-cell stimulation, PBMCs were pulsed with 1 µg/mL per peptide and cultured for 12 days adding interleukins 4 and 7 (IL-4 and IL-7) on days 0 and 1 as well as IL-2 on days 3, 5, 7, and 9. HLA-A*02 (KLFEKVKEV)- and B*07 (KPSEKIQVL)-restricted control peptides derived from benign tissues (HV-exclusive HLA ligands) served as negative control. Peptide-stimulated PBMCs were analyzed by enzyme-linked immunospot (ELISPOT) assays on day 12. For CD4+ T-cell stimulation, culture was performed as described for CD8+ T cells except for 2 modifications: pulsing was carried out with 10 µg/mL HLA class II peptide and no IL-4 or IL-7 was added.

IFNγ ELISPOT assay

Interferon γ (IFNγ) ELISPOT assays were carried out as described previously.50  In brief, 96-well nitrocellulose plates (Millipore) were coated with 1 mg/mL IFNγ mAb (Mabtech) and incubated overnight at 4°C. Plates were blocked with 10% human serum for 2 hours at 37°C. Prestimulated PBMCs (2.5 × 105 cells per well) were pulsed with 1 µg/mL (HLA class I) or 2.5 µg/mL (HLA class II) peptide and incubated for 24 to 26 hours. Readout was performed according to the manufacturer’s instructions. Phytohemagglutinin was used as positive control. HLA-A*02 (KLFEKVKEV)- and B*07 (KPSEKIQVL)-restricted control peptides derived from benign tissues (HV-exclusive HLA ligands) served as negative control. Spots were counted using an ImmunoSpot S5 analyzer (CTL). T-cell responses were considered to be positive when >10 spots per well were counted and the mean spot count per well was at least threefold higher than the mean number of spots in the negative control wells (according to the cancer immunoguiding program guidelines51 ).

aAPC priming of peptide-specific T cells

For the generation of artificial antigen-presenting cells (aAPCs), 4 × 106 streptavidin-coated polystyrene particles (Bangs Laboratories) per milliliter were resuspended in PBE (phosphate-buffered saline/bovine serum albumin/EDTA; Gibco/Sigma Aldrich/Lonza) containing 200 pM biotinylated MHC-peptide monomer and 20 nM antihuman biotinylated CD28 antibody and incubated at room temperature for 30 minutes. After washing, the aAPCs were stored at 4°C prior to use.52  CD8+ T cells from MM patients and HV were enriched by positive selection using magnetic cell sorting (Miltenyi Biotec). Stimulations were initiated in 96-well plates with 1 × 106 T cells plus 2 × 105 aAPCs in 200 μL of T-cell medium complemented with 5 ng/mL human IL-12 (PromoKine). IL-2 (65 U/μL; R&D Systems) was added on day 5. aAPC stimulation was repeated on day 10, for a total of 3 cycles.

Tetramer staining

The frequency of peptide-specific CD8+ T cells was determined on a FACSCanto II cytometer (BD Bioscience) by staining with anti-CD8 (Biolegend) and HLA:peptide-tetramer-phycoerythrin as described previously.52  Staining with tetramers containing the cytomegalovirus (CMV) pp65 A*02 peptide NLVPMVATV served as positive control, tetramers containing irrelevant; nonprimed A*02-restricted control peptides served as negative controls. Successful priming was considered if frequency of peptide-specific CD8+ T cells was >0.1% of viable cells and at least threefold higher than the frequency of peptide-specific CD8+ T cells in the negative control.

Software and statistical analysis

Flow cytometric data analysis was performed using FlowJo 7.2 (TreeStar). In-house R and Python scripts were used for the generation of virtual ligandomes and definition of virtual TAAs in the analysis of TAA FDRs and for the TAA-plateau regression analysis. The standard R heatmap.2 script was used for the unsupervised cluster analysis of HLA ligand source proteins. GraphPad Prism 6.0 (GraphPad Software) was used for statistical analysis. Statistical analysis of HLA surface expression was based on unpaired t tests.

HLA class I surface expression is elevated on myeloma cells

As loss or downregulation of HLA expression on target cells might severely hamper the effectiveness of T-cell–based immunotherapy, we quantified HLA class I and II surface molecule counts on primary myeloma cells compared with autologous hematopoietic cells and plasma cells derived from the bone marrow of HVs. In MM patients (n = 20), HLA class I expression was found to be heterogeneous with mean expression levels on CD38+CD138+ myeloma cells of 416 000 ± 54 500, which was found to be significantly higher as compared with autologous normal CD19+CD20+ B cells (198 5000 ± 20 500, P = .001), CD3+ T cells (167 500 ± 15 500, P = .0002), and CD34+CD38 HPCs (204 000 ± 32 500, P = .002, Figure 1A). In addition, HLA class I expression on primary MM cells was also found to be significantly higher than that on CD38+CD138+ plasma cells of HVs (n = 15, 291 500 ± 25 500, P < .05; Figure 1C). No significant differences in HLA class I expression were observed when comparing normal B cells, T cells, and HPCs of MM patients to the corresponding cell populations of HVs (supplemental Figure 1, see supplemental Data available on the Blood Web site). HLA-DR expression levels on myeloma cells were generally found to be much lower than HLA class I levels. Mean HLA-DR surface molecule counts on myeloma cells (27 000 ± 7000) showed no significant difference compared with autologous HPCs (35 000 ± 5000) and T cells (18 000 ± 13 000) or plasma cells of HVs (39 500 ± 5000) (Figure 1B,D). HLA-DR expression of MM patient CD19+CD20+ B cells (104 000 ± 7000) was significantly higher compared with myeloma cells (P < .0001). No correlation of HLA surface expression on myeloma cells with patient characteristics such as sex, age, disease stage, risk classification, or prior therapy was observed (supplemental Table 1).

Figure 1

HLA class I and II surface expression on myeloma patient and HV bone marrow cells. Quantification of HLA surface expression was performed using a bead-based flow cytometric assay. (A) HLA class I and (B) HLA-DR expression on CD38+CD138+ primary myeloma cells compared with autologous CD34+CD38 HPCs, CD19+CD20+ B cells, and CD3+ T cells. (C) HLA class I and (D) HLA-DR expression on primary MM cells compared with bone marrow–derived plasma cells of HVs. n.s., not significant; *P < .05; **P < .01; ***P < .001.

Figure 1

HLA class I and II surface expression on myeloma patient and HV bone marrow cells. Quantification of HLA surface expression was performed using a bead-based flow cytometric assay. (A) HLA class I and (B) HLA-DR expression on CD38+CD138+ primary myeloma cells compared with autologous CD34+CD38 HPCs, CD19+CD20+ B cells, and CD3+ T cells. (C) HLA class I and (D) HLA-DR expression on primary MM cells compared with bone marrow–derived plasma cells of HVs. n.s., not significant; *P < .05; **P < .01; ***P < .001.

Close modal

MS identifies naturally presented HLA ligands of primary MM samples and MCLs

Mapping the HLA class I ligandomes of 10 myeloma patients and 5 MCLs (Table 1), we identified a total of 17 583 different peptides representing 7574 source proteins, attaining >80% of the maximum attainable coverage (Figure 2A). The mean number of unique peptide identifications (IDs) was 1059 IDs for primary myeloma samples and 2243 IDs for MCLs (supplemental Table 2). Overall, peptides restricted by 20 different HLA-A and -B allotypes were identified in this study, covering 99.3% of the German population (calculated according to Bui et al53 ). No mutated peptide identifications were obtained when processing the data set against a reference proteome containing recurrent mutations in hematologic malignancies.

Figure 2

Comparative HLA ligandome profiling and identification of myeloma-associated antigens. (A) Saturation analysis of HLA class I ligand source protein identifications in MM patients. Number of unique HLA ligand source protein identifications as a function of cumulative HLA ligand source protein identifications in 10 MM patients. Exponential regression allowed for the robust calculation (R2 = 0.99) of the maximum attainable number of different source protein identifications (dashed line). The dotted line depicts the source proteome coverage achieved in our MM patient cohort. (B) Overlap analysis of HLA class I ligand source proteins of primary MM samples (n = 10), MCLs (n = 5), and HV samples (total n = 45: PBMC [n = 30], BMNC [n = 10], granulocytes [n = 5]). (C) Comparative profiling of HLA ligand source proteins based on the frequency of HLA-restricted presentation in MM and HV ligandomes. Frequencies of MMs/HVs positive for HLA-restricted presentation of the respective source protein (x-axis) are indicated on the y-axis. The box on the left highlights the subset of myeloma-associated antigens showing MM-exclusive presentation in >25% of myeloma samples. (D) Statistical assessment of false-positive myeloma-antigen IDs at different threshold values. The numbers of original TAAs identified based on the analysis of the MM and HV cohorts were compared with random virtual TAAs. Virtual MM and HV samples were generated in silico based on random weighted sampling from the entirety of protein identifications in both original cohorts. These randomized virtual ligandomes of defined size (n = 957 proteins, which is the mean number of protein identifications in all analyzed samples) were used to define TAAs based on simulated cohorts of 15 MM vs 45 HV samples. The process of protein randomization, cohort assembly, and TAA identification was repeated 1000 times and the mean value of resultant virtual TAAs was calculated and plotted for the different threshold values. The corresponding FDRs for any chosen TAA threshold are listed below the x-axis. sum, summary.

Figure 2

Comparative HLA ligandome profiling and identification of myeloma-associated antigens. (A) Saturation analysis of HLA class I ligand source protein identifications in MM patients. Number of unique HLA ligand source protein identifications as a function of cumulative HLA ligand source protein identifications in 10 MM patients. Exponential regression allowed for the robust calculation (R2 = 0.99) of the maximum attainable number of different source protein identifications (dashed line). The dotted line depicts the source proteome coverage achieved in our MM patient cohort. (B) Overlap analysis of HLA class I ligand source proteins of primary MM samples (n = 10), MCLs (n = 5), and HV samples (total n = 45: PBMC [n = 30], BMNC [n = 10], granulocytes [n = 5]). (C) Comparative profiling of HLA ligand source proteins based on the frequency of HLA-restricted presentation in MM and HV ligandomes. Frequencies of MMs/HVs positive for HLA-restricted presentation of the respective source protein (x-axis) are indicated on the y-axis. The box on the left highlights the subset of myeloma-associated antigens showing MM-exclusive presentation in >25% of myeloma samples. (D) Statistical assessment of false-positive myeloma-antigen IDs at different threshold values. The numbers of original TAAs identified based on the analysis of the MM and HV cohorts were compared with random virtual TAAs. Virtual MM and HV samples were generated in silico based on random weighted sampling from the entirety of protein identifications in both original cohorts. These randomized virtual ligandomes of defined size (n = 957 proteins, which is the mean number of protein identifications in all analyzed samples) were used to define TAAs based on simulated cohorts of 15 MM vs 45 HV samples. The process of protein randomization, cohort assembly, and TAA identification was repeated 1000 times and the mean value of resultant virtual TAAs was calculated and plotted for the different threshold values. The corresponding FDRs for any chosen TAA threshold are listed below the x-axis. sum, summary.

Close modal

As controls, we analyzed the HLA class I ligandomes of 45 HV-derived samples (30 PBMC, 10 BMNC, and 5 granulocyte specimens) identifying a total of 20 171 different peptides representing 7729 source proteins (supplemental Table 3). The HLA allotype distribution in the HV cohort covered >80% of HLA-A and -B alleles in the MM sample cohort.54  Analysis of HLA class II ligandomes was performed for 7 MM patients and 5 MCLs. A total of 6076 unique peptides representing 1743 source proteins were identified. The HLA class II HV cohort (13 PBMC, 5 BMNC, 5 granulocyte specimens) yielded 2899 different peptides representing 889 source proteins (supplemental Table 2).

HLA class I ligandome profiling identifies a novel panel of myeloma-associated antigens

To identify myeloma-associated antigens, we comparatively analyzed the HLA ligandomes of the MM sample and HV cohorts at the source protein level. Overlap analysis of HLA ligand source proteins identified 2412 proteins (corresponding to 31.8% of the mapped MM HLA source proteome) to be exclusively represented in the HLA ligandomes of MM samples. Of these MM-exclusive source proteins, 68.3% were solely identified on MCL samples, whereas 13.2% of proteins were found to be presented both on MCLs and primary MM samples. A fraction of 18.5% of myeloma-exclusive source proteins was found to be restricted to primary MM samples (Figure 2B). To identify broadly presented TAAs, myeloma-exclusive source proteins were ranked according to their frequencies of representation in the MM sample cohort (Figure 2C). To statistically assess and optimize the stringency of antigen identification, we simulated randomized virtual ligandomes in silico and calculated the resultant number of false-positive TAAs at different frequencies of representation (Figure 2D). We set the frequency threshold for HLA class I TAA definition to >25% of myeloma-exclusive antigen presentation, yielding 58 TAAs with an estimated FDR of 4.1% (supplemental Table 4). This novel panel of frequently presented myeloma-associated antigens was represented by 197 unique HLA class I ligands and constitutes 0.8% of the mapped myeloma HLA ligand source proteome. Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway analysis55  and functional annotation clustering of these antigens with respect to their biological function (GO Term BP FAT56 ) did not identify any statistically significant overrepresented pathways or functional clusters. Notably, the proto-oncogene multiple myeloma SET domain containing protein (MMSET) was identified as a TAA showing representation in 33% of MM patient ligandomes and was found to be represented by 3 different HLA ligands (ASNPSNPRPSK (HLA-A*30:01), KAMEAASSL (A*02:01), SLLEQGLVEA (A*02:01)). Moreover, MMSET was detected on both MM patients with the oncogenic translocation t(4;14), but only on 1 of 6 patients (17%) without this aberration.

Representation of established myeloma-associated antigens in the HLA class I ligandome

Based on our extensive HLA ligandome data set, we investigated the presentation of established myeloma-associated antigens within the different sample cohorts. We identified 73 different HLA ligands representing 22 of 25 (88%) previously described myeloma antigens.42  We found 9 of the 22 detectable antigens (41%) to be exclusively presented on MM samples, 10 of 22 antigens (45.5%) to be represented both on MM and HV samples, and 3 of 22 (13.6%) exclusively presented on HV-derived samples (Figure 3A). Of note, 7 of 9 MM-exclusive antigens (77.8%) were only detectable on MCLs. Only 2 of 9 (22.2%) of these MM-exclusive antigens HLA ligands were detected on primary MM patient samples (Figure 3B). For reference, only 7 of 58 (12.1%) of the newly defined myeloma antigens showed presentation exclusively on MCLs, whereas the majority (51 of 58 [87.9%]) of antigens were also presented on primary MM patient samples as well, which underlines their potential as clinical target antigens (Figure 3C).

Figure 3

Representation of established myeloma-associated antigens in the HLA ligandomes of MM and HV. (A) Representation of previously described MM-associated antigens in HLA class I ligandomes. Bars indicate relative representation (%) of respective antigens by HLA class I ligands on primary MM samples, MCLs, and HV samples. Dashed lines divide the antigens into 4 groups according to their degree of MM association (MM and MCL-exclusive, MCL-exclusive, mixed presentation, HV-exclusive). (B-C) Distribution of myeloma-exclusive antigen presentation for (B) previously described antigens and (C) ligandome-defined TAAs on MCLs (white) and MM plus MCLs (shaded).

Figure 3

Representation of established myeloma-associated antigens in the HLA ligandomes of MM and HV. (A) Representation of previously described MM-associated antigens in HLA class I ligandomes. Bars indicate relative representation (%) of respective antigens by HLA class I ligands on primary MM samples, MCLs, and HV samples. Dashed lines divide the antigens into 4 groups according to their degree of MM association (MM and MCL-exclusive, MCL-exclusive, mixed presentation, HV-exclusive). (B-C) Distribution of myeloma-exclusive antigen presentation for (B) previously described antigens and (C) ligandome-defined TAAs on MCLs (white) and MM plus MCLs (shaded).

Close modal

Moreover, unsupervised clustering of source protein presentation in the HLA ligandomes revealed the cluster of MCLs to be highly distinct from primary MM samples (supplemental Figure 2).

Analysis of HLA class II ligandomes identifies potentially synergistic vaccine candidates

As the direct involvement of CD4+ T cells in tumor control is established,57  we further aimed to identify HLA class II antigens. Overlap analysis of HLA class II ligand source proteins identified 1135 myeloma exclusive antigens (Figure 4A). Comparative profiling of HLA class II ligandomes identified a single antigen (TFRC) represented by 67 HLA class II ligands showing MM-exclusive presentation at FDR <5% (Figure 4B-C). Functional characterization of the most abundant TFRC peptide (NSVIIVDKNGRLV) by IFNγ ELISPOT revealed memory T-cell responses in 2 of 5 MM patients (supplemental Figure 3A).

Figure 4

Identification of synergistic HLA class II–restricted myeloma-associated antigens. (A) Overlap analysis of HLA class II ligand source proteins of primary MM samples (n = 7), MCLs (n = 5), and HV samples (total, n = 23: PBMC [n = 13], BMNC [n = 5], granulocytes [n = 5]). (B) Statistical analysis of false-positive myeloma-antigen identifications at different threshold values, as described in Figure 2. Randomized virtual ligandome sizes were set to 226 proteins and TAAs were defined based on simulated cohorts of 12 MM vs 23 HV samples. (C) Comparative profiling of HLA class II ligand source proteins based on the frequency of HLA-restricted presentation in MM and HV ligandomes. Frequencies of MMs/HVs positive for HLA-restricted presentation of the respective source protein (x-axis) are indicated on the y-axis. (D) Overlap analysis of HLA class I TAAs (n = 58) and HLA class II MM-exclusive antigens (n = 1135). (E) HLA class I TAAs, which also yield potentially synergistic HLA class II ligands. (F) Overlap analysis comprising the entire HLA class I and II ligand source proteomes of MM samples. sum, summary; rep., representation.

Figure 4

Identification of synergistic HLA class II–restricted myeloma-associated antigens. (A) Overlap analysis of HLA class II ligand source proteins of primary MM samples (n = 7), MCLs (n = 5), and HV samples (total, n = 23: PBMC [n = 13], BMNC [n = 5], granulocytes [n = 5]). (B) Statistical analysis of false-positive myeloma-antigen identifications at different threshold values, as described in Figure 2. Randomized virtual ligandome sizes were set to 226 proteins and TAAs were defined based on simulated cohorts of 12 MM vs 23 HV samples. (C) Comparative profiling of HLA class II ligand source proteins based on the frequency of HLA-restricted presentation in MM and HV ligandomes. Frequencies of MMs/HVs positive for HLA-restricted presentation of the respective source protein (x-axis) are indicated on the y-axis. (D) Overlap analysis of HLA class I TAAs (n = 58) and HLA class II MM-exclusive antigens (n = 1135). (E) HLA class I TAAs, which also yield potentially synergistic HLA class II ligands. (F) Overlap analysis comprising the entire HLA class I and II ligand source proteomes of MM samples. sum, summary; rep., representation.

Close modal

As CD4+ T cells play pivotal roles in the induction and maintenance of antigen-specific CD8+ T-cell responses,58-60  we implemented a second approach to identify potentially synergistic HLA class II–restricted peptides derived from HLA class I TAAs. Overlap analysis of the 58 HLA class I antigens with the 1135 HLA class II presented MM-exclusive proteins identified a panel of 6 class-spanning antigens represented by 31 peptides (Figure 4D-E; supplemental Table 5). Functional characterization of synergistic HLA class II ligands revealed peptide-specific T-cell responses in myeloma patients for 3 of 5 tested peptides (Figure 4E, supplemental Figure 3B).

The overall comparison of the HLA class I and II ligandomes of MM samples revealed 80% (1622) of HLA class II–presented proteins to be also presented on HLA class I (Figure 4F). Functional annotation clustering (GO Term CC clustering using DAVID56 ) was performed on the top 500 most frequently presented proteins in each HLA class to identify the cellular compartments from which these proteins derive. Antigens presented on class I displayed highly enriched clusters for nuclear proteins as well as for ribosomal, cytoskeletal, and vesicle-derived proteins. Notably, this pattern was recapitulated in the clustering of proteins presented on both HLA classes, albeit with a higher ranking and an almost threefold higher enrichment for vesicle-derived proteins. HLA class II–presented antigens showed intermediate enrichment for plasma membrane, vesicle-derived, and lysosomal proteins (supplemental Table 6).

HLA class I TAAs are targeted by spontaneous T-cell responses in myeloma patients

Functional characterization of the novel myeloma antigens was performed in panels of 11 HLA-A*02– and 2 HLA-B*07–restricted peptides, including 2 HLA-A*02 ligands derived from MMSET (Figure 5A). Myeloma-associated peptides were evaluated in 12-day recall IFNγ ELISPOT assays using PBMC obtained from MM patients and HVs. We observed IFNγ secretion for 5 of 11 A*02 ligands and 1 of 2 B*07 ligands in myeloma patients, as shown exemplarily in Figure 5C. Both peptides (P1 and P2) derived from MMSET showed specific T-cell recognition in 2 of 16 (13%) and 1 of 8 (13%) MM patients, respectively. Importantly, no myeloma peptide-specific IFNγ secretion was observed in 10 HLA-matched healthy controls (Figure 5B). Notably, T-cell responses were only observed for myeloma-associated peptides identified on primary myeloma samples (10 of 13), and never for the 3 of 13 peptides identified on MCLs only. The frequencies of peptide-specific T-cell responses detected in MM patients by ELISPOT were generally in the same range as the frequencies of presentation of the respective peptide in allotype-matched ligandomes of MM patients (Figure 5A). Due to limitations in the numbers of cells available for analysis, further controls with target cells expressing the corresponding antigens could not be performed. We therefore cannot exclude that T-cell reactivity is directed against impurities contained in the synthetic peptide batch. Indeed, it is well known that synthetic peptides contain impurities, for example, peptides modified with a protecting group, and that these impurities are immunogenic. However, HLA-A*02– and -B*07–restricted control peptides derived from benign tissues (HV-exclusive HLA ligands) used in all ELISPOTs in the study at hand never resulted in significant IFNγ release (Figure 5C).

Figure 5

Functional characterization of myeloma-associated antigens. (A) Myeloma-associated T-cell epitopes with their corresponding HLA restrictions and frequencies of immune recognition by myeloma patient-derived T cells in IFNγ-ELISPOT assays. (B) Example of myeloma-associated T-cell epitopes evaluated in an IFNγ-ELISPOT using HV PBMCs. An EBV epitope mix containing the frequently recognized peptides BRLF109-117 YVLDHLIVV (A*02) and EBNA3247-255 RPPIFIRRL (B*07 served as positive control). Benign-tissue–derived peptides KLFEKVKEV (HLA-A*02) and KPSEKIQVL (B*07) served as negative control. (C) Examples of myeloma-associated T-cell epitopes evaluated in IFNγ-ELISPOTs using MM patient PBMCs (n = 3). Results are shown only for immunoreactive peptides. An EBV epitope mix containing 5 frequently recognized peptides (BRLF109-117 YVLDHLIVV [A*02], EBNA3471-479 RLRAEAQVK [A*03], EBNA3247-255 RPPIFIRRL [B*07], BZLF1190-197 RAKFKQLL [B*08], EBNA6162-171 AEGGVGWRHW [B*44]) was used as positive control. Benign-tissue–derived peptides KLFEKVKEV (HLA-A*02) and KPSEKIQVL (B*07) served as negative control. (D-E) Tetramer staining of CD8+ T cells after 3 cycles of aAPC-based in vitro primings using T cells derived from (D) a healthy individual and (E) a myeloma patient: Leftmost panels, P2-tetramer staining of CD8+ T cells primed with P2-aAPCs (SLLEQGLVEA, A*02); left middle panels, ex vivo P2-tetramer staining of CD8+ T cells; right middle panels, control staining with A*02-tetramer containing a nonrelevant A*02-restricted control peptide (KAMEAASSL, A*02) on CD8+ T cells derived from the same population as T cells depicted in the left panels. Rightmost panels, Positive control: tetramer staining of CD8+ T cells primed with CMV-aAPCs (NLVPMVATV, A*02). EBV, Epstein-Barr virus; neg., negative; pos., positive; UPN, uniform patient number.

Figure 5

Functional characterization of myeloma-associated antigens. (A) Myeloma-associated T-cell epitopes with their corresponding HLA restrictions and frequencies of immune recognition by myeloma patient-derived T cells in IFNγ-ELISPOT assays. (B) Example of myeloma-associated T-cell epitopes evaluated in an IFNγ-ELISPOT using HV PBMCs. An EBV epitope mix containing the frequently recognized peptides BRLF109-117 YVLDHLIVV (A*02) and EBNA3247-255 RPPIFIRRL (B*07 served as positive control). Benign-tissue–derived peptides KLFEKVKEV (HLA-A*02) and KPSEKIQVL (B*07) served as negative control. (C) Examples of myeloma-associated T-cell epitopes evaluated in IFNγ-ELISPOTs using MM patient PBMCs (n = 3). Results are shown only for immunoreactive peptides. An EBV epitope mix containing 5 frequently recognized peptides (BRLF109-117 YVLDHLIVV [A*02], EBNA3471-479 RLRAEAQVK [A*03], EBNA3247-255 RPPIFIRRL [B*07], BZLF1190-197 RAKFKQLL [B*08], EBNA6162-171 AEGGVGWRHW [B*44]) was used as positive control. Benign-tissue–derived peptides KLFEKVKEV (HLA-A*02) and KPSEKIQVL (B*07) served as negative control. (D-E) Tetramer staining of CD8+ T cells after 3 cycles of aAPC-based in vitro primings using T cells derived from (D) a healthy individual and (E) a myeloma patient: Leftmost panels, P2-tetramer staining of CD8+ T cells primed with P2-aAPCs (SLLEQGLVEA, A*02); left middle panels, ex vivo P2-tetramer staining of CD8+ T cells; right middle panels, control staining with A*02-tetramer containing a nonrelevant A*02-restricted control peptide (KAMEAASSL, A*02) on CD8+ T cells derived from the same population as T cells depicted in the left panels. Rightmost panels, Positive control: tetramer staining of CD8+ T cells primed with CMV-aAPCs (NLVPMVATV, A*02). EBV, Epstein-Barr virus; neg., negative; pos., positive; UPN, uniform patient number.

Close modal

Antigen-specific T cells can be induced in vitro from naive T cells of MM patients or HVs

To assess whether myeloma antigen-specific T-cell responses can be induced from naive T cells in vitro, we isolated CD8+ T cells from 1 healthy individual and 1 MM patient. We performed aAPC primings using the MMSET-derived peptide SLLEQGLVEA (P2). Using HV-derived CD8+ T cells, a population of 0.403% P2-tetramer positive CD8+ T cells was detected after in vitro priming. No tetramer-positive T-cell populations >0.1% were detectable ex vivo. After priming of T cells from an MM patient without previous T-cell reactivity for P2 (as detected by 12-day-recall IFNγ-ELISPOT and ex vivo tetramer staining), we detected the induction of a small population of 0.236% P2-tetramer positive CD8+ T cells (Figure 5E). Importantly, control stainings performed with an A*02-tetramer containing a nonrelevant A*02 control peptide were performed in parallel on cells derived from the same wells as used for the relevant staining and did not yield any specific tetramer-positive T-cell populations (Figure 5D).

With the advent of immune checkpoint modulation, T-cell–based immunotherapy has opened up new avenues for the treatment of a range of solid tumors61-68  and is undergoing clinical evaluation in hematologic malignancies.69,70  These novel therapeutic options may further be improved in combination with antigen-specific immunotherapy, which may help induce and guide specific anticancer immune responses. To this end, the exact knowledge of tumor-associated epitopes which can act as rejection antigens is indispensable.4,8  In myeloma, a multitude of studies have investigated MM-associated antigens, yielding a handful of promising targets.25-35  However, these studies were limited in scope of antigens and HLA allotypes.42  In our previous studies in acute myeloid leukemia and chronic lymphocytic leukemia, we demonstrated that the comprehensive and comparative analysis of the HLA-presented antigenome can identify extensive panels of broadly presented T-cell epitopes covering a multitude of HLA allotypes.12,13  Here, we translated this systematic approach to MM.

Quantification of HLA surface expression on different cell populations in the bone marrow of myeloma patients and healthy volunteers demonstrated that HLA loss or downregulation on malignant plasma cells is of no concern, even in patients who received prior therapy. Comparative analysis of the HLA ligandomes of these cell populations revealed distinct antigenic signatures and identified a panel of myeloma-associated antigens. These antigens did not cluster into any of the major tumor-associated pathways and did not show enrichment for highly prioritized targets characterized by tumor-associated overexpression. This indicates a more complex interplay of underlying mechanisms such as protein turnover71  and antigen processing as the source of altered antigen presentation on tumor cells and further underscores the isolated character of the HLA ligandome and the distorted correlation with its upstream sources.72,73 

Importantly, a substantial proportion of established myeloma antigens was found to be only infrequently presented on primary myelomas or to show suboptimal degrees of myeloma specificity. Of note, the majority of these antigens was selectively detected on myeloma cell lines but not in primary samples, indicating that selection of pathophysiologically relevant antigens should be based on analysis of primary tumor samples.

A notable exception was the established myeloma-associated protein MMSET, which is currently being investigated as a target for the therapy of poor-prognosis t(4;14) myeloma patients.74-77  Although MMSET-derived peptides were frequently identified on t(4;14) myeloma samples, we also detected MMSET peptides in the HLA ligandomes of a t(4;14)-negative patient and 1 t(4;14)-negative MCL (U266). Strikingly, functional characterization by ELISPOT revealed memory T-cell responses targeting these MMSET-derived epitopes exclusively in myeloma patients and not in HV. This suggests myeloma-dependent priming of anti-MMSET T-cell responses in vivo in MM patients, which underscores the pathophysiological relevance of this antigen. In concordance with the HLA ligandomics data, we found these T-cell responses not to be restricted to t(4;14) myeloma patients. Results of in vitro primings suggest that MMSET-specific CD8+ T-cell responses can be induced from naive T cells, both in healthy individuals and, importantly, also in myeloma patients, albeit with limited magnitudes. With the current strategies focusing on inhibition of MMSET by small molecules or small interfering RNAs,78,79  our identification of myeloma-exclusive MMSET-derived T-cell epitopes provides new options for targeting MMSET by T-cell–based immunotherapy. Notably, this therapeutic strategy may not necessarily have to be restricted to t(4;14) myelomas, as we observed MMSET-presentation and immune recognition irrespective of the mutational status. This might be explained by the distorted correlation of gene expression and HLA-restricted antigen presentation as well as by the subclonal distribution of t(4;14) in myeloma cells and genomic plasticity occurring over the course of disease.73,80 

Together, our findings illustrate how antigen identification guided by HLA ligandomics can pinpoint novel MM-associated T-cell epitopes and allow for direct assessment of antigen distribution patterns in patient cohorts. In parallel to our findings with MMSET, our study features an extensive panel of novel antigens previously not associated with myeloma or cancer in general. Analogously to MMSET, we detected preexisting T-cell responses against a substantial proportion of these targets in myeloma patients, indicating a high enrichment for relevant MM-associated antigens. In conclusion, our ligandome-centric study may guide the design of future antigen-specific T-cell immunotherapy in MM.

The online version of this article contains a data supplement.

There is an Inside Blood Commentary on this article in this issue.

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.

The authors thank Claudia Falkenburger, Patricia Hrstic, Nicole Zuschke, Katharina Graf, and Beate Pömmerl for excellent technical support.

This work was supported by the German Cancer Consortium (DKTK), the Deutsche Forschungsgemeinschaft (DFG; SFB 685), and the European Union (EU; ERC AdG 339842 MUTAEDITING).

Contribution: S.W., J.S.S., D.J.K., S.S., and H.-G.R. designed the study; S.W., J.S.S., D.J.K., H.-G.R., L.K., K.W., H.R.S., and S.S. drafted the manuscript; S.W., J.S.S., and H.S. performed flow cytometric analysis and HLA surface expression quantitation; S.W., D.J.K., J.S.S., and S.K. conducted MHC peptide immunoaffinity purification; S.W. and D.J.K. performed MS-based analyses and quantitation; S.W., A.N., T.S., J.S.S., and H.S. conducted in vitro T-cell experiments; J.S.S., L.K., K.W., M.H., and H.R.S. conducted patient data collection and medical evaluation; S.W., J.S.S., D.J.K., O.K., and L.B. performed statistical analyses; and S.S., H.-G.R., and J.S.S. supervised the study.

Conflict-of-interest disclosure: The authors declare no competing financial interests.

Correspondence: Juliane S. Stickel, Department of Hematology and Oncology, University of Tübingen, Otfried-Müller Strasse 10, 72076 Tübingen, Germany; e-mail: juliane.stickel@med.uni-tuebingen.de.

1
Small
 
EJ
Schellhammer
 
PF
Higano
 
CS
et al. 
Placebo-controlled phase III trial of immunologic therapy with sipuleucel-T (APC8015) in patients with metastatic, asymptomatic hormone refractory prostate cancer.
J Clin Oncol
2006
, vol. 
24
 
19
(pg. 
3089
-
3094
)
2
Walter
 
S
Weinschenk
 
T
Stenzl
 
A
et al. 
Multipeptide immune response to cancer vaccine IMA901 after single-dose cyclophosphamide associates with longer patient survival.
Nat Med
2012
, vol. 
18
 
8
(pg. 
1254
-
1261
)
3
Perez-Gracia
 
JL
Labiano
 
S
Rodriguez-Ruiz
 
ME
Sanmamed
 
MF
Melero
 
I
Orchestrating immune check-point blockade for cancer immunotherapy in combinations.
Curr Opin Immunol
2014
, vol. 
27
 (pg. 
89
-
97
)
4
van Rooij
 
N
van Buuren
 
MM
Philips
 
D
et al. 
Tumor exome analysis reveals neoantigen-specific T-cell reactivity in an ipilimumab-responsive melanoma.
J Clin Oncol
2013
, vol. 
31
 
32
(pg. 
e439
-
e442
)
5
Robbins
 
PF
Lu
 
YC
El-Gamil
 
M
et al. 
Mining exomic sequencing data to identify mutated antigens recognized by adoptively transferred tumor-reactive T cells.
Nat Med
2013
, vol. 
19
 
6
(pg. 
747
-
752
)
6
Tran
 
E
Turcotte
 
S
Gros
 
A
et al. 
Cancer immunotherapy based on mutation-specific CD4+ T cells in a patient with epithelial cancer.
Science
2014
, vol. 
344
 
6184
(pg. 
641
-
645
)
7
Schumacher
 
T
Bunse
 
L
Pusch
 
S
et al. 
A vaccine targeting mutant IDH1 induces antitumour immunity.
Nature
2014
, vol. 
512
 
7514
(pg. 
324
-
327
)
8
Snyder
 
A
Makarov
 
V
Merghoub
 
T
et al. 
Genetic basis for clinical response to CTLA-4 blockade in melanoma.
N Engl J Med
2014
, vol. 
371
 
23
(pg. 
2189
-
2199
)
9
Snyder
 
A
Chan
 
TA
Immunogenic peptide discovery in cancer genomes.
Curr Opin Genet Dev
2015
, vol. 
30
 (pg. 
7
-
16
)
10
Rizvi
 
NA
Hellmann
 
MD
Snyder
 
A
et al. 
Cancer immunology. Mutational landscape determines sensitivity to PD-1 blockade in non-small cell lung cancer.
Science
2015
, vol. 
348
 
6230
(pg. 
124
-
128
)
11
Linnemann
 
C
van Buuren
 
MM
Bies
 
L
et al. 
High-throughput epitope discovery reveals frequent recognition of neo-antigens by CD4+ T cells in human melanoma.
Nat Med
2015
, vol. 
21
 
1
(pg. 
81
-
85
)
12
Berlin
 
C
Kowalewski
 
DJ
Schuster
 
H
et al. 
Mapping the HLA ligandome landscape of acute myeloid leukemia: a targeted approach toward peptide-based immunotherapy.
Leukemia
2015
, vol. 
29
 
3
(pg. 
647
-
659
)
13
Kowalewski
 
DJ
Schuster
 
H
Backert
 
L
et al. 
HLA ligandome analysis identifies the underlying specificities of spontaneous antileukemia immune responses in chronic lymphocytic leukemia (CLL).
Proc Natl Acad Sci USA
2015
, vol. 
112
 
2
(pg. 
E166
-
E175
)
14
Kuehl
 
WM
Bergsagel
 
PL
Multiple myeloma: evolving genetic events and host interactions.
Nat Rev Cancer
2002
, vol. 
2
 
3
(pg. 
175
-
187
)
15
Röllig
 
C
Knop
 
S
Bornhauser
 
M
Multiple myeloma.
Lancet
2015
, vol. 
385
 
9983
(pg. 
2197
-
2208
)
16
Barlogie
 
B
Tricot
 
GJ
van Rhee
 
F
et al. 
Long-term outcome results of the first tandem autotransplant trial for multiple myeloma.
Br J Haematol
2006
, vol. 
135
 
2
(pg. 
158
-
164
)
17
Ferrero
 
S
Ladetto
 
M
Drandi
 
D
et al. 
Long-term results of the GIMEMA VEL-03-096 trial in MM patients receiving VTD consolidation after ASCT: MRD kinetics’ impact on survival.
Leukemia
2015
, vol. 
29
 
3
(pg. 
689
-
695
)
18
Martinez-Lopez
 
J
Lahuerta
 
JJ
Pepin
 
F
et al. 
Prognostic value of deep sequencing method for minimal residual disease detection in multiple myeloma.
Blood
2014
, vol. 
123
 
20
(pg. 
3073
-
3079
)
19
Björkstrand
 
B
Iacobelli
 
S
Hegenbart
 
U
et al. 
Tandem autologous/reduced-intensity conditioning allogeneic stem-cell transplantation versus autologous transplantation in myeloma: long-term follow-up.
J Clin Oncol
2011
, vol. 
29
 
22
(pg. 
3016
-
3022
)
20
El-Cheikh
 
J
Crocchiolo
 
R
Furst
 
S
et al. 
Long-term outcome after allogeneic stem-cell transplantation with reduced-intensity conditioning in patients with multiple myeloma.
Am J Hematol
2013
, vol. 
88
 
5
(pg. 
370
-
374
)
21
Koehne
 
G
Giralt
 
S
Allogeneic hematopoietic stem cell transplantation for multiple myeloma: curative but not the standard of care.
Curr Opin Oncol
2012
, vol. 
24
 
6
(pg. 
720
-
726
)
22
Riley
 
JL
Combination checkpoint blockade--taking melanoma immunotherapy to the next level.
N Engl J Med
2013
, vol. 
369
 
2
(pg. 
187
-
189
)
23
Perez
 
SA
von Hofe
 
E
Kallinteris
 
NL
et al. 
A new era in anticancer peptide vaccines.
Cancer
2010
, vol. 
116
 
9
(pg. 
2071
-
2080
)
24
Rosenblatt
 
J
Bar-Natan
 
M
Munshi
 
NC
Avigan
 
DE
Immunotherapy for multiple myeloma.
Expert Rev Hematol
2014
, vol. 
7
 
1
(pg. 
91
-
96
)
25
Brossart
 
P
Schneider
 
A
Dill
 
P
et al. 
The epithelial tumor antigen MUC1 is expressed in hematological malignancies and is recognized by MUC1-specific cytotoxic T-lymphocytes.
Cancer Res
2001
, vol. 
61
 
18
(pg. 
6846
-
6850
)
26
Zhou
 
FL
Meng
 
S
Zhang
 
WG
et al. 
Peptide-based immunotherapy for multiple myeloma: current approaches.
Vaccine
2010
, vol. 
28
 
37
(pg. 
5939
-
5946
)
27
Hundemer
 
M
Schmidt
 
S
Condomines
 
M
et al. 
Identification of a new HLA-A2-restricted T-cell epitope within HM1.24 as immunotherapy target for multiple myeloma.
Exp Hematol
2006
, vol. 
34
 
4
(pg. 
486
-
496
)
28
Jalili
 
A
Ozaki
 
S
Hara
 
T
et al. 
Induction of HM1.24 peptide-specific cytotoxic T lymphocytes by using peripheral-blood stem-cell harvests in patients with multiple myeloma.
Blood
2005
, vol. 
106
 
10
(pg. 
3538
-
3545
)
29
Chiriva-Internati
 
M
Liu
 
Y
Weidanz
 
JA
et al. 
Testing recombinant adeno-associated virus-gene loading of dendritic cells for generating potent cytotoxic T lymphocytes against a prototype self-antigen, multiple myeloma HM1.24.
Blood
2003
, vol. 
102
 
9
(pg. 
3100
-
3107
)
30
Rew
 
SB
Peggs
 
K
Sanjuan
 
I
et al. 
Generation of potent antitumor CTL from patients with multiple myeloma directed against HM1.24.
Clin Cancer Res
2005
, vol. 
11
 
9
(pg. 
3377
-
3384
)
31
van Rhee
 
F
Szmania
 
SM
Zhan
 
F
et al. 
NY-ESO-1 is highly expressed in poor-prognosis multiple myeloma and induces spontaneous humoral and cellular immune responses.
Blood
2005
, vol. 
105
 
10
(pg. 
3939
-
3944
)
32
Schuberth
 
PC
Jakka
 
G
Jensen
 
SM
et al. 
Effector memory and central memory NY-ESO-1-specific re-directed T cells for treatment of multiple myeloma.
Gene Ther
2013
, vol. 
20
 
4
(pg. 
386
-
395
)
33
Bae
 
J
Tai
 
YT
Anderson
 
KC
Munshi
 
NC
Novel epitope evoking CD138 antigen-specific cytotoxic T lymphocytes targeting multiple myeloma and other plasma cell disorders.
Br J Haematol
2011
, vol. 
155
 
3
(pg. 
349
-
361
)
34
Bae
 
J
Carrasco
 
R
Lee
 
AH
et al. 
Identification of novel myeloma-specific XBP1 peptides able to generate cytotoxic T lymphocytes: a potential therapeutic application in multiple myeloma.
Leukemia
2011
, vol. 
25
 
10
(pg. 
1610
-
1619
)
35
Bae
 
J
Smith
 
R
Daley
 
J
et al. 
Myeloma-specific multiple peptides able to generate cytotoxic T lymphocytes: a potential therapeutic application in multiple myeloma and other plasma cell disorders.
Clin Cancer Res
2012
, vol. 
18
 
17
(pg. 
4850
-
4860
)
36
Oka
 
Y
Tsuboi
 
A
Fujiki
 
F
et al. 
WT1 peptide vaccine as a paradigm for “cancer antigen-derived peptide”-based immunotherapy for malignancies: successful induction of anti-cancer effect by vaccination with a single kind of WT1 peptide.
Anticancer Agents Med Chem
2009
, vol. 
9
 
7
(pg. 
787
-
797
)
37
Kuball
 
J
de Boer
 
K
Wagner
 
E
et al. 
Pitfalls of vaccinations with WT1-, Proteinase3- and MUC1-derived peptides in combination with MontanideISA51 and CpG7909.
Cancer Immunol Immunother
2011
, vol. 
60
 
2
(pg. 
161
-
171
)
38
Greiner
 
J
Schmitt
 
A
Giannopoulos
 
K
et al. 
High-dose RHAMM-R3 peptide vaccination for patients with acute myeloid leukemia, myelodysplastic syndrome and multiple myeloma.
Haematologica
2010
, vol. 
95
 
7
(pg. 
1191
-
1197
)
39
Schmitt
 
M
Schmitt
 
A
Rojewski
 
MT
et al. 
RHAMM-R3 peptide vaccination in patients with acute myeloid leukemia, myelodysplastic syndrome, and multiple myeloma elicits immunologic and clinical responses.
Blood
2008
, vol. 
111
 
3
(pg. 
1357
-
1365
)
40
Rapoport
 
AP
Aqui
 
NA
Stadtmauer
 
EA
et al. 
Combination immunotherapy using adoptive T-cell transfer and tumor antigen vaccination on the basis of hTERT and survivin after ASCT for myeloma.
Blood
2011
, vol. 
117
 
3
(pg. 
788
-
797
)
41
Hobo
 
W
Strobbe
 
L
Maas
 
F
et al. 
Immunogenicity of dendritic cells pulsed with MAGE3, Survivin and B-cell maturation antigen mRNA for vaccination of multiple myeloma patients.
Cancer Immunol Immunother
2013
, vol. 
62
 
8
(pg. 
1381
-
1392
)
42
Wang
 
L
Jin
 
N
Schmitt
 
A
et al. 
T cell-based targeted immunotherapies for patients with multiple myeloma.
Int J Cancer
2015
, vol. 
136
 
8
(pg. 
1751
-
1768
)
43
Goswami
 
M
Hensel
 
N
Smith
 
BD
et al. 
Expression of putative targets of immunotherapy in acute myeloid leukemia and healthy tissues.
Leukemia
2014
, vol. 
28
 
5
(pg. 
1167
-
1170
)
44
Chng
 
WJ
Dispenzieri
 
A
Chim
 
CS
et al. 
IMWG consensus on risk stratification in multiple myeloma.
Leukemia
2014
, vol. 
28
 
2
(pg. 
269
-
277
)
45
Kowalewski
 
DJ
Stevanović
 
S
Biochemical large-scale identification of MHC class I ligands.
Methods Mol Biol
2013
, vol. 
960
 (pg. 
145
-
157
)
46
Käll
 
L
Canterbury
 
JD
Weston
 
J
Noble
 
WS
MacCoss
 
MJ
Semi-supervised learning for peptide identification from shotgun proteomics datasets.
Nat Methods
2007
, vol. 
4
 
11
(pg. 
923
-
925
)
47
Schuler
 
MM
Nastke
 
MD
Stevanovikć
 
S
SYFPEITHI: database for searching and T-cell epitope prediction.
Methods Mol Biol
2007
, vol. 
409
 (pg. 
75
-
93
)
48
Sturm
 
T
Leinders-Zufall
 
T
Maček
 
B
et al. 
Mouse urinary peptides provide a molecular basis for genotype discrimination by nasal sensory neurons.
Nat Commun
2013
, vol. 
4
 pg. 
1616
 
49
Garboczi
 
DN
Hung
 
DT
Wiley
 
DC
HLA-A2-peptide complexes: refolding and crystallization of molecules expressed in Escherichia coli and complexed with single antigenic peptides.
Proc Natl Acad Sci USA
1992
, vol. 
89
 
8
(pg. 
3429
-
3433
)
50
Widenmeyer
 
M
Griesemann
 
H
Stevanović
 
S
et al. 
Promiscuous survivin peptide induces robust CD4+ T-cell responses in the majority of vaccinated cancer patients.
Int J Cancer
2012
, vol. 
131
 
1
(pg. 
140
-
149
)
51
Britten
 
CM
Gouttefangeas
 
C
Welters
 
MJ
et al. 
The CIMT-monitoring panel: a two-step approach to harmonize the enumeration of antigen-specific CD8+ T lymphocytes by structural and functional assays.
Cancer Immunol Immunother
2008
, vol. 
57
 
3
(pg. 
289
-
302
)
52
Rudolf
 
D
Silberzahn
 
T
Walter
 
S
et al. 
Potent costimulation of human CD8 T cells by anti-4-1BB and anti-CD28 on synthetic artificial antigen presenting cells.
Cancer Immunol Immunother
2008
, vol. 
57
 
2
(pg. 
175
-
183
)
53
Bui
 
HH
Sidney
 
J
Dinh
 
K
Southwood
 
S
Newman
 
MJ
Sette
 
A
Predicting population coverage of T-cell epitope-based diagnostics and vaccines.
BMC Bioinformatics
2006
, vol. 
7
 pg. 
153
 
54
Schipper
 
RF
van Els
 
CA
D’Amaro
 
J
Oudshoorn
 
M
Minimal phenotype panels. A method for achieving maximum population coverage with a minimum of HLA antigens.
Hum Immunol
1996
, vol. 
51
 
2
(pg. 
95
-
98
)
55
Kanehisa
 
M
Goto
 
S
KEGG: kyoto encyclopedia of genes and genomes.
Nucleic Acids Res
2000
, vol. 
28
 
1
(pg. 
27
-
30
)
56
Huang da
 
W
Sherman
 
BT
Lempicki
 
RA
Systematic and integrative analysis of large gene lists using DAVID bioinformatics resources.
Nat Protoc
2009
, vol. 
4
 
1
(pg. 
44
-
57
)
57
Braumüller
 
H
Wieder
 
T
Brenner
 
E
et al. 
T-helper-1-cell cytokines drive cancer into senescence.
Nature
2013
, vol. 
494
 
7437
(pg. 
361
-
365
)
58
Schoenberger
 
SP
Toes
 
RE
van der Voort
 
EI
Offringa
 
R
Melief
 
CJ
T-cell help for cytotoxic T lymphocytes is mediated by CD40-CD40L interactions.
Nature
1998
, vol. 
393
 
6684
(pg. 
480
-
483
)
59
Janssen
 
EM
Lemmens
 
EE
Wolfe
 
T
Christen
 
U
von Herrath
 
MG
Schoenberger
 
SP
CD4+ T cells are required for secondary expansion and memory in CD8+ T lymphocytes.
Nature
2003
, vol. 
421
 
6925
(pg. 
852
-
856
)
60
Greiner
 
J
Ono
 
Y
Hofmann
 
S
et al. 
Mutated regions of nucleophosmin 1 elicit both CD4(+) and CD8(+) T-cell responses in patients with acute myeloid leukemia.
Blood
2012
, vol. 
120
 
6
(pg. 
1282
-
1289
)
61
Wolchok
 
JD
Kluger
 
H
Callahan
 
MK
et al. 
Nivolumab plus ipilimumab in advanced melanoma.
N Engl J Med
2013
, vol. 
369
 
2
(pg. 
122
-
133
)
62
Topalian
 
SL
Hodi
 
FS
Brahmer
 
JR
et al. 
Safety, activity, and immune correlates of anti-PD-1 antibody in cancer.
N Engl J Med
2012
, vol. 
366
 
26
(pg. 
2443
-
2454
)
63
Hodi
 
FS
O’Day
 
SJ
McDermott
 
DF
et al. 
Improved survival with ipilimumab in patients with metastatic melanoma.
N Engl J Med
2010
, vol. 
363
 
8
(pg. 
711
-
723
)
64
Robert
 
C
Ribas
 
A
Wolchok
 
JD
et al. 
Anti-programmed-death-receptor-1 treatment with pembrolizumab in ipilimumab-refractory advanced melanoma: a randomised dose-comparison cohort of a phase 1 trial.
Lancet
2014
, vol. 
384
 
9948
(pg. 
1109
-
1117
)
65
Brahmer
 
JR
Tykodi
 
SS
Chow
 
LQ
et al. 
Safety and activity of anti-PD-L1 antibody in patients with advanced cancer.
N Engl J Med
2012
, vol. 
366
 
26
(pg. 
2455
-
2465
)
66
Hamid
 
O
Robert
 
C
Daud
 
A
et al. 
Safety and tumor responses with lambrolizumab (anti-PD-1) in melanoma.
N Engl J Med
2013
, vol. 
369
 
2
(pg. 
134
-
144
)
67
Motzer
 
RJ
Rini
 
BI
McDermott
 
DF
et al. 
Nivolumab for metastatic renal cell carcinoma: results of a randomized phase II trial.
J Clin Oncol
2015
, vol. 
33
 
13
(pg. 
1430
-
1437
)
68
Lynch
 
TJ
Bondarenko
 
I
Luft
 
A
et al. 
Ipilimumab in combination with paclitaxel and carboplatin as first-line treatment in stage IIIB/IV non-small-cell lung cancer: results from a randomized, double-blind, multicenter phase II study.
J Clin Oncol
2012
, vol. 
30
 
17
(pg. 
2046
-
2054
)
69
Ansell
 
SM
Hurvitz
 
SA
Koenig
 
PA
et al. 
Phase I study of ipilimumab, an anti-CTLA-4 monoclonal antibody, in patients with relapsed and refractory B-cell non-Hodgkin lymphoma.
Clin Cancer Res
2009
, vol. 
15
 
20
(pg. 
6446
-
6453
)
70
Ansell
 
SM
Lesokhin
 
AM
Borrello
 
I
et al. 
PD-1 blockade with nivolumab in relapsed or refractory Hodgkin’s lymphoma.
N Engl J Med
2015
, vol. 
372
 
4
(pg. 
311
-
319
)
71
Bassani-Sternberg
 
M
Pletscher-Frankild
 
S
Jensen
 
LJ
Mann
 
M
Mass spectrometry of human leukocyte antigen class I peptidomes reveals strong effects of protein abundance and turnover on antigen presentation.
Mol Cell Proteomics
2015
, vol. 
14
 
3
(pg. 
658
-
673
)
72
Stickel
 
JS
Weinzierl
 
AO
Hillen
 
N
et al. 
HLA ligand profiles of primary renal cell carcinoma maintained in metastases.
Cancer Immunol Immunother
2009
, vol. 
58
 
9
(pg. 
1407
-
1417
)
73
Weinzierl
 
AO
Lemmel
 
C
Schoor
 
O
et al. 
Distorted relation between mRNA copy number and corresponding major histocompatibility complex ligand density on the cell surface.
Mol Cell Proteomics
2007
, vol. 
6
 
1
(pg. 
102
-
113
)
74
Min
 
DJ
Ezponda
 
T
Kim
 
MK
et al. 
MMSET stimulates myeloma cell growth through microRNA-mediated modulation of c-MYC.
Leukemia
2013
, vol. 
27
 
3
(pg. 
686
-
694
)
75
Martinez-Garcia
 
E
Popovic
 
R
Min
 
DJ
et al. 
The MMSET histone methyl transferase switches global histone methylation and alters gene expression in t(4;14) multiple myeloma cells.
Blood
2011
, vol. 
117
 
1
(pg. 
211
-
220
)
76
Keats
 
JJ
Maxwell
 
CA
Taylor
 
BJ
et al. 
Overexpression of transcripts originating from the MMSET locus characterizes all t(4;14)(p16;q32)-positive multiple myeloma patients.
Blood
2005
, vol. 
105
 
10
(pg. 
4060
-
4069
)
77
Brito
 
JL
Walker
 
B
Jenner
 
M
et al. 
MMSET deregulation affects cell cycle progression and adhesion regulons in t(4;14) myeloma plasma cells.
Haematologica
2009
, vol. 
94
 
1
(pg. 
78
-
86
)
78
Smith
 
EM
Boyd
 
K
Davies
 
FE
The potential role of epigenetic therapy in multiple myeloma.
Br J Haematol
2010
, vol. 
148
 
5
(pg. 
702
-
713
)
79
Xie
 
Z
Gunaratne
 
J
Cheong
 
LL
et al. 
Plasma membrane proteomics identifies biomarkers associated with MMSET overexpression in T(4;14) multiple myeloma.
Oncotarget
2013
, vol. 
4
 
7
(pg. 
1008
-
1018
)
80
Hébraud
 
B
Caillot
 
D
Corre
 
J
et al. 
The translocation t(4;14) can be present only in minor subclones in multiple myeloma.
Clin Cancer Res
2013
, vol. 
19
 
17
(pg. 
4634
-
4637
)

Author notes

S.W. and J.S.S. contributed equally to this work.

Sign in via your Institution