A transcriptional MAPK Pathway Activity Score (MPAS) is a clinically relevant biomarker in multiple cancer types

Transcriptomic analysis of HEK293A cells with a CRISPR/Cas9-mediated TDP1 knockout

Tyrosyl-DNA phosphodiesterase 1 (TDP1) is a human DNA repair protein. It is a member of the phospholipase D family based on structural similarity. TDP1 is a key enzyme of the repair of stalled topoisomerase 1 (TOP1)-DNA complexes. Previously, with the CRISPR/Cas9 method, we obtained HEK293A cells with a homozygous knockout of the TDP1 gene and used the TDP1 knockout cells as a cellular model for studying mechanisms of action of an anticancer therapy. In the present work, we hypothesized that the TDP1 knockout would alter the expression of DNA repair-related genes. By transcriptomic analysis, we investigated for the first time the effect of the TDP1 gene knockout on genes' expression changes in the human HEK293A cell line. We obtained original data implying a role of TDP1 in other processes besides the repair of the DNA-TOP1 complex. Differentially expressed gene analysis revealed that TDP1 may participate in cell adhesion and communication, spermatogenesis, mitochondrial function, neurodegeneration, a cytokine response, and the MAPK signaling pathway.

Genomic and Transcriptomic Landscape of RET Wild-Type Medullary Thyroid Cancer and Potential Use of Mitogen-Activated Protein Kinase-Targeted Therapy

BACKGROUND

About 75% of medullary thyroid cancers (MTCs) are sporadic with 45% to 70% being driven by a RET mutation. Selpercatinib is an approved treatment for RET-mutated (mut RET ) MTC; however, treatments are needed for wild-type RET MTC (wt RET ). Genomic alterations and transcriptomic signatures of wt RET MTC may reveal new therapeutic insights.

STUDY DESIGN

We did a retrospective analysis of MTC samples submitted for DNA/RNA sequencing and programmed cell death ligand 1 expression using immunohistochemistry at a Clinical Laboratory Improvement Amendments/College of American Pathologists-certified laboratory. Tumor microenvironment immune cell fractions were estimated using RNA deconvolution (quanTIseq). Transcriptomic signatures of inflammation and MAP kinase pathway activation scores were calculated. Mann-Whitney U, chi-square, and Fisher's exact tests were applied (p values adjusted for multiple comparisons).

RESULTS

The 160-patient cohort included 108 mut RET and 52 wt RET MTC samples. wt RET tumors frequently harbored mitogen-activated protein kinase (MAPK) pathway mutations, including HRAS (42.31%), KRAS (15.7%), NF1 (6.7%), and BRAF (2%), whereas only 1 MAPK pathway mutation ( NF1 ) was identified among mut RET MTC. Recurrent mutations seen in wt RET MTC included MGA , VHL, APC , STK11 , and NFE2L2 . Increased transcriptional activation of the MAPK pathway was observed in patients with wt RET harboring mutations in MAPK genes. Although the frequency of programmed cell death ligand 1 expression was similar in wt RET and mut RET (10.2% vs 7%, p = 0.531), wt RET tumors were more often tumor mutational burden high (7.7% vs 0%, p = 0.011), and wt RET MTC exhibited higher expression of immune checkpoint genes.

CONCLUSIONS

We identified molecular alterations and immune-related features that distinguish wt RET from mut RET MTC. Although RET mutation drives MTC in the absence of other alterations, we showed that wt RET MTC frequently harbors MAPK pathway mutations. These findings may indicate a potential basis for MAPK-targeted therapy, possibly in combination with immuno-oncology agents for selected patients with wt RET MTC.

Capmatinib is an effective treatment for MET-fusion driven pediatric high-grade glioma and synergizes with radiotherapy

BACKGROUND

Pediatric-type diffuse high-grade glioma (pHGG) is the most frequent malignant brain tumor in children and can be subclassified into multiple entities. Fusion genes activating the MET receptor tyrosine kinase often occur in infant-type hemispheric glioma (IHG) but also in other pHGG and are associated with devastating morbidity and mortality.

METHODS

To identify new treatment options, we established and characterized two novel orthotopic mouse models harboring distinct MET fusions. These included an immunocompetent, murine allograft model and patient-derived orthotopic xenografts (PDOX) from a MET-fusion IHG patient who failed conventional therapy and targeted therapy with cabozantinib. With these models, we analyzed the efficacy and pharmacokinetic properties of three MET inhibitors, capmatinib, crizotinib and cabozantinib, alone or combined with radiotherapy.

RESULTS

Capmatinib showed superior brain pharmacokinetic properties and greater in vitro and in vivo efficacy than cabozantinib or crizotinib in both models. The PDOX models recapitulated the poor efficacy of cabozantinib experienced by the patient. In contrast, capmatinib extended survival and induced long-term progression-free survival when combined with radiotherapy in two complementary mouse models. Capmatinib treatment increased radiation-induced DNA double-strand breaks and delayed their repair.

CONCLUSIONS

We comprehensively investigated the combination of MET inhibition and radiotherapy as a novel treatment option for MET-driven pHGG. Our seminal preclinical data package includes pharmacokinetic characterization, recapitulation of clinical outcomes, coinciding results from multiple complementing in vivo studies, and insights into molecular mechanism underlying increased efficacy. Taken together, we demonstrate the groundbreaking efficacy of capmatinib and radiation as a highly promising concept for future clinical trials.

Rebound growth of BRAF mutant pediatric glioma cells after MAPKi withdrawal is associated with MAPK reactivation and secretion of microglia-recruiting cytokines

INTRODUCTION

Patients with pediatric low-grade gliomas (pLGGs), the most common primary brain tumors in children, can often benefit from MAPK inhibitor (MAPKi) treatment. However, rapid tumor regrowth, also referred to as rebound growth, may occur once treatment is stopped, constituting a significant clinical challenge.

METHODS

Four patient-derived pediatric glioma models were investigated to model rebound growth in vitro based on viable cell counts in response to MAPKi treatment and withdrawal. A multi-omics dataset (RNA sequencing and LC-MS/MS based phospho-/proteomics) was generated to investigate possible rebound-driving mechanisms. Following in vitro validation, putative rebound-driving mechanisms were validated in vivo using the BT-40 orthotopic xenograft model.

RESULTS

Of the tested models, only a BRAFV600E-driven model (BT-40, with additional CDKN2A/Bdel) showed rebound growth upon MAPKi withdrawal. Using this model, we identified a rapid reactivation of the MAPK pathway upon MAPKi withdrawal in vitro, also confirmed in vivo. Furthermore, transient overactivation of key MAPK molecules at transcriptional (e.g. FOS) and phosphorylation (e.g. pMEK) levels, was observed in vitro. Additionally, we detected increased expression and secretion of cytokines (CCL2, CX3CL1, CXCL10 and CCL7) upon MAPKi treatment, maintained during early withdrawal. While increased cytokine expression did not have tumor cell intrinsic effects, presence of these cytokines in conditioned media led to increased attraction of microglia cells in vitro.

CONCLUSION

Taken together, these data indicate rapid MAPK reactivation upon MAPKi withdrawal as a tumor cell intrinsic rebound-driving mechanism. Furthermore, increased secretion of microglia-recruiting cytokines may play a role in treatment response and rebound growth upon withdrawal, warranting further evaluation.

A phase 1 trial of the MEK inhibitor selumetinib in combination with pembrolizumab for advanced or metastatic solid tumors

MEK inhibitors have immunomodulatory activity and potential for synergistic activity when combined with PD-1 inhibitors. We evaluated selumetinib (inhibitor of MEK1/2) plus pembrolizumab (anti‒PD-1 antibody) in patients with advanced/metastatic solid tumors. In this phase 1b study, adults with previously treated advanced/metastatic solid tumors received pembrolizumab 200 mg intravenously every 3 weeks plus selumetinib on days 1‒14 per 3-week cycle (2 weeks on/1 week off); selumetinib dosing began at 50 mg orally twice daily with escalation in 25 mg increments for ≤ 35 cycles. Primary endpoints were dose-limiting toxicities (DLTs), adverse events (AEs), and treatment discontinuations due to AEs. Thirty-two patients were enrolled. Dose escalation was completed up to selumetinib 125 mg twice daily. The target DLT rate of 30% was not reached at any dose level. In the selumetinib 100 mg group, 2/11 patients (18.2%) experienced DLTs (n = 1 grade 3 diarrhea, n = 1 grade 3 fatigue). In the selumetinib 125 mg group, 3/14 (21.4%) experienced DLTs (n = 1 grade 2 retinal detachment, n = 1 grade 3 retinopathy, n = 1 grade 3 stomatitis). Dose-related changes in pharmacokinetic exposures were observed for selumetinib and N-desmethyl selumetinib up to 100 mg (saturation at 125 mg). Two patients achieved partial responses (1 each with selumetinib 75 mg and 125 mg) for an objective response rate of 6%. The study was stopped early because of insufficient efficacy. Although the target DLT rate was not reached at any dose level and no new safety signals were identified, selumetinib plus pembrolizumab had limited antitumor activity in this population. Trial registration: ClinicalTrials.gov , NCT03833427.

Molecular Characterization and Therapeutic Opportunities in KRAS Wildtype Pancreatic Ductal Adenocarcinoma

PURPOSE

To investigate the molecular characteristics of and potential for precision medicine in KRAS wildtype pancreatic ductal adenocarcinoma (PDAC).

PATIENTS AND METHODS

