SEARCH · Search PubMed

Results for “ScRNA-seq”

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.

SourceResearch connectionFrom yearThrough year

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

At least 73 records · Page 4Linked to original sources

PubMed2026

transFusion: a novel comprehensive platform for integration analysis of single-cell and spatial transcriptomics.

MOTIVATION: Understanding spatial organization, intercellular interactions, and regulatory networks within the spatial context of tissues is crucial for uncovering complex biological processes and disease mechanisms. Spatial transcriptomics technologies have revolutionized this field by enabling the spatially resolved profiling of gene expression. 10× Visium has emerged as the predominant spatial technology, but its low resolution and the complexity of integrating multimodal datasets present significant analytical challenges, particularly for researchers with limited computational and statistical expertise. Current spatial transcriptomics analysis platforms generally fall short of effectively integrating multimodal data and maximizing the utility of spatial information-such as uncovering complex cellular spatial dependencies, multimodal gradient patterns, and spatial coexpression of ligand-receptor pairs and regulatory networks related to disease or biological states-thereby limiting their ability to provide comprehensive end-to-end analytical workflows when analyzing 10× Visium data. RESULTS: To address these limitations, we developed transFusion, a novel, advanced web-based platform specializing in the most comprehensive and effective integration analysis of scRNA-seq and 10× Visium spatial transcriptomics data. transFusion offers 12 key functions, from basic visualization to advanced analyses, including intercellular dependency analysis, ligand-receptor coexpression identification and visualization, and spatial multimodal gradient variation patterns. Two case studies were used to demonstrate transFusion's capabilities in exploring tissue architecture, intercellular communication, dependency networks, and multimodal gradient variation patterns with minimal computational skills and statistical expertise. transFusion provides a flexible and powerful framework for multimodal data integration analysis. AVAILABILITY AND IMPLEMENTATION: transFusion is freely available at https://github.com/WQLin8/transFusion.

Spatial Transcriptomics↗

PubMed2026

Nested co-expression network analysis identifies compact gene clusters in a black box.

