SEARCH · Search PubMed

Results for “Single-cell RNA sequencing (scRNA-seq)”

Search indexed PubMed citations on genomics, clinical trials, systematic reviews and public health. Explore titles, authors and supplied subject terms, then open the PubMed record.

SourceResearch connectionFrom yearThrough year

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

At least 55 records · Page 3Linked 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

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

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

Colorectal Neoplasms↗

PubMed2026

Increased IL4I1 expression predicts poor survival and modulates the immune microenvironment in acute myeloid leukemia.

BACKGROUND: The immunometabolic enzyme Interleukin-4-induced-1 (IL4I1) is implicated in cancer pathogenesis, yet its specific function and clinical relevance in acute myeloid leukemia (AML) remain unclear. METHODS: Comparative analysis of IL4I1 mRNA levels between AML patients and normal controls was performed using the Cancer Genome Atlas (TCGA) and Genotype-Tissue Expression (GTEx) databases. The Kaplan–Meier survival analysis was conducted to evaluate the prognostic value of IL4I1. Functional insights were derived from analyses of differentially expressed genes (DEGs), Gene Set Enrichment Analysis (GSEA), and Gene Ontology (GO)/Kyoto Encyclopedia of Genes and Genomes (KEGG) enrichment. Immune infiltration was evaluated using the ssGSEA, ESTIMATE, quanTIseq and single-cell RNA sequencing (scRNA-seq) analysis. Finally, in vitro and in vivo functional experiments were perfromed to explore the impact of IL4I1 on AML progression and immunoregulation. RESULTS: IL4I1 expression was significantly elevated in AML compared to normal controls (p = 0.0004) and associated with poorer overall survival (p = 0.003). Bioinformatic analysis revealed that IL4I1 was linked to immune-related pathways—including humoral immune response, leukocyte interactions, and chemokine signaling—and to cellular amino acid metabolism. Its expression correlated with immune cell infiltration and checkpoint molecule expression. Experimentally, IL4I1 promoted leukemia cell proliferation in vitro and in vivo (p < 0.05). Furthermore, silencing IL4I1 suppressed M2 macrophage polarization and reduced secretion of inflammatory factors (p < 0.05). CONCLUSIONS: IL4I1 may serve as a potential biomarker for poor prognosis and an attractive target for immune-based therapeutic interventions in AML.

Humans↗

PubMed2026

MX1+ effector T cells hyperactivation at the maternal-fetal interface in unexplained recurrent pregnancy loss.

BACKGROUND: Immune tolerance breakdown at the maternal-fetal interface is implicated in unexplained recurrent pregnancy loss (URPL), but the interplay between T cell hyperactivation and dendritic cells (DCs)-mediated signaling remains poorly defined. METHODS: First-trimester decidual tissues from 5 healthy controls and 6 URPL patients underwent single-cell RNA sequencing (scRNA-seq, 10× Genomics). Computational analyses included clustering (Seurat), trajectory inference (scTour), intercellular communication (CellChat) and metabolic pathway enrichment (Gene Ontology and scMetabolism). Flow cytometry was performed from 11 patients and 11 healthy controls. Spatial validation was performed via multiplex immunohistochemistry and immunohistochemistry on 12 additional controls and 12 URPL cases. Statistical significance was assessed using Student’s t-test. RESULTS: URPL decidua exhibited marked CD3+ T cells and MX1+effector T (Tem) cells infiltration and activation. Flow cytometry analysis confirmed a significant decidua-specific upregulation of T cell activation markers CD25 and CD69 specifically on the MX1+Tem subset in URPL patients compared to controls. MX1+Tem cell subset demonstrated interferon hyperactivation, proliferative hyperactivity and lipid-biased immunometabolism. Pseudotemporal analysis positioned MX1+ Tem cells between classical Tem and exhausted T cell states, suggesting progressive differentiation. CellChat identified DCs as key regulators of MX1+ Tem expansion via aberrant ICOSL signaling, validated by spatial co-localization of ICOSL+ DCs and MX1+ Tem cells in URPL tissues. CONCLUSION: Our findings demonstrate that the aberrant activation and proliferation of MX1+Tem cells as a key immunological feature associated with URPL patients.

