Citations

BACKGROUND Most functionally‐significant single nucleotide polymorphisms (SNPs) are outside of open reading frames. Mechanisms by which these SNPs contribute to gene regulation are poorly understood. Aryl hydrocarbon receptor (AHR) is a nuclear receptor that regulates drug metabolism, xenobiotic response and immune response. Our previous data indicate that SNPs outside response elements can impact transcription factor binding and gene regulation. We hypothesized that SNPs near but outside AHR response elements (AhREs) might impact AHR binding and gene regulation. METHODS We have collected dense SNP‐genotypes across 300 lymphoblastoid cell lines (LCLs) – a powerful tool for generating and testing pharmacogenomic hypotheses. Eight LCLs were used to generate hypotheses about AHR SNP‐dependent transcription utilizing RNA‐ and ChIP‐seq, which were tested in additional LCLs using qRT‐PCR and ChIP after AHR agonist or antagonist treatment. RESULTS 69 genes were highly significantly differentially expressed after agonist or antagonist treatment and 80% of those genes had an AhRE within 50kb. Five SNP‐gene pairs have been functionally validated in additional LCLs (Table 1). CONCLUSION SNPs that are outside of AhREs can alter AHR binding and gene transcription. This is a general phenomenon that may have large‐scale influence on gene transcriptional regulation and pharmacogenomics. This abstract is from the Experimental Biology 2018 Meeting. There is no full text article associated with this abstract published in The FASEB Journal .

Objectives Betaine‐homocysteine methyltransferase (BHMT) catalyzes the remethylation of homocysteine. BHMT in the human liver is responsible for up to 25% of homocysteine metabolism. Knowledge of the relationship between genetic polymorphisms and BHMT phenotypes could enhance our understanding of the role of this reaction in drug metabolism. Methods Levels of BHMT enzyme activity and immunoreactive protein were measured in 268 human liver biopsies. Twelve tag BHMT SNPs from our own gene resequencing and HapMap data were used to genotype DNA from these samples. Genotype‐phenotype association analysis was then performed. BHMT enzyme activity averaged 0.077 ± 0.028 units/μg, with a 12 fold variation. Results BHMT protein levels averaged 6.9 ± 3.3 units/μg, with a 50‐fold variation. Three SNPs (rs41272270, rs16876512, and rs6875201), located 28kb upstream, in the 5′‐UTR and in intron 1 of BHMT , respectively, were significantly correlated with both phenotypes (p<10 −8 for enzyme activity and p<10 −5 for protein level). rs41272270 was associated with altered protein binding during electrophoresis mobility shift assay. Conclusions Three tag SNPs were significantly associated with human hepatic BHMT enzyme activity and immunoreactive protein. These observations may help us understand the regulation of BHMT in the human liver and its role in drug metabolism. Supported by NIH U01 GM61388, R01 GM28157, R01 CA132780.

Abstract Purpose: Patient-derived xenografts (PDX) are a research tool for studying cancer biology and drug response phenotypes. While engraftment rates are higher for tumors with more aggressive characteristics, it is uncertain whether engraftment is prognostic for cancer recurrence. Patients and Methods: In a prospective study of patients with breast cancer treated with neoadjuvant chemotherapy (NAC) with taxane ± trastuzumab followed by anthracycline-based chemotherapy, we report the association between breast cancer events and PDX engraftment using tumors derived from treatment naïve (pre-NAC biopsies from 113 patients) and treatment resistant (post-NAC at surgery from 34 patients). Gray test was used to assess whether the cumulative incidence of a breast cancer event differs with respect to either pre-NAC PDX engraftment or post-NAC PDX engraftment. Results: With a median follow-up of 5.7 years, the cumulative incidence of breast cancer relapse did not differ significantly according to pre-NAC PDX engraftment (5-year rate: 13.6% vs. 13.4%; P = 0.89). However, the incidence of a breast event was greater for patients with post-NAC PDX engraftment (5-year rate: 50.0% vs. 19.6%), but this did not achieve significance (P = 0.11). Conclusions: In treatment-naïve breast cancer receiving standard NAC, PDX engraftment was not prognostic for breast cancer recurrence. Further study is needed to establish whether PDX engraftment in the treatment-resistant setting is prognostic for cancer recurrence.

