Identification and analysis of diverse cell death patterns in diabetic kidney disease using microarray-based transcriptome profiling and single-nucleus RNA sequencing

Abstract

Background: Diabetic kidney disease (DKD) is the most lethal complication of diabetes. Diverse programmed cell death (PCD) has emerged as a crucial disease phenotype that has the potential to serve as an indicator of renal function decline and can be used as a target for researching drugs for DKD.

Methods: Microarray-based transcriptome profiling and single-nucleus transcriptome sequencing (snRNA-seq) related to DKD were retrieved from the Gene Expression Omnibus (GEO) database. 13 PCD-related genes (including alkaliptosis, apoptosis, autophagy-dependent cell death, cuproptosis, disulfidptosis, entotic cell death, ferroptosis, lysosome-dependent cell death, necroptosis, netotic cell death, oxeiptosis, parthanatos, and pyroptosis) were obtained from various public databases and reviews. The gene set variation analysis (GSVA) analysis was used to explore the pathway activity of these 13 PCDs in DKD, and the pathway activity of these PCDs in different renal cells was studied based on DKD-related snRNA-seq data. To identify the core PCDs that play a significant role in DKD, we analyzed the relationships between different types of PCD and immune infiltration, fibrosis-related gene expression levels, glomerular filtration rate (GFR), and diagnostic efficiency in DKD. Using the Weighted Gene Co-expression Network Analysis (WGCNA) algorithm, we screened for core death genes among the core PCDs and constructed a cell death-related signature (CDS) risk score based on the Least Absolute Shrinkage and Selection Operator (LASSO). Finally, we validated the predictive performance of the CDS risk score in an independent validation set.

Results: We identified 4 core PCD pathways, namely entotic cell death, apoptosis, necroptosis, and pyroptosis in DKD, and further applied the WGCNA algorithm to screen 4 core death genes (CASP1, CYBB, PLA2G4A, and CTSS) and constructed a CDS risk score based on these genes. The CDS risk score demonstrated high diagnostic efficiency for DKD patients, and those with higher scores had higher levels of immune cell infiltration and poorer GFR.

Conclusion: Our study sheds light on the fact that multiple PCDs contribute to the progression of DKD, highlighting potential therapeutic targets for treating this disease.

Keywords: Diabetic kidney disease; Programmed cell death; Transcriptome profiling; snRNA-seq.