Targetable ERBB2 mutation status is an independent marker of adverse prognosis in estrogen receptor positive, ERBB2 non-amplified primary lobular breast carcinoma: a retrospective in silico analysis of public datasets ========================================================================================================================================================================================================================= * Sasagu Kurozumi * Mansour Alsaleem * Cíntia J. Monteiro * Kartikeya Bhardwaj * Stacey E. P. Joosten * Takaaki Fujii * Ken Shirabe * Andrew R. Green * Ian O. Ellis * Emad A. Rakha * Nigel P. Mongan * David M. Heery * Wilbert Zwart * Steffi Oesterreich * Simon J. Johnston ## Abstract **Background** Invasive lobular carcinoma (ILC) accounts for 10-15% of primary breast cancers and is typically estrogen receptor alpha positive (ER+) and *ERBB2* non-amplified. Somatic mutations in *ERBB2/3* are emerging as a tractable mechanism underlying enhanced human epidermal growth factor 2 (HER2) activity. We tested the hypothesis that therapeutically targetable *ERBB2/3* mutations in primary ILC of the breast associate with poor survival outcome in large public datasets. **Methods** We performed *in silico* comparison of *ERBB2* non-amplified cases of ER+ stage I-III primary ILC (*N*=279) and invasive ductal carcinoma (IDC, *N*=1,301) using METABRIC, TCGA and MSK-IMPACT information. Activating mutations amenable to HER2-directed therapy with neratinib were identified using existing functional data from *in vitro* cell line and xenograft experiments. Multivariate analysis of 10-year overall survival (OS) with tumor size, grade and lymph node status was performed using a Cox regression model. Differential gene expression analyses by *ERBB2* mutation and amplification status was performed using weighted average differences and an *in-silico* model of response to neratinib derived from breast cancer cell lines. **Results** ILC tumors comprised 17.7% of all cases in the dataset but accounted for 47.1% of *ERBB2*-mutated cases. Mutations in *ERBB2* were enriched in ILC *versus* IDC cases (5.7%, *N*=16 *vs*. 1.4%, *N*=18, *p*<0.0001) and clustered in the tyrosine kinase domain of HER2. *ERBB3* mutations were not enriched in ILC (1.1%, *N*=3 *vs*. 1.8%, *N*=23; *p*=0.604). Median OS for patients with *ERBB2*-mutant ILC tumors was 66 months *versus* 211 months for *ERBB2* wild-type (*p*=0.0001), and 159 *vs*. 166 months (*p*=0.733) for IDC tumors. Targetable *ERBB2* mutational status was an independent prognostic marker of 10-year OS – but only in ILC (hazard ratio, HR=3.7, 95% CI 1.2–11.0; *p*=0.021). Findings were validated using a novel *ERBB2* mutation gene enrichment score (HR for 10-year OS in ILC=2.3, 95% CI 1.04–5.05; *p*=0.040). **Conclusions** Targetable *ERBB2* mutations are enriched in primary ILC and their detection represents an actionable strategy with the potential to improve patient outcomes. Biomarker-led clinical trials of adjuvant HER-targeted therapy are warranted for patients with *ERBB2*-mutated primary ILC. Key Words * Lobular * breast cancer * *ERBB2* * HER2 * mutation * prognosis * therapeutic biomarker * adjuvant ## Background Invasive lobular carcinoma of the breast (ILC) accounts for 10-15% of all breast cancer with an estimated 250,000 cases per year worldwide [1-4]. Nearly all cases of ILC derive from luminal cells that express estrogen receptor alpha but lack E-Cadherin (*CDH1*) expression or *ERBB2* amplification [5]. Clinical evidence suggests that despite favorable prognostic indicators, e.g. ER+ and/or progesterone receptor positive, low Ki67 proliferation index and HER2-status, patients with ILC have similar or worse long-term outcomes compared to those with invasive ductal carcinoma (IDC, otherwise known as invasive carcinoma of no special type) [5, 6]. Adjuvant treatment of ER+ ILC with letrozole, an aromatase inhibitor, may be superior to tamoxifen, and ILC cells demonstrate resistance to tamoxifen *in vitro* [7, 8]. However, patients with ILC are treated according to identical protocols as those with IDC [9, 10]. HER2-targeted therapy is indicated for patients whose tumors are HER2+ by immunohistochemistry (IHC) or, if IHC is equivocal, where *ERBB2* is amplified as detected by *in situ* hybridization (ISH) [11, 12]. Recent evidence indicates that in *ERBB2* non-amplified breast cancer, somatic mutation of *ERBB2* (*ERBB2*mut) and/or *ERBB3* (*ERBB3*mut) may provide an alternative mechanism for upregulation of HER2 activity that is therapeutically tractable using second generation HER2 tyrosine kinase inhibitors such as neratinib [13, 14]. The epidermal growth factor (EGF) family of receptor tyrosine kinases (HER1-4) are activated by ligand-dependent homo-/heterodimerisation and regulate cellular proliferation and tumor progression [15, 16]. In *ERBB2* amplified cells, the oncogenic effect of HER2 is mediated by heterodimerisation with HER3 in a ligand-independent manner [17]. Thus HER3 is necessary for HER2 oncogenic activity, and both HER2 and HER3 are therapeutic targets in *ERBB2* amplified breast cancer [18]. Mutations in *ERBB2* cluster in the tyrosine kinase and extra-cellular domains of HER2 and exert their oncogenic effects by activating tyrosine kinase activity or increasing HER2 dimerisation, respectively [13]. *In vitro* studies of HER2 activity in cell line and xenograft models identified 13 mutations (listed in the Methods) that enhanced proliferation and/or demonstrated growth inhibition with the irreversible HER2/EGFR tyrosine kinase inhibitor, neratinib [13]. Mutations in HER3, a critical binding partner for HER2, have been shown to promote ligand (EGF)-independent transformation of breast epithelial cells only in the presence of kinase-active HER2 [19]. This indicates that known oncogenic mutations in *ERBB3*, e.g. G284R and E928G, may also be therapeutically targetable *via* HER2 inhibition [19]. *ERBB2* mutations have previously been linked with worse prognosis in *CDH1*-altered ILC: a study of ILC cases in the TCGA dataset (*N*=169) found that *ERBB2* mutations (*N*=6) occurred exclusively in *CDH1-*altered tumors (*N*=100) [20]. However, prognostic data on the 6 *ERBB2*mut cases were limited to 2 patient events for both disease-free survival and OS analyses. The overall rate of *ERBB2*mut in this study of primary ILC was 3.6% (*N*=6 out of 169). A study of relapsed *CDH1*-mutated ILC found *ERBB2*mut in 18% of cases (*N*=4 of 22), suggesting further acquisition of *ERBB2* mutations in *CDH1*-altered ILC due to the selective pressure of treatment [21]. To demonstrate the potential clinical benefit of targeting low frequency somatic mutations, prognostic analyses using large clinical datasets are required. In the MA12 trial comprising 328 premenopausal patients with ER+ primary breast cancer of all histological subtypes, non-silent *ERBB2* mutations occured in 5.2% of patients (*N*=17) and were adversely prognostic of OS (*p*=0.0114) [22]. A study of 5,605 cases of relapsed breast cancer found *ERBB2*mut in 2.4% of cases (*N*=138), of which 20% (*N*=27 of 138) were in ILC tumors [23]. However, neither of these larger studies stratified clinical outcome by histological subtype. An ILC-specific study of 630 cases of primary ILC found *ERBB2* and *ERBB3* mutations in 5.1% and 3.6% of tumors respectively [24]. Comparison of cases of ER+, HER2-ILC from the same study (*N*=371) with cases of ER+, HER2-IDC from TCGA (*N*=338) indicated significant enrichment of both *ERBB2* and *ERBB3* mutations in ILC (4.3% and 3.5% in ILC *vs*. 1.5% and 0.6% in IDC) [24]. The study reported limited statistical evidence of a time-dependent effect of *ERBB2* mutational status associated with short-term breast cancer-specific survival. However, confirmation in datasets including patients with long-term follow-up and decoupling of activating from silent mutations is needed. We hypothesized that low-frequency somatic mutations in *ERBB2* and *ERBB3* are enriched in ER+, *ERBB2* non-amplified primary ILC cases and may have demonstrable prognostic effect. We tested these hypotheses by mining a combined dataset of the three largest primary breast cancer series with data on tumor *ERBB2* and *ERBB3* mutational status, gene expression, clinicopathological features and patient survival outcomes. Our overall goal was to determine the association between targetable *ERBB2/3* mutations and survival in ILC, and thereby provide evidence for a clinically actionable strategy to improve outcomes for patients with ILC. ### Patients and methods #### Patients and outcome measures Genomic and clinical outcome data associated with tumor samples from patients with primary breast cancer in TCGA 2015 (*N*=817), METABRIC 2012 and 2016 (*N*=2,509) and MSK-IMPACT 2018 (*N*=918) were accessed online *via* CBioportal [25-27]. From these datasets, ER+ and HER2-cases of stage I-III ILC and IDC with both clinical outcome and mutational data called from next-generation sequencing (NGS) analyses were selected (*N*=1,580). Cases of mixed or non-ILC/IDC histology, ER-negative/undetermined, HER2+/undetermined, carcinoma-*in-situ*, and stage 4/undetermined were excluded. For TCGA and MSK datasets, HER2 status was determined by IHC (positive/negative) or, where IHC was indeterminate, by ISH assessment of *ERBB2* amplification, in line with standard clinical practice [11]. For METABRIC cases, HER2 status was determined using the Affymetrix SNP6 copy number inference pipeline. The primary outcome measure, available in all datasets, was OS. Variables included *ERBB2* and *ERBB3* mutational status. *ERBB2*mut status was subcategorized as oncogenic or uncharacterized by cross-reference with existing data that identified mutations targetable by HER2-inhibition: G309A/E, S310F, L755S, del755-759, S760A, D769H, D769Y, V777L, P780ins, V842I, R896C [13]. Cases were denominated onc*ERBB2*mut if tumors harbored at least one oncogenic *ERBB2* mutation. Clinical and NGS mutation data were integrated with clinicopathological features including histological subtype, lymph node (LN) status and tumor size and grade. Normalized gene expression data were publicly available for METABRIC (Illumina HT12 microarray) and TCGA (RNA-seq) datasets. Baseline clinicopathological characteristics are summarized in Table 1. View this table: [Table 1:](http://medrxiv.org/content/early/2020/01/27/2020.01.24.20018622/T1) Table 1: Baseline clinicopathological characteristics of the combined cohort. *Cases of ILC/IDC histology, stage I-III, ER+ and HER2-status with clinical outcome and mutational data were selected *via* CBioportal. * *at last follow-up. ### Statistical analysis For analysis of binary somatic mutation status, combined cohort analysis was performed. Enrichment for cases with mutations in candidate genes (*ERBB2/3*mut) by histological subtype (ILC *vs*. IDC) was determined by χ2 test for association between categorical variables. For ILC and IDC separately, Kaplan-Meier (KM) survival curves stratified by mutation status were compared using Logrank and generalized Wilcoxon tests. Multivariate analysis of OS was performed using a Cox-regression model. Covariates included tumor grade, size (<20mm or ≥20mm) and LN status (positive or negative). Tumor grade was classified as low (grade 1-2) or high (grade 3). To derive a novel gene expression signature of HER2 activity that accounted for the effect of potentially targetable *ERBB2*mut in *ERBB2* non-amplified tumors, we applied a weighted average difference (WAD) method to gene expression data in cases from the METABRIC 2012 (*N*=1,980) and TCGA 2015 (*N*=817) cohorts [25, 28, 29]. Gene expression in *ERBB2*mut cases (*N*=38, selected by *ERBB2* non-amplified status and patient age>50) was compared with that in a similar number of *ERBB2* wild-type cases (*N*=79, selected by *ERBB2* non-amplified status, grade>1, stage>I, patient age>50). This was repeated for onc*ERBB2*mut cases (*N*=23) using the same comparator and selection criteria. To incorporate the effect of HER2 activity *via ERBB2* amplification, the overlap of differentially expressed genes DEGs shared by both comparisons (*ERBB2*mut and onc*ERBB2*mut *vs. ERBB2* wild-type) with DEGs from a further comparison of *ERBB2* amplified (*N*=247) *versus* non-amplified (*N*=1,733) cases in METABRIC was calculated. Finally, to incorporate the downstream phenotype (HER2 status), the overlap of this list with DEGs from a comparison of clinical HER2+ *versus* HER2-cases in TCGA was calculated. Multiple gene expression signatures of HER2 activity have been derived using cell line models and patient tumors [30-33]. We compared our novel gene signature with the HER2 activity signature established by Desmedt *et al* [31] with respect to its ability to detect potentially targetable *ERBB2*mut cases in our ILC/IDC dataset. This was achieved by multivariate regression modeling of response to neratinib for each gene signature using breast cancer cell line pharmacogenomic data from the BROAD Institute, accessed online *via* the CellMinerCDB portal [34, 35]. The Pearson coefficient for each significantly correlated signature gene was used to calculate normalized signature scores for each METABRIC case in the current dataset. To validate the prognostic effect of *ERBB2* in ILC, cases were then stratified by gene signature score (upper *vs*. lower quartiles) and the signatures compared by histological subtype using a Cox regression model of 10-year OS. ## Results ### Clinicopathological landscape of the combined cohort To detect prognostic effects of low frequency mutations and add to the existing body of literature on ILC-specific mutational drivers, a combined cohort of three public datasets was collated. We first evaluated the implications of combining cases of ILC and IDC from potentially disparate datasets. In the combined cohort, long-term follow-up of patients is dominated by the largest dataset (METABRIC, N=702: mean follow-up 133 months) whilst early events are enriched by the two smaller datasets (TCGA and MSK, N=878: mean follow-up 33 months). Compared to the smaller datasets with shorter follow-up (TCGA and MSK), patients in METABRIC were more likely to be over 50 years of age and have T1 tumours (<20mm diameter), but no significant difference in grade or LN status was found (see Table 2). As the principle skew in the combined dataset is towards longer follow-up in METABRIC cases, OS analysis was limited to 10 years – thus providing a clinically meaningful endpoint for all patients, irrespective of age. View this table: [Table 2:](http://medrxiv.org/content/early/2020/01/27/2020.01.24.20018622/T2) Table 2: Comparison of long *vs*. short follow-up cohorts: METABRIC (largest dataset, long follow-up) *vs*. TCGA and MSK (combined smaller datasets, short follow-up). Significant difference was found with respect to age and tumor size, but not tumor grade or LN status. Table excludes “unknown” cases for each variable. ### *ERBB2* mutations are enriched in primary ILC and cluster in the HER2 tyrosine kinase domain We next assessed the prevalence of *ERBB2/3*mut in our dataset of ILC and IDC (*N*=1,580). Overall prevalence of *ERBB2*mut was 2.2% (*N*=34). *ERBB2*mut was enriched in ILC, with prevalence of 5.7% (*N*=16) *vs*. 1.4% in IDC (*N*=18) (*p*<0.0001). In contrast, prevalence of *ERBB3*mut was lower (1.6% overall, *N*=26), and there was no enrichment of *ERBB3*mut in ILC *vs*. IDC (1.1%, *N*=3 *vs*. 1.8%, *N*=23; *p*=0.604). Due to the small number of cases in ILC, further analysis of *ERBB3*mut was not performed. In ILC, *ERBB2*mut clustered in the tyrosine kinase domain of HER2 (15 out of 16 cases; 93.8%). Of these, the majority have been characterized as oncogenic and potentially targetable using neratinib (onc*ERBB2*mut, 11 out of 15 kinase domain mutations in ILC; 73.3%). In IDC, all kinase domain mutations were known oncogenic (*N*=8) and non-characterized mutations were distributed evenly across protein domains (see Figure 1). The most frequently occurring *ERBB3*mut coded for E928G, which lies in the tyrosine kinase domain of HER3 (*N*=7 out of 26; 27%). ![Figure 1:](http://medrxiv.org/http://medrxiv.stage.highwire.org/content/medrxiv/early/2020/01/27/2020.01.24.20018622/F1.medium.gif) [Figure 1:](http://medrxiv.org/content/early/2020/01/27/2020.01.24.20018622/F1) Figure 1: Rate of *ERBB2* and *ERBB3* mutations and their spatial distribution on HER2/3 in ILC and IDC. *ERBB2*mut were found to be **(A)** enriched in ILC (yellow bar) *vs*. IDC (blue bar) and **(B)** clustered in the tyrosine kinase domain of HER2; **(C)** *ERBB3*mut occurred at lower frequency, with the high-frequency outlier in IDC coding for known oncogenic HER3 kinase domain alteration E928G (*N*=6). Y-axes show the number of cases harboring at least one *ERBB2/3* mutation at a specific amino acid (aa) of HER2/3, shown along the x-axes. Yellow-filled circles indicate onc*ERBB2*mut. Extracellular domains of HER2/3: Receptor L, Furin-like and Growth Factor Receptor IV; intracellular: tyrosine protein kinase. **p*<0.001; n/s = not significant. ### *ERBB2* mutation is an adverse prognostic indicator of survival for patients with primary ILC Using our meta-cohort of 1,580 cases, we next assessed *ERBB2*mut as a prognostic marker of OS in ILC (*N*=279) and IDC (*N*=1,301). Median duration of patient follow-up was 50 months (range 0 – 351 months). In patients with ILC, median OS was significantly shorter if tumors were *ERBB2*mut positive (inclusive of oncogenic and uncharacterized mutations) *versus ERBB2*mut negative (66 *vs*. 211 months, *p*=0.0001). In contrast, there was no significant difference in OS for *ERBB2*mut cases of IDC (159 *vs*. 166 months, *p*=0.733) (see Figure 2). ![Figure 2:](http://medrxiv.org/http://medrxiv.stage.highwire.org/content/medrxiv/early/2020/01/27/2020.01.24.20018622/F2.medium.gif) [Figure 2:](http://medrxiv.org/content/early/2020/01/27/2020.01.24.20018622/F2) Figure 2: OS by *ERBB2* mutational status in ILC (left) and IDC (right). Gray line indicates *ERBB2* wild-type cases; blue line indicates cases with at least one *ERBB2*mut. *ERBB3*mut status was not a significant prognostic indicator of OS in cases of IDC (HR=1.46, 95% CI 0.65 – 3.28; *p*=0.359). There were no events (patient deaths) at ten years of follow-up in the three cases of *ERBB3*mut ILC. ### Targetable *ERBB2* mutation status is an independent adverse prognostic marker of 10-year overall survival in ER+, *ERBB2* non-amplified ILC To test the effect of therapeutically actionable *ERBB2* mutations on a clinically relevant endpoint, we stratified 10-year OS by onc*ERBB2*mut status. In ER+, *ERBB2* non-amplified ILC, onc*ERBB2*mut status was prognostic of 10-year OS independently of LN status, tumor grade and size (HR 3.65, 95% CI 1.21 – 11.00; *p*=0.021, see Figure 3). ![Figure 3.](http://medrxiv.org/http://medrxiv.stage.highwire.org/content/medrxiv/early/2020/01/27/2020.01.24.20018622/F3.medium.gif) [Figure 3.](http://medrxiv.org/content/early/2020/01/27/2020.01.24.20018622/F3) Figure 3. Univariate and multivariate analyses of 10-year OS in *N*=279 cases of ER+, *ERBB2* non-amplified ILC. Gray square dot indicates hazard ratio (HR) in univariate analysis and the gray bar indicates the 95% confidence interval (CI). Significant prognostic variables in univariate analyses (where 95% CI does not span HR=1) are included in multivariate analysis, shown in blue. For each variable, cases with unknown values are excluded from the analysis. Unselected *ERBB2*mut status (including oncogenic and uncharacterized mutations) was also adversely prognostic of 10-year OS in univariate analysis (HR 3.66, 95% CI 1.54 – 8.73; *p*=0.003), but in contrast to characterized onc*ERBB2*mut this prognostic effect was not independent of LN stage, tumor grade and / or size. ### A novel gene signature of targetable HER2 activity using *ERBB2*mut status Since existing gene signatures of HER activity derive from HER2 status (by IHC and/or ISH) we generated a novel gene signature incorporating DEGs in *ERBB2*mut and onc*ERBB2*mut cases, as described in the Methods and outlined in Figure 4. A list of up and down-regulated *ERBB2* “mutant” DEGs (*N*=20) was generated by combining the overlap between DEGs for METABRIC *ERBB2* amplified *versus* non-amplified, *ERBB2* mutated *versus* wild-type (*ERBB2*mut and onc*ERBB2*mut separately) and TCGA HER2+ *versus* HER2-(Figure 4A, and see Additional file 1 for supplemental tables S1-3 with all gene lists). ![Figure 4:](http://medrxiv.org/http://medrxiv.stage.highwire.org/content/medrxiv/early/2020/01/27/2020.01.24.20018622/F4.medium.gif) [Figure 4:](http://medrxiv.org/content/early/2020/01/27/2020.01.24.20018622/F4) Figure 4: A novel gene signature of HER2 activity incorporating *ERBB2*mut, *ERBB2* amplification and clinical HER2 status. **(A)** Generation of a 20-gene signature of HER2 activity. The upper Venn diagrams show the overlap between the top 500 DEGs by WAD score for METABRIC amplified *vs*. non-amplified, *ERBB2* mutant *vs*. wild-type, and TCGA HER2+ *vs*. HER2-. Up-regulated DEGs are shaded blue and down-regulated DEGs yellow. *ERBB2*mut and onc*ERBB2*mut *vs*. wild-type are analysed separately, and the overlap combined in the lower Venn diagrams. **(B)** Comparison with an established 24-gene signature of HER pathway activation [31], using a gene signature (genesig) score derived from multivariate analysis of response to neratinib in breast cancer cell lines. Cases with *ERBB2*mut (gray lines) and onc*ERBB2*mut (blue lines) clustered in the upper quartile of normalized genesig scores for the novel signature but not the established signature. **(C)** 10-year OS analysis of cases in the current study stratified by histological subtype and novel genesig score (upper *vs*. lower quartiles) indicates that *ERBB2*mut-associated DEGs are prognostic in ILC but not IDC. GE = geneset enrichment. As shown in Figure 4B, *ERBB2*mut and onc*ERBB2*mut clustered in the upper quartile of the novel gene signature score in cases of ER+, HER2-ILC/IDC from METABRIC. In contrast, *ERBB2*mut and onc*ERBB2*mut cases did not cluster when scored by the established gene signature of HER2 pathway activation but were evenly distributed across the cohort. Both the novel and established HER2 pathway gene signatures were predictive of response to neratinib in multivariate analysis using a pharmacogenomic breast cancer cell line model (Pearson R=0.9, *p*<1e-14 for both; see Additional file 2 for the original plots from the CellMinerCDB online portal). The regression model used incorporates drug response and gene expression data from *N*=36 breast cancer cell lines [35, 36]. Our findings using this model verify that the novel signature derived from DEGs in *ERBB2*mut *versus* wild-type cases reflects oncogenic HER2 activity due to *ERBB2*mut. This suggests that a gene signature may have value in predicting the presence of *ERBB2*mut in HER2-breast cancer. Finally, Figure 4C demonstrates that the novel gene signature score (stratified into upper *vs*. lower quartiles) was adversely associated with 10-year OS (HR=2.3, 95% CI 1.04–5.05; *p*=0.040), thus providing mRNA level validation of *ERBB2*mut as an adverse prognostic marker. ## Discussion In this study, we mined clinical and NGS data from the largest clinical cohorts with data in the public domain to test whether there was an ILC-specific association of *ERBB2* mutations with OS. In our meta-cohort of ER+, *ERBB2* non-amplified cases of ILC and IDC, we found that *ERBB2* mutations are enriched ILC, cluster in the functional kinase domain of HER2, and robustly associate with adverse clinical outcomes – independently of known prognostic clinicopathological features including LN status and tumor grade. In contrast, there was no ILC-specific enrichment of *ERBB3*mut. Compared to the largest previous study with data on *ERBB2* mutational status in primary ILC (*N*=371) [24], the current study found higher frequency of *ERBB2*mut (5.7% *vs*. 4.3%) but lower frequency of *ERBB3*mut (1.1% *vs*. 3.5%). Hazard proportionality in the study by Desmedt *et al* 2016 suggested that *ERBB2* mutational status had a time-dependent effect associated with short-term risk of breast cancer relapse [24]. The current study adds data on the prognostic effect of specific *ERBB2* mutations that were previously shown in cell line, *in vitro* and xenograft studies to respond to existing clinical compounds such as trastuzumab or neratinib [13]. We demonstrated that the status of these therapeutically targetable *ERBB2* mutations is an independent adverse prognostic marker of overall survival in a large cohort of patients with ILC. Findings from the current study imply that targeted sequencing of ILC tumors for *ERBB2* mutations may be an actionable and viable strategy to improve outcomes for patients with primary ILC by providing a biomarker for HER-targeted therapy in the adjuvant setting. Based on worldwide ILC incidence of 250,000 [1-4], 10% rate of ER- and/or HER2-status in ILC, and 5% rate of potentially targetable *ERBB2*mut in ER+, *ERBB2* non-amplified ILC, we conservatively estimate that 11,250 additional ILC patients per year who have *ERBB2*-mutated tumors may benefit from HER2-targeting therapy. Existing clinical studies of HER2-targeting therapy are underway in advanced breast cancer. For example, neratinib (alone or in combination with fulvestrant) was tested in phase II study for patients with advanced, *ERBB2* non-amplified, *ERBB2*mut breast cancer (clinical trial registration number [NCT01670877](http://medrxiv.org/lookup/external-ref?link_type=CLINTRIALGOV&access_num=NCT01670877&atom=%2Fmedrxiv%2Fearly%2F2020%2F01%2F27%2F2020.01.24.20018622.atom)). Of the 16 patients who had evaluable tumor responses, *N*=2 had partial response by RECIST criteria (13%). Out of 14 patients with tumors harboring onc*ERBB2*mut, *N*=5 (36%) derived clinical benefit, defined as stable disease or partial/complete response by RECIST criteria [37]. Given that the patients in this trial had received multiple prior lines of treatment, it is conceivable that patients with *ERBB2* non-amplified primary tumors harboring onc*ERBB2*mut will also derive benefit from (neo)adjuvant HER2-targeted therapy such as neratinib. In the current cohort of primary ER+, HER2-ILC and IDC, *ERBB2* mutations in ILC tumors clustered almost exclusively in the active HER2 tyrosine kinase domain (*N*=15 of 16). Our retrospective clinical outcome data indicates an ILC-specific adverse survival association for patients with *ERBB2* non-amplified tumors harboring potentially targetable *ERBB2* mutations. This primary finding was validated by existing mRNA data *via* a novel gene signature of HER2 activity, which linked *ERBB2*mut with response to neratinib. Subject to independent external validation, a gene signature of HER2 activity could have clinical application as a biomarker of response to HER2-targeted therapy and may be more cost-effective than targeted NGS in the detection of *ERBB2*mut in primary ILC. Taken together, these new findings suggest that focusing early phase studies of HER2-inhibition (e.g. with neratinib) in patients with *ERBB2* non-amplified, *ERBB2*mut primary ILC may be an effective strategy to demonstrate feasibility and clinical benefit in the (neo)adjuvant setting. Limitations of the current study include cohort size, quality of clinical outcome data and bias in long-term follow-up towards cases from the largest dataset (METABRIC). Statistical power is limited by the small number of patient deaths in cases of ILC and IDC with *ERBB2*mut (*N*=7 and 8 respectively). However, the clear deviation of the KM plot of OS in ILC stratified by *ERBB2*mut (Figure 2) is reflected in a low type I error of one in 10,000 (*p*=0.0001). In contrast, there was no ILC-specific enrichment of *ERBB3*mut in our database. However, in IDC there was enrichment of the known oncogenic HER3 kinase domain mutation, E928G. To demonstrate the prognostic effect of mutations affecting HER3, larger datasets will be required. The effect of incomplete clinical outcome data on the primary endpoint (OS) is difficult to quantify but is accounted for statistically in Cox regression modelling by censoring patients at last follow-up. Patients in METABRIC presented at older age, which is associated with more indolent, less biologically aggressive breast cancer [38]. This implies an underestimation of the adverse impact of *ERBB2* mutations on 10-year survival outcomes in the current study. Since the METABRIC dataset did not have complete HER2 IHC or reverse-phase protein array data, it was not possible to associate activating *ERBB2*mut directly with downstream HER2 protein expression. Instead, we generated a gene signature incorporating clinical HER2 status to validate our findings. ## Conclusions Targetable kinase domain *ERBB2* mutations are enriched in primary ILC and their detection by targeted sequencing or validated gene signature surrogate may provide an actionable strategy to improve patient outcomes. Biomarker-led clinical trials of adjuvant HER2-targeted therapy to treat breast tumors with activating *ERBB2* mutations are warranted in HER2-primary ILC. ## Data Availability The datasets generated and/or analysed during the current study are available via the CBioportal online portal ([http://www.cbioportal.org](http://www.cbioportal.org)) and in the Genome-Phenome Archives EGAS00000000098 and phs000178. [https://ega-archive.org/studies/EGAS00000000098](https://ega-archive.org/studies/EGAS00000000098) [https://ega-archive.org/studies/phs000178](https://ega-archive.org/studies/phs000178) ## Declarations ### Ethics approval and consent to participate Ethics approvals and consent for each cohort (METABRIC, TCGA and MSK-IMPACT) was obtained in the original studies. ### Consent for publication Not applicable ### Availability of data and materials The datasets generated and/or analysed during the current study are available *via* the CBioportal online portal ([http://www.cbioportal.org](http://www.cbioportal.org)) and in the Genome-Phenome Archives EGAS00000000098 and phs000178 ([https://ega-archive.org/studies/EGAS00000000098](https://ega-archive.org/studies/EGAS00000000098) and [https://ega-archive.org/studies/phs000178](https://ega-archive.org/studies/phs000178)). ### Competing interests The authors declare that they have no competing interests. ### Funding The corresponding author (SJJ) is funded by the Wellcome Trust (Royal Academy of Medical Sciences grant number AAM 127669) and the National Institute for Health Research UK. Dr Oesterreich’s work on ILC is supported by the Breast Cancer Research Foundation and a Komen Scholar awards (SAC160073). ### Authors’ contributions SK, DMH, WZ, SO and SJJ conceived, designed and supervised the study. SK, MA, CM, KB, NPM and SJJ analysed and interpreted the data. SK, WZ, SO and SJJ prepared the manuscript. ARG, IOE and EAR had supervisory and administrative roles. SJ, TF and KS had consultative roles in study design. All authors were involved in manuscript preparation and reviewing for submission, gave their final approval and agreed to be accountable for all aspects of the work. ## Additional files **Additional file. 1** Supplemental tables: **(S1)** list of genes in the novel (*N*=20) and established (Desmedt *et al* [31], *N*=24) gene signatures of HER2 activity; **(S2 and S3)** top 500 DEGs by WAD score stratified into genes that are up-regulated (S2) and down-regulated (S3) for METABRIC *ERBB2* amplified *vs*. non-amplified (*ERBB2*amp), METABRIC *ERBB2*mut *vs*. wild-type (*ERBB2*mut) and onc*ERBB2*mut *vs*. wild-type (onc*ERBB2*mut), and TCGA HER2+ *vs*. HER2-(TCGA HER2+). **Additional file 2**. Supplemental figure showing the correlation between the observed response of breast cancer cell lines (*N*=36) to neratinib and the response predicted by expression of (A) novel (*N*=20) and (B) established (Desmedt *et al* [31], *N*=24) gene panels. Plots were generated using CellMinerCDB online portal and cell line data from the BROAD Institute [34, 35]. ## Acknowledgements Not applicable ## List of abbreviations CI : Confidence interval DEG : Differentially expressed gene EGF : Epidermal growth factor ER : Estrogen receptor alpha *ERBB2/3*mut : *ERBB2/3* mutant HER2/3 : Human epidermal growth factor receptor 2 / 3 HR : Hazard ratio IDC : Invasive ductal carcinoma of the breast IHC : Immunohistochemistry ILC : Invasive lobular carcinoma of the breast ISH : *In situ* hybridization KM : Kaplein-Meier LN : Lymph node NGS : Next generation sequencing onc*ERBB2*mut : Oncogenic *ERBB2* mutant OS : Overall survival RECIST : Response evaluation criteria in solid tumors WAD : Weighted average difference * Received January 24, 2020. * Revision received January 24, 2020. * Accepted January 27, 2020. * © 2020, Posted by Cold Spring Harbor Laboratory This pre-print is available under a Creative Commons License (Attribution-NoDerivs 4.0 International), CC BY-ND 4.0, as described at [http://creativecommons.org/licenses/by-nd/4.0/](http://creativecommons.org/licenses/by-nd/4.0/) ## References 1. 1.Martinez V, Azzopardi JG. Invasive lobular carcinoma of the breast: incidence and variants. Histopathology. 1979;3(6):467–88. [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1111/j.1365-2559.1979.tb03029.x&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=229072&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2020%2F01%2F27%2F2020.01.24.20018622.atom) [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=A1979HU94500003&link_type=ISI) 2. 2.Li CI, Anderson BO, Daling JR, Moe RE. Trends in incidence rates of invasive lobular and ductal breast carcinoma. JAMA. 2003;289(11):1421–4. [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1001/jama.289.11.1421&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=12636465&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2020%2F01%2F27%2F2020.01.24.20018622.atom) [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=000181616200034&link_type=ISI) 3. 3.Li CI, Daling JR. Changes in breast cancer incidence rates in the United States by histologic subtype and race/ethnicity, 1995 to 2004. Cancer Epidemiol Biomarkers Prev. 2007;16(12):2773–80. [Abstract/FREE Full Text](http://medrxiv.org/lookup/ijlink/YTozOntzOjQ6InBhdGgiO3M6MTQ6Ii9sb29rdXAvaWpsaW5rIjtzOjU6InF1ZXJ5IjthOjQ6e3M6ODoibGlua1R5cGUiO3M6NDoiQUJTVCI7czoxMToiam91cm5hbENvZGUiO3M6NDoiY2VicCI7czo1OiJyZXNpZCI7czoxMDoiMTYvMTIvMjc3MyI7czo0OiJhdG9tIjtzOjUwOiIvbWVkcnhpdi9lYXJseS8yMDIwLzAxLzI3LzIwMjAuMDEuMjQuMjAwMTg2MjIuYXRvbSI7fXM6ODoiZnJhZ21lbnQiO3M6MDoiIjt9) 4. 4.Ferlay J, Colombet M, Soerjomataram I, Mathers C, Parkin DM, Pineros M, et al. Estimating the global cancer incidence and mortality in 2018: GLOBOCAN sources and methods. Int J Cancer. 2019;144(8):1941–53. [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1002/ijc.31937&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=30350310&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2020%2F01%2F27%2F2020.01.24.20018622.atom) 5. 5.Pestalozzi BC, Zahrieh D, Mallon E, Gusterson BA, Price KN, Gelber RD, et al. Distinct clinical and prognostic features of infiltrating lobular carcinoma of the breast: combined results of 15 International Breast Cancer Study Group clinical trials. J Clin Oncol. 2008;26(18):3006–14. [Abstract/FREE Full Text](http://medrxiv.org/lookup/ijlink/YTozOntzOjQ6InBhdGgiO3M6MTQ6Ii9sb29rdXAvaWpsaW5rIjtzOjU6InF1ZXJ5IjthOjQ6e3M6ODoibGlua1R5cGUiO3M6NDoiQUJTVCI7czoxMToiam91cm5hbENvZGUiO3M6MzoiamNvIjtzOjU6InJlc2lkIjtzOjEwOiIyNi8xOC8zMDA2IjtzOjQ6ImF0b20iO3M6NTA6Ii9tZWRyeGl2L2Vhcmx5LzIwMjAvMDEvMjcvMjAyMC4wMS4yNC4yMDAxODYyMi5hdG9tIjt9czo4OiJmcmFnbWVudCI7czowOiIiO30=) 6. 6.Rakha EA, El-Sayed ME, Powe DG, Green AR, Habashy H, Grainge MJ, et al. Invasive lobular carcinoma of the breast: response to hormonal therapy and outcomes. Eur J Cancer. 2008;44(1):73–83. [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1016/j.ejca.2007.10.009&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=18035533&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2020%2F01%2F27%2F2020.01.24.20018622.atom) [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=000253038100018&link_type=ISI) 7. 7.Metzger Filho O, Giobbie-Hurder A, Mallon E, Gusterson B, Viale G, Winer EP, et al. Relative Effectiveness of Letrozole Compared With Tamoxifen for Patients With Lobular Carcinoma in the BIG 1-98 Trial. J Clin Oncol. 2015;33(25):2772–9. [Abstract/FREE Full Text](http://medrxiv.org/lookup/ijlink/YTozOntzOjQ6InBhdGgiO3M6MTQ6Ii9sb29rdXAvaWpsaW5rIjtzOjU6InF1ZXJ5IjthOjQ6e3M6ODoibGlua1R5cGUiO3M6NDoiQUJTVCI7czoxMToiam91cm5hbENvZGUiO3M6MzoiamNvIjtzOjU6InJlc2lkIjtzOjEwOiIzMy8yNS8yNzcyIjtzOjQ6ImF0b20iO3M6NTA6Ii9tZWRyeGl2L2Vhcmx5LzIwMjAvMDEvMjcvMjAyMC4wMS4yNC4yMDAxODYyMi5hdG9tIjt9czo4OiJmcmFnbWVudCI7czowOiIiO30=) 8. 8.Sikora MJ, Cooper KL, Bahreini A, Luthra S, Wang G, Chandran UR, et al. Invasive lobular carcinoma cell lines are characterized by unique estrogen-mediated gene expression patterns and altered tamoxifen response. Cancer Res. 2014;74(5):1463–74. [Abstract/FREE Full Text](http://medrxiv.org/lookup/ijlink/YTozOntzOjQ6InBhdGgiO3M6MTQ6Ii9sb29rdXAvaWpsaW5rIjtzOjU6InF1ZXJ5IjthOjQ6e3M6ODoibGlua1R5cGUiO3M6NDoiQUJTVCI7czoxMToiam91cm5hbENvZGUiO3M6NjoiY2FucmVzIjtzOjU6InJlc2lkIjtzOjk6Ijc0LzUvMTQ2MyI7czo0OiJhdG9tIjtzOjUwOiIvbWVkcnhpdi9lYXJseS8yMDIwLzAxLzI3LzIwMjAuMDEuMjQuMjAwMTg2MjIuYXRvbSI7fXM6ODoiZnJhZ21lbnQiO3M6MDoiIjt9) 9. 9.NICE. Early and locally advanced breast cancer: diagnosis and management. Online: National Institute for Health and Care Excellence; 2018 [NICE Pathway]. Available from: [https://www.nice.org.uk/guidance/ng101](https://www.nice.org.uk/guidance/ng101). 10. 10.Cardoso F, Kyriakides S, Ohno S, Penault-Llorca F, Poortmans P, Rubio IT, et al. Early breast cancer: ESMO Clinical Practice Guidelines for diagnosis, treatment and follow-updagger. Ann Oncol. 2019;30(8):1194–220. [PubMed](http://medrxiv.org/lookup/external-ref?access_num=http://www.n&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2020%2F01%2F27%2F2020.01.24.20018622.atom) 11. 11.Wolff AC, Hammond MEH, Allison KH, Harvey BE, Mangu PB, Bartlett JMS, et al. Human Epidermal Growth Factor Receptor 2 Testing in Breast Cancer: American Society of Clinical Oncology/College of American Pathologists Clinical Practice Guideline Focused Update. J Clin Oncol. 2018;36(20):2105–22. [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1200/JCO.2018.77.8738&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=http://www.n&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2020%2F01%2F27%2F2020.01.24.20018622.atom) 12. 12.Denduluri N, Chavez-MacGregor M, Telli ML, Eisen A, Graff SL, Hassett MJ, et al. Selection of Optimal Adjuvant Chemotherapy and Targeted Therapy for Early Breast Cancer: ASCO Clinical Practice Guideline Focused Update. J Clin Oncol. 2018;36(23):2433–43. 13. 13.Bose R, Kavuri SM, Searleman AC, Shen W, Shen D, Koboldt DC, et al. Activating HER2 mutations in HER2 gene amplification negative breast cancer. Cancer Discov. 2013;3(2):224–37. [Abstract/FREE Full Text](http://medrxiv.org/lookup/ijlink/YTozOntzOjQ6InBhdGgiO3M6MTQ6Ii9sb29rdXAvaWpsaW5rIjtzOjU6InF1ZXJ5IjthOjQ6e3M6ODoibGlua1R5cGUiO3M6NDoiQUJTVCI7czoxMToiam91cm5hbENvZGUiO3M6NzoiY2FuZGlzYyI7czo1OiJyZXNpZCI7czo3OiIzLzIvMjI0IjtzOjQ6ImF0b20iO3M6NTA6Ii9tZWRyeGl2L2Vhcmx5LzIwMjAvMDEvMjcvMjAyMC4wMS4yNC4yMDAxODYyMi5hdG9tIjt9czo4OiJmcmFnbWVudCI7czowOiIiO30=) 14. 14.Hanker A, Koch J, Ye D, Sliwoski G, Sheehan J, Kinch L, et al. Abstract PD3-05: Co-occurring gain-of-function mutations in HER2 and HER3 cooperate to enhance HER2/HER3 binding, HER- dependent signaling, and breast cancer growth. Cancer Research. 2019;79(4 Supplement):PD3-05-lPD3-. [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1158/1538-7445.SABCS18-PD3-05&link_type=DOI) 15. 15.Eroglu Z, Tagawa T, Somlo G. Human epidermal growth factor receptor family-targeted therapies in the treatment of HER2-overexpressing breast cancer. Oncologist. 2014;19(2):135–50. [Abstract/FREE Full Text](http://medrxiv.org/lookup/ijlink/YTozOntzOjQ6InBhdGgiO3M6MTQ6Ii9sb29rdXAvaWpsaW5rIjtzOjU6InF1ZXJ5IjthOjQ6e3M6ODoibGlua1R5cGUiO3M6NDoiQUJTVCI7czoxMToiam91cm5hbENvZGUiO3M6MTM6InRoZW9uY29sb2dpc3QiO3M6NToicmVzaWQiO3M6ODoiMTkvMi8xMzUiO3M6NDoiYXRvbSI7czo1MDoiL21lZHJ4aXYvZWFybHkvMjAyMC8wMS8yNy8yMDIwLjAxLjI0LjIwMDE4NjIyLmF0b20iO31zOjg6ImZyYWdtZW50IjtzOjA6IiI7fQ==) 16. 16.Hynes NE, MacDonald G. ErbB receptors and signaling pathways in cancer. Curr Opin Cell Biol. 2009;21(2):177–84. [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1016/j.ceb.2008.12.010&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=19208461&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2020%2F01%2F27%2F2020.01.24.20018622.atom) [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=000265773000006&link_type=ISI) 17. 17.Junttila TT, Akita RW, Parsons K, Fields C, Lewis Phillips GD, Friedman LS, et al. Ligand-independent HER2/HER3/PI3K complex is disrupted by trastuzumab and is effectively inhibited by the PI3K inhibitor GDC-0941. Cancer Cell. 2009;15(5):429–40. [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1016/j.ccr.2009.03.020&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=19411071&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2020%2F01%2F27%2F2020.01.24.20018622.atom) [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=000265903200010&link_type=ISI) 18. 18.Mishra R, Patel H, Alanazi S, Yuan L, Garrett JT. HER3 signaling and targeted therapy in cancer. Oncol Rev. 2018;12(1):355. [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.4081/oncol.2018.355&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=http://www.n&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2020%2F01%2F27%2F2020.01.24.20018622.atom) 19. 19.Jaiswal BS, Kljavin NM, Stawiski EW, Chan E, Parikh C, Durinck S, et al. Oncogenic ERBB3 mutations in human cancers. Cancer Cell. 2013;23(5):603–17. [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1016/j.ccr.2013.04.012&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=23680147&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2020%2F01%2F27%2F2020.01.24.20018622.atom) [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=000319103400010&link_type=ISI) 20. 20.Ping Z, Siegal GP, Harada S, Eltoum IE, Youssef M, Shen T, et al. ERBB2 mutation is associated with a worse prognosis in patients with CDH1 altered invasive lobular cancer of the breast. Oncotarget. 2016;7(49):80655–63. [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.18632/oncotarget.13019&link_type=DOI) 21. 21.Ross JS, Wang K, Sheehan CE, Boguniewicz AB, Otto G, Downing SR, et al. Relapsed classic E- cadherin (CDH1)-mutated invasive lobular breast cancer shows a high frequency of HER2 (ERBB2) gene mutations. Clin Cancer Res. 2013;19(10):2668–76. [Abstract/FREE Full Text](http://medrxiv.org/lookup/ijlink/YTozOntzOjQ6InBhdGgiO3M6MTQ6Ii9sb29rdXAvaWpsaW5rIjtzOjU6InF1ZXJ5IjthOjQ6e3M6ODoibGlua1R5cGUiO3M6NDoiQUJTVCI7czoxMToiam91cm5hbENvZGUiO3M6MTA6ImNsaW5jYW5yZXMiO3M6NToicmVzaWQiO3M6MTA6IjE5LzEwLzI2NjgiO3M6NDoiYXRvbSI7czo1MDoiL21lZHJ4aXYvZWFybHkvMjAyMC8wMS8yNy8yMDIwLjAxLjI0LjIwMDE4NjIyLmF0b20iO31zOjg6ImZyYWdtZW50IjtzOjA6IiI7fQ==) 22. 22.Griffith OL, Spies NC, Anurag M, Griffith M, Luo J, Tu D, et al. The prognostic effects of somatic mutations in ER-positive breast cancer. Nat Commun. 2018;9(1):3476. [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1038/s41467-018-05914-x&link_type=DOI) 23. 23.Ross JS, Gay LM, Wang K, Ali SM, Chumsri S, Elvin JA, et al. Nonamplification ERBB2 genomic alterations in 5605 cases of recurrent and metastatic breast cancer: An emerging opportunity for anti-HER2 targeted therapies. Cancer. 2016;122(17):2654–62. 24. 24.Desmedt C, Zoppoli G, Gundem G, Pruneri G, Larsimont D, Fornili M, et al. Genomic Characterization of Primary Invasive Lobular Breast Cancer. J Clin Oncol. 2016;34(16):1872–81. [Abstract/FREE Full Text](http://medrxiv.org/lookup/ijlink/YTozOntzOjQ6InBhdGgiO3M6MTQ6Ii9sb29rdXAvaWpsaW5rIjtzOjU6InF1ZXJ5IjthOjQ6e3M6ODoibGlua1R5cGUiO3M6NDoiQUJTVCI7czoxMToiam91cm5hbENvZGUiO3M6MzoiamNvIjtzOjU6InJlc2lkIjtzOjEwOiIzNC8xNi8xODcyIjtzOjQ6ImF0b20iO3M6NTA6Ii9tZWRyeGl2L2Vhcmx5LzIwMjAvMDEvMjcvMjAyMC4wMS4yNC4yMDAxODYyMi5hdG9tIjt9czo4OiJmcmFnbWVudCI7czowOiIiO30=) 25. 25.Ciriello G, Gatza ML, Beck AH, Wilkerson MD, Rhie SK, Pastore A, et al. Comprehensive Molecular Portraits of Invasive Lobular Breast Cancer. Cell. 2015;163(2):506–19. [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1016/j.cell.2015.09.033&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=26451490&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2020%2F01%2F27%2F2020.01.24.20018622.atom) 26. 26.Pereira B, Chin SF, Rueda OM, Vollan HK, Provenzano E, Bardwell HA, et al. The somatic mutation profiles of 2,433 breast cancers refines their genomic and transcriptomic landscapes. Nat Commun. 2016;7:11479. [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1038/ncomms11479&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=27161491&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2020%2F01%2F27%2F2020.01.24.20018622.atom) 27. 27.Razavi P, Chang MT, Xu G, Bandlamudi C, Ross DS, Vasan N, et al. The Genomic Landscape of Endocrine-Resistant Advanced Breast Cancers. Cancer Cell. 2018;34(3):427–38 e6. [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1016/j.ccell.2018.08.008&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=30205045&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2020%2F01%2F27%2F2020.01.24.20018622.atom) 28. 28.Kadota K, Nakai Y, Shimizu K. A weighted average difference method for detecting differentially expressed genes from microarray data. Algorithms Mol Biol. 2008;3:8. [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1186/1748-7188-3-8&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=18578891&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2020%2F01%2F27%2F2020.01.24.20018622.atom) 29. 29.Curtis C, Shah SP, Chin SF, Turashvili G, Rueda OM, Dunning MJ, et al. The genomic and transcriptomic architecture of 2,000 breast tumours reveals novel subgroups. Nature. 2012;486(7403):346–52. [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1038/nature10983&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=22522925&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2020%2F01%2F27%2F2020.01.24.20018622.atom) [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=000305466800033&link_type=ISI) 30. 30.Ferrari A, Vincent-Salomon A, Pivot X, Sertier AS, Thomas E, Tonon L, et al. A whole-genome sequence and transcriptome perspective on HER2-positive breast cancers. Nat Commun. 2016;7:12222. [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1038/ncomms12222&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=27406316&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2020%2F01%2F27%2F2020.01.24.20018622.atom) 31. 31.Desmedt C, Haibe-Kains B, Wirapati P, Buyse M, Larsimont D, Bontempi G, et al. Biological processes associated with breast cancer clinical outcome depend on the molecular subtypes. Clin Cancer Res. 2008;14(16):5158–65. [Abstract/FREE Full Text](http://medrxiv.org/lookup/ijlink/YTozOntzOjQ6InBhdGgiO3M6MTQ6Ii9sb29rdXAvaWpsaW5rIjtzOjU6InF1ZXJ5IjthOjQ6e3M6ODoibGlua1R5cGUiO3M6NDoiQUJTVCI7czoxMToiam91cm5hbENvZGUiO3M6MTA6ImNsaW5jYW5yZXMiO3M6NToicmVzaWQiO3M6MTA6IjE0LzE2LzUxNTgiO3M6NDoiYXRvbSI7czo1MDoiL21lZHJ4aXYvZWFybHkvMjAyMC8wMS8yNy8yMDIwLjAxLjI0LjIwMDE4NjIyLmF0b20iO31zOjg6ImZyYWdtZW50IjtzOjA6IiI7fQ==) 32. 32.Sareyeldin RM, Gupta I, Al-Hashimi I, Al-Thawadi HA, Al Farsi HF, Vranic S, et al. Gene Expression and miRNAs Profiling: Function and Regulation in Human Epidermal Growth Factor Receptor 2 (HER2)-Positive Breast Cancer. Cancers (Basel). 2019;11(5). 33. 33.Creighton CJ, Hilger AM, Murthy S, Rae JM, Chinnaiyan AM, El-Ashry D. Activation of mitogen-activated protein kinase in estrogen receptor alpha-positive breast cancer cells in vitro induces an in vivo molecular phenotype of estrogen receptor alpha-negative human breast tumors. Cancer Res. 2006;66(7):3903–11. [Abstract/FREE Full Text](http://medrxiv.org/lookup/ijlink/YTozOntzOjQ6InBhdGgiO3M6MTQ6Ii9sb29rdXAvaWpsaW5rIjtzOjU6InF1ZXJ5IjthOjQ6e3M6ODoibGlua1R5cGUiO3M6NDoiQUJTVCI7czoxMToiam91cm5hbENvZGUiO3M6NjoiY2FucmVzIjtzOjU6InJlc2lkIjtzOjk6IjY2LzcvMzkwMyI7czo0OiJhdG9tIjtzOjUwOiIvbWVkcnhpdi9lYXJseS8yMDIwLzAxLzI3LzIwMjAuMDEuMjQuMjAwMTg2MjIuYXRvbSI7fXM6ODoiZnJhZ21lbnQiO3M6MDoiIjt9) 34. 34.Rees MG, Seashore-Ludlow B, Cheah JH, Adams DJ, Price EV, Gill S, et al. Correlating chemical sensitivity and basal gene expression reveals mechanism of action. Nat Chem Biol. 2016;12(2):109–16. [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1038/nchembio.1986&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=26656090&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2020%2F01%2F27%2F2020.01.24.20018622.atom) 35. 35.Rajapakse VN, Luna A, Yamade M, Loman L, Varma S, Sunshine M, et al. CellMinerCDB for Integrative Cross-Database Genomics and Pharmacogenomics Analyses of Cancer Cell Lines. iScience. 2018;10:247–64. 36. 36.Basu A, Bodycombe NE, Cheah JH, Price EV, Liu K, Schaefer GI, et al. An interactive resource to identify cancer genetic and lineage dependencies targeted by small molecules. Cell. 2013;154(5):1151–61. [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1016/j.cell.2013.08.003&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=23993102&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2020%2F01%2F27%2F2020.01.24.20018622.atom) [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=000323767300026&link_type=ISI) 37. 37.Ma CX, Bose R, Gao F, Freedman RA, Pegram MD, Blackwell K, et al. Phase II trial of neratinib for HER2 mutated, non-amplified metastatic breast cancer (HER2(mut) MBC). J Clin Oncol. 2016;34(15). 38. 38.Syed BM, Green AR, Paish EC, Soria D, Garibaldi J, Morgan L, et al. Biology of primary breast cancer in older women treated by surgery: with correlation with long-term clinical outcome and comparison with their younger counterparts. Br J Cancer. 2013;108(5):1042–51. [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1038/bjc.2012.601&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=23462719&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2020%2F01%2F27%2F2020.01.24.20018622.atom)