Carter, Jodi M.

Breast cancer is a heterogenous disease with variability in tumor cells and in the surrounding tumor microenvironment (TME). Understanding the molecular diversity in breast cancer is critical for improving prediction of therapeutic response and prognostication. High-plex spatial profiling of tumors enables characterization of heterogeneity in the breast TME, which can holistically illuminate the biology of tumor growth, dissemination and, ultimately, response to therapy. The GeoMx Digital Spatial Profiler (DSP) enables researchers to spatially resolve and quantify proteins and RNA transcripts from tissue sections. The platform is compatible with both formalin-fixed paraffin-embedded and frozen tissues. RNA profiling was developed at the whole transcriptome level for human and mouse samples and protein profiling of 100-plex for human samples. Tissue can be optically segmented for analysis of regions of interest or cell populations to study biology-directed tissue characterization. The GeoMx Breast Cancer Consortium (GBCC) is composed of breast cancer researchers who are developing innovative approaches for spatial profiling to accelerate biomarker discovery. Here, the GBCC presents best practices for GeoMx profiling to promote the collection of high-quality data, optimization of data analysis and integration of datasets to advance collaboration and meta-analyses. Although the capabilities of the platform are presented in the context of breast cancer research, they can be generalized to a variety of other tumor types that are characterized by high heterogeneity.

Abstract Benign breast disease (BBD) is an established risk factor for development of breast cancer (BC). However, little is known about the prevalence of somatic genetic alterations at the time of BBD diagnosis and relationship to subsequent BC. To characterize DNA variants in BBD and investigate association with BC risk, we isolated DNA from formalin-fixed paraffin-embedded (FFPE) BBD biopsies and assessed mutations using a panel of 93 breast cancer predisposition genes (Qiagen GeneRead). The study set consisted of 120 BBD patients from three age-matched groups defined based on 16 years of BC follow-up: BBD-controls (cancer-free after 16 years), BBD-ER+, or BBD-ER- breast cancer (within 16 years). Bioinformatics quality checks included comparisons of paired fresh-frozen and FFPE replicates to reduce false variant calls induced by FFPE artifacts. Gene-level variant burden differences among study groups were assessed using 12 statistical approaches, adjusted for patient age, year of biopsy, and histology. To investigate possible immune implications of increased variant burden, extent of leukocytes in the BBD biopsies was assessed using CD45 immunohistochemistry. To assess variant profile changes during the progression from BBD to BC, we also performed targeted sequencing of 12 BBD-matched subsequent tumor samples (six ER+ and six ER-) developed within 5 years of BBD diagnosis. Gene-level association results showed consistently skewed distributions suggesting a higher variant burden in BBD controls: among the top 10 genes associated with case-status (p<0.05), all had BC-risk odds ratio estimates less than 1.0. Immune infiltrates, assessed by CD45 H-scores, were positively associated with overall variant burden (r=0.42, p=0.0031), and noticeably higher in cancer-free controls (r=0.50, p=0.005). Comparison of mutation variants in 12 BBD cases that developed cancer within 5 years of diagnosis and their subsequent tumors revealed that more than half of the non-SNP variants were shared between the BBD and the subsequent tumor (54.6%), whereas less than one quarter were unique to the tumor (20.6%) or unique to the BBD (24.9%). Through BBD-BC paired variant analysis, three subjects were found with newly emergent clones uniquely appearing in the tumors and four subjects had loss of clones in the tumor that were present in the BBD. These results show that somatic variant burden is lower in patients with BBD that progress to BC, mutation profiles of BC progression are heterogeneous, and BBD tissues with reduced variant burden also show reduced immune infiltrate. This implicates a potential protective role for immune activation in BBD. Future directions will define the characteristics of the potentially protective immune response and its role in variant heterogeneity in a larger sample set, integrated with paired RNA sequencing to better understand the pathways involved and the relationship to immune infiltration. Citation Format: Stacey J. Winham, Julie M. Cunningham, Yuanhang Liu, Aditya V. Bhagwate, Ethan P. Heinzen, Samantha J. McDonough, Melody L. Stallings-Mann, Robert A. Vierkant, Jodi M. Carter, Mark E. Sherman, Derek C. Radisky, Amy C. Degnim, Chen Wang. Somatic variant burdenin benign breast disease and association with risk of subsequent breast cancer [abstract]. In: Proceedings of the American Association for Cancer Research Annual Meeting 2021; 2021 Apr 10-15 and May 17-21. Philadelphia (PA): AACR; Cancer Res 2021;81(13_Suppl):Abstract nr 2274.

