Genomic landscape of ductal carcinoma in situ and association with progression.

PURPOSE: The detection rate of breast ductal carcinoma in situ (DCIS) has increased significantly, raising the concern that DCIS is overdiagnosed and overtreated. Therefore, there is an unmet clinical need to better predict the risk of progression among DCIS patients. Our hypothesis is that by combining molecular signatures with clinicopathologic features, we can elucidate the biology of breast cancer progression, and risk-stratify patients with DCIS. METHODS: Targeted exon sequencing with a custom panel of 223 genes/regions was performed for 125 DCIS cases. Among them, 60 were from cases having concurrent or subsequent invasive breast cancer (IBC) (DCIS + IBC group), and 65 from cases with no IBC development over a median follow-up of 13 years (DCIS-only group). Copy number alterations in chromosome 1q32, 8q24, and 11q13 were analyzed using fluorescence in situ hybridization (FISH). Multivariable logistic regression models were fit to the outcome of DCIS progression to IBC as functions of demographic and clinical features. RESULTS: We observed recurrent variants of known IBC-related mutations, and the most commonly mutated genes in DCIS were PIK3CA (34.4%) and TP53 (18.4%). There was an inverse association between PIK3CA kinase domain mutations and progression (Odds Ratio [OR] 10.2, p < 0.05). Copy number variations in 1q32 and 8q24 were associated with progression (OR 9.3 and 46, respectively; both p < 0.05). CONCLUSIONS: PIK3CA kinase domain mutations and the absence of copy number gains in DCIS are protective against progression to IBC. These results may guide efforts to distinguish low-risk from high-risk DCIS.

Phylogenetic analyses of melanoma reveal complex patterns of metastatic dissemination.

Melanoma is difficult to treat once it becomes metastatic. However, the precise ancestral relationship between primary tumors and their metastases is not well understood. We performed whole-exome sequencing of primary melanomas and multiple matched metastases from eight patients to elucidate their phylogenetic relationships. In six of eight patients, we found that genetically distinct cell populations in the primary tumor metastasized in parallel to different anatomic sites, rather than sequentially from one site to the next. In five of these six patients, the metastasizing cells had themselves arisen from a common parental subpopulation in the primary, indicating that the ability to establish metastases is a late-evolving trait. Interestingly, we discovered that individual metastases were sometimes founded by multiple cell populations of the primary that were genetically distinct. Such establishment of metastases by multiple tumor subpopulations could help explain why identical resistance variants are identified in different sites after initial response to systemic therapy. One primary tumor harbored two subclones with different oncogenic mutations in CTNNB1, which were both propagated to the same metastasis, raising the possibility that activation of wingless-type mouse mammary tumor virus integration site (WNT) signaling may be involved, as has been suggested by experimental models.

Genome co-amplification upregulates a mitotic gene network activity that predicts outcome and response to mitotic protein inhibitors in breast cancer.

BACKGROUND: High mitotic activity is associated with the genesis and progression of many cancers. Small molecule inhibitors of mitotic apparatus proteins are now being developed and evaluated clinically as anticancer agents. With clinical trials of several of these experimental compounds underway, it is important to understand the molecular mechanisms that determine high mitotic activity, identify tumor subtypes that carry molecular aberrations that confer high mitotic activity, and to develop molecular markers that distinguish which tumors will be most responsive to mitotic apparatus inhibitors. METHODS: We identified a coordinately regulated mitotic apparatus network by analyzing gene expression profiles for 53 malignant and non-malignant human breast cancer cell lines and two separate primary breast tumor datasets. We defined the mitotic network activity index (MNAI) as the sum of the transcriptional levels of the 54 coordinately regulated mitotic apparatus genes. The effect of those genes on cell growth was evaluated by small interfering RNA (siRNA). RESULTS: High MNAI was enriched in basal-like breast tumors and was associated with reduced survival duration and preferential sensitivity to inhibitors of the mitotic apparatus proteins, polo-like kinase, centromere associated protein E and aurora kinase designated GSK462364, GSK923295 and GSK1070916, respectively. Co-amplification of regions of chromosomes 8q24, 10p15-p12, 12p13, and 17q24-q25 was associated with the transcriptional upregulation of this network of 54 mitotic apparatus genes, and we identify transcription factors that localize to these regions and putatively regulate mitotic activity. Knockdown of the mitotic network by siRNA identified 22 genes that might be considered as additional therapeutic targets for this clinically relevant patient subgroup. CONCLUSIONS: We define a molecular signature which may guide therapeutic approaches for tumors with high mitotic network activity.

