Complex Structural Variation Associated with Enhancer Hijacking and Loss of Tumor Suppressors in Mantle Cell Lymphoma
Bibliographic record
Abstract
Abstract Mantle cell lymphoma (MCL) is a rare, incurable mature B-cell lymphoma that can have either an aggressive or indolent clinical course. The hallmark aberration of MCL is the t(11;14)(q13;q32) translocation that places CCND1 under control of the immunoglobulin heavy chain (IGH) locus resulting in its constitutive expression. CCND1/IGH translocation occurs in nearly all MCL cases and arises during VDJ recombination. Whole genome sequencing (WGS) of MCLs has shown the prevalence of additional structural variations (SV), particularly in tumors with poor outcome. Complex rearrangements, such as chromothripsis and chromoplexy have been observed in MCL but their role in lymphomagenesis has not been determined. We hypothesized that some of these complex rearrangements afford a selective advantage to the tumor by disrupting tumor suppressor genes or by placing oncogenes in proximity to regulatory elements in cis, thereby resulting in ectopic expression. We performed WGS on tumor DNA extracted from 106 RCHOP-treated MCL patients. Matching ribosomal-depleted RNA-seq data was available for all tumor samples. Gene expression count values were obtained using Salmon and counts were normalized using the DESeq2 method. Structural variants were identified using a consensus between GRIDSS and Manta and these were analyzed to identify the topology of complex rearrangements using JaBbA. This allows the annotation of rearrangements such as chromothripsis and chromoplexy and enables the detection of genes near new regulatory elements through complex (multi-breakpoint) rearrangements. We supplemented this analysis by including 57 published classical and non-nodal leukemic MCL genomes. Chromothripsis and chromoplexy were observed in 8 and 12 tumors, respectively. The majority of the genomes were associated with isolated structural variations such as large deletions, inversions and reciprocal translocations were more common. Six genes (BCL10, TRAF6, TRAF3, MAP3K7, BTK, RELB) coding for canonical and non-canonical NFκB signaling proteins were found rearranged to within 1Mbp of a naïve B-cell super-enhancer region. The TRAF6 rearrangement was part of a chromothripsis and chromoplexy event involving both chr9 and chr11. A separate chromothripsis example involving chr1 placed each of BCL10 and NOTCH2 within 100-200kbp of super-enhancers. A 1Mbp region containing RELB was amplified and inserted into chr19p approximately 150kbp downstream a super-enhancer. An 800kbp deletion brought MAP3K7 to within 400kbp of a super-enhancer. DAZAP1, a gene known to be recurrently mutated in MCL, was translocated upstream to within 300kbp of a super-enhancer by t(5;19)(p13.3;p15.33). TRAF3 was translocated 400kbp upstream a super-enhancer by t(14;20)(q32.32;q13.13). None of the aforementioned SVs were associated with a detectable increase in expression of the affected gene. In contrast, we identified one genome in which the MYC oncogene was relocated 500kbp upstream of a super-enhancer via an unbalanced t(4;8)(q21.23;q24.21) translocation. In this case, MYC expression was in the 95 th percentile of MYC expression across the cohort. Focal deletions and amplifications were also found affecting lymphoma driver genes. Focal amplifications affecting 3q (34 tumors) and 5p (9 tumors), were among the most common recurrent events. These respectively affect the TERC and TERT genes, both involved in telomerase function. All tumors with TERT amplifications also showed TERC amplifications. TERC and TERT were expressed higher on average in amplified tumors than in unamplified tumors. Two tumors showed focal deletions affecting the 3' end of BIRC3. The focal deletion in one tumor was found to span to the 3' end of BIRC2 resulting in a BIRC2/BIRC3 fusion. Aberrant splicing across the two genes was evident in matching tumor RNA-seq data. Complex rearrangements in MCL have been found to link distant super-enhancer elements with a variety of lymphoma oncogenes. We noted a recurrence of such events affecting known regulators of NFκB signaling. We are using nanopore-based long-read PromethION sequencing to validate the structure of the derivative chromosome in these cases. Although these genes were not detectably overexpressed, deregulation of genes may be occurring by other means. The full extent of deregulation of NFκB and other oncogenic pathways will be revealed as complex rearrangements are studied in additional MCL tumors. Disclosures Coyle: Allakos, Inc.: Consultancy. Melnick: Epizyme: Consultancy; Daiichi Sankyo: Research Funding; Sanofi: Research Funding; Janssen Pharmaceuticals: Research Funding; Constellation: Consultancy; KDAC Pharma: Membership on an entity's Board of Directors or advisory committees. Scott: BC Cancer: Patents & Royalties: Patent describing assigning DLBCL COO by gene expression profiling--licensed to NanoString Technologies. Patent describing measuring the proliferation signature in MCL using gene expression profiling. ; AstraZeneca: Consultancy; Abbvie: Consultancy; NanoString Technologies: Patents & Royalties: Patent describing measuring the proliferation signature in MCL using gene expression profiling.; Rich/Genentech: Research Funding; Celgene: Consultancy; Incyte: Consultancy; Janssen: Consultancy, Research Funding. Morin: Epizyme: Patents & Royalties; Celgene: Consultancy; Foundation for Burkitt Lymphoma Research: Membership on an entity's Board of Directors or advisory committees.
Fetched live from OpenAlex and de-inverted. Abstracts are not stored in this database: the inverted indexes are 8.6 GB of the frame’s 9.3 GB of text, and the host has 13 GB free.
How this classification was reachedexpand
Full frame machine prediction
Teacher imitationNot calibrated prevalence, not ground truth. Human validation pending. The Gemma side is a direct model label for every work in the frame, read from the title-only record. The Codex side is a classifier learned from the 10,348 direct Codex labels and calibrated to design-weighted sample rates; fields without enough sample support carry no Codex call. Candidate is the union of the two sides; consensus is their intersection. These outputs are machine_predicted_unvalidated and are not human labels.
Distilled classifier scores by category (both heads)
| Category | Codex | Gemma |
|---|---|---|
| Metaresearch | 0.000 | 0.000 |
| Meta-epidemiology (narrow) | 0.000 | 0.000 |
| Meta-epidemiology (broad) | 0.000 | 0.000 |
| Bibliometrics | 0.001 | 0.000 |
| Science and technology studies | 0.000 | 0.000 |
| Scholarly communication | 0.000 | 0.000 |
| Open science | 0.000 | 0.000 |
| Research integrity | 0.000 | 0.000 |
| Insufficient payload (model declined to judge) | 0.002 | 0.000 |
Machine scores (provisional)
The two teacher heads of the student model, read on this work. A score orders the frame for review; it never asserts a category, and the validation status ships verbatim with every row.
Baseline scores from an immature model (maturity gate not passed, 7 training rounds). Scores rank; they never assert a category.
score_only:v0-immature-baseline · verbatim from the scoring run: score_only means the number may rank works, and no category label ships from itClassification
machine, unvalidatedMachine predicted; a candidate call from one source (direct Gemma or distilled Codex), not a consensus.
How this classification was reached, model by model and score by score, is at the end of the page under "How this classification was reached".