Functional delineation of the luminal epithelial microenvironment in breast using cell-based screening in combinatorial microenvironments.

Microenvironment signals are potent determinants of cell fate and arbiters of tissue homeostasis, however understanding how different microenvironment factors coordinately regulate cellular phenotype has been experimentally challenging. Here we used a high-throughput microenvironment microarray comprised of 2640 unique pairwise signals to identify factors that support proliferation and maintenance of primary human mammary luminal epithelial cells. Multiple microenvironment factors that modulated luminal cell number were identified, including: HGF, NRG1, BMP2, CXCL1, TGFB1, FGF2, PDGFB, RANKL, WNT3A, SPP1, HA, VTN, and OMD. All of these factors were previously shown to modulate luminal cell numbers in painstaking mouse genetics experiments, or were shown to have a role in breast cancer, demonstrating the relevance and power of our high-dimensional approach to dissect key microenvironmental signals. RNA-sequencing of primary epithelial and stromal cell lineages identified the cell types that express these signals and the cognate receptors in vivo. Cell-based functional studies confirmed which effects from microenvironment factors were reproducible and robust to individual variation. Hepatocyte growth factor (HGF) was the factor most robust to individual variation and drove expansion of luminal cells via cKit+ progenitor cells, which expressed abundant MET receptor. Luminal cells from women who are genetically high risk for breast cancer had significantly more MET receptor and may explain the characteristic expansion of the luminal lineage in those women. In ensemble, our approach provides proof of principle that microenvironment signals that control specific cellular states can be dissected with high-dimensional cell-based approaches.

Computational drug repositioning for the identification of new agents to sensitize drug-resistant breast tumors across treatments and receptor subtypes.

Introduction: Drug resistance is a major obstacle in cancer treatment and can involve a variety of different factors. Identifying effective therapies for drug resistant tumors is integral for improving patient outcomes. Methods: In this study, we applied a computational drug repositioning approach to identify potential agents to sensitize primary drug resistant breast cancers. We extracted drug resistance profiles from the I-SPY 2 TRIAL, a neoadjuvant trial for early stage breast cancer, by comparing gene expression profiles of responder and non-responder patients stratified into treatments within HR/HER2 receptor subtypes, yielding 17 treatment-subtype pairs. We then used a rank-based pattern-matching strategy to identify compounds in the Connectivity Map, a database of cell line derived drug perturbation profiles, that can reverse these signatures in a breast cancer cell line. We hypothesize that reversing these drug resistance signatures will sensitize tumors to treatment and prolong survival. Results: We found that few individual genes are shared among the drug resistance profiles of different agents. At the pathway level, however, we found enrichment of immune pathways in the responders in 8 treatments within the HR+HER2+, HR+HER2-, and HR-HER2- receptor subtypes. We also found enrichment of estrogen response pathways in the non-responders in 10 treatments primarily within the hormone receptor positive subtypes. Although most of our drug predictions are unique to treatment arms and receptor subtypes, our drug repositioning pipeline identified the estrogen receptor antagonist fulvestrant as a compound that can potentially reverse resistance across 13/17 of the treatments and receptor subtypes including HR+ and triple negative. While fulvestrant showed limited efficacy when tested in a panel of 5 paclitaxel resistant breast cancer cell lines, it did increase drug response in combination with paclitaxel in HCC-1937, a triple negative breast cancer cell line. Conclusion: We applied a computational drug repurposing approach to identify potential agents to sensitize drug resistant breast cancers in the I-SPY 2 TRIAL. We identified fulvestrant as a potential drug hit and showed that it increased response in a paclitaxel-resistant triple negative breast cancer cell line, HCC-1937, when treated in combination with paclitaxel.

Urinary Comprehensive Genomic Profiling Correlates Urothelial Carcinoma Mutations with Clinical Risk and Efficacy of Intervention.

