SEARCH · PubMed Health

Results for “scRNA-seq data”

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 55 records · Page 3Linked to original sources

PubMed2026

Spatial Transcriptomics Identifies Characteristic Immunological Niches in Atopic Dermatitis.

BACKGROUND: Atopic dermatitis (AD) is primarily driven by a Type 2 immune response, with T helper (TH2) cells producing IL-4 and IL-13, thereby promoting inflammation, itch, and a compromised skin barrier. Yet, the spatial organization of pathogenic immune cells and their interactions with stromal and epithelial compartments in human AD skin remain incompletely understood. METHODS: We performed 10× Genomics Visium spatial transcriptomics on FFPE skin biopsies from patients with AD (n = 6), psoriasis (n = 2), and healthy controls (n = 5). Data were integrated with AD single-cell RNA sequencing (scRNA-seq) datasets and complemented by imaging mass cytometry (IMC) and multiplex immunofluorescence (IF) to validate the spatial localization of immune cells. Cell-cell communication analysis revealed putative signaling interactions within immune niches. RESULTS: Spatial clustering resolved tissue compartments and demonstrated transcriptional dysregulation in keratinocytes in AD and psoriasis. AD lesions showed a conserved spatial organization of immune aggregates within the superficial dermis. Integration of scRNA-seq signatures revealed spatially organized co-localization of T cells and mature migratory dendritic cells (mmDCs). We developed a ring-based neighborhood analysis to characterize the cellular organization of the immune-stromal niches, revealing T cell-enriched regions surrounded by inflammatory fibroblasts and activated keratinocytes. Intercellular communication analysis further identified putative signaling within mmDC-T cell niches that may promote pathogenic T cell recruitment and activation. Application of tertiary lymphoid structure (TLS) signatures indicated the presence of TLS-like regions. IMC and IF validated the close spatial proximity between activated TH2 cells and mmDCs. CONCLUSION: AD lesions contain spatially organized TLS-like immune niches at the dermal-epidermal junction, characterized by the close association of T cells and mmDCs and coordinated interactions with surrounding stromal and epithelial compartments. These mmDC-T cell niches may represent potential targets for future therapeutic strategies aimed at disrupting persistent local inflammatory pathways and improving long-term disease control.

atopic dermatitis↗

PubMed2026

A novel lactylation-related gene signature deciphers the immunosuppressive microenvironment and stratifies precision therapy in colorectal cancer.

BACKGROUND: Colorectal cancer (CRC) remains a leading cause of cancer mortality, largely due to the heterogeneity of the tumor microenvironment (TME) and the limited efficacy of immunotherapy in microsatellite stable (MSS) tumors. Histone lactylation, a post-translational modification derived from the Warburg effect, serves as a critical bridge linking metabolic reprogramming to gene regulation and immune evasion; however, its specific prognostic value and clinical implications in CRC remain to be fully elucidated. METHODS: In this study, we systematically analyzed transcriptome profiling data from The Cancer Genome Atlas (TCGA) and Gene Expression Omnibus (GEO) cohorts, supplemented by single-cell RNA sequencing (scRNA-seq) analysis and Human Protein Atlas (HPA) protein-level validation. By integrating univariate Cox regression, Least Absolute Shrinkage and Selection Operator (LASSO) analysis, and multivariate Cox regression, we constructed a novel lactylation-related gene (LRG) risk signature. We extensively evaluated the association between this risk signature and patient prognosis, immune infiltration patterns, somatic mutations, and therapeutic sensitivity. RESULTS: A robust 9-gene prognostic signature (DHRS7, SPR, MBD2, RBM17, CSRP2, S100A4, TMSB4X, TKT, COPS4) was identified and corroborated at the protein level. Patients with high risk scores exhibited significantly worse overall survival (OS) across the training and two independent validation cohorts. Immunogenomic and scRNA-seq analyses revealed that high-risk tumors were characterized by an immunosuppressive and stromal-dense microenvironment-with stromal cells exhibiting the highest lactylation risk scores-enriched with regulatory T cells (Tregs), and frequently harbored PIK3CA mutations. Differential expression analysis indicated that this immune exclusion is structurally maintained by enriched extracellular matrix (ECM) organization and TGF-β signaling. Conversely, low-risk tumors displayed an inflamed phenotype with active antitumor immunity. Pharmacogenomic prediction identified distinct therapeutic stratifications: low-risk patients exhibited significant sensitivity to standard chemotherapeutics (fluorouracil, oxaliplatin) and EGFR/HER2 inhibitors (e.g., lapatinib, erlotinib). In contrast, high-risk patients showed specific vulnerabilities to novel targeted agents, including PI3K pathway inhibitors (TG-100-115, XL765), microenvironment-modulating agents (sildenafil, GANT-61), and epigenetic inhibitors (UNC0638). CONCLUSION: We established a novel lactylation-related risk signature that effectively stratifies CRC patients by prognosis and TME characteristics. By elucidating the crosstalk between metabolic dysregulation, stromal barriers, and immune exclusion, this study provides potential biomarkers and stratified therapeutic strategies-ranging from standard chemotherapy to targeted metabolic and stromal interventions-to optimize precision medicine for CRC patients.

