Search

Cancer of unknown primary (CUP) tumors are biologically very heterogeneous, which complicates stratification of patients for treatment. Consequently, these patients face limited treatment options and a poor prognosis. With this study, we aim to expand on the current knowledge of CUP biology by analyzing two cohorts: a well-characterized cohort of 44 CUP patients, and 213 metastatic patients with known primary. These cohorts were treated at the same institution and characterized by identical molecular assessments. Through comparative analysis of genomic and transcriptomic data, we found that CUP tumors were characterized by high expression of immune-related genes and pathways compared to other metastatic tumors. Moreover, CUP tumors uniformly demonstrated high levels of tumor-infiltrating leukocytes and circulating T cells, indicating a strong immune response. Finally, the genetic landscape of CUP tumors resembled that of other metastatic cancers and demonstrated mutations in established cancer genes. In conclusion, CUP tumors possess a distinct immunophenotype that distinguishes them from other metastatic cancers. These results may suggest an immune response in CUP that facilitates metastatic tumor growth while limiting growth of the primary tumor., (© 2024 The Authors. Molecular Oncology published by John Wiley & Sons Ltd on behalf of Federation of European Biochemical Societies.)

Background: Circulating tumor DNA (ctDNA) has emerged as a promising tool for early cancer detection and minimal residual disease monitoring. However, the biology underlying ctDNA release and its variation across cancer types and histologies remains poorly understood. This study investigated the biology behind ctDNA shedding in colorectal cancer., Methods: The study included a local cohort of 747 stage I-III colorectal cancer patients. All patients had ctDNA measurement prior to treatment and extensive clinical data. Primary tumor RNA sequencing and whole exome sequencing was performed in 95 and 652 patients respectively. Additionally, the study evaluated 89 non-small cell lung cancer patients from the TRACERx cohort, comprising primary tumor RNA sequencing and ctDNA measurement., Results: We found tumor size and proliferative capacity to be key factors associated with ctDNA shedding in colorectal cancer. Furthermore, we found that the secretory and CMS3 colorectal cancer subtypes exhibited lower ctDNA shedding, while microsatellite instability (MSI) tumors had higher levels of ctDNA. Mutational analysis did not reveal any genes or pathways associated with ctDNA shedding in colorectal cancer. A comparison of transcriptomic profiles across multiple cancer types demonstrated that colorectal cancer and lung squamous cell carcinoma tumors shared a high-proliferative ctDNA shedding phenotype, while lung adenocarcinoma tumors displayed a distinct low-proliferative subgroup. Additionally, proliferation levels correlated with ctDNA detection sensitivity across multiple cancer types., Conclusion: These findings suggest that tumor size and proliferative capacity are drivers of ctDNA release in colorectal cancer and provide insights into the biology of ctDNA shedding on a pan-cancer level., Competing Interests: Declaration of Competing Interest The authors declare that they have no known competing financial interests or personal relationships that could have appeared to influence the work reported in this paper., (Copyright © 2024 The Authors. Published by Elsevier Ltd.. All rights reserved.)

Purpose: Cancer continues to be a major global health challenge, affecting millions of individuals and placing substantial burdens on healthcare systems worldwide. Recent research suggests a complex relationship between obesity and cancer, with obesity increasing the risk of various cancers while potentially improving outcomes for diagnosed patients, a phenomenon termed the "obesity paradox". In this study, we used a cohort of 1781 patients to investigate the impact of obesity on tumor characteristics, including gene expression, pathway dysfunction, genetic alterations and immune infiltration., Methods: Patient samples spanned 10 different cancer types, and were obtained from the Cancer Genome Atlas, with annotations for body mass index (BMI), age, sex, tumor size and tumor gene expression data., Results: When we compared the proportion of large (T3-T4) to small tumors (T1-T2) between obese and non-obese patients, we found that obese patients tended to present with smaller, less invasive tumors and exhibited distinct gene expression profiles, particularly in metabolic and proliferative pathways. Moreover, smaller tumors in obese patients show higher immune cell infiltration and increased T cell diversity, suggesting enhanced immune activity., Conclusion: Taken together, these findings highlight the influence of obesity on tumor biology, with implications for personalized treatment strategies that consider patient physiology alongside tumor characteristics., (© 2024. The Author(s).)

Not available.

Intratumor heterogeneity provides the fuel for the evolution and selection of subclonal tumor cell populations. However, accurate inference of tumor subclonal architecture and reconstruction of tumor evolutionary histories from bulk DNA sequencing data remains challenging. Frequently, sequencing and alignment artifacts are not fully filtered out from cancer somatic mutations, and errors in the identification of copy number alterations or complex evolutionary events (e.g., mutation losses) affect the estimated cellular prevalence of mutations. Together, such errors propagate into the analysis of mutation clustering and phylogenetic reconstruction. In this Protocol, we present a new computational framework, CONIPHER (COrrecting Noise In PHylogenetic Evaluation and Reconstruction), that accurately infers subclonal structure and phylogenetic relationships from multisample tumor sequencing, accounting for both copy number alterations and mutation errors. CONIPHER has been used to reconstruct subclonal architecture and tumor phylogeny from multisample tumors with high-depth whole-exome sequencing from the TRACERx421 dataset, as well as matched primary-metastatic cases. CONIPHER outperforms similar methods on simulated datasets, and in particular scales to a large number of tumor samples and clones, while completing in under 1.5 h on average. CONIPHER enables automated phylogenetic analysis that can be effectively applied to large sequencing datasets generated with different technologies. CONIPHER can be run with a basic knowledge of bioinformatics and R and bash scripting languages., (© 2023. Springer Nature Limited.)

