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.

Cross-species oncogenomics offers insight into human muscle-invasive bladder cancer

BACKGROUND

In humans, muscle-invasive bladder cancer (MIBC) is highly aggressive and associated with a poor prognosis. With a high mutation load and large number of altered genes, strategies to delineate key driver events are necessary. Dogs and cats develop urothelial carcinoma (UC) with histological and clinical similarities to human MIBC. Cattle that graze on bracken fern also develop UC, associated with exposure to the carcinogen ptaquiloside. These species may represent relevant animal models of spontaneous and carcinogen-induced UC that can provide insight into human MIBC.

RESULTS

Whole-exome sequencing of domestic canine (n = 87) and feline (n = 23) UC, and comparative analysis with human MIBC reveals a lower mutation rate in animal cases and the absence of APOBEC mutational signatures. A convergence of driver genes (ARID1A, KDM6A, TP53, FAT1, and NRAS) is discovered, along with common focally amplified and deleted genes involved in regulation of the cell cycle and chromatin remodelling. We identify mismatch repair deficiency in a subset of canine and feline UCs with biallelic inactivation of MSH2. Bovine UC (n = 8) is distinctly different; we identify novel mutational signatures which are recapitulated in vitro in human urinary bladder UC cells treated with bracken fern extracts or purified ptaquiloside.

CONCLUSION

Canine and feline urinary bladder UC represent relevant models of MIBC in humans, and cross-species analysis can identify evolutionarily conserved driver genes. We characterize mutational signatures in bovine UC associated with bracken fern and ptaquiloside exposure, a human-linked cancer exposure. Our work demonstrates the relevance of cross-species comparative analysis in understanding both human and animal UC.

Next Generation Sequencing Workshop at the Royal Society of Medicine (London, May 2022): how genomics is on the path to modernizing genetic toxicology

The use of error-corrected Next Generation Sequencing (ecNG) to determine mutagenicity has been a subject of growing interest and potentially a disruptive technology that could supplement, and in time, replace current testing paradigms in preclinical safety assessment. Considering this, a Next Generation Sequencing Workshop was held at the Royal Society of Medicine in London in May 2022, supported by the United Kingdom Environmental Mutagen Society (UKEMS) and TwinStrand Biosciences (WA, USA), to discuss progress and future applications of this technology. In this meeting report, the invited speakers provide an overview of the Workshop topics covered and identify future directions for research. In the area of somatic mutagenesis, several speakers reviewed recent progress made with correlating ecNGS to classic in vivo transgenic rodent mutation assays as well as exploring the use of this technology directly in humans and animals, and in complex organoid models. Additionally, ecNGS has been used for detecting off-target effects of gene editing tools and emerging data suggest ecNGS potential to measure clonal expansion of cells carrying mutations in cancer driver genes as an early marker of carcinogenic potential and for direct human biomonitoring. As such, the workshop demonstrated the importance of raising awareness and support for advancing the science of ecNGS for mutagenesis, gene editing, and carcinogenesis research. Furthermore, the potential of this new technology to contribute to advances in drug and product development and improve safety assessment was extensively explored.

Abstract 218: Somatic mutations in single-cell derived alveolar organoids

The diversity in somatic variation across human tissues outshining human cancers has only recently been started to be appreciated. In particular, healthy individuals carry large numbers of cells and clonal expansions with mutations in driver genes commonly observed in cancer. We previously demonstrated that the majority of epithelial cells in the proximal airway accumulate DNA damage proportionally to an individual's smoking habits. However, the analysis of whole genomes indicated the abundance of a population of cells with a near-normal mutation burden, mirroring cellular populations found in never-smoking individuals. In contrast to the proximal airway, the genomics of cells in the distal airway of healthy individuals with diverse smoking histories, and thus the foundation of lung adenocarcinoma (LUAD), has not been investigated yet. To assess the influence of smoking on the distal airway and to compare the genomics of proximal and distal airway in the lung, we established an experimental procedure to create single-cell derived alveolar type II (AT2) organoids from primary human tissue of 9 patients. In particular, we focus on AT2 cells, given their hypothesised role as a potential cell of origin for LUAD.

Whole genome sequencing (WGS) of >400 organoids was utilised to assess the landscape of somatic variation, including single nucleotide variants (SNVs), insertions and deletions (indels), copy number and structural variation (CNVs, SVs). Focusing on SNVs, we found that the mutation burden of never-smokers is significantly lower compared to ex- and current-smokers. Emphasising the effect of tobacco smoking on the distal airway, this difference was driven by the presence of mutations related to SBS4 and SBS92, mutational signatures associated with tobacco smoking. Interestingly, and in contrast to the proximal airway, we did not find a difference in the mutation burden of ex- and current-smokers, suggesting that smoking damage to alveolar cells is more severe and long-term. Lastly, to characterise the driver mutation landscape in depth, we used targeted single-molecule sequencing of 250 known cancer genes to profile 15 AT2 single-cell suspensions as well as 22 frozen parenchyma biopsies. To date, our preliminary results indicate positive selection on lineage defining mutations including SFTPB and SFTPC, coding for surfactant as well as TP53. As we begin to explore these data further, we expect this work to provide unprecedented insights into the contribution of somatic and driver mutations to the earliest stages of lung cancer development.

