Vanderwalde, Ari M.

Abstract Natural killer (NK) cells represent key effectors of antitumor immunity, yet emerging evidence highlights populations with distinct roles in cancer. Despite such expanded diversity within the NK cell repertoire, we lack an understanding of how this heterogeneity affects immune responses and downstream clinical outcomes. Using single-cell RNA sequencing, we systematically profiled NK cells across cancer and uncovered a dichotomous phenotypic and functional landscape of tumor-infiltrating NK cells shaped by opposing intrinsic signaling programs that drive the expression of IFNG or TGFB1. These divergent programs are associated with distinct transcription factor circuits that integrate cues within the tumor microenvironment and skew NK cells toward proinflammatory or suppressive functions. We found that the capacity for NK cells to engage in either functional direction is intrinsically linked to their phenotypic identity. Canonical NK cells recruited from circulation predominantly directed suppressive TGFB1 signals toward effector CD8+ T cells in tumors. Of note, these subsets exhibited higher TGFB1 expression than intratumoral myeloid cells across tumor types. In contrast, a tissue-resident (TR) adaptive subset exhibited exclusively proinflammatory IFNG-driven profiles and was associated with prolonged survival in both primary and metastatic tumor settings. Moreover, these TR adaptive NK cells, but not other subsets, were linked to response to immune checkpoint blockade. Collectively, our study reveals a previously unrecognized regulatory axis in NK cells that shapes NK cell diversity and augments broader antitumor immune responses.

e14555 Background: Natural killer (NK) cells are innate lymphocytes that engage in direct cytotoxic responses against malignant cells and promote pro-inflammatory tumor microenvironments. While studies have predominantly focused on conventional NK cell subsets, recent findings uncovered adaptive and tissue-resident (TR) subsets with unique functional properties. The influence of distinct NK cell subsets on tumor immunobiology and clinical outcomes remains unknown. Methods: Single-cell RNA-sequencing (scRNA-seq) of primary, therapy-naïve tumors from 276 patients across 14 cancer types, as well as 39 patients treated with immune checkpoint blockade (ICB), was accessed through publicly available datasets. We employed systematic computational analyses to phenotype and examine intra-tumoral NK cells. Gene signatures for NK cell subsets were applied to bulk RNA-sequencing from 212,105 tumor samples to delineate their associations with clinical outcomes using the Caris CODEai database. Results: Six main subsets of NK cells were identified in solid tumors, including conventional CD56 dim and CD56 bright NK cells, alongside adaptive, TR adaptive, TR CD56 bright , and TR CD56 bright /adaptive subsets. Across cancer, TR adaptive NK cells emerged as the most prognostic subset, with high infiltration conferring significantly longer overall survival (HR[95% CI] = 0.68[0.66-0.69], p < 0.001) than conventional CD56 dim and CD56 bright subsets. This was also observed in cohorts of head and neck, non-small cell lung, and esophageal cancers, particularly those treated with ICB. Notably, elevated TR adaptive NK cell infiltration was more prognostic than high tumor mutational burden in these settings. This was most apparent in ICB-treated head and neck cancers (HR = 0.48[0.41-0.57], p < 0.001 vs HR = 0.90[0.76-1.04], p = 0.16). Additionally, TR adaptive NK cells were strongly correlated to the infiltration of CD8 + T cells and B cells in tumors. Underlying these observations, our scRNA-seq analyses demonstrated that TR adaptive NK cells exhibit the highest immunomodulatory activity of all subsets, expressing high levels of pro-inflammatory cytokines ( CCL3 , CXCL13 , and IFNG ) and low levels of immunosuppressive cytokines ( TGFB1 ). Mechanistically, this activity was linked to the upregulation of transcription factors associated with NK cell activation, function, and persistence and downregulation of genes related to stress response and exhaustion. Interestingly, scRNA-seq from paired pre- and post-ICB-treated tumors revealed a unique expansion of cytotoxic GZMB + TR adaptive NK cells after treatment. Conclusions: TR adaptive NK cells represent a novel subset of NK cells that may drive anti-tumor immune responses, particularly in the context of ICB treatment. Features of this unique population may be harnessed to design more effective NK cell-based therapies against solid tumors.

