SEARCH · PubMed Health

Results for “Single cell transcriptomics”

Explore indexed PubMed citations for clinical trials, systematic reviews and public health research. Read source abstracts and follow each citation to its original PubMed record.

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

At least 37 records · Page 2Linked to original sources

PubMed2026

From wild to domestic: Single-cell transcriptomic perspectives on hippocampal regulation and evolution.

How domestication shapes brain evolution remains an open question. In this study, we integrated single-nucleus RNA sequencing (snRNA-seq), population genomics, and machine learning to investigate the hippocampal evolution under domestication. Across-species comparisons revealed that hippocampal cell type profiles are largely conserved across vertebrate species, while supporting the presence of adult hippocampal neurogenesis in birds. We further found that domestication and selective breeding likely influence the cellular composition and molecular regulation of the hippocampus. Our findings provide cellular evidence supporting the hypothesis that domestication affects adult hippocampal neurogenesis. Additionally, we showed that genes associated with neural progenitor cells (NPC) states and cell-marker programs are enriched for signatures of selection. Many of these genes function as regulators of neurogenesis and pathways mediating stress and fear reduction. Specifically, we identified selection at the FKBP5 promoter that may influence its expression in the NPC lineage, potentially contributing to stress-response regulation during domestication. Collectively, these results suggest that domestication is associated with hippocampal remodeling as part of an adaptive response to human-managed environments. This study provides a cellular and genetic perspective on how domestication reshapes the brain and offers a basis for further investigation into the mechanisms of neural evolution within the context of microevolution.

Animals↗

PubMed2024

Genetic association between epilepsy and gliomas: Insights from Mendelian randomization and single-cell transcriptomic analyses.

BACKGROUND: Seizures are prevalent in glioma patients, especially in those with low-grade gliomas. The interaction between gliomas and epilepsy involves complex biological mechanisms that are not fully understood. METHODS: We collected Genome-Wide Association Study data for epilepsy and gliomas, performed differential expression analysis, and conducted Gene Ontology (GO) enrichment analysis on the identified genes. Single-cell RNA sequencing data (scRNA-seq) from GSE221534 dataset in Gene Expression Omnibus (GEO) were used to analyze cell-cell interactions within glioma samples from patients with and without epilepsy. RESULTS: Mendelian Randomization (MR) analysis revealed significant associations between genetic variants related to epilepsy and glioma risk, suggesting a potential causal relationship, especially in astrocytomas. Differential expression analysis identified epilepsy-related genes that were significantly upregulated in astrocytoma tissues compared to normal brain tissues. GO enrichment analysis indicated that these genes are involved in critical biological processes such as neurogenesis and cellular signaling. The scRNA-seq analysis showed, compared to non-epileptic samples, glioma stem cells, microglia, and NK cells are increased in the core regions of astrocytomas in epileptic patients. Additionally, intercellular communication between tumor cells and other non-tumor cells is markedly enhanced in astrocytoma samples from epileptic patients. CONCLUSION: This study provides evidence of a genetic association between epilepsy and gliomas and elucidates the biological mechanisms through which epilepsy may influence glioma progression.

Humans↗

PubMed2026

REACTOR: REgulon Activity analysis and Comparison Tool for single-cell transcriptOmics Research.

SUMMARY: We introduce REACTOR, a computational tool designed to detect differential activity of transcriptional regulators and their target genes (regulons) in single-cell RNA-sequencing data. It expands the currently available framework for regulon analysis by introducing a robust statistical test to detect differential regulon activity between conditions, such as disease versus control, with multiple replicates. By contrasting different conditions, REACTOR enables identification of key condition- and cell type-specific regulons. To demonstrate the use of REACTOR, we illustrate its performance in a publicly available COVID-19 dataset. AVAILABILITY: REACTOR R-package together with an implementation vignette are available at https://www.github.com/elolab/REACTOR.

Regulon↗

PubMed2025

A Multi-omics Exploration Revealing SLIT2 as a Prime Therapeutic Target for Peripheral Facial Paralysis: Integrating Single-Cell Transcriptomics and Plasma Proteome Data.

