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 379 records · Page 21Linked to original sources

PubMed2025

BCLAF1 restrains stress responses in hematopoietic stem cells to support expansion and repopulation.

Hematopoietic stem cells (HSCs) rapidly expand during fetal development and after stress. Here, we identify B-cell lymphoma-2-associated factor 1 (BCLAF1) as a regulator of HSC repopulation activity, with roles in the expansion of fetal HSCs and hematopoietic reconstitution after stem cell transplantation. Using mice with hematopoietic-specific and inducible deletion of Bclaf1, we find that BCLAF1 promotes fetal HSC development but is dispensable for the maintenance of adult HSCs at steady state. Loss of BCLAF1 in either fetal or adult HSCs significantly impairs their self-renewal and multilineage reconstitution activity after stem cell transplantation. Single-cell RNA sequencing of fetal hematopoietic progenitors reveals that loss of BCLAF1 reduces long-term HSCs and restrains the expression of stress response genes. BCLAF1 associates with chromatin throughout the genome of fetal and adult hematopoietic cells, likely through indirect mechanisms, to regulate transcriptional programs. These results establish a novel function for the transcriptional regulator BCLAF1 in limiting stress responses in HSCs, thereby preserving HSC development during embryogenesis and repopulation function after stem cell transplant.

Hematopoietic Stem Cells↗

PubMed2025

The MYC/TXNIP axis mediates NCL-Suppressed CD8+T cell immune response in lung adenocarcinoma.

BACKGROUND: Lung adenocarcinoma is a deadly malignancy with immune evasion playing a key role in tumor progression. Glucose metabolism is crucial for T cell function, and the nucleolar protein NCL may influence T cell glucose metabolism. This study aims to investigate NCL's role in T cell glucose metabolism and immune evasion by lung adenocarcinoma cells. METHODS: Utilizing single-cell RNA sequencing (scRNA-seq) data from the Gene Expression Omnibus (GEO) and The Cancer Genome Atlas (TCGA), we analyzed cell clustering, annotation, and prognosis. In vitro experiments involved manipulating NCL expression in CD8+ T cells to study immune function and glucose metabolism. In vivo studies using an orthotopic transplant mouse model monitored NCL's impact on CD8+ T cell glucose metabolism and anti-tumor immune function. RESULTS: NCL was associated with T cell dysfunction and glucose metabolism. NCL silencing enhanced CD8+ T cell glucose metabolism, cytotoxicity, and infiltration, while NCL overexpression had the opposite effect. NCL overexpression relieved MYC-mediated transcriptional repression of TXNIP, reducing CD8+ T cell glucose metabolism. In vivo, NCL inhibited CD8+ T cell glucose metabolism through the MYC/TXNIP axis, hindering anti-tumor immune function. CONCLUSIONS: NCL overexpression suppresses CD8+ T cell glucose metabolism and anti-tumor immune function, promoting lung adenocarcinoma progression via the MYC/TXNIP axis.

CD8-Positive T-Lymphocytes↗

PubMed2025

Transcriptome-wide analysis reveals potential roles of CFD and ANGPTL4 in fibroblasts regulating B cell lineage for extracellular matrix-driven clustering and novel avenues for immunotherapy in breast cancer.

