Haplotype-specific assembly of shattered chromosomes in esophageal adenocarcinomas.

The epigenetic landscape of cancer is regulated by many factors, but primarily it derives from the underlying genome sequence. Chromothripsis is a catastrophic localized genome shattering event that drives, and often initiates, cancer evolution. We characterized five esophageal adenocarcinoma organoids with chromothripsis using long-read sequencing and transcriptome and epigenome profiling. Complex structural variation and subclonal variants meant that haplotype-aware de novo methods were required to generate contiguous cancer genome assemblies. Chromosomes were assembled separately and scaffolded using haplotype-resolved Hi-C reads, producing accurate assemblies even with up to 900 structural rearrangements. There were widespread differences between the chromothriptic and wild-type copies of chromosomes in topologically associated domains, chromatin accessibility, histone modifications, and gene expression. Differential epigenome peaks were most enriched within 10 kb of chromothriptic structural variants. Alterations in transcriptome and higher-order chromosome organization frequently occurred near differential epigenetic marks. Overall, chromothripsis reshapes gene regulation, causing coordinated changes in epigenetic landscape, transcription, and chromosome conformation.

Serum testosterone and prostate cancer in men with germline BRCA1/2 pathogenic variants.

Objectives: The relation of serum androgens and the development of prostate cancer (PCa) is subject of debate. Lower total testosterone (TT) levels have been associated with increased PCa detection and worse pathological features after treatment. However, data from the Reduction by Dutasteride of Prostate Cancer Events (REDUCE) and Prostate Cancer Prevention (PCPT) trial groups indicate no association. The aim of this study is to investigate the association of serum androgen levels and PCa detection in a prospective screening study of men at higher genetic risk of aggressive PCa due to BRCA1/2 pathogenic variants (PVs), the IMPACT study. Methods: Men enrolled in the IMPACT study provided serum samples during regular visits. Hormonal levels were calculated using immunoassays. Free testosterone (FT) was calculated from TT and sex hormone binding globulin (SHBG) using the Sodergard mass equation. Age, body mass index (BMI), prostate-specific antigen (PSA) and hormonal concentrations were compared between genetic cohorts. We also explored associations between age and TT, SHBG, FT and PCa, in the whole subset and stratified by BRCA1/2 PVs status. Results: A total of 777 participants in the IMPACT study had TT and SHBG measurements in serum samples at annual visits, giving 3940 prospective androgen levels, from 266 BRCA1 PVs carriers, 313 BRCA2 PVs carriers and 198 non-carriers. The median number of visits per patient was 5. There was no difference in TT, SHBG and FT between carriers and non-carriers. In a univariate analysis, androgen levels were not associated with PCa. In the analysis stratified by carrier status, no significant association was found between hormonal levels and PCa in non-carriers, BRCA1 or BRCA2 PVs carriers. Conclusions: Male BRCA1/2 PVs carriers have a similar androgen profile to non-carriers. Hormonal levels were not associated with PCa in men with and without BRCA1/2 PVs. Mechanisms related to the particularly aggressive phenotype of PCa in BRCA2 PVs carriers may therefore not be linked with circulating hormonal levels.

Diagnostic utility of whole genome sequencing in adults with B-other acute lymphoblastic leukemia.

Genomic profiling during the diagnosis of B-cell precursor acute lymphoblastic leukemia (BCP-ALL) in adults is used to guide disease classification, risk stratification, and treatment decisions. Patients for whom diagnostic screening fails to identify disease-defining or risk-stratifying lesions are classified as having B-other ALL. We screened a cohort of 652 BCP-ALL cases enrolled in UKALL14 to identify and perform whole genome sequencing (WGS) of paired tumor-normal samples. For 52 patients with B-other, we compared the WGS findings with data from clinical and research cytogenetics. WGS identified a cancer-associated event in 51 of 52 patients, including an established subtype defining genetic alterations that were previously missed with standard-of-care (SoC) genetics in 5 of them. Of the 47 true B-other ALL, we identified a recurrent driver in 87% (41). A complex karyotype via cytogenetics emerges as a heterogeneous group, including distinct genetic alterations associated with either favorable (DUX4-r) or poor outcomes (MEF2D-r and IGK::BCL2). For a subset of 31 cases, we integrated the findings from RNA sequencing (RNA-seq) analysis to include fusion gene detection and classification based on gene expression. Compared with RNA-seq, WGS was sufficient to detect and resolve recurrent genetic subtypes; however, RNA-seq can provide orthogonal validation of findings. In conclusion, we demonstrated that WGS can identify clinically relevant genetic abnormalities missed with SoC testing as well as identify leukemia driver events in virtually all cases of B-other ALL.