Peripheral facial paralysis (PFP) is a common neurological disorder characterized by facial-nerve dysfunction. Identifying therapeutic targets and understanding the molecular and cellular mechanisms underlying PFP are crucial for developing effective treatment strategies. This study combined Mendelian randomization (MR) analysis and single-cell RNA sequencing (scRNA-seq) to explore potential therapeutic candidates and their roles in PFP pathophysiology. The MR analysis included 1925 publicly available plasma protein cis-heritability instruments. Instrumental variables were selected for MR analysis to identify plasma proteins associated with PFP, followed by colocalization analysis to evaluate shared genetic variants between the identified proteins and PFP. After the initial identification of plasma proteins associated with Bell's palsy using MR analysis, a rat model of facial-nerve injury was established to further dissect underlying mechanisms at cellular and molecular levels. Using scRNA-seq technology, we delved deeply into cellular Heterogeneity and dynamic changes in gene expression in the facial-nerve nucleus tissues under both injured and control conditions, thereby achieving a systematic study ranging from macroscopic genetic associations to microscopic cellular functions. Finally, expression patterns were preliminarily validated by performing in vitro immunofluorescence analysis on the facial-nerve nucleus samples of SD rats. The MR analysis results identified 30 plasma proteins significantly associated with PFP, with nine target genes showing differential expression in the scRNA-seq data. Colocalization analysis demonstrated that slit guidance Ligand 2 (SLIT2), semaphorin 4D (SEMA4D), EGF containing fibulin extracellular matrix protein 1 (EFEMP1), and sprouty related EVH1 domain containing 2 (SPRED2) shared causal variants with PFP. SLIT2 was highly expressed in the microglia and inhibitory neurons in the experimental group, whereas SEMA4D showed elevated expression across multiple glial cell types in the same group. In contrast, EFEMP1 and SPRED2 showed distinct expression patterns in fibroblasts and oligodendrocytes. The role of SLIT2 has been previously well-documented in many central nervous system diseases. However, for the first time, this study detected SLIT2 alteration after facial-nerve injury. Altered intercellular signaling, particularly enhanced SLIT2-ROBO signaling between neurons and glial cells, was observed in the PFP group. Pseudotime analysis revealed dynamic SLIT2 expression during microglia and inhibitory neuron differentiation, mirroring changes in ROBO1 expression. Immunofluorescence analysis of rat facial-nerve nucleus samples verified that SLIT2 protein levels were significantly increased in the facial-nerve nuclei of injured samples. In conclusion, despite the fact that this study is primarily founded on animal models and despite notable differences existing between animals and humans in terms of the facial motor nucleus, this study successfully identified SLIT2 as potential therapeutic targets for PFP. The SLIT2-ROBO axis stands out as a particularly promising candidate. SLIT2 may play a role in modulating neuroimmune interactions and promoting nerve repair. These findings provide a foundation for future clinical studies and targeted interventions to enhance recovery from PFP. Future research should focus on human sample validation to enhance clinical translation.

Animals↗

PubMed2026

Single-cell transcriptomics reveals distinct microglial state remodeling associated with the (R)-nicotine/diosmetin combination and galantamine in LPS-challenged BV2 cells.

BACKGROUND: Microglial neuroinflammation contributes to the progression of neurodegenerative diseases, yet it remains challenging to attenuate inflammatory responses while preserving cellular function. The effects of (R)-nicotine, diosmetin, their combined administration, and galantamine on heterogeneous BV2 transcriptional states have not been compared at single-cell resolution. METHODS: LPS-stimulated BV2 microglial-like cells were treated with (R)-nicotine, diosmetin, their combination (DR), or galantamine. Single-cell RNA sequencing was performed with three biological replicates per group and integrated with RNA velocity and SCENIC regulon inference to characterize treatment-associated state redistribution, inferred local transcriptional directionality and regulon-activity patterns. Functional validation included CCK-8 metabolic activity assays, multiplex cytokine ELISA, BDNF/GDNF quantification, qPCR, and high-content immunofluorescence analysis of iNOS and Arg1 at single-cell resolution. RESULTS: LPS decreased the relative abundance of the Itgae+/Plk4+ cluster while increasing the Nmur1+/Limk2+ cluster and inflammatory effector programs. DR treatment suppressed pro-inflammatory cytokine release without significantly reducing CCK-8-assessed metabolic activity, increased the Itgae+/Plk4+ cluster proportion, and increased BDNF/GDNF relative to LPS. RNA velocity and SCENIC analyses suggested that DR and galantamine showed distinct transcriptional and regulatory patterns: DR attenuated Batf-associated inflammatory regulons and was associated with increased Atf3-linked stress-response activity, whereas galantamine preferentially engaged DNA repair and genome-maintenance programs. CONCLUSION: These findings indicate that the DR condition was associated with remodeling of LPS-challenged BV2 microglial-like states, attenuation of inflammatory programs, and increased neurotrophic outputs relative to LPS, without significantly reducing CCK-8-assessed metabolic activity. This study provides a single-cell characterization of distinct treatment-associated responses to (R)-nicotine, diosmetin, their combined administration, and galantamine.