MOTIVATION: Digital analysis of biological systems requires methods capable of identifying both broad and nested gene modules reflecting complex biological processes. Existing transcriptomic methods often miss compact gene sets corresponding to subprocesses in specialized cell types, limiting insights into functional heterogeneity. RESULTS: We present Nested-WGCNA, a two-stage unsupervised network analysis algorithm designed to identify coarse-grained and fine-grained gene modules. Applied to bulk RNA-Seq data, Nested-WGCNA reveals stable modules reproducible across datasets. When validated against scRNA-Seq data, these modules correspond to both major and minor immune cell subtypes. Application to immunotherapy response datasets uncovers predictive and prognostic biomarkers, highlighting its utility in treatment stratification and biomarker discovery. AVAILABILITY: The NestedWGCNA source code and analysis pipeline are available on GitHub (https://github.com/ilyada/NestedWGCNA) and archived on Zenodo (https://doi.org/10.5281/zenodo.18959244).

Algorithms↗

PubMed2026

Integrated multi-omics strategies for identifying novel therapies in psoriasis.

MOTIVATION: Psoriasis is a chronic, immune-mediated disorder with an unmet need for effective treatments. To systematically prioritize therapeutic targets, we integrated proteome-wide Mendelian randomization (MR) with expression validation in blood/skin, genetic susceptibility analysis, differential gene expression (DGE) from bulk and single-cell RNA sequencing (scRNA-seq), colocalization, pathway enrichment, and protein-protein interaction analyses. RESULTS: Proteome-wide MR identified 29 candidate protein targets (Bonferroni-corrected), all replicated in independent datasets. Fifteen targets showed significant expression associations in blood or skin. Eleven proteins-UBLCP1, IL23A, ASF1A, RARRES2, ICAM1, PRSS53, ICAM5, GCA, IL2RA, DBI, and NFKB1-exhibited consistent directional effects with their genes. Genetic susceptibility analysis confirmed 20 target-specific polygenic scores for psoriasis and five for psoriatic arthritis. DGE analysis identified 13 targets in bulk and 13 in scRNA-seq-primarily in keratinocytes and immune cells-with IL2RA, COMP, and A2ML1 dysregulated across both. Colocalization analysis implicated shared causal variants for psoriasis in ASF1A, CD8A, CTF1, IL7R, MMP12, RARRES2, XCL2, DBI, IL23A, IL2RA, SGSH, and TIMD4. Enrichment analyses highlighted involvement in cytotoxicity, immune regulation, and JAK-STAT signaling. Eighteen targets interacted with approved anti-psoriasis drugs. Notably, drugs targeting IL2RA, IL7R, CTF1, ICAM1, MMP12, NFKB1, CD8A, DDX58, IL12A, SGSH, and FAP are approved or in trials for other diseases, suggesting repurposing potential. Our integrative multi-omics approach prioritized 29 high-confidence targets, including 13 novel candidates (RARRES2, ASF1A, CTF1, DBI, B3GNT2, CD8A, TIMD4, CRTAM, SGSH, XCL2, DAPK2, A2ML1, and FAP). Several high-priority targets-such as IL2RA, IL23, MMP12, RARRES2, IL7R, and ICAM1-were supported across analytical layers. These findings provide a robust foundation for psoriasis drug development. AVAILABILITY AND IMPLEMENTATION: The code used for the analyses in this manuscript has been archived in Zenodo at [DOI: 10.5281/zenodo.19692128].

Psoriasis↗

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

TMEM176B is co-expressed with TMEM176A and upregulated in peripheral blood monocytes of patients with primary Sjögren's syndrome.

OBJECTIVES: This study aims to determine the role of acid‑sensitive nonspecific cation channels transmembrane protein 176A (TMEM176A) and TMEM176B in autoimmune diseases, with a focus on primary Sjögren's Syndrome (pSS). METHODS: We examined the expression patterns of TMEM176A and TMEM176B across tissues and cells utilizing bulk RNA-seq and scRNA-seq datasets. Immunophenotyping analysis was performed by flow cytometry to compare CD62L expression between TMEM176B⁺ and TMEM176B⁻ monocytes. The proportion of TMEM176B+ cells in monocytes was interrogated in both pSS patients and healthy controls. Clinical correlations of TMEM176B with anti‑SSB antibody and complement C4 levels were also evaluated. RESULTS: TMEM176A and TMEM176B showed conserved co‑expression and were significantly upregulated in autoimmune diseases. TMEM176B⁺ monocytes displayed higher CD62L positivity rate than TMEM176B- monocytes. In pSS patients, the proportion of TMEM176B⁺ monocytes was elevated in total monocytes, classical monocytes (cMo) and intermediate monocytes (iMo). The proportion of TMEM176B+ monocytes positively correlated with anti‑SSB levels, while several TMEM176B-associated monocyte subset markers inversely correlated with C4. CONCLUSION: TMEM176A and TMEM176B are highly correlated. TMEM176B expression in monocytes is linked to pSS and may serve as a novel auxiliary diagnostic biomarker.

Humans↗

PubMed2026

Spatial proteomic mapping of the human and mouse retina using IBEX.

We generated a comparative spatial proteomic atlas of the human and mouse retina using a highly multiplexed immunohistochemistry technique called iterative bleaching extends multiplexity (IBEX). We refined the IBEX workflow by integrating an antibody dissociation option alongside chemical bleaching. This dual strategy enabled removal of the entire antibody complex, permitting the flexible use of antibodies from the same host species across iterative cycles. We coupled this workflow with super-resolution imaging via deconvolution and applied it to the retina of healthy humans and WT mice and the Crb1rd8 mouse model. We successfully imaged over 25 protein markers on human and mouse tissue sections, generating spatial atlases of the major retinal cell populations. Cross-species protein expression was compared to scRNA-seq datasets to identify protein and transcript disparities. Super-resolution IBEX delineated the ultrastructural features of the outer limiting membrane (OLM), identifying CD44 as a core structural component tightly colocalized with a highly organized F-actin belt within Müller glial endfeet. Using the Crb1rd8 mouse model, disruption of this complex was spatially associated with rosette formation and OLM structural failure. In summary, spatial proteomic atlases of the human and mouse retina were used to reveal insights into the arrangement of major retinal cell populations and OLM structure.

Animals↗

PubMed2025

From genetic causality to druggable targets: A multiomics framework identifies ZSCAN16 in gout pathogenesis.

ObjectiveGout is a prevalent form of inflammatory arthritis in which many patients respond suboptimally to current therapies. Drug development is hampered by a lack of genetically validated targets, leading to high clinical trial attrition. This study aimed to systematically identify and prioritize novel, druggable targets for gout via a multilayered genetic and functional genomics approach.MethodsWe performed two-sample Mendelian randomization (MR) using cis-expression quantitative trait locus (cis-eQTL) data and dual independent gout genome-wide association study (GWAS) cohorts (openGWAS and FinnGen). The candidate genes were subjected to a rigorous validation pipeline including Bayesian colocalization, phenome-wide association studies (PheWASs) to assess pleiotropy and on-target safety, and single-cell RNA sequencing (scRNA-seq) to delineate the cellular context. Molecular docking was used to evaluate the structural druggability of prioritized targets.ResultsMR analysis revealed 15 genes causally associated with gout. Colocalization analysis (PPH4 > 0.8) prioritized two targets: ZSCAN16 (risk-increasing, OR = 1.04, 95% CI [1.02-1.06]) and TRIM10 (protective, OR = 0.96, 95% CI [0.94-0.98]). Crucially, PheWAS revealed that ZSCAN16 is highly specific to gout, whereas TRIM10 exhibited extensive pleiotropy with hematological and cardiometabolic traits, indicating significant safety risks. Single-cell analysis provided orthogonal validation, demonstrating flare-specific upregulation of ZSCAN16 in cytotoxic T/NK cells. Molecular docking confirmed ZSCAN16 as a structurally druggable target, showing high-affinity binding with known compounds (e.g. digoxin, binding energy = -9.6 kcal/mol).ConclusionsOur study identifies ZSCAN16 as a high-potential, druggable therapeutic target for gout, highlighting its genetic influence on specific immune cell activities during acute flares. Conversely, TRIM10 was deprioritized owing to substantial pleiotropic liabilities and poor chemical tractability. These findings suggest that ZSCAN16 could play a crucial role in the pathogenesis of gout and may provide a valuable lead for future drug discovery efforts.

Humans↗

PubMed2026

Epigenetic and immunological alterations in umbilical cord blood of overweight/obese women with gestational diabetes mellitus: insights into DNA methylation signatures and immune cell dysregulation.

BACKGROUND: Gestational diabetes mellitus (GDM) is a common pregnancy complication associated with adverse maternal and neonatal outcomes. Epigenetic modifications may reflect intrauterine metabolic exposure and contribute to immune and metabolic alterations. This study aimed to explore DNA methylation profiles in umbilical cord blood from overweight and obese women with and without GDM. METHODS: Umbilical cord blood samples from 30 overweight/obese pregnant women (with and without GDM) were analyzed using the Illumina 850 K methylation array to identify differentially methylated positions (DMPs) and regions (DMRs). Gene Ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) enrichment analyses were performed to assess the functional relevance of methylation changes. Immune cell composition was estimated using deconvolution analysis and further examined in an independent single-cell RNA sequencing (scRNA-seq) cohort. Lasso regression was applied to identify CpG sites associated with GDM status and construct a preliminary methylation-based classification model. RESULTS: A total of 23,331 hypermethylated and 29,501 hypomethylated DMPs were identified between women with and without GDM, with hypomethylation predominating. Enrichment analyses indicated associations with neurodevelopmental pathways, metabolic processes, immune regulation, and epigenetic modification. Immune deconvolution analysis suggested reduced proportions of CD4+ T cells (p < 0.05) and a trend toward decreased NK cells in the GDM group, alongside increased CD8+ T cells and neutrophils. Seven CpG sites were selected for model construction and demonstrated strong discriminatory performance within this cohort. CONCLUSION: This exploratory study identifies distinct cord blood DNA methylation patterns associated with GDM in overweight/obese pregnancies. The findings suggest potential links between epigenetic alterations and immune cell composition in GDM-exposed offspring. The identified CpG signature warrants further validation in larger, prospective cohorts to determine its clinical applicability.

