Skip to main content

ORIGINAL RESEARCH article

Front. Endocrinol., 30 June 2021
Sec. Clinical Diabetes

Clinical Characteristics and Risk of Diabetic Complications in Data-Driven Clusters Among Type 2 Diabetes

Lin Xing&#x;Lin Xing1†Fangyu Peng&#x;Fangyu Peng1†Qian Liang&#x;Qian Liang2†Xiaoshuang DaiXiaoshuang Dai1Junli RenJunli Ren1Han WuHan Wu2Shufen YangShufen Yang2Yaxin ZhuYaxin Zhu2Lijing Jia*Lijing Jia2*Shancen Zhao*Shancen Zhao1*
  • 1BGI Institute of Applied Agriculture, BGI-Shenzhen, Shenzhen, China
  • 2Endocrinology Department, Shenzhen People’s Hospital, Shenzhen, China

Background: This study aimed to cluster newly diagnosed patients and patients with long-term diabetes and to explore the clinical characteristics, risk of diabetes complications, and medication treatment related to each cluster.

Research Design and Methods: K-means clustering analysis was performed on 1,060 Chinese patients with type 2 diabetes based on five variables (HbA1c, age at diagnosis, BMI, HOMA2-IR, and HOMA2-B). The clinical features, risk of diabetic complications, and the utilization of elven types of medications agents related to each cluster were evaluated with the chi-square test and the Tukey–Kramer method.

Results: Four replicable clusters were identified, severe insulin-resistant diabetes (SIRD), severe insulin-deficient diabetes (SIDD), mild obesity-related diabetes (MOD), and mild age-related diabetes (MARD). In terms of clinical characteristics, there were significant differences in blood pressure, renal function, and lipids among clusters. Furthermore, individuals in SIRD had the highest prevalence of stages 2 and 3 chronic kidney disease (CKD) (57%) and diabetic peripheral neuropathy (DPN) (67%), while individuals in SIDD had the highest risk of diabetic retinopathy (32%), albuminuria (31%) and lower extremity arterial disease (LEAD) (13%). Additionally, the difference in medication treatment of clusters were observed in metformin (p = 0.012), α-glucosidase inhibitor (AGI) (p = 0.006), dipeptidyl peptidase 4 inhibitor (DPP-4) (p = 0.017), glucagon-like peptide-1 (GLP-1) (p <0.001), insulin (p <0.001), and statins (p = 0.006).

Conclusions: The newly diagnosed patients and patients with long-term diabetes can be consistently clustered into featured clusters. Each cluster had significantly different patient characteristics, risk of diabetic complications, and medication treatment.

Introduction

Diabetes is a chronic disease, not only has caused heavy social and economic burdens, but also prone to leading to multiple complications, which have profound impacts on the life quality of patients and may potentially cause death in severe cases. The prevalence of diabetes is rapidly increasing worldwide, so effectively preventing and managing diabetes has become an important topic at this stage (1, 2).

Diabetes is characterized by hyperglycemia, the causes of which are highly heterogeneous (3). Based on current classification criteria, diabetes is currently divided into two major subtypes, type 1 diabetes (T1D) and type 2 diabetes (T2D) which is approximately 85% (2, 4, 5). This classification relies on the age of disease diagnosis; however, it may not be enough to characterize complications and outcomes for subtypes. Individuals with diabetes have a different natural course of hyperglycemia, therefore theoretically should be treated with different clinical strategies better fitted their metabolic characteristics (6, 7). A novel approach to detailed characterize the diabetes population and explore the clinical features can be very beneficial to aid with the treatment of diabetes patients.

In recent years, novel stratifications of diabetes have been attempted worldwide. Three subgroups of T2D were identified using a topological analysis based on patient–patient networks (8). It is a valuable attempt to classify the patients, however, because the approach required genotype data from patients, this can be difficult to implement in clinical settings. Moreover, in Ahlqvist and colleagues’ study, five replicable clusters of diabetes based on six common clinical variables were found, which included glutamic acid decarboxylase autoantibody (GADA), HbA1c, BMI, age at onset of diabetes, and homeostasis model estimates of β-cell function (HOMA2B) and insulin resistance (HOMA2IR) (9). The five diabetes clusters were cluster 1, severe autoimmune diabetes (SAID); cluster 2, severe insulin-deficient diabetes (SIDD); cluster 3, severe insulin-resistant diabetes (SIRD); cluster 4, mild obesity-related diabetes (MOD); and cluster 5, mild age-related diabetes (MARD). Cluster 1 was characterized by the presence of GADA, being similar to T1D, and the other four clusters were T2D with the absence of GADA positivity and the other five variables. Safai et al. (10) used similar routine clinical markers to sub-group the patients into five clusters and reported the difference in probability of diabetes complications, such as cardiovascular disease, nephropathy and neuropathy. In Zaharia et al. (11), patients in German with newly diagnosed type 1 or type 2 diabetes were grouped into the same five clusters, validating such techniques in different population. In addition, the five-year follow-up research also reported the different prevalence in subgroups in terms of non-alcoholic fatty liver disease and diabetic neuropathy. In Dennis et al. (12), the five clusters could be replicated and the difference in glycemic progression were identified. Ahlqvist and colleagues in 2020 revisited the sub-grouping technique in several populations and reported the difference among groups for diabetic complications, such as retinopathy, neuropathy, kidney disease and fatty liver (13). To alleviate the requirement on clinical markers, in Kahkoska et al. (14), three global trails (DEVOTE, LEADER, and SUSTAIN-6) with recent-onset diabetes were tested for the clustering technique based on three variables, age at diabetes diagnosis, baseline glycated hemoglobin (HbA1c) and body mass index (BMI). The four T2D clusters could be fully replicated and the risk of major adverse cardiovascular events and death differed significantly in the follow-up duration. In addition, the risk of nephropathy differed in clusters. With clinical variables, HbA1c, BMI, age at onset of diabetes, HOMA2B and HOMA2IR, there were attempts that clustering the newly diagnosed diabetes patients in the United States and China into four subgroups and confirmed that different ethnic groups can also be clustered by the same variables (15). For many previous researches that had already proven the promising application of such sub-grouping techniques in precision medicine, the clustering target were usually newly diagnosed patients. Since the applications of sub-grouping technique on different populations were proven to be quite robust, we would like to make some exploration of the clustering techniques on both newly diagnosed and long-term diabetes patients, aiming to facilitating the application domain as many admitted patients were in various stages.