Foundation models in deep learning are characterized by a single large-scale model trained on vast amounts of data serving as the foundation for various downstream tasks. Foundation models are generally trained using self-supervised learning and excel in reducing the demand for training samples in downstream applications. This is especially important in medicine, where large labelled datasets are often scarce. Here, we developed a foundation model for cancer imaging biomarker discovery by training a convolutional encoder through self-supervised learning using a comprehensive dataset of 11,467 radiographic lesions. The foundation model was evaluated in distinct and clinically relevant applications of cancer imaging-based biomarkers. We found that it facilitated better and more efficient learning of imaging biomarkers and yielded task-specific models that significantly outperformed conventional supervised and other state-of-the-art pretrained implementations on downstream tasks, especially when training dataset sizes were very limited. Furthermore, the foundation model was more stable to input variations and showed strong associations with underlying biology. Our results demonstrate the tremendous potential of foundation models in discovering new imaging biomarkers that may extend to other clinical use cases and can accelerate the widespread translation of imaging biomarkers into clinical settings., Competing Interests: Competing interestsThe authors declare no competing interests., (© The Author(s) 2024.)

Immunotherapy has revolutionized treatment of patients diagnosed with metastatic melanoma, where nearly half of patients receive clinical benefit. However, immunotherapy is also associated with immune-related adverse events, which may be severe and persistent. It is therefore important to identify patients that do not benefit from therapy early. Currently, regularly scheduled CT scans are used to investigate size changes in target lesions to evaluate progression and therapy response. This study aims to explore if panel-based analysis of circulating tumor DNA (ctDNA) taken at 3-week intervals may provide a window into the growing cancer, can be used to identify nonresponding patients early, and determine genomic alterations associated with acquired resistance to checkpoint immunotherapy without analysis of tumor tissue biopsies. We designed a gene panel for ctDNA analysis and sequenced 4-6 serial plasma samples from 24 patients with unresectable stage III or IV melanoma and treated with first-line checkpoint inhibitors enrolled at the Department of Oncology at Aarhus University Hospital, Denmark. TERT was the most mutated gene found in ctDNA and associated with a poor prognosis. We detected more ctDNA in patients with high metastatic load, which indicates that more aggressive tumors release more ctDNA into the bloodstream. Although we did not find evidence of specific mutations associated with acquired resistance, we did demonstrate in this limited cohort of 24 patients that untargeted, panel-based ctDNA analysis has the potential to be used as a minimally invasive tool in clinical practice to identify patients where the benefits of immunotherapy outweigh the drawbacks., (Copyright © 2023 The Author(s). Published by Wolters Kluwer Health, Inc.)

The application of next-generation sequencing (NGS) has transformed cancer research. As costs have decreased, NGS has increasingly been applied to generate multiple layers of molecular data from the same samples, covering genomics, transcriptomics, and methylomics. Integrating these types of multi-omics data in a combined analysis is now becoming a common issue with no obvious solution, often handled on an ad hoc basis, with multi-omics data arriving in a tabular format and analyzed using computationally intensive statistical methods. These methods particularly ignore the spatial orientation of the genome and often apply stringent p-value corrections that likely result in the loss of true positive associations. Here, we present GENIUS (GEnome traNsformatIon and spatial representation of mUltiomicS data), a framework for integrating multi-omics data using deep learning models developed for advanced image analysis. The GENIUS framework is able to transform multi-omics data into images with genes displayed as spatially connected pixels and successfully extract relevant information with respect to the desired output. We demonstrate the utility of GENIUS by applying the framework to multi-omics datasets from the Cancer Genome Atlas. Our results are focused on predicting the development of metastatic cancer from primary tumors, and demonstrate how through model inference, we are able to extract the genes which are driving the model prediction and are likely associated with metastatic disease progression. We anticipate our framework to be a starting point and strong proof of concept for multi-omics data transformation and analysis without the need for statistical correction., Competing Interests: MS, AK, LD, BH, HJ, NB No competing interests declared, (© 2023, Sokač et al.)

Metastatic disease is responsible for the majority of cancer-related deaths 1 . We report the longitudinal evolutionary analysis of 126 non-small cell lung cancer (NSCLC) tumours from 421 prospectively recruited patients in TRACERx who developed metastatic disease, compared with a control cohort of 144 non-metastatic tumours. In 25% of cases, metastases diverged early, before the last clonal sweep in the primary tumour, and early divergence was enriched for patients who were smokers at the time of initial diagnosis. Simulations suggested that early metastatic divergence more frequently occurred at smaller tumour diameters (less than 8 mm). Single-region primary tumour sampling resulted in 83% of late divergence cases being misclassified as early, highlighting the importance of extensive primary tumour sampling. Polyclonal dissemination, which was associated with extrathoracic disease recurrence, was found in 32% of cases. Primary lymph node disease contributed to metastatic relapse in less than 20% of cases, representing a hallmark of metastatic potential rather than a route to subsequent recurrences/disease progression. Metastasis-seeding subclones exhibited subclonal expansions within primary tumours, probably reflecting positive selection. Our findings highlight the importance of selection in metastatic clone evolution within untreated primary tumours, the distinction between monoclonal versus polyclonal seeding in dictating site of recurrence, the limitations of current radiological screening approaches for early diverging tumours and the need to develop strategies to target metastasis-seeding subclones before relapse., (© 2023. The Author(s).)

Lung cancer is the leading cause of cancer-associated mortality worldwide 1 . 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., (© 2023. The Author(s).)

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 elusive 1,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 adenocarcinoma 3 . 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., (© 2023. The Author(s).)

