SEARCH · PubMed Health

Results for “gene set enrichment analysis (GSEA)”

Explore indexed PubMed citations for clinical trials, systematic reviews and public health research. Read source abstracts and follow each citation to its original PubMed record.

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

At least 19 recordsLinked to original sources

PubMed2026

High-fat diet-responsive DNM1 promotes hepatocellular carcinoma progression and predicts poor prognosis in viral-associated patients.

Hepatocellular carcinoma (HCC) arises from diverse etiologies, among which metabolic dysfunction-associated liver disease and chronic viral hepatitis are the two major drivers worldwide. However, the molecular mechanisms linking metabolic stress to HCC progression remain incompletely understood. Dynamin-1 (DNM1), primarily known for its role in vesicular trafficking, has emerged as a potential oncogene, yet its prognostic and functional significance in HCC remains largely unexplored. Here, we investigated the role of DNM1 in high-fat diet (HFD)-associated hepatocarcinogenesis. Transcriptomic profiling was conducted to identify differentially expressed genes between normal and high-fat diet murine models, with human orthologs mapped. Clinical relevance was validated using The Cancer Genome Atlas (TCGA-LIHC) dataset. Survival analysis, GSEA (Gene Set Enrichment Analysis), and subgroup stratifications based on viral hepatitis status were performed. In vitro, loss-of-function assays (shRNA knockdown) were executed in HepG2 and SK-Hep1 cell lines to assess cell viability and migration. DNM1 was significantly upregulated in high-fat diet models. In the TCGA-LIHC cohort, high DNM1 expression was an independent risk factor for poor overall survival (HR=1.44, P=0.039) and correlated with advanced tumor stages (Stage III+IV, P=0.010). In vitro knockdown of DNM1 profoundly impaired cell proliferation and migration in HCC cell lines. Strikingly, DNM1 expression was further elevated in patients with concurrent viral hepatitis (P=0.009). GSEA revealed that high DNM1 expression was positively associated with viral infection pathways and negatively correlated with critical immune responses, including interferon-alpha/gamma responses and host immune cytolysis. Survival analysis stratified by four subgroups demonstrated that patients with both viral infection and high DNM1 expression exhibited the worst prognosis (Overall Log-rank P < 0.001). Our findings identify DNM1 as a high-fat diet-responsive regulator that links metabolic stress to hepatocellular carcinoma progression. Elevated DNM1 expression promotes malignant phenotypes in HCC and identifies a subgroup of viral-associated patients with particularly poor prognosis, highlighting DNM1 as a potential prognostic biomarker and therapeutic target.

Hepatocellular carcinoma (HCC)↗

PubMed2026

SERPINE1-centric inflammatory signature associates with treatment resistance and survival in laryngeal squamous cell carcinoma.

BACKGROUND: Laryngeal squamous cell carcinoma (LSCC) prognosis remains poor despite treatment advances. More accurate prognostic assessment models can help guide individualized treatment and improve prognosis. Chronic inflammation contributes to tumorigenesis, yet inflammatory response-related genes (IRGs) in LSCC prognosis are underexplored. This study aimed to construct an IRG prognostic signature for LSCC and further dissect core IRG-mediated mechanisms of immune escape and chemoresistance. METHODS: Transcriptional profiles and clinical data from LSCC patients were retrieved from The Cancer Genome Atlas (TCGA). IRGs were sourced from Gene Set Enrichment Analysis (GSEA) hallmark gene set. We identified differentially expressed IRGs linked to survival outcomes in LSCC. Key IRGs were subsequently selected using least absolute shrinkage and selection operator (LASSO) Cox regression analysis to establish an inflammatory risk score model. This model underwent internal validation within the TCGA cohort and external validation using independent Gene Expression Omnibus (GEO) datasets. We further assessed the model's association with the tumor immune microenvironment and the impact of IRGs on chemotherapy response. Finally, the functional roles of interested signature IRG were experimentally validated in LSCC cell lines. RESULTS: Four significant IRGs (AQP9, ITGA5, LCK, SERPINE1) were identified to build the risk score model. The model stratified LSCC patients into distinct prognostic groups: TCGA cohort: 5-year area under the curve (AUC) =0.836, P<0.001; GSE25727 cohort: 5-year AUC =0.706, P=0.02; GSE27020 cohort: 5-year AUC =0.798, P<0.01. Multivariate analysis confirmed the risk score as an independent prognostic factor (P<0.05). High-risk patients showed reduced immune cell infiltration (CD8+ T cells, dendritic cells) and suppressed immune pathways. Multi-algorithm immune analysis further revealed defective antigen presentation and reduced anti-tumor immune infiltration in high-risk LSCC, promoting tumor immune escape. GSEA/Gene Ontology (GO) enrichment combined with drug sensitivity prediction further revealed that high-risk tumors activate invasive signaling and acquire broad chemoresistance alongside impaired anti-tumor immunity. SERPINE1 might be associated with chemotherapy resistance and exhibited the highest alteration frequency (predominantly amplification) and overexpression in LSCC tissues. Its knockdown significantly suppressed proliferation, migration, invasion and chemoresistance in LSCC cells. Immunohistochemistry (IHC) confirmed tumor SERPINE1 overexpression (P=0.002 vs. normal tissues), correlating with poor survival (P<0.001). CONCLUSIONS: The 4-IRG risk signature is a reliable prognostic indicator reflecting immune dysfunction in LSCC. SERPINE1 is validated as a therapeutic target and biomarker, enriching our understanding of gene regulation dynamics in LSCC.