Current diabetes treatments mainly focus on controlling blood glucose levels, but with many researches elucidated specific characteristics with subgroups of T2D, more informative treatments on diabetes related issues, such as kidney, cardiovascular and cerebrovascular diseases, became promising. We believe a precise characterization of T2D patient populations in the clinical settings would be beneficial for the understanding of T2D pathophysiology and improvement of clinical management. Therefore, the objectives of this study are to cluster newly diagnosed and previously known Chinese T2D patients, and to explore the clinical characteristics, the risks of diabetic complications, and medication treatment in each cluster. We aim to identify different subgroups of T2D patients through the K-means clustering method with five commonly used clinical variables including HbA1c, BMI, age at diagnosis, HOMA2-B, and HOMA2-IR, and then compare clinical characteristics, identify individuals with increased risk of complications, and distinguish different medication treatment in each subgroup.

Materials and Methods

Study Population

This was a cross-sectional study conducted from January 2018 to November 2019 at the No.1 Shenzhen People’s Hospital. Medical records of 1,240 participants were collected on a first come first chosen basis within a one-year time window, which include anthropometric measurements, laboratory tests, complication diagnostic information, and medication regime. All participants, diagnosed with type 2 diabetes, involved in this study were aged 18 years and above. This study was approved by the medical research ethics committee at Shenzhen People’s Hospital. Informed consent was obtained from the participants subjected to anonymous information utilization in medical research.

Measurements

The height and weights of participants were measured using an automatic anthropometer and the body mass index (BMI) was calculated as body weight/height (kg/m2). Blood pressure was measured by trained nurses with a blood pressure monitor. Laboratory measurements were taken in a fasting state following the standardized procedures during the health examination. Biochemical indices, such as fasting blood glucose (FPG), urine acid (UA), total cholesterol (TC), triglycerides (TG), high-density lipoprotein (HDL) cholesterol and low-density lipoprotein (LDL) cholesterol, were measured by the hexokinase method and C-peptide (CP) concentrations were measured by radioimmunoassay. The FPG and CP were used to calculate homeostasis model assessment 2 estimates of insulin resistance (HOMA2IR) and homeostasis model assessment 2 estimates of β-cell function (HOMA2B) with the HOMA2 calculator v2.2.3 at www.dtu.ox.ac.uk (16).

Definitions of Diabetes and Diabetic Complications

The criteria to diagnose participants to have diabetes followed the internationally adopted standards set by the World Health Organization (WHO) Diabetes Expert Committee in 1999 (17). Estimated glomerular filtration rate (eGFR) was calculated using the chronic kidney disease epidemiology collaboration (CKD-EPI) equation, which was used to classify kidney function as normal (stage 1, eGFR >90 ml/min per 1.73 m2), abnormal (stage 2, eGFR 60–90 ml/min per 1.73 m2 or stage 3, eGFR <60 ml/min per 1.73 m2) (18, 19). The range of urinary albumin to creatinine ratio (UACR) was used to describe the Albuminuria progression. The UACR less than 30, between 30 and 300, and above 300, were defined as normal, microalbuminuria, and macroalbuminuria, respectively. Ankle-brachial index (ABI) measurements were taken after a five-minute break with the supine position, which were used to identify the lower extremity arterial disease (LEAD) including hardened vessels and arterial occlusion (20). The ABI values between 0.9 and 1.3, larger than 1.3 or less than 0.9 were considered as normal and abnormal arterial, respectively (21). Both diabetic retinopathy and diabetic peripheral neuropathy (DPN) were defined following the American Diabetes Association’s criteria (22, 23). The diagnosis of DPN was based on the multiple symptoms and diabetes history. For the impact on the small fibers, the symptoms usually involved pain and dysesthesia, which can be assessed with the pinprick and temperature sensation tests. As the impact developed on large fibers, the symptoms usually involved numbness and loss of protective sensation, which can be confirmed by the vibration perception and 10-g monofilament tests. For diabetic retinopathy, the diagnosis was based on an initial dilated and comprehensive eye examination performed by an ophthalmologist or optometrist. Patients were arranged for diagnosis examinations during the first visit, based on which the initial status of diabetic retinopathy was determined (22, 23).