Microglia↗

PubMed2026

A latent activated olfactory stem cell state revealed by single-cell transcriptomic and epigenomic profiling.

The olfactory epithelium is one of the few regions of the nervous system that sustains neurogenesis throughout life. Its experimental accessibility makes it especially tractable for studying molecular mechanisms that drive neural regeneration in response to injury. In this study, we used single-cell sequencing to identify transcriptional and epigenetic processes involved in determining olfactory epithelial stem cell fate during injury-induced regeneration. By combining gene expression and accessible chromatin profiles of individual lineage-traced olfactory stem cells, we identified transcriptional heterogeneity among activated stem cells at a stage when cell fates are being specified. We further identified a subset of resting cells that appears poised for activation, characterized by accessible chromatin around silent genes prior to their expression in response to injury. These results provide evidence for a latent activated stem cell state in which a subset of quiescent olfactory epithelial stem cells are epigenetically primed to support injury-induced regeneration.

Animals↗

PubMed2025

Unsupervised multiscale clustering of single-cell transcriptomes to identify hierarchical structures of cell subtypes.

BACKGROUND: Cell clustering is an essential step in uncovering cellular architectures in single-cell RNA sequencing (scRNA-seq) data. However, the existing cell clustering approaches are not well designed to dissect complex structures of cellular landscapes at a finer resolution. RESULTS: Here, we develop a multiscale clustering (MSC) approach to construct a sparse cell-cell correlation network for unsupervised identification of de novo cell types and subtypes across multiple resolutions. Based upon simulated silver- and gold-standard data as well as real scRNA-seq data in diseases, MSC demonstrates significantly improved performance compared to established benchmark methods and reveals a biologically meaningful cell hierarchy to facilitate the discovery of novel disease-associated cell subtypes and mechanisms. CONCLUSIONS: We present MSC as a new single-cell multiscale clustering framework as a powerful tool for advancing discoveries in disease-associated cell populations using single-cell sequencing data.

Single-Cell Analysis↗

PubMed2026

Dynamic evolution of chaperone-mediated autophagy is associated with tumor microenvironment remodeling and prognostic stratification in lung adenocarcinoma: insights from single-cell transcriptomics, ensemble machine learning, and experimental validation.

BACKGROUND: Lung adenocarcinoma (LUAD) shows prognostic heterogeneity, and tumor-node-metastasis (TNM) staging is limited for individualized management. Chaperone-mediated autophagy (CMA) maintains proteostasis, but its role during adenocarcinoma in situ (AIS)-minimally invasive adenocarcinoma (MIA)-invasive adenocarcinoma (IAC) progression remains unclear. METHODS: Single-cell RNA sequencing (scRNA-seq) data from GSE189357 and bulk transcriptomes from The Cancer Genome Atlas (TCGA)-LUAD and Gene Expression Omnibus (GEO) cohorts were integrated. CMA activity, cell-cell communication, weighted gene co-expression network analysis (WGCNA), tumor-normal differential expression, machine-learning survival modeling, tumor microenvironment (TME) features, drug sensitivity, and EPC1 function were analyzed. RESULTS: CMA-high tumor epithelial cells increased from AIS (58.1%) to MIA (65.7%) but declined in IAC (44.4%; p < 0.001). CMA-low cells preferentially received fibroblast-derived extracellular matrix cues. A CMA-negatively correlated module identified 69 core genes. Random survival forest (RSF) performed best among 117 machine-learning combinations (mean concordance index > 0.873). High-risk patients had worse survival across cohorts, and the risk score was independently associated with overall survival (hazard ratio = 16.013, 95% confidence interval: 9.579-26.768, p < 0.001). High-risk tumors showed proliferative activation and M0 macrophage enrichment, whereas low-risk tumors showed stronger immune-related signaling. EPC1 overexpression suppressed malignant phenotypes in A549 cells. CONCLUSION: CMA dynamics are associated with stromal and immune remodeling during LUAD progression. A CMA-based model provides robust prognostic stratification and may offer a basis for future TME-guided studies.