We investigated 27 patients with KRASWT PDAC at our institution. Clinical data were obtained via chart review. Tumor specimens for each subject were interrogated for somatic single nucleotide variants, insertion and deletions, and copy number variants by DNA sequencing. Gene fusions were detected from RNA-seq. A patient-derived organoid (PDO) was developed from a patient with a MET translocation and expanded ex vivo to predict therapeutic sensitivity prior to enrollment in a phase 2 clinical trial.

RESULTS

Transcriptomic analysis showed our cohort may be stratified by the relative gene expression of the KRAS signaling cascade. The PDO derived from our patient harboring a TFG-MET rearrangement was found to have in vitro sensitivity to the multi-tyrosine kinase inhibitor crizotinib. The patient was enrolled in the phase 2 SPARTA clinical trial and received monotherapy with vebrelitinib, a c-MET inhibitor, and achieved a partial and durable response.

CONCLUSIONS

KRASWT PDAC is molecularly distinct from KRASMUT and enriched with potentially actionable genetic variants. In our study, transcriptomic profiling revealed that the KRAS signaling cascade may play a key role in KRASWT PDAC. Our report of a KRASWT PDAC patient with TFG-MET rearrangement who responded to a cMET inhibitor further supports the pursuit of precision oncology in this sub-population. Identification of targetable mutations, perhaps through approaches like RNA-seq, can help enable precision-driven approaches to select optimal treatment based on tumor characteristics.

Specific oncogene activation of the cell of origin in mucosal melanoma

Mucosal melanoma (MM) is a deadly cancer derived from mucosal melanocytes. To test the consequences of MM genetics, we developed a zebrafish model in which all melanocytes experienced CCND1 expression and loss of PTEN and TP53. Surprisingly, melanoma only developed from melanocytes lining internal organs, analogous to the location of patient MM. We found that zebrafish MMs had a unique chromatin landscape from cutaneous melanoma. Internal melanocytes could be labeled using a MM-specific transcriptional enhancer. Normal zebrafish internal melanocytes shared a gene expression signature with MMs. Patient and zebrafish MMs have increased migratory neural crest gene and decreased antigen presentation gene expression, consistent with the increased metastatic behavior and decreased immunotherapy sensitivity of MM. Our work suggests the cell state of the originating melanocyte influences the behavior of derived melanomas. Our animal model phenotypically and transcriptionally mimics patient tumors, allowing this model to be used for MM therapeutic discovery.

Direct and selective pharmacological disruption of the YAP-TEAD interface by IAG933 inhibits Hippo-dependent and RAS-MAPK-altered cancers

The YAP-TEAD protein-protein interaction mediates YAP oncogenic functions downstream of the Hippo pathway. To date, available YAP-TEAD pharmacologic agents bind into the lipid pocket of TEAD, targeting the interaction indirectly via allosteric changes. However, the consequences of a direct pharmacological disruption of the interface between YAP and TEADs remain largely unexplored. Here, we present IAG933 and its analogs as potent first-in-class and selective disruptors of the YAP-TEAD protein-protein interaction with suitable properties to enter clinical trials. Pharmacologic abrogation of the interaction with all four TEAD paralogs resulted in YAP eviction from chromatin and reduced Hippo-mediated transcription and induction of cell death. In vivo, deep tumor regression was observed in Hippo-driven mesothelioma xenografts at tolerated doses in animal models as well as in Hippo-altered cancer models outside mesothelioma. Importantly this also extended to larger tumor indications, such as lung, pancreatic and colorectal cancer, in combination with RTK, KRAS-mutant selective and MAPK inhibitors, leading to more efficacious and durable responses. Clinical evaluation of IAG933 is underway.

Predictive value of CDC37 gene expression for targeted therapy in metastatic colorectal cancer

BACKGROUND

CDC37 is a key determinant of client kinase recruitment to the HSP90 chaperoning system. We hypothesized that kinase-specific dependency on CDC37 alters the efficacy of targeted therapies for metastatic colorectal cancer (mCRC).

MATERIAL AND METHODS

Two independent mCRC cohorts were analyzed to compare the survival outcomes between CDC37-high and CDC37-low patients (stratified by the median cutoff values): the CALGB/SWOG 80405 trial (226 and 207 patients receiving first-line bevacizumab- and cetuximab-containing chemotherapies, respectively) and Japanese retrospective (50 refractory patients receiving regorafenib) cohorts. A dataset of specimens submitted to a commercial CLIA-certified laboratory was utilized to characterize molecular profiles of CDC37-high (top quartile, N = 5055) and CDC37-low (bottom quartile, N = 5055) CRCs.

RESULTS

In the bevacizumab-treated group, CDC37-high patients showed significantly better progression-free survival (PFS) (median 13.3 vs 9.6 months, hazard ratio [HR] 0.59, 95% confidence interval [CI] 0.44-0.79, p < 0.01) than CDC37-low patients. In the cetuximab-treated group, CDC37-high and CDC37-low patients had similar outcomes. In the regorafenib-treated group, CDC37-high patients showed significantly better overall survival (median 11.3 vs 6.0 months, HR 0.24, 95% CI 0.11-0.54, p < 0.01) and PFS (median 3.5 vs 1.9 months, HR 0.51, 95% CI 0.28-0.94, p = 0.03). Comprehensive molecular profiling revealed that CDC37-high CRCs were associated with higher VEGFA, FLT1, and KDR expressions and activated hypoxia signature.

CONCLUSIONS

CDC37-high mCRC patients derived more benefit from anti-VEGF therapies, including bevacizumab and regorafenib, but not from cetuximab. Molecular profiles suggested that such tumors were dependent on angiogenesis-relating pathways.

Epidermal growth factor receptor (EGFR) is a target of the tumor-suppressor E3 ligase FBXW7

FBXW7 is an E3 ubiquitin ligase that targets proteins for proteasome-mediated degradation and is mutated in various cancer types. Here, we use CRISPR base editors to introduce different FBXW7 hotspot mutations in human colon organoids. Functionally, FBXW7 mutation reduces EGF dependency of organoid growth by ~10,000-fold. Combined transcriptomic and proteomic analyses revealed increased EGFR protein stability in FBXW7 mutants. Two distinct phosphodegron motifs reside in the cytoplasmic tail of EGFR. Mutations in these phosphodegron motifs occur in human cancer. CRISPR-mediated disruption of the phosphodegron motif at T693 reduced EGFR degradation and EGF growth factor dependency. FBXW7 mutant organoids showed reduced sensitivity to EGFR-MAPK inhibitors. These observations were further strengthened in CRC-derived organoid lines and validated in a cohort of patients treated with panitumumab. Our data imply that FBXW7 mutations reduce EGF dependency by disabling EGFR turnover.

Liprin-α1 contributes to oncogenic MAPK signaling by counteracting ERK activity

PTPRF interacting protein alpha 1 (PPFIA1) encodes for liprin-α1, a member of the leukocyte common antigen-related protein tyrosine phosphatase (LAR-RPTPs)-interacting protein family. Liprin-α1 localizes to adhesive and invasive structures in the periphery of cancer cells, where it modulates migration and invasion in head and neck squamous cell carcinoma (HNSCC) and breast cancer. To study the possible role of liprin-α1 in anticancer drug responses, we screened a library of oncology compounds in cell lines with high endogenous PPFIA1 expression. The compounds with the highest differential responses between high PPFIA1-expressing and silenced cells across cell lines were inhibitors targeting mitogen-activated protein kinase kinase (MEK)/extracellular signal-regulated kinases (ERK) signaling. KRAS proto-oncogene, GTPase (KRAS)-mutated MDA-MB-231 cells were more resistant to trametinib upon PPFIA1 knockdown compared with control cells. In contrast, liprin-α1-depleted HNSCC cells with low RAS activity showed a context-dependent response to MEK/ERK inhibitors. Importantly, we showed that liprin-α1 depletion leads to increased p-ERK1/2 levels in all our studied cell lines independent of KRAS mutational status, suggesting a role of liprin-α1 in the regulation of MAPK oncogenic signaling. Furthermore, liprin-α1 depletion led to more pronounced redistribution of RAS proteins to the cell membrane. Our data suggest that liprin-α1 is an important contributor to oncogenic RAS/MAPK signaling, and the status of liprin-α1 may assist in predicting drug responses in cancer cells in a context-dependent manner.

First-Line Ipatasertib, Atezolizumab, and Taxane Triplet for Metastatic Triple-Negative Breast Cancer: Clinical and Biomarker Results

PURPOSE

To evaluate a triplet regimen combining immune checkpoint blockade, AKT pathway inhibition, and (nab-) paclitaxel as first-line therapy for locally advanced/metastatic triple-negative breast cancer (mTNBC).

PATIENTS AND METHODS

The single-arm CO40151 phase Ib study (NCT03800836), the single-arm signal-seeking cohort of IPATunity130 (NCT03337724), and the randomized phase III IPATunity170 trial (NCT04177108) enrolled patients with previously untreated mTNBC. Triplet therapy comprised intravenous atezolizumab 840 mg (days 1 and 15), oral ipatasertib 400 mg/day (days 1-21), and intravenous paclitaxel 80 mg/m2 (or nab-paclitaxel 100 mg/m2; days 1, 8, and 15) every 28 days. Exploratory translational research aimed to elucidate mechanisms and molecular markers of sensitivity and resistance.

RESULTS

Among 317 patients treated with the triplet, efficacy ranged across studies as follows: median progression-free survival (PFS) 5.4 to 7.4 months, objective response rate 44% to 63%, median duration of response 5.6 to 11.1 months, and median overall survival 15.7 to 28.3 months. The safety profile was consistent with the known toxicities of each agent. Grade ≥3 adverse events were more frequent with the triplet than with doublets or single-agent paclitaxel. Patients with PFS >10 months were characterized by NF1, CCND3, and PIK3CA alterations and increased immune pathway activity. PFS <5 months was associated with CDKN2A/CDKN2B/MTAP alterations and lower predicted phosphorylated AKT-S473 levels.

