Evolution of Metastases in Space and Time under Immune Selection.

We examined how the immune microenvironment molds tumor evolution at different metastatic organs in a longitudinal dataset of colorectal cancer. Through multiplexed analyses, we showed that clonal evolution patterns during metastatic progression depend on the immune contexture at the metastatic site. Genetic evidence of neoantigen depletion was observed in the sites with high Immunoscore and spatial proximity between Ki67+ tumor cells and CD3+ cells. The immunoedited tumor clones were eliminated and did not recur, while progressing clones were immune privileged, despite the presence of tumor-infiltrating lymphocytes. Characterization of immune-privileged metastases revealed tumor-intrinsic and tumor-extrinsic mechanisms of escape. The lowest recurrence risk was associated with high Immunoscore, occurrence of immunoediting, and low tumor burden. We propose a parallel selection model of metastatic progression, where branched evolution could be traced back to immune-escaping clones. The findings could inform the understanding of cancer dissemination and the development of immunotherapeutics.

The evolution of lung cancer and impact of subclonal selection in TRACERx.

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.

Integrative Analyses of Colorectal Cancer Show Immunoscore Is a Stronger Predictor of Patient Survival Than Microsatellite Instability.

Microsatellite instability in colorectal cancer predicts favorable outcomes. However, the mechanistic relationship between microsatellite instability, tumor-infiltrating immune cells, Immunoscore, and their impact on patient survival remains to be elucidated. We found significant differences in mutational patterns, chromosomal instability, and gene expression that correlated with patient microsatellite instability status. A prominent immune gene expression was observed in microsatellite-instable (MSI) tumors, as well as in a subgroup of microsatellite-stable (MSS) tumors. MSI tumors had increased frameshift mutations, showed genetic evidence of immunoediting, had higher densities of Th1, effector-memory T cells, in situ proliferating T cells, and inhibitory PD1-PDL1 cells, had high Immunoscores, and were infiltrated with mutation-specific cytotoxic T cells. Multivariate analysis revealed that Immunoscore was superior to microsatellite instability in predicting patients' disease-specific recurrence and survival. These findings indicate that assessment of the immune status via Immunoscore provides a potent indicator of tumor recurrence beyond microsatellite-instability staging that could be an important guide for immunotherapy strategies.

Immune evasion before tumour invasion in early lung squamous carcinogenesis.

Early detection and treatment are critical for improving the outcome of patients with cancer1. Understanding the largely uncharted biology of carcinogenesis requires deciphering molecular processes in premalignant lesions, and revealing the determinants of the intralesional immune reaction during cancer development. The adaptive immune response within tumours has previously been shown to be strongest at the earliest stage of carcinoma2,3. Here we show that immune activation and immune escape occur before tumour invasion, and reveal the relevant immune biomarkers of the pre-invasive stages of carcinogenesis in the lung. We used gene-expression profiling and multispectral imaging to analyse a dataset of 9 morphological stages of the development of lung squamous cell carcinoma, which includes 122 well-annotated biopsies from 77 patients. We identified evolutionary trajectories of cancer and immune pathways that comprise (1) a linear increase in proliferation and DNA repair from normal to cancerous tissue; (2) a transitory increase of metabolism and early immune sensing, through the activation of resident immune cells, in low-grade pre-invasive lesions; (3) the activation of immune responses and immune escape through immune checkpoints and suppressive interleukins from high-grade pre-invasive lesions; and, ultimately, (4) the activation of the epithelial-mesenchymal transition in the invasive stage of cancer. We propose that carcinogenesis in the lung involves a dynamic co-evolution of pre-invasive bronchial cells and the immune response. These findings highlight the need to develop immune biomarkers for early detection as well as immunotherapy-based chemopreventive approaches for individuals who are at high risk of developing lung cancer.

Endosome and Golgi-associated degradation (EGAD) of membrane proteins regulates sphingolipid metabolism.

