Data-driven long-term glycaemic control trajectories and their associated health and economic outcomes in Finnish patients with incident type 2 diabetes

Abstract

Background

Treatments should be customized to patients to improve patients’ health outcomes and maximize the treatment benefits. We aimed to identify meaningful data-driven trajectories of incident type 2 diabetes patients with similarities in glycated haemoglobin (HbA1c) patterns since diagnosis and to examine their clinical and economic relevance.

Materials and methods

A cohort of 1540 patients diagnosed in 2011–2012 was retrieved from electronic health records covering primary and specialized healthcare in the North Karelia region, Finland. EHRs data were compiled with medication purchase data. Average HbA1c levels, use of medications, and incidence of micro- and macrovascular complications and deaths were measured annually for seven years since T2D diagnosis. Trajectories were identified applying latent class growth models. Differences in 4-year cumulative healthcare costs with 95% confidence intervals (CIs) were estimated with non-parametric bootstrapping.

Results

Four distinct trajectories of HbA1c development during 7 years after T2D diagnosis were extracted: patients with “Stable, adequate” (66.1%), “Slowly deteriorating” (24.3%), and “Rapidly deteriorating” glycaemic control (6.2%) as well as “Late diagnosed” patients (3.4%). During the same period, 2.2 (95% CI 1.9–2.6) deaths per 100 person-years occurred in the “Stable, adequate” trajectory increasing to 3.2 (2.4–4.0) in the “Slowly deteriorating”, 4.7 (3.1–6.9) in the “Rapidly deteriorating” and 5.2 (2.9–8.7) in the “Late diagnosed” trajectory. Similarly, 3.5 (95% CI 3.0–4.0) micro- and macrovascular complications per 100 person-years occurred in the “Stable, adequate” trajectory increasing to 5.1 (4.1–6.2) in the “Slowly deteriorating”, 5.5 (3.6–8.1) in the “Rapidly deteriorating” and 7.3 (4.3–11.8) in the “Late diagnosed” trajectory. Patients in the “Stable, adequate” trajectory had lower accumulated 4-year medication costs than other patients.

Conclusions

Data-driven patient trajectories have clinical and economic relevance and could be utilized as a step towards personalized medicine instead of the common “one-fits-for-all” treatment practices.

Introduction

Main target of type 2 diabetes (T2D) care is to achieve and maintain good treatment balance to avoid complications and premature deaths [1–4]. Typically, glycated haemoglobin A1c (HbA1c) is used as a measure of long-term blood sugar level with values <53 mmol/mol (7.0%) indicating good treatment balance. However, individualized, more-stringent goals for young T2D patients with short disease duration and less-rigorous, personalized treatment targets for older and the most comorbid patients are recommended [1, 4]. Treatments should be customized to patients to improve patients’ health outcomes and maximize the treatment benefits.

While the need of individual tailoring of treatments due to heterogeneity of T2D patients is recognized, majority of patients still receive standardized “one-size-fits-all” care due to limitations in available resources for individualized treatment. However, identification of patient groups with similarities in care needs and treatment outcomes, taking also into account budget constraints and cost consequences, is a step towards personalized care. It is also economically rational to identify subgroups where the treatment is cost-effective [5] and to offer care accordingly. Clustering of treatment outcomes, namely HbA1c, of T2D patients provides tools for treatment targeting and service design, and, further, optimizing the allocation of limited economic resources. Patients in different trajectories are likely to prefer different configurations of diabetes care and support. For example, patients with a deteriorating treatment balance profile should receive more intensive disease management supported by health professionals whereas adequate self-management support could be well enough for patients within stable treatment balance.

The objective of this study was to identify meaningful data-driven trajectories of incident T2D patients with similarities in HbA1c patterns since diagnosis. Further, we aimed to describe incidence of diabetes complications, all-cause mortality, and accumulation of healthcare costs by the trajectories to validate clinical applicability of the results in optimization and supporting effectiveness of treatment.

Materials and methods

Study setting

In Finland, most health services are public meaning all residents have the right to use and equal access to them. Municipalities are responsible for organizing and financing public health and social care services for their residents. Siun sote—Joint municipal authority for North Karelia social and health services was established in the region of North Karelia at the beginning of 2017 [6]. Instead of one municipality, Siun sote joins healthcare organisations of 13 municipalities with around 166 000 inhabitants. Before the integration to Siun sote, all the municipalities in the North Karelia region had introduced an identical electronic health register (EHR) called Mediatri in 2010–2011, and its use has been continued also after the integration. EHR data from the Mediatri cover both public primary and specialized healthcare. However, these EHRs do not cover patients utilizing only private healthcare services that produce about one quarter of all the health and social care services.

