Recurrent Copy Number Aberrations Research Articles

Genomic, Transcriptional, and Immunological Validation of Distinct Molecular Subtypes of Classic Hodgkin Lymphoma through Tissue-Based and Noninvasive Methods

Introduction: The scarcity of malignant Hodgkin and Reed-Sternberg (HRS) cells has hampered comprehensive genomic profiling of classic Hodgkin lymphoma (cHL) from tumor tissue. Multiple recent studies have demonstrated that plasma cfDNA profiling facilitates molecular characterization of cHL. Leveraging noninvasive genotypes and Latent Dirichlet Allocation, we recently defined 2 genetic cHL subtypes in a large international cohort comprising 366 patients including pediatric and adult patients of all ages (Alig et al., ASH 2022 and ICML 2023). Cluster H1 comprised ~2/3 of cases and was dominated by a high somatic mutational burden, and non-silent mutations in genes canonically involved in NF-κB, JAK/STAT and PI3K signaling as well as in B2M. Conversely, cluster H2 (~1/3 of cases) was primarily characterized by recurrent somatic copy number aberrations as well as mutations in TP53 and KMT2D. Herein, we validate these previously identified genetic subtypes in external datasets as well as through orthogonal methods including tissue-based and noninvasive transcriptional and immune profiling. Methods: To validate genetic cHL subtypes, we first leveraged public cHL genotypes from 61 patients obtained through whole genome/exome sequencing (WGS/WES) of flow-sorted HRS cells (Maura et al, Blood Cancer Disc. 2023). Using our previously locked down probabilistic classifier, we assigned the H1/H2 subtype, and then correlated cluster assignments with clinical variables. To explore transcriptional differences between genetic subtypes, we profiled baseline plasma samples (n=113) from the larger plasma genotyping cohort (n=293) using EPIC-Seq (Esfahani et al, Nat Biotechnol 2022), which allows for noninvasive gene expression profiling from cfDNA fragmentation patterns at transcription start sites. Further, we used SABER (Sworder, Cancer Cell 2023) to enumerate T-cell receptor (TCR) rearrangements in cfDNA (cfTCR, n=292). Finally, to assess subtype-specific immune infiltrate patterns, we profiled cHL tumor specimens using RNA-Seq (n=64), and applied CIBERSORTx. Results: Similar to observations in our cfDNA discovery cohort, H1 was found to be the more prevalent subtype in the external validation set comprising 56% of tumors, while 44% were classified as H2. Recurrence frequencies of genetic features were comparable to and significantly correlated with those from our plasma discovery cohort ( R S=0.59 [H1] and R S=0.63 [H2], P<0.001 each). Of note, when considering the whole genome space, the higher mutational burden of H1 tumors could be confirmed ( P<0.01), and this association was independent of the tumor EBV status. In agreement with the discovery cohort, the bimodal age distribution and increased EBV positivity of the H2 subtype could also be recapitulated (31% vs 6% EBV+, P<0.01). To explore transcriptional differences between H1/H2 subtypes, we took advantage of the plasma enrichment of HRS cell ctDNA and utilized EPIC-Seq to noninvasively infer expression of 1,676 targeted genes. Tracking signatures derived from scRNA-Sequencing in plasma samples, we found that both HRS cells and the cHL tumor microenvironment can be successfully profiled by EPIC-Seq. Strikingly, we found substantial enrichment of a cytokine response signature in H1 tumors, while T-cell activation was among the top upregulated signatures in H2 tumors ( Fig. A). Importantly, the T-cell signature found in H2 was accompanied by a higher abundance of T-cell clones as quantified by SABER in baseline plasma samples ( P<0.001, Fig. B). Notably, cfTCR fragment length profiles resembled the mutant ctDNA profiles, strongly suggesting a tumor origin of the TCR rearrangements detected in plasma. Lastly, immune cell deconvolution of bulk RNA-Seq specimens indicated a higher abundance of CD8+ T-cells in H2 tumors ( P<0.01), further confirming our prior observations. Conclusions: Collectively, these results serve to validate H1 and H2 as distinct cHL subtypes, to confirm the characteristic genotypes defining them, and recapitulate their distinctive associations with key clinical and pathological variables including age and EBV status. We further validate the subtypes using orthogonal methods revealing dominant cytokine driven signaling in H1. Conversely, H2 tumors, which largely lack B2M mutations, despite their lower mutational burden, are rather immunogenic triggering T-cell responses.

Abstract 6067: Characterizing the landscape of rare structural variant events in newly diagnosed multiple myeloma