CONCLUSIONS

In patients with mTNBC receiving an ipatasertib/atezolizumab/taxane triplet regimen, molecular characteristics may identify those with particularly favorable or unfavorable outcomes, potentially guiding future research efforts.

Lineage-specific intolerance to oncogenic drivers restricts histological transformation

Lung adenocarcinoma (LUAD) and small cell lung cancer (SCLC) are thought to originate from different epithelial cell types in the lung. Intriguingly, LUAD can histologically transform into SCLC after treatment with targeted therapies. In this study, we designed models to follow the conversion of LUAD to SCLC and found that the barrier to histological transformation converges on tolerance to Myc, which we implicate as a lineage-specific driver of the pulmonary neuroendocrine cell. Histological transformations are frequently accompanied by activation of the Akt pathway. Manipulating this pathway permitted tolerance to Myc as an oncogenic driver, producing rare, stem-like cells that transcriptionally resemble the pulmonary basal lineage. These findings suggest that histological transformation may require the plasticity inherent to the basal stem cell, enabling tolerance to previously incompatible oncogenic driver programs.

Compartment-specific multiomic profiling identifies SRC and GNAS as candidate drivers of epithelial-to-mesenchymal transition in ovarian carcinosarcoma

BACKGROUND

Ovarian carcinosarcoma (OCS) is an exceptionally aggressive and understudied ovarian cancer type harbouring distinct carcinomatous and sarcomatous compartments. Here, we seek to identify shared and compartment-specific events that may represent potential therapeutic targets and candidate drivers of sarcomatous compartment formation through epithelial-to-mesenchymal transition (EMT).

METHODS

We performed multiomic profiling (exome sequencing, RNA-sequencing, microRNA profiling) of paired carcinomatous and sarcomatous components in 12 OCS cases.

RESULTS

While paired sarcomatous and carcinomatous compartments demonstrate substantial genomic similarities, multiple loci are recurrently copy number-altered between components; regions containing GNAS and SRC are recurrently gained within the sarcomatous compartment. CCNE1 gain is a common event in OCS, occurring more frequently than in high grade serous ovarian carcinoma (HGSOC). Transcriptomic analysis suggests increased MAPK activity and subtype switching toward poor prognosis HGSOC-derived transcriptomic subtypes within the sarcomatous component. The two compartments show global differences in microRNA profiles, with differentially expressed microRNAs targeting EMT-related genes (SIRT1, ZEB2) and regulators of pro-tumourigenic pathways (TGFβ, NOTCH); chrX is a highly enriched target of these microRNAs and is also frequently deleted across samples. The sarcomatous component harbours significantly fewer CD8-positive cells, suggesting poorer immune engagement.

CONCLUSION

CCNE1 gain and chrX loss are frequent in OCS. SRC gain, increased GNAS expression and microRNA dysregulation represent potential mechanisms driving sarcomatous compartment formation.

Mesothelin expression correlates with elevated inhibitory immune activity in patients with colorectal cancer

The expression of the protein Mesothelin (MSLN) is highly variable in several malignancies including colorectal cancer (CRC) and high levels are associated with aggressive clinicopathological features and worse patient survival. CRC is both a common and deadly cancer; being the third most common in incidence and second most common cause of cancer related death. While systemic therapy remains the primary therapeutic option for most patients with stage IV (metastatic; m) CRC, their disease eventually becomes treatment refractory, and 85% succumb within 5 years. Microsatellite-stable (MSS) CRC tumors, which affect more than 90% of patients with mCRC, are generally refractory to immunotherapeutic interventions. In our current work, we characterize MSLN levels in CRC, specifically correlating expression with clinical outcomes in relevant CRC subtypes and explore how MSLN expression impacts the status of immune activation and suppression in the peritumoral microenvironment. High MSLN expression is highly prevalent in CMS1 and CMS4 CRC subtypes as well as in mCRC tissue and correlates with higher gene mutation rates across the patient cohorts. Further, MSLN-high patients exhibit increased M1/M2 macrophage infiltration, PD-L1 staining, immune-inhibitory gene expression, enrichment in inflammatory, TGF-β, IL6/JAK/STAT3, IL2/STAT5 signaling pathways and mutation in KRAS and FBXW7. Together, these results suggest MSLN protein is a potential target for antigen-specific therapy and supports investigation into its tumorigenic effects to identify possible therapeutic interventions for patients with high MSLN expressing MSS CRC.

SOS1 inhibition enhances the efficacy of and delays resistance to G12C inhibitors in lung adenocarcinoma

Clinical effectiveness of KRAS G12C inhibitors (G12Cis) is limited both by intrinsic and acquired resistance, necessitating the development of combination approaches. We found that targeting proximal receptor tyrosine kinase (RTK) signaling using the SOS1 inhibitor (SOS1i) BI-3406 both enhanced the potency of and delayed resistance to G12Ci treatment, but the extent of SOS1i effectiveness was modulated by both SOS2 expression and the specific mutational landscape. SOS1i enhanced the efficacy of G12Ci and limited rebound RTK/ERK signaling to overcome intrinsic/adaptive resistance, but this effect was modulated by SOS2 protein levels. Survival of drug-tolerant persister (DTP) cells within the heterogeneous tumor population and/or acquired mutations that reactivate RTK/RAS signaling can lead to outgrowth of tumor initiating cells (TICs) that drive therapeutic resistance. G12Ci drug tolerant persister cells showed a 2-3-fold enrichment of TICs, suggesting that these could be a sanctuary population of G12Ci resistant cells. SOS1i re-sensitized DTPs to G12Ci and inhibited G12C-induced TIC enrichment. Co-mutation of the tumor suppressor KEAP1 limits the clinical effectiveness of G12Cis, and KEAP1 and STK11 deletion increased TIC frequency and accelerated the development of acquired resistance to G12Ci in situ. SOS1i both delayed acquired G12Ci resistance and limited the total number of resistant colonies regardless of KEAP1 and STK11 mutational status. These data suggest that SOS1i could be an effective strategy to both enhance G12Ci efficacy and prevent G12Ci resistance regardless of co-mutations.

Integrative whole-genome and transcriptome analysis of HER2-amplified metastatic breast cancer

BACKGROUND

In breast cancer, the advent of anti-HER2 therapies has made HER2+ tumors a highly relevant subgroup. However, the exact characteristics which prohibit clinical response to anti-HER2 therapies and drive disease progression are not yet fully known. Integrative whole-genome and transcriptomic sequencing data from both primary and metastatic HER2-positive breast cancer will enhance our understanding of underlying biological processes.

METHODS

Here, we used WGS and RNA sequencing data of 700 metastatic breast tumors, of which 68 being HER2+, to search for specific genomic features of HER2+ disease and therapy resistance. Furthermore, we integrated results with transcriptomic data to associate tumors exhibiting a HER2+-specific gene expression profile with ERBB2 mutation status, prior therapy and relevant gene expression signatures.

RESULTS

Overall genomic profiles of primary and metastatic HER2+ breast cancers were similar, and no specific acquired genomics traits connected to prior anti-HER2 treatment were observed. However, specific genomic features were predictive of progression-free survival on post-biopsy anti-HER2 treatment. Furthermore, a HER2-driven expression profile grouped HER2-amplified tumors with ERBB2-mutated cases and cases without HER2 alterations. The latter were reported as ER positive in primary disease, but the metastatic biopsy showed low ESR1 expression and upregulation of the MAPK pathway, suggesting transformation to ER independence.

CONCLUSIONS

In summary, although the quantity of variants increased throughout HER2-positive breast cancer progression, the genomic composition remained largely consistent, thus yielding no new major processes beside those already operational in primary disease. Our results suggest that integrated genomic and transcriptomic analyses may be key in establishing therapeutic options.

A short-term three dimensional culture-based drug sensitivity test is feasible for malignant bone tumors

The feasibility of a short-term, three-dimensional (3D) culture-based drug sensitivity test (DST) for surgically resected malignant bone tumors, including osteosarcoma (OS), was evaluated utilizing two OS cell line (KCS8 or KCS9)-derived xenograft (CDX) models. Twenty-three (KCS8) or 39 (KCS9) of 60 tested drugs were likely effective in OS cells derived from a cell line before xenografting. Fewer drugs (19: KCS8, 26: KCS9) were selected as effective drugs in cells derived from a CDX tumor, although the drug sensitivities of 60 drugs significantly correlated between both types of samples. The drug sensitivity of a CDX tumor was not significantly altered after the depletion of non-tumorous components in the sample. In a surgically resected metastatic tumor obtained from a patient with OS, for whom a cancer genome profiling test detected a pathogenic PIK3CA mutation, DST identified mTOR and AKT inhibitors as effective drugs. Of two CDX and six clinical samples of OS and Ewing's sarcoma, DST identified proteasome inhibitors (bortezomib, carfilzomib) and CEP-701 as potentially effective drugs in common. This unique method of in vitro drug testing using 3D-cell cultures is feasible in surgically resected tissues of metastatic malignant bone tumors.

Expressed Barcoding Enables High-Resolution Tracking of the Evolution of Drug Tolerance