BACKGROUND: The remodeling of the extracellular matrix (ECM) plays a pivotal role in tumor progression and drug resistance. However, the compositional patterns of ECM in breast cancer and their underlying biological functions remain elusive. METHODS: Transcriptome and genome data of breast cancer patients from TCGA database was downloaded. Patients were classified into different clusters by using non-negative matrix factorization (NMF) based on signatures of ECM components and regulators. Weighted Gene Co-expression Network Analysis (WGCNA) was used to identify core genes related to ECM clusters. Additional 10 independent public cohorts including Metabric, SCAN_B, GSE12276, GSE16446, GSE19615, GSE20685, GSE21653, GSE58644, GSE58812, and GSE88770 were collected to construct Training or Testing cohort, following machine learning calculating ECM correlated index (ECI) for survival analysis. Pathway enrichment and correlation analysis were used to explore the relationship among ECM clusters, ECI and TME. Single-cell transcriptome data from GSE161529 was processed for uncovering the differences among ECM clusters. RESULTS: Using NMF, we identified three ECM clusters in the TCGA database: C1 (Neuron), C2 (ECM), and C3 (Immune). Subsequently, WGCNA was employed to pinpoint cluster-specific genes and develop a prognostic model. This model demonstrated robust predictive power for breast cancer patient survival in both the Training cohort (n = 5,392, AUC = 0.861) and the Testing cohort (n = 1,344, AUC = 0.711). Upon analyzing the tumor microenvironment (TME), we discovered that fibroblasts and B cell lineage were the core cell types associated with the ECM cluster phenotypes. Single-cell RNA sequencing data further revealed that angiopoietin like 4 (ANGPTL4)+ fibroblasts were specifically linked to the C2 phenotype, while complement factor D (CFD)+ fibroblasts characterized the other ECM clusters. CellChat analysis indicated that ANGPTL4+ and CFD+ fibroblasts regulate B cell lineage via distinct signaling pathways. Additionally, analysis using the Kaplan-Meier Plotter website showed that CFD was favorable for immunotherapy response, whereas ANGPTL4 negatively impacted the outcomes of cancer patients receiving immunotherapy. CONCLUSION: We identified distinct ECM clusters in breast cancer patients, irrespective of molecular subtypes. Additionally, we constructed an effective prognostic model based on these ECM clusters and recognized ANGPTL4+ and CFD+ fibroblasts as potential biomarkers for immunotherapy in breast cancer.

Humans↗

PubMed2026

Epigenetic and immunological alterations in umbilical cord blood of overweight/obese women with gestational diabetes mellitus: insights into DNA methylation signatures and immune cell dysregulation.

BACKGROUND: Gestational diabetes mellitus (GDM) is a common pregnancy complication associated with adverse maternal and neonatal outcomes. Epigenetic modifications may reflect intrauterine metabolic exposure and contribute to immune and metabolic alterations. This study aimed to explore DNA methylation profiles in umbilical cord blood from overweight and obese women with and without GDM. METHODS: Umbilical cord blood samples from 30 overweight/obese pregnant women (with and without GDM) were analyzed using the Illumina 850 K methylation array to identify differentially methylated positions (DMPs) and regions (DMRs). Gene Ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) enrichment analyses were performed to assess the functional relevance of methylation changes. Immune cell composition was estimated using deconvolution analysis and further examined in an independent single-cell RNA sequencing (scRNA-seq) cohort. Lasso regression was applied to identify CpG sites associated with GDM status and construct a preliminary methylation-based classification model. RESULTS: A total of 23,331 hypermethylated and 29,501 hypomethylated DMPs were identified between women with and without GDM, with hypomethylation predominating. Enrichment analyses indicated associations with neurodevelopmental pathways, metabolic processes, immune regulation, and epigenetic modification. Immune deconvolution analysis suggested reduced proportions of CD4+ T cells (p < 0.05) and a trend toward decreased NK cells in the GDM group, alongside increased CD8+ T cells and neutrophils. Seven CpG sites were selected for model construction and demonstrated strong discriminatory performance within this cohort. CONCLUSION: This exploratory study identifies distinct cord blood DNA methylation patterns associated with GDM in overweight/obese pregnancies. The findings suggest potential links between epigenetic alterations and immune cell composition in GDM-exposed offspring. The identified CpG signature warrants further validation in larger, prospective cohorts to determine its clinical applicability.

Humans↗

PubMed2026

Thyroid-stimulating hormone receptor mediates peripheral-central neuroimmune crosstalk in autoimmune thyroid diseases.

