Indulkar, Shreya

Endocrine therapy remains the cornerstone of treatment for estrogen receptor alpha (ERα)-positive breast cancer, yet the relationship between individual genetic background and molecular response to ERα-targeted therapy remains incompletely understood. Using a well-characterized panel of lymphoblastoid cell lines (LCLs), we performed genome-wide pharmacogenomic expression quantitative trait locus (PGx-eQTL) analysis to identify estradiol (E2)- and tamoxifen (TAM)-induced SNP-gene pairs. PGx-eQTL signals were integrated with previously published breast cancer genome-wide association study datasets to examine their association with clinically relevant breast cancer phenotypes. We identified two ER-mediated PGx-eQTL SNP-gene pairs associated with breast cancer prognosis post-treatment with E2 or TAM. Notable loci included E2-regulated MRPL15 and TAM-regulated SYCP3, which have effects in a genotype-dependent manner, with genotype-dependent survival outcomes, from worse to better relapse-free survival. Similar endocrine-therapy effects on patient survival and breast cancer risk were observed in SIK2, post TAM-treatment, and in LSM4 with E2-treatment. Moreover, qRT-PCR validation in an independent LCL panel confirmed the genotype-specificity of these signals. Overall, we identified ER-mediated PGx-eQTL SNP-gene pairs which represent potential pharmacogenomic tools for identifying patients likely to benefit from ERα-targeted endocrine therapy, offering a foundation for more genotype-informed individualized treatment decisions in breast cancer.

Abstract Poly (ADP-ribose) polymerase (PARP) inhibitors (PARPi) are a mainstay therapy for homologous recombination (HR)-deficient cancers; however, resistance remains a major clinical challenge. Previously, through a genome wide CRISPR screen, we identified ZNF251 haploinsufficiency as a novel driver of PARPi resistance. In BRCA1 -mutant ( BRCA1 mut ) cells, ZNF251 deficiency led to HR hyperactivation, conferring PARPi resistance that could be reversed by RAD51 inhibition. In this study, we further show that ZNF251 deficiency induces replication stress and a pseudo-senescence state, in which cells exhibit molecular and phenotypic markers of senescence while retaining proliferative capacity. Because senomorphic and senolytic therapies can target senescent cells, we tested whether these approaches could overcome PARPi resistance in ZNF251 -deficient breast cancer cells. Critically, targeting this senescence-like state with either senomorphic agents, such as cytokine inhibitors, or senolytic agents, such as BCL-2 and BCL-XL inhibitors, overcame PARPi resistance ex vivo and in vivo . Importantly, pseudo-senescence was also observed in other PARPi-resistant contexts driven by HR hyperactivation, including 53BP1 - and Shieldin -mutant cells, suggesting that it may represent a broader mechanism underlying PARPi resistance in breast cancer. Furthermore, in two olaparib-resistant, BRCA -mutant triple-negative breast cancer organoid models, treatment with DT2216, a BCL-XL–targeting PROTAC, sensitized both models to olaparib. Together, our work defines a novel pathway linking HR hyperactivation, replication stress, and pseudo-senescence, and positions both senomorphic and senolytic therapies as promising strategies to overcome PARPi resistance in BRCA1 mut breast cancer. Highlights ZNF251 deficiency drives PARP inhibitor resistance through HR hyperactivation, replication stress, and pseudo-senescence in BRCA1 -mutant breast cancer. Senomorphic and senolytic therapies overcome PARP inhibitor resistance in vitro , in vivo , and in patient-derived organoid models. Pseudo-senescence represents a shared vulnerability of HR-hyperactivated PARPi-resistant cancers and can be therapeutically targeted.

