Search PubMedSearch

SEARCH · Search PubMed

Results for “cell state transition”

Search indexed PubMed citations on genomics, clinical trials, systematic reviews and public health. Explore titles, authors and supplied subject terms, then open the PubMed record.

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 19 recordsLinked to original sources

Identifying fate-determining transcription factors with single-cell omics.

Single-cell sequencing enables the systematic discovery of cell fate-determining transcription factors (TFs), or key TFs, that define cellular identity or drive cell state transitions. A wide range of computational methods have been developed for this goal, but they differ substantially in the input data and the biological questions they address. In this article, we systematically review computational approaches for key TF identification and organize them from three perspectives: whether they identify TFs defining cell state identity or driving state transitions, whether transitions are modeled as discrete or continuous processes, and whether TFs act individually or combinatorially. We summarize key features and application scenarios of relevant methods to guide tool selection and discuss emerging trends in this field toward programmable and active control of cell fate.

Transcription Factors

Trajectory inference from single-cell genomics data with a process time model.

Single-cell transcriptomics experiments provide gene expression snapshots of heterogeneous cell populations across cell states. These snapshots have been used to infer trajectories and dynamic information even without intensive, time-series data by ordering cells according to gene expression similarity. However, while single-cell snapshots sometimes offer valuable insights into dynamic processes, current methods for ordering cells are limited by descriptive notions of "pseudotime" that lack intrinsic physical meaning. Instead of pseudotime, we propose inference of "process time" via a principled modeling approach to formulating trajectories and inferring latent variables corresponding to timing of cells subject to a biophysical process. Our implementation of this approach, called Chronocell, provides a biophysical formulation of trajectories built on cell state transitions. The Chronocell model is identifiable, making parameter inference meaningful. Furthermore, Chronocell can interpolate between trajectory inference, when cell states lie on a continuum, and clustering, when cells cluster into discrete states. By using a variety of datasets ranging from cluster-like to continuous, we show that Chronocell enables us to assess the suitability of datasets and reveals distinct cellular distributions along process time that are consistent with biological process times. We also compare our parameter estimates of degradation rates to those derived from metabolic labeling datasets, thereby showcasing the biophysical utility of Chronocell. Nevertheless, based on performance characterization on simulations, we find that process time inference can be challenging, highlighting the importance of dataset quality and careful model assessment.

Single-Cell Analysis

Deciphering Cell Fate and Clonal Dynamics via Integrative Single-Cell Lineage Modeling.

Through natural or synthetic lineage barcodes, single-cell technologies now enable the joint measurement of molecular states and clonal identities, providing an unprecedented opportunity to study cell fate and dynamics. Yet, most computational methods for inferring cell development and differentiation rely exclusively on transcriptional similarity, overlooking the lineage information encoded by lineage barcodes. This limitation is exemplified by T cells, where subtle transcriptional differences mark divergent fates with distinct biological activity. Single-cell RNA and matched TCR sequencing is now ubiquitous in the analysis of clinical samples, where the TCR sequence provides an endogenous clonal barcode and could reveal clonal T cell responses. We present Clonotrace, a computational framework that jointly models gene expression and clonotype information to infer cell state transitions and fate biases with higher fidelity. While motivated by challenges in analyzing T cell populations, especially in the tumor microenvironment and immunotherapy settings, Clonotrace is broadly applicable to any lineage-barcoded single-cell dataset. Across diverse systems including T cells, hematopoietic differentiation, and cancer therapy resistance models, Clonotrace reveals differentiation hierarchies, distinguishes unipotent from multipotent states, and identifies candidate fate-determining genes driving lineage commitment.

Journal Article

Deciphering Cell Fate and Clonal Dynamics via Integrative Single-Cell Lineage Modeling.