Aromatase inhibitors (AIs) are the treatment of choice for hormone receptor–positive early breast cancer in postmenopausal women. None of the third‐generation AIs are superior to the others in terms of efficacy. We attempted to identify genetic factors that could differentiate between the effectiveness of adjuvant anastrozole and exemestane by examining single‐nucleotide polymorphism (SNP)‐treatment interaction in 4,465 patients. A group of SNPs were found to be differentially associated between anastrozole and exemestane regarding outcomes. However, they showed no association with outcome in the combined analysis. We followed up common SNPs near LY75 and GPR160 that could differentiate anastrozole from exemestane efficacy. LY75 and GPR160 participate in epithelial‐to‐mesenchymal transition and growth pathways, in both cases with SNP‐dependent variation in regulation. Collectively, these studies identified SNPs that differentiate the efficacy of anastrozole and exemestane and they suggest additional genetic biomarkers for possible use in selecting an AI for a given patient.

Objectives We previously discovered that the single nucleotide polymorphisms (SNP) rs9940645 in the ZNF423 gene regulate ZNF423 expression and serve as a potential biomarker for response to selective estrogen receptor modulators (SERMs). Here we explored pathways involved in ZNF423-mediated SERMs response and drugs that potentially sensitize SERMs. Methods RNA sequencing and label-free quantitative proteomics were performed to identify genes and pathways that are regulated by ZNF423 and the ZNF423 SNP. Both cultured cells and mouse xenograft models with different ZNF423 SNP genotypes were used to study the cellular responses to metformin. Results We identified ribosome and AMP-activated protein kinase (AMPK) signaling as potential pathways regulated by ZNF423 or ZNF423 rs9940645 SNP. Moreover, using clustered regularly interspaced short palindromic repeats/Cas9-engineered ZR75-1 breast cancer cells with different ZNF423 SNP genotypes, striking differences in cellular responses to metformin, either alone or in the combination of tamoxifen, were observed in both cell culture and the mouse xenograft model. Conclusions We found that AMPK signaling is modulated by the ZNF423 rs9940645 SNP in estrogen and SERM-dependent fashion. The ZNF423 rs9940645 SNP affects metformin response in breast cancer and could be a potential biomarker for tailoring the metformin treatment.