Cluster Analysis

The data cleaning process followed four steps. Step-1, twenty-two participants without type 2 diabetes were removed from the dataset as they were not part of the target group. Step-2, twenty-three individuals with missing information in variables, such as BMI, HbA1c, and so on, were removed. Step-3, To improve the group clustering quality, one hundred thirty-four participants with values beyond the defined range of HOMA2 calculators were removed before the feature engineering. Step-4, by checking the variables individually, one extreme outlier was identified and removed. After the data cleaning procedures, 1,060 participants were included in cluster analysis, with the result of which two additional analyses focused on medication usage and diabetes complications were carried out. For the comparison of diabetes complications, 1,060 participants were included in the analysis. For the medication usage comparison, after the removal of participants that did not have the required information, 486 participants were included in medication analysis (Figure 1).

FIGURE 1
www.frontiersin.org

Figure 1 Flowchart for the selection of participants in the different analysis section.

Statistical Analysis

The k-means was used to cluster the data according to five variables, age at diagnosis, BMI, HbA1c, HOMA2B, and HOMA2IR. All data were scaled to mean zero and unit variance before clustering. K-means clustering was performed using the Hartigan and Wong algorithm implemented in the ‘k-means’ package and the optimal number of clusters was determined by elbow method from ‘NbClust’ package in R. To use the elbow method in the given data set, the number of clusters were plotted against the total explained variance. With a straight line going across the different number of clusters, the point where the increase of variance explained became slow was the target. The number of clusters corresponding to the point was the reasonable number of clusters. Based on the characteristics of clusters described by Ahlqvist et al. (9), patients were assigned to explainable clusters, severe insulin-deficient diabetes (SIDD), severe insulin-resistant diabetes (SIRD), mild obesity-related diabetes (MOD), or mild age-related diabetes (MARD). For descriptive statistics between subgroups, the chi-square test was used for categorical data. Skewed data were log-transformed before analysis. To account for the impact of age, the linear model incorporated age variable as correction factor when comparing the clinical features in subgroups. To understand the impact of gender on the diabetes-related complication of subgroups, odds ratios (ORs) were modeled using logistic regression. P-values less than 0.05 were considered statistically significant. Statistical analyses were done with R version 3.5.3.

Results

Cluster Analysis

General characteristics of the population and clinical data for these patients are shown in Table 1. In this study, the age of participants ranged from 24 to 99, with most individuals being male (61%). According to the onset of the patient’s diabetes, around 16% of the patients were newly diagnosed diabetes, and 84% of the patients were long-term diabetes (Figure S1).

TABLE 1
www.frontiersin.org

Table 1 Characteristics of the participants in allocated clusters.

One thousand sixty patients were classified into four clusters, each of which had distinctive clinical features (Figure 2 and S2). Cluster 1, including 21% of the patients characterized by high age, low BMI, low insulin secretion (low HOMA2B index), and poor metabolic control (the highest HbA1c and Fasting blood glucose level), was identified as severe insulin-deficient diabetes (SIDD). Cluster 2 included 21% of patients who were labeled as severe insulin-resistant diabetes (SIRD), which was characterized by insulin resistance (high level of HOMA2-B and HOMA2-IR), high BMI, and good metabolic control. Cluster 3, including 25% of patients, was labeled as mild obesity-related diabetes (MOD) characterized by obesity, low age, average β-cell function, and insulin resistance. Cluster 4, the largest subgroup (33%), was labeled as mild age-related diabetes (MARD), which include older participants with modest metabolic control, insulin resistance, β-cell function, and the lowest BMI.

FIGURE 2
www.frontiersin.org

Figure 2 Characteristics of clusters and distribution of participants. (A–E) The distribution of HbA1c, BMI, Age, HOMA2B, and HOMA2IR in SIDD, SIRD, MOD, and MARD clusters. (F) Distribution of patients (n = 1060) in subgroups clustered by k-means approach. SIDD, severe insulin-deficient diabetes; SIRD, severe insulin-resistant diabetes; MOD, mild obesity-related diabetes; MARD, mild age-related diabetes; HOMA2B, homeostatic model assessment 2 estimates of β-cell function; HOMA2IR, homeostatic model assessment 2 estimates of insulin resistance.

Comparison of Clinical Features in Subgroups

The clinical characteristics of subgroups can be found in Table 1. In this section, the study mainly focused on the characteristics of subgroups in blood pressure, renal function, and lipids. Blood pressure, systolic blood pressure (SBP), and diastolic blood pressure (DBP) were following a similar trend, being highest in SIDD subgroup. Besides, in terms of DBP levels, there was no significant difference between SIRD and MARD (p = 0.104) (Tables S1, S2).