BACKGROUND: Organ-specific autoimmune diseases, particularly Graves' disease (GD) and its extrathyroidal manifestation, Graves' orbitopathy (GO), are characterized by systemic autoimmunity that may extend its impact to the central nervous system (CNS). While thyroid-stimulating hormone receptor (TSHR) is the primary driver of pathological remodeling in the thyroid and orbital tissues, emerging evidence suggests it is also expressed in the brain and may participate in neuroimmune signaling. However, the molecular mechanisms linking peripheral TSHR-driven autoimmunity to these extended systemic features remain unclear. Thus, GD and GO provide a unique window to investigate how peripheral autoantibodies influence CNS involvement as part of its broader pathological spectrum. METHODS: Genome-wide association studies (GWAS) and post-GWAS analyses were integrated with bulk RNA sequencing, single-cell and spatial transcriptomics, and brain imaging phenotypes to comprehensively characterize peripheral and central alterations in GD and GO. Mendelian randomization was applied to test causal relationships between genetic variants and brain signatures. Structural biology analyses were further conducted including protein-protein docking, small-molecule docking, and normal mode dynamics to identify prospective modulators of TSHR. Immunofluorescence staining was performed in a GO mouse model to validate the colocalization of potential interacted proteins in the specific brain region. RESULTS: Brain imaging-derived phenotypes (IDPs) alterations in GO and GO were systematically analyzed to identify neuroanatomical and functional alterations. TSHR was further identified as a shared genetic driver across peripheral and central compartments. TSHR was expressed in spiny projection neurons, microglia, and peripheral T cells, with cell-cell communication analyses highlighting TSHR-mediated interactions among neurons, endothelial cells, and microglia. Immunofluorescence staining in a GO mouse model confirmed the colocalization of TSHR with FN1 and GNAS in the basal ganglia, providing tissue-level validation of the computationally predicted ligand-receptor interactions. Immune profiling further showed immune alterations in GD and GO. Structural modeling supported plausible physical interfaces between TSHR and interacting proteins, and small-molecule screening identified three repurposable compounds - venetoclax, irinotecan, and dutasteride - with predicted favorable docking scores and stable binding poses in our simulations. CONCLUSIONS: These findings demonstrate that TSHR acts as a molecular hub mediating peripheral-central neuroimmune crosstalk in GD and GO. The results support a broader "disease-molecule axis" framework that links genetic susceptibility with multi-level immune and neural mechanisms. This work provides mechanistic insights relevant to the development of TSHR-targeted therapies, with implications for both peripheral immune modulation and central regulation. However, the limited sample size, lack of longitudinal follow-up, and absence of in vivo validation warrant cautious interpretation and further investigation.

Receptors, Thyrotropin↗

PubMed2024

Integrated bioinformatics analysis and experimental validation reveal the relationship between ALOX5AP and the prognosis and immune microenvironment in glioma.

BACKGROUND: Treatment of gliomas, the most prevalent primary malignant neoplasm of the central nervous system, is challenging. Arachidonate 5-lipoxygenase activating protein (ALOX5AP) is crucial for converting arachidonic acid into leukotrienes and is associated with poor prognosis in multiple cancers. Nevertheless, its relationship with the prognosis and the immune microenvironment of gliomas remains incompletely understood. METHODS: The differential expression of ALOX5AP was evaluated based on public Databases. Kaplan-Meier, multivariate Cox proportional hazards regression analysis, time-dependent receiver operating characteristic, and nomogram were used to estimate the prognostic value of ALOX5AP. The relationship between ALOX5AP and immune infiltration was calculated using ESTIMATE and CIBERSORT algorithms. Relationships between ALOX5AP and human leukocyte antigen molecules, immune checkpoints, tumor mutation burden, TIDE score, and immunophenoscore were calculated to evaluate glioma immunotherapy response. Single gene GSEA and co-expression network-based GO and KEGG enrichment analysis were performed to explore the potential function of ALOX5AP. ALOX5AP expression was verified using multiplex immunofluorescence staining and its prognostic effects were confirmed using a glioma tissue microarray. RESULT: ALOX5AP was highly expressed in gliomas, and the expression level was related to World Health Organization (WHO) grade, age, sex, IDH mutation status, 1p19q co-deletion status, MGMTp methylation status, and poor prognosis. Single-cell RNA sequencing showed that ALOX5AP was expressed in macrophages, monocytes, and T cells but not in tumor cells. ALOX5AP expression positively correlated with M2 macrophage infiltration and poor immunotherapy response. Immunofluorescence staining demonstrated that ALOX5AP was upregulated in WHO higher-grade gliomas, localizing to M2 macrophages. Glioma tissue microarray confirmed the adverse effect of ALOX5AP in the prognosis of glioma. CONCLUSION: ALOX5AP is highly expressed in M2 macrophages and may act as a potential biomarker for predicting prognosis and immunotherapy response in patients with glioma.

Humans↗

PubMed2025

Genetic targets related to aging for the treatment of coronary artery disease.

