Mutational Landscape Research Articles

Genomic progression for local invasion of cutaneous squamous cell carcinoma from the superficial to the deep portion.

Cutaneous squamous cell carcinoma (cSCC) is the second most common skin cancer. While many treatments exist, our understanding of its genomic progression, especially from the epidermis to the deep dermis, remains limited. This study aims to identify genetic mutations associated with the progression of cSCC into the deep dermis, providing insights into its aggressive behavior and high-risk features. We performed high-depth whole-exome sequencing on 12 cSCC tissues, along with paired normal tissues from six patients, using microdissection techniques. The mutational analysis focused on identifying alterations enriched during cSCC progression. Gene Ontology enrichment analysis, immunohistochemical assays, and external single-cell RNA data were utilized for validation. A total of 8863 non-synonymous somatic mutations were identified in 4092 genes across the superficial and deep portions of cSCCs. Analysis of deep portion mutations revealed a significant correlation with gene ontology biological processes, particularly cell junction organization, and cell-cell adhesion. Clonal mutations in these processes were more prevalent in the deep portions, indicating their impact on the cSCC mutation landscape. Genetic evolution analysis identified 29 causal genes associated with dermal invasion in cSCC. We highlight somatic mutations in cSCC, revealing heterogeneity between superficial and deep regions. Altered genes in cell junction organization and cell-cell adhesion emerged as pivotal in dermal invasion. We identified 29 causal genes primarily in deep tumor regions. Our findings emphasize analyzing multiple tumor regions to capture varied mutational landscapes. These insights advance our understanding of cSCC progression, emphasizing genetic and cellular changes during tumor evolution.

Integrated single-cell and bulk RNA-Seq analysis enhances prognostic accuracy of PD-1/PD-L1 immunotherapy response in lung adenocarcinoma through necroptotic anoikis gene signatures

In addition to presenting significant diagnostic and treatment challenges, lung adenocarcinoma (LUAD) is the most common form of lung cancer. Using scRNA-Seq and bulk RNA-Seq data, we identify three genes referred to as HMR, FAM83A, and KRT6A these genes are related to necroptotic anoikis-related gene expression. Initial validation, conducted on the GSE50081 dataset, demonstrated the model's ability to categorize LUAD patients into high-risk and low-risk groups with significant survival differences. This model was further applied to predict responses to PD-1/PD-L1 blockade therapies, utilizing the IMvigor210 and GSE78220 cohorts, and showed strong correlation with patient outcomes, highlighting its potential in personalized immunotherapy. Further, LUAD cell lines were analyzed using quantitative PCR (qPCR) and Western blot analysis to confirm their expression levels, further corroborating the model's relevance in LUAD pathophysiology. The mutation landscape of these genes was also explored, revealing their broad implication in various cancer types through a pan-cancer analysis. The study also delved into molecular subclustering, revealing distinct expression profiles and associations with different survival outcomes, emphasizing the model’s utility in precision oncology. Moreover, the diversity of immune cell infiltration, analyzed in relation to the necroptotic anoikis signature, suggested significant implications for immune evasion mechanisms in LUAD. While the findings present a promising stride towards personalized LUAD treatment, especially in immunotherapy, limitations such as the retrospective nature of the datasets and the need for larger sample sizes are acknowledged. Prospective clinical trials and further experimental research are essential to validate these findings and enhance the clinical applicability of our prognostic model.

Biomimetic hydrogels with mesoscale collagen architecture for patient-derived tumor organoids culture

Patient-derived tumor organoids (PDTOs) shows great potential as a preclinical model. However, the current methods for establishing PDTOs primarily focus on modulating local properties, such as sub-micrometer topographies. Nevertheless, they neglect to capture the global millimeter or intermediate mesoscale architecture that have been demonstrated to influence tumor response to therapeutic treatment and tumor progression. In this study, we present a rapid technique for generating collagen bundles with an average length of 90 ± 27 μm and a mean diameter of 5 ± 1.5 μm from tumor tissue debris that underwent mechanical agitation following enzymatic digestion. The collagen bundles were subsequently utilized for the fabrication of biomimetic hydrogels, incorporating microbial transglutaminase (mTG) crosslinked gelatin. These biomimetic hydrogels, referred to as MC-gel, were specifically designed for patient-derived tumor organoids. The lung cancer organoids cultured in MC-gel exhibited larger diameters and higher cell viability compared to those cultured in gels lacking the mesoscale collagen bundle; moreover, their irregular morphology more closely resembled that observed in vivo. The MC-gel-based lung cancer organoids effectively replicated the histology and mutational landscapes observed in the original donor patient’s tumor tissue. Additionally, these lung cancer organoids showed a remarkable similarity in their gene expression and drug response across different matrices. This recently developed model holds great potential for investigating the occurrence, progression, metastasis, and management of tumors, thereby offering opportunities for personalized medicine and customized treatment options.

Comprehensive Characterization of the Integrin Family Across 32 Cancer Types

