Comprehensive mutational scanning of EGFR reveals TKI sensitivities of extracellular domain mutants.

The epidermal growth factor receptor, EGFR, is frequently activated in lung cancer and glioblastoma by genomic alterations including missense mutations. The different mutation spectra in these diseases are reflected in divergent responses to EGFR inhibition: significant patient benefit in lung cancer, but limited in glioblastoma. Here, we report a comprehensive mutational analysis of EGFR function. We perform saturation mutagenesis of EGFR and assess function of ~22,500 variants in a human EGFR-dependent lung cancer cell line. This approach reveals enrichment of erlotinib-insensitive variants of known and unknown significance in the dimerization, transmembrane, and kinase domains. Multiple EGFR extracellular domain variants, not associated with approved targeted therapies, are sensitive to afatinib and dacomitinib in vitro. Two glioblastoma patients with somatic EGFR G598V dimerization domain mutations show responses to dacomitinib treatment followed by within-pathway resistance mutation in one case. In summary, this comprehensive screen expands the landscape of functional EGFR variants and suggests broader clinical investigation of EGFR inhibition for cancers harboring extracellular domain mutations.

RAS/RAF Comutation and ERBB2 Copy Number Modulates HER2 Heterogeneity and Responsiveness to HER2-directed Therapy in Colorectal Cancer.

PURPOSE: ERBB2-amplified colorectal cancer is a distinct molecular subtype with expanding treatments. Implications of concurrent oncogenic RAS/RAF alterations are not known. EXPERIMENTAL DESIGN: Dana-Farber and Foundation Medicine Inc. Colorectal cancer cohorts with genomic profiling were used to identify ERBB2-amplified cases [Dana-Farber, n = 47/2,729 (1.7%); FMI, n = 1857/49,839 (3.7%)]. Outcomes of patients receiving HER2-directed therapies are reported (Dana-Farber, n = 9; Flatiron Health-Foundation Medicine clinicogenomic database, FH-FMI CGDB, n = 38). Multisite HER2 IHC and genomic profiling were performed to understand HER2 intratumoral and interlesional heterogeneity. The impact of concurrent RAS comutations on the effectiveness of HER2-directed therapies were studied in isogenic colorectal cancer cell lines and xenografts. RESULTS: ERBB2 amplifications are enriched in left-sided colorectal cancer. Twenty percent of ERBB2-amplified colorectal cancers have co-occurring oncogenic RAS/RAF alterations. While RAS/RAF WT colorectal cancers typically have clonal ERBB2 amplification, colorectal cancers with co-occurring RAS/RAF alterations have lower level ERRB2 amplification, higher intratumoral heterogeneity, and interlesional ERBB2 discordance. These distinct genomic patterns lead to differential responsiveness and patterns of resistance to HER2-directed therapy. ERBB2-amplified colorectal cancer with RAS/RAF alterations are resistant to trastuzumab-based combinations, such as trastuzumab/tucatinib, but retain sensitivity to trastuzumab deruxtecan in in vitro and murine models. Trastuzumab deruxtecan shows clinical efficacy in cases with high-level ERBB2-amplified RAS/RAF coaltered colorectal cancer. CONCLUSIONS: Co-occurring RAS/RAF alterations define a unique subtype of ERBB2-amplified colorectal cancer that has increased intratumoral heterogeneity, interlesional discordance, and resistance to trastuzumab-based combinations. Further examination of trastuzumab deruxtecan in this previously understudied cohort of ERBB2-amplified colorectal cancer is warranted.

XRN1 deletion induces PKR-dependent cell lethality in interferon-activated cancer cells.

Emerging data suggest that induction of viral mimicry responses through activation of double-stranded RNA (dsRNA) sensors in cancer cells is a promising therapeutic strategy. One approach to induce viral mimicry is to target molecular regulators of dsRNA sensing pathways. Here, we show that the exoribonuclease XRN1 is a negative regulator of the dsRNA sensor protein kinase R (PKR) in cancer cells with high interferon-stimulated gene expression. XRN1 deletion causes PKR pathway activation and consequent cancer cell lethality. Disruption of interferon signaling with the JAK1/2 inhibitor ruxolitinib can decrease cellular PKR levels and rescue sensitivity to XRN1 deletion. Conversely, interferon-ß stimulation can increase PKR levels and induce sensitivity to XRN1 inactivation. Lastly, XRN1 deletion causes accumulation of endogenous complementary sense/anti-sense RNAs, which may represent candidate PKR ligands. Our data demonstrate how XRN1 regulates PKR and how this interaction creates a vulnerability in cancer cells with an activated interferon cell state.