Abstract Background: Approximately 70% of breast cancers are estrogen receptor-positive (ER+), and endocrine therapies, including selective estrogen receptor modulators (SERMs) and aromatase inhibitors (AIs), are standard-of-care. However, interpatient variability in treatment response and relapse remains substantial, suggesting underlying germline pharmacogenomic factors that modify ER signaling. Here, we mapped ER pharmacogenomic eQTLs to identify germline variants that modulate transcriptional responses to estradiol (E2) and tamoxifen (Tam). Methods: We treated 30 genotyped lymphoblastoid cell lines (LCLs) with vehicle, E2, Tam, or E2+Tam and performed RNA-seq to quantify ligand-induced expression fold changes. Using MatrixEQTL, we identified SNP-gene PGx-eQTLs (cis, ±200 kb), requiring evidence of ER transcriptional specificity via (i) reversal by dual treatment, and (ii) overlap with ER ChIP-seq binding annotations. Clinical relevance was assessed by linkage to breast cancer GWAS datasets, including NSABP, MA.27, and M3. Chromatin regulatory states were annotated (ChromHMM), and transcription factor motif enrichment was performed from SNP-centered ±200 bp regions. Results: Distinct PGx-eQTL architectures were identified for E2 and Tam, with partial overlap in ligand-responsive loci. ER PGx-eQTL variants were predominantly enriched in enhancer chromatin and co-localized with canonical ER-binding motifs. Integrative analysis with GWAS of breast cancer outcomes highlighted multiple drug-response-associated loci, including rs311392-MRPL15 and rs7953325-SYCP3. High MRPL15 expression was linked to shorter relapse-free survival, whereas high SYCP3 expression correlated with improved outcomes. Additional loci derived from a GWAS catalog for cancer and hormone phenotypes included rs2850247-SIK2 and rs11673604-LSM4, where higher SIK2 expression was associated with improved relapse-free survival and higher LSM4 expression with poorer outcomes. Conclusions: We identified ligand-dependent ER PGx-eQTLs that mechanistically connect germline regulatory variants to transcriptional response heterogeneity and clinical endocrine therapy outcomes. These SNP-gene relationships represent potential biomarkers for precision endocrine therapy selection and may reveal endocrine-dependent therapeutic vulnerabilities. Citation Format: Martin Meng, Huanyao Gao, August John, Arnab Ghosh, Shreya Indulkar, Liewei Wang, . Estrogen receptor-mediated pharmacogenomic eQTLs as predictors of endocrine therapy response in breast cancer [abstract]. In: Proceedings of the American Association for Cancer Research Annual Meeting 2026; Part 1 (Regular Abstracts); 2026 Apr 17-22; San Diego, CA. Philadelphia (PA): AACR; Cancer Res 2026;86(7 Suppl):Abstract nr 3610.

Abstract Background: Second generation androgen receptor signaling inhibitors (ARSIs), such as abiraterone, have demonstrated improved survival in metastatic castration-resistant prostate cancer (mCRPC) patients, however, acquired resistance to these agents often arises. Novel and alternative therapies are needed to combat the ARSI resistant tumor. In the last two decades patient-derived xenograft (PDX) and PDX organoids have become pivotal models in anticancer drug development and biomarker discovery due to their ability to faithfully recapitulate the complex and heterogeneous tumor biology. Methods: We analyzed RNAseq data from mCRPC patients and corresponding PDX models (PROMOTE study, NCT #01953640) and employed L1000 Chem Perturbation down/up signatures database to find candidate drugs to overcome abiraterone resistance by reverse signature search. We then performed proliferation assay using 49 compounds in three PROMOTE-PDX organoid models (2 abiraterone-resistant and 1 abiraterone-sensitive) using our semi high-throughput drug screening platform. Results: Overall, we identified a total of 32 compounds with IC50<10µM in at least one of the models, among which, 5 agents are active in all three models and 5 in at least two of the models. (Table 1). Within active compounds we find selective kinase inhibitors (BRAF, JAK2, AKT, FLT3) and histone deacetylase inhibitors (belinostat and pracinostat). Particularly interesting are two drugs that are only active in the abiraterone resistant setting but not the sensitive PDX. Conclusions: In summary, we leveraged RNAseq data from mCRPC patients and PDX models and performed a high-throughput drug screen in PDX-derived organoid models to identify potential drugs that could overcome abiraterone-resistance in lethal advanced prostate cancer. These data highlight potential drug combinations or subsequent therapeutic strategies for these patients. Citation Format: Irene Marín-Goñi, Huanyao Gao, Joachim L. Petit, Shreya Indulkar, Adam M. Kase, Cassandra N. Moore, Krishna R. Kalari, Michael T. Barrett, Richard Weinshilboum, Winston Tan, John A. Copland, Liewei Wang. Overcoming resistance in mCRPC: drug discovery using PDX models and high-throughput drug screening [abstract]. In: Proceedings of the American Association for Cancer Research Annual Meeting 2025; Part 1 (Regular Abstracts); 2025 Apr 25-30; Chicago, IL. Philadelphia (PA): AACR; Cancer Res 2025;85(8_Suppl_1):Abstract nr 5658.