Cellular homeostasis requires the ubiquitin-dependent degradation of membrane proteins. This was assumed to be mediated exclusively either by endoplasmic reticulum-associated degradation (ERAD) or by endosomal sorting complexes required for transport (ESCRT)-dependent lysosomal degradation. We identified in Saccharomyces cerevisiae an additional pathway that selectively extracts membrane proteins at Golgi and endosomes for degradation by cytosolic proteasomes. One endogenous substrate of this endosome and Golgi-associated degradation pathway (EGAD) is the ER-resident membrane protein Orm2, a negative regulator of sphingolipid biosynthesis. Orm2 degradation is initiated by phosphorylation, which triggers its ER export. Once on Golgi and endosomes, Orm2 is poly-ubiquitinated by the membrane-embedded "Defective in SREBP cleavage" (Dsc) ubiquitin ligase complex. Cdc48/VCP then extracts ubiquitinated Orm2 from membranes, which is tightly coupled to the proteasomal degradation of Orm2. Thereby, EGAD prevents the accumulation of Orm2 at the ER and in post-ER compartments and promotes the controlled de-repression of sphingolipid biosynthesis. Thus, the selective degradation of membrane proteins by EGAD contributes to proteostasis and lipid homeostasis in eukaryotic cells.

TOR complex 2 (TORC2) signaling and the ESCRT machinery cooperate in the protection of plasma membrane integrity in yeast.

The endosomal sorting complexes required for transport (ESCRT) mediate evolutionarily conserved membrane remodeling processes. Here, we used budding yeast (Saccharomyces cerevisiae) to explore how the ESCRT machinery contributes to plasma membrane (PM) homeostasis. We found that in response to reduced membrane tension and inhibition of TOR complex 2 (TORC2), ESCRT-III/Vps4 assemblies form at the PM and help maintain membrane integrity. In turn, the growth of ESCRT mutants strongly depended on TORC2-mediated homeostatic regulation of sphingolipid (SL) metabolism. This was caused by calcineurin-dependent dephosphorylation of Orm2, a repressor of SL biosynthesis. Calcineurin activity impaired Orm2 export from the endoplasmic reticulum (ER) and thereby hampered its subsequent endosome and Golgi-associated degradation (EGAD). The ensuing accumulation of Orm2 at the ER in ESCRT mutants necessitated TORC2 signaling through its downstream kinase Ypk1, which repressed Orm2 and prevented a detrimental imbalance of SL metabolism. Our findings reveal compensatory cross-talk between the ESCRT machinery, calcineurin/TORC2 signaling, and the EGAD pathway important for the regulation of SL biosynthesis and the maintenance of PM homeostasis.

Deep cell phenotyping and spatial analysis of multiplexed imaging with TRACERx-PHLEX.

The growing scale and dimensionality of multiplexed imaging require reproducible and comprehensive yet user-friendly computational pipelines. TRACERx-PHLEX performs deep learning-based cell segmentation (deep-imcyto), automated cell-type annotation (TYPEx) and interpretable spatial analysis (Spatial-PHLEX) as three independent but interoperable modules. PHLEX generates single-cell identities, cell densities within tissue compartments, marker positivity calls and spatial metrics such as cellular barrier scores, along with summary graphs and spatial visualisations. PHLEX was developed using imaging mass cytometry (IMC) in the TRACERx study, validated using published Co-detection by indexing (CODEX), IMC and orthogonal data and benchmarked against state-of-the-art approaches. We evaluated its use on different tissue types, tissue fixation conditions, image sizes and antibody panels. As PHLEX is an automated and containerised Nextflow pipeline, manual assessment, programming skills or pathology expertise are not essential. PHLEX offers an end-to-end solution in a growing field of highly multiplexed data and provides clinically relevant insights.

Spatial Architecture of Myeloid and T Cells Orchestrates Immune Evasion and Clinical Outcome in Lung Cancer.