The clinical standard of care for urothelial carcinoma (UC) relies on invasive procedures with suboptimal performance. To enhance UC treatment, we developed a urinary comprehensive genomic profiling (uCGP) test, UroAmplitude, that measures mutations from tumor DNA present in urine. In this study, we performed a blinded, prospective validation of technical sensitivity and positive predictive value (PPV) using reference standards, and found at 1% allele frequency, mutation detection performs at 97.4% sensitivity and 80.4% PPV. We then prospectively compared the mutation profiles of urine-extracted DNA to those of matched tumor tissue to validate clinical performance. Here, we found tumor single-nucleotide variants were observed in the urine with a median concordance of 91.7% and uCGP revealed distinct patterns of genomic lesions enriched in low- and high-grade disease. Finally, we retrospectively explored longitudinal case studies to quantify residual disease following bladder-sparing treatments, and found uCGP detected residual disease in patients receiving bladder-sparing treatment and predicted recurrence and disease progression. These findings demonstrate the potential of the UroAmplitude platform to reliably identify and track mutations associated with UC at each stage of disease: diagnosis, treatment, and surveillance. Multiple case studies demonstrate utility for patient risk classification to guide both surgical and therapeutic interventions.

A multiplex implantable microdevice assay identifies synergistic combinations of cancer immunotherapies and conventional drugs.

Systematically identifying synergistic combinations of targeted agents and immunotherapies for cancer treatments remains difficult. In this study, we integrated high-throughput and high-content techniques-an implantable microdevice to administer multiple drugs into different sites in tumors at nanodoses and multiplexed imaging of tumor microenvironmental states-to investigate the tumor cell and immunological response signatures to different treatment regimens. Using a mouse model of breast cancer, we identified effective combinations from among numerous agents within days. In vivo studies in three immunocompetent mammary carcinoma models demonstrated that the predicted combinations synergistically increased therapeutic efficacy. We identified at least five promising treatment strategies, of which the panobinostat, venetoclax and anti-CD40 triple therapy was the most effective in inducing complete tumor remission across models. Successful drug combinations increased spatial association of cancer stem cells with dendritic cells during immunogenic cell death, suggesting this as an important mechanism of action in long-term breast cancer control.

INHBA is a mediator of aggressive tumor behavior in HER2+ basal breast cancer.

BACKGROUND: Resistance to HER2-targeted therapeutics remains a significant clinical problem in HER2+ breast cancer patients with advanced disease. This may be particularly true for HER2+ patients with basal subtype disease, as recent evidence suggests they receive limited benefit from standard of care HER2-targeted therapies. Identification of drivers of resistance and aggressive disease that can be targeted clinically has the potential to impact patient outcomes. METHODS: We performed siRNA knockdown screens of genes differentially expressed between lapatinib-responsive and -resistant HER2+ breast cancer cells, which corresponded largely to luminal versus basal subtypes. We then validated hits in 2-d and 3-d cell culture systems. RESULTS: Knockdown of one of the genes, INHBA, significantly slowed growth and increased sensitivity to lapatinib in multiple basal HER2+ cell lines in both 2-d and 3-d cultures, but had no effect in luminal HER2+ cells. Loss of INHBA altered metabolism, eliciting a shift from glycolytic to oxidative phosphorylative metabolism, which was also associated with a decrease in tumor invasiveness. Analysis of breast cancer datasets showed that patients with HER2+ breast cancer and high levels of INHBA expression had worse outcomes than patients with low levels of INHBA expression. CONCLUSIONS: Our data suggest that INHBA is associated with aggressiveness of the basal subtype of HER2+ tumors, resulting in poor response to HER2-targeted therapy and an invasive phenotype. We hypothesize that targeting this pathway could be an effective therapeutic strategy to reduce invasiveness of tumor cells and to improve therapeutic response.

Ferrous iron-activatable drug conjugate achieves potent MAPK blockade in KRAS-driven tumors.

KRAS mutations drive a quarter of cancer mortality, and most are undruggable. Several inhibitors of the MAPK pathway are FDA approved but poorly tolerated at the doses needed to adequately extinguish RAS/RAF/MAPK signaling in the tumor cell. We found that oncogenic KRAS signaling induced ferrous iron (Fe2+) accumulation early in and throughout mutant KRAS-mediated transformation. We converted an FDA-approved MEK inhibitor into a ferrous iron-activatable drug conjugate (FeADC) and achieved potent MAPK blockade in tumor cells while sparing normal tissues. This innovation allowed sustainable, effective treatment of tumor-bearing animals, with tumor-selective drug activation, producing superior systemic tolerability. Ferrous iron accumulation is an exploitable feature of KRAS transformation, and FeADCs hold promise for improving the treatment of KRAS-driven solid tumors.

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.

Cellular senescence in cancer: from mechanisms to detection.