Abstract Background: Structural variants (SV) are known to play a critical role in the pathogenesis of multiple cancer types. Using whole genome sequencing (WGS), we recently characterized the SV landscape of 752 newly diagnosed multiple myeloma (MM) patients, identifying 68 SV hotspots and 152 recurrent copy number aberrations (CNA; Rustad et al. Blood Cancer Discovery 2020). Despite comprehensive annotation, more than half of SVs were not linked to any known MM genomic driver. The biological impact of these SV events, here defined as rare SV, occurring in 93% (702/752) of patients, is unknown. Methods: To study the biological impact of rare SVs, we interrogated WGS (n=752) and RNAseq (n=591) in the CoMMpass trial. Recurrent SVs identified by involvement in canonical Ig translocations, recurrent CNAs, or SV hotspots were excluded. All SVs within an event must not involve a recurrent region to be defined rare. To determine SV class-specific gene relationships, breakpoint enrichment was compared against a permuted background model for each SV class and gene expression direction, up to 1 Mb. Genes were considered affected if expression was above a gene specific outlier Z-score of +/- 2. Lastly, we modeled breakpoint density to the nearest MM superenhancer up to 10 MB, and compared to permuted background rates. Results: Of the total 8,942 SVs, 4,959 (55%) were identified as rare. 201 (34%) patients had at least 1 enriched rare SV event associated with gene expression outliers. Amongst over-expressed gene outliers, rare templated insertions and duplications were enriched within the gene body and up to 1 MB away. Rare inversions were enriched in genes 100kb and 1MB away, and rare translocations were associated with outliers 1 MB away. Amongst under-expressed gene outliers, rare complex SVs were enriched within the gene body, while deletion SVs were enriched in the gene body and up to 1 MB away. Rare duplications, translocations and templated insertions were enriched up to 1 MB of superenhancers. Rare templated insertions were significantly enriched against the background model (p < 0.001, Fisher Exact). Overall, 82% (104/126) of gene outliers affected by rare templated insertions were associated with superenhancers, (95 over-expressed, e.g. IRF6 and 9 under-expressed), 54% (130/237) by rare translocations, (105 overexpressed, e.g. FAM46A and 25 under-expressed), and 55% (96/172) by rare duplication events, (93 overexpressed, e.g. MAPK13, and 7 under-expressed). In addition, among the 853 outlier genes affected by enriched rare SVs, at least 15 are involved with B-cell development, suggesting a potential driver role in myeloma pathogenesis. Conclusion: In summary, leveraging WGS and RNA-seq of clinical samples, we demonstrate that rare SVs are frequently associated with aberrant gene expression, expanding our understanding of their potential role in heterogenous clinical response in patients diagnosed with MM. Citation Format: Monika Chojnacka, Benjamin Diamond, Bachisio Ziccheddu, Even Rustad, Kylee Maclachlan, Marios Papadimitriou, Eileen Boyle, Patrick Blaney, Saad Usmani, Gareth Morgan, Ola Landgren, Francesco Maura. Characterizing the landscape of rare structural variant events in newly diagnosed multiple myeloma [abstract]. In: Proceedings of the American Association for Cancer Research Annual Meeting 2023; Part 1 (Regular and Invited Abstracts); 2023 Apr 14-19; Orlando, FL. Philadelphia (PA): AACR; Cancer Res 2023;83(7_Suppl):Abstract nr 6067.

Improved Uveal Melanoma Copy Number Subtypes Including an Ultra-High-Risk Group.

PurposeTo evaluate the clinical relevance of low-frequency copy number aberrations (CNAs) in uveal melanoma (UM) and to discern residual genomic and clinical heterogeneity within established molecular subtypes based on genome-wide CNA profiling of 921 primary tumors.DesignRetrospective single-center case series.ParticipantsPatients with primary UM referred for genetic testing between 2008 and 2016 (n = 921). The Cancer Genome Atlas cohort with clinical outcome data available (n = 70) was used to validate findings.MethodsGenome-wide CNAs were generated for primary tumors from 921 patients and for 19 metastatic UM (mUM) in the liver. Of the 921 patients, metastatic outcome was known for 678 patients with a median time to metastasis of 4.5 years. The primary tumors were processed on the Affymetrix arrays SNP-5.0 (n = 140), SNP-6.0 (n = 359), or CytoScanHD (n = 422), and the metastatic tumors on the CytoScanHD array (n = 19). Recurrent CNAs were identified, and the prognostic effect of individual CNAs and multiple CNA clustering strategies, including more specific molecular subgroups with rare CNAs, were evaluated.Main Outcome MeasuresCNA recurrence, and effect of CNAs and derived molecular subtypes on metastatic-free survival.ResultsGenomic profiling revealed CNAs associated with risk of metastasis and demonstrated a strong association between chromosomal instability and patient prognosis. Using standard prognostic CNAs, 6 clusters were detected, and inclusion of chromosome 16q deletion revealed an additional cluster. Of these 7 genomic clusters, 5 patient groups showed distinct rates of metastasis, indicating that different genomic patterns can have similar patient outcomes. A small group of patients with a significantly higher rate of metastasis was characterized by monosomy 3, 8q amplification, and deletion of 1p or 16q. Although this ultra–high-risk group accounts for only 7% of this cohort, 88% demonstrated metastasis within 4 years, compared with 45% in the second-highest risk group.ConclusionsThese results suggest that 1p and 16q deletion should be incorporated in clinical assays to assess prognosis at diagnosis and to guide enrollment in clinical trials for adjuvant therapies.

Low Mutational Burden of Extra Nodal Marginal Zone Lymphoma of Mucosa-Associated Lymphoid Tissue in Patients with Primary Sjogren's Syndrome

Background: Primary Sjogren's syndrome (pSS) is an autoimmune disease characterized by chronic inflammation of the exocrine glands. Patients with pSS are at increased risk of developing extra nodal marginal zone lymphoma of the mucosa-associated lymphoid tissue (MALT) type, predominantly in the salivary glands. The presence of cryoglobulinemia, low C4 levels and salivary gland enlargement, are important risk factors for lymphoma development in pSS. Although lymphoma associated mutations have been described in circulating B-cells of patients with pSS associated cryoglobulinemia, it is unknown which somatic aberrations underlie the development of pSS-associated MALT lymphomas. The aim of the current study was to define the genomic landscape of pSS salivary gland MALT lymphomas. Identification of specific recurrent mutations or copy number aberrations and defining the affected pathways would help to understand the mechanisms by which intra-epithelial B-cells undergo malignant transformation.Methods: Whole exome sequencing was performed on 14 fresh frozen and 3 paraffin embedded MALT tissue samples of pSS patients at the University Medical Center Groningen (UMCG). Matched peripheral white blood cell samples were available for 12 patients and were used as germline controls. Fluorescence in situ hybridization was performed for the detection of MALT1 translocations.Results: Presence of a clonal B-cell population, indicative of a MALT lymphoma was confirmed by IgH PCR (BIOMED-2) for all patients. Whole exome sequencing resulted in a median target coverage of 130x, ranging from 84-163. More than 90% of the target regions had a coverage of >50x. In total we identified 232 somatic mutations in 184 genes. Three cases had a relatively high mutational load (>25 / case), while the median number of mutations in the remaining 14 cases was 7. Across the 17 cases there were 18 recurrently mutated genes. PRKD1, associated with malignant epithelial tumor of the salivary glands, was the most frequently mutated gene (six cases). Besides PRKD1, five additional recurrently mutated genes are also involved in epithelial surface and/or extracellular matrix (MAMDC4, COL14A1, CAMSAP3, TMEM2 and MUC4). Five genes, all with mutations in two cases, were shown to be associated with lymphomagenesis (ID3, TBL1XR1, PAX5, IGLL5, APC). A total of 18 copy number alterations (CNAs) were detected in 8 of the 14 evaluable cases, with gains of part of chromosomes 2 and 19, and loss of chromosome 6 each being detected in two cases. With respect to outcome, only 2 cases with high mutational load relapsed outside of the salivary glands, suggesting a more advanced stage of lymphoma.Conclusion: The low mutational load and lack of a clear lymphoma related gene signature suggests that localized pSS MALT lymphomas are genetically stable and most likely depend on the inflammatory state of the micro-environment. Only those cases characterized by higher mutational load showed a relapse-remitting disease, as is typical for indolent lymphoma. These observations could be translated to the clinical setting and potentially have value as biomarker for identification of patients with a more aggressive disease. Confirmation in larger cohorts are required to confirm this potential association. This study also paves way for treatment strategies targeting the micro-environment. DisclosuresNo relevant conflicts of interest to declare.

