Skip to main content

The effect of comorbidities on glycemic control among Colombian adults with diabetes mellitus: a longitudinal approach with real-world data



Achieving an optimal glycemic control has been described to reduce the incidence of diabetes mellitus (DM) related complications. The association between comorbidities and glycemic control remains unclear. Our aim is to evaluate the effect of comorbidities on glycemic control in people living with DM.


A retrospective longitudinal study on data from the National Registry of Chronic Kidney Disease from 2014 to 2019 in Colombia. The outcome was poor glycemic control (PGC = HbA1c ≥7.0%). The association between each comorbidity (hypertension (HTN), chronic kidney disease (CKD) or obesity) and PGC was evaluated through multivariate mixed effects logistic regression models. The measures of effect were odds ratios (OR) and their 95% confidence intervals (CI). We also evaluated the main associations stratified by gender, insurance, and early onset diabetes as well as statistical interaction between each comorbidity and ethnicity.


From 969,531 people at baseline, 85% had at least one comorbidity; they were older and mostly female. In people living with DM and CKD, the odds of having a PGC were 78% (OR: 1.78, CI 95%: 1.55-2.05) higher than those without CKD. Same pattern was observed in obese for whom the odds were 52% (OR: 1.52, CI 95%: 1.31-1.75) higher than in non-obese. Non-significant association was found between HTN and PGC. We found statistical interaction between comorbidities and ethnicity (afro descendant) as well as effect modification by health insurance and early onset DM.


Prevalence of comorbidities was high in adults living with DM. Patients with concomitant CKD or obesity had significantly higher odds of having a PGC.

Peer Review reports


Diabetes Mellitus (DM) prevalence is increasing worldwide [1]. The International Diabetes Federation estimated that 9.3% (463 million people) of the global adult population had DM in 2019; with an unsettling projection of 10.2% (578 million people) and 10.9% (700 million people) for 2030 and 2045, respectively [1]. The same phenomena was projected for the South and Central American region; with a prevalence up to 8.5% (31.6 million people) and a prediction up to 9.9% (49.1 million people) in 2045 [1]. The aforementioned region includes Colombia, a country where DM prevalence ranges from 1.8-11.2% depending on the report consulted [2]. This is alarming given the effects DM has on morbidity and the sustainability of health systems [3].

Achieving an optimal glycemic control has been described to decrease the incidence of DM related complications [4,5,6,7]. As for microvascular complications, in a recent meta-analysis a more rigorous control was found to reduce the relative risk by 20 and 13% for kidney and eye outcomes, respectively [5]. On the other hand, an effect on macrovascular complications was also described with a reduction in cardiovascular events when intensive glycemic control is employed [6, 7]. In spite of this, only few patients reach the established target; a meta-analysis described that approximately 36% of patients in randomized control trials and 34% in cross-sectional studies achieved the hemoglobin A1c (HbA1c) target (≤7.0% or < 7.0%, respectively) [8].

Different sociodemographic and clinical factors appear to be associated with glycemic control. Regarding age and gender, female and younger patients were more commonly associated with poor control [9]. Taking into account race and ethnicity, studies have demonstrated significant differences between ethnic groups in achieving glycemic target [10, 11]. An additional determinant of inadequate glycemic control is a longer time since the diagnosis of DM [9, 12].

Comorbid conditions, mainly hypertension (HTN), obesity, hyperlipidemia, chronic kidney disease (CKD) and cardiovascular disease, are common in patients living with DM [13]. Studies evaluating the relationship between the presence of comorbidities and glycemic control have showed contradictory findings depending on the analysis. When the total number of comorbidities was used, studies have found no relationship [14] or an inverse relationship [15] between comorbidities and HbA1c. Research using comorbidity scores also established no relationship with HbA1c levels [16]. Whereas in a study evaluating specific comorbidities, subjects with HTN and dyslipidemia were found to have lower odds of achieving glycemic target while no significant association was observed with CKD [17]. Lastly, patients living with DM and overweight or obesity were more commonly associated with poor glycemic levels [18].

Considering the detrimental effects of suboptimal glycemic control, the identification of patients at an increased risk of PGC is crucial for health systems and policy makers as an aid to develop strategies to capture at-risk individuals. Therefore, we aimed to assess the longitudinal effect that common DM-related comorbidities (HTN, CKD, and obesity) have on glycemic control among Colombian adults diagnosed with DM from 2014 to 2019.


Data sources

We conducted a retrospective longitudinal analysis on data of people diagnosed with DM who were reported by health insurers and providers to the National Registry of Chronic Kidney Disease (NRCKD) from July 1st, 2013 to June 30th, 2019 in Colombia. The NRCKD is managed by the High-Cost Diseases Fund (“Cuenta de Alto Costo” in Spanish). It has been operating since 2008 by a resolution from the Ministry of Health of Colombia [19]. The NRCKD aims to evaluate the burden of CKD and its most common precursors (HTN and DM) as well as the effective access to health services related to prevention, diagnosis and control across the country.

The NRCKD is an administrative and passive registry with a national scope because approximately 98% of the population is affiliated to the national healthcare system [20], and their insurers are required to report patients living with HTN, DM or CKD to the registry [19]. A unique identification number is assigned to each individual protecting their personal information and allowing future follow up. When a new case enters the NRCKD, a complete registration is done; for old cases, data are updated every year.

