Clonal status of actionable driver events and the timing of mutational processes in cancer evolution

The evolutionary history of metastatic pancreatic neuroendocrine tumours reveals a therapy driven route to high-grade transformation

Tumour evolution with acquisition of more aggressive disease characteristics is a hallmark of disseminated cancer. Metastatic pancreatic neuroendocrine tumours (PanNETs) in particular may progress from a low/intermediate to a high-grade disease. The aim of this work was to understand the molecular mechanisms underlying metastatic progression as well as PanNET transformation from a low/intermediate to a high-grade disease. We performed multi-omics analysis (genome/exome sequencing, total RNA-sequencing and methylation array) of 32 longitudinal samples from six patients with metastatic low/intermediate grade PanNET. The clonal composition of tumour lesions and underlying phylogeny of each patient were determined with bioinformatics analyses. Findings were validated in post-alkylating chemotherapy samples from 24 patients with PanNET using targeted next generation sequencing. We validate the current PanNET evolutionary model with MEN1 inactivation that occurs very early in tumourigenesis. This was followed by pronounced genetic diversity on both spatial and temporal levels, with parallel and convergent tumour evolution involving the ATRX/DAXX and mechanistic target of the rapamycin (mTOR) pathways. Following alkylating chemotherapy treatment, some PanNETs developed mismatch repair deficiency and acquired a hypermutational phenotype. This was validated among 16 patients with PanNET who had high-grade progression after alkylating chemotherapy, of whom eight had a tumour mutational burden >50 (50%). In comparison, among the eight patients who did not show high-grade progression, 0 had a tumour mutational burden >50 (0%; odds ratio 'infinite', 95% confidence interval 1.8 to 'infinite', p = 0.02). Our findings contribute to broaden the understanding of metastatic/high-grade PanNETs and suggests that therapy driven disease evolution is an important hallmark of this disease. © 2024 The Author(s). The Journal of Pathology published by John Wiley & Sons Ltd on behalf of The Pathological Society of Great Britain and Ireland. This article has been contributed to by U.S. Government employees and their work is in the public domain in the USA.

Molecular profiling reveals novel therapeutic targets and clonal evolution in ovarian clear cell carcinoma

BACKGROUND

Ovarian clear cell carcinoma (OCCC) has a disproportionately high incidence among women in East Asia. Patients diagnosed with OCCC tend to experience worse clinical outcomes than those with high-grade serous carcinoma (HGSC) at advanced stages. The unfavorable prognosis of OCCC can be partly attributed to its frequent resistance to conventional chemotherapy. Within a precision medicine framework, we sought to provide a comprehensive molecular characterization of OCCC using whole-exome sequencing to uncover potential molecular targets that may inform novel therapeutic strategies.

METHODS

We performed whole-exome sequencing analysis on tumor-normal paired samples from 102 OCCC patients. This comprehensive genomic characterization of a substantial cohort of OCCC specimens was coupled with an analysis of clonal progression.

RESULTS

On analyzing 102 OCCC samples, ARID1A (67%) and PIK3CA (49%) emerged as the most frequently mutated driver genes. We identified tier 1 or 2 clinically actionable molecular targets in 40% of cases. This included DNA mismatch repair deficiency (n = 1), as well as BRCA2 (n = 1), PIK3CA (n = 36), KRASG12C (n = 1), and ATM (n = 4) mutations. Furthermore, 45% of OCCC samples displayed ARID1A biallelic loss. Interestingly, we identified previously unreported mutations in the 5' untranslated region of the TERT gene that harbored an adverse prognostic significance. Clock-like mutational processes and activated APOBECs were major drivers of somatic point mutations. Mutations arising from DNA mismatch repair deficiency were uncommon. Reconstruction of clonal evolution revealed that early genetic events likely driving tumorigenesis included mutations in the ARID1A, PIK3CA, TERT, KRAS, and TP53 genes.

CONCLUSIONS

Our study provides a comprehensive characterization of the genomic landscape and clonal evolution in OCCC within a substantial cohort. These findings unveil potentially actionable molecular alterations that could be leveraged to develop targeted therapies.

Genomic and transcriptomic landscape of human gastrointestinal stromal tumors

Gastrointestinal stromal tumor (GISTs) are clinically heterogenous exhibiting varying degrees of disease aggressiveness in individual patients. We comprehensively describe the genomic and transcriptomic landscape of a cohort of 117 GISTs including 31 low-risk, 18 intermediate-risk, 29 high-risk, 34 metastatic and 5 neoadjuvant GISTs from 105 patients. GISTs have notably low tumor mutation burden but widespread copy number variations. Aggressive GISTs harbor remarkably more genomic aberrations than low-/intermediate-risk GISTs. Complex genomic alterations, chromothripsis and kataegis, occur selectively in aggressive GISTs. Despite the paucity of mutations, recurrent inactivating YLPM1 mutations are identified (10.3%, 7 of 68 patients), enriched in high-risk/metastatic GIST and functional study further demonstrates YLPM1 inactivation promotes GIST proliferation, growth and oxidative phosphorylation. Spatially and temporally separated GISTs from individual patients demonstrate complex tumor heterogeneity in metastatic GISTs. Finally, four prominent subtypes are proposed with different genomic features, expression profiles, immune characteristics, clinical characteristics and subtype-specific treatment strategies. This large-scale analysis depicts the landscape and provides further insights into GIST pathogenesis and precise treatment.

Multiomics Profiling Distinguishes Sebaceous Carcinoma from Benign Sebaceous Neoplasms and Provides Insight into the Genetic Evolution of Sebaceous Carcinogenesis

PURPOSE

Sebaceous carcinoma is the third most common nonkeratinocyte skin cancer in the United States with 1,000 cases per year. The clinicopathologic features of sebaceous carcinoma and benign sebaceous neoplasms (adenomas, sebaceomas) can overlap, highlighting the need for molecular biomarkers to improve classification. This study describes the genomic and transcriptomic landscape of sebaceous neoplasms in order to understand tumor etiology and biomarkers relevant for diagnosis and treatment.

EXPERIMENTAL DESIGN

We performed whole-genome sequencing (WGS) and whole-transcriptome sequencing (WTS) of sebaceous neoplasms from six academic and two federal healthcare facilities in the United States diagnosed between January 1, 1999, and December 31, 2021.

RESULTS

We evaluated 98 sebaceous neoplasms: 64 tumors (32 adenomas, 2 sebaceomas, 5 atypical sebaceous neoplasms, 25 carcinomas) had sufficient material for WGS, 96 tumors (42 adenomas, 11 sebaceomas, 8 atypical sebaceous neoplasms, 35 carcinomas) had sufficient material for WTS, and 62 tumors (31 adenomas, 2 sebaceomas, 5 atypical sebaceous neoplasms, 24 carcinomas) had sufficient material for combined WGS and WTS. Overall, we found decreased cholesterol biosynthesis and increased TP53 mutations, copy number gains (chromosome 6, 8q, and/or 18), and tumor mutation burden-high (>10 mutations/MB) in carcinomas compared to adenomas. Although diminished compared to adenomas, most carcinomas still had higher cholesterol biosynthesis than nonmalignant skin. Multiomics profiling also supported a precancerous model of tumor evolution with sebaceomas and atypical sebaceous neoplasms being likely intermediate lesions.

CONCLUSIONS

The study findings highlight key diagnostic biomarkers for sebaceous carcinoma and suggest that immunotherapy and modulation of cholesterol biosynthesis could be effective treatment strategies.

Whole-exome profiles of inflammatory breast cancer and pathological response to neoadjuvant chemotherapy

BACKGROUND

Neoadjuvant chemotherapy (NACT) became a standard treatment strategy for patients with inflammatory breast cancer (IBC) because of high disease aggressiveness. However, given the heterogeneity of IBC, no molecular feature reliably predicts the response to chemotherapy. Whole-exome sequencing (WES) of clinical tumor samples provides an opportunity to identify genomic alterations associated with chemosensitivity.

METHODS

We retrospectively applied WES to 44 untreated IBC primary tumor samples and matched normal DNA. The pathological response to NACT, assessed on operative specimen, distinguished the patients with versus without pathological complete response (pCR versus no-pCR respectively). We compared the mutational profiles, spectra and signatures, pathway mutations, copy number alterations (CNAs), HRD, and heterogeneity scores between pCR versus no-pCR patients.

RESULTS

The TMB, HRD, and mutational spectra were not different between the complete (N = 13) versus non-complete (N = 31) responders. The two most frequently mutated genes were TP53 and PIK3CA. They were more frequently mutated in the complete responders, but the difference was not significant. Only two genes, NLRP3 and SLC9B1, were significantly more frequently mutated in the complete responders (23% vs. 0%). By contrast, several biological pathways involved in protein translation, PI3K pathway, and signal transduction showed significantly higher mutation frequency in the patients with pCR. We observed a higher abundance of COSMIC signature 7 (due to ultraviolet light exposure) in tumors from complete responders. The comparison of CNAs of the 3808 genes included in the GISTIC regions between both patients' groups identified 234 genes as differentially altered. The CIN signatures were not differentially represented between the complete versus non-complete responders. Based on the H-index, the patients with heterogeneous tumors displayed a lower pCR rate (11%) than those with less heterogeneous tumors (35%).

CONCLUSIONS