Understanding the role of the tumor microenvironment (TME) in lung cancer is critical to improving patient outcomes. We identified four histology-independent archetype TMEs in treatment-naïve early-stage lung cancer using imaging mass cytometry in the TRACERx study (n = 81 patients/198 samples/2.3 million cells). In immune-hot adenocarcinomas, spatial niches of T cells and macrophages increased with clonal neoantigen burden, whereas such an increase was observed for niches of plasma and B cells in immune-excluded squamous cell carcinomas (LUSC). Immune-low TMEs were associated with fibroblast barriers to immune infiltration. The fourth archetype, characterized by sparse lymphocytes and high tumor-associated neutrophil (TAN) infiltration, had tumor cells spatially separated from vasculature and exhibited low spatial intratumor heterogeneity. TAN-high LUSC had frequent PIK3CA mutations. TAN-high tumors harbored recently expanded and metastasis-seeding subclones and had a shorter disease-free survival independent of stage. These findings delineate genomic, immune, and physical barriers to immune surveillance and implicate neutrophil-rich TMEs in metastasis. SIGNIFICANCE: This study provides novel insights into the spatial organization of the lung cancer TME in the context of tumor immunogenicity, tumor heterogeneity, and cancer evolution. Pairing the tumor evolutionary history with the spatially resolved TME suggests mechanistic hypotheses for tumor progression and metastasis with implications for patient outcome and treatment. This article is featured in Selected Articles from This Issue, p. 897.

The role of APOBEC3B in lung tumor evolution and targeted cancer therapy resistance.

In this study, the impact of the apolipoprotein B mRNA-editing catalytic subunit-like (APOBEC) enzyme APOBEC3B (A3B) on epidermal growth factor receptor (EGFR)-driven lung cancer was assessed. A3B expression in EGFR mutant (EGFRmut) non-small-cell lung cancer (NSCLC) mouse models constrained tumorigenesis, while A3B expression in tumors treated with EGFR-targeted cancer therapy was associated with treatment resistance. Analyses of human NSCLC models treated with EGFR-targeted therapy showed upregulation of A3B and revealed therapy-induced activation of nuclear factor kappa B (NF-κB) as an inducer of A3B expression. Significantly reduced viability was observed with A3B deficiency, and A3B was required for the enrichment of APOBEC mutation signatures, in targeted therapy-treated human NSCLC preclinical models. Upregulation of A3B was confirmed in patients with NSCLC treated with EGFR-targeted therapy. This study uncovers the multifaceted roles of A3B in NSCLC and identifies A3B as a potential target for more durable responses to targeted cancer therapy.

Mixed responses to targeted therapy driven by chromosomal instability through p53 dysfunction and genome doubling.

The phenomenon of mixed/heterogenous treatment responses to cancer therapies within an individual patient presents a challenging clinical scenario. Furthermore, the molecular basis of mixed intra-patient tumor responses remains unclear. Here, we show that patients with metastatic lung adenocarcinoma harbouring co-mutations of EGFR and TP53, are more likely to have mixed intra-patient tumor responses to EGFR tyrosine kinase inhibition (TKI), compared to those with an EGFR mutation alone. The combined presence of whole genome doubling (WGD) and TP53 co-mutations leads to increased genome instability and genomic copy number aberrations in genes implicated in EGFR TKI resistance. Using mouse models and an in vitro isogenic p53-mutant model system, we provide evidence that WGD provides diverse routes to drug resistance by increasing the probability of acquiring copy-number gains or losses relative to non-WGD cells. These data provide a molecular basis for mixed tumor responses to targeted therapy, within an individual patient, with implications for therapeutic strategies.

APOBEC3 as a driver of genetic intratumor heterogeneity.

Our recent study revealed that APOBEC3B is upregulated during the preinvasive stages of non-small cell lung cancer and breast cancer. In addition to its role in mediating single nucleotide variants, we propose that APOBEC3 promotes copy number intratumor heterogeneity prior to invasion, providing a substrate for cancer evolution.

Quantifying the impact of immunotherapy on RNA dynamics in cancer.

