SEARCH · PubMed Health

Results for “single cell RNA sequencing”

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 91 records · Page 5Linked to original sources

PubMed2025

Colorectal Liver Metastasis Pathomics Model: Integrating Single-Cell and Spatial Transcriptome Analysis With Pathomics for Predicting Liver Metastasis in Colorectal Cancer.

The liver is the primary target organ for hematologic metastasis of colorectal cancer (CRC), and CRC liver metastasis (CRLM) often precludes radical resection, making it the leading cause of death in patients with CRC. To improve the identification and prediction of liver metastasis risk, we identified a cell type of liver metastasis--triggering malignant cells (LMTMCs) through integrating single-cell RNA sequencing and spatial transcriptome analysis. Multiomics cell communication analysis indicated that the interaction between fibroblasts and LMTMCs through the COL1A1-CD44/SDC4 and LAMA4-CD44 signaling axes could promote CRLM. By applying the one-class logistic regression algorithm, we developed a CRLM scoring system in the bulk RNA-sequencing data according to the abundance of LMTMCs in each individual. Using the grouping labels derived from the CRLM scoring system in the bulk data and the corresponding whole-slide images without any manual annotations at the region or pixel level, processed via slide-level weakly supervised learning, a deep-learning model based on the ResNet18 architecture, called Colorectal Liver Metastasis Pathomics Model, was developed to predict the risk of liver metastasis in patients with CRC. The Colorectal Liver Metastasis Pathomics Model achieved an area under the curve of 0.84 at the internal test set of The Cancer Genome Atlas-CRC histology images. In the external independent validation sets, namely the Affiliated Hospital of Southwest Medical University and the Affiliated Traditional Chinese Medicine Hospital of Southwest Medical University cohorts, the areas under the curve were 0.89 and 0.72, respectively, indicating effective classification performances. This study provided new insights and tools for the early identification of CRLM and demonstrated the potential of combining multiomics with deep learning-based pathomics in cancer research.

Humans↗

PubMed2026

Single-cell transcriptomics reveals heterogeneous stress responses and Mg2+-mediated survival mechanisms in Lactobacillus delbrueckii subsp. bulgaricus during freeze-drying and storage.

Maintaining the viability of lactic acid bacteria during dehydration and subsequent storage remains a significant challenge. Here, we employed single-cell RNA sequencing to reveal the heterogeneous stress responses of Lactobacillus delbrueckii subsp. bulgaricus, identifying seven distinct transcriptional clusters across the liquid culture, freeze-drying, and storage phases. The dominant clusters in the freeze-drying and storage were not completely consistent, showing significant functional differentiation. Genomic stability may be important for survival during freeze-drying and storage, while intracellular energy homeostasis appears important for viability during storage. The magnesium transporter mgtB was highly expressed in clusters tolerant to freeze-drying and storage, suggesting a critical role for Mg2+ homeostasis. Further experimental validation confirmed that Mg2+ treatment significantly bolstered stress resistance, increasing immediate post-freeze-drying survival by over 2-fold (up to 92.90%) and post-storage survival by over 5-fold (up to 5.98%). Proteomic data indicated that Mg2+ supplementation correlated with the maintenance of several biological functions potentially relevant to bacterial survival during freeze-drying and storage, including DNA repair, translation, and central carbon metabolism. These findings provide a map of microbial stress resistance through population heterogeneity and offer a potential strategy that may be adapted for enhancing the stability of other industrial lactic acid bacteria products.

Freeze Drying↗

PubMed2026

scGPA: an LLM-assisted workflow for directional virtual gene perturbation analysis from single-cell transcriptomes.