This is the first study aiming at identifying correlations between the WES data of IBC samples and the achievement of pCR to NACT. Our results, obtained in this 44-sample series, suggest a few subtle genomic alterations associated with pathological response. Additional investigations are required in larger series.

Comprehensive pan-cancer analysis of FUTs family as prognostic and immunity markers based on multi-omics data

BACKGROUND

The dysregulation of fucosyltransferases (FUTs) contributes to alterations in fucosylated epitope expression, which serve as distinctive features of cancer cells. Nonetheless, a comprehensive elucidation of the prognostic biological marker and therapeutic target of the FUTs family in pan-cancer remains elusive.

METHODS

Over 10,000 individuals' profiling information was examined, including information on 750 small molecule drugs, 33 types of cancer, and 24 types of immune cells. We focused on POFUT2's function and applied GSVA (Gene Set Variation Analysis) to calculate the FUT score. Survival and cancer pathways were found to be correlated with this score. After deriving a signature via univariate Cox and LASSO regression, we generated and analyzed the ROC curve and developed a nomogram.

RESULTS

Our comprehensive analysis revealed epigenetic, genomic, and immunogenomic changes in FUTs, particularly POFUT2, resulting in aberrant expression. Elevated frequencies of CNV (Copy number variation), SNV (Single Nucleotide Variant), and hypermethylation were observed in FUTs. Additionally, the survival of patients with various types of cancers may be predicted by FUT expression. Immune response and prognosis in numerous types of cancer were found to be strongly linked to aberrant POFUT2 expression. Pathway analysis unveiled the role of FUTs in apoptosis, epithelial-to-mesenchymal transition (EMT), cell cycle, DNA damage response, RAS/MAPK, TSC/mTOR, PI3K/AKT, AR, ER, and RTK. A prognostic index for patients diagnosed with adrenocortical carcinoma (ACC) was established by applying a risk model incorporating nine FUTs and based on the findings of the GSVA.

CONCLUSIONS

FUTs, particularly POFUT2, emerge as candidate targets for improving the outcomes of immune therapy. The significance of aberrant MUC12 expression, cancer immune therapy, and patient survival in the context of diverse malignancies is enhanced by the strong correlation observed among these factors. Our five-gene risk signature provides patients with ACC with an independent prognostic indicator, emphasizing the critical function of these genes in inhibiting the immune system's response in ACC.

Longitudinal Comparative Analysis of Circulating Tumor DNA and Matched Tumor Tissue DNA in Patients with Metastatic Colorectal Cancer Receiving Palliative First-Line Systemic Anti-Cancer Therapy

PURPOSE

This study aimed to compare tumor tissue DNA (ttDNA) and circulating tumor DNA (ctDNA) to explore the clinical applicability of ctDNA and to better understand clonal evolution in patients with metastatic colorectal cancer undergoing palliative first-line systemic therapy.

MATERIALS AND METHODS

We performed targeted sequencing analysis of 88 cancer-associated genes using germline DNA, ctDNA at baseline (baseline-ctDNA), and ctDNA at progressive disease (PD-ctDNA). The results were compared with ttDNA data.

RESULTS

Among 208 consecutively enrolled patients, we selected 84 (41 males; median age, 59 years; range, 35 to 90 years) with all four sample types available. A total of 202 driver mutations were found in 34 genes. ttDNA exhibited the highest mutation frequency (n=232), followed by baseline-ctDNA (n=155) and PD-ctDNA (n=117). Sequencing ctDNA alongside ttDNA revealed additional mutations in 40 patients (47.6%). PD-ctDNA detected 13 novel mutations in 10 patients (11.9%) compared to ttDNA and baseline-ctDNA. Notably, seven mutations in five patients (6.0%) were missense or nonsense mutations in APC, TP53, SMAD4, and CDH1 genes. In baseline-ctDNA, higher maximal variant allele frequency (VAF) values (p=0.010) and higher VAF values of APC (p=0.012), TP53 (p=0.012), and KRAS (p=0.005) mutations were significantly associated with worse overall survival.

CONCLUSION

While ttDNA remains more sensitive than ctDNA, our ctDNA platform demonstrated validity and potential value when ttDNA was unavailable. Post-treatment analysis of PD-ctDNA unveiled new pathogenic mutations, signifying cancer's clonal evolution. Additionally, baseline-ctDNA's VAF values were prognostic after treatment.

A unique interplay of access and selection shapes peritoneal metastasis evolution in colorectal cancer

Whether metastasis in humans can be accomplished by most primary tumor cells or requires the evolution of a specialized trait remains an open question. To evaluate whether metastases are founded by non-random subsets of primary tumor lineages requires extensive, difficult-to-implement sampling. We have realized an unusually dense multi-region sampling scheme in a cohort of 26 colorectal cancer patients with peritoneal metastases, reconstructing the evolutionary history of on average 28.8 tissue samples per patient with a microsatellite-based fingerprinting assay. To assess metastatic randomness, we evaluate inter- and intra-metastatic heterogeneity relative to the primary tumor and find that peritoneal metastases are more heterogeneous than liver metastases but less diverse than locoregional metastases. Metachronous peritoneal metastases exposed to systemic chemotherapy show significantly higher inter-lesion diversity than synchronous, untreated metastases. Projection of peritoneal metastasis origins onto a spatial map of the primary tumor reveals that they often originate at the deep-invading edge, in contrast to liver and lymph node metastases which exhibit no such preference. Furthermore, peritoneal metastases typically do not share a common subclonal origin with distant metastases in more remote organs. Synthesizing these insights into an evolutionary portrait of peritoneal metastases, we conclude that the peritoneal-metastatic process imposes milder selective pressures onto disseminating cancer cells than the liver-metastatic process. Peritoneal metastases' unique evolutionary features have potential implications for staging and treatment.

Exploring evolutionary trajectories in ovarian cancer patients by longitudinal analysis of ctDNA

OBJECTIVES

We analysed whether temporal heterogeneity of ctDNA encodes evolutionary patterns in ovarian cancer.

METHODS

Targeted sequencing of 275 cancer-associated genes was performed in a primary tumor biopsy and in ctDNA of six longitudinal plasma samples from 15 patients, using the Illumina platform.

RESULTS

While there was low overall concordance between the mutational spectrum of the primary tumor biopsies vs. ctDNA, TP53 variants were the most commonly shared somatic alterations. Up to three variant clusters were detected in each tumor biopsy, likely representing predominant clones of the primary tumor, most of them harbouring a TP53 variant. By tracing these clusters in ctDNA, we propose that liquid biopsy may allow to assess the contribution of ancestral clones of the tumor to relapsed abdominal masses, revealing two evolutionary patterns. In pattern#1, clusters detected in the primary tumor biopsy were likely relapse seeding clones, as they contributed a major share to ctDNA at relapse. In pattern#2, similar clusters were present in tumors and ctDNA; however, they were entirely cleared from liquid biopsy after chemotherapy and were undetectable at relapse. ctDNA private variants were present among both patterns, with some of them mirroring subclonal expansions after chemotherapy.

CONCLUSIONS

We demonstrate that tracing the temporal heterogeneity of ctDNA, even below exome scale resolution, deciphers evolutionary trajectories in ovarian cancer. Furthermore, we describe two evolutionary patterns that may help to identify relapse seeding clones for targeted therapy.

Multi-regional genomic and transcriptomic characterization of a melanoma-associated oral cavity cancer provide evidence for CASP8 alteration-mediated field cancerization

BACKGROUND

Precancerous and malignant tumours arise within the oral cavity from a predisposed "field" of epithelial cells upon exposure to carcinogenic stimulus. This phenomenon is known as "Field Cancerization". The molecular genomic and transcriptomic alterations that lead to field cancerization and tumour progression is unknown in Indian Oral squamous cell carcinoma (OSCC) patients.

METHODS

We have performed whole exome sequencing, copy-number variation array and whole transcriptome sequencing from five tumours and dysplastic lesions (sampled from distinct anatomical subsites - one each from buccal anterior and posterior alveolus, dorsum of tongue-mucosal melanoma, lip and left buccal mucosa) and blood from a rare OSCC patient with field cancerization.

RESULTS

A missense CASP8 gene mutation (p.S375F) was observed to be the initiating event in oral tumour field development. APOBEC mutation signatures, arm-level copy number alterations, depletion of CD8 + T cells and activated NK cells and enrichment of pro-inflammatory mast cells were features of early-originating tumours. Pharmacological inhibition of CASP8 protein in a CASP8-wild type OSCC cell line showed enhanced levels of cellular migration and viability.

CONCLUSION

CASP8 alterations are the earliest driving events in oral field carcinogenesis, whereas additional somatic mutational, copy number and transcriptomic alterations ultimately lead to OSCC tumour formation and progression.

Microsatellite Instability, Tumor Mutational Burden, and Response to Immune Checkpoint Blockade in Patients with Prostate Cancer

PURPOSE

Patients with microsatellite instability-high/mismatch repair-deficient (MSI-H/dMMR) and high tumor mutational burden (TMB-H) prostate cancers are candidates for pembrolizumab. We define the genomic features, clinical course, and response to immune checkpoint blockade (ICB) in patients with MSI-H/dMMR and TMB-H prostate cancers without MSI [TMB-H/microsatellite stable (MSS)].

EXPERIMENTAL DESIGN

