Sensitivity to targeted therapy differs between HER2-amplified breast cancer cells harboring kinase and helical domain mutations in PIK3CA.

BACKGROUND: HER2-amplified breast cancer is a clinically defined subtype of breast cancer for which there are multiple viable targeted therapies. Resistance to these targeted therapies is a common problem, but the mechanisms by which resistance occurs remain incompletely defined. One mechanism that has been proposed is through mutation of genes in the PI3-kinase pathway. Intracellular signaling from the HER2 pathway can occur through PI3-kinase, and mutations of the encoding gene PIK3CA are known to be oncogenic. Mutations in PIK3CA co-occur with HER2-amplification in ~ 20% of cases within the HER2-amplified subtype. METHODS: We generated isogenic knockin mutants of each PIK3CA hotspot mutation in HER2-amplified breast cancer cells using adeno-associated virus-mediated gene targeting. Isogenic clones were analyzed using a combinatorial drug screen to determine differential responses to HER2-targeted therapy. Western blot analysis and immunofluorescence uncovered unique intracellular signaling dynamics in cells resistant to HER2-targeted therapy. Subsequent combinatorial drug screens were used to explore neuregulin-1-mediated resistance to HER2-targeted therapy. Finally, results from in vitro experiments were extrapolated to publicly available datasets. RESULTS: Treatment with HER2-targeted therapy reveals that mutations in the kinase domain (H1047R) but not the helical domain (E545K) increase resistance to lapatinib. Mechanistically, sustained AKT signaling drives lapatinib resistance in cells with the kinase domain mutation, as demonstrated by staining for the intracellular product of PI3-kinase, PIP3. This resistance can be overcome by co-treatment with an inhibitor to the downstream kinase AKT. Additionally, knockout of the PIP3 phosphatase, PTEN, phenocopies this result. We also show that neuregulin-1, a ligand for HER-family receptors, confers resistance to cells harboring either hotspot mutation and modulates response to combinatorial therapy. Finally, we show clinical evidence that the hotspot mutations have distinct expression profiles related to therapeutic resistance through analysis of TCGA and METABRIC data cohorts. CONCLUSION: Our results demonstrate unique intracellular signaling differences depending on which mutation in PIK3CA the cell harbors. Only mutations in the kinase domain fully activate the PI3-kinase signaling pathway and maintain downstream signaling in the presence of HER2 inhibition. Moreover, we show there is potentially clinical importance in understanding both the PIK3CA mutational status and levels of neuregulin-1 expression in patients with HER2-amplified breast cancer treated with targeted therapy and that these problems warrant further pre-clinical and clinical testing.

GRB7 dependent proliferation of basal-like, HER-2 positive human breast cancer cell lines is mediated in part by HER-1 signaling.

GRB7 gene encodes a multi-domain signal transduction molecule and is part of the core of the HER-2 amplicon. GRB7 is commonly co-amplified and overexpressed with HER-2 in human breast cancer. This study addresses the role of GRB7 in HER-2 positive human breast cancers resistant to HER-2 targeted therapy. HCC1954, 21MT1, and JIMT1 are basal like HER-2 positive breast cancer cell lines based on expression profiling. These three cell lines are resistant to trastuzumab and lapatinib treatment. Knockdown of GRB7 protein expression with siRNA transfection as well as lentiviral vector mediated shRNA over-expression decreased the growth of HCC1954, 21MT1, and JIMT1 cells in vitro and the growth of tumor xenografts these cells formed in animal models. When assayed by ki-67 staining and TUNEL assay, the mechanism of reduced tumor xenograft growth appeared to be distinct. Reduced proliferation and increased apoptosis were seen in 21MT1 cells, while reduced proliferation was seen in HCC1954 cells and increased apoptosis in JIMT1 cells. Phospho-proteome profiling found HER-1 tyrosine phosphorylation was reduced with GRB7 knock down in JIMT1 cells. Immuno-blotting and immuno-precipitation experiments found HER-1 phosphorylation was reduced with GRB7 knock down in all three cell lines. HER-1 knock down via siRNA transient transfection as well as blocking HER-1 function with panitumumab decreased proliferation of all three cell lines in vitro. Our study finds that GRB7 has an essential growth promoting function which is mediated in part by HER-1 activation. The potential of HER-1 targeting in therapy resistant HER-2 positive breast cancer merits further study.

Proteomics advances for precision therapy in ovarian cancer.

