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.

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 preclinical 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.

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. This article is featured in Selected Articles from This Issue, p. 695.

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.

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.

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.

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.

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.