We sequenced 3,244 tumors from 2,257 patients with prostate cancer. MSI-H/dMMR prostate cancer was defined as an MSIsensor score ≥10 or MSIsensor score ≥3 and <10 with a deleterious MMR alteration. TMB-H was defined as ≥10 mutations/megabase. PSA50 and RECIST responses were assigned. Overall survival and radiographic progression-free survival (rPFS) were compared using log-rank test.

RESULTS

Sixty-three (2.8%) men had MSI-H/dMMR, and 33 (1.5%) had TMB-H/MSS prostate cancers. Patients with MSI-H/dMMR and TMB-H/MSS tumors more commonly presented with grade group 5 and metastatic disease at diagnosis. MSI-H/dMMR tumors had higher TMB, indel, and neoantigen burden compared with TMB-H/MSS. Twenty-seven patients with MSI-H/dMMR and 8 patients with TMB-H/MSS tumors received ICB, none of whom harbored polymerase epsilon (polE) catalytic subunit mutations. About 45% of patients with MSI-H/dMMR had a RECIST response, and 65% had a PSA50 response. No patient with TMB-H/MSS had a RECIST response, and 50% had a PSA50 response. rPFS tended to be longer in patients with MSI-H/dMMR than in patients with TMB-H/MSS who received immunotherapy. Pronounced differences in genomics, TMB, or MSIsensor score were not detected between MSI-H/dMMR responders and nonresponders.

CONCLUSIONS

MSI-H/dMMR prostate cancers have greater TMB, indel, and neoantigen burden than TMB-H/MSS prostate cancers, and these differences may contribute to profound and durable responses to ICB.

Regulation, functional impact, and therapeutic targeting of APOBEC3A in cancer

Enzymes of the apolipoprotein B mRNA editing catalytic polypeptide like (APOBEC) family are cytosine deaminases that convert cytosine to uracil in DNA and RNA. Among these proteins, APOBEC3 sub-family members, APOBEC3A (A3A) and APOBEC3B (A3B), are prominent sources of mutagenesis in cancer cells. The aberrant expression of A3A and A3B in cancer cells leads to accumulation of mutations with specific single-base substitution (SBS) signatures, characterized by C→T and C→G changes, in a number of tumor types. In addition to fueling mutagenesis, A3A and A3B, particularly A3A, induce DNA replication stress, DNA damage, and chromosomal instability through their catalytic activities, triggering a range of cellular responses. Thus, A3A/B have emerged as key drivers of genome evolution during cancer development, contributing to tumorigenesis, tumor heterogeneity, and therapeutic resistance. Yet, the expression of A3A/B in cancer cells presents a cancer vulnerability that can be exploited therapeutically. In this review, we discuss the recent studies that shed light on the mechanisms regulating A3A expression and the impact of A3A in cancer. We also review recent advances in the development of A3A inhibitors and provide perspectives on the future directions of A3A research.

The genomic trajectory of ovarian high-grade serous carcinoma can be observed in STIC lesions

Ovarian high-grade serous carcinoma (HGSC) originates in the fallopian tube, with secretory cells carrying a TP53 mutation, known as p53 signatures, identified as potential precursors. p53 signatures evolve into serous tubal intraepithelial carcinoma (STIC) lesions, which in turn progress into invasive HGSC, which readily spreads to the ovary and disseminates around the peritoneal cavity. We recently investigated the genomic landscape of early- and late-stage HGSC and found higher ploidy in late-stage (median 3.1) than early-stage (median 2.0) samples. Here, to explore whether the high ploidy and possible whole-genome duplication (WGD) observed in late-stage disease were determined early in the evolution of HGSC, we analysed archival formalin-fixed paraffin-embedded (FFPE) samples from five HGSC patients. p53 signatures and STIC lesions were laser-capture microdissected and sequenced using shallow whole-genome sequencing (sWGS), while invasive ovarian/fallopian tube and metastatic carcinoma samples underwent macrodissection and were profiled using both sWGS and targeted next-generation sequencing. Results showed highly similar patterns of global copy number change between STIC lesions and invasive carcinoma samples within each patient. Ploidy changes were evident in STIC lesions, but not p53 signatures, and there was a strong correlation between ploidy in STIC lesions and invasive ovarian/fallopian tube and metastatic samples in each patient. The reconstruction of sample phylogeny for each patient from relative copy number indicated that high ploidy, when present, occurred early in the evolution of HGSC, which was further validated by copy number signatures in ovarian and metastatic tumours. These findings suggest that aberrant ploidy, suggestive of WGD, arises early in HGSC and is detected in STIC lesions, implying that the trajectory of HGSC may be determined at the earliest stages of tumour development. © 2024 The Author(s). The Journal of Pathology published by John Wiley & Sons Ltd on behalf of The Pathological Society of Great Britain and Ireland.

Replication timing alterations are associated with mutation acquisition during breast and lung cancer evolution

During each cell cycle, the process of DNA replication timing is tightly regulated to ensure the accurate duplication of the genome. The extent and significance of alterations in this process during malignant transformation have not been extensively explored. Here, we assess the impact of altered replication timing (ART) on cancer evolution by analysing replication-timing sequencing of cancer and normal cell lines and 952 whole-genome sequenced lung and breast tumours. We find that 6%-18% of the cancer genome exhibits ART, with regions with a change from early to late replication displaying an increased mutation rate and distinct mutational signatures. Whereas regions changing from late to early replication contain genes with increased expression and present a preponderance of APOBEC3-mediated mutation clusters and associated driver mutations. We demonstrate that ART occurs relatively early during cancer evolution and that ART may have a stronger correlation with mutation acquisition than alterations in chromatin structure.

Spatial intra-tumour heterogeneity and treatment-induced genomic evolution in oesophageal adenocarcinoma: implications for prognosis and therapy

BACKGROUND

Oesophageal adenocarcinoma (OAC) is a highly heterogeneous cancer with poor survival. Standard curative treatment is chemotherapy with or without radiotherapy followed by oesophagectomy. Genomic heterogeneity is a feature of OAC and has been linked to treatment resistance.

METHODS

Whole-genome sequencing data from 59 treatment-naïve and 18 post-treatment samples from 29 OAC patients was analysed. Twenty-seven of these were enrolled in the DOCTOR trial, sponsored by the Australasian Gastro-Intestinal Trials Group. Two biopsies from each treatment-naïve tumour were assessed to define 'shared' (between both samples) and 'private' (present in one sample) mutations.

RESULTS

Mutational signatures SBS2/13 (APOBEC) and SBS3 (BRCA) were almost exclusively detected in private mutation populations of treatment-naïve tumours. Patients presenting these signatures had significantly worse disease specific survival. Furthermore, mutational signatures associated with platinum-based chemotherapy treatment as well as high platinum enrichment scores were only detected in post-treatment samples. Additionally, clones with high putative neoantigen binding scores were detected in some treatment-naïve samples suggesting immunoediting of clones.

CONCLUSIONS

This study demonstrates the high intra-tumour heterogeneity in OAC, as well as indicators for treatment-induced changes during tumour evolution. Intra-tumour heterogeneity remains a problem for successful treatment strategies in OAC.

Genetic and epigenetic instability as an underlying driver of progression and aggressive behavior in IDH-mutant astrocytoma

In recent years, the classification of adult-type diffuse gliomas has undergone a revolution, wherein specific molecular features now represent defining diagnostic criteria of IDH-wild-type glioblastomas, IDH-mutant astrocytomas, and IDH-mutant 1p/19q-codeleted oligodendrogliomas. With the introduction of the 2021 WHO CNS classification, additional molecular alterations are now integrated into the grading of these tumors, given equal weight to traditional histologic features. However, there remains a great deal of heterogeneity in patient outcome even within these established tumor subclassifications that is unexplained by currently codified molecular alterations, particularly in the IDH-mutant astrocytoma category. There is also significant intercellular genetic and epigenetic heterogeneity and plasticity with resulting phenotypic heterogeneity, making these tumors remarkably adaptable and robust, and presenting a significant barrier to the design of effective therapeutics. Herein, we review the mechanisms and consequences of genetic and epigenetic instability, including chromosomal instability (CIN), microsatellite instability (MSI)/mismatch repair (MMR) deficits, and epigenetic instability, in the underlying biology, tumorigenesis, and progression of IDH-mutant astrocytomas. We also discuss the contribution of recent high-resolution transcriptomics studies toward defining tumor heterogeneity with single-cell resolution. While intratumoral heterogeneity is a well-known feature of diffuse gliomas, the contribution of these various processes has only recently been considered as a potential driver of tumor aggressiveness. CIN has an independent, adverse effect on patient survival, similar to the effect of histologic grade and homozygous CDKN2A deletion, while MMR mutation is only associated with poor overall survival in univariate analysis but is highly correlated with higher histologic/molecular grade and other aggressive features. These forms of genomic instability, which may significantly affect the natural progression of these tumors, response to therapy, and ultimately clinical outcome for patients, are potentially measurable features which could aid in diagnosis, grading, prognosis, and development of personalized therapeutics.

Evaluating the expression pattern of ATXN1 and CDC42EP1 genes and related long noncoding RNAs in oral squamous cell carcinoma

BACKGROUND