Abstract Purpose: Programmed death ligand 1 [PD-(L)1]-targeted therapies have shown modest survival benefit in triple-negative breast cancer (TNBC). PD-L1+ microenvironments in TNBC are not well characterized and may inform combinatorial immune therapies. Herein, we characterized clinicopathologic features, RNA-based immune signatures, and spatially defined protein-based tumor–immune microenvironments (TIME) in early-stage PD-L1+ and PD-L1− TNBC. Experimental Design: From a large cohort of chemotherapy-naïve TNBC, clinicopathologic features, deconvoluted RNA immune signatures, and intraepithelial and stromal TIME (Nanostring GeoMX) were identified in subsets of PD-L1+ and PD-L1− TNBC, as defined by FDA-approved PD-L1 companion assays. Results: 228 of 499 (46%) TNBC were PD-L1+ (SP142: ≥1% immune cells-positive). Using PD-L1 22C3, 46% had combined positive score (CPS) ≥ 1 and 16% had CPS ≥10. PD-L1+ TNBC were higher grade with higher tumor-infiltrating lymphocytes (TIL; P < 0.05). PD-L1 was not associated with improved survival following adjustment for TILs and other variables. RNA profiles of PD-L1+ TNBC had increased dendritic cell, macrophage, and T/B cell subset features; and decreased myeloid-derived suppressor cells. PD-L1+ stromal and intraepithelial TIMEs were highly enriched in IDO-1, HLA-DR, CD40, and CD163 compared with PD-L1-TIME, with spatially specific alterations in CTLA-4, Stimulator of Interferon Genes (STING), and fibronectin. Macrophage- and antigen presentation–related proteins correlated most strongly with PD-L1 protein. Conclusions: In this early-stage TNBC cohort, nearly 50% were PD-L1+ (SP142 companion assay) while 16% were PD-L1+ with the 22C3 companion assay. PD-L1+ TNBC had specific myeloid-derived and lymphoid features. Spatially defined PD-L1+ TIME were enriched in several clinically actionable immune proteins. These data may inform future studies on combinatorial immunotherapies for patients with PD-L1+ TNBC. See related commentary by Symmans, p. 5446

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.

548 Background: sTILs are a well-established prognostic and predictive biomarker in patients with operable TNBC receiving pre or postoperative systemic therapy. We1 and others2,3 have also shown that sTILs are prognostic in patients who did not receive adjuvant chemotherapy. Here, we detail the outcomes of systemically untreated patients with N0 TNBC according to sTIL score. We focused on the N0 subset as a group of patients who may be candidates for future prospective therapy de-escalation trials. Methods: From a clinically annotated cohort of 605 patients with centrally confirmed TNBC (ER/PR < 1% and HER2 negative) with long-term outcomes data, we identified 182 patients treated with locoregional therapy only (breast surgery +/- radiation therapy and no chemotherapy). The clinicopathological characteristics of this cohort have previously been published1. In this analysis, we report the 5- and 10-year invasive disease-free survival (iDFS) and overall survival (OS) rates of patients with N0 TNBC according to sTIL levels. IDFS and OS were defined as per the STEEP classification and estimated using the Kaplan–Meier method. Comparisons of the survival distributions between groups were assessed by the log-rank test. sTILs were assessed as a continuous parameter according to the International TIL Working Group guidelines. For comparisons of outcomes between groups, tumors were classified as lymphocyte-predominant TNBC (defined as containing ≥50% sTILs) vs non-lymphocyte-predominant ( < 50% sTILs). Results: Of 182 systemically untreated patients, 149 (82%) were N0 and most (78%) were post-menopausal. T stage distribution was T1: 68%, T2: 28%, T3/4: 4%. Among N0 patients, 31 (21%) had lymphocyte-predominant TNBC, and in this group the 5-year iDFS and OS were 89% (95% CI 76-100) and 96% (95% CI 89-100), while the 10-year iDFS and OS were 89% (95% CI 76-100) and 87% (95% CI 73-100), respectively. In contrast, outcomes for patients with non-lymphocyte predominant TNBC were significantly worse. For this group, 5-year iDFS and OS were 62% (95% CI 53-73) and 78% (95% CI 71-86) while the 10-year iDFS and OS were 45% (95% CI 36-58) and 66% (95% CI 68-76), respectively ( log-rank p = 0.02 for iDFS and log-rank p = 0.03 for OS). Conclusions: sTIL quantification identifies a subset of patients with early-stage N0 TNBC with an exceedingly good prognosis, even in the absence of adjuvant chemotherapy. These data provide support for the evaluation of sTILs as part of prospective investigation of systemic therapy de-escalation strategies in N0 TNBC. References:1Leon-Ferre et al, Breast Cancer Res Treat (2018) 167:89-99 2Park et al, Ann Oncol (2019) 12:1941-1949 3De Jong et al, ESMO 2020

