Hoadley, Katherine A.

Abstract Introduction: With emerging antibody drug conjugates (ADC) as treatment of metastatic breast cancer (mBC), it is key to understand the mechanisms that drive resistance to inform optimal ADC selection and sequencing. The AURORA (AUR) US Network was established to characterize paired primary and metastatic (met) samples via multiplatform profiling, with a prospective phase of met tissue collection that included pre and post-ADC samples. Methods: Fresh frozen or FFPE met and archival primary tumors were analyzed across 3 platforms: DNA low pass whole genome and exome sequencing (WES), DNA methylation arrays (Methyl) and RNA sequencing (RNA seq). The molecular cohort included pts with data from ≥1 platform. Pairwise gene expression comparisons were performed in R using t-tests, and pathway analysis with GSEA. Results: A total of 168 pts comprised the molecular cohort: WES, 168; Methyl, 163; RNA seq, 168. Median age at mBC diagnosis was 56 (28-88). 21 pts (12.5%) had known germline BRCA1/2 mutations; 37 (22.0%) presented with de novo mBC. Subtype at mBC (if not available, at primary) diagnosis was: 108 (64.3%) HR+/HER2-, 22 (13.1%) HER2+, 33 (19.6%) TNBC and 5 (3.0%) unknown. Among 139 pts with known subtype in the first AUR met collected, 95 (68.3%) were HR+/HER2-, 17 (12.2%) HER2+ and 27 (19.4%) TNBC. 72 pts received ≥1 ADC with 1 (0-4) median line of chemotherapy for mBC prior to ADC1. As ADC1, 42 pts received HER2 ADC (35 trastuzumab deruxtecan, T-DXd; 7 trastuzumab emtansine), 29 TROP2 ADC (28 sacituzumab govitecan, SG; 1 datopotamab deruxtecan) and 1 NECTIN4 ADC. 22 pts received ≥2 ADC (16 immediate sequence; 6 intervening therapy); 11 T-DXd and 11 SG as ADC2. Median duration of ADC1 was 141 days (6-522); 81 (7-425) for ADC2. Progressive disease was the most common reason for ADC discontinuation, denoting resistance in post-ADC samples. 99 samples were collected pre-ADC1 and 24 post-ADC1, of which 5 were pre-ADC2. ERBB2 expression did not differ between samples after HER2 (n=9) vs non-HER2 ADC (n=8) (p=0.55), nor did TACSTD2 after TROP2 (n=8) vs non-TROP2 ADC (n=9) (p=0.16). In pts with HER2+ subtype at mBC diagnosis or AUR collection who received HER2 ADC, ERBB2 expression significantly decreased in post (n=4) vs pre-treatment (n=7) samples (p=0.043). ERBB2 downregulation was observed in 5 pts with paired pre vs post-HER2 ADC samples. In pts who received TROP2 ADC, TACSTD2 expression did not differ in post (n=10) vs pre-treatment (n=47) samples (p=0.63). GSEA of post-ADC samples revealed significant pathway differences by payload, including actin cytoskeleton in post-microtubule inhibitor (MTi) vs non-MTi ADC samples (p=0.001) and ATM-dependent DNA damage response after topoisomerase I inhibitor (TOP1i) vs non-TOP1i ADCs (p=0.043). Lysosomal vesicle biogenesis (p=0.014) and calcium release (p=0.05) pathways significantly differed between ADC naïve vs treated metastases. ERBB2 hypermethylation was not detected in the cohort. TACST2 promoter methylation was observed in 10 samples (β value >0.3) without association with TROP2 ADC exposure. ERBB2 or TACSTD2 mutations were not identified in the ADC-treated cohort. A novel loss-of-function (LOF) mutation in MED12 was detected in one pt immediately post-T-DXd. MED12 LOF has been reported to promote PARP inhibitor and platinum resistance in BRCA-deficient tumors and may represent a novel resistance mechanism to TOP1i. Conclusions: Multiplatform molecular characterization of mBC treated with ADC revealed potential mechanisms of resistance related to downregulation of target expression, payload (e.g., microtubule formation, DNA damage repair) and lysosomal processing. Additional analyses integrating these multimodal data are ongoing. Strategies to overcome resistance, e.g. modifying ADC target, are being evaluated in clinical trials. Citation Format: A. C. Garrido-Castro, W. C. Nenad, J. M. Balko, T. Hinoue, G. L. Wheeler, B. J. Kelly, M. Singha, B. M. Felsheim, A. R. Michmerhuizen, U. Chandran, A. V. Lee, L. A. Carey, T. A. King, E. Mardis, P. W. Laird, K. V. Giridhar, K. A. Hoadley, C. M. Perou. Molecular characterization of resistance to antibody drug conjugates in metastatic breast cancer: a prospective analysis from the AURORA US Network [abstract]. In: Proceedings of the San Antonio Breast Cancer Symposium 2025; 2025 Dec 9-12; San Antonio, TX. Philadelphia (PA): AACR; Clin Cancer Res 2026;32(4 Suppl):Abstract nr PD3-09.

