Barman, Poulami

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.

Abstract Purpose: To determine genetic variants that could differentiate aromatase inhibitors (AIs) efficacy in early stage breast cancer. Experimental Design: We performed a stratified cox-proportional hazards analysis utilizing stratification factors and other covariates examining Breast Cancer Free Interval (BCFI) SNP-treatment interaction in Caucasian patients entered on the MA.27 trial. This analysis involved 4465 patients (2226 on anastrozole and 2239 on exemestane arm) including 245 (121 on anastrozole and 124 events on exemestane arm) with a breast event. Preclinical laboratory studies included luciferase activity, chromatin immunoprecipitation (ChIP) assay, and cell migration assays. Results: We identified 887 SNPs with a p-value <1E-4 that could differentiate anastrozole from exemestane efficacy. We next used GTEx databases to determine whether these SNPs might be cis-eQTLs with nearby genes, and found that 95 SNPs were eQTL with 14 genes. Functional validation of SNP effect in these 14 genes on response to anastrozole and exemestane revealed that 3 SNPs showed genotype-dependent differences between anastrozole and exemestane. Two of the three SNPs, rs1877193 and rs6735923 located upstream of LY75 gene, were associated with higher LY75 gene expression. Our drug-SNP interaction GWAS showed that both SNPs were associated with better BCFI for exemestane compared to anastrozole (HR= 0.447, 0.458). The SNP rs62293499 located downstream of the GPR160 gene was also associated with better BCFI for exemestane compared to anastrozole (HR=0.433). Interestingly, all 3 SNPs were associated with worse outcome (shorter BCFI) in anastrozole treated patients (HR=1.39~1.58), but with longer BCFI in the exemestane treated patients (HR=0.58~0.71) based on our GWAS results in MA.27 trial. However, all 3 SNPs showed no association with BCFI if the two treatment arms, anastrozole and exemestane, were combined (HR=0.92~1.08). Consistently, cells with Ly75 or GPR160 SNP variant genotypes were more sensitive to exemestane compared to anastrozole. Mechanistically, the 3 SNPs regulated estrogen receptor-dependant LY75 and GPR160 expression. LY75 suppression induced epithelial-to-mesenchymal transition (EMT) in breast cancer cell lines, accompanied by increased migratory capacity in vitro. LY75 knockdown also resulted in predominant downregulation of functional pathways such as cell proliferation, while pathways associated with mesenchymal stimulation were generally increased. GPR160 knockdown also resulted in downregulation of cell proliferation pathways. Conclusions: This SNP genotype and AI treatment interaction clinical study revealed unique genetic variants that differentiate anastrozole and exemestane efficacy. The signals were lost in the GWAS analysis when anastrozole and exemestane were combined. Preclinical laboratory studies revealed novel functions of LY75 and GPR160 in breast cancer. These findings represent potential steps towards individualized AI therapy. Citation Format: Junmei Cairns, James N. Ingle, Krishna R. Kalari, Lois E. Shepherd, Matthew J. Ellis, Paul E. Goss, Poulami Barman, Erin E. Carlson, Matthew P. Goetz, Richard M. Weinshilboum, Liewei Wang. The interaction between SNP genotype and aromatase inhibitor treatment response in early breast cancer [abstract]. In: Proceedings of the Annual Meeting of the American Association for Cancer Research 2020; 2020 Apr 27-28 and Jun 22-24. Philadelphia (PA): AACR; Cancer Res 2020;80(16 Suppl):Abstract nr LB-103.

Objectives Based on our previous findings that postmenopausal women with estrone (E1) and estradiol (E2) concentrations at or above 1.3 pg/ml and 0.5 pg/ml, respectively, after 6 months of adjuvant anastrozole therapy had a three-fold risk of recurrence, we aimed to identify a single-nucleotide polymorphism (SNP)-based model that would predict elevated E1 and E2 and then validate it in an independent dataset. Patients and methods The test set consisted of 322 women from the M3 study and the validation set consisted of 152 patients from MA.27. All patients were treated with adjuvant anastrozole, had on-anastrozole E1 and E2 concentrations and genome-wide genotyping. Results SNPs were identified from the M3 genome-wide association study. The best model to predict the E1-E2 phenotype with high balanced accuracy was a support vector machine model using clinical factors plus 46 SNPs. We did not have an independent cohort that is similar to the M3 study with clinical, E1-E2 phenotypes and genotype data to test our model. Hence, we chose a nested matched case-control cohort (MA.27 study) for testing. Our E1-E2 model was not validated but we found the MA.27 validation cohort was both clinically and genomically different. Conclusions We identified a SNP-based model that had excellent performance characteristics for predicting the phenotype of elevated E1 and E2 in women treated with anastrozole. This model was not validated in an independent dataset but that dataset was clinically and genomically substantially different. The model will need validation in a prospective study.

