SEARCH · PubMed Health

Results for “single-cell RNA-seq”

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

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

At least 37 records · Page 2Linked to original sources

PubMed2026

Multimodal Analysis Reveals Aberrant Expression of SUMO2 and Its Significant Association With Key Mechanisms of Metabolic Pathways in Hepatocellular Carcinoma.

BACKGROUND: Hepatocellular carcinoma (HCC) is the third leading cause of cancer-related deaths worldwide. However, the role of small ubiquitin-like modifier 2 (SUMO2), a core member of the small ubiquitin-like modifier (SUMO) family, regarding its expression patterns and metabolism-related functions in HCC remains inadequately understood. METHODS: A multidimensional analytical framework was applied, integrating immunohistochemistry (153 HCC vs. 21 non-HCC samples), proteomics (159 paired samples), bulk transcriptomics (3240 HCC vs. 2267 non-HCC samples), single-cell RNA sequencing (RNA-seq) (10 HCC vs. 8 non-HCC samples), spatial transcriptomics, and external CRISPR/Cas9 functional genomics data. Systematic analyses included standardized mean difference (SMD), pathway enrichment, pseudotime trajectory inference, in silico knockout, cell-cell communication, metabolic flux scoring, immune infiltration, clinical correlation, drug sensitivity prediction, and molecular docking. RESULTS: At the protein level, immunohistochemistry (nuclear positivity) and external proteomic data collectively demonstrated consistent SUMO2 overexpression in HCC. Consistent upregulation was also observed at the mRNA level across large-scale cohorts. Single-cell RNA-seq and spatial transcriptomics localized SUMO2 enrichment to malignant hepatocytes and tumor-dominant regions. CRISPR-mediated SUMO2 knockout suppressed proliferation in multiple HCC cell lines. Mechanistically, high SUMO2 expression was significantly associated with metabolic reprogramming involving glycolysis/gluconeogenesis, pyruvate metabolism, and the tricarboxylic acid cycle. SUMO2-high malignant hepatocyte subpopulations exhibited enhanced activity of the macrophage migration inhibitory factor signaling axis and enhanced iron-sensor interactions. Further, the immune infiltration analysis revealed a negative correlation between SUMO2 expression and M1 macrophages and a positive correlation between follicular helper T cells and regulatory T cells. Clinically, elevated SUMO2 levels were found to be associated with adverse prognostic features. Furthermore, high SUMO2 expression was associated with increased sensitivity to dasatinib, and molecular docking simulations predicted potential binding between SUMO2 and dasatinib, with a Vina score of -8.5 kcal/mol. CONCLUSIONS: SUMO2 is aberrantly expressed at the protein, mRNA, single-cell, and spatial transcriptomic levels in HCC and is significantly associated with metabolic reprogramming and altered migration inhibitory factor (MIF)-mediated intercellular communication, suggesting its potential as a novel biomarker for diagnosis and treatment.

Humans↗

PubMed2026

SSB deficiency-induced R-loop accumulation triggers podocyte inflammation in DKD.

INTRODUCTION: Diabetic kidney disease (DKD) is fundamentally a podocytopathy in which sterile inflammation plays a central pathogenic role, yet the upstream triggers that initiate inflammatory cascades in podocytes remain elusive. R-loops are critical regulators of genomic stability, and their pathological accumulation triggers DNA damage and innate immune activation. Whether R-loop dysregulation contributes to podocyte-driven inflammation in DKD is unknown. METHODS: We integrated single-cell transcriptomic profiling, dual machine learning algorithms, and functional experiments to dissect the R-loop regulatory network in the diabetic kidney. RESULTS: Integrated analysis of human diabetic kidney single-cell RNA-seq data revealed a globally compromised R-loop regulatory network selectively within podocytes. Intersection of podocyte-specific transcriptomic shifts with validated R-loop regulators identified 93 candidate genes, from which dual machine learning algorithms pinpointed SSB (Sjögren syndrome antigen B) as the principal podocyte-selective R-loop resolver and a superior diagnostic biomarker (AUC = 0.983). SSB expression was selectively downregulated in diabetic podocytes and showed the strongest positive correlation with the R-loop resolution module. Mechanistically, SSB loss impaired RNA splicing and stability pathways, leading to aberrant R-loop accumulation that activated the cGAS-dependent inflammatory signaling in podocytes. In two murine DKD models and high glucose-challenged podocytes, SSB was markedly reduced. Remarkably, SSB knockdown in podocytes alone sufficed to trigger R-loop accumulation and pro-inflammatory cytokine expression, whereas both RNase H1-mediated R-loop removal and cGAS co-depletion blunted this response. DISCUSSION: These findings suggest that an SSB-governed R-loop -cGAS -inflammatory signaling axis may link genomic instability to podocyte inflammation and contribute to DKD progression, nominating R-loop homeostasis as a previously unrecognized potential therapeutic target.