Subtype and pathway specific responses to anticancer compounds in breast cancer.

Breast cancers are comprised of molecularly distinct subtypes that may respond differently to pathway-targeted therapies now under development. Collections of breast cancer cell lines mirror many of the molecular subtypes and pathways found in tumors, suggesting that treatment of cell lines with candidate therapeutic compounds can guide identification of associations between molecular subtypes, pathways, and drug response. In a test of 77 therapeutic compounds, nearly all drugs showed differential responses across these cell lines, and approximately one third showed subtype-, pathway-, and/or genomic aberration-specific responses. These observations suggest mechanisms of response and resistance and may inform efforts to develop molecular assays that predict clinical response.

Loss-of-function mutations in Notch receptors in cutaneous and lung squamous cell carcinoma.

Squamous cell carcinomas (SCCs) are one of the most frequent forms of human malignancy, but, other than TP53 mutations, few causative somatic aberrations have been identified. We identified NOTCH1 or NOTCH2 mutations in ~75% of cutaneous SCCs and in a lesser fraction of lung SCCs, defining a spectrum for the most prevalent tumor suppressor specific to these epithelial malignancies. Notch receptors normally transduce signals in response to ligands on neighboring cells, regulating metazoan lineage selection and developmental patterning. Our findings therefore illustrate a central role for disruption of microenvironmental communication in cancer progression. NOTCH aberrations include frameshift and nonsense mutations leading to receptor truncations as well as point substitutions in key functional domains that abrogate signaling in cell-based assays. Oncogenic gain-of-function mutations in NOTCH1 commonly occur in human T-cell lymphoblastic leukemia/lymphoma and B-cell chronic lymphocytic leukemia. The bifunctional role of Notch in human cancer thus emphasizes the context dependency of signaling outcomes and suggests that targeted inhibition of the Notch pathway may induce squamous epithelial malignancies.

Cross-platform pathway-based analysis identifies markers of response to the PARP inhibitor olaparib.

Poly(ADP-ribose) polymerase (PARP) is an enzyme involved in DNA repair. PARP inhibitors can act as chemosensitizers, or operate on the principle of synthetic lethality when used as single agent. Clinical trials have shown drugs in this class to be promising for BRCA mutation carriers. We postulated that inability to demonstrate response in non-BRCA carriers in which BRCA is inactivated by other mechanisms or with deficiency in homologous recombination for DNA repair is due to lack of molecular markers that define a responding subpopulation. We identified candidate markers for this purpose for olaparib (AstraZeneca) by measuring inhibitory effects of nine concentrations of olaparib in 22 breast cancer cell lines and identifying features in transcriptional and genome copy number profiles that were significantly correlated with response. We emphasized in this discovery process genes involved in DNA repair. We found that the cell lines that were sensitive to olaparib had a significant lower copy number of BRCA1 compared to the resistant cell lines (p value 0.012). In addition, we discovered seven genes from DNA repair pathways whose transcriptional levels were associated with response. These included five genes (BRCA1, MRE11A, NBS1, TDG, and XPA) whose transcript levels were associated with resistance and two genes (CHEK2 and MK2) whose transcript levels were associated with sensitivity. We developed an algorithm to predict response using the seven-gene transcription levels and applied it to 1,846 invasive breast cancer samples from 8 U133A/plus 2 (Affymetrix) data sets and found that 8-21 % of patients would be predicted to be responsive to olaparib. A similar response frequency was predicted in 536 samples analyzed on an Agilent platform. Importantly, tumors predicted to respond were enriched in basal subtype tumors. Our studies support clinical evaluation of the utility of our seven-gene signature as a predictor of response to olaparib.

Xenografts faithfully recapitulate breast cancer-specific gene expression patterns of parent primary breast tumors.

