13q12.2 deletions in acute lymphoblastic leukemia lead to upregulation of FLT3 through enhancer hijacking

Developmental origins shape the paediatric cancer genome

In the past two decades, technological advances have brought unprecedented insights into the paediatric cancer genome revealing characteristics distinct from those of adult cancer. Originating from developing tissues, paediatric cancers generally have low mutation burden and are driven by variants that disrupt the transcriptional activity, chromatin state, non-coding cis-regulatory regions and other biological functions. Within each tumour, there are multiple populations of cells with varying states, and the lineages of some can be tracked to their fetal origins. Genome-wide genetic screening has identified vulnerabilities associated with both the cell of origin and transcription deregulation in paediatric cancer, which have become a valuable resource for designing new therapeutic approaches including those for small molecules, immunotherapy and targeted protein degradation. In this Review, we present recent findings on these facets of paediatric cancer from a pan-cancer perspective and provide an outlook on future investigations.

Characterisation of FLT3 alterations in childhood acute lymphoblastic leukaemia

BACKGROUND

Alterations of FLT3 are among the most common driver events in acute leukaemia with important clinical implications, since it allows patient classification into prognostic groups and the possibility of personalising therapy thanks to the availability of FLT3 inhibitors. Most of the knowledge on FLT3 implications comes from the study of acute myeloid leukaemia and so far, few studies have been performed in other leukaemias.

METHODS

A comprehensive genomic (DNA-seq in 267 patients) and transcriptomic (RNA-seq in 160 patients) analysis of FLT3 in 342 childhood acute lymphoblastic leukaemia (ALL) patients was performed. Mutations were functionally characterised by in vitro experiments.

RESULTS

Point mutations (PM) and internal tandem duplications (ITD) were detected in 4.3% and 2.7% of the patients, respectively. A new activating mutation of the TKD, G846D, conferred oncogenic properties and sorafenib resistance. Moreover, a novel alteration involving the circularisation of read-through transcripts (rt-circRNAs) was observed in 10% of the cases. Patients presenting FLT3 alterations exhibited higher levels of the receptor. In addition, patients with ZNF384- and MLL/KMT2A-rearranged ALL, as well as hyperdiploid subtype, overexpressed FLT3.

DISCUSSION

Our results suggest that specific ALL subgroups may also benefit from a deeper understanding of the biology of FLT3 alterations and their clinical implications.

Disruption of the TP53 locus in osteosarcoma leads to TP53 promoter gene fusions and restoration of parts of the TP53 signalling pathway

TP53 is the most frequently mutated gene in human cancer. This gene shows not only loss-of-function mutations but also recurrent missense mutations with gain-of-function activity. We have studied the primary bone malignancy osteosarcoma, which harbours one of the most rearranged genomes of all cancers. This is odd since it primarily affects children and adolescents who have not lived the long life thought necessary to accumulate massive numbers of mutations. In osteosarcoma, TP53 is often disrupted by structural variants. Here, we show through combined whole-genome and transcriptome analyses of 148 osteosarcomas that TP53 structural variants commonly result in loss of coding parts of the gene while simultaneously preserving and relocating the promoter region. The transferred TP53 promoter region is fused to genes previously implicated in cancer development. Paradoxically, these erroneously upregulated genes are significantly associated with the TP53 signalling pathway itself. This suggests that while the classical tumour suppressor activities of TP53 are lost, certain parts of the TP53 signalling pathway that are necessary for cancer cell survival and proliferation are retained. In line with this, our data suggest that transposition of the TP53 promoter is an early event that allows for a new normal state of genome-wide rearrangements in osteosarcoma. © 2023 The Authors. The Journal of Pathology published by John Wiley & Sons Ltd on behalf of The Pathological Society of Great Britain and Ireland.

Cancer Evolution: A Multifaceted Affair

Cancer cells adapt and survive through the acquisition and selection of molecular modifications. This process defines cancer evolution. Building on a theoretical framework based on heritable genetic changes has provided insights into the mechanisms supporting cancer evolution. However, cancer hallmarks also emerge via heritable nongenetic mechanisms, including epigenetic and chromatin topological changes, and interactions between tumor cells and the tumor microenvironment. Recent findings on tumor evolutionary mechanisms draw a multifaceted picture where heterogeneous forces interact and influence each other while shaping tumor progression. A comprehensive characterization of the cancer evolutionary toolkit is required to improve personalized medicine and biomarker discovery.

SIGNIFICANCE

Tumor evolution is fueled by multiple enabling mechanisms. Importantly, genetic instability, epigenetic reprogramming, and interactions with the tumor microenvironment are neither alternative nor independent evolutionary mechanisms. As demonstrated by findings highlighted in this perspective, experimental and theoretical approaches must account for multiple evolutionary mechanisms and their interactions to ultimately understand, predict, and steer tumor evolution.

Spatial promoter-enhancer hubs in cancer: organization, regulation, and function