Humans↗

PubMed2026

Hepatic metabolic adaptation to endurance exercise: temporal and sex differences by multiomics integration and validation.

BACKGROUND: Although endurance exercise benefits liver health, sex-specific adaptive trajectories remain unclear. This study mapped dynamic liver adaptation in males and females during prolonged training and identified underlying molecular programs. METHODS: Using publicly available time-resolved liver multi-omics data generated by the Molecular Transducers of Physical Activity Consortium (MoTrPAC), we established a computational pipeline for differential analysis of transcriptomic, proteomic, phosphoproteomic, and metabolomic data with FDR correction, followed by FGSEA pathway enrichment. Kinase activities were inferred through ortholog mapping and PhosphoSitePlus. Cross-omics co-expression networks were constructed using WGCNA and topological overlap to link omics features with physiological phenotypes. For experimental validation, liver tissues were collected from endurance-trained Sprague-Dawley rats, and key nodes were confirmed by Western blotting, qRT-PCR, and immunofluorescence/immunohistochemical staining. Public scRNA-seq data were further integrated to map multi-omics signals to single-cell resolution and assess functional changes in specific cell types. RESULTS: The hepatic response to exercise stress was stage-specific, shifting from early transcriptional activation to later proteomic and metabolic remodeling. Multi-omics integration revealed distinct sex-associated adaptive trajectories: males were more strongly associated with energy metabolism, redox-related programs, and amino acid/organic acid catabolism, whereas females showed prominent membrane lipid remodeling, proteostasis -related programs, and mitochondrial/ribosomal translational features. Single-cell analysis showed that tissue remodeling occurred without major lineage turnover, instead involving altered communication among pre-existing cell communities. Validation of PPP1R3G identified a protein-dominant exercise-responsive marker, supporting the contribution of post-transcriptional or protein-level regulation. CONCLUSIONS: Hepatic adaptation to endurance stress follows a cross-omics evolutionary pattern with sex-specific reprogramming of energy supply and homeostatic maintenance. This time-resolved framework clarifies how exercise improves liver function and supports sex-oriented metabolic interventions and therapeutic target discovery.