Humans↗

PubMed2026

Integrated single-cell transcriptomics, Mendelian randomization, and machine learning identify CEBPZ as an immune-related biomarker in oral lichen planus.

BACKGROUND: Oral lichen planus (OLP) is a chronic, immune-mediated oral mucosal disease with complex pathophysiology and potential for malignant transformation. Understanding its molecular basis is critical for the development of precise diagnostic and therapeutic strategies. OBJECTIVES: We aimed to identify key immune-related biomarkers and characterize cellular dynamics in OLP, with a particular focus on the role of CEBPZ in disease pathogenesis. MATERIAL AND METHODS: We analyzed single-cell RNA sequencing (scRNA-seq) data from OLP lamina propria samples (GSE211630) to identify disease-specific T-cell subpopulations using high-dimensional weighted gene co-expression network analysis (hdWGCNA) for oxidative stress-related gene modules.-data-based Mendelian randomization (SMR) integrated FinnGen genome-wide association study (GWAS; 342,499 Europeans) data with Genotype-Tissue Expression (GTEx) expression quantitative trait loci (eQTL) data to identify causal genes. Machine learning (ML) models (least absolute shrinkage and selection operator (LASSO) and convolutional neural network (CNN)) were developed using bulk RNA-seq datasets (GSE52130 and GSE38616) for diagnostic purposes. RESULTS: We identified OLP-specific T-cell populations (clusters 0, 3, 5, 7, 13, and 15) with enhanced migration inhibition factor (MIF) pathway signaling toward B cells and monocytes. Two oxidative stress-associated modules contained hub genes, including CEBPZ. Summary-data-based Mendelian randomization analysis identified 231 OLP-associated genes, with CEBPZ uniquely intersecting LASSO-selected markers (odds ratio (OR) = 1.057, 95% confidence interval (95% CI) = 1.013-1.102, p = 0.010). Machine learning models achieved area under the curve (AUC) values ranging from 0.653 to 0.745, with the CNN model reaching a validation accuracy of 0.735. CEBPZ showed elevated expression in OLP T cells and correlated with enhanced MIF-(CD74+CXCR4) signaling. CONCLUSIONS: This integrative approach identifies CEBPZ as a pivotal biomarker linking genetic susceptibility, oxidative stress, and immune dysregulation in OLP. Our diagnostic models offer promising tools for OLP management.

CEBPZ↗

PubMed2026

Non-structural maintenance of chromosome condensin I complex subunit H knockdown suppresses malignant progression of esophageal squamous cell carcinoma via the Wnt/β-catenin signaling pathway.

BACKGROUND: Esophageal squamous cell carcinoma (ESCC) remains a major cause of cancer-related mortality, and effective therapeutic targets are still limited. Non-structural maintenance of chromosome condensin I complex subunit H (NCAPH) has been implicated in tumorigenesis; however, its clinical relevance, functional roles, and underlying mechanisms in ESCC are not fully defined. We aimed to characterize the expression pattern, prognostic value, biological functions, and mechanistic basis of NCAPH in ESCC. METHODS: Public datasets from The Cancer Genome Atlas (TCGA) and Gene Expression Omnibus (GEO) were analyzed to evaluate NCAPH expression and clinical associations. Single-cell RNA sequencing (scRNA-seq) data were used to map cell-type-specific distribution of NCAPH in tumor and adjacent tissues. NCAPH was silenced in KYSE150 and KYSE510 cells using lentiviral short hairpin RNAs (shRNAs), followed by Cell Counting Kit-8 (CCK-8), colony formation, wound-healing, and Transwell migration/invasion assays. A nude mouse xenograft model was established to assess the effect of NCAPH knockdown in vivo. RNA sequencing (RNA-seq), quantitative polymerase chain reaction (qPCR), western blotting, and enzyme-linked immunosorbent assay (ELISA) were performed to explore potential mechanisms. RESULTS: NCAPH was consistently upregulated in ESCC across multiple cohorts and was associated with unfavorable clinicopathological features and poorer survival. Functional assays demonstrated that NCAPH knockdown significantly inhibited ESCC cell proliferation, migration, invasion, and clonogenic growth. In vivo, NCAPH silencing suppressed xenograft tumor growth. Mechanistically, transcriptomic profiling and molecular validation indicated attenuation of Wnt/β-catenin signaling following NCAPH depletion, accompanied by reduced β-catenin and downstream targets. CONCLUSIONS: NCAPH promotes malignant progression of ESCC, at least in part through activation of the Wnt/β-catenin pathway, and may serve as a potential biomarker and therapeutic target.