Copy Number Aberration Analysis to Predict Response to Neoadjuvant Anti-HER2 Therapy: Results from the NeoALTTO Phase III Clinical Trial.

The heterogeneity of response to anti-HER2 agents represents a major challenge in patients with HER2-positive breast cancer. To better understand the sensitivity and resistance to trastuzumab and lapatinib, we investigated the role of copy number aberrations (CNA) in predicting pathologic complete response (pCR) and survival outcomes in the NeoALTTO trial. The neoadjuvant phase III NeoALTTO trial enrolled 455 patients with HER2-positive early-stage breast cancer. DNA samples from 269 patients were assessed for genome-wide copy number profiling. Recurrent CNAs were found with GISTIC2.0. CNA estimates were obtained for 184 patients included in NeoALTTO. Among those, matched transcriptome and whole-exome data were available for 154 and 181 patients, respectively. A significant association between gene copy number and pCR was demonstrated for ERBB2 amplification. Nevertheless, ERBB2 amplification ceased to be predictive once ERBB2 expression level was considered. GISTIC2.0 analysis revealed 159 recurrent CNA regions. Lower copy number levels of the 6q23-24 locus predicted absence of pCR in the whole cohort and in the estrogen receptor-positive subgroup. 6q23-24 deletion was significantly more frequent in TP53 wild-type (WT) compared with TP53-mutated, resulting in copy number levels significantly associated with lack of pCR only in the TP53 WT subgroup. Interestingly, a gene-ontology analysis highlighted several immune processes correlated to 6q23-24 copy number. Our analysis identified ERBB2 copy number as well as 6q23-24 CNAs as predictors of response to anti-HER2-based treatment. ERBB2 expression outperformed ERBB2 amplification. The complexity of the 6q23-24 region warrants further investigation.

Detection of new drivers of frequent B-cell lymphoid neoplasms using an integrated analysis of whole genomes.

B-cell lymphoproliferative disorders exhibit a diverse spectrum of diagnostic entities with heterogeneous behaviour. Multiple efforts have focused on the determination of the genomic drivers of B-cell lymphoma subtypes. In the meantime, the aggregation of diverse tumors in pan-cancer genomic studies has become a useful tool to detect new driver genes, while enabling the comparison of mutational patterns across tumors. Here we present an integrated analysis of 354 B-cell lymphoid disorders. 112 recurrently mutated genes were discovered, of which KMT2D, CREBBP, IGLL5 and BCL2 were the most frequent, and 31 genes were putative new drivers. Mutations in CREBBP, TNFRSF14 and KMT2D predominated in follicular lymphoma, whereas those in BTG2, HTA-A and PIM1 were more frequent in diffuse large B-cell lymphoma. Additionally, we discovered 31 significantly mutated protein networks, reinforcing the role of genes such as CREBBP, EEF1A1, STAT6, GNA13 and TP53, but also pointing towards a myriad of infrequent players in lymphomagenesis. Finally, we report aberrant expression of oncogenes and tumor suppressors associated with novel noncoding mutations (DTX1 and S1PR2), and new recurrent copy number aberrations affecting immune check-point regulators (CD83, PVR) and B-cell specific genes (TNFRSF13C). Our analysis expands the number of mutational drivers of B-cell lymphoid neoplasms, and identifies several differential somatic events between disease subtypes.

The landscape of copy number variations in classical Hodgkin lymphoma: a joint KU Leuven and LYSA study on cell-free DNA

The low abundance of Hodgkin/Reed-Sternberg (HRS) cells in lymph node biopsies in classical Hodgkin lymphoma (cHL) complicates the analysis of somatic genetic alterations in HRS cells. As circulating cell-free DNA (cfDNA) contains circulating tumor DNA (ctDNA) from HRS cells, we prospectively collected cfDNA from 177 patients with newly diagnosed, mostly early-stage cHL in a monocentric study at Leuven, Belgium (n = 59) and the multicentric BREACH study by Lymphoma Study Association (n = 118). To catalog the patterns and frequencies of genomic copy number aberrations (CNAs), cfDNA was sequenced at low coverage (0.26×), and data were analyzed with ichorCNA to yield read depth-based copy number profiles and estimated clonal fractions in cfDNA. At diagnosis, the cfDNA concentration, estimated clonal fraction, and ctDNA concentration were significantly higher in cHL cases than controls. More than 90% of patients exhibited CNAs in cfDNA. The most frequent gains encompassed 2p16 (69%), 5p14 (50%), 12q13 (50%), 9p24 (50%), 5q (44%), 17q (43%), 2q (41%). Losses mostly affected 13q (57%), 6q25-q27 (55%), 4q35 (50%), 11q23 (44%), 8p21 (43%). In addition, we identified loss of 3p13-p26 and of 12q21-q24 and gain of 15q21-q26 as novel recurrent CNAs in cHL. At diagnosis, ctDNA concentration was associated with advanced disease, male sex, extensive nodal disease, elevated erythrocyte sedimentation rate, metabolic tumor volume, and HRS cell burden. CNAs and ctDNA rapidly diminished upon treatment initiation, and persistence of CNAs was associated with increased probability of relapse. This study endorses the development of ctDNA as gateway to the HRS genome and substrate for early disease response evaluation.