Abstract Background: Androgen receptor (AR) is expressed in 60%-90% of breast cancers. The role of AR in breast cancer largely depends on estrogen receptor (ER) status and remains controversial. In ER+ cancers, AR expression is associated with improved prognosis. Enzalutamide (ENZ), an AR antagonist, impairs AR signaling, inhibits ER+/AR+ breast cancer cell proliferation and has been shown to mitigate resistance to anti-estrogen therapies. A recent study identified that RAD140, an AR agonist, suppressed the growth of ER+/AR+ breast cancer via stimulating AR signaling, resulting in down regulation of ERα. Both AR agonists and antagonists are effective in treating ER+ breast cancers. Identifying the subgroup of patients who most likely will benefit from an individual drug is important in clinical practice. Here, we sought to characterize the relative AR to ER levels and their relationships to drug response in ER+ breast cancers.Methods: We evaluated AR and ER expression levels among 68 samples with ER+ breast cancers. On each tumor, IHC staining for AR and ER were performed and results were scored by multiplying the percentage of positive cells by the intensity. TCGA-RPPA (reverse phase protein array) dataset was used to verify AR protein distribution in ER+ cancers. We created MCF7 and T47D cells with doxycycline-inducible AR and ER expression system to manipulate the relative AR to ER ratios and determined the antitumor activity of ENZ and RAD140 by cell proliferation assays. The levels of AR and ER in cells were measured by western blot. Linear regression was used to test the association between dose-response area under curve (AUC) and AR, ER levels. Drug preference was modeled against AR/ER expression levels (AR/ER ratios) using logistic regression.Results: Among 68 ER+ breast cancers, 69.12% were AR-positive. More than 2/3 cases (48 of 68) had more ER than AR expressed. The AR to ER ratios varied from 0 to 6. In the TCGA cohort, 84.15% of 347 ER+ patients had AR/ER ratios less than 1.The range of AR to ER ratios in TCGA dataset were comparable with our data (0.004 to 4.210). In our cell line models, the AR/ER ratios were controlled between 0.19 to 4.05. We found that the AUC of RAD140 was negatively associated with AR protein levels (P=0.008) and AR/ER ratios (P=0.013), and not significantly associated with ER expression levels. On the other hand, the AUC of ENZ was significantly negatively associated with ER protein (P=0.0016) but postiviely associated with AR/ER ratio (p=0.037). Preferred treatment comparing efficacy of the two drug can be best determined at extremes of AR/ER ratios. Using clinically relevant dosages, our model predicted that ENZ (10µM) would be a preferred treatment choice and have a better treatment efficacy compared with RAD140 when the AR/ER ratio was ≤ 0.42, whereas RAD140 (1µM) would be the preferred choice with a better treatment efficacy compared with ENZ when AR/ER ratios were ≥ 3.1 The efficacy preference of the two drugs are equivocal for AR/ER ratios between 0.42 and 3.10. Conclusion: We developed preclinical models using AR and ER expression levels to predict AR-targeting drug response. The results support the use of RAD140 in AR high patients and those with an AR/ER ratio &gt;3.10, and enzalutamide in AR low patients and those with an AR/ER ratio &lt;0.42. Equipoise on choice of drug was found for AR/ER ratios of 0.42-3.10. RAD140 and enzalutamide are compelling candidates for monotherapy or combination with anti-estrogen therapies in ER+/AR+ breast cancer. Future clinical validation of the models and therapeutic effect is warranted. Citation Format: Lixuan Wei, Huanyao Gao, Jia Yu, Duan Liu, Huan Zhang, Thanh Nguyen, Marie R Passow, Jodi M Carter, Richard M Weinshilboum, James N Ingle, Liewei Wang. Computational modeling of androgen receptor (AR) and estrogen receptor as predictive biomarkers of response to AR agonists and antagonists [abstract]. In: Proceedings of the 2020 San Antonio Breast Cancer Virtual Symposium; 2020 Dec 8-11; San Antonio, TX. Philadelphia (PA): AACR; Cancer Res 2021;81(4 Suppl):Abstract nr PS5-36.