Transcriptional dysregulation is a hallmark of cancer and can be driven by altered enhancer landscapes. Recent studies in genome organization have revealed that multiple enhancers and promoters can spatially coalesce to form dynamic topological assemblies, known as promoter-enhancer hubs, which strongly correlate with elevated gene expression. In this review, we discuss the structure and complexity of promoter-enhancer hubs recently identified in multiple cancer types. We further discuss underlying mechanisms driving dysregulation of promoter-enhancer hubs and speculate on their functional role in pathogenesis. Understanding the role of promoter-enhancer hubs in transcriptional dysregulation can provide insight into new therapeutic approaches to target these complex features of genome organization.

Systematic characterization of the HOXA9 downstream targets in MLL-r leukemia by noncoding CRISPR screens

Accumulating evidence indicates that HOXA9 dysregulation is necessary and sufficient for leukemic transformation and maintenance. However, it remains largely unknown how HOXA9, as a homeobox transcriptional factor, binds to noncoding regulatory sequences and controls the downstream genes. Here, we conduct dropout CRISPR screens against 229 HOXA9-bound peaks identified by ChIP-seq. Integrative data analysis identifies reproducible noncoding hits, including those located in the distal enhancer of FLT3 and intron of CDK6. The Cas9-editing and dCas9-KRAB silencing of the HOXA9-bound sites significantly reduce corresponding gene transcription and impair cell proliferation in vitro, and in vivo by transplantation into NSG female mice. In addition, RNA-seq, Q-PCR analysis, chromatin accessibility change, and chromatin conformation evaluation uncover the noncoding regulation mechanism of HOXA9 and its functional downstream genes. In summary, our work improves our understanding of how HOXA9-associated transcription programs reconstruct the regulatory network specifying MLL-r dependency.

Optical Genome Mapping Identifies Novel Recurrent Structural Alterations in Childhood ETV6::RUNX1+ and High Hyperdiploid Acute Lymphoblastic Leukemia

The mutational landscape of B-cell precursor acute lymphoblastic leukemia (BCP-ALL), the most common pediatric cancer, is not fully described partially because commonly applied short-read next generation sequencing has a limited ability to identify structural variations. By combining comprehensive analysis of structural variants (SVs), single-nucleotide variants (SNVs), and small insertions-deletions, new subtype-defining and therapeutic targets may be detected. We analyzed the landscape of somatic alterations in 60 pediatric patients diagnosed with the most common BCP-ALL subtypes, ETV6::RUNX1+ and classical hyperdiploid (HD), using conventional cytogenetics, single nucleotide polymorphism (SNP) array, whole exome sequencing (WES), and the novel optical genome mapping (OGM) technique. Ninety-five percent of SVs detected by cytogenetics and SNP-array were verified by OGM. OGM detected an additional 677 SVs not identified using the conventional methods, including (subclonal) IKZF1 deletions. Based on OGM, ETV6::RUNX1+ BCP-ALL harbored 2.7 times more SVs than HD BCP-ALL, mainly focal deletions. Besides SVs in known leukemia development genes (ETV6, PAX5, BTG1, CDKN2A), we identified 19 novel recurrently altered regions (in n ≥ 3) including 9p21.3 (FOCAD/HACD4), 8p11.21 (IKBKB), 1p34.3 (ZMYM1), 4q24 (MANBA), 8p23.1 (MSRA), and 10p14 (SFMBT2), as well as ETV6::RUNX1+ subtype-specific SVs (12p13.1 (GPRC5A), 12q24.21 (MED13L), 18q11.2 (MIB1), 20q11.22 (NCOA6)). We detected 3 novel fusion genes (SFMBT2::DGKD, PDS5B::STAG2, and TDRD5::LPCAT2), for which the sequence and expression were validated by long-read and whole transcriptome sequencing, respectively. OGM and WES identified double hits of SVs and SNVs (ETV6, BTG1, STAG2, MANBA, TBL1XR1, NSD2) in the same patient demonstrating the power of the combined approach to define the landscape of genomic alterations in BCP-ALL.

Transcription factor networks link B-lymphocyte development and malignant transformation in leukemia

Rapid advances in genomics have opened unprecedented possibilities to explore the mutational landscapes in malignant diseases, such as B-cell acute lymphoblastic leukemia (B-ALL). This disease is manifested as a severe defect in the production of normal blood cells due to the uncontrolled expansion of transformed B-lymphocyte progenitors in the bone marrow. Even though classical genetics identified translocations of transcription factor-coding genes in B-ALL, the extent of the targeting of regulatory networks in malignant transformation was not evident until the emergence of large-scale genomic analyses. There is now evidence that many B-ALL cases present with mutations in genes that encode transcription factors with critical roles in normal B-lymphocyte development. These include PAX5, IKZF1, EBF1, and TCF3, all of which are targeted by translocations or, more commonly, partial inactivation in cases of B-ALL. Even though there is support for the notion that germline polymorphisms in the PAX5 and IKZF1 genes predispose for B-ALL, the majority of leukemias present with somatic mutations in transcription factor-encoding genes. These genetic aberrations are often found in combination with mutations in genes that encode components of the pre-B-cell receptor or the IL-7/TSLP signaling pathways, all of which are important for early B-cell development. This review provides an overview of our current understanding of the molecular interplay that occurs between transcription factors and signaling events during normal and malignant B-lymphocyte development.