The NRCKD undergoes a data auditing process to ensure the veracity of the information. The first step involves the use of an algorithm to identify mistakes in the reporting process. After this, an experienced team compares the reported information with health clinical records by a well-stablished data monitoring process in a representative stratified sample of cases with HTN and/or DM, with or without CKD. In case of any inconsistency is identified, the real data on clinical records is captured.

Eligibility of participants

A total of 1,081,863 people aged 18 years or more, diagnosed with DM were reported to the NRCKD from 2014 to 2019. 112,332 were excluded because they had no information available on the measurement of HbA1c. Indeed, 969,531 people met the inclusion criteria and were analyzed. Figure 1 shows the eligibility flow of the population studied.

Fig. 1
figure 1

Eligibility flow of studied population from the National Registry of Chronic Kidney Disease, Colombia 2014-2019. Abbreviations: DM: diabetes mellitus, NRCKD: National Registry of Chronic Kidney Disease, HbA1c: glycated hemoglobin

Exposure and outcome variables

The presence of comorbidities (HTN, CKD, obesity), reported by health insurers to the NRCKD was considered as the exposure of interest. HTN diagnosis was reported to the NRCKD as its presence or absence. This diagnosis is defined by the treating physician based on his/her clinical judgement [21, 22]. CKD diagnosis was defined by the patients’ treating physician and reported to the NRCKD by their insurers as follows: yes/no/undetermined/not studied for CKD during the reporting period. CKD diagnosis was verified by a data monitoring process following an algorithm based on the Colombian CKD guideline that defines the presence of CKD as abnormalities in the structure or function (GFR < 60 ml/min/1.73 m2) for more than 3 months [23]. People classified as undetermined were those with early kidney impairment that require evidence of functional or structural injury for more than 3 months. Lastly, obesity was determined using the World Health Organization definition (body mass index [BMI] ≥30 kg/m2) [24]. Exhaustiveness and accuracy of exposures were verified in clinical records by a well-established data monitoring process.

A poor glycemic control (PGC) was the dependent variable. It was defined by the measurement of HbA1c, reported during the last 6 months of each follow up period. PGC was treated as a dichotomous variable for each point in time using the cutoff of HbA1c ≥ 7.0%, according to the standard defined by the American Diabetes Association [25].


Demographic information included age, gender, health insurance (third payer, state insurance, special insurance, exception insurance or uninsured), and region of residence following the Administrative Department for National Statistics (“Departamento Administrativo Nacional de Estadística” in Spanish) classification (Bogotá, D.C., Central, Eastern, Pacific, Caribbean and Amazonian) [26] (Additional File 1: Table S1). Ethnicity was self-reported and classified in indigenous, afro descendant and other (white, mestizo, gypsy or not classified). In respect of clinical data, age at diagnosis of DM and comorbidities were included. Reported information on weight (kg) and height (meters) was used to calculate BMI (kg/m2). Otherwise, glomerular filtration rate (GFR) was estimated using the CKD-EPI method and expressed as mL/min/1.73m2 [27]. On the other hand, CKD stages were defined based on the GFR as follows: 1: GFR > 90 with clinical evidence of kidney damage; 2: GFR: 60-90 with clinical evidence of kidney damage; 3: GFR: 30-60, 4: GFR 15-30, 5: GFR: < 15 [28]. Lastly, DM duration was calculated using the date of last contact (last HbA1C test) and the date of diagnosis and informed in years.

Statistical analysis

Descriptive analysis

Baseline demographic and clinical characteristics of the cohort are presented as means or medians and their standard deviation (SD) or interquartile range (IQR), according to their distribution, evaluated by statistical and graphical methods. Otherwise, qualitative data is informed as absolute values (proportions). Baseline characteristics were compared across comorbidities (HTN, CKD and obesity) by using a X2 test for categorical data and Wilcoxon test was used for numeric variables. We also calculated absolute standardized differences (ASD) between groups of comorbidities to avoid limitations of hypothesis tests taking into account the effect of the sample size. In the case of continuous variables with skewed distribution, rank statistics were used instead the mean [29]. ASD > 10% were considered significant.

Assessment of HbA1c trend

Furthermore, the trend of variation on HbA1c levels during the follow-up for patients with different baseline comorbidity profiles was evaluated by estimating generalized linear mixed models assuming random intercepts for each participant. An interaction term between HbA1c levels and time was included to take into account variations over time, adjusting by each comorbidity. From each model, the predicted Hb1Ac values were plotted over time.

Evaluation of the association between comorbidities and poor glycemic control

The association of interest was estimated using a mixed effects logistic regression model with aleatory intercept, which modelates the odds of PGC for each point in time, taking into account the effect of repeated measures and considering the parameters related to the intercept differentially between each participant, but with a constant slope. The final model for each comorbidity was adjusted by age, gender, ethnicity, health insurance, duration of diabetes. In case of HTN and CKD we also included the BMI (time variant) as a confounder in the final models. Effect measures were the odds ratios and their 95% confidence intervals. We also repeated the final models stratifying by gender, health insurance, and early onset DM (diagnosis < 45 years old [30]) to evaluate a potential effect modification.

Estimation of interaction models

