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 19 recordsLinked to original sources

PubMed2026

A contextual activity score (CAS) for inferring ADAR-associated transcriptional activity across RNA-seq, single-cell, and spatial transcriptomics.

BACKGROUND AND OBJECTIVE: Adenosine-to-inosine RNA editing, catalyzed by Adenosine Deaminases Acting on RNA (ADARs), is a widespread modification involved in neural function, immune regulation, and cancer. The Alu Editing Index (AEI) is the standard metric to estimate ADAR activity but requires raw sequencing reads and is poorly suited for single-cell and spatial transcriptomic data. This study aimed to develop an alternative framework for inferring ADAR-associated transcriptional activity from gene expression data across diverse transcriptomic technologies. METHODS: We developed the Contextual Activity Score (CAS), a framework based on transcriptional signatures from ADAR perturbation experiments. Context-specific signatures were generated for human neurons, mouse neurons, and cancer models to infer ADAR1 and ADAR2 activity. CAS was computed from normalized gene expression matrices using regulon-based enrichment analysis. Performance was evaluated by comparing with the Alu Editing Index across bulk RNA sequencing datasets, simulated sequencing depths, and library preparation protocols. RESULTS: CAS showed strong concordance with the Alu Editing Index across multiple datasets, while remaining robust to reduced sequencing depth and different library protocols. Unlike the Alu Editing Index, CAS can be applied to single-cell and spatial transcriptomic data and enables the independent assessment of ADAR2 activity. In cancer and neuronal contexts, CAS captured biologically meaningful variations in ADAR-associated transcriptional activity at sample, cell-type, and spatial levels. CONCLUSION: CAS provides a scalable approach applicable across multiple RNA-seq protocols for estimating ADAR-associated transcriptional activity using gene expression data. This method, implemented in an open-source R package for broad adoption, expands the ability to study ADAR-associated transcriptional activity across transcriptomic modalities where direct editing quantification is challenging, such as single-cell and spatial transcriptomics.

Adenosine Deaminase↗

PubMed2026

A Multi-omics Regulated Cell Death Framework Defines Immune Phenotypes and Guides Precision Therapy in Colorectal Cancer.

Colorectal cancer (CRC) is molecularly and immunologically heterogeneous, contributing to variable treatment response. Because regulated cell death (RCD) intersects with tumor metabolism, immune regulation, and therapeutic susceptibility, we built an RCD-centered framework for CRC stratification. Multi-cohort transcriptomic data were used to infer RCD subtypes with non-negative matrix factorization (NMF) and non-negative least squares (NNLS). Genomic, bulk RNA-seq, single-cell RNA-seq, and spatial transcriptomic datasets were integrated to characterize subtype-associated biology. Machine-learning models were developed for immunotherapy response and survival-risk estimation. Candidate compounds were screened by GDSC2-based drug-sensitivity modeling and molecular docking, and FSTL3 was functionally assessed in vitro. The framework separated CRC samples into two RCD-related phenotypes resembling immune-hot and immune-cold states. RCD1 showed immune activation and higher mutational burden, whereas RCD2 showed immune-suppressed features, intratumoral heterogeneity, and aggressive biology. RCD-associated signatures showed potential for predicting immunotherapy response and survival risk. Dasatinib was prioritized for immune-cold, high-risk tumors, with preliminary evidence supporting its activity in CRC cells, while functional assays suggested a role for FSTL3 in growth, invasion, epithelial-mesenchymal transition, and apoptosis regulation. These findings suggest that RCD-based multi-omics analysis may refine CRC stratification and help generate therapeutic hypotheses.

Colorectal cancer↗

PubMed2026

Integrated Genomic and Proteomic Analysis Reveals T-B Lymphocyte Signatures in the MYCN Driven "Immune Desert" of Specific Neuroblastoma Subtypes.