Through natural or synthetic lineage barcodes, single-cell technologies now enable the joint measurement of molecular states and clonal identities, providing an unprecedented opportunity to study cell fate and dynamics. Yet, most computational methods for inferring cell development and differentiation rely exclusively on transcriptional similarity, overlooking the lineage information encoded by lineage barcodes. This limitation is exemplified by T cells, where subtle transcriptional differences mark divergent fates with distinct biological activity. Single-cell RNA and matched TCR sequencing is now ubiquitous in the analysis of clinical samples, where the TCR sequence provides an endogenous clonal barcode and could reveal clonal T cell responses. We present Clonotrace, a computational framework that jointly models gene expression and clonotype information to infer cell state transitions and fate biases with higher fidelity. While motivated by challenges in analyzing T cell populations, especially in the tumor microenvironment and immunotherapy settings, Clonotrace is broadly applicable to any lineage-barcoded single-cell dataset. Across diverse systems including T cells, hematopoietic differentiation, and cancer therapy resistance models, Clonotrace reveals differentiation hierarchies, distinguishes unipotent from multipotent states, and identifies candidate fate-determining genes driving lineage commitment.

Journal Article

Spatial proximity sequencing maps developmental dynamics in the germinal center.

Spatial profiling of proteins and protein interactions facilitates understanding of cell functions within tissues and is essential for studies in signaling, immunity, and cancer. We present spatial proximity sequencing (Sprox-seq) for simultaneous profiling of surface proteins, protein complexes, and mRNAs, recording the tissue location of each molecule. Sprox-seq profiled 32 proteins, 528 pairwise interactions, and thousands of mRNAs with spatial resolution across human tonsils and germinal centers. Mapping tissue-wide protein interactions recapitulated RNA-defined tissue architecture but also revealed higher interaction complexity in the light zone. Protein-interaction trajectories uncovered a B cell state transition distinct from that inferred by RNA. Integrated protein-complex and mRNA analysis related spatially enriched complexes with mitotic pathways. Sprox-seq captured cell-cell interactions, such as B cell-follicular dendritic cell interactions mediated by the receptor complex VLA-4-VCAM1. Sprox-seq provides a spatially resolved multi-modal view of cell states and an integrated study of protein and cellular interactions across tissues.

Humans

Cell-type-specific defense priming and NHP-dependent systemic immunity against Pectobacterium in Chinese cabbage.

Defense priming is an effective strategy for mounting the defensive capacity of plants. Primed plants undergo minimal changes, enabling rapid, robust responses to subsequent pathogen attack. The cell-type-specific mechanisms underlying priming and systemic immunity remain unclear. Using single-nucleus RNA sequencing (snRNA-seq), we reveal how primed Chinese cabbage coordinates N-hydroxypipecolic acid (NHP)-dependent systemic immunity against Pectobacterium carotovorum. NHP accumulates as a phloem-mobile priming signal. snRNA-seq identified a priming-specific epidermal cell state transition and transcriptional cascade: BrWRKY18-1 activates BrWRKY33-1 expression during priming, initiating NHP-dependent systemic immunity. At post-challenge primed state, pathogenesis-related genes (PRs), such as BrPR3, BrPR4-1, and BrPR4-2, are predominantly expressed in systemic tissues to inhibit infection. Notably, the functionality of these PR genes in immunity necessitates their expression in all cells to effectively enhance defensive resistance. This NHP-transcription factor-PR axis couples mobile signaling to systemic immunity, elucidating cell-type-specific defense priming.

Plant Immunity

Identification of Drug-resistant Cell Subpopulations in Colorectal Cancer Through Single-cell Analysis and Exploration of Potential Therapeutic Strategies.

INTRODUCTION: The therapeutic efficacy of Colorectal Cancer (CRC) is often compromised by resistance to the standard chemotherapy agent oxaliplatin. METHODS: This study obtained single-cell RNA sequencing (scRNA-seq) data from the Gene Expression Omnibus (GEO) database. Differentially Expressed Genes (DEGs) between resistant and sensitive epithelial subpopulations were identified, followed by enrichment analysis. Pseudotemporal trajectory and cell-cell communication were analyzed using Monocle2 and CellChat, respectively. The candidate drug was predicted by Connectivity Map (cMAP) analysis. External validation included assessment of the EpC2 signature in an oxaliplatin-resistant cell line dataset (GSE76092), survival analysis using The Cancer Genome Atlas (TCGA) cohorts, and re-analysis of the GSE179784 dataset to assess the reproducibility of EpC2-like subpopulations and their DNA Damage Repair (DDR) scores. RESULTS: Cell subpopulations were divided into 10 clusters. Among them, epithelial cells comprised 5 subpopulations, with EPC2 identified as a potential oxaliplatin-resistant subset. DEGs were enriched in the TNF and IL-17 pathways. External validation confirmed the enrichment of EpC2 in resistant cell lines and its association with poor survival. Pseudotemporal trajectory revealed that epithelial cells underwent state transitions, forming two distinct branches. The resistant group exhibited enrichment in RNA splicing and NF-κB pathways. Cell-cell communication analysis revealed interactions involving MDK- NCL and PPIA-BSG. Dasatinib was predicted as a candidate drug. DISCUSSION: We identified an oxaliplatin-resistant subpopulation of Epithelial Cells (EpC2) in CRC, elucidated its multi-layered resistance mechanisms, and integrated multi- omics and cMAP database analyses to predict a potential intervention drug. CONCLUSION: This study provided potential therapeutic possibilities for oxaliplatin resistance, contributing to CRC treatment.