Podocytes↗

PubMed2026

SHMT2: a Metabolic and Immune Biomarker of Aggressive Lung Adenocarcinoma.

Serine/glycine-one-carbon (SGOC) metabolism is frequently altered in lung adenocarcinoma (LUAD), but its relationship to tumor behavior and predicted immunotherapy responsiveness remains incompletely defined. Metabolomic profiling of 23 paired LUAD and adjacent normal lung tissues was performed using internal extractive electrospray ionization mass spectrometry. Transcriptomic and clinical data from The Cancer Genome Atlas LUAD cohort (TCGA-LUAD) were analyzed to assess SHMT2 expression, prognosis, differentially expressed genes, and immune-related features. Predicted response to immune checkpoint blockade was evaluated using Tumor Immune Dysfunction and Exclusion (TIDE) and The Cancer Immunome Atlas (TCIA), and drug sensitivity was inferred using oncoPredict. Single-cell RNA-seq data were used to examine the cellular distribution of SHMT2. Experimental validation included quantitative reverse-transcription PCR (RT-qPCR), western blotting, Human Protein Atlas (HPA) immunohistochemistry, and short hairpin RNA (shRNA)-mediated SHMT2 knockdown followed by proliferation, wound-healing and colony formation assays. Metabolomic analysis identified glycine, serine, and threonine metabolism as a prominently altered pathway in LUAD. SHMT2 was upregulated in LUAD and associated with worse overall survival and adverse clinicopathological features. SHMT2-high tumors displayed enrichment of cell-cycle and SGOC-related transcriptional programs, lower immune and stromal scores, and reduced predicted responsiveness to immunotherapy. Single-cell analysis showed relative enrichment of SHMT2 expression in B cell populations. In vitro, SHMT2 was overexpressed in LUAD cells, and its knockdown suppressed proliferation, migration, and clonogenic growth. Collectively, SHMT2 is associated with SGOC metabolic reprogramming, aggressive tumor phenotypes, and an immune-disadvantaged state in LUAD, supporting its potential relevance as a biomarker; therapeutic targeting requires additional pharmacologic and in vivo validation.

Humans↗

PubMed2026

A single-cell meta-analysis evidences transposable element dysregulation in sex-based differences in Parkinson's disease.

Transposable elements (TEs) (mobile genetic elements comprising ∼45% of the human genome) have recently emerged as potential contributors to Parkinson's disease (PD); however their role and sex-specific impact remain poorly understood. Here, we present the first integrative meta-analysis of TE expression across 4 substantia nigra single-nucleus RNA-seq datasets, comprising a total of 66 donors, generating a cell-type-resolved atlas of TE dysregulation in PD. We identified widespread TE activation across major brain cell types (i.e. neurons, astrocytes, oligodendrocytes and microglia), with marked upregulation of L1s in neurons and HERVs in oligodendrocytes. Sex-stratified analyses revealed distinct male- and female-biased TE signatures, indicating regulatory programs uniquely affected in each sex, including MIR elements in microglia and Alu subfamilies in neurons. Correlation and genomic proximity analyses also uncovered TE-gene associations linked to important PD pathways such as neuroinflammation or myelination. Collectively, our study positions TEs as potential sex-modulated contributors to PD pathology and also provides a public web resource (PATOSS) to explore PD-associated TE transcriptional deregulation.

Parkinson's disease↗

PubMed2026

