From DNA Copy Number Gains and Tumor Dependencies to Novel Therapeutic Targets for High-Risk Neuroblastoma

Neuroblastoma is a pediatric tumor arising from the sympatho-adrenal lineage and a worldwide leading cause of childhood cancer-related deaths. About half of high-risk patients die from the disease while survivors suffer from multiple therapy-related side-effects. While neuroblastomas present with a low mutational burden, focal and large segmental DNA copy number aberrations are highly recurrent and associated with poor survival. It can be assumed that the affected chromosomal regions contain critical genes implicated in neuroblastoma biology and behavior. More specifically, evidence has emerged that several of these genes are implicated in tumor dependencies thus potentially providing novel therapeutic entry points. In this review, we briefly review the current status of recurrent DNA copy number aberrations in neuroblastoma and provide an overview of the genes affected by these genomic variants for which a direct role in neuroblastoma has been established. Several of these genes are implicated in networks that positively regulate MYCN expression or stability as well as cell cycle control and apoptosis. Finally, we summarize alternative approaches to identify and prioritize candidate copy-number driven dependency genes for neuroblastoma offering novel therapeutic opportunities.

Genetic alterations associated with malignant transformation of sporadic vestibular schwannoma

Introduction Malignant peripheral nerve sheath tumor of the vestibulocochlear nerve (VN-MPNST) is exceedingly rare and carries a poor prognosis. Little is known about its underlying genetics and in particular the process of malignant transformation. There is an ongoing debate on whether the transformation is initiated by ionizing radiation. We present here the analysis and comparison of two post-radiation VN-MPNST and one undergoing spontaneous transformation. Methods Four tumors from three patients (radiation-naïve vestibular schwannoma before (VS) and after (VN-MPNST) malignant transformation in addition to two post-radiation VN-MPNST) were subjected to DNA whole-genome microarray and whole-exome sequencing and tumor-specific mutations were called. Mutational signatures were characterized using MuSiCa. Results The tumor genomes were characterized predominantly by copy-number aberrations with 36–81% of the genome affected. Even the VS genome was grossly aberrated. The spontaneous malignant transformation was characterized by a near-total whole-genome doubling, disappearance of NF2 mutation and new mutations in three cancer-related genes (GNAQ, FOXO4 and PDGFRB). All tumors had homozygous loss of the tumor suppressor CDKN2A. Neither mutational signature nor copy number profile was associated with ionizing radiation. Conclusion The VN-MPNST genome in our cases is characterized by large copy-number aberrations and homozygous deletion of CDKN2A. Our study demonstrates a VS with genetic alterations similar to its malignant counterpart, suggesting the existence of premalignant VS. No consistent mutational signature was associated with ionizing radiation.

The Genomic Landscape of Waldenström Macroglobulinemia Reveals Sustained Germinal Center Activity and Late-Developing Copy Number Aberrations