For a majority of patients with non-small cell lung cancer with EGFR mutations, treatment with EGFR inhibitors (EGFRi) induces a clinical response. Despite this initial reduction in tumor size, residual disease persists that leads to disease relapse. Elucidating the preexisting biological differences between sensitive cells and surviving drug-tolerant persister cells and deciphering how drug-tolerant cells evolve in response to treatment could help identify strategies to improve the efficacy of EGFRi. In this study, we tracked the origins and clonal evolution of drug-tolerant cells at a high resolution by using an expressed barcoding system coupled with single-cell RNA sequencing. This platform enabled longitudinal profiling of gene expression and drug sensitivity in response to EGFRi across a large number of clones. Drug-tolerant cells had higher expression of key survival pathways such as YAP and EMT at baseline and could also differentially adapt their gene expression following EGFRi treatment compared with sensitive cells. In addition, drug combinations targeting common downstream components (MAPK) or orthogonal factors (chemotherapy) showed greater efficacy than EGFRi alone, which is attributable to broader targeting of the heterogeneous EGFRi-tolerance mechanisms present in tumors. Overall, this approach facilitates thorough examination of clonal evolution in response to therapy that could inform the development of improved diagnostic approaches and treatment strategies for targeting drug-tolerant cells.

SIGNIFICANCE

The evolution and heterogeneity of EGFR inhibitor tolerance are identified in a large number of clones at enhanced cellular and temporal resolution using an expressed barcode technology coupled with single-cell RNA sequencing.

Comparative Genomic Analysis of Pancreatic Acinar Cell Carcinoma (PACC) and Pancreatic Ductal Adenocarcinoma (PDAC) Unveils New Actionable Genomic Aberrations in PACC

PURPOSE

Pure pancreatic acinar cell carcinomas (PACC) are rare malignancies with no established treatment. PACC demonstrates significant genetic intertumoral heterogeneity with multiple pathways involved, suggesting using targeted cancer therapeutics to treat this disease. We aggregated one of the largest datasets of pure PACC to examine the genomic variability and explore patient-specific therapeutic targets.

EXPERIMENTAL DESIGN

PACC specimens (n = 51) underwent next-generation sequencing of DNA (n = 29) or whole exome (n = 22) and RNA (whole transcriptome, n = 29) at a commercial laboratory. We performed comparative analyses of a genomic cohort of pancreatic ductal adenocarcinomas (PDAC; n = 4,205). In parallel, we conducted a retrospective review of patients with PACC treated at Huntsman Cancer Institute (HCI).

RESULTS

The real-world dataset included samples from 51 patients with PACC. We found key molecular differences between pure PACC and PDAC, highlighting the unique characteristics of pure PACC. Major differences in PACC include lower MAPK signaling and less stromal cell abundance compared with PDAC. Pure PACC showed genomic loss-of-heterozygosity to largely coincide with mutations in BRCA1, BRCA2, and PALB2. Of the 7 patients treated at HCI, one had a tumor that harbored a BRAF-V600E mutation. Leveraging precision oncology, this patient is being treated with encorafenib plus binimetinib, achieving an exceptionally durable and ongoing complete response of more than 3 years.

CONCLUSIONS

There are major differences between PACC and PDAC, including downregulation of the MAPK signaling pathway, and less stromal cell abundance. In addition, genomic characterization of pure PACC revealed frequent targetable alterations, which can guide patient treatment.

Combined CDK4/6 and ERK1/2 Inhibition Enhances Antitumor Activity in NF1-Associated Plexiform Neurofibroma

PURPOSE

Plexiform neurofibromas (PNF) are peripheral nerve sheath tumors that cause significant morbidity in persons with neurofibromatosis type 1 (NF1), yet treatment options remain limited. To identify novel therapeutic targets for PNF, we applied an integrated multi-omic approach to quantitatively profile kinome enrichment in a mouse model that has predicted therapeutic responses in clinical trials for NF1-associated PNF with high fidelity.

EXPERIMENTAL DESIGN

Utilizing RNA sequencing combined with chemical proteomic profiling of the functionally enriched kinome using multiplexed inhibitor beads coupled with mass spectrometry, we identified molecular signatures predictive of response to CDK4/6 and RAS/MAPK pathway inhibition in PNF. Informed by these results, we evaluated the efficacy of the CDK4/6 inhibitor, abemaciclib, and the ERK1/2 inhibitor, LY3214996, alone and in combination in reducing PNF tumor burden in Nf1flox/flox;PostnCre mice.

RESULTS

Converging signatures of CDK4/6 and RAS/MAPK pathway activation were identified within the transcriptome and kinome that were conserved in both murine and human PNF. We observed robust additivity of the CDK4/6 inhibitor, abemaciclib, in combination with the ERK1/2 inhibitor, LY3214996, in murine and human NF1(Nf1) mutant Schwann cells. Consistent with these findings, the combination of abemaciclib (CDK4/6i) and LY3214996 (ERK1/2i) synergized to suppress molecular signatures of MAPK activation and exhibited enhanced antitumor activity in Nf1flox/flox;PostnCre mice in vivo.

CONCLUSIONS

These findings provide rationale for the clinical translation of CDK4/6 inhibitors alone and in combination with therapies targeting the RAS/MAPK pathway for the treatment of PNF and other peripheral nerve sheath tumors in persons with NF1.

ASCL1 is activated downstream of the ROR2/CREB signaling pathway to support lineage plasticity in prostate cancer

Lineage plasticity is a form of therapy-induced drug resistance. In prostate cancer, androgen receptor (AR) pathway inhibitors potentially lead to the accretion of tumor relapse with loss of AR signaling and a shift from a luminal state to an alternate program. However, the molecular and signaling mechanisms orchestrating the development of lineage plasticity under the pressure of AR-targeted therapies are not fully understood. Here, a survey of receptor tyrosine kinases (RTKs) identifies ROR2 as the top upregulated RTK following AR pathway inhibition, which feeds into lineage plasticity by promoting stem-cell-like and neuronal networks. Mechanistically, ROR2 activates the ERK/CREB signaling pathway to modulate the expression of the lineage commitment transcription factor ASCL1. Collectively, our findings nominate ROR2 as a potential therapeutic target to reverse the ENZ-induced plastic phenotype and potentially re-sensitize tumors to AR pathway inhibitors.

ICAM1 antibody drug conjugates exert potent antitumor activity in papillary and anaplastic thyroid carcinoma

Treatment options for anaplastic thyroid cancer (ATC) and refractory papillary thyroid carcinoma (PTC) are limited and outcomes remain poor. In this study, we determined via bioinformatic expression analyses and immunohistochemistry staining that intercellular adhesion molecule-1(ICAM1) is an attractive target for ATC and PTC. We designed and engineered two ICAM1-directed antibody-drug conjugate (I1-MMAE and I1-DXd), both of which potently and selectively ablate multiple human ATC and PTC cell lines without affecting non-plastic cells in vitro. Furthermore, I1-MMAE and I1-DXd mediated a potent tumor regression in ATC and PTC xenograft models. To develop a precision medicine, we also explored magnetic resonance imaging (MRI) as a non-invasive biomarker detection method to quantitatively map ICAM1 antigen expression in heterogeneous thyroid tumors. Taken together, this study provides a strong rationale for the further development of I1-MMAE and I1-DXd as promising therapeutic candidates to treat advanced PTC and ATC.

Integrated analysis of canine soft tissue sarcomas identifies recurrent mutations in TP53, KMT genes and PDGFB fusions

Soft tissue sarcomas (STS) are a heterogenous group of mesenchymal tumors representing over 50 distinct types with overlapping histological features and non-specific anatomical locations. Currently, localized sarcomas are treated with surgery + / -  radiation in both humans and dogs with few molecularly targeted therapeutic options. However, to improve precision-based cancer therapy through trials in pet dogs with naturally occurring STS tumors, knowledge of genomic profiling and molecular drivers in both species is essential. To this purpose, we sought to characterize the transcriptomic and genomic mutation profiles of canine STS subtypes (fibrosarcoma, undifferentiated pleomorphic sarcoma, and peripheral nerve sheath tumors), by leveraging RNAseq, whole exome sequencing, immunohistochemistry, and drug assays. The most common driver mutations were in cell cycle/DNA repair (31%, TP53-21%) and chromatin organization/binding (41%, KMT2D-21%) genes. Similar to a subset of human sarcomas, we identified fusion transcripts of platelet derived growth factor B and collagen genes that predict sensitivity to PDGFR inhibitors. Transcriptomic profiling grouped these canine STS tumors into 4 clusters, one PNST group (H1), and 3 FSA groups selectively enriched for extracellular matrix interactions and PDFGB fusions (H2), homeobox transcription factors (H3), and elevated T-cell infiltration (H4). This multi-omics approach provides insights into canine STS sub-types at a molecular level for comparison to their human counterparts, to improve diagnosis, and may provide additional targets for chemo- and immuno-therapy.

Antitumor Efficacy of Dual Blockade with Encorafenib + Cetuximab in Combination with Chemotherapy in Human BRAFV600E-Mutant Colorectal Cancer

PURPOSE

Encorafenib + cetuximab (E+C) is an effective therapeutic option in chemorefractory BRAFV600E metastatic colorectal cancer (mCRC). However, there is a need to improve the efficacy of this molecular-targeted therapy and evaluate regimens suitable for untreated BRAFV600E in patients with mCRC.

EXPERIMENTAL DESIGN

We performed a series of in vivo studies using BRAFV600E mCRC tumor xenografts. Mice were randomized to receive 5-fluoruracil (5-FU), irinotecan, or oxaliplatin regimens (FOLFIRI or FOLFOX), (E+C) or the combination. Patients received long-term treatment until disease progression, with deescalation strategies used to mimic maintenance therapy. Transcriptomic changes after progression on cytotoxic chemotherapy or targeted therapy were assessed.