Laryngeal cancer↗

PubMed2024

Uncovering essential anesthetics-induced exosomal miRNAs related to hepatocellular carcinoma progression: a bioinformatic investigation.

BACKGROUND: Anesthetic drugs may alter exosomal microRNA (miRNA) contents and mediate cancer progression and tumor microenvironment remodeling. Our study aims to explore how the anesthetics (sevoflurane and propofol) impact the miRNA makeup within exosomes in hepatocellular carcinoma (HCC), alongside the interconnected signaling pathways linked to the tumor immune microenvironment. METHODS: In this prospective study, we collected plasma exosomes from two groups of HCC patients (n = 5 each) treated with either propofol or sevoflurane, both before anesthesia and after hepatectomy. Exosomal miRNA profiles were assessed using next-generation sequencing (NGS). Furthermore, the expression data from The Cancer Genome Atlas-Liver Hepatocellular Carcinoma (TCGA-LIHC) was used to pinpoint the differentially expressed exosomal miRNAs (DEmiRNAs) attributed to the influence of propofol or sevoflurane in the context of HCC. Gene set enrichment analysis (GSEA) and gene set variation analysis (GSVA) were used to dissect the signaling pathways and biological activities associated with the identified DEmiRNAs and their corresponding target genes. RESULTS: A total of 35 distinct DEmiRNAs were exclusively regulated by either propofol (n = 9) or sevoflurane (n = 26). Through TCGA-LIHC database analysis, 8 DEmiRNAs were associated with HCC. These included propofol-triggered miR-452-5p and let-7c-5p, as well as sevoflurane-induced miR-24-1-5p, miR-122-5p, miR-200a-3p, miR-4686, miR-214-3p, and miR-511-5p. Analyses revealed that among these 8 DEmiRNAs, the upregulation of miR-24-1-5p consistently demonstrated a significant association with lower histological grades (p < 0.0001), early-stage tumors (p < 0.05) and higher survival (p = 0.029). Further analyses using GSEA and GSVA indicated that miR-24-1-5p, along with its target genes, were involved in governing the tumor immune microenvironment and potentially inhibiting tumor progression in HCC. CONCLUSIONS: This study provided bioinformatics evidence suggesting that sevoflurane-induced plasma exosomal miRNAs may have a potential impact on the immune microenvironment of HCC. These findings established a foundation for future research into mechanistic outcomes in cancer patients.

Carcinoma, Hepatocellular↗

PubMed2025

Bioinformatics identification and validation of pyroptosis-related gene for ischemic stroke.

