DrugMap: A quantitative pan-cancer analysis of cysteine ligandability

Cysteine-focused chemical proteomic platforms have accelerated the clinical development of covalent inhibitors for a wide range of targets in cancer. However, how different oncogenic contexts influence cysteine targeting remains unknown. To address this question, we have developed "DrugMap," an atlas of cysteine ligandability compiled across 416 cancer cell lines. We unexpectedly find that cysteine ligandability varies across cancer cell lines, and we attribute this to differences in cellular redox states, protein conformational changes, and genetic mutations. Leveraging these findings, we identify actionable cysteines in NF-κB1 and SOX10 and develop corresponding covalent ligands that block the activity of these transcription factors. We demonstrate that the NF-κB1 probe blocks DNA binding, whereas the SOX10 ligand increases SOX10-SOX10 interactions and disrupts melanoma transcriptional signaling. Our findings reveal heterogeneity in cysteine ligandability across cancers, pinpoint cell-intrinsic features driving cysteine targeting, and illustrate the use of covalent probes to disrupt oncogenic transcription-factor activity.

Patterns in the tapestry of chromatin-bound RB

The retinoblastoma protein (RB)-mediated regulation of E2F is a component of a highly conserved cell cycle machine. However, RB's tumor suppressor activity, like RB's requirement in animal development, is tissue-specific, context-specific, and sometimes appears uncoupled from cell proliferation. Detailed new information about RB's genomic distribution provides a new perspective on the complexity of RB function, suggesting that some of its functional specificity results from context-specific RB association with chromatin. Here we summarize recent evidence showing that RB targets different types of chromatin regulatory elements at different cell cycle stages. RB controls traditional RB/E2F targets prior to S-phase, but, when cells proliferate, RB redistributes to cell type-specific chromatin loci. We discuss the broad implications of the new data for RB research.

Distinguishing preferences of human APOBEC3A and APOBEC3B for cytosines in hairpin loops, and reflection of these preferences in APOBEC-signature cancer genome mutations

The APOBEC3 enzymes convert cytosines in single-stranded DNA to uracils to protect against viruses and retrotransposons but can contribute to mutations that diversify tumors. To understand the mechanism of mutagenesis, we map the uracils resulting from expression of APOBEC3B or its catalytic carboxy-terminal domain (CTD) in Escherichia coli. Like APOBEC3A, the uracilomes of A3B and A3B-CTD show a preference to deaminate cytosines near transcription start sites and the lagging-strand replication templates and in hairpin loops. Both biochemical activities of the enzymes and genomic uracil distribution show that A3A prefers 3 nt loops the best, while A3B prefers 4 nt loops. Reanalysis of hairpin loop mutations in human tumors finds intrinsic characteristics of both the enzymes, with a much stronger contribution from A3A. We apply Hairpin Signatures 1 and 2, which define A3A and A3B preferences respectively and are orthogonal to published methods, to evaluate their contribution to human tumor mutations.

Mesoscale DNA features impact APOBEC3A and APOBEC3B deaminase activity and shape tumor mutational landscapes

Antiviral DNA cytosine deaminases APOBEC3A and APOBEC3B are major sources of mutations in cancer by catalyzing cytosine-to-uracil deamination. APOBEC3A preferentially targets single-stranded DNAs, with a noted affinity for DNA regions that adopt stem-loop secondary structures. However, the detailed substrate preferences of APOBEC3A and APOBEC3B have not been fully established, and the specific influence of the DNA sequence on APOBEC3A and APOBEC3B deaminase activity remains to be investigated. Here, we find that APOBEC3B also selectively targets DNA stem-loop structures, and they are distinct from those subjected to deamination by APOBEC3A. We develop Oligo-seq, an in vitro sequencing-based method to identify specific sequence contexts promoting APOBEC3A and APOBEC3B activity. Through this approach, we demonstrate that APOBEC3A and APOBEC3B deaminase activity is strongly regulated by specific sequences surrounding the targeted cytosine. Moreover, we identify the structural features of APOBEC3B and APOBEC3A responsible for their substrate preferences. Importantly, we determine that APOBEC3B-induced mutations in hairpin-forming sequences within tumor genomes differ from the DNA stem-loop sequences mutated by APOBEC3A. Together, our study provides evidence that APOBEC3A and APOBEC3B can generate distinct mutation landscapes in cancer genomes, driven by their unique substrate selectivity.

Tumor cell-based liquid biopsy using high-throughput microfluidic enrichment of entire leukapheresis product

Circulating Tumor Cells (CTCs), interrogated by sampling blood from patients with cancer, contain multiple analytes, including intact RNA, high molecular weight DNA, proteins, and metabolic markers. However, the clinical utility of tumor cell-based liquid biopsy has been limited since CTCs are very rare, and current technologies cannot process the blood volumes required to isolate a sufficient number of tumor cells for in-depth assays. We previously described a high-throughput microfluidic prototype utilizing high-flow channels and amplification of cell sorting forces through magnetic lenses. Here, we apply this technology to analyze patient-derived leukapheresis products, interrogating a mean blood volume of 5.83 liters from patients with metastatic cancer, with a median of 2,799 CTCs purified per patient. Isolation of many CTCs from individual patients enables characterization of their morphological and molecular heterogeneity, including cell and nuclear size and RNA expression. It also allows robust detection of gene copy number variation, a definitive cancer marker with potential diagnostic applications. High-volume microfluidic enrichment of CTCs constitutes a new dimension in liquid biopsies.

Acquired Cross-resistance in Small Cell Lung Cancer due to Extrachromosomal DNA Amplification of MYC paralogs

Small cell lung cancer (SCLC) presents as a highly chemosensitive malignancy but acquires cross-resistance after relapse. This transformation is nearly inevitable in patients but has been difficult to capture in laboratory models. Here, we present a pre-clinical system that recapitulates acquired cross-resistance, developed from 51 patient-derived xenograft (PDX) models. Each model was tested in vivo against three clinical regimens: cisplatin plus etoposide, olaparib plus temozolomide, and topotecan. These drug-response profiles captured hallmark clinical features of SCLC, such as the emergence of treatment-refractory disease after early relapse. For one patient, serial PDX models revealed that cross-resistance was acquired through MYC amplification on extrachromosomal DNA (ecDNA). Genomic and transcriptional profiles of the full PDX panel revealed that MYC paralog amplifications on ecDNAs were recurrent in relapsed cross-resistant SCLC, and this was corroborated in tumor biopsies from relapsed patients. We conclude that ecDNAs with MYC paralogs are recurrent drivers of cross-resistance in SCLC.

Blood-based screening for HPV-associated cancers

Background

HPV-associated oropharyngeal cancer (HPV+OPSCC) is the most common HPV-associated cancer in the United States yet unlike cervical cancer lacks a screening test. HPV+OPSCCs are presumed to start developing 10-15 years prior to clinical diagnosis. Circulating tumor HPV DNA (ctHPVDNA) is a sensitive and specific biomarker for HPV+OPSCC. Taken together, blood-based screening for HPV+OPSCC may be feasible years prior to diagnosis.

Methods

We developed an HPV whole genome sequencing assay, HPV-DeepSeek, with 99% sensitivity and specificity at clinical diagnosis. 28 plasma samples from HPV+OPSCC patients collected 1.3-10.8 years prior to diagnosis along with 1:1 age and gender-matched controls were run on HPV-DeepSeek and an HPV serology assay.

Results

22/28 (79%) of cases and 0/28 controls screened positive for HPV+OPSCC with 100% detection within four years of diagnosis and a maximum lead time of 7.8 years. We next applied a machine learning model classifying 27/28 cases (96%) with 100% detection within 10 years. Plasma-based PIK3CA gene mutations, viral genome integration events and HPV serology were used to orthogonally validate cancer detection with 68% (19/28) of the cohort having multiple cancer signals detected. Molecular fingerprinting of HPV genomes was performed across patients demonstrating that each viral genome was unique, ruling out contamination. In patients with tumor blocks from diagnosis (15/28), molecular fingerprinting was performed within patients confirming the same viral genome across time.

