SEARCH · Search PubMed

Results for “single cell sequencing”

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 307 records · Page 17Linked to original sources

PubMed2026

A transition zone enriched WIF1+ basal cell subtype is associated with benign prostatic hyperplasia.

The cellular composition and disease susceptibilities of the distinct zones of the human prostate remain incompletely understood. Benign prostatic hyperplasia (BPH) is a common condition that causes widespread morbidity and is nearly exclusively localized to the transition zone (TZ). Through extensive single-cell RNA sequencing (scRNA-seq) of benign regions from prostatectomy specimens, we identified a basal cell population expressing WIF1, VCAN, and NRG1, among other genes, that was significantly enriched in the TZ. Analysis of previously published scRNA-seq datasets further confirmed that WIF1+ basal cells were significantly enriched in BPH compared with normal prostate. Pathway and cell-cell communication analyses revealed that this basal subtype is associated with programs related to cell proliferation, epithelial-mesenchymal transition, immune regulation, angiogenesis, and hormone response. Together, the molecular signature, zonal distribution, and pathway enrichment suggest that TZ-enriched WIF1+ basal cells may contribute to BPH pathogenesis by promoting epithelial and stromal remodeling. © 2026 The Author(s). The Journal of Pathology published by John Wiley & Sons Ltd on behalf of The Pathological Society of Great Britain and Ireland.

NRG1↗

PubMed2026

Peripheral pain threshold, glycaemic status, and LAMP3 genetic variation: A community-based analysis.

Diabetic polyneuropathy is a common complication of diabetes, yet substantial inter-individual variation in peripheral pain perception suggests underlying genetic influences. This population-based study investigated clinical, metabolic, and genetic determinants of pain threshold using intraepidermal electrical stimulation in 906 participants from the Iwaki Health Promotion Project 2017. Genome-wide association analysis identified 12 loci showing suggestive associations, among which a missense variant in LAMP3 (rs482912) was prioritized as a biologically plausible candidate. Phenotype-stratified analyses showed that individuals carrying the CT or CC genotypes had lower PINT indices than those with the TT genotype, indicating reduced pain thresholds. Notably, the CC genotype retained an association with lower pain threshold using intraepidermal electrical stimulation under conditions of metabolic stress, including impaired glucose tolerance, elevated HbA1c, and obesity, whereas this association was attenuated in the presence of hypertension. Single-cell RNA sequencing analysis of human skin revealed that LAMP3-positive mature dendritic cells, enriched in immunoregulatory molecules, exhibited transcriptional enrichment of inflammatory, antigen-presenting, and nociception-related pathways, including NF-κB, JAK-STAT, cytokine signaling, and neuroimmune sensitization cascades. Autopsy-based skin analysis further demonstrated genotype-associated differences in dermal LAMP3-positive cell infiltration and CD8-positive T-cell abundance, while CD4-positive T-cell abundance and intraepidermal nerve fiber density remained unchanged across genotypes. Taken together, these findings suggest a potential association between LAMP3 variation and individual differences in peripheral pain threshold and provide biological context supporting a role for neuroimmune interactions in early sensory modulation under metabolic stress. Given the suggestive genetic evidence and indirect mechanistic data, these observations should be interpreted as exploratory and hypothesis-generating.

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

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↗

PubMed2026

Integration of single-cell transcriptomics and genomic mutation analysis identifies an immunotherapy-resistant tumor subcluster and validates ARNTL2 as a malignant driver in lung adenocarcinoma.

BACKGROUND: Immunotherapy resistance in lung adenocarcinoma (LUAD) remains a critical clinical challenge, and the mechanisms underlying resistance-associated intratumoral heterogeneity are poorly characterized. METHODS: We performed single-cell RNA sequencing of LUAD patients receiving neoadjuvant immunotherapy (responders vs. non-responders), integrating inferCNV, GSVA, and differential expression analyses. Cluster-specific genes were validated across seven independent cohorts (TCGA-LUAD, GSE13213, GSE26939, GSE29016, GSE30219, GSE31210, GSE42127). A multi-algorithm machine learning framework was used to construct a prognostic model, and the immune microenvironment was characterized using TCIA scoring, seven infiltration algorithms, and ESTIMATE. ARNTL2 function was assessed by CCK-8 and Transwell assays in A549 and H1299 cells. RESULTS: Non-responders showed significant enrichment of epithelial cells, depletion of cytotoxic T/NK cells, and elevated copy number variation burden versus responders (p < 0.0001). A resistance-enriched malignant subcluster (Cluster 2) exhibited hyperproliferative and metabolic reprogramming signatures with upregulated KRT17, S100A2, and CST6, which showed tumor-specific overexpression, adverse prognostic value, and genomic amplification across cohorts. CoxBoost combined with survivalSVM achieved optimal predictive performance (C-index = 0.686), yielding robust risk stratification (HR: 2.54-10.51, all p < 0.05). Low-risk patients showed greater immune infiltration and higher TCIA immunophenoscores. ARNTL2 was an independent prognostic factor (HR: 2.07-4.64) strongly correlated with risk score (r = 0.69), and its knockdown suppressed proliferation and invasion in both LUAD cell lines (all p < 0.05). CONCLUSION: This study identifies a resistance-associated malignant subcluster in LUAD, constructs a validated CoxBoost + survivalSVM prognostic model with robust immune stratification, and establishes ARNTL2 as a core oncogenic driver and therapeutic target.