Genomic and Immunophenotypic Landscape of Acquired Resistance to PD-(L)1 Blockade in Non-Small-Cell Lung Cancer.

PURPOSE: Although immune checkpoint inhibitors (ICI) have extended survival in patients with non-small-cell lung cancer (NSCLC), acquired resistance (AR) to ICI frequently develops after an initial benefit. However, the mechanisms of AR to ICI in NSCLC are largely unknown. METHODS: Comprehensive tumor genomic profiling, machine learning-based assessment of tumor-infiltrating lymphocytes, multiplexed immunofluorescence, and/or HLA-I immunohistochemistry (IHC) were performed on matched pre- and post-ICI tumor biopsies from patients with NSCLC treated with ICI at the Dana-Farber Cancer Institute who developed AR to ICI. Two additional cohorts of patients with intervening chemotherapy or targeted therapies between biopsies were included as controls. RESULTS: We performed comprehensive genomic profiling and immunophenotypic characterization on samples from 82 patients with NSCLC and matched pre- and post-ICI biopsies and compared findings with a control cohort of patients with non-ICI intervening therapies between biopsies (chemotherapy, N = 32; targeted therapies, N = 89; both, N = 17). Putative resistance mutations were identified in 27.8% of immunotherapy-treated cases and included acquired loss-of-function mutations in STK11, B2M, APC, MTOR, KEAP1, and JAK1/2; these acquired alterations were not observed in the control groups. Immunophenotyping of matched pre- and post-ICI samples demonstrated significant decreases in intratumoral lymphocytes, CD3e+ and CD8a+ T cells, and PD-L1-PD1 engagement, as well as increased distance between tumor cells and CD8+PD-1+ T cells. There was a significant decrease in HLA class I expression in the immunotherapy cohort at the time of AR compared with the chemotherapy (P = .005) and the targeted therapy (P = .01) cohorts. CONCLUSION: These findings highlight the genomic and immunophenotypic heterogeneity of ICI resistance in NSCLC, which will need to be considered when developing novel therapeutic strategies aimed at overcoming resistance.

Comprehensive genomic characterization of HER2-low and HER2-0 breast cancer.

The molecular underpinnings of HER2-low and HER2-0 (IHC 0) breast tumors remain poorly defined. Using genomic findings from 1039 patients with HER2-negative metastatic breast cancer undergoing next-generation sequencing from 7/2013-12/2020, we compare results between HER2-low (n = 487, 47%) and HER2-0 tumors (n = 552, 53%). A significantly higher number of ERBB2 alleles (median copy count: 2.05) are observed among HER2-low tumors compared to HER2-0 (median copy count: 1.79; P = 2.36e-6), with HER2-0 tumors harboring a higher rate of ERBB2 hemideletions (31.1% vs. 14.5%). No other genomic alteration reaches significance after accounting for multiple hypothesis testing, and no significant differences in tumor mutational burden are observed between HER2-low and HER2-0 tumors (median: 7.26 mutations/megabase vs. 7.60 mutations/megabase, p = 0.24). Here, we show that the genomic landscape of HER2-low and HER2-0 tumors does not differ significantly, apart from a higher ERBB2 copy count among HER2-low tumors, and a higher rate of ERBB2 hemideletions in HER2-0 tumors.

XRN1 deletion induces PKR-dependent cell lethality in interferon-activated cancer cells.

Emerging data suggest that induction of viral mimicry responses through activation of double-stranded RNA (dsRNA) sensors in cancer cells is a promising therapeutic strategy. One approach to induce viral mimicry is to target molecular regulators of dsRNA sensing pathways. Here, we show that the exoribonuclease XRN1 is a negative regulator of the dsRNA sensor protein kinase R (PKR) in cancer cells with high interferon-stimulated gene (ISG) expression. XRN1 deletion causes PKR activation and consequent cancer cell lethality. Disruption of interferon signaling with the JAK1/2 inhibitor ruxolitinib can decrease cellular PKR levels and rescue sensitivity to XRN1 deletion. Conversely, interferon-ß stimulation can increase PKR levels and induce sensitivity to XRN1 inactivation. Lastly, XRN1 deletion causes accumulation of endogenous complementary sense/anti-sense RNAs, which may represent candidate PKR ligands. Our data demonstrate how XRN1 regulates PKR and nominate XRN1 as a potential therapeutic target in cancer cells with an activated interferon cell state.