BACKGROUND: Coronary Artery Disease (CAD) is the most common cardiovascular disease worldwide, threatening human health, quality of life and longevity. Aging is a dominant risk factor for CAD. This study aims to investigate the potential mechanisms of aging-related genes and CAD, and to make molecular drug predictions that will contribute to the diagnosis and treatment. METHODS: We downloaded the gene expression profile of circulating leukocytes in CAD patients (GSE12288) from Gene Expression Omnibus database, obtained differentially expressed aging genes through "limma" package and GenaCards database, and tested their biological functions. Further screening of aging related characteristic genes (ARCGs) using least absolute shrinkage and selection operator and random forest, generating nomogram charts and ROC curves for evaluating diagnostic efficacy. Immune cells were estimated by ssGSEA, and then combine ARCGs with immune cells and clinical indicators based on Pearson correlation analysis. Unsupervised cluster analysis was used to construct molecular clusters based on ARCGs and to assess functional characteristics between clusters. The DSigDB database was employed to explore the potential targeted drugs of ARCGs, and the molecular docking was carried out through Autodock Vina. Finally, single-cell data (GSE159677) of arterial intima was used to further explore the expression of aging signature genes in different cell subpopulations. RESULTS: We identified 8 ARCGs associated with CAD, in which HIF1A and FGFR3 were up while NOX4, TCF7L2, HK3, CDK18, TFAP4, and ITPK1 were down in CAD patients. Based on this, CAD patients can be divided into two molecular clusters, among which cluster A mainly involves functional pathways such as ECM receptor interaction and focal adhesion; cluster B mainly involves functional pathways such as amimo sugar and nucleotide sugar metabolism and pyrimidine metabolism. In addition, the molecular docking results showed that retinoic acid and resveratrol had good binding affinity with targets genes. Further single-cell analysis results showed that NOX4, TCF7L2, ITPK1, and HIF1A were specifically expressed in different types of cells in atherosclerotic tissues. CONCLUSION: Our study identified several ARCGs that may be involved in the pathogenesis and progression of CAD. Further, retinoic acid and resveratrol were potential candidate molecule drugs for inhibiting these targets.

Humans↗

PubMed2025

IGF1R deficiency mitigates acute lung injury by promoting anti-inflammatory transcriptional profiles.

BACKGROUND: Acute lung injury (ALI), acute respiratory distress syndrome (ARDS) and COVID-19 are characterized by hyperinflammation, commonly referred to as "cytokine storm". The insulin-like growth factor (IGF) pathway, particularly the type 1 receptor (IGF1R), plays a critical role in lung homeostasis and has been implicated in the pathogenesis of pulmonary inflammatory diseases. In mice, widespread Igf1r deficiency attenuates lung inflammation and alveolar damage in bleomycin (BLM)-induced ALI. METHODS: We analyzed single-cell RNA sequencing datasets from lung tissue of COVID-19 cases and control donors as well as mouse lungs to determine Igf1r and IGF family expression across pulmonary cell types. Furthermore, we conducted bulk RNA sequencing on lungs from Igf1r-deficient mice three days after BLM or saline instillation, followed by differential expression and functional enrichment analyses. Findings were further tested through protein detection, assessment of DNA damage and methylation in lung tissues, and functional assays using Igf1r-deficient primary mouse embryonic fibroblasts (MEFs). RESULTS: IGF1R was broadly expressed across multiple cell types in both human and mouse lungs under normal and pathological conditions. Other IGF family members showed cell-type-specific expression, which was modulated by lung injury. Transcriptomic profiling revealed differentially expressed genes between BLM-challenged and control mouse lungs, detecting biological processes and signaling pathways involved in ALI pathobiology. Igf1r deficiency in BLM-challenged mice reversed a large fraction of the transcriptional changes triggered by BLM, including "cytokine storm"-related gene expression. Functional enrichment analysis additionally revealed significant modulation of pathways related to DNA damage, metabolic reprogramming, mitochondrial homeostasis, and epigenetic regulation. In vitro, Igf1r-deficient MEFs exhibited decreased mitochondrial respiration and glycolysis, protection against BLM-induced nuclear damage and mitochondrial accumulation, and decreased histone H3 acetylation. Moreover, Igf1r-deficient mouse lungs displayed increased global DNA methylation following BLM challenge. CONCLUSIONS: IGF1R is a key modulator of the inflammatory and molecular response to ALI pathogenesis. IGF1R deficiency dampens the "cytokine storm", modifies transcriptional and epigenetic profiles and promotes protective cellular responses. These findings highlight IGF1R signaling as a potential therapeutic target in ARDS and related lung injuries.

Animals↗

PubMed2026

Necroptosis in alveolar epithelium orchestrates lung ischemia-reperfusion injury: a multi-omics study.