Furthermore, the interaction between each comorbidity and intermediate variables such us ethnicity was explored by the inclusion of a multiplicative term in the final models. p values < 0.05 (two-tailed) were considered statistically significant, except for the effect modification model (p-value < 0.20).

General considerations

Trend and association models were estimated using a case complete approach including people with 6-year HbA1c measurement (n = 107,091). All statistical analyses were performed in Stata version 13 (StataCorp LP, College Station, Texas, USA) and graphs were created in R version 4.0.2.

Sensitivity analysis

We performed sensitivity analysis including a BMI ≥25 kg/m2 instead of ≥30 kg/m2 as the exposition to verify the consistency of the association obtained in the main model. Furthermore, we evaluated the association in people with CKD by using a different cut off to define the outcome (HbA1c ≥ 8.0%). Taking into account the potential effect of early onset on PGC, it was tested as a covariate in the final models.


Demographic and clinical characteristics at baseline

969,531 people were analyzed. Mean age at baseline was 61.4 years (SD 13.2), 58.1% were female and 75.7% had third payer insurance. Most people (29.5%) lived in the Central region, followed by Bogotá, D.C. (20.3%). Around 96.0% were self-recognized as “other” ethnicity and 3.5% were afro descendant. Mean age at DM diagnosis was 59.2 years (SD 13.0), median of DM duration was 2.76 years (IQR 1.0-6.0), 56.8% had DM evolution < 5 years and 9.9% had early onset DM. Prevalence of HTN, CKD and obesity were 71.5, 30.6 and 35.2%, respectively. Majority of the people with CKD were classified in stages 2 (43.9%) and 3 (28.9%).

Table 1 shows a characterization by comorbidities. Only 14.9% did not have an associated comorbidity at baseline. Prevalence of comorbidities was higher in women. People with comorbidities were older, except for those with obesity. Age at DM diagnosis was also higher in people with comorbidities. ASD were significant for DM duration in all groups, except for the obese. When comparing the prevalence of obesity by comorbidities, it was also higher in people with HTN and CKD.

Table 1 Baseline characteristics of the population studied according to comorbidities in Colombia 2014-2019

Trends of longitudinal variation on HbA1c by comorbidities

After a 6-year follow up, people living with DM and HTN or CKD had a better glycemic control (lower HbA1c levels) than those with obesity. In both DM subjects with HTN or CKD, HbA1c levels were lower in those who had the condition. These findings were contrary to the obese/not obese subjects, where we found higher HbA1c levels when obesity was present. A higher gap in the HbA1c levels between the presence or not of comorbidities was observed in people living with DM and HTN/no HTN than in the other groups. When estimating the longitudinal variation on HbA1c when changing the cut-off from BMI ≥ 30 kg/m2 to 25 kg/m2, the trend was opposite to the observed in obese, with a better glycemic control in those with a BMI ≥ 25 kg/m2 (Fig. 2).

Fig. 2
figure 2

Longitudinal trend of glycated hemoglobin by comorbidities. Colombia 2014-2019. Figure 2 shows the variation on HbA1c over time by comorbidities. In general, there was a pattern on HbA1c levels in people with diabetes and HTN (a) or CKD (b), characterized by higher HbA1c in those without the comorbidities. However, an inverse trend was observed in obese (c). Abbreviations: HTN: hypertension, CKD: chronic kidney disease, HbA1c: glycated hemoglobin

Longitudinal association between comorbidities and glycemic control

The mean number of HbA1c measurements was 3.18 (SD 1.7). Only 11.0% of the patients had 6-year HbA1c measurements and were included in the final models. Global prevalence of PGC was 52.7% (95% CI: 52.6 a 52.8). When analyzing the effect of comorbidities on PGC (HbA1c ≥7.0), we observed no significant association between HTN and PGC. Nevertheless, CKD and obesity were significant and directly associated with PGC in all models. In subjects living with DM and CKD, the odds of having a PGC were 78.0% (OR: 1.78, CI 95%: 1.55-2.05) higher than those without CKD. We also observed a strong direct association among people living with DM and obesity, in whom the odds were 52.0% (OR: 1.52, CI 95%: 1.31-1.75) higher than in non-obese (Table 2).

Table 2 Multivariate-adjusted odds of poor glycemic control by comorbidities in Colombian adults with diabetes mellitus, 2014-2019

In stratified models we did not observe a significant difference by gender. Regarding health insurance, people cover by the state were more likely to achieve a better glycemic control compared with those insured to the third payer, but this effect was only significantly different in CKD. People with obesity and early onset DM had higher odds of having PGC and this effect was significantly different than observed in those diagnosed ≥45 years (Table 2).

When evaluating the interaction between comorbidity and ethnicity, we only identified a statistically significant interaction between each comorbidity and afro-descendants (all p-values were < 0.01). Predictive probabilities of having a PGC were significantly lower in afro-descendant with HTN compared with afro descendant without HTN. Otherwise, an inverse effect was found in afro descendant with CKD who had a significantly higher probability of a PGC. Afro descendant with obesity, had a lower probability of a PGC, following the same pattern observed with HTN, however, the effect was not significant (Fig. 3).

Fig. 3
figure 3