Though xenografts are used extensively for drug development in breast cancer, how well xenografts reflect the breadth of primary breast tumor subtypes has not been well characterized. Moreover, few studies have compared the gene expression of xenograft tumors to the primary tumors from which they were derived. Here we investigate whether the ability of human breast tumors (n = 20) to create xenografts in immune-deficient mice is associated with breast cancer immunohistochemical (IHC) and intrinsic subtype. We also characterize how precisely the gene expression of xenografts reprises that of parent breast tumors, using hierarchical clustering and other correlation-based techniques applied to Agilent 44K gene expression data from 16 samples including four matched primary tumor-xenograft pairs. Of the breast tumors studied, 25 % (5/20) generated xenografts. Receptor and intrinsic subtype were significant predictors of xenograft success, with all (4/4) triple-negative (TN) tumors and no (0/12) HR+Her2- tumors forming xenografts (P = 0.0005). Tumor cell expression of ALDH1, a stem cell marker, trended toward successful engraftment (P = 0.14), though CDK5/6, a basal marker, did not. Though hierarchical clustering across the 500 most variable genes segregated human breast tumors from xenograft tumors, when clustering was performed over the PAM50 gene set the primary tumor-xenograft pairs clustered together, with all IHC subtypes clustered in distinct groups. Greater similarity between primary tumor-xenograft pairs relative to random pairings was confirmed by calculation of the within-pair between-pair scatter ratio (WPBPSR) distribution (P = 0.0269), though there was a shift in the xenografts toward more aggressive features including higher proliferation scores relative to the primary. Triple-negative breast tumors demonstrate superior ability to create xenografts compared to HR+ tumors, which may reflect higher proliferation or relatively stroma-independent growth of this subtype. Xenograft tumors' gene expression faithfully resembles that of their parent tumors, yet also demonstrates a shift toward more aggressive molecular features.

Comparative analyses of gene copy number and mRNA expression in glioblastoma multiforme tumors and xenografts.

Development of model systems that recapitulate the molecular heterogeneity observed among glioblastoma multiforme (GBM) tumors will expedite the testing of targeted molecular therapeutic strategies for GBM treatment. In this study, we profiled DNA copy number and mRNA expression in 21 independent GBM tumor lines maintained as subcutaneous xenografts (GBMX), and compared GBMX molecular signatures to those observed in GBM clinical specimens derived from the Cancer Genome Atlas (TCGA). The predominant copy number signature in both tumor groups was defined by chromosome-7 gain/chromosome-10 loss, a poor-prognosis genetic signature. We also observed, at frequencies similar to that detected in TCGA GBM tumors, genomic amplification and overexpression of known GBM oncogenes, such as EGFR, MDM2, CDK6, and MYCN, and novel genes, including NUP107, SLC35E3, MMP1, MMP13, and DDX1. The transcriptional signature of GBMX tumors, which was stable over multiple subcutaneous passages, was defined by overexpression of genes involved in M phase, DNA replication, and chromosome organization (MRC) and was highly similar to the poor-prognosis mitosis and cell-cycle module (MCM) in GBM. Assessment of gene expression in TCGA-derived GBMs revealed overexpression of MRC cancer genes AURKB, BIRC5, CCNB1, CCNB2, CDC2, CDK2, and FOXM1, which form a transcriptional network important for G2/M progression and/or checkpoint activation. Our study supports propagation of GBM tumors as subcutaneous xenografts as a useful approach for sustaining key molecular characteristics of patient tumors, and highlights therapeutic opportunities conferred by this GBMX tumor panel for testing targeted therapeutic strategies for GBM treatment.

A systems analysis of the chemosensitivity of breast cancer cells to the polyamine analogue PG-11047.

BACKGROUND: Polyamines regulate important cellular functions and polyamine dysregulation frequently occurs in cancer. The objective of this study was to use a systems approach to study the relative effects of PG-11047, a polyamine analogue, across breast cancer cells derived from different patients and to identify genetic markers associated with differential cytotoxicity. METHODS: A panel of 48 breast cell lines that mirror many transcriptional and genomic features present in primary human breast tumours were used to study the antiproliferative activity of PG-11047. Sensitive cell lines were further examined for cell cycle distribution and apoptotic response. Cell line responses, quantified by the GI50 (dose required for 50% relative growth inhibition) were correlated with the omic profiles of the cell lines to identify markers that predict response and cellular functions associated with drug sensitivity. RESULTS: The concentrations of PG-11047 needed to inhibit growth of members of the panel of breast cell lines varied over a wide range, with basal-like cell lines being inhibited at lower concentrations than the luminal cell lines. Sensitive cell lines showed a significant decrease in S phase fraction at doses that produced little apoptosis. Correlation of the GI50 values with the omic profiles of the cell lines identified genomic, transcriptional and proteomic variables associated with response. CONCLUSIONS: A 13-gene transcriptional marker set was developed as a predictor of response to PG-11047 that warrants clinical evaluation. Analyses of the pathways, networks and genes associated with response to PG-11047 suggest that response may be influenced by interferon signalling and differential inhibition of aspects of motility and epithelial to mesenchymal transition.