Conclusions

We demonstrate accurate blood-based detection of HPV-associated cancers with lead times up to 10 years before clinical cancer diagnosis and in doing so, highlight the enormous potential of ctDNA-based cancer screening.

APOBEC3A induces DNA gaps through PRIMPOL and confers gap-associated therapeutic vulnerability

Mutation signatures associated with apolipoprotein B mRNA editing catalytic polypeptide-like 3A/B (APOBEC3A/B) cytidine deaminases are prevalent across cancers, implying their roles as mutagenic drivers during tumorigenesis and tumor evolution. APOBEC3A (A3A) expression induces DNA replication stress and increases the cellular dependency on the ataxia telangiectasia and Rad3-related (ATR) kinase for survival. Nonetheless, how A3A induces DNA replication stress remains unclear. We show that A3A induces replication stress without slowing replication forks. We find that A3A induces single-stranded DNA (ssDNA) gaps through PrimPol-mediated repriming. A3A-induced ssDNA gaps are repaired by multiple pathways involving ATR, RAD51, and translesion synthesis. Both ATR inhibition and trapping of poly(ADP-ribose) polymerase (PARP) on DNA by PARP inhibitor impair the repair of A3A-induced gaps, preferentially killing A3A-expressing cells. When used in combination, PARP and ATR inhibitors selectively kill A3A-expressing cells synergistically in a manner dependent on PrimPol-generated gaps. Thus, A3A-induced replication stress arises from PrimPol-generated ssDNA gaps, which confer a therapeutic vulnerability to gap-targeted DNA repair inhibitors.

DrugMap: A quantitative pan-cancer analysis of cysteine ligandability

Cysteine-focused chemical proteomic platforms have accelerated the clinical development of covalent inhibitors of a wide-range of targets in cancer. However, how different oncogenic contexts influence cysteine targeting remains unknown. To address this question, we have developed DrugMap , an atlas of cysteine ligandability compiled across 416 cancer cell lines. We unexpectedly find that cysteine ligandability varies across cancer cell lines, and we attribute this to differences in cellular redox states, protein conformational changes, and genetic mutations. Leveraging these findings, we identify actionable cysteines in NFκB1 and SOX10 and develop corresponding covalent ligands that block the activity of these transcription factors. We demonstrate that the NFκB1 probe blocks DNA binding, whereas the SOX10 ligand increases SOX10-SOX10 interactions and disrupts melanoma transcriptional signaling. Our findings reveal heterogeneity in cysteine ligandability across cancers, pinpoint cell-intrinsic features driving cysteine targeting, and illustrate the use of covalent probes to disrupt oncogenic transcription factor activity.

ATR inhibition induces synthetic lethality in mismatch repair-deficient cells and augments immunotherapy

The mismatch repair (MMR) deficiency of cancer cells drives mutagenesis and offers a useful biomarker for immunotherapy. However, many MMR-deficient (MMR-d) tumors do not respond to immunotherapy, highlighting the need for alternative approaches to target MMR-d cancer cells. Here, we show that inhibition of the ATR kinase preferentially kills MMR-d cancer cells. Mechanistically, ATR inhibitor (ATRi) imposes synthetic lethality on MMR-d cells by inducing DNA damage in a replication- and MUS81 nuclease-dependent manner. The DNA damage induced by ATRi is colocalized with both MSH2 and PCNA, suggesting that it arises from DNA structures recognized by MMR proteins during replication. In syngeneic mouse models, ATRi effectively reduces the growth of MMR-d tumors. Interestingly, the antitumor effects of ATRi are partially due to CD8+ T cells. In MMR-d cells, ATRi stimulates the accumulation of nascent DNA fragments in the cytoplasm, activating the cGAS-mediated interferon response. The combination of ATRi and anti-PD-1 antibody reduces the growth of MMR-d tumors more efficiently than ATRi or anti-PD-1 alone, showing the ability of ATRi to augment the immunotherapy of MMR-d tumors. Thus, ATRi selectively targets MMR-d tumor cells by inducing synthetic lethality and enhancing antitumor immunity, providing a promising strategy to complement and augment MMR deficiency-guided immunotherapy.

The RNA m5C modification in R-loops as an off switch of Alt-NHEJ

The roles of R-loops and RNA modifications in homologous recombination (HR) and other DNA double-stranded break (DSB) repair pathways remain poorly understood. Here, we find that DNA damage-induced RNA methyl-5-cytosine (m5C) modification in R-loops plays a crucial role to regulate PARP1-mediated poly ADP-ribosylation (PARylation) and the choice of DSB repair pathways at sites of R-loops. Through bisulfite sequencing, we discover that the methyltransferase TRDMT1 preferentially generates m5C after DNA damage in R-loops across the genome. In the absence of m5C, R-loops activate PARP1-mediated PARylation both in vitro and in cells. Concurrently, m5C promotes transcription-coupled HR (TC-HR) while suppressing PARP1-dependent alternative non-homologous end joining (Alt-NHEJ), favoring TC-HR over Alt-NHEJ in transcribed regions as the preferred repair pathway. Importantly, simultaneous disruption of both TC-HR and Alt-NHEJ with TRDMT1 and PARP or Polymerase θ inhibitors prevents alternative DSB repair and exhibits synergistic cytotoxic effects on cancer cells, suggesting an effective strategy to exploit genomic instability in cancer therapy.

Mesoscale DNA Features Impact APOBEC3A and APOBEC3B Deaminase Activity and Shape Tumor Mutational Landscapes

Antiviral DNA cytosine deaminases APOBEC3A and APOBEC3B are major sources of mutations in cancer by catalyzing cytosine-to-uracil deamination. APOBEC3A preferentially targets singlestranded DNAs, with a noted affinity for DNA regions that adopt stem-loop secondary structures. However, the detailed substrate preferences of APOBEC3A and APOBEC3B have been fully established, and the specific influence of the DNA sequence on APOBEC3A APOBEC3B deaminase activity remains to be investigated. Here, we find that APOBEC3B selectively targets DNA stem-loop structures, and they are distinct from those subjected deamination by APOBEC3A. We develop Oligo-seq, a novel in vitro sequencing-based to identify specific sequence contexts promoting APOBEC3A and APOBEC3B activity. Through this approach, we demonstrate that APOBEC3A an APOBEC3B deaminase activity is strongly regulated by specific sequences surrounding the targeted cytosine. Moreover, we identify structural features of APOBEC3B and APOBEC3A responsible for their substrate preferences. Importantly, we determine that APOBEC3B-induced mutations in hairpin-forming sequences within tumor genomes differ from the DNA stem-loop sequences mutated by APOBEC3A. Together, our study provides evidence that APOBEC3A and APOBEC3B can generate mutation landscapes in cancer genomes, driven by their unique substrate selectivity.

Distinguishing preferences of human APOBEC3A and APOBEC3B for cytosines in hairpin loops, and reflection of these preferences in APOBEC-signature cancer genome mutations

The APOBEC3 family of enzymes convert cytosines in single-stranded DNA to uracils thereby causing mutations. These enzymes protect human cells against viruses and retrotransposons, but in many cancers they contribute to mutations that diversify the tumors and help them escape anticancer drug treatments. To understand the mechanism of mutagenesis by APOBEC3B, we expressed the complete enzyme or its catalytic carboxy-terminal domain (CTD) in repair-deficient Eschericia coli and mapped the resulting uracils using uracil pull-down and sequencing technology. The uracilomes of A3B-full and A3B-CTD showed peaks in many of the same regions where APOBEC3A also created uracilation peaks. Like A3A, the two A3B enzymes also preferred to deaminate cytosines near transcription start sites and in the lagging-strand template at replication forks. In contrast to an earlier report that A3B does not favor hairpin loops over linear DNA, we found that both A3B-full and A3B-CTD showed a strong preference for cytosines in hairpin loops. The major difference between A3A and A3B was that while the former enzyme prefers 3 nt loops the best, A3B prefers loops of 4 nt over those of other lengths. Furthermore, within 5 nt loops, A3A prefers cytosine to be in the penultimate position, while A3B prefers it to be at the 3' end of the loop. We confirmed these loop size and sequence preferences experimentally using purified A3A and A3B-CTD proteins. Reanalysis of hairpin loop mutations in human tumors using the size, sequence and position preferences of the two enzymes found that the tumors displayed mutations with intrinsic characteristics of both the enzymes with a stronger contribution from A3A.