Introduction: Due to the relatively low mutation rate and high frequency of copy number variation, finding actionable genetic drivers of high-grade serous carcinoma (HGSC) is a challenging task. Furthermore, emerging studies show that genetic alterations are frequently poorly represented at the protein level adding a layer of complexity. With improvements in large-scale proteomic technologies, proteomics studies have the potential to provide robust analysis of the pathways driving high HGSC behavior. Areas covered: This review summarizes recent large-scale proteomics findings across adequately sized ovarian cancer sample sets. Key words combined with 'ovarian cancer' including 'proteomics', 'proteogenomic', 'reverse-phase protein array', 'mass spectrometry', and 'adaptive response', were used to search PubMed. Expert opinion: Proteomics analysis of HGSC as well as their adaptive responses to therapy can uncover new therapeutic liabilities, which can reduce the emergence of drug resistance and potentially improve patient outcomes. There is a pressing need to better understand how the genomic and epigenomic heterogeneity intrinsic to ovarian cancer is reflected at the protein level and how this information could be used to improve patient outcomes.

MiR-486-5p Downregulation Marks an Early Event in Colorectal Carcinogenesis.

BACKGROUND: MicroRNAs are dysregulated in colorectal cancer and subsets correlated with advanced tumor stage and metastasis. Data are lacking on microRNA dysregulation from early to late-stage disease. OBJECTIVE: The purpose of this study was to identify a microRNA signature associated with the primary tumor and metastatic site in stage IV disease and to examine whether the signature is evident in earlier stages. DESIGN: A microRNA profile was generated and then explored in normal colon tissue (n = 5), early stage (stage I and II; n = 10), and late-stage (stage III and IV; n = 14) colorectal primary tumors via polymerase chain reaction to delineate molecular events that may promote colorectal carcinogenesis. SETTING: Genome-wide microRNA expression profiling was performed. PATIENTS: A total of 14 patient-matched stage IV primary colorectal cancer tumors and corresponding liver metastases were included. MAIN OUTCOME MEASURES: MicroRNA array technology was used to identify microRNA expression-predictive metastatic potential in the primary tumor. RESULTS: A distinct 9-member signature group of microRNAs was concurrent in stage IV primary colorectal cancer and their corresponding liver metastases, when compared with surrounding unaffected colon and liver tissue (microRNA-18b, microRNA-93, microRNA-182, microRNA-183, microRNA21, microRNA-486-5p, microRNA-500a, microRNA-552, and microRNA-941). Of the microRNA panel, only microRNA486-5p was differentially expressed in early stage colorectal cancer samples compared with normal tissue (p = 0.001) and additionally differentially expressed between late-stage colorectal cancer samples and normal tissue (p < 0.01). LIMITATIONS: Our microRNA profile was generated in a small subset of patients and will require validation in more samples. CONCLUSIONS: We identified a distinct microRNA signature in primary colon and matched metastatic disease. On additional investigation, 1 microRNA was differentially expressed in both early and late-stage cancer patient samples, and it may herald an early event in colorectal carcinogenesis. This study warrants additional investigation with a larger patient cohort to better understand the effect of microRNAs in carcinogenesis. See Video Abstract at http://links.lww.com/DCR/A723.

BRIT1 regulates early DNA damage response, chromosomal integrity, and cancer.

BRIT1, initially identified as an hTERT repressor, has additional functions at DNA damage checkpoints. Here, we demonstrate that BRIT1 formed nuclear foci minutes after irradiation. The foci of BRIT1 colocalized with 53BP1, MDC1, NBS1, ATM, RPA, and ATR. BRIT1 was required for activation of these elements, indicating that BRIT1 is a proximal factor in the DNA damage response pathway. Depletion of BRIT1 increased the accumulation of chromosomal aberrations. In addition, decreased levels of BRIT1 were detected in several types of human cancer, with BRIT1 expression being inversely correlated with genomic instability and metastasis. These results identify BRIT1 as a crucial DNA damage regulator in the ATM/ATR pathways and suggest that it functions as a tumor suppressor gene.

A collection of breast cancer cell lines for the study of functionally distinct cancer subtypes.