Citation Format: Moritz Jakob Przybilla, Amany Ammar, Ryan Chuen Khaw, Andrew R. Lawson, Pantelis Nicola, Kate Davies, Zoe Frazer, Kate H. Gowers, Timothy M. Butler, Sarah E. Clarke, Inigo Martincorena, Sam M. Janes, Peter J. Campbell. Somatic mutations in single-cell derived alveolar organoids [abstract]. In: Proceedings of the American Association for Cancer Research Annual Meeting 2023; Part 1 (Regular and Invited Abstracts); 2023 Apr 14-19; Orlando, FL. Philadelphia (PA): AACR; Cancer Res 2023;83(7_Suppl):Abstract nr 218.

Substantial somatic genomic variation and selection for BCOR mutations in human induced pluripotent stem cells

We explored human induced pluripotent stem cells (hiPSCs) derived from different tissues to gain insights into genomic integrity at single-nucleotide resolution. We used genome sequencing data from two large hiPSC repositories involving 696 hiPSCs and daughter subclones. We find ultraviolet light (UV)-related damage in ~72% of skin fibroblast-derived hiPSCs (F-hiPSCs), occasionally resulting in substantial mutagenesis (up to 15 mutations per megabase). We demonstrate remarkable genomic heterogeneity between independent F-hiPSC clones derived during the same round of reprogramming due to oligoclonal fibroblast populations. In contrast, blood-derived hiPSCs (B-hiPSCs) had fewer mutations and no UV damage but a high prevalence of acquired BCOR mutations (26.9% of lines). We reveal strong selection pressure for BCOR mutations in F-hiPSCs and B-hiPSCs and provide evidence that they arise in vitro. Directed differentiation of hiPSCs and RNA sequencing showed that BCOR mutations have functional consequences. Our work strongly suggests that detailed nucleotide-resolution characterization is essential before using hiPSCs.

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.

Stage-stratified molecular profiling of non-muscle-invasive bladder cancer enhances biological, clinical, and therapeutic insight

Understanding the molecular determinants that underpin the clinical heterogeneity of non-muscle-invasive bladder cancer (NMIBC) is essential for prognostication and therapy development. Stage T1 disease in particular presents a high risk of progression and requires improved understanding. We present a detailed multi-omics study containing gene expression, copy number, and mutational profiles that show relationships to immune infiltration, disease recurrence, and progression to muscle invasion. We compare expression and genomic subtypes derived from all NMIBCs with those derived from the individual disease stages Ta and T1. We show that sufficient molecular heterogeneity exists within the separate stages to allow subclassification and that this is more clinically meaningful for stage T1 disease than that derived from all NMIBCs. This provides improved biological understanding and identifies subtypes of T1 tumors that may benefit from chemo- or immunotherapy.

Genomic reconstruction of the SARS-CoV-2 epidemic in England

The evolution of the severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) virus leads to new variants that warrant timely epidemiological characterization. Here we use the dense genomic surveillance data generated by the COVID-19 Genomics UK Consortium to reconstruct the dynamics of 71 different lineages in each of 315 English local authorities between September 2020 and June 2021. This analysis reveals a series of subepidemics that peaked in early autumn 2020, followed by a jump in transmissibility of the B.1.1.7/Alpha lineage. The Alpha variant grew when other lineages declined during the second national lockdown and regionally tiered restrictions between November and December 2020. A third more stringent national lockdown suppressed the Alpha variant and eliminated nearly all other lineages in early 2021. Yet a series of variants (most of which contained the spike E484K mutation) defied these trends and persisted at moderately increasing proportions. However, by accounting for sustained introductions, we found that the transmissibility of these variants is unlikely to have exceeded the transmissibility of the Alpha variant. Finally, B.1.617.2/Delta was repeatedly introduced in England and grew rapidly in early summer 2021, constituting approximately 98% of sampled SARS-CoV-2 genomes on 26 June 2021.

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 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.

Patterns of within-host genetic diversity in SARS-CoV-2

Monitoring the spread of SARS-CoV-2 and reconstructing transmission chains has become a major public health focus for many governments around the world. The modest mutation rate and rapid transmission of SARS-CoV-2 prevents the reconstruction of transmission chains from consensus genome sequences, but within-host genetic diversity could theoretically help identify close contacts. Here we describe the patterns of within-host diversity in 1181 SARS-CoV-2 samples sequenced to high depth in duplicate. 95.1% of samples show within-host mutations at detectable allele frequencies. Analyses of the mutational spectra revealed strong strand asymmetries suggestive of damage or RNA editing of the plus strand, rather than replication errors, dominating the accumulation of mutations during the SARS-CoV-2 pandemic. Within- and between-host diversity show strong purifying selection, particularly against nonsense mutations. Recurrent within-host mutations, many of which coincide with known phylogenetic homoplasies, display a spectrum and patterns of purifying selection more suggestive of mutational hotspots than recombination or convergent evolution. While allele frequencies suggest that most samples result from infection by a single lineage, we identify multiple putative examples of co-infection. Integrating these results into an epidemiological inference framework, we find that while sharing of within-host variants between samples could help the reconstruction of transmission chains, mutational hotspots and rare cases of superinfection can confound these analyses.