Loss of heterozygosity does not occur in BRCA1/2 mutant pediatric solid and central nervous system tumors.

Utilization of tumor-only sequencing has expanded in pediatric cancer patients, which can lead to identification of pathogenic variants in genes that may be germline and/or have uncertain relevance to the tumor in question, such as the homologous recombination (HR) pathway genes BRCA1/2. We identified patients with pathogenic BRCA1/2 mutations from somatic tumor sequencing, and performed additional germline sequencing to assess for the presence of loss of heterozygosity (LOH). Of seven patients identified, four (57.1%) mutations were found in the germline and none had associated LOH. Our data suggest that BRCA1/2 mutations identified in this context are likely incidental findings.

Oncogenic Drivers and Therapeutic Vulnerabilities in KRAS Wild-Type Pancreatic Cancer.

PURPOSE: Approximately 8% to 10% of pancreatic ductal adenocarcinomas (PDAC) do not harbor mutations in KRAS. Understanding the unique molecular and clinical features of this subset of pancreatic cancer is important to guide patient stratification for clinical trials of molecularly targeted agents. EXPERIMENTAL DESIGN: We analyzed a single-institution cohort of 795 exocrine pancreatic cancer cases (including 785 PDAC cases) with a targeted multigene sequencing panel and identified 73 patients (9.2%) with KRAS wild-type (WT) pancreatic cancer. RESULTS: Overall, 43.8% (32/73) of KRAS WT cases had evidence of an alternative driver of the MAPK pathway, including BRAF mutations and in-frame deletions and receptor tyrosine kinase fusions. Conversely, 56.2% of cases did not harbor a clear MAPK driver alteration, but 29.3% of these MAPK-negative KRAS WT cases (12/41) demonstrated activating alterations in other oncogenic drivers, such as GNAS, MYC, PIK3CA, and CTNNB1. We demonstrate potent efficacy of pan-RAF and MEK inhibition in patient-derived organoid models carrying BRAF in-frame deletions. Moreover, we demonstrate durable clinical benefit of targeted therapy in a patient harboring a KRAS WT tumor with a ROS1 fusion. Clinically, patients with KRAS WT tumors were significantly younger in age of onset (median age: 62.6 vs. 65.7 years; P = 0.037). SMAD4 mutations were associated with a particularly poor prognosis in KRAS WT cases. CONCLUSIONS: This study defines the genomic underpinnings of KRAS WT pancreatic cancer and highlights potential therapeutic avenues for future investigation in molecularly directed clinical trials. See related commentary by Kato et al., p. 4527.

Cancer aneuploidies are shaped primarily by effects on tumour fitness.

Aneuploidies-whole-chromosome or whole-arm imbalances-are the most prevalent alteration in cancer genomes1,2. However, it is still debated whether their prevalence is due to selection or ease of generation as passenger events1,2. Here we developed a method, BISCUT, that identifies loci subject to fitness advantages or disadvantages by interrogating length distributions of telomere- or centromere-bounded copy-number events. These loci were significantly enriched for known cancer driver genes, including genes not detected through analysis of focal copy-number events, and were often lineage specific. BISCUT identified the helicase-encoding gene WRN as a haploinsufficient tumour-suppressor gene on chromosome 8p, which is supported by several lines of evidence. We also formally quantified the role of selection and mechanical biases in driving aneuploidy, finding that rates of arm-level copy-number alterations are most highly correlated with their effects on cellular fitness1,2. These results provide insight into the driving forces behind aneuploidy and its contribution to tumorigenesis.

Impact of Aneuploidy and Chromosome 9p Loss on Tumor Immune Microenvironment and Immune Checkpoint Inhibitor Efficacy in NSCLC.