Recent studies suggest that thousands of genes may contribute to breast cancer pathophysiologies when deregulated by genomic or epigenomic events. Here, we describe a model "system" to appraise the functional contributions of these genes to breast cancer subsets. In general, the recurrent genomic and transcriptional characteristics of 51 breast cancer cell lines mirror those of 145 primary breast tumors, although some significant differences are documented. The cell lines that comprise the system also exhibit the substantial genomic, transcriptional, and biological heterogeneity found in primary tumors. We show, using Trastuzumab (Herceptin) monotherapy as an example, that the system can be used to identify molecular features that predict or indicate response to targeted therapies or other physiological perturbations.

Genomic and transcriptional aberrations linked to breast cancer pathophysiologies.

This study explores the roles of genome copy number abnormalities (CNAs) in breast cancer pathophysiology by identifying associations between recurrent CNAs, gene expression, and clinical outcome in a set of aggressively treated early-stage breast tumors. It shows that the recurrent CNAs differ between tumor subtypes defined by expression pattern and that stratification of patients according to outcome can be improved by measuring both expression and copy number, especially high-level amplification. Sixty-six genes deregulated by the high-level amplifications are potential therapeutic targets. Nine of these (FGFR1, IKBKB, ERBB2, PROCC, ADAM9, FNTA, ACACA, PNMT, and NR1D1) are considered druggable. Low-level CNAs appear to contribute to cancer progression by altering RNA and cellular metabolism.

FOXP3-positive regulatory T lymphocytes and epithelial FOXP3 expression in synchronous normal, ductal carcinoma in situ, and invasive cancer of the breast.

FOXP3-expressing T regulatory lymphocytes (Tregs) have been described as putative mediators of immune tolerance, and thus facilitators of tumor growth. When found in association with various malignancies, Tregs are generally markers of poor clinical outcome. However, it is unknown whether they are also associated with cancer progression. We evaluated quantitative FOXP3 expression in lymphocytes as well as in epithelial cells in a set of thirty-two breast tumors with synchronous normal epithelium, ductal carcinoma in situ (DCIS), and invasive ductal carcinoma (IDC) components. Tumors were stained for FOXP3 and CD3 expression and Tregs quantified by determining the ratio of colocalized FOXP3 and CD3 relative to 1) total CD3-expressing lymphocytes and 2) to FOXP3-expressing epithelial cells. The median proportion of FOXP3-expressing CD3 cells significantly increased with malignant progression from normal to DCIS to IDC components (0.005, 0.019 and 0.030, respectively; p ≤ 0.0001 for normal vs. IDC and p = 0.004 for DCIS vs. IDC). The median intensity of epithelial FOXP3 expression was also increased with invasive progression and most markedly augmented between normal and DCIS components (0.130 vs. 0.175, p ≤ 0.0001). Both Treg infiltration and epithelial FOXP3 expression were higher in grade 3 vs. grade 1 tumors (p = 0.014 for Tregs, p = 0.038 for epithelial FOXP3), but did not vary significantly with hormone receptor status, size of invasive tumor, lymph node status, or disease stage. Notably, Treg infiltration significantly correlated with epithelial up-regulation of FOXP3 expression (p = 0.013 for normal, p = 0.001 for IDC). These findings implicate both Treg infiltration and up-regulated epithelial FOXP3 expression in breast cancer progression.

In situ analyses of genome instability in breast cancer.

Transition through telomere crisis is thought to be a crucial event in the development of most breast carcinomas. Our goal in this study was to determine where this occurs in the context of histologically defined breast cancer progression. To this end, we assessed genome instability (using fluorescence in situ hybridization) and other features associated with telomere crisis in normal ductal epithelium, usual ductal hyperplasia, ductal carcinoma in situ and invasive cancer. We modeled this process in vitro by measuring these same features in human mammary epithelial cell cultures during ZNF217-mediated transition through telomere crisis and immortalization. Taken together, the data suggest that transition through telomere crisis and immortalization in breast cancer occurs during progression from usual ductal hyperplasia to ductal carcinoma in situ.

Robust biomarker discovery through multiplatform multiplex image analysis of breast cancer clinical cohorts.

