Thank you for visiting nature.com. You are using a browser version with limited support for CSS. To obtain the best experience, we recommend you use a more up to date browser (or turn off compatibility mode in Internet Explorer). In the meantime, to ensure continued support, we are displaying the site without styles and JavaScript.
Advertisement
Scientific Reports volume 12, Article number: 21411 (2022)
2
Metrics details
In view of the alarming increase in the burden of diabetes mellitus (DM) today, a rising number of patients with diabetic kidney disease (DKD) is forecasted. Current DKD predictive models often lack reliable biomarkers and perform poorly. In this regard, serum myoglobin (Mb) identified by machine learning (ML) may become a potential DKD indicator. We aimed to elucidate the significance of serum Mb in the pathogenesis of DKD. Electronic health record data from a total of 728 hospitalized patients with DM (286 DKD vs. 442 non-DKD) were used. We developed DKD ML models incorporating serum Mb and metabolic syndrome (MetS) components (insulin resistance and β-cell function, glucose, lipid) while using SHapley Additive exPlanation (SHAP) to interpret features. Restricted cubic spline (RCS) models were applied to evaluate the relationship between serum Mb and DKD. Serum Mb-mediated renal function impairment induced by MetS components was verified by causal mediation effect analysis. The area under the receiver operating characteristic curve of the DKD machine learning models incorporating serum Mb and MetS components reached 0.85. Feature importance analysis and SHAP showed that serum Mb and MetS components were important features. Further RCS models of DKD showed that the odds ratio was greater than 1 when serum Mb was > 80. Serum Mb showed a significant indirect effect in renal function impairment when using MetS components such as HOMA-IR, HGI and HDL-C/TC as a reason. Moderately elevated serum Mb is associated with the risk of DKD. Serum Mb may mediate MetS component-caused renal function impairment.
Diabetic kidney disease (DKD) is a common microvascular complication of diabetes1,2, and 50% of patients with type 2 diabetic mellitus (T2DM) will eventually develop DKD3. Novel biomarker discovery and accurate prediction models will aid DKD intervention and reduce disease burden4,5. However, currently available predictive models often lack reliable biomarkers and perform poorly6,7,8,9.
Components of metabolic syndrome (MetS), a cluster of risk indicators reflecting conditions such as insulin resistance (IR), dyslipidemia, and comorbidities (i.e. prothrombotic and proinflammatory states)10, have been widely shown to be associated with cardiovascular disease events11,12. MetS components are also used for DKD risk prediction, but their predictive power is not ideal. For example, based on the ADVANCE study (2012) of 18,517 T2DM patients, the COX proportional hazards regression approach indicated that baseline blood glucose was not a good predictor of DKD, whereas baseline estimated glomerular filtration rate (eGFR) and urinary albumin-creatinine ratio (ACR) dominated risk prediction13.
To identify any potential DKD-related indicators, we performed machine learning (ML)-based data analysis to evaluate all relevant indicators in 51,866 hospitalized patients with T2DM. Interestingly, we found serum myoglobin (Mb), a 17 kDa monomeric hemoglobin14 that present in human myocardial and skeletal muscle tissues to facilitate the diffusion of O2 from the blood to muscle mitochondria15, as an important feature when excluding renal function indicators (Fig. S1). Since this is only a preliminary analysis finding and the influences of missing values and related diseases are not considered, further in-depth analysis is needed. From a molecular mechanism perspective, subclinical elevations in serum Mb have recently been observed in DM and chronic kidney disease (CKD)16, which may be explained by diabetic muscle damage and abnormal glucose and lipid metabolism, suggesting that MetS components may be contributing to elevated serum Mb levels17,18,19,20. Further, tissue damage caused by Mb has been confirmed at the molecular level. For example, studies have shown that Mb interacts with glucose to form end-of-glycosylation end products (AGEs), which ultimately lead to intramolecular and/or intermolecular cross-linking of long-lived proteins, causing chronic inflammation21,22,23. However, the importance of serum Mb in DKD has not been evaluated, nor has serum Mb been investigated as a risk factor in models of DKD. Moreover, the pathological role of serum Mb as a mediator in MetS component-induced renal function impairment has not been reported.
Therefore, to clarify the possible pathophysiological significance of serum Mb in DKD, we used a cross-sectional study design and obtained a total of 51,866 electronic health records (EHRs) of hospitalized patients with T2DM in the Third Xiangya Hospital of Central South University. We established DKD ML models incorporating serum Mb and simultaneously evaluated the models’ performance and feature importance to verify the importance of serum Mb in DKD. Logistic regression-based restricted cubic spline (RCS) models were applied to further assess the relationship between serum Mb and DKD. We explored the causal relationship between MetS components, serum Mb, and DKD using causal mediation analysis, aiming to test the hypothesis that serum Mb mediates the induction of renal function impairment by MetS components.
After obtaining the data access authorization from the Information Network Center of the Third Xiangya Hospital, we retrospectively collected the EHRs of 51,866 hospitalized patients from the Hospital Information System (HIS). The Third Xiangya Hospital of Central South University is located in Changsha City, Hunan Province. It is a grade A tertiary general hospital serving the middle reaches of the Yangtze River. This study was approved by the IRB of the Third Xiangya Hospital of Central South University (No: 21112). All data had been completely anonymized prior to our access to the HIS. Because the patients’ personally identifiable information had been effectively protected, the IRB of the Third Xiangya Hospital of Central South University waived the requirement for informed consent. The present study was conducted in accordance with the guidelines of the Declaration of Helsinki.
The EHRs of 51,866 hospitalized patients with T2DM in the Third Xiangya Hospital from January 1, 2013 to December 31, 2020 were collected. Exclusion criteria: (1) under 18 years of age (n = 30), (2) diseases affecting Mb (n = 7457), (3) infectious diseases (n = 9258), (4) organ dysfunction (n = 4169), (5) special body status (n = 4810), and (6) missing urine microalbumin (UMALB), eGFR and/or serum Mb data (n = 25,387). Finally, 728 patients (DKD = 286, non-DKD = 442) were eligible for the study. Eligible patients were divided into training set (60%, earlier than March 9, 2016) and testing set (40%, later than March 9, 2016), including 437 (DKD = 158, non-DKD = 279) and 291 patients (DKD = 128, non-DKD = 163), respectively, based on hospitalization date. The inclusion and exclusion criteria of the study population are presented in Fig. 1, and the demographic and clinical features are available in Table 1.
Inclusion and exclusion criteria of study population. Electronic health records (EHRs) of 51,866 hospitalized patients with type 2 diabetic mellitus (T2DM) were collected. Patients meeting the following exclusion criteria were removed: (1) under 18 years of age (n = 30), (2) diseases affecting serum myoglobin (Mb) (n = 7457), (3) infectious diseases (n = 9258), (4) organ dysfunction (n = 4169), (5) special body status (n = 4810), and (6) missing urine microalbumin (UMALB), estimated glomerular filtration rate (eGFR) and/or serum Mb data (n = 25,387). Finally, 728 patients were eligible for our study. The training set and testing set sample sizes were 437 and 291, respectively.
From the EHRs of T2DM, we extracted a total of 157 features in the following 7 domains: demographic information, hospitalization-related information, general information, laboratory tests, diagnosis, medication, and others. The data domain categories and descriptions are listed in Table 2. The details of feature calculation are shown in Table S1.
Diagnoses were organized using the ICD -9 and ICD-10 hierarchies. These diseases were diagnosed by experienced doctors who were licensed to practice medicine. T2DM was diagnosed according to the World Health Organization (WHO) criteria (1999)24. DKD was diagnosed by the NKF/KDOQI (2007) classification25. Serum Mb was tested by Mb-LATEX CN from Denka Seiken (Reference range 0–70 ng/mL).
The components of MetS have the following three sets of features: IR and β-cell function were evaluated by fasting C-peptide, 120 min C-peptide, Homeostatic Model Assessment of Insulin Resistance26 (HOMA-IR), Gutt index27, triglyceride-glucose (TyG) index28, fasting insulin to glucose ratio29 (I/G 0 min), 120 min insulin to glucose ratio29 (I/G 120 min), insulinogenic index30 (IGI), IGI to HOMA-IR ratio31 (IGI/HOMA-IR), and Homeostasis Model of Assessment for beta-cell function32 (HOMA-BETA). Glucose metabolism indicators included fasting glucose, 120 min glucose, glycated hemoglobin A1c (GHBA1C), and hemoglobin glycation index33 (HGI). Lipid metabolism was assessed by total cholesterol (TC), triglycerides (TG), low-density lipoprotein cholesterol (LDL-C), high-density lipoprotein cholesterol (HDL-C), HDL-C/TC, and TG/HDL-C ratio34.
Two popular ML methods, extreme gradient boosting (XGBoost) and random forest (RF), were used to build DKD classification models based on EHR data. RF, first proposed by Breiman in 200135, is an ensemble of unpruned classification or regression trees created using bootstrap samples of training data and random feature selection in tree induction. XGBoost is an efficient and scalable implementation of the gradient boosting framework36. It is used to develop models in a sequential stagewise fashion like other boosting methods and generalize them by allowing optimization of arbitrary differentiable loss functions. RF and XGBoost were applied using the randomForest (4.6–14)37 and xgboost (1.3.2.1)38 R packages, respectively. All ML methods used grid search and fivefold cross-validation for parameter adjustment. The final hyperparameters for the DKD models are listed in Table S2.
The evaluation indicators of the confusion matrix, including sensitivity (SE), specificity (SP), positive predictive value (PPV), negative predictive value (NPV), recall, F1, and balanced accuracy (BA), were used to analyze the discrimination ability of DKD ML models. The probability scores output by the models were used to calculate the area under the receiver operating characteristic (ROC) curve (AUC). These statistical parameters are defined as follows:
where true positive (TP) is the number of correctly classified objects, true negative (TN) is the number of correctly misclassified objects, false positive (FP) is the number of incorrectly classified objects, and false negative (FN) is the number of incorrectly misclassified objects.
Feature importance of RF was estimated by (1) mean decrease in prediction accuracy without the feature in the model and (2) mean decrease in the Gini index, a measure of impurity of the dataset (i.e. risk of misclassifying data), by including the feature. For both parameters, the higher the score, the more important the feature.
We evaluated feature contributions to model predictions using an interpretability method called SHapley Additive exPlanation (SHAP). SHAP is based on the Shapley value, a concept introduced in the 1950s to measure each player’s contribution to a collaborative game39. In 2017, Lundberg and Lee proposed a broadly applicable SHAP method to interpret a variety of models, including regression and classification40,41. For binary classification in XGBoost, the output of the model is the log odds ratio, that is, the SHAP values sum up to the log loss of the model for each sample. The function is defined as follows
where (z) is the log odds ratio, ({x}^{^{prime}}) is a simplified binary feature vector containing values 1 or 0 for observed and missing features, M is the number of simplified input features, and ({varphi }_{i}) is the attribution value of the ith molecular descriptor or the SHAP value. In probability theory, (z) is the corresponding value of probability p after a certain transformation, and the transformation function is defined as follows
where (p) varies monotonically in the range (0, 1) and (z) varies from infinitesimal to infinity. A feature is predicted to be a risk factor if its (z)-value is greater than 0. The larger the z-value, the more likely it is a risk factor rather than a protective factor. Therefore, the SHAP method can help us understand the model more intuitively. The SHAPforxgboost (0.1.1) and fastshap (0.0.5) R packages were applied in the analysis.
We used DCA to assess the clinical utility of ML models. A decision curve is generated by plotting net benefit against a range of threshold probabilities. Decision curves include benefits and harms on the same scale, so they can be directly compared, supporting comparisons between models. The net benefit is calculated as
where N is the total sample size and ({P}_{t}) represents the threshold probability. Two extreme strategies, classifying for all and classifying for none, were added for reference. The R package rmda (1.6) was used.
The association between serum Mb and DKD was evaluated on a continuous scale with RCS curves based on logistic regression models. To balance best fit and overfitting, the number of knots was chosen empirically to be 4. The serum Mb associated with the lowest risk of DKD was the concentration with the lowest odds ratio on the spline curve. Odds ratios were adjusted for gender, age, BMI, hyperlipidemia, hypertension, age-adjusted Charlson comorbidity index (ACCI), and hospitalization date.
Causal mediation analysis identifies potential pathways to explain the observed association between a given exposure and a particular outcome42, examining how a third intermediate variable, the mediator, is related to the observed outcome relationship. By using a counterfactual framework, causal mediation analysis allows us to decompose the total effect (TE) into direct effects (DE) and indirect effects (IE). In our hypothetical model, we defined TE as the effect of MetS components on renal function impairment, DE as the effect of MetS components on renal function impairment after removing the individual effect of Mb on renal function impairment, and IE as the effect of MetS components on renal function impairment specially associated with Mb. A directed acyclic graph (DAG) illustrating mediation process is shown in Fig. 5A. Through the causal mediation analysis, we estimated the mediation effect of serum Mb based on logistic regression in all eligible patients, DKD or non-DKD subgroups. According to previous DKD investigators43,44, we adjusted and standardized the MetS components, Mb, UMALB, and eGFR for 7 variables: gender, age, BMI, hyperlipidemia, hypertension, ACCI, and hospitalization date. The proportion of serum Mb-mediated renal function impairment due to the MetS components was calculated as
Causal mediation analysis was performed using the R package mediation (4.5.0).
R software (version 4.0.3) was used as a statistical analysis tool in this study. Data are presented as median [interquartile range] for continuous variables or number (%) for categorical variables. Differences between groups were based on the Wilcoxon test for continuous variables and the chi-square test for categorical variables. Relationships between MetS components, serum Mb, and DKD were assessed with the Spearman’s correlation coefficient (rs), controlling for 7 variables: gender, age, BMI, hyperlipidemia, hypertension, ACCI, and admission date. All statistical significance was defined as P < 0.1. Missingness was addressed by setting missing categorical data to 0, whereas for numerical data, imputation was completed using the predictive mean matching (PMM) strategy.
Observer bias may arise from the data on patients hospitalized in a grade A tertiary general hospital, as these patients may have more severe than average DM or DKD. Besides, selection bias caused by exclusion of UMALB missing individuals brings DKD rate elevation but comparable cardiovascular risks, leading to limited representation (Table S3). However, the significance and possible mediator role of serum Mb in DKD, as a common rule, was mainly focused in our study, and we believe that the bias did not greatly affect our conclusion.
The clinical characteristics of the participants are presented in Table 1. There were no significant between-group differences in age, sex, and BMI. However, the DKD group showed significantly more diabetes complications such as diabetic retinopathy and diabetic neuropathy. Importantly, serum Mb was significantly elevated in the DKD group, consistent with previous studies45.
To verify and evaluate the significance of serum Mb in DKD, we built DKD ML models incorporating serum Mb, taking advantage of its freedom from collinearity in processing high-dimensional EHR data on disease phenotypes and clinical outcomes46,47,48, aiming to discover important features overlooked by traditional methods. We first accessed the model performance to ensure that the following feature importance analysis and model interpretation are accurate and reliable.
After removing features with a missing ratio of more than 40% or zero variance, a total of 88 important features (including renal function) were available for DKD modeling to obtain high-accuracy DKD classification models. Two popular ML algorithms, XGBoost and RF algorithms, have been applied, named XGBoost-88 and RF-88, respectively. The model performances in the testing set are summarized in Table 3 and the ROC curves are shown in Fig. 2A. Overall, both approaches performed satisfactorily (e.g. all AUCs > 0.85). Clearly, the best classifier based on the RF algorithm was able to achieve AUC values of 0.884, BA values of 0.798, and PPV values of 0.854 in the testing set.
Receiver operating characteristic (ROC) curves and decision curve analysis (DCA) results in the testing set. (A,B) ROC curves of XGBoost and RF models based on 88 features and 83 features, respectively. (C,D) DCA for classifying DKD in T2DM patients by machine learning (ML) models based on 88 features and 83 features, respectively. The net benefit and 95% confidence interval for the ML models are shown with blue and red dashed lines. Two extreme strategies, classifying for all and classifying for none, are represented by solid grey and solid black lines.
We also applied XGBoost and RF algorithms to renal function exclusion dataset (contained 83 features and removed UMALB, serum creatinine (Scr), eGFR, blood urea nitrogen (BUN), and serum uric acid (SUA)), aiming to find the features hidden by renal function features. XGBoost-83 and RF-83 were named respectively. Model performances in the testing set are summarized in Table 3 and ROC curves are shown in Fig. 2B. On the whole, the performances of both approaches were relatively satisfactory (e.g. all AUCs > 0.80). Compared with the model including renal function features, BA decreased by about 0.1, SP did not change much, and SE decreased by 0.2 in the renal-function-excluded-model, indicating that renal function indicators are crucial for identifying DKD, while other features are important for ruling out DKD.SE sensitivity, SP specificity, PPV positive predictive value, NPV negative predictive value, BA balanced accuracy, AUC area under curve.
To assess the clinical utility of several ML models we built, DCA was applied to these ML models (Fig. 2C,D). For both XGBoost/RF-88 and XGBoost/RF-83, the net benefits of two models were similar across the range of threshold probabilities. When the threshold probability was between 0.4 and 0.6, both models outperformed the reference strategies. These results of DCA show the consistency and efficiency of our models and are consistent with our model performance data (Table 3).
Next, to identify important features hidden from indicators of renal function and attest the importance of serum Mb for DKD, feature importance analysis and model interpretation were applied, the reliability of which has been proved by model performance and DCA.
The top 30 important variables for RF-88 and RF-83 are presented in Fig. 3A,B, and two estimators, mean decrease accuracy and mean decrease Gini, have been applied to measure feature importance. In RF-88, UMALB, eGFR, and Scr were listed as top-ranked features. It’s worth noting that serum Mb consistently ranked in the front. Meanwhile, MetS components such as I/G 120 min, I/G 0 min, HOMA-BETA, HOMA-IR, TyG index, Gutt index, 120 min C-peptide (reflecting IR and β-cell function), HGI, fasting glucose (reflecting glucose metabolism), and LDL-C, TG/HDL ratio (reflecting lipid metabolism) were on the list. Intriguingly, in RF-83, serum Mb emerged as the most important feature under all estimators. Meanwhile, MetS components such as 120 min C-peptide, fasting C-peptide, IGI, IGI/HOMA-IR, HOMA-BETA, HOMA-IR, I/G 0 min, TyG index (reflecting IR and β-cell function), fasting glucose (reflecting glucose metabolism), LDL-C, HDL-C (reflecting lipid metabolism) were relatively important features.
Feature importance analysis and model interpretation results. (A,B) Feature importance was estimated by mean decrease accuracy (left) and mean decreased Gini index (right) based on RF-88 and RF-83, respectively. Serum Mb and metabolic syndrome (MetS) components are indicated by red and orange dots, respectively. (C) Global and local importance for top 30 important features detecting by SHapley Additive exPlanation (SHAP).
To overcome the black-box nature of XGBoost and RF algorithms and better understand the prediction models, model interpretation based on XGBoost-88 correct classifications was carried out using the SHAP approach. The top 30 important features identified by SHAP are shown in Fig. 3C. We found that the renal function features UMALB, Scr, and eGFR remained dominant. Serum Mb (ranked 13th) was recognized as an important feature. MetS components, including HOMA-IR, fasting C-peptide, 120 min C-peptide, IGI/HOMA-IR, HOMA-BETA, I/G 0 min, Gutt index (reflecting IR and β-cell function), 120 min glucose, fasting glucose, HGI (reflecting glucose metabolism), and HDL-C (reflecting lipid metabolism) were important features. These results are consistent with the feature importance analysis described above. We noted that higher serum Mb, worsening renal functions (i.e. higher UMALB and Scr and lower eGFR), more diabetes comorbidities (i.e. diabetic neuropathy and diabetic retinopathy), worsening IR and β-cell function (higher IGI/HOMA-IR, AC, HOMA-BETA, and I/G 0 min), and worse glucose control (e.g. higher HGI and AGLU) were associated with DKD.
Our results suggest that in addition to renal function features, serum Mb is a hidden important feature in DKD, and that elevated serum Mb is associated with DKD classification. Meanwhile, MetS components also play a key role in helping to build ML decision-making algorithms in DKD.
By using ML algorithm-based variable importance analysis and SHAP, we verified the importance of serum Mb for DKD discrimination, supporting serum Mb as a potential important feature of DKD. Next, we further evaluated the association between serum Mb and DKD using RCS on a continuous scale (Fig. 4A,B). We found that the odds ratio was greater than 1 when serum Mb was greater than about 80, indicating that moderately elevated serum Mb is associated with DKD risk. Similar results were seen in male and female patients. However, when adjusting for eGFR, serum Mb was not significant associated with DKD across the range (Fig S2), indicating a strong correlation between serum Mb and eGFR.
Restricted cubic spline models for log odds ratios of DKD versus serum Mb on a continuous scale. Analyses were adjusted for gender, age, BMI, hyperlipidemia, hypertension, ACCI, and hospitalization date. (A) The red solid line is the multivariable adjusted log odds ratio, and the shaded area shows the 95% confidence interval derived from restricted cubic spline regressions with four knots. The red dashed curve represents the distribution density of the population with different serum Mb levels. (B) The orange and blue solid lines are the multivariable adjusted log odds ratio for female and male subgroups, respectively, and the shaded areas show the 95% confidence intervals derived from restricted cubic spline regressions with four knots. The orange and blue dashed curves show the distribution density of populations with different serum Mb levels in female and male subgroups, respectively.
Based on the aforementioned ML variable importance analysis, SHAP approach, and RCS results, we noticed that serum Mb and MetS components were closely related to DKD. However, studies have shown that serum Mb is elevated in the case of metabolism disorders, and Mb may directly or indirectly cause tissue damage. Thus, we systematically evaluated whether serum Mb mediates renal impairment associated with the MetS components (20 indicators), and the most significant results are shown in Fig. 5. The spearman correlation coefficient between MetS components and among MetS components, serum Mb, and renal function indicators are shown in Fig. S3 and Table S4, indicating that each component was irreplaceable. Renal function impairment (or DKD severity) was evaluated by eGFR and UMALB.
Serum Mb mediates the relationship between MetS components and renal function. (A) A directed acyclic graph (DAG) illustrating the mediation process. (B–F) Results of causal mediation analysis of 5 significant indirect effects (IE) for serum Mb in all eligible patients, DKD or non-DKD subgroups. (B,C) significant IE seen in all eligible patients and DKD subgroup using UMALB as an indicator of renal function. (D,E,F) Significant IE seen in all eligible patients, DKD and non-DKD subgroups, respectively, using eGFR as an indicator of renal function.
Significant IE and TE were found in 8 out of 20 causal mediation models when evaluating renal impairment using eGFR across all observations. Six indicators (HGI, 120 min glucose, fasting glucose, HOMA-IR, IGI/HOMA-IR and IGI in descending order of proportion mediated, the same below) reflecting IR, β-cell function, and glucose metabolism in the DKD subgroup and two indicators (HDL-C, HDL-C/TC) reflecting lipid metabolism in the non-DKD subgroup showed significant IE and TE (see Table 4 for significant results and Table S5 for full results).
Significant IE and TE were found in 3 out of 20 causal mediation models when evaluating renal impairment using UMALB across all observations. 5 indicators (Gutt index, I/G 0 min, 120 min glucose, GHBA1C, and HGI) reflecting IR, β-cell function, and glucose metabolism in the DKD subgroup and no indicator in the non-DKD subgroup showed significant IE and TE (see Table 5 for significant results and Table S6 for full results).
Based on the above systematic evaluation, we support that serum Mb is a potential mediator of MetS component-induced renal function impairment. Generally, compared with UMALB, an indicator of renal function impairment, eGFR has more significant IE and TE as well as higher mediated proportion, indicating serum Mb as a mediator mainly participates in the decline of eGFR, that is, the later phase of DKD progression. The results also suggest that lipid metabolism-induced renal function impairment through serum Mb is important in the non-DKD subgroup, whereas IR, β-cell function, and glucose metabolism are important in the DKD subgroup. These may indicate different causes of elevated serum Mb, and therefore different glycemic and/or lipid control strategies should be employed.
Current DKD prediction models based on large cohorts of T2DM patients, such as ADVANCE study (2012)13, Diabetes Cohort Study (DCS) (2013)49, and National Diabetes Register (NDR) (2011)44, are generally limited by lack of biomarkers and poor performance. A recently observed relationship between elevated serum Mb and DKD suggests that serum Mb may be a risk factor for DKD43,45. However, to our knowledge, no studies have evaluated serum Mb as a risk factor for DKD prediction, nor has serum Mb been used as a risk factor for DKD model construction, partially due to limitations in methods to address collinearity.
We first verified the importance of serum Mb in DKD. To find those features that might be overlooked, the ML algorithms were chosen because they are not constrained by collinearity when processing high-dimensional EHR data. Apart from common DKD-related features, such as renal function features and fasting glucose, we also specially developed ML models including serum Mb. Our DKD models had comparable performance to other models, achieving an AUC of 0.85 in the presence of renal function features and 0.80 after removal of renal function features. Further, in our study, feature importance analysis and SHAP have recognized serum Mb as an important feature. Equally important, RCS found that moderately elevated serum Mb was associated with the risk of DKD, which is consistent with previous findings17,18,19,20. Therefore, we believe that serum Mb is a potential risk factor for DKD and deserves more attention and further study.
Current DKD management focuses solely on blood pressure control, glycemic and lipid control, weight loss, and diet50,51,52, in other words, we emphasize the strong association between MetS components and DKD, but more efforts are needed to find appropriate treatment for DKD53,54. Recently, elevated Mb have been observed in DM and DKD. In a 200 controlled study of 188 DM patients, serum Mb levels in DM patients were elevated and were much higher in DM patients with metabolic syndrome16. In a cross-sectional study of DKD, serum Mb was found to be positively correlated with DKD in a model adjusted for age, gender, cardiovascular risk factors, glycemic control, and blood sugar45. In another study, serum Mb was found to be related to CKD, and elevated serum Mb levels were related to advanced stage of CKD43. Interestingly, the elevation of serum Mb cannot be explained by Statins uses55. The cause of elevated serum Mb may be due to abnormal glucose and lipid metabolism and diabetic muscle tissue damage17,18. In diabetic rats, liraglutide and insulin increased carnitine palmitoyl transferase 1 (CPT1) expression, decreased Mb, and improved cardiovascular function19. In a skeletal muscle-specific (Cpt1bm-/-) mouse model, a 25% fat diet resulted in increased intramuscular lipids and serum Mb and decreased mouse vitality20.
Further, at the molecular level, Mb interacts with glucose to form AGEs33,56, which are active contributors to cellular stress and chronic inflammation34,57. Also, oxidative modification of Mb induced by stress-related protein carbonyl compounds is an important cause for the development of microvascular and macrovascular complications in diabetes21,22,23.
These cumulative evidence indicate that serum Mb has a mediator role in DKD. Confirming the role of serum Mb in DKD progression may provide a thorough comprehension of DKD pathogenesis and a promising therapeutic target for DKD. However, as far as we known, the causal mediation relationship between serum Mb and MetS components and renal function impairment has not been reported.
In our study, feature importance analysis, SHAP, and RCS revealed that MetS components and elevated serum Mb were closely related to DKD. In causal mediation analysis, significant IE and TE were found using either UMALB or eGFR as indicators of renal function impairment, indicating serum Mb as a mediator in MetS component-induced renal function impairment. Thus, we believe that more studies on the role of serum Mb in DKD, based on multilevel research designs and emerging strategies, will provide insights into the pathogenesis of DKD and reveal new treatment targets for DKD.
Our findings should be further verified in cohort studies, as the strength of the evidence we generated is limited by the cross-sectional study design. External validity was difficult to perform in our study because serum Mb was not a commonly used indicator in DKD studies. Further, it is necessary to gain insights into the pathophysiological mechanisms at the cellular and molecular levels to reveal the causal relationship between MetS components, serum Mb, and renal function impairment.
The datasets generated during and/or analyzed during the current study are not publicly available usage due to the other ongoing projects but are available from the corresponding author on reasonable request.
Tanabe, K., Wada, J. & Sato, Y. Targeting angiogenesis and lymphangiogenesis in kidney disease. Nat. Rev. Nephrol. 16(5), 289–303. https://doi.org/10.1038/s41581-020-0260-2 (2020).
Article CAS Google Scholar
Nusinovici, S. et al. Retinal microvascular signs and risk of diabetic kidney disease in asian and white populations. Sci. Rep. 11(1), 4898. https://doi.org/10.1038/s41598-021-84464-7 (2021).
Article ADS CAS Google Scholar
Thomas, M. C. et al. Diabetic kidney disease. Nat. Rev. Dis. Primers 1(July), 1–20. https://doi.org/10.1038/nrdp.2015.18 (2015).
Article Google Scholar
Alicic, R. Z., Rooney, M. T. & Tuttle, K. R. Diabetic kidney disease: Challenges, progress, and possibilities. Clin. J. Am. Soc. Nephrol. 12(12), 2032–2045. https://doi.org/10.2215/CJN.11491116 (2017).
Article CAS Google Scholar
Macisaac, R. J., Ekinci, E. I. & Jerums, G. Markers of and risk factors for the development and progression of diabetic kidney disease. Am. J. kidney Dis. Off J. Natl. Kidney Found 63(2 Suppl 2), S39-62. https://doi.org/10.1053/j.ajkd.2013.10.048 (2014).
Article Google Scholar
Tziomalos, K. & Athyros, V. G. Diabetic nephropathy: New risk factors and improvements in diagnosis. Rev. Diabet. Stud. 12(1–2), 110–118. https://doi.org/10.1900/RDS.2015.12.110 (2015).
Article Google Scholar
Colhoun, H. M. & Marcovecchio, M. L. Biomarkers of diabetic kidney disease. Diabetologia 61(5), 996–1011. https://doi.org/10.1007/s00125-018-4567-5 (2018).
Article CAS Google Scholar
Jiang, W. et al. Establishment and validation of a risk prediction model for early diabetic kidney disease based on a systematic review and meta-analysis of 20 cohorts. Diabetes Care 43(4), 925–933. https://doi.org/10.2337/dc19-1897 (2020).
Article CAS Google Scholar
Slieker, R. C. et al. Performance of prediction models for nephropathy in people with type 2 diabetes: Systematic review and external validation study. BMJ 374, n2134. https://doi.org/10.1136/bmj.n2134 (2021).
Article Google Scholar
Cornier, M.-A. et al. The metabolic syndrome. Endocr. Rev. 29(7), 777–822. https://doi.org/10.1210/er.2008-0024 (2008).
Article CAS Google Scholar
Wu, S. H., Liu, Z. & Ho, S. C. Metabolic syndrome and all-cause mortality: A meta-analysis of prospective cohort studies. Eur. J. Epidemiol. 25(6), 375–384. https://doi.org/10.1007/s10654-010-9459-z (2010).
Article Google Scholar
Mottillo, S. et al. The metabolic syndrome and cardiovascular risk a systematic review and meta-analysis. J. Am. Coll. Cardiol. 56(14), 1113–1132. https://doi.org/10.1016/j.jacc.2010.05.034 (2010).
Article Google Scholar
Jardine, M. J. et al. Prediction of kidney-related outcomes in patients with type 2 diabetes. Am. J. kidney Dis. Off J. Natl. Kidney Found 60(5), 770–778. https://doi.org/10.1053/j.ajkd.2012.04.025 (2012).
Article CAS Google Scholar
Berenbrink, M. Myoglobin’s old and new clothes: From molecular structure to integrated function and evolution. J. Exp. Biol. 213(Pt 16), 2711–2712. https://doi.org/10.1242/jeb.048918 (2010).
Article Google Scholar
Wittenberg, B. A. & Wittenberg, J. B. Myoglobin-mediated oxygen delivery to mitochondria of isolated cardiac myocytes. Proc. Natl. Acad. Sci. U. S. A. 84(21), 7503–7507. https://doi.org/10.1073/pnas.84.21.7503 (1987).
Article ADS CAS Google Scholar
Odum, E. P. & Young, E. E. Elevated cardiac troponin I, creatine kinase and myoglobin and their relationship with cardiovascular risk factors in patients with type 2 diabetes. Diabetes Metab. Syndr. 12(2), 141–145. https://doi.org/10.1016/j.dsx.2017.12.001 (2018).
Article Google Scholar
Miyamoto, T., Watanabe, K., Fukuda, K. & Moritani, T. Near-infrared spectroscopy of vastus lateralis muscle during incremental cycling exercise in patients with type 2 diabetes. Phys. Ther. Res. 23(1), 23–30. https://doi.org/10.1298/ptr.E9984 (2020).
Article Google Scholar
Gildea, N. et al. Influence of type 2 diabetes on muscle deoxygenation during ramp incremental cycle exercise. Respir. Physiol. Neurobiol. 269, 103258. https://doi.org/10.1016/j.resp.2019.103258 (2019).
Article CAS Google Scholar
Huang, Q. et al. Incremental effect of liraglutide on traditional insulin injections in rats with type 2 diabetes mellitus by maintaining glycolipid metabolism and cardiovascular function. Exp. Ther. Med. 17(3), 1863–1869. https://doi.org/10.3892/etm.2019.7148 (2019).
Article CAS Google Scholar
Warfel, J. D. et al. A low fat diet ameliorates pathology but retains beneficial effects associated with CPT1b knockout in skeletal muscle. PLoS ONE 12(12), e0188850. https://doi.org/10.1371/journal.pone.0188850 (2017).
Article CAS Google Scholar
Ghelani, H., Razmovski-Naumovski, V., Pragada, R. R. & Nammi, S. Attenuation of glucose-induced myoglobin glycation and the formation of advanced glycation end products (AGEs) by (R)-α-lipoic acid in vitro. Biomolecules https://doi.org/10.3390/biom8010009 (2018).
Article Google Scholar
Raupbach, J., Ott, C., Koenig, J. & Grune, T. Proteasomal degradation of glycated proteins depends on substrate unfolding: Preferred degradation of moderately modified myoglobin. Free Radic. Biol. Med. 152, 516–524. https://doi.org/10.1016/j.freeradbiomed.2019.11.024 (2020).
Article CAS Google Scholar
You, Y., Liu, F., Du, K.-J., Wen, G.-B. & Lin, Y.-W. Structural and functional alterations of myoglobin by glucose-protein interactions. J. Mol. Model. 20(7), 2358. https://doi.org/10.1007/s00894-014-2358-6 (2014).
Article CAS Google Scholar
Alberti, K. G. & Zimmet, P. Z. Definition, diagnosis and classification of diabetes mellitus and its complications. Part 1: Diagnosis and classification of diabetes mellitus provisional report of a WHO consultation. Diabet. Med. 15(7), 539–553. https://doi.org/10.1002/(SICI)1096-9136(199807)15:7%3c539::AID-DIA668%3e3.0.CO;2-S (1998).
<a data-track="click" rel="nofollow noopener" data-track-label="10.1002/(SICI)1096-9136(199807)15:7<539::AID-DIA668>3.0.CO;2-S" data-track-action="article reference" href="https://doi.org/10.1002%2F%28SICI%291096-9136%28199807%2915%3A7%3C539%3A%3AAID-DIA668%3E3.0.CO%3B2-S" aria-label="Article reference 24" data-doi="10.1002/(SICI)1096-9136(199807)15:7Article CAS Google Scholar
KDOQI. KDOQI clinical practice guidelines and clinical practice recommendations for diabetes and chronic kidney disease. Am. J. kidney Dis. Off J. Natl. Kidney Found 49(2 Suppl 2), S12-154. https://doi.org/10.1053/j.ajkd.2006.12.005 (2007).
Article Google Scholar
Mingrone, G. et al. Metabolic surgery versus conventional medical therapy in patients with type 2 diabetes: 10-year follow-up of an open-label, single-centre, randomised controlled trial. Lancet 397(10271), 293–304. https://doi.org/10.1016/S0140-6736(20)32649-0 (2021).
Article Google Scholar
Rynders, C. A. et al. Comparing simple insulin sensitivity indices to the oral minimal model postexercise. Med. Sci. Sports Exerc. 48(1), 66–72. https://doi.org/10.1249/MSS.0000000000000728 (2016).
Article CAS Google Scholar
Xie, Y. et al. Temporal relationship between body mass index and triglyceride-glucose index and its impact on the incident of hypertension. Nutr. Metab. Cardiovasc. Dis. 29(11), 1220–1229. https://doi.org/10.1016/j.numecd.2019.07.003 (2019).
Article CAS Google Scholar
von Post-Skagegård, M., Vessby, B. & Karlström, B. Glucose and insulin responses in healthy women after intake of composite meals containing cod-, milk-, and soy protein. Eur. J. Clin. Nutr. 60(8), 949–954. https://doi.org/10.1038/sj.ejcn.1602404 (2006).
Article CAS Google Scholar
Lemieux, P. et al. Effects of 6-month vitamin D supplementation on insulin sensitivity and secretion: A randomised, placebo-controlled trial. Eur. J. Endocrinol. 181(3), 287–299. https://doi.org/10.1530/EJE-19-0156 (2019).
Article CAS Google Scholar
Monti, L. D. et al. Decreased diabetes risk over 9 year after 18-month oral L-arginine treatment in middle-aged subjects with impaired glucose tolerance and metabolic syndrome (extension evaluation of L-arginine study). Eur. J. Nutr. 57(8), 2805–2817. https://doi.org/10.1007/s00394-017-1548-2 (2018).
Article CAS Google Scholar
McKay, D. L., Eliasziw, M., Chen, C. Y. O. & Blumberg, J. B. A pecan-rich diet improves cardiometabolic risk factors in overweight and obese adults: A randomized controlled trial. Nutrients https://doi.org/10.3390/nu10030339 (2018).
Article Google Scholar
Hempe, J. M. et al. The hemoglobin glycation index identifies subpopulations with harms or benefits from intensive treatment in the ACCORD trial. Diabetes Care 38(6), 1067–1074. https://doi.org/10.2337/dc14-1844 (2015).
Article CAS Google Scholar
Block, G. et al. Diabetes prevention and weight loss with a fully automated behavioral intervention by email, web, and mobile phone: A randomized controlled trial among persons with prediabetes. J. Med. Internet Res. 17(10), e240. https://doi.org/10.2196/jmir.4897 (2015).
Article Google Scholar
Breiman, L. Random forests. Mach. Learn. 45(1), 5–32. https://doi.org/10.1023/A:1010933404324 (2001).
Article MATH Google Scholar
KDD ’16: Proceedings of the 22nd ACM SIGKDD International Conference on Knowledge Discovery and Data Mining (Association for Computing Machinery, 2016).
Liaw, A. & Wiener, M. Classification and regression by random forest. Forest. 23 (2001).
Chen T, Guestrin C. XGBoost: A Scalable Tree Boosting System. In: Proceedings of the 22nd ACM SIGKDD International Conference on Knowledge Discovery and Data Mining 785–794 (Association for Computing Machinery, 2016). doi:https://doi.org/10.1145/2939672.2939785.
Eisenman, R. L. A profit-sharing interpretation of Shapley value for N-person games. Behav. Sci. 12(5), 396–398. https://doi.org/10.1002/bs.3830120506 (1967).
Article CAS Google Scholar
Lundberg, S. M. & Lee, S. I. A unified approach to interpreting model predictions. Adv. Neural Inf. Process. Syst. 30(2), 4766–4775 (2017).
Google Scholar
Lundberg, S.M., Lee, S.-I. Consistent feature attribution for tree ensembles. http://arxiv.org/abs/1706.06060 (2017).
Richiardi, L., Bellocco, R. & Zugna, D. Mediation analysis in epidemiology: Methods, interpretation and bias. Int. J. Epidemiol. 42(5), 1511–1519. https://doi.org/10.1093/ije/dyt127 (2013).
Article Google Scholar
Lenglet, A. et al. Prognostic implications of plasma myoglobin levels in patients with chronic kidney disease. Int. J. Artif. Organs 35(11), 959–968. https://doi.org/10.5301/ijao.5000124 (2012).
Article CAS Google Scholar
Afghahi, H. et al. Risk factors for the development of albuminuria and renal impairment in type 2 diabetes–the Swedish National Diabetes Register (NDR). Nephrol. Dial. Transplant. Off Publ. Eur. Dial. Transpl. Assoc.—Eur. Ren. Assoc. 26(4), 1236–1243. https://doi.org/10.1093/ndt/gfq535 (2011).
Article Google Scholar
Carlsson, A. C. et al. Growth differentiation factor 15 (GDF-15) is a potential biomarker of both diabetic kidney disease and future cardiovascular events in cohorts of individuals with type 2 diabetes: A proteomics approach. Ups. J. Med. Sci. 125(1), 37–43. https://doi.org/10.1080/03009734.2019.1696430 (2020).
Article Google Scholar
Lemley, K. V. Machine learning comes to nephrology. J. Am. Soc. Nephrol. 30(10), 1780–1781. https://doi.org/10.1681/ASN.2019070664 (2019).
Article Google Scholar
Dimopoulos, A. C. et al. Machine learning methodologies versus cardiovascular risk scores, in predicting disease risk. BMC Med. Res. Methodol. 18(1), 179. https://doi.org/10.1186/s12874-018-0644-1 (2018).
Article Google Scholar
Cowie, M. R. et al. Electronic health records to facilitate clinical research. Clin. Res. Cardiol. 106(1), 1–9. https://doi.org/10.1007/s00392-016-1025-6 (2017).
Article Google Scholar
Elley, C. R. et al. Derivation and validation of a renal risk score for people with type 2 diabetes. Diabetes Care 36(10), 3113–3120. https://doi.org/10.2337/dc13-0190 (2013).
Article CAS Google Scholar
Liew, A. et al. Asian pacific society of nephrology clinical practice guideline on diabetic kidney disease. Nephrology 25(Suppl 2), 12–45. https://doi.org/10.1111/nep.13785 (2020).
Article Google Scholar
Doshi, S. M. & Friedman, A. N. Diagnosis and management of type 2 diabetic kidney disease. Clin. J. Am. Soc. Nephrol. 12(8), 1366–1373. https://doi.org/10.2215/CJN.11111016 (2017).
Article CAS Google Scholar
Barrett, E. J. et al. Diabetic microvascular disease: An endocrine society scientific statement. J. Clin. Endocrinol. Metab. 102(12), 4343–4410. https://doi.org/10.1210/jc.2017-01922 (2017).
Article Google Scholar
Stanton, R. C. Clinical challenges in diagnosis and management of diabetic kidney disease. Am. J. kidney Dis. Off J. Natl. Kidney Found 63(2 Suppl 2), S3-21. https://doi.org/10.1053/j.ajkd.2013.10.050 (2014).
Article Google Scholar
Lytvyn, Y., Bjornstad, P., van Raalte, D. H., Heerspink, H. L. & Cherney, D. Z. I. The new biology of diabetic kidney disease-mechanisms and therapeutic implications. Endocr. Rev. 41(2), 202–231. https://doi.org/10.1210/endrev/bnz010 (2020).
Article Google Scholar
Yasuda, G. et al. Safety and efficacy of fluvastatin in hyperlipidemic patients with chronic renal disease. Ren. Fail. 26(4), 411–418. https://doi.org/10.1081/jdi-120039826 (2004).
Article CAS Google Scholar
Banerjee, S., Maity, S. & Chakraborti, A. S. Methylglyoxal-induced modification causes aggregation of myoglobin. Spectrochim. Acta A Mol. Biomol. Spectrosc. 155, 1–10. https://doi.org/10.1016/j.saa.2015.10.022 (2016).
Article ADS CAS Google Scholar
Goh, S.-Y. & Cooper, M. E. Clinical review: The role of advanced glycation end products in progression and complications of diabetes. J. Clin. Endocrinol. Metab. 93(4), 1143–1152. https://doi.org/10.1210/jc.2007-1817 (2008).
Article CAS Google Scholar
Download references
We would like to thank the Information Network Center of the Third Xiangya Hospital for providing the clinical data.
This work was supported by the National Natural Science Foundation of China (grant no. 81870352, 81970252), the Key Research and Development Project of Hunan Province (grant no.2020SK2087, 2019SK2041), the Undergraduate Training Program for Innovation and Entrepreneurship of Hunan Province (S2020105330657).
Health Management Center, The Third Xiangya Hospital of Central South University, Changsha, 410013, Hunan, People’s Republic of China
Ruoru Wu & Hongwei Lu
Xiangya School of Medicine, Central South University, Changsha, 410013, Hunan, People’s Republic of China
Ruoru Wu, Fei Zou, Saolai Chan & Yaxian Hu
Department of Cardiology, The Third Xiangya Hospital, Central South University, No. 138, Tongzi, Changsha, 410013, Hunan, People’s Republic of China
Zhihao Shu, Shaoli Zhao & Hongwei Lu
Center for Experimental Medicine, The Third Xiangya Hospital of Central South University, Changsha, 410013, Hunan, People’s Republic of China
Hong Xiang
Department of Biochemistry, School of Life Sciences, Central South University, Changsha, 410013, Hunan, People’s Republic of China
Shuhua Chen
Xiangya School of Pharmaceutical Sciences, Central South University, Changsha, 410013, People’s Republic of China
Li Fu & Dongsheng Cao
You can also search for this author in PubMed Google Scholar
You can also search for this author in PubMed Google Scholar
You can also search for this author in PubMed Google Scholar
You can also search for this author in PubMed Google Scholar
You can also search for this author in PubMed Google Scholar
You can also search for this author in PubMed Google Scholar
You can also search for this author in PubMed Google Scholar
You can also search for this author in PubMed Google Scholar
You can also search for this author in PubMed Google Scholar
You can also search for this author in PubMed Google Scholar
You can also search for this author in PubMed Google Scholar
R.W., Z.S., H.L., D.C., F.Z., S.Z., S.C. conceived and designed the study. R.W., Z.S., F.Z., S.Z., S.C., Y.H., H.X. coordinated data collection. R.W., F.Z., Z.S., S.Z., L.F. conducted the analyses. R.W., Z.S., F.Z., H.L., L.F. wrote the manuscript. All co-authors critically revised the manuscript.
Correspondence to Hongwei Lu.
The authors declare no competing interests.
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article’s Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article’s Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http://creativecommons.org/licenses/by/4.0/.
Reprints and Permissions
Wu, R., Shu, Z., Zou, F. et al. Identifying myoglobin as a mediator of diabetic kidney disease: a machine learning-based cross-sectional study. Sci Rep 12, 21411 (2022). https://doi.org/10.1038/s41598-022-25299-8
Download citation
Received: 18 March 2022
Accepted: 28 November 2022
Published: 10 December 2022
DOI: https://doi.org/10.1038/s41598-022-25299-8
Anyone you share the following link with will be able to read this content:
Sorry, a shareable link is not currently available for this article.
Provided by the Springer Nature SharedIt content-sharing initiative
By submitting a comment you agree to abide by our Terms and Community Guidelines. If you find something abusive or that does not comply with our terms or guidelines please flag it as inappropriate.
Advertisement
© 2022 Springer Nature Limited
Sign up for the Nature Briefing newsletter — what matters in science, free to your inbox daily.