Multiple forms of atypical rearrangements generating supernumerary derivative chromosome 15.

BACKGROUND: Maternally-derived duplications that include the imprinted region on the proximal long arm of chromosome 15 underlie a complex neurobehavioral disorder characterized by cognitive impairment, seizures and a substantial risk for autism spectrum disorders1. The duplications most often take the form of a supernumerary pseudodicentric derivative chromosome 15 [der(15)] that has been called inverted duplication 15 or isodicentric 15 [idic(15)], although interstitial rearrangements also occur. Similar to the deletions found in most cases of Angelman and Prader Willi syndrome, the duplications appear to be mediated by unequal homologous recombination involving low copy repeats (LCR) that are found clustered in the region. Five recurrent breakpoints have been described in most cases of segmental aneuploidy of chromosome 15q11-q13 and previous studies have shown that most idic(15) chromosomes arise through BP3:BP3 or BP4:BP5 recombination events. RESULTS: Here we describe four duplication chromosomes that show evidence of atypical recombination events that involve regions outside the common breakpoints. Additionally, in one patient with a mosaic complex der(15), we examined homologous pairing of chromosome 15q11-q13 alleles by FISH in a region of frontal cortex, which identified mosaicism in this tissue and also demonstrated pairing of the signals from the der(15) and the normal homologues. CONCLUSION: Involvement of atypical BP in the generation of idic(15) chromosomes can lead to considerable structural heterogeneity.

Differentiation-state plasticity is a targetable resistance mechanism in basal-like breast cancer.

Intratumoral heterogeneity in cancers arises from genomic instability and epigenomic plasticity and is associated with resistance to cytotoxic and targeted therapies. We show here that cell-state heterogeneity, defined by differentiation-state marker expression, is high in triple-negative and basal-like breast cancer subtypes, and that drug tolerant persister (DTP) cell populations with altered marker expression emerge during treatment with a wide range of pathway-targeted therapeutic compounds. We show that MEK and PI3K/mTOR inhibitor-driven DTP states arise through distinct cell-state transitions rather than by Darwinian selection of preexisting subpopulations, and that these transitions involve dynamic remodeling of open chromatin architecture. Increased activity of many chromatin modifier enzymes, including BRD4, is observed in DTP cells. Co-treatment with the PI3K/mTOR inhibitor BEZ235 and the BET inhibitor JQ1 prevents changes to the open chromatin architecture, inhibits the acquisition of a DTP state, and results in robust cell death in vitro and xenograft regression in vivo.

Quantification of sensitivity and resistance of breast cancer cell lines to anti-cancer drugs using GR metrics.

Traditional means for scoring the effects of anti-cancer drugs on the growth and survival of cell lines is based on relative cell number in drug-treated and control samples and is seriously confounded by unequal division rates arising from natural biological variation and differences in culture conditions. This problem can be overcome by computing drug sensitivity on a per-division basis. The normalized growth rate inhibition (GR) approach yields per-division metrics for drug potency (GR50) and efficacy (GRmax) that are analogous to the more familiar IC50 and Emax values. In this work, we report GR-based, proliferation-corrected, drug sensitivity metrics for ~4,700 pairs of breast cancer cell lines and perturbagens. Such data are broadly useful in understanding the molecular basis of therapeutic response and resistance. Here, we use them to investigate the relationship between different measures of drug sensitivity and conclude that drug potency and efficacy exhibit high variation that is only weakly correlated. To facilitate further use of these data, computed GR curves and metrics can be browsed interactively at http://www.GRbrowser.org/.

Integrative molecular network analysis identifies emergent enzalutamide resistance mechanisms in prostate cancer.

Recent work demonstrates that castration-resistant prostate cancer (CRPC) tumors harbor countless genomic aberrations that control many hallmarks of cancer. While some specific mutations in CRPC may be actionable, many others are not. We hypothesized that genomic aberrations in cancer may operate in concert to promote drug resistance and tumor progression, and that organization of these genomic aberrations into therapeutically targetable pathways may improve our ability to treat CRPC. To identify the molecular underpinnings of enzalutamide-resistant CRPC, we performed transcriptional and copy number profiling studies using paired enzalutamide-sensitive and resistant LNCaP prostate cancer cell lines. Gene networks associated with enzalutamide resistance were revealed by performing an integrative genomic analysis with the PAthway Representation and Analysis by Direct Reference on Graphical Models (PARADIGM) tool. Amongst the pathways enriched in the enzalutamide-resistant cells were those associated with MEK, EGFR, RAS, and NFKB. Functional validation studies of 64 genes identified 10 candidate genes whose suppression led to greater effects on cell viability in enzalutamide-resistant cells as compared to sensitive parental cells. Examination of a patient cohort demonstrated that several of our functionally-validated gene hits are deregulated in metastatic CRPC tumor samples, suggesting that they may be clinically relevant therapeutic targets for patients with enzalutamide-resistant CRPC. Altogether, our approach demonstrates the potential of integrative genomic analyses to clarify determinants of drug resistance and rational co-targeting strategies to overcome resistance.