Chaperone-mediated autophagy↗

PubMed2025

Identification of a prognostic signature consisting of three macrophage-related genes for glioblastoma based on bulk and single-cell transcriptomes analyses.

BACKGROUND: Tumor-associated macrophages have been implicated in the progression and treatment resistance of glioblastoma (GBM). This study aimed to identify macrophage-related genes associated with prognosis and therapeutic response in GBM. MATERIALS AND METHODS: Bulk RNA-seq data from 533 patients with GBM were downloaded from the Cancer Genome Atlas (TCGA) and Chinese Glioma Genome Atlas (CGGA) databases. Bioinformatic tools were used to detect the co-expression gene modules associated with the infiltration of immune cells, identify a prognostic macrophage-related gene signature, and explore their association with sensitivity to chemotherapeutic drugs and immune checkpoint blockade. Single-cell RNA-seq data and multiplexed immunofluorescence were used to validate ISG20 expression (a member of the identified gene signature) in macrophages. RESULTS: We detected gene modules associated with macrophages and identified a signature consisting of three macrophage-related genes (ISG20, PARP12 and IFIT5) in the discovery set (TCGA-GBM, n = 159), and validated its prognostic value in the validation set (CGGA-GBM, n = 374). This gene signature demonstrated favorable accuracy in predicting prognosis and resistance of immuno- and chemo-therapy. The co-expression of ISG20 and PD-1 in macrophages was verified by single-cell RNA-seq data and multiplex immunofluorescence. CONCLUSIONS: This study presents a macrophage-related gene signature to predict prognosis and therapeutic response in GBM. ISG20, PARP12 and IFIT5 are interferon-stimulated genes, and further investigations may provide new insights into the interplay between macrophages and interferon signaling in GBM.

Humans↗

PubMed2025

A single-cell transcriptomic atlas of the pigtail macaque placenta in late gestation.

The placenta is a complex organ with multiple immune and non-immune cell types that promote fetal tolerance and facilitate the transfer of nutrients and oxygen. The nonhuman primate (NHP) is a key experimental model for studying human pregnancy complications, in part due to similarities in placental structure, which makes it essential to understand how single-cell populations compare across the human and NHP maternal-fetal interface. We constructed a single-cell RNA-Seq (scRNA-Seq) atlas of the placenta from the pigtail macaque ( Macaca nemestrina ) in the third trimester, comprising three different tissues at the maternal-fetal interface: the chorionic villi (placental disc), chorioamniotic membranes, and the maternal decidua. Each tissue was separately dissociated into single cells and processed through the 10X Genomics and Seurat pipeline, followed by aggregation, unsupervised clustering, and cluster annotation. Next, we determined the maternal-fetal origins of cell populations and analyzed single-cell RNA trajectory, Gene Ontology enrichment, and cell-cell communication. Single-cell populations in the pigtail macaque were strikingly similar in their identity and frequency to those found in the human placenta, including cells from trophoblast, stromal cell, immune, and macrophage lineages. An advantage of our approach was the deep sequencing of three tissues at the maternal-fetal interface, which yielded a rich diversity of common and rare single-cell populations. The third-trimester pigtail macaque single-cell atlas enables the identification of cellular subclusters analogous to those in humans and provides a powerful resource for understanding experimental perturbations on the NHP placenta.

Journal Article↗

PubMed2026