Characterizing genetic intra-tumor heterogeneity across 2,658 human cancer genomes

Intra-tumor heterogeneity (ITH) is a mechanism of therapeutic resistance and therefore an important clinical challenge. However, the extent, origin, and drivers of ITH across cancer types are poorly understood. To address this, we extensively characterize ITH across whole-genome sequences of 2,658 cancer samples spanning 38 cancer types. Nearly all informative samples (95.1%) contain evidence of distinct subclonal expansions with frequent branching relationships between subclones. We observe positive selection of subclonal driver mutations across most cancer types and identify cancer type-specific subclonal patterns of driver gene mutations, fusions, structural variants, and copy number alterations as well as dynamic changes in mutational processes between subclonal expansions. Our results underline the importance of ITH and its drivers in tumor evolution and provide a pan-cancer resource of comprehensively annotated subclonal events from whole-genome sequencing data.

Reliable detection of somatic mutations in solid tissues by laser-capture microdissection and low-input DNA sequencing

Somatic mutations accumulate in healthy tissues as we age, giving rise to cancer and potentially contributing to ageing. To study somatic mutations in non-neoplastic tissues, we developed a series of protocols to sequence the genomes of small populations of cells isolated from histological sections. Here, we describe a complete workflow that combines laser-capture microdissection (LCM) with low-input genome sequencing, while circumventing the use of whole-genome amplification (WGA). The protocol is subdivided broadly into four steps: tissue processing, LCM, low-input library generation and mutation calling and filtering. The tissue processing and LCM steps are provided as general guidelines that might require tailoring based on the specific requirements of the study at hand. Our protocol for low-input library generation uses enzymatic rather than acoustic fragmentation to generate WGA-free whole-genome libraries. Finally, the mutation calling and filtering strategy has been adapted from previously published protocols to account for artifacts introduced via library creation. To date, we have used this workflow to perform targeted and whole-genome sequencing of small populations of cells (typically 100-1,000 cells) in thousands of microbiopsies from a wide range of human tissues. The low-input DNA protocol is designed to be compatible with liquid handling platforms and make use of equipment and expertise standard to any core sequencing facility. However, obtaining low-input DNA material via LCM requires specialized equipment and expertise. The entire protocol from tissue reception through whole-genome library generation can be accomplished in as little as 1 week, although 2-3 weeks would be a more typical turnaround time.

Abstract 2342: Tobacco exposure and somatic mutations in normal bronchial epithelia

Background

Lung cancer is the leading cause of cancer death, of which 80-90% are attributed to tobacco smoking. Our understanding of how tobacco exposure affects mutational burden, mutational signatures, driver mutations and clonal dynamics in normal lung tissue is limited. Also, genetic differences between normal cells and cancer cells in lung has not been fully elucidated.

Methods

To access the landscape of somatic mutations in normal bronchial epithelia, we sequenced whole genomes of 632 single-cell-derived colonies of lung basal cells from 16 patients, including 3 children, 4 never-smokers, 6 ex-smokers and 3 current smokers. Five patients had squamous cell carcinomas or carcinoma in situ, which we also sequenced to compare the genetic alteration between normal and cancer cells.

Results

A positive correlation between the number of base substitutions and age was observed, and 22 somatic mutations accumulated per year (95% confidence interval:20-25; P=10−8). Previous or current smoking significantly increased mutational burden: 2330 substitutions in ex-smokers and 5300 in current smokers. In addition, tobacco smoking is massively increasing both between-subject and within-subject variance of mutation burden. A population of cells in subjects with smoking history had mutation burdens equivalent to that expected for never-smokers: these cells had less damage from tobacco-specific mutational processes, and were four-fold more frequent in ex-smokers than current smokers. Signature analysis revealed that the same mutational signatures seen in lung cancers operate both in patients with and without smoking history: endogeneous signatures, including COSMIC signatures 1 and 5, and APOBEC-related signatures. Three mutational signatures were largely restricted to current or ex-smokers, including known smoking-related COSMIC signature 4 characterized by C>A transversions, signature 16 characterized by T>C mutations with extremely strong transcription strand bias and a new signature characterized by T>A and T>C mutations. dN/dS method identified driver genes in normal bronchial epithelium, including NOTCH1, TP53 and FAT1, which were overlapped with those seen in squamous cell lung cancers and other normal squamous tissues such as esophagus and skin. Driver mutations increased in frequency with age, affecting 4-14% of cells in middle-aged never-smokers. In current smokers, ≥25% of cells carried driver mutations and 0-6% cells had 2 or even 3 drivers. Layering driver mutations onto phylogenetic trees revealed that driver mutations occurred in early life. Compared to the normal bronchial epithelial cells, lung cancers and precancerous lesions were characterized by extensive copy number changes or structural variants and distinct set of driver mutations.