Correction: Decoupling of the PI3K Pathway via Mutation Necessitates Combinatorial Treatment in HER2+ Breast Cancer.

[This corrects the article DOI: 10.1371/journal.pone.0133219.].

Patient-specific factors influence somatic variation patterns in von Hippel-Lindau disease renal tumours.

Cancer development is presumed to be an evolutionary process that is influenced by genetic background and environment. In laboratory animals, genetics and environment are variables that can largely be held constant. In humans, it is possible to compare independent tumours that have developed in the same patient, effectively constraining genetic and environmental variation and leaving only stochastic processes. Patients affected with von Hippel-Lindau disease are at risk of developing multiple independent clear cell renal carcinomas. Here we perform whole-genome sequencing on 40 tumours from six von Hippel-Lindau patients. We confirm that the tumours are clonally independent, having distinct somatic single-nucleotide variants. Although tumours from the same patient show many differences, within-patient patterns are discernible. Single-nucleotide substitution type rates are significantly different between patients and show biases in trinucleotide mutation context. We also observe biases in chromosome copy number aberrations. These results show that genetic background and/or environment can influence the types of mutations that occur.

ARTIK-52 induces replication-dependent DNA damage and p53 activation exclusively in cells of prostate and breast cancer origin.

The realization, that the androgen receptor (AR) is essential for prostate cancer (PC) even after relapse following androgen deprivation therapy motivated the search for novel types of AR inhibitors. We proposed that targeting AR expression versus its function would work in cells having either wild type or mutant AR as well as be independent of androgen synthesis pathways. Previously, using a phenotypic screen in androgen-independent PC cells we identified a small molecule inhibitor of AR, ARTIK-52. Treatment with ARTIK-52 caused the loss of AR protein and death of AR-positive, but not AR-negative, PC cells. Here we present data that ARTIK-52 induces degradation of AR mRNA through a mechanism that we were unable to establish. However, we found that ARTIK-52 is toxic to breast cancer (BC) cells expressing AR, although they were not sensitive to AR knockdown, suggesting an AR-independent mechanism of toxicity. Using different approaches we detected that ARTIK-52 induces replication-dependent double strand DNA breaks exclusively in cancer cells of prostate and breast origin, while not causing DNA damage, or any toxicity, in normal cells, as well as in non-PC and non-BC tumor cells, independent of their proliferation status. This amazing specificity, combined with such a basic mechanism of toxicity, makes ARTIK-52 a potentially useful tool to discover novel attractive targets for the treatment of BC and PC. Thus, phenotypic screening allowed us to identify a compound, whose properties cannot be predicted based on existing knowledge and moreover, uncover a barely known link between AR and DNA damage response in PC and BC epithelial cells.

Cellular androgen content influences enzalutamide agonism of F877L mutant androgen receptor.

Prostate cancer is the most commonly diagnosed and second-most lethal cancer among men in the United States. The vast majority of prostate cancer deaths are due to castration-resistant prostate cancer (CRPC) - the lethal form of the disease that has progressed despite therapies that interfere with activation of androgen receptor (AR) signaling. One emergent resistance mechanism to medical castration is synthesis of intratumoral androgens that activate the AR. This insight led to the development of the AR antagonist enzalutamide. However, resistance to enzalutamide invariably develops, and disease progression is nearly universal. One mechanism of resistance to enzalutamide is an F877L mutation in the AR ligand-binding domain that can convert enzalutamide to an agonist of AR activity. However, mechanisms that contribute to the agonist switch had not been fully clarified, and there were no therapies to block AR F877L. Using cell line models of castration-resistant prostate cancer (CRPC), we determined that cellular androgen content influences enzalutamide agonism of mutant F877L AR. Further, enzalutamide treatment of AR F877L-expressing cell lines recapitulated the effects of androgen activation of F877L AR or wild-type AR. Because the BET bromodomain inhibitor JQ-1 was previously shown to block androgen activation of wild-type AR, we tested JQ-1 in AR F877L-expressing CRPC models. We determined that JQ-1 suppressed androgen or enzalutamide activation of mutant F877L AR and suppressed growth of mutant F877L AR CRPC tumors in vivo, demonstrating a new strategy to treat tumors harboring this mutation.

