Modeling differentiation-state transitions linked to therapeutic escape in triple-negative breast cancer

Exploring the resistance mechanism of triple-negative breast cancer to paclitaxel through the scRNA-seq analysis

BACKGROUND

The triple negative breast cancer (TNBC) is the most malignant subtype of breast cancer with high aggressiveness. Although paclitaxel-based chemotherapy scenario present the mainstay in TNBC treatment, paclitaxel resistance is still a striking obstacle for cancer cure. So it is imperative to probe new therapeutic targets through illustrating the mechanisms underlying paclitaxel chemoresistance.

METHODS

The Single cell RNA sequencing (scRNA-seq) data of TNBC cells treated with paclitaxel at different points were downloaded from the Gene Expression Omnibus (GEO) database. The Seurat R package was used to filter and integrate the scRNA-seq expression matrix. Cells were further clustered by the FindClusters function, and the gene marker of each subset was defined by FindAllMarkers function. Then, the hallmark score of each cell was calculated by AUCell R package, the biological function of the highly expressed interest genes was analyzed by the DAVID database. Subsequently, we performed pseudotime analysis to explore the change patterns of drug resistance genes and SCENIC analysis to identify the key transcription factors (TFs). Finally, the inhibitors of which were also analyzed by the CTD database.

RESULTS

We finally obtained 6 cell subsets from 2798 cells, which were marked as AKR1C3+, WNT7A+, FAM72B+, RERG+, IDO1+ and HEY1+HCC1143 cell subsets, among which the AKR1C3+, IDO1+ and HEY1+ cell subsets proportions increased with increasing treatment time, and then were regarded as paclitaxel resistance subsets. Hallmark score and pseudotime analysis showed that these paclitaxel resistance subsets were associated with the inflammatory response, virus and interferon response activation. In addition, the gene regulatory networks (GRNs) indicated that 3 key TFs (STAT1, CEBPB and IRF7) played vital role in promoting resistance development, and five common inhibitors targeted these TFs as potential combination therapies of paclitaxel were identified.

CONCLUSION

In this study, we identified 3 paclitaxel resistance relevant IFs and their inhibitors, which offers essential molecular basis for paclitaxel resistance and beneficial guidance for the combination of paclitaxel in clinical TNBC therapy.

Quantifying cancer cell plasticity with gene regulatory networks and single-cell dynamics

Phenotypic plasticity of cancer cells can lead to complex cell state dynamics during tumor progression and acquired resistance. Highly plastic stem-like states may be inherently drug-resistant. Moreover, cell state dynamics in response to therapy allow a tumor to evade treatment. In both scenarios, quantifying plasticity is essential for identifying high-plasticity states or elucidating transition paths between states. Currently, methods to quantify plasticity tend to focus on 1) quantification of quasi-potential based on the underlying gene regulatory network dynamics of the system; or 2) inference of cell potency based on trajectory inference or lineage tracing in single-cell dynamics. Here, we explore both of these approaches and associated computational tools. We then discuss implications of each approach to plasticity metrics, and relevance to cancer treatment strategies.

Morphodynamical cell state description via live-cell imaging trajectory embedding

Time-lapse imaging is a powerful approach to gain insight into the dynamic responses of cells, but the quantitative analysis of morphological changes over time remains challenging. Here, we exploit the concept of "trajectory embedding" to analyze cellular behavior using morphological feature trajectory histories-that is, multiple time points simultaneously, rather than the more common practice of examining morphological feature time courses in single timepoint (snapshot) morphological features. We apply this approach to analyze live-cell images of MCF10A mammary epithelial cells after treatment with a panel of microenvironmental perturbagens that strongly modulate cell motility, morphology, and cell cycle behavior. Our morphodynamical trajectory embedding analysis constructs a shared cell state landscape revealing ligand-specific regulation of cell state transitions and enables quantitative and descriptive models of single-cell trajectories. Additionally, we show that incorporation of trajectories into single-cell morphological analysis enables (i) systematic characterization of cell state trajectories, (ii) better separation of phenotypes, and (iii) more descriptive models of ligand-induced differences as compared to snapshot-based analysis. This morphodynamical trajectory embedding is broadly applicable to the quantitative analysis of cell responses via live-cell imaging across many biological and biomedical applications.

Mathematical characterization of population dynamics in breast cancer cells treated with doxorubicin

The development of chemoresistance remains a significant cause of treatment failure in breast cancer. We posit that a mathematical understanding of chemoresistance could assist in developing successful treatment strategies. Towards that end, we have developed a model that describes the cytotoxic effects of the standard chemotherapeutic drug doxorubicin on the MCF-7 breast cancer cell line. We assume that treatment with doxorubicin induces a compartmentalization of the breast cancer cell population into surviving cells, which continue proliferating after treatment, and irreversibly damaged cells, which gradually transition from proliferating to treatment-induced death. The model is fit to experimental data including variations in drug concentration, inter-treatment interval, and number of doses. Our model recapitulates tumor cell dynamics in all these scenarios (as quantified by the concordance correlation coefficient, CCC > 0.95). In particular, superior tumor control is observed with higher doxorubicin concentrations, shorter inter-treatment intervals, and a higher number of doses (p < 0.05). Longer inter-treatment intervals require adapting the model parameterization after each doxorubicin dose, suggesting the promotion of chemoresistance. Additionally, we propose promising empirical formulas to describe the variation of model parameters as functions of doxorubicin concentration (CCC > 0.78). Thus, we conclude that our mathematical model could deepen our understanding of the cytotoxic effects of doxorubicin and could be used to explore practical drug regimens achieving optimal tumor control.

