Genome sequencing of pediatric medulloblastoma links catastrophic DNA rearrangements with TP53 mutations.

Genomic rearrangements are thought to occur progressively during tumor development. Recent findings, however, suggest an alternative mechanism, involving massive chromosome rearrangements in a one-step catastrophic event termed chromothripsis. We report the whole-genome sequencing-based analysis of a Sonic-Hedgehog medulloblastoma (SHH-MB) brain tumor from a patient with a germline TP53 mutation (Li-Fraumeni syndrome), uncovering massive, complex chromosome rearrangements. Integrating TP53 status with microarray and deep sequencing-based DNA rearrangement data in additional patients reveals a striking association between TP53 mutation and chromothripsis in SHH-MBs. Analysis of additional tumor entities substantiates a link between TP53 mutation and chromothripsis, and indicates a context-specific role for p53 in catastrophic DNA rearrangements. Among these, we observed a strong association between somatic TP53 mutations and chromothripsis in acute myeloid leukemia. These findings connect p53 status and chromothripsis in specific tumor types, providing a genetic basis for understanding particularly aggressive subtypes of cancer.

Patterns of somatic structural variation in human cancer genomes.

A key mutational process in cancer is structural variation, in which rearrangements delete, amplify or reorder genomic segments that range in size from kilobases to whole chromosomes1-7. Here we develop methods to group, classify and describe somatic structural variants, using data from the Pan-Cancer Analysis of Whole Genomes (PCAWG) Consortium of the International Cancer Genome Consortium (ICGC) and The Cancer Genome Atlas (TCGA), which aggregated whole-genome sequencing data from 2,658 cancers across 38 tumour types8. Sixteen signatures of structural variation emerged. Deletions have a multimodal size distribution, assort unevenly across tumour types and patients, are enriched in late-replicating regions and correlate with inversions. Tandem duplications also have a multimodal size distribution, but are enriched in early-replicating regions-as are unbalanced translocations. Replication-based mechanisms of rearrangement generate varied chromosomal structures with low-level copy-number gains and frequent inverted rearrangements. One prominent structure consists of 2-7 templates copied from distinct regions of the genome strung together within one locus. Such cycles of templated insertions correlate with tandem duplications, and-in liver cancer-frequently activate the telomerase gene TERT. A wide variety of rearrangement processes are active in cancer, which generate complex configurations of the genome upon which selection can act.

Somatic structural variant formation is guided by and influences genome architecture.

The occurrence and formation of genomic structural variants (SVs) is known to be influenced by the 3D chromatin architecture, but the extent and magnitude have been challenging to study. Here, we apply Hi-C to study chromatin organization before and after induction of chromothripsis in human cells. We use Hi-C to manually assemble the derivative chromosomes following the occurrence of massive complex rearrangements, which allows us to study the sources of SV formation and their consequences on gene regulation. We observe an action-reaction interplay whereby the 3D chromatin architecture directly impacts the location and formation of SVs. In turn, the SVs reshape the chromatin organization to alter the local topologies, replication timing, and gene regulation in cis We show that SVs have a strong tendency to occur between similar chromatin compartments and replication timing regions. Moreover, we find that SVs frequently occur at 3D loop anchors, that SVs can cause a switch in chromatin compartments and replication timing, and that this is a major source of SV-mediated effects on nearby gene expression changes. Finally, we provide evidence for a general mechanistic bias of the 3D chromatin on SV occurrence using data from more than 2700 patient-derived cancer genomes.

Somatic CpG hypermutation is associated with mismatch repair deficiency in cancer.

Somatic hypermutation in cancer has gained momentum with the increased use of tumour mutation burden as a biomarker for immune checkpoint inhibitors. Spontaneous deamination of 5-methylcytosine to thymine at CpG dinucleotides is one of the most ubiquitous endogenous mutational processes in normal and cancer cells. Here, we performed a systematic investigation of somatic CpG hypermutation at a pan-cancer level. We studied 30,191 cancer patients and 103 cancer types and developed an algorithm to identify somatic CpG hypermutation. Across cancer types, we observed the highest prevalence in paediatric leukaemia (3.5%), paediatric high-grade glioma (1.7%), and colorectal cancer (1%). We discovered germline variants and somatic mutations in the mismatch repair complex MutSα (MSH2-MSH6) as genetic drivers of somatic CpG hypermutation in cancer, which frequently converged on CpG sites and TP53 driver mutations. We further observe an association between somatic CpG hypermutation and response to immune checkpoint inhibitors. Overall, our study identified novel cancer types that display somatic CpG hypermutation, strong association with MutSα-deficiency, and potential utility in cancer immunotherapy.