Abstract Purpose: To determine if the degree of estrogen suppression with aromatase inhibitors (AI: anastrozole, exemestane, letrozole) is associated with efficacy in early-stage breast cancer, and to examine for differences in the mechanism of action between the three AIs. Experimental Design: Matched case–control studies [247 matched sets from MA.27 (anastrozole vs. exemestane) and PreFace (letrozole) trials] were undertaken to assess whether estrone (E1) or estradiol (E2) concentrations after 6 months of adjuvant therapy were associated with risk of an early breast cancer event (EBCE). Preclinical laboratory studies included luciferase activity, cell proliferation, radio-labeled ligand estrogen receptor binding, surface plasmon resonance ligand receptor binding, and nuclear magnetic resonance assays. Results: Women with E1 ≥1.3 pg/mL and E2 ≥0.5 pg/mL after 6 months of AI treatment had a 2.2-fold increase in risk (P = 0.0005) of an EBCE, and in the anastrozole subgroup, the increase in risk of an EBCE was 3.0-fold (P = 0.001). Preclinical laboratory studies examined mechanisms of action in addition to aromatase inhibition and showed that only anastrozole could directly bind to estrogen receptor α (ERα), activate estrogen response element-dependent transcription, and stimulate growth of an aromatase-deficient CYP19A1−/− T47D breast cancer cell line. Conclusions: This matched case–control clinical study revealed that levels of estrone and estradiol above identified thresholds after 6 months of adjuvant anastrozole treatment were associated with increased risk of an EBCE. Preclinical laboratory studies revealed that anastrozole, but not exemestane or letrozole, is a ligand for ERα. These findings represent potential steps towards individualized anastrozole therapy.

Abstract Background: The LAR subtype is a genomically distinct subset of TNBC. Using a large cohort of non-metastatic TNBC patients (pts) with long term follow-up, we sought to further characterize the clinicopathologic features and outcomes of LAR vs non-LAR TNBC. Methods: From a cohort of 9982 women with surgically-treated non-metastatic breast cancer, 605 met criteria for TNBC (ER/PR<1% and HER2-negative) by central pathology. RNA extracted from 304 FFPE tumor specimens using the HighPure RNA extraction kit was subjected to TruSeq RNA Access library preparation and sequencing on a HiSeq2500. Adequate RNA was available for 283 pts. Tumors were classified as LAR or non-LAR using a shrunken centroid model, CABAL (Clustering Among BAsal and Luminal androgen receptor). In addition to previously described analyses [Leon-Ferre et al, Breast Cancer Res Treat 2017], immunohistochemical (IHC) androgen receptor (AR) staining was performed and the impact of various parameters on invasive disease-free survival (IDFS) and overall survival (OS) was assessed using Cox proportional hazards models. Results: 58 (20%) tumors were classified as LAR and 225 (80%) as non-LAR. Compared to non-LAR, LAR pts were older (mean age 65 vs 54) and more often postmenopausal (79%vs53%), both p=0.01. Apocrine histology was more common among LAR tumors (21%vs0%), which were also lower grade (grade3: 69%vs95%) and had lower Ki-67 (Ki-67>15%: 64%vs82%), all p<0.01. Additionally, LAR tumors had lower median stromal tumor infiltrating lymphocytes (TILs, 20%vs25%) and were less frequently lymphocyte-predominant [≥50% stromal or intratumoral TILs (19%vs32%)], although neither reached statistical significance. AR IHC was available for 223 of 283 tumors. Median AR IHC score in LAR was 65% (range 0-100%) vs 0% (range 0-90%) in non-LAR. T/N stage, surgery type, and receipt of adjuvant chemotherapy (AdjCT) or radiotherapy were similar between LAR and non-LAR. LAR pts had shorter IDFS and OS compared to non-LAR (5.6 vs 11.8 yrs and 10.8 vs 20.8 yrs, respectively), although this did not reach statistical significance. Test of proportional hazard assumption was not significant for IDFS or OS (p = 0.30 and 0.09). IDFS estimates were numerically higher in LAR vs non-LAR (80.2%vs70.5%,p = 0.92) at 3yrs post-diagnosis; whereas the opposite was true (40.9%vs55.6%,p = 0.07) after 10yrs. OS estimates at 3 and 5yrs were similar between LAR and non-LAR, but at 10yrs OS was inferior in LAR (40.9%vs66.4%,p = 0.24). In a univariate analysis including both LAR and non-LAR, older age, higher N stage, lower TILs and absence of AdjCT were associated with poorer IDFS and OS. In a multivariate analysis, higher N stage and absence of AdjCT remained associated with both poorer IDFS and OS; while lower stromal TILs were associated with poorer IDFS (p=0.01), and with a trend towards poorer OS (p=0.07). Conclusions: LAR TNBCs occurred in older women, were lower grade, and had lower TIL density than nonLAR tumors. While significant differences in IDFS or OS were not demonstrated, LAR pts exhibited a numerically lower risk of a disease event at 3yrs, but higher risk by 10yrs compared to nonLAR pts. In the entire cohort, higher N stage, absence of AdjCT and lower TILs were independently associated with poorer outcomes. Citation Format: Leon-Ferre RA, Polley M-Y, Liu H, Kalari KR, Boughey JC, Liu MC, Cafourek V, Negron V, Ingle JN, Thompson KJ, Tang X, Barman P, Carlson E, Visscher DW, Carter JC, Couch FJ, Goetz MP. Characteristics, outcomes and prognostic factors of luminal androgen receptor (LAR) triple-negative breast cancer (TNBC) [abstract]. In: Proceedings of the 2018 San Antonio Breast Cancer Symposium; 2018 Dec 4-8; San Antonio, TX. Philadelphia (PA): AACR; Cancer Res 2019;79(4 Suppl):Abstract nr P3-08-01.