Clonal diversification and histogenesis of malignant germ cell tumours.

Germ cell tumours (GCTs) are a collection of benign and malignant neoplasms derived from primordial germ cells. They are uniquely able to recapitulate embryonic and extraembryonic tissues, which carries prognostic and therapeutic significance. The developmental pathways underpinning GCT initiation and histogenesis are incompletely understood. Here, we study the relationship of histogenesis and clonal diversification in GCTs by analysing the genomes and transcriptomes of 547 microdissected histological units. We find no correlation between genomic and histological heterogeneity. However, we identify unifying features including the retention of fetal developmental transcripts across tissues, expression changes on chromosome 12p, and a conserved somatic evolutionary sequence of whole genome duplication followed by clonal diversification. While this pattern is preserved across all GCTs, the developmental timing of the duplication varies between prepubertal and postpubertal cases. In addition, tumours of younger children exhibit distinct substitution signatures which may lend themselves as potential biomarkers for risk stratification. Our findings portray the extensive diversification of GCT tissues and genetic subclones as randomly distributed, while identifying overarching transcriptional and genomic features.

Unified classification and risk-stratification in Acute Myeloid Leukemia.

Clinical recommendations for Acute Myeloid Leukemia (AML) classification and risk-stratification remain heavily reliant on cytogenetic findings at diagnosis, which are present in <50% of patients. Using comprehensive molecular profiling data from 3,653 patients we characterize and validate 16 molecular classes describing 100% of AML patients. Each class represents diverse biological AML subgroups, and is associated with distinct clinical presentation, likelihood of response to induction chemotherapy, risk of relapse and death over time. Secondary AML-2, emerges as the second largest class (24%), associates with high-risk disease, poor prognosis irrespective of flow Minimal Residual Disease (MRD) negativity, and derives significant benefit from transplantation. Guided by class membership we derive a 3-tier risk-stratification score that re-stratifies 26% of patients as compared to standard of care. This results in a unified framework for disease classification and risk-stratification in AML that relies on information from cytogenetics and 32 genes. Last, we develop an open-access patient-tailored clinical decision support tool.

Diverse mutational landscapes in human lymphocytes.

The lymphocyte genome is prone to many threats, including programmed mutation during differentiation1, antigen-driven proliferation and residency in diverse microenvironments. Here, after developing protocols for expansion of single-cell lymphocyte cultures, we sequenced whole genomes from 717 normal naive and memory B and T cells and haematopoietic stem cells. All lymphocyte subsets carried more point mutations and structural variants than haematopoietic stem cells, with higher burdens in memory cells than in naive cells, and with T cells accumulating mutations at a higher rate throughout life. Off-target effects of immunological diversification accounted for approximately half of the additional differentiation-associated mutations in lymphocytes. Memory B cells acquired, on average, 18 off-target mutations genome-wide for every on-target IGHV mutation during the germinal centre reaction. Structural variation was 16-fold higher in lymphocytes than in stem cells, with around 15% of deletions being attributable to off-target recombinase-activating gene activity. DNA damage from ultraviolet light exposure and other sporadic mutational processes generated hundreds to thousands of mutations in some memory cells. The mutation burden and signatures of normal B cells were broadly similar to those seen in many B-cell cancers, suggesting that malignant transformation of lymphocytes arises from the same mutational processes that are active across normal ontogeny. The mutational landscape of normal lymphocytes chronicles the off-target effects of programmed genome engineering during immunological diversification and the consequences of differentiation, proliferation and residency in diverse microenvironments.