Abstract Background: Patients with metastatic breast cancer (MBC) often have short survival, and clinical disease management is challenging due to tumor heterogeneity, treatment resistance, and unique contributions of organ-specific metastatic microenvironments. Methods: The retrospective and prospective phases of the AURORA US Metastasis Project have profiled samples from 220 MBC patients using RNA-seq, DNA-seq, and DNA methylation arrays. RNA-seq has been performed on 392 specimens from 143 primary tumors and 249 metastases. In total, there are 162 primary-metastasis pairs from 134 patients (13 patients have 2+ metastases profiled). The 249 metastatic samples include 116 liver, 29 lymph node, 17 lung, 14 bone, 10 brain, and 8 skin metastases. Each specimen was assayed to determine its clinical and molecular subtypes. Here we investigate the molecular causes of clinical and PAM50 subtype switching in primary tumor-metastasis pairs using multi-platform data. Results: Of the patients with paired data in the AURORA dataset, 38/134 (28%) had clinical subtype (HR+/HER2+, HR-/HER2+, HR+/HER2-, TNBC) switches between the primary tumor and paired metastasis. The most common clinical subtype switch was loss of hormone receptor expression (16/38 patients) in the metastatic tumor. Many of these observations have been confirmed by concomitant RNA and/or DNA changes. In addition, 59/134 (44%) patients had PAM50 molecular subtype switches. 26/59 (44%) patients had a Luminal A primary tumor which switched either into Luminal B (16/26) or HER2-Enriched (7/26) subtypes, suggesting recurrent progression mechanisms during metastasis. For example, patient ARRP had a Luminal A primary tumor and a HER2-Enriched bone metastasis. Using DNA clonality data, we identified a clonal mutation in the DNA binding domain of FOXA1 (H247L) as well as mutations in CDH1 (frameshift) and NLRC3 (Q650H) that originated in the primary tumor. A distinct, dominant clone arose in the metastasis, harboring oncogenic activating mutations in PIK3CA (H1047R) and ERBB2 (L869R). There was also an increase in RNA expression of ERBB2 and FGFR4 but no change in ERBB2 copy number. This is in alignment with previous publications demonstrating a shift to the HER2-Enriched subtype, often without gain of HER2 amplification. Interestingly, there were 21 patients that had both clinical and PAM50 subtype switches. Patient ANFC had a HR+/HER2- primary tumor (Luminal A) and a TNBC liver metastasis (HER2-Enriched). Notably, the metastasis lost clinical expression of ER and PR (by IHC), and this was also observed through a decrease in ESR1 and PGR mRNA expression. Using DNA clonality data, we identified 5 new clones in the metastasis that were not present in the primary tumor. Finally, we observed that basal-like primary tumors with a PAM50 subtype switch in the paired metastasis had a significantly lower correlation to the basal-like centroid compared to primary tumors with a basal-like paired metastasis (p=0.0053). This suggests that either these tumors had significant intra-tumoral heterogeneity in the primary tumor and the metastasis originated from a non-basal-like portion of the tumor, or the switches from a basal-like primary tumor are a byproduct of the PAM50 bioinformatics algorithm (i.e. the tumor is very close to 2 subtype centroids). Ongoing analyses will precisely determine which subtype switches can be attributed to genetic or epigenetic causes, organ site biases, and/or bioinformatic methodologies. Conclusions: Using multi-platform molecular analyses, we identified cases of clinical and molecular subtype switching in primary-metastasis pairs. These findings begin to elucidate the underlying molecular drivers of tumor progression and subtype switching in patients with MBC, which may have implications for clinical treatment. Citation Format: A. R. Michmerhuizen, B. M. Felsheim, G. L. Wheeler, B. J. Kelly, M. Singha, T. Hinoue, W. C. Nenad, A. C. Garrido-Castro, J. M. Balko, K. Giridhar, U. Chandran, A. V. Lee, L. A. Carey, T. A. King, K. A. Hoadley, P. W. Laird, E. R. Mardis, C. M. Perou, AURORA US Metastasis Network. Molecular causes of subtype switching in primary and metastatic breast tumor pairs from the AURORA US Project [abstract]. In: Proceedings of the San Antonio Breast Cancer Symposium 2025; 2025 Dec 9-12; San Antonio, TX. Philadelphia (PA): AACR; Clin Cancer Res 2026;32(4 Suppl):Abstract nr PD2-06.