Spatial profiling of tissues promises to elucidate tumor-microenvironment interactions and enable development of spatial biomarkers to predict patient response to immunotherapy and other therapeutics. However, spatial biomarker discovery is often carried out on a single patient cohort or imaging technology, limiting statistical power and increasing the likelihood of technical artifacts. In order to analyze multiple patient cohorts profiled on different platforms, we developed methods for comparative data analysis from three disparate multiplex imaging technologies: 1) cyclic immunofluorescence data we generated from 102 breast cancer patients with clinical follow-up, in addition to publicly available 2) imaging mass cytometry and 3) multiplex ion-beam imaging data. We demonstrate similar single-cell phenotyping results across breast cancer patient cohorts imaged with these three technologies and identify cellular abundance and proximity-based biomarkers with prognostic value across platforms. In multiple platforms, we identified lymphocyte infiltration as independently associated with longer survival in triple negative and high-proliferation breast tumors. Then, a comparison of nine spatial analysis methods revealed robust spatial biomarkers. In estrogen receptor-positive disease, quiescent stromal cells close to tumor were more abundant in good prognosis tumors while tumor neighborhoods of mixed fibroblast phenotypes were enriched in poor prognosis tumors. In triple-negative breast cancer (TNBC), macrophage proximity to tumor and B cell proximity to T cells were greater in good prognosis tumors, while tumor neighborhoods of vimentin-positive fibroblasts were enriched in poor prognosis tumors. We also tested previously published spatial biomarkers in our ensemble cohort, reproducing the positive prognostic value of isolated lymphocytes and lymphocyte occupancy and failing to reproduce the prognostic value of tumor-immune mixing score in TNBC. In conclusion, we demonstrate assembly of larger clinical cohorts from diverse platforms to aid in prognostic spatial biomarker identification and validation.

Micelle-Formulated Juglone Effectively Targets Pancreatic Cancer and Remodels the Tumor Microenvironment.

Pancreatic cancer remains a formidable challenge due to limited treatment options and its aggressive nature. In recent years, the naturally occurring anticancer compound juglone has emerged as a potential therapeutic candidate, showing promising results in inhibiting tumor growth and inducing cancer cell apoptosis. However, concerns over its toxicity have hampered juglone's clinical application. To address this issue, we have explored the use of polymeric micelles as a delivery system for juglone in pancreatic cancer treatment. These micelles, formulated using Poloxamer 407 and D-α-Tocopherol polyethylene glycol 1000 succinate, offer an innovative solution to enhance juglone's therapeutic potential while minimizing toxicity. In-vitro studies have demonstrated that micelle-formulated juglone (JM) effectively decreases proliferation and migration and increases apoptosis in pancreatic cancer cell lines. Importantly, in-vivo, JM exhibited no toxicity, allowing for increased dosing frequency compared to free drug administration. In mice, JM significantly reduced tumor growth in subcutaneous xenograft and orthotopic pancreatic cancer models. Beyond its direct antitumor effects, JM treatment also influenced the tumor microenvironment. In immunocompetent mice, JM increased immune cell infiltration and decreased stromal deposition and activation markers, suggesting an immunomodulatory role. To understand JM's mechanism of action, we conducted RNA sequencing and subsequent differential expression analysis on tumors that were treated with JM. The administration of JM treatment reduced the expression levels of the oncogenic protein MYC, thereby emphasizing its potential as a focused, therapeutic intervention. In conclusion, the polymeric micelles-mediated delivery of juglone holds excellent promise in pancreatic cancer therapy. This approach offers improved drug delivery, reduced toxicity, and enhanced therapeutic efficacy.

MYC Deregulation and PTEN Loss Model Tumor and Stromal Heterogeneity of Aggressive Triple-Negative Breast Cancer.

Triple-negative breast cancer (TNBC) patients have a poor prognosis and few treatment options. Mouse models of TNBC are important for development of new therapies, however, few mouse models represent the complexity of TNBC. Here, we develop a female TNBC murine model by mimicking two common TNBC mutations with high co-occurrence: amplification of the oncogene MYC and deletion of the tumor suppressor PTEN. This Myc;Ptenfl model develops heterogeneous triple-negative mammary tumors that display histological and molecular features commonly found in human TNBC. Our research involves deep molecular and spatial analyses on Myc;Ptenfl tumors including bulk and single-cell RNA-sequencing, and multiplex tissue-imaging. Through comparison with human TNBC, we demonstrate that this genetic mouse model develops mammary tumors with differential survival and therapeutic responses that closely resemble the inter- and intra-tumoral and microenvironmental heterogeneity of human TNBC, providing a pre-clinical tool for assessing the spectrum of patient TNBC biology and drug response.

Chemotherapy response and recurrence-free survival in neoadjuvant breast cancer depends on biomarker profiles: results from the I-SPY 1 TRIAL (CALGB 150007/150012; ACRIN 6657).