Single-Cell Transcriptome-Wide Mendelian Randomization and Colocalization Uncover Potential Immunocytes-Related Therapeutic Targets for Obesity.

Weight-loss treatment is crucial for individuals with obesity to prevent various complications. The role of Immune cells in obesity has been recently recognized, whereas its translation into therapy requires identifying key target genes. We performed Mendelian randomization (MR) analysis to assess causal relationships between expression quantitative trait loci (eQTL) of 14 immune cells and obesity-related traits (obesity, body mass index and body fat percentage), and validated the results in colocalization analysis. For the putative causal genes identified by the MR and colocalization analyses, we conducted pathway enrichment, differential expressed gene (DEG) analysis and search of druggable evidence, and utilized a Tier system to prioritize drug targets for obesity. MR and colocalization evidence was observed for 1630 genes associated with one or more obesity-related traits, mainly expressed in CD4+ naive/central memory T cells and enriched in antigen processing and presentation pathways. Forty-one genes showed causal relationship with all three outcomes, among which 19 genes have not been reported for obesity previously. DEG analysis using single-cell RNA sequencing data of blood or adipose tissue indicated that the differential expression of UBE2Z in monocytes, ZCCHC7 in T cells, and FNBP4 in B cells between lean and obese individuals were consistent with the MR results. By searching drug-gene interaction databases, we found targeted drugs for PYGB and PRUNE1, and PYGB was the top gene ranked in the Tier system. This study provides evidence for the involvement of immune cells in obesity, and the potential cell-specific, immune-related targets for obesity treatment.

Obesity↗

PubMed2026

Full-length single-cell spatial transcriptomics reveals spatial and cell-type-specific transcript isoforms in the primate brain.

The primate brain exhibits complex RNA alternative splicing heterogeneity crucial for functional complexity, yet systematic spatial isoform characterization has been lacking. We developed Fullscope-seq, a full-length single-molecule large field-of-view spatial transcriptomics sequencing method at single-cell resolution, based on programmed concatenation cDNA for multiple long-read sequencing platforms. Applying Fullscope-seq to the macaque brain, we uncovered thousands of genes exhibiting differential transcript usage (DTU) across cortical layers, cell types and brain regions. Fullscope-seq resolved hundreds of major isoform switches across distinct brain regions and identified DTUs between superficial and deep cortical layers. Cortical layer-specific DTUs showed cell-composition dependence, whereas regional DTUs were regulated according to both cellular composition and spatial contexts. These isoform variations showed substantial enrichment for neuropsychiatric disorder-associated genes and were conserved across platforms and species. Our study establishes a scalable framework for spatial isoform analysis and provides a resource for understanding transcriptomic diversity in complex tissues.

Animals↗

PubMed2026

HIF1A+CSF3R+ neutrophils-dominated hypoxic niche induced metabolic reprogramming for neoadjuvant therapy resistance in NSCLC.

BACKGROUND: Non-small cell lung cancer (NSCLC) is one of the frequently occurring cancers characterized by molecular heterogeneity and multiple immune cell infiltration patterns, which are associated with treatment sensitivity and resistance. However, the specific microenvironmental cells and their mechanisms that lead to treatment resistance in patients need to be explored in greater depth. METHODS: On the basis of patients receiving neoadjuvant therapy in our center, a multicenter, multicohort NSCLC spatial transcriptome, single-cell transcriptome, T-cell receptor repertoire sequencing, bulk RNA transcriptome, phosphorylated proteome, genome mutation, and clinical data were included for a comprehensive assessment of the therapeutic and prognostic impact of HIF1A+ CSF3R+ neutrophils in NSCLC. In vitro experiments validated the functional phenotype of HIF1A+ CSF3R+ neutrophils and co-localization interactions with other cellular subpopulations. Gradient boosting machine (GBM) constructed region of interest (ROI) models for evaluation. Computer-aided drug design (CADD) was used to predict targeted small molecule drugs, and in vivo mouse models were constructed to assess the effectiveness of the combination treatment regimen. RESULTS: Centered on HIF1A+ CSF3R+ neutrophils, recruited exhausted T cells and stromal cells form a hypoxic niche within the tumor region, which was enriched in non-response patients. ROI composed of these specific cellular subpopulations, associated with senescence and glycolysis, accurately predicting NSCLC progression, prognosis, and microenvironment composition. CADD analysis identified that platycodin-D2 specifically targeted CSF3R, reducing HIF1A expression and inhibiting neutrophil activity. Combining navitoclax, platycodin-D2 with anti-programmed cell death protein 1 (PD-1) significantly suppressed tumor proliferation and improved the immunosuppressive microenvironment. CONCLUSION: Our study emphasized the role of HIF1A+ CSF3R+ neutrophils in immunotherapeutic resistance of NSCLC, constructed a microenvironmental immune dysregulation network in a hypoxic ecological niche with HIF1A+ CSF3R+ neutrophils as the center. Platycodin-D2 specifically targeted HIF1A+ CSF3R+ neutrophils, enhancing the efficacy of anti-PD-1 therapy in NSCLC.