Regarding renal function, patients in SIRD had the lowest eGFR and highest urine acid (UA) compare to those in other subgroups. Urinary albumin to creatinine ratio (UACR) was highest in MOD cluster whereas there no statistical significance in MARD compared to SIDD (p = 0.174) and SIRD (p = 0.334), respectively (Table S5).

Concerning lipids, total cholesterol (TC) was the highest in patients assigned to the SIDD cluster and being statistically higher than MARD (p <0.001) whereas there was no significant difference between MOD and MARD (Table S6). Patients in MARD subgroup had the lowest triglycerides (TG) and highest HDL-cholesterol (HDL) levels comparable to those in other subgroups (Tables S7–S8). LDL-cholesterol (LDL) was the lowest in patients assigned to the SIRD cluster compared with all other clusters and there was no statistical significance between MOD and MARD clusters (p = 0.373) (Table S9).

Comparison of Diabetes Complications in Subgroups

The differences in the prevalence of diabetes complications between the subgroups were shown in Figure 3. Patients in SIRD had the highest prevalence of stages 2 and 3 CKD in the SIRD cluster (57%) compared to those in SIDD (41%; p = 0.001), MOD (29%; p <0.001) and MARD (43%; p= 0.002) (Figure 3A). The SIDD cluster, being similar to MOD (p = 0.172), had the highest prevalence of albuminuria (31%), compared to those in SIRD (18%; p = 0.002) and MARD (14%; p <0.001) (Figure 3B). Patients in the MOD (48%) cluster had lowest prevalence of DPN (Figure 3C) compared to those in SIDD (60%; p = 0.02), SIRD (67%, p = 0.001) and MARD (63%; p = 0.003). The prevalence of retinopathy was highest in SIDD (32%) and MOD (32%) (Figure 3D), but there was no significant difference between the rest groups (p >0.05). The LEAD had the highest prevalence in SIDD (13%), but not being significantly different from SIRD (8%; p = 0.209) and MARD (8%; p = 0.168) (Figure 3E).

FIGURE 3
www.frontiersin.org

Figure 3 Characteristics of diabetes-related complications in subgroups. (A) chronic kidney disease; (B) albuminuria; (C) diabetic peripheral neuropathy; (D) diabetic retinopathy; (E) Lower-extremity arterial disease diabetes-related complications. SIDD, severe insulin-deficient diabetes; SIRD, severe insulin-resistant diabetes; MOD, mild obesity-related diabetes; MARD, mild age-related diabetes.

Separate analyses stratified by gender were performed to illustrate the risk trends of diabetes-related complications in subgroups. For male participants, there was no significantly different (p >0.05) in albuminuria and retinopathy among subgroups (Table 2). In other complications, the MARD consistently had a higher risk of disease onset compared to other subgroups (ORs >1; p <0.05). For female participants, the risk of diabetes complication, retinopathy, was similar among subgroups (p >0.05), whereas, in other subgroups, the risk trends varied (Table 2).

TABLE 2
www.frontiersin.org

Table 2 Association between diabetes-related complications and clusters in different gender stratification.

Comparison of Medication Application in Subgroups

The detailed comparison of nine types of glucose-lowering drugs, statins, and antihypertensive in different subgroups were shown in Table 3. Significant differences in treatment between the subgroups were observed in metformin (p = 0.012), α-glucosidase inhibitor (AGI) (p = 0.006), dipeptidyl peptidase 4 inhibitor (DPP-4) (p = 0.017), glucagon-like peptide-1 (GLP-1) (p <0.001), insulin (p <0.001), statins (p = 0.006).

TABLE 3
www.frontiersin.org

Table 3 Characteristics of medication application by allocated clusters.

Patients in the SIDD cluster had the highest proportion of Insulin (76%), AGI (57%), and statins (81%) treatments compared to those in other subgroups. The MOD cluster had the highest proportion of metformin treatment (86%), compared with SIDD (72%), SIRD (74%), and MARD (69%). Patients in MOD cluster also had the highest percentage of GLP-1 treatment (16%), compared with those in SIDD (1%), SAID (4%), and MARD (1%). Moreover, DPP-4 treatment was most frequently used in MARD subgroup (64%).

Discussion

In this study, the data-driven approach to distinguish the status of diabetes was reproducible and the distribution of patients was similar to that of the Swedish cohort (9). To improve the comprehensive application of such clustering approach in clinical settings, this study tested on population included both newly diagnosed patients and patients with long-term diabetes, which is different from previous studies that only covered newly diagnosed individuals (9, 11, 12, 15). Through the exploration on mixed population, we hope to extend the approach closer to clinical practices as many admitted patients were not newly diagnosed. To evaluate the performance of the clustering technique on data set, we checked the mathematical clustering stability and interpretation of the grouped clusters. Overall, the data driven approach could be used on the mixed population and yielded decent stability. Notably, the clinical characteristics of clusters in this study were similar to previous researches, with a slightly higher proportion of patients in the severe insulin-deficient diabetes (SIDD) cluster (21%) and severe insulin-resistant diabetes (SIRD) cluster (21%). This discrepancy may result from the higher proportion of hospitalized and severe participants in our study than other research population.