Abstract Integrin genes are widely involved in tumorigenesis. Yet, a comprehensive characterization of integrin family members and their interactome at the pan-cancer level is lacking. Here, we systematically analyzed integrin family in approximately 10,000 tumors across 32 cancer types. Globally, integrins represent a frequently altered and misexpressed pathway, with alteration and dysregulation overall being protumorigenic. Expression dysregulation, better than mutational landscape, of integrin family successfully identifies a subgroup of aggressive tumors with a high level of proliferation and stemness. The results showed that several molecular mechanisms collectively regulate integrin expression in a context-dependent manner. For potential clinical usage, we constructed a weighted scoring system, integrinScore, to measure integrin signaling patterns in individual tumors. Remarkably, integrinScore consistently correlated with predefined molecular subtypes in multiple cancers, with integrinScore high tumors being more aggressive. Importantly, integrinScore was cancer-dependently and closely associated with proliferation, stemness, tumor microenvironment, metastasis, and immune signatures. IntegrinScore also predicted patient’s response to immunotherapy. By mining drug databases, we unraveled an array of compounds that may modulate integrin signaling. Finally, we built a user-friendly database, Pan-cancer Integrin Explorer (PIExplorer) (http://computationalbiology.cn/PIExplorer), to facilitate researchers to explore integrin-related knowledge. Collectively, we provide a comprehensive characterization of integrins across cancers and offer gene-specific and cancer-specific rationales for developing integrin-targeted therapy.

Molecular alterations in claudin 18 suppressed and non-suppressed gastric adenocarcinomas to guide targeted therapies

ABSTRACT Background Gastric adenocarcinoma represents an aggressive type of cancer and an important cause of cancer mortality. Progress in gastric cancer therapeutics has resulted from a better understanding of the molecular pathogenesis of the disease and introduction of targeted therapies, but most gastric cancer patients still rely on non-targeted chemotherapy as the mainstay of treatment for advanced disease. Methods An analysis of publicly available series from The Cancer Genome Atlas (TCGA) gastric cancer cohort was undertaken to delineate the clinical and genomic landscape of gastric cancers with suppressed expression of claudin 18 compared with cancers with non-suppressed claudin 18. Claudin 18 suppressed cancers were defined as having an mRNA expression z-score relative to normal samples (log RNA Seq V2) of less than −1. Claudin 18 non-suppressed cancers were defined as having an mRNA expression z-score relative to normal samples (log RNA Seq V2) above 0.5. Results Gastric cancers with claudin 18 mRNA suppression represented 7.7% of the gastric adenocarcinomas of TCGA cohort, while non-suppressed cancers represented 46.6% of the cases. The two groups did not differ in clinical and genomic characteristics, such as mean age, histology, grade, and stage. The mutation landscape of claudin 18 suppressed cases included high mutation rates of TP53, of genes of the WNT/β-catenin pathway and of ubiquitin ligase FBXW7. Moreover, a subset of both claudin 18 suppressed and non-suppressed cancers displayed mutations in Mismatch Repair (MMR) associated genes or a high tumor mutation burden (TMB). At the mRNA expression level, claudin 18 suppressed gastric cancers showed up-regulation of EMT core transcription factor Snail 2 and down-regulation of genes of HLA cluster. The survival of gastric cancer patients with claudin 18 mRNA suppression was not significantly different compared with patients with non-suppressed claudin 18. Conclusion Sub-sets of gastric cancers with claudin 18 mRNA suppression displayed characteristics of potential therapeutic interest, such as mutations in WNT and PI3K pathways and MMR defects. These may guide the development of alternative targeted therapies, in this sub-set of gastric cancers which are not candidates for claudin 18 targeting therapies.

Unraveling the immunogenic cell death pathways in gastric adenocarcinoma: A multi-omics study.

Gastric cancer (GC) is a prevalent malignant tumor of the gastrointestinal (GI) system. However, the lack of reliable biomarkers has made its diagnosis, prognosis, and treatment challenging. Immunogenic cell death (ICD) is a type of programmed cell death that is strongly related to the immune system. However, its function in GC requires further investigation. We used multi-omics and multi-angle approaches to comprehensively explore the prognostic features of ICD in patients with stomach adenocarcinoma (STAD). At the single-cell level, we screened genes associated with ICD at the transcriptome level, selected prognostic genes related to ICD using weighted gene co-expression network analysis (WGCNA) and machine learning, and constructed a prognostic model. In addition, we constructed nomograms that incorporated pertinent clinical features and provided effective tools for prognostic prediction in clinical settings. We also investigated the sensitivity of the risk subgroups to both immunotherapy and drugs. Finally, in addition to quantitative real-time polymerase chain reaction, immunofluorescence was used to validate the expression of ICD-linked genes. Based on single-cell and transcriptome WGCNA analyses, we identified 34 ICD-related genes, of which 11 were related to prognosis. We established a prognostic model using the least absolute shrinkage and selection operator (LASSO) algorithm and identified dissimilarities in overall survival (OS) and progression-free survival (PFS) in risk subgroups. The nomograms associated with the ICD-related signature (ICDRS) demonstrated a good predictive value for clinical applications. Moreover, we detected changes in the tumor microenvironment (TME), including biological functions, mutation landscapes, and immune cell infiltration, between the high- and low-risk groups. We constructed an ICD-related prognostic model that incorporated features related to cell death. This model can serve as a useful tool for predicting the prognosis of GC, targeted prevention, and personalized medicine.

The Association between Mutational Signatures and Clinical Outcomes among Patients with Early-Onset Breast Cancer

Early-onset breast cancer (EoBC), defined by a diagnosis <40 years of age, is associated with poor prognosis. This study investigated the mutational landscape of non-metastatic EoBC and the prognostic relevance of mutational signatures using 100 tumour samples from Alberta, Canada. The MutationalPatterns package in R/Bioconductor was used to extract de novo single-base substitution (SBS) and insertion–deletion (indel) mutational signatures and to fit COSMIC SBS and indel signatures. We assessed associations between these signatures and clinical characteristics of disease, in addition to recurrence-free (RFS) and overall survival (OS). Five SBS and two indel signatures were extracted. The SBS13-like signature had higher relative contributions in the HER2-enriched subtype. Patients with higher than median contribution tended to have better RFS after adjustment for other prognostic factors (HR = 0.29; 95% CI: 0.08–1.06). An unsupervised clustering algorithm based on absolute contribution revealed three clusters of fitted COSMIC SBS signatures, but cluster membership was not associated with clinical variables or survival outcomes. The results of this exploratory study reveal various SBS and indel signatures may be associated with clinical features of disease and prognosis. Future studies with larger samples are required to better understand the mechanistic underpinnings of disease progression and treatment response in EoBC.

Mutation characteristics and molecular evolution of ovarian metastasis from gastric cancer and potential biomarkers for paclitaxel treatment

Ovarian metastasis is one of the major causes of treatment failure in patients with gastric cancer (GC). However, the genomic characteristics of ovarian metastasis in GC remain poorly understood. In this study, we enroll 74 GC patients with ovarian metastasis, with 64 having matched primary and metastatic samples. Here, we show a characterization of the mutation landscape of this disease, alongside an investigation into the molecular heterogeneity and pathway mutation enrichments between synchronous and metachronous metastasis. We classify patients into distinct clonal evolution patterns based on the distribution of mutations in paired samples. Notably, the parallel evolution group exhibits the most favorable prognosis. Additionally, by analyzing the differential response to chemotherapy, we identify potential biomarkers, including SALL4, CCDC105, and CLDN18, for predicting the efficacy of paclitaxel treatment. Furthermore, we validate that CLDN18 fusion mutations improve tumor response to paclitaxel treatment in GC with ovarian metastasis in vitro and vivo.

Abstract PS12-02: Exploration of ctDNA Dynamics in the PACE Trial: A Randomized Phase II Study of Fulvestrant, Palbociclib, and Avelumab for HR+/HER2- Metastatic Breast Cancer

Abstract Background The PACE trial for HR+/HER2- metastatic breast cancer (MBC) prospectively evaluated whether continuation of the CKD4/6 inhibitor (CDK4/6i) palbociclib (P) beyond progression on prior CDK4/6i and aromatase inhibitor (AI), with a change in endocrine therapy (ET) to fulvestrant (F), improved outcomes beyond change to F alone, as well as explored the activity of a palbociclib, fulvestrant, and avelumab (P+F+A) triplet. The primary analysis did not show improvement in progression-free survival (PFS) with F+P versus F alone. In this analysis, an evaluation of serial translational PACE samples was performed to gain insights into the mutational landscape of HR+/HER2- MBC after prolonged exposure to CDK4/6i, evaluate early mutational dynamics with ongoing therapy, and identify potential genetic markers predictive of treatment response. Additionally, the identification of acquired mutations at the time of tumor progression could provide valuable information about mechanisms of treatment resistance. Methods PACE is a multicenter randomized open-label investigator-initiated phase II trial. Eligible patients (pts) had HR+/HER2- evaluable MBC with prior progression on AI and any CDK4/6i after > 6 months (mo) of therapy in the MBC setting, or during/within 12 mo in the adjuvant setting, with no more than 1 prior line of chemotherapy for MBC. Pts were randomized 1:2:1 to F alone, F+P, or F+P+A, with tumor assessments every 8 weeks. Correlative samples including circulating tumor DNA (ctDNA) were collected at baseline (BL), times of tumor assessments (C3D1), and end of study/progression (EOT). Targeted NGS (Guardant360) to evaluate mutations (somatic single nucleotide variants and indels) was performed. Kaplan-Meier survival curves were used for assessing PFS and the log-rank test for comparisons. McNemar’s test with continuity correction was used to test the changes in the frequency of mutations. Results Of the 220 randomized pts, 211 contributed ctDNA samples (200 at BL, 129 at C3D1, 146 at EOT). Characteristics of the correlative science population were similar to the intention to treat population. At BL, the most common mutations observed were ESR1 (57%), TP53 (38%), PIK3CA (36%), GATA3 (20%) ATM (12%) and Rb1 (12%). The presence of a PIK3CA mutation vs WT-PIK3CA at BL was associated with shorter PFS in the F alone arm (p=0.001) but not in the F+P or F+P+A arms. ESR1 mutations present at baseline were D538G (34%), Y537S (21%), Y537N (15%), E380Q (10%) and L536H (4%). At C3D1, the most common mutations observed were ESR1 (47%), TP53 (33%), PIK3CA (26.1), GATA3 (18.1), ATM (13%) and PTEN (8.6%). Including only BL and C3D1 matched samples (N=124), the prevalence of ESR1 mutations decreased after 2 cycles of treatment (57% at BL vs 44% at C3D1), predominantly observed in the F+P+A arm. For those in any arm with at least 6 months of disease stability on study (N=53), the presence of an ESR1 D538G mutation decreased during therapy, with 38% having a mutation at BL, 24% at C3D1, and no pts gaining ESR1 D538G at C3D1. The prevalence of PIK3CA mutations also decreased on treatment (32% at BL to 24% at C3D1), with the greatest decrease in the F+P arm. A predictive model correlating early genomic dynamics and associations with outcome will be presented. Conclusions The PACE trial offers one of the largest clinical trial cohorts describing the genomic landscape of HR+/HER2- MBC post-CDK4/6 exposure. Results from this correlative study provide a comprehensive mutational landscape post-CDK4/6i, demonstrate associations between mutations and clinical outcomes, and suggest the potential value of early ctDNA dynamics in this setting. Citation Format: Rinath Jeselsohn, Jingxin Fu, Yue Ren, Reshma Mahtani, Cynthia Ma, Angela DeMichele, Massimo Cristofanilli, Jane Meisel, Kathy Miller, Yara Abdou, Elizabeth Riley, Rubina Qamar, Priyanka Sharma, Sonya Reid, Naomi Ko, Yuan Liu, Harold Burstein, Michelle DeMeo, Sara Tolaney, Erica Mayer. Exploration of ctDNA Dynamics in the PACE Trial: A Randomized Phase II Study of Fulvestrant, Palbociclib, and Avelumab for HR+/HER2- Metastatic Breast Cancer [abstract]. In: Proceedings of the 2023 San Antonio Breast Cancer Symposium; 2023 Dec 5-9; San Antonio, TX. Philadelphia (PA): AACR; Cancer Res 2024;84(9 Suppl):Abstract nr PS12-02.

Abstract PO5-27-03: Mutational impact of APOBEC3A and APOBEC3B in a human cell line

Abstract A major source of mutation in cancer is DNA cytosine deamination by APOBEC3 enzymes resulting in C-to-T and C-to-G mutations in TCA and TCT motifs. Here, we develop a selectable system to quantify genomic mutations and compare the mutagenic activities of three leading APOBEC3 candidates - APOBEC3A, APOBEC3B, and APOBEC3H. The human cell line, HAP1, is engineered to express the thymidine kinase (TK) gene of HSV-1, which confers sensitivity to ganciclovir. Expression of APOBEC3A and APOBEC3B, but not catalytic mutant controls or APOBEC3H, triggers increased frequencies of TK mutation and nearly indistinguishable TC-biased cytosine mutation profiles in the selectable TK reporter gene. Whole genome sequences from TK mutant clones enabled an analysis of thousands of single base substitution mutations and extraction of local sequence preferences with APOBEC3A preferring YTCW motifs over 70% of the time and APOBEC3B just under 50% of the time (Y=C/T; W=A/T). Signature comparisons with breast tumor whole genome sequences indicate that most malignancies manifest intermediate percentages of APOBEC3 signature mutations in YTCW motifs, mostly between 50 and 70%, suggesting that both enzymes are contributing in a combinatorial manner to the overall mutation landscape. These studies combine to help resolve a long-standing etiologic debate on the source of APOBEC3 signature mutations in cancer and indicate that future diagnostic and therapeutic efforts should focus on both enzymes. Citation Format: Michael Carpenter, Nuri Temiz, Mahmoud Ibrahim, Matthew Jarvis, Margaret Brown, Prokopios Argyris, William Brown, Douglas Yee, Reuben Harris. Mutational impact of APOBEC3A and APOBEC3B in a human cell line [abstract]. In: Proceedings of the 2023 San Antonio Breast Cancer Symposium; 2023 Dec 5-9; San Antonio, TX. Philadelphia (PA): AACR; Cancer Res 2024;84(9 Suppl):Abstract nr PO5-27-03.

Abstract PO1-15-11: Molecular and clinical disparities in breast cancer among East Asian populations

Abstract Background Breast Cancer (BC) is a complex disease with profound genomic aberrations. However, the underlying molecular disparity influenced by age and ethnicity remains elusive. Method We aimed to investigate the molecular properties of 843 primary and metastatic BC patients enrolled in the K-MASTER program. We explored the unique molecular profiles of breast cancer patients by categorizing them into two distinct age subgroups. The first subgroup comprised young breast cancer (YBC) patients, defined as women under the age of 40, while the second subgroup consisted of older breast cancer (OBC) patients. Our study aimed to characterize the molecular differences between these age-defined groups, providing insights into the underlying mechanisms of breast cancer in different age cohorts. Additionally, we leveraged large-scale genomic data from the TCGA and MSK-IMPACT studies to examine the ethnic-driven molecular and clinical disparities. Results The K-MASTER patients were mainly comprised of triple-negative breast cancer (TNBC) and HER2-positive tumors, while the TCGA and MSK-IMPACT cohorts exhibited a predominance of hormone receptor-positive (HR+) subtype. We observed a high prevalence of PI3KCA mutations in K-MASTER HER2+ tumors, particularly in older patients. Moreover, we identified increased mutation rates in DNA damage response molecules, including ARID1A, MSH6, and MLH1. Importantly, GATA3 mutations were less frequently observed in East Asian patients, which correlated with poor clinical outcomes. Our analysis also revealed specific associations between genomic aberrations, organotropisms, and metastatic breast cancer. In addition to characterizing the molecular disparities, we developed a gradient-boosting multivariable model to identify a new molecular signature that could predict the therapeutic response to platinum-based chemotherapy. Conclusion Our findings collectively provide unprecedented insights into the significance of age and ethnicity on the molecular and clinical characteristics of BC patients. Figure 1. Genomic landscape of KM BRCA. (A) Overall characteristic of KM BRCA sample and clinical data. Clinical features of KM BRCA sample. (B) Genomic landscape of somatic mutations and copy number alterations of BRCA by age group. Figure 2. Genomic landscape of KM BRCA: Significantly mutated genes by age group in KM BRCA. Figure 3. Survival Curve Utilizing a curated Gene Set Developed by a Gradient-Boosting Multivariable Model (A) OS of platinum-based therapy treated KM with curated Gene set. (B) OS of platinum-based therapy treated MSK with curated Gene set. Citation Format: Kim Ju Won, LEE JIWON, Kyong Hwa Park, Seung Pil Jung, Jason Sa, SHIN SANG WON, An Jung Seok. Molecular and clinical disparities in breast cancer among East Asian populations [abstract]. In: Proceedings of the 2023 San Antonio Breast Cancer Symposium; 2023 Dec 5-9; San Antonio, TX. Philadelphia (PA): AACR; Cancer Res 2024;84(9 Suppl):Abstract nr PO1-15-11.

Abstract PO5-14-05: Circulating tumor DNA fraction correlates with residual cancer burden post-neoadjuvant chemotherapy in triple negative breast cancer patients

Abstract Background: Pathologic response after preoperative/neoadjuvant chemotherapy is a continuum that can range from no residual cancer to extensive disease. Viable cancers can shed tumor DNA into the blood and we hypothesized that plasma ctDNA levels can reflect pathologic response as continuum. Patients and methods: Post-neoadjuvant plasma samples were obtained after completion of chemotherapy but before surgery from patients with stage I to III TNBC who received treatment with durvalumab concurrent with weekly nab-paclitaxel followed by dose dense doxorubicin/cyclophosphamide in a single arm neoadjuvant clinical trial (NCT02489448). Up to 3 mL of blood was collected into EDTA- and heparin- containing vacutainer tubes, and plasma cfDNA was extracted and twice purified in the lab at Predicine Inc. Pathologic response was quantified on a continuous scale using the Residual Cancer Burden (RCB) score, and also categorized as RCB-0 (complete pathologic response, n=21), -I (minimal residual cancer, n=7), -II (moderate residual cancer, n=13), and -III (extensive residual cancer, n=3). Median follow-up was 35 months (range 25-70 months). Tumor whole exome sequencing (WES) data was previously generated using HiSeq 4000 at the Yale Center for Genome Analysis (dbGaP: PRJNA558949). A personalized MRD panel of up to 60 personalized variants per patient was designed based on the tumor WES data, and ultra-deep next-generation sequencing (100,000X) of plasma cfDNA was performed using the PredicineBEACON assay for MRD detection by Predicine. CtDNA positivity (MRD) was defined as a weighted score equivalent to 2 or more tumor variants detected in the plasma. cfDNA tumor fraction was measured based on mutant allele fraction of observed variants. Adjusted tumor fraction could be confidently measured for 37 out of 44 patients that had enough confidently called personalized ctDNA variants. We compared tumor fraction versus RCB score as continuous variables, and also compared tumor fraction between RCB categories as ordinal variables. Invasive disease-free survival was examined by Kaplan-Meier analysis. Bioinformatic data analysis was performed using Graphpad Prism and R. Results: 44 of 48 plasma samples yielded successful cfDNA generation and sequencing. Six patients (13.6%) were MRD positive. Among these 6 cases, 3 had RCB-0 response (14.3% of RCB 0 cases), 1 had RCB II response (7.7% of RCB II cases), and 2 had RCB III response (66.7% of RCB III cases). The median cfDNA tumor fraction of patients with pathological complete response to treatment (pCR) was significantly lower than patients with residual disease (RD) (0.06% vs. 0.3%, p=0.02). There was a positive correlation between RCB scores and tumor fraction (r= 0.45, Spearman). Tumor fraction was more strongly correlated with later recurrence than a binary ctDNA positivity call. Tumor fraction was significantly higher in patients that later experienced local or metastatic recurrence (n=8) compared to those with no evidence of disease (n=29) (Wilcoxon test, p=0.02). Both ctDNA positive RCB III cases had metastatic recurrence while the remaining 4 ctDNA positive cases currently remain recurrence-free. Mutational genomic landscape of both tumor WES and plasma ctDNA samples identified hotspot mutations in 19 cancer-associated genes. Novel hot-spot mutations were observed in the plasma of two patients following treatment. TP53 mutations were highly prevalent in tissue and plasma, however no significant association of mutation of any gene with treatment response was detected. Conclusions: PredicineBeacon MRD assay was successfully performed on 3 mL plasma samples collected using from EDTA and heparin-containing tubes. Tumor fraction correlated with residual disease, RCB score and disease recurrence. Citation Format: Naing Lin Shan, Billie Gould, Kim Blenman, Julia Foldi, Pan Du, Frank Zhang, Myles Walsh, Lajos Pusztai. Circulating tumor DNA fraction correlates with residual cancer burden post-neoadjuvant chemotherapy in triple negative breast cancer patients [abstract]. In: Proceedings of the 2023 San Antonio Breast Cancer Symposium; 2023 Dec 5-9; San Antonio, TX. Philadelphia (PA): AACR; Cancer Res 2024;84(9 Suppl):Abstract nr PO5-14-05.

Abstract PO1-13-02: Circulating tumor DNA mutation landscape in HR+/HER2− patients with mBC treated with cyclin-dependent kinase 4/6 inhibitors in the SCRUM-Japan MONSTAR-SCREEN study

Abstract Introduction: Combination of cyclin-dependent kinase 4 and 6 inhibitors (CDK 4/6i) + endocrine therapy (ET) is recommended as first- or second-line treatment (1L or 2L Tx) for hormone receptor-positive, human epidermal growth factor receptor-2-negative metastatic breast cancer (HR+/HER2− mBC). In this study, we assessed cancer-specific genetic alterations in circulating tumor DNA (ctDNA) and tissues in abemaciclib-treated patients with HR+/HER2− mBC in Japan. Methods: MONSTAR-SCREEN is a multicenter study and part of a nation-wide cancer genome screening project for patients with advanced solid tumors (SCRUM Japan). Regardless of the Tx line, the current analysis included abemaciclib-treated patients with HR+/HER2− mBC (N=97). Data were collected between Jan 2020 and Dec 2021. Blood samples collected at ≤21 days before abemaciclib initiation (baseline; n=77) and at disease progression or Tx discontinuation (paired post-Tx; n=33) were tested with the FoundationOne® Liquid Companion Diagnostic panel. Additionally, archival or fresh tissue biopsy samples (n=79) were tested with the FoundationOne® Companion Diagnostic panel. Patient characteristics, baseline genetic alterations, neoplastic burden (measured by shedding rate and maximum variant allele frequency [VAF]), and alterations at clinical progression were reported. Results: The total population (N=97) included patients with next generation sequencing ctDNA and/or tissue data. All were female with median age of 57 years (interquartile range 50, 67). Bone, either alone or with other sites, was the most common metastatic site (61%). 78% of patients received abemaciclib in the first or second line of Tx (LoT). Fulvestrant with or without gonadotropin-releasing hormone (55%), letrozole (24%), and anastrozole (11%) were the most common concomitant drugs. Among patients with baseline ctDNA data (n=77), 30% and 21% had received prior ET and prior chemotherapy, respectively. In baseline ctDNA samples, PIK3CA (37%) and TP53 (28%) alterations were detected most frequently across all LoT; while ESR1 (16%), GATA3 (11%), and FGF3/4/19 (9%) were more frequently detected in samples from later LoT. The frequency of ESR1 alterations (p< 0.01), shedding rate, and maximum VAF (both p< 0.05) were significantly higher in prior ET recipients than those who did not receive prior ET. Frequency of post-Tx ESR1 alterations numerically decreased from baseline in patients on fulvestrant (23% to 18%) and numerically increased in patients on aromatase inhibitors (27% to 45%). Other genes showed either a numeric increase or similar alteration frequency between the baseline and post-tx. Newly acquired post-Tx alterations were detected in FGF3/4/19 (18%); PIK3CA, TP53, and RB1 (all 15%) and ESR1 (12%). The table shows PIK3CA and ESR1 alterations detected at baseline and post-Tx. In the tissue samples (n=79), TP53 (32% vs 20%) and ESR1 (14% vs 7%) alterations were numerically more frequent in metastatic versus primary lesions. Other gene alteration frequencies were similar. Conclusion: PIK3CA alterations were the most frequently detected genetic alterations in both pre- and post- abemaciclib Tx patients with HR+/HER2− mBC. Patients with prior ET had more frequent ESR1 alterations and higher neoplastic burden. Evidence from this study provides insight into the ctDNA dynamics and potential resistance mechanisms in patients with HR+/HER2− mBC. PIK3CA and ESR1 alterations detected in circulating tumor DNA at baseline and post abemaciclib treatment Citation Format: Masaya Hattori, Victoria Serelli-Lee, Yoichi Naito, Takashi Yamanaka, Hiroyuki Yasojima, Rikiya Nakamura, Takao Fujisawa, Mitsuho Imai, Yoshiaki Nakamura, Hideaki Bando, Tsutomu Kawaguchi, Takayuki Yoshino, Hiroji Iwata. Circulating tumor DNA mutation landscape in HR+/HER2− patients with mBC treated with cyclin-dependent kinase 4/6 inhibitors in the SCRUM-Japan MONSTAR-SCREEN study [abstract]. In: Proceedings of the 2023 San Antonio Breast Cancer Symposium; 2023 Dec 5-9; San Antonio, TX. Philadelphia (PA): AACR; Cancer Res 2024;84(9 Suppl):Abstract nr PO1-13-02.

Abstract PO3-23-11: Multiomics analysis of ER+ breast cancer with matched primaries, lymph node metastases and recurrent cancers on adjuvant endocrine therapy

Abstract Background: 80% of all breast cancers (BCs) are ER-positive (ER+). A significant number develop endocrine therapy resistance (ETR) . Clones develop in lymph nodes (LNs) which differ from the primary and these are implicated in ETR. A multiomics analysis of primary ER+ breast cancers, matched LN metastasis (LN+), if present, and matched recurrences (primary, nodal or distant) has been performed. Patients and Methods: The patients included are: Cohort A: 95 patients with matched primary and recurrent ER+ BC on different endocrine therapies. 35% had involved LNs at diagnosis. Recurrences: local in 59/95, concurrent local and nodal in 33/95, LN-only in 3/95. Median time to recurrence 4.5 years. Cohort B: 18 patients with ER+ LN+ BC who developed multiple recurrences (local and/or nodal, >3 consecutive recurrences/patient) on different lines of sequential endocrine therapy. Cohort C: 7 patients with ER+ breast cancer with multiple involved LNs treated with 3-6 months of neoadjuvant ET before surgery. Biopsies of primary and matched LNs performed at diagnosis and throughout treatment (3-4 timepoints/patient). Clinical response by serial ultrasound. Primary and LN samples (mean 12 samples/patient) were analysed together with normal DNA. Targeted RNAseq and DNA exome sequencing and whole-genome expression analysis was performed. Somatic mutations and copy number alterations (CNA) were determined, and differential gene expression analysis was performed. Validation of pathways implicated in ETR was by spatial NanoString GeoMx. All patients have long term follow-up. Results: Cohort A: A mean of 6 potential driver mutations identified by DNA sequencing in each primary BC (range 1-16). Most commonly mutated genes were PIK3CA, TP53, ERBB4, MAP3K1 and ERBB3. Multiple aberrations and a highly diverse mutational landscape were observed in all recurrences. In the majority of recurrences significant changes in variant allele frequency (VAF) were identified across several ETR-implicated genes and in over 50% new gains/losses were identified. Most frequently, losses were seen in ERBB3, ARID1A, HRAS and FOXA1 with gains in CDH1, NF1, PIK3CA, ESR1, TBX3 and GATA3 in recurrences vs primaries. Transcriptomic analysis revealed differentially expressed pathways including ER, HER2, GATA3, AKT, RAS and p63 signalling. Integrated analysis revealed several driver mutations with corresponding changes in upstream pathway activity. 33 patients had concurrent local and nodal recurrences and 52% had similar somatic mutation profiles in both with only minor variance in VAF. In 48%, 1 or more new potential driver mutations was identified in axillary nodes compared to recurrent primary cancers with losses in AKT1, CDH1 and TP53 being most common. Multiomics profiling of Cohorts B and C is complete and analysis underway. Discussion: This multiomics study is the largest cohort to-date of matched early and recurrent ER+/HER2- cancers. It sheds new light on the complex somatic and transcriptomic changes in local and nodal recurrent disease. The mechanisms of ETR are diverse and underlying mechanisms and clinically meaningful biomarkers including potentially actionable mutations and targets have been identified. Cohort C will uniquely capture and compare clonal evolution of breast cancer in primary and multiple matched lymph during neoadjuvant ET. Citation Format: Carlos Martinez-Perez, Charlene Kay, Rebecca Swan, Lorna Renshaw, J Michael Dixon, Arran K Turnbull. Multiomics analysis of ER+ breast cancer with matched primaries, lymph node metastases and recurrent cancers on adjuvant endocrine therapy [abstract]. In: Proceedings of the 2023 San Antonio Breast Cancer Symposium; 2023 Dec 5-9; San Antonio, TX. Philadelphia (PA): AACR; Cancer Res 2024;84(9 Suppl):Abstract nr PO3-23-11.

Abstract PO1-05-07: Circulating tumor DNA dynamics in acelERA Breast Cancer: a Phase II study of giredestrant for estrogen receptor-positive, HER2-negative, previously treated advanced breast cancer

Abstract BACKGROUND Giredestrant (GIR) is a highly potent, oral, selective estrogen receptor antagonist and degrader (SERD) that exhibits robust estrogen receptor (ER) occupancy. The Phase II, randomized, open-label acelERA Breast Cancer (BC) study (NCT04576455) assessed GIR vs physician’s choice of endocrine therapy (PCET) in second- or third-line ER-positive, HER2-negative advanced BC (ER+, HER2– aBC). The study did not reach statistical significance for the primary endpoint of investigator-assessed progression-free survival (PFS); however, the benefit of GIR was of larger magnitude among patients (pts) with ESR1-mutated tumors (ESR1m; a common cause of acquired resistance to endocrine therapy [ET]). We present an exploratory biomarker analysis of circulating tumor (ct)DNA dynamics. METHODS Pts (n = 303) were randomized 1:1 to GIR or PCET (75% of pts had fulvestrant [FUL]; 25%, an aromatase inhibitor). Partial response (PR), stable disease (SD), and progressive disease (PD) were categorized by RECIST v1.1. Plasma samples taken at Cycle 1, Day 1 (C1D1; n = 229), while on treatment (tx) at C2D1 (n = 220), and at the end of therapy (EOT; n = 155) were evaluated with the FoundationOne Liquid CDx next-generation sequencing assay. Gene mutations were defined as variants with known or likely impact on protein function. Composite tumor fraction (cTF) was defined as the total estimated tumor ctDNA content in each sample, and was only evaluable for a subset of samples (137 C1D1; 136 C2D1; 99 EOT). cTF or ESR1 mutant allele frequency (MAF) changes were calculated as a percent change from C1D1. ESR1m clonality was estimated as variant tumor fraction/estimated tumor fraction. Statistical significance was evaluated using the Mann–Whitney test. RESULTS The overall mutation landscape remained relatively unchanged with GIR tx at C2D1 or EOT, except for ESR1m prevalence, which declined from 43% at C1D1 to 26% on tx and 27% at EOT. With PCET, ESR1m prevalence decreased from 34% at C1D1 to 29% on tx, but increased to 43% by EOT. There were no concordant changes for other top mutated genes in this study, including PIK3CA, TP53, or DNMT3A. cTF decreased on tx in 38/65 (59%) and 32/59 (54%) of evaluable pts with GIR and PCET, respectively; the median change from C1D1 was –23% with GIR compared with –5% with PCET or –15% with FUL. The median change in cTF levels on tx with GIR was significantly greater in pts with a PR (–88%) vs PD (+6%, p = 0.002) or SD (–19%, p = 0.004). The degree of cTF decline on tx with GIR was significantly higher in pts with baseline ESR1m (median change –58%) vs no ESR1m detected (+5%, p = 0.012). Pts with a clonal ESR1m had higher levels of cTF decline on tx vs pts with subclonal ESR1m, suggesting cTF dynamics are a function of tumor heterogeneity. ESR1 MAF decreased on tx in 93% of pts with baseline ESR1m on GIR, vs 60% with PCET and 70% with FUL. ESR1 MAF decline on tx was significantly higher with GIR vs PCET at both C2D1 and EOT (p < 0.0001). The average ESR1 MAF decline on tx with GIR was greater in pts with a PR (–97%) vs PD (–54%, p = 0.025) or SD (–48%, p = 0.051), and all seven GIR-treated pts with a PR showed a near or total loss of ESR1 MAF by C2D1. The ESR1m variants D538G and Y537X had a significantly greater decline in MAF on tx with GIR vs PCET or FUL (p < 0.0001). CONCLUSIONS Data show that cTF and ESR1m ctDNA dynamics were associated with clinical response to GIR in ER+, HER2– aBC. ESR1 MAF decline was significantly greater with GIR vs PCET or FUL. ESR1 MAF declined to a greater degree with GIR compared with cTF, which is consistent with the larger magnitude of PFS benefit seen with GIR in pts with ESR1m tumors. The specific ESR1 variants D538G and Y537X showed greater sensitivity to GIR compared with PCET or FUL. Citation Format: Ann Collier, Aditya Bardia, Joo Hyuk Sohn, Elgene Lim, Marianna Chavez, Miguel Martín, Jorge Martinalbo, Pablo Perez-Moreno, Heather Moore. Circulating tumor DNA dynamics in acelERA Breast Cancer: a Phase II study of giredestrant for estrogen receptor-positive, HER2-negative, previously treated advanced breast cancer [abstract]. In: Proceedings of the 2023 San Antonio Breast Cancer Symposium; 2023 Dec 5-9; San Antonio, TX. Philadelphia (PA): AACR; Cancer Res 2024;84(9 Suppl):Abstract nr PO1-05-07.

Abstract PO1-13-05: Genomic characteristics and its therapeutic implications in breast cancer patients with detectable molecular residual disease

Abstract Background: Breast cancer constitutes approximately 30% of cancers in women, with a mortality-to-incidence ratio of 15%. While early and middle-stage breast cancer patients can undergo radical surgical resection, postoperative recurrence remains a significant challenge for clinicians. Molecular residual disease (MRD) is the main cause of postoperative recurrence of breast cancer. However, the baseline tumor genomic characteristics and therapeutic implications of breast cancer patients with detectable MRD are still unknown. Methods: In this study, we enrolled 80 patients with breast cancer who underwent next-generation sequencing (NGS)-based genetic testing of 1,021 cancer-related genes performed on baseline tumor and postoperative plasma from Army Specialty Medical Center between June 2017 and January 2023. The criteria for the patient’s enrollment were: (1) pathological diagnosis as breast cancer, (2) no distal metastasis, (3) radical surgery and R0 resection, (4) 1021 cancer-related gene testing of the tumor, (5) plasma MRD testing performed within 3 months after surgery and before adjuvant therapy. Results: A total of 18 patients in our cohort had detectable MRD. Baseline clinical characteristics found that patients with higher clinical stages were more likely to have detectable MRD. Analysis of single nucleotide variations and small insertions/deletions in baseline tumors showed that somatic mutations in MAP3K1, ATM, FLT1, GNAS, POLD1, SPEN, and WWP2 were significantly enriched in patients with postoperative MRD. Oncogenic signaling pathway analysis revealed that alteration of the Cell cycle pathway was more likely to occur in patients with detectable MRD (p=0.0125). Mutational signature analysis showed that defective DNA mismatch repair and activation-induced cytidine deaminase (AID) mediated somatic hypermutation (SHM) were associated with detectable MRD. Somatic copy number variation analysis found that amplification of CDKN2A, HOXB13, PPM1D, MPL, and VHL was significantly enriched in patients with detectable MRD. According to the OncoKB database, 77.8% (14/18) of patients with detectable MRD had U.S. Food and Drug Administration-approved mutational biomarkers and targeted therapy. Conclusion: For the first time, our study reports genomic characteristics of breast cancer patients with detectable MRD. The tumor biology is probably the main driver of detectable MRD. We also found the vast majority of patients with detectable MRD have the opportunity to access targeted therapy. Keywords: Breast cancer, MRD, Genomic character, Cell cycle pathway, Defective DNA mismatch repair, Activation-induced cytidine deaminase Figure 1. Somatic mutational landscape and clinical actionability of breast cancer patients with detectable MRD A. Oncoplot of top 20 genes altered in patients with detectable MRD. B. Forest plot of different mutant genes between patients with detectable and undetectable MRD. OR means odds ratio. Inf means infinity. ＊ means P-value < 0.05. C. The ten cancer-related signaling pathways were affected by somatic mutations both in patients with detectable and undetectable MRD. ＊ means P-value < 0.05. D&E. The fraction and etiology of each signature in patients with detectable (D) or undetectable (E) MRD. F. The number of patients with detectable MRD in different actionable alterations. gBRCA1 means germline BRCA1 mutation. sBRCA1 means somatic BRCA1 mutation. sBRCA2 means somatic BRCA2 mutation. TMB-H means TMB-High. Citation Format: Xu Yan, Shu Zhang, Lu zhou, Yan Jiang, Jing Xu, Gang Zhang, Lu Shen. Genomic characteristics and its therapeutic implications in breast cancer patients with detectable molecular residual disease [abstract]. In: Proceedings of the 2023 San Antonio Breast Cancer Symposium; 2023 Dec 5-9; San Antonio, TX. Philadelphia (PA): AACR; Cancer Res 2024;84(9 Suppl):Abstract nr PO1-13-05.

An entropy-based study on the mutational landscape of SARS-CoV-2 in USA: comparing different variants and revealing co-mutational behavior of proteins

COVID-19 emergency has pushed the international scientific community to use every resource to combat the spread of the virus, to understand its biology and predict its possible evolution in terms of new variants. Since the first SARS-CoV-2 virus nucleotide and amino acid sequences were made available, information theory was used to study how viral information content was changing over time and then trace the evolution of its mutational landscape. In this work we analyzed SARS-CoV-2 sequences collected mainly in the USA in a period from March 2020 until December 2022 and computed mutation profiles of viral proteins over time through an entropy-based approach using Shannon Entropy and Hellinger distance. This representation allows an at-a-glance view of the mutational landscape of viral proteins over time and can provide new insights on the evolution of the virus from different points of view. Non-structural proteins typically showed flat mutation profiles, characterized by a very low Average mutation Entropy, while accessory and structural proteins showed mostly non uniform and high mutation profiles, often coupled with the predominance of variants. Interestingly NSP2 protein, whose function is currently still debated, falls in the same branch of NSP14 and NSP10 in the phylogenetic tree of mutations constructed through correlations of mutation profiles, suggesting a co-evolution of those proteins and a possible functional link with each other. To the best of our knowledge this is the first study based on a massive amount of data (n=107,939,973) that analyzes from an entropy point of view the mutational landscape of SARS-CoV-2 over time and depicts a mutational temporal profile of each protein of the virus.

Mutational landscape of oral mucosal melanoma based on comprehensive cancer genomic profiling tests in a Japanese cohort

ObjectivesOral mucosal melanoma (OMM) is a rare but aggressive melanoma subtype. Due to its rarity, the genomic landscape of OMM remains unknown despite a relatively thorough understanding of the genetic profile of cutaneous melanoma (CM). In this study, we analyzed the genomic mutational profiles of Japanese patients with OMM and compared them with those of patients with nose/sinuses mucosal melanoma (NMM) and CM to identify potential therapeutic targets. Materials and MethodsWe extracted clinical and genomic information of patients with OMM (n = 15), NMM (n = 63), and CM (n = 413) who underwent comprehensive genomic profiling tests under the National Health Insurance between June 2019 and November 2023 from the Center for Cancer Genomics and Therapeutics database. ResultsThe most frequent genomic alteration identified in OMM was RICTOR (40%) followed by CDK4 (33.3%), MDM2 (33.3%), KDR (30%), KIT (26.7%), and NF1 (26.7%). CDK4 and MDM2 were co-amplified. Gene alterations in MYC and NRAS were the highest in patients with NMM, followed by those with CM, and no MYC alteration was observed in patients with OMM. BRAF V600 mutation, which is frequently observed in patients with CM (23.2%) were only present in 1.6% of patients with NMM and none in patients with OMM. ConclusionThis study clarified the genetic differences between OMM and NMM, and the first to report the frequent occurrence of RICTOR amplification in OMM. This analysis offers insights into the development of personalized therapeutics for OMM.

Clinical and genomic landscape of RAS mutations in gynecologic cancers.

We aimed to describe RAS mutations in gynecologic cancers as they relate to clinicopathologic and genomic features, survival, and therapeutic implications. Gynecologic cancers with available somatic molecular profiling data at our institution between February 2010 and August 2022 were included and grouped by RAS mutation status. Overall survival was estimated by Kaplan-Meier method, and multivariable analysis was performed using Cox proportional-hazards model. Of 3328 gynecologic cancers, 523 (15.7%) showed any RAS mutation. Patients with RAS-mutated tumors were younger (57 vs 60 years non-mutated), had higher prevalence of endometriosis (27.3% vs 16.9%), and lower grades (grade 1/2, 43.2% vs 8.1%, all p<0.0001). Highest prevalence of KRAS mutation was in mesonephric-like endometrial (100%, n=9/9), mesonephric-like ovarian (83.3%, n=5/6), mucinous ovarian (60.4%), and low-grade serous ovarian (44.4%) cancers. After adjustment for age, cancer type, and grade, RAS mutation was associated with worse overall survival (HR=1.3, p=0.001). Specific mutations were in KRAS (13.5%), NRAS (2.0%), and HRAS (0.51%), most commonly KRAS G12D (28.4%) and G12V (26.1%). Common co-mutations were PIK3CA (30.9%), PTEN(28.8%), ARID1A (28.0%), and TP53 (27.9%), of which 64.7% were actionable. RAS+MAPK pathway-targeted therapies were administered to 62 patients with RAS-mutated cancers. While overall survival was significantly higher with therapy (8.4 years [95%CI 5.5-12.0] vs 5.5 years [95%CI 4.6-6.6], HR=0.67, p=0.031), this effect did not persist in multivariable analysis. RAS mutations in gynecologic cancers have a distinct histopathologic distribution and may impact overall survival. PIK3CA, PTEN, and ARID1A are potentially actionable co-alterations. RAS pathway-targeted therapy should be considered.

A comparative analysis of the clinical and genetic profiles of blast phaseBCR::ABL1-negative myeloproliferative neoplasm and acute myeloid leukemia, myelodysplasia-related.

The classic Philadelphia chromosome-negative myeloproliferative neoplasms (Ph (-) MPNs), have variable potential for progression to the blast phase (MPN-BP) of the disease. Except initiated by distinct driver mutations, MPN-BP frequently carry similar genetic abnormalities defining acute myeloid leukemia myelodysplasia-related (AML-MR). Because of dissimilar initial pathogenesis, MPN-BP and AML-MR are retained under different disease categories. To determine if separately classifying these entities is justified, we compare MPN-BP with AML-MR patients based on mutational landscape and clinical parameters. 104 MPN-BP patients and 145 AML-MR patients were identified with available clinical, cytogenetic, and genetic data. AML-MR patients presented with a higher blast count (median, 51% vs. 30%) while MPN-BP patients had higher WBC counts, platelet counts and bone marrow cellularity (all p<0.0001). Patients with MPN-BP showed similar genetic mutations with similar mutation pattern (functional domain, hotspot and locus involved by the mutations) but a different mutation rate from AML-MR, with more frequent JAK2, CALR, MPL, ASXL1, IDH2, SETBP1 and SRSF2 mutations and less frequent TP53 and DNMT3A mutations. The overall survival (OS) of MPN-BP (OS post-BP-progression) is comparable to that of AML-MR (median OS, 9.5 months vs. 13.1 months, p=0.20). In addition, the subgroups of MPN-BP show similar OS as AML-MR. When harboring certain mutation such as TP53, ASXL1, DNMT3A, TET2, RUNX1, IDH1, IDH2, EZH2, U2AF1, BCOR and SRSF2, MPN-BP and AML-MR patients carrying the same somatic mutation show no difference in OS. MPN-BP and AML-MR harbor similar somatic mutations and clinical outcomes, suggesting a unified clinical disease entity.