Esophageal squamous cell carcinoma (ESCC)↗

PubMed2026

Construction of molecular signatures based on the co-expression network of NECSO-related gene TRPM4 and its prognostic value in hepatocellular carcinoma.

BACKGROUND: Hepatocellular carcinoma (HCC) demonstrates significant prognostic variability that is not entirely accounted for by traditional staging systems. Necrosis by sodium overload (NECSO) is an emerging programmed cell death pathway, but its clinical relevance in HCC remains undefined. Therefore, this study aimed to identify TRPM4-associated core genes, develop and validate a prognostic signature, and investigate its relationship with the tumor immune microenvironment, tumor mutational burden, and single-cell expression patterns in HCC. METHODS: We integrated transcriptomic, clinical, and mutational datasets from The Cancer Genome Atlas-Liver Hepatocellular Carcinoma (TCGA-LIHC) (n=421) and Gene Expression Omnibus (GEO) cohorts (n=115) to identify genes co-expressed with TRPM4-a key NECSO mediator-and those differentially expressed in HCC. A prognostic signature was developed using least absolute shrinkage and selection operator (LASSO)-Cox regression and validated through survival analysis, time-dependent receiver operating characteristic (ROC) curves, and multivariate Cox regression analysis. The immune landscape was characterized using CIBERSORT, somatic mutation data were used to calculate tumor mutational burden (TMB) and assess its correlation with the risk score, and single-cell RNA sequencing (scRNA-seq) resolved cell-type-specific expression patterns. RESULTS: From 294 TRPM4-associated core genes, we identified an 11-gene signature (BRSK1, MMP1, GRIN2D, GP6, MYOM2, N4BP3, CCDC112, TSEN54, MAP3K9, SPP1, B3GNT4) that independently predicted overall survival (OS) (hazard ratio =5.419, P<0.001) with areas under the curve (AUCs) of 0.779, 0.693, and 0.701 at 1, 3, and 5 years. These values were superior or comparable to conventional clinicopathologic variables after direct comparison. High-risk patients exhibited an immunosuppressive microenvironment, characterized by enrichment of M0 macrophage, a higher M2/M1 ratio (P<0.001) and distinct immune checkpoint profiles. When integrated with TMB, the prognostic stratification was further refined: high-TMB/high-risk patients had poorest outcomes (median OS, 15.3 months), while low-TMB/low-risk patients had the most favorable survival (median OS, 68.7 months). Single-cell analysis revealed that MMP1 was induced in cancer-associated fibroblasts (CAFs) and SPP1 was downregulated in macrophages, single-cell risk scores confirmed TAFs and macrophages as the main contributors to the prognostic model. CONCLUSIONS: The TRPM4-centered 11-gene signature provides robust and independent prognostic stratification in HCC by integrating immune, mutational, and single-cell features. This signature serves as a potential tool for prognostic evaluation and may help inform immunotherapeutic strategies for HCC.

Hepatocellular carcinoma (HCC)↗

PubMed2026

Identification of Drug-resistant Cell Subpopulations in Colorectal Cancer Through Single-cell Analysis and Exploration of Potential Therapeutic Strategies.