Abstract Background: Cyclin dependent 4/6 kinase inhibitors (CDK4/6i) and endocrine therapy (ET) have improved progression-free survival (PFS) and overall survival in hormone-receptor (HR)-positive metastatic breast cancer (MBC), but progression of disease is inevitable. Serum thymidine kinase-1 (TK1) is a secreted marker of proliferation that is prognostic in patients (pts) with HR-positive, HER2-negative MBC and may be predictive of ET and CDK 4/6i response. PROMISE (NCT0281902) is a prospective study enrolling women with HR-positive MBC starting palbociclib (Pb) + letrozole (L) (1st line) or Pb + fulvestrant (2nd line). We undertook a comprehensive “omic” assessment of blood, tumor, urine and the fecal microbiome in order to identify novel genomic variants and pathways associated with an early decline in TK1 (measured after 2 months) and PFS. Additionally, patient derived xenografts/organoids were generated at baseline and progression to test new therapeutic approaches to overcome resistance to CDK4/6i and ET. We report the initial association between the baseline genomic landscape and baseline TK1 levels. Methods: FFPE tumor biopsies were obtained for DNA/RNA sequencing (TempusTM) and blood samples for TK1 (Divitum® assay, Biovica) were collected pretreatment (pre-Pb) and after 2 cycles of Pb + ET (post-Pb2). Both whole-exome (exome capture) sequencing (WES) and RNA-Seq used the Integrated DNA Technologies xGen Exome Research Panel v1.0 capture kit. TK1+ disease was defined as &gt; 200 Du/L and TK1- disease as below limit of detection up to 200 Du/L. We tested the association between genomic and transcriptomic characteristics with baseline TK1 data in pretreatment samples where both WES and RNA-seq and TK1 was available. The data were analyzed using bioinformatics pipelines for somatic and germline mutations and copy number alterations. The current analysis focuses on baseline 1st-line pre-Pb omics data in conjunction with baseline TK1 levels. The database was locked for analysis on 5/29/2020. Results: Thirty-three pts (median age: 59 yrs.) were evaluable, with paired samples for TK1 in 32. Six pts had TK1+ disease pre-Pb and post-Pb2. Twenty-two pts had TK1- disease pre-Pb and post-Pb2. Four pts had a decrease in TK1 after 2 cycles of treatment that altered the classification from TK1+ to TK1-. Both baseline RNA seq and serum TK1 (n=16) were available for 4 TK1+ and 12 TK1- pts. In this group, 476 genes were differentially regulated (398 upregulated; 78 downregulated). Pathway analysis demonstrated enrichment in complement and coagulation cascade pathway, PPAR signaling pathway, and metabolism-related pathways related to up-regulation of CYP and UGT gene families. Further testing for the association of WES data with baseline TK1+ (n=8) and TK1- (n=16) disease demonstrated somatic copy number variations on chromosomes 6, 11, 12 and 15. CDK4 copy number gains were observed in 3/8 TK1+ pts and 0/16 TK1- pts. We also observed that somatic mutations (LOH, copy number and/or SNV/INDELs) were more prevalent in the TK1+ compared to the TK1- pre-Pb group in several cancer-associated genes (FAS [p=0.06] PTEN, PIK3CB, NAB2, SOX9 and FAT1 [p=0.08], TP53, and MAP2K4 [p=0.22]). Conversely, we also noted that 6/7 pts with GATA3 mutations had TK1- disease (p=0.23). Conclusions: Using a comprehensive “omics” approach, our data suggest that a secreted biomarker of proliferation (TK1) obtained prior to initiating CDK4/6i and ET for the first line treatment of HR+ MBC is associated with established and novel genes and pathways associated with prognosis of pts receiving ET and CDK 4/6i. Analysis of on-treatment (after 2 cycles) tumor RNA seq and its association with change in TK1 as well as data from the 2nd-line cohort will be presented at the meeting. Citation Format: Ciara C O Sullivan, Krishna R Kalari, Vera J Suman, Peter T Vedell, Ann Moyer, Abraham D Eyman Casey, Jason Sinnwell, Xiaojia Tang, Kevin Thompson, Alvaro Moreno-Aspitia, Donald W Northfelt, Minetta C Liu, Tufia C Haddad, Saranya Chumsri, Prema Peethambaram, Kathryn J Ruddy, Karthik V Giridhar, Roberto A Leon-Ferre, Adrian Nordmark, Mattias Bergqvist, Brendan P McMenomy, Richard M Weinshilboum, Liewei Wang, Matthew P Goetz. Novel genomic variants and pathways associated with baseline serum thymidine kinase 1 levels in HR-positive HER2-negative MBC patients commencing palbociclib and letrozole [abstract]. In: Proceedings of the 2020 San Antonio Breast Cancer Virtual Symposium; 2020 Dec 8-11; San Antonio, TX. Philadelphia (PA): AACR; Cancer Res 2021;81(4 Suppl):Abstract nr PS5-24.