Cancer-associated cachexia (CAC) is a major contributor to morbidity and mortality in individuals with non-small cell lung cancer. Key features of CAC include alterations in body composition and body weight. Here, we explore the association between body composition and body weight with survival and delineate potential biological processes and mediators that contribute to the development of CAC. Computed tomography-based body composition analysis of 651 individuals in the TRACERx (TRAcking non-small cell lung Cancer Evolution through therapy (Rx)) study suggested that individuals in the bottom 20th percentile of the distribution of skeletal muscle or adipose tissue area at the time of lung cancer diagnosis, had significantly shorter lung cancer-specific survival and overall survival. This finding was validated in 420 individuals in the independent Boston Lung Cancer Study. Individuals classified as having developed CAC according to one or more features at relapse encompassing loss of adipose or muscle tissue, or body mass index-adjusted weight loss were found to have distinct tumor genomic and transcriptomic profiles compared with individuals who did not develop such features. Primary non-small cell lung cancers from individuals who developed CAC were characterized by enrichment of inflammatory signaling and epithelial-mesenchymal transitional pathways, and differentially expressed genes upregulated in these tumors included cancer-testis antigen MAGEA6 and matrix metalloproteinases, such as ADAMTS3. In an exploratory proteomic analysis of circulating putative mediators of cachexia performed in a subset of 110 individuals from TRACERx, a significant association between circulating GDF15 and loss of body weight, skeletal muscle and adipose tissue was identified at relapse, supporting the potential therapeutic relevance of targeting GDF15 in the management of CAC., (© 2023. The Author(s) under exclusive license to Springer Nature America, Inc.)

Circulating tumour DNA (ctDNA) can be used to detect and profile residual tumour cells persisting after curative intent therapy 1 . The study of large patient cohorts incorporating longitudinal plasma sampling and extended follow-up is required to determine the role of ctDNA as a phylogenetic biomarker of relapse in early-stage non-small-cell lung cancer (NSCLC). Here we developed ctDNA methods tracking a median of 200 mutations identified in resected NSCLC tissue across 1,069 plasma samples collected from 197 patients enrolled in the TRACERx study 2 . A lack of preoperative ctDNA detection distinguished biologically indolent lung adenocarcinoma with good clinical outcome. Postoperative plasma analyses were interpreted within the context of standard-of-care radiological surveillance and administration of cytotoxic adjuvant therapy. Landmark analyses of plasma samples collected within 120 days after surgery revealed ctDNA detection in 25% of patients, including 49% of all patients who experienced clinical relapse; 3 to 6 monthly ctDNA surveillance identified impending disease relapse in an additional 20% of landmark-negative patients. We developed a bioinformatic tool (ECLIPSE) for non-invasive tracking of subclonal architecture at low ctDNA levels. ECLIPSE identified patients with polyclonal metastatic dissemination, which was associated with a poor clinical outcome. By measuring subclone cancer cell fractions in preoperative plasma, we found that subclones seeding future metastases were significantly more expanded compared with non-metastatic subclones. Our findings will support (neo)adjuvant trial advances and provide insights into the process of metastatic dissemination using low-ctDNA-level liquid biopsy., (© 2023. The Author(s).)

Intratumour heterogeneity (ITH) fuels lung cancer evolution, which leads to immune evasion and resistance to therapy 1 . Here, using paired whole-exome and RNA sequencing data, we investigate intratumour transcriptomic diversity in 354 non-small cell lung cancer tumours from 347 out of the first 421 patients prospectively recruited into the TRACERx study 2,3 . Analyses of 947 tumour regions, representing both primary and metastatic disease, alongside 96 tumour-adjacent normal tissue samples implicate the transcriptome as a major source of phenotypic variation. Gene expression levels and ITH relate to patterns of positive and negative selection during tumour evolution. We observe frequent copy number-independent allele-specific expression that is linked to epigenomic dysfunction. Allele-specific expression can also result in genomic-transcriptomic parallel evolution, which converges on cancer gene disruption. We extract signatures of RNA single-base substitutions and link their aetiology to the activity of the RNA-editing enzymes ADAR and APOBEC3A, thereby revealing otherwise undetected ongoing APOBEC activity in tumours. Characterizing the transcriptomes of primary-metastatic tumour pairs, we combine multiple machine-learning approaches that leverage genomic and transcriptomic variables to link metastasis-seeding potential to the evolutionary context of mutations and increased proliferation within primary tumour regions. These results highlight the interplay between the genome and transcriptome in influencing ITH, lung cancer evolution and metastasis., (© 2023. The Author(s).)

Lung adenocarcinomas (LUADs) display a broad histological spectrum from low-grade lepidic tumors through to mid-grade acinar and papillary and high-grade solid, cribriform and micropapillary tumors. How morphology reflects tumor evolution and disease progression is poorly understood. Whole-exome sequencing data generated from 805 primary tumor regions and 121 paired metastatic samples across 248 LUADs from the TRACERx 421 cohort, together with RNA-sequencing data from 463 primary tumor regions, were integrated with detailed whole-tumor and regional histopathological analysis. Tumors with predominantly high-grade patterns showed increased chromosomal complexity, with higher burden of loss of heterozygosity and subclonal somatic copy number alterations. Individual regions in predominantly high-grade pattern tumors exhibited higher proliferation and lower clonal diversity, potentially reflecting large recent subclonal expansions. Co-occurrence of truncal loss of chromosomes 3p and 3q was enriched in predominantly low-/mid-grade tumors, while purely undifferentiated solid-pattern tumors had a higher frequency of truncal arm or focal 3q gains and SMARCA4 gene alterations compared with mixed-pattern tumors with a solid component, suggesting distinct evolutionary trajectories. Clonal evolution analysis revealed that tumors tend to evolve toward higher-grade patterns. The presence of micropapillary pattern and 'tumor spread through air spaces' were associated with intrathoracic recurrence, in contrast to the presence of solid/cribriform patterns, necrosis and preoperative circulating tumor DNA detection, which were associated with extra-thoracic recurrence. These data provide insights into the relationship between LUAD morphology, the underlying evolutionary genomic landscape, and clinical and anatomical relapse risk., (© 2023. The Author(s) under exclusive license to Springer Nature America, Inc.)