Distinct genomic features across cytolytic subgroups in skin melanoma

BackgroundSkin melanoma is a highly immunogenic cancer. The intratumoral immune cytolytic activity (CYT) reflects the ability of cytotoxic T and NK cells to eliminate cancer cells, and is associated with improved patient survival. Despite the enthusiastic clinical results seen in advanced-stage metastatic melanoma patients treated with immune checkpoint inhibitors, a subgroup of them will later relapse and develop acquired resistance. We questioned whether CYT associates with different genomic profiles and thus, patient outcome, in skin melanoma.MethodsWe explored the TCGA-SKCM dataset and stratified patients to distinct subgroups of cytolytic activity. The tumor immune contexture, somatic mutations and recurrent copy number aberrations were calculated using quanTIseq, MutSigCV and GISTIC2. Chromothriptic events were explored using CTLPScanner and cancer neoepitopes were predicted with antigen garnish. Each tumor's immunophenoscore was calculated using Immunophenogram. Mutational signatures and kataegis were explored using SigProfiler and compared to the known single or doublet base substitution signatures from COSMIC.ResultsMetastatic skin melanomas had significantly higher CYT levels compared to primary tumors. We assessed enrichment for immune-related gene sets within CYT-high tumors, whereas, CYT-low tumors were enriched for non-immune related gene sets. In addition, distinct mutational and neoantigen loads, primarily composed of C > T transitions, along with specific types of copy number aberrations, characterized each cytolytic subgroup. We found a broader pattern of chromothripsis across CYT-low tumors, where chromosomal regions harboring chromothriptic events, contained a higher number of cancer genes. SBS7a/b, SBS5 and SBS1 were the most prevalent mutational signatures across both cytolytic subgroups, but SBS1 differed significantly between them. SBS7a/b was mutually exclusive with SBS5 and SBS1 in both CYT subgroups. CYT-high patients had markedly higher immunophenoscore, suggesting that they should display a clinical benefit upon treatment with immune checkpoint inhibition therapy, compared to CYT-low patients.ConclusionsOverall, our data highlight the existence of distinct genomic features across cytolytic subgroups in skin melanoma, which might affect the patients' relapse rate or their acquisition of resistance to immune checkpoint inhibition therapies.

Somatic Copy Number Alterations and Associated Genes in Clear-Cell Renal-Cell Carcinoma in Brazilian Patients.

Somatic copy number aberrations (CNAs) have been associated with clear-cell renal carcinoma (ccRCC) pathogenesis and are a potential source of new diagnostic, prognostic and therapeutic biomarkers. Recurrent CNAs include loss of chromosome arms 3p, 14q, 9p, and gains of 5q and 8q. Some of these regional CNAs are suspected of altering gene expression and could influence clinical outcomes. Despite many studies of CNAs in RCC, there are currently no descriptions of genomic copy number alterations in a Brazilian ccRCC cohort. This study was designed to evaluate the chromosomal profile of CNAs in Brazilian ccRCC tumors and explore clinical associations. A total of 92 ccRCC Brazilian patients that underwent nephrectomy at Barretos Cancer Hospital were analyzed for CNAs by array comparative genomic hybridization. Most patients in the cohort had early-stage localized disease. The most significant alterations were loss of 3p (87.3%), 14q (35.8%), 6q (29.3%), 9p (28.6%) and 10q (25.0%), and gains of 5q (59.7%), 7p (29.3%) and 16q (20.6%). Bioinformatics analysis revealed 19 genes mapping to CNA significant regions, including SETD2, BAP1, FLT4, PTEN, FGFR4 and NSD1. Moreover, gain of 5q34-q35.3 (FLT4 and NSD1) and loss of 6q23.2-q23.3 (MYB) and 9p21.3 (MLLT3) had gene expression levels that correlated with TCGA data and was also associated with advanced disease features, such as larger tumors, Fuhrman 3, metastasis at diagnosis and death. The loss of region 14q22.1 which encompasses the NIN gene was associated with poor overall survival. Overall, this study provides the first CNA landscape of Brazilian patients and pinpoints genomic regions and specific genes worthy of more detailed investigations. Our results highlight important genes that are associated with copy number changes involving large chromosomal regions that are potentially related to ccRCC tumorigenesis and disease biology for future clinical investigations.

TAMI-59. METABOLIC REGULATION OF THE EPIGENOME DRIVES LETHAL INFANTILE EPENDYMOMA