Predicted probabilities of poor glycemic control according to comorbidities and ethnicity, Colombia 2014-2019. Figure 3 shows the predicted probabilities of poor glycemic control according to the presence of comorbidities and ethnicity estimated in a multivariate model adjusted by for age, gender, BMI (only for HTN and CKD), health insurance and diabetes duration. Predicted probabilities of poor glycemic control were significantly lower in afro descendant with HTN (a). An opposite association was found in afro descendant with CKD (b). In afro descendant with obesity, the association was similar to those with HTN although it was not statistically significant (c). Abbreviations: HTN: hypertension, CKD: chronic kidney disease, HbA1c: glycated hemoglobin, BMI: body mass index

Sensitivity analysis

The associations of interest were consistent with the main models (Additional File 2: Table S2). When analyzing the effect of a BMI ≥ 25 kg/m2 on the outcome, we observed the same direction and significance but a lower magnitude in respect of the final model including obesity (OR: 1.26, CI 95%: 1.06-1.49). Similarly, direction and statistical significance remained, although the magnitude decreased when the cut-off point of HbA1c was slightly higher to define the PGC (OR: 1.30, CI 95%: 1.11-1.52). Final models adjusted by early onset DM did not show substantial changes, except for HTN in which direction of the association was inverse although it remained non-significant. Therefore, we hypothesized a potential modifying effect and estimated stratified models as previously described.


This study evaluated the effect that common diabetes associated comorbidities (CKD, HTN and obesity) have on glycemic control among a Colombian adult population living with DM. Longitudinal, non-adjusted trend of HbA1c over time showed a better glycemic control in people with HTN, CKD and non-obese. However, when estimating multivariate models, an inverse effect was found, and people living with diabetes and CKD or obesity had significantly higher odds of having a PGC compared to those without the comorbidities. Whereas a non-significant association was found between HTN and a PGC.

Our study showed that CKD was associated with 78% increased odds of having a PGC. Taking into account that a less strict goal may be recommended in patients living with DM and advanced comorbid diseases [25, 31], we conducted a sensitivity analysis with a less strict threshold (PGC as HbA1c ≥8.0%). Even with this flexible target, the findings were consistent. Patients living with DM and CKD had an increase likelihood of PGC when compared to those without CKD even with a less stringent goal.

Interestingly, we found statistical heterogeneity in the effect of CKD on PGC by health insurance with 84% decreased odds of PGC in people under state insurance whereas third payer affiliates maintained the same trend of the general findings. Previous studies have evaluated the effect of health insurance on glycemic control [32, 33]. A cross-sectional study in Mexico showed that public health affiliates were more likely to have a better glycemic control than uninsured subjects [32]. Similar findings were described by Zhang et al. in the USA where uninsured subjects were more likely to have HbA1c > 9% [33]. In Colombia, both types of insurance have similar access to healthcare services and medications [34], complicating comparisons with other health systems. A possible explanation to our findings could be a higher proportion of stage 5 CKD in the state insurance compared with third payer affiliates (10% versus 1.81%) since the precision of HbA1c reduces with advanced CKD stages and this can lead to underestimation of glycemic levels [31, 35]. Furthermore, when evaluating the interaction between CKD and ethnicity, we found that afro descendants living with DM and CKD also had an increased probability of PGC compared to those without CKD.

Our overall findings are consistent to those found by De Cosmo et al. where CKD (defined as low GFR, high albuminuria or both) was associated with a lower probability of achieving the HbA1c target established in a large Italian population living with DM [36]. However, others have reported that people living with DM and CKD had lower odds of having suboptimal glycemic control [37]. It is important to note that despite the previously described limitations, HbA1c is still considered the standard for glycemic follow-up in patients with CKD [31, 38] thus giving relevance to our results. These findings should increase the concern of clinicians taking care of patients with kidney dysfunction given the consequences of a PGC have on DM patients with CKD. A meta-analysis of patients living with DM treated with hemodialysis described an increased risk of mortality in people with HbA1c ≥8.5% when compared to lower levels (6.5-7.4%) [39]. Also, Kuo et al. observed that HbA1c levels > 9% were associated with an increased risk of end stage renal disease and mortality in patients living with diabetes and CKD stages 3-4 [40].

Similar to the CKD effect on PGC, we observed that obese people were 52% more likely to have a PGC than those without obesity. When a lower threshold was used (≥25 kg/m2), the results followed a similar trend, but magnitude decreased to 26%. Together, these results indicate that an elevated BMI increases the risk of PGC in patients living with DM. Similar conclusions have been described in previous studies [37, 41]. A cross sectional study in Iranian population living with diabetes found an association between obesity and poor glycemic control [41]. Furthermore, Bae et al. describe an increased probability of suboptimal glycemic control when overweight or obesity was present in patients living with DM in the USA [37]. Nevertheless, Chetoui et al. did not found an association between overweight/obesity and poor glycemic control in a Moroccan population living with type 2 DM [12]. In a therapeutically point of view, Gummesson et al. stated that weight reduction in overweight and obese patients living with type 2 DM is consistently followed by a reduction on HbA1c levels [42].

Moreover, when evaluating the effect of comorbidities on PGC stratified by age of diagnosis (< 45 years and ≥ 45 years), we found that people with obesity and early onset DM had higher odds of having PGC. These findings are of interest since previous studies have identified a high prevalence of obesity in patients with early onset DM [43, 44]. Additionally, Gopalan et al. described that patients with a diagnosis at 21-44 years had higher initial HbA1c levels and lower odds of achieving glycemic control after a 1 year follow up compared with those diagnosed at 45-64 years [45]. Therefore, these findings support the importance of early-onset DM to further implement specific cost-benefit strategies within this high-risk group to avoid the long-term deleterious complications of PGC.

