Correlation Between Serum 25-Hydroxyvitamin D Levels in Albuminuria Progression of Diabetic Kidney Disease and Underlying Mechanisms By Bioinformatics Analysis

Abstract

AimThis study aimed to assess the correlation between serum concentration of 25-hydroxyvitamin D and albuminuria progression of diabetic kidney disease (DKD) and to use bioinformatics methods to determine the potential mechanism in the pathological process of advanced DKD.MethodsA total of 178 type 1 diabetes mellitus (T1DM) patients with microalbuminuria complications who were hospitalized at least twice (with an interval > 24 months) in the Department of Endocrinology of The First Affiliated Hospital of USTC were included in this study. According to the urinary albumin creatinine ratio (UACR), we classified DKD stages as follows: microalbuminuria (UACR, 30-300 mg/g), and macroalbuminuria (UACR, >300 mg/g). We divided the patients into DKD progression (N=44) and stable group (N=134) on account of urinary albumin-to-creatinine ratio (UACR) by at least two randomized measurements. Stable group was defined as UACR between 30 and 300 mg/g, whereas progression group was defined as UACR >300 mg/g at the end of follow-up. Data were obtained from participants’ medical records, and the 25-hydroxyvitamin D level was categorized into three groups as follows: G1 (N=45), <10 ng/mL; G2 (N=80), 10-20 ng/ml; and G3 (N=53), ≥20 ng/mL. The Nephroseq database (http://v5.nephroseq.org) was used to identify VDR expression in diabetic nephropathy. The dataset GSE142025 from GEO (http://www.ncbi.nlm.nih.gov/geo) was downloaded. After stratification by the median-centered log2 VDR expression value, the 21 advanced DKD samples were divided into two groups (low VDR expression group and high VDR expression group). Gene set enrichment analysis (GSEA) (http://software.broadinstitute.org/gsea/index.jsp). Differentially expressed genes (DEGs) were screened by the limma package (adjusted p < 0.05, |logFC| > 1). The Gene Ontology (GO; http://www.geneontology.org/) database and pathway analysis within the Kyoto Encyclopedia of Genes and Genomes (KEGG; https://www.kegg.jp/) were performed using the R package ClusterProfile. The CIBERSORT (Cell type Identification By Estimating Relative Subsets Of known RNA Transcripts) algorithm was utilized for calculating the infiltrated immune cells in advanced kidney tissues.Results1) A multivariate Cox regression analysis revealed that DR (diabetic retinopathy), eGFR (estimated glomerular filtration rate), and 25-hydroxyvitamin D were significant independent predictors of DKD progression (HR: 2.57, 95% CI: 1.44.4.24, p=0.007; HR: 2.13, 95% CI: 1.58.3.79, p = 0.011; HR: 0.732, 95% CI: 0.232–0.816, p = 0.023, respectively). 2) Kaplan–Meier survival curves of DKD progression by serum 25-hydroxyvitamin D stratification showed that the G2 and G3 groups were significantly different when compared with the G1 group (log-rank χ2 = 14.69, p <0.001; χ2 = 28.26, p <0.001, respectively). 3) There was a weak negative correlation between 25-hydroxyvitamin D level and UACR at baseline,and the overall mean rate of change in eGFR was 1.121 ± 0.19 ml/min/1.73 m2/year. Neither crude nor adjusted rate of decline in eGFR was significantly different among patients classified according to baseline serum 25-hydroxyvitamin D levels (all p<0.05). 4) The high expression of VDR group was most positively correlated with enriched gene sets like reactome innate immune system and reactome G alpha I signaling events when compared with the low expression of VDR group. 5) The CIBERSORT algorithm showed decreased M2 macrophage infiltration in advanced kidneys in comparison to low VDR expression and high VDR expression.ConclusionThis study concluded that low 25-hydroxyvitamin D levels can predict an increased risk of DKD albuminuria progression and eGFR decline. Decreased M2 macrophage infiltration may be a potential mechanism involved in this pathogenesis.

Identification of a copper metabolism‐related gene signature for predicting prognosis and immune response in glioma

“…The “DESeq2” package was used to obtain the DEGs between the high‐ and low‐risk groups in the TCGA set. Then, GSEA 32 was performed using the “clusterProfiler” package. 33 …”

Identification of key biomarkers of the glomerulus in focal segmental glomerulosclerosis and their relationship with immune cell infiltration based on WGCNA and the LASSO algorithm

“…The volcano map of DEGs was generated using the ‘ggplot2’ package of R software. Data standardization was performed, and a p value <0.05 and log fold change > 1 were used as thresholds [ 11 ].…”

Identification of diagnostic markers related to oxidative stress and inflammatory response in diabetic kidney disease by machine learning algorithms: Evidence from human transcriptomic data and mouse experiments

“…For example, Ma used a bioinformatic approach to analyze gene expression profiles and underlying functional networks in cardiac tissue from patients with dilated cardiomyopathy ( 13 ). Huang analyzed the correlation of serum 25-hydroxyvitamin D levels in the progression of proteinuria in DKD and its underlying mechanisms ( 14 ). Yang explored seven immune-related genes that can predict the progression of atherosclerotic plaque based on machine learning ( 15 ).…”