Abstract Human tissue samples commonly preserved as formalin-fixed paraffin-embedded (FFPE) tissues after diagnostic or surgical procedures in the clinic represent an invaluable source of clinical specimens for in-depth characterization of signaling networks to assess therapeutic options. Tyrosine phosphorylation (pTyr) plays a fundamental role in cellular processes and is commonly dysregulated in cancer but has not been studied to date in FFPE samples. In addition, pTyr analysis that may otherwise inform therapeutic interventions for patients has been limited by the requirement for large amounts of frozen tissue. Here we describe a method for highly sensitive, quantitative analysis of pTyr signaling networks, with hundreds of sites quantified from one to two 10-μm sections of FFPE tissue specimens. A combination of optimized magnetic bead–based sample processing, optimized pTyr enrichment strategies, and tandem mass tag multiplexing enabled in-depth coverage of pTyr signaling networks from small amounts of input material. Phosphotyrosine profiles of flash-frozen and FFPE tissues derived from the same tumors suggested that FFPE tissues preserve pTyr signaling characteristics in patient-derived xenografts and archived clinical specimens. pTyr analysis of FFPE tissue sections from breast cancer tumors as well as lung cancer tumors highlighted patient-specific oncogenic driving kinases, indicating potential targeted therapies for each patient. These data suggest the capability for direct translational insight from pTyr analysis of small amounts of FFPE tumor tissue specimens. Significance: This study reports a highly sensitive method utilizing FFPE tissues to identify dysregulated signaling networks in patient tumors, opening the door for direct translational insights from FFPE tumor tissue banks in hospitals.

Abstract Background: Benign breast disease (BBD) is a risk factor for breast cancer (BC); however, little is known about the genetic alterations present at the time of BBD diagnosis, and how these relate to risk of incident BC. Methods: A subset of a long-term BBD cohort was selected to examine DNA variation across three BBD groups (42 future estrogen receptor-positive (ER+) BC, 36 future estrogen receptor-negative (ER-) BC and 42 controls cancer-free for at least 16 years post-BBD). DNA extracted from archival formalin fixed, paraffin-embedded (FFPE) tissue blocks was tested and analyzed for presence of DNA alterations using a targeted panel of 93 BC-associated genes. To address artifacts frequently observed in FFPE tissues (e.g., C>T changes), we applied three filtering strategies based on alternative allele frequencies and nucleotide substitution context. Gene-level associations were performed using two burden tests, SKAT-O and logistic regression, and adjusted for clinical and technical covariates including epithelial percentage, histologic impression, patient age and year of BBD biopsy. Results: After filtering, the variant frequency of SNPs in our samples was highly consistent with population allele frequencies reported in 1KG/ExAC (0.986, p <1e-16). Ten genes were found to be associated with later cancer status by four of 12 association methods, with nominal p-values <0.05: MED12, MSH2, BRIP1, PMS1, GATA3, MUC16, FAM175A, EXT2, MLH1 and TGFB1. Additionally, analysis revealed 11 and 19 gene-level associations with p-values <0.05 and OR<1 for all cancer cases and ER- cases, respectively. Comparing between the three case groups, BBD ER+ cases were closer to controls in mutation profile, while BBD ER- cases were distinct. Notably, the variant burden was significantly higher in controls than in either ER+ or ER- cases. CD45 expression was associated with mutational burden (p=0.003).Conclusions: Somatic mutations were more frequent in benign breast tissue from women who did not develop cancer, opening questions of clonal diversity or immune-mediated restraint on future cancer development. CD45 expression was positively associated with mutational burden, most strongly in controls. Further studies in both normal and premalignant tissues are needed to better understand the role of somatic gene mutations and their contribution to future cancer development.