BACKGROUND: Checkpoint inhibitor (CPI) immunotherapies have provided durable clinical responses across a range of solid tumor types for some patients with cancer. Nonetheless, response rates to CPI vary greatly between cancer types. Resolving intratumor transcriptomic changes induced by CPI may improve our understanding of the mechanisms of sensitivity and resistance. METHODS: We assembled a cohort of longitudinal pre-therapy and on-therapy samples from 174 patients treated with CPI across six cancer types by leveraging transcriptomic sequencing data from five studies. RESULTS: Meta-analyses of published RNA markers revealed an on-therapy pattern of immune reinvigoration in patients with breast cancer, which was not discernible pre-therapy, providing biological insight into the impact of CPI on the breast cancer immune microenvironment. We identified 98 breast cancer-specific correlates of CPI response, including 13 genes which are known IO targets, such as toll-like receptors TLR1, TLR4, and TLR8, that could hold potential as combination targets for patients with breast cancer receiving CPI treatment. Furthermore, we demonstrate that a subset of response genes identified in breast cancer are already highly expressed pre-therapy in melanoma, and additionally we establish divergent RNA dynamics between breast cancer and melanoma following CPI treatment, which may suggest distinct immune microenvironments between the two cancer types. CONCLUSIONS: Overall, delineating longitudinal RNA dynamics following CPI therapy sheds light on the mechanisms underlying diverging response trajectories, and identifies putative targets for combination therapy.

Bioinformatics for cancer immunology and immunotherapy.

Recent mechanistic insights obtained from preclinical studies and the approval of the first immunotherapies has motivated increasing number of academic investigators and pharmaceutical/biotech companies to further elucidate the role of immunity in tumor pathogenesis and to reconsider the role of immunotherapy. Additionally, technological advances (e.g., next-generation sequencing) are providing unprecedented opportunities to draw a comprehensive picture of the tumor genomics landscape and ultimately enable individualized treatment. However, the increasing complexity of the generated data and the plethora of bioinformatics methods and tools pose considerable challenges to both tumor immunologists and clinical oncologists. In this review, we describe current concepts and future challenges for the management and analysis of data for cancer immunology and immunotherapy. We first highlight publicly available databases with specific focus on cancer immunology including databases for somatic mutations and epitope databases. We then give an overview of the bioinformatics methods for the analysis of next-generation sequencing data (whole-genome and exome sequencing), epitope prediction tools as well as methods for integrative data analysis and network modeling. Mathematical models are powerful tools that can predict and explain important patterns in the genetic and clinical progression of cancer. Therefore, a survey of mathematical models for tumor evolution and tumor-immune cell interaction is included. Finally, we discuss future challenges for individualized immunotherapy and suggest how a combined computational/experimental approaches can lead to new insights into the molecular mechanisms of cancer, improved diagnosis, and prognosis of the disease and pinpoint novel therapeutic targets.

A local human Vδ1 T cell population is associated with survival in nonsmall-cell lung cancer.

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.

Clonal somatic copy number altered driver events inform drug sensitivity in high-grade serous ovarian cancer.

Chromosomal instability is a major challenge to patient stratification and targeted drug development for high-grade serous ovarian carcinoma (HGSOC). Here we show that somatic copy number alterations (SCNAs) in frequently amplified HGSOC cancer genes significantly correlate with gene expression and methylation status. We identify five prevalent clonal driver SCNAs (chromosomal amplifications encompassing MYC, PIK3CA, CCNE1, KRAS and TERT) from multi-regional HGSOC data and reason that their strong selection should prioritise them as key biomarkers for targeted therapies. We use primary HGSOC spheroid models to test interactions between in vitro targeted therapy and SCNAs. MYC chromosomal copy number is associated with in-vitro and clinical response to paclitaxel and in-vitro response to mTORC1/2 inhibition. Activation of the mTOR survival pathway in the context of MYC-amplified HGSOC is statistically associated with increased prevalence of SCNAs in genes from the PI3K pathway. Co-occurrence of amplifications in MYC and genes from the PI3K pathway is independently observed in squamous lung cancer and triple negative breast cancer. In this work, we show that identifying co-occurrence of clonal driver SCNA genes could be used to tailor therapeutics for precision medicine.

Evasion before invasion: Pre-cancer immunosurveillance.