INTRODUCTION: Although gene-level copy number alterations have been studied as a potential biomarker of immunotherapy efficacy in NSCLC, the impact of aneuploidy burden and chromosomal arm-level events on immune checkpoint inhibitor (ICI) efficacy in NSCLC is uncertain. METHODS: Patients who received programmed cell death protein 1 or programmed death-ligand 1 (PD-L1) inhibitor at two academic centers were included. Across all 22 chromosomes analyzed, an arm was considered altered if at least 70% of its territory was either gained or deleted. Among nonsquamous NSCLCs which underwent targeted next-generation sequencing, we retrospectively quantified aneuploidy using the adjusted fraction of chromosomal arm alterations (FAA), defined as the number of altered chromosome arms divided by the number of chromosome arms assessed, adjusted for tumor purity. RESULTS: Among 2293 nonsquamous NSCLCs identified, the median FAA increased with more advanced cancer stage and decreased with higher PD-L1 tumor proportion score (TPS) levels (median FAA in TPS < 1%: 0.09, TPS 1%-49%: 0.08, TPS ≥ 50%: 0.05, p < 0.0001). There was a very weak correlation between FAA and tumor mutational burden when taken as continuous variables (R: 0.07, p = 0.0005). A total of 765 advanced nonsquamous NSCLCs with available FAA values were treated with ICIs. With decreasing FAA tertiles, there was a progressive improvement in objective response rate (ORR 15.1% in upper tertile versus 23.2% in middle tertile versus 28.4% in lowest tertile, p = 0.001), median progression-free survival (mPFS 2.5 versus 3.3 versus 4.1 mo, p < 0.0001), and median overall survival (mOS 12.5 versus 13.9 versus 16.4 mo, p = 0.006), respectively. In the arm-level enrichment analysis, chromosome 9p loss (OR = 0.22, Q = 0.0002) and chromosome 1q gain (OR = 0.43, Q = 0.002) were significantly enriched in ICI nonresponders after false discovery rate adjustment. Compared with NSCLCs without chromosome 9p loss (n = 452), those with 9p loss (n = 154) had a lower ORR (28.1% versus 7.8%, p < 0.0001), a shorter mPFS (4.1 versus 2.3 mo, p < 0.0001), and a shorter mOS (18.0 versus 9.6 mo, p < 0.0001) to immunotherapy. In addition, among NSCLCs with high PD-L1 expression (TPS ≥ 50%), chromosome 9p loss was associated with lower ORR (43% versus 6%, p < 0.0001), shorter mPFS (6.4 versus 2.6 mo, p = 0.0006), and shorter mOS (30.2 versus 14.3 mo, p = 0.0008) to immunotherapy compared with NSCLCs without 9p loss. In multivariable analysis, adjusting for key variables including FAA, chromosome 9p loss, but not 1q gain, retained a significant impact on ORR (hazard ratio [HR] = 0.25, p < 0.001), mPFS (HR = 1.49, p = 0.001), and mOS (HR = 1.47, p = 0.003). Multiplexed immunofluorescence and computational deconvolution of RNA sequencing data revealed that tumors with either high FAA levels or chromosome 9p loss had significantly fewer tumor-associated cytotoxic immune cells. CONCLUSIONS: Nonsquamous NSCLCs with high aneuploidy and chromosome 9p loss have a distinct tumor immune microenvironment and less favorable outcomes to ICIs.

Clinicopathologic, Genomic, and Immunophenotypic Landscape of ATM Mutations in Non-Small Cell Lung Cancer.

PURPOSE: ATM is the most commonly mutated DNA damage and repair gene in non-small cell lung cancer (NSCLC); however, limited characterization has been pursued. EXPERIMENTAL DESIGN: Clinicopathologic, genomic, and treatment data were collected for 5,172 patients with NSCLC tumors which underwent genomic profiling. ATM IHC was performed on 182 NSCLCs with ATM mutations. Multiplexed immunofluorescence was performed on a subset of 535 samples to examine tumor-infiltrating immune cell subsets. RESULTS: A total of 562 deleterious ATM mutations were identified in 9.7% of NSCLC samples. ATM-mutant (ATMMUT) NSCLC was significantly associated with female sex (P = 0.02), ever smoking status (P < 0.001), non-squamous histology (P = 0.004), and higher tumor mutational burden (DFCI, P < 0.0001; MSK, P < 0.0001) compared with ATM-wild-type (ATMWT) cases. Among 3,687 NSCLCs with comprehensive genomic profiling, co-occurring KRAS, STK11, and ARID2 oncogenic mutations were significantly enriched among ATMMUT NSCLCs (Q < 0.05), while TP53 and EGFR mutations were enriched in ATMWT NSCLCs. Among 182 ATMMUT samples with ATM IHC, tumors with nonsense, insertions/deletions, or splice site mutations were significantly more likely to display ATM loss by IHC (71.4% vs. 28.6%; P < 0.0001) compared with tumors with only predicted pathogenic missense mutations. Clinical outcomes to PD-(L)1 monotherapy (N = 1,522) and chemo-immunotherapy (N = 951) were similar between ATMMUT and ATMWT NSCLCs. Patients with concurrent ATM/TP53 mutations had significantly improved response rate and progression-free survival with PD-(L)1 monotherapy. CONCLUSIONS: Deleterious ATM mutations defined a subset of NSCLC with unique clinicopathologic, genomic, and immunophenotypic features. Our data may serve as resource to guide interpretation of specific ATM mutations in NSCLC.