Objective: Circulating tumor DNA (ctDNA) is a candidate biomarker of cancer with practice-changing potential in the detection of both early and residual disease. Disease stage and tumor size affect the probability of ctDNA detection, whereas little is known about the influence of other tumor characteristics on ctDNA detection. This study investigates the impact of tumor cell whole-genome doubling (WGD) on the detection of ctDNA in plasma collected preoperatively from newly diagnosed colorectal cancer (CRC) patients., Methods: WGD was estimated from copy numbers derived from whole-exome sequencing (WES) data of matched tumor and normal DNA from 833 Danish CRC patients. To explore if tumor WGD status impacts ctDNA detection, we applied tumor-informed ctDNA analysis to preoperative plasma samples from all patients., Results: Patients with WGD+ tumors had 53% increased odds of being ctDNA positive (OR = 1.53, 95%CI: 1.12-2.09). After stratification for UICC stage, the association persisted for Stage I (OR = 2.44, 95%CI: 1.22-5.03) and Stage II (OR = 1.76, 95%CI: 1.11-2.81) but not for Stage III (OR = 0.83, 95%CI: 0.44-1.53) patients., Conclusion: The presence of WGD significantly increases the probability of detecting ctDNA, particularly for early-stage disease. In patients with more advanced disease, the benefit of WGD on ctDNA detection is less pronounced, consistent with increased DNA shedding from these tumors, making ctDNA detection less dependent on the amount of ctDNA released per tumor cell.

Immunotherapy has revolutionised cancer treatment. However, not all cancer patients benefit, and current stratification strategies based primarily on PD1 status and mutation burden are far from perfect. We hypothesised that high activation of an innate response relative to the adaptive response may prevent proper tumour neoantigen identification and decrease the specific anticancer response, both in the presence and absence of immunotherapy. To investigate this, we obtained transcriptomic data from three large publicly available cancer datasets, the Cancer Genome Atlas (TCGA), the Hartwig Medical Foundation (HMF), and a recently published cohort of metastatic bladder cancer patients treated with immunotherapy. To analyse immune infiltration into bulk tumours, we developed an RNAseq-based model based on previously published definitions to estimate the overall level of infiltrating innate and adaptive immune cells from bulk tumour RNAseq data. From these, the adaptive-to-innate immune ratio (A/I ratio) was defined. A meta-analysis of 32 cancer types from TCGA overall showed improved overall survival in patients with an A/I ratio above median (Hazard ratio (HR) females 0.73, HR males 0.86, P < 0.05). Of particular interest, we found that the association was different for males and females for eight cancer types, demonstrating a gender bias in the relative balance of the infiltration of innate and adaptive immune cells. For patients with metastatic disease, we found that responders to immunotherapy had a significantly higher A/I ratio than non-responders in HMF (P = 0.036) and a significantly higher ratio in complete responders in a separate metastatic bladder cancer dataset (P = 0.022). Overall, the adaptive-to-innate immune ratio seems to define separate states of immune activation, likely linked to fundamental immunological reactions to cancer. This ratio was associated with improved prognosis and improved response to immunotherapy, demonstrating potential relevance to patient stratification. Furthermore, by demonstrating a significant difference between males and females that associates with response, we highlight an important gender bias which likely has direct clinical relevance., Competing Interests: The authors have declared that no competing interests exist., (Copyright: © 2023 Ahrenfeldt et al. This is an open access article distributed under the terms of the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original author and source are credited.)

Cancer metastasis is the lethal developmental step in cancer, responsible for the majority of cancer deaths. To metastasise, cancer cells must acquire the ability to disseminate systemically and to escape an activated immune response. Here, we endeavoured to investigate if metastatic dissemination reflects acquisition of genomic traits that are selected for. We acquired mutation and copy number data from 8332 tumours representing 19 cancer types acquired from The Cancer Genome Atlas and the Hartwig Medical Foundation. A total of 827,344 non-synonymous mutations across 8332 tumour samples representing 19 cancer types were timed as early or late relative to copy number alterations, and potential driver events were annotated. We found that metastatic cancers had a significantly higher proportion of clonal mutations and a general enrichment of early mutations in p53 and RTK/KRAS pathways. However, while individual pathways demonstrated a clear time-separated preference for specific events, the relative timing did not vary between primary and metastatic cancers. These results indicate that the selective pressure that drives cancer development does not change dramatically between primary and metastatic cancer on a genomic level, and is mainly focused on alterations that increase proliferation.

Metastasis is the main cause of cancer death, yet the evolutionary processes behind it remain largely unknown. Here, through analysis of large panel-based genomic datasets from the AACR Genomics Evidence Neoplasia Information Exchange project, including 40,979 primary and metastatic tumors across 25 distinct cancer types, we explore how the evolutionary pressure of cancer metastasis shapes the selection of genomic drivers of cancer. The most commonly affected genes were TP53, MYC, and CDKN2A, with no specific pattern associated with metastatic disease. This suggests that, on a driver mutation level, the selective pressure operating in primary and metastatic tumors is similar. The most highly enriched individual driver mutations in metastatic tumors were mutations known to drive resistance to hormone therapies in breast and prostate cancer (ESR1 and AR), anti-EGFR therapy in non-small cell lung cancer (EGFR T790M), and imatinib in gastrointestinal cancer (KIT V654A). Specific mutational signatures were also associated with treatment in three cancer types, supporting clonal selection following anticancer therapy. Overall, this implies that initial acquisition of driver mutations is predominantly shaped by the tissue of origin, where specific mutations define the developing primary tumor and drive growth, immune escape, and tolerance to chromosomal instability. However, acquisition of driver mutations that contribute to metastatic disease is less specific, with the main genomic drivers of metastatic cancer evolution associating with resistance to therapy., Significance: This study leverages large datasets to investigate the evolutionary landscape of established cancer genes to shed new light upon the mystery of cancer dissemination and expand the understanding of metastatic cancer biology., (©2022 American Association for Cancer Research.)