Importance Both tumor-infiltrating lymphocytes (TILs) assessment and immune-related gene expression signatures by RNA profiling predict higher pathologic complete response (pCR) and improved event-free survival (EFS) in patients with early-stage ERBB2/HER2-positive breast cancer. However, whether these 2 measures of immune activation provide similar or additive prognostic value is not known. Objective To examine the prognostic ability of TILs and immune-related gene expression signatures, alone and in combination, to predict pCR and EFS in patients with early-stage ERBB2/HER2-positive breast cancer treated in 2 clinical trials. Design, Setting, and Participants In this prognostic study, a correlative analysis was performed on the Cancer and Leukemia Group B (CALGB) 40601 trial and the PAMELA trial. In the CALGB 40601 trial, 305 patients were randomly assigned to weekly paclitaxel with trastuzumab, lapatinib, or both for 16 weeks. The primary end point was pCR, with a secondary end point of EFS. In the PAMELA trial, 151 patients received neoadjuvant treatment with trastuzumab and lapatinib for 18 weeks. The primary end point was the ability of the HER2-enriched subtype to predict pCR. The studies were conducted from October 2013 to November 2015 (PAMELA) and from December 2008 to February 2012 (CALGB 40601). Data analyses were performed from June 1, 2020, to January 1, 2022. Main Outcomes and Measures Immune-related gene expression profiling by RNA sequencing and TILs were assessed on 230 CALGB 40601 trial pretreatment tumors and 138 PAMELA trial pretreatment tumors. The association of these biomarkers with pCR (CALGB 40601 and PAMELA) and EFS (CALGB 40601) was studied by logistic regression and Cox analyses. Results The median age of the patients was 50 years (IQR, 42-50 years), and 305 (100%) were women. Of 202 immune signatures tested, 166 (82.2%) were significantly correlated with TILs. In both trials combined, TILs were significantly associated with pCR (odds ratio, 1.01; 95% CI, 1.01-1.02; P = .02). In addition to TILs, 36 immune signatures were significantly associated with higher pCR rates. Seven of these signatures outperformed TILs for predicting pCR, 6 of which were B-cell related. In a multivariable Cox model adjusted for clinicopathologic factors, including PAM50 intrinsic tumor subtype, the immunoglobulin G signature, but not TILs, was independently associated with EFS (immunoglobulin G signature–adjusted hazard ratio, 0.63; 95% CI, 0.42-0.93; P = .02; TIL-adjusted hazard ratio, 1.00; 95% CI, 0.98-1.02; P = .99). Conclusions and Relevance Results of this study suggest that multiple B-cell–related signatures were more strongly associated with pCR and EFS than TILs, which largely represent T cells. When both TILs and gene expression are available, the prognostic value of immune-related signatures appears to be superior.