RESULTS

Antitumor activity of either FOLFIRI or E+C was better as first-line treatment as compared with second-line, with partial cross-resistance seen between a cytotoxic regimen and targeted therapy with an average 62% loss of efficacy for FOLFIRI after E+C and a 45% loss of efficacy of E+C after FOLFIRI (P < 0.001 for both). FOLFIRI-treated models had upregulation of epithelial-mesenchymal transition (EMT) and MAPK pathway activation, where E+C treated models had suppressed MAPK signaling. In contrast, with chemotherapy with E+C, EMT and MAPK signaling remained suppressed. FOLFOX or FOLFIRI, each in combination with E+C, were the most active first-line treatments as compared with E+C or to chemotherapy alone. Furthermore, FOLFOX in combination with E+C as first-line induction therapy, followed by E+C ± 5-FU as maintenance therapy, was the most effective strategy for long-term disease control.

CONCLUSIONS

These results support the combination of cytotoxic chemotherapy and molecular-targeted therapy as a promising therapeutic approach in the first-line treatment of BRAFV600E mCRC.

Capicua (CIC) mutations in gliomas in association with MAPK activation for exposing a potential therapeutic target

Gliomas are the most prevalent neurological cancer in the USA and care modalities are not able to effectively combat these aggressive malignancies. Identifying new, more effective treatments require a deep understanding of the complex genetic variations and relevant pathway associations behind these cancers. Drawing connections between gene mutations with a responsive genetic target can help drive therapy selections to enhance patient survival. We have performed extensive molecular profiling of the Capicua gene (CIC), a tumor and transcriptional suppressor gene, and its mutation prevalence in reference to MAPK activation within clinical glioma tissue. CIC mutations occur far more frequently in oligodendroglioma (52.1%) than in low-grade astrocytoma or glioblastoma. CIC-associated mutations were observed across all glioma subtypes, and MAPK-associated mutations were most prevalent in CIC wild-type tissue regardless of the glioma subtype. MAPK activation, however, was enhanced in CIC-mutated oligodendroglioma. The totality of our observations reported supports the use of CIC as a relevant genetic marker for MAPK activation. Identification of CIC mutations, or lack thereof, can assist in selecting, implementing, and developing MEK/MAPK-inhibitory trials to improve patient outcomes potentially.

An allosteric pan-TEAD inhibitor blocks oncogenic YAP/TAZ signaling and overcomes KRAS G12C inhibitor resistance

The Hippo pathway is a key growth control pathway that is conserved across species. The downstream effectors of the Hippo pathway, YAP (Yes-associated protein) and TAZ (transcriptional coactivator with PDZ-binding motif), are frequently activated in cancers to drive proliferation and survival. Based on the premise that sustained interactions between YAP/TAZ and TEADs (transcriptional enhanced associate domain) are central to their transcriptional activities, we discovered a potent small-molecule inhibitor (SMI), GNE-7883, that allosterically blocks the interactions between YAP/TAZ and all human TEAD paralogs through binding to the TEAD lipid pocket. GNE-7883 effectively reduces chromatin accessibility specifically at TEAD motifs, suppresses cell proliferation in a variety of cell line models and achieves strong antitumor efficacy in vivo. Furthermore, we uncovered that GNE-7883 effectively overcomes both intrinsic and acquired resistance to KRAS (Kirsten rat sarcoma viral oncogene homolog) G12C inhibitors in diverse preclinical models through the inhibition of YAP/TAZ activation. Taken together, this work demonstrates the activities of TEAD SMIs in YAP/TAZ-dependent cancers and highlights their potential broad applications in precision oncology and therapy resistance.

Genomic hallmarks and therapeutic implications of G0 cell cycle arrest in cancer

BACKGROUND

Therapy resistance in cancer is often driven by a subpopulation of cells that are temporarily arrested in a non-proliferative G0 state, which is difficult to capture and whose mutational drivers remain largely unknown.

RESULTS

We develop methodology to robustly identify this state from transcriptomic signals and characterise its prevalence and genomic constraints in solid primary tumours. We show that G0 arrest preferentially emerges in the context of more stable, less mutated genomes which maintain TP53 integrity and lack the hallmarks of DNA damage repair deficiency, while presenting increased APOBEC mutagenesis. We employ machine learning to uncover novel genomic dependencies of this process and validate the role of the centrosomal gene CEP89 as a modulator of proliferation and G0 arrest capacity. Lastly, we demonstrate that G0 arrest underlies unfavourable responses to various therapies exploiting cell cycle, kinase signalling and epigenetic mechanisms in single-cell data.

CONCLUSIONS

We propose a G0 arrest transcriptional signature that is linked with therapeutic resistance and can be used to further study and clinically track this state.

BH3 mimetics targeting BCL-XL impact the senescent compartment of pilocytic astrocytoma

BACKGROUND

Pilocytic astrocytoma (PA) is the most common pediatric brain tumor and a mitogen-activated protein kinase (MAPK)-driven disease. Oncogenic MAPK-signaling drives the majority of cells into oncogene-induced senescence (OIS). While OIS induces resistance to antiproliferative therapies, it represents a potential vulnerability exploitable by senolytic agents.

METHODS

We established new patient-derived PA cell lines that preserve molecular features of the primary tumors and can be studied in OIS and proliferation depending on expression or repression of the SV40 large T antigen. We determined expression of anti-apoptotic BCL-2 members in these models and primary PA. Dependence of senescent PA cells on anti-apoptotic BCL-2 members was investigated using a comprehensive set of BH3 mimetics.

RESULTS

Senescent PA cells upregulate BCL-XL upon senescence induction and show dependency on BCL-XL for survival. BH3 mimetics with high affinity for BCL-XL (BCL-XLi) reduce metabolic activity and induce mitochondrial apoptosis in senescent PA cells at nano-molar concentrations. In contrast, BH3 mimetics without BCL-XLi activity, conventional chemotherapy, and MEK inhibitors show no effect.

CONCLUSIONS

Our data demonstrate that BCL-XL is critical for survival of senescent PA tumor cells and provides proof-of-principle for the use of clinically available BCL-XL-dependent senolytics.

HER2+ endometrioid endometrial cancer possesses distinct molecular and immunologic features associated with a more active immune microenvironment and worse prognosis

OBJECTIVE

HER2 status is not routinely evaluated in endometrioid endometrial cancer (E-EMCA), though it is frequently overexpressed or amplified in high grade E-EMCA and uterine serous carcinoma. Defining characteristics and survival outcomes of HER2+ E-EMCA could reveal subsets of patients who may benefit from targeted therapies.

METHODS

2927 E-EMCA tumors from the Caris Life Sciences database were analyzed by next-generation sequencing and whole exome sequencing, whole transcriptome sequencing, and immunohistochemistry for molecular and genomic features in a CLIA/CAP-certified laboratory (Caris Life Sciences, Phoenix, AZ). HER2 status was determined by transcriptomic cutoff extrapolated from uterine serous carcinoma. The relationship between HER2 status and patient outcomes was determined by Kaplan-Meier analysis.

RESULTS

HER2 positivity was detected in 5.47% of E-EMCA. Differences in molecular alterations based on HER2 status were most apparent in microsatellite stable (MSS) tumors, which displayed increased TP53 mutations and loss of heterozygosity (LOH) and decreased PTEN and CTNNB1 mutations. HER2+ tumors had increased immune checkpoint gene expression and immune cell infiltration, particularly among MSS tumors. All HER2+ tumors displayed increased MAPK pathway activation scores (MPAS) and patients with HER2+ tumors experienced worse overall survival.

CONCLUSIONS

HER2 positivity in E-EMCA corresponds with a unique molecular landscape, particularly in MSS tumors. HER2+ tumors are also associated with increased MAPK pathway activation and exhibit features of a more active immune microenvironment. These findings suggest a potential benefit of HER2 and MAPK targeted therapies as well as immunotherapies in this patient population.

Optimisation Models for Pathway Activity Inference in Cancer

BACKGROUND

With advances in high-throughput technologies, there has been an enormous increase in data related to profiling the activity of molecules in disease. While such data provide more comprehensive information on cellular actions, their large volume and complexity pose difficulty in accurate classification of disease phenotypes. Therefore, novel modelling methods that can improve accuracy while offering interpretable means of analysis are required. Biological pathways can be used to incorporate a priori knowledge of biological interactions to decrease data dimensionality and increase the biological interpretability of machine learning models.

METHODOLOGY

A mathematical optimisation model is proposed for pathway activity inference towards precise disease phenotype prediction and is applied to RNA-Seq datasets. The model is based on mixed-integer linear programming (MILP) mathematical optimisation principles and infers pathway activity as the linear combination of pathway member gene expression, multiplying expression values with model-determined gene weights that are optimised to maximise discrimination of phenotype classes and minimise incorrect sample allocation.

RESULTS

The model is evaluated on the transcriptome of breast and colorectal cancer, and exhibits solution results of good optimality as well as good prediction performance on related cancer subtypes. Two baseline pathway activity inference methods and three advanced methods are used for comparison. Sample prediction accuracy, robustness against noise expression data, and survival analysis suggest competitive prediction performance of our model while providing interpretability and insight on key pathways and genes. Overall, our work demonstrates that the flexible nature of mathematical programming lends itself well to developing efficient computational strategies for pathway activity inference and disease subtype prediction.

Clinical results of the EVESOR trial, a multiparameter phase I trial of everolimus and sorafenib combination in solid tumors

PURPOSE

