Search

[Introduction] Proliferative lesions in the breast have been implicated in the development of breast cancer. Previous studies showed that some proliferative lesions and adjacent breast cancers shared common genetic alterations, suggesting that these originated from the same ancestral cell. However, the clonal structure of normal epithelia and their clonal history during evolution to cancer are poorly understood. In this study, we analyzed genetic profiles of normal epithelia and proliferative lesions in the cancer-borne breast to illustrate the clonal evolution of cancer from a normal epithelial cell. [Methods] Single cell-derived organoids (n=47) were established from breast milk of 4 healthy women aged 22–36 and normal breast tissue of 15 breast cancer patients aged 29–83 to evaluate somatic mutation rate in normal epithelial cells. Multiple normal lobules and proliferative lesions together with cancer lesions were collected using laser-capture micro-dissection (LCM) from fresh frozen (n=3) or formalin-fixed paraffin-embedded (n=5) surgical specimens in 9 premenopausal breast cancer patients. Somatic mutations and copy number alterations were evaluated using whole-genome sequencing. [Results] The mutation profile of single cell-derived organoids suggests that somatic mutations accumulate in normal mammary epithelial cells at a constant rate of 19.4/genome/year before menopause, and the mutation rate decreases to 6.9/genome/year after menopause. Parity was negatively associated with mutation number (-49.3 per life birth). In total, we analyzed 143 LCM samples, including those from 72 normal lobules, 43 proliferative lesions, and 19 non-invasive and 9 invasive cancer samples. Five cases showed a large expansion of proliferative lesions sharing a substantial number of somatic mutations with cancer. These lesions expanded over a distance of 35-90 mm, sharing tens to hundreds of mutations including those in breast cancer-related driver genes, such as PIK3CA, AKT1, GATA3, CBFB and PTEN, while harboring private mutations or copy number alterations of their own. Of interest, the cancers in 4 out of these 5 cases was luminal-A type invasive ducal carcinoma or ER-positive HER2-negative ductal carcinoma in situ, and characterized in common by the presence of der(1;16), concurrent whole-arm 1q gain and 16q loss, in both cancer and proliferative lesions. Phylogenetic analysis adapted with the mutation rate in normal cells predicted that der(1;16) had been acquired between puberty and early 20’s, and the common ancestors of non-cancerous and cancerous lesions emerged by early 30’s, >10 years earlier than at the time of cancer diagnosis. By contrast, analysis of non-cancerous lobules unrelated to cancer showed that der(1;16)-negative non-cancer clones that had emerged after puberty stayed within a single lobule or spatially confined to adjacent lobules and rarely expanded to a large area as observed for those carrying der(1;16), even if the clones had acquired mutations in driver genes such as PIK3CA and PIK3R1, which highlighted the role of der(1;16) in wide clonal expansion. [Conclusions] Our results suggest that in some breast cancer cases, particularly in those with der(1;16), a highly recurrent translocation accounting for the major subset of Luminal A breast cancer, the clones with the funder driver alterations expanded macroscopically long before the onset of cancer, in which further clonal evolutions recursively occur multi-focally, giving rise to multiple proliferative lesions and ultimately, invasive cancers. Our findings provide new insight into the early development of breast cancer. Citation Format: Tomomi Nishimura, Nobuyuki Kakiuchi, Kenichi Yoshida, Takaki Sakurai, Tatsuki R. Kataoka, Eiji Kondoh, Yoshitsugu Chigusa, Masahiko Kawai, Morio Sawada, Takuya Inoue, Yasuhide Takeuchi, Hirona Maeda, Satoko Baba, Yusuke Shiozawa, Ryunosuke Saiki, Masahiro M. Nakagawa, Yasuhito Nannya, Yotaro Ochi, Tomonori Hirano, Yukiko Inagaki-Kawata, Kosuke Aoki, Masahiro Hirata, Eiji Suzuki, Masahiro Takada, Masahiro Kawashima, Kosuke Kawaguchi, Kenichi Chiba, Yuichi Shiraishi, Junko Takita, Satoru Miyano, Masaki Mandai, Kengo Takeuchi, Hironori Haga, Masakazu Toi, Seishi Ogawa. Clonal evolution of mammary epithelial cells into breast cancers [abstract]. In: Proceedings of the 2022 San Antonio Breast Cancer Symposium; 2022 Dec 6-10; San Antonio, TX. Philadelphia (PA): AACR; Cancer Res 2023;83(5 Suppl):Abstract nr P5-13-04.

Esophageal squamous cell carcinoma (ESCC) often recurs after chemoradiotherapy, and the prognosis of ESCC after chemoradiotherapy has not improved over the past few decades. The mutation process in chemoradiotherapy-resistant clones and the functional relevance of genetic alterations remain unclear. To address these problems, we performed whole-exome sequencing of 52 tumor samples from 33 patients with ESCC who received radiotherapy combined with 5-fluorouracil/platinum. In multiregion analyses of pretreatment and locally recurrent lesions from five cases, most driver gene-altered clones remained under chemoradiotherapy selection pressure, while few driver gene alterations were acquired at recurrence. The mutation signatures of recurrent ESCC, including increased deletion frequency and platinum dose-dependent base substitution signatures, were substantially different from those of primary ESCC and reflected the iatrogenic impacts of chemoradiotherapy. Single-region analysis of 28 pretreatment tumors indicated that focal copy-number gain at the MYC locus was significantly associated with poor progression-free survival and overall survival after chemoradiotherapy. MYC gain remained throughout the chemoradiotherapy course and potentially contributes to intrinsic resistance to chemoradiotherapy. Consistent with these findings, MYC copy number and mRNA and protein levels in ESCC cell lines correlated positively with resistance to radiotherapy, and MYC knockdown improved sensitivity to radiotherapy. Overall, these data characterize the clonal evolution process induced by chemoradiotherapy and clinically relevant associations for genetic alterations in ESCC. These findings increase our understanding of therapeutic resistance and support the rationale for precision chemoradiotherapy.Significance:Whole-exome sequencing reveals the genetic evolution of ESCC during chemoradiotherapy, highlighting MYC gain in pretreatment tumors as a potential marker of therapy resistance.

table S1 : recombination detection primers and qPCR primers. tableS2 : antibodies for immunohistochemistory and two chamber assay. tableS3 : list of PDAC patients

Supplementary table 1-9