The genomic landscape of Waldenström Macroglobulinemia (WM) is characterized by recurrent somatic mutations in MYD88, with a lower incidence of mutations affecting CXCR4, ARID1A, CD79B and the NFKB signaling pathway (Hunter et. al. Blood 2014). We aimed to characterize the relationship between single base substitutions (SBS), mutational signatures, copy number aberrations (CNA) and structural variants (SV) in WM. We performed whole genome sequencing (WGS) on 14 primary samples from WM patients at various clinical stages, including IgM monoclonal gammopathy (n=1), smoldering (n=5), newly diagnosed (n=7) and relapsed WM (n=1). We identified a median of 2806 clonal SBS per sample (IQR 1870-3079), and 12/14 (85%) samples harbored MYD88 mutations. To investigate which mutational processes are involved in shaping the genomic landscape of WM we performed a mutational signature analysis. Four previously reported SBS signatures were detected: SBS1 and SBS5 (aging), SBS9 (germinal center; GC) and SBS8, with the contribution of age-related signatures SBS1/SBS5 being directly correlated with age at presentation (R 2=0.44, p=0.014). The GC signature SBS9 demonstrated sustained GC activity, as evidenced by the same proportion of mutations attributable to SBS9 at both the clonal and subclonal level (24%). At the immunoglobulin loci, we observed evidence of clustered SBS84 (AID), reflecting somatic hypermutation, with SBS84 accounting for 30% of signature contribution from subclonal mutations. Overall, these data suggest that, similarly to MM and other hematological malignancies, the interaction between WM and the GC is sustained over time. We have previously demonstrated that SV and complex events are critical in the pathogenesis and clinical outcomes of multiple myeloma. In contrast, in this WM WGS cohort, we found a low prevalence of complex SV, with no chromothripsis detected, and a single chromoplexy event found in 3 patients (21%), all of whom had progressed to symptomatic WM. To explore WM CNA features in a larger cohort, we examined the WGS data together with 38 MYD88-mutated WM samples for which targeted sequencing was available (MSK-IMPACT-Heme 400 gene panel). In this combined dataset (n=52), GISTIC analysis identified significantly deleted regions at 6q16.1, 7q34, 17p13.1 (TP53) and 21q22.2, along with significant amplification at 6p22.1 (HLA-A). To better characterize the HLA loci using the loss of heterozygosity in human leukocyte antigen (LOHHLA) tool (McGranahan et. al. Cell 2017) we found the presence of HLA-specific loss of heterozygosity in 1 sample, while 4 samples had HLA CN >2.5 (all from patients who progressed to symptomatic WM). CNA analysis demonstrated that while some samples harbored typical CNA features, others had minimal changes, with MGUS / smoldering WM samples having less CNA compared with those who progressed to symptomatic WM. The 2 MYD88 wild type WGS contained a clonal gain affecting chromosome 12, which is typically an early event in chronic lymphocytic leukemia. Molecular time analysis (the corrected ratio between duplicated and non-duplicated clonal mutations within large chromosomal gains [Maura et al. Nat Comm 2019]) demonstrated that these 2 chromosomal gain events occurred early in cancer development (relative timing <0.5), while multiple other CNA changes occurred later in the disease course (timing >0.5) and tended to be subclonal. This data suggests that, while MYD88-mutations are central to WM clone establishment and can be observed in precursor disease, CNA may contribute to later phases and disease progression. In summary, WGS in WM allows the demonstration that germinal center activity is sustained over time. CNA in WM are not random in distribution, with specific loci being significantly amplified or deleted, and a potential role for HLA CNA. In contrast to MYD88 mutations, which are carried by stable precursor patients, the subclonal status and late molecular time of most CNA changes suggest a late role in cancer progression. Disclosures Kastritis: Pfizer: Consultancy, Honoraria, Research Funding; Takeda: Honoraria; Janssen: Consultancy, Honoraria, Research Funding; Genesis Pharma: Honoraria; Amgen: Consultancy, Honoraria, Research Funding. Diamond: Sanofi: Honoraria; Medscape: Honoraria. Kazandjian: Arcellx: Honoraria, Membership on an entity's Board of Directors or advisory committees; BMS: Honoraria, Membership on an entity's Board of Directors or advisory committees. Papaemmanuil: Isabl Technologies: Divested equity in a private or publicly-traded company in the past 24 months; Kyowa Hakko Kirin Pharma: Consultancy. Dogan: Roche: Consultancy, Research Funding; Seattle Genetics: Consultancy; EUSA Pharma: Consultancy; Peer View: Honoraria; Takeda: Consultancy, Research Funding; Physicians' Education Resource: Honoraria. Lesokhin: Trillium Therapeutics: Consultancy; Serametrix, Inc: Patents & Royalties; Genetech: Research Funding; Iteos: Consultancy; bristol myers squibb: Research Funding; Behringer Ingelheim: Honoraria; pfizer: Consultancy, Research Funding; Janssen: Honoraria, Research Funding. Landgren: Janssen: Other: IDMC; Janssen: Honoraria; Celgene: Research Funding; Amgen: Honoraria; Janssen: Research Funding; Amgen: Research Funding; Takeda: Other: IDMC; GSK: Honoraria. Palomba: Rheos: Honoraria; Pluto: Honoraria; Lygenesis: Honoraria; Ceramedix: Honoraria; Seres: Honoraria, Other: Stock, Patents & Royalties, Research Funding; Nektar: Honoraria; PCYC: Consultancy; Wolters Kluwer: Patents & Royalties; Notch: Honoraria, Other: Stock; Priothera: Honoraria; Kite: Consultancy; Novartis: Consultancy; Magenta: Honoraria; WindMIL: Honoraria; BeiGene: Consultancy; Juno: Patents & Royalties. Maura: OncLive: Honoraria; Medscape: Consultancy, Honoraria. Dimopoulos: Amgen: Honoraria; BMS: Honoraria; Janssen: Honoraria; Takeda: Honoraria; BeiGene: Honoraria.