17p13 (TP53) Deletions Are Associated With an Aggressive Phenotype but Unrelated to Patient Prognosis in Urothelial Bladder Carcinomas.

17p13 deletions including TP53 and other genes represent a common cause for reduced/lost p53 function in tumor cells. In this study, we analyzed the impact of 17p13 (TP53) deletions and p53 expression on tumor aggressiveness and patient prognosis in urothelial carcinoma. The 17p13 copy number status was analyzed by fluorescence in situ hybridization (FISH) on more than 2700 urothelial bladder carcinomas in a tissue microarray format. 17p13 deletion data were compared to p53 expression data measured by immunohistochemistry (IHC) in a previous study. Different types of p53 alterations were compared with tumor phenotype and clinical outcome data. Deletions of 17p13 occurred in 23% of 2185 analyzable carcinomas. The fraction of tumors with 17p13 deletions increased from pTa G2 low (9%) to pTa G3 (24%, p < 0.0001). In muscle-invasive carcinomas, 17p13 deletions were associated with advanced pT stage (p = 0.0246), but unrelated to patient prognosis (p > 0.5). 17p13 deletions were significantly related to p53 immunostaining (p = 0.0375). 17p13 deletions were most common in tumors with complete lack of p53 staining (31%), which supports the concept that many of these tumors have a complete loss of p53 function (p53 null phenotype). 17p13 deletions were also increased in tumors with high p53 staining (25%). In conclusion, 17p13 deletions were most commonly seen in p53 negative cancers, supporting their role as a cause for the p53 null phenotype in urothelial cancer. The association of 17p13 deletions with high grade and advanced pT stage may reflect increasing genomic instability going along with stage and grade progression.

The landscape of genomic alterations across childhood cancers.

Pan-cancer analyses that examine commonalities and differences among various cancer types have emerged as a powerful way to obtain novel insights into cancer biology. Here we present a comprehensive analysis of genetic alterations in a pan-cancer cohort including 961 tumours from children, adolescents, and young adults, comprising 24 distinct molecular types of cancer. Using a standardized workflow, we identified marked differences in terms of mutation frequency and significantly mutated genes in comparison to previously analysed adult cancers. Genetic alterations in 149 putative cancer driver genes separate the tumours into two classes: small mutation and structural/copy-number variant (correlating with germline variants). Structural variants, hyperdiploidy, and chromothripsis are linked to TP53 mutation status and mutational signatures. Our data suggest that 7-8% of the children in this cohort carry an unambiguous predisposing germline variant and that nearly 50% of paediatric neoplasms harbour a potentially druggable event, which is highly relevant for the design of future clinical trials.

Author Correction: The landscape of genomic alterations across childhood cancers.

