Role of Aneuploidy in Transcriptional Regulation and Clinical Prognosis in Relapsed and/or Refractory Multiple Myeloma (RRMM)

• Published in: Blood 2020-11, Vol.136 (Supplement 1), p.45-46
• Main Authors: Su, Christopher T, Chen, Liying, Chen, Jason, Parkin, Brian, Polk, Avery, Kandarpa, Malathi, Cole, Craig E., Campagnaro, Erica, Vo, Josh, Robinson, Dan, Wu, Yi-Mi, Talpaz, Moshe, Yesil, Jennifer, Auclair, Daniel, Bergsagel, P. Leif, Chinnaiyan, Arul, Baladandayuthapani, Veerabhadran, Ye, J Christine
• Format: Article
• Language: English
• ISSN: 0006-4971; 1528-0020
• DOI: 10.1182/blood-2020-134558

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 no