Abstract Introduction Tumor specific somatic mutations can potentially influence patient response to targeted immunotherapies. Since somatic mutations generate neoantigens, novel tumor-specific peptides can be processed, loaded into MHC class I molecules, and presented to T-cells as potentially immunogenic targets. Indeed, both somatic mutation burden and neoantigen burden have been correlated with responsiveness to checkpoint inhibitor blockade in highly mutagenic advanced melanoma and lung cancer. However, correlation of response with mutation and neoantigen burden has not been reported for cancers of lower characteristic mutational burden such as multiple myeloma (MM). Indeed early results of checkpoint inhibitors in MM have been underwhelming. We therefore wished to better understand the mutation and neoantigen landscape of MM patients as it relates to responsiveness to current and future therapies. Materials/Methods Somatic mutation burden, neoantigen load, and response to therapy was obtained using interim data from theCoMMpassstudy (NCT01454297), a study by the Multiple Myeloma Research Foundation that collects longitudinal clinical and molecular data on 1,000 newly diagnosed MM patients from baseline through treatment and disease progression. Next-generation exome sequencing data was used to call somatic mutations andbioinformaticallydetermine patient specific HLA type. Neoantigen binding affinity was predicted for all unique patient specific somatic missense mutation and HLA allele combinations using a well validated neural network and machine learning based algorithm. Binding affinity of the corresponding non-mutated peptide was also predicted to identify mutated neoantigens witha high binding affinity while the wild type counterpart has low predicted binding affinity. Results For 627 patients with available sequencing and clinical data, mutation and neoantigen profiles were created. Within this population, the mean overall canonical somatic mutation load was 405.84 (s=608.55) mutations per patient. These mutations included 38,222 unique missense mutations in 12,994 genes. The mean somatic missense mutation load was 64.51(s=95.29) mutations per patient including immunoglobulin genes (52.70, s=94.85 excluding immunoglobulins). There was a positive linear correlation between mutation burden and neoantigen burden. The mean neoantigen load using thresholds of mutant peptide binding affinity (half maximal inhibitory concentration, IC50) less than 500nM was 69.57 (s=124.76) neoantigens (58.53, s=123.58 excluding immunoglobulins). Limiting binding affinity of wild type peptide sequence to greater than 500nM decreased the mean neoantigen load to 23.84 (s=53.55) neoantigens (20.74, s=53.28 excluding immunoglobulins). Progression free survival (PFS) following induction therapy was used to determine correlates of mutation burden excluding immunoglobulin mutations with response. Kaplan-Meier survival analysis revealed significantly shorter PFS time in patients with greater than average tumor somatic missense mutation load (N=126) relative to those with less than average mutation load (49.1% vs. 73.0% PFS at 2yrs, p=0.0055, Figure 1) Similarly, patients with greater than average neoantigen load (N=164) had significantly shorter PFS relative to those with less than average neoantigen load (48.0% vs. 74.2% PFS at 2yrs, p=0.00259, Figure 2). Moreover, patients with both high mutation and neoantigen load (N=85) had significantly shorter PFS relative to those with less than average mutation and neoantigen load (45.3% vs. 72% PFS at 2yrs, p= 0.0168). Conclusion While mutation and neoantigen load has been correlated with response to targeted immunotherapies in highly mutated solid tumors, we report here a negative correlation between mutation, neoantigen load and MM patient survival. This indicates that a high mutation andneoantigen load negatively impacts survival of myeloma patients under current standard of care. It is hypothesized that this is because current treatments fail to recruit T cell immunity necessary to exploit tumor neoantigens fundamental to immune-mediated killing of cancer cells. This presents an encouraging basis for pursuing targeted and personalized immunotherapies for multiple myeloma, especially for patients with a higher mutational burden. PFS and mutation load. PFS and mutation load. PFS and neoantigen load. PFS and neoantigen load. Disclosures Lonial: BMS: Consultancy; Celgene: Consultancy; Janssen: Consultancy; Onyx: Consultancy; Celgene: Consultancy; Novartis: Consultancy; Janssen: Consultancy; Onyx: Consultancy; Novartis: Consultancy; BMS: Consultancy; Millenium: Consultancy; Merck: Consultancy.
HOXA9 Is a Novel Therapeutic Target in Multiple Myeloma.
Abstract Abstract 832 We hypothesized that new therapeutic targets for multiple myeloma (MM) could be discovered through the integrative computational analysis of genomic data. Accordingly, we generated gene expression profiling and copy number data on 250 clinically-annotated MM patient samples. Utilizing an outlier statistical approach, we identified HOXA9 as the top candidate gene for further investigation. HOXA9 expression was particularly high in patients lacking canonical MM chromosomal translocations, and allele-specific expression analysis suggested that this overexpression was mono-allelic. Indeed, focal copy number amplifications at the HOXA locus were observed in some patients. Outlier HOXA9 expression was further validated in both a collection of 52 MM cell lines and 414 primary patient samples previously described. To further verify the aberrant expression of HOXA9 in MM, we performed quantitative RT-PCR, which confirmed expression in all MM patients and cell lines tested, with high-level expression in a subset. To further investigate the mechanism of aberrant HOXA9 expression, we interrogated the pattern of histone modification at the HOXA locus because HOXA gene expression is particularly regulated by such chromatin marks. Accordingly, immunoprecipitation studies showed an aberrantly low level of histone 3 lysine 27 trimethylation marks (H3K27me3) at the HOXA9 locus. H3K27me3 modification is normally associated with silencing of HOXA9 in normal B-cell development. As such, it appears likely that the aberrant expression of HOXA9 in MM is due at least in part to defects in histone modification at this locus. To determine the functional consequences of HOXA9 expression in MM, we performed RNAi-mediated knock-down experiments in MM cell lines. Seven independent HOXA9 shRNAs that diminished HOXA9 expression resulted in growth inhibition of 12/14 MM cell lines tested. Taken together, these experiments indicate that HOXA9 is essential for survival of MM cells, and that the mechanism of HOXA9 expression relates to aberrant histone modification at the HOXA9 locus. The data thus suggest that HOXA9 is an attractive new therapeutic target for MM. Disclosures: No relevant conflicts of interest to declare.
Precision medicine in multiple myeloma: are we there yet?
Expert Review of Precision Medicine and Drug Development | February 7, 2019 | Volume: 4 | Issue: 2 | Pages: 51-53
Key Points Using the largest set of patients with newly diagnosed myeloma, we identified 63 mutated driver genes. We identified oncogenic dependencies, particularly relating to primary translocations, indicating a nonrandom accumulation of genetic hits.
The Role of Minimal Residual Disease Testing in Myeloma Treatment Selection and Drug Development: Current Value and Future Applications
Clinical Cancer Research | April 20, 2017 | Volume: 23 | Issue: 15 | Pages: 3980-3993
To gain insight into the genomic basis of diffuse large B-cell lymphoma (DLBCL), we performed massively parallel whole-exome sequencing of 55 primary tumor samples from patients with DLBCL and matched normal tissue. We identified recurrent mutations in genes that are well known to be functionally relevant in DLBCL, including MYD88 , CARD11 , EZH2 , and CREBBP . We also identified somatic mutations in genes for which a functional role in DLBCL has not been previously suspected. These genes include MEF2B , MLL2 , BTG1 , GNA13 , ACTB , P2RY8 , PCLO , and TNFRSF14 . Further, we show that BCL2 mutations commonly occur in patients with BCL2 / IgH rearrangements as a result of somatic hypermutation normally occurring at the IgH locus. The BCL2 point mutations are primarily synonymous, and likely caused by activation-induced cytidine deaminase–mediated somatic hypermutation, as shown by comprehensive analysis of enrichment of mutations in WRCY target motifs. Those nonsynonymous mutations that are observed tend to be found outside of the functionally important BH domains of the protein, suggesting that strong negative selection against BCL2 loss-of-function mutations is at play. Last, by using an algorithm designed to identify likely functionally relevant but infrequent mutations, we identify KRAS , BRAF , and NOTCH1 as likely drivers of DLBCL pathogenesis in some patients. Our data provide an unbiased view of the landscape of mutations in DLBCL, and this in turn may point toward new therapeutic strategies for the disease.
Abstract 4822: Analysis of somatic copy number alterations associated with poor prognosis and progression of multiple myeloma in African Americans
Cancer Research | February 25, 2011 | Volume: 70 | Issue: 8_Supplement | Pages: 4822-4822
Abstract The plasma cell malignancy Multiple Myeloma (MM) is a rare but deadly form of cancer with 5 year survival rates of ∼30%. The disease exacts an especially heavy burden on African American (AA) patients. In fact, African Americans are twice a likely to be diagnosed with and die from MM as compared to Caucasian Americans. Although socio-economic factors likely influence this disparity, the role of biology has not been ruled out. We have profiled ∼250 MM tumors for somatic copy number alterations using array-based Comparative Genomic Hybridization (aCGH) and Gene Expression Profiling as part of the Multiple Myeloma Genomics Initiative. To determine if somatic alterations may play a role in the disparity in incidence and/or severity of MM in AA, we compared the frequency of specific chromosomal aberrations in genes and genomic regions previously found to be involved in progression and poor outcome in MM among a cohort of one hundred ninety-six (196) patients including 16 African Americans and 180 Caucasians. The frequency of occurrence (expressed as a percentage) was determined for each of the following aberrations: 1p deletion, homozygous deletion of p18, chromosome 13 deletion, and chromosome 13 monosomy. No significant difference was detected in the frequency of chromosome 13 deletion and chromosome 13 monosomy between the two groups. However, significant differences were observed for chromosome 1p alterations. Homozygous deletion at the p18 locus at 1p occurred in 12.5% (2/16) of tumors from AA patients and 2.8% (5/180) of CA patient tumors (p = 0.045). We further stratified tumors by hyperdiploid or non-hyperdiploid status using a-CGH profiles to determine ploidy status. Most significantly, p18 homozygous deletion was found in 25% (2/8) of AA patients with non-hyperdiploid tumors and in only 1.2% (1/83) of CA patients with non-hyperdiploid tumors (p=3.2 × 10-4). A nominally significant difference among non-hyperdiploid patients was also present for overall 1p deletion, which was observed in 2/8 (25%) of AA non-hyperdiploid tumors and 33/83 (40%) of CA with non-hyperdiploid tumors (p=0.05). No differences were observed for 1p deletion or p18 homozygous deletion among hyperdiploid tumors between AA and CA patients. These results suggest that deletions at 1p are more common among African American MM patients. Furthermore, homozygous deletion encompassing the p18 locus is more significantly common among AA MM patients, especially AA patients with non-hyperdiploid tumors. As 1p deletion and p18 homozygous deletion have been linked to poor outcome, our results suggest that the increased frequency of these alterations in tumors from AA patients may in part explain the disparity in unfavorable outcomes seen among AA MM patients Citation Format: {Authors}. {Abstract title} [abstract]. In: Proceedings of the 101st Annual Meeting of the American Association for Cancer Research; 2010 Apr 17-21; Washington, DC. Philadelphia (PA): AACR; Cancer Res 2010;70(8 Suppl):Abstract nr 4822.
Abstract B64: Analysis of somatic copy number alterations among African American patients with multiple myeloma using array CGH
Cancer Epidemiology, Biomarkers & Prevention | December 6, 2010 | Volume: 19 | Issue: 10_Supplement | Pages: B64-B64
Abstract Multiple myeloma (MM) is a cancer of the plasma cells in the bone marrow. Although rare, MM is a deadly form of cancer with a 5 year survival rate of only 37%. African Americans (AA) are twice as likely to be diagnosed with and die from MM as compared to European Americans (EA). Socioeconomic factors may contribute to this disparity, however, biological factors may also be an influence. We have identified somatic copy number alterations in MM as a part of the Multiple Myeloma Genomics Initiative by profiling 250 MM tumors using array-based comparative genomic hybridization (aCGH). Using this cohort, we identified 16 African American and 180 European American patients using a MMRC database that indicated the racial demographics of each patient. To determine if somatic alterations may play a role in the disparity in incidence and/or severity of MM in AA, we compared the frequency of specific chromosomal aberrations in genes and genomic regions found to be involved in progression and poor outcome. The number of events of occurring was calculated and a univariate analysis was conducted using the Fisher's exact test. P-values were used to determine statistical significance. Results revealed significant differences for chromosome 1 p alterations. Specifically, homozygous deletion at the p18 locus at 1 p occurred in 12.5% (2/16) of tumors from AA patients and 2.8% (5/180) of EA patient tumors (p = 0.045). Further stratification of tumors by hyperdiploid or non-hyperdiploid status revealed significant differences in p18 homozygous deletion. Homozygous deletion of p18 was found in 25% (2/8) of AA patients and only 1.2% (1/83) of EA patients with non-hyperdiploid tumors (p=3.2 × 10-4). A nominally significant difference among non-hyperdiploid patients was also present for overall 1 p deletion, which was observed in 2/8 (25%) of AA non-hyperdiploid tumors and 33/83 (40%) of EA with non-hyperdiploid tumors (p=0.05). No differences were observed for 1 p deletion or p18 homozygous deletion among hyperdiploid tumors between AA and EA patients. We then combined our AA dataset with an independent dataset of 8 AA MM patients generated at the Mayo Clinic, Scottsdale AZ and repeated the analysis. Interestingly, very little or no aberrations were found in the dataset from the Mayo Clinic, however, differences of p18 homozygous deletion and 16q deletion among AA and EA remained statistically significant. To determine the correlation between hemizygous and biallelic loss of p18 and poor prognosis we used the 70 gene signature (Shaughnessy, et al., 2007) as a surrogate to identify patients with poor prognosis within the Multiple Myeloma Genomics Initiative dataset. Genomic identification of significant targets in cancer (GISTIC) was used (Beroukhim et al., 2007) to identify regions of aberrant copy number that were over-represented among samples. A relationship between poor prognosis and biallelic loss was found to be statistically significant (p = 0002) while samples with hemizygous loss showed no relationship (p = 0.57). Although the size of the current AA cohort is small, these results suggest that homozygous deletions at the p18 locus may be more common among African American MM patients with non-hyperdiploid tumors which may help to explain the disparity in AA patients with MM. Array CGH is currently being conducted on an additional 40 tumors from AA MM patients for a more robust analysis. Citation Information: Cancer Epidemiol Biomarkers Prev 2010;19(10 Suppl):B64.
DNA Methylation Analysis Determines the High Frequency of Genic Hypomethylation and Low Frequency of Hypermethylation Events in Plasma Cell Tumors
Cancer Research | August 24, 2010 | Volume: 70 | Issue: 17 | Pages: 6934-6944
AbstractGenetic heterogeneity between individuals confounds the comparison of gene profiling of multiple myeloma (MM) cells versus normal plasma cells (PCs). To overcome this barrier, we compared the gene expression profile of CD138+ MM cells from a patient bone marrow (BM) sample with CD138+ PCs from a genetically identical twin BM sample using microarray profiling. Two hundred and ninety-six genes were up-regulated and 103 genes were down-regulated at least 2-fold in MM cells versus normal twin PCs. Highly expressed genes in MM cells included cell survival pathway genes such as mcl-1, dad-1, caspase 8, and FADD-like apoptosis regulator (FLIP); oncogenes/transcriptional factors such as Jun-D, Xbp-1, calmodulin, Calnexin, and FGFR-3; stress response and ubiquitin/proteasome pathway–related genes and various ribosomal genes reflecting increased metabolic and translational activity. Genes that were down-regulated in MM cells versus healthy twin PCs included RAD51, killer cell immunoglobulin-like receptor protein, and apoptotic protease activating factor. Microarray results were further confirmed by Western blot analyses, immunohistochemistry, fluorescent in situ hybridization (FISH), and functional assays of telomerase activity and bone marrow angiogenesis. This molecular profiling provides potential insights into mechanisms of malignant transformation in MM. For example, FGFR3, xbp-1, and both mcl-1 and dad-1 may mediate transformation, differentiation, and survival, respectively, and may have clinical implications. By identifying genes uniquely altered in MM cells compared with normal PCs in an identical genotypic background, the current study provides the framework to identify novel therapeutic targets.