SEARCH · Search PubMed

Results for “Single-cell RNA sequencing (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 37 records · Page 2Linked to original sources

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

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

SCMO: a deep learning model integrating the single-cell resolution TME ecosystem and multi-omics for survival prediction in CRC patients.

BACKGROUND: Colorectal cancer (CRC) remains a leading cause of global cancer mortality, highlighting the need for precise survival prediction to guide clinical decisions. Although tissue-level multi-omics is widely utilized for survival prediction, its limited resolution cannot capture tumor heterogeneity. Single-cell RNA sequencing (scRNA-seq) enables dissection of the tumor microenvironment (TME) at cellular resolution, supporting personalized prognostic assessment. METHODS: We collected 213 CRC scRNA-seq samples and established a CRC-specific TME atlas comprising 339,060 cells. Using this atlas as a reference, we deconvolved bulk RNA-seq data from TCGA-CRC cohort with the EcoTyper algorithm to reconstruct TME features. Clinical, genomic, and transcriptomic data were obtained from the Xena platform; microbial data were sourced from the BIC database. We integrated TME and multi-omics features through a self-normalizing neural network to construct a deep learning model (single-cell resolution TME ecosystem with multi-omics data [SCMO]) for survival prediction. To enhance interpretability, we utilized the Integrated Gradients algorithm and spatial transcriptomic data to analyze multi-omics and TME features. We performed anticancer drug screening with tumor necrosis factor receptor-associated protein 1 (TRAP1), a critical feature according to the Integrated Gradients algorithm, as a potential target. RESULTS: We identified 13 survival-related TME features from the CRC-specific atlas: 12 cell states and one multi-cellular ecosystem. SCMO, which combined TME and multi-omics features, improved survival prediction and outperformed existing methods, achieving a concordance index of 0.762. The SCMO demonstrated robust performance for long-term predictions, achieving areas under the curve (AUCs) of 0.752, 0.772, and 0.869 for 1-, 3-, and 5-year predictions in the training set, with corresponding test set AUCs of 0.639, 0.756, and 0.772. TME features from the SCMO model revealed that ecosystem density increased with CRC malignancy. Multi-omics features included TRAP1 as a potential drug target. Drug screening identified saikosaponin A as a novel TRAP1 inhibitor, and its anticancer activity was validated in vitro. We developed SCMO-Lite, a simplified model incorporating 12 high-attribution-weight multi-omics features, which demonstrated robust risk stratification. CONCLUSIONS: SCMO combines analytical precision with biological interpretability, offering novel insights for oncology survival prediction.

Humans↗

PubMed2025

Multi-omics analysis reveals distinct spatial compartmentalization of lung repair niches in pediatric ARDS.

BACKGROUND: Pediatric acute respiratory distress syndrome (PARDS), often triggered by viral infections, is a life-threatening condition. Despite its severity, children demonstrate significantly better survival rates and superior lung repair compared to adults. However, the mechanisms underlying this age-specific advantage remain incompletely understood. PATIENTS AND METHODS: We conducted a pilot multi-omics study of influenza-associated PARDS integrating single-cell RNA sequencing (scRNA-seq) of pediatric lung tissue and bronchoalveolar lavage fluid (BALF), spatial transcriptomics, and plasma proteomics. Analyses were harmonized with the Human Lung Cell Atlas (HLCA) reference, reanalysis of public pediatric PARDS airway scRNA-seq, and contextual comparisons to adult lethal COVID-19 lung. RESULTS: Tissue scRNA-seq and spatial data indicated outcome-linked divergence in PARDS. Survivor showed spatially restricted repair with preserved alveolar type II (AT2) cells, AT2-to-alveolar type I (AT1) differentiation signatures, and higher KRT17, whereas fatal case and adults exhibited diffuse immune activation with pro-fibrotic and pro-apoptotic signaling. In BALF, KRT17-positive airway stress–repair epithelial cells (hillock-like) increased from the acute to recovery phase, and plasma proteomics showed higher circulating KRT17 in survivors. HLCA-based label transfer strengthened cell-type definitions and enabled pediatric–adult comparisons suggesting biological and developmental differences; the adult lethal COVID-19 atlas provided a benchmark with attenuated epithelial repair and prominent collagen CTHRC1-pathologic fibroblasts. Fibroblast programs were regionally compartmentalized, with injury-enriched CTHRC1+ states versus alveolar fibroblasts in preserved areas, and showed stronger injury–homeostasis anti-correlation in fatalities. Myeloid remodeling included BALF transitions from FCN1-high inflammatory states toward FABP4-positive resident-like states, consistent with public pediatric datasets showing reduced inflammatory and interferon-stimulated gene (ISG) modules and severity-linked increases in aged neutrophils. CONCLUSIONS: This pilot multi-omics case series outlines putative pediatric lung repair niches in influenza-associated PARDS. KRT17-positive transitional epithelium, preserved AT2 differentiation, and restoration of resident-like macrophages may align with recovery, whereas diffuse immune activation and CTHRC1-enriched fibroblast programs may accompany worse outcomes. HLCA-guided annotations and adult benchmarks indicate possible age-related differences, warranting validation in larger multi-center cohorts.

Humans↗

PubMed2026

Senescent fibroblasts drive CD8+ T cell dysfunction in colorectal cancer via CD36-mediated lipid transfer and peroxidation.

BACKGROUND: Functional exhaustion of tumor-infiltrating CD8+ T cells represents a hallmark of colorectal cancer (CRC) immunosuppression, though its mechanistic drivers remain elusive. Given the established correlation between CRC progression and stromal senescence characterized by pathological lipid accumulation and impaired immunity, we investigated whether and how senescent fibroblasts actively regulate CD8+ T cell dysfunction. METHODS: Single-cell RNA sequencing (scRNA-seq) analysis was conducted to unveil the diverse fibroblast populations and the significant lipid metabolism changes between senescent fibroblasts and non-senescent fibroblasts in human CRC specimens and adjacent normal mucosa. Machine-learning identified senescent fibroblasts with a distinct gene signature. Cell-cell communication analysis was used to evaluate the interactions between senescent fibroblasts and CD8+ T cells in colorectal cancer. Co-culture experiments were conducted among senescent fibroblasts, CD8+ T cells and patient-derived organoids of CRC (CRC-PDOs), with the results evaluated with high-content imaging and propidium iodide/Hoechst 33,342 staining. Flow cytometry, ELISA and lipid pulse-chase with BODIPY FL C16 were performed to detect the alterations of CD8+ T cell cytotoxic function and metabolic status. AOM/DSS-induced CRC mouse model was used to conduct in vivo validation to evaluate whether senolytics could suppress CRC progression. Patients from the Cancer Genome Atlas colorectal cancer cohort were stratified into CD36-high and CD36-low groups by median expression, and drug sensitivity for GDSC2 compounds was predicted computationally using the oncoPredict R package. RESULTS: ScRNA-seq demonstrated the specific cell population presence and divergence of senescent fibroblasts between neoplastic and histologically normal adjacent cell clusters in CRC. Random Forest was employed for cell senescence classification. Feature importance analysis identified five genes as key contributors to the model’s decision process. Cell-cell communication analysis revealed enhanced interactions between senescent fibroblasts and CD8+ T cells in CRC. Co-culture of senescent fibroblasts significantly impaired the cytotoxic functions of CD8+ T cells on CRC-PDOs, which was reflected by the declined proportions of granzyme B (GZMB) + and interferon gamma (IFNγ) + CD8+ T cells and enhanced viability of CRC-PDOs. Mechanistically, the co-culture with senescent fibroblasts promoted the lipid shuttling into CD8+ T cells to induce lipid peroxidation and downstream impairment of cytotoxicity. Furthermore, the inhibition of CD36, the specific scavenger receptor for lipid uptake of CD8+ T cells, effectively suppressed lipid transfer and peroxidation thereby preserving the effector functions of CD8+ T cells and ultimately promoting tumor apoptosis. Complementarily, in vivo senolytic treatment significantly suppressed CRC progression in AOM-DSS CRC mouse models. Top 12 therapeutic agents were identified significantly enhanced predicted efficacy in CD36-high tumors. CONCLUSIONS: Our study identified a substantial population of senescent fibroblasts in human CRC through single cell transcriptomics, machine-learning and clinical biopsies. These senescent fibroblasts impair CD8+ T cell-mediated killing of CRC-PDOs via CD36-dependent lipid transfer, suggesting senolytic targeting of stromal cells as a promising immunotherapeutic strategy for CRC.

Colorectal Neoplasms↗

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

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

Ligand-based directed differentiation to produce granulosa-like cells expressing steroidogenic enzyme genes.

The ovarian granulosa cells are responsible for producing hormones and supporting oocytes through maturation and meiotic resumption. There is a need to generate granulosa-like cells (GLCs) from human induced pluripotent stem cells (hiPSCs) to better model human gonadal development and to test the effects of exogenous or pharmaceutical compounds on the ovary. Here we report a rapid ligand-based protocol for differentiating hiPSCs into cells that express markers of the transient developmental lineages and steroidogenic pathway genes. Single-cell RNA-sequencing (scRNA-seq) analysis identified canonical granulosa cell genes were expressed in a subset of cells and identified new genes of interest that were significantly associated with computationally modeled pseudotime. HSD17B1 was expressed in resulting GLCs but at low levels, suggesting an immature granulosa cell phenotype. The GLCs were produced using a simple culture method that could be augmented for granulosa cell functions such as sustaining oocyte growth. Producing GLCs through protocols such as this one is a first step toward designing large-scale ovarian endocrinology assays and developing personalized cell-based fertility and hormone restoration technologies in the future. This rapid protocol produced cells that express steroidogenic enzyme genes etoc blurb. Kubo and colleagues present a 5-day rapid protocol to generate immature granulosa-like cells from hiPSCs. Cells differentiated with inhibition of DKK1, a WNT signaling target gene, expressed gonadal ridge markers and FOXL2 transcripts and protein. Additionally, steroidogenic enzyme genes were expressed. A small population of differentiated cells were identified as expressing early-stage granulosa cell genes by single-cell RNA-seq.

Female↗

PubMed2026

shinyDeepGxP: a user-friendly R shiny app for predicting surface protein abundance from scRNA-seq expression using deep learning in blood cells.

MOTIVATION: Understanding accurate immune cell heterogeneity and function in single-cell datasets requires access to protein-level information, which is often unavailable due to experimental limitations. RESULTS: We present shinyDeepGxP, an interactive web application featuring our deep learning model, DeepGxP, for predicting surface protein abundance from single-cell RNA-sequencing (scRNA-seq) data. This platform makes DeepGxP accessible to researchers without programming skills. Users can upload scRNA-seq count matrices and use "Predict Protein" to predict the abundance of 224 biologically relevant surface proteins. shinyDeepGxP provides visualizations to help identify distinct cell populations based on predicted protein profiles. Moreover, users can choose "Explore Model" to reveal key RNA predictors and their associated biological pathways for each protein. Overall, shinyDeepGxP is a user-friendly, freely available web tool that provides protein-level detail for RNA-only single-cell datasets, enabling multimodal discovery without additional experiments. AVAILABILITY AND IMPLEMENTATION: shinyDeepGxP can be launched on https://shiny.crc.pitt.edu/deepgxp/.

Journal Article↗

PubMed2026

Deciphering novel targets in salivary gland pleomorphic adenoma by integrating plasma proteomics and parotid transcriptomics analyses.

BACKGROUND/PURPOSE: Pleomorphic adenoma (PA) is the most common salivary gland benign tumor, with its molecular drivers elusive due to a lack of experimental models. This study aimed to decipher novel targets in PA by systematically integrating plasma protein quantitative trait loci (pQTL)-based Mendelian randomization (MR) with multi-omics profiling of parotid gland tissues. MATERIALS AND METHODS: We performed two-sample MR using 5450 plasma pQTLs and genome-wide association study summary for benign or broader salivary gland diseases from FinnGen consortium. Bulk RNA-sequencing (RNA-seq) and single-cell RNA-seq (scRNA-seq) comparing PA to normal tissue were used for transcriptomic validation. Immunohistochemistry (IHC) was applied for protein-level validation in human PA, adenoid cystic carcinoma (ACC), and murine inflammatory lesions. RESULTS: MR identified 12 plasma proteins associated with benign salivary gland tumor risk. Transmembrane serine protease 6 (TMPRSS6) was the only protein significantly risk-increasing for both benign and broader salivary gland diseases. Strikingly, mitogen-activated protein kinase kinase 4 (MAP2K4) showed opposite MR effects between benign and all-lesion outcomes. Bulk RNA-seq showed limited concordance with MR findings, while scRNA-seq revealed a unique plastic epithelium and partially validated candidates at cellular resolution. Critically, IHC confirmed MAP2K4 protein overexpression specifically in human PA, but not in ACC or inflammatory lesions, while TMPRSS6 was downregulated in established pathologies despite its genetic risk association. CONCLUSION: By integrating plasma proteome-based causal inference with parotid tissue multi-omics, this study unveils MAP2K4 as a potential PA-specific driver. This integrative framework provides novel, context-specific targets for further functional investigation in salivary gland tumorigenesis.

Gene expression profiling↗

PubMed2026

GeneExt: a gene model extension tool for enhanced single-cell RNA-seq analysis.

MOTIVATION: Incomplete gene models negatively impact single-cell gene expression quantification. This is particularly true in non-model species where often gene 3' ends are inaccurately annotated, while most scRNA-seq methods only capture the 3' transcript region. This results in many genes being incorrectly quantified or not detected. RESULTS: GeneExt leverages scRNA-seq data to refine gene annotations. We exemplify GeneExt usage and its impact on the gene expression quantification of eight non-model organism single-cell atlases. By extending and homogenizing gene annotations, our tool will help improve biological interpretation and cross-species comparisons of cell type expression atlases. AVAILABILITY: GeneExt is available at https://github.com/sebepedroslab/GeneExt (DOI: https://doi.org/10.5281/zenodo.18712940) under a GNU General Public license, together with test data and usage instructions.

Software↗

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

Predicting gene-specific regulation with transcriptomic and epigenetic single-cell data.

MOTIVATION: Analysis of single cell ATAC-seq and RNA-seq data has allowed to gain unprecedented insights into gene regulation by allowing to define cell type-specific regulatory regions and their effects on gene expression. While powerful, such analysis is challenging due to the inherent sparsity of single cell data. RESULTS: We present a new approach, MetaFR, to learn gene-specific models that link open-chromatin variation from scATAC-seq data to gene expression from scRNA-seq. Using efficient regression trees, we illustrate that accurate expression prediction models can be learned on the single-cell or meta-cell level. Validation was done using fine-mapped eQTLs. Meta-cell models were found to outperform single-cell models for most genes. Comparison to the SOTA method SCARlink revealed advantages of MetaFR in terms of runtime and prediction performance. MetaFR thus allows time-efficient analysis and obtains reliable models of gene expression prediction, which can be used to study gene regulation in any organism for which scRNA-seq and scATAC-seq data is available. AVAILABILITY AND IMPLEMENTATION: MetaFR is available under https://github.com/SchulzLab/MetaFR.

Single-Cell Analysis↗

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↗

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↗