INTRODUCTION: The therapeutic efficacy of Colorectal Cancer (CRC) is often compromised by resistance to the standard chemotherapy agent oxaliplatin. METHODS: This study obtained single-cell RNA sequencing (scRNA-seq) data from the Gene Expression Omnibus (GEO) database. Differentially Expressed Genes (DEGs) between resistant and sensitive epithelial subpopulations were identified, followed by enrichment analysis. Pseudotemporal trajectory and cell-cell communication were analyzed using Monocle2 and CellChat, respectively. The candidate drug was predicted by Connectivity Map (cMAP) analysis. External validation included assessment of the EpC2 signature in an oxaliplatin-resistant cell line dataset (GSE76092), survival analysis using The Cancer Genome Atlas (TCGA) cohorts, and re-analysis of the GSE179784 dataset to assess the reproducibility of EpC2-like subpopulations and their DNA Damage Repair (DDR) scores. RESULTS: Cell subpopulations were divided into 10 clusters. Among them, epithelial cells comprised 5 subpopulations, with EPC2 identified as a potential oxaliplatin-resistant subset. DEGs were enriched in the TNF and IL-17 pathways. External validation confirmed the enrichment of EpC2 in resistant cell lines and its association with poor survival. Pseudotemporal trajectory revealed that epithelial cells underwent state transitions, forming two distinct branches. The resistant group exhibited enrichment in RNA splicing and NF-κB pathways. Cell-cell communication analysis revealed interactions involving MDK- NCL and PPIA-BSG. Dasatinib was predicted as a candidate drug. DISCUSSION: We identified an oxaliplatin-resistant subpopulation of Epithelial Cells (EpC2) in CRC, elucidated its multi-layered resistance mechanisms, and integrated multi- omics and cMAP database analyses to predict a potential intervention drug. CONCLUSION: This study provided potential therapeutic possibilities for oxaliplatin resistance, contributing to CRC treatment.

Humans↗

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

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

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

Female↗

PubMed2024

Genetic association between epilepsy and gliomas: Insights from Mendelian randomization and single-cell transcriptomic analyses.

BACKGROUND: Seizures are prevalent in glioma patients, especially in those with low-grade gliomas. The interaction between gliomas and epilepsy involves complex biological mechanisms that are not fully understood. METHODS: We collected Genome-Wide Association Study data for epilepsy and gliomas, performed differential expression analysis, and conducted Gene Ontology (GO) enrichment analysis on the identified genes. Single-cell RNA sequencing data (scRNA-seq) from GSE221534 dataset in Gene Expression Omnibus (GEO) were used to analyze cell-cell interactions within glioma samples from patients with and without epilepsy. RESULTS: Mendelian Randomization (MR) analysis revealed significant associations between genetic variants related to epilepsy and glioma risk, suggesting a potential causal relationship, especially in astrocytomas. Differential expression analysis identified epilepsy-related genes that were significantly upregulated in astrocytoma tissues compared to normal brain tissues. GO enrichment analysis indicated that these genes are involved in critical biological processes such as neurogenesis and cellular signaling. The scRNA-seq analysis showed, compared to non-epileptic samples, glioma stem cells, microglia, and NK cells are increased in the core regions of astrocytomas in epileptic patients. Additionally, intercellular communication between tumor cells and other non-tumor cells is markedly enhanced in astrocytoma samples from epileptic patients. CONCLUSION: This study provides evidence of a genetic association between epilepsy and gliomas and elucidates the biological mechanisms through which epilepsy may influence glioma progression.

Humans↗

PubMed2026

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

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

Journal Article↗

PubMed2024

Single cell RNA sequencing provides novel cellular transcriptional profiles and underlying pathogenesis of presbycusis.

Age-related hearing loss (ARHL) or presbycusis is associated with irreversible progressive damage in the inner ear, where the sound is transduced into electrical signal; but the detailed mechanism remains unclear. Here, we sought to determine the potential molecular mechanism involved in the pathogeneses of ARHL with bioinformatics methods. A single-cell transcriptome sequencing study was performed on the cochlear samples from young and aged mice. Detection of identified cell type marker allowed us to screen 18 transcriptional clusters, including myeloid cells, epithelial cells, B cells, endothelial cells, fibroblasts, T cells, inner pillar cells, neurons, inner phalangeal cells, and red blood cells. Cell-cell communications were analyzed between young and aged cochlear tissue samples by using the latest integration algorithms Cellchat. A total of 56 differentially expressed genes were screened between the two groups. Functional enrichment analysis showed these genes were mainly involved in immune, oxidative stress, apoptosis, and metabolic processes. The expression levels of crucial genes in cochlear tissues were further verified by immunohistochemistry. Overall, this study provides new theoretical support for the development of clinical therapeutic drugs.