Anastrozole is a widely prescribed aromatase inhibitor for the therapy of estrogen receptor positive (ER+) breast cancer. We performed a genome‐wide association study (GWAS) for plasma anastrozole concentrations in 687 postmenopausal women with ER+ breast cancer. The top single‐nucleotide polymorphism (SNP) signal mapped across SLC38A7 (rs11648166, P = 2.3E‐08), which we showed to encode an anastrozole influx transporter. The second most significant signal (rs28845026, P = 5.4E‐08) mapped near ALPPL2 and displayed epistasis with the SLC38A7 signal. Both of these SNPs were cis expression quantitative trait loci (eQTL)s for these genes, and patients homozygous for variant genotypes for both SNPs had the highest drug concentrations, the highest SLC38A7 expression, and the lowest ALPPL2 expression. In summary, our GWAS identified a novel gene encoding an anastrozole transporter, SLC38A7, as well as epistatic interaction between SNPs in that gene and SNPs near ALPPL2 that influenced both the expression of the transporter and anastrozole plasma concentrations.

5053 Background: Genomic aberrations associated with resistance/response to AA/P are not known. In a prospective study we assessed whole-exome/RNA-seq based aberrations in CRPC metastatic biopsies for identifying molecular markers associated with primary resistance and response duration. Methods: Sequencing of metastatic biopsies was performed for analyzing molecular aberrations that predict primary resistance (defined as progression at 12-weeks of therapy (non-responders) using PSA, RECIST, bone scan criteria per PCWG2). Gene network analysis was performed in genes mutated more frequently in non-responders and in genes differentially expressed between non-responders and responders using a “risk ratio” (RR) of ≥2. Cox regression models with multiple gene network pathways were used for determining association with time to treatment change (TTTC). Results: Of 92 enrolled pts 82 had complete whole-exome, RNA-seq & 12-week outcome data available for analysis. At 12-weeks 33/82 had progressed. Using a RR of ≥2, 113 genes were more frequently mutated in non-responders & 292 in responders. In non-responders, gene network analysis revealed frequent mutations in Wnt/β-catenin pathway genes; frequent deletion of negative regulators of Wnt pathway ( DKK4, SFRP2, LRP6). Gene expression analyses revealed significantly reduced expression levels of Wnt/β-catenin pathway inhibitors and increased expression levels of cell cycle proliferation (CCP) genes in non-responders. Median study follow up was 32 months during which time 58/82 pts progressed and switched treatments. Median TTTC was 10.1 months (IQR:4.4-24.1). In multivariate analysis CCP scores of ≥50 predicted shorter TTTC (HR = 2.11, 95% CI: 1.17-3.80; p = 0.01). Conclusions: In metastases Wnt/β-catenin pathway activation is associated with primary AA/P resistance and increased CCP with acquired drug resistance. These findings offer molecular based predictive biomarkers in CRPC stage treatment. Clinical trial information: NCT#01953640.