Humans

Genome-wide chromatin recording resolves dynamic cell state changes.

Understanding how the chromatin state of a cell influences its future behavior is a major challenge throughout biology. However, most chromatin profiling methods are limited to endpoint assays. Here, we present LagTag, a method for recovery of earlier and endpoint chromatin states in the same mammalian cells. In this approach, transient activity of bacterial adenine methyltransferase fusions records the DNA-binding profiles of chromatin-associated proteins of interest at earlier time points. Subsequent tagmentation and sequencing recover the earlier chromatin profile from adenine methylation profiles, alongside endpoint profiles of endogenous chromatin-associated proteins. We verified that LagTag profiles aligned with those from established methods in mouse and human cells. We then applied LagTag to record and recover dynamic chromatin state transitions during mouse embryonic stem cell differentiation, capturing transcriptional signatures from pre- and post-differentiation time points within the same cell population. LagTag thus provides a foundation for temporally resolved chromatin profiling. A record of this paper's transparent peer review process is included in the supplemental information.

chromatin organization

Genome-wide chromatin recording resolves dynamic cell state changes.

Understanding how the chromatin state of a cell influences its future behavior is a major challenge throughout biology. However, most chromatin profiling methods are limited to endpoint assays. Here, we present LagTag, a method for recovery of earlier and endpoint chromatin states in the same mammalian cells. In this approach, transient expression of bacterial adenine methyltransferase fusions records the DNA binding profiles of chromatin-associated proteins of interest at earlier timepoints. Subsequent tagmentation and sequencing recovers the earlier chromatin profile from adenine methylation profiles, alongside endpoint profiles of endogenous chromatin-associated proteins. We verified that LagTag profiles aligned with those from established methods in mouse and human cells. We then applied LagTag to record and recover dynamic chromatin state transitions during mouse embryonic stem cell differentiation, capturing transcriptional signatures from pre- and post-differentiation timepoints within the same cell population. LagTag thus provides a foundation for temporally resolved chromatin profiling.

Journal Article

Defining active and repressive chromatin states in neural crest cells using low-input CUT&RUN.

The transition of neural crest cells (NCCs) from a multipotent state to lineage-restricted derivatives, including melanocytes, is governed by tightly regulated epigenetic mechanisms that orchestrate cell type specific gene expression programs. Histone post-translational modifications (PTMs), in particular, play an important role in modulating chromatin accessibility, enhancer activation, and transcription factor occupancy, thereby facilitating dynamic chromatin and transcriptional reprogramming required during development. However, profiling such chromatin states in rare and transient Neural Crest Cell (NCC) populations in vivo remains technically challenging. To address this, we present an optimized low-input Cleavage Under Targets and Release Using Nuclease (CUT&RUN) workflow tailored for fluorescence-activated cell sorting (FACS) isolated NCCs from zebrafish embryos. This approach enables high-resolution and low-background mapping of key histone modifications, including H3K27ac, H3K4me3, and H3K27me3, from limited cell numbers. Collectively, these methodologies provide a robust framework for dissecting chromatin state dynamics in developmental systems and can also offer insights into epigenetic dysregulation associated with disease.

Animals

Mirror worlds: The shared regulatory architecture of cell fate in development and cancer.