All community-dwelling residents of Finland are eligible for medication reimbursement, and reimbursed medication purchases from community pharmacies are registered in the Finnish Prescription Register (FPR) established in 1994 and maintained by the Social Insurance Institution of Finland (SII). However, patients staying at public nursing homes or hospitals >90 consecutive days are not eligible for medication reimbursement and their purchases are not registered during these stays. In addition, the FPR does not include records of over-the-counter medication purchases. The Special Reimbursement Register (SRR), established in 1964 and also maintained by the SII, includes records of entitlements to higher medication reimbursements due to certain chronic diseases, such as T2D.

Study design

This study utilized data provided by Siun sote, where the prevalence of T2D was 6.2% in 2012 [7]. All patients with newly diagnosed T2D in 2011–2012 were identified with ICD-10 code [8] E11 from the EHRs. Collected EHR data contained information on diagnoses, medication prescriptions, as well as laboratory assays during 2010–2017. EHR data were combined with information on medication purchases (such as dispensing date, Anatomical Therapeutic Chemical (ATC) classification code [9], and cost) from the FPR for the years 1995–2017. Entitlements to special refund for diabetes medications prior to 2011 were retrieved from the SRR.

Timing of T2D diagnosis was ascertained with data on diabetes medication purchases (ATC code A10) from the FPR in 1995–2010 and the retrieved entitlements to special refund for diabetes medications since 1964. If a patient had T2D diagnosed by healthcare providers outside of Siun sote, for example, before moving to the North Karelia region, the first record of T2D could be tracked from the national FPR and SRR instead of the first record in the EHRs of Siun sote. T2D diagnosis date was considered to be the first ever record of either diabetes medication purchase in the FPR, entitlement to special refund in the SRR, or T2D diagnosis in the EHRs, and only patients with an ascertained diagnosis recorded in 2011–2012 were included.

Glycaemic control

