Keats, Jonathan J.

Abstract Multiple myeloma (MM) is the second most common hematologic cancer, accounting for roughly 2% of all cancer deaths. Unfortunately, MM remains incurable and nearly all patients experience relapse. Recently, we have identified a novel cIAP1 antagonist, LCL161, with remarkable anti-MM effects. Interestingly, unlike in other tumor types, our previous work found that LCL161 does not directly induce apoptosis of MM cells, but rather drives activation and infiltration of macrophages to the tumor bed and phagocytosis of tumor cells. Furthermore, in vitro experiments demonstrate that the anti-MM activity of the agent relies on the presence of the tumor, with no observed increase in MM phagocytosis when macrophages are treated in isolation. However, LCL161 is a strong and indiscriminate transcriptional activator and thus, the molecular mechanisms underlying this immune response in MM remain poorly understood. In an effort to characterize molecular signatures associated with this response, we collected whole spleens from seven WT mice and eleven mice with VK12598 transplanted myeloma for single cell RNA-seq. Of the eleven tumor-bearing mice, three remained untreated, three were treated with LCL161, three with bortezomib, and three with cyclophosphamide. Additionally, one of the WT mice was treated with LCL161. Using the 10X genomics chromium platform, we measured mRNA levels in over 90,000 individual cells. We identified groups of cells with shared expression profiles, corresponding to known cell types, using graph-based clustering. Consistent with our previous work, we observed a marked reduction in tumor cells and activation of the noncanonical NFkB pathway in mice treated with LCL161. Turning our attention to cells identified as monocytes, with a particular focus on putative macrophages, we sought to identify molecular signatures associated with LCL161 treatment. Indeed, using a binomial model we identified over 500 upregulated genes in LCL161-treated macrophages. However, our previous data suggested that the efficacy of LCL161 was dependent on the presence of tumor cells. Thus, we refined these results by comparing the expression profile of macrophages treated with LCL161 in WT and tumor-bearing mice. These contrasts identified dozens of genes that were significantly upregulated in macrophages from LCL161 treated, tumor-bearing mice. Of note, we identified a small number of genes (including Spp1, Irg1, and Lipg) that were active only in macrophages from mice treated with LCL161 and bearing a transplanted tumor and unexpressed in all other cells. We believe these genes—and the pathways in which they are active—hold promise for the development of novel targeted therapies to better treat MM. Citation Format: Nicholas Banovich, Martin Boateng, Marta Chesi, Leif Bergsagel, Jonathan Keats. Single-cell RNA sequencing identifies macrophage-specific expression signatures associated with phagocytosis of multiple myeloma after treatment with cIAP antagonist [abstract]. In: Proceedings of the American Association for Cancer Research Annual Meeting 2018; 2018 Apr 14-18; Chicago, IL. Philadelphia (PA): AACR; Cancer Res 2018;78(13 Suppl):Abstract nr 3135.

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.

Abstract The gene KDM6A/UTX was recently identified as a putative tumor suppressor (van Haaften et al). The gene is frequently lost in a variety of malignancies but most commonly in esophageal squamous cell carcinomas (8%) and multiple myeloma (10%). The initial characterization of UTX abnormalities in multiple myeloma utilized a mixture of 58 cell lines and patient samples. To establish the frequency of UTX inactivation in multiple myeloma we studied the copy number profiling dataset (n=239) available from the Multiple Myeloma Genomics Initiative (MMGI) and our internal dataset of 68 human myeloma cell lines which were characterized by array-based CGH (Agilent 400k), whole exome sequencing, and mRNA sequencing. We identified homozygous deletions in 2.5% of patients and 26.5% of cell lines. There are no mutations in the available whole genome and exome sequencing data from 38 patients in the MMGI that partially overlaps the 239 copy number samples sequencing. However, in the cell lines another 14.3% of samples showed loss of UTX through a variety of mechanisms (missense, nonsense, frameshifts, single exon deletions and insertions of LINE elements) detected by exome and mRNA sequencing, bringing the total percentage of cell lines with UTX inactivation to 40.8%. The increased frequency of UTX inactivation in cell lines likely suggests that it is associated with a more aggressive disease or promotes the growth of the tumor outside of the normal bone marrow microenvironment. As part of the cell line characterization phase of the project we found that ARD has a homozygous deletion encompassing the entire UTX locus while its isogenic counterparts, ARP-1 and CAG, express wild-type UTX at a level comparable to other cell lines without UTX abnormalities. To ensure these isogenic cell lines represent the same tumor not different tumors from the same individual we cloned the immunoglobulin receptor rearrangements in all there lines, which confirmed these are derived from the same original tumor. To develop this model system we first generated single cell clones from each of the three cell lines. Six clones from each were randomly selected and their DNA content was determined by flow cytometry. In all cases the clones matched the parental population in ploidy (ARD 3.2N, ARP 3.3N, CAG 3.8N). We randomly select one subclone from each cell line to characterize in detail using spectral karyotyping, array-based CGH, whole exome sequencing, mRNA sequencing, and CHIP-seq for H3K4me3, H3K27ac, H3K27me1, H3K27me2, and H3K27me3. The detailed characterization of this isogenic model system will make it an excellent platform for future work aimed at understanding how UTX loss contributes to myelomagenesis. Citation Format: Kristi Allen, Venkata Yellapantula, Quan Zhao, David Edwards, Marta Chesi, Leif Bergsagel, John Carpten, Jonathan Keats. Characterization of an isogenic model system for KDM6A/UTX loss in multiple myeloma. [abstract]. In: Proceedings of the 104th Annual Meeting of the American Association for Cancer Research; 2013 Apr 6-10; Washington, DC. Philadelphia (PA): AACR; Cancer Res 2013;73(8 Suppl):Abstract nr 5209. doi:10.1158/1538-7445.AM2013-5209

