Yu, Jia

Abstract Background: Breast cancer is the most frequently diagnosed cancer and the leading cause of cancer-related death among women worldwide. Triple negative breast cancer (TNBC) represents 15-20% of all breast cancers and is an aggressive subtype. Chemotherapy-based treatments remain the standard of care for TNBC. Unfortunately, chemotherapy resistance is common, and for these patients, outcomes are poor and alternative treatment strategies remain an unmet need. circRNAs are a newly identified class of noncoding RNA molecules with covalently closed circular structures. An increasing number of recent studies including ours have indicated that circRNAs play crucial roles in regulating tumor development and chemoresistance. However, the role of circRNAs in the process of chemotherapy resistance and TNBC progression is not clear. Materials and Methods: As a first step towards identifying circRNAs that participate in the development of chemoresistance in TNBC cells, and to determine if targeting such circRNAs is a novel and efficacious therapeutic strategy, doxorubicin-resistant (Doxo-R), paclitaxel-resistant (PTX-R), and double-resistant (DP-R) cell lines were generated from MDA-MB-231. Human circRNA microarrays were utilized to profile the expression of approximately 14,000 known circRNAs in normal breast tissue, matched patient-derived xenografts (PDX) generated prior to and following neoadjuvant chemotherapy (NAC), and TNBC chemosensitive and chemoresistant cell lines. Top hits were validated using RT-PCR. Results: circRNA microarray profiling identified 429 and 310 transcripts differentially expressed in doxorubicin and paclitaxel resistant cells, respectively, compared to parental chemosensitive cell lines (|FC| ≥ 1.5; p value < 0.05). In comparison to pre-NAC derived xenografts, 1,396 circRNAs were dysregulated among post-NAC PDX models. Further, three circRNAs (hsa_circ_001388, hsa_circ_104652, and hsa_circ_061260) were upregulated in Doxo-R, PTX-R, and post-NAC PDX samples compared to their respective controls. Among these three circRNAs, hsa_circ_001388 (also known as circNSD2) was the only transcript also predicted to be translated into a novel and uncharacterized protein given the presence of a high confidence translation initiation site and IRES sequence. Ongoing studies are aimed at determining the role of circNSD2 protein in breast cancer carcinogenesis, progression, and response to standard of care chemotherapeutics and the mechanistic process by which this protein functions. Additionally, the efficacy of specifically targeting this circRNA as a novel therapeutic approach is being explored. Conclusions: Increasing knowledge of the important functions of circRNAs underlying drug resistance will provide new opportunities for developing efficacious therapeutic strategies and prognostic/predictive biomarkers for TNBC. Citation Format: Xiyin Wang, Michael J. Emch, Xiaojia Tang, Jia Yu, Krishna R. Kalari, Liewei Wang, Matthew P. Goetz, John R. Hawse. The role of circular RNAs in triple negative breast cancer and chemotherapy resistance [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 1689.

Abstract Zinc finger E-box–binding homeobox 1 (ZEB1) is a transcription factor that can promote tumor invasion and metastasis by inducing epithelial-to-mesenchymal transition (EMT). To date, regulation of ZEB1 by RAS/RAF signaling remains unclear, and few studies have examined posttranslation modification of ZEB1, including its ubiquitination. In human colorectal cancer cell lines with RAS/RAF/MEK/ERK activation, an interaction of ZEB1 with the deubiquitinase ubiquitin-specific protease 10 (USP10) was identified whereby USP10 modifies ZEB1 ubiquitination and promotes its proteasomal degradation. Regulation of the USP10–ZEB1 interaction by MEK–ERK signaling was shown whereby constitutive activation of ERK can phosphorylate USP10 at Ser236 to impair its interaction with ZEB1 and enable ZEB1 protein stabilization. Stabilized ZEB1 was shown to promote colorectal cancer metastatic colonization in a mouse tail vein injection model. Conversely, MEK–ERK inhibition blocked USP10 phosphorylation and enhanced the USP10–ZEB1 interaction shown to suppress ZEB1-mediated tumor cell migration and metastasis. In conclusion, we demonstrate a novel function of USP10 in the regulation of ZEB1 protein stability and its ability to mediate tumor metastasis in a preclinical model. Implications: The MEK–ERK-regulated interaction of USP10 with ZEB1 can promote the proteasomal degradation of ZEB1 and thereby suppress its demonstrated ability to mediate tumor metastasis.

Abstract Androgen receptor (AR) is expressed in 80% to 90% of estrogen receptor α–positive (ER+) breast cancers. Accumulated evidence has shown that AR is a tumor suppressor and that its expression is associated with improved prognosis in ER+ breast cancer. However, both a selective AR agonist (RAD140) and an AR inhibitor (enzalutamide, ENZ) have shown a therapeutic effect on ER+ breast cancer, so the potential for clinical application of AR-targeting therapy for ER+ breast cancer is still in dispute. In this study, we evaluated the efficacy of ENZ and RAD140 in vivo and in vitro in AR+/ER+ breast cancer models, characterizing the relationship of AR and ER levels to response to AR-targeting drugs and investigating the alterations of global gene expression and chromatin binding of AR and ERα after ENZ treatment. In the AR-low setting, ENZ directly functioned as an ERα antagonist. Cell growth inhibition by ENZ in breast cancer with low AR expression was independent of AR and instead dependent on ER. In AR-high breast cancer models, AR repressed ERα signaling and ENZ promoted ERα signaling by antagonizing AR. In contrast, RAD140 activated AR signaling and suppressed AR-high tumor growth by deregulating ERα expression and blocking ERα function. Overall, analysis of the dynamic efficacies and outcomes of AR agonist, and antagonist in the presence of different AR and ERα levels reveals regulators of response and supports the clinical investigation of ENZ in selected ER+ tumors with a low AR/ER ratio and AR agonists in tumors with a high AR/ER ratio. Significance: The ratio of androgen receptor to estrogen receptor in breast cancer dictates the response to AR-targeted therapies, providing guidelines for developing AR-directed treatment strategies for patients with breast cancer.

Abstract Understanding the function of non-coding genomic sequence variants represents a challenge for biomedicine. Many diseases are products of gene-by-environment interactions with complex mechanisms. This study addresses these themes by mechanistic characterization of non-coding variants that influence gene expression only after drug or hormone exposure. Using glucocorticoid signaling as a model system, we integrated genomic, transcriptomic, and epigenomic approaches to unravel mechanisms by which variant function could be revealed by hormones or drugs. Specifically, we identified cis-regulatory elements and 3D interactions underlying ligand-dependent associations between variants and gene expression. One-quarter of the glucocorticoid-modulated variants that we identified had already been associated with clinical phenotypes. However, their affected genes were ‘unmasked’ only after glucocorticoid exposure and often with function relevant to the disease phenotypes. These diseases involved glucocorticoids as risk factors or therapeutic agents and included autoimmunity, metabolic and mood disorders, osteoporosis and cancer. For example, we identified a novel breast cancer risk gene, MAST4, with expression that was repressed by glucocorticoids in cells carrying the risk genotype, repression that correlated with MAST4 expression in breast cancer and treatment outcomes. These observations provide a mechanistic framework for understanding non-coding genetic variant-chemical environment interactions and their role in disease risk and drug response.

AbstractAlternative polyadenylation (APA) is a common genetic regulatory mechanism that generates distinct 3′ ends for RNA transcripts. Changes in APA have been associated with multiple biological processes and disease phenotypes. However, the role of hormones and their drug analogs in APA remains largely unknown. In this study, we investigated transcriptome‐wide the impact of glucocorticoids on APA in 30 human B‐lymphoblastoid cell lines. We found that glucocorticoids could regulate APA for a subset of genes, possibly by changing the expression of 142 RNA‐binding proteins, some with known APA‐regulating properties. Interestingly, genes with glucocorticoid‐mediated APA were enriched in viral translation‐related pathways, while genes with glucocorticoid‐mediated expression were enriched in interferon and interleukin pathways, suggesting that glucocorticoid‐mediated APA might result in functional consequences distinct from gene expression. For example, glucocorticoids, a pharmacotherapy for severe COVID‐19, were found to change the APA but not the expression of LY6E, an important antiviral inhibitor in coronavirus diseases. Glucocorticoid‐mediated APA was also cell‐type‐specific, suggesting an action of glucocorticoids that may be unique to immune regulation. We also observed evidence for genotype‐dependent glucocorticoid‐mediated APA (referred to as pharmacogenomic‐alterative polyadenylation quantitative trait loci), providing potential functional mechanisms for a series of common genetic variants that had previously been associated with immune disorders, but without a clear mechanism. In summary, this study reports a series of observations regarding the impact of glucocorticoids on APA, raising the possibility that this mechanism might have implications for both disease pathophysiology and drug therapy.

Neoadjuvant chemotherapy (NAC) remains the cornerstone of the treatment for triple negative breast cancer (TNBC), with the goal of complete eradication of disease. However, for patients with residual disease after NAC, recurrence and mortality rates are high and the identification of novel therapeutic targets is urgently needed. We quantified tyrosine phosphorylation (pTyr)-mediated signaling networks in chemotherapy sensitive (CS) and resistant (CR) TNBC patient-derived xenografts (PDX), to gain novel therapeutic insights. The antitumor activity of SFK inhibition was examined in vivo. Treated tumors were further subjected to phosphoproteomic and RNAseq analysis, to identify the mechanism of actions of the drug. We identified Src Family Kinases (SFKs) as potential therapeutic targets in CR TNBC PDXs. Treatment with dasatinib, an FDA approved SFK inhibitor, led to inhibition of tumor growth in vivo. Further analysis of post-treatment PDXs revealed multiple mechanisms of actions of the drug, confirming the multi-target inhibition of dasatinib. Analysis of pTyr in tumor specimens suggested a low prevalence of SFK-driven tumors, which may provide insight into prior clinical trial results demonstrating a lack of dasatinib antitumor activity in unselected breast cancer patients. Taken together, these results underscore the importance of pTyr characterization of tumors, in identifying new targets, as well as stratifying patients based on their activated signaling networks for therapeutic options. Our data provide a strong rationale for studying SFK inhibitors in biomarker-selected SFK-driven TNBC.

Abstract Background The AMP‐activated protein kinase (AMPK) is a central regulator of energy homeostasis, with deregulation leading to cancer and other diseases. However, how this pathway is dysregulated in cancer has not been well clarified. Methods Using a tandem affinity purification/mass‐spec technique and biochemical analyses, we identified tumor protein D52 (TPD52) as an AMPKα‐interacting molecule. To explore the biological effects of TPD52 in cancers, we conducted biochemical and metabolic assays in vitro and in vivo with cancer cells and TPD52 transgenic mice. Finally, we assessed the clinical significance of TPD52 expression in breast cancer patients using bioinformatics techniques. Results TPD52, initially identified to be overexpressed in many human cancers, was found to form a stable complex with AMPK in cancer cells. TPD52 directly interacts with AMPKα and inhibits AMPKα kinase activity in vitro and in vivo. In TPD52 transgenic mice, overexpression of TPD52 leads to AMPK inhibition and multiple metabolic defects. Clinically, high TPD52 expression predicts poor survival of breast cancer patients. Conclusion The findings revealed that TPD52 is a novel regulator of energy stress‐induced AMPK activation and cell metabolism. These results shed new light on AMPK regulation and understanding of the etiology of cancers with TPD52 overexpression.

Approximately one‐third of patients with metastatic castration‐resistant prostate cancer (CRPC) exhibited primary abiraterone resistance. To identify alternative treatment for abiraterone nonresponders, we performed drug discovery analyses using the L1000 database using differentially expressed genes identified in tumor biopsies and patient‐derived xenograft (PDX) tumors between abiraterone responders and nonresponders enrolled in PROMOTE trial. This approach identified 3 drugs, including topoisomerase II (TOP2) inhibitor mitoxantrone, CDK4/6 inhibitor palbociclib, and pan‐CDK inhibitor PHA‐793887. These drugs significantly suppressed the growth of abiraterone‐resistant cell lines and PDX models. Moreover, we identified 11 genes targeted by all 3 drugs that were associated with worse outcomes in both the PROMOTE and Stand Up To Cancer cohorts. This 11‐gene panel might also function as biomarkers to select the 3 alternative therapies for this subgroup of patients with CRPC, warranting further clinical investigation.

Abstract Background: Androgen receptor (AR) is widely expressed in ER+ breast cancers. Accumulated evidence showed that AR is a tumor suppressor in ERα+ breast cancers. Enzalutamide (ENZ), an AR antagonist, blocks AR signalling, and was reported to inhibit estrogen-driven tumor growth as effectively as tamoxifen. Hence, this study was performed to further examine the mechanism of action of ENZ in ERα+ breast cancers. Methods: The AR gene was knocked out using the CRISPR-Cas9 system. Four crRNAs targeting AR and trans-activating crRNA were co-transfected to Cas9-positive MCF7 and T47D cells, two ERα+ breast cancer cell lines with relatively low AR levels based on RNA seq. ESR1 knock down was performed using a lentivirus system. MCF7 and T47D cells were infected with the ESR1-shRNA lentivirus and selected by puromycin for 2 weeks. Single clones of knockdown or knockout cells were tested using western blot. MCF7 and T47D cells were collected for RNA-seq and ChIP-seq after treated with DMSO or 10µM ENZ for 6 hours in complete medium. All sequences were aligned to the Human Reference Genome (hg38, UCSC). Differential gene expression was performed using limma package and upstream analysis was done in Ingenuity Pathway Analysis software. ChIP peaks were called using MACS2 and overlapping was analyzed using ChIPpeakAnno package in R. Tritium labeled ENZ or estradiol were incubated with ERα overnight. The bound radioactivity was measured using a Beckmann LS 6500 liquidscintillation counter. Results: We found the cell growth inhibition effect of ENZ was little impaired by AR gene knock out in MCF7 and T47D cells as was expected. As ENZ has been shown to inhibit estrogen-driven tumor growth, we performed studies to understand the mechanism of ENZ.We performed RNA-seq and identified 193 differential expressed genes (76 up and 117 down) between DMSO and ENZ-treated MCF7 cells (FC>2 or <0.5, FDR<0.05). Upstream analysis showed that ESR1 might be the upstream regulator of those differentially expressed genes and ESR1 signaling was predicted to be inhibited after ENZ treatment (P=1.69×10-10). To examine the alteration of ERα genomic binding, we performed global ERα ChIP-seq and found 39.58% (8,681/21,933) of ERα-bound sites were abolished by ENZ and the binding intensity was significantly decreased upon ENZ treatment in 12,868 shared peaks. Moreover, radiolabled binding assays showed that ENZ displaced estrogen binding to ERα and directly bound to ERα with a KD= 399.60nM. We further observed a reduction of proliferation inhibition of ENZ in ESR1-KD cells compared with control cells. This suggests that ENZ inhibits ER+ cell growth by displacing E2 from ERα binding sites and subsequently diminishes ERα genomic binding and blocks ERα signaling in ER+ cancers. Conclusion: We found that enzalutamide inhibits cell growth not through AR but by antagonizing ERα activity in MCF7 and T47D ER+ cells. These results reveal ER is the main factor affecting enzalutamide efficacy in these two ERα+ breast cancer cell lines. Further clinical trials are needed to investigate AR antagonists in ERα+ breast cancers. Citation Format: Lixuan Wei, Jia Yu, Huanyao Gao, Huan Zhang, Thanh Nguyen, Yayun Gu, Richard M. Weinshilboum, James N. Ingle, Liewei Wang. Estrogen receptor is a target of enzalutamide in ER+ breast cancer [abstract]. In: Proceedings of the 2021 San Antonio Breast Cancer Symposium; 2021 Dec 7-10; San Antonio, TX. Philadelphia (PA): AACR; Cancer Res 2022;82(4 Suppl):Abstract nr PD1-03.

Abstract Background: Hormone receptor-positive (HR+) breast cancer (BC) accounts for ~60-70% of all BC cases. CDK4/6i combined with endocrine therapy (ET) is the standard of care for first- and second-line therapy of patients with HR+ metastatic breast cancer (MBC) based on substantial improvements in progression free survival (PFS) and overall survival compared to ET alone. Unfortunately, most patients (pts) develop resistance to CDK4/6i. The mechanisms underlying CDK4/6i resistance have not been fully elucidated. Robust in vivo CDK4/6i resistant animal models from pts with de novo and acquired CDK 4/6i resistance are lacking. Patient-derived xenografts (PDXs) are important preclinical tools in oncology to advance our understanding of cancer biology and to evaluate drug response phenotypes. Unfortunately, the establishment of sustainable PDXs from pts with HR+ BC is much more difficult than from HR-negative BC. In order to obtain clinically relevant models that can be utilized to better understand mechanisms of primary and acquired resistance to ET and CDK4/6i, we report the development of PDXs from a prospective study in which percutaneous biopsies were obtained from women receiving palbociclib for the treatment of HR+ advanced or MBC (PROMISE; NCT03281902). Methods: PROMISE (Prospective Study to Evaluate the Role of Tumor Sequencing in Women Receiving Palbociclib for Advanced Hormone Receptor (HR)-Positive Breast Cancer) is a prospective trial which enrolled pts with HR+, HER2-negative MBC who planned to start treatment with palbociclib and letrozole (1st line [FL]) or palbociclib and fulvestrant (2nd line [SL]). For pts enrolled at Mayo Clinic Rochester (MCR), serial tumor biopsies were obtained prior to the start of treatment and at disease progression to establish PDXs. Because baseline take rate was low, we amended the protocol to attempt PDX generation after the 2nd cycle (C2) of therapy. Tissues were received within an hour of excision and were orthotopically implanted in the mammary fat pad of 6~8-week-old nonobese diabetic (NOD)/SCID/IL-2g-receptor null (NSG) mice. Tumor growth was defined as tumor volume ≥ 50 mm3. Tumor take rate was defined as percent of patients with development of at least one stably transplantable (passed at least for two generations) xenograft that was pathologically confirmed as BC. Results: Of the 50 pts enrolled at MCR, 37 were FL pts (2 ineligible) and 13 were SL pts (1 ineligible). Fresh tumor samples obtained pre-treatment from 43 eligible pts (32 FL), after C2 from 17 eligible pts (13 FL), and at disease progression from 10 eligible pts (5 FL) were implanted. Tumor growth was seen in at least one mouse from a pre-treatment specimen in 1 (3.1%) FL pt (confirmed BC) and 4 (36.4%) SL pts (2 confirmed BC); from end of C2 specimens in 2 (15.4%) FL pts and 2 (50%) SL pts (1 confirmed BC); and at progression in 1 (20%) FL pt and 2 (40%) SL pts (1 confirmed BC). The take rate from pre-treatment samples among FL pts was 1/32 (3.1%; 95% CI: 0.08-16.2%). Take rates in other groups are pending histological review. Discussion: In pts receiving palbociclib for HR+ MBC, the take rate using a pre-treatment percutaneous biopsy was low (3.1%) in the FL setting. However, higher growth rates were seen in the SL setting. The observation of higher growth rates at the end of C2 and at progression may relate to the activation of G1/S that occurs with withdrawal of palbociclib, and should be further studied. PDX remains a valuable tool for the evaluation of tumor biology and the development of new therapeutics, especially for those models established from pts with palbociclib and ET resistance. Proteogenomic and sequencing studies are ongoing to fully characterize these PDXs. Citation Format: Jia Yu, Vera J. Suman, Jun He, Jason P. Sinnwell, Ann M. Moyer, Bo Qin, Yayun Gu, Huan Zhang, Jodi M. Carter, Krishna R. Kalari, Brendan P. McMenomy, Minetta C. Liu, Tufia C. Haddad, Kathryn J. Ruddy, Fergus J. Couch, Alvaro Moreno-Aspitia, Donald W. Northfelt, Richard Weinshilboum, Ciara O'Sullivan, Matthew P. Goetz, Liewei Wang. Patient-derived xenografts (PDXs) generated from hormone receptor-positive breast cancer (BC) before and after cyclin-dependent kinase 4 and 6 inhibitors (CDK4/6i) treatment: Initial findings from the PROMISE study [abstract]. In: Proceedings of the 2021 San Antonio Breast Cancer Symposium; 2021 Dec 7-10; San Antonio, TX. Philadelphia (PA): AACR; Cancer Res 2022;82(4 Suppl):Abstract nr P5-01-07.