Senescence refers to a cellular state featuring a stable cell-cycle arrest triggered in response to stress. This response also involves other distinct morphological and intracellular changes including alterations in gene expression and epigenetic modifications, elevated macromolecular damage, metabolism deregulation and a complex pro-inflammatory secretory phenotype. The initial demonstration of oncogene-induced senescence in vitro established senescence as an important tumour-suppressive mechanism, in addition to apoptosis. Senescence not only halts the proliferation of premalignant cells but also facilitates the clearance of affected cells through immunosurveillance. Failure to clear senescent cells owing to deficient immunosurveillance may, however, lead to a state of chronic inflammation that nurtures a pro-tumorigenic microenvironment favouring cancer initiation, migration and metastasis. In addition, senescence is a response to post-therapy genotoxic stress. Therefore, tracking the emergence of senescent cells becomes pivotal to detect potential pro-tumorigenic events. Current protocols for the in vivo detection of senescence require the analysis of fixed or deep-frozen tissues, despite a significant clinical need for real-time bioimaging methods. Accuracy and efficiency of senescence detection are further hampered by a lack of universal and more specific senescence biomarkers. Recently, in an attempt to overcome these hurdles, an assortment of detection tools has been developed. These strategies all have significant potential for clinical utilisation and include flow cytometry combined with histo- or cytochemical approaches, nanoparticle-based targeted delivery of imaging contrast agents, OFF-ON fluorescent senoprobes, positron emission tomography senoprobes and analysis of circulating SASP factors, extracellular vesicles and cell-free nucleic acids isolated from plasma. Here, we highlight the occurrence of senescence in neoplasia and advanced tumours, assess the impact of senescence on tumorigenesis and discuss how the ongoing development of senescence detection tools might improve early detection of multiple cancers and response to therapy in the near future.

Enzalutamide response in a panel of prostate cancer cell lines reveals a role for glucocorticoid receptor in enzalutamide resistant disease.

Representative in vitro model systems that accurately model response to therapy and allow the identification of new targets are important for improving our treatment of prostate cancer. Here we describe molecular characterization and drug testing in a panel of 20 prostate cancer cell lines. The cell lines cluster into distinct subsets based on RNA expression, which is largely driven by functional Androgen Receptor (AR) expression. KLK3, the AR-responsive gene that encodes prostate specific antigen, shows the greatest variability in expression across the cell line panel. Other common prostate cancer associated genes such as TMPRSS2 and ERG show similar expression patterns. Copy number analysis demonstrates that many of the most commonly gained (including regions containing TERC and MYC) and lost regions (including regions containing TP53 and PTEN) that were identified in patient samples by the TCGA are mirrored in the prostate cancer cell lines. Assessment of response to the anti-androgen enzalutamide shows a distinct separation of responders and non-responders, predominantly related to status of wild-type AR. Surprisingly, several AR-null lines responded to enzalutamide. These AR-null, enzalutamide-responsive cells were characterized by high levels of expression of glucocorticoid receptor (GR) encoded by NR3C1. Treatment of these cells with the anti-GR agent mifepristone showed that they were more sensitive to this drug than enzalutamide, as were several of the enzalutamide non-responsive lines. This is consistent with several recent reports that suggest that GR expression is an alternative signaling mechanism that can bypass AR blockade. This study reinforces the utility of large cell line panels for the study of cancer and identifies several cell lines that represent ideal models to study AR-null cells that have upregulated GR to sustain growth.

Transcriptional signatures in histologic structures within glioblastoma tumors may predict personalized drug sensitivity and survival.

BACKGROUND: Glioblastoma is a rapidly fatal brain cancer that exhibits extensive intra- and intertumoral heterogeneity. Improving survival will require the development of personalized treatment strategies that can stratify tumors into subtypes that differ in therapeutic vulnerability and outcomes. Glioblastoma stratification has been hampered by intratumoral heterogeneity, limiting our ability to compare tumors in a consistent manner. Here, we develop methods that mitigate the impact of intratumoral heterogeneity on transcriptomic-based patient stratification. METHODS: We accessed open-source transcriptional profiles of histological structures from 34 human glioblastomas from the Ivy Glioblastoma Atlas Project. Principal component and correlation network analyses were performed to assess sample inter-relationships. Gene set enrichment analysis was used to identify enriched biological processes and classify glioblastoma subtype. For survival models, Cox proportional hazards regression was utilized. Transcriptional profiles from 156 human glioblastomas were accessed from The Cancer Genome Atlas to externally validate the survival model. RESULTS: We showed that intratumoral histologic architecture influences tumor classification when assessing established subtyping and prognostic gene signatures, and that indiscriminate sampling can produce misleading results. We identified the cellular tumor as a glioblastoma structure that can be targeted for transcriptional analysis to more accurately stratify patients by subtype and prognosis. Based on expression from cellular tumor, we created an improved risk stratification gene signature. CONCLUSIONS: Our results highlight that biomarker performance for diagnostics, prognostics, and prediction of therapeutic response can be improved by analyzing transcriptional profiles in pure cellular tumor, which is a critical step toward developing personalized treatment for glioblastoma.