Lineage plasticity has emerged as a central mechanism through which cancer cells adapt to therapeutic pressure, evade immune surveillance, and acquire aggressive phenotypes. Although recognized across tumor types, the regulatory principles governing how cancer cells reprogram cellular identity remain incompletely understood. In this review, we propose that lineage plasticity in cancer reflects the redeployment of regulatory frameworks established during normal development. Rather than representing a stochastic byproduct of genomic instability, cancer plasticity frequently unfolds within gene regulatory architectures that also govern cell fate specification, lineage commitment, and controlled state transitions during embryogenesis and tissue homeostasis. Developmental transcription factors, including members of the SOX family, FOXA1, ASCL1, NKX2-1, and epithelial-mesenchymal transition regulators, function as lineage gatekeepers during development but are repurposed in cancer to destabilize lineage commitment and enable phenotypic switching. Similarly, epigenetic regulators that guide developmental trajectories, including chromatin remodeling complexes, Polycomb group proteins, and DNA methylation machinery, are frequently dysregulated or redistributed in tumors, altering the repression of lineage-stabilizing and alternative lineage programs and thereby weakening epigenetic barriers to lineage transitions. Together, these observations support a model in which development and cancer operate as mirror regulatory systems: one establishing and stabilizing cellular identity, the other exploiting the same regulatory architecture to permit adaptive reprogramming under selective pressure. We further discuss how emerging single-cell and spatial multi-omics technologies, integrated with artificial intelligence-based modeling, enable mapping of cell state landscapes and transitional trajectories, transforming lineage plasticity from a descriptive phenomenon into a measurable and predictable property of tumor evolution.

Humans

Tumor microenvironment governs the prognostic landscape of immunotherapy for head and neck squamous cell carcinoma: A computational model-guided analysis.

Immune checkpoint inhibition (ICI) has emerged as a critical treatment strategy for squamous cell carcinoma of the head and neck (HNSCC) that halts the immune escape of the tumor cells. Increasing evidence suggests that the onset, progression, and lack of/no response of HNSCC to ICI are emergent properties arising from the interactions within the tumor microenvironment (TME). Deciphering how the diversity of cellular and molecular interactions leads to distinct HNSCC TME subtypes subsequently governing the ICI response remains largely unexplored. We developed a cellular-molecular model of the HNSCC TME that incorporates multiple cell types, cellular states, and transitions, and molecularly mediated paracrine interactions. Simulation across the selected parameter space of the HNSCC TME network shows that distinct mechanistic balances within the TME give rise to the five clinically observed TME subtypes such as immune/non-fibrotic, immune/fibrotic, fibrotic only and immune/fibrotic desert. We predict that the cancer-associated fibroblast, beyond a critical proliferation rate, drastically worsens the ICI response by hampering the accessibility of the CD8 + killer T cells to the tumor cells. Our analysis reveals that while an Interleukin-2 (IL-2) + ICI combination therapy may improve response in the immune desert scenario, Osteopontin (OPN) and Leukemia Inhibition Factor (LIF) knockout with ICI yields the best response in a fibro-dominated scenario. Further, we predict Interleukin-8 (IL-8), and lactate can serve as crucial biomarkers for ICI-resistant HNSCC phenotypes. Overall, we provide an integrated quantitative framework that explains a wide range of TME-mediated resistance mechanisms for HNSCC and predicts TME subtype-specific targets that can lead to an improved ICI outcome.

Tumor Microenvironment

Integrative analysis of single-cell sequencing identifies CD8+ TIM3+ CD101+ T cell-associated genes as prognostic biomarkers in breast cancer.