The cGAS-STING pathway serves a critical role in anticancer therapy. Particularly, response to immunotherapy is likely driven by both active cGAS-STING signaling that attracts immune cells, and by the presence of cancer neoantigens that presents as targets for cytotoxic T cells. Chromosomal instability (CIN) is a hallmark of cancer, but also leads to an accumulation of cytosolic DNA that in turn results in increased cGAS-STING signaling. To avoid triggering the cGAS-STING pathway, it is commonly disrupted by cancer cells, either through mutations in the pathway or through transcriptional silencing. Given its effect on the immune system, determining the cGAS-STING activation status prior to treatment initiation is likely of clinical relevance. Here, we used combined expression data from 2,307 tumors from five cancer types from The Cancer Genome Atlas to define a novel cGAS-STING activity score based on eight genes with a known role in the pathway. Using unsupervised clustering, four distinct categories of cGAS-STING activation were identified. In multivariate models, the cGAS-STING active tumors show improved prognosis. Importantly, in an independent bladder cancer immunotherapy-treated cohort, patients with low cGAS-STING expression showed limited response to treatment, while patients with high expression showed improved response and prognosis, particularly among patients with high CIN and more neoantigens. In a multivariate model, a significant interaction was observed between CIN, neoantigens, and cGAS-STING activation. Together, this suggests a potential role of cGAS-STING activity as a predictive biomarker for the application of immunotherapy., Significance: The cGAS-STING pathway is induced by CIN, triggers inflammation and is often deficient in cancer. We provide a tool to evaluate cGAS-STING activity and demonstrate clinical significance in immunotherapy response., Competing Interests: M. Sokač reports grants from Lundbeck Foundation and Aarhus University Research Foundation during the conduct of the study. K. Litchfield reports personal fees from Monopteros therapeutics; grants from Ono pharmaceuticals and Genesis therapeutics outside the submitted work. L. Dyrskjøt reports grants from Ferring, Natera, C2i Genomics, AstraZeneca, Photocure, and personal fees from Ferring outside the submitted work. M.R. Jakobsen reports other from Stipe Therapeutics outside the submitted work. N.J. Birkbak reports grants from The Lundbeck Foundation, Aarhus University Research Foundation, and Danish Cancer Society during the conduct of the study; in addition, N.J. Birkbak has a patent to quantifying homologous recombination deficiency issued, licensed, and with royalties paid and a patent to a prognostic gene expression signature for lung cancer pending. No disclosures were reported by the other authors., (© 2022 The Authors; Published by the American Association for Cancer Research.)

Background: Checkpoint inhibitors have revolutionized the treatment of metastatic melanoma, yielding long-term survival in a considerable proportion of the patients. Yet, 40-60% of patients do not achieve a long-term benefit from such therapy, emphasizing the urgent need to identify biomarkers that can predict response to immunotherapy and guide patients for the best possible treatment. Here, we exploited an unsupervised machine learning approach to identify potential inflammatory cytokine signatures from liquid biopsies, which could predict response to immunotherapy in melanoma., Methods: We studied a cohort of 77 patients diagnosed with unresectable advanced-stage melanoma undergoing treatment with first-line nivolumab plus ipilimumab or pembrolizumab. Baseline and on-treatment plasma samples were tested for levels of PD-1, PD-L1, IFNγ, IFNβ, CCL20, CXCL5, CXCL10, IL6, IL8, IL10, MCP1, and TNFα and analyzed by Uniform Manifold Approximation and Projection (UMAP) dimension reduction method and k-means clustering analysis., Results: Interestingly, using UMAP analysis, we found that treatment-induced cytokine changes measured as a ratio between baseline and on-treatment samples correlated significantly to progression-free survival (PFS). For patients treated with nivolumab plus ipilimumab we identified a group of patients with superior PFS that were characterized by significantly higher baseline-to-on-treatment increments of PD-1, PD-L1, IFNγ, IL10, CXCL10, and TNFα compared to patients with worse PFS. Particularly, a high PD-1 increment was a strong individual predictor for superior PFS (HR = 0.13; 95% CI 0.034-0.49; p = 0.0026). In contrast, decreasing levels of IFNγ and IL6 and increasing levels of CXCL5 were associated with superior PFS in the pembrolizumab group, although none of the cytokines were individually predictors for PFS., Conclusions: In short, our study demonstrates that a high increment of PD-1 is associated with superior PFS in advanced-stage melanoma patients treated with nivolumab plus ipilimumab. In contrast, decreasing levels of IFNγ and IL6, and increasing levels of CXCL5 are associated with response to pembrolizumab. These results suggest that using serial samples to monitor changes in cytokine levels early during treatment is informative for treatment response.

The molecular landscape in non-muscle-invasive bladder cancer (NMIBC) is characterized by large biological heterogeneity with variable clinical outcomes. Here, we perform an integrative multi-omics analysis of patients diagnosed with NMIBC (n = 834). Transcriptomic analysis identifies four classes (1, 2a, 2b and 3) reflecting tumor biology and disease aggressiveness. Both transcriptome-based subtyping and the level of chromosomal instability provide independent prognostic value beyond established prognostic clinicopathological parameters. High chromosomal instability, p53-pathway disruption and APOBEC-related mutations are significantly associated with transcriptomic class 2a and poor outcome. RNA-derived immune cell infiltration is associated with chromosomally unstable tumors and enriched in class 2b. Spatial proteomics analysis confirms the higher infiltration of class 2b tumors and demonstrates an association between higher immune cell infiltration and lower recurrence rates. Finally, the independent prognostic value of the transcriptomic classes is documented in 1228 validation samples using a single sample classification tool. The classifier provides a framework for biomarker discovery and for optimizing treatment and surveillance in next-generation clinical trials.

