Understanding Mechanisms of Response to CAR T-cell Therapy through Single-Cell Sequencing: Insights and Challenges

SUMMARY

Single-cell RNA sequencing has emerged as a powerful technique to understand the molecular features of chimeric antigen receptor (CAR) T cells that associate with clinical outcomes. Here we discuss the common themes that have emerged from across single-cell studies of CAR T-cell therapy, and summarize the challenges in interpreting this complex data type.

Latent human herpesvirus 6 is reactivated in CAR T cells

Cell therapies have yielded durable clinical benefits for patients with cancer, but the risks associated with the development of therapies from manipulated human cells are understudied. For example, we lack a comprehensive understanding of the mechanisms of toxicities observed in patients receiving T cell therapies, including recent reports of encephalitis caused by reactivation of human herpesvirus 6 (HHV-6)1. Here, through petabase-scale viral genomics mining, we examine the landscape of human latent viral reactivation and demonstrate that HHV-6B can become reactivated in cultures of human CD4+ T cells. Using single-cell sequencing, we identify a rare population of HHV-6 'super-expressors' (about 1 in 300-10,000 cells) that possess high viral transcriptional activity, among research-grade allogeneic chimeric antigen receptor (CAR) T cells. By analysing single-cell sequencing data from patients receiving cell therapy products that are approved by the US Food and Drug Administration2 or are in clinical studies3-5, we identify the presence of HHV-6-super-expressor CAR T cells in patients in vivo. Together, the findings of our study demonstrate the utility of comprehensive genomics analyses in implicating cell therapy products as a potential source contributing to the lytic HHV-6 infection that has been reported in clinical trials1,6-8 and may influence the design and production of autologous and allogeneic cell therapies.

Pan-cancer analysis of post-translational modifications reveals shared patterns of protein regulation