The role of stem cells in pituitary tumour formation.

Pituitary tumours are intracranial neoplasms that pose significant clinical challenges due to their potential for recurrence, therapeutic resistance and resultant endocrine dysfunction and mass effects. In the normal anterior pituitary, resident pituitary stem cells (PSCs) contribute to tissue homeostasis and cellular turnover. The extent to which PSCs contribute to tumourigenesis is not known, but an increasing number of studies have been aiming to address this. In this review, we summarise current evidence implicating PSCs and tumour stem-like populations in pituitary tumour biology, including potential roles in tumour initiation, maintenance and progression. We outline practical criteria for defining tumour stem cells and evaluate findings from functional studies of human tumours, emerging single-cell and spatial transcriptomic datasets and murine lineage-tracing models. We also provide a curated overview of published single-cell RNA sequencing studies of pituitary tumours, highlighting reported stem/progenitor populations and transcriptional signatures across tumour subtypes and propose a framework for future genomic analyses. Finally, we discuss the translational implications of these findings, including the potential for targeting stem-like populations and their associated signalling pathways.

Humans↗

PubMed2026

Benchmarking computational decontamination of ambient RNA.

Gene expression profiling of single cells using single-cell and single-nucleus RNA sequencing (sxRNA-seq) enables researchers to characterize cellular heterogeneity and unraveling complex biological processes at unprecedented resolution. However, sxRNA-seq faces challenges due to the presence of ambient RNA, extraneous RNA molecules not originating from the cells of interest. Sample preparation is a major source of ambient RNA, where harsh conditions can lead to cell lysis and the release of intracellular RNA. This inescapable inclusion of ambient RNA can cause erroneous results and hinder downstream analyses. To address this issue, various methodologies have been developed to identify, quantify, and remove ambient RNA. Here, we rigorously evaluate 7 state-of-the-art methodologies for ambient RNA removal using simulated datasets, species-mixing experiments of varying complexities, and genotype-mixing experiments. We find that no single method performs the best across all datasets and metrics, but CellBender, DecontX and SoupX generally perform well.

ambient RNA↗

PubMed2026

Synovial short-lived plasma cells mediate adalimumab resistance in rheumatoid arthritis via MIF-CD74 axis-driven, partially TNF-α-independent inflammation.

OBJECTIVE: Synovial plasma cell infiltration predicts inadequate response to adalimumab in patients with rheumatoid arthritis (RA), yet the cellular and molecular mechanisms underlying this association remain unclear. This study aimed to dissect the functional heterogeneity of synovial plasma cells between adalimumab responders and non-responders at single-cell resolution, and to identify the molecular pathways driving treatment resistance. METHODS: This study was based on a prospective clinical cohort of 101 RA patients receiving adalimumab, from which synovial tissues of 8 patients (4 ACR20 responders and 4 non-responders) were profiled by 10x Genomics single-cell RNA sequencing (66,539 high-quality cells). A systematic ligand-receptor screening was performed to identify candidate signaling axes. Core findings were validated at four levels: an independent single-cell validation cohort (n = 4), external bulk RNA-seq cohorts (GSE15602, GSE47726), multiplex immunofluorescence on synovial tissues (n = 9 per group), and in vitro functional experiments using patient-derived peripheral blood monocyte-derived macrophages stimulated with recombinant human MIF under pharmacological intervention with adalimumab, the MIF inhibitor ISO-1, and an anti-CD74 neutralizing antibody. RESULTS: Plasma cells were significantly enriched in non-responder synovium, with a heterogeneous pattern characterized by quantitative accumulation of long-lived plasma cells (LLPCs) and functional dominance of short-lived plasma cells (SLPCs): SLPCs contributed 58.15% of total ribosomal module activity and preferentially overexpressed MIF. Systematic screening of 145 candidate ligand-receptor pairs identified MIF-CD74 as the only axis satisfying all four independent evidence layers. Tissue-level immunofluorescence confirmed that approximately 95% of synovial CD138+ plasma cells in non-responders co-expressed MIF, compared with approximately 45% in responders. In vitro, rh-MIF upregulated macrophage activation markers (CD74, CD80, CD86, HLA-DR) and induced IL-6 and TNF-α secretion. Adalimumab neutralized supernatant TNF-α but failed to suppress MIF-driven IL-6 and IL-1β activation, whereas ISO-1 and anti-CD74 effectively blocked MIF-induced effects at all levels examined. These findings were replicated in patient-derived PBMC macrophages. CONCLUSION: In adalimumab-resistant RA, a functionally active SLPC subset drives partially TNF-α-independent macrophage inflammation through the MIF-CD74 axis, representing a resistance pathway not fully addressed by anti-TNF therapy. Targeting MIF or CD74 blocked this axis in vitro, supporting MIF-CD74-directed precision intervention.