Exome Sequencing of Cell-Free DNA from Metastatic Cancer Patients Identifies Clinically Actionable Mutations Distinct from Primary Disease.

The identification of the molecular drivers of cancer by sequencing is the backbone of precision medicine and the basis of personalized therapy; however, biopsies of primary tumors provide only a snapshot of the evolution of the disease and may miss potential therapeutic targets, especially in the metastatic setting. A liquid biopsy, in the form of cell-free DNA (cfDNA) sequencing, has the potential to capture the inter- and intra-tumoral heterogeneity present in metastatic disease, and, through serial blood draws, track the evolution of the tumor genome. In order to determine the clinical utility of cfDNA sequencing we performed whole-exome sequencing on cfDNA and tumor DNA from two patients with metastatic disease; only minor modifications to our sequencing and analysis pipelines were required for sequencing and mutation calling of cfDNA. The first patient had metastatic sarcoma and 47 of 48 mutations present in the primary tumor were also found in the cell-free DNA. The second patient had metastatic breast cancer and sequencing identified an ESR1 mutation in the cfDNA and metastatic site, but not in the primary tumor. This likely explains tumor progression on Anastrozole. Significant heterogeneity between the primary and metastatic tumors, with cfDNA reflecting the metastases, suggested separation from the primary lesion early in tumor evolution. This is best illustrated by an activating PIK3CA mutation (H1047R) which was clonal in the primary tumor, but completely absent from either the metastasis or cfDNA. Here we show that cfDNA sequencing supplies clinically actionable information with minimal risks compared to metastatic biopsies. This study demonstrates the utility of whole-exome sequencing of cell-free DNA from patients with metastatic disease. cfDNA sequencing identified an ESR1 mutation, potentially explaining a patient's resistance to aromatase inhibition, and gave insight into how metastatic lesions differ from the primary tumor.

Upregulation of ER Signaling as an Adaptive Mechanism of Cell Survival in HER2-Positive Breast Tumors Treated with Anti-HER2 Therapy.

PURPOSE: To investigate the direct effect and therapeutic consequences of epidermal growth factor receptor 2 (HER2)-targeting therapy on expression of estrogen receptor (ER) and Bcl2 in preclinical models and clinical tumor samples. EXPERIMENTAL DESIGN: Archived xenograft tumors from two preclinical models (UACC812 and MCF7/HER2-18) treated with ER and HER2-targeting therapies and also HER2+ clinical breast cancer specimens collected in a lapatinib neoadjuvant trial (baseline and week 2 posttreatment) were used. Expression levels of ER and Bcl2 were evaluated by immunohistochemistry and Western blot analysis. The effects of Bcl2 and ER inhibition, by ABT-737 and fulvestrant, respectively, were tested in parental versus lapatinib-resistant UACC812 cells in vitro. RESULTS: Expression of ER and Bcl2 was significantly increased in xenograft tumors with acquired resistance to anti-HER2 therapy compared with untreated tumors in both preclinical models (UACC812: ER P = 0.0014; Bcl2 P < 0.001 and MCF7/HER2-18: ER P = 0.0007; Bcl2 P = 0.0306). In the neoadjuvant clinical study, lapatinib treatment for 2 weeks was associated with parallel upregulation of ER and Bcl2 (Spearman coefficient: 0.70; P = 0.0002). Importantly, 18% of tumors originally ER-negative (ER(-)) converted to ER(+) upon anti-HER2 therapy. In ER(-)/HER2(+) MCF7/HER2-18 xenografts, ER reexpression was primarily observed in tumors responding to potent combination of anti-HER2 drugs. Estrogen deprivation added to this anti-HER2 regimen significantly delayed tumor progression (P = 0.018). In the UACC812 cells, fulvestrant, but not ABT-737, was able to completely inhibit anti-HER2-resistant growth (P < 0.0001). CONCLUSIONS: HER2 inhibition can enhance or restore ER expression with parallel Bcl2 upregulation, representing an ER-dependent survival mechanism potentially leading to anti-HER2 resistance.