BACKGROUND: Existing virtual perturbation methods can often infer directional changes by comparing predicted post-perturbation expression profiles with control cells. However, workflows that directly return direction-specific downstream candidate genes together with confidence scores, evidence support and interpretable summaries remain limited. We developed scGPA, an LLM-assisted workflow system for directional single-cell virtual gene perturbation analysis. METHODS: scGPA starts from raw single-cell RNA sequencing data and performs quality control, normalization, dimensionality reduction, clustering and cell-group selection. It then constructs cell-group-specific wild-type regulatory networks using repeated subsampling, principal component regression (PCR)/Ridge-based network inference and CP tensor denoising. Based on these networks, scGPA simulates dose-aware virtual knockdown of the target gene and applies signed perturbation propagation to estimate the magnitude and direction of downstream transcriptional responses. LLM assistance is used for marker-based cell-type annotation, evidence-guided candidate prioritization and user-facing biological summarization. RESULTS: We benchmarked scGPA across five public Perturb-seq datasets and compared its performance with GEARS, scGPT and a random baseline. The overall correct prediction rate of scGPA was 23.0%, exceeding those of GEARS (20.7%), scGPT (15.1%) and the random baseline (13.6%). These results indicate that scGPA achieved a higher correct prediction rate than the two comparator models and the random baseline. We subsequently evaluated scGPA using a public osteosarcoma single-cell dataset and performed qRT-PCR validation in 143B osteosarcoma cells. Among genes with significant experimental changes, scGPA achieved a directional concordance of 76.9%. When all tested downstream genes were counted, 37.0% were directionally correct, 51.9% showed no significant change and 11.1% changed in the opposite direction. CONCLUSIONS: scGPA provides a practical workflow system for predicting and prioritizing direction-specific downstream transcriptional responses after target-gene perturbation. By integrating single-cell regulatory network inference, signed virtual perturbation and LLM-assisted interpretation, scGPA supports target-gene function inference and downstream mechanistic investigation from single-cell transcriptomic data.

Single-Cell Gene Expression Analysis↗

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↗

PubMed2025

Differentiation latency and dormancy signatures define fetal liver hematopoietic stem cells at single-cell resolution.

Decoding the mechanisms governing the self-renewal of hematopoietic stem cells (HSCs) during their expansion in the fetal liver (FL) could unlock novel therapeutic strategies to expand transplantable HSCs, a long-standing challenge. To explore intrinsic and extrinsic regulation of FL-HSC self-renewal at single-cell resolution, we engineered a culture platform replicating the FL endothelial niche that supports the amplification of serially engraftable HSCs. Leveraging this platform together with single-cell index flow cytometry, live imaging, transplantation assays, and single-cell RNA sequencing, we demonstrate that differentiation latency, cell-division symmetry, and transcriptional signatures of biosynthetic dormancy are distinguishing properties of rare FL-HSCs capable of serial multilineage hematopoietic reconstitution. Our findings support a paradigm in which intrinsic programs and niche-derived signals together facilitate the symmetric self-renewal of FL-HSCs while delaying their active participation in hematopoiesis. Our study also provides a resource for future investigations into intrinsic and extrinsic signaling pathways governing FL-HSC self-renewal.

Hematopoietic Stem Cells↗

PubMed2026

Monocytes Defined by Platelet Interactions and Oxidative Stress Signaling Underlie HIV-Associated Atherosclerosis.

BACKGROUND: Monocytes contribute to atherosclerosis by migrating into inflamed endothelium and differentiating into lipid-laden macrophages. In people living with HIV, chronic inflammation increases atherosclerosis risk, yet the role of specific monocyte subsets remains unclear. We investigated how distinct monocyte populations contribute to vascular pathology in early HIV-associated atherosclerosis. METHODS: We profiled 123 965 circulating monocytes using single-cell RNA sequencing and integrated plasma microparticle proteomics in 32 individuals stratified by HIV and atherosclerosis status. Supervised learning identified cluster- and disease-specific signatures, validated by platelet-monocyte cocultures, reverse transcription-quantitative polymerase chain reaction, bulk RNA sequencing, flow cytometry, and ELISA. RESULTS: Seven monocyte clusters were identified, including a subset characterized by platelet-monocyte complexes. Bulk RNA sequencing of platelet-monocyte cocultures revealed platelet-driven upregulation of genes involved in inflammation, lipid metabolism, oxidative stress, and endothelial adhesion. Platelet-monocyte complex-derived macrophages secreted higher levels of TGF-β (transforming growth factor-β) and IL-10 (interleukin-10), displayed decreased CD14 and increased CD80/CD86 while retaining CD36, and promoted endothelial-to-mesenchymal transition (decreased expression of CDH5 and PECAM1; increased expression of S100A4 and markers of vascular inflammation (ICAM1), VCAM), and IL-6 (interleukin-6), and CCL2 (C-C motif chemokine ligand 2) secretion). Additionally, CD14+ monocytes from HIV+ atherosclerosis-negative and HIV+ atherosclerosis-positive groups showed enhanced ROS-NRF2 (reactive oxygen species-nuclear factor erythroid 2-related factor 2) pathway activities, supported by increased basal and H2O2-induced p90RSK phosphorylation, indicating oxidative stress priming. CONCLUSIONS: Two monocyte clusters contribute independently to vascular immune dysregulation in people living with HIV. Platelet-monocyte complex-derived macrophages promote endothelial dysfunction while adopting a profibrotic cytokine profile. CD14+ monocytes show heightened oxidative signaling and stress responses, consistent with vascular activation. Together, these mechanisms may accelerate atherosclerosis development in HIV, even in the absence of traditional cardiovascular risk factors.

