Data availability. The WGS and RNA expression data can be found at the European Genome-phenome Archive (EGA) under accessions EGAD00001004417 and EGAD00001004423, respectively. Code availability. Code associated with the analysis is available upon request. Ethics. The study was registered (UKCRNID 8880), approved by the Institutional Ethics Committees (REC 07/H0305/52 and 10/ H0305/1), and all subjects gave individual informed consent. Reporting summary. Additional information is included in the Life Sciences Reporting Summary, which details exact software and biological materials used and efforts made to ensure reproducibility of results. Author contributions RCF and AMF conceived the overall study. AMF and SJ analyzed the genomic data and performed statistical analyses. RCF, AMF and XL designed the experiments. AMF, XL and JM performed the experiments. GC contributed to the structural variant analysis and data visualization. SK helped compile the clinical data and aided statistical analyses. JP and SA produced and QC'ed the RNA-seq data. EO aided the whole genome sequencing of EAC cell lines. SM and NG coordinated the clinical centres and were responsible for sample collections. ME benchmarked our mutation calling pipelines. MO led the pathological sample QC for sequencing. LB and GD constructed and managed the sequencing alignment and variant calling pipelines. RCF and ST supervised the research. RCF and ST obtained funding. AMF and RCF wrote the manuscript. All authors approved the manuscript.

Background: The incidence of esophageal adenocarcinoma (EAC) is rapidly rising and has a 5-year survival rate of <20%. Beyond TNM (tumorenodeemetastasis) staging, no reliable risk stratification tools exist and no large-scale studies have profiled circulating tumor DNA (ctDNA) at relapse in EAC. Here we analyze the prognostic potential of ctDNA dynamics in EAC, taking into account clonal hematopoiesis with indeterminate potential (CHIP). Patients and methods: A total of 245 samples from 97 patients treated with neoadjuvant chemotherapy and surgery were identified from the prospective national UK Oesophageal Cancer Clinical and Molecular Stratification (OCCAMS) consortium data set. A pan-cancer ctDNA panel comprising 77 genes was used. Plasma and peripheral blood cell samples were sequenced to a mean depth of 7082Â (range 2196-28 524) and ctDNA results correlated with survival. Results: Characteristics of the 97 patients identified were as follows: 83/97 (86%) male, median age 68 years (SD 9.5 years), 100% cT3/T4, 75% cNþ. EAC-specific drivers had higher variant allele fractions than passenger mutations. Using stringent quality criteria 16/79 (20%) were ctDNA positive following resection; recurrence was observed in 12/16 (75%) of these. As much as 78/97 (80%) had CHIP analyses that enabled filtering for CHIP variants, which were found in 18/78 (23%) of cases. When CHIP was excluded, 10/63 (16%) patients were ctDNA positive and 9/10 of these (90%) recurred. With correction for CHIP, median cancer-specific survival for ctDNA-positive patients was 10.0 months versus 29.9 months for ctDNA-negative patients (hazard ratio 5.55, 95% confidence interval 2.42-12.71; P ¼ 0.0003). Similar outcomes were observed for disease-free survival. Conclusions: We demonstrate in a large, national, prospectively collected data set that ctDNA in plasma following surgery for EAC is prognostic for relapse. Inclusion of peripheral blood cell samples can reduce or eliminate false positives from CHIP. In future, post-operative ctDNA could be used to risk stratify patients into high-and low-risk groups for intensification or de-escalation of adjuvant chemotherapy.

B cells are frequently found in the margins of solid tumours as organized follicles in ectopic lymphoid organs called tertiary lymphoid structures (TLS)1,2. Although TLS have been found to correlate with improved patient survival and response to immune checkpoint blockade (ICB), the underlying mechanisms of this association remain elusive1,2. Here we investigate lung-resident B cell responses in patients from the TRACERx 421 (Tracking Non-Small-Cell Lung Cancer Evolution Through Therapy) and other lung cancer cohorts, and in a recently established immunogenic mouse model for lung adenocarcinoma3. We find that both human and mouse lung adenocarcinomas elicit local germinal centre responses and tumour-binding antibodies, and further identify endogenous retrovirus (ERV) envelope glycoproteins as a dominant anti-tumour antibody target. ERV-targeting B cell responses are amplified by ICB in both humans and mice, and by targeted inhibition of KRAS(G12C) in the mouse model. ERV-reactive antibodies exert anti-tumour activity that extends survival in the mouse model, and ERV expression predicts the outcome of ICB in human lung adenocarcinoma. Finally, we find that effective immunotherapy in the mouse model requires CXCL13-dependent TLS formation. Conversely, therapeutic CXCL13 treatment potentiates anti-tumour immunity and synergizes with ICB. Our findings provide a possible mechanistic basis for the association of TLS with immunotherapy response.