The highest prevalence of stages 2 and 3 chronic kidney disease (CKD) was in the SIRD subgroup with the lowest eGFR and highest urine acid. The possible explanation may be that the significant feature of the SIRD group is insulin resistance, which could lead to water-sodium retention, glomerular hypertension, hyperfiltration, and hyperuricemic, thus accelerating the progression of CKD. The trend of risk levels associated with CKD was different in male and female gender (Table 2). There were more subgroups associated with increased risk of CKD in males compared to females, which leads to a possible belief that the prognosis of CKD may depend on gender. This was confirmed by research showing that the female gender was associated with a slower decline in GFR and better patient and renal survival in a 10-year following-up study (24). In 2016, summarized multiple publication papers related to CKD research, Goldberg and Krause reported that the mortality risk of CKD in males was higher than in females (25).

Diabetic retinopathy, presented in approximately 30% of the individuals with diabetes (26), occurred most frequently in SIDD (32%) and mild obesity-related diabetes (MOD) (32%). SIDD, together with MOD, had the highest levels of fasting blood glucose and blood pressure (BP). Multiple studies showed that hyperglycemia and hypertension are risk factors for retinopathy and were different explanations regarding the connections. One possible reason is metabolic pathways triggered by hyperglycemia in diabetes, such as the polyol and the hexosamine pathways, the de novo synthesis of diacyl-glycerol, and advanced glycosylation end products (AGEs), can promote the development of retinopathy (27). Another possible reason is that hypertension increases the expression of pro-inflammatory molecules in the retina (28). The albuminuria, which is a powerful predictor of renal and cardiovascular risk (29), especially microalbuminuria, was highest in the SIDD subgroup. Since hypertension and hypercholesterolemia are causal risk factors of cardiovascular diseases as well (30), it is no surprise to have the SIDD group with a high rate of cardiovascular diseases. Based on the literature, retinopathy can often precede diabetic nephropathy in patients with T2D and this was confirmed by the clustering result that SIDD subgroup has the highest prevalence on both albuminuria and retinopathy. In terms of gender effect, there was a study indicating that the female showed a significantly higher prevalence and the female gender was an independent factor of disease development (31).

The prevalence of lower extremity arterial disease (LEAD) was highest in the SIDD subgroup (13%). Moreover, the retinopathy and albuminuria also had the highest prevalence in SIDD and there was literature indicating that they may be independent risk factors for LEAD (3234). Similar to other researches, this study found the blood pressure and blood lipids were key risk factors for LEAD, being highest in SIDD (34, 35).

Regarding medication strategy, subgroups had different treatments, being consistent with the physiological characteristics of patients. The proportion of insulin, AGI, and statins was significantly higher in SIDD compared with other subgroups. On the other hand, individuals in SIDD hold very low β-cell reserves (the level of HOMA2-B was lowest), thus were treated preferentially with insulin. AGI drugs can reduce the amounts of insulin needed to control postprandial hyperglycemia by slowing down the digestion of complex carbohydrates and sucrose, therefore it is suitable for the insulin-deficient characteristic (36). The possible reason for the high prevalence application of statins in SIDD may result from the treatment of prevailed hyperlipidemia conditions.

In terms of MOD, having the highest BMI and mild diabetic symptoms, metformin, and GLP-1 were the medications used most frequently. Studies have shown that weight loss can effectively control diabetes’s disease course (37). As metformin and GLP-1 have significant effects on weight loss, the medications are suitable for MOD. DDP-4 is characterized as a low risk of hypoglycemia, a high compliance rate, therefore it is suitable for age-related diabetes (MARD) (38). The SIRD subgroup, characterized by the most severe insulin resistance and best β-cell function, was treated most frequently with metformin drugs. As metformin could increase insulin sensitivity, the application of insulin in this subgroup was lowest among all subgroups (39). Overall, the utilization of metformin, as the first-line treatment of type 2 diabetes, had the highest proportion in four subgroups, all above 69.4%.

To extend our study and comprehensively aid with clinical treatment on diabetes, long-term follow-up studies are necessary to explore disease progression and the treatment response. Additionally, newly discovered biomarkers with advanced techniques including genomics, transcriptomics, gut microbiota, could be considered to refine the current stratification strategy. To make the approach being feasible to apply in clinical conditions, a decision-making support system is necessary.

In conclusion, this study testified that the data-driven approach to cluster both newly diagnosed patients and patients with long-term diabetes can yield a consistent result. Each cluster had significantly different patient characteristics, risk of diabetic complications, and medication treatment. These findings have a potential value for clinical trial enrolment and early treatment stratification.

Data Availability Statement

The raw data supporting the conclusions of this article will be made available by the authors, without undue reservation.

Ethics Statement

The studies involving human participants were reviewed and approved by the Medical Research Ethics Committee at Shenzhen People’s Hospital. The patients/participants provided their written informed consent to participate in this study.

Author Contributions

LX and FP performed data analysis, interpretation and manuscript writing.QL supervised the data collection and research collaboration. XD, JR, HW and SY participated in data collection and literature search. LJ and SZ designed the experiment and supervised the overall progress. All authors contributed to the article and approved the submitted version.

Conflict of Interest