Prognostic impact of chromosomal abnormalities and copy number alterations in adult B-cell precursor acute lymphoblastic leukaemia: a UKALL14 study.

Chromosomal abnormalities are established prognostic markers in adult ALL. We assessed the prognostic impact of established chromosomal abnormalities and key copy number alterations (CNA) among 652 patients with B-cell precursor ALL treated on a modern MRD driven protocol. Patients with KMT2A-AFF1, complex karyotype (CK) and low hypodiploidy/near-triploidy (HoTr) had high relapse rates 50%, 60% & 53% and correspondingly poor survival. Patients with BCR-ABL1 had an outcome similar to other patients. JAK-STAT abnormalities (CRLF2, JAK2) occurred in 6% patients and were associated with a high relapse rate (56%). Patients with ABL-class fusions were rare (1%). A small group of patients with ZNF384 fusions (n = 12) had very good survival. CNA affecting IKZF1, CDKN2A/B, PAX5, BTG1, ETV6, EBF1, RB1 and PAR1 were assessed in 436 patients. None of the individual deletions or profiles were associated with survival, either in the cohort overall or within key subgroups. Collectively these data indicate that primary genetic abnormalities are stronger prognostic markers than secondary deletions. We propose a revised UKALL genetic risk classification based on key established chromosomal abnormalities: (1) very high risk: CK, HoTr or JAK-STAT abnormalities; (2) high risk: KMT2A fusions; (3) Tyrosine kinase activating: BCR-ABL1 and ABL-class fusions; (4) standard risk: all other patients.

Genetic and genomic analysis of acute lymphoblastic leukemia in older adults reveals a distinct profile of abnormalities: analysis of 210 patients from the UKALL14 and UKALL60+ clinical trials.

Despite being predominantly a childhood disease, the incidence of acute lymphoblastic leukemia (ALL) has a second peak in adults aged 60 years and over. These older adults fare extremely poorly with existing treatment strategies and very few studies have undertaken a comprehensive genetic and genomic characterization to improve prognosis in this age group. We performed cytogenetic, single nucleotide polymorphism (SNP) array and next-generation sequencing (NGS) analyses on samples from 210 patients aged ≥60 years from the UKALL14 and UKALL60+ clinical trials. BCR-ABL1-positive disease was present in 26% (55/210) of patients, followed by low hypodiploidy/near triploidy in 13% (28/210). Cytogenetically cryptic rearrangements in CRLF2, ZNF384 and MEF2D were detected in 5%, 1% and <1% of patients, respectively. Copy number abnormalities were common and deletions in ALL driver genes were seen in 77% of cases. IKZF1 deletion was present in 51% (40/78) of samples tested and the IKZF1plus profile was identified in over a third (28/77) of cases of B-cell precursor ALL. The genetic good-risk abnormalities high hyperdiploidy (n=2), ETV6-RUNX1 (no cases) and ERG deletion (no cases) were exceptionally rare in this cohort. RAS pathway mutations were seen in 17% (4/23) of screened samples. KDM6A abnormalities, including biallelic deletions, were discovered in 5% (4/78) of SNP arrays and 9% (2/23) of NGS samples, and represent novel, potentially therapeutically actionable lesions using EZH2 inhibitors. Outcome remained poor with 5-year event-free and overall survival rates of 17% and 24%, respectively, across the cohort, indicating a need for novel therapeutic strategies.

Convergent somatic mutations in metabolism genes in chronic liver disease.