Abstract Background: Women with benign breast disease (BBD) are at increased risk of breast cancer (BC). Classic studies based on film-based mammographic screening and pathology diagnosis of surgical biopsies conducted in the 1980s established a hierarchy of increasing BC risk: non-proliferative (NP) BBD, proliferative BBD without atypia (PDWA) and atypical hyperplasia (AH). Given changes in epidemiological BC risk factors and introduction of percutaneous core needle biopsy (CNB) in mid-1990s, and later, digital mammography, we hypothesized that the patient characteristics and relative frequency of BBD diagnoses have changed over time. Accordingly, we performed a longitudinal analysis of the frequency of patient characteristics and BBD diagnoses in the Mayo BBD cohort. Methods: Utilizing the Mayo Clinic Surgical and Pathology Indices, women ages 18 to 85 who had a BBD biopsy between 1/1/67 and 12/31/13 were identified. Breast pathologists reviewed biopsies masked to diagnoses of incident BC diagnosed in follow-up. Demographic characteristics and BC events were obtained by query of institutional data sources and participant surveys. Trends were evaluated for the following eras: 1: pre-mammogram (1967-1981), 2: pre-CNB (1982-1992), 3: CNB Transition (1993-2001), and 4: CNB (2002-2013). Demographics were formally compared across eras using chi-square tests for categorical variables and analyses of variance (ANOVAs) for continuous variables. Results: From 1967-2013, the cohort includes 19,582 unique women with BBD. The frequency of CNB increased from eras 1-4: 0.04%, 0.6%, 51.3 %, and 88.9%, respectively. Mean age at BBD diagnosis was younger in era 1 (48.0 years) vs eras 2-4 (53.2, 52.0, and 51.8, respectively, p&lt;0.001). The percentage of biopsies diagnosed as PDWA increased from era 1-4 (25.7%, 34.3%, 35.2%, 46.2%, p&lt;0.001), as did the percentage with AH (2.4%, 5.1%, 8.6%, 12.3%, p&lt;0.001). Over eras 1-4, the percentages of women with a strong family history of BC increased (9.9%, 12.7%, 17.1%, and 29.0%, p&lt;0.001) as did mean BMI (24.8, 26.4, 27.4, and 28.6, p&lt;0.001). With a median follow-up of 10.9 years, 1,719 breast cancers have developed, with increasing proportion of noninvasive (DCIS-only) disease across eras 1-4: 15.0%, 21.1%, 21.2%, and 33.2%, p&lt; 0.001. Conclusions: Analysis of this large, single institution BBD cohort for the 46 year period 1967-2013 demonstrates that BC risk factors among BBD patients has changed over time, with subjects demonstrating increasing age, BMI, and family history, and that the percentages of BBD classified as PDWA and AH have increased. Impact on BC risk will be further investigated. Citation Format: Amy Degnim, Karthik Ghosh, Jodi M Carter, Robert A Vierkant, Matthew R Jensen, Stacey J Winham, Tanya L Hoskin, Marlene Frost, Teresa M Allers, Denise L Gehling, Mindy J Kern, Laura M Pacheco-Spann, Celine M Vachon, Derek C Radisky, Daniel W Visscher, Mark E Sherman. Benign breast disease: Temporal trends from 1967 to 2013 [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 PS7-11.

Abstract Background: Growing data established the pivotal role of preexisting immune response in triple negative breast cancer (TNBC). Conventionally, preexisting immune response can be evaluated by quantifying tumor infiltrating lymphocytes mainly in the stroma or gene expression analysis from the whole tumor section. Due to technical challenges with these conventional methods, limited data regarding specific subtypes and spatial distribution of these immune infiltrates are currently available. Methods: NanoString IO360 gene expression analysis and Digital Spatial Profiling (DSP) were used. DSP was used to quantify 29 immune-related proteins in stromal and tumor-enriched segments from 44 TNBC samples from the FinXX trial (NCT00114816) and 335 samples from the Mayo Clinic (MC) cohort of centrally reviewed TNBC (Leon-Ferre BCRT 2018). In FinXX trial, 22 patients with recurrence and 22 patients without recurrence were included. In MC cohort, 217/335 patients received adjuvant chemotherapy while 118 patients had surgery only without adjuvant chemotherapy. Regions were segmented based on pancytokeratin staining. The general linear model was used for statistical analysis of differential expression with recurrence free survival (RFS) as a categorical variable (recur yes or no). Kaplan-Meier (KM) estimates and Cox regression models were also used for analysis. Results: In the FinXX trial, there were 12 out of 29 proteins in tumor epithelial segments (intraepithelial) which were significantly expressed at higher levels among patients who were free of recurrence. These proteins include Beta-2 microglobulin, CD11c, CD20, CD40, CD56, CD8, Granzyme B, HLA-DR, ICOS, PD-L1, PD-L2, and TGFB1. In contrast, merely 5 out of 29 proteins in stromal segments were significantly differentially expressed in these 2 groups of patients. Granzyme B, IDO1, PD-L1, and PD-L2 in stroma were significantly higher and SMA was significantly lower in patients without recurrence. Using Cox regression models, intraepithelial CD56, CD40, and HLA-DR were significantly associated with outcome. When comparing between highest and lowest intraepithelial protein expression by tertile, intraepithelial CD56 (HR 0.12, 95%CI 0.03-0.39, p &lt; 0.001), CD40 (HR 0.13, 95%CI 0.04-0.46, p = 0.002), and HLA-DR (HR 0.24, 95%CI 0.06-0.89, p = 0.032) were significantly associated with improved outcome. However, expression of these same proteins in stroma was not associated with outcome. Using KM estimates, intraepithelial CD56 (p &lt; 0.0001), CD40 (p = 0.0006), and HLA-DR (p = 0.013) were also significantly associated with improved outcome. Nonetheless, RNA expression of these proteins by IO360 from whole tumor sections were not significantly associated with outcome (CD56 p = 0.27, CD40 p = 0.21, HLA-DR p = 0.48). Similar findings with DSP were observed in MC TNBC cohort. Comparing between the highest and lowest quartiles, there were significantly fewer patients who developed recurrence with high protein expression of intraepithelial CD56 (p &lt; 0.001), CD40 (p = 0.002), and HLA-DR (p = 0.006). Conclusions: Using an in-depth analysis with spatially defined context, we identify that there were numerically more intraepithelial immune-related proteins associated with outcome compared to proteins in stroma. Specifically, intraepithelial CD56, CD40, and HLA-DR were significantly associated with improved outcome in both FinXX and MC TNBC cohorts. However, neither expression of these proteins in stroma nor RNA expression from whole tumor were associated with outcome. Our study highlights the impact of spatial biology and the importance of evaluating each potential biomarker in a spatially defined manner. Support: W81XWH-15-1-0292, BCRF 19-161, P50CA116201-9, P50CA015083 Citation Format: Saranya Chumsri, Jodi M. Carter, Yaohua Ma, Douglas Hinerfeld, Heather Ann Brauer, Sarah Warren, Torsten O. Nielsen, Karama Asleh, Heikki Joensuu, Edith A. Perez, Roberto A. Leon-Ferre, David W. Hillman, Judy C. Boughey, Minetta C. Liu, James N. Ingle, Krishna R. Kalari, Fergus J. Couch, Keith L. Knutson, Matthew P. Goetz, E. A. Thompson. Spatially defined immune-related proteins and outcome in triple negative breast cancer in the FinXX trial and Mayo Clinic cohort [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 PS6-02.

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.