Early intervention and risk stratification solutions for lung cancer are limited by our understanding of how carcinogenesis transforms the pre-invasive epithelium and its microenvironment before the carcinoma stage. We describe the sequence of molecular and cellular changes leading to cancer formation and the co-evolution of the earliest immune response. We revealed that immune sensing, infiltration and activation of immune cells, immune escape, and microenvironment reorganization occur early in pre-cancer. These findings urge the need for broadening the scope of the established immunotherapy approaches toward prophylactic cancer treatment and preventive intervention. Leveraging the immune contexture and the mechanisms of immune modulation for individuals at risk of developing cancer and further to the general population will allow for early detection, chemoprevention, and risk stratification in the near future.

Characterisation of tumour microenvironment remodelling following oncogene inhibition in preclinical studies with imaging mass cytometry.

Mouse models are critical in pre-clinical studies of cancer therapy, allowing dissection of mechanisms through chemical and genetic manipulations that are not feasible in the clinical setting. In studies of the tumour microenvironment (TME), multiplexed imaging methods can provide a rich source of information. However, the application of such technologies in mouse tissues is still in its infancy. Here we present a workflow for studying the TME using imaging mass cytometry with a panel of 27 antibodies on frozen mouse tissues. We optimise and validate image segmentation strategies and automate the process in a Nextflow-based pipeline (imcyto) that is scalable and portable, allowing for parallelised segmentation of large multi-image datasets. With these methods we interrogate the remodelling of the TME induced by a KRAS G12C inhibitor in an immune competent mouse orthotopic lung cancer model, highlighting the infiltration and activation of antigen presenting cells and effector cells.

Induction of APOBEC3 Exacerbates DNA Replication Stress and Chromosomal Instability in Early Breast and Lung Cancer Evolution.

APOBEC3 enzymes are cytosine deaminases implicated in cancer. Precisely when APOBEC3 expression is induced during cancer development remains to be defined. Here we show that specific APOBEC3 genes are upregulated in breast ductal carcinoma in situ, and in preinvasive lung cancer lesions coincident with cellular proliferation. We observe evidence of APOBEC3-mediated subclonal mutagenesis propagated from TRACERx preinvasive to invasive non-small cell lung cancer (NSCLC) lesions. We find that APOBEC3B exacerbates DNA replication stress and chromosomal instability through incomplete replication of genomic DNA, manifested by accumulation of mitotic ultrafine bridges and 53BP1 nuclear bodies in the G1 phase of the cell cycle. Analysis of TRACERx NSCLC clinical samples and mouse lung cancer models revealed APOBEC3B expression driving replication stress and chromosome missegregation. We propose that APOBEC3 is functionally implicated in the onset of chromosomal instability and somatic mutational heterogeneity in preinvasive disease, providing fuel for selection early in cancer evolution. SIGNIFICANCE: This study reveals the dynamics and drivers of APOBEC3 gene expression in preinvasive disease and the exacerbation of cellular diversity by APOBEC3B through DNA replication stress to promote chromosomal instability early in cancer evolution.This article is highlighted in the In This Issue feature, p. 2355.

LAMTOR/Ragulator regulates lipid metabolism in macrophages and foam cell differentiation.

Late endosomal/lysosomal adaptor and MAPK and mTOR activator (LAMTOR/Ragulator) is a scaffold protein complex that anchors and regulates multiprotein signaling units on late endosomes/lysosomes. To identify LAMTOR-modulated endolysosomal proteins, primary macrophages were derived from bone marrow of conditional knockout mice carrying a specific deletion of LAMTOR2 in the monocyte/macrophage cell lineage. Affymetrix-based transcriptomic analysis and quantitative iTRAQ-based organelle proteomic analysis of endosomes derived from macrophages were performed. Further analyses showed that LAMTOR could be a novel regulator of foam cell differentiation. The lipid droplet formation phenotype observed in macrophages was additionally confirmed in MEFs, where lipidomic analysis identified cholesterol esters as specifically downregulated in LAMTOR2 knockout cells. The data obtained indicate a function of LAMTOR2 in lipid metabolism.