Background: Intra-tumor heterogeneity (ITH) of colorectal cancer (CRC) complicates molecular tumor classification, such as transcriptional subtyping. Differences in cellular states, biopsy cell composition, and tumor microenvironment may all lead to ITH. Here we analyze ITH at the transcriptomic and proteomic levels to ascertain whether subtype discordance between multiregional biopsies reflects relevant biological ITH or lack of classifier robustness. Further, we study the impact of tumor location on ITH., Methods: Multiregional biopsies from stage II and III CRC tumors were analyzed by RNA sequencing (41 biopsies, 14 tumors) and multiplex immune protein analysis (89 biopsies, 29 tumors). CRC subtyping was performed using consensus molecular subtypes (CMS), CRC intrinsic subtypes (CRIS), and TUMOR types. ITH-scores and network maps were defined to determine the origin of heterogeneity. A validation cohort was used with one biopsy per tumor (162 tumors)., Results: Overall, inter-tumor transcriptional variation exceeded ITH, and subtyping calls were frequently concordant between multiregional biopsies. Still, some tumors had high transcriptional ITH and were classified discordantly. Subtyping of proximal MSS tumors were discordant for 50% of the tumors, this ITH was related to differences in the microenvironment. Subtyping of distal MSS tumors were less discordant, here the ITH was more cancer-cell related. The subtype discordancy reflected actual molecular ITH within the tumors. The relevance of the subtypes was reflected at protein level where several inflammation markers were significantly increased in immune related transcriptional subtypes, which was verified in an independent cohort (Wilcoxon rank sum test; p<0.05). Unsupervised hierarchical clustering of the protein data identified large ITH at protein level; as the multiregional biopsies clustered together for only 9 out of 29 tumors., Conclusion: Our transcriptomic and proteomic analyses show that the tumor location along the colorectum influence the ITH of CRC, which again influence the concordance of subtyping., Competing Interests: The authors have declared that no competing interests exist.

Chromosomal instability in cancer consists of dynamic changes to the number and structure of chromosomes 1,2 . The resulting diversity in somatic copy number alterations (SCNAs) may provide the variation necessary for tumour evolution 1,3,4 . Here we use multi-sample phasing and SCNA analysis of 1,421 samples from 394 tumours across 22 tumour types to show that continuous chromosomal instability results in pervasive SCNA heterogeneity. Parallel evolutionary events, which cause disruption in the same genes (such as BCL9, MCL1, ARNT (also known as HIF1B), TERT and MYC) within separate subclones, were present in 37% of tumours. Most recurrent losses probably occurred before whole-genome doubling, that was found as a clonal event in 49% of tumours. However, loss of heterozygosity at the human leukocyte antigen (HLA) locus and loss of chromosome 8p to a single haploid copy recurred at substantial subclonal frequencies, even in tumours with whole-genome doubling, indicating ongoing karyotype remodelling. Focal amplifications that affected chromosomes 1q21 (which encompasses BCL9, MCL1 and ARNT), 5p15.33 (TERT), 11q13.3 (CCND1), 19q12 (CCNE1) and 8q24.1 (MYC) were frequently subclonal yet appeared to be clonal within single samples. Analysis of an independent series of 1,024 metastatic samples revealed that 13 focal SCNAs were enriched in metastatic samples, including gains in chromosome 8q24.1 (encompassing MYC) in clear cell renal cell carcinoma and chromosome 11q13.3 (encompassing CCND1) in HER2 + breast cancer. Chromosomal instability may enable the continuous selection of SCNAs, which are established as ordered events that often occur in parallel, throughout tumour evolution.

Although thousands of solid tumors have been sequenced to date, a fundamental under-sampling bias is inherent in current methodologies. This is caused by a tissue sample input of fixed dimensions (e.g., 6 mm biopsy), which becomes grossly under-powered as tumor volume scales. Here, we demonstrate representative sequencing (Rep-Seq) as a new method to achieve unbiased tumor tissue sampling. Rep-Seq uses fixed residual tumor material, which is homogenized and subjected to next-generation sequencing. Analysis of intratumor tumor mutation burden (TMB) variability shows a high level of misclassification using current single-biopsy methods, with 20% of lung and 52% of bladder tumors having at least one biopsy with high TMB but low clonal TMB overall. Misclassification rates by contrast are reduced to 2% (lung) and 4% (bladder) when a more representative sampling methodology is used. Rep-Seq offers an improved sampling protocol for tumor profiling, with significant potential for improved clinical utility and more accurate deconvolution of clonal structure., Competing Interests: Declaration of Interests S.T. reports grants from Roche Tissue Diagnostics. K.L., S.T., and C.S. report speaker fees from Roche Tissue Diagnostics. K.L., S.S., S.T., and N.R.A. report a patent pending using Rep-Seq profiling to support detection of MRD. N.R.A., S.S., L.G., A.B., K.F.L., and S.H. have patents pending on representative sampling of solid tumors. C.S. receives grant support from Pfizer, AstraZeneca (AZ), Bristol-Myers Squibb (BMS), Roche-Ventana, Boehringer Ingelheim, and Ono. C.S. has consulted for Pfizer, Novartis, GlaxoSmithKline (GSK), Merck Sharpe & Dohme (MSD), BMS, Celgene, AZ, Illumina, Genentech, Roche-Ventana, GRAIL, Medicxi, and the Sarah Cannon Research Institute. C.S. is a shareholder of Apogen Biotechnologies, Epic Bioscience, and GRAIL and has stock options in and is co-founder of Achilles Therapeutics. Outside of the submitted work, K.L., S.T., and C.S. have a patent pending on indel burden and CPI response and a patent pending on targeting of frameshift neoantigens for personalized immunotherapy. J.L. reports institutional research support from BMS, MSD, Novartis, Pfizer, Achilles Therapeutics, Roche, Nektar Therapeutics, Covance, Immunocore, Pharmacyclics, and Aveo and consultancy support from Achilles, AZ, Boston Biomedical, BMS, Eisai, EUSA Pharma, GSK, Ipsen, Imugene, Incyte, iOnctura, Kymab, Merck Serono, MSD, Nektar, Novartis, Pierre Fabre, Pfizer, Roche/Genentech, Secarna, and Vitaccess. S.T. has received speaking fees from Astra Zeneca, Novartis, and Ipsen. S.T., K.L., and C.S. have the following patents filed: Clear Cell Renal Cell Carcinoma BiomarkersP113326GB., (Copyright © 2020 The Authors. Published by Elsevier Inc. All rights reserved.)