AIMS: This study aims to systematically dissect how MYCN amplification shapes the immunosuppressive tumor microenvironment (TME) in high-risk neuroblastoma, elucidating key mechanisms underlying immune evasion. METHODS: We performed an integrated multi-omics analysis of bulk RNA-seq (n = 721), single-cell RNA-seq (n = 9), proteomic data (n = 49) and spatial transcriptomics (Visium, with external validation in melanoma). Analyses included unsupervised clustering, cell-cell communication inference, transcriptional regulatory network reconstruction, and spatial proximity assessment to map the immune landscape. RESULTS: A distinct molecular subtype (Class C), defined by MYCN amplification and poor prognosis, exhibited a comprehensive "immune desert" phenotype characterized by low immune scores and minimal leukocyte infiltration. Single-cell analysis confirmed significant depletion of T and B lymphocytes within the Class C TME. Dysregulated transcriptional networks were identified, including upregulation of REL and EOMES in T cells-with EOMES potentially driving exhaustion via regulation of Transient Receptor Potential (TRP) genes, and REL inhibition enhancing cytotoxic function in vitro. A unique immunosuppressive B-cell subset (B7) engaged in enhanced crosstalk with exhausted T cells and harbored a MYC-centered network linked to cell cycle dysregulation and poor survival. Spatial transcriptomics revealed significant proximity between B7-active regions and Treg/exhaustion-enriched areas, externally validated in melanoma. Proteomic data validated elevated REL expression in MYCN-amplified tumors. CONCLUSION: This work delineates the immunosuppressive architecture of MYCN-driven neuroblastoma, revealing novel regulatory nodes within specific lymphocyte compartments. Integrating single-cell, spatial, and proteomic evidence, we propose REL inhibition as a therapeutic candidate, the EOMES/TRP axis as a bioinformatically supported hypothesis, and the B7/MYC hub as a hypothesis supported by transcriptomic and spatial evidence.

Humans↗

PubMed2026

Human Monocytic Models Reveal Genotype-Dependent Inflammatory Programs in VEXAS Syndrome.

OBJECTIVES: VEXAS syndrome is a severe X-linked autoinflammatory disorder caused by somatic mutations in ubiquitin-like modifier activating enzyme 1 (UBA1), with clinical outcomes that vary by UBA1 genotype. We aimed to elucidate genotype-specific inflammatory programs and identify potential therapeutic targets. METHODS: We conducted longitudinal deep phenotyping, including whole-blood RNA sequencing (RNA-seq) and clinical activity assessment. Peripheral blood samples were analyzed by single-cell RNA-seq. Human monocytic cell lines harboring each major UBA1 mutation (p.Met41Val, p.Met41Thr, or p.Met41Leu) were generated and subjected to transcriptomic and functional analyses. RESULTS: Thirteen patients with VEXAS syndrome contributed a total of 79 RNA-seq samples. Among genes upregulated in VEXAS syndrome, RNASE1 showed the strongest correlation with longitudinal disease activity (r = 0.70, FDR < 0.05) and was upregulated in patients' monocytes. In UBA1-mutant monocytic cell lines, genotype-dependent ubiquitination defects were observed in a graded manner (p.Met41Val > p.Met41Thr > p.Met41Leu), even in the absence of exogenous stimuli. These defects were accompanied by unfolded protein response activation, increased pro-inflammatory cytokine production, progressive cell death, and RNASE1 upregulation, all following the same graded pattern, recapitulating patient genotype-phenotype associations. Transcriptomic analyses demonstrated enrichment of pro-inflammatory, interferon, and necroptosis signatures in more severe genotypes. Notably, inhibition of receptor-interacting protein kinase 3 (RIPK3) markedly attenuated all pathological features, including RNASE1 upregulation. CONCLUSIONS: Our UBA1-mutant monocytic cell-line models, representing three distinct genotypes, recapitulate genotype-dependent inflammatory phenotypes that can be modulated by RIPK3 inhibition, providing a translational platform for mechanistic investigation and precision therapy development in VEXAS syndrome.