Post-translational modifications (PTMs) play key roles in regulating cell signaling and physiology in both normal and cancer cells. Advances in mass spectrometry enable high-throughput, accurate, and sensitive measurement of PTM levels to better understand their role, prevalence, and crosstalk. Here, we analyze the largest collection of proteogenomics data from 1,110 patients with PTM profiles across 11 cancer types (10 from the National Cancer Institute's Clinical Proteomic Tumor Analysis Consortium [CPTAC]). Our study reveals pan-cancer patterns of changes in protein acetylation and phosphorylation involved in hallmark cancer processes. These patterns revealed subsets of tumors, from different cancer types, including those with dysregulated DNA repair driven by phosphorylation, altered metabolic regulation associated with immune response driven by acetylation, affected kinase specificity by crosstalk between acetylation and phosphorylation, and modified histone regulation. Overall, this resource highlights the rich biology governed by PTMs and exposes potential new therapeutic avenues.

Systematic characterization of therapeutic vulnerabilities in Multiple Myeloma with Amp1q reveals increased sensitivity to the combination of MCL1 and PI3K inhibitors

The development of targeted therapy for patients with Multiple Myeloma (MM) is hampered by the low frequency of actionable genetic abnormalities. Gain or amplification of chr1q (Amp1q) is the most frequent arm-level copy number gain in patients with MM, and it is associated with higher risk of progression and death despite recent advances in therapeutics. Thus, developing targeted therapy for patients with MM and Amp1q stands to benefit a large portion of patients in need of more effective management. Here, we employed large-scale dependency screens and drug screens to systematically characterize the therapeutic vulnerabilities of MM with Amp1q and showed increased sensitivity to the combination of MCL1 and PI3K inhibitors. Using single-cell RNA sequencing, we compared subclones with and without Amp1q within the same patient tumors and showed that Amp1q is associated with higher levels of MCL1 and the PI3K pathway. Furthermore, by isolating isogenic clones with different copy number for part of the chr1q arm, we showed increased sensitivity to MCL1 and PI3K inhibitors with arm-level gain. Lastly, we demonstrated synergy between MCL1 and PI3K inhibitors and dissected their mechanism of action in MM with Amp1q.

Abstract 960: Overlapping expression landscape of antibody drug conjugate targets, trophoblast cell surface antigen 2 (Trop-2) & human epidermal growth factor receptor 2 (HER2), in breast cancer

Background Antibody drug conjugates (ADC) are novel drugs linking potent payloads to antibodies targeting antigen-expressing tumors. Sacituzumab govitecan (SG), targeting Trop-2, is approved for metastatic triple negative breast cancer (TNBC); and trastuzumab deruxtecan, targeting HER2, is approved for HER2-positive and HER2-low metastatic breast cancer. To understand the potentially overlapping clinical landscape of Trop-2 and HER2 antigens, we evaluated RNA expression data in breast cancer from The Cancer Genome Atlas (TCGA) project.

Methods TCGA dataset was assessed for Trop-2 and HER2 expression via processed RNA sequencing (RNA-seq) data of the corresponding genes TACSTD2 and ERBB2. Medium/high gene expression was assessed as >100 transcripts per million (TPM). Samples were classified HER2-low per ASCO/CAP guidelines. Gene expression across clinical parameters was assessed via one-way ANOVA.

Results 1076 patients with primary breast cancer were included. The majority (59%) had both high TACSTD2 expression (TACSTD2hi) and high ERBB2 expression (ERBB2hi) (see Table). Median TACSTD2 expression was 572 TPM (IQR 349-666 TPM); median ERBB2 expression was 122 TPM (IQR 73-192 TPM). No significant difference was observed in TACSTD2 or ERBB2 expression among invasive ductal carcinoma, invasive lobular carcinoma, mixed histology, or other (p = 0.07, 0.23). No significant difference in TACSTD2 expression was noted between HER2-low and HER2-negative subtypes (p=0.34).

Conclusions While SG is approved in TNBC, TACSTD2 is expressed across all breast cancer subtypes, including HER2-low, suggesting a broader population may benefit from Trop-2-targeted ADCs. Furthermore, given that over half of breast cancers have high expression of both TACSTD2 and ERBB2, additional studies are needed to understand the optimal sequencing of ADC-based therapies for patients with breast cancer.

Table Patient subsets(HR = hormone receptor) n TACSTD2hi & ERBB2hi TACSTD2hi & ERBB2low TACSTD2low & ERBB2hi TACSTD2low & ERBB2low All 1,076 632 (59%) 394 (36%) 19 (2%) 31 (3%) Histology 843 Invasive ductal carcinoma 507 267 (32%) 215 (25%) 10 (1%) 15 (2%) Invasive lobular carcinoma 130 103 (12%) 24 (3%) 1 (<1%) 2 (<1%) Mixed 91 63 (7%) 24 (3%) 2 (<1%) 2 (<1%) Other 115 65 (8%) 43 (5%) 2 (<1%) 5 (<1%) HER2 Status 345 HR+/HER2-low 254 170 (49%) 77 (22%) 6 (1%) 1 (<1%) TNBC/HER2-low 55 12 (3%) 39 (11%) 1 (<1%) 3 (<1%) HR+/HER2-negative 24 13 (4%) 11 (3%) 0 0 TNBC/HER2-negative 12 1 (<1%) 10 0 1 (<1%)

Citation Format: Arielle J. Medford, Nicholas J. Haradhvala, Neelima Vidula, Rachel Abelman, Laura M. Spring, Leif W. Ellisen, Gad Getz, Aditya Bardia. Overlapping expression landscape of antibody drug conjugate targets, trophoblast cell surface antigen 2 (Trop-2) & human epidermal growth factor receptor 2 (HER2), in breast cancer [abstract]. In: Proceedings of the American Association for Cancer Research Annual Meeting 2023; Part 1 (Regular and Invited Abstracts); 2023 Apr 14-19; Orlando, FL. Philadelphia (PA): AACR; Cancer Res 2023;83(7_Suppl):Abstract nr 960.

Immune biomarkers of response to immunotherapy in patients with high-risk smoldering myeloma

Patients with smoldering multiple myeloma (SMM) are observed until progression, but early treatment may improve outcomes. We conducted a phase II trial of elotuzumab, lenalidomide, and dexamethasone (EloLenDex) in patients with high-risk SMM and performed single-cell RNA and T cell receptor (TCR) sequencing on 149 bone marrow (BM) and peripheral blood (PB) samples from patients and healthy donors (HDs). We find that early treatment with EloLenDex is safe and effective and provide a comprehensive characterization of alterations in immune cell composition and TCR repertoire diversity in patients. We show that the similarity of a patient's immune cell composition to that of HDs may have prognostic relevance at diagnosis and after treatment and that the abundance of granzyme K (GZMK)+ CD8+ effector memory T (TEM) cells may be associated with treatment response. Last, we uncover similarities between immune alterations observed in the BM and PB, suggesting that PB-based immune profiling may have diagnostic and prognostic utility.

Designing sensitive viral diagnostics with machine learning

Design of nucleic acid-based viral diagnostics typically follows heuristic rules and, to contend with viral variation, focuses on a genome's conserved regions. A design process could, instead, directly optimize diagnostic effectiveness using a learned model of sensitivity for targets and their variants. Toward that goal, we screen 19,209 diagnostic-target pairs, concentrated on CRISPR-based diagnostics, and train a deep neural network to accurately predict diagnostic readout. We join this model with combinatorial optimization to maximize sensitivity over the full spectrum of a virus's genomic variation. We introduce Activity-informed Design with All-inclusive Patrolling of Targets (ADAPT), a system for automated design, and use it to design diagnostics for 1,933 vertebrate-infecting viral species within 2 hours for most species and within 24 hours for all but three. We experimentally show that ADAPT's designs are sensitive and specific to the lineage level and permit lower limits of detection, across a virus's variation, than the outputs of standard design techniques. Our strategy could facilitate a proactive resource of assays for detecting pathogens.

Abstract 3575: Differential dynamics of response at single cell resolution between axi-cel and tisa-cel CAR-T therapy in refractory B-cell lymphomas

Chimeric Antigen Receptor (CAR)-T cell therapy has revolutionized the treatment of hematologic malignancies. Approximately half of patients with refractory large B-cell lymphomas achieve durable responses from CD19-targeting CAR-T treatment, across the commercially available CAR-T products with differing designs. Known failure mechanisms such as antigen loss account for only a fraction of cases without durable responses, and this knowledge gap has limited advances in CAR-T engineering and optimal targeting to patients.

We hypothesized that characterization of the transcriptional programs and temporal evolution of CAR-T and host immune cell populations could provide novel insights into the basis of clinical response to CAR-T cell therapy for B cell lymphoma. We performed 10X single-cell RNA sequencing on serial samples collected from 32 individuals with high grade B cell lymphoma treated with the two first FDA-approved CD19 CAR-T products: axicabtagene ciloleucel (axi-cel, utilizing a CD28z costimulatory domain) and tisagenlecleucel (tisa-cel, with a 4-1BB domain). We analyzed 106 samples, including pre-infusion blood samples, infusion product, and post-infusion T cells sorted by flow cytometry into CAR+ and CAR- populations. Analyzing 602,577 single-cell transcriptomes, we discerned major differences in the dynamics of response of the two products. Tisa-cel responders showed dramatic expansion of CD8+ T cells at day 7 after infusion, which represented less than 10% of cells in the product. Conversely, CD8+ T cells in products of non-responders failed to expand to the same degree post-infusion and had a more effector- than memory-like T cell phenotype. In one tisa-cel–treated patient who had no CD8+ T cell expansion after initial infusion and relapsed at 6 months post-infusion, re-treatment with a second dose of the same product led to a durable response and was associated with greater CD8+ T cell expansion as well as a shift in CD4+ T phenotype from cytotoxic to helper. In contrast, axi-cel responders had pre-expanded effector populations distributed more heterogeneously among CD4+ and CD8+ T cells. Finally, we identified nominal elevations in CAR-T regulatory cells (CAR-Tregs) among both axi-cel and tisa-cel non-responders in our dataset, which we confirmed in an external dataset. These small increases in CAR-Tregs were sufficient to uniformly suppress conventional CAR-T cell expansion and drive late relapses in an in vivo mouse model of lymphoma after treatment with CARs with either CD28z or 4-1BB co-stimulatory domains. In summary, this represents the largest CAR-T scRNAseq cohort established thus far and provides important insights into (i) the temporal dynamics of a successful CAR-T response, (ii) the molecular phenotypes of CAR-T cells with different costimulatory domains, and (iii) the capacity for small increases in CAR-Tregs to drive relapse.

Citation Format: Nicholas J. Haradhvala, Mark B. Leick, Katie Maurer, Satyen Gohil, Rebecca C. Larson, Estelle Yao, Matthew J. Frigault, Shuqiang Li, Kenneth J. Livak, Kahn Rhrissorrakrai, Filippo Utro, Chaya Levovitz, Raquel A. Jacobs, Kara Slowik, Brian P. Danysh, Laxmi Parida, Catherine J. Wu, Gad Getz, Marcela V. Maus. Differential dynamics of response at single cell resolution between axi-cel and tisa-cel CAR-T therapy in refractory B-cell lymphomas [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 3575.

CAR T cell killing requires the IFNγR pathway in solid but not liquid tumours

Chimeric antigen receptor (CAR) therapy has had a transformative effect on the treatment of haematologic malignancies^1-6, but it has shown limited efficacy against solid tumours. Solid tumours may have cell-intrinsic resistance mechanisms to CAR T cell cytotoxicity. Here, to systematically identify potential resistance pathways in an unbiased manner, we conducted a genome-wide CRISPR knockout screen in glioblastoma, a disease in which CAR T cells have had limited efficacy^7,8. We found that the loss of genes in the interferon-γ receptor (IFNγR) signalling pathway (IFNGR1, JAK1 or JAK2) rendered glioblastoma and other solid tumours more resistant to killing by CAR T cells both in vitro and in vivo. However, loss of this pathway did not render leukaemia or lymphoma cell lines insensitive to CAR T cells. Using transcriptional profiling, we determined that glioblastoma cells lacking IFNγR1 had lower upregulation of cell-adhesion pathways after exposure to CAR T cells. We found that loss of IFNγR1 in glioblastoma cells reduced overall CAR T cell binding duration and avidity. The critical role of IFNγR signalling in susceptibility of solid tumours to CAR T cells is surprising, given that CAR T cells do not require traditional antigen-presentation pathways. Instead, in glioblastoma tumours, IFNγR signalling was required for sufficient adhesion of CAR T cells to mediate productive cytotoxicity. Our work demonstrates that liquid and solid tumours differ in their interactions with CAR T cells and suggests that enhancing binding interactions between T cells and tumour cells may yield improved responses in solid tumours.

DNA Polymerase and Mismatch Repair Exert Distinct Microsatellite Instability Signatures in Normal and Malignant Human Cells

Although replication repair deficiency, either by mismatch repair deficiency (MMRD) and/or loss of DNA polymerase proofreading, can cause hypermutation in cancer, microsatellite instability (MSI) is considered a hallmark of MMRD alone. By genome-wide analysis of tumors with germline and somatic deficiencies in replication repair, we reveal a novel association between loss of polymerase proofreading and MSI, especially when both components are lost. Analysis of indels in microsatellites (MS-indels) identified five distinct signatures (MS-sigs). MMRD MS-sigs are dominated by multibase losses, whereas mutant-polymerase MS-sigs contain primarily single-base gains. MS deletions in MMRD tumors depend on the original size of the MS and converge to a preferred length, providing mechanistic insight. Finally, we demonstrate that MS-sigs can be a powerful clinical tool for managing individuals with germline MMRD and replication repair-deficient cancers, as they can detect the replication repair deficiency in normal cells and predict their response to immunotherapy. SIGNIFICANCE: Exome- and genome-wide MSI analysis reveals novel signatures that are uniquely attributed to mismatch repair and DNA polymerase. This provides new mechanistic insight into MS maintenance and can be applied clinically for diagnosis of replication repair deficiency and immunotherapy response prediction.This article is highlighted in the In This Issue feature, p. 995.

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.

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.

The repertoire of mutational signatures in human cancer

Somatic mutations in cancer genomes are caused by multiple mutational processes, each of which generates a characteristic mutational signature¹. Here, as part of the Pan-Cancer Analysis of Whole Genomes (PCAWG) Consortium² of the International Cancer Genome Consortium (ICGC) and The Cancer Genome Atlas (TCGA), we characterized mutational signatures using 84,729,690 somatic mutations from 4,645 whole-genome and 19,184 exome sequences that encompass most types of cancer. We identified 49 single-base-substitution, 11 doublet-base-substitution, 4 clustered-base-substitution and 17 small insertion-and-deletion signatures. The substantial size of our dataset, compared with previous analyses^3-15, enabled the discovery of new signatures, the separation of overlapping signatures and the decomposition of signatures into components that may represent associated-but distinct-DNA damage, repair and/or replication mechanisms. By estimating the contribution of each signature to the mutational catalogues of individual cancer genomes, we revealed associations of signatures to exogenous or endogenous exposures, as well as to defective DNA-maintenance processes. However, many signatures are of unknown cause. This analysis provides a systematic perspective on the repertoire of mutational processes that contribute to the development of human cancer.

RNA sequence analysis reveals macroscopic somatic clonal expansion across normal tissues

How somatic mutations accumulate in normal cells is poorly understood. A comprehensive analysis of RNA sequencing data from ~6700 samples across 29 normal tissues revealed multiple somatic variants, demonstrating that macroscopic clones can be found in many normal tissues. We found that sun-exposed skin, esophagus, and lung have a higher mutation burden than other tested tissues, which suggests that environmental factors can promote somatic mosaicism. Mutation burden was associated with both age and tissue-specific cell proliferation rate, highlighting that mutations accumulate over both time and number of cell divisions. Finally, normal tissues were found to harbor mutations in known cancer genes and hotspots. This study provides a broad view of macroscopic clonal expansion in human tissues, thus serving as a foundation for associating clonal expansion with environmental factors, aging, and risk of disease.

Passenger Hotspot Mutations in Cancer

Current statistical models for assessing hotspot significance do not properly account for variation in site-specific mutability, thereby yielding many false-positives. We thus (i) detail a Log-normal-Poisson (LNP) background model that accounts for this variability in a manner consistent with models of mutagenesis; (ii) use it to show that passenger hotspots arise from all common mutational processes; and (iii) apply it to a ∼10,000-patient cohort to nominate driver hotspots with far fewer false-positives compared with conventional methods. Overall, we show that many cancer hotspot mutations recurring at the same genomic site across multiple tumors are actually passenger events, recurring at inherently mutable genomic sites under no positive selection.

Abstract 2727: The premalignant state captured in the landscape of somatic mutations can reveal the cancer cell-of-origin

Despite increasing knowledge of tumorigenesis, the identity of the cancer cell-of-origin, i.e. the normal cell type that acquired the cancer-initiating event, remains largely unknown. Our approach of identifying the cell-of-origin is based on two observations: (1) the chromatin structure is cell-specific; and (2) the density of somatic mutations along the genome is associated with the regional profile of chromatin modifications.

We have previously developed a method that quantifies the ability to predict the mutational distribution along the cancer genome from the profile of epigenetic modifications in different normal cell types. Here we present the largest application of our method using 2,550 whole genomes representing 32 distinct cancer types. To identify the cell-of-origin, we determined the correlation between the observed density of mutations along the genome and the predicted values based on chromatin modifications from 104 different normal tissue types. The normal cell type that showed the strongest correlation with a specific cancer mutational landscape was the candidate cell-of-origin.

We found that in almost all cancer types the cell-of-origin can be characterized solely from DNA sequences. Interestingly, we found that the fallopian tube was the best match for high-grade serous ovarian cancer, providing independent evidence that this is the cancer’s site of origin. For breast cancer we found that the four distinct subtypes best-matched cells from the luminal cell lineage: basal-like breast cancer likely originates from luminal progenitors, whereas all other subtypes from luminal mature cells. This association holds true even when accounting for different alterations in the homologous recombination repair pathway, suggesting that subtypes are more determined by the cell-of-origin than the specific DNA repair defect. In addition, we found that we could identify the cell-of-origin using metastatic samples – a finding that may help in difficult clinical diagnoses. Moreover, we demonstrate that cancer drivers, both germline risk alleles and somatically mutated drivers, reside in active chromatin regions in the respective cell-of-origin.

Taken together, our findings indicate that many of the somatic mutations accumulated while the cells maintained a chromatin structure similar to the cell-of-origin (likely occurring prior to transformation). Therefore, this historical record, captured in the DNA, can be used to identify, the often elusive, cancer cell-of-origin. Our approach can ultimately help better understand the potential of particular normal cell types to transform and initiate cancer, as well as the association of the cell-of-origin with tumor subtypes and sensitivity to treatment.

Citation Format: Kirsten Kubler, Rosa Karlic, Nicholas J. Haradhvala, Kyungsik Ha, Jaegil Kim, Maja Kuzman, Wei Jiao, Sitanshu Gakkhar, Kent W. Mouw, Lior Z. Braunstein, Olivier Elemento, Andrew V. Biankin, Ilse Rooman, Mendy Miller, Christopher D. Nogiec, Edward Curry, Mari Mino-Kenudson, Leif W. Ellisen, Robert Brown, Alexander Gusev, Cristian Tomasetti, Hong-Gee Kim, Hwajin Lee, Kristian Vlahovicek, Charles Sawyers, Katherine A. Hoadley, Edwin Cuppen, Amnon Koren, Peter F. Arndt, David N. Louis, Lincoln Stein, William D. Foulkes, Paz Polak, Gad Getz. The premalignant state captured in the landscape of somatic mutations can reveal the cancer cell-of-origin [abstract]. In: Proceedings of the American Association for Cancer Research Annual Meeting 2019; 2019 Mar 29-Apr 3; Atlanta, GA. Philadelphia (PA): AACR; Cancer Res 2019;79(13 Suppl):Abstract nr 2727.

Abstract 139: Single-cell RNA sequencing reveals compromised immune microenvironment in precursor stages of multiple myeloma

In multiple myeloma (MM), despite well-characterized precursor states such as monoclonal gammopathy of undetermined significance (MGUS) and smoldering multiple myeloma (SMM), there is a lack of sufficient biomarkers to predict disease progression. Most genomic analyses have studied the malignant plasma cells, however, cancers form a complex ecosystem with the immune and stromal microenvironment. To characterize the cellular composition and transcriptional programs of each component of the tumor and microenvironment at different stages of MM progression, we employed single-cell RNA sequencing on 48K plasma and 40.8K immune microenvironmental cells from a cohort of 22 patients with varying stages of disease progression and 9 healthy donors.

Expression profiles of plasma cells revealed clear tumor-specific differences in known oncogenic drivers in MM (MMSET/FGFR3, CCND1 and MAFB) as well as other clonally expressed genes (LAMP5, HIST1H1C, and AREG), distinguishing them from healthy plasma cells. We identified a subset of cycling plasma cells in malignant samples, observing a range of proliferative capacity across disease stages. Furthermore, our approach allowed a unique head-to-head comparison of gene expression changes in normal and malignant plasma cells from the same individual, revealing early alterations in genes related to immune modulation (NKBIA) or controlling transcription and differentiation (EID1). Some alterations were patient-specific, while others, such as MHC I overexpression and CD27 loss, were recurrently observed across subsets of the cohort.

Analysis of the BM microenvironment demonstrated significant infiltration of natural killer cells, non-classical monocytes/macrophages, and T cells, even in the earliest stages of the disease. Further investigation revealed upregulation of MHC II expression at the mRNA level in CD14+ monocytes/macrophages and yet, intriguingly, analysis by CyTOF and immunohistochemistry revealed a shift towards intracellular localization of MHC II in these cells. Co-culture with MM cell lines was sufficient to induce the decrease of extracellular MHC II, providing strong evidence for MM-induced compromised antigen presentation by macrophages, and hinting at a mechanism of immune evasion.

Together, our results provide a comprehensive view at the complex interplay of the immune and malignant cells in different stages of the disease. We demonstrate the immune response beginning in premalignant conditions to be heterogeneous, including compromised antigen presentation as well as alterations in cellular composition and signaling. Consideration of the type of immunological response may prove valuable in determination of progression risk, as well as open up potential strategies for therapy.

Citation Format: Nicholas J. Haradhvala, Oksana Zavidij, Tarek H. Mouhieddine, Romanos Sklavenitis-Pistofidis, Jihye Park, Mairead Reidy, Abdallah Flaifel, Benjamin Ferland, Salomon Manier, Mark Bustoros, Daisy Huynh, Marzia Capelletti, Brianna Berrios, Mahshid Rahmat, Chia-Jen Liu, Meng Xiao He, Esteban Braggio, Rafael Fonseca, Yosef Maruvka, Jennifer Guerriero, Melissa Goldman, Eliezer Van Allen, Steven McCarroll, Jamil Azzi, Gad Getz, Irene M. Ghobrial. Single-cell RNA sequencing reveals compromised immune microenvironment in precursor stages of multiple myeloma [abstract]. In: Proceedings of the American Association for Cancer Research Annual Meeting 2019; 2019 Mar 29-Apr 3; Atlanta, GA. Philadelphia (PA): AACR; Cancer Res 2019;79(13 Suppl):Abstract nr 139.

Quantification of somatic mutation flow across individual cell division events by lineage sequencing

Mutation data reveal the dynamic equilibrium between DNA damage and repair processes in cells and are indispensable to the understanding of age-related diseases, tumor evolution, and the acquisition of drug resistance. However, available genome-wide methods have a limited ability to resolve rare somatic variants and the relationships between these variants. Here, we present lineage sequencing, a new genome sequencing approach that enables somatic event reconstruction by providing quality somatic mutation call sets with resolution as high as the single-cell level in subject lineages. Lineage sequencing entails sampling single cells from a population and sequencing subclonal sample sets derived from these cells such that knowledge of relationships among the cells can be used to jointly call variants across the sample set. This approach integrates data from multiple sequence libraries to support each variant and precisely assigns mutations to lineage segments. We applied lineage sequencing to a human colon cancer cell line with a DNA polymerase epsilon (POLE) proofreading deficiency (HT115) and a human retinal epithelial cell line immortalized by constitutive telomerase expression (RPE1). Cells were cultured under continuous observation to link observed single-cell phenotypes with single-cell mutation data. The high sensitivity, specificity, and resolution of the data provide a unique opportunity for quantitative analysis of variation in mutation rate, spectrum, and correlations among variants. Our data show that mutations arrive with nonuniform probability across sublineages and that DNA lesion dynamics may cause strong correlations between certain mutations.

Structural Alterations Driving Castration-Resistant Prostate Cancer Revealed by Linked-Read Genome Sequencing

Nearly all prostate cancer deaths are from metastatic castration-resistant prostate cancer (mCRPC), but there have been few whole-genome sequencing (WGS) studies of this disease state. We performed linked-read WGS on 23 mCRPC biopsy specimens and analyzed cell-free DNA sequencing data from 86 patients with mCRPC. In addition to frequent rearrangements affecting known prostate cancer genes, we observed complex rearrangements of the AR locus in most cases. Unexpectedly, these rearrangements include highly recurrent tandem duplications involving an upstream enhancer of AR in 70%-87% of cases compared with <2% of primary prostate cancers. A subset of cases displayed AR or MYC enhancer duplication in the context of a genome-wide tandem duplicator phenotype associated with CDK12 inactivation. Our findings highlight the complex genomic structure of mCRPC, nominate alterations that may inform prostate cancer treatment, and suggest that additional recurrent events in the non-coding mCRPC genome remain to be discovered.

Abstract 4443: A post-transcriptional program of chemoresistance regulators in quiescent cancer cells

Quiescent (G0) cells are a clinically relevant fraction in several cancers, which include dormant cancer stem cells, and resist clinical therapy. G0 cells reveal extensive changes in gene expression at the protein and translation levels. We previously identified that the translation mechanism is altered in G0 cancer cells. MicroRNAs, noncoding RNAs that target distinct mRNAs to alter gene expression—were found to associate with an important RNA-binding protein, and enable specialized functions in G0—where they recruit non-canonical translation factors to regulate specific mRNA translation. We find that G0 leukemic cells show similar proteome and translatome to cells isolated post-chemotherapy. These data suggest that specialized post-transcriptional and translational mechanisms in G0 leukemic cells regulate a distinct translatome to mediate chemoresistance.

To understand the role of post-transcriptional and translational regulation in chemoresistance, we compared global RNA, translational and proteome profiling in chemoresistant G0 acute monocytic leukemic (AML) cells. We find that chemotherapy or G0 induction leads to DNA damage responsive ATM and stress signaling, which alter post-transcriptional and translational mechanisms. ATM and stress activated p38 MAPK/MK2 increase AU-rich-element (ARE) bearing pro-inflammatory cytokine and immune gene mRNAs by regulating a key ARE RNA binding protein, and by activating STAT1/interferon pathway to alter canonical translation. AREs are present on 3'UTRs of critical, tightly regulated oncogenes and cytokines to post-transcriptionally control their expression. These changes permit translation of ARE bearing pro-inflammatory cytokine TNFα, and immune and cell-migration modulators that promote survival. Co-inhibiting p38 MAPK and TNFα that promote anti-apoptosis—prior to or alongwith chemotherapy—decreases chemoresistance in AML cell lines, in vivo, and in patient samples, without affecting normal cells. These studies reveal a pro-inflammatory subpopulation in AML that mediates chemoresistance, enabled by DNA damage- and stress-regulated post-transcriptional and translational mechanisms that are mediated by AU-rich-elements and a critical ARE RNA binding protein. Disrupting ARE regulation reduces TNFα and chemoresistance, revealing AREs and an important ARE RNA binding protein as key regulators of inflammation-mediated chemoresistance. These studies reveal the significance of post-transcriptional regulation of inflammation/immune gene-mediated chemoresistance.

Correspondence: vasudevan.shobha@mgh.harvard.edu

Citation Format: Sooncheol Lee, Samuel S. Truesdell, Syed I. Bukhari, Myriam Boukhali, Dongjun Lee, Maria A. Mazzola, Radhika Raheja, Adam Langenbucher, Nicholas J. Haradhvala, Michael Lawrence, Roopali Gandhi, David A. Sweetser, Wilhelm Haas, Shobha Vasudevan. A post-transcriptional program of chemoresistance regulators in quiescent cancer cells [abstract]. In: Proceedings of the American Association for Cancer Research Annual Meeting 2018; 2018 Apr 14-18; Chicago, IL. Philadelphia (PA): AACR; Cancer Res 2018;78(13 Suppl):Abstract nr 4443.

Distinct mutational signatures characterize concurrent loss of polymerase proofreading and mismatch repair

Fidelity of DNA replication is maintained using polymerase proofreading and the mismatch repair pathway. Tumors with loss of function of either mechanism have elevated mutation rates with characteristic mutational signatures. Here we report that tumors with concurrent loss of both polymerase proofreading and mismatch repair function have mutational patterns that are not a simple sum of the signatures of the individual alterations, but correspond to distinct, previously unexplained signatures: COSMIC database signatures 14 and 20. We then demonstrate that in all five cases in which the chronological order of events could be determined, polymerase epsilon proofreading alterations precede the defect in mismatch repair. Overall, we illustrate that multiple distinct mutational signatures can result from different combinations of a smaller number of mutational processes (of either damage or repair), which can influence the interpretation and discovery of mutational signatures.

Polymerase proofreading and the mismatch repair pathway maintain the fidelity of DNA replication. Here the authors show that tumors with concurrent loss of function of both pathways lead to mutation signatures that are not simply a sum of the signatures found in tumors involving single alteration.

A mutational signature reveals alterations underlying deficient homologous recombination repair in breast cancer

Biallelic inactivation of BRCA1 or BRCA2 is associated with a pattern of genome-wide mutations known as signature 3. By analyzing ∼1,000 breast cancer samples, we confirmed this association and established that germline nonsense and frameshift variants in PALB2, but not in ATM or CHEK2, can also give rise to the same signature. We were able to accurately classify missense BRCA1 or BRCA2 variants known to impair homologous recombination (HR) on the basis of this signature. Finally, we show that epigenetic silencing of RAD51C and BRCA1 by promoter methylation is strongly associated with signature 3 and, in our data set, was highly enriched in basal-like breast cancers in young individuals of African descent.

Abstract LB-280: The landscape of somatic microsatellite indels across cancer: detection and identification of driver events

Microsatellites (MSs) are tracts of variable-length repeats of short DNA motifs that are abundant in the human genome and exhibit high rates of mutations in the form of insertions or deletions of the repeated motif (MS indels). Despite their prevalence, the contribution of somatic MS indels to cancer is largely unexplored due to difficulties in detecting them and assessing their significance. Here, we present a comprehensive analysis of MS indels across 20 tumor types. We characterize the overall MS indel landscape and detect genes with candidate driver MS indel events. We present two novel tools: MSMuTect for accurate detection of somatic MS indels and MSMutSig for identifying candidate cancer genes containing events at higher frequency than expected by chance. We observe high variability of the frequency of MS indels across tumors and demonstrate that the number and pattern of MS indels can accurately distinguish microsatellite stable (MSS) tumors from tumors with microsatellite instability (MSI). Applying MSMutSig across 6,788 tumors from 20 different tumor types identified 7 genes with significant MS indel hotspots: ACVR2A, RNF43, DOCK3, MSH3, ESRP1, PRDM2 and JAK1. In the four genes that have been previously implicated in cancer (ACVR2A, RNF43, JAK1 and MSH3), we identified previously unreported MS indels events. Three of the genes with significant loci - DOCK3, PRDM2 and ESRP1- had not been previously listed as cancer genes. MS indels in DOCK3, a negative regulator of the WNT pathway, were mutually exclusive with mutations in CTNNB1. MS indels in ESRP1, an RNA processing gene, correlated with alternative splicing of FGFR2, an event associated with the epithelial-to-mesenchymal transition. Overall, our comprehensive analysis of somatic MS indels across cancer highlights their importance, particularly in

MSI tumors, significantly contributes to the ongoing global efforts to detect cancer genes, and may improve classification of patients into clinically-relevant subgroups.

Citation Format: Yosef E. Maruvka, Kent W. Mouw, Rosa Karlic, Rasanna Parasuraman, Atanas Kamburov, Paz Polak, Nicholas J. Haradhvala, Julian M. Hess, Esther Rheinbay, Yehuda Brody, Lior Z. Braunstein, Alan D’Andrea, Michael S. Lawrence, Adam Bass, Andre Bernards, Franziska Michor, Gad Getz. The landscape of somatic microsatellite indels across cancer: detection and identification of driver events [abstract]. In: Proceedings of the American Association for Cancer Research Annual Meeting 2017; 2017 Apr 1-5; Washington, DC. Philadelphia (PA): AACR; Cancer Res 2017;77(13 Suppl):Abstract nr LB-280. doi:10.1158/1538-7445.AM2017-LB-280

Mutational Strand Asymmetries in Cancer Genomes Reveal Mechanisms of DNA Damage and Repair

Mutational processes constantly shape the somatic genome, leading to immunity, aging, cancer, and other diseases. When cancer is the outcome, we are afforded a glimpse into these processes by the clonal expansion of the malignant cell. Here, we characterize a less explored layer of the mutational landscape of cancer: mutational asymmetries between the two DNA strands. Analyzing whole-genome sequences of 590 tumors from 14 different cancer types, we reveal widespread asymmetries across mutagenic processes, with transcriptional ("T-class") asymmetry dominating UV-, smoking-, and liver-cancer-associated mutations and replicative ("R-class") asymmetry dominating POLE-, APOBEC-, and MSI-associated mutations. We report a striking phenomenon of transcription-coupled damage (TCD) on the non-transcribed DNA strand and provide evidence that APOBEC mutagenesis occurs on the lagging-strand template during DNA replication. As more genomes are sequenced, studying and classifying their asymmetries will illuminate the underlying biological mechanisms of DNA damage and repair.