The progression of chronic liver disease to hepatocellular carcinoma is caused by the acquisition of somatic mutations that affect 20-30 cancer genes1-8. Burdens of somatic mutations are higher and clonal expansions larger in chronic liver disease9-13 than in normal liver13-16, which enables positive selection to shape the genomic landscape9-13. Here we analysed somatic mutations from 1,590 genomes across 34 liver samples, including healthy controls, alcohol-related liver disease and non-alcoholic fatty liver disease. Seven of the 29 patients with liver disease had mutations in FOXO1, the major transcription factor in insulin signalling. These mutations affected a single hotspot within the gene, impairing the insulin-mediated nuclear export of FOXO1. Notably, six of the seven patients with FOXO1S22W hotspot mutations showed convergent evolution, with variants acquired independently by up to nine distinct hepatocyte clones per patient. CIDEB, which regulates lipid droplet metabolism in hepatocytes17-19, and GPAM, which produces storage triacylglycerol from free fatty acids20,21, also had a significant excess of mutations. We again observed frequent convergent evolution: up to fourteen independent clones per patient with CIDEB mutations and up to seven clones per patient with GPAM mutations. Mutations in metabolism genes were distributed across multiple anatomical segments of the liver, increased clone size and were seen in both alcohol-related liver disease and non-alcoholic fatty liver disease, but rarely in hepatocellular carcinoma. Master regulators of metabolic pathways are a frequent target of convergent somatic mutation in alcohol-related and non-alcoholic fatty liver disease.

IKZF1 alterations are not associated with outcome in 498 adults with B-precursor ALL enrolled in the UKALL14 trial.

IKZF1 deletions (ΔIKZF1) are commonly detected in B-precursor acute lymphoblastic leukemia (ALL; B-ALL) and are widely assumed to have a significant impact on outcome. We compared the ability of multiplex ligand-dependent probe amplification (MLPA) and polymerase chain reaction (PCR) to detect ΔIKZF1 and to determine the impact on event-free survival of patients with precursor B-ALL aged 23 to 65 years recruited to the completed trial UKALL14 (ISRCTN 66541317). From 655 recruits with BCR-ABL1+ and BCR-ABL1- B-ALL, all available diagnostic DNA samples (76% of the recruited population) were screened by multiplex end point PCR covering 4 deletions: dominant-negative (DN) Δ4-7 or the loss of function Δ2-7, Δ4-8, and Δ2-8 (n = 498), MLPA (n = 436), or by both (n = 420). Although patients with BCR-ABL1- ΔIKZF1 were more likely to have minimal residual disease at the end of induction, we did not find any impact of ΔIKZF1 (including subgroup analysis for DN or loss-of-function lesions) or the IKZF1plus genotype on event-free, overall survival, or relapse risk by univariable or multivariable analyses. Consistent with the technical approach, MLPA not only detected a wider range of deletions than PCR but also failed to detect some PCR-detected lesions. The main difference between our study and others reporting an association between ΔIKZF1 and outcome is the older age of participants in our population. The impact of ΔIKZF1 in ALL may be less marked in an older population of patients. Our study underscores the need for analyses in large, harmonized data sets. This trial was registered at www.clinicaltrials.gov as #NCT01085617.

The mutational landscape of human somatic and germline cells.

Over the course of an individual's lifetime, normal human cells accumulate mutations1. Here we compare the mutational landscape in 29 cell types from the soma and germline using multiple samples from the same individuals. Two ubiquitous mutational signatures, SBS1 and SBS5/40, accounted for the majority of acquired mutations in most cell types, but their absolute and relative contributions varied substantially. SBS18, which potentially reflects oxidative damage2, and several additional signatures attributed to exogenous and endogenous exposures contributed mutations to subsets of cell types. The rate of mutation was lowest in spermatogonia, the stem cells from which sperm are generated and from which most genetic variation in the human population is thought to originate. This was due to low rates of ubiquitous mutational processes and may be partially attributable to a low rate of cell division in basal spermatogonia. These results highlight similarities and differences in the maintenance of the germline and soma.

Whole-genome sequencing reveals progressive versus stable myeloma precursor conditions as two distinct entities.

Multiple myeloma (MM) is consistently preceded by precursor conditions recognized clinically as monoclonal gammopathy of undetermined significance (MGUS) or smoldering myeloma (SMM). We interrogate the whole genome sequence (WGS) profile of 18 MGUS and compare them with those from 14 SMMs and 80 MMs. We show that cases with a non-progressing, clinically stable myeloma precursor condition (n = 15) are characterized by later initiation in the patient's life and by the absence of myeloma defining genomic events including: chromothripsis, templated insertions, mutations in driver genes, aneuploidy, and canonical APOBEC mutational activity. This data provides evidence that WGS can be used to recognize two biologically and clinically distinct myeloma precursor entities that are either progressive or stable.