Therapy-induced APOBEC3A drives evolution of persistent cancer cells

Acquired drug resistance to anticancer targeted therapies remains an unsolved clinical problem. Although many drivers of acquired drug resistance have been identified1-4, the underlying molecular mechanisms shaping tumour evolution during treatment are incompletely understood. Genomic profiling of patient tumours has implicated apolipoprotein B messenger RNA editing catalytic polypeptide-like (APOBEC) cytidine deaminases in tumour evolution; however, their role during therapy and the development of acquired drug resistance is undefined. Here we report that lung cancer targeted therapies commonly used in the clinic can induce cytidine deaminase APOBEC3A (A3A), leading to sustained mutagenesis in drug-tolerant cancer cells persisting during therapy. Therapy-induced A3A promotes the formation of double-strand DNA breaks, increasing genomic instability in drug-tolerant persisters. Deletion of A3A reduces APOBEC mutations and structural variations in persister cells and delays the development of drug resistance. APOBEC mutational signatures are enriched in tumours from patients with lung cancer who progressed after extended responses to targeted therapies. This study shows that induction of A3A in response to targeted therapies drives evolution of drug-tolerant persister cells, suggesting that suppression of A3A expression or activity may represent a potential therapeutic strategy in the prevention or delay of acquired resistance to lung cancer targeted therapy.

Acquired Cross-resistance in Small Cell Lung Cancer due to Extrachromosomal DNA Amplification of MYC paralogs

Small cell lung cancer (SCLC) presents as a highly chemosensitive malignancy but acquires cross-resistance after relapse. This transformation is nearly inevitable in patients but has been difficult to capture in laboratory models. Here we present a pre-clinical system that recapitulates acquired cross-resistance in SCLC, developed from 51 patient-derived xenografts (PDXs). Each model was tested for in vivo sensitivity to three clinical regimens: cisplatin plus etoposide, olaparib plus temozolomide, and topotecan. These functional profiles captured hallmark clinical features, such as the emergence of treatment-refractory disease after early relapse. Serially derived PDX models from the same patient revealed that cross-resistance was acquired through a MYC amplification on extrachromosomal DNA (ecDNA). Genomic and transcriptional profiles of the full PDX panel revealed that this was not unique to one patient, as MYC paralog amplifications on ecDNAs were recurrent among cross-resistant models derived from patients after relapse. We conclude that ecDNAs with MYC paralogs are recurrent drivers of cross-resistance in SCLC.

SIGNIFICANCE

SCLC is initially chemosensitive, but acquired cross-resistance renders this disease refractory to further treatment and ultimately fatal. The genomic drivers of this transformation are unknown. We use a population of PDX models to discover that amplifications of MYC paralogs on ecDNA are recurrent drivers of acquired cross-resistance in SCLC.

DNA hypomethylation silences anti-tumor immune genes in early prostate cancer and CTCs

Cancer is characterized by hypomethylation-associated silencing of large chromatin domains, whose contribution to tumorigenesis is uncertain. Through high-resolution genome-wide single-cell DNA methylation sequencing, we identify 40 core domains that are uniformly hypomethylated from the earliest detectable stages of prostate malignancy through metastatic circulating tumor cells (CTCs). Nested among these repressive domains are smaller loci with preserved methylation that escape silencing and are enriched for cell proliferation genes. Transcriptionally silenced genes within the core hypomethylated domains are enriched for immune-related genes; prominent among these is a single gene cluster harboring all five CD1 genes that present lipid antigens to NKT cells and four IFI16-related interferon-inducible genes implicated in innate immunity. The re-expression of CD1 or IFI16 murine orthologs in immuno-competent mice abrogates tumorigenesis, accompanied by the activation of anti-tumor immunity. Thus, early epigenetic changes may shape tumorigenesis, targeting co-located genes within defined chromosomal loci. Hypomethylation domains are detectable in blood specimens enriched for CTCs.

Chromatin-bound protein colocalization analysis using bedGraph2Cluster and PanChIP

Computational pipelines for chromatin immunoprecipitation sequencing analysis can neglect colocalization events that occur in a mere subset of the genome. Here, we detail a streamlined approach for assessing colocalization of chromatin-bound proteins using the bedGraph2Cluster and PanChIP algorithms. Using histone modifications as an example, bedGraph2Cluster performs clustering analysis on chromatin binding patterns of target proteins. PanChIP then compares these clusters with a reference library of chromatin binding patterns and measures the overlap in peaks, capturing the heterogeneity in chromatin binding and colocalization patterns. For complete details on the use and execution of this protocol, please refer to Sanidas et al. (2022).1.

TGF-β in the microenvironment induces a physiologically occurring immune-suppressive senescent state

TGF-β induces senescence in embryonic tissues. Whether TGF-β in the hypoxic tumor microenvironment (TME) induces senescence in cancer and how the ensuing senescence-associated secretory phenotype (SASP) remodels the cellular TME to influence immune checkpoint inhibitor (ICI) responses are unknown. We show that TGF-β induces a deeper senescent state under hypoxia than under normoxia; deep senescence correlates with the degree of E2F suppression and is marked by multinucleation, reduced reentry into proliferation, and a distinct 14-gene SASP. Suppressing TGF-β signaling in tumors in an immunocompetent mouse lung cancer model abrogates endogenous senescent cells and suppresses the 14-gene SASP and immune infiltration. Untreated human lung cancers with a high 14-gene SASP display immunosuppressive immune infiltration. In a lung cancer clinical trial of ICIs, elevated 14-gene SASP is associated with increased senescence, TGF-β and hypoxia signaling, and poor progression-free survival. Thus, TME-induced senescence may represent a naturally occurring state in cancer, contributing to an immune-suppressive phenotype associated with immune therapy resistance.

Abstract IA010: Hypomethylation of immune-related loci contribute to prostate tumorigenesis

Cancer-associated alterations in DNA methylation include hypermethylation of focal CpG islands within promoters, as well as hypomethylation of large stretches of gene-poor chromosomal domains, so called Partially Methylated Domains (PMDs). While CpG island hypermethylation is associated with silencing of tumor suppressor genes, the functional consequences of hypomethylation are uncertain. PMDs may be cause or consequence of cancer cell proliferation, and they have been associated with transcriptional de-repression of repetitive sequences, as well as histone modification-triggered silencing of rare protein encoding genes within these domains. In studying Circulating Tumor Cells (CTCs) isolated from blood specimens of patients with prostate cancer, we mapped PMDs at the single cell level, identifying a set of 40 “core PMDs” that are shared across cancer cells, within and across different individuals, as well as in prostate cancer cell lines. By analyzing single nuclei from micro-dissected early-stage prostate tumors, we find that hypomethylation of these core PMDs occurs early during prostate tumorigenesis. Core PMDs are highly enriched for immune-related genes, many of which are co-located within a single chromosomal locus targeted for hypomethylation-associated silencing. Nested between adjacent PMDs, we also identify domains that have retained DNA methylation (Preserved Methylation Islands; PMIs), which have sharply demarcated borders, spanning one or a small number of genes. In contrast to PMDs, PMIs are enriched for cell proliferation-associated genes. Together, these observations suggest that hypomethylation-associated silencing of immune-related genes, with sparing of adjacent proliferation-associated genes, is an early event in prostate cancer development. These large-scale chromosomal modifications are detectable within CTCs and may present diagnostic opportunities.