Colorectal cancer↗

PubMed2025

scnanoseq: an nf-core pipeline for Oxford Nanopore single-cell RNA-sequencing.

MOTIVATION: Recent advancements in long-read single-cell RNA sequencing (scRNA-seq) have facilitated the quantification of full-length transcripts and isoforms at the single-cell level. Historically, long-read data would need to be complemented with short-read single-cell data in order to overcome the higher sequencing errors to correctly identify cellular barcodes and unique molecular identifiers. Improvements in Oxford Nanopore sequencing, and development of novel computational methods have removed this requirement. Though these methods now exist, the limited availability of modular and portable workflows remains a challenge. RESULTS: Here, we present, nf-core/scnanoseq, a secondary analysis pipeline for long-read single-cell and single-nuclei RNA that delivers gene and transcript-level quantification. The scnanoseq pipeline is implemented using Nextflow and is built upon the nf-core framework, enabling portability across computational environments, scalability and reproducibility of results across pipeline runs. The nf-core/scnanoseq workflow follows best practices for analyzing single-cell and single-nuclei data, performing barcode detection and correction, genome and transcriptome read alignment, unique molecular identifier deduplication, gene and transcript quantification, and extensive quality control reporting. AVAILABILITY AND IMPLEMENTATION: The source code, and detailed documentation are freely available at https://github.com/nf-core/scnanoseq and https://nf-co.re/scnanoseq under the MIT License. Documentation for the version of nf-core/scnanoseq used for this paper, including default parameters and descriptions of output files are available at https://nf-co.re/scnanoseq/1.1.0.

Single-Cell Analysis↗

PubMed2025

Spatial transcriptomics-aided localization for single-cell transcriptomics with STALocator.

Single-cell RNA-sequencing (scRNA-seq) techniques can measure gene expression at single-cell resolution but lack spatial information. Spatial transcriptomics (ST) techniques simultaneously provide gene expression data and spatial information. However, the data quality of the spatial resolution or gene coverage is still much lower than the quality of the single-cell transcriptomics data. To this end, we develop a ST-Aided Locator for single-cell transcriptomics (STALocator) to localize single cells to corresponding ST data. Applications on simulated data showed that STALocator performed better than other localization methods. When applied to the human brain and squamous cell carcinoma data, STALocator could robustly reconstruct the relative spatial organization of critical cell populations. Moreover, STALocator could enhance gene expression patterns for Slide-seqV2 data and predict genome-wide gene expression data for fluorescence in situ hybridization (FISH) and Xenium data, leading to the identification of more spatially variable genes and more biologically relevant Gene Ontology (GO) terms compared with the raw data. A record of this paper's transparent peer review process is included in the supplemental information.

Single-Cell Analysis↗

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

Integrated single-cell RNA sequencing and mendelian randomization analysis identifies causal immune-related driver genes in the heart failure inflammatory microenvironment.