Everolimus (EVE) and sorafenib (SOR) combination was associated with synergistic activity in preclinical models. However, previous clinical studies were hampered by cumulative toxicities when both were given continuously. The academic EVESOR trial (NCT01932177) was designed to assess alternative doses and intermittent dosing schedules of EVE and SOR combination therapy to improve the benefit-risk ratio for patients with solid tumors.

METHODS

EVESOR is a multiparameter dose-escalation phase I trial investigating different doses and dosing schedules, with the final objective of generating data for modeling and simulation. Patients were allocated into continuous (A and B) or intermittent (C and D) schedules to determine the recommended phase II dose (RP2D). The clinical outcomes are presented here.

RESULTS

Forty-three patients were included from 2013 to 2019. Most of them had gynecological (25.6%), cholangiocarcinomas (23.2%), colorectal (14.0%), and breast cancers (11.6%). Dose-escalation up to EVE 10 mg QD and SOR 400 mg BID was possible on intermittent schedules. Five dose-limiting toxicities were observed, and dose reductions were required in 39.5% patients, stabilizing at EVE 5 mg and SOR 200 mg BID for 58.1% of them. The overall response rate was 6.3%, and disease control rate was 75.0%. The median progression-free survival (PFS) was 3.6 months. The longest median PFS were observed in cholangiocarcinomas (9.9 months), and gynecological adenocarcinomas (9.2 months).

CONCLUSION

Intermittent arms were associated with improved efficacy/toxicity profiles; and EVE 5 mg QD and SOR 200 mg BID was defined a clinically feasible dose. Strong signs of efficacy were found in cholangiocarcinomas and gynecologic carcinomas.

TRIAL REGISTRATION

ClinicalTrials.gov Identifier: NCT01932177.

Combined KRASG12C and SOS1 inhibition enhances and extends the anti-tumor response in KRASG12C-driven cancers by addressing intrinsic and acquired resistance

Efforts to improve the anti-tumor response to KRASG12C targeted therapy have benefited from leveraging combination approaches. Here, we compare the anti-tumor response induced by the SOS1-KRAS interaction inhibitor, BI-3406, combined with a KRASG12C inhibitor (KRASG12Ci) to those induced by KRASG12Ci alone or combined with SHP2 or EGFR inhibitors. In lung cancer and colorectal cancer (CRC) models, BI-3406 plus KRASG12Ci induces an anti-tumor response stronger than that observed with KRASG12Ci alone and comparable to those by the other combinations. This enhanced anti-tumor response is associated with a stronger and extended suppression of RAS-MAPK signaling. Importantly, BI-3406 plus KRASG12Ci treatment delays the emergence of acquired adagrasib resistance in both CRC and lung cancer models and is associated with re-establishment of anti-proliferative activity in KRASG12Ci-resistant CRC models. Our findings position KRASG12C plus SOS1 inhibition therapy as a promising strategy for treating both KRASG12C-mutated tumors as well as for addressing acquired resistance to KRASG12Ci.

Integrated analysis of canine soft tissue sarcomas identifies recurrent mutations in TP53, KMT genes and PDGFB fusions

Canine soft tissue sarcomas (STS) are a heterogenous group of malignant tumors arising from mesenchymal cells of soft tissues. This simplified collective of tumors most commonly arise from subcutaneous tissues, are treated similar clinically, and conventionally exclude other sarcomas with more definitive anatomical, histological, or biological features. Histologically, canine STS sub-types are difficult to discern at the light microscopic level due to their overlapping features. Thus, genomic, and transcriptomic profiling of canine STS may prove valuable in differentiating the diverse sub-types of mesenchymal neoplasms within this group. To this purpose we sought to characterize the transcript expression and genomic mutation profiles of canine STS. To delineate transcriptomic sub-types, hierarchical clustering was used to identify 4 groups with district expression profiles. Using the RNAseq data, we identified three samples carrying driver fusions of platelet derived growth factor B ( PDGFB ) and collagen genes. Sensitivity to imatinib was evaluated in a canine STS cell line also bearing a PDGFB fusion. Using whole exome sequencing, recurrent driver variants were identified in the cancer genes KMT2D (21% of the samples) and TP53 (21%) along with copy number losses of RB1 and CDKN2A. Gene amplifications and resulting transcript increases were identified in genes on chromosomes 13, 14, and 36. A subset of STS was identified with high T-cell infiltration. This multi-omics approach has defined canine STS sub-types at a molecular level for comparison to their human counterparts, to improve diagnosis, and may provide additional targets for therapy.

High-plex imaging of RNA and proteins at subcellular resolution in fixed tissue by spatial molecular imaging

Resolving the spatial distribution of RNA and protein in tissues at subcellular resolution is a challenge in the field of spatial biology. We describe spatial molecular imaging, a system that measures RNAs and proteins in intact biological samples at subcellular resolution by performing multiple cycles of nucleic acid hybridization of fluorescent molecular barcodes. We demonstrate that spatial molecular imaging has high sensitivity (one or two copies per cell) and very low error rate (0.0092 false calls per cell) and background (~0.04 counts per cell). The imaging system generates three-dimensional, super-resolution localization of analytes at ~2 million cells per sample. Cell segmentation is morphology based using antibodies, compatible with formalin-fixed, paraffin-embedded samples. We measured multiomic data (980 RNAs and 108 proteins) at subcellular resolution in formalin-fixed, paraffin-embedded tissues (nonsmall cell lung and breast cancer) and identified >18 distinct cell types, ten unique tumor microenvironments and 100 pairwise ligand-receptor interactions. Data on >800,000 single cells and ~260 million transcripts can be accessed at http://nanostring.com/CosMx-dataset .

Extensive preclinical validation of combined RMC-4550 and LY3214996 supports clinical investigation for KRAS mutant pancreatic cancer

Over 90% of pancreatic cancers present mutations in KRAS, one of the most common oncogenic drivers overall. Currently, most KRAS mutant isoforms cannot be targeted directly. Moreover, targeting single RAS downstream effectors induces adaptive resistance mechanisms. We report here on the combined inhibition of SHP2, upstream of KRAS, using the allosteric inhibitor RMC-4550 and of ERK, downstream of KRAS, using LY3214996. This combination shows synergistic anti-cancer activity in vitro, superior disruption of the MAPK pathway, and increased apoptosis induction compared with single-agent treatments. In vivo, we demonstrate good tolerability and efficacy of the combination, with significant tumor regression in multiple pancreatic ductal adenocarcinoma (PDAC) mouse models. Finally, we show evidence that ¹⁸F-fluorodeoxyglucose (FDG) positron emission tomography (PET) can be used to assess early drug responses in animal models. Based on these results, we will investigate this drug combination in the SHP2 and ERK inhibition in pancreatic cancer (SHERPA; ClinicalTrials.gov: NCT04916236) clinical trial, enrolling patients with KRAS-mutant PDAC.

ERBB2 as a prognostic biomarker correlates with immune infiltrates in papillary thyroid cancer

Epidermal growth factor receptor 2 (ERBB2) is commonly over-expressed in advanced or metastatic tissues of papillary thyroid cancer (PTC) with poor prognosis, while it remains unknown whether ERBB2 plays a role in the progression of PTC. Thus, we analyzed the data derived from online repositories, including TCGA, KEGG, GO, GeneMANIA, and STRING, to explore the relationship between ERBB2 expression and prognosis, tumor phenotypes of interest, and immune infiltrates in PTC. Compared to normal thyroid tissue, ERBB2 was up-regulated in PTC samples (p < 0.001); In comparison with the group with low expression of ERBB2, the group with high expression of ERBB2 had poorer progression-free interval in stage III/IV patients (p = 0.008) and patients aged >45 years (p = 0.019). The up-regulated ERBB2 was associated with iodine metabolism dysfunction, proliferation, metastasis, angiogenesis, and drug resistance. The expression of ERBB2 negatively correlated with enrichment scores of B cells (r = −0.176, p < 0.001), CD8⁺ T cells (r = −0.160, p < 0.001), cytotoxic cells (r = −0.219, p < 0.001), NK CD56dim cells (r = −0.218, p < 0.001), plasmacytoid dendritic cells (r = −0.267, p < 0.001), T cells (r = −0.164, p < 0.001), T follicular helper cells (r = −0.111, p = 0.012), gamma delta T cells (r = −0.105, p = 0.017), and regulatory T cells (r = −0.125, p = 0.005). In conclusion, ERBB2 may serve as a prognostic biomarker and an immunotherapeutic target in PTC, deserving further exploration.

Optimized Transcriptional Signature for Evaluation of MEK/ERK Pathway Baseline Activity and Long-Term Modulations in Ovarian Cancer

Ovarian cancer is the most aggressive and lethal of all gynecologic malignancies. The high activity of the MEK/ERK signaling pathway is tightly associated with tumor growth, high recurrence rate, and treatment resistance. Several transcriptional signatures were proposed recently for evaluation of MEK/ERK activity in tumor tissue. In the present study, we validated the performance of a robust multi-cancer MPAS 10-gene signature in various experimental models and publicly available sets of ovarian cancer samples. Expression of four MPAS genes (PHLDA1, DUSP4, EPHA2, and SPRY4) displayed reproducible responses to MEK/ERK activity modulations across several experimental models in vitro and in vivo. Levels of PHLDA1, DUSP4, and EPHA2 expression were also significantly associated with baseline levels of MEK/ERK pathway activity in multiple human ovarian cancer cell lines and ovarian cancer patient samples available from the TCGA database. Initial platinum therapy resistance and advanced age at diagnosis were independently associated with poor overall patient survival. Taken together, our results demonstrate that the performance of transcriptional signatures is significantly affected by tissue specificity and aspects of particular experimental models. We therefore propose that gene expression signatures derived from comprehensive multi-cancer studies should be always validated for each cancer type.