Citation Format: Daniel A. Haber, Hongshan Guo, Ben S. Wittner, Joanna A Vuille, Michael S. Lawrence, Chin-Lee Wu, Mehmet Toner, David T. Miyamoto, Shyamala Maheswaran. Hypomethylation of immune-related loci contribute to prostate tumorigenesis [abstract]. In: Proceedings of the AACR Special Conference: Cancer Metastasis; 2022 Nov 14-17; Portland, OR. Philadelphia (PA): AACR; Cancer Res 2022;83(2 Suppl_2):Abstract nr IA010.

Single-cell transcriptomic profiling for inferring tumor origin and mechanisms of therapeutic resistance

Head and Neck Squamous Cell Carcinoma (HNSCC) is an aggressive epithelial cancer with poor overall response rates to checkpoint inhibitor therapy (CPI) despite CPI being the recommended treatment for recurrent or metastatic HNSCC. Mechanisms of resistance to CPI in HNSCC are poorly understood. To identify drivers of response and resistance to CPI in a unique patient who was believed to have developed three separate HNSCCs, we performed single-cell RNA-seq (scRNA-seq) profiling of two responding lesions and one progressive lesion that developed during CPI. Our results not only suggest interferon-induced APOBEC3-mediated acquired resistance as a mechanism of CPI resistance in the progressing lesion but further, that the lesion in question was actually a metastasis as opposed to a new primary tumor, highlighting the immense power of scRNA-seq as a clinical tool for inferring tumor origin and mechanisms of therapeutic resistance.

Chromatin-bound RB targets promoters, enhancers, and CTCF-bound loci and is redistributed by cell-cycle progression

The interaction of RB with chromatin is key to understanding its molecular functions. Here, for first time, we identify the full spectrum of chromatin-bound RB. Rather than exclusively binding promoters, as is often described, RB targets three fundamentally different types of loci (promoters, enhancers, and insulators), which are largely distinguishable by the mutually exclusive presence of E2F1, c-Jun, and CTCF. While E2F/DP facilitates RB association with promoters, AP-1 recruits RB to enhancers. Although phosphorylation in CDK sites is often portrayed as releasing RB from chromatin, we show that the cell cycle redistributes RB so that it enriches at promoters in G1 and at non-promoter sites in cycling cells. RB-bound promoters include the classic E2F-targets and are similar between lineages, but RB-bound enhancers associate with different categories of genes and vary between cell types. Thus, RB has a well-preserved role controlling E2F in G1, and it targets cell-type-specific enhancers and CTCF sites when cells enter S-phase.

Abstract 657: Impact of therapy induced APOBEC3A mutagenesis on tumor evolution in non small cell lung cancer

Acquired drug resistance to even the most effective anti cancer targeted therapies remains an unsolved clinical problem. Although many drivers of acquired drug resistance have been identified, the underlying molecular mechanisms shaping tumor evolution during treatment are incompletely understood. We have seen that lung cancer targeted therapies commonly used in the clinic induce the expression of cytidine deaminase APOBEC3A (A3A), leading to sustained mutagenesis in drug tolerant cancer cells persisting during therapy. Preventing therapy induced A3A mutagenesis by gene deletion delayed the emergence of drug resistance. Here, we show that therapy induced A3A mutagenesis contributes to tumor evolution in NSCLC. Whole genome sequencing revealed that resistant clones that evolved from persistent drug tolerant cells (late clones) harbored more A3A mutations compared to the parental cell population than pre existing resistant clones (early clones). In a subset of NSCLC patients who received targeted therapies, we observed A3A mutations accompanied clonal evolution during treatment. Comparison of APOBEC mutation fractions in short vs long term responders suggests that short responders with acquired resistance mechanisms that evolved from pre existing resistant clones have less accumulation of APOBEC mutations. Collectively, these findings insist that an increase in mutagenic processes drives tumor evolution during targeted therapy treatment and leads to acquired resistance. Thus, suppressing expression or enzymatic activity of A3A may represent a potential therapeutic strategy to halt the evolution of resistant clones and prevent acquired resistance to lung cancer targeted therapy.

Citation Format: Hideko Isozaki, Ammal Abbasi, Naveed Nikpour, Marcello Stanzione, Ramin Sakhtemani, Susanna L. Monroe, Alice T. Shaw, Jessica J. Lin, Lecia V. Sequist, Zofia Piotrowska, Rémi Buisson, Michael S. Lawrence, Aaron N. Hata. Impact of therapy induced APOBEC3A mutagenesis on tumor evolution in non small cell lung cancer [abstract]. In: Proceedings of the American Association for Cancer Research Annual Meeting 2022; 2022 Apr 8-13. Philadelphia (PA): AACR; Cancer Res 2022;82(12_Suppl):Abstract nr 657.

Translesion DNA synthesis mediates acquired resistance to olaparib plus temozolomide in small cell lung cancer

In small cell lung cancer (SCLC), acquired resistance to DNA-damaging therapy is challenging to study because rebiopsy is rarely performed. We used patient-derived xenograft models, established before therapy and after progression, to dissect acquired resistance to olaparib plus temozolomide (OT), a promising experimental therapy for relapsed SCLC. These pairs of serial models reveal alterations in both cell cycle kinetics and DNA replication and demonstrate both inter- and intratumoral heterogeneity in mechanisms of resistance. In one model pair, up-regulation of translesion DNA synthesis (TLS) enabled tolerance of OT-induced damage during DNA replication. TLS inhibitors restored sensitivity to OT both in vitro and in vivo, and similar synergistic effects were seen in additional SCLC cell lines. This represents the first described mechanism of acquired resistance to DNA damage in a patient with SCLC and highlights the potential of the serial model approach to investigate and overcome resistance to therapy in SCLC.

Extracellular matrix proteins regulate NK cell function in peripheral tissues

Natural killer (NK) cells reject major histocompatibility complex class I (MHC-I)-deficient bone marrow through direct cytotoxicity but not solid organ transplants devoid of MHC-I. Here, we demonstrate an immediate switch in NK cell function upon exit from the circulation, characterized by a shift from direct cytotoxicity to chemokine/cytokine production. In the skin transplant paradigm, combining an NK cell-specific activating ligand, m157, with missing self MHC-I resulted in complete graft rejection, which was dependent on NK cells as potential helpers and T cells as effectors. Extracellular matrix proteins, collagen I, collagen III, and elastin, blocked NK cell cytotoxicity and promoted their chemokine/cytokine production. NK cell cytotoxicity against MHC-I-deficient melanoma in the skin was markedly increased by blocking tumor collagen deposition. MHC-I down-regulation occurred in solid human cancers but not leukemias, which could be directly targeted by circulating cytotoxic NK cells. Our findings uncover a fundamental mechanism that restricts direct NK cell cytotoxicity in peripheral tissues.

Active RB causes visible changes in nuclear organization

RB restricts G1/S progression by inhibiting E2F. Here, we show that sustained expression of active RB, and prolonged G1 arrest, causes visible changes in chromosome architecture that are not directly associated with E2F inhibition. Using FISH probes against two euchromatin RB-associated regions, two heterochromatin domains that lack RB-bound loci, and two whole-chromosome probes, we found that constitutively active RB (ΔCDK-RB) promoted a more diffuse, dispersed, and scattered chromatin organization. These changes were RB dependent, were driven by specific isoforms of monophosphorylated RB, and required known RB-associated activities. ΔCDK-RB altered physical interactions between RB-bound genomic loci, but the RB-induced changes in chromosome architecture were unaffected by dominant-negative DP1. The RB-induced changes appeared to be widespread and influenced chromosome localization within nuclei. Gene expression profiles revealed that the dispersion phenotype was associated with an increased autophagy response. We infer that, after cell cycle arrest, RB acts through noncanonical mechanisms to significantly change nuclear organization, and this reorganization correlates with transitions in cellular state.