BACKGROUND: Breast cancer is a prevalent and deadly malignancy that significantly impacts women's quality of life and imposes financial burdens. Despite therapeutic advancements, tumour heterogeneity and frequent relapses remain major challenges. Accordingly, this study aimed to characterize immune features associated with CD8+ TIM3+ CD101+ T cells and develop a prognostic signature for breast cancer. METHODS: This study integrated single-cell and bulk transcriptomic datasets to characterize CD8+ TIM3+ CD101+ T cell (CCT)-related immune features and construct a prognostic signature in breast cancer. Single-cell RNA-seq data were sourced from the Gene Expression Omnibus (GEO) repository, and bulk transcriptomic data were from The Cancer Genome Atlas (TCGA) and GEO databases. Analytical methods included pseudo-time trajectory reconstruction (Monocle2), intercellular signalling analysis (CellChat), functional enrichment (ClusterProfiler), and immune profiling (ssGSEA). Prognostic modeling was conducted using least absolute shrinkage and selection operator (LASSO) Cox regression, with validation via Kaplan-Meier and time-dependent receiver operating characteristic (ROC) analyses. RESULTS: Single-cell analysis identified 17 clusters spanning seven cell types, including T cells, myeloid cells, and epithelial cells. T-cell sub-clustering revealed four subtypes. Pseudotime analysis suggested a potential state-transition relationship between CD8+ CD101- TIM3+ and CD8+ CD101+ TIM3+ T-cell states. A total of 121 differentially expressed genes were enriched in vital biological processes. An 11-gene prognostic model showed strong predictive power across cohorts. Single-cell T-cell reclustering identified a CD8+ CD101+ TIM3+ T-cell subpopulation, which was primarily characterized by the expression of markers such as CD101 and HAVCR2/TIM3. CONCLUSIONS: This study maps cellular heterogeneity and molecular networks in breast cancer, offering insights for targeted therapy and improved prognosis.

Breast invasive carcinoma

Transcription regulation of cell fate plasticity - from embryonic development to tissue regeneration.

Cell fate plasticity refers to the capacity of cells sharing the same genome to alter, reverse, or reconfigure their identity under physiological, pathological, or experimental conditions. This property underlies embryonic development, cellular reprogramming, and tissue regeneration, but becomes progressively restricted as lineage identity is stabilized. Embryonic development represents an intrinsic process of fate transitions, whereas reprogramming and regeneration reveal how differentiated cells can dedifferentiate or transdifferentiate under specific conditions. Across these contexts, plasticity is governed by multilayered regulatory networks involving transcription factors, epigenetic regulators, cofactors, and the core transcription machinery. Robust regulatory programs stabilize cell identity, whereas stochastic fluctuations in gene expression and chromatin state can prime cells for fate transitions, adding a tunable dimension to plasticity control. In this review, we synthesize recent advances in the regulation of cell fate plasticity across development, reprogramming, and regeneration, highlighting how transcription factors, epigenetic modifications, transcriptional cofactors, and core transcription machinery cooperate to control cell fate decisions and plasticity.

Animals

Spatially defined microenvironmental niches are associated with clinical outcome and tumor ecosystem diversity in head and neck cancer.

BACKGROUND: Head and neck squamous cell carcinoma (HNSCC) exhibits substantial biological heterogeneity that is not fully explained by human papillomavirus (HPV) status. The spatial organization of tumor, immune, and stromal cell populations and its relationship to clinical outcome remain incompletely understood. METHODS: We performed single-cell spatial transcriptomic and proteomic profiling of 44 primary HNSCC tumors, generating a spatial atlas of 19,471,501 cells across whole-slide tissue sections. Spatial niches and ecosystem states were identified through integrated computational analyses and evaluated for associations with tumor programs, clinicopathologic features, and patient outcomes. FINDINGS: HPV-negative tumors were enriched for fibroblast-rich, immune-poor niches associated with epithelial-mesenchymal transition and hypometabolic tumor programs, whereas HPV-positive tumors displayed more diverse immune, stromal, and vascular niche combinations and were enriched for immunogenic ecosystem states. Approximately 20% of HPV-positive tumors exhibited fibroblast-rich ecosystem architectures resembling HPV-negative disease and were associated with less favorable outcomes than other HPV-positive tumors of similar stage. In patient-derived co-culture models, extracellular matrix-associated fibroblasts were associated with epithelial-mesenchymal transition (EMT)-like tumor states, CD8+ T cell dysfunction, and chemotherapy resistance-associated phenotypes. CONCLUSIONS: Spatial ecosystem architecture is associated with clinically relevant heterogeneity beyond conventional HPV-based classification. Fibroblast-rich, immune-poor ecosystem states characterize a high-risk subset of HPV-positive tumors and may provide a framework for improved biological classification and risk stratification in HNSCC. FUNDING: This work was supported by the National Institutes of Health (R01CA291607 and R21CA267527-01) and the Feldstein Medical Foundation.

Humans

Toward AI Virtual Cells for Hepatology: Representation, Generation, Dynamics, and Intervention in Single-Cell Models.

