A Mathematical Model of Breast Tumor Progression Based on Immune Infiltration

Breast cancer is the most prominent type of cancer among women. Understanding the microenvironment of breast cancer and the interactions between cells and cytokines will lead to better treatment approaches for patients. In this study, we developed a data-driven mathematical model to investigate the dynamics of key cells and cytokines involved in breast cancer development. We used gene expression profiles of tumors to estimate the relative abundance of each immune cell and group patients based on their immune patterns. Dynamical results show the complex interplay between cells and molecules, and sensitivity analysis emphasizes the direct effects of macrophages and adipocytes on cancer cell growth. In addition, we observed the dual effect of IFN-γ role= presentation style= box-sizing: border-box; max-height: none; display: inline; line-height: normal; word-spacing: normal; overflow-wrap: normal; white-space: nowrap; float: none; direction: ltr; max-width: none; min-width: 0px; min-height: 0px; border: 0px; padding: 0px; margin: 0px; position: relative; \u3eγ on cancer proliferation, either through direct inhibition of cancer cells or by increasing the cytotoxicity of CD8+ T-cells

Panobinostat Induced Spatial In Situ Biomarkers Predictive of Anti-PD-1 Efficacy in Mouse Mammary Carcinoma

Immunotherapies, including anti-PD-1 immune checkpoint blocking (ICB) antibodies, have revolutionized the treatment of many solid malignancies. However, their efficacy in breast cancer has been limited to a subset of patients with triple-negative breast cancer, where ICBs are routinely combined with a range of cytotoxic and targeted agents. Reliable biomarkers predictive of the therapeutic response to ICB in breast cancer are critically missing, though a combination response has been associated with immunogenic cell death (ICD). Here, we utilized a recently developed integrated analytical platform, the multiplex implantable microdevice assay (MIMA), to evaluate the presence and spatial cell relations of literature-based candidate markers predictive of ICB efficacy in luminal mouse mammary carcinoma. MIMA integrates (i) an implantable microdevice for the localized delivery of small amounts of drugs inside the tumor bed with (ii) sequential multiplex immunohistochemistry (mIHC) and spatial cell analysis pipelines to rapidly (within days) describe drug mechanisms of action and find predictive biomarkers in complex tumor tissue. We show that the expression of cleaved caspase-3, ICAM-1, neuropilin-1, myeloperoxidase, calreticulin, galectin-3, and PD-L1 were spatially associated with the efficacy of panobinostat, a pan-HDAC inhibitor that was previously shown to induce immunogenic cell death and synergize with anti-PD-1 in breast cancer. PD-L1 by itself, however, was not a reliable predictor. Instead, ICB efficacy was robustly identified through the in situ hotspot detection of galectin-3-positive non-proliferating tumor zones enriched in cell death and infiltrated by anti-tumor cytotoxic neutrophils positive for ICAM-1 and neuropilin-1. Such hotspots can be specifically detected using distance-based cluster analyses. Single-cell measurements of the functional states in the tumor microenvironment suggest that both qualitative and quantitative effects might drive effective therapy responses. Overall, the presented study provides (i) complementary biological knowledge about the earliest cell events of induced anti-tumor immunity in breast cancer, including the emergence of resistant cancer stem cells, and (ii) newly identified biomarkers in form of specific spatial cell associations. The approach used standard cell-type-, IHC-, and FFPE-based techniques, and therefore the identified spatial clustering of in situ biomarkers can be readily integrated into existing clinical or research workflows, including in luminal breast cancer. Since early drug responses were detected, the biomarkers could be especially applicable to window-of-opportunity clinical trials to rapidly discriminate between responding and resistant patients, thus limiting unnecessary treatment-associated toxicities

Panobinostat Induced Spatial In Situ Biomarkers Predictive of Anti-PD-1 Efficacy in Mouse Mammary Carcinoma

Immunotherapies, including anti-PD-1 immune checkpoint blocking (ICB) antibodies, have revolutionized the treatment of many solid malignancies. However, their efficacy in breast cancer has been limited to a subset of patients with triple-negative breast cancer, where ICBs are routinely combined with a range of cytotoxic and targeted agents. Reliable biomarkers predictive of the therapeutic response to ICB in breast cancer are critically missing, though a combination response has been associated with immunogenic cell death (ICD). Here, we utilized a recently developed integrated analytical platform, the multiplex implantable microdevice assay (MIMA), to evaluate the presence and spatial cell relations of literature-based candidate markers predictive of ICB efficacy in luminal mouse mammary carcinoma. MIMA integrates (i) an implantable microdevice for the localized delivery of small amounts of drugs inside the tumor bed with (ii) sequential multiplex immunohistochemistry (mIHC) and spatial cell analysis pipelines to rapidly (within days) describe drug mechanisms of action and find predictive biomarkers in complex tumor tissue. We show that the expression of cleaved caspase-3, ICAM-1, neuropilin-1, myeloperoxidase, calreticulin, galectin-3, and PD-L1 were spatially associated with the efficacy of panobinostat, a pan-HDAC inhibitor that was previously shown to induce immunogenic cell death and synergize with anti-PD-1 in breast cancer. PD-L1 by itself, however, was not a reliable predictor. Instead, ICB efficacy was robustly identified through the in situ hotspot detection of galectin-3-positive non-proliferating tumor zones enriched in cell death and infiltrated by anti-tumor cytotoxic neutrophils positive for ICAM-1 and neuropilin-1. Such hotspots can be specifically detected using distance-based cluster analyses. Single-cell measurements of the functional states in the tumor microenvironment suggest that both qualitative and quantitative effects might drive effective therapy responses. Overall, the presented study provides (i) complementary biological knowledge about the earliest cell events of induced anti-tumor immunity in breast cancer, including the emergence of resistant cancer stem cells, and (ii) newly identified biomarkers in form of specific spatial cell associations. The approach used standard cell-type-, IHC-, and FFPE-based techniques, and therefore the identified spatial clustering of in situ biomarkers can be readily integrated into existing clinical or research workflows, including in luminal breast cancer. Since early drug responses were detected, the biomarkers could be especially applicable to window-of-opportunity clinical trials to rapidly discriminate between responding and resistant patients, thus limiting unnecessary treatment-associated toxicities

A PDE Model of Breast Tumor Progression in MMTV-PyMT Mice

The evolution of breast tumors greatly depends on the interaction network among different cell types, including immune cells and cancer cells in the tumor. This study takes advantage of newly collected rich spatio-temporal mouse data to develop a data-driven mathematical model of breast tumors that considers cells’ location and key interactions in the tumor. The results show that cancer cells have a minor presence in the area with the most overall immune cells, and the number of activated immune cells in the tumor is depleted over time when there is no influx of immune cells. Interestingly, in the case of the influx of immune cells, the highest concentrations of both T cells and cancer cells are in the boundary of the tumor, as we use the Robin boundary condition to model the influx of immune cells. In other words, the influx of immune cells causes a dominant outward advection for cancer cells. We also investigate the effect of cells’ diffusion and immune cells’ influx rates in the dynamics of cells in the tumor micro-environment. Sensitivity analyses indicate that cancer cells and adipocytes’ diffusion rates are the most sensitive parameters, followed by influx and diffusion rates of cytotoxic T cells, implying that targeting them is a possible treatment strategy for breast cancer

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

SUMMARY 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

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

Abstract 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