Humans↗

PubMed2026

Single-cell analysis of the human retina reveals stage-linked microglial states and neural-immune circuit rewiring in diabetic retinopathy.

Diabetic retinopathy (DR) is a major cause of vision loss worldwide. Here, we conduct single-cell RNA sequencing of twenty human retina samples (from living and post-mortem donors) across non-diabetic, diabetic, and DR states to create a comprehensive transcriptomic atlas. We identify two stable microglial populations-homeostatic and inflammatory-that exist along a functional continuum, plus a neutrophil cluster within C1QA+ myeloid cells with dynamic transitions occurring throughout disease progression. Module-level analysis reveals divergent transcriptional trajectories: homeostatic microglia maintain energetic programs while selectively upregulating stress elements, whereas inflammatory microglia layer additional pro-inflammatory programs onto preserved biosynthetic foundations. Eleven co-expression modules organize into two major axes: an inflammatory-stress axis, and a regulatory/metabolic-motility axis, with a stable translation module persisting across disease stages. Cell communication analysis further highlights sophisticated neural-immune interactions, particularly between photoreceptors and microglia. Our findings provide insights into the complex cellular dynamics of DR progression and suggest potential therapeutic targets for early intervention.

Humans↗

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↗

PubMed2026

MX1+ effector T cells hyperactivation at the maternal-fetal interface in unexplained recurrent pregnancy loss.

BACKGROUND: Immune tolerance breakdown at the maternal-fetal interface is implicated in unexplained recurrent pregnancy loss (URPL), but the interplay between T cell hyperactivation and dendritic cells (DCs)-mediated signaling remains poorly defined. METHODS: First-trimester decidual tissues from 5 healthy controls and 6 URPL patients underwent single-cell RNA sequencing (scRNA-seq, 10× Genomics). Computational analyses included clustering (Seurat), trajectory inference (scTour), intercellular communication (CellChat) and metabolic pathway enrichment (Gene Ontology and scMetabolism). Flow cytometry was performed from 11 patients and 11 healthy controls. Spatial validation was performed via multiplex immunohistochemistry and immunohistochemistry on 12 additional controls and 12 URPL cases. Statistical significance was assessed using Student’s t-test. RESULTS: URPL decidua exhibited marked CD3+ T cells and MX1+effector T (Tem) cells infiltration and activation. Flow cytometry analysis confirmed a significant decidua-specific upregulation of T cell activation markers CD25 and CD69 specifically on the MX1+Tem subset in URPL patients compared to controls. MX1+Tem cell subset demonstrated interferon hyperactivation, proliferative hyperactivity and lipid-biased immunometabolism. Pseudotemporal analysis positioned MX1+ Tem cells between classical Tem and exhausted T cell states, suggesting progressive differentiation. CellChat identified DCs as key regulators of MX1+ Tem expansion via aberrant ICOSL signaling, validated by spatial co-localization of ICOSL+ DCs and MX1+ Tem cells in URPL tissues. CONCLUSION: Our findings demonstrate that the aberrant activation and proliferation of MX1+Tem cells as a key immunological feature associated with URPL patients.

Humans↗

PubMed2026

Single-cell transcriptomic landscape of the southern green stink bug (Nezara viridula) midgut.