175 Background: Genome and transcriptome aberrations associated with primary resistance to abiraterone acetate/prednisone (AA/P) in mCRPC are not known. In a prospective trial (NCT#01953640) we assessed whole-exome and RNA-seq based aberrations in metastases of CRPC stage patients (pts) for identifying markers associated with primary resistance to AA/P. Methods: Whole-exome and RNA-seq of biopsies from metastases was performed followed by analyses for association between genomic aberrations & primary resistance. Primary resistance was defined by progression on AA/P after 12-weeks of therapy(non-responders) using PSA, RECIST, bone scan and symptom criteria (per PCWG2). Gene network analysis was performed in genes mutated more frequently in non-responders, and also in genes that were differentially expressed between non-responders and responders and a “risk ratio” (RR) was calculated thereof. Results: Between 6/2013 & 8/2015, 92 pts were enrolled of which 82 had complete whole-exome, RNA-seq and 12-week outcome data available for analysis. Median age was 72.5 yrs (IQR: 68.5-78); median PSA-18 ng/ml (IQR: 8.1- 46.6). At 12-weeks 33/82 had progressed. Using the RR of 2 as threshold, we identified 113 and 292 genes that were more frequently mutated in non-responders & responders respectively. Among the 113 genes, OBSCN, ADAM21, LPHN3, DOCK10 ( P = 0.08, RR= inf) & USP42 ( P = 0.16, RR = 5.71) were the top candidates.Gene network analysis revealed that in non-responders, genes involved in the Wnt/β-catenin pathway were frequently mutated and negative regulators of Wnt pathway ( DKK4, SFRP2 & LRP6) were also frequently deleted. Gene expression analyses revealed the expression levels of Wnt/β-catenin pathway inhibitors were significantly reduced while expression levels of cell cycle regulatory genes were significantly increased in non-responders. Conclusions: In this study we observed Wnt/β-catenin pathway activation to be associated with primary AA/P resistance. This finding offers a potential for the development of predictive biomarkers and modulation of targeted pathways to overcome AA/P resistance. Clinical trial information: NCT# 01953640.

249 Background: Evaluating tumor specific pathologic and genetic profiles in mCRPC stage is difficult due to the limited availability of mCRPC tissue. We describe findings from a prospective cohort study performed to obtain concomitant histopathology and genetic information in mCRPC. Methods: Patients with mCRPC initiating abiraterone acetate therapy underwent 2 serial (3 months apart) metastatic (met) site needle core biopsies (NCB1/2). Bone lesions were biopsied using 11-13G core needles and 18G core biopsy devices were used for non-osseous masses. Up to 4 cores were obtained at each biopsy with the 1st core (S1) sent for DNA sequencing, the 2nd (S2) for RNA sequencing, and the 3rd/4th (X3, X4) submitted for xenograft implantation. From each, 1-2 mm segments were separated and formalin-fixed for histopathologic examination (HPE). Results: A total of 54 patients, enrolled between June 2013 and July 2014, underwent 94 NCB (54 NCB1, 40 NCB2), rendering a total of 259 samples (94 S1; 75 S2; 59 X3 and 31 X4) of which 85 (32%) were positive (pos) for tumor by HPE. Positivity for tumor in S1, S2, X3, and X4 cores was 42%, 33%, 22%, and 19%, respectively. At least one core pos for met tumor was observed in 62% NCB1 and 52% NCB2 (overall 52/94; 56%). Met sites biopsied include bone (71), lymph nodes (18), liver (3), penile (1) and pelvic (1) soft tissues. HPE revealed met adenocarcinoma in 44/52, and poorly differentiated carcinoma in 8/52. Gleason grade applied to the met ranged from 3+4 to 5+5. 45/85 pos samples (53%) had > 50% tumor cellularity. Pos NCB in bone lesions was observed in 23/71 (32%), compared to 17/23 (73%) of non-bone sites (p=.0004). 85 S1 samples yielded DNA material, of which 66 (77%) had ≥10% tumoral DNA. 24 cases with negative HPE had ≥10% tumoral DNA; 7 cases with pos HPE had <10 tumoral DNA. In all, 76/94 (81%) NCB yielded tumor material either by HPE or DNA analysis. Conclusions: NCB of mCRPC is feasible and provides adequate tissue for pursuing HPE and sequencing studies. HPE of extracted material correlates with DNA sequencing data and provides complementary information on tumor features. Bone lesions yield significantly less tumoral material than non-osseous sites.