Abstract Previous epidemiological studies have suggested that African Americans (AA) are twice as likely to be diagnosed with and to die from Multiple Myeloma (MM) as compared to European Americans (EA). In recent years, however, mortality rates have become more comparable between these groups. Although socioeconomic factors could play a major role to this potential disparity, biological factors may also be responsible for influencing these differences in either incidence and/or mortality. Here we examined a multi-ethnic cohort of MM tumors to assess potential differences in the frequency of molecular events in tumors derived from either African American (AA) or European American (EA) patients. First, we compared the frequency of 14q32 translocation breakpoints at the IgH locus among a series of data from 115 AA MM patients from three studies and EA MM from the Eastern Cooperative Oncology Group (ECOG). Furthermore, we analyzed matching genome-wide copy number and transcriptional data among 47 AA MM tumors and 196 EA MM tumors to determine differences in the frequency of molecular events in tumors from these populations. Univariate analyses were conducted using the Fisher's Exact Test and p-values were calculated to measure statistical significance. We observed statistically significant differences in the frequency of IgH translocations and somatic copy number alterations. Although we did not detect differences for specific translocation subtypes, we saw a statistically significant difference in the overall frequency of IgH translocations between the two groups, which occurred more frequently in EA patients (p = 0.032). When assessing genomic copy number events previously shown to be associated with poor outcome in MM, chromosome 13 monosomy was more frequent in hyperdiploid MM from AA patients and significant after correcting for multiple testing using the Bonferoni correction method (p = 0.001). Moreover, when using the Benjamini-Hochberg multiple testing correction method, chromosome 13 monosomy remained significantly different among hyperdiploid AA patients (p = 0.0121). Chromosome 1q amplification was also significant among hyperdiploid tumors from EA patients when using this correction method (p = 0.036). No differences in the frequency of hyperdiploid status or gene expression-based subtypes were detected when comparing tumors from AA and EA patients. Our study represents the first comprehensive assessment of the frequency and distribution of molecular alterations in MM tumors from both AA and EA patients and could help shed light on possible biological features associated with population differences in incidence and outcome in MM. Citation Format: {Authors}. {Abstract title} [abstract]. In: Proceedings of the 103rd Annual Meeting of the American Association for Cancer Research; 2012 Mar 31-Apr 4; Chicago, IL. Philadelphia (PA): AACR; Cancer Res 2012;72(8 Suppl):Abstract nr 5064. doi:1538-7445.AM2012-5064