Journal Article↗

PubMed2026

Machine learning identifies ac4C-related prognostic signature and TUBA1C as therapeutic target in COAD.

To explore the role of N4-acetylcytidine (ac4C)-related genes (acRGs) in colon adenocarcinoma (COAD) and identify reliable prognostic biomarkers and potential therapeutic targets. Multi-source transcriptomic datasets (TCGA-COAD, GSE39582, GSE17536) and single-cell RNA-seq data were analyzed. Ten machine learning algorithms were integrated to construct an acRG-based prognostic signature (acRGBS). Immune microenvironment (TME) and genomic profiling were performed, with in vitro functional experiments validating TUBA1C's role. acRGBS, comprising four hub genes (SARAF, CDC42SE2, TSPYL2, TUBA1C), effectively stratified COAD patients into high- and low-risk groups with distinct survival outcomes and was an independent prognostic factor. High-risk patients exhibited increased genomic instability and immunosuppressive TME, while low-risk patients had favorable immunotherapy response. TUBA1C was overexpressed in COAD cells, and its knockdown inhibited proliferation/migration and induced apoptosis. The acRGBS is a robust prognostic tool for COAD, and TUBA1C serves as a candidate therapeutic target, providing new insights for personalized COAD management.

Humans↗

PubMed2026

Identification of CD55 as a downstream factor of EP4 receptor signaling in colorectal cancer cells.

Prostaglandin E2 (PGE2) signaling through the E-type prostanoid 4 (EP4) receptor has been implicated in the pathophysiology of colorectal cancer (CRC). We herein identified decay-accelerating factor, also known as CD55, as a novel CRC-associated downstream factor of the EP4 receptor. The integration of transcriptomic profiling of PGE2-stimulated HCA-7 human colon cancer cells with analyses of cancer genomic databases predicted CD55 as a potential EP4 receptor-regulated target. Inhibitor-based experiments showed the induction of CD55 after a PGE2 stimulation required the EP4 receptor and Gi protein in HCA-7 cells, whereas protein kinase A signaling was dispensable. In combination with a toxicogenomic database analysis, p38 mitogen-activated protein kinase (MAPK) was identified as the predominant effector connecting the EP4 receptor to CD55 upregulation. A single-cell RNA-seq re-analysis of human CRC tissues revealed CD55 upregulation and p38 MAPK-related gene set enrichment in epithelial cells expressing the EP4 receptor, suggesting that this induction mechanism may operate in a subset of epithelial cells in clinical specimens. Collectively, these results delineate a PGE2/EP4 receptor/Gi protein/p38 MAPK signaling axis that induces CD55 expression in HCA-7 cells and epithelial tumor cells, provide new mechanistic clues for understanding the regulation of complement regulatory molecule CD55 expression by prostaglandin signaling.

Humans↗

PubMed2026

Single-cell multiomics reveals exosome-mediated reprogramming and clonotypic remodeling of T cells in triple-negative breast cancer.

Triple-negative breast cancer (TNBC) is an aggressive and immunogenic subtype lacking targeted therapies. While tumor-derived exosomes are known to modulate immune function, their direct impact on human T cell plasticity and antigen specificity remains poorly defined. Here, we conducted a comprehensive single-cell multiomic analysis of primary human T cells exposed to exosomes derived from 17 genomically diverse TNBC cell lines and 35 patient samples. Integrating single-cell RNA-seq, V(D)J sequencing, non-coding RNA profiling, bulk and single-cell cytokine analyses, we uncovered conserved and subtype-specific immunomodulatory programs induced by TNBC exosomes. Exosome-treated T cells displayed skewing toward regulatory and dysfunctional phenotypes, including Th17-like, Treg, and PD-1⁺/PD-L1⁺ Tfh cells. Functional profiling revealed suppression of early activation markers and cytokine responses, alongside selective preservation of cytotoxic features in γδ T and NKT subsets. Transcriptomic and miRNA network analyses demonstrated widespread downregulation of immune effector genes (e.g., HBEGF and TNFSF9) mediated by exosome-delivered regulatory miRNAs (has-miR-98-5p). Notably, exosome-stimulated T cells displayed distinct clonotypic expansions, characterized by the emergence of five tumor-specific γδ TCR clonotypes and 30 unique αβ TCR CDR3 sequences that were absent in mock-treated controls, underscoring the role of exosomes in shaping TCR repertoire dynamics.