NR4A1 regulates expression of immediate early genes, suppressing replication stress in cancer

Deregulation of oncogenic signals in cancer triggers replication stress. Immediate early genes (IEGs) are rapidly and transiently expressed following stressful signals, contributing to an integrated response. Here, we find that the orphan nuclear receptor NR4A1 localizes across the gene body and 3' UTR of IEGs, where it inhibits transcriptional elongation by RNA Pol II, generating R-loops and accessible chromatin domains. Acute replication stress causes immediate dissociation of NR4A1 and a burst of transcriptionally poised IEG expression. Ectopic expression of NR4A1 enhances tumorigenesis by breast cancer cells, while its deletion leads to massive chromosomal instability and proliferative failure, driven by deregulated expression of its IEG target, FOS. Approximately half of breast and other primary cancers exhibit accessible chromatin domains at IEG gene bodies, consistent with this stress-regulatory pathway. Cancers that have retained this mechanism in adapting to oncogenic replication stress may be dependent on NR4A1 for their proliferation.

Abstract 2473: NR4A1 suppresses cancer replication stress through R-loop-dependent inhibition of immediate early gene transcriptional elongation

Deregulation of oncogenic proliferative signals triggers replication stress in cancer cells to which they must adapt. Immediate early genes (IEGs), identified by their rapid stress-induced transient bursts of expression, are critical to integrating downstream signaling pathways. In studying tumor initiation by patient-derived breast cancer cells, we observed acquisition of open chromatin domains at the genebody and 3'-UTR of IEGs, uniquely across the genome. Through in vivo and in vitro modeling, we show that the IEG and orphan nuclear receptor NR4A1 localizes across multiple IEG genebodies, where it binds to RNA Pol II, arresting transcriptional elongation and generating extensive R-loops and accessible chromatin. Acute stress promptly removes NR4A1 from IEG genebodies, triggering immediate release of their poised transcripts. In breast cancer cells, NR4A1 overexpression increases tumorigenesis; conversely, its deletion leads to uncompensated replication stress, chromosomal instability and mitotic catastrophe, driven by deregulation of its IEG target FOS. A large fraction of primary breast and other cancers exhibit open genebody chromatin at IEGs, consistent with preserved NR4A1 function. Thus, NR4A1 mediates a novel transcriptional elongation checkpoint, unique to stress-induced genes and required for their rapid bursts of expression. Cancers that have retained this mechanism in adapting to chronic replication stress may be dependent on NR4A1 for proliferation.

Citation Format: Hongshan Guo, Gabriel Golczer, Shyamala Maheswaran, Michael S. Lawrence, Daniel A. Haber. NR4A1 suppresses cancer replication stress through R-loop-dependent inhibition of immediate early gene transcriptional elongation [abstract]. In: Proceedings of the American Association for Cancer Research Annual Meeting 2021; 2021 Apr 10-15 and May 17-21. Philadelphia (PA): AACR; Cancer Res 2021;81(13_Suppl):Abstract nr 2473.

Abstract 39: APOBEC3A drives acquired resistance to targeted therapies in non small cell lung cancer

Acquired drug resistance to even the most effective anti cancer targeted therapies remains an unsolved clinical problem. Although many drivers of acquired drug resistance have been identified, the underlying molecular mechanisms shaping tumor evolution during treatment are incompletely understood. The extent to which therapy actively drives tumor evolution by promoting mutagenic processes or simply provides the selective pressure necessary for the outgrowth of drug resistant clones remains an open question. Here, we report that lung cancer targeted therapies commonly used in the clinic induce the expression of cytidine deaminase APOBEC3A (A3A), leading to sustained mutagenesis in drug tolerant cancer cells persisting during therapy. Induction of A3A facilitated the formation of double strand DNA breaks (DSBs) in cycling drug-treated cells, and fully resistant clones that evolved from drug tolerant intermediates exhibited an elevated burden of chromosomal aberrations such as copy number alterations and structural variations. Preventing therapy induced A3A mutagenesis either by gene deletion or RNAi mediated suppression delayed the emergence of drug resistance. Finally, we observed accumulation of A3A mutations in lung cancer patients who developed drug resistance after treatment with sequential targeted therapies. These data suggest that induction of A3A mutagenesis in response to targeted therapy treatment may facilitate the development of acquired resistance in non small cell lung cancer. Thus, suppressing expression or enzymatic activity of A3A may represent a potential therapeutic strategy to prevent or delay acquired resistance to lung cancer targeted therapy.

Citation Format: Hideko Isozaki, Ammal Abbasi, Naveed Nikpour, Adam Langenbucher, Wenjia Su, Marcello Stanzione, Lecia V. Sequist, Rémi Buisson, Michael S. Lawrence, Aaron N. Hata. APOBEC3A drives acquired resistance to targeted therapies in non small cell lung cancer [abstract]. In: Proceedings of the American Association for Cancer Research Annual Meeting 2021; 2021 Apr 10-15 and May 17-21. Philadelphia (PA): AACR; Cancer Res 2021;81(13_Suppl):Abstract nr 39.

Dissection of PIK3CA Aberration for Cervical Adenocarcinoma Outcomes

Simple Summary

There is limited information about genomic markers, especially for cervical adenocarcinoma treatment decisions. In this prospective study, it was found that nonsynonymous PIK3CA mutation detected in the patient’s circulating DNA collected before treatment or during follow-up was significantly associated with decreased progression-free survival or overall survival. It is the first indication of the predictive power of PIK3CA aberration in cervical adenocarcinoma. The work contributes to the development of liquid biopsies for the prolonged strategy of surveillance and indicates the possibility of tailoring management of this particular women’s cancer.

Abstract

Personalized treatment of genetically stratified subgroups has the potential to improve outcomes in many malignant tumors. This study distills clinically meaningful prognostic/predictive genomic marker for cervical adenocarcinoma using signature genomic aberrations and single-point nonsynonymous mutation-specific droplet digital PCR (ddPCR). Mutations in PIK3CA E542K, E545K, or H1047R were detected in 41.7% of tumors. PIK3CA mutation detected in the patient’s circulating DNA collected before treatment or during follow-up was significantly associated with decreased progression-free survival or overall survival. PIK3CA mutation in the circulating DNA during follow-up after treatment predicted recurrence with 100% sensitivity and 64.29% specificity. It is the first indication of the predictive power of PIK3CA mutations in cervical adenocarcinoma. The work contributes to the development of liquid biopsies for follow up surveillance and a possibility of tailoring management of this particular women’s cancer.

RASAL2 Confers Collateral MEK/EGFR Dependency in Chemoresistant Triple-Negative Breast Cancer

PURPOSE

While chemotherapy remains the standard treatment for triple-negative breast cancer (TNBC), identifying and managing chemoresistant tumors has proven elusive. We sought to discover hallmarks and therapeutically actionable features of refractory TNBC through molecular analysis of primary chemoresistant TNBC specimens.

EXPERIMENTAL DESIGN

We performed transcriptional profiling of tumors from a phase II clinical trial of platinum chemotherapy for advanced TNBC (TBCRC-009), revealing a gene expression signature that identified de novo chemorefractory tumors. We then employed pharmacogenomic data mining, proteomic and other molecular studies to define the therapeutic vulnerabilities of these tumors.

RESULTS