Hyper-active RAS/MAPK introduces cancer-specific mitotic vulnerabilities

Aneuploidy, the incorrect number of whole chromosomes, is a common feature of tumors that contributes to their initiation and evolution. Preventing aneuploidy requires properly functioning kinetochores, which are large protein complexes assembled on centromeric DNA that link mitotic chromosomes to dynamic spindle microtubules and facilitate chromosome segregation. The kinetochore leverages at least two mechanisms to prevent aneuploidy: error correction and the spindle assembly checkpoint (SAC). BubR1, a factor involved in both processes, was identified as a cancer dependency and therapeutic target in multiple tumor types; however, it remains unclear what specific oncogenic pressures drive this enhanced dependency on BubR1 and whether it arises from BubR1's regulation of the SAC or error-correction pathways. Here, we use a genetically controlled transformation model and glioblastoma tumor isolates to show that constitutive signaling by RAS or MAPK is necessary for cancer-specific BubR1 vulnerability. The MAPK pathway enzymatically hyperstimulates a network of kinetochore kinases that compromises chromosome segregation, rendering cells more dependent on two BubR1 activities: counteracting excessive kinetochore-microtubule turnover for error correction and maintaining the SAC. This work expands our understanding of how chromosome segregation adapts to different cellular states and reveals an oncogenic trigger of a cancer-specific defect.

The kinase PLK1 promotes the development of <i>Kras</i>/<i>Tp53</i>-mutant lung adenocarcinoma through transcriptional activation of the receptor RET

Increased abundance of polo-like kinase 1 (PLK1) is observed in various tumor types, particularly in lung adenocarcinoma (LUAD). Here, we found that PLK1 accelerated the progression of LUAD through a mechanism that was independent of its role in mediating mitotic cell division. Analysis of human tumor databases revealed that increased PLK1 abundance in LUAD correlated with mutations in KRAS and p53, with tumor stage, and with reduced survival in patients. In a mouse model of KRAS^G12D-driven, p53-deficient LUAD, PLK1 overexpression increased tumor burden, decreased tumor cell differentiation, and reduced animal survival. PLK1 overexpression in cultured cells and mice indirectly increased the expression of the gene encoding the receptor tyrosine kinase RET by phosphorylating the transcription factor TTF-1. Signaling by RET and mutant KRAS in these tumors converged to activate the mitogen-activated protein kinase (MAPK) pathway. Pharmacological inhibition of the MAPK pathway kinase MEK combined with inhibition of either RET or PLK1 markedly suppressed tumor growth. Our findings show that PLK1 can amplify MAPK signaling and reveal a potential target for stemming progression in lung cancers with high PLK1 abundance.

Sensitivity towards HDAC inhibition is associated with RTK/MAPK pathway activation in gastric cancer

Gastric cancer ranks the fifth most common and third leading cause of cancer-related deaths worldwide. Alterations in the RTK/MAPK, WNT, cell adhesion, TP53, TGFβ, NOTCH, and NFκB signaling pathways could be identified as main oncogenic drivers. A combination of altered pathways can be associated with molecular subtypes of gastric cancer. In order to generate model systems to study the impact of different pathway alterations in a defined genetic background, we generated three murine organoid models: a RAS-activated (Kras^G12D , Tp53^R172H ), a WNT-activated (Apc^fl/fl , Tp53^R172H ), and a diffuse (Cdh1^fl/fl , Apc^fl/fl ) model. These organoid models were morphologically and phenotypically diverse, differed in proteome expression signatures and possessed individual drug sensitivities. A differential vulnerability to RTK/MAPK pathway interference based on the different mitogenic drivers and according to the level of dependence on the pathway could be uncovered. Furthermore, an association between RTK/MAPK pathway activity and susceptibility to HDAC inhibition was observed. This finding was further validated in patient-derived organoids from gastric adenocarcinoma, thus identifying a novel treatment approach for RTK/MAPK pathway altered gastric cancer patients.

Identification of cell type-specific correlations between ERK activity and cell viability upon treatment with ERK1/2 inhibitors

Increased MAPK signaling is a hallmark of various cancers and is a central regulator of cell survival. Direct ERK1/2 inhibition is considered a promising approach to avoid ERK1/2 reactivation caused by upstream kinases BRAF, MEK1/2, and KRAS, as well as by receptor tyrosine kinase inhibitors, but the dynamics and selectivity of ERK1/2 inhibitors are much less studied compared with BRAF or MEK inhibitors. Using ERK1/2 and downstream kinase ELK1 reporter cell lines of lung cancer (H1299; NRAS^Q61K), colon cancer (HCT-116; KRAS^G13D), neuroblastoma (SH-SY5Y), and leukemia (U937), we examined the relationship between ERK inhibition and drug-induced toxicity for five ERK inhibitors: SCH772984, ravoxertinib, LY3214996, ulixertinib, and VX-11e, as well as one MEK inhibitor, PD0325901. Comparing cell viability and ERK inhibition revealed different ERK dependencies for these cell lines. We identify several drugs, such as SCH772984 and VX-11e, which induce excessive toxicity not directly related to ERK1/2 inhibition in specific cell lines. We also show that PD0325901, LY3214996, and ulixertinib are prone to ERK1/2 reactivation over time. We distinguished two types of ERK1/2 reactivation: the first could be reversed by adding a fresh dose of inhibitors, while the second persists even after additional treatments. We also showed that cells that became resistant to the MEK1/2 inhibitor PD0325901 due to ERK1/2 reactivation remained sensitive to ERK1/2 inhibitor ulixertinib. Our data indicate that correlation of ERK inhibition with drug-induced toxicity in multiple cell lines may help to find more selective and effective ERK1/2 inhibitors.

Molecular Characterization of KRAS Wild-type Tumors in Patients with Pancreatic Adenocarcinoma

PURPOSE

KRAS mutation (MT) is a major oncogenic driver in pancreatic ductal adenocarcinoma (PDAC). A small subset of PDACs harbor KRAS wild-type (WT). We aim to characterize the molecular profiles of KRAS WT PDAC to uncover new pathogenic drivers and offer targeted treatments.

EXPERIMENTAL DESIGN

Tumor tissue obtained from surgical or biopsy material was subjected to next-generation DNA/RNA sequencing, microsatellite instability (MSI) and mismatch repair status determination.

RESULTS

Of the 2,483 patients (male 53.7%, median age 66 years) studied, 266 tumors (10.7%) were KRAS WT. The most frequently mutated gene in KRAS WT PDAC was TP53 (44.5%), followed by BRAF (13.0%). Multiple mutations within the DNA-damage repair (BRCA2, ATM, BAP1, RAD50, FANCE, PALB2), chromatin remodeling (ARID1A, PBRM1, ARID2, KMT2D, KMT2C, SMARCA4, SETD2), and cell-cycle control pathways (CDKN2A, CCND1, CCNE1) were detected frequently. There was no statistically significant difference in PD-L1 expression between KRAS WT (15.8%) and MT (17%) tumors. However, KRAS WT PDAC were more likely to be MSI-high (4.7% vs. 0.7%; P < 0.05), tumor mutational burden-high (4.5% vs. 1%; P < 0.05), and exhibit increased infiltration of CD8+ T cells, natural killer cells, and myeloid dendritic cells. KRAS WT PDACs exhibited gene fusions of BRAF (6.6%), FGFR2 (5.2%), ALK (2.6%), RET (1.3%), and NRG1 (1.3%), as well as amplification of FGF3 (3%), ERBB2 (2.2%), FGFR3 (1.8%), NTRK (1.8%), and MET (1.3%). Real-world evidence reveals a survival advantage of KRAS WT patients in overall cohorts as well as in patients treated with gemcitabine/nab-paclitaxel or 5-FU/oxaliplatin.

CONCLUSIONS

KRAS WT PDAC represents 10.7% of PDAC and is enriched with targetable alterations, including immuno-oncologic markers. Identification of KRAS WT patients in clinical practice may expand therapeutic options in a clinically meaningful manner.

MITF deficiency accelerates GNAQ-driven uveal melanoma

Cutaneous melanoma (CM) and uveal melanoma (UM) both originate from the melanocytic lineage but are primarily driven by distinct oncogenic drivers, BRAF/NRAS or GNAQ/GNA11, respectively. The melanocytic master transcriptional regulator, MITF, is essential for both CM development and maintenance, but its role in UM is largely unexplored. Here, we use zebrafish models to dissect the key UM oncogenic signaling events and establish the role of MITF in UM tumors. Using a melanocytic lineage expression system, we showed that patient-derived mutations of GNAQ (GNAQQ209L) or its upstream CYSLTR2 receptor (CYSLTR2L129Q) both drive UM when combined with a cooperating mutation, tp53M214K/M214K. The tumor-initiating potential of the major GNAQ/11 effector pathways, YAP, and phospholipase C-β (PLCβ)–ERK was also investigated in this system and thus showed that while activated YAP (YAPAA) induced UM with high potency, the patient-derived PLCβ4 mutation (PLCB4D630Y) very rarely yielded UM tumors in the tp53M214K/M214K context. Remarkably, mitfa deficiency was profoundly UM promoting, dramatically accelerating the onset and progression of tumors induced by Tg(mitfa:GNAQQ209L);tp53M214K/M214K or Tg(mitfa:CYSLTR2L129Q);tp53M214K/M214K. Moreover, mitfa loss was sufficient to cooperate with GNAQQ209L to drive tp53–wild type UM development and allowed Tg(mitfa:PLCB4D630Y);tp53M214K/M214K melanocyte lineage cells to readily form tumors. Notably, all of the mitfa−/− UM tumors, including those arising in Tg(mitfa:PLCB4D630Y);tp53M214K/M214K;mitfa−/− zebrafish, displayed nuclear YAP while lacking hyperactive ERK indicative of PLCβ signaling. Collectively, these data show that YAP signaling is the major mediator of UM and that MITF acts as a bona fide tumor suppressor in UM in direct opposition to its essential role in CM.