Conclusions

Tobacco smoking increases mutation burden, cell-to-cell heterogeneity and driver mutations. Our data of genetic lesions in normal bronchial cells provides insights into genetic alterations that drive carcinogenesis in lung.

Citation Format: Kenichi Yoshida, Kate HC Gowers, Henry Lee-Six, Deepak P. Chandrasekharan, Tim Coorens, Elizabeth F. Maughan, Kathryn Beal, Andrew Menzies, Fraser R. Millar, Elizabeth Anderson, Sarah E. Clarke, Adam Pennycuick, Ricky M. Thackeray, Colin R. Butler, Nobuyuki Kakiuchi, Tomonori Hirano, Robert E. Hynds, Michael R. Stratton, Inigo Martincorena, Sam M. Janes, Peter J. Campbell. Tobacco exposure and somatic mutations in normal bronchial epithelia [abstract]. In: Proceedings of the Annual Meeting of the American Association for Cancer Research 2020; 2020 Apr 27-28 and Jun 22-24. Philadelphia (PA): AACR; Cancer Res 2020;80(16 Suppl):Abstract nr 2342.

Somatic Evolution in Non-neoplastic IBD-Affected Colon

Inflammatory bowel disease (IBD) is a chronic inflammatory disease associated with increased risk of gastrointestinal cancers. We whole-genome sequenced 446 colonic crypts from 46 IBD patients and compared these to 412 crypts from 41 non-IBD controls from our previous publication on the mutation landscape of the normal colon. The average mutation rate of affected colonic epithelial cells is 2.4-fold that of healthy colon, and this increase is mostly driven by acceleration of mutational processes ubiquitously observed in normal colon. In contrast to the normal colon, where clonal expansions outside the confines of the crypt are rare, we observed widespread millimeter-scale clonal expansions. We discovered non-synonymous mutations in ARID1A, FBXW7, PIGR, ZC3H12A, and genes in the interleukin 17 and Toll-like receptor pathways, under positive selection in IBD. These results suggest distinct selection mechanisms in the colitis-affected colon and that somatic mutations potentially play a causal role in IBD pathogenesis.

Spatial competition shapes the dynamic mutational landscape of normal esophageal epithelium

During aging, progenitor cells acquire mutations, which may generate clones that colonize the surrounding tissue. By middle age, normal human tissues, including the esophageal epithelium (EE), become a patchwork of mutant clones. Despite their relevance for understanding aging and cancer, the processes that underpin mutational selection in normal tissues remain poorly understood. Here, we investigated this issue in the esophageal epithelium of mutagen-treated mice. Deep sequencing identified numerous mutant clones with multiple genes under positive selection, including Notch1, Notch2 and Trp53, which are also selected in human esophageal epithelium. Transgenic lineage tracing revealed strong clonal competition that evolved over time. Clone dynamics were consistent with a simple model in which the proliferative advantage conferred by positively selected mutations depends on the nature of the neighboring cells. When clones with similar competitive fitness collide, mutant cell fate reverts towards homeostasis, a constraint that explains how selection operates in normal-appearing epithelium.

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.

Comprehensive molecular characterization of mitochondrial genomes in human cancers

Mitochondria are essential cellular organelles that play critical roles in cancer. Here, as part of the International Cancer Genome Consortium/The Cancer Genome Atlas Pan-Cancer Analysis of Whole Genomes Consortium, which aggregated whole-genome sequencing data from 2,658 cancers across 38 tumor types, we performed a multidimensional, integrated characterization of mitochondrial genomes and related RNA sequencing data. Our analysis presents the most definitive mutational landscape of mitochondrial genomes and identifies several hypermutated cases. Truncating mutations are markedly enriched in kidney, colorectal and thyroid cancers, suggesting oncogenic effects with the activation of signaling pathways. We find frequent somatic nuclear transfers of mitochondrial DNA, some of which disrupt therapeutic target genes. Mitochondrial copy number varies greatly within and across cancers and correlates with clinical variables. Co-expression analysis highlights the function of mitochondrial genes in oxidative phosphorylation, DNA repair and the cell cycle, and shows their connections with clinically actionable genes. Our study lays a foundation for translating mitochondrial biology into clinical applications.

Long-term expansion, genomic stability and in vivo safety of adult human pancreas organoids

BACKGROUND

Pancreatic organoid systems have recently been described for the in vitro culture of pancreatic ductal cells from mouse and human. Mouse pancreatic organoids exhibit unlimited expansion potential, while previously reported human pancreas organoid (hPO) cultures do not expand efficiently long-term in a chemically defined, serum-free medium. We sought to generate a 3D culture system for long-term expansion of human pancreas ductal cells as hPOs to serve as the basis for studies of human pancreas ductal epithelium, exocrine pancreatic diseases and the development of a genomically stable replacement cell therapy for diabetes mellitus.

RESULTS

