cgpPindel: Identifying Somatically Acquired Insertion and Deletion Events from Paired End Sequencing

ONCOLINER: A new solution for monitoring, improving, and harmonizing somatic variant calling across genomic oncology centers

The characterization of somatic genomic variation associated with the biology of tumors is fundamental for cancer research and personalized medicine, as it guides the reliability and impact of cancer studies and genomic-based decisions in clinical oncology. However, the quality and scope of tumor genome analysis across cancer research centers and hospitals are currently highly heterogeneous, limiting the consistency of tumor diagnoses across hospitals and the possibilities of data sharing and data integration across studies. With the aim of providing users with actionable and personalized recommendations for the overall enhancement and harmonization of somatic variant identification across research and clinical environments, we have developed ONCOLINER. Using specifically designed mosaic and tumorized genomes for the analysis of recall and precision across somatic SNVs, insertions or deletions (indels), and structural variants (SVs), we demonstrate that ONCOLINER is capable of improving and harmonizing genome analysis across three state-of-the-art variant discovery pipelines in genomic oncology.

Germline variants in patients diagnosed with pediatric soft tissue sarcoma

BACKGROUND

While soft tissue sarcomas affect younger patients, few studies have assessed the distribution of underlying pathogenic germline variants.

PATIENTS AND METHODS

We retrospectively identified all pediatric and young adult patients (0-22 years) at Haukeland University Hospital, Norway (1981-2019), through clinical and pathological records. We identified n = 46 eligible patients. From these 46 patients, adequate material representing normal tissue was available for n = 41 cases (n = 24 diagnosed with rhabdomyosarcoma, 9 with synovial sarcomas, 2 with Ewing sarcomas, and 6 without further classification), with matching tumor tissue for n = 40. Normal tissue samples were analyzed for germline pathogenic variants (PVs) by targeted sequencing of 360 cancer genes.

RESULTS

Out of the 41 analyzed cases, we found PVs or likely PVs in 7 (17%). These variants were found in TP53, MUTYH, FANCC, DICER1, FANCA, MYO3A, and MYO5B. Supporting the causality of these PVs, four cases revealed loss of heterozygosity (LOH) of the wild-type allele in the tumor tissue, one patient with a PV in DICER1 had a second somatic variant in DICER1, and a patient with a PV in TP53 had the altered allele amplified in the tumor. For three out of five with available family history, a history of other cancers in relatives was recorded. Among genes with variants of uncertain significance, CHD1L was of particular interest, revealing a stop-gain and a missense variant.

INTERPRETATION

A high fraction of young patients with soft tissue sarcoma harbor PVs. Among the genes affected, we substantiate a potential role of MYO5B and propose a potential role for MYO3A.

Reconstructing phylogenetic trees from genome-wide somatic mutations in clonal samples

Phylogenetic trees are a powerful means to display the evolutionary history of species, pathogens and, more recently, individual cells of the human body. Whole-genome sequencing of laser capture microdissections or expanded stem cells has allowed the discovery of somatic mutations in clones, which can be used as natural barcodes to reconstruct the developmental history of individual cells. Here we describe Sequoia, our pipeline to reconstruct lineage trees from clones of normal cells. Candidate somatic mutations are called against the human reference genome and filtered to exclude germline mutations and artifactual variants. These filtered somatic mutations form the basis for phylogeny reconstruction using a maximum parsimony framework. Lastly, we use a maximum likelihood framework to explicitly map mutations to branches in the phylogenetic tree. The resulting phylogenies can then serve as a basis for many subsequent analyses, including investigating embryonic development, tissue dynamics in health and disease, and mutational signatures. Sequoia can be readily applied to any clonal somatic mutation dataset, including single-cell DNA sequencing datasets, using the commands and scripts provided. Moreover, Sequoia is highly flexible and can be easily customized. Typically, the runtime of the core script ranges from minutes to an hour for datasets with a moderate number (50,000-150,000) of variants. Competent bioinformatic skills, including in-depth knowledge of the R programming language, are required. A high-performance computing cluster (one that is capable of running mutation-calling algorithms and other aspects of the analysis at scale) is also required, especially if handling large datasets.

Genomic characterization and risk stratification of esophageal squamous dysplasia

Objectives

The majority of esophageal squamous dysplasia (ESD) patients progress slowly, while a subset of patients can undergo recurrence rapidly or progress to invasive cancer even after proper treatment. However, the molecular mechanisms underlying these clinical observations are still largely unknown.

Methods

By sequencing the genomic data of 160 clinical samples from 49 tumor-free ESD patients and 88 esophageal squamous cell carcinoma (ESCC) patients, we demonstrated lower somatic mutation and copy number alteration (CNA) burden in ESD compared with ESCC.

Results

Cross-species screening and functional assays identified ACSM5 as a novel driver gene for ESD progression. Furthermore, we revealed that miR-4292 promoted ESD progression and could serve as a non-invasive diagnostic marker for ESD.

Conclusions

These findings largely expanded our understanding of ESD genetics and tumorigenesis, which possessed promising significance for improving early diagnosis, reducing overtreatment, and identifying high-risk ESD patients.

Geographic variation of mutagenic exposures in kidney cancer genomes

International differences in the incidence of many cancer types indicate the existence of carcinogen exposures that have not yet been identified by conventional epidemiology make a substantial contribution to cancer burden1. In clear cell renal cell carcinoma, obesity, hypertension and tobacco smoking are risk factors, but they do not explain the geographical variation in its incidence2. Underlying causes can be inferred by sequencing the genomes of cancers from populations with different incidence rates and detecting differences in patterns of somatic mutations. Here we sequenced 962 clear cell renal cell carcinomas from 11 countries with varying incidence. The somatic mutation profiles differed between countries. In Romania, Serbia and Thailand, mutational signatures characteristic of aristolochic acid compounds were present in most cases, but these were rare elsewhere. In Japan, a mutational signature of unknown cause was found in more than 70% of cases but in less than 2% elsewhere. A further mutational signature of unknown cause was ubiquitous but exhibited higher mutation loads in countries with higher incidence rates of kidney cancer. Known signatures of tobacco smoking correlated with tobacco consumption, but no signature was associated with obesity or hypertension, suggesting that non-mutagenic mechanisms of action underlie these risk factors. The results of this study indicate the existence of multiple, geographically variable, mutagenic exposures that potentially affect tens of millions of people and illustrate the opportunities for new insights into cancer causation through large-scale global cancer genomics.

The chemotherapeutic drug CX-5461 is a potent mutagen in cultured human cells

The chemotherapeutic agent CX-5461, or pidnarulex, has been fast-tracked by the United States Food and Drug Administration for early-stage clinical studies of BRCA1-, BRCA2- and PALB2-mutated cancers. It is under investigation in phase I and II trials. Here, we find that, although CX-5461 exhibits synthetic lethality in BRCA1-/BRCA2-deficient cells, it also causes extensive, nonselective, collateral mutagenesis in all three cell lines tested, to magnitudes that exceed known environmental carcinogens.

Intestinal-Type Adenocarcinoma in Head and Neck: Dissecting Oncogenic Gene Alterations Through Whole Transcriptome and Exome Analysis