Animals↗

PubMed2026

Increased IL4I1 expression predicts poor survival and modulates the immune microenvironment in acute myeloid leukemia.

BACKGROUND: The immunometabolic enzyme Interleukin-4-induced-1 (IL4I1) is implicated in cancer pathogenesis, yet its specific function and clinical relevance in acute myeloid leukemia (AML) remain unclear. METHODS: Comparative analysis of IL4I1 mRNA levels between AML patients and normal controls was performed using the Cancer Genome Atlas (TCGA) and Genotype-Tissue Expression (GTEx) databases. The Kaplan–Meier survival analysis was conducted to evaluate the prognostic value of IL4I1. Functional insights were derived from analyses of differentially expressed genes (DEGs), Gene Set Enrichment Analysis (GSEA), and Gene Ontology (GO)/Kyoto Encyclopedia of Genes and Genomes (KEGG) enrichment. Immune infiltration was evaluated using the ssGSEA, ESTIMATE, quanTIseq and single-cell RNA sequencing (scRNA-seq) analysis. Finally, in vitro and in vivo functional experiments were perfromed to explore the impact of IL4I1 on AML progression and immunoregulation. RESULTS: IL4I1 expression was significantly elevated in AML compared to normal controls (p = 0.0004) and associated with poorer overall survival (p = 0.003). Bioinformatic analysis revealed that IL4I1 was linked to immune-related pathways—including humoral immune response, leukocyte interactions, and chemokine signaling—and to cellular amino acid metabolism. Its expression correlated with immune cell infiltration and checkpoint molecule expression. Experimentally, IL4I1 promoted leukemia cell proliferation in vitro and in vivo (p < 0.05). Furthermore, silencing IL4I1 suppressed M2 macrophage polarization and reduced secretion of inflammatory factors (p < 0.05). CONCLUSIONS: IL4I1 may serve as a potential biomarker for poor prognosis and an attractive target for immune-based therapeutic interventions in AML.