Supplemental Table 1. PanNETs incidence per genotype. Figure S1: Representative images of LacZ staining in the pancreas in Pdx1-Cre;Rosa26R mice; S2: Experimental strategy for generation of Pdx1-Cre;Trp53R172H mice. S3: H&E, immunohistochemical and immunofluorescence staining of islet cells in Pdx1-Cre and Pdx1-Cre;Trp53R172H;Rbf/f mice at 2 months of age.

Histological analysis of ICC and the incidence of HCC developed in AlbcreKP; RosaNotchOE and AlbcreKP; RosaNotchOE; Hes1flox/flox mice at 8 weeks of age.

Lollipop plot of NUP93 comparing EMPD and PanCancer Atlas dataset

Myofibroblasts secrete Cola 1a1 and Cxcl12 in tumor environment.

Leading edge genes in gene set enrichment analysys using HALLMARK gene set of KRAS_SIGNALING_UP.

Intestinal metaplasia (IM) is a risk factor for gastric cancer following infection with Helicobacter pylori. To explore the susceptibility of pure gastric IM to cancer development, we investigated genetic alterations in single IM gastric glands. We isolated 50 single IM or non-IM glands from the inflamed gastric mucosa of 11 patients with intramucosal gastric carcinoma (IGC) and 4 patients without IGC; 19 single glands in the noninflamed gastric mucosa of 11 individuals from our cohort and previous dataset were also included as controls. Whole-exome sequencing of single glands revealed significantly higher accumulation of somatic mutations in various genes within IM glands compared with non-IM glands. Clonal ordering analysis showed that IM glands expanded to form clusters with shared mutations. In addition, targeted-capture deep sequencing and copy number (CN) analyses were performed in 96 clustered IM or non-IM gastric glands from 26 patients with IGC. CN analyses were also performed on 41 IGC samples and The Cancer Genome Atlas-Stomach Adenocarcinoma datasets. These analyses revealed that polyclonally expanded IM commonly acquired CN aberrations (CNA), including amplification of chromosomes 8, 20, and 2. A large portion of clustered IM glands typically consisted of common CNAs rather than other cancer-related mutations. Moreover, the CNA patterns of clustered IM glands were similar to those of IGC, indicative of precancerous conditions. Taken together, these findings suggest that, in the gastric mucosa inflamed with H. pylori infection, IM glands expand via acquisition of CNAs comparable with those of IGC, contributing to field cancerization.Significance:This study contributes to our understanding of gastric intestinal metaplasia as a risk factor for gastric adenocarcinoma via their multifocal expansion and acquisition of CNAs and somatic mutations.

Intrahepatic cholangiocarcinoma (ICC) is frequently driven by aberrant KRAS activation and develops in the liver with chronic inflammation. Although the Notch signaling pathway is critically involved in ICC development, detailed mechanisms of Notch-driven ICC development are still unknown. Here, we use mice whose Notch signaling is genetically engineered to show that the Notch signaling pathway, specifically the Notch/Hes1 axis, plays an essential role in expanding ductular cells in the liver with chronic inflammation or oncogenic Kras activation. Activation of Notch1 enhanced the development of proliferating ductal cells (PDC) in injured livers, while depletion of Hes1 led to suppression. In correlation with PDC expansion, ICC development was also regulated by the Notch/Hes1 axis and suppressed by Hes1 depletion. Lineage-tracing experiments using EpcamcreERT2 mice further confirmed that Hes1 plays a critical role in the induction of PDC and that ICC could originate from PDC. Analysis of human ICC specimens showed PDC in nonneoplastic background tissues, confirming HES1 expression in both PDC and ICC tumor cells. Our findings provide novel direct experimental evidence that Hes1 plays an essential role in the development of ICC via PDC.Significance:This study contributes to the identification of the cells of origin that initiate ICC and suggests that HES1 may represent a therapeutic target in ICC.

Pancreatic neuroendocrine tumors (PanNET) were classified into grades (G) 1 to 3 by the World Health Organization in 2017, but the precise mechanisms of PanNET initiation and progression have remained unclear. In this study, we used a genetically engineered mouse model to investigate the mechanisms of PanNET formation. Although pancreas-specific deletion of the Rb gene (Pdx1-Cre;Rbf/f) in mice did not affect pancreatic exocrine cells, the α-cell/β-cell ratio of islet cells was decreased at 8 months of age. During long-term observation (18–20 months), mice formed well-differentiated PanNET with a Ki67-labeling index of 2.7%. In contrast, pancreas-specific induction of a p53 mutation (Pdx1-Cre;Trp53R172H) had no effect on pancreatic exocrine and endocrine tissues, but simultaneous induction of a p53 mutation with Rb gene deletion (Pdx1-Cre;Trp53R172H;Rb f/f) resulted in the formation of aggressive PanNET with a Ki67-labeling index of 24.7% over the short-term (4 months). In Pdx1-Cre;Trp53R172H;Rbf/f mice, mRNA expression of Pten and Tsc2, negative regulators of the mTOR pathway, significantly decreased in the islet cells, and activation of the mTOR pathway was confirmed in subsequently formed PanNET. Thus, by manipulating Rb and p53 genes, we established a multistep progression model from dysplastic islet to indolent PanNET and aggressive metastatic PanNET in mice. These observations suggest that Rb and p53 have distinct roles in the development of PanNET.Significance:Pancreas-specific manipulation of Rb and p53 genes induced malignant transformation of islet cells, reproducing stepwise progression from microadenomas to indolent (grade 1) and subsequent aggressive PanNETs (grade 2–3).

Supplementary Figure S5 and legend

Supplementary Data from Expansion of Gastric Intestinal Metaplasia with Copy Number Aberrations Contributes to Field Cancerization

Supplementary Table from Expansion of Gastric Intestinal Metaplasia with Copy Number Aberrations Contributes to Field Cancerization

Supplemental Data

Supplementary Tables S1-S13. Supplementary Table S1. Clinical information and the summary of experiments performed on each sample. Supplementary Table S2. List of genes and regions captured for targeted deep sequencing. Supplementary Table S3. dbSNP ID list for which target baits were designed. Supplementary Table S4. Non-silent mutations detected in a validation cohort by targeted deep sequencing. Supplementary Table S5. Validated non-silent somatic mutations by PCR-based deep sequencing. Supplementary Table S6. Validated silent/non-silent somatic mutations used for clonal analysis of multiple samples from PBL001. Supplementary Table S7. Significantly hypo-/hyper-methylated gene sets in PBL compared to normal pancreas. Supplementary Table S8. Fusion genes detected by whole transcriptome sequencing in PBL. Supplementary Table S9. Genes with a significant differential expression between PBL and normal pancreatic samples. Supplementary Table S10. Significant gene set enrichment in comparison between PBL and normal pancreas. Supplementary Table S11. List of differentially expressed genes in each gene cluster identified by hierarchical clustering. Supplementary Table S12. Top 20 significantly enriched pathways in each gene cluster. Supplementary Table S13. Gene coefficients of PC1 and PC2.