BACKGROUND: Ischemic stroke (IS) is one of the common and frequent diseases with extremely high lethality and disability in the world, and there is no effective treatment at present. This study aimed to screen hub genes involved in cerebral ischemia/reperfusion injury (CIRI) and pyroptosis, and explore promising intervention targets. METHODS: CIRI-related genes (GSE202659 and GSE131193) and pyroptosis-related genes (PRGs) in mice were obtained from the Gene Expression Omnibus (GEO) and GeneCards database. We screened for LASSO regression to construct a prognostic model of GSE131193 and PRGs and examined by GSE137482. The functional enrichment analysis of Gene Ontology (GO), Kyoto Encyclopedia of Genes and Genomes (KEGG), Gene Set Enrichment Analysis (GSEA) and Gene Set Variation Analysis (GSVA) were performed on pyroptosis-related differentially expressed genes (PRDEGs) of GSE202659.The key modules for CIRI and pyroptosis were identified by Weight Gene Co-expression Network Analysis (WGCNA). Subsequently, Protein-protein Interaction (PPI) network and the Cytoscape was constructed to screen out hub genes. Used the starBase to predict miRNA interacting with hub genes and constructed mRNA-miRNA-lncRNA interaction networks. CIRI-related Molecular Subtypes were constructed for hub genes. The relationship between immune cells and hub genes was verified via CIBERSORT. Finally, we selected C57BL/6 mice to construct models to confirm hub genes by enzyme linked immunosorbent assay (ELISA), reverse transcription-polymerase chain reaction (RT-PCR), western blot, and Immunofluorescence. RESULTS: A total of 272 PRGs and 35 PRDEGs were screened. An eight-gene risk prediction models were established (AUC = 0.868). GO, KEGG, GSEA and GSVA analyses revealed that PRDEGs were mainly involved in positive regulation of cytokine production, and NOD-like receptor signaling pathway. And then, seven hub genes (Irf1, Icam1, Tlr2, Tnf, Cebpb, Il1rn, and Casp8) were identified by PPI. Icam1, Tnf, Cebpb, Il1rn, and Casp8 had high expression profiles in Cluster2 by hierarchical clustering. The immune infiltration analysis results showed that among the hub genes, Cebpb, Il1rn, and Casp8, showed a significant positive correlation with the degree of NK.Actived, and Icam1 showed a significant negative correlation with B.Cells.Memory. The results of animal experiments significantly demonstrated an upregulation of Irf1, Icam1, Tlr2, Cebpb, and Il1rn. CONCLUSION: Our finding indicated that Irf1, Icam1, Tlr2, Cebpb, and Il1rn are hub genes associated with pyroptosis, and these genes are all associated with different immune cells, so as to provide new targets for the prevention and treatment of IS from the perspective of pyroptosis.

Pyroptosis↗

PubMed2026

ERP44 Is Associated With Poor Prognosis and Promotes Proliferation and Temozolomide Resistance in Lower-grade Glioma.

BACKGROUND/AIM: Endoplasmic reticulum resident protein 44 (ERP44), a protein disulfide isomerase family member, has been implicated in tumor biology, but its role in lower-grade glioma (LGG) remains unclear. This study investigated the prognostic significance and biological function of ERP44 in LGG, focusing on proliferation and temozolomide (TMZ) resistance. MATERIALS AND METHODS: ERP44 expression, clinicopathological associations, and prognostic value were analyzed using The Cancer Genome Atlas (TCGA), Genotype-Tissue Expression (GTEx), and Chinese Glioma Genome Atlas (CGGA) datasets. Time-dependent receiver operating characteristic (ROC) curves, Cox regression, and a prognostic nomogram were constructed. Differential expression, Gene Set Enrichment Analysis (GSEA), Gene Ontology (GO) enrichment, immune infiltration, and drug sensitivity analyses were performed. Functional validation was conducted in SW1088 and SW1783 cells using shRNA-mediated ERP44 knockdown, followed by RT-qPCR, western blotting, CCK-8, colony formation, and TMZ IC50 assays. Subcutaneous xenograft models with or without TMZ treatment were used for in vivo validation. RESULTS: ERP44 was markedly upregulated in LGG and associated with higher WHO grade, IDH wildtype status, 1p/19q non-codeletion, and poor survival in TCGA and CGGA cohorts. ERP44 showed strong prognostic performance and improved risk stratification in a multivariable nomogram. Enrichment analyses linked high ERP44 expression to immune/inflammatory pathways and reduced neuronal functional signatures. ERP44 positively correlated with immune infiltration, proliferation/stemness markers, and predicted TMZ resistance, while its knockdown inhibited proliferation and colony formation, reduced TMZ IC50, suppressed xenograft growth, enhanced TMZ efficacy, and decreased Ki67 positivity. CONCLUSION: ERP44 is a prognostic biomarker that promotes LGG proliferation and TMZ resistance, suggesting its potential as a therapeutic target.