Our chemically defined, serum-free, human pancreas organoid culture medium supports the generation and expansion of hPOs with high efficiency from both fresh and cryopreserved primary tissue. hPOs can be expanded from a single cell, enabling their genetic manipulation and generation of clonal cultures. hPOs expanded for months in vitro maintain their ductal morphology, biomarker expression and chromosomal integrity. Xenografts of hPOs survive long-term in vivo when transplanted into the pancreas of immunodeficient mice. Notably, mouse orthotopic transplants show no signs of tumorigenicity. Crucially, our medium also supports the establishment and expansion of hPOs in a chemically defined, modifiable and scalable, biomimetic hydrogel.

CONCLUSIONS

hPOs can be expanded long-term, from both fresh and cryopreserved human pancreas tissue in a chemically defined, serum-free medium with no detectable tumorigenicity. hPOs can be clonally expanded, genetically manipulated and are amenable to culture in a chemically defined hydrogel. hPOs therefore represent an abundant source of pancreas ductal cells that retain the characteristics of the tissue-of-origin, which opens up avenues for modelling diseases of the ductal epithelium and increasing understanding of human pancreas exocrine biology as well as for potentially producing insulin-secreting cells for the treatment of diabetes.

Genomic evidence supports a clonal diaspora model for metastases of esophageal adenocarcinoma

The poor outcomes in esophageal adenocarcinoma (EAC) prompted us to interrogate the pattern and timing of metastatic spread. Whole-genome sequencing and phylogenetic analysis of 388 samples across 18 individuals with EAC showed, in 90% of patients, that multiple subclones from the primary tumor spread very rapidly from the primary site to form multiple metastases, including lymph nodes and distant tissues-a mode of dissemination that we term 'clonal diaspora'. Metastatic subclones at autopsy were present in tissue and blood samples from earlier time points. These findings have implications for our understanding and clinical evaluation of EAC.

The landscape of somatic mutation in normal colorectal epithelial cells

The colorectal adenoma-carcinoma sequence has provided a paradigmatic framework for understanding the successive somatic genetic changes and consequent clonal expansions that lead to cancer1. However, our understanding of the earliest phases of colorectal neoplastic changes-which may occur in morphologically normal tissue-is comparatively limited, as for most cancer types. Here we use whole-genome sequencing to analyse hundreds of normal crypts from 42 individuals. Signatures of multiple mutational processes were revealed; some of these were ubiquitous and continuous, whereas others were only found in some individuals, in some crypts or during certain periods of life. Probable driver mutations were present in around 1% of normal colorectal crypts in middle-aged individuals, indicating that adenomas and carcinomas are rare outcomes of a pervasive process of neoplastic change across morphologically normal colorectal epithelium. Colorectal cancers exhibit substantially increased mutational burdens relative to normal cells. Sequencing normal colorectal cells provides quantitative insights into the genomic and clonal evolution of cancer.

Somatic mutations and clonal dynamics in healthy and cirrhotic human liver

The most common causes of chronic liver disease are excess alcohol intake, viral hepatitis and non-alcoholic fatty liver disease, with the clinical spectrum ranging in severity from hepatic inflammation to cirrhosis, liver failure or hepatocellular carcinoma (HCC). The genome of HCC exhibits diverse mutational signatures, resulting in recurrent mutations across more than 30 cancer genes1-7. Stem cells from normal livers have a low mutational burden and limited diversity of signatures8, which suggests that the complexity of HCC arises during the progression to chronic liver disease and subsequent malignant transformation. Here, by sequencing whole genomes of 482 microdissections of 100-500 hepatocytes from 5 normal and 9 cirrhotic livers, we show that cirrhotic liver has a higher mutational burden than normal liver. Although rare in normal hepatocytes, structural variants, including chromothripsis, were prominent in cirrhosis. Driver mutations, such as point mutations and structural variants, affected 1-5% of clones. Clonal expansions of millimetres in diameter occurred in cirrhosis, with clones sequestered by the bands of fibrosis that surround regenerative nodules. Some mutational signatures were universal and equally active in both non-malignant hepatocytes and HCCs; some were substantially more active in HCCs than chronic liver disease; and others-arising from exogenous exposures-were present in a subset of patients. The activity of exogenous signatures between adjacent cirrhotic nodules varied by up to tenfold within each patient, as a result of clone-specific and microenvironmental forces. Synchronous HCCs exhibited the same mutational signatures as background cirrhotic liver, but with higher burden. Somatic mutations chronicle the exposures, toxicity, regeneration and clonal structure of liver tissue as it progresses from health to disease.

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.

Abstract 970: The mutational landscape of normal human endometrial epithelium