We reveal the RAS-GTPase-activating protein (RAS-GAP) RASAL2 as an upregulated factor that mediates chemotherapy resistance but also an exquisite collateral sensitivity to combination MAP kinase kinase (MEK1/2) and EGFR inhibitors in TNBC. Mechanistically, RASAL2 GAP activity is required to confer kinase inhibitor sensitivity, as RASAL2-high TNBCs sustain basal RAS activity through suppression of negative feedback regulators SPRY1/2, together with EGFR upregulation. Consequently, RASAL2 expression results in failed feedback compensation upon co-inhibition of MEK1/2 and EGFR that induces synergistic apoptosis in vitro and in vivo. In patients with TNBC, high RASAL2 levels predict clinical chemotherapy response and long-term outcomes, and are associated via direct transcriptional regulation with activated oncogenic Yes-Associated Protein (YAP). Accordingly, chemorefractory patient-derived TNBC models exhibit YAP activation, high RASAL2 expression, and tumor regression in response to MEK/EGFR inhibitor combinations despite well-tolerated intermittent dosing.

CONCLUSIONS

These findings identify RASAL2 as a mediator of TNBC chemoresistance that rewires MAPK feedback and cross-talk to confer profound collateral sensitivity to combination MEK1/2 and EGFR inhibitors.

An extended APOBEC3A mutation signature in cancer

APOBEC mutagenesis, a major driver of cancer evolution, is known for targeting TpC sites in DNA. Recently, we showed that APOBEC3A (A3A) targets DNA hairpin loops. Here, we show that DNA secondary structure is in fact an orthogonal influence on A3A substrate optimality and, surprisingly, can override the TpC sequence preference. VpC (non-TpC) sites in optimal hairpins can outperform TpC sites as mutational hotspots. This expanded understanding of APOBEC mutagenesis illuminates the genomic Twin Paradox, a puzzling pattern of closely spaced mutation hotspots in cancer genomes, in which one is a canonical TpC site but the other is a VpC site, and double mutants are seen only in trans, suggesting a two-hit driver event. Our results clarify this paradox, revealing that both hotspots in these twins are optimal A3A substrates. Our findings reshape the notion of a mutation signature, highlighting the additive roles played by DNA sequence and DNA structure.

Spectrum of Mechanisms of Resistance to Crizotinib and Lorlatinib in ROS1 Fusion-Positive Lung Cancer

PURPOSE

Current standard initial therapy for advanced, ROS proto-oncogene 1, receptor tyrosine kinase fusion (ROS1)-positive (ROS1⁺) non-small cell lung cancer (NSCLC) is crizotinib or entrectinib. Lorlatinib, a next-generation anaplastic lymphoma kinase/ROS1 inhibitor, recently demonstrated efficacy in ROS1⁺ NSCLC, including in crizotinib-pretreated patients. However, mechanisms of lorlatinib resistance in ROS1⁺ disease remain poorly understood. Here, we assessed mechanisms of resistance to crizotinib and lorlatinib.

EXPERIMENTAL DESIGN

Biopsies from patients with ROS1 ⁺ NSCLC progressing on crizotinib or lorlatinib were profiled by genetic sequencing.

RESULTS

From 55 patients, 47 post-crizotinib and 32 post-lorlatinib biopsies were assessed. Among 42 post-crizotinib and 28 post-lorlatinib biopsies analyzed at distinct timepoints, ROS1 mutations were identified in 38% and 46%, respectively. ROS1 G2032R was the most commonly occurring mutation in approximately one third of cases. Additional ROS1 mutations included D2033N (2.4%) and S1986F (2.4%) post-crizotinib and L2086F (3.6%), G2032R/L2086F (3.6%), G2032R/S1986F/L2086F (3.6%), and S1986F/L2000V (3.6%) post-lorlatinib. Structural modeling predicted ROS1^L2086F causes steric interference to lorlatinib, crizotinib, and entrectinib, while it may accommodate cabozantinib. In Ba/F3 models, ROS1^L2086F, ROS1^{G2032R/L2086F}, and ROS1^{S1986F/G2032R/L2086F} were refractory to lorlatinib but sensitive to cabozantinib. A patient with disease progression on crizotinib and lorlatinib and ROS1 L2086F received cabozantinib for nearly 11 months with disease control. Among lorlatinib-resistant biopsies, we also identified MET amplification (4%), KRAS G12C (4%), KRAS amplification (4%), NRAS mutation (4%), and MAP2K1 mutation (4%).

CONCLUSIONS

ROS1 mutations mediate resistance to crizotinib and lorlatinib in more than one third of cases, underscoring the importance of developing next-generation ROS1 inhibitors with potency against these mutations, including G2032R and L2086F. Continued efforts are needed to elucidate ROS1-independent resistance mechanisms.

HPV+ oropharyngeal squamous cell carcinomas from patients with two tumors display synchrony of viral genomes yet discordant mutational profiles and signatures

Human papillomavirus (HPV) positive oropharyngeal squamous cell carcinoma (HPV + OPSCC) is increasing in prevalence in the USA, as are cases of patients with multiple HPV + OPSCCs (mHPV + OPSCC). mHPV + OPSCCs present a unique opportunity to examine HPV + OPSCC mutation acquisition and evolution. We performed sequencing of the viral genome, somatic exome and somatic transcriptome from 8 patients each with 2 spatially distinct HPV + OPSCCs, and 37 'traditional' HPV + OPSCCs to first address if paired tumors are caused by the same viral isolate and next, if acquired alterations, and the underlying processes driving mutagenesis, are shared within pairs. All tumor pairs contained viral genomes from the same HPV type 16 sublineage and differed by 0-2 clonal single nucleotide polymorphisms (SNPs), suggesting infection with the same viral isolate. Despite this, there was significant discordance in expression profiles, mutational burden and mutational profiles between tumors in a pair, with only two pairs sharing any overlapping mutations (3/3343 variants). Within tumor pairs there was a striking discrepancy of mutational signatures, exemplified by no paired tumors sharing high APOBEC mutational burden. Here, leveraging mHPV + OPSCCs as a model system to study mutation acquisition in virally mediated tumors, in which the germline, environmental exposures, immune surveillance and tissue/organ type were internally controlled, we demonstrate that despite infection by the same viral isolate, paired mHPV + OPSCCs develop drastically different somatic alterations and even more strikingly, appear to be driven by disparate underlying mutational processes. Thus, despite a common starting point, HPV + OPSCCs evolve through variable mutational processes with resultant stochastic mutational profiles.

Abstract 4893: Multi-modality sequencing of multiple HPV mediated oropharyngeal squamous cell carcinomas from single patients reveals synchrony in viral genomes yet discordant mutational profiles and signatures

Human Papilloma Virus (HPV) mediated oropharyngeal squamous cell carcinoma (HPVmOPSCC) is increasing in prevalence in the United States, as are cases of patients with multiple, spatially distinct, HPVmOPSCCs (mHPVmOPSCC). Patients with mHPVmOPSCCs present a unique opportunity to study the biology of HPV mediated tumorigenesis in the oropharynx, which currently is poorly understood. We performed whole exome sequencing and RNA-Seq on 16 paired HPVmOPSCCs (eight patients) and whole genome sequencing of the viral genome using a HPV16 NGS assay described previously. Viral sequencing revealed considerable diversity between patients, across multiple viral sublineages. However, all tumor pairs from a given patient contained the same viral genome, with at most 2 clonal SNPs distinguishing the viruses within a pair. In contrast, tumor DNA sequencing revealed striking dissimilarity between tumors from a given patient. Seven of eight tumor pairs shared no overlapping somatic mutations with one pair sharing only six variants out of 277 (2%). Mutational signature analysis demonstrated a dominant APOBEC signature in six tumors. Tumor pairs showed striking discordance of mutational signatures. In 3/8 pairs, one tumor exhibited a nearly 100% pure APOBEC+ mutational profile, while the other was totally devoid of APOBEC-signature mutations. Here, we have demonstrated that spatially independent tumors, originating within the same individual, in the same tissue type, are caused by the same virus. Interesting, despite this shared etiology, the pair of tumors within a patient develop with vastly different mutational profiles and have different mutational processes driving tumorigenesis. This suggests that despite infection by a single oncogenic virus, HPVmOPSCCs may develop via multiple, disparate pathways. Of particular interest is the tumor pairs' clear discordance in APOBEC mutational burden, which challenges existing hypotheses on the role of APOBEC mutagenesis in this disease.