CDKN2A deletion is a frequent event associated with poor outcome in patients with peripheral T-cell lymphoma not otherwise specified (PTCL-NOS).

Nodal peripheral T-cell lymphoma not otherwise specified (PTCL-NOS) remains a diagnosis encompassing a heterogenous group of PTCL cases not fitting criteria for more homogeneous subtypes. They are characterized by a poor clinical outcome when treated with anthracycline-containing regimens. A better understanding of their biology could improve prognostic stratification and foster the development of novel therapeutic approaches. Recent targeted and whole exome sequencing studies have shown recurrent copy number abnormalities (CNAs) with prognostic significance. Here, investigating 5 formalin-fixed, paraffin embedded cases of PTCL-NOS by whole genome sequencing (WGS), we found a high prevalence of structural variants and complex events, such as chromothripsis likely responsible for the observed CNAs. Among them, CDKN2A and PTEN deletions emerged as the most frequent aberration, as confirmed in a final cohort of 143 patients with nodal PTCL. The incidence of CDKN2A and PTEN deletions among PTCL-NOS was 46% and 26%, respectively. Furthermore, we found that co-occurrence of CDKN2A and PTEN deletions is an event associated with PTCL-NOS with absolute specificity. In contrast, these deletions were rare and never co-occurred in angioimmunoblastic and anaplastic lymphomas. CDKN2A deletion was associated with shorter overall survival in multivariate analysis corrected by age, IPI, transplant eligibility and GATA3 expression (adjusted HR =2.53; 95% CI 1.006-6.3; p=0.048). These data suggest that CDKN2A deletions may be relevant for refining the prognosis of PTCL-NOS and their significance should be evaluated in prospective trials.

The mutational landscape of normal human endometrial epithelium.

All normal somatic cells are thought to acquire mutations, but understanding of the rates, patterns, causes and consequences of somatic mutations in normal cells is limited. The uterine endometrium adopts multiple physiological states over a lifetime and is lined by a gland-forming epithelium1,2. Here, using whole-genome sequencing, we show that normal human endometrial glands are clonal cell populations with total mutation burdens that increase at about 29 base substitutions per year and that are many-fold lower than those of endometrial cancers. Normal endometrial glands frequently carry 'driver' mutations in cancer genes, the burden of which increases with age and decreases with parity. Cell clones with drivers often originate during the first decades of life and subsequently progressively colonize the epithelial lining of the endometrium. Our results show that mutational landscapes differ markedly between normal tissues-perhaps shaped by differences in their structure and physiology-and indicate that the procession of neoplastic change that leads to endometrial cancer is initiated early in life.

Timing the initiation of multiple myeloma.

The evolution and progression of multiple myeloma and its precursors over time is poorly understood. Here, we investigate the landscape and timing of mutational processes shaping multiple myeloma evolution in a large cohort of 89 whole genomes and 973 exomes. We identify eight processes, including a mutational signature caused by exposure to melphalan. Reconstructing the chronological activity of each mutational signature, we estimate that the initial transformation of a germinal center B-cell usually occurred during the first 2nd-3rd decades of life. We define four main patterns of activation-induced deaminase (AID) and apolipoprotein B mRNA editing catalytic polypeptide-like (APOBEC) mutagenesis over time, including a subset of patients with evidence of prolonged AID activity during the pre-malignant phase, indicating antigen-responsiveness and germinal center reentry. Our findings provide a framework to study the etiology of multiple myeloma and explore strategies for prevention and early detection.

Type 1 Interferon Responses Underlie Tumor-Selective Replication of Oncolytic Measles Virus.