Human endometrium is a highly dynamic tissue that undergoes hundreds of cycles of breakdown, rapid repair and remodelling in response to the oscillating levels of oestrogen and progesterone during female reproductive years. The marked regenerative capacity of this tissue’s epithelial compartment is maintained by intra-glandular adult stem cells (ASCs) that reside within the stratum basalis which is retained during menstruation. Although the endometrial ASCs were first described over a decade ago, they remain poorly characterised in comparison to their counterparts in other tissues, such as the small and the large intestines. In particular, the size of the stem cell pool within individual glands, the rates of their division, and mutational landscape are largely unknown. In this study, we isolated 215 morphologically normal endometrial glands from women aged 19 to 81 using laser capture microscopy. Analysis of whole-genome sequencing data identified that the overwhelming majority of the glands were clonal cell populations, and thus originating from a single ASC. Somatic mutations were found to accumulate at a linear rate during adult life. Elevated body mass index (BMI), a well-recognised risk factor for endometrial cancer, accelerated the rate of mutation acquisition. Surprisingly, despite the heterogeneity in age, reproductive history and BMI in our cohort, we find relatively homogenous mutational processes within normal endometrium. Comparison with cancer, shows lower somatic mutation burden and fewer operative signatures. Remarkably, we not only identify recurrent acquisitions of certain cancer-associated mutations, particularly those that are advantageous to cell growth, proliferation and migration, but also show that such events occur early in life, potentially even before adolescence. Over time, these mutant ASCs serve as a reservoir for the acquisition of further driver mutations to the extent that in some cases, the entire sampled endometrium becomes ‘neoplastic’ on the genomic level while still retaining the apparently normal phenotype. In older individuals, we observe a shift in the spectrum of acquired cancer-associated mutations, possibly reflecting post-menopausal changes in the levels of sex-steroid hormones and the resultant tissue microenvironment.

Citation Format: Luiza Moore, Daniel Leongamornlert, Tim Coorens, Mathijs Sanders, Peter Ellis, Francesco Maura, Kevin Dawson, Simon F. Brunner, Jyoti Nangalia, Henry Lee-Six, Raheleh Rahbari, Patrick Tarpey, Yvette Hooks, Krishnaa Mahbubani, Christine A. Iacobuzio-Donahue, Jan J. Brosens, Inigo Martincorena, Kourosh Saeb-Parsy, Peter J. Campbell, Michael R. Stratton. The mutational landscape of normal human endometrial epithelium [abstract]. In: Proceedings of the American Association for Cancer Research Annual Meeting 2019; 2019 Mar 29-Apr 3; Atlanta, GA. Philadelphia (PA): AACR; Cancer Res 2019;79(13 Suppl):Abstract nr 970.

Somatic mutation and clonal expansions in human tissues

Recent sequencing studies on healthy skin and esophagus have found that, as we age, these tissues become colonized by mutant clones of cells carrying driver mutations in traditional cancer genes. This comment summarizes these findings and discusses their possible implications for our understanding of cancer, ageing, and other diseases.

Cross-species genomic landscape comparison of human mucosal melanoma with canine oral and equine melanoma

Mucosal melanoma is a rare and poorly characterized subtype of human melanoma. Here we perform a cross-species analysis by sequencing tumor-germline pairs from 46 primary human muscosal, 65 primary canine oral and 28 primary equine melanoma cases from mucosal sites. Analysis of these data reveals recurrently mutated driver genes shared between species such as NRAS, FAT4, PTPRJ, TP53 and PTEN, and pathogenic germline alleles of BRCA1, BRCA2 and TP53. We identify a UV mutation signature in a small number of samples, including human cases from the lip and nasal mucosa. A cross-species comparative analysis of recurrent copy number alterations identifies several candidate drivers including MDM2, B2M, KNSTRN and BUB1B. Comparison of somatic mutations in recurrences and metastases to those in the primary tumor suggests pervasive intra-tumor heterogeneity. Collectively, these studies suggest a convergence of some genetic changes in mucosal melanomas between species but also distinctly different paths to tumorigenesis.

An integrated genomic analysis of anaplastic meningioma identifies prognostic molecular signatures

Anaplastic meningioma is a rare and aggressive brain tumor characterised by intractable recurrences and dismal outcomes. Here, we present an integrated analysis of the whole genome, transcriptome and methylation profiles of primary and recurrent anaplastic meningioma. A key finding was the delineation of distinct molecular subgroups that were associated with diametrically opposed survival outcomes. Relative to lower grade meningiomas, anaplastic tumors harbored frequent driver mutations in SWI/SNF complex genes, which were confined to the poor prognosis subgroup. Aggressive disease was further characterised by transcriptional evidence of increased PRC2 activity, stemness and epithelial-to-mesenchymal transition. Our analyses discern biologically distinct variants of anaplastic meningioma with prognostic and therapeutic significance.

Population dynamics of normal human blood inferred from somatic mutations

Haematopoietic stem cells drive blood production, but their population size and lifetime dynamics have not been quantified directly in humans. Here we identified 129,582 spontaneous, genome-wide somatic mutations in 140 single-cell-derived haematopoietic stem and progenitor colonies from a healthy 59-year-old man and applied population-genetics approaches to reconstruct clonal dynamics. Cell divisions from early embryogenesis were evident in the phylogenetic tree; all blood cells were derived from a common ancestor that preceded gastrulation. The size of the stem cell population grew steadily in early life, reaching a stable plateau by adolescence. We estimate the numbers of haematopoietic stem cells that are actively making white blood cells at any one time to be in the range of 50,000-200,000. We observed adult haematopoietic stem cell clones that generate multilineage outputs, including granulocytes and B lymphocytes. Harnessing naturally occurring mutations to report the clonal architecture of an organ enables the high-resolution reconstruction of somatic cell dynamics in humans.