Oral squamous cell carcinoma (OSCC) is a significant health issue worldwide, and the expression of long non-coding RNAs (lncRNAs) are altered in these malignancies. The present study evaluated the expression level of ATXN1 CDC42EP1 genes and the lncRNAs related to these genes (lnc-ATXN1L, lnc-ATXN1, lnc-ATXN10, and lnc-CDC42EP1) in paraffin blocks of oral and pharyngeal squamous cell carcinoma (SCC) samples from patients referred to Amir Alam Hospital in Tehran, Iran.

METHODS AND RESULTS

This cross-sectional study was conducted on 76 paraffin blocks of oral and pharyngeal squamous cell carcinoma (SCC) samples from patients referred to Amir Alam Hospital in Tehran. The expression levels of ATXN1, CDC42EP1, lnc-ATXN1L, lnc-ATXN1, lnc-ATXN10, and lnc-CDC42EP1 were measured in all samples using a qPCR Master Mix kit. Real-time PCR was used to perform the reactions, and GAPDH was considered the housekeeping gene. Statistical analyses were conducted utilizing the Statistical Package for the Social Sciences (SPSS) version 22.0. The expression of lnc-ATXN1, lnc-ATXN10, and lnc-CDC42EP1 significantly differed between the two groups. All of them were downregulated (p < 0.05), and no significant difference was observed between the SCC samples and the adjacent tissue in other genes (p > 0.05). The expression of genes was not related to age, sex, size, and tumor location (p > 0.05).

CONCLUSIONS

Dysexpression of lnc-ATXN1, lnc-ATXN10, and lnc-CDC42EP1 can be used for diagnosing OSCC.

Passive Smoking-Induced Mutagenesis as a Promoter of Lung Carcinogenesis

INTRODUCTION

The International Agency for Research on Cancer has classified passive smoking (PS) or secondhand smoke exposure as a group 1 carcinogen linked to lung cancer. However, in contrast to active smoking, the mutagenic properties of PS remain unclear.

METHODS

A consecutive cohort of 564 lung adenocarcinoma samples from female never-smokers, who provided detailed information about their exposure to PS during adolescence and in their thirties through a questionnaire, was prepared. Of these, all 291 cases for whom frozen tumor tissues were available were subjected to whole exome sequencing to estimate tumor mutational burden, and the top 84 cases who were exposed daily, or not, to PS during adolescence, in their thirties or in both periods, were further subjected to whole genome sequencing.

RESULTS

A modest yet statistically significant increase in tumor mutational burden was observed in the group exposed to PS compared with the group not exposed to PS (median values = 1.44 versus 1.29 per megabase, respectively; p = 0.020). Instead of inducing driver oncogene mutations, PS-induced substantial subclonal mutations exhibiting APOBEC-type signatures, including SMAD4 and ADGRG6 hotspot mutations. A polymorphic APOBEC3A/3B allele-specific to the Asian population that leads to up-regulated expression of APOBEC3A accentuated the mutational load in individuals exposed daily to PS during adolescence.

CONCLUSIONS

This study reveals that PS-induced mutagenesis can promote lung carcinogenesis. The APOBEC3A/3B polymorphism may serve as a biomarker for identifying passive nonsmoking individuals at high risk of developing lung cancer.

In silico study of heterogeneous tumour-derived organoid response to CAR T-cell therapy

Chimeric antigen receptor (CAR) T-cell therapy is a promising immunotherapy for treating cancers. This method consists in modifying the patients' T-cells to directly target antigen-presenting cancer cells. One of the barriers to the development of this type of therapies, is target antigen heterogeneity. It is thought that intratumour heterogeneity is one of the leading determinants of therapeutic resistance and treatment failure. While understanding antigen heterogeneity is important for effective therapeutics, a good therapy strategy could enhance the therapy efficiency. In this work we introduce an agent-based model (ABM), built upon a previous ABM, to rationalise the outcomes of different CAR T-cells therapies strategies over heterogeneous tumour-derived organoids. We found that one dose of CAR T-cell therapy should be expected to reduce the tumour size as well as its growth rate, however it may not be enough to completely eliminate it. Moreover, the amount of free CAR T-cells (i.e. CAR T-cells that did not kill any cancer cell) increases as we increase the dosage, and so does the risk of side effects. We tested different strategies to enhance smaller dosages, such as enhancing the CAR T-cells long-term persistence and multiple dosing. For both approaches an appropriate dosimetry strategy is necessary to produce "effective yet safe" therapeutic results. Moreover, an interesting emergent phenomenon results from the simulations, namely the formation of a shield-like structure of cells with low antigen expression. This shield turns out to protect cells with high antigen expression. Finally we tested a multi-antigen recognition therapy to overcome antigen escape and heterogeneity. Our studies suggest that larger dosages can completely eliminate the organoid, however the multi-antigen recognition increases the risk of side effects. Therefore, an appropriate small dosages dosimetry strategy is necessary to improve the outcomes. Based on our results, it is clear that a proper therapeutic strategy could enhance the therapies outcomes. In that direction, our computational approach provides a framework to model treatment combinations in different scenarios and to explore the characteristics of successful and unsuccessful treatments.

Genome-wide loss of heterozygosity predicts aggressive, treatment-refractory behavior in pituitary neuroendocrine tumors

Pituitary neuroendocrine tumors (PitNETs) exhibiting aggressive, treatment-refractory behavior are the rare subset that progress after surgery, conventional medical therapies, and an initial course of radiation and are characterized by unrelenting growth and/or metastatic dissemination. Two groups of patients with PitNETs were sequenced: a prospective group of patients (n = 66) who consented to sequencing prior to surgery and a retrospective group (n = 26) comprised of aggressive/higher risk PitNETs. A higher mutational burden and fraction of loss of heterozygosity (LOH) was found in the aggressive, treatment-refractory PitNETs compared to the benign tumors (p = 1.3 × 10-10 and p = 8.5 × 10-9, respectively). Within the corticotroph lineage, a characteristic pattern of recurrent chromosomal LOH in 12 specific chromosomes was associated with treatment-refractoriness (occurring in 11 of 14 treatment-refractory versus 1 of 14 benign corticotroph PitNETs, p = 1.7 × 10-4). Across the cohort, a higher fraction of LOH was identified in tumors with TP53 mutations (p = 3.3 × 10-8). A machine learning approach identified loss of heterozygosity as the most predictive variable for aggressive, treatment-refractory behavior, outperforming the most common gene-level alteration, TP53, with an accuracy of 0.88 (95% CI: 0.70-0.96). Aggressive, treatment-refractory PitNETs are characterized by significant aneuploidy due to widespread chromosomal LOH, most prominently in the corticotroph tumors. This LOH predicts treatment-refractoriness with high accuracy and represents a novel biomarker for this poorly defined PitNET category.

Comprehensive pan-cancer analysis of inflammatory age-clock-related genes as prognostic and immunity markers based on multi-omics data

Inflammatory age (iAge) is a vital concept for understanding the intricate interplay between chronic inflammation and aging in the context of cancer. However, the importance of iAge-clock-related genes (iAge-CRGs) across cancers remains unexplored. This study aimed to explore the mechanisms and applications of these genes across diverse cancer types. We analyzed profiling data from over 10,000 individuals, covering 33 cancer types, 750 small molecule drugs, and 24 immune cell types. We focused on DCBLD2's function at the single-cell level and computed an iAge-CRG score using GSVA. This score was correlated with cancer pathways, immune infiltration, and survival. A signature was then derived using univariate Cox and LASSO regression, followed by ROC curve analysis, nomogram construction, decision curve analysis, and immunocytochemistry. Our comprehensive analysis revealed epigenetic, genomic, and immunogenomic alterations in iAge-CRGs, especially DCBLD2, leading to abnormal expression. Aberrant DCBLD2 expression strongly correlated with cancer-associated fibroblast infiltration and prognosis in multiple cancers. Based on GSVA results, we developed a risk model using five iAge-CRGs, which proved to be an independent prognostic index for uveal melanoma (UVM) patients. We also systematically evaluated the correlation between the iAge-related signature risk score and immune cell infiltration. iAge-CRGs, particularly DCBLD2, emerge as potential targets for enhancing immunotherapy outcomes. The strong correlation between abnormal DCBLD2 expression, cancer-associated fibroblast infiltration, and patient survival across various cancers underscores their significance. Our five-gene risk signature offers an independent prognostic tool for UVM patients, highlighting the crucial role of these genes in suppressing the immune response in UVM.Kindly check and confirm whether the corresponding affiliation is correctly identified.I identified the affiliation is correctly.thank you.Per style, a structured abstract is not allowed so we have changed the structured abstract to an unstructured abstract. Please check and confirm.I confirm the abstract is correctly ,thank you.

Mutational landscape of inflammatory breast cancer

BACKGROUND

Inflammatory breast cancer (IBC) is the most pro-metastatic form of BC. Better understanding of its enigmatic pathophysiology is crucial. We report here the largest whole-exome sequencing (WES) study of clinical IBC samples.

METHODS

We retrospectively applied WES to 54 untreated IBC primary tumor samples and matched normal DNA. The comparator samples were 102 stage-matched non-IBC samples from TCGA. We compared the somatic mutational profiles, spectra and signatures, copy number alterations (CNAs), HRD and heterogeneity scores, and frequencies of actionable genomic alterations (AGAs) between IBCs and non-IBCs. The comparisons were adjusted for the molecular subtypes.