Endoluminal Biopsy for Molecular Profiling of Human Brain Vascular Malformations

BACKGROUND AND OBJECTIVES

Ras-mitogen-activated protein kinase (MAPK) signaling abnormalities occur in most brain arteriovenous malformations (bAVMs). No means exist to molecularly profile bAVMs without open surgery, limiting precision medicine approaches to treatment. Here, we report use of endoluminal biopsy of the vessel lumen of bAVMs to characterize gene expression and blood flow-mediated transcriptional changes in living patients.

METHODS

Endoluminal biopsy and computational fluid dynamic modeling (CFD) were performed in adults with unruptured AVMs with cerebral angiography. Each patient underwent surgical resection and cell sampling from a contiguous arterial segment. Fluorescence-assisted cell sorting enriched endothelial cells, which were sequenced on an Illumina HiSeq 4000 sequencer. Gene expression was quantified with RNA sequencing (RNAseq). Differential gene expression, ontology, and correlative analyses were performed. Results were validated with quantitative reverse transcription PCR (RT-qPCR).

RESULTS

Endoluminal biopsy was successful in 4 patients without complication. Endoluminal biopsy yielded 269.0 ± 79.9 cells per biopsy (control 309.2 ± 86.6 cells, bAVM 228.8 ± 133.4 cells). RNAseq identified 106 differentially expressed genes (DEGs) in bAVMs (false discovery rate ≤0.05). DEGs were enriched for bAVM pathogenic cascades, including Ras-MAPK signaling (p < 0.05), and confirmed with RT-qPCR and a panel predictive of MAPK/extracellular signal-regulated kinase inhibitor response. Compared to patient-matched surgically excised tissues, endoluminal biopsy detected 83.3% of genes, and genome-wide expression strongly correlated (Pearson r = 0.77). Wall shear stress measured by CFD correlated with inflammatory pathway upregulation. Comparison of pre-embolization and postembolization samples confirmed flow-mediated gene expression changes.

DISCUSSION

Endoluminal biopsy allows molecular profiling of bAVMs in living patients. Gene expression profiles are similar to those of tissues acquired with open surgery and identify potentially targetable Ras-MAPK signaling abnormalities in bAVMs. Integration with CFD allows determination of flow-mediated transcriptomic alterations. Endoluminal biopsy may help facilitate trials of precision medicine approaches to bAVMs in humans.

Transcriptional Profiling of Malignant Melanoma Reveals Novel and Potentially Targetable Gene Fusions

Simple Summary

Malignant melanoma is a complex disease that is estimated to claim over 7000 lives in the United States in 2021. Although recent advances in genomic technology have helped with the identification of driver variants, molecular studies and clinical trials have often focused on prevalent alterations, such as the BRAF-V600E mutation. With the inclusion of whole transcriptome sequencing, molecular profiling of melanomas has identified gene fusions and revealed gene expression profiles that are consistent with the activation of signaling pathways by common driver mutations. Patients harboring such fusions may benefit from currently approved targeted therapies and should be considered in the design of future clinical trials to further personalize treatments for patients with malignant melanoma.

Abstract

Invasive melanoma is the deadliest type of skin cancer, with 101,110 expected cases to be diagnosed in 2021. Recurrent BRAF and NRAS mutations are well documented in melanoma. Biologic implications of gene fusions and the efficacy of therapeutically targeting them remains unknown. Retrospective review of patient samples that underwent next-generation sequencing of the exons of 592 cancer-relevant genes and whole transcriptome sequencing for the detection of gene fusion events and gene expression profiling. Expression of PDL1 and ERK1/2 was assessed by immunohistochemistry (IHC). There were 33 (2.6%) cases with oncogenic fusions (14 novel), involving BRAF, RAF1, PRKCA, TERT, AXL, and FGFR3. MAPK pathway-associated genes were over-expressed in BRAF and RAF1 fusion-positive tumors in absence of other driver alterations. Increased expression in tumors with PRKCA and TERT fusions was concurrent with MAPK pathway alterations. For a subset of samples with available tissue, increased phosphorylation of ERK1/2 was observed in BRAF, RAF1, and PRKCA fusion-positive tumors. Oncogenic gene fusions are associated with transcriptional activation of the MAPK pathway, suggesting they could be therapeutic targets with available inhibitors. Additional analyses to fully characterize the oncogenic effects of these fusions may support biomarker driven clinical trials.

Characterizing the molecular and immune landscape of canine bladder cancer

Transitional cell carcinoma (TCC), also known as urothelial carcinoma, is the most common bladder cancer in humans and dogs. Approximately one-quarter of human TCCs are muscle-invasive and associated with a high risk of death from metastasis. Canine TCC (cTCC) tumours are typically high-grade and muscle-invasive. Shared similarities in risk factors, histopathology, and clinical presentation suggest that cTCC may serve as a model for the assessment of novel therapeutics that may inform therapies for human muscle-invasive TCC. The goal of this study was to characterize cTCC at the molecular level to identify drivers of oncogenesis and druggable targets. We performed whole exome sequencing (WES) of 11 cTCC tumours and three matched normal samples, identifying 583 variants in protein-coding genes. The most common variant was a V-to-E missense mutation in BRAF, identified in 4 out of 11 samples (36%) via WES. Sanger sequencing identified BRAF variants in 8 out of the same 11 cTCC samples, as well as in 22 out of 32 formalin-fixed paraffin embedded (FFPE) cTCC samples, suggesting an overall prevalence of 70%. RNA-Seq was performed to compare the gene expression profiles of cTCC tumours to normal bladder tissue. cTCC tumours exhibited up-regulation of genes involved in the cell cycle, DNA repair, and antiviral immunity. We also analysed the immune landscape of cTCC using immune gene signatures and immunohistochemical analysis. A subset of tumours had characteristics of a hot tumour microenvironment and exhibited high expression of signatures associated with complete response to PD-1/PD-L1 blockade in human bladder cancer.

Activation of a cGAS-STING-mediated immune response predicts response to neoadjuvant chemotherapy in early breast cancer

BACKGROUND

The DNA-damage immune-response (DDIR) signature is an immune-driven gene expression signature retrospectively validated as predicting response to anthracycline-based therapy. This feasibility study prospectively evaluates the use of this assay to predict neoadjuvant chemotherapy response in early breast cancer.

METHODS

This feasibility study assessed the integration of a novel biomarker into clinical workflows. Tumour samples were collected from patients receiving standard of care neoadjuvant chemotherapy (FEC + /-taxane and anti-HER2 therapy as appropriate) at baseline, mid- and post-chemotherapy. Baseline DDIR signature scores were correlated with pathological treatment response. RNA sequencing was used to assess chemotherapy/response-related changes in biologically linked gene signatures.

RESULTS

DDIR signature reports were available within 14 days for 97.8% of 46 patients (13 TNBC, 16 HER2 + ve, 27 ER + HER2-ve). Positive scores predicted response to treatment (odds ratio 4.67 for RCB 0-1 disease (95% CI 1.13-15.09, P = 0.032)). DDIR positivity correlated with immune infiltration and upregulated immune-checkpoint gene expression.

CONCLUSIONS

This study validates the DDIR signature as predictive of response to neoadjuvant chemotherapy which can be integrated into clinical workflows, potentially identifying a subgroup with high sensitivity to anthracycline chemotherapy. Transcriptomic data suggest induction with anthracycline-containing regimens in immune restricted, "cold" tumours may be effective for immune priming.

TRIAL REGISTRATION

Not applicable (non-interventional study). CRUK Internal Database Number 14232.

CD97 is associated with mitogenic pathway activation, metabolic reprogramming, and immune microenvironment changes in glioblastoma

Glioblastoma (GBM) is the most common primary brain tumor with a median survival under two years. Using in silico and in vitro techniques, we demonstrate heterogeneous expression of CD97, a leukocyte adhesion marker, in human GBM. Beyond its previous demonstrated role in tumor invasion, we show that CD97 is also associated with upregulation of the mitogen-activated protein kinase/extracellular signal-regulated kinase (MAPK/Erk) and phosphatidylinositol 3-kinase/protein kinase B (PI3K/Akt) pathways in GBM. While CD97 knockout decreased Akt activation, CD97 targeting did not alter MAPK/Erk activation, did not slow GBM cell proliferation in culture, and increased levels of glycolytic and oxidative phosphorylation metabolites. Treatment with a soluble CD97 inhibitor did not alter activation of the MAPK/Erk and PI3K/Akt pathways. Tumors with high CD97 expression were associated with immune microenvironment changes including increased naïve macrophages, regulatory T cells, and resting natural killer (NK) cells. These data suggest that, while CD97 expression is associated with conflicting effects on tumor cell proliferative and metabolic pathways that overall do not affect tumor cell proliferation, CD97 exerts pro-tumoral effects on the tumor immune microenvironment, which along with the pro-invasive effects of CD97 we previously demonstrated, provides impetus to continue exploring CD97 as a therapeutic target in GBM.