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 37 records · Page 2Linked to original sources

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

CeLLTra: aligning cell names with gene expression via a pathway-informed transformer.

MOTIVATION: Single-cell RNA sequencing (scRNA-Seq) technology enables detailed exploration of gene expression at the individual cell level, crucial for annotating cell types and understanding cellular diversity. Traditional methods for cell type annotation often rely on marker genes and manual labeling, posing challenges due to low data quality and incomplete reference datasets. RESULTS: We developed CeLLTra, a novel contrastive learning framework that leverages a Transformer-based model integrating biological pathway information to group genes into super tokens, effectively capturing comprehensive gene expression from scRNA-Seq data. By combining this pathway-informed Transformer with a pretrained domain-specific language model, CeLLTra accurately aligns cell-type annotations with gene expression profiles. Evaluations on a large-scale human scRNA-Seq dataset showed that CeLLTra significantly outperformed state-of-the-art methods in supervised and zero-shot cell-type prediction. Additionally, CeLLTra generalized well to external datasets, improving clustering performance and enabling better characterization of cancerous cell states in tumor-infiltrating myeloid cells from non-small cell lung cancer patients. AVAILABILITY AND IMPLEMENTATION: CeLLTra is freely available on GitHub (https://github.com/WJZheng-group/CeLLTra) and Zenodo (https://doi.org/10.5281/zenodo.17666735). The datasets underlying this article are the following: GSE201333 and GSE127465. All these datasets are publicly available and can be freely accessed on the Gene Expression Omnibus repository.

Humans↗

PubMed2026

Big data analytics for CLEC5A dynamics based on single cell genomics and proteomics reveal its diverse functions in human diseases.

BACKGROUND: CLEC5A (C-type lectin domain family 5 member A) is an innate immune receptor implicated in inflammatory signaling, contributing to hyperinflammatory responses in infections and sterile inflammation. However, CLEC5A dynamics in human diseases remain to be identified. Here, we systematically characterized CLEC5A dynamics in humans across cells, tissues, and disease states, and to explore the functional significance of CLEC5A in macrophage activation based on single-cell genomics. METHODS: With multi-omics (scRNA-seq, proteomics and big data analytics), we analyzed extensive human transcriptomic datasets (>42,000 samples) to profile CLEC5A expression by cell type, tissue, and disease. Single-nucleus RNA-seq (snRNA-seq) from pediatric congenital heart disease and a virtual CLEC5A gene knockout were also performed to characterize CLEC5A dynamics in humans. RESULTS: CLEC5A is highly enriched in innate immune cells, particularly in macrophages and neutrophils. Baseline CLEC5A in most tissues is low, but it is markedly upregulated in inflammatory and infectious diseases. CLEC5A expression has sex-specific differences in certain organs. Single-cell analysis showed that CLEC5A can be considered novel marker of proinflammatory macrophages with elevated cytokine production, antigen presentation, and impaired phagocytosis. Virtual CLEC5A knockout analysis identified coordinated perturbation of immune-regulatory pathways and overlapping genes linking CLEC5A to macrophage activation networks. CONCLUSION: CLEC5A is predominantly expressed in myeloid cells and acts as a key amplifier of inflammation in human diseases. Our findings highlight CLEC5A as a potential biomarker and therapeutic target in myeloid-driven hyperinflammatory conditions, warranting further experimental and translational validation.

Humans↗

PubMed2026

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

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

Spatial Transcriptomics↗

PubMed2026

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

Enhancing pan-cancer spatial transcriptomics at single-cell resolution with stPainter.

Subcellular spatial transcriptomics can resolve tissue architecture at cellular scale, but sparse gene panels and limited detection sensitivity constrain downstream analysis. Existing enhancement methods often require tissue-matched single-cell RNA sequencing (scRNA-seq) references and dataset-specific retraining. Here we show that stPainter, a conditional generative model pretrained on a pan-cancer scRNA-seq atlas, can enhance spatial transcriptomics data without matched references or retraining. Using a latent diffusion architecture guided by Stochastic Differential Equations (SDE), stPainter reconstructs expanded expression profiles from sparse measurements and produces latent representations for clustering and cell-state analysis. When we apply stPainter upon 6 spatial transcriptomics datasets of different cancer types, we demonstrate that our model empowers downstream biological analyses, including fine-grained subpopulation clustering and pathway enrichment. Comparison with spatially resolved proteomics (CODEX) provided independent support for regional agreement between imputed cellular compositions and protein-level tissue organization. These results establish stPainter as a scalable approach for analyzing tumor microenvironments without auxiliary sequencing data.

Spatial Transcriptomics↗

PubMed2026

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

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

CRISPR↗

PubMed2025

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 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

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

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↗