Humans↗

PubMed2026

Genes near tRNAs are enriched in translational machinery.

Transfer RNAs (tRNAs) are known for delivering amino acids to the growing polypeptide chain during translation. They can also influence gene expression, especially in times of nutrient starvation, through differential tRNA expression and modification. Transfer RNAs have a highly consistent cloverleaf structure, but relatively few known regulatory elements govern this conserved structure despite the 20 different standard isotypes. This study examines gene enrichment patterns near tRNA genes across 1149 fungal genomes. Genes enriched in proteasome regulation, ion transport, and rRNA were found to be significantly closer to tRNAs than other pathways. These results were consistent across KEGG overrepresentation analysis (ORA), KEGG gene set enrichment analysis (GSEA), and gene ontology (GO) analysis. Proteasome, ion transport, and RNA are all important aspects of protein production and regulation, suggesting that genes required for the synthesis and quality control of proteins, including tRNAs, are located near each other. Protein regulation is an energetically expensive process, and local co-regulation could increase efficiency and stress impacts on proteins.

RNA, Transfer↗

PubMed2025

Stage-dependent proteomic alterations in aqueous humor of diabetic retinopathy patients based on data-independent acquisition and parallel reaction monitoring.

BACKGROUND: Diabetic retinopathy (DR), a microvascular complication of diabetes mellitus (DM), represents the predominant cause of preventable vision loss in working-age populations globally. While the pathophysiological mechanisms underlying DR progression remain incompletely understood, our study employs comprehensive proteomic profiling of aqueous humor (AH) to identify stage-specific biomarkers and therapeutic targets in type 2 diabetes mellitus (T2DM) patients across DR progression. METHODS: Utilizing data-independent acquisition (DIA) mass spectrometry, we quantified AH proteomes in a discovery cohort comprising 24 subjects: 18 T2DM patients stratified by DR severity [6 non-DR, 6 non-proliferative DR (NPDR), 6 proliferative DR (PDR)] and 6 cataract controls without diabetes (non-DM). Validation cohort analysis (including 10 AH samples in each group) was performed using parallel reaction monitoring (PRM) strategy for verification of target proteins. Comprehensive bioinformatics analyses included gene set enrichment analysis (GSEA), weighted gene co-expression network analysis (WGCNA), Kyoto encyclopedia of genes and genomes (KEGG) enrichment analysis, protein-protein interaction (PPI) network construction, receiver operating characteristic (ROC) curve analysis, and ConnectivityMap (Cmap)-based drug prediction. RESULTS: Proteomic profiling identified 739 quantifiable AH proteins (62% extracellular) with clear separation among the four clinical stages in the discovery cohort. GSEA uncovered altered expression of proteins mainly related to complement and coagulation cascades, folate metabolism, and the selenium micronutrient network in patients with DR. WGCNA-derived protein modules yielded 83 PRM-validated targets, including 5 hub proteins differentiating NPDR from non-DR and 33 hub proteins showed significant upregulation in PDR versus NPDR comparison. Clinical correlation analysis identified F2, FGG, FGB, RBP4, AMBP, VTN, C8A, CPB2, and C2 associated with clinical traits. C6, FAM3C, SPP1, and JCHAIN levels were altered post-anti-VEGF treatment. Pharmacological prediction identified potential therapeutic compounds, including perindopril, triciribine, and XAV-939 for NPDR, and topiramate, triciribine, and vecuronium for PDR. CONCLUSION: This study established a comprehensive AH proteomic signature of DR progression, offering insights into the pathogenesis of DR and highlighting potential biomarkers and novel therapeutic targets.

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