Citation Format: Daniel L. Faden, Adam Langenbucher, Krystle Kuhs, James S. Lewis, Michael S. Lawrence, Lisa Mirabello, Robert L. Ferris. Multi-modality sequencing of multiple HPV mediated oropharyngeal squamous cell carcinomas from single patients reveals synchrony in viral genomes yet discordant mutational profiles and signatures [abstract]. In: Proceedings of the Annual Meeting of the American Association for Cancer Research 2020; 2020 Apr 27-28 and Jun 22-24. Philadelphia (PA): AACR; Cancer Res 2020;80(16 Suppl):Abstract nr 4893.

Small cell transformation of ROS1 fusion-positive lung cancer resistant to ROS1 inhibition

Histologic transformation from non-small cell to small cell lung cancer has been reported as a resistance mechanism to targeted therapy in EGFR-mutant and ALK fusion-positive lung cancers. Whether small cell transformation occurs in other oncogene-driven lung cancers remains unknown. Here we analyzed the genomic landscape of two pre-mortem and 11 post-mortem metastatic tumors collected from an advanced, ROS1 fusion-positive lung cancer patient, who had received sequential ROS1 inhibitors. Evidence of small cell transformation was observed in all metastatic sites at autopsy, with inactivation of RB1 and TP53, and loss of ROS1 fusion expression. Whole-exome sequencing revealed minimal mutational and copy number heterogeneity, suggestive of “hard” clonal sweep. Patient-derived models generated from autopsy retained features consistent with small cell lung cancer and demonstrated resistance to ROS1 inhibitors. This case supports small cell transformation as a recurring resistance mechanism, and underscores the importance of elucidating its biology to expand therapeutic opportunities.

Identification of Somatically Acquired BRCA1/2 Mutations by cfDNA Analysis in Patients with Metastatic Breast Cancer

PURPOSE

Plasma genotyping may identify mutations in potentially "actionable" cancer genes, such as BRCA1/2, but their clinical significance is not well-defined. We evaluated the characteristics of somatically acquired BRCA1/2 mutations in patients with metastatic breast cancer (MBC).

EXPERIMENTAL DESIGN

Patients with MBC undergoing routine cell-free DNA (cfDNA) next-generation sequencing (73-gene panel) before starting a new therapy were included. Somatic BRCA1/2 mutations were classified as known germline pathogenic mutations or novel variants, and linked to clinicopathologic characteristics. The effect of the PARP inhibitor, olaparib, was assessed in vitro, using cultured circulating tumor cells (CTCs) from a patient with a somatically acquired BRCA1 mutation and a second patient with an acquired BRCA2 mutation.

RESULTS

Among 215 patients with MBC, 29 (13.5%) had somatic cfDNA BRCA1/2 mutations [nine (4%) known germline pathogenic and rest (9%) novel variants]. Known germline pathogenic BRCA1/2 mutations were common in younger patients (P = 0.008), those with triple-negative disease (P = 0.022), and they were more likely to be protein-truncating alterations and be associated with TP53 mutations. Functional analysis of a CTC culture harboring a somatic BRCA1 mutation demonstrated high sensitivity to PARP inhibition, while another CTC culture harboring a somatic BRCA2 mutation showed no differential sensitivity. Across the entire cohort, APOBEC mutational signatures (COSMIC Signatures 2 and 13) and the "BRCA" mutational signature (COSMIC Signature 3) were present in BRCA1/2-mutant and wild-type cases, demonstrating the high mutational burden associated with advanced MBC.

CONCLUSIONS

Somatic BRCA1/2 mutations are readily detectable in MBC by cfDNA analysis, and may be present as both known germline pathogenic and novel variants.

MET Alterations Are a Recurring and Actionable Resistance Mechanism in ALK-Positive Lung Cancer

PURPOSE

Most ALK-positive lung cancers will develop ALK-independent resistance after treatment with next-generation ALK inhibitors. MET amplification has been described in patients progressing on ALK inhibitors, but frequency of this event has not been comprehensively assessed.

EXPERIMENTAL DESIGN

We performed FISH and/or next-generation sequencing on 207 posttreatment tissue (n = 101) or plasma (n = 106) specimens from patients with ALK-positive lung cancer to detect MET genetic alterations. We evaluated ALK inhibitor sensitivity in cell lines with MET alterations and assessed antitumor activity of ALK/MET blockade in ALK-positive cell lines and 2 patients with MET-driven resistance.

RESULTS

MET amplification was detected in 15% of tumor biopsies from patients relapsing on next-generation ALK inhibitors, including 12% and 22% of biopsies from patients progressing on second-generation inhibitors or lorlatinib, respectively. Patients treated with a second-generation ALK inhibitor in the first-line setting were more likely to develop MET amplification than those who had received next-generation ALK inhibitors after crizotinib (P = 0.019). Two tumor specimens harbored an identical ST7-MET rearrangement, one of which had concurrent MET amplification. Expressing ST7-MET in the sensitive H3122 ALK-positive cell line induced resistance to ALK inhibitors that was reversed with dual ALK/MET inhibition. MET inhibition resensitized a patient-derived cell line harboring both ST7-MET and MET amplification to ALK inhibitors. Two patients with ALK-positive lung cancer and acquired MET alterations achieved rapid responses to ALK/MET combination therapy.

CONCLUSIONS

Treatment with next-generation ALK inhibitors, particularly in the first-line setting, may lead to MET-driven resistance. Patients with acquired MET alterations may derive clinical benefit from therapies that target both ALK and MET.

A post-transcriptional program of chemoresistance by AU-rich elements and TTP in quiescent leukemic cells

BACKGROUND

Quiescence (G0) is a transient, cell cycle-arrested state. By entering G0, cancer cells survive unfavorable conditions such as chemotherapy and cause relapse. While G0 cells have been studied at the transcriptome level, how post-transcriptional regulation contributes to their chemoresistance remains unknown.

RESULTS

We induce chemoresistant and G0 leukemic cells by serum starvation or chemotherapy treatment. To study post-transcriptional regulation in G0 leukemic cells, we systematically analyzed their transcriptome, translatome, and proteome. We find that our resistant G0 cells recapitulate gene expression profiles of in vivo chemoresistant leukemic and G0 models. In G0 cells, canonical translation initiation is inhibited; yet we find that inflammatory genes are highly translated, indicating alternative post-transcriptional regulation. Importantly, AU-rich elements (AREs) are significantly enriched in the upregulated G0 translatome and transcriptome. Mechanistically, we find the stress-responsive p38 MAPK-MK2 signaling pathway stabilizes ARE mRNAs by phosphorylation and inactivation of mRNA decay factor, Tristetraprolin (TTP) in G0. This permits expression of ARE mRNAs that promote chemoresistance. Conversely, inhibition of TTP phosphorylation by p38 MAPK inhibitors and non-phosphorylatable TTP mutant decreases ARE-bearing TNFα and DUSP1 mRNAs and sensitizes leukemic cells to chemotherapy. Furthermore, co-inhibiting p38 MAPK and TNFα prior to or along with chemotherapy substantially reduces chemoresistance in primary leukemic cells ex vivo and in vivo.

CONCLUSIONS

These studies uncover post-transcriptional regulation underlying chemoresistance in leukemia. Our data reveal the p38 MAPK-MK2-TTP axis as a key regulator of expression of ARE-bearing mRNAs that promote chemoresistance. By disrupting this pathway, we develop an effective combination therapy against chemosurvival.

Pathway and network analysis of more than 2500 whole cancer genomes