Spatial and clonality-resolved 3D cancer genome alterations reveal enhancer-hijacking as a potential prognostic marker for colorectal cancer

The regulatory effect of non-coding large-scale structural variations (SVs) on proto-oncogene activation remains unclear. This study investigated SV-mediated gene dysregulation by profiling 3D cancer genome maps from 40 patients with colorectal cancer (CRC). We developed a machine learning-based method for spatial characterization of the altered 3D cancer genome. This revealed a frequent establishment of "de novo chromatin contacts" that can span multiple topologically associating domains (TADs) in addition to the canonical TAD fusion/shuffle model. Using this information, we precisely identified super-enhancer (SE)-hijacking and its clonal characteristics. Clonal SE-hijacking genes, such as TOP2B, are recurrently associated with cell-cycle/DNA-processing functions, which can potentially be used as CRC prognostic markers. Oncogene activation and increased drug resistance due to SE-hijacking were validated by reconstructing the patient's SV using CRISPR-Cas9. Collectively, the spatial and clonality-resolved analysis of the 3D cancer genome reveals regulatory principles of large-scale SVs in oncogene activation and their clinical implications.

Enhancer Coamplification and Hijacking Promote Oncogene Expression in Liposarcoma

SIGNIFICANCE

Comprehensive profiling of the enhancer landscape and 3D genome structure in liposarcoma identifies extensive enhancer-oncogene coamplification and enhancer hijacking events, deepening the understanding of how oncogenes are regulated in cancer.

Genomic structural variation: A complex but important driver of human evolution

Structural variants (SVs)-including duplications, deletions, and inversions of DNA-can have significant genomic and functional impacts but are technically difficult to identify and assay compared with single-nucleotide variants. With the aid of new genomic technologies, it has become clear that SVs account for significant differences across and within species. This phenomenon is particularly well-documented for humans and other primates due to the wealth of sequence data available. In great apes, SVs affect a larger number of nucleotides than single-nucleotide variants, with many identified SVs exhibiting population and species specificity. In this review, we highlight the importance of SVs in human evolution by (1) how they have shaped great ape genomes resulting in sensitized regions associated with traits and diseases, (2) their impact on gene functions and regulation, which subsequently has played a role in natural selection, and (3) the role of gene duplications in human brain evolution. We further discuss how to incorporate SVs in research, including the strengths and limitations of various genomic approaches. Finally, we propose future considerations in integrating existing data and biospecimens with the ever-expanding SV compendium propelled by biotechnology advancements.

Hyperdiploidy: the longest known, most prevalent, and most enigmatic form of acute lymphoblastic leukemia in children

Hyperdiploidy is the largest genetic entity B-cell precursor acute lymphoblastic leukemia in children. The diagnostic hallmark of its two variants that will be discussed in detail herein is a chromosome count between 52 and 67, respectively. The classical HD form consists of heterozygous di-, tri-, and tetrasomies, whereas the nonclassical one (usually viewed as "duplicated hyperhaploid") contains only disomies and tetrasomies. Despite their apparently different clinical behavior, we show that these two sub-forms can in principle be produced by the same chromosomal maldistribution mechanism. Moreover, their respective array, gene expression, and mutation patterns also indicate that they are biologically more similar than hitherto appreciated. Even though in-depth analyses of the genomic intricacies of classical HD leukemias are indispensable for the elucidation of the disease process, the ensuing results play at present surprisingly little role in treatment stratification, a fact that can be attributed to the overall good prognoses and low relapse rates of the concerned patients and, consequently, their excellent treatment outcome. Irrespective of this underutilization, however, the detailed genetic characterization of HD leukemias may, especially in planned treatment reduction trials, eventually become important for further treatment stratification, patient management, and the clinical elucidation of outcome data. It should therefore become an integral part of all upcoming treatment studies.

Oncogenic lesions and molecular subtypes in adults with B-cell acute lymphoblastic leukemia

B-cell acute lymphoblastic leukemia (B-ALL), a genetically heterogeneous disease, is classified into different molecular subtypes that are defined by recurrent gene rearrangements, gross chromosomal abnormalities, or specific gene mutations. Cells with these genetic alterations acquire a leukemia-initiating ability and show unique expression profiles. The distribution of B-ALL molecular subtypes is greatly dependent on age, which also affects treatment responsiveness and long-term survival, partly accounting for the inferior outcome in adolescents and young adults (AYA) and (older) adults with B-ALL. Recent advances in sequencing technology, especially RNA sequencing and the application of these technologies in large B-ALL cohorts have uncovered B-ALL molecular subtypes prevalent in AYA and adults. These new insights supply more precise estimations of prognoses and targeted therapies informed by sequencing results, as well as a deeper understanding of the genetic basis of AYA/adult B-ALL. This article provides an account of these technological advances and an overview of the recent major findings of B-ALL molecular subtypes in adults.