Adalimumab↗

PubMed2026

ScRNA-seq analysis reveals the effects of nitrite stress on the endocrine system of the eyestalk in Litopenaeus vannamei.

Nitrite is a harmful substance generated in Litopenaeus vannamei farming systems, largely originating from the inadequate breakdown of surplus feed and shrimp feces. Its accumulation in the water can affect the growth and physiological functions of shrimp, damage the immune system, and even cause mass mortality, thus becoming a key environmental factor restricting the green development of the industry. Under nitrite stress, the eyestalk, as an important neuroendocrine regulatory center in crustaceans, participates in the stress adaptation of the organism and exerts a protective effect by regulating energy metabolism and immune function. However, the molecular regulatory mechanism of the eyestalk in response to nitrite stress remains unclear. In this study, single-cell RNA sequencing (scRNA-seq) technology was used to analyze the heterogeneity of eyestalk cells in L. vannamei under nitrite stress. A total of 18, 394 high-quality cells were obtained, and six major cell subpopulations, including Neurosecretory cell, Motor neuron, Sensory neuron, Interneuron, Neurogliocyte, and Support cell, were identified. Differential expression analysis identified 839 differentially expressed genes, and different cell types showed distinct specific responses to nitrite stress. Functional enrichment analysis indicated that pathways such as glycolysis, oxidative phosphorylation, ribosome function, and endoplasmic reticulum protein processing were significantly activated, while signal transduction and DNA repair-related pathways were inhibited. Further analysis revealed that nitrite stress could induce mitochondrial function changes and trigger oxidative stress, thereby affecting the neuroendocrine system function of the eyestalk. This study provided insights into transcriptomic responses of the eyestalk to nitrite stress at the single-cell level, laying a theoretical foundation for the management of aquaculture environments.

Animals↗

PubMed2026

ARCADIA reveals spatially dependent transcriptional programs through integration of scRNA-seq and spatial proteomics.

MOTIVATION: Cellular states are strongly influenced by spatial context, but single-cell RNA sequencing (scRNA-seq) loses information about local tissue organization, while spatial proteomic assays capture limited marker panels that constrain transcriptomic inference. Integrating these modalities can elucidate how spatial niches shape transcriptional programs, yet existing approaches depend on either feature-level correspondence such as gene-protein linkage or cell-level barcode pairing, which is often unavailable. RESULTS: We present ARCADIA (ARchetype-based Clustering and Alignment with Dual Integrative Autoencoders), a generative framework for cross-modal integration that operates without cell barcode pairing and does not assume direct feature-to-feature correspondence. ARCADIA identifies modality-specific archetypes, that is, convex combinations of cells representing extreme phenotypic states, and aligns these anchors across modalities by minimizing the discrepancy between their cell-type composition profiles. The aligned archetypes define a shared coordinate system that anchors dual variational autoencoders (VAEs) trained with cross-modal geometric regularization, preserving archetype structure and spatial neighborhood information while enabling bidirectional translation between modalities. On semi-synthetic CITE-seq data, ARCADIA outperforms existing weak-linkage methods. Applied to independent human tonsil scRNA-seq and CODEX data, ARCADIA reconstructs known tissue architecture and reveals spatially dependent transcriptional programs linking B-cell maturation and T-cell activation or exhaustion to microenvironmental niches. AVAILABILITY AND IMPLEMENTATION: Source code is accessible at https://github.com/azizilab/ARCADIA_public. Reproducibility scripts and data are available at https://github.com/azizilab/arcadia_reproducibility.