Abstract Background: The tamoxifen (TAM) metabolite, ENDX, demonstrated promising antitumor activity in endocrine resistant breast cancer (BC) in both phase I and phase II settings. Furthermore, ENDX resulted in superior in vivo antitumor activity compared to TAM and letrozole in aromatase-expressing aromatase inhibitor-sensitive and resistant MCF7AC1 models. Recently, we identified protein kinase C beta 1 (PKCβ1), which regulates cell proliferation and tumorigenic transformation, as a novel target of ENDX. ENDX-bound PKCβ1 at concentrations achieved in phase I/II ENDX studies (100-300 nM). In contrast, TAM binding to PKCβ1 occurred at concentrations 7-10 folds higher (2 μM) than achievable with TAM 20 mg/day dosing. However, the clinical relevance of targeting PKCβ1 kinase activity is unclear, since drugs that target PKCβ1 (enzastaurin) have been ineffective in BC and other solid tumors. Therefore, we sought to understand how ENDX altered PKCβ1 and to further compare and contrast ENDXs effects to that of PKCβ1 kinase inhibition in ERα+ BC. Methods: The effects of PKCβ1 silencing and ENDX treatment on gene expression was analyzed by RNAseq in MCF7AC1 cells. The impact of PKCβ1-silencing on cell cycle was evaluated by flow cytometry. Protein expression of cell cycle regulators in PKCβ1 and ENDX-treated MCF7AC1 and T47D cells were compared to TAM and enzastaurin in vitro and to letrozole, TAM or control in vivo. The effects of PKCβ1 and drugs on growth were analyzed by cell proliferation assays. PRKCB gene amplification was assessed in primary tumors using TCGA data and in metastatic tumors using whole-exome sequencing data from patients enrolled in the PROMISE study (NCT 03281902). Results: RNAseq analysis revealed E2F targets and G2M checkpoints as the top hallmark genesets significantly downregulated in both PKCβ1-silenced and ENDX-treated MCF7AC1 cells. Flow cytometry demonstrated that PKCβ1 silencing increased G1 and reduced S phases of the cell cycle. Western blot analyses of PKCβ1-silenced MCF7AC1 and T47D cells displayed reduced protein levels of the cell cycle regulators Cyclin D1, Retinoblastoma (Rb), phospho-RbS807/811, CDK4, Chk1 and E2F1 that regulate G1/S transition. While short term ENDX (48 hours) treatment did not alter PKCβ1 levels, prolonged in vitro ENDX treatment profoundly reduced PKCβ1 protein levels and the aforementioned cell cycle regulators, faithfully replicating PKCβ1 silencing effects. In contrast, enzastaurin had no impact on proliferation or cell cycle proteins in either model. Consistent with this finding, ENDX, but not TAM or letrozole, reduced protein levels of ERα and cell cycle regulators in vivo. Overexpression of PKCβ1 induced TAM, but not ENDX, resistance and had little impact on responsiveness to enzastaurin. While PRKCB gene amplification was uncommon in newly diagnosed ERα+/HER2- BC (5%, TCGA), PRKCB was amplified in 40% of metastatic ERα+/HER2- BC (PROMISE study). Conclusion: We have confirmed the relevance of a new ENDX target, PKCβ1, in ERα+/HER2- BC. While targeting PKCβ1 kinase activity elicited no anticancer effects in ERα+ cells, PKCβ1 downregulation, either by siRNA or ENDX, resulted in profound ERα turnover, reduced protein levels of essential cell cycle mediators and profoundly inhibited cell proliferation. Furthermore, PKCβ1 protein expression is associated with TAM, but not ENDX, resistance, a finding whose clinical relevance is further magnified by identification of PRKCB amplification in metastatic ERα+ BC, confirming its potential importance in progression. Efforts are currently underway to elucidate the mechanistic basis for ENDX-induced PKCβ1 and ERα degradation and the contribution of these effects to the superior antitumor activity of ENDX in ERα+ BC. Citation Format: Swaathi Jayaraman, Mary J Kuffel, Krishna R Kalari, Kevin J Thompson, Xiaojia Tang, Vera J Suman, Elizabeth S Bruinsma, Ciara C O'Sullivan, Liewei Wang, Richard Weinshilboum, James N Ingle, John R Hawse, Matthew P Goetz. Antitumor activity of Z-endoxifen (ENDX) is mediated via PKCβ1-dependent ERα loss and cell cycle arrest in ERα-positive breast cancer [abstract]. In: Proceedings of the 2020 San Antonio Breast Cancer Virtual Symposium; 2020 Dec 8-11; San Antonio, TX. Philadelphia (PA): AACR; Cancer Res 2021;81(4 Suppl):Abstract nr PD8-04.