Three-dimensional genome organization in immune cell fate and function

Immune cell development and activation demand the precise and coordinated control of transcriptional programmes. Three-dimensional (3D) organization of the genome has emerged as an important regulator of chromatin state, transcriptional activity and cell identity by facilitating or impeding long-range genomic interactions among regulatory elements and genes. Chromatin folding thus enables cell type-specific and stimulus-specific transcriptional responses to extracellular signals, which are essential for the control of immune cell fate, for inflammatory responses and for generating a diverse repertoire of antigen receptor specificities. Here, we review recent findings connecting 3D genome organization to the control of immune cell differentiation and function, and discuss how alterations in genome folding may lead to immune dysfunction and malignancy.

Epigenetic activation of the FLT3 gene by ZNF384 fusion confers a therapeutic susceptibility in acute lymphoblastic leukemia

FLT3 is an attractive therapeutic target in acute lymphoblastic leukemia (ALL) but the mechanism for its activation in this cancer is incompletely understood. Profiling global gene expression in large ALL cohorts, we identify over-expression of FLT3 in ZNF384-rearranged ALL, consistently across cases harboring different fusion partners with ZNF384. Mechanistically, we discover an intergenic enhancer element at the FLT3 locus that is exclusively activated in ZNF384-rearranged ALL, with the enhancer-promoter looping directly mediated by the fusion protein. There is also a global enrichment of active enhancers within ZNF384 binding sites across the genome in ZNF384-rearranged ALL cells. Downregulation of ZNF384 blunts FLT3 activation and decreases ALL cell sensitivity to FLT3 inhibitor gilteritinib in vitro. In patient-derived xenograft models of ZNF384-rearranged ALL, gilteritinib exhibits significant anti-leukemia efficacy as a monotherapy in vivo. Collectively, our results provide insights into FLT3 regulation in ALL and point to potential genomics-guided targeted therapy for this patient population.

The genomic landscape of pediatric acute lymphoblastic leukemia

Acute lymphoblastic leukemia (ALL) is the most common childhood cancer. Here, using whole-genome, exome and transcriptome sequencing of 2,754 childhood patients with ALL, we find that, despite a generally low mutation burden, ALL cases harbor a median of four putative somatic driver alterations per sample, with 376 putative driver genes identified varying in prevalence across ALL subtypes. Most samples harbor at least one rare gene alteration, including 70 putative cancer driver genes associated with ubiquitination, SUMOylation, noncoding transcripts and other functions. In hyperdiploid B-ALL, chromosomal gains are acquired early and synchronously before ultraviolet-induced mutation. By contrast, ultraviolet-induced mutations precede chromosomal gains in B-ALL cases with intrachromosomal amplification of chromosome 21. We also demonstrate the prognostic significance of genetic alterations within subtypes. Intriguingly, DUX4- and KMT2A-rearranged subtypes separate into CEBPA/FLT3- or NFATC4-expressing subgroups with potential clinical implications. Together, these results deepen understanding of the ALL genomic landscape and associated outcomes.

Enhancer retargeting of CDX2 and UBTF::ATXN7L3 define a subtype of high-risk B-progenitor acute lymphoblastic leukemia

Transcriptome sequencing has identified multiple subtypes of B-progenitor acute lymphoblastic leukemia (B-ALL) of prognostic significance, but a minority of cases lack a known genetic driver. Here, we used integrated whole-genome (WGS) and -transcriptome sequencing (RNA-seq), enhancer mapping, and chromatin topology analysis to identify previously unrecognized genomic drivers in B-ALL. Newly diagnosed (n = 3221) and relapsed (n = 177) B-ALL cases with tumor RNA-seq were studied. WGS was performed to detect mutations, structural variants, and copy number alterations. Integrated analysis of histone 3 lysine 27 acetylation and chromatin looping was performed using HiChIP. We identified a subset of 17 newly diagnosed and 5 relapsed B-ALL cases with a distinct gene expression profile and 2 universal and unique genomic alterations resulting from aberrant recombination-activating gene activation: a focal deletion downstream of PAN3 at 13q12.2 resulting in CDX2 deregulation by the PAN3 enhancer and a focal deletion of exons 18-21 of UBTF at 17q21.31 resulting in a chimeric fusion, UBTF::ATXN7L3. A subset of cases also had rearrangement and increased expression of the PAX5 gene, which is otherwise uncommon in B-ALL. Patients were more commonly female and young adult with median age 35 (range,12-70 years). The immunophenotype was characterized by CD10 negativity and immunoglobulin M positivity. Among 16 patients with known clinical response, 9 (56.3%) had high-risk features including relapse (n = 4) or minimal residual disease >1% at the end of remission induction (n = 5). CDX2-deregulated, UBTF::ATXN7L3 rearranged (CDX2/UBTF) B-ALL is a high-risk subtype of leukemia in young adults for which novel therapeutic approaches are required.