Whole genome copy number analyses reveal a highly aberrant genome in TP53 mutant lung adenocarcinoma tumors

Background Genetic alterations are common in non-small cell lung cancer (NSCLC), and DNA mutations and translocations are targets for therapy. Copy number aberrations occur frequently in NSCLC tumors and may influence gene expression and further alter signaling pathways. In this study we aimed to characterize the genomic architecture of NSCLC tumors and to identify genomic differences between tumors stratified by histology and mutation status. Furthermore, we sought to integrate DNA copy number data with mRNA expression to find genes with expression putatively regulated by copy number aberrations and the oncogenic pathways associated with these affected genes. Methods Copy number data were obtained from 190 resected early-stage NSCLC tumors and gene expression data were available from 113 of the adenocarcinomas. Clinical and histopathological data were known, and EGFR-, KRAS- and TP53 mutation status was determined. Allele-specific copy number profiles were calculated using ASCAT, and regional copy number aberration were subsequently obtained and analyzed jointly with the gene expression data. Results The NSCLC tumors tissue displayed overall complex DNA copy number profiles with numerous recurrent aberrations. Despite histological differences, tissue samples from squamous cell carcinomas and adenocarcinomas had remarkably similar copy number patterns. The TP53-mutated lung adenocarcinomas displayed a highly aberrant genome, with significantly altered copy number profiles including gains, losses and focal complex events. The EGFR-mutant lung adenocarcinomas had specific arm-wise aberrations particularly at chromosome7p and 9q. A large number of genes displayed correlation between copy number and expression level, and the PI(3)K-mTOR pathway was highly enriched for such genes. Conclusions The genomic architecture in NSCLC tumors is complex, and particularly TP53-mutated lung adenocarcinomas displayed highly aberrant copy number profiles. We suggest to always include TP53-mutation status when studying copy number aberrations in NSCLC tumors. Copy number may further impact gene expression and alter cellular signaling pathways.

Sparse canonical correlation to identify breast cancer related genes regulated by copy number aberrations

Background: Copy number aberrations (CNA) have proved to be of clinical and therapeutic significance for many diseases including breast cancer, since they drive numerous key underlying biological processes, by regulating molecular phenotypes like gene expression and others. To comprehensively assess the effect of CNAs, it is not sufficient to only identify significant CNA-gene expression pairs, but also to identify the overall gene networks and regulatory structures that are influenced by CNAs, subsequently producing change in outcomes. Methods: In this article, we adopt a two-step analysis approach to identify CNA regulated genes whose expression levels affect breast cancer related outcomes: (1) we identify gene modules that are regulated by CNAs through sparse canonical correlation analysis (sCCA) which selects a set of closely located CNAs that regulates the expression levels of selected genes. (2) then, we use a using generalized linear model, to identify which genes within the gene modules are associated with breast cancer related outcomes. Results: Analyzing clinical and genomic data on 1904 breast cancer patients from the METABRIC study, we found 14 gene modules to be regulated by groups of proximally located CNA sites. The identification of gene modules was further validated using independent data on individuals in a study of breast invasive carcinoma from The Cancer Genome Atlas (TCGA). Association analysis on 7 different breast cancer related outcomes identified several novel and interpretable regulatory associations which highlights how CNA can impact key biological pathways and process in context of breast cancer. Through downstream analysis of two example outcomes: estrogen receptor status and overall survival, we show that the identified genes were enriched in relevant biological pathways and the key advantage of our method is that we additionally identify the CNA that regulate these genes. Due to the availability of multiple types of outcomes, we further meta-analyzed the results to identify genes that had potentially associations with multiple outcomes. Conclusions: Overall we present a generalizable analysis approach to identify genes associated to different outcomes that are regulated by sets of CNA and can further be used to combine results across various types of outcomes. The results show that our method can identify novel and interpretable associations, by providing mechanistic insights on how the effects of CNA are cascaded via gene expression to impact breast cancer and related outcomes.