Abstract Pulmonary large cell neuroendocrine carcinoma (LCNEC) is a rare, aggressive lung tumor marked by significant molecular heterogeneity. In a study of 590 patients across two independent cohorts, we observed comparable overall survival across treatment regimens (chemotherapy, chemoimmunotherapy, immunotherapy) without unexpected adverse events. Genomic analysis identified distinct NSCLC-like (KEAP1, KRAS, STK11 mutations) and SCLC-like (RB1, TP53 mutations) LCNEC subtypes, with 80% aligning with SCLC transcriptional profiles. Serial sampling revealed stable mutational but shifting transcriptomic landscapes over time. NSCLC-like LCNECs showed elevated FGL-1 (a LAG-3 ligand) and SPINK1 expression, while SCLC-like subtypes expressed higher levels of DLL3. Immunofluorescence confirmed FGL-1 in NSCLC-like LCNECs, and H&E slide analyses indicated fewer tumor-infiltrating lymphocytes in LCNECs versus other lung cancers. These findings highlight LCNEC’s distinct immunogenomic profile, supporting future investigations into LAG-3, SPINK1, and DLL3-targeted therapies.

Abstract B7-H3 (CD276) is a transmembrane glycoprotein of the B7 immune checkpoint superfamily that has emerged as a promising therapeutic target. To better understand the applicability of B7-H3–directed therapies, we analyzed 156,791 samples comprising 50 cancer types to interrogate the clinical, genomic, transcriptomic, and immunologic correlates of B7-H3 mRNA expression. DNA (592-gene/whole-exome) and RNA (whole-transcriptome) sequencing was performed from samples submitted to Caris Life Sciences. B7-H3 high versus low expression was based on top and bottom quartiles for each cancer type. Patients’ overall survival was determined from insurance claims data. Pathway analysis was performed using gene set enrichment analyses. Immune cell fractions were inferred using quanTIseq. B7-H3 is expressed across several human malignancies including prostate, pancreatic, ovarian, and lung cancers. High B7-H3 expression is associated with differences in overall survival, possibly indicating a prognostic role of B7-H3 for some cancers. When examining molecular features across all cancer types, we did not identify recurrent associations between B7-H3 expression and genetic alterations in TP53, RB1, and KRAS. However, we find consistent enrichment of epithelial-to-mesenchymal transition, Wnt, TGFβ, and Notch signaling pathways. In addition, tumors with high B7-H3 expression are associated with greater proportions of M1 macrophages, but lower fractions of CD8+ T cells. We have begun to define the genomic, transcriptomic, clinical, and immunologic features associated with B7-H3 expression in 50 cancer types. We report novel clinical and molecular features of B7-H3–high tumors which may inform how current B7-H3 therapeutics should be deployed and prioritized. Significance: B7-H3–targeting therapeutics have shown promising results in initial clinical trials. In this pan-cancer analysis of B7-H3 mRNA expression, we found that B7-H3 exhibits robust expression in many common cancer types. These results may inform further development of B7-H3–targeting therapeutics and may guide clinical decisions for patients with limited treatment options.

9069 Background: With the emergence of effective targeted-therapies for KRAS G12C, and the development of promising agents for KRAS G12V and G12D among others, the identity of unique KRAS mutations in non-small cell lung cancer (NSCLC) has become increasingly relevant. Acquired KRAS mutations are a known mechanism of resistance observed across driver mutation-positive (DM+) NSCLC. The incidence and diversity of these acquired alterations, and whether they differ from those observed in de novo KRAS-mutated NSCLC, is unknown. Methods: NSCLC molecular profiles were obtained using next-generation sequencing (Caris Life Sciences) with paired whole-transcriptomic sequencing (Illumina NovaSeq) and immunohistochemistry (Caris Life Sciences). Demographic data were abstracted from medical records. KRAS mutation diversity was defined for de novo KRAS-mutated ( KRASmt) NSCLC ( KRAS as only identified driver – KRASmt_dn) and DM+ NSCLC with acquired KRAS mutations (concurrent KRASmt with other known driver – KRASmt_acq). Fisher’s exact test was used to compare the distribution of unique KRAS mutations between the groups. Due to the unique biology of NSCLC with class II/III BRAF mutations, this subset was removed from the KRASmt_acq subgroup for the final analysis. Results: A total of 5932 KRASmt NSCLC samples were identified, among which 5879 were KRASmt_dn and 53 were KRASmt_acq (see table). The distribution of unique KRAS mutations was not significantly different between groups (p = 0.14). Within the KRASmt_dn group, KRAS G12C was most common (40%), followed by G12V (19%), G12D (14%), G12A (7%) and Q61H/G13C (4% each). In the KRASmt_acq group, KRAS G12C was most common (30%), followed by G12D (19%), G12V (17%), Q61H (11%), G12A (7%), and G13C (6%). The most common observed driver mutations in the KRASmt_acq group were EGFR (38%) and MET (30%). In the EGFR group, KRAS G12D (25%), Q61H/G12C (20% each), G12V/G13C (10% each), and G12F/A/R (5% each) were observed. In the MET group, KRAS G12C (25%), G12D/V/A (19% each), Q61H (12%) and G12S (6%) were seen. Among 7 patients with a documented history of smoking in the KRASmt_acq group, KRAS G12D and Q61H were most common (29% each), followed by G12V/F and G13C (14% each). Conclusions: While the distribution of unique KRAS mutations did not differ significantly between KRASmt_dn and KRASmt_acq groups, acquired KRAS mutations were seen across DM+ NSCLC subsets, among which the frequency of observed mutations appeared to vary. The functional and immunological significance of these mutations, and their impact on clinical outcomes, warrants further investigation. [Table: see text]