Construction of three-gene-based prognostic signature and analysis of immune cells infiltration in children and young adults with B-acute lymphoblastic leukemia

Background

Although B‐acute lymphoblastic leukemia (B‐ALL) patients' survival has been improved dramatically, some cases still relapse. This study aimed to explore the prognosis‐related novel differentially expressed genes (DEGs) for predicting the overall survival (OS) of children and young adults (CAYAs) with B‐ALL and analyze the immune‐related factors contributing to poor prognosis.

Methods

GSE48558 and GSE79533 from Gene Expression Omnibus (GEO) and clinical sample information and mRNA‐seq from Therapeutically Applicable Research to Generate Effective Treatments (TARGET) database were retrieved. Prognosis‐related key genes were enrolled to build a Cox proportional model using multivariate Cox regression. Five‐year OS of patients, clinical characteristic relevance and clinical independence were assessed based on the model. The mRNA levels of prognosis‐related genes were validated in our samples and the difference of immune cells composition between high‐risk and low‐risk patients were compared.

Results

One hundred and twelve DEGs between normal B cells and B‐ALL cells were identified based on GSE datasets. They were mainly participated in protein binding and HIF‐1 signaling pathway. One hundred and eighty‐nine clinical samples were enrolled in the study, both Kaplan–Meier (KM) analysis and univariate Cox regression analysis showed that CYBB, BCL2A1, IFI30, and EFNB1 were associated with prognosis, CYBB, BCL2A1, and EFNB1 were used to construct prognostic risk model. Moreover, compared to clinical indicators, the three‐gene signature was an independent prognostic factor for CAYAs with B‐ALL. Finally, the mRNA levels of CYBB, BCL2A1, and EFNB1 were significantly lower in B‐ALL group as compared to controls. The high‐risk group had a significantly higher percentage of infiltrated immune cells.

Conclusion

We constructed a novel three‐gene signature with independent prognostic factor for predicting 5‐year OS of CAYAs with B‐ALL. Additionally, we discovered the difference of immune cells composition between high‐risk and low‐risk groups. This study may help to customize individual treatment and improve prognosis of CAYAs with B‐ALL.

Alterations in transcriptional networks in cancer: the role of noncoding somatic driver mutations

Aberrant gene expression is a cancer hallmark and it is known that almost every tumor acquires somatic mutations in transcription factors, chromatin regulators, or the DNA regulatory elements that are critical for transcriptional control and cell phenotype. While the role of transcription factors and chromatin regulators has been widely studied, relatively few noncoding driver mutations have been identified and functionally characterized to date. Here, we review the current understanding of somatic variants in noncoding regions of the cancer genome and their impact on chromatin architecture and transcriptional networks. We also discuss approaches and ongoing challenges for noncoding driver discovery, and highlight insights gained from recent studies exploring the nature and impact of noncoding drivers on tumor formation.

3D chromatin architecture and transcription regulation in cancer

Chromatin has distinct three-dimensional (3D) architectures important in key biological processes, such as cell cycle, replication, differentiation, and transcription regulation. In turn, aberrant 3D structures play a vital role in developing abnormalities and diseases such as cancer. This review discusses key 3D chromatin structures (topologically associating domain, lamina-associated domain, and enhancer-promoter interactions) and corresponding structural protein elements mediating 3D chromatin interactions [CCCTC-binding factor, polycomb group protein, cohesin, and Brother of the Regulator of Imprinted Sites (BORIS) protein] with a highlight of their associations with cancer. We also summarise the recent development of technologies and bioinformatics approaches to study the 3D chromatin interactions in gene expression regulation, including crosslinking and proximity ligation methods in the bulk cell population (ChIA-PET and HiChIP) or single-molecule resolution (ChIA-drop), and methods other than proximity ligation, such as GAM, SPRITE, and super-resolution microscopy techniques.

An alternative CYB5A transcript is expressed in aneuploid ALL and enriched in relapse

Background

B-cell precursor acute lymphoblastic leukemia (BCP-ALL) is a genetically heterogenous malignancy with poor prognosis in relapsed adult patients. The genetic basis for relapse in aneuploid subtypes such as near haploid (NH) and high hyperdiploid (HeH) BCP-ALL is only poorly understood. Pathogenic genetic alterations remain to be identified. To this end, we investigated the dynamics of genetic alterations in a matched initial diagnosis-relapse (ID-REL) BCP-ALL cohort. Here, we firstly report the identification of the novel genetic alteration CYB5Aalt, an alternative transcript of CYB5A, in two independent cohorts.

Methods

We identified CYB5alt in the RNAseq-analysis of a matched ID-REL BCP-ALL cohort with 50 patients and quantified its expression in various molecular BCP-ALL subtypes. Findings were validated in an independent cohort of 140 first diagnosis samples from adult BCP-ALL patients. Derived from patient material, the alternative open reading frame of CYB5Aalt was cloned (pCYB5Aalt) and pCYB5Aalt or the empty vector were stably overexpressed in NALM-6 cells. RNA sequencing was performed of pCYB5Aalt clones and empty vector controls followed by differential expression analysis, gene set enrichment analysis and complementing cell death and viability assays to determine functional implications of CYB5Aalt.