Tumour mutational burden (TMB) predicts immunotherapy outcome in non-small cell lung cancer (NSCLC), consistent with immune recognition of tumour neoantigens. However, persistent antigen exposure is detrimental for T cell function. How TMB affects CD4 and CD8 T cell differentiation in untreated tumours, and whether this affects patient outcomes is unknown. Here we paired high-dimensional flow cytometry, exome, single-cell and bulk RNA sequencing from patients with resected, untreated NSCLC to examine these relationships. TMB was associated with compartment-wide T cell differentiation skewing, characterized by loss of TCF7-expressing progenitor-like CD4 T cells, and an increased abundance of dysfunctional CD8 and CD4 T cell subsets, with significant phenotypic and transcriptional similarity to neoantigen-reactive CD8 T cells. A gene signature of redistribution from progenitor-like to dysfunctional states associated with poor survival in lung and other cancer cohorts. Single-cell characterization of these populations informs potential strategies for therapeutic manipulation in NSCLC., Competing Interests: Competing interest statement S.A.Q., K.S.P. and C.S. are co-founders of Achilles Therapeutics. C.S. receives grant support from Pfizer, AstraZeneca, BMS, Roche-Ventana and Boehringer-Ingelheim. C.S. has consulted for Pfizer, Novartis, GlaxoSmithKline, MSD, BMS, Celgene, AstraZeneca, Illumina, Genentech, Roche-Ventana, GRAIL, Medicxi, and the Sarah Cannon Research Institute. C.S. is a shareholder of Apogen Biotechnologies, Epic Bioscience, GRAIL, and has stock options in and is co-founder of Achilles Therapeutics. R.R., N.M. and G.A.W. have stock options in and have consulted for Achilles Therapeutics. J.L.R and M.A.B have consulted for Achilles Therapeutics. P.D.B and M.W.S are employees of Achilles Therapeutics.

Whole-genome doubling (WGD) is a prevalent event in cancer, involving a doubling of the entire chromosome complement. However, despite its prevalence and prognostic relevance, the evolutionary selection pressures for WGD in cancer have not been investigated. Here, we combine evolutionary simulations with an analysis of cancer sequencing data to explore WGD during cancer evolution. Simulations suggest that WGD can be selected to mitigate the irreversible, ratchet-like, accumulation of deleterious somatic alterations, provided that they occur at a sufficiently high rate. Consistent with this, we observe an enrichment for WGD in tumor types with extensive loss of heterozygosity, including lung squamous cell carcinoma and triple-negative breast cancers, and we find evidence for negative selection against homozygous loss of essential genes before, but not after, WGD. Finally, we demonstrate that loss of heterozygosity and temporal dissection of mutations can be exploited to identify novel tumor suppressor genes and to obtain a deeper characterization of known cancer genes.

Metastatic cancer is a major cause of death and remains largely incurable. A better understanding of metastasis is therefore desperately needed to improve prognosis for late-stage disease. Here we survey the landscape of studies exploring the genomics of metastatic cancer. We consider evidence for genomic drivers of metastasis and explore studies investigating modes of metastatic spread., (Copyright © 2019 Elsevier Inc. All rights reserved.)

Somatic mutations together with immunoediting drive extensive heterogeneity within non-small-cell lung cancer (NSCLC). Herein we examine heterogeneity of the T cell antigen receptor (TCR) repertoire. The number of TCR sequences selectively expanded in tumors varies within and between tumors and correlates with the number of nonsynonymous mutations. Expanded TCRs can be subdivided into TCRs found in all tumor regions (ubiquitous) and those present in a subset of regions (regional). The number of ubiquitous and regional TCRs correlates with the number of ubiquitous and regional nonsynonymous mutations, respectively. Expanded TCRs form part of clusters of TCRs of similar sequence, suggestive of a spatially constrained antigen-driven process. CD8 + tumor-infiltrating lymphocytes harboring ubiquitous TCRs display a dysfunctional tissue-resident phenotype. Ubiquitous TCRs are preferentially detected in the blood at the time of tumor resection as compared to routine follow-up. These findings highlight a noninvasive method to identify and track relevant tumor-reactive TCRs for use in adoptive T cell immunotherapy.

An aim of molecular biomarkers is to stratify patients with cancer into disease subtypes predictive of outcome, improving diagnostic precision beyond clinical descriptors such as tumor stage 1 . Transcriptomic intratumor heterogeneity (RNA-ITH) has been shown to confound existing expression-based biomarkers across multiple cancer types 2-6 . Here, we analyze multi-region whole-exome and RNA sequencing data for 156 tumor regions from 48 patients enrolled in the TRACERx study to explore and control for RNA-ITH in non-small cell lung cancer. We find that chromosomal instability is a major driver of RNA-ITH, and existing prognostic gene expression signatures are vulnerable to tumor sampling bias. To address this, we identify genes expressed homogeneously within individual tumors that encode expression modules of cancer cell proliferation and are often driven by DNA copy-number gains selected early in tumor evolution. Clonal transcriptomic biomarkers overcome tumor sampling bias, associate with survival independent of clinicopathological risk factors, and may provide a general strategy to refine biomarker design across cancer types.

Introduction: Immunotherapy is regarded as one of the major breakthroughs in cancer treatment. Despite its success, only a subset of patients responds-urging the quest for predictive biomarkers. We hypothesize that artificial intelligence (AI) algorithms can automatically quantify radiographic characteristics that are related to and may therefore act as noninvasive radiomic biomarkers for immunotherapy response., Patients and Methods: In this study, we analyzed 1055 primary and metastatic lesions from 203 patients with advanced melanoma and non-small-cell lung cancer (NSCLC) undergoing anti-PD1 therapy. We carried out an AI-based characterization of each lesion on the pretreatment contrast-enhanced CT imaging data to develop and validate a noninvasive machine learning biomarker capable of distinguishing between immunotherapy responding and nonresponding. To define the biological basis of the radiographic biomarker, we carried out gene set enrichment analysis in an independent dataset of 262 NSCLC patients., Results: The biomarker reached significant performance on NSCLC lesions (up to 0.83 AUC, P < 0.001) and borderline significant for melanoma lymph nodes (0.64 AUC, P = 0.05). Combining these lesion-wide predictions on a patient level, immunotherapy response could be predicted with an AUC of up to 0.76 for both cancer types (P < 0.001), resulting in a 1-year survival difference of 24% (P = 0.02). We found highly significant associations with pathways involved in mitosis, indicating a relationship between increased proliferative potential and preferential response to immunotherapy., Conclusions: These results indicate that radiographic characteristics of lesions on standard-of-care imaging may function as noninvasive biomarkers for response to immunotherapy, and may show utility for improved patient stratification in both neoadjuvant and palliative settings., (© The Author(s) 2019. Published by Oxford University Press on behalf of the European Society for Medical Oncology.)