NR3C1 Modulates Wnt Signalling to Influence the Invasiveness and Immune Features of Nonfunctioning Invasive Pituitary Adenomas.

Pituitary adenomas (PAs) are common intracranial tumours, and invasiveness in nonfunctioning invasive pituitary adenomas (NIPAs) predicts poor prognosis. The molecular mechanisms driving this phenotype remain unclear. This study explored the role of nuclear receptor subfamily 3 group C member 1 (NR3C1) in NIPA invasiveness and its regulation of Wnt signalling. mRNA expression profiles of 32 PA samples were generated by RNA-seq, and proteomic data from 19 samples were obtained by mass spectrometry. Immune-related differentially expressed genes (DEGs) were retrieved from GeneCards. Weighted gene coexpression network analysis identified modules and hub genes linked to invasiveness, while machine learning methods (support vector machine, LASSO, random forest) prioritised key genes. Gene set enrichment analysis (GSEA) assessed pathways associated with candidate gene expression. NR3C1 expression and function were validated by immunohistochemistry, Western blotting and invasion assays. Integration of transcriptomic, proteomic and immune-related datasets yielded 11 overlapping genes, with NR3C1 emerging as the top candidate. NR3C1 was significantly upregulated in NIPAs and demonstrated good discriminatory power by ROC analysis. GSEA associated high NR3C1 expression with Wnt pathway activation. Functional experiments confirmed that NR3C1 overexpression enhances the invasive capacity of PA cells. NR3C1 promotes the invasive phenotype of NIPAs by activating Wnt signalling. These findings suggest NR3C1 as a potential biomarker and therapeutic target for invasive pituitary adenomas.

Humans↗

PubMed2024

Assessment of Gene Set Enrichment Analysis using curated RNA-seq-based benchmarks.

Pathway enrichment analysis is a ubiquitous computational biology method to interpret a list of genes (typically derived from the association of large-scale omics data with phenotypes of interest) in terms of higher-level, predefined gene sets that share biological function, chromosomal location, or other common features. Among many tools developed so far, Gene Set Enrichment Analysis (GSEA) stands out as one of the pioneering and most widely used methods. Although originally developed for microarray data, GSEA is nowadays extensively utilized for RNA-seq data analysis. Here, we quantitatively assessed the performance of a variety of GSEA modalities and provide guidance in the practical use of GSEA in RNA-seq experiments. We leveraged harmonized RNA-seq datasets available from The Cancer Genome Atlas (TCGA) in combination with large, curated pathway collections from the Molecular Signatures Database to obtain cancer-type-specific target pathway lists across multiple cancer types. We carried out a detailed analysis of GSEA performance using both gene-set and phenotype permutations combined with four different choices for the Kolmogorov-Smirnov enrichment statistic. Based on our benchmarks, we conclude that the classic/unweighted gene-set permutation approach offered comparable or better sensitivity-vs-specificity tradeoffs across cancer types compared with other, more complex and computationally intensive permutation methods. Finally, we analyzed other large cohorts for thyroid cancer and hepatocellular carcinoma. We utilized a new consensus metric, the Enrichment Evidence Score (EES), which showed a remarkable agreement between pathways identified in TCGA and those from other sources, despite differences in cancer etiology. This finding suggests an EES-based strategy to identify a core set of pathways that may be complemented by an expanded set of pathways for downstream exploratory analysis. This work fills the existing gap in current guidelines and benchmarks for the use of GSEA with RNA-seq data and provides a framework to enable detailed benchmarking of other RNA-seq-based pathway analysis tools.

Humans↗

PubMed2025

Expression and prognosis of CXCL13 in uterine corpus endometrial carcinoma based on bioinformatics analysis.

