ABSTRACT
Objective The methods currently used to calculate GFR underestimate this measurement in the population of diabetic patients, so there is a need to look for more accurate methods of estimating GFR in this specific population. This study aims to evaluate a predictive model based on the use of HbA1c to estimate the variability of GFR in diabetic patients with or without CKD.
Methods We analyzed data from diabetic patients belonging to a cohort of prospective follow-up of a renal health surveillance program attached to a Peruvian hospital. The following factors were included in the multiple linear regression model: age, sex, diastolic blood pressure (DBP), systolic blood pressure (SBP), body mass index (BMI), cholesterol, triglycerides, HDL, LDL, serum creatinine, Urinary creatinine, microalbuminuria, hemoglobin, basal glycemia and HbA1c.
Results 122 patients were included in the analysis. The final multivariate model, which included variation of HbA1c, age and creatinine variation, was very significant (p <0.0001) with an adjusted R2 of 80%. The other variables analyzed were not significant to predict the variation of the GFR despite showing some correlation.
Conclusions The study shows that HbA1c, age and creatinine variation significantly predict the variation of GFR in diabetic patients with or without CKD and opens the possibility of use as a prognostic tool for this specific population.
INTRODUCTION
Chronic kidney disease (CKD) affects about 50% of the global diabetes population [1], with rates ranging from 40% to 83% in developed countries [2-4]. It is estimated that the prevalence of CKD in patients with diabetes will continue to increase dramatically over the next decade, regardless of their stage [5].
In Peru, up to 59.09% of people with type 2 diabetes are affected by CKD. We also know that it is the main cause of end-stage kidney disease and that it requires dialysis. The national rate of Peruvian prevalence per million population for renal replacement treatment is 335.3, with hemodialysis, peritoneal dialysis, and functional kidney transplant rates of 230.7, 39.1 and 65.5, respectively. Diabetes represents 40% of haemodialysis patients in a public hospital [6]. This requires an adequate capacity to evaluate and monitor the evolution of CKD in this specific population.
The use of glomerular filtration rate (GFR) is undoubtedly the most commonly used tool for detecting, managing and predicting the prognosis of CKD in patients with metabolic diseases such as type 2 diabetes mellitus [7]; however, some studies show that the methods for its estimation such as the formula MDRD (Modification of Diet in Renal Disease) or the formula of the CKD-EPI (Chronic Kidney Disease Epidemiology Collaboration) may be inadequate in diabetic patients due to an underestimation of these values [8,9].
In addition, other studies suggest that the values of Glycosylated Hemoglobin (HbA1c), an essential test for glycemic monitoring and control in diabetic patients, may be related to the variation in their GFR [10, 11], an inverse correlation between the increase in HbA1c and the decrease in GFR [12] was reported, even in the population of non-diabetic patients [13, 14]. A study assessing GFR predictive models using new parameters showed that the use of HbA1c has a high percentage of accuracy for predicting GFR [15].
The use of HbA1c to predict the variation of GFR would be very useful, given the need to find new alternatives for estimating GFR in diabetic patients without the need to increase costs, since these tests are part of the monitoring and control of all diabetic patients.
The aim of this study was to evaluate a predictive model for estimating GFR variability in diabetic patients based on HbA1c variability through a multiple linear regression analysis, incorporating into the analysis other variables that may influence GFR variability such as age, sex, blood pressure, BMI and routine laboratory tests.
MATERIAL AND METHODS
The study is based on a cohort of patients with prospective follow-up belonging to the Adult Surveillance Program with Non-Transmission Chronic Disease, specifically, those who belong to the Renal Health Surveillance Program (VISARE) of the Hospital-I Pacasmayo of Essalud, La Libertad - Peru, a program that is implemented nationally for patients who are treated in Essalud.
This cohort has regular quarterly controls for the follow-up of patients with a diagnosis of Diabetes Mellitus type 2 with or without CKD, having a medical evaluation every 3 months, and monthly checks by the nursing staff for taking AP, BMI and programming of laboratory tests already stipulated as part of the program of Renal Health of the diabetic patient by Essalud, which include monthly evaluation of fasting glycemia, half-yearly evaluation of HbA1c, annual evaluation of lipid profile, hemoglobin, microalbuminurea, urine creatinine and serum creatinine.
The use of data was authorized by the Hospital I Pacasmayo to use data from the VISARE program and to have access to patients’ medical records. Information was subsequently collected from the hospital database for the period from 1 January 2014 to 31 December 2015.
We included patients who presented continuity in their medical controls for this program. The total number of patients identified according to the database obtained was submitted to the study selection criteria.
Selection criteria
Inclusion criteria: Patients over 18 years of age with confirmed diagnosis of Diabetes Mellitus type 2 and belonging to the Renal Health surveillance program with at least an annual control in said program; patients who record age data, sex, laboratory tests at the time of periodic control in the program, including: serum creatinine, urine creatinine, microalbuminuria, hemoglobin, basal glycemia and HbA1c; diastolic blood pressure (DBP), systolic blood pressure (SBP), body mass index (BMI), GFR, use of nephroprotective drugs and comorbidities. Exclusion criteria: Patients with missing or poorly recorded data in any of the annual controls considered for the study.
After screening patients according to the selection criteria, information was collected from medical records. This information was encoded for analysis, which was executed using Stata 15.0 statistical software (StataCorp, College Station, TX, USA), performing a multiple linear regression that included the factors previously mentioned as possible predictors of GFR variability.
The result variable of the linear model was the variation of the TFG in the period of 1 year, which was calculated according to the values reported from the results for the years 2014 and 2015. Although many patients had quarterly controls, annual GFR variability was taken into consideration to better assess the plausibility of the factors included in this study, since evaluating short periods of these factors would have less influence on the variation of the TFG.
The variables blood pressure and basal glycemia were defined as the average of at least 3 measurements in consecutive months prior to the control by the program, due to their greater variability by various factors.
Statistical analysis
The descriptive analysis was performed for each factor studied by calculating proportions for the categorical variables, as well as the mean, standard deviation and median for the numerical variables; additionally, the p value was reported to evaluate the difference between groups in the baseline analysis.
A Pearson correlation analysis was performed for variables of normal distribution and the Spearman test for variables without normal distribution. The analysis to determine the distribution of normality of the variables was performed using the Shapiro-Wilk test.
The inclusion of the factors in the multiple linear regression model was performed after evaluating the correlation between each factor studied and the variation of the GFR, the model was adjusted with the factors that showed to have a significant correlation.
Ethical aspects
The research did not involve the direct participation of human beings (biological samples or interventions), it was worked with a hospital database authorized for use by the management of the institution. The handling of the information collected by the investigators was confidential and no personal data of any kind will be disclosed.
RESULTS
We identified 234 patients, of whom only 122 diabetic patients met the selection criteria and were included in the multiple linear regression analysis to evaluate the predictive capacity of HbA1c, serum creatinine, and age. Of the study participants, the majority were male (55.7%) with no significant differences with respect to females (p= 0.20); the mean age was higher for men (67.4 years, C.I.: 65.2 - 69.6) than for women, there were significant differences between the two genders (p = 0.006). Of all patients, 63 (51.6%) had associated comorbidities, mainly hypertension (HT), with no difference between groups (p= 0.72); 76 (62.3%) used a nephroprotective drug (ACE inhibitor, ARA inhibitor or both), significantly different from those who did not use nephroprotective drugs (p= 0.006).
Regarding the BMI variation, there were no differences between both sexes (p= 0.16), with the male sex showing a variation of 0.02 (I.C.: -0.31, 0.28) and the female sex a variation of 0.33 (I.C.: -0.67, -0.001) (Table 1).
Baseline of the study cohort of diabetic patients.
Variables with normal distribution included Hb variation, serum creatinine variation, and HbA1c variation.
The correlation analysis showed that only HbA1c variation, serum creatinine variation, and microalbuminuria variation showed significant correlation when evaluated for GFR variability (Table 2).
Correlation of the factors analyzed with respect to the variation of the GFR.
Multiple linear regression predictive model for GFR variation
Multivariate Analysis
the final multiple linear regression obtained after selecting the best adjusted model for predicting GFR variability was that which included HbA1c variation, serum creatinine variation, and age (Figures 1, 2, 3). The coefficients for the model were: -1.94 for HbA1c variation (p= 0.006), -108.81 for serum creatinine variation (p < 0.0001) and -0.04 for age (p= 0.009). The final model was very significant (p < 0.0001) with a variation of the HbA1 variance that explains 80% of the GFR (adjusted R2 = 0.80). The determination of the best model led to the elimination of the constant from the final formula as not being significant to the model.
Linear prediction of GFR variation according to HbA1c variation
Linear prediction of GFR variation according to Creatinine variation
Linear prediction of GFR variation by age.
DISCUSSION
The prevalence of diabetes is increasing, for 2014 it was estimated that there are 387 million cases worldwide [16], compared to previous decades [17]. There are no reliable data on the prevalence of diabetes in Peru, the few published reports show important percentages [18, 19]; of these, the percentage of patients with CKD is also important [6]. Diabetes in developed and some developing countries has become the leading cause of CKD. Therefore, it is necessary to accurately evaluate the GFR variation in diabetic patients, whether or not they have CKD.
Our analysis showed that HbA1c variation, serum creatinine variation, and age significantly predict GFR variation in diabetic patients regardless of their CKD condition. In research by Tsuda et al, a new equation was developed that includes HbA1c as a new variable. However, the sample size was very small (40 cases) [20]. The number of patients in our study was much higher, and the performance according to our study was very good.
The multiple regression model included an adequate number of patients, although the sample corresponded to a single center. We must also clarify that the variables, creatinine variability and age, are variables that show collinearity with the formula of Cockcroft - Gault and include them in the model was obviously necessary.
The results showing an explanatory level of 80% for the final model are quite satisfactory for this specific population, the model that proved to be very significant discusses the need to improve prognostic tools for CKD in these patients; such improvements should take into account the underlying disease and therefore try to use tests that are already part of the periodic control of the disease such as HbA1c.
It is important to note the limitations of our study. One of them is having performed the analysis on subjects that were restricted to a single center, under a renal control program. It is also important to mention that the study cohort was not completely uniform for age and nephroproprotection variables, there being a statistical difference when analyzed with respect to the sex of the patient, The average age in male patients is clearly higher than in the female sex, in the same way the proportion of patients who received some nephroproprotective drug is almost double that of those who do not receive any drug. These differences in baseline may have influenced the final model obtained since the relationship of GFR variation with age is known [21, 22], as well as the obvious effect of neproprotective drugs (ACE inhibitors and ARAs) on the evolution of GFR. Despite this, a very significant predictor model was achieved and it is necessary to emphasize that studies such as these are necessary in our environment, since many of the GFR estimation formulas are based on foreign populations with very different realities in terms of the parameters commonly used and the same idiosyncrasy of each individual.
CONCLUSIONS
The optimal study model for predicting GFR variation in diabetic patients with or without CKD included HbA1c, serum creatinine variation, and age, with a high percentage of predictability, however it is important to note that the HbA1c correlation level was low and there were factors that could have affected the results in the study population, such as the differentiated use of nephroproprotective drugs by much of the study cohort and the age difference by sex. New research studies are required that consider a much larger sample and standardize criteria to eliminate confounding factors that could alter the predictive capacity of HbA1c for prognosis of GFR variation. It is possible to further improve the predictive capacity of HbA1c by improving the analysis methodology.
Data Availability
All data produced in the present study are available upon reasonable request to the authors