PURPOSE CALGB 40603 ( NCT00861705 ), a 2 × 2 randomized phase II trial, demonstrated that adding carboplatin or bevacizumab to weekly paclitaxel (wP) followed by doxorubicin and cyclophosphamide significantly increased the pathologic complete response (pCR) rate in stage II-III triple-negative breast cancer. We now report long-term outcomes (LTOs) and correlative science end points. PATIENTS AND METHODS The Kaplan-Meier method was used to estimate LTOs in 443 patients who initiated study treatment. Log-rank tests and Cox proportional hazards models evaluated the impact of clinical characteristics, pathologic response, calculated residual cancer burden (RCB) in patients with residual disease (RD), treatment assignment, and dose delivery during wP on LTOs, including event-free survival (EFS). Genomic predictors of treatment response and outcomes were assessed on pretreatment tumor samples by mRNA sequencing. RESULTS Among baseline characteristics, only the clinical stage was associated with LTOs. At a median follow-up of 7.9 years, LTOs were not significantly improved with either carboplatin or bevacizumab, overall or in patients with basal-like subtype cancers by genomic analysis. Patients with pCR (n = 205, 46.3%) had significantly higher 5-year EFS (85.5% v 56.6%, log-rank P < .0001) and overall survival (87.9% v 63.4%, P < .0001) rates compared with patients with RD, even those with RCB class I. Among clinical and genomic features, evidence of immune activation, including tumor-infiltrating lymphocytes and low B-cell receptor evenness, was associated with pCR and improved EFS. CONCLUSION Despite higher pCR rates, neither carboplatin nor bevacizumab appeared to improve LTOs although the study was not powered to assess these secondary end points. pCR was associated with superior LTOs even when compared with minimal RD. Markers of immune activation in pretreatment tumor biopsies were independently associated with higher pCR rates and improved survival.

Abstract Patients with metastatic breast cancer (MBC) typically have short survival times and their successful treatment represents one of most challenging aspects of patient care. This poor prognostic behavior is in part due to molecular features including increased tumor cell clonal heterogeneity, multiple drug resistance mechanisms, and alterations of the tumor microenvironment. The AURORA US Metastasis Project was established with the goal to identify molecular features specifically associated with metastasis. We therefore collected and molecularly characterized specimens from 55 metastatic breast cancer (BC) patients representing 51 primary cancers and 102 metastases. The 153 unique tumors were assayed using RNAseq, tumor/germline DNA exomes and low pass whole genome sequencing, and global DNA methylation microarrays. We found intrinsic molecular subtype differences between primary tumors and their matched metastases to be rare in triple negative breast cancer (TNBC)/Basal-like subtype tumors. Conversely, tumor subtype changes were relatively frequent in estrogen receptor positive (ER+) cancers where ~30% of Luminal A cases switched to Luminal B or HER2-enriched (HER2E) subtypes. Clonal evolution studies identified changes in expression subtype coincident with DNA clonality shifts, especially involving HER2 amplification and/or the HER2E expression subtype. We further found evidence for ER-mediated downregulation of genes involved in cell-cell adhesion in metastases. Microenvironment differences varied according to tumor subtype where ER+/Luminal metastases had lower fibroblast and endothelial cell content, while TNBC/Basal-like metastases showed a dramatic decrease in B cells and T cells. In 17% of metastatic tumors, we identified DNA hypermethylation and/or focal DNA deletions near HLA-A that were associated with its’ significantly reduced expression, and with lower immune cell infiltrates. We also identified low immune cell features in brain and liver metastases when compared to other metastatic sites, even within the same patient. These findings have direct implications for the treatment of metastatic breast cancer patients with immune- and HER2-targeting therapies and suggest potential novel therapeutic avenues for the improvement of outcomes for some patients with MBC.