Some authors have described the presence of HTN as a predictor of poor glycemic control [17, 41]. However, we did not find a significant association between HTN and PGC. Our results are consistent with those of Badedi et al. describing no association between HTN and HbA1c [46]. Though, this was different when evaluating the interaction between ethnicity and HTN. We found that afro descendants living with DM and HTN had a lower predictive probability of PGC than those without HTN; other ethnicities did not have a significant association. Various authors have described an association between ethnicity and glycemic control [10, 11]. A cross sectional study in the USA, found that non-Hispanic blacks and Mexican Americans when compared to non-Hispanic whites had lower odds of adequate glycemic control [11]. Additionally, the same study described that the proportion of patients with three comorbidities (HTN, Hyperlipidemia and obesity) achieving glycemic control (HbA1c < 7.0%) changed according to their ethnicity: 22.8, 11.4 and 0% of non-Hispanic whites, non-Hispanic blacks and Mexican Americans, respectively [11].

High prevalence of PGC could negatively impact health systems. In people living with DM, higher HbA1c levels were associated with an increase in the rate of DM-related hospital admissions [47]. Following a similar trend, Gil et al. evaluated the effect of uncontrolled DM (HbA1c ≥ 7.5%) on health care utilization and found an increase in the number of general practitioner and specialist visits as well as hospital length of stay among male patients living with DM [48]. Furthermore, a cross-sectional study in Spain described that people with PGC had a higher total health-related costs including medication and hospitalization costs when compared to those controlled [49]. The above suggests that a higher proportion of people with PGC may significantly increase the burden on healthcare systems, with a worse impact in low-and-middle-income countries. Thus, our results should encourage public health decision-making for early identification, increased access to health care, and surveillance in subjects living with DM at a higher risk of PGC due to the coexistence of obesity or CKD.

This study has important strengths. First of all, this report reflects a national scope since data come from the largest registry of people with CKD in the country. By including people with different comorbidities, our results are representative of the Colombian population living with DM who receive services within the framework of the health system. Furthermore, a longitudinal approach was used and the effect of repeated measures during the follow up was considered in the final models. This differs from the commonly used cross sectional analysis in studies evaluating factors influencing glycemic control. On the other hand, the use of HbA1c has been well described to evaluate long term glycemic control and because we only studied people with a new DM diagnosis, changes in glycemic control over time are meaningful for decision making process at a clinical and organizational level.

Some limitations should be considered when interpreting our findings. There is a clear effect of medication and non-pharmacological treatment on prospective glycemic control. Nonetheless, we were unable to adjust for treatment because these data were unavailable so residual confusion cannot be ruled out. Following the same line, although some sociodemographic factors such as educational level, income and urban/rural areas have been associated with PGC, we did not evaluate their effects because these variables are not available in the NRCKD. Moreover, the NRCKD is a passive registry and reporting process is performed by health care insurers; this can lead to underreporting and despite it is minimal, it must be declared. Linked to the above, data on variables are likely subjected to measurement error. Nevertheless, its impact was limited through the data monitoring process. In any case this error would be non-differential and effect estimations tend to the null.


In this adult population living with DM, patients with concomitant CKD or obesity had an increased likelihood of PGC. This supports the importance of weight reduction in diabetes treatment; as well as the increased surveillance and treatment needed in the CKD population. Our findings are based on data from real-world clinical practice thus they should be helpful for public health entities in decision making process. The strategies should be directed in improving access to healthcare, surveillance and treatment in these patients given the detrimental effects of PGC. Whether the presence of HTN modifies the probability of developing a PGC could not be clarified in this study; this effect was only modified in afro descendant ethnicity where HTN decreased the odds of a PGC. Further studies evaluating the effect of pharmacological and lifestyle interventions as well as sociodemographic factors such as educational level, socioeconomic status, or living area (rural versus urban) on the association of interest are needed. The effect of early-onset DM on the increased likelihood of having PGC is noteworthy; hence, further research that evaluates the long-term effect of optimal glycemic levels on the risk of DM-associated complications and mortality is needed.

Availability of data and materials

The datasets generated and/or analysed during the current study are not publicly available due to that they are owned by the Colombian health system but are available from the corresponding author on reasonable request.



Absolute Standard Differences


Body Mass Index


Chronic Kidney Disease


Chronic Kidney Disease Epidemiology Collaboration


Diabetes Mellitus


Glomerular Filtration Rate


Hemoglobin A1c




Interquartile range


National Registry of Chronic Kidney Disease


