Despite remarkable therapeutic advances in recent years, multiple myeloma (MM) remains incurable, hence representing an unmet clinical need. Therapeutic approaches based on induced T cell anti-tumor immunity towards cancer eradication show promising results in many hematologic malignancies, including multiple myeloma (MM), highlighting the need for comprehensive understanding of the implicated mechanisms. Similar to other mature B cell malignancies, immunoglobulin (IG) gene rearrangements in MM lead to the expression of unique, novel peptide sequences that can be used as neoantigens, conceivably representing a potential attractive target for cancer-specific immune responses. Here we sought to explore this possibility by performing ad hoc prediction of putative T-cell class I/II neoepitopes contained within the MM clonotypic IG gene rearrangements. To that end, we investigated 16 newly diagnosed patients with MM from whom we isolated the peripheral blood and/or bone marrow mononuclear cell (PBMC/BMMC) fraction as well as CD138+ myeloma cells, CD4+ and CD8+ T cells. Clonotypic IG heavy and light chain gene rearrangements were RT-PCR amplified on RNA extracted from the PBMC/BMMC fraction using subgroup-specific leader primers for the IGHV/IGK/LV gene and universal primers annealing to the constant domain (IGHG, IGHM, IGKC, IGLC), in order to produce the full-length V-(D)-J gene rearrangement sequence, plus the start of the constant domain. The corresponding sequences were determined by bidirectional Sanger sequencing. Moreover, T cell receptor beta (TRB) chain gene rearrangements were RT-PCR amplified on RNA extracted from CD4+ and CD8+ T cells and then subjected to paired-end next generation sequencing (NGS). The deduced amino acid sequences of the IG heavy and light chain variable/constant domains were subsequently parsed in peptides and subjected to bioinformatics analysis for the identification of putative T-cell class I/II neoepitopes using the NetMHCpan and NetMHIICpan softwares. The rank score was calculated, considering the 4-digit HLA-A, -B, -C and DRB1 typing for each individual patient. High- and medium-binding peptides (rank score <2%) were selected. Prediction of the binding specificity of T Cell Receptors (TR) to MHC-peptide complexes (pMHCs) was performed by the ERGO-II tool. The rank score was calculated, considering the peptides resulted by NetMHCpan/ NetMHIICpan, the HLA, and the clonotypes for each patient (AUC>0.8). Exact matches to germline and/or proteome databases were excluded. Overall, 757,263 TRΒ distinct clonotypes were assessed (median= 18,244 /sample). All patients displayed oligoclonal T cell expansions in both CD4+ and CD8+ T cells, albeit these were significantly (p<0.001) more pronounced in the latter. Both T cell subpopulations displayed skewed TRBV and TRBJ gene repertoires (7 TRBV genes and 3 TRBJ genes accounting for about 40% and 59% of the repertoire, respectively). Overall, 1,219 predicted neoepitopes were identified. All patients had predicted CD4+ and CD8 + T-cell epitopes within the MM clonotypic IG, either in the heavy chain (16/16 pts, n=459 epitopes) or the light chain (16/16 pts, n=760 epitopes). There was no statistically significant difference in the rank score of peptides involving the complementarity determining region 3 (CDR3) vs. all other IG regions. Interestingly, 998/1,219 (82%) peptides with strong binding affinity resulted from IG gene sequence motifs outside the CDR3. Most of these peptides (853/1,219, 75%) resulted from somatic hypermutations (SHM) across the IG variable domain. Overall, 51,574 clonotypes of CD4+ T cells and 35,751 clonotypes of CD8+ T cells were in silico predicted by ERGO to have strong binding affinity. Relevant to mention, 528 clonotypes with high binding score were found to be shared between CD4+ T cells of all patients and 209 clonotypes between CD8+ T cells of all patients and, thus, were deemed as public. In conclusion, in silico prediction identified a significant number of putative T-cell class I/II neoepitopes contained within the clonotypic IG of MM patients. Many of the identified peptides derive from SHM, implying that the SHM which shapes the MM BcR IG repertoire may produce immunogenic CD4+/CD8+ T cell epitopes. Their actual immunogenicity has to be tested in ex vivo studies, currently underway by our group.
Read full abstract