Abstract Background: Tumor infiltrating lymphocytes (TILs) play an important, but incompletely understood role in chemotherapy response and prognosis. In breast cancer, there appear to be distinct immune responses by subtype, but most studies have used limited numbers of protein markers or bulk sequencing of RNA to characterize immune response, in which spatial organization cannot be assessed. Methods: To identify immune phenotypes of Basal-like vs. Luminal breast cancer we used the GeoMx® (NanoString) platform to perform digital spatial profiling (DSP) of immune-related proteins in tumor whole sections and tissue microarrays (TMA). Visualization of CD45, CD68, or pan-Cytokeratin by immunofluorescence was used to select regions of interest (ROIs) in formalin-fixed paraffin embedded tissue sections. 44 antibodies representing stromal markers and multiple immune cell types were applied to quantify the tumor microenvironment. Results: In whole tumor slides, immune hot spots (CD45+) had increased expression of many immune markers, suggesting a diverse and robust immune response. In epithelium-enriched areas, immune signals were also detectable and varied by subtype, with Regulatory T cell (Treg) markers (CD4, CD25, FOXP3) being higher in Basal-like vs. Luminal breast cancer. Extending these findings to TMAs with more patients (n=75), we confirmed subtype-specific immune profiles, including enrichment of Treg markers in Basal-likes. Conclusion: This work demonstrated that immune responses can be detected in epithelium-rich tissue, and that TMAs are a viable approach for obtaining important immunoprofiling data. In addition, we found that immune marker expression is associated with BC subtype, suggesting possible prognostic or targetable differences. Citation Format: Andrea Walens, Linnea T. Olsson, Xiaohua Gao, Alina M. Hamilton, Erin L. Kirk, Stephanie M. Cohen, Bentley R. Midkiff, Yongjuan Xia, Mark E. Sherman, Nana Nikolaishvili-Feinberg, Katherine A. Hoadley, Melissa A. Troester, Benjamin C. Calhoun. Protein-based immune profiles of basal-like vs luminal breast cancers [abstract]. In: Proceedings of the AACR Virtual Special Conference on Tumor Heterogeneity: From Single Cells to Clinical Impact; 2020 Sep 17-18. Philadelphia (PA): AACR; Cancer Res 2020;80(21 Suppl):Abstract nr PO-090.

PURPOSE CALGB 40601 assessed whether dual versus single human epidermal growth factor receptor 2 (HER2) –targeting drugs added to neoadjuvant chemotherapy increased pathologic complete response (pCR). Here, we report relapse-free survival (RFS), overall survival (OS), and gene expression signatures that predict pCR and survival. PATIENTS AND METHODS Three hundred five women with untreated stage II and III HER2-positive breast cancer were randomly assigned to receive weekly paclitaxel combined with trastuzumab plus lapatinib (THL), trastuzumab (TH), or lapatinib (TL). The primary end point was pCR, and secondary end points included RFS, OS, and gene expression analyses. mRNA sequencing was performed on 264 pretreatment samples. RESULTS One hundred eighteen patients were randomly allocated to THL, 120 to TH, and 67 to TL. At more than 7 years of follow-up, THL had significantly better RFS and OS than did TH (RFS hazard ratio, 0.32; 95% CI, 0.14 to 0.71; P = .005; OS hazard ratio, 0.34; 95% CI, 0.12 to 0.94; P = .037), with no difference between TH and TL. Of 688 previously described gene expression signatures, significant associations were found in 215 with pCR, 45 with RFS, and only 22 with both pCR and RFS (3.2%). Specifically, eight immune signatures were significantly correlated with a higher pCR rate and better RFS. Among patients with residual disease, the immunoglobulin G signature was an independent, good prognostic factor, whereas the HER2-enriched signature, which was associated with a higher pCR rate, showed a significantly shorter RFS. CONCLUSION In CALGB 40601, dual HER2-targeting resulted in significant RFS and OS benefits. Integration of intrinsic subtype and immune signatures allowed for the prediction of pCR and RFS, both overall and within the residual disease group. These approaches may provide means for rational escalation and de-escalation treatment strategies in HER2-positive breast cancer.