Clinical Targeted Next-Generation Panel Sequencing Reveals MYC Amplification Is a Poor Prognostic Factor in Osteosarcoma.

PURPOSE: Osteosarcoma risk stratification, on the basis of the presence of metastatic disease at diagnosis and histologic response to chemotherapy, has remained unchanged for four decades, does not include genomic features, and has not facilitated treatment advances. We report on the genomic features of advanced osteosarcoma and provide evidence that genomic alterations can be used for risk stratification. MATERIALS AND METHODS: In a primary analytic patient cohort, 113 tumor and 69 normal samples from 92 patients with high-grade osteosarcoma were sequenced with OncoPanel, a targeted next-generation sequencing assay. In this primary cohort, we assessed the genomic landscape of advanced disease and evaluated the correlation between recurrent genomic events and outcome. We assessed whether prognostic associations identified in the primary cohort were maintained in a validation cohort of 86 patients with localized osteosarcoma tested with MSK-IMPACT. RESULTS: In the primary cohort, 3-year overall survival (OS) was 65%. Metastatic disease, present in 33% of patients at diagnosis, was associated with poor OS (P = .04). The most frequently altered genes in the primary cohort were TP53, RB1, MYC, CCNE1, CCND3, CDKN2A/B, and ATRX. Mutational signature 3 was present in 28% of samples. MYC amplification was associated with a worse 3-year OS in both the primary cohort (P = .015) and the validation cohort (P = .012). CONCLUSION: The most frequently occurring genomic events in advanced osteosarcoma were similar to those described in prior reports. MYC amplification, detected with clinical targeted next-generation sequencing panel tests, is associated with poorer outcomes in two independent cohorts.

Velcrin-induced selective cleavage of tRNALeu(TAA) by SLFN12 causes cancer cell death.

Velcrin compounds kill cancer cells expressing high levels of phosphodiesterase 3A (PDE3A) and Schlafen family member 12 (SLFN12) by inducing complex formation between these two proteins, but the mechanism of cancer cell killing by the PDE3A-SLFN12 complex is not fully understood. Here, we report that the physiological substrate of SLFN12 RNase is tRNALeu(TAA). SLFN12 selectively digests tRNALeu(TAA), and velcrin treatment promotes the cleavage of tRNALeu(TAA) by inducing PDE3A-SLFN12 complex formation in vitro. We found that distinct sequences in the variable loop and acceptor stem of tRNALeu(TAA) are required for substrate digestion. Velcrin treatment of sensitive cells results in downregulation of tRNALeu(TAA), ribosome pausing at Leu-TTA codons and global inhibition of protein synthesis. Velcrin-induced cleavage of tRNALeu(TAA) by SLFN12 and the concomitant global inhibition of protein synthesis thus define a new mechanism of apoptosis initiation.

Rare FGFR Oncogenic Alterations in Sequenced Pediatric Solid and Brain Tumors Suggest FGFR Is a Relevant Molecular Target in Childhood Cancer.