BACKGROUND: The southern green stink bug (SGSB), Nezara viridula, is a globally distributed hemipteran pest that damages many economically important crops. Its midgut supports digestion, defense, symbiosis, and interactions with orally delivered control agents, yet the cellular composition of this tissue remains poorly characterized. We therefore developed a single-cell transcriptomic atlas of the N. viridula midgut. RESULTS: Single-cell RNA sequencing of two biological replicates yielded a quality-filtered data set of 13,763 cells. Unsupervised clustering identified 12 transcriptionally distinct populations with putative annotations, including a stem cell/enteroblast (SC/EB)-like population, seven enterocyte-related populations, goblet-like cells, enteroendocrine cells, visceral muscle cells, and an extracellular-matrix-associated epithelial population. Enterocyte-related populations accounted for more than 77% of recovered cells. Putative annotations were assigned primarily from marker gene enrichment and homology to markers reported in other insects. Gene Ontology and Kyoto Encyclopedia of Genes and Genomes analyses identified population-associated functional enrichment patterns, and pseudotime analysis suggested transcriptional relationships between the SC/EB-like population and several enterocyte- and secretory-associated populations without establishing developmental lineages. Immune- and defense-associated transcripts were preferentially enriched in the pEC2 population, and genes associated with symbiont recognition, insecticide action, xenobiotic transport, and orally delivered double-stranded RNA showed population-biased expression. Descriptive comparisons with published insect midgut data sets identified shared and data-set-specific patterns among annotated populations. CONCLUSION: This atlas provides the first single-cell transcriptomic resource for a stink bug midgut and establishes a descriptive cellular framework for SGSB midgut biology. The dataset prioritizes candidate genes and cell populations for future spatial validation, functional testing, and studies of hemipteran midgut physiology, symbiosis, immunity, and pest-management-relevant traits. © 2026 Society of Chemical Industry.

Nezara viridula↗

PubMed2026

Identification and evaluation of glutamine-related gene characteristics based on multi-omics to predict the prognosis of patients with colorectal cancer.

BACKGROUND: Colorectal cancer (CRC), a prevalent malignancy of the gastrointestinal tract, ranks among the leading causes of cancer-related morbidity and mortality. Its clinical course is marked by high fatality and poor prognosis. Elucidating the mechanisms underlying CRC initiation and recurrence is therefore critical for identifying novel therapeutic targets. METHODS: This study incorporated two datasets, TCGA-CRC and GSE17537. A total of 84 glutamine metabolism-related genes (GMRGs) were identified, and differential expression analysis was conducted using the TCGA-CRC dataset. Weighted Gene Co-expression Network Analysis (WGCNA) was applied to determine gene modules most strongly associated with GMRG scores. Single-cell RNA sequencing (scRNA-seq) was utilized to characterize key cellular clusters and to identify differentially expressed genes (DEGs) between high and low glutamine metabolism (GM) groups. Overlapping GMRGs were visualized using the ggVennDiagram package in R. A CRC risk prediction model was developed through Cox proportional hazards and LASSO regression analyses, with performance evaluated by ROC curves. Cell type enrichment across 64 immune and stromal populations was assessed via xCell, and intergroup differences were tested using the Wilcoxon rank-sum test. TIDE scores were used to estimate immunotherapy responsiveness, while oncoPredict facilitated drug sensitivity profiling. PCOLCE2 expression in CRC was validated by RT-qPCR and Western blotting. Its functional role was examined through CCK-8 assays, invasion and migration tests, flow cytometry, and glutamate quantification. RESULTS: ScRNA-seq analysis identified two key cell populations and 437 DEGs associated with GM status. WGCNA pinpointed the MEgreen module as most significantly correlated with GMRG scores, encompassing 1075 genes. Integration of DEGs, module genes, and GM-related DEGs yielded 60 candidate genes for downstream analysis. A GMRG-based prognostic model comprising six genes (SRPX, CXCL1, GPX3, PCOLCE2, CLU, SEMA3E) demonstrated strong predictive performance. Prognostic gene expression correlated with immune and stromal infiltration patterns, as indicated by Spearman correlation analysis. The high-risk group exhibited diminished predicted response to immunotherapy (TIDE scores). Drug sensitivity analysis identified four compounds—Dasatinib-51, WH-4-023-56, TWS-119-366, and LDN-193189-478—with elevated efficacy in high-risk CRC cases. PCOLCE2 expression was significantly reduced in CRC tissues. Functional assays revealed that PCOLCE2 knockdown did not substantially affect cell proliferation but significantly impaired invasion and migration in CRC cells, increased apoptosis, and suppressed both glutamine uptake and glutamate production—highlighting its oncogenic role. CONCLUSION: Six GMRGs—SRPX, CXCL1, GPX3, PCOLCE2, CLU, and SEMA3E—were identified as key components of a robust prognostic model for CRC. These findings offer valuable insights into CRC pathogenesis and potential therapeutic strategies. Notably, this study provides the first evidence implicating PCOLCE2 as a tumor-promoting factor in CRC.

