Exosome-related gene identification and diagnostic model construction in hepatic ischemia-reperfusion injury

Hepatic ischemia-reperfusion injury (HIRI) may cause severe hepatic impairment, acute hepatic insufficiency, and multiorgan system collapse. Exosomes can alleviate HIRI. Therefore, this study explored the role of exosomal-related genes (ERGs) in HIRI using bioinformatics to determine the underlying molecular mechanisms and novel diagnostic markers for HIRI. We merged the GSE12720, GSE14951, and GSE15480 datasets obtained from the Gene Expression Omnibus (GEO) database into a combined gene dataset (CGD). CGD was used to identify differentially expressed genes (DEGs) based on a comparison of the HIRI and healthy control cohorts. The impact of these DEGs on HIRI was assessed through gene set enrichment analysis (GSEA) and gene set variation analysis (GSVA). ERGs were retrieved from the GeneCards database and prior studies, and overlapped with the identified DEGs to yield the set of exosome-related differentially expressed genes (ERDEGs). Functional annotations and enrichment pathways of these genes were determined using Gene Ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) analyses. Diagnostic models for HIRI were developed using least absolute shrinkage and selection operator (LASSO) regression and support vector machine (SVM) algorithms. Key genes with diagnostic value were identified from the overlap, and single-sample gene-set enrichment analysis (ssGSEA) was conducted to evaluate the immune infiltration characteristics. A molecular regulatory interaction network was established using Cytoscape software to elucidate the intricate regulatory mechanisms of key genes in HIRI. Finally, exosome score (Es) was obtained using ssGSEA and the HIRI group was divided into the Es_High and Es_Low groups based on the median Es. Gene expression was analyzed to understand the impact of all genes in the CGD on HIRI. Finally, the relative expression levels of the five key genes in the hypoxia-reoxygenation (H/R) model were determined using quantitative real-time PCR (qRT-PCR). A total of 3810 DEGs were identified through differential expression analysis of the CGD, and 61 of these ERDEGs were screened. Based on GO and KEGG enrichment analyses, the ERDEGs were mainly enriched in wound healing, MAPK, protein kinase B signaling, and other pathways. GSEA and GSVA revealed that these genes were mainly enriched in the TP53, MAPK, TGF[Formula: see text], JAK-STAT, MAPK, and NFKB pathways. Five key genes (ANXA1, HNRNPA2B1, ICAM1, PTEN, and THBS1) with diagnostic value were screened using the LASSO regression and SVM algorithms and their molecular interaction network was established using Cytoscape software. Based on ssGSEA, substantial variations were found in the expression of 18 immune cell types among the groups (p < 0.05). Finally, the Es of each HIRI patient was calculated. ERDEGs in the Es_High and Es_Low groups were enriched in the IL18, TP53, MAPK, TGF[Formula: see text], and JAK-STAT pathways. The differential expression of these five key genes in the H/R model was verified using qRT-PCR. Herein, five key genes were identified as potential diagnostic markers. Moreover, the potential impact of these genes on pathways and the regulatory mechanisms of their interaction network in HIRI were revealed. Altogether, our findings may serve as a theoretical foundation for enhancing clinical diagnosis and elucidating underlying pathogeneses.