PubMed HealthSearch

SEARCH · PubMed Health

Results for “gene set enrichment analysis”

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

Research on identification of key genes and immune-metabolic mechanisms in atrial fibrillation through integrated multi-cohort transcriptomic analysis and machine learning.

This study aimed to integrate multiple datasets for the identification of atrial fibrillation (AF)-related differentially expressed genes (DEGs), analyze their underlying mechanisms through functional enrichment and machine learning, construct diagnostic models, and explore immune-metabolic interactions to provide novel biomarkers and theoretical foundations. Gene expression datasets were integrated and normalized, with batch effects removed using principal component analysis. Differential expression analysis, functional enrichment analysis (Gene Ontology and Kyoto Encyclopedia of Genes and Genomes pathways), and machine learning-based feature gene selection and model construction were performed. Shapley additive explanations analysis was utilized to interpret the constructed models, while gene set enrichment analysis, gene set variation analysis, and immune cell infiltration analysis were conducted to investigate the associations between feature genes and immune infiltration. After integrating and normalizing gene expression data and eliminating batch effects via principal component analysis, 6 DEGs were identified, including 4 upregulated and 2 down-regulated ones. Functional enrichment analysis showed these DEGs were significantly enriched in neuro-related biological processes and pathways, indicating their key roles in AF pathogenesis. Five key feature genes were selected using LASSO, random forest, and support vector machine-recursive feature elimination algorithms. They had significant expression differences between the AF and control groups (P&#x2005;<&#x2005;.001) and were located on distinct chromosomes. The constructed random forest and support vector machine models performed excellently (area under the curve&#x2005;&#x2265;&#x2005;0.85). Shapley additive explanations analysis revealed TNNI1 contributed most to model prediction, with its expression significantly positively correlated with immune cell infiltration. Gene set enrichment analysis and gene set variation analysis analyses further showed feature genes participated in AF pathogenesis by regulating immune modulation, metabolic pathways, and autophagy. Immune cell infiltration analysis found altered proportions of T-cell subsets and M0 macrophages in the AF group, along with complex links between feature gene expression and immune cell function. This study systematically elucidated the unique gene expression patterns and key regulatory pathways associated with AF, clarifying the crucial roles of feature genes in immune regulation, metabolic imbalance, and cellular dysfunction. These findings provide a theoretical basis and potential therapeutic targets for understanding AF pathogenesis and developing targeted treatment strategies.

Atrial Fibrillation

A multifaceted investigation into the impact of m6A methylation-related genes on pancreatic cancer, integrating insights from various databases and foundational experimental research.

BACKGROUND: Despite advances in surgical techniques, immunotherapy, the mortality rate associated with pancreatic cancer (PC) has been on the rise in recent years. Understanding the importance of RNA N6-methyladenosine (m6A) in PC is critical for prognosis, tumor microenvironment, and immunotherapy efficacy. The study aims to identify m6A methylation regulators that play an important role in the development and progression of PC by mining databases. The effect of insulin-like growth factor-binding protein 3 (IGFBP3) on pancreatic tumors was explored, and the related mechanisms were explored. METHODS: We analyzed the expression of m6A regulators in PC by digging deeper into the datasets of The Cancer Genome Atlas and Gene Expression Omnibus (GEO) databases, and analyzed its relationship with the prognosis of patients with PC, looking for m6A methylation regulators that play an important role in the development and progression of PC. Reuse the ConsensusClusterPlus package, Cox analysis, and unsupervised clustering to delineate three distinct m6A clusters - designated as m6A cluster A, m6A cluster B, and m6A cluster C single-sample gene set enrichment analysis, gene set variation analysis, Gene Ontology, and Kyoto Encyclopedia of Genes and Genomes (KEGG) analyses evaluated the different pathway roles of these clusters in the development and progression of PC. Finally, the cell lines with IGFBP3 overexpression and knockdown were constructed by lentivirus transfection, the transfection effect was identified by WB, and the effects of IGFBP3 overexpression/knockdown on the survival and growth of PC cell lines were verified by cell cloning experiments and cell counting kit-8 experiments, and the possible related pathways were explored by KEGG. RESULTS: Most m6A regulatory factors are highly expressed in PC, and their high expression is negatively correlated with the prognosis of patients with PC. Furthermore, m6A regulatory factors may influence the occurrence and development of PC through metabolic pathways, stroma activation pathways, immune regulatory processes, and the immune microenvironment. Finally, the overexpression of IGFBP3 promoted the growth of PC cells, and vice versa. CONCLUSIONS: Most m6A regulatory factors are differentially expressed in PC and are associated with the prognosis of patients with PC, potentially influencing the occurrence and development of PC through pathways such as the immune microenvironment. The overexpression of IGFBP3 can promote the growth of PC cells and vice versa.

IGFBP3

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

Inflammatory pathways and immune dysregulation in pediatric postoperative septic shock: A study integrating transcriptomics, machine learning and molecular docking.

This study elucidates the molecular and immune regulatory mechanisms of pediatric postoperative septic shock. Transcriptomic data were obtained from the Gene Expression Omnibus database. Differentially expressed genes were identified using the limma package, and gene co-expression modules were constructed using Weighted Gene Co-expression Network Analysis. Functional enrichment was performed via gene set enrichment analysis, Gene Ontology, and Kyoto Encyclopedia of Genes and Genomes analyses. Immune cell infiltration was assessed using ESTIMATE and CIBERSORT. Mendelian randomization was applied to explore causal relationships between gene expression and septic shock. Feature genes were selected using machine learning algorithms, and a diagnostic nomogram model was constructed. Finally, molecular docking analysis was performed to screen and evaluate the binding affinity of traditional Chinese medicine monomers to core target proteins. A total of 1331 differentially expressed genes were identified, and the turquoise module was strongly correlated with septic shock. Enrichment analysis revealed significant activation of IL-6/JAK/STAT3, TNF-&#x3b1;/NF-&#x3ba;B, and PI3K/Akt/mTOR pathways. Immune infiltration analysis indicated suppressed immune scores and imbalances in neutrophils, macrophages, T cells, and B cells. Mendelian randomization confirmed causal associations for 6 genes, including PIM3. The predictive model based on feature genes demonstrated high diagnostic performance. Molecular docking suggested that quercetin and astramembrannin I could stably bind PIM3. This study systematically identified core genes, dysregulated immune pathways, and candidate small-molecule interventions in pediatric septic shock, providing novel insights for early diagnosis and targeted therapy.

Humans

POU2F3 expression in lung squamous cell carcinoma: transcriptomic and immunohistochemical profiling with prognosis.

BACKGROUND: Lung squamous cell carcinoma (LUSC) lacks well-defined molecular targets. This study investigated the clinical and biological relevance of POU class 2 homeobox 3 (POU2F3), a tuft cell-associated transcription factor, in LUSC. METHODS: RNA sequencing data of patients with LUSC from The Cancer Genome Atlas (TCGA cohort, n&#xa0;=&#xa0;190) was analysed and compared to a cohort of surgically resected cases analyzed via immunohistochemistry (IHC cohort, n&#xa0;=&#xa0;137). Prognostic impact was assessed via survival analyses. Transcriptomic features, pathway enrichment, and immune profiles were evaluated via differentially expressed gene analysis, Gene Set Enrichment Analysis, and CIBERSORTx. RESULTS: High POU2F3 expression independently predicted poor overall survival in the TCGA cohort (HR&#xa0;=&#xa0;2.06, 95% CI: 1.04-4.08, P&#xa0;=&#xa0;0.039). In contrast, POU2F3 expression was not prognostic in the IHC cohort (P&#xa0;=&#xa0;0.995). Morphologically, POU2F3-positive tumours were enriched for non-keratinizing and poorly differentiated subtypes. Transcriptomic analysis showed suppression of proliferation and immune-related pathways (FDR&#xa0;<&#xa0;0.001), with suggestive enrichment of the TGF-&#x3b2; (FDR&#xa0;=&#xa0;0.143) and p53 (FDR&#xa0;=&#xa0;0.229) signaling pathways. On immune deconvolution, POU2F3-high tumours showed a nominal increase in activated dendritic cells, which did not withstand multiple testing correction. POU2F3 protein was detected in 12.4% of tumours and was significantly associated with p53 or RB1 abnormalities (single or double) (P&#xa0;=&#xa0;0.028). CONCLUSIONS: POU2F3 marks a transcriptionally distinct, early-stage subtype of LUSC with keratinization-related features. Its prognostic relevance appears context-dependent and requires prospective validation in uniformly treated cohorts.

Humans

Integrative Multi-Omics Analysis of Stem Growth Habit Divergence in Wild Soybean (Glycine soja).

Stem architecture is a major determinant of lodging resistance, biomass accumulation, and harvest efficiency in soybean. However, the molecular features associated with contrasting stem growth habits in wild soybean remain incompletely characterised. Here, we performed an integrated transcriptomic, metabolomic, and epigenomic analysis of stem growth-habit divergence in wild soybean, comparing the wild-type accession ZYD7068 with contrasting vining and erect mutant lines derived from carbon-ion beam mutagenesis. Pairwise transcriptomic comparisons identified between 20&#x2009;311 and 28&#x2009;705 differentially expressed genes per contrast, with a core set of 2672 genes consistently altered across the comparisons. Functional enrichment, gene set variation analysis, and gene set enrichment analysis converged on xylem and phloem pattern formation as a prominent molecular pathway associated with growth-habit divergence. Random forest analysis identified BBR-BPC and ARF transcription factor families as major molecular discriminators, while metabolomic profiling revealed distinct metabolic profiles involving amino-acid-derived and lipid-associated metabolites. Whole-genome bisulfite sequencing revealed context-specific DNA methylation differences, including substantial variation in CHG methylation among erect mutant lines. Integrated network and in silico perturbation analyses prioritised four candidate genes associated with vascular development for future functional validation. Together, these results provide a multi-layer molecular resource for investigating stem growth-habit divergence in G. soja and establish testable candidate pathways and genes for subsequent functional studies and soybean improvement.

glycine soja

Transcriptomic Profiling Reveals NF-&#x3ba;B-Associated Immune Regulatory Signatures Underlying the Regenerative Effects of Hypoxia-Preconditioned Tendon Stem Cell-Derived Extracellular Vesicles.

Remodeling of the immune microenvironment is a critical determinant of tissue regeneration, yet the molecular programs associated with the enhanced therapeutic activity of hypoxia-preconditioned extracellular vesicles remain incompletely defined. In this study, we investigated the regenerative and immunomodulatory effects of hypoxia-preconditioned tendon stem cell-derived extracellular vesicles (Hypo-EVs) and employed transcriptomic profiling to identify molecular signatures associated with their biological activity. The therapeutic effects of Hypo-EVs were evaluated using a rat patellar tendon defect model and lipopolysaccharide-stimulated RAW 264.7 macrophages. Histological analysis, immunostaining, biomechanical testing, and reverse transcription-quantitative polymerase chain reaction were performed to assess tendon healing and macrophage polarization, while RNA sequencing was conducted in macrophages treated with Hypo-EVs or normoxia-derived EVs, followed by Gene Set Enrichment Analysis, Gene Ontology, and Kyoto Encyclopaedia of Genes and Genomes pathway analyses. Hypo-EVs significantly alleviated local inflammatory responses, improved collagen organization and biomechanical properties of repaired tendons, and promoted macrophage polarization toward a reparative M2 phenotype both in&#xa0;vivo and in&#xa0;vitro. Consistent with these biological effects, transcriptomic profiling revealed extensive remodeling of inflammation-related gene expression programs, including significant suppression of NF-&#x3ba;B, TNF, IL-17, and cytokine-cytokine receptor interaction pathways. Integrative bioinformatic analyses identified an NF-&#x3ba;B-associated immune-regulatory signature that distinguished Hypo-EV-treated macrophages from those receiving normoxic EVs. Mechanistically, Hypo-EVs attenuated NF-&#x3ba;B activation, as evidenced by reduced phosphorylation of p65 and I&#x3ba;B&#x3b1;, whereas TNF-&#x3b1;-mediated NF-&#x3ba;B activation partially diminished their macrophage-repolarizing effects. Collectively, these findings demonstrate that hypoxic preconditioning enhances the immunomodulatory and regenerative functions of tendon stem cell-derived EVs. Transcriptomic analyses identified an NF-&#x3ba;B-associated immune-regulatory signature linked to the biological activity of Hypo-EVs, providing a molecular framework for understanding EV-mediated immune modulation and supporting the development of transcriptome-guided molecular signatures for regenerative therapies targeting tendon immune homeostasis.

Animals

New insights into diagnostic values and mechanisms of ferroptosis associated with immune infiltration in diabetic kidney disease.

The pathogenesis of diabetic kidney disease (DKD) is complex and closely related to ferroptosis and immune dysregulation, but the relevance is unclear. The present study investigates the potential mechanisms of ferroptosis-related genes (FRGs) in DKD and their relationship with the immune-inflammatory response. It searches for new diagnostic biomarkers to help diagnose and treat DKD. Four Gene Expression Omnibus (GEO) datasets, GSE30528, GSE30529 and GSE30122 as the test set, and GSE96804 for validation, were analyzed. FRGs were obtained from GeneCards, and 47 ferroptosis-related differentially expressed genes (FRDEGs) were identified by intersecting with DKD-related differentially expressed genes. Functional enrichment analyses, including Gene Ontology, Kyoto Encyclopedia of Genes and Genomes, Gene Set Enrichment Analysis and Gene Set Variation Analysis, revealed that these FRDEGs are primarily associated with ferroptosis, hypoxia response and immune inflammation. Subsequently, the weighted gene co-expression network analysis (WGCNA) was employed to expand the ferroptosis-related gene network, and intersection of the 47 FRDEGs with key WGCNA module genes yielded 10 key genes. Based on the 10 key genes, the least absolute shrinkage and selection operator and support vector machine algorithms identified three hub genes [chemokine ligand 5 (CCL5), forkhead box C1 (FOXC1) and lactotransferrin (LTF)] for DKD diagnosis. Receiver operating characteristic curves confirmed their diagnostic value, with FOXC1 and LTF validated in the independent dataset. Immune infiltration analysis via CIBERSORT revealed eight immune cell types with significantly different infiltration levels between the DKD and control group in the integrated GEO datasets. Notably, both LTF and CCL5 showed a significant positive correlation with gamma delta T cells (&#x3b3;&#x3b4;T). Quantitative PCR results confirmed differential expression of the three hub genes in the DKD group, with elevated expression observed in DKD mice following intervention with rosiglitazone and hyperoside.

bioinformatics analysis

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&#x2009;=&#x2009;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

Gene and pathway analysis of genome-wide genetic associations of bladder cancer.

BACKGROUND: Although genetic variants associated with bladder cancer (BCa) risk have been identified through hypothesis-driven and genome-wide association studies, a systematic understanding of BCa genetic susceptibility at the gene and pathway levels remains to be achieved. MATERIALS AND METHODS: In this 2-stage functional genomics study, we used 5 independent tools for genome-wide gene mapping and ranking based on BCa genome-wide association studies summary statistics, followed by a meta-analysis of gene-level significance p values, to obtain a consensus gene ranking in terms of association with BCa. Subsequently, we performed preranked gene-set enrichment analysis to identify the functional pathways involved in BCa genetic susceptibility. Joint analysis with gene-set enrichment analysis, based on somatic alteration frequency, was performed to explore the pathway-level relationships between genetic susceptibility and somatic alterations in BCa. RESULTS: Other than the well-known BCa genes (such as FGFR3, MYC, TERT, CCNE1, and TP63), we additionally prioritized a set of novel genes likely to be genetically implicated in BCa development, including SETD2, a possible tumor suppressor gene involved in chromatin remodeling. We further demonstrated convergence between genetic associations and somatic alterations at both the gene (eg, FGFR3 and TERT) and pathway levels (eg, cell cycle and chromatin modification), as well as functional ontologies specifically implicated in germline predisposition to BCa (eg, CD8/TCR signaling, immune checkpoints, and cytokine signaling). CONCLUSIONS: We identified several novel genes associated with BCa and demonstrated that genetic variants contribute to the development of BCa by affecting antitumor immunity, response to toxic exposure, and RNA and protein homeostasis and synergizing with somatic alterations in various cancer-related pathways.

Bladder cancer

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

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&#x2009;=&#x2009;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&#x2009;=&#x2009;9) or sevoflurane (n&#x2009;=&#x2009;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&#x2009;<&#x2009;0.0001), early-stage tumors (p&#x2009;<&#x2009;0.05) and higher survival (p&#x2009;=&#x2009;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

Dynamic Alterations in the Blood Transcriptome Characterize Drug Use Behavior and Co-Morbidities in Cocaine Use Disorder: A Preliminary Study.

Individuals with cocaine use disorder (CUD) who attempt abstinence experience craving and relapse that can benefit from multimodal treatment monitoring. Longitudinal studies linking behavioral manifestations in CUD to the blood transcriptome are not only limited but also computationally complex. Therefore, we developed an analytical pipeline to investigate the connection between drug use behaviors during abstinence and change in the blood transcriptome. We conducted a longitudinal study with CUD (n&#x2009;=&#x2009;12 subjects) and collected behavioral metrics and blood RNA-seq at baseline, 3, 6, and 9&#x2009;months. Our analytical pipeline of the high-dimensional data encompasses hierarchical k-means clustering to classify subjects to responder groups based on behavioral scores and abstinence duration, in silico cell deconvolution, differential analysis with correlated multivariate testing over time, gene set enrichment analysis, and gene co-expression with time splines and RNA-seq data. The pipeline captured dynamic changes in behavioral scores and abstinence duration in responder groups. Genes showing differential transcript-level expression were enriched in substance use and cardiovascular disease-associated genetic risk loci in responder groups. Lastly, time-dependent gene co-expression revealed dynamic changes related to immune processes, cell cycle, RNA-protein synthesis, and second messenger signaling for days of abstinence. This is a preliminary investigation, providing an innovative and scalable pipeline for blood-based longitudinal RNA-seq studies in CUD, potentially applicable to other substance use disorders. It outlines a data-driven approach for analyzing composite longitudinal drug use behavioral phenotypes with blood-based transcriptomics. We also demonstrate changes in drug use behaviors and the blood transcriptome during drug abstinence.

Humans

Integrating Genomic and Nongenomic Data to Stratify the Risk of Contralateral Breast Cancer After Radiation Therapy.

PURPOSE: Women treated with radiation therapy (RT) for breast cancer have an increased risk of developing radiation-associated contralateral breast cancer (CBC). Predicting CBC events is challenging because of the complex interplay of genomic, treatment, personal, and clinical factors. This study investigated computational methods that integrate genome-wide single-nucleotide polymorphisms and nongenomic data to develop a risk stratification model for developing CBC in women treated with RT for their first primary breast cancer. METHODS AND MATERIALS: This study used a subset of the population-based Women's Environmental Cancer and Radiation Epidemiology study that included 633 CBC cases and 1253 individually matched unilateral breast cancer controls who were treated with RT and had single-nucleotide polymorphism data available from a genome-wide association study. The study population was split into training, validation, and test sets for rigorous modeling and validation. Three data integration methods were compared in terms of their ability to stratify CBC risk: (1) naive integration; (2) sequential integration; and (3) sequential iterative integration. A biological analysis of the final model was performed using gene set enrichment analysis and protein-protein interaction analysis with gene annotation information informed by the model. RESULTS: The best-performing integration method was the sequential iterative integration equipped with the mixed-effect random forest algorithm. This approach achieved an area under the curve of 0.64 to stratify CBC risk in the test set, representing moderate predictive power. Calibration analysis showed good agreement between the lowest and highest risk bins stratified using sorted predicted values in the test set, resulting in an odds ratio of 3.27 for both predicted and observed CBC occurrence. Gene set enrichment analysis and protein-protein interaction analysis revealed that genes with high importance scores were associated with pathways relevant to lipid and fatty acid metabolism as well as breast cancer sensitivity to tamoxifen. CONCLUSIONS: The mixed-effect random forest approach demonstrated the potential for integrating high-dimensional genomic and low-dimensional nongenomic data to stratify CBC risk.

Humans

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

Identification of multicohort-based predictive signature for NMIBC recurrence reveals SDCBP as a novel oncogene in bladder cancer.

BACKGROUND: Despite surgical and intravesical chemotherapy interventions, non-muscle invasive bladder cancer (NMIBC) poses a high risk of recurrence, which significantly impacts patient survival. Traditional clinical characteristics alone are inadequate for accurately assessing the risk of NMIBC recurrence, necessitating the development of novel predictive tools. METHODS: We analyzed microarray data of NMIBC samples obtained from the ArrayExpress and GEO databases. LASSO regression was utilized to develop the predictive signature. We combined gene signature and clinicopathological factors to construct a clinical nomogram for estimating NMIBC recurrence in a local cohort. Finally. the biological functions and potential mechanisms of SDCBP in bladder cancer were investigated experimentally in vitro and in vivo. RESULTS: An 8-gene signature was developed, and its efficiency for predicting NMIBC recurrence was evaluated using Kaplan-Meier and time-dependent ROC curves in both training and validation datasets. Immunohistochemical testing revealed elevated levels of ACTN4 and SDCBP in recurrent NMIBC tissues. We integrated the two proteins with clinical factors to develop a nomogram model, which showed superior accuracy compared to individual parameters. Gene Set Variation Analysis and Gene Set Enrichment Analysis unveiled SDCBP exerted cancer-promoting biological processes, such as angiogenesis, EMT, metastasis and proliferation. Experimental procedures demonstrated that silencing SDCBP attenuated cell growth, glucose metabolism and extracellular acidification rate, accompanied by decreased expression of p-AKT, p-ERK1/2, LDHA and Vimentin. CONCLUSIONS: The established 8-gene signature holds promise as a tool for predicting NMIBC recurrence, while targeting SDCBP may represent a potential strategy for delaying disease relapse.

Urinary Bladder Neoplasms

Development and validation of a novel risk stratification signature derived from migrasome and tumor microenvironment-related genes for molecular subtyping and improving clinical outcomes in head and neck squamous cell carcinoma.

BACKGROUND: The tumor microenvironment (TME) and migrasomes released by tumor cells significantly influence carcinogenesis and immune evasion. However, our understanding of the prognostic and therapeutic implications of migrasome and tumor microenvironment-related genes (mtmRGs) in head and neck squamous cell carcinoma (HNSCC) remains limited. METHODS: We explored the relationship between mtmRGs and HNSCC prognosis by utilizing The Cancer Genome Atlas (TCGA) and the Gene Expression Omnibus (GEO) databases. Subsequently, we developed an innovative prognostic signature, and assessed its prognostic significance using the Kaplan-Meier method, time-dependent receiver operating characteristic (ROC), and Cox regression analyses. To explore the underlying mechanisms, we conducted gene set variation analysis (GSVA), gene set enrichment analysis (GESA), and immune infiltration analysis. A nomogram was developed to estimate the overall survival (OS) rates for HNSCC patients. Lastly, we chose P4HA1, which was part of the signature, for additional experimental validation in vitro and in vivo. RESULTS: The mtmRGs signature effectively classifies HNSCC patients into two distinct risk subgroups, with the high-risk cohort demonstrating significantly poorer OS. The risk score serves as an independent prognostic factor for HNSCC patients; those with lower risk scores are more likely to exhibit favorable responses to immunotherapy, particularly with CTLA4 inhibitors. Furthermore, a lower risk score is significantly correlated with the sensitivity of HNSCC patients to cyclophosphamide, gemcitabine, and axitinib. CONCLUSION: This study presents an innovative gene signature associated with mtmRGs, which may be utilized both for predicting survival and directing personalized chemotherapy and immunotherapy regiments for patients with HNSCC.

Humans

Differential Proteomic Profiling of Responders and Non-responders to Direct-Acting Antivirals Treatment in Chronic Hepatitis C Virus Infection.

Hepatitis C Virus (HCV), particularly genotype 3 (GT-3), is highly prevalent in India and is associated with faster progression to cirrhosis, hepatocellular carcinoma, and higher treatment failure rates. Although Direct-Acting Antivirals (DAAs) have revolutionized HCV therapy, 5-10% of patients fail to achieve sustained virological response (SVR). This proteomic study aimed to identify changes in the proteomic profile before and after treatment of both responders and non-responders to HCV treatment. Paired plasma samples from HCV GT-3 infected patients were collected before and 12 weeks after initiating DAAs treatment, along with healthy controls. Quantitative proteomic analysis was performed on the paired samples. Differentially expressed proteins (DEPs) were identified and subjected to functional analysis including gene set enrichment analysis (GSEA) and protein-protein interaction (PPI) network analysis. GSEA revealed enrichment in extracellular matrix organization and innate immune pathways. Expression patterns of candidate proteins selected based on fold change and false discovery rate (FDR) criteria were further evaluated in an independent cohort. Western blot confirmed key expression trends of candidate proteins. Proteins linked to extracellular matrix remodeling and angiogenesis showed differential expression patterns. Successful validation of these candidate proteins in large independent cohorts holds potential to predict therapeutic outcomes.

Humans