Adenocarcinomas of the nasal/paranasal sinuses are uncommon, but intestinal-type adenocarcinomas (ITACs) are important. Due to the rarity of these tumors, their molecular profile is not well known. To further investigate the molecular profile and find potential oncogenic drivers, we compared the whole transcriptome and exome of ITACs at different anatomic locations in the head and neck. Twenty-one head and neck adenocarcinomas were used in this study, divided into 10 sinonasal adenocarcinomas (SNT) and 11 extrasinonasal (T) head and neck adenocarcinomas according to anatomic location and histology. Tumor samples along with normal mucosa were microdissected from formalin-fixed, paraffin-embedded samples, and RNA and DNA were subjected to whole-transcriptome and -exome shotgun sequencing. Analysis of ITACs at sinonasal locations showed 410 subtype-specific differentially expressed (DE) genes and noncoding transcripts compared with the group of other anatomic locations, with 2909 subtype-specific DE genes. The groups shared 872 genes, with 17 highly different or opposing DE genes. Whole-exome mutation analysis revealed the gene MLL3 (KMT2C) to be exhibiting the most frequent loss-of-function mutations in all adenocarcinomas investigated. The results suggest that the head and neck ITACs investigated were mainly caused by loss-of-function mutations in MLL3 that disabled chromatin methylation and remodeling of all MLL3-targeted enhancers in the tumors. This changed the activity of multiple genes/gene clusters, supporting oncogenicity mostly via pathways of signaling, dedifferentiation, proliferation, migration, and immune and inflammatory deregulation, indicating a truly epigenetic event as the root cause for the heterogenous diversity of these enteric types of cancer. The data of this study form the basis for understanding cell fate determination and cellular homeostasis in the normal respiratory mucosa at different anatomic sites and show the contribution of different mucosal components to the etiology/molecular pathology of ITAC.

Clonal selection of hematopoietic stem cells after gene therapy for sickle cell disease

Gene therapy (GT) provides a potentially curative treatment option for patients with sickle cell disease (SCD); however, the occurrence of myeloid malignancies in GT clinical trials has prompted concern, with several postulated mechanisms. Here, we used whole-genome sequencing to track hematopoietic stem cells (HSCs) from six patients with SCD at pre- and post-GT time points to map the somatic mutation and clonal landscape of gene-modified and unmodified HSCs. Pre-GT, phylogenetic trees were highly polyclonal and mutation burdens per cell were elevated in some, but not all, patients. Post-GT, no clonal expansions were identified among gene-modified or unmodified cells; however, an increased frequency of potential driver mutations associated with myeloid neoplasms or clonal hematopoiesis (DNMT3A- and EZH2-mutated clones in particular) was observed in both genetically modified and unmodified cells, suggesting positive selection of mutant clones during GT. This work sheds light on HSC clonal dynamics and the mutational landscape after GT in SCD, highlighting the enhanced fitness of some HSCs harboring pre-existing driver mutations. Future studies should define the long-term fate of mutant clones, including any contribution to expansions associated with myeloid neoplasms.

Effects of psoriasis and psoralen exposure on the somatic mutation landscape of the skin

Somatic mutations are hypothesized to play a role in many non-neoplastic diseases. We performed whole-exome sequencing of 1,182 microbiopsies dissected from lesional and nonlesional epidermis from 111 patients with psoriasis to search for evidence that somatic mutations in keratinocytes may influence the disease process. Lesional skin remained highly polyclonal, showing no evidence of large-scale spread of clones carrying potentially pathogenic mutations. The mutation rate of keratinocytes was similarly only modestly affected by the disease. We found evidence of positive selection in previously reported driver genes NOTCH1, NOTCH2, TP53, FAT1 and PPM1D and also identified mutations in four genes (GXYLT1, CHEK2, ZFP36L2 and EEF1A1) that we hypothesize are selected for in squamous epithelium irrespective of disease status. Finally, we describe a mutational signature of psoralens-a class of chemicals previously found in some sunscreens and which are used as part of PUVA (psoralens and ultraviolet-A) photochemotherapy treatment for psoriasis.

Neoantigen identification: Technological advances and challenges

Neoantigens have emerged as promising targets for cutting-edge immunotherapies, such as cancer vaccines and adoptive cell therapy. These neoantigens are unique to tumors and arise exclusively from somatic mutations or non-genomic aberrations in tumor proteins. They encompass a wide range of alterations, including genomic mutations, post-transcriptomic variants, and viral oncoproteins. With the advancements in technology, the identification of immunogenic neoantigens has seen rapid progress, raising new opportunities for enhancing their clinical significance. Prediction of neoantigens necessitates the acquisition of high-quality samples and sequencing data, followed by mutation calling. Subsequently, the pipeline involves integrating various tools that can predict the expression, processing, binding, and recognition potential of neoantigens. However, the continuous improvement of computational tools is constrained by the availability of datasets which contain validated immunogenic neoantigens. This review article aims to provide a comprehensive summary of the current knowledge as well as limitations in neoantigen prediction and validation. Additionally, it delves into the origin and biological role of neoantigens, offering a deeper understanding of their significance in the field of cancer immunotherapy. This article thus seeks to contribute to the ongoing efforts to harness neoantigens as powerful weapons in the fight against cancer.

Fanconi anemia-isogenic head and neck cancer cell line pairs: A basic and translational science resource

Fanconi anemia (FA) is a heritable malformation, bone marrow failure and cancer predisposition syndrome that confers an exceptionally high risk of squamous carcinomas. These carcinomas originate in epithelia lining the mouth, proximal esophagus, vulva and anus: their origins are not understood, and no effective ways have been identified to prevent or delay their appearance. Many FA-associated carcinomas are also therapeutically challenging: they may be multi-focal and stage-advanced at diagnosis, and most individuals with FA cannot tolerate standard-of-care systemic therapies such as DNA cross-linking drugs or ionizing radiation due to constitutional DNA damage hypersensitivity. We developed the Fanconi Anemia Cancer Cell Line Resource (FA-CCLR) to foster new work on the origins, treatment and prevention of FA-associated carcinomas. The FA-CCLR consists of Fanconi-isogenic head and neck squamous cell carcinoma (HNSCC) cell line pairs generated from five individuals with FA-associated HNSCC, and five individuals with sporadic HNSCC. Sporadic, isogenic HNSCC cell line pairs were generated in parallel with FA patient-derived isogenic cell line pairs to provide comparable experimental material to use to identify cell and molecular phenotypes driven by germline or somatic loss of Fanconi pathway function, and the subset of these FA-dependent phenotypes that can be modified, complemented or suppressed. All 10 FANC-isogenic cell line pairs are available to academic, non-profit and industry investigators via the "Fanconi Anemia Research Materials" Resource and Repository at Oregon Health & Sciences University, Portland OR.

A transposase-derived gene required for human brain development

DNA transposable elements and transposase-derived genes are present in most living organisms, including vertebrates, but their function is largely unknown. PiggyBac Transposable Element Derived 5 (PGBD5) is an evolutionarily conserved vertebrate DNA transposase-derived gene with retained nuclease activity in cells. Vertebrate brain development is known to be associated with prominent neuronal cell death and DNA breaks, but their causes and functions are not well understood. Here, we show that PGBD5 contributes to normal brain development in mice and humans, where its deficiency causes disorder of intellectual disability, movement and seizures. In mice, Pgbd5 is required for the developmental induction of post-mitotic DNA breaks and recurrent somatic genome rearrangements in neurons. Together, these studies nominate PGBD5 as the long-hypothesized neuronal DNA nuclease required for brain function in mammals.

