A lineage tree-based hidden Markov model quantifies cellular heterogeneity and plasticity

TNBC response to paclitaxel phenocopies interferon response which reveals cell cycle-associated resistance mechanisms

Paclitaxel is a standard of care neoadjuvant therapy for patients with triple negative breast cancer (TNBC); however, it shows limited benefit for locally advanced or metastatic disease. Here we used a coordinated experimental-computational approach to explore the influence of paclitaxel on the cellular and molecular responses of TNBC cells. We found that escalating doses of paclitaxel resulted in multinucleation, promotion of senescence, and initiation of DNA damage induced apoptosis. Single-cell RNA sequencing (scRNA-seq) of TNBC cells after paclitaxel treatment revealed upregulation of innate immune programs canonically associated with interferon response and downregulation of cell cycle progression programs. Systematic exploration of transcriptional responses to paclitaxel and cancer-associated microenvironmental factors revealed common gene programs induced by paclitaxel, IFNB, and IFNG. Transcription factor (TF) enrichment analysis identified 13 TFs that were both enriched based on activity of downstream targets and also significantly upregulated after paclitaxel treatment. Functional assessment with siRNA knockdown confirmed that the TFs FOSL1, NFE2L2 and ELF3 mediate cellular proliferation and also regulate nuclear structure. We further explored the influence of these TFs on paclitaxel-induced cell cycle behavior via live cell imaging, which revealed altered progression rates through G1, S/G2 and M phases. We found that ELF3 knockdown synergized with paclitaxel treatment to lock cells in a G1 state and prevent cell cycle progression. Analysis of publicly available breast cancer patient data showed that high ELF3 expression was associated with poor prognosis and enrichment in programs associated with cell cycle progression. Together these analyses disentangle the diverse aspects of paclitaxel response and identify ELF3 upregulation as a putative biomarker of paclitaxel resistance in TNBC.

Survival Resonances during Fractional Killing of Cell Populations

Fractional killing in response to drugs is a hallmark of nongenetic cellular heterogeneity. Yet how individual lineages evade drug treatment, as observed in bacteria and cancer cells, is not quantitatively understood. We study a stochastic population model with age-dependent division and death rates, allowing for persistence. In periodic drug environments, we discover peaks in the survival probabilities at division or death times that are multiples of the environment duration. The survival resonances are unseen in unstructured populations and are amplified by persistence.

TNBC response to paclitaxel phenocopies interferon response which reveals cell cycle-associated resistance mechanisms

Paclitaxel is a standard of care neoadjuvant therapy for patients with triple negative breast cancer (TNBC); however, it shows limited benefit for locally advanced or metastatic disease. Here we used a coordinated experimental-computational approach to explore the influence of paclitaxel on the cellular and molecular responses of TNBC cells. We found that escalating doses of paclitaxel resulted in multinucleation, promotion of senescence, and initiation of DNA damage induced apoptosis. Single-cell RNA sequencing (scRNA-seq) of TNBC cells after paclitaxel treatment revealed upregulation of innate immune programs canonically associated with interferon response and downregulation of cell cycle progression programs. Systematic exploration of transcriptional responses to paclitaxel and cancer-associated microenvironmental factors revealed common gene programs induced by paclitaxel, IFNB, and IFNG. Transcription factor (TF) enrichment analysis identified 13 TFs that were both enriched based on activity of downstream targets and also significantly upregulated after paclitaxel treatment. Functional assessment with siRNA knockdown confirmed that the TFs FOSL1, NFE2L2 and ELF3 mediate cellular proliferation and also regulate nuclear structure. We further explored the influence of these TFs on paclitaxel-induced cell cycle behavior via live cell imaging, which revealed altered progression rates through G1, S/G2 and M phases. We found that ELF3 knockdown synergized with paclitaxel treatment to lock cells in a G1 state and prevent cell cycle progression. Analysis of publicly available breast cancer patient data showed that high ELF3 expression was associated with poor prognosis and enrichment programs associated with cell cycle progression. Together these analyses disentangle the diverse aspects of paclitaxel response and identify ELF3 upregulation as a putative biomarker of paclitaxel resistance in TNBC.

Quantifying Dynamic Phenotypic Heterogeneity in Resistant Escherichia coli under Translation-Inhibiting Antibiotics

Understanding the phenotypic heterogeneity of antibiotic-resistant bacteria following treatment and the transitions between different phenotypes is crucial for developing effective infection control strategies. The study expands upon previous work by explicating chloramphenicol-induced phenotypic heterogeneities in growth rate, gene expression, and morphology of resistant Escherichia coli using time-lapse microscopy. Correlating the bacterial growth rate and cspC expression, four interchangeable phenotypic subpopulations across varying antibiotic concentrations are identified, surpassing the previously described growth rate bistability. Notably, bacterial cells exhibiting either fast or slow growth rates can concurrently harbor subpopulations characterized by high and low gene expression levels, respectively. To elucidate the mechanisms behind this enhanced heterogeneity, a concise gene expression network model is proposed and the biological significance of the four phenotypes is further explored. Additionally, by employing Hidden Markov Model fitting and integrating the non-equilibrium landscape and flux theory, the real-time data encompassing diverse bacterial traits are analyzed. This approach reveals dynamic changes and switching kinetics in different cell fates, facilitating the quantification of observable behaviors and the non-equilibrium dynamics and thermodynamics at play. The results highlight the multi-dimensional heterogeneous behaviors of antibiotic-resistant bacteria under antibiotic stress, providing new insights into the compromised antibiotic efficacy, microbial response, and associated evolution processes.