The catalog of cancer driver mutations in protein-coding genes has greatly expanded in the past decade. However, non-coding cancer driver mutations are less well-characterized and only a handful of recurrent non-coding mutations, most notably TERT promoter mutations, have been reported. Here, as part of the ICGC/TCGA Pan-Cancer Analysis of Whole Genomes (PCAWG) Consortium, which aggregated whole genome sequencing data from 2658 cancer across 38 tumor types, we perform multi-faceted pathway and network analyses of non-coding mutations across 2583 whole cancer genomes from 27 tumor types compiled by the ICGC/TCGA PCAWG project that was motivated by the success of pathway and network analyses in prioritizing rare mutations in protein-coding genes. While few non-coding genomic elements are recurrently mutated in this cohort, we identify 93 genes harboring non-coding mutations that cluster into several modules of interacting proteins. Among these are promoter mutations associated with reduced mRNA expression in TP53, TLE4, and TCF4. We find that biological processes had variable proportions of coding and non-coding mutations, with chromatin remodeling and proliferation pathways altered primarily by coding mutations, while developmental pathways, including Wnt and Notch, altered by both coding and non-coding mutations. RNA splicing is primarily altered by non-coding mutations in this cohort, and samples containing non-coding mutations in well-known RNA splicing factors exhibit similar gene expression signatures as samples with coding mutations in these genes. These analyses contribute a new repertoire of possible cancer genes and mechanisms that are altered by non-coding mutations and offer insights into additional cancer vulnerabilities that can be investigated for potential therapeutic treatments.

Pan-cancer analysis of whole genomes

Cancer is driven by genetic change, and the advent of massively parallel sequencing has enabled systematic documentation of this variation at the whole-genome scale1-3. Here we report the integrative analysis of 2,658 whole-cancer genomes and their matching normal tissues across 38 tumour types from the Pan-Cancer Analysis of Whole Genomes (PCAWG) Consortium of the International Cancer Genome Consortium (ICGC) and The Cancer Genome Atlas (TCGA). We describe the generation of the PCAWG resource, facilitated by international data sharing using compute clouds. On average, cancer genomes contained 4-5 driver mutations when combining coding and non-coding genomic elements; however, in around 5% of cases no drivers were identified, suggesting that cancer driver discovery is not yet complete. Chromothripsis, in which many clustered structural variants arise in a single catastrophic event, is frequently an early event in tumour evolution; in acral melanoma, for example, these events precede most somatic point mutations and affect several cancer-associated genes simultaneously. Cancers with abnormal telomere maintenance often originate from tissues with low replicative activity and show several mechanisms of preventing telomere attrition to critical levels. Common and rare germline variants affect patterns of somatic mutation, including point mutations, structural variants and somatic retrotransposition. A collection of papers from the PCAWG Consortium describes non-coding mutations that drive cancer beyond those in the TERT promoter4; identifies new signatures of mutational processes that cause base substitutions, small insertions and deletions and structural variation5,6; analyses timings and patterns of tumour evolution7; describes the diverse transcriptional consequences of somatic mutation on splicing, expression levels, fusion genes and promoter activity8,9; and evaluates a range of more-specialized features of cancer genomes8,10-18.

Analyses of non-coding somatic drivers in 2,658 cancer whole genomes

The discovery of drivers of cancer has traditionally focused on protein-coding genes1-4. Here we present analyses of driver point mutations and structural variants in non-coding regions across 2,658 genomes from the Pan-Cancer Analysis of Whole Genomes (PCAWG) Consortium5 of the International Cancer Genome Consortium (ICGC) and The Cancer Genome Atlas (TCGA). For point mutations, we developed a statistically rigorous strategy for combining significance levels from multiple methods of driver discovery that overcomes the limitations of individual methods. For structural variants, we present two methods of driver discovery, and identify regions that are significantly affected by recurrent breakpoints and recurrent somatic juxtapositions. Our analyses confirm previously reported drivers6,7, raise doubts about others and identify novel candidates, including point mutations in the 5' region of TP53, in the 3' untranslated regions of NFKBIZ and TOB1, focal deletions in BRD4 and rearrangements in the loci of AKR1C genes. We show that although point mutations and structural variants that drive cancer are less frequent in non-coding genes and regulatory sequences than in protein-coding genes, additional examples of these drivers will be found as more cancer genomes become available.

Histone Lysine Methylation Dynamics Control EGFR DNA Copy-Number Amplification

Acquired chromosomal DNA copy gains are a feature of many tumors; however, the mechanisms that underpin oncogene amplification are poorly understood. Recent studies have begun to uncover the importance of epigenetic states and histone lysine methyltransferases (KMT) and demethylases (KDM) in regulating transient site-specific DNA copy-number gains (TSSG). In this study, we reveal a critical interplay between a myriad of lysine methyltransferases and demethylases in modulating H3K4/9/27 methylation balance to control extrachromosomal amplification of the EGFR oncogene. This study further establishes that cellular signals (hypoxia and EGF) are able to directly promote EGFR amplification through modulation of the enzymes controlling EGFR copy gains. Moreover, we demonstrate that chemical inhibitors targeting specific KMTs and KDMs are able to promote or block extrachromosomal EGFR amplification, which identifies potential therapeutic strategies for controlling EGFR copy-number heterogeneity in cancer, and, in turn, drug response. SIGNIFICANCE: This study identifies a network of epigenetic factors and cellular signals that directly control EGFR DNA amplification. We demonstrate that chemical inhibitors targeting enzymes controlling this amplification can be used to rheostat EGFR copy number, which uncovers therapeutic opportunities for controlling EGFR DNA amplification heterogeneity and the associated drug response.This article is highlighted in the In This Issue feature, p. 161.

Combination Olaparib and Temozolomide in Relapsed Small-Cell Lung Cancer

Small-cell lung cancer (SCLC) is an aggressive malignancy in which inhibitors of PARP have modest single-agent activity. We performed a phase I/II trial of combination olaparib tablets and temozolomide (OT) in patients with previously treated SCLC. We established a recommended phase II dose of olaparib 200 mg orally twice daily with temozolomide 75 mg/m² daily, both on days 1 to 7 of a 21-day cycle, and expanded to a total of 50 patients. The confirmed overall response rate was 41.7% (20/48 evaluable); median progression-free survival was 4.2 months [95% confidence interval (CI), 2.8-5.7]; and median overall survival was 8.5 months (95% CI, 5.1-11.3). Patient-derived xenografts (PDX) from trial patients recapitulated clinical OT responses, enabling a 32-PDX coclinical trial. This revealed a correlation between low basal expression of inflammatory-response genes and cross-resistance to both OT and standard first-line chemotherapy (etoposide/platinum). These results demonstrate a promising new therapeutic strategy in SCLC and uncover a molecular signature of those tumors most likely to respond. SIGNIFICANCE: We demonstrate substantial clinical activity of combination olaparib/temozolomide in relapsed SCLC, revealing a promising new therapeutic strategy for this highly recalcitrant malignancy. Through an integrated coclinical trial in PDXs, we then identify a molecular signature predictive of response to OT, and describe the common molecular features of cross-resistant SCLC.See related commentary by Pacheco and Byers, p. 1340.This article is highlighted in the In This Issue feature, p. 1325.

GTP-Dependent Formation of Multimeric G-Quadruplexes

G-Quadruplexes are noncanonical nucleic acid structures made up of stacked guanosine tetrads connected by short loops. They are frequently used building blocks in synthetic biology and thought to play widespread biological roles. Multimerization can change the functional properties of G-quadruplexes, and understanding the factors that modulate this process remains an important goal. Here, we report the discovery of a novel mechanism by which the formation of multimeric G-quadruplexes can be controlled using GTP. We show that GTP likely inhibits multimer formation by becoming incorporated into a tetrad in the monomeric form of the structure and define the sequence requirements of G-quadruplexes that form GTP-dependent structures. These experiments provide new insights into the small molecule control of G-quadruplex multimerization. They also suggest possible roles for GTP-dependent multimeric G-quadruplexes in both synthetic and natural biological systems.