Abstract Ependymomas are malignant glial tumors that occur throughout the central nervous system. Of the nine distinct molecular types of ependymoma, Posterior Fossa A(PFA) ependymomas, found in the hindbrain of infants and young children, are the most prevalent type. They rarely harbor recurrent somatic single nucleotide or copy number aberrations, but rather display epigenetic dysregulations. Transcriptional analyses show that PFAs, but not other molecular ependymoma variants, have elevated hypoxic signaling. We hypothesized that the metabolic environment of the developing human fetal hindbrain contributes to PFA ependymomas through the intermediary metabolic mechanism. Hypoxia microenvironment is essential for PFA survival and upholding their epigenetic dysregulation. Hypoxia blocks methyltransferase activity by upregulating EZHIPs expression and restricting SAM abundance. Fine-tuning the abundance of a-KG and acetyl-CoA, hypoxia fuels demethylase, and acetyltransferase activity which collectively resulting in H3K27 hypomethylation and hyperacetylation. Genome-wide essentially screen further underscore that minimal basal level of H3K27me3 is essential for PFA survival and further attenuation of H3K27me3 decreases the fitness of PFAs. PFA Ependymomas have a unique epigenome, suggesting a model in which they thrive in a narrow Goldilocks zone, with deviation to either increased or decreased H3K27me3 levels leading to diminished cellular fitness. PFAs are maintained under hypoxia, associated with restricted availability of specific metabolites to diminish histone methylation and increase histone demethylation and acetylation at H3K27. PFAs initiate from a cell lineage in the first trimester of human development that resides in restricted oxygen. Unlike other ependymomas, transient exposure of PFA cells to ambient oxygen induces irreversible cellular toxicity. PFA tumors exhibit a low basal level of H3K27me3 and, paradoxically, inhibition of H3K27 methylation specifically disrupts PFA growth. Microenvironmental regulation of PFA ependymoma epigenome appears to play a major role in tumorigenesis. Targeting metabolism and/or the epigenome presents a unique opportunity for rational therapy for infants with PFA ependymoma.

Targeting TRIM37-driven centrosome dysfunction in 17q23-amplified breast cancer

Genomic instability is a hallmark of cancer, and has a central role in the initiation and development of breast cancer1,2. The success of poly-ADP ribose polymerase inhibitors in the treatment of breast cancers that are deficient in homologous recombination exemplifies the utility of synthetically lethal genetic interactions in the treatment of breast cancers that are driven by genomic instability3. Given that defects in homologous recombination are present in only a subset of breast cancers, there is a need to identify additional driver mechanisms for genomic instability and targeted strategies to exploit these defects in the treatment of cancer. Here we show that centrosome depletion induces synthetic lethality in cancer cells that contain the 17q23 amplicon, a recurrent copy number aberration that defines about 9% of all primary breast cancer tumours and is associated with high levels of genomic instability4-6. Specifically, inhibition of polo-like kinase4 (PLK4) using small molecules leads to centrosome depletion, which triggers mitotic catastrophe in cells that exhibit amplicon-directed overexpression of TRIM37. To explain this effect, we identify TRIM37 as a negative regulator of centrosomal pericentriolar material. In 17q23-amplified cells that lack centrosomes, increased levels of TRIM37 block the formation of foci that comprise pericentriolar material-these foci are structures with a microtubule-nucleating capacity that are required for successful cell division in the absence of centrosomes. Finally, we find thatthe overexpression of TRIM37 causes genomic instability by delaying centrosome maturation and separation at mitotic entry, and thereby increases the frequency of mitotic errors. Collectively, these findings highlight TRIM37-dependent genomic instability as a putative driver event in 17q23-amplified breast cancer and provide a rationale for the use of centrosome-targeting therapeutic agents in treating these cancers.

Gamma Knife Radiosurgery does not alter the copy number aberration profile in sporadic vestibular schwannoma

IntroductionIonizing radiation is a known etiologic factor in tumorigenesis and its role in inducing malignancy in the treatment of vestibular schwannoma has been debated. The purpose of this study was to identify a copy number aberration (CNA) profile or specific CNAs associated with radiation exposure which could either implicate an increased risk of malignancy or elucidate a mechanism of treatment resistance.Methods55 sporadic VS, including 18 treated with Gamma Knife Radiosurgery (GKRS), were subjected to DNA whole-genome microarray and/or whole-exome sequencing. CNAs were called and statistical tests were performed to identify any association with radiation exposure. Hierarchical clustering was used to identify CNA profiles associated with radiation exposure.ResultsA median of 7 (0–58) CNAs were identified across the 55 VS. Chromosome 22 aberration was the only recurrent event. A median aberrant cell fraction of 0.59 (0.25–0.94) was observed, indicating several genetic clones in VS. No CNA or CNA profile was associated with GKRS.ConclusionGKRS is not associated with an increase in CNAs or alteration of the CNA profile in VS, lending support to its low risk. This also implies that there is no major issue with GKRS treatment failure being due to CNAs. In agreement with previous studies, chromosome 22 aberration is the only recurrent CNA. VS consist of several genetic clones, addressing the need for further studies on the composition of cells in this tumor.

Pancreatic ductal adenocarcinomas from Mexican patients present a distinct genomic mutational pattern.

Pancreatic ductal adenocarcinoma (PDAC) is one of the deadliest cancers in humans, with less than 5% 5-year survival rate. PDAC is characterized by a small number of recurrent mutations, including KRAS, CDKN2A, TP53, and SMAD4 and a long "tail" of infrequent mutated genes. Most of the studies have been performed in US and European populations, so new studies are needed to describe the mutational landscape of these tumors in other cohorts. The present study analyzed the exome and transcriptome of four PDAC tumors from Mexican patients. We found a paucity of the previously described recurrent mutations, with mutations in only three genes (HERC2, CNTNAP2 and HMCN1) previously reported in PDAC with a frequency > 1%. In addition, we discovered several recurrent putative copy number aberrations in SKP2, BRAF, CSSF1R, FOXE1, JAK2 and MET genes and in genes previously reported as putative drivers in PDAC, including KRAS, SF3B1, BRAF, MYC and MET. Although a larger cohort is needed to validate these findings, our results could be pointing toward potential differences in contributing factors for PDAC in Latin-American populations.

Array Comparative Genomic Hybridization Analysis Reveals Significantly Enriched Pathways in Canine Oral Melanoma.