A living biobank of patient-derived ductal carcinoma in situ mouse-intraductal xenografts identifies risk factors for invasive progression

Ductal carcinoma in situ (DCIS) is a non-obligate precursor of invasive breast cancer (IBC). Due to a lack of biomarkers able to distinguish high- from low-risk cases, DCIS is treated similar to early IBC even though the minority of untreated cases eventually become invasive. Here, we characterized 115 patient-derived mouse-intraductal (MIND) DCIS models reflecting the full spectrum of DCIS observed in patients. Utilizing the possibility to follow the natural progression of DCIS combined with omics and imaging data, we reveal multiple prognostic factors for high-risk DCIS including high grade, HER2 amplification, expansive 3D growth, and high burden of copy number aberrations. In addition, sequential transplantation of xenografts showed minimal phenotypic and genotypic changes over time, indicating that invasive behavior is an intrinsic phenotype of DCIS and supporting a multiclonal evolution model. Moreover, this study provides a collection of 19 distributable DCIS-MIND models spanning all molecular subtypes.

A phase I prospective, non-randomized trial of autologous dendritic cell-based cryoimmunotherapy in patients with metastatic castration-resistant prostate cancer

Metastatic castration-resistant prostate cancer (mCRPC) is an immunologically cold disease with dismal outcomes. Cryoablation destroys cancer tissue, releases tumor-associated antigens and creates a pro-inflammatory microenvironment, while dendritic cells (DCs) activate immune responses through processing of antigens. Immunotherapy combinations could enhance the anti-tumor efficacy. This open-label, single-arm, single-center phase I trial determined the safety and tolerability of combining cryoablation and autologous immature DC, without and with checkpoint inhibitors. Immune responses and clinical outcomes were evaluated. Patients with mCRPC, confirmed metastases and intact prostate gland were included. The first participants underwent prostate cryoablation with intratumoral injection of autologous DCs in a 3 + 3 design. In the second part, patients received cryoablation, the highest acceptable DC dose, and checkpoint inhibition with either ipilimumab or pembrolizumab. Sequentially collected information on adverse events, quality of life, blood values and images were analyzed by standard descriptive statistics. Neither dose-limiting toxicities nor adverse events > grade 3 were observed in the 18 participants. Results indicate antitumor activity through altered T cell receptor repertoires, and 33% durable (> 46 weeks) clinical benefit with median 40.7 months overall survival. Post-treatment pain and fatigue were associated with circulating tumor cell (CTC) presence at inclusion, while CTC responses correlated with clinical outcomes. This trial demonstrates that cryoimmunotherapy in mCRPC is safe and well tolerated, also for the highest DC dose (2.0 × 10⁸) combined with checkpoint inhibitors. Further studies focusing on the biologic indications of antitumor activity and immune system activation could be considered through a phase II trial focusing on treatment responses and immunologic biomarkers.

Cis-regulatory mutations associate with transcriptional and post-transcriptional deregulation of gene regulatory programs in cancers

Most cancer alterations occur in the noncoding portion of the human genome, where regulatory regions control gene expression. The discovery of noncoding mutations altering the cells' regulatory programs has been limited to few examples with high recurrence or high functional impact. Here, we show that transcription factor binding sites (TFBSs) have similar mutation loads to those in protein-coding exons. By combining cancer somatic mutations in TFBSs and expression data for protein-coding and miRNA genes, we evaluate the combined effects of transcriptional and post-transcriptional alterations on the regulatory programs in cancers. The analysis of seven TCGA cohorts culminates with the identification of protein-coding and miRNA genes linked to mutations at TFBSs that are associated with a cascading trans-effect deregulation on the cells' regulatory programs. Our analyses of cis-regulatory mutations associated with miRNAs recurrently predict 12 mature miRNAs (derived from 7 precursors) associated with the deregulation of their target gene networks. The predictions are enriched for cancer-associated protein-coding and miRNA genes and highlight cis-regulatory mutations associated with the dysregulation of key pathways associated with carcinogenesis. By combining transcriptional and post-transcriptional regulation of gene expression, our method predicts cis-regulatory mutations related to the dysregulation of key gene regulatory networks in cancer patients.

Defining the origin, evolution, and immune composition of SDH-deficient renal cell carcinoma

Succinate dehydrogenase (SDH)-deficient renal cell carcinoma represents a rare subtype of hereditary kidney cancer. Clinical diagnosis can be challenging and there is little evidence to guide systemic therapeutic options. We performed genomic profiling of a cohort of tumors through the analysis of whole genomes, transcriptomes, as well as flow cytometry and immunohistochemistry in order to gain a deeper understanding of their molecular biology. We find neutral evolution after early tumor activation with a lack of secondary driver events. We show that these tumors have epithelial derivation, possibly from the macula densa, a specialized paracrine cell of the renal juxtaglomerular apparatus. They subsequently develop into immune excluded tumors. We provide transcriptomic and protein expression evidence of a highly specific tumor marker, PAPPA2. These translational findings have implications for the diagnosis and treatment for this rare tumor subtype.

p53 mutation in normal esophagus promotes multiple stages of carcinogenesis but is constrained by clonal competition

Aging normal human oesophagus accumulates TP53 mutant clones. These are the origin of most oesophageal squamous carcinomas, in which biallelic TP53 disruption is almost universal. However, how p53 mutant clones expand and contribute to cancer development is unclear. Here we show that inducing the p53R245W mutant in single oesophageal progenitor cells in transgenic mice confers a proliferative advantage and clonal expansion but does not disrupt normal epithelial structure. Loss of the remaining p53 allele in mutant cells results in genomically unstable p53R245W/null epithelium with giant polyaneuploid cells and copy number altered clones. In carcinogenesis, p53 mutation does not initiate tumour formation, but tumours developing from areas with p53 mutation and LOH are larger and show extensive chromosomal instability compared to lesions arising in wild type epithelium. We conclude that p53 has distinct functions at different stages of carcinogenesis and that LOH within p53 mutant clones in normal epithelium is a critical step in malignant transformation.

The architecture of clonal expansions in morphologically normal tissue from cancerous and non-cancerous prostates

BACKGROUND

Up to 80% of cases of prostate cancer present with multifocal independent tumour lesions leading to the concept of a field effect present in the normal prostate predisposing to cancer development. In the present study we applied Whole Genome DNA Sequencing (WGS) to a group of morphologically normal tissue (n = 51), including benign prostatic hyperplasia (BPH) and non-BPH samples, from men with and men without prostate cancer. We assess whether the observed genetic changes in morphologically normal tissue are linked to the development of cancer in the prostate.

RESULTS

Single nucleotide variants (P = 7.0 × 10-03, Wilcoxon rank sum test) and small insertions and deletions (indels, P = 8.7 × 10-06) were significantly higher in morphologically normal samples, including BPH, from men with prostate cancer compared to those without. The presence of subclonal expansions under selective pressure, supported by a high level of mutations, were significantly associated with samples from men with prostate cancer (P = 0.035, Fisher exact test). The clonal cell fraction of normal clones was always higher than the proportion of the prostate estimated as epithelial (P = 5.94 × 10-05, paired Wilcoxon signed rank test) which, along with analysis of primary fibroblasts prepared from BPH specimens, suggests a stromal origin. Constructed phylogenies revealed lineages associated with benign tissue that were completely distinct from adjacent tumour clones, but a common lineage between BPH and non-BPH morphologically normal tissues was often observed. Compared to tumours, normal samples have significantly less single nucleotide variants (P = 3.72 × 10-09, paired Wilcoxon signed rank test), have very few rearrangements and a complete lack of copy number alterations.