Strategies to Reduce Long-Term Drug Resistance by Considering Effects of Differential Selective Treatments

Despite great advances in modeling and cancer therapy using optimal control theory, tumor heterogeneity and drug resistance are major obstacles in cancer treatments. Since recent biological studies demonstrated the evidence of tumor heterogeneity and assessed potential biological and clinical implications, tumor heterogeneity should be taken into account in the optimal control problem to improve treatment strategies. Here, first we study the effects of two different treatment strategies (i.e., symmetric and asymmetric) in a minimal two-population model to examine the long-term effects of these treatment methods on the system. Second, by considering tumor adaptation to treatment as a factor of the cost function, the optimal treatment strategy is derived. Numerical examples show that optimal treatment decreases tumor burden for the long-term by decreasing rate of tumor adaptation over time.

Measuring and Modelling the Epithelial- Mesenchymal Hybrid State in Cancer: Clinical Implications

The epithelial-mesenchymal (E/M) hybrid state has emerged as an important mediator of elements of cancer progression, facilitated by epithelial mesenchymal plasticity (EMP). We review here evidence for the presence, prognostic significance, and therapeutic potential of the E/M hybrid state in carcinoma. We further assess modelling predictions and validation studies to demonstrate stabilised E/M hybrid states along the spectrum of EMP, as well as computational approaches for characterising and quantifying EMP phenotypes, with particular attention to the emerging realm of single-cell approaches through RNA sequencing and protein-based techniques.

Changes in Tumor Stem Cell Markers and Epithelial-Mesenchymal Transition Markers in Nonluminal Breast Cancer after Neoadjuvant Chemotherapy and Their Correlation with Contrast-Enhanced Ultrasound

Nonluminal breast cancer has high early metastasis and treatment resistance, and neoadjuvant chemotherapy (NAC) is needed. The presence of cancer stem cells (CSC) and epithelial-mesenchymal transition (EMT) leads to poor prognosis. This study investigated the changes in CSC markers and EMT markers after NAC in nonluminal breast cancer and their correlation with contrast-enhanced ultrasound (CEUS) features and chemotherapy efficacy. Before NAC, the range of nonluminal breast cancer on CEUS was larger than that of two-dimensional ultrasound, but after NAC, it was significantly smaller than that of two-dimensional ultrasound and closer to the postoperative pathological size. After NAC, the enlarged lesions and perfusion defects were significantly less than those before NAC. The time-intensity curve showed the characteristics of slow-in, low enhancement, and low perfusion. Nonluminal breast cancer downregulated the expression of CSC markers and EMT markers after NAC, but the epithelial phenotype of nonluminal breast cancer with good response to chemotherapy was upregulated. In nonluminal breast cancer with poor response to chemotherapy, markers of CSC and EMT were highly expressed before chemotherapy. In conclusion, CEUS is better than conventional ultrasound in estimating NAC efficacy in this mode. CEUS can also predict the prognosis of nonluminal breast cancer before NAC with the characteristics of enhanced enlargement and perfusion defects. The contrast-enhanced time-intensity curve of lesions with relatively poor blood supply may have more CSC and EMT characteristics.

Cancer stem cell and mesenchymal cell cooperative actions in metastasis progression and hormone resistance in prostate cancer: Potential role of androgen and gonadotropin‑releasing hormone receptors (Review)

Prostate cancer (PCa) is the leading cause of male cancer‑associated mortality worldwide. Mortality is associated with metastasis and hormone resistance. Cellular, genetic and molecular mechanisms underlying metastatic progression and hormone resistance are poorly understood. Studies have investigated the local effects of gonadotropin‑releasing hormone (GnRH) analogs (used for androgen deprivation treatments) and the presence of the GnRH receptor (GnRH‑R) on PCa cells. Furthermore, cell subpopulations with stem‑like properties, or cancer stem cells, have been isolated and characterized using a cell culture system derived from explants of human prostate tumors. In addition, the development of preclinical orthotopic models of human PCa in a nonobese diabetic/severe combined immunodeficiency mouse model of compromised immunity has enabled the establishment of a reproducible system of metastatic progression in vivo. There is increasing evidence that metastasis is a complex process involving the cooperative actions of different cancer cell subpopulations, in which cancer stem‑like cells would be responsible for the final step of colonizing premetastatic niches. It has been hypothesized that PCa cells with stemness and mesenchymal signatures act cooperatively in metastatic progression and the inhibition of stemness genes, and that overexpression of androgen receptor (AR) and GnRH‑R decreases the rate the metastasis and sensitizes tumors to hormone therapy. The aim of the present review is to analyze the evidence regarding this cooperative process and the possible influence of stem‑like cell phenotypes, AR and GnRH‑R in metastatic progression and hormone resistance. These aspects may represent an important contribution in the understanding of the mechanisms underlying metastasis and hormone resistance in PCa, and potential routes to blocking these processes, enabling the development of novel therapies that would be particularly relevant for patients with metastatic and castration‑resistant PCa.