OBJECTIVE: The biological significance of the chemokine ligand C-X-C motif chemokine ligand 13 (CXCL13) may play a significant role in the pathogenesis of uterine corpus endometrial carcinoma (UCEC). This study aims to identify and verify CXCL13 with predictive value for prognosis in UCEC. METHODS: CXCL13 mRNA expression differences were analyzed using R software in three independent datasets: one each from The Cancer Genome Atlas (TCGA) and two from the Gene Expression Omnibus (GEO), namely GSE17025 and GSE106191. The correlation between CXCL13 expression and prognosis was evaluated by Kaplan-Meier analysis. Univariate and multivariate Cox analyses were utilized to construct a prognostic nomogram. Tumor Immune Estimation Resource (TIMER) and the Tumor and Immune System Interaction Database (TISIDB) were employed to assess the relationship between CXCL13 and tumor immune infiltration. Coexpressed genes with CXCL13 were identified by the Spearman correlation analysis. A CXCL13 protein-protein interaction (PPI) network was constructed with the STRING website tool and hub genes were screened out. Gene Ontology (GO) and Kyoto Encyclopedia of Genes and Genome (KEGG) analyses were performed with the "clusterProfiler" R package. Gene set enrichment analysis (GSEA) was used to identify underlying biological mechanisms. A drug-gene interaction network was constructed in the Comparative Toxicogenomics Database (CTD). RESULTS: High CXCL13 mRNA expression were validated in UCEC in the above three independent datasets. High CXCL13 expression was associated with favorable prognosis in UCEC. A nomogram for predicting the 1-, 3-, and 5-year survival probability in UCEC was construct based on CXCL13 expression and other clinical parameters. The use of Spearman correlation indicated certain correlation between CXCL13 and immune cells and immune checkpoint (ICP) genes. Seven hub genes were upregulated in UCEC, namely CXCL9, IFNG, CXCL10, CXCL11, GBP5, CCL18, and GZMB. The expression and prognostic relevance of CXCL9, IFNG, GBP5, and GZMB were in accordance with CXCL13. The main biological processes enriched were cytokine-cytokine receptor interaction and chemokine signaling pathway. CONCLUSIONS: The above comprehensive analyses suggest that CXCL13 may serve as a potential prognostic biomarker for UCEC, specifically for early-stage UCEC.

CXCL13↗

PubMed2026

Phenotypic, physiological and transcriptomic analysis of graded salt stress responses in Pyrus betulifolia Bunge and functional characterization of the hub gene PbSTY46.

Pyrus betulifolia Bunge is a salt‑tolerant rootstock for pear, but its salt‑tolerance mechanisms remain largely unknown. In this study, P. betulifolia seedlings were subjected to graded NaCl stress at concentrations of 0 (CK), 50 (T1), 100 (T2), and 200 (T3) mM. We integrated phenotypic observation, physiological assessment, transcriptomic profiling, and functional gene validation to systematically elucidate its salt tolerance mechanisms. Salt stress inhibited seedling growth and root traits in a concentration-dependent manner, and T3 caused the most severe damage. Osmotic solutes responded differentially: soluble sugars peaked under T2, while proline peaked under T3. Antioxidant enzymes showed tissue-specific biphasic responses and declined after prolonged T3 stress. Meanwhile, chlorophyll and photosynthesis decreased, whereas anthocyanin increased, indicating a metabolic shift from photosynthesis to photoprotection. Transcriptome analysis revealed distinct responses depending on stress intensity: mild stress induced membrane lipid remodeling, moderate stress activated circadian rhythm and hormone signaling, and severe stress enhanced phenylpropanoid biosynthesis and thiamine metabolism. Gene Set Enrichment Analysis (GSEA) further highlighted progressive enrichment of phenylpropanoid biosynthesis, heme binding, and oxidoreductase activity. Weighted Gene Co‑expression Network Analysis (WGCNA) identified a blue module significantly positively correlated with root traits, from which the hub gene PbSTY46 was identified. Functional validation via overexpression, loss‑of‑function mutants, and pharmacological interventions (MeJA/DIECA) confirmed that PbSTY46 acts through JA signaling to enhance antioxidant enzyme activities and thereby confer salt tolerance. Collectively, P. betulifolia adopts a "survival‑first" strategy that coordinates growth arrest, osmotic homeostasis, and ROS scavenging. These findings establish PbSTY46 as a key regulator that links JA signaling to antioxidant defense. Thus, PbSTY46 represents a promising candidate for marker‑assisted breeding of salt‑tolerant pear cultivars.

Salt Stress↗

PubMed2026

Transcriptomic analysis identifies novel ferroptosis-related biomarkers and therapeutic targets in pulmonary arterial hypertension.