Proteomics↗

PubMed2026

Single-cell transcriptional profiling identifies the swimming crab Portunus trituberculatus in response to bacterial infection.

Crustaceans rely entirely on innate immunity, yet the cellular composition, functional specialization, and pathogen-induced remodeling of their immune system remain poorly resolved. Here, we generated a high-resolution single-cell transcriptomic atlas of hemocytes from the swimming crab Portunus trituberculatus following Vibrio parahaemolyticus infection using 10× Genomics scRNA-seq. Seven putatively distinct hemocyte clusters were identified, including granulocytes, semigranular hemocytes, prohemocytes, unresolved hemocytes, hyalinocyte-like hemocytes, biosynthetically active secretory hemocytes, and regulatory hemocytes. Although the overall cellular composition remained relatively stable after infection, hemocytes exhibited pronounced cluster-specific transcriptional reprogramming involving Toll/NF-κB signaling, antimicrobial peptide synthesis and metabolic rewiring. By integrating single-cell and bulk transcriptomes, we identified multiple anti-lipopolysaccharide factors (ALFs) as key secretory effectors and experimentally validated their antibacterial activities. FITC-based bacterial engulfment assays and RNA-seq of sorted phagocytes demonstrated that phagocytic capability was shared across multiple hemocyte clusters. Notably, the immunoglobulin superfamily receptor DSCAM displayed extensive alternative splicing and strong infection-induced activation in unresolved hemocytes. Immune-training experiments showed that prior bacterial exposure was associated with altered DSCAM expression and reduced early cumulative mortality upon secondary challenge, suggesting a memory-like immune phenotype. These findings provide a foundational framework for understanding crustacean immunity and advancing disease-resistant breeding in aquaculture.

Antimicrobial peptides↗

PubMed2026

Prophylactic Inhaled Pattern Recognition Receptor Agonists Reprogram Lung Epithelial Response and Prevent Type 2 Allergic Inflammation.

Prophylactic inhalation of the synergistic agents ODN M362 and Pam2CSK4 ("Pam2ODN") protects mice against allergic lung disease, including allergic inflammation caused by house dust mite (HDM). By preventing sensitization, Pam2ODN reduces HDM-induced eosinophilic and lymphocytic inflammation. How Pam2ODN affects interactions among lung epithelial cells, dendritic cells, and T cells to prevent eosinophilic lung inflammation remains unclear. In the present study, we show that a single inhaled dose of Pam2ODN before HDM sensitization reduces airway Th2 polarization without affecting Th1 or Treg responses. Furthermore, Pam2ODN pretreatment inhibits the recruitment of lung monocyte-derived dendritic cells (moDCs) and conventional Type 2 dendritic cells (DC2s), while preventing the HDM-induced decrease in conventional Type 1 dendritic cells (DC1s). Bulk RNA-seq of the whole lung reveals that Pam2ODN pretreatment restricts the expression of proinflammatory transcripts induced by HDM sensitization. This tolerogenic effect is also reflected at the single-cell level in lung epithelial cells, where proinflammatory transcripts, pathways, and chromatin accessibility are inhibited. These results indicate that Pam2ODN reprograms lung epithelial cells to attenuate allergen-induced Th2-promoting cytokines and DCs while maintaining the population of protective DC1s. These findings suggest a strategy to mitigate chronic allergic lung diseases.

Animals↗

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↗

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↗

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

Malignant epithelial states drive immune dysfunction in ampulla of Vater carcinoma.