Authors LX, FP, XD, JR, and SZ were employed by company BGI-Shenzhen.

The remaining authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.

Supplementary Material

The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fendo.2021.617628/full#supplementary-material

References

1. Ogurtsova K, da Rocha Fernandes J, Huang Y, Linnenkamp U, Guariguata L, Cho NH, et al. IDF Diabetes Atlas: Global Estimates for the Prevalence of Diabetes for 2015 and 2040. Diabetes Res Clin Pract (2017) 128:40–50. doi: 10.1016/j.diabres.2017.03.024

PubMed Abstract | CrossRef Full Text | Google Scholar

2. Sahay M, Kalra S, Badani R, Bantwal G, Bhoraskar A, Das A, et al. Diabetes and Anemia: International Diabetes Federation (IDF)–Southeast Asian Region (SEAR) Position Statement. Diabetes Metab Syndrome: Clin Res Rev (2017) 11:S685–95. doi: 10.1016/j.dsx.2017.04.026

CrossRef Full Text | Google Scholar

3. Himsworth HP. Diabetes Mellitus. Its Differentiation Into Insolin-Sensitive and Insulin-Insensitive Types. Lancet (1936) 230:127–30. doi: 10.1016/S0140-6736(01)36134-2

CrossRef Full Text | Google Scholar

4. Association AD. 2. Classification and Diagnosis of Diabetes: Standards of Medical Care in Diabetes—2020. Diabetes Care (2020) 43:S14–31. doi: 10.2337/dc20-S002

PubMed Abstract | CrossRef Full Text | Google Scholar

5. N.D.D. Group. Classification and Diagnosis of Diabetes Mellitus and Other Categories of Glucose Intolerance. Diabetes (1979) 28:1039–57. doi: 10.2337/diab.28.12.1039

PubMed Abstract | CrossRef Full Text | Google Scholar

6. Gloyn AL, Drucker DJ. Precision Medicine in the Management of Type 2 Diabetes. Lancet Diabetes Endocrinol (2018) 6:891–900. doi: 10.1016/S2213-8587(18)30052-4

PubMed Abstract | CrossRef Full Text | Google Scholar

7. Prasad RB, Groop L. Precision Medicine in Type 2 Diabetes. J Internal Med (2019) 285:40–8. doi: 10.1111/joim.12859

PubMed Abstract | CrossRef Full Text | Google Scholar

8. Li L, Cheng W-Y, Glicksberg BS, Gottesman O, Tamler R, Chen R, et al. Identification of Type 2 Diabetes Subgroups Through Topological Analysis of Patient Similarity. Sci Trans Med (2015) 7:311ra174–311ra174. doi: 10.1126/scitranslmed.aaa9364

CrossRef Full Text | Google Scholar

9. Ahlqvist E, Storm P, Käräjämäki A, Martinell M, Dorkhan M, Carlsson A, et al. Novel Subgroups of Adult-Onset Diabetes and Their Association With Outcomes: A Data-Driven Cluster Analysis of Six Variables. Lancet Diabetes Endocrinol (2018) 6:361–9. doi: 10.1016/S2213-8587(18)30051-2

PubMed Abstract | CrossRef Full Text | Google Scholar

10. Safai N, Ali A, Rossing P, Ridderstråle M. Stratification of Type 2 Diabetes Based on Routine Clinical Markers. Diabetes Res Clin Pract (2018) 141:275–83. doi: 10.1016/j.diabres.2018.05.014

PubMed Abstract | CrossRef Full Text | Google Scholar

11. Zaharia OP, Strassburger K, Strom A, Bönhof GJ, Karusheva Y, Antoniou S, et al. Risk of Diabetes-Associated Diseases in Subgroups of Patients With Recent-Onset Diabetes: A 5-Year Follow-Up Study. Lancet Diabetes Endocrinol (2019) 7:684–94. doi: 10.1016/S2213-8587(19)30187-1

PubMed Abstract | CrossRef Full Text | Google Scholar

12. Dennis JM, Shields BM, Henley WE, Jones AG, Hattersley AT. Disease Progression and Treatment Response in Data-Driven Subgroups of Type 2 Diabetes Compared With Models Based on Simple Clinical Features: An Analysis Using Clinical Trial Data. Lancet Diabetes Endocrinol (2019) 7:442–51. doi: 10.1016/S2213-8587(19)30087-7

PubMed Abstract | CrossRef Full Text | Google Scholar

13. Ahlqvist E, Prasad RB, Groop L. Subtypes of Type 2 Diabetes Determined From Clinical Parameters. Diabetes (2020) 69:2086–93. doi: 10.2337/dbi20-0001

PubMed Abstract | CrossRef Full Text | Google Scholar

14. Kahkoska AR, Geybels MS, Klein KR, Kreiner FF, Marx N, Nauck MA, et al. Validation of Distinct Type 2 Diabetes Clusters and Their Association With Diabetes Complications in the DEVOTE, LEADER and SUSTAIN-6 Cardiovascular Outcomes Trials. Diabetes Obes Metab (2020) 22:1537–47. doi: 10.1111/dom.14063

PubMed Abstract | CrossRef Full Text | Google Scholar

