Vo, Josh N.

Summary Cancer genome sequencing is essential for understanding tumor evolution and advancing precision medicine. 1 However, reference gaps and germline variants obscure detection of small and large somatic variants and methylation in repetitive regions. 1–3 It is common for tumor cells to gain or lose chromosome arms due to somatic structural changes that occur inside highly repetitive satellite DNA sequences in the centromeres. 4 To identify the full spectrum of somatic variants, including complex rearrangements, we construct and curate near-complete, haplotype-resolved assemblies of the most recent common ancestor of an early-passage broadly-consented hypodiploid pancreatic cancer cell line and matched normal tissues. The tumor assembly completely recapitulates all 35 tumor chromosomes observed with karyotyping, with multiple translocation-induced hybrid chromosomes. The hybrid chromosomes contain putative functional dicentric and fused centromeres, nested foldback inversions causing 14 breakpoints with a haplotype switch in a single event, and centromeric satellite tandem duplications up to 136 kbp. Direct comparison of tumor and normal assembly haplotypes uncovers >7,000 variants altering >1 Mbp of sequence in repetitive regions that have been hidden by reference gaps and germline variants. 44 % of somatic small variants change representation because they alter germline variants on GRCh38, impacting mutational signatures and kataegis/omikli clusters. Most somatic LINE insertions originate from two hypomethylated non-reference germline LINE insertions, highlighting their impact on insertion mutation burden. These assemblies demonstrate that centromeric, acrocentric, and telomeric regions conventionally excluded from analysis harbor extensive somatic and epigenetic changes. Resolving complete tumor genomes enables a deeper understanding of cancer structural plasticity and the endpoints of breakage-fusion-bridge cycles. These assembled, curated paired normal-tumor benchmarks will serve as a critical foundation for developing future algorithms to characterize the most intractable regions of cancer genomes.

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.