Humans↗

PubMed2026

Cell type resolved MR based on brain single cell eQTLs corroborated by single cell RNA sequencing uncovers neuroimmune and vascular programs in intracerebral hemorrhage.

BACKGROUND: Intracerebral hemorrhage (ICH) lacks effective neuroprotective therapies. We integrated cell type–resolved genetic inference with single-cell profiling to map putative causal programs and multicellular circuitry relevant to ICH. METHODS: Cis-eQTLs from eight human brain cell types were used as instruments for two-sample Mendelian randomization (MR), with an ICH meta-analysis from large biobanks and a stroke consortium as the outcome. Instruments were LD-pruned and restricted to strong variants (F > 10). Inverse-variance weighting (IVW) was the primary estimator, supported by robustness methods, heterogeneity/pleiotropy diagnostics, and false discovery rate control. Experimental validation used mouse collagenase ICH single-cell RNA-seq at 24 h (n = 3 sham; n = 3 ICH) with Seurat integration, composition testing, Slingshot pseudotime, and CellChat. An independent mouse cohort underwent qRT–PCR for selected genes. RESULTS: The ICH meta-analysis showed acceptable genomic control, supporting downstream MR. We identified 524 nominal gene–cell type associations, with a glia-weighted signal landscape. Enrichment implicated autophagy/mitophagy, antigen processing, cytoskeletal and vesicular trafficking, endothelial matrix–adhesion programs, ferroptosis, and myelin stress pathways. In mouse scRNA-seq, disease-associated microglia expanded with reciprocal loss of homeostatic microglia and increased neutrophils and T cells. Prioritized genes showed directional concordance; qRT–PCR confirmed ARPC3 and EIF2AK2 upregulation and TBCK and SPECC1 downregulation in ICH versus sham. Pseudotime supported a shift toward disease-associated microglial states, and CellChat indicated increased network interaction strength with microglia and endothelium as hubs. CONCLUSIONS: Cell type–specific MR combined with single-cell validation highlights neuroimmune and neurovascular programs in ICH and links genetic signals to state transitions and inferred intercellular communication.

Animals↗

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

Non-structural maintenance of chromosome condensin I complex subunit H knockdown suppresses malignant progression of esophageal squamous cell carcinoma via the Wnt/β-catenin signaling pathway.

BACKGROUND: Esophageal squamous cell carcinoma (ESCC) remains a major cause of cancer-related mortality, and effective therapeutic targets are still limited. Non-structural maintenance of chromosome condensin I complex subunit H (NCAPH) has been implicated in tumorigenesis; however, its clinical relevance, functional roles, and underlying mechanisms in ESCC are not fully defined. We aimed to characterize the expression pattern, prognostic value, biological functions, and mechanistic basis of NCAPH in ESCC. METHODS: Public datasets from The Cancer Genome Atlas (TCGA) and Gene Expression Omnibus (GEO) were analyzed to evaluate NCAPH expression and clinical associations. Single-cell RNA sequencing (scRNA-seq) data were used to map cell-type-specific distribution of NCAPH in tumor and adjacent tissues. NCAPH was silenced in KYSE150 and KYSE510 cells using lentiviral short hairpin RNAs (shRNAs), followed by Cell Counting Kit-8 (CCK-8), colony formation, wound-healing, and Transwell migration/invasion assays. A nude mouse xenograft model was established to assess the effect of NCAPH knockdown in vivo. RNA sequencing (RNA-seq), quantitative polymerase chain reaction (qPCR), western blotting, and enzyme-linked immunosorbent assay (ELISA) were performed to explore potential mechanisms. RESULTS: NCAPH was consistently upregulated in ESCC across multiple cohorts and was associated with unfavorable clinicopathological features and poorer survival. Functional assays demonstrated that NCAPH knockdown significantly inhibited ESCC cell proliferation, migration, invasion, and clonogenic growth. In vivo, NCAPH silencing suppressed xenograft tumor growth. Mechanistically, transcriptomic profiling and molecular validation indicated attenuation of Wnt/β-catenin signaling following NCAPH depletion, accompanied by reduced β-catenin and downstream targets. CONCLUSIONS: NCAPH promotes malignant progression of ESCC, at least in part through activation of the Wnt/β-catenin pathway, and may serve as a potential biomarker and therapeutic target.