CONCLUSIONS

Cells within regions of morphologically normal tissue (both BPH and non-BPH) can expand under selective pressure by mechanisms that are distinct from those occurring in adjacent cancer, but that are allied to the presence of cancer. Expansions, which are probably stromal in origin, are characterised by lack of recurrent driver mutations, by almost complete absence of structural variants/copy number alterations, and mutational processes similar to malignant tissue. Our findings have implications for treatment (focal therapy) and early detection approaches.

Detection of early seeding of Richter transformation in chronic lymphocytic leukemia

Richter transformation (RT) is a paradigmatic evolution of chronic lymphocytic leukemia (CLL) into a very aggressive large B cell lymphoma conferring a dismal prognosis. The mechanisms driving RT remain largely unknown. We characterized the whole genome, epigenome and transcriptome, combined with single-cell DNA/RNA-sequencing analyses and functional experiments, of 19 cases of CLL developing RT. Studying 54 longitudinal samples covering up to 19 years of disease course, we uncovered minute subclones carrying genomic, immunogenetic and transcriptomic features of RT cells already at CLL diagnosis, which were dormant for up to 19 years before transformation. We also identified new driver alterations, discovered a new mutational signature (SBS-RT), recognized an oxidative phosphorylation (OXPHOS)^high–B cell receptor (BCR)^low-signaling transcriptional axis in RT and showed that OXPHOS inhibition reduces the proliferation of RT cells. These findings demonstrate the early seeding of subclones driving advanced stages of cancer evolution and uncover potential therapeutic targets for RT.

Single-cell genomic and transcriptomic analyses of longitudinal samples of patients with Richter syndrome reveal the presence and dynamics of clones driving transformation from chronic lymphocytic leukemia years before clinical manifestation

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.

Mechanisms of APOBEC3 mutagenesis in human cancer cells

The APOBEC3 family of cytosine deaminases has been implicated in some of the most prevalent mutational signatures in cancer^1-3. However, a causal link between endogenous APOBEC3 enzymes and mutational signatures in human cancer genomes has not been established, leaving the mechanisms of APOBEC3 mutagenesis poorly understood. Here, to investigate the mechanisms of APOBEC3 mutagenesis, we deleted implicated genes from human cancer cell lines that naturally generate APOBEC3-associated mutational signatures over time⁴. Analysis of non-clustered and clustered signatures across whole-genome sequences from 251 breast, bladder and lymphoma cancer cell line clones revealed that APOBEC3A deletion diminished APOBEC3-associated mutational signatures. Deletion of both APOBEC3A and APOBEC3B further decreased APOBEC3 mutation burdens, without eliminating them. Deletion of APOBEC3B increased APOBEC3A protein levels, activity and APOBEC3A-mediated mutagenesis in some cell lines. The uracil glycosylase UNG was required for APOBEC3-mediated transversions, whereas the loss of the translesion polymerase REV1 decreased overall mutation burdens. Together, these data represent direct evidence that endogenous APOBEC3 deaminases generate prevalent mutational signatures in human cancer cells. Our results identify APOBEC3A as the main driver of these mutations, indicate that APOBEC3B can restrain APOBEC3A-dependent mutagenesis while contributing its own smaller mutation burdens and dissect mechanisms that translate APOBEC3 activities into distinct mutational signatures.

The longitudinal dynamics and natural history of clonal haematopoiesis

Clonal expansions driven by somatic mutations become pervasive across human tissues with age, including in the haematopoietic system, where the phenomenon is termed clonal haematopoiesis1-4. The understanding of how and when clonal haematopoiesis develops, the factors that govern its behaviour, how it interacts with ageing and how these variables relate to malignant progression remains limited5,6. Here we track 697 clonal haematopoiesis clones from 385 individuals 55 years of age or older over a median of 13 years. We find that 92.4% of clones expanded at a stable exponential rate over the study period, with different mutations driving substantially different growth rates, ranging from 5% (DNMT3A and TP53) to more than 50% per year (SRSF2P95H). Growth rates of clones with the same mutation differed by approximately ±5% per year, proportionately affecting slow drivers more substantially. By combining our time-series data with phylogenetic analysis of 1,731 whole-genome sequences of haematopoietic colonies from 7 individuals from an older age group, we reveal distinct patterns of lifelong clonal behaviour. DNMT3A-mutant clones preferentially expanded early in life and displayed slower growth in old age, in the context of an increasingly competitive oligoclonal landscape. By contrast, splicing gene mutations drove expansion only later in life, whereas TET2-mutant clones emerged across all ages. Finally, we show that mutations driving faster clonal growth carry a higher risk of malignant progression. Our findings characterize the lifelong natural history of clonal haematopoiesis and give fundamental insights into the interactions between somatic mutation, ageing and clonal selection.

Clonal dynamics of haematopoiesis across the human lifespan

Age-related change in human haematopoiesis causes reduced regenerative capacity1, cytopenias2, immune dysfunction3 and increased risk of blood cancer4-6, but the reason for such abrupt functional decline after 70 years of age remains unclear. Here we sequenced 3,579 genomes from single cell-derived colonies of haematopoietic cells across 10 human subjects from 0 to 81 years of age. Haematopoietic stem cells or multipotent progenitors (HSC/MPPs) accumulated a mean of 17 mutations per year after birth and lost 30 base pairs per year of telomere length. Haematopoiesis in adults less than 65 years of age was massively polyclonal, with high clonal diversity and a stable population of 20,000-200,000 HSC/MPPs contributing evenly to blood production. By contrast, haematopoiesis in individuals aged over 75 showed profoundly decreased clonal diversity. In each of the older subjects, 30-60% of haematopoiesis was accounted for by 12-18 independent clones, each contributing 1-34% of blood production. Most clones had begun their expansion before the subject was 40 years old, but only 22% had known driver mutations. Genome-wide selection analysis estimated that between 1 in 34 and 1 in 12 non-synonymous mutations were drivers, accruing at constant rates throughout life, affecting more genes than identified in blood cancers. Loss of the Y chromosome conferred selective benefits in males. Simulations of haematopoiesis, with constant stem cell population size and constant acquisition of driver mutations conferring moderate fitness benefits, entirely explained the abrupt change in clonal structure in the elderly. Rapidly decreasing clonal diversity is a universal feature of haematopoiesis in aged humans, underpinned by pervasive positive selection acting on many more genes than currently identified.

Mutational landscape of normal epithelial cells in Lynch Syndrome patients

Lynch Syndrome (LS) is an autosomal dominant disease conferring a high risk of colorectal cancer due to germline heterozygous mutations in a DNA mismatch repair (MMR) gene. Although cancers in LS patients show elevated somatic mutation burdens, information on mutation rates in normal tissues and understanding of the trajectory from normal to cancer cell is limited. Here we whole genome sequence 152 crypts from normal and neoplastic epithelial tissues from 10 LS patients. In normal tissues the repertoire of mutational processes and mutation rates is similar to that found in wild type individuals. A morphologically normal colonic crypt with an increased mutation burden and MMR deficiency-associated mutational signatures is identified, which may represent a very early stage of LS pathogenesis. Phylogenetic trees of tumour crypts indicate that the most recent ancestor cell of each tumour is already MMR deficient and has experienced multiple cycles of clonal evolution. This study demonstrates the genomic stability of epithelial cells with heterozygous germline MMR gene mutations and highlights important differences in the pathogenesis of LS from other colorectal cancer predisposition syndromes.