Human Mucosal Melanoma (hMM) is an aggressive neoplasm of neuroectodermal origin with distinctive features from the more common cutaneous form of malignant melanoma (cMM). At the molecular level, hMMs are characterized by large chromosomal aberrations rather than single-nucleotide mutations, typically observed in cMM. Given the scarcity of available cases, there have been many attempts to establish a reliable animal model. In pet dogs, Canine Oral Melanoma (COM) is the most common malignant tumor of the oral cavity, sharing clinical and histological aspects with hMM. To improve the knowledge about COM's genomic DNA alterations, in the present work, formalin-fixed, paraffin-embedded (FFPE) samples of COM from different European archives were collected to set up an array Comparative Genomic Hybridization (aCGH) analysis to estimate recurrent Copy Number Aberrations (CNAs). DNA was extracted in parallel from tumor and healthy fractions and 19 specimens were successfully submitted to labeling and competitive hybridization. Data were statistically analyzed through GISTIC2.0 and a pathway-enrichment analysis was performed with ClueGO. Recurrent gained regions were detected, affecting chromosomes CFA 10, 13 and 30, while lost regions involved chromosomes CFA 10, 11, 22, and 30. In particular, CFA 13 showed a whole-chromosome gain in 37% of the samples, while CFA 22 showed a whole-chromosome loss in 25%. A distinctive sigmoidal trend was observed in CFA 10 and 30 in 25 and 30% of the samples, respectively. Comparative analysis revealed that COM and hMM share common chromosomal changes in 32 regions. MAPK- and PI3K-related genes were the most frequently involved, while pathway analysis revealed statistically significant perturbation of cancer-related biological processes such as immune response, drug metabolism, melanocytes homeostasis, and neo-angiogenesis. The latter is a new evidence of a significant involvement of neovascularization-related pathways in COMs and can provide the rationale for future application in anti-cancer targeted therapies.

Whole Genome Profiling of Adult B-Other Acute Lymphoblastic Leukaemia on the UKALL14 Trial

B-other ALL represents a working definition for patients with B cell precursor (BCP) ALL without a known primary chromosomal abnormality. In this study we use whole genome sequencing (WGS) to characterize adult B-other cases (age ≥ 25yrs) from the UKALL14 trial (NCT01085617). Of 652 patients aged 25-65yrs enrolled onto UKALL14, 333 (51%) had B-other ALL. Sufficient material was available to screen 156/333 B-other cases for recurrent Ph-like fusion events (CLRF2, JAK2, ABL1, ABL2PDGFRB) using FISH and MLPA (kit P335). This identified 28 (18%) Ph-like fusion events (21 CRLF2, 5 ABL-class fusions and 2 JAK). Of the remaining 128 B-other cases 57 had available samples for tumor normal paired WGS (read depth 60x and 30x respectively). Bioinformatic analysis was performed to determine small somatic mutations (SSMs); single nucleotide variants (SNVs) and insertion/deletions (INDELs) as well as copy number aberrations (CNA) and structural variants (SV). We also undertook de novo motif analysis to identify RAG mediated deletions. We present data for 30/57 cases (median age of diagnosis 44, range 25-65), the remainder are undergoing sequencing. Within this cohort we identified 784 SVs, 49,244 SNVs and 2,881 INDELs, with each case having a median representation of 23 SVs, 1,674 SNVs and 86 INDELs. The Median SSM burden was 0.55 per megabase (range 0.31-0.81), which is in the upper third of previous ALL estimates (median 0.26 range 0.03-2.9) but low compared to most other cancer types (Alexandrov et al. 2013 Nature). Fusion gene analysis identified 11/30 (37%) cases with recurrent rearrangements (2 ZNF384r, 2 PAX5r, 5 DUX4r and FGFR1r); this identified a novel 5' partner for ZNF384r (AKAP8) and MYO18A-FGFR1 a fusion previously seen in a single case of 8p11 myeloproliferative syndrome. A further 7 clonal coding drivers were identified, 3 PAX5alt, 3 Ph-like candidates and a single ETV6-RUNX1-like candidate (ETV6 R399C, a dominant negative mutation). We classified our cohort into 6 driver subgroups as shown in table 1 using Gu et al. as a framework (Gu et al. 2019 Nat Gen). We found no evidence of driver subgroups clustering with age, although the Ph-like candidates were all identified in patients 39yrs or older. There was evidence that known driver subgroups have differing mutation burdens but these were not significant in this preliminary cohort; the single MEF2Dr case had a high SV burden compared to all other known subgroups (137 vs. cohort median 22); Favorable outcome subgroups (ZNF384r and DUX4r n=7) had a lower SNV coding mutation burden (median 14 vs. 20; U = 60, p = 0.057; Mann-Whitney U test); PAX5alt (n=5) cases had a higher INDEL burden (median 174 vs. 82; U = 16.5, p =0.057). As expected we found recurrent CNA in CDKN2A/B (63%; 19/30), PAX5 (33%; 10/30), IKZF1 (27%; 8/30), ETV6 (11%; 3/30) and BTG1 (7%; 2/30) across the cohort. IKZF1 deletions were enriched in the Ph-like candidate (n=4) subgroup compared to all other known groups (75% vs. 12%; p = 0.028; Fisher's Exact). RAG mediated deletions have been established as an oncogenic driver mechanism in childhood ALL, so we sought to ascertain its role in adults. We identified 54% (128/236) of the breakpoint resolved deletions had a RAG heptamer site. There was a significant difference (U = 167.5, p = 0.021) in RAG deletion burden between patients over 40yrs and under 40yrs at diagnosis (median 1.5 vs 3). Within the unknown driver subgroup; three cases carried ZEB2 H1038R mutations, two with concurrent IGH-CEBPA/B fusions, identifying them as likely belonging to the G12 cluster identified by Li et al. (Li et al. 2018 PNAS); one case had a high translocation burden (37 cohort median 2) and a single case had a high RAG deletion burden (45, cohort median 2); of the remaining 7 cases, four had either a NRAS or FLT3 subclonal hotspot mutation. Here we present the WGS analyses of 30 cases classified as B-other on the basis of no cytogenetic findings by standard clinical assays. The landscape of adult B-other ALL is highly heterogeneous but with WGS we were able to find at least one disease defining event in 60% of our cohort. These events often encompass novel fusion partners of established genomic subtypes or a cytogenetically cryptic lesion such as IGH-DUX4. Taken together our findings demonstrate the clinical utility of WGS as a diagnostic assay to inform and improve the management of adult B-other ALL patients in the future. Disclosures Fielding: Pfizer: Consultancy; Incyte: Consultancy; Amgen: Consultancy; Novartis: Consultancy. Papaemmanuil:Celgene: Research Funding.