RESULTS

The number of somatic mutations, TMB, and mutational spectra were not different between IBCs and non-IBCs, and no gene was differentially mutated or showed differential frequency of CNAs. Among the COSMIC signatures, only the age-related signature was more frequent in non-IBCs than in IBCs. We also identified in IBCs two new mutational signatures not associated with any environmental exposure, one of them having been previously related to HIF pathway activation. Overall, the HRD score was not different between both groups, but was higher in TN IBCs than TN non-IBCs. IBCs were less frequently classified as heterogeneous according to heterogeneity H-index than non-IBCs (21% vs 33%), and clonal mutations were more frequent and subclonal mutations less frequent in IBCs. More than 50% of patients with IBC harbored at least one high-level of evidence (LOE) AGA (OncoKB LOE 1-2, ESCAT LOE I-II), similarly to patients with non-IBC.

CONCLUSIONS

We provide the largest mutational landscape of IBC. Only a few subtle differences were identified with non-IBCs. The most clinically relevant one was the higher HRD score in TN IBCs than in TN non-IBCs, whereas the most intriguing one was the smaller intratumor heterogeneity of IBCs.

Genetic interactions reveal distinct biological and therapeutic implications in breast cancer

Co-occurrence and mutual exclusivity of genomic alterations may reflect the existence of genetic interactions, potentially shaping distinct biological phenotypes and impacting therapeutic response in breast cancer. However, our understanding of them remains limited. Herein, we investigate a large-scale multi-omics cohort (n = 873) and a real-world clinical sequencing cohort (n = 4,405) including several clinical trials with detailed treatment outcomes and perform functional validation in patient-derived organoids, tumor fragments, and in vivo models. Through this comprehensive approach, we construct a network comprising co-alterations and mutually exclusive events and characterize their therapeutic potential and underlying biological basis. Notably, we identify associations between TP53mut-AURKAamp and endocrine therapy resistance, germline BRCA1mut-MYCamp and improved sensitivity to PARP inhibitors, and TP53mut-MYBamp and immunotherapy resistance. Furthermore, we reveal that precision treatment strategies informed by co-alterations hold promise to improve patient outcomes. Our study highlights the significance of genetic interactions in guiding genome-informed treatment decisions beyond single driver alterations.

Beyond genetics: driving cancer with the tumour microenvironment behind the wheel

Cancer has long been viewed as a genetic disease of cumulative mutations. This notion is fuelled by studies showing that ageing tissues are often riddled with clones of complex oncogenic backgrounds coexisting in seeming harmony with their normal tissue counterparts. Equally puzzling, however, is how cancer cells harbouring high mutational burden contribute to normal, tumour-free mice when allowed to develop within the confines of healthy embryos. Conversely, recent evidence suggests that adult tissue cells expressing only one or a few oncogenes can, in some contexts, generate tumours exhibiting many of the features of a malignant, invasive cancer. These disparate observations are difficult to reconcile without invoking environmental cues triggering epigenetic changes that can either dampen or drive malignant transformation. In this Review, we focus on how certain oncogenes can launch a two-way dialogue of miscommunication between a stem cell and its environment that can rewire downstream events non-genetically and skew the morphogenetic course of the tissue. We review the cells and molecules of and the physical forces acting in the resulting tumour microenvironments that can profoundly affect the behaviours of transformed cells. Finally, we discuss possible explanations for the remarkable diversity in the relative importance of mutational burden versus tumour microenvironment and its clinical relevance.

Incongruity between T cell receptor recognition of breast cancer hotspot mutations ESR1 Y537S and D538G following exogenous peptide loading versus endogenous antigen processing

T cell receptor engineered T cell (TCR T) therapies have shown recent efficacy against certain types of solid metastatic cancers. However, to extend TCR T therapies to treat more patients across additional cancer types, new TCRs recognizing cancer-specific antigen targets are needed. Driver mutations in AKT1, ESR1, PIK3CA, and TP53 are common in patients with metastatic breast cancer (MBC) and if immunogenic could serve as ideal tumor-specific targets for TCR T therapy to treat this disease. Through IFN-γ ELISpot screening of in vitro expanded neopeptide-stimulated T cell lines from healthy donors and MBC patients, we identified reactivity towards 11 of 13 of the mutations. To identify neopeptide-specific TCRs, we then performed single-cell RNA sequencing of one of the T cell lines following neopeptide stimulation. Here, we identified an ESR1 Y537S specific T cell clone, clonotype 16, and an ESR1 Y537S/D538G dual-specific T cell clone, clonotype 21, which were HLA-B*40:02 and HLA-C*01:02 restricted, respectively. TCR Ts expressing these TCRs recognized and killed target cells pulsed with ESR1 neopeptides with minimal activity against ESR1 WT peptide. However, these TCRs failed to recognize target cells expressing endogenous mutant ESR1. To investigate the basis of this lack of recognition we performed immunopeptidomics analysis of a mutant-overexpressing lymphoblastoid cell line and found that the ESR1 Y537S neopeptide was not endogenously processed, despite binding to HLA-B*40:02 when exogenously pulsed onto the target cell. These results indicate that stimulation of T cells that likely derive from the naïve repertoire with pulsed minimal peptides may lead to the expansion of clones that recognize non-processed peptides, and highlights the importance of using methods that selectively expand T cells with specificity for antigens that are efficiently processed and presented.

Evolutionary trajectories of small cell lung cancer under therapy

The evolutionary processes that underlie the marked sensitivity of small cell lung cancer (SCLC) to chemotherapy and rapid relapse are unknown1-3. Here we determined tumour phylogenies at diagnosis and throughout chemotherapy and immunotherapy by multiregion sequencing of 160 tumours from 65 patients. Treatment-naive SCLC exhibited clonal homogeneity at distinct tumour sites, whereas first-line platinum-based chemotherapy led to a burst in genomic intratumour heterogeneity and spatial clonal diversity. We observed branched evolution and a shift to ancestral clones underlying tumour relapse. Effective radio- or immunotherapy induced a re-expansion of founder clones with acquired genomic damage from first-line chemotherapy. Whereas TP53 and RB1 alterations were exclusively part of the common ancestor, MYC family amplifications were frequently not constituents of the founder clone. At relapse, emerging subclonal mutations affected key genes associated with SCLC biology, and tumours harbouring clonal CREBBP/EP300 alterations underwent genome duplications. Gene-damaging TP53 alterations and co-alterations of TP53 missense mutations with TP73, CREBBP/EP300 or FMN2 were significantly associated with shorter disease relapse following chemotherapy. In summary, we uncover key processes of the genomic evolution of SCLC under therapy, identify the common ancestor as the source of clonal diversity at relapse and show central genomic patterns associated with sensitivity and resistance to chemotherapy.

Variant Allele Frequency Analysis of Circulating Tumor DNA as a Promising Tool in Assessing the Effectiveness of Treatment in Non-Small Cell Lung Carcinoma Patients

Despite the different possible paths of treatment, lung cancer remains one of the leading causes of death in oncological patients. New tools guiding the therapeutic process are under scientific investigation, and one of the promising indicators of the effectiveness of therapy in patients with NSCLC is variant allele frequency (VAF) analysis. VAF is a metric characterized as the measurement of the specific variant allele proportion within a genomic locus, and it can be determined using methods based on NGS or PCR. It can be assessed using not only tissue samples but also ctDNA (circulating tumor DNA) isolated from liquid biopsy. The non-invasive characteristic of liquid biopsy enables a more frequent collection of material and increases the potential of VAF analysis in monitoring therapy. Several studies have been performed on patients with NSCLC to evaluate the possibility of VAF usage. The research carried out so far demonstrates that the evaluation of VAF dynamics may be useful in monitoring tumor progression, remission, and recurrence during or after treatment. Moreover, the use of VAF analysis appears to be beneficial in making treatment decisions. However, several issues require better understanding and standardization before VAF testing can be implemented in clinical practice. In this review, we discuss the difficulties in the application of ctDNA VAF analysis in clinical routine, discussing the diagnostic and methodological challenges in VAF measurement in liquid biopsy. We highlight the possible applications of VAF-based measurements that are under consideration in clinical trials in the monitoring of personalized treatments for patients with NSCLC.

Comparison of genomic alterations in Epstein-Barr virus-positive and Epstein-Barr virus-negative diffuse large B-cell lymphoma

BACKGROUND

Epstein-Barr virus (EBV)-positive diffuse large B-cell lymphoma (EBV-posDLBCL) is an aggressive B-cell lymphoma that often presents similar morphological and immune phenotype features to that of EBV-negative DLBCL (EBV-negDLBCL).

AIMS AND METHODS

To better understand their difference in genomic landscape, we performed whole-exome sequencing (WES) of EBV-posDLBCL and EBV-negDLBCL.

RESULTS