BACKGROUND: Ferroptosis plays a significant role in pulmonary arterial hypertension (PAH), although its underlying mechanisms and key pathogenic genes remain unclear. METHODS: Transcriptomic data from human PAH and control lung tissue were obtained from the Gene Expression Omnibus (GEO) database, whereas ferroptosis-related genes (FRGs) were sourced from the MsigDb and FerrDb databases. Differentially expressed FRGs (DE-FRGs) were identified through the intersection of FRGs with differentially expressed genes (DEGs). Functional enrichment analysis was performed using Gene Ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) pathways. Key hub genes were identified through Least Absolute Shrinkage and Selection Operator (LASSO), support vector machine-recursive feature elimination (SVM-RFE), and weighted correlation network analysis (WGCNA). Gene set enrichment analysis (GSEA) was conducted to explore the functional roles and associated pathways of hub genes. The relationship between hub genes and immune infiltration was investigated. Expression levels of potential biomarkers were validated via Quantitative real-time polymerase chain reaction (qRT-PCR) and immunohistochemistry (IHC) in two PAH animal models (monocrotaline-induced and Sugen5416 plus hypoxia-induced PAH). Finally, molecular docking was employed to screen potential therapeutic compounds. RESULTS: A total of 133 DE-FRGs were identified, with KEGG and GO analyses highlighting their involvement in intracellular iron homeostasis and ferroptosis. Hub genes, notably FZD7 and NFE2, were identified using LASSO, SVM-RFE, and WGCNA. Immune infiltration analysis suggested that monocytes and neutrophils play key roles in PAH pathogenesis. Validation in PAH animal models showed significant upregulation of Fzd7 and downregulation of Nfe2 in lung tissues of both MCT- and SuHx-induced PAH models. Molecular docking identified tetrachlorodibenzodioxin (TCDD) has good binding affinity. CONCLUSION: In summary, we investigated two ferroptosis-related biomarkers, FZD7 and NFE2, in PAH using transcriptomics, offering new insights into molecular mechanisms and potential targeted therapies for the disease.

Ferroptosis↗

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

Targeting RELA and STAT3 regulates TNFRSF10A-mediated apoptosis in a novel apoptosis-based prognostic model for clear cell renal cell carcinoma.

BACKGROUND: Clear cell renal cell carcinoma (ccRCC) is the most common subtype of renal malignancy and remains a major cause of cancer-related mortality worldwide. Although advances in surgery, targeted therapy, and immunotherapy have improved outcomes for patients, reliable biomarkers for predicting prognosis remain limited. Therefore, robust gene-based prognostic models are urgently needed to improve risk stratification and guide individualized treatment strategies. METHODS: We developed a novel prognostic model integrating apoptosis and immune - related genes (AIRGs) to predict overall survival (OS) in patients with ccRCC. RESULT: Using Gene Set Enrichment Analysis (GSEA) combined with least absolute shrinkage and selection operator (LASSO) Cox regression, we identified 7 key prognostic genes, namely, CCR4, TNFRSF10A, TEK, TGFA, CD14, IFITM1, and SEMA3G, that collectively demonstrated strong predictive performance in TCGA cohort with c-index = 0.711. Functional enrichment analyses revealed that apoptosis, immune regulation, and multiple oncogenic signaling pathways were significantly associated with the risk score, highlighting the critical role of the tumor microenvironment in ccRCC progression. Transcription factor binding analysis based on the JASPAR database suggested that RELA and STAT3 with scores of 0.829 and 0.951, respectively are potential upstream regulators within the prognostic network, particularly influencing TNFRSF10A expression. External validation using the International Cancer Genome Consortium (ICGC) dataset confirmed the robustness of the prognostic model with c-index = 0.612 Furthermore, in vitro experiments demonstrated that RELA and STAT3 regulate TNFRSF10A-mediated apoptotic signaling in ccRCC cells, providing mechanistic support for the bioinformatic findings. CONCLUSION: This study establishes a biologically informed and clinically relevant prognostic framework for ccRCC. Our findings highlight the therapeutic potential of targeting the RELA/STAT3-TNFRSF10A axis and contribute to the advancement of precision medicine in ccRCC.

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↗