15. Zou X, Zhou X, Zhu Z, Ji L. Novel Subgroups of Patients With Adult-Onset Diabetes in Chinese and US Populations. Lancet Diabetes Endocrinol (2019) 7:9–11. doi: 10.1016/S2213-8587(18)30316-4

PubMed Abstract | CrossRef Full Text | Google Scholar

16. Levy JC, Matthews DR, Hermans MP. Correct Homeostasis Model Assessment (HOMA) Evaluation Uses the Computer Program. Diabetes Care (1998) 21:2191–2. doi: 10.2337/diacare.21.12.2191

PubMed Abstract | CrossRef Full Text | Google Scholar

17. Alberti KGMM, Zimmet PZ. Definition, Diagnosis and Classification of Diabetes Mellitus and its Complications. Part 1: Diagnosis and Classification of Diabetes Mellitus. Provisional Report of a WHO Consultation. Diabetic Med (1998) 15:539–53. doi: 10.1002/(SICI)1096-9136(199807)15:7<539::AID-DIA668>3.0.CO;2-S

PubMed Abstract | CrossRef Full Text | Google Scholar

18. Levey AS, Coresh J, Greene T, Marsh J, Stevens LA, Kusek JW, et al. Expressing the Modification of Diet in Renal Disease Study Equation for Estimating Glomerular Filtration Rate With Standardized Serum Creatinine Values. Clin Chem (2007) 53:766–72. doi: 10.1373/clinchem.2006.077180

PubMed Abstract | CrossRef Full Text | Google Scholar

19. Levey AS, Stevens LA, Schmid CH, Zhang Y, Castro AF III, Feldman HI, et al. A New Equation to Estimate Glomerular Filtration Rate. Ann Internal Med (2009) 150:604–12. doi: 10.7326/0003-4819-150-9-200905050-00006

CrossRef Full Text | Google Scholar

20. Nativel M, Potier L, Alexandre L, Baillet-Blanco L, Ducasse E, Velho G, et al. Lower Extremity Arterial Disease in Patients With Diabetes: A Contemporary Narrative Review. Cardiovasc Diabetol (2018) 17:1–14. doi: 10.1186/s12933-018-0781-1

PubMed Abstract | CrossRef Full Text | Google Scholar

21. Chevtchouk L, da Silva MHS, do Nascimento OJM. Ankle-Brachial Index and Diabetic Neuropathy: Study of 225 Patients. Arquivos Neuro-psiquiatria (2017) 75:533–8. doi: 10.1590/0004-282x20170084

CrossRef Full Text | Google Scholar

22. Wilkinson CP, Ferris FL III, Klein RE, Lee PP, Agardh CD, Davis M, et al. Proposed International Clinical Diabetic Retinopathy and Diabetic Macular Edema Disease Severity Scales. Ophthalmology (2003) 110:1677–82. doi: 10.1016/S0161-6420(03)00475-5

PubMed Abstract | CrossRef Full Text | Google Scholar

23. Pop-Busui R, Boulton AJ, Feldman EL, Bril V, Freeman R, Malik RA, et al. Diabetic Neuropathy: A Position Statement by the American Diabetes Association. Diabetes Care (2017) 40:136–54. doi: 10.2337/dc16-2042

PubMed Abstract | CrossRef Full Text | Google Scholar

24. Eriksen B, Ingebretsen O. The Progression of Chronic Kidney Disease: A 10-Year Population-Based Study of the Effects of Gender and Age. Kidney Int (2006) 69:375–82. doi: 10.1038/sj.ki.5000058

PubMed Abstract | CrossRef Full Text | Google Scholar

25. Goldberg I, Krause I. The Role of Gender in Chronic Kidney Disease. EMJ (2016) 1:58–64.

Google Scholar

26. Klein R, Klein BE, Moss SE, Davis MD, DeMets DL. The Wisconsin Epidemiologic Study of Diabetic Retinopathy: II. Prevalence and Risk of Diabetic Retinopathy When Age at Diagnosis is Less Than 30 Years. Arch Ophthalmol (1984) 102:520–6. doi: 10.1001/archopht.1984.01040030398010

PubMed Abstract | CrossRef Full Text | Google Scholar

27. Brownlee M. Biochemistry and Molecular Cell Biology of Diabetic Complications. Nature (2001) 414:813–20. doi: 10.1038/414813a

PubMed Abstract | CrossRef Full Text | Google Scholar

28. Silva KC, Pinto CC, Biswas SK, Lopes de Faria JB, Lopes de Faria JM. Hypertension Increases Retinal Inflammation in Experimental Diabetes: A Possible Mechanism for Aggravation of Diabetic Retinopathy by Hypertension. Curr Eye Res (2007) 32:533–41. doi: 10.1080/02713680701435391

PubMed Abstract | CrossRef Full Text | Google Scholar

29. Basi S, Fesler P, Mimran A, Lewis JB. Microalbuminuria in Type 2 Diabetes and Hypertension: A Marker, Treatment Target, or Innocent Bystander? Diabetes Care (2008) 31:S194–201. doi: 10.2337/dc08-s249

PubMed Abstract | CrossRef Full Text | Google Scholar

