GPX3 and GSTT1 as biomarkers related to oxidative stress during renal ischemia reperfusion injuries and their relationship with immune infiltration.

Background: Renal ischemia reperfusion injuries (IRIs) are very common in clinical diagnoses and treatments, which are a common cause of impaired renal functions, worsening pathological damage, affecting disease progression and hindering recovery. Renal IRIs are an inflammatory disease mediated by the adaptive and innate immune system. There is a complex interaction between oxidative stress and immune cell infiltration. Therefore, we aimed to determine biomarkers associated with oxidative stress during renal IRIs and their relationship with immune cell infiltration. Method: A differential gene expression analysis was made based on the GSE148420 dataset from the NCBI Gene Expression Comprehensive Database (GEO) combined with 92 oxidative-stress (OS)-related genes identified in the Molecular Signatures Database. Then we identified differentially-expressed genes (DEOSGs) associated with oxidative stress, which were used for gene ontology (GO) and a Kyoto Encyclopedia of Genomes (KEGG) enrichment analysis. At the same time, we used PPI protein interaction networks and Lasso regression analysis to identify key genes, which were verified by the validation sets GSE58438 and GSE71647, as well as Western Blot detection on rat renal IRI models. At the same time, PAS staining, HE staining and immunohistochemistry were used to detect tissue damage and expression of markers related to oxidative stress during renal ischemia-reperfusion. Single-gene enrichment analysis (GSEA) was used to further clarify the underlying biological functions of key genes. Cibersort was used to analyze the immune cell infiltration during renal IRI and the correlation of key genes with immune cells. At the same time, we constructed a network of transcription-factor (TF)-Hub genes and miRNA-Hub genes. DGIDB was used to predict drugs and molecular compounds that might interact with the Hub genes. Results: Compared with the control group, a total of 5456 differential genes (DEGs) were measured in the renal IRI group, 2486 of which were upregulated and 2970 were down-regulated. Among them, we found 30 DEGs (DEOSGs) associated with oxidative stress. The results of GO and KEGG enrichment analysis showed that these DEOSGs were mainly enriched in glutathione metabolism, the response to oxidative stress stimulation, the regulation of T cell activation and apoptosis signaling pathways. Through a protein interaction network (PPI) and a LASSO regression analysis, a total of two Hub genes were identified, namely GPX3 and GSTT1, which were validated through external validation sets and animal experiments. Through pathological methods, we found that the pathological damage of renal tissue and the expression of oxidative stress markers increased after renal ischemia-reperfusion. The results of GSEA showed that the Hub genes were related to oxidative stress pathways, apoptosis signaling pathways and immune-response-related signaling pathways. An immunoinfiltration correlation analysis showed that genes GPX3 and GSTT1 were significantly positively correlated with plasma cells and macrophage M0, while were negatively correlated with monocytes and macrophages M1 and M2. Using the Strust, Starbase and DGIDB database, we predicted that 81 transcription factors, 49 miRNAs and 13 drug or molecular compounds might interact with the Hub genes. Conclusion: Through a comprehensive analysis of gene expression, our findings may provide new potential biomarkers for the pathogenesis of renal IRIs and a reliable basis for its early diagnosis as well as treatment.