Abstract Introduction: The immune microenvironment of high-risk HR+ breast cancer (BC) is poorly understood. BEAUTY is a prospective neoadjuvant chemotherapy (NAC) study of stage I-III BC patients treated with neoadjuvant weekly taxane followed by anthracycline-based chemotherapy. Among clinical luminal B BC (St. Gallen criteria) from BEAUTY, we used high-plex digital spatial profiling (DSP) to characterize the intra-epithelial tumor and stromal immune landscapes to 1) assess the impact of NAC on these landscapes and 2) identify immune biomarkers predictive of response to NAC. Methods: The tissue set included FFPE sections of resected tumors from 35 patients (median 51y; range: 21-71y) with clinically-defined luminal B BC (ER &gt; 10%/grade 2/Ki67 ≥ 15% or ER &gt; 10%/grade 3), and 16 paired pre-NAC biopsies. Nanostring GeoMX DSP was used to quantitate 58 immune and BC biomarker proteins in intra-epithelial, cytokeratin-positive tumor segments and adjacent stromal (cytokeratin-negative/SYTO13 (nuclear stain)-positive) segments. Data were normalized using the geometric mean of two negative controls (RbIgG and MsIgG1). A general linear model with negative binomial identified differentially-expressed (DE) proteins in pre/post-NAC tumors. Based on DE protein data and biologic function, a targeted protein subset (N=19) was evaluated in pre-NAC biopsies for associations between protein abundance and residual cancer burden (RCB) class (0-II vs. III) (Wilcoxon-rank sum test: p-value &lt; 0.025 considered significant) or RCB ‘breast only’ index (Spearman correlation). Results: Comparing tumor segments in pre-NAC specimens, intra-epithelial segments were predictably enriched in cytokeratin, ER, PR and Ki-67, but also CD127 and NY-ESO-1; whereas stromal segments were enriched in proteins associated with T cell subsets (e.g. CD3, CD4, CD8), macrophages (CD68 and CD163), antigen presentation/dendritic cells (CD11c, HLA-DR) and immune checkpoint proteins (PD-L1, PD-L2, B7-H3, TIM-3, VISTA) among others. While NAC did not alter the spatial distribution of proteins (intra-epithelial vs. stromal segments), it markedly attenuated the immune landscape, with decreased abundance of functionally-diverse immune proteins (e.g. CD45, CD3, CD4, CD127, granzyme B, CTLA4; STING, B7-H3, CD11c, and CD68, log2FC: 0.27-1.52 p &lt; 0.05), including low abundance proteins PD-1, PD-L1, PD-L2, CD20, and OX40L (log 2FC: 0.15-1.16, p &lt; 0.05). Both PR and Ki-67 decreased in post-NAC tumors whereas ER was not significantly altered (p &lt; 0.05). CD8, CD14, CD163, HLA-DR, IDO-1 and TIM-3 were not significantly altered by NAC. In the pre-NAC biopsies (which had subsequent residual tumor burdens of RCB class 0 or I (n=1 each), II (n=6), and III (n=8), and median RCB ‘breast only’ index of 3.39 (range: 0.00-38.02), 19 proteins were used for RCB analysis (CD3, CD4, CD8, CD14, CD34, CD44, CD45, CD68, CD127, CD163, CTLA4, granzyme B, STING, B7-H3, fibronectin, Ki-67, HLA-DR, SMA and TIM3). Among them, stromal CD127 was significantly higher in RCB class III than RCB class 0-II (p=0.021) and RCB ‘breast only’ index was positively correlated with intra-epithelial granzyme B (rho = 0.61; p=0.012), and negatively correlated with intra-epithelial Ki-67 (rho= -0.71; p=0.0022). Conclusion: In this series of clinical luminal B BC, NAC markedly attenuated the tumoral immune landscapes with a small set of “NAC-resistant” immune proteins. Among a targeted set, stromal CD127 was significantly higher in RCB III, and RCB breast-only index was positively correlated with intra-epithelial granzyme B. These data provide insight into the impact of NAC on HR+ BC, and identify potential immune biomarkers to predict response to neoadjuvant chemotherapy. Citation Format: Jodi M Carter, Judy C Boughey, Jun He, Vera J Suman, Xue Wang, Jennifer M Kachergus, Krishna R Kalari, Liewei Wang, Richard Weinshilboum, Ann M Moyer, Sarah A McLaughlin, Alvaro Moreno-Aspitia, Donald W Northfelt, Richard J Gray, James N Ingle, E. Aubrey Thompson, Matthew P Goetz. Neoadjuvant chemotherapy selectively alters spatially-defined immune landscapes in clinical luminal B HR+/HER2- breast cancers: Analysis of the breast cancer genome guided therapy study (BEAUTY) [abstract]. In: Proceedings of the 2020 San Antonio Breast Cancer Virtual Symposium; 2020 Dec 8-11; San Antonio, TX. Philadelphia (PA): AACR; Cancer Res 2021;81(4 Suppl):Abstract nr PD7-05.