In this Article, author Benedikt Brors was erroneously associated with affiliation number '8' (Department of Developmental Neurobiology, St Jude Children's Research Hospital, Memphis, Tennessee, USA); the author's two other affiliations (affiliations '3' and '7', both at the German Cancer Research Center (DKFZ)) were correct. This has been corrected online.

PD-L1 expression in tumor and inflammatory cells is associated with favorable tumor features and favorable prognosis in muscle-invasive urothelial carcinoma of the bladder not treated by immune checkpoint inhibitors.

BACKGROUND: A high level of PD-L1 expression is the most relevant predictive parameter for response to immune checkpoint inhibitor (CPI) therapy in urinary bladder cancer. Existing data on the relationship between PD-L1 expression and the natural course of disease are controversial and sparse. METHODS: To expand our understanding of the relationship between PD-L1 expression and parameters of cancer aggressiveness, PD-L1 was analyzed on tissue microarrays containing 2710 urothelial bladder carcinomas including 512 patients with follow-up data who underwent radical cystectomy and follow-up therapies in the pre-immune checkpoint inhibitor therapy era. RESULTS: Tumor cell positivity in ≥10% of cells were seen in 513 (20%) and an immune cell positivity occurred in 872 (34%) of 2566 interpretable cancers. PD-L1 positivity in tumor cells increased from pTaG2 low grade (0.9% positive) to pTaG3 high grade (4.1%; p = 0.0255) and was even higher in muscle-invasive (pT2-4) carcinomas (29.3%; p < 0.0001). However, within pT2-4 carcinomas, PD-L1 positivity was linked to low pT stage (p = 0.0028), pN0 (p < 0.0001), L0 status (p = 0.0005), and a better prognosis within 512 patients with cystectomy who never received CPIs (p = 0.0073 for tumor cells and p = 0.0086 for inflammatory cells). PD-L1 staining in inflammatory cells was significantly linked to PD-L1 staining in tumor cells (p < 0.0001) and both were linked to a positive p53 immunostaining (p < 0.0001). CONCLUSION: It cannot be fully excluded that the strong statistical link between PD-L1 status and favorable histological tumor features as well as better prognosis could influence the outcome of studies evaluating CPIs in muscle-invasive urothelial carcinoma.

The whole-genome landscape of medulloblastoma subtypes.

Current therapies for medulloblastoma, a highly malignant childhood brain tumour, impose debilitating effects on the developing child, and highlight the need for molecularly targeted treatments with reduced toxicity. Previous studies have been unable to identify the full spectrum of driver genes and molecular processes that operate in medulloblastoma subgroups. Here we analyse the somatic landscape across 491 sequenced medulloblastoma samples and the molecular heterogeneity among 1,256 epigenetically analysed cases, and identify subgroup-specific driver alterations that include previously undiscovered actionable targets. Driver mutations were confidently assigned to most patients belonging to Group 3 and Group 4 medulloblastoma subgroups, greatly enhancing previous knowledge. New molecular subtypes were differentially enriched for specific driver events, including hotspot in-frame insertions that target KBTBD4 and 'enhancer hijacking' events that activate PRDM6. Thus, the application of integrative genomics to an extensive cohort of clinical samples derived from a single childhood cancer entity revealed a series of cancer genes and biologically relevant subtype diversity that represent attractive therapeutic targets for the treatment of patients with medulloblastoma.

Clonal hematopoiesis in elderly twins: concordance, discordance, and mortality.

Clonal hematopoiesis (CH) of indeterminate potential (CHIP) is defined by mutations in myeloid cancer-associated genes with a variant allele frequency of at least 2%. Recent studies have suggested a possible genetic predisposition to CH. To further explore this phenomenon, we conducted a population-based study of 594 twins from 299 pairs aged 73 to 94 years, all with >20 years' follow-up. We sequenced DNA from peripheral blood with a customized 21-gene panel at a median coverage of 6179X. The casewise concordance rates for mutations were calculated to assess genetic predisposition. Mutations were identified in 214 (36%) of the twins. Whereas 20 twin pairs had mutations within the same genes, the exact same mutation was only observed in 2 twin pairs. No significant difference in casewise concordance between monozygotic and dizygotic twins was found for any specific gene, subgroup, or CHIP mutations overall, and no significant heritability could be detected. In pairs discordant for CHIP mutations, we tested if the affected twin died before the unaffected twin, as a direct measurement of the association of having CH when controlling for familial factors. A total of 127 twin pairs were discordant for carrying a mutation, and in 61 (48%) cases, the affected twin died first (P = .72). Overall, we did not find a genetic predisposition to CHIP mutations in this twin study. The previously described negative association of CHIP mutations on survival could not be confirmed in a direct comparison among twin pairs that were discordant for CHIP mutations.

SvABA: genome-wide detection of structural variants and indels by local assembly.

Structural variants (SVs), including small insertion and deletion variants (indels), are challenging to detect through standard alignment-based variant calling methods. Sequence assembly offers a powerful approach to identifying SVs, but is difficult to apply at scale genome-wide for SV detection due to its computational complexity and the difficulty of extracting SVs from assembly contigs. We describe SvABA, an efficient and accurate method for detecting SVs from short-read sequencing data using genome-wide local assembly with low memory and computing requirements. We evaluated SvABA's performance on the NA12878 human genome and in simulated and real cancer genomes. SvABA demonstrates superior sensitivity and specificity across a large spectrum of SVs and substantially improves detection performance for variants in the 20-300 bp range, compared with existing methods. SvABA also identifies complex somatic rearrangements with chains of short (<1000 bp) templated-sequence insertions copied from distant genomic regions. We applied SvABA to 344 cancer genomes from 11 cancer types and found that short templated-sequence insertions occur in ∼4% of all somatic rearrangements. Finally, we demonstrate that SvABA can identify sites of viral integration and cancer driver alterations containing medium-sized (50-300 bp) SVs.

Phenotypic impact of genomic structural variation: insights from and for human disease.

Genomic structural variants have long been implicated in phenotypic diversity and human disease, but dissecting the mechanisms by which they exert their functional impact has proven elusive. Recently however, developments in high-throughput DNA sequencing and chromosomal engineering technology have facilitated the analysis of structural variants in human populations and model systems in unprecedented detail. In this Review, we describe how structural variants can affect molecular and cellular processes, leading to complex organismal phenotypes, including human disease. We further present advances in delineating disease-causing elements that are affected by structural variants, and we discuss future directions for research on the functional consequences of structural variants.

Spectrum and prevalence of genetic predisposition in medulloblastoma: a retrospective genetic study and prospective validation in a clinical trial cohort.

BACKGROUND: Medulloblastoma is associated with rare hereditary cancer predisposition syndromes; however, consensus medulloblastoma predisposition genes have not been defined and screening guidelines for genetic counselling and testing for paediatric patients are not available. We aimed to assess and define these genes to provide evidence for future screening guidelines. METHODS: In this international, multicentre study, we analysed patients with medulloblastoma from retrospective cohorts (International Cancer Genome Consortium [ICGC] PedBrain, Medulloblastoma Advanced Genomics International Consortium [MAGIC], and the CEFALO series) and from prospective cohorts from four clinical studies (SJMB03, SJMB12, SJYC07, and I-HIT-MED). Whole-genome sequences and exome sequences from blood and tumour samples were analysed for rare damaging germline mutations in cancer predisposition genes. DNA methylation profiling was done to determine consensus molecular subgroups: WNT (MBWNT), SHH (MBSHH), group 3 (MBGroup3), and group 4 (MBGroup4). Medulloblastoma predisposition genes were predicted on the basis of rare variant burden tests against controls without a cancer diagnosis from the Exome Aggregation Consortium (ExAC). Previously defined somatic mutational signatures were used to further classify medulloblastoma genomes into two groups, a clock-like group (signatures 1 and 5) and a homologous recombination repair deficiency-like group (signatures 3 and 8), and chromothripsis was investigated using previously established criteria. Progression-free survival and overall survival were modelled for patients with a genetic predisposition to medulloblastoma. FINDINGS: We included a total of 1022 patients with medulloblastoma from the retrospective cohorts (n=673) and the four prospective studies (n=349), from whom blood samples (n=1022) and tumour samples (n=800) were analysed for germline mutations in 110 cancer predisposition genes. In our rare variant burden analysis, we compared these against 53 105 sequenced controls from ExAC and identified APC, BRCA2, PALB2, PTCH1, SUFU, and TP53 as consensus medulloblastoma predisposition genes according to our rare variant burden analysis and estimated that germline mutations accounted for 6% of medulloblastoma diagnoses in the retrospective cohort. The prevalence of genetic predispositions differed between molecular subgroups in the retrospective cohort and was highest for patients in the MBSHH subgroup (20% in the retrospective cohort). These estimates were replicated in the prospective clinical cohort (germline mutations accounted for 5% of medulloblastoma diagnoses, with the highest prevalence [14%] in the MBSHH subgroup). Patients with germline APC mutations developed MBWNT and accounted for most (five [71%] of seven) cases of MBWNT that had no somatic CTNNB1 exon 3 mutations. Patients with germline mutations in SUFU and PTCH1 mostly developed infant MBSHH. Germline TP53 mutations presented only in childhood patients in the MBSHH subgroup and explained more than half (eight [57%] of 14) of all chromothripsis events in this subgroup. Germline mutations in PALB2 and BRCA2 were observed across the MBSHH, MBGroup3, and MBGroup4 molecular subgroups and were associated with mutational signatures typical of homologous recombination repair deficiency. In patients with a genetic predisposition to medulloblastoma, 5-year progression-free survival was 52% (95% CI 40-69) and 5-year overall survival was 65% (95% CI 52-81); these survival estimates differed significantly across patients with germline mutations in different medulloblastoma predisposition genes. INTERPRETATION: Genetic counselling and testing should be used as a standard-of-care procedure in patients with MBWNT and MBSHH because these patients have the highest prevalence of damaging germline mutations in known cancer predisposition genes. We propose criteria for routine genetic screening for patients with medulloblastoma based on clinical and molecular tumour characteristics. FUNDING: German Cancer Aid; German Federal Ministry of Education and Research; German Childhood Cancer Foundation (Deutsche Kinderkrebsstiftung); European Research Council; National Institutes of Health; Canadian Institutes for Health Research; German Cancer Research Center; St Jude Comprehensive Cancer Center; American Lebanese Syrian Associated Charities; Swiss National Science Foundation; European Molecular Biology Organization; Cancer Research UK; Hertie Foundation; Alexander and Margaret Stewart Trust; V Foundation for Cancer Research; Sontag Foundation; Musicians Against Childhood Cancer; BC Cancer Foundation; Swedish Council for Health, Working Life and Welfare; Swedish Research Council; Swedish Cancer Society; the Swedish Radiation Protection Authority; Danish Strategic Research Council; Swiss Federal Office of Public Health; Swiss Research Foundation on Mobile Communication; Masaryk University; Ministry of Health of the Czech Republic; Research Council of Norway; Genome Canada; Genome BC; Terry Fox Research Institute; Ontario Institute for Cancer Research; Pediatric Oncology Group of Ontario; The Family of Kathleen Lorette and the Clark H Smith Brain Tumour Centre; Montreal Children's Hospital Foundation; The Hospital for Sick Children: Sonia and Arthur Labatt Brain Tumour Research Centre, Chief of Research Fund, Cancer Genetics Program, Garron Family Cancer Centre, MDT's Garron Family Endowment; BC Childhood Cancer Parents Association; Cure Search Foundation; Pediatric Brain Tumor Foundation; Brainchild; and the Government of Ontario.

Dissecting the genomic complexity underlying medulloblastoma.

Medulloblastoma is an aggressively growing tumour, arising in the cerebellum or medulla/brain stem. It is the most common malignant brain tumour in children, and shows tremendous biological and clinical heterogeneity. Despite recent treatment advances, approximately 40% of children experience tumour recurrence, and 30% will die from their disease. Those who survive often have a significantly reduced quality of life. Four tumour subgroups with distinct clinical, biological and genetic profiles are currently identified. WNT tumours, showing activated wingless pathway signalling, carry a favourable prognosis under current treatment regimens. SHH tumours show hedgehog pathway activation, and have an intermediate prognosis. Group 3 and 4 tumours are molecularly less well characterized, and also present the greatest clinical challenges. The full repertoire of genetic events driving this distinction, however, remains unclear. Here we describe an integrative deep-sequencing analysis of 125 tumour-normal pairs, conducted as part of the International Cancer Genome Consortium (ICGC) PedBrain Tumor Project. Tetraploidy was identified as a frequent early event in Group 3 and 4 tumours, and a positive correlation between patient age and mutation rate was observed. Several recurrent mutations were identified, both in known medulloblastoma-related genes (CTNNB1, PTCH1, MLL2, SMARCA4) and in genes not previously linked to this tumour (DDX3X, CTDNEP1, KDM6A, TBR1), often in subgroup-specific patterns. RNA sequencing confirmed these alterations, and revealed the expression of what are, to our knowledge, the first medulloblastoma fusion genes identified. Chromatin modifiers were frequently altered across all subgroups. These findings enhance our understanding of the genomic complexity and heterogeneity underlying medulloblastoma, and provide several potential targets for new therapeutics, especially for Group 3 and 4 patients.

Genomes of early onset prostate cancer.

PURPOSE OF REVIEW: Prostate cancer is a disease of the elderly but a clinically relevant subset occurs early in life. In the current review, we discuss recent findings and the current understanding of the molecular underpinnings associated with early-onset prostate cancer (PCa) and the evidence supporting age-specific differences in the cancer genomes. RECENT FINDINGS: Recent surveys of PCa patient cohorts have provided novel age-dependent links between germline and somatic aberrations which points to differences in the molecular cause and treatment options. SUMMARY: Identifying the earliest molecular alterations in PCa can provide insight into the cause of the disease and biomarkers for patient risk stratification. Genomic aberrations of early-onset PCas display several patterns distinct from late-onset PCa genomes, suggesting age-dependent pathomechanisms involving alterations in the androgen receptor pathway.

A cell-based model system links chromothripsis with hyperploidy.

A remarkable observation emerging from recent cancer genome analyses is the identification of chromothripsis as a one-off genomic catastrophe, resulting in massive somatic DNA structural rearrangements (SRs). Largely due to lack of suitable model systems, the mechanistic basis of chromothripsis has remained elusive. We developed an integrative method termed "complex alterations after selection and transformation (CAST)," enabling efficient in vitro generation of complex DNA rearrangements including chromothripsis, using cell perturbations coupled with a strong selection barrier followed by massively parallel sequencing. We employed this methodology to characterize catastrophic SR formation processes, their temporal sequence, and their impact on gene expression and cell division. Our in vitro system uncovered a propensity of chromothripsis to occur in cells with damaged telomeres, and in particular in hyperploid cells. Analysis of primary medulloblastoma cancer genomes verified the link between hyperploidy and chromothripsis in vivo. CAST provides the foundation for mechanistic dissection of complex DNA rearrangement processes.

Heterogeneity of ERG expression in prostate cancer: a large section mapping study of entire prostatectomy specimens from 125 patients.

BACKGROUND: TMPRSS2:ERG fusions are frequent in prostate cancer, and occur predominantly in young patients. Several studies had proposed intratumoral heterogeneity of these fusions. This study was designed to determine frequency and extent of ERG fusion heterogeneity in early-onset prostate cancer (EO-PCA, <50 years) and in elderly patients. METHODS: The prostates from 63 EO-PCA and 62 elderly prostate cancer patients were thoroughly reviewed for presence of cancer foci. All 1592 tumor-containing sections were analyzed by immunohistochemistry for ERG expression. RESULTS: The prostates included in this study contained one tumor focus in 44, two tumor foci in 21, three tumor foci in 32, four tumor foci in 15, and five or more tumor foci in 13 patients. Among 59 cancer foci with ≤3 mm, 19 (32.2 %) were homogeneously ERG positive, 39 66.1 %) were homogeneously ERG negative, and one case (1.7 %) showed a heterogeneous ERG status. The fraction of homogeneously ERG positive cancer foci remained largely constant (14-37 %) with increasing tumor focus diameter but the fraction of heterogeneous ERG findings continuously increased with tumor size and reached 39 % in cancer foci larger than 22 mm. On a patient level, ERG expression was markedly more frequent in EO-PCA than in elderly patients: 13 % of EO-PCA were homogeneously and 62 % were heterogeneously ERG positive. In elderly patients, 3 % of cancers were homogeneously and 57 % were heterogeneously ERG positive (p = 0.0721). CONCLUSION: These data show that about 20-30 % of prostate cancer foci have early ERG fusions. ERG fusions further occur in about 50 % of initially ERG negative cancer foci during cancer progression. The vast majority of cancers are heterogeneous for TMPRSS2:ERG fusions on a patient level, challenging the concept of classifying prostate cancer patients into "fusion type" and "non-fusion type" prostate cancer.