Humans↗

PubMed2026

Stage-specific ROMO1 in rheumatoid arthritis: predictive immune insights into the MIF pathway and HLA-DR/IL2RA axis via integrated GWAS, transcriptomic, single-cell, and spatial profiling.

Emerging evidence links reactive oxygen species modulator 1 (ROMO1), a key mitochondrial ROS regulator, to rheumatoid arthritis (RA) pathogenesis. However, its exact mechanism remains elusive given the conflicting evidence about its specific function. We used a four-level integrative framework combining multi-omics data and literature‑supported mechanistic inference. At the genetic level, Mendelian randomization (MR) was performed to explore potential causal relationships between ROMO1, IL2RA, HLA-DR, MIF, and RA risk, followed by differential expression analysis and machine learning-based feature selection to identify key mROS genes. The temporal expression dynamics of ROMO1 were assessed in RA progression. At the cellular and tissue levels, we integrated single-cell RNA sequencing and spatial transcriptomics to map cell-type-specific expression and synovial localization of ROMO1-related immune cells and pathways. Finally, our multi-omics findings were contextualized with literature-supported mechanistic inference. (1) MR results were consistent with a potential protective effect of ROMO1 on RA (OR = 0.52) and its potential regulation of risk factors IL2RA (OR = 0.46) and HLA-DR (OR = 0.40). Conversely, IL2RA (OR = 1.42), HLA-DR (OR = 1.88), and MIF (OR = 1.17) were positively associated with RA risk. Additionally, ROMO1 was identified as a top candidate diagnostic predictor with stage-specific dynamics: downregulated in the early but upregulated in the late/remission stages. (2) Single-cell RNA sequencing showed ROMO1's cell-specific expression in CD14+ HLA-DR+ CD74+ monocytes and CD4+ IL2RA+ T cells. Cell communication analysis further suggested that these cells may participate in MIF pathway regulation. Spatial transcriptomics subsequently identified that ROMO1-related cells localized to synovial pathological regions, with MIF pathway changes correlated with RA progression. (3) Finally, literature-supported mechanistic inference suggests that ROMO1 may modulate mROS levels to promote anti-inflammatory M2 macrophage polarization, which could theoretically contribute to reduced systemic inflammation and the alleviation of multi-organ decline in RA. This integrated multi-omics investigation, supported by literature-based mechanistic inference, suggests ROMO1 as a stage-dependent biomarker candidate and potential immune regulator in RA.

Humans↗

PubMed2026

A nucleolar stress gene signature enables quantitative scoring across multi-omics contexts.

The nucleolus is essential for ribosome biogenesis and cellular homeostasis, and its dysfunction can induce nucleolar stress, a process implicated in cancer and other diseases. However, nucleolar stress is commonly inferred from morphological changes or a limited set of functional assays, and quantitative approaches based on gene expression profiles remain lacking. Here, we integrate literature curation with multi-dataset screening to define a nucleolar stress gene signature and develop a nucleolar stress score (NuS) applicable to bulk transcriptomics, single-cell transcriptomics, proteomics, and spatial transcriptomics. Using this framework, we show in colorectal cancer models that oxaliplatin induces nucleolar stress, suppresses nascent rRNA synthesis, and activates p53 signaling, whereas these responses are attenuated in oxaliplatin-resistant cells. Combined with a ribosome biogenesis activity score (RiboSis), NuS captures related but distinct dimensions of nucleolar function and stratifies tumors into functional states associated with clinical outcomes. NuS-based analysis of perturbational transcriptomes further prioritizes compounds with putative nucleolar stress-inducing activity. Collectively, this study provides a quantitative framework for evaluating nucleolar stress and illustrates its applications in disease stratification and drug mechanism discovery.