Abstract Background: It has become increasingly clear that effective treatment of metastatic breast cancer (MBC) requires an in-depth understanding of the molecular differences between primary tumors and metastases. The AURORA US Network was established to collect primary breast cancer-metastasis pairs for multi-platform genomic profiling in order to identify the molecular drivers of metastatic disease. AURORA US has both a retrospective and prospective phase. This is the first report of the retrospective phase. Methods: Archived tissue samples from the primary tumor and at least one distant metastasis were retrospectively collected from 83 MBC patients. Following internal quality assessment, samples from 55 pts, including 105 distinct metastatic lesions, were subject to DNA low pass whole genome and exome sequencing, DNA methylation arrays, and RNA sequencing. Early analyses of these multi-platform data include: DNA methylation, tumor gene expression and microenvironmental signatures, somatic and germline variants, DNA copy number changes, and structural variants between breast primaries and matched metastases. Results: Median age at diagnosis was 49 years (25-76); 32 (58%) were Stage I or II at presentation, 27 (49%) had a family history of breast cancer, and 20 (36%) had a second breast primary. Median disease-free interval before developing MBC was 2 years (range 0-36, 5 patients presented with Stage IV). Median overall survival from initial presentation was 4 years (range 0-37). Median survival after developing MBC was 1 year (range 0-13), with a median of three treatments. Primary tissue samples were banked from 1977-2017 and metastases were banked from 1999-2017. Clinical phenotypes of the primaries included 27 HR+ (49%), 15 triple negative (TNBC, 27%), and 11 HER2+ (20%, 12 missing HER2 status). Intrinsic subtype distribution of the primaries included 17 Basal-like (31%), 17 Luminal A (31%), 7 Normal-like (13%), 5 HER2-enriched (9%), and 1 Luminal B, with 8 pending. All metastases from the Basal-like cases remained Basal-like, while metastases from luminal primaries tended to gain HER2-Enriched subtype features (5/18, p = 0.01). Overall, we identified significant metastasis-enriched alterations in metabolism pathways, an increase in proliferation, and the loss of differentiation signatures and immune infiltrates with progression; the latter being the most pronounced in brain metastases. The most frequent somatic mutations in this cohort were in TP53, NCOR1, and RUNX1. Interestingly, ERBB2, EGFR, and ATM were also mutated in ≥10% of the tumors sequenced. In almost all cases, CpG island hypermethylation was clonally present in the primary tumor and persisted stably in the majority of metastatic lesions. Promoter CpG island hypermethylation was also identified in some metastatic lesions at JAM3, an important cellular adhesion molecule,and this was accompanied by reduced mRNA expression. CONCLUSIONS: Collection of banked primary and metastatic tissue pairs identified a young MBC cohort with a high frequency of breast cancer family history and second breast primaries. Molecular characterization of luminal tumor pairs highlighted acquisition of aggressive traits including increased proliferation and loss of differentiation in the metastases. In contrast, basal-like pairs remained relatively unchanged, except for the loss of immune activation. Ongoing analyses to be presented include clonal heterogeneity and phylogeny, novel metastasis signature discovery, gene fusion, and endogenous retrovirus detection. Citation Format: Tari A King, Minetta C Liu, Marni B McClure, Toshinori Hinoue, Benjamin J Kelly, Chad J Creighton, Jay Bowen, Kristen Leraas, Robyn T Burns, Sara Coppens, Salma Rezk, Amy L Garrett, Justin M Balko, Joel S Parker, Ben H Park, Ian Krop, Carey Anders, Katherine A Hoadley, Julie Gastier-Foster, Mothaffar F Rimawi, Rita Nanda, Nancy U Lin, Claudine Isaacs, P. Kelly Marcom, Anna Maria Storniolo, Fergus J Couch, Elaine R Mardis, Adrian V Lee, Uma Chandran, Peter W Laird, Susan G Hilsenbeck, Larry Norton, Andrea L Richardson, W. Fraser Symmans, Lisa A Carey, Antonio C Wolff, Nancy E Davidson, Charles M Perou, the AURORA US Network. Multiplatform analysis of matched primary and metastatic breast tumors from the AURORA US Network [abstract]. In: Proceedings of the 2019 San Antonio Breast Cancer Symposium; 2019 Dec 10-14; San Antonio, TX. Philadelphia (PA): AACR; Cancer Res 2020;80(4 Suppl):Abstract nr GS3-08.