Abstract Background: In order to understand the genetic events associated with disease progression and the development of drug resistance in multiple myeloma, we studied four longitudinally collected samples from a single, high risk, t(4;14) multiple myeloma patient. Samples over a four-year period were available from diagnosis, at time of first progression on lenalidomide, second progression while on bortezomib, at time of relapsed refractory disease and with end stage plasma cell leukemia. To better understand these genetic events, we utilized DNA extracted from germline tissue and CD138 purified bone marrow tumor cells to conduct whole genome sequencing (WGS) utilizing the Life Technologies SOLiD™ platform. Results: We have sequenced over 100 billion bases per sample with whole genome fragment library sequencing providing an average of 30X aligned coverage at all four time points studied. Each sample was aligned to the reference genome (NCBI36/hg18) with BFAST and post-alignment removal of duplicates and recalibration performed with Picard and Genome Analysis Toolkit (GATK) respectively. We then called point mutations with SolSNP and compared the germline sample to each tumor sample utilizing the Paired Mutation Walker. To date, we have identified 224 nonsynonymous point mutations within 217 genes. We then utilized the Sorting Intolerant from Tolerant (SIFT) algorithm and identified 109 potentially damaging point mutations which include 108 novel mutations in coding exons. Of the 224 somatic mutations observed, only 20 of these mutations are common to all tumor time points. Interestingly, the diagnostic and bortezomib refractory (second progression) samples had 11 mutations in common while the lenalidomide refractory (first progression) sample shared only 4 mutations in common with the bortezomib refractory sample. This suggests the presence of multiple clones at diagnosis that, due to treatment selection pressure, may emerge or regress as the dominant clone. Preliminary pathway analysis reveals the presence of mutations in genes regulating NFKB and Akt, whose pathways are often dysregulated in cancer. Conclusions: Whole genome sequencing has revealed evidence of mutational differences in exons between tumor and germline DNA indicative of the presence of multiple clones at diagnosis and evolution of these clones over time. This data supports the use of multiple drug therapies that target multiple clones and challenges the long held assumption that a patient who becomes refractory to a particular treatment will always be refractory to that treatment. Additional evaluation in a larger population is needed to confirm these findings. Citation Format: {Authors}. {Abstract title} [abstract]. In: Proceedings of the 102nd Annual Meeting of the American Association for Cancer Research; 2011 Apr 2-6; Orlando, FL. Philadelphia (PA): AACR; Cancer Res 2011;71(8 Suppl):Abstract nr 4838. doi:10.1158/1538-7445.AM2011-4838

Abstract Multiple Myeloma (MM) is an incurable plasma-cell disorder progressing from indolent monoclonal gammopathy of undetermined significance (MGUS) and smoldering myeloma (SMM) to relapse MM and finally plasma cell leukemia (PCL). Alterations affecting TP53 are among the principal genetic progression events associated with disease progression. We analyzed results from gene expression profiling (GEP), array-based comparative genomic hybridization (aCGH), fluorescent in situ hybridization (FISH) and mutational data on the largest cohort published to date of 806 patients from all MM stages and 48 human myeloma cell lines (HMCLs). The prevalence of monoallelic TP53 deletions by aCGH is 1.4%, 3.7%, 8.9% and 25%, in MGUS, SMM, newly diagnosed and relapsed MM respectively. It reaches 53% in HMCLs and 75% in secondary PCL. No biallelic TP53 deletion were observed in patients and only 8% in HMCLs. In relapsed MM, the prevalence of TP53 deletions increases from 12/75 patients (16%) at first, to 9/31 patients (29%) in second and 5/7 patients (71%) in third relapse or more. We studied by FISH the status of TP53 in 113 patients from early stages progressing to relapse. Overall 15/113 patients had TP53 deletion, we confirmed that 10/15 patients acquired TP53 deletion at late stages. Next, we investigated the mutational status of TP53 in 59 relapsed MM patients, 5/59 samples (8.5%) had TP53 mutations, 4/5 mutations were located in the DNA binding domain. Presence of TP53 deletion conferred poor overall survival in relapsed MM patients (4.2 vs. 37.8 months, p = 0.015). A cohort of 239 patients had GEP-aCGH paired analysis we used a validated 70-gene model to stratify patients based on calculations obtained by GEP. The cohort was divided into quartiles and we observed the prevalence of TP53 in the low and high-risk groups. TP53 deletion was observed in 33/239 patients (14%) by aCGH, 2/33 (6%) had GEP values in the low-risk and 13/33 (39%) in the high-risk group. Next, we used a previously proposed cut-off value predicting TP53 deletion when GEP values are less than 733. We tried to validate this value as a surrogate for direct detection of TP53 deletion, and used the positive predictive value (PPV) of such level as a potential clinically useful marker. The fact that aCGH is the gold standard approach to define copy-number alterations, the proposed cut-off value failed to predict TP53 deletion as the PPV was low at 18% and the sensitivity of 69%. Finally, we studied by aCGH the status of MDM2 and CDKN2A, two key regulators of TP53. MDM2 gain and CDKN2A loss were infrequent events in patients even in advanced disease. In 48 HMCLs, only one had MDM2 copy gain (2%) while 14 cell lines had CDKN2A loss (29%). In MM, emergence of TP53 deletion/mutations in MM and its increase in aggressive stages of disease remain the overriding genetic factor for poor prognosis determination and a marker of progressive genetic events. Based on GEP, the 733 cut-off value did not accurately predict TP53 deletion. Citation Format: {Authors}. {Abstract title} [abstract]. In: Proceedings of the 102nd Annual Meeting of the American Association for Cancer Research; 2011 Apr 2-6; Orlando, FL. Philadelphia (PA): AACR; Cancer Res 2011;71(8 Suppl):Abstract nr 219. doi:10.1158/1538-7445.AM2011-219