ARNTL2↗

PubMed2026

Single-cell profiling reveals a novel CAF subpopulation linking stromal heterogeneity to immune suppression in breast cancer subtypes.

BACKGROUND: The tumor microenvironment critically influences breast cancer (BC) progression, immune surveillance, and therapeutic response. Cancer-associated fibroblasts (CAFs), a heterogeneous stromal population, are key regulators of these processes, yet their subtype-specific contributions in BC remain insufficiently defined. METHODS: We integrated three single-cell RNA sequencing datasets from 29 BC patients to characterize stromal populations. Bulk RNA-seq data from The Cancer Genome Atlas (TCGA) were analyzed to assess correlations between CAF subsets and immune infiltration. Gene signatures were derived to identify subtype-specific CAF-immune interactions, prognostic markers, and potential predictors of chemotherapy response. RESULTS: Three conserved stromal populations (iCAFs, myCAFs, and pericytes) were identified, along with a previously unrecognized subset, the cluster 3 (CL3) CAF-like cells, referred as metabolic stressed CAF (msCAF). msCAF cells displayed transcriptional programs associated with antigen presentation, stress response, glycolysis, and extracellular matrix remodeling. Their abundance was inversely correlated with T-cell infiltration and function, in a subtype-specific manner: triple negative breast cancer (TNBC) was enriched for msCAFs in immune-infiltrated but functionally constrained microenvironments, whereas Luminal A tumors exhibited weaker immune infiltration with heterogeneous CAF-immune associations. msCAFs were characterized by a conserved gene signature (HLA-A, HLA-C, IL32, EMP3) and subtype-specific genes related to T-cell exhaustion. Several genes demonstrated prognostic relevance with distinct patterns in Luminal A (IER3, TIMP1, TBX3, SEC61G) and TNBC (ADM, C4orf3, LDHA) tumors, as well as shared biomarkers (FN1, LOXL2, P4HA1). Multiple msCAF genes also predicted chemotherapy response, suggesting utility as treatment stratification biomarkers. CONCLUSION: msCAFs represent a clinically relevant CAF subset that drives immune suppression, impacts subtype-specific prognosis, and influences therapy response in BC. These findings highlight msCAFs as promising targets for enhancing immunotherapy and personalizing treatment strategies.

Humans↗

PubMed2025

Characterization of non-crossover recombination spectrum by single-microspore sequencing in maize and rice.

Meiotic DNA double-strand breaks (DSB) are crucial for chromosome recombination. The repair of DSB gives two outcomes: crossover (CO) and non-crossover (NCO). CO involves the bidirectional exchange between homologous chromosomes, whereas NCO refers to the unidirectional transfer of chromosome fragments. NCO can be categorized into NCO with gene conversion and NCO without gene conversion. Due to technological constraints, previous studies have focused more on CO than on NCO. In this study, we isolated single microspores from meiotic tetrads of maize (Zea mays) and rice (Oryza sativa) and conducted deep single-microspore genome sequencing to characterize NCO gene conversion (NCO-GC). Under highly stringent conditions, 101 CO and 902 NCO-GC tracts were identified in four maize tetrads, while 173 CO and 279 NCO-GC tracts were identified in six rice tetrads. In both maize and rice, NCO-GC was more prone to occur in the upstream and downstream of genes, as well as the introns. It also had a significant distribution in transposon regions. A common A-rich motif was enriched in the NCO-GC tracts of maize and rice. GC-biased gene conversion (gBGC) likely contributed to the bimodality of the GC content at the third codon position (GC3), and we discovered a significant proportional relationship between the number of DSBs and the GC content. These findings provide evidence that NCO-GC exhibits a distinct pattern compared with CO and may play an important role in gene and genome evolution.

Oryza↗

PubMed2026

Localized PD-1 CAR T therapy reprograms neuroinflammation.