This analysis revealed a new mutational signature 17 (unknown) and signature 29 (smoking) in EBV-posDLBCL as well as a specific mutational signature 24 (associated with aflatoxin) in EBV-negDLBCL. Compared with EBV-negDLBCL, more somatic copy number alterations (CNAs) and deletions were detected in EBV-posDLBCL (p = 0.01). The most frequent CNAs specifically detected in EBV-posDLBCL were gains at 9p24.1 (PDL1 and JAK2), 8q22.2-q24.23 (DEPTOR and MYC), and 7q31.31-q32.2 (MET), which were validated in additional EBV-posDLBCL cases. Overall, 53.7% (22/41) and 62.9% (22/35) of the cases expressed PD-L1 and c-MET, respectively, in neoplastic cells, whereas only 15.4% (4/26) expressed c-MYC. Neoplastic c-MET expression was positively correlated with PD-L1 (p < 0.001) and MYC expression (p = 0.016). However, EBV-posDLBCL cases did not show any differences in overall survival between PD-L1-, c-MET-, or c-MYC-positive and -negative cases or between age-related groups. Analysis of the association between somatic mutation load and EBV status showed no difference in the distribution of tumor mutant burden between the two lymphomas (p = 0.41). Recurrent mutations in EBV-posDLBCL implicated several genes, including DCAF8L1, KLF2, and NOL9, while in EBV-negDLBCL, ANK2, BPTF, and CNIH3 were more frequently mutated. Additionally, PIM1 is the most altered gene in all the WES-detected cases.

CONCLUSIONS

Our results confirm that genomic alteration differs significantly between EBV-posDLBCL and EBV-negDLBCL, and reveal new genetic alterations in EBV-posDLBCL. The positive correlation of c-MET and PD-L1/c-Myc expression may be involved in the pathogenesis of EBV-posDLBCL, which is should be explored prospectively in trials involving MET-directed therapies.

Circulating Tumor DNA Dynamics Fail to Predict Efficacy of Poly(ADP-ribose) Polymerase/VEGFR Inhibition in Patients With Heavily Pretreated Advanced Solid Tumors

PURPOSE

Cell-free circulating tumor DNA (ctDNA) has shown its potential as a quantitative biomarker for longitudinal monitoring of response to anticancer therapies. However, ctDNA dynamics have not been studied in patients with heavily pretreated, advanced solid tumors, for whom therapeutic responses can be weak. We investigated whether changes in ctDNA could predict clinical outcomes in such a cohort treated with combined poly(ADP-ribose) polymerase/vascular endothelial growth factor receptor inhibitor therapy.

MATERIALS AND METHODS

Patients with metastatic pancreatic ductal adenocarcinoma (PDAC), triple-negative breast cancer (TNBC), small-cell lung cancer (SCLC), or non-small-cell lung cancer (NSCLC) received up to 7 days of cediranib 30 mg orally once daily monotherapy lead-in followed by addition of olaparib 200 mg orally twice daily. Patients had progressed on a median of three previous lines of therapy. Plasma samples were collected before and after cediranib monotherapy lead-in and on combination therapy at 7 days, 28 days, and every 28 days thereafter. ctDNA was quantified from plasma samples using a multigene mutation-based assay. Radiographic assessment was performed every 8 weeks.

RESULTS

ctDNA measurements were evaluable in 63 patients. The median baseline ctDNA variant allele fractions (VAFs) were 20%, 28%, 27%, and 34% for PDAC, TNBC, SCLC, and NSCLC, respectively. No association was observed between baseline VAF and radiographic response, progression-free survival, or overall survival (OS). Similarly, no association was found between ctDNA decline and radiographic response or survival. However, an increase in ctDNA at 56 days of combination therapy was associated with disease progression and inferior OS in a landmark analysis.

CONCLUSION

ctDNA levels or dynamics did not correlate with radiographic response or survival outcomes in patients with advanced metastatic malignancies treated with olaparib and cediranib.

APOBEC Mutagenesis in Cancer Development and Susceptibility

APOBEC cytosine deaminases are prominent mutators in cancer, mediating mutations in over 50% of cancers. APOBEC mutagenesis has been linked to tumor heterogeneity, persistent cell evolution, and therapy responses. While emerging evidence supports the impact of APOBEC mutagenesis on cancer progression, the understanding of its contribution to cancer susceptibility and malignant transformation is limited. We examine the existing evidence for the role of APOBEC mutagenesis in carcinogenesis on the basis of the reported associations between germline polymorphisms in genes encoding APOBEC enzymes and cancer risk, insights into APOBEC activities from sequencing efforts of both malignant and non-malignant human tissues, and in vivo studies. We discuss key knowledge gaps and highlight possible ways to gain a deeper understanding of the contribution of APOBEC mutagenesis to cancer development.

CONIPHER: a computational framework for scalable phylogenetic reconstruction with error correction

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.

ZNF689 deficiency promotes intratumor heterogeneity and immunotherapy resistance in triple-negative breast cancer

Triple-negative breast cancer (TNBC) is an aggressive disease characterized by remarkable intratumor heterogeneity (ITH), which poses therapeutic challenges. However, the clinical relevance and key determinant of ITH in TNBC are poorly understood. Here, we comprehensively characterized ITH levels using multi-omics data across our center's cohort (n = 260), The Cancer Genome Atlas cohort (n = 134), and four immunotherapy-treated cohorts (n = 109). Our results revealed that high ITH was associated with poor patient survival and immunotherapy resistance. Importantly, we identified zinc finger protein 689 (ZNF689) deficiency as a crucial determinant of ITH formation. Mechanistically, the ZNF689-TRIM28 complex was found to directly bind to the promoter of long interspersed element-1 (LINE-1), inducing H3K9me3-mediated transcriptional silencing. ZNF689 deficiency reactivated LINE-1 retrotransposition to exacerbate genomic instability, which fostered ITH. Single-cell RNA sequencing, spatially resolved transcriptomics and flow cytometry analysis confirmed that ZNF689 deficiency-induced ITH inhibited antigen presentation and T-cell activation, conferring immunotherapy resistance. Pharmacological inhibition of LINE-1 significantly reduced ITH, enhanced antitumor immunity, and eventually sensitized ZNF689-deficient tumors to immunotherapy in vivo. Consistently, ZNF689 expression positively correlated with favorable prognosis and immunotherapy response in clinical samples. Altogether, our study uncovers a previously unrecognized mechanism underlying ZNF689 deficiency-induced ITH and suggests LINE-1 inhibition combined with immunotherapy as a novel treatment strategy for TNBC.

Oral tumor heterogeneity, its implications for patient monitoring and designing anti-cancer strategies

Oral cancer tumors occur in the mouth and are mainly derived from oral mucosa linings. It is one of the most common and fatal malignant diseases worldwide. The intratumor heterogeneity (ITH) of oral cancerous tumor is vast, so it is challenging to study and interpret. Due to environmental selection pressures, ITH arises through diverse genetic, epigenetic, and metabolic alterations. The ITH also talks about peri-tumoral vascular/ lymphatic growth, perineural permeation, tumor necrosis, invasion, and clonal expansion/ the coexistence of multiple subclones in a single tumor. The heterogeneity offers tumors the adaptability to survive, induce growth/ metastasis, and, most importantly, escape antitumor therapy. Unfortunately, the ITH is prioritized less in determining disease pathology than the traditional TNM classifications or tumor grade. Understanding ITH is challenging, but with the advancement of technology, this ITH can be decoded. Tumor genomics, proteomics, metabolomics, and other modern analyses can provide vast information. This information in clinics can assist in understanding a tumor's severity and be used for diagnostic, prognostic, and therapeutic decision-making. Lastly, the oral tumor ITH can lead to individualized, targeted therapy strategies fighting against OC.

Neoadjuvant pyrotinib plus trastuzumab and chemotherapy for HER2-positive breast cancer: a prospective cohort study

BACKGROUND

This prospective study aims to investigate the efficacy and safety of pyrotinib (P) combined with 4 cycles of epirubicin and cyclophosphamide followed by 4 cycles of taxane and trastuzumab (P + EC-TH) regimen as neoadjuvant therapy for human epidermal growth factor receptor 2 (HER2) positive breast cancer and to investigate the predictive value of p53, p63, and epidermal growth factor receptor (EGFR) status for neoadjuvant efficacy.

METHODS

A total of 138 HER2-positive breast cancer patients who received neoadjuvant therapy and underwent surgery were included. Case group: 55 patients received P + EC-TH regimen.

CONTROL GROUP

83 patients received EC-TH regimen. The chi-square test, Fisher's exact test, and logistic regression analysis were applied. The primary endpoint was total pathologic complete response (tpCR), and the secondary endpoints were breast pathologic complete response (bpCR), overall response rate (ORR), and adverse events (AEs).

RESULTS

In the case group, the tpCR rate was 63.64% (35/55), the bpCR rate was 69.09% (38/55), and the ORR was 100.00% (55/55). In the control group, the tpCR rate was 39.76% (33/83), the bpCR rate was 44.58% (37/83), and the ORR was 95.18% (79/83). The case group had significantly higher tpCR and bpCR rates than those of the control group (P < 0.05), but there was no significant difference in ORR (P > 0.05). The tpCR was associated with the status of estrogen receptor (ER), progesterone receptor (PR), and androgen receptor (AR), and the patients with any negative ER, PR, AR, or combined, were more likely to achieve tpCR than those with positive results (P < 0.05). The p53-positive patients were more likely to achieve tpCR and bpCR than p53-negative patients (P < 0.05). The incidence of hypokalemia and diarrhea in the case group was higher than that in the control group (P < 0.05). The AEs developed were all manageable, and no treatment-related death occurred.

CONCLUSION