Abstract The longitudinal evolution of a myeloma genome from diagnosis to plasma cell leukemia has not previously been reported. We used whole-genome sequencing (WGS) on 4 purified tumor samples and patient germline DNA drawn over a 5-year period in a t(4;14) multiple myeloma patient. Tumor samples were acquired at diagnosis, first relapse, second relapse, and end-stage secondary plasma cell leukemia (sPCL). In addition to the t(4;14), all tumor time points also shared 10 common single-nucleotide variants (SNVs) on WGS comprising shared initiating events. Interestingly, we observed genomic sequence variants that waxed and waned with time in progressive tumors, suggesting the presence of multiple independent, yet related, clones at diagnosis that rose and fell in dominance. Five newly acquired SNVs, including truncating mutations of RB1 and ZKSCAN3, were observed only in the final sPCL sample suggesting leukemic transformation events. This longitudinal WGS characterization of the natural history of a high-risk myeloma patient demonstrated tumor heterogeneity at diagnosis with shifting dominance of tumor clones over time and has also identified potential mutations contributing to myelomagenesis as well as transformation from myeloma to overt extramedullary disease such as sPCL.

Abstract Emerging evidence indicates that tumors can follow several evolutionary paths over a patient's disease course. With the use of serial genomic analysis of samples collected at different points during the disease course of 28 patients with multiple myeloma, we found that the genomes of standard-risk patients show few changes over time, whereas those of cytogenetically high-risk patients show significantly more changes over time. The results indicate the existence of 3 temporal tumor types, which can either be genetically stable, linearly evolving, or heterogeneous clonal mixtures with shifting predominant clones. A detailed analysis of one high-risk patient sampled at 7 time points over the entire disease course identified 2 competing subclones that alternate in a back and forth manner for dominance with therapy until one clone underwent a dramatic linear evolution. With the use of the Vk*MYC genetically engineered mouse model of myeloma we modeled this competition between subclones for predominance occurring spontaneously and with therapeutic selection.

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 Microarray technologies are powerful assays, which can identify copy-number abnormalities (CNA) underlying the pathogenesis of cancer. Several commercial arrays are available but it is unclear which is the best platform. We did a side-by-side comparison between two leading arrays; the Affymetrix SNP6 and Agilent 1M, to compare data compression, noise, gene and exon coverage, and CNA detection. Five samples were run on both platforms and analyzed as recommended by the manufacturer. We first performed an in silico comparison to define the number of probes per serial 5k interval, number of probes per gene and exon, and the coefficient of variation of different smoothing windows. Comparing the CV of the two arrays showed a similar level of variation between a 3 and 5 probe averaging on the 1M and SNP6, respectively. Using these smoothing levels the 1M interrogates 81.2% of the genes in the genome compared to 69.3% on the SNP6. At the exon level the 1M has at least one probe in 64.8% compared to 41.2% on the SNP6. We observed a similar level of data compression (∼10%) on the two arrays but the SNP6 had nearly twice the noise (SD 0.38 versus 0.21). Finally, there was a marked increment in the number of CNA identified by 1M (mean of 78.5) than SNP6.0 (59.5). In this comparison the Agilent 1M out performed the Affymetrix SNP6. Even with half the number of probes, the 1M was superior in CNA detection, gene/exon coverage and showed less signal variation. A comparison with previous platforms showed that these two arrays did not identify new gene amplifications but did detect new deletions, however, neither platform detected several known deletions detected by other methods. Therefore, we designed a custom, exon-centered, 2×400K array (Agilent) with the goal of identifying additional deletions. The design outline was: a) all the probes from the Agilent CGH 44K array; b) 3 probes per exon across the human exome; c) high resolution tiling of 1450 cancer implicated genes; and d) tiling the 4kb flanking all MIRs. To improve performance, we only included probes with performance scores higher than 0.8 and 0.85 in queries b) and c), respectively. Lastly, 15000 probes from the 1M array that map to exons missed by our design parameters were also included. Our custom array has half the probes (415K) of the 1M (960K) and 4.5x less than the SNP6 (1.85M) but shows the best gene and exon coverage. When the minimum number of probes needed for calculating a CNA is considered, 86% of genes are interrogated compared to 81.2% (1M) and 69.3% (SNP6). The difference is even more marked at the exon level with 84% covered by at least one probe and at least 63.2% interrogated by two probes compared to 5.3% (1M) and 18.0% (SNP6). The promising in silico comparison of this exon-centered array, its more stringent design parameters in combination with significantly reducing the cost per run should make this a very promising array design for CNA analysis of tumor samples in the future 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 4957.