PURPOSE: Multiple FGFR inhibitors are currently in clinical trials enrolling adults with different solid tumors, while very few enroll pediatric patients. We determined the types and frequency of FGFR alterations (FGFR1-4) in pediatric cancers to inform future clinical trial design. METHODS: Tumors with FGFR alterations were identified from two large cohorts of pediatric solid tumors subjected to targeted DNA sequencing: The Dana-Farber/Boston Children's Profile Study (n = 888) and the multi-institution GAIN/iCAT2 (Genomic Assessment Improves Novel Therapy) Study (n = 571). Data from the combined patient population of 1,395 cases (64 patients were enrolled in both studies) were reviewed and cases in which an FGFR alteration was identified by OncoPanel sequencing were further assessed. RESULTS: We identified 41 patients with tumors harboring an oncogenic FGFR alteration. Median age at diagnosis was 8 years (range, 6 months-26 years). Diagnoses included 11 rhabdomyosarcomas, nine low-grade gliomas, and 17 other tumor types. Alterations included gain-of-function sequence variants (n = 19), amplifications (n = 10), oncogenic fusions (FGFR3::TACC3 [n = 3], FGFR1::TACC1 [n = 1], FGFR1::EBF2 [n = 1], FGFR1::CLIP2 [n = 1], and FGFR2::CTNNA3 [n = 1]), pathogenic-leaning variants of uncertain significance (n = 4), and amplification in combination with a pathogenic-leaning variant of uncertain significance (n = 1). Two novel FGFR1 fusions in two different patients were identified in this cohort, one of whom showed a response to an FGFR inhibitor. CONCLUSION: In summary, activating FGFR alterations were found in approximately 3% (41/1,395) of pediatric solid tumors, identifying a population of children with cancer who may be eligible and good candidates for trials evaluating FGFR-targeted therapy. Importantly, the genomic and clinical data from this study can help inform drug development in accordance with the Research to Accelerate Cures and Equity for Children Act.

Discovery and Structure-Based Design of Potent Covalent PPARγ Inverse-Agonists BAY-4931 and BAY-0069.

The ligand-activated nuclear receptor peroxisome-proliferator-activated receptor-γ (PPARG or PPARγ) represents a potential target for a new generation of cancer therapeutics, especially in muscle-invasive luminal bladder cancer where PPARγ is a critical lineage driver. Here we disclose the discovery of a series of chloro-nitro-arene covalent inverse-agonists of PPARγ that exploit a benzoxazole core to improve interactions with corepressors NCOR1 and NCOR2. In vitro treatment of sensitive cell lines with these compounds results in the robust regulation of PPARγ target genes and antiproliferative effects. Despite their imperfect physicochemical properties, the compounds showed modest pharmacodynamic target regulation in vivo. Improvements to the in vitro potency and efficacy of BAY-4931 and BAY-0069 compared to those of previously described PPARγ inverse-agonists show that these compounds are novel tools for probing the in vitro biology of PPARγ inverse-agonism.

Tangent normalization for somatic copy-number inference in cancer genome analysis.

MOTIVATION: Somatic copy-number alterations (SCNAs) play an important role in cancer development. Systematic noise in sequencing and array data present a significant challenge to the inference of SCNAs for cancer genome analyses. As part of The Cancer Genome Atlas, the Broad Institute Genome Characterization Center developed the Tangent normalization method to generate copy-number profiles using data from single-nucleotide polymorphism (SNP) arrays and whole-exome sequencing (WES) technologies for over 10 000 pairs of tumors and matched normal samples. Here, we describe the Tangent method, which uses a unique linear combination of normal samples as a reference for each tumor sample, to subtract systematic errors that vary across samples. We also describe a modification of Tangent, called Pseudo-Tangent, which enables denoising through comparisons between tumor profiles when few normal samples are available. RESULTS: Tangent normalization substantially increases signal-to-noise ratios (SNRs) compared to conventional normalization methods in both SNP array and WES analyses. Tangent and Pseudo-Tangent normalizations improve the SNR by reducing noise with minimal effect on signal and exceed the contribution of other steps in the analysis such as choice of segmentation algorithm. Tangent and Pseudo-Tangent are broadly applicable and enable more accurate inference of SCNAs from DNA sequencing and array data. AVAILABILITY AND IMPLEMENTATION: Tangent is available at https://github.com/broadinstitute/tangent and as a Docker image (https://hub.docker.com/r/broadinstitute/tangent). Tangent is also the normalization method for the copy-number pipeline in Genome Analysis Toolkit 4 (GATK4). SUPPLEMENTARY INFORMATION: Supplementary data are available at Bioinformatics online.