New Recurrent Structural Aberrations in the Genome of Chronic Lymphocytic Leukemia Based on Exome-Sequencing Data.

Chronic lymphocytic leukemia (CLL) is the most frequent lymphoproliferative syndrome in Western countries, and it is characterized by recurrent large genomic rearrangements. During the last decades, array techniques have expanded our knowledge about CLL’s karyotypic aberrations. The advent of large sequencing databases expanded our knowledge cancer genomics to an unprecedented resolution and enabled the detection of small-scale structural aberrations in the cancer genome. In this study, we have performed exome-sequencing-based copy number aberration (CNA) and loss of heterozygosity (LOH) analysis in order to detect new recurrent structural aberrations. We describe 54 recurrent focal CNAs enriched in cancer-related pathways, and their association with gene expression and clinical evolution. Furthermore, we discovered recurrent large copy number neutral LOH events affecting key driver genes, and we recapitulate most of the large CNAs that characterize the CLL genome. These results provide “proof-of-concept” evidence supporting the existence of new genes involved in the pathogenesis of CLL.

Mytype: A Capture Based Sequencing Approach to Detect Somatic Mutations, Copy Number Changes and IGH Translocations in Multiple Myeloma

IntroductionAt diagnosis, Multiple Myeloma (MM) is traditionally classified into two clinical and prognostic subgroups groups on the basis of initiating cytogenetic abnormalities: IGH translocations and hyperdiploidy. Currently, these events are clinically ascertained by Fluorescent In-Situ Hybridization (FISH). In recent years, comprehensive genome profiling studies have shown that MM pathogenesis is defined by a spectrum of acquired somatic lesions, many of which are biologically and clinically relevant. To this effect, targeted gene sequencing approaches are becoming routine in the upfront diagnostic settings. Here we present myTYPE, a MM-specific targeted next generation sequencing panel to identify germline and somatic substitutions, indels, Copy Number Aberrations (CNA) and IGH translocations.MethodsA multiplex bait panel was designed to capture the exons of 120 genes implicated in MM pathogenesis, entire IGH locus as well as genome wide representation of single nucleotide polymorphisms (SNPs) (1 in 3Mb) to enable detection of arm level copy number events and recurrent focal events. These 120 genes were selected on the basis of 1) frequently mutated and driver genes in MM 2) genes in important signaling pathways, e.g the NFKB pathway 3) treatment targets and candidate genes for drug resistance, e.g. cereblon.To validate the efficacy of the assay, 16 constitutional bone marrow samples and 18 tumor samples were sequenced using myTYPE. For validation, 6/18 tumor/normal pairs sequenced using myTYPE were subject to WGS and remaining 12/18 tumor samples were subject to FISH. After sequencing, we obtained an overall median target coverage of 815x.ResultsAfter alignment, substitutions and indels were called using Caveman, Pindel and Strelka. CNAs were identified using Facets and IGH translocations were identified using Delly along with a modified version of BRASS. Below is a description of the genomic abnormalities captured by the myTYPE assay.SNVs and IndelsFor the 6 tumor/normal pairs sequenced using myTYPE and WGS, we obtained a total of 21 (median = 3) non-synonymous mutations using myTYPE. When limiting the WGS calls to myTYPE targets, we recovered 20/21 non-synonymous mutations identified by myTYPE. These involved SNVs and indels in key MM related drivers including NRAS, KRAS, FAM46C and TP53 among others. For the mutations identified by both myTYPE and WGS, there was a high correlation between the variant VAFs, R2 = 0.99 and as expected is better in capturing subclonal mutations.IGH rearrangements and Copy Number Aberrations (CNA)Next we compared myTYPE and WGS results for recurrent CNAs in MM. We specifically looked at deletions of 1p, 13p, 16q, 17p and gains of 1q, 11q and found a 100% concordance of these aberrations identified by both assays.The remaining 12 samples sequenced using myTYPE also had orthogonal FISH. myTYPE identified a total of 7 IGH rearrangements, 4 of which are also reported by FISH. Three additional t(11;14) translocations were uniquely identified by myTYPE in cases that remained clinically uncharacterized. FISH was also used to probe deletions in 17q, 13q, 1p and 1q gain. All aberrations identified by FISH were also identified in myType. Additionally, 13q- in four samples and 1p- in one sample were uniquely identified by myTYPE.ConclusionIn summary, we present a targeted assay capable of identifying somatic mutations, CNAs and IGH translocations of prognostic and diagnostic relevance in MM. When compared to conventional assays currently used in clinical practice, myTYPE identified at least one disease defining alterations in all samples screened. Evaluation of sensitivity and specificity will require larger clinical cohorts. Importantly, myTYPE enables comprehensive profiling, large sample multiplexing and short turn around times which renders it as an optimal assay for utilisation in the upfront clinical setting. DisclosuresKorde:Amgen: Research Funding. Mailankody:Juno: Research Funding; Janssen: Research Funding; Takeda: Research Funding; Physician Education Resource: Honoraria. Landgren:Pfizer: Consultancy; Amgen: Consultancy, Research Funding; Celgene: Consultancy, Research Funding; Janssen: Consultancy, Membership on an entity's Board of Directors or advisory committees, Research Funding; Merck: Membership on an entity's Board of Directors or advisory committees; Takeda: Consultancy, Membership on an entity's Board of Directors or advisory committees, Research Funding; Karyopharm: Consultancy.

Unraveling the Landscape of Copy Number Aberrations in Hodgkin Lymphoma: A Joint KU Leuven and Lysa Study on Circulating Cell Free DNA