Environmental carcinogenic exposures are major contributors to global disease burden yet how they promote cancer is unclear. Over 70 years ago, the concept of tumour promoting agents driving latent clones to expand was rst proposed. In support of this model, recent evidence suggests that human tissue contains a patchwork of mutant clones, some of which harbour oncogenic mutations, and many environmental carcinogens lack a clear mutational signature. We hypothesised that the environmental carcinogen, <2.5μm particulate matter (PM2.5), might promote lung cancer promotion through nonmutagenic mechanisms by acting on pre-existing mutant clones within normal tissues in patients with lung cancer who have never smoked, a disease with a high frequency of EGFR activating mutations. We analysed PM2.5 levels and cancer incidence reported by UK Biobank, Public Health England, Taiwan Chang Gung Memorial Hospital (CGMH) and Korean Samsung Medical Centre (SMC) from a total of 463,679 individuals between 2006-2018. We report associations between PM2.5 levels and the incidence of several cancers, including EGFR mutant lung cancer. We nd that pollution on a background of EGFR mutant lung epithelium promotes a progenitor-like cell state and demonstrate that PM accelerates lung cancer progression in EGFR and Kras mutant mouse lung cancer models. Through parallel exposure studies in mouse and human participants, we nd evidence that in ammatory mediators, such as interleukin-1 , may act upon EGFR mutant clones to drive expansion of progenitor cells. Ultradeep mutational pro ling of histologically normal lung tissue from 247 individuals across 3 clinical cohorts revealed oncogenic EGFR and KRAS driver mutations in 18% and 33% of normal tissue samples, respectively. These results support a tumour-promoting role for PM acting on latent mutant clones in normal lung tissue and add to evidence providing an urgent mandate to address air pollution in urban areas.

Lung cancer is the leading cause of cancer-associated mortality worldwide1. Here we analysed 1,644 tumour regions sampled at surgery or during follow-up from the first 421 patients with non-small cell lung cancer prospectively enrolled into the TRACERx study. This project aims to decipher lung cancer evolution and address the primary study endpoint: determining the relationship between intratumour heterogeneity and clinical outcome. In lung adenocarcinoma, mutations in 22 out of 40 common cancer genes were under significant subclonal selection, including classical tumour initiators such as TP53 and KRAS. We defined evolutionary dependencies between drivers, mutational processes and whole genome doubling (WGD) events. Despite patients having a history of smoking, 8% of lung adenocarcinomas lacked evidence of tobacco-induced mutagenesis. These tumours also had similar detection rates for EGFR mutations and for RET, ROS1, ALK and MET oncogenic isoforms compared with tumours in never-smokers, which suggests that they have a similar aetiology and pathogenesis. Large subclonal expansions were associated with positive subclonal selection. Patients with tumours harbouring recent subclonal expansions, on the terminus of a phylogenetic branch, had significantly shorter disease-free survival. Subclonal WGD was detected in 19% of tumours, and 10% of tumours harboured multiple subclonal WGDs in parallel. Subclonal, but not truncal, WGD was associated with shorter disease-free survival. Copy number heterogeneity was associated with extrathoracic relapse within 1 year after surgery. These data demonstrate the importance of clonal expansion, WGD and copy number instability in determining the timing and patterns of relapse in non-small cell lung cancer and provide a comprehensive clinical cancer evolutionary data resource.

This is a repository copy of Tracking early lung cancer metastatic dissemination in TRACERx using ctDNA.