Pancreatoblastoma is a rare pediatric pancreatic malignancy for which the molecular pathogenesis is not understood. In this study, we report the findings of an integrated multiomics study of whole-exome and RNA sequencing as well as genome-wide copy number and methylation analyses of ten pancreatoblastoma cases. The pancreatoblastoma genome was characterized by a high frequency of aberrant activation of the Wnt signaling pathway, either via somatic mutations of CTNNB1 (90%) and copy-neutral loss of heterozygosity (CN-LOH) of APC (10%). In addition, imprinting dysregulation of IGF2 as a consequence of CN-LOH (80%), gain of paternal allele (10%), and gain of methylation (10%) was universally detected. At the transcriptome level, pancreatoblastoma exhibited an expression profile characteristic of early pancreas progenitor-like cells along with upregulation of the R-spondin/LGR5/RNF43 module. Our results offer a comprehensive description of the molecular basis for pancreatoblastoma and highlight rational therapeutic targets for its treatment.Significance: Molecular genetic analysis of a rare untreatable pediatric tumor reveals Wnt/IGF2 aberrations and features of early pancreas progenitor-like cells, suggesting cellular origins and rational strategies for therapeutic targeting. Cancer Res; 78(4); 865–76. ©2017 AACR.

Supplementary Table S1

Supplemental Figure legends and Table Legend

Purpose:Extramammary Paget disease (EMPD) is an uncommon skin malignancy whose genetic alterations are poorly characterized. Previous reports identified mutations in chromatin remodeling genes and PIK3CA. In order to unambiguously determine driver mutations in EMPD, we analyzed 87 EMPD samples using exome sequencing in combination with targeted sequencing.Experimental Design:First, we analyzed 37 EMPD samples that were surgically resected using whole-exome sequencing. Based on several in silico analysis, we built a custom capture panel of putative driver genes and analyzed 50 additional formalin-fixed, paraffin-embedded samples using target sequencing. ERBB2 expression was evaluated by HER2 immunohisotochemistry. Select samples were further analyzed by fluorescence in situ hybridization.Results:A median of 92 mutations/sample was identified in exome analysis. A union of driver detection algorithms identified ERBB2, ERBB3, KMT2C, TP53, PIK3CA, NUP93, AFDN, and CUX1 as likely driver mutations. Copy-number alteration analysis showed regions spanning CDKN2A as recurrently deleted, and ERBB2 as recurrently amplified. ERBB2, ERBB3, and FGFR1 amplification/mutation showed tendency toward mutual exclusivity. Copy-number alteration load was associated with likelihood to recur. Mutational signatures were dominated by aging and APOBEC activation and lacked evidence of ultraviolet radiation. HER2 IHC/fluorescence in situ analysis validated ERBB2 amplification but was underpowered to detect mutations. Tumor heterogeneity in terms of ERBB2 amplification status was observed in some cases.Conclusions:Our comprehensive, unbiased analysis shows EMPD is characterized by alterations involving the PI3K–AKT pathway. EMPD is distinct from other skin cancers in both molecular pathways altered and etiology behind mutagenesis.

Supplementary Methods

Supplementary Figures S1-S10. Supplementary Figure S1. Depths and coverages of targeted deep sequencing. Supplementary Figure S2. Sequencing coverages and numbers of mutations detected by WES. Supplementary Figure S3. Tumor subclonal populations in multiple samples from PBL001 estimated by PyClone. Supplementary Figure S4. Mutational signatures of primary and metastatic PBL. Supplementary Figure S5. SNP array-based copy number analysis of PBL. Supplementary Figure S6. Targeted deep sequencing based allele-specific copy number analysis of chromosome 11. Supplementary Figure S7. Bisulfite Sanger sequencing of ICR1 on chromosome 11p15.5. Supplementary Figure S8. Relative expression levels of IGF2 in CN-LOH-negative samples. Supplementary Figure S9. Aberrant activation of Wnt signaling pathway in PBL009. Supplementary Figure S10. Global DNA methylation profile of PBL.

