Multi-omics Analyses Identify AKR1A1 as a Biomarker for Diabetic Kidney Disease

Diabetic kidney disease (DKD) is the leading cause of end-stage kidney disease. Because many genes associate with DKD, multiomics approaches were used to narrow the list of functional genes, gene products, and related pathways providing insights into the pathophysiological mechanisms of DKD. The Kidney Precision Medicine Project human kidney single-cell RNA-sequencing (scRNA-seq) data set and Mendeley Data on human kidney cortex biopsy proteomics were used. The R package Seurat was used to analyze scRNA-seq data and data from a subset of proximal tubule cells. PathfindR was applied for pathway analysis in cell type-specific differentially expressed genes and the R limma package was used to analyze differential protein expression in kidney cortex. A total of 790 differentially expressed genes were identified in proximal tubule cells, including 530 upregulated and 260 downregulated transcripts. Compared with differentially expressed proteins, 24 genes or proteins were in common. An integrated analysis combining protein quantitative trait loci, genome-wide association study hits (namely, estimated glomerular filtration rate), and a plasma metabolomics analysis was performed using baseline metabolites predictive of DKD progression in our longitudinal Diabetes Heart Study samples. The aldo-keto reductase family 1 member A1 gene (AKR1A1) was revealed as a potential molecular hub for DKD cellular dysfunction in several cross-linked pathways featured by deficiency of this enzyme.