Glutamine↗

PubMed2026

Integrated single-cell and bulk transcriptomic analysis identifies a novel senescent fibroblast subtype associated with poor prognosis in acral melanoma.

BACKGROUND: Acral melanoma (AM) exhibits significant intratumoral heterogeneity, but its tumor microenvironment (TME) and immune regulation remain unclear. This study aims to dissect TME heterogeneity and establish a prognostic model based on key cell subpopulations. METHODS: We collected AM single-cell RNA sequencing (scRNA-seq) and bulk RNA-seq data from the Gene Expression Omnibus (GEO) and the Cancer Genome Atlas (TCGA). Unsupervised clustering, CellChat, and Scissor analysis were performed to characterize cellular heterogeneity, cell-cell communication, and prognosis-related cell subpopulations. Kaplan-Meier analysis was used to assess the prognostic value of key genes, which were further validated by multiplex immunohistochemistry (mIHC). RESULTS: In AM, Mel_C2, C7, and C9 with high SEMA6A and KIT expression were strongly linked to poor prognosis. We further identified a senescent fibroblast subpopulation (sCAF_CDKN2A) characterized by high fibroblast senescence signature (FSS) scores. Integrating Scissor analysis of fibroblast subtypes with bulk prognostic data, we identified COL3A1, VCAN, and KIT as prognosis-associated genes upregulated in poor-outcome-related fibroblast subsets. Cell-cell communication analysis revealed that sCAF_CDKN2A engages in an immunosuppressive network, interacting with regulatory T cells (Tregs) via MIF signaling and receiving signals from exhausted CD8+ T cells through PPIA-BSG interactions. Using transcription factor expression patterns from these fibroblast subtypes, we constructed a prognostic model that effectively stratified patients into distinct risk groups with significant differences in overall survival (OS). mIHC confirmed significantly higher protein levels of SEMA6A and COL3A1 in tumor tissues compared to matched normal tissues. CONCLUSIONS: We established a novel prognostic model for AM and identified sCAF_CDKN2A as an immunosuppressive senescent fibroblast subpopulation driving poor prognosis.

Acral melanoma↗

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↗

PubMed2024

Quantification of escape from X chromosome inactivation with single-cell omics data reveals heterogeneity across cell types and tissues.

Several X-linked genes escape from X chromosome inactivation (XCI), while differences in escape across cell types and tissues are still poorly characterized. Here, we developed scLinaX for directly quantifying relative gene expression from the inactivated X chromosome with droplet-based single-cell RNA sequencing (scRNA-seq) data. The scLinaX and differentially expressed gene analyses with large-scale blood scRNA-seq datasets consistently identified the stronger escape in lymphocytes than in myeloid cells. An extension of scLinaX to a 10x multiome dataset (scLinaX-multi) suggested a stronger escape in lymphocytes than in myeloid cells at the chromatin-accessibility level. The scLinaX analysis of human multiple-organ scRNA-seq datasets also identified the relatively strong degree of escape from XCI in lymphoid tissues and lymphocytes. Finally, effect size comparisons of genome-wide association studies between sexes suggested the underlying impact of escape on the genotype-phenotype association. Overall, scLinaX and the quantified escape catalog identified the heterogeneity of escape across cell types and tissues.

X Chromosome Inactivation↗

PubMed2026

FIERCE: reconstructing dynamic trajectories from the differentiation potency of single cells.

MOTIVATION: Since the introduction of single-cell RNA sequencing (scRNA-seq), numerous computational approaches have been developed to reconstruct dynamic cellular processes from static transcriptional profiles. These methods order cells along continuous trajectories by assessing their similarity in the gene-expression space. However, they rely on several assumptions, such as prior knowledge of the structure and directionality of the expected genealogy. These assumptions can limit their application to complex cellular systems with poorly understood developmental paths. RESULTS: To address this challenge, we introduce FIERCE (Framework for InfERence of the veloCity of Entropy), a novel computational pipeline designed to predict the changes in the differentiation potency of single cells during dynamic processes. Through a fully unsupervised approach, FIERCE enables the inference of cell lineages directly on the differentiation landscape of the biological system, thus eliminating the need for prior specification of developmental parameters. We demonstrate the efficacy of FIERCE by reconstructing three well-known mouse differentiation systems and by quantifying its accuracy on simulated data. AVAILABILITY AND IMPLEMENTATION: The FIERCE R package is available on GitHub at https://github.com/bicciatolab/FIERCE.