Whole Genome Copy Number Analyses Reveal a Highly Aberrant Genome in TP53 Mutant Lung Adenocarcinoma Tumors

Background: Genetic alterations are common in non-small cell lung cancer (NSCLC), and DNA mutations and translocations are targets for therapy. Copy number aberrations occur frequently in NSCLC tumors and may influence gene expression and further alter signaling pathways. In this study we aimed to characterize the genomic architecture of NSCLC tumors and to identify genomic differences between tumors stratified by histology and mutation status. Furthermore, we sought to integrate DNA copy number data with mRNA expression to find genes with expression putatively regulated by copy number aberrations and the oncogenic pathways associated with these affected genes. Methods: Copy number data were obtained from 190 resected early-stage NSCLC tumors and gene expression data were available from 113 of the adenocarcinomas. Clinical and histopathological data were known, and EGFR-, KRAS- and TP53 mutation status was determined. Allele-specific copy number profiles were calculated using ASCAT, and regional copy number aberration were subsequently obtained and analyzed jointly with the gene expression data.Results:The NSCLC tumors tissue displayed overall complex DNA copy number profiles with numerous recurrent aberrations. Despite histological differences, tissue samples from squamous cell carcinomas and adenocarcinomas had remarkably similar copy number patterns. The TP53-mutated lung adenocarcinomas displayed a highly aberrant genome, with significantly altered copy number profiles including gains, losses and focal complex events. The EGFR-mutant lung adenocarcinomas had specific arm-wise aberrations particularly at chromosome7p and 9q. A large number of genes displayed correlation between copy number and expression level, and the PI(3)K-mTOR pathway was highly enriched for such genes.Conclusions: The genomic architecture in NSCLC tumors is complex, and particularly TP53-mutated lung adenocarcinomas displayed highly aberrant copy number profiles. We suggest to always include TP53-mutation status when studying copy number aberrations in NSCLC tumors. Copy number may further impact gene expression and alter cellular signaling pathways.

Signatures of Discriminative Copy Number Aberrations in 31 Cancer Subtypes

Copy number aberrations (CNA) are one of the most important classes of genomic mutations related to oncogenetic effects. In the past three decades, a vast amount of CNA data has been generated by molecular-cytogenetic and genome sequencing based methods. While this data has been instrumental in the identification of cancer-related genes and promoted research into the relation between CNA and histo-pathologically defined cancer types, the heterogeneity of source data and derived CNV profiles pose great challenges for data integration and comparative analysis. Furthermore, a majority of existing studies have been focused on the association of CNA to pre-selected “driver” genes with limited application to rare drivers and other genomic elements. In this study, we developed a bioinformatics pipeline to integrate a collection of 44,988 high-quality CNA profiles of high diversity. Using a hybrid model of neural networks and attention algorithm, we generated the CNA signatures of 31 cancer subtypes, depicting the uniqueness of their respective CNA landscapes. Finally, we constructed a multi-label classifier to identify the cancer type and the organ of origin from copy number profiling data. The investigation of the signatures suggested common patterns, not only of physiologically related cancer types but also of clinico-pathologically distant cancer types such as different cancers originating from the neural crest. Further experiments of classification models confirmed the effectiveness of the signatures in distinguishing different cancer types and demonstrated their potential in tumor classification.