Bayesian networks elucidate complex genomic landscapes in cancer

Bayesian networks (BNs) are disciplined, explainable Artificial Intelligence models that can describe structured joint probability spaces. In the context of understanding complex relations between a number of variables in biological settings, they can be constructed from observed data and can provide a guiding, graphical tool in exploring such relations. Here we propose BNs for elucidating the relations between driver events in large cancer genomic datasets. We present a methodology that is specifically tailored to biologists and clinicians as they are the main producers of such datasets. We achieve this by using an optimal BN learning algorithm based on well established likelihood functions and by utilising just two tuning parameters, both of which are easy to set and have intuitive readings. To enhance value to clinicians, we introduce (a) the use of heatmaps for families in each network, and (b) visualising pairwise co-occurrence statistics on the network. For binary data, an optional step of fitting logic gates can be employed. We show how our methodology enhances pairwise testing and how biologists and clinicians can use BNs for discussing the main relations among driver events in large genomic cohorts. We demonstrate the utility of our methodology by applying it to 5 cancer datasets revealing complex genomic landscapes. Our networks identify central patterns in all datasets including a central 4-way mutual exclusivity between HDR, t(4,14), t(11,14) and t(14,16) in myeloma, and a 3-way mutual exclusivity of three major players: CALR, JAK2 and MPL, in myeloproliferative neoplasms. These analyses demonstrate that our methodology can play a central role in the study of large genomic cancer datasets.

Bayesian network inference on several blood and solid cancer genomic datasets provides more accessible ways to explore driver events in cancer.

A suspension technique for efficient large-scale cancer organoid culturing and perturbation screens

Organoid cell culture methodologies are enabling the generation of cell models from healthy and diseased tissue. Patient-derived cancer organoids that recapitulate the genetic and histopathological diversity of patient tumours are being systematically generated, providing an opportunity to investigate new cancer biology and therapeutic approaches. The use of organoid cultures for many applications, including genetic and chemical perturbation screens, is limited due to the technical demands and cost associated with their handling and propagation. Here we report and benchmark a suspension culture technique for cancer organoids which allows for the expansion of models to tens of millions of cells with increased efficiency in comparison to standard organoid culturing protocols. Using whole-genome DNA and RNA sequencing analyses, as well as medium-throughput drug sensitivity testing and genome-wide CRISPR-Cas9 screening, we demonstrate that cancer organoids grown as a suspension culture are genetically and phenotypically similar to their counterparts grown in standard conditions. This culture technique simplifies organoid cell culture and extends the range of organoid applications, including for routine use in large-scale perturbation screens.

Identification of neoantigens for individualized therapeutic cancer vaccines

Somatic mutations in cancer cells can generate tumour-specific neoepitopes, which are recognized by autologous T cells in the host. As neoepitopes are not subject to central immune tolerance and are not expressed in healthy tissues, they are attractive targets for therapeutic cancer vaccines. Because the vast majority of cancer mutations are unique to the individual patient, harnessing the full potential of this rich source of targets requires individualized treatment approaches. Many computational algorithms and machine-learning tools have been developed to identify mutations in sequence data, to prioritize those that are more likely to be recognized by T cells and to design tailored vaccines for every patient. In this Review, we fill the gaps between the understanding of basic mechanisms of T cell recognition of neoantigens and the computational approaches for discovery of somatic mutations and neoantigen prediction for cancer immunotherapy. We present a new classification of neoantigens, distinguishing between guarding, restrained and ignored neoantigens, based on how they confer proficient antitumour immunity in a given clinical context. Such context-based differentiation will contribute to a framework that connects neoantigen biology to the clinical setting and medical peculiarities of cancer, and will enable future neoantigen-based therapies to provide greater clinical benefit.

Somatic mutation rates scale with lifespan across mammals

The rates and patterns of somatic mutation in normal tissues are largely unknown outside of humans1-7. Comparative analyses can shed light on the diversity of mutagenesis across species, and on long-standing hypotheses about the evolution of somatic mutation rates and their role in cancer and ageing. Here we performed whole-genome sequencing of 208 intestinal crypts from 56 individuals to study the landscape of somatic mutation across 16 mammalian species. We found that somatic mutagenesis was dominated by seemingly endogenous mutational processes in all species, including 5-methylcytosine deamination and oxidative damage. With some differences, mutational signatures in other species resembled those described in humans8, although the relative contribution of each signature varied across species. Notably, the somatic mutation rate per year varied greatly across species and exhibited a strong inverse relationship with species lifespan, with no other life-history trait studied showing a comparable association. Despite widely different life histories among the species we examined-including variation of around 30-fold in lifespan and around 40,000-fold in body mass-the somatic mutation burden at the end of lifespan varied only by a factor of around 3. These data unveil common mutational processes across mammals, and suggest that somatic mutation rates are evolutionarily constrained and may be a contributing factor in ageing.

Functional patient-derived organoid screenings identify MCLA-158 as a therapeutic EGFR × LGR5 bispecific antibody with efficacy in epithelial tumors

Patient-derived organoids (PDOs) recapitulate tumor architecture, contain cancer stem cells and have predictive value supporting personalized medicine. Here we describe a large-scale functional screen of dual-targeting bispecific antibodies (bAbs) on a heterogeneous colorectal cancer PDO biobank and paired healthy colonic mucosa samples. More than 500 therapeutic bAbs generated against Wingless-related integration site (WNT) and receptor tyrosine kinase (RTK) targets were functionally evaluated by high-content imaging to capture the complexity of PDO responses. Our drug discovery strategy resulted in the generation of MCLA-158, a bAb that specifically triggers epidermal growth factor receptor degradation in leucine-rich repeat-containing G-protein-coupled receptor 5-positive (LGR5+) cancer stem cells but shows minimal toxicity toward healthy LGR5+ colon stem cells. MCLA-158 exhibits therapeutic properties such as growth inhibition of KRAS-mutant colorectal cancers, blockade of metastasis initiation and suppression of tumor outgrowth in preclinical models for several epithelial cancer types.

The molecular characteristics of high-grade gastroenteropancreatic neuroendocrine neoplasms