AbstractBackground and ObjectivesPublished data comparing peritoneal metastases from appendiceal cancers (pAC) and colorectal cancers (pCRC) remain sparse. We compared pAC and pCRC using comprehensive tumor profiling (CTP).MethodsCTP was performed, including next‐generation sequencing and analysis of copy number variation (CNV), microsatellite instability (MSI) and tumor mutational burden (TMB).ResultsOne hundred thirty‐six pAC and 348 pCRC samples underwent CTP. The cohorts' age and gender were similar. pCRC demonstrated increased pathogenic variants (PATHs) in APC (48% vs. 3%, p < 0.01), ARID1A (12% vs. 2%, p < 0.01), BRAF (12% vs. 2%, p < 0.01), FBXW7 (7% vs. 2%, p < 0.01), KRAS (52% vs. 41%, p < 0.05), PIK3CA (15% vs. 2%, p < 0.01), and TP53 (53% vs. 23%, p < 0.01), and decreased PATHs in GNAS (8% vs. 31%, p < 0.01). There was no difference in CNV, fusion rate, or MSI. Median TMB was higher in pCRC (5.8 vs. 5.0 mutations per megabase, p = 0.0007). Rates of TMB‐high tumors were similar (pAC 2.1% vs. pCRC 9.0%, p = 0.1957). pCRC had significantly more TMB‐high tumors at lower thresholds.ConclusionsDespite a reduced overall TMB, pAC demonstrated mutations distinct from those seen in pCRC. These may serve as discrete biomarkers for future study.

Abstract Purpose: Mutations in KRAS/NRAS (RAS) predict lack of anti-EGFR efficacy in metastatic colorectal cancer (mCRC). However, it is unclear if all RAS mutations have similar impact, and atypical mutations beyond those in standard guidelines exist. Experimental Design: We reviewed 7 tissue and 1 cell-free DNA cohorts of 9,485 patients to characterize atypical RAS variants. Using an in vitro cell-based assay (functional annotation for cancer treatment), Ba/F3 transformation, and in vivo xenograft models of transduced isogenic clones, we assessed signaling changes across mutations. Results: KRAS exon 2, extended RAS, and atypical RAS mutations were noted in 37.8%, 9.5%, and 1.2% of patients, respectively. Among atypical variants, KRAS L19F, Q22K, and D33E occurred at prevalence ≥0.1%, whereas no NRAS codon 117/146 and only one NRAS codon 59 mutation was noted. Atypical RAS mutations had worse overall survival than RAS/BRAF wild-type mCRC (HR, 2.90; 95% confidence interval, 1.24–6.80; P = 0.014). We functionally characterized 114 variants with the FACT assay. All KRAS exon 2 and extended RAS mutations appeared activating. Of 57 atypical RAS variants characterized, 18 (31.6%) had signaling below wild-type, 23 (40.4%) had signaling between wild-type and activating control, and 16 (28.1%) were hyperactive beyond the activating control. Ba/F3 transformation (17/18 variants) and xenograft model (7/8 variants) validation was highly concordant with FACT results, and activating atypical variants were those that occurred at highest prevalence in clinical cohorts. Conclusions: We provide best available evidence to guide treatment when atypical RAS variants are identified. KRAS L19F, Q22K, D33E, and T50I are more prevalent than many guideline-included RAS variants and functionally relevant.