Cell Differentiation↗

PubMed2026

Protocol to identify genes required for cardiomyocyte development using Perturb-Seq.

While Perturb-Seq combines CRISPR-based screening with single-cell RNA sequencing (scRNA-seq), large-scale experiments are costly and its application during development is complicated by differentiation heterogeneity. Here, we present a protocol to identify genes required for cardiomyocyte development using Perturb-Seq. We describe steps for sgRNA (single guide RNA) library cloning and infection, cardiomyocyte differentiation, cell hashing, super loading, and scRNA-seq. We then detail procedures for sequencing, mapping, and data analysis. For complete details on the use and execution of this protocol, please refer to Sivakumar et al.1.

CRISPR↗

PubMed2025

Causal Inference for Genomic Data with Multiple Heterogeneous Outcomes.

With the evolution of single-cell RNA sequencing techniques into a standard approach in genomics, it has become possible to conduct cohort-level causal inferences based on single-cell-level measurements. However, the individual gene expression levels of interest are not directly observable; instead, only repeated proxy measurements from each individual's cells are available, providing a derived outcome to estimate the underlying outcome for each of many genes. In this paper, we propose a generic semiparametric inference framework for doubly robust estimation with multiple derived outcomes, which also encompasses the usual setting of multiple outcomes when the response of each unit is available. To reliably quantify the causal effects of heterogeneous outcomes, we specialize the analysis to standardized average treatment effects and quantile treatment effects. Through this, we demonstrate the use of the semiparametric inferential results for doubly robust estimators derived from both Von Mises expansions and estimating equations. A multiple testing procedure based on Gaussian multiplier bootstrap is tailored for doubly robust estimators to control the false discovery exceedance rate. Applications in single-cell CRISPR perturbation analysis and individual-level differential expression analysis demonstrate the utility of the proposed methods and offer insights into the usage of different estimands for causal inference in genomics.

Derived outcomes↗

PubMed2026

FZD5 drives macrophage-mediated immunomodulation and predicts prognosis in glioma: evidence from single-cell sequencing.

BACKGROUND: Gliomas are highly malignant brain tumors characterized by an immunosuppressive microenvironment, which limits therapeutic efficacy and contributes to poor clinical outcomes. The WNT/β-catenin signaling pathway is critically involved in tumor progression, and FZD5, a key receptor within this pathway, may participate in immune regulation. However, its specific role and underlying mechanisms in glioma remain unclear. METHODS: RNA-seq and microarray datasets from the Chinese Glioma Genome Atlas (CGGA) and The Cancer Genome Atlas (TCGA), together with single-cell RNA sequencing (scRNA-seq) datasets from GEO, were comprehensively analyzed. The Seurat package was used to identify macrophage-related clusters and mitophagy-associated pathways. Cox and LASSO regression analyses, along with a prognostic nomogram, were applied to evaluate the prognostic significance of FZD5. Immune infiltration, functional enrichment, and immunotherapy response analyses were conducted, followed by validation using spatial transcriptomics, immunohistochemistry, and in vitro assays. RESULTS: In bulk glioma transcriptomes, FZD5 emerged as an independent predictor of poor prognosis. Crucially, single-cell and spatial analyses revealed that the biologically significant FZD5 signal originated predominantly within tumor-associated macrophages (TAMs), where it colocalized with the M2 marker CD163. Consistently, elevated FZD5 levels correlated with increased myeloid infiltration and an immunosuppressive tumor microenvironment. Functionally, macrophage-expressed FZD5 was associated with mitophagy-related programs and promoted an M2-skewed phenotype, thereby enhancing glioma cell proliferation, migration, and invasion via macrophage-glioma crosstalk. CONCLUSION: FZD5 is a TAM-enriched marker in glioma tissues and a potential regulator of macrophage-associated immunosuppressive programs, supporting its utility as a prognostic biomarker and a candidate target for microenvironment-oriented interventions in glioma.

Humans↗