BACKGROUND: Heart failure (HF) is a major global cause of cardiovascular death and disability. Chronic inflammation and immune dysregulation are critical in its development. The cardiac immune microenvironment, especially macrophages, drives HF progression, yet its molecular mechanisms and prognostic impact are not fully clear. This study aimed to identify causal immune-related driver genes in the HF inflammatory microenvironment. METHODS: We combined single-cell RNA sequencing (scRNA-seq) and Mendelian randomization (MR) to study how the inflammatory immune microenvironment affects HF risk. Using two public scRNA-seq datasets, we identified differentially expressed genes (DEGs) in HF heart tissues and selected 489 candidate genes. Causal relationships between these genes and HF were tested using expression quantitative trait loci (eQTL) data and HF genome-wide association study (GWAS) summary statistics. RESULTS: MR analysis showed that 65 genes were causally linked to HF risk. These genes were enriched in pathways related to cardiomyopathy, leukocyte migration, natural killer (NK) cell cytotoxicity, neutrophil extracellular traps, and NF-κB signaling. HF hearts displayed increased levels of macrophages, T cells, B cells, lymphoid cells, and mast cells, while neutrophils were reduced. CONCLUSIONS: Our integrated analysis reveals the central role of the cardiac inflammatory immune microenvironment in HF and identifies 65 key genes causally associated with HF susceptibility. These genes influence specific immune pathways and cell infiltration, shaping HF progression, and provide a basis for developing new biomarkers and immune-targeted therapies.

Heart failure (HF)↗

PubMed2026

Integrated bulk and single-cell RNA sequencing reveals a prognostic neuro-mimicry signature in papillary thyroid carcinoma.

BACKGROUND: Cancer cells can acquire neuron-like characteristics ("neural mimicry") to promote progression. However, the role of specific ion channel genes in Papillary Thyroid Carcinoma (PTC) and their clinical significance remains unclear. METHODS: We included transcriptomic data from 521 PTC patients in the TCGA cohort. A neuron-specific gene set was used to screen for potential targets. We constructed a prognostic model using LASSO logistic regression. To verify the cellular origin of the signature, we performed single-cell RNA sequencing (scRNA-seq) analysis on the GSE184362 dataset. RESULTS: We established an 8-gene signature involving KCNN4, KCNN1, KCNT2, SNAP25, KCNK16, GABRG1, GABRG2, and GABRB2. The model demonstrated good predictive performance for lymph node metastasis, with an AUC of 0.721 (95% CI 0.677-0.765). Single-cell analysis of seven integrated tumor samples (N = 65,744 cells) confirmed that GABRB2 was specifically enriched in malignant thyrocytes (EPCAM+/KRT18+) at 200-fold higher detection rates than immune cells (20.0% vs. 0.1%, P ≈ 0), supporting tumor-intrinsic neural mimicry. High-risk patients showed immunosuppressive features with altered immune cell infiltration patterns. CONCLUSION: This study identifies a malignant cell-intrinsic signature for predicting PTC prognosis. Validated by single-cell data, our findings suggest that targeting ion channels may represent a potential therapeutic strategy for modulating neuro-immune interactions in thyroid cancer, pending experimental validation.

GABRB2↗

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

Single-cell transcriptomics on FFPE placenta: A novel method for comprehensive exploration of an entire placental section.

INTRODUCTION: The placenta's complex cellular diversity challenges traditional transcriptomic analyses. Single-cell RNA sequencing (scRNA-seq) offers breakthrough capabilities by enabling transcriptome profiling at the single-cell level. However, traditional scRNA-seq relies on fresh or frozen samples, which present practical storage and quality challenges. Applying scRNA-seq to Formalin-Fixed, Paraffin-Embedded (FFPE) placentas could harness archived samples for clinical insights. METHODS: We used 10x Genomics Flex technology to analyze 8 non-pathological placentas ranging from 21 + 6 weeks of gestation (WoG) to 39 + 4 WoG. RESULTS: Our approach identifies diverse cell populations and allows us to discern maternal from fetal cells. Despite sample size limitations, the method yields comparable data to prior fresh/frozen tissue studies and we complete these data by integrating new molecular markers. The potential to correlate single-cell results with histopathology enables us to conduct an in-depth analysis across entire placental sections by concurrently addressing both fetal and maternal cells. We could thus confirm molecular markers like KRT5/6 using immunohistochemistry by revisiting the slide. DISCUSSION: This innovation could aid in understanding focal anomalies observed on standard histology slides, thereby enhancing traditional histopathological assessments. Given its practicality, integrating our method into routine practice is both feasible and promising.