A Multi-center Study on the Reproducibility of Drug-Response Assays in Mammalian Cell Lines.

Evidence that some high-impact biomedical results cannot be repeated has stimulated interest in practices that generate findable, accessible, interoperable, and reusable (FAIR) data. Multiple papers have identified specific examples of irreproducibility, but practical ways to make data more reproducible have not been widely studied. Here, five research centers in the NIH LINCS Program Consortium investigate the reproducibility of a prototypical perturbational assay: quantifying the responsiveness of cultured cells to anti-cancer drugs. Such assays are important for drug development, studying cellular networks, and patient stratification. While many experimental and computational factors impact intra- and inter-center reproducibility, the factors most difficult to identify and control are those with a strong dependency on biological context. These factors often vary in magnitude with the drug being analyzed and with growth conditions. We provide ways to identify such context-sensitive factors, thereby improving both the theory and practice of reproducible cell-based assays.

Using Microarrays to Interrogate Microenvironmental Impact on Cellular Phenotypes in Cancer.

Understanding the impact of the microenvironment on the phenotype of cells is a difficult problem due to the complex mixture of both soluble growth factors and matrix-associated proteins in the microenvironment in vivo. Furthermore, readily available reagents for the modeling of microenvironments in vitro typically utilize complex mixtures of proteins that are incompletely defined and suffer from batch to batch variability. The microenvironment microarray (MEMA) platform allows for the assessment of thousands of simple combinations of microenvironment proteins for their impact on cellular phenotypes in a single assay. The MEMAs are prepared in well plates, which allows the addition of individual ligands to separate wells containing arrayed extracellular matrix (ECM) proteins. The combination of the soluble ligand with each printed ECM forms a unique combination. A typical MEMA assay contains greater than 2,500 unique combinatorial microenvironments that cells are exposed to in a single assay. As a test case, the breast cancer cell line MCF7 was plated on the MEMA platform. Analysis of this assay identified factors that both enhance and inhibit the growth and proliferation of these cells. The MEMA platform is highly flexible and can be extended for use with other biological questions beyond cancer research.

Maintenance of MYC expression promotes de novo resistance to BET bromodomain inhibition in castration-resistant prostate cancer.

The BET bromodomain protein BRD4 is a chromatin reader that regulates transcription, including in cancer. In prostate cancer, specifically, the anti-tumor activity of BET bromodomain inhibition has been principally linked to suppression of androgen receptor (AR) function. MYC is a well-described BRD4 target gene in multiple cancer types, and prior work demonstrates that MYC plays an important role in promoting prostate cancer cell survival. Importantly, several BET bromodomain clinical trials are ongoing, including in prostate cancer. However, there is limited information about pharmacodynamic markers of response or mediators of de novo resistance. Using a panel of prostate cancer cell lines, we demonstrated that MYC suppression-rather than AR suppression-is a key determinant of BET bromodomain inhibitor sensitivity. Importantly, we determined that BRD4 was dispensable for MYC expression in the most resistant cell lines and that MYC RNAi + BET bromodomain inhibition led to additive anti-tumor activity in the most resistant cell lines. Our findings demonstrate that MYC suppression is an important pharmacodynamic marker of BET bromodomain inhibitor response and suggest that targeting MYC may be a promising therapeutic strategy to overcome de novo BET bromodomain inhibitor resistance in prostate cancer.

Modeling Tumor Phenotypes In Vitro with Three-Dimensional Bioprinting.