Neoadjuvant chemotherapy for breast cancer allows individual tumor response to be assessed depending on molecular subtype, and to judge the impact of response to therapy on recurrence-free survival (RFS). The multicenter I-SPY 1 TRIAL evaluated patients with ≥ 3 cm tumors by using early imaging and molecular signatures, with outcomes of pathologic complete response (pCR) and RFS. The current analysis was performed using data from patients who had molecular profiles and did not receive trastuzumab. The various molecular classifiers tested were highly correlated. Categorization of breast cancer by molecular signatures enhanced the ability of pCR to predict improvement in RFS compared to the population as a whole. In multivariate analysis, the molecular signatures that added to the ability of HR and HER2 receptors, clinical stage, and pCR in predicting RFS included 70-gene signature, wound healing signature, p53 mutation signature, and PAM50 risk of recurrence. The low risk signatures were associated with significantly better prognosis, and also identified additional patients with a good prognosis within the no pCR group, primarily in the hormone receptor positive, HER-2 negative subgroup. The I-SPY 1 population is enriched for tumors with a poor prognosis but is still heterogeneous in terms of rates of pCR and RFS. The ability of pCR to predict RFS is better by subset than it is for the whole group. Molecular markers improve prediction of RFS by identifying additional patients with excellent prognosis within the no pCR group.

Absence of telomerase and shortened telomeres have minimal effects on skin and pancreatic carcinogenesis elicited by viral oncogenes.

The telomere-stabilizing enzyme telomerase is induced in tumors and functionally associated with unlimited replicative potential. To further explore its necessity, transgenic mice expressing SV40 or HPV16 oncogenes, which elicit carcinomas in pancreas and skin, respectively, were rendered telomerase-deficient. Absence of telomerase had minimal impact on tumorigenesis, even in terc(-/-) generations (G5-7) exhibiting shortened telomeres and phenotypic abnormalities in multiple organs. Analyses of chromosomal aberrations were not indicative of telomere dysfunction or increased genomic instability in tumors. Quantitative image analysis of telomere repeat intensities comparing biopsies of skin hyperplasia, dysplasia, and carcinoma revealed that telomere numbers and relative lengths were maintained during progression, implicating a means for preserving telomere repeats and functionality in the absence of telomerase.

A framework for multiplex imaging optimization and reproducible analysis.

Multiplex imaging technologies are increasingly used for single-cell phenotyping and spatial characterization of tissues; however, transparent methods are needed for comparing the performance of platforms, protocols and analytical pipelines. We developed a python software, mplexable, for reproducible image processing and utilize Jupyter notebooks to share our optimization of signal removal, antibody specificity, background correction and batch normalization of the multiplex imaging with a focus on cyclic immunofluorescence (CyCIF). Our work both improves the CyCIF methodology and provides a framework for multiplexed image analytics that can be easily shared and reproduced.

An omic and multidimensional spatial atlas from serial biopsies of an evolving metastatic breast cancer.

Mechanisms of therapeutic resistance and vulnerability evolve in metastatic cancers as tumor cells and extrinsic microenvironmental influences change during treatment. To support the development of methods for identifying these mechanisms in individual people, here we present an omic and multidimensional spatial (OMS) atlas generated from four serial biopsies of an individual with metastatic breast cancer during 3.5 years of therapy. This resource links detailed, longitudinal clinical metadata that includes treatment times and doses, anatomic imaging, and blood-based response measurements to clinical and exploratory analyses, which includes comprehensive DNA, RNA, and protein profiles; images of multiplexed immunostaining; and 2- and 3-dimensional scanning electron micrographs. These data report aspects of heterogeneity and evolution of the cancer genome, signaling pathways, immune microenvironment, cellular composition and organization, and ultrastructure. We present illustrative examples of how integrative analyses of these data reveal potential mechanisms of response and resistance and suggest novel therapeutic vulnerabilities.

Toward reproducible, scalable, and robust data analysis across multiplex tissue imaging platforms.

The emergence of megascale single-cell multiplex tissue imaging (MTI) datasets necessitates reproducible, scalable, and robust tools for cell phenotyping and spatial analysis. We developed open-source, graphics processing unit (GPU)-accelerated tools for intensity normalization, phenotyping, and microenvironment characterization. We deploy the toolkit on a human breast cancer (BC) tissue microarray stained by cyclic immunofluorescence and present the first cross-validation of breast cancer cell phenotypes derived by using two different MTI platforms. Finally, we demonstrate an integrative phenotypic and spatial analysis revealing BC subtype-specific features.