B cell-depleting therapies are effective in multiple sclerosis (MS), yet some patients relapse, underscoring the need for more precise interventions. To identify new therapeutic targets, we generated a single-cell RNA sequencing (scRNA-seq) atlas of cerebrospinal fluid (CSF), brain, and blood from non-inflammatory controls and patients with MS or other neuroinflammatory diseases. We found disease-associated enrichment of class-switched immunoglobulin G+ (IgG+) B cells and plasma cells in MS CSF. Unbiased analysis identified a rare disease-enriched subset of activated, T cell receptor (TCR)-restricted, PD-1+ T follicular helper-like cells with B cell-recruiting features. To target this population, we developed PD-1-directed chimeric antigen receptor (CAR) T cells that selectively depleted pathogenic PD-1+ CD4 T cells and locally released IL-10. This strategy attenuated central nervous system (CNS) inflammation, reprogrammed the local immune milieu, and improved clinical outcomes across murine neuroinflammation models. These findings define a CNS-localized adaptive immune circuit in MS and nominate programmable PD-1 CAR T cells as a strategy to disrupt it.

Animals↗

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↗

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↗

PubMed2001

Regulation of immunoglobulin light chain gene rearrangements during early B cell development in the human.

Southern blot analyses of immunoglobulin light chain gene rearrangements in human leukemias and myelomas indicated that lambda loci in kappa-producing cells are largely unrearranged while kappa loci in lambda producers are often rearranged and inactivated by rearrangements of the kappa-deleting element (KDE). For a systematic analysis of the regulation of light chain rearrangements during early B cell development in normal human B cells also considering functionality of the rearrangements, we used FACS-sorted single naive kappa- and lambda-expressing B cells from peripheral blood of healthy humans. V(kappa)J(kappa) and V(lambda)J(lambda) joints and rearrangements involving the KDE were amplified simultaneously from single cells and sequenced. Whereas only 2 - 3 % of kappa-expressing cells carry V(lambda)J(lambda) joints, nearly all lambda-expressing cells have rearranged kappa loci and indeed carry V(kappa)J(kappa) joints. The V(kappa)J(kappa) joints in lambda-expressing cells exhibit preferential J(kappa)4 and J(kappa)5 over J(kappa)1 and J(kappa)2 usage compared to kappa-expressing cells. Thirty percent of the V(kappa)J(kappa) joints in lambda producers are rearranged in-frame. These data indicate extensive sequential V(kappa)-J(kappa) rearrangements and inactivation of functional V(kappa)J(kappa) joints in lambda-expressing cells, presumably before V(lambda)J(lambda) joining.

B-Lymphocytes↗

PubMed2026

Clinical, Histopathological, and Molecular Characterization of Pediatric MN1::ZNF341-Associated Cancer.

A lethal round-cell malignancy with an MN1::ZNF341 fusion has recently been reported in three infants. Here, we describe four further tumors, three in newborns (including monozygotic twins), and one in an adolescent. Detailed clinical, radiological, and histopathological data differentiate these tumors from their main mimics, neuroblastoma and round-cell sarcomas. Single-cell RNA sequencing confirms the tumor to be transcriptionally distinct from neuroblastoma, instead exhibiting steroidogenic differentiation. Whole genome and targeted DNA sequencing yield no further driver events. Our work reveals a broader clinicopathological phenotype than previously appreciated and corroborates suggestions that this tumor is a distinct and aggressive childhood cancer.

Humans↗

PubMed2026

Plasma Proteomics Identifies Thousand-and-One-Amino Acid Kinase 3 as a Potential Biomarker of Rheumatoid Arthritis Activity and a Novel Therapeutic Target.

OBJECTIVE: Bone destruction associated with active rheumatoid arthritis (RA) remains a major therapeutic challenge, with a lack of reliable molecular markers reflecting bone injury. This study aims to identify novel biomarkers linked to bone destruction in active RA through proteomic analysis, providing new strategies for precise monitoring and targeted therapy. METHODS: Data-independent acquisition mass spectrometry was used for proteomic quantification and bioinformatic analysis on plasma samples from 160 patients with RA and 40 healthy controls. Key proteins associated with bone destruction were screened by integrating Sharp scores with synovial single-cell RNA sequencing data and subsequently validated in two independent cohorts (N1 = 50 and N2 = 10) using enzyme-linked immunosorbent assay and multiplex immunohistochemistry. Functional studies were conducted using fibroblast-like synoviocytes (FLSs) in vitro and a collagen-induced arthritis (CIA) mouse model in vivo. RESULTS: A total of 4,998 plasma proteins were identified, with 506 showing significant differential expression between active and remitted RA. Thousand-and-one-amino acid kinase 3 (TAOK3) levels were positively associated with Sharp scores and markedly elevated in patients with active RA. Combining TAOK3 with C-reactive protein improved diagnostic accuracy for active RA (area under the curve = 0.915). High TAOK3 expression was also associated with increased relapse frequency. Functional studies showed that TAOK3 knockdown suppressed the tumor-like phenotype of FLSs and down-regulated matrix metalloproteinase 1/2/3 and cathepsin K, whereas TAOK3 overexpression promoted pannus cell-mediated bone erosion, mitigated by TAOK3-targeted inhibitor. In vivo, its inhibition showed therapeutic effects in CIA mice. CONCLUSION: TAOK3 serves as a potential biomarker for bone destruction in active RA and as a therapeutic target for precision monitoring and intervention.