Introduction Minimal residual disease (MRD) detection in solid tumors describes isolation of circulating tumor DNA (ctDNA) molecules in plasma following definitive treatment of a cancer. Detection of MRD following surgical tumor excision categorizes patients as high risk for disease recurrence. Establishing an MRD approach to treating early-stage NSCLC will facilitate escalation of standard of care (SoC) treatment only in patients destined to relapse from their cancer and overcome challenges associated with conventional adjuvant drug-trial design. Here, we present data from the lung TRACERx study where patients with early-stage NSCLC underwent phylogenetic ctDNA profiling following resection. Methods Patient specific anchored-multiplex PCR (AMP) enrichment panels were generated for 78 lung TRACERx patients who underwent surgery for stage I-III NSCLC; 608 plasma samples were analyzed. Extensive patient-specific cfDNA enrichment panels targeted a median of 196 (range 72 to 482) clonal and subclonal variants detected in primary tumor tissue by multi-region exome sequencing. A novel MRD-caller controlled and estimated background sequencing error to maximize ctDNA detection at low mutant allele frequencies (MAFs). Analytical validation experiments benchmarked assay performance. Results Analytical validation of a 50-variant AMP-MRD assay demonstrated a sensitivity of 89% for mutant DNA at a MAF of 0.008% (with 25ng of DNA input into the assay), specificity was 100% experimentally and 99.9% (95% CI: 99.67 to 99.99%) modelled in-silico. 45 patients suffered relapse of their primary NSCLC; ctDNA was detected at or before clinical relapse in 37 of 45 patients. In these 37 patients the median ctDNA lead-time (time from ctDNA detection to clinical relapse) was 151 days (range 0 to 984 days) and the median time to relapse from surgery was 413 days (range 41 to 1242 days). In 10 of 10 patients who developed second primary cancers during follow-up no ctDNA was detected, reflecting specificity of the MRD assay toward the primary tumor. In 23 patients who remained relapse-free during a median of 1184 days of study follow-up, ctDNA was detected in 1 of 199 time-points analyzed. Analysis of SoC adjuvant surveillance imaging (CT, PET-CT or MRI, 220 encounters) revealed examples of MRD positive patients where SoC radiological surveillance was negative for impending relapse. Through application of large cfDNA enrichment panels targeting up to 483 variants per patient we observed dynamic changes in clonal composition and copy-number status prior to NSCLC relapse, categorized relapse as monoclonal or polyclonal and identified distinct subclonal dynamics during systemic intervention for disease recurrence. Conclusions ctDNA is an adjuvant biomarker capable of both detecting MRD following surgery and defining the clonality of relapsing disease. These data pave the way for clinical trials predicated on escalation of adjuvant standard of care in NSCLC patients who exhibit MRD positive status following surgery. Citation Format: Chris Abbosh, Alexander Frankell, Aaron Garnett, Thomas Harrison, Morgan Weichert, Abel Licon, Selvaraju Veeriah, Bob Daber, Mike Moreau, Adrian Chesh, Kevin Litchfield, Emilia Lim, Daniel Cooke, Clare Puttick, Maise Al Bakir, Fabio Gomes, Akshay Patel, Lizi Manzano, Ariana Huebner, Nicolas Carey, Joan Riley, Paula Roberts, Todd Druley, Jacqui A. Shaw, Nicholas McGranahan, Mariam Jamal-Hanjani, Nicolai Birkbak, Josh Stahl, Charles Swanton, Lung TRACERx consortium. Phylogenetic tracking and minimal residual disease detection using ctDNA in early-stage NSCLC: A lung TRACERx study [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 CT023.

Murine tissues harbor signature γδ T cell compartments with profound yet differential impacts on carcinogenesis. Conversely, human tissue-resident γδ cells are less well defined. In the present study, we show that human lung tissues harbor a resident Vδ1 γδ T cell population. Moreover, we demonstrate that Vδ1 T cells with resident memory and effector memory phenotypes were enriched in lung tumors compared with nontumor lung tissues. Intratumoral Vδ1 T cells possessed stem-like features and were skewed toward cytolysis and helper T cell type 1 function, akin to intratumoral natural killer and CD8+ T cells considered beneficial to the patient. Indeed, ongoing remission post-surgery was significantly associated with the numbers of CD45RA−CD27− effector memory Vδ1 T cells in tumors and, most strikingly, with the numbers of CD103+ tissue-resident Vδ1 T cells in nonmalignant lung tissues. Our findings offer basic insights into human body surface immunology that collectively support integrating Vδ1 T cell biology into immunotherapeutic strategies for nonsmall cell lung cancer.