PubMed2026

Evaluation of the subtype-specific epigenetic prognostic association of HELLS in non-small cell lung cancer: integrated clinical and molecular insights.

BACKGROUND: Helicase, lymphoid-specific (HELLS) is an epigenetic chromatin remodeler implicated in several cancers, but its prognostic role in non-small cell lung cancer (NSCLC) subtypes remains unclear. We investigated the expression, prognostic significance, and subtype-specific associations of HELLS in lung adenocarcinoma (LUAD) and lung squamous cell carcinoma (LUSC). METHODS: The Cancer Genome Atlas (TCGA) and independent Gene Expression Omnibus (GEO) datasets were analyzed. HELLS expression was compared between tumor and normal tissues, survival was evaluated separately in LUAD and LUSC, and gene set enrichment analysis (GSEA) was performed. Multivariable analyses were used to assess associations between HELLS and selected oncogenic and immune-related genes after adjustment for clinical variables. RESULTS: HELLS was significantly upregulated in both LUAD and LUSC compared with normal lung tissues (P<0.001). High HELLS expression was associated with shorter overall survival (OS) in LUAD (log-rank P=0.001) and in the TCGA-LUSC cohort (log-rank P=0.002); however, external validation in GSE42127 (LUSC, n=43) was not significant [log-rank P=0.12; hazard ratio (HR) =0.49, 95% confidence interval (CI): 0.20-1.22, P=0.13]. HELLS-high LUAD tumors showed enrichment trends enriched in proliferation-related pathways, whereas HELLS-low LUSC tumors were enriched in inflammatory and apoptotic pathways. HELLS expression remained associated with KRAS, BRAF, and CD274 in LUAD after adjustment for age, sex, and stage, while only limited associations were observed in LUSC. CONCLUSIONS: HELLS shows a subtype-dependent prognostic and molecular association in NSCLC, with the strongest and most reproducible signal in LUAD; however, its prognostic value is attenuated after multivariable adjustment and is not consistently reproduced across external cohorts.

Helicase, lymphoid-specific (HELLS)↗

PubMed2026

Transcriptomic Profiling of Canine Testicular Leydig Cell Tumors Uncovers Key Upregulated Gene Pathways.

Total RNA was isolated from sections of healthy testes and Leydig cell tumors of mixed-breed dogs using TMA Master II device. The RNA-seq libraries were sequenced on the Illumina platform. Following differential expression analysis, Gene Ontology (GO), Kyoto Encyclopedia of Genes and Genomes (KEGG), and Gene Set Enrichment Analysis (GSEA) were applied with quality control obtained using FastQC and Trimmomatic. This analysis revealed 1500 transcripts, including 928 upregulated and 168 downregulated genes. The results demonstrated that a significant proportion of these differentially expressed genes are directly involved in the control of sex steroid production (CYP11A1, STAR, and 3β-HSD3B1) or tube formation, angiogenesis, and extracellular matrix remodeling in interstitial cells (ESM1, FGG, and VEGFA). Moreover, we identified the upregulation of transcripts responsible for neurotransmitter or neuroendocrine signaling (SLC6A4, GRIN2C, GABRB3) and cholesterol metabolism and its regulation (GPX3, MSMO1, DHCR24). These genes were strongly associated with the phosphatidylinositol-3-kinase (PI3K)-Protein Kinase B (Akt) cascade and extracellular matrix interactions, features shared with various malignancies. Alterations in estrogen and relaxin signaling appear to be distinctive, understudied mechanisms specific to canine Leydig cell tumors. Concurrently, downregulated genes (e.g., DMRTC2, SEMA3C, ALOX12) were linked with cell differentiation, signaling and immunoregulatory pathway suppression involved in tumorigenesis. A complex transcriptomic profile of canine Leydig cell tumors was developed, revealing a conserved oncogenic core shared in some aspects with human malignancies alongside unique species-specific alterations. Findings seem to be useful for identifying novel diagnostic biomarkers and targeted therapies in veterinary oncology, establishing canine reproductive tissues as a valuable comparative biomedical model for research in human.

Leydig cell tumor↗