Glycaemic control was measured with HbA1c with the turbidimetric inhibition immunoassay method (TINIA) in the Eastern Finland laboratory (ISLAB, https://www.islab.fi). The method and its analysis were standardized across the North Karelia region. Values were standardized to the International Federation of Clinical Chemistry (IFCC) units. EHRs data on measurements from 3 months prior to T2D diagnosis until a maximum of 89 months after it were used. Patients with ≥1 HbA1c measurement during the study period were included. Baseline HbA1c (mmol/mol) level was defined to be the one measured ±3 months from T2D diagnosis. Follow-up measures were categorized as occurring 6–17 (1 year), 18–29 (2 years), 30–41 (3 years), 42–53 (4 years), 54–65 (5 years), 66–77 (6 years), and 78–89 months (7 years) after diagnosis. If several measurements occurred within a period, the latest one was selected. Sensitivity analyses were conducted against the average of the HbA1c measurements within a period and the results were comparable. If HbA1c value was unmeasured within a certain period, it was considered as a missing value. HbA1c values were measured until the date of death, moving outside of the study area, or 31^st December 2017, whichever occurred first.

Deaths and diabetes complications

Information on the dates of deaths and complications occurring after T2D diagnosis were retrieved from the EHRs. Micro- and macrovascular complications were identified with ICD-10 codes which have been in use in Finland since 1996 (S1 Table). Complications were examined as a composite of micro- and macrovascular complications as well as separately. Incidence of all-cause deaths and diabetes complications were followed up until the date of the specific event, moving outside of the study area, or 31^st December 2017, whichever occurred first.

Healthcare costs

Healthcare costs were calculated as a cumulative sum of primary and specialized healthcare costs, nursing home and home care costs, and outpatient medication costs during 2014–2017. EHR data included all hospitalizations (outpatient visits, inpatient admissions) and their costs derived from the hospitals’ cost accounting systems in 2014–2017. Restriction in cost analyses to years 2014–2017 was applied because data on outpatient visits were not complete and, therefore, comparable between municipalities during the transitional phase to the Mediatri before 2014. The individual-level resource use data for primary care were also obtained from the EHRs (with diagnosis and activity information) and grouped using the Ambulatory and Primary Care-Related Patient Groups (APR) grouper software, a grouping system equivalent to DRG used in hospital care. The unit cost estimates for grouped primary care contacts were obtained from the national standard price lists for primary care encounters [10]. Outpatient medication costs were calculated based on the information on the total costs of all reimbursed as well as reimbursed antidiabetic (ATC code A10) medications only from the FPR. In addition to the combined total costs, accumulated social- and healthcare costs and medication costs were examined separately. Costs were estimated from the payer’s perspective, including all direct health and social care costs.

Baseline characteristics

Patients’ age and gender were obtained from the EHRs. Concordant and discordant diseases that occurred prior to the T2D diagnosis date were defined from the EHRs (S2 Table). Micro- and macrovascular comorbidities occurring prior to T2D diagnosis were identified utilizing the same definitions as for the outcome (S1 Table). Statin (ATC code C10A) use was identified from the FPR.

For body mass index (BMI), blood pressure, fasting plasma glucose (FPG), and creatinine, measurements occurring most closely to the T2D diagnosis date measured before or within 3 months from the T2D diagnosis, and for low-density lipoprotein (LDL), measurements occurring most closely to the T2D diagnosis date measured before or within 1 month from the diagnosis were extracted from the EHRs. Based on these definitions, indicators of having LDL, BMI, FPG, blood pressure, creatinine and HbA1c measured at the baseline were created. In addition, estimated glomerular filtration rate (eGFR) was calculated with the Chronic Kidney Disease Epidemiology Collaboration (CKD-EPI) equation [11] based on creatinine, age, and gender. eGFR was further classified into normal/high (stage I, >90 ml/min/1.73 m²), mildly decreased (stage II, 60–90 ml/min/1.73 m²), mildly to severely decreased (stage III, 30–59 ml/min/1.73 m²) and severely decreased/kidney failure (stage IV, <30 ml/min/1.73 m²).

Diabetes medication use during the follow-up

Diabetes medication use was examined from the FPR for 7 years since T2D diagnosis. The following subgroups of medication use were formed: users of metformin (ATC code A10BA02) only, users of metformin and other oral antidiabetic medications (A10BA02 + other A10B), users of only other diabetes medications than insulin or metformin (A10B excluding A10BA02 and A10A), users of insulin and oral antidiabetic medication (including metformin) (A10A + A10B), and users of insulin (A10A) only. Patients who died were taken into account in the calculation of prevalence of diabetes medication use.

Statistical analyses

Patients were clustered into distinct trajectories of glycaemic control with latent class growth analysis (LCGA) [12–15]. Models were fitted iteratively with 1–5 classes and varying shapes for the trajectories (linear, quadratic, cubic). LCGAs were estimated with full-information maximum likelihood and missing data were not imputed but all available data were used. The method assumes no variation within the trajectories. Selection of the best model was done utilizing information from fit indices and classification performance of the model as well as clinical interpretation of the trajectories. Bayesian Information Criteria (BIC) was utilized to compare models with lower values indicating better fit. Low-Mendel-Rubin likelihood ratio test with P value <0.05 was considered to indicate better fit for a n-class model than for n–1-class model. Entropy was used to guide in the classification accuracy of the model with higher values indicating better classification. Small classes with <2% of population were not accepted to have large enough groups of at least 30 patients. LCGAs were performed using Mplus [16] Version 8.

Differences in baseline characteristics between the identified trajectories were examined with one-way ANOVA or nonparametric Kruskal-Wallis test for continuous variables and chi squared test for categorical ones. Association of baseline characteristics with trajectory membership was examined with multinomial logistic regression model. Cox regression models were utilized to describe incidence of deaths and micro- and macrovascular complications by the trajectories. Models were adjusted for age, sex, concordant and discordant diseases, and microvascular and macrovascular comorbidities at baseline. Calendar year- and trajectory-specific healthcare costs were weighted with Kaplan-Meier sample average estimator [17] to account for differences between trajectories in censoring due to death in calculation of adjusted cumulative costs during 2014–2017. Further, non-parametric bootstrapping generating 1000 bootstrap samples with replacement was applied to estimate mean differences in adjusted cumulative healthcare costs with 95% confidence intervals (CIs). Analyses were conducted with SAS version 9.4 (SAS Institute Inc., Cary, North-Carolina, USA). Two-tailed P values < 0.05 were considered statistically significant.

Ethics statement

Use of the data was approved by the Ethics Committee of the Northern Savonia Hospital District (diary number 81/2012) that considers all applications within its university hospital catchment area. The study protocol was also approved by the register administrator, Siun sote—Joint municipal authority for North Karelia social and health services. A separate permission to link data on medication purchases and special reimbursements was achieved from the Social Insurance Institute (diary number 110/522/2018). We utilized only de-identified register-based data and individuals in the registers were not contacted. Thus, no written consent from the patients was required.

Results

In total, 1647 T2D patients diagnosed in 2011–2012 were identified. After removing patients without HbA1c measurements during the follow-up, 1540 (93.5%) patients remained with 14,173 HbA1c measurements. Mean (SD) age of the patients was 64.9 (12.7) years and 46.5% were female (Table 1). 935 (60.7%) had HbA1c measured ±3 months of T2D diagnosis with an average HbA1c of 48.5 (16.5) mmol/mol.

[Figure omitted. See PDF.]

Of the fitted LCGAs (S3 Table), four-class cubic model was selected as the best fitting with clinically interpretable results (Fig 1). The extracted trajectories were: patients with “Stable, adequate” (66.1%), “Slowly deteriorating” (24.3%), and “Rapidly deteriorating” glycaemic control (6.2%) as well as “Late diagnosed” patients (3.4%).

[Figure omitted. See PDF.]

Patients in the”Rapidly deteriorating” trajectory were more often male and younger than patients in the”Stable, adequate” trajectory (Table 2). Having HbA1c measured at the time of the T2D diagnosis increased and having LDL measured decreased the odds of belonging to the “Slowly deteriorating” and “Late diagnosed” trajectories compared with the”Stable, adequate” trajectory. Having BMI measured decreased the odds of belonging to the “Rapidly deteriorating” and “Late diagnosed” trajectories compared with the”Stable, adequate” trajectory. Patients in the”Late diagnosed” trajectory had less concordant diseases and patients in the “Slowly deteriorating” trajectory used less statins than patients in the”Stable, adequate” trajectory. Average number of HbA1c measurements during the study period differed between the trajectories (Table 1).

[Figure omitted. See PDF.]

Stable, adequate glycaemic control trajectory (n = 1058) as a reference.

Diabetes medication use differed between the trajectories and the content of medication structure changed during the 5-year follow-up (Fig 2, S4 Table). In the “Stable, adequate” trajectory, half of the patients did not purchase diabetes medications during the follow-up and metformin use was most frequent. Among patients in the “Rapidly deteriorating” trajectory, treatments were clearly intensified throughout the follow-up. Use of insulin was the most frequent in the “Late diagnosed” trajectory with almost half of the patients using it, remaining stable over the follow-up. Non-use of diabetes medications decreased in every trajectory during the follow-up.

[Figure omitted. See PDF.]

OAD, oral antidiabetic drugs or GLP-1 analogues (incl. metformin, sulfonylureas, combinations of oral blood glucose lowering drugs, glitazones, DPP-4 inhibitors, glinides, GLP-1 analogues, and SGLT2 inhibitors).

In total, 226 deaths and 318 complications (125 microvascular and 237 macrovascular events) occurred during the follow-up (Fig 3, S1 and S5 Tables). Incidence rate of death after T2D diagnosis was 2.2 (95% CI 1.9–2.6) events per 100 person-years among “Stable, adequate” trajectory, increasing to 3.2 (2.4–4.0) in “Slowly deteriorating”, 4.7 (3.1–6.9) in “Rapidly deteriorating” and 5.2 (2.9–8.7) in the “Late diagnosed” trajectories. Hazard ratio for death was 2.86 (95% CI 1.60–5.12) in the “Late diagnosed” and 2.49 (95% CI 1.61–3.87) in the “Rapidly deteriorating” trajectory compared with the “Stable, adequate” trajectory (S6 Table). Similarly, incidence rate of micro- and macrovascular complications after T2D diagnosis was 3.5 (95% CI 3.0–4.0) events per 100 person-years among “Stable, adequate” trajectory, increasing to 5.1 (4.1–6.2) in “Slowly deteriorating”, 5.5 (3.6–8.1) in “Rapidly deteriorating” and 7.3 (4.3–11.8) in the “Late diagnosed” trajectories. Hazard ratio for complications was 3.35 (95% CI 1.96–5.72) in the “Late diagnosed”, 2.28 (95% CI 1.51–3.46) in the “Rapidly deteriorating”, and 1.39 (95% CI 1.08–1.79) in the “Slowly deteriorating” compared with the “Stable, adequate” trajectory.

[Figure omitted. See PDF.]

(A) Cumulative incidence of death since T2D diagnosis. (B) Cumulative incidence of combined event of micro- and macrovascular complications since T2D diagnosis.

Accumulated total medication costs differed between the trajectories while no difference was seen in accumulated total costs or accumulated social- and healthcare costs (Table 3, S6 Table). However, patients in the “Late diagnosed” trajectory had lower outpatient visit costs in both primary and specialized healthcare than patients in the “Stable, adequate” trajectory. In addition, patients in the “Stable, adequate” trajectory had lower accumulated diabetes medication costs than patients in other trajectories.

[Figure omitted. See PDF.]

Mean differences in accumulated 4-year per-person healthcare costs between HbA1c trajectories with bootstrapped 95% confidence intervals (Stable, adequate glycaemic control trajectory as a reference). All costs are Kaplan-Meier sample average adjusted and calculated for 2014–2017.

Discussion

Based our study results, four trajectories among incident Finnish T2D patients were identified with data-driven technique to classify patients following similar HbA1c patterns over time. Identified trajectories had clinical relevance in terms of diabetes medication use, incidence of diabetes complications and deaths as well as accumulation of medication costs during the study. Monitoring of HbA1c and LDL at the time of T2D diagnosis differed between the trajectories. To our knowledge, novelty of the present study is in the examination of direct long-term costs of social- and healthcare by the trajectories among incident T2D patients in Finland. Due to the restricted healthcare resources, it is crucial to understand the costs associated with the trajectories to give reasons for the most effective care of T2D to avoid costly complications.

Trajectories estimated in this study resemble those identified in previous publications. Typically, 2–5 trajectories have been identified among incident T2D populations [18–21]. More studies have been conducted among prevalent T2D patients [22–27]. Generally, comparison with these previous studies is complicated because of differences in statistical methods used to identify trajectories, follow-up times as well as calendar times: new medications have entered the markets between the studies intensifying the treatment of T2D. In addition, as care and treatment of diabetes varies between geographical areas, and even within one country depending on healthcare providers in addition to differences between populations, results may not be universally generalizable. To note, also clustering based on patient characteristics (age, HbA1c, BMI, homeostasis model assessment (HOMA2) estimates of β-cell function and insulin resistance) at the time of diagnosis in general diabetic population has been examined [28] instead of longitudinal clustering based on HbA1c solely.

Patients in the “Late diagnosed” trajectory had less concordant diseases and measurements of LDL than patients in the “Stable, adequate” trajectory at the time of T2D diagnosis. In addition, having HbA1c measured closely to the T2D diagnosis date increased the odds to be late diagnosed T2D patient compared with the “Stable, adequate” trajectory. These indicate that T2D diagnosis may be delayed due to nonexistence of other diseases or healthcare problems requiring regular controls. In the previous studies [18–21], younger age, male sex, and higher BMI have been found to be associated with poorer glycaemic control trajectories similarly to our study.

To examine clinical relevance of the identified trajectories, we described incidences of all-cause deaths and diabetes complications as well as accumulation of healthcare costs by the trajectories. Our results are in line with a systematic review [29] of studies utilizing longitudinal HbA1c clustering and examining their associations with outcomes reporting poor glycaemic control to be associated with the risk of micro- and macrovascular events and mortality. HbA1c may not reflect only treatment response but also the underlying severity of the disease [28]. Late diagnosed patients used most often insulin reflecting more severe disease, seen also as higher incidence of complications during the 7-year follow-up despite of reaching good glycaemic control 2–3 years after the diagnosis. We observed that accumulated medication costs differed between the trajectories and poorer glycaemic control was associated with increased costs which could be considered as a natural consequence of more severe disease, treatment intensification and incidence of complications.

Medical treatments have a central role in T2D care. Newer diabetes medications, i.e., SGLT2 inhibitors and GLP-1 receptor agonists, were reported to lower the risk of cardiovascular events and subsequent mortality during the 2010’s [3, 4]. Use of these medications, especially SGLT2 inhibitors, started to rise in 2016 in Finland. In our cohort, these medications were most often used by patients from the “Rapidly deteriorating” trajectory (S4 Table). In addition, increased use was observed within the “Late diagnosed” and “Slowly deteriorating” trajectories. When examining the risk of micro- and macrovascular complications and mortality, it seems that our analyses may, therefore, underestimate the underlying risk caused by poorer glycaemic control when compared with the “Stable, adequate” trajectory. Of the complications studied in this study, chronic kidney disease may have been somewhat underdiagnosed leading to underestimation of microvascular diseases at baseline and during the follow-up. As the disease is highly related to the mortality risk, residual confounding in mortality analyses is possible although the analyses accounted for end-stage renal disease.

Strengths of our study are inclusion of all patients with T2D diagnosed in 2011–2012 in the North Karelia region. In addition, all the municipalities in North Karelia use the same regional laboratory and the same standardized methods for HbA1c testing. Due to the Finnish public healthcare system available for all inhabitants, we were able to continuously follow patients’ public healthcare service use until death, moving outside of the study area, or administrative study end of the study. Utilization of register-based data allowed us to avoid recall bias, too. Validation studies of the national Finnish Care Register for Health Care, which also includes EHRs from the North Karelia region, have reported positive predictive values of 87–97% for hospital discharges of vascular diseases [30]. Compared with the previous studies, we were also able to examine economic outcomes in addition to health outcomes.

Our study has some limitations. We had a relatively small sample of incident T2D patients stratified to trajectories which limits the conclusions we can draw from the results due to wide confidence intervals especially in the smallest trajectories, although point estimates show some interesting tendencies. For this, we were also unable to conduct an internal validation study. Although we ascertained the timing of T2D diagnosis observed in the EHRs with data from the SII, the date of diagnosis may not be accurate. We had scarce data on laboratory measurements at the baseline meaning that we could not examine the effect of some important factors, such as LDL or BMI, as such on trajectory probability. In addition, according to our previous study, younger patients with T2D were poorly monitored compared with older patients in the North Karelia region [31]. We did not incorporate information on the diabetic medications in use to our analyses but focused solely on HbA1c levels and associated outcomes. However, differences are reported between medications in their efficacy to manage HbA1c [3, 4].

Interpretation of costs should be cautious. Cumulative healthcare costs did not include costs from the first 1–2 years of follow-up since diagnosis. However, as depicted from Fig 3B, differences in the incidence of micro- and macrovascular complications between the trajectories were already apparent at that time possibly causing some bias, especially to specialized healthcare cost analyses. Medication costs did not include costs of medications used during a long-term (>90 days) hospital stays because they are not registered in the FPR. Furthermore, we did not have information on patients using private healthcare services. However, as persons utilizing these services are not likely to use only private healthcare services in the care of T2D due to, for example, economic issues, this is not a big concern. Finally, we aimed to describe the incidence of complications and deaths by the trajectories. Thus, the results presented do not imply causal relationships since trajectories and incidence of the events were measured simultaneously.

Data-driven HbA1c trajectories among incident T2D patients may be utilized as a first step towards individualized care. Repeated measurements of HbA1c are available from EHRs enabling monitoring of patients’ HbA1c levels during a longer period than one measurement point that is often used as a basis for treatment decisions. With the results of this study, clinicians are provided with more information to detect and treat patients whose diabetes care may need more attention. In addition, as patients have a central role in managing their diabetes care, demonstrating their future scenarios including risk of complications, death and healthcare costs may motivate them to engage more intensively to diabetes care. These aspects enable more personalized treatment of T2D patients. Further, this study shows the importance of proper diabetes care and glycaemic management for the society via the association between glycaemic control and medication costs. According to the results, special effort should be allocated to early detection of glycaemic control of late diagnosed T2D patients with high HbA1c values at the time of diagnosis and those with rapid deterioration in treatment balance to decrease the incidence of complications and medication costs. In future, this study should be replicated with a larger sample of T2D patients and longer follow-up to examine causal relationships.

Supporting information

S1 Table. Definitions of micro- and macrovascular complications and numbers of outcome events by type of analysis (first event occurrence after type 2 diabetes diagnosis).

https://doi.org/10.1371/journal.pone.0269245.s001

S2 Table. Definitions of concordant and discordant diseases at baseline.

https://doi.org/10.1371/journal.pone.0269245.s002

S3 Table. Unadjusted latent class growth analyses for HbA1c values measured annually since T2D diagnosis among patients diagnosed in 2011 or 2012 (n = 1540).

https://doi.org/10.1371/journal.pone.0269245.s003

S4 Table. Antidiabetic medication use (%) by the trajectories after the T2D diagnosis.

https://doi.org/10.1371/journal.pone.0269245.s004

S5 Table. Associations between the estimated HbA1c trajectories and incidence of deaths and diabetes complications.

https://doi.org/10.1371/journal.pone.0269245.s005

S6 Table. Average accumulated 4-year per-person healthcare costs (euros, €) by HbA1c trajectories in 2014–2017.

All the costs are adjusted for censoring due to death utilizing Kaplan-Meier sample average estimator.

https://doi.org/10.1371/journal.pone.0269245.s006

(PDF)

Citation: Lavikainen P, Aarnio E, Linna M, Jalkanen K, Tirkkonen H, Rautiainen P, et al. (2022) Data-driven long-term glycaemic control trajectories and their associated health and economic outcomes in Finnish patients with incident type 2 diabetes. PLoS ONE 17(6): e0269245. https://doi.org/10.1371/journal.pone.0269245

About the Authors:

Piia Lavikainen

Roles: Conceptualization, Formal analysis, Methodology, Software, Visualization, Writing – original draft, Writing – review & editing

E-mail: [email protected]

Affiliation: School of Pharmacy, University of Eastern Finland, Kuopio, Finland

https://orcid.org/0000-0003-3560-8860

Emma Aarnio

Roles: Conceptualization, Writing – review & editing

Affiliation: School of Pharmacy, University of Eastern Finland, Kuopio, Finland

Miika Linna

Roles: Formal analysis, Resources, Writing – review & editing

Affiliation: Aalto University, Helsinki, Finland

Kari Jalkanen

Roles: Conceptualization, Writing – review & editing

Affiliation: School of Pharmacy, University of Eastern Finland, Kuopio, Finland

Hilkka Tirkkonen

Roles: Supervision, Writing – review & editing

Affiliation: Siun Sote – Joint Municipal Authority for North Karelia Social and Health Services, Joensuu, Finland

Päivi Rautiainen

Roles: Supervision, Writing – review & editing

Affiliation: Siun Sote – Joint Municipal Authority for North Karelia Social and Health Services, Joensuu, Finland

Tiina Laatikainen

Roles: Conceptualization, Funding acquisition, Project administration, Resources, Supervision, Writing – review & editing

Affiliations Siun Sote – Joint Municipal Authority for North Karelia Social and Health Services, Joensuu, Finland, Department of Public Health Solutions, National Institute for Health and Welfare, Helsinki, Finland, Institute of Public Health and Clinical Nutrition, University of Eastern Finland, Kuopio, Finland

Janne Martikainen

Roles: Conceptualization, Project administration, Supervision, Writing – review & editing

Affiliation: School of Pharmacy, University of Eastern Finland, Kuopio, Finland

References

1. Type 2 diabetes. Current Care Guidelines. Working group set up by the Finnish Medical Society Duodecim, the Finnish Society of Internal Medicine and the Medical Advisory Board of the Finnish Diabetes Society. Helsinki: The Finnish Medical Society Duodecim, 2020 (referred August 6, 2020). https://www.kaypahoito.fi/hoi50056#K1 (22 April 2022, date last accessed).

2. Arnett DK, Blumenthal RS, Albert MA, Buroker AB, Goldberger ZD, Hahn EJ, et al. 2019 ACC/AHA Guideline on the Primary Prevention of Cardiovascular Disease: A Report of the American College of Cardiology/American Heart Association Task Force on Clinical Practice Guidelines [published correction appears in Circulation. 2019;140(11): e649-e650] [published correction appears in Circulation. 2020;141(4): e60]. Circulation. 2019;140: e596–e646. pmid:30879355

3. Cosentino F, Grant PJ, Aboyans V, Bailey CJ, Ceriello A, Delgado V, et al. 2019 ESC Guidelines on diabetes, pre-diabetes, and cardiovascular diseases developed in collaboration with the EASD. Eur Heart J. 2020;41: 255–323. pmid:31497854

4. Buse JB, Wexler DJ, Tsapas A, Rossing P, Mingrone G, Mathieu C, et al. 2019 update to: Management of hyperglycaemia in type 2 diabetes, 2018. A consensus report by the American Diabetes Association (ADA) and the European Association for the Study of Diabetes (EASD). Diabetologia. 2020;63: 221–228. pmid:31853556

5. Coyle D, Buxton MJ, O’Brien BJ. Stratified cost-effectiveness analysis: a framework for establishing efficient limited use criteria. Health Econ. 2003;12: 421–427. pmid:12720259

6. Wikström K, Toivakka M, Rautianen P, Tirkkonen H, Kivinen P, Laatikainen T. The effect of the integration of health services on health care usage among patients with type 2 diabetes in North Karelia, Finland. BMC Health Serv Res. 2021:21; 65. pmid:33441132

7. Tirkkonen H, Sikiö M, Kekäläinen P, Laatikainen T. Results of treatment of type 2 diabetes in primary health care (in Finnish with English abstract). Finnish Med J. 2014;34: 2027–2032.

8. WHO. International Statistical Classification of Diseases and Related Health Problems 10th revision. http://apps.who.int/classifications/icd10/browse/2010/en. Accessed October 10, 2019.

9. WHO Collaborating Centre for Drug Statistics Methodology, ATC classification index with DDDs, 2019. Oslo, Norway 2018. http://www.whocc.no/atc_ddd_index/. Accessed November 19, 2019.

10. Kapiainen S, Väisänen A, Haula T. Terveyden- ja sosiaalihuollon yksikkökustannukset Suomessa vuonna 2011 (in Finnish). Finnish Institute for Health and Welfare, report 3/2014. http://www.urn.fi/URN:ISBN:978-952-302-079-5. (20 May 2020, date last accessed).

11. Levey AS, Stevens LA, Schmid CH, Zhang YL, Castro AF 3rd, Feldman HI, et al. CKD-EPI (Chronic Kidney Disease Epidemiology Collaboration). A new equation to estimate glomerular filtration rate. Ann Intern Med 2009;150: 604–612.

12. Laird NM, Ware JH. Random-effects models for longitudinal data. Biometrics. 1982;38: 963–974. pmid:7168798

13. Muthén B, Shedden K. Finite mixture modelling with mixture outcomes using the EM algorithm. Biometrics. 1999;55: 463–469. pmid:11318201

14. Jung T, Wickrama KAS. An introduction to latent class growth analysis and growth mixture modeling. Soc Personal Psychol Compass. 2008;2/1: 302–317.

15. Nagin DS, Odgers CL. Group-based trajectory modelling in clinical research. Annu Rev Clin Psychol. 2010;6: 109–138. pmid:20192788

16. Muthén LK, Muthén BO. Mplus User’s Guide, Seventh Edition ed. Muthén & Muthén: Los Angeles, CA; 1998–2015.

17. Gray AM, Clarke PM, Wolstenholme JL, Wordsworth S. Applied Methods of Cost-effectiveness Analysis in Healthcare. Oxford University Press, New York, United States, 2011.

18. Rathmann W, Schwandt A, Hermann JM, Kuss O, Roden M, Laubner K, et al. Distinct trajectories of HbA1c in newly diagnosed Type 2 diabetes from the DPV registry using a longitudinal group-based modelling approach. Diabet Med. 2019;36: 1468–1477. pmid:31392761

19. Hertroijs DFL, Elissen AMJ, Brouwers MCGJ, Schaper NC, Köhler S, Popa MC, et al. A risk score including body mass index, glycated haemoglobin and triglycerides predicts future glycaemic control in people with type 2 diabetes. Diabetes Obes Metab. 2018;20(3): 681–688. pmid:29095564

20. Laiteerapong N, Karter AJ, Moffet HH, Cooper JM, Gibbons RD, Liu JY, et al. Ten-year hemoglobin A1c trajectories and outcomes in type 2 diabetes mellitus: The Diabetes & Aging Study. J Diabetes Complications. 2017;31: 94–100. pmid:27503405

21. Hemo B, Geva D, Shahar DR, Golan R, Heymann AD. Distinct trajectories in HbA1c are associated with different all-cause mortality and morbidity in newly diagnosed patients with type 2 diabetes. Prim Care Diabetes. 2019;10: S1751-9918(19)30350-X. pmid:31836307

22. Luo M, Tan CS, Lim WY, Chia KS, Tang WE, Tai ES, et al. Association of diabetes treatment with long-term glycemic patterns in patients with type 2 diabetes mellitus: A prospective cohort study. Diabetes Metab Res Rev. 2019;35: e3122. pmid:30600922

23. Luo M, Lim WY, Tan CS, Ning Y, Chia KS, van Dam RM, et al. Longitudinal trends in HbA1c and associations with comorbidity and all-cause mortality in Asian patients with type 2 diabetes: A cohort study. Diabetes Res Clin Pract. 2017;133: 69–77. pmid:28898713

24. Walraven I, Mast MR, Hoekstra T, Jansen AP, van der Heijden AA, Rauh SP, et al. Distinct HbA1c trajectories in a type 2 diabetes cohort. Acta Diabetol. 2015;52: 267–275. pmid:25287012

25. Ravona-Springer R, Heymann A, Schmeidler J, Moshier E, Godbold J, Sano M, et al. Trajectories in glycemic control over time are associated with cognitive performance in elderly subjects with type 2 diabetes. PLoS One. 2014;9: e97384 pmid:24887092

26. Chang HY, Wahlqvist ML, Liu WL, Lee MS, Shin SJ, Li YS, et al. Management trajectories in the type 2 diabetes Integrated Delivery System project in Taiwan: accounting for behavioural therapy, nutrition education and therapeutics. Asia Pac J Clin Nutr. 2014;23: 592–606. pmid:25516317

27. Karpati T, Leventer-Roberts M, Feldman B, Cohen-Stavi C, Raz I, Balicer R. Patient clusters based on HbA1c trajectories: A step toward individualized medicine in type 2 diabetes. PLoS One. 2018;13: e0207096. pmid:30427908

28. 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(5): 361–369. pmid:29503172

29. Luo M, Tan KHX, Tan CS, Lim WY, Tai ES, Venkataraman K. Longitudinal trends in HbA1c patterns and association with outcomes: A systematic review. Diabetes Metab Res Rev. 2018;34: e3015. pmid:29663623

30. Sund R. Quality of the Finnish Hospital Discharge Register: A systematic review. Scand J Public Health. 2012;40: 505–515. pmid:22899561

31. Nazu NA, Lindström J, Rautiainen P, Tirkkonen H, Wikström K, Repo T, et al. Maintenance of good glycaemic control is challenging—A cohort study of type 2 diabetes patient in North Karelia, Finland. Int J Clin Pract. 2019;73: e13313. pmid:30664318