A mechanism for epithelial-mesenchymal heterogeneity in a population of cancer cells

Epithelial-mesenchymal heterogeneity implies that cells within the same tumor can exhibit different phenotypes-epithelial, mesenchymal, or one or more hybrid epithelial-mesenchymal phenotypes. This behavior has been reported across cancer types, both in vitro and in vivo, and implicated in multiple processes associated with metastatic aggressiveness including immune evasion, collective dissemination of tumor cells, and emergence of cancer cell subpopulations with stem cell-like properties. However, the ability of a population of cancer cells to generate, maintain, and propagate this heterogeneity has remained a mystifying feature. Here, we used a computational modeling approach to show that epithelial-mesenchymal heterogeneity can emerge from the noise in the partitioning of biomolecules (such as RNAs and proteins) among daughter cells during the division of a cancer cell. Our model captures the experimentally observed temporal changes in the fractions of different phenotypes in a population of murine prostate cancer cells, and describes the hysteresis in the population-level dynamics of epithelial-mesenchymal plasticity. The model is further able to predict how factors known to promote a hybrid epithelial-mesenchymal phenotype can alter the phenotypic composition of a population. Finally, we used the model to probe the implications of phenotypic heterogeneity and plasticity for different therapeutic regimens and found that co-targeting of epithelial and mesenchymal cells is likely to be the most effective strategy for restricting tumor growth. By connecting the dynamics of an intracellular circuit to the phenotypic composition of a population, our study serves as a first step towards understanding the generation and maintenance of non-genetic heterogeneity in a population of cancer cells, and towards the therapeutic targeting of phenotypic heterogeneity and plasticity in cancer cell populations.

The Physics of Cellular Decision Making During Epithelial-Mesenchymal Transition

The epithelial-mesenchymal transition (EMT) is a process by which cells lose epithelial traits, such as cell-cell adhesion and apico-basal polarity, and acquire migratory and invasive traits. EMT is crucial to embryonic development and wound healing. Misregulated EMT has been implicated in processes associated with cancer aggressiveness, including metastasis. Recent experimental advances such as single-cell analysis and temporal phenotypic characterization have established that EMT is a multistable process wherein cells exhibit and switch among multiple phenotypic states. This is in contrast to the classical perception of EMT as leading to a binary choice. Mathematical modeling has been at the forefront of this transformation for the field, not only providing a conceptual framework to integrate and analyze experimental data, but also making testable predictions. In this article, we review the key features and characteristics of EMT dynamics, with a focus on the mathematical modeling approaches that have been instrumental to obtaining various useful insights.

Systems biology approaches to measure and model phenotypic heterogeneity in cancer

The recent wide-spread adoption of single cell profiling technologies has revealed that individual cancers are not homogenous collections of deregulated cells, but instead are comprised of multiple genetically and phenotypically distinct cell subpopulations that exhibit a wide range of responses to extracellular signals and therapeutic insult. Such observations point to the urgent need to understand cancer as a complex, adaptive system. Cancer systems biology studies seek to develop the experimental and theoretical methods required to understand how biological components work together to determine how cancer cells function. Ultimately, such approaches will lead to improvements in how cancer is managed and treated. In this review, we discuss recent advances in cancer systems biology approaches to quantify, model, and elucidate mechanisms of heterogeneity.

Interrogation of Phenotypic Plasticity between Epithelial and Mesenchymal States in Breast Cancer

Dynamic interconversions between transitional epithelial and mesenchymal states underpin the epithelial mesenchymal plasticity (EMP) seen in some carcinoma cell systems. We have delineated epithelial and mesenchymal subpopulations existing within the PMC42-LA breast cancer cell line by their EpCAM expression. These purified but phenotypically plastic states, EpCAMHigh (epithelial) and EpCAMLow (mesenchymal), have the ability to regain the phenotypic equilibrium of the parental population (i.e., 80% epithelial and 20% mesenchymal) over time, although the rate of reversion in the mesenchymal direction (epithelial-mesenchymal transition; EMT) is higher than that in the epithelial direction (mesenchymal-epithelial transition; MET). Single-cell clonal propagation was implemented to delineate the molecular and cellular features of this intrinsic heterogeneity with respect to EMP flux. The dynamics of the phenotypic proportions of epithelial and mesenchymal states in single-cell generated clones revealed clonal diversity and intrinsic plasticity. Single cell-derived clonal progenies displayed differences in their functional attributes of proliferation, stemness marker (CD44/CD24), migration, invasion and chemo-sensitivity. Interrogation of genomic copy number variations (CNV) with whole exome sequencing (WES) in the context of chromosome count from metaphase spread indicated that chromosomal instability was not influential in driving intrinsic phenotypic plasticity. Overall, these findings reveal the stochastic nature of both the epithelial and mesenchymal subpopulations, and the single cell-derived clones for differential functional attributes.