Poor Glycemic Control SD: Standard deviation


  1. Saeedi P, Petersohn I, Salpea P, Malanda B, Karuranga S, Unwin N, et al. Global and regional diabetes prevalence estimates for 2019 and projections for 2030 and 2045: Results from the International Diabetes Federation Diabetes Atlas, 9th edition. Diabetes Res Clin Pract. 2019;157:107843.

    Article  PubMed  Google Scholar 

  2. Vargas-Uricoechea H, Casas-Figueroa LÁ. An epidemiologic analysis of diabetes in Colombia. Ann Glob Health. 2015;81:742–53.

    Article  PubMed  Google Scholar 

  3. Harding JL, Pavkov ME, Magliano DJ, Shaw JE, Gregg EW. Global trends in diabetes complications: a review of current evidence. Diabetologia. 2019;62:3–16.

    Article  PubMed  Google Scholar 

  4. Herman WH, Braffett BH, Kuo S, Lee JM, Brandle M, Jacobson AM, et al. What are the clinical, quality-of-life, and cost consequences of 30 years of excellent vs. poor glycemic control in type 1 diabetes? J Diabetes Complicat. 2018;32:911–5.

    Article  Google Scholar 

  5. Zoungas S, Arima H, Gerstein HC, Holman RR, Woodward M, Reaven P, et al. Effects of intensive glucose control on microvascular outcomes in patients with type 2 diabetes: a meta-analysis of individual participant data from randomised controlled trials. Lancet Diabetes Endocrinol. 2017;5:431–7.

    Article  PubMed  Google Scholar 

  6. Ray KK, Seshasai SRK, Wijesuriya S, Sivakumaran R, Nethercott S, Preiss D, et al. Effect of intensive control of glucose on cardiovascular outcomes and death in patients with diabetes mellitus: a meta-analysis of randomised controlled trials. Lancet. 2009;373:1765–72.

    Article  CAS  PubMed  Google Scholar 

  7. Hayward RA, Reaven PD, Wiitala WL, Bahn GD, Reda DJ, Ge L, et al. Follow-up of glycemic control and cardiovascular outcomes in type 2 diabetes. N Engl J Med. 2015;372:2197–206.

    Article  CAS  PubMed  Google Scholar 

  8. Mannucci E, Monami M, Dicembrini I, Piselli A, Porta M. Achieving HbA1c targets in clinical trials and in the real world: a systematic review and meta-analysis. J Endocrinol Investig. 2014;37:477–95.

    Article  CAS  Google Scholar 

  9. De Pablos-Velasco P, Parhofer KG, Bradley C, Eschwège E, Gönder-Frederick L, Maheux P, et al. Current level of glycaemic control and its associated factors in patients with type 2 diabetes across Europe: Data from the PANORAMA study. Clin Endocrinol (Oxf). 2014;80:47–56.

    Article  CAS  Google Scholar 

  10. Smalls BL, Ritchwood TD, Bishu KG, Egede LE. Racial/ethnic differences in glycemic control in older adults with type 2 diabetes: United States 2003–2014. Int J Environ Res Public Health. 2020;17:950.

    Article  PubMed Central  Google Scholar 

  11. Suh DC, Choi IS, Plauschinat C, Kwon J, Baron M. Impact of comorbid conditions and race/ethnicity on glycemic control among the US population with type 2 diabetes, 1988-1994 to 1999-2004. J Diabetes Complicat. 2010;24:382–91.

    Article  Google Scholar 

  12. Chetoui A, Kaoutar K, Elmoussaoui S, Boutahar K, El Kardoudi A, Chigr F, et al. Prevalence and determinants of poor glycaemic control: a cross-sectional study among Moroccan type 2 diabetes patients. Int Health. 2020;00:1–8.

    Google Scholar 

  13. Iglay K, Hannachi H, Howie PJ, Xu J, Li X, Engel SS, et al. Prevalence and co-prevalence of comorbidities among patients with type 2 diabetes mellitus. Curr Med Res Opin. 2016;32:1243–52.

    Article  PubMed  Google Scholar 

  14. Luijks H, Biermans M, Bor H, van Weel C, Lagro-Janssen T, de Grauw W, et al. The effect of comorbidity on glycemic control and systolic blood pressure in type 2 diabetes: a cohort study with 5 year follow-up in primary care. PLoS One. 2015;10:e0138662.

    Article  PubMed  PubMed Central  CAS  Google Scholar 

  15. Bralić Lang V, Bergman MB. Prevalence of comorbidity in primary care patients with type 2 diabetes and its association with elevated HbA1c: a cross-sectional study in Croatia. Scand J Prim Health Care. 2016;34:66–72.

    Article  PubMed  PubMed Central  Google Scholar 

  16. Hudon C. Comorbidity and glycemia control among patients with type 2 diabetes in primary care. Diabetes Metab Syndr Obes Targets Ther. 2008;1:33–7.

    Article  CAS  Google Scholar 

  17. Almetwazi M, Alwhaibi M, Balkhi B, Almohaini H, Alturki H, Alhawassi T, et al. Factors associated with glycemic control in type 2 diabetic patients in Saudi Arabia. Saudi Pharm J. 2019;27:384–8.

    Article  PubMed  Google Scholar 

  18. Alzaheb R, Altemani A. The prevalence and determinants of poor glycemic control among adults with type 2 diabetes mellitus in Saudi Arabia. Diabetes Metab Syndr Obes Targets Ther. 2018;11:15–21.

    Article  Google Scholar 

  19. Ministerio de Salud y Proteccion Social - Republica de Colombia. Resolucion 2643 del 2014; 2014. p. 1–13. Accessed 9 Oct 2020

    Google Scholar 

  20. Ministerio de Salud y Proteccion Social - Republica de Colombia. Cifras de aseguramiento en salud con corte diciembre 2020. 2020. Accessed 19 Jan 2021.

    Google Scholar 

  21. Ministerio de Salud y Protección Social de Colombia. Guía de práctica clínica Hipertensión arterial primaria (HTA). 2013. Accessed 22 Apr 2021.

    Google Scholar 

  22. Ministerio de Salud y Protección Social de Colombia. Guía de Práctica Clínica para el manejo de la hipertensión arterial primaria (HTA). 2017; 2da edición. Guía No.18. Accessed 22 Apr 2021.

  23. Ministerio de Salud y Protección Social de Colombia. Guía de Práctica Clínica para el diagnóstico y tratamiento de la Enfermedad Renal Crónica (adopción). 2016. Accessed 22 Apr 2021.

    Google Scholar 

  24. WHO Consultation on Obesity (1999:Geneva Switzerland) & World Health Organization. Obesity: Preventing and managing the global epidemic: Report of a WHO Consultation (WHO Technical Report Series 894). 2000:1-252. Accessed 9 Oct 2020.

  25. American Diabetes Association. 6. Glycemic Targets: Standards of Medical Care in Diabetes-2020. Diabetes Care. 2020;43(Suppl. 1):S66–76.

    Article  Google Scholar 

  26. Departamento Administrativo Nacional de Estadistica (DANE). Cuentas departamentales, Producto Interno Bruto por departamento 2019 preliminar; 2020. p. 7–8. Accessed 10 Nov 2020

    Google Scholar 

  27. Levey AS, Stevens LA, Schmid CH, Zhang Y (Lucy), Castro AF, Feldman HI, et al. A new equation to estimate glomerular filtration rate. Ann Intern Med. 2009;150:604–12.

    Article  PubMed  PubMed Central  Google Scholar 

  28. Kidney Disease: Improving Global Outcomes (KDIGO) CKD Work Group. KDIGO 2012 clinical practice guideline for the evaluation and management of chronic kidney disease. Kidney Int Suppl. 2013;3:1–150.

    Google Scholar 

  29. Yang D, Dalton JE. A unified approach to measuring the effect size between two groups using SAS®. SAS Glob Forum Stat Data Anal. 2012;335:1–6.

    Google Scholar 

  30. Wilmot E, Idris I. Early onset type 2 diabetes: risk factors, clinical impact and management. Ther Adv Chronic Dis. 2014;5:234–44.

    Article  PubMed  PubMed Central  Google Scholar 

  31. Kidney Disease: Improving Global Outcomes (KDIGO) Diabetes Work Group. KDIGO 2020 clinical practice guideline for diabetes management in chronic kidney disease. Kidney Int. 2020;98:S1–115.

    Article  Google Scholar 

  32. Sosa-Rubí SG, Galárraga O, López-Ridaura R. Diabetes treatment and control: the effect of public health insurance for the poor in Mexico. Bull World Health Organ. 2009;87:512–9.

    Article  PubMed  PubMed Central  Google Scholar 

  33. Zhang X, McKeever Bullard K, Gregg EW, Beckles GL, Williams DE, Barker LE, et al. Access to health care and control of ABCS of diabetes. Diabetes Care. 2012;35:1566–71.

    Article  PubMed  PubMed Central  Google Scholar 

  34. Ministerio de Salud y Protección Social de Colombia. Todo lo que usted debe saber sobre el plan de beneficios POS. 2014. Accessed 7 May 2021.

    Google Scholar 

  35. Copur S, Onal EM, Afsar B, Ortiz A, van Raalte DH, Cherney DZ, et al. Diabetes mellitus in chronic kidney disease: biomarkers beyond HbA1c to estimate glycemic control and diabetes-dependent morbidity and mortality. J Diabetes Complicat. 2020;34:107707.

    Article  Google Scholar 

  36. De Cosmo S, Viazzi F, Pacilli A, Giorda C, Ceriello A, Gentile S, et al. Achievement of therapeutic targets in patients with diabetes and chronic kidney disease: insights from the Associazione Medici Diabetologi annals initiative. Nephrol Dial Transplant. 2015;30:1526–33.

    Article  PubMed  CAS  Google Scholar 

  37. Bae JP, Lage MJ, Mo D, Nelson DR, Hoogwerf BJ. Obesity and glycemic control in patients with diabetes mellitus: analysis of physician electronic health records in the US from 2009-2011. J Diabetes Complicat. 2016;30:212–20.

    Article  CAS  Google Scholar 

  38. Speeckaert M, Van Biesen W, Delanghe J, Slingerland R, Wiecek A, Heaf J, et al. Are there better alternatives than haemoglobin A1c to estimate glycaemic control in the chronic kidney disease population? Nephrol Dial Transplant. 2014;29:2167–77.

    Article  CAS  PubMed  Google Scholar 

  39. Hill CJ, Maxwell AP, Cardwell CR, Freedman BI, Tonelli M, Emoto M, et al. Glycated hemoglobin and risk of death in diabetic patients treated with hemodialysis: a meta-analysis. Am J Kidney Dis. 2014;63:84–94.

    Article  CAS  PubMed  Google Scholar 

  40. Kuo I-C, Lin HY-H, Niu S-W, Hwang D-Y, Lee J-J, Tsai J-C, et al. Glycated hemoglobin and outcomes in patients with advanced diabetic chronic kidney Disease. Sci Rep. 2016;6:20028.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  41. Esteghamati A, Ismail-Beigi F, Khaloo P, Moosaie F, Alemi H, Mansournia MA, et al. Determinants of glycemic control: phase 2 analysis from nationwide diabetes report of National Program for prevention and control of diabetes (NPPCD-2018). Prim Care Diabetes. 2020;14:222–31.

    Article  PubMed  Google Scholar 

  42. Gummesson A, Nyman E, Knutsson M, Karpefors M. Effect of weight reduction on glycated haemoglobin in weight loss trials in patients with type 2 diabetes. Diabetes Obes Metab. 2017;19:1295–305.

    Article  CAS  PubMed  Google Scholar 

  43. Lerman-Garber I, Aguilar-Salinas C, Tusié-Luna T, Velásquez D, Lobato-Valverde M, Osornio-Flores M, et al. Diabetes tipo 2 de inicio temprano. Experiencia en una institución de tercer nivel de atención médica. Gac Med Mex. 2010;146:179–84.

    PubMed  Google Scholar 

  44. Bo A, Thomsen RW, Nielsen JS, Nicolaisen SK, Beck-Nielsen H, Rungby J, et al. Early-onset type 2 diabetes: age gradient in clinical and behavioural risk factors in 5115 persons with newly diagnosed type 2 diabetes-results from the DD2 study. Diabetes Metab Res Rev. 2018;34:e2968.

    Article  CAS  Google Scholar 

  45. Gopalan A, Mishra P, Alexeeff SE, Blatchins MA, Kim E, Man A, et al. Initial glycemic control and care among younger adults diagnosed with type 2 diabetes. Diabetes Care. 2020;43:975–81.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  46. Badedi M, Solan Y, Darraj H, Sabai A, Mahfouz M, Alamodi S, et al. Factors associated with long-term control of type 2 diabetes mellitus. J Diabetes Res. 2016;2016:2109542.

    Article  PubMed  PubMed Central  CAS  Google Scholar 

  47. Menzin J, Korn JR, Cohen J, Lobo F, Zhang B, Friedman M, et al. Relationship between glycemic control and diabetes-related hospital costs in patients with type 1 or type 2 diabetes mellitus. J Manag Care Pharm. 2010;16:264–75.

    PubMed  Google Scholar 

  48. Gil J, Sicras-Mainar A, Zucchelli E. Uncontrolled diabetes and health care utilisation: panel data evidence from Spain. Eur J Health Econ. 2018;19:785–95.

    Article  PubMed  Google Scholar 

  49. Mata-Cases M, Rodríguez-Sánchez B, Mauricio D, Real J, Vlacho B, Franch-Nadal J, et al. The association between poor glycemic control and health care costs in people with diabetes: a population-based study. Diabetes Care. 2020;43:751–8.

    Article  PubMed  Google Scholar 