BACKGROUND: Ampulla of Vater (AoV) carcinoma is a rare malignancy arising at the junction of intestinal and pancreatobiliary epithelium. Its heterogeneous clinical behavior and histological diversity have hindered therapeutic advances, and the cellular basis of this heterogeneity remains unclear. We aimed to construct a single-cell transcriptomic atlas of AoV carcinoma, with a focus on identifying epithelial subtypes and their interactions with the tumor microenvironment (TME). METHODS: We performed single-cell RNA sequencing on eight primary AoV tumors and four matched normal tissues. Comprehensive clustering and transcriptomic analyses identified cell-type composition, epithelial heterogeneity, and tumor-immune interactions. Findings were validated using deconvolution of bulk RNA-seq data from 62 AoV carcinoma patients. Results Malignant epithelial cells were categorized into four distinct subtypes: Int-Wnt, PB-KRAS, Int-Hypoxia, and Cycling stage. PB-KRAS cells exhibited stem-like transcriptional programs and high genomic instability. Deconvolution analysis of bulk RNA-seq data from the independent AoV cohort revealed that enrichment of the PB-KRAS subtype correlated with tumor recurrence and poor survival. Our immune profiling analysis discovered a significant association between PB-KRAS subtype and GZMK+ CD8+ T cells, which are in a pre-dysfunctional state, alongside SPP1+ macrophages exhibiting immunosuppressive traits. Spatial transcriptome data further supports the immunosuppressive natures of TME around PB-KRAS subtype malignant epithelial cells in AoV carcinoma. CONCLUSIONS: Our study presents a single-cell atlas of AoV carcinoma, highlighting the molecular diversity of malignant epithelium and its association with the immune microenvironment. The PB-KRAS subtype emerges as a stem-like, immunosuppressive tumor state associated with poor prognosis, providing insights for future therapeutic targeting.

Ampulla of Vater carcinoma↗

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

FANCI promotes esophageal squamous cell carcinoma progression and cell cycle regulation and interacts with FANCD2.

BACKGROUND: Esophageal squamous cell carcinoma (ESCC) is an aggressive malignancy with poor clinical outcomes, and reliable molecular biomarkers and therapeutic targets remain limited. Fanconi anemia group I protein (FANCI) is a core component of the Fanconi anemia (FA) pathway, but its expression pattern, clinical significance, and functional role in ESCC have not been comprehensively defined. This study aimed to investigate FANCI expression and prognostic value in ESCC, assess its effects on malignant cellular phenotypes and tumor growth, and explore its potential mechanistic relationship with Fanconi anemia group D2 protein (FANCD2) and cell-cycle regulation. METHODS: Multi-cohort analyses were performed using The Cancer Genome Atlas (TCGA) and Gene Expression Omnibus (GEO) datasets, together with ESCC single-cell RNA sequencing (RNA-seq) data. FANCI functions were assessed by bidirectional gain- and loss-of-function experiments in vitro (proliferation, colony formation, migration, invasion, apoptosis, and cell-cycle assays) and by xenograft models in vivo. Mechanistic studies included protein-protein interaction (PPI) analyses, co-immunoprecipitation (Co-IP), and immunofluorescence (IF) colocalization. RESULTS: FANCI was consistently upregulated in ESCC across bulk transcriptomic datasets and was further supported by quantitative polymerase chain reaction (qPCR), Western blotting, and immunohistochemistry (IHC). FANCI discriminated ESCC from normal tissues in TCGA-ESCC and was independently validated in GSE53624 [area under the curve (AUC) =0.940 and 0.975, respectively]. FANCI was associated with poorer overall survival (OS) and shorter disease-free interval (DFI), and these findings were validated in an independent GEO cohort. Functionally, FANCI promoted ESCC cell proliferation, migration, and invasion, while inhibiting apoptosis; FANCI knockdown suppressed tumor growth in vivo and induced G2/M cell-cycle arrest. Mechanistically, FANCI physically interacted with FANCD2, colocalized with FANCD2 in the nucleus, and was associated with altered FANCD2 protein abundance, consistent with cell-cycle and DNA repair-related programs. Single-cell analysis indicated that FANCI was enriched in epithelial cells and associated with higher activity of malignant functional programs. In TCGA-ESCC, FANCI-high tumors showed distinct mutation profiles, a trend toward increased tumor mutation burden (TMB), and altered immune-associated signatures. CONCLUSIONS: FANCI is upregulated in ESCC and is associated with diagnostic and prognostic value. It promotes malignant phenotypes and tumor growth, potentially through a FANCI-FANCD2-linked cell-cycle/DNA repair program, supporting FANCI as a candidate biomarker and therapeutic target in ESCC.

Esophageal squamous cell carcinoma (ESCC)↗