Yesil, Jennifer

BACKGROUND Aneuploidy, defined by abnormal copy number changes of chromosomes, contributes to genome instability in multiple myeloma and has potential prognostic impact. Previous research has explored transcriptional pathways affected by aneuploidy. We aim to evaluate aneuploidy, clinical prognosis, and gene expression in RRMM patients (pts) who participated in MMRF (Multiple Myeloma Research Foundation) sequencing study at University of Michigan. We further used gene set analysis to identify enrichment and variation of genetic pathways associated with aneuploidy and overall survival (OS). This was a pilot study in view of applying similar methods to larger populations. METHODS DNA and RNA materials were obtained from 51 RRMM pts at the time of disease relapse. Targeted sequencing was performed with the Onco1700 panel and RNA sequencing was performed with a capture protocol using Agilent SureSelect All Exon V4, followed by paired end sequencing. Copy number variation (CNV) was estimated using an in-house pipeline using matched normal samples. Arm-level aneuploidy was defined as copy number status (gain, loss or neutral) with maximal proportion for each chromosomal arm. Aneuploidy score was determined by total number of arm-level CNV aberrations, with the median score defining high and low aneuploidy groups. Survival analysis was performed using Kaplan-Meier (KM) and Cox regression. RNA-seq libraries were aligned with STAR aligner to the hg38 reference and read quantification were performed with featureCounts. We used gene set enrichment analysis (GSEA) and gene set variation analysis (GSVA) to obtain genetic pathway-level summaries for 50 hallmark pathways representative of cancer from the GSEA Molecular Signatures Database (https://www.gsea-msigdb.org/). Supervised GSEA was performed, based on differential analyses comparing gene expression profiles between high versus low aneuploidy groups to identify the differences in genetic pathways in RRMM pts. Unsupervised GSVA was applied to RNAseq count data of RRMM and newly-diagnosed multiple myeloma (NDMM) pts from the MMRF CoMMpass study (797 pts). Univariate Cox proportional hazard model was used to identify pathways having significant association with OS. RESULTS Arm scale aneuploidy analysis revealed high frequency of gains and losses in multiple chromosomes (Figure 1). RRMM pts with high aneuploidy scores had worse OS (median 15.9 months since study enrollment) compared to those with low scores (median not reached, p=0.027). Multivariate analysis demonstrated that high aneuploidy score is persistently associated with poor OS (hazard ratio (HR): 3.9, p=0.006), adjusting for t(4;14) and t(11;14). Hyperdiploidy (any whole chromosome carrying >2 copies) was identified in chromosomes 3, 5, 7, 9, 11, 19, 21, with frequency >20%. Pts with hyperdiploidy had better OS via KM (p=0.04). Using RNAseq gene expression data and GSEA analysis, we found enrichment of gene markers involving the pathways for targets of E2F transcription factors (Normalized Enrichment Score (NES) = 2.07), G2-M checkpoint in the cell division cycle (NES = 1.87), MYC proto-oncogene (NES = 1.93 and 1.85), and DNA repair (NES = 1.68) in RRMM pts with high aneuploidy compared against those with low scores (Benjamini-Hochberg adjusted p-value <0.001 for all pathways, Figure 2). GSVA analysis revealed 21 genetic pathways associated with survival in NDMM and RRMM pts, with 8, 2, and 11 pathways associated with OS in NDMM only, RRMM only, and both, respectively (Figure 3). Enrichment of several biological pathways involved with cellular growth, cell cycle regulation, and stress response (E2F, G2-M, MYC, DNA repair, mTORC1 signaling, UV radiation, mitotic spindle assembly, unfolded protein response, and peroxisome) are associated with poor survival (HR>1) in both NDMM and RRMM pts. Downregulation of KRAS signaling pathway is associated with better survival (HR<1) in both groups. CONCLUSION Aneuploidy, adjusted for chromosomal translocations, is an independent adverse prognostic factor for RRMM. Aneuploidy is correlated with dysregulation of major biological pathways involving cell cycle regulation, cellular division, oncogenes, and DNA damage repair. The majority of common pathways shared between RRMM and NDMM are linked to or reflective of genomic instability, and may explain association with poor survival. Disclosures Talpaz: IMAGO: Consultancy; Takeda: Research Funding; Novartis: Research Funding; Constellation Pharmaceuticals: Membership on an entity's Board of Directors or advisory committees; BMS: Membership on an entity's Board of Directors or advisory committees. Ye:Janssen: Research Funding; Portola: Research Funding; Millennium: Research Funding; Celgene: Research Funding; Sanofi: Research Funding; Karyopharm: Research Funding; Nektar: Research Funding; AbbVie: Research Funding.

Objective To identify the association between aneuploidy and clinical outcome in patients with relapsed and/or refractory multiple myeloma (RRMM) who participated in MMRF (Multiple Myeloma Research Foundation) sequencing study at University of Michigan. Background: Aneuploidy, defined by abnormal copy number changes of chromosomes, is one of the hallmarks in cancer, reflecting and also contributing to genome instability (Ye, Regan et al. 2018). Approximately 90% of cancers have gained or lost one or both arms of at least one chromosome (Taylor, Shih et al. 2018). In recent years, large scale sequencing efforts have extended aneuploidy study, identifying DNA somatic copy number alterations (SCNAs). Published data suggest aneuploidy has a stronger impact on prognosis than gene mutations in multiple myeloma (Walker, Boyle et al. 2015; Bolli, Biancon et al. 2018; Jamal-Haniani, Wilson et al. 2017). Methods: Fifty one RRMM patients from our institute participated in Clinical-Grade Molecular Profiling of Patients with Multiple Myeloma and Related Plasma Cell Malignancies (MMRF-002) from Multiple Myeloma Research Foundation (MMRF) between March 2016 and November 2018. Genomic DNA was obtained from CD138+ sorted myeloma cells and a peripheral blood sample from each patient. Capture exome sequencing on a targeted panel of 1500 genes was performed by the Illumina HiSeq 2500 (2x115bp paired-end reads, average 520X coverage). Copy number changes, loss of heterozygosity (LOH) and tumor purity were jointly estimated using an in-house pipeline for matched tumor/normal libraries. We assessed aneuploidy using chromosomal and arm level SCNAs which are determined by the median of SCNAs and summation of gain and loss of SCNAs for a given chromosomal arm, which are then used for a time-to-event analysis of overall survival (OS) of the patients. The differences between Kaplan-Meier overall survival curves were tested using the log-rank test. Hazard ratios (HR) were estimated from Cox proportional hazard regression. A threshold of significance was taken as p<0.05. Results Patient demographic data is summarized in Table 1. All the patients in this study harbor DNA somatic copy-number abnormalities (SCNAs) on the arm level for at least one chromosome (Figure 1). Gain of a chromosome or a whole arm mostly occurs on 1q, 3, 5, 7, 9, 11, 15, 19, 21 while loss was noted in 1p, 8p, 13, 16q, and 17p. At a global-level, we find that whole genomic aneuploidies (genomic-level across all chromosomes) have a significant association with inferior OS. Specifically, patients with a high frequency of genomic aneuploidies (defined as number of SCNA >8, the first quartile) had a significantly worse prognosis than those with low number of SCNA events (<=8) (p=0.037, Figure 2). Arm-level aneuploidies with high incidences include 11q (54.9%), 13q (52.9%), 15q (51.0%), 1q (51.0%), 9q (51.0%), a finding consistent with previous studies (Walker, Leone et al. 2010). Aneuploidies in 10p (11.8%) and 10q (7.8%) are uncommon in agreement with previous studies (Tricot G, Sawyer et al. 2017) as myeloma structural or numerical abnormalities rarely occur on chromosome 10. Using univariate analysis, we identified five arm-level aneuploidies which are significantly associated with worse OS: 10p, 10q, 11p, 18q and 20q. It is noteworthy that aneuploidies in either arm of chromosome 10 (10p and 10q) had significant negative OS prognostic effect: 10p (median survival for diploid 28.9 vs. aneuploid 4.1 mo; p<0.001) and 10q (28.9 vs 3.8 mo; p=0.015). Similar results with a negative impact on OS were found from 11p (not reached vs 13.8 mo; p=0.036); 18q (28.9 vs 8.5 mo; p=0.025); and 20q (28.9 vs 3.1 mo; p=0.009). Given the small sample size we did not conduct a multivariable analyses and multiplicity adjustment in this preliminary analyses. Conclusion Aneuploidy measured by SCNAs is correlated with unfavorable survival in relapsed/refractory myeloma. The higher the occurrence of aneuploidies, the worse the overall survival, illustrating of the impact of cancer genome instability. Disclosures Ye: Janssen: Research Funding; Karyopharm: Research Funding; Portola: Research Funding; MingSight: Research Funding; Sanofi: Research Funding. Talpaz:Imago BioSciences: Consultancy, Research Funding; Celgene: Consultancy, Research Funding; CTI BioPharma: Research Funding; Constellation: Research Funding; Incyte: Research Funding; Novartis: Research Funding; Samus Therapeutics: Research Funding. Bergsagel:Celgene: Consultancy; Ionis Pharmaceuticals: Consultancy; Janssen Pharmaceuticals: Consultancy.