Abstract RNA sequencing detects estrogen receptor alpha gene (ESR1) fusion transcripts in estrogen receptor positive (ER+) breast cancer but their role in disease pathogenesis remains unclear. Herein we examined multiple in-frame and out-of-frame ESR1 fusions and found only two, both identified in advanced endocrine treatment resistant disease, encoded stable and functional in-frame fusion proteins. In both examples, ESR1-e6&gt;YAP1 and ESR1-e6&gt;PCDH11X, the N-terminal, DNA binding and dimerization motifs encoded by exons 2-6 were fused to C-terminal sequences from the partner gene. Functional properties included estrogen-independent growth, constitutive expression of ER target genes, anti-estrogen resistance, induction of cellular motility in vitro and the development of lung metastasis in vivo. Chromatin immunoprecipitation and RNA sequencing experiments showed both fusions uniquely activated a metastasis-associated transcriptional program. ESR1-e6&gt;YAP1 and ESR1-e6&gt;PCDH11X-induced growth remained sensitive to a CDK4/6 inhibitor, palbociclib, and a patient-derived xenograft (PDX) expressing the ESR1-e6&gt;YAP1 fusion was also responsive. Transcriptionally active ESR1 fusions therefore trigger both endocrine therapy resistance and metastatic progression explaining the association with fatal disease progression, although CDK4/6 inhibitor treatment is predicted to be effective. Citation Format: Jonathan T. Lei, Jieya Shao, Jin Zhang, Michael Iglesia, Doug W. Chan, Jin Cao, Meenakshi Anurag, Purba Singh, Xiaping He, Yoshimasa Kosaka, Ryoichi Matsunuma, Robert Crowder, Jeremy Hoog, Chanpheng Phommaly, Rodrigo Goncalves, Susana Romalho, Raquel M. Peres, Nindo Punturi, Cheryl Schmidt, Alex Bartram, Eric Jou, W V. Lai, Oliver Hampton, Anna Rogers, Ethan Tobias, Poojan Parikh, Sherri R. Davies, Shunqiang Li, Cynthia X. Ma, Vera Suman, Kelly K. Hunt, Mark A. Watson, Katherine A. Hoadley, E A. Thompson, Xi Chen, Shyam M. Kavuri, Chad J. Creighton, Christopher A. Maher, Charles M. Perou, Svasti Haricharan, Matthew J. Ellis. Functional and therapeutic significance of ESR1 gene fusions in breast cancer [abstract]. In: Proceedings of the American Association for Cancer Research Annual Meeting 2018; 2018 Apr 14-18; Chicago, IL. Philadelphia (PA): AACR; Cancer Res 2018;78(13 Suppl):Abstract nr 5240.