Genomic patterns of progression in smoldering multiple myeloma

We analyzed whole genomes of unique paired samples from smoldering multiple myeloma (SMM) patients progressing to multiple myeloma (MM). We report that the genomic landscape, including mutational profile and structural rearrangements at the smoldering stage is very similar to MM. Paired sample analysis shows two different patterns of progression: a "static progression model", where the subclonal architecture is retained as the disease progressed to MM suggesting that progression solely reflects the time needed to accumulate a sufficient disease burden; and a "spontaneous evolution model", where a change in the subclonal composition is observed. We also observe that activation-induced cytidine deaminase plays a major role in shaping the mutational landscape of early subclinical phases, while progression is driven by APOBEC cytidine deaminases. These results provide a unique insight into myelomagenesis with potential implications for the definition of smoldering disease and timing of treatment initiation.

Prediction of acute myeloid leukaemia risk in healthy individuals

The incidence of acute myeloid leukaemia (AML) increases with age and mortality exceeds 90% when diagnosed after age 65. Most cases arise without any detectable early symptoms and patients usually present with the acute complications of bone marrow failure1. The onset of such de novo AML cases is typically preceded by the accumulation of somatic mutations in preleukaemic haematopoietic stem and progenitor cells (HSPCs) that undergo clonal expansion2,3. However, recurrent AML mutations also accumulate in HSPCs during ageing of healthy individuals who do not develop AML, a phenomenon referred to as age-related clonal haematopoiesis (ARCH)4-8. Here we use deep sequencing to analyse genes that are recurrently mutated in AML to distinguish between individuals who have a high risk of developing AML and those with benign ARCH. We analysed peripheral blood cells from 95 individuals that were obtained on average 6.3 years before AML diagnosis (pre-AML group), together with 414 unselected age- and gender-matched individuals (control group). Pre-AML cases were distinct from controls and had more mutations per sample, higher variant allele frequencies, indicating greater clonal expansion, and showed enrichment of mutations in specific genes. Genetic parameters were used to derive a model that accurately predicted AML-free survival; this model was validated in an independent cohort of 29 pre-AML cases and 262 controls. Because AML is rare, we also developed an AML predictive model using a large electronic health record database that identified individuals at greater risk. Collectively our findings provide proof-of-concept that it is possible to discriminate ARCH from pre-AML many years before malignant transformation. This could in future enable earlier detection and monitoring, and may help to inform intervention.

Sequencing of prostate cancers identifies new cancer genes, routes of progression and drug targets

Prostate cancer represents a substantial clinical challenge because it is difficult to predict outcome and advanced disease is often fatal. We sequenced the whole genomes of 112 primary and metastatic prostate cancer samples. From joint analysis of these cancers with those from previous studies (930 cancers in total), we found evidence for 22 previously unidentified putative driver genes harboring coding mutations, as well as evidence for NEAT1 and FOXA1 acting as drivers through noncoding mutations. Through the temporal dissection of aberrations, we identified driver mutations specifically associated with steps in the progression of prostate cancer, establishing, for example, loss of CHD1 and BRCA2 as early events in cancer development of ETS fusion-negative cancers. Computational chemogenomic (canSAR) analysis of prostate cancer mutations identified 11 targets of approved drugs, 7 targets of investigational drugs, and 62 targets of compounds that may be active and should be considered candidates for future clinical trials.

Timing the Landmark Events in the Evolution of Clear Cell Renal Cell Cancer: TRACERx Renal

Clear cell renal cell carcinoma (ccRCC) is characterized by near-universal loss of the short arm of chromosome 3, deleting several tumor suppressor genes. We analyzed whole genomes from 95 biopsies across 33 patients with clear cell renal cell carcinoma. We find hotspots of point mutations in the 5' UTR of TERT, targeting a MYC-MAX-MAD1 repressor associated with telomere lengthening. The most common structural abnormality generates simultaneous 3p loss and 5q gain (36% patients), typically through chromothripsis. This event occurs in childhood or adolescence, generally as the initiating event that precedes emergence of the tumor's most recent common ancestor by years to decades. Similar genomic changes drive inherited ccRCC. Modeling differences in age incidence between inherited and sporadic cancers suggests that the number of cells with 3p loss capable of initiating sporadic tumors is no more than a few hundred. Early development of ccRCC follows well-defined evolutionary trajectories, offering opportunity for early intervention.

The driver landscape of sporadic chordoma

Chordoma is a malignant, often incurable bone tumour showing notochordal differentiation. Here, we defined the somatic driver landscape of 104 cases of sporadic chordoma. We reveal somatic duplications of the notochordal transcription factor brachyury (T) in up to 27% of cases. These variants recapitulate the rearrangement architecture of the pathogenic germline duplications of T that underlie familial chordoma. In addition, we find potentially clinically actionable PI3K signalling mutations in 16% of cases. Intriguingly, one of the most frequently altered genes, mutated exclusively by inactivating mutation, was LYST (10%), which may represent a novel cancer gene in chordoma.Chordoma is a rare often incurable malignant bone tumour. Here, the authors investigate driver mutations of sporadic chordoma in 104 cases, revealing duplications in notochordal transcription factor brachyury (T), PI3K signalling mutations, and mutations in LYST, a potential novel cancer gene in chordoma.