Pediatric T-cell lymphoblastic leukemia evolves into relapse by clonal selection, acquisition of mutations and promoter hypomethylation.

Relapsed precursor T-cell acute lymphoblastic leukemia is characterized by resistance against chemotherapy and is frequently fatal. We aimed at understanding the molecular mechanisms resulting in relapse of T-cell acute lymphoblastic leukemia and analyzed 13 patients at first diagnosis, remission and relapse by whole exome sequencing, targeted ultra-deep sequencing, multiplex ligation dependent probe amplification and DNA methylation array. Compared to primary T-cell acute lymphoblastic leukemia, in relapse the number of single nucleotide variants and small insertions and deletions approximately doubled from 11.5 to 26. Targeted ultra-deep sequencing sensitively detected subclones that were selected for in relapse. The mutational pattern defined two types of relapses. While both are characterized by selection of subclones and acquisition of novel mutations, 'type 1' relapse derives from the primary leukemia whereas 'type 2' relapse originates from a common pre-leukemic ancestor. Relapse-specific changes included activation of the nucleotidase NT5C2 resulting in resistance to chemotherapy and mutations of epigenetic modulators, exemplified by SUZ12, WHSC1 and SMARCA4. While mutations present in primary leukemia and in relapse were enriched for known drivers of leukemia, relapse-specific changes revealed an association with general cancer-promoting mechanisms. This study thus identifies mechanisms that drive progression of pediatric T-cell acute lymphoblastic leukemia to relapse and may explain the characteristic treatment resistance of this condition.