Oligonucleotide conjugated antibody strategies for cyclic immunostaining.

A number of highly multiplexed immunostaining and imaging methods have advanced spatial proteomics of cancer for improved treatment strategies. While a variety of methods have been developed, the most widely used methods are limited by harmful signal removal techniques, difficulties with reagent production and antigen sensitivity. Multiplexed immunostaining employing oligonucleotide (oligos)-barcoded antibodies is an alternative approach that is growing in popularity. However, challenges remain in consistent conjugation of oligos to antibodies with maintained antigenicity as well as non-destructive, robust and cost-effective signal removal methods. Herein, a variety of oligo conjugation and signal removal methods were evaluated in the development of a robust oligo conjugated antibody cyclic immunofluorescence (Ab-oligo cyCIF) methodology. Both non- and site-specific conjugation strategies were assessed to label antibodies, where site-specific conjugation resulted in higher retained binding affinity and antigen-specific staining. A variety of fluorescence signal removal methods were also evaluated, where incorporation of a photocleavable link (PCL) resulted in full fluorescence signal removal with minimal tissue disruption. In summary, this work resulted in an optimized Ab-oligo cyCIF platform capable of generating high dimensional images to characterize the spatial proteomics of the hallmarks of cancer.

Genomic Alterations during the In Situ to Invasive Ductal Breast Carcinoma Transition Shaped by the Immune System.

The drivers of ductal carcinoma in situ (DCIS) to invasive ductal carcinoma (IDC) transition are poorly understood. Here, we conducted an integrated genomic, transcriptomic, and whole-slide image analysis to evaluate changes in copy-number profiles, mutational profiles, expression, neoantigen load, and topology in 6 cases of matched pure DCIS and recurrent IDC. We demonstrate through combined copy-number and mutational analysis that recurrent IDC can be genetically related to its pure DCIS despite long latency periods and therapeutic interventions. Immune "hot" and "cold" tumors can arise as early as DCIS and are subtype-specific. Topologic analysis showed a similar degree of pan-leukocyte-tumor mixing in both DCIS and IDC but differ when assessing specific immune subpopulations such as CD4 T cells and CD68 macrophages. Tumor-specific copy-number aberrations in MHC-I presentation machinery and losses in 3p, 4q, and 5p are associated with differences in immune signaling in estrogen receptor (ER)-negative IDC. Common oncogenic hotspot mutations in genes including TP53 and PIK3CA are predicted to be neoantigens yet are paradoxically conserved during the DCIS-to-IDC transition, and are associated with differences in immune signaling. We highlight both tumor and immune-specific changes in the transition of pure DCIS to IDC, including genetic changes in tumor cells that may have a role in modulating immune function and assist in immune escape, driving the transition to IDC. IMPLICATIONS: We demonstrate that the in situ to IDC evolutionary bottleneck is shaped by both tumor and immune cells.

The impact of tumor epithelial and microenvironmental heterogeneity on treatment responses in HER2+ breast cancer.

Despite the availability of multiple human epidermal growth factor receptor 2-targeted (HER2-targeted) treatments, therapeutic resistance in HER2+ breast cancer remains a clinical challenge. Intratumor heterogeneity for HER2 and resistance-conferring mutations in the PIK3CA gene (encoding PI3K catalytic subunit α) have been investigated in response and resistance to HER2-targeting agents, while the role of divergent cellular phenotypes and tumor epithelial-stromal cell interactions is less well understood. Here, we assessed the effect of intratumor cellular genetic heterogeneity for ERBB2 (encoding HER2) copy number and PIK3CA mutation on different types of neoadjuvant HER2-targeting therapies and clinical outcome in HER2+ breast cancer. We found that the frequency of cells lacking HER2 was a better predictor of response to HER2-targeted treatment than intratumor heterogeneity. We also compared the efficacy of different therapies in the same tumor using patient-derived xenograft models of heterogeneous HER2+ breast cancer and single-cell approaches. Stromal determinants were better predictors of response than tumor epithelial cells, and we identified alveolar epithelial and fibroblastic reticular cells as well as lymphatic vessel endothelial hyaluronan receptor 1-positive (Lyve1+) macrophages as putative drivers of therapeutic resistance. Our results demonstrate that both preexisting and acquired resistance to HER2-targeting agents involve multiple mechanisms including the tumor microenvironment. Furthermore, our data suggest that intratumor heterogeneity for HER2 should be incorporated into treatment design.