Differentially expressed genes (DEG)↗

PubMed2026

Multi-omics and spatial transcriptomics reveal that S100A10 drives CD8+ T-cell exhaustion and immune evasion in hepatocellular carcinoma through cPLA2-5-LOX-mediated arachidonic acid metabolism and ferroptosis.

Immune evasion in hepatocellular carcinoma (HCC) represents a major biological barrier limiting the efficacy of immunotherapy, yet its molecular basis remains incompletely understood. Increasing evidence indicates that tumor metabolic reprogramming and ferroptosis-related signaling play critical roles in shaping an immunosuppressive tumor microenvironment (TME); however, the specific regulatory factors involved remain unclear. This study aims to systematically elucidate the functional role of S100 calcium-binding protein A10 (S100A10) in immune evasion in HCC, with a particular focus on the molecular mechanisms by which S100A10 regulates CD8+ T-cell exhaustion through arachidonic acid (AA) metabolism and ferroptosis, as well as its potential therapeutic implications. To this end, data from The Cancer Genome Atlas Liver Hepatocellular Carcinoma (TCGA-LIHC) cohort are integrated to analyze the expression patterns of S100A10, its prognostic value, and its association with the immune microenvironment. S100A10 overexpression and knockout models are established in HCCLM3 and MHCC97L cell lines, and S100A10-mediated metabolic pathway reprogramming is characterized using transcriptomic profiling, untargeted metabolomics, and ferroptosis-related functional assays. In parallel, single-cell RNA sequencing (scRNA-seq) and spatial transcriptomics are employed to delineate the cell-type specificity and spatial distribution of S100A10. Furthermore, human CD8+ T-cell co-culture systems and orthotopic mouse HCC models are used to evaluate the impact of S100A10 on immune function and responsiveness to anti-programmed cell death protein 1 (anti-PD-1) therapy. The results demonstrate that S100A10 is significantly upregulated in HCC and is closely associated with poor prognosis and an immunosuppressive state. Mechanistically, S100A10 activates cytosolic phospholipase A2-arachidonate 5-lipoxygenase (cPLA2-5-LOX)-mediated AA oxidative metabolism, leading to the accumulation of lipid peroxidation products and ferroptosis-associated signals, thereby driving CD8+ T-cell exhaustion and promoting immune evasion. Significantly, inhibition of S100A10 reshapes the tumor immune microenvironment (TIME) and enhances the therapeutic efficacy of anti-PD-1 treatment. Collectively, these findings identify S100A10 as a critical regulator of metabolic-immune coupling in HCC and provide a theoretical basis for combinatorial strategies targeting metabolism and immunotherapy.

Arachidonic acid metabolism↗

PubMed2026

Integrative Multi-Omics Analysis Identifies Thrombosis-Associated Molecular Features Linked to Germline Susceptibility and Immune Cell Communication in Gastric Cancer.

Emerging evidence indicates that coagulation-related molecular programs are associated with thrombosis, tumor progression, and molecular dysregulation in gastric cancer (GC). However, thrombosis-associated molecular features in GC and their potential links to inherited susceptibility remain insufficiently understood. Integrated analyses of transcriptomic data from The Cancer Genome Atlas (TCGA) and Gene Expression Omnibus (GEO) datasets were performed to identify thrombosis-associated genes and establish a machine learning-based prognostic signature. Genome-wide association study (GWAS), expression quantitative trait loci (eQTL), transcriptome-wide association study (TWAS), and Mendelian randomization (MR) analyses were conducted to investigate susceptibility-associated transcriptional programs in GC. Functional assays were used to evaluate candidate genes associated with malignant phenotypes. Single-cell RNA sequencing (scRNA-seq) and cell-cell communication analyses were further performed to characterize cell-type-specific expression patterns and potential intercellular interactions. A total of 22 differentially expressed thrombosis-associated genes were identified, and a prognostic signature comprising 14 genes was established. The signature stratified patients into high- and low-risk groups and showed prognostic performance in both the training and validation cohorts. Integrative GWAS, eQTL, and TWAS analyses identified susceptibility-associated transcriptional programs that were positively correlated with the thrombosis-associated risk score. Silencing ACTN2 and CRYAB significantly reduced GC cell migration and invasion. scRNA-seq analysis revealed relatively high CRYAB expression in neutrophils, and CellChat analysis suggested potential neutrophil-B cell interactions involving COLLAGEN-related signaling. This integrative multi-omics study identified a thrombosis-associated molecular signature linked to prognosis and germline susceptibility-associated transcriptional programs in GC. ACTN2 and CRYAB may represent candidate genes associated with GC cell migration and invasion, while single-cell analysis suggested potential immune-related communication features.

