Abstract
Rejection seriously affects the success of kidney transplantations. However, the molecular mechanisms underlying this rejection remain unclear. The GSE21374 and GSE36059 datasets were downloaded from the Gene Expression Omnibus (GEO) database. Next, the Cell-type Identification by Estimating Relative Subsets of RNA Transcripts (CIBERSORT) algorithm was used to infer the proportions of 22 immune cells. Moreover, infiltrating immune cell-related genes were identified using weighted gene co-expression network analysis (WGCNA), and enrichment analysis was conducted to observe their biological functions. Extreme Gradient Boosting (XGBoost) and Least Absolute Shrinkage and Selection Operator (LASSO) logistic regression algorithms were used to screen hub genes. Quantitative real-time PCR was conducted to verify the number of immune cells and hub gene expression levels. The rejection and non-rejection groups showed significantly different distributions (P < 0.05) of eight immune cells (B cell memory, Plasma cells, mast cells, follicular helper T cells, T CD8 cells, Macrophages M1, T Cells CD4 memory activated, and gamma delta T cells). Subsequently, CD8A, CRTAM, GBP2, WARS, and VAMP5 were screened as hub genes using the XGBoost and LASSO algorithms and could be used as diagnostic biomarkers. Finally, differential analysis and quantitative real-time PCR suggested that CD8A, CRTAM, GBP2, WARS, and VAMP5 were upregulated in rejection samples compared to non-rejection samples. The present study identified five key infiltrating immune cell-related genes (CD8A, CRTAM, GBP2,WARS, and VAMP5) involved in kidney transplant rejection, which may explain the molecular mechanism of rejection in kidney transplantation development.