The mechanism of tumor-selective replication of oncolytic measles virus (MV) is poorly understood. Using a stepwise model of cellular transformation, in which oncogenic hits were additively expressed in human bone marrow-derived mesenchymal stromal cells, we show that MV-induced oncolysis increased progressively with transformation. The type 1 interferon (IFN) response to MV infection was significantly reduced and delayed, in accordance with the level of transformation. Consistently, we observed delayed and reduced signal transducer and activator of transcription (STAT1) phosphorylation in the fully transformed cells. Pre-treatment with IFNß restored resistance to MV-mediated oncolysis. Gene expression profiling to identify the genetic correlates of susceptibility to MV oncolysis revealed a dampened basal level of immune-related genes in the fully transformed cells compared to their normal counterparts. IFN-induced transmembrane protein 1 (IFITM1) was the foremost basally downregulated immune gene. Stable IFITM1 overexpression in MV-susceptible cells resulted in a 50% increase in cell viability and a significant reduction in viral replication at 24 h after MV infection. Overall, our data indicate that the basal reduction in functions of the type 1 IFN pathway is a major contributor to the oncolytic selectivity of MV. In particular, we have identified IFITM1 as a restriction factor for oncolytic MV, acting at early stages of infection.

Revealing the impact of structural variants in multiple myeloma.

The landscape of structural variants (SVs) in multiple myeloma remains poorly understood. Here, we performed comprehensive analysis of SVs in a large cohort of 752 multiple myeloma patients by low coverage long-insert whole genome sequencing. We identified 68 SV hotspots involving 17 new candidate driver genes, including the therapeutic targets BCMA (TNFRSF17), SLAMF and MCL1. Catastrophic complex rearrangements termed chromothripsis were present in 24% of patients and independently associated with poor clinical outcomes. Templated insertions were the second most frequent complex event (19%), mostly involved in super-enhancer hijacking and activation of oncogenes such as CCND1 and MYC. Importantly, in 31% of patients two or more seemingly independent putative driver events were caused by a single structural event, demonstrating that the complex genomic landscape of multiple myeloma can be acquired through few key events during tumor evolutionary history. Overall, this study reveals the critical role of SVs in multiple myeloma pathogenesis.

Genomic landscape and chronological reconstruction of driver events in multiple myeloma.

The multiple myeloma (MM) genome is heterogeneous and evolves through preclinical and post-diagnosis phases. Here we report a catalog and hierarchy of driver lesions using sequences from 67 MM genomes serially collected from 30 patients together with public exome datasets. Bayesian clustering defines at least 7 genomic subgroups with distinct sets of co-operating events. Focusing on whole genome sequencing data, complex structural events emerge as major drivers, including chromothripsis and a novel replication-based mechanism of templated insertions, which typically occur early. Hyperdiploidy also occurs early, with individual trisomies often acquired in different chronological windows during evolution, and with a preferred order of acquisition. Conversely, positively selected point mutations, whole genome duplication and chromoplexy events occur in later disease phases. Thus, initiating driver events, drawn from a limited repertoire of structural and numerical chromosomal changes, shape preferred trajectories of evolution that are biologically relevant but heterogeneous across patients.

A practical guide for mutational signature analysis in hematological malignancies.

Analysis of mutational signatures is becoming routine in cancer genomics, with implications for pathogenesis, classification, prognosis, and even treatment decisions. However, the field lacks a consensus on analysis and result interpretation. Using whole-genome sequencing of multiple myeloma (MM), chronic lymphocytic leukemia (CLL) and acute myeloid leukemia, we compare the performance of public signature analysis tools. We describe caveats and pitfalls of de novo signature extraction and fitting approaches, reporting on common inaccuracies: erroneous signature assignment, identification of localized hyper-mutational processes, overcalling of signatures. We provide reproducible solutions to solve these issues and use orthogonal approaches to validate our results. We show how a comprehensive mutational signature analysis may provide relevant biological insights, reporting evidence of c-AID activity among unmutated CLL cases or the absence of BRCA1/BRCA2-mediated homologous recombination deficiency in a MM cohort. Finally, we propose a general analysis framework to ensure production of accurate and reproducible mutational signature data.

Author Correction: A practical guide for mutational signature analysis in hematological malignancies.

An amendment to this paper has been published and can be accessed via a link at the top of the paper.