Abstract:Single-cell RNA-sequencing (scRNA-seq) technology has provided significant insights into cancer drug resistance at the single-cell level. However, understanding dynamic cell transitions at the molecular systems level remains limited, requiring a systems biology approach. We present an approach that combines mathematical modeling with a pseudotime analysis using time-series scRNA-seq data obtained from the breast cancer cell line MCF-7 treated with tamoxifen. Our single-cell analysis identified five distinct subpopulations, including tamoxifen-sensitive and -resistant groups. Using a single-gene mathematical model, we discovered approximately 560–680 genes out of 6000 exhibiting multistable expression states in each subpopulation, including key estrogen-receptor-positive breast cancer cell survival genes, such as RPS6KB1. A bifurcation analysis elucidated their regulatory mechanisms, and we mapped these genes into a molecular network associated with cell survival and metastasis-related pathways. Our modeling approach comprehensively identifies key regulatory genes for drug resistance acquisition, enhancing our understanding of potential drug targets in breast cancer.

What problem does this paper attempt to address?

This paper attempts to address the problem of identifying key regulatory genes that lead to drug resistance in breast cancer. Specifically, the researchers aim to understand the dynamic molecular mechanisms of breast cancer cells during drug treatment, particularly how they transition from a drug-sensitive state to a drug-resistant state, by combining mathematical modeling and pseudotime analysis using single-cell time-series transcriptome data. The focus of the study is to identify key genes that exhibit multistable expression patterns during this process and to explore the regulatory mechanisms of these genes. ### Main Research Content: 1. **Background**: - Breast cancer is one of the most common cancers in women, especially ER-positive breast cancer (which relies on estrogen receptor growth). - Although endocrine therapies like tamoxifen have significantly improved survival rates, a considerable proportion of patients experience recurrence within 5-15 years, leading to metastasis and death. - The current understanding of drug resistance mechanisms is still limited, necessitating a systems biology approach for in-depth study. 2. **Methods**: - Single-cell RNA sequencing (scRNA-seq) technology was used to obtain time-series data of the MCF-7 cell line (an ER-positive breast cancer model) under continuous tamoxifen treatment. - Pseudotime analysis and mathematical modeling were employed to identify potential drug-sensitive and drug-resistant cell subpopulations. - The expression pattern of each gene was modeled to determine its stable states over time, thereby identifying multistable genes. - A transcriptional network model was constructed, encompassing pathways such as cell survival, adhesion, and epithelial-mesenchymal transition (EMT). 3. **Results**: - Approximately 560-680 genes were identified to exhibit multistable expression patterns in different subpopulations, including RPS6KB1, a gene associated with poor prognosis, cell survival, and growth. - Branching analysis revealed the regulatory mechanisms of these key regulatory genes and mapped them to molecular networks related to cell survival and metastasis. - It was found that the expression of RPS6KB1 in early pseudotime is mainly controlled by positive regulatory inhibitors, providing potential targets for intervention in drug resistance. 4. **Significance**: - This study provides a powerful tool for understanding drug resistance mechanisms in human diseases, particularly breast cancer. - The identified key regulatory genes can serve as potential drug targets, providing a basis for developing new therapeutic strategies. ### Summary: By combining mathematical modeling and pseudotime analysis, this paper successfully identifies key regulatory genes and their regulatory mechanisms that lead to drug resistance in breast cancer cells during drug treatment using single-cell transcriptome data. This approach not only deepens the understanding of drug resistance mechanisms but also offers new directions for future drug design and therapeutic strategies.

Identifying Key Regulatory Genes in Drug Resistance Acquisition: Modeling Pseudotime Trajectories of Breast Cancer Single-Cell Transcriptome

Identifying key regulatory genes in drug resistance acquisition: Modeling pseudotime trajectories of single-cell transcriptome

Modeling Time-Dependent Transcription Effects Of Her2 Oncogene And Discovery Of A Role For E2f2 In Breast Cancer Cell-Matrix Adhesion

Transformer-based modeling of Clonal Selection and Expression Dynamics (TraCSED) reveals resistance signatures in breast cancer

Integrating transcriptomics and bulk time course data into a mathematical framework to describe and predict therapeutic resistance in cancer

Differential Regulatory Network-Based Quantification and Prioritization of Key Genes Underlying Cancer Drug Resistance Based on Time-Course RNA-seq Data.

Modeling Expression Plasticity of Genes That Differentiate Drug-Sensitive from Drug-Resistant Cells to Chemotherapeutic Treatment

A Novel Mathematical Approach for Analysis of Integrated Cell–Patient Data Uncovers a 6-Gene Signature Linked to Endocrine Therapy Resistance

TGF-β Induced Reprogramming and Drug Resistance in Triple-Negative Breast Cells

Decoding gene regulatory circuitry underlying TNBC chemoresistance reveals biomarkers for therapy response and therapeutic targets

The prognostic genes model of breast cancer drug resistance based on single-cell sequencing analysis and transcriptome analysis

Identification and Pharmacological Targeting of Treatment-Resistant, Stem-like Breast Cancer Cells for Combination Therapy

Differential Gene Expression in Tamoxifen-Resistant Breast Cancer Cells Revealed by a New Analytical Model of RNA-Seq Data

Defining the Regulatory Logic of Breast Cancer Using Single-Cell Epigenetic and Transcriptome Profiling

Mathematical modeling and dynamic analysis for cancer resistance incorporating persister cells

Identification of New Chemoresistance-Associated Genes in Triple-Negative Breast Cancer by Single-Cell Transcriptomic Analysis

A mechanistic model captures the emergence and implications of non-genetic heterogeneity and reversible drug resistance in ER+ breast cancer cells

Single-Cell Image-Based Analysis Reveals Chromatin Changes During the Acquisition of Tamoxifen Drug Resistance.

Mathematical Modeling of Therapy-induced Cancer Drug Resistance: Connecting Cancer Mechanisms to Population Survival Rates

A Kinetic Model of Multiple Phenotypic States for Breast Cancer Cells.

Pseudo-Temporal Dynamics of Chemo-Resistant Triple-Negative Breast Cancer Cells Reveal EGFR/HER2 Inhibition As Synthetic Lethal During Mid-neoadjuvant Chemotherapy