Download references


We would like to thank Dr. Oscar Leonel Rueda and Dr. Lyda Rojas for their recommendations and comments on the manuscript.


This research did not receive any sources of funding.

Author information

Authors and Affiliations



MUJ, LJHP, JAHV, AMVG, LAM, MUT were involved in the conception and design of the study. LJHP and JAHV conducted the statistical analyses and tables/figures. AMVG supervised the statistical analyses. MUJ, LJHP, JAHV, AMVG, LAM, MUT contributed to the interpretation of the analyses. MUJ and JAHV wrote the initial draft of the paper. MUJ, LJHP, JAHV, AMVG, LAM, MUT contributed to the revision of the manuscript. All authors read and approved the final manuscript.

Corresponding author

Correspondence to Miguel Urina-Triana.

Ethics declarations

Ethics approval and consent to participate

All the methods were carried out following national (Resolution 8430 of 1993, stated by The Colombian Health Ministry) and international (The Declaration of Helsinki) standards. Approval by the Ethics Committee of Fundación del Caribe para la Investigación Biomédica (Record 0211 from October 17th, 2020) was obtained. According to the national regulations for conducting human research (Resolution 8430 of 1993, stated by The Colombian Health Ministry), this investigation has no risk for participants. Our research was based on secondary data sources which was anonymized and did not include any private information that could make any subject identifiable. Taking into account the above and the retrospective nature of our research, the Ethics Committee of Fundación del Caribe para la Investigación Biomédica waived the need for informed consent.

Consent for publication

Not applicable.

Competing interests

The authors declare that they have no competing interests.

Additional information

Publisher’s Note

Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Supplementary Information

Additional file 1. Table S1.

Distribution of Colombian departments by region*. 

Additional file 2. Table S2.

Multivariate-adjusted odds of poor glycemic control according to sensitivity analysis criteria, Colombia 2014-2019.

Rights and permissions

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 The Creative Commons Public Domain Dedication waiver ( applies to the data made available in this article, unless otherwise stated in a credit line to the data.

Reprints and permissions

About this article

Check for updates. Verify currency and authenticity via CrossMark

Cite this article

Urina-Jassir, M., Herrera-Parra, L.J., Hernández Vargas, J.A. et al. The effect of comorbidities on glycemic control among Colombian adults with diabetes mellitus: a longitudinal approach with real-world data. BMC Endocr Disord 21, 128 (2021).

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI: