Kannan, Nagarajan

Abstract Cancer can recur when a subset of tumor cells, denoted here as persister cells, are able to survive therapy and re-enter the cell cycle. The precise mechanisms that confer the persister state and whether it is characteristic of a subgroup of cells or arises from multiple cellular lineages remain poorly understood. We hypothesize that an epigenetic signature underlies the drug-tolerant persister state, characterized by transcriptional and chromatin accessibility changes that promote survival of residual cancer following chemotherapy. To identify clinically relevant features of persister cells in untreated tumors and residual disease, we performed single-cell multiomic profiling (snRNA+snATAC) on a cohort of non-malignant fallopian tube, treatment-naïve, and neoadjuvant chemotherapy (NACT)-treated high-grade serous ovarian cancer (HGSOC) samples. We identified differences in gene expression and open chromatin between naïve and residual patient tumors following chemotherapy. Although only a small proportion of the differentially expressed genes enriched in residual HGSOC overlapped with established gene sets for chemo-response and patient prognosis, the epigenomic analysis revealed activity of several DNA-binding factors that are both enriched upon chemotherapy and also high in resistant tumors prior to treatment. From this analysis, we identified an epigenetic signature that precedes expression and defines the persister state. This epigenetic signature also correlated with chemotherapy sensitivity and resistance using patient-derived xenograft models of HGSOC. Gene regulatory networks driven by the persister signature are involved in the activation of oncogenic pathways, including changes to the cell cycle promoting quiescence and stress response. Further study of the persister cells identified by this epigenetic signature may increase understanding of the mechanisms underlying persister cell survival and reveal new vulnerabilities that could be exploited to delay or prevent cancer recurrence.

ABSTRACT In our efforts to enhance sensitivity to PARP inhibitors, we identified clofarabine (CLF) as a potential therapy for drug-resistant ovarian cancer and nuclear trafficking of Cathepsin L (CTSL) as a treatment- responsive biomarker. Using PARP inhibitor-sensitive and -resistant OC cell lines, ex vivo cultures of patient-derived ovarian ascites (OVA), primary ovarian tumors, and xenografts (PDX), we found that CLF monotherapy induces nuclear CTSL (nCTSL) in CLF-responsive cells (CLF-r) and sensitizes them to PARP inhibitors olaparib and rucaparib. In CLF non-responsive cells (CLF-nr), a combination of CLF with olaparib is necessary for nCTSL trafficking and synergy. CLF+olaparib synergy was observed in 47% of CLF-r and 24% of CLF-nr OVA samples. Drug-induced nCTSL is crucial for DNA damage response, including cell cycle arrest and apoptosis. Knockdown of CTSL in both CLF-r and CLF-nr cells conferred resistance to the CLF+olaparib combination, emphasizing nCTSL’s role in the DNA damage response pathway (DDR). Mechanistically, CLF facilitates CTSL nuclear import via KPNB1 in CLF-r cells. In CLF-nr cells, both olaparib and CLF are needed to facilitate CTSL nuclear import. Additionally, CLF downregulates the nuclear export protein CRM1 (XPO1) in both cohorts. Interestingly, CLF does not downregulate CRM1 in a subset of OVAs (29%), and they were classified as CLF-resistant (CLF- Res). In these samples, inhibiting CRM1 with KPT8602 restored synergy between CLF and PARP inhibitors. In vivo, CLF-r and CLF-nr PDX models exhibited enhanced DDR, reduced tumor burden, and prolonged survival with the CLF+olaparib combination. These findings suggest the CLF+olaparib combination is a promising therapeutic strategy for drug-resistant OC by inducing DDR through CTSL nuclear localization.

Multiple myeloma (MM) is an incurable plasma cell (PC) malignancy. Despite therapeutic advances that have vastly improved survival, all patients are expected to relapse. Current standard of care includes autologous stem cell transplant (ASCT) with high-dose melphalan (HDM) administered as a conditioning agent. Cytotoxic chemotherapies such as melphalan have been shown to drive therapy-induced senescence (TIS) in several tissues. We hypothesize that HDM activates TIS pathways in surviving MM cells and that this may correlate with patient outcomes post-ASCT. We developed an in vitro HDM model using 5TGM1 mouse MM cells treated with vehicle (Veh) or HDM (10µM, IC90) for 6h, followed by co-culture with primary mouse bone marrow stromal cells. Image analysis was performed to quantify proliferation by cumulative population doublings (CPDs) over 10 days, after which 5TGM1 cells were evaluated for cell size, telomere-associated DNA damage foci (TAFs), sensitivity to senolysis, and senescence gene expression (scRNA-seq and RT-qPCR). The effect of HDM on CPDs and cell size was also evaluated in human MM cell lines RPMI-8226 and H929. scRNA-seq was used to analyze CD138+ PCs from MM patients at various timepoints post-ASCT. Longitudinal MM patient iliac crest biopsies (diagnosis, post-ASCT; N=20) were stained for CD138+ PCs; PC burden was quantified by AI-assisted histology. Lastly, we established an in vivo model of induction therapy (bortezomib) followed by HDM/SCT in Vk*MYC mouse MM engrafted C57BL/6 mice. Tumor burden was tracked by serum protein electrophoresis for M-spike using submandibular bleeds. Bone marrow was then flushed from long bones and analyzed by flow cytometry and immunofluorescence. HDM-5TGM1 showed decreased CPDs (HDM 0.32±0.31, Veh 6.72±0.69, p<0.0001, N=3), increased size (median pixel size: HDM 413.2±99.98, Veh 207±8.88, p<0.05, N=3), and increased TAF+ cells (Percent TAF+ cells: HDM 34.77%±10.74%, Veh 11.41%±1.58%, p<0.05, N=3) consistent with known TIS features. Preliminary assessment (N=1) of human MM cell lines showed similar growth arrest (RPMI-8226 CPDs: HDM 0.058±0.101, Veh 4.185±0.358; H929 CPDs: HDM 0±0, Veh 4.61±0.243) and size increase (RPMI-8226 mean pixel size: HDM 804.1±6.048, Veh 526.9±2.115; H929 mean pixel size: HDM 371.2±26.95, Veh 130.3±41.43), suggesting a translatable mechanism. HDM-5TGM1 cells were more sensitive to senolysis by dasatinib+quercetin (viability: HDM 54.9%±0.09%, Veh 76.0%±0.05%, N=2). scRNA-seq analysis showed that HDM-5TGM1 were enriched for gene sets related to senescence pathways, including the published set SenMayo (FDR<0.25). HDM-5TGM1 cells also exhibited myeloid gene expression consistent with previous publications on MM dormancy. RT-qPCR confirmed significant increases in HDM-5TGM1 (p<0.001, N=3) in senescence markers (Cdkn1a, Cdkn1c, Glb1), senescent cell anti-apoptosis pathway genes (SCAPs, Bcl2l1), senescence associated secretory phenotype (SASP) genes (Ccl5, Icam1, Mmp13), and myeloid markers (Axl, Fcer1g, Mpeg1). scRNA-seq of patient CD138+ PCs confirmed a myeloid-like cluster that gained prominence in relapsed patients with resistant disease. CD138+ PCs were detected in patient bone biopsies post-ASCT. Interestingly, patients who relapsed (≤3 years, N=9) had a greater percent reduction in PC burden from diagnosis to post-ASCT when compared to patients that did not relapse (N=11) in that time (p=0.02). Further, CD138+ PC burden post-ASCT positively correlated with time in remission in relapsed patients (R=0.697, p=0.03). Bortezomib reduced serum M-spike in Vk*MYC-engrafted mice (p=0.0003); M-spike was further reduced by HDM/SCT compared to Veh-treated mice (p=0.0007). Flow cytometry confirmed a decrease in Vk*MYC cells after HDM/SCT (p=0.0043), and the percentage of Vk*MYC cells with greater size and complexity was increased compared to Veh-treated mice (p=0.0043), consistent with our in vitro data. Vk*MYC cells from HDM-treated mice also trended for more TAF+ cells (>3 TAFs) compared to Veh-treated mice (p=0.0714). Our findings suggest that HDM can induce a TIS and myeloid signature in surviving MM cells. Greater MM burden post-ASCT appears to be associated with longer durable response, again suggesting that these cells may be in growth arrest consistent with TIS. Finally, senolytics may be a novel approach to eliminate surviving MM cells and prevent relapse.

Abstract Multiple myeloma (MM) is an incurable plasma cell (PC) cancer in which all patients are expected to relapse. Cytotoxic agents such as high-dose melphalan (HDM), which is used as a myeloablative agent prior to autologous stem cell transplant (ASCT) in MM, drive therapy-induced senescence (TIS) in non-cancerous tissues. We hypothesized that HDM similarly activates TIS pathways in surviving MM cells and this response may correlate with patient outcomes post-ASCT. To test this, we developed an in vitro model in which 5TGM1 mouse MM cells are cultured with 10uM melphalan (HDM) or vehicle (Veh) for 6 hours, followed by co-culture in normal media with primary mouse bone marrow stromal cells. Cultures were maintained for 10 days and imaged to quantify cumulative population doublings (CPDs). 5TGM1 were then isolated by FACS to measure telomere-associated DNA damage foci (TAFs), sensitivity to senolytic drugs, and senescence gene expression (RT-qPCR and scRNA-seq). To translate these studies, scRNA-seq was used to analyze CD138+ PCs from MM patient bone marrow post-ASCT. Additionally, longitudinal bone biopsies from MM patients at diagnosis and post-ASCT were stained for CD138+ PCs; PC burden was quantified using artificial intelligence assisted histology (HALO). HDM-5TGM1 exhibited stable growth arrest (decreased CPDs) and sustained DNA damage (increased TAFs) compared to Veh-5TGM1. HDM-5TGM1 were also significantly more sensitive to senolytic ablation (Dasatinib+Quercetin or Fisetin). HDM-5TGM1 showed increased senescence (Cdkn1a, Cdkn1c, Glb1), anti-apoptosis (Bcl2l1), and senescence associated secretory phenotype (Ccl5, Icam1, Mmp13) genes, as well as myeloid markers found to be increased in MM dormancy (Axl, Fcer1g, Mpeg1). Analysis by scRNA-seq revealed that HDM-5TGM1 were G0/G1 arrested, and gene set enrichment analysis (GSEA) confirmed upregulation of senescence pathways with downregulation of DNA replication and cell cycle pathways. Consistent with RT-qPCR results, these cells exhibited myeloid gene expression. Analysis of patient CD138+ PCs post-ASCT by scRNA-seq also revealed a myeloid-like cluster. Of interest, of the 491 genes upregulated in this cluster, 268 were also upregulated in HDM-5TGM1. CD138+ PCs could also be detected in patient bone biopsies post-ASCT, although PC burden was significantly reduced post-ASCT. Surprisingly, the percent reduction in PC burden post-ASCT was greater in relapsed (≤3 years) vs non-relapsed patients. Further, in relapsed patients, PC burden post-ASCT positively correlated with time in remission. Altogether, these findings demonstrate that HDM can induce senescence features and a myeloid-dormancy signature in surviving MM cells. This phenotype is associated with longer durable response, suggesting that activation of TIS may be a strategy to extend PFS. Further, senolytic therapy may be a novel approach to eliminate dormant MM cells and prevent relapse. Citation Format: Angelo Jose Guilatco, Marta Diaz-del-Castillo, Gabriel Alvares Borges, Neal I. Sannuli, Megan L. Ritting, Syed Mohammed Aalam, Nagarajan Kannan, Tamar Tchkonia, James L. Kirkland, Yi Lin, Taxiarchis Kourelis, Matthew T. Drake, Thomas Levin Andersen, Megan Weivoda. Melphalan-treated multiple myeloma cells exhibit a senescent-like dormancy phenotype [abstract]. In: Proceedings of the American Association for Cancer Research Annual Meeting 2024; Part 1 (Regular Abstracts); 2024 Apr 5-10; San Diego, CA. Philadelphia (PA): AACR; Cancer Res 2024;84(6_Suppl):Abstract nr 2958.

Abstract Objective: Prolonged use of oral contraceptives (OC) has been associated with a reduced risk of ovarian cancer, yet the underlying mechanism remains elusive. Long-lived self-renewing fallopian tube epithelial stem cells (SCs) capable of cilio-secretory differentiation are considered a potential cell-of-origin for ovarian cancers. Using innovative stem cell methodologies, our objective is to gain insights into the mechanism by which OC exert their protective action. Methods: This research is supported by Mayo Clinic’s clinically annotated regenerative fallopian tube biobank. RNA was isolated from the SC-enriched cell fraction and differentiated cilio-secretory organoids. Absolute mRNA copies of OC targets estrogen receptor-alpha (ESR1; ER) and progesterone receptor (PGR; PR) were determined using optimum RT-ddPCR assay. OC effects on SC-enriched growth in 2D was measured using live-cell imaging. Our primary tubal cells and organoids as well as public single-cell transcriptomics datasets were analyzed for ESR1/PGR expression across tubal cell states. Results: The panel of OCs tested exhibited no growth response in patient-derived cultures enriched with stem cells, and there was an absence of ESR1/PGR expression. Clonal organoid differentiation exhibited high ESR1 levels and lack of PGR, while growth factor cocktail-dependent PGR induction was achieved. Single-cell transcriptomics identified molecular signatures associated with ER/PR-positive OC sensor cells. Conclusions: Based on our findings, we propose a model wherein ER/PR-negative SC differentiate into ER/PR-positive OC-sensors and control SC development through paracrine signals. Our in-vitro conditions to differentiate SC into OC-sensors provide an opportunity to study complex paracrine signaling and explore targets for prevention of ovarian cancers. Citation Format: Megan Ritting, Syed Mohammed Musheer Aalam, Manny Bettencourt, Kinnera Jalagam, Scott Kaufmann, John Weroha, Mark Sherman, Jamie Bakkum-Gamez, Nagarajan Kannan. Lack of ER/PR expression and unresponsiveness to oral contraceptives in fallopian tube stem/progenitor cell states: reevaluating the mechanism of cancer prevention [abstract]. In: Proceedings of the AACR Special Conference on Ovarian Cancer; 2023 Oct 5-7; Boston, Massachusetts. Philadelphia (PA): AACR; Cancer Res 2024;84(5 Suppl_2):Abstract nr B046.

Abstract Ovarian cancer is one of the most common gynecologic malignancies, with high-grade serous carcinoma (HGSC) as the most frequent and deadliest type of ovarian cancer (OC). With the recent advances in cancer epigenomic profiling, it has become clear that intra-tumor heterogeneity driven by genetic and epigenetic factors may be the basis of resistance, where current treatments could efficiently target sub-populations of cells while a sub-group may be resistant, or prone to develop resistance. With epigenomic profiling and single cell sequencing technology, we aim to characterize cell signatures that translate as prognostic markers and novel drug targets. We have mapped the epigenomic landscape of HGSC in 12 patient-derived xenografts (PDX), and their metastatic counterparts (circulating tumor cells and ascites present in the PDX models) using an assay for transposase-accessible chromatin (or ATAC-seq) and compared it with normal fallopian tube organoids to identify tumor and metastatic-specific regulatory elements. We also performed RNA-seq to validate potential transcriptional dependencies associated with resistance to chemotherapy. In order to identify potentially chemotherapy-resistant subpopulations and further validate the activity of some of the transcription factors that can be driving resistance, we performed single cell multiomic profiling using the 10x Multiome sequencing platform (scRNA+scATAC from the same nuclei) in HGSC patient samples and PDX. Finally, we performed immunohistochemistry of patient tissue slides for target markers and correlated it with overall survival in OC patients. As previously suggested, a subset of cancer clonal cells can drive treatment resistance and metastasis and result in poor patient outcomes. Data derived from PDX ATAC-seq showed a cancer stem cell signature driven by SOX2 in metastatic samples compared to the paired PDX primary tumor, suggesting a role for SOX2 in disease progression. Interestingly, SOX2 expression in epithelial cells in the Multiome analysis is more patient-specific than other stem cell markers, like CD133. We also observed that SOX2 expression negatively correlated with prognosis. These data, taken together, point out the promising role of SOX2 as a prognostic and targetable marker for HGSC. Citation Format: Leticia Assad Maia Sandoval, Wazim Mohammed Ismail, Amelia Mazzone, Amik Munankarmy, Jagneet Kaur, Tsering Tashi, Xiaonan Hou, Nagarajan Kannan, Ann Oberg, Scott H. Kaufmann, Saravut John Weroha, Alexandre Gaspar-Maia. Single cell multiomic profiling of high-grade serous carcinoma. [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 4724.

Abstract A problematic feature of many human cancers is a lack of understanding of mechanisms controlling organ-specific patterns of metastasis, despite recent progress in identifying many mutations and transcriptional programs shown to confer this potential. To address this gap, we developed a methodology that enables different aspects of the metastatic process to be comprehensively characterized at a clonal resolution. Our approach exploits the application of a computational pipeline to analyze and visualize clonal data obtained from transplant experiments in which a cellular DNA barcoding strategy is used to distinguish the separate clonal contributions of two or more competing cell populations. To illustrate the power of this methodology, we demonstrate its ability to discriminate the metastatic behavior in immunodeficient mice of a well-established human metastatic cancer cell line and its co-transplanted LRRC15 knockdown derivative. We also show how the use of machine learning to quantify clone-initiating cell (CIC) numbers and their subsequent metastatic progeny generated in different sites can reveal previously unknown relationships between different cellular genotypes and their initial sites of implantation with their subsequent respective dissemination patterns. These findings underscore the potential of such combined genomic and computational methodologies to identify new clonally-relevant drivers of site-specific patterns of metastasis.

Abstract Purpose Breast cancer risk is elevated in pathogenic germline BRCA 1/2 mutation carriers due to compromised DNA quality control and elevated mutagenic load. We hypothesized that if immunosurveillance promotes tumor suppression, then normal breast lobules from BRCA carriers may demonstrate higher immune cell densities. Methods We assessed immune cell composition in normal breast lobules from age-matched women with progressively increased breast cancer risk, including 1) 19 women who donated normal breast tissue to the Komen Tissue Bank (KTB) at Indiana University Simon Cancer Center, 2) 15 women with biopsy-identified benign breast disease (BBD), and 3) 19 prophylactic mastectomies from women with germline mutations in BRCA1/2 genes. We performed immunohistochemical stains and analysis to quantitate immune cell densities from digital images in up to 10 representative lobules per sample. Median cell counts per mm2 were compared between groups using Wilcoxon rank-sum tests. Results Compared to BBD tissues, normal breast lobules from BRCA carriers had significantly higher densities of CD8 + cells (median 354.4 vs 200.0 cells/mm2, p = 0.01). Compared to KTB normal donor breast tissues, normal breast lobules from BRCA carriers also had significantly higher densities of CD4+ (median 116.3 vs 17.7 cells/mm2), CD8+ (354.4 vs 150.9 cells/mm2), CD11c+ (3.5% vs 0.4% pixels positive), and CD68+ (237.5 vs 57.8 cells/mm2) immune cells (each p < 0.001). Conclusions These preliminary analyses show that normal/benign breast lobules of BRCA mutation carriers contain increased immune cells compared with BBD and normal donated tissues. Future studies will be helpful to define relationships of mutations and immune responses.

Abstract Polycystic ovarian syndrome (PCOS) is a complex health condition associated with metabolic disturbances and infertility. Recent data suggest that the prevalence of PCOS is increasing among women globally, although the etiology of these trends is undefined. Consequently, preclinical models that better reflect the biology of PCOS are urgently needed to facilitate research that can lead to the discovery of prevention strategies or improved management. The existing animal models have several limitations as they do not reflect all the PCOS features metabolically and/or phenotypically. Therefore, there is no clear consensus on the use of appropriate animal model and selection of the most appropriate PCOS-inducing agent. To that end, we have established a Swiss albino mouse model of PCOS based on 3 weeks of daily treatment with letrozole (50 μg/day; intraperitoneal) and dehydroepiandrosterone (DHEA, 6 mg/100 g body weight; subcutaneous) in 5-week-old female mice fed on normal or high-fat diet (HFD). Mice were regularly assessed for body weight, blood glucose, and estrous cycle. Three weeks after drug administration, mice were sacrificed and assessed for blood-based metabolic parameters as well as ovarian function. Our results indicate that DHEA combined with HFD produces changes mimicking those of clinical PCOS, including elevated serum testosterone and luteinizing hormone, dyslipidemia, poor ovarian microenvironment, and development of multiple ovarian cysts, recapitulating cardinal features of PCOS. In comparison, normal diet and/or letrozole produced fewer features of PCOS. The data from the experimental models presented here can improve our understanding of PCOS, a growing concern in women’s health.