Arthritis, Rheumatoid↗

PubMed2025

Balancing LncRNA H19 and miR-675 Bioconversion as a Key Regulator of Embryonic Myogenesis Under Maternal Obesity.

BACKGROUND: Maternal obesity (MO) impairs fetal skeletal muscle development, but the underlying mechanisms remain poorly defined. The regulatory roles of lncRNA H19 and its first exon derived microRNA675 (miR675) in prenatal muscle development remain to be examined. H19/Igf2 are in the same imprinting cluster with H19 expressed from the maternal allele while Igf2 expresses paternally. H19 contains a G-rich loop, and KH-type splicing regulatory protein (KHSRP) mediates the biogenesis of pre-miRNAs containing G-rich loops, which depends on its phosphorylation by AKT, a key mediator of IGF2 signalling. This study aims to depict the elusive function of these regulators that are affected by MO during embryonic myogenesis. METHODS: Single-cell transcriptomic sequencing and GeoMx spatial RNA sequencing were performed to identify the differentially expressed genes between embryos from MO and control (CT) mice. Both E11.5 and E13.5 embryos were collected and analysed to validate the sequencing data. The roles of H19 and miR657 in myogenesis were further analysed in P19 embryonic cells via CRISPR/dCas9-mediated H19 activation and inhibition. The epigenetic changes of H19 were analysed by methylated DNA immunoprecipitation, and allele-targeted analysis of H19 was performed by crossing C57BL/6J and CAST/EiJ mice. RESULTS: Transcriptomic analysis showed that MO embryos contained less differentiated myocytes (1.34%) than CT embryos (2.86%). Myogenesis-related GO biological processes were down-regulated in the MO embryonic myotome region. MO embryos showed lower expression of myogenic transcription factors such as Myf5, Myod1, Myog, Mef2c and Myh3 (p < 0.05). MO altered epigenetic modifications of the H19 genomic cluster, showing a decreased methylation level in H19 imprinting control region (p < 0.05) and a diallelic expression pattern of H19, which elevated its expression in MO embryos. Overexpression of H19 inhibited myogenesis in P19 cells, but miR675 promoted myogenesis, suggesting the critical regulatory roles of bioconversion of H19 to miR675. A KHSRP mediates the biogenesis of miR675, a process that relies on its phosphorylation by IGF2/AKT signalling. Knocking-down of KHSRP and inhibition of AKT abolished miR675 biogenesis. MO suppressed IGF2/AKT signalling and blocked KHSRP-dependent miR675 biogenesis in embryos. CONCLUSIONS: We found differential effects of H19 and miR675 on embryonic myogenesis. MO up-regulates H19 but blocks its miR675 bioconversion via suppressing IGF2/AKT/KHSRP signalling axis. Myogenesis in MO embryos was impeded due to the highly accumulated H19 and blocked miR675 biogenesis.

RNA, Long Noncoding↗

PubMed2026

NFS1 activates PI3K/AKT/mTOR signaling to upregulate GPX4 expression and enhance ferroptosis resistance in osteosarcoma.