BACKGROUND: Lung ischemia-reperfusion injury (LIRI) is a leading cause of early morbidity and mortality following lung transplantation and other cardiopulmonary procedures. It is characterized by acute sterile inflammation driven by regulated cell death (RCD). While various RCD modalities, including apoptosis, necroptosis, pyroptosis, and ferroptosis, have been implicated in lung injury, their relative contributions and distinct activation patterns in LIRI remain poorly defined. METHODS: We employed an integrated multi-omics approach combining transcriptomics and proteomics with histological and functional validations in a murine hilar clamping model of LIRI. Key findings were further corroborated using single-cell RNA sequencing (scRNA-seq) data from human lung transplant recipients. The functional role of necroptosis was validated using pharmacological inhibitors (Nec-1, GSK'872) and Mlkl-deficient (Mlkl-/-) mice. RESULTS: LIRI triggered acute, time-dependent lung injury peaking within 24 h of reperfusion. Although transcriptomic profiling suggested broad activation of multiple RCD pathways, proteomic and biochemical analyses revealed a distinct landscape in our experimental setting: markers of apoptosis, pyroptosis, and ferroptosis were either downregulated or showed no significant positive correlation with injury severity and inflammatory peaks. In contrast, the necroptotic pathway emerged as a highly activated modality. Specifically, necroptosis, marked by phosphorylated RIPK1, RIPK3, and MLKL, was localized primarily in alveolar epithelial cells, correlated strongly with cytokine release and histological lung injury, and preceded the inflammatory response. Pharmacological inhibition or genetic ablation of necroptosis significantly attenuated tissue damage and inflammation. This pronounced necroptotic signature appeared distinct from the broad multi-pathway activation observed in lipopolysaccharide (LPS)-induced lung injury. Translational analysis of human scRNA-seq data further confirmed the selective upregulation of necroptosis signatures in alveolar type 2 (AT2) cells following lung transplantation. CONCLUSION: Our multi-omics analysis identifies necroptosis, particularly in alveolar epithelial cells, as a critical driver of sterile inflammation and tissue injury in the early phase of LIRI. Targeting alveolar epithelial necroptosis may represent a precise and promising therapeutic strategy for lung transplantation and ischemia-reperfusion-associated pulmonary disorders.

Animals↗

PubMed2026

Ex vivo long-term expansion of human hematopoietic stem and progenitor cells as a tool for modeling vector integration sites and clonality.

BACKGROUND: Gene therapy (GT) using retroviral vectors (RVs) is efficacious in treating monogenic diseases. However, there is an inherent risk for severe adverse effects due to insertional mutagenesis. Preclinical safety assessment and patient monitoring are inevitable in GT. To assess the genotoxic risk of novel RV vectors, mainly murine hematopoietic stem and progenitor cells (HPSCs) are routinely used, because human HSPCs cannot be immortalized in vitro using mutagenic vectors. In this study, we aim to identify early signs of clonal outgrowth by performing integration site analyses (ISA). METHODS: The small molecules A83-01, pomalidomide, and UM171 (APU) were used for the ex vivo expansion, lentiviral transduction, and long-term cultivation of umbilical cord blood-derived HSPCs. We determined the influence of APU on the stemness of HSPCs and their differentiation capacity via single-cell RNA sequencing (scRNA seq) and in xenotransplantation studies. To track vector insertion site dynamics, we transduced 7-day expanded HSPCs with a mutagenic or a safer RV. ISA was conducted in human HSPCs over a 5-week cultivation in vitro and compared to the bone marrow of xenotransplanted mice to assess clonal skewings. RESULTS: APU supported the expansion of CD34+CD38-CD45RA-CD90+EPCR+ HSPCs. scRNA seq confirmed the enrichment of HSC signature genes in APU-expanded HSPCs compared to the clinically used medium SFT3 (SCF, FLT3-L, TPO, IL-3). After RV transduction, APU still maintained around 30% of CD34+ cells for 5 more weeks. Without the compounds, already 2 weeks post-transduction, less than 10% of cells were CD34+. The long-term culture allowed the detection of high-risk integrations of the mutagenic SIN-LV.SF in MEIS1 or SUSD6 due to their increasing abundance over time. Bone marrow of xenotransplanted mice was less clonal but did not support the outgrowth of insertional mutants. Overall, APU increased clonal diversity. CONCLUSIONS: Our findings propose that long-term cultivation of transduced HSPC in APU allows for outgrowth of clonal integration sites. The decrease of clonality has been observed in gene therapy patient's years after treatment. Thus, the in vitro model could be used to develop novel human HSPC-based genotoxicity assays that predict insertional mutagenesis, in addition to existing preclinical biosafety assays.

Humans↗

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

Analysis and validation of abnormal signaling pathways and immune cell infiltration characteristics in digestive system cancers based on peroxisome-related genes.

BACKGROUND: Although emerging evidence suggests a role for peroxisomes in tumorigenesis, their functions in digestive cancers remain unclear. This study aims to investigate the association between peroxisomes and digestive tract tumors. METHODS: To systematically investigate peroxisomal functions in digestive cancers, we first constructed and validated tumor-specific prognostic signatures based on peroxisome-related genes (PRGs) through univariate Cox, least absolute shrinkage and selection operator (LASSO), and multivariate Cox regression analyses. We then characterized the tumor immune microenvironment (TIME) with CIBERSORT, X-CELL, and EPIC algorithms, and identified tumor-specific and common signalings via Gene Ontology (GO), Kyoto Encyclopedia of Genes and Genomes (KEGG), and gene set enrichment analysis (GSEA). Focusing on hepatocellular carcinoma (HCC), we experimentally validated peroxisome-related therapeutic responses by profiling signature genes in radioresistant cells and an orthotopic transarterial chemoembolization (TACE) rat model. PEX13 knockdown further assessed peroxisomal role in radiosensitivity and targeted therapy response. Clinical relevance of PEX13 was evaluated in HCC cohort. Single-cell RNA sequencing dataset and lipidomics further revealed peroxisomal mechanisms in HCC progression. Finally, peroxisomal function in colorectal cancer (CRC) was validated in vitro. RESULTS: Novel peroxisome-related prognostic signatures demonstrated strong predictive power in HCC, colon adenocarcinoma, rectal adenocarcinoma, pancreatic adenocarcinoma, gastric adenocarcinoma, esophageal adenocarcinoma, esophageal squamous cell carcinoma, and cholangiocarcinoma. High-risk patients displayed an immunosuppressive microenvironment, characterized by increased infiltration of regulatory T cells, M2 macrophages, Th2 cells, or cancer-associated fibroblasts, or Th1 cells' reduction. Peroxisomes engaged in several distinct yet convergent pathways, most notably "positive regulation of response to stimuli". HCC prognostic genes were dynamically regulated in response to therapeutic stimuli, including radiotherapy, targeted therapy, and TACE. Clinically, the expression of PEX13 was markedly upregulated in tumor tissues from therapy-resistant HCC patients. Mechanistically, peroxisomal dysfunction induced by silencing PEX13 in HCC or UBE2D2 in CRC may overcome therapeutic resistance (radiotherapy/ lenvatinib resistance in HCC, radioresistance in CRC) through reprogramming lipid metabolism. CONCLUSIONS: Peroxisomes act as pivotal regulators of digestive cancer progression by modulating signaling pathways, the TIME, therapeutic resistance, and lipid metabolism. Targeting peroxisomal function, particularly in high-risk subgroups of HCC and CRC, warrants further exploration as a promising therapeutic strategy.

Peroxisomes↗

PubMed2025

Understanding disease-associated metabolic changes in human colonic epithelial cells using the iColonEpithelium metabolic reconstruction.

The colonic epithelium plays a key role in the host-microbiome interactions, allowing uptake of various nutrients and driving important metabolic processes. To unravel detailed metabolic activities in the human colonic epithelium, our present study focuses on the generation of the first cell-type-specific genome-scale metabolic model (GEM) of human colonic epithelial cells, named iColonEpithelium. GEMs are powerful tools for exploring reactions and metabolites at the systems level and predicting the flux distributions at steady state. Our cell-type-specific iColonEpithelium metabolic reconstruction captures genes specifically expressed in the human colonic epithelial cells. iColonEpithelium is also capable of performing metabolic tasks specific to the colonic epithelium. A unique transport reaction compartment has been included to allow for the simulation of metabolic interactions with the gut microbiome. We used iColonEpithelium to identify metabolic signatures associated with inflammatory bowel disease. We used single-cell RNA sequencing data from Crohn's Diseases (CD) and ulcerative colitis (UC) samples to build disease-specific iColonEpithelium metabolic networks in order to predict metabolic signatures of colonocytes in both healthy and disease states. We identified reactions in nucleotide interconversion, fatty acid synthesis and tryptophan metabolism were differentially regulated in CD and UC conditions, relative to healthy control, which were in accordance with experimental results. The iColonEpithelium metabolic network can be used to identify mechanisms at the cellular level, and we show an initial proof-of-concept for how our tool can be leveraged to explore the metabolic interactions between host and gut microbiota.

Humans↗

PubMed2026

Regnase-1-mediated regulation of neutrophils modulates SARS-CoV-2 pneumonia.

The innate immune response to viral infection needs to be tightly regulated to ensure effective pathogen clearance while avoiding excessive immune activation. During SARS-CoV-2 infection, however, the immune system often fails to elicit appropriate responses, resulting in cytokine-release syndrome in patients with COVID-19. In this study, we show that reduced expression of Regnase-1, an RNase that negatively regulates immune cell activation, confers resistance to infection with the mouse-adapted SARS-CoV-2 MA10 strain. In Regnase-1+/- mice, altered neutrophil function contributed to the amelioration of MA10-induced pneumonia. Single-cell RNA sequencing of lung tissue during MA10 infection revealed four distinct neutrophil subsets, and among these, a subset characterized by an interferon-stimulated gene (ISG) signature was decreased in Regnase-1+/- mice. Furthermore, Regnase-1+/- neutrophils exhibited reduced ISG expression without corresponding changes in proinflammatory gene expression. Regnase-1 was found to repress the expression of Tsc22d3, a gene involved in the negative regulation of interferon responses, through its 3' untranslated region. Collectively, these findings suggest that Regnase-1 attenuates resistance to SARS-CoV-2 MA10 infection by promoting excessive interferon responses in neutrophils.

Animals↗

PubMed2025

Divergent PTEN-p53 interaction upon DNA damage in a human thyroid organoid model with germline PTEN mutations.

Germline mutations in the tumor suppressor phosphatase and tensin homolog (PTEN) cause PTEN hamartoma tumor syndrome (PHTS). PHTS is characterized by an elevated lifetime risk of differentiated thyroid cancer (DTC), 30 times higher than the general population. However, only 1 in 3 PHTS patients develop DTC, and it remains unknown whether specific PTEN variants are associated with an increased risk of DTC. PTEN antagonizes the phosphatidylinositol 3-kinase (PI3K)-AKT signaling pathway, a frequently affected pathway in sporadic DTC. PTEN also acts as a guardian of the genome by interacting with other tumor suppressors. Here, we report how ionizing radiation, an environmental tumorigenic contributor, modifies the DNA damage response based on the type of germline PTEN variants. We hypothesized that certain PTEN variants associated with DTC create a pro-oncogenic molecular signature upon radiation-induced DNA damage. DTC-associated (PTEN M134R ) or DTC-non-associated (PTEN G132D ) germline PTEN mutant alleles were introduced into a human induced pluripotent cell (hiPSC) line derived from a healthy donor utilizing CRISPR-Cas9 gene editing technology. We determined radiation-induced transcriptomic changes in functional thyroid organoids induced from wild-type and both heterozygous PTEN mutant hiPSCs. Both bulk and single-cell RNA sequencing data indicated that radiation upregulated the p53 network more potently in the thyroid organoids with PTEN WT/G132D than those with PTEN WT/M134R , which could be mediated by AKT-dependent MDM2 inactivation and PTEN-p53 physical interaction. Our data suggest that the lack of p53 pathway activation through PTEN-p53 network interactions explains why PTEN M134R is a DTC-susceptible variant.

Humans↗

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

Pan-cancer characterization of HMGA1 reveals its oncogenic role in tumor microenvironment and stemness: functional validation in pancreatic cancer migration and invasion.

BACKGROUND: HMGA1 is a chromatin-associated oncogenic factor implicated in tumor progression, epithelial-mesenchymal transition (EMT), stemness, and metastasis. However, its pan-cancer expression and prognostic patterns, epigenetic activation, and relationship with malignant-cell stemness/plasticity and tumor microenvironment (TME) remodeling in pancreatic adenocarcinoma (PAAD) remain incompletely defined. This study aimed to systematically characterize HMGA1 across cancers and clarify its clinical and biological relevance in PAAD. METHODS: Pan-cancer transcriptomic, clinical, genetic, methylation, immune, and stemness data were integrated from multiple public databases. PAAD single-cell RNA sequencing data were analyzed to localize HMGA1 expression, infer malignant-cell pseudotime, calculate a stemness module score, and assess ligand-receptor communication using CellChat. Public HMGA1-knockdown RNA sequencing data were reanalyzed to evaluate transcriptional remodeling. The Cancer Genome Atlas (TCGA)-PAAD expression and methylation data were used to assess TME-remodeling, immune-suppression, stemness/plasticity, cytokine/chemokine, checkpoint, and promoter-methylation features. HMGA1 expression and function were further examined using immunohistochemistry (IHC), quantitative real-time polymerase chain reaction, Western blotting, wound-healing assays, and Transwell migration and invasion assays. RESULTS: HMGA1 was upregulated in most tumor types, and high expression was associated with unfavorable survival in multiple cancers, including PAAD. In PAAD, HMGA1 was enriched in malignant epithelial cells and positively correlated with pseudotime (Spearman's rho =0.594), while the stemness module score increased along pseudotime (rho =0.748). HMGA1-high malignant cells showed markedly stronger CellChat-inferred outgoing communication, predominantly involving extracellular matrix (ECM)-receptor, adhesion-related, and selected immunomodulatory ligand-receptor axes. HMGA1 knockdown was associated with broad remodeling of EMT, TGF-β, Hedgehog, IL6/JAK/STAT3, KRAS, and cancer stem cell/stemness-related programs rather than uniform suppression of these programs. HMGA1 promoter methylation was inversely correlated with HMGA1 expression (rho =-0.633) and the TME-remodeling score (rho =-0.347). HMGA1 was associated with selected mediators, including PPIA, PLAU, ANXA1, LGALS9, TGFB1, CD276, and CD47, but not with a generalized checkpoint-high phenotype. Functionally, HMGA1 knockdown significantly reduced pancreatic cancer (PC) cell migration and invasion. CONCLUSIONS: These findings support an association-based model in which promoter hypomethylation-associated HMGA1 activation is linked to malignant epithelial stemness/plasticity, ECM/adhesion-dominant TME remodeling, selected immunomodulatory programs, and aggressive PAAD phenotypes. Further mechanistic and clinical validation is required before HMGA1 can be used for therapeutic stratification or immunotherapy-response prediction.

HMGA1↗

PubMed2026

Identification of a novel signature for prognostic stratification and integrative analyses in lung adenocarcinoma.

BACKGROUND: Recently, research has revealed that the Golgi apparatus is involved in the development process of cancer; however, the specific effect of Golgi apparatus-related genes (GAGs) in lung adenocarcinoma (LUAD) remains unclear. This study aims to construct a more concise and practical risk model in LUAD using GAG. METHODS: The gene expression profiles of patients with LUAD were downloaded from The Cancer Genome Atlas (TCGA) and Gene Expression Omnibus (GEO) databases, and GAGs were downloaded from the Gene Set Enrichment Analysis (GSEA) database. Univariate Cox and least absolute shrinkage and selection operator (LASSO) analyses were performed to identify the prognostic GAG signature. Kaplan-Meier and receiver operating characteristic (ROC) curves were plotted to validate the predictive effect of the prognostic signatures. The correlation between the risk model and the immune landscape was examined using CIBERSORT and TIDE analyses. Also, the genes in the signature were assessed by single-cell RNA sequencing (scRNA-seq). RESULTS: A prognostic signature comprising 5 GAG genes (GNPNAT1, RGS20, CAV3, NTSR1, and FURIN) was established after LASSO and multi-Cox analyses. Both the Kaplan-Meier analysis and the ROC curves supported the strong predictive utility of the risk model. Specifically, the former yielded significant stratification in all three validation datasets (P=1.2001e-05, P=0.006, and P=0.04), while the latter provided further evidence of its predictive precision through the area under the curve. In addition, we found that the low-risk group responded better to immunotherapy than the high-risk group (P<0.0001). scRNA-seq analysis revealed the distribution patterns of the 5 GAG genes in cells. Finally, we assessed the situation of tumor mutation burden (TMB) and performed functional analysis based on the risk model of GAGs. CONCLUSIONS: The risk model based on GAGs can effectively stratify the prognosis of patients and predict immunotherapy responses in LUAD.

Golgi apparatus↗