Abstract:Background: There is a growing public concern about diabetic kidney disease (DKD), which poses a severe threat to human health and life. It is important to discover noninvasive and sensitive immune-associated biomarkers that can be used to predict DKD development. ScRNA-seq and transcriptome sequencing were performed here to identify cell types and key genes associated with DKD. Methods: Here, this study conducted the analysis through five microarray datasets of DKD (GSE131882, GSE1009, GSE30528, GSE96804, and GSE104948) from gene expression omnibus (GEO). We performed single-cell RNA sequencing analysis (GSE131882) by using CellMarker and CellPhoneDB on public datasets to identify the specific cell types and cell-cell interaction networks related to DKD. DEGs were identified from four datasets (GSE1009, GSE30528, GSE96804, and GSE104948). The regulatory relationship between DKD-related characters and genes was evaluated by using WGCNA analysis. Gene Ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) datasets were applied to define the enrichment of each term. Subsequently, immune cell infiltration between DKD and the control group was identified by using the "pheatmap" package, and the connection Matrix between the core genes and immune cell or function was illuminated through the "corrplot" package. Furthermore, RcisTarget and GSEA were conducted on public datasets for the analysis of the regulation relationship of key genes and it revealed the correlation between 3 key genes and top the 20 genetic factors involved in DKD. Finally, the expression of key genes between patients with 35 DKD and 35 healthy controls were examined by ELISA, and the relationship between the development of DKD rate and hub gene plasma levels was assessed in a cohort of 35 DKD patients. In addition, we carried out immunohistochemistry and western blot to verify the expression of three key genes in the kidney tissue samples we obtained. Results: There were 8 cell types between DKD and the control group, and the number of connections between macrophages and other cells was higher than that of the other seven cell groups. We identified 356 different expression genes (DEGs) from the RNA-seq, which are enriched in urogenital system development, kidney development, platelet alpha granule, and glycosaminoglycan binding pathways. And WGCNA was conducted to construct 13 gene modules. The highest correlations module is related to the regulation of cell adhesion, positive regulation of locomotion, PI3K-Akt, gamma response, epithelial-mesenchymal transition, and E2F target signaling pathway. Then we overlapped the DEGs, WGCNA, and scRNA-seq, SLIT3, PDE1A and CFH were screened as the closely related genes to DKD. In addition, the findings of immunological infiltration revealed a remarkable positive link between T cells gamma delta, Macrophages M2, resting mast cells, and the three critical genes SLIT3, PDE1A, and CFH. Neutrophils were considerably negatively connected with the three key genes. Comparatively to healthy controls, DKD patients showed high levels of SLIT3, PDE1A, and CFH. Despite this, higher SLIT3, PDE1A, and CFH were associated with an end point rate based on a median follow-up of 2.6 years. And with the gradual deterioration of DKD, the expression of SLIT3, PDE1A, and CFH gradually increased. Conclusions: The 3 immune-associated genes could be used as diagnostic markers and therapeutic targets of DKD. Additionally, we found new pathogenic mechanisms associated with immune cells in DKD, which might lead to therapeutic targets against these cells.

Small RNA sequencing reveals snoRNAs and piRNA-019825 as novel players in diabetic kidney disease

The Role of Small Extracellular Vesicle Non-Coding RNAs in Kidney Diseases

Noncoding RNAs and diabetic kidney disease

Dysregulated Urinary Extracellular Vesicle Small RNAs in Diabetic Nephropathy: Implications for Diagnosis and Therapy

Single-cell RNA sequencing in diabetic kidney disease: a literature review

Selected miRNAs in urinary extracellular vesicles show promise for non-invasive diagnostics of diabetic kidney disease

Single‐cell RNA‐seq with spatial transcriptomics to create an atlas of human diabetic kidney disease

Noncoding RNAs as therapeutic targets in early stage diabetic kidney disease

Non-Coding RNAs as Biomarkers and Therapeutic Targets for Diabetic Kidney Disease

The role of long non‐coding RNAs in the development of diabetic kidney disease and the involved clinical application

Abstract P124: Human Small Artery MicroRNA Expression Profiles: Changes in Type 2 Diabetes and Associations with Endothelial Function

Identification of urinary extracellular vesicles differentially expressed RNAs in diabetic nephropathy via whole-transcriptome integrated analysis

Bioinformatics prediction and experimental verification of key biomarkers for diabetic kidney disease based on transcriptome sequencing in mice

RNA-sequencing reveals genome-wide long non-coding RNAs profiling associated with early development of diabetic nephropathy

Analysis of potential biomarkers for diabetic kidney disease based on single-cell RNA-sequencing integrated with a single-cell sequencing assay for transposase-accessible chromatin

Single-cell RNA and transcriptome sequencing profiles identify immune-associated key genes in the development of diabetic kidney disease

The advance of single cell transcriptome to study kidney immune cells in diabetic kidney disease

Identification of differentially expressed genes in diabetic kidney disease by RNA-Seq analysis of venous blood platelets

Single-Cell RNA Sequencing Profiles Identify Important Pathophysiologic Factors in the Progression of Diabetic Nephropathy

Kidney microRNA Expression Pattern in Type 2 Diabetic Nephropathy in BTBR Ob/Ob Mice

Long noncoding RNAs may impact podocytes and proximal tubule function through modulating miRNAs expression in Early Diabetic Kidney Disease of Type 2 Diabetes Mellitus patients