The efficacy and safety of the P + EC-TH regimen were verified by this study. The HER2-positive breast cancer patients treated with the EC-TH neoadjuvant regimen were more likely to achieve tpCR or bpCR if pyrotinib was administered simultaneously.

Time-resolved, integrated analysis of clonally evolving genomes

Clonal genome evolution is a key feature of asexually reproducing species and human cancer development. While many studies have described the landscapes of clonal genome evolution in cancer, few determine the underlying evolutionary parameters from molecular data, and even fewer integrate theory with data. We derived theoretical results linking mutation rate, time, expansion dynamics, and biological/clinical parameters. Subsequently, we inferred time-resolved estimates of evolutionary parameters from mutation accumulation, mutational signatures and selection. We then applied this framework to predict the time of speciation of the marbled crayfish, an enigmatic, globally invasive parthenogenetic freshwater crayfish. The results predict that speciation occurred between 1986 and 1990, which is consistent with biological records. We also used our framework to analyze whole-genome sequencing datasets from primary and relapsed glioblastoma, an aggressive brain tumor. The results identified evolutionary subgroups and showed that tumor cell survival could be inferred from genomic data that was generated during the resection of the primary tumor. In conclusion, our framework allowed a time-resolved, integrated analysis of key parameters in clonally evolving genomes, and provided novel insights into the evolutionary age of marbled crayfish and the progression of glioblastoma.

Acute expression of human APOBEC3B in mice results in RNA editing and lethality

BACKGROUND

RNA editing has been described as promoting genetic heterogeneity, leading to the development of multiple disorders, including cancer. The cytosine deaminase APOBEC3B is implicated in tumor evolution through DNA mutation, but whether it also functions as an RNA editing enzyme has not been studied.

RESULTS

Here, we engineer a novel doxycycline-inducible mouse model of human APOBEC3B-overexpression to understand the impact of this enzyme in tissue homeostasis and address a potential role in C-to-U RNA editing. Elevated and sustained levels of APOBEC3B lead to rapid alteration of cellular fitness, major organ dysfunction, and ultimately lethality in mice. Importantly, RNA-sequencing of mouse tissues expressing high levels of APOBEC3B identifies frequent UCC-to-UUC RNA editing events that are not evident in the corresponding genomic DNA.

CONCLUSIONS

This work identifies, for the first time, a new deaminase-dependent function for APOBEC3B in RNA editing and presents a preclinical tool to help understand the emerging role of APOBEC3B as a driver of carcinogenesis.

Genomic profiling of subcutaneous patient-derived xenografts reveals immune constraints on tumor evolution in childhood solid cancer

Subcutaneous patient-derived xenografts (PDXs) are an important tool for childhood cancer research. Here, we describe a resource of 68 early passage PDXs established from 65 pediatric solid tumor patients. Through genomic profiling of paired PDXs and patient tumors (PTs), we observe low mutational similarity in about 30% of the PT/PDX pairs. Clonal analysis in these pairs show an aggressive PT minor subclone seeds the major clone in the PDX. We show evidence that this subclone is more immunogenic and is likely suppressed by immune responses in the PT. These results suggest interplay between intratumoral heterogeneity and antitumor immunity may underlie the genetic disparity between PTs and PDXs. We further show that PDXs generally recapitulate PTs in copy number and transcriptomic profiles. Finally, we report a gene fusion LRPAP1-PDGFRA. In summary, we report a childhood cancer PDX resource and our study highlights the role of immune constraints on tumor evolution.

Comprehensive genomic characterization of sporadic synchronous colorectal cancer: Implications for treatment optimization and clinical outcome

Sporadic synchronous colorectal cancer (SCRC) refers to multiple primary CRC tumors detected simultaneously in an individual without predisposing hereditary conditions, which accounts for the majority of multiple CRCs while lacking a profound understanding of the genomic landscape and evolutionary dynamics to optimize its treatment. In this study, 103 primary tumor samples from 51 patients with SCRC undergo whole-exome sequencing. The germline and somatic mutations and evolutionary and clinical features are comprehensively investigated. Somatic genetic events are largely inconsistent between paired tumors. Compared with solitary CRC, SCRCs have higher prevalence of tumor mutation burden high (TMB-H; 33.3%) and microsatellite-instability high (MSI-H; 29.4%) and different mutation frequencies in oncogenic signaling pathways. Moreover, neutrally evolving SCRC tumors are associated with higher intratumoral heterogeneity and better prognosis. These findings unveil special molecular features, carcinogenesis, and prognosis of sporadic SCRC. Strategies for targeted therapy and immunotherapy should be optimized accordingly.

The nuclear factor erythroid 2-related factor 2/p53 axis in breast cancer

One of the most important factors involved in the response to oxidative stress (OS) is the nuclear factor erythroid 2-related factor 2 (Nrf2), which regulates the expression of components such as antioxidative stress proteins and enzymes. Under normal conditions, Kelch-like ECH-associated protein 1 (Keap1) keeps Nrf2 in the cytoplasm, thus preventing its translocation to the nucleus and inhibiting its role. It has been established that Nrf2 has a dual function; on the one hand, it promotes angiogenesis and cancer cell metastasis while causing resistance to drugs and chemotherapy. On the other hand, Nrf2 increases expression and proliferation of glutathione to protect cells against OS. p53 is a tumour suppressor that activates the apoptosis pathway in aging and cancer cells in addition to stimulating the glutaminolysis and antioxidant pathways. Cancer cells use the antioxidant ability of p53 against OS. Therefore, in the present study, we discussed function of Nrf2 and p53 in breast cancer (BC) cells to elucidate their role in protection or destruction of cancer cells as well as their drug resistance or antioxidant properties.

Combined cytotoxic and immune-stimulatory gene therapy using Ad-TK and Ad-Flt3L: Translational developments from rodents to glioma patients

Gliomas are the most prevalent and devastating primary malignant brain tumors in adults. Despite substantial advances in understanding glioma biology, there have been no regulatory drug approvals in the US since bevacizumab in 2009 and tumor treating fields in 2011. Recent phase III clinical trials have failed to meet their prespecified therapeutic primary endpoints, highlighting the need for novel therapies. The poor prognosis of glioma patients, resistance to chemo-radiotherapy, and the immunosuppressive tumor microenvironment underscore the need for the development of novel therapies. Gene therapy-based immunotherapeutic strategies that couple the ability of the host immune system to specifically kill glioma cells and develop immunological memory have shown remarkable progress. Two adenoviral vectors expressing Ad-HSV1-TK/GCV and Ad-Flt3L have shown promising preclinical data, leading to FDA approval of a non-randomized, phase I open-label, first in human trial to test safety, cytotoxicity, and immune-stimulatory efficiency in high-grade glioma patients (NCT01811992). This review provides a thorough overview of immune-stimulatory gene therapy highlighting recent advancements, potential drawbacks, future directions, and recommendations for future implementation of clinical trials.

Precision medicine in the era of multi-omics: can the data tsunami guide rational treatment decision?

Precision medicine for cancer is rapidly moving to an approach that integrates multiple dimensions of the biology in order to model mechanisms of cancer progression in each patient. The discovery of multiple drivers per tumor challenges medical decision that faces several treatment options. Drug sensitivity depends on the actionability of the target, its clonal or subclonal origin and coexisting genomic alterations. Sequencing has revealed a large diversity of drivers emerging at treatment failure, which are potential targets for clinical trials or drug repurposing. To effectively prioritize therapies, it is essential to rank genomic alterations based on their proven actionability. Moving beyond primary drivers, the future of precision medicine necessitates acknowledging the intricate spatial and temporal heterogeneity inherent in cancer. The advent of abundant complex biological data will make artificial intelligence algorithms indispensable for thorough analysis. Here, we will discuss the advancements brought by the use of high-throughput genomics, the advantages and limitations of precision medicine studies and future perspectives in this field.

KRAS G12C in advanced NSCLC: Prevalence, co-mutations, and testing

KRAS is the most commonly mutated oncogene in advanced, non-squamous, non-small cell lung cancer (NSCLC) in Western countries. Of the various KRAS mutants, KRAS G12C is the most common variant (~40%), representing 10-13% of advanced non-squamous NSCLC. Recent regulatory approvals of the KRASG12C-selective inhibitors sotorasib and adagrasib for patients with advanced or metastatic NSCLC harboring KRASG12C have transformed KRAS into a druggable target. In this review, we explore the evolving role of KRAS from a prognostic to a predictive biomarker in advanced NSCLC, discussing KRAS G12C biology, real-world prevalence, clinical relevance of co-mutations, and approaches to molecular testing. Real-world evidence demonstrates significant geographic differences in KRAS G12C prevalence (8.9-19.5% in the US, 9.3-18.4% in Europe, 6.9-9.0% in Latin America, and 1.4-4.3% in Asia) in advanced NSCLC. Additionally, the body of clinical data pertaining to KRAS G12C co-mutations such as STK11, KEAP1, and TP53 is increasing. In real-world evidence, KRAS G12C-mutant NSCLC was associated with STK11, KEAP1, and TP53 co-mutations in 10.3-28.0%, 6.3-23.0%, and 17.8-50.0% of patients, respectively. Whilst sotorasib and adagrasib are currently approved for use in the second-line setting and beyond for patients with advanced/metastatic NSCLC, testing and reporting of the KRAS G12C variant should be included in routine biomarker testing prior to first-line therapy. KRAS G12C test results should be clearly documented in patients' health records for actionability at progression. Where available, next-generation sequencing is recommended to facilitate simultaneous testing of potentially actionable biomarkers in a single run to conserve tissue. Results from molecular testing should inform clinical decisions in treating patients with KRAS G12C-mutated advanced NSCLC.