High-grade (HG) gastroenteropancreatic (GEP) neuroendocrine neoplasms (NEN) are rare but have a very poor prognosis and represent a severely understudied class of tumours. Molecular data for HG GEP-NEN are limited, and treatment strategies for the carcinoma subgroup (HG GEP-NEC) are extrapolated from small-cell lung cancer (SCLC). After pathological re-evaluation, we analysed DNA from tumours and matched blood samples from 181 HG GEP-NEN patients; 152 neuroendocrine carcinomas (NEC) and 29 neuroendocrine tumours (NET G3). Based on the sequencing of 360 cancer-related genes, we assessed mutations and copy number alterations (CNA). For NEC, frequently mutated genes were TP53 (64%), APC (28%), KRAS (22%) and BRAF (20%). RB1 was only mutated in 14%, but CNAs affecting RB1 were seen in 34%. Other frequent copy number losses were ARID1A (35%), ESR1 (25%) and ATM (31%). Frequent amplifications/gains were found in MYC (51%) and KDM5A (45%). While these molecular features had limited similarities with SCLC, we found potentially targetable alterations in 66% of the NEC samples. Mutations and CNA varied according to primary tumour site with BRAF mutations mainly seen in colon (49%), and FBXW7 mutations mainly seen in rectal cancers (25%). Eight out of 152 (5.3%) NEC were microsatellite instable (MSI). NET G3 had frequent mutations in MEN1 (21%), ATRX (17%), DAXX, SETD2 and TP53 (each 14%). We show molecular differences in HG GEP-NEN, related to morphological differentiation and site of origin. Limited similarities to SCLC and a high fraction of targetable alterations indicate a high potential for better-personalized treatments.

Patterns of genomic change in residual disease after neoadjuvant chemotherapy for estrogen receptor-positive and HER2-negative breast cancer

BACKGROUND

Treatment of patients with residual disease after neoadjuvant chemotherapy for breast cancer is an unmet clinical need. We hypothesised that tumour subclones showing expansion in residual disease after chemotherapy would contain mutations conferring drug resistance.

METHODS

We studied oestrogen receptor and/or progesterone receptor-positive, HER2-negative tumours from 42 patients in the EORTC 10994/BIG 00-01 trial who failed to achieve a pathological complete response. Genes commonly mutated in breast cancer were sequenced in pre and post-treatment samples.

RESULTS

Oncogenic driver mutations were commonest in PIK3CA (38% of tumours), GATA3 (29%), CDH1 (17%), TP53 (17%) and CBFB (12%); and amplification was commonest for CCND1 (26% of tumours) and FGFR1 (26%). The variant allele fraction frequently changed after treatment, indicating that subclones had expanded and contracted, but there were changes in both directions for all of the commonly mutated genes.

CONCLUSIONS

We found no evidence that expansion of clones containing recurrent oncogenic driver mutations is responsible for resistance to neoadjuvant chemotherapy. The persistence of classic oncogenic mutations in pathways for which targeted therapies are now available highlights their importance as drug targets in patients who have failed chemotherapy but provides no support for a direct role of driver oncogenes in resistance to chemotherapy. CLINICALTRIALS.GOV: EORTC 10994/BIG 1-00 Trial registration number NCT00017095.

Mutational signatures in esophageal squamous cell carcinoma from eight countries with varying incidence

Esophageal squamous cell carcinoma (ESCC) shows remarkable variation in incidence that is not fully explained by known lifestyle and environmental risk factors. It has been speculated that an unknown exogenous exposure(s) could be responsible. Here we combine the fields of mutational signature analysis with cancer epidemiology to study 552 ESCC genomes from eight countries with varying incidence rates. Mutational profiles were similar across all countries studied. Associations between specific mutational signatures and ESCC risk factors were identified for tobacco, alcohol, opium and germline variants, with modest impacts on mutation burden. We find no evidence of a mutational signature indicative of an exogenous exposure capable of explaining differences in ESCC incidence. Apolipoprotein B mRNA-editing enzyme, catalytic polypeptide-like (APOBEC)-associated mutational signatures single-base substitution (SBS)2 and SBS13 were present in 88% and 91% of cases, respectively, and accounted for 25% of the mutation burden on average, indicating that APOBEC activation is a crucial step in ESCC tumor development.

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.

Increased somatic mutation burdens in normal human cells due to defective DNA polymerases

Mutation accumulation in somatic cells contributes to cancer development and is proposed as a cause of aging. DNA polymerases Pol ε and Pol δ replicate DNA during cell division. However, in some cancers, defective proofreading due to acquired POLE/POLD1 exonuclease domain mutations causes markedly elevated somatic mutation burdens with distinctive mutational signatures. Germline POLE/POLD1 mutations cause familial cancer predisposition. Here, we sequenced normal tissue and tumor DNA from individuals with germline POLE/POLD1 mutations. Increased mutation burdens with characteristic mutational signatures were found in normal adult somatic cell types, during early embryogenesis and in sperm. Thus human physiology can tolerate ubiquitously elevated mutation burdens. Except for increased cancer risk, individuals with germline POLE/POLD1 mutations do not exhibit overt features of premature aging. These results do not support a model in which all features of aging are attributable to widespread cell malfunction directly resulting from somatic mutation burdens accrued during life.

CloneSig can jointly infer intra-tumor heterogeneity and mutational signature activity in bulk tumor sequencing data

Systematic DNA sequencing of cancer samples has highlighted the importance of two aspects of cancer genomics: intra-tumor heterogeneity (ITH) and mutational processes. These two aspects may not always be independent, as different mutational processes could be involved in different stages or regions of the tumor, but existing computational approaches to study them largely ignore this potential dependency. Here, we present CloneSig, a computational method to jointly infer ITH and mutational processes in a tumor from bulk-sequencing data. Extensive simulations show that CloneSig outperforms current methods for ITH inference and detection of mutational processes when the distribution of mutational signatures changes between clones. Applied to a large cohort of 8,951 tumors with whole-exome sequencing data from The Cancer Genome Atlas, and on a pan-cancer dataset of 2,632 whole-genome sequencing tumor samples from the Pan-Cancer Analysis of Whole Genomes initiative, CloneSig obtains results overall coherent with previous studies.

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.

Development, maturation, and maintenance of human prostate inferred from somatic mutations

Clonal dynamics and mutation burden in healthy human prostate epithelium are relevant to prostate cancer. We sequenced whole genomes from 409 microdissections of normal prostate epithelium across 8 donors, using phylogenetic reconstruction with spatial mapping in a 59-year-old man's prostate to reconstruct tissue dynamics across the lifespan. Somatic mutations accumulate steadily at ∼16 mutations/year/clone, with higher rates in peripheral than peri-urethral regions. The 24-30 independent glandular subunits are established as rudimentary ductal structures during fetal development by 5-10 embryonic cells each. Puberty induces formation of further side and terminal branches by local stem cells disseminated throughout the rudimentary ducts during development. During adult tissue maintenance, clonal expansions have limited geographic scope and minimal migration. Driver mutations are rare in aging prostate epithelium, but the one driver we did observe generated a sizable intraepithelial clonal expansion. Leveraging unbiased clock-like mutations, we define prostate stem cell dynamics through fetal development, puberty, and aging.

Lineage tracing of human development through somatic mutations

The ontogeny of the human haematopoietic system during fetal development has previously been characterized mainly through careful microscopic observations1. Here we reconstruct a phylogenetic tree of blood development using whole-genome sequencing of 511 single-cell-derived haematopoietic colonies from healthy human fetuses at 8 and 18 weeks after conception, coupled with deep targeted sequencing of tissues of known embryonic origin. We found that, in healthy fetuses, individual haematopoietic progenitors acquire tens of somatic mutations by 18 weeks after conception. We used these mutations as barcodes and timed the divergence of embryonic and extra-embryonic tissues during development, and estimated the number of blood antecedents at different stages of embryonic development. Our data support a hypoblast origin of the extra-embryonic mesoderm and primitive blood in humans.

Comprehensive Molecular Profiling of Desmoplastic Small Round Cell Tumor