Introduction. For decades, the study of genomic alterations in Hodgkin lymphoma (HL) has been hampered by the low abundance of the malignant Hodgkin Reed-Sternberg (HRS) cells in HL lymph node biopsies. Laser microdissection or flow cytometric cell sorting allow the study of purified HRS cells, but the use of these applications is restricted to specialized research centers. We have recently demonstrated proof-of-principle that copy number aberrations (CNA) in HRS cells can be retrieved in circulating cell-free DNA (cfDNA) of HL patients (Vandenberghe et al., Lancet Haematol. 2015). Taking advantage of the presence of HRS cell-derived DNA (ctDNA) in plasma, the aim of this study was to catalogue CNA in HL in a large series of prospectively recruited HL patients.Methods. We analyzed plasma collected from 177 patients (median age 29, range 3-86) with newly diagnosed HL. 60 cases were diagnosed at our institution, including all disease subtypes and stages, with a majority of nodular sclerosis. 118 patients were recruited in the context of the BREACH study, a multicentric Phase 2B study for unfavorable early classical HL (NCT02292979). After cfDNA extraction, samples were low-pass sequenced (0.1x coverage) and analyzed using ichorCNA, an algorithm that produces read depth-based log2 CNA profiles and quantifies the cfDNA tumor content.Results. At diagnosis, 86 % (152/177) of patients showed obvious genomic imbalances in cfDNA, in early-stage (85 % (131/155)) as well as in advanced cases (95 % (21/22)). Among the abnormal profiles, gains encompassing 2p16, 5p15, 9p24, 12q13 and 19q13 were the most frequent aberrations, occurring in 79 %, 53 %, 57 %, 63 % and 80 % respectively. Losses most frequently affected regions 4q34, 6q23, 11q22 and 13q13, in 53 %, 63 %, 49 % and 59 % of abnormal profiles respectively (fig. 1A). All these CNAs have been previously described with varying frequencies in smaller series of 10-53 patients, using arrayCGH or whole exome sequencing on microdissected HRS cells or even HRS cell-derived cell lines as input material. The observed pattern of CNAs is distinct from the pattern we observed in cfDNA obtained from other hematological malignances, e.g. multiple myeloma and diffuse large B-cell lymphoma (fig. 1B, C).Analysis of follow-up cfDNA samples revealed that the majority of patients rapidly clear their profiles during the first two treatment cycles. Of 123 samples analysed at d15 of cycle 1, 89 samples (72 %) showed a normal profile without CNAs.This suggests a relationship between the fraction of ctDNA in cfDNA and the disease burden. However, we could not find a correlation between the ctDNA fraction as calculated by ichorCNA and the metabolic tumor volume as computed from PET/CT scans. This could be due to several factors, e.g. the small fraction of HRS cells in the metabolic tumor volume. We are currently investigating whether a correlation does exist between the HRS content in the lymph node biopsy and the ctDNA load.Conclusions. In this largest study of CNA in HL to date, we provide a comprehensive catalogue of the types of CNAs, as well as their frequencies and patterns in HL. In this series of 177 patients, gains affecting 2p, 12q and 19q and loss of 6q and 13q emerge as the most commonly recurrent CNAs in HL, across all HL subtypes and stages. The data confirm and extend our previous findings that the majority of HL patients, including those with early-stage disease, have detectable CNAs in their cfDNA at diagnosis. We demonstrate a broad spectrum of aberrations, gains and losses, some of which recur at higher frequencies than gains of 9p24, harboring PD-L1. This warrants further studies on how these CNAs are implicated in the pathogenesis of HL. It further endorses the use of ctDNA as an alternative gateway to the genome of HRS cells, and as a substrate for the evaluation of early disease response. [Display omitted] DisclosuresNo relevant conflicts of interest to declare.

HetRCNA: A Novel Method to Identify Recurrent Copy Number Alternations from Heterogeneous Tumor Samples Based on Matrix Decomposition Framework.

A common strategy to discovering cancer associated copy number aberrations (CNAs) from a cohort of cancer samples is to detect recurrent CNAs (RCNAs). Although the previous methods can successfully identify communal RCNAs shared by nearly all tumor samples, detecting subgroup-specific RCNAs and their related subgroup samples from cancer samples with heterogeneity is still invalid for these existing approaches. In this paper, we introduce a novel integrated method called HetRCNA, which can identify statistically significant subgroup-specific RCNAs and their related subgroup samples. Based on matrix decomposition framework with weight constraint, HetRCNA can successfully measure the subgroup samples by coefficients of left vectors with weight constraint and subgroup-specific RCNAs by coefficients of the right vectors and significance test. When we evaluate HetRCNA on simulated dataset, the results show that HetRCNA gives the best performances among the competing methods and is robust to the noise factors of the simulated data. When HetRCNA is applied on a real breast cancer dataset, our approach successfully identifies a bunch of RCNA regions and the result is highly correlated with the results of the other two investigated approaches. Notably, the genomic regions identified by HetRCNA harbor many breast cancer related genes reported by previous researches.

Analysis of high-throughput biological data using their rank values

High-throughput biological technologies are routinely used to generate gene expression profiling or cytogenetics data. To achieve high performance, methods available in the literature become more specialized and often require high computational resources. Here, we propose a new versatile method based on the data-ordering rank values. We use linear algebra, the Perron-Frobenius theorem and also extend a method presented earlier for searching differentially expressed genes for the detection of recurrent copy number aberration. A result derived from the proposed method is a one-sample Student's t-test based on rank values. The proposed method is to our knowledge the only that applies to gene expression profiling and to cytogenetics data sets. This new method is fast, deterministic, and requires a low computational load. Probabilities are associated with genes to allow a statistically significant subset selection in the data set. Stability scores are also introduced as quality parameters. The performance and comparative analyses were carried out using real data sets. The proposed method can be accessed through an R package available from the CRAN (Comprehensive R Archive Network) website: https://cran.r-project.org/web/packages/fcros .