The tumor microenvironment plays a critical role in tumor growth, progression, and therapeutic resistance, but interrogating the role of specific tumor-stromal interactions on tumorigenic phenotypes is challenging within in vivo tissues. Here, we tested whether three-dimensional (3D) bioprinting could improve in vitro models by incorporating multiple cell types into scaffold-free tumor tissues with defined architecture. We generated tumor tissues from distinct subtypes of breast or pancreatic cancer in relevant microenvironments and demonstrate that this technique can model patient-specific tumors by using primary patient tissue. We assess intrinsic, extrinsic, and spatial tumorigenic phenotypes in bioprinted tissues and find that cellular proliferation, extracellular matrix deposition, and cellular migration are altered in response to extrinsic signals or therapies. Together, this work demonstrates that multi-cell-type bioprinted tissues can recapitulate aspects of in vivo neoplastic tissues and provide a manipulable system for the interrogation of multiple tumorigenic endpoints in the context of distinct tumor microenvironments.

Kinome rewiring reveals AURKA limits PI3K-pathway inhibitor efficacy in breast cancer.

Dysregulation of the PI3K-AKT-mTOR signaling network is a prominent feature of breast cancers. However, clinical responses to drugs targeting this pathway have been modest, possibly because of dynamic changes in cellular signaling that drive resistance and limit drug efficacy. Using a quantitative chemoproteomics approach, we mapped kinome dynamics in response to inhibitors of this pathway and identified signaling changes that correlate with drug sensitivity. Maintenance of AURKA after drug treatment was associated with resistance in breast cancer models. Incomplete inhibition of AURKA was a common source of therapy failure, and combinations of PI3K, AKT or mTOR inhibitors with the AURKA inhibitor MLN8237 were highly synergistic and durably suppressed mTOR signaling, resulting in apoptosis and tumor regression in vivo. This signaling map identifies survival factors whose presence limits the efficacy of targeted therapies and reveals new drug combinations that may unlock the full potential of PI3K-AKT-mTOR pathway inhibitors in breast cancer.

Microenvironment-Mediated Mechanisms of Resistance to HER2 Inhibitors Differ between HER2+ Breast Cancer Subtypes.

Extrinsic signals are implicated in breast cancer resistance to HER2-targeted tyrosine kinase inhibitors (TKIs). To examine how microenvironmental signals influence resistance, we monitored TKI-treated breast cancer cell lines grown on microenvironment microarrays composed of printed extracellular matrix proteins supplemented with soluble proteins. We tested ∼2,500 combinations of 56 soluble and 46 matrix microenvironmental proteins on basal-like HER2+ (HER2E) or luminal-like HER2+ (L-HER2+) cells treated with the TKIs lapatinib or neratinib. In HER2E cells, hepatocyte growth factor, a ligand for MET, induced resistance that could be reversed with crizotinib, an inhibitor of MET. In L-HER2+ cells, neuregulin1-ß1 (NRG1ß), a ligand for HER3, induced resistance that could be reversed with pertuzumab, an inhibitor of HER2-HER3 heterodimerization. The subtype-specific responses were also observed in 3D cultures and murine xenografts. These results, along with bioinformatic pathway analysis and siRNA knockdown experiments, suggest different mechanisms of resistance specific to each HER2+ subtype: MET signaling for HER2E and HER2-HER3 heterodimerization for L-HER2+ cells.

CCR5 Governs DNA Damage Repair and Breast Cancer Stem Cell Expansion.

The functional significance of the chemokine receptor CCR5 in human breast cancer epithelial cells is poorly understood. Here, we report that CCR5 expression in human breast cancer correlates with poor outcome. CCR5+ breast cancer epithelial cells formed mammospheres and initiated tumors with >60-fold greater efficiency in mice. Reintroduction of CCR5 expression into CCR5-negative breast cancer cells promoted tumor metastases and induced DNA repair gene expression and activity. CCR5 antagonists Maraviroc and Vicriviroc dramatically enhanced cell killing mediated by DNA-damaging chemotherapeutic agents. Single-cell analysis revealed CCR5 governs PI3K/Akt, ribosomal biogenesis, and cell survival signaling. As CCR5 augments DNA repair and is reexpressed selectively on cancerous, but not normal breast epithelial cells, CCR5 inhibitors may enhance the tumor-specific activities of DNA damage response-based treatments, allowing a dose reduction of standard chemotherapy and radiation.Significance: This study offers a preclinical rationale to reposition CCR5 inhibitors to improve the treatment of breast cancer, based on their ability to enhance the tumor-specific activities of DNA-damaging chemotherapies administered in that disease. Cancer Res; 78(7); 1657-71. ©2018 AACR.

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/.

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.].