Humans↗

PubMed2026

Clinical translation of senescence-related pan-cancer multi-omics: tools for assessment and immunotherapy prediction.

Cellular senescence (CS) exerts dual roles in tumorigenesis, yet its pan-cancer molecular characteristics and clinical value remain unclear, hindering its translation to oncology and personalized therapy. To address the lack of specific and universal tools for senescence assessment and immunotherapy response prediction, this study systematically analyzed 1259 CS-related genes from the CellAge database across 31 cancer types by integrating multi-omics data, including bulk RNA-seq, single-cell/spatial transcriptomics, and CRISPR screening. We developed a rank-based algorithm SenScoreR (publicly available at https://gxhub.shinyapps.io/SenScoreR/ ) for senescence quantification, validated with 10 independent datasets, and constructed a machine learning-based predictive model CS.Sig for immunotherapy response. Results showed that tumors had significantly lower Rank-based Senescence Score (RSS) than normal tissues across 31 cancers (average diagnostic AUC = 0.895), with low RSS linked to poor survival; high RSS correlated with reduced genomic instability, enriched CD8⁺ T/NK cell/macrophage infiltration, upregulated PD-L1 expression, and elevated immune cytolytic activity. CS.Sig demonstrated robust performance in predicting ICI response (AUC = 0.716 across 10 cohorts), outperforming 13 existing signatures, while CRISPR screening identified 17 senescence-related targets (e.g., CEP55, PPP1CC) whose knockout enhanced anti-tumor immunity. Our findings clarify CS's role in maintaining tumor genomic stability and shaping immune microenvironments, and the developed SenScoreR, CS.Sig, and identified targets bridge basic CS research with clinical oncology, providing a translational resource and hypothesis basis for future experimental and clinical validation.

Journal Article↗

PubMed2026

PPRC1 is a prognostic biomarker and key regulator of mitochondrial oxidative phosphorylation in multiple myeloma.

BACKGROUND: Multiple myeloma (MM) remains an incurable haematological malignancy, underscoring the need for novel prognostic biomarkers and therapeutic targets. This study aimed to investigate the clinical and biological significance of peroxisome proliferator-activated receptor gamma coactivator-related protein 1 (PPRC1) in MM. METHODS: Expression and clinical data were obtained from public databases and an independent local cohort. Kaplan-Meier and Cox regression analyses were performed to evaluate prognostic value. Differential expression analysis, pathway enrichment analysis and single-cell RNA-seq data analysis were used to explore biological functions. PPRC1 was silenced in MM cell lines using siRNA to assess its effects on cell survival and oxidative phosphorylation. RESULTS: PPRC1 was significantly upregulated in MM and was associated with advanced disease stage and poor overall survival. Multivariate Cox analysis identified PPRC1 as an independent prognostic factor. A nomogram incorporating PPRC1 and revised-ISS improved survival prediction. Functional analyses revealed that PPRC1 was positively correlated with oxidative phosphorylation and oncogenic signalling pathways. A potential connection between PPRC1 expression and immune cell infiltration was observed. PPRC1 knockdown inhibited cell proliferation, induced cell cycle arrest and apoptosis and impaired oxidative phosphorylation in MM. CONCLUSIONS: PPRC1 acts as a prognostic biomarker and metabolic regulator in MM by sustaining mitochondrial oxidative phosphorylation. These findings highlight PPRC1 as a potential therapeutic target in MM.

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

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

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

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

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