Esophageal squamous cell carcinoma (ESCC)↗

PubMed2026

Construction of molecular signatures based on the co-expression network of NECSO-related gene TRPM4 and its prognostic value in hepatocellular carcinoma.

BACKGROUND: Hepatocellular carcinoma (HCC) demonstrates significant prognostic variability that is not entirely accounted for by traditional staging systems. Necrosis by sodium overload (NECSO) is an emerging programmed cell death pathway, but its clinical relevance in HCC remains undefined. Therefore, this study aimed to identify TRPM4-associated core genes, develop and validate a prognostic signature, and investigate its relationship with the tumor immune microenvironment, tumor mutational burden, and single-cell expression patterns in HCC. METHODS: We integrated transcriptomic, clinical, and mutational datasets from The Cancer Genome Atlas-Liver Hepatocellular Carcinoma (TCGA-LIHC) (n=421) and Gene Expression Omnibus (GEO) cohorts (n=115) to identify genes co-expressed with TRPM4-a key NECSO mediator-and those differentially expressed in HCC. A prognostic signature was developed using least absolute shrinkage and selection operator (LASSO)-Cox regression and validated through survival analysis, time-dependent receiver operating characteristic (ROC) curves, and multivariate Cox regression analysis. The immune landscape was characterized using CIBERSORT, somatic mutation data were used to calculate tumor mutational burden (TMB) and assess its correlation with the risk score, and single-cell RNA sequencing (scRNA-seq) resolved cell-type-specific expression patterns. RESULTS: From 294 TRPM4-associated core genes, we identified an 11-gene signature (BRSK1, MMP1, GRIN2D, GP6, MYOM2, N4BP3, CCDC112, TSEN54, MAP3K9, SPP1, B3GNT4) that independently predicted overall survival (OS) (hazard ratio =5.419, P<0.001) with areas under the curve (AUCs) of 0.779, 0.693, and 0.701 at 1, 3, and 5 years. These values were superior or comparable to conventional clinicopathologic variables after direct comparison. High-risk patients exhibited an immunosuppressive microenvironment, characterized by enrichment of M0 macrophage, a higher M2/M1 ratio (P<0.001) and distinct immune checkpoint profiles. When integrated with TMB, the prognostic stratification was further refined: high-TMB/high-risk patients had poorest outcomes (median OS, 15.3 months), while low-TMB/low-risk patients had the most favorable survival (median OS, 68.7 months). Single-cell analysis revealed that MMP1 was induced in cancer-associated fibroblasts (CAFs) and SPP1 was downregulated in macrophages, single-cell risk scores confirmed TAFs and macrophages as the main contributors to the prognostic model. CONCLUSIONS: The TRPM4-centered 11-gene signature provides robust and independent prognostic stratification in HCC by integrating immune, mutational, and single-cell features. This signature serves as a potential tool for prognostic evaluation and may help inform immunotherapeutic strategies for HCC.

Hepatocellular carcinoma (HCC)↗

PubMed2026

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↗

PubMed2026

Single-nucleus profiling of postmortem diffuse midline gliomas identifies mitochondrial biogenesis as a resistance mechanism to imipridone therapy.