Results

RNAseq data analysis revealed non-canonical exon usage of CYB5Aalt starting from a previously undescribed transcription start site. CYB5Aalt expression was increased in relapsed BCP-ALL and its occurrence was specific towards the shared gene expression cluster of NH and HeH BCP-ALL in independent cohorts. Overexpression of pCYB5Aalt in NALM-6 cells induced a distinct transcriptional program compared to empty vector controls with downregulation of pathways related to reported functions of CYB5A wildtype. Interestingly, CYB5A wildtype expression was decreased in CYB5Aalt samples in silico and in vitro. Additionally, pCYB5Aalt NALM-6 elicited a more resistant drug response.

Conclusions

Across all age groups, CYB5Aalt was the most frequent secondary genetic event in relapsed NH and HeH BCP-ALL. In addition to its high subgroup specificity, CYB5Aalt is a novel candidate to be potentially implicated in therapy resistance in NH and HeH BCP-ALL. This is underlined by overexpressing CYB5Aalt providing first evidence for a functional role in BCL2-mediated apoptosis.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12863-022-01041-1.

Enhancer RNAs (eRNAs) in Cancer: The Jacks of All Trades

Simple Summary

This review focuses on eRNAs and the several mechanisms by which they can regulate gene expression. In particular we describe here the most recent examples of eRNAs dysregulated in cancer or involved in the immune escape of tumor cells.

Abstract

Enhancer RNAs (eRNAs) are non-coding RNAs (ncRNAs) transcribed in enhancer regions. They play an important role in transcriptional regulation, mainly during cellular differentiation. eRNAs are tightly tissue- and cell-type specific and are induced by specific stimuli, activating promoters of target genes in turn. eRNAs usually have a very short half-life but in some cases, once activated, they can be stably expressed and acquire additional functions. Due to their critical role, eRNAs are often dysregulated in cancer and growing number of interactions with chromatin modifiers, transcription factors, and splicing machinery have been described. Enhancer activation and eRNA transcription have particular relevance also in inflammatory response, placing the eRNAs at the interplay between cancer and immune cells. Here, we summarize all the possible molecular mechanisms recently reported in association with eRNAs activity.

ProbC: joint modeling of epigenome and transcriptome effects in 3D genome

Background

Hi-C and its high nucleosome resolution variant Micro-C provide a window into the spatial packing of a genome in 3D within the cell. Even though both techniques do not directly depend on the binding of specific antibodies, previous work has revealed enriched interactions and domain structures around multiple chromatin marks; epigenetic modifications and transcription factor binding sites. However, the joint impact of chromatin marks in Hi-C and Micro-C interactions have not been globally characterized, which limits our understanding of 3D genome characteristics. An emerging question is whether it is possible to deduce 3D genome characteristics and interactions by integrative analysis of multiple chromatin marks and associate interactions to functionality of the interacting loci.

Result

We come up with a probabilistic method ProbC to decompose Hi-C and Micro-C interactions by known chromatin marks. ProbC is based on convex likelihood optimization, which can directly take into account both interaction existence and nonexistence. Through ProbC, we discover histone modifications (H3K27ac, H3K9me3, H3K4me3, H3K4me1) and CTCF as particularly predictive of Hi-C and Micro-C contacts across cell types and species. Moreover, histone modifications are more effective than transcription factor binding sites in explaining the genome’s 3D shape through these interactions. ProbC can successfully predict Hi-C and Micro-C interactions in given species, while it is trained on different cell types or species. For instance, it can predict missing nucleosome resolution Micro-C interactions in human ES cells trained on mouse ES cells only from these 5 chromatin marks with above 0.75 AUC. Additionally, ProbC outperforms the existing methods in predicting interactions across almost all chromosomes.

Conclusion

Via our proposed method, we optimally decompose Hi-C interactions in terms of these chromatin marks at genome and chromosome levels. We find a subset of histone modifications and transcription factor binding sites to be predictive of both Hi-C and Micro-C interactions and TADs across human, mouse, and different cell types. Through learned models, we can predict interactions on species just from chromatin marks for which Hi-C data may be limited.

Concurrent CDX2 cis-deregulation and UBTF-ATXN7L3 fusion define a novel high-risk subtype of B-cell ALL

CDX2 cis-deregulation and UBTF::ATXN7L3 fusion driven by focal deletions define a novel subtype of B-ALL.

CDX2/UBTF::ATXN7L3 is a high-risk B-ALL subtype in young adults, which warrants improved therapeutic strategies.