30. Sowers JR, Epstein M, Frohlich ED. Diabetes, Hypertension, and Cardiovascular Disease: An Update. Hypertension (2001) 37:1053–9. doi: 10.1161/01.HYP.37.4.1053

PubMed Abstract | CrossRef Full Text | Google Scholar

31. Kajiwara A, Miyagawa H, Saruwatari J, Kita A, Sakata M, Kawata Y, et al. Gender Differences in the Incidence and Progression of Diabetic Retinopathy Among Japanese Patients With Type 2 Diabetes Mellitus: A Clinic-Based Retrospective Longitudinal Study. Diabetes Res Clin Pract (2014) 103:e7–10. doi: 10.1016/j.diabres.2013.12.043

PubMed Abstract | CrossRef Full Text | Google Scholar

32. Althouse AD, Abbott JD, Forker AD, Bertolet M, Barinas-Mitchell E, Thurston RC, et al. Risk Factors for Incident Peripheral Arterial Disease in Type 2 Diabetes: Results From the Bypass Angioplasty Revascularization Investigation in Type 2 Diabetes (BARI 2d) Trial. Diabetes Care (2014) 37:1346–52. doi: 10.2337/dc13-2303

PubMed Abstract | CrossRef Full Text | Google Scholar

33. Matsushita K, Ballew SH, Coresh J, Arima H, Ärnlöv J, Cirillo M, et al. Measures of Chronic Kidney Disease and Risk of Incident Peripheral Artery Disease: A Collaborative Meta-Analysis of Individual Participant Data. Lancet Diabetes Endocrinol (2017) 5:718–28. doi: 10.1016/S2213-8587(17)30183-3

PubMed Abstract | CrossRef Full Text | Google Scholar

34. Mohammedi K, Woodward M, Hirakawa Y, Zoungas S, Williams B, Lisheng L, et al. Microvascular and Macrovascular Disease and Risk for Major Peripheral Arterial Disease in Patients With Type 2 Diabetes. Diabetes Care (2016) 39:1796–803. doi: 10.2337/dc16-0588

PubMed Abstract | CrossRef Full Text | Google Scholar

35. Adler AI, Stevens RJ, Neil A, Stratton IM, Boulton AJ, Holman RR. UKPDS 59: Hyperglycemia and Other Potentially Modifiable Risk Factors for Peripheral Vascular Disease in Type 2 Diabetes. Diabetes Care (2002) 25:894–9. doi: 10.2337/diacare.25.5.894

PubMed Abstract | CrossRef Full Text | Google Scholar

36. Kumar S, Narwal S, Kumar V, Prakash O. α-Glucosidase Inhibitors From Plants: A Natural Approach to Treat Diabetes. Pharmacognosy Rev (2011) 5:19. doi: 10.4103/0973-7847.79096

CrossRef Full Text | Google Scholar

37. Ross SA, Dzida G, Vora J, Khunti K, Kaiser M, Ligthelm RJ. Impact of Weight Gain on Outcomes in Type 2 Diabetes. Curr Med Res Opin (2011) 27:1431–8. doi: 10.1185/03007995.2011.585396

PubMed Abstract | CrossRef Full Text | Google Scholar

38. Gomez-Peralta F, Abreu C, Gomez-Rodriguez S, Barranco RJ, Umpierrez GE. Safety and Efficacy of DPP4 Inhibitor and Basal Insulin in Type 2 Diabetes: An Updated Review and Challenging Clinical Scenarios. Diabetes Ther (2018) 9:1775–89. doi: 10.1007/s13300-018-0488-z

PubMed Abstract | CrossRef Full Text | Google Scholar

39. Giannarelli R, Aragona M, Coppelli A, Del Prato S. Reducing Insulin Resistance With Metformin: The Evidence Today. Diabetes Metab (2003) 29:6S28–35. doi: 10.1016/S1262-3636(03)72785-2

PubMed Abstract | CrossRef Full Text | Google Scholar

Keywords: diabetic complication, clinical characteristics, k-means, medication treatment, type 2 diabetes

Citation: Xing L, Peng F, Liang Q, Dai X, Ren J, Wu H, Yang S, Zhu Y, Jia L and Zhao S (2021) Clinical Characteristics and Risk of Diabetic Complications in Data-Driven Clusters Among Type 2 Diabetes. Front. Endocrinol. 12:617628. doi: 10.3389/fendo.2021.617628

Received: 15 October 2020; Accepted: 11 June 2021;
Published: 30 June 2021.

Edited by:

Hans Ulrich Häring, Tübingen University Hospital, Germany

Reviewed by:

Jean-Henri Calvet, Impeto Medical, Inc., United States
Stefan Kabisch, German Institute of Human Nutrition Potsdam-Rehbruecke (DIfE), Germany

Copyright © 2021 Xing, Peng, Liang, Dai, Ren, Wu, Yang, Zhu, Jia and Zhao. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.

*Correspondence: Shancen Zhao, zhaoshancen@genomics.cn; Lijing Jia, jialijing2012@126.com

These authors have contributed equally to this work

Disclaimer: All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article or claim that may be made by its manufacturer is not guaranteed or endorsed by the publisher.