Genomic Evolution of Breast Cancer Metastasis and Relapse

Patterns of genomic evolution between primary and metastatic breast cancer have not been studied in large numbers, despite patients with metastatic breast cancer having dismal survival. We sequenced whole genomes or a panel of 365 genes on 299 samples from 170 patients with locally relapsed or metastatic breast cancer. Several lines of analysis indicate that clones seeding metastasis or relapse disseminate late from primary tumors, but continue to acquire mutations, mostly accessing the same mutational processes active in the primary tumor. Most distant metastases acquired driver mutations not seen in the primary tumor, drawing from a wider repertoire of cancer genes than early drivers. These include a number of clinically actionable alterations and mutations inactivating SWI-SNF and JAK2-STAT3 pathways.

Recurrent mutation of IGF signalling genes and distinct patterns of genomic rearrangement in osteosarcoma

Osteosarcoma is a primary malignancy of bone that affects children and adults. Here, we present the largest sequencing study of osteosarcoma to date, comprising 112 childhood and adult tumours encompassing all major histological subtypes. A key finding of our study is the identification of mutations in insulin-like growth factor (IGF) signalling genes in 8/112 (7%) of cases. We validate this observation using fluorescence in situ hybridization (FISH) in an additional 87 osteosarcomas, with IGF1 receptor (IGF1R) amplification observed in 14% of tumours. These findings may inform patient selection in future trials of IGF1R inhibitors in osteosarcoma. Analysing patterns of mutation, we identify distinct rearrangement profiles including a process characterized by chromothripsis and amplification. This process operates recurrently at discrete genomic regions and generates driver mutations. It may represent an age-independent mutational mechanism that contributes to the development of osteosarcoma in children and adults alike.

GOTHiC, a probabilistic model to resolve complex biases and to identify real interactions in Hi-C data

Hi-C is one of the main methods for investigating spatial co-localisation of DNA in the nucleus. However, the raw sequencing data obtained from Hi-C experiments suffer from large biases and spurious contacts, making it difficult to identify true interactions. Existing methods use complex models to account for biases and do not provide a significance threshold for detecting interactions. Here we introduce a simple binomial probabilistic model that resolves complex biases and distinguishes between true and false interactions. The model corrects biases of known and unknown origin and yields a p-value for each interaction, providing a reliable threshold based on significance. We demonstrate this experimentally by testing the method against a random ligation dataset. Our method outperforms previous methods and provides a statistical framework for further data analysis, such as comparisons of Hi-C interactions between different conditions. GOTHiC is available as a BioConductor package (http://www.bioconductor.org/packages/release/bioc/html/GOTHiC.html).

Genome-wide chemical mutagenesis screens allow unbiased saturation of the cancer genome and identification of drug resistance mutations

Drug resistance is an almost inevitable consequence of cancer therapy and ultimately proves fatal for the majority of patients. In many cases, this is the consequence of specific gene mutations that have the potential to be targeted to resensitize the tumor. The ability to uniformly saturate the genome with point mutations without chromosome or nucleotide sequence context bias would open the door to identify all putative drug resistance mutations in cancer models. Here, we describe such a method for elucidating drug resistance mechanisms using genome-wide chemical mutagenesis allied to next-generation sequencing. We show that chemically mutagenizing the genome of cancer cells dramatically increases the number of drug-resistant clones and allows the detection of both known and novel drug resistance mutations. We used an efficient computational process that allows for the rapid identification of involved pathways and druggable targets. Such a priori knowledge would greatly empower serial monitoring strategies for drug resistance in the clinic as well as the development of trials for drug-resistant patients.

Somatic mutations reveal asymmetric cellular dynamics in the early human embryo

Somatic cells acquire mutations throughout the course of an individual's life. Mutations occurring early in embryogenesis are often present in a substantial proportion of, but not all, cells in postnatal humans and thus have particular characteristics and effects. Depending on their location in the genome and the proportion of cells they are present in, these mosaic mutations can cause a wide range of genetic disease syndromes and predispose carriers to cancer. They have a high chance of being transmitted to offspring as de novo germline mutations and, in principle, can provide insights into early human embryonic cell lineages and their contributions to adult tissues. Although it is known that gross chromosomal abnormalities are remarkably common in early human embryos, our understanding of early embryonic somatic mutations is very limited. Here we use whole-genome sequences of normal blood from 241 adults to identify 163 early embryonic mutations. We estimate that approximately three base substitution mutations occur per cell per cell-doubling event in early human embryogenesis and these are mainly attributable to two known mutational signatures. We used the mutations to reconstruct developmental lineages of adult cells and demonstrate that the two daughter cells of many early embryonic cell-doubling events contribute asymmetrically to adult blood at an approximately 2:1 ratio. This study therefore provides insights into the mutation rates, mutational processes and developmental outcomes of cell dynamics that operate during early human embryogenesis.