Oncogenic alterations underlying B-cell acute lymphoblastic leukemia (B-ALL) in adults remain incompletely elucidated. To uncover novel oncogenic drivers, we performed RNA sequencing and whole-genome analyses in a large cohort of unresolved B-ALL. We identified a novel subtype characterized by a distinct gene expression signature and the unique association of 2 genomic microdeletions. The 17q21.31 microdeletion resulted in a UBTF::ATXN7L3 fusion transcript encoding a chimeric protein. The 13q12.2 deletion resulted in monoallelic ectopic expression of the homeobox transcription factor CDX2, located 138 kb in cis from the deletion. Using 4C-sequencing and CRISPR interference experiments, we elucidated the mechanism of CDX2 cis-deregulation, involving PAN3 enhancer hijacking. CDX2/UBTF ALL (n = 26) harbored a distinct pattern of additional alterations including 1q gain and CXCR4 activating mutations. Within adult patients with Ph⁻ B-ALL enrolled in GRAALL trials, patients with CDX2/UBTF ALL (n = 17/723, 2.4%) were young (median age, 31 years) and dramatically enriched in females (male/female ratio, 0.2, P = .002). They commonly presented with a pro-B phenotype ALL and moderate blast cell infiltration. They had poor response to treatment including a higher risk of failure to first induction course (19% vs 3%, P = .017) and higher post-induction minimal residual disease (MRD) levels (MRD ≥ 10⁻⁴, 93% vs 46%, P < .001). This early resistance to treatment translated into a significantly higher cumulative incidence of relapse (75.0% vs 32.4%, P = .004) in univariate and multivariate analyses. In conclusion, we discovered a novel B-ALL entity defined by the unique combination of CDX2 cis-deregulation and UBTF::ATXN7L3 fusion, representing a high-risk disease in young adults.

Redefining the biological basis of lineage-ambiguous leukemia through genomics: BCL11B deregulation in acute leukemias of ambiguous lineage

Acute leukemias of ambiguous lineage (ALAL), including mixed phenotype acute leukemia (MPAL) and related entities such as early T-cell precursor acute leukemia (ETP-ALL), remain diagnostic and clinical challenges due to limited understanding of pathogenesis, reliance of immunophenotyping to classify disease, and the lack of a rational approach to guide selection of appropriate therapy. Recent studies utilizing genomic sequencing and complementary approaches have provided key insights that are changing the way in which such leukemias are classified, and potentially, treated. Several recurrent genomic alterations define leukemias that straddle immunophenotypic entities, such as ZNF384-rearranged childhood B-ALL and B/myeloid MPAL, and BCL11B-rearranged T/myeloid MPAL, ETP-ALL and AML. In contrast, some cases of MPAL represent canonical ALL/AML entities exhibiting lineage aberrancy. For many cases of ALAL, experimental approaches indicate lineage aberrancy arises from acquisition of a founding genetic alteration into a hematopoietic stem or progenitor cell. Determination of optimal therapeutic approach requires genomic characterization of uniformly treated ALAL patients in prospective studies, but several approaches, including kinase inhibitors and BH3 mimetics may be efficacious in subsets of ALAL.

Experimental strategies to achieve efficient targeted knock-in via tandem paired nicking

Tandem paired nicking (TPN) is a method of genome editing that enables precise and relatively efficient targeted knock-in without appreciable restraint by p53-mediated DNA damage response. TPN is initiated by introducing two site-specific nicks on the same DNA strand using Cas9 nickases in such a way that the nicks encompass the knock-in site and are located within a homologous region between a donor DNA and the genome. This nicking design results in the creation of two nicks on the donor DNA and two in the genome, leading to relatively efficient homology-directed recombination between these DNA fragments. In this study, we sought to identify the optimal design of TPN experiments that would improve the efficiency of targeted knock-in, using multiple reporter systems based on exogenous and endogenous genes. We found that efficient targeted knock-in via TPN is supported by the use of 1700-2000-bp donor DNAs, exactly 20-nt-long spacers predicted to be efficient in on-target cleavage, and tandem-paired Cas9 nickases nicking at positions close to each other. These findings will help establish a methodology for efficient and precise targeted knock-in based on TPN, which could broaden the applicability of targeted knock-in to various fields of life science.

Infectious triggers and novel therapeutic opportunities in childhood B cell leukaemia

B cell acute lymphoblastic leukaemia (B-ALL) is the most common form of childhood cancer. Although treatment has advanced remarkably in the past 50 years, it still fails in ~20% of patients. Recent studies revealed that more than 5% of healthy newborns carry preleukaemic clones that originate in utero, but only a small percentage of these carriers will progress to overt B-ALL. The drivers of progression are unclear, but B-ALL incidence seems to be increasing in parallel with the adoption of modern lifestyles. Emerging evidence shows that a major driver for the conversion from the preleukaemic state to the B-ALL state is exposure to immune stressors, such as infection. Here, we discuss our current understanding of the environmental triggers and genetic predispositions that may lead to B-ALL, highlighting lessons from epidemiology, the clinic and animal models, and identifying priority areas for future research.

Genome-wide detection of enhancer-hijacking events from chromatin interaction data in rearranged genomes