Humans↗

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↗

PubMed2024

An isoform-resolution transcriptomic atlas of colorectal cancer from long-read single-cell sequencing.

Colorectal cancer (CRC) ranks as the second leading cause of cancer deaths globally. In recent years, short-read single-cell RNA sequencing (scRNA-seq) has been instrumental in deciphering tumor heterogeneities. However, these studies only enable gene-level quantification but neglect alterations in transcript structures arising from alternative end processing or splicing. In this study, we integrated short- and long-read scRNA-seq of CRC samples to build an isoform-resolution CRC transcriptomic atlas. We identified 394 dysregulated transcript structures in tumor epithelial cells, including 299 resulting from various combinations of splicing events. Second, we characterized genes and isoforms associated with epithelial lineages and subpopulations exhibiting distinct prognoses. Among 31,935 isoforms with novel junctions, 330 were supported by The Cancer Genome Atlas RNA-seq and mass spectrometry data. Finally, we built an algorithm that integrated novel peptides derived from open reading frames of recurrent tumor-specific transcripts with mass spectrometry data and identified recurring neoepitopes that may aid the development of cancer vaccines.

Colorectal Neoplasms↗

PubMed2026

Distinct cochlear cell types associated with genetic susceptibility to sensory and metabolic hearing loss in older adults.

Hearing loss is a heterogeneous condition that can be classified into different subtypes with diverse genetic and cellular components. To investigate the cochlear cell types underlying the genetic basis of sensory and metabolic components of age-related hearing loss (ARHL), we integrated human genome-wide association study data with mouse cochlear single-cell RNA sequencing data using the single-cell disease relevance score tool. These analyses revealed that genes associated with the sensory component of ARHL were most highly expressed in hair cells, while genes associated with the metabolic component of ARHL were most highly expressed in spiral ganglion neurons. We also investigated whether ARHL-associated gene expression patterns differed across subpopulations of the same cell type. Sensory hearing loss-associated genes showed differential expression across supporting cell subpopulations in younger mice, whereas metabolic hearing loss-associated genes exhibited differences across intermediate cell subpopulations of the stria vascularis in older mice. These findings provide evidence for the role of distinct genetic and cellular risk profiles for different ARHL subtypes, suggesting that prevention and therapeutic strategies may require targeting specific cell populations at different life stages.

ARHL↗

PubMed2026

Integrating single-cell transcriptomics to construct an oncogene-driven prognostic model and elucidate metabolic-immune crosstalk in hepatocellular carcinoma.

Hepatocellular carcinoma (HCC) is a leading cause of cancer-related deaths, its progression and treatment heterogeneity are mainly influenced by driver gene and tumor micro-environment (TME) interactions. Nevertheless, the mechanisms of this process at the single-cell level remain unclear. This study integrated TCGA and multi-center single-cell transcriptome data to identify a 575 genes HCC-specific core set, developing a single-cell "oncogene scoring" system to quantify individual carcinogenic activity. This score is significantly elevated in malignant and proliferative T cells and is closely associated with metabolic reprogramming, aberrant cell‒cell communication, and immunosuppressive phenotypes. Based on these characteristics, we constructed a machine learning-based Random Survival Forest (RSF) prognostic model validated in multiple independent cohorts, which classifies patients into distinct risk subtypes. The high-risk group exhibits genomic instability, increased tumor stemness, and immune evasion, while the low-risk group was more sensitive to drugs such as sorafenib. This study highlights the potential pathways by which high oncogenic activity is associated with HCC progression, suggesting a profound link with single-cell metabolic‒immune crosstalk. The constructed RSF model offers a promising computational framework for risk stratification and provides hypothesis-generating insights that may inform future personalized treatment strategies for HCC patients.