Desmoplastic small round cell tumor (DSRCT) is characterized by the EWSR1-WT1 t(11;22) (p13:q12) translocation. Few additional putative drivers have been identified, and research has suffered from a lack of model systems. Next-generation sequencing (NGS) data from 68 matched tumor-normal samples, whole-genome sequencing data from 10 samples, transcriptomic and affymetrix array data, and a bank of DSRCT patient-derived xenograft (PDX) are presented. EWSR1-WT1 fusions were noted to be simple, balanced events. Recurrent mutations were uncommon, but were noted in TERT (3%), ARID1A (6%), HRAS (5%), and TP53 (3%), and recurrent loss of heterozygosity (LOH) at 11p, 11q, and 16q was identified in 18%, 22%, and 34% of samples, respectively. Comparison of tumor-normal matched versus unmatched analysis suggests overcalling of somatic mutations in prior publications of DSRCT NGS data. Alterations in fibroblast growth factor receptor 4 (FGFR4) were identified in 5 of 68 (7%) of tumor samples, whereas differential overexpression of FGFR4 was confirmed orthogonally using 2 platforms. PDX models harbored the pathognomic EWSR1-WT1 fusion and were highly representative of corresponding tumors. Our analyses confirm DSRCT as a genomically quiet cancer defined by the balanced translocation, t(11;22)(p13:q12), characterized by a paucity of secondary mutations but a significant number of copy number alterations. Against this genomically quiet background, recurrent activating alterations of FGFR4 stood out, and suggest that this receptor tyrosine kinase, also noted to be highly expressed in DSRCT, should be further investigated. Future studies of DSRCT biology and preclinical therapeutic strategies should benefit from the PDX models characterized in this study. IMPLICATIONS: These data describe the general quiescence of the desmoplastic small round cell tumor (DSRCT) genome, present the first available bank of DSRCT model systems, and nominate FGFR4 as a key receptor tyrosine kinase in DSRCT, based on high expression, recurrent amplification, and recurrent activating mutations.

A systematic CRISPR screen defines mutational mechanisms underpinning signatures caused by replication errors and endogenous DNA damage

Mutational signatures are imprints of pathophysiological processes arising through tumorigenesis. We generated isogenic CRISPR-Cas9 knockouts (Δ) of 43 genes in human induced pluripotent stem cells, cultured them in the absence of added DNA damage, and performed whole-genome sequencing of 173 subclones. ΔOGG1, ΔUNG, ΔEXO1, ΔRNF168, ΔMLH1, ΔMSH2, ΔMSH6, ΔPMS1, and ΔPMS2 produced marked mutational signatures indicative of being critical mitigators of endogenous DNA modifications. Detailed analyses revealed mutational mechanistic insights, including how 8-oxo-dG elimination is sequence-context-specific while uracil clearance is sequence-context-independent. Mismatch repair (MMR) deficiency signatures are engendered by oxidative damage (C>A transversions), differential misincorporation by replicative polymerases (T>C and C>T transitions), and we propose a 'reverse template slippage' model for T>A transversions. ΔMLH1, ΔMSH6, and ΔMSH2 signatures were similar to each other but distinct from ΔPMS2. Finally, we developed a classifier, MMRDetect, where application to 7,695 WGS cancers showed enhanced detection of MMR-deficient tumors, with implications for responsiveness to immunotherapies.

Somatic mutation landscapes at single-molecule resolution

Somatic mutations drive the development of cancer and may contribute to ageing and other diseases1,2. Despite their importance, the difficulty of detecting mutations that are only present in single cells or small clones has limited our knowledge of somatic mutagenesis to a minority of tissues. Here, to overcome these limitations, we developed nanorate sequencing (NanoSeq), a duplex sequencing protocol with error rates of less than five errors per billion base pairs in single DNA molecules from cell populations. This rate is two orders of magnitude lower than typical somatic mutation loads, enabling the study of somatic mutations in any tissue independently of clonality. We used this single-molecule sensitivity to study somatic mutations in non-dividing cells across several tissues, comparing stem cells to differentiated cells and studying mutagenesis in the absence of cell division. Differentiated cells in blood and colon displayed remarkably similar mutation loads and signatures to their corresponding stem cells, despite mature blood cells having undergone considerably more divisions. We then characterized the mutational landscape of post-mitotic neurons and polyclonal smooth muscle, confirming that neurons accumulate somatic mutations at a constant rate throughout life without cell division, with similar rates to mitotically active tissues. Together, our results suggest that mutational processes that are independent of cell division are important contributors to somatic mutagenesis. We anticipate that the ability to reliably detect mutations in single DNA molecules could transform our understanding of somatic mutagenesis and enable non-invasive studies on large-scale cohorts.

Selection of Oncogenic Mutant Clones in Normal Human Skin Varies with Body Site

Skin cancer risk varies substantially across the body, yet how this relates to the mutations found in normal skin is unknown. Here we mapped mutant clones in skin from high- and low-risk sites. The density of mutations varied by location. The prevalence of NOTCH1 and FAT1 mutations in forearm, trunk, and leg skin was similar to that in keratinocyte cancers. Most mutations were caused by ultraviolet light, but mutational signature analysis suggested differences in DNA-repair processes between sites. Eleven mutant genes were under positive selection, with TP53 preferentially selected in the head and FAT1 in the leg. Fine-scale mapping revealed 10% of clones had copy-number alterations. Analysis of hair follicles showed mutations in the upper follicle resembled adjacent skin, but the lower follicle was sparsely mutated. Normal skin is a dense patchwork of mutant clones arising from competitive selection that varies by location. SIGNIFICANCE: Mapping mutant clones across the body reveals normal skin is a dense patchwork of mutant cells. The variation in cancer risk between sites substantially exceeds that in mutant clone density. More generally, mutant genes cannot be assigned as cancer drivers until their prevalence in normal tissue is known.See related commentary by De Dominici and DeGregori, p. 227.This article is highlighted in the In This Issue feature, p. 211.

Proteogenomic Analysis of Salivary Adenoid Cystic Carcinomas Defines Molecular Subtypes and Identifies Therapeutic Targets

PURPOSE

Salivary gland adenoid cystic carcinoma (ACC) has heterogeneous clinical behavior. Currently, all patients are treated uniformly, and no standard-of-care systemic therapy exists for metastatic ACC. We conducted an integrated proteogenomic analyses of ACC tumors to identify dysregulated pathways and propose a classification with therapeutic implications.

EXPERIMENTAL DESIGN

RNA/DNA sequencing of 54 flash-frozen salivary ACCs and reverse phase protein array (RPPA) in 38 specimens were performed, with validation by Western blotting and/or IHC. Three independent ACC cohorts were used for validation.

RESULTS

Both unbiased RNA sequencing (RNA-seq) and RPPA analysis revealed two molecular subtypes: ACC-I (37%) and ACC-II (63%). ACC-I had strong upregulation of MYC, MYC target genes, and mRNA splicing, enrichment of NOTCH-activating mutations, and dramatically worse prognosis. ACC-II exhibited upregulation of TP63 and receptor tyrosine kinases (AXL, MET, and EGFR) and less aggressive clinical course. TP63 and MYC were sufficient to assign tumors to ACC subtypes, which was validated in one independent cohort by IHC and two additional independent cohorts by RNA-seq. Furthermore, IHC staining for MYC and P63 protein levels can be used to identify ACC subtypes, enabling rapid clinical deployment to guide therapeutic decisions. Our data suggest a model in which ACC-I is driven by MYC signaling through either NOTCH mutations or direct amplification, which in turn suppress P63 signaling observed in ACC-II, producing unique therapeutic vulnerabilities for each subtype.