The interplay between an evolving cancer and a dynamic immune microenvironment remains unclear. Here we analyse 258 regions from 88 early-stage, untreated non-small-cell lung cancers using RNA sequencing and histopathology-assessed tumour-infiltrating lymphocyte estimates. Immune infiltration varied both between and within tumours, with different mechanisms of neoantigen presentation dysfunction enriched in distinct immune microenvironments. Sparsely infiltrated tumours exhibited a waning of neoantigen editing during tumour evolution, indicative of historical immune editing, or copy-number loss of previously clonal neoantigens. Immune-infiltrated tumour regions exhibited ongoing immunoediting, with either loss of heterozygosity in human leukocyte antigens or depletion of expressed neoantigens. We identified promoter hypermethylation of genes that contain neoantigenic mutations as an epigenetic mechanism of immunoediting. Our results suggest that the immune microenvironment exerts a strong selection pressure in early-stage, untreated non-small-cell lung cancers that produces multiple routes to immune evasion, which are clinically relevant and forecast poor disease-free survival.

Judgement, as one of the core tenets of medicine, relies upon the integration of multilayered data with nuanced decision making. Cancer offers a unique context for medical decisions given not only its variegated forms with evolution of disease but also the need to take into account the individual condition of patients, their ability to receive treatment, and their responses to treatment. Challenges remain in the accurate detection, characterization, and monitoring of cancers despite improved technologies. Radiographic assessment of disease most commonly relies upon visual evaluations, the interpretations of which may be augmented by advanced computational analyses. In particular, artificial intelligence (AI) promises to make great strides in the qualitative interpretation of cancer imaging by expert clinicians, including volumetric delineation of tumors over time, extrapolation of the tumor genotype and biological course from its radiographic phenotype, prediction of clinical outcome, and assessment of the impact of disease and treatment on adjacent organs. AI may automate processes in the initial interpretation of images and shift the clinical workflow of radiographic detection, management decisions on whether or not to administer an intervention, and subsequent observation to a yet to be envisioned paradigm. Here, the authors review the current state of AI as applied to medical imaging of cancer and describe advances in 4 tumor types (lung, brain, breast, and prostate) to illustrate how common clinical problems are being addressed. Although most studies evaluating AI applications in oncology to date have not been vigorously validated for reproducibility and generalizability, the results do highlight increasingly concerted efforts in pushing AI technology to clinical use and to impact future directions in cancer care., (© 2019 American Cancer Society.)

Circulating tumour DNA (ctDNA) refers to the fraction of cell-free DNA in a patient's blood that originates from a tumour. Advances in DNA sequencing technologies and our understanding of the molecular biology of tumours have resulted in increased interest in exploiting ctDNA as a tool to facilitate earlier detection of cancer and thereby improve therapeutic outcomes by enabling early intervention. ctDNA analysis might also have utility in the adjuvant therapeutic setting by enabling the identification of patients at a high risk of disease recurrence on the basis of the detection of post-surgical minimal (or molecular) residual disease (MRD). This approach could provide the capability to adapt clinical trials in the adjuvant setting in order to optimize risk stratification, and we argue that this objective is achievable with current technologies. Herein, we evaluate contemporary next-generation sequencing (NGS) approaches to ctDNA detection with a focus on non-small-cell lung cancer. We explain the technical and analytical challenges to low-frequency mutation detection using NGS-based ctDNA profiling and evaluate the feasibility of ctDNA profiling in both screening and MRD assessment contexts.

The original version of this article incorrectly omitted an affiliation of Patricia A. J. Muller: 'Cancer Research UK Manchester Institute, The University of Manchester | Alderley Park, Manchester, SK10 4TG, UK'. This has been corrected in both the PDF and HTML versions of the Article.

Cell-in-cell (CIC) structures are commonly seen in tumours. Their biological significance remains unclear, although they have been associated with more aggressive tumours. Here we report that mutant p53 promotes CIC via live cell engulfment. Engulfed cells physically interfere in cell divisions of host cells and for cells without p53 this leads to host cell death. In contrast, mutant p53 host cells survive, display aberrant divisions, multinucleation and tripolar mitoses. In xenograft studies, CIC-rich p53 mutant/null co-cultures show enhanced tumour growth. Furthermore, our results show that CIC is common within lung adenocarcinomas, is an independent predictor of poor outcome and disease recurrence, is associated with mutant p53 expression and correlated to measures of heterogeneity and genomic instability. These findings suggest that pro-tumorigenic entotic engulfment activity is associated with mutant p53 expression, and the two combined are a key factor in genomic instability.

The first genomic scar-based homologous recombination deficiency (HRD) measures were produced using SNP arrays. As array-based technology has been largely replaced by next generation sequencing approaches, it has become important to develop algorithms that derive the same type of genomic scar scores from next generation sequencing (whole exome "WXS", whole genome "WGS") data. In order to perform this analysis, we introduce here the scarHRD R package and show that using this method the SNP array-based and next generation sequencing-based derivation of HRD scores show good correlation (Pearson correlation between 0.73 and 0.87 depending on the actual HRD measure) and that the NGS-based HRD scores distinguish similarly well between BRCA mutant and BRCA wild-type cases in a cohort of triple-negative breast cancer patients of the TCGA data set., Competing Interests: N.J.B., A.C.E., and Zo.S are listed as co-inventors on a patent on telomeric allelic imbalance, which is owned by Children’s Hospital Boston and licensed to Myriad Genetics. The remaining authors declare no competing interests.