Hepatocellular carcinoma↗

PubMed2026

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

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

Acral melanoma↗

PubMed2025

dbscATAC: a resource of single-cell super-enhancers/enhancers and gene markers derived from scATAC-seq data.

MOTIVATION: scATAC-seq enables high-resolution mapping of cis-regulatory elements. It has been widely applied to uncover cell-type-specific regulatory networks and complement scRNA-seq analysis in numerous studies. However, a large number of datasets generated by scATAC-seq remain underutilized due to limited exploration of super-enhancers/typical enhancers and gene markers. A comprehensive resource enabling cell-type-specific annotation of cis-regulatory elements and their dynamic enhancer-gene linkages remains an urgent unmet need for scATAC-seq. RESULTS: We present dbscATAC, a specialized single-cell database for annotating super-enhancers, gene markers, and enhancer-gene interactions derived from scATAC-seq data. Using improved machine learning algorithms, we identified 213 835 super-enhancers across 520 tissue/cell types from three species, as well as 347 484 gene markers, 13 470 526 enhancers, and 10 402 346 enhancer-gene interactions derived from 1 668 076 single cells spanning 1028 tissue/cell types in 13 species. An easy-to-use online platform with multiple analytic modules and hierarchical query options was developed for searching, browsing and visualizing single-cell super-enhancers, enhancers, and gene markers. dbscATAC provides a comprehensive resource to facilitate the exploration of enhancer landscapes, gene regulation, and cell-type-specific characteristics in single-cell epigenomics. AVAILABILITY AND IMPLEMENTATION: The database with all the super-enhancer/enhancer annotation data is available at http://singlecelldb.com/dbscATAC/index.php. And the source code of dbscATAC for prediction of SEs, enhancers, and gene markers are available at https://github.com/EvansGao/dbscATAC. The source code, tissue/cell type description, and data summary can be downloaded at DOI: 10.6084/m9.figshare.28706414.scATAC-seq, Database, Super-enhancers/enhancers, Gene markers.

Enhancer Elements, Genetic↗

PubMed2026

Amaranth: enhanced single-cell transcript assembly via discriminative modelling of UMI reads and internal reads.

MOTIVATION: Single-cell RNA sequencing (scRNA-seq) has transformed transcriptome profiling at cellular resolution, yet accurate reconstruction of full-length transcripts for individual cells remains a central challenge. Emerging scRNA-seq protocols can produce reads that span entire transcripts, enabling isoform-level expression analysis. For example, Smart-seq protocols combine unique molecular identifier (UMI)-linked reads that index and stitch together multiple reads from the same molecule, with internal reads filling coverage gaps. We demonstrate that these read types exhibit markedly different biological and statistical properties in strandness, 5'/3' coverage bias, and genomic locality. Existing assemblers fail to leverage these distinctions, yielding suboptimal assembly. RESULTS: We developed Amaranth, a novel single-cell assembler that discriminatively models UMI and internal reads. Amaranth implements heuristics specifically designed to address the distinct biases of UMI-linked and internal reads, enabling accurate strandness assignment for internal reads, reliable splicing graph refinement, and precise transcript start site determination. We also developed Amaranth-meta, which integrates information across cells to enhance individual cell assemblies. Benchmarked on Smart-seq3 datasets from human HEK293T and mouse fibroblast cells, Amaranth outperformed other state-of-the-art assemblers in assembling individual cells and in meta-assembly. Amaranth advances isoform-level analysis in single-cell transcriptomics, facilitating detailed studies at cellular resolution. AVAILABILITY AND IMPLEMENTATION: Amaranth is implemented in C++ and is freely available at https://github.com/Shao-Group/amaranth under the BSD-3-Clause license. Scripts, documentation, and data for reproducing experiments in this manuscript are available at https://github.com/Shao-Group/amaranth-test.

Single-Cell Gene Expression Analysis↗