CONCLUSIONS

Cooccurrence of multiple actionable protein/pathways alterations in each subtype indicates unique therapeutic vulnerabilities and opportunities for optimal combination therapy for this understudied and heterogeneous disease.

The mutational signature profile of known and suspected human carcinogens in mice

Epidemiological studies have identified many environmental agents that appear to significantly increase cancer risk in human populations. By analyzing tumor genomes from mice chronically exposed to 1 of 20 known or suspected human carcinogens, we reveal that most agents do not generate distinct mutational signatures or increase mutation burden, with most mutations, including driver mutations, resulting from tissue-specific endogenous processes. We identify signatures resulting from exposure to cobalt and vinylidene chloride and link distinct human signatures (SBS19 and SBS42) with 1,2,3-trichloropropane, a haloalkane and pollutant of drinking water, and find these and other signatures in human tumor genomes. We define the cross-species genomic landscape of tumors induced by an important compendium of agents with relevance to human health.

Extensive heterogeneity in somatic mutation and selection in the human bladder

The extent of somatic mutation and clonal selection in the human bladder remains unknown. We sequenced 2097 bladder microbiopsies from 20 individuals using targeted (n = 1914 microbiopsies), whole-exome (n = 655), and whole-genome (n = 88) sequencing. We found widespread positive selection in 17 genes. Chromatin remodeling genes were frequently mutated, whereas mutations were absent in several major bladder cancer genes. There was extensive interindividual variation in selection, with different driver genes dominating the clonal landscape across individuals. Mutational signatures were heterogeneous across clones and individuals, which suggests differential exposure to mutagens in the urine. Evidence of APOBEC mutagenesis was found in 22% of the microbiopsies. Sequencing multiple microbiopsies from five patients with bladder cancer enabled comparisons with cancer-free individuals and across histological features. This study reveals a rich landscape of mutational processes and selection in normal urothelium with large heterogeneity across clones and individuals.

Immune Surveillance in Clinical Regression of Preinvasive Squamous Cell Lung Cancer

Before squamous cell lung cancer develops, precancerous lesions can be found in the airways. From longitudinal monitoring, we know that only half of such lesions become cancer, whereas a third spontaneously regress. Although recent studies have described the presence of an active immune response in high-grade lesions, the mechanisms underpinning clinical regression of precancerous lesions remain unknown. Here, we show that host immune surveillance is strongly implicated in lesion regression. Using bronchoscopic biopsies from human subjects, we find that regressive carcinoma in situ lesions harbor more infiltrating immune cells than those that progress to cancer. Moreover, molecular profiling of these lesions identifies potential immune escape mechanisms specifically in those that progress to cancer: antigen presentation is impaired by genomic and epigenetic changes, CCL27-CCR10 signaling is upregulated, and the immunomodulator TNFSF9 is downregulated. Changes appear intrinsic to the carcinoma in situ lesions, as the adjacent stroma of progressive and regressive lesions are transcriptomically similar. SIGNIFICANCE: Immune evasion is a hallmark of cancer. For the first time, this study identifies mechanisms by which precancerous lesions evade immune detection during the earliest stages of carcinogenesis and forms a basis for new therapeutic strategies that treat or prevent early-stage lung cancer.See related commentary by Krysan et al., p. 1442.This article is highlighted in the In This Issue feature, p. 1426.

Accelerated single cell seeding in relapsed multiple myeloma

Multiple myeloma (MM) progression is characterized by the seeding of cancer cells in different anatomic sites. To characterize this evolutionary process, we interrogated, by whole genome sequencing, 25 samples collected at autopsy from 4 patients with relapsed MM and an additional set of 125 whole exomes collected from 51 patients. Mutational signatures analysis showed how cytotoxic agents introduce hundreds of unique mutations in each surviving cancer cell, detectable by bulk sequencing only in cases of clonal expansion of a single cancer cell bearing the mutational signature. Thus, a unique, single-cell genomic barcode can link chemotherapy exposure to a discrete time window in a patient′s life. We leveraged this concept to show that MM systemic seeding is accelerated at relapse and appears to be driven by the survival and subsequent expansion of a single myeloma cell following treatment with high-dose melphalan therapy and autologous stem cell transplant.

In multiple myeloma, disease progresses via seeding to different anatomic sites and clonal expansion. Here, utilising autopsy material, the authors show that systemic seeding accelerates at relapse following treatment.

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.

In situ CRISPR-Cas9 base editing for the development of genetically engineered mouse models of breast cancer

Genetically engineered mouse models (GEMMs) of cancer have proven to be of great value for basic and translational research. Although CRISPR-based gene disruption offers a fast-track approach for perturbing gene function and circumvents certain limitations of standard GEMM development, it does not provide a flexible platform for recapitulating clinically relevant missense mutations in vivo. To this end, we generated knock-in mice with Cre-conditional expression of a cytidine base editor and tested their utility for precise somatic engineering of missense mutations in key cancer drivers. Upon intraductal delivery of sgRNA-encoding vectors, we could install point mutations with high efficiency in one or multiple endogenous genes in situ and assess the effect of defined allelic variants on mammary tumorigenesis. While the system also produces bystander insertions and deletions that can stochastically be selected for when targeting a tumor suppressor gene, we could effectively recapitulate oncogenic nonsense mutations. We successfully applied this system in a model of triple-negative breast cancer, providing the proof of concept for extending this flexible somatic base editing platform to other tissues and tumor types.

Pan-cancer analysis of whole genomes

Cancer is driven by genetic change, and the advent of massively parallel sequencing has enabled systematic documentation of this variation at the whole-genome scale1-3. Here we report the integrative analysis of 2,658 whole-cancer genomes and their matching normal tissues across 38 tumour types from the Pan-Cancer Analysis of Whole Genomes (PCAWG) Consortium of the International Cancer Genome Consortium (ICGC) and The Cancer Genome Atlas (TCGA). We describe the generation of the PCAWG resource, facilitated by international data sharing using compute clouds. On average, cancer genomes contained 4-5 driver mutations when combining coding and non-coding genomic elements; however, in around 5% of cases no drivers were identified, suggesting that cancer driver discovery is not yet complete. Chromothripsis, in which many clustered structural variants arise in a single catastrophic event, is frequently an early event in tumour evolution; in acral melanoma, for example, these events precede most somatic point mutations and affect several cancer-associated genes simultaneously. Cancers with abnormal telomere maintenance often originate from tissues with low replicative activity and show several mechanisms of preventing telomere attrition to critical levels. Common and rare germline variants affect patterns of somatic mutation, including point mutations, structural variants and somatic retrotransposition. A collection of papers from the PCAWG Consortium describes non-coding mutations that drive cancer beyond those in the TERT promoter4; identifies new signatures of mutational processes that cause base substitutions, small insertions and deletions and structural variation5,6; analyses timings and patterns of tumour evolution7; describes the diverse transcriptional consequences of somatic mutation on splicing, expression levels, fusion genes and promoter activity8,9; and evaluates a range of more-specialized features of cancer genomes8,10-18.