Predictive modelling of response to neoadjuvant therapy in HER2+ breast cancer

HER2-positive (HER2+) breast cancer accounts for 20-25% of all breast cancers. Predictive biomarkers of neoadjuvant therapy response are needed to better identify patients with early stage disease who may benefit from tailored treatments in the adjuvant setting. As part of the TCHL phase-II clinical trial (ICORG10-05/NCT01485926) whole exome DNA sequencing was carried out on normal-tumour pairs collected from 22 patients. Here we report predictive modelling of neoadjuvant therapy response using clinicopathological and genomic features of pre-treatment tumour biopsies identified age, estrogen receptor (ER) status and level of immune cell infiltration may together be important for predicting response. Clonal evolution analysis of longitudinally collected tumour samples show subclonal diversity and dynamics are evident with potential therapy resistant subclones detected. The sources of greater pre-treatment immunogenicity associated with a pathological complete response is largely unexplored in HER2+ tumours. However, here we point to the possibility of APOBEC associated mutagenesis, specifically in the ER-neg/HER2+ subtype as a potential mediator of this immunogenic phenotype.

Proteomic Profiling and Tumor Microenvironment Characterization Reveal Molecular and Immunological Hallmarks of Left-Sided and Right-Sided Colon Cancer Tumorigenesis

Left-sided and right-sided colon cancer (LSCC and RSCC) display different biological and clinical characteristics. However, the differences in their tumorigenesis and tumor microenvironment remain unclear. In this study, we profiled the proteomic landscapes of LSCC and RSCC with data-independent acquisition mass spectrometry (DIA-MS) using fresh tumor and adjacent normal tissues from 24 patients. A total of 7403 proteingroups were primarily identified with DIA-MS. After quality control, 7212 proteingroups were used for further analysis. Through comparing the difference in proteomic profiles between LSCC and RSCC samples, 2556 commonly and 1982 region-type-specific regulated proteingroups were characterized. During the development of LSCC and RSCC, metabolic, growth, cell division, cell adhesion, and migration pathways were found to be significantly dysregulated (P < 0.05), which was further confirmed by transcriptome data from TCGA. Compared to RSCC, most parts of the immune-related signatures, immune cell infiltration scores, and overall immune scores of LSCC were higher. The systematic elucidation of proteomic and transcriptomic profiles in this work improves our understanding of tumorigenesis and immune microenvironment characteristics of LSCC and RSCC.

Genetic inhibitors of APOBEC3B-induced mutagenesis

The cytidine deaminases APOBEC3A (A3A) and APOBEC3B (A3B) are prominent mutators of human cancer genomes. However, tumor-specific genetic modulators of APOBEC-induced mutagenesis are poorly defined. Here, we used a screen to identify 61 gene deletions that increase A3B-induced mutations in yeast. We also determined whether each deletion was epistatic with Ung1 loss, which indicated whether the encoded factors participate in the homologous recombination (HR)-dependent bypass of A3B/Ung1-dependent abasic sites or suppress A3B-catalyzed deamination by protecting against aberrant formation of single-stranded DNA (ssDNA). We found that the mutation spectra of A3B-induced mutations revealed genotype-specific patterns of strand-specific ssDNA formation and nucleotide incorporation across APOBEC-induced lesions. Combining these three metrics, we were able to establish a multifactorial signature of APOBEC-induced mutations specific to (1) failure to remove H3K56 acetylation, (2) defective CTF18-RFC complex function, and (3) defective HR-mediated bypass of APOBEC-induced lesions. We extended these results by analyzing mutation data for human tumors and found BRCA1/2-deficient breast cancers display three- to fourfold more APOBEC-induced mutations. Mirroring our results in yeast, Rev1-mediated C-to-G substitutions are mainly responsible for increased APOBEC-signature mutations in BRCA1/2-deficient tumors, and these mutations associate with lagging strand synthesis during replication. These results identify important factors that influence DNA replication dynamics and likely the abundance of APOBEC-induced mutation during tumor progression. They also highlight a novel role for BRCA1/2 during HR-dependent lesion bypass of APOBEC-induced lesions during cancer cell replication.

Oncogene-like addiction to aneuploidy in human cancers

Most cancers exhibit aneuploidy, but its functional significance in tumor development is controversial. Here, we describe ReDACT (Restoring Disomy in Aneuploid cells using CRISPR Targeting), a set of chromosome engineering tools that allow us to eliminate specific aneuploidies from cancer genomes. Using ReDACT, we created a panel of isogenic cells that have or lack common aneuploidies, and we demonstrate that trisomy of chromosome 1q is required for malignant growth in cancers harboring this alteration. Mechanistically, gaining chromosome 1q increases the expression of MDM4 and suppresses p53 signaling, and we show that TP53 mutations are mutually exclusive with 1q aneuploidy in human cancers. Thus, tumor cells can be dependent on specific aneuploidies, raising the possibility that these "aneuploidy addictions" could be targeted as a therapeutic strategy.

Resolving the spatial architecture of myeloma and its microenvironment at the single-cell level

In multiple myeloma spatial differences in the subclonal architecture, molecular signatures and composition of the microenvironment remain poorly characterized. To address this shortcoming, we perform multi-region sequencing on paired random bone marrow and focal lesion samples from 17 newly diagnosed patients. Using single-cell RNA- and ATAC-seq we find a median of 6 tumor subclones per patient and unique subclones in focal lesions. Genetically identical subclones display different levels of spatial transcriptional plasticity, including nearly identical profiles and pronounced heterogeneity at different sites, which can include differential expression of immunotherapy targets, such as CD20 and CD38. Macrophages are significantly depleted in the microenvironment of focal lesions. We observe proportional changes in the T-cell repertoire but no site-specific expansion of T-cell clones in intramedullary lesions. In conclusion, our results demonstrate the relevance of considering spatial heterogeneity in multiple myeloma with potential implications for models of cell-cell interactions and disease progression.

Dysregulated Signalling Pathways Driving Anticancer Drug Resistance

One of the leading causes of death worldwide, in both men and women, is cancer. Despite the significant development in therapeutic strategies, the inevitable emergence of drug resistance limits the success and impedes the curative outcome. Intrinsic and acquired resistance are common mechanisms responsible for cancer relapse. Several factors crucially regulate tumourigenesis and resistance, including physical barriers, tumour microenvironment (TME), heterogeneity, genetic and epigenetic alterations, the immune system, tumour burden, growth kinetics and undruggable targets. Moreover, transforming growth factor-beta (TGF-β), Notch, epidermal growth factor receptor (EGFR), integrin-extracellular matrix (ECM), nuclear factor kappa-light-chain-enhancer of activated B cells (NF-κB), phosphoinositol-3-kinase/protein kinase B/mammalian target of rapamycin (PI3K/Akt/mTOR), wingless-related integration site (Wnt/β-catenin), Janus kinase/signal transducers and activators of transcription (JAK/STAT) and RAS/RAF/mitogen-activated protein kinase (MAPK) signalling pathways are some of the key players that have a pivotal role in drug resistance mechanisms. To guide future cancer treatments and improve results, a deeper comprehension of drug resistance pathways is necessary. This review covers both intrinsic and acquired resistance and gives a comprehensive overview of recent research on mechanisms that enable cancer cells to bypass barriers put up by treatments, and, like "satellite navigation", find alternative routes by which to carry on their "journey" to cancer progression.

The landscape of T cell antigens for cancer immunotherapy

The remarkable capacity of immunotherapies to induce durable regression in some patients with metastatic cancer relies heavily on T cell recognition of tumor-presented antigens. As checkpoint-blockade therapy has limited efficacy, tumor antigens have the potential to be exploited for complementary treatments, many of which are already in clinical trials. The surge of interest in this topic has led to the expansion of the tumor antigen landscape with the emergence of new antigen categories. Nonetheless, how different antigens compare in their ability to elicit efficient and safe clinical responses remains largely unknown. Here, we review known cancer peptide antigens, their attributes and the relevant clinical data and discuss future directions.

Clinical trial designs for evaluating and exploiting cancer evolution

The evolution of drug-resistant cell subpopulations causes cancer treatment failure. Current preclinical evidence shows that it is possible to model herding of clonal evolution and collateral sensitivity where an initial treatment could favourably influence the response to a subsequent one. Novel therapy strategies exploiting this understanding are being considered, and clinical trial designs for steering cancer evolution are needed. Furthermore, preclinical evidence suggests that different subsets of drug-sensitive and resistant clones could compete between themselves for nutrients/blood supply, and clones that populate a tumour do so at the expense of other clones. Treatment paradigms based on this clinical application of exploiting cell-cell competition include intermittent dosing regimens or cycling different treatments before progression. This will require clinical trial designs different from the conventional practice of evaluating responses to individual therapy regimens. Next-generation sequencing to assess clonal dynamics longitudinally will improve current radiological assessment of clinical response/resistance and be incorporated into trials exploiting evolution. Furthermore, if understood, clonal evolution can be used to therapeutic advantage, improving patient outcomes based on a new generation of clinical trials.