Osteosarcoma continues to exhibit poor survival outcomes due to chemoresistance and metastasis, with metabolic reprogramming and ferroptosis resistance being key features of tumor heterogeneity, yet their upstream regulators remain poorly defined. NFS1, a cysteine desulfurase essential for iron-sulfur cluster biogenesis, protects multiple cancers from ferroptosis, but its role in osteosarcoma is unknown. In this study, we performed a transcriptomic meta-analysis and found that NFS1 expression was significantly upregulated in osteosarcoma tissues, with further elevation in metastatic lesions, and high NFS1 expression correlated with poor overall survival. Genome‑wide CRISPR screening data revealed a marked NFS1 dependency in osteosarcoma cell lines. Functionally, NFS1 promoted cell proliferation, migration, and invasion, whereas its knockdown suppressed these phenotypes. Using single‑cell RNA sequencing data from 27 osteosarcoma specimens, we applied a multi‑algorithm glycolytic scoring framework and observed NFS1 enrichment in highly glycolytic malignant cells, along with an association with PI3K/AKT/mTOR pathway activation. Mechanistically, NFS1 selectively enhanced PI3K, AKT, and mTOR phosphorylation without altering total protein levels, and upregulated GPX4, a central ferroptosis suppressor, leading to elevated ferroptosis resistance scores in NFS1‑high malignant cells. Collectively, these findings identify a previously unrecognized NFS1-PI3K/AKT/mTOR-GPX4 regulatory axis in osteosarcoma, linking metabolic reprogramming to ferroptosis resistance, and suggest that NFS1 functions as an oncogenic driver, as well as a promising prognostic biomarker and therapeutic target in osteosarcoma.

Humans↗

PubMed2026

Role of IFIT1 and IFIT3 in systemic lupus erythematosus: modeling a diagnosis and exploring immune regulation.

Systemic lupus erythematosus (SLE) is a complex autoimmune disorder characterized by multi-organ involvement and a protracted clinical course. Current diagnostic strategies, which rely heavily on clinical symptoms and serology, are often insufficient for early detection. Therefore, highly accurate diagnostic biomarkers are urgently needed to facilitate early intervention and optimize personalized treatment strategies. D atasets GSE61635 and GSE135779 were integrated to identify differentially expressed genes. Weighted gene co-expression network analysis (WGCNA) was performed to isolate the module with the strongest clinical relevance. Mendelian randomization and single‑cell RNA‑seq were used to identify key disease‑relevant genes. A diagnostic model was then constructed, and gene set variation analysis (GSVA), along with gene set enrichment analysis (GSEA), was conducted to elucidate the underlying molecular pathways. IFIT1 and IFIT3 were identified as 2 core genes highly expressed in monocytes and T cells of SLE patients. Functional enrichment analysis revealed that these genes were enriched in immune-related pathways, metabolic pathways related to inflammation and genomic stability. The diagnostic model showed good accuracy, with an area under the curve (AUC) of 0.974 on the training set and 0.912 on the validation set. IFIT1 and IFIT3 represent promising biomarkers for diagnosing SLE and appear to mediate key immune and metabolic disturbances. Furthermore, the developed model serves as an accurate and reliable instrument for early diagnosis and personalized therapy. Large-scale clinical studies are warranted to further validate these findings and evaluate their clinical application.

Humans↗

PubMed2026

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

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

Colorectal cancer↗

PubMed2026

Immune subtyping of colorectal adenoma identifies a subtype with activated adaptive immunity ahead of progressing to cancer.

BACKGROUND: Colorectal adenomas (CRA) represent precursor lesions with varying risks of malignant transformation. However, molecular subtyping, particularly immune-related classification, remains underexplored in adenomas. This study aims to characterize the immune landscape of CRA through immune subtyping and evaluate its association with cancer progression, gene expression signatures, and functional pathways. METHODS: We conducted a retrospective analysis of transcriptomic data from multiple cohorts of CRA samples. Immune subtypes were identified using non-negative matrix factorization (NMF) based on immune-related genes. Diverse deconvolution algorithms were used to estimate immune cell infiltration. The immune status alteration in premalignant lesion was further consolidated by single-cell transcriptome data. Differential gene expression analysis was performed between subtypes, followed by functional enrichment analyses (Gene Ontology [GO] and Kyoto Encyclopedia of Genes and Genomes [KEGG]). RESULTS: Two distinct immune subtypes were identified: an immune-enriched subtype characterized by high lymphocyte infiltration and elevated expression of immune-related genes, and an immune-deficient subtype with suppressed immune activity. Differential expression analysis revealed significant upregulation of immune response genes (e.g., CD4, CD86, HLA-DRA) in the immune-enriched subtype. GO and KEGG analyses highlighted enrichments in leukocyte transendothelial migration, chemokine signaling, and antigen processing and presentation pathways. Single-cell result revealed an early occurrence of TIGIT activation and exhausted CD8 T cell features in adenoma when compared to normal tissue. CONCLUSION: This study delineates distinct immune subtypes within CRAs. The immune-enriched subtype demonstrates activated adaptive immunity and may reflect a higher potential for immune surveillance, while the immune-deficient subtype exhibits stromal features suggestive of progressive transformation. These findings provide insights into early immune microenvironment alterations and may inform strategies for risk stratification and immunoprevention in colorectal carcinogenesis.

Colorectal adenoma↗