Azacitidine is a mainstay of therapy for MDS-related diseases. The purpose of our study is to elucidate the effect of gene mutations on hematological response and overall survival (OS), particularly focusing on their post-treatment clone size. We enrolled a total of 449 patients with MDS or related myeloid neoplasms. They were analyzed for gene mutations in pre- (n=449) and post- (n=289) treatment bone marrow samples using targeted-capture sequencing to assess the impact of gene mutations and their post-treatment clone size on treatment outcomes. In Cox proportional hazard modeling, multi-hit TP53 mutation (HR, 2.03; 95% CI, 1.42-2.91; P

Intestinal metaplasia (IM) is a risk factor for gastric cancer following infection with Helicobacter pylori. To explore the susceptibility of pure gastric IM to cancer development, we investigated genetic alterations in single IM gastric glands. We isolated 50 single IM or non-IM glands from the inflamed gastric mucosa of 11 patients with intramucosal gastric carcinoma (IGC) and 4 patients without IGC; 19 single glands in the noninflamed gastric mucosa of 11 individuals from our cohort and previous dataset were also included as controls. Whole-exome sequencing of single glands revealed significantly higher accumulation of somatic mutations in various genes within IM glands compared with non-IM glands. Clonal ordering analysis showed that IM glands expanded to form clusters with shared mutations. In addition, targeted-capture deep sequencing and copy number (CN) analyses were performed in 96 clustered IM or non-IM gastric glands from 26 patients with IGC. CN analyses were also performed on 41 IGC samples and The Cancer Genome Atlas-Stomach Adenocarcinoma datasets. These analyses revealed that polyclonally expanded IM commonly acquired CN aberrations (CNA), including amplification of chromosomes 8, 20, and 2. A large portion of clustered IM glands typically consisted of common CNAs rather than other cancer-related mutations. Moreover, the CNA patterns of clustered IM glands were similar to those of IGC, indicative of precancerous conditions. Taken together, these findings suggest that, in the gastric mucosa inflamed with H. pylori infection, IM glands expand via acquisition of CNAs comparable with those of IGC, contributing to field cancerization. Significance: This study contributes to our understanding of gastric intestinal metaplasia as a risk factor for gastric adenocarcinoma via their multifocal expansion and acquisition of CNAs and somatic mutations.

Adult T-cell leukemia/lymphoma (ATL) is an aggressive neoplasm immunophenotypically resembling regulatory T cells, associated with human T-cell leukemia virus type-1. Here, we performed whole-genome sequencing (WGS) of 150 ATL cases to reveal the overarching landscape of genetic alterations in ATL. We discovered frequent (33%) loss-of-function alterations preferentially targeting the CIC long isoform, which were overlooked by previous exome-centric studies of various cancer types. Long but not short isoform–specific inactivation of Cic selectively increased CD4+CD25+Foxp3+ T cells in vivo. We also found recurrent (13%) 3′-truncations of REL, which induce transcriptional upregulation and generate gain-of-function proteins. More importantly, REL truncations are also common in diffuse large B-cell lymphoma, especially in germinal center B-cell–like subtype (12%). In the non-coding genome, we identified recurrent mutations in regulatory elements, particularly splice sites, of several driver genes. In addition, we characterized the different mutational processes operative in clustered hypermutation sites within and outside immunoglobulin/T-cell receptor genes and identified the mutational enrichment at the binding sites of host and viral transcription factors, suggesting their activities in ATL. By combining the analyses for coding and noncoding mutations, structural variations, and copy number alterations, we discovered 56 recurrently altered driver genes, including 11 novel ones. Finally, ATL cases were classified into 2 molecular groups with distinct clinical and genetic characteristics based on the driver alteration profile. Our findings not only help to improve diagnostic and therapeutic strategies in ATL, but also provide insights into T-cell biology and have implications for genome-wide cancer driver discovery.

Germ line DDX41 variants have been implicated in late-onset myeloid neoplasms (MNs). Despite an increasing number of publications, many important features of DDX41-mutated MNs remain to be elucidated. Here we performed a comprehensive characterization of DDX41-mutated MNs, enrolling a total of 346 patients with DDX41 pathogenic/likely-pathogenic (P/LP) germ line variants and/or somatic mutations from 9082 MN patients, together with 525 first-degree relatives of DDX41-mutated and wild-type (WT) patients. P/LP DDX41 germ line variants explained ∼80% of known germ line predisposition to MNs in adults. These risk variants were 10-fold more enriched in Japanese MN cases (n = 4461) compared with the general population of Japan (n = 20 238). This enrichment of DDX41 risk alleles was much more prominent in male than female (20.7 vs 5.0). P/LP DDX41 variants conferred a large risk of developing MNs, which was negligible until 40 years of age but rapidly increased to 49% by 90 years of age. Patients with myelodysplastic syndromes (MDS) along with a DDX41-mutation rapidly progressed to acute myeloid leukemia (AML), which was however, confined to those having truncating variants. Comutation patterns at diagnosis and at progression to AML were substantially different between DDX41-mutated and WT cases, in which none of the comutations affected clinical outcomes. Even TP53 mutations made no exceptions and their dismal effect, including multihit allelic status, on survival was almost completely mitigated by the presence of DDX41 mutations. Finally, outcomes were not affected by the conventional risk stratifications including the revised/molecular International Prognostic Scoring System. Our findings establish that MDS with DDX41-mutation defines a unique subtype of MNs that is distinct from other MNs.

Esophageal squamous cell carcinoma (ESCC) often recurs after chemoradiotherapy, and the prognosis of ESCC after chemoradiotherapy has not improved over the past few decades. The mutation process in chemoradiotherapy-resistant clones and the functional relevance of genetic alterations remain unclear. To address these problems, we performed whole-exome sequencing of 52 tumor samples from 33 patients with ESCC who received radiotherapy combined with 5-fluorouracil/platinum. In multiregion analyses of pretreatment and locally recurrent lesions from five cases, most driver gene-altered clones remained under chemoradiotherapy selection pressure, while few driver gene alterations were acquired at recurrence. The mutation signatures of recurrent ESCC, including increased deletion frequency and platinum dose-dependent base substitution signatures, were substantially different from those of primary ESCC and reflected the iatrogenic impacts of chemoradiotherapy. Single-region analysis of 28 pretreatment tumors indicated that focal copy-number gain at the MYC locus was significantly associated with poor progression-free survival and overall survival after chemoradiotherapy. MYC gain remained throughout the chemoradiotherapy course and potentially contributes to intrinsic resistance to chemoradiotherapy. Consistent with these findings, MYC copy number and mRNA and protein levels in ESCC cell lines correlated positively with resistance to radiotherapy, and MYC knockdown improved sensitivity to radiotherapy. Overall, these data characterize the clonal evolution process induced by chemoradiotherapy and clinically relevant associations for genetic alterations in ESCC. These findings increase our understanding of therapeutic resistance and support the rationale for precision chemoradiotherapy. Significance: Whole-exome sequencing reveals the genetic evolution of ESCC during chemoradiotherapy, highlighting MYC gain in pretreatment tumors as a potential marker of therapy resistance.

Notch/Hes1 signaling has been shown to play a role in determining the fate of pancreatic progenitor cells. However, its function in postnatal pancreatic maturation is not fully elucidated. We generated conditional Hes1 knockout and/or Notch intracellular domain (NICD) overexpression mice in Ptf1a- or Pdx1-positive pancreatic progenitor cells and analyzed pancreatic tissues. Both Ptf1acre/+; Hes1f/f and Ptf1acre/+; Rosa26NICD mice showed normal pancreatic development at P0. However, exocrine tissue of the pancreatic tail in Ptf1acre/+; Hes1f/f mice atrophied and was replaced by fat tissue by 4 weeks of age, with increased apoptotic cells and fewer centroacinar cells. This impaired exocrine development was completely rescued by NICD overexpression in Ptf1acre/+; Hes1f/f; Rosa26NICD mice, suggesting compensation by a Notch signaling pathway other than Hes1. Conversely, Pdx1-Cre; Hes1f/f mice showed impaired postnatal exocrine development in both the pancreatic head and tail, revealing that the timing and distribution of embryonic Hes1 expression affects postnatal exocrine tissue development. Notch signaling has an essential role in pancreatic progenitor cells for the postnatal maturation of exocrine tissue, partly through the formation of centroacinar cells.

In Helicobacter pylori (Hp)-uninfected individuals, diffuse-type gastric cancer (DGC) was reported as the most common type of cancer. However, the carcinogenic mechanism of Hp-uninfected sporadic DGC is largely unknown. We performed whole-exome sequencing of Hp-uninfected DGCs and Hp-uninfected normal gastric mucosa. For advanced DGCs, external datasets were also analyzed. Eighteen patients (aged 29–78 years) with DGCs and nine normal subjects (28–77 years) were examined. The mutation burden in intramucosal DGCs (10–66 mutations per exome) from individuals aged 29–73 years was not very different from that in the normal gastric glands, which showed a constant mutation accumulation rate (0.33 mutations/exome/year). Unbiased dN/dS analysis showed that CDH1 somatic mutation was a driver mutation for intramucosal DGC. CDH1 mutation was more frequent in intramucosal DGCs (67%) than in advanced DGCs (27%). In contrast, TP53 mutation was more frequent in advanced DGCs (52%) than in intramucosal DGCs (0%). This discrepancy in mutations suggests that CDH1-mutated intramucosal DGCs make a relatively small contribution to advanced DGC formation. Among the 16 intramucosal DGCs (median size, 6.5 mm), 15 DGCs were pure signet ring cell carcinoma (SRCC) with reduced E-cadherin expression and a low proliferative capacity (median Ki-67 index, 2.4%). Five SRCCs reviewed endoscopically over 2–5 years showed no progression. Impaired E-cadherin function due to CDH1 mutation was considered as an early carcinogenic event of Hp-uninfected intramucosal SRCC. Genetic and clinical analyses suggest that Hp-uninfected intramucosal SRCCs may be less likely to develop into advanced DGCs.

Background Clonal hematopoiesis (CH) in apparently healthy individuals has been drawing an increasing attention of recent years due to its implication in the risk of hematopoietic malignancies and also in cardiovascular complications. However, our knowledge on CH has been largely based on genetic studies, while few functional analyses have been performed using human materials with most studies being confined to artificial models using mice. The major challenge here is the difficulty of isolating CH clones from wild-type (WT) cells in human bone marrow (BM) to elucidate the effect of CH-mutations. Methods To investigate cellular phenotypes of mutated and WT cells separately in CH, we developed a Fluidigm C1-based single cell-sequencing platform for simultaneous genotyping and gene expression analysis. We analyzed a total of 10,178 hematopoietic stem/progenitors (HSPCs) derived from BM of patients with (n=11) and without (n=17) CH. Results In the analysis of HSPCs from CH(−) elderly individuals, we found a significant positive correlation between age and expression of gene sets implicated in inflammatory responses including TGFβ signalling and IL-6. We next analysed CH(+) samples, including those with mutations in TET2, DNMT3A, SF3B1 and IDH1 in which we investigated the pathways enriched in differentially expressed genes between mutated and unmutated cells. Regardless of mutation type, mutant cells showed an upregulation of genes implicated in an enhanced cell proliferation, while genes related to inflammatory responses were significantly downregulated. These results suggest that mutant HSPCs show an enhanced cell proliferation and an attenuated response to an inflammatory microenvironment in aged BM, compared with endogenous WT counterparts. Based on these observations, we further investigated the role of the BM microenvironment in CH, in which we compared the phenotype of WT cells in age-matched CH(+) and CH(−) cases. Compared with those from CH(−) cases, WT cells from these CH(+) cases showed an enhanced response to proinflammatory cytokines, including IL-6, interferons and TNF, suggesting a possibility that CH(+) BM might be characterized by a more enhanced inflammatory microenvironment, compared with CH(−) BM. To understand the different phenotype of WT cells between CH(+) and CH(−) cases, we next investigated a possible effect of mutant cells on endogenous WT cells. Such an effect was first documented in an IDH1-mutated CH case. IDH1-mutated cells showed a gene expression profile suggestive of an enhanced cell proliferation compared with endogenous IDH1-WT cells (WT IDH1). However, on the basis of comparison of gene expression between WT cells from CH(−) cases, it was better explained by suppressed proliferation of endogenous WT cells, rather than enhanced proliferation of mutant cells. In agreement with this, mouse BM cells treated with 2-hydroxyglutalate (2HG), an oncometabolite produced by IDH1-mutant cells, mimicked the endogenous WT cells in IDH1-mutated cases, including downregulated E2F target genes and upregulated inflammation-related genes, compared with control BM cells. A non-cell autonomous effect of mutations was also seen in cases with TET2-mutated CH, in which not only TET2-mutated cells but also endogenous WT cells exhibited a significantly different gene expression profile, compared with those from CH(−) cases. Of interest, when BM cells from WT Ly5.1/5.2 mice were co-transplanted with those from heterozygous conditional Tet2 knock-out (Ly5.2) or WT (Ly5.2) mice into lethally irradiated mice (Ly5.1) and flow-sorted Lin − WT donor competitor cells were analysed using single cell sequencing, the WT HSC-like cells co-transplanted with Tet2-mutant competitors exhibited enhanced cell proliferation and IFNα and IFNγ pathway genes, compared with those co-transplanted with WT competitors. Conclusions Taken together, these results suggest that mutant cells in CH(+) BM have non-cell autonomous effects on endogenous WT cells, which might be responsible for an accelerated inflammatory microenvironment of aged BM, favor positive selection of CH-clones, and also affect the phenotype of endogenous WT cells, contributing to the pathogenesis of CH. Disclosures Nakagawa: Sumitomo Dainippon Pharma Oncology, Inc.: Research Funding. Inagaki: Sumitomo Dainippon Pharma Oncology, Inc.: Current Employment. Nannya: Otsuka Pharmaceutical Co., Ltd.: Consultancy, Speakers Bureau; Astellas: Speakers Bureau. Yoda: Chordia Therapeutics Inc.: Research Funding. Ogawa: Kan Research Laboratory, Inc.: Consultancy, Research Funding; Dainippon-Sumitomo Pharmaceutical, Inc.: Research Funding; ChordiaTherapeutics, Inc.: Consultancy, Research Funding; Ashahi Genomics: Current holder of individual stocks in a privately-held company; Otsuka Pharmaceutical Co., Ltd.: Research Funding; Eisai Co., Ltd.: Research Funding.

A 39-year-old man had an intracranial tumour without infiltration into the surrounding cerebral tissue. The tumour recurred seven times in 11 years but maintained a well-demarcated character. Histopathological examination of the 4th surgical specimens showed nests of tumour cells surrounding small blood vessels. The tumour cells harboured amphophilic cytoplasm and small round nuclei with fine chromatin, and perinuclear haloes and clear borders were frequently observed, which was unclassifiable histology. By the Deutsches Krebsforschungszentrum methylation classifier, the tumour was not classified into any of the methylation classes. mRNA sequencing identified a novel SREBF1::NACC1 gene fusion. This intracranial tumour could be a novel tumour entity with NACC1 rearrangement showing characteristic histological and diagnostic imaging findings.

Context Pheochromocytoma and paraganglioma (PPGL) may appear as a complication of cyanotic congenital heart disease (CCHD-PPGL) with frequent EPAS1 mutations, suggesting a close link between EPAS1 mutations and tissue hypoxia in CCHD-PPGL pathogenesis. Objective Our aim is to further investigate the role of EPAS1 mutations in the hypoxia-driven mechanism of CCHD-PPGL pathogenesis, particularly focusing on metachronous and/or multifocal CCHD-PPGL tumors. Methods We performed whole-exome sequencing (WES) for somatic and germline mutations in 15 PPGL samples from 7 CCHD patients, including 3 patients with metachronous and/or multifocal tumors, together with an adrenal medullary hyperplasia (AMH) sample. Results We detected EPAS1 mutations in 15 out of 16 PPGL/AMH samples from 7 cases. Conspicuously, all EPAS1 mutations in each of 3 cases with multifocal or metachronous tumors were mutually independent and typical examples of parallel evolution, which is suggestive of strong positive selection of EPAS1-mutated clones. Compared to 165 The Cancer Genome Atlas non–CCHD-PPGL samples, CCHD-PPGL/AMH samples were enriched for 11p deletions (13/16) and 2p amplifications (4/16). Of particular note, the multiple metachronous PPGL tumors with additional copy number abnormalities developed 18 to 23 years after the resolution of hypoxemia, suggesting that CCHD-induced hypoxic environments are critical for positive selection of EPAS1 mutants in early life, but may no longer be required for development of PPGL in later life. Conclusion Our results highlight a key role of activated hypoxia-inducible factor 2α due to mutated EPAS1 in positive selection under hypoxic environments, although hypoxemia itself may not necessarily be required for the EPAS1-mutated clones to progress to PPGL.

Cancer is a clonal disorder derived from a single ancestor cell and its progenies that are positively selected by acquisition of ‘driver mutations’. However, the evolution of positively selected clones does not necessarily imply the presence of cancer. On the contrary, it has become clear that expansion of these clones in phenotypically normal or non-cancer tissues is commonly seen in association with ageing and/or in response to environmental insults and chronic inflammation. Recent studies have reported expansion of clones harbouring mutations in cancer driver genes in the blood, skin, oesophagus, bronchus, liver, endometrium and bladder, where the expansion could be so extensive that tissues undergo remodelling of an almost entire tissue. The presence of common cancer driver mutations in normal tissues suggests a strong link to cancer development, providing an opportunity to understand early carcinogenic processes. Nevertheless, some driver mutations are unique to normal tissues or have a mutation frequency that is much higher in normal tissue than in cancer, indicating that the respective clones may not necessarily be destined for evolution to cancer but even negatively selected for carcinogenesis depending on the mutated gene. Moreover, tissues that are remodelled by genetically altered clones might define functionalities of aged tissues or modified inflammatory processes. In this Review, we provide an overview of major findings on clonal expansion in phenotypically normal or non-cancer tissues and discuss their biological significance not only in cancer development but also in ageing and inflammatory diseases. Clonal expansion in phenotypically normal or non-cancer tissues is commonly seen in association with ageing and/or in response to environmental insults and chronic inflammation, but does not necessarily indicate cancer development. This Review discusses recent findings on clonal expansion in these tissues and their biological significance in cancer development, ageing and inflammatory diseases.

Background: Chronic kidney disease is frequently associated with persistent inflammation, which results in fibrosis in the stroma, a reduced number of renal tubules, the formation of multiple cystic lesions, ultimately terminating in renal failure and hemodialysis. Given the high incidence of renal cell carcinoma (RCC) in hemodialysis patients, suggesting a relationship between tissue remodeling by chronic inflammation and carcinogenesis. However, little is known about the genetic background of cancer development from remaining tubules and cystic lesions in hemodialysis patients. Method: We enrolled 5 patients under hemodialysis who were accompanied by acquired cystic kidney disease and underwent radical nephrectomy for RCC. Surgical specimens were fixed with alcohol-based solution and paraffin-embedded, and after H&E staining, subjected to laser capture microdissection (LCM) to collect remaining tubules and cysts containing approximately 200 cells. DNA was extracted and analyzed for somatic mutations and copy number alteration by whole-exome sequencing. Result: In total, we collected 161 LCM samples, including 118 from proximal tubules, 17 from collecting ducts, and 26 from cysts. Median variant allele frequencies of detected mutations were 0.237 in proximal tubules, 0.133 in collecting duct, and 0.381 in cysts, indicating larger clonal expansion in proximal tubules and cysts than in collecting ducts. Proximal tubules and cysts contained recurrent mutations in FAT1 (11% and 8%, respectively), STAG2 (5% and 8%), and PTEN (3% and 4%), whereas no recurrent driver mutations were identified in collecting ducts. In copy number analysis, 26% of proximal tubules, 35% of collecting ducts, and 80% of cysts had copy number alterations. Of note, samples from cysts had more copy number alterations compared to remaining tubules (on average, 2.3 vs. 1.1). Proximal tubules had recurrent copy number gains of chromosomes 3, 7, 10, 18, and 20, and loss of chromosome 22 and samples from collecting ducts showed gains of chromosome 7 and loss of chromosome 18. In cysts, gains of chromosomes 2, 3, 7, 10, 12, and 16 and loss of chromosomes 15, 16, 21, and 22 were frequently observed. Copy number profile in cysts was similar to that of papillary RCC and acquired cystic disease-associated RCC. Conclusion: In the end-stage cystic kidney, proximal tubules and cysts showed an enrichment of driver mutations commonly found in papillary RCC and RCC with acquired cystic disease, reflecting the fact that the incidence of these two types of RCC increases as the duration of dialysis becomes longer. In addition, clear cell RCC drivers, such as VHL mutation and loss of chromosome 3, were not observed, possibly explaining the rare occurrence of clear cell RCC in hemodialysis patients. Since cysts had more frequent copy number alterations than remaining tubules, cysts in the end-stage kidney might be precursor lesions of RCC. Citation Format: Kosuke Ieiri, Nobuyuki Kakiuchi, Tomonori Hirano, Tomomi Nishimura, Koichi Watanabe, Hiroko Tanaka, Satoru Miyano, Dai Takamatsu, Keisuke Monji, Eiji Kashiwagi, Masaki Shiota, Junichi Inokuchi, Masatoshi Eto, Seishi Ogawa. Genomic analysis of end-stage renal disease [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 6066.

Introduction: The majority of patients with pancreatic cancer are inoperable at the time of diagnosis and even those cases with successful tumor resection for early-stage tumors frequently had second tumors within the remaining pancreas shortly after the surgery, contributing to a poor clinical outcome. Moreover, synchronous multiple pancreatic cancers were often reported. One of the long-standing issues about such synchronous and/or metachronous multiple pancreatic cancers is their clonal origins. In this study, we aimed to reveal the origin of multiple pancreatic cancers using an unbiased detection of somatic mutations in primary and metachronous cancers as well as adjacent precursor lesions. Methods: Serially obtained formalin-fixed paraffin-embedded surgical specimens from 18 patients who had undergone curative surgery for an early-stage pancreatic cancer were subjected to laser microdissection for the enrichment of tumor and precancerous lesions, from which DNA was extracted and analyzed for somatic mutations using whole-exome sequencing (WES) with matched normal DNA. Based on shared and private mutations across different samples, we interrogated history of clonal evolution of these lesions. Results: We analyzed 18 patients with multiple pancreatic cancers of 2 metachronous and synchronous, 14 metachronous, and 2 synchronous tumors. In metachronous cases, pathology for primary cancer were margin-negative in all patients, and the median interval between the initial and second surgery was 38.9 months (7 - 85 months). In addition, a total of 20 pancreatic intraepithelial neoplasia (PanIN) from 5 of these cases were analyzed. We identified a median of 66 (range: 31-232) and 20 (14-42) somatic mutations in cancer and PanIN lesions, respectively. None of the patients have known pathogenic germline variants. All samples had one or more driver mutations. In most cases (N=17), sharing many somatic mutations, multiple tumors had a common clonal origin. In another synchronous case, the synchronous tumors shared only one mutation, suggesting that they were clonally independent tumors. By contrast, none of the mutations other than hotspot KRAS mutations were shared between precursor and cancer lesions. While multiple independent clones were observed in precancerous lesions, most of multiple cancers originated from a common ancestor. Moreover, negative pathology of the margins at the initial surgery suggested that those multiple lesions arose from distant dissemination or metastasis, rather than contiguous, intraductal invasion. Conclusions: In conclusion, our study suggests that even early pancreatic cancer might be disseminated within the pancreas and contribute to synchronous and metachronous cancers. Citation Format: Tomonori Hirano, Nobuyuki Kakiuchi, Yasuhide Takeuchi, Tomomi Nishimura, Toshihiko Masui, Kazuyuki Nagai, Takayuki Anazawa, Sachiko Minamigushi, Hironori Haga, Norimitsu Uza, Hiroshi Seno, Yuzo Kodama, Atsuhiro Masuda, Takeshi Tanaka, Seishi Ogawa, Hiroko Tanaka, Satoru Miyano. Origin of synchronous or metachronous multiple pancreatic cancers [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 1511.

Introduction: The epithelial cells of the endometrium undergo rapid proliferation and shedding each month under the influence of estrogen and progesterone. Recent studies have demonstrated pervasive driver mutations in the endometrium with ageing. However, previous observations have been limited to the microscopic scale, and it is still unclear how clonal expansion occurs over the entire endometrium and what driver mutations are responsible therein. To address these issues, we performed multisampling of the endometrium with extensive spatial mapping. Methods: We collected multiple endometrial glands from normal endometrium in a 5mm lattice in 19 patients who underwent hysterectomy (3-31 samples per case), followed by DNA extraction, detected somatic mutations by whole exome sequencing, to identify large clones that spread across samples based on shared mutations. Results: In total, 216 bulk endometrial glands were analyzed. The number of detected somatic mutations per sample ranged from 6 to 98 (mean±SD: 37.8±16.0), and most of the samples had at least one known driver gene mutation (mean±SD: 4.38±2.27). The mutational signature analysis revealed the predominance of the age-related SBS1 and SBS5 signatures. dN/dS analysis identified 21 genes that were positively selected in normal endometrium, including PIK3CA, PIK3R1, KRAS, PPP2R1A, ARHGAP35, and FBXW7, most of which are also known as driver genes for endometrial cancer. The number of mutations and their frequency for each driver gene did not differ by background diseases: endometriosis, endometriosis-associated ovarian cancer and other non-endometrial related diseases. We found clones showing a large expansion involving multiple samples in 8 cases, most of which harbored one or more driver mutations, such as PIK3CA and KRAS. The size of the largest clone carrying a KRAS mutation, which was found in a 32-year-old woman with cervical cancer, spanned as long as 10 mm in diameter. This suggests that clonal expansion in the endometrium could be very rapid even in healthy women, which contrasts to the observations that clonal expansion in other organs is usually seen in aged individuals. Conclusion: We confirmed the high frequency of driver mutations in the normal endometrium, most of which overlapped to those found in endometrial cancers, suggesting that these driver genes are involved in the early phase of the development in the endometrial cancer. Some of the driver-mutated clones expanded in macroscopic size. Our findings provide an interesting insight into the clonal structure and underlying genetic variation in the normal endometrium. Citation Format: Koichi Watanabe, Nobuyuki Kakiuchi, Shiro Takamatsu, Sachiko Kitamura, Mana Taki, Koji Yamanoi, Ryusuke Murakami, Ken Yamaguchi, Junzo Hamanishi, Hiroko Tanaka, Satoru Miyano, Masaki Mandai, Seishi Ogawa. Clonal expansion with driver mutations in normal human endometrium [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 128.

Introduction: Adjuvant chemotherapy (ACT) has a key role in the treatment of high-risk colorectal cancer (CRC) patients to prevent post surgical recurrence. However, the efficacy of ACT is offset by the development of complication such as peripheral neuropathy, which leads to compromised quality of life, while ACT is expected to prevent tumor recurrence in only 5% and 20% in Stage2 and 3 patients, respectively. Thus, there is an urgent need to develop a biomarker that accurately predicts the effect of ACT and therefore helps guide ACT eligibility. Methods: We enrolled a total of 3,023 cases with CRC, for whom targeted-capture sequencing was performed that covered 169 known or putative driver genes and 1,660 polymorphoic sites to detect somatic mutations and copy number alterations (CNAs), respectively. RNA sequencing was also performed to detect driver gene fusions and to profile gene expression. We identified, through univariate analysis, gene mutations/CNAs as well as clinical parameters using those cases who did not received ACT to establish a model that predict time to recurrence (TTR) using the Cox proportional hazard modelling. Thereafter, the model was applied to those who did receive ACT to assess the impact of the model in the prediction of TTR. Results: We identified a total of 32145 mutations and/or 44511 CNAs in 2970 (98%) out of 3,023 cases. Mutations were found in 55 driver genes, including 15 previously unknown drivers, such as MAP2K1, AKT1, and TCF7. According to the number of mutations, 278 casese were classified as hypermutated cases. The remaining ‘non-hypermutated cases’ (n=2692) were classified into APC-mutated, RNF43/ZNRF3-mutated, CTNNB1-mutated, RSPO-fusion(+) and other Wnt-mutation(−) cases. The TTR prediction model was constructed by repeating Cox regression 1000 times using randomly selected 80% of the non-hypermutated cases and variables identified in univariate analysis. Adopting those variables selected in >70% models, we established the final model including both genetic (mutations/CNAs) and clinical variables, based on which we classified patients into 4 risk groups, low/intermediate/high/very-high. Finally, by applying the classification to those who received ACT, we demonstrated that 18% of Stage2 and the majority of Stage3 cases with high/very-high risks were predicted to benefit from ACT, whereas most of the Stage2 and as many as 19% of Stage3 cases received ACT with no significant improvement of TTR. Approximately 5% of Stage1 cases belonged to the high-risk category and potentially benefit from ACT. Conclusions: On the basis of large-scale genotyping of clinically well-annotated CRC cases, we proposed molecular classification of CRC and established a model that enabled prediction of TTR and risk stratification, which was successfully used to identify cases who benefited (within Stage 2/3 cases) or will benefit (in Stage1 cases) from ACT to prolong TTR. Citation Format: Yoshikage Inoue, Nobuyuki Kakiuchi, Kenichi Yoshida, Yasuhito Nanya, Yasuhide Takeuchi, Yuichi Shiraishi, Kenichi Chiba, Tetusichi Yoshizato, Hiroko Tanaka, Satoshi Nagayama, Kazutaka Obama, Satoru Miyano, Seishi Ogawa. Genomic characterization and risk stratification of colorectal cancer [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 978.

Intrahepatic cholangiocarcinoma (ICC) is frequently driven by aberrant KRAS activation and develops in the liver with chronic inflammation. Although the Notch signaling pathway is critically involved in ICC development, detailed mechanisms of Notch-driven ICC development are still unknown. Here, we use mice whose Notch signaling is genetically engineered to show that the Notch signaling pathway, specifically the Notch/Hes1 axis, plays an essential role in expanding ductular cells in the liver with chronic inflammation or oncogenic Kras activation. Activation of Notch1 enhanced the development of proliferating ductal cells (PDC) in injured livers, while depletion of Hes1 led to suppression. In correlation with PDC expansion, ICC development was also regulated by the Notch/Hes1 axis and suppressed by Hes1 depletion. Lineage-tracing experiments using EpcamcreERT2 mice further confirmed that Hes1 plays a critical role in the induction of PDC and that ICC could originate from PDC. Analysis of human ICC specimens showed PDC in nonneoplastic background tissues, confirming HES1 expression in both PDC and ICC tumor cells. Our findings provide novel direct experimental evidence that Hes1 plays an essential role in the development of ICC via PDC. Significance: This study contributes to the identification of the cells of origin that initiate ICC and suggests that HES1 may represent a therapeutic target in ICC.

Background Accurate preoperative assessment of the longitudinal extension of perihilar cholangiocarcinoma (PHCC) is essential for treatment planning. Mapping biopsies for PHCC remain challenging owing to technical difficulties and insufficient sample amounts. The aim of this study was to investigate the usefulness of a novel technique for mapping biopsies of PHCC. Methods Our novel method focused on a biliary stent delivery system for mapping biopsies. Fifty patients with PHCC undergoing endoscopic transpapillary mapping biopsy using the novel method were reviewed from August 2015 to June 2019. Results The median number of biopsy samples was six (range 1 – 17), and the rate of adequate sampling was 91.4 % (266 /291). Biopsy from the intrahepatic bile duct was possible in 82.0 % of patients (41 /50), and negative margins were confirmed in the resected specimens from 34 /39 patients who underwent surgery (87.2 %). None of the patients had post-endoscopic retrograde cholangiopancreatography pancreatitis. Conclusions With our novel method, accurate assessment of the longitudinal extension of PHCC might be expected with minimal trauma to the duodenal papilla.

Tobacco smoking causes lung cancer1–3, a process that is driven by more than 60 carcinogens in cigarette smoke that directly damage and mutate DNA4,5. The profound effects of tobacco on the genome of lung cancer cells are well-documented6–10, but equivalent data for normal bronchial cells are lacking. Here we sequenced whole genomes of 632 colonies derived from single bronchial epithelial cells across 16 subjects. Tobacco smoking was the major influence on mutational burden, typically adding from 1,000 to 10,000 mutations per cell; massively increasing the variance both within and between subjects; and generating several distinct mutational signatures of substitutions and of insertions and deletions. A population of cells in individuals with a history of smoking had mutational burdens that were equivalent to those expected for people who had never smoked: these cells had less damage from tobacco-specific mutational processes, were fourfold more frequent in ex-smokers than current smokers and had considerably longer telomeres than their more-mutated counterparts. Driver mutations increased in frequency with age, affecting 4–14% of cells in middle-aged subjects who had never smoked. In current smokers, at least 25% of cells carried driver mutations and 0–6% of cells had two or even three drivers. Thus, tobacco smoking increases mutational burden, cell-to-cell heterogeneity and driver mutations, but quitting promotes replenishment of the bronchial epithelium from mitotically quiescent cells that have avoided tobacco mutagenesis. Whole-genome sequencing of normal bronchial epithelium from 16 individuals shows that tobacco smoking increases genomic heterogeneity, mutational burden and driver mutations, whereas stopping smoking promotes replenishment of the epithelium with near-normal cells.