606 Background: Between 3-5% of CRC patients have BRCA1/2 pathogenic mutations. This study aims to identify associations between BRCA1 and BRCA2 mutations and clinical characteristics in CRC. Methods: A total of 6396 CRC tumor samples were tested with Next-Generation Sequencing (NGS) on a 592-gene panel, pathogenic or presumed pathogenic variants were counted as mutations (mt). Microsatellite instability (MSI) and tumor mutational burden (TMB) were tested by NGS. Statistical correlations were investigated using ANOVA, Chi-square and t-test. Results: Among tumors sampled, 53% derived from male patients and median age was 60 years. BRCA1 mt were detected in 1.1% (n = 72) of tumors, while BRCA2 in 2.8% (n = 179). BRCA1 mt were more frequent in women (W;65%) than men (M;35%) (p = 0.0019) while no relationship with sex was seen for BRCA2 mt (42% F vs. 58% M). No significant associations with age were noticed. Majority of pathogenic mt in BRCA1 (52%; n = 34) and BRCA2 (62%; n = 103) occurred in MSI-High (MSI-H) cases. MSI-H pts had more frameshift mt in both BRCA1/2 than MSS pts. MSS cases had lower rates of BRCA1 and 2 pathogenic mt (44% and 37%, respectively). Right-sided tumors were significantly associated with BRCA1 (p = 0.0056) and BRCA2 (p < 0.0001) mt in MSI-H cases only. BRCA1/2 mt were associated with higher TMB in all CRCs, including MSI-H and MSS cases (p < 0.001). POLE mt (n = 31) were associated with higher BRCA1/2 mt rates (9.6%, 55% respectively). Among MSS cases with POLE wild-type status, BRCA1 (p = 0.0269) and BRCA2 (p = 0.0151) mt were associated with high TMB and combining both BRCA1/2 mutations led to an even higher TMB (3.6%; p = 0.001). Conclusions: This is the first study to show that BRCA1/2 mutations are more frequent in MSI-H, and independently associated with higher TMB, pathogenic POLE mutations, and right-sided tumors in MSI-H CRCs. Given their relationship with TMB, the presence of BRCA1/2 mutations may be potential predictive biomarkers for checkpoint or PARP inhibitors in CRC, a finding that should be prospectively evaluated.

60 Background: DDR mutations are associated with higher neoantigen load and tumor infiltrating lymphocytes, and are a potential biomarker for immunotherapy. We investigated the association of DDR mutations in gastric (GC), esophageal (EC), and gastroesophageal junction (GEJ) cancers with MSI, PD-L1, and TMB, known predictors for immune checkpoint inhibitors. Methods: 20 DDR mutations were tested by Next-Generation Sequencing (NGS) with a 592-gene panel on a total of 1935 (709 EC; 831 GC; 355 GEJ) cancers. TMB was assessed by NGS, MSI by NGS or fragment analysis, and PD-L1 by IHC (22c3 for CPS or SP142). Results: GC had the highest DDR mutation rate compared to EC and GEJ (27% vs. 20%, p = 0.0005 and 17%, p = 0.0002, respectively). MSI-High (MSI-H) was significantly more common in the DDR mutated cohort (DDR-M) compared to non-mutated cancers (18% vs. 1%; p < 0.0001). TMB-High (≥ 10 mutations/megabase [mt/MB]) was higher in DDR-M (35% vs. 21%; p < 0.0001); in DDR-M cohort, GC had the highest TMB compared to DDR-M EC and GEJ (mean: 13.8 vs. 9.4 vs. 10 mt/MB, respectively; p < 0.0001). DDR mutations were more frequent in the PD-L1 combined positive score (CPS) ≥ 50 group than CPS 0 (42.9% vs. 24%; p = 0.037) and CPS 1-9 (42.9% vs. 20.6%; p = 0.005). ARID1A, ATRX, BRCA2, and PTEN were the most prevalent DDR mutations in MSI-H (87%, 31%, 25%, 24%, respectively); ARID1A, ATRX, BRCA2, and PTEN in TMB-High (47%, 7.7%, 6.7%, 6.8%) and ARID1A, BRCA2, RAD50, and WRN in PD-L1 high (CPS ≥ 10) (48.5% vs. 5.2% vs. 2.5% vs. 3.4%). Conclusions: MSI-H, TMB-high and high PD-L1 expression were significantly more prevalent in the DDR-M cohort compared to non-DDR-mutated cancers, most pronounced in GC. Alterations in ARID1A, ATRX, BRCA2, and PTEN were correlated with MSI-H and TMB-high while ARID1A, BRCA2, RAD50, and WRN were correlated with increased PD-L1 expression. Our findings may help identify patients for tailored immunotherapy approaches in future clinical trials.