BACKGROUND: Imipridone ONC201 is the first FDA-approved therapy for H3K27-altered diffuse midline glioma; however, clinical responses remain limited. Defining tumor-intrinsic determinants and microenvironmental, extrinsic factors that shape sensitivity or resistance to imipridones will identify actionable therapeutic opportunities and inform improved clinical strategies. METHODS: To identify mechanisms of imipridone resistance, we obtained postmortem brain tissue from DMG patients who had received imipridones and/or standard care. Single-nucleus RNA and open-chromatin sequencing were performed on N = 22 cases. Immunofluorescence-based myeloid phenotyping was performed on N = 46 cases. Mitochondrial copy-number analysis was performed on N = 19 cases. Validation of imipridone sensitivity, its effect on mitochondrial density, and its synergy with inhibition of mitochondrial biogenesis were assessed in DMG primary cells. RESULTS: We established a single-cell RNA/open-chromatin atlas from postmortem DMG cases and found imipridone treatment resulting in regressed mesenchymal transition, reduced myeloid-derived suppressive cells, and reversed aberrant H3K27-altered enhancer activity. Resistant tumors showed increased mitochondrial density, turnover, and membrane potential. Mitochondrial biogenesis and PPARGC1A emerged as resistance biomarkers and actionable targets. CONCLUSIONS: These studies implicate mitochondrial biogenesis as a biomarker of imipridone resistance and a focus for the development of combinatorial strategies to provide effective therapeutic options for a challenging pediatric brain tumor.

Humans↗

PubMed2025

Cellular imbalance in proximal and distal lung of CFTR-/- sheep in utero and at birth.

BACKGROUND: The Lung is the major focus of therapeutic approaches for the inherited disorder cystic fibrosis (CF) as without treatment lung disease is life-limiting. However, the initiating events that predispose the CF lung to cycles of infection, inflammation and resultant tissue damage are still unclear. Inflammation may occur in the CF lung prior to birth in human and several large animal models suggesting an in utero origin for the disease and encouraging further studies prior to birth. METHODS: Here we used the sheep model of CF (CFTR-/-) and age-matched wild-type (WT) sheep of the same breed to investigate the single cell transcriptomes of proximal and distal lung tissue at 80 days and 120 days of gestation and at term (147 days). Single cell RNA-seq was performed on tissues from 4 to 7 animals of each genotype (WT and CFTR-/-) at each time point. RESULTS: At term, FOXJ1-expressing ciliated cells are overrepresented in both lung regions from CFTR-/- lambs, while secretory epithelial and basal cells are underrepresented in proximal lung, as are T cells and monocytes in distal lung. The imbalance in ciliated and basal cells was confirmed by immunohistochemistry. At 120 days of gestation, lymphoid cells are slightly more abundant in proximal and distal lung from CFTR-/- animals compared to WT, consistent with the transient CF-associated inflammatory response in utero. At 80 days of gestation, T and B cells are underrepresented in both lung regions. CONCLUSIONS: The differences in epithelial cell abundance observed in the CFTR-/- lambs at term may reflect sequelae from the loss of CFTR on lung development and differentiation in utero. These findings provide novel insights into the cellular mechanisms of pathology and may be relevant to the design of new therapeutic approaches for CF lung disease.

Animals↗

PubMed2024

Single cell RNA sequencing provides novel cellular transcriptional profiles and underlying pathogenesis of presbycusis.

Age-related hearing loss (ARHL) or presbycusis is associated with irreversible progressive damage in the inner ear, where the sound is transduced into electrical signal; but the detailed mechanism remains unclear. Here, we sought to determine the potential molecular mechanism involved in the pathogeneses of ARHL with bioinformatics methods. A single-cell transcriptome sequencing study was performed on the cochlear samples from young and aged mice. Detection of identified cell type marker allowed us to screen 18 transcriptional clusters, including myeloid cells, epithelial cells, B cells, endothelial cells, fibroblasts, T cells, inner pillar cells, neurons, inner phalangeal cells, and red blood cells. Cell-cell communications were analyzed between young and aged cochlear tissue samples by using the latest integration algorithms Cellchat. A total of 56 differentially expressed genes were screened between the two groups. Functional enrichment analysis showed these genes were mainly involved in immune, oxidative stress, apoptosis, and metabolic processes. The expression levels of crucial genes in cochlear tissues were further verified by immunohistochemistry. Overall, this study provides new theoretical support for the development of clinical therapeutic drugs.

Animals↗