Cell Nucleolus↗

PubMed2025

A pluripotent stem cell atlas of multilineage differentiation.

Human pluripotent stem cells offer a scalable platform to study genetic and signalling mechanisms governing cell lineage decisions during differentiation. Genome-wide and single-cell transcriptomics technologies likewise offer high-throughput analysis of heterogeneous cell differentiation states. While in vivo development has been extensively characterised using these technologies, there remains a need for comprehensive single-cell transcriptomic profiling of stem cell differentiation from pluripotency. Understanding gene expression changes governing differentiation in vitro is key to developing high fidelity differentiation protocols and understanding fundamental mechanisms of development. We generated a single-cell RNA sequencing time course to study the role of developmental signalling pathways on multilineage diversification from pluripotency in vitro. The combined dataset of over 60,000 cells spans cell types from a time course of differentiation across all germ layers, ranging from gastrulation cell states to progenitor and committed cell types. These data provide a diverse benchmarking reference point to compare against in vivo development and advance understanding of signalling regulation of differentiation, providing insights into protocol development, drug screening, and regenerative medicine applications.

Pluripotent Stem Cells↗

PubMed2024

A time-resolved single-cell roadmap of the logic driving anterior neural crest diversification from neural border to migration stages.

Neural crest cells exemplify cellular diversification from a multipotent progenitor population. However, the full sequence of early molecular choices orchestrating the emergence of neural crest heterogeneity from the embryonic ectoderm remains elusive. Gene-regulatory-networks (GRN) govern early development and cell specification toward definitive neural crest. Here, we combine ultradense single-cell transcriptomes with machine-learning and large-scale transcriptomic and epigenomic experimental validation of selected trajectories, to provide the general principles and highlight specific features of the GRN underlying neural crest fate diversification from induction to early migration stages using Xenopus frog embryos as a model. During gastrulation, a transient neural border zone state precedes the choice between neural crest and placodes which includes multiple converging gene programs. During neurulation, transcription factor connectome, and bifurcation analyses demonstrate the early emergence of neural crest fates at the neural plate stage, alongside an unbiased multipotent-like lineage persisting until epithelial-mesenchymal transition stage. We also decipher circuits driving cranial and vagal neural crest formation and provide a broadly applicable high-throughput validation strategy for investigating single-cell transcriptomes in vertebrate GRNs in development, evolution, and disease.

Animals↗

PubMed2025

An intermediate activation state primes Langerhans cell migration from the epidermis.

Langerhans cells (LCs) are a specialized subset of dendritic cells in the epidermis, forming a dense network that acts as a frontline defense through immune surveillance. Upon antigen uptake, LCs become activated and orchestrate subsequent immune responses by migrating to lymphatics. However, how transcriptional programs are regulated during activation and how LCs behave in vivo during this transition remain poorly understood. Here, we combine single-cell transcriptomic analysis and intravital imaging to reconstruct the activation trajectory of epidermal LCs. We present a high-resolution single-cell transcriptomic dataset of over 22,000 high-quality epidermal LCs in both homeostatic and injured conditions. We define specific LC subpopulations representing sequential activation stages, characterized at the level of pathways and transcription factors. Notably, we identify a distinct intermediate state that precedes their migration. Integrating our data with an external dataset from homeostatic and injured skin reveals that wound-specific, WNT-modulated fibroblasts are the primary source of C3, the central component of the complement cascade. Intravital imaging of C3-deficient mice demonstrated that C3 is essential for effective recruitment of activated LCs to wound sites. Together, our findings uncover a novel population of activated epidermal LCs and highlight complement signaling as a critical mediator of LC recruitment during skin injury.

Activation↗