Abstract Background. Dysregulation of the estrogen receptor gene (ESR1) is an established mechanism of inducing endocrine therapy resistance. We previously discovered a chromosomal translocation event generating an estrogen receptor gene fused in-frame to C-terminal sequences of YAP1 (ESR1-YAP1) that contributed to endocrine therapy resistance in estrogen receptor positive (ER+) breast cancer models. This study compares functional, transcriptional, and pharmacological properties of additional ESR1 gene fusion events of both early stage (ESR1-NOP2) late stage (ESR1-YAP1 and ESR1-PCDH11x) breast cancers to gain a better understanding of therapeutic resistance and metastasis. Understanding the role of ESR1 fusions in inducing metastasis is critical, since the primary cause of death in breast cancer patients is through metastasis to distant sites. Methods. RNA-seq screens identified ESR1 fusions from early and late stage, endocrine therapy resistant breast tumor samples. Functional experiments were conducted using ER+ breast cancer cell lines, xenograft, and PDX models to test the ability of ESR1 fusions to induce therapeutic resistance and metastasis. ChIP-seq and RNA-seq were performed to examine transcriptional properties and differential gene expression induced by the fusions which directed subsequent pharmacological experiments with a CDK4/6 inhibitor. Results. ESR1-YAP1 and ESR1-PCDH11x promoted estrogen-independent and fulvestrant-resistant growth in vitro and induced greater tumor growth and increased metastatic capacity to the lungs of xenografted mice. In contrast, the ESR1-NOP2 fusion was sensitive to low estrogen conditions in vitro, and did not promote tumor growth. RNA-seq profiling revealed E2F targets pathway as the most highly enriched pathway induced by the ESR1 fusions. IHC revealed higher levels of pRb in ESR1-YAP1 and ESR1-PCDH11x xenograft tumors and subsequent CDK4/6 inhibition completely blocked tumor growth in an ESR1-YAP1 PDX model. Integrating RNA-seq with ChIP-seq data, we discovered a set of EMT and metastasis genes bound by all ESR1 fusions and WT-ER, but whose expression was strongly and uniquely up-regulated only by the ESR1-YAP1 and ESR1-PCDH11x fusions. These studies also revealed gained sites bound only by the ESR1-YAP1 and ESR1-PCDH11x fusions, not bound by WT-ER nor ESR1-NOP2. Genes mapping to these sites have a role in metastatic biology and were highly up-regulated by the YAP1 and PCDH11x fusions, potentially mediated by long range transcriptional activation. Conclusion. ESR1-YAP1 and ESR1-PCDH11x are driver fusions that occur in drug-resistant, advanced stage breast cancer and are a new class of recurrent somatic mutation that can cause acquired endocrine therapy resistance, yet can be treated with CDK4/6 inhibition. These driver fusions also confer increased metastatic ability through their ability to drive expression of genes that contribute to EMT and metastasis. In contrast, ESR1-NOP2 did not produce functional protein and appears to be a passenger event. These studies may provide pre-clinical rationale for targeting ESR1 translocated breast tumors, since the presence of an ESR1 driver fusion places a patient in a therapeutic category where none of the currently available endocrine therapies are likely to be effective. Citation Format: Lei JT, Shao J, Zhang J, Iglesia M, Chan DW, Cao J, Anurag M, Singh P, Haricharan S, Kavuri SM, Matsunuma R, Schmidt C, Kosaka Y, Crowder R, Hoog J, Phommaly C, Goncalves R, Ramalho S, Rodrigues-Peres RM, Lai W-C, Hampton O, Rogers A, Tobias E, Parikh P, Davies S, Ma C, Suman V, Hunt K, Watson M, Hoadley KA, Thompson A, Perou CM, Creighton CJ, Maher C, Ellis MJ. ESR1 gene fusions drive endocrine therapy resistance and metastasis in breast cancer [abstract]. In: Proceedings of the 2017 San Antonio Breast Cancer Symposium; 2017 Dec 5-9; San Antonio, TX. Philadelphia (PA): AACR; Cancer Res 2018;78(4 Suppl):Abstract nr PD8-03.