``Single-cell and spatial atlases describe the healthy and diseased liver at high resolution, including lobular hepatocyte zonation, fibrotic macrophage-stellate niches, cholangiocyte reactions, immune remodeling, and hepatocellular carcinoma ecosystems. These maps show where cell states occur but do not, by themselves, predict whether liver injury will progress or how the liver will respond to an untested drug, toxicant, or genetic perturbation. In this review, we organize current approaches toward an AI Virtual Cell (AIVC) for the liver into three complementary modeling routes. Generative models represent cell states, dynamics and transport models infer state transitions, and pretrained or foundation models test whether learned representations transfer across donors, etiologies, disease stages, and platforms. Perturbation-response prediction serves as a cross-cutting assessment of whether these layers can predict responses to untested genetic, chemical, inflammatory, or metabolic interventions. Available evidence can be categorized as direct liver validation, liver-included benchmarks, general single-cell evidence, and conceptual applications. Published models demonstrate individual components, including atlas integration, inferred trajectories, transferable representations, and retrospective response programs. However, these models do not constitute a prospectively validated liver simulator. At minimum, evaluation should include donor-, etiology-, stage-, platform-, and perturbation-level hold-outs. Model performance should be reported using response direction, recovery of differentially expressed genes and rare states, and calibrated uncertainty. Claims about tissue- or function-level prediction additionally require independent spatial, histologic, metabolic, and functional readouts. Near-term use should prioritize experiment selection and hypothesis generation, whereas clinical decision support remains a longer-term objective.

AI Virtual Cell

Uncovering the signaling networks of disseminated glioblastoma cells in vivo with INSIGHT.

Dysregulation of intracellular signaling networks underpins cancer. Yet, resolving signaling networks within distinct or rare cell types in cancer in vivo has been unattainable. Here we develop INSIGHT by integrating cell sorting with mass spectrometry to enable quantitative phosphoproteomics and proteomics of discrete cell types from fixed tissues. Using INSIGHT, we map the signaling network within disseminating glioblastoma cells from patient-derived xenografts implanted in mice. Disseminating tumor cells undergo a proteome-wide shift from proliferative to mesenchymal, neural progenitor-like cell states. In parallel, signaling network and global kinase activity are rewired, transitioning from cell cycle-associated circuitries to those governing synaptic function, neuronal migration, and ion channel activity. Changes begin at the tumor margin and persist in distant brain parenchyma. Hornerin and phosphorylation of Ca²⁺-permeable GluA2 at Y876 were identified as mediators of glioblastoma progression. INSIGHT enables systems-level dissection of cell-type-specific signaling circuitries in vivo across wide range of biological systems.

Glioblastoma

Single-cell RNA sequencing of peripheral blood defines two immunological subtypes of Sjögren's disease distinguished by anti-SSA antibodies and aberrant B cell populations.

OBJECTIVES: Sjögren's disease (SjD) is a heterogeneous autoimmune disorder characterized by substantial clinical and molecular diversity. This heterogeneity raises key questions regarding the existence of distinct pathogenic mechanisms underlying disease subtypes. The objective of this study was to comprehensively characterize peripheral immune cell states associated with SjD and to identify features that could enable better patient stratification for targeted treatments. METHODS: We performed single-cell RNA sequencing with surface protein profiling on 1.5 million peripheral blood mononuclear cells (PBMCs) from 333 participants. Individuals were stratified by SjD diagnosis and anti-SSA status to enable comparative analyses between disease subgroups and controls. RESULTS: Our analysis identified two immunological endotypes of SjD, with SSA-positive participants exhibiting a dominant and persistent IFN-I signature that was also associated with altered immune cell composition. Transitional B cells were particularly affected, displaying altered developmental states, reduced BCR diversity, shorter CDR3 regions, and increased predicted interactions with activated immune cell populations, findings consistent with perturbations of early B-cell selection processes. By contrast, SSA-negative SjD participants exhibited limited transcriptional differences compared with symptomatic non-SjD controls, highlighting substantial biological heterogeneity within SjD. CONCLUSIONS: These findings support a two-disease model of SjD and highlight transitional B cells as both a key biomarker and a therapeutic target.

Journal Article