Recent efforts have shown that structural variations (SVs) can disrupt three-dimensional genome organization and induce enhancer hijacking, yet no computational tools exist to identify such events from chromatin interaction data. Here, we develop NeoLoopFinder, a computational framework to identify the chromatin interactions induced by SVs, including interchromosomal translocations, large deletions and inversions. Our framework can automatically resolve complex SVs, reconstruct local Hi-C maps surrounding the breakpoints, normalize copy number variation and allele effects and predict chromatin loops induced by SVs. We applied NeoLoopFinder in Hi-C data from 50 cancer cell lines and primary tumors and identified tens of recurrent genes associated with enhancer hijacking. To experimentally validate NeoLoopFinder, we deleted the hijacked enhancers in prostate adenocarcinoma cells using CRISPR-Cas9, which significantly reduced expression of the target oncogene. In summary, NeoLoopFinder enables identification of critical oncogenic regulatory elements that can potentially reveal therapeutic targets.

Deciphering the Complexity of 3D Chromatin Organization Driving Lymphopoiesis and Lymphoid Malignancies

Proper lymphopoiesis and immune responses depend on the spatiotemporal control of multiple processes, including gene expression, DNA recombination and cell fate decisions. High-order 3D chromatin organization is increasingly appreciated as an important regulator of these processes and dysregulation of genomic architecture has been linked to various immune disorders, including lymphoid malignancies. In this review, we present the general principles of the 3D chromatin topology and its dynamic reorganization during various steps of B and T lymphocyte development and activation. We also discuss functional interconnections between architectural, epigenetic and transcriptional changes and introduce major key players of genomic organization in B/T lymphocytes. Finally, we present how alterations in architectural factors and/or 3D genome organization are linked to dysregulation of the lymphopoietic transcriptional program and ultimately to hematological malignancies.

Epigenetic Control of Infant B Cell Precursor Acute Lymphoblastic Leukemia

B-cell precursor acute lymphoblastic leukemia (BCP-ALL) is a highly aggressive malignancy, with poorer prognosis in infants than in adults. A genetic signature has been associated with this outcome but, remarkably, leukemogenesis is commonly triggered by genetic alterations of embryonic origin that involve the deregulation of chromatin remodelers. This review considers in depth how the alteration of epigenetic profiles (at DNA and histone levels) induces an aberrant phenotype in B lymphocyte progenitors by modulating the oncogenic drivers and tumor suppressors involved in key cancer hallmarks. DNA methylation patterns have been widely studied in BCP-ALL and their correlation with survival has been established. However, the effect of methylation on histone residues can be very different. For instance, methyltransferase KMT2A gene participates in chromosomal rearrangements with several partners, imposing an altered pattern of methylated H3K4 and H3K79 residues, enhancing oncogene promoter activation, and conferring a worse outcome on affected infants. In parallel, acetylation processes provide an additional layer of epigenetic regulation and can alter the chromatin conformation, enabling the binding of regulatory factors. Therefore, an integrated knowledge of all epigenetic disorders is essential to understand the molecular basis of BCP-ALL and to identify novel entry points that can be exploited to improve therapeutic options and disease prognosis.

Exploration of Coding and Non-coding Variants in Cancer Using GenomePaint

GenomePaint (https://genomepaint.stjude.cloud/) is an interactive visualization platform for whole-genome, whole-exome, transcriptome, and epigenomic data of tumor samples. Its design captures the inter-relatedness between DNA variations and RNA expression, supporting in-depth exploration of both individual cancer genomes and full cohorts. Regulatory non-coding variants can be inspected and analyzed along with coding variants, and their functional impact further explored by examining 3D genome data from cancer cell lines. Further, GenomePaint correlates mutation and expression patterns with patient outcomes, and supports custom data upload. We used GenomePaint to unveil aberrant splicing that disrupts the RING domain of CREBBP, discover cis activation of the MYC oncogene by duplication of the NOTCH1-MYC enhancer in B-lineage acute lymphoblastic leukemia, and explore the inter- and intra-tumor heterogeneity at EGFR in adult glioblastomas. These examples demonstrate that deep multi-omics exploration of individual cancer genomes enabled by GenomePaint can lead to biological insights for follow-up validation.

Determinants and role of chromatin organization in acute leukemia

DNA is compacted into higher order structures that have major implications in gene regulation. These structures allow for long-range interactions of DNA elements, such as the association of promoters with their cognate enhancers. In recent years, mutations in genes that control these structures, including the cohesin-complex and the insulator-binding protein CTCF, have been found in a spectrum of hematologic disorders, and especially in acute leukemias. Cohesin and CTCF are critical for mediating looping and establishing boundaries within chromatin. Cells that harbor mutations in these genes display aberrant chromatin architecture and resulting differences in gene expression that contribute to leukemia initiation and progression. Here, we provide detailed discussion of the nature of 3D interactions and the way that they are disrupted in acute leukemia. Continued research in this area will provide new insights into the mechanisms of leukemogenesis and may shed light on novel treatment strategies.