Animals↗

PubMed2027

Teratoma Formation and Genomic Profiling Using Multi-Omics Approaches.

Teratoma formation is the gold standard assay for evaluating the developmental pluripotency of human and mouse embryonic stem cells (ESCs) and induced pluripotent stem cells (iPSCs). Following subcutaneous injection into immunodeficient mice, pluripotent stem cells spontaneously differentiate into derivatives representing all three embryonic germ layers-ectoderm, mesoderm, and endoderm. Beyond serving as a functional assay for pluripotency, teratomas provide a unique three-dimensional model system for studying early human development and lineage specification in vivo. This chapter describes comprehensive protocols for teratoma formation in immunodeficient mice, tissue processing for multiple downstream genomic applications, and multi-omics profiling approaches. We detail methods for embryonic stem cell culture, teratoma generation via subcutaneous injection, tissue dissection and processing for chromatin immunoprecipitation followed by sequencing (ChIP-Seq), RNA sequencing (RNA-Seq), single-cell multiome profiling combining chromatin accessibility (ATAC-Seq) and gene expression (scRNA-Seq), and histological analysis using hematoxylin and eosin (H&E) staining. Additionally, we provide bioinformatics workflows for analyzing the resulting genomic datasets to characterize the epigenetic and transcriptional landscapes of teratoma-derived tissues. These methods enable comprehensive molecular characterization of developmental processes and provide valuable resources for stem cell biologists studying pluripotency, differentiation, and early embryonic development.

Teratoma↗

PubMed2026

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

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

Gene expression profiling↗

PubMed2025

Progressive T cell exhaustion and predominance of aging tissue associated macrophages with advancing disease stage in penile squamous cell carcinoma.

Penile squamous cell carcinoma (PSCC) is a rare malignancy with limited understanding of the tumor immune microenvironment (TIME). The interplay between PSCC and the immune system across disease progression and HPV infection status remains poorly characterized. This study aims to assess the TIME changes from localized to advanced disease and between HPV-positive versus negative tumors to identify potential immune evasion mechanisms in advanced PSCC. scRNA-seq was performed on ten PSCC tissue samples from penile, lymph node and distant metastatic sites with four matched penile and lymph node samples to understand the cellular heterogeneity within PSCC tumors. Analysis of immune cell populations and transcriptional hallmarks were performed stratified by localized (pT1-3, N0) versus advanced (N1-3, M0 or any N, M1) disease states and HPV infection status. We observed significant differences in immune cell infiltration between localized and advanced PSCC disease states and by HPV status. Advanced disease states demonstrated an exhausted immune phenotype, characterized by terminally exhausted CD8+ T cells, M2-like macrophages and hypoxic signature, while localized disease states demonstrated an active innate immune system characterized by increased DCs. HPV-negative tumors displayed low immune cell infiltration while HPV-positive tumors demonstrated an immune exhausted phenotype. These findings offer valuable insights into the evolving PSCC immune landscape, paving the way for the development of potential therapeutic approaches for advanced PSCC.

Humans↗

PubMed2026

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

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

Software↗

PubMed2026

Nested co-expression network analysis identifies compact gene clusters in a black box.

MOTIVATION: Digital analysis of biological systems requires methods capable of identifying both broad and nested gene modules reflecting complex biological processes. Existing transcriptomic methods often miss compact gene sets corresponding to subprocesses in specialized cell types, limiting insights into functional heterogeneity. RESULTS: We present Nested-WGCNA, a two-stage unsupervised network analysis algorithm designed to identify coarse-grained and fine-grained gene modules. Applied to bulk RNA-Seq data, Nested-WGCNA reveals stable modules reproducible across datasets. When validated against scRNA-Seq data, these modules correspond to both major and minor immune cell subtypes. Application to immunotherapy response datasets uncovers predictive and prognostic biomarkers, highlighting its utility in treatment stratification and biomarker discovery. AVAILABILITY: The NestedWGCNA source code and analysis pipeline are available on GitHub (https://github.com/ilyada/NestedWGCNA) and archived on Zenodo (https://doi.org/10.5281/zenodo.18959244).

Algorithms↗